Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “plastic crystals”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
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 implementations of a dislocation density-based continuum model

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, which is referred to as advection. Large strain crystal plasticity models are based on the multiplicative decomposition of the total deformation gradient. Further, the stress is evaluated as a function of the thermoelastic strain, temperature, and other state variables. Advection of tensor quantities, such as the strain, is coupled with possibly significant advection errors. In an effort to reduce the advection errors, we develop two rate forms of an established dislocation density-based continuum model. To that end, we replace the multiplicative decomposition of the deformation gradient with the additive decomposition of the velocity gradient, and define the stress rate instead of the total stress. The Eulerian implementation is compared with Lagrangian calculations, and two numerical examples with severe deformation levels are presented.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Surface smoothing for laser powder-bed Ti-6Al-4V by a transient liquid phase

Surface roughness is the primary driver of fatigue for additively manufactured metals. To address surface roughness, this work introduces a new method to smooth features beyond line-of-sight without material removal. The method applies a coating that triggers local surface remelting by activating a eutectic reaction during heat treatment. The associated liquid phase then wets and isothermally solidifies into a smoother surface. For Ti-6Al-4V fabricated with laser powder bed fusion, samples with and without TLP smoothing (using a Cu coating) were characterized with a suite of techniques, including mechanical testing, electron backscatter diffraction, synchrotron X-ray tomography, and fractography. TLP smoothing reduced surface roughness by 80% and amplified compressive residual stress at the surface by about 50%. With statistically equivalent virtual microstructures, crystal plasticity scrutinized the roles of phases, porosity, and surface roughness. Although the tensile strain-to-failure was reduced to 1% strain, the TLP smoothing process increased high-cycle fatigue strength by about 20% compared to control samples, pointing to future opportunities to optimize the new process through various coating compositions and heat treatment schedules. Overall, this work establishes a new paradigm for treating surfaces of materials for smoothness and compressive residual stress.

Additive manufacturing↗
DOE OSTI2022

Creep anisotropy modeling and uncertainty quantification of an additively manufactured Ni-based superalloy

The advantages offered by additive manufacturing over traditional processes has driven a great deal of industrial and academic interest in recent years. However, the process is relatively new and requires additional investigation to become sufficiently mature for wide scale industrial adoption. Electron beam melting powder bed fusion is one technology that has shown promise for fabricating high temperature resistant materials such as nickel based superalloys. The resulting microstructures typically exhibit a strong fiber texture in the build direction giving rise to anisotropic time-dependent deformation behavior. In order to accelerate the qualification of these materials for industrial adoption accurate numerical models are needed for simulating their behavior. In this work a crystal plasticity model including non-Schmid effects is presented for capturing creep anisotropy observed in additively manufactured IN738LC. The model is calibrated via a probabilistic framework where model parameters are treated as random variables. An iterative sequential design strategy is utilized to efficiently identify the probability density of the unknown model parameters. As a case study the model is utilized to investigate the behavior of randomly oriented equiaxed grain clusters sometimes observed embedded in the additively manufactured columnar structure. A synthetic realization is simulated and uncertainty is propagated through to the full-field response. Results indicate that these features are the source of significant creep relaxation and strain accumulation which partially explains observed grain boundary decohesion at these locations.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Grain size effects on slip band development

Crystallographic slip localizations, known as slip bands, concentrate stress in polycrystals, often leading to the nucleation of damage. Slip band development has been experimentally shown to be sensitive to grain size, tending to develop more frequently and with a greater intensity in large grains. In this work, we investigate the influence of grain size on the propensity for crystallographic slip band development. To this end, we employ the slip band-fast Fourier transform method (SB-FFT). SB-FFT is a 3D, full-field crystal plasticity model that permits the incremental development of discrete crystallographic slip bands according to microstructure and material properties. We present a model Inconel 718 tricrystal to isolate the effect of grain size. Our findings show that slip bands in large grains develop at lower applied strain levels and at a faster rate than slip bands in small grains. The grain size effect is due to a backstress produced by the interaction of the slip band and its neighboring grain. The backstress is most intense at small grain sizes, impeding slip activity within a developing slip band and immediately surrounding the slip band.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Micromechanical Fields Associated with Irregular Deformation Twins in Magnesium

Understanding and controlling the development of deformation twins is paramount for engineering strong and stable hexagonal close-packed (HCP) Mg alloys. Actual twins are often irregular in boundary morphology and twin crystallography, deviating from the classical picture commonly used in theory and simulation. In this work, the elastic strains and stresses around irregular twins are examined both experimentally and computationally to gain insight into how twins develop and the microstructural features that influence their development. A nanoprecession electron diffraction (N-PED) technique is used to measure the elastic strains within and around a $\left\{ {10\overline{1}2} \right\}$ tensile twin in AZ31B Mg alloy with nm scale resolution. A full-field elasto-viscoplastic fast Fourier transform (EVP-FFT) crystal plasticity model of the same sub-grain and irregular twin structure is employed to understand and interpret the measured elastic strain fields. The calculations predict spatially resolved elastic strain fields in good agreement with the measurement, as well as all the stress components and the dislocation density fields generated by the twin, which are not easily obtainable from the experiment. The model calculations find that neighboring twins, several twin thicknesses apart, have little influence on the twin-tip micromechanical fields. Furthermore, this work reveals that irregularity in the twin-tip shape has a negligible effect on the development of the elastic strains around and inside the twin. Importantly, the major contributor to these micromechanical fields is the alignment of the twinning shear direction with the twin boundary.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Development of Predictive Model for Accurate Rupture Time from Multi-Axial Creep in Alloy 709 with Physics-Based Simulations

A physics-based model is developed to predict multiaxial creep behavior in Alloy 709 (A709), an advanced austenitic stainless steel intended for high-temperature applications such as Sodium Fast Reactors (SFRs). Compared to conventional stainless steels like 316H, A709 offers superior high-temperature performance; however, comprehensive data on its multiaxial creep response remain limited. To address this gap, a crystal plasticity finite element (CPFE) framework is used to simulate the deformation and failure mechanisms of A709 under multiaxial loading conditions. The model incorporates an extended Hu-Cocks dislocation creep formulation that accounts for precipitation effects, along with the Sham–Needleman model to capture grain boundary cavitation-driven failure. These advanced constitutive models enable a detailed understanding of the interplay between microstructural evolution and macroscopic creep response. Furthermore, the study evaluates the predictive accuracy of various effective stress measures in estimating creep rupture life, leveraging simulated multiaxial creep data. The findings provide critical insights into the applicability of different stress measures for engineering design and life prediction of A709 components operating under complex loading conditions. This work contributes to improving the reliability of high-temperature structural components by advancing predictive modeling capabilities for advanced austenitic steels.

Alloy 709↗
DOE OSTI2022

The effects of microstructure on deformation twinning in Mg WE43

The interplay between microstructure and deformation twinning in a WE43-T6 Mg alloy under uniaxial compression was investigated using a combination of scanning electron microscopy with digital image correlation (SEM-DIC) and crystal plasticity finite element (CPFE) simulation. To improve understanding of the statistical characteristics of deformation twin formation, microstructural effects were characterized in over 1000 grains through metrics including the nominal Schmid Factor, grain size, geometric compatibility factor (m'), residual Burgers vector, and the strain accommodated by neighboring grains. There was a strong correlation between the nominal Schmid Factor and both twin activation and variant selection, but this was not fully deterministic. Deformation twinning also exhibited a strong dependence on existing slip and twinning in the neighboring grain and on the m' value. Within the range of grain sizes present in this condition, grain size was determined to have minimal effect on deformation twinning. Finally, statistical analysis of CPFE simulations was used to further investigate microstructural effects on twinning, and qualitatively captured the effect of the nominal Schmid Factor.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Solving Stochastic Inverse Problems for Property–Structure Linkages Using Data-Consistent Inversion and Machine Learning

Determining process–structure–property linkages is one of the key objectives in material science, and uncertainty quantification plays a critical role in understanding both process–structure and structure–property linkages. In this work, we seek to learn a distribution of microstructure parameters that are consistent in the sense that the forward propagation of this distribution through a crystal plasticity finite element model matches a target distribution on materials properties. This stochastic inversion formulation infers a distribution of acceptable/consistent microstructures, as opposed to a deterministic solution, which expands the range of feasible designs in a probabilistic manner. Furthermore, to solve this stochastic inverse problem, we employ a recently developed uncertainty quantification framework based on push-forward probability measures, which combines techniques from measure theory and Bayes’ rule to define a unique and numerically stable solution. This approach requires making an initial prediction using an initial guess for the distribution on model inputs and solving a stochastic forward problem. To reduce the computational burden in solving both stochastic forward and stochastic inverse problems, we combine this approach with a machine learning Bayesian regression model based on Gaussian processes and demonstrate the proposed methodology on two representative case studies in structure–property linkages.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Grain boundary slip – twin transmission in titanium

Here, using a combination of mechanical testing, scanning electron microscopy, and a unified crystal plasticity framework for discrete intragranular shear localization, we investigate intense, localized slip bands on prismatic planes and {10$\bar{1}$2} $\langle$$\bar{1}$011$\rangle$ tensile twins, and their transmission across the grain boundaries (GBs) in commercially pure titanium. The analyses show that the orientation and curvature of the GB influence the local stress fields at the GBs, and consequently, the slip/twin transmission across the boundary. In addition to host grain properties, neighboring grain properties, such as active slip systems and instances of heterogeneity like slip bands and twins, heavily affect the deformation mechanisms in each grain. Finally, the applicability of geometric factors to predict the transmission in the experimentally observed co-located pairs is discussed. The local stress field calculated by the discrete slip and twin band model is shown to be capable of determining whether a transmission has occurred in an observed co-located pair, and also the direction of the transmission.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Synergistic effect of microstructure and defects on the initiation of fatigue cracks in additively manufactured Inconel 718

Fatigue cracks in additively manufactured (AMed) Inconel 718 (IN-718) in machined surface condition often initiate from persistent slip bands (PSBs) unlike other popular AM alloys such as 17-4 PH stainless steel or Ti-6Al-4V, where fatigue crack initiation is exclusively from volumetric defects; therefore, a competition between PSB- vs. defect- mediated crack initiation clearly exists. To shed light on the factors governing the competition, this study investigates the characteristics of cyclic strain localization, PSB formation, and crack initiation via crystal plasticity (CP) modeling of cyclic loading on polycrystalline aggregates which are then validated by experiments. In this work, a physics-based, free slip distance (FSD) dependent slip strength evolution law is proposed, which is shown to enable the CP model to simulate the heterogeneous strain distribution in IN-718. Implementing a crack initiation criterion based on strain contrast, the locations and lives for crack initiation can also be calculated. It is shown that both FSD and resolve shear stress influence the strain localization and crack initiation behaviors. The distribution of a localization parameter calculated based on the multiplication of FSD and Schmid factor within a grain is found to correlate well with the locations of PSBs. The maximum values of the localization parameters within a microstructure are shown to correlate well with the experimentally obtained crack initiation lives. The presence of volumetric defects in IN-718 generally do not impact strain localization behavior unless their size is large compared to the grain size.

42 ENGINEERING↗
DOE OSTI2021

Simulated effects of sample size and grain neighborhood on the modeling of extreme value fatigue response

Assessing the size of representative volume elements (RVEs) for fatigue-related applications is challenging. A RVE relevant to random microstructure requires a volume of material that is sufficiently large to capture the grain/phase heterogeneity that captures all statistical moments of the distribution of the driving force for fatigue crack formation at “hot spot” grains. Consequently, the large size of a microstructure RVE required to study fatigue phenomena is largely computationally intractable and difficult to explore. A more realistic objective in this work is to systematically study, as a function of the size of a statistical sample of microstructure, trends towards convergence of the simulated distribution of driving force for fatigue crack formation. Our present work accordingly leverages the recently developed open-source PRISMS-Fatigue framework to examine the trends in convergence of extreme value distributions (EVD) of Fatigue Indicator Parameters (FIPs) in progressively larger polycrystalline microstructure realizations of FCC Al alloy 7075-T6 using crystal plasticity finite element method simulations. The results are compared to the traditional method in which ensembles of statistical volume elements (SVEs) are simulated to build up statistics intended to approximate those associated with a larger volume of material. The convergence of EVDs with increase of size of a SVE of microstructure is closely related to the extent of grain nearest neighbor (NN) interactions. Accordingly, the sensitivity of the local micromechanical response at hot spot grains is quantitatively investigated by systematically varying the orientations of NN grains. Results indicate that SVEs with cubic crystallographic texture tend towards convergence of the EVD of FIPs with tens of thousands of grains while the random and rolled textures require larger volumes. Simple relationships based on microstructure parameters (e.g., Schmid Factor, grain size, NN misorientation) do not completely correlate to fatigue hot spot grains. Finally, the sensitivity of the extreme value fatigue response at hot spot grains extends to the 3rd NN when a single neighborhood grain orientation is altered.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Learning macroscopic internal variables and history dependence from microscopic models

This paper concerns the study of history dependent phenomena in heterogeneous materials in a two-scale setting where the material is specified at a fine microscopic scale of heterogeneities that is much smaller than the coarse macroscopic scale of application. Here, we specifically study a polycrystalline medium where each grain is governed by crystal plasticity while the solid is subjected to macroscopic dynamic loads. The theory of homogenization allows us to solve the macroscale problem directly with a constitutive relation that is defined implicitly by the solution of the microscale problem. However, the homogenization leads to a highly complex history dependence at the macroscale, one that can be quite different from that at the microscale. In this paper, we examine the use of machine-learning, and especially deep neural networks, to harness data generated by repeatedly solving the finer scale model to: (i) gain insights into the history dependence and the macroscopic internal variables that govern the overall response; and (ii) to create a computationally efficient surrogate of its solution operator, that can directly be used at the coarser scale with no further modeling. We do so by introducing a recurrent neural operator (RNO), and show that: (i) the architecture and the learned internal variables can provide insight into the physics of the macroscopic problem; and (ii) that the RNO can provide multiscale, specifically FE 2 , accuracy at a cost comparable to a conventional empirical constitutive relation.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Slip localization behavior at triple junctions in nickel-base superalloys

Incipient slip localization in the vicinity of hundreds of grain boundary triple junctions (TJs) in a lightly deformed nickel-base superalloy IN718 is studied using a combination of three-dimensional (3D) crystal plasticity finite element (CPFE) modeling, high resolution digital image correlation (HR-DIC) and 3D electron back-scatter diffraction tomography (3D EBSD). A 3D reconstruction method enables identification of thousands of TJs and correspondence of any observed slip bands with their originating TJ lines below the specimen surface. Here, we present a large-scale CPFE model of the experimental 3D microstructure composed of high-fidelity representation of the TJ lines and the boundaries and interiors of the parent grains and use it to calculate the local micromechanical response and slip activity of all TJs at the onset of macroscopic yielding. Statistical analysis of the calculated quantities reveal TJs develop larger stress concentration and grain-average re-orientation than grain interiors and grain boundaries, however no substantial differences in cumulative slip were found among these microstructural regions. We find that TJs with observed slip bands generate lower grain-average re-orientation, fewer active slip systems, and more localized slip on a single system than those without. The distinctions in the reorientation and slip activity are stronger in TJs that experience more intense slip.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A general kinetic framework for dislocation mobility derived from probabilistic cellular automaton simulations of the kink-pair mechanism

Dislocation mobility laws are essential components of dislocation-density-based crystal plasticity models. For dislocations governed by the kink-pair mechanism, however, existing formulations are typically limited to specific regimes due to the com plex interplay between stochastic kink-pair nucleation and lateral kink migration. In this work, we develop a general kinetic framework that expresses the average dislocation velocity as a function of mechanism-level variables: positive/negative kink pair nucleation rates, kink migration velocity, dislocation segment length, critical kink-pair width, and kink height. Probabilis tic cellular automaton simulations are used to capture the behavior of conceptual dislocation segments between the limiting conditions of migration outpacing nucleation on the one end and nucleation outpacing migration on the other. An elemen tary functional form that captures the system dynamics is suggested and fitted against the simulation results. This framework remains valid for arbitrary combinations of the six variables and is, therefore, compatible with any admissible constitutive re lations that describe their stress and temperature dependence. Comparisons with established approaches and experimental results confirm the robustness and physical consistency of the formulation, making it broadly applicable to material systems in which dislocation motion is governed by the kink-pair mechanism.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Simulated effect of defect volume and location on very high cycle fatigue of laser beam powder bed fused AlSi10Mg

This study quantifies the interaction between volumetric defect location and size on the very high cycle fatigue (VHCF) of laser beam powder bed fused (LB-PBF) AlSi10Mg. Crystal plasticity finite element method (CPFEM) simulations were used to investigate the effects of defect location and size on the driving force for crack initiation. The CPFEM model was calibrated against uniaxial and cyclic experimental data of LB-PBF AlSi10Mg. Defect characteristics were informed by experimental data from the specimens produced in various geometries to create realistic representative volume elements (RVEs) with equivalent volume fractions of defects. By embedding defects of varying sizes and locations within the RVEs, fatigue indicator parameters (FIPs) were calculated to analyze the impact of defects’ characteristics on fatigue performance. Different combinations of defect volume and locations were generated for various microstructure instantiations, providing insight into extreme value fatigue responses. Larger defect volumes located on free surfaces consistently generated the highest FIPs, suggesting defect size and boundary proximity intensify stress concentration effects. RVEs with multiple smaller defects produced lower FIPs than those with single large critical defects. These findings underscore the critical role of defect characteristics on fatigue life, providing a foundation for future predictive modeling in fatigue-sensitive AM applications.

AlSi10Mg↗
DOE OSTI2024

Stress localization investigation of additively manufactured GRCop-42 thin-wall structure

A full-field crystal plasticity (CP) framework is presented for the GRCop-42 alloy to study microscopic mechanical behavior and local stress heterogeneities. The microstructures of additively manufactured (AM) materials are often unique relative to conventionally processed materials, and the local thermal histories drive these differences during the build process. These thermal histories depend on the process parameters (laser power, scan speed, and scan strategy) and the part geometry. Prior research has shown that the mechanical properties of thin-walled structures can vary significantly with wall thickness due to changes in the thermal boundary conditions during manufacturing. It is, therefore, desirable to perform CP simulations based on the phenomenological constitutive model to predict the local mechanical responses induced by microstructural heterogeneities. This work generates representative microstructures based on experimentally collected grain information (i.e., texture) for grain scale stress analysis, and the material constitutive parameters are calibrated using the experimental mechanical testing data. Here, we specifically investigated the effect of crystallographic texture and grain morphologies on the size-dependent mechanical properties of AM GRCop-42. The selection of appropriate material properties for implementing an effective free surface boundary condition and the influence of adjacent buffer layers are also discussed. Analysis of local field results reveals a strong correlation between stress localization and the initial grain orientation. However, no significant relationship between the misorientation of the individual adjacent grains and the average misorientation is observed.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Prediction of softening kinetics and recrystallization texture in non-isothermally annealed bulged tubes using CPFEM and CA models

A hierarchically coupled cellular automata (CA) model, crystal plasticity finite element method (CPFEM), and thermal finite element (FE) model is developed to predict the softening kinetics of the bulged steel tube during non-isothermal annealing. Through the developed model, the kinetics of softening mechanisms including static recovery (SRV) and static recrystallization (SRX), as well as the recrystallization texture are predicted. Later, the Johnson-Mehl-Avrami-Kohnogorov (JMAK) model based on the predicted SRX data is developed to interpret the recrystallization behavior of the material. To perform this study, diverse experimental tests including tube hydroforming (THF), annealing, uniaxial tensile test, hardness test, as well as microstructure observations through optical microscopy and Electron Backscatter Diffraction (EBSD) tests on steel tube are performed. The obtained experimental data are utilized to calibrate and verify the implemented CPFEM model for simulation of THF process, thermal FE model for prediction of the local temperature over annealing time, and CA algorithm for modeling of the softening kinetics and texture evolution throughout the annealing process. The study shows that the predicted deformation characteristics, softening kinetics, recrystallization texture and temperature profile during non-isothermal annealing are in good agreement with experimental data. During the annealing process, a total of four stages for the kinetics of softening mechanisms is observed: No softening; SRV only; SRV dominant; and SRX dominant. During the progress of SRX, the behavior of recrystallization is abruptly changed, confirming that two different mechanisms are controlling the kinetics of transformation.

36 MATERIALS SCIENCE↗