Untuk membedah secara ilmiah fenomena struktural yang terjadi di lapangan, saya menyusun alur perhitungan mekanika tumbukan yang terjadi antara struktur apung dan struktur tegar. Langkah awal yang paling krusial sebelum menyentuh persamaan numerik adalah mendefinisikan karakteristik fisik dari objek yang bergerak. Ketika sebuah tongkang batubara berukuran standar 300 kaki bergerak menyusuri perairan Sungai Mahakam dengan muatan penuh, bobot mati atau massa kapal diestimasikan berada pada angka kisaran 7.000 ton atau setara dengan 7.000.000 kg. Namun, di dalam dunia teknik perkapalan dan hidrodinamika, perhitungan momentum tidak boleh hanya mengandalkan massa kering kapal itu sendiri. Ketika lambung kapal membelah air dan melaju, ada sejumlah volume fluida di sekeliling lambung yang ikut menempel, bergerak, dan terdorong searah dengan pergerakan kapal. Massa fluida yang ikut terakselerasi ini disebut sebagai Massa Tambah atau Added Mass. Berdasarkan literatur mekanika fluida perkapalan untuk bentuk lambung kapal tongkang yang cenderung kotak dengan nilai block coefficient tinggi, nilai koefisien massa tambah diestimasikan sebesar 20% dari total massa aktualnya. Jika faktor ini diabaikan, maka gaya benturan yang dihasilkan akan jauh lebih kecil dari kondisi riil di lapangan, sehingga berisiko menciptakan kegagalan desain pelindung pilar jembatan. Melalui pertimbangan hidrodinamika tersebut, perhitungan total massa dinamis (M) yang terlibat langsung dalam proses transfer momentum dijabarkan melalui penjumlahan M = M_tongkang + M_added, yang menghasilkan angka 7.000.000 kg + 1.400.000 kg = 8.400.000 kg. Tongkang ini diasumsikan hanyut atau melaju akibat kombinasi tenaga mesin dan dorongan kecepatan arus sungai yang deras, sehingga menghasilkan kecepatan total (v) sebesar 2 m/s atau kurang lebih sekitar 4 knot.
Ketika lambung baja tongkang yang tebal melakukan kontak pertama kali dengan struktur beton pilar jembatan yang masif, terjadi proses penghentian momentum secara mendadak. Berdasarkan Hukum Newton II (F = m * delta_v / delta_t) dan Hukum Momentum-Impuls, seluruh energi kinetik dari massa dinamis sebesar 8.400.000 kg tersebut harus diredam hingga kecepanannya turun menjadi nol dalam rentang waktu yang sangat singkat. Durasi waktu total persentuhan dari awal tumbukan hingga kapal berhenti total diestimasikan berjalan selama delta_t = 1,2 detik. Secara teoritis, total efek hantaman ini didefinisikan sebagai Impuls, yang merupakan luasan daerah di bawah kurva gaya terhadap waktu dengan persamaan dasar:
Impuls = Integral dari t=0 sampai t=1,2 untuk fungsi F(t) dt
Tantangan utama rekayasa ini muncul karena gaya benturan F(t) tidak mungkin konstan. Karakteristik gaya benturan dinamis akan mengikuti fase deformasi material, di mana gaya bernilai nol pada saat awal sentuhan (t = 0), kemudian merangkak naik secara ekstrem seiring rusaknya material lambung and pilar hingga mencapai titik beban puncak tertinggi (F_max) tepat di tengah-tengah durasi benturan (t = 0,6 detik), dan akhirnya meluncur turun kembali ke angka nol ketika seluruh energi kinetik kapal telah habis terserap dan kapal berhenti bergerak (t = 1,2 detik). Untuk merepresentasikan perilaku fisik non-linear yang melengkung tajam ini secara matematis ke dalam perhitungan, saya menggunakan pendekatan fungsi setengah sinus dengan rumus gaya dinamis F(t) = F_max * sin((pi * t) / 1,2).
Karena persamaan fungsi sinus tersebut melibatkan variabel gaya yang terus berubah terhadap waktu secara kontinu, penggunaan matematika analitik biasa akan memakan waktu dan sulit dikembangkan jika nanti ada variasi data lapangan yang diskrit. Di sinilah saya menerapkan Metode Integrasi Numerik Simpson 1/3 untuk memecahkan masalah tersebut. Langkah pertama metode numerik ini adalah melakukan diskretisasi domain waktu, di mana durasi total benturan selama 1,2 detik saya bagi menjadi n = 6 segmen interval yang sama besar. Syarat wajib yang harus dipenuhi dalam metode Simpson 1/3 adalah jumlah segmen harus bernilai genap. Jarak antartitik data atau lebar langkah interval (h) dihitung melalui rumus h = (t_akhir – t_awal) / n, sehingga didapatkan nilai h = (1,2 – 0) / 6 = 0,2 detik. Pembagian domain ini berhasil mendirikan 7 buah titik koordinat waktu diskrit atau grid yang berurutan, dimulai dari t_0 = 0,0 detik, t_1 = 0,2 detik, t_2 = 0,4 detik, t_3 = 0,6 detik, t_4 = 0,8 detik, t_5 = 1,0 detik, hingga t_6 = 1,2 detik.
Dengan menetapkan nilai estimasi gaya puncak batas atas sebesar F_max = 27.500.000 Newton berdasarkan pendekatan konversi energi kinetik total, saya menghitung nilai gaya sesaat F(t) pada masing-masing dari ke-7 titik grid waktu tersebut. Aturan Simpson 1/3 bekerja secara cerdas dengan memberikan bobot pengali atau koefisien spesifik pada deret datanya, yaitu mengikuti pola koefisien (1, 4, 2, 4, 2, 4, 1). Pola bobot ini didasarkan pada integrasi polinomial kuadratik atau pendekatan kurva parabola yang jauh lebih akurat mengikuti lekukan grafik sinus dibandingkan metode trapesium linier biasa. Pada titik awal kontak yaitu t_0 = 0,0 detik, gaya sesaat dihitung melalui rumus F(0,0) = 27.500.000 * sin((3,14159 * 0,0) / 1,2) yang menghasilkan 0 Newton, lalu dikalikan dengan bobot pertama sehingga hasilnya tetap 0. Pada titik penetrasi awal yaitu t_1 = 0,2 detik, nilai gayanya adalah F(0,2) = 27.500.000 * sin((3,14159 * 0,2) / 1,2) yang menghasilkan 13.750.000 Newton setelah dikonversi ke sudut 30 derajat atau 0,5236 radian, kemudian nilai ini dikalikan dengan bobot kedua menjadi 55.000.000. Memasuki fase kenaikan gaya deformasi pada t_2 = 0,4 detik, nilai gayanya meningkat menjadi F(0,4) = 27.500.000 * sin((3,14159 * 0,4) / 1,2) atau setara dengan 23.815.723 Newton, yang setelah dikalikan bobot ketiga menghasilkan angka 47.631.446. Tepat pada titik beban impak maksimum atau peak force di waktu t_3 = 0,6 detik, kurva mencapai puncaknya di mana F(0,6) = 27.500.000 * sin((3,14159 * 0,6) / 1,2) menghasilkan nilai penuh 27.500.000 Newton, dan setelah dikalikan bobot keempat menghasilkan nilai sebesar 110.000.000. Setelah melewati titik puncak, gaya memasuki fase redaman atau pelepasan energi pada t_4 = 0,8 detik dengan nilai F(0,8) = 27.500.000 * sin((3,14159 * 0,8) / 1,2) yang kembali ke angka 23.815.723 Newton, lalu dikalikan bobot kelima menghasilkan 47.631.446. Mendekati fase penghentian kapal pada t_5 = 1,0 detik, nilainya menurun menjadi F(1,0) = 27.500.000 * sin((3,14159 * 1,0) / 1,2) yang menghasilkan 13.750.000 Newton, lalu dikalikan bobot keenam menjadi 55.000.000. Pada akhirnya, saat kapal berhenti sempurna di titik t_6 = 1,2 detik, nilai F(1,2) = 27.500.000 * sin((3,14159 * 1,2) / 1,2) kembali menyentuh angka 0 Newton, dan hasil perkalian bobot terakhirnya adalah 0. Setelah mengeksekusi seluruh baris perhitungan tersebut menggunakan bantuan rumus di dalam spreadsheet, saya melakukan akumulasi atau penjumlahan total terhadap seluruh hasil perkalian bobot dari titik awal hingga akhir, yang menghasilkan nilai akumulatif TOTAL = 0 + 55.000.000 + 47.631.446 + 110.000.000 + 47.631.446 + 55.000.000 + 0, sehingga didapatkan angka mutlak sebesar 315.262.892.
Tahap akhir dari pengoperasian Algoritma Simpson 1/3 ini adalah melakukan finalisasi perhitungan dengan mengalikan angka TOTAL akumulasi deret tersebut dengan faktor pengali sepertiga dari lebar interval waktu atau h/3. Langkah akhir ini akan memberikan nilai luasan definitif di bawah kurva gaya yang merepresentasikan parameter total Impuls Tumbukan secara presisi. Melalui substitusi angka ke dalam rumus Impuls = (h / 3) * TOTAL, perhitungan berjalan menjadi Impuls = (0,2 / 3) * 315.262.892, yang setelah disederhanakan menggunakan perkalian desimal 0,06666667 * 315.262.892 menghasilkan nilai akhir Impuls sebesar 21.017.526 N.s.
Pengembangan Akademik: Pengenalan Ekosistem CAE dan Dinamika Fluida Komputasi
Proses pemahaman saya tidak berhenti pada perhitungan numerik manual berbasis spreadsheet saja. Pada tahap ini, saya memperluas wawasan rekayasa dengan mempelajari bagaimana industri maritim modern menyelesaikan masalah hidrodinamika yang jauh lebih kompleks menggunakan perangkat lunak Dinamika Fluida Komputasi atau CFDSOF yang berbasis metode numerik volume hingga (Finite Volume Method). Di dalam perkuliahan, saya diajarkan bahwa proses rekayasa kapal modern terintegrasi ke dalam tiga pilar utama, yaitu Computer Aided Design (CAD) untuk pemodelan bentuk geometris tiga dimensi dari tongkang dan jembatan, Computer Aided Engineering (CAE) untuk melakukan analisis kekuatan struktur serta simulasi aliran fluida, dan Computer Aided Manufacture (CAM) untuk menerjemahkan hasil desain menjadi proses produksi di galangan. Kasus benturan pilar jembatan ini berada penuh di dalam domain CAE, di mana pendekatan fisika dinamis tidak lagi hanya diselesaikan dengan Hukum Newton linier sederhana, melainkan harus diturunkan ke dalam persamaan diferensial parsial (Partial Differential Equations / PDE) yang mengatur pergerakan fluida di sekitar lambung kapal.
Melalui pemahaman fenomena hidrodinamika ini, pergerakan air di sekitar pilar akibat dorongan kapal dapat dimodelkan secara matematis ke dalam ruang dua dimensi (2D) menggunakan prinsip kekekalan momentum fluida. Persamaan dinamis tersebut diturunkan untuk komponen kecepatan arah sumbu-x (u) dan komponen kecepatan arah sumbu-y (v) secara simultan. Untuk pergerakan fluida pada arah horizontal atau sumbu-x, gaya-gaya yang bekerja (F_x) dipengaruhi oleh interaksi kecepatan fluida, viskositas, gradien tekanan, dan gaya gravitasi lokal, yang dirumuskan melalui persamaan diferensial parsial:
u * (dU/dx) + v * (dU/dy) = mu * (d2U/dy2) – (dP/dx) + g_x
Secara analogi, pergerakan fluida pada arah vertikal atau sumbu-y (F_y) dikontrol oleh kombinasi gaya dinamis yang serupa pada komponen kecepatan v, dengan rumusan:
u * (dV/dx) + v * (dV/dy) = mu * (d2V/dx2) – (dP/dx) + g_y
Pembelajaran mengenai persamaan diferensial parsial ini membuka mata saya bahwa volume air yang terjepit di antara lambung tongkang dan pilar jembatan mengalami perubahan tekanan (dP) dan geseran viskositas (mu) yang sangat ekstrem dalam elemen volume yang sangat kecil (infinitesimal). Integrasi teori matematika tingkat lanjut ini memperkuat validitas analisis saya, karena menunjukkan bahwa di balik penyederhanaan metode integrasi Simpson yang saya lakukan di spreadsheet, terdapat kompleksitas pergerakan fluida riil yang diatur oleh hukum-hukum mekanika kontinu.
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 serta pemahaman mendalam mengenai integrasi CAD-CAE menggunakan persamaan Navier-Stokes 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.
