A bar explosion model based on continuous–discontinuous element method and its application in blasting excavation

Abstract

Blasting excavation is widely used in engineering, often involving complex whole layouts. However, the small size of the blasthole and the large size of the three-dimensional (3D) numerical model lead to the large calculation scale of the 3D blasting numerical simulation, which requires considerable calculation time. Typically, a 3D numerical model is simplified into a two-dimensional (2D) numerical model, and a 1/2- or 1/4-scale model can be adopted to reduce the calculation scale. To solve this problem, a one-dimensional bar explosion model is adopted to replace the traditional solid explosion model under the framework of the continuous–discontinuous element method. The detonation pressure is directly distributed to the elements penetrated by the bar, and the volume expansion of the elements is used to calculate the volume expansion attenuation detonation pressure at each stage of detonation, thus avoiding the problem of local mesh refinement. Compared with the solid explosion case, the reliability of the bar explosion model is verified by the propagation of the explosion stress wave, peak explosion pressure, and damage nephogram. In combination with the engineering background, three blasting conditions are simulated, and the optimal one is evaluated based on fracture degree and blast fragment size pass rate.

Keywords

Bar explosion model / Continuous–discontinuous element method (CDEM) / Blasting excavation / Numerical simulation

References

[1] Asadi, P., Fakhimi, A., & Javad Ashrafi, M. (2023). Dynamic tensile strength of rock specimens with different defect lengths. Engineering Fracture Mechanics, 284, 109245.

[2] Durand, M., Marin, P., Faure, F., & Raffin, B . (2012). DEM-based simulation of concrete structures on GPU. European Journal of Environmental and Civil Engineering, 16(9), 1102-1114.

[3] Fan, H. F., Bergel, G. L., & Li, S. F. (2016). A hybrid peridynamics-SPH simulation of soil fragmentation by blast loads of buried explosive. International Journal of Impact Engineering, 87, 14-27.

[4] Feng, C., Li, S. H., Liu, X. Y., & Zhang, Y. N. (2014). A semi-spring and semi-edge combined contact model in CDEM and its application to analysis of Jiweishan landslide. Journal of Rock Mechanics and Geotechnical Engineering, 6(1), 26-35.

[5] Feng, C., Zhu, X., and Wang, X. (2022). GDEM-Block Dyna (v1.0 ed.). https://support.gdem-tech.com.

[6] Fukuda, D., Mohammadnejad, M., Liu, H. Y., Dehkhoda, S., Chan, A., Cho, S. H., Min, G. J., Han, H. Y., Kodama, J. I., & Fujii, Y. (2019). Development of a GPGPU-parallelized hybrid finite-discrete element method for modeling rock fracture. International Journal for Numerical and Analytical Methods in Geomechanics, 43(10), 1797-1824.

[7] Han, H. Y., Fukuda, D., Liu, H. Y., Fathi Salmi, E., Sellers, E., Liu, T. J., & Chan, A. (2020). FDEM simulation of rock damage evolution induced by contour blasting in the bench of tunnel at deep depth. Tunnelling and Underground Space Technology, 103, 103495.

[8] Henrych, J . (1979). The dynamics of explosion and its use. Elsevier.

[9] Huang, X. P., Zhu, B., & Chen, Y. M. (2025). A coupled and parallel peridynamics-SPH modeling and simulation of buried explosion induced soil fragmentation and cratering. Computers and Geotechnics, 178, 106942.

[10] Li, S. H., Feng, C., & Zhou, D. (2018). Mechanics methods in landslide research. Science Press (in Chinese).

[11] Li, S. H., Zhou, D., & Liu, T. P. (2013). Risk analysis method of accumulated landslide based on fracture degree. Chinese Journal of Rock Mechanics and Engineering, 32(S2), 3909-3917 (in Chinese).

[12] Li, Y. P., Feng, C., Ding, C. X., & Zhang, Y. M. (2022). A novel continuous-discontinuous multi-field numerical model for rock blasting. Applied Sciences, 12(21), 11123.

[13] Li, Y. P., Feng, C., & Zhang, Y. M. (2023). Numerical analysis of directional rock blasting with continuous-discontinuous element method. KSCE Journal of Civil Engineering, 27(8), 3591-3598.

[14] Li, Y. P., Li, J., Feng, C., Wen, M. J., & Zhang, Y. M. (2024). An interface constitutive model of plastic tensile-compressive damage under impact loading based on continuous-discontinuous framework. Computers and Geotechnics, 173, 106502.

[15] Liu, Q. S., Wang, W. Q., & Ma, H. (2020). Parallelized combined finite-discrete element (FDEM) procedure using multi-GPU with CUDA. International Journal for Numerical and Analytical Methods in Geomechanics, 44(2), 208-238.

[16] Lu, A., Yan, P., Lu, W. B., Li, X. F., Liu, X., Luo, S., Huang, S. L., & Grasselli, G. (2024). Crack propagation mechanism of smooth blasting holes for tunnel excavation under high in-situ stress. Engineering Fracture Mechanics, 304, 110144.

[17] Munjiza, A. (2004). Chapter 2, processing of contact interaction in the combined finite discrete element method. In The Combined Finite-Discrete Element Method (pp. 35-72). John Wiley & Sons Inc..

[18] Munjiza, A., Lei, Z., Divic, V., & Peros, B. (2013). Fracture and fragmentation of thin shells using the combined finite-discrete element method. International Journal for Numerical Methods in Engineering, 95(6), 478-498.

[19] Murali, G., Wong, L. S., & Abid, S. R. (2024). A comprehensive review of drop weight impact testing: Evaluating the pros and cons in fiber-reinforced concrete performance assessment. Journal of Building Engineering, 94, 109934.

[20] Park, D., Jeon, B., & Jeon, S. (2009). A numerical study on the screening of blast-induced waves for reducing ground vibration. Rock Mechanics and Rock Engineering, 42(3), 449-473.

[21] Saifi, F., Javaid, M., Haleem, A., & Anas, S. M. (2024). Computational fluid dynamics approach for predicting pipeline response to various blast scenarios: A numerical modeling study. Computer Modeling in Engineering and Sciences, 140(3), 2747-2777.

[22] Sharafisafa, M., Aliabadian, Z., & Shen, L. M. (2025). Numerical study of presplit blasting in rock masses with a closed and filled joint using coupled finite-discrete element method. Simulation Modelling Practice and Theory, 144, 103199.

[23] Xu, S. X., Wu, B., Zhang, H. L., Qi, S. X., Bian, H. B., & Wang, J. J. (2025). Directional crack propagation and optimization strategies for multi-hole shaped charge blasting in tunnel construction. Structures, 72, 108268.

[24] Yin, Y., Esmaeili, K., Sun, Q., & Cao, J. (2025). Numerical investigation of rock damage induced by bilateral-groove-slot shaped charge blasting under the influence of in-situ stresses. Computers and Geotechnics, 180, 107070.

[25] Zhou, H. X., Gao, Q. D., Leng, Z. D., Wang, Y. Q., Fan, Y., & Liu, G. F. (2024). Formation mechanisms of different kinds of blast-induced cracks and their extension characteristics in rock mass. Computers and Geotechnics, 176, 106747.

[26] Zhu, X. G., Feng, C., Cheng, P. D., Wang, X. Q., & Li, S. H. (2021). A novel three-dimensional hydraulic fracturing model based on continuum-discontinuum element method. Computer Methods in Applied Mechanics and Engineering, 383, 113887.