Jbptitbpp Gdl Rudiasnann 30326 3 2008ta 2

Bab 2 Teori Dasar

BAB 2 DASAR TEORI

Gelombang air merupakan manifestasi dari suatu rambatan energi yang memiliki frekuensi dan periode. Gelombang air yang terjadi di laut dapat disebabkan oleh angin, gerakan kapal, gempa atau gaya gravitasi dari matahari dan bulan. Salah satu penyebab gelombang yang utama di laut adalah angin. Setelah gelombang terbentuk, gelombang tersebut akan merambat di permukaan laut. Dalam mempelajari perambatan gelombang ini digunakanlah teori gelombang, diantaranya adalah teori gelombang linier. Teori ini digunakan sebagai salah satu pendekatan dalam memodelkan perambatan gelombang secara matematis. Dalam teori ini digunakan beberapa asumsi untuk menyederhanakan permasalahan. Jika gelombang yang merambat mengenai struktur yang berada di perairan akan dikenai gaya yang ditimbulkan oleh gelombang. Gaya yang mengenai struktur tersebut diformulasikan dalam persamaan Morison, persamaan Froude Krylov dan teori difraksi.

2.1 TEORI GELOMBANG LINIER Pada pengembangan teori ini diasumsikan bahwa fluida dengan aliran stedi, tak mampu mampat (incompressible), dan irotasional. Selain itu berlaku juga hukum kekekalan massa yang merupakan dasar persamaan differensial gerak gelombang. Parameter terpenting dari gelombang adalah panjang gelombang (L) (jarak dari satu puncak gelombang ke puncak gelombang lainnya yang berurutan), tinggi gelombang H ( jarak vertikal dari lembah gelombang ke puncak gelombang), periode gelombang T (waktu yang ditempuh untuk satu lintasan gelombang), dan kedalaman perairan tempat gelombang menjalar. Untuk perumusan teori gelombang digunakan koordinat kartesian, dimana x searah dengan gerak gelombang, z diukur dari permukaan laut rata‐rata keatas secara vertikal, dan y tegak lurus terhadap x dan z. Diasumsikan gelombang dalam dua dimensi pada bidang x‐z, bergerak dalam arah x positif serta bergerak di permukaan laut yang datar dengan kedalaman yang konstan dimana selama bergerak bentuk gelombang tidak berubah. Selain itu fluidanya diasumsikan sebagai fluida seragam dan tidak mampu mampat (incompressible), sehingga kerapatan fluida

ρ

konstan serta pada

permukaan bebas tekanan seragam dan konstan. Asumsi lainnya adalah tegangan permukaan diabaikan. Di bawah ini akan digambarkan profil gelombang dua dimensi yang bergerak dalam arah x. Elevasi muka air

η merupakan fungsi

dari ruang dan waktu (x,t)

Laporan Tugas Akhir “Respon Dinamik Struktur Terapung”

2‐1

Bab 2 Teori Dasar

Dean Dean & Dal Dalrr m le, le, 199 1991 1 Gambar 2. 1 Karakteristik gelombang.

Asumsi dasar yang diterapkan adalah fluida yang tak mampu mampat (incompressible) dan tidak kental (inviscid), selain itu tidak terjadi gerak berputar fluida (irrotational motion), ada potensial kecepatan (velocity potential) yang memenuhi persamaan kontinuitas.

(incompressible)

(2.1)

Atau

(2.2)

(irrotational)

(2.3)

Dengan = vektor kecepatan = potensial kecepatan Persamaan (2.2) ditulis dalam persamaan Laplace sebagai berikut

(2.4) Untuk dua dimensi x dan z, persaman Laplace ditulis menjadi

(2.5)


2‐2

Bab 2 Teori Dasar

Persamaan Laplace merupakan persamaan pengatur BVP (Boundary Value Problem) pada Gambar 2.2 dengan syarat batas sebagai berikut: 1. Syarat batas dasar perairan (the bottom boundary condition, BBC)

pada z = ‐h

(2.6)

2. Syarat batas permukaan:

•

Syarat batas kinematis (kinematic free surface boundary condition, KFSBC)

pada z =

•

(2.7)

Syarat batas dinamis (dynamic free surface boundary condition, DFSBC)

pada z= 3.

η (x,t)

η (x,t)

(2.8)

Syarat batas periodik: ( x,t ) =

( x+L,t )

( x,t ) =

( x,t+T ) (2.9)

Gambar 2. 2 Sketsa definisi masalah nilai batas untuk teori gelombang linier.

Solusi yang tepat dari persamaan diatas sangat sulit untuk ditentukan karena syarat batas di permukaan memiliki suku‐suku tak linier serta kondisi awal di permukaan, z =

η (x,t),

tidak diketahui.

Oleh karenanya dilakukan penyederhanaan dengan melinierkan suku‐suku tak linier. Pelinieran dilakukan dengan membuat asumsi tinggi gelombang H jauh lebih kecil jika dibandingkan dengan Laporan Tugas Akhir “Respon Dinamik Struktur Terapung”

2‐3

Bab 2 Teori Dasar

panjang gelombang L dan kedalaman h, jadi H<
•

Syarat batas kinematis (kinematic free surface boundary condition, KFSBC)

pada z = 0

•

(2.10)

Syarat batas dinamis (dynamic free surface boundary condition, DFSBC)

pada z= 0

(2.11)

Dengan memanfaatkan syarat batas yang baru tersebut, kita dapat menyelesaikan persamaan Laplace di atas dengan menggunakan metoda pemisahan variabel (separation of variables method) sehingga untuk gelombang berjalan didapat potensial kecepatan.

(2.12) Dari syarat batas dinamis, dengan membuat rata‐rata dari η=0 maka C(t) =0 sehingga

pada z=0

(2.13)

Atau

(2.14) Sementara itu dari syarat batas kinematis dengan mensubsitusikan persamaan (2.12) dan (2.14) ke persamaan tersebut, didapat persamaan baru yang disebut persamaan dispersi, yaitu

(2.15) Gelombang berjalan membutuhkan satu perioda T untuk menempuh satu panjang gelombang L, dengan

ω=

2π/T serta k = 2π/L maka cepat rambat gelombang dapat dituliskan dalam

persamaan berikut


2‐4

Bab 2 Teori Dasar

(2.16) Atau

(2.17) 2

Pada laut dalam, kh besar dan tanh 2πh/L = 1, sehingga L= L0 =gT /2π dimana L0 adalah panjang gelombang di laut dalam. Secara umum ditulis sebagai berikut:

(2.18) Jadi panjang gelombang berkurang dengan berkurangnya kedalaman untuk perioda yang konstan.

2.2 KOMPONEN KECEPATAN DAN PERCEPATAN PARTIKEL Kecepatan partikel dapat dihitung dengan diketahuinya potensial kecepatan, dimana kecepatan partikel air merupakan turunan pertama dari potensial kecepatan. Kecepatan horizontal partikel u adalah:

(2.19) Atau

(2.20) Dari kecepatan dapat dihitung percepatan partikel air, dimana percepatan partikel air merupakan turunan pertama kecepatan terhadap waktu. Percepatan lokal dalam arah horizontal adalah

(2.21) Dan kecepatan vertikal w serta percepatannya adalah

(2.22)


2‐5

Bab 2 Teori Dasar

(2.23) Kecepatan dan percepatan merupakan fungsi dari posisi sehingga terdapat beda fase sebesar 0

90 . Kecepatan horizontal akan mempunyai nilai yang ekstrim pada saat posisi fasenya (kx‐ωt) = 0,π,… atau di bawah puncak dan lembah gelombang sementara percepatan mempunyai nilai ekstrim pada saat π/2,3π/2,… atau pada waktu perpindahan muka air adalah nol.

2.3 GAYA GELOMBANG Setiap struktur yang terdapat dilaut akan dibebani oleh gaya gelombang, untuk menghitung gaya tersebut maka diperlukan penerapan mekanika gelombang. Struktur yang sering digunakan di daerah pantai dan lepas pantai adalah struktur yang memanfaatkan tiang sebagai penyangga. Dalam perhitungan gaya gelombang dapat digunakan beberapa rumusan diantaranya adalah persamaan Morison, persamaan Froude Krylov dan teori difraksi. Persamaan‐persamaan ini dapat digunakan berdasarkan ukuran dari struktur yang dilewati oleh gelombang tersebut. Untuk struktur yang berukuran kecil dimana D (diameter tiang penopang dari struktur tersebut) memiliki ukuran yang jauh lebih kecil dari L (panjang gelombangnya) maka perhitungan gaya gelombangnya dapat dihitung dengan persamaan Morison dan persamaan Froude Krylov. Sedangkan untuk struktur yang memiliki ukuran yang besar dimana diameter tiang penopang dari struktur tersebut memiliki ukuran yang tidak berbeda jauh dari panjang gelombangnya maka perhitungan gaya gelombangnya dapat dihitung dengan teori difraksi. Berikut ini akan dijelaskan penggunaan persamaan perhitungan gaya‐ gaya tersebut:

•

Untuk struktur yang berukuran kecil “small body” (D/L

0.2 ) dapat digunakan persamaan

Morison.

• •

Jika struktur (0.2 < D/L Sedangkan jika struktur

0.5) maka dapat digunakan persamaan Froude Krylov. memiliki ukuran yang besar “large body” (D/L > 0.5) dapat

digunakan teori difraksi. Teori difraksi digunakan apabila struktur yang akan dilalui cukup besar sehingga akan mempengaruhi gelombang. Struktur yang besar ini mengakibatkan gelombang akan menghambur ketika melewati struktur sehingga terjadi perubahan pola aliran dari gelombang. Gaya gelombang Morison dan Froude Krylov mensyaratkan bahwa gelombang tidak pecah saat melewati struktur. Selain itu, pengaruh keberadaan tiang terhadap aliran akibat gelombang diabaikan karena tidak merubah pola aliran. Asumsi ini dipakai karena diameter tiang (D) lebih kecil dibandingkan panjang gelombang (L). Dua koefisien gaya, koefisien inersia (CM) dan koefisien seret (CD), digunakan untuk menentukan hubungan gaya dan gelombang.

2.3.1

Persamaan Morison

Morison menyatakan bahwa gaya gelombang dapat diekspresikan sebagai penjumlahan dari gaya seret (drag force) akibat kecepatan partikel air saat melewati struktur dan gaya inersia (inertia force) akibat percepatan partikel air.


2‐6

Bab 2 Teori Dasar

(2.24)

(2.25) Persamaan Morison dibatasi oleh syarat yang dinyatakan sebagai berikut:

Dimana D = diameter struktur L = Panjang gelombang Hal ini berarti bahwa struktur lebih kecil dari panjang gelombang. Untuk menentukan gaya total (F), maka persamaan Morison perlu diintegrasikan sepanjang struktur yang terendam.

(2.26)

(2.27) dimana F

: gaya hidrodinamik per satuan panjang

CD

: koefisien seret

CM

: koefisien inersia : volume struktur yang terkena gaya

A

: luas proyeksi dari struktur arah normal aliran : kecepatan pertikel air, tegak lurus terhadap sumbu struktur : percepatan partikel air, tegak lurus terhadap elemen struktur


2‐7

Bab 2 Teori Dasar

2.3.2

Persamaan Froude Krylov

Persamaan Froude Krylov berlaku untuk menghitung gaya gelombang yang bekerja pada struktur yang cukup besar, namun keberadaan struktur ini tidak memberikan perubahan yang berarti pada pola aliran. Persamaan Froude Krylov dibatasi oleh syarat yang dinyatakan sebagai berikut:

Dimana D = diameter struktur L = Panjang gelombang Persamaan Froude Krylov dinyatakan dalam persamaan matematis sebagai berikut:

(2.28) Dimana

(2.29) Maka persamaan 2.28 dapat ditulis menjadi:

(2.30) dimana F

: gaya hidrodinamik per satuan panjang : massa jenis massa jenis fluida

g

: percepatan gravitasi

H

: tinggi gelombang

k

: bilangan gelombang

n

: unit vektor normal


2‐8

Bab 2 Teori Dasar

2.3.3

Teori difraksi

Teori difraksi dapat digunakan dalam menghitung gaya gelombang yang bekerja pada struktur yang cukup besar. Akibat dari adanya struktur tersebut maka gelombang akan terdifraksi dan terjadi perubahan pada pola aliran. Akibat terjadinya perubahan pola aliran maka terdapat potensial kecepatan yang baru. Potensial kecepatan yang baru didapat dengan menggunakan prinsip super posisi, berikut ini akan akan disajikan dalam bentuk matematis:


2‐9

Bab 2 Teori Dasar

(2.31) dimana : potensial kecepatan total : potensial kecepatan gelombang datang : potensial kecepatan yang menyebar (scattered) Maka gaya gelombang yang bekerja pada struktur dapat dihitung dengan mensubsitusikan persamaan 2.28 dan 2.29 dengan persamaan 2.31, sehingga didapatkan pesamaan berikut:

(2.32)

2.4 STRUKTUR TERAPUNG Setiap struktur yang tedapat dilaut akan dipengaruhi oleh gelombang. Respon dari struktur tersebut akibat adanya gelombang dapat ditentukan dengan menganalisa struktur berdasarkan sudut pandang hidrodinamika. Dalam kasus ini struktur terapung diasumsikan berada pada daerah yang memiliki gelombang yang teratur (regular wave) yang menyerupai bentuk sinusoidal. Untuk mendapatkan total gaya hidrodinamika yang bekerja pada struktur terapung maka perlu kita ketahui komponen penyusun dari gaya hidrodinamika tersebut, diantaranya adalah:

• •

Beban eksitasi (exciation load) Beban radiasi (radiation load)

Ilustrasi dari total gaya hirodinamika tersebut dapat dilihat pada gambar 2.3 berikut ini:

Gambar 2. 3 Ilustrasi total gaya hirodinamika.


2‐10

Bab 2 Teori Dasar

Dalam bentuk matematis total gaya hidrodinamika pada struktur terapung dapat ditulis menjadi:

Untuk ( j = 1,2, … ,6)

(2.33)

Dimana :

“Wave excitatioan forces” yang terdiri dari gaya gelombang Froude Krylov dan gaya gelombang difraksi

:

“Radiation forces” yang terdiri dari “addede mass” dan potensial damping

:

“Viscous forces” yang terdiri dari “skin friction” dan “eddy making”

:

“Hirostatic forces” yang terdiri dari gaya gravitasi dan gaya apung (buoyancy)

Setiap komponen gaya memiliki potensial kecepatan yang berbeda‐beda, penjumlahan dari masing‐masing potensial kecepatan akan menghasilkan sebuah potensial kecepatan total. Dimana masing‐masing potensial kecepatan memenuhi syarat batas yang berlaku. Pengaruh gaya yang ditimbulkan akibat kekentalan fluida (viscous forces) dalam kasus ini dapat diabaikan, hal ini dikarenakan bilangan reynold (persamaan 2.34) memiliki nilai yang besar, sehingga kekentalan fluida tidak berpengaruh.

(2.34) Dimana u

: kecepatan aliran fluida

D

: diameter struktur : kekentalan fluida

Bilangan reynold memiliki nilai yang besar karena dalam kasus ini struktur yang dilewati oleh gelombang memiliki ukuran yang besar. Sedangkan gaya hidrostatik akan berpengaruh dari kekakuan dari stuktur. Untuk struktur terapung total potensial kecepatan merupakan penjumlahan dari potensial kecepatan gelombang datang

, potensial kecepatan difraksi

, yang merupakan

akibat dari gaya luar tanpa adanya pengaruh dari gerakaan struktur, dan juga potensial kecepatan radiasi

, yang merupakan hasil dari gerakan struktur. Dalam bentuk matematis total

potensial kecepatan dapat ditulis sebagai berikut: )


2‐11

Bab 2 Teori Dasar

(2.35)

Dalam bentuk kompleks persamaan 2.34 dapat ditulis menjadi:

(2.36) Dimana a adalah amplitudo gelombang dan

adalah frekuensi gelombang datang. Bentuk

kompleks dari potensial gelombang datang dinyatakan kecepatan difraksi dinyatakan radiasi dinyatakan terapung

2.4.1

. Bentuk kompleks dari potensial

. Sedangkan bentuk kompleks dari potensial kecepatan , dimana

adalah bentuk kompleks dari arah gerakan dari struktur

.

Beban eksitasi (excitation load)

Beban eksitasi merupakan beban luar yang pada struktur. Dalam permasalahan ini beban eksitasi diakibatkan oleh gelombang. Beban eksitasi terdiri dari beban akibat gelombang datang dan beban difraksi. Dimana besar gaya yang ditimbulkan oleh beban eksitasi merupakan penjumlahan dari gaya akibat gelombang datang dan gaya akibat difraksi, dalam rumusan matematis dapat ditulis menjadi:

(2.37) Sehingga didapatkan gaya eksitasi sebagai berikut

(2.38)

2.4.2

Beban radiasi (radiation load)

Beban radiasi muncul akibat

pergerakan dari struktur terapung sehingga menimbulkan

potensial kecepatan baru. Potensial kecepatan ini dikenal sebagai potensial kecepatan radiasi. Gaya yang ditimbulkan oleh beban radiasi ini dinyatakan dalam bentuk matematis sebagai berikut:

(2.39) Potensial kecepatan radiasi dapat dinyatakan dalam persamaan matematis sebagai berikut:


2‐12

Bab 2 Teori Dasar

(2.40) Dimana

adalah arah gerakan dari strukturuntuk arah j, sedangkan

adalah potensial

kecepatan radiasi untuk arah j. arah j.

Dari persamaan 2.39 bentuk riil dari persamaan ini disebut sebagai koefisien massa tambah sedangkan bentuk imajiner dari persamaan ini disebut sebagai koefisien redaman.

(2.41)

(2.42) Sehingga didapatkan bentuk lain dari persamaan gaya yang ditimbulkan oleh beban radiasi adalah:

(2.41) Dimana

merupakan koefisien massa tambah, sedangkan

adalah koefisien redaman.

2.5 STRIP THEORY

Strip theory merupakan suatu metode perhitungan masalah dinamika dari suatu struktur

terapung, dimana gaya‐gaya yang bekerja dan respon suatu struktur tiga dimensi tersebut dapat ditentukan dengan menggunakan hasil dari teori potensial dua dimensi. Teori ini mengasumsikan bahwa struktur ramping atau dengan kata lain ukuran panjang struktur jauh struktur jauh lebih panjang dari lebar dan draft struktur terapung tersebut. Prinsip dari teori ini adalah dengan membagi struktur terapung menjadi beberapa bagian ( 20 bagian), kemudian mengkombinasikan seluruh koefisien dari setiap bagian struktur tersebut untuk memperoleh koefisien “added mass” dan koefisien “damping”.


2‐13

Bab 2 Teori Dasar

Gambar 2. 4 Potongan 2 dimensi dari kapal (struktur terapung).

Koefisien‐koefisien tersebut dapat dibentuk dalam matriks yang berukuran 6

6, dimana

ukuran matriks ini menandakan struktur bergerak dalam 6 derajat kebebasan. Berikut ini adalah rumusan untuk menghitung masing‐masing koefisien pembentuk matriks persamaan gerak dinamik: Untuk arah vertikal:

A11 = a11dx

∫

(2.42)

B11 = b11dx

∫

(2.43)

∫

(2.44)

∫

(2.45)

A13 = A31 = a13dx B13 = B31 = b13dx A15

= − ∫ a13dx −

U0

B13

(2.46)

B15 = − xb13dx + U0 A13

(2.47)

ωe2

∫

A 51

= − ∫ a13dx +

U0

B13

(2.48)

B15 = − xb13dx − U0 A13

(2.49)

ωe2

∫

A33 = a33dx

∫

(2.50)

∫

(2.51)

B33 = b33dx


2‐14

Bab 2 Teori Dasar

A 35

= − ∫ xa33dx −

U0

ωe2

B 33

(2.52)

∫

B35 = − xb33dx − U0 A33 A 53

= − ∫ xa33dx +

U0

ωe2

(2.53)

B33

(2.54)

∫

B53 = − xb33dx − U0 A33

A 55

= −∫ x

B55

= − ∫ x b33dx +

2

U02

a33dx +

(2.55)

A 33

(2.56)

B33

(2.57)

C33 = c33dx = ρg B(x)dx

(2.58)

ωe2 U02

2

∫

ωe2

∫

∫

∫

C35 = C53 = xc33dx = −ρg x B(x)dx

∫

∫

C55 = x2c33dx = ρg x2B(x)dx

(2.59) (2.60)

Untuk arah horizontal:

A22 = a22dx

∫

(2.61)

∫

(2.62)

B22 = b22dx A24

= A 42 = ∫ a24dx

(2.63)

B24

= B42 = ∫ b24dx

(2.64)

A26

= − ∫ xa22dx +

B26

= −∫ xb22dx − U0 A22

A 62

= − ∫ xa22dx −

U0

ωe2

U0

ωe2

B22

B22

∫

B62 = − xb22dx + U0 A22


(2.65)

(2.66)

(2.67)

(2.68)

2‐15

Bab 2 Teori Dasar

A44

= ∫ a44dx

(2.69)

B44

= ∫ b44dx

(2.70)

A 46

= ∫ xa24dx +

B46

= ∫ xb24dx − U0 A24

A 64

= ∫ xa24dx −

B64

= ∫ xb24dx + U0 A24

A 66

= ∫x

B66

= ∫ x b66 dx +

C44

≈ ρg∇GMT

2

U0

ωe2

U0

ωe2

a22dx +

2

B24

(2.71)

(2.72)

B24

U0 2

ωe2 U02

ωe2

(2.73)

(2.74)

A22

(2.75)

B22

(2.76)

(2.77)

Sebagai contoh penggunaan teori strip dalam menentukan massa tambah dan koefisien redaman. Dimisalkan sebuah silinder dengan panjang L dan memiliki jari memiliki jari‐ jari sebesar R. lihat gambar 2.5.

Gambar 2. 5 Penampang melintang buoy berbentuk silinder.

Pusat koordinat berada tepat di pusat berat dari silinder dan permukaan air berada tepat pada sumbu x dari silinder tersebut. Untuk mendapatkan masa tambah pada arah pitch A55, dengan meninjau potongan sebesar dx yang terpengaruh oleh percepatan dalam arah vertikal ‐x akan terdapat massa tambah pada potongan tersebut yang besarnya

. Maka

. Gaya yang

ditimbulkan oleh massa tambah ini akan menimbulkan momen dalam arah y. Dengan mengintegrasikan panjang total dari silinder maka akan didapatkan momen total dalam arah pitch


2‐16

Bab 2 Teori Dasar

sebesar

(2.47) Berdasarkan persamaan 2.41 maka momen total dalam arah pitch ini sama dengan ‐A55 Hal ini berarti

.

. Sedangkan koefisien redamannya bernilai nol karena tidak ada

permukaan bebas.


2‐17

Bab 2 Teori Dasar

Contents BAB 2 ............................................... .............................................................. .......................................... ................. ......................... 1 DASAR TEORI ......................................................... ............................................................. ..................... ........... .......... 1 2.1

TEORI GELOMBANG LINIER .................................................................................................... 1

2.2

KOMPONEN KECEPATAN DAN PERCEPATAN PARTIKEL ....................................................... ... 5

2.3

GAYA GELOMBANG GELOMBANG ..................................................... ..................................................... ....... 6

2.3.1

Persamaan Morison ............................................................................................................ 6

2.3.2

Persamaan Froude Krylov ................................................................................................... 8

2.3.3

Teori difraksi........................................................................................................................ 9

2.4

STRUKTUR TERAPUNG........................................................................................................... 10

2.4.1

Beban eksitasi (excitation load) ........................................................................................ 12

2.4.2

Beban radiasi (radiation load) ................................................... ........................................ 12

2.5

....................................................................................................................... 13 STRIP THEORY .........................................................................................................................

Gambar 2. 1 Karakteristik gelombang. .................................................................................................... 2 Gambar 2. 2 Sketsa definisi masalah nilai batas untuk teori gelombang linier. ..................................... 3 Gambar 2. 3 Ilustrasi total gaya hirodinamika. ................................................ ..................................... 10 Gambar 2. 4 Potongan 2 dimensi dari kapal (struktur terapung). ....................................................... . 14 Gambar 2. 5 Penampang melintang buoy berbentuk silinder. ..................................................... ........ 16


2‐18

Jbptitbpp Gdl Rudiasnann 30326 3 2008ta 2

Recommend Documents