Search NASA⌕ Search

SEARCH · Search NASA

Results for “stress-strain model”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Unified Stress-Strain Model for Plasticity to the Structural Instability

A unified model for the work hardening Θσ Θ(σ) and stress-strain σε σ(ε) behavior is presented that accounts for deformation under tensile loading, from the onset of yielding at the proportional limit up to the ultimate strength as defined at the structural instability. The origin of this approach is based on a negative exponential formulation for an asymptotic-curvilinear work-hardening model that accounts for the rapid strengthening of metals as well as the continuation of steady-state deformation to the instability.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Work Hardening Model of Structure in Hall-Petch Strengthening

The microstructural length scale of metals changes by orders of magnitude under extreme processing conditions producing a concurrent wide range of mechanical strength and plasticity behaviors. A unified stress-strain σε σ(ε) model is formulated that’s based on superposing the components of asymptotic-curvilinear work hardening Θσ Θ(σ) to qualify and quantify these mechanical behaviors. This approach accounts for the rapid strengthening of metals beyond the initial yield point, through extended steady-state deformation, to the structural instability. The relationship between the softening coefficients cbi c bi of the work hardening formulation Θσ Θ(σ) and strength are found to reveal the microstructural scale in the material. Specifically, the rapid decrease in the slope of the Θσ Θ(σ) curve provides a measure for microstructural size consistent with a functional Hall-Petch relationship of strength. A successful application is shown for the tensile behavior of pure aluminum microstructures that result from extreme plastic deformation by equal-channel angle pressing.

Hall-Petch strengthening↗
DOE OSTI2024

Brittle failure analysis and modeling of high-burnup PWR fuel cladding alloys

The aim of this research is the development of methods for predicting mechanical behavior and identification of limiting conditions to prevent brittle failure of high-burnup (HBU) pressure water reactor (PWR) fuel cladding alloys. A finite element (FE) model of the ring compression test (RCT) was created to analyze the failure behavior of zirconium-based alloys with radial hydrides during the RCT. An elastic-plastic material model describes the zirconium alloy. The stress-strain curve needed for the elastic-plastic material model was derived by inverse finite element analyses. Cohesive zone modeling is used to reproduce sudden load drops during RCT loading. Based on the failure mechanism in non-irradiated ZIRLO (R) claddings, a micro-mechanical model was developed that distinguishes between brittle failure along hydrides and ductile failure of the zirconium matrix. Two different cohesive laws representing these types of failure are present in the same cohesive interface. The key differences between these constitutive laws are the cohesive strength, the stress at which damage initiates, and the cohesive energy, which is the damage energy dissipated by the cohesive zone. Statistically generated matrix-hydride distributions were mapped onto the cohesive elements and simulations with focus on the first load drop were performed. Computational results are in good agreement with the RCT results conducted on high-burnup M5 (R) samples. It could be shown that crack initiation and propagation strongly depend on the specific configuration of hydrides and matrix material in the fracture area.

Simbruner, Kai↗
DOE OSTI2025

Physically constrained 3D diffusion for inverse design of fiber-reinforced polymer composite materials

Designing fiber-reinforced polymer composites (FRPCs) with a tailored nonlinear stress-strain response is crucial for applications such as energy absorption in crash structures, flexible robotics, and impact-resistant protective gear. However, the inherent complexities of composite materials and the multitude of parameters involved, render traditional design and optimization methods inadequate for achieving effective inverse design of composites. In this paper, we present an AI-based inverse design framework that effectively and efficiently generates FRPCs with targeted nonlinear stress-strain responses. We introduce a physically constrained diffusion model (PC3D_Diffusion) capable of managing the complexities of composite materials and producing detailed, high-quality designs. We propose a loss-guided, learning-free approach to generate physically feasible microstructure designs by explicitly enforcing physical constraints during the generation process. For training purposes, 1.35 million FRPC samples were created, and their corresponding stress-strain curves were computed using established physics-based computational models. The results show that PC3D_Diffusion consistently generates high-quality designs with tailored mechanical behaviors, while guaranteeing compliance with the physical constraints. PC3D_Diffusion advances FRPC inverse design and may facilitate the inverse design of other 3D materials, offering potential applications in industries reliant on materials with custom mechanical properties.

Xu, Pei [Clemson Univ., SC (United States)]↗
DOE OSTI2025

Inverse design of cellular structures with the targeted nonlinear mechanical response

Advanced additive manufacturing capabilities have enabled a transformational ability to create sophisticated cellular structures using diverse materials. By altering the topology of the unit cell, the mechanical behavior, such as the stress-strain response during compression, can be modulated. Nevertheless, identifying a printable topology within an enormous design space that would precisely deliver the targeted nonlinear material response is challenging. We propose a data-driven generative framework based on a conditional variational autoencoder (cVAE) architecture that can inverse design the cellular structure based on the intended nonlinear stress-strain response. Trained on a dataset of structure-property pairs, the cVAE learns a compact and expressive latent space that enables efficient mapping from targets to feasible geometries. Two inference modes are explored: (1) decoder-only generation, which enables the exploration of diverse designs conditioned solely on the desired mechanical response, and (2) encoder-decoder generation, which further allows for the incorporation of desired topologies, ensuring the generated structure conforms to both mechanical properties and to desired-topology constraints. The results demonstrate that the model can generate structurally plausible and mechanically accurate designs, with the predicted stress-strain curves closely matching the targets. Even under joint conditioning, the model effectively balances geometric fidelity and functional performance.

36 MATERIALS SCIENCE↗
DOE OSTI2026

A finite viscoelastic constitutive model for low to high strain rate response of elastomers with application of strain rate-induced glass transition

Amorphous elastomers exhibit significant rate-stiffening and unique viscous flow characteristics across a wide range of strain rates, often undergoing glass transition above a strain rate threshold. We have developed a thermodynamically-consistent and micromechanically-inspired constitutive model for soft elastomeric materials to capture the rate-dependent stress-strain behavior and hysteresis when subjected to low to high strain rates. Here, our proposed constitutive model encapsulates the viscous flow of materials through molecular motion at low strain rates and local rearrangement and alignment of the molecules trying to overcome the intermolecular resistance at high strain rates, essentially covering the glass transition. We applied our constitutive model to uniaxial compression experiments performed at low and high strain rates for polyborosiloxane (PBS) to identify the material parameters, and subsequently, performed numerical simulations of single and multi-cycle compression, stress relaxation, and small amplitude oscillatory tension-compression. Our analyses indicate that the model predicts higher total energy dissipation with increasing strain rate; however, dissipation associated with molecular relaxation decreases (forming a cusp) because, beyond a crossover strain rate, intermolecular rearrangement and alignment become dominant, which is consistent with the onset of the glass transition. For cyclic loading-unloading, we observed that dissipation over a cycle remains constant at low strain rates but decreases non-monotonically at high strain rates before becoming constant, with the peak stress over the cycle becoming higher, which can be interpreted as more loading being carried elastically by the polymer network as the intermolecular rearrangement process occurs. Additionally, our model was able to predict the qualitative nature of the storage modulus and loss modulus in the limit of small strain over a wide range of frequency sweeps.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Machine-Learning-Based Multiscale Methods for 3D Modelling of Granular Materials by Incorporating History-Dependent State Variables

Over the past decades, the prevalence of machine learning (ML) methods has made the development of ML-based constitutive models for granular materials undoubtedly a popular subject. Numerous studies have been made to feature the loading path or history-dependent stress-strain response of granular media using neural networks. In this work, a novel finite element method (FEM)–ML multiscale approach was developed by incorporating internal variables to improve the simulation accuracy of 3D history-dependent granular materials for the first time. To this end, a surrogate constitutive model based on the single-step-based multi-layer perceptron (MLP) neural network was used to replace representative volume element (RVE) simulations conducted by the discrete element method (DEM) in the multiscale FEM–DEM approach. Although the prediction principle of the MLP aligns with the FEM algorithm, artificially added internal variables are required to differentiate the loading history. To address this issue, history variables associated with the Frobenius norm are proposed to be fed into the MLP coupled with the strain tensor to extract the history-dependent behaviour of granular assemblies. The developed FEM–ML approach was demonstrated in 3D conventional triaxial compression (CTC) simulations. Compared to the multiscale FEM–DEM approach, the proposed FEM–ML method exhibits a significantly improved computational efficiency.

granular materials↗
DOE OSTI2025

In-Silico Design of Next Generation Cellulose-Derived Packaging Materials (CRADA Final Report)

Developing sustainable solutions for single-use packaging is an important objective to combat the environmental crisis of plastics pollution. Most embodiments of cellulose-based packaging materials, including CellophaneTM, are completely biodegradable in both terrestrial and marine environments. However, petroleum-derived alternatives offer some performance advantages for metrics such as moisture barriers and mechanical properties. This project leverages molecular dynamics simulation to investigate how molecular modifications to cellulose-based polymer assemblies impact their material properties. An important performance criterion for the modified materials was to retain biodegradability; thus, modifications by naturally occurring, biodegradable additives were the focus of this study. Specifically, we developed models with xylan and lignin of varying monomeric compositions into the cellulose matrix. The mechanical properties were investigated by performing stress-strain simulations, and the water barrier and hydrophobicity were investigated by simulating the water contact angle. Our findings indicate that the incorporation of xylan into the cellulose matrix tends to increase the mechanical properties with an optimal loading of ~27 wt%. We also predict that orienting the nanoscale directionality of the xylan chains such that they are perpendicular to the cellulose fibrils will dramatically increase mechanical strength. In contrast, the incorporation of lignin tends to weaken the composite at all loadings investigated. Simulations of water contact angle predicted that coating polymers on the surface of the cellulose assembly creates a more hydrophobic surface than incorporating them throughout the matrix. Of the coatings investigated, lignin resulted in the most hydrophobic surface, followed by pectin and keratin, which both imparted modest increases in hydrophobicity. Future experimental work done by Futamura will focus on designing material prototypes to capitalize on the predictions of performance enhancement obtained from molecular modeling. While substantial progress was made by the simulations performed in this project, there still exists a vast parameter space that we were unable to investigate, including branching, functional group decoration, and degree of polymerization of polymer additives. However, the methods developed in this initial investigation will facilitate more rapid evaluation of the impact of molecular characteristics on the performance of biopolymer composite materials and thereby accelerate future materials discovery efforts in this area.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Step-loaded creep testing of Zircaloy-4 cladding at higher temperatures in the α-phase

A refined understanding of zirconium-based cladding thermomechanical performance during rapid transients is essential for enhancing the safety and operation of light-water reactors. Traditional models for zirconium alloys under accident conditions generally assume that creep dominates fuel cladding performance. Here, these historic models have largely remained unchanged and serve as the basis for safety criteria development. As the U.S. nuclear industry pursues higher burnup levels, the increased release of fission gases during transients raises the risk of cladding failure in the low-temperature hcp α-phase, making the fidelity of these models of greater importance. Creep testing was conducted from 550–700°C with 25–120 MPa applied hoop stresses to investigate Zircaloy-4 deformation at accident-relevant temperatures in the α-phase. Step-loading was employed to capture creep behavior across a wide stress range from a single sample. The stress-strain rate data at higher temperatures (650 and 700°C) were well-described by isotropic versions of the Erbacher and Kaddour models, while the lower temperature data (550 and 600°C) were underpredicted by both anisotropic and isotropic model variants. Greater strain rates during the initial loading step at 650 and 700°C were attributed to recrystallization and grain growth of sub-micron crystallites. Yet, texture analysis revealed the basal split texture remained after testing. These observations produced results suggesting Zircaloy-4 claddings experience higher creep rates across the α-phase than previously thought, possibly related to dynamic anisotropy due to temperature dependent activation of deformation mechanisms, effects of biaxial loading, and variation in material condition between the current testing used in previous model development.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Mechanical Characterization of Electrolyzer Membranes and Components Under Compression

Proton-exchange membrane (PEM) water electrolysis is a promising technology for producing clean hydrogen by electrochemically splitting water when paired with renewable energy sources. A major roadblock to improving electrolyzer durability is the mechanical degradation of the cell components, which requires an understanding of their mechanical response under device-relevant conditions. However, there is a lack of studies on the mechanical characterization of the PEM and other components, as well as and their interactions. This study aims to address this gap by using a custom-designed testing apparatus to investigate the mechanics of electrolyzer components in uniaxial compression at 25 and 80 °C. Findings show stress-strain response of components have a varying degree of nonlinearity owing to their distinct deformation mechanisms and morphologies, from porous structures to polymers. These results are used to develop an expression for compressive stress-strain response of Nafion membranes and then analyze the deformation of components under applied pressure by using a 1-D spring network model of cell assembly. This work provides a new understanding of mechanical responses of the electrolyzer membrane and cell components, which can help assess material design and cell assembly strategies for improved electrolyzer durability.

08 HYDROGEN↗
DOE OSTI2025

Grain2Mesh: Mesh Generation for Grain-Scale Nonlinear Elasticity Modeling

The nonlinear hysteretic behavior of rocks under cyclic loading is a crucial area of study in geomechanics. The macroscopic response of a variety of materials has been found to be contingent upon the behavior of the micro-scale structure. This project aims to develop a functional and maintainable software package for generating a multi-phase numerical mesh and accompanying simulation files for finite element modeling used in computational mechanics solvers. Meshes generated from images often lack key preprocessing that reduces noise and prevents mesh element distortion that can increase computational cost. By incorporating user feedback throughout, grain2mesh ensures a high-fidelity mesh that can be used to model grain-scale interactions such as shearing, crack propagation, and interfacial material contrast. Scientific applications of this software include material fracturing, stress-strain analysis for natural and engineered materials, and nonlinear meso-scale analysis.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Anomalous softening of 3D printed elastomeric foam irradiated under compressive strain

Elastomeric foam is an essential component in many industrial and technological settings, primarily as thermal insulators and as positional/mechanical support cushions. In particular, silicone foam is utilized in harsh environments due to exceptional thermal and chemical stability. Under service conditions within certain applications such material gets exposed to a high dosage of gamma radiation, which can permanently alter the material’s structural and mechanical response properties. Most studies on gamma-exposure under inert or oxidative atmosphere indicate hardening of silicone foam, which is attributed to radiation-induced enhancement in chemical cross-linking. Here we report two contrasting effects depending on whether (non-oxidative) radiation exposure is carried out with the foam under zero or finite compressive strain. While in the former case we observe radiation-hardening consistent with previous studies, in the latter case (50% porous foam under 30% uniaxial compression) we see a monotonic decrease in Young’s modulus with increasing dosage, although solvent swelling experiments on the constituent rubber indicate a net increase in cross-link density independent of the state of strain. We quantitatively model all dose-dependent data using the Ogden Hyperfoam strain-energy function within the framework of Tobolsky two-network scheme and attribute the above anomaly to a combined effect of radiation-induced thickness change (compression set) and inherent nonlinearity in the foam’s stress-strain response.

Coarse-grained models↗
DOE OSTI2025

Micro-tensile characteristics of As-fabricated and irradiated AGR-2 TRISO fuel particle buffer, IPyC, and buffer-IPyC interlayer regions

A recently developed micro-tensile sample preparation technique was implemented to evaluate the tensile strengths of the buffer, IPyC, and buffer-IPyC interlayer regions of the unirradiated and irradiated AGR-2 TRISO fuel particles. Understanding the mechanical properties of the buffer-IPyC interlayer is essential for developing thermomechanical models of buffer-IPyC separation, yet there is a lack of experimental data on its micro-tensile properties. TEM analysis was conducted on these regions to determine the microstructural changes relevant to the samples' tensile properties. In the unirradiated TRISO particle samples, the buffer layer demonstrated the weakest tensile strength, while the IPyC layer exhibited the highest. Conversely, in the irradiated TRISO particle samples, the buffer-IPyC interlayer region showed the lowest tensile strength, with the IPyC layer being the strongest. Fractures in the samples from the buffer-IPyC region predominantly occurred either in the buffer layer or at the buffer-IPyC interface. However, some buffer-IPyC interlayer samples displayed stress-strain and fracture behaviors more akin to the IPyC layer than the buffer layer. Analysis of diffraction patterns suggests that irradiation may have increased anisotropy in the three regions tested. Despite this suggested increase in anisotropy, there was no evidence that it affected the measured strengths. The irradiated TRISO particles demonstrated a considerable increase in void space and a decrease in ultimate tensile strength within the buffer-IPyC interlayer region due to the densification and contraction of the buffer layer. Minor variations in diffraction ring patterns were also observed. These changes, coupled with a significant reduction in the Weibull modulus/shape parameter, imply that irradiation-induced densification leads to tearing between the buffer and IPyC layers at locations of elevated porosity in the buffer-IPyC interlayer region.

Tristructural isotropic (TRISO)↗
DOE OSTI2025

Fused filament fabrication of thermoplastic polyurethane composites with microencapsulated phase-change material

Here, the present study examines the thermal energy storage (TES) effectiveness and printability of microencapsulated phase-change material (MEPCM) combined with thermoplastic polyurethane (TPU) for fused filament fabrication (FFF). Two formulations were assessed: 24D MEPCM, which changes phase at 24 ° C, compounded with TPU pellets and 43D MEPCM, which changes phase at 43 ° C, integrated with TPU powder. These combinations are designed to evaluate the effectiveness of the form of the TPU (pellets versus powder) in the FFF process. The investigation includes a comprehensive analysis of thermal characteristics, encompassing phase-change temperature, latent heat of fusion, thermal conductivity, and thermal decomposition, which are assessed through differential scanning calorimetry (DSC) and thermogravimetric analysis (TGA). Additionally, mechanical properties, including stress-strain behavior, are examined to evaluate material suitability for TES applications, while microstructural visualization is used to provide deeper insights into material performance, structural integrity, and the quality of printed components. The 24D MEPCM and TPU pellets formulation experienced a significant loss of approximately 39.6% of PCM during filament extrusion and printing, resulting in a reduced effective latent heat. Therefore, further characterization of the pellet formulation was discontinued due to excessive leakage. In contrast, the 43D MEPCM and TPU powder formulation demonstrated minimal PCM loss, with the 60 wt.% composition achieving an effective latent heat of 132 J/g. This value represents the highest effective latent heat currently documented in the literature for PCM-polymer-composite materials produced using an FFF-based additive manufacturing process.

25 ENERGY STORAGE↗
DOE OSTI2025

Uncertainty Quantification for Smooth Functional Data with Application to Material Properties

This document outlines a method for processing functional output (i.e., curves) for the ultimate purpose of sampling curves under specified input conditions for use in modeling and simulation uncertainty quantification (UQ) studies. A set of benchmark curves sufficiently representative of the relevant scenario(s) being simulated are provided to the process and formatted as described in Section 1. Principal Component Analysis (PCA) is utilized to discover the components of uncertainty in the benchmark curves and is outlined in Section 2. Section 3 describes the application of uncertainty quantification to the PCA results for the purpose of sampling curves to be used in UQ analysis. Section 4 applies these techniques to an example benchmark dataset. Concluding remarks are provided in the final section.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Measurement of the Static Nonlinear Third-Order Elastic Moduli of Rocks: Problems and Applicability

The third-order elastic (TOE) model has been used to describe the widely observed nonlinear mechanical behaviors of earth materials. In addition to linear elastic constants ( λ , μ ), three nonlinear elastic moduli ( A , B , C ) are required for isotropic rocks. Contrary to previous research on dynamic TOE moduli, this study followed the protocol to measure strain and stress under uniaxial and hydrostatic compressive tests statically, which were later used to invert for the full set of TOE moduli for four standard rock types with differing pore structures of a porous oolitic limestone, quartz-rich sandstones, and a dense crystalline basalt. The applicability of the TOE model to characterize nonlinearity depends on the fulfillment of path-independence and small-strain assumptions. Using the measured static TOE moduli, the finite element model demonstrates that the stress in the vicinity of the wellbore is more amplified than the stress in the linear elastic case, which leads to a wider zone of rock failure around the wellbore. Due to the long-standing discrepancy between static and dynamic moduli, the rarely reported full set of static TOE moduli is necessary and will benefit future research in understanding the effect of rock nonlinearity on geophysical and geomechanical applications, such as long-term safe storage of CO 2 and generating process of geohazards.

58 GEOSCIENCES↗
DOE OSTI2026

Controlling Deformation in Al/Ti: How Interface Roughness and Orientation Drive Bimetal Mechanics

The microstructural characteristics and morphology of interfaces in metals can be crucial in governing the initiation of plasticity and early deformation mechanisms under extreme conditions. During high-strain-rate deformation, these interfaces significantly affect dislocation nucleation, twinning, and other mechanisms that directly impact material strength and failure. Despite their importance, a substantial knowledge gap remains between the observed macroscopic material behavior and the underlying role of bimetal interfaces in plasticity initiation. Here, to address this gap, large-scale molecular dynamics simulations are performed on Al/Ti bimetal structures to examine the effect of interface characteristics on the onset of plasticity under uniaxial compression. Specifically, this study investigates how interface roughness (flat vs waveform interfaces) modifies the initiation of plastic events in the microstructure. Atomistic simulations indicate that interface roughness (a microscopic behavior) reduces the stress required for dislocation nucleation, thereby reducing the peak stress relative to a flat interface. For the square interface, the step height strongly influences plasticity initiation by setting the separation of locally flat regions. Varying the interface rotation relative to the loading direction (macroscopic behavior) reveals that the peak stress for both flat and waveform interfaces initially decreases and then increases with rotation, accompanied by a shift in dominant mechanisms—from Al twinning (flat) to interface sliding and Ti-dominated phase transformation and twinning at lower angles. The change in the stress-strain slope during initial compression reflects how rotation alters the resolved shear stress and activates different slip systems. Overall, the study provides valuable insights into the role of bimetal interfaces in controlling plasticity initiation and early deformation pathways for Al/Ti under extreme loading conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2026

ACT University Collaboration Proposal: Role of manufacturing defects on material failure under dynamic loading for developing enhanced failure models and theories (Final Report)

Research on the tensile behavior of additively manufactured 316L stainless steel coupon specimens at increasing strain rates was conducted over the past year at Penn State. Dynamic loading rates into the 1000 1/s loading rates were performed on specimens with nearly 1.32” gage lengths. Stress-strain plots show ductile and plastic behavior well beyond 12% strain with necked specimens, having smaller effective gage lengths showing up to 56% ultimate strain. Additional tests performed on compact tension specimens helped with simulation work to understand the deformation behavior. Using finite elements, it was possible to determine the feasibility of a comprehensive experimental validation study towards an improved model for AM material failure – specifically, the Bai-Wierzbicki approach which accounts for lode angle and triaxiality. A non-significant number of tests are projected for eight specific failure nodes, with additional replicates to provide strain-rate capabilities to the existing formulation. One such correction factor is explored for quasi-static, notched specimens. Finally, stress intensity factor was explored using a set of compact-tension specimens which would also provide useful validation data for any simulation work.

36 MATERIALS SCIENCE↗