Open access

Prediction method of soil horizontal displacement caused by non-uniform distribution of disturbance force in shield construction

Authors: Bao-xin Jia, Zong-xian Gao [email protected], and Hong-yan HuaAuthors Info & Affiliations

Publication: FACETS

9 March 2023

https://doi.org/10.1139/facets-2022-0197

Abstract

Based on the Mindlin solution, this paper considers the influence of factors such as the non-uniform distribution of additional thrust of the cutter head influenced by lateral earth pressure in shield excavation, the non-uniform distribution of friction of shield shell influenced by soil softening and slurry spreading, and the non-uniform distribution of grouting pressure influenced by slurry spreading on the horizontal displacement of soil. The existing prediction formula is revised and verified by engineering examples. It is found that: affected by the shield construction disturbance force, the horizontal displacement behind the excavation surface is larger than that in front of the excavation surface, and the peak value of the horizontal displacement appears around the tunnel axis; through the verification of engineering case, when calculating the horizontal displacement in front of the excavation surface, the calculation results of both the modified formula and the original formula are in good agreement with the measured values, which can reflect the change trend of the measured horizontal displacement; when calculating the horizontal displacement behind the excavation surface, the calculation result of the existing formula has a great error due to the assumption of uniform distribution of disturbance force, which is different from the law of the measured result; the calculation result of the modified prediction formula is obviously in better agreement with the measured deformation value, and the error is smaller, which is more in line with the engineering reality.

1. Introduction

With the rapid development of urban rail transit, shield construction has become the main method of subway tunnel construction. The shield method has the advantages of high efficiency and little impact on the surrounding environment (Jin et al. 2022; Ye et al. 2022). At present, there have been relatively abundant studies on the prediction and influence of soil settlement caused by shield construction, but there are few studies on the horizontal displacement of soil. Horizontal displacement of deep soil is inevitably caused during shield construction, which may not only affect the stability of the tunnel working surface (Cheng et al. 2022; Cheng et al. 2021b), but also may adversely affect the integrity of the adjacent underground infrastructure (Burd et al. 2022; Klar 2022), so it is of great importance to be able to predict the horizontal displacement of the soil caused by the shield machine construction.

At present, many scholars have studied the horizontal displacement of soil caused by shield construction. Jiang et al. (2005, 2011) studied the variation law of the horizontal displacement of the stratum caused by the shield construction process under different stratum conditions. Standing and Selemetas (2013) studied the influence of existing underground structures on soil settlement, horizontal displacement, and pore pressure changes caused by shield excavation. Wu et al. (2018) studied the evolution law of the vertical displacement and horizontal displacement of the surrounding soil caused by double tunnel excavation. Qi et al. (2020) studied the correction effect of grouting on the horizontal displacement of soil in tunnel construction. Sheng et al. (2021) studied the spatial and temporal characteristics of the vertical displacement and horizontal displacement of the double shield tunnel on the existing pile foundation and soil. Zhu et al. (2021) discussed the displacement distribution law of surface and underground soil caused by underwater shield tunnel construction.

However, the above literatures all used numerical simulation and field monitoring methods to study the horizontal displacement of soil. To propose a reliable theoretical prediction method, Loganathan and Poulos (1998), Chou and Bobet (2002), and Cao et al. (2020) studied the effect of shield tail clearance on soil settlement and horizontal deformation through the analysis of a large number of engineering cases and deduced the prediction formula of settlement and horizontal deformation in the soil stability stage. Liang et al. (2015) derived the prediction formulas for the vertical displacement of the surface and the horizontal displacement of the deep soil during shield construction based on the Mindlin solution. Considering the incompressibility of rock and soil mass, based on the principle of plane strain, Tian et al. (2017) derived the calculation formula of the horizontal displacement of the soil along the cross section of the tunnel by using the Peck settlement formula.

However, the above studies either ignored the influence of disturbance forces on the horizontal displacement of soil during shield construction or simply assumed the disturbance forces as uniform distribution for calculation. And in the actual construction process, the disturbance forces such as additional thrust, friction force, and grouting pressure show uneven distribution due to the influence of different construction factors. Based on the Mindlin solution, this paper considers the influence of factors such as the non-uniform distribution of the additional thrust of the cutter head due to the lateral earth pressure of the stratum during shield excavation, the non-uniform distribution of shield shell friction due to the soil softening and slurry spreading, and the non-uniform distribution of grouting pressure due to the grout diffusion on the horizontal displacement of soil. The existing prediction formula is modified and verified by engineering examples.

2. Non-uniform distribution of shield disturbance force

2.1. Assumptions for shield construction calculations

The disturbance force on the surrounding soil during the shield construction process can be divided into three types: the additional thrust at the shield excavation surface, the friction between the shield shell and the surrounding soil, and the grouting pressure at the shield tail. The disturbance force is unevenly distributed under the influence of the lateral pressure of the formation, the softening of the soil, and the spreading of the slurry. For convenience, subsequent calculations satisfy the following assumptions:

The soil is an isotropic linear elastic semi-infinite body.

The shield machine is driven horizontally, regardless of the change of the driving posture.

Only the stratum deformation during the shield construction period is considered, and the stratum deformation caused by disturbed soil consolidation is not included in the calculation in this paper.

Shield propulsion only considers changes in space and does not consider changes in time.

2.2. Additional thrust distribution considering formation lateral pressure

In shield tunneling, there must be a pressure difference between the support pressure applied to the excavation face and the earth and water pressure of the stratum, that is, the shield frontal additional thrust. At present, when many studies analyze the influence of the frontal thrust of shield machine, for earth pressure balance (EPB) shield and slurry shield, the pressure value at the center of the cutter at the excavation face is generally selected for the study, and it is assumed that the cutter frontal thrust and slurry pressure are uniformly distributed on the excavation face, or the additional frontal thrust is simply assumed to be 20 kPa (Deng et al. 2022 ; Jia and Gao 2022; Mei et al. 2022). However, when the diameter of the shield machine is small, it can be assumed that the cutter frontal thrust and slurry pressure are uniformly distributed on the excavation surface, but as the diameter of shield increases, the pressure difference of lateral ground pressure formed at the top and bottom of the excavation face cannot be ignored (Chen et al. 2019; Cheng et al. 2021a). To ensure the stability of the excavation surface, the frontal thrust of the cutter head and the slurry pressure applied by the large-diameter shield machine are non-uniformly distributed under the influence of the differential pressure and the self-weight of the lateral earth pressure of the stratum, and the distribution of support force of shield excavation face is shown in Fig. 1.

Fig. 1.

To ensure the stability of the excavation surface, the frontal thrust and slurry pressure are generally greater than the lateral pressure of the formation. Therefore, the additional thrust that takes into account the uneven distribution of formation lateral pressure effects is as follows:

q = q_{1} - q_{2}

(1)

q_{2} = K_{0} γ (z - H)

(2)

where q is the additional thrust generated by the cutter head of the shield machine; q₁ is the frontal thrust of the cutter and mud pressure; q₂ is the lateral earth and water pressure; K₀ is the lateral pressure coefficient of the soil layer; γ is the soil layer weight; and H is the buried depth of the tunnel.

When there are no actual monitoring data, the additional thrust value can be obtained by the following methods according to Wang et al (Jiang et al. 2018; Wang 2009), when the EPB shield is used,

q = \frac{10.13 (1 - μ) E_{u} π v {(1 - ξ)}^{2}}{(1 + μ) (3 - 4 μ) D k_{u} w} + Δ q'

(3)

When the slurry shield is used,

q = q_{1} - q_{2} - k_{s l} γ_{s l} (z - H)

(4)

where Δq′ is the squeezing pressure generated by the incision cut into the soil, generally taken as 10∼25kP; μ is the Poisson’s ratio of the soil; E_u is the undrained elastic modulus (kPa) of the soil; v is the shield tunneling speed (cm/min); w is the rotational speed of the cutter; k_u is the amplitude of the closed part of the cutter; D is the diameter of the cutter (m); ξ is the opening rate of the cutter (%); γ_sl is the weight of slurry; and k_sl is the lateral pressure coefficient of slurry.

2.3. Friction distribution considering soil softening and slurry spreading

In the process of shield tunneling construction, the friction between the shield shell and the surrounding soil will cause the horizontal displacement of the soil. In most of the current researches, it is generally assumed that the friction force is uniformly distributed on the shield casing, which is inconsistent with the actual situation. Liang et al. (2015) pointed out that the friction between the shield machine and the surrounding soil is similar to the friction characteristics of the pile–soil interface, and the surrounding soil softens during shield tunneling, and the friction is residual friction. According to the shear force formula at the pile–soil interface given by Alonso et al. (2009), the friction force f at any position on the shield shell can be obtained.

{\begin{cases} f = β_{s} σ_{θ} \tan δ_{f} \\ σ_{θ} = σ_{v} s i n^{2} θ + σ_{h} c o s^{2} θ \\ σ_{v} = γ H - γ R s i n θ \\ σ_{h} = K_{0} σ_{v} \end{cases}

(5)

In the formula, β_s is the ratio of residual frictional resistance to ultimate frictional resistance; σ_θ is the radial normal stress received at any point of the shield machine casing; δ_f is the interface friction angle between the casing and the soil; and σ_v and σ_h are the pressure in vertical and horizontal directions.

In previous research, it is usually assumed that the shield machine is a cylinder with a constant diameter. However, in fact, the diameter of the cutter head of the shield machine is larger than the diameter of the shield tail. This difference is very small, generally about 0.4% of the diameter of the shield machine. During the construction of the shield machine, there is a joint gap s between the tail of the shield machine and the soil due to the difference in the diameter of the front and rear of the shield machine, as shown in Fig. 2. Bezuijen (2009) and Nagel and Meschke (2011) found that the grouting slurry will spread to the joint gap of the shield tail during shield construction, and the friction force of the shield tail will be less affected by the slurry. In the calculation, the friction force at the shield tail needs to be multiplied by the reduction factor λ. The grouting spread distance a at any position of the shield tail can be obtained by changing the grouting pressure decay formula of the shield tail.

Δ σ = 2 \frac{Δ R'}{R'} G

(6)

Δ P = 2 \frac{Δ a}{s} τ_{p}

(7)

Fig. 2.

In the formula, Δσ is the change of grouting pressure; R′ is the radius of the shield machine and the thickness of the slurry; ΔR′ is the change of the radius; G is the shear modulus of the soil; ΔP is the change of the grouting pressure after the slurry flows; Δa is the length increment along the tunneling direction; and τ₀ is the shear stress of shield tail grouting.

The calculation method is: divide the soil around the tunnel into independent areas of width Δa and use eq. (6) to calculate the joint width at the tail of the shield machine. This seam width can be used in eq. (7) to calculate the pressure change when grouting flows to the front of the shield machine at distance Δa. The resulting grouting pressure (original grouting pressure minus pressure change) is used to calculate changes in the next area and joint width, etc. For each calculation, it is checked whether the calculated seam width is still positive. When the joint width is 0, the soil is in contact with the shield machine, and the actual diameter of the shield machine in this area is used. For the specific calculation method of the grouting pressure at any position θ of the shield tail, see Section 2.4.

2.4. Grouting pressure distribution considering grouting diffusion

The grouting pressure of the shield tail has an important influence on the horizontal displacement of the soil. Most of the current studies have not considered the spatial non-uniform distribution of the grouting pressure. The synchronous grouting of shield tail is divided into two independent processes: formation and dissipation. The formation of grouting pressure occurs in the stage where the shield tail and the segment are pulled out. At this stage, the grouting hole releases the slurry, and in a very short time the slurry flows along the gap between the shield tail and the segment is filled in the circumferential direction. Due to the influence of the grouting pressure distribution in the grouting hole, the slurry density, and the slurry diffusion, the circumferential grouting pressure distribution is usually a non-uniform distribution pattern with small upper and big lower. The process of grouting pressure dissipation occurs within a certain range behind the shield tail. At this time, the new grouting liquid in the circumferential direction squeezes the grouted liquid, and the extrusion effect gradually decays with the increase of the distance from the shield tail until the grouting pressure reaches a balance with the water and soil pressure of the surrounding stratum. Therefore, the grouting pressure is generated at the shield tail and decays within a certain range behind the shield tail. This is quite different from the simplified treatment of grouting pressure in most existing studies. According to the analysis of the circumferential and longitudinal filling diffusion mechanism of the slurry during the simultaneous grouting process of the shield tail, the following basic assumptions are put forward to construct the theoretical model of its circumferential and longitudinal filling diffusion:

① The slurry is an incompressible and isotropic homogeneous fluid, ignoring the loss and property change of the slurry during the diffusion process;

② After the slurry flows out of the grouting pipe, it first fills the gap in the circumferential direction and then flows in the axial direction of the tunnel in the gap of the shield tail;

③ The axial pressure of the slurry at the same point is equal to the radial pressure;

④ Since the hoop filling stage usually takes only tens of seconds to complete, therefore, regardless of the time spent in the hoop filling stage, the moment when the slurry filling is completed along the hoop is set as the start time of grouting.

Circumferential distribution of grouting pressure

According to the characteristics of the grouting slurry and referring to the existing research, this paper assumes that the grouting slurry is Bingham fluid and analyzes it. The annular distribution formula of the grouting pressure P in the shield tail is

P = P_{i} + ρ_{p} g R_{g} (c o s α_{i} - c o s α) \pm Y (α_{i} - α)

(8)

where i is the number of the grouting hole; P_i is the grouting pressure at the number; “±” is “+” when the slurry spreads downward, and “–” otherwise; ρ_p is the density of the slurry; g is the acceleration of gravity; R_g is the outer diameter of the pipe; G is the angle between the grouting hole at the number and the vertical direction; and α_i is the slurry diffusion angle, when the slurry spreads upward, max α = (α_i + α_i–1)/2 and when the slurry spreads downward, max α = (α_i + α_i+1)/2.

Y^{3} - (\frac{3 R_{g} τ_{0}}{δ p} + \frac{12 R_{g} μ_{p} q_{p}}{b δ_{p}^{3}}) Y^{2} - \frac{4 {R_{g}}^{3} τ_{0}^{3}}{δ_{p}^{3}} = 0

(9)

In the formula, τ₀ is the static shear force of the grouting slurry; μ_p is the plastic viscosity coefficient of the slurry; b is the gap size of the shield tail; δ_p is the thickness of the slurry flow pancake, δ_p = v t; v is the shield tunneling speed; t is the time required for the slurry to fill the gap when the shield tail is synchronously grouted, generally 30∼100 s; and q_p is the interface flow, q_p ≈ πR_gbvm/n.

At present, 4 or 6 grouting holes are generally set for synchronous grouting of shield machine, and the applied grouting pressure is left and right symmetrical. The circumferential distribution of grouting pressure with different numbers of grouting holes is shown in Fig. 3.

Fig. 3.

Longitudinal distribution of grouting pressure

The longitudinal diffusion model of grouting slurry is shown in Fig. 4. Taking the longitudinal slurry-micro-element as the research object, and analyzing the force of the micro-element, the force balance equation can be obtained:

P_{s} d y - (P_{s} + d P) d y - τ d s + (τ + d τ) d s = 0

(10)

Fig. 4.

Let dP_s/ds = B, according to the boundary conditions, at the tail of the shield, S = 0, P_s = P, then the pressure distribution equation of the slurry along the longitudinal direction of the tunnel is

P_{s} = P - B S

(11)

In the formula, P_s is the longitudinal grouting pressure at the distance S from the shield tail, and it will not attenuate when the grouting pressure attenuates to equilibrium with the formation pressure; S is the distance from the shield tail, 0 < S < l; l is the grouting slurry diffusion length; B is the attenuation strength of the grouting pressure along the longitudinal direction; and τ is the shear force of the slurry.

B can be solved by Bingham’s fluid continuity equation, and the exact expression after the solution is

B = \sqrt[3]{\frac{4 τ_{0}^{3}}{{(\frac{12 μ_{p} q_{p}}{b} + 3 τ_{0} h^{2} - h^{3})}^{3}}} \approx \frac{12 π R_{g}^{2} μ_{p} v m}{δ_{p}^{3} n}

(12)

In the formula, h is the filling width of the slurry; m is the injection rate of the slurry, generally 1.1∼1.8; and n is the number of grouting holes.

Spatial distribution of grouting pressure

Combining eq. (1) and eq. (3), the formula for calculating the spatial distribution of grouting pressure can be obtained:

P_{s} = P_{i} + ρ_{p} g R_{g} (c o s α_{i} - c o s α) \pm Y (α_{i} - α) - B S

(13)

3. Calculation of horizontal displacement caused by shield construction

3.1. Mindlin basic solution

Mindlin (1936) deduced the calculation formulas of the vertical displacement and the horizontal displacement of any point (x′, y′, z′) in the space when a point (0, 0, c) in a semi-infinite elastic space is affected by the vertical concentration force P_v and the horizontal concentration force P_h. Under the action of the vertical concentration force P_v and the horizontal concentration force P_h, the horizontal displacements v_v and v_h of any point in the soil are shown in eqs. (14) and (15).

v_{v} = \frac{P_{v} \sqrt{{x'}^{2} + {y'}^{2}} \cos ϕ}{16 π G (1 - μ)} [\frac{z' - c}{{R_{1}}^{3}} + \frac{(3 - 4 μ) (z' - c)}{{R_{2}}^{3}} - \frac{4 (1 - μ) (1 - 2 μ)}{R_{2} (R_{2} + z' + c)} + \frac{6 c z' {(z' + c)}^{2}}{{R_{2}}^{5}}]

(14)

v_{h} = \frac{P_{h}}{16 π G (1 - μ)} {\frac{3 - 4 μ}{R_{1}} + \frac{1}{R_{2}} + \frac{{y'}^{2}}{{R_{1}}^{3}} + \frac{(3 - 4 μ) {y'}^{2}}{{R_{2}}^{3}} + \frac{2 z' c}{{R_{2}}^{3}} (1 - \frac{3 {y'}^{2}}{{R_{2}}^{2}}) + \frac{4 (1 - μ) (1 - 2 μ)}{(R_{2} + z' + c)} [1 - \frac{{y'}^{2}}{R_{2} (R_{2} + z' + c)}]}

(15)

where G is the shear modulus of the soil (kPa), and μ is the Poisson’s ratio.

{\begin{cases} R_{1} = \sqrt{{x'}^{2} + {y'}^{2} + {(z' - c)}^{2}} \\ R_{2} = \sqrt{{x'}^{2} + {y'}^{2} + {(z' + c)}^{2}} \\ c = H - r s i n θ \end{cases}

(16)

Set the space coordinate xyz, and the local coordinate x′y′z′, the corresponding coordinate axes of the two coordinate systems are parallel to each other, and the offset distances of the local coordinate system are l, m, and n, respectively.

{\begin{cases} x' = x - l \\ y' = y - m \\ z' = z - n \end{cases}

(17)

Substituting eq. (17) into eqs. (14)–(16) can obtain the vertical displacement and horizontal displacement at any point (x, y, z) in the spatial coordinates of the soil.

3.2. Horizontal displacement caused by uneven distribution of additional thrust on the excavation face

When the shield machine is excavating, the calculation model of soil deformation caused by the additional thrust of the excavation surface is shown in Fig. 5. The area of any element in the excavation surface is dA = rdrdθ, where r is the distance from the element to the center of the excavation surface and θ is the angle between the element and the horizontal plane of the center of the excavation surface, l = 0, m = rcosθ, n = 0.

Fig. 5.

Calculation formula of horizontal displacement caused by uneven distribution of additional thrust on excavation face during shield tunneling construction is as follows:

v_{h} = \int_{0}^{R} \int_{0}^{2 π} \frac{qrx' y'}{16 π G (1 - μ)} [\frac{1}{{R_{1}}^{3}} + \frac{3 - 4 μ}{{R_{2}}^{3}} - \frac{4 (1 - μ) (1 - 2 μ)}{R_{2} {(R_{2} + z' + c)}^{2}} - \frac{6 c z'}{{R_{2}}^{5}}] d r d θ

(18)

3.3. Horizontal displacement caused by uneven distribution of shield friction

The calculation model of the formation deformation caused by the friction force f between the shield machine shell and the surrounding soil is shown in Fig. 6. Let the length of the shield from the excavation cutter head to the shield tail be L, and the area of any element on the surface of the shell is dA = Rdθdl, where l is the axial distance from the element to the excavation surface, and the concentrated force is dP_h = f Rdθdl.

Fig. 6.

The horizontal displacement caused by the frictional force of the shield shell during the shield tunneling construction can be calculated by the integral of eq. (14).

d v_{f} = \int_{0}^{L} \int_{0}^{2 π} \frac{f R x' y'}{16 π G (1 - μ)} [\frac{1}{{R_{1}}^{3}} + \frac{3 - 4 μ}{{R_{2}}^{3}} - \frac{4 (1 - μ) (1 - 2 μ)}{R_{2} {(R_{2} + z' + c)}^{2}} - \frac{6 c z'}{{R_{2}}^{5}}] d l d θ

(19)

From the analysis in Section 2.2, it can be seen that due to the influence of the shield tail slurry spreading, the friction force needs to be calculated in different regions.

v_{f} = \int_{0}^{L - c} \int_{0}^{2 π} d v_{f_{1}} + \int_{L - c}^{L} \int_{0}^{2 π} d v_{f_{2}}

(20)

3.4. Horizontal displacement caused by uneven distribution of grouting pressure in shield tail

As shown in Fig. 7, the length of the shield from the shield tail to the grouting end is S, the grouting pressure at the shield tail is p, the grouting pressure decays along the grouting length, and the grouting pressure at the distance S away from the shield tail will decrease to 0, any element dA = Rdθds, and its concentrated force dp can be decomposed into horizontal component force dp_h and vertical component force dp_v. Because the influence of vertical component force is small and can be ignored, this paper only considers formation deformation caused by horizontal component force.

Fig. 7.

Set l = −L−s, m = Rcosθ, n = 0, and substitute them into eq. (15) to transform the coordinate system. The horizontal displacement caused by grouting pressure in shield tunneling construction can be calculated by the integral calculation of eq. (15).

v_{p} = \int_{0}^{S} \int_{0}^{2 π} \frac{(P \cos θ) R}{16 π G (1 - μ)} {\frac{3 - 4 μ}{R_{1}} + \frac{{(y - R \cos θ)}^{2}}{{R_{1}}^{3}} + \frac{1}{R_{2}} + \frac{(3 - 4 μ) {(y - R \cos θ)}^{2}}{{R_{2}}^{3}} + \frac{2 z (H - R \sin θ)}{{R_{2}}^{3}} [1 - \frac{3 {(y - R \cos θ)}^{2}}{{R_{2}}^{2}}] + \frac{4 (1 - μ) (1 - 2 μ)}{(R_{2} + z' + c)} [1 - \frac{{(y - R \cos θ)}^{2}}{R_{2} (R_{2} + z' + c)}]} d s d θ

(21)

3.5. Horizontal displacement caused by formation loss

Horizontal displacement formula of the soil above caused by tunnel construction is as follows:

v_{s} = - V_{l} R^{2} y {\frac{1}{y^{2} + {(H - z)}^{2}} + \frac{3 - 4 μ}{y^{2} + {(H + z)}^{2}} - \frac{4 z (H + z)}{{[y^{2} + {(H + z)}^{2}]}^{2}}} e x p {- [\frac{1.38 y^{2}}{{(H + R)}^{2}} + \frac{0.69 z^{2}}{H^{2}}]}

(22)

where V_l is the formation loss rate. A three-dimensional calculation formula for the formation loss along the direction of tunnel excavation is as follows:

V_{l} (x) = \frac{V_{l}}{2} (1 - \frac{x}{\sqrt{x^{2} + y^{2} + H^{2}}})

(23)

Since the stratum loss at the excavation face is basically less than half of the maximum stratum loss at the rear, based on this consideration, a formula for the horizontal displacement of soil mass caused by soil mass loss under three-dimensional conditions is given.

v_{s} = - \frac{V_{l} R^{2} y}{2} (1 - \frac{x}{\sqrt{x^{2} + y^{2} + H^{2}}}) {\frac{1}{y^{2} + {(H - z)}^{2}} + \frac{3 - 4 μ}{y^{2} + {(H + z)}^{2}} - \frac{4 z (H + z)}{{[y^{2} + {(H + z)}^{2}]}^{2}}} e x p {- [\frac{1.38 y^{2}}{{(H + R)}^{2}} + \frac{0.69 z^{2}}{H^{2}}]}

(24)

3.6. Total displacement of stratum caused by shield advancement

The above displacement formulas are superimposed to obtain the solution of the horizontal displacement of the soil caused by the advancement of the shield.

v = v_{q} + v_{f} + v_{p} + v_{s}

(25)

4. Engineering case analysis

To verify the validity of the proposed calculation method, this paper selects two engineering cases and compares the calculation results of different calculation methods to study the influence of uneven distribution of disturbance force on the horizontal displacement of soil caused by shield construction under the influence of multiple factors. Positive values indicate that the soil moves in the direction away from the tunnel, otherwise vice versa.

4.1. Engineering Case 1

Taking a shield tunnel project of Shanghai Metro Line 2 as an example, the stratum in this area has the characteristics of sensitivity, saturation, and low permeability, and the EPB shield machine is used for construction. This section passes through the prosperous commercial area of Shanghai and has strict requirements on stratum deformation, and shield tunneling is relatively cautious. For the specific soil profile and parameters, please refer to the literature Lee et al. (1999) and Liang et al. (2015). The required calculation parameters are as follows: H = 15 m; γ = 17.8 KN/m³; K₀ = 0.57; G = 3.5 MPa; μ = 0.35; v = 12 m/d; L = 6.24 m; R = 3.17 m; R_g = 3.15 m; S = 1 m; n = 4; q = 58.9 kPa; β_s = 0.85; δ_f = 22°; λ = 0.85; p₁ = p₄ = 400 kPa; p₂ = p₃ = 100 kPa; ρ_p = 2190 kg/m³; μ_p = 1.8 Pa·s; m = 1.2; τ₀ = 1.8 × 10⁻² kPa; t = 100 s; and V_s = 0.14%.

Two cross-section positions 1D in front and 1D behind the excavation of the shield machine are selected for research. The horizontal distance from the monitoring point to the tunnel axis is 6.2 m. The horizontal deformation of the formation is calculated for both sections. Compare the calculated value with the field measured value and the calculation result of the prediction formula in the literature, as shown in Figs. 8–9. The horizontal deformation of the soil behind the excavation face is not affected by the additional thrust, and the horizontal displacement value is positive when it is away from the tunnel side.

Affected by the disturbance force of shield construction, the horizontal displacement direction of the deep soil body is the direction away from the shield tunnel, and the horizontal displacement behind the excavation face is greater than that in front of the excavation face, and the peak value of the horizontal displacement appears near the tunnel axis.

It can be found from Fig. 8 that when calculating the horizontal displacement in front of the excavation face, the calculated results considering the non-uniform distribution of disturbance forces agree well with the measured values, and the error between the calculated results assuming uniform distribution of disturbance forces and the measured values is large, but within acceptable limits. The horizontal displacement of the modified prediction formula converges faster on both sides of the tunnel axis.

It can be seen from Fig. 9 that when calculating the horizontal displacement behind the excavation face, the existing prediction formula assumes that the disturbance force is uniformly distributed, and the calculation result has a large error. The final calculation result is only about half of the measured value, and the peak value of the horizontal displacement appears at the tunnel axis, while the peak value of the measured value is above the axis, which is different from the measured result. The correction formula takes into account the uneven distribution of different disturbance forces under the influence of multiple factors. The final calculated results are in better agreement with the measured deformation values, and the error is smaller. When the distance is far from the tunnel axis, the error between the calculation curve of the existing prediction formula and the modified prediction formula and the measured value is large, which may be the influence of the diffusion of the tunnel grouting slurry in the soil. When the mechanical properties of the upper and lower soil layers are quite different, the horizontal displacement of the soil mass will also be affected.

It can be found from Figs. 8–9 that the effect of the modified additional thrust on the horizontal displacement of the soil is greater. Considering the influence of grouting spread, the final calculation result of friction force is smaller and more realistic. The spatiotemporal uneven distribution of grouting pressure has little effect on the horizontal displacement of the soil in front of the shield, which is basically the same as the calculation result of the existing formula, and has a great influence on the soil deformation behind the excavation face. When the space-time is not uniformly distributed, it will cause large error to the final result.

Fig. 8.

Fig. 9.

4.2. Engineering Case 2

Taking the Luoxi–Nanzhou section of Guangzhou Rail Transit Line No. 2 and No. 8 as an example, the tunnel is constructed with two earth pressure shield machines. The landform of the project site belongs to the alluvial plain of the Pearl River Delta, and the ground is flat. For details of soil profile and parameters, please refer to Jiang et al. (2011). The required calculation parameters are as follows: H = 29.3 m; γ = 16.14 KN/m³; K₀ = 0.53; G = 1.8 MPa; E = 4.71 MPa; μ = 0.28; v = 18 m/d; L = 7.5 m; R = 3.13 m; R_g = 3.05 m; S = 1.5 m; n = 4; q = 50 kPa; p₁ = p₂ = p₃ = p₄ = 350 kPa; δ_f = 4.5°; β_s = 0.9; λ = 0.85; ρ_p = 2150 kg/m³; μ_p = 1.7 Pa·s; m = 1.2; τ₀ = 1.8 × 10⁻² kPa; t = 100 s; and V_s = 0.99%.

Two cross-section positions 4.5 m in front and 18 m behind the excavation face of the shield machine are selected for research. The horizontal distance from the monitoring point to the tunnel axis is 6 m. Use the revised prediction formula and the existing prediction formula to calculate the horizontal displacement of the stratum of the two sections, and compare the calculation results with the field measured values, as shown in Figs. 10–11.

The horizontal displacement law of soil caused by shield construction is consistent with Case 1. When calculating the horizontal displacement in front of the excavation, the calculation results considering the non-uniform distribution of the disturbance force are in better agreement with the measured values, which can reflect the change trend of the measured horizontal displacement.

It can be seen from Fig. 11 that when calculating the horizontal displacement behind the excavation face, the calculation results of the existing prediction formula have great error. The calculation results show that the soil moves to the side of the tunnel, which is completely opposite to the measured value. In the calculation of the modified prediction formula, the uneven distribution of different disturbance forces under the influence of multiple factors is considered, and the final calculation result is in better agreement with the measured deformation value, and the error is smaller.

Fig. 10.

Fig. 11.

5. Conclusion

Based on the Mindlin solution theory, this paper considers the influence of the non-uniform distribution of disturbance force on the horizontal displacement of soil and modifies the existing prediction formula. It is verified by engineering cases, which proves the validity and reliability of the calculation method in this paper. The main conclusions are as follows:

Affected by the disturbance force of shield construction, the horizontal displacement direction of the deep soil body is away from the direction of the shield tunnel, and the horizontal displacement behind the excavation face is greater than that in front of the excavation face, and the peak value of the horizontal displacement occurs near the tunnel axis or obliquely above.

After the verification of engineering cases, when calculating the horizontal displacement in front of the excavation face, the calculation results of the revised formula and the original formula are in good agreement with the measured values, which can reflect the change trend of the horizontal displacement of the soil, and the calculation results of the modified formula are more accurate. When calculating the horizontal displacement behind the excavation face, the existing formula assumes that the disturbance force is uniformly distributed, and the error of the calculation result is very large, which is different from the law of the measured result. The calculation result of the revised prediction formula is in better agreement with the measured deformation value, the error is smaller, and it is more in line with the actual engineering.

The uneven distribution of the additional thrust has little effect on the horizontal displacement of the soil. The spatial and temporal uneven distribution of grouting pressure has little influence on the horizontal displacement of the soil in front of the shield but has a great influence on the deformation of the soil behind the excavation face.

Acknowledgements

This work was financially supported by the Natural Science Foundation of China (51774173) and the Liaoning Province “rejuvenating Liaoning talents plan” project (XLYC207163).

References

Alonso E, Josa A, and Villalba A. 2009. Negative skin friction on piles: A simplified analysis and prediction procedure. Géotechnique, 34(3): 341–357.

LOGIN TO YOUR ACCOUNT

Create a new account

Request Username

Change Password

Your password must have 8 characters or more and contain 3 of the following:

Password Changed Successfully

Verify Phone

Congrats!

Abstract

1. Introduction

2. Non-uniform distribution of shield disturbance force

2.1. Assumptions for shield construction calculations

2.2. Additional thrust distribution considering formation lateral pressure

2.3. Friction distribution considering soil softening and slurry spreading

2.4. Grouting pressure distribution considering grouting diffusion

3. Calculation of horizontal displacement caused by shield construction

3.1. Mindlin basic solution

3.2. Horizontal displacement caused by uneven distribution of additional thrust on the excavation face

3.3. Horizontal displacement caused by uneven distribution of shield friction

3.4. Horizontal displacement caused by uneven distribution of grouting pressure in shield tail

3.5. Horizontal displacement caused by formation loss

3.6. Total displacement of stratum caused by shield advancement

4. Engineering case analysis

4.1. Engineering Case 1

4.2. Engineering Case 2

5. Conclusion

Acknowledgements

References

Information

Published In

History

Copyright

Data Availability Statement

Key Words

Sections

Subjects

Plain Language Summary

Authors

Affiliations

Author Contributions

Competing Interests

Metrics

Other Metrics

Citations

Cite As

Export Citations

View options

PDF

Get Access

Media

Other

Share

Share the article link

Share on social media