Search NASA⌕ Search

SEARCH · Search NASA

Results for “viscoplastic”

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

A comparison of surrogate constitutive models for viscoplastic creep simulation of HT-9 steel

Mechanistic microstructure-informed constitutive models for the mechanical response of polycrystals are a cornerstone of computational materials science. However, as these models become increasingly more complex – often involving coupled differential equations describing the effect of specific deformation modes – their associated computational costs can become prohibitive, particularly in optimization or uncertainty quantification tasks that require numerous model evaluations. To address this challenge, surrogate constitutive models that balance accuracy and computational efficiency are highly desirable. Data-driven surrogate models, that learn the constitutive relation directly from data, have emerged as a promising solution. In this work, we develop two local surrogate models for the viscoplastic response of a steel: a piecewise response surface method and a mixture of experts model. These surrogates are designed to adapt to complex material behavior, which may vary with material parameters or operating conditions. The surrogate constitutive models are applied to creep simulations of HT-9 steel, an alloy of considerable interest to the nuclear energy sector due to its high tolerance to radiation damage, using training data generated from viscoplastic self-consistent (VPSC) simulations. In conclusion, we define a set of test metrics to numerically assess the accuracy of our surrogate models for predicting viscoplastic material behavior, and show that the mixture of experts model outperforms the piecewise response surface method in terms of accuracy.

36 MATERIALS SCIENCE↗

Data-Driven Surrogate Modeling with Microstructure-Sensitivity of Viscoplastic Creep in Grade 91 Steel

Abstract To support the development of advanced steel alloys tailored to withstand extreme conditions, it is imperative to account for the mechanical performance of components, while considering the influence of local microstructure on the macroscopic response. To this end, this study focuses on the development of microstructure-sensitive constitutive models for the mechanical response of Grade 91 steel exposed to extreme thermo-mechanical environments. Polynomial chaos expansion (PCE) surrogates are used to emulate high-fidelity polycrystal simulations of the viscoplastic response of Grade 91 steel as a function of the microstructure fingerprint (e.g., dislocations and precipitates). To cover a wide temperature–stress domain, two separate PCE surrogates—one that captures softening and the other that captures hardening behavior—are combined using another (sparse) Gaussian process regression model. The resulting constitutive creep surrogate model is integrated within the MOOSE finite element framework to simulate the intricate effects of microstructure, in particular MX-phase precipitates, on a component with a graded microstructure. Surrogate sensitivity analysis is applied to quantify the relevant impact of spatially varying microstructure on the creep response in a test-case involving a Grade 91 alloy with a prototypical weld.

36 MATERIALS SCIENCE↗

Predicting the viscoplastic response of a crystallizing fluoropolymer using transient network theory

We employ a molecular theory of dynamic polymer networks to describe the viscoplastic response of rubbery FK-800, a thermoplastic copolymer of chlorotrifluoroethylene and vinylidene fluoride, over a broad range of thermal histories. The kinetics of crystallization at different annealing temperatures was modeled using a modified Avrami equation, whose parameters were found to evolve through simple relationships over the full temperature range of the rubbery state. By fitting experimental compression data, we discovered predictable trends for the physical parameters in our mechanical model over its full range of crystallinities (up to ≈20%) and provided insights based on molecular-level physics to justify them. Using this, an end-to-end model was developed to predict the yielding and post-yield behavior of rubbery FK-800 for arbitrary thermal histories. The model successfully predicted the highly nonlinear evolution of characteristic mechanical signatures (stiffness, yield point, post-yield drop) throughout the crystallization process. A statistical analysis of variance test was employed to determine that the measured variations in the mechanical behavior of rubbery FK-800 are primarily dictated by its fractional crystallinity, regardless of its exact thermal history.

36 MATERIALS SCIENCE↗

Elasto-viscoplastic fast Fourier transform modeling framework for assessing microstructural effects on stress intensity factors characterizing fracture toughness

A large-strain elasto-viscoplastic fast Fourier transform (LS-EVPFFT) model with non-periodic (NP) velocity-based boundary conditions is adapted to simulate the sensitivity of stress intensity factors on microstructure for 304L stainless steel. The material was characterized via electron backscattered diffraction (EBSD) serial-sectioning to obtain a measured 3-D microstructural cell to perform simulations. The NP-LS-EVPFFT model, including the simulation setup and boundary conditions, was verified using a crystal plasticity finite element (CPFE) model. To this end, the generation of meshes of notched specimens was developed, which involved creating Python scripts for mesh “cutting” in Abaqus, and Sculpt scripts in Cubit for meshing of the measured microstructural cell processed with DREAM.3D. The complexity of the mesh preparation highlighted the advantages of the FFT-based model, which circumvents the mesh generation process. Given the efficiency of the FFT-based model, statistical distribution of stress intensity factors in function of crystal orientation at the crack tip, grain structure, and crystallographic texture surrounding the crack tip were predicted. Further, the distributions reveal about 10% variation of stress intensity factors with microstructure with the most significant sensitivity found to be the crystal orientation at the crack tip. The methodology developed in this work is discussed as a practical simulation tool for predicting the sensitivity of stress intensity factors on microstructural variability in metallic materials.

36 MATERIALS SCIENCE↗

The Viscoplastic Behavior of Natural Hydrate‐Bearing Sandy‐Silts Under Uniaxial Strain Compression ( K 0 Loading)

Abstract The in‐situ stress state and geomechanical properties of hydrate‐bearing sediments impact hydrate formation and gas production strategies. We explore the uniaxial strain compression and stress evolution of natural hydrate‐bearing sandy‐silts from Green Canyon Block 955 in the deep‐water Gulf of Mexico. We performed constant rate of strain uniaxial strain experiments, interrupted by periods where we held the axial stress constant, to explore the vertical deformation and the evolution of the ratio of lateral to axial effective stress ( K 0 ) with time. The hydrate‐bearing sandy‐silt is stiffer and has a larger K 0 than the equivalent hydrate‐free sediment upon loading. During stress holds, the void ratio decreases sigmoidally with the log of time, and K 0 converges to isotropic conditions. We interpret that during loading, the hydrate bears the load and deforms. With time, the hydrate redistributes the load and K 0 increases. We used a viscoelastic model to describe the behavior. The model accurately captures deformation and K 0 trends but does not reproduce all the complex interactions of the hydrate with the porous skeleton. We anticipate that viscous effects within hydrate sediments will impact reservoir compression and stresses during production (hours to days), result in isotropic stress state over geological timescales, and explain the creeping movement in submarine landslides.

58 GEOSCIENCES↗

Viscoplastic selfconsistent (VPSC) code (version 8, open source)

VPSC is a mean-field polycrystal plasticity code for the prediction of the mechanical response and microstructure evolution of polycrystalline aggregates. The VPSC code is the computational realization of the visco-plastic self-consistent theory of polycrystal deformation. Both the theory and details of its numerical implementation were originally reported by R.A. Lebensohn and C.N. Tomé: "A self-consistent approach for the simulation of plastic deformation and texture development of polycrystals: application to Zirconium alloys", Acta Metallurgica et Materialia 41, 2611 (1993). Since then, the VPSC code has experienced several improvements and extensions and it is nowadays extensively used to simulate plastic deformation of polycrystalline aggregates and to interpret experimental evidence on metals, minerals and polymers.

Lebensohn, Ricardo↗

Preliminary Experimental Results in Support of the Development of EPP and SMT Design Methods and Viscoplastic Model for A709

The ASME code qualification effort for Alloy 709 (A709) is currently underway to qualify it for Class A construction in the American Society of Mechanical Engineers (ASME) Boiler and Pressure Vessel Code, Section III, Division 5. The United States (U.S.) Department of Energy (DOE) national laboratories are collaborating in the advanced materials development initiative to investigate the mechanical performance of A709 in support of its code qualification. As part of the A709 code qualification effort, this report summarizes ORNL’s initial experimental findings that support the integration of A709 into the elastic-perfectly plastic (EPP) strain limits code case. It also covers thermomechanical fatigue testing conducted to develop viscoelastic material models, along with the preliminary results of creep-fatigue experiments at 816°C using the Simplified Model Test (SMT) method.

36 MATERIALS SCIENCE↗

The physics of desiccation cracks 2: Modeling and prediction of the crack patterns

Here, this paper extends the experimental results of our companion paper by modeling and predicting the onset and pattern formation of desiccation cracks in geomaterials (Ruoyu et al., 2023). Thin-layer silt samples in controlled atmospheric conditions were tested to obtain the surface strain maps with the digital image correlation (DIC) method during dehydration. Support experiments, including consolidation and displacement-controlled triaxial tests, were conducted for the properties of geomaterials. These experimental results were used to validate a viscoplastic theoretical model by comparing cracks and singularities locations that distribute following the Cnoidal wave pattern. A critical value exists in the viscoplastic model that determines the number of singularities, which accurately predicts the number of cracks in the experiment. The critical values contain two crucial parameters: rate sensitivity and λ. Rate sensitivity describes the rate-dependent stress–strain relation, while λ shows the ability of pore pressure redistribution under the external mechanical loading rate. These results provide a new view to analyze the desiccation cracks considering the rate-dependent viscoplasticity.

42 ENGINEERING↗

A phase-field fracture formulation for generalized standard materials: The interplay between thermomechanics and damage

Accurately modeling fracture of ductile materials poses open challenges in the field of computational mechanics due to the multiphysics nature of their failure processes. Integrating the interplay between thermodynamics and damage into ductile fracture models is vital for predicting critical failure modes. Here, in this paper, we develop a versatile phase-field (PF) framework for modeling ductile fracture, taking into account finite-strain elasto-plasticity. The framework stems from a variational formulation of constitutive relations for generalized standard materials (GSMs), whose response is described by a Helmholtz free energy and a dissipation pseudo-potential. Its variational structure is based on a minimum principle for a functional that expresses the sum of power densities for reversible and irreversible processes. By minimizing this functional with a constraint on a von Mises yield function, we derive the evolution equation for the equivalent plastic strain and an associative flow rule. This constrained optimization problem is analytically solved for a wide class of thermo-viscoplasticity models. The key innovations of the current work include (i) a cubic plastic degradation function that accounts for a non-vanishing damage-dependent yield stress, (ii) closed-form expressions of the Helmholtz free energy and dissipation pseudo-potential for three thermo-viscoplasticity models, (iii) an extended Johnson–Cook plasticity model with a nonlinear hardening law, and (iv) a plastic work heat source that depends on the plastic degradation function and a variable Taylor–Quinney (TQ) coefficient. The capabilities of the proposed framework are tested with the aid of four ductile fracture problems, including the Sandia Fracture Challenge. In each of these problems, we examine the evolution of relevant field variables such as the PF order parameter, the equivalent plastic strain, the temperature, and the internal power dissipation density, in addition to the overall structural response quantified by the force–displacement curve. These numerical studies demonstrate that the proposed framework effectively represents ductile fracture, yielding computational results that exhibit good agreement with experimental data.

36 MATERIALS SCIENCE↗

The influence of kinematics of deformation on polycrystalline halite dynamic recrystallization: Full-field simulation of simple shear versus pure shear

Rock salt, composed mainly of halite, flows viscoplastically over a wide range of geological conditions, strongly impacting the dynamic evolution of sedimentary basins and orogens. Understanding how dislocation creep, which involves dislocation glide, intracrystalline recovery and dynamic recrystallization, influences the microstructure and rheology of halite under various deformation kinematics and temperatures is crucial for enhancing knowledge of salt flow dynamics. Here, this study employs a full-field numerical simulation method to compare the viscoplastic deformation of polycrystalline halite under simple shear and pure shear conditions up to a natural strain of ε = 1.5 at temperatures ranging from 100 °C to 300 °C. The results are presented in terms of crystallographic preferred orientation (CPO), grain shape preferred orientation (SPO), subgrain boundary direction, grain size and strain rate distribution. The results indicate that the crystallographic anisotropy of individual halite crystals is transferred to the polycrystalline scale, resulting in strain localization, particularly in simple shear simulations. The kinematics of deformation affect the evolution and distribution of high strain-rate bands, determining the direction of intragranular substructures and the morphology of strain-induced grain boundaries, with minimal impact on grain size. The intensity of grain boundary migration increases with temperature, significantly influencing grain morphology and size, thereby obscuring strain localization, while having little effect on CPOs. At low strain (ε < 1.0), CPOs relative to both the maximum shortening direction and the grain SPO are similar regardless of the deformation kinematics. At high strain (ε > 1.0), simple shear CPOs exhibit three stronger {100} maxima with a monoclinic symmetry relative to the grain SPO compared to the six {100} maxima with an orthotropic symmetry relative to the grain SPO generated under pure shear. Therefore, microstructures and CPOs can serve as indicators of the strain path in polycrystalline halite under various conditions, aiding in determining the shear sense and elucidating the deformation kinematics of salt structures.

58 GEOSCIENCES↗

Surrogate Model Integration with MOOSE XFEM for Creep Crack Growth

Ferritic-martensitic steels are key structural materials for advanced reactors but experience time-dependent deformation and damage under prolonged high temperature and irradiation, leading to creep-driven crack initiation and growth. High-fidelity models—crystal plasticity with irradiation mechanisms, phase-field for microstructural evolution, and continuum-damage viscoplasticity—capture the underlying physics but are too computationally intensive for broad design-space exploration and uncertainty quantification. This milestone advances a scalable alternative by integrating a microstructure-sensitive surrogate creep model into the Multiphysics Object-Oriented Simulation Environment (MOOSE) finite element framework and extending it to fracture via the extended finite element method (XFEM). The surrogate model, developed with collaborators at Sandia and Los Alamos National Laboratories, maps relevant microstructural descriptors to the viscoplastic response of HT9. We embed this surrogate within a coupled deformation-damage workflow in MOOSE/XFEM to simulate creep-driven crack initiation and propagation. Implementation enhancements include updates to the material interface, a plastic correction phase involving microstructure evolution, and fracture criteria to ensure numerical robustness and compatibility with the surrogate structure. Demonstrations on canonical creep benchmarks spanning uniaxial and multiaxial states show that the surrogate reproduces key trends of high-fidelity models while substantially reducing computational cost. The resulting capability bridges physics fidelity and performance, providing a practical path to a predictive, microstructure-aware assessment of creep and fracture in reactor materials.

36 - MATERIALS SCIENCE↗

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↗

Bayesian calibration and uncertainty quantification of a rate-dependent cohesive zone model for polymer interfaces

In this work we present a rate-dependent cohesive zone model for the fracture of polymeric interfaces and performs a Bayesian calibration, an uncertainty quantification, and a sensitivity analysis for the model. The proposed cohesive zone model accounts for both reversible elastic and irreversible rate-dependent separation sliding deformation at the interface. The viscous dissipation due to the irreversible opening at the interface is modeled using elastic-viscoplastic kinematics that incorporates the effects of strain rate. Inverse calibration of parameters for such complex models through trial and error is challenging due to the large number of parameters of the model. Moreover, the calibrated parameter values are often non-unique and uncertain when the available experimental data is limited. To tackle this challenge, we employ a Bayesian calibration approach to identify parameters from experimental data, the resulting parameters significantly enhance the accuracy of the model. To quantify the uncertainty associated with the inverse parameter estimation, a modular Bayesian approach is employed to calibrate the unknown model parameters, accounting for the parameter uncertainty of the cohesive zone model. The advantages of the Bayesian calibration over a deterministic parameter fit are demonstrated. Further, to quantify the model uncertainties, such as incorrect assumptions or missing physics, a discrepancy function is introduced, which significantly improves the model’s prediction. Finally, the total uncertainty of the model is quantified in a predictive setting. A sensitivity analysis is performed to assess how changes in the input variables of the model affect the peak load, facilitating the identification of a concise set of highly influential parameters. The present approach can be used for calibration and uncertainty quantification for other complex computational mechanics models. It should also facilitate the designing of interface materials under uncertainty.

42 ENGINEERING↗

Assessment of thermal and radiation induced creep in the dual cooled lead lithium blanket

The creep in recently designed blanket for the fusion nuclear science facility is described by decomposing it into thermal, irradiation and cavity swelling creep. Initial estimates of thermal creep using elastic high temperature rules showed that the recent heat transfer improvements at the first wall (FW) is feasible, and the FW can withstand heat flux of 0.35 MW/m 2 at 3000 hrs without rupture. A viscoplastic model that is based on Norton's power law and relies on Multiphysics coupling of solid mechanics and heat transfer modules is used to capture the inelastic deformation in the blanket. Results showed that the design peak heat flux of 0.25 MW/m 2 produced maximum thermal creep strain of 0.45 % and the relaxation of the thermal stress at the FW. Irradiation creep is prescribed to be proportional to the displacement damage dose and applied stress. The displacements from irradiation creep radially decrease from the FW to the back wall, but the maximum deformation is found at the back wall that is connected to the Helium manifold due to the high stresses at the region. In conclusion, cavity swelling creep is dominant at steady state and combined with radiation creep to produce displacement of about 8 mm at the FW.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Role of grain boundary crystallography on void growth in FCC metals

The growth of an isolated spherical void at and near a grain boundary or a triple junction is studied by means of full-field crystal plasticity simulations with explicit representation of the void. To examine the fundamental aspects of void/grain boundary interactions in these cases, a three-dimensional large-strain elasto-viscoplastic fast-Fourier transform (LS-EVP-FFT) model with axisymmetric tensile periodic boundary conditions was developed. We employ this model to reveal the role of crystallography (grain orientations, grain boundary misorientation and inclination angle with respect to the largest component of the applied stress) and type of loading on void growth. It was found that grain orientation was the primary crystallographic feature governing void growth in strain-rate controlled simulations. For strains small enough to not significantly alter the local geometry, the overall growth of the void in bicrystal boundaries was not strongly affected by grain boundary inclination angle. The relative amount of void growth experienced within each grain in a bicrystal or tri-crystal was strongly dependent on the grain boundary inclination angle. We found a subset of orientations that exhibit faster void growth in bicrystal boundaries than inside the constituent single crystals. Last, we show that when placed primarily in one grain but near the grain boundary, void growth was similar to void growth when placed at the grain boundary. Furthermore, the growing void remained in the original grain and did not grow into the other grain, regardless of crystallographic character.

36 MATERIALS SCIENCE↗

Effect of surface defects on damage and failure of dynamically driven cylindrical shells

Here, in this work, expanding cylinder experiments are performed using a gas gun on the aluminum 6061 alloy to understand the dynamic deformation and failure behavior of this material in the presence of surface defects in the form of scratches. The experiments are then modeled with the modified Tepla (for tensile plasticity) model, which is equipped with the damage yield function of Gurson [“Continuum theory of ductile rupture by void nucleation and growth: Part I—Yield criteria and flow rules for porous ductile media,” (1977)] and Tvergaard and Needleman [Acta Metall. 32(1), 157–169 (1984)], over-stress viscosity [F. Addessio and J. Johnson, J. Appl. Phys. 74(3), 1640–1648 (1993)], and the underlying viscoplastic model of Preston et al. [J. Appl. Phys. 93(1), 211–220 (2003)] (Preston–Tonks–Wallace model) for the dense material. The Tepla model parameters were calibrated using the Bayesian approach against experimentally measured free-surface velocity data obtained from plate impact experiments on the aluminum 6061 alloy. The evolution of radial expansion and shapes of the deformed cylinder from the simulations on a pristine cylinder wall and a cylinder wall with a longitudinal surface defect shows excellent agreement with the corresponding experimental results. Finally, the role of surface defect size and cylinder wall geometry on the dynamic strength of the cylinder wall was also investigated through a series of simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗