Face Stability Analysis of Tunnel with Pipe

Roof Reinforcement Based on Limit Analysis

Haitao Wang
PhD candidate, State Key Lab. of Coastal and Offshore Eng.,
Dalian University of Technology, Dalian, China
e-mail:[email protected].

Jinqing Jia
Professor, State Key Lab. of Coastal and Offshore Eng.,
Dalian University of Technology, Dalian, China
e-mail: [email protected]

ABSTRACT
Based on the kinematic method of limit analysis and the shear strength reduction technique,
the three-dimensional model for expressing the tunnel face stability with pipe roof
reinforcement was established. For a typical example, the solutions computed by the
proposed approach were compared with the results given by wedge model, trapezoid wedge
model and centrifugal-model test. The proposed approach was also employed to study how
the cover depth of tunnel, tunnel diameter and soil parameters affect the support pressure
and safety factor of tunnel face. It is shown that the calculation results of limit analysis is
smaller when compared with that by wedge model and trapezoid wedge model, and also
more closer to the results of centrifugal-model test. It is also found that for the relative
depth investigated ( / 0.5 ~ 8 H D = ), cover depth has a little influence on the stability of
tunnel face for a given diameter, but a good linear relationship exits between safety factor
and tunnel diameter. The linear and quadratic formulas are able to account for the influence
of cohesion and friction angle of soil on limited support pressure of tunnel face,
respectively. The analytical results of tunnel face stability also reveal that the safety factor
of tunnel face with pipe roof reinforcement is about 1.05~1.20 times of that without pipe
roof reinforcement. The effect of pipe roof reinforcement on tunnel face stability is not
significant.
KEYWORDS: pipe roof reinforcement; tunnel face stability; limit analysis;
strength reduction technique
Vol. 14 [2009], Bund. G 2
INTRODUCTION
The construction of tunnel in weak ground can give rise to some problems including
instability of the tunnel, damage to already existing structures, and so forth. In order to
overcome these problems, some auxiliary techniques are adopted in tunneling. One of such
methods is the pipe roof reinforcement method. This method consists on installing, prior to
the excavation of a length of tunnel, a series of pipes, either parallel to the tunnel axis or at a
certain angle with it. By injecting grout through the pipes, the ground in between the pipes is
stiffened and the pipes are connected, creating a kind of umbrella above the area to be
excavated. This arrangement creates a stiff layer of ground and allows safe excavation even in
poor ground conditions.
In recent years, there have been numerous studies into the behaviors of pipe roof
reinforcement. Ibrahim (2008) reported the field monitoring results from a tunneling project
involving pipe roof reinforcement technique. They found that with the use of pipe roof
reinforcement, the stability of the tunnel face in a poor geotechnical condition is improved
and the ground movement can be reduced effectively. Chang (1999) implemented the
analytical design method for the pipe roof reinforcement based on Winkler elastic foundation
beam theory. Considering the delay effect of initial lining and revising the Winkler elastic
foundation model, Wang and Jia (2008) investigated the mechanical behaviors of pipe roof
reinforcement based on Pasternak elastic foundation beam theory. The calculation results
prove that the Pasternak model is a better way to understand the reinforcement mechanism
and improve design practice. Yoo (2002) conducted a parametric study on the effect of
reinforcing layouts on the deformation behavior of the tunnel face and drew a conclusion that
there existed an optimum reinforcing layout to reduce the deformation of the tunnel for a
given tunnel geometry and ground condition. Hisatake and Ohno (2008) conducted a series of
centrifuge tests to clarify the effects of pipe roof supports and the excavation method on the
displacements above a tunnel face. The results of the model tests show the quantitative effects
of the ring-cut excavation method and clarify that the method decreases the ground
displacements remarkably. The model tests also indicate that the maximum settlement of a
ground excavated by the full-face excavation method with pipe roof supports is one fourth of
that without them. Recently, Shin et al. (2008) performed a large scale model testing for the
pipe-reinforced tunnel heading in a granular soil to understand reinforcing mechanism and
improve design practice. They found that the pipe reinforcement of heading increases
longitudinal load transfer to an unexcavated area, and consequently decreases deformation
and increases face stability.
Overview of the previous studies indicates that there are few attempts to evaluate the
tunnel face stability with pipe roof reinforcement quantitatively. In this study, based on the
kinematic method of limit analysis and the shear strength reduction, the three-dimensional
model for expressing the tunnel face stability with pipe roof reinforcement is established. For
Vol. 14 [2009], Bund. G 3
a typical example, the solutions computed by the proposed approach are compared with the
results given by wedge model, trapezoid wedge model and centrifugal-model test to verify the
reasonability of the method. The proposed approach is also employed to study how the cover
depth of tunnel, tunnel diameter and soil parameters affect the support pressure and safety
factor.
LIMIT ANALYSIS FOR TUNNEL FACE STABILITY
Three-dimensional analytical model
Leca and Dormieux (1990) used the limit analysis concept to evaluate the stability of a
tunnel face driven in frictional soil and compared these results with centrifuge tests performed
by Chambon and Corte (1994). A reasonable agreement was found between the theoretical
upper bound estimates and the face pressures measured at failure from the tests. Therefore, by
modifying the upper bound solution suggested by Leca and Dormieux (1990), the face
stability of tunnel with pipe roof reinforcement was evaluated quantitatively in this study. To
establish the three-dimensional analytical model, some key assumptions were necessary and
they were listed as follows:
1. The soil is modeled as a Mohr-Coulomb material, characterized by its cohesion c and
its friction angle throughout this study;
2. Failure mechanism was assumed to involve the movement of solid conical blocks with
circular cross-sections as shown in Fig. 1. Two cones were considered,
1
C (apex
1
, axis
1
)
of which the base
1
is in the same plane as the tunnel face; and
2
C (apex
2
, axis
2
) of
which the base
12
is in plane (
1
C and
2
C have the same geometric properties). The first
moving block
1
B corresponds to the portion of
1
C located below plane and
2
B is the portion
of
2
C located below the lower surface of the pipe roof reinforcement.
Vol. 14 [2009], Bund. G 4

Figure 1: Analytical model for face stability of tunnel with pipe roof reinforcement

3. According to different overburden pressure above the tunnel face, two kinds of the
analytical models were established for tunnel with pipe roof reinforcement and without pipe
roof reinforcement. The overburden pressure above the tunnel face for both kinds was
calculated as follows:
Tunnel without pipe roof reinforcement: The overburden pressure above the tunnel face
can be calculated by Terzaghi earth pressure theory. The mechanical model was shown in
Fig.2 (a);
Tunnel with pipe roof reinforcement: The overburden pressure above the tunnel face is
the counterforce of foundation, which can be obtained from reference 3. The mechanical
model was shown in Fig.2 (b)
Upper-bound solutions for the face stability of tunnel
In order to satisfy the stability condition of the tunnel face according to the upper bound
theorem, the following relation was considered:

e v
P P
(1)
Where P
e

represents the power of the external loads, which is the sum of three
components, P
T
, P
S
, and P

.

P
T
is the power of the support pressure
T
, P
S
is the power of
the surcharge
S
and P

is the power of the soil weight .

v
P denotes the dissipation
power. The derivations for each component of the upper bound theorem have been omitted.
Vol. 14 [2009], Bund. G 5
The upper bound solution for the collapse mechanism of two conical blocks is described
following Leca and Dormieux (1990) as:
2 3 2 2
2 2 2 2 2 2
R R cos ( ) cos cos cos cos( )
sin
4 cos cos ( ) R R 2sin sin( ) 2sin cos ( ) 3
e T S
E A B E
S A T C
C C
P P P P
R R D D
R R V

= + +
( | | + +
= + + ( |
|
+ + +
( \ .
(2)
where

cos( ) cos( )
cos
A
R

+
=
(3)

cos( ) cos( )
sin 2
B
R

+
=
(4)

1/ 2
cos( ) sin( )
cos sin( )
C
R


( +
=
(
+

(5)

sin
sin sin( )
D
R


=
+
(6)

sin( ) cos( )
sin( )
E
R


+
=
+
(7)

(a)Tunnel without pipe roof reinforcement (b) Tunnel with pipe roof reinforcement
Figure 2: Calculating chart of overburden pressure above tunnel face
The energy is dissipated along the lateral surfaces of the moving block
1
B and
2
B and the
base
12
. The dissipation power can be expressed as equation (8).
Vol. 14 [2009], Bund. G 6

1 2 12
2 2 2
2 2 2 2
cos ( ) cos
cos
4 cos sin sin cos ( )
v v v v
A E
C
P P P P
R R D
c V
R

= + +
( +
=
(
+

(8)
where,
1v
P ,
2v
P ,
12v
P are the respective contributions of
1
B and
2
B and
12
.
The upper bound solution for the face stability of the tunnel with pipe roof reinforcement
is obtained by substituting equation (2) and equation (8) in equation (1), resulting in:

( 1) 1 ( 1) ( 1) 1
S T
S P P P
C C C
D
N K N K K

(
+ + +
(

(9)
where
C
is unconfined compression strength of soil,
P
K is Rankines passive earth
pressure coefficient,
S
N and N

are weighting coefficients.

cos
2
1 sin
C
c

(10)

1 sin
1 sin
P
K

+
=

(11)

2
2
1 sin( )
cos cos sin( )
E
S
A
R
N
R

=
+
(12)
3 3
2
1 cos cos( ) 1 sin( )
tan
3 2sin sin( ) 2sin cos cos sin( )
C E
B
A A
R R
N R
R R

( +
= +
(
+ +

(13)
S
N and N

are the functions of the friction angle . Leca and Dormieux (1990) found that
the optimal values
c
s
N and
c
N

were always obtained for essentially the same value

c
of
parameter : 49 / 2
c
=

. This means that in the plane of symmetry of the tunnel, the

angle of the critical failure surface to the horizontal is 49 / 2 +

,which is larger than the

angle of active failure in plane strain conditions, 45 / 2 +

. Therefore the area in front of the

tunnel which is influenced by the collapse is more limited than in the case of a long open cut.
This could be seen as a stabilizing effect due to three-dimensional equilibrium conditions
Vol. 14 [2009], Bund. G 7
around the tunnel face. Since plane is chosen such that
2
is vertical, we have between and
the relation: 2 / 2 = .
Safety factor of tunnel face
The strength reduction technique for finite element analysis of slope stability has been
applied successfully by many authors (Huang, 2009; Wei, 2009). This method allowed
finding the safety factor of a slope by initiating a systematic reduction sequence for the
available shear strength parameters c and to just cause the slope to fail. The reduction
values of shear strength parameters
m
c and
m
were defined as:

m
c
c
K
=
,
tan
tan
m
K

=
(14)
in which the
K
is the strength reduction factor. The safety factor of the slope,
S
K
, is the
value of K to bring the slope to failure.
In this study, the assumption was that when
T
equal to zero, the tunnel face was in limit
equilibrium status. As
T
is the non-linear function of the c and , using the iterative method
to reduce the shear strength parameters until the
T
to zero, the safety factor of the tunnel
face,
S
K , is the value of strength reduction factor.
CALCULATION CASE
There have been many methods to calculate the support pressure of tunnel face and to
evaluate the tunnel face stability. Chambon and Corte (1994) carried out centrifuge tests to
study the face stability of tunnels in sands. Qin (2005) using wedge model to calculate the
support pressure of tunnel face. Wei and He (2007) modified the wedge model and assumed
the up side of sliding block as a trapezoid prism, using the Terzaghi loose earth pressure
theory to calculate limited support pressure of tunnel face based on the principle of sliding
block entire force balance. In this study, the solutions computed by the proposed approach
were compared with the results given by wedge model, trapezoid wedge model and
Vol. 14 [2009], Bund. G 8
centrifugal-model test to verify the reasonability of the method. The results obtained from the
four methods were summarized in Table 1.
Table 1: Comparison of the limited support pressures

D/
m
H/
m
c/
kPa
/
kNm-3
Limited support pressures (kPa)
Centrifuge
test
Wedge
model
Trapezoid wedge
model
Limit analysis
approach
5 2.5 0 16.1 3.6 9.262 7.766 6.524
5 5 0 16.1 3.3 10.623 8.287 6.805
5 10 0 16.1 4 11.280 8.404 7.266
10 10 0 16.0 7.4 21.115 16.472 13.526
10 20 0 16.0 8.0 22.420 16.705 13.561
10 40 0 16.0 8.2 22.606 16.713 13.629
As seen from Table 1, it is apparent that the limited support pressures calculated by
theoretical approach is larger than the results of centrifuge tests, but the calculation results of
limit analysis are in close agreement with test results. In sandy soil, the progression of failure
of tunnel face is a gradual process. When the cover depth is large enough, failure does not
reach the ground surface, so the earth pressure calculated by Terzaghi earth pressure theory is
usually larger than actual value, this maybe the reason to explain why the support pressures
calculated by theoretical approach is larger than the results of centrifuge tests.
Other similarities between upper bound solutions and experimental results were shown in
Fig.3, in which the failure zone observed in the centrifuge along the tunnel centre plane was
represented for the case of a loose sand for H/D =1.0. The critical geometry associated with
the upper bound solution was shown with dashed lines. Even though it does not extend in the
vertical direction as much as the actual failure area, it coincides almost perfectly with the
observed surface in front of the tunnel.

Figure 3: Comparison between theoretical critical surface and observed failure area
Vol. 14 [2009], Bund. G 9
PARAMETRIC STUDY
In this study, to optimize design and construction of tunnel with auxiliary techniques, the
proposed approach was employed to study how the tunnel cover depth, tunnel diameter and
soil parameters affect the support pressure and safety factor. Two kinds of the analytical
models were considered, tunnel with pipe roof reinforcement and without pipe roof
reinforcement. Steel pipes used for pipe roof reinforcement have outer diameter of 108mm,
length of 30m and thickness of 6 mm. The distance between the pipes is 40 cm and overlap
length of pipe is 1.5m. A unit weight of 18 = kN/m
3
was assumed for soil with the cohesion
of 20kPa and the friction angle of 35

.
Influence of cover depth of tunnel
It is assumed that cover depth of tunnel ranges from 2.5m to 40m. In this study, tunnel
diameter is 5m, therefore, / H Dranges from 0.5 to 8. The relation between face support
pressure (
T
), safety factor (
S
K ) and cover depth of tunnel (H) were shown in Fig.4 and
Fig.5, respectively.

Figure 4: Relation between
T
and H
Vol. 14 [2009], Bund. G 10

Figure 5: Relation between K
s
and H
It can be seen in Fig. 4 that the support pressures of tunnel face are minus. These mean
that the tunnel face can maintain stability and doesnt need any auxiliary techniques. For the
relative depth investigated ( / 0.5 ~ 8 H D = ), depth does not have a large influence on the
limited support pressure of tunnel face and safety factor for a given diameter. Such conclusion
is well supported by the results of centrifuge tests reported by Chambon and Corte (1994).
The effect of pipe roof reinforcement on the stability of tunnel face can also be seen from
Fig.4 and Fig.5. When H is 2.5m, the
S
K of tunnel face for reinforced case is 2.137, and 1.835
for unreinforced case; When H is 40m,
S
K of tunnel face for reinforced case and
unreinforced case are 1.797 and 1.561, respectively. It can be understood that the safety factor
of tunnel face with pipe roof reinforcement is about 1.15 times of that without pipe roof
reinforcement.
Influence of tunnel diameter
It is assumed that tunnel diameter ranges from 5m to 15m. In this study, the cover depth
of tunnel was 20m. The relation between face support pressure (
T
), safety factor (
S
K ) and
tunnel diameter (D) were shown in Fig.6 and Fig.7, respectively.
Vol. 14 [2009], Bund. G 11

Figure 6: Relation between
T
and D

Figure 7: Relation between
S
K and D
From Fig.6 and Fig.7, it appears that the limited support pressure of tunnel face increase
linearly and safety factor decrease linearly in according with the tunnel diameter. The
relationship between
T
,
S
K and D can be described as following formula:

0.098 25.581
T
D =
(15)

0.003 1.782
s
K D = +
(16)
When D is 5m, the K
s

of tunnel face for reinforced case is 1.858, and 1.585 for
unreinforced case; When D is 15m, K
s
of tunnel face for reinforced case and unreinforced
case are 1.055 and 0.955, respectively. These indicate that although pipe roof reinforcement
Vol. 14 [2009], Bund. G 12
can be expected to support the soil pressure induced by the loosening area ahead of the face, it
may have limited value on face stability in large excavation span conditions.
Influence of soil parameters
In order to evaluate the tunnel face stability with different soil parameters, a systematic
parametric study was conducted. It is assumed that cohesion of soil ranges from 5kPa to
80kPa, friction angle of soil ranges from 20
o
to 40
o
. The relation between safety factor (K
s
)
and soil parameters were listed in Table 2, and the relation between face support pressures (
T
) and soil parameters were shown in Fig.8 and Fig.9, respectively.

Table 2: Relation between
S
K
and soil parameters
S
K

/ c kPa /


5 10 20 50 80 20

25

30

35

40

Unreinforced
Reinforced
0.824
0.866
1.113
1.231
1.585
1.858
3.271
4.021
6.801
7.607
0.978
1.238
1.168
1.439
1.368
1.643
1.585
1.858
1.817
2.075

From Table 2, it can be noted that the stability of the tunnel face is improved with
increasing c and . It indicates that the improvements of soil conditions can greatly
improve the stability of tunnel face. It was also found that the safety factor of tunnel face with
pipe roof reinforcement is about 1.05~1.20 times of that without pipe roof reinforcement. The
effect of pipe roof reinforcement on tunnel face stability is not significant. As presented in
Fig.8, it was found that a good linear relationship exists between cohesion of soil (c) and face
support pressure (
T
). From Fig. 9, it appears that the quadratic formula is able to account
for the influence of friction angle of soil ( ) on face support pressure (
T
).The relationship
between c , and
T
can be described as following formula:

1.425 11.237
T
c = +
(17)

2
0.084 5.927 86.415
T
= +
(18)

Vol. 14 [2009], Bund. G 13

Figure 8: Relation between
T
andc

Figure 9: Relation between
T
and

CONCLUSIONS
(1) Based on the kinematic method of limit analysis and the shear strength reduction
technique, the three-dimensional model for expressing the tunnel face stability with pipe roof
reinforcement is established. This method can be employed to define the safety factor and its
corresponding critical failure mechanism for a given tunnel.
(2) For a typical example, the solutions computed by the proposed approach were
compared with the results given by wedge model, trapezoid wedge model and
centrifugal-model test. The results indicate that limit analysis approach can provide
reasonable estimates of tunnel face stability.
Vol. 14 [2009], Bund. G 14
(3) It is found that for the relative depth investigated ( / 0.5 ~ 8 H D = ), depth does not
have a large influence on the safety factor of tunnel face for a given diameter, but a good
linear relationship exists between safety factor of tunnel face and tunnel diameter. It is also
found that the linear and quadratic formulas are able to account for the influence of cohesion
and friction angle of soil on limited support pressure of tunnel face, respectively.
(4) The safety factor of tunnel face with pipe roof reinforcement is about 1.05~1.20 times
of that without pipe roof reinforcement. The effect of pipe roof reinforcement on tunnel face
stability is not significant.
REFERENCES
1. Ibrahim, O. (2008) Control of surface settlements with umbrella arch method in
second stage excavations of Istanbul Metro, Tunnelling and Underground Space
Technology, 23(6), pp 674681.
2. Chang, S.D.(1999) Research on pre-brace mechanism of pipe umbrella method,
Chengdu: Southwest Jiaotong University, (in Chinese)
3. Wang, H.T., Jia J.Q. (2008) Analytical method for mechanical behaviors of pipe
roof reinforcement, International Conference on Information Management,
Innovation Management and Industrial Engineering. Taipei: IEEE Computer society,
pp 352357.
4. Yoo, C. S. (2002) Finite-element analysis of tunnel face reinforced by longitudinal
pipes, Computers and Geotechnics, 29(1), pp73-94.
5. Hisatake, M., Ohno, S. (2008) Effects of pipe roof supports and the excavation
method on the displacements above a tunnel face, Tunnelling and Underground
Space Technology, 23(2), pp 120-127.
6. Shin, J. H., Choi, Y. K., Kwon, O. Y., Lee, S. D. (2008) Model testing for
pipe-reinforced tunnel heading in a granular soil, Tunnelling and Underground
Space Technology, 23(3), pp 241-250.
7. Leca, E., Dormieux, L. (1990) Upper and lower bound solutions for the face stability
of shallow circular tunnels in frictional material, Geotechnique, 40 (4), pp 581606.
8. Chambon, P., Corte, J. F. (1994) Shallow tunnels in cohesionless soil: Stability of
tunnel face, Journal of Geotechnical Engineering, 120 (7), pp 11481165.
9. Huang M.S., Jia C. Q. (2009) Strength reduction FEM in stability analysis of soil
slopes subjected to transient unsaturated seepage, Computers and Geotechnics,
36(1), PP 93-101.
Vol. 14 [2009], Bund. G 15
10. Wei, W. B., Cheng, Y. M., Li, L. (2009) Three-dimensional slope failure analysis by
the strength reduction and limit equilibrium methods, Computers and Geotechnics,
36(1), pp70-80.
11. Qin, J. S. (2005) Study on face movement and collapse of earth pressure shield
tunnel, Nanjing: Hohai University. (in Chinese)
12. Wei, Gang., He, F. (2007) Calculation of minimal support pressure acting on shield
face during pipe jacking in sandy soil, Chinese Journal of Underground Space and
Engineering, 3(5): pp903-908. (in Chinese)

2009 ejge