Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “elasto-plasticity”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Effect of heat treatment on the tensile behavior of selective laser melted Ti-6Al-4V by in situ X-ray characterization

Selective Laser Melted Ti-6Al-4V (as-SLMed) exhibits decreased yield strength, increased work hardening, and increased ductility after heat treatment at 730 degrees C (HT-730) or 900 degrees C (HT-900) for 2 h. To understand the change of mechanical properties, in situ high energy X-ray diffraction (HEXRD) is used to examine the phase composition, load partitioning, slip system activity, and dislocation density evolution in all three specimens. The as-SLMed specimen is dominated by martensitic α'. After heat treatment, α' partly or fully decomposes into α+β, reducing the yield strength. In HT-730, beta precipitates with confined size show much higher lattice strain than the α'/α matrix during deformation; in HT-900, the lattice strain difference is mostly eliminated. This is a key reason for the increased ductility in HT-900. From the anisotropic lattice strain development, basal slip is identified as the easiest slip system in α'/α. Using an elasto-plastic self-consistent (EPSC) model, the critical resolved shear stress ratio between prismatic slip and basal slip (CRSS prism /CRSS basal ) is estimated to be 1.31 and 1.16 in the as -SLMed and the HT-900 specimens, respectively. Further, the alpha phase in HT-900 is able to activate multiple slip systems and accumulate more dislocations during plastic deformation. This explains why HT-900 has better ductility and higher work hardening rate than the other two specimens.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Interface morphology and dislocation-mediated processes during rapid solidification of thin films

Rapid solidification experiments have, in recent years, revealed a wealth of new microstructural phenomena that suggest a strong connection between the kinetics of solidification and the crystalline structures that emerge as a result. In this work, we investigate the interplay between interface morphology and defect-mediated processes during rapid solidification conditions using a Phase Field Crystal (PFC) model, enabling us to simultaneously and efficiently explore the physics of solidification and elasto-plasticity in the formalism of a single-field theory. We predict that there are two mechanisms by which dislocations emitted directly from the solid–liquid interface induce orientation gradients as well as the formation of subgrain boundaries within a single solidifying cell. We relate these mechanisms to the morphology of the moving solid–liquid interface and identify a suitable control parameter in the PFC model with which we can go between said morphologies by effectively changing the relative strength of the capillary length and kinetic coefficients of the solid–liquid interface. Thus, we are able to provide mechanistic explanations for several microstructural features (with an emphasis on orientation gradients and subgrain boundaries) observed during the rapid solidification of pure materials. We also provide a simple explanation for the formation of “jagged” subgrain boundaries, which is consistent with our experimental observations in rapidly solidified samples of Aluminum, whose mechanisms have thus far been unknown.

Interface morphology↗
DOE OSTI2021

Identification of crystal plasticity model parameters by multi-objective optimization integrating microstructural evolution and mechanical data

Crystal plasticity models evolve a polycrystalline yield surface using meso-scale descriptions of deformation mechanisms. The activation of deformation mechanisms is governed by crystallography and a set of model parameters, which are typically calibrated through the fitting of mechanical data such as stress–strain curves and elastic lattice strains. Microstructural data such as phase fractions and texture evolution are used for verifying crystal plasticity parameters. In this study, we use a multi-objective genetic algorithm to identify hardening parameters from flow stress curves with an option to incorporate texture into the optimization approach. Robust, generalized objective functions are developed and used to identify sets of parameters pertaining to dislocation density-based hardening laws in visco-plastic and elasto-plastic self-consistent (VPSC and EPSC) homogenization models. First, the parameters are identified for pure Nb directly from texture using an objective function based on generalized spherical harmonics. Since texture evolution is driven by the relative contribution of active slip systems, the parameters governing the evolution of slip resistance ratios can be recovered from fitting discrete textures at a series of strains. Next, a comprehensive set of load reversal data for dual phase (DP) 780 steel is used to fit a hardening law and a back-stress law in EPSC. Finally, parameters pertaining to a complex hardening law for the evolution of slip and twinning in pure α-Ti are identified. Remarkably, using texture as an objective in combination with stress–strain objectives constrains the model of Ti to fully reproduce not only stress–strain and texture evolution but also hierarchical twinning measurements as a function of initial grain size and texture. Furthermore, given an appropriate model fit to representative experimental texture evolution, underlying twin volume fractions contributing to texture evolution can be predicted.

42 ENGINEERING↗
DOE OSTI2024

Poromechanical cohesive interface element with combined Mode I-II cohesive zone elastoplasticity for simulating fracture in fluid-saturated porous media

A combined Mode I-II cohesive zone (CZ) elasto-plastic constitutive model, and a two-dimensional (2D) cohesive interface element (CIE) are formulated and implemented at small strain within an ABAQUS User Element (UEL) for simulating 2D crack nucleation and propagation in fluid-saturated porous media. Here, the CZ model mitigates problems of convergence for the global Newton-Raphson solver within ABAQUS, which when combined with a viscous stabilization procedure allows for simulation of post-peak response under load control for coupled poromechanical finite element analysis, such as concrete gravity dam stability analysis. Verification examples are presented, along with a more complex ambient limestone-concrete wedge fracture experiment, water-pressurized concrete wedge experiment, and concrete gravity dam stability analyses. A calibration procedure for estimating the CZ parameters is demonstrated with the limestone-concrete wedge fracture process. For the water-pressurized concrete wedge fracture experiment it is shown that the inherent time-dependence of the poromechanical CIE analysis provides a good match with experimental force versus displacement results at various crack mouth opening rates, yet misses the pore water pressure evolution ahead of the crack tip propagation. This is likely a result of the concrete being partially-saturated in the experiment, whereas the finite element analysis assumes fully water saturated concrete. For the concrete gravity dam analysis, it is shown that base crack opening and associated water uplift pressure leads to a reduced Factor of Safety, which is confirmed by separate analytical calculations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A mean field homogenization model for the mechanical response of ceramic matrix composites

SiC/SiC composites offer exceptional mechanical stability at high temperatures and under irradiation. These ceramic matrix composites are therefore strong candidate materials for future nuclear energy applications. Their mechanical response, which exhibits pseudo-plasticity, is mediated by matrix cracking, fiber debonding, and fiber pull-out due to slip. Here, this study introduces a mechanistic model for the behavior of unidirectionally reinforced SiC/SiC composites. Specifically a mean field homogenization approach is proposed to account for all deformation and degradation modes during mechanical deformation. The homogenization scheme relies on a Mori Tanaka method that is extended to consider the effects of the coating’s elasto-plastic response on the development of micromechanical fields. Further, the model proposed introduces a method to effectively account for the role of localized damage (i.e., cracks) on mechanical fields within both the fiber and the matrix. Upon validating the model against experimental data, the roles of interface sliding, coating dimensions and intrinsic elastic response, as well as of microstructure (e.g. porosity, fiber volume fraction) are discussed.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Thermal shock resistance of lightweight cements developed for geothermal conditions

Cements are a critical component in well construction, as they act to prevent well fluid and gas escape, prevent corrosion of the casing, and strengthen the wellbore to prevent deformation. Under the high temperature/pressure conditions common in geothermal systems, the injection of cold water for energy production is expected to induce cyclic damage to the borehole cement through the rapid temperature fluctuations. These “thermal shocks” are expected to cause casing shrinkage, annulus formation, and cement tensile stresses. To understand the effect of cold water injection on the wellbore environment, a set of rock-cement-steel samples were created to simulate the structure of a geothermal well. Lightweight thermally-insulating cement blends were tested under thermal shock conditions in this study. In each test, the samples were pressurized to an effective pressure of ~3.5 MPa and placed at high temperatures. Thermal shocks were performed by injecting cold water (~10-15 °C) through the samples at a constant rate while keeping the samples at high temperatures until the sample temperature stopped decreasing and deformation ceased. Eight thermal shock tests were conducted with each sample – two at 100 °C and six at 200 °C. Post-tests analysis was then conducted by cutting open each sample to examine the damage in each component of the simulated wellbore. Experimental results suggest that all samples experienced similar degrees of axial and lateral contraction during cold water injection, but for the most part this contraction is recoverable when injection halts. Post-test analysis revealed that fly ash cenosphere pre-treatment had the best effect on improving thermal shock resistance in the cement blends. Thermomechanical modeling of likely stress paths experienced by the cements during heating/cooling cycles shows that elasto-plastic cement constitutive behavior results in most plastic strain occurring during the initial heating steps, with mostly elastic strain occurring during the thermal shock cycles. In conclusion, this agrees with experimental evidence, suggesting that cement damage from shocking occurs via other mechanisms such as chemical alteration, corrosion, and fatigue.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Influence of cementation on the yield surface of rocks numerically determined from digital microstructures

Digital Rock Physics has reached a level of maturity on the characterisation of primary properties that depend on the microstructure – such as porosity, permeability or elastic moduli – by numerically solving field equations on μCT scan images of rock. After small deformations or at depth though, most rocks eventually reach their limit of elasticity and the complementary plastic properties are needed to describe the full mechanical behaviour. Currently, determination of a rock’s yield surface from its microstructure is often restricted to semi-analytical criteria derived by limit analysis or numerical simulations performed on idealised geometries. Such simplification lacks representativeness, particularly for processes that affect directly the pore-grain interface such as the cementation phenomenon, happening during diagenesis. Eventually, only direct numerical simulation of elasto-plasticity performed on digitalised microstructures can be used to assess the strength of different cemented materials and its evolution with the alteration of the microstructure. In this study, we provide a comprehensive parametric study on the impact of cementation on rock strength for real microstructures of cemented granular materials. Compared to most previous studies, the whole yield surface is determined numerically (using Finite Element Method) in order to assess the influence of cementation for different stress-paths. The previously known tendency of rock to strengthen with increasing cementation volume is verified. New results on the influence of cement property namely Young’s modulus, friction and cohesion on the rock’s yield surface are explored. The envelopes obtained are compared to the ones obtained by experimental data and existing models. Finally, the framework presented in this study showcases the wider possibility of determining any rock’s or porous material’s yield surface from its microstructure.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Crystal plasticity modeling of strain-induced martensitic transformations to predict strain rate and temperature sensitive behavior of 304 L steels: Applications to tension, compression, torsion, and impact

This paper advances crystallographically-based Olson-Cohen (direct γ → α’) and deformation mechanism (indirect γ→ε→α’) phase transformation models for predicting strain-induced austenite to martensite transformation. Here, the advanced transformation models enable predictions of not only strain-path sensitive, but also of strain-rate and temperature sensitive deformation of polycrystalline stainless steels (SSs). The deformation of constituent grains in SSs is modeled as a combination of anisotropic elasticity, crystallographic slip, and phase transformation, while the hardening is based on the evolution of dislocation density and explicit shifts in phase fractions. Such grain-scale deformation is implemented within the meso-scale elasto-plastic self-consistent (EPSC) homogenization model, which is coupled with the implicit finite element (FE) method to provide a constitutive response at each FE integration point for solving boundary value problems at the macro-scale. Parameters pertaining to the hardening and transformation models within FEEPSC are calibrated and validated on a suite of data including flow curves and phase fractions for monotonic compression, tension, and torsion as a function of strain-rate and temperature for wrought and additively manufactured (AM) SS304L. To illustrate the potential and accuracy of the integrated multi-level FE-EPSC simulation framework, geometry, mechanical response, phase fractions, and texture evolution are simulated during gas-gun impact deformation of a cylinder and quasi-static tension of a notched specimen made of AM SS304L. Details of the simulation framework, comparison between experimental and simulation results, and insights from the results are presented and discussed.

304L steels↗
DOE OSTI2021

Learning constitutive relations using symmetric positive definite neural networks

In this work, we present a new neural-network architecture, called the Cholesky-factored symmetric positive definite neural network (SPD-NN), for modeling constitutive relations in computational mechanics. Instead of directly predicting the stress of the material, the SPD-NN trains a neural network to predict the Cholesky factor of the tangent stiffness matrix, based on which the stress is calculated in incremental form. As a result of this special structure, SPD-NN weakly imposes convexity on the strain energy function, satisfies the second order work criterion (Hill's criterion) and time consistency for path-dependent materials, and therefore improves numerical stability, especially when the SPD-NN is used in finite element simulations. Depending on the types of available data, we propose two training methods, namely direct training for strain and stress pairs and indirect training for loads and displacement pairs. We demonstrate the effectiveness of SPD-NN on hyperelastic, elasto-plastic, and multiscale fiber-reinforced plate problems from solid mechanics. The generality and robustness of SPD-NN make it a promising tool for a wide range of constitutive modeling applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

A phase-field fracture formulation for generalized standard materials: The interplay between thermomechanics and damage

Accurately modeling fracture of ductile materials poses open challenges in the field of computational mechanics due to the multiphysics nature of their failure processes. Integrating the interplay between thermodynamics and damage into ductile fracture models is vital for predicting critical failure modes. Here, in this paper, we develop a versatile phase-field (PF) framework for modeling ductile fracture, taking into account finite-strain elasto-plasticity. The framework stems from a variational formulation of constitutive relations for generalized standard materials (GSMs), whose response is described by a Helmholtz free energy and a dissipation pseudo-potential. Its variational structure is based on a minimum principle for a functional that expresses the sum of power densities for reversible and irreversible processes. By minimizing this functional with a constraint on a von Mises yield function, we derive the evolution equation for the equivalent plastic strain and an associative flow rule. This constrained optimization problem is analytically solved for a wide class of thermo-viscoplasticity models. The key innovations of the current work include (i) a cubic plastic degradation function that accounts for a non-vanishing damage-dependent yield stress, (ii) closed-form expressions of the Helmholtz free energy and dissipation pseudo-potential for three thermo-viscoplasticity models, (iii) an extended Johnson–Cook plasticity model with a nonlinear hardening law, and (iv) a plastic work heat source that depends on the plastic degradation function and a variable Taylor–Quinney (TQ) coefficient. The capabilities of the proposed framework are tested with the aid of four ductile fracture problems, including the Sandia Fracture Challenge. In each of these problems, we examine the evolution of relevant field variables such as the PF order parameter, the equivalent plastic strain, the temperature, and the internal power dissipation density, in addition to the overall structural response quantified by the force–displacement curve. These numerical studies demonstrate that the proposed framework effectively represents ductile fracture, yielding computational results that exhibit good agreement with experimental data.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Origins of high ductility exhibited by an extruded magnesium alloy Mg-1.8Zn-0.2Ca: Experiments and crystal plasticity modeling

Low ductility and strength are major bottlenecks against Mg alloys' wide applications. In this work, we systematically design the composition and fabrication process for a low-alloyed Mg-Zn-Ca alloy, showing that it can be extruded at low temperatures (~ 250 degrees C) and high speeds (~2 mm/s). After the extrusion, this alloy exhibits a substantially weakened basal texture, relatively small grain size, very high tensile elongation (~ 30%), and good strength. The origin of the considerably improved ductility was studied using a combination of three-dimensional atom probe tomography (3D-APT), transmission electron microscopy (TEM), electron backscattered diffraction (EBSD) in conjunction with surface slip trace analysis, in-situ synchrotron X-ray diffraction, and elasto-plastic self-consistent (EPSC) modeling. Co-segregation of Zn and Ca atoms at a grain boundary is observed and associated with texture weakening and grain boundary mediated plasticity, both improving the ductility. While basal slip and prismatic slip are identified as the dominant deformation systems in the alloy, the ratio between their slip resistances is substantially reduced relative to pure Mg and most other Mg alloys, significantly contributing to the improved ductility of the alloy. This Mg-Zn-Ca alloy exhibiting excellent mechanical properties and low fabrication cost is a promising candidate for industrial productions.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Size effect in stainless steel thin wires under tension

Miniaturization of modern devices drives the exploration into how materials behave in small dimensions. However, the mechanical behaviors of micro-scale specimens under varied loading conditions are still challenging to characterize and predict. In this work, specimen size-dependent weakening effect was experimentally revealed in 316L stainless steel thin wires under tension with wire diameter decreasing from 200 μm down to 61 μm. The underlying mechanisms were primarily associated with surface grain softening and strain localization, which were largely dependent on the number of grains across wire diameter. Through an analog composite model with surface and interior grains that follows the Hall-Petch relation, the extent of surface grain softening was found to coincide with the flow stress reduction due to lack of boundary strengthening. Further using elasto-plastic self-consistent polycrystal modeling, the stress-strain behaviors of interior and surface grains were well described in analogy to those of bulk materials with similar-size grains and extremely large grains, respectively. Subsequently, the analog composite model well simulated the specimen size-dependent weakening effect under tension. This article provided more insights into the size effect exploration and proposed a practical model for describing size-dependent mechanical behaviors of small-scale materials to be applied to miniaturized devices.

316L stainless steel↗
DOE OSTI2025

Absorbing boundary conditions in material point method adopting perfectly matched layer theory

This study focuses on solving the numerical challenges of imposing absorbing boundary conditions for dynamic simulations in the material point method (MPM). To attenuate elastic waves leaving the computational domain, the current work integrates the Perfectly Matched Layer (PML) theory into the implicit MPM framework. The proposed approach introduces absorbing particles surrounding the computational domain that efficiently absorb outgoing waves and reduce reflections, allowing for accurate modeling of wave propagation and its further impact on geotechnical slope stability analysis. The study also includes several benchmark tests to validate the effectiveness of the proposed method, such as several types of impulse loading and symmetric and asymmetric base shaking. The conducted numerical tests also demonstrate the ability to handle large deformation problems, including the failure of elasto-plastic soils under gravity and dynamic excitations. The findings extend the capability of MPM in simulating continuous analysis of earthquake-induced landslides, from shaking to failure.

58 GEOSCIENCES↗
DOE OSTI2020

Fracture Failure Analysis for MQXFA Magnet Aluminum Shells

The High-Luminosity LHC Accelerator Upgrade Project (HL-LHC AUP) is approaching the production phase of the US-contributed Q1 and Q3 Interaction Region Quadrupoles (MQXFA). The structures for the MQXFA prototypes were design and inspected by the US-LARP (LHC Accelerator Research Program), AUP developed criteria, which will be used for the pre-series structures. As the first two full-length prototypes with 4.2 m magnetic length, MQXFAP1 and MQXFAP2, were designed and assembled at Lawrence Berkeley National Laboratory (LBNL), and tested at Brookhaven National Laboratory (BNL). The end aluminum short shell of MQXFAP2 was fractured along the shell length during the test, and tests were stopped. Analytical and Finite Element analysis were performed in light of the graded procedure defined in the Structure Design Criteria to investigate the fracture failure for MQXFAP2. In this paper, we report the fracture analysis of the current shell design, including the elasto-plastic simulations with submodel technique, and calculations with Linear Elastic Fracture Mechanics (LEFM). Test material properties are also presented. The results of this analysis explain why the end shell of MQXFAP2 failed, and suggest fillets on the end shell notches to meet the margin specified in the Structural Design Criteria.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

A crystal plasticity finite element method modeling of Zircaloy with hydride phases based on scanning electron micrographs

This work developed a microstructure-based finite element model to predict the stress state of alloys with second phase inclusions. Quantitative microstructural details extracted from scanning electron microscopy (SEM) images were used to generate heterogeneous microstructures including the morphology and spatial distribution of hydrides. Generation of digital microstructure was achieved through two steps of tessellations using software Neper and Matlab. The process is demonstrated using an example of Zircaloy material with secondary phases of hydrides dispersed within and stress-strain response of Zircaloy containing hydrides was predicted. The constitutive material model for Zircaloy in this study was based on crystal plasticity theory which considers the hexagonal close-packed (HCP) atomic structure of Zircaloy material. The hydrides were modeled as isotropic elasto-plastic material. A parametric study had been conducted to understand the effect of volume fraction, and lamellae thickness of the hydride phase on the mechanical properties of the overall material. Results can help designers to alter the manufacturing process to obtain the enhanced mechanical properties for components used in nuclear applications made by Zircaloy material.

Kulkarni, Shank S.↗
DOE OSTI2024

Tardigrade, Version 1.1x

Tardigrade is a collection of functions, libraries, and stand-alone code which enables the simulation of materials using the higher-order micromorphic framework. The tools present are the interface to INL's MOOSE code which is used as the underlying FEA solver, several constitutive equations (linear elasticity, a pressure-sensitive elasto-plastic model, crystal plasticity, and additional models for geo-materials) useful for demonstration purposes and the simulation of some materials, the micromorphic filter which converts DNS information to micromorphic stress-deformation quantities, and the overlap coupling framework which enables multi-scale simulations of a coupled higher-order macroscale and a classically defined microscale.

Miller, Nathan↗
DOE OSTI2025

Green's function fast multipole method for continuum mechanics (SM-FMM)

Solid Mechanics Fast Multipole Method based on elastic Green's function and accelerated using FFTs The code calculates the mechanical fields (stress and strain) for a heterogeneous elasto-plastic material unit cell under quasi-static conditions (no dynamic effects). The fields are calculated using a Green's function method, where the strain field is given by a discrete convolution of Green's operator with an auxiliary stress field. The convolution is calculated using the fast multipole method.

Zecevic, Miroslav↗
DOE OSTI2021

A novel physics-based and data-supported microstructure model for part-scale simulation of laser powder bed fusion of Ti-6Al-4V

The elasto-plastic material behavior, material strength and failure modes of metals fabricated by additive manufacturing technologies are significantly determined by the underlying process-specific microstructure evolution. In this work a novel physics-based and data-supported phenomenological microstructure model for Ti-6Al-4V is proposed that is suitable for the part-scale simulation of laser powder bed fusion processes. The model predicts spatially homogenized phase fractions of the most relevant microstructural species, namely the stable β-phase, the stable α s -phase as well as the metastable Martensite α m -phase, in a physically consistent manner. In particular, the modeled microstructure evolution, in form of diffusion-based and non-diffusional transformations, is a pure consequence of energy and mobility competitions among the different species, without the need for heuristic transformation criteria as often applied in existing models. The mathematically consistent formulation of the evolution equations in rate form renders the model suitable for the practically relevant scenario of temperature- or time-dependent diffusion coefficients, arbitrary temperature profiles, and multiple coexisting phases. Due to its physically motivated foundation, the proposed model requires only a minimal number of free parameters, which are determined in an inverse identification process considering a broad experimental data basis in form of time-temperature transformation diagrams. Subsequently, the predictive ability of the model is demonstrated by means of continuous cooling transformation diagrams, showing that experimentally observed characteristics such as critical cooling rates emerge naturally from the proposed microstructure model, instead of being enforced as heuristic transformation criteria. Eventually, the proposed model is exploited to predict the microstructure evolution for a realistic selective laser melting application scenario and for the cooling/quenching process of a Ti-6Al-4V cube of practically relevant size. Numerical results confirm experimental observations that Martensite is the dominating microstructure species in regimes of high cooling rates, e.g., due to highly localized heat sources or in near-surface domains, while a proper manipulation of the temperature field, e.g., by preheating the base-plate in selective laser melting, can suppress the formation of this metastable phase.

Inverse parameter identification↗