Abstract
The ultimate limit state of stability by equilibrium bifurcation, the limit states for stress and strain resulting from this condition were evaluated for an extremely slender real structure of reinforced concrete, with geometry varying along its length. The aspects related to nonlinearities of the material were considered through the recommendations on NBR 6118:2014, from the Brazilian Association of Technical Standards (ABNT). In the analytical solution, developed for stability analysis, all elements of the structural dynamics present in the system were taken into account, including the column selfweight. The critical buckling load was then dynamically defined to different instants of time. Reductions of 70% for the modulus of elasticity and 59% for the critical buckling force were found in analyses performed from zero and five thousand days. It was also possible to obtain the induced stresses on the homogenized crosssections and those transferred to reinforcement steel bars.
Keywords:
Limit states; critical buckling load; analytical solution; NBR 6118:2014
Graphical Abstract
Keywords:
Limit states; critical buckling load; analytical solution; NBR 6118:2014
1 INTRODUCTION
Limit state design is a concept for designing structures, based on a semiprobabilistic consideration, that has been used in place of the admissible stress method and has come to be incorporated by the normalization codes. The principle behind limit state design is to consider the possibility of simultaneous occurrence of action, admitting several combinations in which one of the variable actions is taken as the main and other as a secondary. Dead loads are always weighted with their values as permanent actions. As predicted by NBR 6118:2014 (ABNT, 2014), in the structural analysis, the influence of all actions that may have significant effects on the safety of the structure under consideration must be considered, taking into account the possible ultimate and service limit states. The limit state method is based on the assessment of certain parameters of design in relation to structural use and safety. Therefore, a limit state defines the impropriety for the use of the structure, for reasons of safety, functionality, aesthetics, or performance outside the standards specified for its normal use or interruption of operation due to the failure of one or more of its components. A threshold state separates a desired state from the adverse state (failure), as pointed out Geiker et al. (2019Geiker, M. R.; Michel, A.; Stang, H.; Lepech, M. D., (2019). Limit states for sustainable reinforced concrete structures. Cement Concrete Res. 112: 189195. DOI:10.1016/j.cemconres.2019.04.013
https://doi.org/10.1016/j.cemconres.2019...
).
The evaluation of failure in terms of the load capacity of a column can be accomplished by different methods. For slender elements, the load capacity by the socalled critical buckling load should be analyzed, and not only the capacity to resist stresses. This means a slender column reaches the Ultimate Limit State (ULS), defined by loss of stability, without reaching the maximum stress capacity of the crosssections, i.e., a condition of plastic behavior or proper failure. In case of concrete structures, the problem is aggravated by the possible occurrence of cracks, which decrease the crosssection inertia, and the rheological behavior of the material due to the creep, phenomena that are inherent and typical in brittle and viscoelastic materials, respectively.
The cracking, due to the tension stress induced on structural elements (Shaikh, 2018Shaikh, F. U. A., (2018). Effect of cracking on corrosion of steel in concrete. Int. J. Concr. Struct. Mater 12: 3. DOI:10.1186/s400690180234y.
https://doi.org/10.1186/s400690180234...
), is considered almost inevitable, even when under serviceability conditions, because of the low strength of concrete to this type of load, as described by NBR 6118:2014 (ABNT, 2014). For the creep specifically, which represents the gradual increase of strain over time, its consideration becomes necessary for stability evaluation of slender structures, since they can have their stiffness modified by the material’s own rheology (Metha et al., 1994Metha, P. K.; Monteiro, P. J. M, Filho, A. C., (1994). Concreto: estrutura, propriedades e materiais. Ed. PINI, São Paulo, Brasil.).
The first one to investigate the critical buckling load was Euler in 1774, defining formulations based on statics, and later having his studies complemented by Greenhill (1881Greenhill, A. G., (1881). Determination of the greatest height consistent with stability that a vertical pole or mast can be made, and the greatest height to which a tree of given proportions can grow. Proc. Cambridge Philos. Soc. 4: 6573.), who included the selfweight in the study of stability of poles and masts. According to Timoshenko and Gere (1961Timoshenko, S. P.; Gere, J. M., (1961). Theory of elastic stability, 2nd ed. McGrawHill Book Company, New York.), Euler buckling occurs from the region of elastic behavior of the material, which is defined as the event in which a structural element loses its stability or is bent under the action of an axial compressive force sufficiently large to take it off its straight initial configuration, this might be the most wellknown evaluation criterion for the safety of slender columns. The first postulations have been solved with statics by the equilibrium of forces that develop in the most stressed section of a column, after assuming some bent configuration. On the other hand, the buckling solution must be considered, at the same time, in the field of structural dynamics, which involves the concept of vibration of mechanical systems, for which slenderness and mass constitute the research base.
In the dynamics field of study, particular problems of mechanical vibrations are found, whether damped or not. The undamped free vibrations constitute the modal analysis. According to Shiki et al. (2014Shiki, S. B.; Lopes Jr., V.; da Silva, S., (2014). Identification of nonlinear structures using discretetime Volterra series. J. Braz. Soc. Mech. Sci. Eng., 36(3): 523532. DOI:10.1007/s4043001300889.
https://doi.org/10.1007/s404300130088...
), the modal analyses lead to important information on systems or structures, allowing to know specific characteristics of its dynamics, i.e., the natural frequencies and modes of vibration, having, thus, relevant practical application. Based on modal analysis concepts, Reis et al. (2018Reis, D. G. ; Siqueira, G. H. ; Vieira Jr, L. C. M. ; Ziemian, R. D. (2018). Simplified Approach Based on the Natural Period of Vibration for Considering SecondOrder Effects on Reinforced Concrete Frames. International Journal of Structural Stability and Dynamics. 18 (05): 18500741185007420. DOI: 10.1142/S0219455418500748
https://doi.org/10.1142/S021945541850074...
) and Leitão et al. (2019Leitão, F. F.; Siqueira, G. H.; Vieira Jr., L. C. M.; Almeida, S. J. C., (2019). Analysis of the global secondorder effects on irregular reinforced concrete structures using the natural period of vibration. Rev. IBRACON Estrut. Mater. 12(2): 408428. DOI:10.1590/s198341952019000200012.
https://doi.org/10.1590/s19834195201900...
) evaluated the fundamental period of reinforced concrete buildings, with regular and irregular geometry, to estimate the secondorder effects due to nonlinearities present in the system that reflect into columns of the building.
A reinforced concrete slender column may be described as a continuum system, therefore, with infinite degrees of freedom, submitted to axial compression loads, including its selfweight. This system, from the structural dynamics point of view, can be associated to another similar system, also a continuum, but possessing only one degree of freedom. This way, the first mode of vibration remains restricted to a previously set configuration for a mathematical function that describes the motion. Therefore, the stiffness and masses may be expressed as a function of generalized coordinates, previous and properly defined within the geometry of the system in a way that uniquely represents its vibration mode. Rayleigh (1877Rayleigh (1877). Theory of sound. Dover Publications, New York. Reissued in 1945.) applied this concept to the study of prismatic system vibration, establishing a valid and differentiable function throughout its domain. To find an analytical solution to present problem, a similar approach is taken in this paper.
However, for real structure cases, where the properties of the structural elements vary lengthwise, the developed formulation for the computation of stiffness and mass must be solved by addressing the defined intervals in the system geometry. For these cases, the integral equations obtained by the analytical solution must be solved within the established limits for each segment, in other words, the generalized properties must be solved for each structure segment, as defined in its geometry, and algebraically summed at the end. This procedure must be applied to solve the frequency as well as to determine the critical buckling load, and this analysis should also take into account all parts of the structural stiffness, such as: conventional stiffness, which depends on the material behavior regarding to elasticity, viscoelasticity or even plasticity (Babu et al., 2006Babu, R. R.; Benipal, G. S.; Singh, A. K., (2006). Plasticitybased constitutive model for concrete in stress space. Lat. Am. J. Solids Struct. 3: 417441.; Suter and Benipal, 2007Suter, M.; Benipal, G. S., (2007). Timedependent behaviour of reacting viscoelastic concrete. Lat. Am. J. Solids Struct. 4(2): 103120.); and geometric stiffness, which depends on the normal force acting on the structure, and which also should include the selfweight (Wahrhaftig et al., 2013Wahrhaftig, A. M.; Brasil, R. M. L. R. F.; Balthazar, J. M., (2013). The first frequency of cantilever bars with geometric effect: a mathematical and experimental evaluation. J. Braz. Soc. Mech. Sci. Eng.; 35(4): 457467. DOI:10.1007/s4043001300439.
https://doi.org/10.1007/s404300130043...
; Wahrhaftig and Brasil, 2016, 2017), as well as the soil contribution, given by the soilstructure interaction.
It is important to emphasize that the soilstructure interaction has been considered as an important matter when studying the vibratory motion in real structures (Zaman et al., 1984Zaman, M.; Desai, C. S.; Drumm, E. C., (1984). Interface model for dynamic soilstructure interaction. J. Geotech. Geoenviron. 110(9): 12571273. DOI:10.1061/(ASCE)07339410(1984)110:9(1257).
https://doi.org/10.1061/(ASCE)07339410(...
) since this interaction may significantly influence the behavior of dynamicsbased problems. For this reason, the foundations must be observed in the numerical simulations when elaborating mathematical models, as said by Yazdchi et al. (1999Yazdchi, M.; Khalili, N.; Valliappan, S., (1999). Dynamic soilstructure interaction analysis via coupled finiteelementboundaryelement method. Soil. Dyn. Earthq. Eng. 18(7): 499517. DOI:10.1016/S02677261(99)000196.
https://doi.org/10.1016/S02677261(99)00...
). In this sense, Rosa et al. (2018Rosa, L. M. P.; Danziger, B. R.; Carvalho, E. M. L., (2018). Soilstructure interaction analysis considering concrete creep and shrinkage. Rev. IBRACON Estrut. Mater. 11(3): 564585. DOI:10.1590/s198341952018000300008.
https://doi.org/10.1590/s19834195201800...
), when studying the combined effect of creep and shrinkage, and also of soilstructure interaction, in the evaluation of foundation settlements of a reinforced concrete building, could notice that the mathematical results obtained in the computerassisted simulation completed for that case were in favor of safety when compared to field measurements.
The methodology applied in this paper is aimed at the analysis of any structure that can be modeled as a unidimensional element, including buildings. However, the chosen structure for this study was a slender reinforced concrete (RC) pole, 46 meters long and with variable geometry, whose critical buckling load was dynamically defined. Therefore, the concrete creep and shirking effects were included according to recommendations from NBR 6118:2014 by the Brazilian Association of Technical Standards (ABNT), besides the geometric nonlinearity, which was solved by the geometric stiffness portion. At the end, it was possible to infer results regarding ultimate limit states of strains and calculation of stress on steel verifying the possibility of its yielding. It is worth mentioning that in the analysis of continuum elements, in accordance with this paper, when properly representing the concrete rheology and its interaction with the reinforcement bars, Section 14.2.4 of NBR 6118:2014 (ABNT, 2014) may be used to assess the structural performance in the serviceability case or in failure.
2 SOLUTION OF THE DYNAMICS FOR THE ANALYSIS OF STABILITY ULTIMATE LIMIT STATE BY EQUILIBRUIM BIFURCATION (BUCKLING)
The model presented in Figure 1 represents a column fixed at the bottom and free at the top end, in undamped free vibration, approximating the mathematical model of the structural dynamics field to the practical problem proposed, using modal analysis, where t indicates time, ϕ(x) represents the shape function of the oscillatory motion, L is the total length (or height) of the structure, L _{s} and L _{s−1} are the height at the upper and lower limits of a given segment s, whose length is obtained by the difference between these two positions, and v(t) is the generalized coordinate localized at the free extremity of the column. Equation (1) refers to a trigonometric function taken as a shape function, considered valid at any point in the structural domain and that obeys the boundary conditions of the problem:
where x starts at the base. This model represents a column under axial compressive load N(x) that can have constant or variable properties along the length. These properties include the geometry, elasticity or viscoelasticity and density. Springs k _{s}(x) are applied to represent the lateral soilstructure interaction.
With these conditions, the column is under the influence of gravitational forces, originating from the distributed mass due to the selfweight of the structure and other added masses; and of a concentrated mass at the free top end, to be defined in the event of loss of stability, whose value represents the vertical maximum load that can be applied.
To get the analytical solution of this problem, it was necessary to consider a mathematical function to reduce the problem to a system with only one degree of freedom. For the present case, a trigonometric function given by Equation (1) was adopted. The use of Equation (1) as a valid shape function for a structure with variable geometry was validated by Wahrhaftig (2019aWahrhaftig, A. M., (2019a). Analysis of the first modal shape using two case studies. Int. J. Comput. Meth. 16(6): 1840019. DOI:10.1142/S0219876218400194.
https://doi.org/10.1142/S021987621840019...
) in comparison to computational analyses using the finiteelement method (FEM). The FEM is a technique to discretize the continuum that has been used recently in computer modeling to analyze buckling and vibration of reinforced concrete structures, as mentioned by Rodrigues et al. (2014Rodrigues, E.A.; Manzoli, O.L.; Bitencourt Jr, L.A.G.; Prazeres, P.G.C.; Bittencourt, T.N., (2014). Failure behavior modeling of slender reinforced concrete columns subjected to eccentric load. Lat. Am. J Solids Struct. 12: 520541. DOI:10.1590/167978251224.
https://doi.org/10.1590/167978251224...
), who studied concrete columns eccentrically loaded, considering nonlinearities inherent in the system, including creep and shrinkage, for sections with reinforcement ratio of 1.98% and 3.95%, by Bueno and Loriggio (2016Bueno, J. R.; Loriggio, D. D., (2016). Influence of the flexibility of beams and slabs in static response and dynamic properties. Rev. IBRACON Estrut. Mater. 9(6): 842855. DOI:10.1590/s198341952016000600003.
https://doi.org/10.1590/s19834195201600...
), who evaluated undamped free vibrations of reinforced concrete decks in relation to recommended limit states, and also by Rajasekaran and Khaniki (2018Rajasekaran, S.; Khaniki, H. B., (2018). Bending, buckling and vibration analysis of functionally graded nonuniform nanobeams via finite element method. J. Braz. Soc. Mech. Sci. Eng. 40: 549. DOI:10.1007/s4043001814606.
https://doi.org/10.1007/s404300181460...
), who assessed buckling of beams with functional graded materials.
Applying the principle of virtual work and its derivatives, similarly to that done by Clough and Penzien (1995Clough, R. W.; Penzien, J., (1995). Dynamics of Structures. Computers & Structures, Inc., Berkeley.), the dynamic properties of interest are obtained. The conventional elastic/viscoelastic stiffness is specified as:
with
where, for the segment s of the structure, E _{s}(t) is the viscoelastic module of the material as a function of time, I _{s}(x) is the moment of inertia that varies along the segment about the considered motion, obtained by interpolation of subsequent sections, all homogenized (if constant, it is simply going to be I _{s}), k _{0s}(t) is the term for stiffness variation over time, K _{0}(t) is the final stiffness varying with time, and n is the total number of segments in the analyzed geometry. If the material exclusively shows elastic behavior, the respective equations will not be functions of time anymore.
The geometric stiffness presents as a function of the axial force, including the selfweight contribution, being calculated as follows:
with “j” being just a counter for the summation,
and
where k _{gs}(m _{0}) is the geometric stiffness of segment s, N _{0}(m _{0}) is the concentrated force at the top part of the system, calculated by Equation (5), and K _{g}(m _{0}) is the total geometric stiffness of the structure, given by Equation (6); as can be seen, all of them depend on the mass m _{0} located at the free end of the structure. In Equation (4), N _{s} represents the normal force in the upper segments above the considered segment, which can be obtained from:
where
considering:
with
and
where A _{s}(x) represents the crosssection area and ρ _{s} is the material density for the respective segment s. If the crosssection has a constant area along the interval, A _{s}(x) will be A _{s} , and consequently the mass distribution per unit length will also be constant. In the same way, if the mass m _{0} does not vary, all of its dependent parameters will also be constant.
In order to consider the soil influence on the system vibration, it was necessary to consider it as a series of vertically distributed springs along the foundation that act to introduce a restoring force into the structural system. With k _{s} (x) denoting the spring parameter, the effective stiffness of the soil, as a function of the variable x along the length, can be defined as:
with
and
where the parameter k _{Soil} is an elastic property that consists of the sum of k _{s} (x) along the depth of the foundation, which depends on its geometry D _{s} (x) and on the elastic parameter of the soil S _{s} , considered constant, in this case, on each layer, although it may not necessarily be constant. Considering as positive a compressive normal force, it is possible to obtain the total stiffness of the system as a function with two variables:
Accordingly, the frequency can be calculated, in hertz, with Equation (16), as a function of time and concentrated mass at the free end.
The mathematical method previously described defines the critical buckling load for the nullity of the frequency, representing the occasion at which the structure loses its stiffness. All generalized parameters, such as total stiffness, Equation (15), total mass, Equation (8), as well as the normal force at the free end, Equation (6), were expressed as a function of the concentrated mass at the top part. When the creep is introduced into the analysis, it will be necessary to consider a model that represent the viscoelastic behavior of concrete, and after that, the frequency becomes a timedependent function, since the modulus of elasticity is also timedependent. Thus, by writing the frequency in terms of time and concentrated mass at the free end, Equation (16), the critical buckling load is established for a determined value that makes the frequency becomes null, at any instant of interest, that multiplied by the gravitational acceleration, will define the force corresponding to the collapse of the system.
Taking all the previous considerations into account, and assuming the concentrated mass as an independent variable of the problem, the critical buckling load N
_{bk} can be determined by using the mathematical concept in Equation (17), and the details of the present analytical process can be found in Wahrhaftig et al. (2019bWahrhaftig, A. M.; Silva, M. A.; Brasil, R. M. L. R. F., (2019b). Analytical determination of the vibration frequencies and buckling loads of slender reinforced concrete towers. Lat. Am. J. Solids Struct. 16(5): e196. DOI:10.1590/167978255374.
https://doi.org/10.1590/167978255374...
).
3 CREEP CRITERIA ACCORDING TO THE BRAZILIAN CODE
The model proposed by NBR 6118:2014 (ABNT, 2014) considers the creep strain ε _{cc} to be composed by two terms, one related to the fast strain ε _{cca} and the other one to the slow strain ε _{ccb} . The fast strain is irreversible and happens during the first twentyfour hours after the structure is submitted to loading. On the other hand, the slow strain has two other parts, the slow irreversible strain ε _{ccf} and the slow reversible strain ε _{ccd} . Equation (18) shows the terms used in the determination of concrete creep and the slow strain,
and Equation (19) represents the total strain of concrete ε _{c,tot} which includes an immediate one ε _{c} .
Since the slow strain of concrete ε _{cc} can be expressed in terms of the creep strain coefficient φ, Equation (19) can be replaced by Equation (20),
in which the creep strain coefficient φ is decomposed into the fast strain φ _{a} , the slow irreversible strain φ _{f} , and the slow reversible strain φ _{d} , according to Equation (21),
The fast strain is calculated by the ratio between the concrete strength at the time of the load application and its final strength. The slow irreversible strain is determined as a function of the relative humidity, concrete consistency at placing, assumed thickness of the member, and assumed age of the concrete at the moment of loading. The slow reversible strain depends only on the duration of loading. Its final value and development over time do not depend on the age of concrete at the instant of the load application.
NBR 6118:2014 (ABNT, 2014) introduces some hypotheses that are adopted to define the creep effects when the stresses in the concrete are related to normal service loads. The first one is the consideration of linear variation between the applied stress and the creep strain ε _{cc} . It is also considered that, for the same concrete, the strain curves of the irreversible slow strain φ _{f} as a function of time, and corresponding to different ages of the concrete at the moment of the loading, are obtained one in relation to each other by parallel dislocation of the strain axis, as shown in Figure 2.
NBR 6118:2014 (ABNT, 2014) preconizes that when the determination of the concrete creep is desired to be done for a certain instant t, Equation (18) can be written as a function of time, as can be seen in Equation (22), from which an also timedependent function of the modulus of elasticity is originated, in terms of the creep coefficient φ(t,t _{0} ):
where ε _{cc} (t,t _{0} ) represents the variation over time of the strain due to creep, σ _{c} is the stress in the concrete, E _{c28} is the initial tangent modulus at twentyeight days, and φ(t,t _{0} ) is the creep coefficient as a function of time, given by Equation (23):
where t is the assumed concrete age for the considered instant, in days; t _{0} is the assumed concrete age at the moment of loading, also in days, with φ _{a} calculated by Equation (24); the coefficient φ _{f∞} is the final value for the slow strain given by Equation (25); β _{f} (t) and β _{f} (t _{0} ) are coefficients related to the slow irreversible strain, as a function of the age of the concrete obtained by Equation (26); φ _{d∞} is the final value for the coefficient of the slow reversible strain that is considered equal to 0.4; and β _{d} is the coefficient related to the slow reversible strain, calculated as a function of time after loading (t  t _{0} ), expressed by Equation (27):
Note that the relation between f _{c} (t _{0} ) and f _{c} (t _{∞} ) represents the strength gain curve of concrete; φ _{1c} is the coefficient that depends on the relative humidity, in percentage, and on the consistency of the concrete, and φ _{2c} is the coefficient depending on the assumed thickness of the member, defined by Equation (28); and the A, B, C and D terms are related to the coefficient for slow irreversible strain and can be calculated by the following Equations (29), (30), (31) and (32), respectively:
where h _{fic} and h denote the assumed thickness of the analyzed member, in centimeters and meters, respectively. Therefore, the timedependent function referred to as the modulus of elasticity can be described by Equation (33):
where E _{c0} is the modulus of elasticity of concrete at the instant of loading.
4 EVALUATION OF A CASE STUDY
The analyzed structure in this evaluation is a real pole, extremely slender, made from reinforced concrete, with variable geometry, shown in Figure 3, that is 46 m high, including a 40 m superstructure, with hollow circular crosssection, and the 6 m deep foundation, which is relatively deep, is a belled shaft foundation. The dimensions, heights and structural arrangement of the column and column segments are shown in Figure 4, (a) and (b), where g is the gravitational acceleration; S1, S2, S3, S4, and S5 are the crosssections; D and th indicate the external diameter and wall thickness, respectively; d _{b} represents the diameter of steel bars; n _{b} is the number of bars, c´ is the cover of concrete in the respective crosssections, and “Var” indicates tapered section at the segment.
It is worth mentioning that, although the average wall thickness of the crosssections represents 0.72% of the height of the superstructure, no analyses regarding the local buckling or distortion of the section were considered, since it is a problem that should be treated with another specific approach, as discussed by Pierin and Silva (2015Pierin, I.; Silva, V. P., (2015). Distortional buckling resistance of coldformed steel. J. Braz. Soc. Mech. Sci. Eng. 37: 11631171. DOI:10.1007/s404300140252x.
https://doi.org/10.1007/s404300140252...
). Therefore, the present approximation is restricted to the global instability evaluation, similar to the work developed by Schnabl and Planinc (2017Schnabl, S.; Planinc, I., (2017). Buckling of slender concretefilled steel tubes with compliant interfaces. Lat. Am. J. Solids Struct. 14(10): 18371852. DOI:10.1590/167978253479.
https://doi.org/10.1590/167978253479...
). The moduli of elasticity of concrete at 28 days for the superstructure and foundation, solved with Equation (34) and (35), with f
_{ck} in MPa for both Equations, are, respectively, 34278.92 MPa and 21287.37 MPa, considering nominal strengths f
_{ck} equal to 45 MPa and 20 MPa:
Usually, equipment is attached at the top of the structure, settling a concentrated mass applied at the free end of the column, in which the limit value in relation to the stability loss due to buckling should be determined. Additional devices are also installed along all the structure, adding a distributed mass of 40 kg/m.
The foundation is a belled shaft, with a bell diameter of 140 cm and length of 20 cm; the shaft is 80 cm in diameter and 580 cm long. The lateral soil pressure was represented by an elastic parameter equal to 2669 kN/m^{3}. The physical nonlinearity of the RC was computed following the recommendations from NBR 6118:2014 (ABNT, 2014), which suggests a 50% reduction in the section moment of inertia for analyses under similar circumstances, regardless of the possibility to choose another reduction factor, as discussed by Marin and El Debs (2012Marin, M. C.; El Debs, M. K., (2012). Contribution to assessing the stiffness reduction of structural elements in the global stability analysis of precast concrete multistorey buildings. Rev. IBRACON Estrut. Mater. 5(3): 316342. DOI:10.1590/S198341952012000300005.
https://doi.org/10.1590/S19834195201200...
), but also considering the fact that columns in buildings normally have the same crosssection between floors. The density of reinforced concrete was assumed to be 2600 kg/m^{3} for the superstructure and 2500 kg/m^{3} for the foundation.
To take the steel bars into consideration in the calculation of the crosssection moment of inertia, homogenization factors were used, multiplying the nominal moment of inertia in each section. Using the parallel axis theorem, the values were found to be 1.026, 1.074, 1.050, 1.061, and 1.064 for the sections 1 to 5, respectively, at the commence of loading. The creep due to the viscoelastic behavior of concrete was accounted for in the structure as recommended by NBR 6118:2014 (ABNT, 2014). Studies regarding creep effects transferred to the steel in reinforced concrete columns can be found in Kataoka and Bittencourt (2014Kataoka L. T.; Bittencourt T. N., (2014). Numerical and experimental analysis of timedependent load transfer in reinforced concrete columns. Rev. IBRACON Estrut. Mater. 7(5): 747774. DOI:10.1590/S198341952014000500003.
https://doi.org/10.1590/S19834195201400...
) and Madureira et al. (2013Madureira, E. L.; Siqueira, T. M.; Rodrigues, E. C., (2013). Creep strains on reinforced concrete columns. Rev. IBRACON Estrut. Mater. 6(4): 537560. DOI:10.1590/S198341952013000400003.
https://doi.org/10.1590/S19834195201300...
). Despite the description of the stress transfer to reinforcement steel in the study by Kataoka and Bittencourt (2014), made with short specimens, the numerical evaluation was found to be inconclusive in the study. In their first work, Madureira et al. (2013) describe the transfer mechanism in which forces are carried to the reinforcement steel due to creep of concrete and also indicate that there may be a convergence of the strains after 3000 days of loading. In a second investigation, Madureira and Paiva (2018) discuss the stress transfer to the reinforcement steel due to the creep of concrete in bent RC members.
5 RESULTS AND DISCUSSION
The variation of the critical buckling load and modulus of elasticity over time, as summarized by Equations (17) and (33), can be seen in the graphs in Figure 5, (a) and (b), respectively. The latter provides important information for temporary loading and construction with time of use less than one of the mentioned moments. In order to calculate the variation of the critical load over time, a small computer program was developed, considering a mass increment of 0.1 kg, a time frame of 100 days, and a frequency with 10^{−3} Hz precision.
Besides the arbitrarily chosen intervals, the analysis could be done for any other preferred period of time during the structure’s lifespan, although it is possible to notice that the convergence tends to occur after 5000 days from the beginning of the structure performance, as can be observed in Figure 6(a). Setting the normal force as the same value of the critical load, Figure 6(b) allows the assessment of the variation of the structural frequency over the period considered. The strain analysis in the ultimate limit state (ULS) of a compressed element can be executed by examining the graphs in Figure 7, where it is possible to verify that, at the instants of time defined for the evaluation, the strains are below the limit of 0.002 (0.20%) set by NBR 6118:14 (ABNT, 2014), the value that denotes when the cracks start to occur in the concrete in compression (Tamayo et al.,2015Tamayo, J.L.P.; Morsch, I.B.; Awruch, A.M., (2015). Shorttime numerical analysis of steelconcrete composite beams. J. Braz. Soc. Mech. Sci. Eng. 37(4): 10971109. DOI:10.1007/s4043001402379.
https://doi.org/10.1007/s404300140237...
). The maximum strain value found under the critical buckling load at the most unfavorable instant was 0.12‰, which represents 6.07% of the limit previously set.
These results were obtained by taking into account the total strain of the structure and by computing the strength gain of the concrete as well as the homogenization of the sections, given by Equations (36), (37) and (39). Equation 36 is modified from the original equation proposed in the Eurocode (EN 199211:2004) to estimate the compressive strength of concrete at various ages, introducing the factor of safety γ _{c} , equal to 1.4, according NBR 6118:2014 (ABNT, 2014).
and
with t in day.
The results presented in Figure 7 are compatible with the findings of Fischer et al. (2014Fischer, I.; Pichler, B.; Lach, E.; Terner, C., Barraud, E.; Britz, F., (2014). Compressive strength of cement paste as a function of loading rate: Experiments and engineering mechanics analysis. Cem. Concr. Res. 58: 186200. DOI:10.1016/j.cemconres.2014.01.013.
https://doi.org/10.1016/j.cemconres.2014...
), who experimentally analyzed standard specimens of concrete and verified the existence a linear response of creep to strain below 0.4. The relative stresses, expressed by the ratio between the induced stress in the homogenized section and the strength of concrete, as a function of the force at the top of the structure, can be seen in the graphs on the left in Figure 8, revealing results less than a unit. In fact, the relative stresses decrease over time due to force transfer to the reinforcement steel, influenced by the increase in the area homogenization factor, Equation (40), whose development is presented on the right side in Figure 8. When the normal force at the top equals the critical buckling load, the values found for the relative stresses in sections 3 to 5 were 4.91%, 5.18%, and 5.07%, respectively.
It is important to assess the reduced normal forces ((t) at the crosssections of the structure when the force at the column freeend approaches the value of the critical buckling load, at the respective chosen instants of time. In order to do that, both the homogenization of the section as well as the strength gain of concrete over time were taken into account, using Equation (38),
where f _{c} (t) is the strength of concrete, given by Equation (36), and A _{hom} (t) is the homogenization of the concrete area, obtained by Equation (39)
with
in which A _{cn} is the net area of concrete, and A _{st} and E _{st} are the area of steel and the steel modulus of elasticity, respectively. With the decrease in the concrete modulus of elasticity, and consequent increase in the homogenization factor of the sections ζ(t), it is possible to verify that the internal forces, calculated by equilibrium of forces and compatibility of displacements, occur for an additional part of the loading taken by the reinforcement, compensating for the increase in strain due to creep, an effect that increases together with the reinforcement ratio.
Table 1 summarizes the obtained results, in which the reinforcement ratios (w) for the sections 3 to 5 are presented. The following considerations were made in the analyses of this paper: modulus of elasticity of reinforcement steel equal to 205 GPa, concrete production under standard conditions, relative humidity of 70%, and gravitational acceleration of 9.807 m/s^{2}.
Setting the buckling load as the loading limit that can be applied to each section, with σ(m _{0} ) as the stress acting in the section as a function of the normal top force and σ _{u}(t) = N _{bk} (t)/A _{hom} (t), the normalized stress σ(m _{0} )/σ _{u}(t) becomes a bilinear function composed of inclined and vertical segments, as shown in Figure 9, where the curve’s inflection point at each instant corresponds to the critical buckling load (CBL), from which the structure finds itself on the verge of collapse and any minimal load no longer can be applied. The decreasing of the normalized stress with the time rises with the increasing of the homogenized area of concrete, due the growing of ζ(t) in Equation (40). That means the creep decreases the stress on the sections in terms of the ULS (ultimate limit state) of the capacity resistance under vertical force, Figure 9(left), even though it lifts the strain of the structure, Figure 9(right). In the latter aspect, it is possible to see that the selfweight of the structure imposes an initial strain on the system, and strains are not starting from zero even with no single lumped mass applied.
Normalized stress on the homogenized concrete section in ULS (CBL = Critical buckling load).
As a result of the concrete creep, there is an axial stress transfer from the concrete to the reinforcing bars, without reaching the yield stress of steel, though. In fact, as displayed in Figure 10, the stress in the reinforcement shows a linear growth, reaching the maximum value at the critical buckling load. Calculations that include the weight of upper segments led to the following results: stress of 37.47 MPa (7.49% of f _{y} ) at the critical load for section S3; 42.40 MPa (8.48% of f _{y} ) for section S4; and 30.91 MPa (6.18% of f _{y} ) for section S5, all with steel yield strength f _{y} = 500 MPa. Obviously, the interpretation of these graphs implies understanding their validity until the critical buckling load. Even when admitting a large spectrum of values of the normal force at the top, that overcome the limit value stresses, the reinforcement steel remains in the elastic zone and much below to yield.
5 CONCLUSIONS
The critical buckling load for a real reinforced concrete structure has been determined using an analytical method based on vibration concepts for structural systems. The analyzed structure presented geometric variations along its height, and all necessary parameters for calculations of the dynamics were considered, such as: geometric imperfections, considered by the geometric stiffness; the physical nonlinearity and creep of concrete, that were taken into account according to the recommendations from NBR 6118:2014 by ABNT.
Using the previous considerations in the analytical procedure, an axial compressive force of 252.24 kN at the initial instant and 104.51 kN after 5000 days from the beginning of the structure performance were defined for the buckling of the structure, with both values corresponding to a null first natural vibration frequency, which represents a 59% reduction in the structure’s vertical strength. Even with this reduction, the maximum buckling load the structure can withstand is higher than the loading usually applied to this kind of system. Important to note that no one factor of safety was assigned to the applied normal forces. It is important to clarify that according NBR 6118:2014 (ABNT, 2014) application of a safety coefficient of 1.4 for dead load is waited, being also necessary performing a verification for local second order effects when the slenderness ratio of a structural element is higher than 200.
Due to the creep, the modulus of elasticity exhibited a 70% variation during the considered time frame. The evaluation of induced stresses in the homogenized concrete section revealed a maximum percentage of 5.2% among those presented by different critical buckling loads. For the reinforcement, the maximum relative value found for the stress in the steel, with load transfer due to creep, was 8.5% of the material yield stress (500 MPa), ruling out the possibility of the steel to yield. The strain calculated at the most unfavorable conditions exhibited only 6.1% of the recommended limit. With these results, it has been proved that the critical buckling load determines not only the ultimate limit state to the loss of stability as a deformable solid but also the ultimate limit state for stress and stain for slender structures.
Lastly, it should be noted that the analytical method presented in this paper is free from the application of iterative methods or discretization processes, with a solution directly obtained in the continuum. This method may be adapted to include different geometries, with no, one or some concentrated masses (and forces) along the height, and, still, making it able to solve any problem that may be modeled as unidimensional elements, including buildings.
A comparative analysis between other calculation methods and design criteria from different design codes should be included in future studies, as well as the influences of mechanical and thermalrelated strains, and also the consideration of other boundary conditions.
References
 Associação Brasileira de Normas Técnicas (2014). NBR 6118: Projeto de Estruturas de Concreto  Procedimento. Rio de Janeiro, 2014.
 Babu, R. R.; Benipal, G. S.; Singh, A. K., (2006). Plasticitybased constitutive model for concrete in stress space. Lat. Am. J. Solids Struct. 3: 417441.
 Bueno, J. R.; Loriggio, D. D., (2016). Influence of the flexibility of beams and slabs in static response and dynamic properties. Rev. IBRACON Estrut. Mater. 9(6): 842855. DOI:10.1590/s198341952016000600003.
» https://doi.org/10.1590/s198341952016000600003  Clough, R. W.; Penzien, J., (1995). Dynamics of Structures. Computers & Structures, Inc., Berkeley.
 European Committee for Standardization (2004). EN 199211  Eurocode 2: Design of concrete structures  Part 1: General rules and rules for buildings. Brussels, 2004.
 Fischer, I.; Pichler, B.; Lach, E.; Terner, C., Barraud, E.; Britz, F., (2014). Compressive strength of cement paste as a function of loading rate: Experiments and engineering mechanics analysis. Cem. Concr. Res. 58: 186200. DOI:10.1016/j.cemconres.2014.01.013.
» https://doi.org/10.1016/j.cemconres.2014.01.013  Geiker, M. R.; Michel, A.; Stang, H.; Lepech, M. D., (2019). Limit states for sustainable reinforced concrete structures. Cement Concrete Res. 112: 189195. DOI:10.1016/j.cemconres.2019.04.013
» https://doi.org/10.1016/j.cemconres.2019.04.013  Greenhill, A. G., (1881). Determination of the greatest height consistent with stability that a vertical pole or mast can be made, and the greatest height to which a tree of given proportions can grow. Proc. Cambridge Philos. Soc. 4: 6573.
 Kataoka L. T.; Bittencourt T. N., (2014). Numerical and experimental analysis of timedependent load transfer in reinforced concrete columns. Rev. IBRACON Estrut. Mater. 7(5): 747774. DOI:10.1590/S198341952014000500003.
» https://doi.org/10.1590/S198341952014000500003  Leitão, F. F.; Siqueira, G. H.; Vieira Jr., L. C. M.; Almeida, S. J. C., (2019). Analysis of the global secondorder effects on irregular reinforced concrete structures using the natural period of vibration. Rev. IBRACON Estrut. Mater. 12(2): 408428. DOI:10.1590/s198341952019000200012.
» https://doi.org/10.1590/s198341952019000200012  Madureira, E. L.; Paiva, L. A., (2018). Fresh concrete consistency effect on thinwalled columns creep phenomenon. Rev. IBRACON Estrut. Mater. 11(3): 644651. DOI:10.1590/s198341952018000300012.
» https://doi.org/10.1590/s198341952018000300012  Madureira, E. L.; Siqueira, T. M.; Rodrigues, E. C., (2013). Creep strains on reinforced concrete columns. Rev. IBRACON Estrut. Mater. 6(4): 537560. DOI:10.1590/S198341952013000400003.
» https://doi.org/10.1590/S198341952013000400003  Marin, M. C.; El Debs, M. K., (2012). Contribution to assessing the stiffness reduction of structural elements in the global stability analysis of precast concrete multistorey buildings. Rev. IBRACON Estrut. Mater. 5(3): 316342. DOI:10.1590/S198341952012000300005.
» https://doi.org/10.1590/S198341952012000300005  Metha, P. K.; Monteiro, P. J. M, Filho, A. C., (1994). Concreto: estrutura, propriedades e materiais. Ed. PINI, São Paulo, Brasil.
 NBR 6118:2014 Design of structural concrete  procedure. (2014). Brazilian Association for Standardization (ABNT). Rio de Janeiro, Brazil.
 Pierin, I.; Silva, V. P., (2015). Distortional buckling resistance of coldformed steel. J. Braz. Soc. Mech. Sci. Eng. 37: 11631171. DOI:10.1007/s404300140252x.
» https://doi.org/10.1007/s404300140252x  Rajasekaran, S.; Khaniki, H. B., (2018). Bending, buckling and vibration analysis of functionally graded nonuniform nanobeams via finite element method. J. Braz. Soc. Mech. Sci. Eng. 40: 549. DOI:10.1007/s4043001814606.
» https://doi.org/10.1007/s4043001814606  Rayleigh (1877). Theory of sound. Dover Publications, New York. Reissued in 1945.
 Reis, D. G. ; Siqueira, G. H. ; Vieira Jr, L. C. M. ; Ziemian, R. D. (2018). Simplified Approach Based on the Natural Period of Vibration for Considering SecondOrder Effects on Reinforced Concrete Frames. International Journal of Structural Stability and Dynamics. 18 (05): 18500741185007420. DOI: 10.1142/S0219455418500748
» https://doi.org/10.1142/S0219455418500748  Rodrigues, E.A.; Manzoli, O.L.; Bitencourt Jr, L.A.G.; Prazeres, P.G.C.; Bittencourt, T.N., (2014). Failure behavior modeling of slender reinforced concrete columns subjected to eccentric load. Lat. Am. J Solids Struct. 12: 520541. DOI:10.1590/167978251224.
» https://doi.org/10.1590/167978251224  Rosa, L. M. P.; Danziger, B. R.; Carvalho, E. M. L., (2018). Soilstructure interaction analysis considering concrete creep and shrinkage. Rev. IBRACON Estrut. Mater. 11(3): 564585. DOI:10.1590/s198341952018000300008.
» https://doi.org/10.1590/s198341952018000300008  Schnabl, S.; Planinc, I., (2017). Buckling of slender concretefilled steel tubes with compliant interfaces. Lat. Am. J. Solids Struct. 14(10): 18371852. DOI:10.1590/167978253479.
» https://doi.org/10.1590/167978253479  Shaikh, F. U. A., (2018). Effect of cracking on corrosion of steel in concrete. Int. J. Concr. Struct. Mater 12: 3. DOI:10.1186/s400690180234y.
» https://doi.org/10.1186/s400690180234y  Shiki, S. B.; Lopes Jr., V.; da Silva, S., (2014). Identification of nonlinear structures using discretetime Volterra series. J. Braz. Soc. Mech. Sci. Eng., 36(3): 523532. DOI:10.1007/s4043001300889.
» https://doi.org/10.1007/s4043001300889  Suter, M.; Benipal, G. S., (2007). Timedependent behaviour of reacting viscoelastic concrete. Lat. Am. J. Solids Struct. 4(2): 103120.
 Tamayo, J.L.P.; Morsch, I.B.; Awruch, A.M., (2015). Shorttime numerical analysis of steelconcrete composite beams. J. Braz. Soc. Mech. Sci. Eng. 37(4): 10971109. DOI:10.1007/s4043001402379.
» https://doi.org/10.1007/s4043001402379  Timoshenko, S. P.; Gere, J. M., (1961). Theory of elastic stability, 2nd ed. McGrawHill Book Company, New York.
 Wahrhaftig, A. M., (2019a). Analysis of the first modal shape using two case studies. Int. J. Comput. Meth. 16(6): 1840019. DOI:10.1142/S0219876218400194.
» https://doi.org/10.1142/S0219876218400194  Wahrhaftig, A. M.; Brasil, R. M. L. R. F., (2016). Vibration analysis of mobile phone mast system by Rayleigh method. Appl. Math. Modell. 42: 330345. DOI:10.1016/j.apm.2016.10.020.
» https://doi.org/10.1016/j.apm.2016.10.020  Wahrhaftig, A. M.; Brasil, R. M. L. R. F., (2017). Initial undamped resonant frequency of slender structures considering nonlinear geometric effects: the case of a 60.8 mhigh mobile phone mast. J. Braz. Soc. Mech. Sci. Eng. 39(3): 725735. DOI:10.1007/s4043001605471.
» https://doi.org/10.1007/s4043001605471  Wahrhaftig, A. M.; Brasil, R. M. L. R. F.; Balthazar, J. M., (2013). The first frequency of cantilever bars with geometric effect: a mathematical and experimental evaluation. J. Braz. Soc. Mech. Sci. Eng.; 35(4): 457467. DOI:10.1007/s4043001300439.
» https://doi.org/10.1007/s4043001300439  Wahrhaftig, A. M.; Silva, M. A.; Brasil, R. M. L. R. F., (2019b). Analytical determination of the vibration frequencies and buckling loads of slender reinforced concrete towers. Lat. Am. J. Solids Struct. 16(5): e196. DOI:10.1590/167978255374.
» https://doi.org/10.1590/167978255374  Yazdchi, M.; Khalili, N.; Valliappan, S., (1999). Dynamic soilstructure interaction analysis via coupled finiteelementboundaryelement method. Soil. Dyn. Earthq. Eng. 18(7): 499517. DOI:10.1016/S02677261(99)000196.
» https://doi.org/10.1016/S02677261(99)000196  Zaman, M.; Desai, C. S.; Drumm, E. C., (1984). Interface model for dynamic soilstructure interaction. J. Geotech. Geoenviron. 110(9): 12571273. DOI:10.1061/(ASCE)07339410(1984)110:9(1257).
» https://doi.org/10.1061/(ASCE)07339410(1984)110:9(1257)
Publication Dates

Publication in this collection
03 Feb 2020 
Date of issue
2020
History

Received
29 Aug 2019 
Reviewed
22 Nov 2019 
Accepted
17 Dec 2019 
Published
07 Jan 2020