Effect of initial damage variability on hot-spot nucleation in energetic materials

Abstract

Mechanical insult may be able to produce chemical transformations in solids when the energy is released in highly localized regions. This phenomenon is responsible for the nucleation of hot-spots that are responsible for ignition of energetic materials. The concentration of energy at microstructural defects leads to the probabilistic nature of ignition. The effect of the microstructure of the energetic particles, specifically the influence of the initial crack distribution on the sensitivity to ignition, is studied for a particle embedded in a polymeric matrix at impact velocities 100 m/s and 400 m/s with finite element simulations that couple fracture dynamics and heat transport. A phase field damage model that includes heat sources due to frictional heating at the crack surfaces and heat dissipation during crack propagation is developed and verified. These heat sources are compared and, in the range of impact velocities studied, heat generation due to friction is more important than dissipation due to crack propagation. Hot-spots nucleated at 100 m/s do not reach the critical temperature while conditions consistent with the Lee-Tarver criterion for ignition are observed at 400 m/s impact velocity. The variability observed due to the stochasticity of the initial crack distribution is studied and it increases with a higher impact velocity. In particular, regions of high temperature develop close to cracks intersecting the particle polymer interface. Therefore, controlling the surface quality of the energetic particles may lead to a reduction on the sensitivity uncertainty in polymer-bonded explosives.

Void Collapse in Shocked ‐HMX Single Crystals: Simulations and Experiments

“…The conservation of energy equation including heat flux and heat generation is given by : …”

“…The conservation of energy equation including heat flux and heat generation is given by [7], [26]: where r is a heat source, q 0 ¼ À kr 0 T for assumed Fourier's law, k is the thermal conductivity, E : e is the elastic part of the Lagrange strain rate tensor. The term F mech is the dissipation and it is calculated from the next expression:…”

Comparison of peridynamic and phase-field models for dynamic brittle fracture in glassy materials

“…Notice that this computational phenomenon (unstoppable "thickening" of damage around crack lines) can be observed in many phase-field simulations of dynamic brittle fracture published in the literature (e.g. [90,91], see Fig. 18) even if the topic has not been paid the attention it probably deserves.…”

“…In (a): time evolution of damage field ( from[90]) with white regions in the damage field represent values with ϕ > 0.97. In (b): increasing expansion of damage in the direction perpendicular to the crack growth direction, with time (from[91]) .…”

Plastic Dissipation and Hot‐Spot Formation in HMX‐Based PBXs Subjected to Low Velocity Impact

“…At low impact speed and for isotropic plastic behavior of the crystal, the plasticity developing around the porosities induces a small temperature rise insufficient for initiation [4,16].…”