SEARCH · Search NASA

Results for “elasto-viscoplasticity”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Non-local large-strain FFT-based formulation and its application to interface-dominated plasticity of nano-metallic laminates

Here, this paper presents a novel formulation and its robust numerical implementation of strain-gradient (SG) crystal plasticity within a large-strain (LS) elasto-viscoplastic (EVP) fast Fourier transform (FFT)-based micromechanical model. The resulting non-local SG-LS-EVPFFT formulation is used to model and understand the process of kink band formation during layer-parallel compression of nano-metallic laminates (NMLs). NMLs are layered composites with nanoscale thicknesses, thus requiring consideration of the interaction between dislocations and interfaces within the micromechanical model. The length-scale parameter of the SG model is calibrated by simulating a double pile-up and comparing predictions to analytical solution. This required new expressions for the defect energy, resulting in more accurate double pile-up predictions. The calibrated SG-LS-EVPFFT model is then used to simulate layer-parallel compression of copper-niobium NML. Formation of kink bands is predicted, and the model is used to rationalize the microscopic mechanisms enabling the formation process. It is found that accumulation of dislocations at interfaces leads to activation of layer-parallel slip, which in turn leads to kink band formation.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Coupled chemo-mechanical modeling of point-defect diffusion in a crystal plasticity fast Fourier transform framework

Below the yield strength and at moderate-to-high homologous temperatures, the inelastic deformation of metals is mostly governed/rate-controlled by vacancy diffusion-mediated processes. As a function of grain size, stress, temperature and dislocation content, vacancies (or atoms) can adopt preferential diffusion paths across grain interiors, along grain boundaries, or towards and along dislocations, resulting in climb and self-climb. In the steady state and under constant load, grain boundary and grain bulk vacancy diffusion-mediated plasticity have been described in seminal works by Coble and by Nabarro and Herring, respectively. Yet, the interplay between all aforementioned potential diffusion pathways has not been comprehensively mapped. This work presents a thermodynamically-consistent full-field model integrated within a voxel-based elasto-viscoplastic fast Fourier transform framework, which considers the coupling between the diffusion-mediated plasticity mechanisms. In the proposed approach, the kinetics and kinematics of plastic deformation due to vacancy diffusion along grain boundaries and grain bulk, as well as the exchange between grain boundaries and bulk are described explicitly. A homogenization approach at the voxel level is further introduced to simultaneously consider bulk and grain boundary diffusion in a numerically efficient fashion. The new formulation predicts the expected strain rate dependencies and the scaling of the steady-state creep rate with respect to grain size, temperature, and stress. Finally, the model predicts the transition from grain bulk to grain boundary-dominated diffusion with reduction in grain size, a significant step towards capturing transitions in deformation behavior without any phenomenological or ad-hoc adjustments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Generalized grain boundary constitutive description implemented in a strain-gradient large-strain FFT-based formulation: Application to nano-metallic laminates

Herein this paper presents a general treatment of grain boundary constitutive behavior in the context of strain-gradient (SG) plasticity, and its numerical implementation in a large-strain (LS) elasto-viscoplastic (EVP) fast Fourier transform (FFT)-based micromechanical model. Two novel grain boundary constitutive equations are proposed, allowing for more accurate description of the Burgers vector flow at the grain boundary. The capabilities of the generalized SG-LS-EVPFFT formulation are illustrated for the case of kink-band formation during layer-parallel compression of nano-metallic laminates (NMLs), requiring consideration of the interaction between dislocations and interfaces.

36 MATERIALS SCIENCE↗

DOE OSTI2022

New large-strain FFT-based formulation and its application to model strain localization in nano-metallic laminates and other strongly anisotropic crystalline materials

This paper presents a new robust large-strain (LS) elasto-viscoplastic (EVP) formulation based on Fast Fourier Transforms (FFTs) for the prediction of the micro-mechanical response and microstructure evolution of polycrystalline and multiphase materials, with emphasis on the effect of strong crystallographic and/or morphologic anisotropy on localization of plastic deformation. In this work, the novel LS-EVPFFT formulation allows treatment of complex initial geometries and large deformations considering three grids of material points: a regular grid in the reference configuration, where FFTs can be performed; an irregular grid in the initial configuration, created by applying a stress-free displacement field to the reference regular grid; and an irregular grid in the current configuration, undergoing large strains and rotations as the material is loaded. Further numerical stability of the new formulation also required the use of a novel expression for the discrete modified Green’s operator, which reduces spurious field oscillations. After presenting and validating the new formulation by comparison with preexisting implementations and analytical solutions, LS-EVPFFT is applied to the prediction of slip and kink bands formation in polycrystalline columnar ice, and kink bands in single crystal zinc wires, showing good agreement with classic experiments. Finally, the model is used to study kink band formation during compression of Cu–Nb nano-metallic laminates (NMLs), in which accurate treatment of the complex geometry associated with the tortuosity of interfaces and large deformations become critical, showing consistency with corresponding micropillar experiments.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Twinning pathways enabled by precipitates in $\mathrm{AZ91}$

While precipitates have been shown to substantially strengthen magnesium alloys by blocking the glide of dislocations inside the grain, the interactions between these precipitates and deformation twins, which commonly occur in these alloys, are much less understood. Here, in this work, an elasto-viscoplastic fast-Fourier-transform (EVP-FFT) model is used to study the interactions between plate-shaped basal precipitates and propagating $10\bar{1}2$ tensile twins in AZ91 Mg alloy. The results suggest that while precipitates may impede the propagation and thickening of twins, they can also cause stress localizations that can promote the formation of multiple new twins of the same or different crystallographic variants. We show that the location of the twin-precipitate interaction site, whether precipitate-central or precipitate-edge impingement, and the thickness of the precipitate can influence the propensity for twins to expand around the precipitate or nucleate a new twin on the other side of it. Depending on the twin-precipitate impingement site, we propose multiple twinning pathways that can help explain how twins can proliferate in the magnesium alloys in the presence of precipitates.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Yield strength of CeO 2 measured from static compression in a radial diamond anvil cell

Cerium oxide (ceria, CeO 2 ) is frequently used as a standard in applications such as synchrotron and x-ray free electron lasers for calibrating x-ray wavelengths and offers the potential for understanding the high pressure properties and deformation mechanisms in a wide range of similar face centered cubic (fcc) materials. In this study, the pressure dependence of the strength of ceria was investigated up to 38 GPa using angle dispersive x-ray diffraction in a radial geometry in a diamond anvil cell. In this experiment, the difference in the stress along the axis of compression and perpendicular to the direction of compression can be determined, giving a quantity known as the differential stress. It was found that the differential stress (t), a measure of the lower bound for yield strength, initially increases rapidly from 0.35 ± 0.06 GPa to 2.2 ± 0.4 GPa at pressures of 1.8 and 3.8 GPa, respectively. Above 4 GPa, t increases more slowly to 13.8 ± 2.6 GPa at a pressure of 38 GPa. The changes in the preferred orientation (texture) of CeO 2 with pressure were also measured, allowing for the determination of active deformation mechanisms using an elasto-viscoplastic self-consistent model (EVPSC). It was found that as pressure increased, the [001] direction had a slight preferred orientation along the axis of compression. Our EVPSC model of experimental fiber (cylindrically symmetric) textures and lattice strains were most consistent with dominant slip activity along ${111}$ $\langle$$1\bar10$$\rangle$.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Resolving crystallographic geometrically necessary dislocations in three dimensions in a hexagonal close packed titanium alloy

Abstract Geometrically necessary dislocation (GND) content is measured from mm 3 -scaled Ti7Al three-dimensional (3D) microstructural data using a theory extended for hexagonal close packed crystals, which accounts for basal, prismatic and pyramidal ⟨ c + a ⟩ type dislocation content. The Ti7Al samples have been mechanically pre-strained to two different strain levels, and will then be strained along the same axis in uniaxial tension during simulation. Both inter- and intragranular GNDs across the microstructures have been characterized, with a large contribution of pyramidal ⟨ c + a ⟩ GNDs, consistent with the relative slip activity involved in pre-straining. The spatially resolved crystallographic GND distributions within the 3D microstructures are used to instantiate a microstructure model for forward modeling deformation simulations by a dislocation density hardening elasto-viscoplastic fast Fourier transform framework. Coarsening the voxel resolution during the initial microstructure construction procedure is shown to strongly impact both the magnitude and spatial distribution of the GNDs and in turn the forward deformation response of the pre-strained material. This study indicates that the voxel resolution desired when transferring from measured to model microstructures need not only be proportionally scaled with the microstructure but also sufficiently fine to capture the subgranular orientation gradients that may already be present in the material.

Materials Science↗

DOE OSTI2020

Particle interspacing effects on the mechanical behavior of a Fe–TiB 2 metal matrix composite using FFT-based mesoscopic field dislocation mechanics

This paper presents an application to metal matrix composites (MMCs) of an enhanced elasto-viscoplastic Fast Fourier Transform (EVP-FFT) formulation coupled with a phenomenological continuum Mesoscale Field Dislocation Mechanics (MFDM) theory. Contrary to conventional crystal plasticity, which only accounts for plastic flow and hardening induced by statistically stored dislocations (SSDs), MFDM-EVP-FFT also describes the evolution of polarized geometrically necessary dislocation (GND) density and its effect on both plastic flow and hardening. Numerical results for a Fe–TiB 2 MMC made of a ferrite matrix (α-Fe) and elastic ceramic particles (TiB 2 ) are presented. Full-field simulations are performed using synthetic periodic unit cells representative of the MMC, with single-crystalline and polycrystalline matrix, for different particle interspacing distances. A strong dependence of the predicted equivalent stress, cumulated plastic strain and GND density fields with particle interspacing distance is observed, in contrast with conventional crystal plasticity. Correlations between these mechanical fields and microstructural features, and their influence on local and global mechanical behavior are examined for the different MMC microstructures.

42 ENGINEERING↗

DOE OSTI2021

Deformation and strength of mantle relevant garnets: Implications for the subduction of basaltic-rich crust

Garnet is an important mineral phase in the upper mantle as it is both a key component in bulk mantle rocks, and a primary phase at high pressure within subducted basalt. In this work, we focus on the strength of garnet and the texture that develops within garnet during accommodation of differential deformational strain. We use X-ray diffraction in a radial geometry to analyze texture development in situ in three garnet compositions under pressure at 300 K: a natural garnet (Prp 60 Alm 37 ) to 30 GPa, and two synthetic majorite-bearing compositions (Prp 59 Maj 41 and Prp 42 Maj 58 ) to 44 GPa. All three garnets develop a modest (100) texture at elevated pressure under axial compression. Elasto-viscoplastic self-consistent (EVPSC) modeling suggests that two slip systems are active in the three garnet compositions at all pressures studied: {110}<$1\bar{1}1$> and {001}<110>. We determine a flow strength of ~5 GPa at pressures between 10 to 15 GPa for all three garnets; these values are higher than previously reported yield strengths measured on natural and majoritic garnets. Strengths calculated using the experimental lattice strain differ from the strength generated from those calculated using EVPSC. Prp 67 Alm 33 , Prp 59 Maj 41 , and Prp 42 Maj 58 are of comparable strength to each other at room temperature, which indicates that majorite substitution does not greatly affect the strength of garnets. Additionally, all three garnets are of similar strength as lower mantle phases such as bridgmanite and ferropericlase, suggesting that garnet may not be notably stronger than the surrounding lower mantle/deep upper mantle phases at the base of the upper mantle.

58 GEOSCIENCES↗

DOE OSTI2024

Experimental and Numerical Characterization of High Temperature Deformation Behavior of 347H Stainless Steel

This study investigates how temperature affects the plasticity and thermal creep behavior of 347H stainless steel under uniaxial tension. The research combined experimental testing with advanced computational modeling. Two types of experiments were conducted: uniaxial tensile tests at temperatures from 100°C to 750°C using strain rates of ~10⁻⁴ s⁻¹, and creep tests at temperatures between 600°C and 750°C under various stress levels. These experimental results were used to develop and validate a new integrated mechanistic model that can predict material behavior under any loading condition while accounting for both stress and temperature effects. The model was implemented using a polycrystalline microstructure simulation framework based on elasto-viscoplastic Fast Fourier Transform (EVPFFT). It incorporates three key deformation mechanisms: thermally activated dislocation glide, dislocation climb, and vacancy diffusional creep. The model accounts for internal stress distribution within single crystals and considers how precipitates and solute atoms (both interstitial and substitutional) affect dislocation movement. After validation against experimental data, the model was used to generate Ashby-Weertman deformation mechanism maps for 347H steel, providing new insights into how microstructure influences the activation of different creep mechanisms.

Joy, Jobin K.↗

DOE OSTI2024

Deformation and Transformation Textures in the NaMgF3 Neighborite—Post-Perovskite System

The D″ region of the lower mantle, which lies just above the core–mantle boundary, is distinct from the bulk of the lower mantle in that it exhibits complex seismic heterogeneity and seismic anisotropy. Seismic anisotropy in this region is likely to be largely due to the deformation-induced texture (crystallographic preferred orientation) development of the constituent mineral phases. Thus, seismic anisotropy can provide a marker for deformation processes occurring in this dynamic region of the Earth. Post-perovskite-structured (Mg,Fe)SiO3 is believed to be the dominant mineral phase in many regions of the D”. As such, understanding deformation mechanisms and texture development in post-perovskite is important for the interpretation of observed seismic anisotropy. Here, we report on high-pressure diamond anvil cell deformation experiments on NaMgF3 neighborite (perovskite structure) and post-perovskite. During deformation, neighborite develops a 100 texture, as has been previously observed, both in NaMgF3 and MgSiO3 perovskite. Upon transformation to the post-perovskite phase, an initial texture of {130} at high angles to compression is observed, indicating that the {100} planes of perovskite become the ~{130} planes of post-perovskite. Further compression results in the development of a shoulder towards (001) in the inverse pole figure. Plasticity modeling using the elasto-viscoplastic self-consistent code shows this texture evolution to be most consistent with deformation on (001)[100] with some contribution of glide on (100)[010] and (001)<110> in NaMgF3 post-perovskite. The transformation and deformation mechanisms observed in this study in the NaMgF3 system are consistent with the behavior generally observed in other perovskite–post-perovskite systems, including the MgSiO3 system. This shows that NaMgF3 is a good analog for the mantle bridgmanite and MgSiO3 post-perovskite.

Geochemistry & Geophysics↗

NASA NTRS1982

Structural analysis

Hot section components of aircraft gas turbine engines are subjected to severe thermal-structural loading conditions, especially during the start-up and take-off portions of the engine cycle. The most severe and damaging stresses and strains are those induced by the steep thermal gradients induced during the start-up transient. These transient stresses and strains are also the most difficult to predict, in part because the temperature gradients and distributions are not well known or predictable, and also because the cyclic elasto-viscoplastic behavior of the materials at these extremes of temperature and strain are not well known or predictable. One element of the structures program will develop improved time-varying thermal-mechanical load models for the entire engine mission cycle from start-up to shutdown. The thermal model refinements will be consistent with those required by the structural code including considerations of mesh-point density, strain concentrations, and thermal gradients. Models will be developed for the burner liner, turbine vane and turbine blade.

Johns, R. H.↗

NASA NTRS1987

Thermo-elasto-viscoplastic analysis of problems in extension and shear

The problems of extension and shear behavior of structural elements made of carbon steel and subjected to large thermomechanical loads are investigated. The analysis is based on nonlinear geometric and constitutive relations, and is expressed in a rate form. The material constitutive equations are capable of reproducing all nonisothermal, elasto-viscoplastic characteristics. The results of the test problems show that: (1) the formulation can accommodate very large strains and rotations; (2) the model incorporates the simplification associated with rate-insensitive elastic response without losing the ability to model a rate-temperature dependent yield strength and plasticity; and (3) the formulation does not display oscillatory behavior in the stresses for the simple shear problem.

Riff, R.↗

NASA NTRS1987

The dynamic aspects of thermo-elasto-viscoplastic snap-through and creep buckling phenomena

Use of a mathematical model and solution methodology, to examine dynamic buckling and dynamic postbuckling behavior of shallow arches and spherical caps made of a realistic material and undergoing non-isothermal, elasto-viscoplastic deformation was examined. Thus, geometric as well as material type nonlinearities of higher order are included in this analysis. The dynamic stability problem is studied under impulsive loading and suddenly applied loading with loads of constant magnitude and infinite duration. A finite element model was derived directly from the incrementally formulated nonlinear shell equations, by using a tensor-oriented procedure. As an example of the results, the time history of the midspan displacement of a damped shallow circular arch is presented.

Riff, R.↗

NASA NTRS1987

Non-isothermal elastoviscoplastic snap-through and creep buckling of shallow arches

The problem of buckling of shallow arches under transient thermomechanical loads is investigated. The analysis is based on nonlinear geometric and constitutive relations, and is expressed in a rate form. The material constitutive equations are capable of reproducing all non-isothermal, elasto-viscoplastic characteristics. The solution scheme is capable of predicting response which includes pre and postbuckling with creep and plastic effects. The solution procedure is demonstrated through several examples which include both creep and snap-through behavior.

Simitses, G. J.↗

NASA NTRS1988

Thermo-elasto-viscoplastic analysis of problems in extension and shear

Riff, R.↗

NASA NTRS1988

Non-isothermal buckling behavior of viscoplastic shell structures

Described are the mathematical model and solution methodologies for analyzing the structural response of thin, metallic elasto-viscoplastic shell structures under large thermomechanical loads and their non-isothermal buckling behavior. Among the system responses associated with these loads and conditions are snap-through, buckling, thermal buckling, and creep buckling. This geometric and material nonlinearities (of high order) can be anticipated and are considered in the model and the numerical treatment.

Riff, Richard↗

NASA NTRS1995

Adaptive Finite Element Methods for Continuum Damage Modeling

The paper presents an application of adaptive finite element methods to the modeling of low-cycle continuum damage and life prediction of high-temperature components. The major objective is to provide automated and accurate modeling of damaged zones through adaptive mesh refinement and adaptive time-stepping methods. The damage modeling methodology is implemented in an usual way by embedding damage evolution in the transient nonlinear solution of elasto-viscoplastic deformation problems. This nonlinear boundary-value problem is discretized by adaptive finite element methods. The automated h-adaptive mesh refinements are driven by error indicators, based on selected principal variables in the problem (stresses, non-elastic strains, damage, etc.). In the time domain, adaptive time-stepping is used, combined with a predictor-corrector time marching algorithm. The time selection is controlled by required time accuracy. In order to take into account strong temperature dependency of material parameters, the nonlinear structural solution a coupled with thermal analyses (one-way coupling). Several test examples illustrate the importance and benefits of adaptive mesh refinements in accurate prediction of damage levels and failure time.

Min, J. B.↗