Search

High-fidelity constitutive models are critical for accurate meso-scale continuum modeling and prediction of shock initiation of crystalline energetic materials (EMs). While empirically calibrated or atomistic-guided anisotropic elastoplastic models of EM such as cyclotetramethylene-tetranitramine (HMX) and cyclotrimethylene-trinitramine (RDX) capture important micromechanical phenomena (such as dislocation evolution, slip-resistance, and anisotropic elasticity), the computational cost of using anisotropic single-crystal plasticity models can become prohibitive for meso-scale computations of void-collapse and hotspot formation in microstructures. Thermo-mechanically representative, isotropic, pressure, temperature, and rate-dependent material constitutive models are practical alternatives for meso-scale simulations of the shock response of microstructures. To this end, this work constructs physically consistent isotropic plasticity from anisotropic single-crystal plasticity models for HMX and RDX. State-of-the-art crystal plasticity models for HMX and RDX are used to compute the stress states in single crystals oriented in three different directions relative to shocks generated by impact at velocities ranging from 100 to 1000 m/s. Post-shock von Mises stress fields for the three orientations are then used to calibrate the strain-rate hardening coefficient and the reference strain rate for a rate-dependent Johnson–Cook (JC) yield surface model. We compare the pressures and the post-shock von Mises stresses between the JC and the anisotropic models to show that the isotropic computations closely approximate the averaged deformation response of the three different crystal orientations. We then model the interaction of a shock generated by a 500 m/s impact with a 0.5 μm void and show that the pressures and the deviatoric stresses obtained using the isotropic model closely match those computed from anisotropic models for both HMX and RDX. The resulting isotropic J2 plastic flow model for HMX and RDX can be employed to perform meso-scale simulations for energy localization due to shear bands and void collapse in the two materials. [ABSTRACT FROM AUTHOR]

This paper demonstrates a novel Eulerian computational framework for modeling anisotropic elastoplastic deformations of organic crystalline energetic materials (EM) under shocks. While Eulerian formulations are advantageous for handling large deformations, constitutive laws in such formulations have been limited to isotropic elastoplastic models, which may not fully capture the shock response of crystalline EM. The present Eulerian framework for high-strain rates, large deformation material dynamics of EM incorporates anisotropic isochoric elasticity via a hypo-elastic constitutive law and visco-plastic single-crystal models. The calculations are validated against atomistic calculations and experimental data and benchmarked against Lagrangian (finite element) crystal plasticity computations for shock-propagation in a monoclinic organic crystal, octahydro-1,3,5,7-tetranitro-1,3,5,7 tetrazocine (β-HMX). The Cauchy stress components and the resolved shear stresses calculated using the present Eulerian approach are shown to be in good agreement with the Lagrangian computations for different crystal orientations. The Eulerian framework is then used for computations of shock-induced inert void collapse in β-HMX to study the effects of crystal orientations on hotspot formation under different loading intensities. The computations show that the hotspot temperature distributions and the collapse profiles are sensitive to the crystal orientations at lower impact velocities (viz., 500 m/s); when the impact velocity is increased to 1000 m/s, the collapse is predominantly hydrodynamic and the role of anisotropy is modest. The present methodology will be useful to simulate energy localization in shocked porous energetic material microstructures and other situations where large deformations of single and polycrystals govern the thermomechanical response. [ABSTRACT FROM AUTHOR]

The sensitivity of an energetic material is strongly influenced by its microstructure. This work distinguishes the roles played by the microstructure (i.e., the meso-scale) in the macro-scale shock sensitivity of two different materials: TATB and HMX. To quantify sensitivity, we develop a meso-informed energy deposition model for a porous TATB material, following procedures from the previous work on HMX. Simulations of reactive void collapse in TATB are employed to calculate the rate of initiation and growth of hotspots. These rates are expressed as surrogate models, expressing meso-scale (hotspot) quantities of interest as functions of shock strength P s and void size D v o i d . The hotspot ignition and growth rate surrogates for TATB are compared with those for HMX, providing insights into meso-scale physics underlying shock sensitivity of these two energetic materials. The surrogate models are then used in a meso-informed ignition and growth (MES-IG) model to close macro-scale simulations of the shock response of porous TATB. We also obtain the run-to-detonation distances and generate Pop-plots to quantify macro-scale sensitivity. It is shown that Pop-plots for HMX-based energetic materials accord with behavior observed in experimental studies; however, there is a significant discrepancy between MES-IG predictions and experiments for TATB; the causes for this difference between HMX and TATB are discussed, pointing to areas for future work. [ABSTRACT FROM AUTHOR]

Predictive simulations of shock-to-detonation transitions (SDTs) of energetic materials must contend with uncertainties in the material properties, reactive models, and the microstructures of the material. In this work, we analyze the effects of uncertainties in the run-to-detonation distance h of a pressed energetic (HMX) material due to variabilities in the thermomechanical properties of HMX. The run distances are computed using a recently developed machine-learning based multiscale modeling framework, viz., the Meso-informed Ignition and Growth (MES-IG) model. The input uncertainties are first used in the MES-IG model to quantify the variabilities in the hotspot dynamics at the mesoscale. A Kriging-based Monte Carlo method is used to construct probability density functions (pdfs) for the mesoscale reaction-product formation rates; these are used to propagate the mesoscale uncertainties to the macroscale reaction-progress variables to construct pdfs for the run-to-detonation distance h. We evaluate uncertainties in h due to variabilities in six material properties, viz., specific heat, Grüneisen parameter, bulk modulus, yield strength, thermal expansion coefficient, and the thermal conductivity of the material. Among these six properties, h is found to be most sensitive to the variabilities in the specific heat of the material; the uncertainties in the specific heat amplify exponentially across scales and result in logarithmic pdfs for h. Thus, the paper not only quantifies and propagates uncertainties in material properties across scales in a multiscale model of SDT, but also ranks the properties with respect to the sensitivity of the SDT response of heterogeneous energetic materials on each property. [ABSTRACT FROM AUTHOR]