Abstract

Prediction of tunneling induced soil deformation is an important issue in the design of tunnels constructed in the densely populated urban areas. In this paper, commercial FEM software 3D ABAQUS is adopted to simulate the behavior of soil caused by tunneling under vehicle loads in the coastal reclamation area. A field case study was also carried out to verify the accuracy of the proposed model. A good agreement was achieved. It is also found from the studies that the areas affected by soil deformation can be classified into four zones: the key disturbed zone, the secondary disturbed zone, the general disturbed zone, and minor-disturbed zone. The maximum soil deformation occurs on side of the longitudinal medial axis of the tunnel. The shape of the settlement curve is almost the same as that of the normal distribution curve. The soil deformation in the action zone of vehicle load is greater than that of the nonaction zone of vehicle load.

1. Introduction

The development of the subway system is one of the primary approaches to solve the problem of lacking of space and resource, as well as improving the efficiency of urban land usage and urban center density mitigation. With the increase of technology development on boring machine, tunneling has become one of the major methods for subway construction. During urban tunneling, the tunnel is often required to be constructed under existing structures, such as expressway, railway, sewerage, power transmission tunnels, communication, and other subsurface lifelines [1]. In such situations, soil deformation caused by tunneling may inevitably affect the adjacent structures by generating additional displacements and stresses [2, 3].

Over the last few decades, many approaches, that is, empirical formulas, theoretical analyses, model tests, artificial neural networks and expert intelligent forecasting systems, and numerical simulations, have been performed to investigate the soil deformation with tunneling. Among these methods, the empirical formula method [47] features a clear concept, simple calculation, and ease of use. To a certain extent, it reflects the effect of ground settlement on geological conditions, such as formation, tunnel depth, and excavation diameter. However, the empirical formula method cannot accurately forecast the range of ground settlement as it can only take into account few soil parameters and simple geological conditions. The theoretical analysis method [811] has the advantages of low computational complexity and high precision, but it is hard to determine precisely the boundary conditions in a complex 3D space especially for geotechnical engineering applications. The model test method [12, 13] could more accurately model the influence of the construction factors on ground settlement. Nevertheless, it is costly, is complex, and suffers from poor controllability. The construction parameters are used merely as a reference for theoretical analysis and field monitoring. The artificial neural network and expert intelligent forecasting system [1419] have high nonlinear mapping capability, but they need substantial large amount of in situ monitoring data and errors in knowledge acquisition, parallel inference, adaptive learning, associative inference, and easily occurred fault tolerance.

The numerical simulations [2026] have the advantage of considering different influencing parameters and boundary conditions. For the simulation of the behavior of soil caused by tunneling in the coastal reclamation area, most of the researches in art are based on 2D simulations and few of them on full 3D modeling [1]. Furthermore, the studies on tunneling analysis were focusing on the assessment of the ground surface settlement except little of them which start to pay attention to the interaction between tunneling and existing adjacent structures such as bridges, pile-supported buildings, twin tunnels. Even little research work has been carried out on the interaction between shield tunnels under vehicle loads and soil deformation in a coastal reclamation area.

In this paper, a 3D numerical investigation on soil deformation caused by tunneling under vehicle loads in the coastal reclamation area is carried out. The study aims to identify the fundamental deformation mechanisms involved in the responses of the tunneling under vehicle loads, in terms of the soil deformation and types of vehicles. A field case study was also carried out to verify the accuracy of the proposed model.

2. Description of the Project

2.1. Characteristics of the Tunnel

A typical tunnel section was adopted for the Shenzhen Metro Line 2 construction project. This tunnel was constructed from Houhai Station to Keyuan Station with the earth pressure balance (EPB) shield machine. The EPB shield machine has an external diameter of 6.28 m and a length of 7.58 m. It was launched from the departure shaft at Houhai Station and finally reached the arrival shaft at Keyuan Station (Figure 1). The tunnel of this section is 2473.708 m long, with its mileage from YDK9 + 629.740 to YDK10 + 850.200. At each section, six prefabricated concrete sections were placed to form the tunnel lining and joined with bolts in the longitudinal and circumferential direction. The size of each ring of lining is 0.30 m thick and 1.5 m wide. The concrete tunnel lining has Young’s modulus of 34 500 MPa and Poisson’s ratio of 0.167.

2.2. Characteristics of the Soil

A typical ground condition encountered in the Shenzhen metropolitan area was also considered in this study. Shenzhen lies in the Pearl River Delta plain and has a complex geology of a mostly watery and weak stratum. The soil profile of the city is characterized by high sensitivity, high water content, and high changeability. Moreover, as numerous coastal reclamation projects have been carried out in Shenzhen recently, the city is now confronting with special engineering problems that differ from those in Guangzhou and Shanghai during subway construction in soft soil. Such problems include the following: (1) the soil layer is still under consolidation as the coastal reclamation project was just finished in 2002; (2) the groundwater table is affected by the flood and ebb tides obviously due to the strong hydraulic connection between groundwater and seawater; and (3) there is a great amount of rubble (gravel) with large diameters in the back fill.

A typical geotechnical section is given in Figure 2. The site was reclaimed from the sea by a soft foundation treatment in 2002. In this area, the top layers mainly consist of backfill and muddy clay. A shallow backfill layer with a thickness of 6.0 m is present on top of the muddy clay layer which has a depth of up to 4.0 m. Under the muddy clay layer, there is a thin layer of clay. The layers under this second clay layer are the gravelly sand layer and the weathered rock layer, and the tunnel is mainly embedded in the gravelly sand layer. The groundwater of the working site can be divided into two types: loose soil pore water and bedrock fissure water. The main occurrence of pore water is consistent with the Quaternary System Holocene series marine layer, which includes gravel sand (containing mud) and a proluvium-alluvium gravel sand layer. The bedrock fissure water originates primarily from the bedrock fracture of strongly weathered rock to medium-weathered rock. The two types of groundwater generally have no pressure but similar penetrations. The groundwater table is in the range of 0.70–64.60 m, and the maximum tunnel inflow rate is 7576 m3/d. The site belongs to a very rich water area; water inflow and local quicksand are common occurrences. Therefore, the soft soil in the coastal reclamation area is featured by bad water permeability and big water content. To simplify the calculation, our analysis did not consider the impact of the consolidation and creep of clay.

2.3. Characteristics of the Expressway and Vehicles

The total length of the Binhai Expressway within the territory of Shenzhen is 9.8 km, consisting of 610000 m2 subgrade reclamation from Nanshan District to Luohu District. The Binhai Expressway is an eight-lane urban expressway with heavy traffic and frequent congestion in the morning and evening peak hours. The traffic condition can be assessed in terms of annual average daily traffic (AADT). The AADT is more than 600,000 vehicles per day for each traffic direction. The shield tunnel is from ZDK9 892.500 to ZDK10 45.000 mileage beneath Binhai Expressway and passes through the expressway twice (Figure 1). On this section of the expressway, traffic during the evening peak hours is highly congested, and the vehicles move towards the Binhai Expressway with an average speed of 22.12 km/h, which is much lower than the design speed limit of 80 km/h. The pavement structure sections [27] are illustrated in Table 1. As it can be seen, the pavement structure from the surface to the base can be divided into four layers in terms of its material type.

3. 3D Finite Element Model (FEM) and Methodology

3.1. Numerical Tool and Constitutive Models

In the current study, the finite element (FE) analysis package ABAQUS [28] was used in the numerical analyses to study the effects of soil deformation caused by tunneling under vehicle loads in the coastal reclamation area. The soil was simulated using an elastoplastic constitutive relationship based on the Mohr-Coulomb criterion and the tunnel lining was modeled as an isotropic linear elastic material.

3.2. Calculation Parameters

Table 2 shows the calculation parameters. The model parameters of soil layers in Table 2 are obtained based on geotechnical drilling, in situ data, and laboratory testing of Shenzhen Metro projects [2932]. The moisture content of the soil has been determined through oven drying method. The compression modulus is obtained from standard penetration test, static penetration test, and pressuremeter test. The internal friction and cohesion of the soil have been determined through standard penetration test, dynamic penetration test, and vane shear test. The Poisson’s ratio used in the current study was slightly modified from a previous study by Chen et al. [29, 30, 32] and Wang et al. [31]. The water pressures and initial stresses are calculated separately. The burial depth of groundwater is 0.70–4.60 m, and the water level elevation is −1.40 to 2.81 m. The groundwater of the working site mainly is loose soil pore water. The site belongs to a very rich water area. This paper adopts the method of joint calculation of water and earth pressures. Soil bulk density below the groundwater is calculated using the effective unit weight, which is 10 kN/m3 in this numerical analysis. The initial stresses are calculated by means of the procedure. The value of is based on Jaky’s formula: , in which is the effective angle of the internal friction of the host soil [33]. However, it should be noted that the internal friction angles of gravelly sand and rock are too low (20~24°) for six reasons: mineral composition, grain gradation, original density, water content, structure, and precharge pressure. Because this project is located in coastal reclamation area, the internal friction angle is lower than normal consolidation soil.

3.3. 3D Finite Element Model (FEM)

Considering the symmetry of the geometry and loading conditions, only half of the tunnel (Figure 3(a)) is represented in the 3D FEM model. Figure 3(b) shows the soil profile in this study. All relevant components of the construction process, including the treatment of the soil and ground water, tunneling, and tail void grouting, are taken into account as separate components in the model. The effect of tail void injection is modeled by a pressure diagram applied on a 5-meter distance behind the tail [34]. A partial filling of the void due to earth pressure balance shield machine (EPBS) diameter reduction is assumed by taking into account only 25% of the actual diameter reduction. The calculation range of the model is 36 m × 45 m × 33 m along the , , and axes, respectively. The positive direction of the axis is the direction of the shield tunneling, while the axis is in the upward positive and the axis is in the upward right positive. In this model, the average depth of the tunnel is 16 m (measured from ground surface to tunnel crown) and the excavation diameter of tunnel is 6.28 m. As indicated by displacement boundary conditions, the upper boundary of the model is considered as a free boundary while the bottom of the tunnel is considered as a fixed boundary of which the horizontal and vertical movements are limited. The shield machine is modeled by a plate element, and the solid elements are assigned for lining with an elastic behavior. All the analyses for tunnel excavation are based on the undrained condition. The finite element mesh includes 82524 eight-node isoparametric hexahedral elements with 88874 nodes. In our study, the goal is to compare the calculated values obtained using 3D calculations to the experimental ones. So, for the sake of numerical modeling, four critical measurement points (monitoring points 1, 2, 3, and 4) and two critical monitoring sections (monitoring section 1 and monitoring section 2) in the tunnel were set (Figure 3(a)). They are located at the ground surface. The monitoring section 1 is located 15 m from the beginning of the axes. The monitoring section 2 is located 30 m from the beginning of the axes. Monitoring points 1, 2, 3, and 4 are located at the center of tunnel (resp., 7.5 m, 15.0 m, 22.5 m, and 30.0 m from the beginning of the axes). Note that the stress path analysis was performed for each incremental advance of the tunnel face using an iterative solver, as reported by Chen et al. [32]. Therefore, in the numerical analysis, the influence of stress path was ignored. For more on the influence of stress path on tunnel excavation response, please refer to [32, 35].

3.4. Vehicle Loads

According to vehicle loads investigations, the vehicle loads of 200 kN and 550 kN have the most significant influence on the soil deformation on Binhai Expressway. Based on the characteristic of vehicle loads and the research results [3638], vehicle loads are simplified as quasi-static loads which distributed centrally at the places of the axles. The value of the load applied at each axle follows the standard of technical standard of highway engineering (JTG-2003) [39] (Figures 4(a) and 4(b) and Table 3). At the same time, the force of wheels acting on a road is also assumed to be concentrated force (Figures 5(a) and 5(b)). The vehicle loads of 200 kN and 550 kN are adopted for calculation under extreme condition with symmetrical dynamic load; that is, four vehicles act on a calculation section (Figure 6). The rear axle load is adopted to calculate the size of the vehicle loads according to the literatures [37, 38].

3.5. Calculation Conditions

To analyze the effect of shield construction on soil deformation in the study site comprehensively while taking the disturbance from long-term cyclic dynamic loading of vehicles on Binhai Expressway into consideration, it is necessary to categorize the tunnel space, buried depth, construction parameters, time, and other factors into two stages for calculation: (1) stage I: the vehicle load is 200 kN during the shield beneath Binhai Expressway; and (2) stage II: the vehicle load is 550 kN during the shield beneath Binhai Expressway. The simulation and calculation of each vertical excavation were implemented with an interval of 1-segment rings for a group, after which comparative analysis of the two stages was conducted.

3.6. Excavation Construction Procedure

The numerical simulation has been carried out by modeling the excavation of 45 m of tunnel, following a model change method (1.5 m per step) [28]. The excavation and installation of the liner are modeled in four analysis steps; that is, the initial stress state is applied and the liner elements are removed, the tunnel excavation begins by reducing the concentrated loads on the tunnel surface, the loads are reduced by 25% in this step, and after that the liner is installed, and then the surface load is applied, and the excavation is completed by removing the remainder of the load on the tunnel perimeter. Such a simplified modeling approach adopted in the current study was based on a previous study by Yoo [2], Li et al. [40], Wong et al. [41], and Wang et al. [31]. The tunneling process, consisting of a series of excavation and shotcrete lining installation stages, was closely followed by removing and adding corresponding elements at designated steps assuming full face excavation with an advance length of 1.5 m. In the numerical analysis, a constant pressure which equals the at-rest earth pressure at the centre of tunnel face was assumed following Wong et al. [41] and Yoo [2]. The tunnel excavation was then simulated in 30 steps.

4. Results and Discussions

Four critical measurement points and two critical monitoring sections in the tunnel were set: shield tunneling at 5th, 10th, 15th, and 20th rings (Figure 3(a)). To illustrate soil deformation, we take monitoring section 1 as an example.

4.1. Analysis of Transverse Settlement

The transverse settlement distances curves of stages I and II in the shield beneath Binhai Expressway are shown in Figures 7 and 8 and Table 4. The maximum soil deformation occurs near the longitudinal medial axis of the tunnel, whereas the minimum cumulative settlement occurs at a distance away from the farthest point of the midline. The settlement is mainly distributed within the area of four times of the tunnel diameter, and the shape of the settlement curve is highly similar to the normal distribution curve. The results indicate that the closer the distance from tunnel face, the greater the disturbance from vehicle load and tunnel excavation to the soil deformation of stratum.

The areas affected by soil deformation can be divided into four disturbed zones, namely, the key disturbed zone, the secondary disturbed zone, the general disturbed zone, and the minor-disturbed zone. The definitions of the key disturbed zone, the secondary disturbed zone, the general disturbed zone, are the minor-disturbed zone are the property of the soft soil which allows it to produce deformation and breakage within 4 meters, 4 to 10 meters, 10 to 25 meters, and 25 to 36 meters from the tunnel centerline (settlement trough width), respectively, under the interaction of shield tunneling and vehicle loads. The locations of these four zones are 4 m, 4–10 m, 10–25 m, and 25–36 m away from the lateral midline of the tunnel, respectively. Among the four zones, the key disturbed zone, which exhibits disturbance of soil above the tunnel, is affected by shield construction most remarkably. Table 4 shows the soil deformation of the shield beneath Binhai Expressway, in which the vehicle load 550 kN exhibits more substantial deformation than the 200 kN does. The maximum soil deformations are −16.21 mm, −6.36 mm, +1.72 mm, and +1.07 mm for the key disturbed zone, secondary disturbed zone, general disturbed zone, and minor-disturbed zone, respectively. The soil layer in the coastal reclamation area is loose. The shield beneath Binhai Expressway passes mainly through a sandy (gravel) cohesive soil layer which has good permeability and abundant water. As the layer has unevenly thickness and poor self-stability, it is easy to lose its strength when affected by groundwater. During the initial shield construction, time, excessively large established earth pressure, grouting volume, grouting pressure, long-term loading of ground vehicles, and other factors will impose a considerable effect on original formation and stress field and cause the stress redistribution across the tunnel. Therefore, strict controlling of the earth pressure, driving speed, grouting pressure, grouting volume, and tunneling direction is recommended. The tunnel excavation shows that the vehicle loads have influences on soil deformation, as reported by Wang et al. [31]. The original structure and the state of limit equilibrium of soil are completely destructed during tunneling. The stress distribution around the tunnel was changed due to tunnel construction. But the soft soil in the coastal reclamation area is characterized by high sensitivity, high water content, and high changeability. Tunneling in the coastal reclamation area is very difficult, because the soft soil’ self-stability is very poor. The uplifts of the ground surface have occurred due to the effect of excavation unloading, the friction between the shield shell and soil mass, and the effective synchronous grouting. However, the flexural deflection of road has occurred under vehicle loads, and it is a tensile stress upwards component, with the properties uploading, counteracting the effects of partial vehicle load. Then, the downwards components of pressure stress with the properties loading have prevented the uplifts of the soil. The effect of vehicle loading is more significant than that of excavation which is affected synthetically by loading, uploading of vehicle loads, and tunnel excavation.

4.2. Analysis of Longitudinal Settlement

The longitudinal settlements of stages I and II are shown in Figures 9 and 10. As it can be seen, the soil deformation in the action zone of vehicle load is greater than that of the nonaction zone of vehicle load. With the increase of vehicle loads, the soil deformation increases gradually. The maximum soil deformations of stage I are −6.52 mm, −6.93 mm, −7.5 mm, and −7.6 mm for the shield tunneling at 5th, 10th, 15th, 20th rings, respectively. The maximum soil deformations of stage II are −14.81 mm, −15.34 mm, −16.05 mm, and −16.21 mm for the shield tunneling at 5th, 10th, 15th, and 20th rings, respectively.

5. Comparison between the Numerical Predictions and Experimental Results

5.1. Shield Construction Monitoring

To verify the accuracy of the finite element simulation results, two-point settlement observation sections were placed in the tunnel. In this case, we established a shield tunneling to 10th, 20th ring and selected the shield tunneling to 10th ring sections from the five observation sections for monitoring (Figure 11). A settlement measurement point and a cross-section measurement point were designed along the longitudinal line of the ground every 5 m and 30 m along the advancing direction of the shield, respectively.

5.2. Comparative Analysis of Transverse Deformation

The curves of the soil deformation durations (shield beneath Binhai Expressway) for stages I and II are shown in Figure 12. The comparative analysis indicates that the following.(1)The soil deformation of the shield beneath Binhai Expressway in stage II is larger than that observed in stage I, indicating the complex geology of the coastal reclamation area. The calculation conditions of construction, changes of the original formation and stress filed, and excavation stress impose a strong effect on soil deformation. The cumulative maximum soil deformation occurs at the longitudinal medial axis of the tunnel, whereas the minimum cumulative settlement is observed at a distance away from the farthest point of the midline. Compared with the observed data between stage I and stage II, the simulated results fit well with the field observation. At monitoring section 1, the maximum simulated cumulative settlement for the stage II is −15.34 mm while the maximum and the minimum monitoring cumulative settlements are −12.68 mm and −0.29 mm, respectively. This is because the tunnel width is small, and the formation is characterized primarily by sandy (gravel) cohesive soil ⑧2, ⑧3. The soil exhibits plasticity to low plasticity and comprises 10%–30% gravel volume. A large boulder group formed in the prophase sea reclamation exists on top of the tunnel. The rubble has a large diameter (Figure 13), which will possibly reduce driving speed, slow down cutter rotational speed, and increase thrust in shield construction. Therefore, the comprehensive analysis for controlling soil deformation should consider prediction and monitoring data.(2)Ground settlement is distributed in the diameter zone with a fourfold diameter. The centerline of this zone is the symmetrical axis. The key ground settlement occurs above the tunnel, which is the secondary ground settlement zone, with a 0.5–2.6 times diameter, on the centerline of the tunnel. The minor-ground settlement zone is beyond the 4 times diameter and can be disregarded. This conclusion is consistent with the prediction result on soil deformation.

5.3. Comparative Analysis of Longitudinal Deformation

Figure 14 shows the cumulative settlements of the longitudinal of stages I and II in the shield beneath Binhai Expressway. The comparative analysis indicates the following.(1)The overall control of the surface settlement caused by the tunnel is good, indicating that controlling soil deformation on the basis of the FE analysis and site monitoring is a reasonable approach. The maximum cumulative surface settlement of stage I is −6.93 mm, located at 14.3 m away from the ring number 0, whereas the minimum settlement is −0.12 mm, located at 45 m way from the ring number 0. The surface settlement of stage II is larger than that of stage I. The ground settlement reaches −15.34 mm, especially before the shield beneath Binhai Expressway. However, this issue is effectively controlled after using a reasonable control technology. The average settlement of the action zone of vehicle loads is −13.21 mm. These results can be used to forecast the soil deformation and accordingly a measurement for shield construction control can be proposed. In this measurement, an accurate earth pressure is selected, driving speed, mixed proportions of grouting material, and the amount of oil in the shield tail are controlled, monitoring and measurement are strengthened, and information management is implemented. The overall control of surface settlement ensures construction security and tunnel stability. It also reflects the importance of the prediction model ABAQUS.(2)Disturbance from the long-term cyclic dynamic loading of vehicles on Binhai Expressway and the Binhai Expressway overpass near the tunnel also influence soil deformation. The mileage of 9.3 to 20.5 m is located below the main carriageway of Binhai Expressway. The cumulative settlement is large; the maximum is −15.34 mm; and the minimum is −8.46 mm. This result indicates the positive effect of controlling surface settlement.

6. Conclusions

In this paper, 3D ABAQUS is adopted to simulate the soil deformation due to tunneling under vehicle loads in the coastal reclamation area. A field case study was also carried out to verify the FEM model. A good agreement between the numerical results and the experimental results was achieved. Other main conclusions derived from the study are summarized as follows.(1)The areas affected by soil deformation can be classified into the key disturbed zone, the secondary disturbed zone, the general disturbed zone, and the minor-disturbed zone. Among the four zones, the key disturbed zone that exhibits disturbance of the soil above the tunnel is affected most remarkably by the soil deformation under tunneling. The four zones are located at 4 m, 4–10 m, 10–25 m, and 25–36 m away from the midline of tunnel along the lateral direction, respectively.(2)The maximum soil deformation occurs near the longitudinal medial axis of the tunnel, whereas the minimum cumulative settlement occurs at the farthest point of the midline. The settlement distributed in a range of four times of the tunnel diameter and the shape of the settlement curve is highly similar to the normal distribution curve.(3)The soil deformation in the action zone of vehicle load is greater than that of the nonaction zone of vehicle load. The closer the distance from the tunnel face, the greater the disturbance from vehicle load and tunnel excavation to the soil deformation.

In the numerical analysis, the interaction of vehicle loads, road, soil, and tunnel was not considered. In addition, the influence of strut installation was ignored to simplify the problem. In the future, the effects of additional factors and other affecting parameters, such as vehicle types, vehicle frequencies, the distribution pattern of the vehicle loads, and the rheological behavior of the soft soil on the distribution and the magnitude of soil deformation under shield construction, are suggested to be studied through numerical methods, theoretical analysis, and the physical models.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The financial supports for this research project by the National Natural Science Foundation of China (no. 41072205), the Foundation of Zhejiang University of Science and Technology (no. F702104E03 and no. F703104D01), the Program for New Century Excellent Talents in University (no. NCET-11-1080), the first level talent of “151 Talents Project” in Zhejiang Province (no. 12-1-076), the Natural Science Foundation of Zhejiang Province, China (no. LQ13D020001), and the Foundation of China Railway No. 2 Engineering Group Co. Ltd. (no. 201218) are all gratefully acknowledged. The authors wish to thank Professor Zhao Zhiye of the School of Civil and Environmental Engineering, Nanyang Technological University (NTU), for his kind suggestion and the discussions on the relevant topics of this research work. The authors are grateful to Dr. Nie Wen, Dr. He Lei, Dr. An Xinmei, and Mr. Xu Zhipeng from NTU and Dr. Yang Dong from Tongji University for providing useful information and discussion in this paper.