Contoh Soal MDOF Rev
Contoh Soal MDOF Rev
Contoh Soal MDOF Rev
Jawab:
𝑁 1
𝐸𝐼 = (2 × 104 2
) ( × 150 × 2003 𝑚𝑚4 ) = 2 × 1012 𝑁𝑚𝑚2
𝑚𝑚 12
Asumsi: Untuk tiap lantai, massa dijadikan lumped mass dengan rigid beam dan deformasi
aksial diabaikan.
12𝐸𝐼
𝑘= (𝑓𝑖𝑥𝑒𝑑 − 𝑟𝑖𝑔𝑖𝑑 𝑏𝑒𝑎𝑚)
𝐿3
Maka,
• Untuk massa 1:
𝑚1 𝑢̈ 1 + 3𝑘1 𝑢1 − 3𝑘2 (𝑢2 − 𝑢1 ) = 0
𝑚1 𝑢̈ 1 + (3𝑘1 + 3𝑘2 )𝑢1 + (−3𝑘2 )𝑢2 = 0
• Untuk massa 2:
𝑚2 𝑢̈ 2 + 3𝑘2 (𝑢2 − 𝑢1 ) + 𝑘3 𝑢2 + (−2𝑘2 )(𝑢3 − 𝑢2 ) = 0
𝑚2 𝑢̈ 2 + (−3𝑘2 )𝑢1 + (5𝑘2 + 𝑘3 )𝑢2 + (−2𝑘2 )𝑢3 = 0
• Untuk massa 3:
𝑚3 𝑢̈ 3 + 2𝑘2 (𝑢3 − 𝑢2 ) = 0
𝑚3 𝑢̈ 3 + (−2𝑘2 )𝑢2 + 2𝑘2 𝑢3 = 0
𝑀𝑈̈(𝑡) + 𝐾𝑈(𝑡) = 0
Selanjutnya, cari matriks massa dan kekakuan struktur dari persamaan di atas
𝑚1 0 0 200 0 0
→𝑀=[ 0 𝑚2 0 ]=[ 0 100 0 ] 𝑘𝑔
0 0 𝑚3 0 0 200
2804300 −1679300 0
→ 𝐾 = [−1679300 3358601 −1119534] 𝑁/𝑚
0 −1119534 1119534
Atau,
2,505 −1,5 0
→ 𝐾 = 1119534 [ −1,5 3 −1] 𝑁/𝑚
0 −1 1
poin (a):
Analisis Eigenvalue:
Asumsikan bahwa
Sehingga,
(𝐾 − 𝑀𝜔2 )𝜙𝑠𝑖𝑛𝜔𝑡 = 0
|𝐾 − 𝑀𝜔2 | = 0
Misal
100𝜔2
𝐵=
1119534
Maka,
2,505 − 2𝐵 −1,5 0
| −1,5 3−𝐵 −1 | = 0
0 −1 1 − 2𝐵
Dengan Kalkulator:
𝐵1 = 3,632
𝐵2 = 0,912
𝐵3 = 0,208
Maka,
1119534
𝜔12 = = 3082,416 → 𝜔1 = 55,52 𝑟𝑎𝑑/𝑠
100(3,632)
1119534
𝜔22 = = 12275,59 → 𝜔2 = 110,80 𝑟𝑎𝑑/𝑠
100(0,912)
1119534
𝜔32 = = 53823,73 → 𝜔3 = 232 𝑟𝑎𝑑/𝑠
100(0,208)
poin (b):
Analisis Eigenvector:
Mode 1:
(𝐾 − 𝑀𝜔2 )𝜙 = 0
𝜙11 1
𝜙1 = [𝜙21 ] = [1,303]
𝜙31 2,899
Mode 2:
(𝐾 − 𝑀𝜔2 )𝜙 = 0
2804300 − 200𝜔2 −1679300 0 𝜙12
[ −1679300 3358601 − 100𝜔2 −1119534 ] [𝜙22 ] = 0
2 𝜙
0 −1119534 1119534 − 200𝜔 32
𝜙12 1
𝜙
𝜙2 = [ 22 ] = [ 0,208 ]
𝜙32 −0,174
Mode 3:
(𝐾 − 𝑀𝜔2 )𝜙 = 0
𝜙13 1
𝜙3 = [𝜙23 ] = [−4,74]
𝜙13 0,55
1 1 1
𝜙 = [𝜙1 𝜙2 𝜙3 ] = [1,303 0,208 −4,74]
2,899 −0,174 0,55
Contoh Soal Analisis MDOF: Menentukan Gaya Dalam pada Shear Building
Jawab:
Dengan nilai-nilai yang sudah didapat pada nomor 3, normalized mode shapes dapat dicari
dengan menentukan nilai 𝛼 terlebih dahulu.
1
𝛼𝑗 =
𝑇
√∑ 𝑥𝑖𝑗 𝑚𝑖𝑗 𝑥𝑖𝑗
• Untuk mode 1:
Diketahui:
𝑥11 1
𝑥
𝑥1 = [ 21 ] = [1,303]
𝑥31 2,899
Maka,
1
𝛼1 = = 0,0221
√(1)2 (200) + (1,303)2 (100) + (2,899)2 (200)
𝑥11 1 0,0221
𝑥
𝜙1 = 𝛼 [ 21 ] = 0,0221 [1,303] = [0,0288]
𝑥31 2,899 0,0640
• Untuk mode 2:
Diketahui:
𝑥12 1
𝑥2 = [𝑥22 ] = [ 0,208 ]
𝑥32 −0,174
Maka,
1
𝛼2 = = 0,0689
√(1)2 (200) + (0,208)2 (100) + (−0,174)2 (200)
𝑥11 1 0,0689
𝜙1 = 𝛼 [𝑥21 ] = 0,0689 [ 0,208 ] = [ 0,0143 ]
𝑥31 −0,174 −0,012
• Untuk mode 3:
Diketahui:
𝑥13 1
𝑥3 = [𝑥23 ] = [−4,74]
𝑥13 0,55
Maka,
1
𝛼2 = = 0,02
√(1)2 (200) + (−4,74)2 (100) + (0,55)2 (200)
𝑥11 1 0,02
𝜙3 = 𝛼 [𝑥21 ] = 0,02 [−4,74] = [−0,0947]
𝑥31 0,55 0,011
Dan
Besarnya storey displacement, u(t) bisa dicari melalui persamaan gerak berikut.
𝑀Φ𝑌̈(𝑡) + 𝐶Φ𝑌
̇ (𝑡) + 𝐾Φ𝑌(𝑡) = 𝐹(𝑡)
Φ𝑇 𝑀Φ𝑌̈(𝑡) + Φ𝑇 𝐶Φ𝑌
̇ (𝑡) + Φ𝑇 𝐾Φ𝑌(𝑡) = Φ𝑇 𝐹(𝑡)
Φ𝑇 𝑀Φ = 𝐼
Φ𝑇 𝐾Φ = 𝑑𝑖𝑎𝑔(𝜔2 )
1 0 0
Φ𝑇 𝑀Φ = 𝐼 = [0 1 0]
0 0 1
3082,416 0 0
𝑇 02 12275,59 0
Φ 𝐾Φ = 𝑑𝑖𝑎𝑔(𝜔 ) = [ ]
0 0 53823,73
64,021
𝑇 ∗ (𝑡) −12,016
Φ 𝐹(𝑡) = 𝑓 =[ ] sin(2𝜋𝑡) 𝑁
10,988
Maka, persamaan gerak uncoupled untuk setiap mode dapat dituliskan sebagai berikut.
• Untuk mode 1: 1𝑦̈ 1 (𝑡) + 2𝜁𝑚1∗ 𝜔12 𝑦̇ 1 (𝑡) + 3082,416𝑦1 (𝑡) = 64,021sin (2𝜋𝑡)
• Untuk mode 2: 1𝑦̈ 2 (𝑡) + 2𝜁𝑚2∗ 𝜔22 𝑦̇ 2 (𝑡) + 12275,59𝑦2 (𝑡) = −12,016sin (2𝜋𝑡)
• Untuk mode 3: 1𝑦̈ 3 (𝑡) + 2𝜁𝑚3∗ 𝜔32 𝑦̇ 3 (𝑡) + 52823,73𝑦3 (𝑡) = 10,988sin (2𝜋𝑡)
Dimana persamaan-persamaan gerak di atas termasuk dalam getaran Damped SDOF Harmonic
Loading dimana
𝑃𝑜 /𝑘
𝑦(𝑡) = sin (𝜔𝑡 − 𝜙)
√(1 − 𝛽 2 )2 + (2𝜁𝛽)2
Dimana:
2𝜁𝛽
𝜙 = arctan [ ]
1 − 𝛽2
𝜔
𝛽=
𝜔𝑛
Maka,
• Untuk mode 1:
𝜔 2𝜋
𝛽= = = 0,113
𝜔𝑛 √3082,416/1
2(5%)(0,113)
𝜙 = arctan [ ] = 0,0115 𝑟𝑎𝑑
1 − (0,113)2
64,021/3082,416
𝑦1 (𝑡) = sin (2𝜋𝑡 − 0,0115)
√(1 − 0,1132 )2 + (2 × 5% × 0,113)2
• Untuk mode 2:
𝜔 2𝜋
𝛽= = = 0,057
𝜔𝑛 √12275,59/1
2(5%)(0,057)
𝜙 = arctan [ ] = 0,0057 𝑟𝑎𝑑
1 − (0,057)2
−12,016/12275,59
𝑦2 (𝑡) = sin (2𝜋𝑡 − 0,0057)
√(1 − 0,0572 )2 + (2 × 5% × 0,057)2
• Untuk mode 3:
𝜔 2𝜋
𝛽= = = 0,027
𝜔𝑛 √53823,73/1
2(5%)(0,027)
𝜙 = arctan [ ] = 0,0027 𝑟𝑎𝑑
1 − (0,027)2
10,988/53823,73
𝑦3 (𝑡) = sin (2𝜋𝑡 − 0,0027)
√(1 − 0,00272 )2 + (2 × 5% × 0,0027)2
Gaya geser pada tiap lumped mass dapat dicari dengan menggunakan persamaan
𝑉(𝑡) = 𝐾𝑈(𝑡)
dimana K pada persamaan diatas merupakan matriks kekakuan struktur yang didapatkan pada
Nomor 3.
Gaya geser pada lantai dasar bangunan dapat dicari dengan melakukan perkalian gaya geser
yang bekerja pada tiap lantai dengan transpose matriks 1 yang dimensinya NDOF x 1.
1 𝑇
𝑉𝑏𝑎𝑠𝑒 (𝑡) = [1] 𝑉(𝑡)
1
1 𝑇 284,535 −165,947
𝑉𝑏𝑎𝑠𝑒 (𝑡) = [1] {[−253,685] sin (2𝜋𝑡 − 0,0115) + [ 53,114 ] 𝑠in (2𝜋𝑡 − 0,0057)
1 824,981 28,924
43,899
+ [−74,287] sin (2𝜋𝑡 − 0,0027)}
24,154
𝑽𝒃𝒂𝒔𝒆 (𝒕) = 𝟖𝟓𝟓, 𝟖𝟑𝟏𝐬𝐢𝐧 (𝟐𝝅𝒕 − 𝟎, 𝟎𝟏𝟏𝟓) − 𝟖𝟑, 𝟗𝟎𝟖 𝐬𝐢𝐧 𝒔𝐢𝐧 (𝟐𝝅𝒕 − 𝟎, 𝟎𝟎𝟓𝟕)
− 𝟔, 𝟐𝟑𝟑𝐬𝐢𝐧 (𝟐𝝅𝒕 − 𝟎, 𝟎𝟎𝟐𝟕) [𝑵𝒆𝒘𝒕𝒐𝒏]
Gaya geser pada lantai dasar bangunan dapat dicari dengan melakukan perkalian gaya geser
yang bekerja pada tiap lantai dengan transpose matriks h yang merupakan ketinggian masing2
lumped mass (lantai) diukur dari dasar bangunan yang dimensinya NDOF x 1.
4 4
ℎ = [4 + 3,5] 𝑚 = [7,5] 𝑚
4+7 11
Sehingga,
4 𝑇 284,535 −165,947
𝑀𝑜𝑣 (𝑡) = [7,5] {[−253,685] sin (2𝜋𝑡 − 0,0115) + [ 53,114 ] 𝑠in (2𝜋𝑡 − 0,0057)
11 824,981 28,924
43,899
+ [−74,287] sin (2𝜋𝑡 − 0,0027)}
24,154
𝑴𝒐𝒗 (𝒕) = 𝟖𝟑𝟏𝟎, 𝟐𝟗 𝐬𝐢𝐧 (𝟐𝝅𝒕 − 𝟎, 𝟎𝟏𝟏𝟓) +𝟓𝟐, 𝟕𝟑𝟒 𝐬𝐢𝐧 (𝟐𝝅𝒕 − 𝟎, 𝟎𝟎𝟓𝟕)
− 𝟏𝟏𝟓, 𝟖𝟓𝟖𝐬𝐢𝐧 (𝟐𝝅𝒕 − 𝟎, 𝟎𝟎𝟐𝟕) [𝑵𝒆𝒘𝒕𝒐𝒏. 𝒎𝒆𝒕𝒆𝒓]
Contoh Soal Analisis MDOF: Menentukan Maximum Response
Jawab:
Metode SRSS
2 2
𝑈𝑚𝑎𝑥 = √𝑈𝑚𝑎𝑥,1 + 𝑈𝑚𝑎𝑥,2
Metode CQC
𝑛 𝑛
𝑈𝑚𝑎𝑥 = √∑ ∑ 𝑈𝑖 𝑃𝑖𝑗 𝑈𝑗
𝑖=1 𝑗=1
Dengan:
3
8𝜉 2 (1 + 𝑟)𝑟 2
𝑃𝑖𝑗 =
(1 − 𝑟 2 )2 + 4𝜉 2 𝑟(1 + 𝑟)2
𝜉 = 𝑑𝑎𝑚𝑝𝑖𝑛𝑔 𝑟𝑎𝑡𝑖𝑜
𝜔𝑖
𝑟 = , 𝜔𝑖 < 𝜔𝑗
𝜔𝑗
Ketentuan lainnya:
𝑃𝑖𝑗 = 𝑃𝑗𝑖
𝑖 = 𝑗, 𝑃𝑖𝑗 = 1
Contoh perhitungan:
i j ω r Pij Ui Uj Rij Ri U max
0,00046 0,00046 2,116E-07 2,11,E-07 0,000464
1 1 1 0,0006 0,0006 0,00000036 3,60,E-07 0,000600
0,00134 0,00134 1,7956E-06 1,80,E-06 0,001340
1 55,52
0,00046 -6,8E-05 -5,82163E-10
2 0,501083 0,018611 0,0006 -1,4E-05 -1,56335E-10
0,00134 0,000012 2,9927E-10
-6,8E-05 0,00046 -5,82163E-10 4,04,E-09
1 0,501083 0,018611 -1,4E-05 0,0006 -1,56335E-10 3,97,E-11
0,000012 0,00134 2,9927E-10 4,43,E-10
2 110,8
-6,8E-05 -6,8E-05 4,624E-09
2 1 1 -1,4E-05 -1,4E-05 1,96E-10
0,000012 0,000012 1,44E-10
𝑖 = 1, 𝑗 = 2
0,00046 −0,000068
𝑈1 = [ 0,0006 ] , 𝑈2 = [−0,000014]
0,00134 0,000012
55,52
𝑟= = 0,50108
110,8
3
8 × 5%2 × (1 + 0,50108)0,501082
𝑃12 = = 0,018611
(1 − 0,501082 )2 + 4 × 5%2 × 0,50108(1 + 0,50108)2
Hitung Rij
𝑅𝑖𝑗 = 𝑈𝑖 𝑃𝑖𝑗 𝑈𝑗
Hitung Ri
𝑅1 = 𝑅11 + 𝑅12
2,116 × 10−7 −5,821 × 10−10 2,11 × 10−7
[ 3,6 × 10−7 ] + [−1,563 × 10−10 ] = [ 3,6 × 10−7 ]
1,796 × 10−6 2,993 × 10−10 1,8 × 10−6
Hitung Umax
𝑈𝑚𝑎𝑥 = √𝑅1 + 𝑅2
Metode SRSS
2 2
𝑉𝑚𝑎𝑥 = √𝑉𝑚𝑎𝑥,1 + 𝑉𝑚𝑎𝑥,2
√284,5352 + 165,947 2
𝟑𝟐𝟗, 𝟑𝟗𝟏
2 2 = [𝟐𝟓𝟗, 𝟏𝟖𝟔] (𝑵𝒆𝒘𝒕𝒐𝒏)
√253,685 + 53,114
𝟖𝟐𝟓, 𝟒𝟖𝟖
√824,9812 + 28,9242
[ ]
Metode CQC
i j ω r Pij Vi Vj Rij Ri V max
284,535 284,535 80960,166 80081,381 326,712
1 1 1 -253,685 -253,685 64356,079 64105,306 258,216
824,981 824,981 680593,650 681037,749 826,026
1 55,520
284,535 -165,947 -878,785
2 0,501083032 0,018611 -253,685 53,114 -250,773
824,981 28,924 444,099
-165,947 284,535 -878,785 26659,622
1 0,501083032 0,018611 53,114 -253,685 -250,773 2570,324
28,924 824,981 444,099 1280,697
2 110,800
-165,947 -165,947 27538,407
2 1 1 53,114 53,114 2821,097
28,924 28,924 836,598
c. Base Shear (Vbase)
Nilai Vbase max dari persamaan:
1𝑇
𝑉𝑏𝑎𝑠𝑒,𝑚𝑎𝑥 = [1] 𝑉𝑚𝑎𝑥
1
Metode ABS
1𝑇 450,482
𝑉𝑏𝑎𝑠𝑒,𝑚𝑎𝑥 = [1] × [306,799] = 𝟏𝟔𝟏𝟏, 𝟏𝟖𝟔 (𝑵𝒆𝒘𝒕𝒐𝒏)
1 853,906
Metode SRSS
1𝑇 329,391
𝑉𝑏𝑎𝑠𝑒,𝑚𝑎𝑥 = [1] × [259,186] = 𝟏𝟒𝟏𝟒, 𝟎𝟔𝟒𝟕 (𝑵𝒆𝒘𝒕𝒐𝒏)
1 825,488
Metode CQC
1𝑇 326,712
𝑉𝑏𝑎𝑠𝑒,𝑚𝑎𝑥 = [1] × [258,216] = 𝟏𝟒𝟏𝟎, 𝟗𝟓𝟒𝟑 (𝑵𝒆𝒘𝒕𝒐𝒏)
1 826,026
Metode SRSS
4 𝑇 329,391
𝑉𝑏𝑎𝑠𝑒,𝑚𝑎𝑥 = [7,5] × [259,186] = 𝟏𝟐𝟑𝟒𝟏, 𝟖𝟐𝟑𝟔 (𝑵𝒆𝒘𝒕𝒐𝒏. 𝒎𝒆𝒕𝒆𝒓)
11 825,488
Metode CQC
4 𝑇 326,712
𝑉𝑏𝑎𝑠𝑒,𝑚𝑎𝑥 = [7,5] × [258,216] = 𝟏𝟐𝟑𝟐𝟗, 𝟕𝟓 (𝑵𝒆𝒘𝒕𝒐𝒏. 𝒎𝒆𝒕𝒆𝒓)
11 826,026
Perbandingan hasil
Metode U max (m) V storey (N) V base (N) M ov (Nm)
0,00053 450,482
ABS 0,00061 306,799 1611,186 13495,876
0,00135 853,905
0,00046 329,391
SRSS 0,00060 259,186 1414,065 12341,824
0,00134 825,488
0,00046 326,712
CQC 0,00060 258,216 1410,954 12329,754
0,00134 826,026