Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “polycrystal 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.

At least 37 records · Page 2

Microstructure Clones

Microstructure drives component behavior. Contemporary crystal plasticity studies compare strain measurements of polycrystal specimens to models. Because each specimen is unique, it is impossible to know which differences are significant. In this project, we invented microstructure clones and explored their use in understanding crystal plasticity. Microstructure clones are specimens with nearly identical microstructures, which allows for multiple destructive tests of a microstructure, insight into how a specimen will deform, variability quantification, and the ability to measure the effects of microstructural changes. Several sets of microstructure clones, pure nickel tensile bars, were tested. The techniques of digital image correlation, crystal plasticity finite element analysis, high resolution electron backscatter diffraction, transmission electron microscopy, and dislocation dynamics were used to understand the structural behavior of these microstructures. This work reshapes the fields of crystal plasticity and structure-property relationships by providing a technique to control for specific variables, quantify microstructural stochasticity, and replicate experiments.

36 MATERIALS SCIENCE↗

The effects of free surfaces on deformation twinning in HCP metals

Deformation twinning is a predominant mode of plastic deformation for hexagonal close packed metals, like Mg and Ti. The heterogenous microstructure and the local stresses associated with twinning play a key role in their mechanical response and fracture. Surface analyses, like electron microscopy, are frequently employed to spatially map microstructure and micromechanical fields in order to study twinning behavior. However, these measurements are inherently influenced by the vicinity of the free surface. In this work, an elasto-visco-plastic fast-Fourier-transform (EVP-FFT) polycrystal modeling approach is employed to investigate the effects of free surfaces on twin development before and after loading. We compare calculated micromechanical fields on free surfaces with those calculated inside the bulk and, in some cases, experimental surface measurements. The results indicate that the creation of free surfaces can promote twin propagation and growth and can influence twin morphology by causing a twin lamella to become larger, more blunted and irregular. The structure along the twin boundaries are also affected, due to the higher driving stresses that extend prismatic-basal and basal-prismatic facets. Furthermore, free surfaces invoke different slip activities in the twin and the surrounding parent crystal by enhancing basal, prismatic and pyramidal slip in some localized regions, while reducing slip in others. We demonstrate that the simulated free-surface effects lead to better qualitative and quantitative agreement with experimental measurements from scanning electron microscopy and digital image correlation.

36 MATERIALS SCIENCE↗

Microstructure Clones

Background: A material’s microstructure drives its material performance. Contemporary crystal plasticity experiments compare full-field strain measurements of polycrystal specimens to models. Because each specimen is unique, it is impossible to know which features of the observed deformation are deterministic vs statistical; thus, differences between model and experiment may or may not be significant. Objective: This paper introduces the invention of microstructure clones. Microstructure clones are 2D oligocrystal specimens that have nearly identical microstructures to remedy the aforementioned experimental limitations. Having specimens with nearly identical microstructures will allow for multiple destructive tests of a microstructure (either as repeats or intentionally different experiments), an ability to “see the future” by providing insight into how a specimen will deform, variability quantification, and experimental investigations of response to small microstructural changes. Methods: This work introduces microstructure clones. Repeatability of these clones is demonstrated in tensile bars of pure nickel. Local strain measurements from digital image correlation are compared between clone specimens and compared to results from a crystal plasticity finite element model. Results: Two sets of microstructure clones were tested in this study and displayed very consistent deformation responses within each clone set. Small observed differences in deformation invite investigation into microstructure stochasticity and the effect of small microstructural and loading differences. Conclusions: Microstructure clones represent a significant shift in understanding structure–property relationships. This work reshapes experimental crystal plasticity to allow for experiments that control for specific variables, quantification of microstructural stochasticity (and other sources of stochasticity), and opportunities for replicating experiments.

Crystal Plasticity↗

A generalizable machine learning-assisted fast Fourier transform algorithm to simulate the large strain phenomena in polycrystalline materials

Machine learning methods have shown initial promise in constitutive modeling for single crystals or homogenized polycrystals, delivering notable computational efficiency. However, existing machine learning-based constitutive models often lack generalizability, limiting their application across diverse boundary value problems. This study introduces a thermodynamics-informed artificial neural network model to accelerate rate-tangent crystal plasticity fast Fourier transform simulations for cross-scale deformation behaviors of polycrystals under complex loading. Our model integrates microstructural variability and local interactions effectively. To address local effects in each grain, we employ K-means clustering to group Gauss points within the microstructure into clusters assumed to be in similar mechanical states. This approach, based on self-clustering analysis, extends model scope from macroscopic stress response to the granular level, capturing mechanical responses and orientation evolution across grains. This reduces the number of nonlinear problems to solve, with cluster responses propagated throughout each group. The thermodynamics-based artificial neural network-extracted features are further processed using local material state clusters to account for history-dependent deformation and evolving microstructures. Additionally, representative volume element simulations with rate-tangent crystal plasticity fast Fourier transform provide reliable datasets for model training. The proposed model demonstrates high efficiency, accuracy, self-consistency, and enhanced generalizability in predicting strain–stress responses and orientation evolution at both individual grain and aggregate scales under complex loading conditions, such as biaxial tension and arbitrary loading scenarios.

36 MATERIALS SCIENCE↗

Evaluating the grain-scale deformation behavior of a single-phase FCC high entropy alloy using synchrotron high energy diffraction microscopy

We report although the deformation behavior of high-entropy alloys (HEAs) has been extensively studied at the macroscale, many important properties have yet to be explored for these alloys at the microscale, thus hampering accurate prediction of damage and failure processes. Synchrotron high-energy diffraction microscopy (HEDM) and fast-Fourier transform-based crystal plasticity modeling was conducted to investigate the three-dimensional (3D) grain-resolved micromechanical response for approximately 1,900 constituent grains within a single-phase FCC HEA up to 1% applied strain. The evolution of grain-resolved elastic strains, lattice reorientations, and maximum resolved shear stresses (mRSS) were evaluated to quantify elastic, yield, and fully plastic behavior. Overall, the initial critical resolved shear stress (CRSS), determined via in situ HEDM and companion modeling, was found to be > 20% higher than estimated using the classical polycrystalline Taylor factor (M = 3.06). However, a descriptive parameter based on the average grain-resolved Taylor factor ($\overline{M}$) was found to show excellent agreement with plastic yielding behavior observed within HEDM datasets. Noticeable deviations in HEDM lattice reorientations compared to both EVP-FFT simulations and classical predictions for FCC polycrystals were discovered, highlighting the complexity in correlating local lattice reorientations, Taylor, and Schmid factors with plastic response for this material at the grain-scale. Therefore, it is anticipated that the overall trends and parameter identification of 3D grain-resolved properties in this study can serve as an important foundation for continued mesoscale investigation on both well-established and newly developed Cantor-like HEAs.

36 MATERIALS SCIENCE↗

Comparative analysis of plasticity-based GND density estimation methods in crystal plasticity finite element models

In crystal plasticity finite element (CPFE) simulations, accurately quantifying geometrically necessary dislocations (GNDs) is critical for capturing strain gradients in polycrystals. We compare different methods for quantifying GNDs, all of which originate from the Nye tensor, which is computed as the curl of the plastic deformation gradient. The projection technique directly decomposes the Nye tensor onto individual screw and edge dislocation components to compute GNDs. This approach requires converting a nine-component Nye tensor into densities for a larger number of dislocation systems, a fundamentally underdetermined (non-unique) process, which is resolved using L2 minimization. In contrast, when employing CPFE analysis, one could directly compute dislocation densities on each slip system using shear gradients. Projection and slip gradient methods are compared with respect to their prediction of GNDs with changing grain size, strain, and grain neighborhoods, including multigrain junctions. Although these techniques match analytical GND densities for single slip, single crystal deformation, and are consistent with anticipated overall GND trends, we find that the GND densities from projection techniques are significantly lower than those predicted from CPFE-based slip gradients in polycrystals. A suggested improvement of only using the active dislocation systems in the projection technique almost entirely resolved this mismatch.

Crystal plasticity↗

Simulation of dynamic crystal plasticity with a Lagrangian discontinuous Galerkin hydrodynamic method

Here we present a new Lagrangian modal discontinuous Galerkin (DG) hydrodynamic method that supports a dynamic dislocation based crystal plasticity model for simulating the mechanical behavior of crystallographic materials, both single crystal and polycrystalline, under dynamic conditions. A modal DG approach is used to evolve fields relevant to conservation laws. These fields are approximated by Taylor series polynomials of varying degree. These polynomials describe macro-scale hydrodynamic behavior while their evolution is determined by evaluating the dynamic crystal plasticity model at material points within the element. The dynamic crystal plasticity model is sensitive to the time increment size, with too large of time increments leading to instability in the model. To mitigate this, the temporal evolution of the dynamic crystal plasticity model is achieved with the combination of a sub-incrementing scheme with Heun’s third-order time integration scheme, which is also used to temporally evolve the governing equations. The implementation of the dynamic crystal plasticity model within the DG framework is tested using a 2D approximation of the Taylor impact test with a single crystal material, using quadratic elements that have faces that can bend. In addition to the standard continuous material modeling, we propose a new simulation method that would represent the heterogeneous behavior of polycrystalline microstructures within an element by varying the position and material properties of the material points within the element. This method is demonstrated using random orientation distributions on materials points that are arranged in both structured and random configurations.

42 ENGINEERING↗

Demonstrate new plasticity models for doped UO 2 that capture dislocation mechanisms

In light water reactors, fuel vendors are investigating the use of dopants to modify the properties of UO 2 pellets, with the goal of improving pellet-cladding mechanical interactions during operation. Dopants are expected to ‘soften’ the pellets; that is, the doped pellets have higher plastic deformation than conventional UO 2 . This leads to a reduction in the severity of mechanical pellet-cladding interactions, helping to reduce the hoop strain on the cladding. By minimizing the strain exerted by the pellet on the cladding, it is anticipated that cladding performance under accident conditions can be enhanced (i.e., lowering the risk of burst during a LOCA). Dopants such as chromium (Cr) promote grain growth during pellet fabrication, leading to larger grains; therefore, understanding the link between chemistry, microstructure and mechanical deformation (enhanced creep rates) behavior of UO 2 is critical to helping operators further substantiate the benefits of doping UO 2 . Historically, the nuclear energy industry has relied on empirical models to make assessments of performance. Compared to empirical models, mechanistic physics-based models provide benefits, such as, fewer data points for validation and better extrapolation where experimental data is scarce or non-existent. In this report, Bayesian inference techniques have been applied to a previously developed lower length-scale-informed diffusional creep model. The objective is to i) infer lower-length-scale parameter distributions from available experiment and then ii) determine the uncertainties in the measurable quantity (in this case creep rates) after propagating the inferred lower length scale parameter uncertainties. The approach requires many evaluations of the model, which becomes computationally insurmountable; therefore, a neural-network model is trained to data obtained by sampling the full model over the most important parameters. This neural-network is then used in the Bayesian inference approach to determine probability distributions in the parameter values that represent the uncertainty in the model given what is known from the experiments (posterior). A significant reduction compared to conservative initial (prior) uncertainties is achieved through inference against the experimental data, demonstrating the efficacy of this approach. Furthermore, by accounting for uncertainties in the experimental conditions and sample non-stoichiometry, it is possible to resolve apparent discrepancies in experimental measurements within a self-consistent grain boundary (Coble) creep model that is sensitive to chemistry. This work has been written up and submitted to Nuclear Technology for a special issue on accelerated fuel qualification (AFQ). This uncertainty quantification (UQ) work not only improves the diffusional model, while accounting for uncertainty, but also establishes a framework which can readily be applied to the mechanistic models of dislocation deformation developed in this study. The most likely values from the Bayesian analysis are incorporated into our UO 2 diffusional creep model and a lower length scale-informed irradiation UO 2 creep mechanistic model to generate a dataset. This dataset has been provided to our INL collaborators for training an artificial neural network surrogate model, which will be implemented in the BISON fuel performance code to assess how the results differ from those currently obtained using a fully empirical model and that of using the nominal (uncalibrated) atomic scale parameters in our mechanistic model. Plastic deformation (creep and glide) in UO 2 is a complex phenomenon, governed by multiple underlying processes such as local defect concentrations, applied stresses, and microstructural characteristics. Consequently, there is a need for a meso-scale model with polycrystalline resolution capable of extrapolating to large grain sizes applicable to doped UO 2 , where data is limited and the model can help bridge the knowledge gap. By integrating atomistic data into the polycrystal LApx code, it becomes possible to predict dislocation climb and glide plasticity that simple analytical models cannot accurately represent. The application of atomic-scale data within LApx demonstrated the importance of climb and glide mechanisms in reproducing high-stress UO 2 behavior. Behaviors such as this are crucial to capture and implement in BISON, as parts of the fuel pellet can reach temperatures where glide can occur before pellet cracking. This model which captures dislocation based mechanisms for UO 2 is then used to stand up the doped model accounting for larger grain sizes. It was found that larger grain sizes can lead to enhanced deformation rates in the glide regime, and therefore can help with the pellet cladding mechanical interaction. Therefore if the fuel pellet reaches conditions (stress/temperature) where glide is active, the enhanced creep rates for larger grains in the glide regime (doped UO 2 ) can help with pellet cladding mechanical interactions. Plastic deformation in UO 2 involves multiple mechanisms, including diffusional creep, dislocation climb, and glide. This milestone contains two parts: (1) UQ of a pre-existing lower length scale informed mechanistic diffusional creep model, and (2) development of a new LApx based model for dislocation-mediated creep mechanisms in UO 2 , with application to large-grain doped UO 2 .

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Micromechanics of slip bands on a free surface

A micromechanics analysis for the formation and propagation of slip bands on the free surface of a polycrystal under monotonic loading is presented. For the growth of slip bands, the analysis satisfies the conditions of both equilibrium and displacement continuity, as well as the relation between slip and the resolved shear stress throughout the polycrystal. Numerical calculations show how the microstress field causes the concentration of plastic deformation in discrete sliding bands and give results which are in good qualitative agreement with known slip band observations on aluminum single crystals.

Lin, S. R.↗

Multiscale Modeling of Structurally-Graded Materials Using Discrete Dislocation Plasticity Models and Continuum Crystal Plasticity Models

A multiscale modeling methodology that combines the predictive capability of discrete dislocation plasticity and the computational efficiency of continuum crystal plasticity is developed. Single crystal configurations of different grain sizes modeled with periodic boundary conditions are analyzed using discrete dislocation plasticity (DD) to obtain grain size-dependent stress-strain predictions. These relationships are mapped into crystal plasticity parameters to develop a multiscale DD/CP model for continuum level simulations. A polycrystal model of a structurally-graded microstructure is developed, analyzed and used as a benchmark for comparison between the multiscale DD/CP model and the DD predictions. The multiscale DD/CP model follows the DD predictions closely up to an initial peak stress and then follows a strain hardening path that is parallel but somewhat offset from the DD predictions. The difference is believed to be from a combination of the strain rate in the DD simulation and the inability of the DD/CP model to represent non-monotonic material response.

Saether, Erik↗

Multitask graph neural networks for elastoplastic response prediction in dual-phase polycrystals

Microstructure-sensitive prediction of elastoplastic response remains a recurring bottleneck in multiscale damage and fatigue modeling, where large ensembles of statistically distinct polycrystals are required to quantify variability and extreme-value behavior. In this work, we develop a multitask graph neural network (GNN) surrogate that maps dual-phase ferrite–martensite polycrystal microstructures to Statistical Volume Element (SVE)-level elastoplastic Quantities of Interest (QoIs). Each SVE is represented as a grain-adjacency graph, with node features encoding phase, geometry, and crystallographic orientation, and edge features encoding relative misorientation. A message-passing graph convolution generates node embeddings, which are pooled into a graph representation and passed to a multitask regression head that jointly predicts 10 scalar QoIs and vector-valued stress–strain responses in orthogonal loading directions across multiple martensite volume fractions and SVE sizes. Results show high accuracy for scalar QoIs and strong agreement for full stress–strain trajectories, with population envelopes reproducing both median behavior and finite-SVE variability across compositions and partition scales. A unified model trained on pooled volume-fraction data preserves most within-regime accuracy relative to regime-specific models while also capturing the broader cross-regime variation reflected in the pooled test set. Distributional comparisons further demonstrate that the surrogate preserves heterogeneity under SVE partitioning, enabling statistically consistent block-wise random-field construction for mesoscale analyses. Overall, the proposed grain-graph surrogate provides a practical pathway to accelerate ensemble-based studies of SVE-level constitutive variability in dual-phase polycrystals.

Crystal plasticity↗

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↗

Phase-field framework with constraints and its applications to ductile fracture in polycrystals and fatigue

Abstract Modeling of ductile fracture in polycrystalline structures is challenging, since it requires integrated modeling of cracks, crystal plasticity, and grains. Here we extend the typical phase-field framework to the situations with constraints on the order parameters, and formulate two types of phase-field models on ductile fracture. The Type-I model incorporates three sets of order parameters, which describe the distributions of cracks, plastic strain, and grains, respectively. Crystal plasticity is employed within grain interiors accommodated by J 2 plasticity at grain boundaries. The applications of the Type-I model to single crystals and bicrystals demonstrate the influences of grain orientations and grain boundaries on crack growth. In the Type-II model, J 2 plasticity is assumed for the whole system and grain structures are neglected. Taking advantage of the efficiency of the fast Fourier transform, our Type-II model is employed to study low cycle fatigue. Crack closure and striation-like patterning of plastic strain are observed in the simulations. Crack growth rate is analyzed as a function of the J -integral, and the simulated fatigue life as a function of plastic strain agrees with the Coffin–Manson relation without a priori assumption.

36 MATERIALS SCIENCE↗

Statistical evaluation of microscale stress conditions leading to void nucleation in the weak shock regime

Here, we investigate the heterogeneity of the stress state driven by anisotropic deformation response at the single crystal level through five statistical volume element (SVE) calculations of polycrystalline BCC tantalum. This work focuses on grain boundaries as a prominent material defect type prone to void nucleation based upon experimental observations of predominantly intergranular void nucleation in this material. The SVEs are constructed to be statistically representative of larger volumes of material and are meshed such that mean and standard deviation of grain size and orientation information is reconstructed. The computational meshes feature hexahedral (brick) elements and smooth conformal grain boundaries where significant stress concentration is known to occur, a tail effect of interest in the extreme events process of dynamic ductile damage. An existing micromechanical crystallographic plasticity model shown to capture the single crystal behavior of BCC tantalum well is used to perform the polycrystal calculations. The model includes representation of the non-Schmid effect of non-planar screw dislocation kinetics in tantalum. A three-dimensional stress state time profile predicted by damage modeling of a flyer plate impact experiment is applied as boundary conditions to each SVE. Resulting grain boundary stress state statistics are strongly non-Gaussian. Significant structural evolution is observed within the compressive hold before unloading into tension in the stress profile. Strong angular dependence of grain boundary traction magnitude with shock direction is observed. Non-Schmid effects continue to suggest their influence on propensity of microstructural defect types to nucleate voids. A general void nucleation criterion is proposed using probability theory. The general framework is specified to polycrystalline BCC tantalum in the weak shock regime to include the SVE calculations and literature molecular dynamics calculations of grain boundary void nucleation strength. Probability density functions (PDFs) are used to describe the interaction between the local stress state heterogeneity and the distributed grain boundary void nucleation strength state. A causation entropy maximization procedure removes the requirement for ad hoc selection of a PDF functional form and provides a rigorous procedure for data-based PDF determination. The resulting physically informed PDF describes the spatial appearance frequency of nucleated voids as a function of applied macroscale pressure. Lower length scale physics are thus packaged in a precise and computationally efficient way to provide computational plasticity insight to macroscale dynamic ductile damage models.

36 MATERIALS SCIENCE↗

Atomistic and phase field simulations of three dimensional interactions of {$10\bar12$} twins with grain boundaries in Mg: twin transmission and dislocation emission

In polycrystals, the interaction of dislocations and twins with grain boundaries (GBs) plays a role in hardening and formability during plastic deformation. While dislocation-GB interactions are relatively well-understood, twin-GB interactions remain mostly unknown. In this work, an approach using molecular dynamics and phase-field simulations is followed to study the forward and lateral interactions between {$10\bar12$} twins and tilt grain boundaries in Mg. Molecular dynamics results show that the resolved shear stress on slip/twinning modes of the neighboring grain, not the geometric alignment, is the dominant factor in determining the outcome of the twin-GB interactions. For some lateral interaction configurations, as the misorientation angle increases, the resolved shear stress on the same {$10\bar12$} twin variant of the neighboring grain reduces while it increases for slip or I 2 stacking fault emissions or other twin modes such as {$11\bar12$} and {$10\bar11$}, explaining why twin transmission is not seen at high misorientation angles. Furthermore, lateral and forward interactions of the twin with tilt grain boundaries whose misorientation axes are normal to the coherent twin boundary show significantly different outcomes. For the forward interaction, the twin is absorbed and stacking faults are emitted when interacting for low misorientation angles (up to 30°) while the lateral interaction results in twin transmission, nucleation of a {$11\bar12$} twin, and emission of I 2 stacking faults. Finally, comparisons between twin interactions with symmetric and asymmetric tilt GBs with different GB structures show similar outcomes.

36 MATERIALS SCIENCE↗

Chapter 4: Physically informed deep learning networks for simulating microstructure evolution of 3D polycrystals

As discussed in the previous chapter, high energy diffraction microscopy (HEDM) is used to study the micromechanical evolution of a material during in situ loading. HEDM experiments have been used to verify crystal plasticity (CP) simulations [119, 91, 90, 120], for experimental planning, material design, and to further analyze experimental results. However, Fast Fourier transform-based CP (CP-FFT) or finite element-based CP (CP-FE) methods are often too slow to be used in real-time during an experiment. CP-FFT is faster than CP-FE simulations due to the absence of meshing, but can still take hours to simulate the response of a single volume depending on the size and number of strain steps [127]. Reducing computation time would create a larger exploration space in planning and design, and enable faster analysis of experimental results and real-time feedback during an experiment. This research expands upon previous works to develop a workflow for predicting the full-field evolution of a 3D polycrystal. The workflow is simplified from previous works to predict only orientation and elastic strain tensors (from which stress tensors are calculated). The network is physically informed through loss functions and network architecture for a more robust model. The orientation predictions are informed about the cubic crystal symmetry of the material by incorporating disorientation and misorientation information into the network architecture and loss. The Von Mises stress is used to enforce the correct stress-strain trends in the strain tensor predictions. Additional total strain steps from the elastic and elastoplastic region are included to better capture the stress-strain evolution at smaller total strain steps. Material and hardening parameters are additional inputs into the networks to further inform the network and to study the network’s ability to predict different materials other than those used for training.

36 MATERIALS SCIENCE↗

LANL end of the year report (2023)

This report summarizes the advances made by LANL in the XMAT-H2 program over FY23. In summary LANL was tasked with leading efforts related to (1) the development of models predicting hydrogen (H) redistribution within polycrystals (H redistribution model), (2) the development of models for the mechanical performance of metals/steels under H environments (Model for H effects on plasticity), The development of a framework to quantify the uncertainty in model parameters (Uncertainty in model parameters), (3) small scale characterization of the mechanical response of 347H (Microstructure characterization).

08 HYDROGEN↗

Experiment and non-local crystal plasticity finite element study of nanoindentation on Al-8Ce-10Mg alloy

Cerium and magnesium strengthened aluminum alloys, or Al-Ce-Mg, is a recently developed alloy family that exhibits good mechanical properties at elevated temperatures (~300 °C). To examine the single-crystal properties of Al-Ce-Mg alloys, nanoindentation experiments are conducted in this study. A crystal plasticity finite element model (CPFEM) with the evolution of geometrically necessary dislocations (GNDs) is applied to simulate the indentation in individual grains. A parametric study is carried out to investigate the sensitivity of each crystal plasticity model parameter to the indentation behavior. The highly sensitive parameters are calibrated by matching the indentation load-depth curves, while the rest parameters are obtained from bulk polycrystal uniaxial tension tests. Overall, satisfactory matching between experiment and simulation is obtained for each individual grain. The calculated hardness, as determined from the experiment, shows the dependence on indent depth, which is captured by the GND model. Furthermore, the effect of grain orientation and neighboring grains to nanoindentation behavior have been discussed with the comparison between the simulation and experiments.

36 MATERIALS SCIENCE↗