SEARCH · Engineering Papers

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 361 records · Page 20

All sources NASA NTRS DOE OSTI

NASA NTRS2003

An Approximate Dissipation Function for Large Strain Rubber Thermo-Mechanical Analyses

Mechanically induced viscoelastic dissipation is difficult to compute. When the constitutive model is defined by history integrals, the formula for dissipation is a double convolution integral. Since double convolution integrals are difficult to approximate, coupled thermo-mechanical analyses of highly viscous rubber-like materials cannot be made with most commercial finite element software. In this study, we present a method to approximate the dissipation for history integral constitutive models that represent Maxwell-like materials without approximating the double convolution integral. The method requires that the total stress can be separated into elastic and viscous components, and that the relaxation form of the constitutive law is defined with a Prony series. Numerical data is provided to demonstrate the limitations of this approximate method for determining dissipation. Rubber cylinders with imbedded steel disks and with an imbedded steel ball are dynamically loaded, and the nonuniform heating within the cylinders is computed.

Johnson, Arthur R.↗

NASA NTRS1988

Life prediction of thermomechanical fatigue using total strain version of strainrange partitioning (SRP): A proposal

A method is proposed (without experimental verification) for extending the total strain version of Strainrange Partitioning (TS-SRP) to predict the lives of thermomechanical fatigue (TMF) cycles. The principal feature of TS SRP is the determination of the time-temperature-waveshape dependent elastic strainrange versus life lines that are added subsequently to the classical inelastic strainrange versus life lines to form the total strainrange versus life relations. The procedure is based on a derived relation between failure and flow behavior. Failure behavior is represented by conventional SRP inelastic strainrange versus cyclic life relations, while flow behavior is captured in terms of the cyclic stress-strain response characteristics. Stress-strain response is calculated from simple equations developed from approximations to more complex cyclic constitutive models. For applications to TMF life prediction, a new testing technique, bithermal cycling, is proposed as a means for generating the inelastic strainrange versus life relations. Flow relations for use in predicting TMF lives would normally be obtained from approximations to complex thermomechanical constitutive models. Bithermal flow testing is also proposed as an alternative to thermomechanical flow testing at low strainranges where the hysteresis loop is difficult to analyze.

Saltsman, James F.↗

NASA NTRS1986

Estimation of high temperature low cycle fatigue on the basis of inelastic strain and strainrate

Fatigue life at elevated temperature can be predicted by introducing parametric values obtained from monotonic constitutive behavior into the Universal-Slopes Equation. For directionally solidified MAR-M200+HF at 975 C, these parameters are the maximum stress achievable under entirely plastic (time-independent) and purely creep (time-dependent) conditions and the corresponding inelastic strains, as well as the elastic modulus. For materials which exhibit plasticity/creep interaction, two more pairs of monotonic parameters must be evaluated for fatigue life prediction. This life-prediction method based on the Universal-Slopes Equation, resulted from a constitutive model characterizing monotonic and cyclic data as inelastic strainrate as a function of inelastic strain. Characterizing monotonic data is this way, permitted distinction between different material responses such as strain-hardening, strain-softening, and dynamic recovery effects. Understanding and defining the region of influence of each of these effects facilitated formulation of the constitutive model in relation to the mechanical and microstructural processes occurring in the material under cyclic loading.

Berkovits, A.↗

NASA NTRS2013

PowderSim: Lagrangian Discrete and Mesh-Free Continuum Simulation Code for Cohesive Soils

PowderSim is a calculation tool that combines a discrete-element method (DEM) module, including calibrated interparticle-interaction relationships, with a mesh-free, continuum, SPH (smoothed-particle hydrodynamics) based module that utilizes enhanced, calibrated, constitutive models capable of mimicking both large deformations and the flow behavior of regolith simulants and lunar regolith under conditions anticipated during in situ resource utilization (ISRU) operations. The major innovation introduced in PowderSim is to use a mesh-free method (SPH-based) with a calibrated and slightly modified critical-state soil mechanics constitutive model to extend the ability of the simulation tool to also address full-scale engineering systems in the continuum sense. The PowderSim software maintains the ability to address particle-scale problems, like size segregation, in selected regions with a traditional DEM module, which has improved contact physics and electrostatic interaction models.

Johnson, Scott↗

NASA NTRS1993

Evolution of damage and plasticity in titanium-based, fiber-reinforced composites

The inelastic deformation mechanisms were evaluated for a model titanium-based, fiber-reinforced composite: a beta titanium alloy (Ti-15V-3Al-3Cr-3Sn) reinforced with SiC (SCS-6) fibers. The primary emphasis of this article is to illustrate the sequence in which damage and plasticity evolved for this system. The mechanical responses and the results of detailed microstructural evaluations for the 0(8), 90(8), and +/- 45(2s) line oriented laminates are provided. It is shown that the characteristics of the reaction zone around the fiber play a very important role in the way damage and plasticity evolve, particularly in the microyield regime of deformation, and must be included in any realistic constitutive model. Fiber-matrix debonding was a major damage mode for the off-axis systems. The tension test results are also compared with the predictions of a few constitutive models.

Majumdar, B. S.↗

NASA NTRS1989

Predictions of Poisson's ratio in cross-ply laminates containing matrix cracks and delaminations

A damage-dependent constitutive model for laminated composites has been developed for the combined damage modes of matrix cracks and delaminations. The model is based on the concept of continuum damage mechanics and uses second-order tensor valued internal state variables to represent each mode of damage. The internal state variables are defined as the local volume average of the relative crack face displacements. Since the local volume for delaminations is specified at the laminate level, the constitutive model takes the form of laminate analysis equations modified by the internal state variables. Model implementation is demonstrated for the laminate engineering modulus E(x) and Poisson's ratio nu(xy) of quasi-isotropic and cross-ply laminates. The model predictions are in close agreement to experimental results obtained for graphite/epoxy laminates.

Harris, Charles E.↗

NASA NTRS1989

Computer simulation of the mathematical modeling involved in constitutive equation development: Via symbolic computations

Development of new material models for describing the high temperature constitutive behavior of real materials represents an important area of research in engineering disciplines. Derivation of mathematical expressions (constitutive equations) which describe this high temperature material behavior can be quite time consuming, involved and error prone; thus intelligent application of symbolic systems to facilitate this tedious process can be of significant benefit. A computerized procedure (SDICE) capable of efficiently deriving potential based constitutive models, in analytical form is presented. This package, running under MACSYMA, has the following features: partial differentiation, tensor computations, automatic grouping and labeling of common factors, expression substitution and simplification, back substitution of invariant and tensorial relations and a relational data base. Also limited aspects of invariant theory were incorporated into SDICE due to the utilization of potentials as a starting point and the desire for these potentials to be frame invariant (objective). Finally not only calculation of flow and/or evolutionary laws were accomplished but also the determination of history independent nonphysical coefficients in terms of physically measurable parameters, e.g., Young's modulus, was achieved. The uniqueness of SDICE resides in its ability to manipulate expressions in a general yet predefined order and simplify expressions so as to limit expression growth. Results are displayed when applicable utilizing index notation.

Arnold, S. M.↗

DOE OSTI2024

Dynamic Modeling and Characterization of Nuclear-grade Graphite

Idaho National Labs serves as the spearhead for many innovative energy solutions to the world's energy crisis. One such solution is the INL's Microreactor which is designed to deploy to extreme/remote environments where other sources of power are either unavailable or unreliable. In order to best design these energy solutions for their operational environments, it is crucial to understand how the design, components, and materials will respond to the environmental conditions. One key material in these innovative designs is a nuclear-grade graphite known as PCEA. This study examines the behavior of PCEA graphite under dynamic loading, similar to that which may occur in extreme environments. The objective is to characterize the dynamic behavior and produce an accurate, reliable constitutive model suitable for use in simulation tools such as INL's MOOSE. Graphite specimens were tested using a Split Hopkinson Pressure Bar (SHPB) to administer the dynamic compressive load. The SHPB was charged at various pressures to produce a range of strain rates on the material in compression. Data was acquired via strain gauges on the SHPB setup, from which stress, strain, and time data were collected. Analysis revealed the stress-strain behavior of the material as well as insights into the material behavior's relationship to strain rate. Further work must continue to characterize the various other dynamic behaviors of the material which will combine to create a substantially trustworthy constitutive model for this grade of nuclear-grade graphite. Ultimately, this will allow for realistic simulation of the material in reactor designs, allowing for prediction of design weaknesses and leading to improved designs for increased resilience, security, and reliability.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS1988

Constitutive response of Rene 80 under thermal mechanical loads

The applicability of a classical constitutive model for stress-strain analysis of a nickel base superalloy, Rene' 80, in the gas turbine thermomechanical fatigue (TMF) environment is examined. A variety of tests were conducted to generate basic material data and to investigate the material response under cyclic thermomechanical loading. Isothermal stress-strain data were acquired at a variety of strain rates over the TMF temperature range. Creep curves were examined at 2 temperature ranges, 871 to 982 C and 760 to 871 C. The results provide optimism on the ability of the classical constitutive model for high temperature applications.

Kim, K. S.↗

DOE OSTI2024

KOMPASS-II: Compaction of Crushed salt for Safe Containment – Phase 2

Long-term stable sealing elements are a basic component in the safety concept for a possible repository for heat-emitting radioactive waste in rock salt. The sealing elements will be part of the closure concept for drifts and shafts. They will be made from a welldefinied crushed salt in employ a specific manufacturing process. The use of crushed salt as geotechnical barrier as required by the German Site Selection Act from 2017 /STA 17/ represents a paradigm change in the safety function of crushed salt, since this material was formerly only considered as stabilizing backfill for the host rock. The demonstration of the long-term stability and impermeability of crushed salt is crucial for its use as a geotechnical barrier. The KOMPASS-II project, is a follow-up of the KOMPASS-I project and continues the work with focus on improving the understanding of the thermal-hydraulic-mechanical (THM) coupled processes in crushed salt compaction with the objective to enhance the scientific competence for using crushed salt for the long-term isolation of high-level nuclear waste within rock salt repositories. The project strives for an adequate characterization of the compaction process and the essential influencing parameters, as well as a robust and reliable long-term prognosis using validated constitutive models. For this purpose, experimental studies on long-term compaction tests are combined with microstructural investigations and numerical modeling. The long-term compaction tests in this project focused on the effect of mean stress, deviatoric stress and temperature on the compaction behavior of crushed salt. A laboratory benchmark was performed identifying a variability in compaction behavior. Microstructural investigations were executed with the objective to characterize the influence of pre-compaction procedure, humidity content and grain size/grain size distribution on the overall compaction process of crushed salt with respect to the deformation mechanisms. The created database was used for benchmark calculations aiming for improvement and optimization of a large number of constitutive models available for crushed salt. The models were calibrated, and the improvement process was made visible applying the virtual demonstrator.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2026

Improving the Performance of NEML2 with Modern Graph Compilation Backends

NEML2 vectorizes constitutive-model evaluation for large-scale multiphysics simulation, using PyTorch as its tensor backend so that a batch of material-point updates runs on CPU or GPU through a single implementation. In the two prior reports in this series it was a C++-native library, deployed through TorchScript tracing and just-in-time (JIT) compilation; it has since been rewritten from the ground up into a Python-native library deployed through Ahead-of-Time Inductor (AOTInductor), a modern PyTorch graph-compilation backend. The rewrite is driven by a persistent tension, not a language preference: NEML2 composes constitutive models at runtime from a registry of small, independently-authored pieces, and that flexibility is difficult to reconcile with the compile-time knowledge an efficient GPU kernel needs. This report documents the rewrite and the investment that accompanied it: the AOTInductor export pipeline that turns a Python-authored model into a portable, Python-free compiled artifact loadable from pure C++; the eager and compiled runtimes and the new implicit solver layer built on them; a head-to-head benchmark of legacy JIT against AOTInductor; the physics-model catalog and its worked examples; the developer tooling; and the corresponding overhaul of MOOSE’s NEML2 integration that lets MOOSE consume it. A central objective is to examine whether modern PyTorch graph-compilation backends are effective for MOOSE GPU integration. The benchmark answers directly: AOTInductor outperforms legacy JIT on every GPU scenario measured, by 1.0–4.5×. Modern graph-compilation backends are effective for MOOSE GPU integration, and AOTInductor specifically – not compilation in the abstract – is why.

Hu, Gary (Tianchen) [Argonne National Laboratory (↗

DOE OSTI2024

Enabling Parallel Performance and Portability of Solid Mechanics Simulations Across CPU and GPU Architectures

Efficiently simulating solid mechanics is vital across various engineering applications. As constitutive models grow more complex and simulations scale up in size, harnessing the capabilities of modern computer architectures has become essential for achieving timely results. This paper presents advancements in running parallel simulations of solid mechanics on multi-core CPUs and GPUs using a single-code implementation. This portability is made possible by the C++ matrix and array (MATAR) library, which interfaces with the C++ Kokkos library, enabling the selection of fine-grained parallelism backends (e.g., CUDA, HIP, OpenMP, pthreads, etc.) at compile time. MATAR simplifies the transition from Fortran to C++ and Kokkos, making it easier to modernize legacy solid mechanics codes. We applied this approach to modernize a suite of constitutive models and to demonstrate substantial performance improvements across different computer architectures. This paper includes comparative performance studies using multi-core CPUs along with AMD and NVIDIA GPUs. Results are presented using a hypoelastic–plastic model, a crystal plasticity model, and the viscoplastic self-consistent generalized material model (VPSC-GMM). The results underscore the potential of using the MATAR library and modern computer architectures to accelerate solid mechanics simulations.

Morgan, Nathaniel (ORCID:0000000276118449)↗

DOE OSTI2025

Fully implicit crystal plasticity models representing orientations with modified Rodrigues parameters

Here, this work describes a crystal plasticity formulation combining several mathematical, numerical, and implementation choices to produce a highly efficient model. Specifically, the key choices in the implementation are (1) representing orientations with modified Rodrigues parameters, (2) implementing a fully coupled implicit time integration for the elastic stretch, the crystal orientations, and the model internal variables, (3) implementing the model in the NEML2 constitutive modeling framework, based on PyTorch, to vectorize the calculations and port the computation to GPUs and other hardware accelerators, and (4) an exact implementation of the consistent tangent matrix, even for arbitrary coupling to other field variables beyond the displacements, like temperature, neutron fluence, etc. The first two features of the model are, to our knowledge, novel. The paper considers each of these choices individually as well as the final model as a whole. This includes a full description of modified Rodrigues parameters, their advantages over other representations of orientations, the mathematical formulae and tools required to implement a model with modified Rodrigues parameters, and a detailed description of the geometry of the space of modified Rodrigues parameters (in an appendix). It also includes a description of a fully implicit time integration scheme for the orientations and the advantages in representing orientations with modified Rodrigues parameters in implementing such a model. The work then assess, via numerical examples, the advantages of fully coupled implicit time integration versus more common decoupled and explicit time integration schemes. These studies demonstrate the computational advantages of fully coupled integration versus other time integration algorithms, though the performance of the competing models depends on the complexity of the underlying single crystal model. The study concludes by demonstrating that the choice of time integration method affects the sharpness of the predicted texture, with explicit methods for integrating the orientations overestimating texture sharpness and implicit methods underestimating texture sharpness.

Crystal plasticity↗

NASA NTRS2016

Analysis and Characterization of Damage Utilizing an Orthotropic Generalized Composite Material Model Suitable for Use in Impact Problems

The need for accurate material models to simulate the deformation, damage and failure of polymer matrix composites under impact conditions is becoming critical as these materials are gaining increased usage in the aerospace and automotive communities. In order to address a series of issues identified by the aerospace community as being desirable to include in a next generation composite impact model, an orthotropic, macroscopic constitutive model incorporating both plasticity and damage suitable for implementation within the commercial LS-DYNA computer code is being developed. The plasticity model is based on extending the Tsai-Wu composite failure model into a strain hardening-based orthotropic plasticity model with a non-associative flow rule. The evolution of the yield surface is determined based on tabulated stress-strain curves in the various normal and shear directions and is tracked using the effective plastic strain. To compute the evolution of damage, a strain equivalent semi-coupled formulation is used in which a load in one direction results in a stiffness reduction in multiple material coordinate directions. A detailed analysis is carried out to ensure that the strain equivalence assumption is appropriate for the derived plasticity and damage formulations that are employed in the current model. Procedures to develop the appropriate input curves for the damage model are presented and the process required to develop an appropriate characterization test matrix is discussed.

Impact↗

NASA NTRS2016

Analysis and Characterization of Damage Utilizing an Orthotropic Generalized Composite Material Model Suitable for Use in Impact Problems

Impact↗

NASA NTRS2018

High Strain Rate Response of Adhesively Bonded Fiber-Reinforced Composite Joints – A Computational Study to Guide Experimental Design

Adhesively bonded carbon fiber-reinforced epoxy composite laminates are widely used in aerospace applications. During a high energy impact event, these laminates are often subjected to high strain rate loading. However, the influence of high strain rate loading on the response of these composite joints is not well understood. Computational finite element (FE) modeling and simulations are conducted to guide the design of high strain rate experiments. Two different experimental designs based on split Hopkinson bar were numerically modeled to simulate Mode I and Mode II types loading in the composite. In addition, the computational approach adopted in this study helps in understanding the high strain rate response of adhesively bonded composite joints subjected to nominally Mode I and Mode II loading. The modeling approach consists of a ply-level 3D FE model, a progressive damage constitutive model for the composite material behavior and a cohesive tie-break contact element for interlaminar delamination.

Ravindran, Suraj↗

NASA NTRS1987

Structural analysis demonstration of constitutive and life models

The overall objective of this program is to demonstrate the applicability of NASA-developed advanced constitutive and life damage models for calculating cyclic structural response and crack initiation in selected components of reusable space propulsion systems. The computer model resulting from this program will enable the user to produce an accurate life prediction of hot gas path, life limiting components of propulsion systems such as the space shuttle main engine (SSME). Previously developed computer models addressing constitutive modeling and life damage will be combined in an advanced finite element analysis to generate a sophisticated baseline life prediction program. A material data base will be established for the constitutive and life models parametrically involving temperature, strain range, strain rate, mean strain/stress, and dwell time. The verified computer program will be used to accomplish the life predictions of three SSME critical components as evidence of the model functionality.

Bechtel, G.↗

NASA NTRS1992

A viscoplastic model for single crystals

A viscoplastic constitutive model is described in which deformation behavior is postulated on representative slip systems and the behavior of the entire crystal is determined by summing the slip on the active slip systems. By building in the slip geometry known from the metallurgical literature, it is possible to predict the anisotropic deformation behavior and to model in a straightforward manner other phenomena which have been described by metallurgists in crystallographic terms. Elevated temperature tension-torsion tests were run and used to verify the model's predictive abilities. Ratchetting behavior under thermomechanical loading conditions is specifically addressed.

Jordan, E. H.↗