Research Article  Open Access
Chongxiang Zhang, Nan Jiang, "Effects of EquipmentStructureSoil Interaction on Seismic Response of Equipment and Structure via RealTime Dynamic Substructuring Shaking Table Testing", Shock and Vibration, vol. 2017, Article ID 1291265, 13 pages, 2017. https://doi.org/10.1155/2017/1291265
Effects of EquipmentStructureSoil Interaction on Seismic Response of Equipment and Structure via RealTime Dynamic Substructuring Shaking Table Testing
Abstract
Equation of motion for an equipmentstructuresoil (ESS) interaction system was derived using the branch substructure method. After rearrangement, this equation was applied to realtime dynamic substructuring shaking table (RTDSST) testing of the ESS system. This method adopts the equipmentstructure (ES) subsystem as the experimental substructure and a modal reduced soil model as the numerical substructure: the former is tested via the shaking table, and the latter is numerically simulated, while realtime data communication occurs between the two substructures during testing. A scale model of the ESS system was designed and underwent an RTDSST test. The experimental data were found to be consistent with the numerical calculation results, which corroborates the reliability and validity of the proposed testing method. A comparison of the experimental results from different earthquake stages implies that seismic responses of the equipment and structure decreased in general due to the intervention of soil, but the soil effect weakened with earthquake intensity.
1. Introduction
With increasing number of nonstructural members being incorporated into modern building systems, damage to precision equipment and other nonstructural members may result in huge economic losses and serious safety hazards. Recent research has shown that disregard of equipmentstructure interaction in some situations has an adverse impact on the seismic design of structures and equipment [1, 2]. Most studies of the ES subsystem were based on the rigid foundation assumption with no account taken of the soil effect, which is at odds with engineering reality and impairs the seismic performance research validity. The soilstructure interaction problem had also been investigated by using the shaking table test [3–5] and numerical analysis of soilstructure interaction [6–8]. These studies treated soil and structure as a fully coupled interaction system to elucidate the interaction effects.
Jiang et al. [9, 10] performed shaking table testing with an account of soilstructure interaction, which was found to have a strong effect on the upper structure response. The soil contribution can prolong the period of vibration of the entire model, deteriorate the upper structure actual seismic performance, and even introduce errors into the calculation results. Therefore, the overall ESS interaction system performance requires the experimental corroboration. However, in the conventional shaking table testing, construction of the entire ESS system would be too complicated and expensive. This necessitates the development of a new test method that would improve the efficiency and accuracy of such testing and provide a more reliable basis for the seismic design of structures and equipment.
The implementation of such largescale complex tests became possible with the development of the following approach. In 1992, Nakashima et al. [11] introduced the realtime dynamic substructuring (RTDS) test method. This is a method of dynamically testing a structure without experimentally testing a physical model of the entire system. Instead, the structure can be split into two coupled parts, the region of particular interest, which is tested experimentally, and the remainder, which is tested numerically. The force and displacement at the interface of numerical and physical substructures are transferred via a transfer system using a servohydraulic actuator or a shaking table.
This method can effectively reduce the model scale and testing difficulties. In the RTDSST tests with an account of the soil effect [12–14], the soil and upper structure are mostly treated as the numerical and experimental substructures, respectively, while the dynamic response of the overall coupled system investigation is based on the realtime data communication between the substructures. Given the realtime requirement of testing, calculations of the numerical substructure are often reformed via simplified models, such as the lumped parameter soil model proposed by Wang et al. [12], the modified Penzien model adopted by Yan et al. [13], and the simplified singledegreeoffreedom (DOF) model introduced by Tang et al. [14]. The use of simplified computational models can ensure efficient calculation of numerical substructure in substructure testing but may hurt the accuracy of calculation when applied to complex soil models. Moreover, the direct application of complex finite element models does not allow one to implement realtime substructure testing. Jiang et al. [15, 16] employed a branch substructure method in an indepth study of soilstructure interaction, wherein the calculation time was substantially reduced by the soil modal reduction. The motion equation of equipstructuresoil dynamic interaction system can be deduced using the method of branch substructure. After rearrangement, this equation was applied to realtime dynamic substructuring shaking table test.
In this study, the ESS interaction system equation of motion was derived using the branch substructure method presented elsewhere in the literature [15]. The equation was then transformed and applied to the RTDSST testing of an ESS system. This testing method adopts the ES subsystem as the experimental substructure and the soil as the numerical substructure. The former is subjected to loading and control imposed by the shaking table, while the latter is numerically simulated using a respective software program. A scale model of the ESS interaction system was designed, where the equipment was reduced to a singleDOF model and the structure was modeled as a fourstory steel frame. An experimental system was constructed to perform an RTDSST test on the model, to study the behavior of the ESS interaction system under seismic excitations. After the reliability of the testing method proposed in this study was verified, a comparative analysis of the seismic performance of the equipment and structure under the soil effect was performed.
2. The RTDSST Testing Method Based on the Branch Substructure Method
An ESS interaction system consists of three components: equipment, structure, and soil. The testing method proposed in this study treats the ES subsystem as the experimental substructure and the soil as the numerical substructure. The former is subjected to the loading and control applied by the shaking table, while the latter is computed and analyzed by the simulation software. Realtime data communication occurs between the two substructures during testing. Based on the branch substructure method, the ESS interaction system is subdivided into branch corresponding to the rigid ES system on elastic soil and branch corresponding to the deformed ES system on rigid soil, as is shown in Figures 1(b) and 1(c), respectively.
The modal transformation matrix consisting of th order modes for branch , denoted by , can be obtained as follows: where and are stiffness and mass matrices of branch , respectively. The calculation matrices after modal transformation of branch are as follows:where and are physical displacements and modal coordinates of branch , respectively; denotes the damping matrix of branch , which can be calculated using the general damping theory; and is the load matrix of branch . As the ES subsystem has a small number of DOFs, there is no need for modal reduction, and the complete calculation matrices are retained. The displacement of the ES subsystem, , is composed of two parts, rigidbody displacement of the ES subsystem (Figure 1(b)) caused by soil displacement and the intrinsic displacement of the subsystem (Figure 1(c)), and can be written in the following form: where can be calculated from the rigidbody displacement of the ES subsystem induced by the deformation of soil in branch . The relationship between the displacements of the ES subsystem and soil can be expressed in the matrix form:
Combining and rearranging the above equations by the coordinate relationship of deformation gives the equation of motion for the entire ESS interaction system:where , , , and are the mass, damping, stiffness, and load matrices of the ES subsystem in branch , respectively. The coupling term in the equation above acts as a link between the ES subsystem and soil. By rearranging the couplingrelated terms to the righthand side of this equation, the equations of motion of the soil and the ES subsystem are obtained. These equations are as follows:
In the shaking table testing substructure used in this study, the equipment in the ESS system was simulated with a singleDOF model, and the fourstory steel frame was reduced to a shear model with four DOFs. Thus, the equipment and structure formed a computational model with five DOFs, so can be presented as a matrix ( is the total number of DOFs of the soil) as follows:As seen from Figure 2, which depicts the external and coupling loads on the righthand side of the equation of motion for the ES subsystem in (7), the effect of foundation rotational acceleration, , on each DOF correlates with height. To realize the foundation rotational acceleration in the horizontal shaking table test, is converted into an equivalent horizontal load based on the principle of equivalence. The respective conversion is given by
The procedure of the RTSST testing of the ESS interaction system is as follows:(a)Assume that the external load and coupling load acting on the soil at the th time step are known.(b)Calculate the soil acceleration response at the ()th time step.(c)Obtain the external and coupling loads acting on the ES subsystem.(d)Calculate the equivalent horizontal acceleration based on the principle of equivalence and use the results as a new input data to control the shaking table and, thus, apply the equivalent load to the experimental substructure.(e)Calculate the force generated by the experimental substructure from the measured acceleration of the ES subsystem and seismic excitations and then transfer this force to the soil for use in the next numerical simulation. In this way, the interactive communication of experimental and numerical data between the shaking table and simulation software is realized until the test termination.
The stiffness matrix of equipmentstructuresoil motion equation will be decoupled using the branch substructure method so that the motion equation can be directly applied to a shaking table realtime substructure experiment. The branch substructure method also provides a suitable solution for other complex interaction systems.
3. Simulation and Test Procedures
A onefifth scale model of an ESS interaction system was elaborated in this study. In this model, the equipment was a singleDOF model built on the top of the fourstory steel frame structure. The ES subsystem was used as the experimental substructure, which would be subjected to loading and control imposed by the shaking table during testing. The soil was modeled as a type III site, and it was adopted as a numerical substructure in the RTDSST testing. Because the finite element model of soil which contained too many DOFs may not satisfy the realtime requirement of shaking table substructure experiment. The Ritz vector method was used to reduce DOFs. The testing was performed on a 3 m × 3 m shaking table at the Beijing University of Technology. The El Centro, TianJin, and PerSon ground motions, which represent different loading patterns applicable for type III sites, were selected as the seismic inputs and the respective acceleration time history curves (Figure 3). The El Centro and TianJin ground motions were seismic records waves, while the PerSon ground motion was an artificial seismic wave. The strain, acceleration, and displacement responses of the total system were measured during testing.
(a) El Centro
(b) TianJin
(c) PerSon
3.1. EquipmentStructure Model
The structure of the ES subsystem model in Figure 4 consisted of Hbeam steel sections. The first story of the structure was 0.68 m high, and each of the other stories was 0.63 m high. The structure length and width were both equal to 1.6 m. Each of the lower three stories weighed 1,700 kg, while the fourth story weighed 1,540 kg. The equipment in the model consisted of circular steel pipes and had a height of 0.5 m and a weight of 90 kg. The elastic modulus and yield strength of the Hbeam steel sections were 202.0 GPa and 339.6 MPa, respectively, while those of the circular steel pipes were 192.0 GPa and 421.4 MPa, respectively. A steel plate was attached to the center of the structure top, while the equipment was anchored to the plate using anchor bolts during testing.
3.2. Soil Model
Table 1 shows the material parameters of the soil model. The overall dimensions of the soil model were 30 m × 15 m × 15 m, whereas the foundation embedded in it had dimensions of 2.2 m × 2.2 m × 0.4 m. A soil finite element model (FEM) (Figure 5) was developed using the ANSYS commercial finite element software by the experimental design requirements [17]. An eightnode solid element was used to model the soil, with each node having three DOFs, that is, translation along the , , and directions. The element discretization should satisfy the relationship , where is the wavelength, which ensures that the element width does not exceed fivefold values of the element height. The motion of the experimental substructure model is dominated by translation along the axis, and the soil primarily underwent translational displacement along the axis. Therefore, a constraint equation of the FEM was applied for the problem simplification. Regarding boundary conditions, the soil bottom was fixed, while its other sides were modeled by viscoelastic boundaries [18].

As the soil FEM had a large number of DOFs, it could not be directly used in the RTDSST testing. Therefore, modal reduction of the soil was performed. First, calculation matrices for the soil were exported via the ANSYS software command stream and then imported into MATLAB for the modal reduction. Given the discrete distribution of the dominant soil mode, Li and Jiang [19] analyzed the relationship between the external load on soil and its vibration behavior and calculated the soil modal vector using the Ritz vector method. It was found that this method could tune off the vibration modes that had a relatively feeble effect on the soil vibration. In the current study, this method was employed for the realization of the Ritz vector analysis with the soil vibration mode and frequency assessment at each order. Table 2 compares the calculation results with those from the complete modal model of soil. The standard requires the total ratio of participating mass to exceed 90%. The mass participation factors for various orders of the Ritz vectors were calculated for the soil. According to the results obtained, the sum of the mass participation factors for the 1st through 60th orders reached 0.95. The modal transformation matrix for the soil model, , was determined by the above method. Then the calculation matrices for the reduced soil FEM were determined using (2).

The equipmentstructurenonlinear and equipmentstructurereduced soil models were elaborated, and the calculation results from the two models were compared to evaluate the modal reduction feasibility. The translational acceleration at the center of the foundation was used as the evaluation criterion. The ES subsystem can be reduced to a 5DOF computational model, and the corresponding calculation matrices can be derived from the design parameters of the subsystem [20]. The design parameters imply that the each of the lower three stories of the structure weighed 1,700 kg, that is, 1,700 kg, and the mass of fourth story, , was 1,540 kg. Based on this, the stiffness () and damping coefficients () of these stories were calculated: the 1st story: 2,423,080 N/m, 4,463 N/(m/s); the 2–4th stories: 3,833,120 N/m, 7,061 N/(m/s). The mass, stiffness, and damping coefficients of the equipment are as follows: kg, 93,046 N/m, and N/(m/s). The design parameters of the nonlinear and reduced soil models are presented in Tables 1 and 2, respectively. Based on the design parameters, an equipmentstructurenonlinear soil model was constructed via the ANSYS, and an equipmentstructurereduced one was created in MATLAB using (5). Then the equation of motion was solved via the Newmarkβ direct integration method.
Taking the TianJin ground motion during the minor, medium, and major earthquakes stages as an example, the acceleration responses at the foundation center were calculated by the two computational models and depicted in Figure 6. The foundation center peak acceleration during the minor, medium, and major earthquakes calculated by the reduced soil model exceeded those predicted by the nonlinear one by 2.4%, 3.7%, and 11.5%, respectively. The acceleration evolutions derived by the two models were consistent for the minor and medium earthquakes but differed dramatically in case of a major earthquake. This was primarily because the soil behavior was linear or exhibited a small nonlinear deformation during the minor and medium earthquakes and thus could be quite accurately described by the reduced model, while during a major earthquake the soil underwent significant nonlinear deformation, which was not taken into account by the reduced model. This proves that the reduced soil model is not applicable to the case of a major earthquake. Therefore, the study adopted the reduced soil model for the minor and medium (up to 8th degree) earthquakes in the RTDSST testing, in order to ensure the testing feasibility.
(a) Minor earthquake
(b) Medium earthquake
(c) Major earthquake
In this way, the modal reduction can effectively reduce the number of DOFs of the soil model, while ensuring the accuracy of calculations, so that adopting the reduced soil model as the numerical substructure in the analysis will provide a more efficient calculation of the numerical substructure.
3.3. Test Procedure
Figure 7 shows the composition of the RTDSST test system. The test procedure involves three major steps: () constructing a realtime test system; () designing substructure controllers; and () calculating the numerical substructure. A specialized computer for engineering control was used as the server in the calculation of the numerical substructure, whereas the SIMULINK simulation software was implemented for use in the substructure controller design and calculation of the numerical substructure. With the SIMULINK Desktop RealTime [21], the test module can be compiled into a C code that is executable in a realtime environment. Then a computer equipped with a compiler and an I/O device can act as a realtime system and communicate with the external equipment via the I/O device. In this study, a closedloop realtime substructure test system was set up using this method, which offered a costeffective data communication during the realtime substructure testing process.
The dynamic response of the numerical substructure was calculated using the SIMULINK state equations. The calculation modules are shown in Figure 7. The state variable is composed of the modal displacement and modal velocity of the soil. The input variable is composed of the external and coupling loads acting on the soil. The output variable is the physical acceleration of the soil after coordinate transformations. Then, the translational and rotational acceleration of the foundation were selected and incorporated via the selector module, and its equivalent horizontal acceleration was calculated using (9). The actual seismic excitations experienced by the experimental substructure were then determined from the acceleration and used as the driving commands for the shaking table. The matrices in the state equations, denoted by , , , and , are given in
To track the testing setup performance more accurately, the compensation control was implemented in the RTDSST test system. Timedelay compensation normally involves the loading equipment reduction to a timedelay system, giving a constant value to the system time delay , and sending output test commands seconds in advance by different prediction approaches. This procedure can ensure a consistency of the loading system input and output signals. Various prediction methods, such as singlestep prediction [22, 23] and multistep prediction [24], have been proposed and applied to the realtime substructure testing. Insofar as a shaking table has a complex dynamic behavior, its reduction to a timedelay system will not completely compensate the errors in the loading system. In the inverse dynamics compensation method of the shaking table [25], the relationship between the input signal and output signal can be expressed as . If an accurate mathematical model of the shaking table is constructed and an inverse solution is derived from it, then , and the amplitude and phase errors in the loading system would be eliminated. However, the complex dynamic behavior of the testing setup, combined with the influence of some external factors, often hinders the achievement of desired compensation effect by this method.
Tagawa et al. introduced a virtual object control concept and proposed an inverse dynamics compensation via the simulation scheme [26]. Further studies [27, 28] have corroborated that this scheme application allowed one to avoid interference from noise and other interfering signals from the external sources and verified its effectiveness through the practical application to nonlinear control. Guo et al. [29] further explored a strategy for inverse dynamics compensation and developed a new control methodology called fullstate compensation via simulation (FSCS). Given various uncertainties, a control method combining the closedloop control and realtime inverse dynamics compensation via simulation was applied to the shaking table control in this study. Figure 8 illustrates the principle of this method, which involves separate compensations of the shaking table predictable and unpredictable behaviors. For dynamic behavior that could be described by mathematical models, the inverse dynamics compensation via simulation was implemented. For uncertain dynamics that could not be described by mathematical models, a physical closedloop control was adopted.
Based on the measured acceleration of shaking table and the reference signal within control computer, the dynamic characteristics of the shaking table system were calculated, while the fourthorder transfer function of the shaking table was identified using the MATLAB function “fminsearch.” It can be expressed as follows:
The amplitudes and phases obtained from the shaking table testing and the systematic identification results are plotted in Figure 9 to describe the shaking table dynamic behavior, where the transfer function fitted by the systematic identification is capable of describing the actual vibrational motion of the shaking table and is thus applicable to the inverse dynamics compensation of the numerically controlled part of the test system. The physical part is driven by the shaking table controller, which can further reduce the experimental error and provide more accurate tracking and control of the shaking table trajectory.
4. Results
The RTDSST of the ESS interaction system was carried out using the constructed testing platform established in Section 3.3. Ground motions with degree VIII intensity were input to observe the system responses during minor and medium earthquakes. The equipment top acceleration was used as an indicator to assess the reliability of the proposed RTDSST testing method. The results from the RTDSST testing were compared with those calculated by the integrated numerical calculation. To evaluate the soil effect on seismic responses of the equipment and structure, the shaking table tests on the ESS system and ES subsystem were performed, whereas ground motions were directly applied to the latter.
4.1. Experimental Verification
The validity of the proposed RTDSST testing method was verified by a comparison of the experimental and numerical acceleration of the equipment top. The numerical acceleration was calculated by the equipmentstructurereduced soil model implemented in MATLAB. The comparison between the experimental acceleration from the RTDSST testing and the numerical results from the integrated numerical calculation is presented in Figure 10. The comparative analysis between experimental and the numerical accelleration responses of the equipment top reveals that, during a minor earthquake with degree VIII intensity, the experimental peak acceleration of the equipment top in response to the El Centro, TianJin, and PerSon ground motions was decreased by 6.0%, 4.5%, and 3.8%, respectively, as compared to the corresponding numerical values. The equipment top acceleration evolutions obtained by the two methods exhibited a close correlation for the three ground motions. The satisfactory overall accuracy of the testing proves that the proposed RTDSST testing method is reliable and valid. The deviations between the experimental and numerical results can be attributed to the limited performance of the shaking table controller and the differences between the computational model and the actual model of the ES subsystem.
(a) El Centro
(b) TianJin
(c) PerSon
4.2. Variation Characteristics
Before the earthquakesimulating excitation test, a sine sweep test was conducted separately for the equipment and structure fixed on the shaking table, to investigate their vibration characteristics. The sine sweep signals with frequencies from 0 Hz to 40 Hz were used to excite the shaking table. Then acceleration responses were measured at different locations of the equipment and structure. Next, the fast Fourier transform (FFT) was applied to the acceleration signals. The fundamental frequencies of the equipment and structure were found to be 5.12 Hz and 2.38 Hz, respectively. After that, the equipment was connected to the structure to form an ES subsystem, which was then mounted on the shaking table. The results of sine sweep for the equipment top indicated that the initial first and secondorder frequencies of the subsystem were 2.35 Hz and 4.97 Hz, respectively. The former was close to the basic structure frequency, and the latter approached the basic equipment frequency.
After the minor and medium earthquake stages of the RTDSST test were completed, sine sweep was performed for the ES subsystem. The vibration characteristics of the ES subsystem were obtained through the FTT of the acceleration signals from the equipment top and plotted in Figure 11. As seen from Figure 11, after the minor earthquake stage was finished, the first and secondorder frequencies of the subsystem did not differ from its natural ones. When the medium earthquakesimulating stage was over, the subsystem first and secondorder frequencies were 2.28 Hz and 4.69 Hz, respectively, which indicated their reduction by 3.0% and 5.6%, respectively, from their corresponding initial values. These results strongly suggest that the ES subsystem entered into the nonlinear state in the medium earthquake stage but fail to specify in which particular cases the plastic deformation of the equipment or structure took place.
(a) Minor earthquake
(b) Medium earthquake
4.3. Strain Response
Strain gauges were mounted to the equipment and structure bottom parts to measure their strain responses. The variations in peak strain can reflect the occurrence of deformation of the equipment or structure. The strain measurements of the equipment and structure indicated that the equipment behavior was linear for the El Centro and PerSon ground motions at minor and medium earthquake intensities, as well as for the TianJin ground motion at the minor earthquake intensity. During the medium earthquake stage of the TianJin ground motion, the equipment bottom in the ESS system and the ES subsystem exhibited peak strain values of 2,204 με and 2,392 με, respectively, while the maximal peak strain of the structure bottom in the ESS system and the ES subsystem reached 1,856 με and 2,081 με, respectively, which exceeded the elastic limits of the equipment and structure. This implies that the hardware and structure both exhibited a nonlinear behavior at this stage. The strain measurements also revealed that the states of the equipment and structure varied with ground motion types. During the medium earthquake, both the equipment and structure behaviors were nonlinear for the TianJin ground motion and linear for other cases.
4.4. Acceleration Response
Variations in the acceleration responses of the equipment and structure under different ground motions were analyzed to assess the soil effect on their seismic performance. Tables 3 and 4 summarize the peak acceleration measured at the tops of the equipment and structure in the ES subsystem and the ESS system, respectively. A comparative analysis shows that, due to the soil effect, the equipment top peak acceleration under the El Centro, TianJin, and PerSon ground motions decreased by 9.2%, 18.3%, and 12.0%, respectively, during the minor earthquake stage. Meanwhile, the structure top peak acceleration dropped by 15.1%, 23.6%, and 10.9%, respectively. During the medium earthquake stage, the interaction between the ES subsystem and the soil resulted in the reduction of equipment top peak acceleration under the El Centro, TianJin, and PerSon ground motions by 9.4%, 3.7%, and 11.5%, respectively. The corresponding reductions of the structure top acceleration were 14.5%, 16.8%, and 11.4%, respectively.


The comparative analysis about the acceleration responses of the equipment and structure proves that, in the presence of the soil, the ESS subsystem experienced a reduction of the equipment top and structure top acceleration responses, and the rate of this reduction closely correlated with the type of ground motion and seismic intensity. For instance, under the TianJin ground motion, the equipment top acceleration decreased by 18.3% during the minor earthquake, which was much higher than the reduction by 3.7%, which was observed during the medium earthquake simulation.
4.5. Shear Forces
The shear forces of the ES subsystem under different conditions can be calculated from the acceleration responses. Figure 12 compares the shear force responses of the ES subsystem and the ESS system during the minor earthquake: the ESS system shear forces are much lower than those of the ES subsystem. In case of the soil effect account, the peak values of the ESS system shear force under the El Centro, TianJin, and PerSon ground motions decreased by 20.3%, 17.7%, and 28.7%, respectively. The shear force responses of the ES subsystem and the ESS system during the medium earthquake are illustrated in Figure 13. Under the El Centro, TianJin, and PerSon ground motions, the shear force peak values of the ESS system dropped by 20.5%, 4.3%, and 28.2%, respectively, as compared to the ES subsystem.
The comparison of shear force responses between ES subsystem and ESS system shows that the shear forces of the ESS system tended to decrease under the soil effect because the interaction with the soil mitigated the seismic responses of the equipment and structure and thereby reduced the shear forces of the subsystem. However, under the TianJin ground motion, the rate of shear force reduction varied significantly between the two earthquake stages. This can be attributed to the fact that the ES subsystem behavior became nonlinear at the medium earthquake stage.
4.6. Displacement Response
The displacement responses of all stories of the equipment and structure were measured. Figure 14 shows the peak displacements of different stories of the ES subsystem and the ESS system during the minor earthquake. The measurements from the first through fourth stories represent the peak displacements of the four stories of the structure, and those from the fifth story are the equipment top peak displacements. One can see from Figure 14 that the peak displacements of each story of the ESS system are much smaller than the corresponding values of the ES subsystem. Under the El Centro, TianJin, and PerSon ground motions, the peak displacements of the equipment top in the ESS system decreased by 19.1%, 26.5%, and 19.0%, respectively, as compared to those measured in the ES system. In case of the soil effect account, the respective structure top peak displacements exhibited the reduction by 16.2%, 20.1%, and 25.5%. The peak displacement responses during the medium earthquake are shown in Figure 15. In case of the soil effect account, the peak displacements of the equipment top under the three types of ground motions dropped by 18.6%, 8.1%, and 18.9%, respectively, while those of the structure top decreased by 15.9%, 9.5%, and 26.3%, respectively.
As described in this section about the displacement responses of the equipment and structure, the peak displacements of the equipment and structure were substantially reduced due to the soil effect, mainly because the soil damping increased the energy dissipation. During a minor earthquake, the maximum drop of the equipment peak displacement occurred under the TianJin ground motion, while the maximum decrease in the structure peak displacement occurred under the PerSon ground motion. Under the TianJin ground motion, the tops of the two components underwent greater variations in peak displacement during a medium earthquake than during a minor one, primarily due to nonlinear behavior of both the equipment and structure during a medium earthquake. This implies that, as the seismic intensity increased, the selfdeformation states of the equipment or structure tended to exert a significant influence on their seismic responses.
In this study, the integrated testing of an ESS interaction system was successfully implemented using a realtime shaking table testing method based on the branch substructure method. The proposed method adopts the equipmentstructure subsystem as the experimental substructure and effectively combines the numerical calculation with experimental testing. Moreover, the modal reduction of the numerical model of soil improved the efficiency of the realtime substructure testing with no deterioration of the calculation accuracy. Therefore, the proposed method is quite precise and applicable to practical experiments, which makes it quite lucrative and instrumental in the evaluation of seismic performance of various interactive systems.
5. Conclusions
In this study, the branch substructure method was used to derive the equation of motion of an equipmentstructuresoil interaction system, which was then subjected to the RTDSST test. The results obtained make it possible to draw the following conclusions:(1)After rearrangement, the formulated equation of motion was successfully applied to the RTDSST testing of the ESS interaction system. This testing method has a clear concept and is convenient to use.(2)The foundation rotational acceleration was converted to the equivalent horizontal acceleration applied to the experimental substructure based on the principle of equivalence. In this way, the rotation of foundation was taken into account in the translational shaking table testing and, thus, the accuracy of test results was improved.(3)A comparison of the experimental and numerical results revealed their good fit, which confirms the reliability and validity of the RTDSST testing method.(4)The soil effect on the equipment and structure is manifested by the overall reduction of their seismic responses. With an increase in the seismic intensity, the selfdeformation states of the equipment and structure increasingly affect their seismic responses.
The modal reduction involved in the RTDSST method, which is proposed in this study, quite accurately accounts for the effect of soil with a linear behavior, but further research is needed to investigate more complex cases of nonlinear soil behavior.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This work was financially supported by the National Science Foundation Project (Research Project no. 51478312), National Youth Science Foundation Project (Research Project no. 51208356), and the National Science Foundation Project (Research Project no. 51278335), China.
References
 J. Li, H. Chen, and Z. S. Sun, “Shaking table tests on spatial structureequipment model systems,” Engineering Mechanics, vol. 20, no. 1, pp. 157–161, 2003. View at: Google Scholar
 E. Pantoli, X. Wang, M. Chen, T. Hutchinson, B. Meacham, and H. J. Park, “Shake table testing of a fullscale fivestory building: performance of the major nonstructural componentsegress and facades,” in Proceedings of the Structures Congress 2013: Bridging Your Passion with Your Profession, pp. 1447–1459, Pittsburgh, Pennsylvania, PA, USA, May 2013. View at: Google Scholar
 G. Biondi, M. R. Massimino, and M. Maugeri, “Experimental study in the shaking table of the input motion characteristics in the dynamic SSI of a SDOF model,” Bulletin of Earthquake Engineering, vol. 13, no. 6, pp. 1835–1869, 2015. View at: Publisher Site  Google Scholar
 X. Gao, X.Z. Ling, L. Tang, and P.J. Xu, “Soilpilebridge structure interaction in liquefying ground using shake table testing,” Soil Dynamics and Earthquake Engineering, vol. 31, no. 7, pp. 1009–1017, 2011. View at: Publisher Site  Google Scholar
 R. Paolucci, M. Shirato, and M. T. Yilmaz, “Seismic behaviour of shallow foundations: Shaking table experiments vs numerical modelling,” Earthquake Engineering & Structural Dynamics, vol. 37, no. 4, pp. 577–595, 2008. View at: Publisher Site  Google Scholar
 G. Abate, M. R. Massimino, and M. Maugeri, “Numerical modelling of centrifuge tests on tunnel–soil systems,” Bulletin of Earthquake Engineering, vol. 13, no. 7, pp. 1–25, 2015. View at: Publisher Site  Google Scholar
 R. Figini, R. Paolucci, and C. Chatzigogos, “A macroelement model for nonlinear soilshallow foundationstructure interaction under seismic loads: Theoretical development and experimental validation on large scale tests,” Earthquake Engineering & Structural Dynamics, vol. 41, no. 3, pp. 475–493, 2012. View at: Publisher Site  Google Scholar
 P. Galvín and A. Romero, “A MATLAB toolbox for soilstructure interaction analysis with finite and boundary elements,” Soil Dynamics and Earthquake Engineering, vol. 57, pp. 10–14, 2014. View at: Publisher Site  Google Scholar
 X. L. Jiang and B. W. Xu, “Analysis of shaking table test of largescale soilpilecomplex structure interaction,” China Civil Engineering Journal, vol. 43, no. 10, pp. 98–105, 2010. View at: Google Scholar
 X. Jiang and Y. Li, “Shaking table test of soilpileeccentric structure interaction system,” Journal of Building Structures, vol. 31, no. 8, pp. 106–111, 2010. View at: Google Scholar
 M. Nakashima, H. Kato, and E. Takaoka, “Development of real‐time pseudo dynamic testing,” Earthquake Engineering & Structural Dynamics, vol. 21, no. 1, pp. 79–92, 1992. View at: Publisher Site  Google Scholar
 Q. Wang, J.T. Wang, F. Jin, F.D. Chi, and C.H. Zhang, “Realtime dynamic hybrid testing for soilstructure interaction analysis,” Soil Dynamics and Earthquake Engineering, vol. 31, no. 12, pp. 1690–1702, 2011. View at: Publisher Site  Google Scholar
 X.Y. Yan, Z.X. Li, Q. Han, and X.L. Du, “Shaking tables test on a longspan rigidframed bridge considering soilstructure interaction,” Engineering Mechanics, vol. 31, no. 2, pp. 58–65, 2014. View at: Publisher Site  Google Scholar
 Z. Tang, H. Ma, J. Guo, and Z. Li, “Effect of SoilStructure Interaction on Seismic Performance of LongSpan Bridge Tested by Dynamic Substructuring Method,” Shock and Vibration, vol. 2017, Article ID 4358081, 12 pages, 2017. View at: Publisher Site  Google Scholar
 X. L. Jiang and S. C. Yan, “Two step method of the branched modal and its use for liquidstructurepilesoil interaction system,” Journal of Vibration Engineering, vol. 7, no. 4, pp. 346–351, 1994. View at: Google Scholar
 F. Wang and X. L. Jiang, “Analysis of soilstructure interaction system based on mixed branch mode and constrained mode twostep method,” Journal of Earthquake Engineering and Engineering Vibration, vol. 30, no. 4, pp. 24–30, 2010. View at: Google Scholar
 F. Wang, Dynamic substructure method for soilstructuremultistory buildings interaction system, Tianjin University, 2010.
 X. L. Jiang and H. S. Zhang, “Mixed constraint modal method for analysis of soilstructure system under complex conditions,” Journal of Hohai university (Natural sciences), vol. 44, no. 4, pp. 283–290, 2016. View at: Google Scholar
 Y. Li and X. Jiang, “Parametric analysis of eccentric structuresoil interaction system based on branch mode decoupling method,” Soil Dynamics and Earthquake Engineering, vol. 48, pp. 63–70, 2013. View at: Publisher Site  Google Scholar
 S.K. Lee, E. C. Park, K.W. Min, and J.H. Park, “Realtime substructuring technique for the shaking table test of upper substructures,” Engineering Structures, vol. 29, no. 9, pp. 2219–2232, 2007. View at: Publisher Site  Google Scholar
 R. Grepl, “Realtime control prototyping in MATLAB/simulink: review of tools for research and education in mechatronics,” in Proceedings of the IEEE International Conference on Mechatronics (ICM 2011), pp. 881–886, Turkey, April 2011. View at: Publisher Site  Google Scholar
 T. Horiuchi, M. Inoue, T. Konno, and Y. Namita, “Realtime hybrid experimental system with actuator delay compensation and its application to a piping system with energy absorber,” Earthquake Engineering & Structural Dynamics, vol. 28, no. 10, pp. 1121–1141, 1999. View at: Publisher Site  Google Scholar
 F.D. Chi, J.T. Wang, Q. Wang, and F. Jin, “Delaydependent stability analysis of MDOF realtime dynamic hybrid testing considering compensation,” Engineering Mechanics, vol. 28, no. 4, pp. 200–207, 2011. View at: Google Scholar
 M. I. Wallace, D. J. Wagg, and S. A. Neild, “An adaptive polynomial based forward prediction algorithm for multiactuator realtime dynamic substructuring,” Proceedings of the Royal Society A Mathematical, Physical and Engineering Sciences, vol. 461, no. 2064, pp. 3807–3826, 2005. View at: Publisher Site  Google Scholar  MathSciNet
 J. E. Carrion, B. F. Spencer Jr., and B. M. Phillips, “Realtime hybrid simulation for structural control performance assessment,” Earthquake Engineering and Engineering Vibration, vol. 8, no. 4, pp. 481–492, 2010. View at: Publisher Site  Google Scholar
 Y. Tagawa, J.Y. Tu, and D. P. Stoten, “Inverse dynamics compensation via 'simulation of feedback control systems',” Proceedings of the Institution of Mechanical Engineers, Part I: Journal of Systems and Control Engineering, vol. 225, no. 1, pp. 137–153, 2011. View at: Publisher Site  Google Scholar
 X. Ji, K. Kajiwara, T. Nagae, R. Enokida, and M. Nakashima, “A substructure shaking table test for reproduction of earthquake responses of highrise buildings,” Earthquake Engineering & Structural Dynamics, vol. 38, no. 12, pp. 1381–1399, 2009. View at: Publisher Site  Google Scholar
 J. Y. Tu, P. Y. Lin, D. P. Stoten, and G. Li, “Testing of dynamically substructured, baseisolated systems using adaptive control techniques,” Earthquake Engineering & Structural Dynamics, vol. 39, no. 6, pp. 661–681, 2010. View at: Publisher Site  Google Scholar
 J. Guo, Z. Tang, S. Chen, and Z. Li, “Control strategy for the substructuring testing systems to simulate soilstructure interaction,” Smart Structures and Systems, vol. 18, no. 6, pp. 1169–1188, 2016. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2017 Chongxiang Zhang and Nan Jiang. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.