SEISMIC DYNAMIC COMPUTATION

I. Data :
Modulus of Elasticity of Steel Es := 200000  MPa

Modulus of Elsaticity of Concrete Ec := 21000  MPa

m
gravitational constant g = 9.807
2
s

kN
Concrete Unit weight γc := 23.5 
3
m

Building width Bw := 10m

Building length BL := 14m

Slab thickness tf := 0.15m

Sizes of Structural Members

Level Column height Column sizes

C1h := 3m Lc := C1h C1b := 500  mm C1a := 350  mm

Level 1 1

Level 2 C2h := 3m Lc := C2h C2b := 400  mm C2a := 300  mm

2

Level 3 C3h := 3m Lc := C3h C3b := 300  mm C3a := 300  mm

3

Beam sizes

Roof deck b4 := 250mm h4 := 350mm

Floor 3 b3 := 250mm h3 := 450mm

Floor 2 b2 := 300mm h2 := 500mm

Slab Weight

2
Floor Area: Af := Bw  BL Af = 140 m
 Floor weight Ws := Af  tf  γc Ws = 493.5  kN

Beam Weight
Total beam length TBL := 4  Bw + 3  BL TBL = 82 m

Roof Deck WBR := b4  ( h4)  TBL  γc WBR = 168613 N

Floor 5 WB3 := b3  ( h3)  TBL  γc WB3 = 216788 N

Floor 4 WB2 := b2  ( h2)  TBL  γc WB2 = 289050 N

Note : for T-beam - slab

system,
deduct the thickness of the
slab WBRe := b4  ( h4 - tf )  TBL γc
from the depth of the beam.
WBRe = 96350 N

Column Weight
Number of columns Nc := 12

 C3h 
Roof Deck WcR := Nc   C3b  C3a  γc WcR = 38.07  kN
 2 
 C3h C2h 
3rd floor Wc3 := Nc    C3b  C3a + C2b  C2a  γc Wc3 = 88.83  kN
 2 2 
 C2h C1h 
2nd floor Wc2 := Nc    C2b  C2a + C1b  C1a  γc Wc2 = 124.79  kN
 2 2 

Wall Weight

Total wall length WL := 3  BL + 4  Bw WL = 82 m

Wall load γw := 2.314  kPa

Roof deck parapet wall height Rh := 1.5m

Roof level WRh := Rh WwR := WRh  WL  γw WwR = 284.622  kN

Floor 3 W3h := C3h Ww3 := W3h  WL  γw Ww3 = 569.24  kN

Floor 2 W2h := C2h Ww2 := W2h  WL  γw Ww2 = 569  kN

 II. Total Weight per Floor Level

Roof Level WR := Ws + WBR + WcR + WwR WR = 985  kN

Floor 3 W3 := Ws + WB3 + Wc3 + Ww3 W3 = 1368  kN

Floor 2 W2 := Ws + WB2 + Wc2 + Ww2 W2 = 1477  kN

Type this command in order to set the first element of any matrix or vector with index as 1.

ORIGIN  1

Diagonal Mass Matrix

W2  150569 
  1  
Ma := W3   Ma = 139534  kg
  g  
WR  100422 

Ma1 0 0 
  150569 0 0 
   
M :=  0 Ma2 0  M= 0 139534 0  kg
   
 0 0 Ma   0 0 100422 
 3

III. Calculation of Column Stiffeness

Moment of Inertia of each column per floor

1   3
I :=  C3a   C3b  8 4
Floor Level 3 3 12 I = 6.75  10  mm
3
1   3
Floor Level 2 I :=  C2a   C2b  9 4
2 12 I = 1.6  10  mm
2
1   3
Floor Level 1 I :=  C1a   C1b  9 4
1 12 I = 3.646  10  mm
1

Individual Column Stiffness per floor Total column stiffness per

 floor

12  Ec  I
3 kN
third floor kc := k := Nc  kc k = 75600 
3 3 3 3 3 m
( C3h)

12  Ec  I
2 kN
kc := k := Nc  kc k = 179200 
second floor 2 3 2 2 2 m
( C2h)

12  Ec  I
1 kN
kc := k := Nc  kc k = 408333 
ground floor 1 3 1 1 1 m
( C1h)

The Stiffness in vector form

k1   5
  4.083  10 
    kN
k := k2  k = 1.792  105  
    m
k3  4
 7.56  10 
   

Assembly of stiffness matrix resisting lateral loads

 k1 -k 0 
 2 
 
K := -k2 k2 + k 3 -k3 
 
0 -k
3
k 
3 


 4.083  105 -1.792  105 

 0 
  kN
K = -1.792  105 2.548  105 -7.56  104  
  m
 4 4 
 0 -7.56  10 7.56  10 

IV. Eigenvalue Problem -

Calculating the dynamic properties of the building:

1 ) Circular frequency, ω
  3631 
 K   rad 2
Ω2 := eigenvals   Ω2 = 1362.5  
 M   s2
 297.4 

thus the circular frequencies are

60.257 
  rad
ω := Ω2 ω = 36.912  
  s
17.246 

sort the frequencies from lowest to highest

17.246 
  rad
ω := sort ( ω) ω = 36.912  
  s
60.257 

Consequently, the periods of virbation are

0.364 
2 π  
T := T =  0.17   sec
ω  
0.104 

V. Calculate the Mode Shapes

Calculate the mode shape for each frequency, using the eigenvec command for each eigenvalue

1) Mode shape 1 is the eigenvector, ϕ1, of the eigenvalue ω12

0.324 
K 2  
ϕ1 := eigenvec  ,  ω   ϕ1 = 0.609 
1
M     
0.724 

2) Mode shape 2 is the eigenvector, ϕ2, of the eigenvalue ω22

  -0.58 
K 2  
ϕ2 := eigenvec  , ω
 ϕ2 = -0.609 
2
M    
 0.541 

3) Mode shape 3 is the eigenvector, ϕ3, of the eigenvalue ω32

-0.808 
K 2  
ϕ3 := eigenvec  ,  ω   ϕ3 =  0.578 
3
M     
-0.109 

ϕ := augment ( ϕ1 , ϕ2 , ϕ3)

0.324 -0.58 -0.808 

 
ϕ = 0.609 -0.609 0.578 
 
0.724 0.541 -0.109 

Normalizing the mode shapes to take advanage of the orthogonal property

rows of matrix n := 3
columns of matrix i := 1 .. n j := i

ϕ 0.00093 -0.0016 -0.00211 

i,j  
Φ :=
i,j Φ = 0.00176 -0.00168 0.00151 
n  M 
2 k , k  
0.00209 0.00149 -0.00028 

k = 1
 ϕk , j  kg 
 

Similarly, the matrix of the eignevvectors can be disassembled into vectors of mode shapes

1 2 3

Φ1 := Φ Φ2 := Φ Φ3 := Φ
 0.00093   -0.0016  -0.002114 
     
Φ1 = 0.00176  Φ2 = -0.00168  Φ3 =  0.001513 
     
0.00209   0.00149  -0.000285 

VI. Calculate Modal Particiapation Factor

1 
influence coefficient  
R := 1 
 
1 

1) Modal participation factor for mode shape 1

T
Φ1  M  R
Γ := Γ = 595.48
1 T 1
Φ1  M  Φ1

2) Modal participation factor for mode shape 2

T
Φ2  M  R
Γ := Γ = -324.8
2 T 2
Φ2  M  Φ2
3) Modal participation factor for mode shape 3
T
Φ3  M  R
Γ := Γ = -135.794
3 T 3
Φ3  M  Φ3

The modal particiaption factor in vector form

 595.48 
 
Γ = -324.80 
Γ := Γ
 
-135.79 

VII. Calculate Effective Modal Mass and Mass Participation

a) Effective modal mass

 354596 
2  
Me := Γ
Me = 105495 
 
 18440 

total effective mass Met :=

 Me Met = 478531

b) Percent participation of each effective mass

74.10 
Me  
percent :=  100 percent = 22.05 
Met  
 3.85 

Total mass participation

i  74.1 
 
TMP :=
i  percent
k
TMP = 96.15 
 
k = 1  100 

VIII. Applying the Ground Motion determine the corresponding

Spectral Acceleration for each Period

From NSCP Table Na := 1.0 Nv := 1.0

Ca := 0.44  Na Cv := 0.64  Nv

Ca = 0.44 Cv = 0.64

Cv  s
control period Ts := Ts = 0.582 s
2.5  Ca

Corresponding spectral acceleration for each periods

 0.364 
 
T =  0.17  s
 
0.104 

Cv  s m
T = 0.364 s Sa := g Sa = 17.227
1 1 T 1 2
1 s
m
T = 0.17 s Sa := 1.1  g Sa = 10.787
2 2 2 2
s
m
T = 0.104 s Sa := 1.1  g Sa = 10.787
3 3 3 2
s

IX. Calculate the Lateral Forces for each mode and the Total Base Shear

Importance Factor I := 1.0

R := 8.5
Response Modification Factor
Mode shape 1:

170 
I  
F1 :=  M   Φ1  Γ  Sa  F1 = 296   kN
R  1 1
 
253 

Mode shape 2:

 99 
I  
F2 :=  M   Φ2  Γ  Sa  F2 =  97   kN
R  2 2
 
-62 

Mode shape 3:

 55 
I  
F3 :=  M   Φ3  Γ  Sa  F3 = -36   kN
R  3 3
 
5 
 Calculate the maximum response in terms of lateral forces
Using the modal combination SRSS (square-root of the sum of the squares)

2 2 2
Fx := F1 + F2 + F3

204 
 
Maximum Lateral forces Fx = 313   kN
 
261 

Dynamic Base Shear

Vd :=
 Fx Vd = 778  kN

X. Calculate the Lateral Displacements for each Mode and the

Total Maximum Displacements

Mode shape 1

Sa 3.79 
I 1  
X1 :=  Φ1  Γ  X1 = 7.13   mm
R 1 2
 ω1  
  8.48 

Mode shape 2

Sa  0.48 
I 2  
X2 :=  Φ2  Γ  X2 =  0.51   mm
R 2 2
 ω2  
  -0.45 

Mode shape 3

Sa  0.1 
I 3  
X3 :=  Φ3  Γ  X3 = -0.07   mm
R 3 2
 ω3  
   0.01 
 Maximum displacement response using the SRSS (square-root-of the sum-of-the-squares)

2 2 2
X := X1 + X2 + X3

3.822 
 
Maximum Floor Displacements X = 7.144   mm
 
 8.49 

Calculate the Interstory Drift

Interstory Drifts
Dx := X
1 1

Dx := X - X
2 2 1

Dx := X - X
3 3 2

3.82 
 
Dx = 3.32   mm
 
1.35 

Calculating the Inelastic Drift using the NSCP

22.742 
 
Δi := 0.7  R  Dx Δi = 19.766   mm
 
 8.005 

NSCP drift limits

75 
 
Dr := 0.025  Lc Dr = 75   mm
 
75 

XI. Equivalent Lateral Force Procedure

W := Ma  g
i i
 1477 
 
Weight per floor W = 1368   kN
 
 985 

Total weight
Wt :=
W Wt = 3830  kN

Height per floor from ground level

4 
 
Hn :=  7 
 
10 

hn := Hn hn = 10
3

Zone = 4 Z := 0.4

Importance Factor I := 1.0

Response Modification Factor R := 8.5

Distance from nearest Fault Line D := 50

(km)
Soil Type , Sd
Ca = 0.44 Cv = 0.64

Fundamental Period of Vibration

hn
for RC structures Ct := 0.0731 hn :=
m
3
4 1
T := Ct  hn T = 0.411
0.75
m

Base shear

Cv  I  Wt
Vsb := 1
R T Vsb = 701 m  kN
 Upper limit

2.5  Ca I
Vupper :=  Wt Vupper = 496  kN
R

Lower limits

Vlow1 := 0.11  Ca  I  Wt Vlow1 = 185  kN

0.8  Z  Nv  I
Vlow2 :=  Wt Vlow2 = 144  kN
R

Governing Static Base Shear

Vs := 862kN

Distribution of base shear to floor level forces

Vs  W  Hn
i i
Fxs :=
i n


k = 1
 W k  Hnk
 

201  4 
   
Fxs = 326   kN Hn =  7 
   
335  10 

Comparison of Static versus Dynamic Base Shear

The code specifiies that the dynamic base shear must not be less than 90% of the static
base shear. Therefore scale up the dynamic base shear before it is distributed to
equivalent lateral forces.

Vd
Rv := Rv = 0.902 this is less than 90%
Vs
 0.9  Vs
SF := SF = 0.997
Vd

The adjusted Dynamic Base shear to be distributed is

Vda := SF  Vd Vda = 776  kN

 203.5 
 
Fxa := SF  Fx Fxa = 312.34   kN
 
259.96 