## PEMODELAN NUMERIK KEKUATAN PINTU AIR BAJA

## THE NUMERICAL MODEL OF THE STEEL FLOOD GATE STRENGTH
Fitriansyah *1 , Miftahul Iman 2 , Aminullah 3
1 Dosen, Teknik Sipil, Fakultas Teknik, Universitas Achmad Yani Banjarmasin
2 Dosen, Teknik Sipil, Fakultas Teknik, Universitas Borneo Tarakan
3 Dosen, Teknik Sipil, Fakultas Teknik, Universitas Nahdlatul Ulama Kalimantan Selatan Korespondensi: fitriansyah@uay.ac.id

## ABSTRAK
Bangunan  intake  untuk  Embung  Sei  Bubu  telah  dirancang  untuk  menjadi  saluran  terbuka  dengan  gerbang banjir.  Gerbang  banjir  telah  dirancang  dari  baja  dengan  lebar  bersih  1,35  m  dan  ketinggian  bukaan  pintu 0,0076  m.  Korosi  merupakan  salah  satu  ancaman  berbahaya  bagi  kekuatan  dan  daya  tahan  gerbang  banjir. Penelitian  ini  secara  numerik  memodelkan  gerbang  banjir  yang  diserang  oleh  korosi  lubang.  Korosi  lubang telah  dimodelkan  dalam  beberapa  lubang  kecil  secara  acak  didistribusikan  di  permukaan  pintu  banjir.  Efek dari  level  air  telah  dimasukkan.  Pemodelan  numerik  telah  dilakukan  dalam  metode  elemen  hingga memanfaatkan  program  komputer  seperti  Abaqus.  Hasil  penelitian  menunjukkan  ada  penurunan  kapasitas gerbang banjir baja karena lubang meningkat. Penurunan kapasitas tegangan telah ditunjukkan oleh tegangan yang terjadi di sekitar korosi lubang. Penurunan tegangan terjadi dengan perubahan persentase area distribusi korosi  pitting  sebesar  85,47%  (70,62  MPa),  91,95%  (127,42  MPa)  dan  82%  (57,04  MPa)  untuk  2/3H. Perpindahan maksimum terjadi di IFM, sub model 2 dalam 16,39 mm.
Kata Kunci: Baja, Korosi, Lubang, Pintu air, Tegangan

## ABSTRACT
The intake building for Embung Sei Bubu had been designed to be an open channel with the flood gate. The flood  gate  had  been  designed  from  steel  with  a  net  width  of  1.35  m  and  the  height  of  the  door  opening  is 0.0076  m. The  corrosion  is  one  of  the hazardous  threat  to  the  strengthness  and durability  of  the  flood  gate. This research numerically models flood gate that was attacked by pitting corrosion. The pitting corrosion had been modelled in several small holes randomly were distributed on the surface of the flood gate, precisely on the  surface  of  the  water.    The  numerical  modeling  had  been  performed  in  finite  element  method  utilized computer programs such Abaqus. The results showed there was a reduction in the capacity of the steel flood gate due the hole increasing. The reduction in stress capacity had been indicated by the stress concentration that was occured around the pitting corrosion. The stress reduction occured with the change in the percentage of pitting corrosion distribution area of 10% (225 MPa), 20% (175 MPa) and 30% (120 MPa).
Keywords: Corrosion, Flood gate. Hole, Steel, Stress

## MEDIA ILMIAH TEKNIK SIPIL
Volume 9 Nomor 2 Juni 2021 Hal. 116-123

## PENDAHULUAN
Embung  Sei  Bubu  telah  dirancang  untuk memenuhi  kebutuhan  air  dasar  untuk  kebutuhan domestik,  serta  air  irigasi  untuk  meningkatkan intensitas penanaman. Daerah irigasi yang potensial  untuk  dilayani  adalah  Daerah  Irigasi  Sei Bubu dengan total area pengembangan 200 Ha. Sei Bubu adalah salah satu sungai di Kabupaten Tanah Bumbu.  Sungai  ini  membentang  dari  utara  ke selatan dengan panjang hingga muara adalah  1,5 km.  Bangunan  intake  untuk  reservoir  Sei  Bubu telah dirancang di saluran terbuka dengan gerbang banjir.

## TINJAUAN PUSTAKA
Dimensi gerbang banjir ditentukan berdasarkan  ketinggian  air  di  reservoir  dan persyaratan  maksimum  untuk  air  irigasi  dan  air baku.
H g C Q A . . 2 . =
(1)
dengan:
A : Luas pintu air
Q : Debit
C : Koefisien debit untuk saluran terbuka (= 0.60)
g : Percepatan gravitasi
H : Tinggi kehilangan energy pintu air
Alex  dkk.  (2017)  menyajikan  pemanfaatan model numerik untuk menentukan gaya hidrodinamik  yang  bekerja  pada  pintu  air  gerobak datar,  yang  digunakan  dalam  karya  pengalihan PLTA.  Triatmodjo  (1990)  menjelaskan  bahwa tekanan hidrostatik yang dilakukan pada penampang  pintu  air  dapat  dihitung  dengan menggunakan Persamaan 2.
gh P  = (2)
dengan P adalah tekanan hidrostatik (N/m 2 ),  adalah  massa  jenis  air  (kg/m 3 ), g adalah percepatan  gravitasi  (m/s 2 ),  dan h adalah ketinggian air (m).
James  (2015)  Banjir  yang  terjadi  di  Jakarta sudah  menjadi  hal  biasa  bagi  masyarakat  yang tinggal  di  sana.  Hampir  setiap  tahun  banjir menggenangi  jalan  utama  dan  daerah  perumahan.
Fenomena  ini  tentu  menjadi  perhatian  banyak orang  karena  memengaruhi  kehidupan  ekonomi masyarakat sekitar. Salah satu penyebabnya adalah berkurangnya kapasitas saluran sungai Ciliwung.
Daniel  (2010)  Pemodelan  pergerakan  air banjir  dapat  dicapai  dengan  menggunakan  banyak metode  berbeda  dengan  berbagai  tingkat  detail fisik.  Model  numerik  yang  menggunakan  metode perutean sederhana atau versi persamaan Navier-Stokes yang disederhanakan dapat digunakan  untuk  meningkatkan  kesiapsiagaan banjir publik.
Miftahul  dkk.  (2018)  Korosi  pitting  sering terdapat berupa lubang kecil pada baja tubular dari sebuah  struktur  platform  ketika  lapisan  pelindung korosi  (coating)  telah  rusak.  Sebuah  lubang tunggal  pada  elemen  tekan  yang  langsing  dapat menyebabkan pengurangan yang signifikan terhadap  kapasitas  tekuk  batang.  Penelitian  ini dilakukan  untuk  mempelajari  secara  numerik pengaruh  lubang  bundar  terhadap  kapasitas  tekuk baja tubular langsing.

## METODE
Penelitian ini adalah studi numerik berdasarkan metode elemen hingga yang memodelkan  saluran  irigasi  (Gambar.1).  Model dalam  penelitian  ini  terdiri  dari  dua  model,  model pintu  air  tanpa  korosi  lubang  dan  model  pintu  air dengan  korosi  lubang.  Korosi  pitting  pada penelitian  ini  dimodelkan  sebagai  lubang-lubang kecil  yang  terdistribusi  pada  permukaan pelat  baja gerbang  banjir.  Model  pintu  banjir  dengan  lubang telah  dikelompokkan  menjadi  3  sub  model berturut-turut,  yaitu:  sub  model  1,  sub  model  2, dan sub model 3. Distribusi korosi lubang di mana terletak  di  bagian  atas,  tepi  dan  tengah  pelat. Jumlah  lubang  telah  ditentukan  sebanyak  5  buah per distribusi.

## Perfect model (PFM)
Pemodelan  numerik  pintu  banjir  (Gbr.2) menggunakan  kombinasi  elemen  heksahedron padat  dan  tetrahedron  yang  tersedia  di  Abaqus. Model  ini  terdiri  dari  3  bagian,  yaitu  profil  UNP 120 X 55 x 3,8;
Profil  L75  x  75  x  7  dan  pelat  tebal  10  mm terbuat  dari  baja.  Tiga  bagian  dihubungkan  oleh kendala  dasi.  Kondisi  batas  diidealkan  dengan dukungan sederhana.
Beban  pada  model  didefinisikan  sebagai tindakan tekanan hidrostatik pada permukaan pelat baja gerbang banjir. Tekanan hidrostatik diperoleh dengan  rumus  hidrostatik  (Persamaan  2)  dimana dilakukan pada 0,50 H dan 2 / 3H (H = ketinggian pintu  air)  pada  permukaan  pelat  baja  model gerbang banjir.

## Gambar 2. Flood gate model
Gambar. 2 menunjukkan bahwa model pintu banjir dikonfigurasi dalam model padat. Model ini menggunakan elemen hexahedron dengan perkiraan  ukuran  global  adalah  25  (Gbr.  3).  Pelat baja diberi penyangga bersama dengan perpindahan  arah  vertikal  (U3)  tidak  sama  dengan nol. Ini dilakukan karena pelat baja dimungkinkan bergerak  secara  vertikal.  Pelat  baja  ditempatkan pada balok beton berukuran 10 cm x 10 cm x 135 cm sebagai model pelat beton untuk dasar saluran. Kekuatan  tekan  beton  ditentukan  (f)  20  MPa  dan rasio Poisson (  ) 0,2.

## Gambar 3. Flood gate steel plate model
Balok beton telah dibatasi dengan pelat baja yang digunakan  pembatas  di  Abaqus. Oleh  karena itu  pelat  baja  dan  balok  beton  dapat  merespons beban  tekanan  secara  bersamaan.  Permukaan  air telah  dievaluasi  berdasarkan  ketinggian  pelat  baja (H  =  1350  mm).  Ketinggian  air  telah  ditentukan sebelumnya:  1  /  3H;  1  /  2H;  2  /  3H,  di  mana  H adalah ketinggian pelat baja gerbang banjir.

## Gambar 4. Hydrostatic pressure load
Gambar.  4  menunjukkan  bahwa  beban bertindak  sebagai  beban  tekanan  hidrostatik. Tekanan hidrostatik telah tergantung pada ketinggian  air.  Dalam  penelitian  ini  ketinggian  air dikelompokkan menjadi 3 jenis. Ada 0,33H, 0,5H, dan 0,67H, di mana H adalah ketinggian pintu air.

## Imperfect Model
Model  tidak  sempurna  terdiri  dari  tiga sub  model:  sub  model  1,  sub  model  2,  dan  sub model  3.  Model-model  tersebut  telah  dibedakan berdasarkan lokasi korosi pitting, mereka berada di kedua  sisi  pelat  baja,  di  tengah-tengah  pelat  baja, dan di pelat baja atas selaras vertikal. Setiap model memiliki ketinggian air yang berbeda, 0,33 H, 0,50 H, dan 0,67 H di mana H adalah ketinggian air.

## Sub Model 1
Sub model 1 memiliki korosi pitting di mana terletak  di  kedua  sisi  pelat  baja  (Gbr.  5).  Model telah  dibangun  dalam  tetrahedron  dan  elemen wedge  dalam  diskritisasi  25  elemen  ukuran  dalam idealisasi solid 3D.

## (a) (b)                    (c)

## Gambar 5. The meshing of the sub model
Meshing  di  sekitar  lubang ditentukan dalam elemen  tetrahedron  (Gambar  6).  Lubang  telah diidealkan  oleh  kecil  (diameter  =  5  mm,  dan kedalaman  lubang  3,5  mm).  Itu  artinya  lubang tidak bisa menembus pelat baja.

## Gambar 6. Meshing around the hole

## HASIL DAN PEMBAHASAN
Hasil  penelitian  telah  dibagi  menjadi  dua model berturut-turut yaitu model pintu banjir tanpa efek  korosi  pitting  (model  pintu  banjir  sempurna, PFM)  dan  model  pintu  air  dengan  efek  korosi lubang  (model  imperfect  flood  gate  model,  IFM). Model  IFM  terdiri  dari  tiga  sub-model  yang dibedakan  oleh  efek  ketinggian  airembahasan berisi  tentang  hasil  penelitian  yang  telah dibandingkan  dengan  data-data  yang  diperoleh peneliti sebelumnya.

## Perfect Model (PFM)
Output  yang  diperoleh  dari  PFM  adalah stres  dan  perpindahan.  Tegangan  yang  dievaluasi adalah  kriteria  stres  maksimum  (absolut)  dan  von mises. Evaluasi perpindahan telah difokuskan pada perpindahan  tegak  lurus  (U2)  dari  pelat  baja gerbang banjir.
Gambar. 7 menunjukkan kontur tegangan karena tekanan hidrostatik dilakukan pada permukaan  pelat  model  gerbang  banjir  tanpa  efek korosi  pitting  (PFM).  Tegangan  yang  terjadi  pada pelat berkisar antara 0,02-20,5 MPa dalam 2 / 3H, di  mana  H  adalah  kedalaman  pelat  baja  gerbang banjir.
Tabel  1  menunjukkan  bahwa  tegangan rata-rata  dan  perpindahan  terjadi  dalam  model gerbang  banjir  sempurna  (tanpa  penurunan kualitas).  Tabel  menunjukkan  bahwa  tegangan maksimum terjadi ketika level air dalam 1 / 2H dan 2/3H  (10,67  MPa).  Sebaliknya  mereka  memiliki perpindahan terendah dalam 2 arah (0,41 mm).
Table 1. Tegangan dan perpindahan (PFM) Tegangan rerata von mises (MPa) Perpindahan (mm) 0,33H 4,25 -5,98 0,50H 10,27 0,41 0,67H 10,27 0,41
Tegangan  terendah  pada  model  sempurna (PFM)  terjadi  ketika  permukaan  air  dalam  1  /  3H (4,25  MPa),  di  mana  H  adalah  ketinggian  pelat baja.  Perpindahan  maksimum  telah  ditunjukkan pada  Tabel  1  ketika  permukaan  air  pada  1  /  3H (-5,98).  Perpindahan  telah  dipertimbangkan  dalam 2 arah (U2).

## Gambar 7. Stress contour on PFM
Kontur  tegangan  (von  mises)  pada  pelat baja  telah  dievaluasi  ketika  ketinggian  air  pada 0,3H adalah 4,22 MPa (Gambar. 7 (a)). Gambar. 7 (b)  dan  Gambar.  7  (c)  menunjukkan  bahwa  pelat baja ketika tingkat air pada 0,5H dan 0,67H adalah 10,67 MPa.

## Gambar 8. Stress contour (PFM, pressure)
Kontur tegangan (tekanan) telah ditunjukkan pada  Gambar.  8  dengan  ketinggian  air  yang
berbeda.  Tegangan  maksimum  ditunjukkan  pada 0,95  MPa  ketika  level  air  dalam  1  /  3H.  Kontur pemindahan (U2) telah disajikan pada Gambar.9.
Gambar 9. Displacement contour (U 2 , PFM)

## Imperfect Model (IFM)
Model  pintu  banjir  yang  tidak  sempurna terdiri  dari  tiga  submodel  yang  dibedakan berdasarkan  tingkat  distribusi  korosi  lubang.  Sub model  pertama  adalah  pelat  di  mana  distribusi korosi  lubang  ditemukan  di  kedua  tepi  pelat  pintu banjir  (Sub  model  1).  Sub  model  kedua  adalah pelat  yang  memiliki  distribusi  korosi  pitting  di tengah pelat pintu banjir (Sub model 2). Sub model ketiga  dari  pelat  yang  memiliki  distribusi  korosi lubang di sisi atas pelat pintu banjir (Sub model 3).

## Sub model 1
Sub  model  1  memiliki  tegangan  luntur  di sekitar  lubang  (von  mises)  mulai  dari  4,22  MPa hingga 10,27 MPa (Gambar. 11).

## Gambar 10. Stress contour on IFM
Gambar  10  menunjukkan  bahwa  kontur tegangan  (von  mises)  berbeda  pada  tiga  kondisi keadaan permukaan air. Gambar. 10 (a) menunjukkan bahwa kontur tegangan terjadi ketika ketinggian  air  pada  0,3H.  Gambar.  10  (b)  dan Gambar.  10  (c)  menunjukkan  kontur  tegangan dengan ketinggian air yang sesuai 0,5H dan 0,67H.

## Gambar 11. Stress concentration at IFM 1
Gambar.  11  menunjukkan  konsentrasi  tegangan yang terjadi pada mengelilingi lubang di kedua sisi pelat baja model pintu banjir.
Gambar 12. IFM 1 stress (pressure)
Gambar  12  menunjukkan  bahwa  kontur  tegangan (tekanan)  berbeda  tiga  kondisi  keadaan  ketinggian air.  Gambar.  12(a)  menunjukkan  bahwa  kontur tegangan  (tekanan)  terjadi  ketika  ketinggian  air pada  0,3H.  Gambar.  12(b)  dan  Gambar.  12(c) menunjukkan  kontur  tegangan  dengan  ketinggian air yang sesuai 0,5H dan 0,67H.
Gambar 13. IFM 1 stress (pressure)
Gambar 13 menunjukkan bahwa kontur perpindahan  (von  mises)  berbeda  tiga  kondisi tingkat permukaan air. Gambar. 12 (a) menunjukkan bahwa kontur perpindahan (tekanan) terjadi ketika ketinggian air pada 0,3H. Gambar. 12 (b)  dan  Gambar.  12  (c)  menunjukkan  kontur perpindahan  dengan  ketinggian  air  yang  sesuai 0,5H dan 0,67H.
Table 2 . Tegangan dan perpindahan(IFM 1)
Tegangan  rerata vm (MPa) Perpindahan (mm) 0,33H 50,04 2,11 0,50H 40,41 1,40 0,67H 70,62 2,46

## Sub model 2
Sub  model  2  memiliki  tegangan  di  sekitar  lubang mulai  dari  145  hingga  175  MPa  (von  mises). Konsentrasi  tegangan  yang  mengelilingi  lubang telah dievaluasi dalam 35 MPa hingga 75 MPa.
Gambar 14. Stress contour on IFM 2 Von Mises
Gambar. 14 menunjukkan bahwa kontur stres (von mises)  berbeda  tiga  kondisi  keadaan  permukaan air.  Gambar.  14  (a)  menunjukkan  bahwa  kontur tegangan  terjadi  ketika  ketinggian  air  pada  0,3H. Gambar.  14  (b)  dan  Gambar  14  (c)  menunjukkan kontur tegangan dengan ketinggian air yang sesuai 0,5H dan 0,67H.
Gambar 15. Stress concentration at IFM 2
Gambar.  15  menunjukkan  konsentrasi  tegangan yang terjadi pada mengelilingi lubang di kedua sisi pelat baja model gerbang banjir.
Gambar 16. Stress contour (pressure, submodel 2)
Gambar.  16  menunjukkan  bahwa  kontur  tegangan (tekanan)  berbeda  tiga  kondisi  keadaan  ketinggian air.  Gambar.  16  (a)  menunjukkan  bahwa  kontur tegangan  (tekanan)  terjadi  ketika  ketinggian  air pada  0,3H.  Gambar.  16  (b)  dan  Gambar.  16  (c) menunjukkan  kontur  tegangan  dengan  ketinggian air yang sesuai 0,5H dan 0,67H.
Table 3 Stress and Displacement (IFM 2) Average stress vm (MPa) Displacement (mm) 0,33H 46,85 5,50 0,50H 97,57 11,99 0,67H 127,42 16,39
Tabel 3. menunjukkan bahwa tegangan maksimum dan  perpindahan  pada  submodel  2  ketika ketinggian  air  pada  2  /  3H.  Tegangan  minimum dan  dispelement  terjadi  ketika  level  air  pada  1  / 3H.
Gambar 17. Displacement contour (U2) Gambar. 17 menunjukkan bahwa kontur
perpindahan  (von  mises)  berbeda  tiga  kondisi keadaan ketinggian air. Gambar 17 (a) menunjukkan bahwa kontur perpindahan (tekanan) terjadi ketika ketinggian air pada 0,3H. Gambar. 17 (b)  dan  Gambar.  17  (c)  menunjukkan  kontur perpindahan  dengan  ketinggian  air  yang  sesuai 0,5H dan 0,67H.

## Sub model 3
Sub  model  3  memiliki  tegangan  di  sekitar  lubang mulai  dari  75  MPa  hingga  125  MPa  (von  mises, Gambar.8 (a)).

## Gambar 18. Stress contour IFM Von Mises
Gambar. 18 menunjukkan bahwa kontur stres (von mises)  berbeda  tiga  kondisi  keadaan  permukaan air.  Gambar  18  (a)  menunjukkan  bahwa  kontur tegangan terjadi ketika ketinggian air pada 0,3 jam. Gambar.  18  (b)  dan  Gambar.  18  (c)  menunjukkan kontur tegangan dengan ketinggian air yang sesuai 0,5H dan 0,67H.

## Gambar 19. Stress contour at IFM, submodel 3
Gambar.  19  menunjukkan  konsentrasi  tegangan yang terjadi pada mengelilingi lubang di kedua sisi pelat baja model pintu air.

## Gambar 20. Stress contour (pressure) sub model 3
Gambar  20  menunjukkan  bahwa  kontur  tegangan (tekanan)  berbeda  pada  tiga  kondisi  keadaan ketinggian air. Gambar 20 (a) menunjukkan bahwa kontur tegangan (tekanan) terjadi ketika ketinggian air  pada  0,3H.  Gambar  20  (b) dan  Gambar. 20  (c)
menunjukkan  kontur  tegangan  dengan  ketinggian air yang sesuai 0,5H dan 0,67H.
Gambar 21 . Displacement contour (U 2 , IFM 3)
Gambar. 21 menunjukkan bahwa kontur perpindahan  (von  mises)  berbeda  tiga  kondisi keadaan ketinggian air. Gambar. 21 (a) menunjukkan bahwa kontur perpindahan (tekanan) terjadi ketika ketinggian air pada 0,3H. Gambar. 21 (b)  dan  Gambar.  21  (c)  menunjukkan  kontur perpindahan  dengan  ketinggian  air  yang  sesuai 0,5H dan 0,67H.
Table 3. Stress and Displacement (IFM 3)
Average stress vm (MPa) Displacement (mm) 0,33H 38,25 2,25 0,50H 39,73 1,62 0,67H 57,04 1,04
Tabel  3  menunjukkan  tegangan  rata-rata  (von mises)  yang  terjadi  pada  permukaan  pelat  dengan ketinggian  air  yang  berbeda.  Tabel  menunjukkan bahwa sub model 3 dengan ketinggian air dalam 2 /  3H  adalah  tegangan  maksimum  (57,04  MPa) tetapi memiliki perpindahan terendah (1,04 mm).
Gambar 22 . Stress – water elevation curve
Gambar.  22  menunjukkan  kurva  elevasi air-stres  dari  model  oleh  efek  elevasi  air.
Penurunan  tegangan  terjadi  dengan  perubahan persentase  area  distribusi  korosi  pitting  sebesar 85,47%  (70,62  MPa),  91,95%  (127,42  MPa)  dan 82%  (57,04  MPa)  dalam  2/3H  .  Ini  berarti keberadaan cacat pada permukaan pelat baja secara signifikan  mempengaruhi  kinerja  dan  daya  tahan gerbang banjir.

## PENUTUP

## Kesimpulan
Hasil  dari  hasil  ini  telah  disimpulkan sebagai berikut:
1.  Hasil  penelitian  menunjukkan  ada  penurunan kapasitas  pintu  air  baja  karena  lubang meningkat.
2.  Penurunan kapasitas tegangan telah ditunjukkan oleh  konsentrasi  tegangan  yang  terjadi  di sekitar lubang.
3.  Penurunan  tegangan  terjadi  dengan  perubahan persentase  area  distribusi  korosi  pitting  sebesar 85,47%  (70,62  MPa),  91,95%  (127,42  MPa) dan 82% (57,04 MPa) dalam 2 / 3H .
4.  Perpindahan  maksimum  terjadi  di  IFM,  sub model 2 dalam 16,39 mm.
5.  Lokasi korosi pengutipan di tengah pelat (IFM,
sub  model  2)  lebih  berbahaya  daripada  lokasi lain.  kesimpulan  dapat  dibuat  dalam  bentuk paragraf atau dapat juga diuraikan per poin.

## DAFTAR PUSTAKA
Alex,  C.  R,  Luiz,  F.  M,  Jairo,  A.  M,  (2017), Numerical  Simulation  of  Drainage  in  the Floodgate of a Diversion Tunnel of
Hydroelectric  Plants,  Revista  Espacios,  Vol 38, Spanyol Daniel,  W.  G,  (2010),  Application  of  Numerical Models for Improvement of Flood Preparedness,  Nama  Jurnal,  Nomor  terbit, Iowa James, Z., (2015), Efektifitas Hidraulik
Penambahan  Pintu  Air  Melalui  Uji  Model Fisik  3D  dan  Model  Numerik  1D  (Studi Kasus:  Pintu  Air  Manggarai), Jurnal  Teknik Hidraulik , Vol 6 No. 01, Indonesia
Miftahul,  I,  Bambang,  S,  Hrc  Priyosulistyo  and Muslikh,  (2018),  Applied  Mechanic  and Materials,  Trans  Tech  Publications,  vol 881 UK
Triatmodjo,  B,  (1991), HIDRAULIKA  1 ,  1,  Beta Offset Yogyakarta