Document Type : Original Research Paper


1 Department of Mining Engineering, Faculty of Engineering, Tarbiat Modares University, Tehran, Iran

2 Department of Mining Engineering, Faculty of Engineering, University of Kurdistan, Sanandaj, Iran


Determining the appropriate blasting pattern is important to prevent any damage to the tunnel perimeter in conventional tunneling by blasting operation in hard rocks. In this research work, the LS-DYNA software and numerical finite element method (FEM) are used for simulation of the blasting process in the Miyaneh-Ardabil railway tunnel. For this aim, the strong explosive model and nonlinear kinematic plastic material model are considered. Furthermore, the parameters required for the Johnson-Holmquist behavioral model are based on the Johnson-Holmquist-Ceramic material model relationships and are determined for the andesitic rock mass around studied tunnel. The model geometry is designed using AUTOCAD software and Hyper-mesh software is applied for meshing simulation. After introducing elements properties and material behavioral models and applying control and output parameters in LS-PrePost software, the modeling process is performed by LS-DYNA software. Different patterns of blastholes including 66, 23, and 19 holes, with diameters of 40 and 51 mm, and depths of 3 to 3.8 m are investigated by three-dimensional FEM. The borehole pressure caused by the ammonium nitrate-fuel oil (ANFO) detonation is considered based on the Jones-Wilkins-Lee (JWL) equation of state in the LS-DYNA software. The outer boundaries of the model are considered non-reflective to prevent the wave’s return. The results showed that LS-DYNA software can efficiently simulate the blasting process. Moreover, the post-failure rate of the blasting is reduced by more than 30% using the main charge with less explosive power and reducing the distance and diameter of contour holes.


Main Subjects

[1]. Holmberg, R. & Persson, P. A. (1980). Design of Tunnel perimeter blasthole patterns to prevent rock damage. Transactions Institutions of Mining and Metallurgy, 89, 37-40.
[2]. Saiang, D. (2008). Behavior of Blast-Induced Damaged Zone around Underground Excavations in Hard Rock Mass, Department of Civil, Mining and Environmental Engineering, Luleå University of Technology, Doctoral Thesis.
[3]. Saiang, D. & Nordlund, E. (2007). Failure Mechanisms around shallow tunnels in Brittle rock. In: R.E. Sousa, Olalla and Grossman (Editors), 11th International Society of Rock Mechanics Symposium. Talyor and Francis, Lisbon, Portugal, 883-890.
[4]. Saiang, D. & Nordlund, E. (2009). Numerical Analyses of the Influence of Blast-Induced Damaged Rock Around Shallow Tunnels in Brittle Rock. Rock Mechanics and Rock Engineering, 42(3), 421-448.
[5]. Iverson, S.R., Hustrulid, W.A., Johnson, J.C., Tesarik, D. & Akbarzadeh, Y. (2009). The extent of blast damage from a fully coupled explosive charge. Proceedings of the ninth international symposium on rock fragmentation by blasting. Granada, Spain, Taylor & Francis, 459–68.
[6]. Wang, Z.L., Konietzky, H. & Shen, R.F. (2009). Coupled finite element and discrete element method for underground blast in faulted rock masses. Soil Dynamics and Earthquake Engineering, 29, 939–945.
[7]. Li, H., Xiang, X., Li, J., Zhao, J., Liu, B. & Liu, Y. (2011). Rock damage control in bedrock blasting excavation for a nuclear power plant. International Journal of Rock Mechanics and Mining Sciences48(2), 210 – 218. 
[8]. Onederra, I. A., et al. (2013). Modelling Blast Induced Damage from a Fully Coupled Explosive Charge. International Journal of Rock Mechanics and Mining Sciences, 58, 73–84.
[9]. Yang, J., Lu, W., Hu, Y. et al. (2015). Numerical Simulation of Rock Mass Damage Evolution during Deep-Buried Tunnel Excavation by Drill and Blast. Rock Mechanics Rock Engineering, 48, 2045–2059.
[10]. Yingjie, L., Dingli, Z., Qian, F., Qingchun, Y. & Lu, X. (2014). A physical and numerical investigation of the failure mechanism of weak rocks surrounding tunnels. Computers and Geotechnics, 61, 292-307.
[11]. Xie, L. X., Lu, W., Zhang, Q. B., Jiang, Q. H., Wang, G. H. & Zhao, J. (2016). Damage evolution mechanisms of rock in deep tunnels induced by cut blasting. Tunneling and Underground Space Technology58, 257-270. 
[12]. Eslami, M. & Goshtasbi, K. (2018). Blasting Damage Predictions by Numerical Modeling in Siahbishe Pumped Storage Powerhouse. Journal of The Institution of Engineers (India) Series D, 99 (8), 133–146.
[13]. Dang, V.K., Dias, D., Do, N. A. & Vo, T.H. (2021). Impact of blasting at tunnel face on an existing adjacent tunnel. GEOMATE Journal, 15(47), 22-31.
[14]. Daraei, A. & Shokrollah, Z. (2018). Prediction of overbreak depth in Ghalaje road tunnel using strength factor. International Journal of Mining Science and Technology, 28 (4), 679-684.
[15]. Verma, H.K., Samadhiya, N.K., Singh, M., Goel, R.K. & Singh, P.K. (2018). Blast induced rock mass damage around tunnels. Tunneling and Underground Space Technology, 71, 149–158.
[16]. Villalobos, S., Cacciari, P. & Futai, M. (2020). Stability assessment around a railway tunnel using terrestrial laser scanner data and finite element analysis. Revista ingeniería de construcción, 35(1), 21-33.
[17]. Zou, B., Xu, Zh., Wang, J., Luo, Zh. & Hu, L. (2020). Numerical Investigation on Influential Factors for Quality of Smooth Blasting in Rock Tunnels. Advances in Civil Engineering, Article ID 9854313.
[18]. Ali dadi, G., Ahmadi, M. & Nejati, H.R. (2021). An equation for determining the thickness of the damaged area due to unneling using numerical and statistical methods. Iranian journal of rock mechanics, 5(3), 33-43. (In Persian).
[19]. Chinaei, F., Ahangari, K. & Shirinabadi, R. (2023). Determination of Hoek and Brown damage factor due to explosion around tunnels using numerical modeling. Journal of Analytical and Numerical Methods in Mining Engineering, 13(34), 15-26.
[20]. Iran Oston Consulting Engineers, (2017). General, Structural, and Engineering Geology Report of Ardebil-Miyaneh Railway Tunnel (In Persian).
[21]. Bakhshandeh, H. & Nejat Dehkordi, F. (2016). Simulation of airblast and introducing an appropriate model according to the conditions of the study zone. 4th International Mine & Mining Industries Congress & 6th Iranian Mining Engineering Conference, Tehran.
[22]. Sanchidrián, J. A., Castedo, R., López, L. M., Segarra, P. & Santos, A. P. (2015). Determination of the JWL Constants for ANFO and Emulsion Explosives from Cylinder Test Data. Central European Journal of Energetic Materials, 12(2), 177-194.
[23]. Wang, J., Yin, Y. & Luo, C. (2018). Johnson–Holmquist-II (JH-2) Constitutive Model for Rock Materials: Parameter Determination and Application in Tunnel Smooth Blasting. Applied Sciences, 8(9), 1675.
[24]. Livermore Software Technology Corporation (LSTC), (2017). LS-DYNA Keyword user’s manual. Volume II,  R10.0, 1577 pp.
[25]. Saleh Asheghabadi, M. & Cheng, X. (2020). Analysis of Undrained Seismic Behavior of Shallow Tunnels in Soft Clay Using Nonlinear Kinematic Hardening Model. Applied Sciences, 10(8), 2834.
[26]. Yang, R., Bawden, W.F. & Katsabanis, P.D. (1996). A new constitutive model for blast damage. International Journal of Rock Mechanics and Mining Sciences, 33, 245–254.
[27]. Hallquist, John O. (2006). Ls-Dyna Theory Manual. Livermore Software Technology Corporation.
[28]. Murray, Y.D., Reid, J.D., Faller, R.K., Bielenberg, B.W. & Paulsen, T.J. (2007). Manual for LS-DYNA wood material model 143 Report No. FHWA-HRT-04-096. US Federal Highway Administration, Turner-Fairbank Highway Research Center.
[29]. Sun, Ch. (2013). Damage zone prediction for rock blasting. Ph.D. dissertation, Department of Mining Engineering, University of Utah, 225 pp.