Cumulative Damage Evaluation of RC Bridge Piers subjected to Multiple Earthquakes

Successive earthquakes of Kocaeli and Duzce within three months indicated that even the survived lifeline structures such as bridges under the former event may have damage or collapse potential under the latter event due to their possible stiffness degradation. It is thus important that a rigorous seismic analysis of such structures should account for the effect of prior earthquake damage. For this purpose, nonlinear seismic analysis of a reinforced concrete bridge structure has been carried out under both single and multiple earthquake ground motions. Behavior and response evaluation of the bridge piers subjected to such motions have been discussed in terms of using both flexure-axial and flexure-shear-axial interaction models. Analytical results show that the stiffness degradation under multiple earthquake ground motions is more pronounced than that under single earthquake ground motion. In addition, comparison of the response without and with shear demonstrates that shear deformation is of significance. .e response with shear exhibits the increase in displacement demand and decrease in lateral force carrying capacity, leading to a decrease in energy dissipation capacity. It is concluded that seismic analysis of reinforced concrete bridge structure should account for the effect of multiple earthquake ground motions to assess the demand on such structure properly.


Introduction
Most of the progress in earthquake-resistant design has been achieved in terms of the observation of damage inflicted by earthquakes. e perception of the way in which bridge responds to an earthquake was dramatically changed by the damage observed after the 1971 San Fernando earthquake. Bridge damage induced by that earthquake clearly demonstrated that the seismic design provisions at that time were inadequate [1]. Subsequent to the San Fernando earthquake, much research effort has been dedicated to improve the seismic resistance of bridges. Although significant advances have been achieved, bridges, in particular those constructed from reinforced concrete, have not performed as expected.
Whereas most of the damage analyses for such structures have been conducted in terms of the single ground motion of all components, a few studies have been reported in the literature regarding the seismic response analyses of such structures subjected to multiple earthquake ground motions. e earthquake of Kocaeli on 17 August 1999 caused severe damage to many reinforced concrete structures [2]. Following the event (three months later), on the evening of 12 November 1999, another earthquake stroke the area of Duzce and Bolu in the northwestern part of Turkey, which is a nearby region of Kocaeli (Izmit) [3]. is implies that the survived lifeline structures such as bridges under the former event may have damage or collapse potential under the latter event due to their possible stiffness degradation. erefore, it is important that a rigorous seismic analysis of such structures should account for the effect of multiple earthquake ground motions. Work by Aschheim and Black [4] was among those a few studies that took into account the effect of prior earthquake damage. However, this study was limited since the prior damage was modeled as an assumed reduction in initial stiffness.
Studies considering repeated sequence effects were conducted in terms of using single-degree-of-freedom systems [5,6]. Recently, Hatzigeorgiou and Beskos [7] developed a methodology for the quantification of the seismic sequence effect directly onto displacement demands. After a year later, Hatzigeorgiou [8] proposed a new procedure for the ductility control of single-degree-of-freedom systems under repeated near-fault earthquakes. Subsequently, elastic and inelastic spectra were computed for various viscous damping ratios, postyield stiffness ratios, and ductility levels. Yaghmaaeri-Sabegh and Ruiz-Garcia [9] also investigated the ground motion effects of the Varzaghan-Ahar accelerogram sequences recorded during Iran earthquake of 2012 on the nonlinear response of elastic-perfectly plastic singledegree-of-freedom systems.
In addition, multiple earthquake effect was also carried out for reinforced concrete and steel structures. Fragiacomo et al. [10] investigated the damage accumulation produced on multi-degree-of-freedom steel frames by repeated earthquakes, and they compared responses between multidegree-of-freedom and equivalent single-degree-of-freedom systems. Li and Ellingwood [11] evaluated the potentials for aftershocks to cause additional damage to steel moment frame buildings and provided a probabilistic description of structural damage states prior to and following the aftershocks. Ruiz-Garcia and Negrete-Manriquez [12] conducted an analytical study aiming at providing an understanding on the influence of aftershocks on drift demands in regular existing moment-resisting frame buildings. Sarno [13] investigated the effects of multiple earthquakes on reinforced concrete structures. Comprehensive parametric response spectra were evaluated with degrading (stiffness and/or strength) and nondegrading hysteretic models to account for the response of reinforced concrete structures. Raghunandan and Liel [14] quantified the influence of ground motion duration on the predicted collapse response of concrete frame structures. Incremental dynamic analysis was carried out on a set of archetypical reinforced concrete buildings subjected to various ground motions with varying duration. Abdelnaby and Elnashai [15] conducted an assessment of reinforced concrete frames subjected to repeated strong motions in terms of using numerical models that capture the degrading behavior of reinforced concrete frame. Hosseinpour and Abdelnaby [16] also investigated the effects of repeated shaking on the behavior of two eightstory reinforced concrete buildings (both regular and irregular in height) prone to multiple earthquakes.
However, most of the previous studies considering repeated earthquakes are related to single-degree-of-freedom systems and building structures, and very few studies are conducted for reinforced concrete bridge structures. is situation has motivated the present work; the main objective of which is to investigate the effects of multiple earthquake ground motions on the stiffness degradation of reinforced concrete bridges, particularly piers.
In addition, axial forces on bridge piers may vary during earthquake excitations due to vertical ground motion. Vertical ground motion components combined with horizontal components may significantly affect the axial force variation [17]. In particular, taking into account that bridge piers are in general subjected to a high level of axial force due to the weight of the superstructure, the axial force variation can be of significance for such members. Although the axial force variation during earthquake excitations is qualitatively acknowledged, only a few studies include the varying axial force effects on reinforced concrete frame structures [18][19][20].
Galal and Ghobarah [21] suggested an element for the modeling of the biaxial flexural and shear behavior of reinforced concrete columns subjected to various levels of axial forces. Prakash et al. [22] conducted experimental work on reinforced concrete circular columns under various levels of combined cyclic bending, shear, and torsion, together with finite element analyses. However, these studies are somewhat limited to simulate the response of such columns subjected to continually varying axial forces. In addition, Setzler and Sezen [23] developed a model that can predict the monotonic lateral force-displacement relationship for reinforced concrete columns subjected to lateral loading. Subsequently, Sezen and Chowdhury [24] also proposed a macromodel to simulate the cyclic lateral load-deformation response of reinforced concrete columns. Cyclic responses corresponding to flexural, longitudinal bar slip and shear deformations were calculated using the model. However, those models do not likely seem to simulate the interactive behaviour of flexure-shear-axial force variation.
In line with the above, cumulative earthquake damage due to multiple earthquakes was investigated on a reinforced concrete bridge structure. For this purpose, nonlinear inelastic response analysis of a reinforced concrete bridge was carried out under both single and multiple earthquake ground motions. In addition, hysteretic shear-axial interaction representation was employed to simulate flexureshear-axial interaction. Subsequently, behavior and response with and without shear for reinforced concrete bridge piers are discussed.

Analysis Program and Shear Representation
2.1. Analysis Program ZeusNL. ZeusNL [25] has been developed for the nonlinear inelastic analysis of two-and three-dimensional steel, reinforced concrete, and composite structures under static and dynamic loading, taking into account the effects of both geometric and material nonlinearities. ZeusNL [25] provides a series of solution strategies for nonlinear analysis as described below.
In static analysis, both force and displacement loading can be applied with independent values or constrained to vary in proportional ratios. In addition, displacement and acceleration time histories can be applied at the supports. Also availed of is static adaptive strategy where the applied load pattern is varying through the different steps of the procedure.
is allows the stiffness degradation and the period elongation of the structures to be accurately described. e solution procedure can be full or modified Newton-Raphson method. Automatic load-step reduction is employed to provide an optimum efficiency, and convergence criteria can be defined in terms of either displacement or force.
For dynamic analysis, the Lanczos algorithm is used for eigenvalue analysis to obtain the required natural frequencies and mode shapes. Time-history analysis is performed through numerically integrating the equation of motion using either the unconditionally stable Newmark family of algorithm or the Hilber-Hughes-Taylor algorithm [26] within automatic time-step adjustment during analysis. Equilibrium is ensured for each time step using the same iterative strategies as employed in static analysis.
A variety of cross section types are available including steel rectangular solid, hollow, and I-sections, reinforced concrete columns (rectangular and circular) and T-beam sections, and both fully and partially encased composite sections over which a number of monitoring areas are divided, in order to account accurately for the inelastic response of structural members. e layout of steel and both confined and unconfined concrete within the cross section can then be modeled, with the computation of stress resultants being performed automatically.

Hysteretic Shear Representation.
Inelastic deformations generated during seismic response are not limited to flexural deformation. Saatciouglu and Ozcebe [33] indicated from reinforced concrete column tests that inelastic shear deformation can significantly affect total deformation in local areas. Moreover, a greater shear capacity in comparison with that corresponding to flexural capacity does not always ensure elastic shear behavior. erefore, there is a need for a simple yet accurate model for predicting the flexure and shear interactive behavior of reinforced concrete columns, combined with axial force variation.
In view of the above, the hysteretic shear models for variable axial forces were developed and implemented in ZeusNL [25]. e new formulations were composed of two parts, the envelope curve and a set of hysteresis rules. e envelope curve is derived from the modified compression field theory (referred to as MCFT [34] hereafter). e MCFT [34] was widely used for the analysis of reinforced concrete members and gave a good correlation ( [35] amongst many others). However, the stress-strain relationship for concrete in compression used in the MCFT does not seem to provide an effective confinement to the core concrete confined with circular hoops or spirals. us, the stress-strain relationship for concrete in compression is modified to consider the confinement effect [36]. e envelope curve is represented by a quatrilinear symmetric curve as shown in Figure 1(a) which accounts for the cracking, yielding, and ultimate states denoted by "cr," "y," and "ult," respectively. e cracking force is defined as the force at which the principal tensile stress is equal to 0. 33 ��� f co (f co is a maximum stress of concrete), and yield force is defined as the point obtained by extrapolating a straight line from the origin through 0.75 of the ultimate force to meet the horizontal line passing through the ultimate force.
To describe the inelastic hysteretic shear response, a set of hysteresis rules proposed by Ozcebe and Saatcioglu [37] was adopted for the evaluation of loading and unloading (total and partial). However, since most of the expressions derived in the development of the hysteretic shear model were based on the limited experimental data, direct use of the model in inelastic time-history analysis was uncertain.
us, modifications and developments were made with regard to the hysteresis rules [38]. Primary developments are the relationship between axial force and its effect on pinching for various levels of axial forces, and formulation between the number of cycles at a constant level of deformation and its effect on stiffness degradation. Detailed description can be found in references [36,38].
Based on the envelope and a set of shear hysteresis rules, hysteretic shear-axial interaction representation was developed. e basic concept in including the effect of varying axial forces is that the stiffness in the current time step is calculated by introducing appropriate shifts corresponding to the current level of axial force between series of envelope curves derived for constant levels of axial force. is is an equivalent stepwise linear approach. ese transitions represent either hardening (increase in stiffness) or softening (decrease in stiffness) of the member due to variation of axial force. Cracking, yielding, and ultimate levels are also shifted in accordance with the axial force variation. e envelop curves are derived for several different levels of constant axial force, which are defined by the user and incorporate three levels in compressive axial force, zero axial force, and two levels in tensile axial force. Such an arrangement enables the user with a proper selection of axial force range of interest. e envelope curve corresponding to an axial force between these prespecified levels is established using extrapolation. e graphical representation of the shear stiffness transition due to varying axial forces is illustrated in Figure 1(b). As representative cases, the analytical predictions obtained with and without the new formulations are illustrated in Figure 2. As observed, while a good agreement was achieved between experiment and analysis with shear, a significant difference was obtained between experiment and analysis without shear. A more detailed comparison can be found elsewhere [36,38,39].

Analytical Model and Input Ground Motion
In order to investigate the effect of multiple earthquake ground motions on the inelastic response of a bridge, a reinforced concrete bridge structure, severely damaged by the Northridge earthquake of 17 January 1994, is selected. e bridge under consideration is a ramp structure (Collector Distributor 36), which continues on a line close to that of the main freeway, La Cienega-Venice Boulevard sector of the I-10 [40]. e choice of the bridge structure is intuitive. It is appreciated that the chosen bridge structure is not as conclusive as a damaged bridge structure under the two Turkey earthquakes. However, since very few literature studies are available on the damaged reinforced concrete bridge structure under the two successive Turkey earthquakes, the bridge is selected for comparative analyses.  over three single-column bents 6, 7, and 8, and over the pier wall of bent 9. In the deck, a movement joint forming a structural hinge is placed just under 5 m away from the bent 6.

Analytical Model of Santa
e presence of this movement joint may allow both relative displacement and rotation to occur between the ends of the deck. e columns of all bents consist of 1219 mm diameter reinforced concrete circular sections. Column longitudinal reinforcement is identical for piers 6, 7, and 8, while less longitudinal reinforcement is employed in the columns of bent 5. Cross sections of three-celled box girder and cross sections of columns and pier wall are shown in Figures 3(a) and 3(b), respectively. e general layout of the ramp structure is shown in Figure 3(c), which represents the flexure-axial interactive model. In this model, five cubic inelastic elements are employed in the piers, with shorter elements at the base and top of the piers and longer elements toward the center. is arrangement allows plastic hinges to be captured accurately, and hence, the inelastic behavior. To account for the effect of shear, each pier of the ramp structure is modeled by a combination of cubic inelastic elements with a single joint element representing the developed axial-shear interaction at the bottom of the piers.
us, flexure-shear-axial interaction will be simulated through this model. e graphical representation of the model is illustrated in Figure 3(d). e longitudinal reinforcement is continuous to the bottom of the footing, providing moment resistance. e connection between the piers and the deck is also effectively continuous. Detailed description regarding the bridge structure can be found in reference [41].
Use is made of MCFT [34] to define the input parameters of the new shear element. For this purpose, each pier was analyzed for several levels of constant axial force. ese levels are 10%, 20%, and 30% of compressive axial force capacity, zero axial force, and 10% and 30% of the tensile axial force capacity (where both capacity values were calculated according to ACI 318-14 [42]). Since the exceedance of 15% of compressive axial force capacity is not common for reinforced concrete bridge piers [43], the selection of the above axial force range can be considered as an upper bound.

Input Ground Motion.
Two sets of accelerograms (four in total) are selected for comparative dynamic analyses. One set of accelerograms is Kocaeli (K) and Duzce (D) earthquakes, recorded at the same station, Duzce under the two Turkey earthquakes, Kocaeli of August 1999 and Duzce of November 1999, respectively. e other set of accelerograms is San Fernando (F) and Northridge (N) earthquakes, recorded at Pacoima Dam, San Fernando earthquake of 1971 and the City Hall grounds, Northridge earthquake of 1997, respectively. Whereas the former set may represent the effect of damage accumulation on the structure due to their successive nature in the same area, the latter set can be considered as a benchmark since the two stations are not identical although they are geographically proximate. Nonetheless, it is worth noting that the chosen ground motions are not sufficient enough and rather intuitive. However, since the purpose of the present study is to assess the response of the bridge structure experiencing multiple main shocks recorded at almost identical stations with a time interval, the study can be considered as exploratory comparative analyses on the damage accumulation with and without shear. e peak ground accelerations of all components for the four earthquake ground motions are given in Table 1.

Static Analysis for Piers.
Static analysis of each pier was carried out using the MCFT [34] for preliminary evaluation of the pier capacities. Gravity loads assessed from the cross-sectional areas of the box girder and the piers were applied as axial force at the top of each pier. Figure 4 shows the lateral force-displacement response of each pier in the transverse direction. Comparison of static response of piers shows that while the response of piers 6 and 8 is nearly   identical, the response of pier 7 displays much less stiffness than that of piers 6 and 8. It is also noteworthy that the pier 7 experiences the greatest displacement. is can be attributed to different heights between the piers, being the tallest for pier 7. Static response characteristics are summarized in Table 2.
As observed, pier 8 carries a greater axial force due to adjacent long spans. While the yield displacement is identical (32 mm) for piers 6 and 8, that is much greater (62 mm) for pier 7.         show comparisons of displacement ductility demand for piers in the transverse direction, subjected to single (D) and multiple earthquakes (K + D) under biaxial and triaxial input ground motions, respectively. In general, the analyses under multiple earthquakes show higher displacement ductility demand than those under single earthquake, being approximately 10% increase. Also depicted in the figure are comparisons between with and without shear. As observed, the analyses with shear exhibit higher displacement ductility demand than those without shear, up to approximately 30% increase for pier 6. Shown in Figures 6(a) and 6(b) are comparisons of displacement ductility demand for piers in the transverse direction, subjected to single (N) and multiple earthquakes (F + N) under biaxial and triaxial input ground motions, respectively. Whereas the same trend of displacement ductility demand increase is observed in Figure 5, the discrepancy between single and multiple earthquakes is rather significant, being approximately 50% increased ductility demand for both with and without shear. It is therefore noteworthy that displacement response can be significantly affected by the applied ground motion characteristics. However, a minor effect of vertical ground motion on the ductility demand is observed from the figures. e maximum displacement in the transverse direction is summarized in Table 3. As observed in general, the maximum displacement with shear is greater than that without shear.

Nonlinear Inelastic
Further investigation was carried out for the effect of multiple earthquakes in terms of the hysteretic response of piers in the transverse direction. Figures 7 and 8 show comparisons of transverse hysteretic response of piers for flexure-axial only model cases between single and multiple earthquakes of biaxial input ground motions, subjected to a series of Turkey and USA earthquakes, respectively. For both cases in general, the response under multiple earthquake motions exhibits more pronounced increased inelastic displacement than that under single earthquake motion. In addition, the former experiences a number of inelastic cycles as well as low cycles. is implies that the pier stiffness will decrease with the number of cycles experienced, even if the maximum ductility discrepancy between single and multiple earthquake motions is small (10% difference in the former case), as shown in Figure 5. In other words, the damage to the structure subjected to multiple earthquake motions is more than that under single earthquake motion. Hence, the best measure of damage evaluation should include both the maximum displacement experienced and the stiffness due to the number of cycles. is also stresses the importance of response analysis under multiple earthquake motions.
In order to investigate the effect of multiple earthquakes on the response with shear, comparisons of hysteretic response are made for piers 6 and 8 of the flexure-shear-axial interaction model. Figures 9 and 10 show the transverse hysteretic response for piers 6 and 8 of the flexure-shearaxial interaction model cases between single and multiple earthquakes of triaxial input ground motions, subjected to series of Turkey and USA earthquakes, respectively.
As observed, the response under multiple earthquakes experiences a greater displacement and shows more pronounced stiffness degradation. It is important to note from the response with shear under multiple earthquakes that once single cycle of maximum displacement occurs, the reloading stiffness is significantly reduced in comparison with the stiffness of the first cycle.
is indicates that multiple earthquake motions combined with shear affect the energy dissipation capacities of piers, leading to the reduction of stiffness of piers and hence the damage mode.
Comparisons of the response with and without shear are made to investigate the effect of shear on the damage mode of the bridge structure subjected to multiple earthquakes. Figures 11 and 12 illustrate the transverse displacement response without and with shear subjected to biaxial input ground motions of multiple earthquakes, K + D and F + N, respectively. e increase in total displacement with the inclusion of shear is clearly observed in the figures. is indicates that piers are subjected to higher displacement demands due to shear (approximately 20% higher in general), particularly for piers 6 and 8. is may be attributed to the irregular bridge structure under consideration, where piers have different heights, and thus, the ductility demand may be concentrated on a particular member. e breakdown of total displacement into its components shows that shear displacement is pronounced. e contribution of shear displacement to total displacement reaches a significant level, approximately 45% and 40% for piers 6 and 8, and a lower value of 30% for pier 7. e maximum transverse response with reference to displacement components is summarized in Table 4. Figures 13 and 14 show comparisons of the corresponding transverse lateral force response without and with shear subjected to biaxial input ground motions of multiple earthquakes, K + D and F + N, respectively. Reduction in lateral force carrying capacity for the response with shear is obvious for piers 6 and 8, although vibration cycles are nearly identical. is explains that those piers are expected to experience significant stiffness degradation due to both increase in displacement and decrease in lateral force carrying capacity. is is supported further by the evaluation of hysteretic response. Figure 15 illustrates the transverse hysteretic response for piers 6 and 8 subjected to biaxial input ground motions of multiple earthquakes, F + N. While the response without shear shows relatively stable hysteresis loops, that with shear exhibits severe stiffness degradation, increased displacement, and sudden drop in lateral force carrying capacity. In  addition, the amount of energy absorbed and dissipated by piers can be significantly affected by the presence of shear. Particular emphasis is placed on the reloading stiffness, once the first cycle of maximum displacement is attained.
Significant reduction in stiffness (or pinching) is pronounced for the response with shear. is also suggests that shear under multiple earthquake motions affects the damage mode of piers. Table 5 shows the maximum transverse lateral   force response without and with shear. Considerable reduction in lateral force carrying capacity for the response with shear is clearly demonstrated.

Conclusions
Analytical study has been undertaken for the nonlinear inelastic seismic response of a reinforced concrete bridge under single and multiple earthquake ground motions. e analytical results show that the maximum displacement ductility demand imposed on the bridge piers depends primarily on the applied ground motion characteristics. Although the ductility demand under multiple earthquake motions is greater than that under single motion, the difference between the two cases is not significantly considerable. However, the stiffness degradation under multiple earthquake ground motions is more pronounced than that under single earthquake motion. is is supported by the hysteretic response of piers. e response shows that once the first cycle of maximum displacement is attained, loading stiffness of the second cycle under multiple earthquake ground motions is significantly reduced, in conjunction with a number of inelastic cycles, leading to experience more damage. e effect of multiple earthquake ground motions on the response with shear is also investigated. Comparisons of the response without and with shear demonstrate that shear deformation is important. e response with shear exhibits the increase in displacement demand, decrease in lateral  force carrying capacity, and energy absorption and dissipation capacity of piers. Moreover, shear deformation reaches a significant level (up to 45% to total deformation in this study). In short, the multiple earthquake ground motions can considerably affect both flexural and shear stiffness degradation and deformation capacity in conjunction with shear. In addition, the response of bridge piers can be affected by the applied earthquake input motion sequence. us, the seismic analysis and design of bridges piers should account for the effect of multiple earthquake ground motions in order to assess the demand on such members properly.

Data Availability
is study is a numerical study, and the numerical study results may be sent upon requests.

Disclosure
Parts of figures and analytical results were employed from the paper by the first author, Do Hyung Lee. e paper is "Damage potential of RC bridge piers due to multiple earthquakes," 13 th World Conference on Earthquake Engineering, Paper no. 1020, Vancouver, B.C., Canada, 2004.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.