Numerical simulation of explosive fracturing with smoothed particle hydrodynamics

In this paper we study explosive fracturing with smoothed particle hydrodynamics (SPH). As a particle based Lagrangian method, SPH is particularly suited to the analysis of fracture due to its full Lagrangian frame and capacity to model large deformation. We adopt the Jones-Wilkins-Lee equation as equation of state of the trinitrotoluene (TNT) explosive and a continuum elasto-damage model to predict the fracture of the rock. We predict the evolution of damage using the strain history of each particle. To strengthen the interaction of coupling interfaces we use a penalty function to avoid penetration between different material particles