Response Displacement Method for Seismic Calculation of Subway Station Complex Structure of TOD Mode

: The TOD mode with rail transit stations as the center has become an important direction of future central city construction. However, at present, there is still a lack of simpliﬁed calculation methods for the seismic design of a subway station complex structure of TOD mode. Based on the load-structure model of the classic response displacement method and the seismic deformation characteristics of the underground complex structure, an improved response displacement method for seismic calculation of the subway station complex structure considering the inﬂuence of the upper frame structure was proposed in this paper. Then, based on the actual engineering case, the applicability of the improved reaction displacement method was veriﬁed through investigating the inﬂuencing factors such as seismic wave type, ground motion intensity, building position, structure form, structure stiffness, and stratum stiffness. Furthermore, a series of numerical simulation experiments were conducted to verify the simpliﬁed method and further evaluate its computational accuracy. The result shows that the error in calculating the internal force and deformation response of a station complex structure by using the improved reaction displacement method can be controlled at about 10%. The improved response deformation method is proved to be a highly practical pseudo-static method.


Introduction
In recent years, with the rapid development of urbanization in China, the urban development model of high-density and the land use mode of intensity have gradually become the main development direction of future urban construction [1].Under this background, the TOD mode appears based on the rail transit subway station as the core, integrating transportation, offices, a shopping mall, residences and other functions in the urban complex structure.The subway station complex structure of TOD mode was of a dual characteristic-both common underground subway station and aboveground frame structure.Currently, the research on its seismic response characteristics and failure mechanism is still insufficient, and there is a lack of simplified methods to directly guide its seismic design [2,3].
Wang G B [4] studied the influence of the upper frame structure on the seismic response of the underground subway station via numerical simulation.The analysis shows that the influence of the lighter upper structure on the dynamic response of the structural system is limited.Zhang et al. [5] established a three-dimensional finite element model of a structure system of a soil-subway station and its upper cover structure by using ABAQUS software, then the seismic performance of the structure system was analyzed, and the effects of different vertical ground motion and beam stiffness on the system were obtained.An [6] analyzed the seismic response of a large-chassis subway station under an upper cover single tower frame structure by using the finite element software Midas GTS NX, and the calculation results show that the existence of the upper cover structure increases the internal force of the subway station structure.Han [7,8] established a three-dimensional finite element numerical model of an underground subway station-soilground building integrated structural system, and analyzed the seismic response, seismic damage evolution process, and failure mechanism of the system under different influencing factors.The results show that the existence of surface buildings changes the internal force and dynamic characteristics of the original station structure.Liu [9] established the dynamic analysis model of the integrated structural system by using the SATWE module in the design software PKPM and the finite element software MIDAS.The elastic analysis under frequent earthquakes and the elastic-plastic analysis under rare earthquakes were analyzed, and the seismic performance of the structural system was evaluated.Currently, the most widely used method in the seismic design of underground structures is the response displacement method [10].The traditional response displacement method was improved by many scholars aiming at the limitation of it, such as the integral response displacement method [11], the inertial force-displacement method [12], the generalized response displacement method [13,14] and so on.Among them, Qiu [15] proposed the seismic design method of subway station structure in the vicinity of the existence of ground buildings, and the influence of the size of the numerical model on the accuracy of the simplified method was analyzed.The results showed that the influence of the size effect on the calculation results can be ignored when the distance between the model boundary and the structure exceeds four times the length of the structure.However, previous studies have focused more on the improvement of the traditional response displacement method with errors, and have not considered the seismic response characteristics of the underground complex structure under the upper structure, which is not suitable for direct use.
Therefore, this paper proposes an improved response displacement method suitable for the seismic analysis of the subway station complex structure of the TOD mode considering the influence of the additional seismic load generated by the upper frame structure on the subway station complex structure based on the load-structure model of the traditional response displacement method.The calculation results of time-history analysis are used as exact solutions to discuss the applicability of the improved response displacement method.

Proposition and Implementation of Improved Response Displacement Method 2.1. Traditional Response Displacement Method and Its Limitations
The response displacement method is a common method for seismic design of underground structures.It uses beam element modeling and soil springs to consider the dynamic response and interaction between underground structures and soil.The seismic load is applied to the structure in the form of relative displacement of the soil layer in the free field, the shear force of soil around the structure, and the peak acceleration of ground.Finally, the internal force and deformation of the structure under seismic load are obtained via calculation.The simplified mechanical model is shown in Figure 1.

Improved Response Displacement Method
Compared with common subway stations, the TOD structure model is an integrated structural system of ground and underground.The upper cover structure is directly em- The response displacement method has a rigorous theory and high calculation accuracy.The upper structure has a great influence on the internal force and deformation of the subway station [6].Applying the load-structure model of the traditional response displacement method to the seismic analysis of the subway station structure under the upper structure will lead to two problems.On the one hand, the unfavorable factors of the upper cover structure are not considered, resulting in a smaller calculation result.On the other hand, it cannot reflect the influence of changes in the internal force and deformation of the upper-cover structure of the subway station.Therefore, the calculation results of the traditional response displacement method do not meet the requirements for the integrated subway station and its upper cover structure.

Improved Response Displacement Method
Compared with common subway stations, the TOD structure model is an integrated structural system of ground and underground.The upper cover structure is directly embedded in the roof of the subway station complex.Under a horizontal earthquake, the base of the upper structure will generate a large bottom bending moment, torque, and shear force, etc.Meanwhile, as the upper structure is often a high-rise building, the vertical inertia force generated under the vertical earthquake is not negligible.Here, the additional internal force of the foundation base and vertical inertial force generated by the upper cover structure are unified as the additional seismic load of the upper cover structure.Based on the load-structure model of the traditional response displacement method and the seismic response law of the underground complex structure of the TOD mode, a simplified mechanical model for the seismic calculation of the subway station complex structure considering the influence of the upper cover structure frame structure is proposed, as shown in Figure 2.

Improved Response Displacement Method
Compared with common subway stations, the TOD structure model is an integrated structural system of ground and underground.The upper cover structure is directly embedded in the roof of the subway station complex.Under a horizontal earthquake, the base of the upper structure will generate a large bottom bending moment, torque, and shear force, etc.Meanwhile, as the upper structure is often a high-rise building, the vertical inertia force generated under the vertical earthquake is not negligible.Here, the additional internal force of the foundation base and vertical inertial force generated by the upper cover structure are unified as the additional seismic load of the upper cover structure.Based on the load-structure model of the traditional response displacement method and the seismic response law of the underground complex structure of the TOD mode, a simplified mechanical model for the seismic calculation of the subway station complex structure considering the influence of the upper cover structure frame structure is proposed, as shown in Figure 2.

Additional Seismic Load of Superstructure
Both the equivalent base shear method and the mode-superposition response spectrum method are based on the basic assumption that the structure is in linear elasticity.The traditional mode-superposition response spectrum method and the equivalent base shear method assume that the bottom is a fixed end in seismic calculation.The premise of

Additional Seismic Load of Superstructure
Both the equivalent base shear method and the mode-superposition response spectrum method are based on the basic assumption that the structure is in linear elasticity.The traditional mode-superposition response spectrum method and the equivalent base shear method assume that the bottom is a fixed end in seismic calculation.The premise of this assumption is that the foundation stiffness is very large.For the seismic design of the underground complex structure of the TOD mode, the upper structure is directly embedded in the roof of the lower large underground complex structure, and the stiffness of the underground structure is greater than that of the general foundation.Therefore, it is more in line with the assumption of considering setting the bottom of the upper cover structure as a fixed-end constraint to construct a calculation diagram.It is feasible to use the traditional 'gourd string' model to calculate the seismic action of the superstructure.
(1) Bottom shear The calculation of the bottom shear force of the upper cover structure can be obtained by the mode-superposition response spectrum method.The calculation formula is: where G i is the weight of the simplified particle i of the structural floor; α j is the jth vibration mode seismic influence coefficient calculated according to the jth order period of the system; γ j is the participation coefficient of the j-mode; φ ji is the mode displacement of j vibration mode at i particle.
When the height of the upper structure does not exceed 40 m and it is a frame structure, the bottom shear force of the building can also be obtained by the equivalent base shear method.The calculation formula is: where F EK is the bottom shear force of the structure; G eq is the structural equivalent total mass load; α 1 is the seismic influence coefficient of the first mode.
(2) Bottom bending moment Under the action of the horizontal earthquake, the bending moment generated by the bottom of the upper-structure column directly acts on the top plate of the subway station structure, which can be approximately solved by the anti-bending point method.It is assumed that the contra flexure point of the bottom column is at a height of 2/3 column from the top plate of the station.
The shear force of the bottom columns of the upper cover structure can be calculated using the following formula: where i k is the linear stiffness of the kth column in the bottom floor; ∑ m k =1 i k is the sum of the linear stiffnesses of all the columns in the bottom floor; V k is the shear force of the kth column in the bottom floor; and V is the sum of the shear forces of the columns in the bottom floor.
After obtaining the shear force of each column, the column bottom bending moment of each column can be obtained according to the position of the assumed contra flexure point.
(3) Vertical inertia force The equivalent base shear method can be used to calculate the vertical inertia force of the upper building.It is an equivalent horizontal seismic action method.The calculation diagram is shown in Figure 3.The calculation formula is: where F EVk is the standard value of the total vertical seismic action of the structure; F vi is the standard value of vertical seismic action of particle i; α max is the maximum horizontal seismic influence coefficient; G eg is the structural equivalent total gravity load; G i is the gravity load of the layer i; H i is the height of the particle from the ground.
where  is the standard value of the total vertical seismic action of the structure;  is the standard value of vertical seismic action of particle i;  is the maximum horizontal seismic influence coefficient;  is the structural equivalent total gravity load;  is the gravity load of the layer i;  is the height of the particle from the ground.The vertical seismic load is calculated according to the distribution ratio of the gravity load borne by each component multiplied by the dynamic response increasing coefficient of 1.5.

Improved Response Displacement Method Implementation Steps
(1) One-dimensional free-field analysis.Free Field Analysis, a program included in Midas GTS NX, is used to conduct the 1D free-field analysis, find the moment when the maximum relative displacement of the soil occurs at the top and bottom plate positions, and record the relative displacement of the soil layer from the top plate to the bottom plate, the soil shear force at the top and bottom plate positions, and the horizontal acceleration of the ground soil layer at this moment.
(2) Solve for the foundation spring stiffness.The convergent shear modulus of each soil layer calculated using the one-dimensional free field is used to establish a finite element model of the soil layer, as shown in Figure 4.The horizontal and vertical uniformly distributed loads are applied at the location of the model structure to obtain the deformation in both directions and to derive the coefficient of foundation, which can also be derived from the inverse by applying forced displacements.Finally, the foundation spring stiffness is obtained by the formula k = Kld, where k is the foundation spring stiffness; K is the bed coefficient; l is the concentrated spring spacing of the foundation; d is the calculated length of the stratum along the longitudinal direction of the underground structure.The vertical seismic load is calculated according to the distribution ratio of the gravity load borne by each component multiplied by the dynamic response increasing coefficient of 1.5.

Improved Response Displacement Method Implementation Steps
(1) One-dimensional free-field analysis.Free Field Analysis, a program included in Midas GTS NX, is used to conduct the 1D free-field analysis, find the moment when the maximum relative displacement of the soil occurs at the top and bottom plate positions, and record the relative displacement of the soil layer from the top plate to the bottom plate, the soil shear force at the top and bottom plate positions, and the horizontal acceleration of the ground soil layer at this moment.
(2) Solve for the foundation spring stiffness.The convergent shear modulus of each soil layer calculated using the one-dimensional free field is used to establish a finite element model of the soil layer, as shown in Figure 4.The horizontal and vertical uniformly distributed loads are applied at the location of the model structure to obtain the deformation in both directions and to derive the coefficient of foundation, which can also be derived from the inverse by applying forced displacements.Finally, the foundation spring stiffness is obtained by the formula k = Kld, where k is the foundation spring stiffness; K is the bed coefficient; l is the concentrated spring spacing of the foundation; d is the calculated length of the stratum along the longitudinal direction of the underground structure.

𝐹
0.65 ,  1, … , where  is the standard value of the total vertical seismic action of the structure;  is the standard value of vertical seismic action of particle i;  is the maximum horizontal seismic influence coefficient;  is the structural equivalent total gravity load;  is the gravity load of the layer i;  is the height of the particle from the ground.The vertical seismic load is calculated according to the distribution ratio of the gravity load borne by each component multiplied by the dynamic response increasing coefficient of 1.5.

Improved Response Displacement Method Implementation Steps
(1) One-dimensional free-field analysis.Free Field Analysis, a program included in Midas GTS NX, is used to conduct the 1D free-field analysis, find the moment when the maximum relative displacement of the soil occurs at the top and bottom plate positions, and record the relative displacement of the soil layer from the top plate to the bottom plate, the soil shear force at the top and bottom plate positions, and the horizontal acceleration of the ground soil layer at this moment.
(2) Solve for the foundation spring stiffness.The convergent shear modulus of each soil layer calculated using the one-dimensional free field is used to establish a finite element model of the soil layer, as shown in Figure 4.The horizontal and vertical uniformly distributed loads are applied at the location of the model structure to obtain the deformation in both directions and to derive the coefficient of foundation, which can also be derived from the inverse by applying forced displacements.Finally, the foundation spring stiffness is obtained by the formula k = Kld, where k is the foundation spring stiffness; K is the bed coefficient; l is the concentrated spring spacing of the foundation; d is the calculated length of the stratum along the longitudinal direction of the underground structure.(3) Additional seismic load on the upper structure.The calculation process refers to Section 2.2.1 to find out the key parameters such as bottom shear force, bottom bending moment, and vertical inertia force for the seismic load of the upper structure.
(4) Establishment of improved load-structure model of the response displacement method.Ground movement, ground shear force, horizontal and vertical acceleration of the structure, and additional seismic load of the upper cover structure are applied at the corresponding position of the structure: where the stratum displacement is imposed at the end of the horizontal soil spring on the side walls and bedplate of the structure; the stratigraphic shear force is applied on the side wall and the bottom plate of the structure; the horizontal acceleration applied to the structure is the horizontal relative acceleration of each soil layer along the height in one-dimensional free-field analysis; and the additional seismic load is applied at the location where the upper structure column is connected to the top slab of the subway station.

Calculation Model and Parameters
The calculation model comes from a project in Tongzhou, Beijing, where the underground subway station complex is a three-story and ten-span box frame structure, the above-tower building is a frame structure, and the cross section of the structural system and its corresponding dimensions are shown in Figure 5.The cross-sectional dimension of the column in the underground metro station is 0.8 m × 1.2 m, using C45 concrete; the side walls are 1.1 m thick, and the thickness of the top, middle, and bottom slabs are 1.0 m, 0.4 m, and 1.2 m, respectively, using C40 concrete; the slab of the upper structure is 0.15 m thick, and the cross-sectional dimension of the column is 0.7 m × 0.7 m, and the upper structure is all made of C35 concrete.The parameters of each soil layer are derived from the geological survey report of the project site.The soil layers with similar physical and mechanical properties are merged and simplified into eight soil layers.The main parameters of the soil layer are shown in Table 1.
(4) Establishment of improved load-structure model of the response displacement method.Ground movement, ground shear force, horizontal and vertical acceleration of the structure, and additional seismic load of the upper cover structure are applied at the corresponding position of the structure: where the stratum displacement is imposed at the end of the horizontal soil spring on the side walls and bedplate of the structure; the stratigraphic shear force is applied on the side wall and the bottom plate of the structure; the horizontal acceleration applied to the structure is the horizontal relative acceleration of each soil layer along the height in one-dimensional free-field analysis; and the additional seismic load is applied at the location where the upper structure column is connected to the top slab of the subway station.

Calculation Model and Parameters
The calculation model comes from a project in Tongzhou, Beijing, where the underground subway station complex is a three-story and ten-span box frame structure, the above-tower building is a frame structure, and the cross section of the structural system and its corresponding dimensions are shown in Figure 5.The cross-sectional dimension of the column in the underground metro station is 0.8 m × 1.2 m, using C45 concrete; the side walls are 1.1 m thick, and the thickness of the top, middle, and bottom slabs are 1.0 m, 0.4 m, and 1.2 m, respectively, using C40 concrete; the slab of the upper structure is 0.15 m thick, and the cross-sectional dimension of the column is 0.7 m × 0.7 m, and the upper structure is all made of C35 concrete.The parameters of each soil layer are derived from the geological survey report of the project site.The soil layers with similar physical and mechanical properties are merged and simplified into eight soil layers.The main parameters of the soil layer are shown in Table 1.Based on the project, the improved response displacement method model for seismic calculation of underground structures is shown in Figure 6.The load-structure model is modeled by the equivalent frame method [16], that is, the middle column of the subway station is taken as the prototype size of 0.8 m × 1.2 m, the lateral width of the side wall and the top and bottom plates is taken as the actual width, and the longitudinal length is taken as the sum of the half column spacing before and after the middle column, that is 8.0 m.At the same time, the calculation accuracy of the improved response displacement method was evaluated, and the calculation results of the time-history analysis method were used as the benchmark.The finite element software Midas GTS NX was used to establish a three-dimensional numerical model of soil-structure interaction with a model length of 560 m, a width of 80 m, and a soil thickness of 69 m, as shown in Figure 7.The modified Moore-Coulomb constitutive model and elasticity constitutive model are adopted for the model soil and structure, respectively, and the model system adopts Rayleigh damping with the damping ratio taken as 5%.The bottom boundary of the numerical model is a fixed constraint, and the free field boundary is used all around.The structural dimensions, material parameters, and seismic waves and their peak accelerations of the numerical model are kept consistent with the pseudo-static model.For the structural system under the TOD mode to carry out the internal force and deformation analysis under seismic action, the structural system is designed in accordance with the elastic behavior, and it is assumed that the structure and members are in the elastic working state, and the internal force and deformation analysis adopts the linear dynamic method.
station is taken as the prototype size of 0.8 m × 1.2 m, the lateral width of the side wall and the top and bottom plates is taken as the actual width, and the longitudinal length is taken as the sum of the half column spacing before and after the middle column, that is 8.0 m.At the same time, the calculation accuracy of the improved response displacement method was evaluated, and the calculation results of the time-history analysis method were used as the benchmark.The finite element software Midas GTS NX was used to establish a three-dimensional numerical model of soil-structure interaction with a model length of 560 m, a width of 80 m, and a soil thickness of 69 m, as shown in Figure 7.The modified Moore-Coulomb constitutive model and elasticity constitutive model are adopted for the model soil and structure, respectively, and the model system adopts Rayleigh damping with the damping ratio taken as 5%.The bottom boundary of the numerical model is a fixed constraint, and the free field boundary is used all around.The structural dimensions, material parameters, and seismic waves and their peak accelerations of the numerical model are kept consistent with the pseudo-static model.For the structural system under the TOD mode to carry out the internal force and deformation analysis under seismic action, the structural system is designed in accordance with the elastic behavior, and it is assumed that the structure and members are in the elastic working state, and the internal force and deformation analysis adopts the linear dynamic method.

Input Ground Motion and Seismic Response of Underground Structure
In order to analyze the applicability of the improved response displacement method under the different seismic waves, the Beijing artificial wave and the Kobe wave, the nearfield seismic wave Loma Prieta wave and the far-field seismic wave Landers wave were selected as the input seismic bedrock motion, and the acceleration-time-course curves of seismic waves are shown in Figure 8.The acceleration response spectra of the three seismic components of the seismic waves are shown in Figure 9.In addition, the analysis of the influencing factors of the improved response displacement method is based on the calculated results of the structural reaction under the Beijing artificial wave, and the conditions under other seismic waves are compared with this condition.

Input Ground Motion and Seismic Response of Underground Structure
In order to analyze the applicability of the improved response displacement method under the different seismic waves, the Beijing artificial wave and the Kobe wave, the nearfield seismic wave Loma Prieta wave and the far-field seismic wave Landers wave were selected as the input seismic bedrock motion, and the acceleration-time-course curves of seismic waves are shown in Figure 8.The acceleration response spectra of the three seismic components of the seismic waves are shown in Figure 9.In addition, the analysis of the influencing factors of the improved response displacement method is based on the calculated results of the structural reaction under the Beijing artificial wave, and the conditions under other seismic waves are compared with this condition.
field seismic wave Loma Prieta wave and the far-field seismic wave Landers wave were selected as the input seismic bedrock motion, and the acceleration-time-course curves of seismic waves are shown in Figure 8.The acceleration response spectra of the three seismic components of the seismic waves are shown in Figure 9.In addition, the analysis of the influencing factors of the improved response displacement method is based on the calculated results of the structural reaction under the Beijing artificial wave, and the conditions under other seismic waves are compared with this condition.
(a) (  For the convenience of analysis, the middle columns in the subway station are numbered column 1 to column 9, whereas the columns in the subway stations with upper cover structures are numbered column 3 to column 7, as shown in Figure 10.The calculated results of the bending moment and deformation of the middle column of the subway station complex and its superstructure system under the action of the Beijing artificial wave with a peak acceleration of 0.2 g are shown in Figure 10.As can be seen from the bending moment cloud atlas in Figure 10, the top of the negative first floor and the bottom of the negative third floor of each middle column of the subway station have significantly larger bending moments compared to the rest of the column.At the same time, the internal force and deformation of the middle column in the area with and without the upper cover structure are also different.Comparing the top bending moment of the middle column of the negative layer of the subway station at different locations, it is found that the bending moment of the building with upper cover is obviously larger than that of the building without upper cover, and the deformation of the underground structure with upper cover is still shearing type deformation, which is consistent with the conclusion obtained in Reference [6].In addition, the horizontal deformation of the middle column of the underground structure is slightly larger than that of the side wall, and that of the middle column at the position directly connected with the upper building is the largest.For the convenience of analysis, the middle columns in the subway station are numbered column 1 to column 9, whereas the columns in the subway stations with upper cover structures are numbered column 3 to column 7, as shown in Figure 10.The calculated results of the bending moment and deformation of the middle column of the subway station complex and its superstructure system under the action of the Beijing artificial wave with a peak acceleration of 0.2 g are shown in Figure 10.As can be seen from the bending moment cloud atlas in Figure 10, the top of the negative first floor and the bottom of the negative third floor of each middle column of the subway station have significantly larger bending moments compared to the rest of the column.At the same time, the internal force and deformation of the middle column in the area with and without the upper cover structure are also different.Comparing the top bending moment of the middle column of the negative layer of the subway station at different locations, it is found that the bending moment of the building with upper cover is obviously larger than that of the building without upper Buildings 2023, 13, 2987 9 of 20 cover, and the deformation of the underground structure with upper cover is still shearing type deformation, which is consistent with the conclusion obtained in Reference [6].In addition, the horizontal deformation of the middle column of the underground structure is slightly larger than that of the side wall, and that of the middle column at the position directly connected with the upper building is the largest.
with a peak acceleration of 0.2 g are shown in Figure 10.As can be seen from the bending moment cloud atlas in Figure 10, the top of the negative first floor and the bottom of th negative third floor of each middle column of the subway station have significantly large bending moments compared to the rest of the column.At the same time, the internal forc and deformation of the middle column in the area with and without the upper cover struc ture are also different.Comparing the top bending moment of the middle column of th negative layer of the subway station at different locations, it is found that the bending moment of the building with upper cover is obviously larger than that of the building without upper cover, and the deformation of the underground structure with upper cove is still shearing type deformation, which is consistent with the conclusion obtained in Ref erence [6].In addition, the horizontal deformation of the middle column of the under ground structure is slightly larger than that of the side wall, and that of the middle column at the position directly connected with the upper building is the largest.As the top of the middle column in the negative layer, the bottom of the middle col umn in the negative three layers, and the top and bottom of the side wall in the under ground structure have the largest internal force and deformation, it is the key part to b considered in seismic design.Because of the limited space, only the calculation results o As the top of the middle column in the negative layer, the bottom of the middle column in the negative three layers, and the top and bottom of the side wall in the underground structure have the largest internal force and deformation, it is the key part to be considered in seismic design.Because of the limited space, only the calculation results of these key parts are concerned when the quasi-static seismic calculation of the underground structure is carried out.The monitoring position of the structure when taking the improved response displacement method proposed in this paper for seismic calculation is shown in Figure 11.Where, A is the top of the side wall, B is the bottom of the side wall, C is the top of the middle column of the building without upper cover, D is the bottom of the middle column of the building without upper cover, E is the top of the middle column of the building with upper cover, and F is the bottom of the middle column of the building with upper cover.

One-Dimensional Free-Field Analysis and Bottom Shear Method Calculation Results
The displacements, accelerations, and shear stresses at different locations of the onedimensional free-field soil layers were obtained via Free Field Analysis, the bed coefficients in the horizontal and vertical directions were obtained by the finite element method, and the additional seismic loads on the upper structure, including shear force, bending moment and vertical inertia force, were obtained via the bottom shear method, and the data are shown in Table 2.

One-Dimensional Free-Field Analysis and Bottom Shear Method Calculation Results
The displacements, accelerations, and shear stresses at different locations of the onedimensional free-field soil layers were obtained via Free Field Analysis, the bed coefficients in the horizontal and vertical directions were obtained by the finite element method, and the additional seismic loads on the upper structure, including shear force, bending moment and vertical inertia force, were obtained via the bottom shear method, and the data are shown in Table 2.

Seismic Wave Types and Different Peak Acceleration
The bending moment, deformation, and error of the key parts of the subway station complex and its upper structure system calculated via the traditional response displacement method, the time-history method and the improved response displacement method, under the action of the Kobe wave, the Loma Prieta wave, the Landers wave with peak acceleration of 0.2 g, and the Beijing artificial wave with peak acceleration of 0.05 g, 0.1 g, 0.2 g, and 0.4 g, respectively, are given in Tables 3-6.The deformation error is defined as the ratio of the difference between the numerical results calculated via the response displacement method and the time-history method to the numerical results calculated via the time-history method, the result is a positive number, same as below.It can be clearly seen from Tables 3 and 4 that under the action of different seismic waves, the calculation results of the traditional response displacement method and the timehistory analysis method are quite different, and the calculation results of the traditional response displacement method are generally smaller than the results of the time-history analysis method.The maximum error of structural deformation is 25%, and the bending moment error at the middle column E is greater.Because the improved response displacement method considers the additional seismic load of the upper structure under the seismic action, the calculation results are closer to the real situation, the results of the improved response displacement method and the time-history analysis method are closer, with a maximum error of about 10%.
Meanwhile, under the action of near-field seismic waves and far-field seismic waves, the calculation results of the improved response displacement method are close to the results of the time-history analysis method, so the improved response displacement method can also provide high accuracy for the action of different types of ground shaking.
The internal force and deformation of the side wall and the middle column calculated via the improved response displacement method are increased compared with the traditional response displacement method, and the growth law of the internal force and deformation of the structure are basically the same under different seismic waves.It can be seen from Tables 3 and 5 that the top bending moment of the side wall and the middle column calculated via the improved response displacement method is greatly increased compared with the traditional response displacement method, with an increase of about 10 to 40%.The increase of the bottom bending moment of the side wall and the middle column is small, and the maximum increase is not more than 20%.The bending moment of the middle column with upper cover structure increases more than that without upper cover structure.It can be seen from Tables 4 and 6 that the internal forces of the side wall and the middle column in the improved response displacement method are similar to those in the time-history analysis method, but the displacement of the former is slightly larger than that of the latter.This is due to the fact that the top bending moment of the structural side walls and middle column increases more than the bottom in the improved response displacement method compared to the traditional response displacement method, the significant increase in internal force at the top of the structure makes the horizontal deformation at the top of the structure increase.
From Tables 5 and 6, it can be seen that the change pattern in the calculation results obtained via the three calculation methods is consistent for the subway station complex structure of the TOD mode under the action of the Beijing artificial waves with different ground vibration intensities, and the internal forces and deformations of the subway station structures increase with the peak seismic acceleration.

Analysis of Different Upper Cover Structure Position Forms
To investigate the general applicability of the improved response displacement method, the calculation conditions where the subway station and the upper cover structure are in different relative positions are set up, as shown in Figure 12.Tables 7 and 8 show the internal forces and deformations of the structures of the subway station complexes obtained via the three calculation methods, which differ in their upper cover structure with respect to their relative positions.

Analysis of Different Upper Cover Structure Position Forms
To investigate the general applicability of the improved response displacement method, the calculation conditions where the subway station and the upper cover structure are in different relative positions are set up, as shown in Figure 12.Tables 7 and 8 show the internal forces and deformations of the structures of the subway station complexes obtained via the three calculation methods, which differ in their upper cover structure with respect to their relative positions.It can be seen from Table 7 that with the change in the relative position between the upper cover structure and the lower subway station structure, the bending moment at the side wall A and the middle column C and D decreases first and then increases.The bending moment at the side wall B does not change significantly, and the bending moment at the middle column E and F decreases first, then increases and then decreases.The variation law in the bending moment at different positions of the subway station structure calculated via the improved response displacement method is basically consistent with that calculated by the time-history analysis method, and the maximum error is 12%.However, the change in internal force and deformation of the subway station complex is not obvious, and some of the calculation results obtained via the traditional response displacement method differ greatly from the time analysis method, with the maximum error reaching 206%, which is very unreasonable, during the change in the position of the superstructure frame structure, since the traditional reaction displacement method cannot consider the additional seismic load of the superstructure.

Different Cover Structures
In order to explore the influence of different upper cover structure forms on the calculation accuracy of the improved response displacement method, three working conditions of underground station complex superstructure frame structure and frame-symmetric shear wall structure, and frame-asymmetric shear wall structure are analyzed.
The plan layout of the three upper cover structures forms is shown in Figure 13, where the frame column cross-section size of the frame-shear wall is 0.7 m × 0.7 m, the thickness of the shear wall is 0.5 m, and the beam has a cross-section size of 0.4 m × 0.8 m.It can be seen from Table 7 that with the change in the relative position between the upper cover structure and the lower subway station structure, the bending moment at the side wall A and the middle column C and D decreases first and then increases.The bending moment at the side wall B does not change significantly, and the bending moment at the middle column E and F decreases first, then increases and then decreases.The variation law in the bending moment at different positions of the subway station structure calculated via the improved response displacement method is basically consistent with that calculated by the time-history analysis method, and the maximum error is 12%.However, the change in internal force and deformation of the subway station complex is not obvious, and some of the calculation results obtained via the traditional response displacement method differ greatly from the time analysis method, with the maximum error reaching 206%, which is very unreasonable, during the change in the position of the superstructure frame structure, since the traditional reaction displacement method cannot consider the additional seismic load of the superstructure.

Different Cover Structures
In order to explore the influence of different upper cover structure forms on the calculation accuracy of the improved response displacement method, three working conditions of underground station complex superstructure frame structure and frame-symmetric shear wall structure, and frame-asymmetric shear wall structure are analyzed.
The plan layout of the three upper cover structures forms is shown in Figure 13, where the frame column cross-section size of the frame-shear wall is 0.7 m × 0.7 m, the thickness of the shear wall is 0.5 m, and the beam has a cross-section size of 0. The improved response displacement method model established when the upper structure is a frame-shear wall structure, the parameters of shear wall are determined according to the equivalent stiffness, and that of the frame columns remain the same as the upper frame structure.
The internal forces and deformations of the subway station obtained via different calculation methods for the three structural forms are given in Tables 9 and 10 and Figure 14.The improved response displacement method model established when the upper structure is a frame-shear wall structure, the parameters of shear wall are determined according to the equivalent stiffness, and that of the frame columns remain the same as the upper frame structure.
The internal forces and deformations of the subway station obtained via different calculation methods for the three structural forms are given in Tables 9 and 10 and Figure 14.influence on the bending moment of the sidewall of the subway station, but has more influence on the bending moment at middle columns E and F. This is because the superstructure form of central column E and central column F has changed from a frame structure to a shear wall structure, which changes the form of load transmission, resulting in a large change in the bending moment of some middle columns in the subway station, and the variation law in bending moment is different between the two upper structure forms.In addition, the three-dimensional numerical model based on time-history analysis shows that the bending moment of the column in the subway station complex structure changes along the cross section of the station and the longitudinal section of the station.The models based on the traditional response displacement method and the improved response displacement method are both two-dimensional quasi-static models, and only the variation in the middle column bending moment can be reflected on the cross section of the subway station.Therefore, the calculation results based on the traditional response displacement method and the improved response displacement method are analyzed in Figure 14c; the calculation results of the bending moment at column E and column F of the subway station both show large errors, indicating that the improved response displacement method is not applicable to this upper structure frame-shear wall structure with asymmetric shear wall arrangement.
(   It can be seen that the internal force and deformation of the underground structure are related to the form of the upper structure from Tables 9 and 10 and Figure 14.By comparing the calculation results obtained from the three different upper structure forms using the time-history analysis method, it is found that the upper structure form has less influence on the bending moment of the sidewall of the subway station, but has more influence on the bending moment at middle columns E and F. This is because the superstructure form of central column E and central column F has changed from a frame structure to a shear wall structure, which changes the form of load transmission, resulting in a large change in the bending moment of some middle columns in the subway station, and the variation law in bending moment is different between the two upper structure forms.In addition, the three-dimensional numerical model based on time-history analysis shows that the bending moment of the column in the subway station complex structure changes along the cross section of the station and the longitudinal section of the station.The models based on the traditional response displacement method and the improved response displacement method are both two-dimensional quasi-static models, and only the variation in the middle column bending moment can be reflected on the cross section of the subway station.Therefore, the calculation results based on the traditional response displacement method and the improved response displacement method are analyzed in Figure 14c; the calculation results of the bending moment at column E and column F of the subway station both show large errors, indicating that the improved response displacement method is not applicable to this upper structure frame-shear wall structure with asymmetric shear wall arrangement.

Analysis of Different Relative Positions of Underground Stations to the Ground
Compared with the common subway station, the roof of subway stations with upper structures is closer to the ground, but not necessarily at the same height as the ground.In order to investigate the applicability of the improved response displacement method under different relative positions of underground stations and the ground, three calculation conditions were set up for the top plate of the subway station 4 m above the ground, equal to the ground and 4 m below the ground.The working condition that the top plate of subway station is located 4 m below the ground is that the bottom of the upper structure has been buried into the soil, and the external wall of the upper structure buried in the soil is 0.8 m thick steel and concrete out-wall.Tables 11 and 12 and Figure 15 show the internal forces and deformations of the subway station structure obtained via different calculation methods for three working conditions: 4 m above ground level, equal to ground level, and 4 m below ground level.to the ground and 4 m below the ground.The working condition that the top plate of subway station is located 4 m below the ground is that the bottom of the upper structure has been buried into the soil, and the external wall of the upper structure buried in the soil is 0.8 m thick steel and concrete out-wall.Tables 11 and 12 and Figure 15 show the internal forces and deformations of the subway station structure obtained via different calculation methods for three working conditions: 4 m above ground level, equal to ground level, and 4 m below ground level.
(a) (b) (c)   The calculation results of Tables 11 and 12 and Figure 15 show that in the process of the location of the roof of the subway station from above the ground to level with the ground, and then to below the ground, the bending moments at the A and B positions of the side wall of the subway station and the C and D positions of the middle column become larger, the bending moments at the E and F positions of the middle column increase first and then decrease, and the deformation of the side wall and the middle column decreases.The internal force and deformation of the side wall calculated based on the traditional response displacement method can clearly reflect the above variation law, but the calculation results are small.In addition, the variation law in bending moment at E and F of the middle column is obviously different from that obtained via the time-history analysis method, and the calculation results are also quite different such as the maximum error of deformation is as high as 47%.The deformation of the structure calculated via the modified response displacement method is slightly smaller than that obtained by the time-history analysis method for the working condition in which the top plate position is 4 m above the ground level, unlike the other two working conditions.However, the internal force calculation results obtained by the improved response displacement method are similar to those obtained via the time-history analysis method, which can well reflect the variation in the internal force of the side walls and the middle column.

Different Station Complex Structure Stiffness
In order to analyze the influence of structural stiffness on the applicability of the improved response displacement method, the structural stiffness of the subway station complex is adjusted to 1/2, 1, 2, and 4 times of the original.Tables 13 and 14 show the internal force and deformation results of the subway station structure calculated via different methods under four working conditions.
Analysis of the data in Tables 13 and 14 shows that the increase in structural stiffness leads to an increase in internal forces at different locations.Therefore, it is not appropriate to resist external loads by increasing the structural stiffness only in seismic design.However, the increase in structural stiffness makes the deformation of underground structure decrease, and the structural stiffness can be appropriately increased for the underground structure with strict control of deformation.Compared with the results of the four working conditions obtained via the time-history analysis method, the maximum error of the deformation calculation results of the improved response displacement method is only 11%, and the greater the stiffness, the smaller the error; while the calculation error of the traditional response displacement method is mostly above 20%.

Different Soil Stiffness
Tables 15 and 16 and Figure 16 show the results of internal force and deformation of the structure calculated via different methods with the soil stiffness is adjusted to 1/4, 1/2, 1, and 2 times of the original.Analysis of the data in Tables 13 and 14 shows that the increase in structural stiffness leads to an increase in internal forces at different locations.Therefore, it is not appropriate to resist external loads by increasing the structural stiffness only in seismic design.However, the increase in structural stiffness makes the deformation of underground structure decrease, and the structural stiffness can be appropriately increased for the underground structure with strict control of deformation.Compared with the results of the four working conditions obtained via the time-history analysis method, the maximum error of the deformation calculation results of the improved response displacement method is only 11%, and the greater the stiffness, the smaller the error; while the calculation error of the traditional response displacement method is mostly above 20%.

Different Soil Stiffness
Tables 15 and 16 and Figure 16 show the results of internal force and deformation of the structure calculated via different methods with the soil stiffness is adjusted to 1/4, 1/2, 1, and 2 times of the original.Analysis of Tables 15 and 16 and Figure 16 shows that as the stiffness of the soil increases, the internal force and deformation of the structure tend to become smaller, indicating that the greater the stiffness of the soil, the smaller the external load effect on the underground structure.In engineering, measures to increase the stiffness of the surrounding rock can be used to improve the seismic performance of the structure.The calculation results of the improved response displacement method are slightly smaller than those of the time-history analysis method when the soil stiffness is small, and that of the improved response displacement method are slightly larger when the soil stiffness is larger.The calculation results of internal force and deformation of traditional response displacement method are always small, and the maximum error reaches 26%.Overall, the error of the calculation obtained by the pseudo-static method to calculate the underground structure tend to decrease with the increase in soil stiffness.

Conclusions
An improved response displacement method for the seismic calculation of a subway station complex structure considering the influence of the upper frame structure was proposed based on the theory of the traditional response displacement method, and verified by case analysis, and the main conclusions are as follows: (1) For seismic calculations, both the equivalent base shear method and the modesuperposition response spectrum method are calculated assuming that the structural system is in a linear elastic state, and the structural system in the TOD mode is designed according to the elastic behavior, and the analysis of the internal forces and deformations is carried out using the linear dynamic method.(2) Compared with the calculation results of the time-history analysis method, the internal force and deformation error of a subway station complex calculated via the traditional response displacement method is larger, and the overall value is smaller, so the design is dangerous.The error of internal force and deformation calculated by the improved response displacement method is small, the maximum error is only about 10%, and the variation law is basically consistent with the results in the time-history analysis method.
(3) The near-field seismic wave and far-field seismic wave are selected as the input seismic bedrock motion, respectively, and compared with the calculation results of the time-history analysis method, the improved response displacement method can still provide better accuracy for different types of ground shaking.(4) The traditional response displacement method cannot consider the additional seismic load and vertical inertia force of the superstructure, some of the calculation results obtained via the traditional reaction displacement method are not in conformity with the time-range analysis, and the maximum error is 206%, which makes the calculation results obviously distorted in the process of changing the position of the superstructure frame structure.While the improved reaction displacement method in the calculation results of the underground station structure at different locations of the bending moment change rule and the time-history analysis method is basically consistent, the maximum error is 12%.(5) The deformation of the subway station complex structure under strong seismic action is shear-type deformation in the presence of a upper cover single tower frame structure, and the upper cover structure significantly increases the internal force and deformation of the columns in the underground structure connected to it.The top and bottom of the middle columns and side walls in the structure of the subway station complex of TOD mode, especially the middle column members directly connected to the upper structure, are the weak links in seismic design.(6) The internal forces and deformations of the subway station complex structure calculated via the time-history analysis method, the traditional response displacement method, and the improved response displacement method under different ground vibration intensities increase with the increase of peak acceleration of ground vibration, respectively.

Figure 3 .
Figure 3. Vertical seismic action of upper building.

Figure 4 .
Figure 4. Sketch of calculating the foundation coefficient by the static finite element method.

Figure 3 .
Figure 3. Vertical seismic action of upper building.

Figure 3 .
Figure 3. Vertical seismic action of upper building.

Figure 4 .
Figure 4. Sketch of calculating the foundation coefficient by the static finite element method.

Figure 4 .
Figure 4. Sketch of calculating the foundation coefficient by the static finite element method.

Figure 5 .
Figure 5. Cross section of subway station complex structure with upper frame structure (Unit: m).

Table 1 .Figure 5 .
Figure 5. Cross section of subway station complex structure with upper frame structure (Unit: m).

Figure 7 .
Figure 7. Numerical model for time-history analysis of soil-structure interaction.

Figure 7 .
Figure 7. Numerical model for time-history analysis of soil-structure interaction.

Figure 10 .
Figure 10.Cloud diagram of columns in underground structure.(a) Middle column bending moment cloud diagram.(b) Column deformation cloud (unit: mm).

Figure 10 .
Figure 10.Cloud diagram of columns in underground structure.(a) Middle column bending moment cloud diagram.(b) Column deformation cloud (unit: mm).

Buildings 2023 ,
13,  x FOR PEER REVIEW 10 of 21 these key parts are concerned when the quasi-static seismic calculation of the underground structure is carried out.The monitoring position of the structure when taking the improved response displacement method proposed in this paper for seismic calculation is shown in Figure11.Where, A is the top of the side wall, B is the bottom of the side wall, C is the top of the middle column of the building without upper cover, D is the bottom of the middle column of the building without upper cover, E is the top of the middle column of the building with upper cover, and F is the bottom of the middle column of the building with upper cover.

Figure 11 .
Figure 11.Location of the monitoring point of the structural section.

Figure 11 .
Figure 11.Location of the monitoring point of the structural section.

Figure 13 .
Figure 13.Layout plan of different superstructure structures.(a) Upper frame structure.(b) Upper frame-symmetric shear wall structure.(c) Upper frame-asymmetric shear wall structure.

Figure 13 .
Figure 13.Layout plan of different superstructure structures.(a) Upper frame structure.(b) Upper frame-symmetric shear wall structure.(c) Upper frame-asymmetric shear wall structure.

Figure 14 .
Figure 14.Bending moment calculation results of different upper cover structures.(a) Frame structure.(b) Frame-symmetric shear wall structure.(c) Frame-asymmetric shear wall structure.

Figure 14 .
Figure 14.Bending moment calculation results of different upper cover structures.(a) Frame structure.(b) Frame-symmetric shear wall structure.(c) Frame-asymmetric shear wall structure.

Figure 15 .
Figure 15.Bending moment calculation results of different relative height of subway stations.(a) 4 m above the ground.(b) Equal to the ground.(c) 4 m below the ground.

Figure 15 .
Figure 15.Bending moment calculation results of different relative height of subway stations.(a) 4 m above the ground.(b) Equal to the ground.(c) 4 m below the ground.

Table 3 .
Bending moment of different seismic waves conditions (unit: kN•m).
Note: In this table, TRDM represents the traditional response displacement method, THAM represents the timehistory analysis method, and IRDM represents the improved response displacement method.The representation method in the article is the same below.

Table 4 .
Structural deformation of different seismic wave conditions (unit: mm).

Table 5 .
Bending moment of different earthquake intensities conditions (unit: kN•m).

Table 7 .
Bending moment calculation results of different upper cover structure locations (unit: kN•m).

Table 8 .
Deformation calculation results of different superstructure locations (unit: mm).

Table 7 .
Bending moment calculation results of different upper cover structure locations (unit: kN•m).

Table 8 .
Deformation calculation results of different superstructure locations (unit: mm).

Table 9 .
Bending moment calculation results of different upper cover structures (unit: kN•m).

Table 10 .
Deformation calculation results of different upper cover structures (unit: mm).

Table 9 .
Bending moment calculation results of different upper cover structures (unit: kN•m).

Table 10 .
Deformation calculation results of different upper cover structures (unit: mm).

Table 11 .
Bending moment calculation results of different relative height of subway stations (unit: kN•m).

Table 12 .
Deformation calculation results of different relative height of subway stations (unit: mm).

Table 11 .
Bending moment calculation results of different relative height of subway stations (unit: kN•m).

Table 12 .
Deformation calculation results of different relative height of subway stations (unit: mm).

Table 13 .
Bending moment calculation results of different structural stiffness (unit: kN•m).

Table 14 .
Deformation calculation results of different structural stiffness (unit: mm).

Table 15 .
Bending moment calculation results of different soil layer stiffness (unit: kN•m).

Table 16 .
Deformation calculation results of different soil layer stiffness (unit: mm).