Numerical analysis of blast-induced wave propagation and spalling damage in a rock plate

The paper focuses on the numerical analysis of blast-induced stress wave propagation and related spalling damage in a rock plate or wall. Firstly, a large-scale contact explosion is simplified as a one-dimensional (1-D) strain problem, and a cross-format centered finite-difference scheme is presented. Secondly, the finite-difference code is used to explore the wave propagation in a rock plate, and compared with the commercial software LS-DYNA. Thirdly, a continuum damage constitutive model for rock is introduced and successfully incorporated into the LS-DYNA through its user-defined subroutine. Coupling with the erosion technique, the improved LS-DYNA is used to simulate the blast-induced spalling damage at the back-side of rock plate. Numerical results show that the 1-D finite difference code is convenient and accurate. Also, the user-defined subroutine can well describe the process of spalling damage.