Next Article in Journal
3D Suspended Polymeric Microfluidics (SPMF3) with Flow Orthogonal to Bending (FOB) for Fluid Analysis through Kinematic Viscosity
Next Article in Special Issue
Research on the Rational Yield Ratio of Isolation System and Its Application to the Design of Seismically Isolated Reinforced Concrete Frame-Core Tube Tall Buildings
Previous Article in Journal
From Point Clouds to Building Information Models: 3D Semi-Automatic Reconstruction of Indoors of Existing Buildings
Previous Article in Special Issue
Experimental Study on the Performance of Polyurethane-Steel Sandwich Structure under Debris Flow
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Soil–Structure–Equipment Interaction and Influence Factors in an Underground Electrical Substation under Seismic Loads

1
School of Civil Engineering, Xi’an University of Architecture & Technology, Xi’an 710055, China
2
School of Humanities, Xi’an Polytechnic University, Xi’an 710048, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2017, 7(10), 1044; https://doi.org/10.3390/app7101044
Submission received: 25 July 2017 / Revised: 27 August 2017 / Accepted: 25 September 2017 / Published: 12 October 2017

Abstract

:
Underground electrical substations play an increasingly significant role in urban economic development for the power supply of subways. However, in recent years, there have been few studies on the seismic performance of underground electrical substations involving the interaction of soil–structure–equipment. To conduct the study, three-dimensional finite element models of an underground substation are established. The implicit dynamic numerical simulation analysis is performed by changing earthquake input motions, soil characteristics, electrical equipment type and structure depths. According to a seismic response analysis, acceleration amplification coefficients, displacements, stresses and internal forces are obtained and analyzed. It is found that (1) as a boundary condition of soil–structure, the coupling boundary is feasible in the seismic response of an underground substation; (2) the seismic response of an underground substation is sensitive to burial depth and elastic modulus; (3) the oblique incidence of input motion has a slight influence on the horizontal seismic response, but has a significant impact on the vertical seismic response; and (4) the bottom of the side wall is the seismic weak part of an underground substation, so it is necessary to increase the stiffness of this area.

1. Introduction

In order to meet the requirements of the rapid growth of urban power load under the constraint of densely populated urban land, the number of underground electrical substations has been increasing in recent years. The disaster prevention and mitigation issue of underground structures, including with regard to underground electrical substations, is receiving increasing attention. Therefore, many scholars have focused on the interaction of structure–equipment and the seismic performance of underground structures.
Suarezn and Singh (1989) [1] first developed the interaction principle of primary–secondary structure and obtained response spectrum curves based on the different floors by a simple model involving equipment–structure interaction. Pires (1996) [2] established a simulation model of a San Francisco substation for seismic reliability analysis and compared the numerical substation damages with the real damages caused by Loma Prieta earthquake in 1989. Bazán-Zurita (2009) [3] illustrated concepts interpreting and supplementing the seismic provisions of Manual 113 (Kempner, 2008 [4]) which defined seismic design spectra as simultaneous occurrence of ground motion in three mutually perpendicular directions and estimation of deflections under seismic loading. Poor seismic performance of these system components could result in extended power outages that might propagate far beyond the local epicentral region. Knight and Kempner, Jr. (2009) [5] proposed a Seismic Vulnerability Reduction Program (SVRP) design method to improve the seismic weak parts of substation and present several suggestions on the seismic design of a transformer substation. The international atomic energy agency (IAEA, [6]) provided seismic evaluation criteria on the interaction of structure–equipment in a nuclear power plant. Furthermore, from the study of the interaction between the main plant and electrical equipment (EE) in transformer substations, it was found that the seismic response of the main structure and equipment was not usually synchronized and equipment had a negative effect on the seismic ability of the substation (Wen and Niu, 2013 [7]; Wen et al., 2015 [8]).
To date, many researchers have focused on the interaction of soil–structure. Assuming that interaction force is transformed by a spring and damper, Yuan (1970, [9]) introduced the soil–structure interaction principle. Rodriguez (1985, [10]) analyzed the static interaction between soil and the tunnel lining and found that the force between soil and structure had a significant relationship with contact surface roughness. Penzien (2000, [11]) developed an analytic method involving shear interaction between soil and structure. Based on Penzien’s studies, Huo (2000, [12]) established a hysteretic non-linear soil model to investigate the load transfer mechanisms between underground structure and surrounding soil. Therefore, dynamic finite element analyses were conducted to investigate the load transfer mechanisms between underground structure and surrounding soil. Gong (2002, [13]) studied the responses of an underground structure under vertical seismic excitation through the shaking table test, analyzing the dynamic interaction of soil and underground structure. Based on the background of metro construction in Nanjing city, China, a shaking table test was conducted to model the interaction of soil–metro structure with a liquefiable soil effect. Chen and Zhuang (2006, 2007, [14,15]) studied the dynamic damage characteristics of underground structures in far-field earthquake motions and near-field earthquake motions. According to the static and dynamic coupling effects on the underground structure and surrounding rock system, Du (2009, [16]) proposed stress viscoelastic artificial boundary conditions to establish an underground structure–soil computational model. Tao (2012, [17]) calculated the side wall deformation and structural internal force of a single-story frame of Zhijing line subway in Beijing with different geological conditions and cover depths. In order to understand the distribution of internal force of the underground shield tunnel under earthquake motions, Geng (2013, [18]) investigated the seismic response of underground tunnels in different site types using a quasi-static method and indicated that such a method was more suitable to underground tunnels in soft foundation. Debiasi (2013, [19]) proposed a nonlinear static analysis of a rectangular tunnel by using several different simplified numerical methods, and the qualitative relationship between seismic dynamic strength and structure geometry, overburden thickness and soil stiffness were obtained. Liu (2014, [20]) put forward a Pushover analysis method with good accuracy and reliability, which applied to substructure seismic analysis and design. Kang (2014, [21]) studied the influence of foundation liquefaction on the substructure and conducted an experimental study and numerical analysis on the seismic response of a circular subsurface structure using a modified simple calculating method. Kyriazis (2014, [22]) performed a numerical analysis on an arched tunnel using finite element software and found that some given parameters such as the sectional area of arch, and the thickness of surrounding rock and soil characteristic, had a significant influence on the tunnel seismic response.
From the above research, it is known that underground structure types become more complicated and there are few studies on the seismic performance of underground substations involving the interaction of soil–structure–equipment. In past studies, such interaction in underground substations was usually ignored, and the seismic response caused by characteristics of electrical equipment and soil were neglected. For these reasons, some positive research is proposed in this paper. A three-dimensional finite element model of an underground substation is established to conduct dynamic numerical simulation analysis. By changing some parameters such as ground motions, properties of surrounding soils, electrical equipment and structure depth, the dynamic results of acceleration amplification coefficients, displacements, stresses and internal forces under different earthquake motions are obtained and analyzed.

2. Dynamic Interaction Mechanism of Soil–Structure–Equipment

It is well known that the dynamic equation of structure under earthquake motions can be shown as Equation (1) (Clough and Penzien, 1993, [23]).
[ M ] { u } + [ C ] { u } + [ K ] { u } = { F }
where M is the mass matrix, u is the acceleration vector, C is the damping matrix, u is the velocity vector, K is the stiffness matrix, u is the displacement vector, and F is the load vector.
For the underground substation model, it usually includes three parts: soil, structure and electrical equipment. Similar to Equation (1), the dynamic equation of an underground substation can be described as Equation (2), in which subscripts s, c, and j are denoted as soil, structure, and electrical equipment respectively.
[ M s M c M j ] { u } + [ C s s C s c C c s C c c C j s C j c C s j C c j C j j ] { u } + [ K s s K s c K c s K c c K j s K j c K s j K c j K j j ] { u } = { F s F c F j }
The approximate solution is solved by the Newton Method [24], as shown in Equations (3)–(5).
F N ( u M ) = 0
F N ( u i M + c i + 1 M ) = 0
u i + 1 M = u i M + c i + 1 M
where u M is the exact solution of Equation (2) and u i M is the approximate solution of Equation (2). The dynamic equation numerical solution is solved when a small enough number c i + 1 M exists in Equation (5).

3. Numerical Simulations

3.1. Finite Element Models

The prototype substation used in this study is selected from Urban underground substation design rule (2005, China) [25]. The main plant of the underground substation is a three-story, four-bay reinforced concrete shear wall structure, as shown in Figure 1, Figure 2, Figure 3 and Figure 4. The dimensions of the main plant are 50 m by 25 m, and floor heights from the bottom to top are 4.4 m, 4.8 m, and 5.2 m respectively.
The exterior wall of the main plant structure is made of a reinforced concrete wall, 0.45 m in width. The dimensions of the main beams are 600 × 1300 mm2, 500 × 1400 mm2, 500 × 1000 mm2 separately, and the cross-sectional dimension of the transverse beams is 400 × 1000 mm2. The thicknesses of the first, second, third, and roof slabs are 900 mm, 400 mm, 400 mm, and 600 mm respectively. The total weight of the plant is approximately 27,530 kN.
Ten sets of electrical equipment, shown in Figure 5, are installed on the second floor, the dimensions of which are 6.0 m (length), 2.0 m (width) and 3.5 m (height). The total weight of the equipment is approximately 964 kN, which corresponds to 3.5% of the plant weight.
The plant was designed by Chinese loading Code for Design of Building Structures (2012, [26]), in which gravity loads are selected as 1.8 kN/m2 for floors and 3.5 kN/m2 for the roof. Live loads are selected as 2.5 kN/m2 for floors and 2.0 kN/m2 for the roof. The gravity load of equipment is 8.0 kN/m2. The seismic design of the substation is allowed according to Chinese Code for Seismic Design of Buildings (2010, [27]) and Chinese Code for Seismic Design of Electrical Installations (2013, [28]). The underground substation is assumed to be located in Xi’an, China, where the design basis acceleration associated with 10% probability exceedance in 50 years is 0.2 g (g denotes the coefficient of gravitational acceleration). The compressive strength of concrete and the yield strength of steel are taken as fc = 30 MPa, and fy = 400 MPa, respectively. According to the first eigen-frequency, the factor β is shown in Equation (6).
β = 2 ξ 1 ω 1
Based on a natural frequency extraction analysis of the interaction system, the first eigen-frequency is 0.39 rad/s. Therefore, β is 0.179 s.
The models were developed and analyzed in ABAQUS (2006, [29]). The beams and columns are modeled in three-dimensional beam elements, the floors and concrete walls are modeled in shell elements, and the equipment is modeled in pipe elements, as shown in Figure 6. The three-dimensional finite element electrical equipment model is shown in Figure 7. The concrete plastic-damage model is used for the constitutive properties of concrete material and the double linear dynamic reinforcement model is used for rebar. Druker–Prager, a generalized Mises yield criterion, is adopted in soil dilation. The effect of hydrostatic pressure on yield and strength is used as the constitutive model of soil (Liu et al., 2006, [30]). The physical and mechanical parameters of soil, as determined by the penetration test and consolidation test, are listed in Table 1. Based on the artificial constrained boundary, the boundary element is used to simulate the infinite field of soil, and the horizontal reaction forces on the lateral boundary are determined by static analysis. According to the literature (Liu et al., 2006 [31]; Huang et al., 2010 [32]), the coupling analysis method is used in the soil boundary condition. Meanwhile, it is suggested from literature (Lou et al., 1999 [33]; Lou et al., 2000 [34]) that five times the width of the structure can be defined as the width of the computational area of soil in order to reduce the motion reflection on the boundary. Following this rule, five times the width of the structure is selected as the soil width in the finite element model. Therefore, the soil size of the finite element is 250 m × 132.5 m × 50 m and the size of the entire soil is 500 m × 265 m × 50 m. The arbitrary Lagrangian Eulerian (ALE) adaptive meshing method is used to maintain a high-qualified meshing system for soil, which allows the meshes to move independently when large deformation of liquefied soils occurs. The soil element meshes are shown in Figure 8, in which the finite element of soil is simulated by unit C3D8R. In order to simulate the most realistic condition of soil and simplify some less salient factors in this paper, the following basic assumptions are advised in analysis (Chen [35], 2015; Novak [36], 1972).
(1)
Each layer of soil is homogeneous, and extends infinitely in the horizontal direction;
(2)
There is no relative sliding among different layers of soil;
(3)
Earthquake excitations are supposed from the bottom of the foundation;
(4)
The movement of each point on the bottom of the foundation is consistent without a traveling motion effect.
According to the classification of site type and seismic fortification intensity, three earthquake records—El-Centro record, Taft record and Lanzhou record—were derived from the registered seismic data. Moreover, in order to obtain earthquake motions, three earthquake records were back-calculated to the bottom of the foundation level by the SHAKE91 procedure (Idriss et al., 1992 [37]). According to the finite element division of the soil layer, soil properties such as density, shear velocity, shear modulus of elasticity, etc., were taken into account in earthquake motion inversion. Time-histories and spectra of these records after inversion are shown in Figure 9, Figure 10, Figure 11, Figure 12, Figure 13 and Figure 14.

3.2. Accuracy Verification of the Finite Element Model

The soil boundary condition of the coupling analysis method of the finite element and infinite element, which reasonably reflects the real boundary condition, gives a unified solving format, and makes the infinite element of the soil section become part of the finite element model (Liang et al., 2010 [38]; Meng et al., 2012 [39]). Compared with the viscoelastic boundary condition, the infinite element boundary condition can not only simulate the far field boundary to absorb seismic energy, but also correctly simulate the boundary conditions whose infinity displacement is zero. Therefore, the coupling analysis method of the finite element and infinite element is used as the soil boundary condition to simulate semi-infinite domain soil.
For the three-dimensional finite element model, the distance value r is first calculated from the random boundary node to the structure center. Secondly, normal spring stiffness, normal damping coefficient, tangential spring stiffness, and tangential damping coefficient of the node on the boundary are calculated by Equations (7) and (8). Finally, all of the boundary nodes are calculated cyclically in order to obtain a viscoelastic artificial boundary.
K l i = λ + 2 G r   C l i = ρ c P
K l i = G r   C l i = ρ c S
where cp is the longitudinal velocity of the medium, cs is the longitudinal velocity of the medium, ρ is density, λ and G are lame constants, and r is the distance from the artificial boundary point to the scattering source.
The soil surface is supposed to be a free boundary, and the horizontal reaction forces at the nodes on the lateral boundary and bottom boundary are calculated according to the artificially constrained boundary assumption (Liu et al., 2006 [31]; Huang et al., 2010 [32]). Therefore, the fixed boundary model, viscoelastic boundary model and coupling analysis model are set up to verify different soil boundary conditions, which are simplified as follows:
  • Model 1 (M1): Coupling boundary between finite elements and infinite elements
  • Model 2 (M2): Viscoelastic boundary
  • Model 3 (M3): Fixed boundary
The acceleration time–history curves at the soil surface midpoint are shown in Figure 15.
From Figure 15, it can be inferred that the curves are almost the same in three different conditions before 4 s. After 4 s, only curves of Model 1 and Model 2 are similar to each other, while the acceleration values of Model 3 gradually increase compared to those of the other two curves. The reason is that seismic energy is absorbed and the seismic wave is rarely reflected when the motions pass through the infinite element boundary or viscoelastic artificial boundary, while seismic energy is greatly reflected when the motions pass through the fixed boundary, which leads to the accumulation of acceleration values on the free boundary. Therefore, compared with the viscoelastic boundary condition, the infinite element boundary condition does not need to be involved in the analytical solution expression in ABAQUS software. For this reason, Model 1 is selected as the soil boundary condition in this paper.

4. Numerical Simulation Analysis

The underground substation model involving interaction of soil–structure–equipment is established for numerical simulation analysis. Based on finite element analysis, the influence of some parameters on seismic response such as oblique incidence of input motion, properties of surrounding soils, and burial depth of the underground substation are analyzed.

4.1. The Seismic Response under Different Ground Motions

According to Chinese Code for Seismic Design of Buildings (2010, China [27]), the peak ground acceleration (PGA) is given in Table 2, in which the scales of PGA in X, Y and Z directions are adjusted to 0.20 g, 0.17 g, and 0.13 g, respectively.

4.1.1. Acceleration Responses

The acceleration amplification factors under 0.20 g are shown in Table 3, in which the acceleration amplification factor is defined as the ratio of the peak acceleration of the layer over the peak acceleration of input motion.
The following can be inferred from Table 3:
(1)
The acceleration amplification factors in the X-direction are larger than those in the Y-direction in M1 and M3, which means that the seismic response of the underground substation in the X-direction is more severe than that in the Y-direction. The seismic response including acceleration responses, displacement responses and internal force responses of 0.17 g PAG in the X-direction and 0.2 g PAG in the Y-direction would be reduced but the response rule is uniform.
(2)
The acceleration amplification factors in the Z-direction of M3 are larger than those in M2. The reason is that the seismic response is magnified under multi-direction earthquake motions compared with that under the one-direction earthquake input motion and the horizontal earthquake motions have a reverse effect on the seismic response in the Z-direction.
(3)
The peak vertical and horizontal acceleration amplification factors are 1.526 and 0.760, respectively, which means that the oblique incidence of input motion has a slight influence on the horizontal seismic response, but has a significant impact on the vertical seismic response.

4.1.2. Displacement Responses

The horizontal deformation of the side wall causes eccentric compression and a large bending moment, which can lead to damage of the side wall. Therefore, the horizontal layer drift of the underground substation is the major factor of side wall damage, and it is necessary to analyze the horizontal relative deformation of the underground substation under different ground motions.
The relative displacement time–history curves of the top floor under El-Centro motion and Lanzhou motion are shown in Figure 16, in which the relative displacement is defined as horizontal displacement values, and layer drifts in M1 and M3 under 0.20 g are shown in Table 4.
It can be seen from Figure 16 and Table 4 that layer drifts in the Y-direction are larger than those in the X-direction, which coincides with the feature that stiffness in X-direction is larger than that in the Y-direction. There is only a 2% difference in the value of the drifts between one-way input motions and multi-way input motions in the X-direction, which indicates that there is no obvious difference in torsion coupling in the X-direction. On the other hand, in the Y-direction, the horizontal drifts with multi-way input motions are about 1.4 times larger than those with one-way input earthquake motion. Therefore, the drifts of the corner column under multi-way input motions are larger than those under one-way input motions, which leads to severe damage of the corner columns. It is essential to improve the constraints of corner columns and increase the number of stirrups in these columns, which can strengthen the deformation ability of such components.

4.1.3. Internal Force Responses

The shear wall is an important component under seismic action. As the first defense of energy dissipation under earthquake motions, coupling beams cannot bear the internal force value owing to the limited sectional dimension. Thus, the main internal force value is borne by the shear wall after internal force redistribution. There are six edges of concrete walls represented by Q-2, Q-7, Q-12, Q-17, Q-22 and Q-27 on surface-6, as shown in Figure 17. Furthermore, there is a similar seismic response among different sections. Therefore, the surface-6 section is selected to analyze the internal force response of the side walls. The peak internal forces in each section of surface-6 under 0.20 g are given in Table 5 and Figure 18.
The following can be concluded from Figure 18 and Table 5:
(1)
Different input motions have different influences on the internal force of side walls. The shear forces of components are similar between horizontal direction motion and coupling input motion. Furthermore, the vertical pressure and bending moments of components under vertical seismic motion are similar to those in coupling input motion. Oblique incidence motion makes the vertical pressure and bending moments increase obviously.
(2)
In comparing the peak value of internal force in different sections, it can be seen that the shear force and vertical pressure of the bottom of the side wall on the third floor are, at most, about 1.03 × 106 N/m and 3.67 × 106 N/m respectively, values which are significantly greater than those on the other floors. So, the bottom of the side walls is the seismic weak part of the structure.
(3)
The bending moments at joint parts of the side wall on the second floor are larger than the rest of the structure, because the extra inertia force caused by the electrical equipment leads to the second floor being in a complex bend torsion condition. The bending moments of the second floor under coupling seismic action increase by 20% compared to those under horizontal seismic action. So, it can be speculated that the vertical pressure and bending moments increase in vertical seismic motion.

4.2. Soil Properties in the Seismic Response of the Underground Substation

In order to analyze the seismic response of soil and the underground substation, M1 (1 time soil elasticity modulus model), M2 (2 times soil elasticity modulus model) and M3 (3 times soil elasticity modulus model) were set up separately, the seismic response results of which are shown in Table 6.

4.2.1. Displacement Responses

The relative displacements under 0.20 g are shown in Figure 19. The layer drifts and layer drift angles under 0.20 g are shown in Table 7.
(1) From Figure 19, it can be inferred that the displacements of the side wall decrease with the increasing of the elastic modulus of soil and the horizontal drift diminishes mostly at the top of the side wall, and the horizontal relative displacement in the Y-direction decreases by about 73.9% in Lanzhou motion. The reason is that, with the increasing of the soil elastic modulus, the structure constraint increases and structural deformation decreases, which usually lead to the diminishing of horizontal drift.
(2) From Table 7, the drift angle is inversely proportional to the elastic modulus of surrounding soil. With the increase in the soil elastic modulus, the drift angle of each layer significantly reduces; the minimum value decreases sharply from 1/1611 to 1/7059. Therefore, the soil elastic modulus has a very important influence on the deformation of the underground substation.

4.2.2. Stress Responses

The stress responses of the underground substation are similar under El-Centro, Taft and Lanzhou earthquakes; therefore, the stress response of the structure under Lanzhou earthquake is taken as an example. The stress nephograms of the underground substation under 0.20 g are shown in Figure 20. The peak stresses and stress amplitude in different locations under 0.20 g are shown in Table 8.
With the increase in the elastic modulus of soil, stress amplitudes of layers and side walls do not simply decrease with the decrease of layer drift, but firstly decrease and then increase. The second-floor stress amplitude maximum value is 46.23% when the soil elastic modulus increases from two- to three-fold. The reason is that with the increase in the elastic modulus of soil, the deformation of soil reduces and the corresponding displacement of the structure reduces as well, which leads to a decrease in structure stress to some extent. However, some adverse factors, such as non-uniform stiffness, big holes in floors, large weight and the volume of electrical equipment, can lead to a nonlinear response in the internal forces.

4.2.3. Internal Force Responses

The peak internal forces of each section on surface-6 under different seismic motions are shown in Table 9 and the peak internal forces of different points on surface-6 are shown in Figure 21.
From Table 9 and Figure 21, it can be inferred that the shear forces of the underground substation gradually increase with the increase in the elastic modulus of the surrounding soil, while vertical pressure and bending moments change little under the same condition. The reason is that soil shear deformation reduced with the increase in the soil elastic modulus, and the corresponding displacements of the structure decreased as well.
The change of the soil’s physical properties has a great effect on the seismic response of the underground substation. Therefore, the interaction between the underground substation and soil should be taken into account in the seismic design.

4.3. Seismic Response Influence of Structure–Electrical Equipment Interaction

Due to its heavy weight and huge size, electrical equipment fixed on the floors of the substation usually produced additional inertial force, which had an additional effect on the structure under seismic load and bending torsional load.
In previous seismic designs of underground substations, the interaction between structure and equipment was usually ignored. When the size and weight of electrical equipment is low, such a method is acceptable. Nowadays, the size and weight of electrical equipment is larger than before, and such an equivalent load method would result in some error. For this reason, Model 1 (M1) of the underground substation with the interaction between electrical equipment and structure and Model 2 (M2) without the interaction between them were established. According to a comparison and analysis of the internal force and seismic response in these two models, the seismic performance effect of electrical equipment on the underground substation can be determined.

4.3.1. Acceleration Responses

The acceleration amplification factors under 0.40 g, with the interaction between electrical equipment and structure, are shown in Table 10, in which the acceleration amplification factor is defined as the ratio between the peak acceleration of both the layer and input seismic motion. Similarly, the soil acceleration amplification factor is defined as the value at the same level of the structure.
(1) Under the same direction and earthquake motion, acceleration amplification factors of soil and structure gradually increase; acceleration amplification factors in the horizontal X-direction are greater than those in the horizontal Y-direction. Meanwhile, the peak acceleration amplification factors of structure and soil are 1.070 and 1.772 (in the horizontal X-direction under Taft motion) respectively, which means that seismic responses are more severe in these cases than in other cases.
(2) Under different earthquake motions, acceleration amplification factors at the top of electrical equipment are different from each other. The peak acceleration amplification factor on the second floor under Taft motion is 3.52 (in the horizontal or the X-direction) which is much greater than 2.0 given in Code for the seismic design of electrical installations (2013, China). The acceleration amplification factors of electrical equipment in the X-direction are obviously greater than those on the same floor, which suggests that the acceleration response of the structure and electrical equipment is not synchronous. Therefore, the interaction between structure and equipment should be taken into account in the seismic design of an underground substation.

4.3.2. Displacement Responses

Time–history curves of displacement on the top floor in M1 and M2 under El-Centro motion and Lanzhou motion are shown in Figure 22. Layer drift and layer drift angle of each layer under 0.40 g are shown in Table 11.
From Figure 22 and Table 11, it can be seen that layer drift and layer drift angles of M1 are larger than those of M2 under 0.40 g. The reason is mainly that electrical equipment produces additional inertial force and extra relative displacement, which leads to structure displacement increasing at the same time. Therefore, the interaction of structure–equipment should be involved in seismic design. Only when the effect of electrical equipment is taken into account is the seismic design of an underground substation safe.

4.3.3. Internal Force Responses

Shear force and bending moment time–history curves in the X-direction at the bottom of side wall 4 under 0.40 g are shown in Figure 23. The peak internal forces of the side wall under 0.40 g are shown in Table 12.
From Figure 23 and Table 12, it is shown that the difference of internal force is obvious in M1 and M2, in which the peak shear amplification is 18.99% from M2 to M1 and the peak bending moment difference increased by 19.93%. Such data reaffirm that the interaction between electrical equipment and structure could not be ignored.

4.4. Influence of Burial Depth

Burial depth is one of the important factors in the seismic response of the underground structure. The same underground structures may have different seismic responses in different burial depths. To study the influence of different burial depths on the seismic response of an underground substation in the horizontal direction, five conditions including different burial depth at 1.2 m (M1), 3 m (M2), 5 m (M3), 7 m (M4) and 9 m (M5) were chosen for the samples.

4.4.1. Displacement Responses

Layer drifts and drift angles under 0.20 g are shown in Table 13 under El-Centro motion.
(1) With the increase in buried depth, the relative horizontal displacements in the X-direction first increase and then decrease, and then reach the maximum values when the buried depth is 7.0 m. Simultaneously, the relative horizontal displacements in the Y-direction first decrease and then increase and reach the maximum values when the buried depth is 1.2 m. The layer drift of each floor has the same trend and the maximum layer drift angle value in the X-direction reaches 1/809 (the third underground floor when the buried depth is 7.0 m), which is close to the limit value of 1/800 stated in Code for seismic design of buildings (2010, China).
(2) It is usually believed that the deeper the structures are buried, the safer they are under earthquake motions. However, the rule is unsuitable for some special industrial structures such as an underground substation. The reasons are that equipment of high quality and stiffness with complex shapes makes displacements of structure complicated, and overlaying soil reduces the overall displacement of the structure with the increase in buried depth. The layer drift is critical in the Y-direction when the buried depth is shallow, while it is critical in the X-direction with the increase in buried depth. In order to prevent the danger caused by excessive displacement, buried depths which can cause a large seismic response should be avoided in underground substation seismic design.

4.4.2. Internal Force Responses

It has been previously demonstrated that the side wall is the weak area in the structure; therefore, shear force of side wall 3 is chosen to be analyzed and the results are shown in Table 14, in which the static values are defined as static computational results under the effect of gravity, and dynamic values are defined as dynamical computational results under earthquake motions.
Comparing shear force amplifications at the top and bottom of side walls in different buried depths, it can be concluded that dynamic load has a greater influence on the top of side walls when the buried depth is 1.2 m when the shear force amplifications of the top of the side wall are larger than those of bottom. With the increase of buried depths, dynamic loads have a greater effect on the base of side walls, and the shear force amplifications on the top of the side wall are less than those of the bottom, which causes shear forces within the bottom of the wall to reach the maximum values under static and dynamic loads. Thus, the bottom of the side wall is the weak part of the underground substation. Therefore, it should be given more attention during the seismic design of such a system.

5. Conclusions

Three-dimension finite element models involving the interaction of soil–structure–equipment are established, and dynamic numerical simulation analysis of such models is performed. By changing some parameters such as ground motion input motions, properties of surrounding soils, burial depth, etc., the seismic performance of underground substations is proposed, and simulated results, such as acceleration amplification coefficient, displacement, stress and internal forces, are obtained from research. The findings are summarized as follows:
(1)
The coupling boundary is selected as the boundary condition of soil–structure, which is feasible in a seismic response of an underground substation.
(2)
The seismic response of an underground substation is more sensitive to vertical earthquake motion than the seismic response of normal structures. Therefore, vertical earthquake motions should be taken into account in the seismic design of such underground structures.
(3)
Burial depth and elastic modulus are the main factors for the seismic performance of the underground electrical substation. With the increase in burial depth, layer drifts of an underground substation first increase and then decrease. With the increase in the elastic modulus of soil, the constraint of soil on structure increases and the deformation of structure decreases.
(4)
The acceleration amplification factors of electrical equipment are obviously greater than those of the same layer of the structure, and the peak drift of equipment is not synchronous with structure. Therefore, the interaction between structure and equipment should be taken into account in the seismic design of underground substations and some seismic measures should be proposed to control the dynamic response of electrical equipment.
(5)
The bottom of the side walls is the weak part of an underground substation. Therefore, it is necessary to increase the stiffness of side walls and strengthen the connection between the bottom floor and side walls.

Acknowledgments

The authors wish to express their gratitude towards the Major Program of the National Natural Science Foundation (No. 51590914) (in China), the National Natural Science Foundation Project (No. 51578450, No. 51378225) (in China).

Author Contributions

Bo Wen and Lu Zhang conceived and designed the models; Lu Zhang analyzed the data; Ditao Niu contributed analysis tools; Bo Wen wrote the paper; and Muhua Zhang searched some references.

Conflicts of Interest

The authors declare no conflict of interest. And the founding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, and in the decision to publish the results.

References

  1. Suarez, L.E.; Singh, M.P. Floor Spectra with Equipment Structure Equipment Interaction Effects. J. Eng. Mech. 1989, 115, 247–264. [Google Scholar] [CrossRef]
  2. Pires, J.A.; Ang, A.H.-S.; Villaverde, R. Seismic reliability of electrical power transmission systems. Nucl. Eng. Des. 1996, 160, 427–439. [Google Scholar] [CrossRef]
  3. Brzan-Zurita, E.; Bielak, J.; Digioia, A.M. Seismic Design of Substation Structures. In Proceedings of the Electrical Transmission and Substation Structures Conference ASCE, Fort Worth, TX, USA, 8–12 November 2009; pp. 1–12. [Google Scholar]
  4. Kempner, L.; Reston, V.A., Jr. (Eds.) Substation Design Guide; ASCE: Oakland, CA, USA, 2008. [Google Scholar]
  5. Knight, B.T.; Kempner, L. Vulnerabilities and Retrofit of High-Voltage Electrical Substation Facilities. In Proceedings of the Technical Council on Lifeline Earthquake Engineering Conference, Oakland, CA, USA, 28 June–1 July 2009; pp. 232–243. [Google Scholar]
  6. International Atomic Energy Agency (IAEA). Seismic Evaluation of Existing Nuclear Installations; IAEA: Vienna, Austria, 2008. [Google Scholar]
  7. Wen, B.; Niu, D.T. Seismic vulnerability analysis for the main building of the large substation. China Civ. Eng. J. 2013, 46, 19–23. (In Chinese) [Google Scholar]
  8. Wen, B.; Taciroglu, E.; Niu, D.T. Shake table testing and numerical analysis of transformer substations including main plant and electrical equipment interaction. Adv. Struct. Eng. 2015, 18, 1959–1980. [Google Scholar] [CrossRef]
  9. Yuan, H.F.; Walker, R.E. The Investigation of a Simple Soil-Structure Interaction Model. In Dynamic Waves in Civil Engineering; Howells, D.A., Haigh, I.P., Taylor, C., Eds.; Wiley: New York, NY, UAS, 1970; pp. 247–266. [Google Scholar]
  10. Rodriguez, B. Numerical Analysis of an Experimental Tunnel. Proceeding of the 11th International Conference on Soil Mechanics and Foundation Engineering San Francisco, San Francisco, CA, USA, 12–16 August 1985; pp. 789–792. [Google Scholar]
  11. Penzien, J. Seismically Induced Racking of Tunnel Linings. Earthq. Eng. Struct. Dyn. 2000, 29, 683–691. [Google Scholar] [CrossRef]
  12. Huo, H. Seismic Design and Analysis of Rectangular Underground Structures. Tunneling and Underground Space Tecnology United Kingdom. Ph.D. Thesis, Purdue University, West Lafayette, IN, USA, 2000. [Google Scholar]
  13. Gong, B.N.; Zhao, D.P. Experimental Study on Dynamic Interaction of Underground Structure and Soil. Undergr. Space 2002, 22, 320–324. (In Chinese) [Google Scholar]
  14. Zhuang, H.Y. Study on Dynamic Contact Properties of Soil-Subway Underground Structure Interaction System; Nanjing University of Technology: Nanjing, China, 2006. (In Chinese) [Google Scholar]
  15. Chen, G.X.; Zhuang, H.Y.; Du, X.L. Analysis of large-scale shaking table test of dynamic soil-subway station interaction. Earthq. Eng. Eng. Vib. 2007, 27, 171–176. (In Chinese) [Google Scholar]
  16. Li, L.Y.; Du, X.L.; Li, L. Review on Mechanical Performanical of Soil-Structure Contact. Adv. Mech. 2009, 39, 588–597. (In Chinese) [Google Scholar]
  17. Tao, L.J.; Wang, W.P.; Zhang, B. Difference law study of seismic design methods for subway structures. China Civ. Eng. J. 2012, 45, 170–176. (In Chinese) [Google Scholar]
  18. Geng, P.; He, C.; Tang, J.; Yan, Q.; Xu, Y. Appropriate Seismic Analysis Method of Subway Shield Tunnels in Soft Ground. Int. Efforts Lifeline Earthq. Eng. 2013, 283–290. [Google Scholar] [CrossRef]
  19. Debiasi, E.; Gajo, A.; Zonta, D. On the Seismic Response of Shallow-Buried Rectangular Structures. Tunn. Undergr. Space Technol. 2013, 38, 99–113. [Google Scholar] [CrossRef]
  20. Liu, J.B.; Wang, W.H.; Dasgupta, G. Pushover analysis of underground structures: Method and Application. Sci. China Technol. Sci. 2014, 57, 423–437. (In Chinese) [Google Scholar] [CrossRef]
  21. Kang, G.H.; Tobita, T.; Lai, S. Seismic simulation of liquefaction-induced uplift behavior of a hollow cylinder structure buried in shallow ground. Soil Dyn. Earthq. Eng. 2014, 64, 85–94. [Google Scholar] [CrossRef]
  22. Pitilakis, K.; Tsinidis, G.; Leanza, A.; Maugeri, M. Seismic behavior of circular tunnels accounting for above ground structures interaction effects. Soil Dyn. Earthq. Eng. 2014, 67, 1–15. [Google Scholar] [CrossRef]
  23. Clough, R.; Penzien, J. Dynamics of Structures; McGraw-Hill: New York, NY, USA, 1993. [Google Scholar]
  24. Zhang, F.; Cheng, C.S. A modified Newton method for radial distribution system power flow analysis. IEEE Trans. Power Syst. 2002, 12, 389–397. [Google Scholar] [CrossRef]
  25. 35kV–220kV Urban Underground Substation Design Rules; China Electric Power Press: Beijing, China, 2005. (In Chinese)
  26. Chinese Loading Code for Design of Building Structures; GB50009-2012; Architecture & Building Press: Beijing, China, 2012. (In Chinese)
  27. Chinese Code for Seismic Design of Buildings; GB50011-2010; Architecture & Building Press: Beijing, China, 2010. (In Chinese)
  28. Chinese Code for Seismic Design of Electrical Installations; GB50260-2013; Architecture & Building Press: Beijing, China, 2013. (In Chinese)
  29. Hibbitt, K.; Sorensen, I. ABAQUS User Subroutines Reference Manual; Elsevier Science Publishers: Philadelphia, PA, USA, 2006. [Google Scholar]
  30. Liu, J.L.; Luan, M.T.; Xu, C.S. Study on parametric characters of Drucker-prager criterion. Chin. J. Rock Mech. Eng. 2006, 25, 4009–4015. (In Chinese) [Google Scholar]
  31. Liu, J.B.; Gu, Y.; Du, Y.X. Consistent viscous-spring artificial boundaries and viscous-spring boundary elements. Chin. J. Geotech. Eng. 2006, 28, 1070–1075. (In Chinese) [Google Scholar]
  32. Huang, S.; Chen, W.Z.; Wu, G.J. Study of method of earthquake input in a seismic analysis for underground engineering. Chin. J. Rock Mech. Eng. 2010, 29, 1254–1262. [Google Scholar]
  33. Lou, M.L.; Chen, H.J. Research on the Effects of Lateral Boundary of Pile Seismic Response; Tongji University Press: Shanghai, China, 1999. (In Chinese) [Google Scholar]
  34. Lou, M.L.; Wang, W.J.; Zhu, T. Soil lateral boundary effect in shaking table model test of soil-structure system. Earthq. Eng. Eng. Vib. 2000, 20, 30–36. (In Chinese) [Google Scholar]
  35. Chen, G.X.; Chen, S.; Zuo, X. Shaking-table tests and numerical simulations on a subway structure in soft soil. Soil Dyn. Earthq. Eng. 2015, 76, 13–28. [Google Scholar]
  36. Novak, M.; Bereduge, Y.O. Vertical Vibration of Embedded Footings. J. Soil Mech. Found. Div ASCE 1972, 98, 1291–1310. [Google Scholar]
  37. Idriss, I.M.; Sun, J. User’s Manual for SHAKE91. Center Geotech. Model. 1992, 388, 279–360. [Google Scholar]
  38. Liang, X.W.; Qian, L.; Tan, L.N. The Concrete Damage Plastic Constitutive Relation Research Based on ABAQUS. In Proceedings of the 8th Session of the National Earthquake Engineering, Chongqing, China, 24–26 December 2010; Volume 32, pp. 646–648. (In Chinese). [Google Scholar]
  39. Meng, W.Y.; Wang, J.F.; Zhang, R. Comparison Analysis of Constitutive Models of Reinforced Concrete Structures Based on ABAQUS. J. North China Inst. Water Conserv. Hydroelectr. Power 2012, 33, 40–42. (In Chinese) [Google Scholar]
Figure 1. Plan of the second floor.
Figure 1. Plan of the second floor.
Applsci 07 01044 g001
Figure 2. Dimensions of the floor and wall.
Figure 2. Dimensions of the floor and wall.
Applsci 07 01044 g002
Figure 3. Long cross-section.
Figure 3. Long cross-section.
Applsci 07 01044 g003
Figure 4. Short cross-section.
Figure 4. Short cross-section.
Applsci 07 01044 g004
Figure 5. Electrical equipment.
Figure 5. Electrical equipment.
Applsci 07 01044 g005
Figure 6. Finite element model of a substation.
Figure 6. Finite element model of a substation.
Applsci 07 01044 g006
Figure 7. Finite element model of equipment.
Figure 7. Finite element model of equipment.
Applsci 07 01044 g007
Figure 8. Finite element model of soil.
Figure 8. Finite element model of soil.
Applsci 07 01044 g008
Figure 9. El-Centro motion in the horizontal direction after inversion.
Figure 9. El-Centro motion in the horizontal direction after inversion.
Applsci 07 01044 g009
Figure 10. El-Centro motion in the vertical direction after inversion.
Figure 10. El-Centro motion in the vertical direction after inversion.
Applsci 07 01044 g010
Figure 11. Taft motion in the horizontal direction after inversion.
Figure 11. Taft motion in the horizontal direction after inversion.
Applsci 07 01044 g011
Figure 12. Taft motion in the vertical direction after inversion.
Figure 12. Taft motion in the vertical direction after inversion.
Applsci 07 01044 g012
Figure 13. Lanzhou motion in the horizontal direction after inversion.
Figure 13. Lanzhou motion in the horizontal direction after inversion.
Applsci 07 01044 g013
Figure 14. Lanzhou motion in the vertical direction after inversion.
Figure 14. Lanzhou motion in the vertical direction after inversion.
Applsci 07 01044 g014
Figure 15. Time–history curves of the midpoint’s acceleration of soil.
Figure 15. Time–history curves of the midpoint’s acceleration of soil.
Applsci 07 01044 g015
Figure 16. Relative displacement curves of the top floor.
Figure 16. Relative displacement curves of the top floor.
Applsci 07 01044 g016aApplsci 07 01044 g016b
Figure 17. Position of the observation surface and point.
Figure 17. Position of the observation surface and point.
Applsci 07 01044 g017
Figure 18. The peak internal forces of surface-6 under El-Centro motion.
Figure 18. The peak internal forces of surface-6 under El-Centro motion.
Applsci 07 01044 g018
Figure 19. Relative horizontal displacement of the side wall.
Figure 19. Relative horizontal displacement of the side wall.
Applsci 07 01044 g019
Figure 20. Stress nephograms of the underground substation.
Figure 20. Stress nephograms of the underground substation.
Applsci 07 01044 g020
Figure 21. The peak internal force of each section on surface-6 under El-Centro motion.
Figure 21. The peak internal force of each section on surface-6 under El-Centro motion.
Applsci 07 01044 g021
Figure 22. Time–history curves of displacement on the top floor.
Figure 22. Time–history curves of displacement on the top floor.
Applsci 07 01044 g022
Figure 23. Shear force and bending moment time–history curves in the X-direction at the bottom of side wall 4.
Figure 23. Shear force and bending moment time–history curves in the X-direction at the bottom of side wall 4.
Applsci 07 01044 g023
Table 1. The physical and mechanical soil parameters.
Table 1. The physical and mechanical soil parameters.
NumberName of SoilDensity (kg/m3)Elastic Modulus (MPa)Poisson RatioCohesion(kPa)Internal Friction Angle (°)Soil Depth (m)
1Plain fill17205.00.3922154.30
2New loess160012.00.3028186.20
3Silty clay195015.00.3045151.70
4Pebble bed225055.00.150452.10
5Silty clay195015.00.30451527.0
6Coarse ands192040.00.260354.50
Table 2. The conditions of peak ground acceleration (PGA).
Table 2. The conditions of peak ground acceleration (PGA).
ModelDirectionAcceleration Input Value
M1Horizontal0.2 g:0.17 g (X:Y)
M2Vertical0.13 g (Z)
M3Coupling of horizontal and vertical0.2 g:0.17 g:0.13 g (X:Y:Z)
Table 3. Acceleration amplification factors under 0.20 g.
Table 3. Acceleration amplification factors under 0.20 g.
SectionEL-Centro MotionLanzhou Motion
M1M2M3M1M2M3
RoofX-direction0.671/0.6830.683/0.714
Y-direction0.468/0.4800.489/0.502
Z-direction/1.4761.526/1.1081.163
Floor 3X-direction0.677/0.6950.699/0.726
Y-direction0.573/0.5820.503/0.520
Z-direction/1.4761.525/1.1081.162
Floor 2X-direction0.679/0.7190.722/0.742
Y-direction0.595/0.6230.536/0.565
Z-direction/1.4751.523/1.0161.159
Floor 1X-direction0.708/0.7210.746/0.760
Y-direction0.681/0.6920.568/0.590
Z-direction/1.4751.522/1.1051.158
Table 4. Layer drifts (mm).
Table 4. Layer drifts (mm).
SectionEl-Centro MotionLanzhou Motion
M1M3M1M3
X-DirectionY-DirectionX-DirectionY-DirectionX-DirectionY-DirectionX-DirectionY-Direction
Floor 31.652.471.764.031.622.181.773.14
Floor 21.562.301.653.811.532.001.652.98
Floor 11.702.501.793.831.682.101.772.99
Table 5. The peak internal forces of each section on surface-6.
Table 5. The peak internal forces of each section on surface-6.
SectionVertical Pressure (105 N/m)Shear Force (105 N/m)Bending Moment (103 N·m/m)
M1M2M3M1M2M3M1M2M3
Floor 3El-CentroTop3.744.324.435.414.745.427.198.228.28
Bottom4.875.205.476.154.826.176.546.947.16
LanzhouTop3.753.934.015.614.365.647.307.878.01
Bottom4.514.684.866.144.576.236.907.217.47
Floor 2El-CentroTop6.266.767.147.495.127.506.337.677.85
Bottom9.3910.1310.317.834.867.929.4111.5011.62
LanzhouTop5.756.076.247.044.827.226.336.976.97
Bottom8.699.059.447.004.537.359.0810.0710.13
Floor 1El-CentroTop11.3912.5212.776.623.146.646.327.958.19
Bottom31.4236.6636.7010.25.7610.35.235.705.81
LanzhouTop10.0311.1611.555.653.005.866.516.906.91
Bottom30.2032.6633.118.895.399.225.455.936.00
Table 6. Soil elasticity modulus (MPa).
Table 6. Soil elasticity modulus (MPa).
NumberName of SoilM1M2M3
1Plain fill5.010.015.0
2New loess12.024.036.0
3Silty clay15.030.045.0
4Pebble bed55.0110.0165.0
5Silty clay15.030.045.0
Table 7. Layer drifts (mm) and layer drift angles (rad).
Table 7. Layer drifts (mm) and layer drift angles (rad).
SectionEl-Centro MotionLanzhou Motion
M1M2M3M1M2M3
Floor 3layer drift in the X-direction1.761.260.81.771.151.04
layer drift angle in the X-direction1/28981/40481/63751/28811/44351/4904
layer drift in the Y-direction4.032.371.913.141.830.76
layer drift angle in the Y-direction1/12661/20251/25131/15291/26231/6316
Floor 2layer drift in the X-direction1.651.230.861.651.120.78
layer drift angle in the X-direction1/29091/39021/55811/29091/42861/6154
layer drift in the Y-direction3.812.341.852.981.760.68
layer drift angle in the Y-direction1/12601/20511/25951/16111/27271/7059
Floor 1layer drift in the X-direction1.791.360.961.771.231.23
layer drift angle in the X-direction1/26821/35291/50001/27121/39021/3902
layer drift in the Y-direction3.832.361.872.991.750.93
layer drift angle in the Y-direction1/12531/20341/25671/16051/27431/5161
Table 8. The peak stresses and stress amplitudes under Lanzhou motion.
Table 8. The peak stresses and stress amplitudes under Lanzhou motion.
Section M1M2M3
Static (MPa)Dynamic (MPa)AmplificationDynamic (MPa)AmplificationDynamic (MPa)Amplification
−1 Floor3.974.9524.664.8020.895.0827.89
−2 Floor3.614.2216.904.1314.405.2846.23
−3 Floor3.904.4313.744.6519.385.5141.39
Side wall 1 *1.992.2714.092.105.442.2311.98
Side wall 32.933.3213.383.2611.263.5822.16
Side wall 45.436.5620.806.1112.596.5119.97
Side wall 57.728.5410.708.368.339.0417.08
Side wall 64.154.7013.144.6512.365.4130.23
Side wall A4.144.7113.884.406.484.7314.49
Side wall B4.094.141.094.100.254.223.07
Side wall C3.735.1939.084.7928.285.1437.69
Side wall D2.482.604.702.594.432.8113.19
* Side wall 1 means the number of the wall, i.e., located at the cross of axis 1. The rest have similar meanings.
Table 9. The peak internal forces of each section on surface-6 under 0.20 g.
Table 9. The peak internal forces of each section on surface-6 under 0.20 g.
SectionPeak Vertical Pressure (105 N/m)Peak Shear Force (105 N/m)Peak Bending Moment (103 N·m/m)
M1M2M3M1M2M3M1M2M3
Floor 3El-Centro motionTop4.434.284.405.425.675.818.288.288.36
Bottom5.575.325.396.176.806.857.167.677.87
Lanzhou motionTop4.614.144.505.445.475.618.018.288.61
Bottom5.304.835.276.136.266.597.177.187.36
Floor 2El-Centro motionTop7.146.696.817.508.728.947.457.777.96
Bottom10.519.809.927.929.409.799.3210.0110.33
Lanzhou motionTop6.946.196.707.227.518.206.977.628.02
Bottom9.849.209.717.357.578.449.9310.1011.06
Floor 1El-Centro motionTop12.7711.9112.116.648.008.515.396.476.75
Bottom36.7032.7932.4810.3312.2713.225.215.465.67
Lanzhou motionTop12.2511.3212.005.265.556.626.916.987.22
Bottom35.1131.2734.198.228.479.915.265.525.69
Table 10. Acceleration amplification factors.
Table 10. Acceleration amplification factors.
Acceleration Amplification FactorSoil SurfaceRoofFloor 3Floor 2Floor 1Equipment Top
El-Centro motionStructureX-direction/0.8170.820.8230.8242.600
Y-direction/0.4980.5820.6610.7381.194
SoilX-direction0.7520.8100.9851.0641.219/
Y-direction0.5940.7800.8170.8410.864/
Lanzhou motionStructureX-direction/0.8050.8280.8450.8652.554
Y-direction/0.5020.5150.5450.5991.096
SoilX-direction0.8121.2871.3271.3191.457/
Y-direction0.7880.8220.8360.8290.845/
Taft motionStructureX-direction/0. 900.9241.0471.0703.518
Y-direction/0.690.7590.8070.9541.494
SoilX-direction1.2671.3421.5511.5831.772/
Y-direction1.0211.0361.0491.0931.177/
Table 11. Layer drift (mm) and layer drift angle (rad).
Table 11. Layer drift (mm) and layer drift angle (rad).
SectionEl-Centro MotionLanzhou Motion
M1M2M1M2
Floor 3layer drift in the X-direction3.1102.7372.9492.589
layer drift angle in the X-direction1/16401/18631/17291/1970
layer drift in the Y-direction8.0746.8035.3734.611
layer drift angle in the Y-direction1/6321/7501/9491/1106
Floor 2layer drift in the X-direction2.9842.6232.8082.463
layer drift angle in the X-direction1/16091/18301/17091/1949
layer drift in the Y-direction7.6606.4545.1854.510
layer drift angle in the Y-direction1/6271/7441/9261/1064
Floor 1layer drift in the X-direction3.1352.7542.9522.610
layer drift angle in the X-direction1/15311/17431/16261/1839
layer drift in the Y-direction8.1326.7965.4024.621
layer drift angle in the Y-direction1/5901/7061/8891/1039
Table 12. Peak internal forces of the side wall.
Table 12. Peak internal forces of the side wall.
SectionShear Force (105 N/m)Bending Moment (103 N·m/m)
M1M2Increase (%)M1M2Increase (%)
Side wall 1El-Centro motion6.535.8012.5913512210.66
Lanzhou motion6.015.499.4713011711.11
Side wall 3El-Centro motion10.989.7912.1694785011.41
Lanzhou motion9.768.6812.448938129.98
Side wall 4El-Centro motion11.089.4217.6265054219.93
Lanzhou motion9.438.1016.4281868319.77
Side wall 5El-Centro motion11.7010.5710.69123411329.01
Lanzhou motion10.028.8313.485665179.48
Side wall 6El-Centro motion3.763.1618.9912610915.60
Lanzhou motion3.893.3117.52998713.79
Table 13. Layer drift (mm) and drift angle (rad) under El-Centro motion.
Table 13. Layer drift (mm) and drift angle (rad) under El-Centro motion.
SectionM1M2M3M4M5
Floor 3layer drift in the X-direction1.762.252.423.062.97
layer drift angle in the X-direction1/28981/22671/21071/16671/1717
layer drift in the Y-direction4.032.883.824.013.91
layer drift angle in the Y-direction1/12661/17701/13351/12711/1304
Floor 2layer drift in the X-direction1.652.102.232.691.81
layer drift angle in the X-direction1/29091/22861/21521/17841/2652
layer drift in the Y-direction3.812.743.603.763.62
layer drift angle in the Y-direction1/12601/17521/13331/12771/1326
Floor 1layer drift in the X-direction1.792.252.365.932.83
layer drift angle in the X-direction1/26821/21331/20341/8091/1696
layer drift in the Y-direction3.832.743.593.763.73
layer drift angle in the Y-direction1/12531/17521/13371/12771/1287
Table 14. Shear forces of internal side wall 3 in different buried depths.
Table 14. Shear forces of internal side wall 3 in different buried depths.
Top of Floor 3Bottom of Floor 3Top of Floor 2Bottom of Floor 2Top of Floor 1Bottom of Floor 1
M1Static ((105 N/m)3.534.205.598.2011.1816.88
Dynamic (105 N/m)4.575.467.2910.5814.4121.56
Amplification (%)29.5130.0230.4328.9428.8527.71
M2Static (105 N/m)6.196.427.7210.4213.6120.15
Dynamic (105 N/m)7.748.129.7713.1817.1625.41
Amplification (%)24.8826.4626.5626.4826.1126.10
M3Static (105 N/m)8.749.039.9412.7016.0623.43
Dynamic (105 N/m)11.0611.2612.6116.1720.4029.78
Amplification (%)26.6224.6826.8627.3626.9827.09
M4Static (105 N/m)11.3912.4612.4915.4018.9828.43
Dynamic (105 N/m)14.2415.0015.7419.6824.2935.85
Amplification (%)25.0620.3525.9627.7727.9926.11
M5Static (105 N/m)13.5414.9614.5317.4121.0930.13
Dynamic (105 N/m)17.2218.5018.5622.3726.9138.52
Amplification (%)27.1823.6927.7028.4627.6127.85

Share and Cite

MDPI and ACS Style

Wen, B.; Zhang, L.; Niu, D.; Zhang, M. Soil–Structure–Equipment Interaction and Influence Factors in an Underground Electrical Substation under Seismic Loads. Appl. Sci. 2017, 7, 1044. https://doi.org/10.3390/app7101044

AMA Style

Wen B, Zhang L, Niu D, Zhang M. Soil–Structure–Equipment Interaction and Influence Factors in an Underground Electrical Substation under Seismic Loads. Applied Sciences. 2017; 7(10):1044. https://doi.org/10.3390/app7101044

Chicago/Turabian Style

Wen, Bo, Lu Zhang, Ditao Niu, and Muhua Zhang. 2017. "Soil–Structure–Equipment Interaction and Influence Factors in an Underground Electrical Substation under Seismic Loads" Applied Sciences 7, no. 10: 1044. https://doi.org/10.3390/app7101044

APA Style

Wen, B., Zhang, L., Niu, D., & Zhang, M. (2017). Soil–Structure–Equipment Interaction and Influence Factors in an Underground Electrical Substation under Seismic Loads. Applied Sciences, 7(10), 1044. https://doi.org/10.3390/app7101044

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop