Abstract

To investigate the seismic response of large-scale tunnel in compacted clay and effect of shock absorbing layer to the tunnel, a series of three dimensional (3D) shaking table model tests were carried out. The similarity ratio of the model is 1 : 8 and the size of the model container is 9.3 m (length) × 3.7 m (width) × 2.5 m (height). The cross-sectional diameter of the model tunnel is 0.9 m, and the thickness of the tunnel lining is 0.06 m. To simulate the clay soil surrounding condition, the container was filled with clay soil. During the tests, the concrete strain, acceleration, and dynamic soil pressure on the surface of the model tunnel were measured. The results show the existence of tunnel can decrease the maximum acceleration of the model in the X direction; the shock absorbing layer can further decrease the maximum acceleration, however, cannot change the dominant frequency of the ground motion. The longitudinal and hoop strain of the model tunnel with excitation of the input motion is mainly in tension state and the maximum hoop deformation of the model tunnel is located at the conjugate 45°. In addition, the shock absorbing layer has an effect on the strain and dynamic earth pressure of the model tunnel.

1. Introduction

Tunnel behaved better than above ground structure during the earthquake loading. However, earthquake events in recent years show that the underground structure has serious damage in seismically active area; for example, the collapse of the Daikai subway metro station in Kobe, in Japan (1995) [1]; failure of Bolu tunnel in Turkey (1999) [2]; damage of LongXi tunnel in WenChuan (2008) [3]. These failures revealed that the current practices of seismic design for underground structures cannot ensure the seismic performance of tunneling structures. Therefore, how to improve safety of tunnels during earthquakes has become one of the major concerns for the engineering.

Many scholars have investigated the seismic response of tunnel during the earthquake loading through theoretical analysis [411], which has a certain significance of engineering guiding in seismic of tunnel. However, the results of these methods may substantially deviate, even under the same assumptions. Therefore, it is a better way to find the weakness parts of tunnels and explore the structure of the failure mechanism of tunnels through the model test. Cheng et al. [12] found that soil liquefaction significantly influenced the propagation of seismic motions and the dynamic responses of the tunnel. Ulas et al. [13] carried out dynamic centrifuge modeling to circular tunnel in dry clay (scale 1 : 50); largest axial dynamic lining forces were measured near shoulders. Tsinidis et al. [14] found the horizontal acceleration recorded at several locations was generally amplified towards the soil surface of square tunnel in dry clay, while the tunnel presence affected this amplification (scale 1 : 50). Tsinidis et al. [15] revealed a minor effect of the interface properties on the tunnels racking deformation (scale 1 : 40). Jiang et al. [16] found that the soil pressures were symmetrically distributed on the left and right side walls. Cilingir and Madabhushi [17] conducted centrifuge modeling of square tunnel to investigate the effect of different shallow depths of clay deposits. The results show that the depth of the tunnel does not affect the deformation pattern of the tunnel significantly during an earthquake event (scale 1 : 45). The above literature of model tests has a common feature, that is, the small scale model testing, which has a certain value for investigating the seismic response of tunnel. However, these small-scale model tests can oversimplify the complexity of the detail and response characteristic of a tunnel structure. It is even more difficult to show the realistic response to the lining force and failure characteristics of the tunnel. Moreover, Andrew [18] found that as the difference in the scales of the model and prototype increases, the difference in the three-dimensional and time-varying properties also increases. Therefore, based on reducing the difference brought by the small scale and the influence of boundary effects, in order to investigate the dynamic response of a tunnel and the soil-tunnel interaction, it is necessary to carry out the shake table test of large-scale (1 : 8) shield tunneling models in a 9.3 m-long, 3.7 m-wide, and 2.5 m-high rigid box.

The objectives of this study are (a) to provide experimental data on the seismic response of circular tunnels in compacted clay with ground motions excitation; (b) to determine the differences in the measured accelerations, strains, and dynamic earth pressures between the cross sections with and without a shock absorbing layer (SAL); (c) to evaluate the efficiency of the SAL during earthquake loadings. In this paper, a series of shake table tests on large-scaled circular model tunnel embedded in clay soil are initially presented. The tests were conducted at Shake table laboratory of National Engineering Laboratory for Construction Technology of High Speed Rail in China. Following a detailed description of the experimental set up, a detailed description of the experimental setup is provided, focusing on the design of the soil container, clay soil tested and the scaled model tunnel, instrumentation, and test procedure. The main results obtained from the tests are summarized and discussed experiment measurements in terms of concrete strain, dynamic earth pressure, and acceleration. The collected experimental data are primarily to be used to better understand the seismic behavior of circular tunnel and further optimizing seismic design methods for tunnels.

2. Test Setup and Model Preparation

2.1. Shake Table System

The shake table tests were performed using the Laboratory Multifunction Shaking Table System (LMSTS) in the National Engineering Laboratory of High Speed Railway Construction Technology in Central South University, China. The LMSTS has an electrohydraulic servo drive, and it consists of three 4 m × 4 m 6-DOF mobile stations, in which can be controlled independently to simulate single or multiple three-directional ground motions. The LMSTS was built in 2009 and equipped with the advanced instruments high-performance dynamic signal acquisition system, which carry out independent sampling of 320 input channels for accelerometers, strain gauges, and earth pressure sensors. Figure 1 shows the layout of LMSTS and the horizontal and vertical actuators. For the tests conducted in this study, only Tables A and B were used. Main parameters of shake table facility are summarized in Table 1.

2.2. The Soil Container
2.2.1. Design of Soil Container

The soil container is a rigid box with the dimensions of 9.3 m (length) × 3.7 m (width) × 2.5 m (height), as shown in Figure 2. The soil container includes the box wall, the inner framework, and the outer framework. The box wall was made by steel plates and the thickness of the plates is 6 mm. The inner framework consists of steel strips welded onto the box wall to ensure the box is rigid; the outer framework consists of steel angles welded between steel plates. Two transverse ribs were welded onto the top of the model container. The bottom of the box is a 5 mm thick steel plate welded on to the frame. The inner bottom of the box was divided into 1.0 m2 baffle plates using welded steel angles. The total weight of the empty box is 3,800 kg, approximately. The soil container was bolted firmly onto tables A and B using a total of 64 high-strength bolts. The length/height ratio of the container is 3.72, and this design can reduce the boundary effect of the rigid end wall [19, 20].

2.2.2. Boundary Condition

The boundary condition is crucial for the result of the shake table tests and a wave-absorbing boundary condition is needed. The flexible materials are generally installed inside the rigid container to reduce the boundary effect of the box wall during shake table tests and these flexible materials are usually polystyrene foam boards and sponge rubber [2124]. In this study, a 10 cm-thick polystyrene foam board was attached to the container walls; meanwhile, a 5 cm-thick and a 2.5 cm-thick sponge rubber boards were attached to the end walls and the side walls, respectively.

2.3. Model Tunnel and Materials

The model tunnel is scaled at 1 : 8. The prototype tunnel has a cross-sectional diameter of 7.2 m and a lining thickness of 0.48 m. According to the similarity law [25], the cross-sectional diameter of the model tunnel is 0.9 m and the thickness of the tunnel lining is 0.06 m. Table 2 shows the similarity ratios of the model tunnel.

The steps of how to install the model tunnel and tamper the clay soil are as follow:

Firstly, at the base of the soil container, a layer of 2-3 cm-diameter aggregate was placed to increase the frictional resistance between the model clay and the floor to reduce potential slippage.

Secondly, the clay soil was placed into soil container and the thickness of rammed layer is 20 cm. When the clay soil was buried to a height of 30 cm, nine model tunnels were installed in the middle of the soil container along the X direction. The clay soil was compacted to the height of 191 cm, which mainly simulated the dynamic response of the shield tunnel in clay stratum. The clay soil was brought from YueLu Mountain, and its properties were tested using samples taken directly from the soil container [26, 27]. Table 3 shows the physical properties of the compacted clay.

As shown in Figure 3, on the basis of equivalent stiffness principle, the tunnel lining is made of microconcrete with an inside diameter of 0.78 m and a thickness of 0.06 m. It is reinforced with steel mesh (6 mm diameter) at the radial ring spacing of 10 mm. Table 4 shows the typical mechanical properties of the model tunnel used in the tests. A 5 cm-thick foam board was used to wrap a section of the model tunnel (Figure 4) to investigate the effect of the SAL on the tunnel response under seismic shaking, and previous dynamic finite difference method revealed that 5 cm-thick foam board can significantly reduce the seismic dynamic response of the second lining of tunnels. The shake table tests also revealed that the strains of lining are reduced by the SAL [2830]. Table 5 shows the mechanical properties of the foam board.

2.4. Measuring Devices and Instrumentation Layout

Three sets of sensors were installed on the model tunnel to monitor the model tunnel seismic response during the tests; the measuring rang and minimum precision of the capacitive micromachined accelerometers (model 2210 by Vigor,T., Ltd.) are ±2 g and 0.001 g, respectively, and the measuring rang of microresistance strain earth pressure cells (Zhuochuan Electronic Technology of Beijing Co., Ltd) is ±300 kPa and the minimum precision is 0.001 kPa; the resistance strain gauges (model BX120-100AA, Taizhou Technology Co., Ltd) has a measuring rang of about 0∼12000 με and the minimum presion is 0.001 με. Figure 5 shows the layout of the instrumentation.

As shown in Figure 5(a), two main cross-sections were designed for monitoring in the soil container. One is the cross section # 1 which has one model tunnel, the other is the cross section # 2 which has one model tunnel wrapped by a 5 cm-thick foam board. As shown in Figure 5(b)-5(c), three sets of sensors have been installed in the rigid box. The first set of accelerometers is used to measure the accelerations in the X, Y, and Z directions of tunnel lining. The second set is used to measure the accelerations of the model clay, and the third set is used to measure the accelerations on the shake table. The dynamic earth pressure cells are installed to measure the earth pressure at the interface of the clay, and the strain gauges are installed to measure the deformation of the inner and outer surfaces of the model tunnel.

2.5. Experimental Program and Test Procedure

During the tests, white noise and harmonic motions frequency sweeps and the recorded earthquake input motions were used. White noise sweep can determine the approximate range of the natural frequency of the model casing. Harmonic motions sweep can accurately confirm the natural frequency of the soil container. In order to prevent any adverse effect of damage accumulation on the test results, the sweep motion tests of white noise and Harmonic motions have been conducted with 0.05 g peak acceleration. The results of the tests are as follows. The natural frequency of the empty soil container in the X direction is 35.31 Hz. However, the natural frequency of the clay soil in the model in the X direction is 4.62 Hz. This is an indication that the shear stiffness of the soil container is much higher than that of the clay in the model, which can prevent the resonance effect between the soil container and the clay with the ground motion excitation.

Two real earthquake records from the Kobe earthquake measured in 1995 at OKA Station and the Chichi earthquake measured in 1999 at CHK station were used as input motions in the shake table tests, which came from the similar geological conditions of PEER Strong motion Database. Duration and frequency of the earthquake records were scaled for model input motions by a scale factor of 8, as shown in Figures 6 and 7. The input motions were scaled up to different peak accelerations during the testing program. Four cases of peak acceleration were considered and scaled to 0.05 g, 0.1 g, 0.2 and 0.4 g, as tabulated in Table 6.

3. Test Results and Discussion

3.1. Profiles of Maximum Acceleration with Depth

Figure 8 shows the profiles of the maximum horizontal acceleration in the X direction with model tunnel and the surrounding soil for events EQ1, EQ3, EQ5, and EQ7. The profiles of these maximum accelerations from the base of soil container to the surface of the clay (i.e., Sensors A1, A5, A7, and A4). Figure 9 shows the profiles of the maximum horizontal acceleration in the X direction with model tunnel and the surrounding soil for events EQ2, EQ4, EQ6, and EQ8. It can be seen that the maximum acceleration increases gradually from the shake table to the surface of the model clay in “free field” vertical array. However, the maximum acceleration of model tunnel firstly decreases from the shake table to the base of soil container; the reason is that the weight of the model tunnel is more than 300 kg and it influences the clay motion below the bottom of model tunnel, and then, these accelerations gradually increase from the bottom of model tunnel to surface of the clay in the model tunnel array. Then, a remarkable increase in acceleration takes place from the bottom to the crown of model tunnel compared to the same depth “free field,” which mainly attributes to the damping of wave propagation in the macroconcrete is lower than the clay soil.

Due to space limitations, only the results for excitation of EQ1, EQ3, EQ5, and EQ7 are summarized in Figure 10, where they are plotted in terms of the ratio between the peak accelerations measured by from the base to the surface of clay for vertical array of the tunnel and the surrounding soil. This ratio is the conventional amplification factor and can be used to investigate the influence of the tunnel on the vertical propagation of shear motions. The amplification factor at the base of clay (A1) in tunnel vertical array is lower than unity for all earthquakes, compared with same position (A15) in free field vertical array, which was attributed to the influence of large-scale model tunnel. The amplification factor at the bottom of model tunnel also is lower than unity, while the amplification factor at the crown of model tunnel (A7) is larger than the same position (A17) in the surrounding soil. This could be an explanation for the nonuniformity behavior of the clay, which is due to the arching effects, as well as the influence of the cavity boundary, reflecting down some of the shear waves propagating upward.

3.2. Effect of SAL on Tunnel Acceleration

Figures 11 and 12 show a comparison between acceleration time history and corresponding spectrum of with SAL and without tests for the bottom (A5) and crown (A7) of the model tunnel subjected to the Kobe with peak acceleration levels (0.05 g, 0.1 g, 0.2 g, and 0.4 g), respectively. With the different peak acceleration levels excitation, the maximum accelerations and the maximum amplitude of frequency with the SAL are still lower than those of the model tunnel without SAL, which means that the SAL can reduce the acceleration of the model tunnel. On the other hand, the main frequency of the model tunnel without SAL is the same as the model tunnel with SAL and the difference is the size of the amplitude of frequency, which means that the SAL cannot change the dominant frequency of ground motion.

Figure 13 shows a comparison of the maximum acceleration development for the tunnel with and without the SAL. These accelerations were recorded using the X-direction accelerometers. It can be seen that the maximum acceleration firstly decreases from the shake table to the bottom of model tunnel; then, these accelerations gradually increase from the bottom of model tunnel to ground surface of the clay in the model tunnel array of section # 2. It is interesting that the point of minimum acceleration is located at the bottom of model tunnel with SAL. Meanwhile, the maximum accelerations at the crown and at the bottom (i.e., location A7 and A5) of the model tunnel with the SAL are generally lower than those of the model tunnel without, as the peak acceleration of input motion increases, this phenomenon becomes more pronounced.

Figure 14 shows the amplification factor of the bottom and crown of the model tunnel with and without SAL. It can be seen that at the bottom (at location A5) and at the crown (at location A7) of the model tunnel, the amplification factor of the model tunnel with the SAL is no more than 0.7; however, the amplification factor of the model tunnel is no less than 0.8, which means the SAL further reduces the amplification factor of the model tunnel. It can also be seen that there is a tendency for the peak acceleration of the tunnel to attenuate from the bottom to the crown of model tunnel with the SAL. This shows that the foam board has positive effect on the reduction of the acceleration.

Figure 15 shows the reduction factor of the SAL with different peak acceleration of input motions. The calculation of the reduction factor is shown in Equation (1). It can be seen that with small peak accelerations, the reduction factor of the SAL at the bottom of model tunnel is the same as the crown of model tunnel. The reduction factor of the SAL to the bottom and crown of the model tunnel is more than 50% with the 0.1 g excitation. However, as the peak accelerations of input motion increase (i.e., 0.2 g and 0.4 g), the reduction factor of the SAL gradually decreases.where, is the maximum acceleration of model tunnel with SAL and without; is the maximum acceleration of model tunnel without SAL; is the maximum acceleration of model tunnel with SAL; and is the maximum acceleration of shaking table.

3.3. Tunnel Deformation

Because of the variation of soil deformation over the depth, the deformation of the cross section of the tunnel also changes (as shown in Figure 16).

Figure 17 shows the maximum longitudinal and hoop strains at the inner and outer surfaces of the tunnel under excitation in the X direction (i.e., for events EQ1, EQ3, EQ5, and EQ7).

Figure 18 shows the maximum longitudinal and hoop strains at the inner and outer surfaces of the tunnel under excitation in the X direction (i.e., for events EQ2, EQ4, EQ6, and EQ8). The deformation characteristics of the model tunnel are described using the absolute maximum strain in the longitudinal and hoop direction. A positive strain value indicates a “tension” state of model tunnel, and a negative strain value indicates a “compression” state of model tunnel.

Figures 17(a)17(c) show the longitudinal strain at the outer and inner surfaces the model tunnel with the Kobe motions excitation, respectively. It can be seen that at the small peak acceleration (0.05 g), the longitudinal strain of the model tunnel is in tension and no more than 5 με. As the peak acceleration of input motion increases, the tension strain gradually increases, the maximum tension strain of inner and outer surface of the model tunnel are located at the 0° and 270° of the model tunnel, respectively. Figures 18(a)18(c) also show the longitudinal strain at the outer and inner surfaces the model tunnel with the Chichi waves excitation, respectively. It reveals that at the small peak acceleration (0.05 g), the longitudinal strain of the model tunnel is in tension and no more than 10 με. As the peak acceleration of input motion increases, the tension strain gradually increases, the maximum tension strain of inner and outer surfaces of the model tunnel are located at the 315° and 270° of the model tunnel, respectively. During the tests, the longitudinal strain of the model tunnel along the X direction is in tension with the different seismic motions, and the maximum longitudinal strain is no more than 50 με. Figures 17(b)17(d) show the hoop strain of the outer and inner surface of the model tunnel with different peak acceleration Kobe motions. Under the EQ1 excitation, the hoop strain of the model tunnel is in tension and no more than 10 με. With the increasing of input motions, the maximum hoop strain of the outer and inner surface is located at 225° of the model tunnel. The value of maximum hoop strain of the inner and outer surface is 39.12 με and 49.56 με, respectively. Figures 18(b)18(d) also show the hoop strain of the outer and inner surfaces of the model tunnel with different peak acceleration Chichi motions. Under the EQ2 excitation, the hoop strain of the model tunnel is in tension and no more than 10 με. With the increasing of input motions, the maximum hoop strain of the inner and outer surface are located at 45° and 315° of the model tunnel, respectively. The value of maximum hoop strain of the inner and outer surface is 57.11 με and 37.93 με, respectively. As can be seen, the hoop strain of the model tunnel is in tension, and then, the maximum deformation of the model tunnel is located at the 45° conjugate angle.

The deformation of the tunnel can be strongly influenced by the ground deformation. Under the ground motion, the ground deformation developed from underground to the ground surface, as shown in Figure 16. The development of shear strain in the clay due to base excitation caused a shear force at the tunnel.

3.4. Effect of SAL to Concrete Strains

Figure 19 shows the maximum longitudinal and hoop strains at the inner and outer surfaces of the model tunnel with the SAL under excitation for events EQ1, EQ3 EQ5, and EQ7.

Figures 19(a)19(c) show the longitudinal strain at the outer and inner surface the model tunnel, respectively. It shows that as the peak acceleration of input motion increases, the tension strain gradually increases. However, the maximum tension strain of inner and outer surface of the model tunnel is smaller compared to that of model tunnel without SAL, and the same phenomenon was observed in the Figures 19(b)19(d).

The maximum strain of the model tunnel with SAL and without is summarized in Table 7. In view of difference of the strain of the model tunnel, the reduction factor of the SAL to the strain of model tunnel was shown by the ratios of the difference between the strain of model tunnel and the strain of the model tunnel with SAL. The longitudinal strain of the model tunnel for the SAL can reduce more than 20% when compared with the model tunnel without SAL; however, the hoop strain of the model tunnel with the SAL can reduce more than 20%.

3.5. Soil-Structure Interaction

In this study, the soil-structure interaction is investigated through the dynamic earth pressure between the model tunnel and the clay. The dynamic earth pressure cells have measured the dynamic earth pressures at four symmetrical points of the model tunnel (Figure 5(b)). Figure 20 shows the maximum dynamic earth pressures for different peak acceleration levels (0.05 g, 0.1 g, 0.2 g, and 0.4 g). It can be seen that the maximum earth pressure took place at the bottom of the model tunnel and the reason is the weight of model tunnel is about 300 kg, and the earth pressure at the bottom of the model tunnel is caused by the inertia of model tunnel weight and clay-model tunnel interaction with excitation of ground motion. The earth pressures at the crown of model tunnel are higher than those on the side of model tunnel.

3.6. Effect of the SAL to Soil-Structure Interaction

Figure 21 shows a comparison of the maximum dynamic earth pressures of the model tunnel with and without the SAL under the excitation of E1, EQ3, EQ5, and EQ7. The earth pressures of the tunnel with the SAL were lower than those of the tunnel without SAL. Figure 22 shows the reduction factor of the SAL to the dynamic earth pressure subjected to different peak acceleration levels. The calculation of the reduction factor is shown in Equation (2). It can be seen that the reduction factor of the SAL firstly decreases from 0.05 g to 0.1 g of input motions; then, with the increase of the input motions, the reduction factor of the SAL also increases. The SAL can reduce 10%–50% dynamic earth pressure the of the tunnel model without SAL. The reduction factor of the SAL to earth pressure is minimum when the peak acceleration is 0.1 g. is the earth pressure of model tunnel without SAL; is the earth pressure of model tunnel with SAL.

4. Conclusion

This paper presented experimental results obtained from a series of shake table tests on a circle model tunnel embedded in compacted clay. A soil container filled with clay was used. Using strain gauges, accelerometers, and dynamic earth pressure cells, the seismic response of the tunnel during shake table experiment was measured. The dynamic response of the model tunnel with and without the SAL material has also been compared and analyzed. From the results, we can summarize the following:(i)The peak accelerations of the X direction increase gradually from the shake table to the surface of the model clay in “free field” vertical array; the peak accelerations in the X direction firstly decrease and then increase from the base of soil container to the surface of the model clay. In the model tunnel vertical array, due to the weight of model tunnel and the difference of wave propagation in the macroconcrete and clay.(ii)The SAL can reduce the acceleration of model tunnel, and the horizontal acceleration of the model tunnel with SAL in X direction was just 30–60% of that for tunnel lining without a SAL. The reduction factor of SAL to the acceleration firstly increases; however, with the increasing of peak accelerations, the reduction factor of the SAL decreases instead.(iii)The deformation characteristic of the model tunnel is mainly in tension, and as the peak accelerations of input motion increase, the tension strain of the model tunnel in longitudinal and hoop direction also increases; however, the maximum tension strain of model tunnel is located at the 45° conjugate angle due to the shear force of input motion.(iv)The longitudinal strain of the model tunnel for the SAL can reduce more than 20% that of the model tunnel without with the main X excitation; however, the hoop strain of the model tunnel with the SAL also can reduce more than 20% that of the model tunnel without, which means the SAL can effectively reduce the strain of the model tunnel in longitudinal and hoop direction.(v)The dynamic earth pressures were larger at the bottom of the model tunnel compared to the other position of model tunnel(vi)The SAL can reduce 10%–50% dynamic earth pressure of the tunnel model without SAL. The reduction factor of the SAL firstly decreases, and then, as the peak accelerations of input motion increase, the reduction factor of the SAL also increases.

Data Availability

All data generated or analyzed during the study are included in this paper.

Conflicts of Interest

The authors declare that they have no conflicts of interest regarding the publication of this paper.

Acknowledgments

This research was supported by the China Railway Corporation Key Project (grant no. 2014G010-G) and National Engineering Laboratory for Construction Technology of High Speed Railway. All supports are gratefully acknowledged.