Analysis of the global second-order effects on irregular reinforced concrete structures using the natural period of vibration 2

Resumo The χT parameter, a simplified method recently presented, allows to estimate the global second-order effects on reinforced concrete frames using the natural period of vibration. This parameter was developed based on the fact that both natural period of vibration and global second-order effects depend essentially on the stiffness and mass matrices of the structure, being thus related. In this paper, numerical analyses are conducted on nine models with different patterns of irregularity in terms of geometry in plan and stiffness. The main purpose of these analyses is to evaluate the applicability of the χT parameter in asymmetric structures as well as that can present torsional modes as the fundamental mode of vibration. In addition, different hypotheses are tested in order to verify the influence of the different modes of vibration in the structural sensitivity to global second-order effects. Results of the simplified analyses were compared to the final bending moment values obtained through a nonlinear numerical analysis considering the P-Δ effect. It is observed that the parameter χT is a promising indicator for a simplified estimation of the global secondorder effects for concrete frames, especially when higher modes of vibration are taken account in the analysis.


Introduction
In the latest years, computational advancement and improvement of the material technology lead to the design of slender buildings more susceptible to global second order effects.Enhancement in computer's processors also enabled more complex structural analysis that were not assessed in the past because of limited computational resources.Along with this, using high-strength concretes allows the construction of higher buildings, which are consequently, more flexible [1].Taller structures are subjected to higher horizontal forces, due to wind load.The combination of these forces with vertical loadings results in higher displacements of structural nodes and, consequently, significant increase of bending moment in its elements.Therefore, the structure must be assessed in its deformed shape.The additional forces from this analysis are known as the second-order moments [2].In order to ensure structural safety, this geometrical nonlinearity must be considered in the structural stability analysis.In addition, reinforced concrete frames have a nonlinear material behavior.Concrete cracking, concrete creep and, the rebar yielding are known as material physical nonlinearities [3].Excessive lateral displacement must be verified considering both geometrical and material nonlinearities.An option for the second-order analysis is to adapt the stiffness matrix (K) to include second-order effects [4].In this method, the stiffness matrix (K) is composed by the elastic stiffness matrix (K E ) and geometric stiffness matrix (K G ); this problem are usually solved through iterative incremental methods [5].In structures with multiple degrees of freedom, however, incremental method has high computational cost.Therefore, direct and simplified methods are preferred to the second-order analysis.These methods adopt approximations and simplifications to describe the structural behavior; their efficiency relies on the selection of the most important variables and its influence on the structural behavior.Although these procedures are simplified, they can be quickly applied and provide satisfactory results [6].Rutenberg [7] proposed a direct approach to obtain the global second-order effects.The stiffness matrix is modified through the insertion of fictitious diagonal braces and columns with negative cross section area and flexural stiffness, respectively.The fictitious elements simulate the geometric stiffness matrix, which includes second-order effects while equilibrium is verified.Thus, a linear analysis can be carried out and no iterations are needed.Wilson and Habilullah [8], White et al. [9], [10] and Vandepitte [11] are other examples of direct approaches.Simplified methods are also included in many countries design standards.The American standard AISC/LRFD [12] introduces B1 and B2 coefficients to evaluate second-order effects in steel structures.The second-order forces are obtained through linear analysis and, are vastly used by steel structural design engineers.To define the requirement of second-order analysis, ABNT NBR 6118:2014 [2] adopts the simplified parameters α and γ z .Parameter α classify the structure as nonsway or sway frames.Global second-order effects can be neglected for the latter.
The γ z coefficient [13], besides providing a classification regarding structural susceptibility to second order effects, it is also used to estimate second-order forces through first-order analysis.ABNT NBR 6118:2014 [2] allows material nonlinearity to be considered applying reduction factors to the element stiffness.It estimates the structural elements stiffness change and assumes different values for each type of element.The coefficient γ z is calculated assuming that lack of stability will only occur laterally without rotational component.Structures, which the center of mass does not correspond to the center of rigidity, however, may present as the first stability mode the torsional mode.In this case, stability analysis with this coefficient may not be sufficient to assess structural integrity [14].Therefore, to apply γ z coefficient, stability modes must also be assessed.Buildings global stability is mainly influenced by lateral stiffness and vertical loads applied [15].Based on this, Statler et al. [16] have demonstrated that it is possible to evaluate the structure susceptibility to global second order effects using the natural period of vibration in steel structures.Recently, Reis et al. [17] presented the development of the simplified χ T parameter.This parameter uses the natural period of vibration to estimate bending moments considering non-linear effects.To account for the global second-order effects, first-order moments are amplified in reinforced concrete frames.Analyses carried out by Reis et al. [17], however, applied the χ T parameter only in structures with floor plan symmetry and no stiffness variation along its height.These characteristics diverge completely from most constructions.In Brazil, floor plan designs are usually asymmetric, which leads to: (i) asymmetrical structures, that have eccentricity between the center of mass and center of rigidity and makes the structure vulnerable to significant torsional forces, (ii) vertical stiffness variation along the building's height caused by the sudden change in geometry between two floors or change in the cross-section area, (iii) bracing elements non-parallel or non-symmetric to principal axes, and, (iv) structures with multiple floor plan that do not have the same projection.Thus, it is essential, to validate the parameter χ T, that a greater number of analysis are carried out considering structural models with similar characteristics to existing buildings.Furthermore, results obtained by Reis et al. [17] presented conservative values for the critical columns on the assessed models.Therefore, an improvement in moment amplification methods due to global secondorder effects must be investigated.In the field of structural dynamics, modal superposition is commonly used.The total structure's response is obtained by superposing the response for each mode shape.Therefore, parameter χ T can be improved considering higher vibration modes.In this context, the current work aims to assess the application of the χ T parameter, developed by Reis et al. [17], in reinforced concrete structures with different irregularities on the floor plan and, along its height.The study presented herein is intended to cover the major structures designed in Brazil and, to understand the limitation of the methods.Two new hypotheses are developed and verified.The importance of the mode shape direction in the application of the simplified parameter is assessed, as well as, the influence of higher frequency modes.
Analysis of the global second-order effects on irregular reinforced concrete structures using the natural period of vibration

χ T Parameter
The χ T parameter was developed by Reis et al. [17] aiming to evaluate reinforced concrete structures susceptibility to global secondorder effects; in addition, it is possible to estimate the final internal forces through a simplified first-order analysis.The parameter was developed applying Rayleigh's method [18], in which coordinate generalization and the principle of virtual work were used.Rayleigh's method [18] is based on the principle of energy conservation.If the physical properties are known, the natural period of vibration of an element can be estimated by the selection of a displacement function (ψ) that respects the essential and/or natural boundary conditions.In Rayleigh method [18], a satisfactory approximation for the structure fundamental period of vibration is obtained adopting any displacement function that satisfies the essential (kinematic) boundary conditions.Adopting, however, a shape function that does not satisfies the natural boundary conditions, impose additional restrains and external constrained to the system [19].These restrictions increase the system stiffness and, consequently, influence significantly on the natural period of vibration estimated by the method.The χ T parameter theoretical basis relies on the application of generalized coordinates to represent the dynamic displacement of a cantilever beam-column, with mass and stiffness equivalent to a regular structure.The column is subjected to a uniform distributed horizontal loading F(x) along its height.The floor masses (m p ) are concentrated along its length with constant spacing of 3,0 m.The column mass is uniformly distributed along the column height.The model representation is shown in Figure [1].The χ T parameter complete original formulation obtained by Reis et al. [17] is presented in Equation [1]: Onde, M 1 -First-order moment; M 2 -Final second-order moment; H -Building's height; T -Natural period of vibration; n -Number of floors; g -Acceleration of gravity; κ pav -Ratio between the sum of all floor weight and the structure total weight.A simplification was developed based on the influence of each A simplification was developed based on the influence of each polynomial term of the original equation and adopting κ pav equal to 0,80.The complete and simplified equation presented a maximum divergence of 1,5%.The equation details, as well as simplification, can be found in Reis et al. [17].The χ T parameter final version is presented in Equation [2]: Equation [2] depends essentially on the natural period of vibration, building's height and number of floors.The generalized system adopted by Reis et al. [17] for the development of the χ T parameter suppose that the fundamental mode of vibration is flexural.Thus, the natural period of vibration adopted to global second-order analysis must be related to the first flexural mode of vibration in each assessed direction.Then, the natural period of vibration adopted in the χ T parameter equation for moment amplification (M 2 /M 1 ) must be related to this vibration mode.Reis et al. [17] assessed the χ T parameter for five different structure models with one axis of symmetry at least.The moment amplification results were compared with the results of a P-D second-order analysis.It was shown that the simplified χ T parameter represents well the second-order effects in reinforced concrete structures.The authors, however, commented that more studies must be carried out to adjust and validate the simplified methodology, mainly on the application of irregular structures which can be vulnerable to torsional effects.In addition, the χ T parameter values presented by the models in Reis et al. [17] were conservative for the critical columns.The adopted displacement function does not satisfy the natural boundary conditions, which leads to a conservative moment amplification for the critical columns.Mamone et al. [20] extended χ T parameter using a displacement function that satisfies simultaneously essential and natural conditions, which lead to a more flexible system.Considering this adaptation, the χ T parameter application was assessed in two simplified models.The first model is a cantilever beam-column and the second model is a simple frame with regular geometry.The results obtained for the moment amplification were closer to the values obtained for the critical columns; the models studied, however, present floor plan symmetry and the stiffness does not vary along its height.The χ T parameter equation can also lead to greater values when considering a modified period of vibration with the superposition of kinematic mode shapes of the structure (vibration modes).Thus, the moment amplification due to global second-order effect can be estimated considering natural period of vibration related to superior vibration modes.The current study intends to contribute to the verification of the method's limit of application and, to its development.The parameter is applied: (i) in new structural models with different patterns of irregularities on the floor plan and (ii) stiffness variation along the building's height, and (iii) in structures vulnerable to torsional forces, which the first mode of vibration is torsional.Two new hypotheses are developed and tested to consider higher vibration modes in the χ T parameter equation.The hypotheses present criteria that, based on the effective modal mass, select the natural period of vibration to be applied in the χ T parameter equation.These hypotheses are explained later in the following section.

Model description
Spatial frames with different geometries and irregularities on its floor plan or along the building's height were modelled herein and, are presented in Figures [2] to [6].The models have multimodal characteristics, which assist the assessment of the influence of higher vibration modes on the proposed simplified parameter.Models 7 and 8 are real structures designed in Brazil, aimed to allow a realistic study of the parameter χ T. Observing the geometry of both structures (Figure [4] and

Figure 4 Model 7 geometry
Analysis of the global second-order effects on irregular reinforced concrete structures using the natural period of vibration [5]), one can note that there are many stiffness variations between floors.These geometry variations characterize a stiffness variation along the building's height.Thus, models 7 and 8 will contribute to evaluate the parameter χ T application in structures with this type of irregularity.All models were generated with the software S4 [21] used in the structural design of reinforced concrete structures.All geometries were assessed varying the number of floors until the γ z coefficient reached the limit value of 1.3.This is the maximum value allowed in the application of the simplified method suggested by ABNT NBR 6118:2014 [2].Model 4 was an exception to this γ z limitation, a limit of 1.2 was imposed since the building height respective to this value was already considered too high.The models' characteristics are summarized in Table [1].The number of floors, maximum height, material mechanical properties, vertical loading per floor and, γ z value for both orthogonal directions in the plan (x and y) are presented.

Design Considerations
After defining the maximum number of floors, the software S4 [21] model was exported to SAP2000 [22] software.SAP2000 [22] is a finite element software extensively used by structural designers worldwide because of its user-friendly interface to carry out linear, nonlinear and modal analysis.In SAP2000 [22], columns and beams are modelled as unidimensional elements (beam element).These elements have six degree of freedom per node; three translations and, three rotations.Due to slabs' high horizontal stiffness, they are modelled as rigid diaphragms.Thus, the slabs behave similarly to a rigid body that does not deform axially [23].

Figure 6 Model 9 geometry
Analysis of the global second-order effects on irregular reinforced concrete structures using the natural period of vibration Load cases and loading combinations were also exported from S4 software [21].Loading combinations include the effects of simultaneous application of vertical (dead and live loads) and horizontal (wind load) loadings.For each wind loading direction, a different loading combination was created.The wind load was considered in the structure's principal directions 0º(Fx+), 90º(Fy-), 180º(Fx-) and 270º(Fy+), according to standard ABNT NBR 6123:1988 [24].Figure [7] illustrates wind load directions considered in the models.
Wind load magnitude were divided by pressure ranges according to each structure height, aimed to more realistic represent the wind pressure acting on the structure.To consider the rigid diaphragm effect, however, these loadings were modeled as concentrated horizontal forces in each floor level.All floor nodes displacements equivalent to the diaphragm center of mass displacement.Thus, concentrated horizontal loadings were applied in these nodes on each floor level.The material nonlinearity was considered approximately using the reduction factors prescribed by ABNT NBR 6118:2014 [2].The columns were considered encastred on the foundations.For all loading combinations, two analyses types were carried out in SAP2000 [22]: considering only the overturning moment due to horizontal wind load, in order to obtain the first-order bending moments (M 1 ); and, a second analysis considering the P-Δ effect to determine the second-order bending moments (M 2 ).The obtained ratio between second-order and first order bending moments (M 2 /M 1 ) from both analyses express the moment amplification due to global second-order effects.P-Δ effect are taken in account in SAP2000 [22] considering a modified geometrical stiffness matrix K G in addition to the linear-elastic stiffness matrix K E .Simultaneously, modal analyses were carried out to determine the natural periods of vibration and its respective mode shapes, as well as the effective modal mass participation factor.Modal analysis, in SAP2000 [22], consists of a generalized eigenvalue problem, which the eigenvectors are orthogonal to mass and stiffness matrices.Free undamped systems are considered.The eigenvalues represent the natural period of vibration and their respective eigenvector are the vibration modes.For each geometry studied, after defining the maximum height and the maximum number of floors, all analyses were repeated remov-   In models 1 to 8, the same structural configuration was adopted, independently of the number of floors.Model 9 has the same geometry of Model 2. For each analysis, the building structural elements were redesigned once a floor level was removed.Model 9 crosssection variation is depicted in Figure [6], in which the ground floor geometry for a 21-storey and 4-storey structures are presented.Therefore, Model 9, due to element redesign, also contributed to evaluate the χ T parameter in realistic buildings.More details about the evaluated hypotheses are presented as follows.

Hypotheses
The original χ T parameter hypothesis, presented by Reis et al. [17], shall consider the first flexural vibration mode according to the assessed direction.In order to verify the influence of different vibration modes to the global second-order effects, two hypotheses were proposed.Each one has specific criteria to select the natural period of vibration that will be used in the simplified χ T parameter.
Results are compared to the final moment amplification factors of the structures columns.Table [2] summarizes the hypotheses for the natural period of vibration selection, including the original hypothesis (Hypothesis I), adopted by Reis et al. [17].

Hypothesis I (original from REIS et al. [17]) -1 st flexural mode is considered direction to assess second-order effects
Hypothesis I is the original hypothesis presented by Reis et al. [17].
The selected natural period of vibration is computed to determine the χ T parameter; only first flexural vibration mode is taken in account to consider global second-order effects.Reis et al. [17] have used the effective modal mass coefficient as the value used to determine the first flexural mode of vibration.According to Reis et al. [17] the natural period of vibration must be chosen if the respective mode shape presents an effective modal mass coefficient greater than 35% in the considered direction of the global second-order analysis.
Cases, in which none of the first vibration modes present greater effective modal mass than 35% in any direction, are hardly found.If this is the situation, it must be considered the vibration mode which depict the greatest effective modal coefficient in the assessed direction when compared to the orthogonal directions.This hypoth-esis contributes to further development of Reis et al. [17] studies by verification of its applicability in structural models with several irregularity patterns.

Hypothesis II -First mode shape of vibration (fundamental vibration mode)
In this hypothesis, the selected vibration mode is the fundamental mode shape, independently from the assessed direction and the vibration mode shape.This analysis intends to verify the consequences of applying the simplified parameter without the concern regarding the vibration mode shape (flexural, torsional or flexuraltorsional).In addition, this hypothesis will help understanding the application of the fundamental vibration mode as an indicator of the global second-order effects in three-dimensional frames.

Hypothesis III -Average of the natural period of vibration weighted by the modal participation factor
The selected period of vibration in this hypothesis is the composition of natural period of vibration weighted by the effective modal mass in the direction considered the greatest second-order effects impact.Herein, the weighted natural period of vibration can be obtained by Equation [3]: In which, T n is the natural period of vibration of the n-th mode and, U n the modal participation factor of the n-th mode in the considered direction.This hypothesis aims to verify the influence of the higher modes of vibration in the structure sensitiveness to the global second-order effects.
To illustrate hypothesis III natural period selection method, the first vibration modes of a 21-storey building (Model 2 geometry) are shown in Table [3].Weighted by the effective modal mass, the weighted natural period of vibration in the assessed direction is calculated.In this example, it is intended to calculate the moment amplification (M 2 /M 1 ) in the y-direction using the χ T parameter.
From the presented effective modal mass in Table [3], the weighted period of vibration of hypothesis III is 5.49 s.The definition of the number of vibration modes to be used in this hypothesis was performed through a comparative study, in which the number of selected vibration modes were varied, according to the total effective modal mass participation.The total effective modal mass varied from 75% to 90% in the assessed direction.Analysis of the global second-order effects on irregular reinforced concrete structures using the natural period of vibration The goal is to define the minimum number of vibration modes to be considered on Hypothesis III.Analyses were carried out for the models from Figures [2] to [6] with 10-storey, 20-storey and the maximum number of for each model.The results of this comparison are presented in Table [4].Difference between the results obtained considering 75% and 90% of the total effective modal mass is small.The maximum difference, considering all the structures, was 2.4%.Therefore, the contribution can be disregarded and, it is recommended to adopt the number of natural periods that adds to 75% of total effective modal mass in the assessed direction (x or y).

Results and discussion
Modal analysis and second-order analysis results of each model are depicted in Figures [8]- [12].From the modal analyses, curves were plotted referent to the effective modal mass percentage of the first vibration mode as a function of the number of floors for each model.From the global second-order analysis, the bending mo-ment amplification factor (M 2 /M 1 ) of all columns in each direction were plotted for every node in the model (gray dots), as a function of the number of storey.Moment amplification factors lower than an unity were disregarded.Columns located at the base were defined as critical columns and, its moment amplification factors are highlighted in the graph, as an analysis criterion.
Adopting the criteria of each hypothesis presented above, the natural period of vibration was selected and used to calculate the simplified χ T parameter.The dashed vertical line in each graph indicates a change in the first vibration mode type (flexural, torsional, flexural-torsional, etc.) of the structure as a function of the number of storey.Table [5] presents a descriptive legend aimed to summarize the analyses discussed herein; each hypothesis is represented by a different type and line color.The simplified γ z parameter, presented in the Brazilian standard ABNT NBR 6118:2014 [2], is also plotted for comparison and validation of the obtained results.Following, a detailed discussion is presented for the modal analysis and the global second-order effects analysis through the simplified χ T parameter.

Table 5
Informative legend for the result graphs

Table 4
Comparative analysis between χ T values obtained with hypothesis III, varying the number of modes of vibration in function of total mobilized mass -75% and 90%

Modal analysis
The structure's dynamic characteristics are investigated through the modal analysis.The first vibration mode shapes and the respective natural period of vibration were determined.Figures [8] and [9] present the effective modal mass percentage for the fundamental vibration mode in the x-direction, y-direction and rotation along the z axis as a function of each assessed geometry.The graphs depict the fundamental vibration mode according to the number of storey.Therefore, the direction with more significant effective modal mass values indicate the structure's type of the vibration mode shape.
In Figure [8], for instance, the y-direction curve for Model 1 pres-

Figure 8
Effective modal mass of the first mode of vibration for models 1, 2, 3, 4, 5 and 6 Analysis of the global second-order effects on irregular reinforced concrete structures using the natural period of vibration ents an increase on the effective modal mass for the models with 14 storey or more in comparison to the other directions.Below 14 storey, the predominant mode shape is torsional mode shape.
Table [6] summarizes the fundamental vibration mode type of all models as a function of the number of storey.
Table [6] depicts also a great number of geometries with torsional vibration mode, which, indirectly, represents structures vulnerable to torsional effects.It must be pointed out that this behavior is usually avoided in the structural design because of the requirement of special rebar detailing related to torsional forces and additional lateral displacements that can be caused due to the combined action of torsion and bending moments.These structures, however, are studied to verify the simplified χ T parameter capacity in evaluating the global second-order effects in structures with significant eccentricity between the rigidity center and mass center.

Global second-order effects analysis ( χ T parameter)
The global second-order effect was evaluated through a P-Δ analysis.Its bending moment amplification results for each column are plotted for every node.It is intended to compare the bending moment amplification from a more accurate analysis, the P-Δ analysis, and the moment amplification from the χ T parameter.

Figure 9
Effective modal mass of the first mode of vibration for models 7, 8 and 9  Analysis of the global second-order effects on irregular reinforced concrete structures using the natural period of vibration Figure [10], [11] and [12] show the results of the global secondorder effects in the x and y-direction, where all hypotheses results are compared, for each model and direction.The moment amplification obtained with the P-Δ analyses are presented as a function of the number of storey.
It must be emphasized that, the comparison intends to verify the simplified χ T parameter accuracy to evaluate global second-order effects in regular and irregular three-dimensional frames.Then, it is expected to evaluate limitations in the use of the χ T parameter.
For instance, using the fundamental natural period of vibration, independently of the direction, the χ T parameter is evaluated in structures with torsional and flexural-torsional vibration modes.Also, the influence of higher vibration modes is studied.The bending moment amplification factor calculated using the simplified χ T parameter (Hypothesis I) and the bending moments in the critical columns increase with the number of storey in each model, as mentioned by Reis et al. [17].This systemic increase reinforces that the natural period of vibration and moment amplification due to global second-order effects are related.
As expected, the most flexible direction presented the greatest values for the moment amplification factor.However, it is noticed that when there is a change in the first vibration mode shape, the results for the χ T parameter and the critical columns are significantly different.One can clearly note that, observing the results for Model 5 in the y-direction.At the point that the first mode of vibration changes from torsional to flexural in y-direction, the parameter curve has a pronounced shift.In this model, for structures with 23 storey or less, the first flexural mode of vibration in the y-direction is the second global vibration mode.Thus, in hypothesis I, it is the second natural period of vibration that must be used (Figure [11]).For structures, with more than 23 storey, the fundamental vibration mode (first vibration mode) is flexural.For the 24-storey building, specifically, the fundamental mode of vibration in y-direction has an effective modal mass less than 35%.Therefore, in Hypothesis I, the second natural period of vibration was adopted for the simplified parameter calculation.Only in structures with more than 25 storey, the first vibration mode in y-direction present effective modal mass participation greater than 35%, which explains the shift in Hypothesis I curves on Figure [11].
The fundamental mode of vibration meets the hypothesis criteria and, therefore, the natural period of vibration that must be used is the fundamental natural period of vibration, which increases significantly the χ T parameter.It can be observed that Hypothesis I, adopted originally by Reis et al. [17], reaches conservatives values for the moment amplification factor, when compared only with the moment amplification factor for the critical columns.It represents a upper bound for the other structure columns [17].Comparing the χ T parameter (hypothesis I) with γ z parameter, Model 1 presented values similar, with maximum difference of 0,71% between both parameters.This behavior was also noted in Model 6 for the x-direction.The maximum difference was 2,60% in the 30-storey structure.If this comparison is carried out for the other models, however, the differences were pronounced, indicating conservative values for the χ T parameter in relation to γ z parameter.One can note that the torsional mode is the fundamental vibration mode as shown in frames with 14 storey or less for Model 1, frames with 10 storey or less for Model 2 and, frames with 23 storey or less for Model 5.
For these structures, it is observed that the χ T parameter given the first vibration mode (fundamental period -hypothesis II) is similar to the χ T parameter from the original hypothesis proposed by Reis et al. [17] (hypothesis I), when considering the direction with smaller stiffness.Although the fundamental vibration mode is torsional, the results similarity is due to the second mode shape being flexural; the maximum difference being of 8,0%.Table [7] presents the first three vibration modes characteristics for Models 1, 2, 5 and 9, that showed similar results.For Model 1, the maximum difference between hypotheses I and II curves was 0,66% in a 14-storey building.It can be observed in Table [7], that the first vibration mode of Model 1 was torsional, and the second vibration mode was flexural in the y-direction.Thus, the χ T parameter calculated with hypothesis II (first mode) considers the first natural period of vibration (3,39 s), while hypothesis I (first flexural mode of vibration in the considered direction) considers the value corresponding to the second period of vibration (3,38 s).The natural periods of vibration are very similar, which explains the small difference between the curves.One can also note this observation in Model 2 in the y-direction.The maximum difference  In Table [7] shall be noted that the first natural period of vibration (2,95 s) and the first flexural natural period of vibration in the ydirection (2,84 s) have close values.The difference between the natural period of vibration values is 3,87%, which leads the simplified χ T parameter to similar values in both hypotheses.
When considering these structures flexibility in the x-direction (natural period of vibration presented in Table [7], there is considerably difference between the value of the moment amplification factor for the fundamental period of vibration (hypothesis II) and the other hypotheses.If the fundamental period of vibration is used in threedimensional frames, independently from the assessed direction, the moment amplification for the global second-order effects are conservative when applied to the direction with greatest stiffness.Besides that, once the simplified method is used, the structural engineer is led to verify the building vulnerability to torsion and the authors understand that is an important consideration to be taken in account during the structural design and structural system conception.On the other hand, the effective modal mass participation factor, aforementioned, is capable to accurately assess the moment amplification factor given that the structure is susceptible to vibration modes with torsional components (flexural-torsional mode shapes).However, considering the value of the natural period of vibration corresponding to the flexural vibration mode, the simplified parameter approaches the moment amplification for the critical columns as an upper bound, particularly in the direction with lower stiffness.
The objective of the hypothesis III development is to verify the influence of higher vibration mode shapes in the χ T parameter behavior.To assess this hypothesis results, the period of vibration adopted in the parameter is calculated by weighted average of the natural period of vibration, considering the effective modal mass participation in the studied direction.
Hypothesis III results presented values similar to the moment amplification for the critical columns in each direction.Figure [13] presents a comparison between the maximum differences obtained using the χ T parameter -hypothesis III and the moment amplification factor for the critical columns of each model.The maximum difference between hypothesis III curves and the moment amplification for the critical columns was 4.65%.This difference was observed considering all geometries from the nine models studied herein, varying the number of storey and the assessed direction for the lateral load (x and y).
The authors would like to emphasize the importance in considering higher modes of vibration in the simplified χ T parameter analyses, with the equation presented by Reis et al. [17].For Model 9, which were obtained with the redesign of Model 2 as storey were introduced, hypothesis III results are similar to the moment amplification in the critical columns.Note that this structure presents stiffness variation along the building's height and floor-plan asymmetry; these considerations are more realistic situation in structural design.One can note that, for this model, there are cases in which the moment amplification for the critical columns reduces, as the number of floors increases.This reduction can indicate a case in which the critical column has its participation in the structural global stiffness reduced because of the new design.Or, a case in which, after the redesign, another column was considered critical; in other words, presented higher internal forces.Figure [14] shows the changes of the elements' dimensions in the ground floor for 12-storey and 13-storey structures.Both models are based on Model 9 geometry.Columns P9 and P17 are highlighted in Figure [14]; these are the columns with the greatest bending moment in the 12-storey and 13-storey buildings, respectively.The overturning moment values of columns P9 and P17 for the 12-storey and 13-storey buildings of Model 9 is presented in Table [8].Columns' section and moment of inertia is also presented in Table [8].P9 is the critical column in the 12-storey building.The overturning moment is 144,9 kN.m and the moment amplification factor (M 2 /M 1 ) is equal to 1,16.For the 13-storey building, however, column P17 presents an overturning base moment of 205,5 kN.m, which is greater than the moment in the  Analysis of the global second-order effects on irregular reinforced concrete structures using the natural period of vibration P9 column (156,5 kN.m) and the moment amplification factor is 1,14.One can note that the moment amplification factor was reduced, even with an increase in the number of storey, from 12 to 13 storey.The reduction is explained by the fact that the critical column is different in both analyses.Hypothesis III curves also presented shifts in the same points in which these changes were noticed.Therefore, the simplified method was also sensitive to these structural behavior changes.
In general, applying the original equation (Hypothesis I) resulted in conservative moment amplification factor in relation to the critical columns.This can be explained by adoption of a displacement function that only satisfies the essential boundary conditions.The weighted natural period of vibration improved significantly the accuracy of the moment amplification factor.The weighting is a correction for the system's stiffening arising from the original formulation.Once the structural system becomes more flexible, there is improvement of the simplified analysis representation for the critical columns.

Conclusions
This paper presented the applications and limitations of a simplified parameter, based on the natural period of vibration, used to evaluate the global second-order effects on irregular reinforced concrete structures.The objective was to evaluate the χ T parameter applicability to irregular structures susceptible to torsion.Nine models with different geometries were assessed.The models presented different patterns of irregularity on the floor plan and stiffness variations.In addition, two new hypotheses for the χ T parameter were proposed.It was intended to verify the importance of the vibration mode direction on the χ T parameter and, also verify the contribution of the higher vibration modes in the global second-order effects.
The obtained results have demonstrated that the χ T parameter is a promising estimator for the global second-order effects, even in irregular reinforced concrete structures.More accurate results are found when considering higher vibration modes.Most of the models studied herein consisted of structures vulnerable to torsion.In other words, structures in which the first natural vibration mode is torsional or, structures with high effective modal mass for the rotation along the vertical axis (z-axis).From the results, it is clear that the simplified method proposed herein is able to capture the moment amplification factor for mixed vibration mode (flexural torsional) or predominantly torsional.Moreover, the analyses indicate the importance of observing the structure fundamental vibration modes during structural design.Structures with torsional vibration mode as the fundamental mode indicates great eccentricities between the center of mass and center of rigidity, which leads to the conclusion that the structure is vulnerable to torsion and, the effects shall be carefully considered in the structural analysis.
In order to avoid this behavior on the design of reinforced concrete frames, some corrective measures can be taken.
In addition, the simplified methodologies usually used to evaluate the global second-order effects have some restrictions on structures with torsional behavior; these methodologies do not present any indication of the fundamental vibration mode or information about the direction of the critical mode.Therefore, the simplified analysis with the χ T parameter become an advantageous tool for the structural engineer.The coupled analy-sis considering the modal parameters (period of vibration, vibration mode and, effective modal mass) lead to a better understanding of the type and direction of the structure vibration mode shape.Furthermore, this observation will provide to the structural engineer basis for any necessary precautions related to the structure vulnerability to torsional effects.For the application of the simplified parameter, it must be assured that the parameter is calculated with the periods of vibration associated with the flexural vibration mode for each direction.
The global second order analysis applying the χ T parameter with hypothesis II (fundamental vibration mode) presented excessively conservative results in the direction of the higher structural stiffness.Therefore, it is important to use the χ T parameter considering each orthogonal direction in the floor plan.The effective modal mass must be greater than 35% for the first flexural vibration mode in the assessed direction.
Hypothesis III, which considers an average natural period of vibration weighted by the effective modal mass, presented moment amplification similar to the moment amplification obtained for the critical columns in the design.The latter moment amplification was obtained through a P-D analysis.The maximum difference for the amplification curves of this hypothesis was 4,65% when compared with the moment amplification of the critical columns.
It must be emphasized that Model 7 and 8 are built structures in Brazil.And, Model 9 structures were redesigned as new storey were introduced, thus, representing real structural systems.These structures have stiffness variation along the building's height.The χ T parameter with hypothesis III was able to satisfactorily represent the moment amplification for the critical.Therefore, even if the stiffness changes along the building's height, as the modal analysis considers the global stiffness and mass matrices, the simplified parameter is capable of determining the moment amplification for the critical columns (critical in the structural design).Finally, it must be emphasized that the χ T parameter equation used herein was developed considering a uniformly distributed lateral loading, which represents the wind load.New studies, however, may be carried out to adequate the parameter for wind loading variation along the building's height.In order to evaluate the global second-order effects in structures vulnerable to torsional effects, methodologies that consider the dynamic loading effects must be used.For instance, the synthetic wind methodology can be used in further studies.

Figure 7
Figure 7Principal directions considered for the lateral wind loading I -1 st Flexural Mode Reis et al. (2017)) Hypothesis I (1 st flexural mode) χ T(II -1 st Mode of Vibration) Hypothesis II (fundamental mode of vibration) χ T(III -Weighted Average) Hypothesis III (weighted average) Columns Columns Amplification M 2 /M 1 γ Z Coefficient γ Z values Critical columns Critical columns amplification M 2 /M 1

F 6 F
Figure 10 Moment amplification M 2 /M 1 in function of the natural period of vibration in the assessed direction -Models 1, 2, 3 Torsional (rotation along z-axis) 2º 3.44 Flexural in x-direction 3º 3.34 Flexural in y-direction F. F. LEITÃO | G. H. SIQUEIRA | L. C. M. VIEIRA JR. | S. J. C. ALMEIDA between using hypothesis II (first mode) and the hypothesis I (first flexural mode -Reis et al. [17]) was 1,29% in a 10-storey building.

Figure 13
Figure 13Maximum differences between hypothesis III values and moment amplification for the critical column in each direction of the assessed models

Figure 14 F
Figure 14Cross section of P9 and P17 columns in the ground floor for 12-storey and 13-storey structures -Model 9

Table 1
Assessed models characteristics

Table 2
Summary of each hypothesis for the natural period of vibration selection [8]flexural mode in the considered direction -original hypothesis (REIS et al.[8])II1 st mode of vibration (fundamental mode) III Average of the natural period of vibration weighted by effective modal mass

Table 3
Example of the first 12 modes of vibration of a 21-storey structure generated with Model 2

Table 6
Fundamental mode of vibration direction as function of the number of storey for the assessed models

Table 7
Structures modes of vibration of Models 1, 2, 5 and 9

Table 8
Overturning base moment of critical columns for 12-storey and 13-storey structures -Model 9