Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “STRAIN HARDENING”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Mechanical properties and deformation mechanisms of single crystal Mg micropillars subjected to high-strain-rate C-axis compression

Here, the mechanical properties and deformation mechanisms of single crystal magnesium under c-axis quasi-static and high-strain rate compressions are investigated through in situ scanning electron microscope (SEM) experiments and post-mortem transmission electron microscope (TEM) characterization. The findings revealed that ductility and high rates of hardening are preserved for pillars as large as 15 μm. Furthermore, rate effects result in a mild increase in flow stress with plastic deformations controlled primarily by the slip of type dislocations. Importantly and in contrast to other literature reports, plastic deformation occurs in the absence of twining. As the strain increases and plastic deformation exceeds about 4%, crystal rotation activates basal slip, <$\mathrm{a}$> type dislocations, resulting in a more rate independent flow stress. TEM observation on micropillars compressed at a strain rate of 250/s, revealed the activation of {${11}$$\bar{2}$$\bar{2}$} < $\bar{1}$$\bar{1}23$ > slip systems and high mobility of screw dislocations as major contributors to plastic strains in excess of 10% without fracture. These findings are relevant to the design of lightweight materials used in transportation systems, e.g., selection of material grain size. Moreover, the experimental data here reported provides the materials science community with a unique opportunity to validate discrete dislocation dynamics (DDD) formulations employed in multiscale design of materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Unraveling kinking: A plasticity enhancing failure mode in high strength nano metallic laminates

Kinking is an important and plasticity-enhancing deformation/failure mode in numerous mechanically anisotropic materials including high-strength nano metallic laminates (NMLs). However, our current limited understanding of the mechanics of kinking and its dependence on microstructural attributes is insufficient for thoroughly comprehending and eventually being able to control failure behaviors of materials. In this study, we investigate kinking dependencies on microstructural attributes in NMLs via in situ micropillar compression, multiscale microstructure characterization, dislocation dynamic simulations, and crystal plasticity modeling. Additionally, by examining several NML systems (Cu/Fe, Ag/Fe, Al-4Mg/Fe), we demonstrate that the development of internal stresses during loading activates local layer-parallel glide triggering kinking in NMLs. Furthermore, this work reveals the effect of key microstructural features including layer thickness, layer waviness, interface barrier strength, and work hardening capacity on kink band formation in NMLs. More broadly, our efforts represent a generically applicable approach for probing large-strain deformation behavior of complex materials via synergetic modeling and experimental efforts.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Martensitic transformation induced strength-ductility synergy in additively manufactured maraging 250 steel by thermal history engineering

Maraging steels are known for their exceptional strength but suffer from limited work hardening and ductility. Here, in this study, we report an intermittent printing strategy to tailor the microstructure and mechanical properties of maraging 250 steel via tuning the thermal history during wire-arc directed energy deposition. By introducing a dwell time between adjacent layers, the maraging 250 steel is cooled below the martensite start temperature, triggering thermally-driven martensitic transformation during the printing process. Thermal cycling during subsequent layer deposition results in the formation of reverted austenite which shows a refined microstructure and induces elemental segregation between martensite and reverted austenite. The Ni enrichment in the austenite promotes stabilization of the reverted austenite upon cooling to room temperature. The reverted austenite is metastable during deformation, leading to strain-induced martensitic transformation under loading. Specifically, a 3 min interlayer dwell time produces a maraging 250 steel with approximately 8% reverted austenite, resulting in improved work hardening via martensitic transformation induced plasticity during deformation. Meanwhile, the higher cooling rate and refined prior austenite grains lead to substantially refined martensitic grains (by approximately fivefold) together with an increased dislocation density. With 3 min interlayer dwell time, the yield strength of the printed maraging 250 steel increases from 836 MPa to 990 MPa, and the uniform elongation is doubled from 3.2% to 6.5%. This intermittent deposition strategy demonstrates the potential to tune the microstructure of maraging steels for achieving strength-ductility synergy by engineering the thermal history during additive manufacturing.

Additive manufacturing↗
DOE OSTI2020

Proppant embedment in coal and shale: Impacts of stress hardening and sorption

During methane production in CBM reservoirs, the influence of proppant embedment and permeability damage cannot be neglected – especially where the wall-rock is soft. Effective stresses are elevated during methane recovery, increasing both normal loading stress and confinement and simultaneously overprinting sorption-induced volumetric strains. Experiments and analytic modeling are conducted to define key mechanisms controlling these competitive effects. We independently measure overall sample compaction (external LVDT) and local strain (strain gauge) in the matrix to deconvolve proppant embedment in a propped fracture for different conditions of confining stress. The results show symptomatic behaviors of elastic (shale) and elastoplastic (coal) responses of embedment. Different from shale, the evolution of embedment is convex upwards with increased stress where indented depth increases more rapidly as loading stress increases under constant confinement. In addition, a stress-hardening effect is found to play a pivotal role in determining the characteristics of indentation, which are examined in terms of evolution profiles, deformation regimes, embedment slopes, curvatures, yield points and irreversible indentations. Based on the experimental observations a semianalytical model predicts indentation and the evolution of propped permeability under recreated in-situ stress conditions. A simplified case study is conducted to further illustrate the evolution of aperture and permeability of a propped fracture in CBM reservoirs. The modeling results suggest that proppant embedment is significantly overestimated if the variable stress-hardening (VSH) effect is neglected, especially when effective stress is large. Moreover, a decrease in indentation depth possibly occurs during late stage methane production, resulting in a reversal/recovery in fracture closure. This is because desorption-induced shrinkage becomes the predominant effect, causing an increase in aperture and a reduction in the indented volume of proppant. The resulting recovery in permeability implies that the propped coal fracture has the potential to optimally facilitate methane production as a pathway, even at high closure stresses generated by methane drainage.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2021

On the computational solution of vector-density based continuum dislocation dynamics models: A comparison of two plastic distortion and stress update algorithms

Continuum dislocation dynamics models of mesoscale plasticity consist of dislocation transport-reaction equations coupled with crystal mechanics equations. The coupling between these two sets of equations is such that dislocation transport gives rise to the evolution of plastic distortion (strain), while the evolution of the latter fixes the stress from which the dislocation velocity field is found via a mobility law. Earlier solutions of these equations employed a staggered solution scheme for the two sets of equations in which the plastic distortion was updated via time integration of its rate, as found from Orowan's law. In this work, we show that such a direct time integration scheme can suffer from accumulation of numerical errors. We introduce an alternative scheme based on field dislocation mechanics that ensures consistency between the plastic distortion and the dislocation content in the crystal. The new scheme is based on calculating the compatible and incompatible parts of the plastic distortion separately, and the incompatible part is calculated from the current dislocation density field. Stress field and dislocation transport calculations were implemented within a finite element based discretization of the governing equations, with the crystal mechanics part solved by a conventional Galerkin method and the dislocation transport equations by the least squares method. A simple test was first performed to show the accuracy of the two schemes for updating the plastic distortion, which shows that the solution method based on field dislocation mechanics is more accurate. This method then was used to simulate an austenitic steel crystal under uniaxial loading and multiple slip conditions. By considering dislocation interactions caused by junctions, a hardening rate similar to discrete dislocation dynamics simulation results was obtained. Finally, the simulations show that dislocations exhibit some self-organized structures as the strain is increased.

36 MATERIALS SCIENCE↗
DOE OSTI2026

TEM Approaches for Microstructure-Informed Prediction of Mechanical Properties in Structural Alloys

Predicting the mechanical performance of structural alloys from their evolving microstructure remains a major challenge in materials science, particularly for nuclear structural materials, where irradiation-induced defects span multiple types and length scales and interact through complex mechanisms. The dispersed barrier hardening (DBH) [1] and Friedel–Kroupa–Hirsch (FKH) [2,3] models have been widely used to evaluate the hardening contributions of individual obstacles and to estimate tensile strength from quantified microstructures; however, when multiple size-dependent obstacles coexist and evolve, predicting temperature-dependent tensile strength becomes significantly more complex, and a fully consistent hardening model is still lacking. Transmission electron microscopy (TEM) plays a central role in refining hardening models and enabling predictive assessments of tensile strength evolution by providing quantitative characterization of dislocations, irradiation-induced defects (e.g., dislocation loops and cavities), precipitates, and grain structure (Fig. 1.). These experimentally measured defect densities are incorporated into physically based hardening models with size- and shape- dependent obstacle strengths [4], using root-sum-square superposition for obstacles of comparable strength and linear superposition for dissimilar ones [5]. In addition, recent advances in TEM [6-8], including high-resolution imaging, 4D-STEM strain mapping, EDS/EELS elemental analysis, and flash-polishing-based TEM specimen preparation and extraction-replica methods (Fig. 2), further improve the accuracy of microstructural quantification. By comparison with prior studies as well as our own results, we show that when TEM-derived microstructural information is carefully integrated with physically grounded hardening models, yield strength (or irradiation-induced hardening) measured at room temperature can be predicted with good quantitative agreement across multiple alloy classes. In-situ TEM combined with high-temperature mechanical testing represents an important next step for refining hardening models by directly probing dislocation–obstacle interactions across varying irradiation doses and temperatures [9]. Because the barrier strength factor (α) depends on both temperature and obstacle size, it should not be treated as a constant fitting parameter; rather, it must be explicitly evaluated to achieve physically meaningful predictions of mechanical behaviour at operating temperatures. This presentation therefore discusses why all strengthening contributions (e.g., Peierls stress, solid-solution strengthening, voids, bubbles, dislocation loops, dislocation lines, and grain boundaries) must be considered collectively, why appropriate superposition methods are essential when obstacles possess different barrier strength factors, how hardness measurements can be meaningfully related to tensile properties, and how TEM-derived microstructural information can be systematically incorporated into hardening models. More broadly, it outlines a pathway toward microstructure-informed prediction of mechanical properties and supports the goal of establishing science-based tools for evaluating structural materials in extreme environments [10].

Lin, Yan-Ru [ORNL] (ORCID:0000000339991473)↗
DOE OSTI2021

Shear-based deformation processing of age-hardened aluminum alloy for single-step sheet production

Shear-based deformation processing by hybrid cutting-extrusion and free machining are used to make continuous strip, of thickness up to one millimeter, from low-workability AA6013-T6 in a single deformation step. The intense shear can impose effective strains as large as 2 in the strip without pre-heating of the workpiece. The creation of strip in a single step is facilitated by three factors inherent to the cutting deformation zone: highly confined shear deformation, in situ plastic deformation-induced heating and high hydrostatic pressure. The hybrid cutting-extrusion, which employs a second die located across from the primary cutting tool to constrain the chip geometry, is found to produce strip with smooth surfaces (S a < 0.4 μm) that is similar to cold-rolled strip. The strips show an elongated grain microstructure that is inclined to the strip surfaces – a shear texture – that is quite different from rolled sheet. Furthermore, this shear texture (inclination) angle is determined by the deformation path. Through control of the deformation parameters such as strain and temperature, a range of microstructures and strengths could be achieved in the strip. When the cutting-based deformation was done at room temperature, without workpiece pre-heating, the starting T6 material was further strengthened by as much as 30% in a single step. In elevated-temperature cutting-extrusion, dynamic recrystallization was observed, resulting in a refined grain size in the strip. Implications for deformation processing of age-hardenable Al alloys into sheet form, and microstructure control therein, are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Additive friction stir deposition of SS316: Effect of process parameters on microstructure evolution

Solid state nature of additive friction stir deposition (AFSD) additive manufacturing process is very advantageous in terms of defect formation and microstructural refinement in the material. Current study presents the process optimization, microstructural evolution and kinetics of recrystallization for AFSD deposited low stacking fault energy material - SS316. As deposited microstructure shows equiaxed ultra fine grains with an average grain size of ~5.0 ± 0.5 μm. Shear deformation at high temperature during processing leads to the operation of restoration mechanisms. Observation of necklace type microstructure in the as deposited SS316 is attributed to discontinuous dynamic recrystallization during processing. Recrystallization kinetics of the AFSD SS316 is characterized using Johnson-Mehl-Avarami-Kolmogorov (JMAK) model. Deformation – thermal cycling during AFSD process resulted in inconsistent recrystallization kinetics. Variation in strain, strain rate and temperature during processing with processing parameters result in varying microstructure and tool wear. Further, high strength-ductility combination and sustained work hardening in as deposited SS316 appear to arise from transformation and twinning during deformation, leading to the formation of hierarchical twins and martensitic phase after deformation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

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

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

Crystal plasticity finite element method↗
DOE OSTI2021

Synergic grain boundary segregation and precipitation in W- and W-Mo-containing high-entropy borides

The structures of W- and W-Mo-containing high-entropy borides (HEBs) are systematically studied by combining atomic-resolution transmission electron microscopy imaging, electron diffraction, and chemical analysis. Here, we reveal that W or W-Mo addition in HEBs leads to segregation of these elements to the grain boundaries (GBs). In the meantime, W- or W-Mo-rich precipitates also form along the GBs. Crystallographic analysis and atomic-scale imaging show that the GB precipitates in both W- and W-Mo-containing HEBs have a cube-on-cube orientation relationship with the matrix. With further strain analysis, the coherency of the precipitate/matrix interface is validated. Nanoindentation tests show that the simultaneous GB segregation and coherent precipitation, as a supplement to the grain hardening, provide additional hardening of the HEBs. Our work provides an in-depth understanding of the GB segregation and precipitation behaviors of HEBs. It suggests that GB engineering could be potentially used for optimizing the performance of high-entropy ceramics.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2022

Damage Detection and Electrical Performance Impact of Flat-Flexible Circuits

Flat-flexible circuits or ‘flex circuits’ have varied applications in electronics and in explosive detonation systems. They pose advantages over conventional coaxial cables in their ability to fit into small mechanical envelopes and routings. A typical flex circuit is composed of a laminated assembly of multiple layers of conductor and dielectric materials. Copper, because of its high electrical conductivity, is one of the primary materials used in flex circuit conductors. The conductor must maintain structural integrity and electrical continuity when the flexible circuit is bent during processing or to fit an application envelope. Under some strain conditions a cable may have degraded electrical performance even though visible damage may not be evident. Dislocation density can increase by mechanisms such as dislocation forest hardening, where gliding dislocations in primary slip planes get pinned at obstacles and form Frank-Read sources. Dislocation density increases within the conductor layers may impact electrical resistivity but would not necessarily result in visible mechanical failure. In some detonation applications, such as for building implosion for demolition, the time for signal transmission and voltage used for firing detonators can be critical. If the resistance in a single flex cable within a building demolition array is too high, the resultant demolition shape may be asymmetric, leading to unplanned damage to nearby structures and potential hazards to public safety. This study will apply non-destructive evaluation techniques to evaluate the detectability of damage accumulation in flat, flexible copper-conductor cables and discuss potential implications to their practical use.

42 ENGINEERING↗
DOE OSTI2022

Study of the interplay between lower-order and higher-order energetic strain-gradient effects in polycrystal plasticity

in this report strain-gradient (SG) plasticity refers to a class of non-local theories in which gradients of plastic slip determine the storage of geometrically necessary dislocations, introducing a length-scale dependence in the mechanical behavior of crystalline materials, which is otherwise lacking in local theories. In this work, we incorporate lower-order (LO) and higher-order energetic (HOE) strain-gradient effects into a crystal plasticity fast Fourier transform (FFT)-based formulation to investigate the interplay of the length scale that each strain-gradient term introduces at the microscale, and the mechanical properties that result at the macroscale. For an applicable range of length scales, we consider two systems: a 1-D two-phase face centered cubic (FCC) laminate and a 3-D FCC polycrystal, and two uniaxial deformation modes: monotonic tension and cyclic tension–compression. We show that increases in the individual LO and HOE length scales increase the hardening rate and strength of the material, respectively. When combined, the strong LO hardening is less pronounced than the effect alone due to the lowering of the gradients due to the HOE microstress. We demonstrate that the LO and HOE hardening manifest as “isotropic” (yield surface expansion) and “kinematic” (yield surface shift) effects, respectively, consistent with their theoretical origins. We show that in cyclic loading, the Bauschinger effect emerges in both local and non-local calculations and link its origins and severity to the behavior in the strain field, slip-system rates, and the HOE microforce.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Unveiling the origins of work-hardening enhancement and mechanical instability in laser shock peened titanium

Laser shock-peening (LSP) produces complex gradients in microstructures and residual stresses, which enables fabrication of metallic engineering components with superior mechanical properties. In this work, we use non-destructive and spatially resolved high-energy X-ray diffraction (HE-XRD) to investigate the influences of LSP on microstructure, surface topography, and residual stress for high purity titanium plates. LSP is found to produce large compressive in-plane residual stresses near the peened surface which monotonically decay to zero about 2.5 mm below the surface. These properties were also tracked during in-situ tensile loading, allowing stress partitioning and work-hardening rates to be measured as a function of depth. The surface region is found to have the highest work-hardening rate and remain mechanically stable up to sample failure. In addition, a crystal rigid rotation along a transversal axis of about 27 & nbsp; was observed at large applied strain near the surface region, which is attributed to the formation of denser dispersed shear bands. Furthermore, a sudden drop in tensile stress was observed at a depth of similar to 160 mu m from the LSP surface, indicating the formation of a localized mechanical instability zone (LMIZ). The origin of this LMIZ could be attributed to the increase in shear stress induced by a load transfer from LSP surface due to the formation of denser dispersed shear bands. These observations indicate a localized transition in deformation behavior has occurred from dispersed shear banding to LMIZ and homogeneous deformation.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Dose and dose rate effects on the microstructural and mechanical stability of long-range ordered precipitates in Inconel 718

Nickel-based alloys offer exceptional mechanical properties and corrosion resistance for various applications in extreme environments such as high neutron flux, temperature, and stress. Precipitates, resulting from thermal treatments, increase the strength and creep resistance of these alloys. Here, heavy ion irradiation is shown to destabilize the γ' and γ" precipitates in Inconel 718 nickel-based alloy by disordering and dissolution processes. Transmission electron microscopy and atom probe tomography reveal that the disordering and dissolution rates of the precipitates are proportional to both irradiation dose and dose rate. Irradiation induced disordering occurs at a higher rate than the dissolution. From analyzing the hardening contributions, microscale mechanical results suggest radiation-induced disordering is responsible for the majority of the softening response at low dose. In addition, low dose radiation-induced disordering results in strain localization and is believed to be the primary cause for the drastic degradation of uniform elongation and toughness in irradiated Inconel 718.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Examining deformation localization of irradiated tungsten under uniaxial compression with crystal plasticity

In this work, the deformation behavior of irradiated tungsten at elevated service temperatures is examined using mechanism-based crystal plasticity finite element framework by considering multiple mechanisms, including thermal softening and irradiation hardening. The previously developed mechanism-based crystal plasticity model is first utilized to assess the yield stress of unirradiated and irradiated tungsten upon low-temperature irradiation. Additionally, the engineering stress-strain response under uniaxial compression is then predicted by considering the effect of irradiation-induced defect clusters on the motion of deformation-induced dislocations, and the predicted results agree well with the experimental measurements. More importantly, the formation and evolution of plastic strain localization zones in the form of shear bands are qualitatively captured, and the maximal principal strain distribution is quantitatively linked to the dose of irradiation damage. Finally, the deformation-induced shear band as a result of plastic instability in the compressed tungsten specimen is discussed.

36 MATERIALS SCIENCE↗
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 OSTI2022

Strengthening Mechanisms and Deformation Behavior of Industrially-Cast and Lab-Cast Dual-Phase High Entropy Alloy

An iron-rich high entropy alloy Al 0.65 CoCrFe 2 Ni comprising of both FCC and BCC phases was designed and fabricated using induction melting at a large industrial scale and suction casting at a small lab scale. Spinodally decomposed interdendritic regions were uniformly distributed in the industrially melted alloy, and a finer dendritic structure with smaller sized spinodal structures was observed in the alloy suction cast at lab scale. The suction cast alloy fabricated at the lab scale exhibited significantly higher compressive yield strength compared to the industrially cast alloy. The relative strengthening contributions in the alloys due to solid solution, grain boundary, dislocation interactions and interphase boundary were analyzed. Dislocation hardening was the major contributor responsible for the higher strength of the SC alloy. The deformation behavior in the alloys was examined with the help of orientation image mapping. Finally, intense slip band formation at large strains in the industrially melted alloy indicates that cold working can be an effective technique to enhance the strength of the dual-phase high entropy alloy processed at a large scale.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Gas Bubble Evolution in Polycrystalline UMo Fuels Under Elastic-Plastic Deformation: A Phase-Field Model With Crystal-Plasticity

In monolithic UMo fuels, the interaction between the Al cladding and large gas bubble volumetric swelling causes both elastic-plastic and creep deformation. In this work, a phase-field model of gas bubble evolution in polycrystalline UMo under elastic-plastic deformation was developed for studying the dynamic interaction between evolving gas bubble/voids and deformation. A crystal plasticity model, which assumes that the plastic strain rate is proportional to resolved shear stresses of dislocation slip systems on their slip planes, was used to describe plastic deformation in polycrystalline UMo. Xe diffusion and gas bubble evolution are driven by the minimization of chemical and deformation energies in the phase-field model, while evolving gas bubble structure was used to update the mechanical properties in the crystal plasticity model. With the developed model, we simulated the effect of gas bubble structures (different volume fractions and internal gas pressures) on stress-strain curves and the effect of local stresses on gas bubble evolution. The results show that 1) the effective Young’s modulus and yield stress decrease with the increase of gas bubble volume fraction; 2) the hardening coefficient increases with the increase of gas bubble volume fraction, especially for gas bubbles with higher internal pressure; and 3) the pressure dependence of Xe thermodynamic and kinetic properties in addition to the local stress state determine gas bubble growth or shrinkage. The simulated results can serve as a guide to improve material property models for macroscale fuel performance modeling.

36 MATERIALS SCIENCE↗