# Semianalytical Solution and Parameters Sensitivity Analysis of Shallow Shield Tunneling-Induced Ground Settlement

Authors: Jifeng Liu; Huizhi Zhang
Journal: Advances in Materials Science and Engineering
 (2017)

---

## Abstract

The influence of boundary soil properties on tunneling-induced ground settlement is generally not considered in current analytic solutions, and the hypothesis of equal initial stress in vertical and horizontal makes the application of the above solutions limited. Based on the homogeneous half-plane hypothesis, by defining the boundary condition according to the ground loss pattern in shallow tunnel, and with the use of Mohr-Coulomb plastic yielding criteria and classic Lame and Kiersch elastic equations by separating the nonuniform stress field to uniform and single-direction stress field, a semiempirical solution for ground settlement induced by single shallow circular tunnel is presented and sensitivity to the ground parameters is analyzed. The methods of settlement control are offered by influence factors analysis of semiempirical solution. A case study in Beijing Metro tunnel shows that the semiempirical solution agrees well with the in situ measured results.


---

## Body

## 1. Introduction

With the rapid development of China’s urban infrastructure construction, the shield tunneling is widely used as a safe and efficient tunnel construction method in urban subway, underground pipeline, and so on. At present, there are still a series of problems needing to be solved in the metro shield tunneling construction; one of the problems that should be urgently solved is the prediction and control of the shield tunneling construction induced ground surface settlement [1, 2]. The research methods of this problem include the observational method, numerical simulations, and analytical methods [3].In recent years, a large number of analytical solutions for axisymmetric circular tunnels in strain plane state have been presented [3–11], but those analytic methods are seldom considering the influence of ground soil properties on shield tunnel induced ground settlement, and some methods are based on the assumption of homogeneous semi-infinite plane; the calculated results are suitable for the deep buried tunnel surface settlement calculation but are not fully applicable for most shallow-buried tunnel, which limits the use of those methods. Therefore, it is an urgent problem to find suitable methods to calculate shallow shield tunnel surface settlement.


## 2. Semianalytical Analysis of Shallow Single Tunnel

As shown in Figure1, the assumption of shallow tunnel surrounding rock analytic analysis is as follows: (1) the rock is homogeneous and isotropic; (2) tunnel excavation is a plane strain axisymmetry problem; (3) the original stress affecting the surrounding rock is not equal in the vertical and the horizontal direction; the static earth pressure coefficient is 0.5. The other parameters are tunnel buried depth of H (m), average overburden soil weight γ (kN/m3), excavation radius a (m), plastic radius R (m), and face support reaction stress σ0 (kPa).Figure 1
Initial stress field of shallow tunnel.Equilibrium differential equation for axisymmetric problem is(1)dσrdr+σr-σθr=0.During the tunnel excavation, the soil deformation is elastic; then the expression of the stress-strain relationship of axisymmetric elasticity is(2)εr=1+μE1-μσr-μσθεθ=1+μE1-μσθ-μσr,where μ is the Poisson ratio of surrounding rock of tunnel and E is Young’s modulus of surrounding rock of tunnel.Two-way nonisobaric state could be divided into the external pressure asσh uniform stress field and the external pressure as (σv-σh) one-way stress field. For uniform stress field and one-way stress field, the Lame and Kiersch formula is, respectively, used to calculate the stress field; the elastic zone can be obtained by summing the above calculated uniform stress field and one-way stress field of the surrounding rock: (3)σr=-σ0R2r2-121+K0σv1-R2r2+121-K0σv1-R2r21-2R2r2cos⁡2θσθ=σ0R2r2-121+K0σv1+R2r2+121-K0σv1+3R4r4cos⁡2θ.By (2) and (3), and considering the boundary conditions and symmetry, the elastic zone displacement u can be obtained, and the tunnel excavation induced relative displacement can be obtained by subtracting the initial stress generated initial displacement from the calculated elastic zone displacement u, as shown in(4)u′=12GR2rσ0-121+K0σv-121-K0×σvR2r2-41-μcos⁡2θ,where G is shear modulus of soil.The Mohr-Coulomb theory strength criterion is adopted to judge whether the soil is in plastic state. The yield condition is(5)σθ-σr=2sin⁡φ1-sin⁡φσr+ccot⁡φ,where c is soil cohesion, kPa, which can be obtained from the shear test and φ is soil internal friction angle, °, which can be obtained from the shear test.The stress expression of the plastic zone can be obtained by (1) and (5), shown as(6)σr=σ0+ccot⁡φra2sin⁡φ/1-sin⁡φ-ccot⁡φσθ=1+sin⁡φ1-sin⁡φσ0+ccot⁡φra2sin⁡φ/1-sin⁡φ-ccot⁡φ,where c is soil cohesion, kPa; φ is soil internal friction angle, °; a is tunnel excavating radius, m; σa is face supporting pressure, kPa; and r is the radius of the calculating stress point, m.Along the boundary between the elastic and plastic zone, the stress boundary conditions should be met at the same time; that is, there isσθPlasticr=R=σθElasticr=R when the radius r=R; the plastic radius R is shown as follows:(7)R=a1-sin⁡φ1/21+K0σv+1-K0σvcos⁡2θ+ccot⁡φσ0+ccot⁡φ2sin⁡φ/1-sin⁡φ.Whena≤r≤R, the Mohr-Coulomb yield criterion nonassociative plastic flow method is used to solve the displacement [12]. The tunnel excavation is unloading process; (8) can be obtained by nonassociated flow rule.(8)εrεθ=-1+sin⁡ψ1-sin⁡ψ,where ψ is dilatancy angle of soil.Substituting (2) and (3) into (8), (9) can be available.(9)εr+1+sin⁡ψ1-sin⁡ψεθ=1m1-1+sin⁡ψ1-sin⁡ψ×μ1-μσr+1+sin⁡ψ1-sin⁡ψ-μ1-μσθ,where m=E/1-μ2.For the soil or rock, the value ofm is very large, the value to the right of (9) is very small, and then (9) can be simplified to (10) by assuming no deformation in the tunnel axis direction and neglecting the elastic strain in the plastic zone [6]. (10)εr+1+sin⁡ψ1-sin⁡ψεθ=0.From the geometric equationεr=du′/dr, εθ=u′/r, and (10), the relative displacement of the plastic zone can be solved as follows: (11)u′=12GR2/1-sin⁡ψr1+sin⁡ψ/1-sin⁡ψσ0-121+K0σv-12×1-K0σv4μ-3cos⁡2θsin⁡φ.The deformation of (11) is radial uniform, which is different from the actual situation. Park [13] presented that the ground deformation model for shallow-buried tunnel displacement boundary conditions is adopted to modify (11) so that it can better predict the subsidence of shallow-buried tunnels. The displacement boundary expression of Park is shown as (12), where u0 is the uniform radial deformation boundary, which is the value of (11) at the r=a position; the meaning of θ is shown in Figure 1, sin⁡θ=-(z-H)/r, cos⁡θ=x/r, and z is the vertical distance of the forecast point to the ground surface. The deformation pattern of (12) is shown in Figure 2.(12)urr=a=-u045+3sin⁡θ-3cos2⁡θ(13)u′=R2/1-sin⁡ψr1+sin⁡ψ/1-sin⁡ψσ0-121+K0σv-12×1-K0σv4μ-3cos⁡2θ×18Gsin⁡φ5+3sin⁡θ-3cos2⁡θ.Figure 2
Boundary condition of displacement.


## 3. Parameters Sensitivity Analysis of Semianalytical Solution

For a given depth and radius of the tunnel, taking the soil horizontal pressure coefficientK0=0.5, it is shown that the influence parameters of the ground deformation are K0, Young’s modulus E, Poisson’s ratio μ, soil cohesion c, internal friction angle φ, and dilatancy ψ. The calculated results of (13) by Matlab software are shown as Figures 3–7. Figure 3 shows the influence of Young’s modulus on the settlement when c=30 kPa, ψ=5°, φ=20°, and μ=0.2; Young’s modulus is inversely proportional to the maximum subsidence of the surface. Figure 4 shows the influence of Poisson’s ratio on the ground settlement, which shows that the maximum settlement increases with increase of Poisson’s ratio μ.Figure 3
Influence ofE to surface settlement.Figure 4
Influence ofμ on surface settlement.Figure 5
Influence ofφ on surface settlement.Figure 6
Influence ofc on surface settlement.Figure 7
Influence ofψ on surface settlement.Figure5 is the influence of soil internal friction angle φ on surface settlement when E=15.0 MPa, ψ=5°, c=30 kPa, and μ=0.2, Figure 6 is the influence of cohesion c on surface settlement when E=15.0 MPa, ψ=5°, μ=0.2, and φ=20°, and Figure 7 is the influence of dilatancy ψ on surface settlement when E=15.0 MPa, c=30 kPa, μ=0.2, and φ=20°. Figures 5–7 show that the maximum surface settlement decreases with E and c increase and increases with ψ, μ, and φ increase. In comparison, the influence of E, μ, ψ, and φ on the maximum surface settlement is larger than c.The above analysis for the influence of surrounding soil properties on the ground subsidence shows that there are two ways to control the shield tunnel induced ground settlement: one way is the active way, which is to reduce the surrounding soil properties weakening degree by optimizing the construction parameters of shield tunneling, such as choosing the rationale supporting face force and excavating speed or shield tail grouting reinforcement; the other way is the passive protection measures, which is to improve soil physical parameters by engineering methods to reduce the surface settlement of shield construction, the feasible methods including bolt reinforcement, partition wall, pile foundation underpinning, or the affected building body reinforcement [14].


## 4. Case Study

The 11th Project of Beijing Metro Line 10 is underground twin tunnels and is constructed by shield tunnel; the 6 m diameter tunnels are buried 16 m deep with nonsynchronous construction; the left tunnel is excavated more than 200 m ahead of the right tunnel. The ground elevation is 38.15 m. The shield tunnels are mainly through ⑤ the medium-coarse sand, ⑥ silty clay, ⑥1 clay, and ⑥2 silt. The perched groundwater is 2.74–3.58 m deep and the diving groundwater is 7.26–7.52 m deep. The layered and main mechanical parameters of foundation soil are found in Table 1.Table 1
Physical and mechanical parameters of foundation soils.
Number
Name
Weight/ kN/m3
Sr
e
I
p
Standard penetration/N
Fast shear test
UU
ModulusEs/MPa
Poisson’s ratioμ
Layer height/m
Cohesionc/kPa
Internal friction angle/°
Cohesionc/kPa
Internal friction angle/°
①
Miscellaneous fill
16.6
/
/
/
/
10
8
10
5
/
/
2.4
③1~③
Silt-silty clay
19.3
78.12–96.3
0.61–0.77
7.2–11.8
13–10
30.8
15.7
26.7
17.2
6.6–12.7
0.30
5.4
④
Silty clay
20.1
96.5
0.59
11.9
10
35
9
30
15
5.5–10
0.31
4.0
⑥
Silty clay
19.8
96.2
0.71
13.3
15
46
15
30
15
8.3–11.9
0.29
4.4
⑥
Silt
20.2
95.3
0.58
8.4
17
33
28
15
15
9.3–17.1
0.30
6.0
⑦1~⑦
Silty sand
19.9
54–48
0
32.8
30–40
0.28
3.6
⑦
Rounded gravel
21.1
0
40
60
0.25
1.4
⑧
Silty clay
19.7
96.4
0.73
12.9
25
54
5
9.2–13.6
0.30
5.0
⑨
Silty fine sand
20.6
52
0
30
35
0.27
1.6
⑩
Silty clay
19.7
20
20
8.7–13.1
0.29
6.2As shown in Figure8, the left shield tunnel machine is nearby underpassing the North Street 8# and South Street 8# residential building; in this underpassing process, the shield machine is just in the transition from circular curve to straight line, and the two tunnels center distance is 11.15 m~15.28 m, which puts great challenge to the control of strata settlement caused by the shield construction. In order to carry out the information construction and to reasonably control the ground surface subsidence for ensuring the safety of the affected building group, the corresponding settlement measurement point is laid out before the construction.Figure 8
Plan of tunnel line, buildings, and layout of surveying points.According to the relevant provisions of the national code of “code for design of building foundation (GB50007-2011),” taking into account the importance, the longtime service, and the poor deformation coordination of the upper panel type structure of the tunnels excavation affected buildings, the foundation deformation limit should be more strict than the national code; the previous studies suggest that the whole building tilt value should be controlled within 0.002, but the building tilt value of North Street 8# and South Street 8# residential building should be controlled within 0.001 and the maximum settlement under the wall which is nearby the street should be less than 10 mm because those buildings have tiled about 0.001 [15]. In order to test the impact of shield construction on the settlement of buildings, according to the mechanical properties of various strata and construction parameters, and considering the shield tunneling straightening on the surrounding soil disturbance, the following parameters are used for analytical analysis of the shield tunnel induced ground surface settlement, that is, the tunnel buried depth H=16 m, the static earth pressure coefficient K0=0.5, excavation radius a=3 m, weight of soil σv=310 kPa, face support reaction stress σ0=120 kPa, cohesion c=25 kPa, internal friction angle φ=20°, Young’s modulus E=15 MPa, Poisson ratio μ=0.24, and dilatancy angle ψ=3°. After the shield tunnel straightened, the shield tunneling causing surrounding soil properties weakening is relatively small, so the above-mentioned parameters were adjusted to c=30 kPa, φ=25°, E=20 MPa, μ=0.25, and ψ=5°. The two sets of parameters are taken into (13), and the Matlab7.0 is used for programming; the calculated results are shown in Figure 9. The maximum ground surface settlement is 23.8 mm when shield tunnel is in straightening phase, and the maximum ground surface settlement is 16.4 mm after posture adjustment.Figure 9
Surface settlement calculation curves.Figure9 shows that the maximum settlement of the foundation under the external wall nearby street will be more than 3 mm due to the left-line tunnel construction if the corresponding engineering measures are not taken. With the same construction technology and considering the influence of left tunnel excavation, the right tunnel construction induced maximum settlement of the foundation under the external wall nearby street will exceed 10 mm; it is difficult to ensure the safety of nearby buildings. So, the double-row high-pressure jet grouting cement piles were constructed between the external wall nearby street and the right tunnel; before the shield tunneling, the North Street 8# building pile diameter is 500 mm with 1000 mm center distance, and the South Street 8# building pile diameter is 450 mm with 950 mm center interval. Pile length is more than 1 m at the bottom of the shield tunnel, the water-cement ratio of pile is 0.8 ~1.2, the P.O.32.5 Portland cement is used, and the I20~I30 steel I-beam is inserted in piles. In the construction process, the inertia synchronous grouting slurry is replaced as hard synchronous grouting slurry, and the amount of cement is increased, the secondary grouting pressure is 0.3~0.4 MPa, and the grouting amount is 1.5 m3 per liner segment.Figure10 is settlement monitoring results of some sections; it shows that the left tunnel excavation induced about 5 mm ground surface subsidence on the right tunnel axis which is 11 m away from left tunnel axis and induced about 2 mm ground surface subsidence on the building which is 18 m away from left tunnel axis; the measured results of shield tunneling causing ground surface settlement are accordant with the analytical calculation results. It is revealed that the foundation reinforcement and grouting slurry adjustment measures are reasonable. By adopting active and passive engineering measures, the twin tunnels successfully underpassed the adjacent buildings, the shield tunneling-induced ground settlement is in standard control, and the nearby buildings safety is ensured.Figure 10
Surveying curves of some cross-sections.


## 5. Conclusions

Based on the assumption of homogeneous semi-infinite space, the initial two-way nonisobaric stress field of shallow-buried shield tunnel is divided into the uniform stress field and the one-way stress field; the semiempirical solution for shallow tunnel surface subsidence prediction is presented by using the analytical mechanic theory and definition the surrounding rock displacement boundary conditions of shield tunnel. The Matlab7.0 is chosen for programming to analyze the influence of the surrounding rock and soil mass properties on the ground surface settlement, and the corresponding surface settlement control method is proposed. The results of this research are used in the 11th Project of Beijing Metro Line 10 to predict and control the settlement of adjacent residential buildings affected by shield tunneling; the results show that the semiempirical solution is in good agreement with the in situ measured results; it can be used in the control of shield tunnel induced ground subsidence.



---

*Source: 3748658-2017-10-19.xml*