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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Sound velocities and thermal equation of state of fcc -iron-nickel alloys at high pressure and high temperature: Implications for the cores of Moon and several planets

Fcc-Fe-Ni alloy is believed to be the most dominant solid constitute of moderate-sized terrestrial planetary cores. Investigating the physical properties, especially the density and sound velocity of Fe-Ni alloys and comparing them with seismic observations is an indispensable approach to constructing compositional models for planetary interiors. In this study, we conducted sound velocity measurements on Fe-Ni alloys with 10 wt.% and 20 wt.% Ni up to ∼13.5 GPa and 1073 K, using the ultrasonic interferometry technique in a multi-anvil apparatus in conjunction with synchrotron radiation. By fitting the experimental data to finite strain equations, the bulk and shear moduli and their pressure and temperature derivatives are derived, yielding K S0 =145.8(14) GPa, G 0 = 73.2(7) GPa, K S0 ’ = 5.89(24), G 0 ’ = 2.89(8), (∂K S /∂T) P = -0.0181(12) GPa/K and (∂G/∂T) P = -0.0393(10) GPa/K for fcc-Fe 80 Ni 20 . An examination of the density-velocity relationship shows that compressional wave velocity is insensitive to temperature within the current pressure and temperature range, while shear wave velocity exhibits a large reduction with increasing temperature. Here, extrapolation of the sound velocities following the finite strain theories suggests that much slower Vs should be expected at pressure and temperature conditions corresponding to those of the lunar core. Possible core density and velocity profiles for other moderate planets and satellites, such as Mars, Mercury, and Ganymede are also calculated.

Equation of state↗
DOE OSTI2020

The Elastic Properties of β-Mg 2 SiO 4 Containing 0.73 wt.% of H 2 O to 10 GPa and 600 K by Ultrasonic Interferometry with Synchrotron X-Radiation

Here, we measured the elastic velocities of a synthetic polycrystalline β-Mg 2 SiO 4 containing 0.73 wt.% H 2 O to 10 GPa and 600 K using ultrasonic interferometry combined with synchrotron X-radiation. Third-order Eulerian finite strain analysis of the high P and T data set yielded K so = 161.5(2) GPa, G o = 101.6(1) GPa, and (∂K s /∂P) T = 4.84(4), (∂G/∂P) T = 1.68(2) indistinguishable from K so = 161.1(3) GPa, G o = 101.4(1) GPa, and (∂K s /∂P) T = 4.93(4), (∂G/∂P) T = 1.73(2) from the linear fit. The hydration of the wadsleyite by 0.73 wt.% decreases K s and G moduli by 5.3% and 8.6%, respectively, but no measurable effect was noted for (∂K s /∂P) T and (∂G/∂P) T . The temperature derivatives of the K s and G moduli from the finite strain analysis (∂K S /∂T) P = –0.013(2) GPaK –1 , (∂G/∂T) P = –0.015(0.4) GPaK –1 , and the linear fit (∂K S /∂T) P = –0.015(1) GPaK –1 , (∂G/∂T) P = –0.016(1) GPaK –1 are in agreement, and both data sets indicating the |(∂G/∂T) P | to be greater than |(∂K S /∂T) P |. Calculations yield ΔV p(α-β) = 9.88% and ΔV S(α-β) = 8.70% for the hydrous β-Mg 2 SiO 4 and hydrous α-Mg 2 SiO 4 , implying 46–52% olivine volume content in the Earth’s mantle to satisfy the seismic velocity contrast ΔV s = ΔV P = 4.6% at the 410 km depth.

58 GEOSCIENCES↗
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↗
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↗
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↗
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↗
DOE OSTI2022

Marmot V2

MARMOT is a robust numerical tool for mesoscale modeling of fuel performance developed under the NEAMS Fuels technical area to predict the coevolution of microstructure and properties in fuel and cladding materials. MARMOT accomplishes this using the phase field method coupled with finite strain mechanics and heat conduction. MARMOT is based on the open source Multiphysics Object-Oriented Simulation Environment (MOOSE) and solves the coupled partial differential equations defining the physics using the finite element method. MARMOT is being developed in order to facilitate the development of improved materials models for fuel performance, but it is also being developed as a powerful tool in and of itself for the simulation of mesoscale fuel performance.

Aagesen, LarryK.↗
DOE OSTI2025

Modeling low cycle fatigue (LCF) of additively manufactured Hastelloy X using An accelerated crystal plasticity fatigue damage model

This paper presents a microstructure-based model for low cycle fatigue (LCF) behavior and life of Nickel-based alloy Hastelloy X manufactured using laser-powder bed fusion (L-PBF) additive manufacturing (AM). AM Hastelloy X, a solution-strengthened alloy, is tested at elevated temperature under fully reversed LCF conditions at different strain levels. A generalized plane strain finite element model is generated from electron backscatter diffraction (EBSD) characterization. The constitutive behavior of the material under fatigue is modeled using crystal plasticity and calibrated with both monotonic tensile and cyclic stress–strain data. The fatigue micro-crack initiation and propagation in the microstructure is modeled using a modified Chaboche fatigue damage model. An embedded boundary condition with a homogenous medium is used to apply the cyclic deformation and prevent numerically introduced over-constraints during fatigue simulation. A ‘cycle-jump’ method is used to accelerate the fatigue simulation and reduce the computational cost. The simulation results are compared to LCF experiments, showing satisfactory matches in cyclic stress behavior and number of cycles to macro-crack initiation for all applied strain ranges. In addition, the model illustrates the potential for quantifying microscale fatigue life impacting factors such as microstructure and surface roughness, which is needed to accurately quantify the reliability of AM components in service.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Calibrating constitutive models with full‐field data via physics informed neural networks

Abstract The calibration of solid constitutive models with full‐field experimental data is a long‐standing challenge, especially in materials that undergo large deformations. In this paper, we propose a physics‐informed deep‐learning framework for the discovery of hyperelastic constitutive model parameterizations given full‐field surface displacement data and global force‐displacement data. Contrary to the majority of recent literature in this field, we work with the weak form of the governing equations rather than the strong form to impose physical constraints upon the neural network predictions. The approach presented in this paper is computationally efficient, suitable for irregular geometric domains, and readily ingests displacement data without the need for interpolation onto a computational grid. A selection of canonical hyperelastic material models suitable for different material classes is considered including the Neo–Hookean, Gent, and Blatz–Ko constitutive models as exemplars for general non‐linear elastic behaviour, elastomer behaviour with finite strain lock‐up, and compressible foam behaviour, respectively. We demonstrate that physics informed machine learning is an enabling technology and may shift the paradigm of how full‐field experimental data are utilized to calibrate constitutive models under finite deformations.

Hamel, Craig M.↗
DOE OSTI2025

A fully-integrated lattice Boltzmann method for fluid–structure interaction

Here we present a fully-integrated lattice Boltzmann (LB) method for fluid–structure interaction (FSI) simulations that efficiently models deformable solids in complex suspensions and active systems. Our Eulerian method (LBRMT) couples finite-strain solids to the LB fluid on the same fixed computational grid with the reference map technique (RMT). An integral part of the LBRMT is a new LB boundary condition for moving deformable interfaces across different densities. With this fully Eulerian solid–fluid coupling, the LBRMT is well-suited for parallelization and simulating multi-body contact without remeshing or extra meshes. We validate its accuracy via a benchmark of a deformable solid in a lid-driven cavity, then showcase its versatility through examples of soft solids rotating and settling. The LBRMT achieves a spatial convergence rate between first-order and second-order for FSI simulations and is designed for low to intermediate Reynolds number flows with finite inertia at small Mach numbers. With simulations of complex suspensions mixing, we highlight the potential of the LBRMT for studying collective behavior in soft matter and biofluid dynamics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Response tailoring of elasto-plastic trusses

Abstract In this work we tailor the response of trusses loaded beyond their yield limit. The truss structures are modeled using finite strain theory and rate-independent elasto-plasticity. We design trusses with a tailored mechanical response that is between “elastic” and “elastic-ideal-plastic” subject to the volume constraint. The design updates are generated by the gradient-based Method of Moving Asymptotes (MMA) solver and the sensitivities of the response functions are computed using a path-dependent adjoint sensitivity analysis. The computations are performed in Matlab.

42 ENGINEERING↗
DOE OSTI2020

Inelastic equation of state for solids

In this study, a complete inelastic equation of state (IEOS) for solids is developed based on a superposition of thermodynamic energy potentials. The IEOS allows for a tensorial stress state by including an isochoric hyperelastic Helmholtz potential in addition to the zero-kelvin isotherm and lattice vibration energy contributions. Inelasticity is introduced through the nonlinear equations of finite strain plasticity which utilize the temperature dependent Johnson–Cook yield model. Material failure is incorporated into the model by a coupling of the damage history variable to the energy potentials. The numerical evaluation of the IEOS requires a nonlinear solution of stress, temperature and history variables associated with elastic trial states for stress and temperature. The model is implemented into the ALEGRA shock and multi-physics code and the applications presented include single element deformation paths, the Taylor anvil problem and an energetically driven thermo-mechanical problem.

42 ENGINEERING↗
DOE OSTI2021

Simulations of crack extensions in arc-shaped tension specimens of uncharged and hydrogen-charged 21-6-9 austenitic stainless steels using cohesive zone modeling with varying cohesive parameters

Crack extensions in uncharged and hydrogen-charged side-grooved A(T) specimens of conventionally forged 21-6-9 austenitic stainless steels are simulated using the cohesive zone modeling approach. Two-dimensional plane strain finite element analyses with fixed cohesive parameters are first conducted to fit the experimental load-displacement curves. Similar analyses using varying cohesive parameters as functions of the crack extension are then conducted to fit the experimental load-crack extension and crack extension-displacement curves. The computational results with varying cohesive parameters can fit very well the experimental data. The computational results also indicate that the average cohesive energy for the hydrogen-charged A(T) specimen is lower than that for the uncharged A(T) specimen.

42 ENGINEERING↗
DOE OSTI2021

Anomalous elastic behavior of tantalum at high pressures: Experimental and theoretical studies

Elastic wave velocities of polycrystalline tantalum have been measured up to 13.6 GPa at room temperature using ultrasonic interferometry technique in a multi-anvil apparatus. The bulk and shear moduli, as well as their pressure derivatives are obtained using Eulerian 3rd order finite strain equations, yielding K S0 = 193.9(29) GPa, G 0 = 69.3(10) GPa, K S0' = 3.18(5) and G 0' = 0.88(1). Here the Steinberg-Guinan yield stress model is examined using the current shear modulus and its pressure derivative, where a slightly lower A 0 = 0.0127 GPa –1 compared to previous studies is obtained. First-principles calculations on the elasticity of tantalum have also been performed up to 200 GPa. Anomalous softening of the elastic shear constant C 44 and shear wave velocity V S is observed between 80– 180 GPa, whereas the anisotropy shows a softening-stiffening behavior with increasing pressure and reaches a minimum at around 150 GPa. Such anomalies might commonly originate from the nesting of Fermi surface for Group VB metals.

36 MATERIALS SCIENCE↗