Search NASA⌕ Search

SEARCH · Search NASA

Results for “crystal plasticity simulation”

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

Anisotropic strength behavior of single-crystal TATB

High-rate strength behavior plays an important role in the shock initiation of high explosives, with plastic deformation serving to localize heat into hot spots and as a mechanochemical means to enhance reactivity. Recent simulations predict that detonation-like shocks produce highly reactive nanoscale shear bands in the layered crystalline explosive TATB (1,3,5-triamino-2,4,6-trinitrobenzene), but the thresholds leading to this response are poorly understood. We utilize molecular dynamics to simulate the high-rate compressive stress–strain response of TATB, with a focus on understanding flow behavior. Here, the dependence of strength on pressure and loading axis (crystal orientation) is explored. The deformation mechanisms fall broadly into two categories, with compression along crystal layers activating a buckling/twinning mode and compression normal to the layers producing nanoscale shear bands. Despite the complexity of the underlying mechanisms, the crystal exhibits relatively straightforward stress–strain curves. Most of the crystal orientations studied show rapid strain softening following the onset of yielding, which settles to a steady flow state. Trajectories are analyzed using five metrics for local states and structural order, but most of these metrics yield similar distributions for these deformation mechanisms. On the other hand, a recently proposed measure of intramolecular strain energy is found to most cleanly distinguish between these mechanisms, while also providing a plausible connection with mechanochemically accelerated decomposition kinetics. Localization of intramolecular strain energy is found to depend strongly on crystal orientation and pressure.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Cross-slip and easy-glide CRSS in titanium: Theoretical predictions and in-situ TEM measurements

This study investigates the mechanics of prismatic and first-order pyramidal $\langle$a$\rangle$ slip in titanium (Ti), elucidating the physics of easy-glide and cross-slip through a combination of theory and experiments. Screw-character prismatic (Pr) dislocations in Ti are of particular interest because their complex cores can be stable or unstable, leading to activation by either cross-slip or planar glide. To investigate these mechanisms, site-specific micro tensile samples were prepared using focused ion beam (FIB) milling and mounted on a push-to-pull (PTP) device for in-situ transmission electron microscopy (TEM) tensile testing. The in-situ experiments provide direct observations of the onset of dislocation motion and the precise determination of the critical resolved shear stress (CRSS) for the activated mechanisms, and its evolution with load cycling. A comprehensive theory has been developed to predict the CRSS values for easy glide, cross-slip, and multiplication of dislocations. Predicted critical stresses for pyramidal (π)-to-Pr and reverse cross-slip agree closely with the experimental measurements. The latter cross-slip stress is a factor of two higher than that of unobstructed planar slip. The model accounts for overlapping dislocation cores and employs a Wigner-Seitz based cell to evaluate misfit energies. By combining ab initio density functional theory (DFT) with anisotropic elasticity, the framework identifies minimum energy pathways for dislocation glide, which can be intermittent and zig-zag. A simplified expression utilizing (π) and Pr Schmid factor ratios is proposed for critical stress corresponding to (π)-to-Pr cross-slip transition. The results are strongly dependent on crystal orientation, underscoring non-Schmid behavior. Overall, this study explores key critical stress parameters essential for informing higher-scale simulations of plasticity in Ti.

Critical stress↗
DOE OSTI2024

Designing Ti-6Al-4V microstructure for strain delocalization using neural networks

Abstract The deformation behavior of Ti-6Al-4V titanium alloy is significantly influenced by slip localized within crystallographic slip bands. Experimental observations reveal that intense slip bands in Ti-6Al-4V form at strains well below the macroscopic yield strain and may serially propagate across grain boundaries, resulting in long-range localization that percolates through the microstructure. These connected, localized slip bands serve as potential sites for crack initiation. Although slip localization in Ti-6Al-4V is known to be influenced by various factors, an investigation of optimal microstructures that limit localization remains lacking. In this work, we develop a novel strategy that integrates an explicit slip band crystal plasticity technique, graph networks, and neural network models to identify Ti-6Al-4V microstructures that reduce the propensity for strain localization. Simulations are conducted on a dataset of 3D polycrystals, each represented as a graph to account for grain neighborhood and connectivity. The results are then used to train neural network surrogate models that accurately predict localization-based properties of a polycrystal, given its microstructure. These properties include the ratio of slip accumulated in the band to that in the matrix, fraction of total applied strain accommodated by slip bands, and spatial connectivity of slip bands throughout the microstructure. The initial dataset is enriched by synthetic data generated by the surrogate models, and a grid search optimization is subsequently performed to find optimal microstructures. Describing a 3D polycrystal with only a few features and a combination of graph and neural network models offer robustness compared to the alternative approaches without compromising accuracy. We show that while each material property is optimized through a unique microstructure solution, elongated grain shape emerges as a recurring feature among all optimal microstructures. This finding suggests that designing microstructures with elongated grains could potentially mitigate strain localization without compromising strength.

Ahmadikia, Behnam↗
DOE OSTI2025

Fatigue life prediction of powder bed fused–laser beam AlSi10Mg: Incorporating critical defects via crystal plasticity modelling

The current study provides a microstructurally-based computational framework to predict the fatigue life of additive manufactured (AM), i.e., powder bed fused–laser beam (PBF-LB), AlSi10Mg specimens using the crystal plasticity finite element method (CPFEM). The fractography analysis, electron backscatter diffraction (EBSD), uniaxial and cyclic responses, and fatigue life of specimens were used to inform the computational framework. CPFE simulation was used to compute fatigue indicator parameters (FIPs) as fatigue driving forces. A new fatigue criterion is introduced based on FIPs, which was calibrated using experimental fatigue data. The proposed fatigue measure was evaluated versus the specimens with critical defects of various sizes and locations subjected to different stress amplitudes. The results show that the developed framework can capture the fatigue life of samples with different critical defect locations and sizes along with different stress amplitudes for both high-cycle fatigue (HCF) and very high-cycle fatigue (VHCF) regimes.

Additive manufacturing↗
DOE OSTI2023

Predicting Failure Using Deep Learning SAND Report

Accurate prediction of ductile failure is critical to Sandia’s NW mission, but the models are computationally heavy. The costs of including high-fidelity physics and mechanics that are germane to the failure mechanisms are often too burdensome for analysts either because of the person-hours it requires to input them or because of the additional computational time, or both. In an effort to deliver analysts a tool for representing these phenomena with minimal impact to their existing workflow, our project sought to develop modern data-driven methods that would add microstructural information to business-as-usual calculations and expedite failure predictions. The goal is a tool that receives as input a structural model with stress and strain fields, as well as a machine-learned model, and output predictions of structural response in time, including failure. As such, our project spent substantial time performing high-fidelity, three-dimensional experiments to elucidate materials mechanisms of void nucleation and evolution. We developed crystal-plasticity finite-element models from the experimental observations to enrich the findings with fields not readily measured. We developed engineering length-scale simulations of replicated test specimens to understand how the engineering fields evolve in the presence of fine-scale defects. Finally, we developed deep learning convolutional neural networks, and graph-based neural networks to encode the findings of the experiments and simulations and make forward predictions in time for structural performance. This project demonstrated the power of data-driven methods for model development, which have the potential to vastly increase both the accuracy and speed of failure predictions. These benefits and the methods necessary to develop them are highlighted in this report. However, many challenges remain to implementing these in real applications, and these are discussed along with potential methods for overcoming them.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Tensile deformation behavior of twist grain boundaries in CoCrFeMnNi high entropy alloy bicrystals

High entropy alloys (HEA) are a class of materials that consist of multiple elemental species in similar concentrations. The use of elements in far from dilute concentrations introduces a multi-dimensional composition design space by which the properties of metallic systems can be tailored. While the mechanical behavior of HEAs has been the subject of active research recently, the role of grain boundaries (GBs) in their deformation behavior remains poorly understood. Motivated by recent experiments on HEAs demonstrating that GBs act as nucleation sites for deformation twins, herein, we leverage atomistic simulations to construct a series of equiatomic CoCrFeMnNi HEA bicrystals with $\langle 110 \rangle$ and $\langle 111 \rangle$ symmetric twist GBs and examine their tensile behavior and underlying deformation mechanisms at 77 K. Simulation results reveal that plastic deformation proceeds by the nucleation of partial dislocations from GBs, which then grow with further loading by bowing into the bulk crystals leaving behind stacking faults. Variations in the nucleation stress exist as function of GB character, defined in this work by the twist angle. Our results provide future avenues to explore GBs as a microstructure design tool to develop HEAs with tailored properties.

36 MATERIALS SCIENCE↗
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 OSTI2022

Effect of loading path on grain misorientation and geometrically necessary dislocation density in polycrystalline aluminum under reciprocating shear

Solid phase processing (SPP) is a promising alloy fabrication technique to produce fine and homogeneous grain structures for high-performance alloys. However, there is very limited modeling capability to understand and predict the grain refinement during SPP. In this work, the crystal plasticity theory was used to study elastic-plastic deformation in polycrystalline aluminums under large shear deformation. Two approaches, kernel averaged misorientation (KAM) and grain reference orientation deviation (GROD), were used to assess the grain misorientations. The geometrically necessary dislocation (GND) density was computed with the plastic strain rate. The deformation simulations were carried out under two loading conditions to investigate the effect of loading paths on the evolutions of grain misorientation and GND density. The results show that the regions with high misorientation and GND density first appear near grain boundaries. These regions then extend toward interior grains. The loading path affects dislocation system activation and dislocation recovery, hence dislocation evolution and misorientation. In conclusion, both two- and three-dimensional simulations showed that the spatial and temporal evolutions of GROD, KAM, and GND density in are closely correlated, which indicates they all can be used as criteria of grain refinement or recrystallization.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Eulerian finite element simulations of the drop weight impact test with a dislocation Density-based continuum model

During a drop weight impact test, the kinetic energy of the falling weight is transferred into the sample resting on the anvil. The plastic deformation in the sample is an important mechanism for the dissipation of the input kinetic energy. We use Eulerian finite element analysis to simulate the deformation and temperature evolution in a 1,3,5-trinitro-1,3,5-triazine (RDX) sample consisting of multiple crystals. In Eulerian finite element simulations, the mesh moves relative to the material. After every change of position between the mesh and the material, the state variables are interpolated to the new mesh position, i.e., advection. In an effort to reduce the advection errors, we use a rate form of a dislocation density-based continuum model by Luscher et al. Here, the simulations predict localization of plastic deformation, and plastic dissipation as a significant source of heat generation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Effect of neighboring grain orientation on strain localization in slip bands in HCP materials

Particularly in plastically anisotropic crystals, such as hexagonal close packed (HCP) materials, plastic deformation is realized by slip acting in the small volumes within individual crystals. Here we extend a full field fast Fourier transform (FFT)-based elasto-viscoplastic formulation to simulate the development of a single slip band on either prismatic or basal planes spanning a crystal. Calculations of the strain and stress fields induced locally within the band and parent crystal, and ahead of the band/grain boundary junction in the neighboring crystal are analyzed as the slip band intensifies under increasing applied strain. We report a substantial influence of the crystallographic orientation of the nearest neighboring grain on the rate of slip band localization. Performing the analysis on two materials, CP-Ti and Mg, indicates that the strength of the material affects the rate of localization, with stronger materials tending to localize more easily. A slip band tip stress-based criterion is proposed for identifying the nearest neighbor orientations in which slip band transmission is possible and the likely slip system for which it occurs. This indicator is validated against experimental studies on commercially pure Ti, an Mg–Y alloy, and Ti–6Al–4V. Finally, we show that for low GB misorientations, the slip band is likely to transmit into another slip band of the same type in the neighbor grain, while for high GB misorientations, it is likely to transmit into one of a different type or to not transmit at all.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Elastic Properties and Deformation Mechanisms in the van der Waals Single-Crystalline Indium Selenide

The bulk van der Waals (vdW) single-crystalline indium selenide (β-InSe) displays exceptional plasticity in a layered crystalline form at both micro- and macroscale. However, the nanoscale origin of plasticity remains unclear. Herein, an atomic-level study on the deformation mechanisms of InSe by using first-principles calculations is reported. Remarkable anisotropic elasticity is observed in the vdW InSe layered crystal, and the stiffness is dramatically softened because of the vdW gap in the layered structure. The simulations capture the distinct fracture modes in the uniaxial tensile deformation, depending on the loading directions—brittle fracture in the [100] and [110] directions while ductile failure in the [001] direction. The InSe layered crystal structure exhibits superplastic deformability under uniaxial compression. Different transition pathways, including interlayer tangling, amorphization, and cross-linking, are tracked along respective deformation directions. The unprecedented plasticity of InSe layered crystals can be attributed to the phase transition coupled with interlayer gliding and cross-layer dislocation slipping. Here, this study deepens our understanding of the deformation mechanisms of layered materials at the atomistic level and provides insights into tailoring material properties for low-dimensional material design based on its deformation mechanisms.

2D materials↗
DOE OSTI2025

Calibration of RAFM Micromechanical Model for Creep Using Bayesian Optimization for Functional Output

A Bayesian optimization procedure is presented for calibrating a multimechanism micromechanical model for creep to experimental data of F82H steel. Reduced activation ferritic martensitic (RAFM) steels based on Fe(8–9)%Cr are the most promising candidates for some fusion reactor structures. Although there are indications that RAFM steel could be viable for fusion applications at temperatures up to 600°C, the maximum operating temperature will be determined by the creep properties of the structural material and the breeder material compatibility with the structural material. Due to the relative paucity of available creep data on F82H steel compared to other alloys such as Grade 91 steel, micromechanical models are sought for simulating creep based on relevant deformation mechanisms. As a point of departure, this work recalibrates a model form that was previously proposed for Grade 91 steel to match creep curves for F82H steel. Due to the large number of parameters (9) and cost of the nonlinear simulations, an automated approach for tuning the parameters is pursued using a recently developed Bayesian optimization for functional output (BOFO) framework (Huang et al., 2021, “Bayesian optimization of functional output in inverse problems,” Optim. Eng., 22, pp. 2553–2574). Incorporating extensions such as batch sequencing and weighted experimental load cases into BOFO, a reasonably small error between experimental and simulated creep curves at two load levels is achieved in a reasonable number of iterations. In conclusion, validation with an additional creep curve provides confidence in the fitted parameters obtained from the automated calibration procedure to describe the creep behavior of F82H steel.

42 ENGINEERING↗
DOE OSTI2020

Particle interspacing effects on the mechanical behavior of a Fe–TiB 2 metal matrix composite using FFT-based mesoscopic field dislocation mechanics

This paper presents an application to metal matrix composites (MMCs) of an enhanced elasto-viscoplastic Fast Fourier Transform (EVP-FFT) formulation coupled with a phenomenological continuum Mesoscale Field Dislocation Mechanics (MFDM) theory. Contrary to conventional crystal plasticity, which only accounts for plastic flow and hardening induced by statistically stored dislocations (SSDs), MFDM-EVP-FFT also describes the evolution of polarized geometrically necessary dislocation (GND) density and its effect on both plastic flow and hardening. Numerical results for a Fe–TiB 2 MMC made of a ferrite matrix (α-Fe) and elastic ceramic particles (TiB 2 ) are presented. Full-field simulations are performed using synthetic periodic unit cells representative of the MMC, with single-crystalline and polycrystalline matrix, for different particle interspacing distances. A strong dependence of the predicted equivalent stress, cumulated plastic strain and GND density fields with particle interspacing distance is observed, in contrast with conventional crystal plasticity. Correlations between these mechanical fields and microstructural features, and their influence on local and global mechanical behavior are examined for the different MMC microstructures.

42 ENGINEERING↗
DOE OSTI2025

Anomalous entropy-driven kinetics of dislocation nucleation

The kinetics of dislocation reactions, such as dislocation multiplication, controls the plastic deformation in crystals beyond their elastic limit, therefore critical mechanisms in a number of applications in materials science. We present a series of large-scale molecular dynamics simulations that shows that one such type of reactions, the nucleation of dislocation at free surfaces, exhibit unconventional kinetics, including unexpectedly large nucleation rates under compression, very strong entropic stabilization under tension, as well as strong non-Arrhenius behavior. These unusual kinetics are quantitatively rationalized using a variational transition state theory approach coupled with an efficient numerical scheme for the estimation of vibrational entropy changes. These results highlight the need for a variational treatment of the kinetics to quantitatively capture dislocation reaction kinetics, especially at low-to-moderate strains where large deformations are required to activate reactions. These observations suggest possible explanations to previously observed unconventional deformation kinetics in both molecular dynamics simulations and experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Analytic model of the remobilization of pinned glide dislocations: Including dislocation drag from phonon wind

Here, we discuss the effect of a non-constant dislocation drag coefficient on the very high strain rate regime within an analytic model describing mobile-immobile dislocation intersections applicable to fcc polycrystals. Based on previous work on dislocation drag, we estimate its temperature and pressure dependence and its effects on stress–strain rate relations. In the high temperature regime, we show that drag can remain the dominating effect even down to intermediate strain rates. We also discuss the consequences of having a limiting dislocation velocity, a feature which is typically predicted by analytic models of dislocation drag, but which is somewhat under debate because a number of MD simulations predict supersonic dislocations.

42 ENGINEERING↗
DOE OSTI2020

A predictive strain-gradient model with no undetermined constants or length scales

A general meso–scale (GM) crystal plasticity (CP) model was developed that accounts for lower-order (strain hardening) and higher-order (internal stress) effects of geometrically necessary dislocations (GNDs). It is predictive: no arbitrary parameters or length scales were invoked and no ad hoc numerical techniques were employed. It uses general stress field equations for GND content and a novel harmonization technique to enforce consistency of elastic long-range singular defect fields with applied elastic-plastic fields. The model facilitates implementation in commercial finite element programs without requiring special elements, special boundary conditions, or access to element shape functions. GM simulations confirmed, with improved accuracy, previously published predictions of the Hall-Petch effect, Bauschinger effect, and anelasticity. Previously unpredicted phenomena were also predicted: anelasticity and hysteresis for single Ta crystals and strain-hardening stagnation. The internal stresses (higher-order effect) dominate at large length scales, while at small length scales, the GND density hardening (lower-order effect) dominates. As a result, GM predicts that strain heterogeneity and consequent GND internal stresses are important factors in anelasticity.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Deformation, dislocation evolution and the non-Schmid effect in body-centered-cubic single- and polycrystal tantalum

A physically-informed continuum crystal plasticity model is presented to elucidate deformation mechanisms, dislocation evolution and the non-Schmid effect in body-centered-cubic (bcc) tantalum widely used as a key structural material for mechanical and thermal extremes. We show the unified structural modeling framework informed by mesoscopic dislocation dynamics simulations is capable of capturing salient features of the large inelastic behavior of tantalum at quasi-static (10 -3 s -1 ) to extreme strain rates (5000 s -1 ) and at low (77 K) to high temperatures (873 K) at both single- and polycrystal levels. Here we also present predictive capabilities of the model for microstructural evolution in the material. To this end, we investigate the effects of dislocation interactions on slip activities, instability and the non-Schmid behavior at the single crystal level. Furthermore, ex situ measurements on crystallographic texture evolution and dislocation density growth are carried out for polycrystal tantalum specimens at increasing strains. Numerical simulation results also support that the modeling framework is capable of capturing the main features of the polycrystal behavior over a wide range of strains, strain rates and temperatures. The theoretical, experimental and numerical results at both single- and polycrystal levels provide critical insight into the underlying physical pictures for micro- and macroscopic responses and their relations in this important class of refractory bcc materials undergoing large inelastic deformations.

36 MATERIALS SCIENCE↗