Next Article in Journal
Evaluating the Impact of External Support on Green Building Construction Cost: A Hybrid Mathematical and Machine Learning Prediction Approach
Previous Article in Journal
Application of Glass Structures in Architectural Shaping of All-Glass Pavilions, Extensions, and Links
Previous Article in Special Issue
Study on Damping Performance of Hyperboloid Damper with SMA-Negative Stiffness
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Experimental Hybrid Simulation of Severe Aftershocks Chains on Buildings Equipped with Curved Surface Slider Devices

1
Civil Engineering and Architecture Deptartment, University of Pavia, 27100 Pavia, Italy
2
EUCENTRE Foundation, 27100 Pavia, Italy
*
Author to whom correspondence should be addressed.
Buildings 2022, 12(8), 1255; https://doi.org/10.3390/buildings12081255
Submission received: 28 June 2022 / Revised: 4 August 2022 / Accepted: 11 August 2022 / Published: 16 August 2022
(This article belongs to the Special Issue Building Structures and Seismic Technologies)

Abstract

:
In this research work the outcomes of a hybrid experimental campaign are analyzed, in order to evaluate the influence of aftershock events on the frictional response of sliding-based isolation devices for buildings. To achieve this, a hybrid testing framework was accordingly defined, by considering a numerical substructure, in terms of a simplified analytical model of a case study structure, and a physical substructure, as a full-scale Curved Surface Slider device, tested within the Bearing Tester System of the EUCENTRE Foundation Laboratory in Pavia (Italy). The tested isolator was equipped with a special sliding material, made up of a Poly-Tetra-Fluoro-Ethylene-based compound (PTFE), filled with carbon fibers and with a solid lubrication. The hybrid tests were performed, in terms of earthquake simulations, and the response of the base-isolated structural system was computed, by applying single-events, rather than aftershock chains. Results lead to a better understanding of the behavior of sliding-based seismic isolation systems, characterized by medium-to-high tribological properties, in terms of peak and residual displacements for both the single-event and the mean responses. Specifically, this work provides hybrid experimental evidence of the influence of an initial displacement offset on the overall behavior of the considered structural system.

1. Introduction

The knowledge about the effects induced by earthquake excitations on structural systems has grown exponentially in last decades, also as a result of a number of natural disasters occurring in several locations worldwide. Human losses and economic exposure in countries that experienced severe seismic events led to the need to understand all the possible strategies for the most efficient reduction of the vulnerability of building and bridge structures, in order to limit the overall seismic risk level. One of the most innovative solutions for the mitigation of the vulnerability of structural systems is represented by seismic isolation, which consists of the design of special devices that introduce a low-stiffness layer between the ground and the considered structure, which provide a period shift and a consequent reduction of the overall base shear and internal forces [1,2,3]. At the same time, the resulting high displacements are limited by a certain dissipative capacity of the installed isolators, which is fairly represented by the hysteresis of the force displacement response. Among the possible technical solutions, Curved Surface Slider devices are able to provide sufficiently low force responses, together with a high energy dissipation, through the frictional motion [4,5,6,7,8,9,10]. Research activities from both the numerical and experimental perspectives have led to a high level of accuracy in predicting the performance of such isolation systems, and more and more realistic results can be achieved [11,12,13,14,15]. Concerning experimental techniques, shake table tests, which generally provide the most realistic behavior of a structural system, can not be easily defined for base-isolated structures, without a proper scaling procedure, which my lead to unrealistic stress distributions inside the devices. Consequently hybrid simulations are currently widely used for the evaluation of the seismic response of not only base isolated buildings, but generally for any kind of structural and non-structural systems, by considering numerical and physical substructures [16,17,18,19,20,21]: consequently, much more realistic results can be achieved, extremely close to what could be obtained through full-scale shake table tests. Although sliding isolators have significantly good performance under seismic loading, residual displacements at the end of a single event could represent an important issue, due to the eventual high-frictional properties, specifically designed for some applications. Earthquakes are commonly considered as the mainshock, which is generally the event with the highest peak ground acceleration of a global seismic sequence. However, in some cases, aftershock events can have amplitudes more or less comparable to the corresponding mainshock value, and possible cumulative residual displacements can occur and the collapse of some devices can be achieved. As highlighted in recent research applications, the evaluation of residual displacements at the end of subsequent seismic events can lead to significantly high values [22,23,24]. Such residual deformations represent an initial offset for the device, with a non-negligible reduction of the maximum stroke along one direction. Numerical results show that if multiple earthquake excitations are subsequently applied, as a seismic chain of events, the residual displacement at the end of the considered ground motions gains with the number of the considered aftershocks [25,26]. Nonetheless, in the scientific literature, results are available for numerical simulations only, and no experimental data are publicly available for full-scale structural systems. The reason is mainly due to a number of issues related to the actual possibility of performing shake table tests on full-scale base-isolated systems, since the maximum geometrical characteristics of specimens do not allow to obtain realistic loading conditions for the implemented isolators, and, consequently, alternative testing strategies are needed [27,28].
The present research aims at investigating the influence of residual displacements on the overall response of a case study structure, evaluated through experimental hybrid tests. The numerical substructure consists of a simplified analytical model of the building, whereas a full-scale isolation device was considered as physical substructure, tested within the Bearing Tester System of the EUCENTRE Foundation Laboratory in Pavia (Italy). Results are provided by comparing the response of the system when single events and aftershock chains are applied, and the effects of possible initial displacement offset are discussed.

2. Case Study Structure

In this work, a case study structure was studied in order to evaluate the influence of severe aftershock events on the overall response of the isolation system. More specifically, a Reinforced Concrete frame structure was considered, laid on a reinforced concrete slab, which represents the actual interface between the superstructure and the isolation system. In agreement with the Italian Building Code, external loads were combined, in order to consider the seismic combination, by considering permanent structural (3.65 kN/m2) and non-structural (1.55 kN/m2) self-weights of the flooring system, 30% of the live load related to residential buildings (30% of 2.0 kN/m2–0.6 kN/m2), and the self weight of beams columns plates, which is automatically computed by the software. The concrete class is C40/50, which is associated with a characteristic value of cylindrical compression strength of 40 MPa, whereas the reinforcement steel grade is B450C, with a characteristic value of yielding strength of 450 MPa. All the building elements were modeled by means of linear-elastic frame elements, since the response of the superstructure is expected to be within the elastic range, and the slab interface was implemented by means of a mesh of shell elements. In Figure 1a 3D rendering of the overall structural system is shown.
Along both x and y directions four spans are designed, with equal span length (6 m each), with a constant interstorey height of 3 m. The reinforced concrete slab at the base of the building has thickness of 500 mm, and plan dimensions were obtained, by considering the plan development of the superstructure, with an additional offset of 1.5 m along both directions x and y directions (consequently, 27 m × 27 m). In order to consider the effective flexural stiffness of all the frame elements, moment of inertia of all the implemented sections was reduced through ad hoc reduction factors, which are higher for columns, in order to take into account the effect of the axial load. As a consequence, the linear-elastic branch of the bi-linear approximation of the capacity curve, returned by distributed-plasticity modeling, results aligned with respect to the capacity curve of the linear-elastic model, as reported in Figure 2.
The non-linear capacity curve is represented as the normalized base-shear with respect to the total mass of the building in fixed-base configuration vs. building total drift, as the ratio between the roof displacement and the total height: results were computed by means of a pushover analysis, which was performed with the software SeismoStruct [29]. In details, all frame elements were implemented with force-based inelastic frames, with distributed plasticity, defined according to specific non-linear constitutive laws for materials (Mander’s model for Concrete and Menegotto-Pinto’s model for steel reinforcement [30,31]); in addition, the inelastic frame elements related to concrete columns were defined with the proper confinement coefficients, which provide the correct influence of lateral reinforcements in the increase of both compressive strength and ductility capacities of the concrete core, aiming at obtaining the most realistic response for the building. More specifically, for both columns and beams Force-Based inelastic frame elements were implemented, represented by five Gauss’ integration points. The overall force-based formulation provide a strong form of the solution, and consequently a single frame element can be adopted for all the elements, and no longitudinal mesh has to be defined (which is fundamental for the proper approximation of Displacement-Based weak formulations). For all sections, a large number of fibers was set up, aiming at obtaining the most realistic representation of both concrete and steel rebar contributions. The capacity curve related to the fixed-base configuration of the building was returned by push-over analysis under displacement control, for the proper evaluation of possible softening behaviors after the achievement of the peak strength [32]. The adopted isolators consist of Curved Surface Slider devices, which can provide high dissipative capacity through the frictional response, together with a certain recentering behavior, induced by the spherical shape of the sliding surfaces [33,34,35,36,37].

3. Seismic Input

In order to perform earthquake simulations by using the Hybrid testing technique, a suite of seven one-directional natural seismic events was considered [38]. More specifically, the ground acceleration time series were selected, as ruled by the Italian Building Code [39], in agreement with the seismic hazard level defined for the construction site. For the adopted case study structure, the construction site is assumed to be located at L’Aquila (Italy), and Soil class B and topography category T1 were considered. Since this work focuses the attention on the behavior of the isolation devices, the Collapse Limit State was studied, and consequently 5% probability of exceedance in the reference life of the structural system was assumed. Thus, with a reference life period of 50 years, a return period of 975 years is obtained. The spectrum-compatibility of the adopted suite of natural records was studied, in agreement with the Italian Building Code. Hence, the mean response spectrum of the records selection, obtained by computing the average spectral coordinate at all period values, was ensured to be bounded between the lower and upper limits, corresponding to 90% and 130% respectively of the target value, within a period range between 0.15 s and 3.0 s. In Table 1 the selected seismic events are listed, together with some of the main characteristics.
It has to be noted that the selection of seismic events contains seven individual signals, even though some of them are related to the same earthquake; nonetheless, they were recorded at different stations, thus are associated with different stations, rather than along opposite directions. In order to be consistent with the characteristics of the construction site, the suite of natural seismic events was selected in the available databases of the software REXEL by limiting soil classification to the same properties of the target spectrum (i.e., Soil Class B). This is commonly performed in real applications, in order to consider earthquake excitations which are more comparable to what is expected at the construction site of the case study structure under investigation. Furthermore, the disaggregation analysis was considered for the proper selection of the seismic events by setting specific ranges for magnitude and epicentral distance which led to the highest contributions in the definition of the overall seismic hazard. It has to be noted that all these limitations and boundaries increase the probability to not achieve convergence, and maybe no selection can be obtained in some cases. In addition, records were scaled in order to better achieve the spectrum compatibility, also from a single-event perspective; moreover, all scale factors were limited between 0.5 and 2.0, in order not to obtain unrealistic signals, characterized by inconsistent balance between the frequency content and the considered Peak Ground Acceleration (PGA) values.
In Figure 3 graphical results for the spectrum compatibility study are provided in terms of single-event and mean response spectra, in comparison to the target, together with the code-based lower and upper bound.
Thanks to the aforementioned spectrum-compatibility study, the single-event discrepancy between the related response spectrum and the target spectral coordinates is very limited in the period range of interests for the designed structural system (values higher than 2.0 s). The selected records were applied both as single events and as three different chains of aftershocks earthquakes. It is highlighted that the adopted approach is the most demanding, since individual mainshocks events were applied consecutively: results are expected to provide a big picture of what happens in the very worst case.

4. Hybrid Simulation Framework

This section provides the description of the framework designed for experimental hybrid simulations of the previously described case study structure, performed at the Laboratory of EUCENTRE Foundation in Pavia (Italy). More specifically, the reinforced concrete building was numerically modeled through a statically condensed multi degree of freedom oscillator, by computing the actual stiffness and mass matrices from the full three-dimensional Finite Element Model (F.E.M.). On the other hand, the base-isolation system was represented by the response of a single physical full-scale isolator, which is assumed to be representative of the whole isolation system: in order to obtain the actual force response of the isolation system, the force feedback returned by the testing equipment is scaled by the ratio between the total structural weight and the applied vertical load to the physical device. Consequently, purely translational motions are considered, with no torsional movements of the superstructure [40,41].
The hybrid simulation was set up as an equivalent time history analysis, by applying at each time instant the selected ground motion time series to the dynamic system of the Multi Degree of Freedom (MDOF) system [42] (i.e., the numerical substructure); the resulting value of displacement at the isolation level is then applied to the full-scale device (i.e., the physical substructure), tested within the Bearing Tester System (BTS), and the consequent force feedback value is used for the next iteration in the dynamic system. The stepwise interaction between the numerical and the physical substructures leads to the final solution of the time integration of the dynamic systems, which consists of the displacement time series at all levels of the building, together with all the feedback signals recorded by the testing equipment of the bearing.
In order to obtain the most realistic results possible, a time scale factor equal to 8.0 was adopted, which has provided structural responses significantly close to real-time. Consequently, the motion applied to the physical device through the Bearing Tester System is represented by the same displacement amplitudes in comparison to the real-time numerical simulation, even though velocities values are scaled by 8.0. The assumption of the proper scale factor for experimental hybrid simulations could be crucial, especially for rate-dependent devices, such as the ones adopted in this study, which present the commonly known “velocity effect”. Nonetheless, thanks to the low scale factor, and given the high values of sliding velocity achieved during earthquake simulations, the device is subjected to motions which imply a constant friction coefficient, the peak velocity value being greater than the value correspondent to the end of the increasing transition branch of the frictional properties.

4.1. Physical Substructure

A Double Curved Surface Slider isolator [5] was considered as the physical substructure, being one of the most typical configuration of sliding-based devices adopted for real practice applications. In Figure 4, a section of the 3D rendering of the device is provided, which shows all the internal components. Such an isolator is made up of two individual spherical backing plates, which were realized with the same radius of curvature and with structural steel material S355JR. Within counterbore gaps of 5 mm stainless steel sliding surfaces are provided, with a thickness of 2 mm, and both are polished to mirror finish in order to achieve a roughness index Ra of 0.2 [40]. In the middle of the sliding surfaces, a non-articulated slider is installed, which is a unique steel block (material S355JR). The internal slider houses two circular sliding pads, made up of a special sliding material with the same diameter (160 mm), in order to consider the same frictional response at both the sliding interfaces. The overall equivalent radius of curvature of the device can be obtained by the summation of the radii of curvature of the spherical surfaces (1600 mm each), and by subtracting the height of the internal slider (120 mm): thus, the equivalent radius of curvature of the isolator is equal to 3080 mm.
Both the backing plates of the isolator have an external diameter equal to 548 mm, whereas the inner slider diameter is equal to 170 mm: consequently, the maximum displacement allowed by the geometrical configuration of the device is equal to 360 mm. A special sliding material was adopted within the sliding interfaces of the device, which consists of pigmented graded Polytetrafluoroethylene (PTFE) filled with carbon fibers, with a solid lubrication: such a sliding material showed significantly good frictional properties, together with a negligible break-away and stick slip friction coefficients and with low damage after several dynamic tests. During the hybrid simulations the physical device was subjected to vertical load equal to 400kN. In order to evaluate the tribological properties of the device, a number of characterization tests were performed, as ruled by the European Standard Code for Anti-Seismic devices UNI:EN15129:2009 [43], and consequently the reference friction coefficient was computed:
μ = E D C 4 D m a x W
EDC is the Energy Dissipated per Cycle, obtained as the integral of the force-displacement hysteretic loop, Dmax is the maximum displacement of the test and W is the applied vertical load. According to the assumed vertical load, a friction coefficient value equal to 0.10 (10%) was computed.
Such frictional characteristics are generally provided by medium to high tribological properties of the sliding material commonly adopted in real applications and is expected to lead to interesting results, since in this work special attention is focused on the influence of initial offset and residual displacements on the response of aftershocks events.

4.2. Numerical Substructure

The case study structure was numerically modeled as an equivalent Multi Degree of Freedom oscillator (MDOF), thanks to a special static condensation procedure, which was applied to the full 3D Finite Element Model (F.E.M.) of the structure ([41]; Figure 5).
The resulting mass and stiffness matrices allow to obtain dynamic properties significantly close to the ones returned by the three-dimensional F.E.M. model, in terms of periods, modal shapes and modal participating mass ratios, which were implemented within the commercial software SAP2000 [44]. The full model of the building was developed by considering linear-elastic frame elements for both beams and columns, by considering 30 GPa and 0.25 as Young’s modulus and Poisson’s ratio respectively. It has to be noted that a linear elastic model can fairly represent the response of the building when the isolation system is installed, since all standard codes allow to linearly model the superstructure, thanks to the vulnerability reduction induced by the isolators. The high level of accuracy of the adopted numerical model, despite the simplification of the overall system, is a direct consequence of the static condensation procedure, which actually returns the stiffness and mass characteristics of the structure as a function of the effective flexural behavior of both beams and columns. More specifically, in order to compute the i-th column of the condensed stiffness matrix, the horizontal translational degree of freedom at the i-th level of the system is unrestrained and a unit horizontal force is lumped at the same location. Such a reference point for the i-th floor corresponds to the location of the center of mass, which is connected to the related story by means of a rigid diaphragm constraint, so that the out-of-plane stiffness of beams is considered. On the other hand, the remaining degrees of freedom (centers of mass of the remaining floors) are restrained. Thus, the i-th column of the matrix can be computed by dividing the vector containing the reactions and the lumped force by the horizontal displacement at the i-th unrestrained floor, originated by the applied lumped load.
Such a stiffness matrix is computed without considering any force response coming from the isolation level, since the hysteretic behavior is considered as the horizontal force feedback returned by the testing equipment, within the time integration algorithm which rules the experimental hybrid simulation. Thus, the resulting stiffness matrix is representative of the linear elastic response of the superstructure, by accounting also for the horizontal rigid body translational motion. With the adopted procedure, similar dynamic properties of the full 3D model of the considered structural system can be obtained, since all the implemented structural elements of the model are able to deform, and no restrained behavior is assigned. In agreement with all the aforementioned assumptions, the overall structural system was experimentally evaluated, by performing hybrid earthquake simulation, according to the following dynamic system of motion equations:
M ¯ ¯ ( u ¨ 0 _ u ¨ 1 u ¨ 2 u ¨ 3 ) + C ¯ ¯ ( u ˙ 0 _ u ˙ 1 u ˙ 2 u ˙ 3 ) + K ¯ ¯ ( u 0 _ u 1 u 2 u 3 ) + F i s ( 1 _ 0 0 0 ) = M ¯ ¯ ( 1 _ 1 1 1 ) x ¨ g
where:
  • M ¯ ¯ is the mass matrix, represented by a diagonal matrix can be obtained, which has the components of the main diagonal equal to the summation of the assembled mass values of all the points of a given storey;
  • C ¯ ¯ is the viscous damping matrix, defined according to a Rayleigh formulation, with a low damping ratio, namely 2%, in order not to underestimate the overall response;
  • K ¯ ¯ is the stiffness matrix returned by the applied static condensation procedure;
  • u i , u ˙ i and u ¨ i are the translational degrees of freedom and the time derivatives at the i-th level, relative to the ground;
  • F i s is the force response feedback, which is stepwise returned by the Bearing Tester System;
  • x ¨ g is the applied ground motion.
Since earthquake simulations were performed, the structure was assumed to be in equilibrium at the beginning of the motion; thus, zero initial conditions were considered for both displacement and velocity values at all degrees of freedom. Thus, the solution of the dynamic system was returned by the time-integration algorithm, which considers a stepwise interaction between the numerical substructure (i.e., the MDOF oscillator representing the building) and the physical substructure (i.e., the full-scale device tested within the Bearing Tester System), by considering a scaled time axis (time scale factor: 8.0 -Figure 6).
The presented hybrid simulation framework was implemented within the MATLAB Environment & Simulink [45], and allows to time-integrate the overall dynamic system, by considering the stepwise interaction of both the numerical and physical substructures. Since hybrid tests are performed to simulate earthquake excitation, the presented system is analyzed through a special time-integration algorithm, in order to perform an equivalent Non-Linear Time History Analysis. The non-linearity of all the implemented isolators is fairly and realistically represented by the experimental response of the physical substructure. Consequently, the ground acceleration time series represent the input signals for the dynamic system of the overall base-isolated building, and the displacement response at all levels of the building can be computed through the time-integration of the whole set of dynamic equilibrium equations. More specifically, at each time step, the displacement at the isolation layer becomes the control command for the Bearing Tester System (Figure 7), which is able to apply such a deformation value to the physical isolator, by considering a displacement-control method.
The force feedback response returned by the testing equipment, scaled by the effective ratio between the total weight of the structural system and the vertical load applied to the device, provides an updated force contribution for the dynamic equilibrium equation at the isolation level, which can be further time-integrated. Thanks to the presented algorithm the response of the whole base-isolated building can be obtained.

5. Experimental Results and Discussion

In this section the experimental results are shown for all the hybrid earthquake excitations. Special attention is focused on the global hysteretic response of the single-event cases, in order to experimentally check the tribological properties of the adopted sliding material. Then, the comparison of the displacement time series at the isolation level returned by the single-event excitations are graphically compared to three different aftershock chains, which were defined by considering three different orderings of the adopted earthquakes. The same comparison was studied by analyzing peak displacements and residual displacements as a direct consequence of initial displacement offset when aftershocks are applied.
In Figure 8, results of the adopted single events are provided in terms of hysteretic responses of the isolation layer.
More specifically, on the horizontal axis, the isolation displacement is reported, whereas on the vertical axis the isolation force response is considered, normalized with respect to the total weight of the overall structural system. This representation of the hysteretic response allows to notice the resulting frictional properties, which corresponds to the numerical value of normalized force at zero-displacement configuration: it should be noted that the desired friction coefficient was achieved, approximately equal to 10% in all cases. The mean displacement demand achieves values of about 83 mm, with a coefficient of variation of 43%. Concerning the sliding velocity, a mean value of 46 mm/s was found, related to the experimental response of the physical substructure, which corresponds to approximately 370 mm/s in the numerical real-time simulation, with a coefficient of variation of 35%. With the adopted sliding material being a PTFE-based compound, results of previous characterization testing campaigns have shown that 50 mm/s can be considered as a threshold value related to the achievement of a constant asymptotic value of friction coefficient, regardless of the applied vertical load at the physical device. Consequently, it could be assessed that negligible differences in the presented results could be obtained also for lower time scale factors, which correspond to simulations closer to real time (i.e., time scale factors 4.0, 2.0 or 1.0). In Figure 9, Figure 10 and Figure 11, graphical results are provided, in terms of experimental displacement time series at the isolation level, for both aftershocks chains and single events, by overlapping the obtained signals. It can be noted that all single-events (colored lines) start from the zero-displacement configuration; on the other hand, all the performed aftershock chains (black lines) presents an initial displacement offset, which causes an increase of the peak displacement response. This aspect is amplified by the numerical value of the experimental friction coefficient (10%), which generally leads to non-negligible residual displacements at the end of all the events, even though the recentering contribution induced by the spherical shape of the sliding surfaces provides a certain limitation.
Furthermore, after the strong motion range of all the seismic events, a smoothly decreasing branch of the residual displacement can be detected: this behavior could be attributed to the response of the building, which is subjected to negligible ground acceleration values, and consequently a free-vibrations phase could be considered; such a vibrating motion provides an additional recentering contribution, which gradually decreases the residual displacements for both the single events and the aftershocks chains. The presented results provide evidence that the residual displacement at the end of all aftershocks does not increase significantly, even though all events start with an initial displacement offset.
In order to better understand the overall behavior of the considered structural system, further analyses were computed, by considering the initial offset, the peak and residual displacements and the drift response of the building.
In Figure 12 the building total drift and peak normalized base shear values were analyzed for all the simulations, in order to ensure if the superstructure was subjected to limited plastic deformations, and consequently to prove that the linear elastic modeling assumption for the building is reasonable.
As can be noted all drift values are averagely equal to the yielding reference value (9%) among the selected seismic events of the adopted suite of natural earthquakes; only for events #5 and #7 would some limited plastic deformations be considered, even though the variation percentages between such values and the target elastic drift do not exceed 20% in the worst cases. On the other hand, normalized values for building base shear are lower than the corresponding target value at the yielding point of the capacity curve in all cases. Thus, the overall linear elastic modeling assumption for the building can be considered reasonable for the presented study.
In Figure 13, the ratio between the initial offset and the peak displacement of the chain aftershocks are provided (left), together with the percentage variation between the peak displacement response of the single event and the aftershock chains.
Results show that the initial offset at the beginning of all the aftershock events can be a significantly high percentage of the peak displacement response of the isolation system. Approximately the same behavior can be noticed in most cases, regardless of the actual sequence of seismic events, which seems to suggest that the initial offset is strictly related to the overall characteristics of the selected ground motions, with no influence of their position within the aftershocks chains. This is not expected to hold for cases such as the MW 7.3 earthquake on 9 March 2011 that was followed by the stronger MW 9.0 Tohoku earthquake on 11 March 2011 in Japan (additional details can be found in [46,47,48] which provide natural time analysis of seismicity). In the worst case, for Event #4, the initial offset is approximately 60% of the peak displacement response of the subsequent event. Nonetheless, as previously assessed for the displacement time series, the initial offset does not increase as the aftershock events occur. As a direct consequence, the ratio between the single event and the aftershock chains peak displacement responses does not increase as well, and very similar variations can be obtained, regardless of the actual position of each event within the performed chains. The initial offset lead to variations in the peak displacement response bounded between −20% and −10% in all cases, which provide evidence that the initial offset leads to an increasing lateral deformation of the isolation layer, which does not gain as all the aftershocks occur.
In Figure 14 results are presented in terms of ratio between the residual displacement for single event and for aftershock chains for all the selected earthquakes.
Additionally for residual displacement response, no significant influence of the actual position of all events within the considered aftershock chains can be averagely detected. Furthermore, as for the peak displacement response, the residual deformation returned by the application of single events is lower than the corresponding value related obtained in the aftershock chain. On the other hand, it can be noted that for some events no relevant differences can be found from the considered approaches. Such an issue seems to highlight again that the main characteristics of the suite of selected ground motions, and no significant influence of the actual sequence of events can be clearly assessed.
Finally, peak displacement and interstory drift responses were analyzed for all the levels of the considered case study structure, as shown in Figure 15 and Figure 16, respectively, in terms of variation percentages between the single event and the aftershock chains values.
Approximately same conclusions drawn for the isolation layer (Level #0) can be considered also for all the levels of the superstructure, which shows that the single-event response is lower in comparison to the application of the same events within aftershock chains.
In addition, it can be assessed that the variation percentage looks to decrease, if higher levels of the building are considered, and the maximum variation at the isolation level (Level #0: −20%) significantly decreases (Level #3: −10%). Concerning the drift response, much lower variation percentages can be computed, averagely bounded between ±5%, and no clear trend can be defined. Consequently, the influence on the building response of the application of a single event, rather than aftershock chains, seems to be negligible, according to the presented experimental results.

6. Conclusions

The presented research provides an investigation of the behavior of friction-based seismic isolators when subjected to subsequent earthquake excitations. To this aim, a special framework for experimental hybrid tests was defined in order to obtain the most realistic results possible. The hybrid tests were set up as seismic simulation of a case study structure, base-isolated by means of Curved Surface Slider devices. More specifically, a numerical substructure was considered within the experimental hybrid algorithm, consisting of a simplified lumped-masses oscillator with the same dynamic properties of the full 3D F.E.M. model; on the other hand, a full-scale device was tested as physical substructure. The time-integration method implemented within the hybrid testing algorithm considers a time scale factor equal to 8.0, which is considerably close to real-time simulations and realistic results can be obtained. The isolator was equipped with an innovative sliding material, filled with carbon fibers and solid lubrication, which provides 10% friction coefficient at the applied vertical load. The results were obtained, by performing hybrid earthquake simulations on both single-events and with aftershock chains. By analyzing the results, that the following can be concluded:
  • The outcomes of the displacement response at all levels showed that the application of a single event generally leads to lower peak values, in comparison to aftershock sequences of earthquakes, due to the initial displacement offset.
  • Both the initial offset and the residual displacements do not increase as all the subsequent aftershock events occur, but numerical values remain limited and do not gain over the application of the seismic sequence.
  • The variation percentage between the single event and the aftershock chains peak displacement responses at all levels is averagely equal to −20% at the isolation system, and decreases to −10% as upper floors are considered.
  • This behavior can be assessed, regardless of the location of a single-event within the adopted aftershock chains.
  • Concerning the building response, limited variations were computed for interstory drifts at all levels, bounded between ±5%, which suggests that the superstructure is not significantly affected by the application of sequences of seismic events rather than single ground motions.
The outcomes of the present work seem to suggest that initial displacement offsets could lead to a certain unexpected discrepancy between the single-event and seismic chain responses, with consequent limitations in the total stroke of the isolation devices. Nonetheless, this issue would need further investigations, especially from the experimental standpoint, by also analyzing also other seismic isolation systems commonly adopted in real applications, together with several case study structures, in order to come up with generalized conclusions. A better understanding of the recentering behavior of isolation devices would be valuable for the revisions of the main standard codes for both anti-seismic devices and structural analysis.

Author Contributions

A.P. Funding acquisition and supervision; I.L. Software and methodology; M.F. Project administration, investigation, writing-original draft and writing review and editing. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by Italian Civil Protection (Convenzione DPC-EUCENTRE 2018: ReLUIS project 2019-2021).

Acknowledgments

Part of the current work was carried out under the financial support of the Italian Civil Protection, within the frameworks of the Executive Project 2018 (Project 3–Assessment of the seismic isolation of building structures through hybrid tests with numerical substructuring) and the national Research Project DPC–ReLUIS (National Network of Laboratories of Seismic Engineering) 2019-2021, WP15–Isolation and Dissipation.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. De Luca, A.; Guidi, L.G. State of art in the worldwide evolution of base isolation design. Soil Dyn. Earthq. Eng. 2019, 125, 105722. [Google Scholar] [CrossRef]
  2. Dall’Asta, A.; Leoni, G.; Micozzi, F.; Gioiella, L.; Ragni, L. A Resilience and Robustness Oriented Design of Base-Isolated Structures: The New Camerino University Research Center. Front. Built Environ. 2020, 6, 50. [Google Scholar] [CrossRef]
  3. Tsiavos, A.; Sextos, A.; Stavridis, A.; Dietz, M.; Dihoru, L.; Di Michele, F.; Nicholas, A. Low-cost hybrid design of masonry structures for developing countries: Shaking table tests. Soil Dyn. Earthq. Eng. 2021, 146, 106675. [Google Scholar] [CrossRef]
  4. Tsiavos, A.; Schlatter, D.; Markic, T.; Stojadinovic, B. Shaking table investigation of inelastic deformation demand for a structure isolated using friction-pendulum sliding bearings. Structures 2021, 31, 1041–1052. [Google Scholar] [CrossRef]
  5. Constantinou, M.; Mokha, A.; Reinhorn, A. Teflon bearings in base isolation II: Modeling. J. Structral Eng. 1990, 116, 455–474. [Google Scholar] [CrossRef]
  6. Constantinou, M.C.; Tsopelas, P.; Kasalanati, A.; Wolff, E. Property Modification Factors for Seismic Isolation Bearings; Report No.99-0012; Multidisciplinary Center for Earthquake Engineering Research: Buffalo, NY, USA, 1999. [Google Scholar]
  7. Calvi, P.M.; Calvi, G.M. Historical development of friction-based seismic isolation systems. Soil Dyn. Earthq. Eng. 2018, 106, 14–30. [Google Scholar] [CrossRef]
  8. Messina, A.; Miranda, E. Characterization of the Frictional Behavior of Steel-Polymer Interfaces with Pronounced Stick-Slip Effect for Use in Seismic Isolation. Bull Earth. Eng. 2022. [Google Scholar] [CrossRef]
  9. Maegawa, S.; Liu, X.; Itoigawa, F. Discussion of Stick-Slip Dynamics of 2DOF Sliding Systems Based on Dynamic Vibration Absorbers Analysis. Lubricants 2022, 10, 113. [Google Scholar] [CrossRef]
  10. O’Reilly, G.J.; Yasumoto, H.; Suzuki, Y.; Calvi, G.M.; Nakashima, M. Risk-based seismic design of base-isolated structures with single surface friction sliders. Earth. Eng. Struct. Dyn. 2022. [Google Scholar] [CrossRef]
  11. De Domenico, D.; Gandelli, E.; Quaglini, V. Adaptive isolation system combining low-friction sliding pendulum bearings and SMA-based gap dampers. Eng. Struct. 2020, 212, 110536. [Google Scholar] [CrossRef]
  12. Losanno, D.; De Domenico, D.; Madera Sierra, I.E. Experimental testing of full-scale fiber reinforced elastomeric isolators (FREIs) in unbounded configuration. Eng. Struct. 2022, 260, 114234. [Google Scholar] [CrossRef]
  13. Dolce, M.; Cardone, D.; Croatto, F. Frictional behavior of steel-PTFE interfaces for seismic isolation. Bull. Earthq. Eng. 2005, 3, 75–99. [Google Scholar] [CrossRef]
  14. Fenz, D.; Constantinou, M.C. Behaviour of the double concave friction pendulum bearing. Earthq. Eng. Struct. Dyn. 2006, 35, 1403–1424. [Google Scholar] [CrossRef]
  15. Kumar, M.; Whittaker, A.S.; Constantinou, M.C. Characterizing friction in sliding isolation bearings. Earthq. Eng. Struct. Dyn. 2015, 44, 1409–1425. [Google Scholar] [CrossRef]
  16. Park, J.; Strepelias, E.; Stathas, N.; Kwon, O.S.; Bousias, S. Application of hybrid simulation method for seismic performance evaluation of RC coupling beams subjected to realistic boundary condition. Earth. Eng. Struct. Dyn. 2021, 50, 375–393. [Google Scholar] [CrossRef]
  17. Schellenberg, A.H.; Becker, T.C.; Mahin, S.A. Hybrid shake table testing method: Theory, implementation and application to midlevel isolation. Struct. Contr. Health Monit. 2017, 24, e1915. [Google Scholar] [CrossRef]
  18. Zhang, Z.; Park, J.; Kwon, O.S.; Sextos, A.; Strepelias, E.; Stathas, N.; Bousias, S. Hybrid Simulation of Structure-Pipe-Structure Interaction within a Gas Processing Plant. J. Pipe. Syst. Eng. Pract. 2021, 12, 16. [Google Scholar] [CrossRef]
  19. Cheng, M.; Becker, T.C. Performance of unscented Kalman filter for model updating with experimental data. Earth. Eng. Struct. Dyn. 2021, 50, 1948–1966. [Google Scholar] [CrossRef]
  20. Tsokanas, N.; Zhu, X.; Abbiati, G.; Marelli, S.; Sudret, B.; Stojadinović, B. A Global Sensitivity Analysis Framework for Hybrid Simulation with Stochastic Substructures. Front. Built Env. 2021, 71. [Google Scholar] [CrossRef]
  21. Abbiati, G.; Lanese, I.; Cazzador, E.; Bursi, O.S.; Pavese, A. A computational framework for fast-time hybrid simulation based on partitioned time integration and state-space modeling. Struct. Contr. Health Monit. 2019, 26, e2419. [Google Scholar] [CrossRef]
  22. Quaglini, V.; Gandelli, E.; Dubini, P.; Limongelli, M.P. Total displacement of curved surface sliders under non seismic and seismic actions: A parametric study. Struct. Control Health Monit. 2017, 24, e2031. [Google Scholar] [CrossRef]
  23. Berton, S.; Infanti, S.; Castellano, M.G.; Hikosaka, H. Self-centring capacity of seismic isolation systems. Struct. Contr. Health Monit. 2007. [Google Scholar] [CrossRef]
  24. Katsaras, C.P.; Panagiotakos, T.B.; Kolias, B. Restoring capability of bilinear hysteretic seismic isolation systems. Earth. Eng. Struct. Dyn. 2008. [Google Scholar] [CrossRef]
  25. Christopoulos, C.; Filiatrault, A.; Folz, B. Seismic response of self-centring hysteretic SDOF systems. Earth. Eng. Struct. Dyn. 2002. [Google Scholar] [CrossRef]
  26. Macrae, G.A.; Kawashima, K. Post-Earthquake Residual Displacements of Bilinear Oscillators. Earth. Eng. Struct. Dyn. 1998. [Google Scholar] [CrossRef]
  27. Quaglini, V.; Gandelli, E.; Dubini, P. Experimental investigation of the re-centring capability of curved surface sliders. Struct. Contr. Health Monit. 2017. [Google Scholar] [CrossRef]
  28. Quaglini, V.; Dubini, P.; Gandelli, E. Effect of Friction on the Re-Centring Capability of Sliding Bearings with Curved Surfaces. In Proceedings of the 7AESE-7th International Conference on Advances in Experimental Structural Engineering, Pavia, Italy, 6–8 September 2017. [Google Scholar] [CrossRef]
  29. Seismosoft. SeismoStruct 2018–A Computer Program for Static and Dynamic Nonlinear Analysis of Framed Structures. 2018. Available online: http://www.seismosoft.com (accessed on 28 June 2022).
  30. Mander, J.B.; Priestley, M.J.N.; Park, R. Theoretical stress-strain model for confined concrete. J. Struct. Eng. 1988, 114, 1804–1826. [Google Scholar] [CrossRef]
  31. Menegotto, M.; Pinto, P.E. Method of analysis for cyclically loaded R.C. plane frames including changes in geometry and non-elastic behaviour of el-ements under combined normal force and bending. In Proceedings of the Symposium on the Resistance and Ultimate Deformability of Structures Acted on by Well Defined Repeated Loads, International Association for Bridge and Structural Engineering, Zurich, Switzerland, 28 June 1973; pp. 15–22. [Google Scholar]
  32. Antoniou, S.; Pinho, R. Advantages and Limitations of Force-based Adaptive and Non-Adaptive Pushover Procedures. J. Earth. Eng. 2004, 8, 497–522. [Google Scholar] [CrossRef]
  33. Mazza, F.; Mazza, M. Sensitivity to modelling and design of curved surface sliding bearings in the nonlinear seismic analysis of base-isolated r.c. framed buildings. Soil Dyn. Earthq. Eng. 2017, 79, 951–970. [Google Scholar] [CrossRef]
  34. Mokha, A.; Contantinou, M.C.; Reinhorn, A.M. Teflon Bearings in Aseismic Base Isolation: Experimental Studies and Mathematical Modeling; Report No. 88-0038, MCEER; National Center for Earthquake Engineering Research: Buffalo, NY, USA, 1988. [Google Scholar]
  35. Castaldo, P.; Amendola, G. Optimal DCFP bearing properties and seismic performance assessment in nondimensional form for isolated bridges. Earth. Eng. Struct. Dyn. 2021, 50, 2442–2461. [Google Scholar] [CrossRef]
  36. Castaldo, P.; Amendola, G. Optimal sliding friction coefficients for isolated viaducts and bridges: A comparison study. Struct. Contr. Health Mon. 2021, 28, e2838. [Google Scholar] [CrossRef]
  37. Yang, T.; Calvi, P.M.; Wiebe, R. Numerical implementation of variable friction sliding base isolators and preliminary experimental results. Earthq. Spectra 2020, 36, 767–787. [Google Scholar] [CrossRef]
  38. Iervolino, I.; Galasso, C.; Cosenza, E. REXEL: Computer aided record selection for code-based seismic structural analysis. Bull. Earthq. Eng. 2009, 8, 339–362. [Google Scholar] [CrossRef]
  39. Ministero delle Infrastrutture e dei Trasporti. Norme Tecniche per le Costruzioni, Decreto Ministeriale del 17 Gennaio 2018; Italian Building Code; Ministero delle Infrastrutture e dei Trasporti: Rome, Italy, 2018.
  40. Furinghetti, M.; Pavese, A.; Quaglini, V.; Dubini, P. Experimental Investigation of The Cyclic Response of Double Curved Surface Sliders Subjected to Radial and Bidirectional Sliding Motions. Soil Dyn. Earthq. Eng. 2019, 117, 190–202. [Google Scholar] [CrossRef]
  41. Furinghetti, M.; Lanese, I.; Pavese, A. Experimental assessment of the seismic response of a base isolated building through hybrid simulation technique. Front. Built Environ. 2020, 6, 33. [Google Scholar] [CrossRef]
  42. Chopra, A.K. Dynamics of Structures Theory and Applications to Earthquake Engineering; Prentice Hall: Upper Saddle River, NJ, USA, 1995. [Google Scholar]
  43. CEN. Comité Européen de Normalisation TC 340, European Code UNI EN 15129:2009 Anti-Seismic Devices; European Committee for Standardization: Brussels, Belgium, 2009. [Google Scholar]
  44. Computers and Structures Inc. SAP2000 Ultimate: Structural Analysis Program. Version 21.2.0. 2021. Available online: http://www.https://www.csi-italia.eu/ (accessed on 26 April 2022).
  45. MathWorks. MATLAB Webpage. MathWorks. Available online: https://ch.mathworks.com/ (accessed on 7 February 2018).
  46. Varotsos, P.A.; Skordas, E.S.; Sarlis, N.V. Fluctuations of the entropy change under time reversal: Further investigations on identifying the occurrence time of an impending major earthquake. EPL 2020, 130. [Google Scholar] [CrossRef]
  47. Varotsos, P.A.; Sarlis, N.V.; Skordas, E.S. Self-organized criticality and earthquake predictability: A long-standing question in the light of natural time analysis. EPL 2020, 132. [Google Scholar] [CrossRef]
  48. Sarlis, N.V.; Skordas, E.S.; Varotsos, P.A. Order parameter fluctuations of seismicity in natural time before and after mainshocks. EPL 2010, 91. [Google Scholar] [CrossRef]
Figure 1. 3D rendering of the case study structure.
Figure 1. 3D rendering of the case study structure.
Buildings 12 01255 g001
Figure 2. Capacity curve of the building.
Figure 2. Capacity curve of the building.
Buildings 12 01255 g002
Figure 3. Spectrum-compatibility of the selected records.
Figure 3. Spectrum-compatibility of the selected records.
Buildings 12 01255 g003
Figure 4. Photographic and 3D representations of the tested physical device.
Figure 4. Photographic and 3D representations of the tested physical device.
Buildings 12 01255 g004
Figure 5. Reference scheme for the dynamic system of motion.
Figure 5. Reference scheme for the dynamic system of motion.
Buildings 12 01255 g005
Figure 6. Testing procedure for experimental hybrid simulations.
Figure 6. Testing procedure for experimental hybrid simulations.
Buildings 12 01255 g006
Figure 7. Testing setup for the physical substructure of the SHAKE Lab at EUCENTRE Foundation in Pavia (Italy).
Figure 7. Testing setup for the physical substructure of the SHAKE Lab at EUCENTRE Foundation in Pavia (Italy).
Buildings 12 01255 g007
Figure 8. Single-event hysteretic response of the isolation system.
Figure 8. Single-event hysteretic response of the isolation system.
Buildings 12 01255 g008
Figure 9. Displacement results of aftershock chain #1.
Figure 9. Displacement results of aftershock chain #1.
Buildings 12 01255 g009
Figure 10. Displacement results of aftershock chain #2.
Figure 10. Displacement results of aftershock chain #2.
Buildings 12 01255 g010
Figure 11. Displacement results of aftershock chain #3.
Figure 11. Displacement results of aftershock chain #3.
Buildings 12 01255 g011
Figure 12. Building total drift (left) and base shear (right) responses.
Figure 12. Building total drift (left) and base shear (right) responses.
Buildings 12 01255 g012
Figure 13. Initial offset ratio (left) and peak displacement variation (right) for aftershock chains.
Figure 13. Initial offset ratio (left) and peak displacement variation (right) for aftershock chains.
Buildings 12 01255 g013
Figure 14. Ratio between single-event and aftershock chain residual displacements.
Figure 14. Ratio between single-event and aftershock chain residual displacements.
Buildings 12 01255 g014
Figure 15. Peak displacement ratio between single event and aftershock chain at all levels.
Figure 15. Peak displacement ratio between single event and aftershock chain at all levels.
Buildings 12 01255 g015
Figure 16. Peak interstory drift ratio between single event and aftershock chains at all levels.
Figure 16. Peak interstory drift ratio between single event and aftershock chains at all levels.
Buildings 12 01255 g016
Table 1. Selection of natural seismic events.
Table 1. Selection of natural seismic events.
Event [#] 1Station IDEarthquake NameComponentMwFault MechanismEpicentral Distance [km]Original PGA [g]Scaled PGA [g]Scale Factor [#]
1ST_58065Loma Prietay6.9oblique27.590.5120.3790.74
2HECHector Minex7.1strike-slip28.610.3360.4201.25
3ST_58065Loma Prietax6.9oblique27.590.3240.3561.10
4AQGL’Aquila mainshockx6.3normal4.390.4450.7581.70
5LGPCLoma Prietax6.9oblique18.750.5860.5861.00
6LGPCLoma Prietay6.9oblique18.750.9650.4820.50
7KARGazlix6.7reverse12.780.7170.7891.10
1 [#] units correspond to dimensionless parameters.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Furinghetti, M.; Lanese, I.; Pavese, A. Experimental Hybrid Simulation of Severe Aftershocks Chains on Buildings Equipped with Curved Surface Slider Devices. Buildings 2022, 12, 1255. https://doi.org/10.3390/buildings12081255

AMA Style

Furinghetti M, Lanese I, Pavese A. Experimental Hybrid Simulation of Severe Aftershocks Chains on Buildings Equipped with Curved Surface Slider Devices. Buildings. 2022; 12(8):1255. https://doi.org/10.3390/buildings12081255

Chicago/Turabian Style

Furinghetti, Marco, Igor Lanese, and Alberto Pavese. 2022. "Experimental Hybrid Simulation of Severe Aftershocks Chains on Buildings Equipped with Curved Surface Slider Devices" Buildings 12, no. 8: 1255. https://doi.org/10.3390/buildings12081255

APA Style

Furinghetti, M., Lanese, I., & Pavese, A. (2022). Experimental Hybrid Simulation of Severe Aftershocks Chains on Buildings Equipped with Curved Surface Slider Devices. Buildings, 12(8), 1255. https://doi.org/10.3390/buildings12081255

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop