Next Article in Journal
A Novel Heart Rate Robust Method for Short-Term Electrocardiogram Biometric Identification
Previous Article in Journal
Capacity Credit Evaluation of Correlated Wind Resources Using Vine Copula and Improved Importance Sampling
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Local, Story, and Global Ductility Evaluation for Complex 2D Steel Buildings: Pushover and Dynamic Analysis

by
Mario D. Llanes-Tizoc
,
Alfredo Reyes-Salazar
*,
Eden Bojorquez
,
Juan Bojorquez
,
Arturo Lopez-Barraza
,
J. Luz Rivera-Salas
and
Jose R. Gaxiola-Camacho
Facultad de Ingeniería, Universidad Autónoma de Sinaloa, Culiacán Rosales, Sinaloa 80013, Mexico
*
Author to whom correspondence should be addressed.
Appl. Sci. 2019, 9(1), 200; https://doi.org/10.3390/app9010200
Submission received: 17 November 2018 / Revised: 31 December 2018 / Accepted: 2 January 2019 / Published: 8 January 2019

Abstract

:
A numerical investigation regarding ductility evaluation of steel buildings with moment resisting steel frames is conducted. Bending (µ) and tension (µ) local ductilities as well as story (µS) and global ductilities are studied. Global ductility is calculated as the mean values of story ductilities (µGS) and as the ratio of the maximum inelastic to yielding top displacements (µGt). The ductility capacity is associated to drifts of about 5%. Ductility values significantly may vary with the strong motion, ductility definition, structural element, story number, type of analysis, and model. µ is much larger for beams than for columns. Even though the demands of µ are considered an important issue they are less relevant than µ. µS is much smaller than µ for beams. µGS for dynamic analysis give reasonable values, but µGt does not. µ, µS and µGS obtained from pushover are larger than those obtained from dynamic analysis and unlike the case of dynamic analysis, µ tend to increase with the story number showing an opposite trend. Considering that: µGt for dynamic analysis results in unreasonable values, pushover analysis does not consider energy dissipation, the strong column–weak beam (SCWB) concept was followed in the model designs, and µ is not relevant in framed steel buildings, the ratio (RLG) of global to local ductility capacity is calculated as the ratio of µGS to µ of beams, for dynamic analysis. A value of 1/3 is proposed. Thus, if bending local ductility capacity is stated as the basis for the design, the global ductility capacity can be easily estimated.

1. Introduction

The seismic behavior of a building under the action of a strong motion represents a very complicated phenomenon, particularly when the building is deformed into the inelastic range. In spite of this most of building codes around the world permit the use of simple elastic procedures to determine the seismic demands on steel buildings either for small or large deformations. Due to their relatively simplicity in their application, simplified methods like the Static Equivalent Lateral Force (SELF) procedure, are broadly used. For example, The International Building Code (IBC) [1], The National Building Code of Canada (NBCC) [2], The Mexico Federal District Code (MFDC) [3], and The Eurocode 8 (EC) [4], permit the use of the mentioned procedure for regular buildings with relatively short periods (low- and medium-rise). FEMA-273 and ATC-40 also permit to use nonlinear static procedure or pushover analysis. In the procedure, static analysis of the buildings under the action of equivalent lateral forces, which are related to the properties of the structure and the seismicity of the region, provides the design forces; then a serviceability revision is performed. Thus, it can be said that conventional seismic design considered in many seismic codes is essentially force-based with a final check on displacements.
In the above mentioned procedure, the ductility parameter (µ) plays an important role in the determination of the design seismic forces of building structural systems, it represents the capacity of a structure to dissipate energy, allowing for a reduction of the elastic strength demands; the larger the ductility, the smaller the design seismic forces. It is particularly important for steel structures since there are many sources of ductility and of energy dissipation. 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 or unanimity on the profession on how to define it [5,6]. In fact, the magnitude of the reduction of the elastic design seismic forces directly depends on the force reduction factor (R), which significant depends on an associated parameter called ductility reduction factor (Rµ) [7,8]. The estimation of these parameters represents one of the most controversial issues in the SELF procedure.
In IBC the reduction factor is called ‘response modification factor’ (R); it is stated that it depends on many parameters including the ductility capacity and inelastic performance of structural material. In NBCC the reduction parameter is called ‘force modification factor’ (Rd), where it is explicitly stated that it depends on the ductility capacity and materials as well as on the structural overstrength. In MFDC this parameter is called ‘seismic reduction factor’ (Q′) which depends on the structural material, the structural system and detailing. In EC the factor is called ‘seismic reduction factor’ (q′) and depends on an overstrength factor and on factors which in turn depends on structural materials and the structural system. It is implicitly assumed in these codes that ductility represents the capacity of the structure to dissipate energy.
It is concluded that, it is essential to establish a measure of ductility. In this regard, the various types of ductility involved in a building must be considered [9]; local ductility must be differentiated from story or from global ductility. For the case of steel buildings, local ductility may be associated to the rotational capacity of a member under bending or to the longitudinal deformation capacity of a member under tension. Story ductility is essentially associated to the relative story displacements while global ductility is normally expressed in terms of story ductility or in terms of the absolute displacements of the roof. It is generally accepted than the 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 is the ratio of the maximum interstory lateral displacement of a structure during the application of the seismic loading to the corresponding displacement when first yielding occurs at any member of the story, while ductility capacity is the ratio of the maximum permissible inelastic lateral displacement to the displacement when first yield occurs. Ductility capacity is usually obtained from experimental results of individual members (local ductility). Therefore, it is important to relate it 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. Moreover, local ductility cannot be exceeded; otherwise, the ductility corresponding to the collapse mechanism will not be the ductility capacity. For that reason, some researchers [10] suggest using local ductility as the basis for design because there are numerous laboratory studies on ductility for members. In this regard, as stated above, some relationships need to be established between local, story and global ductility.
The central objective of this paper is to evaluate the ductility parameter for steel buildings with typical welded moment-resisting frames. Different types of local ductility as well as story and global ductility are calculated according to both, nonlinear dynamic analysis and nonlinear static analysis (pushover). The relationship between local and global ductility is calculated. Due to the advancements in the computer technology, the computational capabilities have significantly increased in the recent years allowing us to estimate the nonlinear seismic response by modeling structures as complex MDOF systems with hundred and even thousands of degrees of freedoms and applying the seismic loadings in time domain as realistically as possible. Responses obtained in this way may represent the best estimate of the seismic responses. Then the values of ductility demands for the different definitions can be properly estimated.

2. Literature Review

Studying the µ parameter and the associated ductility reduction factor (Rμ for the case of the IBC code) of steel buildings has been an important research topic during the last decades. There have been a significant number of studies for single degree of freedom (SDOF) systems based on analytical and/or empirical observations and considerations. The Rμ factor was first introduced in ATC-3-06 [11] in the late 1970s. Other of the first investigations was conducted by Newmark and Hall [9]; they proposed a procedure to relate Rμ and µ by constructing the inelastic response spectra from the basic elastic design spectra. Hadjian [12] studied the reduction of the spectral accelerations to account for the inelastic behavior of structures. Miranda and Bertero [13] proposed simplified expressions to estimate the inelastic design spectra. Tiwari and Gupta [14] proposed a preliminary scaling model to estimate the ductility reduction factors of horizontal ground motions. Significant contributions regarding the evaluation of the ductility and ductility reduction parameters for SDOF systems can be found in other publications [15,16,17,18,19,20]. More recently Zhai et al. [21] investigated the strength reduction factor of SDOF systems with constant ductility performance subjected to the mainshock–aftershock sequence-type ground motions. Ghods et al. [22], by using the finite element method (FEM), investigated the forming of plastic hinges, distribution of stresses, and ductility and stiffness of steel systems composed of reinforced concrete column to steel beam connections.
There are also several studies regarding the evaluation of R (or Rµ) and µ factors for multi degree of freedom (MDOF) systems. Nassar and Krawinkler [23] studied the relationship between force reduction factors and ductility for simplified (three-story single-bay) MDOF systems. Santa-Ana and Miranda [24] studied the strength reductions factors for several steel frames modeled as plane MDOF systems. Moghaddam and Mohammadi [25] 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 [26] investigated the relationship between the lateral capacity, the design force reduction factor, the ductility factor and the overstrength factor for reinforced-concrete buildings. Medina and Krawinkler [27] presented an evaluation on drift demands for regular moment resisting frame structures subjected to ordinary ground motions. In another study Medina and Krawinkler [28] studied the strength demands relevant for the seismic design of moment-resisting frames. Important results, regarding ductility, ductility reduction factor and other related parameters for structures modeled as MDOF systems can be found in some research reports and papers [29,30,31,32,33,34,35,36,37,38].
More recently Reyes-Salazar et al. [39], 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. Valenzuela-Beltran et al. [40] proposed a reliability-based criterion including two simplified mathematical expressions, which depends on the ductility of the structural system, to estimate strength amplification factors for buildings with asymmetric yielding. Fanaie and Shamlou [41] studied the seismic behavior, in terms of response modification factor and ductility factor, of mixed structures. Vuran and Aydınoğlu [42] developed simple capacity and ductility demand estimation tools for coupled core wall systems. Gómez-Martínez et al. [43] analytically studied the local and global ductility of wide-beam reinforced concrete moment resisting frames. Wang et al. [44] studied the seismic performances of steel braced truss-RC column hybrid structure. Liu et al. [45] 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. Kang and Mory [47] proposed simplified procedure to estimate the peak inter-story drift ratios of steel frames with hysteretic dampers for SDOF and MDOF systems where the energy dissipated by hysteretic behavior and the involved ductility were explicitly considered.
The abovementioned studies represent a significant contribution regarding the evaluation of ductility or force reduction factors, however, in most of them SDOF systems, plane shear buildings, or a limited level of inelastic deformation were considered. Therefore, they did not explicitly consider the inelastic behavior and energy dissipation of the structural elements existing in actual systems. It has been shown [48,49,50] that ductility demands as well as the ductility reduction factors depend on the amount of dissipated energy, which in turn depends on the plastic mechanism 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 to the ductility capacity. Moreover, local ductility taking into account the maximum inelastic curvature and tensional elongations as well as relationship between local and global ductility, estimated by dynamic and pushover analysis, have not been considered.

3. Objectives

The central objective of this paper is to evaluate ductility demands of steel buildings; local, story and global ductility as well as relationship among them are considered. To this aim, the nonlinear seismic responses of low- and medium-rise steel building models with typical welded connections, idealized as complex plane MDOF systems, are calculated. The ductility values are compared with those resulting from pushover analysis. The specific objectives are:
(1)
Calculate local ductility for individual structural elements (beams and columns) associated to the maximum inelastic curvatures and axial elongations according to nonlinear seismic analysis. Several intensities of the strong motions are considered. The larger intensity of the earthquakes will correspond to a deformation state close to a structural collapse mechanism (maximum drifts of about 5%).
(2)
Calculate story and global ductility for the same cases above mentioned.
(3)
Compare the ductility values obtained from nonlinear dynamic analysis with those of pushover analysis.
(4)
Propose a relationship between local and global ductility.

4. Procedure and Structural Models

Two steel buildings, modeled as MDOF plane frames, and twenty strong seismic motions are used in the study. Local, story and global ductility, according to several definitions, are calculated. The strong motions are scaled down and up to get different levels of structural deformations. The maximum deformation level is very close to the formation of a collapse mechanism; therefore, it is associated to the ductility capacity. In some experimental studies 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 stable manner [51,52,53,54]. This is approximately the maximum deformation level considered in this study. Thus, the ductility values obtained for this deformation, according to the different definitions, are assumed to be the ductility capacity.

4.1. Steel Building Models

Several steel model buildings were designed, as part of the SAC (Structural Engineers of California, Applied Technology Council, and Consortium of Universities for Research in Earthquake Engineering) steel project [55], by three consulting firms. The buildings are supposed to satisfy all code requirements existed at the time of evaluation for the following three cities: Los Angeles [56], Seattle [56], and Boston [57]. The buildings were designed to be case studies to evaluate many aspects related to the seismic behavior of steel buildings with moment resisting frames [55]. Even though the three sets of buildings were used in the SAC project only the perimeter moment resisting frames (PMRF) of the 3- and 10-level buildings of the project, which are assumed to be located in the Los Angeles area, are considered in this study to address the issues discussed earlier. Isometric views for the 3- and 10-level buildings, which are referred hereafter as Models 1 and 2, respectively, are given in Figure 1 and Figure 2, where the PMRF (exterior) and the interior gravity frames are clearly identified. Sizes of structural elements, namely beams (horizontal) and columns (vertical), as reported, are given in Table 1 for the two models. The first three translational periods of the plane frames associated to 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 building models are assumed to have typical welded connections. The structural members, which were modeled as one-dimensional beam-column elements, were designed following the strong column–weak beam (SCWB) concept. Each column is represented by one element and each girder by two elements, having a node at the mid-span. Each node is considered to have three degrees of freedom. The damping is considered to be 3% of critical damping. The RUAUMOKO Computer Program [58] is used for the time history nonlinear dynamic analysis. The frames are modeled as complex 2D MDOF systems. The Newmark constant average acceleration method together with the Newton–Raphson iteration scheme are used to solve the differential equation systems. For a given load increment convergence occurs (it is the iterative process is terminated) when the ratio of the incremental displacements at the current iteration to the current estimate of the displacements is small enough (about 0.00001). In addition lumped mass matrix, Rayleigh damping as well as large displacement effects, are considered in the nonlinear dynamics 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, as well as elastic and inelastic section properties are given or read within the Ruaumoko computer program environment. No strength degradation member, bilinear behavior with 5% of the initial stiffness in the second region and concentrated plasticity are assumed. The interaction axial load-bending moment is given by the yield interaction surface proposed by Chen and Atsuta [59].
It is worth mentioning that the numerical investigation was developed on the base of the above-described simplified FE models, where the mechanical, geometrical, and dynamics properties of structural elements are roughly, but efficiently, described. It has been shown [60,61,62,63] that this type of FE formulation results in a good approximation of the structural response in parametric studies like that presented in this paper, as long as it is provided a realistic modeling of the structure and of the cyclic load deformation characteristics of its structural elements.

4.2. Earthquake Loading

The responses of a structure excited by different strong motions, even when they are normalized with respect to the same response parameter, are expected to be different reflecting the influence of their different frequency content. Thus to get meaningful results, the models under consideration are excited by twenty strong motions in time domain. The characteristics of the records are given in Table 2 for the NS direction.
The predominant periods of the earthquakes vary from 0.11 to 0.62 s, which are defined as the period where the largest peak in the pseudo-acceleration (Sa) elastic response spectrum occurs. The earthquake time histories were obtained from the Data Sets of the National Strong Motion Program (NSMP) of the United States Geological Surveys (USGS). One horizontal seismic component at a time as well as the vertical component and the gravity loads are applied to the models. For a given direction, half of the seismic loading and gravity loading are applied to the corresponding PMRF.
The models behave essentially elastic under the action of any of the strong motions. In order to have different levels of deformation as well as moderate and significant inelastic behavior, the strong motions are scaled in terms of Sa evaluated in 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.8 g for the 10-level model, with increments of 0.2 g. The maximum considered values (Sa = 1.2 g and Sa = 0.8 g) correspond to a deformation state very close to a collapse mechanism where, as stated above, interstory drifts of about 5% were developed for many strong motions.

4.3. Gravity Loading

In addition to the seismic loading, the following gravity loads [55] were used in the analysis: (a) the floor dead load for weight calculations was 4.5 kN/m2; (b) the floor dead load for mass calculations was 4.04 5 kN/m2; (c) the roof dead load was 3.9 kN/m2; (d) the reduced live load per floor and for roof was 0.94 kN/m2. The seismic mass for the entire structure was: (a) for the roof of the 3-level building it was 1023.09 kN-s2/m; (b) for floors 2 of the 3-level building it was 945.6 kN-s2/m; (c) for the roof of the 10-level building it was 1054.83 kN-s2/m; (d) for floor 2 of the 10-level building it was 996.25 kN-s2/m; (e) for floors 3 to 9 of the 10-level building it was 979.22 kN-s2/m.

5. Ductility Definitions

The main objective of this paper is to evaluate the ductility parameter (local, story, and global) associated to several levels of deformations of steel buildings under the action of seismic loading. From a conceptual point of view, in the context of seismic analysis of SDOF systems, 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 total load 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 calculate these two parameters (Dmax and Dy); many alternatives are possible. The implication of this is that there is no unanimity in the profession on how to calculate the ductility parameter for MDOF systems. As it will be additionally discussed below, global ductility for MDOF systems, typically is expressed as the ratio of the maximum absolute lateral displacement of the roof after the complete application of the loading to the absolute lateral displacements of the roof when yielding at any member occurs by the first time. However, since global ductility should represent the overall structural inelastic deformation, some researchers suggest defining it in terms of relative lateral displacements [5,7,54]. In addition, it is necessary to properly define local, story and global ductility and establish appropriate relationships among these parameters. In this paper the following definitions of ductility are considered.

5.1. Local Ductility

Definition 1.
The local ductility of a flexural member (µ) for a given joint is defined 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
μ L ϕ = ϕ m a x ϕ y
Thus in the case of nonlinear time history analysis, 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 member under consideration.
Definition 2.
The local ductility of a tensile member (µ) is defined as the ratio of the maximum inelastic axial deformation that the member undergoes during the total time of excitation (δmax) to the axial deformation of the member when it yields for the first time (δy). It is expressed as
μ L δ = δ m a x δ y

5.2. Story Ductility

Definition 3.
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
μ S = Δ m a x Δ y
The Δy parameter in Equation (3) needs additional discussion at this stage: it is assumed that for a given story of a given frame, the beams on the story and the columns connecting beneath it are part of the story. For example, for the 3-story model, going from the top to the bottom (see PMRF in Figure 1), 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.

5.3. Global Ductility

Definition 4.
Global ductility is defined as the mean value of the story ductilities; it is calculated as
μ G S = 1 n i = 1 n ( μ S ) i
where n is the number of stories.
Definition 5.
Global ductility is calculated as the ratio of the maximum inelastic displacement at the top of the building during the total time of excitation (Dmax,t) to the top displacement when any member of the building yields for the first time (Dy,t). It can be expressed as:
μ G t = D max , t D y , t

6. Objective 1: Local Ductility, Dynamic Analysis

The local ductility parameter associated to bending, as defined by Equation (1), is calculated for all the structural members of both models for several intensities of the 20 strong motions as well as for the NS and EW directions. First, the models are subjected to the simultaneous action of the horizontal seismic component oriented in the EW direction, the vertical seismic component and the gravity loads. Then, the models are subjected to a similar set of loads, but the other horizontal component (NS direction) is considered instead. It is important to mention that even though, as stated earlier, the strong column–weak beam concept was considered in building design calculations, column hinging occurred in many cases.
For a given story, the µ values are averaged, first over all the beams, and then over all the columns. The resulting averages for the beams of the 3-level building and the EW horizontal component, for seismic intensities Sa = 0.4 g, 0.8 g and 1.2 g, are presented in Figure 3a–c, respectively, while those associated to columns are given in Figure 3d–f, for the same strong motion intensities. The corresponding results for the 10-level building are presented in Figure 4, but in this case for seismic intensities of 0.2 g, 0.6 g and 0.8 g. In these figures, the symbol “ST” stands for the story level. It is worth to mention that moderate yielding occurs 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 deformations are 1.2 g and 0.8 g, for the 3- and 10-level models, respectively. In fact, even though it is not shown in the paper, the drifts for these significant levels of deformation were about 5% for some of the strong motions; the corresponding pattern of plastic hinges were very close to define a collapse mechanism. For this reason, the mentioned levels of maximum deformation are assumed to be associated the structural capacity of the models.
Results of the figures indicate that, for beams, the maximum values of µ (bending local ductility capacity) are about 15 and 20 for the 3- and 10-level buildings, respectively. There are some numerical and experimental [9,64,65] evidence that the bending local ductility capacity can reach values larger than 20; however, it is for monotonic loading and individual members. It is also observed that, for a given value of Sa, the magnitude of µ significantly varies from one seismic motion to another even though the seismic demands normalized according to Sa for each seismic motion was the same.
For example, for the case of beams of Story 3 of the 3-level building and Sa = 1.2 g, values as small as 4 and as larger as 15 are observed for the twentieth and twelfth strong motions, respectively. Such broad variation reflects the influence of the strong motion frequencies and the contribution of several vibration modes to the structural response. Results also indicate that the µ values increase as the seismic intensity increases and that the values are much larger for beams than for columns, as expected. Plots similar to those given in of Figure 3 and Figure 4 were also developed for the NS direction. In total 20 and 16 plots were developed for the 3- and 10-level buildings, respectively, but only the fundamental statistics in terms of the mean values (MV) and coefficients of variation (CV) are given for all cases; the results are presented in Table 3 and Table 4.
It can be observed from the tables 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, as observed from individual graphs, 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. Unlike what observed for beams, the mean values of µ for columns tend to decrease with the story number. For both, beams and columns, the mean values tend to increase with the seismic intensity and are larger for the NS that for the EW direction; the uncertainty in the estimation is moderate in most of the cases.
The statistics for the 10-level building (Table 4) resemble those of the 3-level building in the sense that the mean values of µ are much larger for beams than for columns and that the uncertainty in the estimation is, in general, moderate. 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 µ mean values tend to decrease with the story number for columns; no yielding occurs in columns in most of the cases for low 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 are observed to be larger for the 10-level building.
Local ductility values associated to axial elongations (µ) were also calculated. It is worth to mention that since axial deformations of beams are negligible in comparison with those of bending, µ is not calculated for beams. In other words, plasticization of beams (formation of plastic hinges) is essentially produced by bending. In addition, it was observed that, due to their location, the axial loads at interior columns are smaller than that of exterior columns in such a way that plasticization of interior columns are essentially due to the action of bending moments too. Thus, that the µ parameter is calculated only for exterior columns
Similar to the case of µ, for a given story and building, the µ values are averaged (but in this case only over the two exterior columns) and graphs for individual members are developed; only the fundamental statistics for the 3-level building are presented (Table 5), however. The main observations that can be made are that the mean values of µ are comparable to those of µ for columns but much smaller than those of µ for beams. However, unlike the case of µ for beams or columns, the maximum values occur for the upper story.
The previous discussion clearly indicates that the ductility demands associated to column axial deformations of steel framed structures are much smaller than those associated to bending of beams; for example, for the 3-level building and the deformation state close to collapse (Sa = 1.2 g), the average ductility demand associated to bending of beams ranges from 7.8 to 11.86, while that associated to axial elongation on columns ranges from 1.00 to 3.57. Thus, even though tensional ductility demands are considered as an important issue in some research reports [9] and steel structural members may have a high ductility capacity in tension, these types of ductility demands are less relevant for the type of the structural system under consideration.

7. Objective 2: Story and Global Ductility, Dynamic Analysis

7.1. Story Ductility

The µS values for the NS direction of the 3-level building, calculated according to Equation (3), are shown in Figure 5a–c, for seismic intensities of 0.4 g, 0.8 g, and 1.2 g, respectively; while those of the 10-level building for the same direction are given in Figure 5d–f, for seismic intensities of 0.2 g, 0.6 g, and 0.8 g, respectively. These results resemble those of µ in the sense that they do not show any trend with the strong motions; the additional observation that can be made is that bending local ductility values are, in general, larger than those of story ductility. The fundamental statistics, for all the seismic intensities and structural directions under consideration, are presented in Table 6 and Table 7, for the 3- and 10-level buildings, respectively. It can be observed that, as for the case of µ, 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, unlike the case of µ the values do not tend to decrease with high.
By comparing the results of Table 6 and Table 7 with those of Table 4 and Table 5, it is concluded that, as observed from individual plots, µS is much smaller than µ of beams; the 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. As for µ for beams, the mean values and the uncertainty in the estimation tend to increase with the seismic intensity, however the uncertainty in the estimation is larger for µS.

7.2. Global Ductility

The global ductility values, calculated according to the two definitions under consideration (µGS and µGt), are now discussed. The results are given in Table 8 and Table 9, for the 3- and 10-level buildings, respectively. It can be observed that the values of µGS tend to increase (as expected) with the strong motion intensity and that they are larger for the 3- than for the 10-level building; for the ductility demands (ductility capacity) associated to the maximum deformations µGS takes values of 3.95 and 3.89, for the NS and EW direction, respectively, for the 3-level building, while the corresponding values are 3.82 and 3.31 for the 10-level building. By comparing the results of Table 8 and Table 9 with those of µ for beams (Table 3 and Table 4) it is noted that, as for the case µS, the values are significantly smaller for µGS. Most of the observations made µGS also apply to µGt; however, unreasonable large values are observed for the latter, values of 36.34 and 38.78 were obtained for the 10-level buildings. The reason for this is not that very large drifts were obtained, but the top displacements when first yielding occurs in the buildings were very small indicating an important participation of the higher modes in the response. The implication of this is that the ductility definition based on the building top displacement, even though used in the profession, is not appropriate for the case of dynamic analysis, particularly for high buildings where significant structural contributions of the higher modes are expected.

8. Objective 3: Ductility According to Pushover Analysis

It is accepted that nonlinear time history analysis is the most accurate and reliable analysis procedure as long as realistic modeling of the structure and the cyclic load deformation characteristics of its structural elements are provided. The latter characteristic cannot be considered in pushover analysis implying that the effect of the dissipation of energy on the seismic response is neglected. Despite this, nonlinear static procedures are broadly used to estimate seismic responses in terms of different parameter for low- and medium-rise buildings. In this section of the paper, the different types of ductility under consideration are calculated by using pushover analysis and compared to those of dynamic analysis. The pattern of loads is the classical one where a triangular distribution is used; gravity loads were also considered. In order to have a reasonable comparison, the maximum drifts considered in pushover analysis (about 5%) is quite similar to those observed in dynamic analysis for many strong motions for the case of maximum seismic intensity (Sa = 1.2 g and 0.8 g for the 3- and 10-level buildings).
Similar to the case of dynamic analysis, as the lateral loads are gradually increased, the structural deformations are gradually increased too. As soon the first plastic hinge is developed in a particular member, in any member of a story, or in any member of the structure, the corresponding deformations ϕy, Δy and Dy,t as indicated in Equations (1), (3), and (5), respectively, are recorded. Then, by considering the associated maximum deformations (ϕmax, Δmax and Dmax,t) the values of local, story and global ductility are calculated.
The results for bending local ductility are shown in last column of Table 3 for the 3-level building. As for the case of dynamic results, the values are much larger for beams than for columns. It is observed that for beams, the values of µ of pushover analysis, which range from 13.01 to 16.64, can be much larger than the µ mean values of dynamic analysis (ranging from 7.81 to 11.86). There are several reasons for this: (a) as stated earlier, the maximum drift was approximately the same for both type of analyses, however, it was observed for many but not for all the strong motions in the case of dynamic analysis, for many of them the maximum drift were about 4% for the maximum seismic intensity, (b) for the case of pushover, the drifts were very close to 5% for all stories, but for dynamic analysis values of 3% and even smaller were observed for some stories, (c) the drift at which first yielding occurred for the stories were smaller for pushover than for of dynamic analysis. Results also indicate that for the case of columns, µ values may be larger or smaller than those of dynamic analysis; however, the values tend to increase and decrease with the story number for pushover and dynamic analysis, respectively, showing an opposite trend.
The µ pushover results for the 10-level building are shown in last column of Table 4. The major observations made before for the 3-level building apply to this building: (a) the pushover values are larger for beams than for columns (in both cases tend to increase with the story number), (b) for the case of beams, the maximum pushover values are larger than the maximum dynamic mean values, and (c) for beams and columns, the values tend to increase and decrease with the story number for pushover and dynamic analysis, respectively. The additional observation that can be made is that for the case of beams, the maximum pushover or dynamic values are larger for the 10- than for the 3-level buildings. The µS values for pushover analysis are given in last columns of Table 6 and Table 7 for the 3- and 10-level models, respectively. It can be observed that, as for bending ductility demands of beams, µS in general is larger for pushover than for dynamic analysis. It is also shown that the values are larger for the 3- than for the 10-level building and that for the 3-level building the values tend to slightly decrease with the story number which is not observed for the 10-level building.
The global ductility values, as for the case of dynamic analysis, are calculated in terms of interstory displacements (µGS) and of top displacements (µGt). The results are shown in last columns of Table 8 and Table 9. As for µS, µGS is larger for pushover (7.72 and 5.08) than for dynamic analysis (3.95 and 3.32) for both models. µGt is larger for pushover than for dynamic analysis for the 3-level building, but unlike the case of dynamic analysis reasonable large values are observed for the case of the 10-level building.

9. Objective 4: Relationship between Local and Global Ductility

In experimental studies the ductility capacity is usually obtained for individual members (local ductility), for this reason, as stated earlier, it is suggested [10] considering local ductility as the basis for design. Thus it results convenient relate local to global ductility. In this section the ratio of global to local ductility, denoted by the RLG parameter is presented and discussed. The ratio is calculated only for bending local ductility (µ) and global ductility, according to the two definitions under consideration (µGS y µGt), for dynamics and pushover analysis. The results are summarized in Table 10.
It is observed that the RLG values associated to µGt and dynamic analysis give unreasonable larger values particularly for the 10-story building; the reason for this is that, as mentioned in Section 7, the normalizing quantity Dy,t in Equation (5) may be very small. These values are larger than that of pushover analysis which in turn are larger than those of µGS and dynamic analysis. Considering that: (a) µGt of dynamic analysis fails addressing the normalizing quantity (Dy,t in Equation (5)) giving unreasonable large values of µGt, (b) pushover does not take into account dissipation of energy, (c) SCWB concept was followed in the building design, and (d) the local axial ductility is not important in framed steel building structures, the ratio of global to local ductility capacity, proposed in this study, is calculated as the ratio of µGS to µ of beams for dynamic analysis. Results from the table indicate that, for this definition, the mean values of RLG and the uncertainty in the estimation are larger for the 3- than for the 10-level building and larger for the EW than for the NS direction. Based on the values obtained in this study (0.41, 0.34, 0.35 and 0.32) a value of 1/3 is proposed for the RLG ratio. Thus, if local ductility capacity is stated as the basis for the design, say 15 or 12, the global ductility capacity can be estimated as 5 or 4.

10. Conclusions

A numerical investigation regarding ductility evaluation of steel buildings with moment resisting steel frames is conducted in this paper. Some steel model buildings and some strong motions used in the SAC steel project are used in the study. Bending local (µ) and tension local (µ) ductilities, story ductility (µS) as well as global ductility are studied; the most appropriate definition of global ductility is identified. Global ductility is calculated as the mean values of story ductilities (µGS) and as the ratio of the maximum inelastic top displacement to the top displacement when yielding occurs for the first time (µGt). The ductility values calculated according to dynamic analysis are compared to those of pushover analyses. Results of the study indicate that the ductility values significantly may vary with the strong motion, the ductility definition, the structural element, the story number, the type of analysis, and the structural model. The main conclusions are:
(1)
For beams and columns and the maximum seismic intensity, the mean values of µ tend to decrease with the story number; however, as expected, the beam values are much larger than those of columns (10.22 and 14.69 against 4.8 and 6.51). The µ values are only significant for exterior columns and much smaller than that of µ for beams. Even though tensional ductility demands are considered as an important issue in some research reports, and steel structural members may have a high ductility capacity in tension, these type of ductility is less relevant for the structural system under consideration.
(2)
The µS values for the maximum seismic intensity range from 3.84 to 4.28 and from 2.49 to 5.19 for the low- and medium-rise models, respectively. The values of µGS for dynamic analysis associated to the maximum seismic intensity (ductility capacity) are larger for the low- than for the medium-rise building (3.95 and 3.89 against 3.82 and 3.31). Unreasonable large values are obtained for µGt, values of 36.34 and 38.78 were obtained for the 10-level buildings. The implication of this is that the ductility definition based on the building top displacement, even though used in the profession, is not appropriate for the case of dynamic analysis, particularly for high buildings where significant structural contributions of the higher modes are expected.
(3)
As for the case of dynamic results, the ductility values obtained from pushover analysis are much larger for beams than for columns. For beams, the values of µ, µS and µGS obtained from pushover can be much larger than the corresponding values obtained from dynamic analysis.
(4)
The ratio (RLG) of global to local ductility (µ) is also calculated. Considering that: (a) µGt for dynamic analysis results in unreasonable large values, (b) pushover analysis does not take into account for dissipation of energy, (c) SCWB concept was followed in the building design, and (d) the local axial ductility is not important in framed steel building structures, the RLG ratio of dynamic analysis associated to bending local ductility is assumed to be the most appropriated definition. A value of 1/3 is proposed for this ratio. Thus, if local ductility capacity is stated as the basis for the design, say 15 or 12, the global ductility capacity can be estimated as 5 or 4.

Author Contributions

M.D.L.-T., A.R.-S. and E.B. performed the dynamic and push-over analyses for the 3-level model and built the plots and tables of the results obtained. J.B. and A.L.B. performed the dynamic analyses for the 9-level model and built the plots and tables of the results obtained. J.L.R.-S. and J. R.G.-C. performed the push-over analyses for the 9-level model and built the plots and tables of the results obtained. M.D.L.-T. and A.R.-S. interpreted, reviewed results and wrote the paper.

Acknowledgments

The research presented in this paper was financially supported by La Universidad Autónoma de Sinaloa (UAS) under grant PROFAPI-2015/235. 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 no conflicts of interest.

References

  1. International Code Council. International Building Code; International Code Council: Falls Church, VA, USA, 2009. [Google Scholar]
  2. National Research Council of Canada. National Building Code of Canada (NBCC); National Council of Canada: Ottawa, ON, Canada, 2005.
  3. Mexico City Building Code (MCBC). Complementary Technical Norms for Earthquake Resistant Design; MCBC: Mexico City, Mexico, 2004. [Google Scholar]
  4. Committee Euripean of Normalisation. Design of Structures for Earthquake Resistance; CEN: Brussels, Belgium, 2004. [Google Scholar]
  5. SAC. Steel Moment Frame Connections; Advisory No. 3 D-146; Structural Engineers Associated of California, Applied Technology Council and California University for Research in Earthquake Engineering: Sacramento, CA, USA, 1995. [Google Scholar]
  6. Reyes-Salazar, A. Ductility and Ductility Reduction Factor for MDOF systems. Struct. Eng. Mech. Int. J. 2002, 13, 369–385. [Google Scholar] [CrossRef]
  7. Uang, C.M. Establishing R (or Rw) and Cd Factors for Building Seismic Provisions. J. Struct. Eng. ASCE 1991, 117, 19–28. [Google Scholar] [CrossRef]
  8. Uang, C.M. Structural Overstrength and Limit State Philosophy in Seismic Design Provisions; Department of Civil Engineering, Northeastern University: Boston, MA, USA, 1991. [Google Scholar]
  9. Newmark, N.M.; Hall, W.J. Earthquake Spectra and Design Monograph Series; Earthquake Engineering Research Institute: Berkeley, CA, USA, 1982. [Google Scholar]
  10. Osteraas, J.D.; Krawinkler, H. Strength and Ductility Considerations in Seismic Design; Blume Earthquake Engineering Center, Stanford University: Stanford, CA, USA, 1990. [Google Scholar]
  11. ATC (Applied Technology Council). Tentative Provisions for the Development of Seismic Regulation Buildings; Report No. ATC-3-06; ATC: Redwood City, CA, USA, 1978. [Google Scholar]
  12. Hadjian, A.H. An Evaluation of the Ductility Reduction Factor Q in the 1976 Regulations for the Federal District of Mexico. Earthq. Eng. Struct. Dyn. 1989, 18, 217–231. [Google Scholar] [CrossRef]
  13. Miranda, E.; Bertero, V. Evaluation of Strength Reduction Factors for Earthquake-Resistant Design. Earthq. Spectra 1994, 10, 357–379. [Google Scholar] [CrossRef]
  14. Tiwari, A.K.; Gupta, V.K. Scaling of ductility and damage-based strength reduction factors for horizontal motions. Earthq. Eng. Struct. Dyn. 2000, 29, 969–987. [Google Scholar] [CrossRef]
  15. Arroyo-Espinoza, D.; Teran-Gilmore, A. Strength reduction factors for ductile structures with passive energy dissipating devices. J. Earthq. Eng. 2003, 7, 297–325. [Google Scholar] [CrossRef]
  16. Levy, R.; Rutenberg, A.; Qadi, K.H. Equivalent linearization applied to earthquake excitations and the R-µ-T0 relationships. Eng. Struct. 2006, 28, 216–228. [Google Scholar] [CrossRef]
  17. Karmakar, D.; Gupta, V.K. A parametric study of strength reduction factors for elasto-plastic oscillators. Sādhanā 2006, 31, 343–357. [Google Scholar] [CrossRef] [Green Version]
  18. Rupakhety, R.; Sigbjörnsson, R. Ground-motion prediction equations (GMPEs) for inelastic response and structural behavior factors. Bull. Earthq. Eng. 2009, 7, 637–659. [Google Scholar] [CrossRef]
  19. Fanaie, N.; Dizaj, E.A. Response modification factor of the frames braced with reduced yielding segment BRB. Struct. Eng. Mech. 2014, 50, 1–17. [Google Scholar] [CrossRef]
  20. Zahrai, S.M.; Jalali, M. Experimental and analytical investigations on seismic behavior of ductile steel knee braced frames. Steel Compos. Struct. 2014, 16, 1–21. [Google Scholar] [CrossRef]
  21. Zhai, C.H.; Wen, W.P.; Li, S.; Xie, L.L. The ductility-based strength reduction factor for the mainshock–aftershock sequence-type ground motions. Bull. Earthq. Eng. 2015, 13, 2893–2914. [Google Scholar] [CrossRef]
  22. Ghods, S.; Kheyroddin, A.; Nazeryan, M.; Mirtaheri, S.M.; Gholhaki, M. Nonlinear behavior of connections in RCS frames with bracing and steel plate shear wall. Steel Compos. Struct. 2016, 22, 915–935. [Google Scholar] [CrossRef]
  23. Nassar, A.; Krawinkler, H. Seismic Demands of SDOF and MDOF Systems; John, A., Ed.; Blume Earthquake Engineering Center, Standford University: Stanford, CA, USA, 1991. [Google Scholar]
  24. Santa-Ana, P.; Miranda, E. Strength Reduction Factors of Multi-Degree-of-Freedom Systems. In Proceedings of the 12th World Conference on Earthquake Engineering, Auckland, New Zealand, 30 January–4 February 2000. [Google Scholar]
  25. Moghaddam, H.; Mohammadi, R.K. Ductility reduction factor of MDOF shear-building structures. J. Earthq. Eng. 2001, 5, 425–440. [Google Scholar] [CrossRef]
  26. Elnashai, A.S.; Mwafy, A.M. Overstrength and force reduction factors of multistorey reinforced-concrete buildings. Struct. Des. Tall Build. 2002, 11, 329–351. [Google Scholar] [CrossRef]
  27. Medina, R.; Krawinkler, H. Evaluation of drift demands for the seismic performance assessment of frames. J. Struct. Eng. ASCE 2005, 131, 1003–1013. [Google Scholar] [CrossRef]
  28. Medina, R.; Krawinkler, H. Strength demand issues relevant for the seismic design of moment-resisting frames. Earthq. Spectra 2005, 21, 415–439. [Google Scholar] [CrossRef]
  29. De Stefano, M.; Marino, E.M.; Rossi, P.P. Effect of Overstrength on the Seismic Behaviour of Multi-Storey Regularly Asymmetric Buildings. Bull. Earthq. Eng. 2006, 4, 23–42. [Google Scholar] [CrossRef]
  30. Cai, J.; Zhou, J.; Fang, X. Seismic Ductility Reduction Factors for Multi-Degree-of-Freedom Systems. Adv. Struct. Eng. 2006, 9, 591–601. [Google Scholar] [CrossRef]
  31. Chopra, A.K. Dynamics of Structures; Prentice Hall: Upper Saddle River, NJ, USA, 2007. [Google Scholar]
  32. Mollaioli, F.; Bruno, S. Influence of side effects on inelastic displacement ratios for SDOF and MDOF systems. Comput. Math. Appl. 2008, 55, 184–207. [Google Scholar] [CrossRef]
  33. Vielma, J.C.; Barbat, A.H.; Oller, S. Seismic safety of low ductility structures used in Spain. Bull. Earthq. Eng. 2010, 8, 135–155. [Google Scholar] [CrossRef]
  34. Ceylan, M.; Arslan, M.H.; Ceylan, R.; Kaltakci, M.Y.; Ozbay, Y. A new application area of ANN and ANFIS: Determination of earthquake load reduction factor of prefabricated industrial buildings. Civ. Eng. Environ. Syst. 2010, 27, 53–69. [Google Scholar] [CrossRef]
  35. Ganjavi, B.; Hao, H. Effect of structural characteristics distribution of strength demand and ductility reduction factor of MDOF systems considering soil-structure interaction. Earthq. Eng. Eng. Vib. 2012, 11, 205–220. [Google Scholar] [CrossRef]
  36. Kumar, M.; Stafford, P.J.; Elghazouli, A.Y. Seismic shear demands in multi-storey steel frames designed according to Eurocode 8. Eng. Struct. 2013, 52, 69–87. [Google Scholar] [CrossRef]
  37. Kumar, M.; Stafford, P.J.; Elghazouli, A.Y. Influence of ground motions characteristics on drift demands in steel moment frames designed according to Eurocode 8. Eng. Struct. 2013, 52, 502–517. [Google Scholar] [CrossRef]
  38. Serror, M.H.; Diab, R.A.; Mourad, S.A. Seismic force reduction factor for steel moment resisting frames with supplemental viscous dampers. Earthq. Struct. 2014, 7, 1171–1186. [Google Scholar] [CrossRef]
  39. Reyes-Salazar, A.; Bojórquez, E.; Velazquez-Dimas, J.I.; López-Barraza, A.; Rivera-Salas, J.L. Ductility and ductility reduction factors for steel buildings considering different structural representations. Bull. Earthq. Eng. 2015, 13, 1749–1771. [Google Scholar] [CrossRef]
  40. Valenzuela-Beltrán, F.; Ruiz, S.E.; Reyes-Salazar, A.; Gaxiola-Camacho, J.R. On the Seismic Design of Structures with Tilting Located within a Seismic Region. Appl. Sci. 2017, 7, 1146. [Google Scholar] [CrossRef]
  41. Fanaie, N.; Shamlou, S.O. Response modification factor of mixed structures. Steel Compos. Struct. 2015, 19, 1449–1466. [Google Scholar] [CrossRef]
  42. Vuran, E.; Aydınoğlu, M.N. Capacity and ductility demand estimation procedures for preliminary design of coupled core wall systems of tall building. Bull. Earthq. Eng. 2016, 14, 721–745. [Google Scholar] [CrossRef]
  43. Gómez-Martínez, F.; Alonso-Durán, A.; De Luca, F.; Verderame, G.M. Ductility of wide-beam RC frames as lateral resisting system. Bull. Earthq. Eng. 2016, 14, 1545–1569. [Google Scholar] [CrossRef] [Green Version]
  44. Wang, B.; Dai, H.; Wu, T.; Bai, G.; Bai, Y. Experimental Investigation on Seismic Behavior of Steel Truss-RC Column Hybrid Structure with Steel Diagonal Braces. Appl. Sci. 2018, 8, 131. [Google Scholar] [CrossRef]
  45. Liu, G.; Lian, J.; Liang, C.; Zhao, M. Structural response analysis in time and frequency domain considering both ductility and strain rate effects under uniform and multiple-support earthquake excitations. Earthq. Struct. 2016, 10, 989–1012. [Google Scholar] [CrossRef]
  46. Hashemi, S.S.; Sadeghi, K.; Vaghefi, M.; Shayan, K. Evaluation of ductility and response modification factor in moment-resisting steel frames with CFT columns. Earthq. Struct. 2017, 12, 643–652. [Google Scholar] [CrossRef]
  47. Kang, J.D.; Mori, Y. Evaluation of a Simplified Method to Estimate the Peak Inter-Story Drift Ratio of Steel Frames with Hysteretic Dampers. Appl. Sci. 2017, 7, 449. [Google Scholar] [CrossRef]
  48. Reyes-Salazar, A.; Haldar, A. Dissipation of Energy in Steel Frames With PR Connections. Struct. Eng. Mech. Int. J. 2000, 9, 241–256. [Google Scholar] [CrossRef]
  49. Reyes-Salazar, A.; Haldar, A. Energy Dissipation at PR Frames Under Seismic Loading. J. Struct. Eng. ASCE 2001, 127, 588–593. [Google Scholar] [CrossRef]
  50. Reyes-Salazar, A.; Haldar, A. Seismic Response and Energy Dissipation in Partially Restrained and Fully Restrained Steel Frames: An Analytical Study. Steel Compos. Struct. Int. J. 2001, 1, 459–480. [Google Scholar] [CrossRef]
  51. Roeder, C.W.; Scheiner, S.P.; Carpenter, J.E. Seismic Behavior of Moment-Resisting Steel Frames: Analytical Study. J. Struct. Eng. ASCE 1993, 119, 1866–1884. [Google Scholar] [CrossRef]
  52. Leon, R.T.; Shin, K.J. Performance of Semi-Rigid Frames. In Proceedings of the Structure XIII Congress, Boston, MA, USA, 2–5 April 1995; pp. 1020–1035. [Google Scholar]
  53. Nader, M.N.; Astaneh, A. Dynamic Behavior of Flexible, Semi-rigid and Rigid Frames. J. Constr. Steel Res. 1991, 18, 179–192. [Google Scholar] [CrossRef]
  54. Osman, A.; Ghobarah, A.; Korol, R.M. Implications of Design Philosophies of Seismic Response of Steel Moments Frame. Earthq. Eng. Struct. Dyn. 1995, 24, 127–143. [Google Scholar] [CrossRef]
  55. FEMA. State of the Art Report on Systems Performance of Steel Moment Frames Subject to Earthquake Ground Shaking, Prepared for the SAC joint Venture Partnership by Helmut Krawinkler; Report No. FEMA 355-C; FEMA: Washington, DC, USA, 2000.
  56. Uniform Building Code (UBC). International Conference of Building Officials (ICBO), Whittier, CA, USA. 1997. [Google Scholar]
  57. BOCA. 12th Edition Building Officials & Code Administration International; BOCA: Country Club Hills, IL, USA, 1993.
  58. Carr, A. “RUAUMOKO” Inelastic Dynamic Analysis Program; University of Cantenbury, Department of Civil Engineering: Christchurch, New Zealand, 2011. [Google Scholar]
  59. Chen, W.F.; Atsuta, T. Interaction Equations for Biaxially Loaded Sections; Fritz Laboratory Report (72-9); Leigh University: Bethlehem, PA, USA, 1972. [Google Scholar]
  60. Aboosaber, M.; Hines, E.M. Modeling Reserve System Performance for Low-Ductility Braced Frames; Technical Report TUSSR-2011/1; Tufts University Structural Systems Research (TUSSR), Department of Civil and Environmental Engineering: Medford, MA, USA.
  61. Amadio, C.; Chiara Bedon, CH.; Fasan, M.; Pecce, M.R. Refined numerical modelling for the structural assessment of steel-concrete composite beam-to-column joints under seismic loads. Eng. Struct. 2017, 138, 394–409. [Google Scholar] [CrossRef]
  62. Christidis, A.A.; Dimitroudi, E.G.; Hatzigeorgiou, G.D.; Beskos, D.E. Maximum seismic displacements evaluation of Steel frames from their post-earthquake residual deformation. Bull. Earthq. Eng. 2013, 11, 2233–2248. [Google Scholar] [CrossRef]
  63. Amadio, C.; Chiara Bedon, C.H.; Fasan, M. Numerical assessment of slab-interaction effects on the behaviour of steel-concrete composite joints. J. Construct. Steel Res. 2017, 139, 397–410. [Google Scholar] [CrossRef]
  64. Liew, A.; Gardner, L.; Block, P. Moment-Curvature-Thrust Relationships for Beam-Columns. Structures 2017, 11, 146–154. [Google Scholar] [CrossRef] [Green Version]
  65. Galambos, T.V. Structural Members and Frames; Prentice Hall: Englewood Cliffs, NJ, USA, 2016. [Google Scholar]
Figure 1. Isometric view of the 3-level building.
Figure 1. Isometric view of the 3-level building.
Applsci 09 00200 g001
Figure 2. Isometric view of the 9-level building.
Figure 2. Isometric view of the 9-level building.
Applsci 09 00200 g002
Figure 3. Bending local ductility (µ), EW direction, 3-level building; beams: (a) Sa = 0.4 g, (b) Sa = 0.8 g, (c) Sa = 1.2 g; columns: (d) Sa = 0.4 g, (e) Sa = 0.8 g, (f) Sa = 1.2 g.
Figure 3. Bending local ductility (µ), EW direction, 3-level building; beams: (a) Sa = 0.4 g, (b) Sa = 0.8 g, (c) Sa = 1.2 g; columns: (d) Sa = 0.4 g, (e) Sa = 0.8 g, (f) Sa = 1.2 g.
Applsci 09 00200 g003
Figure 4. Bending local ductility (µ), EW direction, 10-level building; beams: (a) Sa = 0.2 g, (b) Sa = 0.6 g, (c) Sa = 0.8 g; columns: (d) Sa = 0.2 g, (e) Sa = 0.6 g, (f) Sa = 0.8 g.
Figure 4. Bending local ductility (µ), EW direction, 10-level building; beams: (a) Sa = 0.2 g, (b) Sa = 0.6 g, (c) Sa = 0.8 g; columns: (d) Sa = 0.2 g, (e) Sa = 0.6 g, (f) Sa = 0.8 g.
Applsci 09 00200 g004
Figure 5. Story ductility (µS) values; NS direction; 3-level building: (a) Sa = 0.4 g, (b) Sa = 0.8 g, (c) Sa = 1.2 g; 10-level buildings: (d) Sa = 0.2 g, (e) Sa = 0.6 g, (f) Sa = 0.8 g.
Figure 5. Story ductility (µS) values; NS direction; 3-level building: (a) Sa = 0.4 g, (b) Sa = 0.8 g, (c) Sa = 1.2 g; 10-level buildings: (d) Sa = 0.2 g, (e) Sa = 0.6 g, (f) Sa = 0.8 g.
Applsci 09 00200 g005
Table 1. Beam and columns sections for Models 1 and 2.
Table 1. Beam and columns sections for Models 1 and 2.
ModelStoryColumnsGirder
ExteriorInterior
11/2W14X257W14X311W33X118
2/3W14X257W14X312W30X116
3/ROOFW14X257W14X313W24X68
2−1/1W14X370W14X500W36X160
1/2W14X370W14X500W36X160
2/3W14X370W14X500, W14X455W36X160
3/4W14X370W14X455W36X135
4/5W14X370, W14X283W14X455, W14X370W36X135
5/6W14X283W14X370W36X135
6/7W14X283, W14X257W14X370, W14X283W36X135
7/8W14X257W14X283W30X99
8/9W14X257, W14X233W14X283, W14X257W27X84
9/ROOFW14X233W14X257W24X68
Table 2. Earthquake records, NS component. Tn: Elastic period; ED: epicentral distance; M: Magnitude moment; PGA: Peak ground acceleration.
Table 2. Earthquake records, NS component. Tn: Elastic period; ED: epicentral distance; M: Magnitude moment; PGA: Peak ground acceleration.
No.PlaceDate (Day/Month/Year)StationTn (s)ED (km)MPGA (cm/s2)
1Landers, California28/06/1992Fun Valley, Reservoir 3610.11317.3213
2MammothLakes, California27/05/1980Convict Creek0.16126.3316
3Victoria09/06/1980Cerro Prieto0.16376.1613
4Parkfield, California28/09/2004Parkfield; JoaquinCanyon0.17156609
5PugetSound, Washington29/04/1965Olympia Hwy Test Lab0.17896.5216
6Long Beach, California10/03/1933UtilitiesBldg, Long Beach0.2296.3219
7Sierra El Mayor, Mexico04/04/2010El centro, California0.21777.2544
8Petrolia/Cape Mendocino, California25/04/1992Centerville Beach, Naval Facility0.21227.2471
9Morgan Hill24/04/1984GilroyArraySta #40.22386.2395
10Western Washington13/04/1949Olympia Hwy Test Lab0.22397.1295
11San Fernando09/02/1971Castaic-Old Ridge Route0.23246.6328
12MammothLakes, California25/05/1980Long Valley Dam0.24136.5418
13El Centro18/05/1940El Centro-ImpVallIrrDist0.27127350
14Loma Prieta, California18/10/1989Palo Alto0.29476.9378
15Santa Barbara, California13/08/1978UCSB Goleta FF0.36145.1361
16Coalinga, California02/05/1983ParkfieldFaultZone 140.39386.2269
17Imperial Valley, California15/10/1979Chihuahua0.4196.5262
18Northridge, California17/01/1994Canoga Park, Santa Susana0.6166.7602
19Offshore Northern, California10/01/2010Ferndale, California0.61436.5431
20Joshua Tree, California23/04/1992Indio, Jackson Road0.62266.1400
Table 3. Statistics for µ, 3-level building.
Table 3. Statistics for µ, 3-level building.
Type of MemberStoryDynamic AnalysisPUSH
NS DirectionEW Direction
Sa/g ValuesSa/g Values
0.40.60.81.01.20.40.60.81.01.2
MVCVMVCVMVCVMVCVMVCVMVCVMVCVMVCVMVCVMVCV
BEAMS12.610.304.270.355.990.408.230.3610.220.372.190.213.560.204.560.236.330.177.810.2013.01
22.480.284.600.337.170.269.570.2511.860.262.260.224.150.255.870.267.690.269.320.2716.64
32.010.334.310.376.710.308.940.2611.200.281.990.373.730.365.590.367.670.329.300.3016.32
COLUMNS11.060.091.320.281.950.442.650.473.700.511.010.021.250.191.720.242.220.322.870.391.00
21.000.051.000.071.090.101.250.291.410.291.000.071.010.021.030.051.210.231.290.191.40
31.000.041.000.061.060.051.120.151.230.231.000.031.000.041.020.021.130.141.260.248.71
Table 4. Statistics for µ, 10-level building.
Table 4. Statistics for µ, 10-level building.
Type of MemberStoryDynamic AnalysisPUSH
NS DirectionEW Direction
Sa/g ValuesSa/g Values
0.20.40.60.80.20.40.60.8
MVCVMVCVMVCVMVCVMVCVMVCVMVCVMVCV
BEAMS21.930.345.600.409.910.4114.690.351.600.334.650.398.100.3911.890.371.60
31.650.445.230.338.940.3212.820.291.290.314.380.407.510.3910.640.384.16
41.830.385.470.288.610.2812.610.231.410.264.580.417.580.4010.800.368.02
51.890.485.180.278.180.2011.360.201.410.324.210.366.970.379.960.3311.06
61.560.484.110.326.460.238.680.201.200.243.440.335.490.357.710.3615.62
71.420.543.750.405.910.257.770.231.260.303.240.355.140.347.210.2920.75
82.610.414.960.296.580.218.360.222.530.304.700.256.330.227.970.2323.15
93.010.555.040.346.480.267.840.212.660.354.820.226.220.197.880.2022.38
102.110.623.740.455.090.366.320.241.490.533.270.354.850.266.340.2521.15
COLUMNS21.000.071.710.263.610.666.510.631.000.041.410.492.870.474.800.481.00
31.000.051.000.041.570.722.540.741.000.041.000.081.900.431.540.461.00
41.000.061.000.041.000.051.870.671.000.091.000.071.000.051.370.411.00
51.000.031.000.031.000.041.640.491.000.031.000.041.000.061.310.291.00
61.000.011.000.021.000.031.460.471.000.051.000.041.000.041.260.321.00
71.000.011.000.021.000.021.000.041.000.021.000.051.000.031.000.041.00
81.000.031.000.011.000.011.000.031.000.031.000.031.000.051.000.061.60
91.000.021.000.011.000.011.000.021.000.011.000.031.000.021.000.043.60
101.000.021.000.001.000.011.000.011.000.001.000.021.000.011.000.018.73
Table 5. Statistics for (µ) for columns of the 3-level building.
Table 5. Statistics for (µ) for columns of the 3-level building.
StoryNS DirectionEW Direction
Sa/g ValuesSa/g Values
0.40.60.81.01.20.40.60.81.01.2
MVCVMVCVMVCVMVCVMVCVMVCVMVCVMVCVMVCVMVCV
11.000.001.000.001.000.001.000.001.000.001.000.001.000.001.000.001.000.001.000.00
21.000.001.000.001.000.001.030.131.260.621.000.001.000.001.070.291.110.441.230.57
31.390.462.540.852.340.792.200.612.580.591.150.221.960.763.940.734.480.733.570.54
Table 6. Statistics for µS, 3-level building.
Table 6. Statistics for µS, 3-level building.
StoryDynamic analysisPUSH
NS directionEW direction
Sa/g valuesSa/g values
0.40.60.81.01.20.40.60.81.01.2
MVCVMVCVMVCVMVCVMVCVMVCVMVCVMVCVMVCVMVCV
11.610.241.990.342.810.403.440.513.950.551.480.201.900.312.800.372.890.344.280.698.63
21.700.182.200.292.860.463.310.423.990.501.650.172.150.242.580.313.090.323.540.457.66
31.530.161.930.202.690.272.990.233.900.331.520.242.100.242.720.293.190.373.840.426.87
Table 7. Statistics for µS, 10-level building.
Table 7. Statistics for µS, 10-level building.
StoryDynamic AnalysisPUSHEr
NS DirectionEW Direction
Sa/g ValuesSa/g Values
0.20.40.60.80.20.40.60.8
MVCVMVCVMVCVMVCVMVCVMVCVMVCVMVCV
21.210.232.230.513.770.565.190.591.040.161.730.322.660.423.950.515.90
31.180.252.110.373.570.514.400.501.050.201.900.382.570.443.630.456.59
41.290.242.350.453.520.464.790.481.240.272.110.362.910.413.880.437.62
51.260.272.130.343.100.384.120.451.130.241.960.332.740.363.410.396.88
61.150.291.880.382.560.363.380.431.040.161.700.312.540.352.930.355.07
71.120.311.770.402.320.333.140.381.080.191.660.312.060.342.770.403.97
81.640.542.230.472.880.273.550.351.440.242.210.242.810.283.140.353.78
91.760.482.320.302.540.283.310.451.570.262.320.312.710.223.400.342.43
101.330.361.740.342.180.362.490.351.170.261.600.282.170.292.630.335.90
Table 8. Global ductility (µGS y µGt), 3-level building.
Table 8. Global ductility (µGS y µGt), 3-level building.
DEFINDynamic AnalysisPUSH
NS DirectionEW Direction
Sa/g ValuesSa/g Values
0.40.60.81.01.20.40.60.81.01.2
µGS1.612.042.793.253.951.552.052.703.053.897.72
µGt2.453.214.315.846.912.573.304.884.886.508.37
Table 9. Global ductility (µGS y µGt), 10-level building.
Table 9. Global ductility (µGS y µGt), 10-level building.
DEFINDynamic AnalysisPUSH
NS DirectionEW Direction
Sa/g ValuesSa/g Values
0.20.40.60.80.20.40.60.8
µGS1.272.082.943.821.181.912.573.315.08
µGt8.4011.7022.4036.344.8913.5225.5138.787.98
Table 10. Statistics of the RLG ratio.
Table 10. Statistics of the RLG ratio.
DefinitionModelDynamic AnalysisPushover
EONS
MVCVMVCV
µGS3-LEVEL0.410.450.340.270.50
10-LEVEL0.350.330.320.190.36
µGt3-LEVEL0.79NA0.59NA0.55
10-LEVEL2.87NA2.59NA0.56
NA = not applicable.

Share and Cite

MDPI and ACS Style

Llanes-Tizoc, M.D.; Reyes-Salazar, A.; Bojorquez, E.; Bojorquez, J.; Lopez-Barraza, A.; Rivera-Salas, J.L.; Gaxiola-Camacho, J.R. Local, Story, and Global Ductility Evaluation for Complex 2D Steel Buildings: Pushover and Dynamic Analysis. Appl. Sci. 2019, 9, 200. https://0-doi-org.brum.beds.ac.uk/10.3390/app9010200

AMA Style

Llanes-Tizoc MD, Reyes-Salazar A, Bojorquez E, Bojorquez J, Lopez-Barraza A, Rivera-Salas JL, Gaxiola-Camacho JR. Local, Story, and Global Ductility Evaluation for Complex 2D Steel Buildings: Pushover and Dynamic Analysis. Applied Sciences. 2019; 9(1):200. https://0-doi-org.brum.beds.ac.uk/10.3390/app9010200

Chicago/Turabian Style

Llanes-Tizoc, Mario D., Alfredo Reyes-Salazar, Eden Bojorquez, Juan Bojorquez, Arturo Lopez-Barraza, J. Luz Rivera-Salas, and Jose R. Gaxiola-Camacho. 2019. "Local, Story, and Global Ductility Evaluation for Complex 2D Steel Buildings: Pushover and Dynamic Analysis" Applied Sciences 9, no. 1: 200. https://0-doi-org.brum.beds.ac.uk/10.3390/app9010200

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop