SeismicResponse Studyof TunnelsRunningunderneath a Subway Station in Parallel

A tunnel passing below a metro station is taken as the object of our study, and a two-dimensional plane model is established to study the effects of their dynamic interaction to seismic excitation. Comparative analysis is used to obtain the influence law between the underground structures and on the soil. .e results show that (1) the influence of the underground structure on the soil response is related to structure depth..e range of influence of the station structure on the surface is approximately five times the width of the station, and the surface response is obviously significant within this range. (2) .e existence of the tunnel is conducive to reducing the acceleration of the column in the station. It increases the displacement difference between the stations, but the impact on the internal force of the station is not significant. (3).e influence of the station on the dynamic response of the tunnel is consistent with regard to tunnel acceleration, difference in displacement between the top and bottom of the tunnel, and internal force. .e presence of the station will reduce the dynamic response of the tunnel.


Introduction
With increased development and utilization of urban underground space and the expansion and transformation of existing underground structures, underground engineering structures inevitably appear to pass each other at close distances, and cross-shaped and X-shaped traverse have emerged. As underground structures play an important role in the daily operation of the city, it is very critical to study the seismic performance of closed underground tunneling structures. Because an underground structure is buried in the soil, it is bound by the surrounding soil during an earthquake. e earthquake resistance mechanism of an underground structure is inevitably different from that of a surface structure: the seismic response of the surface structure is mainly controlled by the inertial force, whereas that of the underground structure is affected by the surrounding soil and mainly controlled by soil deformation. It is generally believed that the seismic performance of underground structures is superior to that of surface structures; therefore, the seismic problems of underground structures are often ignored. However, numerous underground structures have been seriously damaged in a series of large earthquakes around the world in the recent years. For example, the destruction of subway stations and tunnels caused by the Hanshin earthquake and Wenchuan earthquake indicates that underground structures are not absolutely safe, and their seismic performance needs to be studied [1][2][3].
Presently, there are many research studies on seismic resistance of underground structures, which have yielded constructive results. In terms of theoretical research, Lee and Trifunac [4] studied the dynamic response of tunnels during the early stage of an earthquake. e analytical solution of the tunnel dynamic response during an earthquake was derived through theoretical analysis and reasonable assumptions. In addition, the effects of soil-structure stiffness and tunnel depth on the dynamic response of the tunnel were analyzed. Bourouaiah et al. [5] revealed that soil properties have a great impact on soil-structure interactions, and the behavior at the interface between the soil and the structure depends directly on the soil properties, but the horizontal and vertical directions are affected differently. In terms of numerical simulation, Anastasopoulos et al. [6] simulated the immersed tunnel segment joints with nonlinear superelastic elements and studied the seismic response law of multisegment long-span tunnels under nonuniform seismic excitation. e analysis considered that the tunnel joint can effectively "absorb" a part of the deformation, and increasing the thickness of the joint gasket or reducing the length of the joint makes the structure more resistant to earthquakes. Pitilakis and Tsinidis [7] and Pitilakis et al. [8] analyzed the deformation behavior of rectangular underground structures during nonuniform earthquakes and pointed out that the special characteristics of underground structures make their seismic performance different from that of aboveground structures. Maugeri et al. [9] studied the performance of the shallow tunnel lining under seismic loading through numerical simulation and metaheuristic techniques. e results showed that genetic algorithm can be used to evaluate the safety performance of the tunnel lining. Wen et al. [10] studied the plastic failure behavior of tunnels, and the results showed that the plastic failure of the upper cross section of the tunnel was more serious than that of the lower part. Presently, there are numerous research results on the lateral response of structures, and the longitudinal response of the structure is gradually receiving attention. Park et al. [11] performed pseudostatic three-dimensional finite element analysis of tunnel structures under the action of spatially varied ground motion. It was pointed out that the spatial ground motion will cause the tunnel to bend longitudinally and increase the longitudinal axial force of the tunnel. Chen et al. [12] carried out a model shaking table test of the tunnel structure. It was considered that the dynamic response of the tunnel structure under nonuniform seismic excitation was greater than that under uniform seismic excitation. erefore, the spatial effect of seismic excitation should be considered in the tunnel design. A related test by Cilingir and Gopal Madabhushi [13,14] found that the ground motion response of the tunnel structure depends on the peak acceleration of the seismic excitation, and the change in the buried depth does not affect the tunnel deformation mode; it mainly affects the tunnel acceleration and internal force. Moghadam and Baziar [15] studied the mode of influence of the tunnel on the ground dynamic response during an earthquake by using a combination of the shaking table test and numerical simulation. Hashash et al. [16] reviewed the current state of seismic analysis and design of underground structures, introduced the quantitative methods used by engineers to assess the seismic action on underground structures, and expounded the deterministic and probabilistic methods for seismic hazard analysis. Numerical analysis tools such as the finite element method or finite difference method are used to analyze the dynamic soil-structure interactions. e effects of the relative structural flexibility of soil on the structural deformation and other properties have been studied, and calculation methods have been proposed. Abate and Massimino [17] studied the effects of the tunnel depth, aboveground structure position, and input seismic waveform on the dynamic interaction of the tunnel-soil-ground structure system based on an actual case in Italy. e system was fully coupled and analyzed.
irty different recorded accelerations were analyzed based on the time domain, frequency domain, seismic bending moment, and axial force on the tunnel lining.
Most of the existing related researches are directed to a single simple structure, such as a rectangular subway station or circle tunnel, and the influencing factors (seismic excitation amplitude and type, structural depth, structure size, and form) are relatively simple. ere is a lack of research on multiple structural systems, such as tunnels crossing each other, subway stations crossing each other, and tunnels crossing subway stations from below or above. In addition, their interactions have not been sufficiently explored, and it is difficult to draw a systematic general law. e purpose of this paper is to establish a two-dimensional model to study the influence of different positional relationships between structures on the dynamic response of the structures in a traversing system and draw some general conclusions.

Calculation Models
ere are many engineering cases of tunnels passing through subway stations in China. In order to reflect the generality, the research object of this article is not a specific engineering case, but an abstract general model from actual projects, which can obtain a general research result.

Calculation Range and Boundary Conditions.
is paper studies the seismic response of tunnels running underneath a metro station structure. Because the tunnel and station structure are not changed along the longitudinal section, the problem can be simplified to a plane problem under the action of seismic waves perpendicular to the longitudinal direction of the tunnel and the metro station. erefore, this paper establishes a two-dimensional plane model, and the calculation software used is ABAQUS. e subway station is a two-story and two-span frame structure with a cross-sectional dimension of 24 m × 14 m, side wall width of 0.8 m, center pillar width of 0.6 m, top plate and bottom plates are both of 0.9 m, and a middle floor slab thickness of 0.4 m. e outer diameter of the tunnel is 6 m, inner diameter is 5.3 m, thickness of the lining is 0.35 m, and the distance between the centers of the two tunnels is considered as 0.5D, D, 2D, and 3D (D is tunnel diameter). e tunnel depth is 17.5 m or 26 m, corresponding to shallow buried tunnels and deep buried tunnels, respectively. Based on existing research results [18], the width of the calculation model is 7 times the width of the station (three times the width of each station on the left and right sides), that is, 168 m, and the depth is 50 m to the bedrock surface.
e established station and tunnel models are shown in Figure 1. e finite element model includes a free-field model, spatially varying ground motion soil-station model, soil-double hole tunnel model, and soil-tunnel-station model, as shown in Figure 2. e size requirement of the grid unit is L ≤ λ/10 � v/10f, where λ is the shortest wavelength, v is the equivalent shear wave velocity of the soil, and f is the maximum frequency within the range of the ground motion energy to be considered. According to the soil material parameters, the equivalent shear wave velocity of the soil is approximately 220 m/s, and the ground motion energy is mainly concentrated within 10 Hz. Considering f � 10 Hz, we obtain L ≤ 2.2 m. erefore, the size of the grid unit is 1 m.
Because the modal analysis is performed first, the vibration characteristics of the system are determined, and then, the transient analysis is performed. e boundary conditions of the modal analysis are as follows: the bottom of the model is set as the fixed boundary, the side is set to the vertical displacement constraint (i.e., the lateral slip boundary), and the top is free. In the transient analysis, in order to input the horizontal transverse seismic wave at the bottom of the model, the horizontal constraint of the bottom of the model is removed; the other aspects of the analysis are the same as in the modal analysis.

Constitutive Model.
In this paper, the equivalent linear model, called the Davidenkov model, is used to consider the nonlinear characteristics of soil. Chazelas et al. [19] realized the secondary development of this model based on ABA-QUS. e Davidenkov model can be described as where A and B are fitting constants; c r is the reference shear strain; c d is the instantaneous dynamic shear strain; G d and λ are the instantaneous dynamic shear modulus and damping ratio; and G max and λ max are the maximum dynamic shear modulus and maximum damping ratio. When A � 1, B � 0.5, the Davidenkov model degenerates into the common Hardin-Drnevich model. Based on the existing dynamic test data of soil, different types of soils are selected with different A and B values, and the initial damping ratio is 0.003. e soil body adopts the solid element, and the unit type is CPE4R [20]. Because the equivalent nonlinear Davidenkov model is employed, the damping characteristics are considered in the model. e tunnel structure and station structure adopt the beam unit of unit type B21 [20], and the tunnel and the station adopt the elastic model. e calculation parameters of each material are shown in Table 1. e soil and station and the soil and tunnel are in tie contact, that is, the relative displacement of the contact interface is not considered. e famous El-Centro wave is selected firstly in the paper, and the amplitude of the input seismic wave was adjusted to 0.1 g (g is the gravitational acceleration) according to Chinese Code. en, in order to compare the seismic response of the tunnel-soil-station system under different seismic excitations, Kobe wave is selected for comparison, whose time history and spectrum curves are shown in Figure 3. e calculation uses the display algorithm; the step time of the calculation is taken as 0.002 s, and the calculation time is taken uniformly as 30 s.

Test Conditions.
El-Centro wave (near-field wave) ground motion is selected in the seismic response analysis. e acceleration time-history curve of the seismic wave and its Fourier spectrum curve are shown in Figure 3. e seismic response law of the structural system under the action of seismic waves is analyzed. e specific working conditions are shown in Table 2.
ere are 18 working conditions. In Tables 2, S represents the net spacing between the two tunnels, and H represents the net distance between the top of the tunnel and the bottom of the station.

Test Points Program of Models.
To monitor the seismic response of the soil and structure, the measuring points are arranged at equal intervals on the surface and in the depth direction of the soil, and the acceleration, velocity, and displacement response of the soil and structure are mainly monitored. e specific arrangement of the test points is shown in Figure 4. Due to the symmetrical structure, the surface test points are only on the semistructure and are arranged as one point every 4 m; the arrangement in the depth direction is one point every 2 m. e station measurement point is on the middle column, and the tunnel measurement points are at the top and bottom.

Verification of Numerical Analysis
In order to ensure the accuracy of the numerical simulation, this section verifies the tunnel shaking table test in reference [15].

Verification of Modal Analysis.
In order to verify the calculation model thoroughly, two models are both used in the paper. According to the actual size of the model box in reference [15], a 1 m × 1.6 m model is established (Figure 5(a)), and a prototype with a size of 160 m × 32 m is established according to the similarity relationship ( Figure 5(b)). e calculation model is established as discussed in Section 2, and the soil parameters are obtained according to the values in the reference. e first-order frequencies of the model and the prototype are 36.948 Hz and 2.748 Hz, respectively.
is shows a good agreement with the theoretical calculations of 37 Hz and 2.74 Hz; the relative errors are 0.14% and 0.29%, respectively, indicating the correctness of the modal analysis process.

Verification of Dynamic Analysis.
Because the input seismic excitation selected in the shaking table test in the literature cannot be obtained, the classical El-Centro wave is used to verify the relationship diagram shown in Figure 21 ( Figure 6(a) in this paper) in reference [15]. A free-field model (FF model) and a model containing an underground tunnel (SF model) are established, and the amplification factor at V S � 100 m/s is obtained, as shown in Figure 6(b).
It can be seen from Figure 6 that the simulation results in this paper are similar to those in reference [15]. When the dimensionless distance X/a � 0, the amplification factor is the smallest; when X/a � 3 (or 2), the amplification factor is the largest; and when X/a > 6, the amplification factor tends to 1. Despite the influence of different factors such as input seismic waves, the calculation results in this paper and the experimental data in [15] show only a slight difference in amplitude. ey maintain a good agreement on the regular trend, and the differences are within the acceptable error range.
erefore, the dynamic analysis model and the method followed in this paper can be considered reasonable and feasible and can be used for subsequent parameter analysis.

Modal Analysis.
e modal analysis of the established model is carried out to determine the vibration frequency characteristics of each working condition. e first 3 frequencies of the free field, soil-station model, soil-shallow buried tunnel model, and the corresponding soil-stationshallow buried tunnel model are shown in Table 3.
It can be seen from the Table 3 that (1) f � 1.1 Hz is calculated based on the theoretical formula f � V/4H for the self-vibration characteristics of the free-field soil. It is closer to the frequency of the model in the table, indicating the correctness of the model. (2) e fundamental frequencies of the working conditions are similar, and the underground structure has no significant influence on the dynamic characteristics of the system. Its influence is considered negligible from the engineering point of view.
Free-field FF19 Free-field Kobe wave (0.1 g) SS20 Soil-station Buried depth is 3  Figure 7 shows the acceleration-time history diagram of the midpoint of the surface of the free field under the action of 0.1 g El-Centro wave. e magnitudes of the acceleration at the midpoint of the ground are listed in Table 4 and compared with those of the free field. To facilitate quantitative analysis, the influence coefficient of the underground structure on soil is defined as where A 0 represents the positive acceleration amplitude of the center of the free-field table, and A CS represents the    Shock and Vibration 7 positive acceleration amplitude of the center of the surface when the underground structure exists. It can be seen from the Table 4 that the station structure has small buried depth and large cross-sectional size, which have a great influence on the surface acceleration. e station structure significantly increases the acceleration of the surface center, and the impact amplitude is approximately 25%. However, owing to the large depth and small size of the tunnel, the impact on the surface acceleration is small, which can be neglected.

Influence of Underground Structure on Soil Response.
Taking the deep buried tunnel as an example, the effects of the underground structure on the surface acceleration response and soil displacement along the depth of the soil are explored. e surface acceleration amplitude curves of the four systems ( Figure 8) and the soil displacement map in the depth direction (Figure 9) are presented.
It can be seen from the figure that the degree of influence of the station on the ground surface acceleration is greater than that of the tunnel; the influence of the station on the surface acceleration is approximately five times the station width (60 m on both sides of the center of the earth's surface). In this range, the surface acceleration is obviously larger than the free field, and the surface acceleration outside this range is different much to free field. Based on the displacement of the soil in the depth direction, the existence of the underground structure is beneficial in reducing the lateral displacement of the soil, and the influence of the station structure is greater than the influence of the tunnel structure. e influence of the station structure plays a major role when the station and the tunnel exist simultaneously.

Influence of the Seismic Wave Type.
e soil and structural responses may be different for different types of seismic waves. e soil response is taken as an example to  analyze the effects of different types of seismic waves and the surface acceleration of the free field, and the soil-station model ( Figure 10) and soil displacement map are obtained in the depth direction ( Figure 11) under different seismic waves. It can be seen from the figure that the influence law and degree of influence of different seismic waves on the soil response are different. Based on the surface acceleration response and soil displacement response, the influence of the El-Centro wave is larger than that of the Kobe wave.

Analysis of the Influence of Tunnel on Station Response.
is section compares the different aspects of the seismic response of the station with and without tunnels. To facilitate the analysis of the problem, the influence coefficient of the tunnel on the station is defined as where C C represents the response of the station in the earthstation system, and C CS represents the response of the station in the earth-tunnel-station system.

Acceleration Analysis of the Top of the Middle Column of the Station.
In order to explore the influence of the acceleration response at the top of the column in the station with the tunnel passing underneath, the acceleration response of the station center column is given for the case when only the station structure is present ( Figure 12). Table 5 lists the acceleration amplitudes for the station structures with and without the tunnel.
It can be seen from the Table 5 that the existence of the tunnel can reduce the acceleration of the column in the station by approximately 10%. Generally, the influence of the net spacing between the tunnels is not obvious. e greater the distance between the tunnel and the station, the smaller the acceleration of the column in the station.

Analysis of the Displacement Angle between Floors of the Station.
ere are three measuring points on the column in the station, which are located at the top, middle, and bottom of the column. e displacement response of these points can provide the difference in the interlayer displacement between the upper and lower layers and the  e time-history map of the displacement angle between the stations is shown in Figure 13 when only the station structure is present. e maximum displacement angle between the stations in other systems is given in Table 6.
It can be seen from the graph that the interlayer displacement angle of the lower layer of the station structure is large; the existence of the tunnel structure increases the interlayer displacement angle of the station structure. e smaller the spacing between the tunnel and the station and the smaller the spacing between the tunnels, the greater the increase in the displacement angle between the stations. When the distance between the tunnel and the station is large, the change in the net spacing between the tunnels will not cause a large change in the displacement angle between the stations.

Internal Force Analysis of the Station.
By observing the bending moment cloud diagram of the station structure, it is found that the maximum bending moment always appears at the corner of the station at 2.180 s. e bending moment cloud at 2.18 s is shown in Figure 14. e maximum values (absolute value) of the bending moment of the left corner of the station for other conditions are listed in Table 7.
It can be seen from the graph that the bending moment at the corner of the station structure is large and stress concentration can easily occur. When the distance between the tunnel and the station is small, the existence of the tunnel increases the bending moment of the station, but the overall impact is not significant.

Analysis of the Influence of the Station on Tunnel Response.
is section analyzes the response of the tunnel and discusses the impact of the station on the dynamic response of the tunnel. To facilitate the analysis of the problem, the influence coefficient of the station on the tunnel is defined as where S S represents the response of the tunnel in the soiltunnel system and S CS represents the response of the tunnel in the corresponding soil-tunnel-station system.

Analysis of Acceleration at the Top of the Tunnel.
Taking the soil-shallow buried tunnel with a tunnel spacing of 3 m as an example, the acceleration diagram of the top of the tunnel is shown in Figure 15. e maximum accelerations of the top of the tunnel under other conditions are shown in Table 8. It can be seen from the chart that stations influence shallow buried tunnels and deep buried tunnels differently. In the case of a shallow buried tunnel (distance between the shallow buried tunnel and station is 0.5 m), the presence of a station slightly enlarges the acceleration of the tunnel top. e amplification is between 2% and 4%;      the larger the net spacing between the tunnels, the more obvious the amplification effect. In the case of a deep buried tunnel, the presence of the station weakens the acceleration response at the top of the tunnel by 5%-9%; the degree of attenuation decreases as the spacing between the tunnels increases.

Analysis of the Displacement Difference between the
Top and Bottom of the Tunnel. Taking the soil-shallow buried tunnel with a tunnel spacing of 3 m as an example, the displacement difference between the top and bottom of the tunnel is shown in Figure 16. e maximum displacement differences between the top and bottom of the tunnel for various conditions are shown in Table 9. It can be seen from the chart that the station reduces the difference in the horizontal displacement between the top and bottom of the tunnel, and the degree of attenuation decreases with an increase in the spacing between the tunnels. However, the impact amplitudes of the shallow tunnel and deep tunnel are different. e attenuation for the shallow tunnel is 4%-15%, whereas that for the deep tunnel is 16%-27%.

Internal Force
Analysis of the Tunnel. By observing the bending moment cloud diagram of the tunnel structure, it is found that the maximum bending moment appears in the 45°d irection of the tunnel structure. e bending moment cloud diagram of the tunnel structure in the soil-shallow buried tunnel system (when the tunnel spacing is 3 m) is shown in Figure 17. e maximum value of the bending moment of the working tunnel structure is shown in Table 10.
It can be seen from the graph that the maximum bending moment of the tunnel structure appears in the 45°direction of the tunnel section. e larger the spacing between the tunnels, the larger the bending moment on the tunnel; the influence of the station structure on the tunnel bending moment is inconsistent. In general, the presence of a station can reduce the bending moment of the tunnel structure.

Conclusions
In this paper, a two-hole tunnel underpass station model was established. e ABAQUS finite element software was used to analyze the seismic response. e seismic response law of the cross-underground structure was discussed. e following rules were obtained from the calculation and analysis: (1) Influence of the underground structure on the dynamic response of soil: the station structure has a great impact whereas the tunnel structure has less impact on the surface response. e range within which the station structure has a great impact on the surface response is five times the width of the station. Within this range, the surface acceleration response increases significantly, whereas the influence outside this range is small. e tunnel has a large buried depth, and its influence on the surface response is basically negligible. e station depth is small because of which it has a great influence on the surface response. When the station and the tunnel structure exist simultaneously, the influence of the station structure plays a major role. Overall, the presence of underground structures reduces the horizontal displacement of the soil.
(2) Impact of the tunnel on the dynamic response of the station: the existence of the tunnel can reduce the acceleration of the column in the station by approximately 10%; the impact of the net spacing between the tunnels is not obvious. When the distance between the tunnel and the station is large, the acceleration of the column in the station is greatly weakened. e existence of the tunnel structure increases the interlayer displacement angle of the station structure; the smaller the distance between the tunnel and the station and the smaller the spacing between the tunnels, the greater the increase in the displacement between the stations. e impact of the tunnel on the internal force response of the station is not significant.
(3) Impact of the station on the dynamic response of the tunnel: the station influences the acceleration of the shallow buried tunnel and the deep buried tunnel differently. In the case of the shallow buried tunnel, the existence of the station slightly amplifies the acceleration response at the top of the tunnel. In the case of the deep buried tunnel, the presence of the station weakens the acceleration response at the top of the tunnel. As the spacing between the tunnels increases, the degree of attenuation decreases. e presence of the station reduces the difference in displacement between the top and bottom of the tunnel, and with an increase in the spacing between the tunnels, the degree of attenuation decreases. e presence of the station is conducive to reducing the internal force of the tunnel.
At present, the research on seismic resistance of complex underground structure is limited, and the related research still needs further improvement. Underground structures tend to be large in size and traverse different strata. Moreover, they do not exist independently and have a complex environment. erefore, it is 14 Shock and Vibration necessary to consider more influencing factors, such as structural characteristics, changes in soil parameters, seismic excitation characteristics, spatial effects of ground motion effects, and aboveground buildings.

Data Availability
No data were used to support the findings of the study.

Conflicts of Interest
e authors declare that they have no conflicts of interest.

Authors' Contributions
Fuxue Sun conducted investigation on existing research results extensively and deeply and supervised the development of research work. Wang Guo-bo conceptualized the actual project, determined the reasonable modeling method, and acquired funds to ensure the success of the research work. Xiang-jun Peng was responsible for modelling and simulating actual engineering problems through software and data analysis. Zhou-zhou Jin verified the correctness and reliability of the results and wrote the original draft. Xiaochun Li administrated the actual project and provided necessary resources for the research. Jie-lin Zhao was responsible for reviewing and editing the original draft to improve overall English language usage, readability, clarity, grammar, and spelling.