Next Article in Journal
An Easy and Effective Technique for Silicone Facial Implant Insertion and Fixation to Periosteum
Previous Article in Journal
The Protective Effect of Static Magnetic Fields with Different Magnetic Inductions against Fluoride Toxicity Is Related to the NRF2 Signaling Pathway
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Seismic Risk Assessment of a Novel Self-Centering Precast Concrete Frame under Near-Fault Ground Motions

1
Nanjing Institute of Technology, Nanjing 211176, China
2
Tongji Architectural Design (Group) Co., Ltd., Shanghai 200092, China
3
Key Laboratory of Concrete and Prestressed Concrete Structures of Ministry of Education, Southeast University, Nanjing 210096, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2020, 10(18), 6510; https://doi.org/10.3390/app10186510
Submission received: 29 July 2020 / Revised: 11 September 2020 / Accepted: 12 September 2020 / Published: 18 September 2020
(This article belongs to the Section Civil Engineering)

Abstract

:

Featured Application

The novel self-centering precast concrete (SCPC) frame with hysteretic dampers has favorable seismic behavior and achieves the desirable demands at the design basis earthquake (DBE) and the maximum considered earthquake (MCE) seismic levels under the near-fault earthquake waves. In addition, the novel connection shows desirable energy consumption and self-recovering capacities, which is applicable to the precast structure in near-fault seismic region.

Abstract

The damage to structures caused by the velocity pulse effect of near-fault earthquake waves cannot be ignored, yet there are few studies on the risk assessment of seismic performance for precast concrete frame under near-fault earthquake waves. A novel self-centering precast concrete (SCPC) frame with hysteretic dampers is proposed to obtain great self-recovering and energy consumption characteristics. To accurately assess the seismic behaviors of the novel SCPC frame under the near-fault earthquake waves, a prototype structure is modelled and elastoplastic dynamic analysis is conducted at the design basis earthquake (DBE) and the maximum considered earthquake (MCE) seismic levels. Incremental dynamic analysis and the vulnerability analysis are performed. Annual and 50-year exceeding probabilities of the novel SCPC frame are calculated afterwards. In addition, the reinforced concrete (RC)frame and the traditional SCPC frame are also modelled, whose section sizes, reinforcements arrangement and seismic intensity are consistent with the novel SCPC frame. The dynamic time-history analysis at the two seismic levels are also carried out for two types of frames. The analysis results demonstrate that the novel SCPC frame has great seismic performance and low seismic risk possibility under the near-fault earthquakes loading.

1. Introduction

At present, the velocity pulse effect caused by near-fault earthquake waves has attracted the extensive attention of researchers. Although the Imperial Valley earthquake in 1979, the Northridge earthquake in 1994, the Kobe earthquake in 1995, the Kocaeli earthquake and Chi-Chi earthquake in 1999 are all moderate earthquakes, they caused serious damage and loss of life and damage to structures [1]. From the existing strong earthquake records, the predominant velocity pulse emerges from either forward directivity or fling step effects [2,3,4,5]. Forward directivity effects show that the component perpendicular to the fault has a strong velocity pulse and fling step effects show residual and irrecoverable displacement due to the relative movements of ruptures [6,7,8]. Hence, the velocity pulse effect of near-fault earthquake waves can increase the seismic energy and cause serious damage to the structure and leave a permanent displacement. Different from the common far-field earthquake waves, short-duration and high-energy characteristics of the near-fault earthquake waves are a great potential destructive threat and thus increase the seismic requirement of structures [9].
The unique velocity pulse effect has been studied by many researchers on different structural systems [10]. Fibio et al. studied the influence of near-fault earthquake waves on the isolation structure with high damping rubber bearing [11]. At the same time, the comprehensive effects of horizontal and vertical ground motion components are considered and the numerical analysis of the bi-reinforcement concrete buildings is carried out according to the European seismic code. Mortezaei et al. introduced the influence on the seismic performance of the original reinforced concrete (RC)structure and RC structure strengthened by Fiber Reinforced Polymer (FRP)at the near-fault earthquake waves loading [12]. The results showed that the strengthening technology of FPR can improve the rigidity, strength and lateral displacement capacity of the structure. Babak et al. analyzed the important influence of the velocity pulse effect of near-fault earthquake waves on the seismic performance of elastic frame and inelastic frame [13]. The results showed that the early yielding occurs at the higher storied of the structure when the natural vibration period is longer than the pulse period. With the increase of the seismic intensity, the ductility demands move to the bottom stories. The change rule is opposite when the natural vibration period is shorter than the pulse period. Chin-Hsiung Lou et al. studied the influence of near-fault earthquake waves from the Chi-Chi earthquake on the seismic performance of the structure [14]. Based on the test study, the response of a single degree of the freedom inelastic system is calculated by the improved nonlinear hysteretic model. According to the nonlinear elastic response spectrum, the variation of base shear coefficient is discussed. Reza Sehhati et al. analyzed the impulse effect of near-fault earthquake waves on the multi-story building and 54 near-fault earthquake wave tests were carried out for three kind of multi-story building [15]. Besides, when the ratio of pulse period and the basic period are in the range of 0.5~2.5, the equivalent pulse model can better reflect the influence of forward directionality on the structure. From the aforementioned study, seismic risk assessments of the precast concrete structures under near-fault earthquake waves using the fragility analysis has been little reported in the literature [16].
Since the velocity pulse effect of near-fault earthquake waves can generate more significant inelastic deformation, many structural systems with high ductility and effective seismic energy dissipation have been developed in recent years. The self-centering reinforced concrete frame using the unbonded prestressed tendons have been proposed, which can maintain good self-restoring ability and retain minimal residual deformation [17,18,19]. Priestley and Tao proposed the concept of assembling precast components with unbonded tendons and the tests of connection and frame were also carried out afterwards [20]. The results showed the residual deformation of the joint was negligible and shear capacity of the joints were obviously improved. EI-Sheikh et al. established the analysis models of some 6-story and 4-span precast concrete frames using the elaborate finite elements in light of the different seismic regions and sites [21]. Pushover analysis and elastoplastic dynamic analysis based on 16 seismic waves were also carried out. The results displayed frames were in nonlinear elastic state, the damage was small and collapse could be avoided. Cheok et al. evaluated the seismic behaviors of a precast hybrid frame using the mix beam-column connections, the seismic energy consumption capacity of joints was improved through the yielding of ordinary reinforcements. A versatile hysteretic model simulating the nonlinear behavior of the joint area was proposed [22]. Morgen et al. used a new kind of friction damper installed on the beam-column joint with unboned tendons and proposed the seismic design method of the frame with the novel connections [23]. In addition, pushover analysis of the multi-story frame with the damper was carried out under the low cyclic loadings. The results showed that the frictions can provide supplemental energy dissipation, the residual deformation of the connection was minimal and the friction device was easy to replace. Geoffery et al. conducted the low cyclic load tests of two internal and external beam-column connections specimens with the energy dissipation devices (HF2V). The results showed the HF2V device can effectively enhance the energy dissipation and self-recovering characteristics when HF2V devices were reliably anchored [24]. In addition, Geoffery et al. established the mechanical model of precast beam-column connection using the HF2V devices under low cyclic repeated load [25]. Lu et al. designed the two story two way and three story three way scaled self-centering precast concrete frames and conducted a shake table test [26]. The experimental results showed that the seismic behavior and self-restoring capacity of the frame were excellent. They recommended that seismic energy dissipation devices can be added to enhance the seismic behavior of the precast frame.
To further enhance the energy consumption capacity, maintain the self-recovering capacity and reduce the local concrete damage; this study develops a new kind of hybrid self-centering precast (SCPC) frame joint with hysteretic dampers. In this novel connection, Post-tensioned (PT) tendons undertake moment and maintain the self-recovering capability after the ground motion loadings. The hysteretic dampers are arranged at the beam end to dissipate seismic energy and have no interference with the floor slab. To investigate the seismic risk assessment of the novel SCPC frame under near fault earthquake waves, a prototype frame with the novel SCPC connection is modelled and dynamic time-history analysis is performed at the design basis earthquake (DBE) and the maximum considered earthquake (MCE) seismic levels. Incremental dynamic analysis and vulnerability analysis are carried out afterwards and annual and 50-year exceeding probabilities of two limit states are obtained. In addition, the elastoplastic dynamic analysis of the RC frame and the traditional SCPC frame are also performed. Considering the near-fault’s strong velocity pulse effect, the analysis results validate the great seismic performance and low seismic risk possibility of the novel SCPC frame.

2. Basic Characteristics of Novel SCPC Connection

2.1. Construction of Connection

The novel SCPC hybrid joint is a new ductile connection, the deformations of precast beams and columns are within the elastic range and the inelastic damage of beams and column is minimal. The inelastic damage is usually concentrated in an easily replaceable energy dissipation device. Figure 1a illustrates the detailed construction of the novel joint. One duct is arranged at the central axis of the precast beam section and unbonded tendons are arranged along the central axis of the beam horizontally. PT tendons not only assemble the precast beam and column in the construction stage but also resist the bending moment in the use stage under the seismic actions. If the moment is large beyond the critical opening moment, the gap–opening of the hybrid joint has occurred. The length and the stresses of PT tendons increase subsequently. The residual deformation is eliminated after the earthquake due to the characteristic of the self-recovering of prestressed tendons provided. Shear force at the beam end is undertaken by the friction at the beam-column interface caused by tendons and the friction bolts. The steel plate is shop fabricated in the column through the embedded part to reduce the local concrete damage. As shown in Figure 1a, the hysteretic dampers are installed at the upper and lower part of the beam end and welded between the steel connectors and the steel plates. The yield point of steel used in hysteretic dampers is low and the yield point is only 235 N/mm2. In addition, the slot area of long strip on the steel plate will firstly produce in-plane plastic deformation. Hence, the hysteretic damper can enter the yielding state faster under the seismic loadings and then dissipate seismic energy by generating large plastic deformation after yielding. Compared with the out-of-plane buckling hysteretic damper, the in-plane buckling hysteretic damper has the advantage of large initial stiffness, high bearing capacity and good energy dissipation capability. The four steel connectors are bolted on the two sides of the beam and every steel connector uses two high strength friction bolts. The other side of the steel connector is welded with the hysteretic damper, in order to fix the damper and limit the damper’s displacement. Therefore, the mechanical principle of the novel connection is—shear force at the beam end is taken by friction and the moment is undertaken by unbonded prestressed tendons and hysteretic dampers, PT tendons provide the self-recovering and hysteretic dampers improve energy consumption characteristics.

2.2. Mechanical Behavior Analysis of Hysteretic Dampers

The hysteretic dampers improve the seismic energy dissipation capability of the novel beam-column connection due to their good plastic hysteresis performance, which generates large stresses and large deformations after yielding. Figure 2 shows the construction of the hysteretic damper. For the novel frame, the high initial stiffness of the hysteretic damper can increase structural stiffness to limit the structural displacement under the minor earthquake and generate significant plastic deformation after yielding to dissipate the seismic energy under the major earthquake. Therefore, the hysteretic damper must have relatively large initial stiffness and the good seismic energy capability. The section size of the hysteretic damper is determined according to the required stiffness and yield strength. The steel plate of the hysteretic dampers with a rectangular slot only occur in-plane deformation and can effectively avoid out of plane bulking. Three-dimensional solid element models are established using ABAQUS software to clearly determine the hysteretic behavior of the novel SCPC joint. The eight-node, hexahedral linear reduced integral solid elements C3D8R are adopted to accurately simulate a hysteretic damper. The 2-node linear 3D truss elements T3D2 are used to simulate the reinforcement and unbonded tendons. Figure 3a displays the three-dimensional numerical simulation model of the new joint in ABAQUS and the hysteretic curve of new joint with a hysteretic damper under the cyclic loading, respectively. Energy consumption coefficient β E is an important performance index to measure the energy consumption capacity of novel SCPC connection [27]. β E is the ratio of the hysteresis loop area of the SCPC joint and that of the two line elastoplastic system with the same bending strength. From Figure 3b we can see that the hysteretic curve M- θ r is a typical double flag-shape, where θ r is the connection rotation at the beam-column interface. The hysteretic loops become much plumper and the energy consumption coefficient β E is 0.42, which indicates that the energy consumption capacity of the new connection with hysteretic dampers is sufficient.
To validate the numerical model and study the seismic behavior of the novel joint, the cyclic loading test of a full-sized novel SCPC connection specimen with hysteretic dampers was conducted. The details of the cyclic loading tests setup is shown in Figure 4. The hinge support was installed on the strong floor and the precast column was hinged. The displacement of the column top was constrained by the support beam and reaction wall. A 50t MTS actuator was fixed on the strong floor to apply the cyclic loading on the beam end. The horizontal displacement transducers were installed on the upper and lower surfaces of the beam end to measure the deformation of the beam end. M was the moment at the beam end. The connection rotation θr was calculated by the absolute length of opening and compression between the upper and lower of beam end and the distance of the horizontal transducers. Figure 5 shows the comparison of hysteresis curve between the numerical simulation and test results. The numerical simulation results agree with the experimental data on the whole, such as gap opening and closing, energy consumption and self-centering characteristics.
For the novel beam-column connection, the simplified hysteretic model of connection under cyclic loads is shown as Figure 6. As shown in Figure 6a, the theoretical hysteresis curve of the self-centering system taken by PT tendons is modelled with an idealized bilinear elastic spring. The analysis demonstrates hysteretic damper has ideally a much higher initial stiffness, it can be modelled as elastic plastic model. As shown in Figure 6b, the theoretical hysteresis curve of hysteretic dampers is an elastic-plastic model and the stiffness after yielding is zero. Hence, combined with the PT tendons and hysteretic dampers, the novel connection can be modelled as the bilinear elastic spring taken parallel with elastic plastic model. The theoretical hysteresis curve of novel hybrid joint is modelled as a typical double flag-shape hysteretic model.
The mechanical principle of the novel hybrid joint with hysteretic dampers is defined as the gap opening and closing alternately under the near-fault earthquake wave actions. Figure 7 shows the mechanical principle of novel connection composites when it bears the positive moment under seismic actions—where T is the resultant tension of all the unbonded tendons; Nt and Nb are the horizontal resultant forces of the upper and lower hysteretic dampers, when the beam-column joint subjected to the positive moment under the loadings, the upper hysteretic damper bear compressive force and the lower hysteretic damper bear tensile force; V c is the shear force taken by the fiction, reflecting the fictions caused by prestressed tendons and high strength bolts, respectively; d1 is the distance between the top edge of the beam and the resultant point of PT tendons; d t and d b are the distances between the center point of the top and bottom hysteretic dampers and the top edge, respectively; The moment M at the beam end includes the moment provided by unbonded tendons and the hysteretic dampers. N P T is the pre-pressure of PT tendons at the end beam, μ 1 is the fiction coefficient of the beam-column interface, N v b is the shear bearing capacity of single high strength bolt on the steel connection, the formula is:
M = T d 1   -   N t d t   +   N b d b
V c = F f = N P T × μ 1 + 4 × 2 × N b v .
The value of initial resultant force of unbonded tendons is T 0 , when the bending moment at beam-column full contact interface exceed the critical opening moment M I O under the external load, the joint gap has opened, the formula of M I O is as follows:
M I O = T 0 d 1   -   N t d t +   N b d b

3. Structural Model

3.1. Prototype Frame Description

To assess the seismic risk of the novel SCPC building under the near-fault earthquake waves, a five-story four-span prototype precast concrete frame with hysteretic dampers is designed for analysis. According to the design requirements of near-fault ground motion, the distance between the fault fracture surface and the site is within 15 km. Near-fault earthquake waves are chosen from the PEER (Pacific Earthquake Engineering Research Center) record database, all the near-fault earthquake waves only considered the velocity pulse effect generated by the forward-directivity effect. Based on Chinese Standard GB 50010-2010 [28] and GB 50011-2010 [29], the site class of the prototype structures is defined as group II, the characteristic site period Tg is 0.40 s and the seismic intensity is categorized as 7.0. The original structural layout is regular, with dimensions of 24 m × 24 m. The total height of the frame is 32.8 m including 4.8 m for the first floor and 4.0 m for the upper floors. In the light of Chinese Standard Code [30], the dead load is 7.0 kN/m2 for the roof while it is 5.0 kN/m2 for other floor. Live loads for all the stories are 2.0 kN/m2. The calculation method of the representative value of gravity load of the novel frame is dead loads plus 0.5 times of live loads for the novel SCPC frame. The P-Δ effect of frame caused by gravity load is fully considered (the additional internal force and deformation produced by the precast concrete columns).
Figure 8b displays the configuration and the reinforcing details. The grade of concrete strength is C40. The longitudinal reinforcement, stirrup and waist bars use HRB400, the corresponding elastic modulus is 2 × 105 MPa. The nominal sectional area of one unbonded tendon is 140 mm2 and six tendons are arranged in one duct. The yield strength of PT strands is fpy = 1675 MPa and the ultimate strength of PT strands is fpy = 1860 MPa, the corresponding modulus of elasticity is 1.95 × 105 MPa. The steel used for the hysteretic dampers and steel plated are Q235 and Q335, the corresponding yield strengths are 235 MPa and 335 MPa, respectively. For the beam, four longitudinal reinforcements of 18 mm are arranged at the bottom and the bottom, respectively. Two waist bars of 10 mm are arranged and tie bars of 10 mm are arranged with a spacing of 250 mm. For the column, the diameters of longitudinal reinforcements arranged in the column cross sections are 20 mm. The diameter of stirrups is 8 mm diameter and the spacing of the stirrups is 100 mm in the stirrup densification area. One duct use the polyvinyl chloride (PVC)pipe with the diameter of 40 mm arranged at the beam center and four unbonded tendons do not grout in the duct. The configuration of the steel connector is shown in Figure 8c. The steel connectors are welded between the steel plate and the hysteretic dampers and bolted on the two sides of the beam using high strength friction bolts to limit the dampers’ displacements. For the hysteretic dampers and steel connectors, the thickness of all the steel plates is 10 mm.
This study uses performance-based seismic design methods to assess the seismic performance of a novel structure from elasticity to partial or overall collapse. According to the Federal Emergency Management Agency (FEMA)356 and characteristics of the novel SCPC frame, the design approach uses two performance state limits—Immediate occupancy (IO) and Repairable (RE) [31]. In addition, two seismic levels are considered—the design basis earthquake (DBE) and the maximum considered earthquake (MCE), whose earthquake return period are 475 years and 2475 years, respectively. The maximum horizontal seismic influence coefficient at the DBE level is about half of that at the MCE level. To assess the damage degree of the novel frame accurately, the displacement evaluation index, including the peak inter-story drift ratio and the peak residual inter-story drift ratio, are used. Therefore, the design target of structural limit states at two performance levels displays as follows.
Under the DBE earthquake level, the novel frame must fit the IO limit state. The structure damage level is very small; it can be reused without repair. When the beam-column connection gap is open, the concrete of the connection is in a state of eliminating pressure, the stress at the beam end declines to zero under the seismic actions. Minor hairline cracking appeared in a few places in precast columns and slight yielding phenomena occurred in the column. The allowable peak inter-story drift ratio cannot exceed in 1.0% and the peak residual inter-story drift ratio cannot exceed in 0.2% [30].
Under the MCE earthquake level, the novel frame must fit the RE limit state. The extensive cracks and plastic hinges are appeared in the column elements. Minimal yielding at the column bases is allowed, PT tendons are always elastic and do not yield. If the PT tendons yield, the connection will lose the restoring function and generate large residual deformation. The overall damage level is moderate and can be reused after repairing. The allowable peak inter-story drift ratio cannot exceed 2% and the residual story drift ratios cannot exceed 0.5%. The analysis from the existing earthquake damage reports shows—if the peak residual inter-story drift ratio exceed 0.5%, it does not make sense to repair the damage components from the viewpoint of economic efficiency [31].

3.2. Selected Near-Fault Ground Motions

The near-fault earthquake waves usually mean the distance between the rupture fault and building site is less than 20 km, epicenter intensity is less than 6.5 degrees and peak acceleration is more than 0.15 g [32]. To accurately assess the seismic performance and energy consumption capacity of the novel SCPC frame, 21 near-fault earthquake waves which are collected in the 12 earthquakes from PEER (Pacific Earthquake Engineering Research Center) record database were chosen for the elastoplastic dynamic analysis. The substantial observed seismic data demonstrate the near-fault ground motions are obviously different from the far-field ground motions. Because the rupture is close to building site, high-frequency damping is minimal and the velocity and displacement time-history of near-fault earthquake waves have the large amplitude and short duration. However, the velocity pulse effect of the near-fault earthquake waves means forward-directional effect and fling step effect. The directional effect can be divided into forward directional effect, backward directional effect and neutral directional effect and directional effect mainly refers to forward directional effect. Forward directional effect is the main cause of pulse effect and the damage of structures will be aggravated by forward-directional effect. Directional effect can make seismic energy accumulate in a short time and cause impact ground motion, which is reflected in large peak value, obvious pulse waveform and short duration of earthquake on time-history curve.
The fling step effect is generated by the relative movement of two fault rupture, which raises a single direction half pulse or one side pulse form in the velocity time-history curve and generates residual deformation on the ground [33]. Since the earthquake wave records in the PEER database are processed, the displacement is zero at the end of the record, so there is no record reflecting the slip effect. At present, structural damage caused by the fling step effect can only be solved by the method of seismic conceptual design. Near-fault earthquake waves usually only considered the forward directivity effect. The information of 21 ground motions is listed in Table 1, which includes the peak ground acceleration (PGA), peak ground velocity (PGV), peak ground displacement (PGD). The distance from the rupture fault to the structure site is less than 11 km and magnitudes Mw ranges from 6.19 to 7.51. The ground acceleration, displacement and velocity time-histories of near-fault ground motions KJM090 and H-E07230 are shown in Figure 9. The adjusted peak ground acceleration of each near-fault ground motion is equivalent to the two seismic levels. From the aforementioned research, the Rayleigh-type 3% damping ratios of the novel SCPC frame is based on the first and second modes of the vibration [34]. The corresponding acceleration response spectrum of each near-fault earthquake wave is shown in Figure 10 and the mean spectrum and the target spectrum are also shown in Figure 10.

4. Finite Element Study of Novel SCPC Joint

To evaluate the seismic behaviors of the novel SCPC frame, a 5-story prototype building is modelled as a 2D plane numerical model using the OpenSees software and the nonlinear time-history analysis was performed [35]. When the gap-opening occurred under the load, the sections do not fit the plane section assumption; the mechanical behavior of the novel connection is complex. The connection analysis model is based on the following basic assumptions—(1) the joint only considers the load in the plane and has the deformation in the plane without torsion; (2) shear failure of beam and column can be avoided through the reasonable stirrup arrangement; (3) the hysteretic dampers, beam and column do not slip in the whole loading process.
Figure 11 displays the accurate finite element model of the novel SCPC joint, which mainly consists of nonlinear beam column elements, zero length elements and truss elements. Two different types of constitutive relationships were used for concrete materials of the connection—(1) unconfined concrete in protective cover (without stirrup restraint); (2) concrete in core area (stirrup restraint). They are simulated by Concrete 01 material model consisting of three parts—nonlinear ascending section, linear descent section and horizontal section [36]. Concrete 01 material model do not consider the tensile strength of concrete and the concrete crushing, residual strength is accounted for. Considering the amplification effect of stirrup on confined concrete, there are various confined concrete model to calculate the confined concrete strength enhancement for both normal-strength and high-strength concrete [37,38,39,40]. The amplification factor of compressive strength and strain is 0.2 according to the parameters such as stirrup ratio and stirrup spacing. The steel fiber including steel bars and PT tendons are assigned to the Steel 02 material model, which are based on Giuffre-Menegotto-Pinto model [41]. Steel 02 material model has an initial stress, the Bauschinger effect and isotropic reinforcement effect can be reflected in the cyclic loading. As shown in Figure 12, all the beams and columns are modelled by Nonlinear Beam Column element, which consider the spread of plasticity along the element through a number of integration points. Each element is divided by many fiber cross sections by integration points and each fiber was given a uniaxial stress-strain relationship. The hysteretic dampers arranged at the upper and lower at the beam-column interface are simulated by zero length spring elements between node 14, 15 and 16, 17 [42]. The gap opening and closing behavior is also modelled by two pairs of zero length elements. The element is assigned to Elastic PP Gap material model with compression-only material properties. Although the gap is opened, the hysteretic dampers and beam-column interface do not slip. When the gap is opened, the sum of horizontal relative displacements of the upper and lower rotation points is Δgap, the vertical distance of the two points is h, so the relative connection rotation is θ r = Δ g a p / h . Three degrees of freedom of node 15 and node 17 are coupled with node 9, three degrees of freedom of node 14 and node 16 are coupled with node 4 and node 3. The rigid link elements (which is used to construct a single multipoint constraint) are selected to reflect the coupling characteristics between the above nodes. In the modelling of the beam-column interface, the effect of critical region length was considered. The critical region length is the distance between node 9 and node 10, where the most severe damage will take place [43,44]. The constitutive relationships of the concrete and steel materials model in the critical region are assigned to compression only material properties in the critical region. Truss element model between node 7 and node 12 are used to simulate the unbonded tendons at the beam center. Steel 02 material model with initial stress is adopted to simulate the constitutive relationship of the truss element due to prestressed tendons have the initial force. Truss element only consider the axial force and axial deformation, ignore geometric distortion in the vertical direction. In addition, the three degrees of node 10 and node 11 are coupled with node 13 respectively and use the equal Degree of freedom (DOF)constraints [45].
The finite element model method of the whole new SCPC frame with hysteretic dampers is the same as the model method of the new SCPC frame beam-column joint described in the above study. The model analyze is firstly conducted and the fundamental period of the novel SCPC frame is 1.04 s.

5. Seismic Risk Assessment of Novel SCPC Frame

5.1. Overall Responses to Near-Fault Earthquakes

Nonlinear time-history analysis is applied to the five-story four-span novel SCPC frame numerical model at the DBE and MCE earthquake levels using OpenSees software. The horizontal seismic influence coefficient of the prototype frame at the DBE and MCE levels are 0.117 and 0.234 calculated by the seismic intensity, site category, design earthquake grouping, the fundamental period and damping ratio according to Chinese Standard Code [29]. The gravity load of each story is evenly distributed at the central node of the core area of the beam-column connection. In order to calculate the residual deformations of the novel frame, 20 s last-time free vibration analysis of the structure continued after applying each near-fault ground motion action. From the aforementioned analysis results, the peak inter-story drift ratios and peak residual inter-story drift ratios are selected as the main evaluation index to define IO and RE limit state under the near-fault earthquake wave loadings.
Figure 13 displays the top floor displacement time-history curves of the near-fault earthquake waves NO.10-LGP090, NO.13-H-E04230, NO.18-CPM090 respectively. The peak value of top floor displacements are 23.5 mm, 20.8 mm and 32.4 mm under the DBE level and the peak value of top floor displacements are 54.2 mm, 46.7 mm and 69.8 mm under MCE level, respectively. The peak value of top floor displacements occur about in 11 s, 7 s and 4 s for the above three waves, which indicates the displacement responses of the frame are relatively large at the beginning and the displacements reduce significantly because of the self-centering capability for the novel SCPC frame. Hence, the final residual displacements for the novel frame are minimal under each ground motion. The final residual top floor displacements for the precast frame are 0.23 mm, 0.15 mm and 0.41 mm for NO.10-LGP090, NO.13-H-E04230, NO.18-CPM090 under the MCE level, respectively. The maximum top floor displacements, the times the maximum response occurs and the residual top displacements under 21 earthquake waves are all listed on Table 2. The above analysis results present the novel SCPC frame has greater lateral stiffness and good self-centering capability, the top floor displacement responses are in a reasonable range under DEB and MCE levels.
The peak inter-story drift ratios for each earthquake wave and for each story subjected to DBE and MCE levels under the near-fault ground motions are displayed in Figure 14. The maximum seismic responses for each floor, such as inter-story drift ratio, residual inter-story ratio, inter-story shear force and connection relative rotation, are the average values of the peak responses under 21 earthquake waves. The peak inter-story drift ratios of the novel frame are 0.72% and 2.03% corresponding for 21 near-fault earthquake waves at the DBE and MCE levels, respectively. The mean value of peak inter-story drift ratios of the frame for 21 earthquake waves are 0.34% and 0.81% subjected to the two levels, which are below the limiting requirements. The peak inter-story drift ratios for each story happen on the second story under the two seismic levels. The maximum responses from the third and the fifth story are equal to 0.29%, 0.22%, 0.15% at the DBE level and 0.68%, 0.51% and 0.29% at the MCE level, respectively. The peak inter-story drift ratios for each story at MCE level are about twice that at DBE level. It indicates that the degree of gap-opening at the beam-column interface has significantly increased at the MCE level, the bending stiffness of the connection has reduced.
Figure 15 shows the maximum residual inter-story ratio for each story at the DBE and MCE levels for the 21 near-fault earthquake waves. The peak residual inter-story drift ratios of the novel frame for all the stories are minimal, satisfy the limiting value which cannot exceed in 0.2% for IO limit state and 0.5% for RE limit state. The maximum residual inter-story ratios both occur on the first floor, which are only 3.85 × 10−5 rad and 1.87 × 10−4 rad at DBE and MCE level and the peak residual inter-story drift ratios of the above stories are similar and close to zero. It indicates the novel frame with hysteretic dampers still has excellent self-centering ability under the near-fault seismic wave actions.
Figure 16a and b show the maximum inter-story shear force for each earthquake wave and for each story under the two levels for the 21 near-fault earthquake waves, respectively. The peak inter-story shear force are 631 kN and 930 kN corresponding to 21 earthquake waves at the DBE and MCE levels, respectively. The peak inter-story shear force for each story have occurred on the first story under the two levels. The maximum responses of the first story are 365 kN and 650 kN at the DBE and MCE levels, respectively and the disparities between them are less than twice. From the second floor to the fifth floor, the inter-story shear forces decrease gradually. The maximum responses form the second to fifth story are 343 kN, 302 kN, 256 kN and 202 kN under the DBE level and 562 kN, 515 kN, 401 kN and 305 kN under the MCE level, respectively. It indicates that the novel SCPC frame has sufficient lateral stiffness to undertake more loadings.
Figure 17 displays the maximum connection relative rotation and the maximum standard PT tendons force for each story at the two seismic levels dunder the near-fault earthquake waves loading. TN defines as the standardized PT force, which is the ratio of PT force and PT yielding force T y . Figure 17a indicates—the maximum connection relative rotation of the novel SCPC frame for each story is minimal. The peak value is 8.5 × 10−4 rad for the second floor under the DBE level and it means the beam-column interface do not open basically. The peak connection relative rotations increase obviously at the MCE level, the peak value is 3 × 10−3 rad occurred at the first floor. The maximum connection rotations decrease gradually and the value of the top floor is only 1.3 × 10−3 rad. Figure 17b indicates—the maximum PT tendons force is relatively small for each story under DBE level and increase significantly at the MCE seismic level. In order to avoid great residual displacement in the core zones of beam-column connection, the PT tendons should not yield. The peak values of the maximum standardized PT tendon forces are about 0.48 times and 0.6 times of the yielding tendon force at the DBE and MCE levels. From the analysis results, unbonded tendons force for each story is in the elastic stage, far from reaching its yield strength under the MCE level.
The hysteretic relationship of the beam end moment M and connection relative rotation θ r is the important index to assess the energy consumption characteristics of the novel SCPC structure. Figure 18 shows the M- θ r curve of the joint for the earthquake wave KJM090 loading under the MCE level. The joint is located in the right of the second span of the second story. The figure shows the peak connection rotation is only 0.003 and the peak connection moment is 620 KN. In addition, the beam end M- θ r relation curves of the novel joint in the middle column is stable and shows typical double-flag shape. This sufficiently indicates that the novel frame joint with hysteretic dampers has effective energy consumption and self-recovering characteristics.

5.2. Incremental Dynamic Analysis (IDA)

Incremental dynamic analysis (IDA) is one of potential parameter analysis methods used in performance-based seismic engineering to assess the seismic behavior of the structure with the increase of the ground motion intensity. The technique includes a series of nonlinear dynamic analysis using one or more seismic ground motions multiplied by a series of proportional coefficients, to cover a whole dynamic performance of the structures with multiple intensities from elastic state to collapse state. To plot the IDA curves, the intensity measure (IM) of earthquake waves and the engineering demand parameters (EDP) are required. In this paper, 21 near-fault ground motions recommended by PEER record database are selected to conduct the IDA and assess the seismic behavior of novel SCPC structure. Select the proportional strength measurement index such as the acceleration response spectrum value corresponding to the first model period Sa (T1) as the intensity measure (IM), compared with other intensity measures such as the PGA and PGV. The EDP select the maximum inter-story drift ratio θmax and the maximum residual inter-story drift ratio θr,max to reflect seismic response at the near-fault earthquake waves. The acceleration response spectrum value corresponding to the first model period Sa (T1) change from 0.05 g to 0.6 g, at an interval 0.05 g. When applying the nonlinear dynamic history analysis on the novel SCPC frame, we can obtain the EDP index θmax and θr,max for each near-fault earthquake wave. The relationship between the EDP θmax and θr,max and the acceleration response spectrum value Sa (T1) are displayed in Figure 18. The figures illustrate the peak inter-story drift ratio and the peak residual inter-story drift ratio for each earthquake wave has a certain degree of dispersion and the dispersion increases with the increase of the seismic intensity.
For the probabilistic estimation of demand, the relationship between engineering demand parameters (EDP) and intensity measure Sa (T1) shows as follows:
D = a ( S I ) b ε ,
where a and b denote the regression coefficient of linear regression analysis, SI refers to acceleration response spectrum value corresponding to the first model period Sa (T1), ε is lognormal distribution random variable with a median of 1.0 and logarithmic standard deviation of βd.
According to the Formula (4), regression analysis is conducted for the structural demands (θmax and θr,max) and the regression curves can be obtained shown as in Figure 19. The figures also reflect the linear regression results together with 16th percentile and 84th percentile lines. It can be observed from the regression curves that both EDP index θmax and θr,max of the novel SCPC structure nonlinear increased approximately with the seismic intensity of near-fault earthquake waves. However, the dispersion of the peak residual inter-story drift ratio is significantly larger than the peak inter-story drift ratio. The logarithmic standard deviation of the maximum inter-story drift ratio and the maximum residual inter-story drift ratio are 0.3767 and 0.9172. It also presents the mean values of the peak story drift ratio of the novel SCPC frame are 0.34% and 0.83% at the two levels and the mean values of peak residual story drift ratios of the novel SCPC frame are 0.027% and 0.012% at the two seismic levels.
The seismic vulnerability of structures is a characteristic of the structure, which means the possibility of the structure reaching a certain limit state or performance level with the increasing of seismic intensities. The vulnerability analysis usually uses fragility curves to evaluate the seismic performance. The vulnerability function are defined as follows:
F ( x ) = x P [ L S | S I = x ] .
F(x) refers to the seismic vulnerability of the structure, value x refers to the value of certain seismic intensity corresponding to the SI. The structural fragility function can be described by a lognormal model corresponding the limit state shown as follows [46]:
F ( x ) = Φ [ ln x ln m R β d 2 + β c 2 ] ,
where Φ is a standard normal distribution function, m R is median value of the fragility function, β d , β c are assigned to logarithmic standard deviation of the structural demand, capacity respectively. β d can carried out by linear regression of the structural response under different intensities. β c means 0.3 to construct the fragility curves in this paper [47]. The seismic hazard probability model H(x) represents the annual exceeding probability of the certain seismic intensity SI at the certain site, which is related to the magnitude, epicenter distance and the attenuation regulation of ground motion [48]. If the seismic hazard probability is small, the simplified seismic hazard model can be used, the model is expressed as:
H ( x ) = P [ S I > x ] = k 0 x k ,
where k and k0 are shape coefficient of seismic hazard curves, the values of k and k0 are subjected to four parameters—① annual exceeding probability with 10% exceeding probability in 50 years; ② annual exceeding probability with 2% exceeding probability in 50 years; ③ Sa (T1) with 10% exceeding probability in 50 years; ④ Sa (T1) with 2% exceeding probability in 50 years.
Compared with the seismic hazard probability model, the seismic hazard curves of the novel SCPC frame under near-fault ground motions is shown in Figure 20. From seismic hazard curves and Chinese Standard Code [29], Sa (T1) and annual exceeding probability corresponding to the DBE level are 0.1176 and 1/475, respectively, Sa (T1) and annual exceeding probability corresponding to the MCE level are 0.235 and 1/2475, respectively.
Based on the fragility curves and the seismic hazard curves, the seismic risk assessment of a structure can be expressed by annual and 50-year exceeding probability of a certain limit state. The annual probability of exceeding a certain limit state can be calculated by the formula:
P L S = F ( x ) d H ( x ) ,
combined with Formulas (6) and (7), the annual probability can be rewritten as follows:
P L S = ( k 0 m R k ) exp [ ( k β R ) 2 2 ]
the 50-year exceeding probability corresponding to a certain limit state can be written as follows:
P L S / 50 = 1 ( 1 P L S ) 50 .
Based on the theory of fragility analysis and the results of IDA analysis under near-fault earthquake wave, the fragility curves of the novel SCPC structure are plotted in Figure 21. The black and red curves represent the limit state of IO and RE for the engineering demand parameters including θmax and θr,max, respectively. From the fragility curves, the vertical lines represent the corresponding spectral acceleration Sa (T1) under the DBE and MCE levels, which are 0.117 g and 0.234 g calculated with 3% damping ratio of the novel SCPC frame. The vulnerability curves for the seismic limit state from IO to RE become smooth and move significantly toward the right, which means under the same Sa (T1), the exceeding probabilities of the corresponding θmax and θr,max drop obviously. The exceeding probabilities of each performance limit state for the novel frame are also shown in Table 3. At the DBE level, the exceeding probability is 5.5% for IO limit state and 0.58% for RE limit state with θmax as evaluating indicator and the exceeding probability is 1.90% for IO performance limit state and 0.19% for RE limit state with θr,max as evaluating indicator. At the MCE level, the exceeding probability is 36.7% for IO limit state and 2.5% for RE limit state with θmax as evaluating indicator and the exceeding probability is 6.65% for IO performance limit state and 0.45% for RE limit state with θr,max as evaluating indicator. Hence, the novel frame system has good seismic performance and less damage possibility under the IO and RE limit states.
Based on the Chinese Standard Code [29], the DBE and MCE levels have 10% and 2% 50-year exceeding probabilities, respectively. The design objective of the novel SCPC frame can be also summarized as—the 50-year exceeding probability of IO performance limit state is less than 10% and the 50-year exceeding probability of RE performance limit state is less than 2%. Based on the structural fragility analysis and seismic hazard analysis theory, the annual and 50-year exceeding probabilities of different limit states are shown in Figure 22. When the EDP θmax is taken as the evaluating indicator, the annual exceeding probabilities of the novel frame are 0.0628% and 0.0159% for the IO and RE limit state and the 50-year exceeding probabilities of the novel frame are 3.09% and 0.79% for the IO and RE limit state, respectively. When the EDP θr,max is taken as the evaluating indicator, the annual exceeding probabilities of the novel frame are 0.0263% and 0.0045% for the IO and RE limit state, the 50-year exceeding probabilities are 1.31% and 0.23% for the IO and RE limit state, respectively. Therefore, according to the seismic risk assessment of different performance limit states based on structural fragility and seismic hazard analysis theory, the novel SCPC frame can satisfy the seismic limit state and has great structural reliability under the near-fault ground motions.

6. Further Discussion

To verify the great seismic behavior of the novel SCPC frame under near-fault earthquake waves and the accuracy of seismic risk assessment of numerical model, a five-story four-span cast-in-place RC frame (a reinforced concrete frame) and the traditional SCPC frame (a traditional self-centering precast concrete frame without hysteretic dampers) have been modelled by OpenSees. The section sizes and reinforcements arrangement of the two frames are consistent with the prototype structures.
Figure 23a displays the analytical connection model of the RC frame; all the structure members are simulated by the Nonlinear Beam Column element. The Joint2D element is adopted in the modelling of the connection core area and is idealized as a parallelogram shaped shear panel with adjacent elements connected to its mid-points. The Joint2D element includes a spring element inside the joint to represent the shear behavior of the core area and four spring elements located at the beam column interface to simulate the sliding the pulling out of reinforcements. The hysteretic curves of force and deformation of each spring can be simulated by the Pinching4 material model [49]. The Pinching4 material model also reflects the pinching effect, degradation effect of stiffness and strength under cyclic loading. The analytical model of the traditional SCPC frame (shown in Figure 23b) and the modelling method of the novel joint are close to the aforementioned novel SCPC beam-column connection. The concrete is also assigned to the Concrete 01 material model and the enhancement factor of the stirrup restraint effect is 0.2. The steel bars and PT tendons are also assigned to the Steel 02 material model, which are also based on the Giuffre-Menegotto-Pinto model. Because there are no hysteretic dampers to dissipate energy, the two pairs of zeroLength spring elements simulate the gap of connection opening and closing at the end of the beam under earthquake actions. The shear slip at the beam-column interface is avoided when gap has opened. The rigidlink elements are selected to reflect the coupling characteristics between the node of gaps and the relevant nodes of structural components. The Truss element is also used to simulate the unbonded PT tendons between node 7 and node 12 for connection of the traditional SCPC frame.
The site class and seismic intensity are similar to the prototype structures, 21 near-fault earthquake waves from PEER record database are also chosen for the elastoplastic dynamic analysis for the RC frame and the traditional SCPC frame. The peak ground acceleration of selected ground motions are all satisfy the DBE and MCE levels through linearly scaled. The damping ratio of the Rayleigh damper is 5% for the RC structure and 3% for the traditional SCPC structure, respectively. Tenety second last-time free vibration analysis of the frames are also continued after the applying each near-fault ground motion to calculate the residual deformation.
The peak inter-story drift ratios θmax of the three frames for each near-fault ground motion and for each story at the DBE and MCE levels are displayed in Figure 24a,b. Figure 24a indicates that the peak inter-story drift ratios of the traditional SCPC frame are the largest and the novel SCPC frame are the smallest subjected to each earthquake wave at the DBE and MCE levels. The maximum difference between the maximum response of the RC frame and the novel SCPC frame subjected to all the earthquake waves are 11.3% and 12.6% at the DBE and MCE level. The maximum difference between the maximum response of the traditional SCPC frame and the novel SCPC frame subjected to all the earthquake waves is 26.4% and 13.8% at the DBE and MCE levels, respectively. The mean value of the peak story drift ratios of the RC structure subjected to all the earthquake waves are 0.36% and 0.88% subjected to the two levels. And the mean value of peak story drift ratios of the traditional SCPC structure corresponding to all the earthquake waves are 0.40% and 0.91% subjected to the two levels.
Figure 24b shows that the peak inter-story drift ratios of the traditional SCPC structure are the largest and the novel SCPC frame are the smallest for each story under the two levels. The maximum responses of all the stories all happen on the second story for the three structures. The peak responses for the RC frame are equal to 0.72%, 0.54%, 0.31% from the third to top story at the MCE level and the maximum responses for the traditional SCPC frame are 0.76%, 0.58%, 0.33% from the third to fifth story at the MCE level. From the results of the three frames, the extent of gap opening of the novel SCPC frame is small at the DBE level, in addition to the initial stiffness provided by the hysteretic dampers, the peak inter-story drift ratios of the novel frame are relatively small. When the gap opening becomes larger and the bending stiffness significantly decreases under the MCE level, the frame has been in a plastic state. The hysteretic dampers dissipate seismic energy after yielding; the peak responses of the novel frame are still the smallest of the three frames. It contributes the high initial stiffness and better dissipates seismic energy capacity after the yielding provided by the hysteretic dampers.
Figure 25 displays the peak residual inter-story drift ratios θr,max of the three frames for each story at the DBE and MCE levels subjected to the near-fault earthquake waves. The peak residual inter-story drift ratios of the three frames for each story are far below the limiting value of 0.2% and 0.5% at the DBE and MCE levels. The maximum residual responses of the RC frame for each story are significantly larger than the novel SCPC frame and the traditional SCPC frame. The peak residual story drift ratios of the RC frame for all the stories are 1.28 × 10−4 rad and 6.54 × 10−4 rad at the DBE and MCE levels and happen on the first floor. The peak residual story drift ratios of the novel SCPC frame and traditional SCPC frame are similar and close to zero, which attributes the good self-centering capability of the two self-centering precast concrete frames.
The maximum inter-story shear forces of the three frames for each earthquake wave and for each story at the DBE and MCE level are displayed in Figure 26a and Figure 26b, respectively. It shows that the maximum inter-story shear force of the RC structure is the largest of the three frames. The peak inter-story shear force of the RC frame and the traditional SCPC frame are 1034 kN and 1001 kN, corresponding to 21 earthquake waves at the MCE level, respectively. The peak responses of the RC structure and the traditional SCPC structure for all the stories are 710 kN and 687 kN under the MCE level and both happen on the first floor. Due to the different construction mechanisms of the beam-column connections of different concrete frames, horizontal earthquake action will weaken with the increase of natural periods of the novel SCPC frame and the traditional SCPC frame, in addition to the better energy dissipation provided by the hysteretic dampers. Therefore, the inter-story shear force of the novel SCPC frame is the smallest at the MCE level.
The maximum connection rotations and PT tendons forces of the traditional SCPC structure and novel SCPC structure at the MCE level are also drawn in Figure 27. It illustrates—(1) the peak connection rotation of the traditional SCPC frame for each story is obviously larger than that of the novel SCPC frame. The maximum connection rotation of the traditional SCPC frame for all the stories is 3.38 × 10−3 rad and happens on the first story, then decreases gradually. The maximum difference between the maximum rotations of the two SCPC frames for all the stories are 26.7% under the MCE level; (2) The PT tendons forces of the traditional SCPC frame are also significantly larger than the novel frame at the MCE level and the change regulation of PT force in each story of the two frames are relatively homogeneous. The maximum tendon force of the traditional SCPC frame happens on the second floor, which is 0.67 times the yielding tendon force at the MCE level. It reflects that the PT forces of the two SCPC frames are still in the elastic stage and the good energy consumption characteristic of the novel SCPC frame at the MCE level.

7. Conclusions

This study develops a new self-centering precast concrete (SCPC) frame with hysteretic dampers. The configuration and working principles of beam-column connections are illustrated. A prototype structure is modelled and the elastoplastic dynamic responses at the DBE and MCE levels are studied. IDA and seismic fragility analysis are also carried out. In addition, the elastoplastic dynamic responses of the RC structure and the traditional SCPC structure are also studied. From the aforementioned research, the main conclusions are summarized:
  • The novel SCPC frame can satisfy the seismic target under near-fault earthquake waves. The peak and residual inter-story drift ratios all achieve the desirable demands at the DBE and MCE levels. The maximum inter-story shear forces indicate the novel frame can undertake more load due to the large stiffness. The maximum connection relative rotations are minimal and PT forces are all below the yielding tendon force at the MCE level to maintain the self-recovering characteristic. Besides, the beam end M- θ r relation curves of the novel joint are stable and display a typical double-flag shape.
  • The results of IDA indicate the peak and residual inter-story drift ratios both nonlinearly increase approximately with the increasing intensity under near-fault earthquake waves. The dispersion degree of the peak residual inter-story drift ratios is significantly higher than that of the peak inter-story ratios and the dispersion increases with increasing intensity levels.
  • The fragility curves for the seismic limit state from IO to RE become smooth and move significantly toward the right, which means under the same Sa (T1) the exceeding probabilities of the corresponding θmax and θr,max drop obviously. Based on the fragility analysis and seismic hazard analysis, the 50-year exceeding probabilities can be obtained and satisfy the seismic design objective.
  • The peak inter-story drift ratios of the traditional SCPC frame are the largest of the three frames under near-fault earthquake waves loading. The peak residual inter-story drift ratios and the peak inter-story shear forces of the RC frame are the largest; the peak residual inter-story drift ratios of two SCPC frames are minimal. The maximum connection rotations and PT tendons forces of the traditional SCPC frame are significantly larger than those of the novel SCPC frame. The results show that the SCPC frame with hysteretic dampers has great seismic behavior, energy consumption and self-recovering capacities under near-fault earthquake waves loading.

Author Contributions

Original ideas F.G. and Y.D.; data curation, F.G., H.W. and K.Y; funding acquisition, Y.D.; methodology, F.G. and H.W.; writing—original draft preparation F.G.; writing—review and editing, F.G. and K.Y. All authors have read and agreed to the published version of the manuscript.

Funding

The authors gratefully appreciate the support of Program of National Natural Science Foundation (Nos.51608258, 51978154), the Distinguished Young Scientists of Jiangsu Province (No. BK20190013), the Jiangsu Government Scholarship for Overseas Studies (grant number: JS-2018-074), the Program of China Postdoctoral Foundation (No. 2019M651595), Six Peak High Talent Projects of Jiangsu Province (No. JZ-032), the Fundamental Research Fund for Tongji Architectural Design (Group) (No. 2017KY13).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mavroeidis, G.P.; Papageorgiou, A.S. A mathematical representation of near-fault ground motions. Bull. Seismol. Soc. Am. 2003, 93, 1099–1131. [Google Scholar] [CrossRef]
  2. Adanur, S.; Altunişik, A.C.; Bayraktar, A.; Akköse, M. Comparison of near-fault and far-fault ground motion effects on geometrically nonlinear earthquake behavior of suspension bridges. Nat. Hazards 2012, 64, 593–614. [Google Scholar] [CrossRef]
  3. Somerville, P.G. Magnitude scaling of the near fault rupture directivity pulse. Phys. Earth Planet. Inter. 2003, 137, 201–212. [Google Scholar] [CrossRef]
  4. Alavi, B.; Krawinkler, H. Strengthening of moment-resisting frame structures against near-fault ground motion effects. Earthq. Eng. Struct. Dyn. 2004, 33, 707–722. [Google Scholar] [CrossRef]
  5. Yadav, K.K.; Gupta, V.K. Near-fault fling-step ground motions: Characteristics and simulation. Soil Dyn. Earthq. Eng. 2017, 101, 90–104. [Google Scholar] [CrossRef]
  6. Phan, V.; Saiidi, M.S.; Anderson, J.; Ghasemi, H. Near-fault ground motion effects on reinforced concrete bridge columns. J. Struct. Eng. 2007, 133, 982–989. [Google Scholar] [CrossRef]
  7. Takewaki, I.; Tsujimoto, H. Scaling of design earthquake ground motions for tall buildings based on drift and input energy demands. Earthq. Struct. 2011, 2, 171–187. [Google Scholar] [CrossRef]
  8. Chioccarelli, E.; Iervolino, I. Near-source seismic hazard and design scenarios. Earthq. Eng. Struct. Dyn. 2013, 42, 603–622. [Google Scholar] [CrossRef]
  9. Tothong, P.; Cornell, C.A. Structural performance assessment under near-source pulse-like ground motions using advanced ground motion intensity measures. Earthq. Eng. Struct. Dyn. 2008, 37, 1013–1037. [Google Scholar] [CrossRef]
  10. Champion, C.; Liel, A. The effect of near-fault directivity on building seismic collapse risk. Earthq. Eng. Struct. Dyn. 2012, 41, 1391–1409. [Google Scholar] [CrossRef]
  11. Mazza, F.; Vulcano, A. Effects of near-fault ground motions on the nonlinear dynamic response of base-isolated rc framed buildings. Earthq. Eng. Struct. Dyn. 2012, 41, 211–232. [Google Scholar] [CrossRef]
  12. Mortezaei, A.; Ronagh, H.R.; Kheyroddin, A. Seismic evaluation of FRP strengthened RC buildings subjected to near-fault ground motions having fling step. Compos. Struct. 2010, 92, 1200–1211. [Google Scholar] [CrossRef]
  13. Alavi, B.; Krawinkler, H. Behavior of moment-resisting frame structures subjected to near-fault ground motions. Earthq. Eng. Struct. Dyn. 2004, 33, 687–706. [Google Scholar] [CrossRef]
  14. Loh, C.H.; Wan, S.; Liao, W.I. Effects of hysteretic model on seismic demands: Consideration of near-fault ground motions. Struct. Des. Tall Build. 2002, 11, 155–169. [Google Scholar] [CrossRef]
  15. Sehhati, R.; Rodriguez-Marek, A.; ElGawady, M.; Cofer, W.F. Effects of near-fault ground motions and equivalent pulses on multi-story structures. Eng. Struct. 2011, 33, 767–779. [Google Scholar] [CrossRef]
  16. Soleimani Amiri, F.; Ghodrati Amiri, G.; Razeghi, H. Estimation of seismic demands of steel frames subjected to near-fault earthquakes having forward directivity and comparing with pushover analysis results. Struct. Des. Tall Spec. Build. 2013, 22, 975–988. [Google Scholar] [CrossRef]
  17. Ricles, J.M.; Sause, R.; Garlock, M.M.; Zhao, C. Posttensioned seismic-resistant connections for steel frames. J. Struct. Eng. 2001, 127, 113–121. [Google Scholar] [CrossRef]
  18. Rojas, P.; Ricles, J.M.; Sause, R. Seismic performance of post-tensioned steel moment resisting frames with friction devices. J. Struct. Eng. 2005, 131, 529–540. [Google Scholar] [CrossRef]
  19. Lin, Y.C.; Sause, R.; Ricles, J. Seismic performance of a large-scale steel self-centering moment-resisting frame: MCE hybrid simulations and quasi-static pushover tests. J. Struct. Eng. 2013, 139, 1227–1236. [Google Scholar] [CrossRef]
  20. Priestley, M.N.; Tao, J.R. Seismic response of precast prestressed concrete frames with partially debonded tendons. PCI J. 1993, 38, 58–69. [Google Scholar] [CrossRef]
  21. El-Sheikh, M.; Pessiki, S.; Sause, R.; Lu, L.W. Moment rotation behavior of unbonded post-tensioned precast concrete beam-column connections. Struct. J. 2000, 97, 122–131. [Google Scholar]
  22. Cheok, G.S.; Stone, W.C.; Kunnath, S.K. Seismic response of precast concrete frames with hybrid connections. Struct. J. 1998, 95, 527–539. [Google Scholar]
  23. Morgen, B.G.; Kurama, Y.C. Seismic design of friction-damped precast concrete frame structures. J. Struct. Eng. 2007, 133, 1501–1511. [Google Scholar] [CrossRef]
  24. Rodgers, G.W.; Mander, J.B.; Geoffrey Chase, J. Modeling cyclic loading behavior of jointed precast concrete connections including effects of friction, tendon yielding and dampers. Earthq. Eng. Struct. Dyn. 2012, 41, 2215–2233. [Google Scholar] [CrossRef]
  25. Rodgers, G.W.; Solberg, K.M.; Mander, J.B.; Chase, J.G.; Bradley, B.A.; Dhakal, R.P. High-force-to-volume seismic dissipators embedded in a jointed precast concrete frame. J. Struct. Eng. 2012, 138, 375–386. [Google Scholar] [CrossRef] [Green Version]
  26. Lu, X.; Cui, Y.; Liu, J.; Gao, W. Shaking table test and numerical simulation of a 1/2-scale self-centering reinforced concrete frame. Earthq. Eng. Struct. Dyn. 2015, 44, 1899–1917. [Google Scholar] [CrossRef]
  27. Seo, C.Y.; Sause, R. Ductility demands on self-centering systems under earthquake loading. Aci Struct. J. 2005, 102, 275. [Google Scholar]
  28. MOHURD. Code for Design of Concrete Structures; MOHURD: Beijing, China, 2010; GB50010-2010.
  29. MOHURD. Code for Seismic Design of Buildings; MOHURD: Beijing, China, 2010; GB50011-2010.
  30. MOHURD. Load Code for the Design of Building Structures; MOHURD: Beijing, China, 2012; GB50009-2012.
  31. FEMA (Federal Emergency Management Agency). Commentary for the Seismic Rehabilitation of Buildings (FEMA-356); FEMA: Washington, DC, USA, 2000.
  32. Wang, G.Q.; Zhou, X.Y.; Zhang, P.Z.; Igel, H. Characteristics of amplitude and duration for near fault strong ground motion from the 1999 Chi-Chi, Taiwan earthquake. Soil Dyn. Earthq. Eng. 2002, 22, 73–96. [Google Scholar] [CrossRef]
  33. Burks, L.S.; Baker, J.W. A predictive model for fling-step in near-fault ground motions based on recordings and simulations. Soil Dyn. Earthq. Eng. 2016, 80, 119–126. [Google Scholar] [CrossRef]
  34. MOHURD. Specification for Seismic Design of Prestressed Concrete Structures; JGJ 140-2004; MOHURD: Beijing, China, 2012.
  35. Mazzoni, S.; Mckenna, F.; Scott, M.H.; Fenves, G.L. Opensees Command Lang. Munual; Pacific Earthquake Engineering Research (PEER) Center: Berkeley, CA, USA, 2007. [Google Scholar]
  36. Kent, D.C.; Park, R. Flexural members with confined concrete. J. Struct. Div. 1971, 97, 1969–1990. [Google Scholar]
  37. Lai, M.H.; Liang, Y.W.; Wang, Q.; Ren, F.M.; Chen, M.T.; Ho, J.C. A stress-path dependent stress-strain model for FRP-confined concrete. Eng. Struct. 2020, 203, 109824. [Google Scholar] [CrossRef]
  38. Lai, M.H.; Chen, M.T.; Ren, F.M.; Ho, J.C. Uni-axial behaviour of externally confined UHSCFST columns. Thin Walled Struct. 2019, 142, 19–36. [Google Scholar] [CrossRef]
  39. Lai, M.; Hanzic, L.; Ho, J.C. Fillers to improve passing ability of concrete. Struct. Concr. 2019, 20, 185–197. [Google Scholar] [CrossRef] [Green Version]
  40. Lai, M.H.; Ho, J.C. Confinement effect of ring-confined concrete-filled-steel-tube columns under uni-axial load. Eng. Struct. 2014, 67, 123–141. [Google Scholar] [CrossRef]
  41. Menegotto MPinto, P.E. Method of analysis for cyclically loaded reinforced concrete plane frames including changes in geometry and non-elastic behavior of elemets under combined normal force and bending. In Proceedings of the IABSE Symposium, Lisbon, Protugal, 27–29 March 1973; Volume 13, pp. 15–22. [Google Scholar]
  42. Iyama, J.; Seo, C.; Ricles, J.M.; Sause, R. Self-centering moment resisting frames MRFs with bottom flange friction devices under earthquake loading. J. Constr. Steel Res. 2009, 65, 314–325. [Google Scholar] [CrossRef]
  43. Pam, H.J.; Ho, J.C. Length of critical region for confinement steel in limited ductility high-strength reinforced concrete columns. Eng. Struct. 2009, 31, 2896–2908. [Google Scholar] [CrossRef]
  44. Ho, J.C.; Pam, H.J. Inelastic design of low-axially loaded high-strength reinforced concrete columns. Eng. Struct. 2003, 25, 1083–1096. [Google Scholar] [CrossRef]
  45. Feng, D.C.; Wu, G.; Lu, Y. Finite element modelling approach for precast reinforced concrete beam-to-column connections under cyclic loading. Eng. Struct. 2018, 174, 49–66. [Google Scholar] [CrossRef] [Green Version]
  46. Ellingwood, B.R.; Kinali, K. Quantifying and communicating uncertainty in seismic risk assessment. Struct. Saf. 2009, 31, 179–187. [Google Scholar] [CrossRef]
  47. Celik, O.C.; Ellingwood, B.R. Seismic fragilities for non-ductile reinforced concrete frames–Role of aleatoric and epistemic uncertainties. Struct. Saf. 2010, 32, 1–2. [Google Scholar] [CrossRef]
  48. Cornell, C.A.; Jalayer, F.; Hamburger, R.O.; Foutch, D.A. Probabilistic basis for 2000 SAC federal emergency management agency steel moment frame guidelines. J. Struct. Eng. 2002, 128, 526–533. [Google Scholar] [CrossRef] [Green Version]
  49. Celik, O.C.; Ellingwood, B.R. Modeling beam-column joints in fragility assessment of gravity load designed reinforced concrete frames. J. Earthq. Eng. 2008, 12, 357–381. [Google Scholar] [CrossRef]
Figure 1. Details of construction of the novel self-centering precast concrete (SCPC) connection.
Figure 1. Details of construction of the novel self-centering precast concrete (SCPC) connection.
Applsci 10 06510 g001
Figure 2. Details of the hysteretic dampers.
Figure 2. Details of the hysteretic dampers.
Applsci 10 06510 g002
Figure 3. Three-dimensional model of hysteretic dampers and hysteretic curve of the novel joint.
Figure 3. Three-dimensional model of hysteretic dampers and hysteretic curve of the novel joint.
Applsci 10 06510 g003
Figure 4. Cyclic loading test set up.
Figure 4. Cyclic loading test set up.
Applsci 10 06510 g004
Figure 5. Comparisons of hysteretic behavior of numerical simulation and test results.
Figure 5. Comparisons of hysteretic behavior of numerical simulation and test results.
Applsci 10 06510 g005
Figure 6. Theoretical hysteretic model of the novel connection.
Figure 6. Theoretical hysteretic model of the novel connection.
Applsci 10 06510 g006
Figure 7. Theoretical hysteretic model of the novel SCPC connection.
Figure 7. Theoretical hysteretic model of the novel SCPC connection.
Applsci 10 06510 g007
Figure 8. Details dimensions of the prototype structure and steel connector (mm).
Figure 8. Details dimensions of the prototype structure and steel connector (mm).
Applsci 10 06510 g008
Figure 9. The time-histories of near-fault earthquake waves KJM090 and H-E07230.
Figure 9. The time-histories of near-fault earthquake waves KJM090 and H-E07230.
Applsci 10 06510 g009
Figure 10. Acceleration response spectrum of near-fault earthquake waves.
Figure 10. Acceleration response spectrum of near-fault earthquake waves.
Applsci 10 06510 g010
Figure 11. Analytical connection model of the novel S CPC frame with Opensees.
Figure 11. Analytical connection model of the novel S CPC frame with Opensees.
Applsci 10 06510 g011
Figure 12. The cross section of Nonlinear Beam Column elements.
Figure 12. The cross section of Nonlinear Beam Column elements.
Applsci 10 06510 g012
Figure 13. Time-history curves of top floor displacement-time history curves under design basis earthquake (DBE) and maximum considered earthquake (MCE).
Figure 13. Time-history curves of top floor displacement-time history curves under design basis earthquake (DBE) and maximum considered earthquake (MCE).
Applsci 10 06510 g013
Figure 14. Maximum inter-story drift ratio for each near-fault ground motion and for each story.
Figure 14. Maximum inter-story drift ratio for each near-fault ground motion and for each story.
Applsci 10 06510 g014
Figure 15. Maximum residual inter-story drift ratio for each story.
Figure 15. Maximum residual inter-story drift ratio for each story.
Applsci 10 06510 g015
Figure 16. Maximum inter-story shear force for each near-fault ground motion and for each story.
Figure 16. Maximum inter-story shear force for each near-fault ground motion and for each story.
Applsci 10 06510 g016
Figure 17. Maximum connection rotation and standardized PT force for each story.
Figure 17. Maximum connection rotation and standardized PT force for each story.
Applsci 10 06510 g017
Figure 18. M- θ r hysteretic relationships of the beam end joint for the new frame.
Figure 18. M- θ r hysteretic relationships of the beam end joint for the new frame.
Applsci 10 06510 g018
Figure 19. The IDA curves of novel SCPC frame.
Figure 19. The IDA curves of novel SCPC frame.
Applsci 10 06510 g019
Figure 20. Seismic hazard curves of the novel SCPC frame.
Figure 20. Seismic hazard curves of the novel SCPC frame.
Applsci 10 06510 g020
Figure 21. Fragility curves of novel SCPC frames (engineering demand parameters (EDP) of θmax and θr,max).
Figure 21. Fragility curves of novel SCPC frames (engineering demand parameters (EDP) of θmax and θr,max).
Applsci 10 06510 g021
Figure 22. Annual and 50-year exceeding probabilities of novel SCPC frame under near-fault earthquake waves.
Figure 22. Annual and 50-year exceeding probabilities of novel SCPC frame under near-fault earthquake waves.
Applsci 10 06510 g022
Figure 23. Analytical connection model of RC frame and traditional SCPC frame.
Figure 23. Analytical connection model of RC frame and traditional SCPC frame.
Applsci 10 06510 g023
Figure 24. Maximum inter-story drift ratios of three frames at DBE and MCE level.
Figure 24. Maximum inter-story drift ratios of three frames at DBE and MCE level.
Applsci 10 06510 g024
Figure 25. Maximum residual inter-story drift ratios of three frames at DBE and MCE level.
Figure 25. Maximum residual inter-story drift ratios of three frames at DBE and MCE level.
Applsci 10 06510 g025
Figure 26. Maximum inter-story shear force of three frames at MCE level.
Figure 26. Maximum inter-story shear force of three frames at MCE level.
Applsci 10 06510 g026
Figure 27. Maximum connection rotation and standardized PT force of two SCPC frames.
Figure 27. Maximum connection rotation and standardized PT force of two SCPC frames.
Applsci 10 06510 g027
Table 1. Information of selected near-fault ground motions.
Table 1. Information of selected near-fault ground motions.
NumberSeismic NameGround Motion NameMagnitude (Mw)Distance (km)PGA (g)PGV (cm/s)PGD (cm)
1Kocaeli, Turkey 1999GBZ0007.5110.920.24450.342.74
2Kobe 1995TAZ0006.90.270.69368.326.65
3Kobe 1995TAZ0906.90.270.69485.316.75
4Kobe 1995KJM0006.90.960.82181.317.68
5Kobe 1995KJM0906.90.960.59974.319.95
6Loma Prieta 1989LGP0006.933.880.56394.841.18
7Loma Prieta 1989LGP0906.933.880.6055111.5
8Morgan Hill 1984HVR2406.193.480.31239.47.66
9Imperial Valley-06 1979H-E042306.537.050.3676.659.02
10Imperial Valley-06 1982H-E072306.530.560.463109.344.74
11Cape Mendocino 1992CPM0007.016.961.497127.441.01
12Cape Mendocino 1992CPM0907.016.961.0394212.39
13Cape Mendocino 1992PET0007.018.180.5948.421.74
14Cape Mendocino 1992PET0907.018.180.66289.729.55
15Duzce, Turkey 1999DZC2707.146.580.53583.551.59
16Northridge 1994JEN0226.695.430.424106.243.06
17Northridge 1994NWH3606.695.920.5997.238.05
18Northridge 1994RRS2286.696.50.838166.128.7
19Northridge 1994SCS1426.75.350.897102.846.99
20Northridge 1994SPV2706.78.440.75384.818.68
21Northridge 1994SYL3606.75.30.843129.632.68
Table 2. Maximum top floor displacements under 21 ground motions.
Table 2. Maximum top floor displacements under 21 ground motions.
NumberGround Motion NameMax Top Displacement (DBE) (mm)
mm((mm)
Time Occur (s)Max Top Displacement (MCE) (mm)Time Occur (s)Residual Top Displacements (mm)
1GBZ00025.17.247.37.20.23
2TAZ00031.85.970.55.440.40
3TAZ09034.66.3875.15.980.42
4KJM00051.83.38112.43.50.51
5KJM09042.49.8490.29.920.45
6LGP00030.212.2868.311.840.38
7LGP09023.511.154.211.140.23
8HVR24033.513.1273.713.20.42
9H-E0423020.87.846.77.840.15
10H-E0723045.47.3895.37.560.47
11CPM00056.08.08113.28.120.53
12CPM09032.43.669.83.60.41
13PET00036.32.9685.12.980.43
14PET09022.13.0846.23.10.18
15DZC27025.83.7255.13.720.22
16JEN02221.43.645.23.60.16
17NWH36031.36.9257.26.980.39
18RRS22835.37.0681.740.41
19SCS14244.35.8291.25.860.47
20SPV27038.53.585.530.44
21SYL36025.37.355.17.40.23
Table 3. Exceeding probabilities of the two performance limit states.
Table 3. Exceeding probabilities of the two performance limit states.
Earthquake LevelExceeding Probability (%,θmax)Exceeding Probability (%,θr,max)
IOREIORE
DBE level5.500.581.900.19
MCE level36.702.506.650.45

Share and Cite

MDPI and ACS Style

Geng, F.; Ding, Y.; Wu, H.; Yang, K. Seismic Risk Assessment of a Novel Self-Centering Precast Concrete Frame under Near-Fault Ground Motions. Appl. Sci. 2020, 10, 6510. https://doi.org/10.3390/app10186510

AMA Style

Geng F, Ding Y, Wu H, Yang K. Seismic Risk Assessment of a Novel Self-Centering Precast Concrete Frame under Near-Fault Ground Motions. Applied Sciences. 2020; 10(18):6510. https://doi.org/10.3390/app10186510

Chicago/Turabian Style

Geng, Fangfang, Youliang Ding, Honglei Wu, and Kang Yang. 2020. "Seismic Risk Assessment of a Novel Self-Centering Precast Concrete Frame under Near-Fault Ground Motions" Applied Sciences 10, no. 18: 6510. https://doi.org/10.3390/app10186510

APA Style

Geng, F., Ding, Y., Wu, H., & Yang, K. (2020). Seismic Risk Assessment of a Novel Self-Centering Precast Concrete Frame under Near-Fault Ground Motions. Applied Sciences, 10(18), 6510. https://doi.org/10.3390/app10186510

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