The theoretical formula of the maximum internal forces for circular tunnel lining structure under impact loads of the underground is deduced in this paper. The internal force calculation formula under different equivalent forms of impact pseudostatic loads is obtained. Furthermore, by comparing the theoretical solution with the measured data of the top blasting model test of circular formula under different equivalent forms of impact pseudostatic loads are obtained. Furthermore, by comparing the theoretical solution with the measured data of the top blasting model test of circular tunnel, it is found that the proposed theoretical results accord with the experimental values well. The corresponding equivalent impact pseudostatic triangular load is the most realistic pattern of all test equivalent forms. The equivalent impact pseudostatic load model and maximum solution of the internal force for tunnel lining structure are partially verified.

1. Introduction

In recent years, the nuclear leakage events have sparked a new national self-examination. The safety and stability problems of the underground structure caused by the impact or blasting load must be paid special attention to. The blasting demolition at the top of underground, chemical explosion devised by terrorism, gas pressure extrusion, and release of shallow layer containing gas and so forth will result in instant huge impact load on the structure of tunnel lining. The structure crack or damage of the tunnel under the impact load will lead directly to the serious nuclear leakage events, with Japan’s Fukushima nuclear power plant explosion causing the serious nuclear leakage from 2011. The former Soviet Union Chernobyl nuclear power plants exploded at 1:23 in the morning of April 26, 1986, with more than 8 tons of highly radioactive materials mixed with graphite fragments and nuclear fuel burning debris spewing out. The radiation pollution caused by the nuclear leakage accident is equivalent to 100 times of the radioactive pollution caused by the atomic bomb explosion in Hiroshima, Japan. Even 20 days after the accident, the temperature of the center of the nuclear reactor is still as high as 270 degrees Celsius. It caused 10 times the number of cancer deaths caused by the accident in the United Nations official estimates, with a global total of 2 billion people affected by the Chernobyl accident; 270 thousand people suffer from cancer, which killed more than 93 thousand people. Experts estimate that eliminating the effects of this catastrophe will take at least 800 years.

The underground structure of tunnel lining locates in the enclosure space; the security problem and the precise evaluation of the dynamic response to the blasting or impact loading becomes an important engineering issue and there is a need to give mathematical model and solution of internal force [16].

First, determining the strength and internal force of structure response under the blast loading is the key factor to solve the dynamic response problem for tunnel lining. By now, the researches of stability of the structure for tunnel lining under impact load are more concentrated in blasting load or gas outburst. Generalized impact action blocks and changes the motion of a moving object. The difference between impact load and impact pseudostatic load is the energy form related to the instantaneous time effect or not. Generally, the quantity of explosive and detonation pressure of explosive shock can be determined according to the energy release and the size of the gas pressure to determine the impact of the air pressure. Blasting seismic wave attenuation can be equivalent to the dynamic load effect of gas outburst on tunnel lining structure. Currently, the theoretical study of the tunnel lining and the internal force response under impact load is relatively little. This paper introduces the free deformation method as the theoretical basis. Derivation of impact load equivalent pseudostatic model is given. Load pattern on the circular tunnel lining under the instantaneous maximum internal force calculation formula is compared and analyzed. By comparing the theoretical analysis to the test data and the numerical simulation results, it intends to get the internal force response of the failure mechanism and the optimal load pattern and the theoretical results can be preliminarily verified by the back analysis [1].

2. The Internal Force Formula of the Circular Tunnel Lining under the Blasting Loading

Considering significantly attenuated effect of cover layer to the impact load and wave, it is supposed that impact loaded on structure of tunnel lining is equal to the earth surface. The additional stresses during elastic wave propagation are ignored. To simulate the stress and deformation of tunneling lining under impact or blasting dynamic load, engineers commonly adopt the simple and safe simplified method, which means the impact pseudostatic load multiplied by a dynamic load factor to estimate the maximum impact load. The core idea of the proposed method is first to find out the maximum equivalent impact pseudostatic stress mode under the impact or blasting load.

The theory formula of the equivalent model for impact load on circular tunnel lining is derived based on using the free deformation method, which is the classical theory method for underground structure [7].

2.1. The Equivalent Impact Load Model

The equivalent impact load model should be given and the equivalent form caused by blasting seismic wave loads on tunnel lining structure should be determined. Left-side blasting load coefficient of dynamic stress distribution for surrounding rock is given by Figure 1. Figure 1(a) gives the mechanical model of left-side blasting load on the circular underground structure; the radial movement of surrounding rock stress concentration factor is given as in Figure 1(b). Upper blasting load coefficient of dynamic stress distribution of surrounding rock is shown in Figure 2. Figure 2(a) gives the mechanical model of the upper-side blasting load on the circular underground structure; the force model of the tunnel lining structure and the surrounding rock of the dynamic stress concentration factor distribution map is given as in Figure 2(b). Figure 1 can give the maximum acceleration stress distribution map of the left-side blasting load in surrounding rock for lining. Figure 2 can give the maximum acceleration stress distribution map of the upper-side blasting load surrounding rock for lining. From Figures 1 and 2, it is shown that the blasting response of the circular lining is according to the following rules: the maximum point of impact load generally is the maximum stress response point, mostly on angle of 0 degrees near the load; the smaller load is distributed on both sides of the lining; the maximum stress zone is concentrated in the range of degrees around both sides of the impact load; on the far side from the impact load, the stress on both ends of the lining is obviously higher than that in the middle of the lining.

The dynamic stress concentration factor and the map of the maximum acceleration distribution obtained by the two modes of lateral explosion and upper explosion of Figures 1 and 2 present typical symmetrical patterns of butterfly distribution.

The maximum internal force response formula of arbitrary angle on lining is derived under any angle impact load.

First consider the special case of the impacted load on the upper side of the lining; the impact load can be expressed by arbitrary rotating angle; the additional force of lining under impact load is divided into three parts.

The first part is the equivalent impact load, with the blasting load applied on the lining; the second part is the equivalent impact reaction force, which is the reverse force of lining and formation, equal to the equivalent impact load for balance effect; the third part is the proposed rock resistance force on both sides of lining; it is caused by the relative displacement mode between the lining and ground.

We give the case study of the triangle resistance distribution; the load model is shown as in Figure 3. The effective counter force which is under the influence of impact in order to maintain the balance of force provided by the formation of lining is equal to the action of the equivalent load and impact load. The distribution of load model can be traced back to the Japanese triangle resistance method [7].

According to the stress characteristics of structure response under impact load, the force model diagram of the circular tunnel lining under the upper blasting impact is given. Two main factors should be considered for affecting the load and the response of the lining under the impact load: load form and load value. Figure 3(a) shows the force model of the circular lining under the condition of impact load. Figure 3(b) is the mechanical model of additional load for blasting shock.

The equivalent action form of impact load is discussed first. There are two loading modes of simulating the blasting and impact loading of rock: one is calculating the explosion hole pressure by explosive burst detonation theory and then loading the calculating blast action on hole wall directly; second is using the empirical formula to calculate the dynamical peak value of the load and then imposing it on the boundary according to triangular pulse wave form. The former method needs to introduce the state equations of the explosive detonation and the rock mass, which is used for single hole blasting or centralized charge blasting. For the impact of porous blasting, the blasting source distribution area is larger; the simulation is not directly loaded with blasting hole pressure according to the equivalent amount of explosives to concentrate on loading and explosive loading. The simulation of blasting vibration effect is based on triangle pulse wave loading [7]. Suppose that the effect of gas outburst, seismic wave, and blasting attenuation on lining is equivalent. The equivalent additional load caused by impact load is assumed to consist of three parts: the upper part load and the two-side reaction force and the bottom reaction force. The influence of the different distribution forms on the internal force of the lining is discussed in three parts. The assumed distribution of the blasting load is shown in Figure 4.

Assume that the impact load produces the same total equivalent value ; the concentration of Figure 4 in four kinds of load on the upper part of the form is (a) < (b) < (c) < (d); Figure 4(d) is the limit form of concentrated load; the impact load can act as the concentrated load. Stress concentration in the middle of the lining is the most significant. The concentration degree of subgrade load under the same condition of the equivalent impact load is ① < ②.

2.2. The Formula Deduction of Lining Internal Force under Impact Load

The internal force formula of lining response under various forms of equivalent impact load is derived. The force in the process is the total load value of the impact load. Because the computation processes of internal force calculation for different impact loads combination are similar, the paper gives a detailed derivation of the upper curve load for reference, and the rest of the cases only give the result.

Suppose the upper impact equivalent load is sinusoidal and the maximum value is ; the sine load value is in the angle of , and take the range of as the load calculation of the equivalent load. When , two load diagrams are shown in Figure 5; the bending moment at caused by is , of which is the load value at , is the action scope, and is the action distance.

While ,The load diagram is shown as in Figure 5, while : can be calculated after the calculation formula of is gotten asThenPut the values of and into .

While ,While ,The relationship between the equivalent blasting force and isThe axial force calculation is (8), while :While ,Bring the values of and into .

While ,While ,The calculation formula of shear force is as follows: while ,While ,Bring and into .

While ,While ,

The maximum value of instantaneous explosive load and maximum internal force can be obtained by bringing maximum into the formula.

The expression of maximum explosion pressure is given in 1956 by Brown in the equivalent value of blasting load [1] as

1971 Sassa gives the following expressions:

The detonation wave spread in the rock mass gives the maximum pressure in the contact interface of explosive and rock. The maximum pressure generated on the contact interface between the detonation wave and the rock mass is related to the rock physical characteristic. The maximum pressure and the maximum explosion pressure relation can be approximately expressed as follows: is rock density; is wave velocity propagated in rocks; is explosive detonation velocity; 0 is explosive density.

The maximum equivalent load is

The relationship between the peak value of blasting load and the distance can be expressed as is peak value of blasting load; is the distance from the calculate point to blasting hole; is the diameter of the contact surface ().

The Statfield dynamic pressure with the time history can be used aswhere is the load constant, having an amount equivalent to the dynamic pressure generated by the charge of each 1 kg explosive [1].

The function of equivalent total load value on the lining with the time from formulas (16)–(21) iswhere is the maximum value of blasting load on the contact lining surface.

The internal force calculated by the different combination of load form is shown in the Appendix.

3. Test Verification of the Model

3.1. Comparison with the Blasting Model Test

In order to further verify the applicability of the theoretical calculation formula of the instantaneous maximum internal force of the circular lining (Figure 4) under the impact load, the calculation and comparison of the microvibration centrifuge model test results on the top of lining are carried out. The top vibration model test diagram is shown in Figure 6.

The parameters of impact load model test on the upper part of the lining are shown in Table 1 [9]. The comparison and analysis of the test results are carried out as follows.

The bending moment of the structure under impact load can be obtained by converting the microstrain model test results. The model tests show that the mix ratio of sand to gypsum material is gypsum plaster material : sand : water = 1 : 0.8 : 0.5; elastic modulus , which is the reduction elasticity modulus of the prototype concrete . Compare two kinds of thickness with 7.5 mm and 12.5 mm; the elastic modulus of aluminum alloy model test is ; the thickness is 3.8 mm.

The similarity criterion of bending deformation of the model centrifuge test is as follows: is the number of centrifugal accelerations; is Poisson’s ratio of the model; is Poisson’s ratio of the prototype; is Young’s modulus of the model; is Young’s modulus of the prototype; is the thickness of the model lining; is the thickness of the prototype lining.

The next formula is available by the formula of (23): and are suitable for the sand gypsum material; the moment can be calculated by directly substituting the experimental microstrain into the prototype; for the aluminum alloy material, it can be calculated byAccording to the formula of , the instantaneous maximum bending moment of the prototype can be obtained. The elasticity modulus of the prototype concrete is . The results of converted moment value of the measuring point are shown in Table 2.

The equivalent blasting load values are calculated by the third groups of experiments, and the parameters of the test are shown in Table 3 [9]. The theoretical solution of maximum internal force to test model is calculated.

ThenThe attenuated equivalent load is by giving according to the test conditions; the conversion coefficient of black powder and TNT is 0.4; the centrifugal acceleration is 50 g; then 1.25 g black powder is equivalent to TNT of 62.5 kg [9]; thenThat means the equivalent blasting load on the lining of third groups experimental results is 778 kN, considering that the size effect of model test will lead to error; the selection range of calculation load is from 600 kN to 800 kN.

The internal force of the maximum impact load for the circular lining under the condition of the load combination is calculated as shown in Figure 7. The equivalent form of the upper impact load is compared with the uniform load, the curve load, the triangle load, and the concentrated load. The concentrated load represents the extreme case where the impact load directly acts on the lining without attenuation through the surrounding rock.

Select the different calculation combination of blasting load form in the top and bottom; the results are compared as shown in Tables 49; take . The comparison of the scatter plots is shown in Figure 7.

The following can be concluded through comparison of a①, b①, c①, and d①:(1)The influence of the upper load form on the bending moment of the measuring point 1 (top) is obvious, the difference range can reach 50%–100%, and the comparison of the load distribution is very necessary.(2)The characteristics of each group showed a clear trend of convergence; the measured moment value of point 1 represents the relationship a① < b① < c① < d①. This shows that with more concentration in the upper part of the load form the moment value of measuring point 1 will be greater.(3)Compared with each test point data under all cases of combination, it can be known that the measuring points of the bending moment values have changed considerably when the upper load form changed from the uniform load to the curve load. The rule is that the measuring point 1 and point 3 (the top and bottom points of the lining) bending moment value increases and the measuring point 2 and point 4 (the point on left-right sides of the lining) bending moment value decreases.(4)Compared with each test point data of 2, 3, and 4 under all cases of combination of b①, c①, and d①, it can be known that the increase concentration of the upper load and the linear increase effect of the internal force response are decreased. The effect of concentration on both sides and the bottom of the lining is very weak (within 1%).

Compared with the test combination of b① and b②, the result showed the following:(1)The form influence of the ground reaction force on the moment of the measuring point 1 (bottom) is obvious, and the difference range is 90%–50%.(2)The difference range of the combination of b① and b② is larger (range: 10%–50%) for the measuring point 1, point 2, and point 4. This shows that the concentration increase of the subgrade force in the bottom of lining will only affect the bending moment of the action position and have less effect on the other points’ bending moment outside of this range.

Comparing the data calculation in the four groups, it is found that the combination of upper triangular load and subgrade counterforce of the outer triangle reaction model has good coincidence degree with the experimental values. This case also has the certain safety reserve; comparatively, this load combination is the most reasonable one. The internal force calculation results (bending moment, shear force, and axial force) of the upper triangular load and the subgrade counterforce of outer triangle reaction model are shown in Figure 8.

Validating the lining safety theoretical results of the third test groups preliminarily, the calculation results of triangular equivalent loads are selected to verify the results, while the maximum value of the calculated moment is 367.2 kN m, the position is S1 measuring point, and the axial force is 247.8 kN.

According to the maximum compressive classical stress formula of the eccentric compression member,Bring the test parameters and calculation results; the following can be concluded:The compressive strain is considered as the critical value. The lining structure can be considered in a safe condition as less than the concrete elastic ultimate compressive strain of 0.002.

The compression plane strain of the test results isIt can be considered that the lining structure is in a safe state under the case study of blasting load.

3.2. Comparison to Model Test of Pneumatic Impact Load

Compare the pneumatic impact load mode to the model test. The test model box is shown as in Figure 9(a) [10]. The air tube is located in the three (left, right, and down) directions of the model box. The air pressure is equivalent to the triangular load. The additional load is caused by air pressure and the equivalent loading is shown in Figure 9(b). The upper uniformly distributed load is the subgrade counterforce provided by the equilibrium pressure from the lower pressure layer.

The comparison of the additional bending moment and the theoretical calculation results of the lining structure under releasing and applying air pressure to the test is shown in Table 10.

It is found that the calculation results of the additional bending moment for the circular lining under the impact pressure were coincident with the model test data. Four points’ precision accuracy is within 15%, with the highest accuracy reaching 95%. The results for the theory of structure response (Figure 3) moment under impact load calculation results are in good agreement with the experimental values. It can provide support for theoretical calculation of tunnel lining structure response under impact loading.

4. Numerical Simulation to Microvibration Test

MIDAS GTS was used to simulate the microvibration model test. After the model feature was calculated, the response simulation of lining (Figure 3) under blasting load was carried out by using the time history analysis.

First, the material parameters are given in Tables 11 and 12 [10].

Using MIDAS GTS to simulate the blasting load, manually input the blasting load function. The dynamic pressure time history functions are shown in Figure 10. The blasting load is applied on the interface in form of the surface pressure.

Table 13 [10] shows the maximum microstrain simulation value of the lining structure under blasting loads in consistency with the test value very well; most of the gap is less than 30%. Most difference of the measuring points is less than 30%. The numerical simulation results not only verify the experimental results but also provide the basis for the correctness of the theoretical calculation results. Compared with the results of numerical simulation, the theoretical formula of the lining maximum internal force under impact loading has good applicability in the case study.

Using the GTS, simulate and validate test results; the impact pseudostatic load simulations are carried out by constrained lining structure of tunnel lining by the spring boundary. Theoretically, the simulation and experimental results are compared with the results shown in Table 14. The three-coincidence degree is good through comparison; most difference of the points is less than 20%. The applicability of the theoretical calculation results is verified further by the results of impact pseudostatic load simulation.

5. Conclusions

Impact pseudostatic load equivalent model and the maximum internal force solution for circular tunnel lining structure were deduced from classical theory of the free deformation method. The equivalent model and the maximum internal force solution of the underground lining are given. The equivalent load forms of different impact loads are compared preliminarily. The following conclusions and recommendations of the maximum internal force of lining structure under impact load may also be drawn by comparing the theoretical solution, experimental data, and numerical simulation results.

(1) The equivalent form of impact load has a remarkable influence on the calculation results. It was reflected in the internal forces distribution pattern of the lining structure in the load range of impact pseudostatic action.

Two equivalent forms of impact loading on the lining structure of tunnel lining were compared by the curve distributed shape and triangle distributed shape. The simplified triangular distributed load is more convenient and feasible to use with the premise of closing the test results.

(2) The additional load caused by the impact load mainly has two types, the ground counterforce and lateral deformation resistance, which can be calculated by the classical Japanese triangle resistance method. The equivalent mode of formation reaction has two modes, which are external triangle distribution and uniform distribution. For the case of side impact load, the distribution of the outer triangle is obviously better than that of the uniform distribution load distribution, whether from the dynamic stress concentration degree of the surrounding ground or from the matching degree to test data.

(3) The influence of soil parameters on the structural lining response under blasting shock load is the following: the soil density influences the transmission and attenuation of shock wave in surrounding rock and the formation parameters for the effect of structural response under blasting load: soil density influences the detonation wave in the soil layer of propagation and attenuation. The characteristics show that the lower the density, the softer the soil and the more significant attenuation of impact loads. The distance between blasting spot and lining structure is similar to the factor of the soil layer density; the farther the distance, the faster the impact load decay rate. The magnitude of the impact load directly affects the response of the lining; the strain of the structure lining presents a nonlinear increasing relationship under the nonlinear interaction to the surrounding ground.

(4) Comparing the theoretical formula and dynamic load simulation to the microvibration test data, the theoretical results of the impact pseudostatic load equivalent model are verified again by the numerical simulation of the equivalent static load.

The results show that the consistency of the test results, theoretical results, and the simulation results of impact pseudostatic load is very high. Most of the difference point is within 20%, which shows that the theoretical calculation results of the proposed model have good applicability.

The maximum microstrain simulation by the dynamic load on lining is in agreement with the test values very well, with most of the difference points being within the range of 30%. The comparison between the experimental results and the theoretical calculation results of the pressure shock load model is carried out; the difference of the closest point is 2.3%.

In conclusion, the impact pseudostatic load equivalent model is proposed for dynamic response analysis of tunnel lining structure. Good agreement among the theoretical analysis, test data, and the simulation is achieved, which preliminarily validates the present model and method.


See Tables 15, 16, and 17.

Competing Interests

The authors declare that there are no competing interests regarding the publication of this article.


This work is supported by the National Natural Science Foundation of China (no. 51378051 and no. 51678038), the Fok Ying Tung Education Foundation (no. 122009), and the Fundamental Research Funds for the Central Universities (no. 2011JBZ008).