Search NASA⌕ Search

SEARCH · Search NASA

Results for “Constitutive modeling”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Calibration of CDM-Based Creep Constitutive Model Using Accelerated Creep Test (ACT) Data

In conventional creep testing (CCT) a specimen is subject to constant load and temperature for a long-duration until creep rupture occurs. Conventional testing can be costly when considering the number of experiments needed to characterize the creep response of a material over a range of stress and temperature. To predict long-term creep-rupture properties, the time-temperature-stress superposition principle (TTSSP) approach has been employed where stress and/or temperature is applied at an elevated level; the result of which are extrapolated down to low stress and/or temperature conditions. These methods have been successful in predicting minimum-creep-strain-rate (MCSR) and stress-rupture (SR) but suffer from an inability to predict the creep deformation curve or account for changes in deformation mechanisms or aging that occurs at long-duration. An accelerated technique, termed the Stepped isostress method (SSM) allows the accelerated testing of materials to determine their creep deformation response. Unlike TTSSP tests, the SSM test employs a single specimen where the stress is periodically step increased until rupture. The SSM creep deformation curve is processed (time and strain shifted) to produce an accelerated creep deformation curve that represent the creep deformation curve at the initial stress level in SSM. A processing procedure for metals has yet to be developed. The research objective of this study is to develop a processing procedure for SSM test data using a creep-damage constitutive model. Triplicate SSM tests were conducted on Ni-based superalloy Inconel 718 at 650°C with stress being periodically increased until rupture. Triplicate CCT tests were conducted at the initial stress level of the SSM tests. The Sine-hyperbolic (Sinh) creep-damage model was employed in this study. The Sinh creep-damage constitutive model is based on coupled creep strain rate and damage evolution equations; where both rates are dependent on the current state of damage. Calibration is two-step: analytical and numerical optimization. Each stepped creep deformation curve is tackled quasi-analytically to determine MCSR and SR related material constants and accumulated damage. The damage accumulated at the end of each step was then passed onto subsequent steps to calibrate the MCSR, rupture prediction, and damage evolution. Numerical optimization was applied to optimize model constants involved in the creep strain constitutive equations in order to generate best-fitted Sinh creep deformation curves. The Sinh model predictions were compared to the SSM and CCT data. The Sinh model satisfactorily predicts the SSM data and thus the calibrated material constants provides a good estimate of rupture found in the CCT data. Calibration using SSM data reduces the number of tests needed to calibrate a model; significantly reducing costs. A single SSM test replaces numerous creep tests at different stresses.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A constitutive model for glass–ceramic materials

Glass–ceramics have received recent attention for use in glass–ceramic to metal hermetic seals. Due to their heterogeneous microstructure, these materials exhibit a number of advantageous responses over conventional glass based seals. Key amongst them is the possibility of a controllable thermal strain response and apparent coefficient of thermal expansion which may be used to minimize thermally induced residual stresses for aforementioned seals. Furthermore, these behaviors result from an inorganic glass matrix and variety of crystalline ceramic phases including silica polymorph(s) that may undergo reversible solid-to-solid transformations with associated inelastic strain. Correspondingly, these materials exhibit complex thermomechanical responses associated with multiple inelastic mechanisms (viscoelasticity and phase transformation). While modeling these behaviors is essential for developing and analyzing the corresponding applications, no such model exists. Therefore, in this work a three-dimensional continuum constitutive model for glass–ceramic materials combining these various inelastic mechanisms is developed via an internal state variable approach. A corresponding fully implicit three dimensional numerical formulation is also proposed and implemented. The model is used to simulate existing experiments and validate the proposed formalism. As an example, the simple seal problem of a glass–ceramic seal inside a concentric metal shell is explored. Finally, the impact of the cooling rates, viscoelastic shift factors, and inelastic strain on final residual stress state are all investigated and the differing contributions highlighted.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A constitutive model for sheared dense fiber suspensions

Here we propose a constitutive model to predict the viscosity of fiber suspensions, which undergoes shear thinning, at various volume fractions, aspect ratios, and shear stresses/rates. We calibrate the model using the data from direct numerical simulation and prove the accuracy by predicting experimental measurements from the literature. We use a friction coefficient decreasing with the normal load between the fibers to quantitatively reproduce the experimentally observed shear thinning in fiber suspensions. In this model, the effective normal contact force, which is directly proportional to the bulk shear stress, determines the effective friction coefficient. A rise in the shear stress reduces the effective friction coefficient in the suspension. As a result, the jamming volume fraction increases with the shear stress, resulting in a shear thinning in the suspension viscosity. Moreover, we extend the model to quantify the effects of fiber volume fraction and aspect ratio in the suspension. We calibrate this model using the data from numerical simulations for the rate-controlled shear flow. Once calibrated, we show that the model can be used to predict the relative viscosity for different volume fractions, shear stresses, and aspect ratios. The model predictions are in excellent agreement with the available experimental measurements from the literature. The findings of this study can potentially be used to tune the fiber size and volume fraction for designing the suspension rheology in various applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Theory and Implementation of the Spectacular Nonlinear Viscoelastic Constitutive Model

This report is a comprehensive guide to the nonlinear viscoelastic Spectacular model, which is an isotropic, thermo-rheologically simple constitutive model for glass-forming materials, such as amorphous polymers. Spectacular is intermediate in complexity to the previous PEC and SPEC models (Potential Energy Clock and Simplified Potential Energy Clock models, respectively). The model form consists of two parts: a Helmholtz free energy functional and a nonlinear material clock that controls the rate of viscoelastic relaxation. The Helmholtz free energy is derived from a series expansion about a reference state. Expressions for the stress and entropy functionals are derived from the Helmholtz free energy following the Rational Mechanics approach. The material clock depends on a simplified expression for the potential energy, which itself is a functional of the temperature and strain histories. This report describes the thermo-mechanical theory of Spectacular, the numerical methods for time-integrating the model, model verification for its implementation in LAMÉ, a user guide for its implementation in LAMÉ, and ideas for future work. A number of appendices provide supplementary mathematical details and a description of the procedure used to derive the simplified potential energy from the full expression for the potential energy. The goal of this report is create a convenient point-of-entry for engineers who wish to learn more about Spectacular, but also to serve as a reference manual for advanced users of the model.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Initial High Temperature Inelastic Constitutive Model for Alloy 617

This report describes progress on a Alloy 617 inelastic constitutive model slated for inclusion in the ASME Boiler & Pressure Vessel Code in Code Case N-898 “Use of Alloy 617 (UNS N06617) for Class A Elevated Temperature Construction, Section III, Division 5” as part of an appendix on inelastic modeling. This appendix will provide a reference material model for A617 for use in conjunction with the ASME design by inelastic analysis criteria. This report focuses on the implementation and evaluation of an existing Alloy 617 model developed to represent the data collected at Idaho National Laboratory (INL) and Oak Ridge National Laboratory (ORNL) as part of the Alloy 617 Code Case experimental campaign. While the model accurately captures the cyclic stress/strain behavior of Alloy 617, it fails to capture the long-term creep and lower temperature tensile data. In addition, the model form is complicated and has over 50 configurable parameters. As such, we recommend the development of a simplified model, incorporating key features of the ORNL model but reducing the model complexity.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Complementing a continuum thermodynamic approach to constitutive modeling with symbolic regression

Here, an interpretable machine learning method, physics-informed genetic programming-based symbolic regression (P-GPSR), is integrated into a continuum thermodynamic approach to developing constitutive models. The proposed strategy for combining a thermodynamic analysis with P-GPSR is demonstrated by generating a yield function for an idealized material with voids, i.e., the Gurson yield function. First, a thermodynamic-based analysis is used to derive model requirements that are exploited in a custom P-GPSR implementation as fitness criteria or are strongly enforced in the solution. The P-GPSR implementation improved accuracy, generalizability, and training time compared to the same GPSR code without physics-informed fitness criteria. The yield function generated through the P-GPSR framework is in the form of a composite function that describes a class of materials and is characteristically more interpretable than GPSR-derived equations. The physical significance of the input functions learned by P-GPSR within the composite function is acquired from the thermodynamic analysis. Fundamental explanations of why the implemented P-GPSR capabilities improve results over a conventional GPSR algorithm are provided.

42 ENGINEERING↗
DOE OSTI2021

Formulation and calibration of two-dimensional constitutive models for composite structures based on panel tests

Steel Plate Concrete (SC) composite members have been widely adopted because of its cost-efficiency and enhanced structural behavior. While researchers have attempted to study its in-plane shear behavior in the past twenty years, very limited number of large-scale pure shear tests were performed due to the challenge of experimental set-up and the availability of facilities. In this paper, a series of uniaxial loading tests and two full-scale pure shear panel tests of SC members were reported, on which the “mechanics-based Membrane Model of SC elements (MM-SC)” is developed. The MM-SC model is based on the fixed-angle crack formulation and the smeared-crack formulation, in which the experimental-based uniaxial constitutive models are implemented, considering the local buckling of faceplate, the tension stiffening of steel plate, the strength degradation of cracked concrete and the confinement effect of concrete. The proposed MM-SC model is subsequently incorporated into the object-oriented software OpenSEES. Finally, the simulation results of proposed model well predict the SC test observations in terms of critical branch points and structural behaviors, including initial stiffness, cracking strength, post-crack stiffness, yield strength, maximum strength, and failure modes.

42 ENGINEERING↗
DOE OSTI2025

Dispatch Manager for NEML2 Constitutive Model Calculations Embedded in MOOSE

This report describes the extended capabilities of the NEML2 constitutive modeling library, including a flexible and efficient work dispatching system designed to leverage both CPU and GPU resources. This enhancement addresses one of the primary computational challenges in large-scale simulations: the ability to distribute and execute batches of material model evaluations across heterogeneous computing devices. The new dispatch system introduces a modular set of dispatcher and scheduler classes that coordinate the flow of data and execution between devices. The dispatcher is responsible for efficiently packaging work, managing device-specific memory operations, and synchronizing results. This modularity allows for extensibility, making it straightforward to integrate additional computing backends in the future. From an implementation standpoint, the dispatcher system interfaces seamlessly with NEML2's existing models. They handle device-aware tensor operations, optimize memory transfers, and support asynchronous execution when applicable. This design ensures that batches of material points can be evaluated concurrently, substantially improving throughput compared to previous single-device or serial implementations. These improvements not only enhance the raw performance of NEML2 but also improve its usability in multiscale and high-fidelity simulations, where the simultaneous evaluation of large material point batches is critical. Benchmarks included in the report demonstrate the system’s scalability, highlighting its effectiveness when leveraging modern GPU architectures.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Progress towards a fully 3D, thermodynamically-consistent, constitutive model for volume changes during crystallization and melting of semicrystalline polymers

A thermodynamically-consistent constitutive model for semicrystalline polymers has been developed that is capable of representing volume strains from crystallization and melting. The new model is targeted at representing the thermal strain behavior of polymeric photo-voltaic module polymer encapsulation which may undergo crystallization and melting in their service environments. In the model, a homogenized material point is represented by a Gibbs free energy incorporating contributions from the crystalline phase, the amorphous phase, and mixing terms.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A constitutive model incorporating grain refinement strengthening on metallic alloys

Surface nano-crystallization techniques have been recently developed as one of the most efficient ways to optimize materials’ structure, and therefore develop the local and global mechanical behavior as to increase strength without compromising ductility. Here, in this work, we present a constitutive model incorporating grain refinement hardening to simulate the nano-crystallization technique, specifically, surface mechanical attrition treatment. The computation is implemented using user-defined VUMAT subroutines. As an example of its application, a geometry model with full coverage of random impacts are employed. The results show that the model has rather precise predictability of grain size evolution during plastic deformation. The readily embedded with a computational code of material dynamics enables this novel model to be a promising tool to study the dynamic evolution of microstructures under plastic deformation.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Towards data-driven constitutive modelling for granular materials via micromechanics-informed deep learning

The analytical description of path-dependent elastic-plastic responses of a granular system is highly complicated because of continuously evolving microstructures and strain localisation within the system undergoing deformation. This study offers an alternative to the current analytical paradigm by developing micromechanics-informed machine-learning based constitutive modelling approaches for granular materials. A set of critical variables associated with the constitutive behaviour of granular materials are identified through an incremental stress-strain relationship analysis. Depending on the strategy to exploit the priori micromechanical knowledge, three different training strategies are explored. The first model uses only the measurable external variables to make stress predictions; the second model utilises a directed graph to link all the external strain sequences and internal microstructural evolution variables into a single prediction model comprised of a series of sub-mappings, and the third model explicitly integrates the physically important non-temporal properties with external strain paths into training through an enhanced Gated Recurrent Unit (GRU). These three models show satisfactory agreement with unseen test specimens based on multi-directional loading cases. The basic features and potential applications of each model are explained. Lastly, the key factors for constitutive training and limitations of the current work are also discussed in detail.

36 MATERIALS SCIENCE↗
DOE OSTI2026

A wide-range thermodynamically consistent constitutive model for dynamic loading of partially saturated porous media

SUMMARY This paper presents a novel thermodynamically consistent constitutive model for partially saturated porous rocks across a wide range of conditions. The material states generated behind the shock wave from an explosive source can vary significantly, ranging from crushed and melted rock near the source to a poroelastic medium in the far field. In the model, rock strength is determined by the effective pressure, which is calculated using two independent equations of state: one for the solid rock and another for the pore fluid. The model accounts for shock-induced liquefaction resulting from fluid pressure buildup in the pore spaces near the explosive source. Simultaneously, it describes the increase in wave propagation speed due to elastic pore contraction in both dry and partially saturated rocks. This model is applied to investigate how fluid saturation affects the amplitude and shape of the generated waves, as well as the residual stress surrounding the cavity formed by spherical explosions.

Vorobiev, Oleg [Lawrence Livermore National Labora↗
DOE OSTI2023

A constitutive model for elastomers tailored by ionic bonds and entanglements

Over the past decade or two, the concept has emerged of using multiple types of weak interactions simultaneously to enhance the mechanical properties of elastomers. These weak interactions include physical entanglements, hydrogen bonds, metal-coordination bonds, dynamic covalent bonds, and ionic bonds. The combination of entanglements and ionic bonding has been minimally explored and is particularly exciting because of the broad application space for polyelectrolytes. In this work, a constitutive model framework is developed to describe the response of elastomers with both ionic bonds and entanglements. Here we formulate a micromechanical model that couples together chain stretching, ionic bond slipping, and entanglement evolution. The ionic bonds provide toughness by enabling plastic deformation in comparison to covalently crosslinked material and add strength compared to a linear polymer. Evolution of the entanglement density is taken as a key mechanism that can govern stiffness, toughness, and self-recovery in elastomers. The model is used to match bulk polyelectrolytes with different fractions of ionic components under a variety of loading histories. The variations in material parameters are then used to help understand the relative importance of different governing mechanisms in the bulk polymers. We show that the theoretical framework can explain our experimental uniaxial tensile experimental results for polyelectrolytes. This model can help to design better material with high stiffness and toughness. We expect that our model can be extended to explain the mechanical behavior of other polyelectrolytes and other soft materials with a wide range of dynamic bonds.

42 ENGINEERING↗
DOE OSTI2025

Modified 316H constitutive model updated with high temperature stress relaxation test data

This report describes the calibration of a new high temperature constitutive model for 316H stainless steel, suitable for use with the ASME Boiler & Pressure Vessel Section III, Division 5, Class A rules for design by inelastic analysis. The model retains the same mathematical form used by the reference model included in Nonmandatory Appendix Z of the Code, but refits the model to an expanded dataset including all the data used to fit the original model plus seven new stress relaxation tests. The addition of these high temperature stress relaxation tests improves the model's accuracy in predicting relaxation at temperatures greater than 700 ⁰C, without compromising the accuracy of the model versus the original calibration data.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Constitutive Model Development for Aging Polymer Encapsulants (ASC P&EM FY2021 L2 Milestone 7836)

This SAND report fulfills the completion requirements for the ASC Physics and Engineering Modeling Level 2 Milestone 7836 during Fiscal Year 2021. The Sandia Simplified potential energy clock (SPEC) non-linear viscoelastic constitutive model was developed to predict a whole host of polymer glass physical behaviors in order to provide a tool to assess the effects of stress on these materials over their lifecycle. Polymer glasses are used extensively in applications such as electronics packaging, where encapsulants and adhesives can be critical to device performance. In this work, the focus is on assessing the performance of the model in predicting material evolution associated with long-term physical aging, an area that the model has not been fully vetted in. These predictions are key to utilizing models to help demonstrate electronics packaging component reliability over decades long service lives, a task that is very costly and time consuming to execute experimentally. The initiating hypothesis for the work was that a model calibration process can be defined that enables confidence in physical aging predictions under ND relevant environments and timescales without sacrificing other predictive capabilities. To test the hypothesis, an extensive suite of calibration and aging data was assembled from a combination of prior work and collaborating projects (Aging and Lifetimes as well as the DoD Joint Munitions Program) for two mission relevant epoxy encapsulants, 828DGEBA/DEA and 828DGEBA/T403. Multiple model calibration processes were developed and evaluated against the entire set of data for each material. A qualitative assessment of each calibration's ability to predict the wide range of aging responses was key to ranking the calibrations against each other. During this evaluation, predictions that were identified as non-physical, i.e., demonstrated something that was qualitatively different than known material behavior, were heavily weighted against the calibration performance. Thus, unphysical predictions for one aspect of aging response could generate a lower overall rating for a calibration process even if that process generated better quantitative predictions for another aspect of aging response. This insurance that all predictions are qualitatively correct is important to the overall aim of utilizing the model to predict residual stress evolution, which will depend on the interplay amongst the different material aging responses. The DSC-focused calibration procedure generated the best all-around aging predictions for both materials, demonstrating material models that can qualitatively predict the whole host of different physical aging responses that have been measured. This step forward in predictive capability comes from an unanticipated source, utilization of calorimetry measurements to specify model parameters. The DSC-focused calibration technique performed better than compression-focused techniques that more heavily weigh measurements more closely related to the structural responses to be predicted. Indeed, the DSC-focused calibration procedure was only possible due to recent incorporation of the enthalpy and heat capacity features into SPEC that was newly verified during this L2 milestone. Fundamentally similar aspects of the two material model calibrations as well as parametric studies to assess sensitives of the aging predictions are discussed within the report. A perspective on the next steps to the overall goal of residual stress evolution predictions under stockpile conditions closes the report.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Experimental Determination of Constitutive Response of Alloy 800H in Support of Viscoplastic Constitutive Model Development for ASME Section III, Division 5, Class A Applications

In FY 2022, a development effort was initiated at the US Department of Energy’s Argonne National Laboratory (Argonne) and Oak Ridge National Laboratory (ORNL) to develop the Alloy 800H viscoplastic constitutive model so that inelastic analysis methods can support high-temperature gas-cooled reactor (HTGR) designs. The viscoplastic constitutive equations will be developed at Argonne. The ORNL effort focuses on generating test data from various high-temperature experiments with specially designed loading histories to probe Alloy 800H’s constitutive response.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Reference constitutive model for Alloy 617 and 316H stainless steel for use with the ASME Division 5 design by inelastic analysis rules

This report describes the development, calibration, and validation of an inelastic constitutive model describing the monotonic, cyclic, and creep behavior of Alloy 617 from room temperature to 983 °C . The model is intended for incorporation into a new Nonmandatory Appendix to Section III, Division 5, Subsection HB, Subpart B of the ASME Boiler & Pressure Vessel Code providing guidance on developing models for use with the ASME design by inelastic analysis provisions as well as a listing of acceptable material models for each Class A material that designers can use without further validation. The report described the development of the A617 model as well as providing an updated draft appendix, ready for ASME ballot, incorporating the A617 and and the previously-developed 316H model. The development of this appendix and the related material models will promote the widespread application of the generally more efficient inelastic analysis rules rules.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗