Abstract

Ductility plays a central role in seismic analysis and design of steel buildings. A numerical investigation regarding the evaluation of energy dissipation, ductility, and ductility reduction factors for local, story, and global structural levels is conducted. Some steel buildings and strong motions, which were part of the SAC Steel Project, are used. Bending local ductility capacity (µLϕ) of beams can reach values of up to 20, as shown in experimental investigations. The values are larger for medium than for low-rise buildings, reflecting the effect of the structural complexity on µLϕ. Most of the dissipated energy occurs on beams; however, resultant stresses at columns are also significantly reduced by beam yielding. A value of 1/3 is proposed for the ratio of global to local ductility; thus, if local ductility capacity is stated as the basis for the design, global ductility capacity can be calculated by using this ratio. It is implicitly assumed in seismic codes that the magnitude of the global ductility reduction factor is about 4; according to the results found in this paper, it is not justified; a value of 3 is observed to be more reasonable. According to the well-known ratio of the ductility reduction factor to ductility, this ratio should be unity for the models under consideration; the results of this study indicate that, for global response parameters, a value of 3/4 is more appropriate and that, for local response parameters, values larger than 2 can be reached; the implication of this is that using simplified methods like the static equivalent lateral force may result in nonconservative designs from a global structural point of view, but in conservative designs from a local point of view. A value of 8 is proposed for the ratio of the global ductility reduction factor to the global normalized energy.

1. Introduction

Even though building structures undergo significant nonlinear deformations when subjected to strong earthquakes, simple elastic procedures are still used to determine the seismic demands (International Building Code (IBC) [1], National Building Code of Canada (NBCC) [2], Mexico Federal District Code (MFDC) [3], and Eurocode 8 (EC) [4]). Simplified methods like the Static Equivalent Lateral Force (SELF) are broadly used for regular buildings with relatively short periods (low- and medium-rise). However, while using this procedure in steel buildings, it is not possible to properly capture the effects of nonlinearity introduced by large deformations, by the connections, and by nonlinear geometry. In addition, dissipation of energy due to yielding of the material is considered in a very crude way. According to the method, static analysis of the buildings subjected to equivalent lateral forces, which are related to the properties of the structure and the seismicity of the region, provides the design forces. In the procedure, the ductility parameter (µ) plays an important role in the determination of the design seismic forces since it is associated with the energy dissipation structural capacity produced by nonlinear behavior, allowing for a reduction of the elastic strength demands; the larger the ductility, the smaller the design seismic forces. In fact, the magnitude of the reduction of the elastic design seismic forces directly depend on a parameter that, in general, can be called “the reduction factor,” which in turn greatly depends on an associated parameter called the ductility reduction factor (Rµ) [5, 6]. It is particularly important for steel structures since the beneficial effect of ductility and energy dissipation is supposed to come from different sources. However, there is not unanimity on the profession on how to define it; it is argued that this parameter is constantly used in the profession in an indirect way to estimate the building seismic design forces, but there is no engineering definition of it in our specifications [7, 8]. Defining the magnitude of these parameters represents one of the most controversial issues in the SELF procedure.

The reduction factor receives the name of the response modification factor (R), force modification factor (Rd), seismic reduction factor (Q′), and seismic reduction factor (q′) in the IBC, NBCC, MFDC, and EC codes, respectively. It is implicitly or explicitly assumed in these codes that ductility represents the capacity of the structure to dissipate energy and that the reduction of the elastic demands importantly depends on the ductility capacity.

Thus, the ductility parameter and the ductility reduction factors play a central role in the seismic design of steel buildings; for that reason, it has been an important research topic during the last recent decades. However, there are many aspects that need additional attention; some of them are addressed in this research. As it is further elaborated below, the central objective of this paper is to evaluate the ductility parameter, the ductility reduction factor (Rµ), for steel buildings with moment-resisting frames (MRF) for different structural levels (local, story, and global) relating them with the dissipated energy. The evaluation of the µ and Rµ parameters according to sophisticated nonlinear analysis procedures to capture the effects of the mentioned sources of nonlinearity is needed. It is accepted that nonlinear time history analysis is the most accurate and reliable analysis procedure providing a realistic modeling of the structure as well as of the cyclic load deformation characteristics of its structural elements. Extensive seismic nonlinear time-history analyses are performed in this investigation to reach the objectives of the study. The results are compared with those specified in the codes.

2. Literature Review and Objectives

Investigations regarding nonlinear analysis of buildings under the action of earthquakes, following different objectives, have been conducted by many researchers during the last decades. In this regard, since this dissipation of energy allows for a reduction of the elastic seismic forces, the quantification of the ductility demand, the ductility reduction factor, and the force (or seismic) reduction factor is of particular interest. There have been many studies considering concrete or steel buildings modeled as SDOF or simple systems. Introduced first in ATC [9] at the end of the 70s, the modification factor was used to reduce the elastic response (base shear) obtained from a 5% damped acceleration response spectra. Among the first investigations, we can also find those of Newmark and Hall [10], and they proposed a procedure to relate Rμ and μ based on the basic elastic design spectra. Hadjian [11] calculated the spectral accelerations considering the nonlinear deformation of structures. Miranda and Bertero [12] proposed simplified expressions to estimate the inelastic design spectra as a function of the maximum tolerable ductility. Whitaker et al. [13] proposed calculating R as the product of factors accounting for viscous damping, ductility, and overstrength. Arroyo-Espinoza and Teran-Gilmore [14] from the dynamic response of SDOF systems proposed expressions to calculate Rµ. Levy et al. [15] derived approximate harmonic equivalent stiffness and damping for bilinear systems in the context of earthquake resistant. Karmakar and Gupta [16] by using elastoplastic oscillators performed a parametric study to estimate the dependence of strength reduction factors on strong motion duration, earthquake magnitude, geological site conditions, and epicentral distance. Rupakhety and Sigbjörnsson [17] presented ground-motion prediction equations which describe constant-ductility inelastic spectral ordinates and structural behavior factors to be applied within the framework of Eurocode 8. Sanchez-Ricart [18] reviewed the backgrounds that support the values of the reduction factor in the United States, Europe, and Japan. Halabian and Kabiri [19] evaluated the effect of the foundation flexibility on the ductility reduction factors of reinforced concrete stack-like structures. AlHamaydeh et al. [20] investigated the force reduction factors for reinforce concrete frame buildings designed according to the IBC code. Naimi et al. [21] studied the reduced beam section approach via the introduction of multilongitudinal voids in the beam web for various beam depths, where the ANSYS finite element program was used in the numerical simulation. An improvement in the connection ductility was observed since the input energy was dissipated uniformly along the beam length. Azimi et al. [22] proposed a new beam to column connection and numerically and experimentally showed that the new connection could improve the ultimate structural lateral resistance, ductility, and energy dissipation capacity. Fanaie and Shamlou [23] calculated the overstrength, ductility, and response modification factors for frames braced with a different type of buckling restrained braces. Zhai et al. [24] investigated the strength reduction factor of SDOF systems with constant ductility performance subjected to the mainshock-aftershock sequence-type ground motions. Wang et al. [25] studied the nonlinear behavior of prestressed steel reinforced concrete beams by using finite element analysis. Based on the model, they analyzed the effect of prestressed force to the stiffness, the ultimate bearing capacity, and ductility of the beams. Cho et al. [26] investigated the shear design equations for prestressed hollow-core slabs and examined the magnitude of the strength reduction factors based on the structural reliability theory. The results showed that the shear strengths of the members with the heights of more than 315 mm are excessively safe, whereas some members with low depths did not satisfy the target reliability index.

There are also several studies regarding the evaluation Rμ and μ factors for multidegree of freedom (MDOF) systems. Nassar and Krawinkler [27] studied the relationship between force reduction factors and ductility for SDOF and simplified (three-story single-bay) MDOF systems. Moghaddam and Mohammadi [28] introduced a modification to the response modification factor and proposed an approach to evaluate the seismic strength and ductility demands of MDOF structures. Elnashai and Mwafy [29] investigated the relationship between the lateral capacity, the design force reduction factor, the ductility factor, and the overstrength factor for reinforced-concrete buildings. Reyes-Salazar [8] studied the ductility capacity of plane steel moment-resisting frames. Medina and Krawinkler [30] presented an evaluation on drift demands for regular moment-resisting frame structures subjected to ordinary ground motions considering the uncertainty due to differences in the frequency content of the ground motions. De Stefano et al. [31] studied the effects of the overstrength on the seismic behavior of multistory asymmetric buildings. Cai et al. [32] estimated ductility reduction factors for MDOF systems by modifying ductility reduction factors of SDOF systems through a modification factor. Chopra [33] studied the ductility reduction factors for MDOF systems modeled as shear buildings and their corresponding equivalent SDOF systems. Mollaioli and Bruno [34] developed constant ductility spectra for SDOF and MDOF systems. Vielma et al. [35] studied the seismic damage and safety assessment of buildings with low ductility used in Spain by using pushover and dynamics analysis. Ceylan et al. [36] estimated the strength reduction factor for prefabricated industrial structures having a single story, one, and two bays. Honglue el al. [37] estimated the accuracy of strength reduction factors proposed by other investigators by using a suite of near-fault earthquake records with directivity-induced pulses. Ganjavi and Hao [38] studied the seismic response of linear and nonlinear MDOF systems subjected to a group of earthquakes recorded on alluvium and soft soils. Serror et al. [39] evaluated the values of both damping and ductility reduction factors for steel moment-resisting frames with supplemental linear viscous dampers. Bojorquez et al. [40] assessed the effect of cumulative damage on the strength requirements of degrading structures through the evaluation of the target ductility and corresponding strength reduction factors of simple degrading structures. The results provide insight into all relevant parameters that should be considered during seismic design of earthquake-resistant structures, and some recommendations to evaluate the effect of cumulative damage on seismic design are suggested. Chaulagain et al. [41] by using pushover analysis studied the effects of overstrength on the ductility factor of reinforced concrete buildings. Reyes-Salazar et al. [42] studied the ductility reduction factor for buildings with moment-resisting steel frames (MRSF) which were modeled as complex MDOF systems, considering an intermediate level of inelastic structural deformation. Vuran and Aydınoğlu [43] developed simple capacity and ductility demand estimation tools for coupled core wall systems to be implemented during the preliminary design stage of such structural system. Gómez-Martínez et al. [44] analytically studied the local and global ductility of wide-beam reinforced concrete moment-resisting frames. Liu et al. [45], for the SDOF system and considering spatially varying ground motions, developed a new response spectrum method by incorporating the ductility factor and strain rate into the conventional response spectrum method. Hashemi et al. [46] presented the results of studies on two important seismic parameters, namely, ductility and response modification factor for moment-resisting frames with concrete-filled steel tube columns.

The abovementioned studies represent a significant contribution regarding the evaluation of ductility or ductility reduction factors; however, in most of them, SDOF systems or a limited level of inelastic deformation were considered. Therefore, they did not explicitly consider the energy dissipation associated with nonlinear behavior of the structural elements existing in actual systems. It has been shown [8,4749] that ductility demands and the ductility reduction factors depend on the amount of dissipated energy, which in turn depends on the pattern of plastic hinges formed in the frames as well as on the loading, unloading, and reloading process at plastic hinges. In addition, a limited level of inelastic deformation is not associated with the ductility capacity. Moreover, local, story, and global ductility, as well as local, story and global ductility reduction factors, as well as relationships among them taking into account the dissipated energy, have not been studied.

The objectives of this research are as follows:(1)Calculate local ductility demands (in terms of curvatures) and capacities for individual structural elements (beams and columns) as well as story and global ductility demands and capacities. The ratio of local to global ductility is also estimated.(2)Calculate the ductility reduction factors as well as the ratio of the ductility reduction factor to ductility, for the three structural levels under consideration, and compare them with those specified in the codes.(3)Estimate the energy demands for local, story, and global levels as well as the ratio of the ductility reduction factor to dissipated energy.

3. Methods, Procedure, and Structural Models

3.1. Steel Building Models

As part of the SAC Steel Project [50], several steel building models were designed, which are supposed to satisfy the code requirements at the time the project started, for the following cities: Los Angeles (IBC), Seattle (IBC), and Boston (BOCA) [51]. The perimeter moment-resisting frames (PMRF) of the 3- and 10-level buildings used in the project located in the Los Angeles area are considered in this study to address the issues discussed earlier. Isometric views of the buildings are shown in Figures 1 and 2, where the PMRF can be easily identified (exterior frames); they are referred hereafter as Models 1 and 2, respectively. The beams and columns of the models were designed following the strong column-weak beam (SC-WB) concept; sizes of beams and columns for the PMRF, as reported, are given in Table 1. The first three translational periods of the plane frames (lateral vibrations) are 1.03 s, 0.30 s, and 0.15 s, for the 3-level model; the corresponding values for the 10-level model are 2.41 s, 0.89 s, and 0.5 s. The damping is considered to be 3% of the critical. Each column is represented by one element and each girder by two elements, having a node at the midspan. The RUAUMOKO software [52] is used to perform the required step-by-step nonlinear seismic analysis.

It is worth to mention that some results of the SAC Steel Project are given in a research report [50] where many aspects of the seismic performance of the mentioned steel frame models subjected to earthquake ground motions are presented. The strong motions were scaled according to several return periods; story drifts of about 1%, 2%, 3%, and 4% and a little larger (between 5% and 6%) were developed in the models.

The frames are modeled as complex 2D MDOF systems, having three degrees of freedom per node. The Newmark constant average acceleration method, lumped mass matrix, Rayleigh damping, and large displacement effects are considered within the Ruaumoko Computer Program environment while performing the required nonlinear seismic analysis; the time increment in the analysis was 0.01 s. The panel zone was considered to be rigid. Typical input data as ground accelerations, boundary conditions, node coordinates, and elastic and inelastic section properties are given or read within the computer program. No strength degradation member, bilinear behavior with 5% of the initial stiffness in the second zone, and concentrated plasticity were assumed. The interaction axial load-bending moment is given by the yield interaction surface proposed by Chen and Atsuta [53].

3.2. Earthquake Loading

When a structure is subjected to the action of two different strong motions, even when they are normalized with respect to the same peak ground acceleration or with respect to any other parameter, it is expected to respond differently, reflecting the influence of the frequency content of the motions and of the structural vibration modes. Thus, to get meaningful results, the models under consideration are excited by twenty strong motions in time domain whose characteristics are given in Table 2. The predominant periods of the strong motions range from 0.11 to 0.62 sec. The earthquake time histories were obtained from the Data Sets of the National Strong Motion Program (NSMP) of the United States Geological Surveys (USGS). For a given direction, half of the seismic loading and the tributary gravity loading are assigned to the corresponding PMRF. In order to have different levels of structural deformation as well as moderate and significant nonlinear behavior, the strong motions are scaled in terms of Sa evaluated at the fundamental lateral vibration period (Sa(T1)) ranging from 0.4 g to 1.2 g for the 3-level model and from 0.2 g to 0.6 g for the 10-level model, with increments of 0.2 g. Therefore, the maximum values of the seismic intensities are Sa = 1.2 g and 0.6 g for the 3- and 10-level models, respectively. Thus, considering two models, twenty strong motions, two horizontal seismic components, and six seismic intensities, in addition to the elastic analyses, more than one thousand of nonlinear analyses of MDOF systems were performed.

The abovementioned maximum levels of seismic intensities produce a deformation state which is very close to that of a collapse mechanism, where interstory drifts of about 5% were developed for some strong motions; therefore, they are associated with the structural and ductility capacity. This is concordant with the results of some experimental studies where it has been shown that moment-resisting steel frames may undergo interstory displacements of up to 5% (and even larger) and still be able to vibrate in a stable manner [5458]. Thus, the ductility values obtained for this level of deformation, according to the different definitions, are assumed to be the ductility capacity.

It is important to note that UBC-1994 in Sections 1629.1 and 1629.2 states: “Dynamic analysis procedures, when used, shall conform to the criteria established in this section. The analysis shall be based on an appropriate ground motion representation … .The ground motions representation shall, as a minimum, be one having a 10% probability of being exceeded in 50% … .” The expression “as a minimum” implies that larger intensities (and larger drifts) of the strong motions can be used. In addition, according to the particular objectives stated in our paper, we need a deformation state close to that of a collapse mechanism which is developed for drifts of about 5%, for some strong motions.

In addition to the seismic loading, the following gravity loads [50] were used in the analysis: (a) the floor dead load for weight calculations was 96 psf; (b) the floor dead load for mass calculations was 86 psf; (c) the roof dead load was 83 psf; (d) the reduced live load per floor and for roof was 20 psf. The seismic mass for the entire structure was as follows: (a) for the roof of the 3-level building, it was 70.90 kips-sec2/ft; (b) for floor 2 of the 3-level building, it was 65.53 kips-sec2/ft; (c) for the roof of the 10-level building, it was 73.10 kips-sec2/ft; (d) for floor 2 of the 10-level building, it was 69.04 kips-sec2/ft; (e) for floors 3 to 9 of the 10-level building, it was 67.86 kips-sec2/ft.

3.3. Ductility Definitions

The discussions made in Sections 1 and 2 clearly indicate that the parameter ductility is a concept of central importance in seismic analysis and design of steel buildings. In this regard, the different types of ductility existing in a steel building must be considered [10]. Within a SDOF system context, ductility is defined as the ratio of the maximum inelastic displacement (Dmax) to the yield displacement (Dy). Dmax is calculated as the maximum displacement that the system undergoes during the application of the seismic loading and Dy as the displacement of the system when yielding occurs on it for the first time. For MDOF systems, however, it is not clearly stated how to define these two parameters (Dmax and Dy); in addition, and as commented above, it is necessary to properly consider the different levels of ductility (local, story, and global). It is generally accepted that local ductility is larger than story ductility, which in turn is larger than global ductility. It must be noted that ductility demand is different from ductility capacity. For example, for the case of a particular story, story ductility demand can be defined as the ratio of the maximum relative lateral displacement (drift) of the story during the application of the seismic loading to the corresponding drift when first yielding occurs at any member of the story, while ductility capacity is the ratio of the maximum permissible inelastic drift to the drift when first yielding occurs. Ductility capacity is usually obtained from experimental results for individual members (local ductility). For that reason, some researchers [59] suggest using local ductility as the basis for design because there are numerous laboratory studies on ductility for members. Therefore, it is important to relate the local to the story, or to the overall structural ductility (global ductility).

Theoretically, ductility capacity should be reached when a collapse mechanism develops in the structure. To obtain this, it needs to be guaranteed that plastic moments are reached at positions of maximum moments before failure due to instability, namely, local buckling or lateral torsional buckling, in a member or in a connection occurs. For the case of steel buildings, local ductility (µLϕ) for a flexural member will be associated with the rotational deformation of the member. It is defined, for a given joint, as the ratio of the maximum inelastic curvature that the joint undergoes during the total time of excitation (ϕmax) to the curvature of the joint when it yields for the first time (ϕy). Mathematically, it is expressed as

Thus, as soon as any of the joints of a given member yields for the first time, the corresponding curvature is identified as ϕy for that particular member. In a similar manner, the curvature is calculated at each time increment of the analysis, and the largest one is identified as ϕmax for the node and member under consideration. This type of ductility demand will be calculated for beams as well as for those columns where plasticization is produced mainly by bending moments; as it will be further discussed below this is the case for most columns, particularly for the 10-level building.

Story ductility is defined in terms of lateral drifts. The ductility of a story (µS) is defined as the ratio of the maximum inelastic drift of the story during the total time of excitation (Δmax) to the drift of the story when any of its members yields for the first time (Δy). Mathematically, we have

Regarding the Δy parameter in Equation (2), it is assumed that, for a given story, the beams and columns connecting beneath it are part of the story. For example, for the 3-story model (top story), going from the top to the bottom, the first three beams and four columns are considered to be part of the third story; in the same manner, the second set of three beams and four columns are considered to be part of Story 2, and so on. Then, in the µS definition, the expression “the drift of the story when any of its members yields for the first time” refers to first yielding of any beam or column that is part of the story under consideration.

Since global ductility should represent the overall structural inelastic deformation, some researchers suggest defining it in terms of relative lateral displacements [5, 10, 58]. Here, it is defined as the mean value of the story ductilities; mathematically, we havewhere n is the number of stories.

3.4. Ductility Reduction Factor

As for the ductility parameter, the ductility reduction factors are estimated for different structural levels. This parameter, in general, can be expressed aswhere Re represents the peak value of a given, local, story, or global response parameter, obtained from elastic analyses (without considering dissipation of energy) and Ri represents the same, but nonlinear analyses (considering dissipation of energy) are considered instead. For the case of bending moments, which is a local response parameter, Equation (4) takes the following particular form:where RµLϕ is the bending local ductility reduction factor and Me and Mi are the elastic and inelastic peak bending moments, respectively, at a given joint of a given member. Equation (5) will be used to calculate the reduction in bending moments for both beam and columns. Reductions in axial loads on columns are not calculated since, as stated above, the number of cases where dissipation of energy due to plasticization by the axial load occurs is very small in comparison with that of bending moments.

The story ductility reduction factor (RµS) is calculated aswhere Vse and Vsi are the elastic and inelastic peak values of the interstory shears, respectively. The global ductility reduction factors (RµG) are calculated as the mean values of RµS, and it is

3.5. Dissipated Energy

The dissipated energy is also calculated for local, story, and global levels. The normalized dissipated energy at a given joint (ELϕ) of a given beam is calculated aswhere EDϕ is the energy demand and ECϕ is the energy capacity of a member under the action of only bending moment. ECϕ in Equation (8) is estimated by the following equation [60, 61]:where MP and θpa in Equation (9) in turn represent the plastic moment and the cumulative plastic rotation capacity of a member under bending, respectively. Even though experiments to estimate the θpa parameter are not common, there is some experimental evidence to provide a reasonable value. A value of 0.23 for θpa will be used in this paper [6065].

For all the beams of a story, the average normalized dissipated energy per joint ((ELϕ)SB) can be calculated aswhere q is the number of beams in the story.

The normalized dissipated energy by bending at a column is calculated aswhere is the reduced energy capacity of a column due to the presence of axial loads. Taking into account that (a) the SC-WB concept was followed in the design of the structural models implying that most of the hysteretic dissipated energy will occur in beams by the action of bending moments, (b) experimental studies regarding the cumulative plastic rotation capacity of a member under the combined action of bending moment and axial load are rare, and (c) the level of average axial load demands observed in this study implies a reduction of the bending moment capacity in columns of about 30–50%, it is conservatively assumed in this research that the cumulative plastic rotation capacity of a member under the action of bending and axial load is one third of that of θpa. Then,

For all the columns of a story, it is obtained aswhere is the normalized energy per joint and r is the number of columns. The average normalized dissipated energy per joint in a given story (ES) is calculated as

Finally, the normalized dissipated energy for the whole frame (EG) is calculated as

All the parameters in Equation (15) were defined earlier.

4. Results and Discussion

4.1. Objective 1: Results in Terms of Ductility
4.1.1. Local Ductility

The bending local ductility parameter (Equation (1)) is calculated for each individual beam and columns of the two models for the seismic intensities of each strong motion mentioned in Section 3.2. First, the models are subjected to the simultaneous action of the horizontal seismic component oriented in the NS direction, the vertical seismic component, and the gravity loads. Then, the models are subjected to a similar set of loads, but the other horizontal (EW) seismic component is applied instead. It is assumed in this research that structural members have no problems of lateral torsional and local or shear buckling in the panel zone. As stated above in Section 3.3, bending local ductilities are expressed in terms of curvatures in this paper. Local ductilities expressed in terms of rotations are expected to be smaller.

For a given story, the µLϕ values are averaged, first over all the beams (Equations (10)), and then over all the columns (Equations (12)). It must be noted that even though the SC-WB concept was followed in the model design, plastic hinges were developed in some columns, particularly for the 10-level model and the largest intensity of the strong motions. Plots for the resulting averages are developed for each strong motion, but they are not presented, about 30 plots were developed. Only the fundamental statistics in terms of the mean values (MV), calculated over all the strong motions, are given. The results are presented in Table 3 for the 3-level model and in Figures 3(a)3(d) for the 10-level model. The term “ST” in the table and in figures stands for the story level. Results indicate that, for the beams of the 3-level building, the maximum bending ductility demands occur, in general, for the second story. For the case of columns, unlike what observed for beams, the mean values of µLϕ tend to decrease with the story number, but the mean values are much smaller than those of beams; in fact, for the two lowest intensities (Sa = 0.4 g and 0.6 g), they are essentially equal to unity for the two upper stories implying no yielding. For both, beams and columns, the mean values tend to increase with the seismic intensity and are larger for the NS than that for the EW direction. The largest observed values for beams are 11.86 and 9.32 for the NS and EW direction, respectively.

The results for the 10-level building resemble those of the 3-level building in the sense that the mean values of µLϕ are much larger for beams than for columns. However, unlike the 3-level building, the mean values tend to decrease with the story number for beams, and as for the 3-story building, the µLϕ mean values tend to decrease with the story number for columns. No yielding occurs in columns in most of the cases for the two lowest intensities of the strong motions (Sa = 0.2 g and 0.4 g). The only additional observation that can be made is that, for the larger strong motion intensities, the maximum mean values of the bending local ductility demands for beams are larger for the 10-level building than that for the 3-level building (14.3 against 11.86).

It is worth to mention that moderate yielding occurred for seismic intensities of 0.4 g and 0.2 g for the 3- and 10-level structural models, respectively; the corresponding seismic intensities for significant deformation, as stated in Section 3.2, are 1.2 g and 0.6 g, for the 3- and 10-level models, respectively. Even though, it is not shown in the paper, the drifts for these significant levels of deformation were about 5% for many of the strong motions. Thus the mean values of µLϕ shown in Table 3, or in Figures 3(a)3(d), are associated with the structural capacity and are assumed to be the bending local ductility capacity. Even though, as stated earlier, results for individual strong motions are not included in the paper, it was observed that the maximum values of µLϕ were close to 20 for some particular earthquake and stories. There is some evidence in the literature [10,6668] that the bending local ductility capacity can reach values larger than 20; however, it was for monotonic loading and individual members. It was also observed that, for a given value of Sa, the magnitude of µLϕ significantly varies from one seismic motion to another, in spite of the structural deformation, in terms of Sa for each seismic motion was the same, reflecting the effect of the frequency content on the values of the µLϕ parameter. The coefficients of variation (COV) are not shown either, but it can be said that the uncertainty in the estimation of µLϕ is moderated in most of the cases (COV∼0.20–0.45).

4.1.2. Story Ductility

Similar to the µLϕ parameter, the µS values as defined by Equation (2) are calculated for both structural models and horizontal directions, and plots for individual strong motions are developed but are not presented. The mean values for all seismic intensities and structural directions are presented in Table 4 and Figures 4(a) and 4(b), for the 3- and 10-level buildings, respectively. It can be observed that, similar to µLϕ for beams, the mean values of µS do not present any tendency with the story number for the case of the 3-level building; for the 10-level building, however, unlike the case of µLϕ, the values do not tend to decrease with height. By comparing the results of Table 3 and Figures 3(a)3(d) with those of Table 4 and Figures 4(a) and 4(b), it is observed that, for beams, µS is much smaller than µLϕ; the mean values for the maximum seismic intensity range from 3.84 to 4.28 and from 2.49 to 5.19, for the 3- and 10-level models, respectively. Both, µLϕ and µS mean values for beams tend to increase with the seismic intensity.

4.1.3. Global Ductility

The global ductility values (µG) are now discussed. The results are given in Table 5 for the two buildings, the two horizontal directions, and all seismic intensities. It can be observed that the mean values of µG tend to increase (as expected) with the strong motion intensity and that they are quite similar for the two buildings; for the demands associated with the maximum deformations (ductility capacity), µG takes values of 3.95 and 3.89, for the NS and EW directions, respectively, for the 3-level building; the corresponding values are 4.11 and 3.60 for the 10-level building. By comparing the results of Table 5 with those of µLϕ for beams (Table 3 and Figures 3(a) and 3(b)), it is noted that, as for the case µS, the values are significantly smaller for µG.

4.1.4. Ratio of Local to Global Ductility

As stated above, the estimation of ductility capacity is commonly based on experimental studies of individual members. For this reason, it is suggested [59] considering local ductility as the basis for design. In this regard, it is important to find a relationship (Q = µG/µLϕ) between global and local ductility. The µG values given in Table 5 and the µLϕ values, averaged over all the beams of the frames, are used to calculate the Q ratio. The results are summarized in the second row of Table 5. It is observed that the ratio values tend to decrease with the intensity of the earthquake. Based on the results obtained for the maximum deformation under consideration (0.36, 0.44, 0.37, and 0.35), a value of 1/3 is proposed for the Q ratio.

4.2. Objective 2: Results for Ductility Reduction Factor
4.2.1. Local Ductility Reduction Factor

The local, story, and global ductility reduction factors, as given by Equations (5)–(7), respectively, are now calculated. Bending local ductility reduction factors (RµLϕ) is discussed first. As for the µLϕ parameter, for a given story, the RµLϕ values are averaged, first over all the beams, and then over all the columns. Plots for the resulting averages, as for the other parameters, are developed for individual strong motions, but they are not presented. Only the mean values are given and discussed. The results are presented in Table 3 for the 3-level model, while those of the 10-level building are given in Figures 3(e)3(h).

It can be seen that, for beams of the 3-level building, the mean values of RµLϕ tend to decrease with the story number but tend to increase with the seismic intensity, implying that the magnitude of the reduction of the bending moments of beams produced by yielding decreases with the story number, but it increases with the strong motion intensity. This variation with the story number is not observed for column moments; in fact, the maximum reductions occur for the upper level. It is also shown that the mean values of RµLϕ are larger for beams than that for columns, and the values larger than 4 are observed for beams. However, for some cases of columns, RµLϕ can be significant even if plastic hinges are not developed. For example, for columns of Story 3 (ST = 3), NS direction and Sa = 0.6 g, RµLϕ equals 1.95; this is an interesting point since, although yielding did not occur on these columns (i.e. µLϕ = 1.00), bending moments are considerably reduced implying that yielding on beams reduce not only their bending moments but also those of columns. From the results of the 10-level building, it can be observed that the maximum values of RµLϕ occur for the two stories below the roof (the values larger than 5 are observed for Sa = 0.6 g) and that these values are larger than those of the 3-level building.

4.2.2. Story Ductility Reduction Factor

The story ductility reduction factor (RµS), calculated according to Equation (6), is now discussed. Even though, as for the case of story ductility demands, RµS is calculated for several levels of structural deformation, the value of this parameter associated with a deformation state close to that of a collapse mechanism is that considered in the seismic guidelines [69]. The results are given in Table 4 and in Figures 4(c) and 4(d), for the 3- and 10-level models, respectively. For the 3-level building, no trend is observed between the mean values and the story number. For the case of the maximum deformation, the mean values are 3.10 and 3.06, for the NS and EW directions, respectively. As for the case of local and story ductility, the values are larger for RµLϕ than for RµS and, however, the differences are larger for the case of ductility demands. For the 10-level building, unlike the 3-story building, the mean values of RµS, even though not in a perfect way, tend to increase with the story number. The values seem to be larger for the 10-level building than that for the 3-level building, values larger than 4 are observed in some cases.

4.2.3. Global Ductility Reduction Factor

The mean values of global ductility reduction factors (RµG) are presented in Table 5 for all cases. The results indicate that, for maximum deformation, RµG are quite similar for both structural models, and the values range from 2.87 to 3.09. It is implicitly assumed in the seismic codes mentioned in Section 1 that the magnitude of the global ductility reduction factor is about 4. According to the results found in this paper, this value is not justified; a value of 3 seems to be more reasonable.

4.2.4. Ratio of Ductility Reduction Factor to Ductility

The ratios of the ductility reduction factor to ductility, for local (QL), story (QS), and global (QG) structural levels, are presented in this section. They are expressed as

The values of QL are given in last rows of Table 3 and in Figure 5, for the 3- and 10-level buildings, respectively. The results indicate that this ratio significantly varies with the type of the structural element (beam or column), story number, strong motion intensity, and building height. For beams of the two models, the QL values tend to decrease with the strong motion intensity, but they tend to increase with the story number for the 3-level building. For the case or maximum deformation, QL is quite similar for the two models; the values range from 0.32 to 0.6. Most of the observations made for beams apply to the case of columns, and however, an important difference is that the QL values are significantly larger for columns; the values of up to 2.4 and 3.0 are observed for the 3- and 10-level buildings, respectively. This result has an important implication: while the well-known ratio between ductility reduction factor to ductility [10, 33] indicates that, for the structural models under consideration (fundamental lateral vibration periods of 1.03 and 2.41), the maximum reduction in the elastic forces should be equal to ductility capacity, and the results in this paper show that the reduction for local response parameters like bending moment in columns can be much larger than that. The implication of this is that the design of columns will be too conservative if the reduction is performed according to the mentioned well-known ratio.

The results for QS are presented in the last rows of Table 4 and in Figure 6, for the 3- and 10-level buildings, respectively, while those of QG are given in the last row of Table 5. The QS values resemble those of QL in the sense that they decrease with the strong motion intensity and that they are quite similar for the two buildings. However, they can be significantly larger for the case of beams; for the maximum deformation, they range from 0.68 to 0.86 and from 0.60 to 1.10 for the 3- and the 10-level buildings, respectively. Regarding the QG parameter, the results in Table 5 indicate that, as for the local and story ratios, QG decreases with the intensity of the earthquakes; for the level of deformation associated with ductility capacity (Sa = 1.2 g and 0.6 g), it takes values of 0.73, 0.76, 0.72, and 0.86. Thus, from an overall point of view, the ratio of the global ductility reduction factor to global ductility is about 3/4, and it is the value proposed in this study. This value is smaller than the value of 1.00 specified by the well-known ratio [10, 33] abovementioned. The implication of this is that if the value of 1.00 is used to calculate the global ductility reduction factor, from a global point of view, nonconservative designs may result.

4.3. Objective 3: Dissipated Energy
4.3.1. Local, Story, and Global Normalized Dissipated Energy per Joint

As for the ductility demand and ductility reduction factor parameters, the dissipated energy is calculated for three levels. The mean values of (ELϕ)SB, SC, ES, and EG are calculated according to Equations (10), (13)–(15), respectively. Plots for individual strong motions were developed for each beam, column, story, and the whole structure; however, only the results in terms of the mean values, averaged over all the strong motions, are given below. For the case of the 3-level building, the mean values of (ELϕ)SB and SC are given in Table 6, while those of ES and EG are given in Tables 7 and 8, respectively. The corresponding results for the 10-level building are given in Figures 7 and 8 and in Table 8. In Table 8, the symbols “B+C” or “B” indicate that EG was calculated considering beam and columns, or only beams, respectively.

Results for the 3-level building indicate that (ELϕ)SB, SC, and ES tend to decrease with height and that (ELϕ)SB is much larger than SC, as expected. For the 10-level building, (ELϕ)SB andSC tend to increase and decrease with the story number, respectively; however, ES does not show a clear tendency. From a global point of view (Table 8), it is shown that the values of the global normalized dissipated energy (EG), even for the maximum deformation, are relatively small when beams and columns are considered; the maximum observed values are 0.129 and 0.218 for the 3- and 10-level buildings, respectively. These relatively small values of EG are due to the fact that the normalized energy demands in columns are quite small so the average values obtained from Equations (14) and (15) for the whole frame are expected to be small too. The EG parameter is also calculated by assuming that the plasticization of the steel frames, perfectly follows the strong column-weak beam concept; in other words, the dissipation of energy in columns is not considered. It is observed from Table 8 in this case that the maximum values of EG are now 0.324 and 0.475, for the 3- and 10-level models, respectively.

4.3.2. Ratio of Ductility Reduction Factor to Normalized Dissipated Energy

It is clear that the reduction of the response from the elastic to the inelastic case due to yielding of the material, which is quantified by the ductility reduction factor, is produced by the dissipated energy. In this regard, it may be of interest to study the ratio of the ductility reduction factor to dissipated energy. In this section of the paper, the ratio of the story ductility reduction factor to story normalized dissipated energy (PS) and the ratio of global ductility reduction factor to global normalized dissipated energy (PG) are given and discussed. They are calculated as

The results for PS are given in Table 9 and in Figure 9, for the 3- and 10-level models, respectively. The results for PG are given in Table 8 for both models. In the calculation of ES, only the dissipated energy at beams is considered. A value of PS = RµS, or PG = RµG, would mean that ES = 1.0, or EG = 1.0, which in turn would correspond to the idealized hypothetical case where all structural elements simultaneously reach, for every strong motion, the energy capacity. In other words, the energy demand equals the energy capacity in all cases. Obviously, it is not the case and, as shown below, values much larger than RµS or RµG are obtained for PS or PG, respectively.

It is observed from Table 9 and Figure 9 that the PS values significantly vary with the model height, the story number, and the strong motion intensity. Since energy dissipation produces the response reduction, one could expect a similar variation of RµS and ES which would approximately imply a constant value for PS; however, these two parameters vary in a different proportion. The results clearly indicate that PS tends to significantly increase with the story number for the 3-level building; this trend is not observed for the 10-level building. For both buildings, the PS values significantly decrease with the strong motion intensity in most of the cases. One reason for this is that, as stated before, the seismic response significantly varies from one ground motion to another even though the level of deformation, in terms of Sa, is approximately the same for any of the strong motions. This variation is observed even for two strong motions with a similar predominant period. These differences are due, in part, to the inherent variability existing in ground motions characterized by their frequency contents; it causes that the response be very sensitive to the particular characteristics of a given story or of a given building model and to the ground motion under consideration. This sensitivity is quite different from one response parameter to another and seems to be much larger for dissipation of energy than for the ductility reduction factor.

Results in Table 8 indicate that some of the observations made before for PS apply the PG Ratio. Based on the results associated with the maximum structural deformation, a value of 8 is proposed for this ratio. Thus, if the global normalized dissipated energy is estimated for a steel building with the structural system under consideration, the global ductility reduction factor can be estimated as eight times EG.

5. Conclusions

The ductility parameter plays a central role in seismic analysis and design of steel buildings. However, it is still used in an indirect way. A numerical investigation regarding the evaluation of local (µLϕ), story (µS), and global (µG) ductility, local (RµLϕ), story (RµS), and global (RµG) ductility reduction factors, as well as local (ELϕ), story (ES), and global (EG) normalized dissipated energy, for steel buildings with moment-resisting frames, which were designed following the strong column-weak beam concept, is presented in this paper. The ratios of µG to µLϕ (Q), RµLϕ to µLϕ (QL), RµS to µS (QS), RµG to µG (QG), RµS to ES (PS), and RµG to EG (PG) are also calculated. 3- and 10-level buildings and some strong motions used in the SAC Steel Project are used in the study. The mentioned parameters are calculated for several intensities of the strong motions. Results of the study indicate that any of the abovementioned parameters significantly may vary with the strong motion, the seismic intensity, the structural element, the story number, and the structural model. The main findings are as follows:(1)Bending local ductility capacity (µLϕ) of beams can reach values of up to 20 for some individual strong motions, as observed in some experimental investigations. On an average basis, µLϕ ranges from 7 to 14. The values are larger for the 10-level building than for the 3-level building, reflecting the effect of the structural complexity on the µLϕ parameter. Many of the observations made for µLϕ apply to µS, but the values are much smaller than those of µLϕ; they range from 3 to 5. The corresponding values of µG range from 3.6 to 3.95. A value of 1/3 is proposed for the Q ratio. Thus, if local ductility capacity is stated as the basis for the design, global ductility capacity can be calculated by using this ratio.(2)The mean values of RµLϕ for beams range between 3.5 and 5, which are larger than those of columns. Even though yielding did not occur in columns for many cases, a considerable reduction is observed in bending moments implying that yielding on beams reduce not only their bending moments but also those of columns. The maximum mean values of RµS is about 4.0. The RµG mean values range from 2.87 to 3.09. It is implicitly assumed in the seismic codes that the magnitude of the global ductility reduction factor is about 4. According to the results found in this paper, this value is not justified; a value of 3 seems to be more reasonable and is proposed in this study.(3)The maximum values of QL range from 0.32 to 0.6 for beams; for columns, however, they are significantly larger; values larger than 2 are observed in many cases. The maximum values of QS range from 0.60 to 1.10 while those of QG from 0.73 to 0.86. These values are smaller than the specified well-known ratio of 1.00. The implication of this is that if the value of 1.00 is used to calculate the global ductility reduction factor, from a global point of view nonconservative designs may result, while from a local point of view, as is the case of columns, very conservative designs may result. A value of QG of 3/4 is proposed in this study.(4)The ELϕ values for beams for maximum deformation range from 0.15 to 0.7, which are much larger than those of columns and larger than those of ES. The maximum values of EG range from 0.25 to 0.47. The values of the PS ratio significantly decrease with the strong motion intensity and significantly increase with the story number, in most of the cases. Based on the results associated with the maximum deformation, a value of 8 is proposed for the PG ratio.

Data Availability

Data used to support the findings of this study are included within the article, and data supporting this study are from previously reported studies, which have been cited.

Disclosure

Any opinions, findings, conclusions, or recommendations expressed in this publication are those of the authors and do not necessarily reflect the views of the sponsors.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

The research presented in this paper was financially supported by La Universidad Autónoma de Sinaloa (UAS) under grant PROFAPI-2015/235.