DINAMIKA STRUKTUR
i
DINAMIKA STRUKTUR
i
KATA PENGANTAR Puji syukur kehadirat Allah SWT yang telah melimpahkan rahmatNya sehingga kami dapat menyelesaikan Buku Ajar Mata Kuliah Dinamika Struktur ini. Buku ajar ini merupakan bagian dari media bahan ajar yang dimaksudkan untuk meningkatkan pemahaman mahasiswa terhadap materi perkuliahan yang disampaikan, khususnya mata kuliah Dinamika Struktur, Jurusan Teknik Sipil, Fakultas Teknik, Universitas Brawijaya, Malang. Buku ajar ini disusun dalam enam bab. Bab I memperkenalkan konsepkonsep dasar mengenai dinamika struktur, respon struktur terhadap beban dinamik, analisa dinamis pada struktur, serta derajat kebebasan. Bab II membahas sistem berderajat kebebasan tunggal (SDOF) yang meliputi pemodelan parameter, pemodelan matematis, free body diagram, diagram, dan persamaan gerak dari suatu struktur. Getaran bebas sistem SDOF untuk kondisi tak teredam dan teredam dibahas pada bab III. Selain itu juga dijelaskan mengenai eksperimen penentuan frekuensi alami dasar dan faktor damping, serta getaran bebas dengan coulomb damping dari sebuah sistem SDOF. Sistem SDOF terhadap gerak harmonis untuk sistem tak teredam dan sistem dengan redaman viskous dijelaskan pada bab IV. Bab V membahas respon sistem SDOF terhadap bentuk spasial dari eksitasi, meliputi respon sistem redaman viskous untuk step input ideal, respon sistem tak teredam pada rectangular pulse dan pembebanan ram, serta impuls dengan durasi pendek, unit respon impuls. Bab VI dibahas tentang respon sistem SDOF pada eksitasi dinamis dengan metode integral duhamel. Akhirnya, pada bab VIII dan IX membahas mengenai sistem berderajat kebebasan banyak (MDOF). Kami menyadari bahwa dalam penyusunan buku ajar ini masih terdapat banyak kekurangan. Oleh karena karena itu, kritik dan saran yang konstruktif sangat kami harapkan. Semoga buku ini dapat memberikan manfaat kepada siapapun yang ingin mengkaji dinamika struktur.
Hormat kami,
Penyusun
DINAMIKA STRUKTUR
ii
DAFTAR ISI
KATA PENGANTAR ......................................... .............................................................. ........................................... ............................... ......... i DAFTAR ISI .................................................... .......................................................................... ........................................... ................................. ............ ii BAB I PENDAHULUAN .......................................... ............................................................... ........................................... ........................ 1 1.1
Pendahuluan mengenai Dinamika Struktur ...................................................... ......................................................1
1.2
Analisa Dinamis pada Struktur Struk tur ............... ........ ............... ............... .............. ............... ............... ............... ............... ............. ......2
1.3
Derajat Kebebasan (Degrees of Freedom) ....................................................... .......................................................4
BAB II SISTEM BERDERAJAT KEBEBASAN TUNGGAL TUNGGAL (SDOF) ............... 6 2.1
Pemodelan Parameter.................................................................................. ..................................................................................... ...6
2.2
Pemodelan Matematis ............................................................................. ..................................................................................... ........7
2.3
Free Body Diagram ............................................................ .......................................................................................... ..............................9
2.4
Persamaan Gerak (Equation of Motion) ......................................................... .........................................................10
2.4.1 Aplikasi Aplikas i dari Hukum Newton Newt on Pada Model-model Lumped Parameter .......... ........ ..10 2.4.2 Prinsip D’Alembert ............................................................... 12 ...................................................................................... .......................12 2.4.3 Solusi Persamaan Gerak SDOF Tak Teredam (Undamped)............................. 15 2.4.4 Persamaan Gerak G erak SDOF Teredam Te redam (Damped) (Dam ped) .................................................. 20
BAB III GETARAN BEBAS SISTEM SDOF ........................................... .................................................... ......... 21 3.1
2 Pendahuluan ................................................................. ................................................................................................ ................................21 .1
3.2
Getaran Bebas Pada Pad a Sistem SDOF Tak Teredam Tereda m (Undamped) (Undamped ).......................... 22
3.3
Getaran Bebas Pada Sistem SDOF Teredam (Damped) (Dam ped) .................................... 23
3.4
Eksperimen Penentuan dari Frekuensi Alami Das ar dan Faktor Damping dari sebuah sistem sist em SDOF SDO F ......................................................... .......................................................................................... ................................. 26
3.5
Getaran Bebas dari sebuah sistem SDOF dengan
Dampin g ............... ............... Coloumb Damping
32
BAB IV RESPON SISTEM SDOF TERHADAP GERAK HARMONIS............ 35 4.1
Respon Sistem SDOF Tak Teredam Terhadap Terhada p Gerakan Harmonis .................... 35
4.2
Respon Sistem SDOF SD OF Redaman Viskous Viskou s Terhadap Gerakan Harmonis ............. ............. 39
BAB V Respon Sistem SDOF Terhadap Bentuk Spesial Dari Eksitasi ................ 44 5.1
Respon Dari Sebuah Viscous-Damped System SD OF Untuk Sebuah Step Input 44 yang Ideal................................ ................................................ ................................. ................................. ............................ ............44
5.2
Persamaan Respon dari sebuah Sistem Undamped SDOF pada Rectangular Pulse dan Pembeban Pe mbebanan an Ram ............................................................................ ............................................................................45
DINAMIKA STRUKTUR 5.3
ii
Respon Dari Sistem SDOF Tak Teredam untuk Impuls dengan Durasi Pendek, Unit Respon Impuls ........................................................................................49
BAB VI Respon System SDOF pada Eksitasi Dinamis ........................................ 53 6.1
Metode Integral Duhamel ...................................................................... 53
BAB VII Respons Spektrum ................................................................................. 62 7.1
Bentuk Respons Spektrum ............................................................................62
7.2
Respons Spektrum pada Pondasi yang Bergerak .............................................. 65
7.3
Besaran- Besaran Respons Spektrum................................................................ 66
7.4
Respons Spektrum untuk Perencanaan Elastis.................................................. 68
BAB VIII SISTEM BERDERAJAT KEBEBASAN BANYAK
(MDOF)...... 70
8.1
Sistem MDOF Sederhana................................................................................... 70
8.2
Hukum Newton Kedua pada Sistem MD OF ....................................................70
8.3
Prinsip D’Alembert’s pada Sistem MDOF .......................................................71
8.4
Sistem Massa – Pegas – Redaman .................................................................72
8.5
Koefisien Kekakuan .......................................................................................74
BAB IX GETARAN BEBAS UNTUK SISTEM MDOF .................................... 77 9.1
Sistem MDOF Tak Teredam ...........................................................................77
9.2
Frekuensi Natural dan Pola Normal ................................................................78
9.3
Sifat Ortogonalitas dari Pola Normal ..............................................................79
9.4
Solusi Persamaan Getaran Bebas pada Sistem Tak teredam ............................ 83
9.5
Respon Pada Gedung Akibat Gempa .............................................................84
DINAMIKA STRUKTUR
1
BAB I PENDAHULUAN
Pendahuluan mengenai Dinamika Struktur Secara sederhana dinamik dapat diartikan sebagai variasi atau perubahan terhadap waktu dalam konteks gaya yang bekerja (eksitasi) pada struktur. Beban dinamis dapat berupa variasi besarannya (magnitude), arahnya (direction) atau posisinya (point of application) berubah terhadap waktu. Demikian pula respons struktur terhadap beban dinamik, yaitu lendutan dan tegangan yang dihasilkan juga perubahan-waktu, atau bersifat dinamik. 1.1
(a)
(b)
Gambar 1.1. Balok kantilever dengan (a) beban statis dan (b) beban dinamis.
Pada gambar diatas terlihat balok kantilever dengan dua jenis pembebanan berbeda yaitu beban statis dan dinamis. a. gambar 1.1 (a) menunjukan balok kantilever dengan beban s tatis, responnya dipengaruhi oleh beban P. b. gambar 1.1 (b) menunjukan balok kantilever dengan beban dinamis atau beban yang bervariasi terhadap waktu P(t). Lendutan dan tegangan internal yang timbul dalam kasus beban statis hanya ditimbulkan langsung oleh beban P, sedangkan dalam kasus beban dinamis, percepatan yang dialami oleh balok akibat P(t) menimbulkan gaya inersia yang terdistribusi pada seluruh bagian balok. Lendutan dan tegangan pada balok sangat dipengaruhi pula oleh gaya inersia yang ditimbulkan oleh massa balok ketika mengalami percepatan. Jika pengaruh gaya inersia yang terjadi sangat signifikan, maka perlu dilakukan analisa dinamis. Perbedaan respon untuk beban statis dan dinamis juga dapat dilihat pada gambar 1.2 berikut.
DINAMIKA STRUKTUR
STATIS P
2
DINAMIS P(t)
Gambar 1.2. Balok dengan (a) beban statis dan (b) beban dinamis
Analisa Dinamis pada Struktur Dapat dikatakan bahwa langkah yang paling diperlukan dalam sebuah analisa dinamis adalah pemodelan matematis. Namun secara keseluruhan langkahlangkah dalam analisa dinamis dapat dilihat pada gambar berikut. 1.2
Gambar 1.3. Langkah-langkah dalam analisa dinamis.
DINAMIKA STRUKTUR
3
Model analitis terdiri dari: a. Asumsi sederhana yang dibuat untuk menyederhanakan suatu sistem. b. Gambar dari model analitis tersebut. c. Daftar parameter desain. Model analitis terbagi dalam dua kategori dasar : a. Model berkesinambungan (continues model) b. Model diskrit (discrete-parameter model) Model berkesinambungan (continues model) mempunyai jumlah derajat kebebasan (number of DOF) tak berhingga. Namun dengan proses idealisasi, sebuah model matematis dapat mereduksi jumlah derajat kebebasan menjadi suatu jumlah diskrit.
(a)
(b)
(c)
Gambar 1.4. Model analitis berkesinambungan ( continues) dan diskrit (discrete-parameter ) pada sebuah balok kantilever.
Model berkesinambungan (continues model) pada gambar 1.4(a) menunjukan jumlah derajat kebebasan tak berhingga, model diskrit pada gambar 1.4 (b) dan (c) ditunjukan dengan model massa terkelompok (lumped-mass model) dimana massa terbagi rata dari sistem dianggap sebagai massa titik atau partikel.
DINAMIKA STRUKTUR
1.3
4
Derajat Kebebasan (Degrees of Freedom) Jumlah koordinat bebas yang menetapkan susunan atau posisi sistem pada setiap saat.
Model Struktur
Model SDOF
Model MDOF
Model Struktur
Model SDOF
Model MDOF
Model Struktur
Model SDOF
Model MDOF
DINAMIKA STRUKTUR
Model Struktur
Model Struktur
Model SDOF
Model SDOF
Model MDOF
Model MDOF
Gambar 1.5. Beberapa model struktur dengan derajat kebebasan SDOF (Single Degree of Freedom) dan MDOF (Multiple Degree of Freedom).
5
DINAMIKA STRUKTUR
6
BAB II SISTEM BERDERAJAT KEBEBASAN TUNGGAL (SDOF)
Pemodelan Parameter Komponen-komponen yang merupakan pemodelan himpunan parameter dari sebuah struktur adalah sesuatu yang menghubungkan gaya dengan perpindahan, kecepatan, dan percepatan. Komponen yang menghubungkan gaya dengan perpindahan disebut pegas. Gambar 2.1 menunjukkan idealisasi pegas tak bermassa dan plot gaya dari pegas terhadap regangan. Gaya pegas selalu bekerja sepanjang garis hubung kedua ujung pegas. Hubungan linier antara gaya dan regangan dinyatakan : ……………(2.1) f s = k e 2.1
dimana, k adalah konstanta pegas. Besaran k adalah pound/inc (lb/in) atau N/m. Energi tegangan dinyatakan dengan ……………(2.2) V = ½ (k e 2 )
Gambar 2.1. Gaya-deformasi pada pegas.
dimana energi tegangan dinyatakan sebagai area dibawah kurva f s terhadap e. Model analitis yang paling umum dari redaman dalam analisa dinamika struktur adalah model tahanan dashpot, yang dapat diilustrasikan pada gambar 2.2.
Gambar 2.2. Model tahanan dashpot.
Gaya redaman f D dinyatakan : f D
c(u 2 u
1
)
……………(2.3)
Dari fungsi linear dari kecepatan relatif antara dua ujung dashpot.
DINAMIKA STRUKTUR
7
Konstanta c disebut koefisien viscositas redaman dan besarannya adalah pond/inc/detik atau N/m/detik. Dalam menulis persamaan gerak dari partikel, hukum kedua dari Newton digunakan,
F ma
……………(2.4)
dimana m adalah massa dan a adalah percepatan relatif dari suatu bidang referensi inersia. Besaran massa adalah lb.det/in atau N.det/in. Untuk permasalahan dinamika struktur seringkali sangat berguna untuk memperkenalkan gaya inersia. f l ma
……………(2.5)
Kemudian persamaan 2.4 bisa ditulis sebagai persamaan dinamik yang semisal :
F ' f F 0 l
……………(2.6)
dengan resultan gaya inersia yang ditambahkan pada resultan gaya lain yang bekerja pada partikel.
Pemodelan Matematis Model matematis dalam analisa dinamika struktur mempunyai beberapa elemen sebagai berikut: 2.2
massa m menyatakan massa dan sifat inersia dari struktur pegas k menyatakan gaya balik elastic dan kapasitas ener gy potensial dari struktur redaman c menyatakan sifat geseran dan kehilangan energy dari struktur
gaya pengaruh F(t) menyatakan gaya luar yang bekerja pada sistem struktur sebagai fungsi dari waktu. Namun dalam pembahasan dinamika struktur dengan analisa sederhana pada sistem berderajat kebebasan tunggal, redaman c diabaikan. Beberapa contoh model matematis pada struktur dapat dilihat pada gambar berikut.
DINAMIKA STRUKTUR
8
m K
K EI
m
y
Model Struktur
Model SDOF
m
P(t)
P(t)
Model Matematis
K m
K1
K2
K
Model Struktur
Model SDOF
P(t)
Model Matematis
Gambar 2.3. Model matematis sistem berderajat kebebasan tunggal.
Pada model diatas, massa m dihambat oleh pegas k dan bergerak menurut garis lurus sepanjang satu sumber koordinat. Karakteristik mekanis pegas digambarkan antara gaya Fs pada ujung pegas dan hasil perpindahan y dapat dilihat pada gambar 2.4 (a) sedangkan tiga jenis pegas ditunjukan secara grafis pada gambar 2.4 (b). Fs (gaya) hard spring linier spring soft spring
y Fs
(a)
y (perpindahan)
(b) Gambar 2.4. Hubungan gaya dan perpindahan pada pegas.
Lengkungan pada pegas kuat (hard spring) menyatakan sifat dimana gaya harus memberikan pengaruh lebih besar untuk suatu perpindahan yang diisyaratkan seiring dengan terdeformasinya pegas. Karakteristik garis lurus pada pegas liniear (linear spring) menggambarkan deformasi yang selaras dengan gaya. Konstanta keselarasan antara gaya dan perpindahan dari pegas linier disebut konstanta pegas (spring constant) k . Sedangkan pada pegas lemah (soft spring),
DINAMIKA STRUKTUR
9
pertambahan gaya untuk memperbesar perpindahan cenderung mengecil pada saat deformasi pegas menjadi makin besar. Jika suatu pegas terpasang secara paralel atau seri, maka diperlukan penentuan konstanta pegas ekivalen dari sistem tersebut.
K1
K1
K2
K2
y m
1 y
k e
1
k 1
1 k 2
P
k e k 1 k 2 (a) (b) Gambar 2.5. Kombinasi pegas (a) pegas paralel (b) pegas seri.
Untuk n pegas yang dipasang parallel, konstanta pegas ekivalennya: n
k e
k i i1
……………(2.7)
Sedangkan untuk n pegas yang terpasang seri :
1 k e
n
1
i1
k i
……………(2.8)
Free Body Diagram Salah satu aspek yang penting dalam analisis dinamis adalah menggambar sebuah diagram free body dari sistem yang memungkinkan penulisan besaran matematik dari sistem tersebut. Free Body Diagram (FBD) adalah suatu sketsa dari benda yang dipisahkan dari benda lainnya, dimana semua gaya luar pada benda terlihat jelas. Sebagai contoh dapat dilihat pada gambar berikut: 2.3
K m
P(t)
f s
P(t)
I Gambar 2.6. Free Body Diagram dari sebuah sistem berderajat kebebasan tunggal.
DINAMIKA STRUKTUR
10
Dari gambar free body diagram diatas, menunjukan bahwa massa m yang dipindahkan dengan adanya gaya luar sebesar P(t), dan memberikan gaya pegas sebesar F s=ky serta gaya inersia I . Persamaan Gerak (Equation of Motion) Pada bagian ini persamaan gerakan dari beberapa model lumped parameter akan diturunkan dengan menggunakan hukum Newton atau yang ekivalen, metode gaya D’Alembert. Hal ini akan berlaku sebagai review atas pelajaran sebelumnya pada dinamika dan juga memperkenalkan prosedur yang digunakan dalam menentukan model matematis dari sistem SDOF. 2.4
2.4.1 Aplikasi dari Hukum Newton Pada Model-model Lumped Parameter Untuk menentukan gerak pada sebuah sistem, yaitu mempelajari perpindahan atau kecepatan massa m pada saat t untuk kondisi awal pada saat
t 0 . Hubungan antara perpindahan dan waktu diberikan oleh Hukum Newton Kedua untuk gerak yang ditulis pada persamaan (2.4), dimana F adalah resultan gaya yang bekerja pada partikel massa m dan a adalah resultan percepatan. Persamaan diatas merupakan persamaan vector yang dapat ditulis dalam bentuk ekivalen, dalam besaran komponennya menurut sumbu koordinat.
F ma F ma F ma x
x
y
y
z
z
……………(2.9)
Contoh 2.1
Gunakan hukum Newton untuk menurunkan persamaan gerakan dari sistem pegas sederhana dan dashpot massa di bawah ini. Asumsikan hanya ada gerakan vertikal. Dan asumsikan bahwa pegas linier dengan konstanta pegas k . Abaikan gesekan udara, massa pegas, dan redaman dalam pegas. P(t) adalah gaya yang bekerja pada massa dari luar.
DINAMIKA STRUKTUR
11
Solusi: Tentukan bidang referensi dan koordinat perpindahan. Pilih sumbu x sepanjang garis pergerakan dan tentukan titik acuan awal (misal x = 0) pada lokasi dimana pegas tidak teregang. u adalah perpindahan pada arah x.
Gambar diagram free body dari partikel.
Gunakan hukum Newton yang kedua
Fx mu
……………(2.10)
(catatan : tanda + menunjukkan arah ke bawah dimana u adalah positif untuk arah ke bawah). Dari diagram free body, tentukan gaya-gaya pada bagian kanan persamaan (2.10) ……………(2.11) p fs fd W mu Hubungkan gaya dengan sistem variabel gerakan fs ke ku fd ce
……………(2.12) ……………(2.13)
cu Gabungkan dan susunlah variabel yang tidak diketahui di bagian kanan pada persamaan ……………(2.14) mu cu ku W p(t ) (Catat bahwa ini adalah persamaan diferensial ordiner ordo dua, linier, non homogen dengan koefisien konstan). Persamaan ini bisa disederhanakan dengan pertimbangan sebagai berikut. Perpindahan statis dari bobot w dinyatakan sebagai perpindahan dari massa terukur berhubungan dengan posisi setimbang, statis sebagai u r sehingga
u ur u st
……………(2.15) dimana ust adalah konstan, persamaan (2.14) bisa ditulis sebagai :
DINAMIKA STRUKTUR
12
……………(2.16) mu r cu r kur p(t ) Persamaan (2.16) pada contoh 2.1 bisa dipertimbangkan sebagai persamaan dasar pada dinamika struktur dan teori getaran linier. Akan diperlukan waktu yang lama untuk menetukan solusinya dan aplikasinya pada soal-soal dinamika struktur, baik sistem SDOF maupun MDOF. Pada contoh 2.1, hukum Newton yang kedua digunakan langsung, sehingga tidak ada gaya inersia yang diperlihatkan pada diagram free body.
2.4.2 Prinsip D’Alembert Alternatif pendekatan untuk mendapatkan persamaan gerak adalah penggunaan Prinsip D’Alembert yang menyatakan bahwa sebuah sistem dapat dibuat dalam keadaan keseimbangan dinamis dengan menambahkan sebuah gaya fiktif pada gaya-gaya luar yang disebut sebagai gaya inersia.
mg
y K m
f s = ky
m
I m y
N (b) (a) Gambar 2.7. Sistem berderajat kebebasan tunggal, (a) model matematis dan (b) diagram Free Body.
Penggunaan Prinsip D’Alembert memungkinkan pemakaian persamaan keseimbangan untuk mendapatkan persamaan gerak. Pada gambar free body diagram diatas dapat dilihat bahwa jumlah gaya-gaya pada arah y memberikan persamaan
H f s m y
Dengan:
0 I 0 ky 0 → Persamaan gerak (Equation of Motion)
y = simpangan 2 2 y d y dt = percepatan m = massa k = kekakuan elemen
DINAMIKA STRUKTUR
Satuan: lb
k m
w
g
in lb in sec2
g 386 in
K
sec2
K1
K1 yo m
y m
I W
A
B
Keterangan: A
: Pegas belum dibebani
B
: Pegas dibebani (statis)
C
: Pegas dibebani (dinamis)
Kondisi (B) Statis f s
V 0
fs W 0
m
W fs W k . y o W
C
13
DINAMIKA STRUKTUR
14
Kondisi (C) Dinamis
f s
V 0 fs I W
fs k y o y
m
I
I m. y W k . y o
fs I W
W
k y o y m. y k . y o k . yo k . y m. y
m. y 0
k . y
k .yo
→ Persamaan gerak (Equation of Motion)
Untuk menunjukkan kegunaan gaya inersia dan juga mengilustrasikan fungsi utama eksitasi terdukung atau gerakan dasar, seperti struktur gedung yang akan mengalaminya selama gempa bumi, dapat dilihat pada contoh 2.2 . Contoh 2.2
Gunakan metode gaya D’Alembert untuk menentukan persamaan gerakan dari massa m, asumsikan bahwa gaya redaman pada sistem bisa diwakili dengan viskous dashpot linier seperti yang diperlihatkan pada gambar di bawah. Asumsikan bahwa eksitasi terdukung z(t) diketahui. ketika u = z = 0, pegas belum diregangkan.
Solusi: Gambarkan diagram free body dari massa termasuk gaya inersia bersama dengan gaya sesungguhnya.
DINAMIKA STRUKTUR
15
Tulis persamaan kesetimbangan dinamis
F ' x 0
……………(2.17)
Dari diagram freebody didapat p fs fd mu 0
……………(2.18)
Hubungkan gaya dengan variable gerakan dan sederhanakan mu c(u z ) k (u z ) p
……………(2.19)
Ingat bahwa gaya redaman dan gaya pegas yang dihubungkan dengan gerakan dari massa mempunyai hubungan dengan gerakan yang terdukung. Persamaan (2.19) bisa dituliskan dengan semua nilai yang diketahui dari bagian kanan persamaan. ……………(2.20) mu cu ku cz kz p Persamaan (2.20) adalah persamaan dari gerakan dari perpindahan aktual dari massa yang berada dalam kerangka acuan inersia yakni untuk u(t) ……………(2.21) w u z Dengan mengalihkan m z pada persamaan (2.19) dan menggunakan persamaan (2.21) bisa didapatkan persamaan berikut : ……………(2.22) mw cw kw p m z 2.4.3 Solusi Persamaan Gerak SDOF Tak Teredam (Undamped) Persamaan gerak untuk sistem berderajat kebebasan tunggal tak teredam adalah
m. y 0
k . y
……………(2.23)
Misal solusi: y Acos t
……………(2.24) ……………(2.25)
y Asin t
Kita menganggap bahwa solusi pada persamaan (2.23) adalah persamaan (2.24) y Acos t
A sin t A 2
y
……………(2.26)
y cos t
Substitusikan persamaan (2.24) dan (2.26) kedalam persamaan (2.23) m. y k . y 0 m 2 k 0 0 2 m A. Cos t k ACos t 0 ACos t m 2 k ACos t 0
Sehingga: m 2 k 0 k 2 m
DINAMIKA STRUKTUR
k
16
→ Frekuensi Alami Struktur [rad/dt] ……………(2.27)
m Sebenarnya persamaan (2.25) juga solusi, maka solusi umumnya adalah: ……………(2.28) y ACos t B Sin t ……………(2.29) y A Sin t BCos t Jika dimasukkan masalah kondisi awal (t = 0) yaitu: Perpindahan awal
: y t y 0
y
……………(2.30)
Kecepatan awal
: y t y 0 V o
……………(2.31)
o
Maka substitusi persamaan (2.30) ke dalam persamaan (2.28) didapat: ……………(2.32) A yo Substitusi persamaan (2.31) dan (2.32) ke dalam persamaan (2.29), maka didapat: V ……………(2.33) B o Substitusi persamaan (2.32) dan (2.33) ke dalam persamaan (2.29), maka didapat: V o Cos t Sin t y yo → Solusi Gerak Respons atau
y C Sin t
dengan: 2
C
2 o
tan
V o o
V o
Gambar 2.8. Respon getaran bebas tak teredam.
f
1 T
T
2
2
→
Frekuensi Alami [Siklus/dt]
→
Periode Getar
DINAMIKA STRUKTUR
17
Contoh 2.3 SDOF 200 lb/ft
Model Struktur :
m F(t)
F(t)
E = 30.10 6 psi I = 82,5 in4 W8x24
15 ft
W = 200 x 25 = 5000 lb g = 386 ft/dt 2
y
Model Matematis
K
m
F(t)
m
f s
Free Body Diagram
F(t)
I
Persamaan Kesetimbangan:
I fs F t m. y k . y
F t
12.30.10 2.82,5 6
12 E 2 I K L3 m
W g
(Equation of Motion)
k
15.123
5000 386 10,185.386
m 1
f 2
10,185 lb / in
2
5000 10,185.386 5000
rad / dt
4.46 sps
DINAMIKA STRUKTUR
Latihan.
Jika:
0,001 ft
Simpangan awal y 0
Kecepatan awal y 0 0,1 ft/dt Gaya luar F(t) Gambarkan Respons Struktur!!
18
DINAMIKA STRUKTUR
Konstanta Pegas (Konstanta Elastis) P K . yo P K yo
yo
P
3
P
PL
yo
48 EI P K P 3 PL y
EI yo
48 EI 3 L
o
48 EI
yo P 3
yo
Ph
12 EI P K P 3 Ph y
EI h
12 EI h3
o
12 EI
P yo
EI
yo
K
Pl
3 EI P P
yo L
3
Pl
3
3 EI
3 EI 3
L
19
DINAMIKA STRUKTUR
20
P yo
yo
Ph
EA P P EA K
h
yo
Ph EA
h
2.4.4 Persamaan Gerak SDOF Teredam (Damped) Pada pembahasan sebelumnya telah dijelaskan beberapa cara untuk memeperoleh persamaan gerak untuk SDOF teredam. Struktur yang dimodelisasikan sebagai sistem sederhana dengan redaman-liat (viscousdamping), seperti pada gambar berikut: m
P(t)
P(t)
m
y
K1
K2
K,c K
m I
c
(b) (c) (a) Gambar 2.9. Sistem SDOF teredam, (a) model struktur, (b) model SDOF, dan (c) model matematis.
Free Body Diagram f s f d
I
P(t)
H 0 I f d f s P (t ) I m y my cy ky f d cy P t
f s
ky
→ solusi persamaan gerak
P(t)
DINAMIKA STRUKTUR
21
BAB III GETARAN BEBAS SISTEM SDOF Pendahuluan Pada semua kasus, persamaan gerak sistem linier berderajat kebebasan tunggal mempunyai bentuk ……………(3.1) mu cu ku p(t ) 3.1
Perpindahan dan kecepatan pada saat t = 0 adalah
u(0) uo ,
u (0) u o
……………(3.2)
dimana, uo dan u o adalah perpindahan awal dan kecepatan awal. Persamaan (3.1) dapat ditulis kembali menjadi 2 2 u 2 nu n u n p(t ) k
……………(3.3)
dimana 2
n
k m
dan c
ccr 2k
2m n 2
ccr
km
n
Untuk getaran bebas →P(t) = 0, maka persamaan (3.1) dan (3.3) menjadi: ……………(3.4) mu cu ku 0 2 2 u n u n
u
0
……………(3.5)
n adalah frekuensi alami sudut tak teredam (rad/s), adalah faktor redaman liat dan ccr adalah koefisien redaman kritis. Respon total: u(t )
u p (t) uc (t )
……………(3.6)
Di dalam istilah matematika, penyelesaian umum dari persamaan diferensial terdiri dari penyelesaian sesungguhnya u p(t) dan penyelesaian komplemen/pelengkap uc(t). Untuk memenuhi persamaan (3.4) dan (3.5), maka digunakan asumsi st ……………(3.7) u Ce Dengan mensubstitusikan persamaan (3.7) kedalam (3.5), maka diperoleh
s
2
2 n s n 2 Ce st 0
……………(3.8)
Agar persamaan (3.8) valid untuk semua nilai t , kita harus menentukan s 2 2 n s
n 2 0
……………(3.9)
DINAMIKA STRUKTUR
22
Getaran Bebas Pada Sistem SDOF Tak Teredam (Undamped) Persamaan gerak untuk sistem berderajat kebebasan tunggal (SDOF) tak teredam adalah ……………(3.10) u 2u 0 3.2
n
Dan persamaan karakteristik yang sesuai adalah s2 n
2
0
……………(3.11)
Akar dari persamaan adalah s1,2 i n
i
dimana
-1
……………(3.12)
Jika akar-akar tersebut di substitusikan ke persamaan (3.7), kita mendapat penyelesaian umum ……………(3.13) u C 1ei nt C 2ei nt dengan memperkenalkan persamaman Euler : ……………(3.14) i e cos i sin kita dapat menulis ulang persamaan (3.13) dalam bentuk fungsi trigonometri, yaitu ……………(3.15) u A cos t A sin t 1
n
2
n
dimana A1 dan A2 adalah konstanta real untuk ditentukan dari kondisi awal yaitu persamaan 3.2. Persamaan 3.2 dan 3.15 3.15 mengacu pada u(0) uo A1 u (0) u
o
……………(3.16)
jadi, A2 n ……………(3.17)
u o u uo cos nt
sin nt
n
adalah respon getaran bebas dari sistem SDOF tak teredam. Pertama-tama dengan mempertimbangkan kasus dari sebuah sistem yang menggantikan dari posisinya yang seimbang dengan jumlah u o dan dibebaskan. Kemudian u(0) = 0 , jadi ……………(3.18) u u cos t o
n
.
Gambar 3.1. Getaran bebas dari sistem SDOF tak teredam dengan
u(0) 0 .
DINAMIKA STRUKTUR
23
Dapat dilihat bahwa gerakan hasil adalah merupakan gerakan harmonik sederhana dengan amplitudo uo, periode alami dari sistem tak teredam (undamped natural period) yaitu ……………(3.19) 2 Tn (s) n dan frekuensi alami dari sistem tak teredam (undamped natural frequency) adalah 1 f n ……………(3.20) (Hz) n T n 2
Gambar 3.2. Respon getaran bebas secara umum dari sistem SDOF tak teredam.
Gambar diatas menunjukkan sebuah plot dari persamaan (3.17) apabila uo ataupun u adalah 0 (nol). Hal ini tetap merupakan gerakan harmonik sederhana o
dengan periode T n u(t) dapat diekspresikan dengan persamaan (3.17) ata u dengan persamaan ) U cos ……………(3.21) u(t ) U cos( nt n 1 n Getaran Bebas Pada Sistem SDOF Teredam (Damped) Persamaan (3.5) ditulis kembali disini :
3.3
2 nu 2u
u
0
……………(3.22)
n
Mengasumsi kembali sebuah solusi dari bentuk : uCe
st
……………(3.23)
sC e st
u
u s C e dan kita akan mendapatkan persamaan karakteristik : 2
s
2
st
2 s 2 n s 0 n
……………(3.24)
nilai s1 dan s2 adalah s1,2
n n 2 1
……………(3.25)
DINAMIKA STRUKTUR
24
Besarnya faktor "damping" ( ) dapat digunakan untuk membedakan 3 kasus, yaitu underdamped (0 < < 1), critically damped ( = 1), dan overdamped ( 1). Respon pada sistem SDOF teredam dengan beberapa variasi nilai redaman dapat dilihat pada gambar berikut.
Gambar 3.3. Respon dari sistem SDOF dengan redaman viskous dan variasi tingkat redaman.
Kasus Underdamped Underdamped (ζ<1) (redaman subkritis) Untuk < 1, lebih mudah bila menulis persamaan (3.25) dalam bentuk s1,2 n i d
……………(3.26)
dimana i 1 adalah unit imajiner dan d d adalah frekuensi alami "damped circular" yang diberikan oleh
d n 1 2
……………(3.27)
periode redaman (T (T d d) adalah 2 Td d
……………(3.28)
Dengan bantuan dari formula Euler, penyelesaian umum, u(t), u(t), dapat ditulis dalam bentuk t u(t ) e n ( A cos cos t A sin t ) ……………(3.29) 1
d
2
d
u0 dan u 0 digunakan untuk mengevaluasi A1 dan A2 , dengan ha sil : u o nuo t ……………(3.30)
n
u(t ) e
uo cos d t
d
A2 sin d t )
persamaan (3.30) dapat ditulis dalam bentuk: u (t ) Ue n cos( d t ) t
u n u0 2 U u0 0 d
……………(3.31)
2
……………(3.32)
jika harga ζ=20%, maka pada persamaan (3.27) (3.27)
d 0,98 n d n
……………(3.33)
DINAMIKA STRUKTUR
25
Substitusi persamaan (3.33) ke dalam persamaan (3.30), maka solusi gerak dapat digambarkan sebagai berikut
Gambar 3.4. Respon getaran bebas dari sistem redaman subkritis.
Gambar berikut menunjukkan perbandingan antara respon-respon dari sistem-sistem SDOF mempunyai level-level yang berbeda dalam redaman subkritis. Dalam tiap kasus, karena uo = 0, respon yang didapat adalah
Gambar 3.5. Pengaruh dari tingkat redaman pada getaran bebas.
Walaupun nilai dari mempunyai efek pada frekuensi, d , efek yang paling berat dari damping adalah pada angka pada saat gerakan menyusut, yaitu pada waktu e- dt . Efek ini akan dibahas lebih lanjut pada bagian 3.4, yang membahas ukuran dari damping.
DINAMIKA STRUKTUR
26
Kasus Critically-damped (ζ=1) (redaman kritis) Ketika ζ=1 maka persamaan (3.25) menjadi s1,2 n
……………(3.34)
Sehingga responnya menjadi: ……………(3.35)
nt
u(t ) (C 1 C2t )e Ketika kondisi adalah:
awal diperhitungkan, maka respon dari sistem redaman kritis
u(t ) [uo (u o un o)t ]e
nt
……………(3.36)
Gambar 3.6. Respon getaran bebas pada redaman kritis.
Kasus Overdamped (ζ>1) (redaman superkritis) Pada sistem redaman superkritis, koefisien redamannya lebih besar dari koefisien redaman kritis yaitu c ……………(3.37) 1 ccr
Eksperimen Penentuan dari Frekuensi Alami Dasar dan Faktor Damping dari sebuah sistem SDOF Metode eksperimen biasa dipakai untuk variabel dinamis pada suatu sistem (misal: frekuensi alami dan faktor redaman). Nilai konstanta pegas (k ) dan massa (m) dari sistem SDOF sederhana dapat diukur secara langsung. Namun nilai faktor redaman sering berubah sehingga perlu pengukuran yang lebih teliti. Bila faktor redaman diketahui, maka koefisien redaman bisa dihitung menggunakan persamaan faktor redaman. Frekuensi alami dari sistem SDOF tak teredam dapat ditentukan secara langsung melalui pengukuran statis. Contoh perhitungannya seperti pada contoh 3.2 berikut. 3.4
DINAMIKA STRUKTUR
27
Contoh 3.1 Tentukan frekuensi natural dari sistem pegas-massa dengan menggunakan pengukuran perpindahan secara statis.
Solusi:
k
Lo
k
fs=kust ust
w w Dari persamaan frekuensi alami struktur, diperoleh persamaan (1)
n2 k m Persamaan keseimbangan massa yang tergantung pada pegas adalah
F 0
(2)
atau
W f s 0
(3)
Dari persamaan gaya-perpindahan pada pegas
f s ku st Kombinasi persamaan 3 dan 4 2 g n u st
(4) (5)
apabila redaman dalam sistem kecil ( < 0.2), persamaan 3.32 menunjukkan bahwa nilai d kurang lebih sama dengan n. Sedangkan dari contoh 3.3 dapat diketahui bagaimana sebuah eksperimen getaran bebas dapat digunakan untuk menentukan frekuensi alami dari sebuah sistem SDOF. Contoh 3.2
Frekuensi natural dari balok kantilever dengan massa lumped (terpusat) bergerak dinamis. Massa bergerak dengan amplitudo A = 1 in kemudian dilepaskan. Gerakan yang terjadi ditunjukkan gambar di bawah yang
DINAMIKA STRUKTUR
28
mengindikasikan bahwa redaman pada struktur sangat kecil. Hitung frekuensi natural dalam radian per detik dan hertz. Berapa periodenya?
Solusi: Pada titik a, mass telah bergetar sepanjang 1,25 putaran. 1.25 putaran f n 3.125 Hz 0.4 s
n 2 f n (6.28)(3.125) 19.6 rad/s T n
1 f n
1 3.125
0.32 s
Terdapat dua metode yang hampir sama untuk menentukan faktor redaman ( ) menggunakan rekaman melemahnya getaran bebas dari sistem SDOF, yaitu metode pengurangan logaritmik dan metode setengah amplit udo dimana keduanya didasarkan pada persamaan 3.31.
Gambar 3.7. Rekaman melemahnya respon pada sistem teredam.
Pada metode setengah amplitudo, gerakan amplitudo (up) pada permulaan putaran dan amplitudo (uQ) pada akhir putaran diperkirakan besarnya. Pada akhir
DINAMIKA STRUKTUR
29
periode (satu putaran) nilai c os d t kembali lagi ke nilai pada a wal putaran. Dari persamaan 3.31, didapatkan rumus: u P e T uQ n d
……………(3.38)
Persamaan pengurangan logaritmik adalah: u P
ln u T n d
Q
……………(3.39)
Dimana Td adalah periode natural teredam yang dirumuskan sebagai berikut.
Td
2
d
2
2 n 1
……………(3.40)
Dari persamaan 3.38 dan 3.39 didapatkan
T nd
2 1
……………(3.41) 2
Untuk faktor redaman kecil ( < 0,2), persamaan persamaan pengurangan logaritmik mendekati nilai ……………(3.42)
2
Sehingga faktor redaman dapat diketahui juga menggunakan persamaan 1 U P ……………(3.43)
ln Q U
2
Prosedur yang sama juga dapat diterapkan pada metode setengah amplitudo, yang menghasilkan perhitungan lebih sederhana untuk faktor redaman. Metode setengah amplitudo didasarkan pada amplitudo dari kurva envelope. t uˆ(t ) Ue n
……………(3.44) Pada dua titik P dan R dimana: uˆ R
uˆ P
……………(3.45)
2
Titik-titik tersebut sejarak periode redaman N, dimana N bukan sebuah bilangan bulat. Selanjutnya, uˆ P
e NT 2
……………(3.46)
nd
uˆ R Dari persamaan 3.40 dan 3.46
2 N 1 2
ln(2)
……………(3.47)
DINAMIKA STRUKTUR
30
Gambar 3.8 menunjukkan hubungan antara dan N.
Gambar 3.8. Faktor redaman vs. jumlah putaran untuk mengurangi ampitudo sebesar 50%. 2
Untuk nilai faktor redaman yang kecil, << 1, persamaan 3.47 menjadi: 2 N ln(2) ……………(3.48) Sehingga, 0.11 ……………(3.49)
N
Persamaan 3.49 memberikan cara yang mudah untuk memperkirakan redaman dalam sebuah sistem yang teredam secara ringan ( < 0.1, misal N > 1) Contoh 3.3 Sebuah sistem bergetar terdiri dari berat W = 10 lb dan pegas dengan kekakuan K = 20 lb/in. Akibat redaman viskous (liat) sehingga terjadi amplitudo puncak 1,0 dan 0,85. Tentukan:
a) n
b) Pengurangan logaritma ln
y1 y2
c) d) c e) D Solusi: a) n
K m
K = 20 lb/in W 10 lb
m
g
n
20 10 386
27,78 rad sec
386 in/sec
2
DINAMIKA STRUKTUR
f
2
27,78
4,42
2
31
SPS
b) ln y1 y2
y1 = 1,00 y2 = 0,85
ln 1,0 0,165 0,85 c ) (untuk ξ kecil) 2 0,163 0,026 2 d )
c ccr
ccr
2
k m 2 10 20 386
c ccr
0,026 2 10 20 386 lb dt 0,037 in
Contoh 3.4
Gunakan metode setengah amplitudo untuk memperkirakan redaman dari sebuah sistem yang gerakannya terekam dalam gambar 3.10
Solusi:
Gambar sketsa dari kurva envelope (terdapat pada gambar) Ambil titik P pada salah satu puncak dan ukur u p
u p = 0,44 in
DINAMIKA STRUKTUR
32
Tempatkan titik R, dimana amplitudo dari kurvanya adalah u p / 2 = 0,22 in Perkirakan jumlah putaran antara P dan R
N = 2,25 putaran
Gunakan persamaan 3.49 untuk memperkirakan : 0.11 0.049 2.25 Level redaman dalam suatu sistem juga tercermin dalam konstanta waktu, , yang didefinisikan sebagai waktu yang diperlukan amplitudo untuk berkurang sejumlah faktor 1/e. Persamaan untuk menghitung konstanta waktu dapat menggunakan langkah yang sama dengan langkah yang dipakai untuk penurunan persamaan pada metode setengah amplitudo. Gunakan kurva envelope pada gambar 3.7 lagi, tentukan titik S dimana: u P uS
u P u P (1/ e)
e
……………(3.50)
Jadi, u P uS
U exp( nt P ) U exp[ n (t P )]
e
……………(3.51)
Atau, e n e
……………(3.52)
Dengan menggunakan logaitma pada kedua sisi, kita dapatkan:
n 1 Selanjutnya konstanta waktu, , didapat dengan persamaan: 1 T n n 2
……………(3.53)
……………(3.54)
Dari persamaan 1/e = 1 / 2,718 = 0,368. Maka, konstanta waktu, ,adalah waktu yang diperlukan amplitudo gerakan untuk berkurang sekitar 63 %. Getaran Bebas dari sebuah sistem SDOF dengan Coloumb Damping Struktur dengan redaman couloumb mempunyai persamaan gerakan diferensial linier sehingga menjadi lebih mudah diselesaikan untuk kasus respon getaran bebas ataupun respon akibat adanya gaya luar. Dalam praktek, redaman ini biasanya terjadi akibat hilangnya sambungan, gesekan antar komponen dan redaman dari material yang semuanya menyebabkan perilaku struktur menjadi nonlinier. Gambar 3.12 menunjukkan sebuah massa meluncur pada permukaan kasar yang menghasilkan gaya gesekan. 3.5
DINAMIKA STRUKTUR
33
Gambar 3.9. Sistem SDOF dengan redaman.
f D k N k mg
……………(3.55) adalah koefisien gesek kinetik atau koefisien gesekan luncur. Gaya
Dimana k
gesek selalu berlawanan arah dengan gerakan gaya u Newton II, kita peroleh:
. Menggunakan hukum
f s f D mu
……………(3.56)
Sedangkan fs = k . u dan Selanjutnya, mu mu
ku k mg , ku k
mg , Dengan u f D
u u
f D k mg sgn(u )
……………(3.57)
0 0
……………(3.58)
1 k g D k 2 n
……………(3.59)
Persamaan 3.58 dan 3.59 dapat digabungkan untuk mendapat: u
n 2u n Du
u
0
u 0
n 2u n Du
u
……………(3.60)
Gambar 3.10. Respon getaran bebas sistem dengan redaman Couloumb.
DINAMIKA STRUKTUR
34
Gerakan yang dihasilkan kemudian diplot dalam gambar 3.10. Yang perlu dicatat pada gambar 3.10 adalah bahwa sistem redaman couloumb berlaku seperti sistem SDOF tak teredam yang posisi seimbangnya berubah pada tiap akhir dari setengah putaran. Tampilan yang beda dari gambar 3.9 adalah amplitudo berkurang secara linier terhadap waktu, tidak secara eksponen seperti pada kasus redaman viskous.
DINAMIKA STRUKTUR
BAB IV RESPON SISTEM HARMONIS
SDOF
TERHADAP
35
GERAK
Pada bab ini, dibahas respon sistem SDOF baik yang tidak teredam maupun dengan redaman viskous terhadap gaya luar, dalam bentuk gerakan harmonis, yaitu struktur yang dibebani oleh gaya atau perpindahan yang besarnya dinyatakan oleh fungsi sinus atau cosinus dari waktu (p(t) = sin Ωt atau p(t) = cos Ωt). Contoh gerakan harmonis adalah gerakan mesin-mesin rotasi yang menghasilkan pengaruh harmonis akibat adanya eksentrisitas massa yang berotasi. Respon Sistem SDOF Tak Teredam Terhadap Gerakan Harmonis Respon total dari sistem linier terdiri dari superposisi respon akibat gerakan gaya luar dan respon dari gerakan natural. Sedangkan pada gerakan harmonis, gaya luarnya berupa respon steady-state. Berdasarkan gambar 4.1 yang menunjukkan Sistem SDOF tak teredam, diasumsikan bahwa sistem linier, amplitudo p 0 dan frekuensi gerakan Ω, persamaan gerakan adalah: ……………(4.1) mu ku p0 cos t 4.1
Nilai dari gaya luar (respon steady-state) berbentuk: u p U cos t
……………(4.2)
Untuk menentukan amplitudo, U, persamaan (4.2) disubstitusikan ke dalam persamaan (4.1): p0 ……………(4.3) U k m 2 Terlihat bahwa k m U 0
2
0 , maka defleksi statis: ……………(4.4)
0
k Kombinasi dari persamaan 4.3 dan 4.4 menghasilkan persamaan fungsi respon
frekuensi: p0 U k m 1 2 k U0 U m 2 1
m k
1 2 n
k H ()
1 1 r
2
,r 1
……………(4.5)
DINAMIKA STRUKTUR
36
dimana: r
……………(4.6)
2 n
dan U H () U0 r H(Ω)
……………(4.7)
= rasio rekuensi = fungsi respon frekuensi
Gambar 4.1. Gerak harmonis dari sistem SDOF tak teredam.
Fungsi respon frekuensi adalah fungsi yang memberikan penambahan atau pembesaran pada gerakan steady-state dalam bentuk nilai absolut dari fungsi respon frekuensi. Faktor pembesaran respon steady-state dirumuskan sebagai berikut: D s H ()
……………(4.8)
Dari gabungan persamaan (4.2) dn (4.5) memberikan persamaan respon steadystate sebagaiU berikut: 0 cos t , r 1 u ……………(4.9) p
1 r
2
Gambar 4.2. Faktor pembesaran untuk sistem SDOF tak teredam
(p(t) = po sin Ωt).
DINAMIKA STRUKTUR
37
Jika r < 1, maka responnya sefase / terdapat di dalam fase gerakan karena (1-r 2) bernilai positif. Jika r > 1, maka responnya 180° diluar fase / tidak sefase dengan gerakan, sehingga u p dapat ditulis: U 0 cos t u p
1 r
2
……………(4.10)
Persamaan respon total terdri dari solusi komplementer (uc) yang memenuhi persamaan homogen dan solusi partikulir (u p) yang memenuhi persamaan differensial nonhomogen. u u p uc
uc A1 cos nt A2 sin nt U 0 cos t A cos u t A sin t 2 1 r 1 n 2 n
……………(4.11)
Persamaan 4.9 dan 4.11 tidak dapat digunakan bila r = 1 atau
Ω = n yang
disebut dengan keadaan resonansi. Dari gambar 4.2 terlihat bahwa frekuensi gerakan yang berada dekat dengan resonansi, responnya menjadi sangat besar karena amplitudonya bernilai tak hingga. Oleh karena itu, memperhitungkan respon struktur terhadap gerakan harmonis sangat penting untuk menghindari kondisi resonansi dimana terjadi nilai amplitudo yang sangat besar. Namun biasanya bahan yang dipakai untuk struktur mempunyai limit kekuatan dan pada kondisi sebenarnya struktur akan runtuh jauh sebelum dicapainya amplitudo maksimum.
Contoh 4.1
Sistem pada gambar 4.1 mempunyai k = 40 lb/in, dan berat benda 38,6 lb.
Jika
uo
u
o
0 dan gaya luar P(t) = 10 cos (10t), tentukan pe rsamaan
gerakannya dan sketsa hasil pergerakannya.
DINAMIKA STRUKTUR
Solusi: Dari persamaan U 0 4.11, respon total adalah: u cos t A cos t A sin t
1 r
2
1
n
2
n
Selanjutnya, persamaan gerakan diturunkan untuk mendapatkan persamaan kecepatan: U 0 u sin t A sin t A cos t
1 r 2 Persamaan 3.4a
k n m
1
1
2
1
kg W
n
2
n
40(386)
2
n
n
20 rad/s
(38.6)
Persamaan 4.4: p 10 0.25 in. U 0 0 k 40 Persamaan 4.7: 10 r 0.5 n 20 Sehingga, U 0 2 1 r
0.25 1 (0.5)
2
0.25 1 0.25
0.33 in
Gunakan kondisi awal untuk menghitung A1 dan A2. U 0 A1 u(0) 0 2 1 r Maka: U 0 0.33 in A1 1 r 2
u (0) 0 A2 n Jadi, A2 = 0 u = 0,33[cos (10t) – cos (20t)] in
38
DINAMIKA STRUKTUR
39
Persamaan yang diperoleh kemudian digambarkan pada kurva di bawah ini.
Dari respon yang digambarkan pada contoh 4.1, maka:
Respon steady-state mempunyai frekuensi yang sama dengan gerakan dan berada di dalam fase gerakan karena r < 1.
Gerakan gaya dan gerakan natural saling memperkuat dan menghilangkan, menghasilkan fenomena tumbukan. Jadi, respon total bukan merupakan gerak harmonis sederhana.
Respon total maksimum (u = -0,66 in pada t = π/10 s) lebih besar pada pembesarannya daripada respon steady-state (u p = 0,33 in pada t = 0).
Total faktor pembesaran dinamis didefinisikan sebagai: D max u(t ) t U 0
……………(4.12)
Jika r = 1, maka asumsi yang digunakan pada persamaan (4.2) adalah : ……………(4.13) u Ct sin t , p
n
Kemudian, dengan mnsubstitusikan persamaan (4.13) ke persamaan (4.1), didapat: C Atau: u
p0
……………(4.14)
2m n
1 (U t ) sin t
p
2
0
n
……………(4.15)
n
Respon Sistem SDOF Redaman Viskous Terhadap Gerakan Harmonis Model analisis klasik dari sistem SDOF adalah model pegas-massa-dashpot (gambar 3.1). Ketika sistem tersebut dikenakan gerakan harmonis (p 0 cos Ωt), maka persamaan gerakannya menjadi: ……………(4.16) mu cu ku p0 cos 4.2
t
DINAMIKA STRUKTUR
Gambar 4.3. Respon u p(t) saat resonansi, Ω=
40
n .
Akibat adanya redaman pada persamaan (4.16), respon steady-state tidak akan berada dalam satu fase dengan respon steady-state: ……………(4.17) u U cos(t ) p
Dimana U adalah amplitudo steady-state dan α adalah sudut fase respon steadystate. Penentuan nilai U dan α dapat dilakukan dngan menggunakan putaran vektor. Kecepatan dan percepatan dirumuskan sebagai berikut: u ) u
p
U sin( t p
……………(4.18)
U cos( t ) 2
Vektor posisi dari gaya luar, perpindahan, kecepatan dan percepatan terlihat pada gambar 4.4.
Gambar 4.4. Vektor gaya, perpindahan, kecepatan dan percepatan.
Gambar 4.5. Poligon vektor gaya.
DINAMIKA STRUKTUR
41
Persamaan 4.17 dan 4.18 disubstitusikan ke persamaan 4.16, menghasilkan persamaan:
m2U cost cU sint kU cost p 0 cos t (4.19) Persamaan di atas diperoleh dari poligon vektor gaya dimana masingmasing variabel gaya menggambarkan gaya yang bekerja pada suatu massa. m U < kU yang berarti < 2
Gambar 4.5 menunjukkan kasus
n . Proyeksi
vektor dengan garis putus-putus pada gambar tersebut ditulis pada bagian kiri persamaan 4.19. Sedangkan proyeksi vektor dengan garis penuh ditulis pada bagian kanan persamaan 4.19. Dari gambar 4.5 juga bias diperoleh hubungan persamaan sebagai berikut:
(kU (kU m2U )2 (c ( cU )2 c tan k m 2 2
p0
…………(4.20a) …………(4.20b)
Sehingga nilai faktor pembesaran steady-state pembesaran steady-state dirumuskan dengan persamaan: D s
U U 0
tan α
1
1 r
2 r
2 2
2
1 2
…………(4.21a)
2ζ r
…………(4.21b) 2 1 r Kombinasi dari amplitude dan fase disebut respon frekuensi. Hubungan
antara rasio frekuensi dan faktor pembesaran steady-state digambarkan pada kurva gambar 4.6.
Gambar 4.6 . Kurva faktor pembesaran vs rasio frekuensi untuk berbagai nilai redaman.
DINAMIKA STRUKTUR
Contoh 4.2 Jika = 0,2 ditambahkan pada sistem contoh 4.1, dengan kondisi dan
perlakuan yang sama, tentukan persamaan gerakannya. gerakannya. Sketsa pergerakannya.
Solusi: Fungsi total respon didapat dari: t u U cos(t ) e n ( A cos t A sin t ) 1
d
2
Dimana:
U
U 0
(1 r )
2 2
(2 r )
2
1 2
n , ud dan r dapat ditemukan dari contoh 4.1 1
k 2 20 rad/s n m p 10 U 0 0.25 in 0 k 40 10 r 0.5 n 20 n (0.2)(20) 4 rad/s Oleh karenanya: U
0.25
1 0.5 2(0.2)(0.5) 2
2
2
2 r 2(0.2)(5) tan 2 0.267 2 1 r 1 (0.5) α = 0,26 rad
1 2
0.32 in
d
42
DINAMIKA STRUKTUR
43
Dari persamaan 3.31a,
d n 1 2 20 1 (0.2)2 19.6 rad/s Hasil diferensial total respon dari waktu: u
cos t A A sin t
t U sin( U sin( t ) e A A A n
2
d
1
n
d
1
d
2
n
Maka,
u(0) 0 0.32cos(0.26) A1 Sehingga:
A1 0.32cos( 0.26) 0.31 in A2 0.11 in Oleh karenanya,
u 0.32cos(10t 0.32cos(10 t 0.26)
[0.31cos(19.6 t ) 0.11sin(19.6t 0.11sin(19.6t )] )] in e4t [0.31cos(19.6t
d
DINAMIKA STRUKTUR
44
BAB V Respon Sistem SDOF Terhadap Bentuk Spesial Dari Eksitasi Pada berbagai situasi riil di lapangan, eksitasi dinamik yang terjadi tidaklah harmonik maupun periodik. Oleh karena itu, pada bab ini akan dibahas respon dinamik dari suatu sistem SDOF terhadap eksitasi. Respon Dari Sebuah Viscous-Damped System SDOF Untuk Sebuah Step Input yang Ideal Prototipe sistem SDOF yang ditunjukkan dalam gambar 3.1 merupakan bentuk subjek untuk sebuah step input yang ideal seperti ditunjukkan pada gambar 5.1. Dari gambar di bawah dapat dilihat bahwa sebuah gaya bekerja secara tiba tiba dari gaya nol (0) sampai dengan P o, selanjutnya nilainya konstan sebesar po. 5.1
P(t)
Po
t
Gambar 5.1. Sebuah stop input yang ideal.
Persamaan dari gerakan diberikan oleh persamaan 3.1, yaitu:
mu p0 ,
cu ku
t 0
……………(5.1)
Dianggap sistem berhenti pada t = 0 (kondisi awal), sehingga: u(0) u (0) 0
……………(5.2)
Penyelesaian dari persamaan 5.1 memuat sebuah bagian penyelesaian (a particular solution) dari persamaan 5.1, yang dapat ditulis menjadi: u p
0
……………(5.3)
k
Dan sebuah penyelesaian pelengkap (a complementary solution) diberikan (untuk
<1) oleh persamaan 3.32 sehingga: u(t )
p0
e t ( A cos t A sin t ) n
1
d
2
k Gunakan kondisi awal untuk evaluasi A 1 dan A2, kita peroleh: p0 n t u(t ) k
1 e
n
……………(5.4)
d
sin d t cos d t d
……………(5.5)
DINAMIKA STRUKTUR
45
Suatu cara yang berguna untuk menentukan respon dinamis suatu sistem adalah dengan memperhitungkan rasio respon atau suatu faktor beban dinamik, R(t ) , yang didefinisikan oleh: ku(t ) R(t ) pmax
……………(5.6)
Suatu faktor beban dinamik adalah rasio dari respon dinamis terhadap deformasi statis. Untuk step input ideal, R(t ) diberikan oleh: n t R(t ) 1 e
n
sin d t cos d t d
……………(5.7)
Suatu faktor beban dinamik yang sejenis diilustrasikan pada gambar 5.2. Pada rasio respon plot R(t ) 1 sesuai dengan posisi dari perpindahan statis. Karena beban diberikan secara langsung, terdapat overshoot, kemudian sistem akan tetap bertahan pada nilai statis yaitu 1 setelah melalui sejumlah gerakan bolak-balik yang teredam.
Gambar 5.2. Plot dari faktor beban dinamik untuk sebuah step input.
Untuk sebuah sistem tak teredam (undamped), persamaan 5.5 menjadi: p u(t ) 0 (1 cos t ) ……………(5.8) n k dan Rmax 2 Persamaan Respon dari sebuah Sistem Undamped SDOF pada Rectangular Pulse dan Pembebanan Ram Pada bagian rectangular pulse akan dibahas efek dari hilangnya beban setelah durasi t d . Gambar 5.3 menunjukkan sebuah input rectangular pulse dan rasio respon untuk sebuah sistem tak teredam untuk 2 kasus, yaitu T n T n t , t (a) d (b) d 2 2 dimana td adalah durasi dari rectangular pulse. 5.2
DINAMIKA STRUKTUR
46
Gambar 5.3. Respon dari sebuah input rectangular pulse, (a) rectangular pulse dan (b) rasio reaksi.
Dari gambar di atas terlihat bahwa ketika ketika
t d T n 2 , maka respon
maksimum terjadi sepanjang Force-vibration era, sedangkan jika t d T n 2 , maka respon maksimum terjadi di Residual-vibration era, dimana nilai maksimumnya dapat ditentukan pada tiap kasusnya. a. Kasus 1 : F orced-vibration era (0 t t d) Gambar 5.3 (b) memperlihatkan R(t) untuk sebuah pulse dengan durasi pembebanan sebesar
t d 5 4T n , dimana R(max) terjadi selama force-
vibration era. Untuk kasus ini, R(t) adalah sama untuk sebuah step ideal yang nilainya diperoleh dari persamaan (5.5) untuk sistem undamped, dimana: ……………(5.9) 0 t t R (t ) 1 cos t , 1
d
n
Nilai maksimumnya adalah ( R1 ) max
b.
T n R 2 1 2
……………(5.10)
Kasus 2 : Residual-vibration era (t d < t) Gambar 5.3(b) menunjukkan R(t) untuk sebuah pulse selama durasi
t d T n 8 . Rmaks terjadi selama residual-vibration era. Karena respon untuk t > td adalah vibrasi bebas dengan kondisi awal “initial condition”
R1 (t d )
dan R 1 (t d ) , maka pe rsamaan (3.17) dapat digunakan dalam bentuk R1(t d ) (t t ) sin (t t ) ………(5.11) R (t ) R (t ) cos 2
1 d
n
d
n
n
d
untuk t t d , dimana R1 (t d ) dan R 1 (t d ) diperoleh dari pe rsamaan 5.9.
DINAMIKA STRUKTUR
47
Gambar 5.4. Rotasi vektor yang merepresentasikan vibrasi bebas tak teredam.
Dari gambar 5.4 dapat dilihat bahwa amplitude U , dan sudut pada persamaan 3.21 ditentukan dengan 2 u U 2 u o2 o n dan tan
u
n
o
uo Persamaan amplitude U tersebut dapat digunakan untuk menentukan amplitude dan respon ini. R (t 2 21 ……………(5.12) 2 R R (t ) 1 d )
1
2 max
d
n
Persamaan di atas dapat ditulis sebagai berikut: t d ( R2 )max
2sin T n
……………(5.13)
Untuk memperhitungkan pengaruh dari durasi pembebanan pada respon maksimum, maka selanjutnya akan dibahas mengenai pengaruh dari peningkatan waktu pembebanan. Gambar 5.4 memperlihatkan hubungan antara beban ramp dengan peningkatan waktu t r yang diterapkan pada sistem undamped SDOF.
DINAMIKA STRUKTUR
48
P(t)
Po
t r
t
Gambar 5.5. Fungsi input ramp.
Persamaan gerakan pada kondisi awal adalah t P 0 t t r mu ku r o
t
P o u0 u 0 0
t t
.......................(5.14a)
.......................(5.14b)
r
……………(5.15)
untuk 0 ≤ t ≤ tr, solusi khususnya adalah u
t p0
……………(5.16)
t r k
p
Kemudian,
t p0 u
……………(5.17)
A1 cos nt A2 sin nt t r k
Dengan menggunakan kondisi awal dari persamaan 5.15, kita dapatkan
p0 t
1
……………(5.18)
sin nt t r n r
u k t
Untuk t t r , persamaan 5.14b dapat diselesaikan menjadi p0 1 u k 1 t
n
sin n t t r sin n t r
……………(5.19)
Gambar 5.6a memperlihatkan respon sebuah masukan dengan
t r T n serta
t r T n . Gambar 5.6b menggambarkan pengaruh dari kenaikan waktu pada respon maksimum.
DINAMIKA STRUKTUR
49
(a)
(b) Gambar 5.6 . Respon dari sistem SDOF tak teredam terhadap input ramp. (a) Respon terhadap input ramp. (b) Respon maksimum terhadap input ramp.
Dari gambar 5.6, dapat dilihat bahwa respon maksimum, Rmax 2 , terjadi pada step input ideal (misalkan untuk t r = 0). Untuk ramp dengan t r >> Tn akan terjadi sedikit overshoot dan sistem mengalami sedikit getaran bolak-balik atas kurva defleksi statis semu (pseudostatic deflection curve).
t p0 , t r k
……………(5.20)
u pseudostaitc
0 t t r
Respon Dari Sistem SDOF Tak Teredam untuk Impuls dengan Durasi Pendek, Unit Respon Impuls Pembebanan impuls adalah pembebanan yang berlangsung dalam selang waktu yang singkat. Impuls pada pembebanan ini didefinisikan sebagai perkalian dari gaya dan selang waktu bekerjanya gaya tersebut. Mengingat sistem SDOF tak teredam menyebabkan gaya dari durasi t d << Tn menghasilkan sebuah impuls 5.3
t d
I p(t )dt 0
……………(5.21)
DINAMIKA STRUKTUR
50
Gambar 5.7. Sistem SDOF tak teredam yang menerima impuls dengan durasi pendek.
Persamaan gerakan dan kondisi awal adalah mu
pt ku 0
u(0) u (0) 0
0 t t d
...................(5.22a)
t d t
...................(5.22b) ……………(5.23)
Dengan mengintegralkan persamaan 5.22a yang mempertimbangkan waktu dan menggabungkan kondisi awal, kita dapatkan mu (t d ) kuavg t d I
……………(5.24)
Dimana uavg adalah perpindahan rata rata(kecil) dalam interval waktu 0 < t ≤ td. Untuk td → 0, yaitu td ≤ Tn, bagian kedua pada persamaan 5.24 dapat diabaikan. ……………(5.25) mu (0 )
I Akan tetapi, sebuah impuls yang terdiri atas sebuah gaya yang besar dan bekerja pada waktu yang singkat berefek memberikan kecepatan awal yang besar, yaitu I u(0 ) …………(5.26a) m Akan tetapi abaikan dengan diganti perpindahan awal, yakni
u(0 ) 0
…………(5.26b) Hal ini bisa digunakan sebagai “awal” kondisi bagi permasalahan getaran bebas pada persamaan 5.22b. Dengan menggunakan persamaan 3.17, kita dapatkan respon dari impuls I ……………(5.27) u(t )
sin nt
mn
Fungsi respon dari impuls untuk sistem SDOF tak teredam didapatkan dari persamaan 5.27 dengan I = 1. Secara konvensi, unit fungsi respon dari impuls seringkali disebut h(t), sehingga 1 ……………(5.28) h(t ) sin nt mn
DINAMIKA STRUKTUR
51
Untuk sistem SDOF teredam-kental (viscous damped) dengan < 1, fungsi respon dari impuls dapat diperlihatkan sebagai
t sin d t e m d
u(t )
……………(5.29)
n
dan kesetaraan fungsi respon dari impuls didapat
h(t ) 1 e nt sin d t m d
……………(5.30)
Contoh 5.1
Asumsikan bahwa impuls I p(t )dt berasal dari gaya konstan po yang terjadi pada interval waktu 0 < t ≤ td hingga sistem SDOF tak teredam berada pada t = 0. Buktikan bahwa untuk td ≤ Tn, persamaan 5.11 dapat diringkas menjadi persamaan 5.27. Solusi: a.
Tentukan u t d dan u t d dari persamaan 5.8 p0 (1 cos t ) u(t ) d
u (t )
k n p0
d
k
...................(1)
n d
sin t
...................(2)
n d
Dengan nT n 2 dan t d << T n, nt d << 2 (1 cos ntd )
1 2
2 n dt
sin n t d n t d
td, sehingga Karena I = po p 1 0 ( t 2 I t d u(t ) ) d
u (t ) t
k
2 m
n d
(
d
b.
k n p 0
2
)
……………(3)
I
n d
m
Tinjau u(t) dari persamaan 5.11 dengan td → 0 I t d I
sin n (t t d ) cos n (t t d ) t 0 m n 2m
u(t ) lim d
...................(4)
DINAMIKA STRUKTUR
Jadi untuk td → 0 I sin nt u(t ) m n Seperti yang terlihat pada persamaan 5.27.
52
...................(5)
DINAMIKA STRUKTUR
53
BAB VI Respon System SDOF pada Eksitasi Dinamis
6.1
Metode Integral Duhamel Metode integral duhamel untuk menentukan respon dari sistem SDOF terhadap eksitasi dinamis secara umum dapat dikembangkan dari persamaan fungsi respon impuls, yang telah dijabarkan pada bab 5.3. Integral duhamel didasarkan pada prinsip superposisi, yang valid hanya untuk sistem linear. Gambar 6.1 menunjukan sistem SDOF tak teredam yang pada awalnya tak terganggu dimana kemudian menjadi pokok inputan P(t). respon dari sistem atas impuls dI= P(τ)dτ yang disebut du(t) dan didapatkan atas persamaan: dI du t sin n t ……………(6.1) mn
Gambar 6.1. Kenaikan respon pada sistem tak teredam.
Total respon pada waktu t adalah jumlah dari respons pada waktu t adalah jumlah dari respons yang berasal dari keseluruhan pias-pias impuls mulai dari waktu awal hingga waktu t, sehingga 1 t u(t ) ……………(6.2) 0 p( )sin n (t )d m n
Atau u(t )
t
0
p( )h(t )d
……………(6.3)
DINAMIKA STRUKTUR
54
dimana h(t- τ) didapat dari persamaan 5.28 untuk sistem tak teredam. Persamaan 6.3 adalah valid untuk sistem teredam jika persamaan 5.30 digunakan untuk mendapatkan h(tuntuk sistem teredam yang mulai dari waktu awal t 1τ). Sehingga, n (t ) u(t ) p( )e sin (t )d ……………(6.4a)
md 0 atau
1
u(t )
m
d
t
I 0
d
(i d n )(t )
p( )e
……………(6.4b)
d
Persamaan 6.2 dan 6.4 menjelaskan pernyataan integral duhamel untuk respons dari sistem SDOF baik yang tak teredam maupun yang teredam. Persamaan 6.3 seringkali menjelaskan sebagai integral konvulasi, dimana bentuk yang lebih umum adalah : x(t )
……………(6.5)
f 1 ( ) f 2 (t )d
Persamaan 6.2 atau 6.4 mungkin akan digunakan untuk menentukan respons dari sistem SDOF hingga eksitasi dinamis secara umum jika sistem ini dimulai dari waktu awal. Jika sistem memiliki kondisi awal tidak sama dengan nol, kemudian respons dari kondisi awal di tentukan dari persamaan 3.17 atau untuk ζ < 1, dari persamaan 3.33. itu, untuk sistem tak teredam u t Oleh 1 karena 0 u(t ) sin ……(6.6) t p( )sin (t )d u cos t
mn 0
0 n
n
n
n
dan untuk sistem di bawah teredam (under damped)
u(t ) 1 t p( )e (t ) sin d (t )d u0 e t cos d t m d 0 1 t ……………(6.7) u 0 n u0 n sin d t e
n
n
d
adalah tepat untuk menggunakan identitas trigonometri ketika mengevaluasi integrasi Duhamel ……………(6.8) sin (t ) sin t cos cos t sin
DINAMIKA STRUKTUR
55
Contoh 6.1 Gunakan integral Duhamel untuk menentukan respons dari sistem SDOF tak teredam dari sebuah beban “ledakan” yang ditentukan oleh pulsa triangular yang diperlihatkan pada gambar di bawah ini.
Menjelaskan pernyataan yang valid untuk t < t d dan untuk t > t d . Sistem tersebut dimulai pada waktu awal. Solusi: Gunakan persamaan 6.2 dengan t p(t ) p o 1 0 t t d t d t d t
p(t ) 0
a. Untuk 0 ≤ t ≤ td 1 u(t ) mn p0
t
o
p0 1 sin n (t )d t d t
k sin n t o 1 cos n d ( n ) t d t p0 k cos n t o 1 sin n d ( n ) t d Gunakan integral parsial, kita dapatkan 1 cos n d ( n ) sin n sin n d ( n ) n 1 sin n cos n n
Juga,
DINAMIKA STRUKTUR
56
1 sin n d ( n ) cos n n sin n Sehingga,
p0
t
1
1
u(t ) k sin nt sin nt t
sin nt cos nt t t d n d n d t 1 sin nt cos nt cos nt 1 cos nt t t d n d
Persamaan tersebut dapat disederhanakan menjadi t 1 sin nt R1 (t ) 1 cos nt t t d n d
b. Untuk td < t
1 t u(t ) mn o
d
Jadi,
p0
po 1 sin n (t )d t d
1
1
cos nt d sin nt t t n d n d 1 sin nt d cos nt 1 t n d
u(t ) k
sehingga,
1 sin nt (1 cos nt d ) cos nt ( nt d sin nt d ) R2 (t ) t n d
Contoh 6.2 Sebuah gedung yang ditujukan untuk mendapatkan gaya ledak dibuat model dengan sistem SDOF. Tentukan gaya ledak maksimum yang dapat ditahan bila perpindahan dibatasi sampai 5 mm dan apabila : (1) t d = 0.4 s, (2) t d = 0.04 s
DINAMIKA STRUKTUR
57
k = 9.0 GN/m, m = 10 Mg
Solusi :
a. Tentukan frekuensi dasar sistem
n f
n
k m
n 2
910
30 rad / s 1010 9
6
4.77 Hz
b. Menentukan rasio reaksi maksimum
Untuk kasus 1, t d = 0.4 s
1.91
f nt d 4.77 0.4
Dari reaksi spectrum, dapat dihitung R max = 1.75
Untuk kasus 2, t d = 0.04 s
0.191
f nt d 4.77 0.04
Dari reaksi spectrum, dapat dihitung R max = 0.58
c. Menentukan perpindahan statis untuk setiap kasus, kemudian tentukan nilai gaya ledak maksimum yang dapat ditahan ( po) umax = 5 mm u R p0 max
max
k
atau po
kumax Rmax
Untuk kasus 1
po 1
910
9
510
1.75 25.7 MN
Untuk kasus 2
3
DINAMIKA STRUKTUR
910
9
58
510 3
po 2
0.58 77.6 MN z(t)
u k
m
c
Gambar 6.2. Prototipe gerakan relatif sistem SDOF.
Sejak banyaknya penerapan dari reaksi spectra yang melibatkan gerakan relatif, ini berguna untuk menentukan jumlah reaksi yang sesuai untuk kasus ini. Gambar 6.2 menunjukkan prototipe dari gerakan relatif sistem SDOF. Seperti pada contoh 2.2, perpindahan relatif menjadi ……………(6.9) w = u – z Kemudian persamaan dari gerakan dapat ditulis ……………(6.10) mw cw kw m z definisi dari rasio reaksi seperti yang telah diberikan pada persamaan 5.6. Rasio reaksi berhubungan pada gerakan relatif hendaknya difenisikan sebagai
R t
w t m z k
2nwt
max
/
Dengan memasukkan
z
……………(6.11)
max
z
z max f a t dan berdasarkan persamaan 6.4 dapat
menunjukan R (t) pada integral Duhamel. Untuk sistem teredam
2 1
ω n
R(t)
ωd 0
n t τ
f a τ e
ζω
w 0
w 0 0
maka untuk
……………(6.12)
sin ωd t τ dτ
dan untuk sistem tak teredam 1
R(t) n f a sin n t d
……………(6.13)
0
Dari persamaan 6.11, perpindahan relatif maksimum diberikan oleh w 1 R z ……………(6.14) max
2n
max
max
DINAMIKA STRUKTUR
59
Jumlah kedua yang menjadi perhatian ialah akselerasi maksimum mutlak, u
max
. Persamaan 6.10 dapat ditulis ……………(6.15)
cw kw 0
mu
Untuk sistem yang tak teredam,
u
max
dapat dengan mudah ditentukan dari
persamaan 6.14 dan 6.15 u
max
n2 wmax
……………(6.16)
atau dari persamaan 6.14 dan 6.16,
u z
max
R max
(c 0)
……………(6.17)
max
Contoh 6.3 Sistem SDOF tak teredam seperti pada gambar 6.3 mengambarkan akselerasi dasar seperti gambar di bawah. Semua kondisi awal nol. Tentukan persamaan untuk wmax dan dan plot log-log dari wmax dengan f n. u
max
Gambar 6.3a. Pemodelan kendaraan yang bergerak di atas lintasan bump.
z dimana
z
f (t)
max a
DINAMIKA STRUKTUR
fa (t)
t 1 t d
0 t 2t d
0
2t d t
60
Solusi :
a. Menentukan reaksi untuk 0 ≤ t ≤ 2td Karena eksitasi untuk 0 ≤ t ≤ 2td memiliki format yang sama pada contoh 6.1, R1(t) akan sama dengan yang diberikan pada persamaan 7 dari contoh 6.1 t 1 sin nt R1 t 1 t cos nt t d n d
b. Menentukan reaksi untuk 2t d < t
R 1 2t d
sin n t 2t d n
R2 t R1 2t d cos n t 2t d
c. Menentukan waktu reaksi maksimum berdasarkan gaya – getaran, 2t d 1 R t 1 t sin t cos t 1 n d n n td
0≤t≤
dimana Ŕ1 = 0
nt m 2 tan1 tn d Pada kasus biasa
nt m 2 tan1 tn d
dengan ( nt m /2) berdasarkan pada kuadran
pertama Pada kasus terakhir
nt m 2tan 1 nt d p dimana p ialah integer terbesar untuk nt m < 2 nt d dan tan-1 ( n t d ) diambil pada kuadran pertama. Kemudian,
DINAMIKA STRUKTUR
nt m
R1 max 1
t
n
61
1
cos nt m sin nt m t d n d
d. Menentukan waktu residual-vibration maksimum yang terjadi, 2t d < t Dari persamaan 2, 2 R 2t 2 R1 2t d
R
1/ 2
1
d
2 max
n
dimana R1(2t d ) dan Ŕ1(2t d ) berdasarkan pada persamaan 1.
e. Menentukan pernyataan untuk wmax dan ϋmax : u wmax 1 dan
ax 2 d
z
t
R
2 n 2 max t d
z
max
R max
max
dimana Rmax adalah reaksi maximax yang diambil nilai terbesar antara (R1 )max dan (R2 )max.
f.
Plot wmax dengan fn menggunakan skala log – log Ini akan lebih mudah untuk menggambar reaksi w
ax
z
ax d
t
2
nondimensional
dengan frekuensi alami nondimensional f t . n d
DINAMIKA STRUKTUR
62
BAB VII Respons Spektrum Respons spektrum (response spectra) adalah plat respons maksimum (perpindahan, kecepatan, percepatan maksimum ataupun besaran yang diinginkan) dari fungsi beban tertentu untuk semua kemungkinan sistem berderajat kebebasan tunggal (SDOF). Perpindahan relatif |y-ys|max
m k
k
Frekuensi Natural f
y s (t)
(b) (a) Gambar 7.1 (a) Sistem SDOF yang dipengaruhi pergerakan tanah, (b) Bentuk spektrum respons.
Gambar 7.1(a) menunjukkan bangunan yang dibebani/dipengaruhi perpindahan tanah yang dinyatakan sebagai fungsi y s(t). Lengkung atau kurva spektrum respons pada gambar 7.1(b) memperlihatkan perpindahan relatif maksimum dari massa m terhadap perpindahan penyokong dari suatu sistem SDOF. Bentuk Respons Spektrum Bentuk grafik spektrum respons dapat dijelaskan dengan menggunakan sebuah osilator tak teredam yang dipengaruhi oleh setengah perioda gaya pengaruh sinusoidal pada gambar 7.2 di bawah ini. 7.1
DINAMIKA STRUKTUR
63
y F(t)
k
F(t)
Fo
m
td
(a) m y
ky
t
(b) F(t )
(c) Gambar 7.2 (a) Osilator sederhana tak teredam yang dipengaruhi beban F(t), (b) Fungsi beban
F t F o sin t 0 t t d
, dan (c) Free body diagram.
Persamaan gerak SDOF tak teredam :
m y F ( t )
ky
……………(7.1)
dengan,
0 F (t )
F sin t
o
untuk t t d
……………(7.2)
untuk 0 t t d
……………(7.3)
t d Solusi umum persamaan (7.1) merupakan superposisi antara solusi komplementer
yc dan solusi partikulir y p , y yc y p
yc Acos t Bsin t y p C sin t dimana
……………(7.4) ……………(7.5) ……………(7.6)
k m adalah frekuensi natural. Solusi khusus untuk selang waktu
0 t t d adalah sebagai berikut. Recall persamaan (7.1).
m y m y
ky F (t ) ky F (t )sin t
……………(7.1a)
Recall persamaan (7.6). y p C sin t y
C cos
t y
C 2 sin t
Susbstitusikan persamaan (7.6a) ke dalam persamaan (7.1a):
……………(7.6a)
DINAMIKA STRUKTUR
m(C sin t )
kC sin t o F Cm 2 Ck F o 2 C (k m ) F o
2
C
64
sin t
F o
……………(7.7)
(k m ) 2
Dengan mengkombinasikan persamaan (7.4) dan persamaan (7.7), maka respon untuk 0 ≤ t ≤ td adalah y yc y p y Acos t B sin t
F o sin t
k m F cos t y A sin t B cos t o 2 k m Untuk kondisi awal y(0) 0 , recall persamaan (7.8). 2
……………(7.8) ……………(7.8a)
0=A+0+0 A = 0 Untuk kondisi awal y (0) 0 , recall persamaan (7.8a). y
0 B
F o k m
2
B 2 k m Masukan harga A dan B kedalam persamaan (8). F sin t F o sin t o 2 y k m 2 k m F o sin t sin t y Fo
F o k y sin t sin t 2 1 k m
2
Untuk memudahkan analisa, persamaan tersebut ditulis F o 2 , , y st T k t d Maka persamaan (7.9) menjadi
……………(7.9)
DINAMIKA STRUKTUR
y y st
65
t T sin 2 t sin t 2t untuk 0 t t d ………(7.10a) T d d T 1 2t d 1
2
Solusi pada persamaan (7.10a) untuk t > td adalah y t t T t cos d sin 2 t d d 2
T 1 2t d
y st
T
T
…………(7.10b)
2T
y y st adalah
Pada persamaan (7.10) terlhat bahwa respon dalam besaran
fungsi dari rasio waktu pulsa (pulse duration) dengan periode natural dari sistem (t d /T ) dan dari waktu yang dinyatakan dengan t/T . Jadi dari harga tertentu parameter t d /T akan diperoleh respon maksimum pada persamaan (7.10). Gambar 7.3 menunjukan respons spektrum dari persamaan (7.10).
y y st T 2 m
k
t d T Gambar 7.3 Respons spektrum untuk setengah gaya sinusoidal dengan selang waktu t d.
Respons Spektrum pada Pondasi yang Bergerak Analisa sistem yang dipengaruhi oleh beban pada perletakan/pondasi merupakan suatu masalah penting yang ada pada dinamika struktur. 7.2
s
m
k
c
k
k
m c
(a)
I m. y
f S I
f D (b)
f S k ( y y s ) f D c( y
y
s
)
DINAMIKA STRUKTUR
66
Gambar 7.4 Osilator sederhana teredam yang dipengaruhi oleh penyokongnya, (b) diagram free body.
Percepatan getaran yang menyebabkan pergerakan pondasi bisa digambarkan sebagai berikut: y s (t )
t
Gambar 7.5 Fungsi percepatan yang memmpengaruhi penyokong dari osilator pada Gambar 7.4.
Persamaan gerak sistem pada diagram free body dapat ditulis
m y 0
c( y y s ) k ( y y s )
……………(7.11)
dengan,
k
m
,
c ccr
,
ccr 2
k .m
maka persamaan (7.11) menjadi: y (t )
2 2 y 2 y s y
(t ) 2 s y
……………(7.12)
yang merupakan persamaan differensial gerak dari osilator teredam dal am besaran gerak absolut. Jika dirumuskan perpindahan relatif u yang didefinisikan sebagai ……………(7.13) u y y s maka persamaan (7.12) dapat ditulis
2 u 2u s y
u
(t )
……………(7.14)
solusi persamaant (7.14) diperoleh dengan menggunakan integral Duhamel sebagai 1 ……………(7.15) u(t )
y s (
(t )
sin (t )d
)e 0
Besaran- Besaran Respons Spektrum Spektrum perpindahan S D adalah perpindahan relatif maksimum yang linear dengan spektrum percepatan Sa yaitu percepatan absolut maksimum. 2 ……………(7.16) S S 7.3
a
S V S D
D
S a
……………(7.17)
DINAMIKA STRUKTUR
67
Suatu contoh respons spektrum perpindahan untuk sistem SDOF yang dipengaruhi oleh gerak penyokong terlihat pada gambar 7.6, yang merupakan respons dari gerakan hasil rekaman percepatan tanah pada gempa di El Centro 1940. Plot dari rekaman percepatan gempa ini terlihat pada Gambar 7.7.
Gambar 7.6 Respons spektrum perpindahan untuk sistem elastis yang dipengaruhi pergerakan tanah akibat gempa di El Centro 1940.
Gambar 7.7 Rekaman percepatan tanah untuk gempa El Centro, California 18 Mei 1940 konponen utara-selatan.
Pada Gambar 7.8, bentuk data yang sama digunakan untuk mendapatkan respons spektrum perpindahan pada Gambar 7.6, yang diplot dalam besaran spektrum kecepatan untuk beberapa harga koefisien redaman, dengan perbedaan absis dan ordinat dalam skala logaritmis.
DINAMIKA STRUKTUR
68
Gambar 7.8 Respons spektrum sistem elastis untuk gempa El Centro 1940.
Untuk menyatakan bentuk dari diagram tiga besaran pada Gambar 7.8, persamaan (7.17) ditulis dalam besaran frekuensi natural f dalam siklus per detik (spd) dan mengambil harga logaritmanya, akan didapat S V S D 2 f S D ……………(7.18) log S V log f log( 2 S D ) S S a S a V 2 f log S V
log f og
l S a
……………(7.19)
2 Respons Spektrum untuk Perencanaan Elastis Gempa bumi terdiri dari suatu seri gerakan tanah yang bersifat acak (random). Biasanya komponen utara-selatan, timur-barat dan komponen vertikal dari percepatan tanah yang diukur. Respons spektrum rencana yang merupakan gabungan spektrum beberapa gempa bumi yang dinyatakan oleh suatu bentuk spektrum “rata-rata” digunakan dalam perencanaan struktur tahan gempa, karena saat ini tidak ada metode yang dapat menduga bentuk gerakan gempa pada suatu lokasi yang akan terjadi. Respons spektru rencana ini diperlihatkan pada Gambar 8.9. 7.4
DINAMIKA STRUKTUR
69
Gambar 7.9 Respons spektrum rencana yang dinormalisasikan untuk 1.0g .
Contoh 7.1 Struktur dengan model sistem berderajat kebebasan tunggal mempunyai perioda alami Tn= 1 detik. Metoda spektrum respons untuk menentukan percepatan absolut maksimum, perpindahan relatif maksimum dan kecepatan (pseudovelocity) relatif maksimum untuk: a) gerakan pondasi yang sama dengan gempa El Centro 1940 b) gempa rencana dengan percepatan tanah maksimum sebesar 0.32g . Dengan anggapan redaman sebesar 10% redaman kritis.
Solusi : a. Dari spectrum respons pada gambar 8.8 dengan f=1/T=1.0 spd dan ξ=0.10, maka S D = 3.3 in S V = 18.5 in/dt S a = 0.30 g b. Dari spectrum dasar rencana pada gambar 8.9, dengan f=1/T=1.0 spd dan ξ=0.10, untuk percepatan tanah maksimum 0.32g, maka S D = 9.5 x 0.32 = 3.04 in,
DINAMIKA STRUKTUR
S V S a
70
= 60 x 0.32 = 19.2 in/det, = 0.95 x 0.32g = 0.304g.
BAB VIII SISTEM BERDERAJAT KEBEBASAN BANYAK (MDOF) Sistem MDOF Sederhana Persamaan gerak untuk sistem MDOF sederhana, dapat diidealisasikan pada struktur portal tingkat dua dengan gaya luar p1(t) dan p2(t) (gambar 8.1). 8.1
Gambar 8.1. (a) Struktur portal tingkat dua (b) gaya yang bekerja pada kedua massa
Pada idealisasi tersebut balok dan lantai adalah kaku. Massa yang terdistribusi pada seluruh gedung. akan diidealisasikan terpusat pada bidang lantai. Asumsi tersebut umumnya sesuai untuk bangunan bertingkat. Pada gambar 8.1a diatas, portal tingkat dua dengan massa terpusat pada setiap lantai memiliki dua DOF : perpindahan lateral u1 dan u2 pada kedua lantai dalam arah x. Gaya-gaya yang bekerja untuk setiap massa lantai m j dapat dilihat pada gambar 8.1b., termasuk gaya luar p j(t), gaya elastic f Sj dan gaya redaman f Dj. Gaya elastis dan redaman menunjukan arah yang berlawan, karena kedua gaya tersebut adalah gaya dalam yang menahan gerakan. Hukum Newton Kedua pada Sistem MDOF Persamaan gerak dari hukum Newton kedua yang diberikan untuk setiap massa adalah 8.2
m j u (t )
j
f D j f Sj p j
……………(8.1)
Persamaan diatas terdiri dari j=1 dan j=2 sehingga dapat ditulis dalam bentuk matrik ; m1 0 u 1 f D1 f S 1 p1 ……………(8.2) (t ) 0 m f p (t ) f u 2 2 D2 S 2 2 atau dapat ditulis ;
mu
f D f S p(t)
……………(8.3)
DINAMIKA STRUKTUR
71
dimana
u1 u u2
m
m1 0
m2 0
f D1 f D2
f D
fS
f f S 1 S2
p1 p p1
Gaya elastis f S berhubungan dengan perpindahan yang terjadi pada setiap lantai u. Oleh karena itu, kekakuan lateral k j untuk setiap lantai ke- j memberikan hubungan geser pada lantai V j terhadap deformasi lantai, Δ j = u j-u j-1. V j k j j
……………(8.4)
Kekakuan pada setiap tingkat atau lantai adalah jumlah kekakuan lateral dari semua kolom di lantai tersebut. Tingkat atau lantai dengan tinggi h dan kolom dengan modulus E dan momen inersia I c maka kekakuan lantai tersebut adalah k j
12 EI h
……………(8.5)
c
3
kolom
Pada gambar 8.1, kita dapat menghubungkan gaya elastis f S1 dan f S2 terhadap u1 dan u2 .Gaya f S1 pada lantai pertama tersusun atas f a dari tingkat atas dan f b dari S 1
tingkat bawah. Oleh karena itu
S 1
f b f a
f S 1
S 1
S 1
f S 1 k 1u1 k 2 (u1 u2 )
……………(8.6a)
Gaya f S2 pada lantai kedua adalah
f S 2 k 2 (u2 u1 )
……………(8.6b)
Persamaan (8.6a) dan (8.6b) dalam bentuk matrik adalah f S 1 k 1 k 2 k 2 u1 atau f ku f k S k u
S2
2
2
……………(8.7)
2
Dengan cara yang sama pada persamaan (8.6), dapat diperoleh
f D1
c1u 1 c2 (u 1 u
2)
dan dalam bentuk matrik adalah f D1 c1 c2
f c D2 2
……………(8.8)
f D2 c2 (u 2 u 1) c2 u 1 atau f
c
2 u2
cu
……………(8.9)
D
Dengan mensubstitusikan persamaan (8.7) dan persamaan (8.9) kedalam persamaan (8.3), maka diperoleh ……………(8.10) mu cu ku p(t) Prinsip D’Alembert’s pada Sistem MDOF Berdasarkan prinsip D’ Alembert’s pada bab sebelumnya, adanya gaya inersia pada kesimbangan dinamis pada sebuah struktur. Untuk dua massa dalam 8.3
DINAMIKA STRUKTUR
72
sistem pada gambar 8.1a, free body diagram dan gaya inersianyanya dapat dilihat pada gambar 8.2, dimana untuk setiap gaya inersia adalah perkalian massa dengan percepatannya.
Gambar 8.2. Free Body Diagram 8.4
Sistem Massa – Pegas – Redaman
Gambar 8.3. (a) Sistem berderajat dua; (b) free body diagram
Persamaan gerak untuk sistem diatas telah ditunjukan oleh persamaan (8.10), sehingga; m1 0 u 1 c1 c2 u 1 k 1 k 2 u1 p1 (t ) ……(8.11) k 2 0 mc2 c p (t ) k c k 2 2 2 u2 2 u2 2 2 u2
Contoh 8.1 Buat persamaan gerak untuk portal dua tingkat dibawah ini.
DINAMIKA STRUKTUR
73
Solusi:
m1 2m k1 2
12( 2 EI c ) h
3
m2 m
48 EI c h
k2 2
3
12(EI c ) h
3
24 EI c h3
Substitusikan ke persamaan (8.2) dan (8.7), sehingga diperoleh matrik massa dan matrik kekakuan: 1 2 0 24 EI c 3 k mm
0
h3
1
Jadi persamaan gerak untuk sistem ini adalah 2 0 u 1 24 EI c m 3
0
1u2
3
h
1
1
1
1u1 p1 (t ) 1 u2 p2 (t ) I2
I1
Contoh 8.2 Buat persamaan gerak untuk portal tiga tingkat (bangunan berlantai tiga) dibawah ini. u3 (t)
(t)
3
3
m3 3 k 3(u u) 32
u2
k 3 (u 3 u2 )
(t)
2(t)
2
m2 2 k 2 (u2 u 1 )
u1
k 2 (u2 u 1 )
(t)
1
(t)
1
m1 1 k 1u1
DINAMIKA STRUKTUR
74
Solusi:
u1 (t)
u2 p3(t)
(t)
1
2
(t)
1
m1 u k1u1
u1 p2(t) k 2(u u )1 2
u3
u2 m2 u
u3 (t)
3
2
m3 3
u) k 3(u 32
Persamaan-persamaan gerak dari masing-masing free body diagram pada se tiap massa, m1u
k 1u1 k 2 (u2 u1 ) p1 (t ) 0 2 k 2 (u 2 u1 ) k 3 (u3 u 2 ) p 2 (t ) 0 3 k 3 (u3 u 2 ) p3 (t ) 0
1
m2u m3u
Sehingga persamaan gerak dalam bentuk matrik dari sistem ini adalah
m1
0
0
m
0
0
2
(k 1 k 2 )
0
u
k 2
0 u1
1
k u2 02 m3 u 3 0
(k k ) 2
k 3
3
k p (t ) 2 3 u2 k u p (t ) 3
p1 (t )
3 3
Koefisien Kekakuan Elemen-elemen dari mtriks kekakuan pada persamaan (8.7) disebut koefisien kekakuan. Dimana pada umumnya koefisien kekakuan k ij didefinisikan sebagai gaya pada koordinat i bila satu perpindahan diberikan pada titik j. Sebagai contoh, koefisien pada baris kesatu dan kolom kesatu dari persamaan (8.7) adalah k 11=k 1+k 2 menyatakan gaya pada lantai kesatu akibat satu perpindahan yang diberikan pada lantai tersebut. 8.5
DINAMIKA STRUKTUR
75
Contoh 8.3 Buat persamaan gerak pada contoh soal 8.1 dengan menggunakan koefisien kekakuan.
Solusi:
Matrik kekakuan Pertama, kita tentukan matriks kekakuan dengan menentukan nilai u1 = 1 dan u2 = 0. Koefisien kekakuan adalah k i1 . Diperlukan gaya pada bagian atas dan bawah untuk setiap lantai atau tingkat untuk menahan perubahan bentuk pada struktur, yang digambarkan oleh kekakuan k 1 dan k 2. Dari contoh 8.1, diperoleh 48 EI c 24 EI c k k 1 2 h3 h3 Dua gaya pada gambar (a) dan (b) diatas, k 11 k 1
k 2
72 EI c h3
k 21 k 2
24 EI c h3
Kedua, kita tentukan matriks kekakuan dengan menentukan nilai u1 = 0 dan u2 = 1. Koefisien kekakuan adalah k i2 . Diperlukan gaya untuk menahan perubahan bentuk yang digambarkan oleh gambar (d). Dua gaya pada gambar (c) dan (d) diatas, k 12
k 2
24 EI c h3
k 22
k 2
24 EI c h3
Dengan koefisien kekakuan yang telah ditentukan, maka matriks kekakuannya adalah k
k 11 k 21
k 12
k 22
24 EI c 3 h3
1
1 1
DINAMIKA STRUKTUR
Sedangkan matrik massa,
2 m m 0 Persamaan gerak adalah 2 0 u m 3
0
1
24 EI c 1u1 p1 (t ) 3 h 1 1 u2 p2 (t ) 0 1u2 1
76
DINAMIKA STRUKTUR
77
BAB IX GETARAN BEBAS UNTUK SISTEM MDOF Sistem MDOF Tak Teredam Persamaan gerak MDOF tak teredam dengan p(t)=0, ……………(9.1) mu ku 0 Terdapat dua kemungkinan gerak harmonis dari struktur sedemikian rupa, dimana semua massa bergerak dengan fasa tertentu pada frekuensi ω1 dan ω2. Setiap karakteristik perubahan bentuk disebut normal atau pola natural dari getaran. Sering disebut dengan pola pertama (first mode) atau pola dasar (fundamental mode) untuk menyatakan pola yang sesuai dengan frekuensi terendah. Pola yang lain disebut pola harmonis atau pola harmonis yang lebih tinggi. Gambar 9.1 dan 9.2 menunjukan getaran bebas pada portal dua tingkat. Kekakuan dan massa yang terpusat dapat dilihat pada gambar 9.1a dan mode getar atau pola getar ditunjukan oleh gambar 9.1b dan 9.2b.Hasil gerak u j pada sistem digambarkan oleh gambar 9.1d dan 9.2d. 9.1
Gambar 9.1. Getaran bebas pada sistem tak teredam dengan pola natural pertama dari getaran (a) Struktur portal tingkat dua; (b) perubahan bentuk struktur pada waktu a,b,c; (c) modal coordinate qn(t) (d) perpindahan
DINAMIKA STRUKTUR
78
Gambar 9.2. Getaran bebas pada sistem tak teredam dengan pola natural kedua dari getaran (a) Struktur portal tingkat dua; (b) perubahan bentuk struktur pada waktu a,b,c; (c) koordinat modal qn(t) (d) perpindahan
Perioda alami dari getaran T n pada sistem MDOF adalah waktu yang diperlukan untuk satu siklus dari gerak harmonis sederhana dalam satu pola natural. Hubungan terhadap frekuensi natural sudut dari getaran adalah ωn dan frekuensi natural adalah f n, 1 2 ……………(9.2) Tn f n n n Gambar 9.1dan 9.2 menunjukan perioda alami T n dan frekuensi natural sudut dari ωn (n=1,2) dari getaran bangunan 2 tingkat dengan pola natural ( )T . n
1n
2n
Frekuensi natural sudut yang lebih kecil diberi notasi ω1 sedangkan yang lebih besar dinotasikan ω2. Sedangkan untuk perioda alami yang lebih panjang dinotasikan T 1 dan yang lebih pendek adalah T 2. Frekuensi Natural dan Pola Normal Getaran bebas pada sistem tak teredam , yang secara grafis telah ditunjukan oleh gambar 9.1 dan 9.3 untuk sistem dua DOF, dapat diuraikan secara matematis adalah 9.2
u(t ) qn (t ) n
……………(9.3)
Variasi waktu pada perpindahan yang terjadi dapat diuraikan dengan fungsi sederhana harmonis ……………(9.4) q (t ) A cos t B sin t n
n
n
n
n
Substitusikan persamaan (9.4) ke (9.3)
u(t ) n ( An cos nt Bn sin nt )
……………(9.5)
DINAMIKA STRUKTUR
79
dimana ωn dan n tidak diketahui. Substitusikan persamaan (9.5) kedalam persamaan (9.1), sehingga didapatkan 2 [ m k ] q (t ) 0 ……………(9.6) n
n
n
n
Persamaan (9.6) dapat diselesaikan dengan satu dari dua cara. Salah satunya, qn(t)=0 yang memberikan nilai u(t)=0 dan tidak adanya gerak pada sistem atau frekuensi natural sudut ωn dan pola perubahan n persamaan aljabar berikut
yang harus memenuhi
2m k 0 n
n
……………(9.7)
n
dimana persamaan ini menunjukan kondisi maksimal. Matriks kekakuan k dan matriks massa m adalah diketahui, masalahnya adalah menentukan nilai skalar 2 dari dan vector dari . Persamaan (9.7) dapat ditulis kembali menjadi n
n
[k
n2m] n 0
……………(9.8)
Persamaan (9.8) adalah masalah matematis yang penting, yang dikenal sebagai “eigenproblem”, yang mempunyai soulusi nontrivial 2 ……………(9.9) det[k nm] 0 Pada umumnya jawaban persamaan (9.9) mempunyai bentuk persamaan polynomial derajat n dalam besaran ω2 yang harus mempunyai n buah harga ω2, yang memenuhi persamaan tersebut atau dikenal sebagai persamaan karakteristik. Sehingga kita dapat menyelesaikan persamaan (9.8). 9.3
Sifat Ortogonalitas dari Pola Normal Kita tinjau kembali persamaan (9.7) , k n
n2m n
……………(9.10)
untuk sistem berderajat kebebsan dua (lihat persamaan 8.7), sehingga k 1 k 2 k 2 1 m1 2 0 1
k 2 2 0 k 2 2 (k 1 k 2 ) 1 k 2 2 m1 2 k k m 2 2 1
2
2
2
2 m2 2
……………(9.11)
2
Digunakan teori Betti yang menyatakan bahwa, pada sebuah struktur yang dibebani oleh dua sistem pembebanan dimana terjadi dua jenis perpindahan, maka kerja yang dilakukan sistem pembebanan pertama sepanjang perpindahan akibat sistem pembebanan kedua, akan sama dengan kerja akibat sistem pembebanan kedua yang bergerak sepanjang perpindahan akibat sistem pembebanan pertama.
DINAMIKA STRUKTUR
(a)
80
(a)
Gambar 9.2. Model sejumlah massa dan perpindahan pada struktur bertingkat dua (a) Sistem I; (b) Sistem II
Kedua sistem pembebanan dan perpindahan yang akan ditinjau adalah Sistem I 2 2 Gaya-gaya m , m dan perpindahannya , 1 1 11
Sistem II 2 Gaya-gaya m
2 1 21
11
21
2 , m dan perpindahannya ,
1 2 12
2 2 22
12
22
Penggunaan teori Betti untuk kedua sistem ini menghasilkan, m 2 m 2 m 2 m 2 1 1 11 12
atau
2 1 21 22
1 2 12 11
2 2 22 21
( 2 2 )(m m 1
2
1 11 12
) 0
……………(9.12)
2 21 22
Jika ω1≠ ω2, maka persamaan 9.12 didapat ……………(9.13)
m1 11 12 m2 21 22 0
Persamaan diatas disebut hubungan ortogonalitas antara pola dasar dari sistem berderajad – kebebasan dua. Untuk sebuah sistem berderajat kebebasan n dimana matriks massa adalah matrik diagonal maka kondisi ortogonalitas antara pola n dan r dapat dinyatakan sebagai T ……………(9.13) nT m r 0 n k r 0 Atau kondisi orthogonal untuk pola n dan r dapat diperoleh melalui penjabaran sebagai berikut 2 [ K ]{ } [ M ]{ } ……………(9.14) n
n 2
n
r
r
r
[ K ]{ } [ M ]{ }
……………(9.15)
Apa bila persamaan (9.14) ditranspose, [ K ]{ } 2 [ M ]{ } sedangkan [ K ] T
n
n
n
dan [M] berbentuk diagonal untuk struktur biasa, oleh sebab itu
[ A] [ A]T matriks simetris Sehingga persamaan (9.14) menjadi T { } [ K ] n
2{ }T [ M] n
n
……………(9.16)
DINAMIKA STRUKTUR
81
Untuk menyelesaikan kedua persamaan matriks diatas maka kalikan persamaan (9.16) dengan { r }dan kalikan persamaan (9.15) dengan { n}T , sehingga T 2 T { } [ K ]{ } { } [ M ]{ } n
r
n
n
r
{ } [ K ]{ } { } [ M ]{ } T
n 2
2
n
r
r
r
2
T
n
r
( ){ } [ M ]{ } 0 T
n
……………(9.17)
r
Jika ωn≠ ωr , maka didapatkan nilai yang sama dengan persamaan sebelumnya (9.13) { n }T [ M ]{ r } 0
juga
{ n }T [ K ]{ r } 0 …………(9.18)
Contoh 9.1 Tentukan frekuensi alami dan pola pada sistem yang ditunjukan gambar dibawah ini.(lihat contoh 8.1)
Solusi: Dari contoh 8.1 diperoleh nilai matriks massa dan matriks kekakuan sebagai berikut 2 0 1 3k k 2m 0 24 EI c 3 24 EI c mm k k 0 3 3 h 1 1 k 1 0 m k h Nilai frekuensi alami ωn dapat diselesaikan dari persamaan (9.9) det[k n m] 0 3k 2m 2 k 2 n det k k m 0 2
n
2m2 5km 4
Akar-akar persamaan diatas adalah k k 2 1 1 2m 2m
2
2k 2 0 2 2
2k
m
2
Pola natural untuk sistem I diperoleh dengan mensubstitusikan persamaan (9.8), sehingga
2k m
ωn = ω1 pada
DINAMIKA STRUKTUR
3k 2m 2 1 k
2k k
k 11
k
k m
k 11 0 0.5k 21
82
0 1 21 2
2
1
2m
Biasanya pola natural atau normal ditentukan dengan menentukan satu satuan harga untuk salah satu pola, jadi ditentukan untuk 21 =1 dan diperoleh nilai 11 =0.5 Pola natural untuk sistem II diperoleh dengan mensubstitusikan persamaan (9.8), sehingga 2 k 12 2k 3k 2m 2 2 2 0 k k m 2 m 2 k k 12 0 k k 22
22
ωn = ω2 pada
Biasanya pola natural atau normal ditentukan dengan menentukan satu satuan harga untuk salah satu pola, jadi ditentukan untuk 22 =1 dan diperoleh nilai 12 =-1 Jadi,
1
11 21
1 2 1
1 1 12 1
kontrol kondisi orthogonal
nT k r 0 1T k 2 0
31
k 1/ 2 1
1 1 1 0 1
T 1 m 2 0
02
m 1/ 2 1
01 1
1 0
Tn m r 0
22
DINAMIKA STRUKTUR
83
Solusi Persamaan Getaran Bebas pada Sistem Tak teredam Solusi umum persamaan gerak, diberikan oleh persamaan (9.5). Sehingga untuk nilai n=1,2,3…,n maka persamaan (9.5) dapat ditulis menjadi 9.4
n
u(t )
( A cos t B sin t ) n
n1 n
n
n
n
……………(9.19)
n
( A sin t B
u (t )
n
n
n
n
……………(9.20)
n
cos n t )
n1
Pada saat t=0 maka persamaan tersebut dapat ditulis
……………(9.21)
n
q (0)
u(0)
n
n1
n
……………(9.22)
n
u (t )
q n
n1
n
(0)
dan saat t=0, persamaan (9.14) dan (9.15) memberikan
……………(9.23)
n
A
u(0)
n1
n
n
An qn (0)
n
u (t )
Jadi, u(t )
n
n
Bn
n1
(0) cos t
(q
n n Bn
n
q n (0)
n
n1
q n (0)
n sin t ) n
n
……………(9.24)
n
q (t )
nn
………(9.25)
n1
Contoh 9.2 Tentukan respon getaran bebas pada portal dua tingkat untuk contoh 9.1. Dengan
nilai q1 (0) 1, q2 (0) 1dan q 1 (0) 0, q 2 (0) 0 Solusi: Dari persamaan (9.20) didapatkan q1 (t ) 1cos 1t
q2 (t ) 1cos 2t
Dengan mensubstitusikan nilai n dari hasil perhitungan contoh 9.1 dan nilai qn (t ) diatas ke persamaan (9.20) u1 (t )
1/ 2 1 cos 1 t cos 2t 1 u2 (t ) 1
DINAMIKA STRUKTUR
9.5
84
Respon Pada Gedung Akibat Gempa Secara umum persamaan geraknya adalah [ M ]{u [ M ]{u }
u g } [C ]{u } [ K ]{u} 0 } [C ]{u } [ K ]{u} [ M ]{u
………(9.26)
g
………(9.27)
{u} { i }{ Ai }
Kalikan persamaan (9.26) dengan { i}T dan susbtitusikan persamaan (9.27) [ M ]{u u g } [C ]{u } [ K ]{u} 0 T T T { } [ M ]{ }{ A } { } [C ]{ }{ A } { } [ K ]{ }{ A } {u [ M ] i
i
i
i
i
i
i
i
i
g
T
}{ }
(9.28)
i
Misal: M i { i} [ M ]{ i } T
Ri { i } [ M ] T
C 2 n [ M ] [ K ] n [ M ] 2
Sehingga persamaan (9.28) menjadi
u
[ M ]{u
g
} [C ]{u } [ K
]{u} 0 M iA i 2 A u
M
2 M
i
A
i
i
i
A i
u
i
R g
i
2 A 2 A
………(9.29)
Ri
i
i
i
i
i
g i
Apabila persamaan (9.29) ditulis D
i
2 i D
i
i2 Di u
Dari persamaan (9.29) dan (9.30) memberikan R Ai i Di D S D pseudodisplecement M i Sehingga nilai perpindahan relatif maksimum umax
({ i} Ai ) 2
0.5
………(9.30)
………(9.31)
………(9.32)
DINAMIKA STRUKTUR
85
Contoh 9.3 Diketahui struktur portal tingkat tiga dengan pembebanan, berat per lantai dan kekakuan kolom seperti tergambar. W1
m1
W2 m2
w1 = 2943 kg, K 1 = 1600 kg/cm w2 = 4414 kg, K 2 = 2000 kg/cm
W3 m3
w3 = 4414 kg, K 3 = 2400 kg/cm
Hitung : 1. Frekuensi alami dan waktu getar alami dari sistem struktur di atas. 2. Gambar mode shape dari masing-masing waktu getar alami yang terjadi. 3. Hitung gaya gempa disetiap lantai dari sistem struktur tersebut jika berada di wilayah gempa 3 dengan jenis tanah lunak SNI .
Solusi: Solusi : a. Menghitung massa beban tiap lantai W 1
W 1 2943 kg
m1
W 2 4414 kg
m2
W 2
W 1 4414 kg
m3
W 3
g
g g
2943 981 4414 981 4414 981
3 4,5 4,5
DINAMIKA STRUKTUR
b. Menyusun matriks kekakuan [ K ] k12 = -k1
k11 = k1
m1
k13 = 0
k23 = -k2
k21 =- k1 m2
k22 = k1+k2
m3
k31 = 0
k32 = -k2
K 12 K 13 K 1 K 11 K 1 K K21 K 22 K 23 K1 K1 K 2 K 31 K 32 K 33 0 K 2 1600 0 1600 1600 1600 2000 2000 0 2000 2000 2400 0 1600 1600 1600 3600 2000 2000 4400 0 Menyusun matriks massa [m] 0 0 3 m1 M 0 m 0 0
0
m3 0
2
0
c. Persamaan Frekuensi
K M 0 2
n
0 K M 0 2
n
0 4,5 0
4,5 0 0
k33 = k2+k3
2 K 2 K 3 0 K
86
DINAMIKA STRUKTUR
87
0 0 0 1600 1600 3 2 2000 n 0 4,5 0 0 1600 3600 0 0 2000 0 4,5 4400 1600 3 2 1600 0 n 3600 4,5 2n 2000 2 0 1600 4400 4,5 0 2000 n 2 1600 3 1600 0 n 2 36004,5 n 2000 2 0 det 1600 0 2000 4400 4,5 n 1600 3 n 2 3600 4,5n 2 4400 4,5 2 160016004400 4,5 2 n n
1600 3n 2 2000 2000 0 6 2311 4 1343210 2 1,26410 8 0 n
n
n
Misalkan n , maka akan diperoleh persamaan 2
1343210 1,264108 0 1 21 116 rad 2 /det 2 1 1 3 2311 2
T1
757 rad
2
2
2
/det
2
2
2
1
2 10,8
2 2
3 23 1438 rad 2 /det 2
116 10,8 rad/det
0,58 det
757
27,5 rad/det
T2
2 2 0,23det
1
3
T3
2 2 0,17 det
27,5
2
1438 37,92 rad/det
37,92
3
Bentuk mode 1 1 1 (vektor sha pe) relatif displacement. 2 K M 0 K 2 M 0
n
1
1
DINAMIKA STRUKTUR
88
1600 0 0 0 11 1600 3 0 n 2 4,5 0 1600 3600 12 4400 0 0 2000 2000 0 4,5 0 13 1600 0 1600 3116 11 0 1600 2000 3600 4,5116 4400 4,5116 12 0 2000 13 0 1252 1600 11 1600 3078 2000 0 12 0 2000 3878 13 1252 11 1600 12 0 1600 11 3078 12 2000 13 0 2000 12 3878 13 0 Bila harga 11 1,0 dan harga-harga yang lain dinyatakan terhadap harga 11, maka diperoleh
1600
0 1252 0,7825 12 1600 2000 12 3878 13 0 20000,7825 3878 13 0
1252 1
12
13
2000 0,7825
0,4036
3878 Sehingga diperoleh mode shape relatif displacement sebagai berikut : 11 1 1 12 0,7825 0,4036
13
Dengan cara yang sama dan dengan menggantikan atau memberikan harga 0 , maka akan harga 2 dan 3 adalam persamaan K 2 M
n
dapat pula :
21 1 0,419 22 2 0,9844 23 1 31 1,696 32 3
33
1,634
DINAMIKA STRUKTUR
W1
m1
W2 m2
W3
1.00
1.00
0.782
- 0.419
0.403
m3
1.00
- 1.696
1.634
- 0.984
T1 0,58det
1 10,8 rad/det
T2 0,23det
2 27,5 rad/det
89
T3 0,17 det
1 37,92 rad/det
Kontrol kondisi orthogonal T M 0 1
2
3
0
1 0,7825 0,4036 0 4,5 0 0
0 1,0 0 0,419 0
4,5 0,984
Kontrol kondisi orthogonal untuk 2 dan 3 → 2T M 3 0
2 T M 3 0
1
3
0,419
0.984 0
0
0 4,5 0
0 1,0 0 1,696 0
4,5
1,634
Dari kontrol orthogonal tersebut di atas menunjukkan harga mode 1, 2, 3 sudah benar. d. Persamaan untuk harga
adalah harga relatif dari simpangan tiap-tiap
lantai, dan bagaimana harga mutlaknya, dapat dijelaskan sebagai berikut :
1 10,8rad/det T1 0,58 det cd 1 0.75 2 27,5rad/det T2 0,23det cd 2 0.75 3 37,92 rad/det T3 0,17 det cd 3 0.55 Dimana harga cd (respon percepatan maksimum dengan satuan “gravitasi”) diperoleh dari grafik koeffisien gempa dasar wilayah zone (SNI).
DINAMIKA STRUKTUR
90
Nilai perpindahan pola, R .c Ai Mi d 2
i S D
ii
dimana
i faktor S D
Ri
partisipasi
M i
pseudodisplacement
cd
i2
Menghitung nilai Ri : R1
1 T M
3
8,33 kg.det 2 /cm
0,403 4,5 4,5
1 0,782
R 2 2,65 kg.det /cm 2
R3 2,91kg.det /cm 2
Menghitung nilai M i : T M M 1
1
1
3
0,403 0 0
1 0,782
0 4,5 0
0 1
7,97 kg.det/cm 0 0,782 4,50,403
Dengan cara yang sama akan diperoleh : M 2 7,00 kg.det /cm 2
M 3 28,03 kg.det /cm 2
Menghitung Ai R .cd 1 A1 12 M
0,75 981 0,316 cm 2 8,33 10,8 7,97
R .c A2 22 Md 2
1
2
1
2
R .c A3 32 Md 3 3
3
2,65 0,75 981 2
0,368 cm 2.91 0,55 981 39,7 28.03 0,036 cm 27,5 7 2