Seismic excitation, which results in large horizontal relative displacements, may cause collisions between two adjacent structures due to insufficient separation distance between them. Such collisions, known as earthquake-induced structural pounding, may induce severe damage. In this paper, the case of pounding between two adjacent buildings is studied by the application of single degree-of-freedom structural models. Impact is numerically simulated with the use of a nonlinear viscoelastic model. Special attention is focused on calculating values of impact forces during collisions which have significant influence of pounding-involved response under ground motions. The results of the study indicate that the impact force time history is much dependent on the earthquake excitation analyzed. Moreover, the peak impact forces during collision depend substantially on such parameters as gap size, coefficient of restitution, impact velocity, and stiffness of impact spring element. The nonlinear viscoelastic model of impact force with the considered relation between the damping coefficient and the coefficient of restitution has also been found to be effective in simulating earthquake-induced structural pounding.
1. Introduction
During ground motions, buildings often collide with each other due to different dynamic characteristics, insufficient gap between them, and out-of-phase vibrations [1]. The phenomenon related to such collisions is often called the earthquake-induced structural pounding. Pounding is expressed by an instance of rapid strong pulsation, which may cause severe damage [2]. Consequently, the probability of structural interactions during earthquakes must be assessed. A number of researchers have studied the problem of pounding for different structural configurations under various ground motions:
Anagnostopoulos [3] was among the first scientists who investigated the earthquake-induced structural pounding by analyzing interactions between buildings in a row. He also described in detail the threats related to such collisions during ground motions [4].
Maison and Kasai [5] presented a formulation and simulated the multi degree-of-freedom equations of motion for floor-to-floor pounding between two 15-storey and 8-storey buildings. The influence of building separation, relative mass, and contact location properties was investigated.
Jankowski [6, 7] and Mahmoud et al. [8, 9] carried out a number of studies concerning structural pounding, related to both experimental and numerical aspects.
Komodromos et al. [10–12] studied pounding between seismically isolated buildings.
Barros and Khatami [13] examined a new model of impact to coordinate of results among numerical and experimental studies. They also suggested an approximate trend to select coefficient of restitution, which becomes equal to impact velocity [14]. Subsequently, a new equation of motion was suggested to simulate impact and figure damping terms out of collision [15].
Furthermore, some more recent numerical analyses have been carried out to study the influence of different parameters in pounding of buildings [16–20].
Nevertheless, there is still a need to investigate different models of structural pounding so as to verify their accuracy in the case of different configurations under different earthquake excitations. This concerns especially the values of impact forces during collisions which are often not studied in the analyses (or the analyses are simplified) since the investigations are rather focused on pounding-involved response under ground motions.
2. Contact Element Methods
The contact element is a special element (usually consisting of a spring and damper) to model impact between two colliding structures, which is widely used to simulate impact force. Impact is parametrically modelled in this way that when relative displacement exceeds the separation distance, the contact element is activated. The general formula for the impact force during collision can be expressed as follows (see [21]):(1)Fimpt=ks·δnt+cimpt·δ˙t,where ks is stiffness of spring, cimp(t) denotes damping of dashpot, and δ(t) and δ˙(t) describe lateral displacement and velocity, respectively. The power of n has been recommended to be 1 or 1.5 depending on the model considered. The damping coefficient cimp(t) is usually related to the coefficient of restitution CR which is defined as the ratio between the postimpact velocity δ˙rebound and the prior-impact velocity δ˙imp [21]:(2)0<CR=δ˙reboundδ˙imp<1.
Anagnostopoulos. The linear viscoelastic model of impact force (with n=1) was considered by Anagnostopoulos [3]. The following formulae were suggested [3, 22]:(3)cimpt=2·ζ·ks·mimjmi+mjζ=lnCRπ2+lnCR2,where mi,mj are the masses of colliding structures. The impact force versus time and lateral displacement for the linear viscoelastic model is shown in Figure 1.
Impact force versus time and lateral displacement for the linear viscoelastic model.
Jankowski. Jankowski [6] considered the nonlinear viscoelastic model when n=1.5. He proposed activating the damping term only during the approach period, in which most of the energy is lost. Therefore, the following formulae were considered [6]:(4)Fimpt=ks·δnt+cimpt·δ˙t⟶δ˙t>0 (approach period)Fimpt=ks·δnt⟶δ˙t≤0 (restitution period)(5)cimpt=2·ζ·ksδt·mimjmi+mj.It was suggested that the appropriate value of impact damping ratio ζ for a specified value of CR can be obtained numerically through iterative simulations in order to satisfy the relation between the postimpact and the prior-impact velocities defined by (2) [6]. The following formula was also proposed [9]:(6)ζ=952·1-CR2CRCR9π-16+16.An alternative approach might also be needed to consider the relation between the damping coefficient cimp(t) and the coefficient of restitution CR in the following form:(7)cimpt=α·CRβδ˙imp1-CRδ¨t·ks·δnt,where α and β are the parameters of the model obtained by fitting the experimental data using the method of the least squares (determined in this way typical values are α=0.01557,β=0.2706). As it can be seen from (7), the impact damping coefficient depends directly on the prior-impact velocity value, which can be obtained for each impact separately. The impact force versus time and lateral displacement for the nonlinear viscoelastic model is shown in Figure 2.
Impact force versus time and lateral displacement for the nonlinear viscoelastic model.
Muthukumar and DesRoches. A different model (Hertzdamp model) was considered by Muthukumar and DesRoches [23]. They presented a formula for the impact force during collision with n=1.5, which depends on three parameters, including coefficient of restitution, stiffness of spring, and prior-impact velocity, as (8)cimpt=ζ·δnt(9)ζ=3·ks1-CR24·δ˙imp.The impact force versus time and lateral displacement for the Hertzdamp model is shown in Figure 3.
Impact force versus time and lateral displacement for the Hertzdamp model.
Ye et al. Ye et al. [24] claimed that (9) is incorrect to calculate the damping ratio for the pounding simulation in order to evaluate impact between two buildings. They numerically indicated that the appropriate relation is [24](10)ζ=8·ks1-CR5·CR·δ˙imp.The impact force versus time and lateral displacement for the modified Hertzdamp model is shown in Figure 4.
Impact force versus time and lateral displacement for the modified Hertzdamp model.
3. Numerical Study
A parametric study has been conducted in order to verify the effectiveness of the nonlinear viscoelastic model of impact force during structural pounding described by (4) and (7). Two single degree-of-freedom systems (see Figure 5) with the gap size of d=2 cm have been used to model the behaviour of adjacent structures under ground motions.
Model of interacting structures.
The dynamic equation of motion for such a model can be expressed as [25](11)m100m2x¨1tx¨2t+C100C2x˙1tx˙2t+K100K2x1tx2t+Fimpt-Fimpt=-m100m2x¨gtx¨gt,where xi(t),x˙i(t),x¨i(t),Ci,Kiare the horizontal displacement, velocity, acceleration, damping coefficient, and stiffness coefficient for structure i(i=1,2), respectively, x¨g(t) stands for the acceleration input ground motion, and Fimp(t) is the pounding force which is equal to zero when δ(t)≤0 and is defined by (4), when δ(t)>0, where δ(t) is defined as(12)δt=x1t-x2t-d.In the numerical analysis, stiffness of each structure has been taken to beK1=K2=740MN/m and storey masses have been assumed to bem1=110 tons andm2=145 tons, respectively. The structural damping ratio of 0.05 has been considered in the analysis.
The dynamic analyses under the Parkfield (1966), San Fernando (1971), Kobe (1995), and El Centro (1940) earthquake records have been performed. These records have different contents of the excitation frequencies, different magnitude of the accelerations, and different time durations. Besides, their place of occurrence and geological conditions close to the epicentre are distinct. San Fernando earthquake had the highest Peak Ground Acceleration (PGA) among the four records discussed. The PGA of the earthquake amounted to 1.164 g, with an epicentre distance less than 12 km. The PGA of the Kobe earthquake was 0.7105 g and it was measured at a distance of 18.3 km. The PGA of the Parkfield earthquake was equal to 0.462 g (measured at a distance of 32 km). Finally, the PGA of the El Centro earthquake was equal to 0.347 g. All mentioned records have been normalized to investigate the effect of earthquake properties on pounding-involved structural response. The examples of the results of the numerical analysis in the form of the lateral displacement time histories under different earthquakes are shown in Figure 6. Additionally, the examples of the impact force time histories for the Kobe and the El Centro earthquakes are presented in Figure 7.
Lateral displacement time histories under different earthquakes.
El Centro
Kobe
San Fernando
Parkfield
Impact force time histories for the Kobe and El Centro earthquakes.
Using four different earthquake excitations, the peak lateral displacements, velocities, and accelerations have also been calculated for different structural periods of colliding structures. The results of the analyses are presented in Figure 8. They indicate that, with the increase in the structural period, the peak lateral displacements show the nonuniform increase trend. Among the ground motions analyzed, the Kobe record gives the maximum lateral peak displacement equal to 3.87 cm, while the minimum peak displacement of about 0.64 cm has been observed for the Parkfield earthquake. In the case of velocity, when the structural period is increased, the peak velocities are nearly the same in the range of 1–4 sec and after that the curves show a sudden decrease trend. Slightly different trend is observed for the El Centro record which demonstrates a slight increase from 5.1 m/s to 7 m/s and, subsequently, shows a sharp decline to 2.64 m/s when the structural period changes its value from 1 to 9 sec. Finally, the peak acceleration curves are quite stable at the beginning of analyzed range of structural period and then they show a substantial increase trend.
Peak lateral displacement, velocity, and acceleration with respect to structural period under different earthquakes.
3.1. Effect of Gap Size
In order to investigate the effect of separation distance between structures, a gap size has been varied from 0 to 8 cm. Figure 9 shows the effect of separation distance on the peak impact force under four different earthquake records. It can be seen from the figure that the curves follow an irregular decrease trend when the gap size increases. In the case of the San Fernando and Kobe ground motions, a sudden decrease is observed after passing a specific gap size value, while a slight declining tendency is visible for two other earthquake records.
Peak impact force with respect to gap size under different earthquakes.
3.2. Effect of Coefficient of Restitution
Different values of coefficient of restitution CR have been considered to investigate the impact forces between structures under different earthquakes. The results of the investigation showing the peak values of impact forces are presented in Figure 10. Similar trend can be observed for all analyzed excitations. The results show a uniform decrease in the force when the coefficient of restitution increases. For instance, the peak impact force for the Kobe earthquake is equal to 341 kN and 92 kN for CR=0.1 and CR=0.9, respectively.
Peak impact force with respect to coefficient of restitution under different earthquakes.
3.3. Effect of Impact Velocity
In order to obtain the responses and compare the results of peak impact forces, different values of impact velocity have been considered from the range 1–25 m/s. The relations between the peak impact force and impact velocity values under different earthquakes are presented in Figure 11. The results show a uniform increase in the peak impact forces with the increase in the impact velocity. The peak impact forces are nearly equal to zero for the velocity of 1 m/s and are as large as 1725 kN, 1578 kN, 1560 kN, and 1405 kN for the 25 m/s impact velocity under the El Centro, Kobe, San Fernando, and Parkfield earthquake, respectively.
Peak impact force with respect to impact velocity under different earthquakes.
3.4. Effect of Stiffness of Impact Spring
Stiffness of impact spring element is considered to be one of the most important parameters when the impact force during collision is calculated. The results of the parametric study showing the peak values of impact forces with respect to stiffness of spring are presented in Figure 12. It can be seen from the figure that the trend for all earthquakes analyzed is similar. The impact forces show nearly linear increase from 210 kN, 214 kN, 273 kN, and 352 kN to 2122 kN, 2242 kN, 2625 kN, and 3342 kN for the Parkfield, San Fernando, Kobe, and El Centro earthquake, respectively.
Peak impact force with respect to impact spring stiffness under different earthquakes.
4. Conclusions
In this paper, earthquake-induced pounding between two adjacent buildings has been studied by the application of single degree-of-freedom structural models. Impact has been numerically simulated with the use of a nonlinear viscoelastic model. Special attention has been focused on calculating values of impact forces during collisions which have significant influence of pounding-involved response under ground motions.
The results of the study indicate that the impact force time history depends substantially on the earthquake excitation analyzed. Moreover, the peak impact force during collision is much dependent on such parameters as gap size, coefficient of restitution, impact velocity, and stiffness of impact spring element. The nonlinear viscoelastic model of impact force with the considered relation between the damping coefficient and the coefficient of restitution has also been found to be effective in simulating pounding between structures during seismic excitations.
The conclusions of the study can be very valuable for the purposes of accurate modelling the phenomenon of earthquake-induced structural pounding. This concerns especially the issue of determination of the precise values of impact forces during collisions which are often not studied in the analyses (or the analyses are simplified) since the investigations are rather focused on pounding-involved response under ground motions. It can be considered as the most significant element of the analysis described in this paper, as compared to other relevant research studies.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
KasaiK.MaisonB. F.Building pounding damage during the 1989 Loma Prieta earthquake199719319520710.1016/s0141-0296(96)00082-x2-s2.0-0031105067JankowskiR.Assessment of damage due to earthquake-induced pounding between the main building and the stairway tower200734733934410.4028/0-87849-444-8.3392-s2.0-34250809549AnagnostopoulosS. A.Pounding of buildings in series during earthquakes198816344345610.1002/eqe.42901603112-s2.0-0023822752AnagnostopoulosS. A.Building pounding re-examined: how serious a problem is it?1996Oxford, UKPergamon, Elsevier ScienceMaisonB. F.KasaiK.Dynamics of pounding when two buildings collide199221977178610.1002/eqe.42902109032-s2.0-0026928014JankowskiR.Non-linear viscoelastic modelling of earthquake-induced structural pounding200534659561110.1002/eqe.4342-s2.0-18744392529JankowskiR.Theoretical and experimental assessment of parameters for the non-linear viscoelastic model of structural pounding2007454931942MahmoudS.JankowskiR.Modified linear viscoelastic model of earthquake-induced structural pounding201135C15162MahmoudS.Abd-ElhamedA.JankowskiR.Earthquake-induced pounding between equal height multi-storey buildings considering soil-structure interaction20131141021104810.1007/s10518-012-9411-62-s2.0-84880628862KomodromosP.PolycarpouP. C.PapaloizouL.PhocasM. C.Response of seismically isolated buildings considering poundings200736121605162210.1002/eqe.6922-s2.0-34648836965PolycarpouP. C.KomodromosP.Numerical investigation of potential mitigation measures for poundings of seismically isolated buildings20112112410.12989/eas.2011.2.1.0012-s2.0-79951879772KomodromosP.PolycarpouP.A nonlinear impact model for simulating the use of rubber shock absorbers for mitigating the effect of structural pounding during earthquake20124281100BarrosR. C.KhatamiS. M.Damping ratios for pounding of adjacent building and their consequence on the evaluation of impact forces by numerical and experimental models201322119131BarrosR. C.NaderpourH.KhatamiS. M.MortezaeiA. R.Influence of seismic pounding on RC buildings with and without base isolation system subject to near-fault ground motions2013113952NaderpourH.BarrosR. C.KhatamiS. M.A new model for calculating the impact force and the energy dissipation based on CR-factor and impact velocity20142214863ColeG.DhakalR.CarrA.BullD.An investigation of the effects of mass distribution on pounding structures201140664165910.1002/eqe.10522-s2.0-79954605552Yaghmaei-SabeghS.Jalali-MilaniN.Pounding force response spectrum for near-field and far-field earthquakes20121951236125010.1016/j.scient.2012.07.0122-s2.0-84868134092TubaldiE.BarbatoM.GhazizadehS.A probabilistic performance-based risk assessment approach for seismic pounding with efficient application to linear systems201236-37142210.1016/j.strusafe.2012.01.0022-s2.0-84857277373ZhaiC.JiangS.LiS.XieL.Dimensional analysis of earthquake-induced pounding between adjacent inelastic MDOF buildings201514229531310.1007/s11803-015-0024-3MadaniB.BehnamfarF.Tajmir RiahiH.Dynamic response of structures subjected to pounding and structure—soil—structure interaction201578466010.1016/j.soildyn.2015.07.002GoldsmithW.19601stLondon, UKEdward ArnoldAnagnostopoulosS. A.Equivalent viscous damping for modeling inelastic impacts in earthquake pounding problems200433889790210.1002/eqe.3772-s2.0-3042642198MuthukumarS.DesRochesR.A Hertz contact model with non-linear damping for pounding simulation200635781182810.1002/eqe.5572-s2.0-33646829138YeK.LiL.ZhuH.A note on the Hertz contact model with nonlinear damping for pounding simulation20093891135114210.1002/eqe.8832-s2.0-68149163529JankowskiR.Impact force spectrum for damage assessment of earthquake-induced structural pounding2005293-29471171810.4028/www.scientific.net/KEM.293-294.7112-s2.0-34249661499