Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “constitutive models”

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

Incorporating bubble growth volume feedback to improve simulation of the response of a structure containing liquid and gas to sudden energy input

The SNS target module is a stainless-steel structure that contains and directs mercury. The mercury is struck with short, intense proton pulses to create neutrons. The pulses deposit energy and cause loads on the target structure and mercury cavitation. Helium bubbles are introduced into the mercury to reduce loads and cavitation. This gas complicates the prediction of the target module’s physical response. The current method used to simulate the structural response to the proton pulse incorporates a simple but effective approach to approximate mercury cavitation effects. This method cannot account for the changes from injected non-condensable gas. The ability to model the target vessel response is a prerequisite for estimating its life. Known approaches for simulating a bubbly fluid mixture’s response are computationally expensive and impractical for detailed engineering models. A constitutive model is proposed for target mercury with injected gas bubbles, with the intent of providing more accurate mercury–vessel response behavior without the computing costs required to track individual bubbles. The model’s assumptions and a computational model embodiment are described. The bubbly mercury computational model is subjected to test problems to assess its potential utility. Here, the model was able to better predict the effect of gas on wave propagation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Sample problems for Section III, Division 5 design by inelastic analysis of Grade 91 components

This report works through two sample design analyses of representative Grade 91 components using the ASME Boiler & Pressure Vessel Code Section III, Division 5, Appendix HBB-T rules for the strain and deformation and creep-fatigue design criteria using a design by inelastic analysis approach. The component geometries and loadings were selected to span a wide variety of potential applications. The sample components are an endplate geometry, which is fairly representative of typical pressure vessel construction, and a core block, representing a critical component in a heat pipe microreactor, which has a complicated geometry and complex thermal stress history. These sample problems serve several functions. First, they are complete worked examples of how to apply the current design by inelastic analysis rules. Second, they serve as test cases to evaluate the current rules and provide feedback on how they could be optimized to deal with complicated component geometries or to be better automated in modern finite element analysis software. Finally, the analysis here uses an inelastic constitutive model for Grade 91 developed by Argonne National Laboratory and proposed for inclusion in the ASME Boiler & Pressure Vessel Code Section III, Division 5 rules as a reference constitutive model. Using this constitutive model in the sample problem analyses tests the model in the context of the ASME Code rules. The final conclusions of this report are that the material model and the current Code rules are satisfactory, but that the Code rules could be optimized to simplify automation and reduce over conservatism. Future work should address the specific optimization topics discussed in this report.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Multi-fidelity microstructure-induced uncertainty quantification by advanced Monte Carlo methods

Quantifying uncertainty associated with the microstructure variation of a material can be a computationally daunting task, especially when dealing with advanced constitutive models and fine mesh resolutions in the crystal plasticity finite element method (CPFEM). Numerous studies have been conducted regarding the sensitivity of material properties and performance to the mesh resolution and choice of constitutive model. However, a unified approach that accounts for various fidelity parameters, such as mesh resolutions, integration time-steps and constitutive models simultaneously is currently lacking. Further, this paper proposes a novel uncertainty quantification (UQ) approach for computing the properties and performance of homogenized materials using CPFEM, that exploits a hierarchy of approximations with different levels of fidelity. In particular, we illustrate how multi-level sampling methods, such as multi-level Monte Carlo (MLMC) and multi-index Monte Carlo (MIMC), can be applied to assess the impact of variations in the microstructure of polycrystalline materials on the predictions of homogenized materials properties. We show that by adaptively exploiting the fidelity hierarchy, we can significantly reduce the number of microstructures required to reach a certain prescribed accuracy. Finally, we show how our approach can be extended to a multi-fidelity framework, where we allow the underlying constitutive model to be chosen from either a phenomenological plasticity model or a dislocation-density-based model.

36 MATERIALS SCIENCE↗
DOE OSTI2025

LAROMance Grade 91 Model Integration in NEML2

New reactor designs are targeting higher operating temperatures for increased thermal efficiency when compared to the current fleet of light water reactors. Designing structural components for these high temperature environments with reliable long-term operations requires material models that can accurately capture the deformation mechanisms active in these environments. The LAROMance surrogate material models are based on a database of mechanistic crystal plasticity simulations for high-temperature conditions. Inputs to the LAROMance models reflect the microstructural pedigree of the material, like dislocation densities and precipitate contents. Based on the evolution of these microstructural features, the LAROMance model provides the engineering scale constitutive model response. The LAROMance model was recently parameterized for Grade 91, a high temperature alloy. In the present work, the Grade 91 LAROMance model is implemented in the New Material Model Library, version 2 (NEML2). NEML2 provides a modular way to build material models from smaller blocks and was developed to vectorize the material update to efficiently run on modern computational architectures with graphics processing unit accelerators. NEML2 constitutive models can be used in simulations based on the multiphysics object-oriented simulation environment (MOOSE). This report provides details on the implementation of the Grade 91 LAROMance model in NEML2 and its verification of engineering scale finite element simulations in MOOSE.

42 - ENGINEERING↗
DOE OSTI2020

An inelastic homogenization framework for layered materials with planes of weakness

Many geologic materials have a composite structure, in which macroscopic mechanical behavior is determined by the properties, shape, and heterogeneous distribution of individual constituents. In particular, sedimentary rocks commonly exhibit a layered microstructure, with distinct bedding planes that can also form planes of weakness. In this work, we present a homogenization framework for modeling inelastic layered media. The proposed constitutive model allows for distinct micro-constitutive laws for each layer, explicit representation of layer distributions, as well as incorporation of imperfect bonding at the interface between adjacent layers. No a priori assumptions are needed regarding the specific constitutive models used for the layers and interfaces, providing significant modeling flexibility. The overall framework provides a simple and physically-motivated way of defining anisotropic material behavior as an emergent property of the layered microstructure. The model is calibrated using triaxial and true-triaxial experimental data to demonstrate its ability to describe anisotropic deformation and multiple modes of failure.

58 GEOSCIENCES↗
DOE OSTI2023

An extended three-field principle to scale-bridge the granular micromechanics of polymer-bonded particulate materials

Observed phenomena of damage, viscoelasticity, and viscoplasticity in polymer-bonded particulate materials (PBPM), like polymer-bonded explosive (PBX), are accommodated in a scale-bridging constitutive model based on granular micromechanics. An extended three-field variational principle with min - max optimization problem is introduced to homogenize the bonded aggregates’ grain-scale contact mechanics. Similarly to prior applications of a three-field principle, volume constraints are invoked to embed microscale-to-macroscale correspondences for kinematic and static volumetric measures, like mean strain and its energy-conjugate stress invariant. Unlike prior applications, our three-field principle realizes a strain-driven model for effective constitutive behavior due to bonded intergranular contact. In conclusion, the resultant micro-sphere integration-based homogenization is demonstrated in comparison to existing publicly-available measurements to capture the extreme tension/compression constitutive asymmetry exhibited by PBPMs (specifically PBX-9501) during material softening.

42 ENGINEERING↗
DOE OSTI2024

Extreme sparsification of physics-augmented neural networks for interpretable model discovery in mechanics

Data-driven constitutive modeling with neural networks has received increased interest in recent years due to its ability to easily incorporate physical and mechanistic constraints and to overcome the challenging and time-consuming task of formulating phenomenological constitutive laws that can accurately capture the observed material response. However, even though neural network-based constitutive laws have been shown to generalize proficiently, the generated representations are not easily interpretable due to their high number of trainable parameters. Sparse regression approaches exist that allow for obtaining interpretable expressions, but the user is tasked with creating a library of model forms which by construction limits their expressiveness to the functional forms provided in the libraries. Here, in this work, we propose to train regularized physics-augmented neural network-based constitutive models utilizing a smoothed version of $L^0$-regularization. This aims to maintain the trustworthiness inherited by the physical constraints, but also enables interpretability which has not been possible thus far on any type of machine learning-based constitutive model where model forms were not assumed a priori but were actually discovered. During the training process, the network simultaneously fits the training data and penalizes the number of active parameters, while also ensuring constitutive constraints such as thermodynamic consistency. We show that the method can reliably obtain interpretable and trustworthy constitutive models for compressible and incompressible hyperelasticity, yield functions, and hardening models for elastoplasticity, using synthetic and experimental data. This work aims to set a new paradigm for interpretable machine learning models in the broad area of solid mechanics where low and limited data is available along with prior knowledge of physical constraints that the learned maps need to obey. This paradigm can potentially be extended to a broader spectrum of scientific exploration.

Data-driven constitutive models↗
DOE OSTI2023

A statistical mechanics framework for polymer chain scission, based on the concepts of distorted bond potential and asymptotic matching

To design increasingly tough, resilient, and fatigue-resistant elastomers and hydrogels, the relationship between controllable network parameters at the molecular level (bond type, non-uniform chain length, entanglement density, etc.) to macroscopic quantities that govern damage and failure must be established. Many of the most successful constitutive models for elastomers have been rooted in statistical mechanical treatments of polymer chains. Typically, such constitutive models have used variants of the freely jointed chain model with rigid links. However, since the free energy state of a polymer chain is dominated by enthalpic bond distortion effects as the chain approaches its rupture point, bond extensibility ought to be accounted for if the model is intended to capture chain rupture. To that end, a new bond potential is supplemented to the freely jointed chain model (as derived in the u FJC framework of Buche and Silberstein (2021) and Buche et al. (2022)), which we have extended to yield a tractable, closed-form model of single chain behavior that should be amenable to continuum-level constitutive model development. Inspired by the asymptotically matched u FJC model response in both the low/intermediate chain force and high chain force regimes, a simple, quasi-polynomial bond potential energy function is derived. This bond potential exhibits harmonic behavior near the equilibrium state and anharmonic behavior for large bond stretches tending to a characteristic energy plateau (akin to the Lennard-Jones and Morse bond potentials). Using this bond potential, approximate yet highly-accurate analytical functions for bond stretch and chain force dependent upon chain stretch are established. Then, using this polymer chain model, a stochastic thermal fluctuation-driven chain rupture framework is developed. This framework is based upon a force-modified tilted bond potential that accounts for distortional bond potential energy, allowing for the derivation and subsequent calculation of the dissipated chain scission energy. Here, the cases of rate-dependent and rate-independent scission are accounted for throughout the rupture framework. The impact of Kuhn segment number on chain rupture behavior is also investigated. The model is fit to single chain mechanical response data collected from atomic force microscopy tensile tests for validation and to glean deeper insight into the molecular physics taking place. Due to their analytical nature, this polymer chain model and the associated rupture framework can in the future be implemented in finite element models accounting for fracture and fatigue in polydisperse elastomer networks.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Hydro-Code Implementation and Testing of a Kinetic Phase Transition Framework

In this report we describe the Kinetic Phase Transition (KPT) framework that has been worked out over the last 10 years (from around 2014) and the implementation of it into three different codes, the one-dimensional hydro- LASLO and the three-dimensional magneto-hydro- ALEGRA, Sandia codes, via subroutines in the LAMBDA Equations of State and constitutive models package, and Flag, an arbitrary Lagrangian-Eulerian multiphysics code developed within the Lagrangian Applications project (LAP) at LANL. We discuss the introduction of phase mass (and/or volume) fractions that are needed in a code for it to be ‘phase aware’, that is, not only the thermodynamic state is known in each point but also the mixture of the materials’ phases in that point. Further we point to the need of a full Equations of State for each phase in a material to achieve phase awareness and we review the equilibrium phase model, where a phase mixture is at its lowest Gibbs free energy state, to make this point clear. Contrasting the kinetic phase transition to this equilibrium model seamlessly introduce us to the KPT framework that is subsequently thoroughly discussed. While the determination of the total state and the states and mass fractions of phases in each point is a problem that can borrow many of its numerical details from Eulerian codes and mixture of materials (not phases), the update of mass fractions with time in a KPT framework needs a new set of considerations. General for any update model is that we need to prevent mass fractions from becoming unphysical (negative or their sum to be larger than one). We have solved this problem by implementing a subdivision of the hydro time step that prevents the phase from being fully present to not present at all in one subdivided time step by limiting the size of the subdivided time step. This scheme also corrects numerical problems from abrupt changes in parameter values, the so called Gibbs phenomena, that gives rise to slushing between phases in the KPT framework. Interspersed throughout the report are discussions on different thermodynamics considerations. EOS validity windows, limitations on the EOS phase space, are needed for the KPT framework and are discussed separately and exemplified. The KPT framework described in this report has been verified by code comparison, but validation is still an active area of research. There is room for improvement in the update model, both in the model for determination of rates and in how to prevent the mass fractions from becoming unphysical. In addition, the parameters in the KPT update model and the placement of the phase boundary in the EOS phase space, and interactions with other constitutive models, are closely related and interfering with each other. One possible way forward is to simultaneously develop KPT parameters, EOS, and constitutive models for each material.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

A continuum damage mechanics (CDM) based Wilshire model for creep deformation, damage, and rupture prediction

In this report, a new constitutive model is derived that combines the Wilshire equations with continuum damage mechanics (CDM) to enable the long-term prediction of creep deformation, damage, and rupture. The Wilshire equations have been demonstrated to accurately extrapolate the stress-rupture, minimum-creep-strain-rate, and time-to-strain of various alloys across decade of life. Recently, the time-to-creep-strain equation has been rearranged to predict creep deformation curves; however, the resulting equation is difficult to calibrate and does not track damage. The CDM-based Sinh constitutive model accurately predicts creep deformation, damage, and stress-rupture for a variety of alloys; however, it lacks an explicit description of stress and temperature co-dependence which limits interpolation and extrapolation ability. In this study, the Wilshire and the CDM-based Sinh equations are combined to create the “WCS” model. The WCS model combines the best features of each model while eliminating their deficiencies. Experimental data for alloy P91 is gathered and the WCS model is calibrated to the data. The WCS model accurately predicts the stress rupture, minimum-creep-strain-rate, creep deformation, and damage of alloy P91 across the experimental stress and temperature range as well as consistency in parametric simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Robust deep learning framework for constitutive relations modeling

Modeling the full-range deformation behaviors of materials under complex loading and materials conditions is a significant challenge for constitutive relations (CRs) modeling. Here, we propose a general encoder-decoder deep learning framework that can model high-dimensional stress-strain data and complex loading histories with robustness and universal capability. The framework employs an encoder to project high-dimensional input information (e.g., loading history, loading conditions, and materials information) to a lower-dimensional hidden space and a decoder to map the hidden representation to the stress of interest. We evaluated various encoder architectures, including gated recurrent unit (GRU), GRU with attention, temporal convolutional network (TCN), and the Transformer encoder, on two complex stress-strain datasets that were designed to include a wide range of complex loading histories and loading conditions. All architectures achieved excellent test results with an root-mean-square error (RMSE) below 1 MPa. Additionally, we analyzed the capability of the different architectures to make predictions on out-of-domain applications, with an uncertainty estimation based on deep ensembles. The proposed approach provides a robust alternative to empirical/semi-empirical models for CRs modeling, offering the potential for more accurate and efficient materials design and optimization.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Modular machine learning-based elastoplasticity: Generalization in the context of limited data

The development of highly accurate constitutive models for materials that undergo path-dependent processes continues to be a complex challenge in computational solid mechanics. Challenges arise both in considering the appropriate model assumptions and from the viewpoint of data availability, verification, and validation. Recently, data-driven modeling approaches have been proposed that aim to establish stress-evolution laws that avoid user-chosen functional forms by relying on machine learning representations and algorithms. However, these approaches not only require a significant amount of data but also need data that probes the full stress space with a variety of complex loading paths. Furthermore, they rarely enforce all necessary thermodynamic principles as hard constraints. Hence, they are in particular not suitable for low-data or limited-data regimes, where the first arises from the cost of obtaining the data and the latter from the experimental limitations of obtaining labeled data, which is commonly the case in engineering applications. In this work, we discuss a hybrid framework that can work on a variable amount of data by relying on the modularity of the elastoplasticity formulation where each component of the model can be chosen to be either a classical phenomenological or a data-driven model depending on the amount of available information and the complexity of the response. The method is tested on synthetic uniaxial data coming from simulations as well as cyclic experimental data for structural materials. The discovered material models are found to not only interpolate well but also allow for accurate extrapolation in a thermodynamically consistent manner far outside the domain of the training data. This ability to extrapolate from limited data was the main reason for the early and continued success of phenomenological models and the main shortcoming in machine learning-enabled constitutive modeling approaches. Training aspects and details of the implementation of these models into Finite Element simulations are discussed and analyzed.

42 ENGINEERING↗
DOE OSTI2022

Effects of internal swelling on residual elasticity of a quasi-brittle material through a composite sphere model

This work describes the development of a micromechanical-based constitutive model accounting for the effect of internal expansion on the residual elasticity of a Hashin composite material. This material is made of spherical inclusions that are subjected to gradual swelling within a quasi-brittle matrix. The main focus of this work is to describe and analyze the material mechanical response, with an additional focus on the internal swelling’s effect on the stress–strain response and residual elasticity. Microstructural features and parameters of major importance for the mechanical responses were identified. The innovative characteristics of the proposed approach are summarized as follows: (1) a full determination of the physics of a complete-damage problem throughout the whole process of inclusion swelling with upscaling techniques, which transfers the microcrack-related properties from the lower scale to upper scale; and (2) an evolution of the mechanical fields and the corresponding residual elasticity for various inclusion swelling levels. Concerning the matrix–inclusion composite, it was hypothesized that only the matrix was susceptible to cracking, with varied degrees of damage, whereas the inclusions behave elastically and the elastic modulus of the expanding inclusions remains constant. It is to emphasize that the gradual swelling of inclusions is modeled by an increasing strain in the current micromechanical-based constitutive model, and the microcracks are represented by a set of randomly oriented penny-shaped microcracks with identical radii (namely, closed cracks). The main contribution of the current research is to establish the exact mathematical solutions for the mechanical fields (stress, strain) caused by the swelling of the inclusions (mechanical loading) and derive the effective residual elasticity of a composite (structural response) subject to internal expansion and quasi-brittle damage. Based on the assumption of closed cracks that could be extended to open cracks in the upcoming work, the results proposed by the present paper help to understand the non-linear mechanical behavior of quasi-brittle materials subject to microcracking and provide a theoretical framework to be used as academical benchmark for numerical simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Stabilized Hyperfoam Modeling of the General Plastics EF4003 (3 PCF) Flexible Foam

Constitutive model parameterizations for the General Plastics EF4003 low density 3 pound per cubic foot are needed for design and qualification purposes in normal and abnormal mechanical simulations. The material is expected to be deformed in two ways: first during preloading, and second under impact conditions of the system (transient dynamic). All analyses are to be performed at room temperature. The goal is to provide the analysis community a robust constitutive model parameterization to represent the compression behavior of the EF4003 foam from small deformations up to massive compressive deformations when the foam is densifying. It is worth noting the EF4003 exhibits anisotropy in its stress-strain behavior between the rise and transverse directions (See figure 2.8c-d) as well as plateau behavior that is very likely to cause material stability issues, due to the buckling transition, (and has historically done so) when using Sandia’s current workhorse models for flexible foams, Hyperfoam and Flex Foam. A Stability-informed Hyperfoam parameterization procedure is developed and executed to calibrate a hyperfoam model for the EF4003 room temperature, transversely loaded data. A rise orientation parameterization was not attempted due to localization in the experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Mechanistic Modeling of Creep and Aging for Time Extrapolation in LPBF 316 SS - LANL

A polycrystal constitutive model was extended to quantify the effects heterogeneity in the microstructure, porosity, and internal stresses on the mechanical response (creep, tensile) of 316H steel processed by additive manufacturing. The starting point is a recently developed and fully calibrated advanced microstructuresensitive mechanistic constitutive model fit to literature data pertaining to conventionally processed steels. The model is embedded within the full field fast Fourier transform framework (LApx) and several hundred simulations were performed. Synthetic microstructures with different grain aspect ratios were generated using the Neper package to study the role of grain morphology on the mechanical responses. The effect of initial cell interior dislocation density, cell-wall density, precipitate content, porosity, and internal stresses on creep rates are quantified for different polycrystalline microstructures. The study allowed to prospectively quantify the role of distinct features of the microstructure on the mechanical response under tension and creep loading scenarios. Critically it is found that minor anisotropies in the tensile response of samples, emanating from the aspect ratio of grains correlated into significant differences in creep rate as a function of loading directions. The study also allows to identify initial microstructure porosities below which the materials response only moderately degraded by the initial microstructure.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A Review of Computational Models for the Flow of Milled Biomass Part II: Continuum-Mechanics Models

The design of efficient material-handling systems for milled lignocellulosic biomass is challenging due to their complex particle morphologies and frictional interactions. Computational modeling, including the discrete element method (DEM) and continuum-based finite-element/volume methods, may offer scientific insight and predictive capabilities for the flow of milled biomass in hoppers and feeders. Herein, this article (Part II) presents a review of current state-of-the-art continuum models for the flow of milled biomass, whereas DEM models are reviewed in a companion article (Part I). Advances of numerical methods to solve the global governing equations are discussed first, followed by a comprehensive review of constitutive models for granular materials, including Drucker–Prager, hypoplastic, Cambridge-type, inertial-rheology, and nonlocal granular fluidity models. Specifically, we provide in-depth discussion on the suitability of those models for milled lignocellulosic biomass materials in terms of nonlinear elasticity, dependence of flow strength on pressure, density and shear rate, and compaction (dilation) associated with hardening (softening). Furthermore, our study shows that, despite the recent advances in continuum granular flow modeling, the most suitable constitutive models still need further development to account for material parametrization, multiflow regimes, and multiscale behavior before they can be reliably used to optimize the design and operation of biomass handling systems.

09 BIOMASS FUELS↗
DOE OSTI2020

Multiscale-Informed Modeling of High Temperature Component Response with Uncertainty Quantification

This report summarizes a joint effort between Argonne National Laboratory, Idaho National Laboratory, and Los Alamos National Laboratory to develop and deploy constitutive models targeted at predicting the life of Grade 91 alloy components subjected to high temperature environments typical of those that structural components in advanced nuclear reactors would experience. Two distinct, but complementary constitutive modeling approaches have been taken here. The first employs a phenomenological viscoplastic model for which parameters have been calibrated based on experimental data for a wide range of Grade 91 alloy that has undergone a variety of processing. A Bayesian approach was used to derive distributions of uncertain parameters for this model based on this data set. The second approach is a reduced order model suitable for engineering-scale analysis that is based on the results of a large set of mesoscale simulations. Mesoscale models allow for the microstructure and composition of a particular alloy to be directly taken into account in the computation of the viscoplastic response, but are computationally expensive, which makes it impractical to directly call those models for the material constitutive response in an engineering-scale simulation. The reduced-order representation of the response of the underlying model used here allows for an engineering-scale model to take into account the characteristics of the underlying microstructure, while only incurring a reasonable computational expense. Both of these approaches have different strengths, and are applicable for different parts of the design/analysis process. The phenomenological models can be readily parameterized based on a set of experimental data for a given class of materials and used for scoping calculations. Once a specific material is chosen and adequately characterized, the reduced order models can accurately predict the response of that specific alloy, and because the models are based on predictive models of the underlying microstructure, they can be used to more confidently predict the response under conditions in regions where there is limited experimental data. Both of these models have been integrated in the Grizzly code, which is used here to perform proof-of-concept uncertainty quantification analyses of a simple component under prototypical conditions. The built- in stochastic analysis capabilities in the MOOSE framework that Grizzly is built on are used here to run large sets of simulations for this uncertainty quantification analysis. As would be expected, because the reduced order models are developed for a much more tightly defined alloy, they predict tighter distributions of the time to failure than the phenomenological models, which are calibrated to a broader set of data. Also important is that these simulations demonstrate that a reduced order modeling approach can be successfully deployed to propagate uncertainties from the material scale to practical engineering-scale component simulations.

42 ENGINEERING↗
DOE OSTI2025

A direct-adjoint approach for material point model calibration with application to plasticity

Here, this paper proposes a new approach for the calibration of material parameters in local elastoplastic constitutive models. The calibration is posed as a constrained optimization problem, where the constitutive model evolution equations for a single material point serve as constraints. The objective function quantifies the mismatch between the stress predicted by the model and corresponding experimental measurements. To improve calibration efficiency, a novel direct-adjoint approach is presented to compute the Hessian of the objective function, which enables the use of second-order optimization algorithms. Automatic differentiation is used for gradient and Hessian computations. Two numerical examples are employed to validate the Hessian matrices and to demonstrate that the Newton–Raphson algorithm consistently outperforms gradient-based algorithms such as L-BFGS-B.

36 MATERIALS SCIENCE↗