Simulating Dynamic Failure of Polymer-Bonded Explosives under Periodic Excitation

Accidental mishandling of explosive materials leads to thousands of injuries in the US every year. Understanding the mechanisms behind the detonation process is crucial to prevent such accidents. In polymer-bonded explosives (PBX), high-frequency mechanical excitation generates thermal energy and can lead to an increase in temperature and vapor pressure, and potentially the initiation of the detonation process. However, the mechanisms behind this energy release, such as the effects of dynamic fracture and friction, are not well understood. Experimental data is difficult to collect due to the different time scales of reactions and vibrations, so research is aided by running simulations to computationally understand experimental results. Using phase-field model of fracture, we simulate the behaviors of various crack orientations in single particles of HMX bonded in a polymer matrix. Larger amplitudes induce higher rates of energy buildup which lead to quicker crack propagation, while higher frequencies generate higher spikes in temperature. However, crack location and orientation with respect to loading also significantly affect damage rates and temperature fluctuations. Cracks perpendicular to the loading vibration wave propagate most readily and appear to generate the most frictional energy, especially along the crystal-polymer interface

Computational Modeling of Energetic Materials Under Impact and Shock Compression

Understanding the fundamental physics involved in the high strain rate deformation of high explosives (HE) is critical for developing more efficient, reliable, and safer energetic materials. When HE are impacted at high velocities, several thermo-mechanical processes are activated, which are responsible for the ignition of these materials. These processes occur at different time and length scales, some of them inaccessible by experimentation. Therefore, computational modeling is an excellent alternative to study multiscale phenomena responsible for the ignition and initiation of HE. This thesis aims to develop a continuum model of β-HMX to study the anisotropic behavior of HE at the mesoscale, including fracture evolution and plastic deformation. This thesis focus on three types of simulations. First, we investigate dynamic fracture and hotspot formation in β-HMX particles embedded in Sylgard binder undergoing high strain rate compression and harmonic excitation. We use the phase field damage model (PFDM) to simulate dynamic fracture. Also, we implement a thermal model to capture temperature increase due to fracture dissipation and friction at both cracks and debonded HMX/Sylgard interface. In our simulations, we observe that crack patterns are strongly dominated by initial defects such as pre-existing cracks and interface debonding. Regions with initial debonding between HMX particles and the polymer are critical sites where cracks nucleate and propagate. Heating due to friction generates in these regions too and caused the formation of critical hotspots. We also run simulations of a β-HMX particle under high-frequency harmonic excitation. As expected, higher frequencies and larger amplitudes lead to an increase in the damage growth rate. The simulations suggest that the intensity of the thermal localization can be controlled more readily by modifying the bonding properties between the particle and the binder rather than reducing the content of bulk defects in the particle. Second, we present simulations of shock compression in β-HMX single crystals. For this purpose, we implemented a constitutive model that simulates the elastoplastic anisotropic response of this type of material. The continuum model includes a rate-dependent crystal plasticity model and the Mie-Grüneisen equation of state to obtain the pressure due to shock. Temperature evolves in the material due to plastic dissipation, shock, and thermoelastic coupling. The model is calibrated with non-reactive atomistic simulations to make sure the model obeys the Rankine-Hugoniot jump conditions. We compare finite element (FE) and molecular dynamic (MD) simulations to study the formation of hot spots during the collapse of nano-size void in a β-HMX energetic crystal. The FE simulations captured the transition from viscoelastic collapse for relatively weak shocks to a hydrodynamic regime for strong shocks. The overall temperature distributions and the rate of pore collapse are similar to MD simulations. We observe that the void collapse rate and temperature field are strongly dependent on the plasticity model, and we quantify these effects. We also studied the collapse of a micron size void in HMX impacted at different crystal orientations and impact velocities. The simulation results of void collapse are in good agreement with a gas gun void collapse experiment. While the void size and crystal orientation do not affect the area ratio rate, they strongly affect the void collapse regime and temperature