Research Article  Open Access
Sahar Mohammadzadeh Osalu, Hamzeh Shakib, "The Effect of Foundation Flexibility on Probabilistic Seismic Performance of PlanAsymmetric Buildings with Different Strength Distributions", Advances in Civil Engineering, vol. 2020, Article ID 5191508, 23 pages, 2020. https://doi.org/10.1155/2020/5191508
The Effect of Foundation Flexibility on Probabilistic Seismic Performance of PlanAsymmetric Buildings with Different Strength Distributions
Abstract
In this research, the probabilistic seismic performance of asymmetric reinforced concrete wallframe buildings with different strength distributions incorporating foundation flexibility effects is examined. By using probabilitybased performance evaluation approach, it is possible to provide a more accurate prediction of the different strength distribution effect on the seismic performance of asymmetric buildings and find the most efficient strength distribution for meeting each performance level. These efficient distributions can be adopted in the performancebased design of asymmetric buildings. For this purpose, first, the regression analysis and the concepts of efficiency and sufficiency were used to determine an optimal intensity measure (IM) for incremental dynamic analysis and evaluating the seismic response of the considered building models. Then, the proper magnitude of interstory drift capacity for this type of buildings in each limit state was estimated using the damage index concept. Finally, the effects of different strength distributions and the flexibility of foundation were studied on the seismic performance of the asymmetric buildings by investigating the mean annual frequencies of exceeding structural performance levels and confidence levels to satisfy performance objectives. It is concluded that irregular distributions of stiffness and strength in the plan of a building highly affect the seismic performance of buildings. Also, the results show that the optimum strength distribution is a function of the objective performance level and these optimum strength distributions are the same for both fixed and flexiblebase conditions. Meanwhile, the flexible effect of foundation increases the mean annual frequencies of exceedance within the range of 10% to 45% and significantly decreases the confidence levels in most cases.
1. Introduction
The experience of past earthquakes has demonstrated that asymmetric buildings are more vulnerable to intensive earthquakes compared to symmetric buildings. Irregular distributions of stiffness and strength among the lateral forceresisting elements (LFREs) of a building can result in coupled lateraltorsional dynamic response which is the main cause of extensive damages in asymmetric buildings [1, 2]. So far, there have been two different concepts for distribution of strength and stiffness among the LFREs of the buildings. Earlier studies traditionally assumed that the stiffness of an element is an independent parameter and can be estimated based on its geometric characteristics. However, in recent years, it has been pointed out that for many commonly used LFREs such as concrete momentresisting frames and shear walls, stiffness is a strengthdependent parameter and will be calculated after the strength allocation among the LFREs [3–5]. These two different types of LFREs are called Ktype [6] and Dtype [7] elements, respectively.
In Dtype elements, the yield displacement is independent of stiffness and strength and can be estimated based on the geometry and material properties. Therefore, it can be said that stiffness and strength of a Dtype element depend on its yield displacement, and both stiffness and strength eccentricities are important parameters that affect the seismic response of structures [5]. Paulay [8] was the first to study the seismic behavior of asymmetric buildings composed of elements having strengthdependent stiffness, and he claimed that the common lateral strength allocation specified in modern seismic building codes is not rational. Tso and Myslimaj [9, 10] proposed a yield displacement distributionbased approach for strength allocation to Dtype LFREs which does not need the information of stiffness distribution prior to the strength allocation. They used a singlestory walltype model and investigated the effects of different arrangements of centers of stiffness (CS) and strength (CR) with respect to the center of mass (CM) on the diaphragm rotation response of the model. They concluded that the optimal arrangement of the centers is one where the centers of strength and stiffness are located on opposite sides of the mass center. This configuration is known as “balanced CSCR location” in the literature. The authors also demonstrated that there is a noticeable change in the trend of the diaphragm rotation at different seismic input levels. At lower earthquake levels, the minimum rotations are obtained in the models with minimum stiffness eccentricity, but there is a shift towards minimum strength eccentricity as the seismic input level increases [9]. Shakib and Ghasemi [11] studied the effects of nearfault and farfault ground motions on structural responses of asymmetric buildings supported by wall elements. Their results indicated that considering Dtype elements, in nearfault ground motions, the balanced configuration of the stiffness and strength centers leads to the minimum rotational response. Aziminejad and Moghadam [12, 13] studied the performance of walltype asymmetric buildings by using nonlinear dynamic analyses. The results were presented in the form of different limitstate fragility curves. They concluded that for interstory drift responses, models with balanced configuration and strength eccentricity equal to 0.25 of the distance between CR and CS perform better. However, in general, the proper configuration of centers differs significantly depending on the selected engineering demand parameter.
Most of the above studies were carried out by employing simple models with shear wall systems, while in most real systems the frames are also involved in lateralload resistance. For this reason, in this study, the seismic performance of planasymmetric reinforced concrete buildings with shear wallframe system has been evaluated. In addition, since there is no general agreement on the optimal configuration of the various centers, an attempt was made to quantify the effects of different strength distributions on the seismic performance of structures and identify a more reliable configuration of centers in each limit state by using a probabilistic performance assessment approach.
On the other hand, the global seismic behavior of a specific structure may be totally different when it is founded on flexible soil compared to a fixed base. This difference is mainly due to the soilstructure interaction (SSI) phenomenon. Several researchers have so far attempted to incorporate the effects of the foundation flexibility on asymmetric buildings [14–16]. Shakib and Fuladgar [17] studied the seismic response of asymmetric structures composed of Ktype elements with soil flexibility. They concluded that SSI effects reduce the lateral and torsional responses of asymmetric buildings. However, as the structural period is increased, the effect of base flexibility is noticeably decreased. They also found that the eccentricity ratio of the asymmetric system has a significant effect on the response of soilstructure interaction. Roy and Chandra Dutta [18] studied the inelastic seismic demand of lowrise asymmetric buildings with soil flexibility and assuming that the stiffness is independent of the strength in the LFREs of the structural models. They concluded that the nonlinear behavior of the asymmetric system compared to the reference symmetric model is not extremely affected by SSI. However, the effect of soilstructure interaction on the nonlinear range of structural behavior may result in increased damages. Shakib and Atefatdoost [19] performed nonlinear analysis of threedimensional dynamic soilstructure interaction of asymmetric buildings with Dtype walls. Their results showed that with increasing base flexibility, the slab rotation of asymmetric buildings decreases, and the minimum slab rotation occurs when CS and CR are on the opposite sides of CM (balanced condition). More recently, Atefatdoost et al. [20] studied the effect of SSI on the seismic response of idealized singlestory asymmetric buildings subject to nearfault pulselike earthquakes. They concluded that for flexiblebase structures, the reduction in the base shear results only in the center of stiffness to control torsional response; meanwhile, for structures with fixed base, both centers of stiffness and strength are controlling.
In the literature, there is no study investigating base flexibility effects on asymmetric buildings from the standpoint of probabilistic performancebased earthquake engineering. Therefore, this study has provided an opportunity to examine the influence of soilfoundation flexibility on the seismic behavior, in addition to strength distribution effects by implementing incremental dynamic analysis (IDA) and reliabilitybased probabilistic approach. For this purpose, first, the efficiency and sufficiency of different intensity measures (IMs) for evaluating the seismic response of the models were investigated and an optimal IM was chosen for performance assessment. Then, the proper magnitude of interstory drift capacity for RC wallframe buildings in each limit state was estimated using the damage index concept. Finally, the effect of different strength distributions was studied on the probabilistic performance of asymmetric buildings in fixed and flexiblebase conditions by evaluating the mean annual frequencies of exceedance and confidence levels to satisfy performance objectives.
2. Characteristics of WallFrame Buildings with DType Elements
For studying the seismic behavior of wallframe buildings with Dtype elements, first, the yield displacements of the elements should be calculated. As previously mentioned, research concerning the concrete momentresisting frames and shear walls has pointed out that the yield displacement of these LFREs depends only on their geometric characteristic and the properties of the material used and is essentially independent of their strength [3, 5]. The yield displacement of a cantilever shear wall, , as recommended by Paulay [3], can be obtained from the following equation:in which is the steel reinforcement yield strain, is the wall thickness, is its height, and is a coefficient which quantifies the lateral force distribution pattern. For example, for shear walls with lateral force applying at the top, is equal to 0.66 [3].
Also, for a concrete momentresisting frame, the following equation can be used to estimate the yield displacement at roof level [5]:where and are the length and depth of the beam, respectively.
According to Reference [3], within reasonable limits, the strength of the elements may be assigned arbitrarily, and so the designer has also the freedom to choose the fractions of base shear force carried by the wall and the frame in wallframe buildings. This arbitrariness can make it possible for the designer to choose the more rational and economical design solution from a number of possible solutions [21]. After calculating the yield displacement of each element, , and assigning its strength, , its stiffness () can be estimated by the following equation:
The global yield displacement of the wallframe system, , can also be obtained fromwhere and are the wall and frame strengths, respectively.
3. Numerical Modeling
3.1. Structural System
A singlestory RC wallframe system is considered in this study as a reference symmetric model, Model 1 (Figure 1). This model consists of a rigid deck of plan dimensions 12 m × 8 m that is supported by two reinforced concrete walls and two momentresisting frames in ydirection, and two walls and four momentresisting frames in xdirection, as shown in the figure. The story height is considered equal to 3 m. The model is assumed to be located on a site with very high seismicity and a sitespecific earthquake acceleration of 0.35 g (Tehran, Iran) according to the Iranian Seismic Design Code [22]. The design gravity load of the model is 150 tonf, and the design base shear, which has been calculated based on Reference [22], is equal to 24 tonf. According to the definition of dual wallframe buildings in this code, the total base shear has been distributed between walls and moment frames with the ratio 3 to 1 in each direction. All walls have a length of 2 m and a thickness of 0.25 m with two layers of distributed reinforcement, φ16@150 mm. Also, the same cross section of 0.3 × 0.3 m^{2} is considered for all beams and columns with 1.5% and 2.5% reinforcement ratios, respectively. A concrete compressive strength of 21 MPa and a yield strength of 300 MPa are assumed for concrete and steel reinforcement, respectively.
The total lateral strengths of all asymmetric models considered in the study (Models 2 to 10) along both directions are assumed to be equal to that of the symmetric model. In order to prepare the asymmetric models, the dimensions of the cross sections and the yield strain of reinforcements in the left and right side elements in the ydirection are changed so that the normalized distance between the yield displacement center (CD) and the mass center (CM) of the models becomes equal to 15% (). To create strength eccentricity, a yield displacement distributionbased approach has been used as suggested by Tso and Myslimaj [9]. Based on this approach, the strength assignment is accomplished by the following two conditions: (a) the strength eccentricity, , is considered linearly proportional to the yield displacement eccentricity, e_{D}, by the relation “,” where is a variable parameter that is selected by the designer; and (b) the radius of gyration for strength distribution is equal to that for the yield displacement distribution. After distributing the strength among the elements, the stiffness of each element can be calculated by equation (3). The general properties of the structural models are presented in Table 1. In this table, is the stiffness eccentricity, also expressed as a percentage of the plan length dimension, like and . , , and represent the first three modal periods (two sway modes and one torsional mode), and and define the lateral sway (corresponding to ) and torsional (corresponding to ) frequencies of the structural models, respectively. The torsional to the lateral frequency ratio () indicates that all considered models are torsionally stiff (). To better clarify the dynamic behavior of the models, the modal mass participation factors for the sway mode in ydirection ( and ) and the torsional mode ( and ) are also shown in Table 1. It is clear that by reducing the parameter from 1.0 to −1.0, the correlation between the lateral and torsional components of displacements increases, consequently, the ratios of and increase, and the corresponding values for and decrease. The mass center of all models has been considered at the geometric center of the rigid diaphragm.

Dynamic nonlinear analyses of the models were performed by OpenSees software [23]. The damping ratio of the models was assumed 5%, and Rayleigh mass and stiffness proportional damping was considered for damping properties [24]. The peakoriented hysteretic model, proposed by Ibarra et al. [25], was used for nonlinear modeling of the elements. Based on this model, nonlinear behavior of beams and columns is modeled by a linear middle part and two lumped plastic hinges at the end of elements, and nonlinear model of walls is composed of a linear part and one plastic hinge at the base [26]. It should be noted that this model considers all degradation modes including the basic and postcapping strength and the loading and accelerated unloading stiffness degradation effects in the cyclic response. A typical backbone curve of the model is shown in Figure 2. In this figure, and represent the yielding and capping moment capacities which can be calculated based on the lateral strength assigned to the element. , , and are the yielding, plastic, and postcapping rotation capacities for the lumped plastic hinges, respectively. In this study, the values of rotation capacities for shear walls and frame elements were adopted from Krawinkler and Zareian [26]. Based on this reference, can be calculated using the geometric and mechanical characteristic of the element. For momentresisting frames, the parameters and are taken equal to 0.03 and 0.15, respectively, and for shear walls, the corresponding parameters are assumed to be 0.02.
3.2. Foundation System
In this study, according to the Iranian Seismic Design Code site classification [22], the underlying soil is assumed to be of Type III () which is equivalent to Soil Class D of ASCE 4113 [27]. The properties of sandy soil material considered for this study are listed in Table 2. The Iranian Foundation Design Code [28] was used for the design of strip foundations. Based on this code, the safety factors of 3.0 and 1.5 were considered against the bearing and sliding capacity failure of the foundations, respectively. Figure 3 shows the arrangement and the dimensions of strip footing elements considered in this study. As it is seen in the figure, the adjacent footing elements were connected to each other to enforce harmonized motions during the seismic excitations.

The wellknown theory of beam on nonlinear Winkler foundation (BNWF) was utilized for modeling the flexibility of the foundation, as used in several studies [29–38]. The simplicity and capability of the BNWF model in accounting for nonlinear behavior of shallow foundations are its two main beneficial features. This model includes elastic beam elements to model the footing of the structure and linear/nonlinear zerolength springs to capture the soil material behavior [30]. The arrangement of these zerolength springs is schematically shown in Figure 4(a). The horizontal springs are known as tx and px springs which capture sliding and passive resistance of the foundation in horizontal movements, respectively. The vertical springs, called qz springs, are used to simulate footing uplift and compression yielding in the soil medium. The nonlinear backbone curves of these three different sets of springs are shown in Figures 4(b)–4(d). The expressions describing the backbone curves of the springs are similar and have been originally mentioned in Reference [38].
(a)
(b)
(c)
(d)
As shown in Figures 4(b)–4(d), each backbone curve has two parts. In the elastic part of the backbone curve, a linear loaddisplacement equation is considered, in which shows the initial elastic stiffness:
According to Reference [27], the initial elastic stiffness of soil springs depends on the rigidity and bending stiffness of the foundation. Strip foundations can be considered rigid when
In this equation, and are the width and length of the footing, respectively, is the elastic modulus of the footing material (concrete), is Poisson's ratio of soil, is the second moment of area of the footing section, and is the shear modulus of the soil material. Note that is expected to decline at large strains due to increasing seismic input levels; this effect can be considered in the modeling by using effective shear modulus ratios. Table 3 summarizes the values of shear modulus ratios recommended by ASCE 4113 [27] for soil type D.

After determining the rigidity of the foundation, the vertical spring coefficient of the foundation can be calculated using the following equation [27]:
Generally, the horizontal spring coefficients for both rigid and flexible foundations are expressed as
In the nonlinear part, the backbone curve is defined by
In this equation, and are the ultimate and the yielding capacities, respectively, is the displacement where 50% of is mobilized, is the displacement corresponding to the yielding pressure, and c and n are constants that control the shape of the nonlinear part of the curve. Raychowdhury and Hutchinson [31] calibrated these two constants from shallow footing tests, and some relations were suggested by Harden et al. [32] for determination of for each set of springs, which were also used in the current study.
For modeling of the soil material properties, the QzSimple1 (Figure 4(b)), PxSimple1 (Figure 4(c)), and TxSimple1 (Figure 4(d)) material models were employed from OpenSees, and the spring elements were distributed along the footings’ length according to the recommendation by Harden et al. [32].
3.3. Ground Motions
A total number of 20 twodirectional ground motions described in Table 4 have been used in this study, selected from the Ground Motion Database of the Pacific Earthquake Engineering Research Center [39]. All ground motion records were recorded on soil type D () with large magnitudes (), and the closest distance to the fault rupture () was considered larger than 10 km with no marks of directivity effects.

4. Incremental Dynamic Analysis
The effects of different strength distributions and soilstructure interaction on the performance of the structural models were examined by using the incremental dynamic analysis (IDA) method. This method is considered as one of the most efficient parametric analysis methods that can accurately predict the complete range of the structural response [40]. In the IDA method, the ground motion intensity measure (IM) is employed as a connector between seismicity properties and structural seismic analysis and the seismic response of the structure is characterized by the engineering demand parameter (EDP), as an observable indicator. The most common EDP for structural buildings is the maximum interstory drift ratio (). The interstory drift capacities at different performance levels for typical lateralload resisting systems of buildings have been defined in different seismic codes [27, 41–43]. However, there is no reference provided for interstory drift capacities for wallframe buildings. In this section, a suitable IM is determined for the considered structural models by using regression analysis and the concepts of efficiency and sufficiency. Then, an attempt was made to estimate interstory drift capacities of wallframe buildings in different limit states.
4.1. The Optimal Intensity Measure
As mentioned at the beginning of this section, an IM serves as an interface parameter between seismic hazard and structural seismic analysis and this implies that to have a more accurate prediction of the structural demand, the correlation between EDPs and the selected IM is highly important in performancebased earthquake engineering. According to Reference [44], the required characteristics for selecting an appropriate IM are efficiency and sufficiency. An efficient IM leads to a relatively small variability of EDP given IM and consequently a small number of ground motion records and dynamic analyses to evaluate the probability of exceedance of an EDP given IM [44]. On the other hand, a sufficient IM is defined as one that results in an EDP that is conditionally independent of earthquake magnitude () and sourcetosite distance () [44]. Using a sufficient IM is favorable because it neutralizes the effects of magnitude and source distance on the estimation of EDP [44]. Many studies are reported in the literature on the specification of efficient and sufficient IMs in estimating seismic demands of structures. The most commonly used IMs for structural buildings are peak ground acceleration (), peak ground velocity (), and spectral acceleration at the fundamental period of structure (). Elenas and Meskouris [45] studied the correlation between alternative IM parameters and damage indices of reinforced concrete frames. They concluded that correlates well to the displacement response EDPs compared to . Akkar and Küçükdoǧan [46] also confirmed that can be considered as a reasonable IM for predicting the peak nonlinear displacement response of structural systems. On the other hand, a number of previous studies have shown that is a more efficient IM than for drift response EDPs [47, 48]. However, for some structures, may produce undesirable results from the efficiency or the sufficiency point of view [44, 49–54]. Cordova et al. [49] proposed a twoparameter seismic intensity measure, , which reflects both spectral intensity and shape. Their results indicated that the proposed IM significantly decreases the number of ground motion records necessary in time history analyses. Vamvatsikos and Cornell [50] studied three different types of building with different fundamental periods and found that when the influence of higher modes is negligible, taking the spectral acceleration at an optimum period higher than the firstmode period provides good efficiency. Bianchini et al. [55] and Pirizadeh and Shakib [56] demonstrated that equal to which is defined as the elastic 5%damped geometric mean of the spectral accelerations over the range of periods from αT_{1}to ηT_{1} with α = 0.5 and η = 1.5 is an appropriate IM for multistory steel momentresisting frame buildings. It was concluded that selecting as the IM can significantly reduce the variability of IDA curves with respect to the elastic 5%damped firstmode spectral acceleration () and . Eads et al. [57] studied the performance of spectral acceleration averaged over a period range () as an IM for estimating the collapse risk of the RC multistory momentresisting frame and also walltype structures. They concluded that , calculated using an appropriate period range, is generally more efficient and sufficient and provides more stable collapse risk estimates than . They also found that a period range from 0.2T_{1} to 3T_{1} is more likely to be an optimal range for considered structures. Tsantaki et al. [52] compared two alternative IMs, and considering Pdelta effects, with the widely used IM, and without Pdelta, and found that both proposed IMs are more efficient than the benchmark IM. In this study, was suggested as an optimal averaged IM for SDOF structures with fundamental periods larger than 0.15 sec. Kohrangi et al. [53, 54] investigated the efficiency and sufficiency of scalar and vector ground motion IMs including and for the seismic response assessment of RC frame buildings. They observed that is more likely to be the most efficient IM in the estimation of structural responses. However, employing vector IMs leads to an improvement in response prediction efficiency. AlBaadani and AlShameri [58] investigated several IMs for evaluation of the seismic response of RC frame minarets and concluded that the optimal IM for these structures is which defines as the geometric mean of the spectral accelerations at and 2T_{1}. Colajanni et al. [59] carried out a correlation study between interstory drift ratios and three different intensity measure parameters, including , , and in irregular RC structures. The results of their study confirmed that is highly correlated to drift responses, and as a result, it is more efficient than and .
Based on the IMs used in the abovementioned studies, 30 different IMs as categorized in Table 5 were considered as potential IMs for predicting the seismic response of the considered structural models in this study. The optimal one was identified based on both efficiency and sufficiency characteristics.

4.1.1. Efficiency
According to Reference [60], a natural logarithmic correlation can be chosen between EDP and IM to utilize a linear regression of the natural logarithms of IM and EDP:
Thus, efficiency can be measured based on the conditional standard deviation of the regression, , which depends on two regression coefficients () and can be calculated by equation (11). The lower values of belong to a more efficient IM [58, 61, 62]:where the total number of dynamic analyses has been denoted by .
The results of regression analysis for all investigated IMs, in terms of dispersion values (σ), are presented in Table 6 for the symmetric model of reference and six asymmetric models with different values as defined in Table 1. It can be seen that in all models, the maximum value for among all IMs belongs to (IM no. 3), which shows that this IM is the least effective one. After , (no. 1), (no. 7), and (no. 30) have less efficiency in most cases. Furthermore, it can be observed that among peak ground responses, (IM no. 2) with lower dispersion values represents a good IM from the efficiency point of view compared to and . The efficiency of this velocitybased IM is even higher than that of some spectral accelerationbased IMs, including , , , , and . Among all IMs, (no. 15) has the smallest dispersion value for the asymmetric models with β = −1.0 and −0.5, which have large stiffness eccentricities. In the other models, (no. 21) with minimum is the most efficient IM.

Figure 5 demonstrates a comparison between the least and the most efficient IMs for the symmetric model ( and ). It is observed that there is a higher dispersion in the diagram and than that in and . The corresponding results are shown in Figure 6 for the asymmetric model with β = −1.0, where it is seen that the scattering about the regression trend line for , as the most efficient IM, is smaller than that for .
(a)
(b)
(a)
(b)
4.1.2. Sufficiency
Using a oneparameter regression analysis of the residuals on the earthquake record characteristics ( and ) simplifies the procedure for quantifying the sufficiency of an IM [44]. The residuals indicate the difference between the EDP values obtained from dynamic analysis and the estimated ones from the regression fit. According to Reference [44], the statistical significance of the regression estimate on and and consequently the sufficiency of IM can be measured by the value, so that a small value for a given IM (less than about 0.05) indicates that the considered IM is insufficient.
The sufficiency results of all considered IMs are presented in Table 7, where the values higher than 0.05 for both and parameters are highlighted. As shown in the table, the sufficiency of (IM no. 15) and (no. 21), which were the most efficient IMs, is not acceptable. Among peak ground responses, only in the symmetric model provides good sufficiency with respect to both and parameters. It can also be seen that (no. 18), (no. 23), (no. 24), (no. 29), and (no. 30) are sufficient IMs in terms of and for all models. On the other hand, according to Table 6, has the least among these five sufficient IMs and hence is the most efficient of the group. Therefore, is selected as the optimal IM for the structural models used in this study. For example, the sufficiency results of this IM are presented in Figures 7 and 8 for the symmetric model and the asymmetric model with β = −1.0, respectively. For improved readability, is renamed in the IDA curves presented in the next sections.
(a)
(b)
(a)
(b)
4.2. Interstory Drift Capacities
In this paper, following the ASCE 4113 [27], three common performance levels, including immediate occupancy (IO), life safety (LS), and collapse prevention (CP) levels are used for the probabilistic performancebased analysis. According to this code, the interstory drift capacities for these performance levels are determined to be 0.5%, 1%, and 2% for RC walltype buildings, and 1%, 2%, and 4% for RC frame buildings, respectively. However, there is no information about the interstory drift capacities for wallframe buildings. Jeong et al. [63] used the limit states of ductile walls suggested in Reference [64] for fragility analysis of wallframe buildings in their study; these interstory drift capacities were 0.7%, 1.5%, and 2.5% for the IO, LS, and CP performance levels, respectively. In this section, an attempt has been made to obtain a more precise estimate of the interstory drift capacities for the wallframe buildings used in this study. For this purpose, a damage index concept was used based on the idea that the state of structural damage due to seismic loads can be quantified by a mathematical indicator called the damage index (). The most commonly used damage index is the Park–Ang model [65] that is expressed by deformation and hysteretic energy due to an earthquake excitation. Kunnath et al. [66] recommended a modification for the Park–Ang model in which the deformation terms of the model were replaced by the rotation ones as presented in the following equation:where is the maximum rotation during the loading time, is the ultimate rotation capacity, is the recoverable rotation when unloading, is the yield moment, and is the dissipated energy of each structural component.
for an individual story can be computed using the weighting factor based on hysteretic energy dissipated at the component () as shown in the following equation:
Park and Ang [65] classified damage levels into five categories defined in Table 8. According to this table, the interstory drift ratios in which the damage indices reach the limits of 0.25, 0.4, and 0.8 represent the interstory drift capacities for IO, LS, and CP performance levels, respectively. So, knowing this fact, in this paper, the following procedure was used to find the interstory drift capacities.

First, IDA curves were generated in the plane for each earthquake record. Then, based on nonlinear analysis data, the damage index in each intensity level () was calculated and the damage index curves versus were obtained for all ground motion records, and the intensity levels () corresponding to DI = 0.25, 0.4, and 0.8 were extracted from damage index curves. Finally, the interstory drift values corresponding to these were determined from IDA curves. For example, this procedure is shown for the symmetric model in Figure 9 and Table 9.
(a)
(b)

The same procedure was used for each structural model, and the average of interstory drift capacities for 20 ground motion records in each performance level was calculated (Table 10). Then, to have one unique criterion for the performance evaluation of all models, the average of the interstory drift capacities of the 10 models, as calculated in the last row of Table 10, was used as the interstory drift capacities for three performance levels. As it can be seen from the table, the interstory drift ratios of 0.9%, 1.6%, and 2.9% were estimated as the IO, LS, and CP limitstate capacities, respectively, for the wallframe structural models used in this study. It is worth noting that these values are between the two corresponding aforementioned interstory drift values proposed by ASCE 4113 [27] for shear walls and RC momentresisting frames in each limit state.

4.3. Median IDA Curves
For investigating the effects of strength distribution on the seismic response of the structures, the median IDA curves of all models (with different values) were derived, as shown in Figure 10. As expected, it is apparent from Figure 10 that the median IM capacities of all asymmetric models are reduced over the whole range of the response compared to the symmetric model. However, in most cases, the reduction effects of the flexiblebase asymmetric models on the IM capacities are slightly higher than the fixedbase ones. In both fixed and flexiblebase conditions, the asymmetric models with β = 0.0, 0.25, and 0.5 have better performance than the other asymmetric models, which is also expected since these models have minimum strength eccentricity and balanced configurations. Furthermore, the largest IM capacity reduction happens in the asymmetric model with β = −1.0, which has maximum and compared to other asymmetric cases, and its stiffness and strength centers are both located on one side of the mass center.
(a)
(b)
The flexible effect of soil can be examined considering the curves shown in Figure 11. The three common performance levels, IO, LS, and CP, were also defined on the median IDA curves of this figure. Based on the results of the previous section, the IO and LS points are placed where reaches 0.9% and 1.6%, respectively, and the CP performance limit is violated at θ_{max} = 2.9%. For a better comparison, the ratios of IM capacities of the models with flexiblebase conditions to those of the corresponding model with fixed base are extracted from Figure 11 and shown in Figure 12. It is clear from this figure that the IM capacities corresponding to all performance levels of the flexiblebase models reduce with regard to those of the fixedbase ones. The maximum reduction rate is about 11% which is observed in the asymmetric models with β = −1.0 and 1.0 for the IO performance level.
5. Evaluation of Structural Performance
The data from the IDA study are employed in the probabilistic seismic performance evaluation of structural models, which characterizes the reliability of the models under the seismic hazard levels [41, 67–69]. Based on this approach, the mean annual frequencies of exceeding structural performance levels and the confidence levels in satisfying performance objectives were used to quantify the effects of different strength distributions and also soil flexibility on the probabilistic performance of the structural models considered in this study. For this purpose, first a seismic hazard curve appropriate for this study is defined, and then the abovementioned quantities are evaluated in the next sections.
5.1. Seismic Hazard Curve
In this study, the seismic hazard curve defined in Figure 13 was used in the evaluation of the seismic performance of the models. This curve has been developed from the information of the base hazard curve of Tehran used in Reference [70]. The latter hazard curve, which has been obtained from the probabilistic seismic hazard analysis (PSHA) method, is for the elastic 5%damped firstmode spectral acceleration (), while we needed a curve corresponding to the optimal IM used in this study (). In order to address this issue, a similar procedure to the one used in the previous paper of the author [56] was applied in the current study as well. Based on this procedure, the “standard design response spectrum” of the Iranian Seismic Design Code [22] for the considered site (Tehran) was used to infer from . According to this code, the design spectral acceleration is defined asin which A is the sitespecific earthquake acceleration which is equal to 0.35 g for Tehran, I is the importance factor which is equal to 1.0 for residential buildings with “moderate importance,” and B is the building response factor which can be defined by the following formulas:where T_{0}, T_{s}, S, and S_{0} are the parameters determined from the soil type and seismicity level.
As mentioned before, the structural models are assumed to be located in Tehran with very high seismicity and the underlying soil is classified as Type III. Therefore, according to the Iranian Seismic Design Code [22], the values of , and are equal to 0.15, 0.7, 1.1, and 1.75, respectively.
N is the spectrum modification factor which is defined as
Based on the above formulas, the “standard design response spectrum” for the considered site is derived and shown in Figure 14.
By using the information of the design spectrum and considering the fact that all points on that have the same probability of occurrence, we can define a conversion factor of 0.835 between and and calculate the value of corresponding to each by equation (17). Final based hazard curve used in the study is illustrated in Figure 13.
5.2. The Mean Annual Frequencies of Different Performance Limit States
In this section, the effect of strength distribution on the mean annual frequency (MAF) of exceeding the IO, LS, and CP performance levels in both fixed and flexiblebase models is evaluated. Here, the equations adopted by the Pacific Earthquake Engineering Research Center [60] are used to estimate the MAFs. In this framework, the general form of MAF is illustrated aswhere is the cumulative distribution function (CDF) of the IM for a limitstate capacity and is the spectral acceleration hazard.
By fitting a powerlaw function between the interstory drift demands and the spectral accelerations (as in equation (10)) and then considering a similar simplification for at each level of IM (), it is possible to provide a closedform analytical solution for equation (18) and derive a straightforward equation for the MAF in the IMbased approach as follows [60]:
In this equation, is the median of the IM capacity, is the logarithmic slope of the locally approximated hazard curve at , and is the standard deviation of the natural logarithm of the IM capacity. In this process, and k are calculated by using the S_{A,avg}based hazard curve of Tehran illustrated in Figure 13.
It is also worth mentioning that equation (19) has so far been used by several researchers (e.g., Cordova et al. [49]; Jalayer and Cornell [69]; Pirizadeh and Shakib [56]; and Shakib and Homaei [37]) as a closedform solution for calculating . Additionally, the accuracy of the closedform solutions has also been examined by Medina [71] and Jalayer and Cornell [69], and it was observed that the results of the closedform expression are reasonably close to those of the numerical integration.
In this section, the mean annual frequencies of exceeding the three common performance levels calculated by equation (19) for the symmetric and asymmetric models with fixed and flexible base are presented in Tables 11 and 12. Furthermore, Figure 15 shows the ratio of MAFs of asymmetric models to the symmetric one. According to the results, the following observations can be made:(i)In most cases, the asymmetric models have larger MAF values in comparison with the symmetric model. Also, in the models with SSI effect, the difference between the MAF of the symmetric and asymmetric models is larger than in the fixedbase condition.(ii)In the IO and LS performance levels, except for models with β = 0.5 and 0.25, the MAF values for the flexiblebase asymmetric models increase by up to 60% relative to the symmetric model. The increase rate is about 40% in the models with fixed base.(iii)The increase of MAF is more remarkable for the CP performance level. At this performance level, the MAF values for the asymmetric model with β = −1.0 (having maximum stiffness eccentricity) are approximately twice the MAF values for the symmetric model. After this case, the model with β = −0.75 has the maximum MAF value among the asymmetric models.(iv)The minimum MAF values among the asymmetric models are highlighted in Tables 11 and 12. It can be seen that for the IO and LS performance levels, the smallest MAFs belong to the asymmetric models with β = 0.5 and 0.25, respectively. As mentioned before, these models have balanced and near balanced configurations of CS and CR, respectively. However, for the CP performance level, the asymmetric model with zero strength eccentricity (β = 0.0) has the smallest MAF value, and its behavior is closer to the symmetric model. It is notable that with regard to optimum strength distributions, the trends of the strength distribution effect are similar for both cases of fixed and flexiblebase conditions, and the same results are obtained from both conditions.
(a)
(b)
The ratio of the MAFs of the flexiblebase models to the fixedbase ones is shown in Figure 16. It is apparent that the flexible effect of soil significantly increases the MAF values of the models in most cases. The highest rate of increase is observed for the CP performance level within the range of 10% to 45%. However, the increase is reduced to up to 20% for the IO and LS performance levels.
5.3. The Confidence Level of Meeting the Performance Levels
In a probabilistic performance evaluation format, it is possible to examine the building performance by estimating a level of confidence of the building about the desired performance objective being met [68]. This format takes into account the randomness and uncertainty due to the seismic hazard and the structural capacity in the analysis. According to the following equations, the confidence level for satisfying a specific performance objective () can be evaluated based on the standard Gaussian approach:where and are the median drift capacity and median drift demand, respectively, which can be derived from the IDA method. and are the factors that characterize the randomness in the prediction of the capacity and demand, respectively. Additionally, is defined as the dispersion parameter which considers the epistemic uncertainty in the estimation of demand and capacity, and and are the standard deviation of the natural logarithmic of the drift capacity and demand, respectively. Finally, is a parameter that reflects the powerlaw relationship between demand and IM.
In this section, the confidence levels for meeting the performance levels are investigated for both the fixed and flexiblebase structural models. For this purpose, the confidence levels for satisfying the performance levels are calculated for different intensity levels () (see Figures 17 and 18). It can be observed that in many cases in both the fixed and flexiblebase conditions, the symmetric model has higher confidence levels than the asymmetric models. For the IO and LS performance levels, the asymmetric model with β = −1.0 (with maximum stiffness eccentricity) has the lowest confidence levels. However, for the CP performance level, the lowest confidence is achieved in the asymmetric model with β = 1.0 (with maximum strength eccentricity and zero stiffness eccentricity).
(a)
(b)
(c)
(a)
(b)
(c)
For a more complete comparison, the IM capacities corresponding to 70% confidence in satisfying each limit state are extracted from the figures and shown in Table 13. The higher IM values, highlighted in the table, reflect better probabilistic performance. Similar to what was seen in the previous section for the MAF values, for both the fixed and flexiblebase conditions, the asymmetric models with β = 0.5 and 0.25 have the best performance in the IO and LS limit state, respectively, but for the CP limit sate, the asymmetric model with no strength eccentricity (β = 0.0) performs better. Furthermore, except for the models with β = 0.5 and 0.25, the IM capacities of the asymmetric models corresponding to CL = 70% are reduced by 2 to 30% relative to that of the symmetric model. It is also observed that the IM capacities of the flexiblebase models corresponding to CL = 70% are decreased by up to 15% in comparison with those of the fixedbase models.
6. Conclusion
This study used the incremental dynamic analysis (IDA) method to investigate the effect of different strength distributions on the seismic performance of singlestory planasymmetric reinforced concrete buildings with wallframe systems in the framework of probabilisticbased earthquake engineering. Two different base conditions (fixed and flexible) were considered for the structural models, and an attempt was made to find the most efficient strength distribution for meeting each performance level in each base condition. Firstly, an optimal intensity measure (IM) was selected for the structures. Then, the interstory drift capacities of the structures for three common performance levels (immediate occupancy, life safety, and collapse prevention) were estimated using the damage index concept. Finally, the probabilistic performance of the structural models was evaluated in terms of the mean annual frequency of exceeding different performance limit states and their confidence level. The following main conclusions were drawn:(1)According to the result derived from the investigation of efficiency and sufficiency of 30 candidate IMs, the average was selected as the optimal IM for the considered structural models.(2)The interstory drift ratios of 0.9%, 1.6%, and 2.9% were estimated as the appropriate IO, LS, and CP limitstate capacities, respectively, for the wallframe structural models of this study.(3)The nonuniform distribution of strength and stiffness properties in the plan has remarkable effects on the seismic performance of the buildings, and by using an appropriate strength distribution, the negative effects of planasymmetry can be reduced.(4)The optimum strength distribution is a function of the performance level objective. In the IO performance level, a desirable strength distribution is the one which results in a balanced configuration of the centers of stiffness and strength, where the eccentricities of both parameters are almost equal (β = 0.5). In the LS performance level, the best probabilistic performance could be achieved in the asymmetric model with a nearly balanced configuration (β = 0.25) in which the strength eccentricity is smaller than the stiffness eccentricity. However, in the CP performance level, the asymmetric model with zero strength eccentricity (β = 0.0) performs better. These optimum strength distributions can be adopted in the performancebased design of asymmetric buildings.(5)The lowest confidence levels of meeting performance levels were observed in the asymmetric models with β = −1.0 and −0.75 (with the center of stiffness being on the same side as the center of strength) and β = 1.0 and 0.75 (having zero or almost zero stiffness eccentricity).(6)Soilstructure interaction effect reduces the IM capacities corresponding to a given confidence level in meeting performance objectives. This implies that ignoring the effect of soil flexibility in the design practice of the asymmetric buildings can result in overestimating the seismic capacity, and this phenomenon should be taken into consideration in the design of this type of buildings.(7)The introduction of soil flexibility in the models increases the mean annual frequency of exceedance of the performance levels. The maximum rate of increase is about 20% for the IO and LS performance levels. However, this rate is observed within the range of 10% to 45% for the CP performance level.(8)The effect of soil flexibility does not influence the optimum strength distributions for improving the seismic performance of wallframe asymmetric models, and the optimum configurations are the same for the structures with fixed and flexible base.
Data Availability
The data used during the study are included in the submitted article, and the software codes used to support the findings are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The authors wish to express their sincere gratitude to Professor Ghyslaine McClure from McGill University for her valuable and constructive comments which remarkably improved the article. Her willingness to give her time so generously has been very much appreciated.
References
 T. Paulay and M. J. N. Priestley, Seismic Design of Reinforced Concrete and Masonry Buildings, John Wiley & Sons, New York, NY, USA, 1992.
 M. N. Fardis, Advances in PerformanceBased Earthquake Engineering (ACES Workshop), Springer, Dordrecht, Netherlands, 2010.
 T. Paulay, “Some design principles relevant to torsional phenomena in ductile buildings,” Journal of Earthquake Engineering, vol. 5, no. 3, pp. 273–308, 2001. View at: Publisher Site  Google Scholar
 M. Aschheim, “Seismic design based on the yield displacement,” Earthquake Spectra, vol. 18, no. 4, pp. 581–600, 2002. View at: Publisher Site  Google Scholar
 M. J. N. Priestley, G. M. Calvi, and M. J. Kowalsky, DisplacementBased Seismic Design of Structures, IUSS Press, Pavia, Italy, 2007.
 W. K. Tso and H. Ying, “Lateral strength distribution specification to limit additional inelastic deformation of torsionally unbalanced structures,” Engineering Structures, vol. 14, no. 4, pp. 263–277, 1992. View at: Publisher Site  Google Scholar
 W. K. Tso and R. S. H. Smith, “Reevaluation of seismic torsional provisions,” Earthquake Engineering & Structural Dynamics, vol. 28, no. 8, pp. 899–917, 1999. View at: Publisher Site  Google Scholar
 T. Paulay, “Seismic torsional effects on ductile structural wall systems,” Journal of Earthquake Engineering, vol. 1, no. 4, pp. 721–745, 1997. View at: Publisher Site  Google Scholar
 W. K. Tso and B. Myslimaj, “A yield displacement distributionbased approach for strength assignment to lateral forceresisting elements having strength dependent stiffness,” Earthquake Engineering & Structural Dynamics, vol. 32, no. 15, pp. 2319–2351, 2003. View at: Publisher Site  Google Scholar
 B. Myslimaj and W. K. Tso, “A designoriented approach to strength distribution in singlestory asymmetric systems with elements having strengthdependent stiffness,” Earthquake Spectra, vol. 21, no. 1, pp. 197–212, 2005. View at: Publisher Site  Google Scholar
 H. Shakib and A. Ghasemi, “Considering different criteria for minimizing torsional response of asymmetric structures under nearfault and farfault excitations,” International Journal of Civil Engineering, vol. 5, no. 4, pp. 247–265, 2007. View at: Google Scholar
 A. Aziminejad and A. S. Moghadam, “Performance of asymmetric multistory shear buildings with different strength distributions,” Journal of Applied Sciences, vol. 6, no. 9, pp. 1082–1089, 2009. View at: Publisher Site  Google Scholar
 A. Aziminejad and A. S. Moghadam, “Fragilitybased performance evaluation of asymmetric singlestory buildings in near field and far field earthquakes,” Journal of Earthquake Engineering, vol. 14, no. 6, pp. 789–816, 2010. View at: Publisher Site  Google Scholar
 T. Balendra, C. W. Tat, and S.L. Lee, “Modal damping for torsionally coupled buildings on elastic foundation,” Earthquake Engineering & Structural Dynamics, vol. 10, no. 5, pp. 735–756, 1982. View at: Publisher Site  Google Scholar
 K. S. Sivakumaran and T. Balendra, “Seismic analysis of asymmetric multistorey buildings including foundation interaction and PΔ effects,” Engineering Structures, vol. 16, no. 8, pp. 609–624, 1994. View at: Publisher Site  Google Scholar
 W.H. Wu, J.F. Wang, and C.C. Lin, “Systematic assessment of irregular buildingsoil interaction using efficient modal analysis,” Earthquake Engineering & Structural Dynamics, vol. 30, no. 4, pp. 573–594, 2001. View at: Publisher Site  Google Scholar
 H. Shakib and A. Fuladgar, “Dynamic soilstructure interaction effects on the seismic response of asymmetric buildings,” Soil Dynamics and Earthquake Engineering, vol. 24, no. 5, pp. 379–388, 2004. View at: Publisher Site  Google Scholar
 R. Roy and S. Chandra Dutta, “Inelastic seismic demand of lowrise buildings with soilflexibility,” International Journal of Nonlinear Mechanics, vol. 45, no. 4, pp. 419–432, 2010. View at: Publisher Site  Google Scholar
 H. Shakib and G. Atefatdoost, “Minimizing the torsional response of asymmetric walltype systems considering soilstructure interaction,” International Journal of Civil Engineering, vol. 12, no. 1A, pp. 14–24, 2014. View at: Google Scholar
 G. R. Atefatdoost, B. JavidSharifi, and H. Shakib, “Effects of foundation flexibility on seismic demands of asymmetric buildings subject to nearfault ground motions,” Structural Engineering and Mechanics, vol. 66, no. 5, pp. 637–648, 2018. View at: Google Scholar
 T. Paulay, “An estimation of displacement limits for ductile systems,” Earthquake Engineering & Structural Dynamics, vol. 31, no. 3, pp. 583–599, 2002. View at: Publisher Site  Google Scholar
 Building and Housing Research Centre, Iranian Code of Practice for Seismic Resistant Design of Buildings (Standard 2800), Building and Housing Research Centre, Tehran, Iran, 4th edition, 2014.
 Pacific Earthquake Engineering Research Center, Open System for Earthquake Engineering Simulation OpenSees, V.2.5.2, Pacific Earthquake Engineering Research Center, Berkeley, CA, USA, 2017.
 A. K. Chopra, Dynamics of Structures: Theory and Applications to Earthquake Engineering, PrenticeHall, Englewood Cliffs, NJ, USA, 2007.
 L. F. Ibarra, R. A. Medina, and H. Krawinkler, “Hysteretic models that incorporate strength and stiffness deterioration,” Earthquake Engineering & Structural Dynamics, vol. 34, no. 12, pp. 1489–1511, 2005. View at: Publisher Site  Google Scholar
 H. Krawinkler and F. Zareian, “Prediction of collapsehow realistic and practical is it, and what can we learn from it?” The Structural Design of Tall and Special BuildingsE, vol. 16, no. 5, pp. 633–653, 2007. View at: Publisher Site  Google Scholar
 American Society of Civil Engineers, Seismic Evaluation and Retrofit of Existing Buildings (ASCE/SEI 4113), American Society of Civil Engineers, Reston, VA, USA, 2013.
 Institute of Building National Code, Iranian National Building Code, Part 7Foundation Design, Institute of Building National Code, Tehran, Iran, 2014.
 F. Behnamfar and M. Banizadeh, “Effects of soilstructure interaction on distribution of seismic vulnerability in RC structures,” Soil Dynamics and Earthquake Engineering, vol. 80, pp. 73–86, 2016. View at: Publisher Site  Google Scholar
 S. Gajan, T. Hutchinson, B. L. Kutter, P. Raychowdhury, and A. Ugalde, “Numerical models for analysis and performancebased design of shallow foundations subjected to seismic loading,” Tech. Rep., Pacific Earthquake Engineering Research Center (PEER), University of California, Berkeley, CA, USA, 2007, Technical Report 2007/04. View at: Google Scholar
 P. Raychowdhury and T. Hutchinson, “Nonlinear material models for Winklerbased shallow foundation response evaluation,” in Proceedings of the Geo Congress, New Orleans, LA, USA, March 2008. View at: Publisher Site  Google Scholar
 C. W. Harden, T. C. Hutchinson, G. R. Martin, and B. L. Kutter, “Numerical modeling of the nonlinear cyclic response of shallow foundations,” Tech. Rep., Pacific Earthquake Engineering Research Center (PEER), University of California, Berkeley, CA, USA, 2005, Technical Report 2005/04. View at: Google Scholar
 P. Raychowdhury, “Seismic response of lowrise steel momentresisting frame (SMRF) buildings incorporating nonlinear soilstructure interaction (SSI),” Engineering Structures, vol. 33, no. 3, pp. 958–967, 2011. View at: Publisher Site  Google Scholar
 M. A. Ghannad and A. H. Jafarieh, “Inelastic displacement ratios for soilstructure systems allowed to uplift,” Earthquake Engineering & Structural Dynamics, vol. 43, no. 9, pp. 1401–1421, 2014. View at: Publisher Site  Google Scholar
 H. Masaeli, F. Khoshnoudian, and R. Ziaei, “Rocking soilstructure systems subjected to nearfault pulses,” Journal of Earthquake Engineering, vol. 19, no. 3, pp. 461–479, 2015. View at: Publisher Site  Google Scholar
 P. Raychowdhury and S. RayChaudhuri, “Seismic response of nonstructural components supported by a 4story SMRF: effect of nonlinear soilstructure interaction,” Structures, vol. 3, pp. 200–210, 2015. View at: Publisher Site  Google Scholar
 H. Shakib and F. Homaei, “Probabilistic seismic performance assessment of the soilstructure interaction effect on seismic response of midrise setback steel buildings,” Bulletin of Earthquake Engineering, vol. 15, no. 7, pp. 2827–2851, 2017. View at: Publisher Site  Google Scholar
 R. W. Boulanger, C. J. Curras, B. L. Kutter, D. W. Wilson, and A. Abghari, “Seismic soilpilestructure interaction experiments and analyses,” Journal of Geotechnical and Geoenvironmental Engineering, vol. 125, no. 9, pp. 750–759, 1999. View at: Publisher Site  Google Scholar
 https://ngawest2.berkeley.edu/.
 D. Vamvatsikos and C. A. Cornell, “Incremental dynamic analysis,” Earthquake Engineering & Structural Dynamics, vol. 31, no. 3, pp. 491–514, 2002. View at: Publisher Site  Google Scholar
 FEMA, “Recommended seismic design criteria for new steel moment frame buildings,” Tech. Rep., Washington, DC, USA, 2000, Report No. FEMA350. View at: Google Scholar
 FEMA, “Prestandard and commentary for the seismic rehabilitation of buildings,” Washington, DC, USA, 2000, Report No. FEMA356. View at: Google Scholar
 FEMA, “Quantification of building seismic performance factors,” Washington, DC, USA, 2009, Report No. FEMAP695. View at: Google Scholar
 N. Luco and C. A. Cornell, “Structurespecific scalar intensity measures for nearsource and ordinary earthquake ground motions,” Earthquake Spectra, vol. 23, no. 2, pp. 357–392, 2007. View at: Publisher Site  Google Scholar
 A. Elenas and K. Meskouris, “Correlation study between seismic acceleration parameters and damage indices of structures,” Engineering Structures, vol. 23, no. 6, pp. 698–704, 2001. View at: Publisher Site  Google Scholar
 S. Akkar and B. Küçükdoğan, “Direct use of PGV for estimating peak nonlinear oscillator displacements,” Earthquake Engineering & Structural Dynamics, vol. 37, no. 12, pp. 1411–1433, 2008. View at: Publisher Site  Google Scholar
 N. Shome, C. A. Cornell, P. Bazzurro, and J. E. Carballo, “Earthquakes, records, and nonlinear responses,” Earthquake Spectra, vol. 14, no. 3, pp. 469–500, 1998. View at: Publisher Site  Google Scholar
 P. Giovenale, C. A. Cornell, and L. Esteva, “Comparing the adequacy of alternative ground motion intensity measures for the estimation of structural responses,” Earthquake Engineering & Structural Dynamics, vol. 33, no. 8, pp. 951–979, 2004. View at: Publisher Site  Google Scholar
 P. P. Cordova, G. G. Deierlein, S. S. F. Mehanny, and C. A. Cornell, “Development of a two parameter seismic intensity measure and probabilistic assessment procedure,” in Proceedings of the 2nd USJapan Workshop on PerformanceBased Earthquake Engineering Methodology for RC Building Structures, Sapporo, Japan, September 2000. View at: Google Scholar
 D. Vamvatsikos and C. A. Cornell, “Developing efficient scalar and vector intensity measures for IDA capacity estimation by incorporating elastic spectral shape information,” Earthquake Engineering & Structural Dynamics, vol. 34, no. 13, pp. 1573–1600, 2005. View at: Publisher Site  Google Scholar
 A. K. Kazantzi and D. Vamvatsikos, “Intensity measure selection for vulnerability studies of building classes,” Earthquake Engineering & Structural Dynamics, vol. 44, no. 15, pp. 2677–2694, 2015. View at: Publisher Site  Google Scholar
 S. Tsantaki, C. Adam, and L. F. Ibarra, “Intensity measures that reduce collapse capacity dispersion of Pdelta vulnerable simple systems,” Bulletin of Earthquake Engineering, vol. 15, no. 3, pp. 1085–1109, 2017. View at: Publisher Site  Google Scholar
 M. Kohrangi, P. Bazzurro, and D. Vamvatsikos, “Vector and scalar IMs in structural response estimation, part I: hazard analysis,” Earthquake Spectra, vol. 32, no. 3, pp. 1507–1524, 2016. View at: Publisher Site  Google Scholar
 M. Kohrangi, P. Bazzurro, and D. Vamvatsikos, “Vector and scalar IMs in structural response estimation, part II: building demand assessment,” Earthquake Spectra, vol. 32, no. 3, pp. 1525–1543, 2016. View at: Publisher Site  Google Scholar
 M. Bianchini, P. P. Diotallevi, and J. W. Baker, “Prediction of inelastic structural response using an average of spectral accelerations,” in Proceedings of the 10th International Conference on Structural Safety and Reliability, Osaka, Japan, September 2009. View at: Google Scholar
 M. Pirizadeh and H. Shakib, “Evaluation of efficient ground motion intensity measure for vertical irregular structures by aims of reducing the dependency of results to record selection,” in Proceedings of the 6th International Conference on Seismology and Earthquake Engineering, Tehran, Iran, May 2011. View at: Google Scholar
 L. Eads, E. Miranda, and D. G. Lignos, “Average spectral acceleration as an intensity measure for collapse risk assessment,” Earthquake Engineering & Structural Dynamics, vol. 44, no. 12, pp. 2057–2073, 2015. View at: Publisher Site  Google Scholar
 M. A. A. AlBaadani and M. AlShameri, “Study of dispersion of the ground motion intensity measure of R/C frame minaret,” MUST Journal of Engineering and Applied Sciences, vol. 1, no. 1, pp. 24–33, 2017. View at: Publisher Site  Google Scholar
 P. Colajanni, S. Pagnotta, and G. Testa, “Comparison of fully nonstationary artificial accelerogram generation methods in reproducing seismicity at a given site,” Soil Dynamics and Earthquake Engineering, vol. 133, Article ID 106135, 2020. View at: Publisher Site  Google Scholar
 C. A. Cornell, F. Jalayer, R. O. Hamburger, and D. A. Foutch, “Probabilistic basis for 2000 SAC federal emergency management agency steel moment frame guidelines,” Journal of Structural Engineering, vol. 128, no. 4, pp. 526–533, 2002. View at: Publisher Site  Google Scholar
 J. E. Padgett and R. DesRoches, “Methodology for the development of analytical fragility curves for retrofitted bridges,” Earthquake Engineering and Structural Dynamics, vol. 37, no. 13, pp. 1157–1174, 2008. View at: Publisher Site  Google Scholar
 H. Shakib and V. Jahangiri, “Intensity measures for the assessment of the seismic response of buried steel pipelines,” Bulletin of Earthquake Engineering, vol. 14, no. 4, pp. 1265–1284, 2016. View at: Publisher Site  Google Scholar
 S.H. Jeong, A. M. Mwafy, and A. S. Elnashai, “Probabilistic seismic performance assessment of codecompliant multistory RC buildings,” Engineering Structures, vol. 34, pp. 527–537, 2012. View at: Publisher Site  Google Scholar
 A. Ghobarah, N. M. Aly, and M. ElAttar, “Seismic reliability assessment of existing reinforced concrete buildings,” Journal of Earthquake Engineering, vol. 2, no. 4, pp. 569–592, 1998. View at: Publisher Site  Google Scholar
 Y. J. Park and A. H. S. Ang, “Mechanistic seismic damage model for reinforced concrete,” Journal of Structural Engineering, vol. 111, no. 4, pp. 722–739, 1985. View at: Publisher Site  Google Scholar
 S. K. Kunnath, A. M. Reinhorn, and R. F. Lobo, “IDARC version 3.0: a program for the inelastic damage analysis of reinforced concrete structures,” Tech. Rep., State University of New York, New York, NY, USA, 1992, NCEER920022, National Center for Earthquake Engineering Research. View at: Google Scholar
 P. Bazzurro and C. A. Cornell, “Seismic hazard analysis of nonlinear structures. I: methodology,” Journal of Structural Engineering, vol. 120, no. 11, pp. 3320–3344, 1994. View at: Publisher Site  Google Scholar
 F. Jalayer and C. A. Cornell, “A technical framework for probabilitybased demand and capacity factor design (DCFD) seismic formats,” Tech. Rep., University of California, Berkeley, CA, USA, 2003, PEER2003/08, Pacific Earthquake Engineering Research Center. View at: Google Scholar
 F. Jalayer and C. A. Cornell, “Alternative nonlinear demand estimation methods for probabilitybased seismic assessments,” Earthquake Engineering & Structural Dynamics, vol. 38, no. 8, pp. 951–972, 2009. View at: Publisher Site  Google Scholar
 B. Asgarian, A. Sadrinezhad, and P. Alanjari, “Seismic performance evaluation of steel moment resisting frames through incremental dynamic analysis,” Journal of Constructional Steel Research, vol. 66, no. 2, pp. 178–190, 2010. View at: Publisher Site  Google Scholar
 R. Medina, “Seismic demands for nondeteriorating frame structures and their dependence on ground motions,” Department of Civil and Environmental Engineering, Stanford, CA, USA, 2002, Ph.D. thesis. View at: Google Scholar
Copyright
Copyright © 2020 Sahar Mohammadzadeh Osalu and Hamzeh Shakib. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.