Research Article  Open Access
Seismic Response Analysis of an Isolated Structure with QZS under NearFault Vertical Earthquakes
Abstract
Seismic isolation devices are usually designed to protect structures from the strong horizontal component of earthquake ground shaking. However, the effect of nearfault (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 quasizero 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 verticaltohorizontal 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 verticaltohorizontal 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 sitetosource distance during the 1994 Northridge earthquake [4].
Vertical seismic isolation of a baseisolated 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 lowrise 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 quasizero 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.
Nearfault (NF) ground motions are characterized by longperiod velocity and displacement pulses [11]. Past researches have revealed that seismic isolators can generate unacceptably large lateral displacement for NF motion due to resonancelike 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 baseisolated 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 selfweight 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 quasizero 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 forcedeflection characteristic for seismic isolated structures with QZS is given by
(a)
(b)
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 asAccording to the nonlinear geometric relationship, is taken as follows:
If the oblique spring yields axial deflection, the relationship of and is calculated aswhich can be written in a nondimensional form aswhere , , , 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 .
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.

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 aswhich can be written in a nondimensional form aswhere , , , , and .
A timestepping integration procedure has been applied to solve the equation of motion (13) numerically.
5. Selection of NearFault 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 shearwave 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 (PEERNGA) 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.

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 quasizero 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 baseisolated structure with different static equilibrium positions are investigated. Numerical studies are carried out using a mathematical model of the baseisolated structure to calculate the response of interest such as the peak absolute acceleration and the peak relative displacement of the isolation structure. A simplified 1DOF model of base isolation is shown in Figure 6. In Figure 6, the static equilibrium equation isThe dynamic equilibrium equation isSubstituting (15) into (16) yieldswhich can be written aswhere is the height of the 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 nearfault 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 nearfault ground motions have similar tendency. The characteristic of the isolator changes from puresoftening characteristic to purehardening characteristic.
(a)
(b)
(c)
(d)
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 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 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)
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:
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^{2}, 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)
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 nearfault ground motions is significant. These studies demonstrate that the seismic isolated system with QZS can provide reasonable isolation effects to control nearfault 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 baseisolated 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 baseisolated structure has a large lengthtowidth 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 baseisolated USC hospital building in Northridge earthquake,” Journal of Structural Engineering, vol. 126, no. 10, pp. 1177–1186, 2000. View at: Google Scholar
 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. View at: Publisher Site  Google Scholar
 S. Furukawa, E. Sato, Y. Shi, T. Becker, and M. Nakashima, “Fullscale shaking table test of a baseisolated medical facility subjected to vertical motions,” Earthquake Engineering & Structural Dynamics, vol. 42, no. 13, pp. 1931–1949, 2013. View at: Google Scholar
 Y. Bozorgnia and K. W. Campbell, “The verticaltohorizontal 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. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 M. Morishita, K. Inoue, and T. Fujita, “Development of threedimensional seismic isolation systems for fast reactor application,” Journal of Japan Association for Earthquake Engineering, vol. 4, no. 3, pp. 305–310, 2004. View at: Google Scholar
 T. Tomizawa, O. Takahashi, H. Aida et al., “Vibration test in a Building named “Chisuikan” using Threedimensional Seismic Isolation System,” in proceedings of the The 15th World Conference on Earthquake Engineering, 2012. View at: Google Scholar
 B. Vu, M. Unal, G. P. Warn, and A. M. Memari, “A distributed flexibility and damping strategy to control vertical accelerations in baseisolated buildings,” Structural Control and Health Monitoring, vol. 21, no. 4, pp. 503–521, 2013. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 Y. Zheng, X. Zhang, Y. Luo, B. Yan, and C. Ma, “Design and experiment of a highstaticlowdynamic stiffness isolator using a negative stiffness magnetic spring,” Journal of Sound and Vibration, vol. 360, pp. 31–52, 2016. View at: Publisher Site  Google Scholar
 M. Dicleli and S. Buddaram, “Equivalent linear analysis of seismicisolated bridges subjected to nearfault ground motions with forward rupture directivity effect,” Engineering Structures, vol. 29, no. 1, pp. 21–32, 2007. View at: Publisher Site  Google Scholar
 R. S. Jangid and J. M. Kelly, “Base isolation for nearfault motions,” Earthquake Engineering & Structural Dynamics, vol. 30, no. 5, pp. 691–707, 2001. View at: Publisher Site  Google Scholar
 L. Y. Lu, P. R. Chen, and K. W. Pong, “Theory and experiment of an inertiatype vertical isolation system for seismic protection of equipment,” Journal of Sound and Vibration, vol. 366, pp. 44–61, 2016. View at: Google Scholar
 S. K. Kunnath, E. Erduran, Y. H. Chai, and M. Yashinsky, “Effect of nearfault vertical ground motions on seismic response of highway overcrossings,” Journal of Bridge Engineering, vol. 13, no. 3, pp. 282–290, 2008. View at: Publisher Site  Google Scholar
 A. Carrella, M. J. Brennan, and T. P. Waters, “Static analysis of a passive vibration isolator with quasizerostiffness characteristic,” Journal of Sound and Vibration, vol. 301, no. 3–5, pp. 678–689, 2007. View at: Publisher Site  Google Scholar
 C. C. Lan, S. A. Yang, and Y. S. Wu, “Design and experiment of a compact quasizerostiffness isolator capable of a wide range of loads,” Journal of Sound and Vibration, vol. 333, no. 20, pp. 4843–4858, 2014. View at: Google Scholar
 B. Chiou, R. Darragh, N. Gregor, and W. Silva, “NGA project strongmotion database,” Earthquake Spectra, vol. 24, no. 1, pp. 23–44, 2008. View at: Publisher Site  Google Scholar
 E. I. Katsanos, A. G. Sextos, and G. D. Manolis, “Selection of earthquake ground motion records: a stateoftheart review from a structural engineering perspective,” Soil Dynamics and Earthquake Engineering, vol. 30, no. 4, pp. 157–169, 2010. View at: Publisher Site  Google Scholar
Copyright
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.