#### 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).