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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Crystal plasticity modeling of strain-induced martensitic transformations to predict strain rate and temperature sensitive behavior of 304 L steels: Applications to tension, compression, torsion, and impact

This paper advances crystallographically-based Olson-Cohen (direct γ → α’) and deformation mechanism (indirect γ→ε→α’) phase transformation models for predicting strain-induced austenite to martensite transformation. Here, the advanced transformation models enable predictions of not only strain-path sensitive, but also of strain-rate and temperature sensitive deformation of polycrystalline stainless steels (SSs). The deformation of constituent grains in SSs is modeled as a combination of anisotropic elasticity, crystallographic slip, and phase transformation, while the hardening is based on the evolution of dislocation density and explicit shifts in phase fractions. Such grain-scale deformation is implemented within the meso-scale elasto-plastic self-consistent (EPSC) homogenization model, which is coupled with the implicit finite element (FE) method to provide a constitutive response at each FE integration point for solving boundary value problems at the macro-scale. Parameters pertaining to the hardening and transformation models within FEEPSC are calibrated and validated on a suite of data including flow curves and phase fractions for monotonic compression, tension, and torsion as a function of strain-rate and temperature for wrought and additively manufactured (AM) SS304L. To illustrate the potential and accuracy of the integrated multi-level FE-EPSC simulation framework, geometry, mechanical response, phase fractions, and texture evolution are simulated during gas-gun impact deformation of a cylinder and quasi-static tension of a notched specimen made of AM SS304L. Details of the simulation framework, comparison between experimental and simulation results, and insights from the results are presented and discussed.

304L steels↗
DOE OSTI2022

Visualization and validation of twin nucleation and early-stage growth in magnesium

The abrupt occurrence of twinning when Mg is deformed leads to a highly anisotropic response, making it too unreliable for structural use and too unpredictable for observation. Here, we describe an in-situ transmission electron microscopy experiment on Mg crystals with strategically designed geometries for visualization of a long-proposed but unverified twinning mechanism. Combining with atomistic simulations and topological analysis, we conclude that twin nucleation occurs through a pure-shuffle mechanism that requires prismatic-basal transformations. Also, we verified a crystal geometry dependent twin growth mechanism, that is the early-stage growth associated with instability of plasticity flow, which can be dominated either by slower movement of prismatic-basal boundary steps, or by faster glide-shuffle along the twinning plane. The fundamental understanding of twinning provides a pathway to understand deformation from a scientific standpoint and the microstructure design principles to engineer metals with enhanced behavior from a technological standpoint.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A viscoelastic viscoplastic damage model (VPVD)

This report documents the development, implementation, and parameter estimation procedures for aviscoelastic-viscoplastic-viscodamage model. This model is being developed to improve the predictive accuracy of simulations involving plastic bonded explosives (PBX). A representative PBX is PBX9501, which is 95% by weight HMX, 2.5% estane binder and 2.5% nitroplasticizer. Mechanical deformation of this material is primarily accommodated by straining of the polymer binder and rearranging the relatively rigid HMX crystals.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2022

An Eulerian crystal plasticity framework for modeling large anisotropic deformations in energetic materials under shocks

Here, this paper demonstrates a novel Eulerian computational framework for modeling anisotropic elastoplastic deformations of organic crystalline energetic materials (EM) under shocks. While Eulerian formulations are advantageous for handling large deformations, constitutive laws in such formulations have been limited to isotropic elastoplastic models, which may not fully capture the shock response of crystalline EM. The present Eulerian framework for high-strain rates, large deformation material dynamics of EM incorporates anisotropic isochoric elasticity via a hypo-elastic constitutive law and visco-plastic single-crystal models. The calculations are validated against atomistic calculations and experimental data and benchmarked against Lagrangian (finite element) crystal plasticity computations for shock-propagation in a monoclinic organic crystal, octahydro-1,3,5,7-tetranitro-1,3,5,7 tetrazocine (β-HMX). The Cauchy stress components and the resolved shear stresses calculated using the present Eulerian approach are shown to be in good agreement with the Lagrangian computations for different crystal orientations. The Eulerian framework is then used for computations of shock-induced inert void collapse in β-HMX to study the effects of crystal orientations on hotspot formation under different loading intensities. The computations show that the hotspot temperature distributions and the collapse profiles are sensitive to the crystal orientations at lower impact velocities (viz., 500 m/s); when the impact velocity is increased to 1000 m/s, the collapse is predominantly hydrodynamic and the role of anisotropy is modest. The present methodology will be useful to simulate energy localization in shocked porous energetic material microstructures and other situations where large deformations of single and polycrystals govern the thermomechanical response.

42 ENGINEERING↗
DOE OSTI2021

Inelastic peridynamic model for molecular crystal particles

The peridynamic theory of solid mechanics is applied to modeling the deformation and fracture of micrometer-sized particles made of organic crystalline material. A new peridynamic material model is proposed to reproduce the elastic–plastic response, creep, and fracture that are observed in experiments. The model is implemented in a three-dimensional, meshless Lagrangian simulation code. In the small deformation, elastic regime, the model agrees well with classical Hertzian contact analysis for a sphere compressed between rigid plates. Under higher load, material and geometrical nonlinearity is predicted, leading to fracture. Finally, the material parameters for the energetic material CL-20 are evaluated from nanoindentation test data on the cyclic compression and failure of micrometer-sized grains.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Processing behavior evolution of recycled polypropylene: An integrated experimental and Computer-Aided engineering simulation study

Polypropylene (PP) comprises 21% of global plastics production and 18% of plastics waste, yet less than 1% of solid-waste PP is recycled in the United States (U.S.), representing significant environmental and economic challenges. Mechanical recycling, the most prevalent recycling method, subject's materials to thermomechanical stresses, which typically degrade polymer properties, affecting the quality of polymer products. This study replicates the impact of mechanical recycling through multiple extrusion cycles to examine the effects on PP's processing behavior. Dynamic scanning calorimetry (DSC) measurements showed stable melting behavior across all processing conditions, while crystallization analysis exhibited consistent shifts in kinetic parameters. Rheological characterization demonstrated progressive viscosity reductions through successive cycles, particularly pronounced at elevated reprocessing temperatures. Here, the integration of this experimental data into injection molding simulations showed that recycled PP maintains viable processing characteristics. Our findings establish quantitative correlations between processing history and material behavior, enabling optimization of processing parameters directly rather than relying on trial-and-error approaches. While these results reflect idealized recycling conditions with minimal contamination, they provide a framework for understanding fundamental property evolution during mechanical recycling.

42 ENGINEERING↗
DOE OSTI2021

A Review: Microstructural and Phase Evolution in Alloys during Extended Plastic Deformation

Materials in service and during processing are often subjected to plastic deformation. For multi-phase metallic alloys, simple geometric models and atomistic simulations suggest that two distinctive regimes in these materials’ evolution during deformation exist. At low strains, evolutions are often dominated by the kinetic roughening of interfaces, which results from the superdiffusive transport of matter in sheared crystals. Furthermore, at high strains and temperatures where thermal diffusion is sluggish, on the other hand, shearing- induced forced atomic mixing dominates these evolutions, resulting in significantly enhanced solubility. Distinguishing these two regimes is shown to provide a convenient framework for rationalizing and analyzing recent experiments and simulations on wear of layered structures in the low strain-regime and nonequilibrium phase co-existence and self-organization in highly immiscible or reactive alloy systems in the high-strain regime.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Parallel Dislocation Simulator

ParaDiS, or Parallel Dislocation Simulator, is a simulation tool that performs direct numerical simulation of dislocation ensembles, the carriers of plasticity, to predict the strength in crystalline materials from the fundamental physics of defect motion, evolution, and interaction. The code has been successfully deployed on high performance computing architectures and used to study the origins of strength and strain hardening for cubic crystals, the strength of micro-pillars, and irradiated materials at LLNL. The ParaDiS code has been successfully deployed on more than one hundred thousand CPU's with over ten million active degrees of freedom.

Bulatov, Vasily↗
DOE OSTI2021

USAMP Low-Cost Mg Sheet Component Development and Demonstration Project

The overarching objective of this USAMP project was to develop and demonstrate door panels made from magnesium (Mg) sheet with a cost penalty over conventional steel stampings of no more than $\$5.50$/kg saved. The technical approach integrated experiments with advanced computational tools based on Integrated Computational Materials Engineering (ICME) methods to develop new alloy chemistries and their thermomechanical processing that promise improved formability and lower forming temperatures. A penultimate task before finally forming the stampings was to incorporate actual microstructure into models that would enable formability simulations. This approach would, for the first time, account for individual magnesium grains moving in an anisotropic fashion unlike that for aluminum or steel that have isotropic properties upon which the current simulation tools are based. In separate activities, new coatings and lubricants to facilitate forming and improved corrosion protection and joining strategies, were developed to ensure that the door could be produced with stated product requirements. A technical cost model, which included parts production, assembly, and paint for a door specifically designed for Mg sheet, showed the cost penalty to be between $\$4.26$ to $\$6.60$/kg saved, which enveloped the project’s cost targets. The cost of the coated Mg sheet was identified as the key driver for the cost penalty. The mass of the Mg-intensive door was 7.9 kg, which was 54% less than the baseline steel door.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Low-Cost Magnesium Alloy Sheet Component Development and Demonstration Project

The overarching objective of this USAMP project was to develop and demonstrate door panels made from magnesium (Mg) sheet with a cost penalty over conventional steel stampings of no more than $5.50/kg saved. The technical approach integrated experiments with advanced computational tools based on Integrated Computational Materials Engineering (ICME) methods to develop new alloy chemistries and their thermomechanical processing that promise improved formability and lower forming temperatures. A penultimate task before finally forming the stampings was to incorporate actual microstructure into models that would enable formability simulations. This approach would, for the first time, account for individual magnesium grains moving in an anisotropic fashion unlike that for aluminum or steel that have isotropic properties upon which the current simulation tools are based. In separate activities, new coatings and lubricants to facilitate forming and improved corrosion protection and joining strategies, were developed to ensure that the door could be produced with stated product requirements. A technical cost model, which included parts production, assembly, and paint for a door specifically designed for Mg sheet, showed the cost penalty to be between 4.26 USD to 6.60 USD/kg saved, which enveloped the project’s cost targets. The cost of the coated Mg sheet was identified as the key driver for the cost penalty. The mass of the Mg-intensive door was 7.9 kg, which was 54% less than the baseline steel door.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Site Disorder Drives Cyanide Dynamics and Fast Ion Transport in Li 6 PS 5 CN

Halide argyrodite solid-state electrolytes of the general formula Li 6 PS 5 X exhibit complex static and dynamic disorder that plays a crucial role in ion transport processes. Here, we unravel the rich interplay between site disorder and dynamics in the plastic crystal argyrodite Li 6 PS 5 CN and the impact on ion diffusion processes through a suite of experimental and computational methodologies, including temperature-dependent synchrotron powder X-ray diffraction, AC electrochemical impedance spectroscopy, 7 Li solid-state NMR, and machine learning-assisted molecular dynamics simulations. Sulfide and (pseudo)halide site disorder between the two anion sublattices unilaterally improves long-range lithium diffusion irrespective of the (pseudo)halide identity, which demonstrates the importance of site disorder in dictating bulk ionic conductivity in the argyrodite family. Furthermore, we find that anion site disorder modulates the presence and time scales of cyanide rotational dynamics. Ordered configurations of anions enable fast, quasi-free rotations of cyanides that occur on time scales of 10 11 Hz at T = 300 K. In contrast, we find that cyanide dynamics are slow or frozen in Li 6 PS 5 CN when site disorder between the cyanide and sulfide sublattices is present at T = 300 K. We rationalize the observed differences in cyanide dynamics in the context of elastic dipole interactions between neighboring cyanide anions and local strain induced by the configurations of site disorder that may impact the energetic landscape for cyanide rotational dynamics. Through this study, we find that anion disorder plays a decisive role in dictating the extent and time scales of both lithium ion and cyanide dynamics in Li 6 PS 5 CN.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Kinetics of HCP-BCC Phase Transition Boundary in Magnesium at High Pressure

Under high pressures, many crystalline metals undergo solid–solid phase transformations. In order to accurately model the behavior of materials under extreme loading conditions, it is essential to understand the kinetics of phase transition. Using molecular dynamics simulations, this work demonstrates the feasibility of characterizing the speeds of a moving phase boundary using atomistic simulations employing a suitable empirical potential for single-crystal magnesium. The model can provide temperature- and tensorial stress-dependent velocity of a moving phase boundary as a rate-limiting contribution to the kinetics of phase transformation in continuum codes. Results demonstrate that a nonlinear interaction exists between plasticity and phase transition, facilitating a jump in the velocity of a moving phase boundary, facilitated by activated plastic deformation mechanisms.

36 MATERIALS SCIENCE↗
DOE OSTI2026

A Solid State Zwitterionic Plastic Crystal with High Static Dielectric Constant

The dielectric data in Figure 3, Figure 4, Figure S6 of the published paper was extracted from 2EOIMTSA-BDS-DATA .txt file. This file can be directly opened using a text file editor. It can also be imported to Excel/ Origin for further plotting and analysis. The G' and G'' in Figure 3 of the publihsed paper was plotted from data in file 2EOImTSA-temperature-sweep.xlsx. This file can be directly opend using Excel. The details of DFT simulations mentioned in Figure 2, Figure 7, and Figure S9 of the published paper are included in the DFT.zip file.

Huang, Zitan [Pennsylvania State University]↗
DOE OSTI2025

Direct sensitivity analysis on the parameterization of crystal plasticity models

Various methods for calibrating crystal plasticity finite element (CPFE) models lead to non-unique input parameter values, which subsequently introduce uncertainty in the predicted mechanical response. Sensitivity analysis (SA) conducted on crystal plasticity models is used to identify how variability in these parameters contribute to output uncertainty. Traditional SA on CPFE parameters uses simplified surrogate models to save computational time. However, the accuracy of the surrogate models depends on the quantity of training data used, and any modeling error can propagate into the SA results, potentially affecting their reliability. In this work, the elementary effects test (EET) method, a global SA technique using direct CPFE simulations was employed, and the results obtained were compared with the First Order Second Moment (FOSM) method. ExaConstit, an open-source GPU-enabled CPFE code, was used to perform the simulations and direct SA. The EET method was accurately able to capture the non-linear effects of all the input parameters on the output and is a valuable approach for reliably attributing parameter sensitivities in CPFE models. Based on the results, efficient strategies to perform future parameter calibration and SA are discussed. Additionally, the SA trends observed in different single crystal orientations closely mirrored the activity of the slip systems.

Elementary Effects Test↗
DOE OSTI2022

Discrete dislocation dynamics for crystal RVEs. Part 1: Periodic network kinematics

A novel implementation of the dislocation flux boundary condition in discrete dislocation dynamics is presented. The continuity of the individual dislocation loops in a periodic representative crystal volume (RVE) is enforced across the boundary of the RVE with the help of a dual topological description for representing dislocation line kinematics in two equivalent spaces representing the deforming crystal, the RVE and the unbounded crystal spaces. The former describes the motion of the dislocations in the simulated crystal RVE whereas the latter represents the motion of dislocations in an infinite space containing all replicas of the RVE. A mapping between the two spaces forms the basis of the implementation of flux boundary condition. The implementation details are discussed in the context of statistical homogeneity of bulk crystals undergoing macroscopically homogeneous plastic deformation. In this case, the boundary nodes associated with dislocation segments bear no relevance in the motion of the dislocations. Finally, some test cases are presented and discussed to establish the proposed approach.

42 ENGINEERING↗
DOE OSTI2021

A constitutive model incorporating grain refinement strengthening on metallic alloys

Surface nano-crystallization techniques have been recently developed as one of the most efficient ways to optimize materials’ structure, and therefore develop the local and global mechanical behavior as to increase strength without compromising ductility. Here, in this work, we present a constitutive model incorporating grain refinement hardening to simulate the nano-crystallization technique, specifically, surface mechanical attrition treatment. The computation is implemented using user-defined VUMAT subroutines. As an example of its application, a geometry model with full coverage of random impacts are employed. The results show that the model has rather precise predictability of grain size evolution during plastic deformation. The readily embedded with a computational code of material dynamics enables this novel model to be a promising tool to study the dynamic evolution of microstructures under plastic deformation.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Monotonic Gaussian Process for Physics-Constrained Machine Learning With Materials Science Applications

Physics-constrained machine learning is emerging as an important topic in the field of machine learning for physics. One of the most significant advantages of incorporating physics constraints into machine learning methods is that the resulting model requires significantly less data to train. By incorporating physical rules into the machine learning formulation itself, the predictions are expected to be physically plausible. Gaussian process (GP) is perhaps one of the most common methods in machine learning for small datasets. In this paper, we investigate the possibility of constraining a GP formulation with monotonicity on three different material datasets, where one experimental and two computational datasets are used. The monotonic GP is compared against the regular GP, where a significant reduction in the posterior variance is observed. The monotonic GP is strictly monotonic in the interpolation regime, but in the extrapolation regime, the monotonic effect starts fading away as one goes beyond the training dataset. Imposing monotonicity on the GP comes at a small accuracy cost, compared to the regular GP. The monotonic GP is perhaps most useful in applications where data are scarce and noisy, and monotonicity is supported by strong physical evidence.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Compressive strength of the plastic-bonded explosive surrogate idoxuridine under quasistatic and dynamic compression

Inert simulant materials, or “mocks”, are often used as surrogates for plastic-bonded explosives (PBX) in non-detonative tests in order to mitigate hazards. Mocks should reproduce as many non-detonative properties of the explosive as possible, including structural behavior in a variety of thermal and mechanical environments. Recently, the molecular crystal idoxuridine (IDOX) has been identified as an ideal mock for the main component in the explosive polymer-matrix composite PBX 9501, and has performed favorably under quasistatic loading conditions. Here, in order to assess robustness over a range of mechanical environments, plastic-bonded IDOX was compression tested from 0.001/s to 2000/s strain rates and compared to PBX 9501 historical data. Plastic-bonded IDOX showed good agreement to PBX 9501 across these strain rates, justifying continued development and production as a mock.

Burch, Alexandra C.↗