Pemodelan Dinamis & Penurunan Persamaan Diferensial Orde Dua

Kajian sistemik aktuator pneumatik kompak bersuspensi aerostatik (ASCPA) untuk penggerindaan robotik presisi: eliminasi gesekan nonlinier sebesar 99,48%, formulasi relasi fisis masukan-keluaran, hingga reduksi matematis menjadi 2 buah persamaan diferensial orde 1 (State-Space).

tune Uji Parameter Interaktif developer_board Skema Fisik & I/O
Reduksi Gesekan 99.48%
Frek. Alami \(\omega_n = 5.0\text{ rad/s}\)
Kutub LHP \(p_{1,2} = -3 \pm j4\)

speed Detail Status Telemetri Aktuator ASCPA

ASCPA-CYL-01
air

Piston Film Aerostatik (\(D_p = 20\,\text{mm}\))

NORMAL OPERATION check_circle AI: Kontrol Linier Valid
Gaya Masukan \(u(t)\) 1.00 N
Perpindahan \(y(t)\) 0.040 m
Tebal Film Udara \(h_0 = 30\,\mu\text{m}\)
Gesekan Coulomb 0.0777 N (-99.48%)
Tekanan Suplai 0.40 MPa
Status Kontak ELASTIS LINIER
Karakteristik LTI Orde 2 KANONIK
Massa \(m\) 1.0 kg
Redaman \(c\) 6.0 N·s/m
Kekakuan \(k\) 25.0 N/m
Rasio redaman \(\zeta = 0.60\) (Underdamped). Respon getaran teredam stabil tanpa fenomena stick-slip.
psychology Asisten State-Space EM-KF-RTS
Reduksi ke 2 Persamaan Orde 1
\(\dot{x}_1 = x_2, \quad \dot{x}_2 = -25x_1 - 6x_2 + u\)

Nilai eigen matriks keadaan: \(\lambda_{1,2} = -3 \pm j4\).

Sistem LTI memenuhi kriteria BIBO & Asimtotik stabil secara matematis.

Overview Paper Internasional Bereputasi (MSSP 2026)

Kajian jurnal ilmiah non-buku teks bereputasi tinggi terindeks Scopus Q1 & Web of Science.

precision_manufacturing

Permesinan Presisi & Kontrol Gaya

Pada manufaktur bilah turbin mesin aero dan lensa optik, end-effector robot membutuhkan aktuator fleksibel yang mampu menjaga gaya kontak normal secara stabil guna mencegah terjadinya retak mikro bawah-permukaan.

error_outline

Kelemahan TPA/LFPA (Stick-Slip)

Aktuator konvensional menggunakan segel karet luncur yang memicu gesekan Coulomb masif (\(15\text{ N}\)). Ini menghasilkan histeresis, osilasi limit cycles, dan perilaku nonlinier parah yang merusak kehalusan penggerindaan.

air

Inovasi Suspensi Aerostatik ASCPA

Dengan menyuplai udara bertekanan hidrostatik (\(P = 0.4\,\text{MPa}\)) melalui \(n=10\) orifis ke celah \(h_0 = 30\,\mu\text{m}\), kontak padat ditiadakan. Gesekan Coulomb turun 99.48% menjadi hanya \(0.0777\text{ N}\), mewujudkan sistem linier murni!

table_chart Tabel Identifikasi Parameter Gesekan Riil Eksperimental (Tabel 2 Paper Yang et al., 2026)

Model LuGre: \(f = (F_c + (F_s - F_c)e^{-(v/v_s)^2})\text{sign}(v) + \sigma_2 v\)
Parameter Model LuGre TPA (Konvensional) LFPA (Low Friction) ASCPA (Suspensi Udara) Tingkat Reduksi
Gesekan Coulomb \(F_c\) 15.00 N 5.40 N 0.0777 N Turun 99.48%
Gesekan Statis Maksimum \(F_s\) 18.45 N 8.70 N 0.0767 N Turun 99.58%
Kecepatan Stribeck \(v_s\) 5.00 mm/s 1.40 mm/s 0.01 mm/s Bebas stick-slip
Koefisien Redaman Viskos \(\sigma_2\) 0.283 N·s/mm 0.010 N·s/mm 0.0012 N·s/mm (\(1.2\,\text{N}\cdot\text{s/m}\)) Dominan Linier

Penjelasan Relasi Fisis, Masukan (Input), dan Keluaran (Output)

Mekanika translasi 1-DOF: dari modulasi tekanan katup servo hingga perpindahan pahat & gaya kontak gerinda.

schema Skema Fisik End-Effector Gerinda ASCPA

Gbr. 1 Laporan
Skema Sistem ASCPA

Menunjukkan rakitan piston bersuspensi film udara hidrostatik tebal \(h_0 = 30\,\mu\text{m}\), batang piston bertumpu pada bantalan udara tanpa kontak, roda gerinda berputar, serta interaksi kontak elastis dengan benda kerja (\(k = 25\,\text{N/m}\)).

input Sinyal Masukan (Input) \(u(t)\) [Satuan: Newton (N)]

Gaya Dorong Pneumatik Bersih pada Piston

\[u(t) = A_p \cdot \Delta P(t) = A_p [P_1(t) - P_2(t)]\]

Asal Fisis: Perbedaan tekanan diferensial \(\Delta P(t)\) antara ruang maju dan mundur silinder yang dimodulasi secara dinamis oleh high-speed proportional servo valve (Festo MPYE-5-1/4-010-B) yang dikontrol tegangan perintah dari Real-Time Control Processor (RCP).

output Sinyal Keluaran (Output) \(y(t)\) [Satuan: Meter (m)]

Perpindahan Aksial Pahat Gerinda & Gaya Kontak Normal

\[y(t) = x(t) \quad \iff \quad F_c(t) = k_e \cdot y(t)\]

Asal Fisis: Perpindahan aksial \(x(t)\) dari rakitan spindel gerinda. Karena antarmuka kontak pahat dengan benda kerja memiliki konstanta kekakuan elastis \(k_e\), gaya penggerindaan normal \(F_c(t)\) berbanding lurus secara linier terhadap \(y(t)\).

Penurunan Lengkap Persamaan Diferensial Orde 2

Dari kesetimbangan translasi Hukum II Newton ke bentuk kanonik sistem getaran LTI.

01

Kesetimbangan Newton

\(m \ddot{y}(t) = \sum F\)

Jumlah gaya mencakup dorong pneumatik \(F_{\text{pneu}}\), gesekan \(F_{\text{fric}}\), dan kontak \(F_{\text{cont}}\).

02

Gaya Konstitutif

\(F_{\text{fric}} = c \dot{y}(t)\)
\(F_{\text{cont}} = k y(t)\)

Lapisan aerostatik mengubah gesekan padat Coulomb menjadi geser laminer linier murni.

03

PDB Orde 2 Gerak

\(m \ddot{y} + c \dot{y} + k y = u\)

Persamaan diferensial biasa (PDB) linier invarian waktu waktu-kontinu yang eksak.

04

Bentuk Kanonik

\(\ddot{y} + 2\zeta\omega_n \dot{y} + \omega_n^2 y = \frac{u}{m}\)

\(\omega_n = \sqrt{k/m} = 5.0\text{ rad/s}\), \(\zeta = \frac{c}{2\sqrt{km}} = 0.60\) (underdamped).

functions Fungsi Transfer Laplace \(G(s)\) & Konstelasi Kutub

Menerapkan Transformasi Laplace unilateral pada kondisi awal nol (\(y(0) = 0, \dot{y}(0) = 0\)):

\[G(s) = \frac{Y(s)}{U(s)} = \frac{1}{m s^2 + c s + k} = \frac{1}{s^2 + 6s + 25} = \frac{1}{(s + 3)^2 + 4^2}\]

Akar persamaan karakteristik menghasilkan sepasang kutub konjugat kompleks: \(p_{1,2} = -3 \pm j4\), dengan laju redaman eksponensial \(\sigma = 3.0\,\text{s}^{-1}\) dan frekuensi sudut teredam \(\omega_d = 4.0\,\text{rad/s}\).

Pemetaan Parameter Riil Paper
Massa Total \(m\) 1.0 kg
Redaman Total \(c\) 6.0 N·s/m
Kekakuan Kontak \(k\) 25.0 N/m
Penguatan DC \(K_{dc}\) 0.040 m/N

Bagaimana Cara Menjadikan PDB Orde 2 ke 2 Buah PDB Orde 1?

Penjelasan langkah demi langkah (step-by-step) reduksi orde ke dalam variabel keadaan (*State-Space Representation*).

help_outline Mengapa Persamaan Diferensial Orde 2 Harus Dipecah Menjadi 2 Buah Persamaan Orde 1?

1. Algoritma Filter Kalman (Paper MSSP) Algoritma EM-KF-RTS pada paper bekerja langsung pada persamaan ruang keadaan matriks orde 1: \(\mathbf{\dot{x}} = \mathbf{A}\mathbf{x} + \mathbf{B}u\).
2. Solver Komputasi Numerik (ODE45/RK4) Komputer tidak dapat mengintegrasikan turunan kedua secara langsung. Solver Runge-Kutta hanya dapat menyelesaikan vektor turunan pertama secara simultan.
3. Teori Kendali Modern (State Feedback) Memungkinkan kontrol posisi dan kecepatan sekaligus menggunakan hukum umpan balik keadaan \(u = -\mathbf{K}\mathbf{x}\).
A

Langkah 1: Definisikan Dua Variabel Keadaan Fisik (State Variables)

Karena sistem mekanik ini memiliki 2 elemen penyimpan energi (energi potensial pegas kontak dan energi kinetik massa), kita mendefinisikan 2 variabel keadaan yang mewakili posisi dan kecepatan:

Variabel Keadaan 1 (Posisi): \[x_1(t) \equiv y(t)\] Perpindahan aksial gerinda dari kesetimbangan [meter]
Variabel Keadaan 2 (Kecepatan): \[x_2(t) \equiv \frac{dy(t)}{dt} = \dot{y}(t)\] Kecepatan luncur aksial rakitan piston [m/s]
B

Langkah 2: Ambil Turunan Pertama Masing-Masing Variabel Keadaan

Cari laju perubahan waktu terhadap masing-masing variabel (\(\dot{x}_1\) dan \(\dot{x}_2\)):

Persamaan Diferensial Orde 1 Pertama:

Berdasarkan definisi kecepatan, turunan posisi adalah kecepatan:

\[\dot{x}_1(t) = \frac{dx_1}{dt} = \frac{dy}{dt} = x_2(t)\]
Persamaan Diferensial Orde 1 Kedua:

Turunan kecepatan adalah percepatan \(\ddot{y}(t)\). Ambil persamaan gerak Hukum Newton:

\[m \ddot{y}(t) + c \dot{y}(t) + k y(t) = u(t) \implies \ddot{y}(t) = -\frac{k}{m} y(t) - \frac{c}{m} \dot{y}(t) + \frac{1}{m} u(t)\]

Substitusikan \(y(t) = x_1(t)\) dan \(\dot{y}(t) = x_2(t)\):

\[\dot{x}_2(t) = -\frac{k}{m} x_1(t) - \frac{c}{m} x_2(t) + \frac{1}{m} u(t)\]
C

Langkah 3: Susun ke dalam Format Matriks Ruang Keadaan (State-Space)

Dua persamaan orde 1 simultan di atas disatukan ke dalam bentuk kanonik matriks:

\[\begin{bmatrix} \dot{x}_1(t) \\ \dot{x}_2(t) \end{bmatrix} = \begin{bmatrix} 0 & 1 \\ -\frac{k}{m} & -\frac{c}{m} \end{bmatrix} \begin{bmatrix} x_1(t) \\ x_2(t) \end{bmatrix} + \begin{bmatrix} 0 \\ \frac{1}{m} \end{bmatrix} u(t)\]

Persamaan keluaran sistem untuk perpindahan \(y(t) = x_1(t)\):

\[y(t) = \begin{bmatrix} 1 & 0 \end{bmatrix} \begin{bmatrix} x_1(t) \\ x_2(t) \end{bmatrix} + [0] u(t)\]
Substitusi Nilai Riil Paper (\(m = 1.0\,\text{kg}, c = 6.0\,\text{N}\cdot\text{s/m}, k = 25.0\,\text{N/m}\)):
\[\mathbf{A} = \begin{bmatrix} 0 & 1 \\ -25 & -6 \end{bmatrix}, \quad \mathbf{B} = \begin{bmatrix} 0 \\ 1 \end{bmatrix}\]
\[\mathbf{C} = \begin{bmatrix} 1 & 0 \end{bmatrix}, \quad D = 0\]
✓ Pembuktian Ekuivalensi: Nilai Eigen Matriks A = Kutub PDB

Persamaan karakteristik matriks: \(\det(s\mathbf{I} - \mathbf{A}) = s(s + 6) - (-1)(25) = s^2 + 6s + 25 = 0\).

\[\lambda_{1,2} = -3 \pm j4\]

Nilai eigen matriks ruang keadaan terbukti identik mutlak dengan kutub fungsi alih Laplace, menjamin integritas fisik pemodelan.

Simulator Dinamis & Parameter Tuning Live

Ubah parameter fisik aktuator dan amati respon transien grafik posisi \(y(t)\) dan kecepatan \(v(t)\) seketika.

tune Slider Parameter Fisik

Massa \(m\): 1.0 kg
Redaman Viskos \(c\): 6.0 N·s/m
Kekakuan Kontak \(k\): 25.0 N/m
Gaya Dorong Masukan \(A\): 1.0 N
Karakteristik Kanonik Terhitung:
\(\omega_n\): 5.00 rad/s
\(\zeta\): 0.600
\(\omega_d\): 4.00 rad/s
Rezim: Underdamped
Kutub Bidang-\(s\): -3.00 ± j4.00

show_chart Respons Waktu-Nyata (Chart.js Engine)

Posisi \(y(t) = x_1\) dan Kecepatan \(v(t) = x_2\)

Posisi Steady-State \(y_{ss} = A/k\): 0.040 m Settling Time (\(2\%\)): \(\approx 1.33\text{ s}\)

Galeri Grafik Validasi Simulasi & Analisis Frekuensi

Grafik resmi hasil eksekusi Python dan MATLAB dari direktori figures/.

Karakteristik Gesekan Gbr. 2 Paper MSSP
Karakteristik Gesekan vs Kecepatan

Membuktikan lenyapnya histeresis gesek padat pada ASCPA berkat pelumasan aerostatik.

S Plane Pole Zero Gbr. 3 Konstelasi Kutub
Konstelasi Kutub Bidang-\(s\)

Sepasang kutub kompleks \(p_{1,2} = -3 \pm j4\) pada LHP dengan sudut \(\theta = 126.87^\circ\).

Diagram Bode Gbr. 4 Diagram Bode
Respon Frekuensi & Diagram Bode

Filter lolos-rendah orde 2 dengan roll-off curam \(-40\,\text{dB/dekade}\) dan fasa hingga \(-180^\circ\).

Respons Pulsa Gbr. 5 Pulsa Persegi
Respons Pulsa Persegi Tunggal

Respon domain waktu untuk durasi \(T = 0.5, 1.0, 2.0\,\text{s}\) (eksitasi paksa & relaksasi bebas).

Spektrum CTFT Gbr. 6 CTFT Spektrum
Spektrum Fourier CTFT

Spektrum keluaran meluruh pada \(\mathcal{O}(1/\omega^3)\) (\(-60\,\text{dB/dekade}\)), memperhalus lonjakan katup.

Impuls dan Undak Gbr. 9 Impuls & Step
Respons Impuls & Undak Satuan

Validasi analitis vs integrasi numerik: overshoot \(9.48\%\), waktu puncak \(t_p = 0.785\text{ s}\).