Engineering Papers⌕ Search

Engineering topics

Bronkhorst, Curt Allan

Publications and source records attributed to Bronkhorst, Curt Allan.

At least 19 records

Anisotropic thermal-conductivity degradation in the phase-field method accounting for crack directionality

Dynamic loads applied to metals may lead to brittle or ductile fracture depending on the imposed strain rates, material properties, and specimen geometry. With an increase in the applied velocity, brittle to ductile failure mode transition may also be observed. Furthermore, at high strain rates, ductile fracture may be preceded by shear bands, which are narrow bands of intense plastic deformation, typically accompanied by a significant rise in temperature. Reliable models are needed to predict the response of metals subject to dynamic loads. In particular, capturing the interplay between heat conduction and crack propagation using the phase-field fracture method is still an open research field. To accurately capture the heat transfer physics across crack surfaces, damage models degrading thermal-conductivity are necessary. While isotropic thermal-conductivity degradation models were proposed, they may lead to errors as they are indifferent to crack directionality. In our current work, a new anisotropic approach is proposed in which thermal-conductivity, which depends on the phase-field gradient, is degraded solely across the crack. Here, it is shown that this approach improves the near-field approximation of temperature and heat flux compared with isotropic degradation, when taking the discontinuous crack solutions as reference. Additionally, the anisotropic degradation approach is implemented in a unified dynamic fracture model and a couple of examples demonstrate the viability of this approach.

42 ENGINEERING↗

Thermomechanical conversion in metals: dislocation plasticity model evaluation of the Taylor-Quinney coefficient

Using a partitioned-energy thermodynamic framework which assigns energy to that of atomic configurational stored energy of cold work and kinetic-vibrational, in this study we derive an important constraint on the Taylor-Quinney coefficient, which quantifies the fraction of plastic work that is converted into heat during plastic deformation. Associated with the two energy contributions are two separate temperatures – the ordinary temperature for the thermal energy and the effective temperature for the configurational energy. We show that the Taylor-Quinney coefficient is a function of the thermodynamically defined effective temperature that measures the atomic configurational disorder in the material. Finite-element analysis of recently published experiments on the aluminum alloy 6016-T4 [1], using the thermodynamic dislocation theory (TDT), shows good agreement between theory and experiment for both stress-strain behavior and temporal evolution of the temperature. The simulations include both conductive and convective thermal energy loss during the experiments, and significant thermal gradients exist within the simulation results. Computed values of the differential Taylor-Quinney coefficient are also presented and suggest a value which differs between materials and increases with increasing strain.

36 MATERIALS SCIENCE↗

Three-dimensional modeling and simulations of single-crystal and bi-crystal titanium for high-strain-rate loading conditions

High purity single crystal titanium (Ti) under shock wave loading is modeled under both one-dimensional and three-dimensional cylindrical conditions. Cylinder sizes of 10 μm and 20 μm radius are both considered in order to assess influence of boundary conditions. A thermodynamically consistent single crystal model for application to shock conditions is presented. The model accounts for the coupled non-linear elastic, dislocation slip, deformation twinning, and structural phase transformation response of the titanium material. Plate impact experiment results using a copper flyer are used to compare against the simulations for crystals oriented in [0001] and [10$\overline1$1] crystallographic directions. The one-dimensional and threedimensional simulations of the two differently oriented single crystals indicate differences between the one-dimensional and three-dimensional representation, especially for the [10$\overline1$1] oriented single crystal. This orientation breaks the relative orientation symmetry between the crystal and cylinder which otherwise exists for the [0001] oriented single crystal. A significant amount of heterogeneity in the field response of the[10$\overline1$1] oriented simulation was demonstrated due to the highly coupled nature of the deformation. A bi-crystal model composed of both the [0001] and [10$\overline1$1] orientations with the boundary between the two along the axis of the cylinder is also considered for a cylinder model size of 10 μm. The results indicate a strong interaction between the two grains that affects the ω phase volume fraction achieved relative to the single crystal calculations.

36 MATERIALS SCIENCE↗

Thermodynamic theory of crystal plasticity: Formulation and application to polycrystal fcc copper

We present a thermodynamic description of crystal plasticity. Our formulation is based on the Langer-Bouchbinder-Lookman thermodynamic dislocation theory (TDT), which asserts the fundamental importance of an effective temperature that describes the state of configurational disorder and therefore the dislocation density of the crystalline material. We extend the TDT description from isotropic plasticity to crystal plasticity with many slip systems. Finite-element simulations show favourable comparison with experiments on polycrystal fcc copper under uniaxial compression, tension, and simple shear. Here, the thermodynamic theory of crystal plasticity thus provides a thermodynamically consistent and physically rigorous description of dislocation motion in crystals. We also discuss new insights about the interaction of dislocations belonging to different slip systems.

36 MATERIALS SCIENCE↗

Three-dimensional explicit finite element formulation for shear localization with global tracking of embedded weak discontinuities

Shear localization is a damage mechanism characterized by the concentration of plastic deformation within thin bands of materials, often observed as a precursor to failure in metals. In this study, we present a three-dimensional explicit finite element formulation with embedded weak discontinuities for the treatment of shear localization under dynamic loading conditions. In this formulation, the onset of localization is detected via a material stability analysis suitable for rate-sensitive materials, and the continuity of propagating shear bands is ensured using a global tracking strategy involving a heat-conduction type boundary value problem, solved for a scalar level set function over the global domain of the problem. In addition, we propose a modified quadrature rule that is used to compute the contribution of an individual element to the global finite element arrays, taking into account the position of the embedded shear band within that element.

42 ENGINEERING↗

Strain localization and dynamic recrystallization in polycrystalline metals: Thermodynamic theory and simulation framework

We describe a theoretical and computational framework for adiabatic shear banding (ASB) and dynamic recrystallization (DRX) in polycrystalline materials. The Langer-Bouchbinder-Lookman (LBL) thermodynamic theory of polycrystalline plasticity, which we recently reformulated to describe DRX via the inclusion of the grain boundary density or the grain size as an internal state variable, provides a convenient and self-consistent way to represent the viscoplastic and thermal behavior of the material, with minimal ad-hoc assumptions regarding the initiation of yielding or onset of shear banding. Here, we implement the LBL-DRX theory in conjunction with a finite-element computational framework. Favorable comparison to experimental measurements on a top-hat AISI 316L stainless steel sample compressed with a split-Hopkinson pressure bar suggests the accuracy and usefulness of the LBL-DRX framework, and demonstrates the crucial role of DRX in strain localization.

36 MATERIALS SCIENCE↗

Structural representation of additively manufactured 316L austenitic stainless steel

Three 316L stainless steel materials are studied and reported upon; wrought, as-built additively manufactured (AM), and heat-treated AM material. The AM material was produced from the laser engineered net shaping (LENS) process. Macroscopic uniaxial compression stress-strain curves were obtained for all three materials. The curves were similar for the wrought and heat-treated AM materials but the as-built AM material demonstrated approximately 1.7 times greater flow stress at any given level of strain than the other two materials. Electron-backscatter diffraction analysis of the materials also showed that the microstructures of the three materials differed; with complex grain morphology for the as-built AM material. The mean grain size of each of the three materials also differed. Furthermore, the initial dislocation density was also measured with neutron diffraction and line-profile analysis for both wrought and as-built AM materials with the density in the AM material approximately 2.5 times greater.

36 MATERIALS SCIENCE↗

Slip transmission of high angle grain boundaries in body-centered cubic metals: Micropillar compression of pure Ta single and bi-crystals

In this study, we seek to probe and understand the mechanical behavior of grain boundaries in pure Ta as a model bcc material using micropillar compression experiments. Three high angle grain boundaries are chosen with varying crystal orientations. Multiple bi-crystal pillars are prepared containing a single, nearly vertical grain boundary in the approximate center of the pillar and compared against their single crystalline pillar counterparts. The main phenomenon of interest was slip transmission or strain transfer in the bi-crystals which was considered to occur when slip traces aligned across the grain boundary. This occurred in two of the three bi-crystals.

36 MATERIALS SCIENCE↗

A computationally efficient ductile damage model accounting for nucleation and micro-inertia at high triaxialities

The computational formulation of a micro-mechanical material model for the dynamic failure of ductile metals is presented in this paper. The statistical nature of porosity initiation is accounted for by introducing an arbitrary probability density function which describes the pores nucleation pressures. Each micropore within the representative volume element is modeled as a thick spherical shell made of plastically incompressible material. The treatment of porosity by a distribution of thick-walled spheres also allows for the inclusion of micro-inertia effects under conditions of shock and dynamic loading. The second order ordinary differential equation governing the microscopic porosity evolution is solved with a robust implicit procedure. A new Chebyshev collocation method is employed to approximate the porosity distribution and remapping is used to optimize memory usage. The adaptive approximation of the porosity distribution leads to a reduction of computational time and memory usage of up to two orders of magnitude. Moreover, the proposed model affords consistent performance: changing the nucleation pressure probability density function and/or the applied strain rate does not reduce accuracy or computational efficiency of the material model. The numerical performance of the model and algorithms presented is tested against three problems for high density tantalum: single void, one-dimensional uniaxial strain, and two-dimensional plate impact. Here, the results using the integration and algorithmic advances suggest a significant improvement in computational efficiency and accuracy over previous treatments for dynamic loading conditions.

42 ENGINEERING↗

Path dependency of phase transformations

At high pressures titanium and zirconium are known to undergo a phase transformation from the hexagonal close packed (HCP), alpha-phase to the simple-hexagonal, omega-phase. Under conditions of shock loading, the high-pressure omega-phase can be retained upon release. It has been shown that temperature, peak shock stress, and texture can influence the transformation. Moreover, under these same loading conditions, plastic processes of slip and twinning are also affected by similar differences in the loading path. To understand this path dependency and it role on microstructural evolution, in-situ velocimetry measurements along with post-mortem metallographic and neutron diffraction characterization of soft recovered specimens have been utilized to qualitatively understand the kinetics of transformation, quantify volume fraction of retained omega-phase and characterize the shocked alpha and omega-phases. Together the work described here can be utilized to understand phase and structure for these metals and lend insight into the partitioning of plastic processes between phases far from equilibrium conditions. Experiments discussed here, reveal that the transformation from alpha to omega phase is more difficult in Ti than in Zr metals. This result is discussed in terms of both differences in shear modulus and free space within the lattice for these systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Experimental Physical Sciences Vistas Performance through Science Winter 2017

This issue of Experimental Physical Sciences Vistas focuses on the integrated science that plays a critical role in Los Alamos National Laboratory’s support of the nation’s nuclear deterrent. I hope you will enjoy reading about these accomplishments, opportunities, and challenges.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

3D Microstructures for Materials and Damage Models

Many challenges exist with regard to understanding and representing complex physical processes involved with ductile damage and failure in polycrystalline metallic materials. Currently, the ability to accurately predict the macroscale ductile damage and failure response of metallic materials is lacking. Research at Los Alamos National Laboratory (LANL) is aimed at building a coupled experimental and computational methodology that supports the development of predictive damage capabilities by: capturing real distributions of microstructural features from real material and implementing them as digitally generated microstructures in damage model development; and, distilling structure-property information to link microstructural details to damage evolution under a multitude of loading states.

36 MATERIALS SCIENCE↗

Modeling and simulation framework for dynamic strain localization in elasto-viscoplastic metallic materials subject to large deformations

This study describes a theoretical and computational framework for the treatment of adiabatic shear band formation in rate-sensitive polycrystalline metallic materials. From a computational perspective, accurate representation of strain localization behavior has been a long-standing challenge. In addition, the underlying physical mechanisms leading to the localization of plastic deformation are still not fully understood. The proposed framework is built around an enhanced-strain finite element formulation, designed to alleviate numerical pathologies known to arise in localization problems, by allowing a localization band of given finite width (weak discontinuity) to be embedded within individual elements. The mechanical threshold strength (MTS) model is used to represent the temperature and strain rate-dependent viscoplastic response of the material. This classical flow stress model employs an internal state variable to quantify the effect of dislocation structure evolution (work hardening and recovery). In light of growing evidence suggesting that the softening effect of dynamic recrystallization may play a significant role, alongside thermal softening, in the process of shear band formation and growth, a simple dynamic recrystallization model is proposed and cast within the context of the MTS model with the aid of the aforementioned internal state variable. An initiation criterion for shear localization in rate and temperature-sensitive materials is introduced and used in the present context of high-rate loading, where material rate-dependence is pronounced and substantial temperature increases are achieved due to the dissipative nature of viscoplastic processes. In addition, explicit time integration is adopted to facilitate treatment of the dynamic problems under consideration, where strain rates in excess of 10 4 s –1 are typically attained. Two series of experiments are conducted on AISI 316L stainless steel, employing the commonly used top-hat sample geometry and the Split-Hopkinson Pressure Bar dynamic test system. Axi-symmetric finite element simulation results are compared to cross-sectional micrographs of recovered samples and experimental load–displacement results, in order to examine the performance of the proposed framework and demonstrate its effectiveness in treating the initiation and growth of adiabatic shear banding in dynamically loaded metallic materials. These comparisons demonstrate that thermal softening alone is insufficient to induce shear localization behaviors observed in some materials, such as stainless steel, and support the hypothesis that dynamic recrystallization and/or other softening mechanisms play an essential role in this process.

36 MATERIALS SCIENCE↗

A High-Rate, Single-Crystal Model including Phase Transformations, Plastic Slip, and Twinning

An anisotropic, rate-­dependent, single-­crystal approach for modeling materials under the conditions of high strain rates and pressures is provided. The model includes the effects of large deformations, nonlinear elasticity, phase transformations, and plastic slip and twinning. It is envisioned that the model may be used to examine these coupled effects on the local deformation of materials that are subjected to ballistic impact or explosive loading. The model is formulated using a multiplicative decomposition of the deformation gradient. A plate impact experiment on a multi-­crystal sample of titanium was conducted. The particle velocities at the back surface of three crystal orientations relative to the direction of impact were measured. Molecular dynamics simulations were conducted to investigate the details of the high-­rate deformation and pursue issues related to the phase transformation for titanium. Simulations using the single crystal model were conducted and compared to the high-­rate experimental data for the impact loaded single crystals. The model was found to capture the features of the experiments.

36 MATERIALS SCIENCE↗

Capturing the Complexity of Additively Manufactured Microstructures

The underlying mechanisms and kinetics controlling damage nucleation and growth as a function of material microstructure and loading paths are discussed. These experiments indicate that structural features such as grain boundaries, grain size distribution, grain morphology crystallographic texture are all factors that influence mechanical behavior.

36 MATERIALS SCIENCE↗