SEARCH · Search NASA
Results for “rate-dependent plasticity”
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.
Bounds on the Rate-Dependent Plastic Flow of Tantalum up to 75 GPa
In this study we report on the continued development of thermodynamics-based analysis of shock waves propagation with the objective of extracting information related to materials strength at high strain rates and pressures. Building on previous results reported for peak stresses of 10 GPa and 25 GPa, we present a series of three-step gas-gun shock experiments designed to explore the pressure and strain rate dependence of plastic flow in polycrystalline tantalum. These experiments at nominal peak stresses of 50 GPa and 75 GPa show the irreversible deformation before pullback to be almost entirely confined to the shock loading, with negligible plastic relaxation on the post-shock plateau. We also add a reverse-ballistics shot at 25 GPa, which was designed to reveal the pullback response with negligible interference from free-surface effects. General thermodynamic considerations allow us to place bounds on the plastic behavior even for parts of the curve that change far too rapidly for the velocimetric time resolution of (conservatively) ~ 5 ns. To analyze the data, we found it necessary to substantially improve the interpolation/extrapolation scheme in order to improve its robustness, flexibility and range of applicability. We describe the new scheme based on splines, as well as an extension of free-surface corrections to the post-shock rarefaction waves. Reanalysis with the new scheme produces results essentially within the error bars previously reported, showing that the known systematic errors associated with free-surface effects are relatively inconsequential for determining thermodynamic paths.
Convergence Analysis of Rate-Dependent J2 Plasticity
End of summer internship write-up/memo
A high-order, localized-artificial-diffusivity method for Eulerian simulation of multi-material elastic-plastic deformation with strain hardening
A high-order method for Eulerian simulation of material undergoing large elastic–plastic deformation is developed. Thermodynamically consistent hyperelastic constitutive relations are assumed, facilitating the treatment of solids, liquids, and gases in a unified manner. Here, the method enables the simulation of multi-material interactions using a diffuse interface approach. Numerical capturing of material interfaces, shock waves, contact surfaces, and elastic-plastic strain discontinuities using high-order compact-difference schemes is assisted by Localized Artificial Diffusivity (LAD). In the new setting involving elastic–plastic deformation, the previously established terms for the artificial properties are verified to effectively regularize normal shocks. Additional LAD terms are introduced to the elastic and plastic kinematic equations to regularize shear shocks and other strain discontinuities, improving solution stability. Other important features of the method that improve robustness include the numerical treatment of compatibility terms in the kinematic equations, and the treatment of rotation. Particular emphasis is focused toward new advancements of the methods for plastic-deformation integration and the associated strain hardening of the material, including rate-dependent plasticity. The method is demonstrated on a variety of test problems, including 1-D impacts, a variant of the Shu-Osher problem, a Taylor impact, and a Richtmyer-Meshkov instability between two elastic–plastic solids with strain hardening.
The physics of desiccation cracks 1: Ductile fracturing and dependence on relative humidity
Ductile deformation is ubiquitously found in the shrinkage of geomaterials. The existence of ductility requires elastoplastic mechanics when analyzing the structure deformation under external stress. Here, in this work, we explore the physics of ductile fracturing based on the results from desiccation experiments and triaxial tests. By using the digital image correlation (DIC) method to generate the strain maps of samples undergoing desiccation cracking under different relative humidities, we obtain results showing the previously postulated Cnoidal Wave ductile failure patterns propagating under atmospheric condition-controlled crack velocities. We then correlate these observations with rate-dependent plasticity models calibrated through triaxial tests undergoing several unloading–reloading cycles and velocity stepping. This work demonstrates the necessity to consider ductility in soil cracking, indicating that the formation of crack patterns in soil desiccation is a slow and predictable process.
A macro-micro approach for identifying crystal plasticity parameters for necking and failure in nickel-based alloy haynes 282
Here, this work develops a two-scales macro-micro approach to address the challenge in calibrating crystal plasticity microstructural models when samples undergo necking prior to fracture. The crystal plasticity models are crucial for predicting the materials’ plastic deformation and failure at the microstructure level, identifying the materials’ intrinsic properties as well as investigating the microstructure-properties relationships. However, after necking occurs, the experimentally measured stress-strain curves fail to reflect the materials ‘true’ stress-strain behavior and cannot be directly fitted into crystal plasticity models. The proposed macro-micro approach employs a top-down strategy to address this challenge, which has been studied with experimental tests on precipitation-strengthened Ni-based superalloy Haynes® 282®. In this approach, a macro rate-dependent anisotropic plasticity model with Voce-type hardening and Rice-Tracey damage law is first utilized to model the deformation and failure of the tensile bar, and calibrated by matching the stress-strain curves, necking strain, and reduction of area. Especially, to match the testing results under different applied strain rates, the rate-sensitivity parameter m and saturation stress in the elasticity model are modified to incorporate dependence on the local strain rate. Then, the ‘true’ stress-strain behaviors are extracted from the necking zone of the macro-model, which are used to calibrate a micro-model with explicit microstructures and governed by an extended crystal plasticity law. The consistency between the micro-model and macro-model are enforced during calibration. The calibration outcomes from the crystal plasticity model elucidate the materials intrinsic properties for slip, hardening, and failure, which is vital for further investigations into the microstructure-properties relationship and for accurate prediction of the material behavior under various test and service conditions.
Strongly nonlinear wave propagation in elasto-plastic metamaterials: Low-order dynamic modeling
Nonlinear elastic metamaterials are known to support a variety of dynamic phenomena that enhance our capacity to manipulate elastic waves. Since these properties stem from complex, subwavelength geometry, full-scale dynamic simulations are often prohibitively expensive at scales of interest. Prior studies have therefore utilized low-order effective medium models, such as discrete mass-spring lattices, to capture essential properties in the long-wavelength limit. While models of this type have been successfully implemented for a wide variety of nonlinear elastic systems, they have predominantly considered dynamics depending only on the instantaneous kinematics of the lattice, neglecting history-dependent effects, such as wear and plasticity. Here, to address this limitation, the present study develops a lattice-based modeling framework for nonlinear elastic metamaterials undergoing plastic deformation. Due to the history- and rate-dependent nature of plasticity, the framework generally yields a system of differential-algebraic equations whose computational cost is significantly greater than an elastic system of comparable size. We demonstrate the method using several models inspired by classical lattice dynamics and continuum plasticity theory and explore means to obtain empirical plasticity models for general geometries, thereby gaining insight into the influence of microstructural plasticity on effective material performance, which can be used to improve the design of nonlinear mechanical metamaterials.
Insights into mixing of non-isothermal multi-polymer melts for complex plastics recycling
Catalytic recycling or upcycling of plastics is often limited not by catalyst performance, but by transport, arising from highly viscous, non-Newtonian polymer melts. In this work, we develop a reactor-scale framework that integrates rheological measurements, constitutive modeling, computational fluid dynamics (CFD), and experiments to quantify mixing, heat transfer, and dispersion in surrogate hydrocarbon melts representing mixed plastics systems. Temperature- and shear rate-dependent viscosity of low-density polyethylene (LDPE) and high-density polyethylene (HDPE) is measured to create two surrogate polymers (PLD and PHD) that capture the dominant shear-thinning flow behavior while neglecting strong elastic effects, enabling tractable simulation of non-isothermal, polymer-melt mixing using a Carreau-Arrhenius generalized Newtonian framework. Three-dimensional CFD simulations are employed to evaluate impeller performance in PLD using mixing time, cavern volume, thermal uniformity, and interfacial area for regimes in which viscoelastic effects are not dominant. We show that magnetic stir bars commonly used in lab-scale studies produce large thermal gradients (~60 °C) and poor mixing, even under idealized power delivery and polymer flow conditions. In contrast, close-clearance anchor impellers achieve near-isothermal operation, reduce mixing times by up to 5×, and provide >90% active circulation volume. We further demonstrate that, at low pseudo-Deborah number (De*), motor power requirements can be predicted directly from shear rate-dependent rheology using the Carreau-Arrhenius framework, enabling rational selection of operating conditions. Extension to surrogate immiscible multi-polymer systems based on PLD and PHD shows that interfacial area is highly sensitive to operating conditions and impeller design, with coaxial anchor-turbine configurations enhancing dispersion by up to 4 × .
Modeling plasticity-mediated void growth at the single crystal scale: A physics-informed machine learning approach
Modeling the evolution of voids during plastic flow as well as their effects on plastic dissipation is critical for both component manufacturing and lifetime estimation purposes. To this end, we propose a rate-dependent constitutive model to homogenize the effects of semi-randomly distributed voids on single crystal plasticity whilst capturing void interaction and plastic anisotropy. Here, this present work focuses on the case of face centered cubic crystals to introduce an anisotropic gauge function applicable within the crystal plasticity formalism. The approach combines analytical methods to describe the micromechanics of the system in combination with symbolic regression to capture analytically intractable mechanisms from data. The hybrid framework uses a physics-informed genetic programming-based symbolic regression algorithm to solve a multiform optimization problem simultaneously producing a new gauge function and a new strain rate equation. This is also a multi-objective optimization problem with many competing objectives. A new search and selection step is introduced to the genetic algorithm that promotes convergence toward a global solution that better satisfies all the objectives. Overall, the symbolic equations produced leverage data-driven methods to achieve greater accuracy than comparable alternatives on an analytically intractable problem while maintaining model transparency.
A rheological model for loose sands with insights from DEM
A rheological model for loose granular media is developed to capture both solid-like and fluid-like responses during shearing. The proposed model is built by following the mathematical structure of an extended Kelvin–Voigt model, where an elastic spring and plastic slider act in parallel to a viscous damper. This arrangement requires the partition of the total stress into rate-independent and rate-dependent stress components. To model the solid-like behavior, a simple frictional plasticity model is adopted without modifications, thus contributing to the rate-independent stress. Instead, the fluid-like or rate-dependent stress is further decomposed into deviatoric and volumetric parts, by proposing a new formulation based on a combination of the μ(I) relation, originally developed under pressure-controlled shear, with a pressure-shear rate relation derived under volume-controlled shear. The proposed formulation allows the model to capture both the increase in the friction coefficient and the enhanced dilation at high shear rates. High-fidelity simulation data, obtained from discrete element method and multiscale modelling, are used to evaluate the performance of the proposed constitutive model. The model provides accurate results under both drained and undrained simple shear paths across a wide range of shear rates. Furthermore, it successfully reproduces at much lower computational cost the flowslide mobility computed through multiscale simulations, which is primarily regulated by the shear rate dependence of the material properties during the dynamic runout stage.
Microstructure characterization and elastic-plastic self-consistent simulation studies of anisotropic deformation of β -tin
Uniaxial compression tests of 99.9% pure polycrystalline $β$-tin (Sn) were conducted at various strain rates (10 -3 /s, 10 -1 /s and 1/s) and temperatures (294 K and 193 K) to understand the effects of these variables on the stress–strain response, microstructure evolution, and stress relaxation behavior. Multiple Sn specimens were subjected to complex compressive loading/unloading/reloading paths at different strain rates. Specimens initially compressed at higher strain rates showed strain rate-dependent texture evolution and more pronounced relaxation upon reloading, as compared to those pre-strained at lower rates at room temperature (RT). Compression tests conducted at low temperature (193 K) revealed increased strength, similar to the enhanced strength observed when the strain rate was increased at RT. Neutron diffraction was employed to characterize the initial and final bulk textures of the RT specimens. Electron backscatter diffraction was utilized to examine the crystallographic grain orientation and morphology, thereby identifying the signatures of dislocation-mediated deformation, grain refinement, recrystallization, and twinning behaviors. Elastic-plastic self-consistent simulations were performed to investigate the deformation modes responsible for the strain rate-dependent macroscopic stress–strain response and texture evolution. Effects of crystallographic orientation on stress relaxation behavior was also examined. In conclusion, the model predictions are in reasonable agreement with experimental observations.
A crystal plasticity finite element model embedding strain-rate sensitivities inherent to deformation mechanisms: Application to alloy AZ31
We report the fundamental power-law relationship representing the flow rule in crystal visco-plasticity ensures uniqueness in the selection of slip systems accommodating imposed plastic strain-rates. The power-law relationship also introduces an artificially high strain-rate sensitivity in crystal plasticity simulations, unless a high value of the power-law exponent is used. However, the use of high values for the exponent is limited by numerical tractability. This paper presents a numerical method implemented in a crystal plasticity finite element (CPFE) model for embedding any value of the power-law exponent reflecting the true material strain-rate sensitivity. Importantly, the method does not increase computation time involved in the simulations. The enhanced CPFE model is used to interpret and predict a complex strain-rate sensitive response and microstructural evolution of AZ31 Mg alloy. Measured values of strain-rate sensitivity for slip and twinning modes are used in the simulations. Calculations show that the model successfully captures the phenomena pertaining to the effect of changing applied strain-rate on the mechanical response including flow stress and evolution of texture and twinning for a broad range of strain-rates ranging from 10 -3 s -1 to 10 3 s -1 and loading orientations in tension and compression. It is shown that such predictions are a consequence of not only relative amounts of slip and twinning activities driven by a set of accurately characterized hardening law parameters but also values of the strain-rate sensitivities inherent to individual deformation mechanisms. Besides, the model validates the measured strain-rate dependency of deformation mechanisms while accurately reproducing the mechanical data. Hence, the model can be used to verify and further refine or infer measured strain-rate sensitivity per deformation mechanism by reproducing experimental data.
Report on FY 2021 testing in support of integrating Alloy 800H and 2.25Cr-1Mo for EPP analysis
The use of simplified Elastic-Perfectly Plastic (EPP) design methods avoids the complexities and limitations of the elastic analysis approach in evaluating strain limits and creep-fatigue damage for Class A metallic coolant boundary components in ASME Boiler and Pressure Vessel Code, Section III, Division 5, Subsection HB, Subpart B. No stress classification and linearization are required in the EPP methods. Further, the strain limits and creep-fatigue damage evaluation procedures based on elastic analysis results were developed for the reactor operating conditions where uncoupled, rate-independent plasticity and stationary creep can be used to approximate the deformation in the component adequately. But at very high temperatures, the deformation is rate-dependent, and it can only be approximated adequately by unified viscoplastic model. Thus, for example, the elastic analysis approach cannot be used to evaluate strain limits and creep-fatigue damage for Alloy 617 above 650 °C. The EPP methods remove the temperature restriction. They can be used to evaluate strain limits and creep-fatigue damage for the full range of temperatures.
2024 Second Half Semi Annual Report: Modeling plasticity-mediated flow in metals with pressurized cavities
The objective is to better predict the bulk-scale mechanical behavior of porous metals that have over pressurized cavities (e.g., irradiated metals with helium bubbles) by quantifying the complex coupling among cavity aspects (e.g., size distribution, inhomogeneous overpressure values, spatial arrangement) and metal properties (e.g., rate-dependency, crystallographic lattice). This requires up-scaling local mechanical fields from the single crystal scale and will be accomplished using a homogenization approach that combines full-field numerical simulations, analytical formalisms, and physics-informed machine learning to produce symbolically-defined constitutive equations (e.g., gauge functions). These equations will satisfy the objective because they enable computationally efficient predictions that approach the accuracy of computationally expensive full-field numerical simulations, abide by theoretical requirements (e.g., conservation of energy, work conjugacy), and retain the transparency of analytical models.
Strain-rate dependent deformation mechanisms in single-layered Cu, Mo, and multilayer Cu/Mo thin films
Here, strain-rate sensitivity and rate-dependent hardness, over a range of 10 -2 to 10 2 s -1 , of sputter-deposited single-layered Cu, Mo, and 5 nm Cu/ 5 nm Mo, and 100 nm Cu/ 100 nm Mo multilayer films with a total film thickness of 5 μm were measured using nanoindentation. The plastic zone underneath the nanoindents was characterized via cross-sectional transmission electron microscopy (XTEM). The multilayer films exhibited enhanced hardness but slightly reduced strain-rate sensitivity with decreasing layer thickness from 100 nm to 5 nm. Only the 5 nm Cu/ 5 nm Mo multilayer film exhibited shear bands underneath the nanoindents, and the size of the shear bands increased with increasing strain rate. In contrast, the 100 nm Cu/ 100 nm Mo multilayer film exhibited material pile-up around the indents and significant nanotwinning within Cu grains. The effect of strain rate and layer thickness on the hardness and strain rate sensitivity of the multilayer thin films is interpreted using a modified confined layer slip (CLS) model. The reduced rate sensitivity at 5 nm as compared to 100 nm correlates with abundant growth nanotwins in the Cu grains in 100 nm and formation of shear bands in 5 nm multilayers. In single layer films, a substructure with a high density of dislocations was observed consistent with the plastic strain gradient in the indent plastic zone. No evidence of deformation twins was noted in any of the samples.
Simulation of dislocation evolution in microparticle impacts over a wide range of impact velocities
We report the dynamic loading environment created by high velocity microparticle impacts is very difficult to model. The extremely high strain rates, large deformations, and high temperatures affect the flow stress of the material in ways that standard flow stress models cannot capture. Plastic deformation and strain hardening in metals is controlled by the motion of dislocations. Dislocations can be nucleated, stored as forest dislocations, or be annihilated as loading progresses. A comprehensive accounting of dislocation density change is needed to accurately describe rate and temperature dependent dislocation glide and evolution across the wide rage of loading conditions present in microparticle impact problems. Therefore, we implement and apply a newly proposed flow stress model (Hunter and Preston, 2015, 2022) that is then coupled with mobile and immobile dislocation density evolution equations. This model is implemented in Los Alamos National Laboratory’s hydrodynamics code, FLAG, to model copper-on-copper microparticle impacts. This new strength model allows for accurate simulation of particle rebound and flattening across a wide range of impact velocities.
Atomistic deformation mechanism of silicon under laser-driven shock compression
Abstract Silicon (Si) is one of the most abundant elements on Earth, and it is the most widely used semiconductor. Despite extensive study, some properties of Si, such as its behaviour under dynamic compression, remain elusive. A detailed understanding of Si deformation is crucial for various fields, ranging from planetary science to materials design. Simulations suggest that in Si the shear stress generated during shock compression is released via a high-pressure phase transition, challenging the classical picture of relaxation via defect-mediated plasticity. However, direct evidence supporting either deformation mechanism remains elusive. Here, we use sub-picosecond, highly-monochromatic x-ray diffraction to study (100)-oriented single-crystal Si under laser-driven shock compression. We provide the first unambiguous, time-resolved picture of Si deformation at ultra-high strain rates, demonstrating the predicted shear release via phase transition. Our results resolve the longstanding controversy on silicon deformation and provide direct proof of strain rate-dependent deformation mechanisms in a non-metallic system.
Digital image correlation and infrared thermography data for seven unique geometries of 304L stainless steel
Material Testing 2.0 (MT2.0) is a paradigm that advocates for the use of rich, full-field data, such as from digital image correlation and infrared thermography, for material identification. By employing heterogeneous, multi-axial data in conjunction with sophisticated inverse calibration techniques such as finite element model updating and the virtual fields method, MT2.0 aims to reduce the number of specimens needed for material identification and to increase confidence in the calibration results. To support continued development, improvement, and validation of such inverse methods—specifically for rate-dependent, temperature-dependent, and anisotropic metal plasticity models—we provide here a thorough experimental data set for 304L stainless steel sheet metal. The data set includes full-field displacement, strain, and temperature data for seven unique specimen geometries tested at different strain rates and in different material orientations. Commensurate extensometer strain data from tensile dog bones is provided as well for comparison. We believe this complete data set will be a valuable contribution to the experimental and computational mechanics communities, supporting continued advances in material identification methods.