Research Article  Open Access
A Study on Seismic Isolation of Shield Tunnel Using QuasiStatic Finite Element Method
Abstract
Using a quasistatic method based on an axisymmetric finite element model for seismic response analysis of seismically isolated tunnels, the seismic isolation effect of the isolation layer is studied, and the seismic isolation mechanism of the isolation layer is clarified. The results show that, along the longitudinal direction of the tunnel, the seismic isolation effect is mainly affected by the shear modulus of the isolation material. The smaller the shear modulus is, the more evident the seismic isolation effect is. This is due to the tunnel being isolated from deformation of its peripheral ground through shear deformation of the isolation layer. However, along the transverse direction of the tunnel, the seismic isolation effect is mainly affected by the shear modulus and Poisson’s ratio of the isolation material. When Poisson’s ratio is close to 0.5, a seismic isolation effect is not evident because the tunnel cannot be isolated from deformation of its peripheral ground through compression deformation of the isolation layer. Finally, a seismic isolation system comprising a shield tunnel in which flexible segments are arranged at both ends of an isolation layer is proposed, and it is proved that the seismic isolation system has significant seismic isolation effects both on the longitudinal direction and on the transverse direction.
1. Introduction
Shield tunnels are widely used for water supply, public transport, communication, sewerage, and other infrastructures. Tunnels do not cause selfexcited vibration under earthquakes but rather are controlled by the surrounding ground deformation. Therefore, it is generally recognized that tunnels are not easily affected by earthquakes. In recent years, with the increase in the number of underground structures and the frequent occurrence of seismic damage to underground structures, the antiseismic issue of underground structures is increasingly attracting high attention from seismologists around the world.
As shown in Figure 1, the seismic isolation layer applied to a tunnel is a technology to reduce the seismic response of the tunnel by isolating the tunnel body from deformation of the ground in the outer periphery of the structure. The effectiveness of the seismic isolation layer has been verified by a series of numerical simulations and experiments [1–4]. Seismic response analysis of tunnels is divided into crosssectional seismic response analysis of tunnels and seismic response analysis of tunnels as a whole. At present, there are many methods for seismic response analysis of tunnels; typical ones include the Bessel function approximate solution, the response displacement method, and the finite element dynamic analysis method. The Bessel function approximate solution is mainly suitable for crosssectional seismic response analysis of tunnels, without dynamic analysis of the surrounding ground, and establishment of a new structural model of the tunnel. The approximate solution of the tunnel lining internal force of a circular cross section can be obtained using a theoretical formula. The analytic target and external force are based on ideal assumption. In the case of a sudden change in ground conditions, the approximate analytical Bessel function solution for general ground cannot give the internal force of the tunnel lining correctly. Under this condition, it is necessary to use other analytic methods [5, 6]. The response displacement method is suitable for crosssectional seismic response analysis of tunnels and seismic responses of tunnels as a whole at the same time. The tunnel lining is simulated as a ring beam for the crosssectional seismic analysis or the straight beam for the whole longitudinal seismic analysis. The ground spring, showing the interaction between the tunnel lining and surrounding ground, is placed around the tunnel lining, and the surrounding ground response displacement and shear stress are input into one end of the ground spring. This method is practical. However, due to the complexity of ground conditions, it is difficult to set up the ground spring, and the beamspring model used for seismic design of underground pipelines does not necessarily give a good approximation for tunnels with large diameters, so it has limitations [7, 8]. In the finite element dynamic nonlinear analytic method, it is needed to consider various nonlinear questions, and the modeling process is complicated. It requires very large calculation resources and time, especially in threedimensional analysis. Meanwhile, it is difficult to provide strict dynamic boundary conditions. Hence, it is not suitable for practical seismic isolation design of tunnels [9, 10].
The behavior of a tunnel subjects to deformations imposed by the surrounding ground. At present, the research on seismic isolation of tunnels mainly focuses on tunnel cross sections. However, as a linear underground structure, the key point is to study the seismic isolation of tunnel as a whole. The whole longitudinal behavior of a tunnel subjected to deformations imposed by the surrounding ground can be divided into two types [11]: (1) compressive and tensile deformations along the longitudinal direction (Figure 2(a)) and (2) bending deformation along the transverse direction (Figure 2(b)). In this paper, a quasistatic method based on an axisymmetric finite element model for seismic response analysis of seismically isolated tunnels is used to research the seismic isolation effect and mechanism of seismic isolation layer of a shield tunnel along both the longitudinal and transverse directions.
2. QuasiStatic Method for Seismic Response Analysis of Seismically Isolated Tunnels
2.1. Outline of the QuasiStatic Method
The quasistatic approach adopting an axisymmetric FEM and the response acceleration method is used for practical design of seismically isolated tunnels; this approach eliminates the troublesome processing of boundary conditions related to reflected waves and modeling complex ground and structural conditions.
It is well known that the response acceleration method provides slightly better evaluations than the finite element dynamic analysis approach; that is, the method presents more reliable and safer evaluation results. This method has been adapted to the seismic design of underground ducts, shafts, and rock caverns. In this method, the surrounding ground is modeled by means of finite elements, and accelerations are calculated in a free field. The calculated accelerations are then applied to the finite element model, including an underground structure [12–14].
Figure 3 illustrates a schematic representation of the axisymmetric modeling. In the upper part of the figure, a tunnel is constructed through a soil deposit, while the lower part represents the axisymmetric modeling of the ground and tunnel conditions shown in the upper part; here, the centerline of the tunnel is set as an axis of symmetry. It is essential that special considerations should be taken when modeling the effects of the ground surface, boundary conditions, and seismic load. Accordingly, if a method can be defined to convert the seismic load in the axisymmetric model, thereby equalizing the ground displacements around the tunnel shown in the upper and lower parts of Figure 3, a simplified procedure can be developed to evaluate the 3D interaction effects around a tunnel.
2.2. Method to Convert the Seismic Load
The fundamental theory regarding the method employed to determine the seismic load in the axisymmetric model will be described in this section. The corresponding underground and loading conditions are shown in Figure 4, where H denotes the thickness of the soil deposits, z denotes the coordinate originating from the bedrock, and h_{c} denotes the height from the bedrock to the tunnel center under actual tunnel and ground conditions. In the axisymmetric model, an axis of symmetry is located at the height of the upper outer surface of the actual tunnel, and the outer radius R denotes the distance from the axis to the bedrock. The outer surface of the axisymmetric model corresponding to the outer surface of the tunnel lining is located at the radius r = r_{0}, where r_{0} denotes the outer radius of the tunnel lining. Then, h_{c} and R can be related by h_{c} = R − r_{0}. The seismic load considered here is an inertial force originating from the ground acceleration due to seismic ground motions. The static loading method, in which the inertial force due to an earthquake is loaded statically while ignoring the damping term in the equations of motion, is employed herein.
Figure 5 illustrates a schematic diagram to describe the method of applying an inertial seismic load on the ground. In the proposed model, the seismic load is divided into two components: an inertial force denoted , which is given in equation (1), acting on the ground beneath the tunnel and a concentrated force S, which is given in equation (2), acting on the outer surface of the tunnel body:
The conversion of the seismic load in the axisymmetric model is conducted as shown in Figure 5.
2.2.1. Seismic Load in the Longitudinal Direction
The conversion of the seismic load in the longitudinal direction is conducted in the following manner [15, 16]. The displacement of the ground at z with the unit width is defined as u when the horizontal force is acting at a height z. The displacement at r (=R − z) of a hollow cylindrical disk with a unit thickness, the outer surface of which is constrained, is defined as when the horizontal force is acting at r. Placing the displacements and equal to each other, can be expressed using , , and as follows:
The concentrated force S in equation (2) can also be replaced by as shown below:
Figures 6(a) and 6(b) represent the relation between the concentrated force and ground shear deformation in the actual ground and in the axisymmetric model, respectively; the case of a soil column is shown in Figure 6(a), and the case of a hollow cylindrical disk is shown in Figure 6(b). The shear displacement values at the location r = r_{0} in Figure 6(a) and the shear displacement values at the location z = h_{c} in Figure 6(b) are coincident with each other. However, the shape of the vertical shear displacement distribution in the hollow cylindrical disk in Figure 6(b) is not straight and is widely different from that in the soil column displayed in Figure 6(a). To make the vertical distribution of the shear displacement in the proposed model equal to that in the actual ground, the following method is adopted.
(a)
(b)
When the concentrated load acts on r = r_{0}, the shear displacement u_{c} is
However, as shown in Figure 7, assuming that the acceleration is uniformly distributed in the vertical direction, the inertial force is
When the inertial force acts on a hollow cylindrical disk with a unit thickness ranging from to , the shear displacement produced by the inertial force is
Then, the shear displacement produced at r = r_{0} is
The uniformly distributed acceleration, which can be derived from equation (8), is given by equation the following equation:
In numerical analysis, earthquake loads acting at are calculated by multiplying the mass matrix by the ground acceleration at . Thus, the acceleration at r used in such calculations can be given by equation (10), which is derived from equations (3) and (9); this acceleration denoted used herein is defined as a modified acceleration coefficient:
Finally, by applying the inertial force, which is the product of the mass and acceleration of each node as the seismic load, the displacement generated in the actual ground can be reproduced. In the analysis of the longitudinal direction, the load acts in the longitudinal direction on each hollow cylindrical disk as an axisymmetric load.
2.2.2. Seismic Load in the Transverse Direction
The conversion of the seismic load in the transverse direction is conducted in the following manner. The displacement of the ground at z with the unit width is defined as u when the horizontal force is acting at a height z. The displacement at r (=R − z) of a hollow cylindrical disk with a unit thickness, the outer surface of which is constrained, is defined as when the horizontal force is acting at . Placing the displacements and equal to each other, can be expressed using , , and as follows:
The concentrated force S in equation (2) can also be replaced by as follows:
Figures 8(a) and 8(b) represent the relation between the concentrated force and ground shear deformation in the actual ground and the axisymmetric model, respectively; the case of a soil column is shown in Figure 8(a), and the case of a hollow cylindrical disk is shown in Figure 8(b). The shear displacement values at the location r = r_{0} in Figure 8(b) and the shear displacement values at the location z = h_{c} in Figure 8(a) are coincident with each other. However, the vertical shear displacement distribution in the hollow cylindrical disk displayed in Figure 8(b) is not straight and is widely different from that in the soil column illustrated in Figure 8(a). To make the vertical distribution of the shear displacement in the proposed model equal to that in the actual subsurface, the following method is adopted.
(a)
(b)
When the concentrated load acts on r = r_{0}, the shear displacement u_{c} is
However, as shown in Figure 9, assuming that the acceleration is uniformly distributed in the vertical direction, the inertial force is
When the inertial force acts on a hollow cylindrical disk with a unit thickness ranging from to , the shear displacement produced by the inertial force is
Then, the shear displacement produced at r = r_{0} is
The uniformly distributed acceleration a, which can be derived from equation (16), is given by the following equation:
In numerical analysis, an earthquake load acting at r is calculated by multiplying the mass matrix by the ground acceleration at r. Thus, the acceleration at r used in such calculations can be given by equation (18), which is derived from equations (11) and (17); the acceleration used herein is defined as a modified acceleration coefficient:
Finally, by applying the inertial force, which is the product of the mass and acceleration at each node as the seismic load, the displacement generated in the actual ground can be reproduced in the axisymmetric model. In the analysis of the transverse direction, the load acts in the transverse direction on each hollow cylindrical disk as an asymmetrical load.
3. Seismic Isolation Effect and Mechanism of the Isolation Layer
3.1. Modeling of Shield Tunnels and Method for Determining Calculation Parameters
A shield tunnel is a structure formed by segments, which are fastened by joints. For reasonably showing the shield tunnel characteristics, joint effects should be included when modeling a shield tunnel. An equivalent stiffness beam model of the shield tunnel is often used in the actual shield tunnel design. As shown in Figure 10, based on the equivalent stiffness beam model, the calculation parameters of the shield tunnel finite element in the axisymmetric finite element model are determined.
The equivalent stiffness of the shield tunnel can be calculated using the following equation:where is the equivalent compression stiffness, is the equivalent tension stiffness, is the equivalent bending stiffness, is the segment width, is the crosssectional area of the segment, is the elastic modulus of the segment, is the sum of the tension stiffness of the spring of ring joints, is the moment of inertia of the segment section, and is the angle corresponding to where the neutral axis of segment section exists.
Therefore, the elastic modulus of the finite elements of shield tunnel lining in the axisymmetric finite element model is as follows:where is the compression elastic modulus, is the tension elastic modulus, and is the bending elastic modulus.
3.2. Research on the Seismic Isolation Effect of Isolation Layers
As shown in Figure 11, to study the seismic isolation effect of the isolation layer, seismic response analyses are conducted on cases in which an isolation layer is applied to a shield tunnel buried in irregularly bounded surface soil deposits, and the results of these analyses are compared in this paper. A shield tunnel with an outer diameter of 5.1 m made of reinforced concrete segments with a thickness of 25 cm constructed in the place where ground conditions change sharply is considered. Because the outer diameter of the shield machine was 5.2 m, the thickness of the isolation layer was set to 10 cm, which equals the thickness of the tail void. Table 1 reports the equivalent stiffness of the shield tunnel. Figure 12 shows the discretization and boundary conditions of the axisymmetric finite element mesh. In the axisymmetric model, the mesh is discretized into 9 layers in the circumferential direction. The lining and isolation layer are each modeled by a layer.

(a)
(b)
As shown in Figure 13, the input earthquake motion for Level 2 earthquake motion was chosen from among standard waves used for seismic design of road bridges [17]. Waves obtained by performing amplitude of the EW component at JR Takatori Station were used as Type II earthquake motion for category II ground, which is a motion originated from a nearfault earthquake.
Seismic isolation materials suitable for underground structures should have certain physicochemical properties, and the most important thing to consider is whether the existence of an isolation layer will affect the static stability of underground structures. A softer isolation layer has a more evident seismic isolation effect; however, considering the uneven settlement of ground, Poisson’s ratio of the isolation material should be close to 0.5. When Poisson’s ratio of the isolation material is close to 0.5, uneven settlement of the ground can be avoided effectively, even if the isolation material is very soft [18]. In this paper, silicone material (SISMO), recommended for seismic isolation design of shield tunnels, is used as the isolation material for numerical simulations. Its shear elastic modulus is relatively small, and its Poisson’s ratio is close to 0.5. It meets requirements for grouting material. After being hardened, there is no harmful ingredient, nor does it pollute underground water. Table 2 reports the parameters for the SISMO silicone isolation material.

Onedimensional site response analysis based on the multireflection theory for layered soil (termed as the multireflection analysis) is conducted using the equivalent linear technique in this paper. The site responses as a free field are calculated, and the equivalent seismic rigidity and damping of the layered soil are obtained. The equivalent seismic rigidity of the soil deposits in Table 3 is used in every seismic analysis hereafter. The calculated accelerations are then applied to the axisymmetric finite element model.

3.2.1. Longitudinal Direction
Figures 14(a) and 14(b) show a comparison of the analysis results for the tunnel lining axial strains and stresses. Evidently, the axial peak strain and stress decrease after implementing the isolation layer. Table 4 indicates a decreased rate of axial strain and stress corresponding to different shear modulus. The lower the shear modulus of the isolation material, the better the isolation effect of the isolation layer.
(a)
(b)

3.2.2. Transverse Direction
Figure 15(a) shows a comparison of the analysis results for the tunnel lining axial strains. After the isolation layer is added, the axial peak tension strain decreases slightly, and the axial peak compression strain increases slightly.
(a)
(b)
Figure 15(b) shows a comparison of the analysis results for the tunnel lining axial stresses. After the isolation layer is added, the axial peak tension stress decreases slightly, and the axial peak compression stress almost remains the same. Therefore, when Poisson’s ratio of isolation material is close to 0.5, the seismic isolation effect of isolation layer is not evident in the transverse direction.
3.3. Research on Seismic Isolation Mechanism of Isolation Layers
A spring is used to simulate the interaction between the seismically isolated tunnel and ground. Figures 16(a) and 16(b) show calculation models of the isolation layer spring stiffness along the longitudinal and transverse directions of the tunnel, respectively. The calculation formula is as follows:wherein is the spring stiffness coefficient of the isolation layer along the longitudinal direction of a tunnel, is the spring stiffness coefficient of the isolation layer along the transverse direction of a tunnel, is the isolation layer outer diameter, is the isolation layer inner diameter, is the shear modulus of the isolation layer, is Poisson’s ratio of the isolation layer, and is the thickness of the isolation layer.
(a)
(b)
Figure 17 shows how the shear modulus of isolation material affects the spring stiffness when Poisson’s ratio is 0.48. As is known, the spring stiffness increase for both transverse and longitudinal tunnels with increased shear modulus. Figure 18 shows how Poisson’s ratio affects the spring stiffness when the shear modulus is 0.3 MPa. When the shear modulus is constant, Poisson’s ratio does not affect the spring stiffness of the isolation layer in the longitudinal direction, but the spring stiffness of the isolation layer in the transverse direction increases with increased Poisson’s ratio. When Poisson’s ratio is greater than 0.4, the increasing amplitude of spring stiffness increases sharply with increased Poisson’s ratio.
For isolation material, the shear modulus should be relatively small, and Poisson’s ratio is close to 0.5, avoiding ground settlement caused by the isolation layer. When the shear modulus is constant, Poisson’s ratio does not affect the spring stiffness of the isolation layer along the longitudinal direction of a tunnel. When an isolation material with a small shear modulus is used, the tunnel is isolated from the deformation of its peripheral ground through the shear deformation of the isolation layer. A good isolation effect along the longitudinal direction of a tunnel is achieved. When Poisson’s ratio is close to 0.5, the spring stiffness is still great when an isolation material with small shear modulus is used, verifying that ground settlement can be avoided effectively when Poisson’s ratio is close to 0.5. However, the tunnel cannot be isolated from deformation of its peripheral ground through the compression deformation of the isolation layer. Therefore, a good isolation effect along the transverse direction of a tunnel cannot be achieved.
4. Research on Seismic Isolation System of Shield Tunnel
Based on the above concept, the material used for the isolation layer in shield tunnels should have relatively low shear modulus, and its Poisson’s ratio should be approximately 0.5. However, the isolation layer cannot provide good isolation effect along the transverse direction of a tunnel. Considering that an earthquake would lead to concentration of stress at shaft joints and the positions at which ground conditions sharply change, flexible segments are correspondingly employed to increase the flexibility of the shield tunnel and improve its bending deformation capacity. However, there exist some difficulties in the layout of flexible segments owing to the uncertainty regarding the strain distribution in a shield tunnel when an earthquake occurs. Thus, a shield tunnel isolation system combining an isolation layer and flexible segments is proposed in this paper. As is shown in Figure 19, flexible segments are set at both ends of an isolation layer.
The seismic responses of the shield tunnel along the longitudinal direction of the tunnel are shown in Figure 20(a). It is observed that when the seismic isolation system is built, the peak tension stress decreases to 47.9% of its original value and the stress distribution becomes smoother along the axial direction of the tunnel. On the contrary, the seismic responses of the shield tunnel along the transverse direction of the tunnel are represented in Figure 20(b). It is observed that when the axial peak tension stress decreases to 16.7% of its original value, the peak compression stress drops to 47.5% of its original value, and the stress distribution becomes smoother along the axial direction of the tunnel. Therefore, the proposed isolation system could have obvious effects on both longitudinal and transverse directions of a shield tunnel.
(a)
(b)
5. Conclusions
In this paper, a quasistatic method based on an axisymmetric finite element model for seismic response analysis of seismically isolated tunnels is used to research the seismic isolation effect and mechanism of seismic isolation layer of a shield tunnel along both longitudinal and transverse directions. The conclusions derived in this paper can be summarized as follows:(1)Along the longitudinal direction of a tunnel, the seismic isolation effect of the elastic isolation layer is mainly affected by the shear modulus. The smaller the shear modulus of the isolation material, the better the isolation effect. This result occurs because the tunnel can be isolated from deformation of its peripheral ground through shear deformation of the isolation layer.(2)Along the transverse direction of a tunnel, the seismic isolation effect of the elastic isolation layer is mainly affected by the shear modulus and Poisson’s ratio. When Poisson’s ratio is approximately 0.5, the seismic isolation effect is not obvious. This result occurs because the tunnel cannot be isolated from the deformation of its peripheral ground through the compression deformation of the isolation layer.(3)A seismic isolation system for a shield tunnel in which flexible segments are arranged at both ends of an isolation layer is proposed. Based on the seismic response analysis, it is concluded that the system could have obvious effects on both the longitudinal and transverse directions of a tunnel.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The authors acknowledge the financial support received from the National Program on Key Basic Research Project of China (973 Program) (Grant no. 2015CB057906) and the National Key Research and Development Program of China (Grant no. 2018YFF01014204).
References
 Z. Y. Chen and H. Shen, “Dynamic centrifuge tests on isolation mechanism of tunnels subjected to seismic shaking,” Tunnelling and Underground Space Technology, vol. 42, pp. 67–77, 2014. View at: Publisher Site  Google Scholar
 W. S. Zhao, W. Z. Chen, X. J. Tan, and S. Huang, “Study on foamed concrete used as seismic isolation material for tunnels in rock,” Materials Research Innovations, vol. 17, no. 7, pp. 465–472, 2014. View at: Publisher Site  Google Scholar
 T. Suzuki and T. Tanaka, “Research and development on the seismic isolation system applied to urban tunnels (part1: development of seismic isolation materials and construction methods),” in Proceedings of the 7th U.S.Japan Workshop on Earthquake Disaster Prevention for Lifeline Systems, pp. 447–459, Seattle, WA, USA, November 1997. View at: Google Scholar
 T. Tanaka and T. Suzuki, “Research and development on the seismic isolation system applied to urban tunnels (part2: effects of seismic isolation and seismic design),” in Proceedings of the 7th U.S.Japan Workshop on Earthquake Disaster Prevention for Lifeline Systems, pp. 461–475, Seattle, WA, USA, November 1997. View at: Google Scholar
 W. Zhao, W. Chen, and D. Yang, “Interaction between strengthdening and isolation layers for tunnels in rock subjected to SH waves,” Tunnelling and Underground Space Technology, vol. 79, pp. 121–133, 2018. View at: Publisher Site  Google Scholar
 M. N. Wang and G. Y. Cui, “Study of the mechanism of shock absorption layer in the supporting system of tunnels in highly seismic areas,” China Civil Engineering Journal, vol. 44, no. 8, pp. 126–131, 2011. View at: Google Scholar
 K. Kawashima, Seismic Design of Underground Structures, Kajima Press, Tokyo, Japan, 1994.
 T. Suzuki, “Comparison between finite elemnt analysis and seismic displacement method on a transmission coefficient of seismic axial strain of shielddriven tunnels,” Journal of Japan Society of Civil Engineers, Ser. A1 (Structural Engineering & Earthquake Engineering (SE/EE)), vol. 65, no. 1, pp. 236–243, 2009. View at: Publisher Site  Google Scholar
 S. Huang, W. Z. Chen, J. P. Yang et al., “Research on earthquakeinduced dynamic responses and aseismic measures for underground engineering,” Chinese Journal of Rock Mechanics and Engineering, vol. 28, no. 3, pp. 483–490, 2009. View at: Google Scholar
 W. S. Zhao, W. Z. Chen, S. S. Ma et al., “Isolation effect of foamed concrete layer on seismic responses of tunnel,” Rock and Soil Mechanics, vol. 39, no. 2, pp. 1027–1036, 2018. View at: Google Scholar
 Y. M. A. Hashash, J. J. Hook, B. Schmidt, and J. IChiang Yao, “Seismic design and analysis of underground structures,” Tunnelling and Underground Space Technology, vol. 16, no. 4, pp. 247–293, 2001. View at: Publisher Site  Google Scholar
 I. Katayama, M. Adachi, Y. Shimada et al., “Proposal of practical quasidynamic analysis method for underground structures named response acceleration method,” Proceedings of the 40th Annual Meeting of JSCE, vol. 369, pp. 737738, 1985. View at: Google Scholar
 I. Katayama, Studies on Fundamental Problems in Seismic Design Analyses of Critical Structures and Facilities, Kyoto University, Kyoto, Japan, 1990, Ph.D. thesis.
 A. Tateishi, “A study on seismic analysis methods in the cross section of underground structures using static finite element method,” Structural Engineering/Earthquake Engineering, vol. 22, no. 1, pp. 41s–54s, 2005. View at: Publisher Site  Google Scholar
 T. Suzuki and M. Maruyama, “The axisymmetric finite element model developed as a measure to evaluate earthquake response of seismically isolated tunnels,” in Proceedings of the 10th Japan Earthquake Engineering Symposium, pp. 2003–2008, Tokyo, Japan, 1998. View at: Google Scholar
 T. Suzuki, “The axisymmetric finite element model developed as a measure to evaluate earthquake response of seismically isolated tunnels,” in Proceedings of the 12th World Conference on Earthquake Engineering, Auckland, New Zealand, September 2000. View at: Google Scholar
 Japan Road Association, Specifications of Highway Bridges (PartV: Seismic Design), Maruzen Press, Tokyo, Japan, 2012.
 T. Suzuki and C. Tamura, “Proposal of a seismic isolation structure for urban tunnels and the method to evaluate its isolation effects,” Doboku Gakkai Ronbunshu, vol. 525/I33, no. 525, pp. 275–285, 1995. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2019 Can Li et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.