# Field Observation and Theoretical Study on an Existing Tunnel Underpassed by New Twin Tunnels.

1. IntroductionThe interaction between new shield construction and existing tunnel has become a common and important issue with the rapid development of underground traffic system, which has been studied in the past using a variety of approaches: field observations, model tests, analytical methods, and finite element modeling. Kim and Liu et al. [1, 2] presented a good summary comparison of the studies and only studies that illustrates the shield underpassing or parallel underpassing the existing tunnel. Yamaguchi et al. [3] presented successively the numerical model and then analyzed three configurations of the twin tunnels in Japan: aligned horizontally, vertically, and inclined. The construction of the upper tunnel at first leads to both higher settlement and bending moment. The maximum soil settlement was obtained for vertical-aligned tunnels, while horizontal-aligned tunnels caused the lowest settlement. Addenbrooke and Potts [4] analyzed the influence of tunnel position, tunnel spacing, rest period, and sequence of excavation on the interaction between the two tunnels. It concentrated on the shape of the settlement profile and the volume loss induced by the two tunnels. The shield underpassing the existing tunnel is highly site specific. The soil condition, the tunnel buried depth, and the relative position of the new tunnel and existing tunnel all affect the response of the existing tunnel. Field observations remain the commonly recognized approach for understanding the interaction behavior between the new tunnel and existing tunnel. Usually, the underpassing of a shield with a loss of soil often causes settlement of the existing tunnel at last [5-7]. Chehade and Shahrour [8] used the finite element method to investigate the influence of the relative position of tunnels and the construction procedure on the soil settlement. The results showed that the settlements of the existing tunnel for the vertical parallel tunnels were larger than those for the horizontal parallel tunnels. Li and Yuan [9] studied the twin tunnels passing under a double-decked tunnel at an angle of 55 in weathered granite gneiss. Only settlements were found, and the horizontal displacement was smaller than the vertical displacement in the existing tunnel. Despite a number of studies having been carried out, current understanding of the interaction between the two tunnels is still limited.

In order to ensure the stability of the existing tunnel, local thickening is needed at the sides of the existing concrete lining. There are a plenty of research about the conventional reinforcement methods such as inner steel plate reinforcement or new reinforcement approaches such as FRC or composite concrete, but they all focus on the performance of reinforcement on the single ring. The effect of reinforcement on the longitudinal behavior of the tunnel is not yet clear in current analysis. The analytical method of longitudinal displacement of the tunnel due to adjacent excavation and multiple tunnelling is researched a lot: (1) the elastic continuum models developed by Vorster et al. [10], (2) the Winkler model [11], and (3) the two-parameter elastic models. When variable stiffness of the existing tunnel needs to be considered, the methods are no longer applicable. Selvadurai [12] divided the beam on foundation into N elements and used the method of initial parameters to analyze the beam's displacement. In the paper, we assume the existing tunnel as a beam resting on a two-parameter foundation, combining the initial parameter method and the transfer matrix method to analyze the reinforcement effect on the longitudinal behavior of the existing tunnel.

The major objectives of this paper are (1) to investigate the influence of the tunnel driving parameters and the relative distance between the shield and the existing tunnel on the existing tunnel based on the interpretation of the field measured data and (2) to study the effects of different stiffness and ranges of reinforcement on the existing tunnel using analytical methods.

2. Project Overview

The location of the Line 4 tunnels and the upline of Line 1 are shown in Figure 1. The upline of Line 1 is the existing tunnel built in 2012. Figure 2 shows a plan view of tunnel alignment and arrangement of the monitoring points in the existing tunnel. The shield of the northbound tunnel of Line 4 starts from Guanhe Station, passes under the existing tunnel with a small angle of 23[degrees], and reaches East railway station at last. Then, the shield machine is reassembled with tunnels in a reversed direction. The northbound tunnel with a length of 329.675 m consists of 275 rings numbered starting with zero from Guanhe Station to East railway station, while the southbound tunnel with a length of 323.834 m consists of 270 rings numbered starting with zero from East railway station to Guanhe Station. The lateral distance between the northbound tunnel and the southbound tunnel is 9.4-15.0 m.

Figure 3 shows the longitudinal profile of soil and tunnels. The vertical distance between the existing tunnel and Line 4 is 2.12 m. The existing tunnel and new tunnels are all built using an articulated shield tunnelling machine with an outer diameter of 6.4 m and a length of 8.5 m. The spokes plus panel-type cutter head are used, and the aperture ratio of cutter head is 40%. Each ring of new tunnels consists of six precast concrete segments. The outer diameter, the inner diameter, and the thickness of the segment are 6.2 m, 5.5 m, and 0.35 m, respectively.

2.1. Soil Condition. The engineering properties of the rock and the soils in this site are very complicated. Figure 4 shows geotechnical parameters ofsoil in this site. The soil layers from top to bottom are (1) fill, (3) silty clay, (4) clay, (6) the muddy silty clay, and (8) soft silty clay and gravel layer. The existing tunnel and new tunnels are located in the muddy silty clay layer at a depth of 20 m and 25 m, respectively. Along the metro line, an extensive geotechnical investigation is carried out. The standard penetration test (SPT), vane shear test (VST), piezocone test (CPTU), flat plate dilatometer test (DMT), and water pumping test (WPT) are conducted along Line 1. The maximum water content of clay and muddy silty clay is about 40%, which is close to the liquid limit value. The average undrained shear strength of undisturbed clay and remolded clay which is obtained by VST is 70.5 kPa and 10.9 kPa, respectively. The coefficient of the at-rest earth pressure (K0) of soft clay is obtained by DMT.

2.2. Hydrological Geology. The shallow ground water is pore phreatic water, mainly found in layers from (1) to (8). The elevation of the water surface is 2.232 m. The laboratory penetration test and field steady flow test can be seen in Table 1. It can be seen that the permeability coefficient of the clay layer is very small, and the permeability coefficient of the gravel layer is 1-2 m/d. The seepage induced from water head difference might influence the stability of tunnels.

3. Shield Driving Parameters

Figure 5 presents advance-time curve of the shield. Before underpassing the existing tunnel, the shield advances with a very slow rate. The tunnel advancing rate is 7-8 rings/day, and the penetration rate of the shield ranges from 20 to 25 mm-min-1.

Figure 6 shows the applied tunnel face pressure. In practice, the tunnel face pressure should be adjusted with the theoretical value as well as the feedback displacement of the existing tunnel. The existing tunnel redistributes the soil stress, and the applied tunnel face pressure should be adjusted by (1) as follows:

[p.sub.2] = [K.sub.0][gamma]'h + [p.sub.w] + [Q.sub.1] - [Q.sub.2] [+ or -] 20 kPa, (1)

where [K.sub.0] is the coefficient of earth pressure at rest; [gamma]' is the effective gravity of soil; h is the thickness of overburden depth (m); [p.sub.w] is the water pressure; Q1 is the weight of the existing tunnel (kN); [Q.sub.2] is the weight of the soil with the same internal volume of the existing tunnel (kN); and 20 kPa is the pressure fluctuation.

The applied tunnel face pressure of the southbound tunnel cannot be calculated by (1) which is due to the redistribution of the earth pressure after the construction of the northbound tunnel nearby and should be adjusted according to the displacement of the existing tunnel.

Figure 7 shows the tail void grouting volume of the shield of the northbound and southbound tunnels of Line 4. Shield parameter adjustments play an important role in the displacement controls, reducing the tail void grouting volume and applied tunnel face pressure when Line 1 heaves. Opposite measures are taken when there is a subsidence in Line 1.

4. Observation Results and Discussions

4.1. Monitoring Arrangement of Line 1. Arrangement of the monitoring rings of the existing tunnel can be seen in Figure 2. The total station Topcon MS05AX fixed on the tunnel sidewall arranged from rings 404 to 559 in the existing tunnel is an automatic real-time measuring system which is used to monitor vertical and horizontal displacement. One monitor section is set every two rings in the most affected zone (from rings 447 to 499 in the existing tunnel), and one monitor section is set every five rings in the rest part. An LECAI D5 hand-held distance finder is used to monitor the converge displacement every 5 rings from 404 to 679 in the existing tunnel. The arrangement of monitoring points at the cross section of the existing tunnel is shown in Figure 8.

4.2. Displacement of the Existing Tunnel. The location of point A as shown in Figure 2 is corresponding to the intersection point of a plan view of the existing tunnel and northbound tunnel of Line 4. The location of point B is corresponding to the intersection point of the plan view of the existing tunnel and southbound tunnel of Line 4. The middle point of points A and B is point C.

4.2.1. Vertical Displacement of the Existing Tunnel. Figure 9 shows time-varying vertical displacement in the monitoring rings of the existing tunnel. The x-axis is the ring number of the existing tunnel. A positive value of the ordinate denotes heave, while a negative value denotes settlement of the existing tunnel. The selected monitoring rings start to heave when the shield face is 0-10 m away from the selected monitoring rings which is mainly due to large applied tunnel face pressure, bulk addictive thrust (Figure 6), and the friction force between the shield shell and soil mass. When the shield tail is far beyond the selected monitoring rings, a reduction of heave in the selected monitoring points is observed which is mainly due to the closure of the tail void. Only small settlements and heaves are measured after the construction of the northbound tunnel. During the 3-month shutdown of the shield, additional settlement ranges from 2 mm to 3 mm can be observed in rings from 430 to 520. Additional settlements ranges from 2 to 4 mm were measured in rings from 460 to 540 during the southbound tunnel's construction. The long-term additional settlements, monitored up to 140 days (from 2014/6/1 to 2014/10/19), range from 2 to 4mm in rings from 460 to 510. The settlements of rings from 430 to 520 range from 2 mm to 12 mm, and the settlement curve of the upline of Line 1 is "U" shaped after the long-term monitoring. The settlement curve is approximately symmetric about the dashed line C after the long-term monitoring. The maximum settlement is 12 mm which is located in ring 487.

4.2.2. Horizontal Displacement of the Existing Tunnel. The horizontal displacement of the existing tunnel with respect to the location of the shield is illustrated in Figure 10. A positive horizontal displacement denotes northward transverse tunnel movement away from the original tunnel centerline, while a negative horizontal displacement denotes southward transverse tunnel movement away from the original tunnel centerline.

There are northward displacement on the left side of point C and southward displacement on the right side of point C. The horizontal displacement curve is approximately symmetric about point C after the long-term monitoring. The maximum northward displacement is 10 mm in the ring near the intersection point A, and the maximum southward displacement is -10.5 mm in the ring near the intersection point B after the completion of the two tunnels' construction. During the construction of the northbound tunnel, rings from 453 to 471 move southward which is likely due to the additional bulkhead additive thrust and the squeezing force provided by the shield shell. When the shield tail leaves ring 459, rings from 439 to 459 move southward slowly. Only 2 to 5 mm additional southward displacement is measured in rings from 487 to 560, and nearly no displacement is observed on the left side of point B during the construction of the southbound tunnel which is mainly due to the northbound tunnel's barrier effect. No change of horizontal displacement is observed in the long-term conditions (from 2014/5/28 to 2014/10/19).

4.2.3. Convergence Displacement of the Existing Tunnel. Figure 11 shows the convergence displacement of the existing tunnel with respect to the locations of the shield. A negative value indicates the reduction horizontal diameter of the existing tunnel, while a positive value indicates the addition horizontal diameter of the existing tunnel. It can be observed that the convergence displacement is not obvious when the shield reaches ring 48 (22 rings away from the intersection point A). During the process of the shield driving from ring 64 to ring 87, a significantly additional increase in the negative convergence displacement in rings from 440 to 490 is observed, and the maximum displacement occurs in the intersection point A. The reason might be that the shield face squeezes one side of the existing tunnel, leading to the reduction of the horizontal diameter of the existing tunnel. With the tunnelling of the southbound tunnel, there are a reduction negative convergence displacement on the left side of point B and an addition positive convergence displacement on the right side of point B. As can be seen from Figures 10 and 11, the shapes of the convergence displacement curve and the horizontal displacement curve are similar because the causes of convergence displacement and horizontal displacement are theoretically the same.

5. Reinforcement Scheme of the Existing Tunnel

Figure 12 shows the reinforcement in the existing tunnel. Radial reinforcement by an arc-shaped supporting steel plate connected to the tunnel segment and longitudinal reinforcement by channel section steel to provide longitudinal tensile stress are conducted in the existing tunnel. Radial steels and longitudinal steels are connected by welding, and so, the reinforcement becomes a whole one. 25 rings are reinforced at first on either side of the intersection point A before the underpassing of the northbound tunnel, and the whole reinforcement is completed after the construction of the northbound tunnel. The reinforcement range is from ring 424 to ring 510 (i.e., 25 rings on either sides of the intersection points A and B).

6. Theoretical Analysis of the Reinforcement Design

The deformed tunnel lining reinforced by inner bonding steel plates has been partially or entirely used in many built tunnels or currently under construction. To better understand the benefits coming from the reinforcement method, the behavior of the tunnel is investigated.

Figure 13 shows a schematic view of the existing tunnel-soil tunnelling interaction. The analysis method demonstrated in this paper can be divided into two steps: firstly, estimating the greenfield displacement induced by the tunnelling. Secondly, calculating the responses of the existing tunnel subjected to the soil displacement. The method of analysis is based on three assumptions: (1) the above existing tunnel does not affect the displacement of soil due to tunnelling, (2) the soil foundation is assumed as the Winkler or the Pasternak foundation, and (3) the soil displacement is calculated by superimposing the independent settlement predicted for each individual tunnelling.

In this paper, the tunnels are considered as an infinite beam on the Winkler foundation, an infinite beam on the Pasternak foundation, and a finite beam on the Pasternak foundation (the new method). Comparing the deflection, rotation angle, normalized bending moment, and shear force of the existing tunnel with constant stiffness based on the three models, the new method is verified. Then, different stiffness of the reinforcement is taken into consideration, and the optimal reinforcement range of the existing tunnel is discussed.

6.1. The Subsurface Soil Displacement due to Tunnelling. For the theoretical analysis, the alignment of the new tunnels and the existing tunnel is assumed to be straight. Figure 14 shows a schematic diagram of the new tunnels and the existing tunnel in this case. The greenfield settlement s(x) due to the tunnelling can be replaced by the equivalent distributed load q(x) acting on the beam on the third assumption as follows:

q(x) = ks(x). (2)

In this study, subsurface settlement s(x) at the depth of z induced by tunnelling is calculated based on closed-form analytical solutions presented by Loganathan and Poulos [13] as follows:

[mathematical expression not reproducible], (3)

where s(x) = subsurface settlement of the soils due to tunnelling; R = tunnel radius; z = depth below ground surface; H = depth of the axis of the new tunnel; v = Poisson's ratio; x = lateral distance from tunnel centerline; and [[epsilon].sub.0] = equivalent ground loss ratio which is defined as

[[epsilon].sub.0] = 4Rg + [g.sup.2]/4[R.sup.2], (4)

g = [G.sub.p] + [u.sup.*.sub.3D] + w [approximately equal to] [G.sub.p], (5)

where g is the gap parameter (Lee et al. [14]); [u.sup.*.sub.3D] is the three-dimensional elastoplastic deformation at the tunnel face; w is the workmanship factor; and [G.sub.p] represents the physical gap between the outer skin of the shield and the lining which is given as

[G.sub.p] = 2[DELTA] + [delta], (6)

where [DELTA] is the thickness of the tail piece and [delta] is the clearance required for erection of the tail piece. The parameter w can be neglected due to considerable experience with the equipment and good tunnelling technique of construction workers based on the study of Rowe and Lee [15].

6.2. Analytical Methods Derived from Prior Studies. Attewell et al. [16] proposed a tunnel-soil interaction model and solved the analytical solution of the longitudinal displacement of an infinite tunnel:

[d.sup.4]w(x)/d[x.sup.4] + 4[[lambda].sup.4] p(x) = 4[[lambda].sup.4]q(x), (7)

where q(x) is the concentrated load acting on the origin point of the beam; EI is the flexural stiffness of a tunnel, which is recommended by Ye et al. [17]; and k is the subgrade modulus that represents the pipe-soil interaction.

If the soil foundation is assumed as the Winkler foundation, the foundation counterforce p(x) can be expressed as

p(x) = kw(x) (8)

The two-parameter model of foundation can capture the shear resistance of soil. Based on the assumption of the plane strain condition, the displacement of ground in the y-direction is neglected. If the soil foundation is assumed as the Pasternak foundation, the foundation counterforce p(x) can be expressed as

p(x) = -G[[nabla].sup.2]w(x) + kw(x), (9)

where G is the coefficient of the shear element in the Pasternak's model with the dimension of force per unit length.

Attewell et al. [16] presented a solution for the concentrated load q on an infinite beam resting on the Winkler foundation and computed the deflection as

w(x) = 1/8EI[[lambda].sup.3] [qe.sup.-[lambda]x] (cos [lambda]x + sin [lambda]x). (10)

It can be concluded that each of the concentrated load ks([tau]) at point [tau] contributes the following amount to the deflection of the existing tunnel on the Winkler foundation:

[mathematical expression not reproducible], (11)

where [L.sub.1] is the distance between the coordinate origin of the existing tunnel and the intersection point between the new tunnel and the existing tunnel and [alpha] is the intersection angle between the new tunnel and existing tunnel. The deflection, rotational angle, bending moment, and shear force of the existing tunnel which is assumed as an infinite beam on the Winkler foundation and the Pasternak foundation are illustrated in Appendix A.

6.3. A Brief Description of the New Method. In this study, we want to figure out the effects of the reinforcement on the existing tunnel, the stiffness varies in this case and the existing tunnel assumed as an infinite beam is not suitable. Assuming the existing tunnel with finite length and introducing (7) to (9), the basic differential equation governing the flexure of the beam resting on the Pasternak foundation can be written in the form with an assumption of G[[lambda].sup.2]/k < 1:

[mathematical expression not reproducible], (12)

where [b.sub.1i] = [b.sub.i](1 + ([square root of [G.sub.i]/[k.sub.i]]/[b.sub.i])), subscript i represents the ith element of the beam, b is the width of the beam, and [b.sub.1i] is the modified width of the beam. Assuming a solution of (12) in the form of w(x) = [e.sup.[lambda]x], the general solution of vertical displacement of the beam is as follows:

w(x) = [4.summation over (i=1)] [C.sub.i][F.sub.i] (x), (13)

where [C.sub.i] (i = 1, 2, 3, 4) are undetermined parameters, which are dependent on the boundary conditions, and [F.sub.i] (x) (i = 1, 2, 3, 4) are the basic functions of the general solution, which are given as follows:

[F.sub.1i] (x) = cos ([a.sub.1i][c.sub.i]x) sinh ([a.sub.2i][c.sub.i]x), (14)

[F.sub.2i] (x) = cos ([a.sub.1i][c.sub.i]x) cosh ([a.sub.2i][c.sub.i]x), (15)

[F.sub.3i] (x) = sin ([a.sub.1i][c.sub.i]x) cosh ([a.sub.2i][c.sub.i]x), (16)

[F.sub.4i] (x) = sin ([a.sub.1i][c.sub.i]x) sinh ([a.sub.2i][c.sub.i]x), (16)

where [mathematical expression not reproducible].

The relationship of the deflection, rotational angle, bending moment, and shear force of the beam proposed by Lancaster and Mitchell [18] is as follows:

[mathematical expression not reproducible], (18)

where [[theta].sub.i](x), [M.sub.i](x), and [Q.sub.i](x) are the amplitudes of rotational angle, bending moment, and shear force of the section area, respectively.

Figure 15 shows that different stiffness is replaced by equivalent stepped stiffness. According to different stiffness and distributed load, the local coordinate system is established as illustrated in Figure 16. Selvadurai [12] suggests an initial parameter method to solve the parameters of the beam. The four parameters at the origin of the coordinates are [w.sub.i] ([0.sub.i]), [[theta].sub.i]([0.sub.i]), [M.sub.i] ([0.sub.i]), and [Q.sub.i] ([0.sub.i]), and four parameters at the end of the element are [w.sub.i]([x.sub.i]), [[theta].sub.i]([x.sub.i]), [M.sub.i]([x.sub.i]), and [Q.sub.i]([x.sub.i]). [C.sub.1], [C.sub.2], [C.sub.3], and [C.sub.4] can be expressed by [w.sub.i]([0.sub.i]), [[theta].sub.i]([0.sub.i]), [M.sub.i] ([0.sub.i]), and [Q.sub.i] ([0.sub.i]) by introducing [x.sub.i] = 0 to (13) and (18).

[w.sub.i] ([x.sub.i]), [[theta].sub.i] ([x.sub.i]), [M.sub.i] ([x.sub.i]), and [Q.sub.i] ([x.sub.i]) can be rewritten in the matrix-array form as

[mathematical expression not reproducible], (19)

where [[k.sub.ij]] is a stiffness submatrix with 5 x 5 order; i, j = 1,2, ..., 5. M and V are, respectively, the bending moment and the generalized shear force. Formula (19) can be simplified as

[N.sub.i] ([x.sub.i]) = [K.sub.ij][N.sub.i]([0.sub.i]), (20)

where [mathematical expression not reproducible].

Figure 17 shows that curved distributed load is replaced by equivalent triangular and square distributed load which can be expressed as

[q.sub.i-1] + [DELTA][q.sub.i] = [q.sub.i-1] + [xi] ([q.sub.i] - [q.sub.i-1]/[L.sub.i]), (21)

where [q.sub.i-1] and q, are, respectively, the load of the starting point and the ending point of the ith element and [xi] is the position of arbitrary load. From (19), we can obtain

[mathematical expression not reproducible]. (22)

In order to avoid numerical errors in the calculation process, the integration of [k.sub.i5] is shown in Appendix B. The rest of [k.sub.ij] can be referred in the work by Selvadurai [12].

Figure 18 shows the tunnel matrix transfer diagram. Assuming that the tunnel is composed of N elements, the x-coordinate of beginning points of the ith element is [x.sub.i](i = 0,1,2, ... n) and the length of each element is [L.sub.i], where [L.sub.i] = [x.sub.i] - [x.sub.i-1].

For the first element,

[N.sub.1] (x) = [K.sub.1] (x) x [N.sub.1] ([0.sub.1]) x [member of] (0, [L.sub.1]). (23)

For the end point of the first element,

[N.sub.1] ([x.sub.1]) = [K.sub.1] ([x.sub.1] - 0) x [N.sub.1] ([0.sub.1]). (24)

For the end point of the ith element,

[mathematical expression not reproducible]. (25)

For x at the ith element section, where x is a global coordinate,

[mathematical expression not reproducible]. (26)

Suppose

[mathematical expression not reproducible]. (27)

Then, (20) can be expressed as

[N.sub.i] (x) = [A.sub.i](x) x [N.sub.1] ([0.sub.1]). (28)

Thus, the relation between the vector [N.sub.n] ([x.sub.n]) at the point [x.sub.n] and the vector [N.sub.1] ([0.sub.1]) at the start point [x.sub.0] is

[N.sub.n]([x.sub.n]) = [A.sub.n]([x.sub.n]) x [N.sub.1] ([0.sub.1]), (29)

where

[mathematical expression not reproducible]. (30)

The boundary condition of the existing tunnel assumed as an infinite long tunnel considering the shear stress of soil can be written as

[mathematical expression not reproducible]. (31)

By applying the boundary conditions to (19), the deflection, rotational angle, bending moment, and shear force of the overall beam can be obtained.

6.4. Case Study

6.4.1. Parameter Selection. To enable a direct comparison, assumed parameters of the soil and existing tunnel are selected in this study. Poisson's ratio of the soil is 0.388. The grouting process is regarded as the inverse process of soil loss due to tunnelling, and g can be expressed as g = (2[DELTA] + [delta])(1 - [eta]), where [eta] is the injection ratio.

For the Pasternak foundation, as the thickness of the elastic layer is assumed to be 2.5 B, the subgrade elastic modulus k and shear stiffness G are suggested by Xu [19] as follows:

[mathematical expression not reproducible], (32)

where [beta] is a modified parameter (3 = 1.29) and [E.sub.0] and [v.sub.0] are, respectively, the elastic modulus and Poisson's ratio for elastic foundation, which are given as follows:

[E.sub.0] = [E.sub.s]/1-[v.sup.2.sub.s], [v.sub.0] = [v.sub.s]/1-[v.sub.s], (33)

where the empirical relation between the Young's modulus [E.sub.s] and the compression modulus [E.sub.s0.1-0.2] of soft soils is built up by Yang and Zhao [20].

[E.sub.0] = (2.5 ~ 3.5)[E.sub.S0.1-0.2]. (34)

As illustrated in Figure 14, points [O.sub.1] and [O.sub.2] are, respectively, the intersection points A and B; in practical, coordinate transformation expressions of the axis x and x" are, respectively,

x' = (x - [L.sub.1])sin [alpha], x" = (x - [L.sub.2])sin [alpha], (35)

where [L.sub.1] and [L.sub.2] are the distance from the beginning point of the existing tunnel to points [O.sub.1] and [O.sub.2], respectively. The length of the existing tunnel is assumed as 250 m, and the ring spacing is 1.2 m. So, the tunnel is divided into 209 sections, and each section is 1.2 m in length.

7. Discussions

The normalized bending moment and shear force are defined as

[mathematical expression not reproducible], (36)

where [[??].sub.i] and [Q.sub.i] are, respectively, the normalized bending moment and shear force; [M.sub.i] and [Q.sub.i] are, respectively, the calculated bending moment and shear force; and [L.sub.0] is the length of the existing tunnel.

Figure 19 shows the comparison of the deflection, rotational angle, bending moment, and shear force of the existing tunnel with a constant stiffness which is based on three kinds of models: an infinite tunnel on the Winkler foundation, an infinite tunnel on the Pasternak foundation, and a finite tunnel on the Winkler foundation. Note that the existing tunnel is only subject to soil deformation force in this case. The displacement and the deflection of beam are particular similar concepts owing to the fact that the existing tunnel is securely attached to the foundation, generally without rigid body motion. Fairly good agreement can be observed between the infinite tunnel on the Pasternak foundation and the finite tunnel on the Pasternak foundation, which shows that the method proposed in this paper is correct and can be used to analyze similar problems. Because of considering soil shear stress, vertical stress and deformation are affected and weakened; the displacement of the existing tunnel on the Pasternak foundation is smaller than that on the Winkler foundation without considering the shear stress of soil. Using the Winkler foundation model is in good agreement with the observed data. It can be seen from (3) that the displacement of the existing tunnel is proportional to the grouting efficiency (i.e., [eta]). Due to the uncertainty of [eta], it cannot be decided which one is the best analytical method. It is clear that the maximum normalized bending moment occurs at the points [O.sub.1] and [O.sub.2], which are the nearest points between the existing tunnel and the new tunnels.

The reinforce range of the existing tunnel is assumed as 25 rings on either side of the intersection point A and the various stiffness is taken into consideration using the method proposed in this paper. Figure 20 shows the influence of different stiffness of the reinforcement on the existing tunnel (i.e., 4EI, EI, and 1/4EI). It is evident that the reinforcement can reduce the displacement, rotational angle, normalized bending moment, and shear force of the existing tunnel, while tunnel damage (reduce the stiffness of reinforcement) has the opposite effect. It is evident that the change of reinforcement stiffness has greater effect on the normalized bending moment and the normalized shear force of the existing tunnel, particularly at the inflection points of the normalized bending moment and normalized shear force curve.

The reinforce stiffness of the reinforcement range is assumed as 4 EI and the various reinforcement ranges is taken into consideration using the method proposed in this paper. Figure 21 shows the influence of different reinforcement ranges on the existing tunnel when the northbound tunnel of Line 4 underpasses. The normalized bending moment on either side of the intersection point is too large when the reinforcement range is 69-89.2 m, and the normalized shear force on either side of the intersection point is relatively too large when the reinforcement range is less than 59.2-99.2 m. So, it is necessary to reinforce larger than the range of 49.2-109.2 m in the existing tunnel in this case, that is, larger than 30 m (25 rings) on either side of the intersection point. There is no big difference in the bending moment and shear force of the existing tunnel when the reinforcement ranges are 49.2-109.2 m, 39.2-119.2 m, and 29.2-129.2 m. In order to reduce the costs of reinforcement, reinforcement in the 30 m range (i.e., 25 rings) on either side of the intersection point is the best choice which verifies the actual design.

8. Conclusions

Based on the abovementioned statements, the response of the existing tunnel is analyzed, and analytical methods are proposed for verification and the reinforcement design. Main conclusions derived from the analysis are as follows:

(1) The displacement of the existing tunnel changes with the relative position of the shield and the existing tunnel. Shield parameters have a major impact on the existing tunnel. The heaves of the existing tunnel might be caused by the large applied face pressure and bulkhead addictive thrust. When the shield tail is driving beyond the selected monitoring rings, the selected monitoring rings settle rapidly due to the closure of the shield tail void.

(2) The southbound tunnel of Line 4 has less effect on the left side of the intersection point A on the existing tunnel because of the northbound tunnel's barrier effect. A variety of measures of displacement prevention of Line 1 such as the control of shield parameters, reinforcement of the existing tunnel can ensure the normal operation of metro Line 1.

(3) The analytical method proposed in the paper is verified by the methods derived from prior studies. It is possible to take the different stiffness and curved loads into account by applying the stepped stiffness and trapezoidal load while using a local coordinate system in the derivation of the new method. The change of reinforcement stiffness has greater effect on the normalized bending moment and the normalized shear force of the existing tunnel, particularly at the inflection points of the normalized bending moment and normalized shear force curve. The analytical method can be applied for the reinforcement design and safety check of the built tunnel.

Appendices

A.1. An Infinite Beam on the Winkler Foundation

[mathematical expression not reproducible], (A.1)

[mathematical expression not reproducible], (A.2)

[mathematical expression not reproducible]. (A.3)

A.2. An Infinite Beam on the Pasternak Foundation

[mathematical expression not reproducible], (A.4)

[mathematical expression not reproducible]. (A.5)

B.1. Coefficients of the Matrix [k.sub.i5]

[mathematical expression not reproducible], (B.1)

[mathematical expression not reproducible], (B.2)

[mathematical expression not reproducible], (B.3)

[mathematical expression not reproducible]. (B.4)

https://doi.org/ 10.1155/2018/1598672

Conflicts of Interest

The author declares that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

The author appreciates the help from the staff of Hangzhou and Hong Kong Tunnel Company Ltd. during instrumentation setup and data acquisition. This study is supported by the National Natural Science Foundation of China (no. 41702313).

References

[1] S. H. Kim, Model testing and analysis of interactions between tunnels in clay [Ph.D. Thesis], Department of Engineering Science, University of Oxford, Oxford, UK, 1996.

[2] H. Y. Liu, J. C. Small, J. P. Carter, and D. J. Williams, "Effects of tunnelling on existing support systems of perpendicularly passing tunnels," Computers and Geotechnics, vol. 36, no. 5, pp. 880-894, 2009.

[3] I. Yamaguchi, I. Yamazaki, and Y. Kiritani, "Study of ground-tunnel interactions of four shield tunnels driven in close proximity, in relation to design and construction of parallel shield tunnels," Tunnelling and Underground Space Technology, vol. 13, no. 3, pp. 289-304, 1998.

[4] T. I. Addenbrooke and D. M. Potts, "Twin tunnel interaction: surface and subsurface effects," International Journal of Geomechanics, vol. 1, no. 2, pp. 249-271, 2001.

[5] J. P. Kimmance, S. Lawrence, O. Hassan, and N. J. Purchase, "Observations of deformations created in existing tunnels by adjacent and cross cutting excavations," in Proceedings of Geotechnical Aspects of Underground Construction in Soft Ground, pp. 707-712, London, UK, July 1996.

[6] E. Leca and B. New, "Settlements induced by tunnelling in soft ground," Tunnelling and Underground Space Technology, vol. 22, no. 2, pp. 119-149, 2007.

[7] S. M. Liao, J. H. Liu, R. L. Wang, and Z. M. Li, "Shield tunnelling and environment protection in Shanghai soft ground," Tunnelling and Underground Space Technology, vol. 24, no. 4, pp. 454-465, 2009.

[8] F. H. Chehade and I. Shahrour, "Numerical analysis of the interaction between twin-tunnels: influence of the relative position and construction procedure," Tunnelling and Underground Space Technology, vol. 23, no. 2, pp. 210-214, 2008.

[9] X. G. Li and D. J. Yuan, "Response of a double-decked metro tunnel to shield driving of twin closely underpasses tunnels," Tunnelling and Underground Space Technology, vol. 28, no. 1, pp. 18-30, 2012.

[10] T. E. Vorster, A. Klar, K. Soga, and R. J. Mair, "Estimating the effects of tunnelling on existing pipelines," Journal of Geotechnical and Geoenvironmental Engineering, vol. 131, no. 11, pp. 1399-1410, 2005.

[11] Z. G. Zhang, M. S. Huang, and W. D. Wang, "Evaluation of deformation response for adjacent tunnels due to soil unloading in excavation engineering," Tunnelling and Underground Space Technology, vol. 38, no. 3, pp. 244-253, 2013.

[12] A. P. S. Selvadurai, Elastic Analysis of Soil-Foundation Interaction, Elsevier Scientific Pub. Co., Amsterdam, Netherlands, 1979.

[13] N. Loganathan and H. G. Poulos, "Analytical prediction for tunnelling-induced ground movements in clays," Journal of Geotechnical and Geoenvironmental Engineering, vol. 124, no. 9, pp. 846-856, 1998.

[14] K. M. Lee, R. K. Rowe, and K. Y. Lo, "Subsidence owing to tunnelling. I. Estimating the gap parameter," Canadian Geotechnical Journal, vol. 29, no. 6, pp. 929-940, 1992.

[15] R. K. Rowe and K. M. Lee, "Subsidence owing to tunnelling. II. Evaluation of a prediction technique," Canadian Geotechnical Journal, vol. 29, no. 6, pp. 941-954, 1992.

[16] P. B. Attewell, J. Yeates, and A. R. Selby, Soil Movements Induced by Tunnelling and their Effects on Pipelines and Structures, Methuen Inc., New York, NY, USA, 1986.

[17] F. Ye, C. He, H. H. Zhu et al., "Analysis of longitudinal equivalent stiffness of shield tunnel considering lateral performance," Chinese Journal of Geotechnical Engineering, vol. 33, no. 12, pp. 1870-1876, 2011, in Chinese.

[18] P. R. Lancaster and D. Mitchell, "Problems in Bending of Bars and Plates [M]," Advanced Solid Mechanics, Macmillan, UK, 1980.

[19] L. Xu, "Study on vertical settlement of shield tunnel in soft soil," Ph.D. dissertation, Tongji University, Shanghai, China, 2005.

[20] M. Yang and X. Zhao, "An approach for a single pile in layered soil," Tongji University. Natural Science, vol. 20, no. 4, pp. 421-428, 1992, in Chinese.

Qiongfang Zhang [ID] (1,2)

(1) MOE Key Laboratory of Soft Soils and Geoenvironmental Engineering, Zhejiang University, Hangzhou 310058, China

(2) Research Center of Coastal and Urban Geotechnical Engineering, Zhejiang University, Hangzhou 310058, China

Correspondence should be addressed to Qiongfang Zhang; yangziduozi@sina.com

Received 10 August 2017; Accepted 16 November 2017; Published 20 March 2018

Academic Editor: Andrea Benedetto

Caption: FIGURE 1: Location of the new tunnels and the existing tunnel in Hangzhou. A: Guanhe Station; B: East railway station.

Caption: FIGURE 2: Plan view of tunnel alignment and location of the monitored rings.

Caption: FIGURE 3: Longitudinal profile of soils and tunnels.

Caption: FIGURE 4: Soil profile and geotechnical parameters.

Caption: FIGURE 5: Advance-time curve of the shield.

Caption: FIGURE 6: The applied tunnel face pressure.

Caption: FIGURE 7: The grouting volume.

Caption: FIGURE 8: Arrangement of the monitoring points at the cross section of the existing tunnel.

Caption: FIGURE 9: Time-varying vertical displacement of the existing tunnel.

Caption: FIGURE 10: Time-varying horizontal displacement of the existing tunnel.

Caption: FIGURE 11: Convergence displacement of the existing tunnel.

Caption: FIGURE 12: Reinforcement in the existing tunnel.

Caption: FIGURE 13: A schematic view of the existing tunnel-soil tunnelling interaction.

Caption: FIGURE 14: A schematic diagram of tunnels for the analytical methods.

Caption: FIGURE 15: Sketch of equivalent stepped stiffness.

Caption: FIGURE 16: Forces and local coordinate system.

Caption: FIGURE 17: Replacement of variable load with trapezoid load.

Caption: FIGURE 18: The matrix transfer diagram.

Caption: FIGURE 19: (a) Comparison of the deflection of the existing tunnel using different models. (b) Comparison of the rotational angle of the existing tunnel using different models. (c) Comparison of the normalized bending moment of the existing tunnel using different models. (d) Comparison of the normalized shear force of the existing tunnel using different models.

Caption: FIGURE 20: (a) The effects of different reinforcement stiffness on the deflection of the existing tunnel. (b) The effects of different reinforcement stiffness on the slope of the existing tunnel. (c) The effects of different reinforcement stiffness on the normalized moment of the existing tunnel. (d) The effects of different reinforcement stiffness on the normalized shear force of the existing tunnel.

Caption: FIGURE 21: (a) The effects of different reinforcement ranges on the deflection of the existing tunnel. (b) The effects of different reinforcement ranges on the slope of the existing tunnel. (c) The effects of different reinforcement ranges on the normalized moment of the existing tunnel. (d) The effects of different reinforcement ranges on the normalized shear force of the existing tunnel.

TABLE 1: Laboratory penetration test and field steady flow test. Laboratory Field steady flow Soil penetration test (cm/s) test (cm/s) layer [K.sub.v] [K.sub.h] K 3-2 4.79 x [10.sup.-5] 7.81 x [10.sup.-5] 4.21 x [10.sup.-3] 3-3 2.55 x [10.sup.-4] 2.30 x [10.sup.-4] 1.57 x [10.sup.-3] 3-5 9.63 x [10.sup.-5] 1.74 x [10.sup.-4] -- 3-6 3.29 x [10.sup.-4] 3.01 x [10.sup.-4] -- 4-3 2.39 x [10.sup.-7] 7.84 x [10.sup.-7] / 6-1 2.06 x [10.sup.-7] 3.03 x [10.sup.-7] -- 6-2 2.62 x [10.sup.-7] 1.28 x [10.sup.-7] -- 8-1 1.88 x [10.sup.-7] 5.11 x [10.sup.-7] --

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Zhang, Qiongfang |

Publication: | Advances in Civil Engineering |

Article Type: | Report |

Geographic Code: | 1USA |

Date: | Jan 1, 2018 |

Words: | 7011 |

Previous Article: | Formation Mechanism and Mechanical Properties of Soil-Rock Mixture Containing Macropore. |

Next Article: | Surficial Failure of Expansive Soil Cutting Slope and Its Flexible Support Treatment Technology. |

Topics: |