Advancements in Analysis and Design of Protective Structures against Extreme Loadings
View this Special IssueResearch Article  Open Access
Xiaojie Zhou, Qinghua Liang, Zhongxian Liu, Ying He, "IBIEM Analysis of Dynamic Response of a Shallowly Buried Lined Tunnel Based on ViscousSlip Interface Model", Advances in Civil Engineering, vol. 2019, Article ID 1025483, 14 pages, 2019. https://doi.org/10.1155/2019/1025483
IBIEM Analysis of Dynamic Response of a Shallowly Buried Lined Tunnel Based on ViscousSlip Interface Model
Abstract
A viscousslip interface model is proposed to simulate the contact state between a tunnel lining structure and the surrounding rock. The boundary integral equation method is adopted to solve the scattering of the plane SV wave by a tunnel lining in an elastic halfspace. We place special emphasis on the dynamic stress concentration of the lining and the amplification effect on the surface displacement near the tunnel. Scattered waves in the lining and halfspace are constructed using the fictitious wave sources close to the lining surfaces based on Green’s functions of cylindrical expansion and the shear wave source. The magnitudes of the fictitious wave sources are determined by viscousslip boundary conditions, and then the total response is obtained by superposition of the free and scattered fields. The slip stiffness and viscosity coefficients at the liningsurrounding rock interface have a significant influence on the dynamic stress distribution and the nearby surface displacement response in the tunnel lining. Their influence is controlled by the incident wave frequency and angle. The hoop stress increases gradually in the inner wall of the lining as sliding stiffness increases under a lowfrequency incident wave. In the highfrequency resonance frequency band, where incident wave frequency is consistent with the natural frequency of the soil column above the tunnel, the dynamic stress concentration effect is more significant when it is smaller. The dynamic stress concentration factor inside the lining decreases gradually as the viscosity coefficient increases. The spatial distribution and the displacement amplitudes of surface displacement near the tunnel change as incident wave frequency and angle increase. The effective dynamic analysis of the underground structure under an actual strong dynamic load should consider the slip effect at the liningsurrounding rock interface.
1. Introduction
Analyses of seismic damage incurred by disasters such as the Kobe earthquake, ChiChi earthquake, and Wenchuan earthquake have shown that underground structures might be severely damaged during strong earthquakes, resulting in massive economic and societal losses [1–5]. As the scope and scale of modern underground structures continually increase, the seismic design grows increasingly complex. In theory, wave scattering and dynamic stress concentration effects should be considered for large underground structures during the seismic wave propagation process. To this effect, it is of great significance to study the seismic response and hazards feature of underground structures.
In general, the calculation methods include the analytical method [6–9], finite element method [10–12], finite difference method [13–15], boundary element method [16–19], boundary integral equation method [20–23], etc. It is worth mentioning that all of these studies assumed that the lining and surrounding rock are completely bonded. In actuality, however, there are varying degrees of slip between the tunnel and surrounding rock—especially in the case of intense dynamic loading. Yi et al. [24] studied the dynamic response of a tunnel lining with a sliding interface under an incident wave based on an interface contact model. Fang and Jin [25, 26] proposed a viscoelastic interface model to solve the dynamic response of a tunnel lining under different interface stiffness and viscosity coefficients with incident and SV waves. Ping et al. [27] calculated the maximum moment and axial force of a circular shield tunnel under interface slip and noslip states based on an equivalent stiffness circle.
It should be noted that the abovementioned research works based on interface contact models were mainly limited to deepburied tunnels, while the response of shallowburied tunnels differs significantly from that of deepburied ones [28]. Yi et al. [29] presented an analytical solution to the outofplane dynamic response of a shallow tunnel lining under the action of a plane SH wave including interface contact stiffness, incident angle, wave frequency, and tunnel depth as these factors affect the dynamic stress concentration of the lining. Fang et al. [30] investigated a lined tunnel in a semiinfinite alluvial valley with an elasticslip interface and analyzed the dynamic stress distribution around the circular tunnel subjected to SH waves.
However, until now, few studies have explored the seismic response of shallow tunnels under incident SV waves with the interfaceslipping model due to the complexity of multimode coupling and hybrid boundary conditions. We used an indirect boundary integral equation method (IBIEM) to solve the scattering of the plane SV wave by a tunnel lining in a halfspace based on the viscousslip interface model. This method had been used effectively to solve the dynamic response of the tunnel structure [31, 32].
This study aims to investigate the dynamic stress concentration effect of the tunnel lining and the surface displacement amplification near the tunnel with a viscousslip interface. We assessed the influence of parameters such as incident wave frequency and angle, viscousslip interface stiffness, and viscosity coefficient on the overall dynamic response of the lining and surrounding rock. This study can provide a theoretical basis for the seismic design of actual underground engineering structures under intense dynamic loads.
2. Calculation Model
As shown in Figure 1, the elastic halfspace contains an infinitely long tunnel lining. Both the lining and the halfspace are linearly isotropic homogeneous elastic media. The viscousslip state between the tunnel and the surrounding rock can be assumed to consist of a series of linear springs and dampers. The parameters related to the halfspace and the tunnel are shown in Table 1.

Let the buried depth of the tunnel be , the inner and outer radii of the lining be and , and the inner and outer boundary surfaces of the lining be and , respectively. Assume that the plane SV wave is incident at an angle from the halfspace.
3. Calculation Method
In this study, we considered the cylindrical shear wave source in the halfspace as the fundamental solution. The indirect boundary integral equation method and viscousslip boundary condition were used to solve the scattering of plane waves by tunnel linings and the dynamic response around the tunnel lining [20].
3.1. Wave Field Analysis
The total halfspace wave field can be viewed as a superposition of a halfspace free field (without tunnel linings) and a scattering field. We first carry out a free field analysis. The shear wave potential function in the halfspace is denoted as (plane strain state), and the plane SV wave with circular frequency is incident at angle . In the Cartesian coordinate system, the incident SV wave potential function can be expressed as follows:
For the sake of simplicity, the time factor is omitted here. Incident plane SV waves generate reflected waves and SV waves on the surface of the halfspace. Their specific expressions are shown in [16].
Scattered fields are generated in the halfspace of a lined tunnel and in the interior of the lining. The fields can be determined by superimposing all the expansion wave and shear wave sources on the virtual wave source surface inside and outside the lining, respectively. Assuming that the scattered field in the halfspace is generated by the virtual source face , the displacement and stress in the halfspace are as follows:where and . and correspond to the density of and SV wave sources at the position on the virtual wave source surface , respectively. and denote Green’s function of displacement and stress in the elastic halfspace, respectively (lower scale index = 1 or 2 corresponding to and SV wave sources ), and function automatically satisfies the wave equation and surface boundary conditions.
The internal scattering field of the lining is obtained by superimposing the action of all the expansion wave sources and shear wave sources on the virtual wave source surfaces and . The internal displacement and stress of the lining are expressed as follows:where , , and . and , respectively, correspond to the density of and SV wave sources at the position on the virtual wave source surface . and correspond to the density of and SV wave sources at the position on the virtual wave source surface . and indicate the displacement and stress Green’s functions in the lining, respectively.
The total displacement and stress fields in the halfspace are obtained by superposition of the scattered wave field and the free field in the halfspace. The internal reactions of the lining are all generated by the scattered fields within it.
3.2. Boundary Conditions and Solutions
We built a viscousslip interface model to determine the influence of the interface effect on the dynamic response. The lining and the halfspace were connected by linear springs and dampers (Figure 1). Spring and damper parameters are represented by stiffness and viscosity coefficients. In this model, the boundary conditions at the interface between the lining and the halfspace are as follows:where superscripts and correspond to the halfspace and the lining, respectively; and are the normal and tangential stiffness coefficients of the viscousslip boundary; and and are the normal and tangential viscosity coefficients of the viscousslip boundary, respectively. To secure a numerical solution to the problem, we first discretely separate the inner and outer surfaces of the lining and the virtual wave source surfaces , , and . The number of discrete points on the inner and outer surfaces of the lining is set to , and the number of discrete points on the virtual wave source surfaces , , and is . The scattered displacement and stress fields in the halfspace can be expressed as follows:where and ( and ). and are the density of wave and SV wave sources at the th discrete point on the virtual source surface . Similarly, the scattering field inside the lining can be constructed from discrete wave sources on and . The scattering displacement and stress fields arewhere , , and (, , and ).
A linear system of equations can be obtained by synthesizing the above formulas:where , , and are the displacement Green’s influence matrix of the discrete points on the outer surface of the lining of the discrete wave source on , , and , respectively; , , and are boundary displacement matrices obtained from boundary conditions; , , and are the stress Green’s influence matrix of the discrete points on the outer surface of the lining of discrete source points on , , and , respectively; and correspond to the stressgreeninfluenced matrix of the discrete points on the inner surface of the lining of the source points of and ; , , and are the virtual wave source density vectors (to be determined) on , , and , respectively; and is the free field vector. System (7) can be solved using the leastsquares method. The virtual wave source density is obtained, and then the scattered field is obtained. The total wave field can be obtained by superimposing the scattering field and the free field, and then the displacement and stress of the halfspace and any point in the lining can be calculated.
4. Numerical Example and Validation
In this section, the ratio of the buried depth of the tunnel to the inner radius of the lining is ; the ratio of the inner and outer radii of the lining is . Poisson’s ratios of the halfspace and the lining material are both 0.25. The material’s hysteretic damping ratio is 0.001, and the density ratio of the halfspace and lining material is . The ratio of the shear wave speed in the halfspace and the lining is . We define as the dimensionless incident frequency . represents the ratio of the inner diameter of the tunnel to the halfspace shear wavelength.
We define the dimensionless dynamic stress concentration factor (DSCF) as . DSCF represents the absolute value of the ratio of the circumferential stress of the lining hole to the amplitude of the incident wave stress in the halfspace; is the dimensionless stiffness factor of the viscousslip interface, and ; and is the dimensionless viscosity factor of the viscousslip interface, and . In this paper, the dimensionless stiffness factor is 1, 5, 10, and 20, respectively. The dimensionless viscosity factor is 0. We found that the interface is essentially in a noslip state when exceeds 20, so when the stiffness factor exceeds 20, we specify the interface is the noslip state.
In order to verify the correctness of the IBIEM method, we first let be 100 and be 0, which is equivalent to the state of consolidation without slip, and then compare it with the state of no slip in [20]. To more fully prove the accuracy, we choose the spectrum analysis of the different locations of the lining and compare them with [20]. The calculation positions are at = 90° (the top of the lining), 45°, 0°, −45°, and −90° (the bottom of the lining), respectively, as shown in Figure 2. The DSCF of = 90° (the top of the lining) and −90° (the bottom of the lining) is 0, so it is not shown in Figure 2. As shown in Figure 2, the results calculated by the IBIEM method are quite consistent with the results in [20]. This proves the correctness of our calculation method.
5. Numerical Analysis
5.1. DSCF of Inner and Outer Lining Wall Surfaces under SingleFrequency Incident Wave
5.1.1. Influence of Stiffness Factor
Figures 3–6 show the distribution of DSCF in the inner and outer wall of a rigid lining under an incident SV wave. Among them, the dimensionless incident frequency is 0.25, 0.5, 1, and 2, respectively. The incidence angle is 0° and 30°, respectively;
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
As shown in Figures 3–6, under an incident plane SV wave, the DSCF distribution curves of inner and outer wall surfaces of the lining are similar. The inner wall stress amplitude is considerably larger than that of the outer wall. When the SV wave is incident at low frequency ( = 0.25, 0.5), the DSCF of the inner lining wall surface increases as the stiffness factor increases while the DSCF of the outer wall decreases as increases. When and , is either 1, 5, 10, or 20; the DSCF of the outer wall surface is, respectively, 23.79, 23.04, 21.45, and 19.91.
When the SV wave is obliquely incident, the increase and decrease amplitudes of the DSCF of the inner and outer wall surfaces is smaller than the normal incidence. When , , and , the DSCF of the inner wall surface is 51.71. It is 58.39 when = 20 marking an increase in amplitude of 13%. When = 1, the DSCF of the outer wall surface is 37.21, corresponding to 16.26 when = 20, and the amplitude of increase is 56%. When , , and , the DSCF of the inner wall surface is 46.74 and it corresponds to 50.34 when = 20, and the increase in amplitude is 8%. When = 1, the DSCF of the outer wall surface is 28.64, corresponding to 20.15 when = 20 at a decrease in amplitude of 30%.
When SV wave is incident with frequency , the DSCF of the inner and outer wall of the lining decreases with the increase of , and the dynamic stress concentration on the inner and outer wall surfaces is very significant when . If the SV wave is obliquely incident, the DSCF of the inner wall surface is to 69.60 when , which is 2.1 times as much as the peak 32.93 under . The peak DSCF of the outer wall surface is 60.36, which is 2.8 times as much as the corresponding peak value 21.92 under . The SV wave is inclined at an angle of 30°; the peak DSCF on the inner wall of the lining is 67.78 when , which is 2.3 times as much as the peak 29.62 under . Section 3.2 discusses this phenomenon in detail.
When the SV wave is incident at high frequency (), the DSCF curve on the inner and outer wall surfaces of the lining oscillates very sharply along the circumference of the hole. There is no obvious relationship between the vibration regularity and . When the SV wave is inclined at an incidence of 30°s and , the DSCF of the inner and outer wall surfaces is minimal: the peak value of the inner wall falls to 12.64. When , the peak values of the hoop stress of the inner and outer wall surfaces are 27.50 and 23.26, respectively.
5.1.2. Influence of Viscosity Factor
Figure 7 shows the distribution of DSCF of the outer lining wall surface under an incident SV wave. The dimensionless incident frequency is 0.25, 0.5, and 1, respectively; the incident angle is 0° and 30°, respectively. The dimensionless stiffness factor , and the viscosity factor is 1, 10, or 100.
(a)
(b)
(c)
(d)
(e)
(f)
As shown in Figure 7, the DSCF of the outer wall of the lining decreases gradually as interfacial viscosity factor increases. However, the influence of the viscosity factor on the circumferential stress distribution of the lining gradually weakens as incident wave frequency increases. For example, when the SV wave is of low frequency () at normal incidence, the DSCF of the outer wall is 27.32 when , that is, 2.3 times the corresponding value of 11.90 when and 2.6 times the corresponding value of 10.43 when . The SV wave is under at normal incidence. The DSCF of the outer wall is 25.56 when , 1.35 times the corresponding value of 18.93 when , and 1.39 times the corresponding value of 18.43 when .
Compared to the normal incidence, the influence of the viscosity factor on the amplitude distribution of the circumferential stress of the lining weakens at oblique incidence. When the SV wave is incident at an angle of 30°, the spatial distribution of the circumferential stress of the lining is relatively gentle along the circumference of the hole. In the lowfrequency region (), when , the DSCF of the outer wall is 23.99, which is equal to 1.3 times the value of 18.23 when . When and , the DSCF of the outer wall of the lining is 27.12, which is 1.1 times the corresponding value of 23.60 when .
5.2. Lining Internal DSCF Spectrum Analysis
Figures 8 and 9 show the influence of the dimensionless stiffness factor on the DSCF at different positions of the inner and outer walls of the tunnel lining when the SV wave is perpendicularly incident in the spectrum state, and the viscosity factor is 0. Calculation positions were same as in Section 3, and the DSCF of = 90° (the top of the lining) and −90° (the bottom of the lining) is not shown in Figures 8 and 9. Similar to our observations in the model without slipping, the stress of the lining is sensitive to changes in frequency. The spectrum curve also has obvious peaks and troughs. The peak frequency here corresponds to the natural frequency in the soil column above the lining [33].
(a)
(b)
(c)
(a)
(b)
(c)
As shown in Figure 8, under an SV wave with normal incidence, the interface slip stiffness factor significantly affects the internal stress spectrum of the lining. When , the stress amplification effect in the resonance reaction section is particularly obvious. At the firstorder resonance frequency (), the stress peak at reaches 70.0; and at the resonance frequency, the stress at has a peak of 77.2. The stress peak in the resonance band gradually decreases as the stiffness factor gradually increases. For example, when (near the noslip state), the stress peaks at the first two resonance frequencies at are approximately 63.4 and 34.0. When the stiffness factor is small, the restraining effect of the lining on the upper soil column weakens and the response amplitude of the liningupper soil column system in the resonant state is rather large, causing an increase in the corresponding stress amplitude.
The resonance frequency point also is offset to a certain extent as the sliding stiffness factor increases. When is 1, 5, and 10, the firstorder resonance frequency is 0.16, 0.18, and 0.20, respectively. This is due to the fact that the overall stiffness of the soil column above the lining increases as slip stiffness factor increases, and the stress spectrums of and are similar. According to the spatial distribution, the stress spectrum curves at several typical observation points markedly differ. If the stress peak at occurs at the firstorder natural frequency, the stress peak at occurs at the secondorder natural frequency.
As shown in Figure 9, the frequency variation rule of the DSCF spectrum curve of the outer wall is similar to that of the inner wall, but the outer wall DSCF spectrum is generally smaller than the inner wall. When is 1, 5, 10, and 20, the peaks of the inner wall DSCF reach 77.24, 65.65, 64.05, and 63.44 while the peaks of the outer wall DSCF are 64.14, 29.14, 24.87, and 22.33, respectively.
Figure 10 shows the DSCF spectrum of the outer wall surface of the lining based on the influence of the interfacial viscosity factor , where the peak of the DSCF spectrum of the outer wall decreases as the viscosity factor increases. When is 1, 10, and 100, the DSCF peaks are 30.43, 21.10, and 20.57, respectively. Increase in the viscosity factor causes greater energy loss during the sliding process, which in turn causes the DSCF of the lining surface to attenuate.
(a)
(b)
(c)
5.3. Displacement Amplitude of Ground Surface
Figures 11 and 12 show the surface displacement amplitude distribution above the tunnel lining asaffected by interface slip stiffness factor under an incident SV wave. The surface displacement amplitude in the figure was standardized according to the displacement amplitude of the incident wave. The dimensionless incident frequency is 0.25 and 0.5; the incident angle is 0° and 30°, and the dimensionless slip stiffness factor is 1, 2, 10, and 20. The interface viscosity factor is 0.
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
Figures 11 and 12 show that when the SV wave has low frequency () at normal incidence, the spatial distribution of surface displacement is basically consistent under different slip stiffness. The variations in the stiffness factor have a considerable influence on the horizontal and vertical displacement amplitudes of the ground surface near the lining. The standard amplitude of the horizontal displacement above the tunnel lining increases, while the standard amplitude as increases. The surface horizontal displacement amplitudes at = 1, 5, 10, and 20 were 2.04, 2.23, 2.29, and 2.32, respectively, while the vertical displacement amplitudes were 1.71, 1.41, 1.27, and 1.16, respectively.
When the SV wave is incident at an angle of 30° and at a low frequency (), any change in slip stiffness factor has little effect on the horizontal displacement above the tunnel lining surface but does influence the spatial distribution and amplitude of the vertical displacement. The influence of gradually increases as the incident frequency of the SV wave increases () while the spatial distribution and amplitude of the lining’s surface displacement also markedly change. The 30° oblique incidence has more significant effects than the normal incidence. Take the vertical incidence of the SV wave (Figure 12(a)) as an example: at the position of the lining directly above the surface (i.e., ), the horizontal displacement amplitude of the surface decreases as increases. The amplitude is 1.47 when = 1, and the corresponding values for = 5, 10, and 20 are 1.16, 1.05, and 0.98, respectively. Just above the two sides, the horizontal displacement amplitude increases as increases.
6. Conclusion
The boundary integral equation method was applied to solve the seismic response of a tunnel lining in an elastic halfspace under incident plane SV waves based on a viscousslip interface model. The effects of key factors such as incident wave frequency and angle, interface slip stiffness, and interfacial viscosity coefficient on the dynamic stress response of the tunnel lining in elastic halfspace and the surface displacement near the tunnel lining were analyzed. Main conclusions can be summarized as follows:(1)The interface slip stiffness factor significantly affects the dynamic stress distribution of the tunnel lining; the response characteristics are controlled by the incident wave frequency. When the slip stiffness is small, the internal stress of the lining changes very radically along the space around the hole, and the spatial distribution of the dynamic stress is highly complex. When the slip stiffness is large (), the dynamic response is close to that of the nonslip model. Under an incident lowfrequency wave, increase in interface slip stiffness causes a gradual increase in the circumferential stress of the inner wall. The dynamic stress concentration is more significant in the noslip state than the slip state. When the SV wave is incident with (close to the highfrequency resonant frequency band), the dynamic stress concentration effect inside the lining is very significant when the interface slip stiffness coefficient is small (). Under highfrequency wave incidence (), the influence of slip stiffness coefficient on the dynamic stress of the lining is more complex and spatial oscillation of the dynamic stress is more severe.(2)The viscosity factor of the viscousslip interface also has a significant influence on the dynamic stress distribution of the tunnel lining. As the viscosity factor increases, the DSCF of the lining outer wall decreases gradually; however, this effect gradually weakens as the incident wave frequency increases. The influence of the viscosity coefficient under a normally incident plane SV wave is greater than that under a wave with oblique incidence.(3)The interface slip stiffness factor has a significant effect on the DSCF spectrum characteristics of the tunnel lining surface. When the slip stiffness is small (), the dynamic stress amplification effect in the highfrequency resonance reaction section is more obvious. The stress peak in the resonance band gradually decreases as the slip stiffness increases. The resonance frequency point is also offset to a certain extent as slip stiffness increases.(4)When the SV wave is incident at a low frequency, the spatial distribution of surface displacements above the lining under different slip stiffness is essentially constant, but the displacement amplitudes are quite different. Any increase in SV wave incident frequency and incidence angle has a significant effect on both the spatial distribution of surface displacement near the lining and the amplitude of said surface displacement.
In this study, we analyzed only the 2D seismic response of a shallow lined tunnel based on the viscousslip interface model of a uniform space in halfspace. Similar seismic response analyses of uneven sites and 3D tunnels merit further research.
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The research described in this paper was financially supported by the National Natural Science Foundation of China under grant nos. 51678389 and 51678390 and the National Basic Research Program of China under grant no. 2015CB058002.
References
 Ö. Aydan, Y. Ohta, M. Geniş, N. Tokashiki, and K. Ohkubo, “Response and stability of underground structures in rock mass during earthquakes,” Rock Mechanics and Rock Engineering, vol. 43, no. 6, pp. 857–875, 2010. View at: Publisher Site  Google Scholar
 Y. M. A. Hashash, J. J. Hook, B. Schmidt, and J. I.C. 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
 K. Uenishi and S. Sakurai, “Characteristic of the vertical seismic waves associated with the 1995 hyogoken nanbu (kobe), Japan earthquake estimated from the failure of the daikai underground station,” Earthquake Engineering and Structural Dynamics, vol. 29, no. 6, pp. 813–821, 2015. View at: Publisher Site  Google Scholar
 I. Arango, “Theme paper: earthquake engineering for tunnels and underground structures,” in Proceedings of a Case History. Geotechnical Earthquake Engineering and Soil Dynamics Congress IV, pp. 1–34, Sacramento, CA, USA, May 2011. View at: Google Scholar
 Y. M. A. Hashash, D. Park, and J. I.C. Yao, “Ovaling deformations of circular tunnels under seismic loading, an update on seismic design and analysis of underground structures,” Tunnelling and Underground Space Technology, vol. 20, no. 5, pp. 435–441, 2005. View at: Publisher Site  Google Scholar
 K. A. Kuo, H. E. M. Hunt, and M. F. M. Hussein, “The effect of a twin tunnel on the propagation of groundborne vibration from an underground railway,” Journal of Sound and Vibration, vol. 330, no. 25, pp. 6203–6222, 2012. View at: Publisher Site  Google Scholar
 J. A. Forrest and H. E. M. Hunt, “A threedimensional tunnel model for calculation of traininduced ground vibration,” Journal of Sound and Vibration, vol. 294, no. 45, pp. 678–705, 2006. View at: Publisher Site  Google Scholar
 Q. Liu and R. Wang, “Dynamic response of twin closelyspaced circular tunnels to harmonic plane waves in a full space,” Tunnelling and Underground Space Technology, vol. 32, no. 6, pp. 212–220, 2012. View at: Publisher Site  Google Scholar
 Y. Gao, D. Dai, N. Zhang et al., “Scattering of plane and cylindrical SH waves by a horseshoe shaped cavity,” Journal of Earthquake and Tsunami, vol. 11, no. 2, Article ID 1650011, 2016. View at: Publisher Site  Google Scholar
 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
 H. Huo, A. Bobet, G. Fernández, and J. Ramírez, “Load transfer mechanisms between underground structure and surrounding ground: evaluation of the failure of the daikai station,” Journal of Geotechnical and Geoenvironmental Engineering, vol. 131, no. 12, pp. 1522–1533, 2005. View at: Publisher Site  Google Scholar
 G. P. Kouretzis, G. D. Bouckovalas, and C. J. Gantes, “3d shell analysis of cylindrical underground structures under seismic shear (s) wave action,” Soil Dynamics and Earthquake Engineering, vol. 26, no. 10, pp. 909–921, 2006. View at: Publisher Site  Google Scholar
 H. Kawabe and K. Kamae, Long Period Ground Motion Prediction of Linked Tonankai and Nankai Subduction Earthquakes Using 3D Finite Difference Method, AGU Fall Meeting, San Francisco, CA, USA, 2005.
 S. Aoi and H. Fujiwara, “3D finite difference method using discontinuous grids,” Bulletin of the Seismological Society of America, vol. 89, no. 4, pp. 918–930, 1999. View at: Google Scholar
 A. Salemi, R. Mikaeil, and S. S. Haghshenas, “Integration of finite difference method and genetic algorithm to seismic analysis of circular shallow tunnels (case study: tabriz urban railway tunnels),” KSCE Journal of Civil Engineering, vol. 22, no. 5, pp. 1978–1990, 2017. View at: Publisher Site  Google Scholar
 J. E. Luco and F. C. P. D. Barros, “Dynamic displacements and stresses in the vicinity of a cylindrical cavity embedded in a half space,” Earthquake Engineering and Structural Dynamics, vol. 23, no. 3, pp. 321–340, 2010. View at: Publisher Site  Google Scholar
 A. A. Stamos and D. E. Beskos, “3D seismic response analysis of long lined tunnels in halfspace,” Soil Dynamics and Earthquake Engineering, vol. 15, no. 2, pp. 111–118, 1996. View at: Publisher Site  Google Scholar
 H. Alielahi, M. Kamalian, and M. Adampira, “Seismic ground amplification by unlined tunnels subjected to vertically propagating SV and P waves using BEM,” Soil Dynamics and Earthquake Engineering, vol. 71, pp. 63–79, 2015. View at: Publisher Site  Google Scholar
 H. Alielahi and M. Adampira, “Seismic effects of twodimensional subsurface cavity on the ground motion by bem: amplification patterns and engineering applications,” International Journal of Civil Engineering, vol. 14, no. 4, pp. 233–251, 2016. View at: Publisher Site  Google Scholar
 Z. X. Liu, J. W. Liang, and H. Zhang, “Scattering of plane P and SV waves by a lined tunnel in elastic half space (i): method,” Journal of Natural Disasters, vol. 19, no. 4, pp. 71–76, 2010. View at: Google Scholar
 M. Bouchon, C. A. Schultz, and M. N. ToksÖ, “A fast implementation of boundary integral equation methods to calculate the propagation of seismic waves in laterally varying layered media,” Bulletin of the Seismological Society of America, vol. 85, no. 6, pp. 1679–1687, 1995. View at: Google Scholar
 M. Bouchon, M. Campillo, and S. Gaffet, “A boundary integral equationdiscrete wavenumber representation method to study wave propagation in multilayered media having irregular interfaces,” Geophysics, vol. 54, no. 9, pp. 1134–1140, 2012. View at: Publisher Site  Google Scholar
 J. Liang and Z. Liu, “Diffraction of plane P waves by a canyon of arbitrary shape in poroelastic halfspace (i): formulation,” Earthquake Science, vol. 22, no. 3, pp. 215–222, 2009. View at: Publisher Site  Google Scholar
 C. Yi, P. Zhang, D. Johansson, and U. Nyberg, “Dynamic response of a circular lined tunnel with an imperfect interface subjected to cylindrical Pwaves,” Computers and Geotechnics, vol. 55, no. 1, pp. 165–171, 2014. View at: Publisher Site  Google Scholar
 X.Q. Fang and H.X. Jin, “Viscoelastic imperfect bonding effect on dynamic response of a noncircular lined tunnel subjected to P and SV waves,” Soil Dynamics and Earthquake Engineering, vol. 88, pp. 1–7, 2016. View at: Publisher Site  Google Scholar
 X.Q. Fang and H.X. Jin, “Dynamic response of a noncircular lined tunnel with viscoelastic imperfect interface in the saturated poroelastic medium,” Computers and Geotechnics, vol. 83, pp. 98–105, 2017. View at: Publisher Site  Google Scholar
 G. Ping, Z. Lv, Q. Yan, and H. E. Chuan, “Internal force solution for a deepburied shield tunnel under the action of a seismic wave,” Modern Tunnelling Technology, vol. 53, no. 1, pp. 45–51, 2016. View at: Google Scholar
 E. Debiasi, A. Gajo, and D. Zonta, “On the seismic response of shallowburied rectangular structures,” Tunnelling and Underground Space Technology, vol. 38, pp. 99–113, 2013. View at: Publisher Site  Google Scholar
 C. Yi, D. Johansson, and U. Nyberg, “Scattering of SHwaves by a shallow circular lined tunnel with an imperfect interface,” in Proceedings of 8th International Symposium on Ground Support in Mining and Underground Construction, E. Nordlund, T. H. Jones, and A. Eitzenberger, Eds., Luleå, Sweden, September 2016. View at: Google Scholar
 X.Q. Fang, T.F. Zhang, and H.Y. Li, “Elasticslip interface effect on dynamic response of a lined tunnel in a semiinfinite alluvial valley under sh waves,” Tunnelling and Underground Space Technology, vol. 74, pp. 96–106, 2018. View at: Publisher Site  Google Scholar
 F. C. P. D. Barros and J. E. Luco, “Diffraction of obliquely incident waves by a cylindrical cavity embedded in a layered viscoelastic halfspace,” Soil Dynamics and Earthquake Engineering, vol. 12, no. 3, pp. 159–171, 1993. View at: Publisher Site  Google Scholar
 Z. Liu, X. Ju, C. Wu, and J. Liang, “Scattering of plane P 1 waves and dynamic stress concentration by a lined tunnel in a fluidsaturated poroelastic halfspace,” Tunnelling and Underground Space Technology, vol. 67, pp. 71–84, 2017. View at: Publisher Site  Google Scholar
 C. Smerzini, J. Avilés, R. Paolucci, and F. J. SánchezSesma, “Effect of underground cavities on surface earthquake ground motion under SH wave propagation,” Earthquake Engineering and Structural Dynamics, vol. 38, no. 12, pp. 1441–1460, 2010. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2019 Xiaojie Zhou 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.