Abstract
With the development of urban tunnels, more and more studies focus on the antiseismic design of tunnels. Setting up a shock absorption layer is a very effective means in the earthquake resistance of an urban tunnel. In this paper, in order to study the influence of the thickness of a shock absorption layer on the antiseismic effect of the urban shallow buried double arch rectangular tunnel, the Tengzhou–Fenghuang tunnel is used as the research background. Firstly, the model is established by finite element software to analyze the tunnel with 50 mm, 100 mm, and 150 mm thick shock absorption layers, and analyze the tunnel displacement, stress, and safety factors of different models under the action of a seismic wave. Finally, the calculation results are compared and analyzed with those of the model without the shock absorption layer to summarize the thickness of the shock absorption layer suitable for urban shallow buried rectangular tunnels. The results show that the displacement and stress of tunnel lining are significantly reduced, and the safety factor is significantly improved after setting the shock absorption layer. Among them, the tunnel with a 100 mm thick shock absorption layer has the greatest reduction in principal stresses, the greatest increase in safety factor, and the best antiseismic effect, while the construction costs are not high. It is recommended that the urban shallow buried double arch rectangular tunnel be constructed with a 100 mm thick shock absorption layer.
1. Introduction
In recent years, many tunnels have been damaged in the earthquake [1, 2]; therefore, there are more and more research studies on the antiseismic technology of tunnels. At present, there are three aspects in the research of improving the tunnel antiseismic performance, which are to improve the strength of tunnel surrounding rock, enhance the toughness and strength of tunnel lining, and construct a shock absorption layer. The main way to improve the strength of surrounding rock is to set anchor bolts to grout the surrounding rock of the tunnel so as to improve the strength of surrounding rock. At present, the research on surrounding rock grouting has been relatively perfect, and the main research achievements include exploring the working mechanism of grouting anchor [3], the interaction between anchor bolts and surrounding rock [4], and evaluating the grouting effect of the surrounding rock [5–7]. Adding steel fiber, glass fiber, and other materials in the concrete can greatly improve the mechanical properties of concrete [8, 9], and using this composite concrete as a lining material in tunnel construction can effectively improve the antiseismic performance of tunnel lining. What is more, there are many research studies on fiber reinforced concrete. Through shaking table test, it is found that polypropylene fiber concrete as lining can reduce the stress and strain of a tunnel under earthquake action [10]. Steel fiber reinforced concrete has better antiseismic performance than ordinary reinforced concrete [11, 12]. Although there are many research studies on fiber reinforced concrete, the application of fiber reinforced concrete in engineering is less, and its antiseismic effect in specific projects is not clear, so further research is needed.
The shock absorption layer is added between the tunnel lining and the surrounding rock, which can reduce the transmission of adverse effects of surrounding rock to the tunnel lining. Xin et al. [13] thought that setting a shock absorption layer is the most simple and effective means to improve the antiseismic performance of a tunnel. The shock absorption layer can reduce the effect of a seismic wave on lining and reduce the pressure from surrounding rock. Huang et al. [14] carried out an antiseismic test on a circular tunnel with different material shock absorption layers through shaking table test, and the results show that the sponge rubber plate has a good antiseismic effect. In addition, the exploration of new damping materials is also carried out. Jiang and Zhang [15] studied the seismic behavior of foam concrete. The results show that foam concrete can prevent lining from being damaged under an earthquake. Su et al. [16] were the first to apply the closed-cell aluminum foam in designing a shock-absorbing layer for a high-speed railway tunnel and studied the seismic behavior of closed-cell aluminum foam. Shaking table test and numerical simulation are two particularly effective methods to study the antiseismic performance of the shock absorbing layer. Many scholars [17–19] have discovered the antiseismic mechanism and antiseismic performance of the shock-absorbing layer through shaking table test and numerical simulation.
Compared with the horseshoe tunnel and the circular tunnel, the rectangular tunnel has uneven stress distribution. Therefore, the rectangular tunnel has more stringent requirements on the antiseismic technology of the tunnel, and there is a lack of research on the antiseismic technology of the rectangular tunnel at present. In order to explore the influence of the thickness of the shock-absorbing layer on the antiseismic performance of an urban shallow buried double arch rectangular tunnel, based on the Fenghuang tunnel, this paper analyzes the displacement, stress, and safety factor of the tunnel lining with three different thickness shock absorption layers of 50 mm, 100 mm, and 150 mm, and combined with the engineering practice, obtains the appropriate antiseismic measures for the urban shallow buried double arch rectangular tunnel.
2. Tunnel Situation
Fenghuang tunnel is located in a flat terrain with plain landform. The upper part of the stratum is clay and the lower part is mainly sandstone (rock quality grade is IV) [20]. The seismic fortification intensity in this area is 6 degrees, and the basic design seismic acceleration value is 0.05 g [21].
The tunnel adopts a closed box structure. There is a layer of 0.2 m cushion between the bottom of the tunnel and the surrounding rock. The cross section of the tunnel is shown in Figure 1. C40 concrete [22] is used for tunnel lining, and C20 concrete is used for cushion. The span of the tunnel is 26.3 m, and the height is 8.2 m. The thickness of the roof, bottom plate, and the middle wall of the tunnel lining are 0.8 m, 1.1 m, and 0.8 m, respectively.

3. Software General Situation and Calculation Condition
3.1. Simulation Method
ABAQUS is a general finite element software, which mainly consists of two analysis modules, ABAQUS/Standard and ABAQUS/Explicit. It can be used to solve linear and nonlinear problems. It is widely used in many scientific research fields and engineering fields such as machinery, water conservancy, materials, and civil engineering. Compared with other software, ABAQUS has several advantages. First, it has powerful functions. This software is a simulation software integrating mechanics, electromagnetism, and thermodynamics, and its operation is very simple. Second, it can solve many kinds of nonlinear problems; in nonlinear analysis, ABAQUS can automatically select corresponding load increments and convergence rates so as to ensure accurate solutions in the analysis process. Third, it has a rich material model library and cell library, which can select different models and elements according to different problems.
At present, many experts have carried out an antiseismic calculation of slope, house, and underground structures by ABAQUS. The antiseismic calculation methods of underground structures mainly include the static method, quasistatic method, and dynamic response analysis method. The error of the static method is large, and now, the main method is the quasistatic method. In the antiseismic calculation of the underground structure, due to the size limitation of the model, it is necessary to transform the semi-infinite boundary into the model with a boundary. The seismic wave will reflect at the boundary of the model, which will interfere with the calculation results. Therefore, appropriate boundary conditions should be set up when the model is established. At present, the commonly used boundary conditions include viscous boundary, viscoelastic boundary, and infinite element boundary. The infinite element boundary can simulate the actual situation of zero displacement of soil at infinity and reduce the calculation amount of the model. Therefore, the infinite element boundary is selected as the boundary condition in this paper.
3.2. Calculation Model
The model is established according to the geological conditions and the tunnel size. When the model is established, the buried depth of the tunnel is 8 m, and the width of the tunnel is 26.3 m. The width of the model is about 10 times of the tunnel, so the width of the model is about 290 m, and the thickness of the model is 20 m. The upper part of the model is silty clay with a thickness of 8 m, the middle part is grade IV surrounding rock with a thickness of 38 m, and the bottom is grade II surrounding rock with a thickness of 20 m. When the model is established, the Mohr Coulomb constitutive model is adopted for surrounding rock, and the linear elastic constitutive model is adopted for the tunnel lining, cushion, and shock absorption layer. In order to make the calculation results accurate and reduce the amount of calculation, the mesh division of tunnel lining is more compact and that of surrounding rock is loose. The boundary condition of the model is that there is no rotation and displacement at the bottom and around, and the infinite element boundary is set at the same time, and the upper boundary is unconstrained. The calculation model is shown in Figure 2.

(a)

(b)
3.3. Calculation Parameters
The surrounding rock conditions are selected according to geological survey data. The material of the shock absorption layer is rubber plate. The calculation parameters are shown in Table 1.
3.4. Calculation Conditions
This paper studies the influence of the thickness of the shock absorption layer on the antiseismic effect of the urban shallow buried double arch rectangular tunnel. There are four groups of calculation conditions, one of which has no shock absorption layer as the control group, and the other three conditions have the shock absorption layer of different thicknesses. The calculation conditions are shown in Table 2.
3.5. Measurement Arrangement
One monitoring section is set every 5 m in the tunnel, and a total of 5 monitoring sections are set, which are S1∼S5, respectively; the section distribution is shown in Figure 3. Nine monitoring points are set up for each monitoring section, and the principal stress, shear stress, and safety factor of each measuring point are analyzed. The arrangement of measuring points is shown in Figure 4.


3.6. Load Input
The seismic fortification intensity of the Tengzhou area is 6 degrees, so the acceleration wave of the Wolong station of the Wenchuan earthquake is selected as the seismic wave, which is standardized according to the intensity of 7 degrees. The peak acceleration of the seismic wave is 0.1 g, and the duration is 15 s. The loading mode of this numerical simulation is conventional dynamic loading, and the seismic waves in three directions (x, y, z) are transmitted to the upper part at the same time through the bottom of the model. The acceleration curve is shown in Figure 5.

(a)

(b)

(c)
3.7. Calculation Method of Antiseismic Effect
In order to check the antiseismic performance of the tunnel, it is usually necessary to calculate the safety factor of the lining. The safety factor is calculated according to the specification [23] as follows:
In the above formula, N represents the axial force; E represents the modulus of elasticity; ε1 and ε2 are the internal and external strains of tunnel lining; h represents section thickness; b is the section width, taken as 1 m.
In the above formula, M represents the bending moment, and others are the same as formula (1).
In formulas (3) and (4), K represents the safety factor of tunnel lining; φ represents the coefficient of longitudinal bending; α represents the influence coefficient of axial eccentricity; Ra represents the ultimate compressive strength; e0 represents the eccentricity; Rl represents the ultimate tensile strength, and the rest is the same as formula (1).
According to the simulation results, the antiseismic effects of displacement, principal stress, and safety factor of different thickness shock absorption layers are calculated. The calculation formula of antiseismic effect is as follows:
In the above formula, ρ represents the antiseismic effect of each simulation result; Qa represents the value of maximum convergence of side wall, principal stress, and safety factor without the shock absorption layer (condition 1); Qb represents the value of maximum convergence of side wall, principal stress, and safety factor when setting different thickness shock absorption layer (condition 2 ∼ condition 4).
4. Numerical Simulation Results and Antiseismic Effect Analysis
4.1. Maximum Convergence Value of Side Wall
The maximum convergence value of the side wall of the tunnel lining under each condition is extracted, and the decreased amplitude of the maximum convergence value of the side wall after setting the shock absorption layer (condition 2 ∼ condition 4) is calculated (compared with condition 1). The calculation results are shown in Table 3.
It can be seen from Table 3 that compared with no shock absorption layer, the maximum convergence value of the side wall decreases after the shock absorption layer is set in the tunnel. The maximum convergence value of the side wall of condition 1 is 1.44 mm and that of condition 2 is 1.20 mm, which is 14.90% less than that of condition 1. The maximum convergence value of the side wall of condition 3 ∼ 4 is also less than that of condition 1, but the reduction is not as large as that of condition 2; among them, the reduction of condition 4 is the smallest. With the increase of the thickness of the shock absorption layer, the decrease of the maximum convergence value of the side wall continues to decrease. This phenomenon is due to the influence of the lateral dislocation between the tunnel lining and the surrounding rock. The influence of lateral dislocation will increase with the increase of the thickness of the shock absorption layer. According to the convergence analysis of the side wall, the effect of a 50 mm thick shock absorption layer is the best.
4.2. Principal Stress
According to the simulation results, the maximum and minimum principal stress cloud charts of tunnel lining under various conditions are obtained, as shown in Figures 6 and 7. At the same time, the maximum and minimum principal stresses of the tunnel lining under various conditions (the principal stress is positive in tension and negative in compression) are extracted, and the effect of different thickness shock absorption layers (condition 2 ∼ condition 4) on the maximum and minimum principal stresses is calculated (compared with condition 1). The calculation results are shown in Table 4.

(a)

(b)

(c)

(d)

(a)

(b)

(c)

(d)
It can be seen from Figures 6, 7, and Table 4 that compared with no shock absorption layer, the maximum and minimum principal stress of each condition decreases after setting the shock absorption layer. In condition 1, the maximum value of maximum and minimum principal stresses is 3.231 × 106 Pa and -5.816 × 106 Pa, respectively. After setting the shock absorption layer, in terms of maximum principal stress, the maximum value of the maximum principal stress of condition 2 is larger than the other conditions, which is 2.780 × 106 Pa and 13.96% lower than that of condition 1; the maximum value of the maximum principal stress of condition 3 is smaller than the other conditions, which is 2.756 × 106 Pa and 14.70% lower than that of condition 1. In terms of minimum principal stress, the maximum value of the minimum principal stress of condition 2 is larger than the other conditions, which is -5.258 × 106 Pa and 9.59% lower than that of condition 1; the maximum value of the maximum principal stress of condition 4 is smaller than the other conditions, which is -5.000 × 106 Pa and 14.03% lower than that of condition 1. From the analysis of the maximum and minimum principal stress, it can be seen that the principal stress of condition 3 and condition 4 decreases greatly, so it is better to set the rubber plate with the thickness of 100 mm or 150 mm.
4.3. Maximum Shear Stress
According to the simulation results, the maximum shear stress cloud chart of tunnel lining under various conditions is obtained, as shown in Figure 8. At the same time, the maximum shear stress of the tunnel lining under various conditions is extracted, and the effect of different thickness shock absorption layers (condition 2 ∼ condition 4) on the maximum shear stress is calculated (compared with condition 1). The calculation results are shown in Table 5.

(a)

(b)

(c)

(d)
It can be seen from Figure 8 and Table 5 that the maximum shear stress of each condition occurs in the middle part of the tunnel. Compared with no shock absorption layer, the maximum shear stress decreases after the shock absorption layer is set. The maximum shear stress of condition 1 is the largest, which is 4.781 × 106 Pa. The maximum shear stress of condition 4 is the smallest, which is 4.486 × 106 Pa, which is 6.17% less than that of condition 1, and the antiseismic effect is the most obvious. The maximum shear stress of condition 2 is 4.627 × 106 Pa, and the decrease is the smallest, which is reduced by 3.22%. The maximum shear stress analysis shows that the thicker the shock absorption layer is, the better the antiseismic effect is.
4.4. Safety Factor
The minimum safety factor of each monitoring section under each condition is shown in Figure 9. The effect of setting different thickness shock absorption layers (condition 2 ∼ condition 4) on the minimum safety factor of tunnel lining is calculated (compared with condition 1), and the results are shown in Table 6.

As can be seen from Figure 9 and Table 6, compared with the safety factor of each section of the tunnel without the shock absorption layer, the safety factor of each section of the tunnel is greatly improved after setting different thicknesses of the shock absorption layer, and the increase of the safety factor in condition 3 and condition 4 is basically the same, which is larger than that in condition 2. The minimum safety factor of the four conditions is all located at S3. The minimum safety factor of the tunnel in condition 1 is 2.74, and the minimum safety factor of the tunnel in condition 2 is 3.45, which is 25.91% higher than that in condition 1. The section with the largest increase in condition 2 is S4, which is 32.13% higher than that in condition 1. The minimum safety factor of the tunnel in condition 3 is 4.12, which is 50.36% higher than that in condition 1. The section with the largest increase in condition 3 is S2, which is 54.45% higher than that in condition 1. The minimum safety factor of the tunnel in condition 4 is 4.19, which is 52.92% higher than that in condition 1. The section with the largest increase in condition 4 is S4, which is 53.07% higher than that in condition 1.
When the thickness of the shock absorption layer increases from 50 mm to 100 mm, the safety factor of each section of the tunnel also increases. When the thickness of the shock absorption layer increases from 100 mm to 150 mm, the safety factor of each section of the tunnel does not increase significantly. Therefore, considering the safety and cost, the effect of setting a 100 mm thick shock absorption layer is the best.
5. Conclusions
(1)According to the analysis of the side wall convergence, the side wall convergence of the tunnel lining is reduced after setting different thicknesses of the shock absorption layer. When a 50 mm thick shock absorption layer is set, the reduction of the maximum convergence of the tunnel side wall is the biggest, which is 14.90%. When the thickness of the 150 mm shock absorption layer is set, the reduction of the maximum convergence of the tunnel side wall is the smallest, which is 4.65%. The reason for this phenomenon is that the tunnel lining and surrounding rock have lateral dislocation. The thicker the shock absorption layer is, the greater the influence of lateral dislocation is.(2)From the analysis of principal stress, it can be seen that the principal stress decreases after the construction of shock absorption layer. When the shock absorption layer is 50 mm thick, the reduction of the maximum and minimum principal stresses is the smallest in the three conditions. In terms of the maximum principal stress, the maximum reduction is 14.70% when the shock absorption layer is 100 mm thick. In terms of the minimum principal stress, the maximum reduction is 14.03% when the 150 mm thick shock absorption layer is set. The results show that the antiseismic effect is better when the thickness of the shock absorption layer is 100 mm or 150 mm.(3)From the analysis of the maximum shear stress, it can be seen that the maximum shear stress is reduced after the construction of the shock absorption layer. The maximum shear stress is reduced by 6.17% after setting a 150 mm thick shock absorption layer, which is the largest of the three conditions. The decrease of the maximum shear stress is small when the shock absorption layer is 50 mm and 100 mm.(4)The safety factor of each section of the tunnel increases after the construction of the shock absorption layer, and the increase is the smallest when the shock absorption layer is 50 mm thick. The minimum safety factor of the tunnel with a a100 mm thick shock absorption layer is 4.12, which is increased by 50.36%; the section with the biggest increase in the safety factor of the tunnel is S2, and the biggest increase is 54.45%. The minimum safety factor of the tunnel with a150 mm thick shock absorption layer is 4.19, increased by 52.92%; the section with the largest increase in the safety factor of the tunnel is S4, and the largest increase is 53.07%. The effect of a 100 mm and 150 mm shock absorption layers is basically the same.(5)When the shock absorption layer is set in the tunnel, the displacement and principal stress are controlled, and the safety factor is increased. When the shock absorption layer is 100 mm or 150 mm thick, the maximum and minimum principal stress and the maximum shear stress are significantly reduced, and the reduction range is basically the same; the safety factor also has a similar increase range. Therefore, the two thickness shock absorption layers can meet the seismic requirements of the tunnel. On this basis, considering the construction cost, a 100 mm thick shock absorption layer is selected. In conclusion, the effect of a 100 mm thick shock absorption layer is better than that of 50 mm and 150 mm.Data Availability
The data used to support the findings of this study are available within the article.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.