A Multimode Adaptive Pushover Procedure for Seismic Assessment of Integral Bridges EhsanMohtashami and

This paper presents a new adaptive pushover procedure to account for the effect of higher modes in order to accurately estimate the seismic response of bridges.The effect of higher modes is considered by introducing a minimum value for the total effective modal mass. The proposed method employs enough number of modes to ensure that the defined total effective modal mass participates in all increments of the pushover loading. An adaptive demand curve is also developed for assessment of the seismic demand. The efficiency and robustness of the proposed method are demonstrated by conducting a parametric study. The analysis includes 18 four-span integral bridges with various heights of piers. The inelastic response history analysis is employed as reference solution in this study. Numerical results indicate excellent accuracy of the proposed method in assessment of the seismic response. For most bridges investigated in this study, the difference between the estimated response of the proposed method and the inelastic response history analysis is less than 25% for displacements and 10% for internal forces. This indicates a very good accuracy compared to available pushover procedures in the literature. The proposed method is therefore recommended to be applied to the seismic performance evaluation of integral bridges for engineering applications.


Introduction
Over the past two decades, there has been an increasing attention to Nonlinear Static Procedures (NSPs) for the seismic assessment and evaluation of bridges.In these methods, a pushover analysis is carried out to estimate the inelastic capacity of the structure as well as its response to different levels of the seismic demand.The NSPs have the simplicity of the linear static methods and the accuracy of the nonlinear dynamic analyses.Therefore, they have become standard tools for the analysis of bridges in many seismic codes and guidelines such as EC8-2 [1] and Caltrans [2].However, the principles of the NSPs are originally developed for buildings.It has been indicated by Isaković and Fischinger [3] that there are major differences between the structural system and the seismic response of buildings and bridges, particularly when the response of a bridge is analyzed in transverse direction.Therefore, the application of NSPs should be modified and adjusted for bridges.
The NSPs for the seismic assessment of bridges are generally based on two classic pushover procedures: the capacity spectrum method (CSM) [4] and the inelastic demand spectrum method also known as the N2 method [5].Both methods consider a single mode and keep the shape of the pushover load distribution constant during the analysis.To overcome these deficiencies, several approaches have been proposed such as Modal Adaptive Nonlinear Static Procedure (MANSP) [6,7] and Incremental Response Spectrum Analysis (IRSA) [8].In these methods, the effect of higher modes is included in the analysis and the pushover load is considered a function of the updated modal characteristics of the structure.
In order to improve the efficiency of adaptive pushover methods for bridges, various studies have been conducted.Isaković and Fischinger [3] employed MANSP and IRSA methods and investigated the effect of higher modes in the seismic assessment of continuous-deck viaducts.Both MANSP and IRSA performed well for regular bridges.However, the accuracy of these methods was reported to decrease as the irregularity of the bridge increases.Casarotti and Pinho [9] developed a displacement-based procedure called Adaptive Capacity Spectrum Method (ACSM).This approach was shown to provide satisfactory predictions in terms of displacements and moments.Pinho et al. [10] employed four pushover procedures and conducted a parametric study for the seismic assessment of 14 continuous-span bridges.They concluded that the accuracy of pushover methods is generally acceptable.However, the dispersion of results was reported to be high for irregular bridges.Therefore, finding a more reliable method requires further investigation.
The main contribution of this paper is to present a refined multimode adaptive capacity spectrum approach in order to enhance the accuracy and reliability of NSPs in the seismic assessment of bridges.To do this, two improvements have been made in estimating the seismic capacity and demand.First, the capability to identify and include important higher modes in the pushover analysis is provided.This is done by introducing a minimum value for the total effective modal mass.The adequate number of modes is then included to ensure that the defined total effective modal mass participates in all increments of the pushover loading.In other words, the proposed method is able to incorporate variable number of modes in the analysis.Second, an adaptive demand curve is developed based on consecutive updating of the seismic demand.The derived adaptive demand curve is then intersected by the adaptive capacity curve to determine the performance point of the structure.To the best of authors' knowledge, the idea of using an adaptive scheme of both capacity and demand curves in the process of seismic assessment is novel and has not been presented yet.The formulation of the proposed method for obtaining the capacity curve of the structure is similar to ACSM [9].However, the proposed method is basically a force-based pushover procedure while the ACSM is a displacement-based approach.Unlike ACSM, the number of the participating modes is not assumed fixed in the proposed method.Moreover, the capacity and demand curves are updated step by step and their intersection is checked in each increment for detection of the performance point.This provides a more straightforward procedure compared to ACSM in which the capacity curve is first obtained and the seismic demand is then applied in an iterative procedure to determine the performance point.
The paper is organized as follows.Section 2 presents the proposed assessment methodology as well as the new improvements implemented in the procedure.To evaluate the accuracy and efficiency of the proposed procedure, a parametric study is conducted in Section 3. Section 4 presents the numerical results of the parametric study and discusses the suitability and adequacy of the proposed method.Section 5 highlights concluding remarks.

The Proposed Seismic Assessment Method
The proposed method is performed in the following main steps: (1) assessment of the capacity curve of the bridge by conducting a multimode adaptive pushover analysis, (2) estimation of the adaptive demand curve based on the updated overdamped elastic response spectrum of the given ground motion, (3) determination of the performance point, and (4) determination of displacements and internal forces of the bridge.Details of these steps are explained as follows.

Assessment of the Adaptive Capacity
Curve.The first step of the proposed procedure is to perform a pushover analysis on the bridge.A fiber based distributed plasticity finite element program is implemented in OpenSees [11] for the seismic assessment of bridges.The following steps are performed at each increment of the pushover loading.
(i) Conduct an eigenvalue analysis to find mode shapes.The proposed method provides the option to define an arbitrary value for the total effective modal mass,  eff , to be participated in the pushover analysis.For instance, if the common value of 90% of the total mass of the structure is selected for  eff , our computer code will include an adequate number of modes in the assembly of the pushover load such that at least 90% of the total mass of the structure participates in the analysis.Since the method is an incremental procedure, this check for  eff is repeated at each load increment.Therefore, the number of the participating modes and consequently the shape of the pushover load vector repeatedly vary during the analysis.As a result, the changes of modal characteristics due to gradual softening of all structural members are well captured and included in the pushover curve.
(ii) Assemble the incremental pushover load vector, {Δ  }, as where Δ base is increment of the base shear defined at the beginning of the analysis,   is the normalized modal force at th degree of freedom (DOF),   is the resultant modal force applied at DOF  based on the Square-Root-of-Sumof-Squares (SRSS) combination rule, m is the number of the participating modes, and   is the modal force of DOF  at mode  obtained as where   is the mass of the th DOF, Φ  is the mode shape of DOF  at mode , Γ  is the modal participation factor of mode , and  , is the spectral pseudo-acceleration of mode  determined from the updated overdamped elastic response spectrum of the given input motion.The procedure for determining this spectrum will be explained in Section 2.2.
(iii) Perform an incremental-iterative analysis by employing the Newton-Raphson numerical method [12] as where {Δ} is the incremental displacement vector of the bridge and [  ] is the tangent stiffness of the structure including geometrical and material nonlinearity.(iv) Derive the adaptive capacity curve of the equivalent single-degree-of-freedom system by calculating the equivalent displacement  ,cap and the equivalent acceleration  ,cap .These values are calculated in kth step of the analysis as [9] ( ,cap where  is the gravity,  base, is the total base shear of the structure at step , and Δ , is the displacement of the mass   at step .Note that the equivalent displacement and acceleration of the structure do not refer to a specific physical location of the bridge and are updated based on the deformed shape of the structure.

Estimation of the Adaptive Demand
Curve.At each step of the proposed procedure, the response spectrum of the ground motion is updated to match the current state of the structure.To do this, the equivalent damping of the structure is determined at each load step by the damping-ductility relationship of the Takeda degrading stiffness hysteretic response [13] as where  0 is the viscous damping of the bridge which is assumed as 5% of the critical damping, and  represents the postyield to preyield stiffness ratio and is assumed to be 0.05 [14].The equivalent ductility of the structure,  eq , is defined as the average ductility of all lateral resistant members, weighted by the shear force of each member.The reason is that the overall ductility of the structure is related to both ductility and stiffness of the members that contribute to the lateral resistance.The equivalent ductility of the bridge is then calculated as where  accounts for the effect of the soilabutment interaction which is explained in Section 3.2.

Determination of the Performance Point.
In each analysis increment, the adaptive capacity curve is compared with the adaptive demand curve for detection of the performance point (see Figure 1).In initial steps of the analysis, the structure is elastic and the demand curve proceeds on the elastic response spectrum.When the first structural member starts to yield, the equivalent damping of the structure begins to increase according to (5) and (6).The demand curve is then shifted from the elastic spectrum to the derived spectrum.As yielding spreads within the structure, the demand curve is shifted from the spectrum of the previous step to the spectrum of the new step.When an intersection point is detected, it is saved in an output file as a potential performance point.However, the analysis continues because more than one intersection between the capacity and demand curves may be found [9].The analysis stops when an ultimate limit state is met.Ultimate limit states are considered to occur either when the structure fails to resist the pushover load or when the columns reach a drift limit of 3% [15].

Determination of the Structural Response.
Once the performance point is detected, the structural response can be determined in terms of displacements and forces.The flowchart of the proposed procedure is schematically summarized in Figure 2. To evaluate the efficiency of the proposed method, a parametric study is conducted in the next section.(This is a possible structural response, because more performance points might be detected during next steps of the analysis.The correct performance point is determined by the user at the end of the analysis.)(I) Assess the adaptive capacity curve (1) Run an eigenvalue analysis for determining mode shapes and the number of important modes (2) Assemble pushover load pattern by SRSS combination rule (3) Apply load increment and conduct nonlinear static analysis (4) Save displacements and internal forces into the database Is an ultimate limit state achieved?
(1) Define input parameters (geometrical and material properties) (2) Define a minimum value for the total effective modal mass (III) Check for performance point Is a performance point detected?
(1) Determine (S d , S a ) demand using the overdamped elastic response spectrum with  eq calculated in the previous analysis step (2) Update  eq for the next step by employing equation (5) (5) Calculate parameters of the equivalent SDOF system: (S d , S a ) capacity

Description of the Parametric Study
The parametric study consists of 180-meter long four-span integral bridges supported on pile foundations.The bridges are primarily designed according to AASHTO LRFD Bridge Design Specifications [16] and Caltrans Seismic Design Criteria [2] for the purpose of this study.Structural details of the investigated bridges are illustrated in Figures 3(a)-3(e).In this study, eighteen bridges with different pier heights are considered.These bridges are denoted by Bijk where , , and  are label numbers that indicate the height of left, middle, and right pier, respectively.The height of each pier can be 6, 12, or 18 meters which will be labeled 1, 2, and 3, respectively.For instance, the bridge in Figure 3(a) is denoted by B123.A fiber-based distributed plasticity finite element program is implemented in OpenSees for the proposed seismic assessment method.The minimum total effective modal mass is defined as 90% of the total mass of the structure.Geometrical nonlinearity is taken into account for all structural elements.The potential of inelasticity is restricted to columns.Nonlinear Winkler springs are employed to represent the reaction of soil on piles and abutments.The accuracy of the proposed method is evaluated by inelastic response history analysis (IRHA).For each of the 18 bridges, the analysis is performed under 10 strong ground motions.This implies that a total number of 180 analyses have been carried out to evaluate the proposed method.The constitutive law of the concrete is described by Mander's confined concrete model [17].The modulus of elasticity of steel reinforcing bars is 200 GPa.The stress-strain relationship of steel is defined by Menegotto-Pinto model [18] with yield stress of 450 MPa.The total mass of the bridge includes the mass of the superstructure, cap beams, and half of piers.The total weights of the investigated bridges are in the range of 50100-50800 kN, depending on the height of their piers.The mass of each element is equally distributed at its end nodes.The bridge deck is a multicell box girder as shown in Figure 3(c).The deck is modeled by 3D beam-column elements fully characterized by elastic sectional properties.A 50% reduction of elastic flexural rigidity (EI) and an 80% reduction of elastic torsional constant (GJ) are considered according to recommendations given in Caltrans [2] for multicell box girder sections.The superstructure has monolithic connections to abutments and piers.
The bridge piers consist of two circular columns as illustrated in Figure 3(c).The cross-section of columns is depicted in Figure 3(d).A rigid end zone of 1.0 m is located at the top of columns to account for the offset between the top of columns and the structural nodes of the superstructure.The columns are modeled by 3D nonlinear beam-column elements with fiber cross-sections.The connection between columns and the pile cap is considered to provide full continuity.The pile cap is supported by 8 steel friction piles, as illustrated in Figure 3(e).The piles are 15 meters long and are modeled by 3D nonlinear beam-column elements.The soil stratum under the whole bridge is assumed to be stiff cohesionless sand with a unit weight of 18 kN/m 3 and an internal friction angle of 40 degrees.All piles lie on the bedrock.

Soil-Structure
Interaction.An important issue in the analysis of integral bridges is how to model the interaction of abutments, foundations, and piles with their surrounding soil.This interaction has significant influence on the seismic response and thus should be carefully included in the structural model of the bridge.In order to realistically model the soil-structure interaction, nonlinear springs are utilized to simulate the horizontal reaction of the soil around abutments and piles.Since the soil-structure interaction is not the main scope of this study, characteristics of soil springs are kept constant for all investigated bridges.The hyperbolic tangent curve recommended by the American Petroleum Institute [19] is employed for the lateral force-displacement relationship (p-y curve) of soil springs around piles.These springs are assigned along the length of piles under abutments and pier foundations.The p-y curves at certain depths are plotted in Figure 4(a).The soil-abutment interaction in the transverse direction also includes the contribution of wingwalls.The effect of wingwalls is represented by a bilinear translational SZO011-NS 0.511 * The first 7 records are selected from [21] and the last three records from [22].
spring as shown in Figure 4(b).This spring is assigned to the top of the abutment in the transverse direction.The stiffness of this spring is determined by the shear capacity of one wingwall [20].The initial stiffness  1 and the postyield stiffness  2 in Figure 4(b) are calculated based on elastic and ultimate shear capacities of the wingwall.

Ground Motion Records.
The employed seismic excitation is defined by a set of ten natural near-fault ground motions from strong and major earthquakes around the world.Seven records are selected from the Pacific Earthquake Engineering Research Center strong-motion database [21] and three records from the Kyoshin strong-motion seismograph network [22].The range of magnitude for the selected records is from 6.4 to 7.6, and the range of epicentral distance is from 2 to 22 km.The characteristics of the input ground motions are summarized in Table 1.The selected earthquake records cover a wide variety of near-fault ground motion parameters.In this study, the PGA of the records is scaled to 0.4 g in order to represent the design earthquake level.The effect of different intensity levels on the accuracy of the proposed method is later investigated in Section 4.4.The ground motions are applied in the transverse direction of the bridge.

Statistical Analysis.
Results of the parametric study are presented in terms of various response parameters among which, deck displacements, pier base displacements, and column base moments are investigated here.In order to evaluate the accuracy of the proposed method, the response parameters are normalized to the response of IRHA for each input ground motion.This provides a good measure for evaluation of the proposed method since all quantities are normalized.
Due to random nature of earthquakes, any seismic response of a structure can be considered as a random variable.This response is a function of the input ground motion as well as the DOF at which it is computed.In order to properly characterize the response of the structure, the behavior of each response parameter has to be studied at different DOFs subjected to various ground motions.In other words, the input ground motion and the DOF should be statistically modeled.As a result, an arbitrary response parameter, say R, can be considered as a random variable.By   , we mean the response parameter  at the th DOF subjected to the jth ground motion.To evaluate the accuracy of our proposed method,   is calculated by our nonlinear static procedure and then compared to the result of the IRHA as the reference solution.This is mathematically represented as where the superscripts NSP and IRHA denote the proposed nonlinear static procedure and the inelastic response history analysis, respectively.Note that we have selected the relative difference,   as the likelihood measure of our proposed method to the reference solution.If the proposed method estimates the seismic response exactly as the IRHA,   is equal to zero.The mean of the random variable  is obtained as where   is the total number of DOFs at which the response parameter is computed and   is the number of input ground motions applied to the bridge.The standard deviation of  indicates the dispersion of the response predicted by the proposed method compared to the IRHA result.It is written as The necessary condition for accepting the accuracy of the proposed method is that the mean of the relative difference should be close to zero.This condition ensures that the proposed procedure is potentially able to properly predict the results obtained by the IRHA.The sufficient condition, on the other hand, is satisfied when the standard deviation of the relative difference is sufficiently small.This condition guarantees that the results estimated by the proposed method are always close to the IRHA results.Therefore, the proposed method is considered accurate if and only if both mean and standard deviation of the relative difference are sufficiently small.
Based on the computed standard deviation, it is useful to consider a confidence range for the relative difference as a measure of accuracy.The 95% confidence is a common value in probabilistic and reliability analysis of structures and is therefore accepted in this study.The 95% confidence range for a normally distributed random variable  is from  − 2  to  + 2  where  and   denote the expected value and the standard deviation of x, respectively.If the dispersion of the relative difference is low, the 95% confidence range is small and the estimation of the response parameter is considered accurate.The normality of the relative difference can be verified by the nonparametric Kolmogorov-Smirnov (K-S) test [23].This test is performed by the widely used statistical software SPSS.All relative differences for the response parameters of this study have successfully passed the K-S test and thus have been verified to follow the normal distribution.Numerical results of the parametric study are presented in the next section.

Discussion on the Numerical Results
The eighteen bridges considered in this study can be geometrically classified as regular and irregular based on the definition of the regularity index proposed by Calvi et al. [24].According to the comprehensive study of Akbari and Maalek [25] on the regularity of bridges, the bridges B111, B121, B131, B222, B232, B323, and B333 of this study are identified as regular, and the other 11 bridges are irregular.The effect of regularity on the response of the proposed method will be investigated in Section 4.2.The accuracy of the proposed method is evaluated by the inelastic response history analysis having a 2% Rayleigh damping assigned to the two modes with the highest effective modal mass.

Results of the Statistical Analysis.
Results of the statistical analysis of Section 3.4 are presented here.The mean and 95% confidence range of the relative difference  are presented in Figures 5, 6, and 7 for deck displacements, pier base displacements, and column base moments, respectively.The standard deviations of the relative difference (  ) are shown in Figure 8. Results for each response parameter are discussed in detail as follows.
The selected DOFs for the deck displacement response are located at every 15 meters along the length of the bridge.This implies a total of 11 DOFs for this response parameter.Considering the 10 input ground motions of Table 1, a total number of 110 realizations are provided for the relative difference .The K-S test shows that this sample follows the normal distribution very well.The mean, , ranges from −8% to +5% which indicates that the proposed method has the    necessary condition to accurately estimate the displacements of the superstructure.With 95% confidence, the relative difference for 12 bridges shown in Figure 5 is less than 25%.This means that for two-thirds of the investigated bridges, the proposed method confidently estimates the deck displacements with less than 25% difference compared to the IRHA results.The standard deviation of the response obtained by the proposed method is in the range from 8% to 17%.As reported by Pinho et al. [10], the standard deviation for deck displacements is up to 40% in different pushover procedures.Therefore, the proposed method satisfies the sufficient condition of accuracy for estimating the deck displacements.It can be concluded then that the proposed method estimates the deck displacement response with satisfactory accuracy.The selected DOFs for the response of each pier base are the displacement in the transverse direction and the rotation about the longitudinal axis of the bridge.Therefore, 6 DOFs are selected for each bridge.Considering the 10 input ground motions, a total number of 60 realizations are included in the sample.The result of the K-S test indicates that the random variable  follows the normal distribution.The mean, , ranges from −4% to +2% which indicates the excellent adequacy of the necessary condition.As can be seen in Figure 8, the standard deviation ranges from 2% to 9%.According to the results of Pinho et al. [10], the standard deviation is low enough to satisfy the sufficient condition of accuracy.It is then concluded that the proposed method excellently estimates the displacement response of the pier base.

Advances in Civil Engineering
For the column base moments, 6 DOFs corresponding to six columns of the bridge are considered.Employing the 10 ground motions of Table 1, a total number of 60 realizations are obtained in the sample.The mean, , ranges from −2% to +3% which indicates that the proposed method provides the necessary condition excellently.As can be seen Advances in Civil Engineering  b  is the period of vibration (in seconds).c  (%) is the effective modal mass divided by the total mass of the structure (given in percentage).
in Figure 7, the proposed method has confidently predicted column base moments in all the eighteen bridges with less than 15% difference compared to the IRHA results.The maximum standard deviation of the relative difference is 7% (see Figure 8).This indicates the remarkable accuracy of the proposed method.Pinho et al. [10] reported the standard deviation of the column shear force (and consequently the base moment) up to 20%.Therefore, the standard deviation obtained by the proposed method well satisfies the sufficient condition of accuracy for estimating the column moment response.Therefore, it can be concluded that the proposed method estimates the column base moment response with desirable accuracy.

Effect of Bridge
Regularity on the Results.The overall effect of bridge regularity on the response of the proposed method is investigated by evaluating the standard deviation of the random variable  as shown in Table 2.The average standard deviation of the seven regular and the eleven irregular bridges is calculated and compared, for the selected response parameters.It is observed that the accuracy of the proposed method in estimating all three response parameters is slightly more reliable for regular bridges.However, this difference is small and the bridge regularity has negligible effect on the accuracy of the proposed procedure.

Case Study Results
. In order to highlight the features of the proposed method, a case study is explained here.
To maintain the generality of discussion, bridge B123 is considered because its piers have different heights and the symmetry does not affect its seismic response.The Northridge record (with the PGA scaled to 0.4 g) is also selected among the records of Table 1. Figure 9 shows the variation of the pushover load pattern acting on the bridge deck.At the beginning of loading, the columns are elastic and their lateral stiffness is considerable.The bridge deck can therefore be considered as a multispan continuous beam in the transverse direction.The intensity of the pushover load in the middle of the bridge is higher than that in the regions near the abutments (see dash lines in Figure 9).As the loading increases, the columns start to yield and their stiffness degrades.Therefore, as seen in Figure 9(a), the intensity of the pushover load shifts from the middle of the bridge to its ends where more stiffness is provided by abutments.As the loading continues and more inelastic deformation is introduced in columns, more load shifts from the middle of the bridge towards abutments as shown in Figure 9(b).
Advances Gradually, the stiffness of abutments degrades and some load is redistributed from the regions near abutments back to the middle of the bridge, as can be seen in Figure 9(c).As the structure approaches the ultimate limit state, columns experience extensive inelastic deformations and their stiffness becomes very small.Therefore, the bridge deck behaves like a single-span beam supported by two abutments.At this point, only the first mode participates in the analysis and, consequently, the pushover load distribution resembles a parabolic shape as shown in Figure 9(d).
Table 3 shows modal characteristics of B123.It is observed that three modes have participated in the analysis.As columns start to yield consecutively, the first mode gradually tends to dominate the vibration of the bridge.The effect of the second mode is negligible, whereas the third mode with the modal mass ratio of about 20% significantly participates in the vibration of the bridge.The period of vibration for the three modes increases as the columns yield.This increase is less significant for higher modes compared to the first mode.The proposed method has been capable of presenting accurate quantitative measure of the bridge modal characteristics.

Effect of Seismic Intensity on the Results.
To investigate the effect of seismic intensity on the accuracy of the proposed method, an appropriate intensity measure should be selected.In this study, PGA is chosen as intensity measure since it has been identified as suitable for the seismic assessment of bridges [26].Here, the PGA of Tabas record is selected and scaled to typical values in the range from 0.2 g to 1.2 g.The records are then applied to one regular bridge (B222) and one irregular bridge (B123).Figure 10 shows the deck displacement pattern for different seismic intensities.It can be seen that the results of the proposed method are almost identical to IRHA results.Only for PGA of 1.2 g, the proposed method underestimates the displacements of the two middle spans, whereas the displacements of the two end spans are still estimated properly.Figure 11 shows column base moments for different intensity levels.The proposed method has been able to estimate the moments very close to the IRHA results in all intensity levels.Therefore, it can be concluded that the intensity of the ground motion does not have a significant effect on the accuracy and robustness of the proposed method.Furthermore, the deformed shape of the bridge and the column base moments are very well captured for both regular and irregular bridges.Therefore, the bridge regularity does not affect the performance of the proposed procedure in different intensity levels.

Conclusions
An efficient nonlinear static procedure for the seismic assessment of integral bridges has been presented in this paper.The proposed method is basically a multimode adaptive capacity spectrum method with improvements in the assessment of capacity and demand curves.The proposed method employs incremental-iterative analysis and updates both capacity and demand curves at each load increment.The performance point is defined as the intersection of the two adaptive curves.The displacements and internal forces are then determined accordingly.
A parametric study has been conducted to evaluate the accuracy and suitability of the proposed method.The inelastic response history analysis has also been selected as the reference solution.Numerical results of the parametric study indicate that the proposed method estimates the displacement and force responses of the simulated bridges close to the results obtained by the reference solution.For most of the bridges investigated in this study, this difference is less than 25% for prediction of displacements.The results are more accurate for column moments with average difference of less than 10%.Besides, the dispersion of results is low such that the maximum standard deviation is 17% for displacements and 7% for column moments.These values are considered quite acceptable compared to the values reported in the literature.
Numerical results also show that the regularity of the bridge has negligible effect on the accuracy of the proposed procedure.The proposed method has been able to predict both displacement and force responses of regular and irregular bridges with similar accuracy.Moreover, the intensity of the ground motion does not affect the accuracy of the proposed method in estimating the force response of the bridge.It is observed, however, that the difference between the displacement response of the proposed method and that of the reference solution slightly increases as the intensity of the ground motion increases.Nevertheless, this increase has not affected the overall accuracy and the proposed method performs satisfactorily for most intensity levels.
It is finally emphasized that the proposed procedure can also be utilized for other types of bridges with small adjustments.Besides, this method provides a robust tool for development of fragility curves and evaluation of the seismic performance of bridges.

Figure 1 :
Figure 1: The proposed procedure for detection of the performance point.
are the displacement ductility of pier  and abutment  at step . pier,  and  abt,  are the shear force of pier  and abutment  at kth step of the analysis.The term  abt, Estimate the adaptive demand curve (IV) Determine the structural response Go back to the database and determine the displacements and internal forces.

Figure 2 :
Figure 2: Flowchart of the proposed seismic assessment procedure.

3. 1 .
Structural Modeling.The three-dimensional (3D) finite element model of a typical simulated bridge is shown in links (to represent the pile cap) Rigid end zone (to represent the cap beam) m = 180 m H 1 (6 m) Pile length = 15 m Abutment height = 4 m H 2 (12 m) H 3 (18 m) (a) Elevation view (b) 3D analytical model in OpenSees HP360 × 174

Figure 3 :
Figure 3: Typical details of the investigated bridges.

Figure 3 (
Figure 3(b).The material of both superstructure and substructure is cast-in-place reinforced concrete.The concrete has Young modulus of 28 GPa and compressive strength of 35 MPa.The constitutive law of the concrete is described by Mander's confined concrete model[17].The modulus of elasticity of steel reinforcing bars is 200 GPa.The stress-strain relationship of steel is defined by Menegotto-Pinto model[18] with yield stress of 450 MPa.The total mass of the bridge includes the mass of the superstructure, cap beams, and half

Figure 4 :
Figure 4: Soil-structure interaction at abutments and foundation of piers.

Figure 6 :
Figure 6: Mean and 95% confidence range for pier base displacements and rotations.

Table 1 :
List of the selected input ground motions.

Table 2 :
Accuracy of the proposed method for regular and irregular bridges.