Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “harden”

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 937 records · Page 52

A predictive discrete-continuum multiscale model of plasticity with quantified uncertainty

Multiscale models of materials, consisting of upscaling discrete simulations to continuum models, are unique in their capability to simulate complex materials behavior. The fundamental limitation in multiscale models is the presence of uncertainty in the computational predictions delivered by them. In this work, a sequential multiscale model has been developed, incorporating discrete dislocation dynamics (DDD) simulations and a strain gradient plasticity (SGP) model to predict the size effect in plastic deformations of metallic micro-pillars. The DDD simulations include uniaxial compression of micro-pillars with different sizes and over a wide range of initial dislocation densities and spatial distributions of dislocations. An SGP model is employed at the continuum level that accounts for the size-dependency of flow stress and hardening rate. Sequences of uncertainty analyses have been performed to assess the predictive capability of the multiscale model. The variance-based global sensitivity analysis determines the effect of parameter uncertainty on the SGP model prediction. The multiscale model is then constructed by calibrating the continuum model using the data furnished by the DDD simulations. A Bayesian calibration method is implemented to quantify the uncertainty due to microstructural randomness in discrete dislocation simulations (density and spatial distribution of dislocations) on the macroscopic continuum model prediction (size effect in plastic deformation). Here, the outcomes of this study indicate that the discrete-continuum multiscale model can accurately simulate the plastic deformation of micro-pillars, despite the significant uncertainty in the DDD results. Additionally, depending on the macroscopic features represented by the DDD simulations, the SGP model can reliably predict the size effect in plasticity responses of the micropillars with below 10% of error.

36 MATERIALS SCIENCE↗

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↗

A crystal plasticity finite element model embedding strain-rate sensitivities inherent to deformation mechanisms: Application to alloy AZ31

We report the fundamental power-law relationship representing the flow rule in crystal visco-plasticity ensures uniqueness in the selection of slip systems accommodating imposed plastic strain-rates. The power-law relationship also introduces an artificially high strain-rate sensitivity in crystal plasticity simulations, unless a high value of the power-law exponent is used. However, the use of high values for the exponent is limited by numerical tractability. This paper presents a numerical method implemented in a crystal plasticity finite element (CPFE) model for embedding any value of the power-law exponent reflecting the true material strain-rate sensitivity. Importantly, the method does not increase computation time involved in the simulations. The enhanced CPFE model is used to interpret and predict a complex strain-rate sensitive response and microstructural evolution of AZ31 Mg alloy. Measured values of strain-rate sensitivity for slip and twinning modes are used in the simulations. Calculations show that the model successfully captures the phenomena pertaining to the effect of changing applied strain-rate on the mechanical response including flow stress and evolution of texture and twinning for a broad range of strain-rates ranging from 10 -3 s -1 to 10 3 s -1 and loading orientations in tension and compression. It is shown that such predictions are a consequence of not only relative amounts of slip and twinning activities driven by a set of accurately characterized hardening law parameters but also values of the strain-rate sensitivities inherent to individual deformation mechanisms. Besides, the model validates the measured strain-rate dependency of deformation mechanisms while accurately reproducing the mechanical data. Hence, the model can be used to verify and further refine or infer measured strain-rate sensitivity per deformation mechanism by reproducing experimental data.

36 MATERIALS SCIENCE↗

Deformation behavior of CP-titanium under strain path changes: Experiment and crystal plasticity modeling

The deformation behavior of commercially pure rolled titanium subjected to strain path changes is studied using experiments and a crystal plasticity model. Four different loading combinations are performed at room temperature to study the activation of slip, twinning, de-twinning and double-twinning in hexagonal closed packed titanium. The strain paths considered are: rolling direction compression (RDC) followed by normal direction compression (NDC), RDC followed by transverse direction compression (TDC), NDC followed by RDC, and NDC followed by TDC. An EBSD-based analysis of the distribution of boundary misorientation angles before and after reload was developed to analyze the evolution of {$10\bar{1}2$} tensile and {$11\bar{2}2$} compression twins. This analysis supports the model results concerning the treatment of twin reorientation. A de-twinning and double-twinning model accounting for back stress effects, an important feature of strain path changes, is implemented within the framework of the visco-plastic self-consistent (VPSC) model along with a dislocation density (DD) based hardening scheme. In the model, plasticity is accommodated by prismatic $\langle$a$\rangle$, basal $\langle$a$\rangle$ and pyramidal $\langle$c+a$\rangle$ slip modes, and {$10\bar{1}2$} tensile and {$11\bar{2}2$} compression twinning modes. The VPSC model predicts the evolution of twinning, de-twinning and double-twinning processes for both tensile and compression twinning modes under strain path change. The model predicts macroscopic stress-strain response, texture evolution, and twin volume fraction that are in agreement with experimental observations. The evolution of texture is investigated in detail by separately analyzing the twinned domains, rather than the evolution of the global texture.

36 MATERIALS SCIENCE↗

Simulation of dislocation evolution in microparticle impacts over a wide range of impact velocities

We report the dynamic loading environment created by high velocity microparticle impacts is very difficult to model. The extremely high strain rates, large deformations, and high temperatures affect the flow stress of the material in ways that standard flow stress models cannot capture. Plastic deformation and strain hardening in metals is controlled by the motion of dislocations. Dislocations can be nucleated, stored as forest dislocations, or be annihilated as loading progresses. A comprehensive accounting of dislocation density change is needed to accurately describe rate and temperature dependent dislocation glide and evolution across the wide rage of loading conditions present in microparticle impact problems. Therefore, we implement and apply a newly proposed flow stress model (Hunter and Preston, 2015, 2022) that is then coupled with mobile and immobile dislocation density evolution equations. This model is implemented in Los Alamos National Laboratory’s hydrodynamics code, FLAG, to model copper-on-copper microparticle impacts. This new strength model allows for accurate simulation of particle rebound and flattening across a wide range of impact velocities.

36 MATERIALS SCIENCE↗

Development and performance of INCONEL® alloy 740H® seam-welded piping

INCONEL® alloy 740H® is an age-hardenable nickel-based alloy approved for pressure vessels and piping within ASME Section I, VIII, and B31.1. Currently, the code applies a weld strength reduction factor (WSRF) of 0.7 to the allowable stresses for longitudinal seam welds in the time dependent creep regime. In this work, a full-scale seam weld was successfully produced using typical industrial practices. The component was solution heat-treated and aged after fabrication to improve the performance of the weldment. Tensile, bend, impact, and low-cycle fatigue tests showed the component met all the specification requirements and elevated temperature properties were within base metal expectations with failures predominately in the weld metal. Long-term creep tests, including large samples more representative of the entire weldment, were fabricated and tested to times in excess of 10,000 h at multiple temperatures. Analysis of the creep data supports a WSRF of 0.9 for the solution annealed + aged weldments in contrast to the WSRF of 0.7 currently applied to welded + aged weldments. Metallurgical analysis shows that the solution annealing causes recrystallization of the weld metal and reduces the chemical compositional and microstructural gradients in the weldments resulting in acceptable short-term performance. In long-term creep, evidence for coarsened zones in the weldments due to discontinuous coarsening reactions were identified as the mechanism leading to accelerated creep damage formation in the weld metal and sample failure.

14 SOLAR ENERGY↗

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↗

Ultra-high temperature testing and performance of L-PBF C103

Additive Manufacturing (AM) of refractory alloys is gaining traction as a materials processing route for components subject to extreme temperature environments. Due to the low oxidation resistance of refractory alloys, novel methods for evaluating their elevated temperature performance must be developed. In this work, a Gleeble® 3800 thermomechanical load frame was modified to evaluate the mechanical properties of laser powder bed fusion (L-PBF) consolidated niobium alloy C103 ranging from room temperature (RT) to 1400 °C. The fixturing and sample geometry were designed to accommodate Joule heating and prevent damage to the test chamber. Oxidation of the samples was minimized via testing in vacuum level of 1E-5 Torr. Ultimate tensile strength (UTS), yield strength (YS), elongation, and strain-hardening behavior were determined as a function of temperature. L-PBF C103 presented an average UTS of ∼650 MPa and over 25 % elongation at RT. Above RT, the UTS and YS dropped then leveled off from 500 °C to 1000 °C with values ranging from ∼400 MPa to ∼460 MPa, which is consistent with dynamic strain aging observed in this class of alloys. The strength rapidly declined after 1200 °C to ∼150 MPa at 1400 °C. Fractography indicated ductile fractures for the C103 at all test temperatures, and Electron Backscatter Diffraction (EBSD) analysis revealed a textured microstructure and the presence of dynamic recrystallization within the necked region of the sample tested at 1400 °C.

33 ADVANCED PROPULSION SYSTEMS↗

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↗

Modeling the viscoplastic behavior of a semicrystalline polymer

In this study, a complex constitutive relation is identified using inverse modeling with the nominal mechanical response as sole experimental input. The methodology is illustrated for a semicrystalline thermoplastic in the presence of strain localization at finite deformations. The experimental database includes cylindrical tensile bars, compression pins and round notched bars loaded at strain rates spanning up to five decades and temperatures below and above T g . The data is organized into a calibration set and a validation set. The response of tensile specimens is determined using finite element analyses and a two-phase constitutive relation for semicrystalline polymers that accounts for temperature- and rate-sensitive plastic flow, pressure-sensitivity, small-strain softening and large-strain orientational hardening of the amorphous phase, along with the evolution of crystallinity. The large number of constitutive parameters is identified using an optimization tool coupled with the finite element solver and the calibration set from experiments. The methodology is shown to be successful in predicting the response of round notched bars and replicating the effects of temperature and strain rate on the severity of necking in tensile bars. The proposed model identification strategy is both simple and effective in comparison with other elaborate methods that attempt to access intrinsic behavior directly from high-fidelity experimental measurements.

36 MATERIALS SCIENCE↗

Optimization of the microstructure and mechanical properties of electron beam welded high-strength medium-entropy alloy (NiCoCr) 94 Al 3 Ti 3

For this work, the weldability and effects of subsequent heat treatment on the microstructure and mechanical properties of the medium-entropy alloy (MEA) (NiCoCr) 94 Al 3 Ti 3 have been explored using electron beam welding (EBW). It was found that EBW joints consist of a supersaturated f.c.c. matrix with coarse, asymmetric columnar grains (17.7–24 μm). The fraction of twins and Σ3 n CSL boundaries in weld joints was only 10% of that in the base metal, a value which increased to 30% after aging at 800 °C. The aging also led to improvements in the hardness, tensile strength, and work hardening rate, a result of the precipitation of a high volume fraction (0.29–0.33) of ultrafine (15–17.4 nm dia.) coherent spherical L1 2 -structured particles. The contributions of the L1 2 nanoparticles and Hall-Petch strengthening to the yield strength were estimated to be 537–573 MPa and 105–123 MPa, respectively. The MEA shows excellent weldability and, after aging, good strength (YS~1098 MPa, UTS~1368 MPa) and elongation to failure (~11%) by ductile fracture.

36 MATERIALS SCIENCE↗

High susceptibility to adiabatic shear banding and high dynamic strength in tungsten heavy alloys with a high-entropy alloy matrix

The susceptibility to adiabatic shear banding (ASB) of tungsten heavy alloys (WHAs) is a key factor to determine their penetration performance for kinetic energy penetrators (KEPs). This work reports a novel WHA with strong susceptibility to ASB that uses a precipitation-hardening high-entropy alloy (HEA) as the binding phase of W particles. The results find that nanoscale L12 precipitates and high dislocation density in the matrix contribute a high dynamic yield strength of ~2300 MPa. Also, a narrow ASB with a width of ~12 μm forms in the novel WHA, which is only 6% of the width of ASB in traditional WHAs like W-NiFe, exhibiting a great susceptibility to ASB. The redissolution of precipitates into the matrix because of temperature rise during shearing could sharply soften the shearing location, which is suggested to be a very important factor of destabilizing mechanism to promote the development of ASB.

36 MATERIALS SCIENCE↗

Investigation of crossed-twin structure formation in magnesium and magnesium alloys

In this work, the effect of alloying addition on the propensity for twin-twin interactions to transform into crossed-twin structures in magnesium alloys is investigated. A full-field elasto-viscoplastic fast Fourier transform (EVP-FFT) framework combined with a discrete twin model and dislocation density-based hardening law for slip strengths is used to calculate the micromechanical fields in the crystalline matrix around the interacting twins. AZ31 and MgLi alloys are selected along with pure Mg to study the influence of plastic anisotropy in connection with alloying elements. These alloys were selected since their plastic anisotropy measure, which is defined as the ratio between the critical resolved shear stress for pyramidal $\langle$c +a$\rangle$ and basal $\langle$a$\rangle$ slip modes, spanned a wide range. To quantify the role of twin thicknesses, we probe a range of impinging twin thicknesses while fixing the recipient twin thickness. The analysis reveals that: (i) the local driving stress for crossed twin structure formation generated from the interaction of the two twins is lower in a low plastically anisotropic material, like a MgLi alloy, than a high plastically anisotropic material like pure Mg and (ii) the critical impinging twin thickness needed to form the crossed twin structure in pure Mg, AZ31 and MgLi alloys is ~0.5, ~0.75 and ~1.5 times the recipient twin thickness. We propose a relationship between the tendency for crossed twin structure formation and the experimentally observed higher ductility in MgLi alloys compared to pure Mg. Finally, one key implication of the findings is that crossed- twin structure formation can be hindered and the ductility of magnesium alloy thereby improved by properly choosing alloying elements that lower the slip strength for pyramidal $\langle$c + a$\rangle$ slip.

36 MATERIALS SCIENCE↗

Gas nitriding behavior of refractory metals and implications for multi-principal element alloy design

Multi-principal element alloys (MPEAs) comprise a large, flexible compositional space that enables tuning of their chemistry, structure, and properties. Here, to facilitate the development of nitriding-based surface-enhancement strategies that harness a broad compositional space, this study examined the gas nitriding behavior of Hf, Mo, Nb, Ta, Ti, and Zr as a function of time, temperature (750 and 1000 °C), and nitriding potential (i.e. ammonia-to-hydrogen ratio). These metals were selected because they have a strong driving force to form nitrides, and appear in many promising refractory MPEA compositions. The nitriding temperatures were selected based on the phase transformation temperature of Ti and Zr, and the nitriding potentials were chosen such that all elements are expected to form nitrides. Mass gain measurements indicate that all six elements follow parabolic kinetics. The microstructure observations and quantitative microchemical analysis show formation of dense and well-adhered compound layers for Mo, Nb, and Ta. Thick diffusion zones appear in Hf, Ta, Ti, and Zr, and diffusion coefficients were fit to the composition profiles. Partial delamination of the compound layer occurred for Ti and Zr. Peak hardness values above 30 GPa are obtained in the dense compound layers, and the solute hardening of the underlying alloy is correlated with the nitrogen content. The results provide insight into the dynamics of nitride compound formation relative to interstitial dissolution of nitrogen, and are discussed in the context of MPEA composition and processing design.

36 MATERIALS SCIENCE↗

Electrochemical Li intercalation in b-As y P 1–y alloys: In-situ Raman spectroscopy study

Structural evolution of b-As y P 1-y alloys with varying As concentration, y was studied during Li-intercalation using in-situ Raman spectroscopy. A monotonic redshift of all the vibrational modes corresponding to P-P bonds (A 1 g , A 2 g , B 2g ), As-As bonds (A 1 g , A 2 g , B 2g ), and As-P bonds of all the samples was observed during the initial stages of the intercalation process due to the softening of each mode caused by the intercalation driven donor-type charge-transfer from Li to b-As y P 1-y . The emergence of a new peak identified as the E g mode of gray As, above an intercalation threshold, is indicative of the presence of an intercalation-driven structural phase segregation process. Further, the A 1g mode of gray As emerging after this intercalation threshold lies in close proximity of the A 2 g Raman mode of b-As y P 1-y . All the peaks beyond the intercalation threshold show an upshift due to the co-existence of gray As with b-As y P 1-y alloys causing strain and thus hardening of the phonon modes. In the sample with the highest As concentration phase segregation takes place during the synthesis process. Computational modeling reveals the co-existence of gray As in the b-As y P 1-y alloys with high As concentrations. It also confirms the existence of a local structural segregation taking place at a critical Li concentration during the intercalation process. Furthermore, this work shows the significance of intercalation on structural changes in the search for new phases of b-As y P 1-y alloys.

2D Materials↗

A moving discontinuous Galerkin finite element method with interface condition enforcement for compressible flows

A variation of moving discontinuous Galerkin finite element method with interface condition enforcement (MDG-ICE) is developed for solving the compressible Euler equations. The MDG-ICE method, originating from the work of Corrigan et al. [1], [2], [3], [4], is based on the space-time DG formulation, where both flow field and grid geometry are considered as independent variables and the conservation laws are enforced both on discrete elements and element interfaces. The element conservation laws are solved in the standard discontinuous solution space to determine conservative quantities, while the interface conservation is enforced using a variational formulation in a continuous space to determine discrete grid geometry. The resulting over-determined system of nonlinear equations arising from the MDG-ICE formulation can then be solved in a least-squares sense, leading to an unconstrained nonlinear least-squares problem that is regularized and solved by Levenberg-Marquardt method. A number of numerical experiments for both 1D unsteady and 2D steady state compressible flow problems are conducted to assess the accuracy and robustness of the MDG-ICE method. Numerical results obtained indicate that the MDG-ICE method is able to implicitly detect and track all types of discontinuities via interface conservation enforcement and satisfy the conservation law on both elements and interfaces via grid movement and grid management, demonstrating that an exponential rate of convergence for Sod and Lax-Harden shock tube problems can be achieved and highly accurate solutions without overheating to both double-rarefaction wave and Noh problems can be obtained.

97 MATHEMATICS AND COMPUTING↗

Helium effects and bubbles formation in irradiated Ti 3 SiC 2

Ti 3 SiC 2 is a potential structural material for nuclear reactor applications. However, He irradiation effects in this material are not well understood, especially at high temperatures. Here, we compare the effects of He irradiation in Ti 3 SiC 2 at room temperature (RT) and at 750 °C. Irradiation at 750 °C was found to lead to extremely elongated He bubbles that are concentrated in the nano-laminate layers of Ti 3 SiC 2 , whereas the overall crystal structure of the material remained intact. In contrast, at RT, the layered structure was significantly damaged and highly disordered after irradiation. Our study reveals that at elevated temperatures, the unique structure of Ti 3 SiC 2 can accommodate large amounts of He atoms in the nano-laminate layer, without compromising the structural stability of the material. The structure and the mechanical tests results show that the irradiation induced swelling and hardening at 750 °C are much smaller than those at RT. Overall, these results indicate that Ti 3 SiC 2 has an excellent resistance to accumulation of radiation-induced He impurities and that it has a considerable tolerance to irradiation-induced degradation of mechanical properties at high temperatures.

36 MATERIALS SCIENCE↗

Design, performance, and cost savings of using GF-PC additively manufactured tooling for stamping of HSS 590 sheet metal

Conventional sheet metal forming tooling in the automotive industry is made up of hardened steel and used for mass-production. Prototype tooling made of metal is durable, but it is only used for a small number of parts despite its high cost, contributing heavily to the vehicle development cost. Additive Manufacturing (AM) offers a low-cost alternative for fabricating tooling suitable for low-volume production. Here, this study investigates the feasibility of using AM polymer composite tooling for the stamping of 1.5mm thickness HSS 590 steel sheets through a two-pronged approach – experimental and numerical analysis. Sheet metal stamping experiments were conducted with glass fiber reinforced polycarbonate (GF-PC) AM polymer tooling and their performance was evaluated based on various metrics such as tool deformation and part accuracy. Finite element simulations of the stamping process were performed, which accurately captured the tool deformations seen in experiments with the use of anisotropic material models. The effective use of simulations in optimizing process parameters to achieve the desired final part geometry is also demonstrated. The GF-PC AM polymer composite tooling was found to be well suited for low volume production applications with potential cost savings.

36 MATERIALS SCIENCE↗