Analysis of the Effect of Mainshock-Aftershock Sequences on the Fragility of RC Bridge Columns

: The mainshock (MS) is often accompanied by a number of aftershocks (AS). The existence of AS may cause the seismic demand to be greater than the MS. In order to better evaluate the impact of AS, this paper takes RC columns as the research object and performs incremental dynamic analysis (IDA) on the actual recorded mainshock-aftershocks (MS-AS). The Park–Ang model and incremental damage index are used to quantify the effect of the MS and AS, respectively. The damage and fragility analysis of the parameters such as reinforcement ratio, axial compression ratio and shear-span ratio are carried out respectively. The results show that the seismic demand of the MS-AS is greater than the MS. Besides, the damage of the column gradually increases with the increase of axial compression ratio and shear-span ratio, and gradually decreases with the increase of the reinforcement ratio. When the seismic design grade is 7, 8, and 9 degree, the maximum increase rate of additional damage caused by aftershocks is 7, 13, and 15% of the MS, respectively. When the column is in a medium damaged and a severely damaged state, the growth rate of additional damage can be estimated to be 12.7 and 11% of the MS, respectively. The fragility of columns in different damage states under the action of MS-AS is greater than that of MS. Reducing the axial pressure ratio can greatly reduce the damage probability of columns in different damage states. The effect of the MS-AS can be comprehensively considered to select appropriate design parameters in the design, and the additional damage caused by the AS can be estimated according to the damage condition of the column.


Introduction
A large amount of historical data indicates that earthquakes often occur in the form of earthquake sequences [1].Although most aftershocks are not as strong as the mainshock, the strength and stiffness of the concrete structure are often degraded due to the effect of the mainshock.Due to the limited time between the mainshock and the aftershock, it is often impossible to repair, and the effect of the aftershock on the structure is even greater than the mainshock, causing further structural damage or even collapse [2,3].The past series of earthquakes also show that the existence of aftershock tends to exacerbate the loss of structures in earthquakes.Quantifying the risk of further damage to the bridge from aftershocks is essential to post mainshock decision-making, functionality, and recovery.The compounding effect of the damage and disruption caused by the earthquake sequences in Chi-Chi (1999), Wenchuan (2008), Tohoku (2011), and Italy (2016) are just a few examples of the implications of mainshock-aftershock event sequences (Atzori et al. 2012 [4]; Kazama and Noda 2012 [5]; Galadini et al., 2017 [6]).Therefore, the damage of aftershock to the structure cannot be ignored.However, the magnitude of aftershock effect on concrete structures and the effect of aftershock in different damage states are not clear.Even the impact of aftershocks is not considered in the existing specification.
In recent years, the research on aftershocks has been beneficially explored by the academic community.Ji et al. [7] selected more than 200 MS-AS sequences from very dense and hard soil profiles, performed a nonlinear response history analysis of the single-degreeof-freedom, and evaluated the effects of structure, strong ground motion, and site-specific parameters on the γ index.Wen et al. [8] studied the effect of rotating mainshock-aftershock sequences in different directions on the seismic performance of the structure.The effect of ground motion rotation is significant for MS-AS sequences and can exceed 25% for all EDPs considered.Amiri S et al. [9] performed MS-AS analysis for a single degree of freedom using residual displacement as an index.The functional relationship between the residual displacement ratio and the elastic vibration period is predicted under different post-yielding stiffness ratio levels and different aftershock intensities.The proposed equations can be used for structural seismic evaluation for MS-AS sequences.Shokrabadi M et al. [10,11] used the MS-AS to evaluate the structural performance of five ductile reinforced concrete frames of different heights using a continuous nonlinear response history analysis method.Zhang et al. [12] proposed a new approach to develop state-dependent fragility curves using real MS-AS records.Results indicated that fragility curves based on residual interstorey drift ratio, peak ground velocity, and maximum inter-storey drift ratio are the best choice to characterize the cumulative damage effect.A summary of the above studies reveals that most of the studies are based on real main aftershock ground vibrations for singledegree-of-freedom structures.However, most of the studies are limited to structures of a particular period, and less research has been conducted on the additional damage caused by aftershocks in different damage states.Therefore, the additional damage caused by aftershocks in different damage states needs to be further explored.
A fragility curve describes the probability of a structure being damaged beyond a specific limit state for various levels of earthquake.An important aspect of structural seismic fragility is to derive the fragility curve to determine the degree of damage under different intensities [13].Chen et al. [14] investigates the seismic fragility and the associated direct financial losses of tall pier bridges subjected to MA sequences using the vector-based intensity measure considering both main-and aftershocks.The analytical results show that while the structural vulnerabilities increase with both the intensities of main-and aftershocks, the role of mainshocks is observed more substantially.Hashemi S V et al. [15] studied the seismic performance of the BRB frame and the effect of SMA on the seismic performance of the frame.The results showed that the use of SMA can reduce the cost of repair and restoration of damaged systems, making building systems more resilient.Ali Massumi et al. [16] presented the effects of aftershock ground motion on the collapse capacity of post-mainshock buildings.The different parameters were selected to evaluate the effect of its characteristics on the collapse capacity of buildings.Irfan Z et al. [17] assessed the seismic performance of existing buildings through an IDA-based vulnerability curve.The type of building damage can be predicted by referring to the vulnerability curve.However, previous studies on bridge seismic risk assessment were less focused on the effect of the mainshocks, and aftershock has often been ignored [18][19][20][21][22][23].Neglecting aftershocks may lead to a severe underestimation of the seismic collapse risk [24][25][26][27][28].An accurate evaluation of highway bridges under the action of seismic sequence is necessary.A summary of the above on vulnerability can be found: IDA is an effective means to study the main-aftershock.Most of the studies of fragility to damage are directed to a specific structure.For bridge structures, the increments of aftershocks on the probability of exceeding different damage states need to be further explored.
Based on the shortcomings of the above studies, the purpose of this paper is to quantify the damage of aftershocks at different earthquake magnitudes and the effects of aftershocks on different damaged structures.Therefore, the study combines the basic structure of the bridge, i.e., the column, to carry out the time history analysis of the MS-AS.The difference in the period of bridge piers can be realized by adjusting the reinforcement ratio and the axial compression ratio.The models of the RC column are established with OpenSees.The nonlinear time history analysis of incremental dynamic analysis (IDA) under multiple working conditions of MS and MS-AS are carried out based on the actually recorded seismic waves.The damage of columns with different parameters is quantified based on the Park-Ang damage model.The additional damage caused by AS under different seismic intensity and the effect of AS under different damage states are obtained.Combined with the theory of vulnerability analysis, the vulnerability curves are used to comprehensively evaluate the damage probability of columns in different damage states.This paper provides guidance for the design of bridges in multi-seismic areas and damaged bridges to resist AS.The flow chart of research content and significance is shown in Figure 1.
shocks on different damaged structures.Therefore, the study combines the basic structure of the bridge, i.e., the column, to carry out the time history analysis of the MS-AS.The difference in the period of bridge piers can be realized by adjusting the reinforcement ratio and the axial compression ratio.The models of the RC column are established with Open-Sees.The nonlinear time history analysis of incremental dynamic analysis (IDA) under multiple working conditions of MS and MS-AS are carried out based on the actually recorded seismic waves.The damage of columns with different parameters is quantified based on the Park-Ang damage model.The additional damage caused by AS under different seismic intensity and the effect of AS under different damage states are obtained.Combined with the theory of vulnerability analysis, the vulnerability curves are used to comprehensively evaluate the damage probability of columns in different damage states.This paper provides guidance for the design of bridges in multi-seismic areas and damaged bridges to resist AS.The flow chart of research content and significance is shown in Figure 1.As can be seen from Figure 1, the seismic magnitude is selected according to the site conditions for site-specific design.Based on the quantitative description of aftershock deterioration damage in this paper, it is verified whether the bridge meets the design requirements under the action of the MS-AS.For the damaged piers, the magnitude of the AS impact on the damaged structure is judged according to the existing damage state so that some necessary measures can be taken to avoid further economic losses.

Damage Model and Damage Index
The Park-Ang damage model was proposed by Park and Ang [29] in 1985.The model used the linear combination of the maximum displacement of the structure and the accumulated hysteretic energy to express the damage of the structure.The calculation formula of the Park-Ang is shown in Equation (1).As can be seen from Figure 1, the seismic magnitude is selected according to the site conditions for site-specific design.Based on the quantitative description of aftershock deterioration damage in this paper, it is verified whether the bridge meets the design requirements under the action of the MS-AS.For the damaged piers, the magnitude of the AS impact on the damaged structure is judged according to the existing damage state so that some necessary measures can be taken to avoid further economic losses.

Damage Model and Damage Index
The Park-Ang damage model was proposed by Park and Ang [29] in 1985.The model used the linear combination of the maximum displacement of the structure and the accumulated hysteretic energy to express the damage of the structure.The calculation formula of the Park-Ang is shown in Equation (1).
In Equation (1): DI is the damage index of the structure; δ max is the maximum displacement of the component under the action of an earthquake; δ u is the ultimate displacement of a component under monotonic loading failure.The deformation corresponding to the strength of the component drops to 85%; F y is the yield shear force of the component; dE is the hysteretic energy absorbed by the component; β is the energy dissipation factor and the value is 0.15 in this study.
The corresponding relationship between the different damage states and the damage index obtained from the Kunnath [30] column test.The research of Hose et al. [31] on the performance level of bridge structures are considered in this study.This paper uses five kinds of levels of damage.The corresponding relationship between different grades and damage index is shown in Table 1.

Medium damage
The cracks developed significantly and the concrete protective layer began to fall off 0.25~0.4

Serious injury
The crack widens sharply, and the partial concrete protective layer falls off 0.4~0.8

Completely destroyed
Stirrup fracture or longitudinal reinforcement Buckling fracture, the core area concrete is crushed ≥0.8

Incremental Damage Evaluation Index
In order to explain the relationship between the additional damage caused by the AS and the damage caused by the MS, Yu [32] proposed the incremental damage index to quantify the damage of AS.The expression is shown in Equation (2).
In Equation (2): δ DI is incremental damage index; DI MA and DI M are the cumulative damage of the structure under the action of the MS-AS and the MS, respectively.

Fragility Analysis Theory
The fragility of structure refers to the probability of exceeding a certain failure state under the action of a certain intensity of ground motion.The expression of beyond probability is shown in Equation (3): In Equation (3): D is the seismic demand of the structure; C is the seismic capacity of the structure; I M is the seismic intensity parameter.
The nonlinear time history analysis of the structure is carried out.The logarithmic discrete points of the seismic intensity parameters and the seismic demand values of the structure are input into the Cartesian coordinate system.The logarithmic discrete point regression analysis is carried out to obtain the relationship between the seismic demand of the structure and the seismic intensity [33,34].The expression is shown in Equation (4): In Equation ( 4): I M represents the ground motion intensity parameter; a, b represent the linear regression coefficient.
The logarithmic standard deviation β D|I M is expressed as Equation (5): In Equation (5), N is the number of nonlinear time history analysis, Sr is the residual sum of squares of regression analysis.
The expression is shown in Equation ( 6): Combining Equations ( 4) and (5) to obtain a probabilistic seismic demand model, as shown in Equation (7): Many scholars had assumed that both the capacity limit state and the seismic response of the component obey the lognormal distribution, and the logarithmic value also obeys the lognormal distribution [35].The lognormal distribution can be expressed as Equation (8): In Equation ( 8): ) In Equation ( 9), σ lnx is the logarithmic standard deviation.In Equation (10), µ lnx is the log mean of the structure response.Equation ( 7) can be expressed as the Equation ( 11): In Equation (11), S c is the mean value of the seismic capacity of the structure, β c and β D are the standard deviations of competence and demand for structure.When the spectral acceleration Sa [36] is used as the independent variable, β 2 D + β 2 C is taken as 0.4; when the peak acceleration PGA is used as the independent variable, β 2 D + β 2 C is taken as 0.5.

Numerical Simulation of Bridge Column
According to the regulation of the bridge in the current "Guideline for Seismic Design of Highway Bridges", the bridge column can usually be simplified as "single column model" in seismic analysis.In this study, the time history analysis of MS-AS is carried out.When the shear-span ratio of the column (L 0 /D) is 10, the axial compression ratio µ is from 0.05 to 0.2.The axial compression ratio in this paper is calculated by the standard value of concrete strength, i.e., f ck = 26.8Mpa.When the axial compression ratio µ is 0.1, the shear span ratio of the column is 6, 10, 14 and 18, respectively.Based on the OpenSees finite element software, different materials are used to simulate the different positions of the section, i.e., the Concrete01 is used to simulate the protective layer and constrained concrete, and the Kent-Park model [37] is used to calculate the stress-strain relationship.The reinforcing steel constitutive model [38] can well consider the yield and fatigue of longitudinal reinforcement and reflect the yield and hardening of reinforcement under earthquake action, which is used to simulate reinforcement.Each node of the column is connected by fiber elements and the cross-section of the column is divided into several fibers, and different fiber materials give different constitutive relationships.The element type of the node connection is the nonlinear beam column.Element deformation takes into account the P-∆ effect.The simplified numerical model of the column is shown in Figure 2. The basic parameters of reinforcement and concrete are shown in Table 2.
tudinal reinforcement and reflect the yield and hardening of reinforcement under earthquake action, which is used to simulate reinforcement.Each node of the column is connected by fiber elements and the cross-section of the column is divided into several fibers, and different fiber materials give different constitutive relationships.The element type of the node connection is the nonlinear beam column.Element deformation takes into account the P-Δ effect.The simplified numerical model of the column is shown in Figure 2. The basic parameters of reinforcement and concrete are shown in Table 2.

Determination of Basic Parameters of Damage Model
In order to obtain the required parameters of the damage model, the piers are subjected to pushover loading.The yield displacement and yield strength of the piers are obtained according to the equivalent energy method.The results are shown in Figure 3. D Y represents the yield displacement; D U represents the ultimate displacement, and F Y represents the yield load. µ represents the axial pressure ratio in Figure 3. Different structures often have different periods, and the basic period of the structure represents the basic characteristics of the structure.Therefore, a basic periodic analysis is performed on the designed piers.The results are shown in Table 3.
jected to pushover loading.The yield displacement and yield strength of the piers are obtained according to the equivalent energy method.The results are shown in Figure 3. DY represents the yield displacement; DU represents the ultimate displacement, and FY represents the yield load. μ represents the axial pressure ratio in Figure 3. Different structures often have different periods, and the basic period of the structure represents the basic characteristics of the structure.Therefore, a basic periodic analysis is performed on the designed piers.The results are shown in Table 3.  Figure 3 indicates that the equivalent yield displacement of the column is almost unchanged with the increase of the reinforcement ratio and axial compression ratio.The ultimate displacement decreases nonlinearly with the increase of the axial compression ratio.The larger the axial compression ratio, the smaller the ultimate displacement.The equivalent yield load of the column increases linearly with the increase of the  Figure 3 indicates that the equivalent yield displacement of the column is almost unchanged with the increase of the reinforcement ratio and axial compression ratio.The ultimate displacement decreases nonlinearly with the increase of the axial compression ratio.The larger the axial compression ratio, the smaller the ultimate displacement.The equivalent yield load of the column increases linearly with the increase of the reinforcement ratio.The larger the axial compression ratio, the greater the equivalent yield load.Table 3 shows that the basic period of the column gradually decreases with the increase of the reinforcement ratio and increases with the increase of the axial compression ratio.
When the axial compression ratio increases, the concrete area in the compression zone increases, and the stress of the tensile steel bars decreases, which can result in a slight increase in the equivalent yield load of the column.The ultimate displacement is the displacement value when the column drops to 85% of the maximum load.Due to the P-∆ effect of the column and eccentric compression, the concrete at the bottom of the column with a larger axial compression ratio is more likely to collapse and reach the ultimate displacement.When the reinforcement ratio increases, the height of the concrete compression zone increases and most of the column section is compressed, and a small eccentric compression failure mode appears.As a result, the concrete in the compression zone is more likely to reach the ultimate compressive stress and reach the ultimate displacement.

Selection and Synthesis of Seismic Waves
Since multiple AS often occur after the MS, this paper selected the largest AS as the AS [39,40].Due to the seismic frequency spectrum characteristics and the effect of duration, the MS-AS seismic waves of different sites are different.According to the shear wave velocity, 26 mainshocks and aftershocks with shear wave velocities between 350 and 500 m/s in the Chi-Chi earthquake were selected.The acceleration response spectra of the MS and AS with damping ratio of 5% are shown in Figures 4 and 5, respectively.The intensity of ground motion [41] can be described by the peak acceleration (PGA), the response spectrum acceleration Sa (T), etc.The research results of Padgett et al. [42] showed that the PGA was a more appropriate index of earthquake intensity.Therefore, this study chooses PGA as the seismic wave intensity index.The MS and AS measured at the same station are connected to form the MS-AS sequences.The MS and AS are separated by 30 s to ensure that the structure has enough free vibration time after MS to simulate real MS-AS sequences [43].
with a larger axial compression ratio is more likely to collapse and reach the ultimate displacement.When the reinforcement ratio increases, the height of the concrete compression zone increases and most of the column section is compressed, and a small eccentric compression failure mode appears.As a result, the concrete in the compression zone is more likely to reach the ultimate compressive stress and reach the ultimate displacement.

Selection and Synthesis of Seismic Waves
Since multiple AS often occur after the MS, this paper selected the largest AS as the AS [39,40].Due to the seismic frequency spectrum characteristics and the effect of duration, the MS-AS seismic waves of different sites are different.According to the shear wave velocity, 26 mainshocks and aftershocks with shear wave velocities between 350 and 500 m/s in the Chi-Chi earthquake were selected.The acceleration response spectra of the MS and AS with damping ratio of 5% are shown in Figure 4 and Figure 5, respectively.The intensity of ground motion [41] can be described by the peak acceleration (PGA), the response spectrum acceleration Sa (T), etc.The research results of Padgett et al. [42] showed that the PGA was a more appropriate index of earthquake intensity.Therefore, this study chooses PGA as the seismic wave intensity index.The MS and AS measured at the same station are connected to form the MS-AS sequences.The MS and AS are separated by 30 s to ensure that the structure has enough free vibration time after MS to simulate real MS-AS sequences [43].

Damage Analysis of RC Column under MS-AS
The Guideline for Seismic Design of Highway Bridges (JTGTB02-01-2008) gives different seismic acceleration peaks to different seismic intensity, i.e., 6 degree (a = 0.05 g), 7 degree

Damage Analysis of RC Column under MS-AS
The Guideline for Seismic Design of Highway Bridges (JTGTB02-01-2008) gives different seismic acceleration peaks to different seismic intensity, i.e., 6 degree (a = 0.05 g), 7 degree (a = 0.1 g or 0.15 g), 8 degree (a = 0.2 g or 0.3 g), 9 degree (a = 0.4 g).This study modulates the selected peak acceleration of ground motion between 0.05 g and 0.4 g at intervals of 0.05 g.A total of 12,896 working conditions is obtained.The maximum displacement at the top of the column and the shear force at the bottom of the column are extracted from the analysis results, and the damage is calculated according to the Park-Ang model.The average value of the damage index under the same peak acceleration is taken as the damage of a certain earthquake intensity, and the effect of aftershock is explained by the incremental damage index.

Effect of Axial Compression Ratio on Damage of RC Column
In order to illustrate the influence of aftershocks on bridge piers with different periods, the cross-combination is carried out in combination with the common reinforcement ratio and axial compression ratio.The MS and MS-AS sequences time history analysis were performed on the designed columns.Due to the limitation of space, only the damage index of partial reinforcement ratio is listed here.The damage of the columns is shown in Figure 6a-d  Figure 6 indicates that the seismic demand of the column under the action of the MS-AS is greater than MS, and this demand increases with the increase of axial compression ratio and decreases with the increase of reinforcement ratio.Besides, the column is in a slight or no damage state when PGA ≤ 0.1 g, and the effect of aftershock can be basically Figure 6 indicates that the seismic demand of the column under the action of the MS-AS is greater than MS, and this demand increases with the increase of axial compression ratio and decreases with the increase of reinforcement ratio.Besides, the column is in a slight or no damage state when PGA ≤ 0.1 g, and the effect of aftershock can be basically ignored.However, the additional damage caused by the aftershock is less when the column damage index is greater than 0.8.The reason is that the average value of the column damage under the same PGA is taken.When the damage of the mainshock is greater than 1, aftershock is not probable to cause additional damage.The number of seismic waves with damage index exceeding 1 caused by mainshock increases with the increase of mainshock damage, resulting in a small proportion of the aftershock damage.
Figure 6a shows that the effect of aftershock under different earthquake intensities decreases with the reduction of mainshock damage.The damage of different columns under the action of the MS-AS is reduced from slight damage state to non-damage state when PGA = 0.2 g, which could possibly improve the durability of the column.Furthermore, increasing the reinforcement ratio can not only significantly reduce the damage of the column caused by the MS, but also reduce the effect of aftershock when PGA ≥ 0.3 g. Figure 6b indicates that the effect of aftershock under different earthquake intensities decreases with the decrease of mainshock damage.Besides, increasing the reinforcement ratio is likely to reduce the damage of the column from medium damage to slight damage when PGA = 0.2 g.Increasing the reinforcement ratio is possible to reduce the column from severe damage to medium damage and avoid further deterioration when PGA = 0.3 g.Moreover, some columns with smaller reinforcement ratio can have collapsed due to the effect of aftershock when PGA = 0.4 g and increasing the reinforcement ratio is likely to avoid the collapse of the bridge.
Figure 6c indicates that increasing the reinforcement ratio can reduce the damage of the column from severe damage to slightly damaged when PGA = 0.2 g, and the columns generally do not collapse.Besides, increasing the reinforcement ratio is possible to reduce the damage index when PGA = 0.3 g, and the column may collapse due to the impact of aftershocks when PGA = 0.4 g and ρ = 3.0%.Figure 6d shows that increasing the reinforcement ratio can reduce the damage of the column from severe damage to medium damage and avoid further deterioration when PGA = 0.2 g.Besides, the column with smaller reinforcement ratio can collapse due to the effect of aftershock, and increasing the reinforcement ratio is probable to prevent the bridge from collapsing when PGA = 0.3 g.Furthermore, it is probable that the columns collapse when PGA = 0.4 g.The collapse of the column is likely to be avoided by adjusting the axial compression ratio.Therefore, the appropriate axial compression ratio and reinforcement ratio should be selected for the high earthquake area.The effect of aftershock should be considered in varying degrees.
In order to study the effect of aftershock, the above Formula ( 2) is used to quantitatively analyze the additional damage caused by aftershock.For the convenience of statistics, the growth rate of the damage of the column caused by aftershock is calculated separately.The maximum damage rate of columns with different reinforcement ratios is regarded as the damage of a certain axial compression ratio.The results are shown in Table 4. Table 4 indicates that the growth rate of additional damage to the columns caused by an aftershock gradually decreases with the increase of the axial compression ratio, and it initially increases and then decreases with the increase of PGA.Besides, the maximum damage growth rate of columns of different magnitudes is not probable to exceed 13%.When the seismic level is 7 degree, the maximum increase rate of damage caused by aftershock is about 7% of the mainshock.When the seismic level is 8 degree, the maximum increase rate of damage caused by aftershock is about 13% of the mainshock.When the seismic level is 9 degree, the maximum growth rate of damage caused by aftershock is about 12% of the mainshock.Due to the damage of the column, a slight aftershock is likely to exceed the bearing capacity limit, causing serious economic losses.Therefore, the additional damage caused by aftershock can be comprehensively considered in the design according to different magnitudes.
The growth rate of additional damage to the column caused by aftershock does not appear at µ = 0.2.The main reasons are: (1) As the axial compression ratio increases, the damage to the column caused by the mainshock gradually increases, and the additional caused by aftershock is relatively reduced.(2) When the damage index of the mainshock is greater than 1, it is probable that the aftershock is not likely to cause additional damage, resulting in a relatively low increase in damage caused by the aftershock.

Effect of Shear-Span Ratio on Damage of RC Column
In order to illustrate the generality of the above research conclusions, a comparative analysis of MS-AS damage of bridge columns with different shear-span ratios is carried out.The basic parameters of the bridge column are ρ = 2.0%, ρ sv = 1.5%, µ = 0.1.Time history and modal analysis of the column with shear-span ratios of 6, 10, 14, and 18 are carried out.Figure 7 shows the damage comparison chart of the bridge column under the action of MS and MS-AS of different seismic intensity.The basic period and additional damage are shown in Table 5.

Effect of Shear-Span Ratio on Damage of RC Column
In order to illustrate the generality of the above research conclusions, a comparative analysis of MS-AS damage of bridge columns with different shear-span ratios is carried out.The basic parameters of the bridge column are ρ = 2.0%, ρsv = 1.5%, μ = 0.1.Time history and modal analysis of the column with shear-span ratios of 6, 10, 14, and 18 are carried out.Figure 7 shows the damage comparison chart of the bridge column under the action of MS and MS-AS of different seismic intensity.The basic period and additional damage are shown in Table 5.     Figure 7 indicates that the MS-AS earthquake demand of columns with different shearspan ratios is greater than the MS.As the shear-span ratio increases, the seismic demand gradually increases.Table 5 shows that the column damage index is less than 0.2 and the growth rate of column damage does not exceed 2.5% when PGA ≤ 0.1 g.Therefore, the effect of aftershock can be ignored.Besides, the columns with different shear-span ratios does not exceed medium damage and the growth rate of additional damage caused by aftershock does not exceed 10% when PGA = 0.2 g.When PGA = 0.3 g and 0.4 g, the damage growth rates of the column with a shear-span ratio of 6 caused by aftershock are 11.71 and 14.79%, respectively.However, the columns are slightly and severely damaged under the action of the mainshock, and the actual damage caused is limited.Furthermore, the growth rate of additional damage to columns caused by aftershock gradually decreases with the increase in shear-span ratio with different magnitudes.This situation is closely related to the additional damage of the column caused by aftershocks less than 1 and the large difference in the period of the column.For the area with a seismic level of 6 degree, the effect of aftershock can be ignored.When the seismic intensity is 7 degree, the additional damage caused by aftershock is about 4%.When the seismic intensity is 8 degree, the additional damage caused by aftershock is about 12%.When the seismic intensity is 9 degree, the additional damage caused by aftershock is about 15%.The damage growth rate of different seismic intensity caused by aftershock is close to the above research.From the above study, it can be seen that the additional damage caused by aftershocks did not exceed 15%, which is smaller than 25% proposed by Wen [8,28].The main reason is that the effect of incidence angle is not considered in this paper, and the damage is additional damage at the specified seismic wave velocity.In addition, the conclusion of this paper is greater than the 12.3% proposed by Zhang [12] because the object of study in reference [12] is wooden structures with a greater capacity to absorb seismic energy under earthquake.

Effect of Aftershock on Column with Different Damage Degree
For bridge structures in practical engineering, the occurrence of MS is often accompanied by AS.Although the intensity of AS is mostly less than that of MS, the structure is often already deteriorated when AS occur, which may lead to further increase of structural damage.Given that most of the existing research is focused on MS and less on AS, Equation ( 2) is used to quantify the additional damage caused by AS in order to make use of the existing research results for MS in structural design.The steps of the analysis in this section are as follows: (1) the IDA of MS and MS-AS is carried out for the pier; (2) the damage states of piers with different reinforcement rates caused by MS are counted and classified according to the damage states selected in Table 1; (3) the damage growth rate of AS in different damage states is counted, and the maximum value is selected as the additional damage caused by AS; (4) For piers with different axial compression ratio, the process of ( 1)-( 3) is repeated to obtain the additional damage under different damage states caused by AS with different axial compression ratio.The additional damage in different damage states is shown in Figure 8a,b.
As can be seen from Figure 8a, when the column is in different damage states, the additional damage caused by aftershocks gradually decreases with the increase of axial compression ratio.Besides, the additional damage to the column caused by aftershock gradually increases with the increase of damage when µ ≥ 0.1.Furthermore, the maximum growth rate of the additional damage to the column caused by aftershocks is 10.4% when the column is slightly damaged, and it gradually decreases as the axial compression ratio increases.However, the damage caused by the mainshock is relatively small, and therefore the actual damage caused by the aftershocks is limited.When the column is in a medium damaged state, the growth rate of additional damage caused by aftershocks is up to 12.7%, and furthermore the axial compression ratio decreases significantly between 0.05 and 0.15.Since the damage index of the medium damage state is between 0.25 and 0.4, the aftershocks have a certain impact on the actual damage of the column.When the column is in a severely damaged state, the growth rate of the additional damage caused by aftershocks decreases slightly with the increase in axial compression ratio, which can be conservatively estimated to be 11% of the mainshock damage.However, the damage caused by aftershocks is relatively large since the damage index of the severe damage state is greater than 0.4.As can be seen from Figure 8b, the additional damage caused by aftershocks gradually decreases with the increase of the shear-span ratio.Besides, when the shear span ratio is between 6 and 10, the additional damage is significantly reduced.However, due to the relatively large mainshock damage index of the columns with a shear-span ratio of 18, the additional damage caused by the aftershocks is relatively small.shocks did not exceed 15%, which is smaller than 25% proposed by Wen [8,28].The main reason is that the effect of incidence angle is not considered in this paper, and the damage is additional damage at the specified seismic wave velocity.In addition, the conclusion of this paper is greater than the 12.3% proposed by Zhang [12] because the object of study in reference 12 is wooden structures with a greater capacity to absorb seismic energy under earthquake.

Effect of Aftershock on Column with Different Damage Degree
For bridge structures in practical engineering, the occurrence of MS is often accompanied by AS.Although the intensity of AS is mostly less than that of MS, the structure is often already deteriorated when AS occur, which may lead to further increase of structural damage.Given that most of the existing research is focused on MS and less on AS, equation (2) is used to quantify the additional damage caused by AS in order to make use of the existing research results for MS in structural design.The steps of the analysis in this section are as follows: (1) the IDA of MS and MS-AS is carried out for the pier; (2) the damage states of piers with different reinforcement rates caused by MS are counted and classified according to the damage states selected in Table 1; (3) the damage growth rate of AS in different damage states is counted, and the maximum value is selected as the additional damage caused by AS; (4) For piers with different axial compression ratio, the process of ( 1)-( 3) is repeated to obtain the additional damage under different damage states caused by AS with different axial compression ratio.The additional damage in different damage states is shown in Figure 8a,b.As can be seen from Figure 8a, when the column is in different damage states, the additional damage caused by aftershocks gradually decreases with the increase of axial compression ratio.Besides, the additional damage to the column caused by aftershock gradually increases with the increase of damage when μ ≥ 0.1.Furthermore, the maximum

Effect of Axial Compression Ratio on Fragility of RC Column
This part combines the results of the IDA time history analysis above and the logarithm of the seismic peak acceleration, and the damage index is plotted as a scatter plot.The least square method is used to linearly fit the scatter plot.Only scatter diagrams with µ = 0.05, ρ = 1.0% are listed due to the limitation of space.The scatter plot is shown in Figure 9.According to the above method, the damage index of each column is linearly fitted.The statistical results of the correlation coefficient of the linear fitting of each column are shown in Table 6.
Table 6 indicates that the correlation of the linear fitting of each column decreases as the axial compression ratio increases and increases as the reinforcement ratio increases.It can be seen that the dispersion of the seismic response of the piers will gradually increase when the piers enter the plastic state from elastic under the seismic load.
According to the linear fitting equation in Table 6, the fragility curves of the columns beyond different damage states under the action of MS and MS-AS are drawn.The difference in median PGA value corresponding to each column is examined.The median value of fragility curves, referred to as the median PGA (g), is the IM corresponding to 50% probability of exceeding a specified limit state.As the median PGA increases, the fragility of the bridge is reduced.Only the fragility curves for the axial compression ratio of 0.05 and 0.2 are listed here.The fragility curves and median values are shown in Figures 10 and 11, respectively.This part combines the results of the IDA time history analysis above and the logarithm of the seismic peak acceleration, and the damage index is plotted as a scatter plot.The least square method is used to linearly fit the scatter plot.Only scatter diagrams with μ = 0.05, ρ = 1.0% are listed due to the limitation of space.The scatter plot is shown in Figure 9.According to the above method, the damage index of each column is linearly fitted.The statistical results of the correlation coefficient of the linear fitting of each column are shown in Table 6.
-3.0 -2.5 -2.0 -1.5 -           Figure 10 indicates that the damage probability of the column gradually decreases with the increase of the reinforcement ratio, and the fragility of the columns under the action of the MS-AS is greater than that of the MS.Increasing the reinforcement ratio can significantly reduce the probability of column failure when μ ≤ 0.10.However, the performance of reducing the fragility of the column by increasing the reinforcement ratio gradually decreases with the increase of the axial compression ratio.The damage probability of the column in different damage states gradually increases with the increase of axial compression ratio; nevertheless, the effect of aftershock is lower in columns with a small axial compression ratio.
Figure 10a shows that the probability of minor damage of columns with different reinforcement ratios exceeds 70% when μ = 0.2 and PGA = 0.1 g.Therefore, the axial compression ratio should be reasonably selected in design to avoid cracks of columns under small earthquakes, which can improve the durability of the structure in special areas.Figure 10b indicates that the probability of medium damage of columns with different reinforcement ratios exceeds 70% when μ = 0.2 and PGA = 0.2 g.It is not ideal to improve the ability of the column by increasing reinforcement ratio to resist medium damage.Besides, reducing the axial compression ratio can effectively improve the probability to resist medium damage.Figure 10c indicates that the probability of serious damage of columns exceeds 70% when μ = 0.2 and PGA = 0.3 g, and the probability of serious damage is reduced to less than 10% by increasing reinforcement ratio.Therefore, it is extremely uneconomical to reduce the fragility of the column by increasing the reinforcement ratio.Besides, the severely damaged column often enters the plastic stage, and the existence of aftershock tends to aggravate the collapse of the structure.Therefore, the effect of aftershock should be considered in the design.Figure 10d shows that increasing the reinforcement ratio to reduce the collapse probability is less than 20% when μ = 0.2.It is not as good as severe or medium damage by increasing the reinforcement ratio to reduce the probability of the column's collapse.Nevertheless, reducing the axial compression ratio can greatly reduce the probability of the collapse of column.
Figure 11 indicates that the median PGA50% of the column gradually increases with the increase of the damage state, and the PGA50% gradually decreases with the increase of the axial compression ratio and increases with the increase of the reinforcement ratio.When the axial compression ratio increases from 0.05 to 0.1, the reduction of PGA50% is particularly prominent.When the axial compression ratio is 0.1 to 0.15, increasing the reinforcement ratio can increase the PGA50% to a certain extent for columns with minor and medium damage.However, this trend flattens out between 0.15 and 0.2.Besides, the Figure 10 indicates that the damage probability of the column gradually decreases with the increase of the reinforcement ratio, and the fragility of the columns under the action of the MS-AS is greater than that of the MS.Increasing the reinforcement ratio can significantly reduce the probability of column failure when µ ≤ 0.10.However, the performance of reducing the fragility of the column by increasing the reinforcement ratio gradually decreases with the increase of the axial compression ratio.The damage probability of the column in different damage states gradually increases with the increase of axial compression ratio; nevertheless, the effect of aftershock is lower in columns with a small axial compression ratio.
Figure 10a shows that the probability of minor damage of columns with different reinforcement ratios exceeds 70% when µ = 0.2 and PGA = 0.1 g.Therefore, the axial compression ratio should be reasonably selected in design to avoid cracks of columns under small earthquakes, which can improve the durability of the structure in special areas.Figure 10b indicates that the probability of medium damage of columns with different reinforcement ratios exceeds 70% when µ = 0.2 and PGA = 0.2 g.It is not ideal to improve the ability of the column by increasing reinforcement ratio to resist medium damage.Besides, reducing the axial compression ratio can effectively improve the probability to resist medium damage.Figure 10c indicates that the probability of serious damage of columns exceeds 70% when µ = 0.2 and PGA = 0.3 g, and the probability of serious damage is reduced to less than 10% by increasing reinforcement ratio.Therefore, it is extremely uneconomical to reduce the fragility of the column by increasing the reinforcement ratio.Besides, the severely damaged column often enters the plastic stage, and the existence of aftershock tends to aggravate the collapse of the structure.Therefore, the effect of aftershock should be considered in the design.Figure 10d shows that increasing the reinforcement ratio to reduce the collapse probability is less than 20% when µ = 0.2.It is not as good as severe or medium damage by increasing the reinforcement ratio to reduce the probability of the column's collapse.Nevertheless, reducing the axial compression ratio can greatly reduce the probability of the collapse of column.
Figure 11 indicates that the median PGA 50% of the column gradually increases with the increase of the damage state, and the PGA 50% gradually decreases with the increase of the axial compression ratio and increases with the increase of the reinforcement ratio.When the axial compression ratio increases from 0.05 to 0.1, the reduction of PGA 50% is particularly prominent.When the axial compression ratio is 0.1 to 0.15, increasing the reinforcement ratio can increase the PGA 50% to a certain extent for columns with minor and medium damage.However, this trend flattens out between 0.15 and 0.2.Besides, the PGA 50% of the MS is larger than the MS-AS, and the PGA 50% caused by the aftershock gradually decreases with the increase of the axial compression ratio.The smaller PGA 50% due to subsequent earthquakes for a certain reinforcement ratio show the importance of aftershock in the progressive damage to the column and further indicates the necessity to capture aftershock effects on fragility assessment of bridges, which has not been considered by current codes to date.The above shows that the optimal design axial compression ratio under the action of the MS-AS should be between 0.05 and 0.1.In this interval, increasing the reinforcement ratio can better improve the PGA 50% of the column, and the effect of aftershock should be considered based on the damage of the mainshock.

Effect of Shear-Span Ratio on Fragility of RC Column
To illustrate the generality of the above research, the fragility analysis of columns with different shear-span ratios is carried out.The scatter diagram is drawn by the above method and fitted linearly.The linear fitting correlation coefficients are shown in Table 7.According to the linear fitting equation in Table 7, the fragility curves of the columns beyond different damage states under the action of MS and MS-AS are drawn, and the median value of fragility curves is obtained.Due to space limitations, only the fragility curves of medium and severe conditions are listed.The fragility curves and the median fragility PGA 50% of columns with different shear-span ratios are shown in Figures 12a-d  PGA50% of the MS is larger than the MS-AS, and the PGA50% caused by the aftershock gradually decreases with the increase of the axial compression ratio.The smaller PGA50% due to subsequent earthquakes for a certain reinforcement ratio show the importance of aftershock in the progressive damage to the column and further indicates the necessity to capture aftershock effects on fragility assessment of bridges, which has not been considered by current codes to date.The above shows that the optimal design axial compression ratio under the action of the MS-AS should be between 0.05 and 0.1.In this interval, increasing the reinforcement ratio can better improve the PGA50% of the column, and the effect of aftershock should be considered based on the damage of the mainshock.

Effect of Shear-Span Ratio on Fragility of RC Column
To illustrate the generality of the above research, the fragility analysis of columns with different shear-span ratios is carried out.The scatter diagram is drawn by the above method and fitted linearly.The linear fitting correlation coefficients are shown in Table 7.According to the linear fitting equation in Table 7, the fragility curves of the columns beyond different damage states under the action of MS and MS-AS are drawn, and the median value of fragility curves is obtained.Due to space limitations, only the fragility curves of medium and severe conditions are listed.The fragility curves and the median fragility PGA50% of columns with different shear-span ratios are shown in Figure 12a Table 7 shows that the logarithmic correlation between the damage index and the acceleration gradually decreases with the increase of the shear-span ratio.Figure 12 shows that the vulnerability of columns under the action of MS-AS is greater than MS.Besides, the fragility of the column for a certain damage state decreases with the decrease of the shear-span ratio, but the effect of aftershock on the fragility gradually increases.As can be seen from Figure 12a, the smaller the shear-span ratio, the greater the effect of aftershock.Figure 12b indicates that the fragility curves with shear-span ratios of 10, 14, and 18 are little different in the severely damaged state.This situation is due to: (1) The period of the columns changes greatly when the shear-span ratio is 10, 14, 18, and the corresponding response spectrum of the period changes greatly.(2) The plastic deformation of the column is irregular and damage index is less than or equal to 1.
As can be seen from Figure 13, the median fragility PGA50% of the column gradually decreases with the increase of the shear-span ratio and gradually increases with the increase of the damage state.Besides, the PGA50% under the MS-AS is smaller than the MS, and the effect of aftershock on the PGA50% decreases as the shear-span ratio increases.Reducing the shear-to-span ratio can greatly reduce the damage probability of bridge piers and increase the ability of piers to resist aftershocks.

Conclusions
Based on the OpenSees analysis platform, the paper selects the mainshocks and aftershocks with shear wave speeds between 350 and 500 m/s for IDA analysis of the bridge columns.The effect of the columns with different axial compression ratios, reinforcement ratios and shear-span ratios is discussed.The conclusions are: 1.The seismic demand of the column under the action of the MS-AS sequences is greater than MS.This demand increases with the increase of axial compression ratio and decreases with the increase of reinforcement ratio.When the column is in a slight or no damage state, the effect of aftershock can be basically ignored.2. The additional damage caused by AS can be comprehensively considered in the design according to different magnitudes.When the seismic level is 7 degree, the maximum increase rate of additional damage is about 7% of MS.When the seismic level Table 7 shows that the logarithmic correlation between the damage index and the acceleration gradually decreases with the increase of the shear-span ratio.Figure 12 shows that the vulnerability of columns under the action of MS-AS is greater than MS.Besides, the fragility of the column for a certain damage state decreases with the decrease of the shear-span ratio, but the effect of aftershock on the fragility gradually increases.As can be seen from Figure 12a, the smaller the shear-span ratio, the greater the effect of aftershock.Figure 12b indicates that the fragility curves with shear-span ratios of 10, 14, and 18 are little different in the severely damaged state.This situation is due to: (1) The period of the columns changes greatly when the shear-span ratio is 10, 14, 18, and the corresponding response spectrum of the period changes greatly.(2) The plastic deformation of the column is irregular and damage index is less than or equal to 1.
As can be seen from Figure 13, the median fragility PGA 50% of the column gradually decreases with the increase of the shear-span ratio and gradually increases with the increase of the damage state.Besides, the PGA 50% under the MS-AS is smaller than the MS, and the effect of aftershock on the PGA 50% decreases as the shear-span ratio increases.Reducing the shear-to-span ratio can greatly reduce the damage probability of bridge piers and increase the ability of piers to resist aftershocks.

Conclusions
Based on the OpenSees analysis platform, the paper selects the mainshocks and aftershocks with shear wave speeds between 350 and 500 m/s for IDA analysis of the bridge columns.The effect of the columns with different axial compression ratios, reinforcement ratios and shear-span ratios is discussed.The conclusions are: 1.
The seismic demand of the column under the action of the MS-AS sequences is greater than MS.This demand increases with the increase of axial compression ratio and decreases with the increase of reinforcement ratio.When the column is in a slight or no damage state, the effect of aftershock can be basically ignored.

2.
The additional damage caused by AS can be comprehensively considered in the design according to different magnitudes.When the seismic level is 7 degree, the maximum increase rate of additional damage is about 7% of MS.When the seismic level is 8 degree, the maximum increase rate of additional damage is about 13% of MS.

Figure 1 .
Figure 1.Research content and significance.

Figure 1 .
Figure 1.Research content and significance.

Figure 3 .
Figure 3.The basic parameters of the column.

Figure 3 .
Figure 3.The basic parameters of the column.

Figure 7 .
Figure 7. Damage index of the column.

Figure 7 .
Figure 7. Damage index of the column.

Figure 9 .
Figure 9. Scatter diagram of the column under MS-AS.

Figure 9 .
Figure 9. Scatter diagram of the column under MS-AS.

Table 1 .
Range of damage index for different damage states.

Table 2 .
Basic parameters of reinforcement.

Table 2 .
Basic parameters of reinforcement.

Table 3 .
Basic period (T) of the columns.

Table 3 .
Basic period (T) of the columns. .
Figure7indicates that the MS-AS earthquake demand of columns with different

Table 6 .
Linear fitting equation of the columns.

Table 7 .
Linear fitting equation of the columns.

Table 7 .
Linear fitting equation of the columns.