Next Article in Journal
COVID-19 Diagnosis from Crowdsourced Cough Sound Data
Next Article in Special Issue
Algorithm for an Effective Ratio of the Transverse Bending Rigidity Based on the Segment Joint Bending Stiffness
Previous Article in Journal
Enhancement of Astaxanthin and Fatty Acid Production in Haematococcus pluvialis Using Strigolactone
Previous Article in Special Issue
Effect of Excavation Blasting in the Arch Cover Method on Adjacent Existing Pipelines in a Subway Station
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Estimating Volume Loss for Shield-Driven Tunnels Based on the Principle of Minimum Total Potential Energy

Key Laboratory of Urban Underground Engineering of Ministry of Education, Beijing Jiaotong University, Beijing 100044, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2022, 12(4), 1794; https://doi.org/10.3390/app12041794
Submission received: 24 December 2021 / Revised: 28 January 2022 / Accepted: 29 January 2022 / Published: 9 February 2022
(This article belongs to the Special Issue Deep Rock Mass Engineering: Excavation, Monitoring, and Control)

Abstract

:
Volume loss is an important method to estimate ground movement during tunnelling. However, volume loss is usually estimated by empirical methods, especially for volume loss at the tunnel face, which is a three-dimensional problem. Based on the principle of minimum total potential energy, we proposed a semi-analytical method to predict the volume loss at the tunnel face and ground surface. The proposed method provides a more direct way of estimating volume loss at the tunnel face from an energy point of view. Moreover, a new deformation mechanism was designed to describe the ground movement before the tunnel face. Based on the proposed method, we investigated the influence of the support pressure, tunnel diameter, and tunnel depth on the volume loss at the tunnel face, and other parameters related to surface subsidence. The volume loss at the tunnel face decreased with the increase in the support pressure ratio and the slurry weight. The volume loss at ground surface was generally smaller than the volume loss at the tunnel face due to the soil compression during ground movement. The bigger the tunnel diameter, the bigger the volume loss at the tunnel face. However, the volume loss at ground surface may not increase with the increase in tunnel depth, because of the soil arching effect. Moreover, the deeper the tunnel, the more obvious the influence of the support pressure ratio on volume loss. Similarly, the bigger the tunnel diameter, the more obvious the influence of the slurry weight on the volume loss.

1. Introduction

The slurry pressure balance shield tunnel is widely adopted in practice, and ground movement is inevitably induced during tunnel excavation. The induced ground movement may cause serious damage to existing buildings and subterranean facilities. Therefore, the estimation of ground movement is important in the design stage of tunnelling. Many researchers have contributed to analyzing the ground movement by different methods. Those methods can be divided into three categories: the analytical method [1,2,3,4,5,6,7,8], the empirical method [9,10,11], and numerical simulation [12,13,14].
Volume loss plays an important role in the analysis of ground movement. The effect of the tail grout on the volume loss was studied [15]. Different ways for determining the volume loss are summarized [11]. Furthermore, many researchers have investigated the volume loss in two-dimensional problems by assuming a zero-volume change in soil (i.e., undrained volume loss). The undrained volume loss is equal to ground loss, which is defined as the ratio of shrinkage deformation to the area of the tunnel (sometimes the undrained volume loss is called “ground loss”, in some literature). The analytical solution for the tunnel excavation as the function of the ground loss and the tunnel ovalization parameter are developed [16]. The deformation was divided into the ground loss, ovalization, and vertical movement [17]. A lot of researchers have investigated undrained volume loss with gap parameters, which was developed by considering the 3D elastoplastic deformation at the tunnel face [18]. The procedure for estimating the gap parameter was supported with 14 case monitoring data [19]. Loganathan and Poulos [3] incorporated the gap parameter into the closed-form solution derived by Verruijt and Booker [16]. Based on the work performed by Loganathan and Poulos [3], A solution for predicting the ground movement in multi-layered clayed soil was proposed [20]. The above research assumed a plane strain condition, which was limited to a two-dimension solution.
Some researchers have devoted their study to investigating the ground deformation and volume loss before the tunnel face. Based on the field monitoring data, Different relationships between the stability number, which is used in the stability analysis for tunnelling, and the volume loss before the tunnel face are established [21,22,23,24].
In recent years, the energy method has been developed rapidly in the analysis of ground deformation due to tunnelling. The two-dimensional ground movement was predicted with the balance between energy and work, ignoring soil volume change [25]. Based on the principle of energy conservation, an innovative method to predict the two-dimensional undrained volume loss were proposed [26]. Plamer and Mair [27] gave a solution for estimating volume loss in the elastic ground on the basis of the reciprocal theorem. As for volume loss at the tunnel face, Klar and Klein [28] estimated the face ‘take’ undrained volume loss with the use of a mobilizable strength design (MSD).
In this paper, we develop a method to predict volume loss at the tunnel face based on the principle of minimum total potential energy. Moreover, a new deformation mechanism is created based on a proposed deformation mechanism [29] to consider the soil volume change. Finally, we analyze the influence of support pressure supplied by slurry, tunnel diameter, and tunnel depth on volume loss. The proposed method provides a more direct way of estimating volume loss at the tunnel face from the energy point of view. The soil volume changes during ground movement are considered in the proposed method.

2. Problem Description

According to the causes of volume loss in slurry pressure balance shield tunnelling [10,24,30], volume loss can be categorized into three components, as shown in Figure 1:
Vlf: Volume loss at the tunnel face. The stability of the working face is achieved by the support pressure provided by a pressurized mixture of bentonite or clay and water (i.e., slurry). However, the support pressure is different from the lateral earth pressure. In general, the soil tends to move towards the excavation chamber and the ground movement will cause the subsidence of the ground surface before the face.
Vls: Volume loss along the shield. To facilitate the advancement of the shield machine, the radius at the head of the shield is bigger than that at the tail. This will cause the overcutting of the soil and the subsidence of the ground surface above the shield.
Vlt: Volume loss at the tail. The voids between the excavation surface and precast segment cause volume loss at the tail. Volume loss at the tail can be compensated by the appropriate grout at the tail.
In this paper, we focus on investigating the volume loss at the tunnel face (VLf) under different support pressure supplied by the slurry. The other components of volume loss are assumed to be zero, and the displacement of the tunnel opening is fixed. The same assumption has been adopted by previous researchers [31,32,33]. VL refers to the volume loss at the ground surface which is defined as the ratio of the integral of the vertical displacement over the three-dimensional settlement trough to the cross-sectional area of the tunnel per unit length. Vlf is defined as the ratio of the integral of the horizontal displacement over the tunnel face to the cross-sectional area of the tunnel per unit length. Due to difficulties in considering the volume changes of soil, most researchers assume that the moved soil has a constant volume and VL is equal to Vlf. In our method, we simply consider the volume change.
As shown in Figure 1, we assumed a circular tunnel of diameter D at finite depth H. The depth of the cover soil is denoted by C. The stress of the soil increases with the depth, and the Buoyant unit weight of the soil is ρ′. The lateral pressure coefficient of the soil is k. The initial horizontal earth pressure at the crown of the tunnel is noted by σu, as shown in Figure 1.
When tunnelling, the soil ahead of the face tends to move into the excavation chamber. The support pressure supplied by the slurry is important in controlling the ground movement. As shown in Figure 1, the support pressure at the tunnel crown is denoted by Pu, and it can be given by Equation (1):
P u = α · σ w = α · k σ u + σ w
where α is the support pressure ratio; σu is the sum of the lateral earth pressure ( σ u ) and the static pore pressure (σw) at the tunnel crown.
Due to the weight of the slurry, the support pressure (P) linearly increases along h and increases to Pb at the tunnel invert (Equation (2)). The support pressure can be given as follows:
P = P u + γ h
where γ is the unit weight of the slurry, and h is the distance from the tunnel crown.

3. Principle of Minimum Total Potential Energy

The principle of minimum total potential energy is widely utilized in different industries, including chemistry, physics, and engineering. Moreover, many numerical simulation methods are based on it. The principle of minimum total potential energy indicates that an unbalanced system will finally translate to a balanced system with the minimum total potential energy. In order to utilize the minimum total potential energy principle in the proposed problem, the start and end of the deformation process first have to be determined. As shown in Figure 2, the soil is stable in the initial unbalanced state. When tunnelling, the support pressure supplied by slurry cannot substitute the original earth pressure, and face stress relief happens. As a result, the soil begins to move into the excavation chamber. In our model, the support pressure was assumed to be constant during the deformation process. At last, the system was stable when the total potential energy reached the minimum potential energy in the final equilibrium state. So, we determined the final equilibrium by finding the state with the minimum total potential energy. In other words, the parameters in the deformation mechanism can be determined by minimizing the total potential energy of the studied system.
In our model, the moved soil is regarded as the studied system in our method. The potential energy of the studied system contains three components: elastic strain energy stored in the deformed system (Us), the potential energy of gravity (Vg), and the potential energy of the support pressure (Vp). Vg and Vp belong to the potential energy of external force associated with the external force (V).
The elastic strain energy can be calculated by the strain of the soil. In elastic mechanics, the strain energy density (strain energy of unit volume) can be calculated as follows:
U s 0 = E 2 ( 1 + μ ) μ 1 2 μ ( ε x + ε y + ε z ) 2 + ( ε x 2 + ε y 2 + ε z 2 ) + 1 2 γ x y 2 + γ x z 2 + γ y z 2
where εx, εy, εz, εxy, εxz, and εyx, are the strain in the Cartesian coordinate system as shown in Figure 2, E is the elastic modulus of the soil, and μ is Poisson’s ratio of the soil. Furthermore, the soil is assumed to behave elastically. The strain energy can be calculated by integrating Us0 over the moved soil:
V P = N U s 0 d n
where N is the volume of the moved soil in the deformation mechanism.
The external force consists of the support pressure and gravity. The potential energy of the external force is generally measured by the work performed by the external force. The increase in the potential energy of the external force is equal to the minus work performed by the external force. So, the potential energy of the support pressure and gravity can be given as follows:
V P = S ( P γ w ) u f d s V g = N g w d n
where uf is the axial displacement on the face, w is the vertical displacement of the soil, and S is the area of the face.
Some parameters for the calculation of potential energy (w, uf, and the strain of soil) are tightly related to the deformation mechanism. Thus, these parameters are introduced in Section 4.2 after the deformation mechanism introduction.

4. Deformation Mechanism

4.1. Deformaiton Mechanism

Two different continuous velocity fields for the collapse of a pressurized tunnel face were constructed [29]. The first velocity field, M1, assumed that the point with maximum velocity was located at the center of the tunnel face. The second velocity field was more suitable for engineer projects, because it assumed that the point with maximum velocity was located below the center of the tunnel face. However, the M2 velocity field was designed for collapse, and ignored the volume change of soil. A new deformation mechanism (Figure 3) was designed based on the modification of the M2 velocity field. The new deformation mechanism contains several parameters which decide the geometry of the new deformation mechanism. The parameters depend on the mechanical properties of the soil.
As shown in Figure 4, the maximum displacement on the tunnel face is Smf and is located at point E. The maximum velocity flow line is a quarter of an ellipse, with respect to the origin O. The semi-major axis of the ellipse is Rf = C + D/2 + L1. L1 is the distance between the center of the face and point E. The semi-minor axis of the ellipse is 0.75 Rf. For a better description of the deformation mechanism, the curvilinear coordinate system (r, θ, β) is adopted.
The new deformation mechanism assumes that the soil within the boundary flows towards the face in a toric manner, and the soil outside the boundary is rigid. The boundary is constructed by the periphery of the different circular cross-section, Πβ, as shown in Figure 5
Rmax is the radius of the circular cross-section, and it can be given by Equation (6):
R max ( θ ,   β ) = D · Ω 1 Ω 2 2 Ω 3
Ω 1 ( β ) = ( 1 + 2 L ¯ 1 ) ( π 2 β ) + 3 β R ¯ f
Ω 2 ( θ , β ) = 2 L ¯ cos θ + 1 2 L ¯ 2 + 2 L ¯ 2 cos 2 θ
Ω 3 = ( 1 + 2 L ¯ 1 ) π
L ¯ = L ¯ 1 l β
where [ · ] ¯ means normalized with respect to tunnel diameter (such as C ¯ = C / D ). L indicates the distance between the point with maximum displacement and the center of the circular cross-section.
On the Πβ plane, the orthoradial displacement component (vr) is assumed to be zero. The axial displacement component (vβ) is assumed to decrease in a parabolic way to vanish at the boundary, as shown in Figure 4. Equation (11) gives the orthoradial displacement of points (r, θ, β) on the plane Πβ:
v β ( r ,   θ ,   β ) = v m ( β ) 1 r 2 R max 2 ( θ ,   β )
where vm(β) is the orthoradial displacement on the maximum velocity flow line. On the face, the displacement can be expressed as follow:
v β ( r ,   θ ,   0 ) = S m f 1 r 2 R max 2 ( θ ,   0 )
where Smf is the maximum displacement on the face. Furthermore, the flux across the tunnel face (qsf) can be acquired by integral Equation (13) over the tunnel face.
q s f = S m f · π D 2 8
In previous deformation mechanism [29], the assumption of zero volume change was made and the flux through any Πβ plane was the same. However, as drainage or volume loss always happens in the sand, we assume that the flux across plane Πβ linearly increases (or decreases) with β:
q s = q s f ( k β + 1 ) = S m f · π D 2 8 ( k β + 1 )
where k is used to describe the extent of the volume change of the soil, and qs denotes the flux through plane Πβ. At the tunnel face, qs is equal to qsf. qs equal to Qs at the ground surface.
q s = 0 2 π 0 R max ( θ , β ) v β ( r ,   θ ,   β ) r d r d θ = 0 2 π 0 R max ( θ , β ) v m ( β ) 1 r 2 R max 2 ( θ ,   β ) r d r d θ = v m ( β ) 0 2 π 0 R max ( θ , β ) 1 r 2 R max 2 ( θ ,   β ) r d r d θ = v m ( β ) π D 2 Ω 1 2 8 Ω 3 2
So, combining Equations (14) and (15), the maximum orthoradial displacement on any plane can be given by Equation (16).
v m ( β ) = S mf ( k β + 1 ) Ω 3 2 Ω 1 2
So, vβ can be derived as follows.
v β ( r ,   θ ,   β ) = S mf ( k β + 1 ) Ω 3 2 Ω 1 2 1 r 2 R max 2 ( θ ,   β )
In solid mechanics, there are two expressions for the volumetric strain (Θ): the ratio of the volume change (∆V) to the original volume (V), and the summation of the normal strains in three mutually perpendicular directions, as shown in Equation (18).
Θ = Δ V V = ε r + ε θ + ε β
Moreover, the volume change (∆V) of the Πβ plane with dβ can be calculated by taking the derivative of qs as shown in Equation (19).
Δ V = q s β · d β = q s f k · d β
The original volume of the Πβ plane with dβ can be given as follows:
V = d β · R f · D s ( β ) 2 π / 4
where DS(β) is the diameter of the Πβ plane, and can be expressed as follows.
D s ( β ) = R max ( θ = 0 ,   β ) + R max ( θ = π 2 ,   β ) = Ω 1 D Ω 3
So, we have:
Θ = S mf · k Ω 3 2 2 R f Ω 1 2
The second expression of the volumetric volume shown in Equation (18) can be calculated by the strain introduced in Section 4.2. The volumetric volume can be also given as follows.
v r r + 1 r v r + 1 R β r cos θ ( υ r cos θ + υ β β ) = Θ
After simplification of Equation (23), we have:
v r r ¯ + ( R ¯ β 2 r ¯ cos θ ) r ¯ ( R ¯ β r ¯ cos θ ) v r = 1 ( R ¯ β r ¯ cos θ ) v β β + S mf · k Ω 3 2 2 R ¯ β Ω 1 2
The derivative of orthoradial displacement vβ with respect to β can be given as follows:
v β β = S mf Ω 3 2 Ω 1 2 k 1 4 Ω 3 2 r ¯ 2 Ω 1 4 Ω 2 2 + ( k β + 1 ) 2 Ω 4 4 Ω 3 2 Ω 5 r ¯ 2 Ω 1 5 Ω 2 3 = S mf Ω 3 2 Ω 1 2 4 Ω 3 2 Ω 1 4 Ω 2 2 k + ( k β + 1 ) Ω 5 Ω 1 Ω 2 r ¯ 2 + k 2 ( k β + 1 ) Ω 4
Ω 5 ( θ , β ) = 4 Ω 2 Ω 4 + 3 Ω 1 cos θ + 2 L ¯ sin 2 θ 1 2 L ¯ 2 + 2 L ¯ 2 cos 2 θ
Equation (24) is an ordinary first-order linear differential equation. In order to simplify the calculation, Equation (24) can be expressed as follows:
v r r ¯ + ( R ¯ β 2 r ¯ cos θ ) r ¯ ( R ¯ β r ¯ cos θ ) v r = 1 ( R ¯ β r ¯ cos θ ) A r ¯ 2 + B + Θ
where A and B can be given as follows:
A = 4 S mf Ω 3 4 Ω 1 6 Ω 2 2 k + ( k β + 1 ) Ω 5 Ω 1 Ω 2 B = S mf Ω 3 2 Ω 1 2 k 2 ( k β + 1 ) Ω 4
Solving it with boundary condition vr(Rmax(θ, β), θ, β) = 0, we have:
v r ( r , θ , β ) = 6 ( B Θ R ¯ β ) ( r ¯ R max ) 2 + 4 cos θ · ( r ¯ R max ) 3 + 3 A ( r ¯ R max ) 4 12 ( R ¯ β r ¯ cos θ )

4.2. Calculation of Potential Energy

(a)
Strain energy
Because the strain of the soil is derived from the deformation mechanism, and the deformation mechanism in this paper is expressed in a curvilinear coordinate system, we have to convert the strain in the curvilinear coordinate system to the strain in the Cartesian coordinate system. Equation (30) gives the relationship between the curvilinear coordinate system and the Cartesian coordinate system.
x = r sin θ y = ( R β r cos θ ) cos β z = ( R β r cos θ ) sin β
R β = cos 2 β / R f 2 + cos 2 β / ( 0.75 R f ) 2 0.5
The conversion can be easily given by tensor calculation, as shown in Equations (32) and (33).
ε x ε x y ε x z ε y x ε y ε y z ε z x ε z y ε z = P ε r ε r y ε r z ε θ r ε θ ε θ β ε β r ε β θ ε β P T
P = sin θ cos β cos θ cos θ sin β cos θ cos β sin θ sin β sin θ 0 sin β cos β
Klar and Klein [28] derived the geometrical equations in the curvilinear coordinate system. Based on the work performed by Klar and Klein [28], we give the geometrical equations modified for the coordinate in our method as follows.
ε r r = υ r r ε r r = 1 2 r υ r θ + r υ θ r υ θ ε r β = 1 2 ( R β r cos θ ) υ β cos θ + υ r β + ( R β r cos θ ) υ β r ε θ θ = 1 r υ r + υ θ θ ε θ β = 1 2 r ( R β r cos θ ) r υ θ β υ β r sin θ + υ r β + ( R β r cos θ ) υ β θ ε β β = 1 R β r cos θ υ θ sin θ υ r cos θ + υ β β
Equation (34) gives the train energy of the soil. Equation (4), the strain energy of the soil, can be calculated by combining Equation (4) with Equation (34).
(b)
Potential energy of external force
According to the deformation mechanism, the vertical displacement can be given as follows:
w = υ r cos θ cos β + υ θ cos β sin θ υ β sin β
Substituting Equation (35) into Equation (6), the potential energy of gravity can be calculated.
When β = 0, the orthoradial displacement vβ(r, θ, β) is equal to the horizontal displacement on the tunnel face uf. Equation (36) gives the horizontal displacement on the tunnel face.
u f = υ β ( r ,   θ ,   0 )
Combining Equation (36) with Equation (5), the potential energy of support pressure can be calculated.

4.3. Volume Loss Calculation

Vlf is the volume loss at the tunnel face and equals the ratio of the flux on the tunnel face (qsf) to the area of the tunnel cross-section.
V L = q s f / π D 2 4 = 1 2 S m f
VL is the volume loss at ground surface and equals the ratio of the flux on the ground surface (Qs) to the area of the tunnel cross-section, as defined in Klar and Klein [28]. Qs equals qs (π/2). So, the volume loss can be given as follows:
V L = q s f ( π 2 k + 1 ) / π D 2 4 = 1 2 S m f ( π 2 k + 1 )

5. Verification with Numerical Simulation

In this section, we perform numerical simulations by a finite difference computer program FLAC3D to validate the accuracy of the proposed method.
As shown in Figure 6, only a half model is built due to the symmetry along the tunnel central line. The tunnel diameter and cover depth are 12 m. To guarantee the accuracy of the numerical model, the dimension size should be large enough to wake the influence of the boundary, and the specific dimension is shown in Figure 6. The mechanical parameters of the soil and support pressure are shown in Table 1.
In this numerical model, the top surface of the model was free, normal velocity on the four vertical surfaces was restricted, and the boundary surface was fixed. Because we focused on the displacement of soil before the tunnel face, a simplified one-step excavation was applied after reaching initial equilibrium, due to gravity. Based on the assumption that the VLs and VLt were zero, the normal displacement of points on the opening was fixed after excavation. Immediately, the gradient support pressure was applied on the face. We monitored the displacement of the points on the face and surface. The calculation was stopped when the displacements of several points were stable.
We extracted the maximum displacement on the tunnel face and ground surface in the numerical simulations, and calculated the volume loss at the tunnel face and the ground surface. The comparison is shown in Figure 7.
The comparison shows that the proposed method has a good agreement with the numerical simulation.

6. Parameter Analysis

6.1. Influence of Support Pressure Ratio

The influence of the support pressure ratio α on the volume loss with the tunnel diameter ranging from 10 m to 14 m is shown in Figure 8. The tunnel depth was 15 m and the slurry weight was 12 kN/m3. The volume loss decreased linearly with the increase in the support pressure ratio, which implied that the increase in the support pressure ratio was favorable to the control of volume loss, including the volume loss at the tunnel face and the volume loss at ground surface. Moreover, the deeper the tunnel, the bigger the volume loss. The volume loss at the tunnel face is generally bigger than the volume loss at ground surface, due to the soil compression during ground movement. The influence of the tunnel diameter on the volume loss at the ground surface is much smaller than that at the tunnel face.
The influence of the support pressure ratio α on the volume loss with the tunnel depth ranging from 15 m to 21 m is shown in Figure 9. The tunnel diameter was 12 m, and the slurry weight was 12 kN/m3. With the increase in support pressure ratio, volume loss decreased linearly. Moreover, the deeper the tunnel, the faster the volume loss decreased. When α = 1, the tunnel depth did not influence the volume loss at the tunnel face. However, the tunnel depth will influence the soil compression, and the volume loss at the ground surface will be affected by the tunnel depth.

6.2. Influence of Slurry Weight

The influence of the slurry weight γ on the volume loss with the tunnel diameter ranging from 10 m to 14 m and tunnel depth ranging from 15 m to 21 m are shown in Figure 10 and Figure 11, respectively. The support pressure ratio was 0.9. The tunnel depth was 15 m in Figure 10, and tunnel diameter was 12 m in Figure 11. The volume loss decreased linearly with the increase in slurry weight, which implies that the increase in the slurry weight was favorable to the control of volume loss, including the volume loss at the tunnel face and the volume loss at the ground surface. The bigger the tunnel diameter, the quicker the volume loss decreased. However, the tunnel diameter did not influence the decrease rapidly. Moreover, the volume loss increased with the increase in the tunnel diameter and tunnel depth.

7. Conclusions

In this paper, a semi-analytical method based on the principle of minimum total potential energy was proposed to solve the volume loss at the tunnel face and the ground face for shield-driven tunnels. Moreover, a new deformation mechanism considering the soil volume change was designed to describe the ground movement before the tunnel face. Based on the proposed method, the influence of different factors on the volume loss was investigated. The main conclusions are given as follow:
The volume loss decreases with the increase in the support pressure ratio and the slurry weight. The increase in the slurry weight is favorable to the control of volume loss during tunnelling
The volume loss at ground surface is smaller than that at the tunnel face due to the soil compression during ground movement.
The bigger the tunnel diameter, the bigger the volume loss. However, the volume loss at ground surface may not increase with the increase in the tunnel depth because of the soil arching effect.
The deeper the tunnel, the quicker the volume loss decreases with the increase in the support pressure ratio. The bigger the tunnel diameter, the quicker the volume loss decreases with the increase in the slurry weight.
The proposed method provides a more direct way for estimating volume loss at the tunnel face from an energy point of view. Most current researchers have assumed zero volume behind the tunnel face (VLt and VLs), and the considered volume loss is limited to the tunnel face (‘tunnel heading’). A more effective way for estimating all volume loss along the shield advancement may be developed by designing multiple deformation mechanisms according to the field data and numerical simulations.

Author Contributions

G.W.: Conceptualization, formal analysis, software. Q.F.: Funding acquisition, Supervision, Validation. J.D.: Investigation, Writing-review & editing, X.Y.: Visualization. J.W.: Formal analysis, Data curation. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Natural Science Foundation of Beijing, China grant number 8202037, and the Key Project of High-speed Rail Joint Fund of National Natural Science Foundation of China grant number U1934210.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within this article.

Acknowledgments

The authors thank the reviewers for their great help on the article during its review progress.

Conflicts of Interest

The authors declare no conflict of interest.

Notations

Psupport pressure supplied by slurry
Pusupport pressure at the tunnel crown
σearth lateral pressure
σuearth lateral pressure at the tunnel crown
Hdepth of centreline of tunnel
Dtunnel diameter
Ccover of tunnel
VLfvolume loss at the tunnel face
αvolume loss at the tunnel face
γslurry weight
Eelastic modulus of soil
μPoisson’s ratio
Parameter of deformation mechanism
L1distance between the centre of the tunnel face and the point with maximum displacement on tunnel face
kparameter describing extent of volume change of soil
Smfmaximum displacement on the tunnel face
Rfsemi-major axis of maximum displacement line
Rβradius of maximum displacement line
Rmaxradius of circular cross-section
vrradial displacement on circular cross-section
vβorthoradial displacement on circular cross-section
vmmaximum displacement on circular cross-section
qsflux through circular cross-section
Smgmaximum settlement on the ground surface

References

  1. Sagaseta, C. Analysis of undrained soil deformation due to ground loss. Geotechnique 1988, 38, 647–649. [Google Scholar] [CrossRef]
  2. Verruijt, A. A Complex Variable Solution for a Deforming Circular Tunnel in an Elastic Half-Plane. Int. J. Numer. Anal. Methods Geomech. 1997, 21, 77–89. [Google Scholar] [CrossRef]
  3. Loganathan, N.; Poulos, H.G. Analytical prediction for tunneling-induced ground movements in clays. J. Geotech. Geoenviron. Eng. 1998, 124, 846–856. [Google Scholar] [CrossRef]
  4. Bobet, A. Analytical Solutions for Shallow Tunnels in Saturated Ground. J. Eng. Mech. 2001, 127, 1258–1266. [Google Scholar] [CrossRef]
  5. Park, K.H. Analytical solution for tunnelling-induced ground movement in clays. Tunn. Undergr. Space Technol. 2005, 20, 249–261. [Google Scholar] [CrossRef]
  6. Li, S.; Li, P.; Zhang, M. Analysis of additional stress for a curved shield tunnel. Tunn. Undergr. Space Technol. 2021, 107, 103675. [Google Scholar] [CrossRef]
  7. Fang, Q.; Wang, G.; Yu, F.C.; Du, J.M. Analytical algorithm for longitudinal deformation profile of a deep tunnel. J. Rock Mech. Geotech. Eng. 2021, 13, 845–854. [Google Scholar] [CrossRef]
  8. Li, P.F.; Wei, Y.J.; Zhang, M.J.; Huang, Q.F.; Wang, F. Influence of non-associated flow rule on passive face instability for shallow shield tunnels. Tunn. Undergr. Space Technol. 2022, 119, 104202. [Google Scholar] [CrossRef]
  9. Fang, Q.; Zhang, D.; Wong, L.N.Y. Shallow tunnelling method (STM) for subway station construction in soft ground. Tunn. Undergr. Space Technol. 2012, 29, 10–30. [Google Scholar] [CrossRef]
  10. Vu, M.N.; Broere, W.; Bosch, J. Volume loss in shallow tunnelling. Tunn. Undergr. Space Technol. 2016, 59, 77–90. [Google Scholar] [CrossRef] [Green Version]
  11. Fang, Q.; Du, J.M.; Li, J.Y.; Zhang, D.L.; Cao, L.Q. Settlement characteristics of large-diameter shield excavation below existing subway in close vicinity. J. Cent. South Univ. 2021, 28, 882–897. [Google Scholar] [CrossRef]
  12. Gui, M.W.; Chen, S.L. Estimation of transverse ground surface settlement induced by DOT shield tunneling. Tunn. Undergr. Space Technol. 2013, 33, 119–130. [Google Scholar] [CrossRef]
  13. Shahin, H.M.; Nakai, T.; Ishii, K.; Iwata, T.; Kuroi, S. Investigation of influence of tunneling on existing building and tunnel: Model tests and numerical simulations. Acta Geotech. 2016, 11, 679–692. [Google Scholar] [CrossRef]
  14. Zheng, H.; Li, P.; Ma, G. Stability analysis of the middle soil pillar for asymmetric parallel tunnels by using model testing and numerical simulations. Tunn. Undergr. Space Technol. 2020, 108, 103686. [Google Scholar] [CrossRef]
  15. Bezuijen, A.; Talmon, A.M. Processes around a TBM. In Geotechnical Aspects of Underground Construction in Soft Ground—Proceedings of the 6th International Symposium, IS-SHANGHAI 2008; CRC Press: Boca Raton, FL, USA, 2009; pp. 3–13. [Google Scholar]
  16. Verruijt, A.; Booker, J.R. Surface settlements due to deformation of a tunnel in an elastic half plane. Geotechnique 1996, 46, 753–756. [Google Scholar] [CrossRef]
  17. González, C.; Sagaseta, C. Patterns of soil deformations around tunnels. Application to the extension of Madrid Metro. Comput. Geotech. 2001, 28, 445–468. [Google Scholar] [CrossRef]
  18. Lee, K.M.; Rowe, R.K.; Lo, K.Y. Subsidence owing to tunnelling. I. Estimating the gap parameter. Can. Geotech. J. 1992, 29, 929–940. [Google Scholar] [CrossRef]
  19. Rowe, R.K.; Lee, K.M. Subsidence owing to tunnelling. II. Evaluation of a prediction technique. Can. Geotech. J. 1992, 29, 941–954. [Google Scholar] [CrossRef]
  20. Cao, L.; Zhang, D.; Fang, Q. Semi-analytical prediction for tunnelling-induced ground movements in multi-layered clayey soils. Tunn. Undergr. Space Technol. 2020, 102, 103446. [Google Scholar] [CrossRef]
  21. Mair, R.J.; Gunn, M.J.; O’Reilly, M.P. Ground movement around shallow tunnels in soft clay. Tunn. Tunn. Int. 1982, 14, 45–48. [Google Scholar]
  22. Attewell, P.B.; Yeates, J.; Selby, A.R. Soil Movements Induced by Tunnelling and their Effects on Pipelines and Structures; Methuen, Inc.: New York, NY, USA, 1986. [Google Scholar]
  23. Macklin, S.R. The prediction of volume loss due to tunnelling in over consolidated clay based on heading geometry and stability number. Ground Eng. 1999, 32, 30–33. [Google Scholar]
  24. Dimmock, P.S.; Mair, R.J. Estimating volume loss for open-face tunnels in London Clay. Proc. Inst. Civ. Eng. Geotech. Eng. 2007, 160, 13–22. [Google Scholar] [CrossRef]
  25. Osman, A.S. Predicting 2D ground movements around tunnels in undrained clay. Geotechnique 2006, 56, 597–604. [Google Scholar] [CrossRef] [Green Version]
  26. Liu, X.; Fang, Q.; Zhang, D.; Liu, Y. Energy-based prediction of volume loss ratio and plastic zone dimension of shallow tunnelling. Comput. Geotech. 2020, 118, 103343. [Google Scholar] [CrossRef]
  27. Palmer, A.C.; Mair, R.J. Ground movements above tunnels: A method for calculating volume loss. Can. Geotech. J. 2011, 48, 451–457. [Google Scholar] [CrossRef]
  28. Klar, A.; Klein, B. Energy-based volume loss prediction for tunnel face advancement in clays. Geotechnique 2014, 64, 776–786. [Google Scholar] [CrossRef]
  29. Mollon, G.; Dias, D.; Soubra, A.-H. Continuous velocity fields for collapse and blowout of a pressurized tunnel face in purely cohesive soil. Int. J. Numer. Anal. Methods Geomech. 2012, 30, 1303–1336. [Google Scholar] [CrossRef] [Green Version]
  30. Attewell, P.B.; Farmer, I.W. Ground disturbance caused by shield tunnel over consolidated clay. Eng. Geol. 1974, 8, 361–381. [Google Scholar] [CrossRef]
  31. Deane, A.P.; Bassett, R.H. The Heathrow Express trial tunnel. Proc. Inst. Civ. Eng. Geotech. Eng. 1995, 113, 144–156. [Google Scholar] [CrossRef]
  32. Li, P.F.; Chen, K.; Wang, F.; Li, Z. An upper-bound analytical model of blow-out for a shallow tunnel in sand considering the partial failure within the face. Tunn. Undergr. Space Technol. 2019, 91, 102989. [Google Scholar] [CrossRef]
  33. Yu, L.; Zhang, D.; Fang, Q.; Cao, L.; Zhang, Y.; Xu, T. Face stability of shallow tunnelling in sandy soil considering unsupported length. Tunn. Undergr. Space Technol. 2020, 102, 103445. [Google Scholar] [CrossRef]
Figure 1. Sketch of proposed problem.
Figure 1. Sketch of proposed problem.
Applsci 12 01794 g001
Figure 2. Schematic diagram of the proposed model.
Figure 2. Schematic diagram of the proposed model.
Applsci 12 01794 g002
Figure 3. Three-dimensional schematic diagram of modified deformation mechanism.
Figure 3. Three-dimensional schematic diagram of modified deformation mechanism.
Applsci 12 01794 g003
Figure 4. Side view of the modified deformation mechanism.
Figure 4. Side view of the modified deformation mechanism.
Applsci 12 01794 g004
Figure 5. Circular cross-section Πβ.
Figure 5. Circular cross-section Πβ.
Applsci 12 01794 g005
Figure 6. Numerical model.
Figure 6. Numerical model.
Applsci 12 01794 g006
Figure 7. Comparison between numerical simulation and proposed method. (a) Volume loss at tunnel face Vlf, (b) volume loss at ground surface Vlg.
Figure 7. Comparison between numerical simulation and proposed method. (a) Volume loss at tunnel face Vlf, (b) volume loss at ground surface Vlg.
Applsci 12 01794 g007
Figure 8. Volume loss versus support pressure ratio α with different tunnel diameter D. (a) volume loss at tunnel face Vlf, (b) volume loss at ground surface Vlg.
Figure 8. Volume loss versus support pressure ratio α with different tunnel diameter D. (a) volume loss at tunnel face Vlf, (b) volume loss at ground surface Vlg.
Applsci 12 01794 g008
Figure 9. Volume loss versus support pressure ratio α with different tunnel depth H. (a) Volume loss at tunnel face Vlf, (b) volume loss at ground surface Vlg.
Figure 9. Volume loss versus support pressure ratio α with different tunnel depth H. (a) Volume loss at tunnel face Vlf, (b) volume loss at ground surface Vlg.
Applsci 12 01794 g009
Figure 10. Volume loss versus slurry weight γ with different tunnel diameter D. (a) Volume loss at tunnel face Vlf, (b) volume loss at ground surface Vlg.
Figure 10. Volume loss versus slurry weight γ with different tunnel diameter D. (a) Volume loss at tunnel face Vlf, (b) volume loss at ground surface Vlg.
Applsci 12 01794 g010
Figure 11. Volume loss versus slurry weight γ with different tunnel depth H. (a) Volume loss at tunnel face Vlf, (b) volume loss at ground surface Vlg.
Figure 11. Volume loss versus slurry weight γ with different tunnel depth H. (a) Volume loss at tunnel face Vlf, (b) volume loss at ground surface Vlg.
Applsci 12 01794 g011
Table 1. Mechanical parameters of soil and support pressure.
Table 1. Mechanical parameters of soil and support pressure.
ParametersValue
Buoyant unit weight of soil (kN/m3)18
Unit weight of water (kN/m3)10
Elastic modulus (MPa)1
Poisson’s ratio0.3
Support pressure ratio0.85, 0.9, 0.95, 1.00
Slurry weight (kN/m3)12, 13, 14
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, G.; Fang, Q.; Du, J.; Yang, X.; Wang, J. Estimating Volume Loss for Shield-Driven Tunnels Based on the Principle of Minimum Total Potential Energy. Appl. Sci. 2022, 12, 1794. https://doi.org/10.3390/app12041794

AMA Style

Wang G, Fang Q, Du J, Yang X, Wang J. Estimating Volume Loss for Shield-Driven Tunnels Based on the Principle of Minimum Total Potential Energy. Applied Sciences. 2022; 12(4):1794. https://doi.org/10.3390/app12041794

Chicago/Turabian Style

Wang, Gan, Qian Fang, Jianming Du, Xiaoxu Yang, and Jun Wang. 2022. "Estimating Volume Loss for Shield-Driven Tunnels Based on the Principle of Minimum Total Potential Energy" Applied Sciences 12, no. 4: 1794. https://doi.org/10.3390/app12041794

APA Style

Wang, G., Fang, Q., Du, J., Yang, X., & Wang, J. (2022). Estimating Volume Loss for Shield-Driven Tunnels Based on the Principle of Minimum Total Potential Energy. Applied Sciences, 12(4), 1794. https://doi.org/10.3390/app12041794

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop