Research Article  Open Access
ThreeDimensional Nonlinear Seismic Response of Immersed Tunnel in Horizontally Layered Site under Obliquely Incident SV Waves
Abstract
A threedimensional (3D) detailed numerical model of an immersed tunnel in a horizontally layered site is established in this study. The 3D seismic response of the immersed tunnel in a horizontally layered site subjected to obliquely incident waves is analyzed based on the precise dynamic stiffness matrix of the soil layer and halfspace via combined viscousspring boundary and equivalent node stress methods. The nonlinear effects of external and internal site conditions on the whole model were determined by equivalent linearization algorithm and Mohr–Coulomb model, respectively. The proposed model was then applied to investigate the nonlinear seismic response of an immersed tunnel in the Haihe River subjected to seismic waves of oblique incidence. The dislocation (opening) of pipe joints in the immersed tunnel were analyzed to determine the response characteristics of the shear keys and overall displacement of the tunnel; the dynamic responses of the immersed tunnel subjected to obliquely incident seismic waves markedly differ from those of vertically incident seismic SV waves. The maximum stress value of shear keys and the maximum dislocation of the pipe joint appear as upon critical angle. The overall displacement of the tunnel increases as incident angle increases. Under severe earthquake conditions, both the pipe corners and midspan section of the roof and floor are likely to produce crack. These areas need careful consideration in the seismic design of immersed tunnel structures.
1. Introduction
Immersed tunnel structures have grown increasingly common in crosssea and crossriver traffic engineering projects. Many of these structures are located in highly seismic regions. The unique characteristics of the immersed tunnel make any structural damage caused by an earthquake particularly difficult to be repaired, so antiseismic design is particularly important. The complexity of the immersed tunnel system makes it difficult to test by shaking table. Numerical analysis methods are an important means to study the seismic performance of immersed tunnels. Immersed tunnels tend to be lengthy and are usually built on weak ground [1, 2], so it is necessary to consider the spatial variation effects of ground motion and the dynamic interaction between the soil and structure in such numerical analyses [3]. Numerical simulation of the whole dynamic analysis model must include both the immersed tunnel and surrounding rock [4].
Recent researchers have conducted many shaking table tests [5–7] and fine numerical simulations [8–10] for investigating the 3D seismic response of immersed tunnels. Hu et al. [11] investigated the mechanical properties of segmental joints in the immersed tunnel of the HongkongZhuhaiMacao Bridge in China; Yu et al. [6] considered the multipoint input of ground motions via experiment. Chen et al. [5] carried out a large immersed tunneljointsite soil shakingtable test to determine the seismic response of immersed tunnels and joints.
It is generally assumed that seismic waves in a deepfocus earthquake strike perpendicularly to the bedrock and transfer it to the structure. For shallowfocus earthquakes with small epicentral distance, the influence of seismic wave incident angle on the seismic response of surrounding structures is an important consideration [12]. Du [13] analyzed the oblique incidence of seismic waves by using the transmission viscoelastic boundary conditions. FCPD Barros and Luco [14] obtained the threedimensional harmonic response in the vicinity of an infinitely long, cylindrical cavity of circular cross section buried in a layered, viscoelastic halfspace. Xiao et al. [15] investigated the shear behavior of an immersion joint subjected to compressive shear loads and found that the capacity of the joint is smaller than the sum of the capacities of all shear keys.
Previous researchers have mainly focused on 2D cross sections. However, it is difficult to accurately describe the nonlinear response state by simply decomposing the immersed tunnel into inplane and outofplane components under the oblique incidence of seismic waves. Unfortunately, it is very challenging to secure an accurate 3D model due to its complexity. Huang et al. [16] proposed a method to realize the SV and SH waves with arbitrary incident angles in ABAQUS and indicated that the nonlinear seismic responses of long lined tunnels are highly affected by the incident angles of S waves. The literature mainly centers around the halfspace uniform field rather than the 3D horizontally layered nonlinear seismic response of the site. Actual construction sites are usually layered; the amplification effect of the layered site is significantly different from that of the equivalent uniform site due to the resonance amplification and filtering effect of the soil layer [17], and this is proved in our paper [18].
In the present study, we established a 3D refined model of horizontal layered site immersed tunnel. We used the dynamic stiffness matrix of the soil layer and halfspace combined with viscousspring boundary and equivalent node stress methods to conduct accurate simulations of 3D seismic responses in a horizontally layered immersed tunnel structure. We applied the equivalent linearization algorithm (selfprogrammed) to determine the nonlinear effects of the external site of the overall model. We also used the Mohr–Coulomb model to determine the nonlinear effects of soil inside the overall model. The nonlinear effects of concrete were determined by a plasticity damage model of concrete in ABAQUS. The nonlinear seismic response of an immersed tunnel in the Haihe River and the seismic behavior of shear keys and pipe joints under obliquely incident SV waves were then determined based on the proposed model.
2. Obliquely Incident Seismic Waves in ABAQUS
2.1. ViscousSpring Boundary
We used the 3D viscousspring boundary set provided by Huang et al. [19] to simulate the stiffness and radiation damping effects of semiinfinite strata in this paper. The viscousspring boundary reflects the propagation of an advancing wave by superimposing the plane wave and the farfield scattered wave. The characteristics of wave attenuation and multiangle incidence are reflected by two empirical parameters; the wave field of the infinite field at the artificial boundary is transformed to stress on artificial boundaries. When the finiteelement integral is stable, the viscousspring boundary remains stable. The spring stiffness and damping coefficient of the 3D viscousspring artificial boundary are shown in equations (1) and (2):where K_{N} and C_{N} are the normal spring stiffness and damping coefficient of the nodes on the artificial boundary, respectively; K_{T} and C_{T} are the tangential spring stiffness and damping coefficient of the nodes on the artificial boundary, respectively; A and B are empirical parameters typically set to 0.8 and 1.1, respectively; G is the shear modulus; is the first Lame constant; r is the distance from the scattered wave source to the artificial boundary, which can be simply taken as the vertical distance from the geometric center of the structure to the artificial boundary; ρ is the soil density of the foundation; and and are S wave velocity and P wave velocity, respectively.
2.2. Earthquake Input Method Based on ViscousSpring Boundary
Here, we combine the finiteelement method with the timedomain method of viscousspring boundaries. The displacement and stress of the artificial boundary caused by equivalent load input on the artificial boundary must be consistent with that of the original free field. In this way, the earthquake input is translated into a freefield motion problem on artificial boundary nodes, whereas freefield motion can be transformed into an equivalent nodal load on artificial boundary nodes [20]. The equivalent nodal load of artificial boundary node l in direction i is calculated as follows:where is the influence area of node l on the artificial boundary surface; is the freefield stress tensor of node l in direction i; is the freefield displacement tensor of node l in direction i; is the freefield velocity vector of node l in direction i; is the viscousspring boundary spring stiffness of node l in direction i; and is the viscousspring boundary damping coefficient of node l in direction i.
2.3. Plane Wave Input under Oblique Incidence in a 3D Layered Site
The schematic diagram of the input wave field when the SV wave is obliquely incident is shown in Figure 1 [20]. The plane SV wave with the incident angle is reflected on the surface of the bedrock, forming a plane SV wave with the reflection angle and a plane P wave with the reflection angle , where and the critical value of the incident angle of the SV wave . When the problem is solved, the input schedule is first transformed to the frequency domain by Fourier transform, and then all the effective frequency points are calculated in the frequency domain. Finally, the timedomain results are obtained by inverse Fourier transformation. The direct stiffness method can be applied to the freefield solution in the frequency domain.
Wolf [21] gives an accurate dynamic stiffness matrix for layered soil and halfspace. The dynamic equilibrium equations of i layer and halfspace in the plane motion (PSV wave) are calculated by equations (4) and (5), respectively:where P_{i}, R_{i}, P_{b}, and R_{b} are horizontal and vertical load amplitudes of the upper and lower soil layers, respectively; P_{o} and R_{o} are horizontal and vertical load amplitudes of the free surface in halfspace; u_{i}, , u_{b}, and are horizontal and vertical displacements of the upper and lower surface of the soil layer; and u_{o} and are horizontal and vertical displacements of the bedrock surface. The stiffness matrix is a symmetric matrix. The expressions of k_{11}–k_{44} and r_{11}–r_{22} are available in our references [21].
The stiffness matrix of each soil layer can be set and adjusted to obtain the overall space stiffness matrix S_{PSV}, which is a (2n + 2) order square matrix where n is the number of soil layers. The overall dynamic equilibrium equation of the site can be written as follows:where P is the load vector, which is only related to the field characteristics, the incident wave amplitude, and the incident angle; the displacement vector u at the layer stratification can be solved accordingly.
The horizontal displacement and vertical displacement at any point within each soil layer can be expressed as equations (7a) and (7b). The velocity at any point can be derived from the displacement to the time as expressed by equations (8a) and (8b). Normal stress amplitude and as well as shear stress amplitude , can be derived from the following constitutive equation:where , , , and represent the incident and reflected amplitude of P and SV waves, respectively; is the wave number; is the phase wave velocity; is the circular frequency; and . For the incident SV wave, the direction cosine ; is the P wave propagation direction cosine.
In this study, we investigated the nonlinear seismic response of the site by the equivalent linearization method. The soil characteristics in the matrix are replaced by the equivalent shear modulus and damping ratio which match the equivalent shear strain:where is the maximum shear modulus; is a reference shear strain during the earthquake (the value can be determined experimentally or empirically); is the dynamic shear strain; and is the maximum damping ratio (which also can be obtained empirically or experimentally).
The coordinate system in the direct stiffness method differs from that in ABAQUS. After coordinate transformation, the solution is plugged into equation (3) to obtain the equivalent nodal load for use in ABAQUS.
3. Input Verification
We established a 3D horizontally layered finiteelement model (Figure 2(a)) and selected a 240 m × 240 m × 120 m cubic computational domain; 5 m × 5 m × 5 m hexagonal solid elements satisfying the precision requirements of the wave finiteelement simulation (equation (11)) were used for discretization, in which A and B are monitoring points. We divided the soil into three layers each 40 m in thickness with the parameters shown in Table 1. The incident plane wave is a pulse wave (equation (12)). The waveform is shown in Figure 3. The SV wave incident angle is . The angle between the incident wave projecting to the ZX plane and the X axis is , and the angle between the incident wave and the Yaxis is (i.e., incident angle in Section 2.3). We used the incident angle of (60°, 30°) for input verification.where is the largest size for the grid, is the minimum shear wave velocity, and is the maximum frequency of fluctuations.where f = 2, 0 ≤ t ≤ 0.5.
(a)
(b)

(a)
(b)
The reflection and refraction waveforms of the SV wave are shown in Figure 4. The viscousspring boundary absorbs the reflected wave very well. The close agreement between the numerical solution and the analytic solution proves that the obliqueincidence plane wave input method yields accurate simulation results in a 3D layered field with viscousspring boundaries.
(a)
(b)
(c)
(d)
4. Seismic Response of Immersed Tunnel under Seismic Wave Oblique Incidence
4.1. Calculation Model and Parameters
The tunnel immersed below the Tianjin Haihe River is a twohole, threepipe gallery structure (Figure 5(a)). The width of the immersed tunnel is 36.6 m, the height of its cross section is 9.65 m, and the thicknesses of the roof, floor, side wall, and partition wall are 1.35 m, 1.4 m, 1.0 m, and 0.6 m, respectively. The span of pipe gallery on both sides is 2.6 m, the span of the middle elements gallery is 2 m, and the span of the dual holes is 12.5 m. The tunnel consists of three pipes (E1–E3); the length of a single pipe is 85 m. The thickness of the crushedstone soil cover above the structure is 1.0 m. The tunnel is composed of C40 concrete. The clay layer in the Haihe Tunnel is defined here as the soil layer, which is divided into three sublayers. The material parameters of the soil and structure are shown in Tables 2 and 3.
(a)
(b)


We established a 3D finiteelement model (300 m × 255 m × 55 m) of the Haihe River immersed tunnel using the parameters discussed above. The overall model includes three pipes (E1–E3), soil layers, and superimposed soil (Figure 5(b)). Around the cutoff boundary and the bottom boundary, a viscousspring artificial boundary approximates an infinite field. We used the Mohr–Coulomb model for the soil and the plastic damage model in ABAQUS for the concrete.
4.2. Seismic Wave Input
A Tianjin wave was input at the bottom of the model with a 20 s duration, 0.02 s interval, and peak acceleration of 3.1 m/s^{2} (Figure 6). Per the incident angle definition above (Figure 2) and the critical incident angle of the SV wave in this model (32.3°), the propagation directions of seismic waves are selected as = 0°, 30°, 60° and = 0°, 15°, 30°; this produced nine angleofincidence combinations. We used the dynamic relaxation method for static application. The seismic wave was input after 1 s of the static calculation; the entire calculation time was 21 s.
4.3. Results
The pipe joint and shear key are significant components in the immersed tunnel—their inherent weakness may result in serious, even irreversible damage to the whole tunnel in the case of an earthquake. We first calculated and analyzed the mechanical properties of the pipe joint and the shear key, then the deformation of the tunnel (including three pipes) under typical earthquake conditions.
4.3.1. Analysis of Shear Keys
The three shear keys on the left side of the E1E2 section interface are numbered A, B, and C. Each shear key was marked with four corners (numbered 1–4) as monitoring points, as shown in Figure 7. The maximum principal stress, minimum principal stress, and shear stress of each monitoring point under the nine different incident angle combinations were extracted as shown in Tables 4–6. Among them, the maximum principal stress, minimum principal stress, and shear stress value of each monitoring point of shear key are maximal when the incidence angle is (30°, 30°), indicating that the angle combination of (30°, 30°) is the most unfavorable (which is close to the critical angle of incidence of SV waves in the halfspace). The secondworst is the angle of incidence (60°, 30°). Figure 8 shows a stress cloud chart of the horizontal shear key A and vertical shear key B (30°, 30°).
(a)
(b)
(c)


 
Note. The angle is expressed as (, ), the angle between the incident wave projecting to the ZX plane and the X axis is , and the angle between the incident wave and the Yaxis is . 
(a)
(b)
(c)
(d)
(e)
(f)
Tables 4–6 and Figure 8 altogether indicate the following.(1)The shear key produces significant compressive stress during an earthquake while the shear stress and tensile stress are relatively low.(2)When the seismic wave incidence transforms from vertical to oblique along the pipe cross section (that is, gradually increases), and the incident angle is small (less than 15°), the stress value of the shear keymonitoring point is fairly stable. With further increase in the incident angle, the stress value increases significantly to a maximum stress of 25.53 MPa, which is 7.1 times greater than the corresponding vertical incidence.(3)When the incidence direction of the seismic wave transforms from the tunnel transverse direction to the axial (that is, gradually increases), the maximum principal stress, minimum principal stress, and shear stress of the horizontal shear key A and vertical shear key B increase at first and then decrease. The maximum principal stress of the vertical shear key C increases at first and then decreases while its minimum principal stress and shear force continue to increase.(4)The horizontal shear stress distribution of compressive stress and shear stress is very uniform, the top corner has stress concentration phenomenon, and the tensile stress is mainly concentrated at the root of the shear key. The stress of the vertical shear keys B and C assumes a trapezoidal distribution, that is, the stress at the lower part of the shear key is greater than the upper stress. Tables 4–6 show that the stress of shear key B on the side wall is greater than that of shear key C on the middle partition wall. Seismic designers, to this effect, are recommended to focus on the root of the shear key which connects the shear key and immersed pipe body as well as the vertical shear key; the choice of reinforced concrete materials for the vertical shear key in the outer partition wall is particularly important.
4.3.2. Dislocation (Opening) of Pipe Joint
As shown in Figure 9, we selected five nodes of the roof and the floor, respectively, at the E1–E2 pipe joint interface and determined the relationship between the dislocation (opening) and the incident angle of the seismic wave. The results are shown in Figure 10.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
The vertical dislocation value can be obtained without indicating the direction of dislocation, as shown in Figure 10. In fact, the vertical dislocation directions at both ends of the cross section of the model are opposite.
The dislocation volume at the joint is significantly affected by the seismic wave incidence angle. The horizontal and vertical dislocations of the roof and floor of the immersed tunnel joints increase as the local seismic wave incidence changes from vertical to oblique along the pipe cross section ( increases) and that of the maximum value is under (30°, 30°). The maximum value of horizontal dislocation at the roof is 10.64 mm, which is 11.9 times greater than the corresponding vertical incidence. The maximum horizontal dislocation at the floor is 2.45 mm, which is 3.7 times greater than the corresponding vertical incidence. The maximum vertical dislocation at the roof is 4.44 mm, which is 6.5 times greater than the corresponding vertical incidence. The maximum vertical dislocation of the floor is 4.47 mm, which is 6.2 times greater than the corresponding vertical incidence.
When the incidence direction of the seismic wave transforms from the tunnel transverse direction to the axial ( increases), the horizontal dislocation at the position of the roof and the floor of the pipe joint interface first increases and then decreases as vertical dislocation continuously increases. The lengthwise opening at the tunnel joint is very small (0.2 mm maximum at the bottom of the joint).
The horizontal dislocations at joints are similar at a certain angle of incidence. The end nodes are larger than the intermediate nodes in the vertical dislocation, and the direction of dislocation is opposite; in other words, the immersed tunnel shows torsional deformation. When the incident angle is (30°, 30°), the Dvalue of the vertical dislocation between the end node and the intermediate node of the floor joint is maximal (0.014°), which is in the scope of security.
4.3.3. Overall Deformation of Immersed Tunnel
The nodes from the immersed tunnel’s central axis including three pipes were used to draw horizontal and vertical dislocation curves at the maximum stress (Figure 11).
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
Figure 11 shows that as the seismic wave incidence changes from vertical to oblique along the pipe cross section (that is, increases), the horizontal and vertical displacement of the immersed tunnel increases. When the incidence direction of seismic wave transforms from the tunnel transverse direction to the axial (that is, increases), the overall horizontal and vertical displacement of the tunnel first increases and then decreases.
Furthermore, when the seismic wave is transversely incident (), the horizontal and vertical displacement distributions at intermediate tunnel joints are larger than at both sides; the displacement of both sides of tunnel joints at the pipe ends is basically the same. The horizontal displacement difference and vertical displacement difference of both sides of the tunnel (E1 and E3) joint increase as α_{1} increases, resulting in axial tensile deformation. When the angle of incidence is (60°, 30°), the horizontal dislocation difference of the end tunnel section reaches its maximum value. The axial tensile strain is approximately 1.73e − 7, which is in the scope of security.
4.3.4. Distribution of Cracks in Pipe Joint Cross Section
The distribution of tensile damage factors of the cross section of the immersed tunnel joint under incident angle (30°, 30°) is shown in Figure 12. This is an approximate representation of the distribution and progression of cross section cracks.
(a)
(b)
(c)
(d)
(e)
As shown in Figure 12, the end of the side wall produces cracks first and is followed by the midspan of the floor, which is connected to horizontal shear key cracks; the midspan of the roof then cracks followed by the midspan of the roof in the middle pipe gallery and four corner points of the tunnel roof. Key parts that should be strengthened in an actual engineering scenario include the tunnel corner, the midspan of the roof, and the floor—that is to say, the connections between the horizontal shear key and floor.
5. Conclusion
The 3D seismic response of an underwater immersed tunnel under obliquely incident seismic waves was investigated in this study based on simulations of the viscousspring artificial boundary and seismic load. The influence of seismic SV wave incident angle on the 3D immersed tunnel joint and overall deformation of the tunnel is discussed by modelling the Tianjin Haihe River tunnel structure. Our main conclusions can be summarized as follows.(1)The oblique incidence of seismic SV waves, unlike the vertical incidence, has a marked effect on the seismic response of the immersed tunnel.(2)Among nine seismic wave combinations, the stress of the shear key is the largest under 30°, 30°; this is the worstcase angle to the site discussed in this paper. The tensile stress of the horizontal shear key is mainly concentrated where the shear key is connected to the main body of the immersed tube. The stress in the lower part of the vertical shear key is larger than that in the upper part, which makes it an important design point.(3)The horizontal and vertical dislocation of the roof and the floor of the immersed tunnel joints increase as the seismic wave incidence angle increases; they are maximal at an incident angle of 30°, 30°. The maximum relative dislocation of tunnel joints observed in this study was about 10 mm, while the average residual compression of GINA water stops is 15 mm–25 mm. In other words, the dislocation is within a safe range. The vertical dislocation of the end section is greater than that of the intermediate node and the dislocation directions are opposite, which suggests that the deformation of the immersed tunnel is torsional.(4)The overall horizontal and vertical displacement of the immersed tunnel increases as the seismic wave incidence angle increases. The displacement difference between the two ends of the same tunnel segment also increases, which can cause axial tensile deformation of the tunnel joints.(5)Under the oblique incidence of SV waves, cracks appear in the ends of the side walls (the corners of the tunnel), the midspan of the roof and the floor, and the part of the horizontal shear key connected to the floor. In an actual engineering scenario, construction measures and anticracking treatment should be targeted to these key areas.
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 there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
The research described in this paper was financially supported by the National Natural Science Foundation of China under grant numbers 51678389 and 51678390.
References
 Y. Yong and H. Yu, “Multishakingtable test of a long immersed tunnel,” in Proceedings of the ITA WTC 2015 Congress, Dubrovnik, Croatia, May 2015. View at: Google Scholar
 H. Yu, Y. Yuan, and A. Bobet, “Seismic analysis of long tunnels: a review of simplified and unified methods,” Underground Space, vol. 2, no. 2, 2017. View at: Publisher Site  Google Scholar
 K. Pitilakis, G. Tsinidis, A. Leanza, and M. Maugeri, “Seismic behaviour of circular tunnels accounting for above ground structures interaction effects,” Soil Dynamics and Earthquake Engineering, vol. 67, no. 67, pp. 1–15, 2014. View at: Publisher Site  Google Scholar
 I. Anastasopoulos, N. Gerolymos, V. Drosos, R. Kourkoulis, T. Georgarakos, and G. Gazetas, “Nonlinear response of deep immersed tunnel to strong seismic shaking,” Journal of Geotechnical and Geoenvironmental Engineering, vol. 133, no. 9, pp. 1067–1090, 2007. View at: Publisher Site  Google Scholar
 H. Chen, X. Li, W. Yan, S. Chen, and X. Zhang, “Shaking table test of immersed tunnel considering the geological condition,” Engineering Geology, vol. 227, pp. 93–107, 2017. View at: Publisher Site  Google Scholar
 H. Yu, Y. Yuan, G. Xu, Q. Su, X. Yan, and C. Li, “Multipoint shaking table test for long tunnels subjected to nonuniform seismic loadingspart II: application to the HZM immersed tunnel,” Soil Dynamics and Earthquake Engineering, vol. 108, pp. 187–195, 2018. View at: Google Scholar
 X. Cheng, L. Jing, J. Cui, Y. Li, and R. Dong, “Shakingtable tests for immersed tunnels at different sites,” Shock and Vibration, vol. 2017, no. 6, pp. 1–11, 2017. View at: Publisher Site  Google Scholar
 Z.n. Hu and Y.l. Xie, “Mechanical and failure characteristics of shear keys on immersed tunnel segment joints under differential settlements,” Procedia Engineering, vol. 166, pp. 373–378, 2016. View at: Publisher Site  Google Scholar
 J.H. Ding, X.L. Jin, Y.Z. Guo, and G.G. Li, “Numerical simulation for largescale seismic response analysis of immersed tunnel,” Engineering Structures, vol. 28, no. 10, pp. 1367–1377, 2006. View at: Publisher Site  Google Scholar
 X. L. Jin, Y. Z. Guo, and J. H. Ding, “Three dimensional numerical simulation of immersed tunnel seismic response based on elasticplastic FEM,” in Proceedings of the AsiaPacific Symposiumon Engineering PLASTICITY and ITS Applications, pp. 661–666, Australia, 2004. View at: Google Scholar
 Z.n. Hu, Y.l. Xie, G.p. Xu et al., “Segmental joint model tests of immersed tunnel on a settlement platform: a case study of the HongkongZhuhaiMacao bridge,” Tunnelling and Underground Space Technology, vol. 78, pp. 188–200, 2018. View at: Publisher Site  Google Scholar
 J. Yang, “Saturation effects on horizontal and vertical motions in a layered soilbedrock system due to inclined SV waves,” Soil Dynamics and Earthquake Engineering, vol. 21, no. 6, pp. 527–536, 2001. View at: Publisher Site  Google Scholar
 X. Du, “A stress artificial boundary in fea for nearfield wave problem,” Chinese Journal of Theoretical & Applied Mechanics, vol. 38, no. 1, pp. 49–56, 2006. View at: 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 & Earthquake Engineering, vol. 12, no. 3, pp. 159–171, 1993. View at: Publisher Site  Google Scholar
 W. Xiao, H. Yu, Y. Yuan, L. Taerwe, and G. Xu, “Compressionshear behavior of a scaled immersion joint with steel shear keys,” Tunnelling and Underground Space Technology, vol. 70, pp. 76–88, 2017. View at: Publisher Site  Google Scholar
 J. Q. Huang, X. L. Du, M. Zhao, and X. Zhao, “Impact of incident angles of earthquake shear (s) waves on 3d nonlinear seismic responses of long lined tunnels,” Engineering Geology, vol. 222, pp. 168–185, 2017. View at: Publisher Site  Google Scholar
 G. Wei and H. J. Qiu, “Summary of model tests and settlement characteristics of base layer in immersed tube tunnel,” Applied Mechanics and Materials, vol. 368–370, no. 1, pp. 1421–1425, 2013. View at: Publisher Site  Google Scholar
 Z. XiaoJie, Z. YueYu, H. E. Ying, and L. ZhongXian, “Seismic response analysis of immersed tunnel under oblique incidence of seismic SV wave,” China Earthquake Engineering Journal, vol. 39, no. 4, pp. 600–608, 2017. View at: Google Scholar
 J. Q. Huang, D. U. XiuLi, Z. Mi, and J. Liu, “A simplified seismic wave input method for nearfield wave analysis,” Zhendong Yu Chongji/Journal of Vibration & Shock, vol. 34, no. 3, pp. 77–82, 2015. View at: Google Scholar
 H. E. Jiantao, M. A. Huaifa, Z. Boyan et al., “Viscoelastic artificial boundary method of earthquake input and implementation,” Journal of Hydraulic Engineering, vol. 41, no. 8, pp. 960–969, 2010, in Chinese. View at: Google Scholar
 J. P. Wolf, Dynamic SoilStructure Interaction, PrenticeHall, Englewood Cliffs, NJ, USA, 1985.
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.