A method to evaluate the stability of tunnel face is proposed in the framework of upper bound theorem. The safety factor which is widely applied in slope stability analysis is introduced to estimate the stability of tunnel face using the upper bound theorem of limit analysis in conjunction with a strength reduction technique. Considering almost all geomaterials following a nonlinear failure criterion, a generalized tangential technique is used to calculate the external work and internal energy dissipation in the kinematically admissible velocity field. The upper bound solution of safety factor is obtained by optimization calculation. To evaluate the validity of the method proposed in this paper, the safety factor is compared with those calculated by limit equilibrium method. The comparison shows the solutions derived from these two methods match each other well, which shows the method proposed in this paper can be considered as effective.
1. Introduction
Shallow tunnels are very common nowadays in municipal engineering construction as they make traveling more convenient and also reduce engineering costs. Though shallow tunnels have a number of advantages, accidents due to instability of tunnels face would occur if the supporting pressure cannot resist the earth pressure. Thus, to avoid the occurrence of such accidents, it is necessary to determine the suitable retaining pressure on the tunnel face. The stability analysis of the face for a tunnel excavated in shallow strata has drawn the attention of many investigators [1–4]. On the other hand, the retaining pressure on the tunnel face is constant in practical engineering when the earth pressure balance of tunnel face is achieved by using compressed air or bentonite slurry. Therefore, it is necessary to establish a method to evaluate the stability of tunnel face when the retaining pressure is known. The safety factor (FOS) is a valid index to evaluate the stability of geotechnical structure which is widely used in geotechnical engineering. Presently, the safety factor used in geotechnical stability analysis is mostly calculated by the finite element method with shear strength reduction technique (SSRFEM) [5–9]. Based on this method, the safety factor and critical failure surface can be obtained by reducing the actual shear strength parameters until geotechnical material fails. So the accuracy of FOS derived from SSRFEM relies strongly on the determination of failure of geotechnical structure. However, the definition of failure for geotechnical structure is a disputed issue. Griffiths and Lane [10] chose the nonconvergence of algorithm as the indicator of failure for homogeneous slope. Snitbhan and Chen [11] regarded bulging or progressive loss of ground along the vertical cut as the failure of vertical slope. As there is divergence of definition for failure in the geotechnical engineering field, the different definition of failure would induce different analysis result, which confines the application of SSRFEM in geotechnical engineering.
Large numbers of theoretical researches and practical projects illustrate that the failure of shallow tunnel face can be divided into collapse and blow-out failure mechanisms. The collapse failure mechanism is used to describe the collapse of soil in front of the tunnel face, which is caused by the insufficient retaining pressure on the tunnel face. Contrarily, blow-out failure mechanism occurs when the retaining pressure is so great that soil is heaved in front of the tunnel face. Based on the failure mode and kinematically admissible velocity characteristic of the shallow tunnel face, Leca and Dormieux [12] proposed a three-dimensional failure mechanism composed of solid conical blocks. Due to the slide between the solid conical blocks and surrounding soil, the plastic flow occurs along the velocity discontinuity surface. Using the energy dissipation rate along the surface and the rate of work caused by external force, Leca and Dormieux [12] derived the upper bound solution of retaining pressure for tunnel face. As the mechanism proposed by Leca and Dormieux [12] is supported by centrifuge model tests and well reflects the mechanical characteristics of failure mode for shallow tunnel face, many scholars used this mechanism to analyze the stability of shallow tunnel face under various conditions [13–16].
These literatures mentioned above all used linear Mohr-Coulomb criterion. However, the stress-strain relation of soils and rocks is nonlinear. This viewpoint has been supported by experiments and some scholars used nonlinear failure criterion to study the stability of tunnel face and other geotechnical structures [17, 18]. Based on the failure mechanism proposed by Leca and Dormieux [12], Huang and Yang [19] calculated the upper bound solution of retaining pressure on the tunnel face using the upper bound theorem in conjunction with nonlinear failure criterion. Senent et al. [20] studied the face stability of circular tunnels excavated in heavily fractured rock masses which are subjected to the nonlinear Hoek-Brown failure criterion. According to their study, the critical retaining pressures computed with limit analysis are very similar to those obtained with the numerical model, which proves their method is valid. Though their studies reflect the influence of nonlinearity on the critical retaining pressures, they failed to propose a method to evaluate the stability of tunnel face when the retaining pressure is known.
In this work, upper bound theorem combined with shear strength reduction technique is used to estimate the face stability of a tunnel excavated in shallow strata following the nonlinear failure criterion. On the basis of upper bound theorem, the rate of external work and the internal energy rate of dissipation for the failure mechanism are calculated. Based on the relationship between the rate of external work and the internal energy rate of dissipation, the convergence of iteration in the strength reduction calculation can be controlled, which avoids the selection of the definition of failure that occurs in the SSRFEM calculation process. To validate the new methodology, the FOS of tunnel face is compared with the result computed by limit equilibrium method. Furthermore, the influence of nonlinear parameter on the stability of tunnel face is investigated.
2. Upper Bound Theorem with Shear Strength Reduction Technique
Shear strength reduction technique was proposed by Bishop [22], whose core content is the reduction of soil shear strength parameters until the soil fails. To achieve this reduction, an important concept is introduced, the shear strength reduction factor Fs. When the actual shear strength parameters c and ϕ are divided by the shear strength reduction factor, the soil strength parameters cf and ϕf used in upper bound analysis are obtained(1)cf=cFs,ϕf=arctantanϕFs.While shear strength reduction factor Fs increases incrementally, newly reduced soil strength parameters are obtained. The iterative process continues until failure occurs.
As mentioned above, the selection of a suitable definition of failure is a problem in SSRFEM. To overcome this difficulty, the upper bound theorem is used to control the convergence of iterative operation in the strength reduction technique. The upper bound theorem states that when the velocity boundary condition is satisfied, the load derived by equating the rate of external work to the rate of the energy dissipation in any kinematically admissible velocity field is no less than the actual collapse load. Therefore, if the reduced soil strength parameters are introduced in the energy dissipation calculation, the shear strength reduction factor can be obtained on the basis of the relationship between the rate of external work and the rate of the energy dissipation.
3. Upper Bound Solution of FOS for Shallow Tunnel Face
In this work, the failure mechanisms proposed by Leca and Dormieux [12] are used to calculate the rate of external work and the rate of the energy dissipation in the framework of the upper bound theorem of limit analysis. The failure mechanism of Leca and Dormieux [12] is composed of collapse and blow-out failure mechanisms, which is shown in Figure 1. For instance, we use failure mechanism III to illustrate the computational process of the upper bound solution of FOS. As Leca and Dormieux [12] have computed the rate of external work Pe for mechanism III, the expression of Pe can be written as(2)Pe=πD24-RBRC2cosϕσs-13sinαRBRC/sinα3-RA/cosα3sin2ϕcosϕγDV,where D is the tunnel diameter, γ is the unit weight of the soil, ϕ is the friction angle of the soil, σs is surcharge loading, V is the velocity of conical block, α is the angle between symmetry axis of conical block and centre line of the tunnel, and the parameters RA, RB, RC are expressed as(3)RA=cosαcosα+ϕcosα-ϕ,RB=sinαsinα+ϕsinα-ϕ,RC=sin2α+2C/D+1sin2ϕcos2ϕ-cos2α,where C is the tunnel depth. Moreover, the rate of the energy dissipation produced in kinematically admissible velocity field is(4)PV=πD24RBRC2-RAcVsinϕ,where c is the cohesion of soil. Based on the upper bound theorem, the expression of retaining pressure σT is obtained by equaling the rate of external work to the rate of the energy dissipation(5)σT=PV-PeπD2RA/4cosϕV.To calculate the upper bound solution of FOS, the initial cohesion c and friction angle ϕ of soil are substituted into (1) to obtain the reduced strength parameters cf and ϕf. Then, the reduced strength parameters cf and ϕf are substituted into (2) and (4) to derive the upper bound solution of retaining pressure σT. Finally, by equaling the practical retaining pressure to the retaining pressure σT expressed in (5), the objective function of safety factor fα,Fs which includes an angle variable is obtained. However, the objective function fα,Fs is just an expression of numerous upper bound solutions. According to the upper bound theorem, the minimum value of objective function fα,Fs is the real upper bound solution of FOS. Therefore, a sequential quadratic programming is employed to search the minimum value of objective function fα,Fs. As some compatibility relations of velocity should be satisfied in the kinematically admissible velocity field, the optimization calculation is achieved when corresponding constraint conditions are satisfied. Therefore, the expression of mathematical planning for the problem can be written as (6)minFs=Fsα,Fs,s.t.Pγ+PT=PV,ϕ<α<π2-ϕ,where PT is the power of practical retaining pressure.
Failure mechanism of tunnel face proposed by Leca and Dormieux [12].
Failure mechanism II
Failure mechanism III
4. Comparison with the Results Computed by Limit Equilibrium Method
To evaluate the validity of the method proposed in this work, the FOS of tunnel face is calculated by limit equilibrium method and the upper bound theorem with shear strength reduction technique. Numerical results for these two methods with different parameters are presented and compared with each other.
Based on the silo theory, a three-dimensional failure mode of the tunnel face composed of wedges, as shown in Figure 2, is proposed by Anagnostou and Kovári [21] in the framework of limit equilibrium method. According to Bishop [22], the FOS for slope can be defined as the ratio of the available shear strength of the soil to that required to maintain equilibrium. Therefore, the FOS of tunnel face derived from the wedge model presented in Figure 2 is expressed as follows:(7)FOS=FrFi,where Fr is the shearing resistance force and Fi is the shearing force acting on the wedge base. The detailed calculating procedure of Fr and Fi can be seen in the appendix.
Wedge stability model of tunnel face proposed by Anagnostou and Kovári [21].
Dawson et al. [23] pointed out that if the actual height of slope is equal to the critical height computed by the upper bound theorem, the FOS of the slope is exactly 1.0. Similarly, if the actual retaining pressure applied on the tunnel face is equal to the upper solution of retaining pressure, the FOS of tunnel face calculated by this method should also be 1.0.
Table 1 presents FOS of tunnel face calculated by LEM and LAM for C/D = 0.4 and γ = 20 kN/m^{3}, with the ϕ varying from 20° to 45°. The retaining pressures σT presented in Table 1 are the upper solutions computed by Leca and Dormieux [12]. It is found from Table 1 that the FOS of tunnel face computed by the method proposed in this paper is almost equal to those using the limit equilibrium method. Therefore, the good agreement between these two methods shows that the method proposed in this paper is an effective method for calculating the FOS of shallow tunnel face.
Comparison of FOS for tunnel face.
Number
C/D
γ (kN/m^{3})
ϕ (°)
σT (Pa)
FOS
LEM
LAM
1
0.4
20
20
24960
1.0442354
1.0085697
2
0.4
20
25
18120
1.0395553
1.0189546
3
0.4
20
30
12960
1.0307114
1.0094596
4
0.4
20
35
10560
1.031479
1.0074599
5
0.4
20
40
7920
1.0314455
1.0098653
6
0.4
20
45
6000
1.0327247
1.0112345
5. Nonlinear Failure Criterion and Generalized Tangential Technique
Zhang and Chen [24] adopted a nonlinear expression to describe the relationship between the normal and shear stress when a plastic flow of geotechnical materials occurs, which is expressed as (8)τ=c01+σnσt1/m,where σn and τ are the normal and shear stresses on the failure surface, respectively, and c0, σt, and m are material parameters determined by geotechnical test. As the strength envelope of this nonlinear failure criterion is curve, the strength parameter of geotechnical materials cannot be determined like linear Mohr-Coulomb failure criterion. To overcome this difficulty, Yang and Yin [25] proposed a generalized tangential technique which uses a tangential line on the nonlinear failure criterion at point M to determine the strength parameter. The tangential line on the curve at the location of tangency point M as shown in Figure 3 can be written as(9)τ=ct+σntanϕt,where ϕt is a tangential frictional angle and ct is the intercept of the straight line on the τ-axis. ct and ϕt can be obtained by the following two expressions:(10)ct=m-1m·c0m·σt·tanϕtc01/1-m+σt·tanϕt,tanϕt=c0mσt1+σnσt1-m/m.Since the tangential line in Figure 3 is random, ϕt is regarded as a variable to calculate the rate of external work and energy dissipation. Using sequential quadratic programming, the upper bound solution of objective function and corresponding value of ϕt are obtained. Obviously, the tangential line is determined by optimization calculation, which indicates that nonlinear failure criterion represented by the tangential line will provide the optimum upper bound of actual load for the geotechnical material.
Generalized tangential technique for a nonlinear failure criterion.
6. Upper Bound Solution of FOS with Nonlinear Failure Criterion
To study the effect of the nonlinear parameter on the FOS of tunnel face, upper bound solutions of FOS are calculated by the method proposed in this paper with different parameters. As nonlinear coefficient m which controls the shape of strength envelope for nonlinear failure criterion has significant influence on the other parameters, this paper focuses on the effect of m on FOS. Using the collapse failure mechanism, the upper bound solutions of FOS with different nonlinear coefficient m are plotted in Figure 4 with σT = 5 kPa, γ = 20 kN/m^{3}, c0 = 10 kPa, σt = 24.73 kPa, and depth ratio C/D = 0.2~1. It can be found from Figure 4 that the FOS decreases with the increase of m and C/D for collapse failure mechanism. Similarly, Figure 5 shows the change law of FOS for varying m and C/D with parameters corresponding to σT = 800 kPa, γ = 20 kN/m^{3}, c0 = 90 kPa, and σt = 247.3 kPa on the basis of blow-out failure mechanism. Figure 5 shows that the values of FOS tend to increase directly with the depth ratio C/D and inversely with m.
Effect of m and C/D on the FOS for collapse failure mechanism (σT = 5 kPa, γ = 20 kN/m^{3}, c0 = 10 kPa, σt = 24.73 kPa, C/D = 0.2~1, m = 1.0~2.0).
Effect of m and C/D on the FOS for blow-out failure mechanism (σT = 800 kPa, γ = 20 kN/m^{3}, c0 = 90 kPa, σt = 247.3 kPa, C/D = 1~3, m = 1.0~2.0).
It can be seen that the change trend of depth ratio C/D on the FOS for collapse failure mechanism is just opposite to that of the blow-out failure mechanism. This phenomenon is caused by the different velocity of solid conical blocks in kinematically admissible velocity field for these two failure mechanisms. In the kinematically admissible velocity field of collapse failure mechanism, the velocity component of solid conical block in front of the tunnel face moves in the opposite direction of the tunnel excavation. On the contrary, in kinematically admissible velocity field of blow-out failure mechanism, induced by the huge retaining pressure on tunnel face, the velocity component of solid conical block in front of the tunnel face moves in the direction of the tunnel excavation. As the velocities of solid conical block in these two failure mechanisms are opposite, the effects of power of the soil weight on the energy dissipation calculation are also opposite. Therefore, when collapse failure occurs on the tunnel face, low value of depth ratio C/D will contribute to the tunnel stability. However, when blow-out failure occurs on the tunnel face, high value of depth ratio C/D will contribute to the tunnel stability. On the other hand, with the increase of nonlinear coefficient m the values of ct and ϕt obtained by generalized tangential technique both decrease. As ct and ϕt are used to calculate the energy dissipation and the energy dissipation is independent of velocity of solid conical blocks, the change laws of nonlinear coefficient m on FOS for these two failure mechanisms are the same.
7. Conclusion
The upper bound theorem combined with shear strength reduction technique is adopted to calculate the FOS of shallow tunnel face in the framework of nonlinear failure criterion. Upper bound solutions of FOS are derived from collapse and blow-out failure mechanisms proposed by Leca and Dormieux [12]. Using generalized tangential technique, the nonlinear failure criterion is introduced in the energy dissipation calculation.
The upper bound solutions of FOS are compared with the results calculated by limit equilibrium method. The solutions of FOS derived in this paper are almost equal to those calculated by limit equilibrium method, which proves that the method proposed in this paper is effective.
Based on the motion features of solid conical block in kinematically admissible velocity field for these two failure mechanisms, the reason for different effect of depth ratio on FOS is explained. By parameter study, it is found that the parameters m and C/D have large influence on the FOS of shallow tunnel face.
Appendix
According to Qiao et al. [26], the detailed calculating procedure of Fr and Fi can be illustrated as follows. On the basis of Terzaghi relation soil pressure theory, the three-dimensional relation soil pressure is (A.1)σv=γ-λ·cλ·k0·tanϕ·1-e-λ·k0·Z·tanϕ+P0·e-λ·k0·Z·tanϕ,where k0 is lateral pressure coefficient, γ is the unit weight of the soil, c is the cohesion of soil, ϕ is the friction angle of the soil, Z is the distance between tunnel roof and ground surface, P0 is surcharge pressure, and λ can be calculated by the following expression:(A.2)λ=2tanα-1/tanβ+1/sinβB1-1/tanα,where α, β, and B are parameters which can be seen in Figure 2. The soil pressure is (A.3)Pv=σv·ACNEF,where ACNEF is the area of CNEF. The weight of the wedge is(A.4)G=D22·cotα·B·γ,where D is tunnel diameter. Using static equilibrium equation, the pressure N is obtained(A.5)N=Psinα+G+Pvcosα,where P is retaining pressure. Based on the Mohr-Coulomb failure criterion, the friction resistance T of slip surface AMEF in the wedge is derived(A.6)T=c·AAMEF+N·tanϕ,where AAMEF is the area of slip surface AMEF in the wedge. Thus, the shearing resistance force Fr and the shearing force Fi can be written as(A.7)Fr=T+P·cosα,Fi=G+Pv·sinα.
Competing Interests
The authors declare that the mentioned received funding in the Acknowledgments did not lead to any competing interests regarding the publication of this paper.
Acknowledgments
The preparation of this paper received financial support from the National Natural Science Foundation of China (nos. 51308072 and 51278071), Educational Commission of Hunan Province of China (no. 15C0052), and Innovation Research Project of Priority Key Disciplines at Changsha University of Science and Technology (no. 15ZDXK13). The financial supports are greatly appreciated.
AhmedM.IskanderM.Evaluation of tunnel face stability by transparent soil modelsChenR. P.LiJ.KongL. G.TangL. J.Experimental study on face instability of shield tunnel in sandKimS. H.TononF.Face stability and required support pressure for TBM driven tunnels with ideal face membrane—Drained caseOresteP.Face stabilization of deep tunnels using longitudinal fibreglass dowelsZhengH.SunG. H.LiuD. F.A practical procedure for searching critical slip surfaces of slopes based on the strength reduction techniqueWangT.WuH. G.LiY.GuiH.ZhouY.ChenM.XiaoX.ZhouW.ZhaoX.Stability analysis of the slope around flood discharge tunnel under inner water exosmosis at Yangqu hydropower stationWeiW. B.ChengY. M.Strength reduction analysis for slope reinforced with one row of pilesHuangM. S.JiaC. Q.Strength reduction FEM in stability analysis of soil slopes subjected to transient unsaturated seepageMarandiS. M.AnvarM.BahramiM.Uncertainty analysis of safety factor of embankment built on stone column improved soft soil using fuzzy logic α-cut techniqueGriffithsD. V.LaneP. A.Slope stability analysis by finite elementsSnitbhanN.ChenW.-F.Elastic-plastic large deformation analysis of soil slopesLecaE.DormieuxL.Upper and lower bound solutions for the face stability of shallow circular tunnels in frictional materialLeeI. M.LeeJ. S.NamS. W.Effect of seepage force on tunnel face stability reinforced with multi-step pipe groutingLiY.EmeriaultF.KastnerR.ZhangZ. X.Stability analysis of large slurry shield-driven tunnel in soft clayLeeI. M.NamS. W.The study of seepage forces acting on the tunnel lining and tunnel face in shallow tunnelsWangH. T.JiaJ. Q.Face stability analysis of tunnel with pipe roof reinforcement based on limit analysisYangX. L.LongZ. X.Seismic and static 3D stability of two-stage rock slope based on Hoek-Brown failure criterionYangX. L.XuJ. S.LiY. X.Collapse mechanism of tunnel roof considering joined influences of nonlinearity and non-associated flow ruleHuangF.YangX. L.Upper bound solutions for the face stability of shallow circular tunnels subjected to nonlinear failure criterionProceedings of the GeoShangai International Conference, Deep and Underground Excavations2010Shangai, China251256SenentS.MollonG.JimenezR.Tunnel face stability in heavily fractured rock masses that follow the Hoek–Brown failure criterionAnagnostouG.KováriK.Face stability conditions with Earth-Pressure-Balanced shieldsBishopA. W.The use of the slip circle in the stability analysis of earth slopesDawsonE. M.RothW. H.DrescherA.Slope stability analysis by strength reductionZhangX. J.ChenW. F.Stability analysis of slopes with general nonlinear failure criterionYangX. L.YinJ. H.Slope stability analysis with nonlinear failure criterionQiaoJ. L.ZhangY. T.GaoJ.Application of strength reduction method to stability analysis of shield tunnel face