Proposing an optimal integralbased intensity measure for seismic collapse capacity assessment of structures under pulselike nearfault ground motions
Masood Yakhchalian^{1} , Ahmad Nicknam^{2} , Gholamreza Ghodrati Amiri^{3}
^{1, 2}School of Civil Engineering, Iran University of Science and Technology, P. O. Box 16765163, Tehran, Iran
^{3}Center of Excellence for Fundamental Studies in Structural Engineering, School of Civil Engineering Iran University of Science and Technology, P. O. Box 16765163, Tehran, Iran
^{3}Corresponding author
Journal of Vibroengineering, Vol. 16, Issue 3, 2014, p. 13601375.
Received 22 January 2014; received in revised form 29 March 2014; accepted 7 April 2014; published 15 May 2014
JVE Conferences
Pulselike nearfault ground motions, which are characterized by the presence of a velocity pulse, can impose large demands in structures, and hence can potentially increase seismic collapse risk. One of the most important parameters affecting the structural demands under nearfault ground motions is the pulse period. The commonly used scalar intensity measure (IM), spectral acceleration at the fundamental period of the structure, $Sa\left({T}_{1}\right)$, is demonstrated to be deficient, and also insufficient with respect to the pulse period, for predicting the structural collapse capacity under pulselike ground motions. Furthermore, it is shown that the recently proposed IM, named ${I}_{Np}$, is not able to fully account for the effect of the pulse period. This study proposes an optimal integralbased IM, named ${I}_{Sa}$, for reliable seismic collapse assessment of structures subjected to pulselike nearfault ground motions. ${I}_{Sa}$ is the integral of the pseudo acceleration response spectrum over an optimal period range that is a function of the fundamental period and ductility of the structure. To propose the new IM, 15 generic frame structures with different fundamental periods and ductility classes were employed. The results indicate that ${I}_{Sa}$ is significantly efficient, which causes reduction in the dispersion of structural collapse capacity, and also sufficient with respect to the pulse period and other ground motion characteristics for seismic collapse capacity assessment.
Keywords: intensity measure, seismic collapse assessment, nearfault ground motions, pulse period, efficiency, sufficiency.
1. Introduction
Nowadays, nonlinear dynamic analysis has become a common method for evaluating the demands on structures under earthquake excitations. Thus, it is necessary to investigate the properties of ground motions that are strongly related to the structural response. A parameter that describes the strength of a ground motion and quantifies its effect on structures is called Intensity Measure (IM). In seismic performance assessment of structures, IMs link the ground motion hazard with the structural response, and using an appropriate IM plays an important role in the prediction of structural response or capacity. Therefore, studies continue to classify the existing IMs and also to propose new optimal IMs (e.g., [14]).
Pulselike nearfault ground motions have two known characteristics: forward directivity and fling step. The propagation of rupture toward a site at a velocity close to the shear wave velocity is called ‘forward directivity’. This phenomenon causes a single large longperiod pulse that occurs at the beginning of the record and contains most of the input energy from the rupture process. Commonly, the longperiod pulse is likely to occur in the faultnormal direction at sites in the proximity of an active fault, where the rupture is propagating toward the site. The permanent static displacement in the faultparallel direction for strikeslip faults or in the faultnormal direction for dipslip faults is called ‘fling step’. Unlike the forward directivity that generates a dynamic twosided velocity pulse, the fling step generates a onesided velocity pulse [57].
Past studies that have investigated the effects of nearfault directivity on structures (e.g., [610]) indicate that pulselike ground motions may cause large displacement and strength demands in structures, and hence can increase the risk of structural collapse compared with high amplitude farfield ground motions. Furthermore, due to the importance of the nearfault effects on structures, some researchers have tried to simulate nearfault ground motions consisting of forward directivity or fling step effects (e.g., [11]). As mentioned by Tothong and Cornell [7], the twosided pulse caused by the forward directivity is more damaging than the onesided pulse caused by the fling step. Hence, this study focused on the forward directivity from the faultnormal component of nearfault ground motions. In these ground motions, the period of the velocity pulse (${T}_{p}$) is an important parameter affecting the structural response [6, 10, 12, 13], and its effect is not properly described by common IMs such as spectral acceleration at the fundamental period of the structure, $Sa\left({T}_{1}\right)$.
In PerformanceBased Earthquake Engineering (PBEE) framework [14], the properties of an optimal IM are as follows: efficiency, that is, the ability of an IM to predict the response or capacity of a structure subjected to ground motion records with small dispersion [15]; sufficiency, which is the ability of an IM to render the structural response or capacity conditionally independent of earthquake characteristics such as magnitude ($M$), sourcetosite distance ($R$), and in nearfault ground motions the pulse period (${T}_{p}$); scaling robustness, which represents sufficiency with respect to scale factor ($SF$); and predictability, that is, having a reliable ground motion prediction equation [2].
Since one of the most important parts of the PBEE is assessing the seismic collapse safety of structures, this issue is an interesting field of study for researchers (e.g., [16, 17]). Recently, several studies have focused on improving the reliability of seismic collapse assessment procedures by employing advanced IMs, and also some advanced ground motion selection methods (e.g., [18, 19]). Incremental Dynamic Analysis (IDA) [20] is a powerful analytical method to obtain the seismic collapse capacity of structures. This method involves subjecting a structural model to several ground motion records, each scaled to multiple levels of intensity, measured by an IM, until it causes collapse. When traditional IMs such as peak ground acceleration (PGA) or $Sa\left({T}_{1}\right)$ are used for collapse capacity prediction, the collapse capacity values can display large recordtorecord variability (dispersion), forcing the use of many records to achieve reliable results [15]. However, by employing an optimal IM, more reliable results with lower dispersion and bias can be achieved. In this study, the performance of the commonly used IM, $Sa\left({T}_{1}\right)$, and the IM recently proposed by Bojórquez and Iervolino [3], ${I}_{Np}$, was investigated for collapse capacity prediction of structures under pulselike nearfault ground motions. To improve the shortcomings of $Sa\left({T}_{1}\right)$ and ${I}_{Np}$, an optimal integralbased IM, ${I}_{Sa}$, which increases the efficiency and decreases the potential bias due to the pulse period, was proposed.
2. Considered intensity measures
In this study, $Sa\left({T}_{1}\right)$ and ${I}_{Np}$ were considered as IMs for collapse capacity prediction of structures. ${I}_{Np}$ is defined as follows [3]:
${N}_{p}=\frac{S{a}_{avg}\left({T}_{1}\dots {T}_{N}\right)}{Sa\left({T}_{1}\right)},$
where $S{a}_{avg}({T}_{1}\dots \text{}{T}_{N})$ is the geometric mean of spectral accelerations over the period range of ${T}_{1}\text{}{T}_{N}$ (${T}_{N}={2T}_{1}$). To propose an optimal IM, the integral of the pseudo acceleration response spectrum over an optimal period range, which is a function of the fundamental period and ductility of the structure, was also considered as a scalar IM:
where $PSA(T,5\%)$ is the $5\%$ damped pseudo spectral acceleration at vibration period $T$, ${T}_{1}$ is the fundamental period of the structure, and $\alpha $ and $\beta $ are coefficients to be determined based on the structural characteristics such as fundamental period and ductility. $PSA(T,5\%)$, which is often shortened to spectral acceleration at vibration period $T$, $Sa\left(T\right)$, can be obtained as follows:
where $Sd(T,5\%)$ is the ordinate of the $5\%$ damped displacement response spectrum at vibration period $T$. Displacement response spectrum is a plot of the maximum displacement response of singledegreeoffreedom (SDOF) systems with various natural periods when subjected to an earthquake ground motion [21]. The damping ratio of the SDOF systems is commonly assumed to be equal to $5\%$.
Since the major aim of proposing an IM is applying it in PBEE, a proposed IM should be predictable (should have a reliable ground motion prediction equation) to be able to link the structural response estimates with the seismic hazard. Tothong and Cornell [7] proposed an IM for seismic performance assessment of structures under nearfault excitations, which needs a firstorder Taylorseries expansion to be predictable. But as mentioned by Bradley et al. [2], using Taylorseries expansion causes inaccuracy due to the large uncertainty in ground motion prediction equations. In this study, to propose ${I}_{Sa}$ as an optimal IM the predictability criterion was considered. Therefore, knowing that ${I}_{Sa}$ is the integral of the pseudo acceleration response spectrum over a specific period range, the ground motion prediction equation for this integralbased IM can be obtained from the ground motion prediction equations for spectral acceleration using the method proposed by Bradley [22].
3. Structural models and method of analyses
Because of the need for generality of the results, and hence, of the proposed IM, the structural models of twodimensional generic onebay frames designed by Medina and Krawinkler [23], and also used by Ibarra and Krawinkler [24] were employed for the analyses. The study performed by Medina and Krawinkler [23] showed that onebay generic frames are generally adequate to capture the global behavior of multibay frames. The considered generic frames have different number of stories, $N=$3, 6, 9, 12 and 15, and fundamental periods, ${T}_{1}$, equal to 0.1$N$. The structural models were created with the OpenSees software [25]. For nonlinear time history analyses, 5 % rayleigh damping was assigned to the first mode and the mode at which the cumulative mass participation exceeds 95 %. The P$\mathrm{\Delta}$ effect was considered for modeling, and nonlinear behavior was modeled by using rotational springs at both beam ends and the bottom end of the firststory columns. Bilinear model [26] was used to specify the nonlinear momentrotation behavior of the rotational springs, and cyclic deterioration was neglected. For each member, the postyield and negative postcapping stiffness ratios of 0.03 and –0.1 were assumed, respectively. To account for the effect of ductility in proposing an optimal IM, three member ductility capacity values, $\delta c/\delta y=$2, 4 and 8, were considered as the representatives of low, intermediate and highductile structures, respectively. Thus, a total number of 15 generic frame structures (5 structures for each ductility class) were considered.
A set consisting of 48 pulselike records rotated to the faultnormal direction with a minimum usable frequency less than 0.2 Hz, extracted from the pulselike nearfault ground motion record set used by Tothong and Cornell [27], was selected for the analyses. The selected ground motion records were taken from the PEER NGA database [28]. To obtain the collapse capacities of the structures, IDAs were performed by using the selected ground motion set. The amplitude of each ground motion in the set, measured by an IM, was scaled to an increasing intensity until it causes collapse. Engineering Demand Parameter (EDP) was considered to be the maximum interstory drift ratio, and collapse was assumed to take place when the dynamic instability occurs and the IDA curve becomes flat [29]. Thus, the collapse capacity corresponding to each ground motion record, $I{M}_{col}$, was obtained. In the IDA approach, the Hunt and Fill algorithm [20] was used to obtain the collapse capacities corresponding to all of the ground motion records. The results of IDAs for the highductile 12story structure considering different IMs are presented in Fig. 1(a)(c). The IDA curves were normalized by the median of the collapse capacity values calculated by considering each IM (more details about the IDA curves are discussed in Section 4).
Fig. 1. IDA curves for the 12story structures with $\delta c/\delta y=$8 with obtained by employing: a) $Sa\left({T}_{1}\right)$, b) ${I}_{Np}$ and c) ${I}_{Sa}$ with the integration period range of 0.74${T}_{1}$–4.3${T}_{1}$
a)
b)
c)
The probabilities of collapse at different IM levels, which can be represented with a collapse fragility curve, can be calculated by fitting a lognormal distribution to the $I{M}_{col}$ values [27]:
where $\mathrm{\Phi}\left(\right)$ is the cumulative distribution function (CDF) of the standard Gaussian distribution; ${\mu}_{\mathrm{l}\mathrm{n}{IM}_{col}}$ and ${\sigma}_{\mathrm{l}\mathrm{n}{IM}_{col}}$ are the estimated logarithmic mean and standard deviation of the collapse capacity, $I{M}_{col}$.
4. Efficiency of the IMs for collapse capacity prediction
It should be noted that when using an IM to estimate the $I{M}_{col}$ values by the IDA approach, the observed collapse capacity dispersion, ${\sigma}_{\mathrm{l}\mathrm{n}\text{}{IM}_{col}}$, is closely connected to the IM employed. An IM is more efficient for collapse capacity prediction if it indicates a lower ${\sigma}_{\mathrm{l}\mathrm{n}{IM}_{col}}$ value [15]. The use of an efficient IM for seismic performance assessment of structures increases the reliability of the assessment. In other words, using an efficient IM can considerably reduce the number of analyses required for the estimation of structural response or capacity with a given accuracy. The values of collapse capacity dispersion, ${\sigma}_{\mathrm{l}\mathrm{n}{IM}_{col}}$, obtained by considering $Sa\left({T}_{1}\right)$ and ${I}_{Np}$ as IMs for collapse capacity prediction of the generic frame structures are presented in Table 1. This table shows that using ${I}_{Np}$ for collapse capacity prediction causes reduction in the dispersion of collapse capacity for all of the structures with different ductility classes, when compared with using $Sa\left({T}_{1}\right)$. Fig. 1(a) and (b) show the IDA curves of the highductile 12story structure obtained by employing $Sa\left({T}_{1}\right)$ and ${I}_{Np}$, respectively. It can be seen that using ${I}_{Np}$ as an IM causes reduction in the dispersion of structural collapse capacity when compared with using $Sa\left({T}_{1}\right)$.
Applying the integral of the pseudo acceleration response spectrum over an optimal period range, ${I}_{Sa}$, as an IM can considerably reduce the dispersion of structural collapse capacity. For example, Fig. 1(c) illustrates that considering the period range of 0.74${T}_{1}$4.3${T}_{1}$ for ${I}_{Sa}$ in the case of the 12story structure with $\delta c/\delta y=$8 leads to a more efficient IM than $Sa\left({T}_{1}\right)$ and ${I}_{Np}$ for collapse capacity assessment of this structure. It is obvious that the optimal period range varies by varying the fundamental period and ductility of the considered structure. In Section 6, the optimal period ranges for ${I}_{Sa}$ considering the generic frame structures with different fundamental periods and ductility classes are obtained.
Table 1. Values of collapse capacity dispersion, ${\sigma}_{\mathrm{l}\mathrm{n}\text{}{IM}_{col}}$, obtained for the generic frames
No. of stories

$IM=Sa\left({T}_{1}\right)$

$IM={I}_{Np}$


$\delta c/\delta y$

$\delta c/\delta y$


2

4

8

2

4

8


3

0.273

0.315

0.380

0.223

0.265

0.323

6

0.366

0.370

0.432

0.274

0.275

0.345

9

0.405

0.428

0.479

0.318

0.332

0.385

12

0.375

0.424

0.475

0.285

0.335

0.378

15

0.329

0.356

0.427

0.232

0.260

0.332

The standard error of the collapse capacity associated with a sample of size $n$ (number of ground motion records) can be expressed as:
where ${\sigma}_{\mathrm{l}\mathrm{n}\text{}{IM}_{col}}$ is the logarithmic standard deviation of the collapse capacity. Based on Eq. (5), by applying a more efficient IM, the number of ground motion records, $n$, can be reduced while the standard error remains the same, which means that lower effort and computational expense are needed for performance assessment.
5. Sufficiency of the IMs
A sufficient IM produces the same distribution of demands and capacities independently of the ground motion selection [15]. The sufficiency of an IM for collapse capacity prediction means that the distribution of collapse capacity, obtained by using the IM, is independent of ground motion characteristics such as earthquake magnitude ($M$) and sourcetosite distance ($R$). In pulselike nearfault ground motions, the pulse period (${T}_{p}$) is the most important ground motion characteristic that can considerably affect the structural response. Neglecting the pulse period effect in seismic performance assessment of structures under pulselike ground motions can make the structural response and collapse capacity estimates biased [10, 13]. Because the distribution of collapse capacity is obtained from the results of a finite number of IDAs, sufficiency is one of the important properties of an optimal IM. If this distribution is dependent on the $M$, $R$ and ${T}_{p}$ values of the ground motions used, then the distribution will be biased if the distribution of $M$, $R$ and ${T}_{p}$ of the ground motions used in the IDAs is not the same as that of the ground motions which will occur at the site in the future [2]. In addition, scaling robustness represents sufficiency with respect to scale factor ($SF$). Scaling robustness results in the unbiased prediction of structural response or capacity when the records are linearly scaled to perform structural analyses. As pointed out by Bojórquez and Iervolino [3], sufficiency of an IM is important because a sufficient IM can be used in the probabilistic structural assessment decoupling the hazard and structural analysis.
To examine the sufficiency of an IM with respect to $M$, $R$, $SF$ and ${T}_{p}$ for predicting the collapse capacity of structures, linear regression can be performed between these parameters and the observed collapse capacities from the IDAs by using the following equation:
where ${\mu}_{\mathrm{ln}I{M}_{col}}$ is the expected value of $\mathrm{l}\mathrm{n}I{M}_{col}$; ${b}_{0}$ and ${b}_{1}$ are coefficients to be estimated from linear regression; and $x$ is one of the parameters $M$, $\mathrm{l}\mathrm{n}R$ (natural logarithm of sourcetosite distance), $\mathrm{l}\mathrm{n}SF$ and ${T}_{p}$. Because the linear regression is based on a finite number of observations, it is essential to use statistical tests to determine the significance of the coefficient ${b}_{1}$. Assuming a Studentt distribution for the coefficient ${b}_{1}$, the Ftest can be used to determine the statistical significance of ${b}_{1}$ [30]. In general, a pvalue less than 0.05 from the Ftest indicates that the slope of the linear regression, ${b}_{1}$, is a statistically significant value, which represents the insufficiency of the IM with respect to $x$ for collapse capacity assessment.
5.1. Sufficiency of the IMs with respect to $\mathit{M}$, $\mathit{R}$, and $\mathit{S}\mathit{F}$ for collapse capacity prediction
The sufficiency of $Sa\left({T}_{1}\right)$ and ${I}_{Np}$ with respect to $M$ and $R$ for collapse capacity prediction of the intermediateductile 3story structure ($\delta c/\delta y=$4) is tested in Fig. 2 and 3, respectively. It can be seen that these IMs are insufficient with respect to M because the pvalues from the Ftest are less than 0.05, whereas they are sufficient with respect to $R$ because the corresponding pvalues are greater than 0.05. In other words, these IMs are unable to fully account for the effect of magnitude, and thus, for reliable seismic collapse assessment of this structure by using $Sa\left({T}_{1}\right)$ or ${I}_{Np}$, the magnitude of ground motions should be considered in the ground motion selection procedure. But the sourcetosite distance of ground motions is not necessary to be considered as a determining factor in ground motion selection. As an alternative, the use of a sufficient IM with respect to $M$ can eliminate the need to consider magnitude as an important criterion in the ground motion selection procedure for seismic analyses.
Fig. 2. Testing the sufficiency of $Sa\left({T}_{1}\right)$ and ${I}_{Np}$ with respect to magnitude for collapse capacity prediction of the 3story structure ($\delta c/\delta y=$ 4): a) $Sa\left({T}_{1}\right)$ and b) ${I}_{Np}$
a)
b)
Fig. 3. Testing the sufficiency of $Sa\left({T}_{1}\right)$ and ${I}_{Np}$ with respect to sourcetosite distance for collapse capacity prediction of the 3story structure ($\delta c/\delta y=$4): a) $Sa\left({T}_{1}\right)$ and b) ${I}_{Np}$
a)
b)
The pvalues obtained from testing the sufficiency of $Sa\left({T}_{1}\right)$ and ${I}_{Np}$ with respect to $M$, $R$ and $SF$ for collapse capacity prediction of the structures are presented in Fig. 4(a)(c), respectively. Fig. 4(a) indicates that $Sa\left({T}_{1}\right)$ and ${I}_{Np}$ in some points (structures), which the pvalues are marginally less than 0.05 (indicated by dashed line), are insufficient with respect to magnitude. In Fig. 4(b), it can be seen that $Sa\left({T}_{1}\right)$ and ${I}_{Np}$ are sufficient with respect to sourcetosite distance for all of the structures. It can therefore be inferred that in the case of structures subjected to nearfault ground motions, neglecting earthquake characteristics such as magnitude and sourcetosite distance often does not cause significant bias in the collapse capacity estimates. In fact, $Sa\left({T}_{1}\right)$ and ${I}_{Np}$ can account for the effects of sourcetosite distance and magnitude (marginally for magnitude). Investigating the sufficiency of $Sa\left({T}_{1}\right)$ and ${I}_{Np}$ with respect to scale factor (see Fig. 4(c)) indicates that in the highductile structures unlike ${I}_{Np}$, $Sa\left({T}_{1}\right)$ may slightly cause insufficiency for predicting the structural collapse capacity.
5.2. Sufficiency of the IMs with respect to ${\mathit{T}}_{\mathit{p}}$ for collapse capacity prediction
A sufficient IM for seismic collapse capacity assessment of structures under pulselike nearfault ground motions should be able to account for the effect of ${T}_{p}$. The pulse period, ${T}_{p}$, is defined as the period associated with the global peak of the velocity response spectrum (${S}_{v}$) [7, 31]. As described earlier, the dependence of the structural collapse capacity on the pulse period can be investigated by linear regression (see Eq. (6)). Therefore, the pvalue for ${b}_{1}$, obtained from the Ftest, can be used as an indicator of the IM sufficiency with respect to ${T}_{p}$. The pvalues obtained from testing the sufficiency of $Sa\left({T}_{1}\right)$ and ${I}_{Np}$ with respect to ${T}_{p}$ for collapse capacity prediction of the structures are presented in Fig 4(d). It can be seen that predicting the collapse capacity for most of the structures, except the 3story structures, is dependent on ${T}_{p}$. In other words, using $Sa\left({T}_{1}\right)$ or ${I}_{Np}$ as an IM causes the seismic collapse capacity assessment of the structures to be biased with respect to ${T}_{p}$.
Fig. 4. Pvalues obtained from testing the sufficiency of $Sa\left({T}_{1}\right)$ and ${I}_{Np}$ with respect to different parameters for collapse capacity prediction of the structures: a) $M$, b) $R$, c) $SF$ and d) ${T}_{p}$
a)
b)
c)
d)
Fig. 5(a)(c) show the collapse capacity values obtained using different IMs, normalized by their median, versus ${T}_{p}/{T}_{1}$ values for the intermediateductile 9story structure. In these figures, for illustrating the dependency of the collapse capacity values on ${T}_{p}$, moving average and also in a simple manner linear regression are applied. The shape of the moving average curve, calculated by averaging the point of interest with the four previous and four subsequent data points, implies the impact of the pulse period variations on the structural collapse capacity. It can be seen that the collapse capacity values obtained using $Sa\left({T}_{1}\right)$ and ${I}_{Np}$ are dependent on ${T}_{p}$. As shown in Fig. 5(a) and (b), when the pulse period is approximately equal to the fundamental period of the structure (${T}_{p}/{T}_{1}\approx 1$), the highest collapse capacities are observed, which means that the structure is least susceptible to seismic collapse. Since the resonance phenomenon for linearelastic structures occurs when the dominant period of the excitation is close to the fundamental period of the structure, the coincidence of the period associated with the peak of the moving average curve with the fundamental period of the structure may appear to be irrational. The reason for this coincidence is that when a structure subjected to a ground motion is near collapse, its period elongates considerably due to severely nonlinear responses. Consequently, the collapse capacity of the structure is dependent on the spectral components in the period range of the nearcollapse period far away from ${T}_{1}$, and hence pulselike records with ${T}_{p}/{T}_{1}>1$ cause more severe damages in the structure, which correspond to lower collapse capacity values. This issue has also been mentioned by Champion and Liel [10].
Fig. 5. Normalized collapse capacity values obtained using different IMs versus ${T}_{p}$/${T}_{1}$ values, and the corresponding moving average curves calculated for the 9story structure ($\delta c\u2044\delta y=$ 4): a) $Sa\left({T}_{1}\right)$, b) ${I}_{Np}$ and c) ${I}_{Sa}$ with the integration period range of 0.82${T}_{1}$3.9${T}_{1}$
a)
b)
c)
In Fig. 5, the area (in absolute values) between the moving average curve and the horizontal dashed line, indicating the median for each IM, can be a representative of the bias due to the pulse period existing in the $I{M}_{col}$ values. Therefore, to propose an optimal IM that satisfies the sufficiency with respect to ${T}_{p}$, the proposed IM should reduce the aforementioned area. It can be seen that the area between the moving average curve and the horizontal dashed line obtained for ${I}_{Np}$ is lower than that of obtained for $Sa\left({T}_{1}\right)$. Although employing ${I}_{Np}$ as an IM for seismic collapse capacity assessment of the selected structure under nearfault ground motions reduces the bias compared with employing $Sa\left({T}_{1}\right)$, the results are still biased. The regression lines in Fig. 5(a) and (b) also imply that $Sa\left({T}_{1}\right)$ and ${I}_{Np}$ are insufficient with respect to ${T}_{p}$, because the pvalues are less than 0.05. Consequently, the collapse capacity estimates obtained by employing these IMs are biased, and depend on the ${T}_{p}$ values of the ground motions used for the analyses.
By comparing Fig. 5(c) with Fig. 5(a) and (b), it can be seen that unlike $Sa\left({T}_{1}\right)$ and ${I}_{Np}$, the integralbased IM, ${I}_{Sa}$, with the integration period range of 0.82${T}_{1}$3.9${T}_{1}$ significantly reduces the bias caused by the pulse period. In fact, a pvalue of 0.92 (corresponding to an approximately horizontal regression line), and also significant reduction in the area between the moving average curve and the horizontal dashed line show the high degree of sufficiency of ${I}_{Sa}$ with respect to ${T}_{p}$ for collapse capacity assessment of the structure.
6. Proposing optimal period ranges for ${\mathit{I}}_{\mathit{S}\mathit{a}}$
It appears that selection of an optimal integration period range for ${I}_{Sa}$, as a function of the fundamental period and ductility of the structure, to satisfy the efficiency, and simultaneously sufficiency with respect to $M$, $R$, $SF$ and ${T}_{p}$ is the most important and challengeable part of proposing ${I}_{Sa}$. To determine the optimal period range of integration, it should be considered that the integration limits may vary with the variations of structural characteristics such as fundamental period and ductility. As mentioned earlier, three ductility classes, member ductility values of $\delta c/\delta y=$2, 4 and 8 as the representatives of low, intermediate and highductile structures, were used that resulted in a total number of 15 generic frame structures. For each structure, 7140 different integration period ranges, with the period step of 0.05 sec and the maximum and minimum periods of 0 and 6 sec, were assumed. The dispersion of structural collapse capacity, ${\sigma}_{\mathrm{l}\mathrm{n}{I}_{Sacol}}$, was calculated considering each integration period range. As an example, Fig. 6 shows the variations of ${\sigma}_{\mathrm{l}\mathrm{n}{I}_{Sacol}}$ considering different period ranges for the intermediateductile 9story structure. In this figure, the region corresponding to the period ranges that lead to lower ${\sigma}_{\mathrm{l}\mathrm{n}{I}_{Sacol}}$ values, shown in black, can be easily located. But in the selection of the optimal period ranges, the sufficiency of the IM should also be considered as a determining criterion. Fig. 7(a)(d) show the pvalues obtained from testing the sufficiency of ${I}_{Sa}$, considering different period ranges, with respect to $M$, $R$, $SF$ and ${T}_{p}$ for collapse capacity prediction of the intermediateductile 9story structure, respectively.
Fig. 6. Variations of ${\sigma}_{\mathrm{l}\mathrm{n}{I}_{Sacol}}$ considering different period ranges ($\alpha {T}_{1}\u2013\beta {T}_{1}$) for the intermediateductile 9story structure
For selecting the optimal period ranges, some filters were used to select them based on the dispersion of structural collapse capacity, ${\sigma}_{\mathrm{l}\mathrm{n}{I}_{Sacol}}$, and the pvalues obtained from testing the sufficiency with respect to $M$, $R$, $SF$ and ${T}_{p}$. For each structure, the optimal period ranges were selected to satisfy the following criteria: the dispersion of ${I}_{Sacol}$ values does not exceed 1.2 times the minimum ${\sigma}_{\mathrm{l}\mathrm{n}{I}_{Sacol}}$ value, and simultaneously the period ranges that lead to acceptable pvalues from testing the sufficiency with respect to all of the considered parameters (i.e., $M$, $R$, $SF$ and ${T}_{p}$).
Fig. 7. Pvalues obtained from testing the sufficiency of ${I}_{Sa}$, considering different period ranges ($\alpha {T}_{1}\u2013\beta {T}_{1}$), with respect to different parameters for collapse capacity prediction of the intermediateductile 9story structure: a) $M$, b) $R$, c) $SF$ and d) ${T}_{p}$
a)
b)
c)
d)
After selection of the optimal period ranges for each structure, ($\alpha {T}_{1}\u2013\beta {T}_{1}$) period ranges, by using the aforementioned filters, the mean values of $\alpha $ and $\beta $ were calculated for each structure. Then, thirdorder polynomial regression was applied to the mean values of $\alpha $ and $\beta $ for the structures with the same class of ductility by employing the following equation:
Fig. 8 shows the $\beta $ values selected for the intermediateductile structures ($\delta c/\delta y=$4), and the fitted regression to obtain $\beta $ as a function of ${T}_{1}$. The regression coefficients calculated by fitting thirdorder polynomial regression to the mean $\beta $ values are presented in Table 2 for the different ductility classes.
Fig. 9(a) and (b) indicate the $\alpha $ and $\beta $ values calculated using the fitted regressions for the different ductility classes. It can be seen that $\alpha $ does not vary considerably with the variation of structural ductility. Thus, the mean of $\alpha $ values obtained from regressions, for different ductility classes, was calculated for the structures with the same fundamental period. Then, thirdorder polynomial regression was applied to the mean $\alpha $ values corresponding to different fundamental periods. The regression coefficients calculated for obtaining $\alpha $ as a function of ${T}_{1}$ are presented in Table 2. Fig. 9(b) indicates that the difference between the $\beta $ values for the high and lowductile structures decreases with increase in the height of the structure. The reason for this is that when a structure becomes taller, the contribution of the P$\mathrm{\Delta}$ effect in the structural collapse becomes more significant, and the P$\mathrm{\Delta}$ effect controls the collapse mode instead of the ductility capacity of the structural members. Furthermore, it is worth noting that under pulselike nearfault ground motions the importance of the P$\mathrm{\Delta}$ effect increases significantly [9].
Table 2. Regression coefficients calculated to obtain $\alpha $ and $\beta $ as functions of ${T}_{1}$
$\alpha $ or $\beta $

$\delta c/\delta y$

$a$

$b$

$c$

$d$

$\beta $

2

–3.423

8.328

–5.812

4.218

4

–7.713

20.543

–16.982

8.194


8

–4.717

11.903

–10.486

7.900


$\alpha $

–

1.445

–4.670

4.725

–0.697

Fig. 8. Obtaining $\beta $ as a function of ${T}_{1}$ for the intermediateductile structures
Fig. 9. $\alpha $ and $\beta $ values calculated using the fitted regressions for the different ductility classes: a) $\alpha $ and b) $\beta $
a)
b)
7. Discussion on the proposed IM, ${\mathit{I}}_{\mathit{S}\mathit{a}}$
The purpose of proposing an optimal IM is to achieve more reliable results. In fact, using an optimal IM for seismic collapse assessment of structures increases the efficiency (decreases the dispersion of structural collapse capacity), and also sufficiency with respect to ground motion characteristics, that leads to more realistic seismic collapse assessment of structures. Fig. 10 indicates the values of collapse capacity dispersion, ${\sigma}_{\mathrm{l}\mathrm{n}\text{}{IM}_{col}}$, obtained by using ${I}_{Sa}$ and $Sa\left({T}_{1}\right)$ for the structures. As shown in this figure, employing ${I}_{Sa}$ for seismic collapse capacity prediction of the structures leads to significant reduction in the dispersion of structural collapse capacity compared with using $Sa\left({T}_{1}\right)$. This reduction varies from 26 to 66 percent for different structures.
Fig. 11(a)(d) show the pvalues obtained from testing the sufficiency of the proposed IM with respect to ${T}_{p}$, $M$, $R$ and $SF$ for collapse capacity prediction of the structures, respectively. The pvalues imply that ${I}_{Sa}$ is able to account for the effect of the ground motion characteristics, and also scale factor, for seismic collapse assessment of the structures. As mentioned earlier, one of the most important parameters in nearfault ground motions is the pulse period, which affects the structural response significantly. Thus, it is expected for an IM, proposed for seismic collapse assessment of structures under pulselike nearfault excitations, to be able to account for the effect of the pulse period. By comparing Fig. 11(a) and 4(d), it can be seen that the pvalues obtained from testing the sufficiency of ${I}_{Sa}$ with respect to the pulse period for collapse capacity prediction of the structures are significantly greater than those obtained by employing $Sa\left({T}_{1}\right)$ and ${I}_{Np}$.
Fig. 10. Values of collapse capacity dispersion, ${\sigma}_{\mathrm{l}\mathrm{n}{IM}_{col}}$, obtained by using $Sa\left({T}_{1}\right)$ and ${I}_{Sa}$ for the structures
Fig. 11. Pvalues obtained from testing the sufficiency of ${I}_{Sa}$ with respect to different parameters for collapse capacity prediction of the structures: a) ${T}_{p}$, b) $M$, c) $R$ and d) $SF$
a)
b)
c)
d)
The pvalue from the Ftest, for the regression coefficient ${b}_{1}$ in Eq. (6), is a simple tool to gauge the bias caused by the pulse period in collapse capacity prediction. Another common method to gauge the bias caused by the pulse period is applying moving average to the collapse capacity values. As described in Section 5.2, the area (in absolute values) between the moving average curve and the horizontal line, representing the median value, can be used as an indicator of the bias due to the pulse period (see Fig. 5(a)(c)). For each structure, the areas, between the moving average curve and the corresponding horizontal line, obtained by employing ${I}_{Sa}$ and ${I}_{Np}$ were normalized to that obtained by employing $Sa\left({T}_{1}\right)$. The values of normalized areas are presented in Fig. 12 for all of the structures. It can be seen that employing ${I}_{Sa}$, considering the integration limits proposed in Section 6, significantly reduces the bias caused by the pulse period compared with using $Sa\left({T}_{1}\right)$ and ${I}_{Np}$. The reduction in the area representing bias achieved by employing ${I}_{Sa}$ instead of $Sa\left({T}_{1}\right)$ varies from 33 to 82 percent for the structures.
Based on the results of the analyses, to achieve more efficiency and sufficiency for collapse capacity prediction of an arbitrary structure under pulselike nearfault ground motions by employing ${I}_{Sa}$, there should be an optimal integration period range, which is a function of the fundamental period and ductility of the structure. The proposed integration limits for the three ductility classes can be good estimates for each arbitrary structure. Furthermore, the optimal period range for a structure can also be determined based on the method applied in this paper.
Fig. 12. Values of normalized areas obtained by employing ${I}_{Sa}$ and ${I}_{Np}$
8. Conclusions
To propose an optimal integralbased IM for reliable seismic collapse assessment of structures under pulselike nearfault ground motions, 15 generic frame structures, with different fundamental periods and ductility classes, were employed. The integral of the pseudo acceleration response spectrum over an optimal period range, which is a function of the fundamental period and ductility of the structure, was proposed as a scalar IM, ${I}_{Sa}$. Efficiency and sufficiency for collapse capacity prediction were considered as two determining criteria to obtain the optimal period range. It was shown that employing ${I}_{Sa}$, with the proposed upper and lower integration limits, results in more efficiency compared with using $Sa\left({T}_{1}\right)$ and ${I}_{Np}$ (2666 percent reduction in the dispersion of structural collapse capacity when compared with $Sa\left({T}_{1}\right)$). The results also indicated that using $Sa\left({T}_{1}\right)$ or ${I}_{Np}$ as an IM for predicting the structural collapse capacity causes significant bias with respect to the pulse period for all of the structures, except the 3story structures. But using ${I}_{Sa}$, considering the proposed integration limits, causes unbiased structural collapse capacity assessments with respect to the pulse period. In other words, ${I}_{Sa}$ can cause 3382 percent reduction in the area representing the bias due the pulse period compared with $Sa\left({T}_{1}\right)$. Furthermore, testing the sufficiency of ${I}_{Sa}$ for collapse capacity assessment revealed that this IM is also sufficient with respect to magnitude, sourcetosite distance and scale factor.
References
 Baker J. W. Vectorvalued ground motion intensity measures for probabilistic seismic demand analysis. Ph.D. Dissertation, Department of Civil and Environmental Engineering, Stanford University, CA, 2005. [Search CrossRef]
 Bradley B. A., Dhakal R. P., MacRae G. A., Cubrinovski M. Prediction of spatially distributed seismic demands in specific structures: ground motion and structural response. Earthquake Engineering and Structural Dynamics, Vol. 39, Issue 5, 2010, p. 501520. [Search CrossRef]
 Bojórquez E., Iervolino I. Spectral shape proxies and nonlinear structural response. Soil Dynamics and Earthquake Engineering, Vol. 31, Issue 7, 2011, p. 9961008. [Search CrossRef]
 Bojórquez E., Iervolino I., ReyesSalazar A., Ruiz S. E. Comparing vectorvalued intensity measures for fragility analysis of steel frames in the case of narrowband ground motions. Engineering Structures, Vol. 45, 2012, p. 472480. [Search CrossRef]
 Somerville P., Smith N., Graves R., Abrahamson N. Modification of empirical strong motion attenuation relations to include the amplitude and duration effects of rupture directivity. Seismological Research Letters, Vol. 68, Issue 1, 1997, p. 199222. [Search CrossRef]
 Kalkan E., Kunnath S. K. Effects of fling step and forward directivity on seismic response of buildings. Earthquake Spectra, Vol. 22, Issue 2, 2006, p. 367390. [Search CrossRef]
 Tothong P., Cornell C. A. Structural performance assessment under nearsource pulselike ground motions using advanced ground motion intensity measures. Earthquake Engineering and Structural Dynamics, Vol. 37, Issue 7, 2008, p. 10131037. [Search CrossRef]
 Luco N., Cornell C. A. Structurespecific scalar intensity measures for nearsource and ordinary earthquake ground motions. Earthquake Spectra, Vol. 23, Issue 2, 2007, p. 357392. [Search CrossRef]
 Sehhati R., RodriguezMarek A., Elgawady M., Cofer W. F. Effects of nearfault ground motions and equivalent pulses on multistory structures. Engineering Structures, Vol. 33, Issue 3, 2011, p. 767779. [Search CrossRef]
 Champion C., Liel A. The effect of nearfault directivity on building seismic collapse risk. Earthquake Engineering and Structural Dynamics, Vol. 41, Issue 10, 2012, p. 13911409. [Search CrossRef]
 Nicknam A., Barkhordari M., Hamidi Jamnani H., Hosseini A. Compatible seismogram simulation at near source site using MultiTaper Spectral Analysis approach (MTSA). Journal of Vibroengineering, Vol. 15, Issue 2, 2013. [Search CrossRef]
 Alavi B., Krawinkler H. Behavior of momentresisting frame structures subjected to nearfault ground motions. Earthquake Engineering and Structural Dynamics, Vol. 33, Issue 6, 2004, p. 687706. [Search CrossRef]
 Baker J. W., Cornell C. A. Vectorvalued intensity measures for pulselike nearfault ground motions. Engineering Structures, Vol. 30, Issue 4, 2008, p. 10481057. [Search CrossRef]
 Deierlein G. G. Overview of a comprehensive framework for performance earthquake assessment. PEER Report 2004/05, Pacific Earthquake Engineering Research Center, University of California, Berkeley, CA, 2004, p. 1526. [Search CrossRef]
 Vamvatsikos D., Cornell C. A. Developing efficient scalar and vector intensity measures for IDA capacity estimation by incorporating elastic spectral shape information. Earthquake Engineering and Structural Dynamics, Vol. 34, Issue 13, 2005, p. 15731600. [Search CrossRef]
 Haselton C. B., Liel A. B., Deierlein G. G., Dean B. S., Chou J. H. Seismic collapse safety of reinforced concrete buildings: I. assessment of ductile moment frames. Journal of Structural Engineering, Vol. 137, Issue 4, 2011, p. 481491. [Search CrossRef]
 Liel A. B., Haselton C. B., Deierlein G. G. Seismic collapse safety of reinforced concrete buildings: II. comparative assessment of nonductile and ductile moment frames. Journal of Structural Engineering, Vol. 137, Issue 4, 2011, p. 492502. [Search CrossRef]
 Haselton C. B., Baker J. W., Liel A. B., Deierlein G. G. Accounting for ground motion spectral shape characteristics in structural collapse assessment through an adjustment for epsilon. Journal of Structural Engineering, Vol. 137, Issue 3, 2011, p. 332344. [Search CrossRef]
 Yakhchalian M., Ghodrati Amiri G., Nicknam A. A new proxy for ground motion selection in seismic collapse assessment of tall buildings. The Structural Design of Tall and Special Buildings, 2013. [Search CrossRef]
 Vamvatsikos D., Cornell C. A. Incremental dynamic analysis. Earthquake Engineering and Structural Dynamics, Vol. 31, Issue 3, 2002, p. 491514. [Search CrossRef]
 Chopra A. K. Dynamics of structures: theory and applications to earthquake engineering. Fourth Edition, Prentice Hall, Upper Saddle River, New Jersey, 2012. [Search CrossRef]
 Bradley B. A. Sitespecific and spatially distributed groundmotion prediction of acceleration spectrum intensity. Bulletin of the Seismological Society of America, Vol. 100, Issue 2, 2010, p. 277285. [Search CrossRef]
 Medina R. A., Krawinkler H. Seismic demands for nondeteriorating frame structures and their dependence on ground motions. PEER Report 2003/15, Pacific Earthquake Engineering Research Center, University of California, Berkeley, CA, 2003. [Search CrossRef]
 Ibarra L. F., Krawinkler H. Global collapse of frame structures under seismic excitations. PEER Report 2005/05, Pacific Earthquake Engineering Research Center, University of California, Berkeley, CA, 2005. [Search CrossRef]
 Open System for Earthquake Engineering Simulation (OpenSees). Pacific Earthquake Engineering Research Center, University of California, Berkeley, CA, 2012, http://opensees.berkeley.edu. [Search CrossRef]
 Ibarra L. F., Medina R. A., Krawinkler H. Hysteretic models that incorporate strength and stiffness deterioration. Earthquake Engineering and Structural Dynamics, Vol. 34, Issue 12, 2005, p. 14891511. [Search CrossRef]
 Tothong P., Cornell C. A. Probabilistic seismic demand analysis using advanced ground motion intensity measures, attenuation relationships, and nearfault effects. PEER Report 2006/11, Pacific Earthquake Engineering Research Center, University of California, Berkeley, CA, 2007. [Search CrossRef]
 Pacific Earthquake Engineering Research Center. PEER Next Generation Attenuation (NGA) Database, 2013, http://peer.berkeley.edu/nga. [Search CrossRef]
 Villaverde R. Methods to assess the seismic collapse capacity of building structures: state of the art. Journal of Structural Engineering, Vol. 133, Issue 1, 2007, p. 5766. [Search CrossRef]
 Ang A. H. S., Tang W. H. Probability concepts in engineering planning and design. Volume I – basic principles. Wiley, Inc., New York, 1975. [Search CrossRef]
 Alavi B., Krawinkler H. Effects of nearfield ground motion on frame structures. Report No. 138, John A. Blume Earthquake Engineering Center, Department of Civil and Environmental Engineering, Stanford University, Stanford, CA, 2001. [Search CrossRef]