1. Introduction
Owing to rapid economic development, logistics, and the vehicle industry, the construction of tunnels with large-span and extra-large-span sections has boomed in recent years and serves as the optimum scheme to improve road capacity and meet the existing travel demand. Compared to a single-hole project with a small span, an extra-large-span-section tunnel consists of a much larger span, a lower high-span ratio, and a greater construction risk. There are no corresponding design codes and construction standards for these projects. Although some construction experience has been obtained, theoretical research is relatively lacking.
To guarantee the rationality of the supporting structure design for super-span tunnels, the first problem to be solved is to clarify the method of predicting the loose pressure in the surrounding rock mass. At present, the empirical formulas commonly used for the loose pressure are divided into three categories: the equivalent thickness of the crown, statistics of collapse law, and surrounding rock classification, as shown in
Figure 1. The equivalent thickness of the crown takes the gravity of the rock and soil mass above the vault as the primary influence for the loose pressure in the surrounding rock. Typical methods include Protodyakonov’s theory [
1], the Bierbaumer theory [
2], the Terzaghi theory [
3], and Xie Jiaxiao’s theory [
4]. The statistics of collapse law method is based on practical engineering and is the leading engineering practice used in design codes in China [
5,
6,
7,
8]. The correlations are relatively simple and consider the influence of the surrounding rock classification. The surrounding rock classification method derives the surrounding rock loose pressure from the RQD, Q, and RMR systems [
9,
10,
11,
12]. The correlations and applicable conditions of various calculation methods are presented in
Table 1.
There are two deficiencies of the above-mentioned empirical formulas: (1) the calculations are derived or statistically obtained under the conditions of a small-span tunnel, and (2) the influencing factors are significantly different. The above two deficiencies restrict the applications of empirical formulas for loose pressure in the extra-large-span cross-section tunnels.
In addition to the empirical formulas, the upper bound theorem of limit analysis with a strictly theoretical basis has also been used recently to calculate the loose pressure in the surrounding rock mass. It focuses on the ultimate failure state of the rock–soil mass without considering the elastic–plastic deformation process. In addition, this method is not restricted by the span. Some scholars have published papers on the surrounding rock pressure of soil tunnels. Atkinson et al. [
13] discussed the stability of a shallow circular tunnel in a cohesionless soil through the upper bound theorem of limit analysis together with results of tests performed in sand soil. Davis et al. [
14] deduced the upper boundary stability solutions for a shallow tunnel excavation at three different positions under undrained conditions and proposed a new stability evaluation method. In order to analyze the stability of an excavation face, Mollon et al. [
15,
16] established the spatial discretization three-dimensional calculation method based on the limit analysis method, and the results of the calculation approached the real situation. Based on the Hoek–Brown failure criterion, Huang et al. [
17] presented an upper bound solution for the collapse shape for a circular tunnel considered the influence of different geological conditions and pore water pressure. Lei et al. [
18] established explicit expressions to calculate the loose pressure on a shallow, buried tunnel based on the upper bound theorem of limit analysis and the Terzaghi failure mode. Based on the spatial discretization model, Pam et al. [
19,
20,
21] studied the impact mechanism of a heterogeneous stratum, pore water pressure, and a non-circular excavation face on the tunnel face stability. Zhang et al. [
22] proposed an analytical method for calculating the ultimate support force of an excavation face in the clay layer and analyzed the local and whole stability of the excavation face for a large-diameter shield tunnel in practical engineering. According to the tunnel face failure mode in purely cohesive soil proposed by Mollon [
23], Huang et al. [
24] established an analytical model considering the shear failure boundary which further improved the calculation accuracy of the upper bound solution for tunnel face stability in non-homogeneous and anisotropic undrained clay. To assess the stability of a tunnel face in inclined strata, Tu et al. [
25] optimized the 3D rotating rigid body collapse mechanism by improving the spatial discretization technique. However, the failure mechanism constructed in the current research is highly subjective and lacks numerical simulation and field test verification. Such flaws can be solved by combining a numerical analysis and an upper bound limit analysis. With its lower cost, shorter period, and better repeatability, a numerical simulation is helpful for predicting the failure position [
26,
27]. Wang et al. [
28,
29] verified the rationality of the analytical solution for the collapse curve via a vertical displacement cloud picture obtained from a simulation.
In view of the shortcomings of the existing research, this paper discussed the impact of rock conditions and cavern size on the range of the loose pressure zone by combining a numerical analysis and an orthogonal test. The investigation is based on the Bifurcation Tunnel. The limit analysis model of the loose pressure in the surrounding rock is established using the fitted boundary function of the loose zone and the upper bound limit analysis. Based on this model, this paper deduces the upper bound solution for a deep tunnel’s loose pressure, which is verified by pressure-monitoring results obtained for the different sections. The results of this research can provide new ideas for the selection of a supporting structure for similar super-section tunnel projects.
2. Project Overview
The Shenzhen Transit Expressway is a major construction project in Shenzhen, Guangdong Province, China. It begins at Luosha Road, connects with the Hong Kong highway network through the port, and ends at the Huiyan expressway. The total length of the route is 32.5 km, while the total investment is about CNY 5.792 billion. It will become a critical traffic artery serving the development of the eastern coastal areas of South China after its completion.
The eastern transit passage of Shenzhen City starts from Luohu Port, connects with the Hong Kong highway network through the port, and ends at the Huiyan expressway. The project is 32.5 km long in total, and it will become a critical traffic artery serving the development of the eastern coastal areas of South China after its completion. The passage planning is shown in
Figure 2.
The Bifurcation Tunnel is a key control project in the entire transit expressway line and has a total length of 1435 m [
30]. It is composed of three parts: the mainline section, the municipal section, and the port connecting section, of which the mainline section is divided into left and right lines. The Bifurcation Section is located at the entrance of the left line. The large-span section includes the maximum section, Gradient Section 1, Gradient Section 2, and Gradient Section 3. The Bifurcation Section is shown in
Figure 3.
The burial depth of the large-span bifurcation section is 78~106 m. The stratum at the tunnel site consists of artificial soil, sandstone, flood alluvial soil, and rhyolite. The rock mass is relatively complete, with only slightly developed fractures to which the degree of development increases with depth. The geological section is shown in
Figure 4, and photos of the drill core are shown in
Figure 5. The geological investigation, combined with the rock mass type, weathering degree, and joint development, provides the GSI (geological strength index) value of the Bifurcation Section, as shown in
Table 2. Furthermore, the maximum excavation span is 30.01 m, and the maximum excavation area is 428.5 m
2; globally, this is a rare extra-large section highway tunnel. Therefore, the key to the construction project’s structural design and safety is determining the influencing factors and distribution pattern of the loose pressure of the surrounding rock for the large-span section.
3. Study of Factors Influencing the Loose Zone in the Surrounding Rock Mass
The formation of the loose zone due to excavation is an inherent characteristic of the surrounding rock mass. The scope of the loose zone directly determines the loose pressure in the surrounding rock mass. Based on the orthogonal test and the finite difference software FLAC3D, this paper designs multiple group model tests to determine the influence of various factors, such as the surrounding rock quality, cavern size, and tunnel burial depth, on the loose zone of an extra-large section tunnel.
3.1. Orthogonal Analysis
The orthogonal experimental design method, proposed by mathematical statistician R.A. Fisher in 1935, has been widely used in mathematics, economics, agriculture, aerospace, and other fields. Compared with the comprehensive test method, it has several advantages, such as fewer experiments, careful consideration of factors, etc. According to the orthogonality, representative test points are selected from all level combinations of all test factors. These test points have uniform dispersion and neat comparability.
For example, if an event has three factors (
), each factor has three levels (
). A comprehensive test must be designed for 3
3 = 27 combination tests, while only nine tests are needed to analyze the degree of influence of each factor (
) via the orthogonal table
. The distribution of the comprehensive and orthogonal test points is shown in
Figure 6.
3.2. Numerical Models and Simulation Scheme
3.2.1. Hoek–Brown (H-B) Criterion
The H-B criterion, which perfectly reflects the non-linear failure characteristics of rocks and rock masses, was proposed by Hoek and Brown in 1980. After years of continuous development, it has become one of the most accepted criteria in the research on rock mass stability analysis and strength prediction [
31]. The expression of the criterion is as follows:
where
and
are the maximum and minimum principal stress of the rock mass, respectively;
is the uniaxial compressive strength of the rock specimen;
,
s, and
a are the material parameters of the rock mass related to lithology and the structural plane, which can be determined by the correlation:
where
is the material constant reflecting the softness and hardness of the rock mass.
D is the parameter reflecting the disturbance of rock mass, and its value ranges from 0 to 1.
3.2.2. Basic Assumptions
Numerical tests were conducted using the FLAC3D software. The following assumptions were made: (1) the models were established based on the plain strain hypothesis; (2) the surrounding rock mass was regarded as an ideal elastic–plastic body, following the H-B failure criterion; (3) the tectonic stresses were neglected in the numerical simulations; (4) the function of the supporting structure was not considered; (5) the model adopted the full-face method without considering the construction process and groundwater.
3.2.3. Simulation Scheme
Based on the traditional research method of loose zones in a surrounding rock mass, the following influencing factors were determined as the primary research objects: tunnel burial depth (H), tunnel span (B), tunnel height–span ratio (λ), lateral pressure coefficient (K), and GSI. There were five influencing factors in total, and each factor had four levels. This paper established 16 numerical models to carry out the orthogonal test in line with the orthogonal table .
According to the geological conditions and cross-sectional dimensions of the large-span section of the bifurcation,
H was 80~110 m,
B was 20~35 m,
λ was 0.6~1.2,
K is 0.5~2.0, and the GSI was 25~40. The orthogonal table and the rock mass mechanical parameters corresponding to GSI are shown in
Table 3 and
Table 4, respectively. The bulk density and Poisson’s ratio of the rock mass were identified by geological exploration reports.
value was mainly obtained from a laboratory test, and
and
can be evaluated via Correlation 2. The GSI-based empirical equations [
32,
33] used to determine the elastic modulus
E and
are proposed as follows:
The numerical calculation model of the unit thickness was set as: height × width = 180 m × 300 m (210 m × 300 m, 240 m × 300 m, 270 m × 300 m). The element mesh and model dimension are shown in
Figure 7.
H represents the vertical distance from the tunnel vault to the ground. The model’s boundary conditions were set so that the upper boundary was a free surface boundary, and both sides and the bottom edges were constrained by the method of normal displacement. The initial in situ stress was the self-weight stress of the rock mass.
3.2.4. Method of Determining the Loose Zone in the Surrounding Rock Mass
The pressure arch in the surrounding rock mass was characterized by an arch structure. The pressure arch body bears its own load and the load of the rock mass above the arch so that the surrounding rock maintains a certain degree of self-stability after excavation. However, the rock mass below the pressure arch has poor stability because it is close to the tunnel excavation outline, which leads to the stress exceeding the rock mass strength and causing damage. Therefore, the loose zone in the surrounding rock mass is defined as the rock mass from the inner boundary of the pressure arch to the tunnel contour. The principal stress vector diagram after the excavation of the No. 16 model is in
Figure 8, showing the pressure arch shape and the ring-shaped distribution of the maximum principal stress.
Based on existing research results [
34], this paper took the connection of the peak point of the surrounding rock’s tangential stress after the excavation of the tunnel as the judgment basis for the pressure arch inner boundary. The inner boundary criterion was embedded into FLAC3D with the Fish programming language to conduct numerical simulations. The correlation of the tangential stress can be calculated as follows:
The simulation was processed as follows: ① The model was established according to the orthogonal test scheme and the initial stress field was calculated. ② The element stress component was extracted after the calculation was balanced and the tangential stress of each component was calculated using the correlation. ③ Values of the tangential stress of each element on the search path and were compared, and the element with the peak tangential stress point was marked. Owing to the model’s symmetry, half of the model was selected for the search. The specific shape of the loose zone was obtained by connecting the inner boundary position of the pressure arch on each search path. ④ The loose zone was smoothed by the least square method to obtain the final scope of the loose zone in the surrounding rock mass, as shown in
Figure 9. It should be noted that the influence of floor heave was not considered in this study.
3.2.5. Analysis of the Factors Influencing the Loose Zone in the Surrounding Rock Mass
Three geometric quantities were selected as evaluation indexes reflecting the scope of the closed annular loose zone, namely
(the span of the loose zone at the tunnel vault),
(the height of the loose zone at the tunnel vault), and
(the area of the loose zone). The geometric quantity diagram of the loose zone is shown in
Figure 10. The data were processed for the 16 groups of working conditions, and the test number of the closed annular loose zone and the evaluation indexes of the loose zone scope were then summarized. The evaluation indexes of the loose area are shown in
Table 5.
The CCR model quantitatively evaluates the factors influencing the loose zone in the surrounding rock mass using a fuzzy restriction on the weight [
35]. This model can realize the evaluation and prediction of multi-index input and multi-index output problems with the help of the linear programming method. The model overcomes the shortcomings of the traditional CCR model, which discards influencing factors (the weight value is 0) when solving for optimal efficiency. The maximum fuzzy membership function constrains the upper and lower limits of the weight values to obtain the model’s optimal weight solution. The CCR model with a fuzzy restriction on the weight can be expressed as:
where
μ is the maximum fuzzy membership function;
is the output value of the
r-th term;
is the input value of the
i-th term;
is the weight value of the
r-th output;
is the weight value of the
i-th input;
and
are the upper limit values of the output weight and input weight, respectively.
When analyzing the factors influencing the loose pressure in the surrounding rock mass of a super-span tunnel, the input factors are
H,
B,
λ,
K, and GSI. The values of each impact factor are shown in
Table 2. The output factors are
L,
, and
S, and their geometric statistics are shown in
Table 6.
The weight values of the influencing factors are summarized in
Table 1. It can be concluded that the order of the weight values of the influencing factors in the loose zone of a super-large section tunnel is
, which indicates that two factors (
B for the tunnel span and GSI for surrounding rock quality) have a high degree of influence on the surrounding rock loose zone, while the other factors have less of an impact.
4. Boundary Function of the Surrounding Rock’s Loose Zone
The boundary curve of the loose zone was fitted by a non-linear function based on the scope of the loose zone in the surrounding rock mass. Owing to the symmetry of the loose zone in the surrounding rock, half of the boundary curve was taken as the research object. The Cartesian rectangular coordinate system was adopted, taking the arch vertex of the loose zone as the coordinate origin. The positive direction on the X-axis was vertically downward, and the positive direction on the Y-axis was horizontal and to the left. Ten characteristic points were selected from the coordinate information on the boundary curve and were fitted by the built-in function in MATLAB according to the changing trend of data points to determine the boundary function of the loose zone. After data sorting and analysis, the quadratic function had the highest degree of fit for the dataset. Therefore, this paper took it as the boundary function with
a,
b, and
c being undetermined coefficients. Taking Model 7 as an example, the fitting effect and relevant information are represented in
Figure 11 and
Table 7.
The following can be concluded from the correlation between the evaluation indexes of the loose zone and the boundary function:
Taking Correlation 5 into the boundary function
, the undetermined coefficients
a,
b, and
c can be approximated:
Using the data in
Table 3,
Table 4 and
Table 5, the multiple regression equations of
and
can be concluded via the process described below:
(1) Based on the principle of the least squares method, the univariate regression equation between the influencing factors of the loose zone with
and
are established, as shown in
Table 8. According to the manual, the correlation coefficient threshold value is 0.83. As illustrated in
Table 8, the correlation coefficients of the five influencing factors all exceed 0.83, indicating that the univariate regression equations of
and
are correct.
(2) The prediction correlations for
and
can be deduced from the univariate regression equations as below:
The correctness of Correlation 7 and Correlation 8 are verified by the variance analysis, as shown in
Table 9. The
F distribution value is
in Correlation 7 and
in Correlation 8. Both are greater than
at the 0.05 significance level, indicating that the accuracy of the regression equation meets the requirements.
It should be noted that the prediction correlations of and were established with , , , , . The parameter values must be defined within the domain.
5. Upper Bound Solution for Loose Pressure in the Surrounding Rock Mass
5.1. Upper Boundary Solution for Loose Pressure for Extra-Large Section Tunnels
5.1.1. Basic Assumptions
Four basic assumptions must be made when the upper bound theory is applied to analyze and calculate the loose pressure of surrounding rock:
① The rock–soil mass is regarded as a perfectly elastic–plastic material that obeys the H-B yield criteria.
② The tunnel arch is under a vertical even load q, and the side walls are under a horizontal even load e. Assuming that q is proportional to e, the proportional coefficient is the side pressure coefficient K, that is,
③ The influences of the construction process and floor heave are not considered.
5.1.2. Description of Failure Mechanism
The basic shape of the loose zone obtained from the orthogonal test shows that the failure forms of the deeply buried tunnel include the tunnel vault subsidence and side wall shrinkage. Based on the determination of the failure zone, the mass of rock and soil in the failure mechanism is divided into rigid blocks that meet the conditions. The velocity and direction of each rigid block are obtained according to the geometric relationship and the sine theorem corresponding to the failure mechanism.
In summary, the failure mechanism of a deeply buried tunnel is represented in
Figure 12. The tunnel section is approximately rectangular,
ABCD, with height
h and span
B.
M and
N are located at the midpoint of the tunnel vault and the midpoint of the loose zone boundary curve, respectively.
O is located on the straight line
MN and is the auxiliary point.
H represents the burial depth of the cavity
represents the distance from
O to the ground plane, and the support reaction forces of the tunnel arch and the side walls are the uniform loads
q and
e, respectively.
Owing to the symmetry of the failure mechanism, this study takes the left half of the failure mechanism for the corresponding description. is extended to intersect the boundary curve at point . is divided into parts with angles for . is divided into n parts with angles . The intersections of each angle and the boundary of the loose zone are recorded as . Connecting , when m and n are large enough, the error between each connecting line and the boundary curve is therefore small enough.
According to the associated flow rule, the velocity vector sum of the rigid slider is zero, and the angle between the relative velocity vector and the discontinuity line should be
.
Figure 13 shows the kinematically admissible velocity field corresponding to the failure mechanism. The velocity of block
is
in the vertically downward direction. The velocity of block
is
in the vertically downward direction. The velocities of the other blocks, from top to bottom, are
and the relative velocities of each rigid slider on the discontinuity line are
. The shape of the blocks is determined by
and
. With reference to Protodyakonov’s theory and the Terzaghi theory, the shapes can be assumed to be
.
5.1.3. Derivation of the Upper Bound Solution for the Loose Pressure
(1) Limit Analysis Based on H-B Strength Criterion
Based on the assumptions mentioned above, the loose zone is regarded as an ideal elastic–plastic body which obeys the H-B criterion. The relation between
and
can be obtained as below using the tangent method (
Figure 14) [
36]:
The non-linear shear strength index in the above correlations serves as an unknown parameter in the upper-bound limit analysis, reflecting the instantaneous variability of the tangent, which is obtained by the optimization algorithm. is obtained by Correlation 10 after is determined.
(2) Calculation of geometric quantities
Owing to the symmetry of the failure mechanism, this paper takes the left part of MN as an example for geometric quantity calculation. The point N is taken as the coordinate origin and a rectangular coordinate system is established. The vertical is the X-axis, where the down is positive, and the horizontal direction is the Y-axis, where the left is positive. According to the loose zone boundary equation, the expression of the curve NC is .
According to the relevant description of failure mechanism and the research results of the loose zone in the surrounding rock mass, it can be known that , , , , , .
It can be concluded from the trigonometric function:
The linear equations of
can be obtained from the coordinate
of point
and the division of
as follows:
The linear equations of
can be obtained from the coordinate
of point
and the division of
as follows:
The coordinates
of
can be obtained by combining the loose zone boundary equation with the linear equations, and the following correlations calculate the distance of each line segment:
can be calculated by the sine theorem:
(3) Calculation of the velocity field
As seen in
Figure 13, the velocity vector relationship of each slider meets the following recurrence correlations:
(4) Calculation of external powers
① Power of gravitational force
Gravitational power generated by each rigid slider:
② Power of Support Reaction
(5) Power of Internal Energy Dissipation
Power of Total Internal Energy Dissipation:
(6) Calculation of support reaction
Based on the virtual work principle, the power produced by the external force is equal to the power dissipated by the internal energy.
The correlation of the loose pressure in the surrounding rock mass is:
To ensure the velocity vector closure, each angle in the failure mechanism must meet the following constraints [
37]:
With the commands of the MATLAB optimization toolbox utilized and the constraints met, a certain number of the loose pressure values in the surrounding rock mass are obtained in which the maximum value is the upper bound solution for the loose pressure in the surrounding rock mass.
The horizontal support reaction can be calculated as:
5.2. Upper Boundary Solution for Loose Pressure for Extra-Large Section Tunnels
Based on the above research results, the process of calculating the loose pressure for a super-span tunnel is represented in
Figure 15. The loose pressure of the large-span section (the maximum Section, Gradient Section 1, and Gradient Section 2) was calculated according to the calculation process outlined above, with the geological conditions of the tunnel taken into consideration. The calculation parameters and the calculation results are represented in
Table 10 and
Table 11, respectively.
5.3. Comparison between the Correlation Results and In Situ Measured Values
Several difficulties and risks existed in the construction of the super-large section. In order to ensure the construction quality and grasp the distribution pattern of the surrounding rock pressure, our research group conducted a half-year on-site monitoring of the large-span section from April 2018 to October 2018. The monitoring items included the displacement of the deep surrounding rock, the contact pressure between the surrounding rock and the structure’s steel, and the stress of steel arch. The Maximum Section, the Gradient Section 1 and the Gradient Section 2 were selected as surrounding rock pressure monitoring sections. The site photos of the earth pressure cell are represented in
Figure 16.
The stable values of the loose pressure of each monitoring section’s surrounding rock were collected to determine the circumferential distribution. They were then compared with the value calculated for the loose pressure in the surrounding rock mass based on the upper bound limit analysis method, as shown in
Figure 17. The figure shows that the distribution law of the surrounding rock’s calculated loose pressure is consistent with the measured results except for the right side of the maximum section. The main reason for monitoring data distortion at the right side of the maximum section was the failure to install pressure sensors in a timely manner after excavation. As the surrounding rock pressure is monitored after the completion of the supporting structure, the monitoring results are the variation relative to the initial stress, and the monitoring values are usually less than the theoretical calculation values. As illustrated in
Figure 17, the calculated results have a good envelope effect on the measured data, which proves the correctness and practicability of the calculation method.
6. Conclusions
(1) With the combination of a numerical analysis and an orthogonal test, the factors affecting the loose zone range of the surrounding rock were studied. According to CCR model with fuzzy restrictions on the weight in statistics, the order of factors influencing the weight value from large to small is .
(2) The loose boundary zone expression of a deeply buried tunnel was determined to be , and a calculation method for the undetermined coefficients of the boundary expression was proposed. The correlations of typical geometric quantities and , representing the loose zone range, were established based on the unitary regression equations between the factors affecting the loose zone and typical geometric quantities.
(3) Based on the boundary function of the loose zone and the upper bound limit analysis method, a limit analysis model of the loose pressure in the surrounding rock mass under a non-linear failure criterion was established via the tangent method and virtual power principle, and the correlations of the loose pressure were deduced. The calculated values of the loose pressure were consistent with the pressure-monitoring results in the different sections of the large-span section, verifying the correctness of the proposed calculation method. The research results can provide new ideas for the selection of a supporting structure for similar super-section tunnel projects.
7. Discussion
(1) The orthogonal test designed in this paper considers several influence factors, and the value range is small (e.g., H is only 80~110 m.). In follow-up research, parameters such as the structural plane of the rock mass will be included in the influencing factors to improve the applicability of the calculation method.
(2) To simplify the analysis process, heterogeneities within the stratum were not considered in the study.
(3) The failure mechanism proposed in this paper did not consider the situation of floor heave, which must be improved in future research.
Author Contributions
Conceptualization, H.G. and X.L.; methodology, W.M.; software, X.L.; validation, H.G.; formal analysis, W.Z.; investigation, J.Z.; resources, W.M.; data curation, J.H.; writing—original draft preparation, H.G.; writing—review and editing, J.H.; visualization, X.L.; supervision, W.M.; project administration, H.G.; funding acquisition, W.M. All authors have read and agreed to the published version of the manuscript.
Funding
This research was funded by China National Railway Group Science and Technology Research Program grant number K2020G034.
Institutional Review Board Statement
Not applicable.
Informed Consent Statement
Not applicable.
Data Availability Statement
The simulation and monitoring data in the article are not freely available due to legal concerns and commercial confidentiality. Nevertheless, all the concepts and procedures are explained in the presented research and parts of the research may be available upon request.
Acknowledgments
The authors of the present work feel grateful and would like to thank the China State Railway Group Co., Ltd., Beijing General Municipal Engineering Design & Research Institute Co., Ltd., and the China Academy of Railway Sciences Co., Ltd. for providing facilities, funds, and support in measurement and monitoring.
Conflicts of Interest
The authors declare no conflict of interest.
References
- Gu, Z.Q.; Peng, S.Z.; Li, Z.K. Underground Cavern Engineering; Tsinghua University Press: Beijing, China, 1994; pp. 50–53. [Google Scholar]
- Bierbaumer, A. Diedimensionerung des Tunnelmauerwerks; Engelmann: Leipzig, Germany, 1913. [Google Scholar]
- Terzaghi, K. Theoretical Soil Mechanics; John Wiley and Sons, Inc.: Hoboken, NJ, USA, 1943. [Google Scholar]
- Xie, J. Earth Pressure on Shallow Burial Tunnel. China Civil Eng. J. 1964, 6, 58–70. [Google Scholar]
- TB10003; Code for Design of Railway Tunnels. China Railway Publishing House: Beijing, China, 2005.
- JTG D70-2010; Code for Design of Road Tunnels. People’s Communications Publishing House: Beijing, China, 2004.
- SL279; Code for Design of Hydraulic Tunnels. China Water & Power Press: Beijing, China, 2002.
- Chongqing Civil Engineering Institute. Underground Construction in Rock; China Architecture & Building Press: Beijing, China, 1982. [Google Scholar]
- Barton, N.; Lien, R.; Lunde, J. Engineering classification of rock masses for the design of tunnel support. Rock Mech. 1974, 6, 183–236. [Google Scholar] [CrossRef]
- Bhasin, R.; Grimstad, E. The Use of Stress-Strengh Relationship in the Assessment of Tunnel Stability. Proc. Recent Adv. Tunn. Technol. 1996, 11, 93–98. [Google Scholar]
- Unal, E. Development of Design Guidelines and Roof-Control Standards for Coal-Mine Roofs; Pennsylvania State University: State College, PA, USA, 1983. [Google Scholar]
- Goel, R.K.; Jethwa, J.L. Prediction of Support Pressure using RMR Classification. In Proceedings of the Indian Geotechnical Conference, Surat, India, 20–22 December 1991. [Google Scholar]
- Atkinson, J.H.; Potts, D.M. Stability of a shallow circular tunnel in cohesion-less soil. Geotechnique 1977, 27, 203–215. [Google Scholar] [CrossRef]
- Davis, E.H.; Gunn, M.J.; Mair, R.J.; Seneviratine, H.N. The stability of shallow tunnels and underground openings in cohesive material. Geotechnique 1980, 30, 397–416. [Google Scholar] [CrossRef]
- Mollon, G.; Dias, D.; Soubra, A.H. Face stability analysis of circular tunnels driven by a pressurized shield. J. Geotech. Geoenvironmental Eng. 2010, 136, 215–229. [Google Scholar] [CrossRef]
- Mollon, G.; Dias, D.; Soubra, A.H. Rotational failure mechanisms for the face stability analysis of tunnels driven by a pressurized shield. Int. J. Numer. Anal. Methods Geomech. 2011, 35, 1363–1388. [Google Scholar] [CrossRef]
- Huang, F.; Yang, X.L. Upper bound limit analysis of collapse shape for circular tunnel subjected to pore pressure based on the Hoek–Brown failure criterion. Tunn. Undergr. Space Technol. Inc. Trenchless Technol. Res. 2011, 26, 614–618. [Google Scholar] [CrossRef]
- Lei, M.; Peng, L.; Shi, C. Calculation of the surrounding rock pressure on a shallow buried tunnel using linear and nonlinear failure criteria. Autom. Constr. 2014, 37, 191–195. [Google Scholar] [CrossRef]
- Pan, Q.; Dias, D. Face stability analysis for a shield-driven tunnel in anisotropic and nonhomogeneous soils by the kinematical approach. Int. J. Geomech. 2016, 16, 04015076. [Google Scholar] [CrossRef]
- Pan, Q.; Dias, D. The effect of pore water pressure on tunnel face stability. Int. J. Numer. Anal. Methods Geomech. 2016, 40, 2123–2136. [Google Scholar] [CrossRef]
- Pan, Q.; Dias, D. Upper-bound analysis on the face stability of a non-circular tunnel. Tunn. Undergr. Space Technol. 2017, 62, 96–102. [Google Scholar] [CrossRef]
- Zhang, F.; Gao, Y.; Wu, Y.; Wang, Z. Face stability analysis of large-diameter slurry shield-driven tunnels with linearly increasing undrained strength. Tunn. Undergr. Space Technol. 2018, 78, 178–187. [Google Scholar] [CrossRef]
- Mollon, G.; Dias, D.; Soubra, A.H. Continuous velocity fields for collapse and blowout of a pressurized tunnel face in purely cohesive soil. Int. J. Numer. Anal. Methods Geomech. 2013, 37, 2061–2083. [Google Scholar] [CrossRef]
- Huang, M.; Tang, Z.; Zhou, W.; Yuan, J. Upper bound solutions for face stability of circular tunnels in non-homogeneous and anisotropic clays. Comput. Geotech. 2018, 98, 189–196. [Google Scholar] [CrossRef]
- Tu, S.; Li, W.; Zhang, C.; Chen, W. Effect of inclined layered soils on face stability in shield tunneling based on limit analysis. Tunn. Undergr. Space Technol. 2023, 131, 104773. [Google Scholar] [CrossRef]
- Wang, H.T.; Liu, P.; Liu, C.; Zhang, X.; Yang, Y.; Liu, L.-Y. Three-dimensional upper bound limit analysis on the collapse of shallow soil tunnels considering roof stratification and pore water pressure. Math. Probl. Eng. 2019, 2019, 8164702. [Google Scholar] [CrossRef]
- Wang, H.; Liu, P.; Wang, L.G.; Liu, C.; Zhang, X.; Liu, L. Three-dimensional collapse analysis for a shallow cavity in layered strata based on upper bound theorem. Comput. Model. Eng. Sci. 2020, 124, 375–391. [Google Scholar] [CrossRef]
- Abate, G.; Grasso, S.; Massimino, M.R. The role of shear wave velocity and non-linearity of soil in the seismic response of a coupled tunnel-soil-above ground building system. Geosciences 2019, 9, 473. [Google Scholar] [CrossRef]
- Castelli, F.; Grasso, S.; Lentini, V.; Sammito, M.S.V. Effects of soil-foundation-interaction on the seismic response of a cooling tower by 3D-FEM analysis. Geosciences 2021, 11, 200. [Google Scholar] [CrossRef]
- Gao, H.; He, P.; Chen, Z.; Li, X. A novel calculation method of process load for extra-large section tunnels. Symmetry 2019, 11, 1228. [Google Scholar] [CrossRef]
- Hoek, E.; Brown, E.T. Empirical strength criterion for rock masses. ASCE J. Geotech. Eng. 1980, 106, 1013–1036. [Google Scholar] [CrossRef]
- Hoek, E.; Diederichs, M.S. Empirical estimation of rock mass modulus. Int. J. Rock Mech. Min. Sci. 2006, 43, 203–215. [Google Scholar] [CrossRef]
- Cui, L.; Sheng, Q.; Zheng, J.-J.; Cui, Z.; Wang, A. Regression model for predicting tunnel strain in strain-softening rock mass for underground openings. Int. J. Rock Mech. Min. Sci. 2019, 119, 81–97. [Google Scholar] [CrossRef]
- Yang, J.H.; Wang, S.R.; Wang, Y.G.; Li, C.L. Analysis of arching mechanism and evolution characteristics of tunnel pressure arch. Jordan J. Civ. Eng. 2015, 9, 125–132. [Google Scholar]
- Wang, S.; Ma, Q.H.; Guan, Z.M. Application of DEA Model with Fuzzy Restriction on Weights. Ind. Eng. Manag. 2006, 5, 103–106. [Google Scholar]
- Yang, X.L.; Li, L.; Yin, J.H. Stability analysis of rock slopes with a modified Hoek-Brown failure criterion. Int. J. Numer. Anal. Methods Geomech. 2004, 28, 181–190. [Google Scholar] [CrossRef]
- Chen, W.F. Limit Analysis and Soil Plasticity; Elsevier Science: Amsterdam, The Netherlands, 2013. [Google Scholar]
Figure 1.
Empirical formula classification.
Figure 1.
Empirical formula classification.
Figure 2.
Route planning of the Greater Bay Area.
Figure 2.
Route planning of the Greater Bay Area.
Figure 3.
Plan view of bifurcation section.
Figure 3.
Plan view of bifurcation section.
Figure 4.
Engineering geologic profile of the project.
Figure 4.
Engineering geologic profile of the project.
Figure 5.
Site photos of the drill core. (a) Small clear-distance section; (b) large-span section.
Figure 5.
Site photos of the drill core. (a) Small clear-distance section; (b) large-span section.
Figure 6.
Distribution of test points. (a) Comprehensive test; (b) orthogonal test.
Figure 6.
Distribution of test points. (a) Comprehensive test; (b) orthogonal test.
Figure 7.
Numerical model.
Figure 7.
Numerical model.
Figure 8.
Maximum principal stress vector cloud diagram (No. 16).
Figure 8.
Maximum principal stress vector cloud diagram (No. 16).
Figure 9.
Search path and loose zone boundary.
Figure 9.
Search path and loose zone boundary.
Figure 10.
Loose zone of surrounding rock mass.
Figure 10.
Loose zone of surrounding rock mass.
Figure 11.
The fitting curve of the surrounding rock’s loose zone.
Figure 11.
The fitting curve of the surrounding rock’s loose zone.
Figure 12.
Failure mechanisms.
Figure 12.
Failure mechanisms.
Figure 13.
The velocity field schematic diagram.
Figure 13.
The velocity field schematic diagram.
Figure 14.
The tangential method.
Figure 14.
The tangential method.
Figure 15.
Flow of calculation for the loose pressure in the surrounding rock mass.
Figure 15.
Flow of calculation for the loose pressure in the surrounding rock mass.
Figure 16.
Installation of the earth pressure cell. (a) Arch crown; (b) haunch.
Figure 16.
Installation of the earth pressure cell. (a) Arch crown; (b) haunch.
Figure 17.
Comparison between the calculated and measured values.
Figure 17.
Comparison between the calculated and measured values.
Table 1.
Different correlations for calculation methods.
Table 1.
Different correlations for calculation methods.
Item | Correlation | Applicable Conditions | Weakness |
---|
Protodyakonov’s theory [1] | | Deep burial condition; small span | Competent coefficient of rock is determined subjectively |
Terzaghi theory [3] | | Shallow burial condition; loose media | Unsuitable for rock stratum |
Bierbaumer theory [2] | | Shallow burial condition; loose media | Unsuitable for rock stratum |
Xie Jiaxiao’s theory [4] | | Shallow burial condition; small span | Unsuitable for deep burial condition |
Code for the Design of Railway Tunnel (TB10003-2016) [5] | | Small span | Unsuitable for conditions with span greater than 15 m |
Specification for the Design of Hydraulic Tunnels (SL279-2016) [7] | | Specific stratum | Neglects the effect of surrounding rock conditions |
Design Code of Artificial Rock Caverns (in Chinese) [8] | | Specific stratum | Neglects the effect of surrounding rock conditions |
Q System | Barton correlation [9] | | Small span | Neglects the effect of geometrical dimension |
Bhasin and Grimstad correlation [10] | | Broken rock condition | Neglects the effect of height |
RMR System | Unal correlation [11] | | Specific stratum | Neglects the effect of buried depth |
Goel and Jethwa correlation [12] | | Deep burial condition | Neglects the effect of height |
Table 2.
Characteristics of the surrounding rock of the project.
Table 2.
Characteristics of the surrounding rock of the project.
Mileage | Length | Surrounding Rock Characteristics | GSI |
---|
0~130 | 130 m | The surrounding rock is slightly weathered sandstone, and the rock mass is fractured with developed joint fissures. Generally, the groundwater is poor, and the permeability is low. | 40 |
130~244 | 114 m | The surrounding rock is a greenish-gray, medium, and slightly weathered sandstone. The rock mass is relatively fractured with developed joint fissures, and the rock core is short, columnar, and massive. The groundwater is poor and does not have the conditions for water inrush. | 28 |
Table 3.
Orthogonal test table.
Table 3.
Orthogonal test table.
Model Number | Buried Depth H/m | Tunnel Span B/m | Height–Span Ratio λ | Lateral Pressure Coefficient K | GSI |
---|
1 | 80 | 20 | 0.6 | 0.5 | 40 |
2 | 80 | 25 | 0.8 | 1.0 | 35 |
3 | 80 | 30 | 1.0 | 1.5 | 30 |
4 | 80 | 35 | 1.2 | 2.0 | 25 |
5 | 90 | 20 | 0.8 | 1.5 | 25 |
6 | 90 | 25 | 0.6 | 2.0 | 30 |
7 | 90 | 30 | 1.2 | 0.5 | 35 |
8 | 90 | 35 | 1.0 | 1.0 | 40 |
9 | 100 | 20 | 1.0 | 2.0 | 35 |
10 | 100 | 25 | 1.2 | 1.5 | 40 |
11 | 100 | 30 | 0.6 | 1.0 | 25 |
12 | 100 | 35 | 0.8 | 0.5 | 30 |
13 | 110 | 20 | 1.2 | 1.0 | 30 |
14 | 110 | 25 | 1.0 | 0.5 | 25 |
15 | 110 | 30 | 0.8 | 2.0 | 40 |
16 | 110 | 35 | 0.6 | 1.5 | 35 |
Table 4.
Mechanical parameters of the rock mass corresponding to different GSI values.
Table 4.
Mechanical parameters of the rock mass corresponding to different GSI values.
GSI | (MPa) | Bulk Density (kN/m3) | Elastic Modulus (GPa) | Poisson’s Ratio | | | |
---|
40 | 70.65 | 23 | 4.73 | 0.28 | 12.11 | 1.136 | 1.018 |
35 | 59.28 | 22 | 3.25 | 0.31 | 10.94 | 0.859 | 0.584 |
30 | 48.89 | 21 | 2.21 | 0.32 | 9.73 | 0.639 | 0.335 |
25 | 39.67 | 20 | 1.49 | 0.34 | 8.48 | 0.466 | 0.192 |
Table 5.
Summary of evaluation indexes of the loose area.
Table 5.
Summary of evaluation indexes of the loose area.
Model Number | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 |
Evaluation index | | 30.46 | 30.56 | 48.34 | N/A | 50.90 | 41.34 | 25.46 | 35.70 |
| 17.25 | 9.77 | 11.19 | N/A | 7.04 | 11.30 | 3.11 | 12.67 |
| 239 | 207 | 336 | N/A | 470 | 203 | 104 | 228 |
Model Number | 9 | 10 | 11 | 12 | 13 | 14 | 15 | 16 |
Evaluation index | | 41.02 | 50.50 | N/A | 61.58 | 45.30 | 58.50 | 39.78 | 65.46 |
| 23.74 | 28.17 | N/A | 24.38 | 25.33 | 9.81 | 17.23 | 11.20 |
| 517 | 597 | N/A | 887 | 552 | 467 | 392 | 772 |
Table 6.
Weight value of influencing factors.
Table 6.
Weight value of influencing factors.
Influencing Factors | H | B | λ | K | GSI |
---|
Weight value | 0.12 | 0.37 | 0.15 | 0.14 | 0.22 |
Table 7.
Results for the fitting curve.
Table 7.
Results for the fitting curve.
Model Number | Formula | Undetermined Coefficient | Correlation Coefficient R2 |
---|
a | b | c |
---|
7 | . | −0.1174 | 2.5181 | −0.1086 | 0.9988 |
Table 8.
Unitary regression equation.
Table 8.
Unitary regression equation.
Factors | Formula | Correlation Coefficient R2 |
---|
H | | 0.9986 |
| 0.9433 |
B | | 0.8477 |
B | 0.8535 |
λ | | 0.9767 |
| 0.9978 |
K | | 0.9991 |
| 0.8659 |
GSI | | 0.8699 |
| 0.8904 |
Table 9.
Variance analysis of multiple regression equation.
Table 9.
Variance analysis of multiple regression equation.
Variance Analysis of Correlation 7 | Variance Analysis of Correlation 8 |
---|
Source of Variance | Sum of Deviation Squares | Freedom | Mean Square | F | Source of Variance | Sum of Deviation Squares | Freedom | Mean Square | F |
---|
Regression | 21391.7 | 6 | 1572.8 | 6.59 | Regression | 23964.4 | 6 | 1160.0 | 7.02 |
Error | 3149.1 | 4 | 251.3 | N/A | Error | 3645.5 | 4 | 269.7 | N/A |
sum | 24646.6 | 10 | N/A | N/A | sum | 29764.1 | 10 | N/A | N/A |
Table 10.
Calculation parameters.
Table 10.
Calculation parameters.
Section Size | GSI | H (m) | B (m) | λ | K |
---|
Maximum section | 40 | 104 | 30.01 | 0.61 | 0.50 |
Gradient Section 1 | 40 | 90 | 27.16 | 0.65 | 0.50 |
Gradient Section 2 | 28 | 82 | 23.62 | 0.70 | 0.50 |
Table 11.
Loose pressure in the surrounding rock mass based on the upper bound analysis.
Table 11.
Loose pressure in the surrounding rock mass based on the upper bound analysis.
Section Size | Evaluation Index | Boundary Function | Surrounding Rock Loose Pressure |
---|
L/m | H1/m | Vault Vertical Pressure (kPa) | Haunch Horizontal Pressure (kPa) |
---|
Maximum Section | 40.54 | 9.8 | | 269.44 | 134.72 |
Gradient Section 1 | 36.72 | 9.2 | | 243.36 | 121.68 |
Gradient Section 2 | 41.26 | 12.6 | | 307.11 | 153.56 |
| Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).