Shock and Vibration
Volume 2018 (2018), Article ID 9149721, 12 pages
https://doi.org/10.1155/2018/9149721

Research Article

Seismic Response Analysis of an Isolated Structure with QZS under Near-Fault Vertical Earthquakes

Dewen Liu,¹ Yang Liu,² Dongfa Sheng,¹ and Wenyuan Liao¹

¹College of Civil Engineering, Southwest Forestry University, Kunming, China
²College of Civil Engineering, Tongji University, Shanghai 200092, China

Correspondence should be addressed to Yang Liu; winner5299@163.com

Received 31 October 2017; Revised 27 January 2018; Accepted 12 February 2018; Published 15 March 2018

Academic Editor: Mahmoud Bayat

Copyright © 2018 Dewen Liu et al. 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.

Abstract.

Seismic isolation devices are usually designed to protect structures from the strong horizontal component of earthquake ground shaking. However, the effect of near-fault (NF) vertical ground motions on seismic responses of buildings has become an important consideration due to the observed building damage caused by vertical excitation. As the structure needs to maintain its load bearing capacity, using the horizontal isolation strategy in vertical seismic isolation will lead to the problem of larger static displacement. In particular, the bearings may generate large deformation responses of isolators for NF vertical ground motions. A seismic isolation system including quasi-zero stiffness (QZS) and vertical damper (VD) is used to control NF vertical earthquakes. The characteristics of vertical seismic isolated structures incorporating QZS and VD are presented. The formula for the maximum bearing capacity of QZS isolation considering the stiffness of vertical spring components is obtained by theoretical derivation. From the static analysis, it is found that the static capacity of the QZS isolation system with vertical seismic isolation components increases when the configurative parameter reduces. Seismic response analyses of the seismic isolated structure model with QZS and VD subjected to NF vertical earthquakes are conducted. The results show that seismic responses of the structure can be controlled by setting the appropriate static equilibrium position, vertical isolation period, and vertical damping ratio. Adding a damping ratio is effective in controlling the vertical large deformation of the isolator.

1. Introduction

Seismic isolation is one of most effective control ways to protect structural systems, nonstructural systems, and content from damage due to horizontal earthquake ground motions [1, 2]. However, most buildings with seismic isolation devices are not designed to reduce the vertical component of earthquake ground motions [3]. In engineering design, the vertical-to-horizontal acceleration () ratios of peak ground acceleration are usually recommended as 2/3. In addition, the shape of the vertical response spectrum is similar to that of the horizontal response spectrum. In recent years, it was found that the areas near the epicenter and faults exert a strong vertical ground motion. The vertical-to-horizontal response spectral ratios are greater than 2/3. The ratios were smaller for long periods than for short periods. In the past researches, it has been found that the response spectral ratios are strongly dependent on period and site-to-source distance during the 1994 Northridge earthquake [4].

Vertical seismic isolation of a base-isolated structure has gained the interest of researchers. Various kinds of vertical seismic isolation devices have been developed to mitigate earthquakes, such as the low shape factor (LSF) elastomeric bearings [5], air springs consisting of cylinders filled with nitrogen gas [6, 7], and vertical flexibility [8]. However, these approaches have their own drawbacks. The LSF approach is used only for low-rise structures. Vertical air springs and vertical flexibility are so costly and hence not likely to be widely applied. Furthermore, in order to achieve the seismic isolation effect, a longer vertical isolation period implies a lower vertical stiffness, which will lead to larger static displacement of the mass. This hinders the development of vertical seismic isolation systems. In recent years, negative stiffness structures (NSS) and quasi-zero stiffness (QZS) systems have been developed to overcome this disadvantage [9, 10]. Nevertheless, the new approaches’ applications are for the isolation of mechanical engineering only, and very limited researches focus on vertical seismic isolation.

Near-fault (NF) ground motions are characterized by long-period velocity and displacement pulses [11]. Past researches have revealed that seismic isolators can generate unacceptably large lateral displacement for NF motion due to resonance-like responses [12]. Accordingly, controlling the isolator’s large displacement can be a challenging problem when developing a vertical isolation system [13]. By far, considering load bearing capacity and seismic response control for the structure under NF vertical ground motions has not been addressed. Therefore, focused on NF vertical seismic protection of structures such as span systems and base-isolated structures [3, 14], the objective of this study is to present a vertical isolation system that has sufficient static vertical rigidity to prevent large static deformations due to the self-weight of supported structures, while having enough dynamic flexibility and energy dissipation to reduce NF vertical seismic responses. Also, this approach is different from the horizontal isolation strategy.

In this paper, the characteristics of static load bearing capacity dynamic response of isolated structures with QZS under NF vertical ground motions are studied. The formula of maximum bearing capacity of quasi-zero stiffness isolation incorporating vertical seismic isolation components is theoretically derived. Subsequently, the effects of the static equilibrium position, vertical period, and damping ratio on the seismic responses for the isolated system with QZS are analyzed. Finally, comparison analyses of the QZS system and the seismic isolated system without QZS are presented to check the isolation effectiveness for the model against NF vertical ground motions.

2. Seismic Isolated Structure Model with QZS

Figure 1 shows the model of a seismic isolated structure with QZS. The system consists of a SDOF isolation system, supported on an intermediate mass , together with a vertical stiffness , in parallel with a viscous damper and two oblique springs with stiffness [15, 16]. Figure 2 shows the configuration of the oblique spring, which is composed of a sleeve, coil spring, and cylindrical bar. One end of the cylindrical bar is fixed to the spring, while the other end is connected to a cover plate by a hinge joint. The sleeve is connected to a base support by a hinge joint. The cylindrical bar can be directly inserted into the hole of the sleeve. Based on the configuration, the supported cover plate is allowed to move in the vertical plane. A force is applied at the mass center of the system. In the initial location, there is an angle between the oblique spring and the horizontal line. According to the static equilibrium condition, the static force-deflection characteristic for seismic isolated structures with QZS is given by

Figure 1: Model of seismic isolated structure with QZS.

(a)

(b)

Figure 2: Model of oblique spring: (a) 3D, (b) cross section.

When it has a downward motion due to a vertical force applied to the mass of the seismic isolated structure, the vertical force can be described as

According to the nonlinear geometric relationship, is taken as follows:

If the oblique spring yields axial deflection, the relationship of and is calculated as

which can be written in a nondimensional form as

where , , , and .

3. Static Performance

Initially, the mass is held in equilibrium by the compression force of the system, which is counterbalanced by the weight of the mass. Therefore, it is necessary to meet the load bearing capacity demand for seismic isolated structures with QZS. The static property should be considered in engineering design. Define as a new parameter; substituting it into (6) yields

It is rational to obtain the maximum force of vertical capacity of a seismic isolated structure with QZS by differentiating (6) with respect to . The formula is expressed as

Setting and substituting it into (8), the following equation can be obtained:

By substituting (9) into , the nondimensional displacement is obtained as follows:

By substituting (10) into (7), the maximum nondimensional force at the position is obtained as follows:

The nondimensional maximum force of vertical capacity of a seismic isolated structure with QZS is plotted in Figure 3 for various values of and . It is found that the maximum force of vertical capacity of a seismic isolated structure with QZS is strongly dependent on , , and .

Figure 3: The nondimensional maximum force of vertical capacity of seismic isolated structures.

With a mass loaded, when the system reaches the horizontal position (, ), the dynamic stiffness is zero. At the QZS state, the system is found to satisfy [15]

where .

In practice, engineers need to choose a convenient angle of the oblique spring to the horizon in degrees to design the seismic isolated structure with QZS. Normally, some special degrees such as 10°, 30°, 45°, and 60° need be used to design the isolation bearing. When and are determined, and QZS position can be preliminarily estimated. Parameters for different angles of the QZS isolator are presented in Table 1.

Table 1: Parameters for different angles.

Degree

10
0.98
24.5
0.17

30
0.87
3.36
0.50

45
0.71
1.22
0.71

60
0.50
0.50
0.87

4. Dynamic Equation of Motion for Seismic Input

Figure 4 shows the state of the seismic isolated structure model under vertical ground motions (). When the mass oscillation reaches the steady state around the static equilibrium position, the dynamic equilibrium can be established as

which can be written in a nondimensional form as

where , , , , and .

Figure 4: Seismic isolated structure model with QZS under vertical ground motion.

A time-stepping integration procedure has been applied to solve the equation of motion (13) numerically.

5. Selection of Near-Fault Vertical Ground Motion Records

Due to the uncertainty of earthquake events, selection of ground motion records has a significant impact on the variability of the structural response. Seismic codes for seismic design of buildings have prescribed general guidelines but they do not provide specifics for selecting the type of earthquake records for nonlinear dynamic analysis. Full agreement has not yet been reached regarding the establishment of commonly accepted selection criteria for earthquake records [17].

Selection criteria for earthquake records are given as follows:

(i)These motions cover a moment magnitude range from 6.0 to 7.6.(ii)The rupture distance (closest distance from site to fault rupture plane) ranges from 0 to 10 km.(iii)Their shear-wave velocities at the uppermost 30 m () soil profile are greater than 200 m/s.

20 different earthquake records (Table 2) are used to examine the effects of the characteristics of the excitation on the seismic response of the seismically isolated building with QZS. These earthquakes are obtained from the Pacific Earthquake Engineering Research Next Generation of Ground Motion Attenuation Models (PEER-NGA) ground motion database [18]. Elastic response spectra for 5% critical damping generated with the vertical components of the recorded ground motions are presented in Figure 5. Also plotted in this figure is the average spectrum. The peak acceleration (PGA) corresponding to earthquakes of major levels is specified to be 0.4 g.

Table 2: Near-fault pulse-like record data.

Number
Event
Station
Moment
magnitude
(km)
Mechanism

Loma Prieta, USA, 1989
Saratoga-Aloha Ave., CMP
6.9
8.50
Reverse
380.89

Kocaeli, Turkey, 1999
Izmit, IZT
7.4
7.21
Strike slip
811.00

Chi-Chi, Taiwan, 1999
TCU089
7.6
8.88
Reverse
671.52

Chi-Chi, Taiwan, 1999
CHY080
7.6
2.69
Reverse
496.21

Chi-Chi, Taiwan, 1999
TCU049
7.6
3.76
Reverse
487.27

Chi-Chi, Taiwan, 1999
TCU051
7.6
7.64
Reverse
350.06

Chi-Chi, Taiwan, 1999
TCU052
7.6
0.66
Reverse
579.1

Chi-Chi, Taiwan, 1999
TCU068
7.6
0.32
Reverse
487.34

Chi-Chi, Taiwan, 1999
TCU075
7.6
0.89
Reverse
573.02

Chi-Chi, Taiwan, 1999
TCU076
7.6
2.74
Reverse
614.98

Imperial Valley, 1979
El Centro #4
6.53
7.05
Strike slip
208.91

Imperial Valley, 1979
El Centro #5
6.53
3.95
Strike slip
205.63

Imperial Valley, 1979
El Centro #7
6.53
0.56
Strike slip
201.51

Imperial Valley, 1979
El Centro #10
6.53
8.6
Strike slip
202.85

Imperial Valley, 1979
El Centro Differential Array
6.53
5.09
Strike slip
202.26

Loma Prieta, USA, 1989
Los Gatos-Lexington Dam
6.9
5.02
Reverse
1070.32

Tabas, Iran, 1978
Tabas
7.35
2.05
Reverse
767.77

Northridge, USA, 1994
LA Dam
6.7
5.92
Reverse
628.99

Northridge, USA, 1994
Pacoima Dam (upper left)
6.7
7.01
Reverse
2016.33

Kobe, Japan, 1995
KJMA
6.9
0.96
Strike slip
312

Figure 5: Elastic earthquake acceleration response spectra for 5% critical damping.

6. Seismic Response Analysis

In the previous researches, the horizontal position was a static equilibrium position. But in actual design, because of the limited deformation ability and the uncertainty static estimation, the position of quasi-zero stiffness and static equilibrium is not always the same. It is essential to study the effect of static equilibrium on the seismic isolation performance of the system. The seismic responses of a base-isolated structure with different static equilibrium positions are investigated. Numerical studies are carried out using a mathematical model of the base-isolated structure to calculate the response of interest such as the peak absolute acceleration and the peak relative displacement of the isolation structure. A simplified 1-DOF model of base isolation is shown in Figure 6. In Figure 6, the static equilibrium equation is

The dynamic equilibrium equation is

Substituting (15) into (16) yields

which can be written as

where is the height of the static equilibrium position.

Figure 6: The dynamic analysis model of a QZS isolator at a static equilibrium position.

The numerical procedure has been programmed by MATLAB software. To evaluate the performance of the system, the acceleration amplification coefficient (AAC) is an index to represent the isolation effect as follows:

where is the nondimensional acceleration of the system and is the maximum nondimensional seismic input acceleration.

Seismic responses of seismic isolated structures with QZS under vertical near-fault earthquakes are conducted to analyze the isolation effects. The comprehensive influences of static equilibrium position, vertical period, and damping ratio are investigated numerically. The length of oblique spring is  m. The angle of the oblique spring to the horizon in degrees is = 30°. We choose = 30° so that the vertical spring is not too lengthy. The initial height of the mass is  m. The resulting value of is 0.87. It is assumed that the mass reaches the position at ( = 0°) which is a QZS position. In the case of QZS position, = 3.36. The elastic displacement versus force response curves of the system for different are shown in Figure 7. It can be seen that when the height of the static equilibrium position increases, the QZS position coordinate also increases. In general, the QZS position coordinate is approximately equal to the static equilibrium position. However, the stiffness is always positive when  m. In addition, with the height of the static equilibrium position increasing, the deformation of the isolator decreases. This indicates that a larger deformation space is needed in the design. The characteristics of the displacement versus force response curves for four near-fault ground motions have similar tendency. The characteristic of the isolator changes from pure-softening characteristic to pure-hardening characteristic.

(a)

(b)

(c)

(d)

Figure 7: Comparison of elastic displacement versus force response curves for a seismic isolated structure with QZS: (a) Loma Prieta-Saratoga-Aloha Ave.; (b) Kocaeli-Izmit; (c) Chi-Chi-TCU089; (d) Chi-Chi-CHY080.

The effect of the static equilibrium position and vertical period on the nondimensional displacement and AACs of the system is shown in Figure 8. It can be seen that the nondimensional displacement of the system decreases when the height of the static equilibrium position increases. The AAC decreases as the height of the static equilibrium position decreases and vertical period increases for the general trend. When vertical periods are relatively large, the AAC is less than 1.0 in the static equilibrium position range ( m), where the isolator is able to mitigate vibrations efficiently.

(a)

(b)

Figure 8: Comparison of nondimensional displacement and AAC for a seismic isolated structure with QZS ().

Figure 9 shows the nondimensional displacement and AAC of the system with a fixed vertical period for different static equilibrium positions and damping ratios. It shows that vertical seismic responses reduce when the damping ratio decreases. Compared with the case of , the AAC is less than 1.0 for other cases. Although the increasing height of the static equilibrium position will result in a detrimental effect on seismic isolation, it is preferable to have relatively high damping to control the response of the system. It is worth noting that the isolation effect is not evident when the damping ratio changes from to . This means that it is essential to choose befitting damper parameters for a good isolator.

(a)

(b)

Figure 9: Comparison of nondimensional displacement for a seismic isolated structure with QZS ( s).

Figure 10 shows the nondimensional displacement and AAC of the system with a fixed static equilibrium position for different vertical periods and damping ratios. It can be observed that the displacement of the system increases when vertical periods increase. A significant decrease in AAC of the system with the increase in vertical periods is observed in Figure 10(b).

(a)

(b)

Figure 10: Comparison of AAC for a seismic isolated structure with QZS ( m).

7. Parameter Design and Performance Evaluation

According to previous static and seismic response results, QZS system parameters design and seismic performance evaluation are conducted. The design chart of the QZS system is shown in Figure 11. Assuming  s and rated load  kN, the vertical spring stiffness can be calculated using given parameters ( N/m). The initial design length of the oblique spring is equal to 1 m, and the initial design angle of the oblique spring to the horizon in degrees is = 30°. The parameter can be determined to be 0.5 m and . It is assumed that the mass reaches the position at  m ( = 0°) which is a QZS position. In the case of QZS position, . At the QZS position, .  kN/m. To meet the static capacity requirement, design parameters should be verified by the following formula:

Figure 11: The design procedures.

Substituting the given parameters into (11) yields = 967 kN. Obviously, the result cannot meet (11). Thus, design parameters need be adjusted to meet the static capacity requirement. Here, changes from 2 s to 1 s, and increases from 1 m to 1.2 m. Then, = 2321 kN . In this case, seismic responses of the QZS system can be obtained by numerical analysis.

In order to compare the performance between the QZS system and the seismic isolated system without QZS, we consider isolating the same mass under the same external excitation conditions. Here, rubber bearings are used for the seismic isolated system without QZS. The vertical stiffness is 1683.2 kN/mm [3]. The vertical nature period of the seismic isolated system without QZS is 0.068 s. Figure 12 shows the results of peak absolute vertical acceleration for the QZS system and the seismic isolated system without QZS. Results from the numerical simulations show that the vertical components of ground motions cause significant amplification in the acceleration of the seismic isolated system without QZS. For example, the peak acceleration response for N2 is 11.19 m/s², and the corresponding AAC is 2.80. The results presented in Figure 12 also suggest that the QZS concept can largely decrease the peak absolute vertical acceleration of the seismic isolated structure.

(a)

(b)

(c)

(d)

Figure 12: Comparison of ACCs for the QZS system and the seismic isolated system without QZS: (a) Loma Prieta-Saratoga-Aloha Ave.; (b) Kocaeli-Izmit; (c) Chi-Chi-TCU089; (d) Chi-Chi-CHY080.

The static and seismic response analysis results show that the seismic isolation system including QZS and vertical damper is a promising strategy in NF vertical ground motions. However, the limitation of the study is that if the rated load of the system exceeds a value, it is hard to realize a good isolation efficiency. Moreover, optimization selection of the complex parameters and superior ability on deformation and load capacity for materials of this system should be taken into account for the isolation system to obtain better isolation efficiency.

8. Conclusions

Vertical static and nonlinear time history analyses were performed to investigate the static capacity, acceleration, and displacement responses of the seismic isolated system with QZS. From the results presented in this paper, the following conclusions are drawn:

(1)The maximum force of the vertical capacity of a seismic isolated structure with QZS is highly related to the configurative parameter, the ratio of spring stiffness, vertical stiffness, and the length of the oblique spring. For fixed vertical stiffness, the length of the oblique spring, and the ratio of spring stiffness, the maximum force of vertical capacity of a seismic isolated structure with QZS increases when the configurative parameter reduces.(2)The effect of static equilibrium position, vertical period, and damping ratio on seismic responses for the seismic isolated system with QZS under near-fault ground motions is significant. These studies demonstrate that the seismic isolated system with QZS can provide reasonable isolation effects to control near-fault earthquakes by designing key model parameters.(3)Although the investigation suggests that the QZS concept might be a promising strategy to achieve vertical seismic protection of base-isolated buildings, further detailed analysis and bearings product are needed when the rated load of the system exceeds a value.(4)It is noted that the analysis results presented in this paper do not consider the effects of horizontal, rocking, and rotational ground motions. These effects are likely important and studies are needed to better quantify them and develop additional methods for their consideration in analysis and design.(5)The analysis model in this paper is assumed to be symmetrical, without considering the eccentricity effect. It is possible to change the critical load capacity and isolation efficiency of the analysis model under vertical and horizontal ground motions simultaneously. Particularly, if the supported structure has an excessive vertical eccentricity for mass and vertical stiffness, there will be a need to control rocking responses. In addition, overturning effects can be considered when a base-isolated structure has a large length-to-width ratio.(6)The feasibility of the seismic isolation system presented in this paper would need to be further explored including specific details of the vertical spring, damper, and oblique springs. Due to the effects of horizontal, rocking, and rotational ground motions, methods to restrain lateral deformation and rocking angle can be further studied.

Notations

: | Mass of a lumped parameter model for a seismic isolated structure
: | The length of the oblique spring
: | Horizontal projection length of the oblique spring
: | The vertical spring stiffness and damping coefficient
: | Vertical displacement of a seismic isolated structure
: | Vertical force applied to the mass of a seismic isolated structure
: | The oblique spring stiffness
: | The angle of the oblique spring to the horizon in degrees
: | The axial deflection of the oblique spring
: | The ratio of spring stiffness
: | The configurative parameter of the system
: | Vertical seismic acceleration.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

The authors gratefully acknowledge the financial support of the National Natural Science Foundation of China (no. 51508414), Yunnan Science Technology Department Fund Project (2016RA079), and Yunnan Provincial Department of Education Fund Project (2015Y298).

References

S. Nagarajaiah and S. Xiaohong, “Response of base-isolated USC hospital building in Northridge earthquake,” Journal of Structural Engineering, vol. 126, no. 10, pp. 1177–1186, 2000.

W.-G. Liu, C. Qin, Y. Liu, and Q.-R. Yang, “Shaking table tests on earthquake response characterization of a complex museum isolated structure in high intensity area,” Shock and Vibration, vol. 2016, Article ID 7974090, 23 pages, 2016.

S. Furukawa, E. Sato, Y. Shi, T. Becker, and M. Nakashima, “Full-scale shaking table test of a base-isolated medical facility subjected to vertical motions,” Earthquake Engineering & Structural Dynamics, vol. 42, no. 13, pp. 1931–1949, 2013.

Y. Bozorgnia and K. W. Campbell, “The vertical-to-horizontal response spectral ratio and tentative procedures for developing simplified V/H and vertical design spectra,” Journal of Earthquake Engineering, vol. 8, no. 2, pp. 175–207, 2004.

G. P. Warn and B. Vu, “Exploring the low shape factor concept to achieve threedimensional seismic isolation,” in Proceedings of the 20th Analysis and Computation Specialty Conference, pp. 1–11, Chicago, Ill, USA, 2012.

M. Morishita, K. Inoue, and T. Fujita, “Development of three-dimensional seismic isolation systems for fast reactor application,” Journal of Japan Association for Earthquake Engineering, vol. 4, no. 3, pp. 305–310, 2004.

T. Tomizawa, O. Takahashi, H. Aida et al., “Vibration test in a Building named “Chisuikan” using Three-dimensional Seismic Isolation System,” in proceedings of the The 15th World Conference on Earthquake Engineering, 2012.

B. Vu, M. Unal, G. P. Warn, and A. M. Memari, “A distributed flexibility and damping strategy to control vertical accelerations in base-isolated buildings,” Structural Control and Health Monitoring, vol. 21, no. 4, pp. 503–521, 2013.

T. D. Le and K. K. Ahn, “Experimental investigation of a vibration isolation system using negative stiffness structure,” International Journal of Mechanical Sciences, vol. 70, pp. 99–112, 2013.

Y. Zheng, X. Zhang, Y. Luo, B. Yan, and C. Ma, “Design and experiment of a high-static-low-dynamic stiffness isolator using a negative stiffness magnetic spring,” Journal of Sound and Vibration, vol. 360, pp. 31–52, 2016.

M. Dicleli and S. Buddaram, “Equivalent linear analysis of seismic-isolated bridges subjected to near-fault ground motions with forward rupture directivity effect,” Engineering Structures, vol. 29, no. 1, pp. 21–32, 2007.

R. S. Jangid and J. M. Kelly, “Base isolation for near-fault motions,” Earthquake Engineering & Structural Dynamics, vol. 30, no. 5, pp. 691–707, 2001.

L. Y. Lu, P. R. Chen, and K. W. Pong, “Theory and experiment of an inertia-type vertical isolation system for seismic protection of equipment,” Journal of Sound and Vibration, vol. 366, pp. 44–61, 2016.

S. K. Kunnath, E. Erduran, Y. H. Chai, and M. Yashinsky, “Effect of near-fault vertical ground motions on seismic response of highway overcrossings,” Journal of Bridge Engineering, vol. 13, no. 3, pp. 282–290, 2008.

A. Carrella, M. J. Brennan, and T. P. Waters, “Static analysis of a passive vibration isolator with quasi-zero-stiffness characteristic,” Journal of Sound and Vibration, vol. 301, no. 3–5, pp. 678–689, 2007.

C. C. Lan, S. A. Yang, and Y. S. Wu, “Design and experiment of a compact quasi-zero-stiffness isolator capable of a wide range of loads,” Journal of Sound and Vibration, vol. 333, no. 20, pp. 4843–4858, 2014.

B. Chiou, R. Darragh, N. Gregor, and W. Silva, “NGA project strong-motion database,” Earthquake Spectra, vol. 24, no. 1, pp. 23–44, 2008.

E. I. Katsanos, A. G. Sextos, and G. D. Manolis, “Selection of earthquake ground motion records: a state-of-the-art review from a structural engineering perspective,” Soil Dynamics and Earthquake Engineering, vol. 30, no. 4, pp. 157–169, 2010.