## HUBUNGAN ANTARA DIMENSI KOLAM OLAK DAN EFEKTIFITAS PEREDAM ENERGI DI HILIR PINTU AIR
Yuliman Ziliwu

## Abstrak
Loncat  air  di  hilir  sluice  gate  sering  menimbulkan  masalah  yaitu  terjadinya local scouring.  Salah  satu  untuk  mengatasi  masalah  tersebut    adalah  menambah  kekasaran apron menggunakan endensill di hilir sluice gate. Tujuan dari penelitian ini adalah untuk mendapatkan hubungan antara dimensi kolam olak dan efektifitas peredam energi di hilir pintu  air.  Pelaksanaan  penelitian  ini  dimulai  dari  kalibrasi  debit  aliran  yang  melalui standard tilting flume dan didapatkan nilai debit yaitu Q= 4830 cm 3 /det yang besarnya disesuaikan dengan kapasitas flume pada saat sluice gate di pasang dengan bukaan pintu terkecil  agar  tidak  meluap.  Dari  hasil  penelitian  :  1.  nilai  frude  number  untuk  model dengan  panjang  kolam  olak  yang  sama  adanya  penambahan  kedalaman  kolam  olak mengakibatkan  LS  semakin  pendek;  2.Untuk  suatu  nilai  Froude  number  yang  sama, panjang  loncat  hidraulik  yang  terjadi  semakin  berkurang  dengan  adanya  penambahan kedalaman  kolam  olak,  penambahan  tailwater  depth  dan  penambahan  tinggi  endsill; 3.Penambahan  kedalaman  dan  panjang  kolam  olak,  penambahan  tailwater  depth  dan penambahan  tinggi  endsill  dapat  memperbesar  pengaruh  peredaman  energi; 4.Penambahan endsill pada ujung hilir kolam olak dapat mengurangi terjadinya osilasi yang menyebabkan terjadinya kerusakan pada tanggul – tanggul dari tanah dan riprap.

## Kata kunci: Loncat air, kolam olak, peredam nenergi

## 1. PENDAHULUAN
Terjadinya loncat hidraulik di hilir sluice gate sering  menimbulkan  terjadinya local scouring.  Berbagai  cara  untuk mengatasi  masalah  tersebut  diantaranya adalah  menambah  kekasaran  apron menggunakan endsill di hilir sluice gate, namun  terkadang  menimbulkan  masalah lain,  yaitu  terjadinya second jump.
Langkah – langkah  tambahan,  seperti misalnya
menurunkan dan memperpanjang  tinggi endsill yang bertujuan  menambah  kekasaran  apron diperlukan agar dapat menjamin perencanaan energi secara memadai. Penelitian ini adalah untuk mendapatkan hubungan antara dimensi kolam olak dan efektifitas  peredaman  energi.  Beberapa
peneliti  telah melakukan  penetitian tentang  pengaruh  letak endsill di  hilir pintu  air  S.WU.  dan  N.  Rajaratnam (1995)  menyelidiki  tentang  pengaruh letak endsill di hilir pintu air pada aliran submerged.  Ranga  Raju  dan  Visaradia (1979)  menyelidiki  tentang  sifat – sifat debit  yang  terjadi  pada  pintu  penguras dengan  kondisi  aliran  bebas.  Mossa  M., Petrillo  A.  dan  chanson  H.  (2002), melakukan penelitian tentang peredaman energi  loncat  hidraulik  dalam kolam olak.

## 2. TINJAUAN PUSTAKA
Ada 3 macam pola loncat hidraulik yang mempunyai kemungkinan terbentuk pada  daerah  hilir sluice  gate (Chow,  V. T., 1959), yaitu :
1. Kasus  1,  apabila  kedalam  hilir (tailwater  depth)  y 3 sama  dengan kedalaman y 2 (sequent depth)  yakni lanjutan  dari  y 1 .  Nilai – nilai  Fr 1 ,  y 1 dan  y 3 (  =y 2 )  akan  memenuhi persamaan 1 :
( ) 1 8 1 2 1 2 1 1 2 − + = Fr y y …(1)
2. Kasus  2,  apabila  kedalam  hilir  y 3 < y 2 . Hal    ini  berarti  kedalaman  hilir pada  kasus  1  berkurang.  Akibatnya, loncat hidraulik akan mundur ke hilir ke  suatu  titik  dimana  persamaan  1 dipenuhi kembali.
3. Kasus  3,  apabila  kedalam  hilir y 3 > y 2 .  Hal    ini  berarti  kedalaman  hilir pada  kasus  1  bertambah  besar. Akibatnya,  loncatan  di  dorong  ke arah  hulu  dan  berubah  menjadi loncatan terendam.

## Persamaan  Momentum  di  Mulut

## Sluice gate
Mossa  M.,  Petrillo  A.  dan  Chanson  H. (2002) menunjukkan bahwa :
Pola aliran =       0 0 1 0 , , y H y y Fr f …(2)
Persamaan dasar aliran untuk penyelesaian  analitis  pada  penelitian tersebut  di  atas  diaplikasikan  dalam penelitian ini.
Gambar 1. Pers. Momentum di mulut sluice gate
Persamaan  momentum  di  mulut sluice
gate (Gambar 1) adalah :
F 0 + F h – F 1 = pQ (U 1 – U 0 ) …(3)
Dengan :
F 0 = ½  x y 0 2 B F h =  H B H y       + 2 0
F 1 = ½  x y 1 2 B

## …(4)
Substitusi  persamaan  (4)  ke  persamaan (3)  dan  dengan  penjabaran  seperlunya, didapatkan :
Fr 0 2 =
                    −         −     −     −     0 1 0 2 0 2 0 1 1 1 2 1 2 y y y H y H y y …(5)

## Kekasaran Saluran
Energi  aliran  superkritis  yang  masuk dalam  kolam  olak  akan  direndam  oleh kekasaran saluran.
Gambar 2. Gaya – gaya yang bekerja pada Loncat Hidraulik
Persamaan momentum pias 0-2 adalah :
 A 0 z 0 -  A  x  h f -  A 2 z 2 =  pQ(U 2 – U 0 ) …(6)
Bila  h f =  Ls  x f S dan f S = 3 4 2 2 2 R A n Q

## …(7)
Substitusi  persamaan  (7)  ke  dalam persamaan  (5)  dan  dengan  penjabaran seperlunya, didapatkan :
Ls = ( ) 3 4 2 2 2 2 0 0 2 2 2 2 2 0 2
R A n Q xy y y y Bg Q y y B     − + − …(8) Persamaan
Momentum Loncat

## Hidraulik
Bila  La adalah panjang  kolam olak,  L  = Ls  +  Lj,  dimana  Ls  adalah  jarak  antara pintu  air  dengan  awal  loncat  hidraulik dan  Lj  adalah  panjang  loncat  hidraulik, maka  posisi  loncat  hidraulik  di  dalam kolam olak, yaitu :
a. Apabila  L/La<1,  maka  persamaan yang berlaku adalah persamaan 1.
b. Apabila  L/La ≥ 1,  maka  fenomena loncat  hidraulik  yang  terjadi  pada peninggian saluran secara mendadak.
Gambar 3. Posisi Loncat Hidraulik dalam Kolam Olak
Persamaan  momentum  di  ujung hilir konlam olak (Gambar 3) adalah :
F 1 – F 2 – F 3 = pQ (U 2 – U 1 ) …(9)
Dengan
F 1 = ½  By 1 2
F 2 = ½ 2 y By 2 2
F 3 = -pU 1 2 BH …(10)
Substutusi  persamaan  (10)  ke  dalam persamaan  (9)  dan  dengan  penjabaran seperlunya, didapatkan :
) 2 1 ( 2 ) 2 1 ( ) 2 1 ( 2 2 2 1 2 2 2 1 1 2 Fr H Fr Fr Fr y y + + + + =
…(11)

## Panjang Loncat Hidraulik
Rajaratnam dan Biro
Reklamasi Amerika  Serikat  (USBR)  mengusulkan hubungan  panjang  loncat  hidraulik  pada
lantai  mendatar  (Chow,  V.T.,  1959), yaitu :
Lj = A(y 2 – y 1 ) …(12)

## 3. METODE PENELITIAN
Pelaksanaan penelitian ini dimulai dengan  kalibrasi debit aliran  yang melalui Standart  Tilting  Flume dan didapatkan  satu  nilai  debit,  yaitu Q  = 4830 cm 3 /det, yang besarnya disesuaikan dengan kapasitas flume pada saat sluice gate dipasang dengan bukaan pintu terkecil agar air tidak meluap.
Tinggi  bukaan  pintu  air  y 0 ,  masing – masing  bukaan 1.25  cm,  1.5  cm dan 2 cm,  agar  didapatkan  angka Froude dengan kisaran 1,0<Fr 1 ,4,5.
Pemberian  nama  pada  model  dengan ketentuan :
a. G1 (model  dengan  bukaan  pintu  2 cm), G2 (model dengan bukaan pintu 1,5  cm)  dan G3 (model  dengan bukaan pintu 1.25 cm) b. K1 (model dengan kedalaman kolam 1 cm), K2 (model dengan kedalaman kolam 2 cm) dan K3 (model dengan kedalaman kolam 3 cm) c. E1 (model  dengan  tinggi endsill 1 cm), E2 (model dengan tinggi endsill 2 cm)
d. D1 (model  dengan  jarak  apron  120
cm), dan  seterusnya  sampai D12
(model dengan jarak apron 10 cm)
Langkah kedua adalah mengatur kedalaman muka air hilir (y 3 ).
Setelah  debit  konstan,  pada  tahapan  ini dilakukan pengukuran ragam kedalaman aliran  yang  terbentuk  di  hilir  pintu  air dengan menggunakan point gauge antara lain kedalaman air di hulu pintu air y hulu pada  daerah vena  contracta y c ,  awal loncat hidraulik  y 1 , tinggi puncak loncat hidraulik  y p ,  kedalaman  aliran  di  atas endsill y edsill dan  kedalaman  aliran  di hilir  loncat  hidraulik  y 3 .  Sejalan  dengan pelaksanaan  tersebut  juga  dilakukan pengukuran jarak antara pintu air dengan awal  loncatan  Ls,  panjang  loncatan  Lj, jarak  antara  awal  loncatan  sampai dengan  puncak  loncatan  Lp  dan  jarak antara awal loncatan dengan letak endsill x.
Pelaksanaan pengukuran model selanjutnya  dilakukan  untuk  setiap perubahan jarak apron, perubahan kedalaman  kolam  loncat  hidraulik, perubahan  tinggi endsill dan  perubahan angka Froude.  Pengukuran  dibatasi hanya pada  kondisi  loncat  hidraulik terbuka  (surface  flow),  tidak  termasuk kondisi loncat hidraulik terendam (submertged flow)

## 4. ANALISIS HASIL DAN PEMBAHASAN

## Kekasaran Saluran
Hasil  pengamatan  secara  fisik  untuk debit  sebesar Q =  4830  cm 3 /det diperoleh kedalaman  normal  sebesar yn =  5,7  cm.  untuk  bukaan  pintu y 0 =  1,25 cm posisi  awal  loncat  hidraulik  adalah Ls  =  191  cm dan  panjang  loncat hidraulik Lj = 31 cm. dengan menggunakan persamaan  8,  didapatkan nilai n = 0,0027.
Pengaruh  Kedalaman  dan  Panjang Kolam Olak (apron) terhadap

## Parameter Loncat Hidraulik
Gambar  4  memperlihatkan  bahwa  untuk suatu  nilai Froude number yang  sama, H/y 1 pada  kedalaman  kolam  1  (H  =  1 cm) < H/y 1 pada kedalaman kolam 2 (H = 2 cm) dengan nilai perbandingan 1 : 2. hal ini menunjukkan bahwa nilai y 1 tidak dipengaruhi oleh kedalaman kolam olak.
Gambar 4. Hubungan ekspermental antara Fr dan H/Y 1 pada kolam Koordinat – koordinat  yang  terbentuk pada gambar 4 menunjukkan bahwa nilai H/y 1 yang  senantiasa  berkurang  setiap dilakukan perubahan panjang La. Hal ini disebabkan perlawanan gaya yang terjadi di  ujung  hilir  kolam  olak  untuk mengimbangi  tekanan  hidrostatis  yang terdistribusi  disepanjang  apron  semakin kecil.  Jarak  untuk  mereduksi  distribusi tekanan hidrostatis yang semakin pendek (La  semakin  pendek),  menyebabkan tekanan  hidrostatis  yang  terjadi  akan semakin kecil (ditandai dengan perubahan  nilai  Fr 1 ). Pada  jarak  La tertentu,  kedalaman  kolam  olak  tidak efektif  lagi  dalam  memecah  energi.  Hal tersebut berkaitan dengan angka Froude, yaitu  pertama  apabila  nilai  Fr 1 terlalu kecil  hingga  melewati  batas vena contracta dan  menimbulkan  loncatan terendam (submerged jump).
Pengaruh  Kedalaman  dan  Panjang
Kolam Olak (apron) Terhadap

## Kedalaman Hilir (Tailwater Depths)
Dengan  menggambarkan  kedalaman (y) berserta  gaya  spesifik  (Fs)  pada  suatu penampang saluran dengan debit tertentu, diperoleh lengkung gaya spesifik. Untuk suatu nilai gaya spesifik, lengkung  gaya  spesifik  memiliki  dua kemungkinan kedalaman yaitu y 1 (initial depth)  dan  y 2 (sequent  depth).  Gaya spesifik  pada  penampang  1  (y 1 )  dan  2 (y 2 )  adalah  sama  (F 1 =  F 2 ),  jika  gaya – gaya luar dan berat air di bagian saluran yang  lurus  di  antara  kedua  penampang dapat diabaikan (Chow, V.T.,1959). Hubungan nilai y 1 dan y 2 hasil penelitian diperlihatkan pada  gambar  5a.  sebagai contoh,  untuk  suatu  model  G 2 K 1 , D 4 = 90 cm; y 1 = 0,0153 m dan y 2 = 0,0515 ; nilai Fs 1 = Fs 2 = 0,000553 m 3 .
Gambar  5b  memperlihatkan  hubungan antara  nilai  y 1 dan  y 3 (tailwater  depth) hasil  penelitian,  dimana  untuk  satu model,  dengan  panjang  apron  yang beragam,  nilai  y 3 tetap.  Pada  model G 1 K1, D 4 = 90 cm, y 1 = 0,0242 m; Fs 1 = 0,000412  m 3 ;  y 3 =  0,050  m.  nilai  y 3 terletak  di  sisi  atas  lengkung  gaya spesifik, yang artinya nilai y 3 >  y 2 . pada model G 2 K 1 , G 3 K 1 , G 2 K 2 , G 3 K 2 memperlihatkan  nilai  y 3 terletak  di  sisi dalam  lengkung  gaya  spesifik,  artinya bahwa nilai y 3 < y 2 .
Gambar  5c memperlihatkan  hubungan antara  nilai  y 1 dan  (y 3 +  H)  hasil penelitian.  Pada  model  G 1 K 1 ,  G 2 K 2 ,
G 3 K 2 ,  dan  sebagian  dari  model  G 2 K 1 ,memperlihatkan  bahwa  nilai  y 3 terletak di  sisi  atas  lengkung  gaya  spesifik (y 3 >y 2 ).  Untuk  model  G 1 K 1 adanya
penambahan  kolam  hanya  berfungsi untuk memperbesar pengaruh pembenaman,  karena  tanpa  adanya penambahan  kolam  nilai  y 3 >y 2 .  untuk model G 3 K 1 dan sebagian yang lain dari model  G 2 K 1 menunjukkan  bahwa  nilai y 3 terletak  di  sisi  dalam  lengkung  gaya spesifik  (y 3 <y 2 ),  untuk  itu  diperlukan penambahan  kedalaman  kolam  agar y 3 ≥ y 2, agar  peredaman  energi  lebih memadai.  Pada  model  G 2 K 1 untuk  nilai y 3 <y 2 terjadi pada apron dengan panjang 70  cm  s/d  120  cm.  untuk  nilai  y 3 =  y 2 terjadi  pada  apron  dengan  panjang  60 cm,  sedangkan  untuk  nilai  y 3 >y 2 terjadi pada  apron  dengan  panjang  20  s/d  50 cm.
Panjang Loncat Hidraulik pada Kolam Olak
Data – data  hasil  percobaan  dapat dipetakan  dengan  absis Froude number
Gambar 5. Nilai y 1 , y 2 & y 3 hasil eksperimen pada lengkung gaya spesifik
dan ordinatnya berupa rasio tak
berdimensi Lj/y 1 . Gambar 6
memperlihatkan  adanya  kesamaan  pola dalam  perubahan.  Semakin  besar  angka Fr,  maka  Lj/y 1 cenderung  bertambah besar.  Untuk  suatu  angka Froude yang sama,  panjang  loncat  hidraulik  yang terjadi semakin berkurang dengan adanya  penambahan  kedalaman  kolam, pengaturan tailwater depth dan penambahan tinggi endsill.
Panjang loncat hidraulik (Lj) pada kolam olak  dengan  H  =  1  cm  (model  G 2 K 1 , gambar  6a)  lebih  besar  dibandingkan kolam  olak  dengan  H=2  cm  (G 2 K 2 , gambar  6b).  untuk  kedalaman  kolam olak  tanpa endsill (model  G 2 K 1 ,  gambar 6a) lebih besar dibandingkan kolam olek dengan  penambahan endsill (G 2 K 1 , gambar  6d).  Begitu  pula  untuk  kolam olak  dengan  kedalaman  kolam  dengan penambahan taiwater depth akan
mengurangi panjang loncat hidraulik.
Gambar 6. Hubungan angka Froude terhadap perubahan panjang loncat hidraulik
Pengaruh  Angka Froude terhadap Kehilangan Energi (Head Losses)
Kehilangan  energi  pada  loncat  hidraulik adalah  sama  dengan  perbedaan  energi
spesifik  sebelum  dan  sesudah  terjadinya loncat  hidraulik.  Hasil  analisa  pada penelitian ini diperoleh hubungan  antara angka Froude yang  terbentuk  pada
bukaan  pintu  (y 0 ) terhadap  besarnya kehilangan  energi  yang  terjadi  dimulai dari  bukaan  pintu  (y 0 )  hingga tailwater depth (y 3 ).
Gambar 7. Perubahan angka Froude terhadap head losses ( ∆ E = y 0 – y 3 )
Gambar  7  memperlihatkan  terjadinya kehilangan  energi  total,  yaitu  dimulai dari hilir sluice  gate (y 0 )  sampai kedalaman  hilir  (y 3 ).  Untuk  suatu  angka Froude yang sama, besarnya head losses dipengaruhi oleh kedalaman kolam yaitu G 2 K 1 dan  G 2 K 2 serta  G 3 K 1 ,  dan  G 3 K 2 . semakin  besar  kedalaman  kolam, maka semakin besar pula energi yang direndam.
Perbandingan  y 2 /y 1 Hasil  Penelitian dengan y 2 /y 1 Hasil Formulasi Analitik Perbandingan  y 2 /y 1 hasil  penelitian dengan  y 2 /y 1 hasil  formulasi  analitik diplotkan dalam bentuk grafik hubungan antara  angka Froude dengan  y 2 /y 1 dan diperlihatkan pada  gambar  8.  posisi
loncat  hidraulik  dalam  kolam  olak mempunyai dua kemungkinan, yaitu : a. Apabila L/La<1 (gambar  3a),  maka persamaan  yang  berlaku  adalah persamaan 1 b. Apabila L/La ≥ 1 (gambar 3b), maka persamaan  yang  berlaku  adalah persamaan 11
Dari  Gambar  8a terlihat  bahwa  kolam olak  dengan  kedalaman  H  =  1  cm sebagian  besar  posisi  loncat  hidraulik terletak pada  L/La ≥ 1,  yang terjadi pada bukaan  pintu  y 0 =  1,25  cm.  hal  ini menunjukkan bahwa angka Froude yang besar menyebabkan  kolam  olak  belum dapat mengendalikan posisi loncat hidraulik.
Adanya  penambahan  kedalaman  kolam seperti  yang  terlihat  pada  gambar  8b menjelaskan  hampir  seluruh  loncat hydraulik  berada  pada posisi  di  dalam kolam olak.
Gambar 8d.  memperlihatkan  bahwa kolam  olak  dengan  penambahan endsill dapat menggeser loncat hidraulik ke arah hulu.  Sehingga  seluruh  posisi  loncat hidraulik berada dalam kolam olak.
Gambar 8. Perbandingan nilai y 2 /y 1 hasil penelitian dengan hasil formulasi analitik
Sedangkan pada gambar 8e. memperlihatkan bahwa adanya penambahan tailwater  depth dan endsill
pada kolam olak hanya akan memperbesar pengaruh pembenaman. Terdapat kecenderungan  yang  sama antara  hasil  penelitian  dengan  hasil  dari formulasi  analitis.  Untuk  kolam  dengan penambahan endsill dan pengaturan hilir antara  hasil  penelitian  dengan  hasil  dari formulasi  analitis  memberikan  hasil yang hampir sama.
Pengaruh  Angka Froude terhadap Posisi Awal Loncat Hidraulik
Froude number dipengaruhi oleh besarnya  debit  dan  tinggi  bukaan sluice gate. Hasil  analisa  pada  penelitian  ini diperoleh hubungan antara angka Froude yang  terbentuk  pada  awal  loncat hidraulik (y 1 ) terhadap posisi awal loncat hidraulik (Ls).
Gambar 9. Hubungan Angka Froude terhadap y 1 /Ls.
Koordinat – koordinat pada model G 3 K 1 , G 2 K 1 ,  dan  G 1 K 1 dengan  panjang  kolam olak yang sama (La sama) dihubungkan, maka  grafik  asitosis  ke  arah  titik  1, begitu diambil suatu angka Froude Fr  = 2,75 untuk model dengan panjang kolam olak yang sama :
a. Pada  model  G 2 K 1 y 1 /Ls  =  0,066 sehingga Ls = y 1 /0,066 b. Pada  model  G 2 K 2 y 1 /Ls  =  0,281 sehingga Ls = y 1 /0,281
Hal  yang sama terjadi pada model G 3 K 1 dan  G 3 K 2 .  Dari  angka – angka  tersebut
dapat disimpulkan bahwa adanya penambahan  kedalaman  kolam  olak mengakibatkan Ls semakin pendek.

## 5. KESIMPULAN
Berdasarkan pembahasan hasil penelitian  yang  telah  dilakukan  maka dapat disimpulkan bahwa :
a. Bila  diambil  suatu  nilai Froude
number untuk model dengan panjang kolam  olak  yang  sama,  adanya penambahan  kedalaman  kolam  olak mengakibatkan Ls semakin pendek
b. Untuk suatu  nilai Froude number yang  sama,  panjang  loncat  hidraulik yang  terjadi  semakin  berkurang dengan adanya penambahan kedalaman  kolam  olak,  penambahan tailwater  depth dan  penambahan tinggi endsill c. Penambahan kedalaman dan panjang kolam  olak,  penambahan tailwater depth dan penambahan tinggi endsill dapat memperbesar pengaruh peredaman energi
d. Penambahan endsill pada  ujung  hilir kolam olak dapat mengurangi terjadinya osilasi yang menyebabkan terjadinya  kerusakan  pada  tanggul – tanggul dari tanah dan riprap.

## 6. DAFTAR PUSTAKA
Bos, M.G., 1978, Discharge Measurement Structure, ICRI,
Netherlands Chow,  V.T.,  1959, Open  Channel
Hydraulics,
McGraw-Hill
Kogakusha Ldt, Tokyo
Hager,  W.H  dan  lin,  D.,  1992, Sill
Controlled  Energy  Dissipator,
Journal  ofHydraulics  Research vol. 30-92, halaman 165-181
Mossa  M.,  Petrillo  A.  dan  Chanson  H.,
2002, Tailwater Levelo Effects on Flow  Conditions  at  an  Abrupt Drop, Journal of Hydraulics
Research vol.40-4
Rangga  Raju,  K.G.,  1986, Aliran Melalui Saluran Terbuka, terjemahan Yan Piter Pangaribuan, Erlangga, Jakarta Sharma, S.K., 1988, Design of Errigation  Structure, S.  Chad  & Company (Pvt) LDT, New Delhi Wu, S. and Rajaratnam, N., 1995, Effect of  Buffles  on  Submerged  Flow, Journal of Hydraulics Engineering vol. 121, no. 9, September, ASCE Biodata :
Yuliman  Ziliwu,  Alumni  (S1)  Teknik Sipil  UTP Surakarta  (1994).  Pasca Sarjana  (S2)  Program  Magister  Teknik Sipil UNDIP Semarang  (2002),  Lahir di Ambukha  Nias  4  Juni  1966.  Dosen  di Program  Studi  Teknik  Sipil  UTP  dari tahun  1996  sampai  sekarang.  Pernah menjabat  sebagai  Pembantu  Dekan  II dari tahun 2002 sampai 2006 dan Wakil Dekan dari tahun 2006 sampai 2010.