Mathematical Problems in Engineering
Volume 2016 (2016), Article ID 2628174, 10 pages
http://dx.doi.org/10.1155/2016/2628174

Research Article

Prediction of Collapse Scope of Deep-Buried Tunnels Using Pressure Arch Theory

Yingchao Wang,1,2 Hongwen Jing,1,2 Qiang Zhang,¹ Ning Luo,¹ and Xin Yin¹

¹State Key Laboratory for Geomechanics and Deep Underground Engineering, China University of Mining and Technology, Xuzhou, Jiangsu 221116, China
²School of Mechanics and Civil Engineering, China University of Mining and Technology, Xuzhou, Jiangsu 221116, China

Received 16 April 2016; Accepted 20 June 2016

Academic Editor: Giovanni Garcea

Copyright © 2016 Yingchao Wang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Tunnel collapse remains a serious problem in practice. Effective prediction methods on tunnel collapse are necessary for tunnel engineering. In this study, systematic study on the pressure arch was presented to predict tunnel collapse. Multiple factors under different conditions were considered. First, the pressure arch was described as a certain scope in comparison with the lowest pressure arch line. Then, a deep-buried circular tunnel was selected as the investigated object. Its collapse scope was analyzed using the lowest pressure arch line. Meanwhile, the main influence from the ground stress field was considered. Different modes of ground stress fields were investigated in detail. The results indicate that the collapse scope varies with different ground stress fields. Determination on the collapse scope is strongly affected by the judgment standard of the pressure arch. Furthermore, a selected case was analyzed with the pressure arch. The area and the height of tunnel collapse were calculated with multiple factors, including ground stress field, judgment standard, and lateral pressure coefficient. Finally, selected results were compared with relevant previous researches, and reasonable results were obtained. The present results are helpful for further understanding of the tunnel collapse and could provide suitable guidance for tunnel projects.

1. Introduction

Roof collapse of underground cavities is a practical problem in geotechnical engineering, which usually arises not only in shallow tunnels but also in deep-buried tunnels. Recently, deeper and longer tunnels have been constructed or are under construction in hydropower and transportation fields in China [1]. During the construction of deep-buried tunnels, collapse frequently occurs, which leads to great economic losses and heavy casualties [2–4]. Many efforts have been made to prevent deep-buried tunnels from collapsing, including monitoring, measurement, and prediction. Meanwhile, collapse prediction for deep-buried tunnels is imperative during tunnel construction.

Deep-buried tunnels present many complicated engineering problems [5]. First, tunneling usually traverses multiple geological conditions, like deep underground strata and thick or extremely thick layered rock strata. Second, the presence of cracks and fractures in rock banks can directly induce tunnel collapse [3, 6]. Third, tunnel collapse is strongly affected by random variation in the mechanical properties of the rock in situ. In addition, tunnel collapse can be caused not only by these natural factors [7] but also by human factors, including blast vibration [8], construction conditions, support design, and unexpected failures at work. As a result, more and more attention has been paid to tunnel collapse, particularly for deep-buried tunnels. In order to better understand tunnel collapse, much research has been done on estimating the instability of tunnel roofs and heading faces, using several approaches, including numerical simulation, experimental methods, and analytical approaches [9–11]. However, it is difficult to investigate tunnel collapse experimentally, since recent measurement techniques are limited by the multiscale complexity of the system. Numerical simulations are usually accompanied by multiple assumptions with a complex calculation process. Due to its feasibility and conveniences, the analytical approach has been widely employed to estimate the stability problems in tunnel engineering.

Previous studies on tunnel collapse have focused on analysis techniques, which have been proposed for evaluating the stability of tunnels and investigating the so-called arching effect for both soft soils and hard rocks [12, 13]. Meanwhile, some techniques depended on the specific prepartitioning of the soil domain with theorems of limit analysis [14]. Three independent methods were used to evaluate the roof collapse of underground rectangular cavities for a range of geometries and rock properties [15]. Moreover, an analytical upper bound method, proposed by Fraldi and Guarracino, was used to predict plastic collapse in circular rock tunnels [6, 16]. According to Hoke-Brown failure criterion and the upper bound theorem of limit analysis, tunnel collapse under seepage forces condition and three-dimensional failure mechanism of a rectangular cavity were studied in detail [9, 17–21]. Senent et al. [22] studied the face stability of circular tunnels, which were excavated in heavily fractured and “low quality” rock masses, and the nonlinear Hoke-Brown (HB) failure criterion was applied. Useful results have been obtained with the abovementioned methods. However, due to the complex features of the engineering conditions, it is difficult for one method to obtain useful results when predicting collapse position and the collapse height. In addition, an accurate method of determining the loading mode has yet to be found, which should be acknowledged and could be applied to a wide range of tunnel collapses. Therefore, further research on predicting tunnel collapses is indeed necessary. To estimate the instability of tunnel roofs, a reliable and feasible method is still needed to calculate the collapse scope more accurately for deep-buried tunnels under different ground stress fields.

In practice, the main issue faced by engineering is the analysis of roof failure occurring in deep-buried tunnels [9]. It is still difficult to accurately calculate the key parameters of pressure in the surrounding rock and the relative arching mechanism is also ambiguous [23]. Meanwhile, the evolution characteristics of the pressure arch for a tunnel are difficult to determine exactly. Generally, there is a lack of theoretical solutions, not only for estimating the instability of tunnel roof, but also for predicting the scope of collapse for deep-buried tunnels, especially under different stress fields. However, there are many methods currently being investigated by engineers and researchers, including the analogue model with experience, the structural-geological model, the constrained model with convergence, and the load model. The load model is regarded as one of the most useful models for calculating the pressure of the surrounding rock in tunnel engineering. In detail, the load model involves the restrictive effect of upper strata on the surrounding rock, which is regarded as the main effect of the pressure arch. Research on the pressure arch has been conducted for a long time, but its application is very limited in some conditions.

In the present work, pressure arch theory was employed as an analysis method to predict the collapse scope in deep-buried tunnels. First, the pressure arch derived from pressure arch theory is described in detail. Meanwhile, the lowest pressure arch curve is presented for prediction. Based on the lowest pressure arch, the presented pressure arch curve was applied to predict the collapse scope for a deep-buried circular tunnel. Then, the main influences of the ground stress field related to tunnel collapses were investigated. Moreover, the area and the height of the collapse scope were calculated and analyzed for the deep-buried circular tunnel. Furthermore, multiple factors, including ground stress field, judgment standard, and lateral pressure coefficient, were considered when calculating the area and height of tunnel collapse to find the tendency and variation. Selected results from the present study were compared with results obtained by other researchers. These results are helpful to estimate the overall burden on the lining and can guide further study of prediction methods for tunnel collapses.

2. The Pressure Arch Theory

The present pressure arch was derived from pressure arch theory but differs from that theory in key ways, which are described as follows.

Before the tunnel excavation, all the surrounding rock was under pressure due to the action of gravity. Thus, the surrounding rock along any inner curve was in a compression state along three axes, known as the triaxial compression state, as shown in Figure 1. For example, A and B are arbitrary curves of the pressure arch, which are stable.

Figure 1: Pressure arch of unexcavated tunnel.

After tunnel excavation, due to the unloading effect, the stress in the surrounding rock around the tunnel contour changes, and the surrounding rock within different scopes might have different stress levels. Thus, there is no guarantee that the surrounding rock along any arbitrary curve is in a compression state. There should be some surrounding rock in a tensile state with different tensile stress levels. Meanwhile, there exists a curve with a tensile stress of zero. This curve is used for partitioning the region of the surrounding rock above the tunnel contour, which can be divided into two main regions.

Above the curve with a tensile stress of zero, the surrounding rock within a certain scope is thought to be in the unidirectional compression state, which seems like an arch and is defined as the pressure arch in the present study. Furthermore, the curve with a tensile stress of zero is selected as the lowest pressure arch curve. Any arbitrary curve above the lowest pressure arch curve, like A in Figure 1, is called the pressure arch. On the pressure arch curve, the surrounding rock is still in a compression state, as shown in Figure 2.

Figure 2: Pressure arch of excavated tunnel.

On the other hand, below the lowest pressure arch line, the surrounding rock within a limited scope above the tunnel contour is in a tension state with instability, since the tensile bearing capacity is low for the cracked surrounding rock. The unstable state of the surrounding rock could cause collapse-induced damage and lead to tunnel collapse. The scope of instability within the surrounding rock is regarded as the collapse scope in the present study.

According to the abovementioned analysis, the lowest pressure arch line is significant for predicting collapse of the surrounding rock, which should be gained in the first place. Then, the collapse would be gained accordingly when the area and the height of collapse scope for tunnels are obtained.

3. Results and Discussion

3.1. The Pressure Arch with Bidirectional and Constant Ground Stress

3.1.1. Results of Different Lateral Pressure Coefficients with Different Judgment Standards

The deep-buried circular tunnel was selected as the investigated object, which was analyzed under bidirectional and constant ground stress, as shown in Figure 3. The parameters to be analyzed are introduced as follows: the vertical component of nature stress for rock is ; the horizontal component is ; the lateral pressure coefficient is defined as ; the radial stress is ; the tangential stress is ; and the minimum principal stress is . Additionally, it was assumed that the compressed state had a positive stress value, while the tension state had a negative stress value. Tangential stress () occurs in the tensile state, which was selected as the judgment standard for the pressure arch. Additionally, the lateral pressure coefficient is . Then, the lowest pressure arch line is calculated, as shown in Figure 4.

Figure 3: Deep tunnel with bidirectional and constant ground stress.

Figure 4: Pressure arch (, , bidirectional and constant ground stress).

In order to calculate the scope of tunnel collapse, a plane rectangular coordinate system was established. As shown in Figure 5, there are three key points in the pressure arch line: , , and . The equation to determine the lowest pressure arch curve can be obtained as follows:

Figure 5: Scope of tunnel collapse (, , bidirectional and constant ground stress).

Equation (1) can also be expressed as , while the equation to determine the tunnel contour curve is presented as

Meanwhile, (2) can also be expressed as .

Therefore, the scope of tunnel collapse can be calculated as follows:

Accordingly, the average height of tunnel collapse can be calculated as follows:

When the tangential stress () in the tensile state is selected as the judgment standard of the pressure arch and the lateral pressure coefficient is , the pressure arch curve appears as shown in Figure 6. According to the abovementioned calculation, the maximum height of tunnel collapse is 1 m, the area of tunnel collapse is 2.05 m², and the average height of tunnel collapse is .

Figure 6: Pressure arch (, , bidirectional and constant ground stress).

When the tangential stress () in the tensile state is selected as the judgment standard of the pressure arch and the lateral pressure coefficient is , the pressure arch curve appears as shown in Figure 7. It is possible for the pressure arch curve and the tunnel contour to coincide with each other.

Figure 7: Pressure arch (, , bidirectional and constant ground stress).

When the minimum principal stress () in a state of tension is selected as the judgment standard of the pressure arch and the lateral pressure coefficient is , the pressure arch curve is the same as in Figure 7, and the surrounding rock is in the compression state. When the lateral pressure coefficient is , the pressure arch curve appears as shown in Figure 8.

Figure 8: Pressure arch (, , bidirectional and constant ground stress).

In order to calculate the area of tunnel collapse, a plane rectangular coordinate system was established. As shown in Figure 9, there are three key points in pressure arch line, , , , and the equation related to the pressure arch line can be obtained as follows:

Figure 9: Calculation of scope of pressure arch (, , bidirectional and constant ground stress).

Equation (5) is changed to , and the equation related to the tunnel contour curve is

In another form, (6) can be presented as . Thus, the collapse area can be calculated as

Thus, the maximum height of tunnel collapse is 10.9m, and the average height of tunnel collapse is .

When the tangential stress in the tensile state is selected as the judgment standard of the pressure arch, the lateral pressure coefficient is , and the pressure arch line is shown in Figure 10. According to the abovementioned calculating process, the maximum height of tunnel collapse is equal to 3.6m, the area of tunnel collapse is equal to 13.38 m², and the average height of tunnel collapse is .

Figure 10: Pressure arch (when , , bidirectional and constant ground stress).

3.1.2. Comparison of Results under Different Conditions

According to the abovementioned analysis, when the vertical stress is three times greater than the horizontal stress (), the tunnel vault would undergo tensile stress. The area undergoing tensile stress increases as the lateral pressure coefficient () decreases.

The effect of the lateral pressure coefficient () on collapse scope is shown in Figure 11, and the effect of the lateral pressure coefficient () on collapse height is shown in Figure 12. As shown in Figure 11, the collapse scope decreases as the lateral pressure coefficient () increases. When the lateral pressure coefficient () is 1/3, the collapse scope is 0. When the tangential stress () in the tensile state is selected as the judgment standard of the pressure arch, the relationship between collapse scope and the lateral pressure coefficient () exhibits a linear tendency. Collapse height increases with the lateral pressure coefficient, as shown in Figure 12.

Figure 11: Effects of the lateral pressure coefficient on the collapse scope.

Figure 12: Effects of the lateral pressure coefficient on the collapse height.

It is important to note that the judgment standard of the pressure arch has a great influence on the pressure arch curve. When the minimum principal stress () in tension is selected as the judgment standard, the collapse scope is much larger than that of tangential stress () in tension.

When the horizontal stress is too high, the definition of can be changed from to , and the rock mass in a tensile state will occur at the waist of the arch.

3.2. The Pressure Arch with the Ground Stress Proportional to the Depth

3.2.1. The Theoretical Results of Ground Stress

As shown in Figure 13, the ground stress is proportional to the depth for a deep-buried circular tunnel. The stress field of surrounding rock with gravity has been solved using functions of complex variables [24]. In this study, the stress field of surrounding rock with any linear ratio was calculated using the elastic mechanics method.

Figure 13: Deep tunnel with the ground stress proportional to the depth.

( 1) Original Stress Field of Surrounding Rock before Tunnel Excavation. The stress field in Cartesian coordinates was converted to polar coordinates, and the results are shown as follows:

where , , and are radial stress, tangential stress, and shear stress for Point 0, respectively. is the buried depth of central tunnel, represents the ratio of the vertical stress to the depth, and represents the ratio between the horizontal stress and the depth. In addition, is the polar radius, and is the polar angle.

( 2) Unloading after Tunnel Excavation. After tunnel excavation, the unloading of the working face, , can be divided into four parts as follows:

(i)The first part of the stress field induced by unloading is as follows:

(ii)The second part of the stress field induced by unloading is as follows:

(iii)The third part of the stress field induced by unloading is as follows:

(iv)The fourth part of the stress field induced by unloading is as follows:

( 3) The Total Stress Field of Surrounding Rock after Tunnel Excavation. The original stress field and the stress field induced by unloading can be added, and the total stress field after excavation can be obtained as follows:

3.2.2. Case Study

In order to analyze the pressure arch in detail, a deep-buried circular tunnel was selected as an example. For the tunnel, the buried depth was 600 m, the radius was 6 m, and the rock mass density was 25 kN/m³ (). The coefficients of lateral pressure with respect to and were assigned values of 0.1 and 0.2, respectively. Poisson’s ratio of rock was 0.2.

When the tangential stress in tension state was selected as the judgment standard of the pressure arch, the typical pressure arch was presented, as shown in Figure 14. There is an obvious difference between tectonic stress fields with and without gradient. For the tectonic stress field with gradient, the value of with tensile stress and tensile scope all decreases obviously when compared with that of the tectonic stress field without gradient.

Figure 14: Pressure arch (, , proportional ground stress).

In order to gain the scope of tunnel collapse, a plane rectangular coordinate system was established. As shown in Figure 15, there are three points in pressure arch line, , , and , and the equation related to the pressure arch line can be obtained as follows:

Figure 15: Scope of tunnel collapse (, , proportional ground stress).

With another form, (15) can be presented as .

The equation relative to tunnel contour line can be described as

With another form, (16) can be presented as .

Then, the area of tunnel collapse can be calculated as

Therefore, the maximum height of tunnel collapse is equal to 0.5 m, the area is equal to 1.27 m², and the average height of tunnel collapse is .

When the minimum principal stress in tension state was selected as the judgment standard of the pressure arch, the pressure arch of is shown in Figure 16, while the pressure arch of is presented in Figure 17.

Figure 16: Pressure arch (, , proportional ground stress).

Figure 17: Pressure arch (, , proportional ground stress).

In order to gain the scope of tunnel collapse, a plane rectangular coordinate system was established. As shown in Figure 18, there are five key points in pressure arch line: , , , , and . The maximum height of tunnel collapse is equal to 2.5 m, and the equation relation of the pressure arch line can be obtained as follows:

Figure 18: Scope of tunnel collapse (, , proportional ground stress).

With another form, (18) can be presented as .

The equation relative to tunnel contour line is

Equation (19) can be changed to another form; that is, . Thus, the scope of tunnel collapse can be calculated as

Therefore, the average height of tunnel collapse is .

4. Comparison of Research Results

4.1. Comparison between Different Ground Stresses

For the area and height of collapse, comparisons between different ground stresses were presented. As shown in Figure 19, there exists an obvious difference between bidirectional ground stress and proportional ground stress. Furthermore, ground stress has an obvious effect on the area and height of collapse. The collapse area is the largest under bidirectional ground stress when the minimum principal stress in tension state is selected as the judgment standard of the pressure arch. However, the scope of collapse is the smallest under proportional ground stress when the tangential stress in tension state is selected as the judgment standard of the pressure arch. Moreover, the variation tendencies for the maximum height and average height are similar to that of the collapse area changing with different ground stresses.

Figure 19: The comparison of different ground stresses ().

Therefore, when the minimum principal stress in tension is selected as the judgment standard, the area of the tensile area is obviously larger than that when tangential stress in tension is selected as the judgment standard. Considering engineering safety, the minimum principal stress in tension is the better judgment standard when analyzing the pressure arch.

4.2. Comparison with Previous Work

In order to verify the proposed method, the results were compared with previous studies. Selected conditions were presented for comparison, and the area of tunnel collapse was calculated based on the pressure arch theory. The corresponding parameters are shown as follows: runnel radius is , rock mass density is kN/m³, the buried depth is , the coefficient of lateral pressure is 0.1, and rock Poisson’s ratio is equal to 0.2. The result is presented in Figure 20.

Figure 20: Comparison between the present result and the relative results of collapse area for deep circular tunnel.

Fraldi and Guarracino [6] have also investigated a similar circular tunnel and obtained the area of collapse scope with varying tunnel radius (Figure 20), which could be used as a comparison to validate the present model. As shown in the figure, the collapse scope in the present study approximates to that of Fraldi and Guarracino. It indicates that the results derived in the present study agree with that of previous studies to predict the collapse scope in deep-buried circular tunnels. The present method could be used to effectively predict tunnel collapse.

5. Conclusion

With the consideration of the serious problem with tunnel collapses, the method with the pressure arch was put forward to predict tunnel collapse in the present work. The method was presented with the lowest pressure arch line, which was used to analyze the collapse scope for a deep-buried circular tunnel. Then, the main influence from ground stress field was considered with different conditions: the ground stress proportional to the depth and the bidirectional and constant ground stress. Moreover, the area and the height of tunnel collapse were calculated in detail with the consideration of multiple factors, like ground stress field, judgment standard, and lateral pressure coefficient. It was found that the scope of tunnel collapse differs with different modes of ground stress field, and the collapse scope is closely related to the judgment standard for the pressure arch. For bidirectional ground stress, the collapse scope decreases with the increasing of the lateral pressure coefficient. For the area and the height of collapse, there exist obvious differences between bidirectional ground stress and proportional ground stress, and the ground stress has an obvious effect on the area and the height of collapse. Furthermore, in order to verify the proposed method, corresponding comparison was presented with the previous findings, and suitable agreements could be gained. The results are significant for further understanding of the essence of tunnel collapse, which could provide more suitable guidance for future tunnel projects.

Nomenclature

: | The average height of tunnel collapse (m)
: | The buried depth of central tunnel (m)
: | The ratio of the vertical stress to the depth
: | The ratio of the horizontal stress to the depth
: | The area of tunnel collapse (m²)
: | The polar radius (m)
: | The rock mass density (kg/m³)
: | The polar angle (°)
: | Rock Poisson’s ratio
: | The lateral pressure coefficient
: | The minimum principal stress (Pa)
: | The horizontal component of nature stress (Pa)
: | The vertical component of nature stress (Pa)
: | The radial stress (Pa)
: | The tangential stress (Pa).

Competing Interests

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

Acknowledgments

Financial supports from the Fundamental Research Funds for the Central Universities (no. 2014QNB37) and the National Natural Science Foundation of China (no. 41572263 and no. 11502282) are sincerely acknowledged.

References

T. H. Ma, C. A. Tang, L. X. Tang, W. D. Zhang, and L. Wang, “Rockburst characteristics and microseismic monitoring of deep-buried tunnels for Jinping II Hydropower Station,” Tunnelling and Underground Space Technology, vol. 49, pp. 345–368, 2015.

A. Aydin, A. Ozbek, and I. Cobanoglu, “Tunnelling in difficult ground: a case study from Dranaz tunnel, Sinop, Turkey,” Engineering Geology, vol. 74, no. 3-4, pp. 293–301, 2004.

H.-S. Shin, Y.-C. Kwon, Y.-S. Jung, G.-J. Bae, and Y.-G. Kim, “Methodology for quantitative hazard assessment for tunnel collapses based on case histories in Korea,” International Journal of Rock Mechanics and Mining Sciences, vol. 46, no. 6, pp. 1072–1087, 2009.

J.-H. Shin, I.-K. Lee, Y.-H. Lee, and H.-S. Shin, “Lessons from serial tunnel collapses during construction of the Seoul subway line 5,” Tunnelling and Underground Space Technology, vol. 21, no. 3-4, pp. 296–297, 2006.

D.-F. Chen, X.-T. Feng, D.-P. Xu, Q. Jiang, C.-X. Yang, and P.-P. Yao, “Use of an improved ANN model to predict collapse depth of thin and extremely thin layered rock strata during tunnelling,” Tunnelling and Underground Space Technology, vol. 51, pp. 372–386, 2016.

M. Fraldi and F. Guarracino, “Evaluation of impending collapse in circular tunnels by analytical and numerical approaches,” Tunnelling and Underground Space Technology, vol. 26, no. 4, pp. 507–516, 2011.

L. P. Li, Q. H. Wang, S. C. Li et al., “Cause analysis of soft and hard rock tunnel collapse and information management,” Polish Journal of Environmental Studies, vol. 23, pp. 1227–1233, 2014.

D. Uystepruyst and F. Monnoyer, “A numerical study of the evolution of the blast wave shape in rectangular tunnels,” Journal of Loss Prevention in the Process Industries, vol. 34, pp. 225–231, 2014.

C. B. Qin, X. L. Yang, Q. J. Pan, Z. B. Sun, L. L. Wang, and T. Miao, “Upper bound analysis of progressive failure mechanism of tunnel roofs in partly weathered stratified Hoek-Brown rock masses,” International Journal of Rock Mechanics and Mining Sciences, vol. 74, pp. 157–162, 2015.

R. Bubbico, B. Mazzarotta, and N. Verdone, “CFD analysis of the dispersion of toxic materials in road tunnels,” Journal of Loss Prevention in the Process Industries, vol. 28, pp. 47–59, 2014.

C. Vianello, B. Fabiano, E. Palazzi, and G. Maschio, “Experimental study on thermal and toxic hazards connected to fire scenarios in road tunnels,” Journal of Loss Prevention in the Process Industries, vol. 25, no. 4, pp. 718–729, 2012.

C. J. Lee, B. R. Wu, H. T. Chen, and K. H. Chiang, “Tunnel stability and arching effects during tunneling in soft clayey soil,” Tunnelling and Underground Space Technology, vol. 21, no. 2, pp. 119–132, 2006.

A. S. Osman, R. J. Mair, and M. D. Bolton, “On the kinematics of 2D tunnel collapse in undrained clay,” Geotechnique, vol. 56, no. 9, pp. 585–595, 2006.

A. Drescher and E. Detournay, “Limit load in translational failure mechanisms for associative and non-associative materials,” Geotechnique, vol. 43, no. 3, pp. 443–456, 1993.

A. M. Suchowerska, R. S. Merifield, J. P. Carter, and J. Clausen, “Prediction of underground cavity roof collapse using the Hoek-Brown failure criterion,” Computers and Geotechnics, vol. 44, pp. 93–103, 2012.

M. Fraldi and F. Guarracino, “Limit analysis of collapse mechanisms in cavities and tunnels according to the Hoek-Brown failure criterion,” International Journal of Rock Mechanics and Mining Sciences, vol. 46, no. 4, pp. 665–673, 2009.

C.-B. Qin, Z.-B. Sun, and Q. Liang, “Limit analysis of roof collapse in tunnels under seepage forces condition with three-dimensional failure mechanism,” Journal of Central South University, vol. 20, no. 8, pp. 2314–2322, 2013.

F. Huang, C.-B. Qin, and S.-C. Li, “Determination of minimum cover depth for shallow tunnel subjected to water pressure,” Journal of Central South University, vol. 20, no. 8, pp. 2307–2313, 2013.

X. L. Yang and F. Huang, “Three-dimensional failure mechanism of a rectangular cavity in a Hoek-Brown rock medium,” International Journal of Rock Mechanics and Mining Sciences, vol. 61, pp. 189–195, 2013.

X. L. Yang, J. S. Xu, Y. X. Li, and R. M. Yan, “Collapse mechanism of tunnel roof considering joined influences of nonlinearity and non-associated flow rule,” Geomechanics and Engineering, vol. 10, no. 1, pp. 21–35, 2016.

X.-L. Yang and J.-H. Yin, “Slope stability analysis with nonlinear failure criterion,” Journal of Engineering Mechanics, vol. 130, no. 3, pp. 267–273, 2004.

S. Senent, G. Mollon, and R. Jimenez, “Tunnel face stability in heavily fractured rock masses that follow the Hoek-Brown failure criterion,” International Journal of Rock Mechanics and Mining Sciences, vol. 60, pp. 440–451, 2013.

J. H. Yang, S. R. Wang, Y. G. Wang, and C. L. Li, “Analysis of arching mechanism and evolution characteristics of tunnel pressure arch,” Jordan Journal of Civil Engineering, vol. 9, no. 1, pp. 125–132, 2015.

G. F. Wang, Tunnel Calculation (The Theory of Joint-Deformation), Chengdu Science & Technology University Press, Chengdu, China, 1992.