## Jurnal APLIKASI
ISSN.1907-753X

## Methode Aplikasi Bangunan Krib Sebagai Pelindung terhadap Bahaya Erosi Tebing Sungai
Suharjoko Staft Pengajar Program Studi D-III Teknik Sipil FTSP ITS email: suharjoko_hydro@ce.its.ac.id

## ABSTRAK
Untuk tujuan disain suatu bangunan krib, perlu dilakukan pengkajian  dengan menggunakan model, sementara membangun model baik model fisik maupun model matematis memerlukan dana yang cukup besar. Oleh karena itu  Suharjoko ( 2001) telah mengkaji dengan melakukan simulasi numerik terhadap 45 kasus yaitu terdiri dari tiga model krib yang di dibedakan  oleh sudut antara krib dengan tebing sungai, setiap modelnya terdapat tiga variasi yang dibedakan oleh ukuran panjang krib dan setiap variasi dikaji terhadap lima kasus yang dibedakan oleh kecepatan arusnya. Tujuan penelitian adalah untuk mendapatkan pedoman perencanaan aplikasi bangunan krib yang sesuai.
Dari hasil simulasi diperoleh bahwa parameter penting yang menjadi penentu adalah lebar sungai (B), Kecepatan arus (v), kedalaman air (h), panjang krib (P), arah krib (  ) itu dipasang dan panjang perlindungan yang dihasilkan (D). Kemudian antara parameter- parameter penting yang berpengaruh dengan cara penyelesaian analisa dimensi dicari hubungannya. Dari penyelesaian tersebut didapat hubungan antara bilangan Froude ( Fr. ) dengan Dh/PB. Besarnya bilangan Froude (Fr) adalah V / gh
Hasil temuan pada studi ini menghasilkan suatu grafik hubungan antara bilangan Froude (Fr.) dengan Dh/PB sehingga memberikan kemudahan dalam perencanaan bangunan Krib, dengan demikian akan sangat menghemat biaya disain karena tahapan studi simulasi tidak lagi perlu dilakukan dalam merencanakan bangunan krib tersebut. 1. PENDAHULUAN

## A. Latar Belakang
Pada sungai dimana kecepatan arusnya tinggi tentu menimbulkan erosi, baik yang tejadi di tebing maupun di dasar sungai, oleh karena itu diperlukan suatu bangunan perlindungan. Bangunan perlindungan tebing dapat berupa plengsengan (lining), tembok, turap, bronjong, krib dan lain-lain. Bangunan Krib adalah salah satu bangunan perlindungan yang dipasang melintang pada tebing sungai, yang tujuannnya untuk memperlambat kecepatan arus di sekitar bangunan krib tersebut sehingga proses erosi akan terhindari bahkan akan terjadi proses sedimentasi. Dengan demikian pada tebing yang semula terjadi proses erosi menjadi terlindungi bahkan terjadi endapan (sedimentasi) akibat adanya bangunan krib tersebut. Namun untuk kebutuhan perencanaan bangunan krib perlu dilakukan analisa efektifitas tata letak bangunan krib tersebut apakah harus tegak lurus atau
menyudut dan berapa panjang krib yang baik, untuk itu perlu dilakukan pengkajian dan analisa terhadap faktor-faktor yang berpengaruh. Untuk melakukan pengkajian tersebut dapat dilakukan dengan suatu model, baik model fisik maupun model matematis sementara membangun model memerlukan dana yang cukup besar. Oleh karena itu Suharjoko, 2001 telah melakukan penelitian dengan membangun model matematik. Model ini merupakan suatu penyelesaian persamaan aliran tak tunak dua Dimensi Horizontal dengan methode karakteristik.
Model matematik tersebut dilakukan untuk mengkaji terhadap 45 kasus model krib yang terdiri dari tiga model yang di yang dibedakan  oleh sudut antara krib tersebut terhadap tebing sungai. Dan setiap modelnya terdapat tiga variasi yang dibedakan oleh ukuran panjang krib dan setiap variasi dikaji terhadap lima kasus yang dibedakan oleh kecepatan arusnya.

## Jurnal APLIKASI
ISSN.1907-753X
Hasil simulasi diperoleh bahwa parameter penting yang menjadi penentu adalah lebar sungai, kecepatan arus, kedalaman air, panjang krib, arah krib itu dipasang dan panjang perlindungan yang dihasilkan. Dari parameter yang dihasilkan tersebut dicari hubungannya sedemikian hingga diperoleh suatu fungsi atau hubungan antar parameter tersebut menjadi hubungan yang general. Dengan demikian diperoleh suatu methode perncanaan bangunan krib yang mudah dalam aplikasi disainnya dan tidak perlu lagi melalui kajian.

## 1.2. Tujuan dan Manfaat
Tujuannya adalah melakukan kajian sedemikian hingga diperoleh hubungan antara parameter penting yang berpengaruh tersebut terhadap perlindungan yang dihasilkan oleh bangunan krib terhadap tebing sungai sehingga diperoleh suatu fungsi atau hubungan antar parameter menjadi hubungan yang general.
Sedangkan manfatnya adalah diperoleh kemudahan dalam merencanakan bangunan krib tersebut sebagai perlindungan sungai sehingga dapat menghemat energi dan biaya yang semula dalam merencanakan bangunan krib harus melalui kajian dengan model menjadi hitungan yang amat sederhana.

## 2. METODOLOGI

## Tahap pertama dilakukan running model
terhadap berbagai kasus dan dilanjutkan analisa terhadap setiap hasil running model yang dihasilkan yakni melakukan penilaian terhadap besaran parameter yang dihasilkan.
Tahap kedua melakukan analisa non-dimensi terhadap parameter penentu untuk mendapatkan  hubungan antar parameter tersebut.
Tahap ketiga menghitung terhadap nilai parameter yang dihasilkan untuk mendapatkan hubungan antar parameter model dan kasus.
Tahap keempat melakukan analisa untuk mendapatkan hubungan antar parameter tersebut.

## 3. TINJAUAN PUSTAKA

## 3.1. Model Matematik
Model ini merupakan suatu penyelesaian persamaan aliran tak tunak (unsteady) dua Dimensi Horizontal dengan methode karakteristik. Persamaan aliran tak tunak dua Dimensi Horizontal adalah terdiri dari persamaan momentum dan persamaan kontinuitas dari Navier Stokes berikut ini
Persamaan kontinuitas :
................. (1) Persamaan momentum : ............... (2a) ............... (2b) dimana   . u & v = Kecepatan rerata vertikal arus arah x.dan arah y, h = Kedalaman sungai,  = Tinggi air terhadap bidang Persamaan atau tinggi pasang surut.
Persamaan di atas diturunkan mengikuti konsep gerakan partikel air dari tempat (x,y,z) menuju ke tempat (x+dx, y+dy, z+dz) dalam selang waktu dt. Pada Gambar 1 berikut ditunjukkan perubahan kecepatan partikel/gerak zat cair arah x,y,z dalam selang waktu dt  yaitu suatu perubahan gerakan pertikel dalam 3-Dimensi.
u x p X y u v x u u t u 2 1                 0 ) ( ) (          y hv x hu t  v y p Y y v v x v u t v 2 1                              2 2 2 2 2 y x

## Jurnal APLIKASI
ISSN.1907-753X
Kemudian dengan mengeliminasi gerakan arah vertical maka persamaan menjadi persamaan gerak aliran 2 Dimensi Horizontal seperti yang ditunjuk  pada persamaan 1 dan 2 diatas.
Selanjutnya persamaan defferensial dari persamaan aliran di atas diselesaikan secara numeris dengan methode karakteristik kompotabel Triatmadja (1990) yang merupakan
Hasil akhir dari penyelesaian tersebut diperoleh suatu skema karakteristik sebagai berikut :
       t hv hu vh uh y x y x h ) (            2 2 2 2 2 ) ( ) / )( ( 2 1   s t v u c   2 ) / ( ( 4 ) ( ) 2 2 2 ( 2 2 ) / ( s t h h gu u u v u c t s                 2 ) / ( ) ( 4 ) ( 4 s t v v hv u u hu      . ............. (3)   t g vu uu u u x y x h       t y v x v t h C v u gv z           ) ( ) ( 2 2 2 2 2 5 , 0 2 2  + t y u x u t h C v u gu z           ) ( ) ( 2 2 2 2 2 5 , 0 2 2  ,  ................ (4)         t g vv uv v v y y x h      2 2 2 2 2 ) / ( ( 4 ) ( ) ( 2 ) / ( s t h h gv v v v u c t s           t y v x v t h C v u gv z           ) ( ) ( 2 2 2 2 2 5 , 0 2 2  , (5)
Sesuai dengan jaring hitungan maka unsur- unsur yang tertulis pada persamaan di atas (Persamaan 3.24, 3.25 dan 3.26), dapat didefinisikan sebagai berikut:
2 2 x u   = (u (i+1,j) -2u (i,j) + u (i-1,j) )/  x 2 ,
2 2 y u   = (u (i,j+1) -2u (i,j) + u (i,j-1) )/  y 2 ,
2 2 x v   = (v (i+1,j) -2v (i,j) + v (i-1,j) )/  x 2 ,
2 2 y v   = (v (i,j+1) -2v (i,j) + v (i,j-1) )/  y 2 ,
h =  + ho,
 s =  x =  y,
c = Cepat rambat gelombang, ( g h (i,j) ) 0,5 ,
C z = Koefisien Chezy  ,
 = Viskositas kinematis = 1x10 -3 .
Persamaan numeric di atas kemudian ditulis dalam bahasa program komputer yang  telah diselesaikan oleh Suharjoko (1999), dengan menggunakan bahasa Program Visual Basic.

## 3.2. Analisa Kasus dan Model
Dengan menggunakan program tersebut telah dikaji (1999) terhadap 45 kasus yaitu terdiri dari tiga model krib yang di dibedakan  oleh sudut antara krib dengan tebing sungai, setiap modelnya terdapat tiga variasi yang dibedakan oleh ukuran panjang krib dan setiap variasi dikaji terhadap lima kasus yang dibedakan oleh kecepatan arusnya.
Hasil model diatas sebagai dicontohkan pada Gambar 2 adalah gambar vektor kecepatan disekitar Krib miring dengan kecepatan di hulu = 0.6 m/dt, Gambar 3 adalah gambar vektor kecepatan disekitar Krib miring dengan kecepatan di hulu = 0.4 m/dt dan Gambar 4. kontour distribusi kecepatan disekitar krib. Dari contoh gambar yang dihasilkan tersebut dapat diketahui bahwa di sekitar tebing sungai yang diharapkan dilindungi dari bangungan krib tersebut telah mencapai distribusi kecepatan kurang dari 0,25 m/dt.
Pengamatan serupa dilakukan pada hasil running model terhadap ke 45 kasus yang yang direncanakan.
dz z u dy y u dx x u u          z x y (x,y,z) V V+dV v u w (x+dx,y+dy,z+dz) dx x u u    dy y u u    dz z u u    dx dz dy Gambar 1 :   Perubahan kecepatan partikel zat cair arah x, y, z dalam selang waktu dt.
Hasil simulasi numeric terhadap 45 kasus model krib yaitu parameter penting yang menjadi penentu adalah lebar sungai, kecepatan arus, kedalaman air, panjang krib, arah krib itu dipasang dan panjang perlindungan yang dihasilkan dicatat dan kemudian di analisa untuk  dicari hubungan pengaruhnya.

## 3.3. Hubungan Antara Parameter Penting
Untuk mendapatkan  hubungan antara parameter-parameter yang berpengaruh tetapi dalam bentuk general maka dilakukan analisa tak berdimensi. Pada Tabel 1. berikut ini ditunjukkan penurunan analisa non-dimensi dengan cara step wise. 10 20 30 40 50 60 70 10 20 Gambar 2 : Gambar vektor kecepatan disekitar Krib miring dengan kecepatan di hulu = 0.6 m/dt,
Gambar 3 : Gambar vektor kecepatan disekitar Krib miring dengan kecepatan di hulu = 0.4 m/dt,
Gambar 4 : Kontur distribusi kecepatan disekitar disekitar Krib miring dengan kecepatan di hulu = 0.4 m/dt

## Jurnal APLIKASI
ISSN.1907-753X
Dengan merujuk pada Gambar 5, parameter- parameter yang akan dicari hubungannya adalah :
dimana  : L : Panjang tebing yang terlindungi  : Kerapatan air u . : kecepatan rata-rata searah dengan alur sungai v . : kecepatan rata-rata tegak lurus dengan alur sungai h : kedalaman air g : percepatan grafitasi υ : viskositas kinematis P : panjang krib B : lebar sungai Α : sudut yang dibentuk oleh krib dengan garis tebing sungai Tabel 1. Tabel penurunan analisa non- dimensi dengan cara step wise
Dari penurunan diatas diperoleh bilangan tak berdimensi sebagai berikut :
Dari hasil penurunan tersebut diatas kemudian dengan mengkombinasikan antara bilangan non-dimensi menghasilkan sebagai berikut :
Fr adalah Froude Number , = u/  gh ,
Dari penyelesaian terebut didapat hubungan antara bilangan Froude ( Fr. ) dengan Dh/PB, yang dapat ditulis sebagai :
Kemudian dari hasil simulasi terhadap 45 model kasus dicari hubungannya dengan melakukan tabulasi dan menghitung harga bilangan  Froude ( Fr. ) dengan Dh/PB seperti yang disajikan pada Tabel 2 berikut ini. Kemudian antara besarnya nilai Froude ( Fr. ) dan Dh/PB di plot dan kemudian dicari garis regrasinya sebagai yang ditunjuk pada Gambar 6 berikut.
Nampak pada gambar tersebut bahwa dari tiga alternatif model tersebut secara umum dapat dikatakan  Model 1 yaitu Krib dengan sudut  = 90 o merupakan pilihan yang paling baik dibanding dengan model lain yang telah diajukan. Oleh karena itu disarankan dipilih bangunan krib  tegak lurus dengan arah aliran.
Fr PB Dh   B D V Krib P Gambar 5  : Parameter – parameter D  u h g  p b tg   M/L 2 ) L M/L 2 L/T L L/T 2 M/LT L L 1 D u h g  p b 1 u  L/T) L L/T L L/T 2 L 2 /T L L 1 D h g/u 2  u p b 1 h(L) L L L -1 L L L 1 D/h gh/u 2  uh p/h b/h 1  ( D/h , gh/u 2 ,  uh , p/h , b/h , tg  ) gh u uh PB Dh     tan

## 4. KESIMPULAN
Hasil temuan pada penelitian diatas yang menghasilkan suatu hubungan antara bilangan Froude ( Fr. ) dengan Dh/PB, sesuai yang ditunjuk dengan Gambar 6, memberikan kemudahan dalam perencanaan bangunan Krib. Dengan demikian akan sangat menghemat biaya disain karena tahapan studi simulasi tidak lagi perlu dilakukan dalam merencanakan bangunan krib tersebut.

## 5. DAFTAR ACUAN
Suharjoko, 1999,  Studi model numeric Aliran “2D Horizontal” pada kasus pemasangan Krib di Estuari, Thesis S2 Program Studi Teknik Sipil Program Pasca Sarjana UGM.
Suharjoko ,  2000, Perilaku aliran di sekitar bangunan Krib di Estuari , Majalah IPTEK,
ISSN 0953-4098, vol 11, No 2, Mei 2000
Suharjoko ,  2001, Model numeric Perilaku Aliran akibat pemasangan Krib di Sungai , Laporan Penelitian – Lemlit ITS .
Suharjoko , 2002, Panjang Perlindungan Pada Tebing Akibat Adanya Bangunan Krib di Sungai , Prosiding Seminar Nasional Sains dan Teknologi – Lemlit ITS, Surabaya 31 juli 2002
0.515 0.520 0.525 0.530 0.535 0.540 0.545 0.550 0.555 0.000 0.025 0.050 0.075 0.100 0.125 0.150 0.175 0.200 0.225 0.250 0.275 0.300 Dh/PB Fr ( F ru od e N u m b er )  = 90 o  = 45 o  = 135 o Gambar 6  : Hubungan antara Froude (Fr,) dengan nilai
Dh/PB Number   u P B h D Dh/PB Fr K - 1 90 0,30 1,25 6 0,030 45,000 11,25 0,046 0,550 K - 2 90 0,40 1,25 6 0,054 35,000 8,75 0,063 0,548 K - 3 90 0,60 1,25 6 0,120 27,000 6,75 0,108 0,553 K - 4 90 0,80 1,25 6 0,230 25,000 6,25 0,191 0,533 K - 5 90 1,00 1,25 6 0,375 23,000 5,75 0,288 0,521 K - 1 90 0,30 2 6 0,030 60,000 15,00 0,038 0,550 K - 2 90 0,40 2 6 0,054 70,000 17,50 0,079 0,548 K - 3 90 0,60 2 6 0,120 40,000 10,00 0,100 0,553 K - 4 90 0,80 2 6 0,230 35,000 8,75 0,168 0,533 K - 5 90 1,00 2 6 0,375 35,000 8,75 0,274 0,521 K - 1 90 0,30 3,25 6 0,030 80,000 20,00 0,031 0,550 K - 2 90 0,40 3,25 6 0,054 80,000 20,00 0,056 0,548 K - 3 90 0,60 3,25 6 0,120 70,000 17,50 0,108 0,553 K - 4 90 0,80 3,25 6 0,230 55,000 13,75 0,162 0,533 K - 5 90 1,00 3,25 6 0,375 50,000 12,50 0,240 0,521 K - 1 135 0,30 1,25 6 0,030 48,000 12,00 0,049 0,550 K - 2 135 0,40 1,25 6 0,054 33,000 8,25 0,060 0,548 K - 3 135 0,60 1,25 6 0,120 18,000 4,50 0,072 0,553 K - 4 135 0,80 1,25 6 0,230 12,000 3,00 0,092 0,533
K - 5 135 1,00 1,25 6 0,375 15,000 3,75 0,188 0,521 K - 1 135 0,30 2 6 0,030 55,000 13,75 0,035 0,550 K - 2 135 0,40 2 6 0,054 58,000 14,50 0,066 0,548 K - 3 135 0,60 2 6 0,120 45,000 11,25 0,113 0,553 K - 4 135 0,80 2 6 0,230 28,000 7,00 0,134 0,533 K - 5 135 1,00 2 6 0,375 25,000 6,25 0,195 0,521 K - 1 135 0,30 3,25 6 0,030 K - 2 135 0,40 3,25 6 0,054 K - 3 135 0,60 3,25 6 0,120 K - 4 135 0,80 3,25 6 0,230 K - 5 135 1,00 3,25 6 0,375 K - 1 45 0,30 1,25 6 0,030 34,000 8,50 0,034 0,550 K - 2 45 0,40 1,25 6 0,054 28,000 7,00 0,051 0,548 K - 3 45 0,60 1,25 6 0,120 23,000 5,75 0,092 0,553 K - 4 45 0,80 1,25 6 0,230 17,000 4,25 0,130 0,533 K - 5 45 1,00 1,25 6 0,375 16,000 4,00 0,200 0,521
K - 1 45 0,30 2 6 0,030 40,000 10,00 0,025 0,550 K - 2 45 0,40 2 6 0,054 35,000 8,75 0,040 0,548 K - 3 45 0,60 2 6 0,120 27,000 6,75 0,068 0,553
K - 4 45 0,80 2 6 0,230 24,000 6,00 0,115 0,533 K - 5 45 1,00 2 6 0,375 23,000 5,75 0,180 0,521 K - 1 45 0,30 3,25 6 0,030 90,000 22,50 K - 2 45 0,40 3,25 6 0,054 90,000 22,50 K - 3 45 0,60 3,25 6 0,120 99,000 24,75 K - 4 45 0,80 3,25 6 0,230 80,000 20,00 K - 5 45 1,00 3,25 6 0,375 90,000 22,50 M odel - 3 Var-1 Var-2 Var-3 M odel - 2 Var-2 Var-3 Var-1 Var-1 Model Krib Variasi Krib Tabel Hubungan antara Dh/PB dengan bilangan Froude. M odel - 1 Kasus Data simulasi dan hasil simulasi Var-2 Var-3 Tabel 2