ccitonline.com

CCIT – Cara Cerdas Ingat Tuhan

| AI-DAI5 | DAI5 AI Agents | NIC | ZWI | | CCITEdu | DAI5 eBook | CFDSOF | Donation | Download | CCIT Corporation | DAI5 | 33 Kriteria Evaluasi Penerapan DAI5 | Search |

Tugas Besar- Pemodelan Dinamika Beban Kendaraan terhadap Jalan dengan Quarter-Car 2 Degree of Freedom Model

ุจูุณู’ู…ู ูฑู„ู„ูŽู‘ูฐู‡ู ูฑู„ุฑูŽู‘ุญู’ู…ูŽูฐู†ู ูฑู„ุฑูŽู‘ุญููŠู…ู

Pemodelan dinamika kendaraan-jalan melalui pendekatan Quarter-Car 2-DOF ini memberikan kerangka analitis yang presisi untuk menghitung beban yang diterima infrastruktur jalan. Beban kendaraan disederhanakan menjadi model Quarter-Car 2-Degree-of-Freedom (2-DOF). Artinya, kita hanya melihat seperempat bagian mobil (satu roda beserta suspensi dan bodi di atasnya) yang bergerak dalam dua derajat kebebasan (gerak vertikal bodi dan gerak vertikal roda).

Perumusan Persamaan Gerak (Hukum Newton II dan Hukum Hooke)

Model Quarter-Car 2-DOF membagi kendaraan menjadi dua massa utama:

  • Massa terpegas (Massa bodi kendaraan) : m_s
  • Massa tak terpegas ( Massa poros/roda) : m_u

Sistem ini dihubungkan oleh pegas dan peredam (damper). Berdasarkan Hukum Hooke, gaya pemulih suspensi proporsional dengan perpindahan, dan untuk peredam proporsional dengan kecepatan:

F=kฮ”zF = k \Delta z
F=cฮ”zห™F = c \Delta \dot{z}

Menggunakan Hukum Newton II, kita dapatkan Persamaan Gerak (Equation of Motion):

โˆ‘F=ma\sum F = ma

Untuk Massa Terpegas (Bodi, m_s):

mszยจs=โˆ’ks(zsโˆ’zu)โˆ’cs(zห™sโˆ’zห™u)m_s \ddot{z}_s = -k_s(z_s – z_u) – c_s(\dot{z}_s – \dot{z}_u)

Untuk Massa Tak Terpegas (Poros/Roda, m_u):

muzยจu=ks(zsโˆ’zu)+cs(zห™sโˆ’zห™u)โˆ’kt(zuโˆ’zr)โˆ’ct(zห™uโˆ’zห™r)m_u \ddot{z}_u = k_s(z_s – z_u) + c_s(\dot{z}_s – \dot{z}_u) – k_t(z_u – z_r) – c_t(\dot{z}_u – \dot{z}_r)

Keterangan:

zs,zu,zrz_s, z_u, z_r
zห™s,zห™u,zห™r\dot{z}_s, \dot{z}_u, \dot{z}_r
zยจs,zยจu\ddot{z}_s, \ddot{z}_u
ks,csk_s, c_s
kt,ct k_t, c_t
  • Perpindahan/posisi vertikal dari bodi, roda, dan profil jalan.
  • Kecepatan vertikal bodi, roda, dan turunan profil jalan.
  • Percepatan vertikal bodi dan roda.
  • Kekakuan dan koefisien redaman suspensi.
  • Kekakuan dan koefisien redaman ban (tire).

Beban Dinamis & Energi pada Jalan (Hukum Newton III)

Berdasarkan Hukum Newton III (Aksi-Reaksi), gaya dinamis (Dynamic Load) yang diberikan oleh roda ke jalan adalah sama dengan gaya kontak dari jalan ke ban. Gaya aksi ini dirumuskan sebagai:

Fjalan(t)=kt(zu(t)โˆ’zr(t))+ct(zห™u(t)โˆ’zห™r(t))F_{jalan}(t) = k_t \big(z_u(t) – z_r(t)\big) + c_t \big(\dot{z}_u(t) – \dot{z}_r(t)\big)

Energi yang Diterima Jalan: Jika Anda ingin menghitung energi ($E$) yang benar-benar ditransfer ke jalan (misalnya untuk analisis kelelahan/kerusakan aspal), ini dihitung dari kerja mekanik. Kerja adalah integral dari gaya terhadap jarak deformasi/defleksi perkerasan jalan.

Ejalan=โˆซ0TFjalan(t)โ‹…ฮดห™r(t)dtE_{jalan} = \int_{0}^{T} F_{jalan}(t) \cdot \dot{\delta}_r(t) \, dt

Reduksi ke Persamaan Diferensial Orde 1 (State-Space)

Metode Runge-Kutta membutuhkan persamaan diferensial orde 1 (ODE orde 1). Kita harus mereduksi dua persamaan ODE orde 2 di atas menjadi empat ODE orde 1. Definisikan variabel state (kondisi) Y:

y1=zs(PosisiBodi)y_1 = z_s (Posisi Bodi)
y2=zห™s(KecepatanBodi)y_2 = \dot{z}_s (Kecepatan Bodi)
y3=zu(PosisiRoda)y_3 = z_u (Posisi Roda)
y4=zห™u(KecepatanRoda)y_4 = \dot{z}_u (Kecepatan Roda)

Vektor turunannya menjadi:

Yห™=f(t,Y)\dot{Y} = f(t, Y)
yห™1=y2\dot{y}_1 = y_2
yห™2=1ms[โˆ’ks(y1โˆ’y3)โˆ’cs(y2โˆ’y4)]\dot{y}_2 = \frac{1}{m_s} \big[ -k_s(y_1 – y_3) – c_s(y_2 – y_4) \big]
yห™3=y4\dot{y}_3 = y_4
yห™4=1mu[ks(y1โˆ’y3)+cs(y2โˆ’y4)โˆ’kt(y3โˆ’zr)โˆ’ct(y4โˆ’zห™r)]\dot{y}_4 = \frac{1}{m_u} \big[ k_s(y_1 – y_3) + c_s(y_2 – y_4) – k_t(y_3 – z_r) – c_t(y_4 – \dot{z}_r) \big]

Penyelesaian Numerik dengan Runge-Kutta Orde 4 (RK4)

Untuk mensimulasikan sistem sepanjang waktu simulasi dengan langkah waktu (time-step) sebesar h, algoritma RK4 menghitung matriks state pada waktu berikutnya (Y_{n+1}) berdasarkan 4 gradien evaluasi (k_1, k_2, k_3, k_4):

Langkah iterasi dalam komputasi:

k1=f(tn,Yn)(Kemiringandiawalinterval)k_1 = f(t_n, Y_n) (Kemiringan di awal interval)
k2=f(tn+h2,Yn+h2k1)(Kemiringandititiktengahberdasarkank1)k_2 = f\left(t_n + \frac{h}{2}, Y_n + \frac{h}{2} k_1\right) (Kemiringan di titik tengah berdasarkan k_1)
k3=f(tn+h2,Yn+h2k2)(Kemiringandititiktengahberdasarkank2)k_3 = f\left(t_n + \frac{h}{2}, Y_n + \frac{h}{2} k_2\right) (Kemiringan di titik tengah berdasarkan k_2)
k4=f(tn+h,Yn+hk3)(Kemiringandiakhirintervalberdasarkank3)k_4 = f(t_n + h, Y_n + h k_3) (Kemiringan di akhir interval berdasarkan k_3)

Pembaharuan Nilai (Update State):

Yn+1=Yn+h6(k1+2k2+2k3+k4)Y_{n+1} = Y_n + \frac{h}{6} (k_1 + 2k_2 + 2k_3 + k_4)

Leave a Reply

Your email address will not be published. Required fields are marked *