Acessibilidade / Reportar erro

Global stability analysis of structures and actions to control their effects

ABSTRACT

In this moment in which civil engineering is undergoing a phase where structural projects have been developed with structural systems composed of different and complex elements, some methods and criteria are used for the purpose of evaluating important aspects with regard to global and local stability. Among them, it is necessary to mention the parameters of instability a and ?z. In this sense, this work has the objective to present the basic concepts of the instability parameters a and ?z in accordance with what is clearly defined in the Brazilian standard ABNT NBR 6118; to present the results of simulations of models in the Brazilian structural software TQS varying the stress of compression in the columns in order to relate these values with the stability parameters.

Keywords:
global stability; gama-z; structural analysis

RESUMO

Nesse momento em que a Engenharia civil vem passando por uma fase onde os projetos estruturais têm sido elaborados com sistemas estruturais compostos por elementos diferenciados e complexos, alguns critérios e métodos são utilizados com a finalidade de avaliar aspectos importantes no que diz respeito à estabilidade global e local. Entre eles, faz-se necessário citar os parâmetros de instabilidade α e Gama-z. Nesse sentido, este artigo tem o objetivo de apresentar os conceitos básicos dos parâmetros de instabilidade α e Gama-z de acordo com o que é claramente definido na norma brasileira ABNT NBR 6118; apresentar os resultados de simulações de modelos no software estrutural Brasileiro TQS variando a tensão de compressão nos pilares com a finalidade de relacionar estes valores com os parâmetros de estabilidade.

Palavras-chave:
estabilidade global; gama-z; análise estrutural

1. Introduction

The analysis of global stability of structures has become very important nowadays. Civil Engineering, as well as other branches of Engineering, has experienced major advances in their fields. This can be noticed by the total amount built or about to be built by the different companies who operate in construction and design. Not only by the significant number of services, this increasing request of Civil Engineering also becomes cause and incentive for the development of new technical production methods aimed at better results.

Major projects are launched to supply the demand, which comes day after day, whether in residential, commercial or industrial construction. Behind this demand, technical production brings an increase for goods and the requirement for the excellence.

Best practices in design and development projects have been considered in the viability procedures of the projects.

One of the advances is evidenced by the improvement design procedures. Factors previously considered only in special buildings (tall buildings or engineering structures), such as the second order effects and global and local stability of the structures, now become necessary in the evaluation in more joint ventures.

These advances are used by most of the current versions of software to simulate with greater precision the effects to which the structure is subjected in different situations.

The Brazilian regulatory standard "NBR 6118:2014 - Concrete Structures Project - Procedure" establishes, in item 15.5, guidelines for dismissal of the global second-order efforts. To analyze stability, the standard uses two parameters to estimate the second-order efforts: (alfa) and (Gama-Z).

The parameter , the simpler of the two, is used to evaluate the global stability in the structure, but is not able to estimate the second-order efforts. This parameter allows classifying a structure as fixed or mobile nodes. Parameter , on the other hand, by means of its formulation, allows estimation of the second-order efforts.

In this context, this paper is presented, aiming to discuss the main concepts involved in the study of global stability of structures, making use of the stability parameters and .

2. Brief history of stability parameters: the beginning

The study on the parameters of stability was started in 1967 by Hurbert Beck and Gert König (1967, cited by Vasconcelos, 1991) after presentation of a doctoral dissertation Germany. This study analyzed the structure of a building with many floors, which were considered the columns braced by rigid walls with a slim structure. Studies led to a complicated solution of the differential equation with variable coefficients, but after a mathematical simplification, the equation was reduced to Bessel functions. Based on the solution of this equation, a parameter called "stability coefficient a" was proposed.

This coefficient was able to provide information of what is the displacement and definitions about fixed and mobile nodes of the structure under analysis. By simplifying displacement calculations of the structure, it states that fixed node structures would be those in which the addition of internal forces (existing moments) was less than 10% of the internal forces already present. For these cases, the value of the coefficient a would be below 0.6.

2.1 Differential equation

Starting from the differential equation of a straight bar of length L, with constant section and made of a linear elastic material, subject to an axial load P at the ends, one has:

where y(x) represents the displacement curve of the bar in the transverse direction and the flexural rigidity. From the equation of the straight bar, it follows that:

where is treated such as a stability coefficient and its value is related to the conditions of support from the bar ends. The value of is low when the bar is embedded and high when the bar is biarticulated or mono crimped.

2.2 The differential equation according Beck e König

The equation proposed by Beck and König (1967, cited by Vasconcelos, 1991) at the time of their studies, was:

when:

In this equation, is the horizontal distributed load across the height L of the building.

After some mathematical operations, the equation of the value of found by Beck and König was:

where p is the distributed load to which the bracing structure is subjected and is the distributed load at the braced structure.

In table 1 the limit values found by Beck and König according with the number of floors are shown.

These results were obtained considering that the bracing structure of the building could be replaced by a structure composed of only

Table 1
Value of α acr found by Beck and König

a single column whose geometric characteristics were equivalent to the original. Thereby, it was decided that the sum of the stiffness of the bracing structure was equal to the rigidity of a single column. To calculate the equivalent rigidity in cases of bracing columns, the horizontal force acting on the spatial frame should be applied, in order to determine the value of the deformation with the horizontal loads. This deformation must be the same as the one obtained from the equivalent column subject to the same horizontal loads.

Once calculated the value, it must be verified to make sure it does not exceed the limit value. In affirmative case, modify the dimensions of the structural elements in order to get below the limit values.

Studies of parameters for stability analysis were continued by other authors with interest in the subject (CEB e ABNT). In 1978, the CEB started using the values described by Beck & König, as well as ABNT NBR 9062 (Brazilian code pre-shaped structure). In 1985, Mario Franco (1985) studied again how to obtain the limit values of the stability coefficient for buildings with more than three floors finding values ranging close to 0.6 in structures with wall-like columns, frames and combination frames and wall-like columns.

Vasconcelos (1985) found for structures with 1, 2 and 3 pavements values as described in Table 2.

In 1986, Augusto C. Vasconcelos took over studies by Beck & König, but this time with the advent of the computer, he could assess the same structure previously studied discretized, instead of continuous structure. The results were not different for up to three floors, and from four floors the difference between and was decreasing with the increase of the amount of floors.

Table 2
Values of α acr in according with the number of floors

This study resulted in the suggestion of a formula for calculating according to the number n of floors. Follows the formula:

In 1990, the CEB decides to drop the considerations on the calculation of stability parameter. In this context, Mario Franco and Augusto C. Vasconcelos (19976 VASCONCELOS, A. C. e FRANÇA, R. C. (1997) - Um método Simplificado e Muito Preciso para Avaliação dos Momentos de 2a Ordem em Edifícios Altos Usuais - Coletânea de Trabalhos sobre Estabilidade Global e Local das Estruturas de Edifícios. São Paulo, 1997., Franco, M. and Vasconcelos, AC, cited Vasconcelos, 1991) introduced the first concept of the amplification coefficient of the moment in the state of Rio de Janeiro, as an auxiliary tool in the study of global stability. Through the P-Delta process, come up to the idea of .

When the structure is requested by the loads, it undergoes a displacement altering its original configuration. From this new configuration, new internal forces are created, which in turn generate new displacements. In this continuous process, it is noted that at every step new strains are generated which, compared to the previous results, can be represented by a geometric progression. This progression has the following ratio:

where:

and is the eccentricity of first order, and it is the first order moment.

Years after the definition of the moment amplification , several studies were conducted on this topic. In the year 1993, the TQS creates the Spatial Frame module in which it begins to use the parameters of stability and . Correlations between and were also formulated as approaches for practical applications. Table 3 has some of the correlations that were presented in articles and studies.

3. Study of global stability of structures

The study of global stability of the structures has gained a high

Table 3
Correlations of α acr and ϒz

Figure 1
Relation moment-curvature (extracted by NBR 6118)

importance position because it has become common to development projects in which buildings have been increasingly high, slender and with features aimed at excellent performance, low consumption of materials, structural systems with speed of execution without giving up the structural safety of the guarantee to the ultimate limit state.

Accordingly, the verification of the global stability can be checked through the global stability parameters, which are described in the standard NBR 6118, in the items 15.5.2 and 15.5.3.

Using these parameters allows us to estimate data on the stability of the structure and estimate the second-order effects in a simplified manner. For this reason, it is necessary to present some concepts related to second-order effects and that should be considered in the project.

3.1 Related concepts to global stability of structures

3.1.1 Physical nonlinearity

According to NBR 6118, in section 15.3, structural analysis considering the effects of second order shall take into account the effects of physical nonlinearity. The nonlinearity is related to the behavior of the material, which in the case of reinforced concrete can be considered a non-linear material. The effects caused by creep, shrinkage and concrete cracking affect the final value of concrete elastic modulus, hence the physical non-linearity of concrete. For consideration of these effects, the rule establishes the concept of moment-curvature, shown in Figure 1.

According Moncayo (20113 MONCAYO, W. J. Z (2011). Análise de Segunda Ordem Global em Edifícios com Estrutura de Concreto Armado. Master's degree dissertation, USP. São Carlos, 2011.), the advantage in using the moment-curvature relationship is due to the fact that is possible to directly obtain the value of rigidity EI, which is used in the structural analysis. In the case of the study of global stability, the NBR 6118 states, in section 15.7.3, that the consideration of the approximated nonlinearity takes as stiffness of structural elements the following values:

Slabs:

Beams:

para

Columns:

where:

is the moment of inertia of the gross concrete section, including, where applicable, the co-operating tables (T section),

A's is the compression reinforcement when double armor,

As is the tensile armor,

Eci is the initial modulus of elasticity of the concrete, given by ;

fck is the characteristic compressive strength of concrete in MPa.

The approximate stiffness values must be used according to the conditions described in paragraph 15.7.3 of NBR 6118.

3.1.2 Geometric nonlinearity

The consideration of the effects of geometric nonlinearity is related to changes that may occur in the geometry of the structural elements. Like the physical nonlinearity, it is also determined by analyzing the deformed structure. Specially for high buildings, the consideration of the effects of geometric nonlinearity becomes important because it generates increased internal forces due to vertical loads and horizontal displacements.

For the consideration of geometric nonlinearity it is common to use the P-? process, where is conducted the study considering the deformed structure.

3.1.3 Stiffness of structural elements

The rigidity of the structural elements has a direct influence on the overall stability of a structure. In a conventional structure, consisting of columns, beams and slabs, the main elements that will work to ensure stability are the columns and beams. This fact can be seen taking the NBR 6118 approximate stiffness values for the slabs, which are lower compared to the stiffness of the beams and columns. In fact, the slabs will work in another property giving overall stability study: The hard diaphragm. The slabs are considered as rigid elements in their own plane, aligning the displacements at all points of this plane.

Another important consideration to be made is the consideration of hard cores in the buildings. Normally they are formed by large columns with "C" shape, and located in the stairs and elevators. This type of structure has a high flexural rigidity in the analysis of stability and horizontal displacement parameters, contribute significantly.

3.2 Global stability parameters

Following studies on the effects of second-order structures, the NBR 6118 code, in chapter 15 - "Instability and second-order effects", provides guidelines on how to apply test methods and analysis of second-order effects.

For their use, it is considered that it is applied to structures with bars subject to flexion-compression efforts that must disregard the twist. The concepts apply to shells structures, walls and wall-beams. The second order effects are those in which the equilibrium analysis is done using the deformed configuration.

The code adopts that to disregard the second-round effects the increases caused by these should not exceed 10% in the reactions and the loads of the structure. It should also ensure that the calculated bearing capacity should not be surpassed.

Figure 2
Equivalence of stiffness - analogy with column in balance

3.2.1 Stability parameters (a)

The structure classification as the lateral displaceability of the nodes is performed by the Stability parameter a . Above a certain limit value, the structure is classified as of mobile nodes, and therefore should be considered second-order analysis. It is emphasized that the parameter is not able to estimate the effects of second order.

Assuming initially that the bracing structure may be represented by a single column clamped at the base and free at the top, an equivalent Ejk stiffness of the sum of the stiffness of the bracing columns of the structure (figure 2) and being constant over the height H of the building, we have the following value for the parameter of instability:

In equation 9, is the sum of all the vertical loads of the structure. It will be a caveat that the amounts to be used are characteristic.

The index "d" indicates design values.

The condition for the structure to be fixed nodes is that , thereby equal to:

The limit value , when used in structures with less than four floors, is generally well applicable to most building structures. However, when the bracing structure is formed exclusively by wall-like columns and frames, the value of should be changed. The new limit values are described below:

- For bracing with wall-like columns;

- For mixed bracing (wall-like columns + frames);

- For bracing in frames.

The equivalent stiffness of the bracing elements (framed systems and mixed systems) of the structure is calculated by using the lateral characteristic load over the entire height of the building and determining its roof top displacement .

Obtained the deformation on top of the building is done with that for a new column equivalent balance, with constant equivalent stiffness, under the action of the same loading, this new column in balance suffer the same displacement at its end cantilevered. Assuming that the load is constant over the entire height of the building, it becomes

where:

is the lateral load characteristic;

is the total height of the building;

is the displacement at the top of the building.

3.2.2 Coefficient ϒz

The Brazilian code gives the coefficient assessing the presence or absence of second-order global efforts in framed structures with at least four floors.

For the determination of its value, it is necessary to make a first order analysis for each loading case. The formula for is:

where:

is the moment of tipping. It is the sum of the moments of all horizontal forces in the considered combination, with design values, in relation to the base structure.

is the sum of the multiplication of the vertical forces acting on the structure, in the combination design values considered, for horizontal displacement of the respective points of application.

The above values are obtained by first order analysis.

For the structure to be considered as a fixed node structure, the condition is that:

The condition above being satisfied, the structural calculation with fixed nodes can be performed on each compressed element alone. The element will be with the end linked to the other structural elements. According to NBR 6118: 2014, in section 15.6 - "nodes fixed structure Analysis", under the action of horizontal forces the structure is always calculated as displaceable.

The fact that the structure can be treated as fixed allows to disregard only the verification of the second order global effects. Local second-order effects must be considered.

The code NBR 6118 says that the global analysis of second order provides only internal forces at the ends of the bars and should be performed an analysis of the local second order effects along the axes of the compressed bars.

On the other hand, when the value of obtained from a first order analysis is higher than , the structure should be considered of mobile nodes. In this case, the effects of physical and geometric nonlinearity should be necessarily be considered. It is therefore necessary to check the local and global second order

Figure 3
Floor plan without wall-like columns

effects.

When the value of is in the range of , NBR 6118: 2014 says that the global second order internal forces can be considered in an approximate way by adopting different values of stiffness for structural elements, which should be the same as shown in item 3.1.1. Finally, the final solution to the approximate consideration of the global second order internal forces is to increase the horizontal forces from the combination of load seen by , obtained after analysis of the first order structure. This process only applies to .

In other cases in which the structure will be of mobile nodes and must be verified the analysis of second-order effects considering refined methods such as .

4. Structural modeling

Figure 4
Structure

To apply the concepts of structural stability, this article will simulate of a building in a structural modeling software. The CAD/TQS software is one of the tools available in the Brazilian market to professionals involved in the preparation of structural design of reinforced concrete buildings and one of the most widely used.

Next will be presented the floor plan of the building used for simulation, variations of each model used, the description of vertical and horizontal loads to which the building is subject, the conditions and considerations adopted for structural analysis and presentation of values of and .

Table 4
Loads used on the structural model

Tabela 5
Types of models in accordance with the stress on the columns

Table 6
Preliminary design of Column P1 - Variation V1

4.1 Structural model adopted

The building consists of 10 floors (Figure 4) in accordance with the floor plan shown in Figures 3 and 5. Were adopted two structural solutions: one without wall-like columns (figure 3) and the other with wall-like columns (Figure 5). This condition was set to compare which are the effects when using wall-like columns.

The CAD / TQS version 17.8 software provides six models for design and structural analysis. In the case of this article "Model VI", also known as Flexibility Model with Slabs, was adopted. This model is the most appropriate when is used to structural analysis of reinforced concrete buildings, because the model considers a three-dimensional space frame consisting of beams, columns and slabs. The model consists of bars, which simulate the slabs, beams and columns.

The structure was modeled considering the following horizontal loads provided in the code: plumb and wind; they were considered automatically by the software during the structural analysis. The wind speed was considered to 30 m/s, factor S1, S2 and S3 are respectively 1, category IV-B and 1.

The vertical loads considered in the floors below according to Table 4. In order to simplify the study and facilitate the comparison of results between each model, we chose to use surface charges for the masonry load.

4.2 Pre-dimensioning of the columns

The principal objective of this article is to vary the compression stress on the columns in order to relate these values with the stability parameters. For this, the pre-sizing of the columns was determined by the method of the areas of influence on each floor of the building. In accordance with what was previously presented by Table 4, the total load surface is 10 kN/m2.

The area of influence of each pavement was determined by the structural calculation software, but could be obtained manually considering that the distribution of loads occurs simply, i.e., a span between two columns in a defined direction, the influence of loads in each of these columns, is given until the middle of this span. Reworking this step for the other directions will be obtained an area. This area is equivalent to the distributed load that this column absorbs.

The characteristic strength of compression adopted for design was:

The variation of compression stress in the columns is related according with the criteria that the reference value for design would be:

Figure 5
Floor plan with wall-like columns

The models studied followed stress variations according to Table 5.

As an example, Table 6 presents how is determined the preliminary design of column P1, without wall-like columns, following V1 variation:

where:

Pav is the floor considered;

A.I. is the Influence Area in m2;

C/Pav is the Full charge adopted on the pavement (given in Table 4), units in kN/m2;

C. Pav is the Load obtained in the pavement for the column P1 (), units in kN; C. Acum: Load accumulated on the floor, where: ;

b and h are the column section dimensions of the considered floor, in cm;

P.P. is the Weight of the column, obtained by multiplying the cross-sectional area of the column by floor height (2,88 m) and the dead weight of the concrete (25 kN/m2);

Column Stress for each floor is a value obtained by , units in kN/m2;

It is observed that the stress in the foundation and first floor do not exceed the limit value for the V1 range: 11.6 MPa = 1.16 kN/m2.

To simplify the presentation of the results, were developed the following models with the variations V1 to V5 (Table 5), where the dimensions of the columns were pre-determined by each floor plan with and without wall-like columns. Ten different models were

Table 7
Types of models

Table 8
Used combinations

determined in accordance with Table 7.

To obtain the results, the load combinations were considered in models presented in Table 8.

The dead and love loads were increased by and wind loads with coefficient . These values were obtained from the Brazilian code NBR 6120: 1980 - Loads for Building Structures Design.

5. Results

The proposed models were processed in the CAD/TQS software, through three-dimensional frame analysis, and the values obtained of e are presented on table 9.

6. Conclusion

After the presentation of the results using Table 9 and Graph 1, it can be seen that the values of obtained for the models calculated without wall-like columns are higher compared to models with wall-like columns, which was expected because the presence of the wall-like columns increases the bending stiffness of the building.

It is also observed that the variation of stress on the columns of the models with wall-like columns results in small differences of the coefficient , with values between 1.064 and 1.161. On models without wall-like columns, the variation of the column stress modify significantly the values of , making them between 1.197 and

Table 9
Results

1.426. It is possible to conclude, therefore, that the variation of the column stress is not the determining factor in acceptable values, but the presence of elements that contribute to increase stiffness bending. In this case, the author has chosen to use wall-like columns, but alternatives can also be used, like beams with high stiffness, walls of reinforced concrete, increase of the columns sections in the direction of less rigidity of the structure.

In this article, the models with wall-like columns have values below the limit, being possible to disregard the second order analysis of the structure. It is only needed to perform first-order analysis and increase the internal forces according to the code recommendations. In models without wall-like columns, all of them with stress

Graph 1
Comparison of obtained results

over 70% of the reference value, second-order analysis is needed.

In interpreting the values of it is observed that most of the models set up with mobile nodes.

Thereby, this article presented the basic concepts of stability parameters a and ; presented the results of models simulations on the structural software CAD/TQS varying the compressive stress in the columns for the purpose of relating these values with the stability parameters and and presented actions that could be taken in buildings to minimize the effects of global instability.

References

  • 1
    ASSOCIAÇÃO BRASILEIRA DE NORMAS TÉCNICAS (2014). Projeto de estruturas de concreto - Procedimento, NBR 6118, Rio de Janeiro, RJ.
  • 2
    ASSOCIAÇÃO BRASILEIRA DE NORMAS TÉCNICAS (1980). Cargas para cálculo de estruturas de edificações, NBR 6120, Rio de Janeiro, RJ.
  • 3
    MONCAYO, W. J. Z (2011). Análise de Segunda Ordem Global em Edifícios com Estrutura de Concreto Armado. Master's degree dissertation, USP. São Carlos, 2011.
  • 4
    VASCONCELOS, A. C. (1997). Origem dos Parâmetros de Estabilidade α e γz - Coletânea de Trabalhos sobre Estabilidade Global e Local das Estruturas de Edifícios. São Paulo, 1997.
  • 5
    VASCONCELOS, A. C. e FRANCO, M. (1997) - Avaliação Prática dos Efeitos de 2a Ordem em Edifícios Altos - Coletânea de Trabalhos sobre Estabilidade Global e Local das Estruturas de Edifícios. São Paulo, 1997.
  • 6
    VASCONCELOS, A. C. e FRANÇA, R. C. (1997) - Um método Simplificado e Muito Preciso para Avaliação dos Momentos de 2a Ordem em Edifícios Altos Usuais - Coletânea de Trabalhos sobre Estabilidade Global e Local das Estruturas de Edifícios. São Paulo, 1997.
  • 7
    CARMO, R.M.S (1995) - Efeitos de Segunda Ordem em Edifícios Usuais de Concreto Armado - Master's degree dissertation. College of Engineering of São Carlos, University of São Paulo. São Carlos - SP, 1995.

Publication Dates

  • Publication in this collection
    Apr 2016

History

  • Received
    10 June 2015
  • Accepted
    18 Sept 2015
IBRACON - Instituto Brasileiro do Concreto Instituto Brasileiro do Concreto (IBRACON), Av. Queiroz Filho, nº 1700 sala 407/408 Torre D, Villa Lobos Office Park, CEP 05319-000, São Paulo, SP - Brasil, Tel. (55 11) 3735-0202, Fax: (55 11) 3733-2190 - São Paulo - SP - Brazil
E-mail: arlene@ibracon.org.br