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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

High‐Pressure and High‐Temperature Single‐Crystal Elasticity of Cr‐Pyrope: Implications for the Density and Seismic Velocity of Subcontinental Lithospheric Mantle

Abstract Single‐crystal X‐ray diffraction and Brillouin spectroscopy experiments were performed on a natural Cr‐pyrope (Prp 71.0 Alm 12.6 Sps 0.7 Grs 3.5 Uvr 12.2 ) at high pressure and high temperature up to 11.0 GPa and 800 K. Fitting the collected data to the third‐order finite strain equation yields bulk modulus ( K S 0 ), shear modulus ( G 0 ), their pressure ((∂ K S /∂ P ) T and (∂ G /∂ P ) T ) and temperature (( ∂K S /∂T ) P and ( ∂G/∂T ) P ) derivatives, K S 0 = 167.7(8) GPa, G 0 = 91.5(5) GPa, (∂ K S /∂ P ) T = 4.3(1), (∂ G /∂ P ) T = 1.4(1), ( ∂K S /∂T ) P = − 0.0175(1) GPa/K and ( ∂G/∂T ) P = − 0.0073(1) GPa/K. Using the obtained results, we examined whether the elastic properties of the Cr‐pyrope can be accurately calculated from those of endmembers including pyrope, almandine, grossular, and uvarovite assuming a linear relationship between elastic properties and composition (end‐member model). The results indicate that the end‐member model provides a sufficient approximation for the elastic properties of Cr‐pyrope in calculating the density and velocity of the subcontinental lithospheric mantle (SCLM). We modeled the densities and velocities of three typical types of SCLM (Archon, Proton, and Tecton) in order to investigate how the variation of chemical composition influences the SCLM. We obtained that the compositional change from the Archon to the Tecton increases the density of the SCLM significantly, which can be an important prerequisite for SCLM delamination. However, the compositional variation only slightly changes the velocity of the SCLM and the change is within the uncertainty of the calculation. Moreover, in comparison to the velocity, ρ / V P and ρ / V S are much more sensitive to the compositional change of the SCLM.

Xu, Jingui↗
DOE OSTI2024

Sound Velocities of Stishovite at Simultaneous High Pressure and High Temperature Suggest an Eclogite‐Rich Layer Beneath the Hawaii Hotspot

Compressional and shear wave velocities of polycrystalline stishovite (SiO 2 ) have been measured at simultaneous high pressures and temperatures up to 14.5 GPa and 800°C. By fitting velocities to the finite strain equations, the elastic moduli and density were determined to be K S0 = 306.6(46) GPa, $K$$^{′}_{S}$ = 4.92(10), ∂K S /∂T = −0.024(1) GPa/K, G 0 = 229.0(34) GPa, G′ = 1.07(10), ∂G/∂T = −0.017(1) GPa/K, ρ 0 = 4.287(2) g/cm 3 . Our modeling suggested that, in the eclogite, coesite-stishovite transition can increase P and S wave velocities by 2.4% and 3.5%, respectively. A comparison between geophysical observations and our model shows that the coesite-stishovite phase transition in the eclogite can potentially be responsible for the occurrence of the X discontinuity beneath Hawaii. In addition, our current results suggest an eclogite-rich layer between 340 and 450 km depth beneath Hawaii. The eclogite concentration at the top and bottom of the layer is 41–55 vol% and >77 vol%, respectively.

36 MATERIALS SCIENCE↗
DOE OSTI2024

In Situ Determination of Thermoelastic Properties of Magnesite at High Pressure and Temperature With Implications to Seismic Detectability of Moderately Carbonated Lithologies in the Earth's Mantle

Abstract Subduction of carbonate‐bearing oceanic plates into Earth's interior recycles carbon from the surface to the deep mantle. The subducted carbonates can significantly affect mantle properties and dynamics. Magnesite is recognized as one of the major potential carbon hosts in the deep mantle because of its stability up to deep lower mantle conditions. However, despite many previous studies on the equation of state and elastic properties of magnesite, large discrepancies still exist for its elastic moduli and their pressure and temperature derivatives. Here we report in situ density and elastic wave velocity measurements on a natural magnesite at simultaneous high pressure‐temperature conditions up to ∼8 GPa–1073 K in a multi‐anvil apparatus using ultrasonic and synchrotron X‐ray techniques. Global fitting of the data set to finite strain equations yieldsK S0 = 114.0 ± 1.2 GPa, = 4.9 ± 0.3, = (−0.019 ± 0.002) GPa/K,G 0 = 68.6 ± 0.4 GPa,G′ = 1.6 ± 0.1, and = (−0.018 ± 0.001) GPa/K. Compared to other major upper mantle and transition zone minerals, magnesite has intermediateV P , the lowestV S , and the lowest density. Thus, magnesite possesses higherV P /V S ratio than other mantle minerals in normal mantle regions, whereas this feature is less pronounced in subduction zone environments. Modeling of the velocity profiles of carbonated lithologies along different geotherms suggests that moderately‐enriched magnesite domains are unlikely to be detected seismically in the Earth's mantle.

Geochemistry & Geophysics↗
DOE OSTI2025

Sound Velocity of Hematite up to 13 GPa

We performed acoustic measurements on polycrystalline hematite using ultrasonic interferometric techniques in a Kawai-type multi-anvil apparatus. The elastic moduli and their pressure derivatives were determined by fitting the experimental data to third-order finite strain equations, yielding: 𝐾 𝑆,0 = 214.8 (20) GPa, 𝐺 0 = 92.0 (2) GPa, 𝐾$^{′}_{S,0}$ = 5.70 (1), and 𝐺$^{′}_{0}$ = 0.88 (4). Modeling of the velocity profiles of hematite along hot and cold slab geotherms shows that compressional (V P ) and shear (V S ) wave velocities are 22%–24% and 53%–55% lower, respectively, compared to the Preliminary Reference Earth Model at the top of the lower mantle. Furthermore, these findings likely offer a potential explanation for the low-velocity anomalies associated with slab graveyards at the base of the lower mantle.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Single Crystal Plasticity Model with Deformation Twinning for the High Rate Deformation of β-HMX

Deformation twinning is an active deformation mechanism during loading along certain crystallographic directions of β-cyclotetramethylene tetranitramine (β-HMX). In this work, a finite strain thermomechanical model is extended to include twinning as a deformation mechanism in addition to plastic slip. The stress is derived from the free energy expression including a term representing the equation of state. The crystal plasticity framework is used to divide the total strain into inelastic and elastic component, where the elastic part is used in the expression for the free energy. The twin systems are treated as pseudo slip systems and the shear rate on the twin systems is evaluated in terms of the twin resistance and the projection of the stress tensor. The model parameters were calibrated against a set of plate impact experiments performed on β-HMX by Dick et al. The remaining plate impact experiments are used to evaluate the predictive capability of the model. The quality of the model fits and predictions is discussed from a physical and modeling perspective. Particularly, the role of twin modeling on results is highlighted. Finally, the model is used to explore the relationship between the propensity for twinning and crystal orientation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Thermoelastic properties of tungsten at simultaneous high pressure and temperature

In this work, the compressional ( P ) and shear wave velocities ( S ) and unit cell volumes (densities) of polycrystalline tungsten (W) have been measured simultaneously up to 10.5 GPa and 1073 K using ultrasonic interferometry in conjunction with x-ray diffraction and x-radiography techniques. Thermoelastic properties of W were derived using different methods. We obtained the isothermal bulk modulus K T 0 = 310.3(1.5) GPa, its pressure derivative K ' T 0 = 4.4(3), its temperature derivative at constant pressure ( ∂ K T / ∂ T ) P = – 0.0138 ( 1 ) GPa K – 1 and at constant volume ( ∂ K T / ∂ T ) V = – 0.0050 GPa K – 1 , the thermal expansion α ( 0, T ) = 1.02(27) × 10 –5 + 7.39(3.2) × 10 –9 T (K –1 ), as well as the pressure derivative of thermal expansion ( ∂ α / ∂ P ) T = – 1.44 ( 1 ) × 10 – 7 K – 1 GP a – 1 based on the high-temperature Birch–Murnaghan equation of state (EOS), the Vinet EOS, and thermal pressure approach. Finite strain analysis allowed us to derive the elastic properties and their pressure/temperature derivatives independent of the choice of pressure scale. A least-squares fitting yielded K S 0 = 314.5(2.5) GPa, K S 0 ' = 4.45(9), (∂ K S /∂ T ) P = – 0.0076(6) GPa K –1 , G 0 = 162.4(9) GPa, G 0 ' = 1.8(1), (∂ G /∂ T ) P = – 0.0175(9) GPa K –1 , and α 298 K = 1.23 × 10 – 5 K – 1 . Fitting current data to the Mie–Grüneisen–Debye EOS with derived θ 0 = 383.4 K yielded γ 0 = 1.81 ( 6 ) and q = 0.3 . The thermoelastic parameters obtained from various approaches are consistent with one another and comparable with previous results within uncertainties. Our current study provides a complete and self-consistent dataset for the thermoelastic properties of tungsten at high P – T conditions, which is important to improve the theoretical modeling of these materials under dynamic conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Elastic constants in monocrystalline tungsten under quasi-hydrostatic pressures to 11.3 GPa

Compressional (P) and shear (S) wave velocities of tungsten single crystals along the [100] and [110] directions were measured using ultrasonic interferometry at room temperature up to 11.3 GPa. Least-squares fitting of $V$$^{[100]}_{P}$, $V$$^{[100]}_{S}$, $V$$^{[110]}_{P}$, and pressure to finite strain (FS) equations yields the elastic constants: C 11 = 523.5(5) GPa, C 12 = 205.1(32) GPa, and C 44 = 160.8(4) GPa, along with their respective pressure derivatives: $C$$^{′}_{11}$ = 6.20(2), $C$$^{′}_{12}$ = 3.35(1), and $C$$^{′}_{44}$ = 1.65(6). Using the Voigt–Reuss–Hill approximation, the elastic moduli were derived as K S0 = 311.2(22) GPa and G 0 = 160.2(7)GPa, along with their respective pressure derivatives: $K$$^{′}_{S0}$ = 4.30(4) and $G$$^{′}_{0}$ = 1.56(1). The Debye temperature of tungsten was determined to be 380.7(8) K, showing good agreement with previous calorimetric measurements. The elastic anisotropy increases slightly from 1.01(1) at ambient pressure to 1.03(1) at 11.3 GPa, suggesting that tungsten remains nearly isotropic under compression. Poisson's ratio slightly increased from 0.281(3) to 0.288(3) with pressure. Additionally, Pugh's ratio decreased from 0.512(4) to 0.494(4), while Pettifor's ratio increased from 0.148(10) to 0.175(10) as pressure reached 11.3 GPa. These results suggest that tungsten is weakly ductile compared to the critical thresholds for ductile behavior (i.e., Pugh's ratio <0.6 and Pettifor's ratio >0, indicating ductile behavior), but that its ductility increases marginally under pressure.

Wang, Ran [Stony Brook University, NY (United Stat↗
DOE OSTI2023

Single-crystal elasticity of (Al,Fe)-bearing bridgmanite up to 82 GPa

Abstract Thermoelastic properties of mantle candidate minerals are essential to our understanding of geophysical phenomena, geochemistry, and geodynamic evolutions of the silicate Earth. However, the lower-mantle mineralogy remains much debated due to the lack of single-crystal elastic moduli (Cij) and aggregate sound velocities of (Al,Fe)-bearing bridgmanite, the most abundant mineral of the planet, at the lower mantle pressure-temperature (P-T) conditions. Here we report single-crystal Cij of (Al,Fe)-bearing bridgmanite, Mg0.88Fe0.1Al0.14Si0.90O3 (Fe10-Al14-Bgm) with Fe3+/ΣFe = ~0.65, up to ~82 GPa using X-ray diffraction (XRD), Brillouin light scattering (BLS), and impulsive stimulated light scattering (ISLS) measurements in diamond-anvil cells (DACs). Two crystal platelets with orientations of (–0.50, 0.05, –0.86) and (0.65, –0.59, 0.48), that are sensitive to deriving all nine Cij, are used for compressional and shear wave velocity (νP and νS) measurements as a function of azimuthal angles over 200° at each experimental pressure. Our results show that all Cij of singe-crystal Fe10-Al14-Bgm increase monotonically with pressure with small uncertainties of 1–2% (±1σ), except C55 and C23, which have uncertainties of 3–4%. Using the third-order Eulerian finite-strain equations to model the elasticity data yields the aggregate adiabatic bulk and shear moduli and respective pressure derivatives at the reference pressure of 25 GPa: KS = 326 ± 4 GPa, µ = 211 ± 2 GPa, KS′ = 3.32 ± 0.04, and µ′ = 1.66 ± 0.02 GPa. The high-pressure aggregate νS and νP of Fe10-Al14-Bgm are 2.6–3.5% and 3.1–4.7% lower than those of MgSiO3 bridgmanite end-member, respectively. These data are used with literature reports on bridgmanite with different Fe and Al contents to quantitatively evaluate pressure and compositional effects on their elastic properties. Comparing with one-dimensional seismic profiles, our modeled velocity profiles of major lower-mantle mineral assemblages at relevant P-T suggest that the lower mantle could likely consist of about 89 vol% (Al,Fe)-bearing bridgmanite. After considering uncertainties, our best-fit model is still indistinguishable from pyrolitic or chondritic models.

Geochemistry & Geophysics↗
DOE OSTI2025

Sound Velocities in Vanadium Reveal Complex Elastic Behavior at High Pressures

Compressional (VP) and shear (VS) wave velocities of polycrystalline vanadium were measured simultaneously up to 11.5 GPa at room temperature using ultrasonic interferometry in a multi-anvil press. Complex softening behavior in VS and resulting shear moduli are discovered, possibly revealing a precursor to the reported phase transition within 30–60 GPa. The current data enables a comprehensive assessment of the elastic and mechanical properties of vanadium at high pressures, including bulk and shear moduli, Young’s modulus, Poisson’s ratio, and Pugh’s ratio. Through fitting to the 3rd-order finite strain equations, the elastic moduli and their pressure derivatives were determined to be K S0 = 151 (2) GPa, G 0 = 46.9 (8) GPa, K$^{’}_{S0}$ = 3.47 (5), and G$^{’}_{0}$ = 0.62 (1). These experimental results allow us to compare with and benchmark the existing Steinberg–Guinan models for extrapolations to extreme pressure and temperature conditions.

36 MATERIALS SCIENCE↗
NASA NTRS1971

Theoretical equations of state

Theoretical equations of state in geophysics, considering systematics approach to laboratory data, seismic velocity profiles, finite strain and atomistic approach

Sammis, C. G.↗
NASA NTRS1972

Wave-front singularities for two-dimensional anisotropic elastic waves.

Wavefront singularities for the displacement functions, associated with the radiation of linear elastic waves from a point source embedded in a finitely strained two-dimensional elastic solid, are examined in detail. It is found that generally the singularities are of order d to the -1/2 power, where d measures distance away from the front. However, in certain exceptional cases singularities of order d to the -n power, where n = 1/4, 2/3, 3/4, may be encountered.

Payton, R. G.↗
NASA NTRS1974

On numerically accurate finite element solutions in the fully plastic range

A general criterion for testing a mesh with topologically similar repeat units is given, and it is shown that only a few conventional element types and arrangements are suitable for computations in the fully plastic range. Further, a new variational principle, which can easily be incorporated into an existing finite element program, is presented. This makes it possible to perform accurate computations even for element designs that would not normally be suitable. An appendix extends the discussion to elastic-plastic computation at finite strain.

Nagtegaal, J. C.↗
NASA NTRS1977

The relationship between crustal tectonics and internal evolution in the moon and Mercury

The relationship between crustal tectonics and thermal evolution is discussed in terms of the moon and Mercury. Finite strain theory and depth and temperature-dependent thermal expansion are used to evaluate previous conclusions about early lunar history. Factors bringing about core differentiation in the first 0.6 b.y. of Mercurian evolution are described. The influence of concentrating radioactive heat sources located in Mercury's crust on the predicted contraction is outlined. The predicted planetary volume change is explored with regard to quantitative limits on the extent of Mercurian core solidification. Lunar and Mercurian thermal stresses involved in thermal evolution are reviewed, noting the history of surface volcanism. It is concluded that surface faulting and volcanism are closely associated with the thermal evolution of the whole planetary volume. As the planet cools or is heated, several types of tectonic and volcanic effects may be produced by thermal stress occurring in the lithosphere.

Solomon, S. C.↗
NASA NTRS1979

Energy absorption mechanisms during crack propagation in metal matrix composites

The stress distributions around individual fibers in a unidirectional boron/aluminum composite material subjected to axial and transverse loadings are being studied utilizing a generalized plane strain finite element analysis. This micromechanics analysis was modified to permit the analysis of longitudinal sections, and also to incorporate crack initiation and propagation. The analysis fully models the elastoplastic response of the aluminum matrix, as well as temperature dependent material properties and thermal stress effects. The micromechanics analysis modifications are described, and numerical results are given for both longitudinal and transverse models loaded into the inelastic range, to first failure. Included are initially cracked fiber models.

Murphy, D. P.↗
NASA NTRS1981

Analysis of crack propagation as an energy absorption mechanism in metal matrix composites

The crack initiation and crack propagation capability was extended to the previously developed generalized plane strain, finite element micromechanics analysis. Also, an axisymmetric analysis was developed, which contains all of the general features of the plane analysis, including elastoplastic material behavior, temperature-dependent material properties, and crack propagation. These analyses were used to generate various example problems demonstrating the inelastic response of, and crack initiation and propagation in, a boron/aluminum composite.

Adams, D. F.↗
NASA NTRS1983

Characteristics of thermally-induced transverse cracks in graphite epoxy composite laminates

The characteristics of thermally induced transverse cracks in T300/5208 graphite-epoxy cross-ply and quasi-isotropic laminates were investigated both experimentally and analytically. The formation of transverse cracks and the subsequent crack spacing present during cool down to -250 F (116K) and thermal cycling between 250 and -250 F (116 and 394K) was investigated. The state of stress in the vicinity of a transverse crack and the influence of transverse cracking on the laminate coefficient of thermal expansion (CTE) was predicted using a generalized plane strain finite element analysis and a modified shear lag analysis. A majority of the cross-ply laminates experienced transverse cracking during the initial cool down to -250 F whereas the quasi-isotropic laminates remained uncracked. The in situ transverse strength of the 90 degree layers was more than 1.9 times greater than the transverse strength of the unidirectional 90 degree material for all laminates investigated.

Adams, D. S.↗