Data-driven plastic anisotropy characterization using crystal plasticity and deep learning models
Explore the source record for details and available documents.
SEARCH · Engineering Papers
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.
Explore the source record for details and available documents.
Dislocation‐based dissipation mechanisms potentially control the viscoelastic response of Earth's upper mantle across a variety of geodynamic contexts, including glacial isostatic adjustment, postseismic creep, and seismic‐wave attenuation. However, there is no consensus on which dislocation‐based, microphysical process controls the viscoelastic behavior of the upper mantle. Although both intergranular (plastic anisotropy) and intragranular (backstress) mechanisms have been proposed, there is currently insufficient laboratory data to discriminate between those mechanisms. Here, we present the results of forced‐oscillation experiments in a deformation‐DIA apparatus at confining pressures of 3–7 GPa and temperatures of 298–1370 K. Our experiments tested the viscoelastic response of polycrystalline olivine—the main constituent of the upper mantle—at stress amplitudes from 70 to 2,800 MPa. Mechanical data are complemented by microstructural analyses of grain size, crystallographic preferred orientation, and dislocation density. We observe amplitude‐ and frequency‐dependent attenuation and modulus relaxation and find that numerical solutions of the backstress model match our results well. Therefore, we argue that interactions among dislocations, rather than intergranular processes (e.g., plastic anisotropy or grain boundary sliding), control the viscoelastic behavior of polycrystalline olivine in our experiments. In addition, we present a linearized version of the constitutive equations of the backstress model and extrapolate it to conditions typical of seismic‐wave propagation in the upper mantle. Our extrapolation demonstrates that the backstress model can explain the magnitude of seismic‐wave attenuation in the upper mantle, although some modification is required to explain the weak frequency dependence of attenuation observed in nature and in previous experimental work.
The change in volume, and therefore the change in mass density, of an aluminum alloy was measured in uniaxial tension using clip-on extensometers. The experimental data do not agree with the assumption of plastic incompressibility found in the classical theories of plasticity. In fact, the elastic and plastic volume changes are of the same order of magnitude. Plastic anisotropy is thought to be the prime cause of this plastic compressibility.
The theory of elasto-plastic flow is developed so that it may accommodate features such as work-hardening, anisotropy, plastic compressibility, non-continuous loading including local or global unloading, and others. A complete theory is given in quasi-linear form; as a result, many useful attributes are accessible. Several integral theorems may be written, finite deformations may be incorporated, and efficient methods for solving problems may be developed; these and other aspects are described in some detail. The theory is reduced to special forms for 2-space, and extensive experience in solving such problems is cited.
Modeling the evolution of voids during plastic flow as well as their effects on plastic dissipation is critical for both component manufacturing and lifetime estimation purposes. To this end, we propose a rate-dependent constitutive model to homogenize the effects of semi-randomly distributed voids on single crystal plasticity whilst capturing void interaction and plastic anisotropy. Here, this present work focuses on the case of face centered cubic crystals to introduce an anisotropic gauge function applicable within the crystal plasticity formalism. The approach combines analytical methods to describe the micromechanics of the system in combination with symbolic regression to capture analytically intractable mechanisms from data. The hybrid framework uses a physics-informed genetic programming-based symbolic regression algorithm to solve a multiform optimization problem simultaneously producing a new gauge function and a new strain rate equation. This is also a multi-objective optimization problem with many competing objectives. A new search and selection step is introduced to the genetic algorithm that promotes convergence toward a global solution that better satisfies all the objectives. Overall, the symbolic equations produced leverage data-driven methods to achieve greater accuracy than comparable alternatives on an analytically intractable problem while maintaining model transparency.
Generative modeling has opened new avenues for solving previously intractable materials design problems. However, these new opportunities are accompanied by a drastic increase in the required amount of training data. This is in stark juxtaposition to the high expense and difficulty in curating such large materials datasets. In this work, we propose a novel framework for integrating generative models within an active learning loop. Further, this enables the training of generative models with datasets significantly smaller than what has previously been demonstrated, providing a direct route for their application in data constrained environments. The functionality of this framework is then demonstrated by addressing the challenge of designing polycrystalline textures associated with target anisotropic mechanical properties. The developed protocol exhibited a cost reduction between 14 to 18 times over a randomly sampled experimental design.
Mechanical properties of fiberglass anisotropic laminated plastics - tensile, compression, impact and shear tests
Crack growth in ductile single crystals (DuSCs) is orientation dependent due to the anisotropies of crystal plasticity and elastic tensor. This study develops a phase-field model incorporating both crystal plasticity and crack growth and proposes a general method to decompose the elastic energy into compressive and tensile parts to prevent crack growth under compression in the phase-field description. The phase-field model, in combination with three Euler angles, is employed to simulate orientation-dependent crack growth in DuSCs. The contributions from crystal plasticity and anisotropic elasticity are compared, and the former is found to dominate in the anisotropy of crack growth in copper single crystals. Furthermore, the simulation results demonstrate that crystal orientation strongly affects the heterogeneous distribution of plastic strain and the interaction between plastic strain and crack growth. High-throughput phase-field simulations are performed with exhaustive crystal orientations, and the results are explained based on the anisotropy of the Taylor factor.
The influence of the internal structure at micrometer length scales on the deformation of polycrystalline materials can be effectively captured using crystal plasticity finite element methods (CPFEM). However, the complexity and nonlinearity of the deformation equations CPFEM solves demand significant computational power and resources to achieve accurate predictions, limiting its broader application. To address this challenge, we have identified a reduced-order representation of the complex data in order to establish a computationally efficient reduced-order models (ROM) and drastically reduce the computational expense of CPFEM. Specifically, in this work, we developed a parametric, data-driven, and non-intrusive ROM framework for CPFEM using proper orthogonal decomposition (POD) and sparse variational Gaussian process (SVGP) regression for single-crystal microstructures under tensile loading conditions. The developed protocol enables one to compress field into a latent/low-dimensional space described by principal component analysis (PCA) via the singular value decomposition (SVD) algorithm. As a result, the high-dimensional data are reduced to a significantly smaller amount of dimensions with POD bases and POD coefficients. Furthermore, we deployed an ensemble of SVGPs—extended from the classical Gaussian process (GP) regression for scalability and handling big data—in a massively parallel manner to train and predict latent POD coefficients using known POD bases from a set of previously obtained simulations results. Lastly, using the predicted POD coefficients, we reconstructed the full-field results and showed reasonable agreement compared with the true values obtained from running CPFEM. The developed framework is validated with a set of CPFEM simulations of a single embedded void in single-crystal aluminum alloy. While the framework is broadly applicable, this work specifically focuses on single-crystal microstructures, a single load case (e.g., tensile), and a specific void geometry (spherical).
Magnetic effects associated with shock induced transformation of fcc antiferromagnetic iron precipitates in polycrystalline copper disks at levels up to 5 GPa in weak magnetic fields (H not greater than 0.5 Oe) were investigated. The demagnetization and anisotropy associated with second order transition, the effects of plastic deformation in imparting magnetic anisotropy and magnetic hardening, and the influence of post shock thermal transients on magnetization associated with recovery, recrystallization and grain growth were studied. It was found that on the microsecond time scale of the shock induced first order transformation, the field sense is recorded in the transformed iron particles. For a given particle size the degree of transformation of fcc iron depends on the level of the shock. For a given shock level the resultant magnetic properties depend on the particle size distribution, with maximum effects noted in specimens with 400 to 600 A particles.
Not Available
Explore the source record for details and available documents.
Advanced 3-D inelastic structural/stress analysis methods and solution strategies for more accurate yet more cost-effective analysis of components subjected to severe thermal gradients and loads in the presence of mechanical loads, with steep stress and strain gradients are being developed. Anisotropy, time and temperature dependent plasticity and creep effects are also addressed. The approach is to develop four different theories, one linear and three higher order theories (polynomial function, special function, general function). The theories are progressively more complex from linear to general function in order to provide streamlined analysis capability with increasing accuracy for each hot section component and for different parts of the same component according to the severity of the local stress, strain and temperature gradients associated with hot spots, cooling holes and surface coating cracks. To further enhance the computational effectiveness, the higher order theories will have embedded singularities (cooling passages, for example) in the generic modeling region. Each of the four theories consists of three formulation models derivable from independent theoretical formulations. These formulation models are based on: (1) mechanics of materials; (2) special finite elements; and (3) an advanced formulation to be recommended by the contractor.
NASA Langley Research Center (LaRC) is applying flow forming technology to aerospace aluminum alloys for the manufacture of integrally stiffened cylinders that meet aircraft fuselage or launch vehicle specifications. High-strength 2xxx series alloys, in various temper/product form conditions, are being compared with medium-strength, very formable AA 6061 as the benchmark. Screening of candidate materials traditionally relies on empirical forming trials, but more recently on metrics derived from constant nominal strain rate tensile testing. In this work, the established hierarchy of formability metrics is adopted as the baseline; tensile area reduction, work hardening exponent, and modulus of resilience. The suite of metrics can be expanded based on phenomenological and mechanistic considerations of plastic flow instabilities to incorporate anisotropy ratio, rate of work hardening, and dynamic strain aging effects. Separate research thrusts at LaRC are aimed at formulating a protocol for scalable production of fully formed, defect-free components: 1) quantifying the extended set of formability metrics from constant and jump strain rate tensile tests with digital image correlation (DIC); 2) using statistically designed experiments on a vertical spin/flow forming facility (WF VUD-600®) to isolate the effect of convenient process variables, such as axial feed rate and thickness reduction; and 3) predicting the effect of feed ratio (roller translation relative to mandrel rotation speeds) on material response via process modeling using the DEFORM® simulation software. Current results will be discussed in terms of correlating intrinsic material properties with the bulk workability characteristics to successfully fabricate cylinders with complex cross-sections.
Metals and their alloys, except for a few intermetallics, are inherently ductile, i.e. plastic deformation precedes fracture in these materials. Therefore, resistance to fracture is directly related to the development of the plastic zone at the crack tip. Recent studies indicate that the fracture toughness of single crystals depends on the crystallographic orientation of the notch as well as the loading direction. In general, the dependence of crack propagation resistance on crystallographic orientation arises from the anisotropy of (i) elastic constants, (ii) plastic deformation (or slip), and (iii) the weakest fracture planes (e.g. cleavage planes). Because of the triaxial stress state at the notch tips, many slip systems that otherwise would not be activated during uniaxial testing, become operational. The plastic zone formation in single crystals has been tackled theoretically by Rice and his co-workers and only limited experimental work has been conducted in this area. The study of the stresses and strains in the vicinity of a FCC single crystal notch tip is of relatively recent origin. We present experimental and numerical investigation of 3D stress fields and evolution of slip sector boundaries near notches in FCC single crystal tension test specimens, and demonstrate that a 3D linear elastic finite element model that includes the effect of material anisotropy is shown to predict active slip planes and sectors accurately. The slip sector boundaries are shown to have complex curved shapes with several slip systems active simultaneously near the notch. Results are presented for surface and mid-plane of the specimens. The results demonstrate that accounting for 3D elastic anisotropy is very important for accurate prediction of slip activation near FCC single crystal notches loaded in tension. Results from the study will help establish guidelines for fatigue damage near single crystal notches.
Chitin nano- and mesoscale structures present in the exoskeleton of crustaceans exhibit exceptional longitudinal stiffness and toughness, rivaling or even exceeding that of many synthetic polymer architectures. Here, we reveal the origin of the asymmetric shear response in chitin multiscale architectures, marked by pronounced anisotropy in deformation. Under shear aligned with the molecular axis, chitin accommodates strain through coherent atomic rearrangements that enable elastic recovery. In contrast, shear applied perpendicular to this axis induces liquid-like plasticity via localized sliding. These results demonstrate the intrinsic mechanical anisotropy of chitin, underpinning a dual function in resisting repetitive loading while dissipating internal stress during high-strain events. Our findings establish the molecular basis of shear elastoplasticity in multiscale chitin structures, wherein axial elasticity supports energy storage in load-bearing regions, whereas transverse plasticity enables controlled energy dissipation. These atomic-scale insights lay a foundation for the predictive design of chitin-based materials with tunable strength-toughness profiles.
Delamination fracture has limited the use of lightweight Al-Li alloys. In the present study, electron backscattered diffraction (EBSD) methods were used to characterize crack paths in Al-Li alloy C458 (2099). Secondary delamination cracks in fracture toughness samples showed a pronounced tendency for fracture between grain variants of the same deformation texture component. These results were analyzed by EBSD mapping methods and simulated with finite element analyses. Simulation procedures include a description of material anisotropy, local grain orientations, and fracture utilizing crystal plasticity and cohesive zone elements. Taylor factors computed for each grain orientation subjected to normal and shear stresses indicated that grain pairs with the largest Taylor factor differences were adjacent to boundaries that failed by delamination. Examination of matching delamination fracture surface pairs revealed pronounced slip bands in only one of the grains bordering the delamination. These results, along with EBSD studies, plasticity simulations, and Auger electron spectroscopy observations support a hypothesis that delamination fracture occurs due to poor slip accommodation along boundaries between grains with greatly differing plastic response.
Seismic anisotropy of the Earth's mantle has been mostly attributed to crystallographic preferred orientation (CPO) generated during subduction and convection of an anhydrous mantle. But some hydrous phases are also stable at mantle conditions. Here we present results from diamond-anvil cell deformation experiments at high pressure and temperature on hydrous phases Egg [AlSiO 3 (OH)], δ [AlO(OH)] and hydrous stishovite [Si 1-n H 4n O 2 ], transformed from the clay mineral kaolinite. They develop strong CPO during axial compression, suggesting that they likely contribute to seismic anisotropy and heterogeneity in the mantle. Comparing experimental results with viscoplastic polycrystal plasticity models suggest that phase Egg deforms dominantly by (001) slip, δ by (010) slip and stishovite by {100} slip which could be incorporated in future models of mantle geodynamics.