Available online at www.sciencedirect.com
ScienceDirect APCBEE Procedia 9 (2014) 339 – 346
ICCEN 2013: December 13-14, Stockholm, Sweden
A Simple Model for Calculating the Fundamental Period of Vibration in Steel Structures Dia Eddin Nassania,* a
Hasan Kalyoncu University, Havalimani yolu - Şahinbey, Gaziantep27410 , Turkey
Abstract Most international codes allow the use of an equivalent static lateral forces method for practical design of real structures to withstand earthquake actions. This method requires the calculation of an approximation for the fundamental mode period of vibrations of these structures and several empirical equations have been suggested. Most of these equations are not considering the effect of semi-rigid connections and provide only a crude approximation to the required fundamental periods. In this paper, a new formula to calculate a satisfactory approximation of fundamental period of the structure is proposed. This formula takes into account the effect of semi-rigid connections, mass and the stiffness of the structure. In order to verify the accuracy of the proposed formula, several examples are presented in which the fundamental mode periods of several structures have been calculated utilizing the proposed formula and the conventional empirical equations. Comparing the obtained results with those obtained from a well known computer program STAAD.PRo has shown that the proposed formula provides a more accurate estimation of the fundamental period of the structures. © 2013 2014 Dia Eddin Nassani. Published by Elsevier B.V. © Published by Elsevier B.V. Selection and/or peer review under responsibility of Asia-Pacific Selection andBiological peer review& under responsibility Engineering of Asia-Pacific Society Chemical, Biological & Environmental Engineering Society Chemical, Environmental Keywords: STEEL STRUCTURE, CONNECTION, FUNDAMENTAL PERIOD
1. Introduction Most structures nowadays are designed to withstand horizontal loading resulting from earthquakes that might strike during the lifetime of the structures. Calculating an approximate estimation of loading resulting from earthquake actions requires the
* Corresponding author. Tel.: +90-342-2118080-1218. E-mail address:
[email protected].
2212-6708 © 2014 Dia Eddin Nassani. Published by Elsevier B.V. Selection and peer review under responsibility of Asia-Pacific Chemical, Biological & Environmental Engineering Society doi:10.1016/j.apcbee.2014.01.060
340
Dia Eddin Nassani / APCBEE Procedia 9 (2014) 339 – 346
calculation of fundamental mode period of vibrations. Since most design methods, recommended by international codes [1- 4], allows the use of the "equivalent static lateral forces method" and this method takes into consideration only the fundamental mode period of vibration to calculate a set of equivalent horizontal static forces to be applied to the structure. The equivalent static lateral forces method is applicable only to regular structures that can be analyzed in two planes and whose responses are not highly dependent on higher modes of vibrations. Several empirical equations have been proposed to calculate an estimation to the fundamental mode period of vibrations of the structures such as the following equations adopted by ASCE 7-05 [1]: Ta
0.0724 hn0.8 for steel moment-resisting frames
(1)
Ta
0.0466 hn0.9 for concrete moment-resisting frames
(2)
Ta
0.0731 hn
0.75
for eccentrically braced steel frame
(3)
Where hn is the total structure height (taken in m) And there are other equations adopted by EC8 [4]: T1
0.075 h t
for framed structures
0.75
Where ht is the total structure height (taken in m) 0.09 h t for framed and wall concrete structures T1 L0.5 Where ht , L are building height and width parallel to the applied forces (in m) T1
§1 2S ¨ ¨g ©
¦w G ¦H G i
2
i
i
i
· ¸ ¸ ¹
(4) (5)
0.5
for all types of structures
(6)
Where wi is the total gravity load of the storey i and δi is the lateral displacement of storey i under the effect of an arbitrary set of horizontal loads Hi applied at each floor level. Equation 6 is known as the Rayleigh equation. Although equations 1, 2, 3, 4 and 5 are simple to use, they usually provide only a crude approximation to fundamental periods of real frames. No mention is given in these equations to the semi-rigid connections, masses and stiffnesses of the structures despite that they are the most important factors which influence the fundamental mode period of vibrations. The Rayleigh equation is known to provide a satisfactory approximation of fundamental periods of structures, but it does not consider the effect of semi-rigid connection and it requires the use of linear elastic analysis to calculate the structure lateral displacement under the effect of a set of horizontal loads, H i , applied at each storey level. This requirement might be considered inconvenient since appropriate sub-assemblages might be used to analyze structures under the effect of either gravity or horizontal loading and complete analysis may never be required during the design of the structures. In this paper a simple equation to calculate the fundamental period of vibration is developed taking in the consideration the effect of semi-rigid connections, mass and the stiffness of the structure. This equation is simple to use and provides a satisfactory approximation of the fundamental mode period of vibrations for real structures. 2. Calculating the fundamental period T It is well known that the Rayleigh-Ritz method is the most popular technique for calculating an approximation to the fundamental periods, T 1, of multi-degree structures. In this technique an approximation
341
Dia Eddin Nassani / APCBEE Procedia 9 (2014) 339 – 346
to the lowest eigenvalue of the general eigenvalue problem >K @I Rayleigh quotient
I a T >K @I a I a T >M @I a
U (I ) a
O>M @I is given by the well-known
O1 U (I ) a O n
and
(7)
Where [K] is the global stiffness matrix of the structure and [M] is the global mass matrix of the structure. (Ia) an assumed approximation to the lowest eigenvector of the problem (I1). In equation 7, if the assumed vector (Ia) is chosen to be identical with the lowest eigenvector of the problem (I1), then the Rayleigh quotient coincides with the lowest eigenvalue ρ(I)1= λ1. The lowest eigenvector is usually unknown in practical cases and an assumed vector is utilized instead. Therefore the Rayleigh quotient will be an upper approximation to the lowest eigenvalue of the problem (see equation 7). The assumed vector (Ia) should be chosen as close as possible to the actual lowest eigenvector of the problem (I1), in order to obtain a good approximation to the lowest eigenvalue of the problem. In a practical analysis the assumed vector (Ia) can be taken as the lateral displacement vector (δ) resulting from subjecting a set of lateral load pattern (H) to the structure, i.e.
>K @G H
, I a
G
(8)
Substituting equation 8 into equation 7 and rearranging gives: n
¦G
G H G T >M @G T
U I a
i
Hi
i 1 n
¦m G i
(9)
2 i
i 1
Where δi is lateral displacement of storey i and Hi is horizontal load applied at storey i and mi is lumped mass of storey i and n is number of stories in the structure. Knowing that λ1=ω12, where ω1 is the lowest circular natural frequency of the structure, hence n
T1
2S
Z
¦mG i
2S
n
2 i
2S
i 1 n
¦H G i
i
i 1
¦wG i
2 i
i 1 n
(10)
g ¦ H iG i i 1
Where g is the gravity acceleration. Equation 10 is known by most international codes [4]. In order to facilitate the use of equation 10, and since the assumed vector (δ) is only an approximation to the real lowest eigenvector of the structure, the applied horizontal loads can be assumed to consist of the total storey gravity load applied horizontally at each storey level, i.e. Hi = wi Hence equation 10 becomes: n
T1
2S
¦w G i
2 i
i 1 n
g ¦ wi G i i 1
To further simplify equation 11, the following assumptions are adopted x The storey gravity loads are uniform in the structure, i.e. wi = cons.
(11)
342
Dia Eddin Nassani / APCBEE Procedia 9 (2014) 339 – 346
x Assuming that the storey displacement δi is linearly proportional with δn the displacement at upper storey hi (i.e. G i uG n ) ht where δn is the horizontal displacement at upper storey of the structure , hi is the storey height taken from ground level and ht is the total height of the structure. x Storey heights are uniform and equal to h (i.e. h=ht /n). Hence equation 11 becomes: wi u h 2 u T1
2S
Gn2 ht
2
g u wi u h u
1
2 2 ..... n 2
Gn
1 2 3 ....n
2
ht
(12)
Rearranging equation 12 mathematically gives; hu T1
Gn
2S
ht
n(n 1)(2n 1) 6 n g u ( n 1) 2 u
(13)
Assuming now that 2n 1 | 2n and rearranging equation 13 gives T1
2S
2G n 3g
(14)
It should noted here that the assumptions mentioned above should not reduce the ability of equation 14 to handle practical problems since a parametric study conducted by the author, not presented in this paper, has shown that these assumptions have only minimum influence on the calculated fundamental period. Moreover these assumptions are often satisfied in practical problems. 3. The proposed method In order to calculate an approximation to G n in equation 14, the model shown in figure 1 is adopted. This model consists of a column connected to upper and lower beams, which for the sake of simplicity, are assumed to be identical. Calculating the end moments of the model and taking anti-clock wise as positive gives: M 21 M 2B
M 12
EI c ª 6Gi º «6T » Lc ¬ Lc ¼
· º EI b ª § 1 «3 ¨ ¸ T» Lb ¬ © 1 3E ¹ ¼
(15) (16)
End moment equation of member supported by an elastic rotational spring of stiffness R from the first side and pinned from other side [5]. Where E young modulus, Hi is horizontal Load, Ic is the column's moment of inertia, Lc is column's height , Ib is beam's moment of inertia, Lb is beam's span, θ is rotation at upper and lower column ends, δi is storey drift and R is connection rotational stiffness.
343
Dia Eddin Nassani / APCBEE Procedia 9 (2014) 339 – 346
EI b LR
E
- Equilibrium at joint 2 requires that
¦
2
§ 6 EI c 3EI b § 6 EI c G i ·· 1 ¨ ¸¸ 0¨ ¨ ¸ ¨ L ¸T L E L Lc 1 3 © ¹¹ c b c ©
M 21 M 2 B
§ 6 EI c ¨ Lc ¨ ¨ 6 EI 3EI b c ¨ ¨ L Lb c ©
T
0
Gi
· ¸ Lc ¸ § ·¸ 1 ¨ ¨ 1 3E ¸ ¸¸ ¸ © ¹¹
(17)
- Considering the equilibrium of the column one obtains:
¦
column
0 M 12 M 21 H i Lc
M
(18)
0
Substituting equations 16 and 17 and accomplishing great amount of manipulation gives: 3
G
HLc 12 EI c Gi
where
Gi
§ · 1 ¨ ¨ 1 3E ¸ ¸ © ¹ 3EI b § · 6EI c 1 ¨ ¸ ¸ Lb ¨ Lc © 1 3E ¹ 3EI b Lb
(19)
Gi
Hi
2
Lc
Ib
B
Ic 1
Ib
A Lb
Fig. 1. Proposed model
In order to use equation 19 for a multistory frame, the following should considered: In sway frames each beam buckles in double curvature. Hence the rotational stiffness of the beam should be multiplied by 2. In real frames, each beam restrains two columns, therefore the beam rotational stiffness should be multiplied by 2. Finally each beam in a real multistory sway frame restrains two columns, i.e. upper and lower columns connected to that beam. Thus the beam rotational stiffness should approximately be divided by 2. As a result of above, the use of equation 19 for real multistory frame is permitted providing that the beam rotational stiffnesses is multiplied by 2, or
¦
Gi
sto rey
6EI b sto rey L b
¦
§ · 1 ¨ ¨ 1 3E ¸ ¸ © ¹ 6EI § · 1 c ¨ ¸ ¨ ¸ Lc © 1 3E ¹
6EI b Lb
¦
sto rey
Ib L sto rey b
¦
§ · 1 ¨ ¨ 1 3E ¸ ¸ © ¹ I § · 1 ¨ ¸ c ¨ ¸ 1 3 E L © ¹ c
Ib Lb
(20)
344
Dia Eddin Nassani / APCBEE Procedia 9 (2014) 339 – 346
It should be noted here that equation 20 underestimates the stiffness of both the first storey, where columns might be fully fixed at base, and topmost storey, where no upper columns are present. This underestimation will not affect the accuracy of equation 20 since the second assumption used to derive equation 14 slightly overestimate the frame stiffness, hence the two effects are likely to be more or less balanced. The factor Gi should be taken as the ratio of the sum of beam stiffnesses of the storey i to the sum of the total storey stiffnesses (beams plus columns). Corresponding to equation 20 the upper storey drift resulting from subjecting the structure to the storey gravity load applied horizontally at each storey level can be calculated as follows: Gn
n
¦G
n
i
i 1
WLc
¦ 12 EI i 1
3
§ n 1 i · ¸ ¨ ¸ ¨ n Gi c ¹ ©
W 12 En
n
¦ i 1
3
Lci I ci
§ n 1 i · ¸ ¨ ¸ ¨ Gi ¹ ©
(21)
Where Ici the sum of the column's moment of inertias of storey i Lci the height of storey i Gi the ratio of beam stiffnesses to column and beam stiffnesses of the storey i (equation 20) W is the total weight of the structure. Therefore T1
2S
W 18Egn
n
¦ i 1
3
Lci I ci
§ n 1 i · ¨ ¸ ¨ ¸ Gi © ¹
(22)
It should be noted that equation 22 is similar to the equation used to calculate the fundamental period of a where m and k are the mass and the stiffness of the single degree of freedom, i.e. T 2S m k
corresponding single degree of freedom structure. 4. Verification of the proposed method In order to verify the accuracy of the proposed method (equation 22) to calculate a good approximation to the fundamental natural periods of structures, the following examples are selected from many examples solved by the author and not included within. 4.1. Example 1 The frame shown in figure (2) which consists of ten multistory two bays steel frame with following properties: Beam section: UB 203x133x30, Column section: UC 254x254x89 Lc = 300 cm, Ac = 113 cm2, Ib = 2896 cm4, Ic = 14270 cm4, σy = 275 N/mm2, E = 210000 N/mm2 Typical extended end plate connection with rotational stiffness R= 10000 kN.m/rad. This value has been chosen from Scerbo [6]. Gravity load for stories (1-9) 20 N/mm and for upper storey 10 N/mm Calculating storey properties for the first storey of the frame produces: E1
EIb Lb R
210000 u102 u 2896 400 u10000 u105
0.152, E2
EI b Lb R
210000 u102 u 2896 600 u10000 u105
2896 § 1 1 · 2896 § · u¨ u¨ ¸ ¸ 600 © 1 3 u 0.102 ¹ 400 © 1 3 u 0.152 ¹ 2896 § 1 1 14270 · 2896 § · 2u u¨ u¨ ¸ ¸ 4u 600 © 1 3 u 0.102 ¹ 400 © 1 3 u 0.152 ¹ 300
0.102
2u
G1
0.061
345
Dia Eddin Nassani / APCBEE Procedia 9 (2014) 339 – 346 300 300 300 300 300 300 300 300 300 300 600
400
600
Fig. 2. The frame of example 1
Repeating same procedure for all stories of the frame provides table 1 Table 1. Calculation of example 1 Storey No. 1 2 3 4 5 6 7 8 9 10
I ci
¦I
3
3
c
storey
x104 mm4 57080 57080 57080 57080 57080 57080 57080 57080 57080 57080
Gi 0.061 0.061 0.061 0.061 0.061 0.061 0.061 0.061 0.061 0.061
Lci I ci
¦
(n 1 i)
1/mm 473 425 378 331 283 236 189 142 94 47
Lci I ci
¦
2598
(n 1 i) Gi 1/mm 7754 6967 6196 5426 4639 3868 3098 2327 1540 770 42585
Employing the well known computer program STAAD.Pro [7] to calculate the fundamental period produces row 2 in table 2. Calculating the fundamental period of the frame utilizing equation (22) provides the value of raw 3 in table2. Employing the conventional method (equation 1) produces the value tabulated in raw 4 table 2. Table 2. Fundamental Period (sec.) for Example 1 Steel structure with semi-rigid connection STAAD.Pro Proposed method Conventional method equation 1
4.1 3.7 1.1
It can be concluded from table 2 that the prediction of the proposed method is more accurate than the conventional method. 4.2. Example 2 Frame consists of three multistory, one bays steel frame with following properties: Beam section: UB 203x102x23, Column section: UC 152x152x23, Lc = 300 cm, Ac = 29.2 cm2, Ib = 2105 cm4, Ic = 1250 cm4, σy = 275 N/mm2, E = 210000 N/mm2, flush end plate connection R= 3800 kN.m/rad [6]. Gravity load for stories (1-2) 30 N/mm and for upper storey 20 N/mm
346
Dia Eddin Nassani / APCBEE Procedia 9 (2014) 339 – 346
Table 3. Fundamental Period (sec.) for Example 2 Steel structure with semi-rigid connection STAAD.Pro Proposed method Conventional method equation 1
1.4 1.7 0.42
It can be concluded from table 3 that the prediction of the proposed method is more accurate than the conventional method. 5. Conclusion The design of real structures to withstand earthquake actions require the calculation of the fundamental mode period of vibration of these structures. Many equations have been suggested by the international codes to calculate the fundamental periods and most of these equations provide only a crude approximation. The proposed formula give a satisfactory approximation of fundamental period of the structure taking into account the mass and the stiffness of the structure including the effect of semi-rigid connections. The accuracy of the proposed formula has been proven by solving several examples to calculate the fundamental mode periods. Comparing the obtained results with those obtained from STAAD.Pro has shown that the proposed formula provides a more accurate estimation of the fundamental period of the structures than the conventional method. References [1] ASCE 7-10. Minimum Design Loads for Buildings and Other Structures. American Society of Civil Engineers. USA; 2010. [2] IBC. International Building Code. International Code Council, INC: USA; 2012. [3] SEAOC. Recommended lateral force requirements. Structural Engineers Association of California. San Francisco; 1996. [4] CEN, Eurocode 8. Design of structures for earthquake resistance. Part 1: general rules, seismic actions and rules for buildings, European Standard EN 1998- 1, 2004, Comité Européen de Normalisation, Brussels: Belgium; 2004. [5] NASSANI DE. Static and Dynamic Behavior of Frames with Semi-Rigid Connection. Ph.D. Dissertation . Aleppo University: Syria; 2011. [6] SCERBO JS. Analysis of steel frames with deformable beam to column connections. M.Sc. Thesis. The University of Manitoba: Canada; 1996. [7] STAAD PRO. Technical reference manual. Research Engineer International; 2012.