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 181 records · Page 10

Quantitative study of the effect of grain boundary parameters on the slip system level Hall-Petch slope for basal slip system in Mg-4Al

Several theoretical studies have reported that the geometry and structure of grain boundaries in polycrystalline materials could impose a significant effect on the Hall-Petch slope. However, experimental observations are primarily limited by the ability of the techniques to accurately quantify the grain boundary strength and validate these theoretical models. Using high-resolution electron backscatter diffraction (HR-EBSD), the local stress tensor ahead of a slip band blocked by a grain boundary was quantified and coupled with a continuum dislocation pile-up model to assess the barrier strength of specific grain boundaries to specific slip systems, referred to as micro-Hall-Petch coefficient. For basal slip system in a deformed Mg-4Al alloy, the micro-Hall-Petch coefficient ($k^{basal}_{μ}$) varied significantly, from 0.054 to 0.184 MPa - m 1/2 for nine different grain boundaries. These results were correlated with geometric descriptors of the respective grain boundaries, with three-dimensional GB profile additionally measured via focused ion beam milling. It was found that the angle between the two slip plane traces on the grain boundary plane was the most sensitive parameter affecting $k^{basal}_{μ}$, followed by the angle between the slip directions. A functional form for calculation of $k^{basal}_{μ}$ depending on these two angles is proposed to augment crystal plasticity constitutive models with slip resistance dependent on some measure of the grain size. The method allows a new pathway to calibrate grain size strengthening parameters in crystal plasticity models, allowing further computational investigations of the interrelationship between texture, grain morphology, and the Hall Petch effect.

36 MATERIALS SCIENCE↗

Developing an experimental-computational framework to investigate the deformation mechanisms and mechanical properties of Al-8Ce-10Mg alloys at micro and macroscales

There is a promising future for the use of aluminum-cerium-magnesium alloys in a broad range of applications, including devices that operate at high temperatures. With cerium currently considered a waste product of rare earth mining validation studies of possible applications are essential to reduce the environmental waste. Here, a computational-experimental framework is developed to investigate the role both the intermetallic and matrix have on the mechanical properties of these alloys. A set of experiments, including SEM/EBSD imaging, nanoindentation, in-situ SEM tensile testing, and in-situ SEM-DIC tests are performed to characterize the microstructure and mechanical properties of these alloys. Furthermore, the elastic, plastic, and failure deformation mechanisms of the microstructure, and their correlation with the bulk scale mechanical properties are investigated. Experimental results are also used to calibrate parameters for a crystal plasticity finite element model, by performing a computational framework that minimizes the error between the computational and experimental results. This model is then utilized to investigate how the area percentage of intermetallics and the crystallographic texture control the mechanical properties of the alloy. Simulation results show that an increase in the percentage of intermetallics increase the strength but decrease the ductility of the alloy. Also, a change in material texture improves strength and reduces damage that leads to material failure. The development of the crystal plasticity model, as discussed in this work, opens opportunities for future investigations of similar aluminum-cerium-magnesium alloys.

36 MATERIALS SCIENCE↗

A new perspective on lone pair dynamics in halide perovskites

Halide perovskites form the foundation of an emerging class of materials for broad application in renewable and sustainable applications, including photocatalysis and solar energy harvesting. These materials exhibit beneficial photophysical properties, including bandgaps suitable for solar energy harvesting and efficient charge screening that underlies efficient charge carrier separation and resistance to defects. For organic–inorganic hybrid perovskites, these benefits are thought to arise, in part, from dipolar molecular cations that can reorient in response to charged particles and defects. In this work, we provide a similar perspective for inorganic metal halide perovskites, which do not contain molecular species with permanent dipoles. We discuss how lone pair electrons lead to dipolar ions that exhibit dynamics in analogy with traditional molecular plastic crystals and hybrid perovskites. We argue that further understanding these electronic plastic crystal motions with first principles simulations and synchrotron scattering can help create a basic understanding of photophysical properties of metal halide perovskites and inform the design of advanced functional materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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)↗

Modeling the non-Schmid crystallographic slip in MAX phases

We present a crystal plasticity constitutive relation for the description of experimentally observed non-Schmid crystallographic slip in a class of ternary carbides and nitrides commonly referred to as MAX phases. In the constitutive relation, we assume that the evolution of the slip system strength in MAX phases has two components – a classical component that depends on the Taylor cumulative shear strain and a non-Schmid component that depends on the stress normal to the slip plane. The non-Schmid crystal plasticity constitutive relation is then used to carry out finite element simulations of micropillar compression of single crystals of two MAX phases, Ti 2 AlC and Ti 3 AlC 2 . The finite element simulations not only quantitatively predict the stress – strain response of a wide range of crystallographic orientations of the micropillars but also rationalize the non-uniform deformation and the deformed shape of the micropillars observed in the experiments for the two materials. As a result, parametric studies are also carried out to quantify the role of the non-Schmid effect and understand the effects of key experimental parameters on the stress – strain response of the micropillars of the two MAX phases.

36 MATERIALS SCIENCE↗

Exploring the Role of Type-II Residual Stresses in A Laser Powder Bed Fusion Nickel-Based Superalloy Using Measurement and Modeling

Far-field high-energy X-ray diffraction microscopy (ff-HEDM) and the crystal plasticity finite element method (CPFEM) are used to investigate the role of grain-scale (Type-II) residual stresses on the fatigue life of additively manufactured (AM) Inconel alloy 625 (IN-625). Grain-averaged orientations, centroids, and residual elastic strain tensors from ff-HEDM data are used to instantiate a crystal plasticity model to simulate the effect of residual stresses at the grain scale. Simulation results indicate that the presence of tensile residual strains increase stress localization and heterogeneity within grains, triggering an earlier onset of plasticity. A microscale fatigue indicator parameter (FIP) is computed to model the impact of these residual strains on the cycles to fatigue crack nucleation. The crack nucleation model, based on the computed FIPs, predicts a significant reduction in the number of cycles for fatigue crack nucleation for mid- and high-cycle fatigue due to the residual strain induced localization, while the residual strains have minimal impact on low-cycle fatigue life.

Inconel↗

Multi-objective surrogate-assisted calibration of CPFEM models using macroscopic response and in situ EBSD measurements of grain reorientation trajectories

Crystal plasticity finite element method (CPFEM) models are widely used to simulate the deformation behaviour of polycrystalline materials, but their calibration is often limited by their high computational cost and the non-convexity of the optimisation landscape. Here, this study develops a multi-objective surrogate-assisted calibration workflow that couples a multi-objective genetic algorithm (MOGA) with an adaptively trained deep neural network (DNN) surrogate model to efficiently identify CPFEM parameters from experimental data. The workflow is demonstrated on three crystal plasticity (CP) formulations of increasing complexity — Voce hardening (VH), two-coefficient latent hardening (LH2), and six-coefficient latent hardening (LH6) — using in situ electron backscatter diffraction (EBSD) measurements of Alloy 617 under uniaxial tensile loading. The CPFEM models are calibrated against the experimentally observed stress–strain response and reorientation trajectories of eight grains, then validated against eight additional trajectories and overall texture evolution. Across the CP formulations, the macroscopic response was reproduced reliably, while differences emerged in the robustness and accuracy of the grain-scale predictions. Including grain reorientation trajectories in the multi-objective calibration improved texture evolution predictions and filtered out physically inconsistent parameter sets that can arise from calibrating against only the stress–strain data. The workflow also demonstrates good transferability of calibrated parameters from a low- to a high-fidelity microstructural model. These results provide practical guidance for integrating in situ microstructural data into CPFEM through efficient, repeatable, and physically meaningful multi-objective calibration.

Crystal plasticity finite element method↗

Effects of Zn and Ca on the deformation texture evolution of Mg–Zn–Ca alloys at elevated temperatures

This work utilized crystal plasticity simulations and experimental thermo-mechanical processing to examine how Zn and Ca content affects texture evolution in Mg–Zn–Ca alloys during simulated hot rolling. Four different compositions were studied: unalloyed Mg and three Mg alloys of ZX0p50 (Mg- 0.5 wt% Zn- 0.1 wt% Ca), ZX30 (Mg-3 wt% Zn- 0.1 wt% Ca), and ZX31 (Mg-3 wt% Zn- 0.3 wt% Ca). Multi-pass Gleeble experiments simulated hot rolling using plane strain conditions. Unalloyed Mg developed a strong basal texture, while the addition of Zn and Ca significantly weakened the texture in ZX30 and ZX31; this reduction was not seen in ZX0p50. Notably, ZX31 exhibited a split basal texture aligned with the rolling direction. Crystal plasticity simulations revealed that the weaker basal texture in ZX30 and ZX31 resulted from reduced activation of basal and twinning modes compared to unalloyed Mg and Zx0p50. Increased Zn and Ca content raised the ratios of basal to pyramidal critical resolved sheer stress (CRSS) and twin to pyramidal CRSS, enhancing pyramidal slip activity and weakening the basal texture. Extension twinning played a critical role in texture development during the first deformation pass. In unalloyed Mg and ZX0p50, twinning led to a strong basal texture, while ZX30 and ZX31 had weak basal textures that increased only slightly with further deformation. The split basal poles simulated in ZX31 were consistent with the experimental findings, highlighting the interplay between alloy composition and texture evolution.

Crystal plasticity finite element↗

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↗

A finite element formulation for deformation twinning induced strain localization in polycrystal magnesium alloys

Deformation twinning induces shear strain localization in hexagonal close-packed crystals and is critical for the material’s ductility and failure. Cracks often occur at twin-twin or twin-grain boundary intersections and propagate along twin bands. However, most crystal plasticity models for deformation twinning are based on a “pseudo-slip” approach and do not capture the localized deformation associated with the formation of each discrete twin band. The few exceptions are discrete twin models that involve very complex numerical algorithms and are often compromised in accuracy due to the numerical convergence. These factors make the discrete twin models hard to adopt. This paper proposes a modification to the conventional finite element weak form, to fully incorporate a twin-induced heterogeneous deformation that does not depend on the “pseudo-slip” assumption. The model starts by splitting the deformation gradient into elastic-slip-twinning components. The twin-induced deformation gradient component is computed separately by solving a microstructural evolution problem and then implemented into finite element weak form by constructing a global “twin-force” vector. The constitutive update (e.g., in the user-defined material subroutine, or UMAT, for ABAQUS) therefore avoids dealing with the twinning and recovers to the form of a regular slip-based crystal plasticity model. The results presented here indicate that the twin-induced strain localization and the associated stress-reversal phenomena near the twin band were naturally captured in the model, which was validated against an in-situ synchrotron X-ray micro-diffraction experiment.

36 MATERIALS SCIENCE↗

Multi-Scale Computational Modeling of Two-Phased Metal Using GMC Method

A multi-scale computational model for determining plastic behavior in two-phased CMSX-4 Ni-based superalloys is developed on a finite element analysis (FEA) framework employing crystal plasticity constitutive model that can capture the microstructural scale stress field. The generalized method of cells (GMC) micromechanics model is used for homogenizing the local field quantities. At first, GMC as stand-alone is validated by analyzing a repeating unit cell (RUC) as a two-phased sample with 72.9% volume fraction of gamma'-precipitate in the gamma-matrix phase and comparing the results with those predicted by finite element analysis (FEA) models incorporating the same crystal plasticity constitutive model. The global stress-strain behavior and the local field quantity distributions predicted by GMC demonstrated good agreement with FEA. High computational saving, at the expense of some accuracy in the components of local tensor field quantities, was obtained with GMC. Finally, the capability of the developed multi-scale model linking FEA and GMC to solve real life sized structures is demonstrated by analyzing an engine disc component and determining the microstructural scale details of the field quantities.

nickel alloys↗

Electron Backscatter Diffraction (EBSD) Analysis and U-Pb Geochronology of the Oldest Lunar Zircon: Constraining Early Lunar Differentiation and Dating Impact-Related Deformation

The evolution of the early moon was dominated by two processes (i) crystallization of the Lunar Magma Ocean (LMO) and differentiation of potassium-rare earth element-phosphorous-rich residual magma reservoir referred to as KREEP, and (ii) an intense meteorite bombardment referred to as lunar cataclysm . The exact timing of these processes is disputed, and resolution relies on collection and interpretation of precise age data. This study examines the microstructure and geochronology of zircon from lunar impact breccias collected during the Apollo 17 mission. A large zircon clast within lunar breccia 72215,195 shows sector zoning in optical microscopy, cathodoluminescence (CL) imaging and Raman mapping, and indicates that it was a relict fragment of a much larger magmatic grain. Sensitive high resolution ion microprobe (SHRIMP) U-Pb analysis of the zircon shows that U and Th concentration correlate with sector zoning, with darkest CL domains corresponding with high-U and Th (approx.150 and approx.100 ppm respectively), and the brightest-CL sectors containing approx.30-50 ppm U and approx.10-20 ppm Th. This indicates that variations in optical CL and Raman properties correspond to differential accumulation of alpha-radiation damage in each sector. Electron backscatter diffraction (EBSD) mapping shows that the quality of electron backscatter patterns (band contrast) varies with sector zoning, with the poorest quality patterns obtained from high-U and Th, dark-CL zones. EBSD mapping also reveals a deformation microstructure that is cryptic in optical, CL and Raman imaging. Two orthogonal sets of straight discrete and gradational low-angle boundaries accommodate approx.12 misorientation across the grain. The deformation bands are parallel to the crystallographic {a}-planes of the zircon, have misorientation axes parallel to the c-axis, and are geometrically consistent with formation by dislocation creep associated with <100>{010} slip. The deformation bands are unlike curved morphology of crystal-plastic microstructures in tectonically deformed terrestrial zircon, and geometrically similar to dislocation microstructures reported in experimentally shocked zircon. We interpret these crystal-plastic deformation microstructures to have resulted from a significant impact, either directly from impact shock, or during ductile flow directly following the impact. The deformation bands appear to continue undeflected through the non-indexed, radiation-damaged areas of the grain, which suggests that the orientation variation predates any significant mechanical weakening from radiation damage in the grain, and therefore occurred early in its history.

Timms, Nick↗

Role of crystallographic orientation on intragranular void growth in polycrystalline FCC materials

In this study, we study the effect of crystallographic orientation and applied triaxiality on the growth of intragranular voids. Two 3D full-field micromechanics methods are used, the dilatational visco-plastic fast-Fourier transform (DVP-FFT) and the crystal plasticity Finite Elements (CP-FE), both of which incorporate a combination of crystalline plasticity and dilatational plasticity. We demonstrate with several select cases that predictions of void growth from both formulations agree qualitatively. With the more computationally efficient DVP-FFT, additional effects of polycrystalline microstructure and the influence of nearest neighborhood are investigated. Crystals bearing a single intracrystalline void are studied in three types of 3D microstructural environments: isolated single crystals, individual equal-sized grains within a regular polycrystal, and individual variable sized grains within a polycrystal with grains and voids randomly located. We show that loading type plays a significant role. In strain-rate controlled conditions, voids in the hardest [111]-crystals grow the fastest in time, whereas in stress-controlled conditions, voids in the softest [100]-crystal grow the fastest in time. The analysis reveals that on average void growth is slower for the same starting orientation in the polycrystal than in the single crystal. We find that at the highest triaxiality tested that the correlation between crystal orientation and void growth rate in the polycrystal strengthens, drawing closer to that seen in the isolated single crystals. These results and model can help guide the microstructural design of polycrystalline materials with high strength and damage-tolerance in high-rate deformation.

36 MATERIALS SCIENCE↗

Sweep-tracing algorithm: in silico slip crystallography and tension-compression asymmetry in BCC metals

Abstract Direct Molecular Dynamics (MD) simulations are being increasingly employed to model dislocation-mediated crystal plasticity with atomic resolution. Thanks to the dislocation extraction algorithm (DXA), dislocation lines can be now accurately detected and positioned in space and their Burgers vector unambiguously identified in silico, while the simulation is being performed. However, DXA extracts static snapshots of dislocation configurations that by themselves present no information on dislocation motion. Referred to as a sweep-tracing algorithm (STA), here we introduce a practical computational method to observe dislocation motion and to accurately quantify its important characteristics such as preferential slip planes (slip crystallography). STA reconnects pairs of successive snapshots extracted by DXA and computes elementary slip facets thus precisely tracing the motion of dislocation segments from one snapshot to the next. As a testbed for our new method, we apply STA to the analysis of dislocation motion in large-scale MD simulations of single crystal plasticity in BCC metals. We observe that, when the crystal is subjected to uniaxial deformation along its [001] axis, dislocation slip predominantly occurs on the {112} maximum resolved shear stress plane under tension, while in compression slip is non-crystallographic (pencil) resulting in asymmetric mechanical response. The marked contrast in the observed slip crystallography is attributed to the twinning/anti-twinning asymmetry of shears in the {112} planes relatively favoring dislocation motion in the twinning sense while hindering dislocations from moving in the anti-twinning directions.

36 MATERIALS SCIENCE↗

(Invited) All-Solid-State Batteries Using Li 7 La 3 Zr 2 O 12 Garnet Electrolyte Framework

All-solid-state batteries (ASSB) using LLZO (variants of Li 7 La 3 Zr 2 O 12 ) garnet solid electrolyte with lithium metal anode potentially offer higher energy density and improved safety. Rational design of cell architecture as well as manufacturing scalability are key aspects to consider as the technology readiness level advances. A thick composite cathode layer attached to a thin solid electrolyte layer is desirable to obtain superior energy densities. Furthermore, the architecture within the composite cathode may be engineered to contain directional conduction paths of lithium ions or electrons for enhanced rate capabilities. Here, we demonstrate a functioning bulk-type LLZO based all-solid-state battery with a practical form factor incorporating the above described design concepts. Freeze-tape-casting (FTC), a scalable and environmentally friendly ceramic processing method, is used to construct 3D porous LLZO scaffolds composed of vertical arrays of LLZO walls. The thin solid-electrolyte layer is fabricated by tape-casting (TC). By sintering the stacks of FTC and TC green tapes, porous/dense bilayers and porous/dense/porous trilayers of LLZO frameworks are obtained. An ASSB was constructed using a porous/dense bilayer by infiltrating LiNi 0.6 Mn 0.2 Co 0.2 O 2 cathode powder and carbon black into the porous layer and adhering lithium metal foil to the dense side. We find it crucial to introduce a plastic crystal based soft solid electrolyte to the porous layer to electrochemically connect all cathode components, obviating the need for co-sintering to establish contact. The soft nature of the plastic crystal based solid-electrolyte may be able to accommodate the volume change of the cathode material as it cycles. In this device, which contained no liquid component, the initial discharge capacities were similar to those observed in a lithium ion battery configuration using the same cathode powder at C/10 rates.

25 ENERGY STORAGE↗

Making Plasticized Polymer Electrolytes Stable Against Sodium Metal for High‐Energy Solid‐State Sodium Batteries

Abstract Solid polymer electrolytes based on plastic crystals are promising for solid‐state sodium metal (Na 0 ) batteries, yet their practicality has been hindered by the notorious Na 0 ‐electrolyte interface instability issue, the underlying cause of which remains poorly understood. Here, by leveraging a model plasticized polymer electrolyte based on conventional succinonitrile plastic crystals, we uncover its failure origin in Na 0 batteries is associated with the formation of a thick and non‐uniform solid electrolyte interphase (SEI) and whiskery Na 0 nucleation/growth. Furthermore, we design a new additive‐embedded plasticized polymer electrolyte to manipulate the Na 0 deposition and SEI formulation. For the first time, we demonstrate that introducing fluoroethylene carbonate (FEC) additive into the succinonitrile‐plasticized polymer electrolyte can effectively protect Na 0 against interfacial corrosion by facilitating the growth of dome‐like Na 0 with thin, amorphous, and fluorine‐rich SEIs, thus enabling significantly improved performances of Na//Na symmetric cells (1,800 h at 0.5 mA cm −2 ) and Na//Na 3 V 2 (PO 4 ) 3 full cells (93.0 % capacity retention after 1,200 cycles at 1 C rate in coin cells and 93.1 % capacity retention after 250 cycles at C/3 in pouch cells at room temperature). Our work provides valuable insights into the interfacial failure of plasticized polymer electrolytes and offers a promising solution to resolving the interfacial instability issue.

Zou, Peichao↗

Making Plasticized Polymer Electrolytes Stable Against Sodium Metal for High‐Energy Solid‐State Sodium Batteries

Solid polymer electrolytes based on plastic crystals are promising for solid-state sodium metal (Na 0 ) batteries, yet their practicality has been hindered by the notorious Na 0 -electrolyte interface instability issue, the underlying cause of which remains poorly understood. Here, in this study, by leveraging a model plasticized polymer electrolyte based on conventional succinonitrile plastic crystals, we uncover its failure origin in Na 0 batteries is associated with the formation of a thick and non-uniform solid electrolyte interphase (SEI) and whiskery Na 0 nucleation/growth. Furthermore, we design a new additive-embedded plasticized polymer electrolyte to manipulate the Na 0 deposition and SEI formulation. For the first time, we demonstrate that introducing fluoroethylene carbonate (FEC) additive into the succinonitrile-plasticized polymer electrolyte can effectively protect Na 0 against interfacial corrosion by facilitating the growth of dome-like Na 0 with thin, amorphous, and fluorine-rich SEIs, thus enabling significantly improved performances of Na//Na symmetric cells (1,800 h at 0.5 mA cm −2 ) and Na//Na 3 V 2 (PO 4 ) 3 full cells (93.0 % capacity retention after 1,200 cycles at 1 C rate in coin cells and 93.1 % capacity retention after 250 cycles at C/3 in pouch cells at room temperature). Our work provides valuable insights into the interfacial failure of plasticized polymer electrolytes and offers a promising solution to resolving the interfacial instability issue.

25 ENERGY STORAGE↗

Experimental verification of a crystal plasticity-based simulation framework for predicting microstructure and geometric shape changes: Application to bending and Taylor impact testing of Zr

This paper is concerned with experimental verification of a recently developed multi-scale simulation framework for plastic deformation of metallic materials from quasi-static to impact deformation conditions. The framework is a visco-plastic self-consistent (VPSC) polycrystalline model embedded in an implicit finite element method (FE-VPSC) to provide a microstructure-sensitive constitutive response at each material point. Each material point of the FEM model is a polycrystalline aggregate with crystallographic deformation mechanisms operating at the single crystal scale with their evolving activity based on a dislocation density-based hardening law and texture. Four beams and three cylinders machined in different orientations from a textured plate of high-purity zirconium are tested quasi-statically in 4-point bending and at speeds of 100 m/s, 170 m/s and 243 m/s during Taylor impact tests, respectively. The variation in dimensional changes resulting from different sample orientations in the plate with respect to loading directions is measured for each sample. Moreover, texture and twinning characterization is performed using electron backscattered diffraction (EBSD). The deformation processes and underlying evolution of microstructure are successfully simulated using the FE-VPSC framework. In doing so, the model parameters are optimized and validated across a broad range of strain rates and temperatures. Simulation results in terms of geometrical changes and microstructural evolution are compared with the experimental measurements. The model predicts anisotropic material flow resulting from the hard-to-deform crystallographic directions, the development of gradients in texture and twinning through the geometries, tension–compression asymmetry, as well as the extent of plasticity under impact.

42 ENGINEERING↗