Abstract

Isolated bridges are commonly designed in the near-fault region to balance excessive displacement and seismic force. Optimal intensity measures (IMs) of probabilistic seismic demand models for isolated bridges subjected to pulse-like ground motions are identified in this study. Four typical isolated girder bridge types with varied pier height (from 4 m to 20 m) are employed to conduct the nonlinear time history analysis. Totally seven structure-independent IMs are considered and compared. Critical engineering demand parameters (EDPs), namely, pier ductility demands and bearing deformation along the longitudinal and transverse directions, are recorded during the process. In general, PGV tends to be the optimal IM for isolated bridges under pulse-like ground motions based on practicality, efficiency, proficiency, and sufficiency criterions. The results can offer effective guidance for the optimal intensity measure selection of the probabilistic seismic demand models (PSDMs) of isolated bridges under pulse-like ground motions.

1. Introduction

It has been observed that ground motion recordings in the near-fault region, which is usually taken as within 20 km from the fault rupture, tend to differ from far-field rivals in three main characteristics, that is, velocity pulse, directivity effect, and large vertical acceleration component [1]. Some of these differences are usually shown in the velocity and displacement time series, which have been attributed to two effects: the rupture directivity effect and the fling step. The rupture directivity denotes the effect of rupture propagation relative to the site [2]. The near-fault directivity often manifests in the fault-normal direction, which is the direction perpendicular to the surface of the fault rupture. A site may experience forward directivity when the rupture propagates towards the site with a velocity almost equal to the shear-wave velocity of the soil. This effect appears in the form of a large, long-period velocity pulse. And the fling step effect is observed in the fault-parallel direction, that is, the direction paralleled to the fault slip, which appears in the form of a permanent displacement. Due to the above features, the near-fault ground motions impose high input energy on structures and result in severe damage. This phenomenon was originally discovered by Bertero et al. [3, 4]. However, it is only after the 1994 Northridge earthquake that the severe effect of near-fault ground motions on the demand of structures and the necessity of incorporating the effect into the design process were recognized [5, 6].

The probabilistic seismic demand models conditioned on a single intensity measure (IM) of ground motions are a common and necessary step in generating the analytical fragility curves. The intensity measures (IMs) represent characteristics of ground motions and play an important role in the probabilistic seismic demand analysis. Tracking and reducing uncertainty associated with the PSDM can be accomplished by selecting an optimal IM upon which the demand models are conditioning [7].

Researchers have done a lot of work on the selection of optimal IMs for far-field ground motions in the recent years. The Modified Mercalli Intensity Scale was chosen by the Applied Technology Council [8] as the preferred IM. These data were adopted in the risk assessment software package HAZUS. The recent versions of HAZUS selected the spectral acceleration at a period of 1 s (Sa−1) and peak ground displacement (PGD) as the IMs. Padgett et al. [7] investigated the seismic response of the multispan simply supported steel girder bridge class and found that peak ground acceleration (PGA) and spectral acceleration at the fundamental period (Sa) tended to be the most optimal IMs for synthetic ground motions and cumulative absolute velocity was also a contender for recorded motions. In addition to the common IMs, more complex parameters, multiparameter, or vector-based IMs have been proposed [912]. Within the bridge engineering community, PGA and Sa are the widely accepted IMs [7, 13].

However, there is relatively less research on the optimal intensity measures of pulse-like ground motions [14]. Yakhchalian et al. [15] proposed that the commonly used scalar intensity measure, namely, the spectral acceleration at the fundamental period of the structure (Sa), is deficient and insufficient to predict the seismic response with respect to the pulse periods of pulse-like ground motions and developed an optimal integral-based IM (I-Sa) with 15 generic frame structures. Baker and Cornell [16] demonstrated that the near-fault pulse-like ground motions are not well described by traditional intensity measures and proposed a vector intensity measure that combines Sa with a measure of the spectral shape. Li et al. [17] proposed two fuzzy-valued (FV) structure-specific intensity measures (one is based on the squared spectral velocity and the other is based on the inelastic spectral displacement) to characterize near-fault pulse-like ground motions. Zhong et al. [1] investigated the response of cable-stayed bridges under synthetic and recorded pulse-like ground motions and found that PGV tends to be the optimal intensity measure compared to other six structure-independent IMs.

In general, the existing IMs can be classified into two categories, that is, structure-dependent IMs and structure-independent IMs. The former combine structure characteristics with intensity measures, while the latter only consider the features of ground motions. However, Padgett et al. [7] proposed that structure-dependent IMs are difficult to obtain in a regional risk assessment for bridge portfolios due to the lack of sufficient information such as fundamental periods of structures. Isolated bridges display totally different seismic response compared to the ordinary short period bridges. Consequently, this study aims to identify the optimal IM for isolated bridges under pulse-like ground motions by evaluating seven existing structure-independent IMs. Four typical isolated bridge types with varied pier height (from 4 m to 20 m) are employed to conduct the nonlinear time history analysis. Efficiency, practicality, proficiency and sufficiency, and hazard computability are selected to evaluate the candidate intensity measures.

2. Fragility Function Methodology

As an emerging tool, the seismic fragility is a useful tool for evaluating the potential seismic damage of bridges during the earthquakes. Seismic fragility refers to the conditional probability that the seismic demand meets or exceeds the seismic capacity of structures subjected to a series of specific ground motions, generally expressed in equation (1):where Pvulnerability is the conditional failure probability, SD and are the median and dispersion of seismic demand, respectively, SC and βc are the corresponding median and dispersion of seismic capacity, respectively, and is the cumulative normal distribution function. In this paper, probabilistic seismic demand analysis (PSDA) is utilized to derive the fragility curves of bridges, which can establish the probabilistic relationship between the engineering demand parameters (EDPs) of components and the ground motion intensity measures (IMs) [18]. During the process, the peak demands (Di) of the critical components, namely, the curvature ductility of the pier and deformation of the bearing, are recorded for the ith ground motion.

If it is assumed that the seismic demand coincides with the lognormal distribution, the relationship between the median of seismic demand (SD) and the intensity measure (IM) can be written in equation (2):where a and b are regression coefficients obtained from the response data of the nonlinear time history analysis. The dispersion is theoretically related to ground motions. But, for the sake of simplification, it is usually assumed that is constant [19] as follows:where n is the number of ground motions.

3. Characteristics of an Optimal IM

Selection of an approximate IM plays a vital role in the probabilistic demand seismic models, which provides the link between EDPs and IMs. An optimal IM should satisfy the following five requirements [1, 7].

3.1. Practicality

Practicality is used to examine the correlation between the demands placed on the structure and the intensity measures of ground motions, which is measured by the regression parameter b in equation (2). A larger value of b indicates a more practical IM.

3.2. Efficiency

Efficiency is known as the common metric that characterizes an optimal IM. An efficient IM can reduce the variation in the estimated demand median, which is represented by the lower regression parameter in equation (3).

3.3. Proficiency: A Composite Metric of Efficiency and Practicality

The conventional selection of an optimal IM based on practicality and efficiency may lead to a challenge in balancing the two parameters. To deal with the problem, Padgett et al. [7] proposed a composite measure ζ to combine practicality and efficiency, which is defined in equation (4). A more proficient IM yields a lower ζ.

3.4. Sufficiency

Sufficiency is an alternative measure to identify an optimal IM of PSDMs [19]. A sufficient IM is statistically independent of ground motion characteristics, such as the magnitude (M), the epicentral distance (R), and fault characteristics. The sufficiency of IMs is assessed with a regression analysis on the residual ε/IM to obtain the corresponding value, from the PSDM to the ground motion characteristics, M or R [1]. A larger value of the regression evaluation indicates a more sufficient IM.

3.5. Hazard Computability

Hazard computability is also a viable measure of an optimal IM. As defined in Giovenale et al. [20], the hazard computability is the level of effort to assess the seismic hazard or derive the hazard curve. PGA or spectral accelerations at the period of 0.2 s or 1 s are readily available in hazard maps; however, other structure-dependent IMs need even more effort and there is difficulty in calculating some IMs from hazard maps.

Therefore, only structure-independent IMs, namely, seven existing IMs in Padgett et al. [7], are considered and compared in this paper. The detailed IMs are listed in Table 1 with their formulas.

4. Characteristics and Modeling of Isolated Bridges

Four typical isolated girder bridge types that are common in transportation networks are adopted in this paper as case studies, namely, multicell box girder (MCBG), multibeam box girder (MBBG), concrete I-girder bridge (CIGB), and steel I-girder bridge (SIGB). To generate PSDMs that are representative of each bridge type, some generalizations of the configurations are necessary to be made. As-built bridge data are collected from the various departments of transportation and the advice of engineers who are experienced in designing and restoring bridges is taken to determine pier details, superstructure geometry, isolation bearing parameters, and so forth. The Guidelines for Seismic Design of Highway Bridges [21] are also referred to in order to ensure that the crucial characteristics (e.g., stiffness, mass, and natural period) are generally consistent with the existing bridges. In particular, the pier dimension is determined to guarantee that the axial compression ratio of piers is about 10% for the ductility demand. In this way, the conclusion obtained in this paper is representative of the isolated girder bridge types rather than a certain bridge (Figure 1).

To demonstrate the effect of structure periods on the optimized results, the pier height varies from 4 m to 20 m to generate a wide period range. The typical longitudinal layout of the bridge is given in Figure 2. The vertical reinforcement ratio of the piers is set as 1.2% uniformly with 16 mm circular ties spaced at 100 mm vertically. The configuration of the pier is also presented in Figure 2, where n is 7 for CIGB and MBBG, 9 for MCBG, and 5 for SIGB and D refers to the side length of the pier. The detailed geometry of each bridge type is summarized in Table 2.

As an effective measure for earthquake resistance and disaster prevention, isolation devices play a vital role in reducing the seismic response of structures. The common types of isolation devices applied in transportation networks include elastomeric bearings (ERB), lead-rubber bearings (LRB), and friction pendulum systems (FPS) [22]. In this paper, LRB is selected to be equipped in the bridge systems. The mechanical properties of LRB can be modeled with a bilinear relationship, which are decided with three parameters, namely, characteristic strength Q, elastic stiffness K1, and postyielding stiffness K2. The specific formulas of these parameters are given in Table 3, where Alead is the cross section area of the lead core, G is the shear modulus of the rubber, A is the area of rubber, and Σtr is the total thickness of rubber layers. Since the parameters rely on the specific site scenarios and there is an available LRB portfolio [23], the bearing parameters are determined directly according to the superstructure gravity, which are listed in Table 4.

The four bridge types considered in this paper include elements that would undergo highly nonlinear behavior under earthquakes, which are incorporated into the three-dimensional detailed nonlinear OpenSees models [24]. The superstructure usually remains elastic subjected to the earthquakes and is simulated with elastic beam-column elements. Due to the constraint effect of stirrups in the piers, we should divide the concrete into cover parts and confined parts, respectively. Therefore, we adopt fiber elements to model the substructures. Each fiber element has its own stress-strain relationship and can simulate the confined concrete, cover concrete, and longitudinal reinforcement, respectively. The response of isolation devices is simulated by a nonlinear bilinear element to account for the yielding of lead. Figure 3 gives the schematic of the analytical model of isolated girder bridges with the reinforcement, concrete, and LRB relationships, where σy, εy, fc, εuu, εcu, and K are yield strength, yield strain, compressive strength, strain at compressive strength for unconfined concrete, strain at compressive strength for confined concrete, and amplification factor due to constraint effect, respectively.

In the model, the yield strength, strain-hardening ratio, and elasticity modulus of steel are 400 MPa, 0.5%, and 2.0 × 105 MPa, respectively. With Concrete01 and Steel01 material models in OpenSees, a modified uniaxial Kent-Scott-Park [25] core concrete model is established considering the degraded unloading/reloading stiffness. The compressive strengths of core concrete and cover concrete are set as 29.8 MPa and 26.8 MPa, respectively, and the crushing strains are 0.016 and 0.0035.

In addition, the damping is taken into account for each sample in this paper. The common damping used for structure dynamic analysis includes viscous damping, modal damping, and structure damping [26, 27]. Although many damping models have been proposed, the ideal damping is still difficult to be determined. In this paper, viscous damping is taken into consideration with Rayleigh damping. Rayleigh damping has the characteristics of algebraic manipulation, uniform mathematic expression, and easy application, which is applied widely in structural dynamic response analysis. It is assumed that the damping matrix [C] is the combination of mass matrix [M] and stiffness matrix [K] [26]:where α0 and α1 are two proportionality coefficients calculated from equation (6):where ξ is the damping ratio (assumed as 0.05 in this paper) and ωi and ωj are the natural vibration frequencies of structures. Since both the longitudinal and transverse seismic responses of bridges are the focus of this paper, the first two-order longitudinal and transverse frequencies are selected as the values of ωi and ωj. The detailed calculation results of Rayleigh damping are given in Table 1.

5. Ground Motion Database

A database consisting of 243 pulse-like ground motions was established by Shahi and Baker [28] in the June 2012 version of the NGA-West2 database by rotating the ground motions and identifying pulses in all orientations. Dabaghi and Kiureghian [29] selected 121 pulse-like ground motions from the database excluding the aftershocks and the records at sites far from the fault (larger than 30 km). Therefore, the remaining pulse-like ground motions database contains 121 pairs in the 2014 version of NGA-West2 database. Figure 4 shows the distribution of pulse-like ground motions with respect to M and R, PGA, PGV, and the pulse period (Tp). The ground motions were recorded from earthquakes, where M ranges from 5.4 to 7.9 and R changes from 0.07 to 28.04 km. The mechanisms of the fault include reverse, strike slip, and reverse oblique. Each record in the pulse-like ground motion database is rotated into the direction containing the largest pulse amplitude and the corresponding orthogonal horizontal direction following the method proposed by Baker [2] and Dabaghi and Kiureghian [29].

In this paper, the nonlinear history analysis is conducted using the abovementioned 121 pairs of pulse-like ground motions. EDPs including curvature ductility demands (μx and μy: the subscripts x and y refer to the longitudinal and transverse directions), bearing deformation at the pier (δx, δy), are tracked and recorded during the process.

With the seismic demands and IMs of ground motions, PSDMs are easily developed. For illustration purpose, PSDMs for curvature ductility μx conditioned on PGA and PGV of MCBG (pier height is 18 m) are given in Figure 5, where R2 is the correlation coefficient. As illustrated in the figure, b values are 0.84 and 1.81, respectively, and are 1.42 and 0.93, which indicates that PGV is more practical and efficient. The corresponding ζ are 1.69 and 0.51, respectively, demonstrating that the more proficient IM is PGV. The sufficiency of PGA is examined with conditional statistical dependence on M and R, which is plotted in Figure 6. As shown in the figure, the values are 0.2527 and 2.3308 × 10−4 with respect to M and R, respectively. This phenomenon reveals that PGV is more independent of M and dependent on R.

To highlight the pulse effect on the seismic demand of structures, the optimal intensity measure selection process is repeated for far-field ground motions case. A set of 80 far-field ground motions in Shafieezadeh et al. [30] is used in this paper as the seismic excitation. Figure 7 shows the distribution of far-field ground motions. The ground motions were recorded from earthquakes with M ranging from 5.5 to 7 and R varying from 12.3 to 60 km.

6. Optimal Intensity Measure under Pulse-Like Ground Motions

6.1. Practicality

The practicality comparison of IMs can be estimated with the regression parameter b. A more practical IM yields a larger b. The practicality comparison of the four bridge types with different pier heights is shown in Figures 8 and 9. As illustrated in the figures, PGV tends to be the most practical IM for the pier curvature ductility in the longitudinal and transverse directions, followed by CAV. In contrast, for the lower height cases (from 4 m to 12 m), Iv tends to be the least practical IM, and for the higher height cases (from 14 m to 20 m), PGA is the least practical IM. For the whole height range, Iv, PGD, CAD, and PGA are less appropriate due to relatively low practicality parameter b. Similarly, for the displacement-related EDPs, namely, the bearing displacements, PGV is the most practical in the longitudinal and transverse directions, followed by CAV; the least practical IM is Iv. For the considered height range, Iv, CAD, and PGD tend to be less appropriate because of the relatively low practicality parameter b. In general, PGV tends to be the most practical IM, followed by CAV. In particular, one should notice that, for the long-period cases (e.g., 16–20 m), PGA is not a practical IM under pulse-like ground motions, which is a commonly used IM to establish the analytical fragility curves [22].

As mentioned above, the optimal intensity measure subjected to far-field ground motions is also identified in this paper. Table 5 summarizes the practicality comparison of the seven structure-independent IMs under far-field ground motions of MCBG. For simplicity, the comparison in the transverse direction is not given because it follows a similar tendency. The practicality comparison results of the other three bridge types are also not listed here due to a similar tendency. As shown in the table, CAV tends to be the most practical IM, followed by PGV; however, Iv tends to be the least practical. To demonstrate the pulse effect on the practicality comparison results graphically, Figure 10 shows the comparison results under pulse-like ground motions and far-field ground motions, where PGA, PGV, Iv, and CAV are chosen as the candidate IMs (conditioning on μy, the other cases follow the similar tendency). As shown in the figure, PGV can be chosen as the appropriate IM to assess the seismic demand under far-field and pulse-like ground motions. In addition, PGA is less appropriate to develop PSDMs with respect to practicality, especially for higher piers (larger than 16 m) under pulse-like ground motions. For example, b values are 0.69 and 0.35 for 8 m and 16 m cases subject to pulse-like ground motions, respectively, and, for far-field ground motions, b varies in a narrow range. This phenomenon reveals that PGA should not be chosen as the IM for long-period structures under pulse-like ground motions.

6.2. Efficiency

As indicated, the efficiency of IMs can be evaluated with the dispersion obtained from the PSDMs. The larger is, the less efficient IM is. Figures 11 and 12 present the efficiency comparison of the four bridge types under pulse-like ground motions. As illustrated in the figures, PGV tends to be the most efficient IM with the pier height ranging from 4 m to 16 m for the pier curvature ductility. However, when the height is 18 m or 20 m, Iv tends to be the most efficient IM. For the bearing deformation, Ia tends to be the most efficient in the longitudinal direction and PGV is the most efficient IM in the transverse direction. In general, PGV tends to be the most efficient IM. An interesting phenomenon observed is that PGA is not efficient similar to the practicality comparison for the long-period structures, which supports the abovementioned assumption strongly.

Similarly, Table 6 lists the efficiency comparison of the seven structure-independent IMs under far-field ground motions. As shown in the table, Iv tends to be the most efficient IM. For the considered height range, PGA and PGD are less efficient because of relatively large efficiency parameter . Similarly, to demonstrate the pulse effect on the comparison results, Figure 13 shows the efficiency comparison of IMs under pulse-like ground motions and far-field ground motions, where PGA, PGV, Iv, and Ia are chosen as the candidate IMs (conditioning on μy, the other cases follow the similar tendency). As shown in the figure, PGV can be chosen as the appropriate IM to predict the seismic response of bridges under either far-field or pulse-like ground motions with respect to efficiency. In addition, PGA is the least efficient IM under either pulse-like or far-field ground motions due to the largest dispersion .

6.3. Proficiency

The composite measure ζ can be utilized to estimate the proficiency of IMs, which combines practicality and efficiency. A more proficient IM has a lower ζ. The proficiency comparison of IMs is given in Figures 14 and 15. As illustrated in the figures, ζ is smaller for PGV conditioning on the curvature ductility and bearing deformation in the longitudinal and transverse directions, which indicates that PGV is the most proficient IM, followed by CAV. For the considered height range, Iv, CAD, PGD, and PGA tend to be less appropriate because of the relatively high proficiency parameters ζ. Similarly, for long-period cases, PGA is not a proficient IM, which is consistent with the practicality and efficiency comparison results. This phenomenon reveals that PGA is not appropriate to be chosen as the IM to develop PSDMs of long-period structures under pulse-like ground motions.

Similarly, Table 7 summarizes the proficiency comparison of the seven structure-independent IMs under far-field ground motions. As shown in the table, CAV tends to be the most proficient IM, followed by PGV. To demonstrate the pulse effect on the comparison results, Figure 16 shows the proficiency comparison of IMs under pulse-like ground motions and far-field ground motions, where PGA, PGV, CAV, and Ia are chosen as the candidate IMs (conditioning on μy, the other cases follow the similar tendency). As shown in the figure, PGV can be chosen as the appropriate IM to assess the seismic demand under both far-field and pulse-like ground motions with respect to proficiency. In addition, PGA is less appropriate to perform the vulnerability analysis under either pulse-like or far-field ground motions because of the relatively large ζ. One should notice that ζ value of PGA increases extremely with the height of pier under pulse-like ground motions, which reveals that PGA is the least proficient under pulse-like ground motions. This phenomenon is consistent with the practicality and efficiency comparison results.

6.4. Sufficiency

The value, which is obtained with a regression analysis on the residual, ε/IM, from the PSDM relative to the ground motion characteristics, that is, M or R, is an indicator of sufficiency of IMs. Generally, the value of 0.1 is regarded as the cutoff of an insufficient IM.

A larger value signifies that the IM is more sufficient. The sufficiency comparison of the four candidate IMs, namely, PGV, PGA, CAV, and Ia, is given in Figures 17 and 18. As indicated in the figures, PGA tends to be the most sufficient IM conditioning on the pier curvature ductility for the lower height range (from 4 m to 6 m) with respect to M and R in the longitudinal and transverse directions. For the higher height cases (larger than 8 m), CAV tends to more sufficient, followed by PGV. For the displacement-related EDP, namely, bearing deformation, PGA tends to be the most sufficient IM for the lower height cases (4 m to 8 m) with respect to M; with height increasing, PGV is more sufficient. However, CAV tends to be the most sufficient IM with respect to R for the lower height range (4 m to 8 m) in the longitudinal and transverse directions. As height increases, CAV and PGV tend to be the most sufficient IMs in the longitudinal and transverse directions, respectively. In general, both PGA and PGV can be selected as the most sufficient IMs for PSDMs of isolated bridges.

Similarly, Table 8 shows the sufficiency comparison of the seven structure-independent IMs under far-field ground motions. As shown in the table, Iv tends to be the most sufficient IM, followed by PGV and CAV. To demonstrate the pulse influence on the sufficiency comparison results, Figure 19 shows the sufficiency comparison of IMs under pulse-like ground motions and far-field ground motions, where PGA, PGV, CAV, and Ia are chosen as the candidate IMs (conditioning on μy, the other cases follow the similar tendency). As shown in the figure, PGV is found to be more sufficient with respect to R and less sufficient with respect to M under pulse-like ground motions. Under far-field ground motions, PGV is appropriately independent of M and R. Similarly, PGA is relatively dependent on M and R, which implies that PGA is less sufficient.

In general, based on the practicality, efficiency, proficiency, and sufficiency comparison results, PGV can be selected as the optimal IM to assess the seismic demand of isolated bridges under pulse-like and far-field ground motions. However, PGA is not an appropriate IM to establish PSDMs of isolated bridges, which is commonly used in the fragility analysis of structures.

7. Conclusion

The intensity measure (IM) plays a vital role in linking the seismic hazard analysis with the structural demand (e.g., curvature ductility and bearing deformation). Traditional ground motion models used in the probabilistic seismic hazard demand analysis do not account for pulse effect and may underestimate the seismic demand of structures located in the near-fault region. In addition, the research of the optimal intensity measure selection targeted at pulse-like ground motions is rare and there is no widely accepted optimal intensity measure for pulse-like ground motions in the bridge engineering community, which emphasizes the necessity of the research of optimal intensity measures targeted at the pulse-like ground motions.

The objective of this paper is to identify the optimal IM of PSDMs for isolated bridges subjected to pulse-like ground motions. For this purpose, we select four common isolated girder bridge types with varied pier height (from 4 m to 20 m), that is, multicell box girder (MCBG), multibeam box girder (MBBG), concrete I-girder bridge (CIGB), and steel I-girder bridge (SIGB), as case studies. Totally seven structure-independent IMs, namely, peak ground acceleration (PGA), peak ground velocity (PGV), peak ground displacement (PGD), Arias intensity (Ia), velocity intensity (Iv), cumulative absolute velocity (CAV), and cumulative absolute displacement (CAD), are selected as the candidate IMs. The probabilistic seismic demand analysis is performed using 121 pairs of pulse-like ground motions. The engineering demand parameters (EDPs), namely, pier curvature ductility (μx, μy) and bearing deformation (δx, δy), are monitored and reordered. Efficiency, practicality, proficiency, sufficiency, and hazard computability are selected as the selection metrics of an optimal intensity measure.

In general, PGV tends to be the optimal IM for conditioning PSDMs of isolated bridges with the considered pier height range based on the five metrics mentioned above. This phenomenon can be explained in that the response of isolated bridges under pulse-like ground motions is dominated by the pulse characteristic. In addition, as a commonly used IM, PGA is not an appropriate choice to establish the fragility curves of isolated bridges under pulse-like ground motions, which signifies that special attention should be paid to the intensity measure selection when the pulse effect cannot be ignored. The results obtained in this paper are intended to offer guidance for probabilistic seismic demand analysis of isolated bridges near the fault.

Data Availability

All data, models, and codes generated or used during the study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This research was supported by the National Natural Science Foundation of China (51608161 and 51778471), the China Postdoctoral Science Foundation (2016M602007), the Fundamental Research Funds for the Central Universities of China (PA2019GDPK0041), and the Ministry of Science and Technology of China (SLDRCE19-B-19). The supports are gratefully acknowledged.