Elastoplastic Analysis of Circular Tunnel Based on Drucker–Prager Criterion
Based on the Drucker–Prager yield criterion, the theoretical solution of stratigraphic deformation in tunnel excavation process is deduced by the cavity expansion theory. In view of soil loosening around the tunnel caused by the tunnel excavation process, the internal friction angle of the surrounding soil is not a constant but a function of normal stress. The piecewise linearization of the nonlinear yield function is used to analyze the elastoplastic solution of the cylindrical hole shrinkage. A comparison is conducted with a plastic zone in which the internal friction angle of the soil remains unchanged. It can be concluded that the radial stress, the tangential stress, the radial strain, and the tangential strain around the inner wall calculated from the former are smaller.
After subway tunnel excavation, the stress of surrounding soil will be redistributed. When the stress on the surrounding soil exceeds the elastic limit of rock mass and enters into the plastic state, elastoplastic analysis of the surrounding soil must be conducted to determine the stability of the subway tunnel and to provide a basis for the quantitative design of roadway support [1,2]. For deep tunnels, to analyze the surrounding soil stress and displacement generated by tunnel excavation and support, it is not necessary to consider the surface boundary conditions and the surrounding soil can be regarded as an infinite domain and can be solved .
There are many examples of prediction of ground displacement caused by tunnel excavation based on the theory of cavity expansion. Mair and Taylor  used the Tresca criterion to predict the deformation of soil around the tunnel in clay and compared them with those displacements measured in Green Park and Regent Park in London. Yu and Rowe  proposed theoretical and semitheoretical solutions to the expansion of the borehole and the expansion of the sphere in the unloading process during drained and undrained conditions. Pinto and Whittle  predicted and expound the deformation of shallowly buried tunnel in the soft soil layer, based on the cavity expansion theory. Mo et al.  deduced the theoretical solution of the concentric regions of two different rock soils and also carried out the finite element simulation. For a long time, the elastoplastic analysis of the tunnel has been conducted by adopting the Mohr–Coulomb criterion [4, 8–10] or Hoek–Brown criterion [11–13], based on the fact that the internal friction angle of the surrounding soil is a constant. However, the internal friction angle is a function of normal stress rather than a constant. Kennedy and Lindberg  derived small strain solutions based on the segmented Mohr–Coulomb yield surface. Florence and Schwer  extended the linear Mohr–Coulomb criterion solution. Recently, Zou et al. [16, 17] focused on employing the elastoplastic theories of circular tunnel to solve geotechnical engineering problems and achieved outstanding results.
In this paper, the Drucker–Prager criterion is used as the yield condition, and the plastic analysis of deep-buried circular subway tunnels is performed by using the piecewise linearization of nonlinear yield function. A comparison is conducted with a plastic zone in which the internal friction angle of the soil remains unchanged.
2. Drucker–Prager Yield Criterion
The Mohr–Coulomb failure criterion can be adopted to describe the yield and failure mechanisms of geomaterials. However, singularity will occur when the directional derivative is on the yield surface at the apex or ridgeline of the hexagonal cone [18, 19]. In order to overcome the shortcoming of the Mohr–Coulomb failure criterion, Drucker and Prager originally proposed the lower limit of the inscribed Mohr–Coulomb yield surface in 1952 , which can be expressed aswhere is the first invariant of the Cauchy stress and is the second invariant of the deviator stress tensor. In the formula, the material constants and can be written aswhere and are the internal friction angle and the cohesion of soil, respectively.
Assuming the compressive stress as positive and the tensile stress as negative in the calculation, we can obtain the following relation on the basis of the cylindrical cavity expansion theory:
Thus, the first invariant of the Cauchy stress and the second invariant of the deviator stress tensor can be simplified as
Based on the positional relationship between the Drucker–Prager criterion and the Mohr–Coulomb criterion on the plane (Figure 1), a variety of expressions for the parameters and can be obtained. The parameters and for various yield criteria are listed in Table 1.
3. Governing Equation
Before analyzing the elastic and plastic zones of the tunnel, the governing equations should be derived firstly. According to the theory of elasticity, the equilibrium equation of the infinitesimal body around the tunnel can be expressed aswith boundary conditions
In addition, strain-displacement relations of the soil can be given by
After the strain components and displacement are obtained, the stress-strain relation involved in plane strain problem can be given bywhere
4. Multisection Drucker–Prager Criterion Solution
In this paper, the elastoplastic closed-form solution of cylindrical hole contraction in the unrestricted Drucker–Prager medium is analyzed by piecewise linearization of nonlinear yield function. The pore radius is set as a, and the hydrostatic pressure is uniformly applied to isotropic soil. Tunnel wall pressure is reduced very slowly, so the power effect can be ignored. As the internal pressure decreases, the inner wall of the hole will fail first. As the pressure of the inner wall of the hole decreases further, the plastic zone will expand outward until it reaches a certain radius . Eventually, the material in the annular area is in the plastic state, and the other area is purely elastic.
4.1. Stress Distribution in the Elastic Zone
As the tunnel pressure p decreases from p0, the deformation of the soil around the tunnel is purely elastic at the beginning. According to the theory of elasticity, the elastic solutions of stress components can be expressed as
With the continuous reduction of the pore pressure, the inner wall will trigger an initial yielding when pore pressure reaches the following condition:
4.2. Stress Distribution in the Plastic Zone
When the pore wall reaches the initial yielding, the plastic area will be formed near the inner wall in the region with decrease of the pressure . As mentioned above, the plastic area is divided into m annular zone, as shown in Figure 2. For different plastic zones, there are different yield conditions represented by different stress functions; hence, each zone can be analyzed separately.
For the plastic zone 1, the yield criterion is given by
By solving the Euler equation, the general solution for stress components in the plastic zone 1 can be obtained:
Additionally, the stress components in the outer elastic zone can be expressed as
The integral constants and can be solved by using the boundary and continuous conditions. Finally, we can obtain
Similarly, the equilibrium equation for the ith plastic zone can be written aswhere
The general solution can be obtained as follows by solving Equation (27).
According to the continuous condition
The stress expressions for other plastic zones can be derived. For the zone , stress components can be obtained as
The boundary radius of each region can be determined by the radial stress continuous conditions at the junction of the zone and the zone , which yield
4.3. Displacement Analysis
This paper mainly discusses the strain and displacement in the elastic and plastic zones. It is convenient to work out the strain expression in the elastic zone by solving the simultaneous Equations (12), (13), (15), and (16). However, the law of plastic flow is required for the calculation strain and displacement in the plastic zone, which will be analyzed by following the associated flow rule.
Firstly, we calculate the strain and displacement distributions of the plastic zone 1 in the region , which is close to the outer elastic area. The relationship between two plastic strain components can be obtained by referring to the associated flow rule:
Equation (33) can be alternatively written aswhere the tangential plastic strain is undetermined. When Equations (34) and (35) are substituted into Equation (11), the governing equation for can be obtained as followswhere
The general solution can be derived as follows by solving Equation (36):where the integral constant is determined by the elastoplastic boundary condition:
Thus, the integral constant can be solved:
Thus, the strain components and displacement in zone 1 are expressed as
Similarly, strain and displacement expressions in other plastic zones can be obtained. The solution of the integral constants can be determined by displacement continuous conditions.
5. Case Study
There is a tunnel with the radius , the cohesion of the surrounding soil C = 5 MPa, the internal friction angle of the surrounding soil , and the initial hydrostatic pressure =28.2 MPa acting uniformly on the soil. After tunnel excavation, the pressure of the inner wall of the tunnel will decrease slowly. When the pressure on the inner wall of the tunnel decreases to 18 MPa, 12.63 MPa, and 3.17 MPa, respectively, the distributions of radial stress and tangential stress (dimensionless) around the tunnel are shown in Figure 3.
It can be seen from Figure 3 that the tunnel is in an elastic state when the pore pressure is greater than 12.63 MPa, and the radial stress increases with the increase of . In contrast, the tangential stress has opposite tendency. Once the inner wall pressure p decreases to 12.63 MPa, plasticity commences around the inner wall of the tunnel, so part of surrounding soil enters the plastic state. The tangential stress in the plastic area increases with the increase of , whereas the tangential stress in the elastic zone decreases with the increase of . In addition, both the radial and tangential stresses far from the tunnel are almost equal to the hydrostatic pressure.
When the tunnel wall pressure p decreases to 3.17 MPa, the radius of the plastic area varies with the internal friction angle and cohesion C of the soil, as shown in Figures 4 and 5, respectively. It can be noted from Figures 4 and 5 that the radius of the plastic area decreases with the increase of the inner friction angle and cohesion C of the surrounding soil. It indicates that it is essential to use varying inner friction angle and cohesion C of the surrounding soil to characterize the mechanical behavior of tunnel excavation.
Notably, it is an efficient way to calculate the elastoplastic solution of tunnel excavation by using piecewise linear yield function to represent the nonlinear yield criterion. Based on the influence of the normal stress on the internal friction angle of the soil, we divided the plastic area into seven zones. The yield criterion for each zone with various internal friction angles can be expressed as
Comparison of stress components between one plastic area with internal friction angle and seven separated plastic zones with various internal friction angles is shown in Figure 6, and the corresponding strain components’ comparison is shown in Figure 7.
In practical engineering, the closer the distance from the tunnel is, the looser the soil becomes during tunnel excavation and the smaller the internal friction angle becomes. Figures 6 and 7 depict that the radius of the plastic area increases when the plastic area is divided into seven zones. The stress and strain components around the inner wall of the tunnel decrease. However, the values of stress and strain components in the elastic area are very close.
To have a better understanding of the stress and strain distributions around the tunnel, finite element software ANSYS is used for numerical simulation . Two-dimensional element PLANE42 is used to simulate surrounding soil. Under symmetrical conditions, a-quarter mode with 1,600 elements and 1,683 nodes is constructed, as depicted in Figure 8. The radius of the surrounding soil is 100 m, and the radius a of the tunnel is equal to 3 m. Comparison between finite element results and analytical results is illustrated in Figures 9 and 10, respectively.
According to Figures 9 and 10, we consider tunnel excavation as the two-dimensional plane strain problem in the theoretical derivation, without considering the influence of strain along the tunnels. This leads to a larger analytical solution of tangential stress in plastic area and a smaller radius in plastic area. Notably, by using the associated flow rule to analyze the plastic strain, the analytical solution of near the inner wall of the tunnel becomes larger. The good agreement in elastic area between finite element results and theoretical results validates the model developed in the current study. Furthermore, our model results are in good agreement with the numerical results proposed by Zhang et al. .
Taking into account loosening of the soil around the tunnel caused by the tunnel excavation, the internal friction angle of rock and soil is not a constant but a function of normal stress. The piecewise linearization of nonlinear yield function was used to analyze the elastoplastic solution of shrinkage of cylindrical pores based on the Drucker–Prager yield criterion. A comparison is conducted with a plastic zone in which the internal friction angle of the soil remains unchanged. Based on this study, the following conclusions can be drawn:(1)The tunnel excavation process leads to decrease of the pressure on the inner wall of the tunnel and increase of the radius of the plastic area. The radius of the plastic area is related to the internal friction angle and cohesion of the soil in the tunnel.(2)The radius of the plastic area decreases with increasing of the inner friction angle and cohesion of the surrounding soil.(3)Compared with only one plastic zone (internal friction angle ), when the plastic area is divided into seven zones, the radius of the plastic area increases and the stress and strain components around the inner wall of tunnel decrease, but the distributions of stress and strain in the elastic area are very close.
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.
The authors gratefully acknowledge the financial support provided by the National Natural Science Foundation of China (nos. 11032005 and 11402099), the Science and Technology Scheme of Guangdong Province (no. 2012A030200003), the Science and Technology Scheme of Guangzhou City (no. 1563000451), and the Fundamental Research Funds for the Central Universities of China (no. 17817004).
Y. Xiang, General Theory of Tunnel Mechanics, Science Press, Beijing, China, 2014, in Chinese.
R. J. Mair and R. N. Taylor, “Prediction of clay behavior around tunnels using plasticity solutions. Predictive Soil Mechanics,” in Proceedings of the Wroth Memorial Symposium, pp. 27–29, Oxford, UK, July 1992.View at: Google Scholar
Z. Qiang, S. Wang, and G. E. Xiurun, “Elastoplastic analysis of circular openings in strain-softening rock masses,” Chinese Journal of Rock Mechanics & Engineering, vol. 29, pp. 1031–1035, 2010.View at: Google Scholar
E. Hoek and E. T. Brown, Underground Excavations in Rock, Institute of Mining and Metallurgy, London, UK, 1980.
T. C. Kennedy and H. E. Lindberg, “Tunnel closure for nonlinear Mohr-Coulomb functions,” Journal of the Engineering Mechanics Division, vol. 104, pp. 1313–1326, 1978.View at: Google Scholar
A. S. Khan and S. Huang, Continuum Theory of Plasticity, John Wiley & Sons, Hoboken, NJ, USA, 1995.
K. M. Thompson and M. J. Thompson, ANSYS Mechanical APDL for Finite Element Analysis, Butterworth-Heinemann, Woburn, MA, USA, 2017.
X. Zhang, G. Zhao, and X. Meng, “Elastoplastic analysis of surrounding rock on circular roadway based on Drucker-Prager yield criterion,” Journal of China Coal Society, vol. 38, pp. 30–37, 2013.View at: Google Scholar