## Analisis Geoteknik Terowongan Batuan Geurutee Aceh Menggunakan Metode Elemen Hingga

## RYAN ACHMAD FADHILLAH, INDRA NOER HAMDHAN
Jurusan Teknik Sipil, Institut Teknologi Nasional, Bandung e-mail: daderaf@gmail.com

## ABSTRAK
Gunung  Geurutee  Aceh  memiliki  kondisi  standar  alinyemen  vertikal  maupun horizontal  masih  belum  terpenuhi,  begitu  juga  waktu  tempuh  efektif.  Timbul wacana pembangunan terowongan di kawasan Gunung Geurutee Aceh. Analisis geoteknik terowongan batuan Geurutee dibutuhkan untuk mengetahui gaya dalam yang  terjadi  pada  terowongn  tersebut  dengan  adanya  sistem  perkuatan sementara, nilai faktor keaamanan (FK) , dan face stability menggunakan metode elemen  hingga.  Metode  konstruksi  terowongan  menggunakan  metode  New Austrian  Tunneling  Methode  (NATM)  dan  Metode  penggalian  mengacu  pada metode  Japanese  Society  of  Civil  Engineers  (JSCE).  Model  yang  dibuat  adalah terowongan  dengan  perkuatan  sementara  dengan  tanpa  perkuatan  sementara, untuk  membandingkan  total  displacement  yang  terjadi  pada  model  tersebut dimana perkuatan sementara yaitu menggunakan rock bolt dan shotcrete. Kedua model  tersebut  menghasilkan  nilai  SF  7,076  untuk  model  dengan  perkuatan sementara dan 4,525 untuk model tanpa perkuatan.
Kata kunci : terowongan batuan,  analisi geoteknik, face stability, metode elemen hingga NATM, JSCE, rock bolt, shotcrete.

## ABSTRACT
Mount Geurutee Aceh has a standard condition horizontal and vertical alignment that  are  still  not  fulfilled,  as  well  as  effective  travel  time.  Arising  discourse construction of tunnels in the area of Mount Geurutee Aceh. Geurutee rock tunnel geotechnical analysis is required to determine the force that occurs in the tunnel with temporary reinforcement system, the value of the safety factor (SF), and face stability using the finite element method. Tunnel construction method using the New Austrian Tunneling Method (NATM) and the excavation method refers to the method  of  Japanese  Society  of  Civil  Engineers  (JSCE).  The  model  created  is temporary  reinforcement  system  tunnels  and  without  temporary  reinforcement system,  to  compare  the  total  displacement  that  occurred  on  that  model  where temporary reinforcement system while using rock bolt and shotcrete. Both models come  SF  value  7.076  for  the  model  with  temporary  reinforcement  system  and 4,525 for the model without temporary reinforcement system.
Keywords :  rock  tunnel,  geotechnical  analysis,  face  stability,  finite  element method NATM, JSCE, rock bolt, shotcrete.

## 1. PENDAHULUAN
Terowongan adalah suatu tempat perpindahan alternatif yang berada di bawah permukaan tanah  atau  gunung . Pembangunan  terowongan  Geurutee  ini  diharapkan  dapat  mengurangi angka  kecelakaan,  mempercepat  arus  lalu  lintas  dan  melancarkan  proses  pengangkutan barang  ke  wilayah  barat  selatan  Aceh . Metode  konstruksi  menggunakan  metode  NATM disebabkan  untuk  mencegah  terjadinya  keruntuhan  di  muka  terowongan  ketika  proses konstruksi. Penggalian dengan metode JSCE memilih cara membagi proses penggalian pada terowongan tersebut untuk mencegah keruntuhan pada saat proses penggalian.
Pada  tugas  akhir  ini  akan  dilakukan  analisis  geoteknik  terowongan  batuan.  Agar  metode konstruksi  dan  penggalian  yang  akan  diterapkan  sesuai  dengan  kondisi  aslinya,  maka dibutuhkan perhitungan yang tepat dengan melakukan pemodelan terlebih dahulu. Pada studi kasus tugas akhir ini digunakan program PLAXIS 3D AE yang berbasis metode elemen hingga .
Manfaat penelitian ini, yaitu dapat menambah pengetahuan untuk menganalisis dan mengolah data tanah dan batuan di suatu tempat juga dapat mendesain terowongan ( tunnel) untuk jalan raya pada kondisi batuan tertentu. Gambar 1 merupakan gambar rencana trase terowongan Geurutee.
Gambar 1. Perencanaan trase Terowongan Geurutee
(Sumber: Balai Geoteknik, 2015)

## Analisis Geoteknik Terowongan Batuan Geurutee Aceh

## 2. TINJAUAN PUSTAKA

## 2.1 Terowongan
Terowongan sebagai prasarana transportasi umumnya tertutup di seluruh sisi kecuali di kedua ujungnya  yang  terbuka  pada  lingkungan  luar.  Beberapa  ahli  teknik  sipil  mendefinisikan terowongan sebagai sebuah tembusan di bawah permukaan yang memiliki panjang minimal 0,1 mil atau 0,161 km, jika lebih pendek dari 0,1 mil disebut underpass.

## 2.2 Metode Penggalian Terowongan
Kondisi  tanah  atau  batuan  dilapangan  dan  kondisi  trase  eksisting  di  lokasi  konstruksi terowongan,  mengakibatkan  dibutuhkannya  metode  penggalian  terowongan  sesuai  dengan kondisi-kondisi  tersebut.  Berikut  di  bawah  ini  merupakan  jenis-jenis  metode  penggalian terowongan:
1. Ledakan dan Bor ( Drill and Blast).
2. Tunnel Boring Machine (TBM).
3. Roadheader.
4. New Austrian Tunneling Methode (NATM) .

## 2.3 Rock Mass Rating (RMR)
Rock Mass Rating (RMR) atau juga dikenal dengan Geomechanics Classification dikembangkan oleh  Bieniawski  pada  tahun  1972-1973.  Metode  klasifikasi  RMR  merupakan  metode  yang sederhana dalam penggunaannya, dan parameter-parameter yang digunan dalam metode ini dapat  diperoleh  baik  dari  data  lubang  bor  maupun  dari  pemetaan  struktur  bawah  tanah. Tujuan  dari  system  RMR  adalah  untuk  mengklasifikasi  kualitas  massa  batuan  dengan menggunakan  data  permukaan,  untuk  memandu  metode  penggalian  dan  juga  untuk memberikan  rekomendasi  penggalian  menggunakan  perkuatan  semantara  dengan  tanpa perkuatan sementara.

## 2.3.1 Klasifikasi Massa Batuan Metode RMR
Dalam  mengklasifikasikan  massa  batuan  berdasarkan  sistem  klasifikasi  RMR,  Bieniawski menggunakan enam parameter utama yang dijumlahkan untuk memperoleh nilai total RMR, yaitu :
1. Uniaxial Compressive Strength (UCS).
2. Rock Quality Designation (RQD).
3. Jarak antar (spasi) kekar ( Spacing of discontinuities).
4. Kondisi kekar ( Condition of discontinuities).
5. Kondisi air tanah ( Groundwater conditions).
6. Orientasi Kekar ( Orientation of discontinuites).

## 2.3.2 Penggunaan Rock Mass Rating (RMR)
Kondisi massa batuan dievaluasi untuk setiap bidang diskontinu yang ada (Bieniawski, 1989). RMR  dapat  digunakan  sebagai  panduan  memilih  penyangga  terowongan.  Panduan  ini tergantung pada beberapa faktor spserti kedalaman lubang bukaan dari permukaan, ukuran dan  bentuk  terowongan  serta  metode  penggalian  yang  akan  di  pakai  (Bieniawski,  1989) dibawah ini Tabel 1 dari RMR (setelah Bieniawski, 1989).

## Ryan Achmad Fadhillah, Indra Noer Hamdhan
Tabel 1. Rock Mass Rating System (after Bieniawski, 1989)
Point load strength index > 10 MPa 4 - 10 MPa 2 - 4 MPa 1 - 2 MPa Uniaxial comp. strength > 250 MPa 100 - 250 MPa 50 - 100 MPa 25 - 50 MPa 5 - 25 MPa 1 - 5 MPa < 1 MPa 15 12 7 4 2 1 0 90% - 100% 75% - 90% 50% - 75% 25% - 50% 20 17 13 8 > 2 m 0,6 - 2 m 200 - 600 mm 60 - 200 mm 20 15 10 8 Very rough surface, not continuous, no separation, unweathered wall rock Slightly rough, surfaces, separation < 1 mm, slightly weathered walls Slightly rough, surfaces, separation < 1 mm, highly weathered walls Slickensided surfaces or Gouge < 5 mm thick or Separation 1-5 mm Continuous 30 25 20 10 Inflow per 10 m tunnel length (l/m) None < 10 10 -- 25 25 - 125 joint water press / major principal 0 < 1,0 0,1 - 0,2 0,2 - 0,5 General conditions Completely Dry Damp Wet Dripping 15 10 7 4 Very favourable Favourable Fair Unfavourable Tunnels & mines 0 -2 -5 -10 Foundations 0 -2 -7 -15 Slopes 0 -5 -25 -50 100 - 81 80 - 61 60 - 41 40 - 21 I II III IV Very good rock Good rock Fair rock Poor rock I II III IV 20 yrs for 15 m span 1 year for 10 span 1 week for 5 m span 10 hrs for 2,5 m span > 400 300 - 400 200 - 300 100 - 200 > 45 35 - 45 25 - 35 15 - 25 < 1 m 1 - 3 m 3 - 10 m 10 - 20 m 6 4 2 1 None < 0,1 mm 0,1 - 1 mm 1 - 5 mm 6 5 4 1 Very rough Rough Slightly rough Smooth 6 5 3 1 None Hard filling < 5 mm Hard filling > 5 mm Soft filling < 5 mm 6 4 2 2 Unweathered Slightly weathered Moderately weathered Highly weathered 6 5 3 1
Drive againts dip - Dip 20 - 45 (Deg) Decomposed 0 F. EFFECT  OF DISCONTINUITY STRIKE AND DIP ORIENTATION IN TUNNELLING Fair
Unfavourable
Dip 45 - 90 (Deg) Dip 20 - 45 (Deg) Very unfavourable Fair
Dip 0 - 20 - Ireespective of strike (Deg) Fair Drive with dip - Dip 45 - 90 (Deg)
Drive with dip - Dip 20 - 45 (Deg)
Very favourable
Favourable
Drive againts dip - Dip 45 - 90 (Deg)
Strike perpendicular to tunnel axis
Strike parallel to tunnel axis
> 20 m 0
> 5 mm
0
Slickensided
0 Roughness Rating Infilling Rating Weathering Rating Soft filling > 5 0

## E. GUIDELINES FOR CLASSIFICATION OF DISCONTINUITY CONDITIONS
Discontinuity length Rating Separation
Rating
D. MEANING OF ROCK CLASSES
Class number Average stand-up time
Cohesion of rock mass (kPa)
Friction angle of rock mass (deg)
V 30 min for 1 m span < 100 < 15 C. ROCK MASS CLASSES DETERMINED FROM TOTAL RATING Rating Class number
Description
< 21 V Very poor rock Rating 0 B. RATING ADJUSMENT FOR DISCONTINUITY ORIENTATIONS (See F) Strike and dip orientations Ratings Very Unvfavourable -12 -25
5 Ground Water > 125 > 0,5 Flowing 4 Condition of discontinuities (See E) Rating Soft gouge > 5 mm thick or Separation > 5 mm Continuous 0 2 Drill core Quality RQD Rating < 25%
3
3 Spacing of discontinuities Rating < 60 mm 5 A. CLASIFICATION PARAMETERS AND THEIR RATINGS Parameter
Range of values
1 Strength of intact rock material Rating For this low range uniaxial compresive test is prefered

## Analisis Geoteknik Terowongan Batuan Geurutee Aceh
Gambar 2. Kolerasi antara 𝐄𝐦 dengan RMR
(Sumber: 2007.-Evaluation-of-the-formation-modulus-of-rock-masses-using-RMR)
Pada Gambar 2 merupakan hubungan antara modulus elastisitas batuan dengan nilai RMR menurut Bienawski dan Sefarin & Preira.
1. Formula Bieniawski (1978) untuk batuan keras (𝑅𝑀𝑅  >  50) , dapat dilihat pada Persamaan 1
𝐸 𝑚 =  2 𝑅𝑀𝑅 –  100 𝐺𝑃𝑎 … (1)
2. Formula Serafin dan Pereira (1983) untuk batuan keras (𝑅𝑀𝑅  <  50) , dapat dilihat pada Persamaan 2
𝐸 𝑚 = 𝜋𝑟 2 =   10 (𝑅𝑀𝑅−10)
40 … (2) 3. Modulus Elastisitas (Ei) , didapat pada formula, dapat dilihat pada Persamaan 3
𝐸 𝑚 =  𝐸𝑖 . 𝑅𝑀𝐹
… (3) 4. Dimana nilai Modulus Reduction Factor (MRF) menurut Nicholson dan Bieniawski (1990), dapat dilihat pada Persamaan 4
𝑀𝑅𝐹  =   𝐸 𝑚 / 𝐸𝑖  =  0,0028 𝑅𝑀𝑅 2 = 𝜋𝑟 2 +  0,9 . 𝑒 (𝑅𝑀𝑅/22,82)
… (4)

## 1.4 Japanese Society of Civil Engineer (JSCE)
Pedoman  metode  perencanaan  penggalian  dan  sistem  perkuatan  terowongan  jalan  pada media  campuran  tanah-batuan  di  Indonesia  umumnya  menggunakan  JSCE  dikarenakan metode ini dapat menjadi acuan untuk perencanaan penggalian terowongan jalan di media campuran tanah-batuan .

## 2.4.1 Metode Penggalian
Tipikal  ilustrasi  penampang  melintang  penggalian  berdasarkan  pembagian  muka  bidang galian,  kondisi  media  yang  sesuai  dan  kategori  batuan/tanah  dapat  dilihat  pada Tabel  2 dengan  tahapan  penggalian  mengacu  pada  penomoran  menggunakan  nomor  terkecil  yang menjadi urutan penggalian pertama yang dilakukan.
Tabel 2. Pemilihan Metode Penggalian dan Karakteristiknya
(Sumber : Japanese Society of Civil Engineer (JSCE) )

## 3. ANALISIS DATA

## 3.1 Data Terowongan
Pada pemodelan terowongan, dibutuhkan data berupa jenis tanah, jenis batuan, parameter tanah, parameter batuan, topografi pegunungan dan data penampang terowongan. Berikut data-data teknis yang di dapat untuk pemodelan terowongan:
1. Jenis batuan : Breksi Granitik 2. Parameter batuan : 𝛾 𝑢𝑛𝑠𝑎𝑡 , 𝛾 𝑠𝑎𝑡 , 𝑐,  , 𝐸’, 𝑣’
3. Parameter shotcrete
: 𝑑, 𝛾, 𝐸, 𝑣
4. Parameter rock bolt : 𝛾 , 𝐸, Ø , pile type

## 3.2 Diagram Alir Metodologi Analisis
Metodologi analisis yang dijadikan dasar studi ini dapat dijabarkan dalam bentuk diagram alir ( flow  chart)  pada Gambar  3 .  Diawali  dari  perumusan  masalah,  studi  pustaka  yaitu  tahap pengumpulan bahan – bahan seperti referensi persiapan, data terowongan seperti yang sudah dijelaskan  pada  sub  bab  sebelumnya,  Parameter  dan  pengolahan  data  diamana  parameter didapat  dari  uji  lab  maupun  lapangan,  pemodelan  terowongan  pada  PLAXIS  3D  yaitu memodelkan geometri bentuk penampang terowongan sesuai dengan gambar rencana, dari hasil  perhitungan  PLAXIS  3D  didapat  dua  jenis  hasil  yaitu  hasil  dari  stabilitas  berupa face stability, faktor keamanan dan hasil dari gaya-gaya yang terjadi berupa deformasi, gaya aksial, gaya geser, dan momen. Terakhir dari diagram alir berupa pembahasan dan kesimpulan dari hasil analisi yang telah dilakukan.

## Analisis Geoteknik Terowongan Batuan Geurutee Aceh
Start Perumusan Masalah Data Terowongan: 1. Jenis Batuan 2. Parameter Batuan, shotcrete, dan rock bolt 3. Topografi Pegunungan 4. Data Penampang Terowongan Mengolah Data dan Menentukan Paremeter yang Digunakan Pemodelan Terowongan Menggunakan PLAXIS 3D dengan Metode Elemen Hingga Analisis Stabilitas ( Face) menggunakan PLAXIS 3D AE Output : 1. Diagram Gaya Aksial 2. Diagram Gaya Momen 3. Diagram Gaya Geser 4. Diagram Deformasi 5. Faktor Keamanan Tahapan Penggalian Analisa Gaya-gaya yang Terjadi Pada Elemen Perkuatan Terowongan Pembahasan dan Kesimpulan Selesai Output : 1. Bidang Longsor ( Face) 2. Faktor Keamanan Tahapan Penggalian Studi Pustaka

## Gambar 3. Diagram alir metodologi analisis
Data parameter batuan, diperoleh dari sampel yang telah di ambil dari tempat yang akan di tinjau. Tabel 3 menunjukan data parameter batuan di gunung Geurutee Aceh.

## Ryan Achmad Fadhillah, Indra Noer Hamdhan
Tabel 3. Parameter Batuan di Titik BT-2 Gunung Geurutee Aceh
Lapisan  Jenis Batuan RQD (%) 𝜸 𝑼𝒏𝒔𝒂𝒕 (kN/m 3 ) 𝜸 𝒔𝒂𝒕 (kN/m 3 ) v’ E (kN/m 2 ) C (kN/m 2 )  (°) 1 Breksi Granitik 24 -73 26,6 26,7 0,15 6,3E5 250 30 2 Breksi Granitik  74 - 98 24,9 25,2 0,15 1,9E6 350 40
(Sumber: Balai Geoteknik, 2015)

## 3.3  Pemodelan Geometri Terowongan Pada PLAXIS 3D AE
Desain bentuk penampang terowongan dimodelkan separuh dari bentuk aslinya dikarenakan simetris dan model materialnya sudah disederhanakan sehingga tidak didesain utuh. Gambar 4 merupakan desain struktur terowongan dengan panjang 19,5 m dan perkuatan sementara pada PLAXIS 3D AE.
Gambar 4. Desain struktur Terowongan Geurutee pada PLAXIS 3D AE

## 3.4  Hasil Perhitungan PLAXIS 3D AE
Hasil  perhitungan  yang  diambil  dari  pemodelan  terowongan  Geurutee  yaitu  nilai Total Displacement, Shear Force, Axial Force, Moments, dan Safety Factor. Nilai Displacement akan dibandingkan  antara  pemodelan  terowongan  dengan  perkuatan  sementara  dan  pemodelan terowongan  tanpa  perkuatan  sementara. Gambar  5 merupakan  hasil  perhitungan  output PLAXIS 3D AE. Dibawah ini merupakan grafik hasil total displacement tiap panjang penggalian 3 m penggalian.

## Analisis Geoteknik Terowongan Batuan Geurutee Aceh

## Gambar 5. Perbandingan deformasi pemodelan terowongan dengan perkuatan sementara dan tanpa perkuatan
Dari Gambar 5 tersebut dapat di lihat terowongan dengan perkuatan sementara lebih kecil terjadi  deformasi  dibandingkan  tanpa  perkuatan  sementara,  dimana  dengan  perkuatan terdeformasi sebesar 7,64E-03 m, sedangkan tanpa perkuatan sebesar 8,60E-3 m. Perbedaan deformasi yang kecil dikarenakan kondisi batuan pada Pegunungan Geurutee terbilang mantap shingga dengan perkuatan sementara pun perbedaan deformasi tidak terlalu signifikan. Hasil perhitungan PLAXIS 3D AE tiap panjang galian dapat dilihat pada Gambar 6 .
Gambar 6. Deformasi hasil perhitungan PLAXIS 3D AE
0.00E+00 1.00E-03 2.00E-03 3.00E-03 4.00E-03 5.00E-03
6.00E-03
7.00E-03
8.00E-03 9.00E-03
1.00E-02 0 3 6 9 12 15 18 21 D efo rmas i (m) Panjang Penggalian (m) Dengan Perkuatan Tanpa Perkuatan
3 m 6 m 9 m 12 m 15 m 18 m 19,5 m 3 m 6 m 9 m 12 m 15 m 18 m 19,5 m Model Terowongan Dengan Perkuatan Model Terowongan Tanpa Perkuatan
Hasil perhitungan selanjutnya berupa nilai gaya geser, gaya aksial, dan momen. Gaya - gaya dalam tersebut dapat dilihat pada struktur perkuatan sementara dan nilai gaya gaya dalam tersebut dapat dilihat pada Tabel 4 .
Tabel 4. Hasil Perhitungan Gaya - Gaya Dalam PLAXIS 3D AE
Tahap Galian (m) Shear Force Axial Force Moments Min (kN/m) Max (kN/m) Min (kN/m) Max (kN/m) Min (kNm/m) Max (kNm/m) 3 -78,65 93,54 -197,2 138,6 -2,788 16,92 6 -136,4 98,92 -361,1 182,4 -4,335 21,46 9 -164,3 94,27 -445,2 183,4 -6,804 22,73 12 -183,7 97,34 -496,2 193,2 -6,559 24,57 15 -196,2 92,96 -526,1 155,1 -6,492 24,85 18 -205,3 98,19 -545 150,9 -6,442 25 19.5 -208,2 65,94 -551,7 142 -6,426 25,32
Dari Tabel 4 dapat dilihat nilai Shear Force terbesar yaitu -208,2 kN/m terjadi pada sistem perkuatan  sementara  ketika  sistem  perkuatan  sementara  terpasang  pada  panjang  19,5  m karena jika penampang lebih luas maka gaya geser yang terjadi pun lebih besar. Sama halnya Axial Force terbesar terjadi pada perkuatan sementara yang terpasang pada panjang 19,5 m sebesar  -551,7  kN/m.  Momen  dominan  terjadi  pada  atap  terowongan  karena  bagian  atas, momen terbesar 25,32 kN m/m ketika perkuatan sementara terpasang pada pajang 19,5 m.
Gambar 7. Nilai SF terhadap tahapan penggalian tanpa perkuatan
Pada Gamba 7 menyimpulkan bahwa pada panjang Penggalian 3 m nilai faktor keamanan (FK) 4,525 lebih besar dari panjang penggalian 19,5 m yaitu 2,644 karena, pada tahap awal penggalian batuan yang tergali tidak sebanyak setelah tergali di 19,5 m sehingga kemungkinan terjadi keruntuhan sangat kecil, sehingga nilai FK di panjang penggalian 3 m lebih besar dari panjang penggalian 19,5 meski pada penggalian tersebut masih aman karena syarat nilai FK untuk batuan dikatakan aman jika > 1,5 . Terakhir adalah pola keruntuhan face stability yang terjadi  pada  pemodelan  terowongan  Geurutee  ketika  tahapan  penggalian,  keruntuhan dominan  terjadi  pada  bagian  atap  terowongan. Gambar  8 merupakan  tampilan  dari  pola keruntuhan yang terjadi pada panjang galian 19,5 m.
0 1 2 3 4 5 6 7 8 0 3 6 9 12 15 18 21 Fak to r Kea m an an Panjang Penggalian (m) Tanpa Perkuatan Dengan Perkuatan

## Analisis Geoteknik Terowongan Batuan Geurutee Aceh
Gambar 8. Face stability pada pemodelan Terowongan Geurutee

## 4. KESIMPULAN
Besarnya deformasi yang terjadi pada terowongan tanpa perkuatan dan dengan  perkuatan didapat  nilai  deformasi  yang  tidak  jauh  berbeda  dikarenakan  kondisi  massa  batuan pegunungan Geurutee terbilang mantap.
Gaya  dalam  yang  terjadi  gaya  aksial,  gaya  geser,  dan  momen  terbesar  pada    perkuatan sementara setelah terpasang 19,5 m karena makin luas penampang maka beban yang diterima pun lebih besar.
Kedua model memiliki nilai SF yang di atas batas aman selama panjang terowongan tersebut 19,5  m  tetapi,  pemodelan  tanpa  perkuatan  sementara  diperkirakan  akan  mengalami penurunan nilai SF sampai < 2 pada penggalian selanjutnya, sehingga pemodelan terowongan tanpa  perkuatan  sementara  tidak  dianjurkan.  Pada  metode  konstruksi  sesungguhnya  di lapangan terowongan tanpa perkuatan tidak dianjurkan karena dari faktor akibat blasting yang memungkinkan terjadinya keruntuhan pada tahap penggalian terowongan tersebut.

## DAFTAR RUJUKAN
Broere, W. (2001). Tunnel Face Stability & New CPT  Application. ISBN 90-407-2215-3 DUP Science. Netherlands.
Han  Kim,  Seung.  (2010).  Large  Tunnels  for  Transportation  Purposes  and  Face  Stability  of Mechanically Driven Tunnels in Soft Ground. Austin: The University of Texas.
Hoek, Evert. (2006). Pratical Rock Engineering. Evert Hoek Consulting Engineer Inc, Columbia. Japan Society of Civil Engineers (JSCE). (2007). Standard Specifications for Tunneling-2006:
Mountain Tunnels. Tokyo: JSCE.
Menteri Pekerjaan Umum dan Perumahan Rakyat.  (2015). Metode  Perencanaan Penggalian dan  Sistem  Perkuatan  Terowongan  Jalan  Pada  Media  Campuran  Tanah-Batuan. Jakarta: Kementerian Pekerjaan Umum dan Perumahan Rakyat.
U.S. Departement of Transportation Federal Highway Administration. (2009). Technical Manual for  Design  and  Construction  of  Road  Tunnels  –  Civil  Elements.  Washington  D.C: National Highway Institute.