Elastic Stress and Plastic Zone Distributions around a Deeply Buried Tunnel under the Nonhydrostatic Pressure

Liu, Shuhong; Fan, Jinlu; Wu, Chenhao; Zhu, Yongquan

doi:https://doi.org/10.1155/2022/2363291

Geofluids

On this page

Abstract Introduction Results and Discussion Conclusions Data Availability Conflicts of Interest Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2022 | Article ID 2363291 | https://doi.org/10.1155/2022/2363291

Elastic Stress and Plastic Zone Distributions around a Deeply Buried Tunnel under the Nonhydrostatic Pressure

Shuhong Liu,¹Jinlu Fan,¹Chenhao Wu,¹and Yongquan Zhu²

Academic Editor: Qibin Lin

Received08 May 2022

Revised12 Jun 2022

Accepted21 Jun 2022

Published04 Jul 2022

Abstract

The elastic stress and the plastic zone are the important mechanical parameters to determine the tunnel support design. Based on Muskhelishvili’s complex variable function, the analytical solution for the elastic stress around a deeply buried noncircular tunnel under the nonhydrostatic pressure is firstly derived. The shape and size of the plastic zone of the surrounding rock mass are then determined by substituting the elastic stresses into the Drucker-Prager yield criterion. Finally, taking a horseshoe-shaped tunnel as an example, the analytical solutions of elastic stress distribution and plastic zone shape around the tunnel under different lateral pressure coefficients are in good agreement with ANSYS numerical solutions. The calculation results show that if the vertical in situ stress exceeds the critical value, with increasing lateral pressure coefficient, the hoop stress at the roof and floor of the tunnel increases significantly and the shape and size of the plastic zone change obviously.

1. Introduction

With the rapid development of modern roads, railways, and mining, the role of tunnels is becoming more and more important. The safety of tunnel construction has attracted attention in both academia and engineering circles. For the plane problem of a hole or a crack in an infinite body, the most effective method is using the complex variable theory. Closed-form solutions for simple shapes, such as an elliptical hole, a crack, and a square hole, have been discussed in detail. The exact solutions for cracks emanating from an elliptical hole were also studied [1–3]. This method has been successfully used to solve the tunnel problem. A deeply buried tunnel can be regarded as a hole in an infinite elastic body. For unlined tunnels, Exadaktylos et al. [4, 5] gave analytical solutions of stress and displacement distributions around the semicircular or notched circular tunnels and compared with numerical results of FLAC2D. Guan et al. [6] and Liu et al. [7, 8] studied horseshoe-shaped tunnels and compared analytical solutions with FLAC3D and ANSYS numerical results, respectively. For a lined straight wall arch tunnel, Kargar et al. [9, 10] and Lu et al. [11] provided analytical solutions for stress according to the different contact condition between lining and surrounding rock; Lu et al. [12] gave analytic solutions for stress and displacement considering support delay. Li and Chen [13] and Liu et al. [14] obtained the analytical solutions for stress and displacement of lined horseshoe-shaped tunnels in isotropic and orthotropic surrounding rock, respectively.

The plastic zone around a tunnel is induced simultaneously when tunnels are excavated in the surrounding rock mass. The extension of the plastic zone may increase pressure of the lining and thus inevitably reduce the safety of the tunnel. Analytical solutions for the elastic–plastic stresses and the plastic zone around a circular tunnel under an axial in situ stress were investigated [15–17]. While the surrounding rock is often under unequal biaxial in situ stresses in actual geological conditions, analytical solutions for the plastic zone in the vicinity of a circular tunnel were also obtained. Behnam et al. [18], Deng et al. [19], and Guo et al. [20] gave the shape of the plastic zone for different lateral pressure coefficient. Guo et al. [21] and Dong et al. [22] verified the correctness of the plastic zone predicted by theoretical calculation through numerical simulation.

Considering the pressure relief effect, Xu et al. [23] proposed a method of grooving slots in the roof and floor of a circular tunnel under the nonuniform stress; then, a circular tunnel actually becomes a circular tunnel with two cracks. Based on the different yield criterion, the shape and size of the crack tip plastic zone were analyzed [24–26]. Li et al. [27] obtained the numerical distribution of the plastic zone in the surrounding rock with a large-span rectangular roadway using FLAC^3D software. Using the Griffith yield criterion, Shi and Bai [28] studied the plastic zone between a horseshoe-shaped tunnel and a concealed cave to ensure rock safety thickness. Guo et al. [29] obtained the explicit form of analytical solution caused by a shallow circular tunnel and analyzed the plastic zone around the tunnel with pile load combined with Mindlin’s solution. Zou et al. [30] considered the effect of the intermediate principal stress and interaction between the surrounding rock and support structure on the plastic zone of the surrounding rock around the circular tunnel. Shi et al. [31] and Zhou and Wu [32] calculated the plastic zones of circular tunnel surrounding rock by different strength criteria. Ma et al. [33, 34] determined the plastic zones around two circular tunnels. Hu et al. [35] investigated the influence of left and right tunnels on the plastic zone shape of the surrounding rock around the center circular tunnel.

The above-described literature shows that many studies focus on analytical solutions for elastic stresses of the surrounding rock with deep-buried tunnels. Additionally, more results of the plastic zone around circular tunnels are presented. The analytical distribution of the plastic zone around noncircular tunnels in engineering is very few due to the complexity of mathematics.

In this paper, based on Muskhelishvili’s complex variable function theory, the analytical solution for the elastic stress of a deeply buried noncircular tunnel under the nonhydrostatic pressure is derived. Subsequently, the elastic stresses are taken into the Drucker-Prager yield criterion; then, the plastic zone can be estimated in the surrounding rock mass. Taking a horseshoe-shaped tunnel as an example, the analytical solution of elastic stress distributions and plastic zone shapes for different lateral pressure coefficients is given. The finite element software ANSYS is used to verify the accuracy of the proposed analytical solutions.

2. Problem and Basic Equations

The schematic diagram of the noncircular tunnel is presented in Figure 1. Assuming that the buried depth is great enough, the effect of gravity can be neglected. The surrounding rock mass is subjected to the nonhydrostatic pressure and acting along - and -axes, respectively. is the lateral pressure coefficient.

Muskhelishvili’s method has obvious advantages in solving plane problems involving holes with complex geometrical shapes. The mapping function transforms the outer region of the tunnel in the physical -plane with into the outer region of a unit circle in the mathematical -plane with , where and are polar coordinates in the plane. The general form of the conformal mapping function can be expressed as where the coefficients and are known constants for the given tunnel shape [36]. In general, only a few terms are sufficient to satisfy the accuracy of mapping.

For plane strain problems, stresses and displacements in the surrounding rock mass are determined by the following stress functions and : where denotes an arbitrary point of the unit circle. The constants , , and are relevant with far-field stresses and :

Substituting equations (1) and (7) into equation (6), can be rewritten as

The stress components , , and in Cartesian coordinates can be represented as where , .

The stress components and in curvilinear coordinates can be represented as

When there is no surface force on the boundary of the tunnel, , then, equation (10) is simplified as

3. Solution for Stress Functions

From equation (1), the following formula may be derived: where the coefficients are expressed by the known coefficients of equation (1).

The analytic function can be expressed as the following series form: where are unknown real coefficients obtained from the following boundary conditions.

Combining equations (12) and (14) yields where

The real coefficients are expressed by the known coefficients of equation (1) and the still unknown coefficients of equation (14). and are not given because they are not used in the following derivation.

Substituting equations (8), (14), and (15) into (4) and calculating the Cauchy integral on the right-hand side of this equation yield

In equation (17), it can be seen that . Because the coefficients of the same negative exponents of on both sides of equation (17) are equal, the following equations are obtained:

The coefficients solved by the equation (18) give the analytic function :

Substituting equations (1), (7), and (19) into (2), the first stress function is given as follows:

Next, another analytic function will be derived. The conjugate of complex numbers in equations (8) and (15) may be written:

Similarly, substituting equation (21) into (5) and calculating the Cauchy integral on the right-hand side of this equation, can be obtained:

Substituting equations (1), (7), and (22) into (3), then, the second stress function is given:

4. Plastic Zone

For the noncircular tunnel under different in situ stresses, the exact closed-form solution for the plastic zone in the surrounding rock mass cannot be obtained. A suitable yield criterion for the surrounding rock mass is introduced to estimate the plastic zone. In this section, the Drucker-Prager (D-P) yield criterion is adopted: where is the first stress invariant and is the second invariant of the deviator stress tensor. and are internal friction angle and adhesion of the surrounding rock mass, respectively. The following formulas may be written:

For plane strain problem, , and μ represents the Poisson’s ratio. Substituting equations (9) and (25) into (24), the relationship between and can be obtained; then, the shape and size of the plastic zone are determined from equation (1).

5. Results and Discussion

In order to verify the derived solution, analytical results of the elastic stress distributions and the plastic zone shapes around a horseshoe-shaped tunnel under different lateral pressure coefficients are compared with numerical results predicted by ANSYS finite element software.

The tunnel is a two-lane highway tunnel located in grade V surrounding rock. The horseshoe-shaped cross-section consists of six successive arcs, as shown in Figure 2. The center position and the radius of each arc are , , , and and , , , and , respectively.

The coefficients in equation (1) which was obtained by searching points on the boundary are as follows [36]: , , , , , , , , , and .

The tunnel cross-section after mapping is adopted for analytical and numerical analyses.

The two-dimensional finite element model of size is much larger than the tunnel. Therefore, the numerical model can be regarded as the infinite surrounding rock mass with a horseshoe-shaped tunnel under far-field loads and , just like the analytic model. The following parameters of the surrounding rock are assumed elastic modulus , Poisson’s ratio , adhesion , and internal friction angle .

5.1. Elastic Stress Distribution

In this section, far-field stress and the lateral pressure coefficient . Figure 3 displays the contours of elastic stresses , , and obtained by the ANSYS finite element software. Figures 4–6 illustrate the distributions of elastic stresses , , and along the tunnel boundary, the -axis and the -axis predicted by the proposed analytical solution and ANSYS finite element software. A positive angle turns counterclockwise from the positive -axis to the positive -axis (Figure 1). In these figures, it can be seen that there is obvious stress concentration around the tunnel. Far away from the tunnel, the stresses , , and tend to the constant values of −3 MPa, −2 MPa, and 0 MPa, which are the values of far-field in situ stresses. The excellent agreement of both predictions may be seen in Figures 4–6, because the same tunnel shape is adopted. These results indicate that the above-derived stress functions are accurate.

(a)

(b)

(c)

Figure 7 illustrates the distribution of hoop stress along the tunnel boundary predicted by the proposed analytical solution for , 0.4, 0.6, 0.8, 1.0, and 1.2 for the case of . An interesting phenomenon is found that there are two identical inflection points of the six curves. With increasing lateral pressure coefficient, the hoop stress at the points and 180° changes significantly.

5.2. Plastic Zone

Figures 8 and 9 illustrate the shape and size of the plastic zone around the tunnel predicted by ANSYS finite element software and the above-derived analytical solution for , 0.4, 0.6, 0.8, 1.0, and 1.2 for the case of . The good agreement of both predictions may be seen from these figures. The plastic zone distribution is symmetrical to the vertical axis, and the boundary of the plastic zone is smooth.

(a)

(b)

(c)

(d)

(e)

(f)

(a)

(b)

(c)

(d)

(e)

(f)

In order to view the variation of the shape and size of the plastic zone change with the lateral pressure coefficient, the six conditions predicted by the analytical solution are plotted in Figure 10. It can be seen that with the lateral pressure coefficient increases, the shape and size of the plastic zone change obviously and the location extends from the tunnel sides to the roof and floor in turn. However, the plastic zone around the tunnel sidewalls shows little change. These results indicate that the plastic zone estimation in Section 4 is accurate. In the study, we also find that if the vertical in situ stress is less than a certain critical value, such as , the location of the plastic zone varies with the lateral pressure coefficient.

6. Conclusions

For the plane strain problem of the surrounding rock mass with a noncircular tunnel under the nonhydrostatic pressure, the analytical solution for the elastic stresses was obtained. Hence, the shape and size of the plastic zone can be estimated. (1)The elastic stress distributions and the plastic zone shapes around the tunnel for different lateral pressure coefficients predicted by the proposed analytical model were in good agreement with those by the ANSYS numerical model(2)There is obvious stress concentration around the tunnel, and the stress values tend to the far-field in situ stresses far away from the tunnel. With increasing lateral pressure coefficient, the hoop stress at the roof and floor of the tunnel changes significantly. For different lateral pressure coefficients, there are two identical inflection points on the distribution curves of hoop stress along the tunnel boundary(3)If the vertical in situ stress exceeds the critical value, with increasing lateral pressure coefficient, the shape and size of the plastic zone change obviously and the location extends from the tunnel sides to the roof and floor in turn. However, the plastic zone around the tunnel sidewalls shows little change

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The work was financially supported by the Natural Science Foundation of China (51778380) and China State Railway Group Co. Ltd. (N2020G040).

References

S. Liu and S. Duan, “Analytical solutions of cracks emanating from an elliptical hole under shear,” Chinese Journal of Aeronautics, vol. 27, no. 4, pp. 829–834, 2014.
View at: Publisher Site | Google Scholar
S. Liu and S. Duan, “Analytical solutions of cracks emanating from an elliptic hole in an infinite plate under tension,” Chinese Journal of Mechanical Engineering, vol. 27, no. 5, pp. 1057–1063, 2014.
View at: Publisher Site | Google Scholar
M. Hajimohamadi and R. Ghajar, “An analytical solution for the stress field and stress intensity factor in an infinite plane containing an elliptical hole with two unequal aligned cracks,” Applied Mathematics and Mechanics (English Edition), vol. 39, no. 8, pp. 1103–1118, 2018.
View at: Publisher Site | Google Scholar
G. Exadaktylos and M. Stavropoulou, “A closed-form elastic solution for stresses and displacements around tunnels,” International Journal of Rock Mechanics and Mining Sciences, vol. 39, no. 7, pp. 905–916, 2002.
View at: Publisher Site | Google Scholar
G. Exadaktylos, P. Liolios, and M. Stavropoulou, “A semi-analytical elastic stress–displacement solution for notched circular openings in rocks,” International Journal of Solids and Structures, vol. 40, no. 5, pp. 1165–1187, 2003.
View at: Publisher Site | Google Scholar
Z. Guan, T. Deng, H. Huang, and Y. Jiang, “Back analysis technique for mountain tunneling based on the complex variable solution,” International Journal of Rock Mechanics and Mining Sciences, vol. 59, pp. 15–21, 2013.
View at: Publisher Site | Google Scholar
S. Liu, Z. Wang, and W. Zhen, “Complex variable function solutions for stress and displacement distributions in surrounding rock of deep buried highway tunnels,” Applied Mathematics and Mechanics, vol. 39, no. 5, pp. 548–557, 2018.
View at: Google Scholar
S. Liu, J. Fan, Z. Chen, and Y. Zhu, “Analytical solution to elastic and viscoelastic displacements for rock mass caused by excavation of deep-buried non-circular tunnel,” Journal of the China Railway Society, vol. 44, no. 3, pp. 142–149, 2022.
View at: Google Scholar
A. Kargar, R. Rahmannejad, and M. Hajabasi, “A semi-analytical elastic solution for stress field of lined non-circular tunnels at great depth using complex variable method,” International Journal of Solids and Structures, vol. 51, no. 6, pp. 1475–1482, 2014.
View at: Publisher Site | Google Scholar
A. Kargar, R. Rahmannejad, and M. Hajabasi, “The stress state around lined non-circular hydraulic tunnels below the water table using complex variable method,” International Journal of Rock Mechanics and Mining Sciences, vol. 78, pp. 207–216, 2015.
View at: Publisher Site | Google Scholar
A. Lu, N. Zhang, and Y. Qin, “Analytical solutions for the stress of a lined non-circular tunnel under full- slip contact conditions,” International Journal of Rock Mechanics and Mining Sciences, vol. 79, pp. 183–192, 2015.
View at: Publisher Site | Google Scholar
A. Lu, N. Zhang, and L. Kuang, “Analytic solutions of stress and displacement for a non-circular tunnel at great depth including support delay,” International Journal of Rock Mechanics and Mining Sciences, vol. 70, pp. 69–81, 2014.
View at: Publisher Site | Google Scholar
Y. Li and S. Chen, “Lining and surrounding rock in non-circular tunnel based on complex variable method,” Journal of Southwest Jiaotong University, vol. 55, no. 2, pp. 265–272, 2020.
View at: Google Scholar
S. Liu, P. Zhang, S. Chen, and Y. Zhu, “Stress and displacement of orthotropic surrounding rock with a deeply-buried non-circular lined tunnel under the complete contact interface condition,” Journal of Mechanical Engineering, vol. 57, no. 23, pp. 149–159, 2021.
View at: Publisher Site | Google Scholar
A. Lu, G. Xu, F. Sun, and W. Sun, “Elasto-plastic analysis of a circular tunnel including the effect of the axial in situ stress,” International Journal of Rock Mechanics and Mining Sciences, vol. 47, no. 1, pp. 50–59, 2010.
View at: Publisher Site | Google Scholar
X. Huang, J. Zhang, L. Yang, S. Yang, and X. Wang, “Elasto-plastic analysis of the surrounding rock mass in circular tunnel based on the generalized nonlinear unified strength theory,” International Journal of Mining Science and Technology, vol. 26, no. 5, pp. 819–823, 2016.
View at: Publisher Site | Google Scholar
R. Mohammand and F. Ahmad, “Analytical solutions for the stresses and deformations of deep tunnels in an elastic-brittle-plastic rock mass considering the damaged zone,” Tunnelling and Underground Space Technology, vol. 58, pp. 186–196, 2016.
View at: Publisher Site | Google Scholar
B. Behnam, S. Fazlollah, and M. Hamid, “Prediction of plastic zone size around circular tunnels in non-hydrostatic stress field,” International Journal of Mining Science and Technology, vol. 24, no. 1, pp. 81–85, 2014.
View at: Publisher Site | Google Scholar
S. Deng, Y. Zheng, L. Feng, L. Tuan, C. Yue, and Z. Wu, “Analysis of plastic zones in surrounding rocks around a circular tunnel considering the effect of intermediate principal stress and heterogeneity,” in International Conference on Advances in Civil Engineering and Materials (ACEM2018) and 1st World Symposium on Sustainable Bio-Composite Materials and Structures (SBMS1), pp. 1–6, Nanjing, China, 2018.
View at: Google Scholar
X. Guo, Z. Zhao, X. Gao, X. Wu, and N. Ma, “Analytical solutions for characteristic radii of circular roadway surrounding rock plastic zone and their application,” International Journal of Mining Science and Technology, vol. 29, no. 2, pp. 263–272, 2019.
View at: Publisher Site | Google Scholar
X. Guo, N. Ma, X. Zhao, Z. Zhao, and Y. Li, “General shapes and criterion for surrounding rock mass plastic zone of round roadway,” Journal of China Coal Society, vol. 41, no. 8, pp. 1871–1877, 2016.
View at: Google Scholar
H. Dong, Q. Gao, Z. Zhang, R. Chen, and N. Liu, “Approximate solution and numerical simulation for plastic zone of roadway surrounding rock under unequal compression in two directions,” Journal of China Coal Society, vol. 44, no. 11, pp. 3360–3368, 2018.
View at: Google Scholar
M. Xu, S. Wu, Y. Gao, J. Ma, and Q. Wu, “Analytical elastic stress solution and plastic zone estimation for a pressure- relief circular tunnel using complex variable methods,” Tunnelling and Underground Space Technology, vol. 84, pp. 381–398, 2019.
View at: Publisher Site | Google Scholar
Y. Zhang, H. Qiang, and Y. Yang, “Unified solutions to mixed mode crack tip under small scale yielding,” Chinese Journal of Mechanical Engineering, vol. 43, no. 2, pp. 50–54, 2007.
View at: Publisher Site | Google Scholar
X. Gao, X. Kang, H. Wang, and Y. Zhang, “Analysis of crack tip plastic zone for composites under small scale yielding,” Journal of Mechanical Engineering, vol. 45, no. 7, pp. 254–258, 2009.
View at: Publisher Site | Google Scholar
X. Gao, H. Wang, X. Kang, and L. Jiang, “Analytic solutions to crack tip plastic zone under various loading conditions,” European Journal of Mechanics - A/Solids, vol. 29, no. 4, pp. 738–745, 2010.
View at: Publisher Site | Google Scholar
C. Li, J. Xu, J. Pan, and C. Ma, “Plastic zone distribution laws and its types of surrounding rock in large-span roadway,” International Journal of Mining Science and Technology, vol. 22, no. 1, pp. 23–28, 2012.
View at: Publisher Site | Google Scholar
H. Shi and M. Bai, “Plastic zone analysis of deep-buried, noncircular tunnel and application on the high-speed railway in the karst area,” Mathematical Problems in Engineering, vol. 2017, Article ID 9523267, 10 pages, 2017.
View at: Publisher Site | Google Scholar
C. Guo, H. Han, H. Kong, and L. Shi, “Explicit form of exact analytical solution for calculating ground displacement and stress induced by shallow tunneling and its application,” Advances in Civil Engineering, vol. 2019, Article ID 5739123, 13 pages, 2019.
View at: Publisher Site | Google Scholar
Z. Zou, J. Yang, Z. Wang, and H. Liu, “An OCaNet model based on octave convolution and attention mechanism for iris recognition,” Mathematical Problems in Engineering, vol. 2021, Article ID 3412060, 10 pages, 2021.
View at: Publisher Site | Google Scholar
H. Shi, Z. Ma, Q. Zhu, J. Shi, and Z. Zhao, “Comparison of shape characteristics of plastic zone around circular tunnel under different strength criteria,” Journal of Mechanics, vol. 36, no. 6, pp. 849–856, 2020.
View at: Publisher Site | Google Scholar
P. Zhou and P. Wu, “Study on strength theory effect of plastic zone distribution of roadway surrounding rock,” Mathematical Problems in Engineering, vol. 2022, Article ID 6165652, 9 pages, 2022.
View at: Publisher Site | Google Scholar
Y. Ma, A. Lu, X. Zeng, and H. Cai, “Analytical solution for determining the plastic zones around twin circular tunnels excavated at great depth,” International Journal of Rock Mechanics and Mining Sciences, vol. 136, p. 104475, 2020.
View at: Publisher Site | Google Scholar
Y. Ma, A. Lu, H. Cai, and X. Zeng, “Analytical solution for determining the plastic zones around two unequal circular tunnels,” Tunneling and Underground Space Technology, vol. 120, p. 104267, 2020.
View at: Publisher Site | Google Scholar
Z. Hu, X. An, F. Li, and Y. Zhang, “The shape characteristics of circular tunnel surrounding rock plastic zone in the complex stress field,” Arabian Journal of Geosciences, vol. 15, no. 2, 2022.
View at: Publisher Site | Google Scholar
Z. Wang, Complex Variable Function Solution of Deep Buried Tunnel Excavation, Shijiazhuang Tiedao University, Shijiazhuang, 2018.

Copyright

Copyright © 2022 Shuhong Liu 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.

PDF Download Citation

Download other formats

Order printed copies

Views

196

Downloads

233

Citations