Rear-Surface Collapse of Finite Thickness Concrete Targets under Internal Explosion

An experimental investigation on the buried internal explosion in finite thickness concrete targets was carried out, with the aim at developing an available criterion for the critical collapse of rear-surface to determine the critical collapse thickness and the critical amount of explosive charge under different depth of buried. It is found, under a certain density and diameter of explosive charge, the critical collapse thickness increases monotonically with the length-to-diameter ratio or the amount of the explosive charge, but the increasing becomes slower down after the length-to-diameter ratio of the explosive charge is larger than about 5, which implies that the geometry of the explosive charge can have much influence on the damage and failure of concrete targets due to different mechanism of energy dissipation. Moreover, by using the dimensional analysis approach, the function relation between the dimensionless critical collapse thickness and the length-to-diameter ratio was obtained, which shows that the dimensionless critical collapse thickness depends on both the amount and the length-to-diameter ratio of the charge.


Keywords:    Internal explosioncritical collapse thicknessdamage and failure spallingconcrete targetlength-to-diameter ratio

Concrete is one of widely used structural material in the civil and military protective engineering. It is associated with the penetration of a missile into and its explosion inside a concrete target, a compressive stress wave is formed and propagates to the rear free surface. When the compressive wave arrived at the rear free surface, a reflected tensile wave comes subsequently into being, which can result in spalling and collapse of the rear free surface under certain critical conditions. The failure fragments with high kinetic energy can endanger considerably the equipment as well as the people inside the protective structure, which makes it of great significance to understand in depth such a rear-surface spalling and collapse phenomena of a finite thickness concrete target under internal explosions and gets more attentions in the realm of protection of the structures.


Over the past decades, much work has been done on the damage and failure of concrete targets under explosion and impact loadings. Fu1, et al., Wang2, et al., and Wu3, et al. investigated the shallow buried explosion in semi-infinite concrete bodies, and the relationship between the geometrical properties of blasting crater and the amount of spherical explosive charge was determined; Xu4, et al., Zhou5,6, et al., Lu7, et al., Tai8, et al., Wang9, et al., and Yi10, et al. studied the dynamic response of concrete thin slabs subjected to the shock wave loadings generated by explosions in air. Rabczuk and Eibl11, Ohtsu12, et al., and Yuan13, et al. concerned mainly with the blast-induced spalling damage at the rear free surface of concrete thin slabs under contact explosion on the front surface of targets. It is worthy of special mention that a spherical explosive charge assumption, even for a cylindrical one, was always made for such a problem in most of previous literature, by which the influence of the geometrical properties of the explosive charge on the material damage and failure are neglected. Liu14, et al. has pointed out experimentally, the length-to-diameter ratio of a cylindrical explosive charge has much effect on the damage and failure of concrete targets. Therefore, in the authors’ opinion, the spherical explosive charge assumption should not be available for the charge in a large length-to-diameter ratio, especially for modern weapons, in which the maximum length-to-diameter ratio of explosive charge can reach high up to 10. However, so far, the rear-surface damage and failure of a finite thickness concrete target under internal explosion is not yet very well understood, and especially less experimental work has addressed such a problem.


In this paper, the rear-surface collapse of a finite thickness concrete target under internal explosions is investigated experimentally, and the geometrical property of explosive charge is taken into account. The amount of explosive charge at the critical collapse state is obtained, and then the relationship between the critical collapse thickness and the critical amount of a cylindrical explosive charge is proposed by using the dimensional analysis approach. The results can serve as a reference for damage analysis of the deep earth-penetrating weapons and the design of deeply buried protective structures.



The concrete used in the study has a uniaxial compressive strength of 35 MPa measured after standard curing of 28-day and the density of 2.5×103 kg/m3. The dimension of each concrete target is 1000 mm × 1000 mm × 700 mm or 1500 mm × 1500 mm × 900 mm. For convenience, as is shown in Fig. 1, several targets are casted into a monolithic one, in which some wood plates are inserted to block off each target as well as to prevent possible crack growth across the boundary between two neighbour targets. To avoid the lateral boundary effect, subsidiary concrete layers of 500 mm in thickness are casted around the monolithic target. A pre-drilled charge hole with its diameter less than one thirtieth of the minimum dimension of the target is located at the center of each of the target. Totally, eight kinds of the charge hole in different depth are pre-drilled for eight different depths of burial of 350 mm, 370 mm, 400 mm, 450 mm, 490 mm, 500 mm, 520 mm, and 525 mm, respectively.


Figure 1.Concrete targets



The casted cylindrical TNT explosive charge is used in all the experiments, with its diameter D = 30 mm and density ρe = 1.58×103 kg/m3, respectively. The charge amount is then determined by the length of the explosive charge. An electrical detonator is fixed at one end of the explosive charge, as shown in Fig. 2, and the positive initiation is adopted.


Figure 2. TNT charge.


It is easy to understand, for an internal explosion under a certain depth of burial, that the rear scab thickness or the depth of the blasting crater increases with the amount of explosive charge, and hence a critical collapse state can be defined as that when the depth of the blasting crater R equals to the residual length of the finite target Rc, as shown in Fig. 3. On the other hand, the residual length Rc can also be taken as a critical collapse thickness for an internal explosion with a given depth of burial, and the corresponding amount of explosive charge Mc can be defined naturally as the critical amount of explosive charge. In the experiments, for a given depth of burial and thus a fixed Rc, the critical collapse state can be acquired through increasing step-by-step the amount of explosive charge, and the critical amount of explosive charge can thus be determined as that at the critical collapse state. Moreover, the charged pre-drilled hole is not plugged up for simulating better the internal explosion of the modern weapons.


Figure 3. TNT charge in a concrete target.



Twenty-six concrete targets are used in the experimental study, and twenty-four tests are completed. Experimental data are listed in Table 1, together with the qualitative descriptions of the rear-surface failure of the concrete targets.


Table 1.Experimental data of concrete targets


A postmortem observation on typical targets of 700 mm in thickness and the critical collapse thickness of 250 mm is presented as shown in Figs. 4(a)-4(d), corresponding to the amount of explosive charge 50 g, 80 g, 100 g, and 110 g. There are no blasting craters in the first two tests, and the depth of the blasting crater for the last two tests is 140 mm, and 310 mm, respectively. In other words, there is R100 < Rc < R110 (the numerical subscripts denote the amount of explosive charge), and thus Mc = 105 g is taken as the critical amount of explosive charge considered the minimum increment of explosive charge of 10 g in the experiment. The critical collapse states corresponding to eight different depths of burial can be determined, as listed in Table 2.





Figure 4. Postmortem observation on the targets with the depth of burial h = 450 mm.


Table 2. Experimental data at the critical collapse states


From the experimental results, the influence of the length-to-diameter ratio on the rear-surface collapse can be estimated. To describe quantitatively the relationship between the critical collapse thickness and the length-to-diameter ratio, the dimensional analysis approach is taken. For such an internal explosion issue, the governing parameters can be divided into two classes that describe explosive charge and concrete, respectively, which are listed in Table 3.


Table 3.Governing parameters in the internal explosion of finite concrete targets


Taken the critical collapse thickness Rc as the quantity being determined, we have



R c =f( ρ e ,γ,L,D,h; ρ c , σ t , σ c ,E,μ ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqipC0xf9vqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbqaq=laadkfada WgaaWcbaGaam4yaaqabaGccqGH9aqpcaWGMbWaaeWaaeaacqaHbpGC daWgaaWcbaGaamyzaaqabaGccaGGSaGaeq4SdCMaaiilaiaadYeaca GGSaGaamiraiaacYcacaWGObGaai4oaiabeg8aYnaaBaaaleaacaWG JbaabeaakiaacYcacqaHdpWCdaWgaaWcbaGaamiDaaqabaGccaGGSa Gaeq4Wdm3aaSbaaSqaaiaadogaaeqaaOGaaiilaiaadweacaGGSaGa eqiVd0gacaGLOaGaayzkaaaaaa@5507@         (1)

It is easy to verify the parameters ρe, D and E as dimensional independent, and according to the Buckingham Π theorem Eqn.(1) can then be transformed into following dimensionless form:



R c D =f( γ, L D , h D ; ρ c ρ e , σ t E , σ c E ,μ ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqipC0xf9vqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaWaaSaaaeaaca WGsbWaaSbaaSqaaiaadogaaeqaaaGcbaGaamiraaaacqGH9aqpcaWG MbWaaeWaaeaacqaHZoWzcaGGSaWaaSaaaeaacaWGmbaabaGaamiraa aacaGGSaWaaSaaaeaacaWGObaabaGaamiraaaacaGG7aWaaSaaaeaa cqaHbpGCdaWgaaWcbaGaam4yaaqabaaakeaacqaHbpGCdaWgaaWcba GaamyzaaqabaaaaOGaaiilamaalaaabaGaeq4Wdm3aaSbaaSqaaiaa dshaaeqaaaGcbaGaamyraaaacaGGSaWaaSaaaeaacqaHdpWCdaWgaa WcbaGaam4yaaqabaaakeaacaWGfbaaaiaacYcacqaH8oqBaiaawIca caGLPaaaaaa@5470@        (2)

In fact, damage and failure at the rear-surface of concrete targets is almost invariable with the same amount and geometrical properties of explosive charge when the depth of burial of explosive charge is large enough, which is satisfied in all these experiments. Dimensionless parameter h/D or the depth of burial of explosive charge h could not be taken into account. Moreover, for a given explosive charge and a target material, all the material parameters as well as corresponding dimensionless parameters are constant. Therefore, Eqn.(2) becomes:


R c D =f( L D ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqipC0xf9vqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaWaaSaaaeaaca WGsbWaaSbaaSqaaiaadogaaeqaaaGcbaGaamiraaaacqGH9aqpcaWG MbWaaeWaaeaadaWcaaqaaiaadYeaaeaacaWGebaaaaGaayjkaiaawM caaaaa@3DC9@        (3)

Combining with the experimental results presented above, the function f in Eqn.(3) can be determined by using the cubic polynomial fitting technique, in which the constant is too small and can be neglected. At the same time, for clarity, the fitted curve as well as the experimental results is depicted in Fig. 5, and the maximum relative error defined by the ratio of the absolute value of the difference between the experimental data and the fitted results to the experimental data is 6.0 per cent as shown in Table 4.


Figure 5. Relationship between Rc/D and L/D.


Table 4. Comparison of fitting results and experimental data



R c D =3.993 L D [ 10.107( L D )+0.004 ( L D ) 2 ] MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqipC0xf9vqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGacaGaaiaabeqaamaabaabaaGcbaWaaSaaaeaaca WGsbWaaSbaaSqaaiaadogaaeqaaaGcbaGaamiraaaacqGH9aqpcaaI ZaGaaiOlaiaaiMdacaaI5aGaaG4mamaalaaabaGaamitaaqaaiaads eaaaWaamWaaeaacaaIXaGaeyOeI0IaaGimaiaac6cacaaIXaGaaGim aiaaiEdadaqadaqaamaalaaabaGaamitaaqaaiaadseaaaaacaGLOa GaayzkaaGaey4kaSIaaGimaiaac6cacaaIWaGaaGimaiaaisdadaqa daqaamaalaaabaGaamitaaqaaiaadseaaaaacaGLOaGaayzkaaWaaW baaSqabeaacaaIYaaaaaGccaGLBbGaayzxaaaaaa@521C@        (4)

which is available in the range of 1.8 ≤ L/D ≤ 9.0. Moreover, it can be seen from Eqn. (4) that the relative critical collapse thickness depends on both the amount and the length-to-diameter ratio of explosive charge, because the amount of explosive charge can be determined by the length-to-diameter ratio of explosive charge under certain density and diameter of explosive charge.


It was observed that relative critical collapse thickness was nonlinear and monotonic increasing function with respect to the length-to-diameter ratio of the explosive charge (or the amount of explosive charge for that the amount and the length-to-diameter ratio of explosive charge are interactional in some sense). And increasing tendency turns to be slower down after the length-to-diameter ratio of explosive charge is larger than 5, which implies that, for relatively larger length-to-diameter ratios, the energy dissipation will exert more and more influence on the relative critical collapse thickness, due to easily release of the explosion gas as well as redistribution of the explosion energy.



The critical collapse condition at the rear-surface of finite thickness targets of concrete under internal explosions was investigated experimentally in this paper. Both the critical collapse thicknesses and the critical amounts of explosive charge under eight different depths of burial were obtained. Combining the obtained experimental data with the dimensional analysis approach, the function relation between the critical collapse thickness and the length-to-diameter ratio of the explosive charge (the critical amount of the explosive charge) was obtained. It was found that relative critical collapse thickness was a nonlinear and monotonic increasing function of the length-to-diameter ratio or the amount of explosive charges and the increasing tendency turns to be slower after the length-to-diameter reaches upto 5. Length-to-diameter ratio of explosive charge has considerable effects on the rear-surface failure of concrete targets in finite thickness.


1. Fu, Yuesheng; Zhang, Qingming & Wang, Haijun. Study on formative mechanism of blasting crater in reinforced concrete under internal blast loading. Key Eng. Mater., 2006, 1645-648.

2. Wang, Z.L.; Konietzky, H.  & Huang, R.Y. Elastic-plastic-hydrodynamic analysis of crater blasting in steel fiber reinforced concrete. Theor. Appl. Fract. Mech., 2009, 52(2), 111-16.

3. Wu, Haijun; Huang, Fenglei; Fu, Yuesheng; Zhang, Qingming & Ma, Ai-e. Numerical simulation of reinforced concrete breakage under internal blast loading. Trans. Beijing Inst. Tech., 2007, 27(3), 200-204. (in Chinese)

4. Xu, K.  & Lu, Y. Numerical simulation study of spallation in reinforced concrete plates subjected to blast loading. Comput. Struct., 2006, 84(6), 431-38.

5. Zhou, X.Q.; Kuznetsov, V.A.; Hao, H.  & Waschl, J. Numerical prediction of concrete slab response to blast loading. Int. J. Impact Eng., 2008, 35(10), 1186-200.

6. Zhou, X.Q.  & Hao, H. Mesoscale modelling and analysis of damage and fragmentation of concrete slab under contact detonation. Int. J. Impact Eng., 2009, 36(12), 1315-326.

7. Lu, Yong & Xu, Kai. Prediction of debris launch velocity of vented concrete structures under internal blast. Int. J. Impact Eng., 2007, 34(11), 1753-767.

8. Tai, Y.S.; Chu, T.L.; Hu, H.T. & Wu, J.Y. Dynamic response of a reinforced concrete slab subjected to air blast load. Theor. Appl. Fract. Mech., 2011, 56(3), 140-47.

9. Wang, Wei; Zhang, Duo; Lu, Fangyun; Wang, Song-Chuan & Tang, Fujing. Experimental study on scaling the explosion resistance of a one-way square reinforced concrete slab under a close-in blast loading. Int. J. Impact Eng., 2012, 49, 158-64.

10. Yi, Na-Hyun; Kim, Jang-Ho Jay; Han, Tong-Seok; Cho, Yun-Gu & Lee, Jang Hwa. Blast-resistant characteristics of ultra-high strength concrete and reactive powder concrete. Constr. Build. Mater., 2012, 28(1), 694-707.

11. Rabczuk, T. & Eibl, J. Modelling dynamic failure of concrete with meshfree methods. Int. J. Impact Eng., 2006, 32(11), 1878-897.

12. Ohtsu, Masayasu; Uddin, Farid A.K.M.; Tong, Weiguang & Murakami, Kiyoshi. Dynamics of spall failure in fiber reinforced concrete due to blasting. Constr. Build. Mater., 2007, 21(3), 511-18.

13. Yuan, Lin; Gong, Shunfeng & Jin, Weiliang. Spallation Mechanism of RC Slabs Under Contact Detonation. Trans. Tianjin Univ., 2008, 14, 464-69.

14.Liu, Yan; Duan, Zhuoping; Huang, Fenglei & Wang, Xinsheng. Damage Effects of Explosion of Shelled Explosive in Concrete. Def. Sci. J., 2010, 60(6), 672-77.


Dr Zhuo-Ping Duan obtained his PhD from Beijing Institute of Technology (BIT), China, in 1994. Presently, he is Professor at State Key Laboratory of Explosion Science and Technology, BIT. His research areas are explosion damage technology and its applications.

Ms Hai-Ying Zhang has received her Bachelor degree in Engineering Mechanics from Taiyuan University of Technology in 2007. Presently, pursuing her PhD in Engineering Mechanics from BIT, China.

Dr Hai-Jun Wu obtained his PhD from Beijing Institute of Technology (BIT), China, in 2003. He is working as an Associate Professor at State Key Laboratory of Explosion Science and Technology, BIT. His research interests are explosion & impact dynamics, dynamic behavior of materials and dynamics numerical simulation.

Dr Yan Liu obtained his PhD from Beijing Institute of Technology (BIT), China, in 2003. He is working as a Professor at State Key Laboratory of Explosion Science and Technology, BIT. His research interests are explosion mechanics, the vulnerability of targets, and damage effects.

Dr. Zhuo-Cheng Ou obtained his PhD (Solid Mechanics) from Xi’an Jiaotong University (XJTU), China, in 2003. Currently working as a Professor in the State Key laboratory of Explosion Science and Technology, BIT, China. His research areas include impact dynamics of solids, dynamic fracture mechanics and fractal fracture mechanics.

Dr Feng-Lei Huang obtained his PhD from Beijing Institute of Technology (BIT), P.R. China in 1992. Presently, he is a professor at State Key Laboratory of Explosion Science and Technology, BIT and his research areas are explosion dynamics, explosion damage and explosion protection.