Search NASA⌕ Search

SEARCH · Search NASA

Results for “material elastic constants”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Superior High-Temperature Strength in a Supersaturated Refractory High-Entropy Alloy

Refractory high-entropy alloys (RHEA) show promising applications at high temperatures. However, achieving high strengths at elevated temperatures above 1,173K is still challenging due to heat softening. Using intrinsic material characteristics as the alloy-design principles, a single-phase body-centered-cubic (BCC) CrMoNbV RHEA with high-temperature strengths (beyond 1,000 MPa at 1,273 K) is designed, superior to other reported RHEAs as well as conventional superalloys. The origin of the high-temperature strength is revealed by in-situ neutron scattering, transmission-electron microscopy, and first-principles calculations. The CrMoNbV’s elevated-temperature strength retention up to 1,273 K arises from its large atomic-size and elastic-modulus mismatches, the insensitive temperature dependence of elastic constants, and the dominance of non-screw character dislocations caused by the strong solute pinning, which makes the solid-solution strengthening pronounced. The alloy-design principles and the insights in this study pave the way to design RHEAs with outstanding high-temperature strength.

36 MATERIALS SCIENCE↗
NASA NTRS1988

Temperature effect on stress concentration around circular hole in a composite material specimen representative of X-29A forward-swept wing aircraft

The theory of anisotropic elasticity was used to evaluate the anisotropic stress concentration factors of a composite laminated plate containing a small circular hole. This advanced composite was used to manufacture the X-29A forward-swept wing. It was found for composite material, that the anisotropic stress concentration is no longer a constant, and that the locations of maximum tangential stress points could shift by changing the fiber orientation with respect to the loading axis. The analysis showed that through the lamination process, the stress concentration factor could be reduced drastically, and therefore the structural performance could be improved. Both the mixture rule approach and the constant strain approach were used to calculate the stress concentration factor of room temperature. The results predicted by the mixture rule approach were about twenty percent deviate from the experimental data. However, the results predicted by the constant strain approach matched the testing data very well. This showed the importance of the inplane shear effect on the evaluation of the stress concentration factor for the X-29A composite plate.

Yeh, Hsien-Yang↗
DOE OSTI2025

DuctGPT: A Generative Transformer for Forward Screening of Ductile Refractory Multi-Principal Element Alloys

Designing ductile materials for extreme environments such as fusion reactors requires a deep understanding of the complex interplay between electronic structure, mechanical stability, and wide compositional space. Here, in this work, we introduce DuctGPT, a physics-informed, GPT-powered machine learning platform that enables rapid and accurate prediction of ductility across a wide range of refractory multi-principal element alloys (MPEAs). Trained on both experimental and high-fidelity computational data, DuctGPT integrates descriptors such as density of states at the Fermi level, elastic constants, and valence electron concentration to capture the fundamental mechanisms governing ductile versus brittle behavior. Using this framework, we screen over 1000 compositions in of body-centered cubic (BCC) MPEAs, including two new alloy classes, i.e., NbTa-rich (NbTa $>$ 50 at.%) NbTa-Ti-V and W-rich ($>$ 50 at.%) W-Ti-V MPEAs, to rapidly identify promising alloy compositions with enhanced ductility. Validation against experimental data confirms the model's ability to predict ductility with high fidelity and low uncertainty. By leveraging conversational AI and robust physical modeling, DuctGPT provides a blueprint for the next generation of alloy design assistants, enabling human-AI collaboration in the accelerated discovery of ductile, high-performance materials for fusion, aerospace, and advanced manufacturing.

AI/ML↗
DOE OSTI2022

Symmetry-correct bonding in density functional theory calculations for delta phase Pu

The long-held conclusion that magnetic order in delta phase Pu makes the structure mechanically unstable in density functional theory calculations is incorrect if the magnetic order is allowed to be non-collinear. The non-collinear 3Q spin structure is intrinsically cubic and makes all structurally equivalent bonds equivalent in their bonding character. Applied to density functional theory calculations on delta phase Pu, the 3Q spin structure results in elastic constants and phonons with the correct symmetry. Finally, the calculated phonon dispersion agrees better with experiment than previous calculations, and the calculated thermal expansion shows the unique behavior seen experimentally: it is negative.

36 MATERIALS SCIENCE↗
NASA NTRS1984

Prediction of fatigue crack-growth patterns and lives in three-dimensional cracked bodies

Fatigue crack growth patterns and lives for surface cracks, surface cracks at holes, and corner cracks at holes in three dimensional bodies were predicted using linear-elastic fracture mechanics concepts that were modified to account for crack-closure behavior. The predictions were made by using stress intensity factor equations for these crack configurations and the fatigue crack-growth (delta K against rate) relationship for the material of interest. The crack configurations were subjected to constant-amplitude fatigue loading under either remote tension or bending loads. The predicted crack growth patterns and crack growth lives for aluminum alloys agreed well with test data from the literature.

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

Prediction of fatigue crack-growth patterns and lives in three-dimensional cracked bodies

Fatigue crack growth patterns and lives for surface cracks, surface cracks at holes, and corner cracks at holes in three dimensional bodies were predicted using linear-elastic fracture mechanics concepts that were modified to account for crack-closure behavior. The predictions were made by using stress intensity factor equations for these crack configurations and the fatigue crack-growth (delta K against rate) relationship for the material of interest. The crack configurations were subjected to constant-amplitude fatigue loading under either remote tension or bending loads. The predicted crack growth patterns and crack growth lives for aluminum alloys agreed well with test data from the literature.

Newman, J. C., Jr.↗
DOE OSTI2022

ORNL_AISD_NiNb

This dataset describes the nickel-niobium solid solution binary alloy, where the two constituent elements nickel (Ni) and niobium (Nb) are randomly placed on an underlying crystal lattice. This dataset for nickel-niobium (Ni-Nb) alloys available includes the formation energy and bulk modulus for each crystal structure. Each atomic sample has a disordered phase which is obtained starting from an initial regular crystal structure of type body-centered cubic (BCC), face-centered cubic (FCC), or hexagonal compact packed (HCP). The geometry optimization ensures that all the alloy samples reached the equilibrium with negative formation energy. We perform geometry optimizations using the LAMMPS simulation package [1], a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. We utilized the embedded atom model (EAM) potential for Ni and Nb developed in a previous study [2]. The potential could describe behaviors of the liquid and solid phases of Ni-Nb alloy. The structural factors and angular distributions of three atoms are well-matched with X-ray and ab initio-based molecular dynamics data. We prepared the three different crystals with different initial lattice parameters (3.52 Ã… for FCC, 3.32 Ã… for BCC, and 3.5 Ã… for HCP). We performed energy minimization in two steps. Firstly, we minimized the structures with an isotropic unit cell to minimize the side effects from our arbitrary lattice parameters for all other compositions. Then, we applied geometry optimization with a triclinic (non-orthogonal) unit cell to fully minimize the stress components to calculate the elastic constants. In this procedure, we chose 10,000 as the maximum number of allowable steps aimed at obtaining fully relaxed atomic geometries. The dataset consists of three sets of crystal structures. The first set contains 46,086 irregular crystal structures, each of them with 54 atoms, obtained through optimization starting from a regular BCC crystal structure. The second set contains 24,543 irregular crystal structures, each of them with 32 atoms, obtained through optimization starting from a regular FCC crystal structure. The third set contains 39,303 irregular crystal structures, each of them with 48 atoms, obtained through optimization starting from a regular HCP crystal structure. The atomic configurations within each set span the possible compositional range. The three sets have been unified in a global dataset, which is extremely heterogeneous in terms of crystal structures, lattice volumes, and atomic configurations. Organization of files inside the dataset: the dataset contains three subdirectories called • BCC_opt • FCC_opt • HCP_opt based on the type of initial regular structure used to start the geometry optimization. Inside each of these folders, every atomic structure is identified by a string “A_B_Câ€, where A denotes the number of Nb in the system, B denotes index of structure with a given Nb number, and C denotes the total number of structures generated with a given Nb number. For each optimized crystal structure identified by the unique string of characters “A_B_Câ€, three files are provided: • A_B_C_opt.xyz: The optimized geometries in xyz format • A_B_C_opt.cfg: The optimized geometries in cfg format. It includes cell information and atomic energy, and forces calculated from LAMMPS. • A_B_C.elastic: Raw data of 21 elastic constants from LAMMPS output. • A_B_C.bulk: Calculated upper and lower bounds of bulk modulus and averaged one based on Voigt-Reuss-Hill approach from *.elastic. References: [1] A. P. Thompson, H. M. Aktulga, R. Berger, D. S. Bolintineanu, W. M. Brown, P. S. Crozier, P. J. in 't Veld, A. Kohlmeyer, S. G. Moore, T. D. Nguyen, R. Shan, M. J. Stevens, J. Tranchida, C. Trott, and S. J. Plimpton. LAMMPS - a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. Comp. Phys. Comm., 271:108171, 2022. [2] Y Zhang, R Ashcraft, MI Mendelev, CZ Wang, and KF Kelton. Experimental and molecular dynamics simulation study of structure of liquid and amorphous ni62nb38 alloy. The Journal of chemical physics, 145(20):204505, 2016.

36 MATERIALS SCIENCE↗
DOE OSTI2025

PyDislocDyn: A Python code for calculating dislocation drag and other crystal properties

PyDislocDyn is a suite of python programs designed to perform various calculations for dislocation dynamics in the continuum limit. In particular, one of its main purposes is to calculate dislocation drag from phonon wind. Additional features include the averaging of elastic constants for polycrystals, the calculation of the dislocation field including its limiting velocities, and the calculation of dislocation self-energy and line tension.

36 MATERIALS SCIENCE↗
DOE OSTI2022

LAVA

We introduce LAVA, a general-purpose python toolkit to provide user-friendly and high throughput calculations for material properties using both LAMMPS and VASP. It contains a set of classes as well as pre-processing and post-processing functions to prepare, execute and extract information from LAMMPS/VASP simulations. An overview of the program structure is provided. The current version contains modules to calculate elastic and mechanical properties such as lattice constant, cohesive energy, cold curve, elastic constants, bulk and shear modulus, volume conserving/non-conserving deformation path, vacancy/interstitial formation energy, surface energy, stacking fault energy, melting point, radial distribution function, and thermal expansion.

Dang, Khanh↗
NASA NTRS1987

A hybrid-stress finite element approach for stress and vibration analysis in linear anisotropic elasticity

A hybrid stress finite element method is developed for accurate stress and vibration analysis of problems in linear anisotropic elasticity. A modified form of the Hellinger-Reissner principle is formulated for dynamic analysis and an algorithm for the determination of the anisotropic elastic and compliance constants from experimental data is developed. These schemes were implemented in a finite element program for static and dynamic analysis of linear anisotropic two dimensional elasticity problems. Specific numerical examples are considered to verify the accuracy of the hybrid stress approach and compare it with that of the standard displacement method, especially for highly anisotropic materials. It is that the hybrid stress approach gives much better results than the displacement method. Preliminary work on extensions of this method to three dimensional elasticity is discussed, and the stress shape functions necessary for this extension are included.

Oden, J. Tinsley↗
DOE OSTI2022

Short- vs . long-range elastic distortion: structural dynamics of a [2 × 2] tetrairon( II ) spin crossover grid complex observed by time-resolved X-Ray crystallography

Spin crossover complexes (SCO) are among the most studied molecular switches due to their potential use in displays, sensors, actuators and memory components. A prerequisite to using these materials is the understanding of the structural changes following the spin transition at out-of-equilibrium conditions. So far, out-of-equilibrium studies in SCO solids have been focused on mononuclear complexes, though a growing number of oligonuclear SCO complexes showing cooperative effects are being reported. Here, we use time-resolved pink Laue crystallography to study the out-of-equilibrium dynamics of a [2 × 2] tetranuclear metallogrid of the form [Fe II 4 L Me 4 ](BF 4 ) 4 ·2MeCN ([L Me ] - = 4-methyl-3,5-bis{6-(2,2'-bipyridyl)}pyrazolate). The out-of-equilibrium spin state switching induced by a ps laser pulse demonstrates that the metallogrid exhibits a multi-step response similar to that reported for mononuclear complexes. Contrary to the mononuclear complexes, the metallogrid shows two types of elastic distortions at different time scales. The first is a short-range distortion that propagates over the entire Fe 4 grid complex during the ps time scale, and it is caused by the rearrangement of the coordination sphere of the photo-switching ion and the constant feedback between strongly linked metal ions. The second is a long-range distortion caused by the anisotropic expansion of the lattice during the ns time scale, observed in mononuclear materials. The structural analysis demonstrates that the long-range prevails over the short-range distortion, inducing the largest deformation of both the entire grid and the coordination sphere of each metal ion. The present study sheds light on the out-of -equilibrium dynamics of a non-cooperative oligonuclear complex.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1990

Inelastic deformation under nonisothermal loading

The objective of this paper is to evaluate, both experimentally and analytically, the appropriate forms of the hardening evolution equations in unified constitutive models for conditions involving nonisothermal loading. Critical experiments were performed for a cast nickel-base superalloy by using variable temperature tensile, creep, and cyclic tests in the 538-982 C temperature range. These experimental results were compared with both isothermal data and predictions of the Bodner-Partom-Bodner, 1987 elastic-viscoplastic theory to assess the effects of thermal history on constitutive behavior. The results indicate that the hardening evolution equations based on isothermal data are applicable for nonisothermal loading of these precipitation strengthened alloys. Additional thermal history effect terms in the hardening evolution equations were not required beyond those accounting for the variation of material constants with temperature. Using material constants determined solely from isothermal data, the inelastic deformation behavior of B1900 + Hf subject to thermomechanical loading were adequately predicted by the Bodner-Partom model.

Chan, K. S.↗
NASA NTRS2005

Nonlinear Visco-Elastic Response of Composites via Micro-Mechanical Models

Micro-mechanical models for a study of nonlinear visco-elastic response of composite laminae are developed and their performance compared. A single integral constitutive law proposed by Schapery and subsequently generalized to multi-axial states of stress is utilized in the study for the matrix material. This is used in conjunction with a computationally facile scheme in which hereditary strains are computed using a recursive relation suggested by Henriksen. Composite response is studied using two competing micro-models, viz. a simplified Square Cell Model (SSCM) and a Finite Element based self-consistent Cylindrical Model (FECM). The algorithm is developed assuming that the material response computations are carried out in a module attached to a general purpose finite element program used for composite structural analysis. It is shown that the SSCM as used in investigations of material nonlinearity can involve significant errors in the prediction of transverse Young's modulus and shear modulus. The errors in the elastic strains thus predicted are of the same order of magnitude as the creep strains accruing due to visco-elasticity. The FECM on the other hand does appear to perform better both in the prediction of elastic constants and the study of creep response.

Gates, Thomas S.↗
NASA NTRS2009

Prediction of Material Properties of Nanostructured Polymer Composites Using Atomistic Simulations

Atomistic models of epoxy polymers were built in order to assess the effect of structure at the nanometer scale on the resulting bulk properties such as elastic modulus and thermal conductivity. Atomistic models of both bulk polymer and carbon nanotube polymer composites were built. For the bulk models, the effect of moisture content and temperature on the resulting elastic constants was calculated. A relatively consistent decrease in modulus was seen with increasing temperature. The dependence of modulus on moisture content was less consistent. This behavior was seen for two different epoxy systems, one containing a difunctional epoxy molecule and the other a tetrafunctional epoxy molecule. Both epoxy structures were crosslinked with diamine curing agents. Multifunctional properties were calculated with the nanocomposite models. Molecular dynamics simulation was used to estimate the interfacial thermal (Kapitza) resistance between the carbon nanotube and the surrounding epoxy matrix. These estimated values were used in a multiscale model in order to predict the thermal conductivity of a nanocomposite as a function of the nanometer scaled molecular structure.

Hinkley, J.A.↗
DOE OSTI2021

A density functional theory and neutron diffraction study of the ambient condition properties of sub-stoichiometric yttrium hydride

Several mechanical and thermophysical properties are required as a function of non-stoichiometry for the successful implementation of YH 2-x for nuclear reactor moderator applications. Density functional theory calculations, in combination with neutron diffraction experiments, were used to study the structural and mechanical properties of YH 2-x . Point defect analysis indicated H occupation primarily at the tetrahedral site within an fcc Y sub-lattice, confirming a fluorite YH 2 structure. The small positive formation energy for H vacancies under Y-rich conditions predicted that hypo-stoichiometry is accommodated by Y+YH 2 at ambient conditions and by H vacancies in the YH 2-x single phase that is relevant to high temperatures. Neutron diffraction studies were used to confirm both the occupation of H on tetrahedral sites and the near-stoichiometric composition of the hydride phase in the two-phase Y+YH 2 region of the phase diagram that dominates at room temperature. Energy minimized special-quasirandom-structures of H vacancies were used to calculate lattice parameters, elastic constants, and several other properties as a function of composition for the single phase YH 2-x . The lattice parameter of YH 2-x decreased by only 0.004 Å with increasing H/Y for the range 1.31 ≤ H / Y ≤ 2.0 indicating a negligible effect on lattice parameters due to vacancy formation. In the two-phase region, however, calculations predicted the density of z Y + ( 1 - z ) Y H 2 - x to increase with decreasing H/Y at lower temperatures due to the increased fraction of high-density Y metal. Additionally, for the high temperature single phase, decreasing H/Y reduced the density as a consequence of the lattice expansion associated with vacancy formation. All elastic constants and moduli increased with increasing hydrogen content in single-phase YH 2-x .

36 MATERIALS SCIENCE↗
NASA NTRS1997

Structure-Property Relationships of Bismaleimides

The purpose of this research was to control and systematically vary the network topology of bismaleimides through cure temperature and chemistry (addition of various coreactants) and subsequently attempt to determine structure-mechanical property relationships. Characterization of the bismaleimide structures by dielectric, rheological, and thermal analyses, and density measurements was subsequently correlated with mechanical properties such as modulus, yield strength, fracture energy, and stress relaxation. The model material used in this investigation was 4,4'-BismaleiMidodIphenyl methane (BMI). BMI was coreacted with either 4,4'-Methylene Dianiline (MDA), o,o'-diallyl bisphenol A (DABA) from Ciba Geigy, or Diamino Diphenyl Sulfone (DDS). Three cure paths were employed: a low- temperature cure of 140 C where chain extension should predominate, a high-temperature cure of 220 C where both chain extension and crosslinking should occur simultaneously, and a low-temperature (140 C) cure followed immediately by a high-temperature (220 C) cure where the chain extension reaction or amine addition precedes BMI homopolymerization or crosslinking. Samples of cured and postcured PMR-15 were also tested to determine the effects of postcuring on the mechanical properties. The low-temperature cure condition of BMI/MDA exhibited the highest modulus values for a given mole fraction of BMI with the modulus decreasing with decreasing concentration of BMI. The higher elastic modulus is the result of steric hindrance by unreacted BMI molecules in the glassy state. The moduli values for the high- and low/high-temperature cure conditions of BMI/MDA decreased as the amount of diamine increased. All the moduli values mimic the yield strength and density trends. For the high-temperature cure condition, the room- temperature modulus remained constant with decreasing mole fraction of BMT for the BMI/DABA and BMI/DDS systems. Postcuring PMR-15 increases the modulus over that of the cured material even though density values of cured and postcured PMR were essentially the same. Preliminary results of a continuous and intermittent stress relaxation experiment for BMI:MDA in a 2:1 molar ratio indicate that crosslinking is occurring when the sample is in the undeformed state. Computer simulation of properties such as density, glass transition temperature, and modulus for the low- temperature cure conditions of BMI/MDA and BMI/DABA were completed. The computer modeling was used to help further understand and confirm the structure characterization results. The simulations correctly predicted the trends of these properties versus mole fraction BMI and were extended to other BMI/diamine systems.

Tenteris-Noebe, Anita D.↗
DOE OSTI2022

The fate of liquid crystal topological defects on chemically patterned surfaces during phase transitions

Controlling topological defects in liquid crystals (LCs) is an essential element in the development of areas such as directed self-assembly and micropatterning materials. However, during the phase transition on confined patterned surfaces, how the morphologies in one liquid crystalline phase change from deformations or defects into another phase is much less known. Here, we examine the fate of defects in a LC confined on a patterned surface during smectic-A–nematic and nematic–isotropic phase transitions, using experiments and simulation analyses. Upon heating from smectic-A to nematic, a Toric focal conical domain (TFCD) melts into a +1 converging boojum defect, which then transitioned into a concentric configuration as temperature increases, attributed to a steeper decrease of the bend and twist modulus compared to splay modulus. During cooling, TFCDs are developed from two distinct pathways depending on the cooling rates. Our continuum simulation recapitulates these transformations and provides elastic constant-based explanations for the two pathways. Although the phase transition pathways of defects are independent of the geometry of the confined patterns, the arrangement of FCDs is highly dependent on the size and shape of the patterns. Finally, taken together, this simple approach offers promising opportunities for tuning the micro- or nano-patterning of topological defects in liquid crystals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1988

Stress concentration around circular hole in a composite material specimen representative of the X-29A forward-swept wing aircraft

The theory of anisotropic elasticity was used to evaluate the anisotropic stress concentration factors of a composite laminated plate containing a small circular hole. This advanced composite material was used to manufacture the X-29A forward swept wing. Observe that the usual isotropic material stress concentration factor is three. However, for composite material, it was found that the anisotropic stress concentration factor is no longer constant, and that the locations of maximum tangential stress points could shift by changing the fiber orientation with respect to the loading axis. The analysis showed that through the lamination process, the stress concentration factor could be drastically reduced, and therefore the structural performance could be improved. Both the mixture rule approach and the constant strain approach were used to calculate the stress concentration factor. The results predicted by the mixture rule approach were about 20 percent deviate from the experimental data. However, the results predicted by the constant strain approach matched the testing data very well. This showed the importance of the inplane shear effect on the evaluation of stress concentration factor for the X-29A composite plate.

Yeh, Hsien-Yang↗