ESTIMASI NUMERIK GAYA IMPAK BENTURAN TONGKANG PADA PILAR JEMBATAN MAHAKAM MENGGUNAKAN METODE INTEGRASI SIMPSON

BAB I: PENDAHULUAN

Latar Belakang

Keselamatan infrastruktur transportasi perairan pedalaman (inland waterways) merupakan salah satu pilar krusial dalam menopang aktivitas ekonomi daerah, khususnya pada wilayah yang mengandalkan jalur sungai sebagai urat nadi distribusi logistik utama. Jembatan Mahakam yang terletak di Samarinda, Kalimantan Timur, merupakan infrastruktur vital yang menghubungkan antarwilayah sekaligus menjadi perlintasan strategis bagi armada industri maritim, terutama tongkang pengangkut batubara. Namun, posisi pilar jembatan yang berada di tengah alur pelayaran aktif menjadikannya sangat rentan terhadap risiko insiden kejut berupa tumbukan hidrodinamis dari sarana apung yang melintas. Keresahan personal mendalam sebagai putra daerah sekaligus mahasiswa teknik perkapalan muncul setiap kali mendengar kabar berulangnya peristiwa pilar Jembatan Mahakam dihantam oleh tongkang batubara berukuran masif. Kesadaran moral dan tanggung jawab keilmuan inilah yang mendasari pentingnya sebuah kajian ilmiah terukur untuk membedah karakteristik mekanis dari fenomena tumbukan tersebut. Melalui pendekatan rekayasa perkapalan, tugas akademis metode numerik ini diamplifikasi menjadi sebuah misi nyata untuk memberikan kontribusi kuantitatif terhadap mitigasi risiko bencana maritim lokal, dengan tujuan akhir merumuskan estimasi gaya rusak tumbukan secara presisi sebagai landasan evaluasi keselamatan struktur publik.

BAB II: METODOLOGI DAN MODEL MATEMATIKA

Definisi Parameter Fisik dan Analisis Hidrodinamika

Pendekatan kuantitatif dimulai dengan menetapkan skenario batas atas (worst-case scenario) berdasarkan dimensi armada tongkang standar yang beroperasi di Sungai Mahakam, yaitu tipe barge 300 kaki (feet). Dalam kondisi bermuatan penuh, bobot mati atau massa kering dari tongkang tersebut diestimasikan sebesar 7.000 ton atau setara dengan 7.000.000 kg. Di dalam ranah teknik perkapalan, perhitungan momentum tumbukan tidak boleh hanya melibatkan massa kering kendaraan apung. Ketika lambung kapal melaju membelah fluida, terdapat sejumlah volume air di sekitar lambung yang ikut menempel dan terakselerasi searah dengan pergerakan kapal. Fenomena fisik ini didefinisikan sebagai Massa Tambah (Added Mass). Mengingat bentuk geometris tongkang yang cenderung menyerupai kotak dengan nilai koefisien blok (block coefficient) yang tinggi, koefisien massa tambah diestimasikan sebesar 20% dari massa aktualnya. Jika efek hidrodinamika air ini diabaikan, maka total energi kinetik yang dihasilkan akan terhitung jauh lebih kecil dari kondisi riil di lapangan, sehingga berisiko menciptakan kegagalan fatal pada desain pelindung pilar. Melalui perumusan total massa dinamis (M = M_tongkang + M_added), didapatkan nilai massa total sebesar 7.000.000 kg + 1.400.000 kg = 8.400.000 kg. Armada ini diasumsikan bergerak hanyut akibat pengaruh kegagalan mesin maupun dorongan arus pasang surut sungai yang deras, menghasilkan kecepatan total ($v$) sebesar 2 m/s atau setara dengan 4 knot.

Pemodelan Fungsi Gaya Non-Linear

Ketika lambung baja tongkang melakukan kontak fisik pertama kali dengan struktur beton tegar pilar jembatan, seluruh momentum dari massa dinamis sebesar 8.400.000 kg tersebut harus diredam hingga kecepatannya turun menjadi nol dalam rentang waktu yang sangat singkat. Durasi waktu total persentuhan dari awal hantaman hingga kapal berhenti total diestimasikan berjalan selama delta_t = 1,2 detik. Berdasarkan Hukum Newton II dan Hukum Momentum-Impuls, efek beban dinamik ini dihitung sebagai Impuls, yang merupakan fungsi integrasi gaya terhadap waktu dengan persamaan dasar:

Tantangan utama dari analisis struktural ini adalah sifat gaya hantaman F(t) yang tidak konstan dari detik ke detik. Gaya bernilai nol pada awal kontak (t = 0), kemudian melonjak secara non-linear mencapai titik beban puncak tertinggi (F_max) di tengah durasi benturan (t = 0,6 detik) seiring dengan terjadinya deformasi plastis material lambung, dan akhirnya meluncur turun kembali ke angka nol ketika seluruh energi kinetik kapal telah habis terserap (t = 1,2 detik). Untuk merepresentasikan perilaku fisik melengkung tajam ini secara matematis, diimplementasikan pendekatan idealisasi fungsi Pulsa Setengah Sinus (Half-Sine Pulse Model) dengan rumus:

BAB III: IMPLEMENTASI METODE NUMERIK DAN KALKULASI

Diskretisasi Domain Waktu

Karena persamaan fungsi sinus tersebut bersifat kontinu, penggunaan matematika analitik murni akan membatasi pengembangan model jika dihadapkan pada variasi data diskrit di lapangan. Oleh karena itu, Metode Integrasi Numerik Simpson 1/3 diterapkan untuk memecahkan masalah integrasi ini. Langkah awal algoritma Simpson 1/3 adalah melakukan diskretisasi domain waktu dengan membagi durasi total 1,2 detik menjadi $n = 6$ segmen interval kecil yang sama besar (memenuhi syarat mutlak aturan Simpson di mana jumlah segmen wajib bernilai genap). Lebar langkah interval (h) dihitung melalui persamaan h = \frac{t_{\text{akhir}} – t_{\text{awal}}}{n}$, sehingga diperoleh $h = \frac{1,2 – 0}{6} = 0,2$ detik. Proses gridding ini mendirikan 7 buah titik koordinat waktu diskrit, dimulai dari t_0 = 0,0 s, t_1 = 0,2 s, t_2 = 0,4 s, t_3 = 0,6 s, t_4 = 0,8 s, t_5 = 1,0 s, hingga t_6 = 1,2 s.

Operasi Tabel Iterasi Spreadsheet

Dengan menetapkan estimasi nilai batas atas gaya puncak sebesar F_max = 27.500.000 Newton berdasarkan pendekatan konversi energi kinetik total, nilai gaya sesaat F(t) pada masing-masing titik grid dihitung menggunakan bantuan spreadsheet. Algoritma Simpson 1/3 mengeksekusi data ini dengan memberikan bobot pengali kuadratik parabola yang spesifik, yaitu mengikuti pola koefisien (1, 4, 2, 4, 2, 4, 1):

• Pada titik t_0 = 0,0 detik, nilai gaya sesaat adalah F(0,0) = 27.500.000 x sin(0 rad) = 0 Newton, dengan hasil perkalian bobot bernilai 0.
• Pada titik t_1 = 0,2 detik, nilai gaya sesaat melonjak menjadi F(0,2) = 27.500.000 x sin(0,5236 rad) = 13.750.000 Newton, dikalikan bobot kedua menghasilkan 55.000.000.
• Pada titik t_2 = 0,4 detik, gaya meningkat ke angka F(0,4) = 27.500.000 x sin(1,0472 rad) = 23.815.723 Newton, dikalikan bobot ketiga menghasilkan 47.631.446.
• Pada titik puncak t_3 = 0,6 detik, kurva mencapai kekuatan penuh di mana F(0,6) = 27.500.000 x sin(1,5708 rad) = 27.500.000 Newton, dikalikan bobot keempat menghasilkan 110.000.000.
• Pada titik t_4 = 0,8 detik, gaya masuk ke fase pelepasan energi menjadi F(0,8) = 27.500.000 x sin(2,0944 rad) = 23.815.723 Newton, dikalikan bobot kelima menghasilkan 47.631.446.
• Pada titik t_5 = 1,0 detik, gaya menurun ke angka F(1,0) = 27.500.000 x sin(2,6180 rad) = 13.750.000 Newton, dikalikan bobot keenam menghasilkan 55.000.000.
• Pada titik akhir t_6 = 1,2 detik, saat kapal berhenti sempurna, gaya kembali menyentuh F(1,2) = 27.500.000 xsin(3,1416 rad) = 0 Newton, dengan perkalian bobot terakhir bernilai 0.

Seluruh baris kalkulasi tersebut diakumulasikan ke bawah melalui fungsi penjumlahan spreadsheet, sehingga diperoleh nilai mutlak TOTAL sebesar:

TOTAL = 0 + 55.000.000 + 47.631.446 + 110.000.000 + 47.631.446 + 55.000.000 + 0 = 315.262.892

Langkah finalisasi dari operasi integrasi numerik ini adalah mengalikan angka TOTAL tersebut dengan faktor pengali sepertiga langkah ($h/3$), guna mematangkan parameter luasan total Impuls tumbukan:

BAB IV: ANALISIS LANJUTAN DAN SIMULASI KOMPUTASI

Integrasi Ekosistem CAD-CAE-CAM dan Teori Fluida Komputasi

Wawasan rekayasa dalam penelitian ini diperluas dengan meninjau bagaimana industri maritim modern menyelesaikan masalah interaksi hidro-struktural yang kompleks melalui integrasi platform Computer Aided Design (CAD) untuk pemodelan lambung 3D, Computer Aided Engineering (CAE) untuk simulasi dinamika fluida, serta Computer Aided Manufacture (CAM) untuk fabrikasi pelindung pilar. Studi kasus benturan ini berada penuh di dalam domain CAE, di mana validasi pergerakan fluida di sekitar pilar jembatan dipelajari melalui perangkat lunak Dinamika Fluida Komputasi (Computational Fluid Dynamics / CFDSOF) yang berbasis Metode Volume Hingga (Finite Volume Method).

Analisis fenomena hidrodinamika air yang terjepit di antara lambung tongkang dan pilar beton dapat diturunkan ke dalam ruang dua dimensi (2D) menggunakan basis matematis Persamaan Diferensial Parsial (Partial Differential Equations / PDE) Navier-Stokes. Persamaan momentum kontinu ini mengatur pergerakan fluida akibat variasi tekanan dan viskositas lokal yang diturunkan untuk komponen arah kecepatan sumbu-x (u) dan sumbu-y (v) secara simultan:

1. Persamaan Momentum Horizontal (Sumbu-x):u * (dU/dx) + v * (dU/dy) = mu * (d2U/dy2) – (dP/dx) + g_x
2. Persamaan Momentum Vertikal (Sumbu-y):u * (dV/dx) + v * (dV/dy) = mu * (d2V/dx2) – (dP/dx) + g_y

Penurunan persamaan di atas membuktikan bahwa volume air di sekitar struktur jembatan mengalami perubahan gradien tekanan (dP/dx) dan gaya geser kinematik viskositas (mu) yang sangat ekstrem dalam hitungan milidetik saat tumbukan berlangsung. Hal ini memperkuat landasan teoretis penelitian, menunjukkan bahwa penyederhanaan kurva setengah sinus pada metode numerik manual yang saya susun di spreadsheet didukung secara ilmiah oleh prinsip mekanika kontinu fluida yang kompleks.

BAB V: KESIMPULAN DAN REKOMENDASI

Kesimpulan

Berdasarkan seluruh hasil perhitungan numerik menggunakan Metode Integrasi Simpson 1/3 yang telah difinalisasi, dapat disimpulkan bahwa interaksi antara tongkang batubara berkapasitas dinamis 8.400.000 kg dengan pilar Jembatan Mahakam pada kecepatan 2 m/s menghasilkan transfer energi kejut yang sangat masif, dengan nilai total Impuls tumbukan sebesar 21.017.526 N.s dalam durasi rentang waktu 1,2 detik. Pola distribusi gaya beban dinamis terbukti membentuk kurva non-linear yang mencapai titik kritis gaya puncak (Peak Force) sebesar 27.500.000 Newton tepat pada detik ke-0,6.

Rekomendasi Teknis

Luaran kuantitatif dari penelitian berbasis metode numerik ini sukses dikembangkan menjadi rekomendasi rekayasa praktis yang nyata. Pertama, hasil estimasi gaya impak ini dapat dijadikan acuan otoritas kesyahbandaran lokal untuk menetapkan zonasi batas kecepatan aman (safe speed) kaku maksimal sebesar 1,5 knot bagi tongkang jumbo yang melintasi kolong jembatan saat kondisi arus sungai deras. Kedua, angka beban puncak sebesar 27,5 Meganewton ini dapat langsung digunakan oleh para perancang struktur sebagai parameter input batas bawah (baseline data) di dalam mendesain sistem struktur pelindung pilar (fendering system) maupun tiang pancang pengarah (dolphin deflector) yang lebih kokoh, guna menjamin proteksi jangka panjang terhadap pilar utama Jembatan Mahakam dari ancaman kegagalan katastrofik di masa depan.

Refleksi Filosofis: Integrasi Perhitungan Numerik dengan Siklus DAI5

Seluruh rangkaian formulasi matematika dan eksekusi numerik di atas sejatinya merupakan manifestasi konkret dari penerapan siklus DAI5 (Deep Awareness of I, Amplify Goals, Intention, Initiative, Instruction), sebuah kerangka berpikir terstruktur yang menuntun jalannya studi kasus ini dari sebuah keresahan personal hingga menjadi solusi rekayasa yang dapat dipertanggungjawabkan.

Perjalanan ini berakar dari tahap Deep Awareness of I, di mana muncul kesadaran mendalam dan keresahan pribadi dalam diri saya sebagai mahasiswa teknik perkapalan yang berasal dari Samarinda. Setiap kali mendengar kabar pilar Jembatan Mahakam kembali dihantam oleh tongkang batubara, nurani saya terusik. Saya menyadari bahwa keilmuan yang saya tekuni di bangku kuliah harus memiliki kebermanfaatan nyata bagi kampung halaman saya, dan dari kesadaran moral itulah dorongan untuk bertindak ini lahir. Kesadaran ini kemudian menuntun saya ke tahap Amplify Goals, di mana saya memperluas dan memperkuat visi jangka panjang dari sekadar menyelesaikan tugas mata kuliah Metode Numerik menjadi sebuah misi optimasi keselamatan maritim dan perlindungan infrastruktur daerah. Saya menetapkan target yang lebih tinggi, yaitu menyajikan sebuah model estimasi gaya benturan yang akurat secara ilmiah agar dapat menjadi landasan mitigasi risiko yang riil di perairan Sungai Mahakam.

Untuk mewujudkan tujuan besar tersebut, saya memasuki tahap Intention, di mana saya menetapkan niat kuat dan memformulasikan strategi pendekatan masalah. Di tahap inilah keputusan teknis krusial diambil, yaitu memilih Metode Integrasi Simpson 1/3 sebagai computational tools utama untuk memecahkan karakteristik beban tumbukan dinamis yang bersifat non-linear dan berubah drastis terhadap waktu. Niat ini bukan sekadar rencana, melainkan pilihan metodologi yang rasional untuk mengejar akurasi tinggi. Komitmen tersebut langsung saya tindak lanjuti pada tahap Initiative, di mana saya secara mandiri mulai mengumpulkan data parameter fisik sekunder—seperti estimasi bobot tongkang 7.000 ton, penentuan faktor added mass sebesar 20%, prediksi kecepatan arus sungai, hingga pemodelan kurva gaya dengan fungsi half-sine pulse. Saya mengambil langkah aktif untuk mendiskretisasi domain waktu menjadi grid-grid numerik yang terstruktur di spreadsheet sebelum masa tenggat tugas mendekat.

Akhirnya, seluruh akumulasi proses berpikir tersebut bermuara pada tahap Instruction, yang ditunjukkan melalui eksekusi teknis yang dijabarkan dalam detail hitungan di atas. Tahap ini adalah momen di mana rumus Simpson benar-benar dioperasikan langkah demi langkah, mengalikan gaya sesaat dengan koefisien bobot (1, 4, 2, 4, 2, 4, 1), melakukan akumulasi angka total 315.262.892, hingga menghasilkan nilai akhir impuls sebesar 21.017.526 N.s. Melalui siklus DAI5 yang utuh ini, proses belajar metode numerik tidak lagi berjalan sebagai rutinitas akademik yang hampa, melainkan bertransformasi menjadi sebuah kerja rekayasa yang terarah, berlandaskan empati, dan menghasilkan keluaran kuantitatif yang valid demi keselamatan infrastruktur publik.