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]

In this article, we trained a machine learning (ML) model to connect microstructural details of an energetic material formulation to its performance for the purpose of guiding the discovery of new explosive formulations. Our hypothesis was that the algorithm would robustly learn the training data and produce an accurate surrogate model. Specifically, the algorithm learned the relationship between details of the void size distribution (VSD), initiating shock pressure, and the energetic material performance. We used realistic constraints on the VSD and a range of cases were ingested by a physically informed reactive flow model working within a hydrodynamic solver running on high-performance computing resources. The ML algorithm produced a surrogate model that accurately predicted known test points around the parameter space. In addition to the utility of the model and the process used for its development, we noted interesting comparisons between what we, the authors—subject matter experts, would heuristically conclude from the training data and the surrogate model predictions. We detected nuanced details that were missed by the surrogate model; however, these details are not important to an energetic material formulator. We concluded that the algorithm did indeed robustly learn the training data and produce an accurate surrogate model. We further concluded that the surrogate model is a powerful tool to guide the formulator in the absence of subject matter experts and limited-access computing resources. [ABSTRACT FROM AUTHOR]

Heterogeneous energetic materials (EMs) contain microstructural defects such as voids, cracks, interfaces, and delaminated zones. Under shock loading, these defects offer potential sites for energy localization, i.e., hotspot formation. In a porous EM, the collapse of one void can generate propagating blast waves and hotspots that can influence the hotspot phenomena at neighboring voids. Such void–void interactions must be accounted for in predictive multi-scale models for the reactive response of a porous EM. To infuse such meso-scale phenomena into a multi-scale framework, a meso-informed ignition and growth model (MES-IG) has been developed, where the influence of void–void interactions is incorporated into the overall reaction rate through a function, f v − v . Previously, MES-IG was applied to predict the sensitivity and reactive response of EM, where f v − v was assumed to be a function of the overall sample porosity alone. This paper performs a deeper analysis to model the strong dependency of f v − v on other factors, such as void size and shock strength. The improved model for void–void interactions produces good agreement with direct numerical simulations of the HE microstructures and, thus, advances the predictive capability of multi-scale models of the shock response and sensitivity of EM. [ABSTRACT FROM AUTHOR]

Accurate simulations of the shock response of heterogeneous energetic (HE) materials require closure models, which account for energy localization in the micro-structure. In a multi-scale framework, closure is provided by reaction rate models that account for ignition and growth of hotspots, allowing for prediction of the overall macro-scale sensitivity of a HE material. In the present meso-informed ignition and growth (MES-IG) model, the reaction rate is expressed as a function of shock pressure and morphology of the void field in a pressed energetic material. In MES-IG, the void morphology is quantified in terms of a limited number of parameters: viz., overall porosity, void size, and shape (aspect ratio and orientation). In this paper, we quantify the effects of arbitrary variations in void shapes on meso-scale energy deposition rates. A collection of voids of arbitrary shapes is extracted from scanning electron microscope (SEM) images of real, pressed HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine) samples and classified into groups based on their similarity in shapes. Direct numerical simulations (DNS) are performed on the highly contorted "real" void shapes, and the calculated hotspot ignition and growth rates are compared with values predicted by the MES-IG. It is found that while the parameterization of complex void morphologies in terms of orientation and aspect ratio gives fairly good agreement between DNS and MES-IG reaction rates, the intricate details of highly complex void shapes impact hotspot characteristics to a significant extent. This work suggests possible improvements for the prediction of reaction rate in the energetic microstructure by adopting a more detailed description of shapes. [ABSTRACT FROM AUTHOR]

Meso-scale simulations of pressed energetic materials are performed using synthetic microstructures generated using deep feature representation, a deep convolutional neural network-based approach. Synthetic microstructures are shown to mimic real microstructures in the statistical representation of global and local features of micro-morphology for three different classes of pressed HMX with distinctive micro-structural characteristics. Direct numerical simulations of shock-loaded synthetic microstructures are performed to calculate the meso-scale reaction rates. For all three classes, the synthetic microstructures capture the effect of morphological uncertainties of real microstructures on the response to shock loading. The calculated reaction rates for different classes also compare well with those of the corresponding real microstructures. Thus, the article demonstrates that machine-generated ensembles of synthetic microstructures can be employed to derive structure–property–performance linkages of a wide class of real pressed energetic materials. The ability to manipulate the synthetic microstructures using deep learning-based approaches then provides an opportunity for material designers to develop and manufacture pressed energetic materials that can yield targeted performance. [ABSTRACT FROM AUTHOR]