Assessment of Seismic Vulnerability of Steel and RC Moment Buildings Using HAZUS and Statistical Methodologies

1Department of Civil Engineering, Birjand University of Technology, Birjand 97175-569, Iran 2Department of Civil and Environmental Engineering, Incheon National University, 12-1 Songdo-dong, Yeonsu-gu, Incheon 406-840, Republic of Korea 3Incheon Disaster Prevention Research Center, Incheon National University, 12-1 Songdo-dong, Yeonsu-gu, Incheon 406-840, Republic of Korea 4Faculty of Engineering, University of Mohaghegh Ardabili, Ardabil 56199-11367, Iran 5TAAT Investment Group, Tehran 18717-13553, Iran


Introduction
Collapse fragility curves give information on how the probability of structural damage increases with increasing peak ground acceleration.Collapse capacities are presumed to follow a lognormal distribution.Literature reviews show that utilization of fragility curves began with nuclear facilities; because the damage to these structures can be so devastating, these structures are considered extremely important.Fragility curves were first designed for nuclear power plants in 1980 [2].Following that, these curves were further developed by Kircher and Martin 1993 [3]; the curves are simple in terms of calculation, and they had been obtained experimentally and with engineering judgment.Moreover, they had considered the vertical and horizontal axes based on emergence probability and the qualitative amount of earth motion, respectively.Recent studies have found the lognormal distribution to provide good results for illustrating the collapse capacity data [4][5][6][7].Several approaches have been suggested for the fragility analysis of structures: empirical, analytical, and hybrid methods [8].Rossetto and Elnashai [9] used the empirical method by using data gathered on structures damaged during past earthquakes.
An investigation on methods to evaluate the seismic collapse of structures by Villaverde [10] and global collapse of frames by Ibarra and Krawinkler [4] as well as the ATC-63 project [11] has reported extensive literature reviews of analytical methods available to evaluate the capacity of structures to withstand seismic collapse.After the Northridge earthquake, more attention was allocated to estimate the amount of structural damage.During a study on California structures, principles of ATC-13 were used to plot the fragility curves in 1994.In that work, fragility curves were plotted for steel moment resisting frames and wooden and reinforced concrete structures.Anagnos et al. [12] undertook further studies based on inserted load distribution in ATC in 1995, and they presented a new model of fragility curves.In this paper, all seismic calculations were performed based on ATC-13, where the probability distribution function had been as assumed to be the normal distribution.Fragility curves were produced for Memphis buildings in 1996, and the displacement among stories was investigated at different levels.Singhal and Kiremidjian estimated fragility curves according to the observations in the one-story buildings.They used the Park-Ang index to assess damage and then expressed the amount by statistical distribution function undergoing various earthquakes [13].Shinozuka [14] presented an investigative paper about the fragility curve of a one-span bridge that, in comparison with former works, was very precise.Because a nonlinear dynamic technique and statistical data had been used, they presented a method for obtaining fragility curves in bridges.Later, Saxena et al. [15] used this procedure for multispan bridges and reinforced concrete, in which the generation of fragility curves for bridges was presented step by step for the first time, and the probability distribution of the data was assumed to be the normal distribution.Tanaka et al. [16] used the lognormal distribution to calibrate fragility curves in 2000.He classified 3683 bridges into 5 groups, defined the amount of damage for 5 levels, and assessed lognormal distribution parameters.
Researchers have also plotted fragility curves for steel, wooden, and concrete structures by using the damage functions and PGV in the Japan Kobe in 2000 [17].They came to the conclusion that the reinforced concrete and steel structures have the minimum and maximum fragility, respectively, in this city.Furthermore, the production procedure of fragility curves generated using the capacity spectrum method was presented in 2001.In this work, researchers   plotted fragility curves for displacement among stories by utilizing the nonlinear static analysis method.Pagni [18] presented a damage model for old reinforced concrete elements in Washington University.He introduced 12 states of damage, which included primary cracking and the common face of beam-column.Moreover, damage states were divided into 2 groups of cracking and crumbling, and the best distribution was selected.Fragility curves were generated for steel moment resisting frame buildings by using Perform software and FEMA codes in 2006.Nonlinear dynamic analysis of 2-, 3-, 4-, 6-, 8-, and 10-story frames was performed.Ozer et al. [19] dealt with developing seismic failure curves for 3-, 5-, 7-, and 9-story buildings that were designed based on the TSC-1998 prevalent code.The structure model was the concrete moment resisting frame which was designed twodimensionally by using Sap 2000 and IDARC-2D software, which were used for analysis and evaluation.Structure frames were classified into poor (with low quality) and superior (with high quality) groups based on the special features, construction process, and observed performance after large earthquakes in Turkey.Failure curves were plotted in different groups according to the parameters affecting the seismic performance, and the number of stories was considered as an important parameter in a building's seismic vulnerability.Barkhordary and Tariverdilo [20] studied an evaluation of the effect of column bar patching and bar slipping on the structure's vulnerability, utilizing concrete moment resisting frame systems in Iran.The structure model was considered two-dimensionally and was analyzed using OpenSees (2006).
The seismic reliabilities of two typical stainless steel legged wine storage tanks have been studied by Colombo and Almazán [21].Seismic fragility of nonlinear connections in large piping systems has been evaluated by Tadinada and Gupta [22].Seismic fragility of RC structures with shape memory alloys has been investigated by Mirtaheri et al. [23].
Finally, the results of additive dynamic analysis were used for generating the fragility curves.Nazri and Saruddin [24] analyzed steel and concrete moment resisting frames under near and far field earthquakes using IDA, and, according to the fragility curve, they presented the possibility of reaching or exceeding damage limits for each structure at 5 levels of performance introduced in FEMA-273.HAZUS (Hazard US) is a well-known loss estimation methodology that defines five damage states, none, slight, moderate, extensive, and complete, using physical   (qualitative) descriptions of damage to building elements.The HAZUS methodology has been applied to various seismic hazard assessment studies by adapting capacity and fragility curves for structures in specific regions [25,26].
In areas with significant seismic activity, choosing the kind of structure has always been a serious challenge for designer engineers.In this regard, development of fragility curves can help designer select a structure that has the least fragility.On the other hand, the question presented is that, among all common methods, which one is more accurate and comprehensive?In this paper, we aim to provide an answer to this and other important questions.

Damage States
To investigate the vulnerability of structures, it is necessary to consider a damage index.In this study, interstory drift is considered as the damage index.Additionally, a series of failure modes is defined and evaluated in terms of seismic vulnerability.In this paper, the limit states used are based on the HAZUS-MH MR5 guidelines [1].According to HAZUS, there are 36 building types, and each type is divided into three categories in terms of number of stories: low-rise, mid-rise, and high-rise, which generally correspond to 1-3, 4-7, and 8+ stories, respectively.The limit and intensity of damage to structural components of a building are described by one of four damage states: slight, moderate, extensive, and complete [1].The values for each damage state are presented in Table 1.Damage occurs when the drift value reaches the limit state.

Building Design and Modeling
According to HAZUS guideline, 4-to 7-story buildings are classified as mid-rise buildings and the buildings with higher 8 stories are classified as tall buildings.Therefore, in this study 5-, 8-, and 12-story buildings are analyzed as a sample of low-rise, border of mid-and high-rise, and high-rise buildings, respectively.The 5-, 8-, and 12-story steel and concrete buildings have the same plan view (Figure 1) and seismic design category; they are designed based on the AISC 360-10, ACI 318, and ASCE 7-10 codes.The height of the first story is 2.7 m (to top of beam), and the height of all the other stories is 3.2 m.The dead and live loads are 460 and 200 kg/m 2 uniformly distributed over each floor.The specimens are analyzed using Perform 3D v4.0.3, which is a nonlinear dynamic analysis program.Beams, columns, and connection regions are modeled individually, and ground acceleration records were applied to the specimen models.All of the structural members for all the buildings can be found in Figures 21 and 22 and Tables 5, 6, 7, 8, 9, and 10.

Ground Motion Database
According to [27], accurate assessment of damage demand can be achieved using 10 to 20 ground motion records.In this study, ten ground motion records were selected for the nonlinear time history analyses.The selected records consist of site classes C and D extracted from the PEER database and scaled from 0.1 g to 1.5 g by steps of 0.1 g.Pertinent information on the ground motion database is presented in Table 2.

Fragility Analysis
To facilitate the assessment of the seismic vulnerability of the structures, fragility functions defining the probability of meeting or exceeding a specific limit state given an earthquake intensity level were developed using the nonlinear time history analyses described above.
The Pacific Earthquake Engineering Center gives a formulation for the fragility analysis and probabilistic seismic assessment framework.This framework proposes a relationship between a certain key decision variable, DV (e.g., annual seismic loss), and a series of conditional probabilities, deaggregating average annual frequency (AAF) of exceeding the decision parameter,  DV , in terms of damage measure, DM, and ground motion intensity measure, IM, as follows: in which (DV | DM) denotes the probability that the DV exceeds specified values given that the engineering DMs are equal to special values, (DM | IM) is the probability that the DMs exceed these values given that the IM is equal special value, and  IM is the AAF of the ground motion IM.The part (DM | IM) in the deaggregated equation is named the seismic fragility function.Seismic fragility relates the probability of achieving or exceeding predefined levels of damage to the severity of ground motion intensity.Analytical fragility curves are derived from the results of numerical models of the buildings subjected to ground motion records.The fragility function is given in where Φ(⋅) stands for the standard normal cumulative distribution function,  is the superior bound for each damage level,  is the median value of demand as the function of IM,  |IM is the logarithmic standard deviation of the demand conditioned on the IM,   is the capacity uncertainty, and   is the modeling uncertainty.
The output of the linear regression analysis, in a logtransformed space and with empirically based limit-state models, may be used to develop closed-form expressions for fragility functions, assuming limit-state models also follow a lognormal distribution: in which the left-hand side is the probability of exceeding the damage limit,  is damage, B is the standard normal cumulative distribution,   is the median estimate of demand (maximum story drift),   is the median estimate of capacity, and   = √ 2 |IM +  2  , where  |IM is the dispersion of the demand conditioned on the intensity measure and   is the dispersion of the capacity.

HAZUS Methodology.
The values of   and   for various building types were taken from the HAZUS guidelines (Table 3).
The parameter   can be calculated as follows: where  and  are regression coefficients which can be obtained using regression analysis of incurred damage versus PGA [1].For example, Figure 2 shows the regression of the 5-story RC frame.After defining and modeling the nonlinear properties of the elements, the scaled ground motion records are applied to structures using Perform 3D software.For each model, the maximum interstory drifts are derived and compared with HAZUS criteria to determine the number meeting or exceeding the limit states.For example, the 5-story RC frame is analyzed 10 times at an acceleration of 0.7 g for 10 ground motion records; the results (Table 4) show that, in 10, 8, 3, and 0 cases, the structure has experienced slight, moderate, extensive, and complete damage states, respectively.
The fragility curves for both steel and RC frames are plotted as shown in Figures 3-8.The curves include four damage states based on HAZUS guidelines.
The curves shifting left mean that the probability of damage has increased.All curves show that, for slight and moderate states, the slope of curve is initially steep and then decreases.Usually, this slope change occurs at PGA > 0.8 g, because until 0.8 g the probability of exceeding the limit states varies between 0 and 100% while after 0.8 g this probability closely approaches 100%.For the extensive and complete damage states, the slope of the curve is nearly constant.The results indicate that in the steel frames the highest probability of extensive damage is 60%, while this value is 80% for RC frames.The probabilities of complete damage for steel and RC frames are 10-15% and 20-30%, respectively.
The fragilities of the steel frames for slight, moderate, extensive, and complete states are compared in Figures 9-12.
Results reveal that, in general, the increase in the number of stories leads to increase in the fragility of both steel and RC structures.This difference is more significant in RC frames.For example, in the moderate damage state, the difference between 5 and 12 steel frames is 15%, while this value for 5 and 12 RC frames is approximately 30%.These comparisons can also be seen in Figures 13-16.
The fragilities of the 12-story steel and RC frames are compared for slight, moderate, extensive, and complete damage states as shown in Figures 3-6.
Figures 3-6 demonstrate that the fragility of a 12-story RC frame is higher compared to that of a 12-story steel frame.This difference is more significant in higher damage states, that is, at the slight damage state there is minimal difference, and at the complete state, there is maximal difference.These figures also show that the curves of the slight and moderate damage states begin to converge with increasing PGA, while the curves of the extensive and complete states begin to diverge with increasing the PGA.
The fragilities of the 8-story steel and RC frames are compared for slight, moderate, extensive, and complete damage states, as shown in Figures 7-10.
Figures 7-10 show that, similar to 12-story models, 8-story RC frames have higher fragility compared with steel frames.These figures also show that the curves of the slight and moderate damage states begin to converge with increasing PGA, while the curves of the extensive and complete states begin to diverge with increasing the PGA.
The fragilities of the 5-story steel and RC frames are compared for the slight, moderate, extensive, and complete states, as shown in Figures 11-14.
Evaluation of  indicates similar results for 12and 8-story frames.Furthermore, in the 5-story frames at the complete damage state, the fragility curves nearly overlap at low PGA and only diverge slightly at high PGA.The most significant difference is seen at the moderate damage state, and the least significant is at the complete damage state.

Statistical Methodology.
As described previously, the fragility curve can be obtained using this equation: In Section 5.1, the required parameters and coefficients have been defined using HAZUS.In this part, these parameters are obtained using nonlinear dynamic analyses and a statistical postprocessing.The ground motion records are scaled from 0.1 g to 1.5 g by a step of 0.1 g. Figure 15 shows the fragility curves for a 5-story steel frame using both HAZUS and the statistical procedure.The figure shows that, generally, the statistical method yields higher fragility.The most significant difference is seen at the extensive damage state and the least significant at the complete damage state.The curves for the slight damage state nearly overlap, and at PGA > 0.7 g the probability approaches 100%.At the moderate damage state, the difference between the statistical and HAZUS methodologies is approximately 10%.
The comparison between statistical and HAZUS methodologies for an 8-story steel frame is shown in Figure 16.
The comparison between statistical and HAZUS methodologies for a 12-story steel frame is shown in Figure 17.
As is evident in Figures 15 to 18, the HAZUS method yields a reduced probability of damage.The difference between the curves is higher compared to the 5-story steel frame.Similar to the 5-story steel frame, in the extensive damage state, the most significant difference (30%) is observed.
Figure 17 shows the fragility curves for a 12-story steel frame.The most significant difference is between the moderate and extensive damage states (20%).This difference is 10% for the slight damage state, and in the complete state the curves are approximately the same.
The fragility curves for RC frames are presented in Figure 18.
Figure 18 shows that the HAZUS procedure yields a reduced probability of damage.At three damage states, slight, moderate, and extensive, the difference between the HAZUS and statistical procedures is 25%, but at the complete state, the curves nearly overlap at low PGA and diverge only slightly at large PGA.A comparison between the 5-story steel and RC frames shows that this difference percentage is higher in the RC case. Figure 19 shows that similar conclusions can be drawn for an 8-story RC frame.
Seismic vulnerability of a 12-story RC frame has been assessed using two approaches (Figure 20).The comparison  shows that, similar to previous conclusions, the HAZUS procedure yields reduced damage probability except at PGA = 0.6 g for the extensive damage state.The most significant difference can be seen in the extensive damage state (statistical method yields 20% higher fragility).In other states, the curves are fairly close together.
Generally, in RC frames, the difference between the statistical and HAZUS procedures is 20%.

Conclusions
This paper investigated the seismic vulnerability of steel and RC frames of various heights by the use of the seismic fragility analysis.Two methodologies, HAZUS and a statistical method, have been employed, and the results were compared.The main conclusions are as follows: (i) Generally, an increase in PGA leads to an increased fragility.(ii) The slope of the fragility curve is larger at lower PGA and smaller at higher PGA.(iii) For both steel and RC frames, the increase in number of stories leads to higher fragility (especially in RC frames).(iv) For RC frames at slight and moderate damage states, the fragility curves of 5-and 8-story frames are very close together, while the 12-story curve shows a significant difference.(v) In steel frames, the difference in fragility for 8-and 12-story frames is not significant (all 4 damage states are similar), while the 5-story frame shows difference from others.(vi) Generally, the fragility of RC frames is higher compared to steel frames.(vii) Comparison between the HAZUS and statistical methodologies shows that the fragility curves obtained from the statistical procedure yield higher fragility; the most difference is approximately 30%.

Figure 3 :Figure 4 :Figure 5 :Figure 6 :
Figure 3: Comparison of the fragility of the 12-story steel and RC frame buildings for the slight damage state.

Figure 7 :
Figure 7: Comparison of the fragility of the 8-story steel and RC frame buildings for the slight damage state.

Figure 8 :Figure 9 :Figure 10 :Figure 11 :Figure 12 :Figure 13 :Figure 14 :Figure 15 :
Figure 8: Comparison of the fragility of the 8-story steel and RC frame buildings for the moderate damage state.

Figure 17 :
Figure 17: The fragility curves for the 12-story steel frame buildings.

Figure 18 :
Figure 18: The fragility curves for the 5-story RC frame buildings.

Figure 19 :
Figure 19: The fragility curves for the 8-story RC frame buildings.

Figure 22 :
Figure 22: Typical plan view of steel columns.

Table 4 :
Number of instances reaching or exceeding the limit states for the 5-story RC frame.

Table 5 :
Sections for the 5-story RC frame.

Table 6 :
Sections for the 8-story RC frame.

Table 7 :
Sections for the 12-story RC frame.

Table 8 :
Sections for the 5-story steel frame.

Table 9 :
Sections for the 8-story steel frame.

Table 10 :
Sections for the 12-story steel frame.