Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “plastic strains”

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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Comparison of three state-of-the-art crystal plasticity based deformation twinning models for magnesium alloys

In magnesium alloys, deformation twinning and its interactions with dislocation slip are responsible for a sigmoidal shape stress–strain behavior and an asymmetrical tension–compression yield strength in magnesium alloys. The sensitivity of twinning to the underlying microstructure renders the crystal plasticity method the most commonly adopted modeling approach for magnesium-twinning. This paper compares three state-of-the-art crystal plasticity-based twinning models from the literature, namely the elastic-viscoplastic self-consistent twinning-detwinning (EVPSC-TDT) model, crystal plasticity finite element model based on enhanced predominate twin reorientation approach (CPFE-ePTR), and the crystal plasticity finite element model based on “discrete twinning” approach (CPFE-DT). A polycrystalline microstructure is simulated with all three methods to compare the resulting stress–strain curves and lattice strains to those from the experimentally measured data. All three methods showed the capability of capturing the experimental results with different levels of accuracy. Additionally, the EVPSC-TDT method avoids solving the finite element matrices and showed the highest computational efficiency. The CPFE-ePTR model shows a higher accuracy in capturing the lattice strain. The CPFE-DT relies on high-resolution finite element mesh and is much slower than the other two methods, but it captured the local deformation concentration and stress reversal phenomena near the twin band, which was not possible with the other two methods. Based on the comparison, guidance for the selection of the appropriate model based on the specific modeling target is provided in this paper.

36 MATERIALS SCIENCE↗
NASA NTRS1974

On numerically accurate finite element

A general criterion for testing a mesh with topologically similar repeat units is given, and the analysis shows that only a few conventional element types and arrangements are, or can be made suitable for computations in the fully plastic range. Further, a new variational principle, which can easily and simply be incorporated into an existing finite element program, is presented. This allows accurate computations to be made even for element designs that would not normally be suitable. Numerical results are given for three plane strain problems, namely pure bending of a beam, a thick-walled tube under pressure, and a deep double edge cracked tensile specimen. The effects of various element designs and of the new variational procedure are illustrated. Elastic-plastic computation at finite strain are discussed.

Nagtegaal, J. C.↗
NASA NTRS1984

Constitutive modeling for isotropic materials (HOST)

The results of the first year of work on a program to validate unified constitutive models for isotropic materials utilized in high temperature regions of gas turbine engines and to demonstrate their usefulness in computing stress-strain-time-temperature histories in complex three-dimensional structural components. The unified theories combine all inelastic strain-rate components in a single term avoiding, for example, treating plasticity and creep as separate response phenomena. An extensive review of existing unified theories is given and numerical methods for integrating these stiff time-temperature-dependent constitutive equations are discussed. Two particular models, those developed by Bodner and Partom and by Walker, were selected for more detailed development and evaluation against experimental tensile, creep and cyclic strain tests on specimens of a cast nickel base alloy, B19000+Hf. Initial results comparing computed and test results for tensile and cyclic straining for temperature from ambient to 982 C and strain rates from 10(exp-7) 10(exp-3) s(exp-1) are given. Some preliminary date correlations are presented also for highly non-proportional biaxial loading which demonstrate an increase in biaxial cyclic hardening rate over uniaxial or proportional loading conditions. Initial work has begun on the implementation of both constitutive models in the MARC finite element computer code.

Lindholm, Ulric S.↗
NASA NTRS1973

Plane elasto-plastic analysis of v-notched plate under bending by boundary integral equation method

A method of solution is presented, which, when applied to the elasto-plastic analysis of plates having a v-notch on one edge and subjected to pure bending, will produce stress and strain fields in much greater detail than presently available. Application of the boundary integral equation method results in two coupled Fredholm-type integral equations, subject to prescribed boundary conditions. These equations are replaced by a system of simultaneous algebraic equations and solved by a successive approximation method employing Prandtl-Reuss incremental plasticity relations. The method is first applied to number of elasto-static problems and the results compared with available solutions. Good agreement is obtained in all cases. The elasto-plastic analysis provides detailed stress and strain distributions for several cases of plates with various notch angles and notch depths. A strain hardening material is assumed and both plane strain and plane stress conditions are considered.

Rzasnicki, W.↗
DOE OSTI2026

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↗
DOE OSTI2022

Predicting extreme anisotropy and shape variations in impact testing of tantalum single crystals

We report recent Taylor cylinder impact tests carried out for Ta single crystals showed strong variations in dimensional changes for different crystallographic directions aligned with the cylindrical axis. In order to capture the effect of crystallography on the deformation characteristics and final shapes of the impacted cylinders, a single crystal material subroutine is adapted and embedded in the solid mechanics/dynamics Finite Element solver Abaqus to simulate the aforementioned single crystal Ta Taylor impact experiments. Details of the coupled model implementation, and insights on the role played by single crystal anisotropic flow on the deformation behavior across a broad range of strain rates and temperatures for different single crystal orientations are presented and discussed. We demonstrate the predictive capability of the adopted crystal plasticity model to capture the significant role played by crystal orientation-induced anisotropy, as well as strain hardening and adiabatic heating, on the dynamic deformation response of crystalline materials. This re-emphasizes the need of microstructure-aware models to improve the accuracy of simulations for high-consequence engineering design.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Crystal plasticity model of BCC metals from large-scale MD simulations

Accurate crystal plasticity models that faithfully capture the behavior of single crystals under a wide range of loading conditions, such as loading direction, strain rate, and temperature, are still lacking. Here we introduce a novel approach in which a crystal plasticity (CP) model is informed directly from and calibrated to large-scale quantum-accurate MD simulations in which single crystal BCC Ta serves as a testbed material. By analyzing our large set of MD simulations several key insights are obtained leading us to modify constitutive assumptions in order to address deficiencies of existing CP models. Importantly, we observe that the standard notion of fixed slip systems – pairs of slip directions and slip planes – is inadequate for describing high-rate plasticity in BCC tantalum at room temperature. Instead, pencil glide defined as dislocation motion in the maximum resolved shear stress planes (MRSSP) of each Burgers vector is fully consistent with our MD simulation data while providing significant simplifications of the constitutive relations. Herein, our resulting new CP model closely matches the behavior of single crystals observed in high-rate MD simulations while being fully consistent with lower rate experimental results.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Report on FY 2021 testing in support of integrating Alloy 800H and 2.25Cr-1Mo for EPP analysis

The use of simplified Elastic-Perfectly Plastic (EPP) design methods avoids the complexities and limitations of the elastic analysis approach in evaluating strain limits and creep-fatigue damage for Class A metallic coolant boundary components in ASME Boiler and Pressure Vessel Code, Section III, Division 5, Subsection HB, Subpart B. No stress classification and linearization are required in the EPP methods. Further, the strain limits and creep-fatigue damage evaluation procedures based on elastic analysis results were developed for the reactor operating conditions where uncoupled, rate-independent plasticity and stationary creep can be used to approximate the deformation in the component adequately. But at very high temperatures, the deformation is rate-dependent, and it can only be approximated adequately by unified viscoplastic model. Thus, for example, the elastic analysis approach cannot be used to evaluate strain limits and creep-fatigue damage for Alloy 617 above 650 °C. The EPP methods remove the temperature restriction. They can be used to evaluate strain limits and creep-fatigue damage for the full range of temperatures.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Plasticity and fracture behavior of Inconel 625 manufactured by laser powder bed fusion: Comparison between as-built and stress relieved conditions

In this work, the influence of stress relief on the plasticity and fracture behavior of Inconel 625 fabricated through laser powder bed fusion additive manufacturing (AM) was investigated. The as-built versus stress relieved microstructures were compared, showing similar grain structures but the presence of ~10 vol % δ phase in the stress relieved condition, and no δ phase in the as-built condition. Mechanical tests under plane strain tension were performed on the stress relieved samples, and an anisotropic plasticity model was calibrated and validated using finite element simulations. Uniaxial and notched tension tests were performed on both as-built and stress relieved samples to probe the effect of stress relief on stress state- and direction-dependent fracture behavior. It was found that on average, the fracture strain of the stress relieved samples along the build direction was 30% higher than that along the perpendicular build direction in the stress state range studied, and the stress relief heat treatment resulted in a 45% decrease in fracture strain. The fracture strain in stress relieved samples was more strongly dependent on stress state than in as-built samples.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A stochastic solver based on the residence time algorithm for crystal plasticity models

Abstract The deformation of crystalline materials by dislocation motion takes place in discrete amounts determined by the Burgers vector. Dislocations may move individually or in bundles, potentially giving rise to intermittent slip. This confers plastic deformation with a certain degree of variability that can be interpreted as being caused by stochastic fluctuations in dislocation behavior. However, crystal plasticity (CP) models are almost always formulated in a continuum sense, assuming that fluctuations average out over large material volumes and/or cancel out due to multi-slip contributions. Nevertheless, plastic fluctuations are known to be important in confined volumes at or below the micron scale, at high temperatures, and under low strain rate/stress deformation conditions. Here, we develop a stochastic solver for CP models based on the residence-time algorithm that naturally captures plastic fluctuations by sampling among the set of active slip systems in the crystal. The method solves the evolution equations of explicit CP formulations, which are recast as stochastic ordinary differential equations and integrated discretely in time. The stochastic CP model is numerically stable by design and naturally breaks the symmetry of plastic slip by sampling among the active plastic shear rates with the correct probability. This can lead to phenomena such as intermittent slip or plastic localization without adding external symmetry-breaking operations to the model. The method is applied to body-centered cubic tungsten single crystals under a variety of temperatures, loading orientations, and imposed strain rates.

42 ENGINEERING↗
DOE OSTI2020

Microstructure, Texture, and Strength Development during High-Pressure Torsion of CrMnFeCoNi High-Entropy Alloy

The equiatomic face-centered cubic high-entropy alloy CrMnFeCoNi was severely deformed at room and liquid nitrogen temperature by high-pressure torsion up to shear strains of about 170. Its microstructure was analyzed by X-ray line profile analysis and transmission electron microscopy and its texture by X-ray microdiffraction. Microhardness measurements, after severe plastic deformation, were done at room temperature. It is shown that at a shear strain of about 20, a steady state grain size of 24 nm, and a dislocation density of the order of 1016 m−2 is reached. The dislocations are mainly screw-type with low dipole character. Mechanical twinning at room temperature is replaced by a martensitic phase transformation at 77 K. The texture developed at room temperature is typical for sheared face-centered cubic nanocrystalline metals, but it is extremely weak and becomes almost random after high-pressure torsion at 77 K. The strength of the nanocrystalline material produced by high-pressure torsion at 77 K is lower than that produced at room temperature. The results are discussed in terms of different mechanisms of deformation, including dislocation generation and propagation, twinning, grain boundary sliding, and phase transformation.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Accounting for the effect of dislocation climb-mediated flow on the anisotropy and texture evolution of Mg alloy, AZ31B

The current study explores the role that dislocation climb has in mediating plasticity a Mg alloy at moderately elevated temperatures. Interrupted tensile tests were performed on samples of Mg alloy, AZ31B, sheet in the O temper condition over a range of strain rates (10 -5 to 10 -1 s -1 ) and temperatures (20–350°C) along the rolling and transverse directions. Experimental measurements of the resulting strain anisotropy and texture evolution were used as constraints during a parametric study employing a new crystal plasticity model (VPSC-CLIMB), which explicitly accounts for the kinematics of dislocation climb. The results reveal that the climb of basal (a) dislocations is not only important for dislocation recovery, but also demonstrate that climb accommodates a significant fraction of the strain in conditions where a power-law creep-type constitutive response prevails. This work does not discredit the notion that non-basal slip of (a) and (c+a) dislocations is important over a wide range of temperatures and strain rates. However, it demonstrates that the activation of dislocation climb as the mechanistic change within the power law regime provides an explanation for a wide range of observations, including the simultaneous reduction in strain anisotropy, slowed texture evolution, and rapid increase in strain rate sensitivity. Finally, it is hypothesized that these conclusions may even apply to cases in which grain boundary sliding and/or dynamic recrystallization are observed.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A better understanding of the mechanics of borehole breakout utilizing a finite strain gradient-enhanced micropolar continuum model

Borehole breakout denotes the failure in rock mass subjected to drilling, caused by stress concentrations exceeding the material strength. Depending on the material properties, the preexisting in situ stress state, and the borehole dimensions, different types of borehole breakout, such as spiral-shaped breakout or v-shaped breakout, are distinguished in the literature. In the present work, we address the influence of the material properties on the predicted borehole breakout mode in a comprehensive finite element study. To this end, herein we employ a gradient-enhanced micropolar damage-plasticity model based on the Mohr–Coulomb strength criterion, formulated in the finite strain regime, which is calibrated based on experimental results from plane strain compression tests on Red Wildmoor sandstone. In the numerical study, the influence of the in situ stress state, the material friction angle, plastic dilation, post peak residual strength, and the inherent material length scale parameters are investigated. Thereby, we demonstrate that depending on the material parameters, substantially different failure modes, characterized by strongly localized shear bands or diffuse failure zones, are predicted. It is shown that in particular the brittleness of the material in the post peak regime has a major influence on the predicted breakout type. Moreover, a statistical validation of the results is obtained by considering different random field distributions of the initial material strength.

58 GEOSCIENCES↗
NASA NTRS1997

Mechanotransduction in bone: osteoblasts are more responsive to fluid forces than mechanical strain

Mechanical force applied to bone produces two localized mechanical signals on the cell: deformation of the extracellular matrix (substrate strain) and extracellular fluid flow. To study the effects of these stimuli on osteoblasts, MC3T3-E1 cells were grown on type I collagen-coated plastic plates and subjected to four-point bending. This technique produces uniform levels of physiological strain and fluid forces on the cells. Each of these parameters can be varied independently. Osteopontin (OPN) mRNA expression was used to assess the anabolic response of MC3T3-E1 cells. When fluid forces were low, neither strain magnitude nor strain rate was correlated with OPN expression. However, higher-magnitude fluid forces significantly increased OPN message levels independently of the strain magnitude or rate. These data indicate that fluid forces, and not mechanical stretch, influence OPN expression in osteoblasts and suggest that fluid forces induced by extracellular fluid flow within the bone matrix may play an important role in bone formation in response to mechanical loading.

Non-NASA Center↗
NASA NTRSDate not supplied

Advances in Integrally Formed Metallic Structures

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.

Aluminum↗
DOE OSTI2020

Creep Deformation and Dynamic Grain Growth in an Interstitial-Free Steel

Dynamic grain growth is demonstrated to be much faster than static grain growth in a body-centered-cubic, interstitial-free steel sheet material at 850°C. Dynamic grain growth occurs during concurrent plastic deformation at elevated temperature, whereas static grain growth occurs during static annealing. Grain growth during steady-state plastic flow in tension at 850°C to a true strain of 0.2 at a true-strain rate of 10 -4 s -1 doubled grain size, while static annealing for the same time produced no increase in grain size. This is described as dynamic normal grain growth (DNGG) because no abnormally large grains were observed. The recrystallized microstructure of the steel demonstrated a log-normal distribution of grain sizes. DNGG produced bimodal grain size distributions that deviate from the theoretical expectation of a simple shift to larger sizes during normal growth. The bimodal distributions contained a remnant of small grains that were not consumed during grain growth. DNGG produced a crystallographic texture that is unique from both the recrystallized material and that produced by lattice rotation alone. DNGG strengthened the {111}⟨110⟩ and {111}⟨112⟩ components of the strong γ-fiber component in the original recrystallization texture. Lattice rotation from tensile deformation, by contrast, strengthened the α-fiber components that intersect the original γ-fiber.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Characterization of a 1D-imaging high-energy x-ray backlighter driven by the National Ignition Facility Advanced Radiographic Capability laser

Plastic deformation of samples compressed to Mbar pressures at high strain rates at the National Ignition Facility (NIF) forms the basis of ongoing material strength experiments in conditions relevant to meteor impacts, geophysics, armor development, and inertial confinement fusion. Hard x-ray radiography is the primary means of measuring the evolution of these samples, typically employing a slit-collimated high-Z microdot driven by the NIF laser to generate >40 keV x rays [E. Gumbrell et al., Rev. Sci. Instrum. 89, 10G118 (2018) and C. M. Huntington et al., Rev. Sci. Instrum. 89, 10G121 (2018)]. Alternatively, a dysprosium “micro-flag” target driven by the Advanced Radiographic Capability laser (~2 kJ, 10 ps) can deliver significantly higher spatiotemporal resolution [M. P. Hill et al., Rev. Sci. Instrum. 92, 033535 (2021)], especially in high-opacity samples. Here, initial experiments revealed problematic brightness and spectral gradients from this source, but by radiographing a set of diamond-turned, 105 µm-thick Pb test objects and supported by simulations using the 3D Monte Carlo code GEANT4, these geometry-dependent gradients across the field of view are quantified and mitigation strategies are assessed. In addition to significantly enhancing the modulation transfer function compared to the existing system, image stacking from multiple layers of image plate is shown to almost double the signal to noise ratio that will reduce uncertainties in future dynamic strength experiments.

47 OTHER INSTRUMENTATION↗