Seismic Response of Tunnel Lining for Shallow-Bias Tunnel with a Small Clear Distance under Wenchuan Earthquake

In order to study the internal force characteristics of shallow-bias tunnel with a small clear distance in earthquake, a large-scale shaking table slope model test was designed, and the geometric scale was 1 :10. In the model test, the Wenchuan (WC) seismic wave was used as the excitation wave.*en, the three-dimensional numerical model was established by usingMIDAS-NX, and the reliability of the numerical model was verified by comparing the acceleration of the test results. *e axial force, bending moment, and shear force of the tunnel cross section and longitudinal direction were calculated by the numerical model under different excitation directions included the horizontal direction (X), the vertical direction (Z), and the horizontal and vertical direction (XZ). *e results show the following. (1)*e internal force of right arch foot of left hole and the left arch foot of right hole is larger than other part of the tunnels because the distance between the two tunnels is smaller and they interact with each other. (2) *e loading direction of single direction loading method is different and the variation trend of tunnel force are different, so the loading direction of seismic wave has a significant influence on the seismic force response of the tunnel. (3) All of the internal force values of tunnel lining under the seismic wave action in bidirection are larger than those in single direction. *e value is not a simple superposition of two directions and has some coupling effect. *e influence of the vertical seismic wave cannot be ignored in dynamic response research. *ese results improve the understanding of the rock slope with small spacing tunnel under seismic action.


Introduction
As an important form of traffic structure, the unsymmetrical loading tunnels usually appear when building roads and railways in hilly, mountainous, sloping areas.Compared with the general tunnel, the bias tunnel has the following characteristics: thin covering rock and soil, asymmetric pressure of surrounding rock, and poor dynamic stability of the whole tunnel.In recent years, earthquakes occur frequently in our country.In Wenchuan earthquake, many tunnels were damaged which including bias tunnels.In all the influence factors of tunnel damage [1], the buried depth is one of the main factors.e tunnel will be constrained greater by surrounding rock with the increase of its depth and will be suffered less damage in earthquake.But in practical engineering, there are many shallow-buried bias tunnels due to topographic conditions, route planning, project cost, and environmental influence.In this condition, the distance between the two tunnels cannot conform to the standard requirements and the pressure on tunnel lining is not even, so many unsymmetrical loading neighborhood tunnels with shallow depth were constructed and, they would be damaged more serious in earthquake.
e dynamic response of tunnels against earthquake has been a subject of intense study after the Wenchuan earthquake and the research, mainly focusing on the experiment and numerical simulation.e shaking table test can simulate the process of earthquake well; it is the most direct and effective method to study the seismic response and failure mechanism of tunnel structure in the laboratory and is one of the important methods in earthquake research.In recent years, many scholars have studied the response of the tunnel under earthquake action by model tests and obtained some valuable results [2][3][4][5][6][7][8][9][10][11][12][13].e longitudinal distances of 100 m in the portal section of Galongla tunnel were selected as the prototype to conduct large-scale shaking table test by Shuping et al. [4]; the results showed that the tunnel vibrates synchronously with the surrounding soil and do not show the natural vibration frequency of itself.In order to research the tunnel structure, seismic responses, and the tunnelsurrounding rock interaction, a large-scale 3D shaking table model test was carried out by Hua et al. [5].Lianjin et al. [6] researched the dynamic response in the portal section of the mountain tunnel with different gradients of upward slope based on large-scale shaking table test.e analysis results showed that an amplified effect of acceleration and displacements of tunnel structures exist at the entrance of the tunnel and there is not much difference in the law of the amplification effect for different gradients of upward slope.In the paper by Feng et al., [7] shaking table tests for seismic response of tunnels were conducted which considering that the tunnel diseases caused by earthquake are different in various depths.In the work of Dong et al., [8] a large-scale shaking table model test is conducted to study the dynamic behavior of entrance slope and its interaction with lining structure of mountain tunnel under earthquake loading.Test results showed that the acceleration response of tunnel entrance slope exhibits obvious amplification effect and surface effect along both vertical and axial directions.
In addition, many scholars also have studied the seismic response of tunnel by numerical simulation.Based on the ABAQUS, Yong [14] studied the effects of tunnel depth, tunnel space, lining stiffness, the properties of surrounding rock, seismic wave incident angle, and seismic intensity on the seismic response of shallow double-hole tunnels.In the paper by Jinsong [15], based on the background of Chongqing Shuangbei tunnel, the three-dimensional numerical model of large section tunnel hole body complex lining was build to research the seismic response calculation and analysis under earthquake.Kunpeng [16] analyzed the seismic dynamic response of unsymmetrically loaded tunnel with the numerical simulation method; it suggested that the unsymmetrical load can increase the entrance length, and the asymmetry of the dynamic response of the tunnel decreased with the increase of the slope soil thickness.e results of numerical calculation are a supplement to the experimental research [16][17][18][19][20][21].
Although numerous investigation and studies on dynamic response of tunnels have been thoroughly accomplished by many researchers, a few studies have been performed about the internal force of shallow-bias tunnel under seismic actions.In this paper, in order to obtain the internal force of shallow-buried tunnels with small distance under seismic action, a shaking table test was carried out.A three-dimensional numerical model was established by using MIDAS-NX, and the reliability of the numerical model was verified by comparing the acceleration of the test results.
is numerical model was used for analyzing the change rule of the internal force of tunnel lining.

Shaking Table Test
2.1.Test Model.In this test, the rigid model box (Figure 1) was used, and its internal dimensions was 3.5 m (length) × 1.5 m (width) × 2.1 m (height).The perspex plates were installed on both sides of the box.e box was opened at one end for filling material, installing sensor, and observing experimental process.
In order to reducing the reflection of seismic waves and the friction on the contact surface of the model and the model box [22,23], the polystyrene foam board was pasted on the inner wall of the left and right sides of the model box, and a layer of PVC film was pasted on the polystyrene foam board.e middle sand and crushed stone were used at the bottom of the model box to increase the friction of the cement mortar and the bottom plate and also reduce the relative slip during the vibration.

Similarity Relation Design.
According to the purpose of this test, the geometric similarity ratio, the acceleration similarity ratio, and the density similarity ratio were used as the main control parameters of the model.e table size, maximum weight capacity, boundary effect of the vibration, and the effect of model material were also comprehensively considered.According to similarity theory, the geometric similarity ratio was 1 : 10, the acceleration similarity ratio was 1 : 1, and the density similarity ratio was 1 : 1; other main similarity constants were deduced by similarity theory and dimensional analysis method.
In this test, it was assumed that the tunnel bias was caused by topographic factors and the gradient of rock slope at the top of tunnel was 1 : 1.5.e rock mass from top to bottom were Weak-weathered rock, Poor rock, and Hard rock, respectively.e rock mass was made of cement mortar which physical and mechanical properties was similar.e mixed proportion from top to bottom is 1 : 6.30 : 1.17 (cement : sand : water) of M7.5 mortar, 1 : 7.25 : 1.35 (cement : sand : water) of M5 mortar, and 1 : 5.58 : 1.04 (cement : sand : water) of M10 mortar, respectively.e tunnel size was designed to be a net height of 60 cm and a net width of 70 cm; the buried depth was 0.9 m and thickness of the middle partition wall was 0.4 m (Figure 2).

Advances in Civil Engineering
e tunnel lining was made of microconcrete which was similar to the physical properties of concrete.Steel bar was simulated with galvanized iron wire.According to the similarity constant, the thickness of tunnel lining was 4 cm.After several times of proportioning test, the optimum tunnel lining material ratio was 1 : 6.9 : 1.3 (cement : sand : water).e strength of lining was converted to 5 MPa in accordance with stress similitude ratio.
e data which included acceleration, dynamic strain, and displacement were collected in the test.Its type was 1221L-002, the range was ±20 m•s −2 , and sensitivity was 2000 mv/g.A total of 20 accelerometer sensors, and 30 strain gauges were arranged in the test.e sensors were arranged on the inner wall of each tunnel lining, separately in the arch foot, arch spandrel, arch vault, because they were the key parts of maximum stress and maximum deformation, as shown in Figure 2. e measuring point number is 1-5 for the right hole and 6-10 for the left hole.

Seismic Wave Loading.
e test was performed in High Speed Railway Construction Technology in Central South University.e vibrating table system consists of one fixed platform and three mobile platforms and these four vibrating tables are in the same line.All of the platforms have six degrees of freedom and the dimensions are 4.0 m (length) × 4.0 m (width); the working frequency of those test systems is 0.1-50 Hz. e maximum displacement of X, Y, Z direction is 250 mm, 250 mm, 160 mm, the maximum acceleration in each direction is ±1.0 g, and the maximum seismic vibration velocity is 1000 mm/s.e Wenchuan artificial wave was input to simulate earthquake.e main frequency was around 5-25 Hz, and each loading time was 20 s, as shown in Figures 3 and 4. Table 1 shows the loading scheme in the test.Seismic wave was applied with peak accelerations of 0.1 g, 0.2 g, 0.4 g, and 0.6 g in X direction and 2/3 of the X direction in Z direction.
e white noise input peak wave was applied before seismic wave loading.e surrounding rocks and linings are modeled using threedimensional solid elements, and the total number of units is 75,979 according to the results, as shown in Figure 5. Physical and mechanical parameters of surrounding rocks and linings are shown in Tables 2 and 3.In numerical calculation, the mechanical behavior of elements obeyed elastic-plastic constitutive relation and Mohr-Coulomb yield criterion.

Verification of Numerical Calculation Model. If
MIDAX/NX numerical model is used to analyze the dynamic response of the tunnel, it must be verified that this model can match the experimental model well.So the seismic wave input is same as the test.After calculation, the peak acceleration at same position of lining is compared.Comparison of acceleration response peak value under WC-Z seismic wave and WC-X seismic wave is shown in Tables 4 and 5, respectively.

Advances in Civil Engineering
e analysis shows that the trend of change for peak acceleration is basically the same and the error is within a reasonable range.e numerical model can be used to calculate and analyze the dynamic response of tunnel under many complex conditions when the experiment cannot be simulated in practical.

Internal Force Analysis of Tunnel Lining Cross Section.
According to the test results and numerical analysis, the acceleration and strain reach the maximum at the peak of seismic wave, so the internal force analysis in this paper refers to the tunnel lining at the peak time of the seismic Advances in Civil Engineering wave.In the following paper, the number of measuring points is L1∼L7 in the left hole and R1∼R7 in the right hole; the detailed measuring points are shown in Figure 6. Figure 7, under the action of horizontal seismic waves, to the left tunnel, the vault and right spandrel are tensiled and the others are compressed.e maximum pressure is on the right arch foot and its value is 499.3 kN.To the right tunnel, the vault is tensiled and the others are compressed; the maximum pressure is on the left arch waist and its value is 488.8 kN.
e change of shear force is similar for both linings, the maximum shear force is on the right arch foot, and the value is −151.2kN for the left hole and 124.3 kN for the right hole.Sign "±" indicates di erent direction of force.
e maximum bending moment is also on the right arch foot; the value is 20,204.6N•m for left hole, and 14,117.8N•m for right hole.Under the action of vertical seismic waves, the axial force of left tunnel lining is pressure, the maximum pressure is on the right arch foot, and its value is 458.8 kN. is is di erent from the horizontal seismic waves.e change of shear force and bending moment is not obvious; the maximum shear force is 137.7 kN and bending moment is 17,600 N•m.To the right hole, the axial force and shear force have similar variation trend as the tunnel lining in the horizontal seismic wave action.But the bending moment of left arch waist changes direction and the inside part is changed to the outside part.e maximum is also on right arch foot and the value is 16,410 N•m.Advances in Civil Engineering rough the force map, the in uence of vertical seismic wave to tunnel lining cannot be ignored; it must be considered in structural design.
e charts also show that all of the internal force values of tunnel lining under the seismic wave action in bidirection (XZ direction) are larger than those in single direction.But the region which contains arch foot of both linings and the left part near the free surface of left tunnel is the largest part of internal force in three seismic wave input conditions.e internal force changes signi cantly in the left side of the left tunnel lining and the arch foot of the right tunnel lining.e reason is that the upper rock layer of the right hole is thick and has good stability, and it is far away from the existing slope and is less a ected by the existing slope.e left hole is close to the existing slope, and the slope has a great in uence on the seismic force response of the right tunnel under the earthquake dynamic action.

Internal Force Changes along the Longitudinal Direction.
Under di erent loading modes, the internal force change along the longitudinal direction is taken as an index to investigate the di erence of seismic force response of tunnel, as shown in Figures 8-10.
As shown in Figure 8, the variation of axial force in the same position is smaller along the longitudinal direction of the tunnel.To the left hole, the axial force reaches the maximum at the arch foot and the value is 502.7 kN at 6 m cross section  Advances in Civil Engineering from the tunnel entrance in the X direction seismic wave.To the right hole, the maximum is left arch foot and the value is 550.6 kN at the 9 m cross section from the entrance in the Z direction seismic wave.rough comparison, it can be found that the axial force at the middle part of tunnel lining is larger but smaller at the entrance and exit of the tunnel.
As shown in Figure 9, it can be found that the shear force is larger at the entrance and exit of the tunnel than the other part.To the left hole, the shear force reaches the maximum at the right arch foot and the value is −149 kN at the cross section 15 m from tunnel entrance in Z direction seismic wave and −148 kN at the cross section of tunnel entrance.To the right hole, the shear force reaches the maximum at the right arch foot and the value is −144 kN at the cross section 15 m from tunnel entrance in the XZ direction seismic wave.From Figure 9, it can be seen that the response of the left hole is more intense than the right hole.e main reason is that in the vertical direction of loading wave, the overlying strata of the right hole is thicker, the action force of the right hole is greater, and the "rise e ect" has a signi cant impact on the shear force response of the right hole.By comparing Figures 8 and 9, it is suggested that the loading direction of single   Advances in Civil Engineering direction loading method is different, the variation trend of tunnel force are different.It shows that the loading direction of seismic wave has a significant influence on the seismic force response of the tunnel.As shown in Figure 10, the bending moment increases significantly at the arch foot of the tunnel in all loading modes.e maximum is located at the right arch foot of the left hole and the value is 20,339 N•m at the cross section of tunnel exit in x direction seismic wave.
By comparing Figures 8-10, it is summarized as follows: Because distance between the two tunnels is closer, the internal force of right arch foot of left hole, and the left arch foot of right hole is larger than other part of the tunnels because they interact with each other.
Because of the influence of the existing slope on the internal force response of the tunnel, the variation trend of the internal force of the left and the right hole never show a symmetrical trend.
It can be seen that under the action of bidirectional wave loading, the variation trend of internal force in is more complex, but this trend is similar to X direction wave loading.e reason may be that the response of the vertical direction is stronger than that of the horizontal direction under the action of the bidirectional wave.

Conclusions
e large-scale shaking table model test was carried out, and the acceleration response of the unsymmetrical loading neighborhood tunnel with shallow depth model in different directions and different acceleration peaks of Whenchuan seismic wave was obtained.e three-dimensional numerical model of the unsymmetrical loading neighborhood tunnel with shallow depth was established by using MIDAS-NX.e acceleration and internal forces of the tunnels under the same working conditions were obtained, and the reliability of the numerical model was verified by comparing the acceleration of the tunnel lining.
e tunnel lining dynamic response under different loading modes (X direction, Z direction, and XZ direction) are studied.
e axial force, shear force, and bending moment of the tunnel cross section were calculated by the numerical model under the earthquake simulation.
e following conclusions are obtained: (1) e internal force of the right arch foot of the left hole and the left arch foot of the right hole is larger than the other part of the tunnels because the distance between the two tunnels is closer and they interact with each other.(2) All of the internal force values of tunnel lining under the seismic wave action in bidirection are larger than those in single direction.e value is not a simple superposition of two directions and has some coupling effect.e influence of the vertical seismic wave cannot be ignored in dynamic response research.

3. 1 .
Numerical Simulation Model.In this paper, MIDAX/NX software is used to analyze the dynamic response of the tunnel.e main parameters of the numerical model are as follows: the total length is 40 m; the height is 17 m; the length of tunnel is 15 m; and the thickness of tunnel lining is 0.4 m.

Figure 8 :
Figure 8: e axial force of tunnel lining along longitudinal direction.

Figure 9 :
Figure 9: e shear force of tunnel lining along longitudinal direction.

Figure 10 :
Figure 10: e bending moment of tunnel lining along longitudinal direction.

Table 1 :
Loading rule of shaking table test.

Table 4 :
Comparison of acceleration response peak value in the vertical direction (unit: m•s −2 ).

Table 5 :
Comparison of acceleration response peak value in the horizontal direction (unit: m•s −2 ).