Modelling the dynamic mechanical response of rock mass under sequential explosive detonation, especially the fracture and energy evolution characteristics, is the key to optimizing the spatial distribution and initiation time interval of explosives. A hybridized continuum-discontinuum element method and energy statistics algorithm are implemented in this study to accurately investigate the rock dynamic response induced by sequential detonation loading. Based on the rock fracture status and the stress–strain curve of bilinear cohesive fracture model, an energy statistics algorithm considering all energy components is first proposed, and its accuracy and robustness are verified. Then, the continuum-discontinuum element method and energy statistics algorithm are adopted to simulate the dynamic response of rock mass under sequential detonation loading. The results indicate that the area of cracked interfaces generated by each explosive first increases and then decreases, which caused by explosive 2 is the largest, accounting for 38.19%. Due to the damage to rock caused by the previous detonation loading, the rock around explosive 1 mainly undergoes shear fracture, while the rock around explosives 2 and 3 mainly undergoes tensile fracture. The change trends of various energy components are different, and the explosion energy is mainly converted into friction energy and element kinetic energy.
Read full abstract