1 Respons Struktur SDOF Akibat Beban Sinusoidal Dengan Metode Integral Duhamel
1 Respons Struktur SDOF Akibat Beban Sinusoidal Dengan Metode Integral Duhamel
1 Respons Struktur SDOF Akibat Beban Sinusoidal Dengan Metode Integral Duhamel
Reni Suryanita
Jurusan Teknik Sipil Fakultas Teknik Universitas Riau
Kampus Bina Widya Jl. H.R. Soebrantas Km. 12,5 Pekanbaru
email : [email protected]
Hendra Sarfika
Jurusan Teknik Sipil Fakultas Teknik Universitas Riau
email : [email protected]
ABSTRAKSI
Dalam bidang teknik, penyederhanaan proses perhitungan sering dilakukan untuk memudahkan
pengolahan data. Beberapa struktur bangunan seperti rumah, tower dan bangunan lainnya
dalam proses perhitungan dapat disederhanakan dengan mengidealisasikan bangunan tersebut
sebagai sistem dengan derajat kebebasan tunggal (SDOF). Dalam sistem SDOF respons
struktur yang terjadi akibat beban dinamik diasumsikan searah horizontal. Untuk mendapatkan
respons yang terjadi disebabkan oleh beban dinamik seperti perpindahan, kecepatan dan
percepatan maksimum dapat digunakan konsep respons struktur. Dalam penelitian ini respons
struktur dihitung secara numerik menggunakan Integral Duhamel.
Langkah perhitungan dimulai dengan menempatkan gaya luar pada sistem struktur dan dihitung
nilai respons yang dihasilkan. Perhitungan dapat diulang untuk beberapa parameter yang
berbeda. Parameter yang divariasikan adalah massa dan kekakuan system. Respons struktur
dihitung dengan tipe pembebanan Sinusiodal untuk sistem teredam dan system tak teredam.
Hasil perhitungan memperlihatkan untuk kekakuan tetap perpindahan maksimum akan
meningkat pada saat massa sistem ditambah dua kali semula, dan akan menurun jika massa
system dikurangi setengah dari massa mula-mula. Hal ini berbanding terbalik terhadap variasi
kekakuan dengan massa sistem tetap.
Kata kunci : respons struktur, Integral Duhamel, pembebanan sinusoidal, massa, kekakuan
ABSTRACT
In engineering field, the simplicity process is often done in order to make easier the processing.
Some buildings e.g houses, tower and others can be idealized as Single Degree of Freedom
System (SDOF) with an assumption that the response caused by dynamical load only taken
place in horizontall direction. To detect the response occurs cause by dynamical load, e.g
displacement, velocity, and acceleration maximum structure, therefore the response of structure
concept can be used. In this case, response of the structure is calculated numerically using
Duhamel’s Integral.
First, the external force put on the system, both on mass and restrain and then calculating the
response that occurred. The calculation can be repeat for any different parameters. The
parameter variation are mass and system stiffness. The response of the structure is calculated
with sinusoidal loading type for dumping system and undumping system.
266
Volume 7 No. 3, Juni 2007 : 266 – 278
The study results showed that the maximum removal will increase at the mass added twice and
decrease at the reducing load a half. The removal is in proportion to the mass of the system and
inversely proportion to the stiffness of the system.
1. PENDAHULUAN
Respons struktur merupakan riwayat waktu dari perpindahan, kecepatan dan
percepatan dari fungsi beban tertentu, untuk struktur dengan derajat kebebasan tunggal dan
banyak. Pada permasalahan beban dinamis seperti beban ledakan, beban angin, beban getaran
mesin dan beban gempa, beban dan respons strukturnya merupakan fungsi dari waktu,
sehingga analisis yang dilakukan harus berdasarkan waktu.
Pembebanan pada struktur akibat beban dinamis dapat terjadi sewaktu-waktu, maka
untuk perencanaan bangunan perlu diperhitungkan pengaruh beban ini. Adakalanya struktur
yang direncanakan harus menerima beban secara berulang-ulang (periodik) yang tidak
diperhitungkan sebelumnya. Pola beban periodik ini dapat berbentuk grafik fungsi sinus
maupun fungsi cosinus, yang dikenal dengan nama sinusoidal. Struktur bangunan yang
mengalami beban ini harus dievaluasi, apakah masih aman untuk ditempati atau harus
diperbaiki dan bahkan mungkin harus dibongkar. Untuk dapat memprediksi kemungkinan-
kemungkinan tersebut, maka perlu ditinjau bagaimana respons yang terjadi pada struktur
bangunan, jika bangunan tersebut menerima beban dinamik dengan pola sinusoidal.
Penulisan artikel ini bertujuan untuk menentukan respons maksimum yang terjadi
pada sistem derajat kebebasan tunggal akibat beban sinusoidal. Respons maksimum yang
didapatkan berguna sebagai acuan untuk merencanakan sistem struktur sehingga dapat
meminimalisir kemungkinan kerusakan pada sistem derajat kebebasan tunggal.
2. PERMASALAHAN
Umumnya beban dinamis sulit dinyatakan dalam suatu bentuk matematis yang
sederhana secara eksak. Untuk itu diperlukan analisis dinamis secara numerik. Simulasi
numerik dalam penelitian ini menggunakan bantuan komputer, agar respons maksimum dari
sistem berderajat kebebasan tunggal yang diinginkan, dapat dihitung lebih akurat, dan
membutuhkan waktu yang relatif lebih singkat.
3. BATASAN MASALAH
Beberapa batasan dalam penelitian ini antara lain:
1) Struktur dimodelkan mempunyai sistem derajat kebebasan tunggal (SDOF)
2) Dimensi balok diasumsikan 30 cm x 40 cm dan kolom dengan ukuran 20 cm x 30 cm.
3) Beban merata (q) yang bekerja pada struktur diasumsikan sebesar 2,5 t/m’. Modulus
elastik beton Eb = 2,2.105 kg/cm2 dan percepatan gravitasi g = 980 cm/dt2.
4) Beban dinamik yang bekerja pada stuktur merupakan fungsi sinusiodal sebesar P(t) =
1500 sin (πt/0,6) selama 1,2 detik.
5) Analisis sistem dilakukan dengan variasi massa dan kekakuan baik untuk sistem
teredam maupun tak teredam dengan pembebanan sinusoidal.
Respons Struktur SDOF Akibat Beban Sinusoidal Dengan Metode Integral Duhamel 267
(Reni Suryanita, Hendra Sarfika)
4. TINJAUAN PUSTAKA
Pada dasarnya setiap struktur mempunyai derajat kebebasan yang tak terhingga
jumlahnya. Menurut Widodo (2000) derajat kebebasan merupakan derajat independensi atau
jumlah koordinat yang diperlukan untuk menyatakan posisi suatu sistem pada setiap saat.
Suatu struktur memiliki frekuensi natural sebanyak derajat kebebasan yang dimilikinya dan
jika beban dinamik yang diterima struktur memiliki frekuensi yang mendekati frekuensi
natural dari struktur maka akan terjadi resonansi yang akan mengakibatkan keruntuhan atau
collapse pada struktur.
Lumantarna (1999), Sistem massa yang berpindah dalam satu arah saja yaitu arah
horizontal dinamakan sistem berderajat kebebasan tunggal (single degree of freedom, SDOF).
Pada sistem SDOF, struktur dimodelkan dengan massa tunggal dan koordinat perpindahan
tunggal.
5. LANDASAN TEORI
Secara sederhana model matematis sistem derajat kebebasan tunggal (SDOF) dapat
dilihat seperti Gambar 2 berikut ini,
268
Volume 7 No. 3, Juni 2007 : 266 – 278
Persamaan gerak sistem derajat kebebasan tunggal menghasilkan gaya luar, Ft berdasarkan
kesetimbangan dinamis free-body :
atau
y (t ) = {A(t ) sin ωt − B(t ) cos ωt}/ mω (2)
dengan
t
A(t ) = ∫ F (τ ) cos ωτdτ
0
t
B(t ) = ∫ F (τ ) sin ωτdτ (3)
0
Perhitungan Integral Duhamel memerlukan evaluasi integral A(t) dan B(t) secara
numerik. Metoda yang sering digunakan adalah aturan Trapesium (trapezoidal rule) dan
aturan Simpson (Simpson’s rule). Untuk mengintegrasi suatu fungsi I(τ)dapat diturunkan
sebagai berikut :
t
A(t) = ∫ I (τ )dτ
0
(4)
Respons Struktur SDOF Akibat Beban Sinusoidal Dengan Metode Integral Duhamel 269
(Reni Suryanita, Hendra Sarfika)
Dalam menggunakan metode ini, fungsi beban F(τ) dianggap mendekati fungsi linier.
Untuk menggambarkan respons secara lengkap integrasi pada Persamaan (7) dapat dinyatakan
dalam bentuk pertambahan (incremental form) berikut.
ti
dengan A(ti) dan B(ti) menyatakan harga dari integral-integral pada Persamaan (7) pada saat ti.
ΔF
F (τ ) = F (t i −1 ) + i (τ − t i −1 ) , t i −1 ≤ τ ≤ t i (8)
Δt i
dengan
ΔFi = F (t i ) − F (t i −1 ) (9)
dan
Δt i = t i − t i −1 (10)
dengan mensubstitusikan Persamaan (8) ke Persamaan (7) diperoleh
⎛ ΔF ⎞⎛ sin ωt i − sin ωt i −1 ⎞
A(t i ) = A(t i −1 ) + ⎜⎜ F (t i −1 ) − t i −1 i ⎟⎟⎜ ⎟
⎝ Δt i ⎠⎝ ω ⎠
(11)
ΔF
+ 2 i {cos ωt i − cos ωt i −1 + ω (t i sin ωt i − t i −1 sin ωt i −1 )}
ω Δt i
dan
⎛ ΔF ⎞⎛ cos ωt i − cos ωt i −1 ⎞
B(t i ) = B(t i −1 ) + ⎜⎜ F (t i −1 ) − t i −1 i ⎟⎟⎜ ⎟
⎝ Δt i ⎠⎝ ω ⎠
(12)
ΔF
+ 2 i {sin ωt i − sin ωt i −1 + ω (t i cos ωt i − t i −1 cos ωt i −1 )}
ω Δt i
F (τ )dτ
dy (t ) = e ξω ( t −τ ) sin ω D (t − τ ) (13)
mω D
dan jumlah respons differensial dari semua waktu pembebanan adalah
t
1
mω D ∫0
y (t ) = F (τ )e −ξω (t −τ ) sin ω D (t − τ )dτ
270
Volume 7 No. 3, Juni 2007 : 266 – 278
ti
AD (t i ) = AD (t i −1 ) + ∫ F (τ )e
ξωτ
cos ω Dτdτ (15)
t i −1
ti
BD (t i ) = BD (t i −1 ) + ∫ F (τ )e
ξωτ
sin ω Dτdτ (16)
t i −1
Untuk semua bagian fungsi pembebanan linier, F(τ) dari Persamaan (8)
disubstitusikan ke Persamaan (15) dan (16) memerlukan evaluasi integral berikut :
e ξωτ
ti
∫ e cos ω Dτdτ =
ξωτ
I1 = (ξω cos ω Dτ + ω D sin ω Dτ ti
(17)
ti −1 (ξω )2 + ω D2 t i −1
e ξωτ
ti
∫ e sin ω Dτdτ =
ξωτ
I2 = (ξω sin ω Dτ − ω D cos ω Dτ ti
(18)
ti −1 (ξω )2 + ω D2 t i −1
ti
⎛ ξω ⎞ ωD
I 3 = ∫ τe ξωτ sin ω Dτdτ = ⎜⎜τ − ⎟I +
2 ⎟ 2
I ti (19)
t i −1 ⎝ (ξω ) + ω D ⎠ (ξω ) + ω D
2 2 2 1 ti −1
ti
⎛ ξω ⎞ ωD
I 4 = ∫ τe ξωτ cos ω Dτdτ = ⎜⎜τ − ⎟I −
2 ⎟ 1
I 2 ttii −1 (20)
ti −1 ⎝ (ξω ) + ω D ⎠ (ξω ) + ω D
2 2 2
⎛ ΔF ⎞ ΔF
AD (t i ) = AD (t i −1 ) + ⎜⎜ F (t i −1 ) − t i −1 i ⎟⎟ I i + i I 4 (21)
⎝ Δt i ⎠ Δt i
⎛ ΔF ⎞ ΔF
BD (t i ) = BD (t i −1 ) + ⎜⎜ F (t i −1 ) − t i −1 i ⎟⎟ I 2 + i I 3 (22)
⎝ Δt i ⎠ Δt i
dari persamaan-persamaan di atas diperoleh perpindahan pada saat ti yaitu,
e −ξωti
y (t i ) = {AD (t i ) sin ω D t i − BD (t i ) cos ω D t i } (23)
mω D
dengan
y(ti) = perpindahan pada titik i (cm)
ξ = rasio redaman
ω = frekuensi sudut (rad/dt)
m = massa (kgdt2/cm)
ωD = frekuensi sudut teredam (rad/dt)
ti = waktu pada titik i (dt)
Respons Struktur SDOF Akibat Beban Sinusoidal Dengan Metode Integral Duhamel 271
(Reni Suryanita, Hendra Sarfika)
Tabel 1. Langkah- langkah perhitungan respons struktur
Flowchart Keterangan
Input data :
NTYPE, indek pengaruh
N, jumlah titik yang ditentukan pada
fungsi pengaruh.
M, massa
K, kekakuan
XI, rasio redaman
e −εωt i Tmax, waktu maksimum dari integrasi
y(i) = AD(t i )sin (ω D t i ) - BD(ti )cos( ωD t i )
mωD DT, selang waktu dari integrasi
y max = y(i) max
INT, indek interpolasi
GR, percepatan gravitasi
T(I), waktu pada titik i
e − εωt i
v(i) = [(ωD BD t i − εωA D t i )sin (ωD t i ) + (ωD A D t i + εωBD t i )cos(ωD t i )] F(I), gaya pada titik i
mωD
v max = v(i) max
Perpindahan :
y(i), perpindahan pada titik i
F(i) - C * v(i) - K * y(i) ymak, perpindahan maksimum
ACC(i) =
m m, massa
ACC max = ACC(i) max
ε, redaman
ω, frekuensi natural
ti, waktu pada titik i
y max , VMAX , ACC MAX
ωD, frekuensi natural teredam
Kecepatan :
v(i), perpindahan pada titik i
vmak, perpindahan maksimum
m, massa
ε, redaman
ω, frekuensi natural
ti, waktu pada titik i
ωD, frekuensi natural teredam
Percepatan :
ACC(i), percepatan pada titik i
ACCmax, percepatan maksimum
C, koefisien redaman liat
K, konstanta pegas
m, massa
Data Perhitungan
Model struktur yang ditinjau mempunyai sistem derajat kebebasan tunggal dengan
dimensi balok 30 cm x 40 cm dan kolom dengan ukuran 20 cm x 30 cm. Struktur diberi
beban merata (q) sebesar 2,5 t/m’. Modulus elastik beton Eb = 2,2.105 kg/cm2 dan percepatan
272
Volume 7 No. 3, Juni 2007 : 266 – 278
gravitasi g = 980 cm/dt2. Struktur dibebani dengan beban dinamik fungsi sinusiodal sebesar
P(t) = 1500 sin (πt/0,6) selama 1,2 detik seperti pada Gambar. 3. di bawah ini,
2000
1000
gaya (kg)
0
0,0 0,2 0,4 0,6 0,8 1,0 1,2
-1000
Respons Struktur SDOF Akibat Beban Sinusoidal Dengan Metode Integral Duhamel 273
(Reni Suryanita, Hendra Sarfika)
Tabel 2. Variasi massa dan kekakuan
Massa Kekakuan
Variasi Keterangan
(kgdt2/cm) (kg/cm)
- 12,8 4027,3 Acuan
1 12,8 2013,6 Massa tetap, kekakuan setengah dari acuan
2 12,8 8054,6 Massa tetap, kekakuan dua kali dari acuan
3 6,4 4027,3 kekakuan tetap, massa setengah dari acuan
4 25,6 4027,3 kekakuan tetap, massa dua kali dari acuan
7.1. Perpindahan
a. Struktur dengan variasi massa dan kekakuan tetap sebesar 4027,3 kg/cm
Sistem dengan massa sebesar 6,4 kgdt2/cm menghasilkan perpindahan maksimum
sebesar 0,44294 cm, massa sebesar 12,8 kgdt2/cm menghasilkan perpindahan maksimum
sebesar 0,51648 cm dan untuk massa sebesar 25,6 kgdt2/cm menghasilkan perpindahan
maksimum sebesar 0,52151 cm.
1 m = 6,3875 kgdt2/cm
m = 12,775 kgdt2/cm
0,8 m = 25,550 kgdt2/cm
0,6
perpindahan (cm)
0,4
0,2
0
-0,2 0 0,2 0,4 0,6 0,8 1
-0,4
-0,6
-0,8 w aktu (detik)
Semakin besar massa maka semakin besar pula perpindahan maksimum yang
dihasilkan atau dengan kata lain massa berbanding lurus dengan perpindahan maksimum,
secara visual dapat dilihat pada Gambar 5.
b. Sistem dengan variasi kekakuan dan massa tetap sebesar 12,8 kgdt2/cm
Sistem dengan kekakuan sebesar 2013,6 kg/cm menghasilkan perpindahan maksimum
sebesar 1,04302 cm, kekakuan sebesar 4027,3 kg/cm menghasilkan perpindahan maksimum
sebesar 0,51648 cm dan untuk kekakuan sebesar 8054,6 kg/cm menghasilkan perpindahan
maksimum sebesar 0,22147 cm.
274
Volume 7 No. 3, Juni 2007 : 266 – 278
k = 2013,64 kg/cm
1,5 k = 4027,27 kg/cm
k = 8054,54 kg/cm
1
perpindahan (cm) 0,5
0
0 0,2 0,4 0,6 0,8 1
-0,5
-1
w aktu (detik)
-1,5
7.2. Kecepatan
a. Sistem dengan variasi massa dan kekakuan tetap sebesar 4027,3 kg/cm
Sistem dengan massa sebesar 6,4 kgdt2/cm menghasilkan kecepatan maksimum
sebesar -4,97423 cm/dt, massa sebesar 12,8 kgdt2/cm menghasilkan kecepatan maksimum
sebesar -6,3269 cm/dt dan untuk massa sebesar 25,6 kgdt2/cm menghasilkan kecepatan
maksimum sebesar -5,78066 cm/dt.
m = 6,3875 kgdt2/cm
10 m = 12,775 kgdt2/cm
m = 25,550 kgdt2/cm
kecepatan (cm/dt)
0
0 0,2 0,4 0,6 0,8 1
-5
b. Sistem dengan variasi kekakuan dan massa tetap sebesar 12,8 kgdt2/cm
Sistem dengan kekakuan sebesar 2013,6 kg/cm menghasilkan kecepatan maksimum
sebesar -11,5613 cm/dt, kekakuan sebesar 4027,3 kg/cm menghasilkan kecepatan maksimum
Respons Struktur SDOF Akibat Beban Sinusoidal Dengan Metode Integral Duhamel 275
(Reni Suryanita, Hendra Sarfika)
sebesar -6,3269 cm/dt dan untuk kekakuan sebesar 8054,6 kg/cm menghasilkan kecepatan
maksimum sebesar -2,48712 cm/dt.
k = 2013,64 kg/cm
15 k = 4027,27 kg/cm
k = 8054,54 kg/cm
10
kecepatan (cm/dt)
0
0 0,2 0,4 0,6 0,8 1
-5
-10
7.3. Percepatan
a. Sistem dengan variasi massa dan kekakuan tetap sebesar 4027,3 kg/cm
Sistem dengan massa sebesar 6,4 kgdt2/cm menghasilkan percepatan maksimum
sebesar -84,1396 cm/dt2, massa sebesar 12,8 kgdt2/cm menghasilkan percepatan maksimum
sebesar 100,7412 cm/s2 dan untuk massa sebesar 25,6 kgdt2/cm menghasilkan percepatan
maksimum sebesar 67,28318 cm/dt2.
276
Volume 7 No. 3, Juni 2007 : 266 – 278
m = 6,3875 kgdt2/cm
180 m = 12,775 kgdt2/cm
m = 25,550 kgdt2/cm
120
percepatan (cm/dt 2)
60
0
0 0,2 0,4 0,6 0,8 1
-60
-120
b. Sistem dengan variasi kekakuan dan massa tetap sebesar 12,8 kgdt2/cm
Sistem dengan kekakuan sebesar 2013,6 kg/cm menghasilkan percepatan maksimum
sebesar 134,5664 cm/dt2, kekakuan sebesar 4027,3 kg/cm menghasilkan percepatan
maksimum sebesar 100,7412 cm/dt2 dan untuk kekakuan sebesar 8054,6 kg/cm menghasilkan
percepatan maksimum sebesar -42,0698 cm/dt2.
k = 2013,64 kg/cm
200 k = 4027,27 kg/cm
150 k = 8054,54 kg/cm
percepatan (cm/dt 2)
100
50
0
-50 0 0,2 0,4 0,6 0,8 1
-100
-150
-200 w aktu (detik)
Respons Struktur SDOF Akibat Beban Sinusoidal Dengan Metode Integral Duhamel 277
(Reni Suryanita, Hendra Sarfika)
8. KESIMPULAN
Dari hasil dan analisis respon struktur SDOF pada penelitian ini dapat diambil
kesimpulan bahwa semakin besar massa maka semakin besar pula perpindahan maksimum
yang dihasilkan atau dengan kata lain massa berbanding lurus dengan perpindahan
maksimum. Lain halnya dengan kekakuan, semakin besar kekakuan maka semakin kecil
perpindahan maksimum yang dihasilkan atau dengan kata lain kekakuan berbanding terbalik
dengan perpindahan maksimum.
Variasi massa dan kekakuan mengasilkan kecepatan maksimum yang semua nilainya
negatif. Hal ini menunjukkan bahwa kecepatan maksimum terjadi pada arah yang berlawanan
dengan arah beban. Nilai kecepatan dan percepatan untuk tiap penambahan massa cenderung
bertambah besar sedangkan nilai kecepatan untuk tiap penambahan kekakuan cenderung akan
semakin kecil.
DAFTAR PUSTAKA
Clough, R.W. & Penzien, J., 1988, Dinamika Struktur, Erlangga, Jakarta.
Dewobroto, W., 2004, Aplikasi Rekayasa Konstruksi dengan SAP2000, PT Elex Komputindo
Kelompok Gramedia, Jakarta.
James,M.L. et al., 1994, Vibration of Mechanical and Structural Systems with Microcomputer
Application. Harper Collins College Publishers, New York.
Liu, J.L, 2002, New Method for Solution of Response of Nonlinear Hysteretic System,
Journal of Vibration and Acoustics, October 2002 , Volume 124, Issue 4, pp. 653-
656
Lumantarna, B.,1999, Pengantar Analisis Dinamis dan Gempa, Andi Offset Yogyakarta.
Paz, M.,1990., Dinamika Struktur. Erlangga, Jakarta.
Widodo.,2000., Respons Dinamik Struktur Elastik. UII Press, Yogyakarta.
278
Volume 7 No. 3, Juni 2007 : 266 – 278