Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Chimney Seismic Analysis

Download as xlsx, pdf, or txt
Download as xlsx, pdf, or txt
You are on page 1of 9

SEISMIC ANALYSIS OF A TALL UNIFORM RCC CHIMNEY: (IS : 1893 - 1984)

DATE: 5/24/2022
Height of chimney, H = 367.5 m
Base diameter, D = 10.5 m
Thickness, t = 0.5 m
H/D ratio = 367.5/10.5
= 35.000

0.5
m

367.5 m

CHIMNEY READY TO RECEIVE THE FLUE DUCT


FROM BOILER HOUSE

Inner diameter, d = 10.5-(2*0.5)


= 9.5 m
Ac = PI()*((10.5^2)-(9.5^2))/4 Ac = Area of concrete cross section
= 15.708 m2
I= PI()*((10.5^4)-(9.5^4))/64 I = Moment of inertia of concrete cross section
= 196.84041 m4
Concrete grade = M 40
E= 5000*SQRT(40) E = Young's modulus of concrete
= 31622.777 N/mm2
Unit weight of concrete, γconc = 25 kN/m3
Seismic zone = V
Foundation type = Raft
Damping factor, β = 5 % (for concrete)
SEISMIC ANALYSIS OF A TALL UNIFORM RCC CHIMNEY: (IS : 1893 - 1984)
DATE: 5/24/2022
As per Table 4,

Importance factor, I = 2
The chimney is divided into 4 parts.

F1 = 1838.901
kN-s2/m 1
91.875 m
F2 = 3677.801 2
kN-s2/m
91.875 m
F3 = 3677.801 3 We have considered equal intervals.
kN-s2/m 91.875 m

F4 = 3677.801 4
kN-s2/m 91.875 m

LUMPED MASS MODEL

Total weight of the chimney = 15.708*367.5*25


= 144316.91 kN
SEISMIC ANALYSIS OF A TALL UNIFORM RCC CHIMNEY: (IS : 1893 - 1984)
DATE: 5/24/2022
Total mass of the chimney = 144316.91/9.81
= 14711.204 kN-s2/m

F 2

a'
4 x

a' = height at which the load is applied


x = height at which deflection needs to be calculated

EIy = (F/6)((a'-x)3+(3a'2x)-(a'3)) if x<a'


= F((a'2x/2)-(a'3/6)) if x>a'

Deflection coefficient = EIy/FH 3

Deflection at
1 (x = 2 (x = 3 (x = 4 (x =
Load at
367.5 275.625 183.75 91.875
m) m) m) m)
1 (a' = 367.5 m) x≥a' x<a' x<a' x<a'
2 (a' = 275.625 m) x≥a' x≥a' x<a' x<a'
3 (a' = 183.75 m) x≥a' x≥a' x≥a' x<a'
4 (a' = 91.875 m) x≥a' x≥a' x≥a' x≥a'

Deflection coefficients:
Deflection at
1 (x = 2 (x = 3 (x = 4 (x = Load value
Load at
91.875 (kN-s /m)
2
367.5 275.625 183.75
m) m) m) m)
1 (a' = 367.5 m) 0.333 0.211 0.104 0.029 1838.901
2 (a' = 275.625 m) 0.211 0.141 0.073 0.021 3677.801
3 (a' = 183.75 m) 0.104 0.073 0.042 0.013 3677.801
4 (a' = 91.875 m) 0.029 0.021 0.013 0.005 3677.801
SEISMIC ANALYSIS OF A TALL UNIFORM RCC CHIMNEY: (IS : 1893 - 1984)
DATE: 5/24/2022

Modal analysis for first mode:

Sample calculation of deflection for first mode shape:

Deflection at 2, δ2 = (H3/EcI)*(Σ(Fr*deflection coefficient))*10 3


= (1838.901*0.211)+(2574.461*0.141)+(1103.340*0.073)+(367.780*0.0(H3/EcI)*103
= 838.041 (H3/EcI)*103

Mr Fr = Mrφra δ (mm) φ r'


φra
Load at (initializati Cycle
(kN-s2/m) (assumed) (kN-s2/m) (H3/EcI)*103 on)
1 1838.901 1.000 1838.901 1281.484 1.000
2 3677.801 0.700 2574.461 838.041 0.654
1st cycle
3 3677.801 0.300 1103.340 430.035 0.336
4 3677.801 0.100 367.780 122.593 0.096
1 1838.901 1.000 1838.901 1258.940 1.000
2 3677.801 0.654 2405.139 823.438 0.654
2nd cycle
3 3677.801 0.336 1234.180 422.932 0.336
4 3677.801 0.096 351.837 120.686 0.096

EcI/H3 = (31622.78*1000)*196.8404/367.5^3
= 125.413 m
= 125412.91 mm

M rφ r' Mr(φr')2 F r φ r'


Load at
(kN-s2/m) (kN-s2/m) (kN-s2/m)
1 1838.901 1838.901 1838.901
2 2405.548 1573.403 1573.135
3 1235.532 415.069 414.614
4 352.567 33.798 33.728
Σ 5832.548 3861.170 3860.378

Frequency, ω1 = ((3860.378/3861.170)*(125412.9/(1258.940*1000)))^0.5
= 0.316 rad/sec
Participation factor, P1 = 5832.548/3861.170
= 1.511
Time period, T1 = 2π/ω1 = 2*PI()/0.316
= 19.909 sec
SEISMIC ANALYSIS OF A TALL UNIFORM RCC CHIMNEY: (IS : 1893 - 1984)
DATE: 5/24/2022
As per Figure 2,

Sa1/g = 0.04
Sa1 = 0.04*9.81 Sa1 = Seismic spectral acceleration for first mode
= 0.392 m/s2

Shear and bending moment calculation (considering only first mode effect):
SEISMIC ANALYSIS OF A TALL UNIFORM RCC CHIMNEY: (IS : 1893 - 1984)
DATE: 5/24/2022
As per Table 3,

β= 1 (for Raft foundation)


As per Table 2,
SEISMIC ANALYSIS OF A TALL UNIFORM RCC CHIMNEY: (IS : 1893 - 1984)
DATE: 5/24/2022

Zone factor, F0 = 0.4 (for Zone V )


Sa1 (mod) = Sa1P1βIF0
= 0.392*1.511*1*2*0.4
= 0.474 m/s2
Fr1 = Mrφr1Sa1 (mod) = inertia force at rth level of the chimney for mode 1

F11

F21

F31

F41

Inertia forces at various levels of chimney

F11 = 1838.901*1.000*0.47 = 872.000 kN


F21 = 3677.801*0.654*0.47 = 1140.703 kN
F31 = 3677.801*0.336*0.47 = 585.885 kN
F41 = 3677.801*0.096*0.47 = 167.186 kN

Base shear = 872.000+1140.703+585.885+167.186


= 2765.774 kN
Bending moment at base = (F11 X H) + (F21 X 0.75H) + (F31 X 0.5H) + (F41 X 0.25H)
(872.000*(91.875+91.875+91.875+91.875))
= +(1140.703*(91.875+91.875+91.875))+(585.885*(91.875+91.875))
+(167.186*91.875)
= 757882.81 kN-m
SEISMIC ANALYSIS OF A TALL UNIFORM RCC CHIMNEY: (IS : 1893 - 1984)
DATE: 5/24/2022

2nd mode analysis:

ΣMr(φr')2 = 3861.170 kN-s2/m ψi = ΣMrφriφra


ΣMr(φr')2

Mr φri (1st φra Mrφriφra ψiφri φraφ = φra- Fr = Mrφraφ δ (mm)


Load at
(kN-s /m)
2 mode) (assumed) (kN-s /m)
2 ψiφri (kN-s /m) (H /EcI)*103
2 3

1 1838.901 1.000 1.000 1838.901 0.351 0.649 1193.138 25.201


2 3677.801 0.654 0.100 240.555 0.230 -0.130 -476.971 -6.686
3 3677.801 0.336 -0.500 -617.766 0.118 -0.618 -2272.780 -21.172
4 3677.801 0.096 -0.300 -105.770 0.034 -0.334 -1227.151 -11.743
Σ= 1355.919 ψ1 = 0.351
1 1838.901 1.000 1.000 1838.901 0.000 1.000 1839.655 36.719
2 3677.801 0.654 -0.265 -638.195 0.000 -0.265 -974.738 -9.983
3 3677.801 0.336 -0.840 -1037.996 0.000 -0.840 -3089.290 -30.478
4 3677.801 0.096 -0.466 -164.294 0.000 -0.466 -1713.685 -16.759
Σ= -1.585 ψ2 = 0.000
φ r'
Load at (initializati Cycle
on)
1 1.000
2 -0.265
1st cycle
3 -0.840
4 -0.466

1 1.000
2 -0.272
2nd cycle
3 -0.830
4 -0.456

M rφ r' Mr(φr')2 F r φ r'


Load at
(kN-s2/m) (kN-s2/m) (kN-s2/m)
1 1838.901 1838.901 1839.655
2 -999.896 271.845 265.005
3 -3052.73 2533.892 2564.239
4 -1678.63 766.164 782.164
Σ -3892.35 5410.801 5451.064

Frequency, ω2 = ((5451.064/5410.801)*(125412.9/(36.719*1000)))^0.5
= 1.855 rad/sec
Participation factor, P2 = -3892.35/5410.801
= -0.719
Time period, T2 = 2π/ω2 = 2*PI()/1.855
= 3.387 sec
SEISMIC ANALYSIS OF A TALL UNIFORM RCC CHIMNEY: (IS : 1893 - 1984)
DATE: 5/24/2022

Sa2 (mod) = Sa1P2βIF0


= 0.392*-0.719*1*2*0.4
= -0.226 m/s2
Fr2 = Mrφr2Sa2 (mod) = inertia force at rth level of the chimney for mode 2

F11

F21

F31

F41

Inertia forces at various levels of chimney

F11 = 1838.901*1.000*-0.226 = -415.268 kN


F21 = 3677.801*-0.272*-0.226 = 225.800 kN
F31 = 3677.801*-0.830*-0.226 = 689.379 kN
F41 = 3677.801*-0.456*-0.226 = 379.075 kN

Base shear = -415.268+225.800+689.379+379.075


= 878.986 kN
Bending moment at base = (F11 X H) + (F21 X 0.75H) + (F31 X 0.5H) + (F41 X 0.25H)
=
(-415.268*(91.875+91.875+91.875+91.875))
+(225.800*(91.875+91.875+91.875))+(689.379*(91.875+91.875))
+(379.075*91.875)
= 71126.234 kN-m

You might also like