Affiliation(s):
1.
Key Laboratory of Urban Security and Disaster Engineering, Ministry of Education, Beijing University of Technology, Beijing 100124, China; moreAffiliation(s): 1.
Key Laboratory of Urban Security and Disaster Engineering, Ministry of Education, Beijing University of Technology, Beijing 100124, China; 2.
Beijing Collaborative Innovation Center for Metropolitan Transportation, Beijing 100124, China; 3.
Key Laboratory for Urban Underground Engineering of the Education Ministry, Beijing Jiaotong University, Beijing 100044, China; less
Peng-fei Li, Qian Fang, Ding-li Zhang. Analytical solutions of stresses and displacements for deep circular tunnels with liners in saturated ground[J]. Journal of Zhejiang University Science A, 2014, 15(6): 395-404.
@article{title="Analytical solutions of stresses and displacements for deep circular tunnels with liners in saturated ground", author="Peng-fei Li, Qian Fang, Ding-li Zhang", journal="Journal of Zhejiang University Science A", volume="15", number="6", pages="395-404", year="2014", publisher="Zhejiang University Press & Springer", doi="10.1631/jzus.A1400023" }
%0 Journal Article %T Analytical solutions of stresses and displacements for deep circular tunnels with liners in saturated ground %A Peng-fei Li %A Qian Fang %A Ding-li Zhang %J Journal of Zhejiang University SCIENCE A %V 15 %N 6 %P 395-404 %@ 1673-565X %D 2014 %I Zhejiang University Press & Springer %DOI 10.1631/jzus.A1400023
TY - JOUR T1 - Analytical solutions of stresses and displacements for deep circular tunnels with liners in saturated ground A1 - Peng-fei Li A1 - Qian Fang A1 - Ding-li Zhang J0 - Journal of Zhejiang University Science A VL - 15 IS - 6 SP - 395 EP - 404 %@ 1673-565X Y1 - 2014 PB - Zhejiang University Press & Springer ER - DOI - 10.1631/jzus.A1400023
Abstract: To investigate the displacement and stress distributions for deep circular tunnels with liners in saturated ground, an analytical model is proposed. For a deep tunnel with drainage conditions, plane strain conditions at any cross-section of the tunnel and the elastic regime of the linear elasticity for the remaining liner are assumed, while the ground is assumed to be linearly elastic and perfectly plastic with a failure surface defined by the Mohr-Coulomb criterion. The post-yield behavior of the ground follows the non-associated flow rule defined by the dilation angle. To solve the proposed problem, two procedures are presented. An axisymmetric model for a deep circular tunnel with a steady-state seepage condition is considered, and then a simple closed-form analytical solution is obtained with a common theoretical framework for the boundary conditions of a constant total head along the tunnel circumference. Assuming that certain ground displacements along the tunnel circumference have occurred before the installation of the liner, analytical solutions of stresses and displacements are derived with particular emphasis on the seepage and the stress release effect induced by tunnelling. The proposed analytical model is validated by numerical simulation.
Darkslateblue:Affiliate; Royal Blue:Author; Turquoise:Article
Article Content
1. Introduction
One of the important issues during the construction of deep tunnels in saturated ground is the appropriate design of the tunnel liners. The convergence-confinement method (CCM) has become very popular for tunnel construction throughout the world due to its technical feasibility, safety, and economic competitiveness (Carranza-Torres and and Fairhurst, 2000; Oreste, 2003a; 2003b; González-Nicieza et al., 2008). According to the concept of the CCM, the tunnel liners should be designed to sustain the loads transmitted from the surrounding ground during excavation as well as the seepage forces due to the water flow towards the tunnel, if any (Lee and Nam, 2001; Bobet, 2003; Nam and Bobet, 2006). The effects of seepage force and in-situ stress on displacements and stress distributions for deep tunnels in saturated ground have been investigated by numerous researchers (Schweiger et al., 1991; Lee and Nam, 2001; Bobet, 2001; 2003; Arjnoi et al., 2009; Li, 2011).
Although there are a lot of numerical procedures capable of estimating the ground displacements and liner stresses, analytical solutions are still desirable in direct applications or in the validation and verification of numerical models for their advantages of conciseness and convenience (Huangfu et al., 2010).
Analytical solutions of the distribution of the pore pressure in the rock mass for steady seepage of deep tunnles are also obtained in previous studies. For instance, based on the theory of complex potentials, the theory of Cauchy integrals, and of singular integral equations, the steady-state fluid flow solution for fractured porous media was obtained by Liolios and Exadaktylos (2006). They also provided a numerical method for the accurate estimation of pore pressure and pore pressure gradient fields due to specified hydraulic pressure or pore pressure gradient acting on the lips of one or multiple non-intersecting curvilinear cracks in a homogeneous and isotropic porous medium, but there is a fly in the ointment for which they could not provide the complete analytical solutions of displacements and stresses for tunnels in saturated ground.
Complete analytical solutions of displacements and stresses for tunnels in saturated ground are obtained by Bobet (2001; 2003) and Nam (2006). These solutions cover different construction processes and ground conditions, including dry or saturated ground, shallow or deep, with or without air pressure, and with or without a gap between the ground and the liner. However, these solutions are restricted to cases where ground displacements are small, since the ground and the liner are assumed to behave elastically.
Based on the assumptions that circular tunnels are excavated in elasto-plastic homogeneous ground under hydrostatic stress conditions, a series of analytical or semi-analytical analyses have been performed (Fenner, 1938; Brown et al., 1983; Carranza-Torres and Fairhurst, 2000; Sharan, 2005; Lee et al., 2000; Fang et al., 2013; Han et al., 2013). All these studies have contributed to the understanding of the surrounding ground-liner interaction mechanism, but the results do not apply to the tunnels in saturated ground.
Considering the effect of seepage forces, ground stresses are obtained by Li et al. (2004) and Liu et al. (2009) for deep circular tunnels in elastic-perfectly plastic material with a Mohr-Coulomb type yield criterion while neglecting the existence of liners. Their studies were complemented by Lu and Xu (2009), who proposed an analytical model to solve the stresses for a subsea tunnel with liner. Note that all these studies are not capable of providing the solutions of ground displacements in the plastic zone.
In this study, elasto-plastic plane strain solutions of the stress and displacement distributions around deep circular lined tunnels in isotropic saturated ground due to uniform ground loads and seepage forces are presented. The derivation process for the expression of displacements in the plastic zone is a highlight in this paper. A gap between the ground and the liner, which represents a certain ground displacement along the tunnel circumference occurred before the installation of the liner is also considered. The proposed model is particularly suitable for the investigation of deep circular tunnels excavated in soft saturated ground.
2. Problem description
The fundamental objective of this paper is to study the stresses and displacements of tunnels with liners in saturated ground considering the effect of both the seepage force and in-situ stress. To derive the analytical solution for the proposed problem, the following assumptions are made: (1) The circular tunnel is located in a fully saturated, homogeneous and isotropic aquifer subjected to uniform in-situ stress. (2) Plane strain conditions are applicable at any tunnel cross-section. (3) The liner behaves in a linear elastic manner. (4) The ground behaves in an elastic-perfectly plastic manner, obeys the linear Mohr-Coulomb yield criterion, and is considered weightless. (5) The water is incompressible, and Darcy’s law is applicable. (6) Before the installation of the liner, certain ground displacements along the tunnel circumference have already occurred. Therefore, there is a constant gap between the ground and the liner along the perimeter of the tunnel (Bobet, 2001), which can represent the displacements due to any physical gap, the unsupported tunnel face, the liner construction lag, or the time effect of shotcrete hardening process. For example, the installation delay of the support will inevitably cause certain ground displacements. Moreover, even though the shotcrete lining is installed immediately after the excavation, certain ground displacements will gradually take place as the shotcrete hardens (Oreste, 2003b).
A schematic model of the proposed problem is shown in Fig. 1. In this model, ri is the inner radius of the liner, rl is the outer radius of the liner, and r0 is the outer boundary of the problem. In this study, r0=20ri. The outer boundary is far enough from the inner tunnel boundary, therefore the presence of the outer boundaries have limited impact on the solutions. This method has been adopted by numerous researchers (Bobet, 2003; Li et al., 2004; Nam and Bobet, 2006; Liu et al., 2009). There are no displacement constraints at the boundaries (Nam and Bobet, 2006). Ambient stresses σi and σ0 are the radial pressures (the tension is positive) at the inner boundary and the outer boundary, respectively. hi and h0 are the water head on the inner boundary and the outer boundary, respectively. kl is the hydraulic conductivity of the liner. kp and ke are the hydraulic conductivities of the plastic zone and the elastic zone, respectively. El and μl are the elasticity modulus and the Poisson’s ratio of the liner, respectively. Ep and μp are the elasticity modulus and the Poisson’s ratio of the ground in the plastic zone, respectively. Ee and μe are the elasticity modulus and the Poisson’s ratio of the ground in the elastic zone, respectively. cp, φp, and ψp are the cohesion, the internal friction angle, and the dilation angle of the ground, respectively.
Fig.1 Schematic illustration of the analytical model
The yield zone radius, rp, the water head on the outer of the liner, hl, the water head on the plastic zone, hp, the radial stress at the ground-liner interface, σl, and the radial stress at the plastic-elastic interface, σp can be derived by using the seepage and mechanical theory.
The radial stress, σr, the tangential stress, σθ, the water head, h, the displacement, u, are the functions of the radial distance from the tunnel center, r, and will be given analytically by a theoretical model in this study.
3. Proposed analytical model
3.1. Analytical solutions of the seepage field
It follows from the symmetry with respect to the axis of the tunnel that the flow is radial, and the discharge vectors at all points are directed toward the center of the tunnel. According to Bear (1972), following the continuity equation, there is , then , where C1 and C2 are two integration constants.
The corresponding boundary conditions for the problem are
Let Ql, Qp, and Qe be the water discharge through the ground in the elastic zone, in the plastic zone, and the liner, respectively. According to Wang et al. (2008), when ri≤r≤rl, there is , analogously,
For Ql=Qe=Q, we can obtain: , , .
3.2. Analytical solutions of stresses and displacements
3.2.1. Basic equations and boundary conditions
The seepage force is a body force (Lambe and Whitman, 1969), of which the direction is radial due to the symmetry of the flow. The seepage force per unit volume is , where γw is the specific weight of the water, ξ is the equivalent surface porosity of the fractured mass. Terzaghi (1936) and Leliavski (1947) determined from field and laboratory studies that ξ, the surface porosity, is generally close to unity and always greater than 0.85.
The solution of a deep tunnel in saturated ground must satisfy the equilibrium equation as follows: , where the subscripts r and θ denote the radial and tangential directions, respectively. Eq. (10) can be solved with the following boundary conditions: , , , , , , where ,, and are the radial stress, the tangential stress, and the radial displacement of the liner, respectively; , , are the radial stress, the tangential stress, and the radial displacement of the ground in the plastic zone, respectively; , , and are the radial stress, the tangential stress, and the radial displacement of the ground in the elastic zone, respectively. w is the gap between the ground and the liner.
3.2.2. Stresses and displacements of liner
The analytical solutions of stresses and displacements acting on the liner due to the ground loads and seepage forces are obtained by a theoretical model based on the solutions of Li (2004): , , , , , , , , , , , .
3.2.3. Stresses and displacements of the ground in the plastic zone
The excavation produces unloading in the radial direction and loading in the tangential direction according to the proposed problem. Therefore, the major principal stress is the tangential stress and the minor principal stress is the radial stress in the plastic zone. At failure, the 2D Mohr-Coulomb criterion can be expressed as f(σr, σθ)=0, i.e., .
According to the equilibrium equation (Eq. (10)), the boundary conditions (Eqs. (11c) and (11d)), and the failure criterion (Eq. (24)), the ground stresses in the plastic zone are:
, , where B is a constant and can be expressed as follows: .
Under the axisymmetric plain strain condition, the strains and the displacements are given by
The compatibility condition is expressed as .
Assuming a small strain deformation analysis, the total strains are decomposed into elastic and plastic components, as follows: where the superscripts e and p denote elastic and plastic strains, respectively.
Hooke’s law is applied to determine the radial and tangential strains in the elastic region under the plane strain condition: where Ee and μe denote the elastic modulus and Poisson’s ratio, respectively.
The elastic strains in the plastic zone can be readily obtained by substituting the stresses in the plastic zone into Hooke’s law (Eq. (31)). The plastic strains in the plastic zone are instead governed by an appropriate flow rule postulated for the yielding behavior. The flow rule of plasticity relating the plastic strain increment dεp to the plastic potential g is governed by ,
where dλ is a non-negative scalar function present throughout the plastic loading process. With the Mohr-Coulomb model, it is often assumed that the plastic potential takes the same form as the yield function while the friction angle is replaced by the dilation angle (which is commonly a smaller angle). The plastic potential can therefore be written as ,
where ψp is the dilation angle.
By substituting the plastic potential (Eq. (33)) into the flow rule (Eq. (32)), radial and tangential plastic strain increments can be obtained: where Nψ=(1+sinψp)/(1−sinψp). Integration for (Eq. (34)), there is
By substituting Eqs. (35) and (31) into Eq. (30), then into Eq. (29), a differential equation can be obtained:
.
The solution is ,
where Kx is a coefficient which can be determined by the plastic strain at the plastic-elastic interface, , as follows: .
By substituting Eq. (38) into Eq. (37), the plastic strain can be written as
The tangential strain in the plastic zone can be derived as follows:
By substituting Eq. (40) into Eq. (28), the solution of displacements Up can be written as
Note that for the proposed problem, the initial radial displacement of ground in the plastic zone u0p induced by the initial stress σ0 should be subtracted from the total displacement in order to calculate the net displacements induced by tunnelling alone. The expression of u0p is given by .
The expression for displacements due to tunnelling in the plastic zone can be derived as follows:
3.2.4. Stresses and displacements of ground in the elastic zone
Note that the initial radial displacement of ground in the elastic zone u0e induced by the initial stress σ0 should be subtracted from the total displacement to calculate the net displacements induced by tunnelling alone (Fang et al., 2013). The expression of u0e is given by .
The ground stresses and displacements in the elastic zone can be derived by using the same methods as calculating the distributions of liner stresses and displacements. Therefore, the solutions can be written as ,
,
,
where K'1–K'7 have the same expressions as K1–K7 but replacing ri, rl, hi, hl, μl, σi, and σl by rp, r0, hp, h0, μe, σp, and σ0, respectively.
3.2.5. Expressions of σri, σrp, and rp
The solutions in this study must satisfy the equilibrium equations, the strain compatibility equations, and the boundary conditions (Eqs. (11c)–(11f)).
According to Eqs. (11c), (11d), and (11f), we can obtain ,
,
where ,
.
By substituting Eqs. (48) and (49) into Eq. (11e), a transcendental equation can be obtained to find the value of rp. By substituting the value of rp into Eq. (6) and Eq. (7), the values of hl and hp can be obtained. Then, (σrl, σθl, ul), (σrp, σθp, up), and (σre, σθe, ue) can be obtained, and so forth.
4. Validation of the analytical solutions
4.1. Parameters used in validation
To validate the analytical model, the results of the analytical model are compared with those obtained from numerical computations. A MATLAB program is used to carry out the computations for the analytical model. The numerical computations are carried out by using a 3D finite difference code FLAC3D (Itasca Consulting Group, Inc., 1997; 2002) which can address 3D continuum problems with coupled groundwater flow. This comparison is intended to validate the solutions. A group of typical parameter values are employed, w=0, 7.5, 15, 23, 30.5, 38, 46, 53.5, 61, 68.5, and 76.5 (unit: mm), varied to simulate different stress releases before liner construction, and the others are shown in Table 1. Only a quarter of the geometry needs to be modeled due to symmetry. The geometry of the mesh and the boundary conditions are shown in Fig. 2. The normal velocities of the gridpoints along the bottom, vertical lateral, front and back boundary planes are fixed to be zero (Fang et al., 2013).
Table 1
Parameters used in the validation
Parameter
value
ri (m)
5
rl (m)
5.3
r0 (m)
100
hi (m)
0
h0 (m)
80
σi (MPa)
0
σ0 (MPa)
2
El (GPa)
10
Ep (GPa)
1
Ee (GPa)
1
ξ
1
μl
0.2
μp
0.35
μe
0.35
φp
40
ψp
0
kl (m/s)
1×10−8
kp (m/s)
1×10−6
ke (m/s)
1×10−6
cp (kPa)
100
Fig.2 Calculation model (a) and boundary conditions (b)
4.2. Relationship between rp−rl and w
The results of the calculations for the relationship between rp−rl and w are given in Fig. 3. It is clear from the plot that the solutions of the study and FLAC3D match well. All illustrative comparisons of the solution and the numerical results show little effect on the FLAC3D boundary constraints along the edges of the mesh. The yield zone radius increases significantly with the increase of w. This shows that the earlier the liner is carried out, the more significant the decreases in the plastic zone.
Fig.3 Relationship between rp−rl and w
4.3. Results of the calculation
Results of the calculation for the stresses and displacements are given in Figs. 4 and 5. The solutions of the study and FLAC3D are matched well.
Fig.4 Relationship between σrl/σ0 and w
Fig.5 Stress distributions around the tunnel when w=15 mm
From the relationship between σrl/σ0 and w in Fig. 4, it is clear that contrary to rp, the effective radial stress at the ground-liner interface, σrl decreases significantly with the increase of w, and σrl≈0 when w=76.5 mm. This shows that the earlier the liner is carried out, the more significant will be the increase in the liner stress. Thus, appropriate time for the liner construction should be selected for a safe and economical design.
Fig. 5 shows stresses distributions around the tunnel when w=15 mm. We can find that: (1) with the increase of distance from the tunnel centerline, the radial stress increases until it ceases to the in site stress, and yet the tangential stress in the ground increases to the maximum equaling two times that of the in site stress and then decreases until it ceases to the in site stress; (2) The tangential stress in the liner decreases with the increase of the distance from the tunnel centerline, and the discontinuity appears at the ground-liner interface because of different elastic moduli between the ground and liner.
5. Conclusions
The analytical solutions presented in this paper expand and complete previous studies, in that we have now obtained complete analytical solutions for a deep tunnel in saturated ground. Although the analytical solutions are valid only for cases under the given assumptions, it appears to be useful for a preliminary design of deep tunnels under high water levels to predict the water pressure, stresses and displacement distributions around the tunnel. Stresses and displacement distributions using the proposed analytical solutions are compared to the numerical solution of FLAC3D, which appears to match well. For more complicated scenarios, a numerical model would be flexible. Nevertheless, the theoretical expressions derived in this study can be used in validating the numerical models.
[1] Arjnoi, P., Jeong, J.H., Kim, C.Y., 2009. Effect of drainage conditions on porewater pressure distributions and lining stressed in drained tunnels. Tunnelling and Underground Space Technology, 24(4):376-389.
[2] Bear, J., 1972. Dynamics of Fluids in Porous Media. American Elsevier Publishing Company,New York :
[3] Bobet, A., 2001. Analytical solutions for shallow tunnels in saturated ground. Journal of Engineering Mechanics, 127(12):1258-1266.
[4] Bobet, A., 2003. Effect of pore water pressure on tunnel support during static and seismic loading. Tunnelling and Underground Space Technology, 18(4):377-393.
[5] Brown, E.T., Bray, J.W., Ladanyi, B., 1983. Ground response curves for rock tunnels. Journal of Geotechnical Engineering, 109(1):15-39.
[6] Carranza-Torres, C., Fairhurst, C., 2000. Application of the convergence-confinement method of tunnel design to rock-masses that satisfy the Hoek-Brown failure criterion. Tunnelling and Underground Space Technology, 15(2):187-213.
[7] Fang, Q., Zhang, D.L., Zhou, P., 2013. Ground reaction curves for deep circular tunnels considering the effect of ground reinforcement. International Journal of Rock Mechanics & Mining Sciences, 60:401-412.
[8] Fenner, R., 1938. Undersuchungen zur Erkenntnis des Gebirgsdruckes. Gluckauf, (in German),74:681-695.
[9] Gonzlez-Nicieza, C., lvarez-Vigil, ., Menndez-Daz, A.E., Gonzlez-Palacio, A.C., 2008. Influence of the depth and shape of a tunnel in the application of the convergence-confinement method. Tunnelling and Underground Space Technology, 23(1):25-37.
[10] Han, J.X., Li, S.C., Li, S.C., 2013. A procedure of strain-softening model for elasto-plastic analysis of a circular opening considering elasto-plastic coupling. Tunnelling and Underground Space Technology, 37:128-134.
[11] Huangfu, M., Wang, M.S., Tan, Z.S., 2010. Analytical solutions for steady seepage into an underwater circular tunnel. Tunnelling and Underground Space Technology, 25(4):391-396.
[12] Itasca Consulting Group, Inc., 1997. Fast Language Analysis of Continua Three Dimensions, :
[13] Itasca Consulting Group, Inc., 2002. FLAC-Fast Lagrangian Analysis of Continua, Version 4.0, Minneapolis :
[15] Lee, I.M., Nam, S.W., 2001. The study of seepage forces acting on the tunnel liner and tunnel face in shallow tunnels. Tunnelling and Underground Space Technology, 16(1):31-40.
[16] Lee, I.M., Kim, J.H., Reddi, L.N., 2000. Permeability reduction model of soil-geotextile system induced by clogging. Journal of the Korean Geotechnical Society, 16(4):107-116.
[17] Li, P.F., 2011. Study on Stability Analysis and Control of Subsea Tunnel Surrounding Rock, PhD Thesis, Beijing Jiaotong University, Beijing, China :
[18] Li, Z.L., Ren, Q.W., Wang, Y.H., 2004. Elasto-plastic analytical solution of deep-buried circle tunnel considering fluid flow field. Chinese Journal of Rock Mechanics and Engineering, 23:1291-1295.
[19] Liolios, P.A., Exadaktylos, G.E., 2006. A solution of steady-state fluid flow in multiply fractures isotropic porous media. International Journal of Solids and Structures, 43(13):3960-3982.
[20] Liu, C.X., Yang, L.D., Li, P., 2009. Elastic-plastic analytical solution of deep buried circle tunnel considering stress redistribution. Engineering Mechanics, 26:16-20.
[21] Lu, X.C., Xu, J.Y., 2009. Elastic-plastic solution for subsea circular tunnel under the influence of seepage field. Engineering Mechanics, 26:216-221.
[22] Nam, S.W., Bobet, A., 2006. Liner stresses in deep tunnels below the water table. Tunnelling and Underground Space Technology, 21(6):626-635.
[23] Oreste, P.P., 2003. Analysis of structural interaction in tunnels using the convergence-confinement approach. Tunnelling and Underground Space Technology, 18(4):347-363.
[24] Oreste, P.P., 2003. A procedure for determining the reaction curve of shotcrete lining considering transient conditions. Rock Mechanics and Rock Engineering, 36(3):209-236.
[25] Schweiger, H.F., Pottler, R.K., Steiner, H., 1991. Effect of seepage forces on the shotcrete liner of a large undersea cavern. Computer Methods and Advances in Geomechanics, :1503-1508.
[26] Sharan, S.K., 2005. Exact and approximate solutions for displacements around circular openings in elastic-brittle-plastic Hoek-Brown rock. International Journal of Rock Mechanics and Mining Sciences, 42(4):542-549.
[28] Wang, X.Y., Tan, Z.S., Wang, M.S., 2008. Theoretical and experimental study of external water pressure on tunnel liner in controlled drainage under high water level. Tunnelling and Underground Space Technology, 23(5):552-560.
Open peer comments: Debate/Discuss/Question/Opinion
Open peer comments: Debate/Discuss/Question/Opinion
<1>