1China Railway 16st Bureau Group Road and Bridge Engineering Co Ltd, Beijing, China
Vibroengineering PROCEDIA, Vol. 42, 2022, p. 125-131.
Received 30 March 2022; received in revised form 14 April 2022; accepted 26 April 2022; published 16 May 2022
56th International Conference on Vibroengineering, Zittau, Germany, May 16, 2022
This work provides new findings on the theory of D–P failure criterion applied to the study of the analytic solution of the radius and stress of tunnel plastic zone and the periphery displacement of hydraulic pressure tunnel under uniform ground stress field. It focuses on elasto-plastic rock masses with the stage of construction and operation. The results that are obtained are compared to other calculation techniques and show that: The elastic-plastic solution calculated based on the theoretical model is larger than that based on the D-P Criterion, indicating safer conditions. The principal stress varies with the internal water pressure, and it firstly decreases to zero and then gradually increase in the plastic zone. With the internal water pressure keeping increasing to a high level, the principal stress in the elastic zone changes from tangential stress to radial stress in the radial distance-enlarging direction, which tends to reach the geostress. In this condition, the radial distance is related to the value of internal water pressure. The calculation results of this study are demonstrated to be in good agreement with those of numerical simulation under the same conditions.
Keywords: hydraulic tunnel, D-P criterion, elastic-plastic solution, the first principal stress.
With the construction of large-scale water conservancy projects such as the Three Gorges Hydropower Station, Jinping Hydropower Station, and the South-to-North Water Diversion Project, a large number of hydraulic tunnels are excavated in China. In these processes, several factors significantly influence the stability and functional use of the tunnels, including the unloading of surrounding rocks, the redistribution of stress, and the emergence of plastic zones. As a result, the stress and the elastic-plastic deformation of hydraulic pressure tunnel have been receiving attention from both the academic and engineering fields.
The studies of axisymmetrical tunnels in perfect elasto-plastic materials and with the Mohr-Coulomb type failure criterion was elaborated by some researchers. When the internal water pressure is extremely high in the hydraulic tunnel, the principal stress changes, producing new plastic zones. In this context, Fenner formula and Kastner formula will be inapplicable. However, there are only a few research achievements about the high-pressure hydraulic tunnel during the operation stage. In the previous theoretical analyses of the elasto-plastic behavior of the rock mass around a circular tunnel [1-8], the circular tunnel is subjected to a hydrostatic in situ stress at infinity, the plane strain assumption are taken into account. For a preliminary design of circular axisymmetric plane strain openings encountered in tunneling problem, an exact and close-form solution with Mohr-Coulomb (M-C) criterion is normally feasible. For the nonlinear yield criteria, the Hoek-Brown criterion [9-11] is also accepted in rock mechanics community, as it provides a reliable tool for predicting the strength of jointed rock mass. However, the D-P failure criterion is widely used in the numerical calculation for the characteristic of no edges and corners. For such cases the D-P failure criterion has been found to be more suitable.
In this paper, the elasto-plastic analysis of rock around a circular pressure tunnel subjected to axial in situ stress and internal water pressure are discussed, using the D-P failure criterion. Effects of the axial in situ stress on the plastic zones, elasto-plastic stresses and periphery displacement around a circular tunnel is analyzed. Influence of the parameters , and on the elasto-plastic stresses, and plastic zones around a circular pressure tunnel is also investigated.
2. Elastic-plastic solution of circular tunnels
The circular tunnel is assumed in the infinite elastic body, with the tunnel excavation radius of , the surrounding rock elastic parameters of and , strength parameters of and , the country rock stress of , the internal water pressure of , and the plastic zone radius of (Fig. 1).
The principal stress of the surrounding rock varies with the working conditions. During the construction stage characterized by low internal water pressure, the principal stress is radial stress, while during the operation stage characterized by high internal water pressure, the principal stress is tangential stress. The radial stress as the first principal stress.
Fig. 1. Calculation model: 1 – surrounding rock; 2 – country rock
2.1. Stress solutions
The excavation of tunnels is usually associated with the release of geostress, the unloading of surrounding rocks, the redistribution of stress, and the emergence of plastic zones. In axisymmetric case regardless of the body stress, the axisymmetric equilibrium equation under the plane strain condition is:
D-P yield criterion can be expressed as:
We substitute that and rearrange the balance equation, getting Eq. (3):
where is the integral constant, determined by the boundary conditions:
When we substitute the and , we can get:
Eq. (5) is the expression of the plastic zone radius in the construction stage when the first principal stress is radial. The plastic zone radius changes with geostress, internal water pressure (or supporting pressure), strength parameters and excavation radius.
2.2. Deformation calculation
In the condition that the deformation is small, considering the planar deformation and the material incompressibility:
where is the integration constant. Therefore, we can get:
Considering the general Hooke law, can be calculated according to Eq. (8), we can get:
Eq. (8) is the displacement expression.
2.3. The tangential stress as the first principle stress
2.3.1. Stress solutions
In the tunnel operation stage, the tunnel wall can withstand certain internal water pressure. Assuming the internal water pressure is unevenly distributed and disregarding the water weight, the operation condition is still axisymmetric, which means the balance Eq. (1) still apply. If we substitute into the yield equation, we can get:
When , stresses are continuous, we obtain:
Eq. (10) expresses the plastic zone radius during the construction stage when the first principal stress is tangential. The radius varies with relevant parameters including geostress, internal water pressure (support pressure), intensity, and excavation radius.
2.3.2. Deformation calculation
can be calculated according to Eq. (10), we can get:
Eq. (12) is the displacement expression.
To sum up, during the construction stage, the first principal stress is radial stress, and the plastic zone radius should be calculated according to Eq. (5); during the operation stage, the first principal stress is tangential stress, and the plastic zone radius should be modified according to Eq. (12).
3. Model validation and example analysis
3.1. Model validation
In order to verify the correctness of the model, two different methods respectively based on the D-P criteria and M-C criteria are implemented to calculate the plastic zone radius and stress, with the model parameter 3 m, normal stress 15 MPa, high geostress 30 MPa, internal friction angle 25°, and cohesion force 6 MPa. The plastic zone radius based on the D-P criteria and M-C criteria under the normal stress ( 15 MPa, 0) are 3.575 m and 3.495 m, respectively, while those under the high geostress ( 30 MPa, 0) are 4.870 m and 4.690 m, respectively.
The variation of elastic zone radius with strength parameters is shown in Figs. 2 and 3. Generally, the two methods agree with each other well. When the cohesion force equals or is smaller than 1.2 MPa, the elastic zone radius quickly increases with the decrease of , which can provide significance guidance for the engineering construction. When the cohesion force is larger than 1.2 MPa, the influence of on is small. Under the high geostress ( 30 MPa), friction angle dramatically influences the plastic zone radius (Fig. 3).
Fig. 2. The variation of with
Fig. 3. The variation of with
The stress distribution of surrounding rocks under different geostress are shown in Figs. 4 and 5. The tangential stress reaches the maximum in the junction part of plastic and elastic zones. Moreover, the radial stress increases with the geostress, changing from 0 to the geostress gradually. It is identified by both the two methods that tangential and radial stresses are close to the country rock stress at the location the radius is .
Fig. 4. Stress distribution under 15 MPa
Fig. 5. Stress distribution under 30 MPa
3.2. Example analysis
According to the critical condition , the plastic zone radius during the construction stage, previously expressed as Eq. (5), can be rewritten as:
If we substitute the known strength parameters into Eq. (13), we can obtain that 3.5585 MPa. It means when 3.5585 MPa, there emerges plastic zones in the surrounding rocks, and the plastic zone radius varies with the water pressure. When 3.5585 MPa, there are no plastic zones in the sounding rocks. Similarly, the critical value calculated based on the M-C criterion is that 3.2229 MPa, slightly smaller than the result of this study.
During the operation stage, the water pressure is extremely high, contributing to the yield of surrounding rocks again. The plastic zone radius increases with the water pressure, which is expressed as Eq. (14). Given the critical condition that , we can obtain:
If we substitute the known strength parameters into Eq. (14), we can obtain that 26.441 MPa. It means when 26.441 MPa, there emerges plastic zones again in the surrounding rocks, and the plastic zone radius varies with the water pressure. When 26.441 MPa, there are no plastic zones in the sounding rocks (shown in Fig. 6).
Fig. 6. The relationship between the plastic zone radius and internal water pressure
In conclusion, the variation of stresses of surrounding rocks with internal water pressure can be divided into three stages. The first stage (I) is the low water pressure stage ( 3.5585 MPa), also referred as the construction stage. The second stage (II) is the medium water pressure stage (3.5585 MPa 26.441 MPa), during which there exists no plastic zone, and the stress distribution can be calculated with the elastic Rummy solutions. The final stage (II) is the high-water pressure stage ( 26.441 MPa), during which the plastic zone area increases with the internal water pressure.
Under the normal stress, the stress distribution during the high-water pressure stage is further analyzed. Figs. 7 and 8 show the stress distribution with the internal water pressure of 30 MPa and 40 MPa, respectively, with corresponding plastic radius of 3.4875 m and 4.9911 m. Moreover, when the tangential or the radial stress is 15 MPa (geostress), the critical radius corresponding to the two different internal water pressure are 6.2863 m and 9.0120 m, respectively. When the radius is larger than the critical value, the first principal stress is tangential again.
Fig. 7. Stress distribution with 30 MPa
Fig. 8. Stress distribution with 40 MPa
Based on the D-P yield criterion, the elastic-plastic analytical solutions of pressure circular tunnels are derived under the plane strain condition considering different working conditions. Finite element models are established under the same condition, which help reach the following conclusions.
1) Elastic-plastic analytical solutions of pressure circular tunnels during the construction stage are verified to be accurate by comparison with results based on M-C criterion. When 1.2 MPa, the plastic zone radius dramatically increases with the decrease of cohesion force. However, when 1.2 MPa, the plastic zone radius gradually decreases with the increase of cohesion force. Under the high geostress ( 30 MPa), friction angle significantly influences the development of the plastic zone. Generally, the plastic zone radius by the model is larger than that based on the M-C criterion, which indicates safer condition.
2) The plastic zone radius decreases with the water pressure when 3.5585 MPa. Then the elastic zone disappears when 3.5585 MPa26.441 MPa. Again, there emerges the plastic zone when 26.441 MPa, with the radius increasing with the water pressure. During the high-water pressure stage, the principle stress of the elastic rock body changes from tangential to radial stress when , and approaches to the geostress. Meanwhile, the radial distance is related to the internal water pressure.
- Y. Wang, “Ground response of circular tunnel in poorly consolidated rock,” Journal of Geotechnical Engineering, Vol. 122, No. 9, pp. 703–708, Sep. 1996, https://doi.org/10.1061/(asce)0733-9410(1996)122:9(703) [Publisher]
- H. Wu, Q. Fang, and Z.-K. Guo, “Zonal disintegration phenomenon in rock mass surrounding deep tunnels,” Journal of China University of Mining and Technology, Vol. 18, No. 2, pp. 187–193, Jun. 2008, https://doi.org/10.1016/s1006-1266(08)60040-1 [Publisher]
- G. Hou and X. Niu, “Perfect elastoplastic solution of axisymmetric circular openings in rock mass based on Levy-Mises constitutive relation and D-P yield criterion,” (in Chinese), Rock and Soil Mechanics, Vol. 30, No. 6, pp. 1555–1562, 2009, https://doi.org/10.16285/j.rsm.2009.06.003 [Search CrossRef]
- Detournay E. and Fairhurst C., “Two-dimensional elastoplastic analysis of a long, cylindrical cavity under non-hydrostatic loading,” International Journal of Rock Mechanics and Mining Sciences and Geomechanics Abstracts, Vol. 25, No. 1, pp. A46–A47, Feb. 1988, https://doi.org/10.1016/0148-9062(88)93094-x [Publisher]
- X. Chen and Cp Tan, “Haberfield CM. Solutions for the deformations and stability of elastoplastic hollow cylinderssubjected to boundary pressures,” International Journal for Numerical and Analytical Methods in Geomechanics, Vol. 23, No. 8, pp. 779–800, 1999, https://doi.org/10.1002/(sici)1096-9853(199907)23:8 [Publisher]
- A. I. Sofianos, “Tunnelling Mohr-Coulomb strength parameters for rock masses satisfying the generalized Hoek-Brown criterion,” International Journal of Rock Mechanics and Mining Sciences, Vol. 40, No. 3, pp. 435–440, Apr. 2003, https://doi.org/10.1016/s1365-1609(03)00017-0 [Publisher]
- A.-Z. Lu, G.-S. Xu, F. Sun, and W.-Q. 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, Jan. 2010, https://doi.org/10.1016/j.ijrmms.2009.07.003 [Publisher]
- J. F. Labuz and A. Zang, “Mohr-Coulomb Failure Criterion,” Rock Mechanics and Rock Engineering, Vol. 45, No. 6, pp. 975–979, Nov. 2012, https://doi.org/10.1007/s00603-012-0281-7 [Publisher]
- X. P. Zhou and J. L. Li, “Hoek-Brown criterion applied to circular tunnel using elastoplasticity and in situ axial stress,” Theoretical and Applied Fracture Mechanics, Vol. 56, No. 2, pp. 95–103, Oct. 2011, https://doi.org/10.1016/j.tafmec.2011.10.005 [Publisher]
- S. K. Sharan, “Elastic-brittle-plastic analysis of circular openings in Hoek-Brown media,” International Journal of Rock Mechanics and Mining Sciences, Vol. 40, No. 6, pp. 817–824, Sep. 2003, https://doi.org/10.1016/s1365-1609(03)00040-6 [Publisher]
- S. K. Sharan, “Exact and approximate solutions for displacements around circular openings in elastic-brittle-plastic Hoek-Brown rock,” International Journal of Rock Mechanics and Mining Sciences, Vol. 42, No. 4, pp. 542–549, Jun. 2005, https://doi.org/10.1016/j.ijrmms.2005.03.019 [Publisher]