Numerical Simulation of Rock Blasting Phenomena in Small Scale Using SPH

"A numerical procedure based on Smoothed Particle Hydrodynamics (SPH) is developed for rock blasting phenomena involving multiple intersecting joints. In this procedure, an elasto-plastic damage constitutive model has been applied in the framework of SPH to simulate yielding response of the surrounding rock medium. In this paper, a square shape rock medium with centrally located emulsion explosive in a borehole is numerically blasted with and without the presence of two diagonal rock joints. The model is simulated by assuming all four sides of the rock medium as free surfaces. The paper elaborates the propagation of stress wave, expansion dynamics of gaseous products, and development of fracture, fragmentation of the medium and dissipation of energy due to damage.INTRODUCTIONDamage prediction of rock mass and its mean fragmentation size are of important to design blast parameters for a given explosive. Current blasting practice solely depends on empirical relationships between fragmentation sizes to those blast design parameters while mathematical descriptions of the physical laws governing the process of rock blasting are limited. A constitutive relationship based on properties of explosives and rock mass is of special interest because of its capabilities in describing the key features of rock fracture and fragmentation, resulting in understanding of the physical process involved in rock blasting. It is well known that blast-induced damage and fragmentation of rock mass occur due to detonation induced stress wave and product gas driven fracture propagation. Nevertheless, the fundamental mechanism of the failure process and the behavior of rock characteristics under these two different loading are not well understood in many cases, which are closely related to growth of fractures and interaction of cracks under different loading conditions.Continuum damage mechanics is considered to be appropriate to describe the failure process under such high-amplitude stress wave. Based on this approach, several theoretical models have been developed to analyze the dynamic damage of brittle material such as rock and concrete with microstructure (Chen, 1999; Taylor et al., 1986; Wu et al., 2004; Yazdchi et al., 1996). These models are mostly based on the theory of damage mechanics or fracture mechanics with some statistical treatment to incorporate random distribution of microcracks. In the application of fracture mechanics theory, sometimes it is difficult to estimate the micromechanical parameters such as crack density and fracture toughness in different modes. In this study, damage in rock material is evaluated according to the earlier model developed by Grady and Kipp (1980), which is a fracture model coupled with a material description for stress wave propagation to predict quantitatively fracture and fragmentation under explosive loading conditions."