## ANALISIS RESPONS DINAMIK PELAT LANTAI BANGUNAN ORTOTROPIK DENGAN DUA PENGAKU YANG DIBERI BEBAN LEDAKAN
Meilani 1 ; Sofia W. Alisjahbana 2
1 Civil Engineering Department, Faculty of Engineering, Binus University Jl. K. H. Syahdan No. 9, Palmerah, Jakarta Barat 11480 meilani@binus.edu
2 Rektor Universitas Bakrie, Jl. H. R.Rasuna Said Kav. C-22, Jakarta 12920 wangsadi@indosat.net.id

## ABSTRACT
Bombings occurred in Indonesia mostly attack high-rise buildings. This fact should be a consideration for the structural engineer when designing high-rise building structures since the design of the structure does not include bomb blast burden in the calculation. This article discusses the dynamic response of buildings with of two plate stiffeners which are then given blast loads. The purpose of this study is to determine the effects of damping and the boom duration of the dynamic response of the plate. The plate is modeled damped orthotropic rectangular with semirigid placement on all sides. A decrease in plate element equilibrium equation is based on moment balance and forces on the plate. Meanwhile, the plate motion equations are derived using Newton's second law. The natural frequency of the plate system is obtained using the Modified Bolotin Method (MBM) as the numerical calculation process is completed with Microsoft Excel 2003 and Mathematica 5.2. Plates will be given damps as much as 0%, 5% and 10%. Each plate with each damp will be given blast loads for 1 ms, 2 ms, 10 ms and 20 ms. The research results shows that the greater the damp is used, the smaller system dynamic response will be. The longer the blast load works, the greater system dynamic response will be.
Keywords: Modified Bolotin Method, orthotropic, stiffener, blast load

## ABSTRAK
Aksi pemboman yang selama ini terjadi di Indonesia kebanyakan menyerang gedung bertingkat tinggi. Fakta ini hendaknya menjadi pertimbangan bagi para insinyur struktur ketika merancang struktur bangunan bertingkat tinggi karena selama ini perancangan struktur tidak memasukkan beban ledakan bom dalam perhitungan. Artikel ini membahas respons dinamik pelat bangunan yang diberi dua buah pengaku berupa balok anak  yang kemudian  diberi beban ledakan. Adapun tujuan dari penelitian ini adalah untuk mengetahui pengaruh redaman dan lamanya ledakan terhadap respons dinamik pelat. Pelat dimodelkan sebagai pelat persegi panjang ortotropik teredam dengan perletakan semirigid pada semua sisinya. Penurunan persamaan keseimbangan elemen pelat dilakukan berdasarkan keseimbangan momen dan gaya yang bekerja pada pelat. Sementara persamaan gerak pelat diturunkan menggunakan hukum Newton II. Frekuensi alami sistem pelat dapat dicari menggunakan Modified Bolotin Method (MBM) dan proses perhitungan numerik diselesaikan dengan bantuan program Microsoft Excel 2003 dan Mathematica 5.2. Pelat tersebut akan diberikan redaman 0%, 5% dan 10 %. Setiap pelat dengan masing-masing redaman akan diberikan beban ledakan selama 1 ms, 2 ms, 10 ms dan 20 ms. Hasil penelitian menunjukkan bahwa semakin besar redaman yang digunakan, respons dinamik system akan semakin kecil dan semakin lama beban ledakan bekerja, respons dinamik system akan semakin besar.
Kata kunci: Modified Bolotin Method, ortotropik, pengaku, beban ledakan

## PENDAHULUAN
Penelitian mengenai pelat sudah dilakukan sejak tahun 1766, berupa penelitian dari  perilaku pelat. Seiring dengan berkembangnya kemajuan jaman penelitian tentang pelat berkembang pesat. Penelitian tentang perilaku pelat bangunan akibat beban dinamik berupa beban ledakan berkembang di Indonesia pasca tahun 2000, di mana aksi pemboman oleh terorisme terjadi di mana-mana. Teroris banyak mengincar gedung bangunan tinggi untuk diledakkan sehingga masalah ini harus dijadikan bahan pertimbangan dalam melakukan desain suatu bangunan. Selama ini yang ramai diperbincangkan adalah ketahanan bangunan terhadap gempa, tetapi menilik situasi dan kondisi di Indonesia akhir- akhir ini, mungkin sudah saatnya kita mulai memperbincangkan ketahanan gedung terhadap beban dinamik yang berupa ledakan. Pelat lantai bangunan pada penelitian ini berupa pelat persegi panjang orthotropik teredam dengan perletakan semirigid pada keempat sisinya yang diperkaku dengan balok anak pada salah satu arah sumbunya. Persamaan gerak pelat persegi panjang orthotropik yang diberi beban ledakan berupa persamaan diferensial orde empat kemudian dapat diselesaikan dengan mencari solusi homogen dan solusi partikuler. Pada pelat tersebut akan diberikan redaman sebesar 0%, 5% dan 10 %. Untuk setiap pelat dengan masing-masing redaman akan diberikan beban ledakan selama 1 ms, 2 ms, 10 ms dan 20 ms. Kemudian akan dicari pengaruh redaman dan lamanya ledakan terhadap respons dinamik pelat.

## METODE

## Permodelan Pelat Lantai Bangunan
Pelat lantai bangunan dimodelkan sebagai pelat orthotropik dengan perletakan semirigid pada keempat sisinya, berukuran 5.05 m x 7.4 m yang diberi dua pengaku pada salah satu arah sumbunya berupa balok anak berukuran 30 cm x 60 cm (Gambar 1).
Gambar 1 . Pelat dengan dua pengaku.
Adapun properties Beton Pelat dan Balok Anak terdapat pada Tabel 1.
Tabel 1 Properties Beton Pelat dan Balok Anak
Adapun persamaan gerak pelat orthotropik dapat dinyatakan sebagai berikut:
(1) (2)
(3) Properties pelat dengan dua pengaku terdapat pada Tabel 2.

## Tabel 2
Properties Pelat dengan Dua Pengaku
Syarat batas yang berlaku untuk pelat persegi panjang dengan perletakan semi rigid di keempat sisinya adalah:
1. Sepanjang x =0 dan x =a:
(4)
2. Sepanjang y =0 dan y =b:
(5)

## Beban ledakan
Beban ledakan akibat ledakan bom dapat dilukiskan sebagai berikut (Gambar 2).
Gambar 2. Hubungan antara tekanan akibat ledakan dengan waktu.
Untuk memperoleh solusi total sistem perlu didefinisikan terlebih dahulu fungsi beban yang bekerja pada sistem tersebut yang dapat dinyatakan sebagai berikut:
(6) di mana:
x (t) = fungsi gerak beban yang bekerja di arah x y (t) = fungsi gerak beban yang bekerja di arah y x (t)= x o dan y (t) = y o = posisi beban pada saat t = t o ( x o, y o) P (t) = beban impuls yang bekerja pada saat t
(7)
Dalam hal ini P max diambil sebesar 1.3 MPa yang bekerja di posisi
.
Bila persamaan (4.3.2) disubstitusikan ke dalam persamaan (4.3.1) akan diperoleh persamaan sebagai berikut:
(8)

## Analisis Umum
Pada getaran bebas pelat ( free vibration ) tanpa redaman dimana harga γ =0 defleksi transversal pelat dapat dinyatakan sebagai berikut:
W ( x , y ,t)=W( x , y ) sin ω t (9)
Dalam persamaan (9) di atas: W( x , y )  = fungsi posisi (spatial function) ω = frekuensi alami sistem (natural frequency system)
Bila persamaan (9) disubstitusikan ke dalam persamaan (1) akan dihasilkan persamaan:
D x ∂ 4 W+2B ∂ 4 W  +D y ∂ 4 W - ρ h ω 2 W = 0 ∂ x 4 ∂ x 2 ∂ y 2 ∂ y 4 (10)
Agar solusi dari persamaan (4) dan persamaan (5) dapat memenuhi syarat batas sesuai dengan persamaan (10), defleksi transversal pelat harus dinyatakan dalam deret Fourier ganda (solusi Navier) yang dapat dinyatakan sebagai berikut:
(11)
di mana:
A mn = koefisien amplitudo yang ditentukan dari kondisi awal (initial condition) m п = wave number di arah x a n п = wave number di arah y b m = indeks integer, pola getar yang terjadi ( mode ) arah x n = indeks integer, pola getar yang terjadi ( mode ) arah y a = panjang pelat di arah x b = panjang pelat di arah y
Jika persamaan (11) disubstitusi ke dalam persamaan (10)  diperoleh frekuensi alami sistem pelat dengan perletakan sendi pada keempat sisinya, yaitu:
(12)
Frekuensi alami sistem untuk pelat persegi panjang dengan perletakan jepit, bebas, semi rigid dapat dicari dengan menganalogikan pelat persegi panjang tersebut sebagai pelat yang memiliki perletakan sendi pada keempat sisinya. Pola getar arah x (m) diganti dengan p dan pola getar di arah y (n) diganti dengan q, dimana p dan q adalah bilangan real yang diperoleh dari persamaan transendetal.
Dengan demikian frekuensi alami sistem dengan perletakan tidak simetris dapat dinyatakan sebagai berikut:
(13)
di mana: p п = wave number di arah x a q п = wave number di arah y b p,q = bilangan real  (untuk perletakan  bukan sendi di kedua sisi yang berhadapan) yang dapat diselesaikan dengan masalah tipe Levy dari dua buah  persamaan auxiliary . Cara ini dikenal juga dengan nama Modified Bolotin Method  (Pevzner, 2000).

## Respons Dinamik Pelat
Respons dinamik pelat dapat ditemukan dengan menggunakan metode pemisahan variabel ( separation of variable method ). Dengan demikian solusi partikuler dari persamaan (1) dapat dinyatakan sebagai:
(14)
di mana T mn (t) = fungsi waktu yang akan diperoleh dengan analisis lebih lanjut.
Persamaan diferensial untuk koefisien fungsi T mn (t) dapat dinyatakan sebagai:
(15)
dimana Q mn adalah faktor normalisasi.
Solusi partikuler dari fungsi T mn (t) bisa dituliskan dalam bentuk integral Duhamel (Paz,1987), yaitu:
Solusi umum sistem akibat beban dinamik transversal dapat dinyatakan sebagai:
(17)

## HASIL DAN PEMBAHASAN

## Frekuensi Alami Sistem
Frekuensi alami sistem dihitung untuk 5 mode pertama di arah x dan 5 mode kedua di arah y . Harga frekuensi alami untuk berbagai mode m dan n dua pengaku tersaji dalam Tabel 3.

## Tabel 3
Harga Frekuensi Alami untuk Berbagai Mode m dan n Dua Pengaku

## Respons Dinamik
Pelat dengan dua pengaku diberi beban ledakan di posisi
dengan lama ledakan td = 1 ms, 2 ms, 10 ms, 20 ms dan diberi redaman 0 %, 5%, 10 %. Hasilnya disajikan dalam Tabel 4 – 6 i bawah ini:

## Tabel 4
Respons Dinamik Pelat Dua Pengaku Dengan Redaman 0%
1 ms 2 ms 10 ms 20 ms Defleksi dinamik 0.0000563048 0.0001988320 0.0029730300 0.0054130400 max (m) 0.0441403000 0.0865300000 0.2749750000 0.3406770000 terhadap fungsi waktu t 8,896.49 32,493.00 140,067.00 298,135.00 terhadap sumbu y 7,188.91 28,452.90 503,627.00 766,515.00 terhadap sumbu x 2,846.10 11,156.80 140,092.00 163,342.00 terhadap fungsi waktu t 2,736,807.57 5,469,148.94 15,564,916.94 19,943,891.79 terhadap sumbu y 1,006,468.71 2,272,659.39 16,761,331.59 19,508,887.43 terhadap sumbu x 985,021.00 2,230,950.15 16,768,857.77 19,521,504.90 terhadap fungsi waktu t 1,262.97 4,667.24 23,334.00 27,974.70 terhadap sumbu y 1,356.30 5,373.01 95,983.00 103,828.00 terhadap sumbu x 402.90 1,583.26 22,417.40 27,959.50 terhadap fungsi waktu t 195,405.00 408,252.00 918,189.00 1,144,038.25 terhadap sumbu y 182,455.00 399,836.00 1,550,306.79 2,442,681.99 terhadap sumbu x 65,163.80 148,546.00 886,625.00 1,075,775.06 terhadap fungsi waktu t 60.34 212.05 3,505.09 6,276.85 terhadap sumbu y 109.42 433.61 8,025.84 13,143.80 terhadap sumbu x 888.96 3,476.07 37,946.10 66,680.20 terhadap fungsi waktu t 51,044.40 100,110.00 319,437.00 395,216.00 terhadap sumbu y 26,649.00 51,582.00 204,425.00 343,315.00 terhadap sumbu x 655,922.00 1,358,931.06 3,013,645.84 3,908,653.86 terhadap fungsi waktu t 793.83 2,740.63 44,981.20 70,575.00 terhadap sumbu y 2,064.73 8,180.27 145,578.00 203,149.00 terhadap sumbu x 257.20 997.07 4,984.62 48,439.00 terhadap fungsi waktu t 128,838.00 261,510.00 689,962.00 961,230.00 terhadap sumbu y 424,016.00 851,171.00 2,791,172.70 3,052,094.23 terhadap sumbu x 112,918.00 232,231.00 447,086.00 507,146.00
Qx (N) Sebelum ledakan Sebelum ledakan Sesudah ledakan 0% 2 td Pengaku Damping ratio Sesudah ledakan Sesudah ledakan Sebelum ledakan Sebelum ledakan Sesudah ledakan Qy (N) Sebelum ledakan Sesudah ledakan Mx (Nm) My (Nm)

## Tabel 5
Respons Dinamik Pelat Dua Pengaku Dengan Redaman 5%
1 ms 2 ms 10 ms 20 ms Defleksi dinamik 0.0000548480 0.0001885210 0.0026619000 0.0046132800 max (m) 0.0408259000 0.0757316000 0.2512600000 0.2540980000 terhadap fungsi waktu t 8,715.99 31,208.20 129,044.00 240,197.00 terhadap sumbu y 7,140.55 28,073.20 474,707.00 715,245.00 terhadap sumbu x 2,818.87 10,945.70 129,066.00 150,332.00 terhadap fungsi waktu t 2,464,585.23 4,582,442.59 8,146,822.50 14,454,993.65 terhadap sumbu y 991,100.00 2,220,735.11 6,825,258.95 11,259,926.53 terhadap sumbu x 970,915.00 2,182,544.13 8,115,879.20 11,681,139.94 terhadap fungsi waktu t 1,239.86 4,502.24 21,682.60 26,056.10 terhadap sumbu y 1,347.30 5,302.08 90,221.30 100,133.00 terhadap sumbu x 399.42 1,556.23 20,937.10 26,041.90 terhadap fungsi waktu t 183,319.00 316,489.00 754,406.00 756,288.00 terhadap sumbu y 177,846.00 385,029.00 1,270,499.32 1,540,132.60 terhadap sumbu x 64,303.80 145,500.00 499,098.00 644,591.00 terhadap fungsi waktu t 58.73 200.69 3,142.28 5,342.73 terhadap sumbu y 108.73 428.19 7,590.99 12,215.70 terhadap sumbu x 879.63 3,403.70 34,321.40 58,586.80 terhadap fungsi waktu t 47,273.00 87,689.30 292,157.00 295,149.00 terhadap sumbu y 26,500.10 46,902.90 107,639.00 257,705.00 terhadap sumbu x 458,671.00 568,060.00 2,934,198.24 3,103,302.94 terhadap fungsi waktu t 771.44 2,583.64 37,965.30 65,609.50 terhadap sumbu y 2,050.96 8,071.60 136,658.00 189,108.00 terhadap sumbu x 253.81 970.87 4,019.42 41,005.50 terhadap fungsi waktu t 100,845.00 169,235.00 539,717.00 559,732.00 terhadap sumbu y 413,803.00 784,230.00 1,763,487.31 2,525,976.51 terhadap sumbu x 75,718.40 94,203.80 250,184.00 307,691.00
Qx (N) Sebelum ledakan Sebelum ledakan Sesudah ledakan 5% 2 td Pengaku Damping ratio Sesudah ledakan Sesudah ledakan Sebelum ledakan Sebelum ledakan Sesudah ledakan Qy (N) Sebelum ledakan Sesudah ledakan Mx (Nm) My (Nm)

## Tabel 6
Respons Dinamik Pelat Dua Pengaku dengan Redaman 10%
1 ms 2 ms 10 ms 20 ms Defleksi dinamik 0.0000534347 0.0001788170 0.0023857500 0.0039706900 max (m) 0.0407047000 0.0730750000 0.0892035000 0.2017350000 terhadap fungsi waktu t 8,540.59 29,994.10 119,343.00 194,730.00 terhadap sumbu y 7,092.69 27,701.20 448,505.00 670,169.00 terhadap sumbu x 2,792.03 10,740.30 119,363.00 138,950.00 terhadap fungsi waktu t 1,401,530.15 2,581,636.65 5,452,933.27 11,046,091.36 terhadap sumbu y 976,093.00 2,170,408.65 4,970,130.21 8,998,345.82 terhadap sumbu x 957,061.00 2,135,344.76 4,623,665.84 6,333,577.69 terhadap fungsi waktu t 1,217.38 4,346.02 20,234.30 24,372.40 terhadap sumbu y 1,338.38 5,232.50 85,017.24 96,742.90 terhadap sumbu x 395.98 1,529.91 19,618.60 24,361.60 terhadap fungsi waktu t 157,110.00 261,310.00 456,811.00 518,031.00 terhadap sumbu y 173,437.00 371,040.00 948,665.00 1,016,337.20 terhadap sumbu x 63,449.70 142,495.00 334,487.00 376,464.00 terhadap fungsi waktu t 57.17 190.00 2,820.33 4,611.67 terhadap sumbu y 108.05 422.87 7,194.34 11,423.80 terhadap sumbu x 870.43 3,333.34 31,166.60 51,338.90 terhadap fungsi waktu t 47,050.80 84,609.20 193,971.00 234,718.00 terhadap sumbu y 26,336.40 45,117.80 84,601.80 86,411.30 terhadap sumbu x 303,144.00 189,900.00 2,395,456.58 2,635,674.09 terhadap fungsi waktu t 749.76 2,436.72 32,302.80 61,263.90 terhadap sumbu y 2,037.31 7,964.97 128,616.00 177,079.00 terhadap sumbu x 250.48 945.48 3,235.21 34,974.90 terhadap fungsi waktu t 78,241.40 129,219.00 228,259.00 364,033.00 terhadap sumbu y 77,656.70 249,600.00 1,404,283.78 1,629,911.19 terhadap sumbu x 53,685.50 60,701.30 153,556.00 177,909.00
Qx (N) Sebelum ledakan Sebelum ledakan Sesudah ledakan 2 td Pengaku Damping ratio 10% Sesudah ledakan Sesudah ledakan Sebelum ledakan Sebelum ledakan Sesudah ledakan Qy (N) Sebelum ledakan Sesudah ledakan Mx (Nm) My (Nm)

## Gaya Dalam Pelat
Berikut adalah gaya dalam pelat dengan dua pengaku diberi beban ledakan di posisi dengan lama ledakan td = 1 ms, 2 ms, 10 ms, 20 ms dan diberi redaman 0 %, 5%, 10 % (Gambar 3 – 5).

## Gaya Dalam Pelat Pada Saat  t = 25 ms
Pada saat 25 ms setelah terjadi ledakan, diamati respons dinamik di tengah bentang pelat dua pengaku redaman 0 % dengan posisi beban berada di x = a /4 dan y = b /4 (Lihat Tabel 7 – 10).

## KESIMPULAN
Berdasarkan analisis respons dinamik lantai pelat bangunan persegi panjang orthotropik dengan bantuan program Mathematica 5.2 dapat disimpulkan dua hal. Pertama, semakin besar redaman yang digunakan, semakin kecil respons dinamik sistem. Kedua, semakin lama beban ledakan bekerja, respons dinamik sistem akan semakin besar.

## DAFTAR PUSTAKA
Alisjahbana, S.W. (2002). Dinamika Struktur Pelat . Jakarta: Jurusan Teknik Sipil Universitas Tarumanagara.
Elishakoff, I. B. (1974). Vibration Analysis of Clamped Square Orthotropic Plate. AIAA Journal, 12 (7), 921-924.
Kadid, Abdelkrim. (2008). Stiffened Plates Subjected To Uniform Blast Loading. Algeria: Journal of Civil Engineering and Managemen t, pp. 155-161.
Meilani. (2010). Respons Dinamik Lantai Bangunan Dengan Pengaku Akibat Beban Ledakan . Tesis Pascasarjana Teknik Sipil. Jakarta: Fakultas Teknik Universitas Tarumanagara
Pevzner, Pavel., Weller, Tanchum & Berkovitz, Abraham. (2000). Further Modification of Bolotin Method in Vibration Analysis of Rectangular Plates. AIAA Journal, 38 (9), 1725-1729.

## APPENDIX
Gambar 3 . Gaya dalam pelat dua pengaku dengan redaman 0%.
Gambar 4 . Gaya dalam pelat dua pengaku dengan redaman 5%.
Gambar 5. Gaya dalam pelat dua pengaku dengan redaman 10%.

## Tabel 7
Momen Arah  X Pelat Dua Pengaku dengan Redaman 0% dengan Posisi Beban di x= a/4 dan y =  b/4

## Tabel 8
Momen Arah Y Pelat Dua Pengaku dengan Redaman 0% dengan Posisi Beban di x= a/4 dan y =  b/4

## Tabel 9
Gaya Geser Arah X Pelat Dua Pengaku dengan Redaman 0% dengan Posisi Beban di x= a/4 dan y =  b/4

## Tabel 10
Gaya Geser Arah Y Pelat Dua Pengaku dengan Redaman 0% dengan Posisi Beban di x= a/4 dan y =  b/4