ccitonline.com

A. Project Title

Analisis Numerik Getaran pada Sistem Suspensi Otomotif: Pendekatan Metode Numerik Berbasis DAI5

B. Author Complete Name

Herdi Agusta

C. Affiliation

Fakultas Teknik, Universitas Indonesia

D. Abstract

Makalah ini membahas analisis getaran pada sistem suspensi otomotif menggunakan metode numerik berbasis kerangka DAI5. Tujuan utama proyek ini adalah memahami karakteristik dinamis sistem suspensi kendaraan dengan menggunakan pendekatan numerik untuk menyelesaikan persamaan gerak. Dengan mengadopsi model sistem massa-pegas-peredam (single degree of freedom) dan menyelesaikannya menggunakan metode Euler dan Runge-Kutta, diperoleh respons sistem terhadap gangguan eksternal. Studi ini menyoroti pentingnya kesadaran mendalam, niat etis, idealisasi yang realistis, dan pendekatan terstruktur dalam memecahkan permasalahan teknik mesin modern.

E. Author Declaration

1. Deep Awareness (of) I

Dalam mengerjakan proyek ini, saya senantiasa mengingat peran diri saya sebagai makhluk ciptaan Tuhan Yang Maha Esa, menyadari keterbatasan dan keharusan untuk bertanggung jawab dalam menghasilkan solusi yang tidak hanya efektif secara teknis, tetapi juga bermanfaat dan bernilai etis.

2. Intention of the Project Activity

Proyek ini bertujuan untuk mengembangkan model numerik yang akurat dan efisien untuk memahami perilaku getaran sistem suspensi otomotif, dengan niat untuk meningkatkan kenyamanan dan keselamatan berkendara serta berkontribusi terhadap perkembangan ilmu teknik yang berkelanjutan.

F. Introduction

Sistem suspensi memainkan peranan penting dalam kendaraan bermotor, bertugas meredam getaran akibat ketidakteraturan permukaan jalan dan mempertahankan stabilitas kendaraan. Pemahaman perilaku dinamis dari sistem ini menjadi kunci untuk merancang kendaraan yang aman dan nyaman.

Initial Thinking (about the Problem):

• Analisis Masalah Secara Sistematis: Sistem suspensi mengalami getaran yang kompleks akibat interaksi antara kendaraan dan permukaan jalan. Getaran ini harus dianalisis untuk menentukan respons dinamis kendaraan.
• Soroti Penelitian Sebelumnya dan Kesenjangan: Studi-studi terdahulu menunjukkan bahwa model matematis sistem suspensi sederhana dapat memberikan gambaran yang cukup akurat, namun pemecahan numeriknya perlu metode yang tepat agar dapat diterapkan dalam berbagai kondisi.
• Mengurai Masalah: Fokus permasalahan diarahkan pada perhitungan respon dinamis gaya luar sederhana terhadap sistem massa-pegas-peredam menggunakan metode numerik.
• Dekonstruksi Prinsip Dasar: Persamaan gerak sistem dapat diperoleh dari hukum Newton II dan diturunkan menjadi persamaan diferensial orde kedua.
• Analisis State-of-the-Art: Teknologi saat ini menggunakan metode numerik canggih dalam simulasi suspensi, namun pendekatan berbasis prinsip idealisasi sederhana masih sangat relevan untuk pembelajaran dan dasar pengembangan lebih lanjut.

G. Methods & Procedures

Idealization:

• Sistem diasumsikan sebagai model single degree of freedom (massa-pegas-peredam) dengan gaya eksitasi sinusoidal.
• Asumsi:
  • Gerak hanya terjadi secara vertikal.
  • Sistem linear, tidak ada efek nonlinier.
  • Pengaruh suhu dan perubahan sifat material diabaikan.

Instruction (Set):

1. Menyusun persamaan gerak: mx¨+cx˙+kx=F(t)m \ddot{x} + c \dot{x} + kx = F(t).
2. Menyusun ulang menjadi sistem persamaan diferensial orde pertama.
3. Implementasi metode numerik:
   • Metode Euler.
   • Metode Runge-Kutta orde 4.
4. Simulasi pada variasi parameter:
   • Massa (m)
   • Konstanta pegas (k)
   • Koefisien redaman (c)
5. Visualisasi hasil:
   • Plot respon displacement terhadap waktu.
6. Analisis performa metode numerik berdasarkan akurasi dan kestabilan.

H. Results & Discussion

Setelah simulasi dilakukan:

• Metode Runge-Kutta menunjukkan hasil yang jauh lebih stabil dan akurat dibandingkan metode Euler, terutama untuk langkah waktu yang lebih besar.
• Sistem dengan redaman rendah menunjukkan respons osilasi lebih panjang dibandingkan sistem dengan redaman tinggi.
• Massa kendaraan mempengaruhi amplitudo respon: semakin besar massa, semakin kecil percepatan respon.

Diskusi:

• Keterbatasan metode Euler yang menghasilkan akumulasi error signifikan.
• Pentingnya pemilihan metode numerik sesuai dengan kebutuhan akurasi dan efisiensi komputasi.
• Simulasi ini berguna untuk pengembangan lebih lanjut seperti kontrol aktif suspensi.

I. Conclusion, Closing Remarks, Recommendations

Studi ini membuktikan bahwa metode numerik sangat efektif untuk menganalisis perilaku dinamis sistem suspensi sederhana. Dengan pendekatan DAI5, tidak hanya solusi teknis yang dihasilkan, namun juga kesadaran akan nilai-nilai etis dan keberlanjutan. Direkomendasikan untuk penelitian lanjutan menggunakan model nonlinier dan studi eksperimen untuk validasi model.

J. Acknowledgments

Penulis mengucapkan terima kasih kepada dosen pengampu mata kuliah Metode Numerik, teman-teman sekelas, dan seluruh pihak yang telah mendukung dalam penyelesaian proyek ini.

K. References

1. Rao, S.S. (2011). Mechanical Vibrations. 5th ed. Pearson Education.
2. Chapra, S.C., Canale, R.P. (2015). Numerical Methods for Engineers. 7th ed. McGraw-Hill.
3. Norton, R.L. (2013). Design of Machinery. McGraw-Hill.

L. Appendices

Appendix A. Kode Simulasi Metode Runge-Kutta

import numpy as np
import matplotlib.pyplot as plt

# Parameter sistem
m = 250  # massa (kg)
k = 15000  # konstanta pegas (N/m)
c = 1200  # konstanta redaman (Ns/m)
F0 = 100  # amplitudo gaya luar (N)
omega = 10  # frekuensi gaya luar (rad/s)

# Persamaan gerak dalam bentuk sistem orde 1
def model(t, y):
    x, v = y
    dxdt = v
    dvdt = (F0*np.sin(omega*t) - c*v - k*x)/m
    return np.array([dxdt, dvdt])

# Metode Runge-Kutta orde 4
def runge_kutta(f, y0, t):
    n = len(t)
    y = np.zeros((n, len(y0)))
    y[0] = y0
    for i in range(n-1):
        dt = t[i+1] - t[i]
        k1 = f(t[i], y[i])
        k2 = f(t[i] + dt/2, y[i] + dt*k1/2)
        k3 = f(t[i] + dt/2, y[i] + dt*k2/2)
        k4 = f(t[i] + dt, y[i] + dt*k3)
        y[i+1] = y[i] + dt*(k1 + 2*k2 + 2*k3 + k4)/6
    return y

# Setup waktu
T = 10  # waktu simulasi (s)
N = 1000  # jumlah langkah
t = np.linspace(0, T, N)

# Inisialisasi kondisi awal
y0 = [0, 0]  # x(0) = 0, v(0) = 0

# Simulasi
solusi = runge_kutta(model, y0, t)

# Plot hasil
plt.figure(figsize=(10,6))
plt.plot(t, solusi[:,0])
plt.title('Respon Displacement Sistem Suspensi')
plt.xlabel('Waktu (s)')
plt.ylabel('Displacement (m)')
plt.grid(True)
plt.show()