Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite strains”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Sound velocities of iron-nickel (Fe90Ni10) alloy up to 8 GPa and 773 K: The effect of nickel on the elastic properties of bcc-iron at high P-T

Sound velocities of iron and iron-based alloys at high pressure and high temperature are crucial for understanding the composition and structure of Earth’s and other telluric planetary cores. In this study, we performed ultrasonic interferometric measurements of both compressional (ν P ) and shear (ν S ) velocities on a polycrystalline body-centered-cubic (bcc)-Fe 90 Ni 10 up to 8 GPa and 773 K. The elastic moduli and their pressure and temperature derivatives are derived from least-square fits to third-order finite strain equations, yielding K S0 = 154.2(8) GPa, G 0 = 73.2(2) GPa, K' S0 = 4.6(2), G' 0 = 1.5(1), ∂K S /∂T = –0.028(1) GPa/K, and ∂G/∂T = –0.023(1) GPa/K. Furthermore, a comparison with literature data on bcc-Fe suggests that nickel not only decreases both P and S wave velocities but also weakens the temperature effects on the elastic moduli of Fe-Ni alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Reassessment of Birch's Law on hcp‐Fe From Ultrasonic Sound Velocity Measurement and Implications on the Velocity Profiles of Earth's Inner Core

Here, we performed in situ X-ray diffraction and ultrasonic sound velocity measurements on hcp-Fe up to 15 GPa, 873 K in a multi-anvil apparatus. The elastic moduli and their pressure and temperature derivatives were determined by fitting the velocity and density data to the third-order finite strain equations, yielding K S0 = 169.0(57) GPa, K$^{′}_{S0}$ = 5.4(6), (∂K S /∂T) P = −0.031(3) GPa/K, G 0 = 104.5(27) GPa, G$^{′}_{0}$ = 1.7(2), (∂G/∂T) P = −0.060(2) GPa/K. Within the experimental P-T range, we find significant temperature effect on the density-velocity (ρ − V P or S ) relations and caution the use of the temperature-independent Birch's law for extrapolation to core conditions. Furthermore, temperature-induced velocity decrease is more significant in V S than in V P , offering a possible explanation for the high Poisson's ratio in the core. Extrapolations based on our results together with previous experimental data suggest that V P of hcp-Fe aligns with PREM at Earth's core conditions, while Vs and density are approximately 10% and 2.7% higher than PREM, respectively.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Exact solutions of burst pressure for thick-walled cylinders in power-law strain hardening steels

Exact solutions of burst pressure for defect-free, thick-walled cylindrical pressure vessels with capped ends are developed using the flow theory of plasticity in terms of the Tresca, von Mises, and Zhu-Leis yield criteria. The material is assumed to obey a power-law strain hardening rule, and the finite strain theory is adopted to describe large plastic deformation. On this basis, internal pressure is obtained as a complex polynomial function of effective strains on the inside and outside surfaces of the pressure vessel with the use of the Bernoulli numbers. At burst failure, the effective strains and stresses on the inside and outside surfaces are first determined, and then three flow solutions of burst pressure are obtained as a power series function of the diameter ratio (D o /D i ), strain hardening exponent (n), and ultimate tensile strength (UTS). The power series solution is confirmed to agree well with the exact flow solution of burst pressure in a closed-form for each yield criterion, and the results show that the von Mises flow solution is an upper bound prediction, the Tresca flow solution is a lower bound prediction, and the Zhu-Leis flow solution is an intermediate prediction of burst pressure. In conclusion, two sets of full-scale burst test data are then utilized to evaluate and validate the proposed flow solutions of burst pressure for both thin-walled pipes and thick-walled cylindrical pressure vessels.

36 MATERIALS SCIENCE↗
NASA NTRS1971

Gravity flow of powder in a lunar environment. Part 2: Analysis of flow initiation

A small displacement-small strain finite element technique utilizing the constant strain triangle and incremental constitutive equations for elasticplastic (media nonhardening and obeying a Coulomb yield condition) was applied to the analysis of gravity flow initiation. This was done in a V-shaped hopper containing a powder under lunar environmental conditions. Three methods of loading were examined. Of the three, the method of computing the initial state of stress in a filled hopper prior to drawdown, by adding material to the hopper layer by layer, was the best. Results of the analysis of a typical hopper problem show that the initial state of stress, the elastic moduli, and the strength parameters have an important influence on material response subsequent to the opening of the hopper outlet.

Pariseau, W. G.↗
DOE OSTI2023

Topology optimization of compliant mechanisms under transient thermal conditions

This work considers multi-material topology optimization of compliant mechanisms under transient thermal and quasi-static mechanical conditions wherein thermally actuated devices are optimized for different operating conditions. The materials are modeled using finite strain thermo-hyperelasticity and a two way coupling between the energy balance and equilibrium equations is investigated. The design updates are generated from the gradient-based method of moving asymptotes optimizer and the sensitivities are computed using the time dependent adjoint sensitivity analysis. Results show the impact of designing for short versus long actuation times.

42 ENGINEERING↗
DOE OSTI2023

Computation of the Strain Induced Critical Current Reduction in the 16 T Nb 3 Sn Test Facility Dipole

A test facility dipole is being developed at LBNL, targeting a 16 T field in a 144 mm wide aperture. The magnet uses a block design, with two double-pancake coils. In order to minimize motion under the large Lorentz forces, the coils are preloaded against a thick aluminum shell and iron yoke using bladder and key technology. It is then crucial to verify that the performance of the magnet is not degraded due to strain induced on the Nb 3 Sn conductor during assembly, cool-down and powering. The critical current of extracted strands was measured in a varying background magnetic field and as a function of the applied longitudinal strain. Finite element analysis was used to extract the strain state inside the superconducting strands during magnet assembly and operation. This strain was then compared to the measurements to evaluate potential reversible and irreversible effects on the magnet performances. Here the results suggest that the magnet can reach 16 T with sufficient margin, with no irreversible degradation in the high field region.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1975

A triangular thin shell finite element: Nonlinear analysis

Aspects of the formulation of a triangular thin shell finite element which pertain to geometrically nonlinear (small strain, finite displacement) behavior are described. The procedure for solution of the resulting nonlinear algebraic equations combines a one-step incremental (tangent stiffness) approach with one iteration in the Newton-Raphson mode. A method is presented which permits a rational estimation of step size in this procedure. Limit points are calculated by means of a superposition scheme coupled to the incremental side of the solution procedure while bifurcation points are calculated through a process of interpolation of the determinants of the tangent-stiffness matrix. Numerical results are obtained for a flat plate and two curved shell problems and are compared with alternative solutions.

Thomas, G. R.↗
DOE OSTI2023

Exact solution of burst pressure for thick-walled pipes using the flow theory of plasticity

This paper develops an exact solution of burst pressure for defect-free, thick-walled pipes using the flow theory of plasticity in terms of the average shear stress yield criterion (or Zhu-Leis yield criterion). The pipe steel is assumed to obey the power-law strain hardening rule, and large plastic deformation is described by the finite strain theory. On this basis, internal pressure is obtained as a power series function of the effective strains on the inside and outside surfaces of the thick-walled pipe with use of the second Bernoulli numbers. At burst failure, the Zhu-Leis flow solution of burst pressure is determined as a power series solution, and the burst effective strains, the burst effective stresses, and the burst pressure are functions of the diameter ratio (D o / D i ), strain hardening exponent (n), and ultimate tensile strength (UTS). Similarly, the Tresca and von Mises flow solutions of burst pressure are also determined as a power series solution in terms of the Tresca and von Mises yield criteria. A general closed-from exact solution of burst pressure is then proposed for the three yield criteria, and the results showed that the proposed exact solution matches well with the power series solution for each yield criterion. Moreover, the von Mises flow solution is an upper bound prediction, the Tresca flow solution is a lower bound prediction, and the Zhu-Leis flow solution is an intermediate prediction that agrees well with the finite element analysis results of burst pressure for thick-walled pipes. Furthermore, two datasets of full-scale burst tests are then utilized to evaluate and validate the proposed flow solutions of burst pressure for both thin and thick-walled pipes.

42 ENGINEERING↗
DOE OSTI2021

Plastic work constrained elastoplastic topology optimization

An elastoplastic topology optimization framework for limiting plastic work generation while maximizing stiffness is presented. The kinematics and constitutive model are based on finite strain linear isotropic hardening plasticity, and the balance laws are solved using a total Lagrangian finite element formulation. Aggregation of the specific plastic work combined with an adaptive normalization scheme efficiently constrains the maximum specific plastic work. The optimization problem is regularized using an augmented partial differential equation filter, and is solved by the method of moving asymptotes where path-dependent sensitivities are derived using the adjoint method. The numerical examples show a clear dependence on the optimized maximum stiffness structures for different levels of constrained specific plastic work. It is also shown that due to the history dependency of the plasticity, the load path significantly influences the structural performance and optimized topology.

42 ENGINEERING↗
NASA NTRS1951

Analysis of plane-plastic stress problems with axial symmetry in strain-hardening range

A simple method is developed for solving plane-plastic-stress problems with axial symmetry in the strain-hardening range which is based on the deformation theory of plasticity employing the finite-strain concept. The equations defining the problems are first reduced to two simultaneous nonlinear differential equations involving two dependent variables: (a) the octahedral shear strain, and (b) a parameter indicating the ratio of principal stresses. By multiplying the load and dividing the radius by an arbitrary constant, it is possible to solve these problems without iteration for any value of the modified load. The constant is determined by the boundary condition. This method is applied to a circular membrane under pressure, a rotating disk without and with a central hole, and an infinite plate with a circular hole. Two materials, inconel x and 16-25-6, the octahedral shear stress-strain relations of which do not follow the power law, are used. Distributions of octahedral shear strain, as well as of principal stresses and strains, are obtained. These results are compared with the results of the same problems in the elastic range.

Wu, M H Lee↗
NASA NTRS1982

Finite-element analysis of initiation, stable crack growth and instability using a crack-tip-opening displacement criterion

An elastic-plastic (incremental and small strain) finite element analysis was used with a crack growth criterion to study crack initiation, stable crack growth, and instability under monotonic loading to failure of metallic materials. The crack growth criterion was a critical crack-tip-opening displacement (CTOD) at a specified distance from the crack tip, or equivalently, a critical crack-tip-opening angle (CTOA). Whenever the CTOD (or CTOA) equaled or exceeded a critical value, the crack was assumed to grow. Single values of critical CTOD were used in the analysis to model crack initiation, stable crack growth, and instability for 7075-T651 and 2024-T351 aluminum alloy compact specimens. Calculated and experimentally measured CTOD values at initiation agreed well for both aluminum alloys. These critical CTOD values were also used to predict failure loads on center-crack tension specimens and a specially-designed three-hole-crack tension specimen made of the two aluminum alloys and of 304 stainless steel. All specimens were 12.7 mm thick. Predicted failure loads for 7075-T651 aluminum alloy and 304 stainless steel specimens were generally within + or - 15 percent of experimental failure loads, whereas the predicted failure loads for 2024-T351 aluminum alloy specimens were generally within + or - 5 percent of the experimental loads.

Newman, J. C., Jr.↗
NASA NTRS1979

Finite dynamic element formulation for a plane triangular element

Higher order correction terms for the stiffness and inertia matrices associated with a triangular plane stress-strain finite dynamic element are developed in detail. Numerical results are presented which indicate that the adoption of these matrices along with a suitable quadratic matrix eigenproblem solver effects a significant economy in the free vibration solution of structure when compared with the analysis based on the usual finite element procedure. Finally, a FORTRAN IV computer program listing of the various relevant matrices is given.

Gupta, K. K.↗
NASA NTRS1986

Space Shuttle filament wound case compressive strength study. II - Analysis

In order to explain the results of the testing conducted in Part I and relate the subscale test specimen performance to case performance during prelaunch loading, a comprehensive analysis program was conducted. Finite element models of both full-scale case and subscale test coupons were made. Validated axisymmetric models were used for the case analyses. Modeling of coupons was accomplished using a generalized plane strain finite element model in most cases, although some three-dimensional finite element modeling was used to investigate the importance of three-dimensional effects on the coupon performance. The large number of plies and complexity of the test coupons used limited the finite element grids to a through-the-thickness discretization of one element per ply. Delaminations observed in the helical ply dropoff region were modeled using frictionless sliders. Comparisons of predicted strains and displacements were made with the test data from both cases and coupons to validate the finite element models. Based on observations of the test results and finite element predictions, a failure criterion for both coupons and cases was postulated and then validated by comparing performance predictions to test results. Finally, the validated failure criterion was used to predict the performance of the first two flight cases under prelaunch loads.

Messick, M. J.↗
NASA NTRS1991

A variational justification of the assumed natural strain formulation of finite elements

The objective is to study the assumed natural strain (ANS) formulation of finite elements from a variational standpoint. The study is based on two hybrid extensions of the Reissner-type functional that uses strains and displacements as independent fields. One of the forms is a genuine variational principle that contains an independent boundary traction field, whereas the other one represents a restricted variational principle. Two procedures for element level elimination of the strain field are discussed, and one of them is shown to be equivalent to the inclusion of incompatible displacement modes. Also, the 4-node C(exp 0) plate bending quadrilateral element is used to illustrate applications of this theory.

Militello, Carmelo↗
NASA NTRS1985

Thermodynamically consistent constitutive equations for nonisothermal large strain, elasto-plastic, creep behavior

The paper is concerned with the development of constitutive relations for large nonisothermal elastic-viscoplastic deformations for metals. The kinematics of elastic-plastic deformation, valid for finite strains and rotations, is presented. The resulting elastic-plastic uncoupled equations for the deformation rate combined with use of the incremental elasticity law permits a precise and purely deductive development of elastic-viscoplastic theory. It is shown that a phenomenological thermodynamic theory in which the elastic deformation and the temperature are state variables, including few internal variables, can be utilized to construct elastic-viscoplastic constitutive equations, which are appropriate for metals. The limiting case of inviscid plasticity is examined.

Riff, R.↗
DOE OSTI2024

A B-spline based gradient-enhanced micropolar implicit material point method for large localized inelastic deformations

The quasi-brittle response of cohesive-frictional materials in numerical simulations is commonly represented by softening plasticity or continuum damage models, either individually or in combination. However, classical models, particularly when coupled with non-associated plasticity, often suffer from ill-posedness and a lack of objectivity in numerical simulations. Moreover, the performance of the finite element method significantly degrades in simulations involving finite strains when mesh distortion reaches excessive levels. This represents a challenge for modeling cohesive-frictional materials, given their tendency to experience strongly localized deformations, such as those occurring during shear band dominated failure. Hence, accurate modeling of the response of cohesive-frictional solids is a demanding task. To address these challenges, we present an extension of the material point method (MPM) for the unified gradient-enhanced micropolar continuum, aiming at the analysis of finite localized inelastic deformations in cohesive-frictional materials. The generalized gradient-enhanced micropolar continuum formulation is employed to tackle challenges related to localization and softening material behavior, while the MPM addresses issues arising from excessive deformations. The method utilizes a B-spline formulation for the rigid background mesh to mitigate the well-known cell crossing errors of the MPM. To demonstrate the performance of the method, 2D and 3D numerical studies on localized failure in sandstone in plane strain compression and triaxial extension tests are presented. A comparison with finite element results confirms the suitability of the formulation. Moreover, an efficient numerical implementation of the formulation is presented, and it is demonstrated that the additional MPM specific overhead is negligible.

B-spline↗
NASA NTRS2000

Simulation of Hypervelocity Impact on Aluminum-Nextel-Kevlar Orbital Debris Shields

An improved hybrid particle-finite element method has been developed for hypervelocity impact simulation. The method combines the general contact-impact capabilities of particle codes with the true Lagrangian kinematics of large strain finite element formulations. Unlike some alternative schemes which couple Lagrangian finite element models with smooth particle hydrodynamics, the present formulation makes no use of slidelines or penalty forces. The method has been implemented in a parallel, three dimensional computer code. Simulations of three dimensional orbital debris impact problems using this parallel hybrid particle-finite element code, show good agreement with experiment and good speedup in parallel computation. The simulations included single and multi-plate shields as well as aluminum and composite shielding materials. at an impact velocity of eleven kilometers per second.

Fahrenthold, Eric P.↗