Next Article in Journal
Inference about the Ratio of the Coefficients of Variation of Two Independent Symmetric or Asymmetric Populations
Previous Article in Journal
The Computer Course Correlation between Learning Satisfaction and Learning Effectiveness of Vocational College in Taiwan
 
 
Erratum published on 7 February 2020, see Symmetry 2020, 12(2), 258.
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analytical Solution of Ground Stress Induced by Shallow Tunneling with Arbitrary Distributed Loads on Ground Surface

1
Beijing Jiaotong University, Beijing 100044, China
2
Beijing Urban Engineering Design & Research Institute Co., Ltd., Beijing 100044, China
3
College of Civil and Transportation Engineering, Shenzhen University, Shenzhen 518060, China
*
Author to whom correspondence should be addressed.
Symmetry 2019, 11(6), 823; https://doi.org/10.3390/sym11060823
Submission received: 22 April 2019 / Revised: 5 June 2019 / Accepted: 13 June 2019 / Published: 21 June 2019

Abstract

:
The impact of shallow tunnel construction on the surrounding environment is often considered as a symmetric half-plane problem with circular holes. In this research, the analytical solutions of the ground stresses and displacements of a shallow circular tunnel in an elastic half-plane under arbitrary distributed loads on ground surface were derived, based on the complex variable method. Then, an application was implemented to analyze the potential plastic zone induced by shallow tunneling adjacent to the ground surface structures. The verification of the results obtained from the proposed analytical prediction model was carried out using the numerical simulations. Additionally, the influences of different boundary condition (different magnitudes and ranges of arbitrary distributed loads and different symmetric boundary conditions of the tunnel perimeter) on the distribution characteristics of the potential plastic zones were analyzed. In general, the results showed that the larger the pile loads and the closer the relative position between the tunnel and distributed loads, the more distinct the coalesced trends of the potential plastic zones around the tunnel and the potential plastic zones around the distributed loads.

1. Introduction

As people pay increasingly more attention to environmental protection and ground transportation becomes busier, the underground excavation method is replacing the cut and cover method to become the most popular construction method [1]. The tunnel construction method in urban areas is divided into two types, namely, the closed face tunneling method (shield method) and the open type face tunneling methods (shallow tunneling method and analysis of controlled deformation of in rocks and soils (ADECO-RS) approach). Many empirical and analytical formulas have been developed to calculate the ground displacements and stresses induced by shallow tunneling for the green-field case, which includes Peck’s empirical formula [2,3], the virtual image technique [4,5,6], the complex variable method [7,8], the general series form stress function in polar coordinates [9,10,11,12], and the stochastic medium theory [13]. Recently, the superposition methods were adopted to calculate the ground displacements and stresses induced by the combination of shallow tunneling and pile foundation loads [14,15].
In this research, the solutions of the ground stresses and displacements of a shallow circular tunnel considering the influence of arbitrary distributed loads on ground surface were derived using the complex variable method. Then, an application was implemented to analyze the potential plastic zone induced by shallow tunneling and loads of the ground surface structures. The results obtained from the proposed analytical prediction model were compared with the results obtained from the numerical simulations to verify the correctness of the theoretical formula. The influences of different boundary condition (different magnitudes and ranges of the arbitrary distributed loads; different symmetric boundary conditions of the tunnel perimeter) on the distribution characteristics of the potential plastic zones were also analyzed.

2. Problem Description

In this study, a symmetric elastic half-plane (y < 0) involving a circular tunnel of radius r at a depth h below the ground surface was considered. Tunneling will produce ground loss, which will cause a change of the deformation and stress of the surrounding strata. The boundary of a shallow circular tunnel undergoes a given distribution of displacements, i.e., uniform radial displacement and ovalization. Moreover, the arbitrary distributed loads are also considered to apply to ground surface to represent the effect of adjacent ground surface structures. The shear modulus and the Poisson’s ratio of the ground are denoted by μ and ν. The geometry of the problem is shown in Figure 1.
The proposed problem was solved in the two following sections: (i) the first partial solution comprises the general formulas of the exact analytical solutions of the ground stresses and the displacements in a symmetric elastic half-plane based on the complex variable method. The solution given in this paper is a generalization of Verruijt’s solution [7,8], which not only considers the ground loss but also includes the effect of arbitrary distributed loads on ground surface (see Section 3). (ii) The second partial solution comprises the derivation and the summary of the recursive relations of the Laurent series coefficients used in the complex variable method obtained from the symmetric boundary conditions of a tunnel and the arbitrary distributed loads on ground surface (see Section 4).

3. General Formulas of the Exact Analytical Solutions of Ground Displacements and Stresses

Figure 2 shows a circular tunnel located in the lower part of a symmetric elastic plane. It was assumed that the outer boundary of the ground surface suffered arbitrary distributed loads and the inner boundary of the tunnel profile underwent a specific symmetric distribution of convergence deformation. According to the complex variable method, the original domain in the z-plane (physical plane) was mapped conformally onto an annular region on the ζ-plane (mapped plane) with the following conformal transformation:
z = ω ( ζ ) = i h 1 α 2 1 + α 2 1 + ζ 1 ζ = i a 1 + ζ 1 ζ ,
where α was obtained by
α = h r ( h r ) 2 1 .
In the z-plane, the solutions of the ground stresses and displacements due to the tunnel construction were donated by two analytic functions, φ1 (z) and ψ1 (z), and they were obtained as follows:
σ x x + σ y y = 2 [ φ 1 ( z ) + φ 1 ( z ) ¯ ] = 4 Re [ φ 1 ( z ) ] ,
σ y y σ x x + 2 i σ x y = 2 [ z ¯ φ 1 ( z ) + ψ 1 ( z ) ] ,
2 μ ( u x + i u y ) = κ φ 1 ( z ) z φ 1 ( z ) ¯ ψ 1 ( z ) ¯ ,
where μ, κ, and i are the shear modulus of the ground material, a related parameter for Poisson’s ratio ν, and the imaginary constant, respectively.
By virtue of the conformal transformation function, the two analytic functions, φ1(z) and ψ1(z), could be transformed to the functions in terms of ζ as follows:
φ 1 ( z ) = φ 1 ( ω ( ζ ) ) = φ ( ζ ) ,
ψ 1 ( z ) = ψ 1 ( ω ( ζ ) ) = ψ ( ζ ) .
Moreover, φ(ζ) and ψ(ζ) would be expanded in a Laurent series in the ζ-space in the following equations due to the characteristics of the analytical function.
φ ( ζ ) = a 0 + k = 1 a k ζ k + k = 1 b k ζ k ,
ψ ( ζ ) = c 0 + k = 1 c k ζ k + k = 1 d k ζ k ,
where the Laurent series coefficients a0, ak, bk, c0, ck, and dk could be calculated by the recursive relations obtained from the boundary conditions.

4. Recursive Relations and Laurent Series Coefficients for the Proposed Model

The solution of displacements and stresses obtained above is a generalization of Verruijt’s solution, which not only considers the ground loss but also includes the effect of arbitrary distributed loads on the ground surface. In this section, the recursive relations from specific boundary conditions are derived, and the Laurent series coefficients are summarized for the boundary deformations of a tunnel and the arbitrary distributed loads on ground surface.

4.1. The Outer Boundary of the Ground Surface

The outer boundary of the ground surface y = 0 underwent arbitrary distributed loads. This gave the following equation:
y = 0 : φ 1 ( z ) + z φ 1 ( z ) ¯ + ψ 1 ( z ) ¯ = i l 0 z ( q x + i q y ) d z = Q 1 ( z ) + C .
This equation can be transformed to the ζ-plane and expressed as the form
| ζ 0 | = 1 : φ ( ζ 0 ) + ω ( ζ 0 ) ω ( ζ 0 ) ¯ φ ( ζ 0 ) ¯ + ψ ( ζ 0 ) ¯ = Q ( ζ 0 ) + C = K = + e k σ k + C .
Then the boundary Condition (7) can be elaborated in terms of Laurent series coefficients as follows:
k = 1 a k σ k + k = 1 b k σ k + 1 2 k = 1 ( k + 1 ) a ¯ k + 1 σ k 1 2 k = 2 ( k 1 ) b ¯ k 1 σ k 1 2 k = 2 ( k 1 ) a ¯ k 1 σ k + 1 2 k = 1 ( k + 1 ) b ¯ k + 1 σ k + a 0 + 1 2 a ¯ 1 + 1 2 b ¯ 1 + c ¯ 0 + k = 1 c ¯ k σ k + k = 1 d ¯ k σ k = K = + e k σ k + C
Therefore, the coefficients c0, ck, and dk satisfy the following recursive relations:
{ c 0 = e 0 ¯ a ¯ 0 1 2 a 1 1 2 b 1 + C ¯ c k = e k b ¯ k + 1 2 ( k 1 ) a k 1 1 2 ( k + 1 ) a k + 1 d k = e ¯ k a ¯ k + 1 2 ( k 1 ) b k 1 1 2 ( k + 1 ) b k + 1 ,
where the related coefficients ek are in Fourier series terms in Equation (9), as seen in Appendix A.

4.2. The Inner Boundary of the Tunnel Profile

In general, the inner boundary of a tunnel profile will undergo convergence deformation. Researchers have proposed a variety of the specific symmetric distributions of the convergence deformation for different soil types or engineering conditions. Verruijt and Booker [5] stated that a shallow tunnel profile undergoes a combination deformation of uniform convergence and ovalization. In this research, considering the buoyancy effect, the vertical translation was incorporated into the convergence deformation of a shallow tunnel [11]. Figure 3 shows the four possible forms of boundary conditions (B.C. 1, B.C.2, B.C.3, and B.C.4) that are more in accord with engineering practice.
The boundary condition for a shallow tunnel profile in the ζ-plane is
| ζ | = α : 2 μ ( u x + i u y ) = κ φ ( ζ 0 ) ω ( ζ 0 ) ω ( ζ 0 ) ¯ φ ( ζ 0 ) ¯ ψ ( ζ 0 ) ¯ = G ( ζ 0 ) ,
( 1 α σ ) [ κ φ ( ζ 0 ) ω ( ζ 0 ) ω ( ζ 0 ) ¯ φ ( ζ 0 ) ¯ ψ ( ζ 0 ) ¯ ] = ( 1 α σ ) G ( ζ 0 ) = G ( ζ 0 ) = G ( α σ ) = k = A k σ k .
After some elaborations, the symmetric boundary Condition (10) now gives
G ( α σ ) = k = A k σ k =    ( κ + 1 ) a 0 + ( 1 α 2 ) a ¯ 1 ( κ + α 2 ) b 1 + ( e 0 + α 2 e ¯ 1 C ) + [ ( κ + 1 ) α 2 a 0 + ( κ α 2 + 1 ) a 1 + ( 1 α 2 ) b ¯ 1 ( α 2 e 0 + e 1 α 2 C ) ] α 1 σ + k = 1 [ ( 1 α 2 ) ( k + 1 ) a ¯ k + 1 ( κ α 2 k + α 2 ) b k + 1 ( 1 α 2 ) k a ¯ k + ( κ α 2 k + 1 ) b k ( e ¯ k α 2 e ¯ k + 1 ) ] α k σ k + k = 2 [ ( κ α 2 k + 1 ) a k + ( 1 α 2 ) k b ¯ k ( κ α 2 k + α 2 ) a k 1 ( 1 α 2 ) ( k 1 ) b ¯ k 1 ( e k α 2 e k 1 ) ] α k σ k
Therefore, the coefficients a0, ak, and bk are calculated using the following equations integrating the related coefficients Ak in Laurent series terms
{ ( 1 α 2 ) a ¯ 1 ( κ + α 2 ) b 1 = A 0 ( κ + 1 ) a 0 + ( e 0 α 2 e ¯ 1 + C ) ( κ α 2 + 1 ) a ¯ 1 + ( 1 α 2 ) b 1 = A ¯ 1 α + ( κ + 1 ) α 2 a ¯ 0 + ( α 2 e ¯ 0 + e ¯ 1 α 2 C ¯ ) ( 1 α 2 ) ( k + 1 ) a ¯ k + 1 ( κ α 2 k + α 2 ) b k + 1 = ( 1 α 2 ) k a ¯ k ( κ α 2 k + 1 ) b k + ( e ¯ k α 2 e ¯ k + 1 ) + A k α k ( κ α 2 k + 2 + 1 ) a ¯ k + 1 + ( 1 α 2 ) ( k + 1 ) b k + 1 = α 2 ( κ α 2 k + 1 ) a ¯ k + ( 1 α 2 ) k b k + ( e ¯ k + 1 α 2 e ¯ k ) + A ¯ k + 1 α k + 1 .
{ ( 1 α 2 ) a ¯ 1 ( κ + α 2 ) b 1 = A 0 ( κ + 1 ) a 0 + ( e 0 α 2 e ¯ 1 + C ) ( κ α 2 + 1 ) a ¯ 1 + ( 1 α 2 ) b 1 = A ¯ 1 α + ( κ + 1 ) α 2 a ¯ 0 + ( α 2 e ¯ 0 + e ¯ 1 α 2 C ¯ ) ,
where the Ak variables are in Laurent series terms in Equations (13) and (14), seen in Appendix B.

5. An Application of the Derived Exact Analytic Solutions

5.1. Prediction of the Distribution Characteristics of the Potential Plastic Zone

It is often unavoidable that urban subway tunnels are constructed adjacent to ground surface structures, e.g., existing buildings. It is important to predict the degree of interactions between the tunnel construction and the ground surface structures. Generally, the degree of the interactions will be evaluated with a prediction of the distribution characteristics of the potential plastic zone. In this paper, the practical problem of the construction of shallow subway tunnels adjacent to surface building is simplified to the mechanics model, as shown in Figure 4.
By combining the exact analytical solutions of the ground stresses obtained from this research and the Mohr–Coulomb yield criterion, an equation for the boundary of the potential plastic zone due to shallow tunneling adjacent to the surface building was obtained. Based on the approach presented above, the boundaries of the potential plastic zones induced by tunneling adjacent to the surface building were drawn using MATLAB software. The assumed parameters in the presented calculations were silt clayey soil with a cohesion c = 30 kPa and an angle of internal friction φ = 30°. The radius of the circular tunnel was 3 m, and the center depth h was 10 m. The uniform convergence u0 was 30 mm. The magnitudes of the distributed loads q were 100, 200, and 300 kN/m, and the ranges of the distributed loads (a0, b0) were (0 m, 7 m), (2 m, 9 m), and (4 m, 11 m).

5.2. Influences of Different Magnitudes and Ranges of Arbitrary Distributed Loads on the Potential Plastic Zones

The influences of the different magnitudes and ranges of arbitrary distributed loads on the distribution characteristics of the potential plastic zones are analyzed in Figure 5 and Figure 6.
The influences of three different magnitudes of the distributed loads q = 100, 200, and 300 kN/m on the distribution characteristics of the potential plastic zones are plotted in Figure 5. The results show that the potential plastic zones around the tunnel and around the ground surface structure coalesced when the magnitude of distributed load q was big enough, while if the magnitude of the distributed load q was relatively small, the potential plastic zones separated into two parts. Moreover, the plastic zone without distributed load obtained from the solution of Verruijt [7,8] is also added in Figure 5 to allow for interpretation of the superiority of the analytical solution proposed in this paper.
Figure 6 presents the influences of three different ranges of distributed loads (a0, b0) = (0 m, 7 m), (2 m, 9 m), (4 m, 11 m) on the distribution characteristics of the potential plastic zone induced by tunneling. It can be observed that when the range of distributed loads (a0, b0) was close enough to the tunnel, the tunneling-induced potential plastic zones around the tunnel and the potential plastic zones around the distributed loads coalesced.
In general, the results show that the larger the pile loads and the closer the relative position between the tunnel and the distributed loads, the more distinct the coalesced trends of the potential plastic zones around the tunnel and the potential plastic zones around the distributed loads.

5.3. Influences of Different Tunnel Boundary Conditions on the Potential Plastic Zones

The influences of different tunnel boundary conditions (B.C.1, B.C.2, B.C.3, and B.C.4) on the distribution characteristics of the potential plastic zones are shown in Figure 7. The results indicated that different tunnel boundary conditions greatly affected the distribution characteristics of the potential plastic zones. More attention should be paid to selecting tunnel boundary conditions because the ovalization and the vertical translation significantly affect the shape of the potential plastic zones.

5.4. Verification of the Potential Plastic Zone Obtained from the Derived Analytical Solutions with the Numerical Simulations

Numerical simulations were performed by using the professional software OptumG2 [16] based on identical parameters. Figure 8a,b indicate the potential plastic zone for the green-field and for the condition with arbitrary distributed loads on ground surface obtained from the numerical simulations. Figure 8c presents the comparison of the potential plastic zone obtained from numerical simulations and analytical solution, which shows that the range of potential plastic zone are similar to each other. It is worth pointing out that there was still a difference between the potential plastic zones obtained from the derived analytical solutions with the numerical simulations. The boundary contour of the potential plastic zone obtained from the derived analytical solutions was obtained by combining the yield criterion with the analytical solution of the elastic hypothesis, while the numerical simulation could directly calculate the stress in the elastic and plastic zones and obtain the potential plastic zone.

6. Conclusions

In this research, the solutions of the ground stresses and displacements of a shallow circular tunnel in a symmetric elastic half-plane under arbitrary distributed loads on ground surface were derived using the complex variable method. Then, an application was implemented to analyze the potential plastic zone caused by shallow tunneling adjacent to the ground surface structures. The influences of different boundary condition (different magnitudes and ranges of the arbitrary distributed loads, different symmetric boundary conditions of the tunnel perimeter) on the distribution characteristics of the potential plastic zones were analyzed. The main conclusions were as follows.
(1)
The potential plastic zones around the tunnel and around the ground surface structure coalesced when the magnitude of the distributed load q was big enough, while if the magnitude of the distributed load q was relatively small, the potential plastic zones separated into two parts.
(2)
If the range of the distributed loads (a0, b0) was close enough to the tunnel, the tunneling-induced potential plastic zones around the tunnel and the potential plastic zones around the distributed loads coalesced.
(3)
The results indicated that different symmetric tunnel boundary conditions greatly affected the distribution characteristics of the potential plastic zone.

Author Contributions

J.W. and K.H. conceived of the idea of using complex variable method to solve the the analytical solutions of the ground stresses and displacements of a shallow circular tunnel in an elastic half-plane under arbitrary distributed loads on ground surface and linked solutions to applications. Z.L. completed most of the details of the calculations.

Funding

The authors acknowledge the financial support provided by the Fundamental Research Funds for the Central Universities of China (Grant no. 2015YJS128).

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

rradius of a circular tunnel
hburial depth of the tunnel below the ground surface
μshear modulus
νPoisson’s ratio
zphysical plane
x, ycoordinate axis in physical plane
ζmapped plane
ξ, ηcoordinate axis in mapped plane
l1, l2, d1, d2range of arbitrary distributed loads on ground surface
qx, qymagnitude of arbitrary distributed loads on ground surface
σxx, σyy, σxycomponents of stress in physical plane
ux, uycomponents of displacement in physical plane
ω(ζ)conformal transformation
αa parameter defines the conformal transformation
φ1(z), ψ1 (z)analytic functions in physical plane
φ(ζ), ψ(ζ)analytic functions in mapped plane
iimaginary constant
a0, ak, bk, c0, ck, dkLaurent series coefficients related to boundary conditions
κrelated parameter to Poisson’s ratio ν

Appendix A. Fourier Coefficients for Arbitrary Distributed Loads

Fourier coefficients for arbitrary distributed loads on ground surface are as follows:
e 0 = q y + i q x 2 π [ σ 1 ( 1 1 2 σ 1 ) σ 2 ( 1 1 2 σ 2 ) ]
e k = q y + i q x 2 k π [ ( σ 1 σ 2 ) sin k + cos ( k σ 1 ) cos ( k σ 2 ) k ] + i [ ( σ 1 σ 2 ) cos k + sin ( k σ 1 ) sin ( k σ 2 ) k ]
e k = e ¯ k

Appendix B. Fourier Coefficients for the Four Different Symmetric Boundary Conditions

Fourier coefficients for the four different symmetric boundary conditions proposed by Park (2004) [11] had been derived by Wang and Li (2009) [17] as follows:
{ A k = 0      k < 0 A 0 = 2 i μ u 0 α A 1 = 2 i μ u 0 A k = 0      k > 1   ( B . C . 1 )
{ A k = 0      k < 0 A 0 = ( 1 + α ) 2 i μ u 0 A 1 = ( 2 + 3 α α 3 ) i μ u 0 A k = ( 1 α 2 ) 2 α k 2 i μ u 0      k > 1   ( B . C . 2 )
{ A k = ( α 2 1 ) 2 4 α k 1 i μ u 0      k < 0 A 0 = ( 1 + α ) 2 i μ u 0 A 1 = ( 3 2 + 3 α + α 2 α 3 1 2 α 4 ) i μ u 0 A k = ( α 2 1 ) 2 ( 4 α + 3 ) + ( α 2 1 ) 3 ( k + 1 ) 4 α k 3 i μ u 0      k > 1   ( B . C . 3 )
{ A k = 3 ( α 2 1 ) 2 8 α k 1 i μ u 0      k < 0 A 0 = ( 3 4 + 5 2 α + 3 4 α 2 ) i μ u 0 A 1 = ( 7 4 + 9 4 α + 3 2 α 2 3 4 α 3 3 4 α 4 ) i μ u 0 A k = 3 ( α 2 1 ) 2 ( k α 2 + α 2 + 2 α + 2 k ) 8 α k 3 i μ u 0      k > 1   ( B . C . 4 )

References

  1. Pinto, F.; Whittle, A.J. Ground movements due to shallow tunnels in soft ground. I: Analytical solutions. J. Geotech. Geoenviron. 2014, 140, 04013040. [Google Scholar] [CrossRef]
  2. Peck, R.B. Deep excavations and tunneling in soft ground. In Proceedings of the 7th International Conference on Soil Mechanics and Foundation Engineering, Mexico City, Mexico, 25–29 August 1969; pp. 225–290. [Google Scholar]
  3. Mair, R.J.; Taylor, R.N. Theme lecture: Bored tunneling in the urban environment. In Proceedings of the 14th International Conference on Soil Mechanics and Foundation Engineering, Hamburg, Germany, 6–12 September 1999; pp. 2353–2385. [Google Scholar]
  4. Sagaseta, C. Analysis of undrained soil deformation due to ground loss. Geotechnique 1987, 37, 301–320. [Google Scholar] [CrossRef]
  5. Verruijt, A.; Booker, J.R. Surface settlements due to deformation of a tunnel in an elastic half plane. Geotechnique 1996, 46, 751–756. [Google Scholar] [CrossRef]
  6. 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]
  7. Verruijt, A. A complex variable solution for a deforming circular tunnel in an elastic half plane. Int. J. Numer. Anal. Met. 1997, 21, 77–89. [Google Scholar] [CrossRef]
  8. Verruijt, A. Deformations of an elastic half plane with a circular cavity. Int. J. Solids Struct. 1998, 35, 2795–2804. [Google Scholar] [CrossRef]
  9. Bobet, A. Analytical solutions for shallow tunnels in saturated ground. ASCE J. Eng. Mech. 2001, 127, 1258–1266. [Google Scholar] [CrossRef]
  10. Chou, W.I.; Bobet, A. Prediction of ground deformations in shallow tunnels in clay. Tunn. Undergr. Space Technol. 2002, 17, 3–19. [Google Scholar] [CrossRef]
  11. Park, K.H. Elastic solution for tunneling-induced Ground Movements in Clays. Int. J. Geomech. 2004, 4, 310–318. [Google Scholar] [CrossRef]
  12. Park, K.H. Analytical solution for tunneling-induced ground movement in clays. Tunn. Undergr. Space Technol. 2005, 20, 249–261. [Google Scholar] [CrossRef]
  13. Yang, J.S.; Liu, B.C.; Wang, M.C. Modeling of tunneling-induced ground surface movements using stochastic medium theory. Tunn. Undergr. Space Technol. 2004, 19, 113–123. [Google Scholar] [CrossRef]
  14. Xiang, Y.Y.; Feng, S.Q. Theoretical prediction of the potential plastic zone of shallow tunneling adjacent to pile foundation in soils. Tunn. Undergr. Space Technol. 2013, 38, 115–121. [Google Scholar] [CrossRef]
  15. Guo, C.X.; Han, K.H.; Kong, H.; Shi, L.L. Explicit form of exact analytical solution for calculating ground displacement and stress induced by shallow tunneling and its application. Adv. Civ. Eng. 2019, 2019, 5739123. [Google Scholar] [CrossRef]
  16. Optum Computational Engineering (OptumCE). Finite-Element Package OptumG2; OptumCE: Copenhagen, Denmark, 2015. [Google Scholar]
  17. Wang, L.Z.; Li, L.L.; Lv, X.J. Complex variable solutions for tunneling-induced ground movement. Int. J. Geomech. 2009, 9, 63–72. [Google Scholar] [CrossRef]
Figure 1. Shallow tunnel under arbitrary distributed loads on ground surface.
Figure 1. Shallow tunnel under arbitrary distributed loads on ground surface.
Symmetry 11 00823 g001
Figure 2. Conformal transformation for a single tunnel.
Figure 2. Conformal transformation for a single tunnel.
Symmetry 11 00823 g002
Figure 3. Summarized boundary conditions of shallow tunnels [11].
Figure 3. Summarized boundary conditions of shallow tunnels [11].
Symmetry 11 00823 g003
Figure 4. Conceptual model of the project practical problem.
Figure 4. Conceptual model of the project practical problem.
Symmetry 11 00823 g004
Figure 5. Boundaries of potential plastic zones for different magnitude of distributed loads q.
Figure 5. Boundaries of potential plastic zones for different magnitude of distributed loads q.
Symmetry 11 00823 g005
Figure 6. Boundaries of potential plastic zones for different range of distributed loads (a0, b0).
Figure 6. Boundaries of potential plastic zones for different range of distributed loads (a0, b0).
Symmetry 11 00823 g006
Figure 7. Boundaries of potential plastic zones for different boundary conditions of shallow tunnel.
Figure 7. Boundaries of potential plastic zones for different boundary conditions of shallow tunnel.
Symmetry 11 00823 g007
Figure 8. Verification of the potential plastic zone obtained from the derived analytical solutions with the numerical simulations.
Figure 8. Verification of the potential plastic zone obtained from the derived analytical solutions with the numerical simulations.
Symmetry 11 00823 g008

Share and Cite

MDPI and ACS Style

Li, Z.; Wang, J.; Han, K. Analytical Solution of Ground Stress Induced by Shallow Tunneling with Arbitrary Distributed Loads on Ground Surface. Symmetry 2019, 11, 823. https://doi.org/10.3390/sym11060823

AMA Style

Li Z, Wang J, Han K. Analytical Solution of Ground Stress Induced by Shallow Tunneling with Arbitrary Distributed Loads on Ground Surface. Symmetry. 2019; 11(6):823. https://doi.org/10.3390/sym11060823

Chicago/Turabian Style

Li, Zejun, Jianchen Wang, and Kaihang Han. 2019. "Analytical Solution of Ground Stress Induced by Shallow Tunneling with Arbitrary Distributed Loads on Ground Surface" Symmetry 11, no. 6: 823. https://doi.org/10.3390/sym11060823

APA Style

Li, Z., Wang, J., & Han, K. (2019). Analytical Solution of Ground Stress Induced by Shallow Tunneling with Arbitrary Distributed Loads on Ground Surface. Symmetry, 11(6), 823. https://doi.org/10.3390/sym11060823

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