Scribd
Upload
182 views

Example-Modal Analysis Mathcad

1) A modal analysis was performed on a 3 degree of freedom system to calculate its natural frequencies, mode shapes, and modal forces. 2) Three natural frequencies were found to be 4.592 Hz, 9.818 Hz, and 14.578 Hz. 3) Corresponding mode shapes and modal forces were determined for each natural frequency/mode of vibration. 4) Total base shear and overturning moments at different floor levels were calculated by combining the modal forces.

Uploaded by

Mad Max
Copyright
© Attribution Non-Commercial (BY-NC)
Available Formats
Download as PDF, TXT or read online on Scribd
182 views

Example-Modal Analysis Mathcad

1) A modal analysis was performed on a 3 degree of freedom system to calculate its natural frequencies, mode shapes, and modal forces. 2) Three natural frequencies were found to be 4.592 Hz, 9.818 Hz, and 14.578 Hz. 3) Corresponding mode shapes and modal forces were determined for each natural frequency/mode of vibration. 4) Total base shear and overturning moments at different floor levels were calculated by combining the modal forces.

Uploaded by

Mad Max
Copyright
© Attribution Non-Commercial (BY-NC)
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 7

ORIGIN := 1

k := 1000 lb

60 60 0
k
K := 60 180 120

in
0 120 300

1
2

1 0 0
2
k sec
M := 0 1.5 0

in
0 0 2
K1 := M

60 60 0
-2
K1 = 40 120 80 sec

0 60 150
evp := eigenvals( K1)

21.088
-2
evp = 96.396 sec

212.516

Freq := sort( evp)

i := 1 .. 3
i :=

Freq

T :=

1 = 4.592 sec

-1

2 = 9.818 sec

-1

21.088
-2
Freq = 96.396 sec

212.516

2
i
T = 1.368 sec
1

3 = 14.578 sec

T = 0.64 sec
2

-1

C:\AAA_Files\Teaching\CE7119
\AAA_Mathcads\Example-modal
analysis.mcd/Example-modal
analysis.mcd

T = 0.431 sec
3

CIVL7119/8119

1/7

Find mode shapes:


v := eigenvecs( K1)

0.813 0.739 0.273


v = 0.527 0.449 0.694

0.246 0.502 0.666


Mode shape 1 corresponds with eigenvalue 21.088. The first column of the v matrix
corresponds with the first evp which was 21.088. The second column corresponds to
the frequency 96.396 (second row of evp). The third column corresponds to the frequency
212.516.

V1 :=

1
v
v

1, 1

V2 :=

2
v
v

1, 2

V3 :=

3
v
v

1, 3

1
V1 = 0.649

0.302
1
V2 = 0.607

0.679
1
V3 = 2.542

2.44

2.531 105
2.531 105

2
-1

K V1 = 2.462 105 lb ft
(1) M V1 = 2.462 105 lb ft-1

5
1.528 10
1.528 10
1.157 106
1.157 106

2
-1
-1
(2) M V2 = 1.053 106 lb ft K V2 =
6
1.053 10 lb ft

6
1.571 10
1.571 10

C:\AAA_Files\Teaching\CE7119
\AAA_Mathcads\Example-modal
analysis.mcd/Example-modal
analysis.mcd

CIVL7119/8119

2/7

Check orthogonality condition:


T

Test1 := V1 M V2

Test1 = 6.46 10

Test2 := V1 M V3

12

) lb sec2 ft-1

Test2 = 1.331 10

11

) lb sec2 ft-1

Test3 := V2 M V3

Test3 = 3.205 10

12

) lb sec2 ft-1

Sv1 := 1.74

1
one := 1

Sv2 := 1.41
Sv3 := 1.20

ft
sec
ft
sec
ft
sec

L1 := V1 M one

L2 := V2 M one

L3 := V3 M one

M 1 := V1 M V1

M 2 := V2 M V2

M 3 := V3 M V3

M 1 = ( 1.813 )
M 2 = ( 2.474 )

k sec
in
k sec

M 3 = ( 22.596 )

L1 = ( 2.577 )

k sec
in

in
k sec

L2 = ( 1.268 )
2

in

C:\AAA_Files\Teaching\CE7119
\AAA_Mathcads\Example-modal
analysis.mcd/Example-modal
analysis.mcd

L3 = ( 2.066 )

k sec

in

k sec

in

CIVL7119/8119

3/7

Y1 := L1 M 1

Y2 := L2 M 2

Y3 := L3 M 3

1 Sv1

Y1 = ( 0.538 ) ft

1 Sv2

Y2 = ( 0.074 ) ft

1 Sv3

Y3 = 7.528 10

) ft

u1 := V1 Y1
u2 := V2 Y2
u3 := V3 Y3

0.074
u2 = 0.045 ft

0.05

0.538
u1 = 0.349 ft

0.163

7.528 10 3

u3 =
ft
0.019
0.018

SRSS Methods

umax1 :=

(u11) + (u21) + (u31)

umax2 :=

(u12) + (u22) + (u32)

umax3 :=

(u13) + (u23) + (u33)

umax1 = 0.543 ft
umax2 = 0.353 ft
umax3 = 0.171 ft

C:\AAA_Files\Teaching\CE7119
\AAA_Mathcads\Example-modal
analysis.mcd/Example-modal
analysis.mcd

CIVL7119/8119

4/7

CQC Procedure
i := 1 .. 3
:= 0.05
j := 1 .. 3

i, j

:=

i, j

=
1

2.138
3.175
0.468
1
1.485

3
2

) (ri, j)

8 1 + r
i , j :=

0.315

0.673

i, j

1 r
+ 4 r ( 1 + r )
i, j
i, j
( i , j)
2

u11 u21 u31


U := u12 u22 u32

u13 u23 u33

1
2

0.538 0.074 7.528 10 3

U = 0.349 0.045
ft
0.019
0.163 0.05
0.018

p := 1 .. 3

Ucqc :=
p

U
p , m m, n p , n

m= 1 n= 1

0.542
Ucqc = 0.353 ft

0.172

Compare this with SRSS results

Sa1 := 1 Sv1
Sa2 := 2 Sv2
Sa3 := 3 Sv3

C:\AAA_Files\Teaching\CE7119
\AAA_Mathcads\Example-modal
analysis.mcd/Example-modal
analysis.mcd

CIVL7119/8119

5/7

Modal Forces
f1 := M V1 L1

f2 := M V2 L2

f3 := M V3 L3

136.254
f1 = 132.549 k

82.257

Sa1
M1

85.134
f2 = 77.464 k

115.609

Sa2
M2

Sa3

19.197
f3 = 73.197 k

93.668

M3

Base Shear for Each Mode


3

V01 :=

f1

V01 = 351.06 k

j= 1

V02 :=

f2

V02 = 107.938 k

j= 1

V03 :=

f3

V03 = 39.668 k

j= 1

Total Base Shear

V0 :=

V01 + V02 + V03

V0 = 369.414 k

C:\AAA_Files\Teaching\CE7119
\AAA_Mathcads\Example-modal
analysis.mcd/Example-modal
analysis.mcd

CIVL7119/8119

6/7

Total Lateral Force


Froof :=

(f11) + (f21) + (f31)

F2nd :=

(f12) + (f22) + (f32)

F1st :=

(f13) + (f23) + (f33)


2

Froof = 161.807 k

F2nd = 170.081 k

F1st = 1.7 10 lb

Overturning Moments at Floor Levels


H1 := 12 ft
H2 := 12 ft
H3 := 12 ft

OMGround := ( H1 + H2 + H3) Froof + ( H2 + H3) F2nd + H3 F1st

OM1st := ( H1 + H2) Froof + H2 F2nd

OM2nd := H1 F1st
4

OMGround = 1.195 10 k ft
3

OM1st = 5.924 10 k ft
3

OM2nd = 2.04 10 k ft

C:\AAA_Files\Teaching\CE7119
\AAA_Mathcads\Example-modal
analysis.mcd/Example-modal
analysis.mcd

CIVL7119/8119

7/7

You might also like