## 3 Overview of Pseudo-Dynamic Testing and Result Analysis
### Pseudo-Dynamic Testing System and Test Methods
Pseudo-dynamic testing differs from classical computer-based structural dynamic simulations, in that the latter depend on experimentally measured restoring forces rather than on restoring forces computed using a mathematical model. Figure 6 is a conceptual illustration of the pseudo-dynamic test system developed in the present study; this figure also shows how specimens are set and tested. As shown in Figure 6, it is a two-degree-of-freedom (TDF) system and composed of the numerical analysis unit, in which the input seismic ground motion is determined by the control computer and the loading test unit in which specimens are actually subjected to testing. The displacement response that has been calculated during the test is actually applied tothe specimen using two hydraulic actuators that are installed in the horizontal direction. The restoring force is experimentally measured during the test, and this measurement is fed to the control computer to calculate the corresponding displacement response. Data conversion is performed by an analog-to-digital/digital-to-analog converter ([DA-16A], [32]), and the seismic response during pseudo-dynamic testing is calculated by the closed-loop control system. Control computer-based numerical analysis was conducted using the Pseudo-Dynamic Testing Program [33]. Based on the restoring force of the specimen against the applied deformation measured with LVDT during the loading test, along with the input seismic acceleration and the measured response in the current step, the amount of response in the next step is calculated through numerical integration using the motion equation, as shown in Equation (1) below. \[M\ddot{y}(t)+C\dot{y}(t)+r(t)[=Ky(t)]=-M\ddot{y_{0}} \tag{1}\]  
where, \(M\), \(C\), and \(K\) refer to the mass, damping, and stiffness matrix of the structure, respectively. y refers to the relevant displacement vector of each layer weight for the foundation. \(r\) is the restoring force vector, while \(\ddot{y_{0}}\) is the input ground acceleration. The numerical integration of the motion equation was conducted using the \(\alpha\)-method [34], and the numerical integration algorithm for pseudo-dynamic testing is shown in Equation (2). \[Ma_{i+1}+(1+\alpha)\text{C}v_{i+1}-\alpha\text{C}v_{i}+(1+\alpha)-\alpha r_{i }=(1+\alpha)f_{i+1}-\alpha f_{i} \tag{2}\]  
\[y_{i+1}=y_{i}+\Delta tv_{i}+\Delta t^{2}\bigg{[}\bigg{(}\frac{1}{2}-\beta \bigg{)}a_{i}+\beta a_{i+1}\bigg{]} \tag{3}\]  
\[v_{i+1}=v_{i}+\Delta t[(1-\gamma)a_{i}+\gamma a_{i+1}] \tag{4}\]  
where, \(y_{i}\), \(v_{i}\), and \(a_{i}\) refer to the joint displacement at the time corresponding to \(i\Delta t\), the velocity, and the acceleration, respectively. \(\Delta t\) is the time interval for integration, \(r_{i}\) is the restoring force vector at the joint displacement, and \(f_{i}\) is the external load vector \((-M\ddot{y_{0}})\). In an elastic structure, \(r_{i}=Ky_{i}\) applies (K is the elastic stiffness matrix of the structure). \(\alpha,~{}\beta\), and \(\gamma\) are parameters that control the numerical characteristics of the algorithm. The conditions \(-5\leq\alpha\leq 0\), \(\beta=\frac{(1-\alpha)^{2}}{4}\), and \(\gamma=\frac{1}{2}-\alpha\) indicate that the system is in a state of unconditional stability. In the next step, the displacement response is calculated based on the stiffness (K), the mass (M), and the coefficient of damping proportional to stiffness (c) using Equations (1)-(4). The damping factor (\(\xi\)) was set to 0.03, which amounts to 3% of the critical damping. As shown in Figure 6, horizontal seismic response deformation is applied to the specimen using two 2000-kN hydraulic MTS actuators on the first and second floors. The horizontal displacement used to calculate the displacement response is measured using a 300-mm linear variable differential transformer (LVDT) installed on each floor.