Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ELASTICITY”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
NASA NTRS2014

Elastic Model Transitions: a Hybrid Approach Utilizing Quadratic Inequality Constrained Least Squares (LSQI) and Direct Shape Mapping (DSM)

A method for transitioning linear time invariant (LTI) models in time varying simulation is proposed that utilizes both quadratically constrained least squares (LSQI) and Direct Shape Mapping (DSM) algorithms to determine physical displacements. This approach is applicable to the simulation of the elastic behavior of launch vehicles and other structures that utilize multiple LTI finite element model (FEM) derived mode sets that are propagated throughout time. The time invariant nature of the elastic data for discrete segments of the launch vehicle trajectory presents a problem of how to properly transition between models while preserving motion across the transition. In addition, energy may vary between flex models when using a truncated mode set. The LSQI-DSM algorithm can accommodate significant changes in energy between FEM models and carries elastic motion across FEM model transitions. Compared with previous approaches, the LSQI-DSM algorithm shows improvements ranging from a significant reduction to a complete removal of transients across FEM model transitions as well as maintaining elastic motion from the prior state.

Jurenko, Robert J.↗
DOE OSTI2022

A New Proof That the Number of Linear Elastic Symmetries in Two Dimensions Is Four

In this work, we present an elementary and self-contained proof that there are exactly four symmetry classes of the elasticity tensor in two dimensions: oblique, rectangular, square, and isotropic. In two dimensions, orthogonal transformations are either reflections or rotations. The proof is based on identification of constraints imposed by reflections and rotations on the elasticity tensor, and it simply employs elementary tools from trigonometry, making the proof accessible to a broad audience. For completeness, we identify the sets of transformations (rotations and reflections) for each symmetry class and report the corresponding equations of motions in classical linear elasticity.

42 ENGINEERING↗
DOE OSTI2017

Temperature dependent elastic properties of γ-phase U – 8 wt% Mo

Polycrystalline elastic moduli and stiffness tensor components of γ-phase U – 8 wt% Mo have been determined by resonant ultrasound spectroscopy in the temperature range of 25-650°C. The ambient temperature elastic properties are compared to results measured via other experimental methods and show reasonable agreement, though there is considerable variation of these properties within the literature at both the U – 8 wt% Mo composition and as a function of Mo concentration. The Young’s modulus of U – 8 wt% Mo measured in this study decreases steadily with temperature at a rate that is slower than trends previously observed at similar Mo concentrations, though the difference is not statistically significant. This first measurement of the temperature dependent elastic stiffness tensor of a polycrystalline U-Mo alloy clarifies that the behavior of the Young’s modulus is due to a strongly weakening C 11 polycrystalline stiffness tensor component, along with milder decreases in C 12 and C 44 . The unique partially auxetic properties recently predicted for singlecrystalline U-Mo are discussed in regard to their possible impact on the polycrystalline behavior of the alloy.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2021

Damage detection in a laboratory-scale wellbore applying Time Reversal and Nonlinear Elastic Wave Spectroscopy (TR NEWS)

We report Time Reversal and Nonlinear Elastic Wave Spectroscopy (TR-NEWS) has been used to focus acoustic energy, and make measurements correlated with damage in a variety of industrial materials. Most studies using TR-NEWS have focused on relatively small objects and may have multiple acoustic sources. In Earth, within energy extraction and carbon sequestration settings, the structure and scale of wellbores makes acoustic focusing challenging, as has been demonstrated through numerical simulations. This paper explores the feasibility of applying TR NEWS for damage detection in wellbores by constructing a laboratory-scale wellbore, and using TR to focus and make dynamic linear and nonlinear elastic measurements on it. During successive cycles of localized mechanical damage induced in the sample, the hysteretic nonlinear elastic parameter α increased with progressive mechanical damage. In addition to these strain-dependent changes in α, the width and peak amplitude of the time reversal focused energy changed near the damage site.

58 GEOSCIENCES↗
DOE OSTI2021

Strain Gradient Elasticity in SrTiO 3 Membranes: Bending versus Stretching

Young’s modulus determines the mechanical loads required to elastically stretch a material and also the loads required to bend it, given that bending stretches one surface while compressing the opposite one. Flexoelectric materials have the additional property of becoming electrically polarized when bent. The associated energy cost can additionally contribute to elasticity via strain gradients, particularly at small length scales where they are geometrically enhanced. Here, we present nanomechanical measurements of freely suspended SrTiO 3 crystalline membrane drumheads. In this work, we observe an unexpected nonmonotonic thickness dependence of Young’s modulus upon small deflections. Furthermore, the modulus inferred from a predominantly bending deformation is three times larger than that of a predominantly stretching deformation for membranes thinner than 20 nm. In this regime we extract a strain gradient elastic coupling of ~2.2 μN, which could be used in new operational regimes of nanoelectro-mechanics.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Elastic wave velocity evolution of shales deformed under uppermost crustal conditions

In this paper, conventional triaxial tests were performed on a series of samples of Tournemire shale along different orientations relative to bedding (0°, 90°). Experiments were carried out up to failure at increasing confining pressures ranging from 2.5 to 80 MPa, and at strain rates ranging between 3 × 10 -7 s -1 and 3 × 10 -5 s -1 . During each experiment, P and S wave elastic velocities were continuously measured along many raypaths with different orientations with respect to bedding and maximum compressive stress. This extensive velocity measurement setup allowed us to highlight the presence of plastic mechanisms such as mineral reorientation during deformation. The evolution of elastic anisotropy was quantified using Thomsen's parameters which were directly inverted from measurement of elastic wave velocity. Brittle failure was preceded by a change in P wave anisotropy, due to both crack growth and mineral reorientation. Anisotropy variations were largest for samples deformed perpendicular to bedding, at the onset of rupture. Anisotropy reversal was observed at the highest confining pressures. For samples deformed parallel to bedding, the P wave anisotropy change is weaker.

58 GEOSCIENCES↗
DOE OSTI2017

Field characterization of elastic properties across a fault zone reactivated by fluid injection

In this paper, we studied the elastic properties of a fault zone intersecting the Opalinus Clay formation at 300 m depth in the Mont Terri Underground Research Laboratory (Switzerland). Four controlled water injection experiments were performed in borehole straddle intervals set at successive locations across the fault zone. A three-component displacement sensor, which allowed capturing the borehole wall movements during injection, was used to estimate the elastic properties of representative locations across the fault zone, from the host rock to the damage zone to the fault core. Young's moduli were estimated by both an analytical approach and numerical finite difference modeling. Results show a decrease in Young's modulus from the host rock to the damage zone by a factor of 5 and from the damage zone to the fault core by a factor of 2. In the host rock, our results are in reasonable agreement with laboratory data showing a strong elastic anisotropy characterized by the direction of the plane of isotropy parallel to the laminar structure of the shale formation. In the fault zone, strong rotations of the direction of anisotropy can be observed. Finally, the plane of isotropy can be oriented either parallel to bedding (when few discontinuities are present), parallel to the direction of the main fracture family intersecting the zone, and possibly oriented parallel or perpendicular to the fractures critically oriented for shear reactivation (when repeated past rupture along this plane has created a zone).

58 GEOSCIENCES↗
DOE OSTI2018

Convergence of the multigrid reduction in time algorithm for the linear elasticity equations

This paper presents some recent advances for parallel-in-time methods applied to linear elasticity. With recent computer architecture changes leading to stagnant clock speeds, but ever increasing numbers of cores, future speedups will be available through increased concurrency. Thus, sequential algorithms, such as time stepping, will suffer a bottleneck. This paper explores multigrid reduction in time (MGRIT) for an important application area, linear elasticity. Previously, efforts at parallel-in-time for elasticity have experienced difficulties, for example, the beating phenomenon. As a result, practical parallel-in-time algorithms for this application area currently do not exist. This paper proposes some solutions made possible by MGRIT (e.g., slow temporal coarsening and FCF-relaxation) and, more importantly, a different formulation of the problem that is more amenable to parallel-in-time methods. Using a recently developed convergence theory for MGRIT and Parareal, we show that the changed formulation of the problem avoids the instability issues and allows the reduction of the error using two temporal grids. We then extend our approach to the multilevel case, where we demonstrate how slow temporal coarsening improves convergence. Finally, the paper ends with supporting numerical results showing a practical algorithm enjoying speedup benefits over the sequential algorithm.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Dependence of the Elastic Stiffness Tensors of PETN, α‐RDX, γ‐RDX, ϵ‐RDX, ϵ‐CL‐20, DAAF, FOX‐7, and β‐HMX on Hydrostatic Compression

Abstract The dependence of the components of the elastic stiffness tensors (or elastic constants) of the organic explosives PETN, RDX, CL‐20, DAAF, FOX‐7, and HMX on hydrostatic pressure up to 10 GPa have been computed using dispersion‐corrected density functional theory. We report the evolution of lattice parameters and the non‐zero stiffnesses for the tetragonal, orthorhombic, and monoclinic crystal symmetries. Linear and quadratic dependencies of the components of the elastic stiffness tensors on volumetric compression and hydrostatic pressure are tabulated for use in single crystal plasticity models.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Simple, accurate surrogate models of the elastic response of three-dimensional open truss micro-architectures with applications to multiscale topology design

Elastic meta-materials are those whose unique properties come from their micro-architecture, rather than, e.g., from their chemistry. The introduction of such architecture, which is increasingly able to be fabricated due to advances in additive manufacturing, expands the design domain and enables improved design, from the most complex multi-physics design problems to the simple compliance design problem that is our focus. Unfortunately, concurrent design of both the micro-scale and the macroscale is computationally very expensive when the former can vary spatially, particularly in three dimensions. Instead, we provide simple, accurate surrogate models of the homogenized linear elastic response of the isotruss, the octet truss, and the ORC truss based on high-fidelity continuum finite element analyses. These surrogate models are relatively accurate over the full range of relative densities, in contrast to analytical models in the literature, which we show lose accuracy as relative density increases. The surrogate models are also simple to implement, which we demonstrate by modifying Sigmund’s 99-line code to solve a three-dimensional, multiscale compliance design problem with spatially varying relative density. Here, we use this code to generate examples in both two and three dimensions that illustrate the advantage of elastic meta-materials over structures with a single length scale, i.e., those without micro-architectures.

42 ENGINEERING↗
DOE OSTI2018

Ab Initio Study of the Electronic Structure, Elastic Properties, Magnetic Feature and Thermodynamic Properties of the Ba{sub 2}NiMoO{sub 6} Material

We report first-principles calculations of the elastic properties, electronic structure and magnetic behavior performed over the Ba{sub 2}NiMoO{sub 6} double perovskite. Calculations are carried out through the full-potential linear augmented plane-wave method within the framework of the Density Functional Theory (DFT) with exchange and correlation effects in the Generalized Gradient and Local Density Approximations, including spin polarization. The elastic properties calculated are bulk modulus (B), the elastic constants (C{sub 11}, C{sub 12} and C{sub 44}), the Zener anisotropy factor (A), the isotropic shear modulus (G), the Young modulus (Y) and the Poisson ratio (υ). Structural parameters, total energies and cohesive properties of the perovskite are studied by means of minimization of internal parameters with the Murnaghan equation, where the structural parameters are in good agreement with experimental data. Furthermore, we have explored different antiferromagnetic configurations in order to describe the magnetic ground state of this compound. The pressure and temperature dependence of specific heat, thermal expansion coefficient, Debye temperature and Grüneisen parameter were calculated by DFT from the state equation using the quasi-harmonic model of Debye. A specific heat behavior C{sub V} ≈ C{sub P} was found at temperatures below T = 400 K, with Dulong–Petit limit values, which is higher than those, reported for simple perovskites.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

High-Pressure and Temperature Dependence of Electronic, Magnetic, Elastic, Thermodynamic, and Transport Properties of Full-Heusler Alloys Co{sub 2}YIn (Y = Nb, Zr)

Structural, elastic, magnetic, and electronic properties of Co{sub 2}YIn (Y = Nb, Zr) have been investigated by means of the full-potential linearized augmented plane wave along with generalized gradient approximation + U/modified Becke-Johnson method in the stable Fm-3m phase. The optimized equilibrium lattice parameter in stable phase is 6.1744 Å for Co{sub 2}NbIn and 6.2824 Å for Co{sub 2}ZrIn. We have calculated the spin-resolved density of states for both the alloys to show that the ferromagnetic behavior of Co{sub 2}NbIn is due to the exchange splitting of Co and Nb atoms and due to Co and Zr in Co{sub 2}ZrIn. In spin-up configuration, there is absence of energy gap while an energy gap can be seen in spin-down state for both the alloys exhibiting a half-metallic character. Through the elastic constants and their discussions, we have predicted the elastic behavior of the materials. The quasi-harmonic approximations have been employed to study the pressure- and temperature-dependent thermodynamic properties. The Boltzmann theory is employed to investigate the electronic transport properties of these alloys. The defined material properties may pay as a prime stand for the application of material in spintronics and thermoelectrics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Elastic Properties of Lithium and Sodium Amides

This paper presents the optimization of the structural parameters of tetragonal and orthorhombic crystals α-LiNH{sub 2} and α-NaNH{sub 2} and computation of elastic constants and flexibility coefficients using the all-electron implementation of the Gaussian and augmented plane wave density functional method. The obtained elastic constants satisfy all conditions of mechanical stability of LiNH{sub 2} and α-NaNH{sub 2} crystals. From the viewpoint of elastic properties, α-LiNH{sub 2} crystal is isotropic, while α-NaNH{sub 2} crystal is rather anisotropic. Anisotropy under compression in α-NaNH{sub 2} crystal is higher than that under shear. The former is brittle, the latter is more plastic.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Recrystallization Behavior and Super-Elasticity of a Metastable β-Type Ti-21Nb-7Mo-4Sn Alloy During Cold Rolling and Annealing

Based on the d-electron alloy design theory, a new metastable β-type titanium alloy for biomedical applications, Ti-21Nb-7Mo-4Sn (wt.%) was designed in this article. This theory predicted this alloy was a metastable single β-type titanium alloy with low elastic modulus at room temperature, and the β-phase transition (β to α″) would be easy to occur during the cold rolling. The evolution of microstructure and mechanical properties of this alloy during cold rolling plus annealing were investigated by means of x-ray diffraction, optical microscope, transmission electron microscope and mechanical properties test. The results indicate that only β phase can be identified in this alloy before cold rolling, while a large amount of lath martensite (α″ phase) appears after cold rolling due to stress-induced martensitic transformation, which is in accordance with the prediction of the d-electron alloy design theory. The recrystallization nucleation occurs preferentially in the martensite lath region during the subsequent annealing process, and this region forms a large number of nano-crystals and microcrystals, showing that cold rolling plus annealing can become a new process for refining grains. Compared with casting and cold rolling samples, the annealing materials at 923 K for 10 min have excellent comprehensive mechanical properties with lower elastic modulus (58 GPa) and higher elastic recovery rate (35.39%).

36 MATERIALS SCIENCE↗
DOE OSTI2023

Ab Initio Modeling on The Thermodynamic and Temperature-Dependent Elastic Properties of Subsystems of The FCC FeNiCoCr Medium Entropy Alloys (MEAs)

The stability, phonon spectra, thermodynamic, and temperature-dependent elasticity of subsystems of the FCC FeNiCoCr MEAs are systematically studied by the ab initio approach. Especially, the quasi-harmonic approximation (QHA) and the innovative Zentropy theory were utilized to predict the thermodynamic properties and elastic properties of FeNi, NiCo, FeNiCo, and FeNiCoCr MEAs with the consideration of magnetic transition. With the ensemble of the partition function based on the multiplicity of each magnetic configuration, the current work successfully reproduced the Curie temperature and the Schottky anomaly of heat capacity of these four MEAs purely based on the ab initio input. Meanwhile, the elastic properties of these alloys at finite temperatures are also successfully predicted with the consideration of magnetic transition. The overall results agree well with the available experimental data and CALPHAD prediction.

36 MATERIALS SCIENCE↗
DOE OSTI2016

On the consistency between nearest-neighbor peridynamic discretizations and discretized classical elasticity models

The peridynamic theory of solid mechanics is a nonlocal reformulation of the classical continuum mechanics theory. At the continuum level, it has been demonstrated that classical (local) elasticity is a special case of peridynamics. Such a connection between these theories has not been extensively explored at the discrete level. This paper investigates the consistency between nearest-neighbor discretizations of linear elastic peridynamic models and finite difference discretizations of the Navier–Cauchy equation of classical elasticity. While nearest-neighbor discretizations in peridynamics have been numerically observed to present grid-dependent crack paths or spurious microcracks, this paper focuses on a different, analytical aspect of such discretizations. We demonstrate that, even in the absence of cracks, such discretizations may be problematic unless a proper selection of weights is used. Specifically, we demonstrate that using the standard meshfree approach in peridynamics, nearest-neighbor discretizations do not reduce, in general, to discretizations of corresponding classical models. We study nodal-based quadratures for the discretization of peridynamic models, and we derive quadrature weights that result in consistency between nearest-neighbor discretizations of peridynamic models and discretized classical models. The quadrature weights that lead to such consistency are, however, model-/discretization-dependent. We motivate the choice of those quadrature weights through a quadratic approximation of displacement fields. The stability of nearest-neighbor peridynamic schemes is demonstrated through a Fourier mode analysis. Finally, an approach based on a normalization of peridynamic constitutive constants at the discrete level is explored. This approach results in the desired consistency for one-dimensional models, but does not work in higher dimensions. The results of the work presented in this paper suggest that even though nearest-neighbor discretizations should be avoided in peridynamic simulations involving cracks, such discretizations are viable, for example for verification or validation purposes, in problems characterized by smooth deformations. Furthermore, we demonstrate that better quadrature rules in peridynamics can be obtained based on the functional form of solutions.

42 ENGINEERING↗
DOE OSTI2018

Unconditionally stable, second-order schemes for gradient-regularized, non-convex, finite-strain elasticity modeling martensitic phase transformations

In the setting of continuum elasticity martensitic phase transformations are characterized by a non-convex free energy density function that possesses multiple wells in strain space and includes higher-order gradient terms for regularization. Metastable martensitic microstructures, defined as solutions that are local minimizers of the total free energy, are of interest and are obtained as steady state solutions to the resulting transient formulation of Toupin’s gradient elasticity at finite strain. This type of problem poses several numerical challenges including stiffness, the need for fine discretization to resolve microstructures, and following solution branches. Accurate time-integration schemes are essential to obtain meaningful solutions at reasonable computational cost. Here in this work we introduce two classes of unconditionally stable second-order time-integration schemes for gradient elasticity, each having relative advantages over the other. Numerical examples are shown highlighting these features.

42 ENGINEERING↗
DOE OSTI2017

Elastic knowledge base of bcc Ti alloys from first-principles calculations and CALPHAD-based modeling

Here, titanium alloys are being investigated as suitable materials for load-bearing implants because of their biocompatibility and mechanical properties. Stress shielding, a common issue with the current load-bearing implant materials, occurs due to a Young’s modulus (E) mismatch between bone (~10–40 GPa) and implants (such as Ti-6Al-4V ~ 110 GPa), which leads to bone dying around the implant and ultimately implant failure. Reducing the Young’s modulus of Ti alloys may overcome the issues of stress shielding and improve implant materials. In the present work, first-principles calculations have been used to predict the single crystal elastic stiffness coefficients (c ij ’s) for the Ti-containing ternary alloys Ti-X-Y(X ≠ Y = Mo, Nb, Sn, Ta, Zr) in the bcc lattice. It is found that the ternary Ti-X-Y (X ≠ Y = Mo, Nb, Ta) alloys behave similarly; so do the ternary Ti-X-Sn (X = Mo, Nb, Ta) alloys and the Ti-X-Zr (X = Mo, Nb, Ta) alloys. This is expected due to the similarity between the Mo, Nb and Ta elements. The results also show that the Ti-Zr-X alloys stabilized the bcc phase at lower alloying concentrations. The polycrystalline aggregate properties are also estimated from the c ij ’s, including bulk modulus, shear modulus and Young’s modulus. The results show that Ti-alloys with compositions close to the bcc stability limit have the lowest E. In combination with previous predictions, a complete elastic database has been established using theCALPHAD (CALculation of PHAse Diagram) based modeling approach. The database results are compared with the E of higher order Ti alloys and shown to be able to predict the E accurately. This complete data-base forms a foundation to tailor Ti alloys for desired elastic properties.

36 MATERIALS SCIENCE↗