Seismic performance assessment of structures is a process in which performance indicators, namely decision variables (i.e. financial losses, downtime and human casualties) are specified and compared with predefined performance objectives. The promising procedure, namely FEMA-58 method , aims to predict the decision variables where various sources of uncertainties are involved through a consistent probabilistic manner. Towards being thoroughly applicable for developing countries, the required data, existing tools and challenges to achieve the final results must be investigated.
According to the force-based seismic design regulations, the earthquake shocks have been modeled as induced inertia forces to the structures . Conversely, the permission of occurrence of structural damages against severe earthquakes, with low probability of exceedance and high consequences, motivates researches to replace seismic-induced displacement with that of inertia force as the design parameter. Performance-based earthquake engineering aims to insure the desired overall structural performance through restricting the seismic-induced forces (for the brittle structural members) and the seismic-induced displacements (for the ductile structural members) against various earthquake hazard levels. Several standards have been developed according to the performance-based earthquake engineering philosophy (namely, CEN-2005 , JGS4001-2004 , AS/NZS-1170-2002  and …). A large amount of seismic-induced economic losses during recent earthquakes revealed limitations of the first generation performance-based regulations . These challenges consist of the structural performance expression using non-comprehensible variables for all stakeholders and the large amount of inherent uncertainties in the predicted performance which have not been involved explicitly in the design process. According to the mentioned shortcomings, the next generation performance-based earthquake engineering guideline was presented. The promising procedure, namely the FEMA P-58, aims to express the structural performance using the comprehensible variables (i.e. seismic-induced economic losses, number of casualties and downtime) for all contributed stakeholders. Furthermore, the inherent uncertainties are involved explicitly. Through the FEMA P-58 procedure, the uncertainties of performance indicators (i.e. Decision Variables) are comprised of strong ground motion Intensity Measure (IM) uncertainty, structural seismic response (Engineering Demand Parameter: (EDP)) uncertainty, and seismic-induced damages (DM). Uncertainty of strong ground motion intensity measure is presented by probabilistic seismic hazard analysis. Probabilistic Seismic Hazard Analysis (PSHA) shows the Mean Annual Frequency of exceedance (MAFE) for the selected IM. Uncertainty of EDP is characterized through Incremental Dynamic Analysis (IDA) procedure, through which the structure is excited by several strong ground motions. Uncertainty of seismic-induced damage is presented by the fragility curves. Fragility curves corresponding to the local and global limit states present the probability of the relevant damage conditioned on the affected EDP . Uncertainty of Decision Variable (DV) (consists of economic losses, number of casualties and downtime) is characterized by loss functions. For each considered damage states, the relevant loss functions show the probability of DV conditioned on the activated damage . Mathematical formulation for estimating the mean annual frequency of decision variables (DV) is presented in Equation (1) which is formulated according to the conditional probability theorem through which the contributed uncertainties due to IM, EDP, DM and DV are involved. The integrations are implemented over Damage Measure (DM); engineering demand parameter (EDP), and intensity of ground motion (IM) variables. In Equation (1) the function G (X|Y) shows the cumulative distribution function of X conditioned on Y. Equation (1) is used to calculate the mean annual frequency of occurrence.
The FEMA-P 58 seismic performance-based methodology is a 10-year-program which aims to develop the seismic design guidelines which will be applicable for new and existing buildings. Furthermore, the methodology is in consistent with the explosion engineering, the fire engineering and the wind engineering regulations. The project is still under development by FEMA team . The seismic performance assessment will be implemented through the intensity-based, the scenario-based and the time-based procedures. The proposed decision variables are repair cost, number of casualties and downtime .
Available tool: Performance Assessment Calculation Tool (PACT software) is one of the tools developed based on FEMA-P 58 formulation. Besides, the necessary supporting data (such as fragility curves and relevant loss functions) to calculate the seismic performance of the buildings are presented in PACT. To use the package, the building general information(such as number of stories, total replacement cost, floor area, story height, building occupancy, etc.), population model (i.e. number of inhabitants), component fragility functions of the building, performance groups hazard curves of the located zone and structural analysis results (i.e. structural seismic response values) are the inputs to the package. The loss curves (which show the probability of loss values) are the output . In this software, it is possible to make the appropriate modifications to the input variables according to the existing features to achieve more accurate results. First of all, performance parameters should be specified. Then, the collapse fragility curve should be specified. This curve is calculated based on Incremental Dynamic Analysis, through which the structure is excited by number of strong ground motions with increasing intensities. The collapse points are characterized as the intensity for which the numerical model will not converge applying any of the existing algorithms. The probability distribution fitted to the collapse points presents the collapse fragility curve. As the next step, the structural analysis results are required. Through this step, vectors of structural demands (such as inter-story drift ratio, floor acceleration, floor velocity, etc.) should be specified against number of strong ground motions scaled to the various intensities. Based on the input demand vectors, the corresponding probabilistic model is generated and hundreds of realizations are produced applying the Monte Carlo simulation method. Using the previously-defined fragility curves of the damage states, the simulated demand vectors present the damage probabilities. Finally, the loss values and the relevant exceedance probabilities are summed up for all damage states  . Applying the collapse fragility curve and the exceedance probability of loss, the final loss value, condition on the affected IM, is estimated based on Equation (2). This equation is used to calculate the total annual losses for both collapsed and non-collapsed conditions using the general probabilistic-statistical laws.
In Equation (2), , represents total loss of the structure (including repair cost, downtime, casualties). The CO and NC stands for collapse and non-collapse states. Therefore, shows the exceedance probability of loss conditioned the affected IM assuming the collapse doesn’t occur. The calculation of this probability is explained in the previous paragraph. is the probability of NC conditioned on IM, which is one minus the collapse fragility curve. is the exceedance probability of loss while the structure collapses. This function can be assumed as the probability of building replacement cost. Finally, is the collapse fragility curve. To obtain the mean annual frequency of exceedance (MAFE) of loss, the calculated loss curve (based on Equation (2)) is integrated over the probabilistic seismic hazard curve, which shows the MAFE for the selected IM (normally Sa(T1)). Towards making the performance assessment calculation tool (PACT) applicable for developing countries, the required input data and the relevant challenges should be investigated.
In the present study, the performance assessment calculation tool (PACT) is applied to study the seismic performance of sampled structures. The sampled structures are representatives of typical structures in developing countries. Diagonal braces, ordinary moment-resisting frame and special moment-resisting frame are selected as typical structures in developing countries. In developing countries, we face to lack lots of information about Seismic design methodology and loss models and seismic hazard. So, in this paper, we try to implement FEMA methodology for these countries. Because of high seismic potential of these countries, the FEMA has been selected. Seismic performances of structures are calculated based on FEMA P-58 methodology and by PACT tool. As the inputs of PACT tool, probabilistic seismic hazard model is presented based on power function whose constants are calibrated based on empirical seismic hazard data. Probabilistic seismic demand model is calculated based on Incremental Dynamic Analysis of the sampled frames against number of earthquake strong ground motions. Fragility functions corresponding to two commonly used limit states (i.e. Immediate Occupancy and Collapse Prevention) are estimated using the IM-based method based on IDA results. Calculated seismic loss curves will present the powerful tool to compare the seismic performance of the sampled structures and to make rational decisions if the reliable prerequisite input data are provided.
2. Modeling and Assumptions
To represent the results of analytical studies, three random structural frames were selected from three separated steel buildings (a braced frame and two special moment resisting frames with high and low ductility. The frames were loaded based on UBC97 code  and designed based on AISC 360  as depicted in Table 1. To validate results, a two-dimensional (2D) model was created as shown in Figure 1. Then, one of the side frames was analyzed as the representative of the buildings to take into account some design parameters, such as biaxial bending for designing of the columns . In this paper, in order to
Table 1. Used sections for braced and special moment frames (with different ductility).
Figure 1. Two dimensional models: (a) braced frame (b) special moment frame (with different ductility) (c) plan view of case study buildings caption).
achieve different ductility, two types of connections were used in existing buildings in Northridge (named as Pre-Northridge and Post-Northridge Connections) to achieve two different ductilities of frames . Common frame connections before Northridge earthquake were mostly welded joints of WUF, but majority of designers preferred to use RBS connections and bolted connections after the earthquake for residential and/or office buildings with medium hazard level .
It should be mentioned that columns with box sections in the moment frames are selected to achieve easier modeling . Moreover, the value of behavior factor for the structure with special moment-resisting frame (SMRF) is equal to 8.5 and 7.5 for buildings with dual system and with special moment frame System. It is worth noting, according to the UBC97 code, buildings with usual moment frame system are not suggested for located in high risk zones. So, three different types of connections were used in modeling of above buildings to achieve different ductility. The first type is center to center non-linear connections used in dual braced frame modeling. Since structural stiffness and strength are two major criteria used to design new buildings and/or used to assess existing buildings. It is recommended that center to center non-linear modeling be used for dual braced frames. Second type of connection, used in the frames with high and low ductility, is non-linear model with spring connection in which beam and column dimensions are regarded as connection spring dimensions. In this type of modeling, a three-linear torsional spring and two double-linear torsional springs can be used in one of the connecting corners to assign a three-linear behavior .
Regarding to convergent bracing frames (CBFs) have numerous applications in buildings as lateral load systems, seismic performance of these bracings depends on their hysteresis cycle. In connection modeling, fiber model was used for braced connections’ modeling and center to center non-linear connections were used for modeling of braced frames modeling.
To model low and high ductilities connections in Open Sees software, plastic hinge length was used as the criteria to make rigid and semi-rigid connections of such buildings (Figure 2). Support columns were used to take into account the effect of gravitational load on interior moment frames and applied demands on exterior moment frames . These columns are, in fact, main rigid members which are connected to the main structure by rigid beams in both sides of the
Figure 2. Schematic model of panel-zone elements: (a) Panel-Zone element used in pre-Northridge buildings (corresponding to WUF connection for low ductility buildings); (b) Panel-zone element used in post-Northridge buildings (corresponding to RBS connection for high ductility buildings); (c) Schematic connection modeling in the braced frame.
joints. Moreover, nodes of beam and column sections are made from joints that do not transfer additional moment resistance to the system.
In order to consider P-delta effects, the weights of internal frame spans (due to the symmetry of half of them) are applied as point loads to the end nodes of rigid beams. When the lateral drift is applied to the buildings, these point loads cause secondary bending moments at the moment-resisting frames. Support columns along with other details are illustrated in Figure 3 based on FEMA-P695 . The analytical period of modeled buildings is also shown in Table 2.
Figure 3. Schematic view of support columns modeling in special moment frames.
Table 2. Obtained periods of buildings.
3. Record Selection
The first step in performance assessment method is drawing IDA curves and preparing a series of earthquake ground motions so that it would express regional seismicity . Usually, due to the lack of a sufficient number of recorded earthquakes in the developing countries under consideration, records from other regions are adopted that may not belong to the same seismic plateau, but are similar in terms of earthquake magnitude, source-to-site distance and local geological and subsoil conditions . Here, a limited number of ground motions with similar conditions in terms of seismic characteristics have been used. Thus, sufficient numbers of ground motions, i.e., 22 records in accordance with FEMA P695 recommendations regarding earthquake event and recording station data for the far-field record set, have been used . Magnitude versus distance of selected ground motions diagram has been shown in Figure 4. Also, we can see mean spectra from selected ground motions in Figure 5.
4. Incremental Dynamic Analysis (IDA)
The response of the structure under the ground movements can be estimated by time history dynamic analysis. The sensitivity of responses to the selected
Figure 4. Magnitude versus distance of selected ground motions.
Figure 5. Mean spectra from selected ground motions.
records is one of the most important problems in applying non-linear dynamic analysis . Representing incremental dynamic analysis and estimating responses based on the probability relations can largely decrease this defect. The method of producing IDA curves is exposing one structural modeling subjected to different levels of seismic intensity of one or more earthquake ground motions. Each accelerometer has been measured to reach a special intensity which not only investigates the seismic behavior of buildings but also gives some information to engineers about structural capacities and estimates the seismic performance of buildings . These methods give some information to engineers about structural behaviors, capacity, and seismic performance of buildings by considering structure response to probable intensities of earthquake and structure behavior to strong and rare earthquakes and natural response of the structure to the gradual increase of earthquake intensity and also by estimating dynamic capacity of structural system. The regular form of non-linear incremental dynamic analysis was conducted first to investigate the performance of buildings. Incremental dynamic analysis is performed differently as a parametric analysis to appropriately assess structure performance due to seismic loads . The gradual increase of seismic loads was performed for the first time to investigate structure performance  . However, it was regulated for the first time to investigate structure performance . PEER has done extensive research on risk analysis and has a close relation with ground science and geological engineering to improve the accuracy of common IMs. Moreover, Cornell investigated different options to seismic intensity values which had the best relation with loss from earthquake . Studies by Shome  showed that spectral acceleration for the first mode of the structure and damping of 5%, Sa (T1, 5%), was considered as appropriate selection for intensity measures (IM). Since this paper focuses on general behavior of the buildings and structural elements, the maximum value of relative story drifts is selected as damage measures (DM) . In the following figures, EDP presents a maximum inter-story drift ratio (max.IDR). It is necessary to use the respective data collected from the records from specified zone to investigate non-linear dynamic analysis of buildings in each site. But these ground motions with wide ranges are not accessible. In this paper, incremental dynamic analysis is used for obtaining IDA curves, in which it is possible to obtain required ground motions by scaling them for analyzing. Obtained curves from this method are shown for three frames in Figure 6 and comparing their percentiles of 50% have been illustrated in Figure 7. As it can be seen, the mean IDA curve of frame with low ductility will become horizontal in high seismic intensity measure level.
Figure 6. Obtained IDA curves for special moment frames: (a) For low ductility frame, (b) For high ductility frame, (c) For braced frame.
Figure 7. Median IDA curves for three frames.
5. Fragility Curves
Fragility curves are used to express structural and non-structural system vulnerabilities as well as performance parameters in buildings . Due to the approach toward earthquake engineering based on performance and due to the importance and necessity of appropriate management of probable crisis such as earthquake, it is necessary to assess vulnerability of structural and non-structural systems. Moreover, an appropriate assessment of probable loss and reasonable decision making becomes more important and necessary. In seismic performance-based design methodology, failure or damage does not happen when the structure is unable to reach defined level of performance requirements   and . Fragility curves of collapse level are the most important tools used in the analysis of damage and direct and indirect economic loss effects of earthquakes on buildings. It is possible to observe some obtained fragility curves for structural and non-structural elements   and . Performance requirements are dependent on two connected relations. One of them is related to the relation of seismic intensity of ground movement and destruction probability (destruction fragility curve). The other one is related to the relation between earthquake risk and earthquake intensity of similar ground movement (earthquake risk curve) . In this article, intensity movement is directly used as destruction probability approach based on seismic intensity (IM-based) . Methods to assess performance values of destruction have been brought in SAC/FEMA guidelines  . Therefore, the fragility curves of a structure are considered as fundamental elements to assess destruction performance, known as destruction probability of a structure for seismic intensity of ground movement. For this case, the limited probability of downtime is calculated according to Equation (3):
In Equation (3), is the cumulative function of capacity probability of seismic intensity. If all input parameters and earthquake effects on buildings are certain, the probable function will be one or zero, but in fact, there are some effects that have inherent nature or cause some changes in the capacity parameter of the structure due to the lack of knowledge. Mathematical view is observable in Equation (4) in which IMC is the defined limited critical case for seismic intensity and is the cumulative probability of structural failure for seismic intensity of imi. IMC is corresponding to IM of unstable dynamic mode.
By using obtained curves, it is possible to calculate the probability of limited case in lieu of each level of IM-provided that IM value is limited to the given level. Using lognormal distribution for points in each IM level, 16th, 50th, and 84th percentiles of IDA curves may be extracted for two damage levels of immediate occupancy (IO) level (1% drift) and collapse prevention (CP) level (10% drift exceeding) . It should be mentioned that, according to FEMA 350 . Obtained Fragility curves for three models and two limit state levels are shown in Figure 8. From this figure, it can be obtained that there is higher damage probability in each seismic intensity measure for frame with high ductility.
6. Hazard Curves
In this section, introduced Seismic hazard Curves derived from Mahdavi et al. research . According to available attenuation relationships, different hazard curves for the seismic zones have been expressed. These curves representing Mean Annual Frequency (MAF) of exceedance from different values of seismic intensity, are obtained for numerous structural periods using the uniform risk map . MAF of exceedance from seismic intensity based on Sa is estimated through a linear relation in Log-Log space. The relation related to this estimation is expressed as Equation (5).
k and t parameters, derived from studies about oscillation periods related to the buildings in this article for the area with high level of risk, are brought in Table 3 . Obtained uniform risk curves are illustrated in Figure 9 for each of three frames. T1 is frame’s period in first shape mode which has obtained from software.
Figure 8. Comparing fragility curves obtained from IM-based method  using IDA curves obtained from the common method for two damage levels corresponding to three modeled buildings.
Table 3. Values of k and t parameters corresponding to Equation (5).
Figure 9. Global comparison between obtained uniform hazard curves in all frames.
7. Mean Annual Frequency
These values are useful as the quantities reflecting the general probabilistic capacity of buildings relying on uncertainties from earthquake. These values can be used as the criteria to measure the structural assurance of available buildings in comparison with other buildings or they can be applied in statutory criteria related to structure designs . Equation (6) is used to calculate mean annual frequency.
In Equation (6), λ(collapse) represents Mean Annual Frequency (MAF) of exceedance for IM in which the quantity in absolute value, hazard gradient of IM, and are damage probability or fragility function value which is obtained using high seismic risk curves (because of building’s site location) and fragility curves and numerical integration of Equation (6). Values related to MAF of limit-state levels (IO and CP level) for buildings are illustrated in Table 4. These values are useful as the quantities reflecting the general probabilistic capacity of buildings relying on uncertainties from earthquake. These values can be used as the criteria to measure the structural assurance of available buildings in comparison with other buildings or they can be applied in statutory criteria related to structure designs . It can be seen in Table 4 that the MAF values are higher in low-ductility frames than high-ductility ones, for all two performance levels.
Table 4. Obtained MAF for three buildings.
8. Loss Evaluation
In order to assess the loss, it is necessary to have functional groups, their positions and their values in the structure. Modeled buildings are assessed for this purpose. The ground story is considered as a parking, and higher stories are categorized as residential building. In each story, there are two units. The least considered performance groups for one story of such buildings have been supposed based on FEMA P-58 instruction and condition of the region  . Then, final loss for three buildings is calculated using PACT software. In the software, incremental dynamic analysis results (IDA curves), fragility analysis results (fragility curves) and hazard analysis results (hazard curves) and total cost of building construction per square meters, have been used as the inputs. To implement PACT software, it is necessary to define cost function. These functions are related to contractor attitudes, approaches and procedures in the region. Incremental dynamic analysis results and PACT software are used to calculate non-downtime mode. In order to obtain results for downtime mode (in which major loss for residential buildings occur), cumulative density function, with average cost equal to three times of probable cost repair from structure damage has been supposed with 5% of distribution. By adding the results for above two modes, loss curves are obtained. Figure 10 and Figure 11 show loss curves and annually probabilistic final loss curves for three case study buildings, respectively.
According to Figure 10 and Figure 11, the braced structure has less annual losses in comparison with other investigated cases, and the structure with high ductility can be considered as the next ones. Performing cost of each building should be considered by contractors. Obviously, the cost of execution of structure with high ductility is more than other buildings. Using the results of this paper can be a common language for decision-making and transfer of information between engineers and employers.
Figure 10. Comparing loss curves obtained from three structure.
Figure 11. Comparing annually probabilistic final loss curves.
The present study makes information about economic feasibility for seismic retrofit of existing buildings. In the performance-based design method, it is possible to estimate the behavior of buildings in the case of earthquake. The most important reason of discussion about seismic design based on performance may be encouraging to use innovation in developing some methods to promote performance. In this article, the general method introduced in FEMA P-58 project has been investigated on three buildings. This method helps engineers make rational decisions regarding the selection of global structural components including selection of lateral-resisting systems and optimum use of materials to achieve the desired performance objectives. Implementing this method in developing countries is probabilistic due to some defects in applying this method in terms of tools and necessary data to develop; moreover, this article shows the application of this method through performance assessment. By applying the introduced method in FEMA P-58 project, it is possible to easily obtain loss curves which lead to better communication and making decision between employers and engineers. Estimating annual loss is considered as a valuable tool to assess insurance value for the buildings. Also, it is possible to say that for the investigated structure, braced structure has less annual loss and the structure with high ductility can be considered as the next option. Costs of performing each of the buildings should be considered by contractors, as the cost of execution a structure with high ductility is more than other buildings; moreover, contractors should regard the ease of windward implementation, its speed and architecture. Thus, there is a better communication between contractors and engineers in comparison with the past which leads to making a better decision about designing a structure.
 Jalayer, F. (2003) Direct Probabilistic Seismic Analysis: Implementation Non-Linear Dynamic Assessments. A Dissertation Submitted to the Department of Civil and Environmental Engineering and the Committee on Graduate Studies of Stanford University in Partial Fulfillment of the Requirements for the Degree of Doctor of Philosophy.
 Tafakori, E., Pourzeynali, S. and Estekanchi, H.E. (2017) Probabilistic Seismic Loss Estimation via Endurance Time Method. Earthquake Engineering & Engineering Vibration, 16, 233-245. https://doi.org/10.1007/s11803-017-0379-8
 Beheshti Aval, S.B., Kouhestani, H. and Mottaghi, L. (2017) Effectiveness of Two Conventional Methods for Seismic Retrofit of Steel and RC Moment Resisting Frames Based on Damage Control Criteria. Earthquake Engineering & Engineering Vibration, 16, 537-555.
 Beheshti Aval, S.B. and Masoumi Verki, A. (2019) Seismic Reliability Assessment of a Steel Moment-Resisting Frame with Two Different Ductility Levels Using a Cloud Analysis Approach. Earthquake Engineering & Engineering Vibration, 18, 171-185.
 Lignos, D. and Krawinkler, H. (2007) A Database in Support of Modeling of Component Deterioration for Collapse Prediction of Steel Frame Structures. Research Frontiers at Structures Congress, Long Beach, 16-19 May 2007, 16-19.
 Masoumi Verki, A. (2012) Performance-Base Design through Implementation of ATC-58 Guideline on a Conventional Structure (Flexural-Bracy Frame) in Iran. Thesis Submitted for the Degree of Master of Seismic Engineering, K. N. Toosi University of Technology (KNTU), Department of Civil Engineering, Tehran.
 Federal Emergency Management Agency (FEMA) (2009) Quantification of Building Seismic Performance Factors. Report No. FEMA-P695, Prepared by the Building Seismic Safety Council for the Federal Emergency Management Agency, Washington DC.
 Khorami, M., Alvansazyazdi, M., Shariati, M., Zandi, Y., Jalali, A. and Tahir, M. (2017) Seismic Performance Evaluation of Buckling Restrained Braced Frames (BRBF) Using Incremental Nonlinear Dynamic Analysis Method (IDA). Earthquakes and Structures, 13, 531-538.
 Li, B., Xie, W. and Pandey, M.D. (2016) Newmark Design Spectra Considering Earthquake Magnitudes and Site Categories. Earthquake Engineering & Engineering Vibration, 15, 519-535. https://doi.org/10.1007/s11803-016-0341-1
 Liang, F., et al. (2017) Accuracy of Three-Dimensional Seismic Ground Response Analysis in Time Domain Using Nonlinear Numerical Simulations. Earthquake Engineering & Engineering Vibration, 16, 487-498.
 Bahramirad, A., Tehranizadeh, M. and Moshref, A. (2015) Equating Incremental Dynamic Analysis with Static Nonlinear Analysis at Near-Field Excitation. Earthquake Engineering & Engineering Vibration, 14, 465-476.
 Mehanny, F. and Deierlein, G. (2000) Modeling of Assessment of Seismic Performance of Composite Frames with Reinforced Concrete Columns and Steel Beams. John A. Blume Earthquake Engineering Research Center, Department of Civil Engineering, Stanford University, Stanford.
 Shome, N. and Cornell, C.A. (1999) Probabilistic Seismic Demand Analysis of Nonlinear Structures. Reliability of Marine Structures, Program Report No. RMS-35, Department of Civil and Environmental Engineering, Stanford University, Stanford.
 Zareian, F., Lignos, D., Krawinkler, H. and Ibarra, L. (2010) Basic Concept and Performance Measures in Prediction of Collapse of Buildings under Earthquake Ground Motions. The Structural Design of Tall and Special Buildings, 19, 167-181.
 Aslani, H. and Miranda, E. (2005) Fragility Assessment of Slab-Column Connections in Existing Non-Ductile Reinforced Concrete Buildings. Journal of Earthquake Engineering, 9, 777-804. https://doi.org/10.1080/13632460509350566
 Beheshti, A., Masoumi Verki, A., Rastegaran, M., et al. (2014) Systematical Approach to Evaluate Collapse Probability of Steel MRF Buildings Based on Engineering Demand and Intensity Measure. International Journal of Structural Analysis & Design, 1, 14-18.
 Federal Emergency Management Agency (FEMA) (2000) State of the Art Report on Performance Prediction and Evaluation of Steel Moment-Frame Buildings. Report No. FEMA-355F, Federal Emergency Management Agency, Washington DC.
 Ibarra, L.F. and Krawinkler, H. (2005) Global Collapse of Frame Buildings under Seismic Excitations. Report No. PEER2005/06, Pacific Earthquake Engineering Research Center, University of California at Berkeley, Berkeley.
 Mahdavi, A., Banazadeh, M. and Deylami, A. (2011) Bayesian Approach for Determination of Drift Hazard Curves for Generic Steel Moment-Resisting Frames in Territory of Tehran. International Journal of Civil Engineering & Structural Engineering, 9, 145-154.
 Mahdavi, A., Banazadeh, M. and Deylami, A. (2010) Determination of Drift Hazard Curves of Steel Moment-Resisting Frames for Territory of Tehran City. 9th US National and 10th Canadian Conference on Earthquake Engineering, Toronto, 25-29 July 2010, 25-29.