Response Modification Factors for Multi-Span Reinforced Concrete Bridges in Pakistan

: In Pakistan, updated codes covering seismic provisions for reinforced concrete bridges do not exist. The majority of the bridge design uses different versions of AASTHO-LRFD provisions. Response modification factors recommended for usage in these codes are primarily for bridges derived from conditions and bridges in the United States of America. This research focuses on the seismic assessment of three real multi-spans simply supported reinforced concrete bridges in Pakistan having multiple bents. This typology of bridges is very common in Pakistan. Non-linear static pushover analysis is performed to derive seismic capacity curves for these bridges, which were used to compute response modification factors. The study results show that response modification factors vary between 4.50 and 5.0 for the bridges in the longitudinal and transverse directions. The results of this work may serve as input in developing the seismic design code of bridges in Pakistan.


Introduction
There are approximately 6000 bridges on the national highways of Pakistan, and 67% of these were built before 1980.The seismic provisions followed for designing bridges in Pakistan are known as the West Pakistan Highway Code (WPHC) [1].The WPHC was adopted from AASTHO (1961).However, since then, the WPHC was never updated.Instead, different versions of AASTHO-LRFD codes were consulted to design bridges.In the WPHC, the design base lateral force equals 2% to 6% of the seismic weight (i.e., 0.02 to 0.06 W).However, it lacks a reliable link between the site seismic hazard and design base lateral force.Kawashima [2] proposed seismic isolation technology for highway bridges.After the Kobe 1995 earthquake event in Japan, the author proposed guidelines and specifications for this technique.Sakellariadis et al. [3] revisited the collapse of the 18-span bridge during the tragic seismic event of Kobe in Japan in 1995.Putra et al. [4] conducted a case study of risk analysis of seismic bridge damage after the Lombok and Palu earthquakes in Indonesia.The authors performed an on-site visual inspection of 38 bridges and reported several damages to the bridges.Memari et al. [5] evaluated the seismic response of five existing bridge piers in the eastern region of the US seismic zone.The AASTHO-LRFD suggests seismic hazard having a return period of 1000 years for calculating seismic actions for bridges.However, such seismic hazard maps are not available to date, which will lead to great destruction because significant energy is released during seismic events [6].In a case study of the Tawarayama Bridge under the 2016 Kumamoto earthquake, Aye et al. [7] reported damage mechanisms based on seismic response analysis induced by an earthquake in a plate girder bridge.The Halabjah earthquake, which struck Iraq in November 2017 and damaged numerous bridge piers and commercial buildings, prompted Al-taie et al. [8] to conduct a case study about the seismic risks associated with Iraqi soil.The seismic hazard map of Pakistan was updated with the advent of the Building Code of Pakistan (BCP-SP-, 2007) [9].However, this suggests values for peak ground acceleration (PGA) having a 10% probability of exceedance in 50 years within 475 years of the return period.Recently, S. Khan and Waseem [10] developed relevant latest probabilistic seismic hazard analysis (PSHA) maps for various return periods to evaluate Pakistan seismic zones.R. Khan et al. [11] conducted a study of the existing nine precast bridges to assess the vulnerability of bridges in Karachi to seismic ground motion.
Many improvements were made at the time by scientists and researchers to deal with seismic hazards, along with seismic site behavior.Numerous studies were conducted to depict structural behavior to dissipate seismic energy by causing adequate ductility to structures, such as in 1998 when AASHTO presented a report titled "Guide Specification for Seismic Isolation Design" [12].Research results led to the establishment of an appropriate response modification factor for a substructure for seismic activity.The response modification factor given in the US codes FEMA1997, and UBC 1997 [13,14] are intended to account for both overstrength and ductility factors (ATC-1978) [15] available in the structural system to respond to design base earthquakes, which can be traced back to the empirical horizontal force factors adopted in the 1959 SEAOC Blue Book [16,17].Inappropriate selection of R-factor may lead to underestimation or overestimation of design forces in the structure.For demonstration, the hazard maps developed by S. Khan and Waseem [10] suggested a maximum PGA of 0.51 g near the city of Quetta for a 975 year return period, rendering an elastic base shear demand of 0.51 W for a bridge having a fundamental time period of 1.0 sec.Considering the maximum design base shear of 0.06 W for the same bridge, as per the suggestion of WPHC, this requires response modification factor R (demand) = 0.51/0.06= 8.50.Such a high value of a response modification factor (demand) may exceed the available response modification factor (capacity) of existing bridges.Nearly two decades ago, in the devastating earthquake that happened in Pakistan (the 8 October 2005 Kashmir earthquake), several bridges were either completely/partially collapsed or experienced extensive damage [1].Many researchers have quantified this factor for appropriate seismic analysis and design.
In addition, many researchers have developed mechanical devices, isolators, and special ductile connections at beam columns or in bridge piers to dissipate seismic energy (R-factor) within the structure, such as in the USA, where many bridges reach their intended life rather than being replaced with new construction [18].The authors proposed an apparatus, "Response Modification Apparatus", which mechanically exhibits and amplifies adequate stiffness and damping of bridges.Gastineau et al. [19] likely further improved the Response Modification (RM) Apparatus for flexibility in member and in overall performance during vibrational effects.Routledge et al. [20] reported the case study of the Wigram-Magdala Bridge in New Zealand, which was difficult to retrofit due to the Christchurch earthquake.They concluded the performance of plastic hinge regions in typical reinforced concrete structures before developing low-damage ductile jointed systems to control damage in plastic hinge regions.Warn et al. [21] analyzed the twenty bilinear bridge isolator responses of various bridges to design a new seismic isolator for the seismic resilience of US bridges.Chavdar et al. [22] reported the response modification of the bridge, which is isolated by a lead rubber bearing pad (LRBs) due to change in low ambient temperature from 20 °C to −30°.Zhang et al. [23] followed up the seismic response of fault crossing simply supported precast bridge by fault rupture.The amplitude of numerical models were evaluated for permanent ground displacement.Correspondingly, Bergami et al. [24] probed the near-fault event for the bridge and evaluated its response by 'Incremental Modal Pushover Analysis' (IMPA).China has rapidly progressed its high-speed railway (HSR) and lies in an earthquake-prone region.The HSR has faced seismic events ever since Guo et al. [25] proposed isolation bearing rather than common bearing to perform better in intense seismic events and working on the principle of "Equivalent energy-based design procedure (EEDP)".
Seismic codes worldwide suggest a different value for ductility factor, e.g., a ductility factor value from 3 to 4 is given in the US seismic codes, a value of 3.0 is given in the seismic code of Japan, a value of 6.0 is given in the seismic code of New Zealand, and a value of 3.50 is given in the seismic code for Europe.The AASHTO recommends a response modification factor of 3.0 and 5.0 for longitudinal and transverse directions, respectively.Kappos et al. [26] analyzed existing reinforced concrete bridges with yielding piers and others having elastomeric bearing pads placed between piers and superstructures.For the first class of bridges, the response modification factor was quantified as the product of ductility and overstrength factors.This factor varied from 4.2 to 10.1 and 3.7 to 11.6 in longitudinal and transverse directions.For bridges with elastomeric bearing pads, the response modification factor was defined as the ratio of spectral acceleration (corresponding to the pertinent predominant period of the bridge) that causes failure to design spectral acceleration.This factor varied from 4.0 to 6.6 and 4.3 to 9.3 in the longitudinal and transverse directions, respectively.They expressed that this factor given in seismic code is possible but, in numerous cases, underestimates the realistic deformability and energy dissipation capacity of the bridges.Zahrai et al. [27] investigated five existing reinforced concrete bridges in Iran with elastomeric and lead rubber bearings and proposed that the average response modification factor for bridges with elastomeric bearing pad isolators is 3.89 and 4.91 in the longitudinal and transverse directions, respectively.For bridges with monolithic deck-column construction, response modification factors of 2.92 and 2.41 were obtained for longitudinal and transverse directions.
In the context of Pakistan, Waseem and Spacone [28] carried out the vulnerability of three RCC bridges in northern Pakistan.They derived the fragility functions of bridge components.However, R-factors were not evaluated for the bridges.S. Ali [29] evaluated the R-factor for the low-cost single pier of the bridges.This study is limited to the pier's response only.Therefore, an attempt is made in this work to study the R-factor numerically for three bridges.This paper attempts to evaluate the R-factor values of bridges with ERB isolators to compare the results with those prescribed by the bridge design codes.According to AASHTO for seismic isolation design, the R-factor for all substructure elements of isolated bridges should be expressed in the AASHTO standard provisions for highway bridges.To investigate the R-factor for the most commonly found bridges in Pakistan, three simply supported spans in real RC bridges equipped with elastomeric bearing pads are analyzed using the static non-linear pushover analysis based on the Uang [30] approach.

Description of Reinforced Concrete Bridges
In Pakistan, the most common type of bridges constructed are simply supported superstructure multi-span RC bridges.In these bridges, the superstructure consists of an RC slab and prestressed girders.The substructure consists of circular single or multiple piers' bent.The superstructure and superstructure are separated by elastomeric bearing pads.Expansion joints are provided with every span and at abutments.Three real bridges are considered in this work as candidates and representatives of the bridges from Pakistan.

Mingora Bridge
Mingora Bridge is located in the Khyber Pakhtunkhwa province in the north of Pakistan.It is a three-span RC bridge, and the two bents consist of three pier shafts and a rectangular cap beam described in Figure 1.(layout plans were provided by "AA Associates" consultant firm).The abutments of the bridge consist of RC walls with a pile cap foundation.The abutment of the bridge has a seating arrangement for the bearing pad.Girder is pre-stressed, I-shaped, and placed over the bearing pads.Expansion joints of 25 mm are provided at every bent and abutment.Shear keys are provided between girders at the top of the cap beam to prevent lateral movement.A gap of 10 mm is provided between the girders and shear key, filled with polystyrene.Girders are connected with a diaphragm, and three diaphragms are provided in each span.The concrete compressive strength of girders is 34.47 MPa (5000 psi), while for all other members are 27.58MPa (4000 psi).The bearing pads used in this bridge are rectangular with a laminated layer of steel.Its dimension is 400 × 500 × 46 mm.The details of the Mingora Bridge are shown in Figure 1.

Shabqadar Bridge Charsadda
Shabqadar Bridge, located in the Charsadda district in the province of Khyber Pakhtunkhwa, is a ten-span, simply supported bridge with four multiple pier shaft bents described in the drawing (provided by "ABM Engineers" consultant firm).The bent possesses a solid circular column and a rectangular transom (cap beam).The abutment consists of an RCC wall supported on a foundation of pile cap and side wing walls.The abutment wall has seating arrangements for the bearing pad.The post-tensioned, I-shaped girders are supported at each end on the rectangular elastomeric bearing pad at a constant distance of 2600 mm.Five rectangular diaphragm beams are connected to hold four girders at each span, having a girder depth of 1800 mm.The symmetry of the members are the same throughout the bridge; each pier shaft is 1200 mm in diameter.The cap beam has a dimension of 1800 mm × 1200 mm.Deck slab, girders at bent, and girders and abutment walls are separated by an expansion joint of 25 mm.The bearing pad is placed on a concrete pad, and a special epoxy grout helps to cement the pad.A shear block on the top of each column bent and the abutment is also built to prevent the lateral load and gap of 10 mm between the girder and shear block filled by polystyrene.The prestressed girder concrete cylinder compressive strength is 34.47 MPa (5000 psi), and the compressive strength of concrete in other members (deck slab, diagphram, and pile) is 27.58 Mpa (4000 psi).Moreover, the compressive strength of the railing post, kerb, plank, and pilecap is 20.68 MPa (3000 psi).The laminated elastomeric bearing pad has internal steel plates in laminated pads according to ASTM A-238, and has a dimension of 400 × 600 × 69 mm.The detail of Shabqadqr Bridge is shown in Figure 2.

Dhoor Bridge Mansehra
The Dhoor Bridge is located in the district of Haripur in the province Khyber Pakhtunkhwa.It is a multi-span bridge with each span having the dimensions 30 × 100 mm, as shown in Figure 3 (layout plan drawings provided by "AA Assosiates" consultant firm).The superstructure consists of a post-tension girder holding deck slab.Four girders are rested on the elastomeric laminated bearing pad.They have nine intermediated bents with four solid circular pier shafts and rectangular transoms (cap beams) 2000 × 1200 mm.Every bent has a pile cap foundation.The bearing pad dimension is 500 × 500 × 82 mm, and pad thickness is 82 mm.It has six 3 mm steel plates embedded in five 12 mm rubber plates and 2 mm at the ends.Three rectangular diaphragm beams are constructed to hold four girders in the transverse direction.Pier shafts have 1200 mm dia.There are four girders in each span, spaced 2500 mm throughout the bridge.Each pier shaft has a 1200 mm dia.The cap beam has a dimension of 1200 mm × 2070 mm.There is a 25 mm expansion joint between the deck slab, girders, and abutment wall.The bearing pad is placed on a concrete pad with the help of a special epoxy grout cement pad.A shear block on the top of each column bent and an abutment is also built to prevent the lateral load and gap of 10 mm between girder and shear block filled by polystyrene.The concrete used in prestressed girders have 34.47MPa (5000 psi) cylinder compressive strength, and other members (deck slab, diagphram, and pile) have The general properties of these three bridges are shown in Table 1.

Computational Models
Meshing is one of the significant components in the computational modeling of structures to find precise outcomes from an FEA numerical model.The elements in the meshing essentially take numerous aspects into account to discretize stress/deformation gradients precisely.Characteristically, the lesser the extent of meshing, the more correct it will be across the real-time problems.The potential bottleneck of this simulation is that the more sophisticated the precision, the larger the simulations' calculations become, and thus solve times are prolonged.There is no worth in outlaying additional hours simulating with a compact mesh if a coarser meshing will give us the desired outcomes.Researchers frequently accomplish convergence studies to achieve the optimum equilib-rium between exactness and time.In this article, we carry out the parametric models of our bridge structures and find out in the convergence study of elements in mesh convergence corresponding to mid-span deflection of the bridge computational model, that increasing the mesh density beyond (1000) elements will not significantly diverge corresponding to deflection, therefore we limit of mesh element density to (1000) in order to optimize our result and reduce the computational time to solving numerical models as shown in Figure 4.
Three-dimensional (3D) models are generated because they can capture the mass and geometric distributions of the structure.The axial load variation and pier response in the transverse direction is captured more accurately in 3D models [31].Models of the bridges are generated in the finite element software package.Frame elements are used to generate the piers, girders, diaphragms, and cap beam components of the bridges.Link elements are used to generate the responses of the bearing pads and shear key [32].Expansion joints are modeled with gap elements and deck slabs with shell elements.The piers, abutments, expansion joint, shear key, and bearing pads are non-linear components of the bridge, while the superstructure is a linear elastic component.Computational modeling of the Mingora, Shabqadar, and Dhoor bridges are shown in Figure 5a-c, respectively.

Material Properties Modeling of Bridge
The ductile behavior of a plastic hinge is significantly affected by the non-linear material property used to define the frame member receiving the hinges.The material strength and stress-strain (σ − ε) relation used for concrete was unconfined and confined using the Mander model and reinforcing the steel model to precisely capture the bridge's capacity and behavior.The properties of cement concrete conformed to Section 3.2 of SDC 2004 [33].The Mander et al. (1988) [34] model represents the uniaxial stress-strain behavior of unconfined and confined concrete, as shown in Figure 6.The tensile stress is measured by ACI 318 as ƒr' = 0.6 ƒc′ (MPa) for concrete having normal weight, with an initial modulus of elasticity Ec according to Section 3.2 of SDC 2004.A moment curvature (M-φ) analysis is to be carried out for the concrete column (see Section 3.4.2); the steel longitudinal and transverse reinforcement properties are to be used according to Sections 3.2.2 and 3.2 of the SDC 2004 Guidelines for Steel ASTM A-706.The steel material model in tension and compression regions assumes an initial elastic behavior to yield strength rather than a yield strength, tracked by a strain-hardening section.According to SDC 2004, Section 3.2, the maximum expected values of yield stress (fye) and ultimate stress (fue) for all bar sizes are to be taken as 448 and 655 MPa, respectively, as shown in Figure 7.

Finite Element Modeling of Bridge
A numerical model of the bridges was developed in the finite element analysis program.The properties used in generating the models are listed in Table 1.The contact (constrain) properties are used for generating the computational model.For example, frames, girders, and piers etc. are modeled by a rigid link element.Shell elements were used to model the slab, while bi-linear link elements were used to model bearing pads and abutments.A gap element was used to model expansion joints.These numerical models of the bridges are shown in Figure 8a-c, respectively.The characterization of the boundary/edge conditions in a structural arrangement is a significant feature in the accumulation of its stiffness matrix, thus affecting both the static and dynamic performance of the structural system.The boundary conditions must be allocated appropriately through simplified and realistic models of the bridge abutments and footing system to properly estimate the ductility capacity and seismic demand on key structural sections.In a dynamic analysis of the bridge, the modal periods and mode shapes, as well as other related properties, are greatly affected by such assignments.In this article, bridge footing is constrained by a fixed support system.Bridge abutment, bearing bad, and shear key ponding is discussed in detail in the later sections.

Bearing Pads
Bearing pads are modeled as link elements with non-linear constituent laws.There are two possible failure modes of the bearing pads: (1) failure in shear between the neoprene and steel reinforcement and (2) failure in sliding between the neoprene and contact surface of the overlaying girder.According to Tortolini and Annunzio, [35] shear failure force varies from 1.2 to 1.8 times G (shear modulus).The seismic design criteria (SDC) of the California State Department of Transportation CALTRANS-2006 [36] recommends the ratio between the normal (dead load) force and horizontal failure force (sliding force) acting on the bearing pad to be 0.40.
In this work, the sliding failure of the bearing is considered, and an elastoplastic force-displacement relationship is assumed for the bearing pads shown in Figure 9. Since the link element has six degrees of freedom, the rotational degree of freedom is assumed as zero.The translational stiffness values are calculated using the relationship given by Equation (1) proposed by [37].Table 2 presents the properties of bridge-bearing pads.
where, Kh = horizontal stiffness of the bearing pad G = shear modulus of neoprene rubber, Ab = area subjected to shear; and H = the total thickness of the pad excluding steel plates.In the vertical direction, the stiffness is assumed as 100 times Kh, shown by Equation (2) as follows: where Kv = horizontal stiffness of the bearing pad, and G = shear modulus considered 1 MPa.Mathematically shear failure force and shear stress in the bearing pad is expressed by Equations ( 3) and ( 4) as follows: where Vf = shearing failures subjected to the area of bearing Ab, and τ = shearing stress.

Abutments
The abutments of all three bridges are seating abutments, and they are considered the non-linear component of the bridge.Their force-displacement relationship is modeled as elastic, perfectly plastic based on the recommendations of Caltrans (2006) [36] as shown in Figure 10a,b.Caltrans proposed the maximum passive pressure (Pbw) given by Equation ( 5).The stiffness of the abutment is computed by Equation ( 6).Ki is the initial stiffness, and its value is taken based on the test of large-scale abutment carried out at the University of California Davis (i.e., 20 kips/in/ft.or 11.5 kN/mm/m proportioned to the height of the back wall) [38].
W = The back wall/diaphragm width for seat and diaphragm abutment, Ki = the initial stiffness value.

Expansion Joint
A gap, ponding, or expansion joint can alter the response of a bridge in a seismic event.We can model in (Sap2000) by connecting link elements by a linear or bilinear spring model as shown in Figure 11a,b.It is challenging to depict the exact behavior of these springs.Susender et.al [39] observed at the bridge's gap and device a method to model it by a stereo mechanical method; during shaking and drift, seismic energy is dissipated and expressed in Equation (7).

Plastic Hinges, Performance Limit and Bridge Piers
The lump plasticity approach is been carried out in this research.Plastic hinges are assigned to piers of bridges.These hinges are provided at the top and bottom of every pier.Hinge lengths (Lp) are computed by Equation (9) given in (Caltrans SDC, 2013) [40].
where, L is the pier; fye denotes the expected yield strength of longitudinal reinforcement, and 'dbl' represents the nominal bar diameter of longitudinal reinforcement.Plastic hinges are assigned at the mid-height of the plastic zone with a segmental length L = Lp as shown in Figure 12.The plastic hinge's rotation directs the capacity of a member to endure inelastic distortion and is used in SAP000 to define column plastic hinge properties.(FEMA-356, 2000) provides a generalized F-D relation model shown in Figure 13 for the static non-linear analysis procedure, and the auto model in SAP2000 (PMM hinge in SAP2000) for the axial load moment hinge.Elastic behavior from initial position "A" to significant yielding point "B" defines the F-D curve.At that point, the stiffness cut from point "B" to "C"."C" point has a confrontation equal to the nominal strength, then an abrupt lessening in the cross-load resistance to the retort at reduced resistance "D" to final loss of resistance "E".The slope of "BC" is generally achieved in "0 and 10%" of the initial condition.The "CD" slope links to the preliminary failure of the section.The "DE" slope depicts the remaining strength of the section.These slopes are referred to (FEMA-356, 2000) in order to limit the "hinge rotation" behavior of RC Sections.The slope lines "B" and "C" represent acceptance criteria for the hinges rotation, which is referred to as immediate occupancy (IO), life safety (LS), and collapse prevention (CP) [41].

Performance Limits
Two categories generally define performance limits or acceptance criteria: (i) global structural performance limits and (ii) local member/elements level performance limits.(FEMA-356, 2000) [26,27].The global-level acceptance criteria depict the capacity of the structure to endure gravity and lateral loading.The lateral load capacity of structural performance should not be degraded by 20% of its overall strength.Lateral displacements should verify against standard limits (ATC-40, 1996) as given in Table 3.The local member/elements level performance limits are often imposed by nonstructural and component damage in structural elements, such as beams and columns.They are constructed on the plastic hinge rotation capacities limit.Table 4, below, states the "deformation limits" referred to by (ATC-40, 1996).Hence, cross-verification of the failure of elements by "flexural demands" and "shear failure" does not befall earlier than these rotation limits are reached.The capacity of all ductile concrete structure sections for non-linear plastic hinge/plastic moments shall be considered by M-ϕ plot dependent on anticipated material properties.M-ϕ plot derives the curvatures linked with a range of moments for an x-section depending on the principles of strain considerations and balance of forces.The moment curvature (M-ϕ) plot can be idealized (bi-linearization) with an elastic, perfectly plastic response to approximate the plastic moment capacity of a member's x-section.The elastic region of the idealized (bi-linear) curve must pass through the point, indicating the first reinforcing bar yield.The idealized plastic moment capacity is attained by corresponding to the areas between the actual and the idealized M-ϕ plots beyond the first reinforcing bar yield point, as shown in Figure 14.Moment curvature analysis for piers is designed and analyzed in SAP2000 as shown in Figures 15-17; maximum axial dead load is carried out to examine the yield moment and consider the yield and ultimate curvature rotation for the three bridges, and in Equations ( 10)-( 12), these parameters are stated for the Mingora, Shabqadar, and Dhoor Bridges, respectively.

Response Modification Factor (R-Factor)
The main objective of earthquake-resistant design is: (1) In the event of minor earthquakes (less damaging), the structural response is in the elastic region.(2) In case of a significant earthquake, the structure can resist and sustain the overall stability.In major earthquake events, structures respond in the non-linear range.Therefore, non-linear analysis is required to describe the behavior of structures.However, due to the simple elastic procedures, traditional analysis procedures are based on linear analysis suggesting the reduced seismic forces.The schematic base shear-displacement response of a traditional structure is shown in Figure 18.The response modification factor (R) is a significant feature in building seismic design.Equivalent statistical analysis (ESA), which is recurrently used to evaluate the seismic response of structures, can be executed by shaping the R-factor.In actuality, R-factor directs the capacity of a structure to dissipate energy in an inelastic manner [42,43].
As mentioned before, R-factor is a significant parameter in designing vulnerable seismic design.During an earthquake event, it considers the level of inelastic behavior of the structure.R-factor is defined by (NEHRP 1988) [44] as the "factor intended to account for both damping and ductility inherent in structural systems at the displacements great enough to approach the maximum displacement of the systems" [44].This definition delivers an intuition into the consideration of the seismic behaviors of the structure and the predictable response of code-based structures in the design earthquake.R-factor decreases the design forces in seismic resilient structures based on the damping, energy dissipation capacity, and overstrength of the structure [45].
This reduced seismic force is extracted by dividing the elastic base shear by (R-factor).In another sense, in the linear analysis [27], the elastic base shear force level (Ve), is reduced to the force level corresponding to the appearance of the first plastic hinge in the structure (Vs).

Uang's Procedure R-Factor Coefficient Method
Uang's procedure [30] is one of the methods used to compute the R-factor.The extreme base shear is considered when the structure remains in the elastic zone.The R-factor comprises a reduction factor due to ductility (Rµ) and an overstrength coefficient (Ω).Rµ is the proportion of base shear at the elastic level (Ve/Vmax) and to the base shear at the level of structural failure (Vy); and Ω is definite as the fraction of the base shear at strength level, which corresponds to the formation of a yield mechanism (Vy) to the base shear of the structure when the primary plastic hinge is formed (Vs).Conferring to Figure 16, the mandatory elastic resistance due to the base shear coefficient is determined by Equation (13).
where, W is the effective structural weight, and Ve is the maximum base shear when the structure is in the elastic range.Regarding the Uang's method, the required equations for calculating R-factor are as follows: The structure ductility, µ, is defined in terms of the maximum structural drift (Δmax) and the displacement corresponding to the idealized yield strength (Δy) as in Equation ( 17): Many researchers have deliberated the two key components of R-factor presented in Equation ( 16).The ductility-dependent component, Rµ, has received substantial consideration.Ductility reduction factor Rµ is a function of the structure's features, including ductility, damping, fundamental period of vibration (T), and the features of earthquake ground motion [46], which presented a relation for Rµ in the following form: where In Equation ( 19), α is the post-yield stiffness given as a percentage of the initial stiffness of the system, and a and b are parameters given as functions of α that can be obtained from Table 5 [46].Alternatively, we can also attain "Rµ" corresponding to "time period" by help of the chart as shown in Figure 19, which is proposed by Nassar and Krawinkler [47] instead of using Equation (18).Pushover analysis is a very popular technique used for seismic assessment of the structure.In this procedure, a gravity load is applied to the structure and a load vector is applied sideways, and the structure is pushed until an ultimate limit is reached.The fundamental output of pushover analysis is a force versus displacement curves, also known as the capacity curve.Capacity curves provide information about the ductility and strength capacity of structures.In this work, the displacement coefficient method pushover analysis documented in (FEMA-356, 2000) [48], is used.
The non-linear static procedure (NLSP) is considered as just an analog of the linear static procedure (LSP) (i.e., that the static later loading is imposed on non-linear models of the structure).The monotonic pushover analysis (POA) is the most important tool of all present nonlinear static procedures (NSP).The primary working principle of pushover analysis (POA) is that after applying a gravity load on the non-linear model of the considered structure, the lateral load pattern (simple distribution of seismic forces) is applied, as shown in Figure 20.This procedure progresses until that structure becomes unstable or deteriorated, or the pre-assumed targeted deformation is touched [49].The key objective of pushover analysis (POA) is to evaluate overall lateral strength, deformation ductility, and the deterioration mechanism of the considered structure with the influences of seismically induced lateral forces on the structure.At the global structural level, the monotonic/pushover or (capacity curve) of the detailed numerical structure model straightaway comprises the predicted nonlinearity of structural members.The maximum seismic displacement and other induced parameters can be obtained via any recognized non-linear static procedure.Then, these pushover curves are idealized (bilinear relation) in any simple form or relation for analysis simplicity.Comparing both linear static procedure (LSP) and nonlinear static procedure (NSP), the latter advances over the former procedure based on its ability to redistribute the elemental internal forces of the structure that constantly undergo incremental seismic lateral loads.This created a vivid understanding of how structures respond inelastically and achieved several limit states when undergoing an inelastic zone.The existing non-linear static procedure (NLSP) is designed on the principle of a detailed structural model of structures converted into an idealized single degree of freedom system (SDOF).The global structure response of the pushover curve can be idealized in such a manner that it depicts the nonlinear force-displacement (F-Δ) relation.This force-displacement (F-Δ) relation is assigned to the idealized (SDOF), which is then predicted to behave as an actual non-linear detailed structural model.This procedure is explained in detail in Figure 21 [49].

The Capacity Spectrum Method (CSM)
The primary nonlinear static procedure (NSP) given in the ATC 40 (1996) [50] report entitled "seismic evaluation and retrofit of concrete buildings", is working on the equivalent linearization approach (the capacity spectrum method, CSM), which was first advanced as a rapid evaluation method for a pilot seismic risk project of the "Puget Sound naval shipyard for the U.S. navy in the early 1970s" [49].In this procedure, the capacity or (NSP) curve is changed into the "capacity spectrum" using point-by-point alteration of the coordinates of the modal spectrum.The first-mode inertia load vector generates this pushover curve.The base shear (Vb) and roof displacement (Droof) of the capacity curves are transformed into the conforming spectral acceleration (Sai) and spectral displacement (Sdi) values on the capacity spectrum using the Equations ( 20) and (21).
where W presented "seismic weight" (i.e., the sum of dead load and the likely live loads), α1 and PF1 are the "modal mass coefficient" and "the modal participation factor for the first vibration mode respectively", and ϕ1,roof is "the roof level amplitude of the first vibration mode shape".Every point on a capacity spectrum curve is associated with a unique spectral acceleration (Sa), spectral velocity (Sv), spectral displacement (Sd), and natural period (T).This reduced response spectrum is also converted to the accelerationdisplacement response spectrum (ADRS) format using Equation (22).
The modified response spectrum in the ADRS format is called the "demand spectrum".This demand spectrum is overlaid with the capacity spectrum to find their intersection point, which corresponds to a condition for which the seismic capacity is equal to the demand imposed on the structure.This point is called the "performance point", which estimates the actual maximum displacement expected during an anticipated earthquake.The overall procedure is shown in Figure 22.

Displacement Coefficient Method (DCM)
This procedure uses the co-efficient to convert displacement of the linear to displacement of the non-linear system.This method is adopted by various codes, e.g., FEMA 356 [48], (FEMA-440) [51], and EC8. Figure 23 provides information about this procedure.This method is carried by modifying the elastic response from the SDOF equivalent with coefficient factors C0, C1, C2, and C3 [48].
where, Te is the effective fundamental period, C0 is the modification factor for the spectral displacement of an equivalent SDOF system to an MDOF system, C1 is the modification factor to relate the expected maximum inelastic displacement to the displacement calculated for the linear elastic response, C2 is the modification factor that represents the effect of pinched hysteretic shape, stiffness degradation, and strength deterioration on the maximum displacement response, C3 is the modification factor that represents increased displacement due to dynamic P−δ effect.Sa is the response spectrum acceleration in g (i.e., acceleration of gravity).Figure 24 shows the performance point according to the displacement coefficient method [41].The first step is to determine the target displacement (δT).The procedure starts by setting an effective fundamental period that considers inelastic conditions.The effective fundamental period represents the linear stiffness of the equivalent SDOF system and corresponds to maximum spectral acceleration (Sd).Then, target displacement can be specified using Equation (23).

Response-Factor Computation
A block diagram of the process for computing the bridge models' R-factors is depicted in Figure 25.The capacity curves obtained by the pushover analysis depicted in Figure 23 and the parameters depicted in Figure 18 are used to find Ve, Vy, and Vs.Based on the elastomeric bearing pad parameters in Table 2, non-linear static analyses performed on the bridge models with the ERB isolator are also computed.The response modification factor is obtained by displacement coefficient method (DCM).To evaluate the force reduction factors of concrete bridges at the ultimate limit state, three bridges that cross the three regions of the KPK province of Pakistan, Mingora, Shabqadar, and Thrace, were selected.The basic characteristics considered in the classification are the type of deck, type of piers, and type of pier-to-deck connections, as mentioned in Table 1.All three of the selected structures belong to the category defined in the previous Section 2, which is a deck supported through elastomeric bearings on elastically responding piers.The selected bridge's main features are given in Figures 1, 2, and 4. The pushover curves derived using analysis with SAP-2000 point hinge models (lumped plasticity), as mentioned in the previous Section 3, were idealized as bilinear curves (Figures 26-31) to define a conventional yield displacement (dy) and ultimate displacement (du) by the method mentioned in a previous Section 4.2.2.The derived overstrength factors for bridges with yielding piers and the ductility factors for the same bridges are given in Tables 6-11.It is well known that both Ω and Rµ range rather sketchily; taking the lowest among the values calculated for the longitudinal and the transverse directions in each bridge, Ω varies from 1.45 to 2.17, and Rµ from 2.96 to 3.74 as mentioned in Tables 12 and 13.It should also be noted that high Rµ values do not essentially agree with high Ω values.The R-factor's calculated values, overstrength, and ductility factors in longitudinal and transverse directions are listed in Tables 4 and 5. R-factor in the longitudinal direction is 4.50, while in the transverse direction, it is 5.18.A comparison of R-factor in longitudinal and transverse directions is shown in Figures 32 and 33, respectively.These values are more than recommended as per AASTHO-LRFD, indicating the greater capacity of the bridges.The higher value of the R-factor in the transverse direction is as in AASTHO-LRFD.In the transverse direction, bridges are more resilient than in the longitudinal direction due to the presence of shear keys or the contribution of the substructure.

Figure 9 .
Figure 9. (a,b) Bearing pad modeling and F-D relation in transverse direction and in (c,d) longitudinal direction.

Figure 13 .
Figure 13.Force deformation relation (back bone curve) of plastic hinge.

Figure 19 .
Figure 19.Strength reduction factor proposed by Nassar and Krawinkler.

Figure 25 .
Figure 25.Flow chart for computing R-factor.

Figure 33 .
Figure 33.Comparison of R-factor in transverse direction.

Table 1 .
Properties of bridges.

Table 2 .
Properties of bridge's bearing pads.

Table 12 .
Computed R-factor in longitudinal direction.

Table 13 .
Computed R-factor in transverse direction.