Effect of Wind Loads on Collapse Performance and Seismic Loss for Steel Ordinary Moment Frames

: The aim of this study is to investigate the effect of wind loads on the seismic collapse performance and seismic loss for steel ordinary moment frames (OMFs). For this purpose, 9-, 12-, 15-, and 18-story steel OMFs are repeatedly designed for (1) gravity load + seismic load, (2) gravity load + seismic load + wind load (wind speed = 44 m/s), and (3) gravity load + seismic load + wind load (wind speed = 55 m/s). The seismic collapse performance and seismic loss of OMFs are evaluated using the procedures in FEMA P695 (FEMA, 2009) and FEMA P58 (FEMA, 2018), respectively. Steel OMFs designed with consideration of wind loads have larger member sections than corresponding steel OMFs designed without consideration of wind loads as expected. Although member sections are increased when wind loads are considered, the growth in the maximum base shear force and lateral stiffness of OMFs are insigniﬁcant. Unlike our expectation, OMFs designed with consideration of wind loads have higher expected annual loss ( EAL ) than corresponding OMFs designed without consideration of wind loads.

The most stringent design and detail requirements are required on SMFs among the three types of moment frames, whereas the least requirements are specified for OMFs [5,6]. According to ANSI/AISC 341-16 [3], connections in SMFs and IMFs should have a rotation capacity of at least 0.04 and 0.02, respectively, whereas OMF connections should have minimal inelastic rotation capacity. The column-to-beam moment strength ratio is not specified for OMFs. Panel zones in OMFs should resist design force calculated from elastic analyses rather than maximum forces transferred from adjoining beams and columns. To compensate for the relaxed design and detail requirements on OMFs, OMFs should be designed with a larger design shear force than IMFs and SMFs.
Many studies have been conducted to assess the seismic performance of steel SMFs [7][8][9][10][11][12][13][14][15][16][17][18]. To obtain the reliable results of seismic performance assessment for steel SMFs, the cyclic behavior of SMF components has been investigated through experimental and numerical studies [19,20]. However, only a limited number of studies have been conducted for steel IMFs, OMFs, and existing moment frames [21][22][23][24][25][26]. According to the commentary of AISC 341-16 [3], the details of OMFs were developed based solely on the judgement of engineers rather than on the results of experimental and numerical studies. Thus, it is important to evaluate the seismic performance of OMFs designed according to current seismic design codes, which was the motivation of this study.
Steel OMFs are permitted only for seismic design categories (SDC) A, B, and C [1]. In such SDCs, wind loads govern design results for OMFs rather than seismic loads, particularly for high-rise buildings [27,28]. However, most studies on steel moment frames have focused on steel SMFs located in high-seismicity regions. In addition, wind loads have often been ignored in seismic design procedures. Thus, it is important to evaluate the seismic performance of OMFs designed according to current seismic design codes, which was the motivation of this study.
Since seismic and wind loads mainly generate lateral force and drift demands on building structures, it is important to assess the effect of wind loads on the seismic performance of the OMFs. Although larger member sections should be used for OMFs when wind loads are considered, the use of the larger member sections may not guarantee better seismic performance of the OMFs. To the best of our knowledge, no study has been conducted to investigate the effect of wind loads on the seismic performance of the OMFs. Assessment of wind load effect was another motivation of this study. The aim of this study is to evaluate the seismic performance of steel OMFs with consideration of wind load effects. The seismic performance was estimated in terms of seismic collapse performance and seismic loss.

Research Methodology
To assess the seismic performance of code-compliant steel OMFs and the effect of wind loads on the seismic performance, collapse performance was estimated according to the methodology provided in FEMA P695 [29]. In ASCE 7 [1], target seismic performance objectives are defined as seismic collapse probabilities.
For this purpose, steel OMFs with four different building heights were designed according to ASCE 7 and reference standards with and without consideration of wind loads. The seismic collapse performance of these OMFs were estimated according to FEMA P695. In addition, seismic loss estimation was conducted for these OMF buildings according to FEMA P58 [30] to investigate the effect of wind loads. The methodologies provided in FEMA P58 had been developed for performance-based earthquake seismic engineering. Figure 1 summarizes the methodology used in this study. Appl. Sci. 2022, 12, x FOR PEER REVIEW 2 of 19 of AISC 341-16 [3], the details of OMFs were developed based solely on the judgement of engineers rather than on the results of experimental and numerical studies. Thus, it is important to evaluate the seismic performance of OMFs designed according to current seismic design codes, which was the motivation of this study. Steel OMFs are permitted only for seismic design categories (SDC) A, B, and C [1]. In such SDCs, wind loads govern design results for OMFs rather than seismic loads, particularly for high-rise buildings [27,28]. However, most studies on steel moment frames have focused on steel SMFs located in high-seismicity regions. In addition, wind loads have often been ignored in seismic design procedures. Thus, it is important to evaluate the seismic performance of OMFs designed according to current seismic design codes, which was the motivation of this study.
Since seismic and wind loads mainly generate lateral force and drift demands on building structures, it is important to assess the effect of wind loads on the seismic performance of the OMFs. Although larger member sections should be used for OMFs when wind loads are considered, the use of the larger member sections may not guarantee better seismic performance of the OMFs. To the best of our knowledge, no study has been conducted to investigate the effect of wind loads on the seismic performance of the OMFs. Assessment of wind load effect was another motivation of this study. The aim of this study is to evaluate the seismic performance of steel OMFs with consideration of wind load effects. The seismic performance was estimated in terms of seismic collapse performance and seismic loss.

Research Methodology
To assess the seismic performance of code-compliant steel OMFs and the effect of wind loads on the seismic performance, collapse performance was estimated according to the methodology provided in FEMA P695 [29]. In ASCE 7 [1], target seismic performance objectives are defined as seismic collapse probabilities.
For this purpose, steel OMFs with four different building heights were designed according to ASCE 7 and reference standards with and without consideration of wind loads. The seismic collapse performance of these OMFs were estimated according to FEMA P695. In addition, seismic loss estimation was conducted for these OMF buildings according to FEMA P58 [30] to investigate the effect of wind loads. The methodologies provided in FEMA P58 had been developed for performance-based earthquake seismic engineering. Figure 1 summarizes the methodology used in this study.

Design of Model Steel Ordinary Moment Frames with Consideration of the Level of Wind Loads
In this study, to evaluate the seismic performance of steel OMFs, model office buildings with different structural heights were designed: 9-, 12-, 15-, and 18-story OMFs. Moment frames are economically effective for buildings up to 25 stories, above which their resistance is costly to control. In addition, member sections in OMFs with more than 18 stories could not be chosen from the section tables provided in the AISC steel construction manual [31]. Thus, this study considered OMFs with 18 or less stories.
The W-shaped sections provided in the AISC steel construction manual [31] are used for the OMFs. The plan and elevation of the 15-story model OMF buildings are shown in Figure 2. OMFs are denoted with thick solid lines, which are placed along the perimeter of the building (Figure 2a). The span length is 6.1 m and the height of the first story is 4.6 m, whereas the story height for other stories is 4.0 m. It is assumed that the column splices are located in the middle of the column for every two stories except the first story, as shown in Figure 2b-d. ETABS software [32] is used for the analyses and design.

Design of Model Steel Ordinary Moment Frames with Consideration of the Level of Wind Loads
In this study, to evaluate the seismic performance of steel OMFs, model office buildings with different structural heights were designed: 9-, 12-, 15-, and 18-story OMFs. Moment frames are economically effective for buildings up to 25 stories, above which their resistance is costly to control. In addition, member sections in OMFs with more than 18 stories could not be chosen from the section tables provided in the AISC steel construction manual [31]. Thus, this study considered OMFs with 18 or less stories.
The W-shaped sections provided in the AISC steel construction manual [31] are used for the OMFs. The plan and elevation of the 15-story model OMF buildings are shown in Figure 2. OMFs are denoted with thick solid lines, which are placed along the perimeter of the building (Figure 2a). The span length is 6.1 m and the height of the first story is 4.6 m, whereas the story height for other stories is 4.0 m. It is assumed that the column splices are located in the middle of the column for every two stories except the first story, as shown in Figure 2b-d. ETABS software [32] is used for the analyses and design. The OMFs were designed according to ASCE/SEI 7-16 [1], ANSI/AISC 360-16 [2], and ANSI/AISC 341-16 [3]. The response modification coefficient (R), overstrength factor (Ω 0 ), and deflection modification factor (C d ) for the steel OMFs are 3.5, 3.0, and 3.0, respectively, according to ASCE/SEI 7-16 [1]. The risk category and importance factor (I e ) are assumed as II and 1.0 according to ASCE 7-16 [1], respectively. Design earthquake spectral response acceleration parameters at short period (S DS ) and 1-s period (S D1 ) are defined as 0.499 and 0.199, respectively. The seismic design category (SDC) of the buildings is classified as C max .
The dead, cladding, and live loads of typical office buildings are 4.31 kN/m 2 , 1.20 kN/m 2 , and 2.40 kN/m 2 , respectively [12]. The steel material used in the design is ASTM A992 Gr.50 steel, which has a specified minimum yield strength (F y ) of 345 MPa and a specified tensile strength (F u ) of 450 MPa. The seismic loads were determined using the equivalent lateral force (ELF) procedure.
Wind loads were determined using basic wind speed with a mean recurrence interval (MRI) of 700 years, according to Figure 26.5-1B of ASCE 7-16. The exposure category of B was used for the OMF buildings because the model buildings were assumed to be located in urban and suburban flat areas. As seen in Figure 26.5-1B of ASCE 7-16, the basic wind speed for most areas in the USA range from 45 m/s to 51 m/s. In this study, the upper (=51 m/s) and lower bound values (=44 m/s) for the range were considered to assess the effect of wind loads on the seismic performance of OMFs.
The exposure category is defined as B, assuming that the model buildings are located in the urban or suburban flat area. Table 1 summarizes the fundamental periods (T 1 ) of the model OMFs. As expected, T 1 tends to increase as the height of OMFs increases. Twelve-, fifteen-, and eighteen-story steel OMFs, designed with consideration of wind loads, have larger member sections than corresponding OMFs designed without consideration of wind loads. The increase in member sections due to wind loads is prominent in lower stories. However, for nine-story OMFs, the influence of wind loads on design results is negligible when the basic design wind speed is 44 m/s. Table 1. Calculated fundamental period (T 1 ) and upper limit of period (C u T a ) according to ASCE 7-16 of the designed steel ordinary moment frames.

Numerical Models of Steel Ordinary Moment Frames
To evaluate the seismic performance of model OMFs, a numerical model is constructed using OpenSees software [33]. Figure 3 presents a nonlinear numerical model used for steel OMFs. The numerical model consists of columns, beams, and panel zones, and a leaning column installed to simulate the P − ∆ effect induced by gravity loads. The leaning column is modeled as a rigid truss element and connected to the ground by a pin. The leaning column is connected to the OMF using pins (Figure 3a). The gravity axial load applied on the leaning column is calculated with effective seismic weight (1.05 DL + 0.25 LL) [29], where DL and LL are the superimposed dead and live loads, respectively. The damping matrix is constructed using the commonly used Rayleigh damping matrix based on a damping ratio of 2% for the first and fifth modes [34][35][36]. It is noted that in NIST GCR 17-917-46v1 [34] and ATC 72-1 [35], it is recommended that the first-mode period and the third or higher period should be used to construct the Rayleigh damping matrix. In this study, the first and fifth modes were used. In order to avoid unrealistic damping forces, stiffness-and mass-proportional damping terms are constructed only using the elastic line elements for beams and columns, and nodes for lumped mass [36]. and a leaning column installed to simulate the − ∆ effect induced by gravity loads. The leaning column is modeled as a rigid truss element and connected to the ground by a pin. The leaning column is connected to the OMF using pins (Figure 3a). The gravity axial load applied on the leaning column is calculated with effective seismic weight (1.05 + 0.25 ) [29], where and are the superimposed dead and live loads, respectively. The damping matrix is constructed using the commonly used Rayleigh damping matrix based on a damping ratio of 2% for the first and fifth modes [34][35][36]. It is noted that in NIST GCR 17-917-46v1 [34] and ATC 72-1 [35], it is recommended that the first-mode period and the third or higher period should be used to construct the Rayleigh damping matrix. In this study, the first and fifth modes were used. In order to avoid unrealistic damping forces, stiffness-and mass-proportional damping terms are constructed only using the elastic line elements for beams and columns, and nodes for lumped mass [36]. Individual beam and column elements are modeled using an elastic line element with rotational plastic hinge elements placed at both ends of the elastic element. The moment capacity of the column section should be reduced due to stress generated by axial force. The modeling parameters of the concentrated plastic hinges for steel columns are determined according to Lignos et al. [37].
The welded unreinforced flange-bolded web (WUF-B) connection is used for OMFs in this study, which is modeled as shown in Figure 3b. A panel-zone consists of four rigid line elements and one rotational spring element simulating the trilinear behavior of the panel zone [38]. To simulate the cyclic behavior of the WUF-B connection, the modified IMK model [39] is used. As shown in Figure 3b, is the drift ratio ( ) at capping strength, is the value at fracture, is the capping moment strength, and is the effective yield strength. The and values for beams are estimated according to Lignos and Krawinkler [19], and the and values are estimated according to ASCE 41 [40]. The residual strength is defined as 20% of the plastic moment strength for WUF-B connections according to Stojadinovic et al. [41] and Han et al. [42]. The accuracy of the connection model was verified by Han et al. [43]. In Figure 3b, empty circles with black, blue, and red colors denote yielding, capping, and fracture, respectively. Panel zones are modeled according to Gupta and Krawinkler [38]. Composite actions [11,[44][45][46]  Individual beam and column elements are modeled using an elastic line element with rotational plastic hinge elements placed at both ends of the elastic element. The moment capacity of the column section should be reduced due to stress generated by axial force. The modeling parameters of the concentrated plastic hinges for steel columns are determined according to Lignos et al. [37].

due to
The welded unreinforced flange-bolded web (WUF-B) connection is used for OMFs in this study, which is modeled as shown in Figure 3b. A panel-zone consists of four rigid line elements and one rotational spring element simulating the trilinear behavior of the panel zone [38]. To simulate the cyclic behavior of the WUF-B connection, the modified IMK model [39] is used. As shown in Figure 3b, θ p is the drift ratio (θ) at capping strength, θ g is the θ value at fracture, M c is the capping moment strength, and M y is the effective yield strength. The M c and M y values for beams are estimated according to Lignos and Krawinkler [19], and the θ p and θ g values are estimated according to ASCE 41 [40]. The residual strength is defined as 20% of the plastic moment strength for WUF-B connections according to Stojadinovic et al. [41] and Han et al. [42]. The accuracy of the connection model was verified by Han et al. [43]. In Figure 3b, empty circles with black, blue, and red colors denote yielding, capping, and fracture, respectively. Panel zones are modeled according to Gupta and Krawinkler [38]. Composite actions [11,[44][45][46] due to slab effects are not considered in this study because of insufficient test data to verify the numerical models.

Seismic Collapse Performance of Steel Ordinary Moment Frames
Nonlinear static analyses (pushover) are performed using lateral forces with a vertical distribution following the fundamental mode shape of a frame according to FEMA P695 [29] methodology. Figure 4 shows the pushover curves of OMFs, from which overstrength factor (Ω) and period-based ductility (µ T ) are determined using the following equations: where V max is the maximum base shear, V design is the design base shear force, δ y,e f f is the effective yield roof displacement, and δ u is the maximum roof displacement measured when the strength decreases to 0.8 V max .
where is the maximum base shear, is the design base shear force, , is the effective yield roof displacement, and is the maximum roof displacement measured when the strength decreases to 0. 8 . Figure 4a,c shows the pushover curves of 9-story and 15-story OMFs. The cyclic curves of OMFs does not vary significantly according to the application of wind loads. The maximum strength and lateral stiffness of W44-and W51-OMFs are similar with those of corresponding W0-OMFs. Figures 4b and 3d show floor displacement normalized by building height ( / ). As mentioned in the previous section, the consideration of wind load in design leads to an increase in member sections, especially for lower stories. Since in upper stories, member sections do not increase unlike in the lower stories, significant damage could occur in the upper stories during earthquakes. This phenomenon is more prominent when OMFs approach collapse state (see Figure 4b,d).   Figures 3d and 4b show floor displacement normalized by building height (δ i /H). As mentioned in the previous section, the consideration of wind load in design leads to an increase in member sections, especially for lower stories. Since in upper stories, member sections do not increase unlike in the lower stories, significant damage could occur in the upper stories during earthquakes. This phenomenon is more prominent when OMFs approach collapse state (see Figure 4b,d). Figure 5 shows the values of Ω and µ T for individual OMFs. As shown in this figure, the Ω value of OMFs does not vary according to structural height and level of wind loads. All steel OMFs have Ω values of about 2.0, which is less than the value of system overstrength factor (Ω 0 = 3.0) specified for steel OMF in ASCE 7 [1]. The µ T value of OMFs tends to decrease gradually with increasing building height. It is also observed that µ T increases with an increase in the level of wind loads. Figure 5 shows the values of Ω and for individual OMFs. As shown in this figure, the Ω value of OMFs does not vary according to structural height and level of wind loads. All steel OMFs have Ω values of about 2.0, which is less than the value of system overstrength factor (Ω = 3.0) specified for steel OMF in ASCE 7 [1]. The value of OMFs tends to decrease gradually with increasing building height. It is also observed that increases with an increase in the level of wind loads.  Figure 6 shows the distribution of damage in the 15-story OMFs designed without and with wind loads (44 m/s and 51 m/s) near collapse. As shown in this figure, damage is concentrated in upper stories for the W44-and W51-OMFs, whereas for the W0-OMF, damage is more widely distributed in the lower stories than that for the W44-, and W51-OMFs. This could be the reason that the maximum shear force of W0-OMF is similar to those of the corresponding W44-, and W51-OMFs (Figure 4a,c), although the member sections for the W0-OMF are less than those of the W44-, and W51-OMFs.  All 15-story OMFs (W0-, W44-, and W51-OMFs) experience yielding and fracture mainly in exterior beams and interior panel zones. The fracture in the interior panel zones can be attributed to the fact that OMF panel zones are designed for the force demand from the load combinations specified in ASCE 7 [1] instead of that developed from the plastic moments of nearby beams. The locations of the most significant damage in the 15-story W0, W44-and W51-OMFs are the fourth, ninth, and twelve stories, respectively, which are the same stories that experience the largest story drift (Figure 4b,d). It is noted that in All 15-story OMFs (W0-, W44-, and W51-OMFs) experience yielding and fracture mainly in exterior beams and interior panel zones. The fracture in the interior panel zones can be attributed to the fact that OMF panel zones are designed for the force demand from the load combinations specified in ASCE 7 [1] instead of that developed from the plastic moments of nearby beams. The locations of the most significant damage in the 15-story W0, W44-and W51-OMFs are the fourth, ninth, and twelve stories, respectively, which are the same stories that experience the largest story drift (Figure 4b,d). It is noted that in all the OMFs in Figure 6, plastic hinges are not developed at the bottom of the first story columns. This indicates that steel OMFs, according to the current design procedure, may not have the maximum base shear capacity intended in seismic design standards. This is attributed to the fact that the Ω values of all OMFs are smaller than overstrength factor (Ω 0 = 3.0) specified for OMFs in ASCE 7 [1] (Figure 5a). For other OMFs with different numbers of stories, similar observations are obtained.
An incremental dynamic analysis (IDA [42]) was also performed for the steel OMFs, from which collapse intensity and fragility curves were obtained. The fragility curve of the model frame was constructed according to the FEMA P695 [29] methodology. An ensemble of 44 far-field ground motions provided by FEMA P695 are used to perform the IDA. The hunt and fill algorithm [47] is used to efficiently conduct the IDA. Due to the limitations of the page, IDA curves for individual ground motions of all frames are not included, and the median IDA curves of each steel OMF are presented in Figure 7. Table 2 summarizes the parameter values obtained from the pushover analysis and the IDA. The collapse intensity (S CT ) of a structure is defined as a 5%-damped pseudo spectral acceleration PSA(T 1 , 5%) at T 1 , which leads to the global dynamic instability (GI) of a structure. When a structure reaches the GI state, the lateral displacement of the structure increases without bound by a slight increase in ground motion intensity.    The median collapse intensity (Ŝ CT ) of W0-OMFs is generally similar with that of corresponding W44-and W51-OMFs, whereŜ CT is the median value of S CT for 44 ground motions (Figure 7). This indicates that the consideration of wind loads in design may not affect the seismic collapse intensity of OMFs. However, the record-to-record dispersion (β RTR ) of W44-and W51-OMFs is larger than that of W0-OMF. With increasing β RTR , the collapse probability of a structure increases. This indicates that the collapse probability of W44-and W51-OMFs is higher than that of corresponding W0-OMF.
According to Eads et al. [48], an increase in β RTR leads to an increase in the collapse probability for low-intensity ground motions, and the effect of β RTR on life-cycle collapse probability is significant. Additionally, with a decrease in T 1 , ground motion intensity demands at a given hazard level increases. The T 1 values of OMFs designed with consideration of wind loads are shorter than those of corresponding OMFs designed without consideration of wind loads. Thus, the collapse performance of W44-and W51-OMFs could be less satisfactory than that of corresponding W0-OMF due to higher seismic intensity demands, similar seismic collapse intensity (capacity), and larger β RTR . Thus, an increase in member sections only in the lower stories of the structure due to the consideration of wind loads may adversely affect the seismic performance of OMFs. This indicates that when evaluating the seismic performance of OMFs, the effect of wind loads should be considered in the evaluation.
In addition, this study inspects the distribution of elastic story drift ratios along the building height for 44 far-field ground motion sets [29]. The story drift ratio (SDR i ) is calculated as the story drift of the ith story divided by story height. Figure 8 shows the distribution of story drift ratios for 15-story W0-and W51-OMFs. A solid gray circle in the figure represents story drift ratios generated by individual ground motions. Median SDR i values for W0-and W51-OMFs are included in this figure.
For the 15-story steel OMFs, the larger elastic SDR i values are obtained in the upper stories of the structure. This is because the effect of higher modes. As seen in the previous section, the W51-OMF has larger SDR i in the upper stories than W0-OMF because of the increase in member sections in the lower stories when wind loads are taken into consideration. The dispersion of SDR i in the upper stories for the W51-OMF is larger than that for the W0-OMF. The larger dispersion of SDR i in the upper stories for the W51-OMF may be attributed to the effect of larger lateral drifts in the upper stories and more significant contributions of higher modes. be considered in the evaluation.
In addition, this study inspects the distribution of elastic story drift ratios along the building height for 44 far-field ground motion sets [29]. The story drift ratio ( ) is calculated as the story drift of the th story divided by story height. Figure 8 shows the distribution of story drift ratios for 15-story W0-and W51-OMFs. A solid gray circle in the figure represents story drift ratios generated by individual ground motions. Median values for W0-and W51-OMFs are included in this figure. For the 15-story steel OMFs, the larger elastic values are obtained in the upper stories of the structure. This is because the effect of higher modes. As seen in the previous section, the W51-OMF has larger in the upper stories than W0-OMF because of the increase in member sections in the lower stories when wind loads are taken into consideration. The dispersion of in the upper stories for the W51-OMF is larger than that for the W0-OMF. The larger dispersion of in the upper stories for the W51-OMF may be attributed to the effect of larger lateral drifts in the upper stories and more significant contributions of higher modes.

Seismic Loss of Steel Ordinary Moment Frames
In this study, the effect of wind loads on the seismic performance of steel OMFs was investigated. To achieve the research objective, the seismic loss estimation was conducted according to FEMA P58 [30]. In seismic loss estimation, the various levels of damage states for structural and nonstructural elements should be considered. As noted that in previous sections, collapse limit state is only considered to evaluate the seismic performance of OMFs, which may overestimate the seismic performance [49].
Three limit states are considered to evaluate the seismic loss: (1) the limit state, in which damage to structural and nonstructural components occurs and results in repair

Seismic Loss of Steel Ordinary Moment Frames
In this study, the effect of wind loads on the seismic performance of steel OMFs was investigated. To achieve the research objective, the seismic loss estimation was conducted according to FEMA P58 [30]. In seismic loss estimation, the various levels of damage states for structural and nonstructural elements should be considered. As noted that in previous sections, collapse limit state is only considered to evaluate the seismic performance of OMFs, which may overestimate the seismic performance [49].
Three limit states are considered to evaluate the seismic loss: (1) the limit state, in which damage to structural and nonstructural components occurs and results in repair costs (Repair); (2) the limit state in which collapse does not occur but requires demolition and reconstruction due to severe damage (Demolition); and (3) the limit state in which collapse occurs due to earthquake (Collapse). Considering each limit state, the total expected seismic loss for a specific ground motion intensity measure (im) in the building (E(L T |im)) can be evaluated based on the following equation proposed by Ramirez and Miranda [50]: where E(L R |R, im), E(L D |D, im), and E(L C |C, im) are the expected seismic loss associated with repair, demolition, and collapse limit state at specific im, respectively. In this equation, P(C|im) is the probability that collapse occurs at a given im; P(NC|im) is the non-collapse probability [1 − P(C|im)]; and P(D|im) is the probability of occurrence of the demolition limit state, which is defined using residual inter-story drift ratio (RIDR) [50]. The seismic loss required to repair structural and nonstructural component damage caused by a ground motion was calculated according to FEMA P58 [30]. A reconstruction procedure is required to restore the structure experiencing collapse. For the demolition limit state, costs are required for the demolition of the building, the removal of debris, and the reconstruction to restore the performance of the building. In this study, the seismic loss associated with demolition limit state is assumed to be equivalent to 125% of the total replacement cost [30]. Therefore, the expected seismic losses associated with demolition and collapse limit states are assigned as 125% and 100% of the total replacement cost, respectively. Individual damage states and repair costs for structural and nonstructural components are defined according to FEMA P58 Vol 3 [51]: supporting electronic materials and background documentation [51]. The total replacement cost of the model buildings was assumed to be 2691 USD per square meter [52]. The quantities and types of structural components considered in this study are briefly summarized in Table 3, and the quantities and types of nonstructural components are defined using normative quantity estimation tool provided by FEMA P58 [51]. The expected annual loss (EAL) can be calculated using E(L T |im) (Equation (4)) and the seismic hazard curves representing the frequency of occurrence of a specific im (λ im ). To construct the seismic hazard curve, the OMFs are assumed to be located at 2393 E Tropicana Ave. in Las Vegas, NV (36.099 • N, −115.119 • W), which used a similar design spectral acceleration as in the seismic design for OMFs. The hazard curve is obtained from the USGS unified hazard tool (https://earthquake.usgs.gov/hazards/interactive/, accessed on 15 February 2022). The EAL value is calculated using Equation (5) as follows: where is the absolute value of the slope of the seismic hazard curve. In this study, nonlinear response analyses were performed for OMF buildings using 44 far-field ground motions scaled to multiple hazard levels [29]. Figure 9 shows the median peak story drift ratio (SDR) of 15-story OMFs designed with and without consideration of wind loads for three different hazard levels. The distribution of SDR for other OMFs is plotted in Appendix A ( Figure A1). Appl In this study, nonlinear response analyses were performed for OMF buildings using 44 far-field ground motions scaled to multiple hazard levels [29]. Figure 9 shows the median peak story drift ratio ( ) of 15-story OMFs designed with and without consideration of wind loads for three different hazard levels. The distribution of for other OMFs is plotted in Appendix A ( Figure A1). As shown in Figure 9, the larger occurs in the upper stories than in the lower stories irrespective of different seismic hazard levels, which is due to the contributions of higher modes. A similar observation is made in elastic range (Figure 8).
The distribution of along the building height affects the damage in the driftcontrolled nonstructural and structural components. For OMFs designed with consideration of wind loads, the concentration of higher in the upper stories becomes more significant than that for the OMF designed without consideration of wind loads. This phenomenon is more apparent when a larger wind load (51 m/s) is used in the design. The W44-OMF has 8.6%, 22.8%, and 11.9% larger than the W0-OMF for the SLE, DBE, and MCE hazard levels, respectively, where is the maximum value among the s of individual stories. The increase in is more prominent for the W51-OMF, As shown in Figure 9, the larger SDR occurs in the upper stories than in the lower stories irrespective of different seismic hazard levels, which is due to the contributions of higher modes. A similar observation is made in elastic range (Figure 8).
The distribution of SDR along the building height affects the damage in the driftcontrolled nonstructural and structural components. For OMFs designed with consideration of wind loads, the concentration of higher SDR in the upper stories becomes more significant than that for the OMF designed without consideration of wind loads. This phenomenon is more apparent when a larger wind load (51 m/s) is used in the design. The W44-OMF has 8.6%, 22.8%, and 11.9% larger SDR max than the W0-OMF for the SLE, DBE, and MCE hazard levels, respectively, where SDR max is the maximum value among the SDRs of individual stories. The increase in SDR max is more prominent for the W51-OMF, which has 27.4%, 42.9%, and 38.6% larger SDR max than the W0-OMF for the SLE, DBE, and MCE hazard levels, respectively. This indicates that the increase in member sections in the lower stories of OMFs due to the consideration of the wind loads may adversely affect not only the collapse risk but also the drift-controlled nonstructural and structural components. A similar observation is made for OMFs with different numbers of stories ( Figure A1). Figures 10 and 11 show repair costs associated with the structural, displacement-based, and acceleration-based nonstructural components in the individual stories for the 15-story OMFs at SLE and DBE hazard levels. The median repair costs normalized by the replacement cost for SLE and DBE hazard levels are only plotted in Figures 10 and 11, where STR, NSTR-D, and NSTR-A represent mean repair cost associated with structural component, drift-controlled nonstructural component, and acceleration-controlled nonstructural component, respectively.  As shown in Figure 10, for the repair cost induced by ground motions for the SLE hazard level, the contribution of NSTR-D is dominant. The W44-and W51-OMFs have 18% and 75% larger expected repair costs for the SLE hazard level than the W0-OMF. Due to damage concentration in the upper stories for the W44-OMF and W51-OMF, the repair costs associated with sharply increases. As mentioned earlier, this phenomenon is attributed to the increase in member sections in the lower stories for these OMFs due to wind loads. However, the repair costs for the lower stories are higher for the W0-OMF than for the W44-OMF and W51-OMF. Figure 11 shows the repair costs for the DBE hazard level. As shown in the figure, the repair costs associated with the NSTR-D is also dominant, as observed for the SLE hazard level. More repair costs for DBE hazard level are required than those for the SLE hazard  As shown in Figure 10, for the repair cost induced by ground motions for the SLE hazard level, the contribution of NSTR-D is dominant. The W44-and W51-OMFs have 18% and 75% larger expected repair costs for the SLE hazard level than the W0-OMF. Due to damage concentration in the upper stories for the W44-OMF and W51-OMF, the repair costs associated with sharply increases. As mentioned earlier, this phenomenon is attributed to the increase in member sections in the lower stories for these OMFs due to wind loads. However, the repair costs for the lower stories are higher for the W0-OMF than for the W44-OMF and W51-OMF. Figure 11 shows the repair costs for the DBE hazard level. As shown in the figure, the repair costs associated with the NSTR-D is also dominant, as observed for the SLE hazard level. More repair costs for DBE hazard level are required than those for the SLE hazard As shown in Figure 10, for the repair cost induced by ground motions for the SLE hazard level, the contribution of NSTR-D is dominant. The W44-and W51-OMFs have 18% and 75% larger expected repair costs for the SLE hazard level than the W0-OMF. Due to damage concentration in the upper stories for the W44-OMF and W51-OMF, the repair costs associated with SDR sharply increases. As mentioned earlier, this phenomenon is attributed to the increase in member sections in the lower stories for these OMFs due to wind loads. However, the repair costs for the lower stories are higher for the W0-OMF than for the W44-OMF and W51-OMF. Figure 11 shows the repair costs for the DBE hazard level. As shown in the figure, the repair costs associated with the NSTR-D is also dominant, as observed for the SLE hazard level. More repair costs for DBE hazard level are required than those for the SLE hazard level. For the DBE hazard level, the W44-and W51-OMFs have 12% and 39% larger repair costs than the W0-OMF, respectively. The contribution of STR on the repair costs for the DBE hazard level becomes more significant than that for the SLE hazard level. It is noted that NSTR-A did not have a significant effect on the repair costs associated with STR and NSTR-D. This is because the peak floor acceleration (PFA) is small in the steel OMFs designed for the SDC level of C max or lower.
To evaluate seismic loss for all limit states according to im, normalized expected seismic loss is plotted in Figure 12, where the horizontal and vertical axes represent the spectral acceleration normalized to the spectral acceleration corresponding to the maximum considered earthquake (S MT ) and the expected seismic loss (EL) normalized to the total replacement cost, respectively. Figure 13 shows the contribution of each limit state (repair, demolition, and collapse limit states) on EL. This figure shows the contributions of each limit state on the expected seismic loss according to im.  As shown in Figure 12, the values of the W0-OMF for individual limit states are less than those of W44-and W51-OMFs irrespective of . The W51-OMF has larger values for individual limit states than W44-OMF. This indicates that the seismic loss could increase when wind loads are considered in the design. As mentioned earlier, the concentrations of lateral displacement in the upper stories for OMFs designed with consideration of wind loads cause a rapid increase in the maximum in the upper stories, which not only accelerates the damage of structural and nonstructural components but also causes an increase in and an increase in the probability of sideway collapse. The of the W44-OMF is at least 7% larger than that of the W0-OMF. The probability of collapse and demolition limit states rapidly increases with an increase in . The value increases up to 25%. For the W51-OMF, the value increases by at least 50% as  As shown in Figure 12, the values of the W0-OMF for individual limit states are less than those of W44-and W51-OMFs irrespective of . The W51-OMF has larger values for individual limit states than W44-OMF. This indicates that the seismic loss could increase when wind loads are considered in the design. As mentioned earlier, the concentrations of lateral displacement in the upper stories for OMFs designed with consideration of wind loads cause a rapid increase in the maximum in the upper stories, which not only accelerates the damage of structural and nonstructural components but also causes an increase in and an increase in the probability of sideway collapse. The of the W44-OMF is at least 7% larger than that of the W0-OMF. The probability of collapse and demolition limit states rapidly increases with an increase in . The value increases up to 25%. For the W51-OMF, the value increases by at least 50% as As shown in Figure 12, the EL values of the W0-OMF for individual limit states are less than those of W44-and W51-OMFs irrespective of im. The W51-OMF has larger EL values for individual limit states than W44-OMF. This indicates that the seismic loss could increase when wind loads are considered in the design. As mentioned earlier, the concentrations of lateral displacement in the upper stories for OMFs designed with consideration of wind loads cause a rapid increase in the maximum SDR in the upper stories, which not only accelerates the damage of structural and nonstructural components but also causes an increase in RIDR and an increase in the probability of sideway collapse. The EL of the W44-OMF is at least 7% larger than that of the W0-OMF. The probability of collapse and demolition limit states rapidly increases with an increase in im. The EL value increases up to 25%. For the W51-OMF, the EL value increases by at least 50% as compared twith that for the W0-OMF. The largest increase in the EL value is 78%. Figure 13 shows the contribution (%) of repair, demolition, and collapse to EL as a function of im. As seen in this figure, seismic loss associated with the repair limit state is dominant on EL at low hazard levels. However, the contribution of seismic loss associated with demolition and collapse grows with an increase in im. As shown in Figure 6, OMFs experience significant damages in few stories and collapse. Therefore, the loss associated with demolition and collapse surpasses that associated with repair when im exceeds the MCE hazard level. For im corresponding to the MCE hazard level (S MT ), seismic loss associated with demolition, and collapse limit states become more than 50% the total EL. The contribution (%) of demolition and collapse on EL is the largest for the W51-OMF and is the least for the W0-OMF. This observation can be expected because OMFs designed with consideration of wind loads have higher seismic risk than corresponding OMFs designed without the consideration of wind loads as described in Section 4. Figure 14 shows the expected annual loss (EAL) of W0-, W44-, and W51-OMFs. As in Figure 13, EAL normalized by replacement cost for repair, demolition, and collapse limit states are marked. Appl Figure 14 shows the expected annual loss ( ) of W0-, W44-, and W51-OMFs. As in Figure 13, normalized by replacement cost for repair, demolition, and collapse limit states are marked. As was observed in , the value also increases when wind loads are considered. The associated with repair also significantly increases for OMFs designed with consideration of wind loads. The value associated with repair is more than 50% of the total . This is because of the frequency of ground motion with low intensity. The total values of 12-, 15-, and 18-story W44-OMFs are 14%, 20%, and 45% larger than those of corresponding W0-OMFs, respectively. An increase in is more significant for W51-OMF. The total s of 12-, 15-, and 18-story W51-OMFs are 48%, 58%, 87%, and 76% larger than those of corresponding in 9-, 12-, 15-, and 18-story W0-OMFs, respectively. In this regard, if seismic performance evaluation is conducted for new buildings designed without consideration of wind loads, unconservative results can be obtained.

Conclusions
In this study, the seismic risk of steel ordinary moment frames (OMFs) designed with and without the consideration of wind loads were evaluated to investigate the effect of wind loads on the seismic performance of the steel OMFs. The seismic collapse intensities and seismic loss were estimated according to FEMA P695 and FEMA P58 methodologies, respectively. Nine-, twelve-, fifteen-, and eighteen-story steel OMFs were considered, and wind load of 0, 44, and 51 m/s were used for the design of OMFs. The following conclusions are obtained from this study.
1. OMFs designed with consideration of gravity, seismic, and wind loads (W44-and W51-OMF) have larger member sections than OMFs designed when only gravity and As was observed in EL, the EAL value also increases when wind loads are considered. The EAL associated with repair also significantly increases for OMFs designed with consideration of wind loads. The EAL value associated with repair is more than 50% of the total EAL. This is because of the frequency of ground motion with low intensity.

Conclusions
In this study, the seismic risk of steel ordinary moment frames (OMFs) designed with and without the consideration of wind loads were evaluated to investigate the effect of wind loads on the seismic performance of the steel OMFs. The seismic collapse intensities and seismic loss were estimated according to FEMA P695 and FEMA P58 methodologies, respectively. Nine-, twelve-, fifteen-, and eighteen-story steel OMFs were considered, and wind load of 0, 44, and 51 m/s were used for the design of OMFs. The following conclusions are obtained from this study.

1.
OMFs designed with consideration of gravity, seismic, and wind loads (W44-and W51-OMF) have larger member sections than OMFs designed when only gravity and seismic loads (W0-OMF) are considered, particularly for interior columns and beams in the lower stories.

2.
The period-based ductility (µ T ) of W44-and W51-OMFs was larger than that of corresponding W0-OMFs. However, values of overstrength factor (Ω) of W44-and W51-OMFs are similar with those of W0-OMFs. All steel OMFs had the Ω values smaller than 3.0, which is specified for steel OMFs in ASCE 7-16 [1]. This is because in OMFs, failure occurred in the upper stories before reaching the intended maximum base shear capacity expected in the design.

3.
Due to the increase in member sections only in the lower stories for W44-and W51-OMFs, lateral drift and damage concentrated in the upper stories. Although W44-and W51-OMFs had larger member sections than corresponding W0-OMFs, those frames had higher collapse risk than corresponding W0-OMFs.

4.
Seismic losses for W51-and W44-OMFs were larger than that for corresponding W0-OMFs at all hazard levels. Seismic losses associated with all limit states for W44and W51-OMFs were also larger than corresponding W0-OMFs, which was more apparent with an increase in building height and the magnitude of wind loads.

5.
The EL values of the W0-OMF for individual limit states are less than those of W44and W51-OMFs irrespective of im. W51-OMFs has larger EL values for individual limit states than W44-OMFs. This indicates that the seismic loss could increase when wind loads are considered in the design. The EL of the W44-OMF is at least 7% larger than that of the W0-OMF. For the W51-OMF, the EL value increases by at least 50% as compared with that for the W0-OMF. The largest increase in the EL value is 78%. 6.
As was observed in EL, the EAL value also increases when wind loads are considered. The total EAL values of 12-, 15-, and 18-story W44-OMFs are 14%, 20%, and 45% larger than those of corresponding W0-OMFs, respectively. An increase in EAL is more significant for W51-OMFs. The total EALs of 12-, 15-, and 18-story W51-OMFs are 48%, 58%, 87%, and 76% larger than those of corresponding W0-OMFs, respectively. In this regard, if seismic performance evaluation is conducted for new buildings designed without consideration of wind loads, unconservative results can be obtained. 7.
Based on this study, it is important to consider the effect of wind loads when evaluating the seismic performance of steel structures. Otherwise, evaluation results could be unconservative.

Conflicts of Interest:
The authors declare no conflict of interest.