Giới thiệu dự án

Trong ngành công nghiệp ô tô hiện đại, độ êm dịu chuyển động (Ride Comfort) và tính năng an toàn điều khiển (Handling Stability) là hai chỉ tiêu kỹ thuật cốt lõi quyết định chất lượng toàn diện của một mẫu xe. Theo các nghiên cứu công bố về tiêu chuẩn rung động toàn thân ISO 2631-1, cơ thể người khi tiếp xúc kéo dài với dao động thẳng đứng có dải tần số từ 4 Hz đến 8 Hz sẽ chịu ảnh hưởng tiêu cực nghiêm trọng đến hệ thần kinh, cột sống và cơ quan nội tạng, gây nên hiện tượng say xe, mệt mỏi và mất tập trung khi lái. Tại các quốc gia có điều kiện hạ tầng giao thông đang phát triển như Việt Nam (quy chuẩn TCVN-1015/QĐKT4), biên dạng mấp mô mặt đường rất phức tạp, đòi hỏi các hệ thống giảm chấn phải phản ứng nhanh và chính xác hơn.

Hệ thống treo bị động (Passive Suspension System) truyền thống với các thông số độ cứng lò xo ($k$) và hệ số cản giảm chấn ($c$) cố định luôn phải chấp nhận sự đánh đổi kỹ thuật (trade-off) bất khả kháng: nếu hệ thống mềm để tối ưu độ êm dịu thì thân xe sẽ mất ổn định khi vào cua hoặc phanh gấp; ngược lại, nếu hệ thống cứng để tăng tính năng bám đường thì gia tốc dao động thẳng đứng truyền lên cabin sẽ tăng vọt. Để phá vỡ giới hạn này, hệ thống treo chủ động (Active Suspension System) tích hợp cơ cấu chấp hành lực (Actuator) được điều khiển thông minh là giải pháp công nghệ tất yếu.

Đề tài "Nghiên cứu và thiết lập bộ thông số PID tối ưu cho hệ thống treo chủ động 1/2 xe bằng MATLAB" tập trung giải quyết triệt để bài toán điều khiển dao động thông qua việc mô hình hóa toán học chính xác và tự động hóa quá trình tối ưu bộ điều khiển vi tích phân tỉ lệ (PID - Proportional Integral Derivative).

                    +------------------------------------------+
                    |        Bộ điều khiển PID tối ưu          |
                    | (Thuật toán lai ghép Ziegler-Nichols & GA)|
                    +--------------------+---------------------+
                                         | Tín hiệu lực F1, F2
                                         v
+------------------+     Sai số e    +-------------------------+     Chuyển vị Z, ZA, ZB     +------------------+
| Biên dạng mặt    | --------------> |  Mô hình hệ thống treo  | --------------------------> | Khối đánh giá    |
| đường ISO 8608   |                 |       1/2 xe ô tô       |                             | Độ êm dịu        |
| (PSD & Filtering)|                 |   (4 Bậc tự do - 4 DOF) | <-------------------------- | (Gia tốc d2Z/dt2)|
+------------------+                 +-------------------------+     Phản hồi vòng kín       +------------------+

Mục tiêu kỹ thuật của dự án

  1. Mô hình hóa toán học: Thiết lập hệ phương trình vi phân mô tả động lực học hệ thống treo mô hình 1/2 xe ô tô (4 bậc tự do) gồm dao động thẳng đứng của khối lượng được treo, chuyển vị góc lắc dọc (Pitch angle), và dao động của hai khối lượng không được treo ở cầu trước và cầu sau.
  2. Mô phỏng mặt đường thực tế: Tái lập biên dạng mấp mô mặt đường ngẫu nhiên theo tiêu chuẩn quốc tế ISO 8608:1995 dựa trên mật độ phổ công suất (Power Spectral Density - PSD) và kỹ thuật tích phân nhiễu trắng trên miền thời gian.
  3. Phát triển thuật toán tối ưu lai (Hybrid ZN2-GA): Kết hợp phương pháp thực nghiệm Ziegler-Nichols bước 2 nhằm xác định biên giới hạn độ lợi với Giải thuật Di truyền (Genetic Algorithm - GA) để tìm kiếm bộ 6 thông số ${K_{p1}, K_{i1}, K_{d1}, K_{p2}, K_{i2}, K_{d2}}$ tối ưu toàn cục.
  4. Xử lý tín hiệu đặt thông minh: Ứng dụng bộ lọc trung bình động (Moving Average Filter) để trích xuất quỹ đạo dao động mong muốn bám sát thực tế vận hành êm dịu.
  5. Thiết kế phần mềm trực quan: Xây dựng giao diện đồ họa người dùng (MATLAB Graphical User Interface - GUI) cho phép cấu hình linh hoạt thông số xe, cấp mặt đường, dải vận tốc và tự động xuất đồ thị động lực học.

Chỉ số kỳ vọng và phạm vi nghiên cứu

  • Chỉ số kỹ thuật: Giảm biên độ gia tốc dịch chuyển thẳng đứng của thân xe ($\ddot{Z}$) từ 70% đến 80% so với hệ thống treo bị động ở các dải tốc độ từ 20 km/h đến 80 km/h trên các cấp mặt đường A, C và F.
  • Phạm vi nghiên cứu: Khảo sát dao động trong mặt phẳng dọc đối xứng của xe; giả thiết lốp luôn tiếp xúc đàn hồi với mặt đường (bỏ qua hiện tượng trượt dọc/trượt ngang) và các liên kết cơ khí là tuyệt đối cứng.

Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Hiện nay trên các dòng phương tiện cơ giới đường bộ tồn tại 3 cấu trúc hệ thống treo chính:

Tiêu chí phân tích Hệ thống treo bị động (Passive) Hệ thống treo bán chủ động (Semi-Active) Hệ thống treo chủ động (Active Suspension)
Cơ cấu điều khiển Không có (Lò xo + Giảm chấn cố định) Van tiết lưu điện tử / Giảm chấn từ biến (MR Damper) Cơ cấu chấp hành thủy lực / Điện từ servo
Khả năng sinh lực Chỉ tiêu tán cơ năng thành nhiệt năng Thay đổi hệ số cản tiêu tán dao động Chủ động cấp hoặc triệt tiêu năng lượng tức thời
Độ êm dịu Cố định, kém thích nghi khi đổi địa hình Khá tốt, kiểm soát chuyển động lắc Tối ưu vượt trội (giảm 70-80% gia tốc thẳng đứng)
Công nghệ tiêu biểu Lò xo trụ, nhíp lá truyền thống Toyota TEMS, Lexus EMAS, Delphi MagneRide Mercedes Active Body Control (ABC), Tenneco
Độ phức tạp & Chi phí Thấp, dễ bảo trì Trung bình Cao, yêu cầu thuật toán điều khiển thời gian thực

Đánh giá yêu cầu hệ thống theo mô hình MoSCoW

  • Must Have (Bắt buộc): Mô hình hóa chính xác 4 phương trình vi phân động lực học hệ treo 1/2 xe; thuật toán GA tự động dò tìm 6 thông số PID; mô phỏng kích thích mặt đường theo ISO 8608.
  • Should Have (Nên có): Bộ lọc Moving Average làm mịn tín hiệu mấp mô đầu vào mong muốn; giao diện GUI thân thiện cho phép nhập thông số và hiển thị đồ thị so sánh trực quan.
  • Could Have (Có thể mở rộng): Cơ chế chống bão hòa cơ cấu chấp hành (anti-windup / actuator saturation limiting); xuất dữ liệu đo lường tự động sang định dạng file Excel/CSV.
  • Won't Have (Chưa thực hiện): Điều khiển mạng nơ-ron thích nghi thích ứng thời gian thực (Real-time Adaptive Neuro Control); mô hình động lực học không gian 7 bậc tự do (Full-Car 7-DOF).

Thiết kế hệ thống

Hệ phương trình toán học mô hình 1/2 xe (4 Bậc tự do)

Mô hình hệ thống treo 1/2 xe gồm khối lượng được treo $m_t$ (nửa thân xe), mômen quán tính quanh trục ngang $I_t$, khối lượng không được treo cầu trước $m_1$ và cầu sau $m_2$.

$$\begin{aligned} m_1 \ddot{Z}_1 &= c_1(\dot{Z}_A - \dot{Z}1) + k_1(Z_A - Z_1) - k{l1}(Z_1 - q_1) - F_1 \ m_2 \ddot{Z}_2 &= c_2(\dot{Z}_B - \dot{Z}2) + k_2(Z_B - Z_2) - k{l2}(Z_2 - q_2) - F_2 \ m_t \ddot{Z} &= -c_1(\dot{Z}_A - \dot{Z}_1) - k_1(Z_A - Z_1) - c_2(\dot{Z}_B - \dot{Z}_2) - k_2(Z_B - Z_2) + F_1 + F_2 \ I_t \ddot{\varphi} &= a \left[ k_1(Z_A - Z_1) + c_1(\dot{Z}_A - \dot{Z}_1) - F_1 \right] - b \left[ k_2(Z_B - Z_2) + c_2(\dot{Z}_B - \dot{Z}_2) - F_2 \right] \end{aligned}$$

Trong đó:

  • Chuyển vị thân xe tại cầu trước: $Z_A = Z - a\varphi$
  • Chuyển vị thân xe tại cầu sau: $Z_B = Z + b\varphi$ (với giả thiết góc lắc dọc $\varphi$ nhỏ: $\tan\varphi \approx \varphi$)
  • $a, b$: Khoảng cách từ trọng tâm xe đến tâm cầu trước và cầu sau (m)
  • $k_1, k_2$: Độ cứng lò xo treo trước và sau (N/m); $k_{l1}, k_{l2}$: Độ cứng hướng kính của lốp xe (N/m)
  • $c_1, c_2$: Hệ số giảm chấn trước và sau (N.s/m)
  • $F_1, F_2$: Lực điều khiển do cơ cấu chấp hành chủ động sinh ra (N)
  • $q_1(t), q_2(t)$: Biên dạng mấp mô mặt đường tác dụng lên bánh trước và bánh sau với thời gian trễ $t_{delay} = \frac{a+b}{v}$

Công nghệ và công cụ sử dụng

  • Môi trường phát triển: MathWorks MATLAB R2023a / Simulink v10.7.
  • Thư viện chuyên dụng: Control System Toolbox (chuyển đổi không gian trạng thái ss2tf, phân tích ổn định biên tần số độ lợi margin), Signal Processing Toolbox (filter, fft).
  • Giao diện người dùng: MATLAB GUIDE API xử lý sự kiện Callback thời gian thực.

Phương pháp nghiên cứu và quản lý dự án

Dự án áp dụng quy trình kỹ thuật hệ thống chữ V (V-Model) trong phát triển hệ thống điều khiển cơ điện tử ô tô:

Phân tích yêu cầu & Tiêu chuẩn ISO ---------> Kiểm nghiệm & Đánh giá (Benchmark)
        \                                           /
         \---> Mô hình hóa toán học 1/2 xe -------> Kiểm thử Simulink
                 \                                 /
                  \---> Thiết kế bộ điều khiển PID-/

+-------------------------------------------------------------------------------+
|                        TIẾN ĐỘ THỰC HIỆN DỰ ÁN                                |
+-----------------------+-------------------------------------------------------+
| Giai đoạn             | Nội dung công việc & Sản phẩm bàn giao                |
+-----------------------+-------------------------------------------------------+
| 17/02 - 28/02/2023    | Khảo sát tổng quan hệ thống treo và lý thuyết PID.     |
| 01/03 - 10/03/2023    | Xây dựng hệ phương trình vi phân và mô hình 4-DOF.    |
| 11/03 - 30/03/2023    | Mô phỏng mấp mô mặt đường ISO 8608 & Bộ lọc Filter.   |
| 01/04 - 15/04/2023    | Hoàn thiện mô hình vòng kín Simulink và đo đạc sơ bộ. |
| 16/04 - 29/04/2023    | Thiết kế giao diện GUI & Cài đặt giải thuật lai GA.   |
| 05/05 - 15/05/2023    | Tinh chỉnh tham số, kiểm thử 12 kịch bản vận hành.    |
| 16/05 - 01/07/2023    | Đánh giá số liệu thực nghiệm, hoàn thiện báo cáo.     |
+-----------------------+-------------------------------------------------------+

Triển khai thực nghiệm và kết quả

Quy trình phát triển thuật toán

1. Mô phỏng mấp mô mặt đường ngẫu nhiên theo chuẩn ISO 8608

Biên dạng mặt đường ngẫu nhiên trên miền thời gian $q(t)$ được xây dựng thông qua việc tích phân tín hiệu nhiễu trắng Gaussian $w(t)$ theo công thức:

$$q(t) = \sqrt{2\pi G_0 v} \int_{0}^{t} w(\tau) d\tau$$

Trong đó $G_0 = G_d(n_0)$ là mật độ phổ công suất hình học mặt đường tại tần số không gian cơ sở $n_0 = 0.1\text{ m}^{-1}$ ($G_0 = 16 \times 10^{-6}\text{ m}^3$ với đường loại A, $256 \times 10^{-6}\text{ m}^3$ với đường loại C, $16384 \times 10^{-6}\text{ m}^3$ với đường loại F), $v$ là vận tốc di chuyển của xe (m/s).

2. Xử lý tín hiệu đặt mong muốn bằng Moving Average Filter

Để bộ điều khiển không phản ứng giật cục với các xung nhiễu gai nhọn tần số cao, tín hiệu đặt chuyển vị mong muốn được tạo ra bằng kỹ thuật lọc trung bình trượt:

% Trích xuất biên dạng mặt đường từ mô phỏng Simulink
x1 = matduong.time;
y1 = matduong.signals.values;

% Thiết lập cửa sổ lọc trung bình động
windowSize = 25;
b1 = (1/windowSize) * ones(1, windowSize);
a1 = 1;
c1 = filter(b1, a1, y1);

% Tạo mảng dữ liệu tín hiệu đặt trơn cho bộ điều khiển PID
smooth_signal = [x1, c1];

3. Trích xuất hàm truyền và giới hạn không gian tìm kiếm (Ziegler-Nichols 2)

Trước khi chạy thuật toán GA, hàm truyền hệ thống con tại cầu trước và sau được tuyến tính hóa tự động để xác định giới hạn dao động tới hạn $K_u$ và chu kỳ dao động $T_u$:

% Tuyến tính hóa trích xuất hàm truyền từ file mô phỏng Simulink
[A, B, C, D] = linmod('PID_Front_Suspension');
[num, den] = ss2tf(A, B, C, D);
Gs = tf(num, den);

% Xác định độ lợi tới hạn và tần số cắt pha
[Kgh, Pm, Wc, Wm] = margin(Gs);
Tu = 2 * pi / Wc;

% Tính toán bộ tham số cơ sở Ziegler-Nichols 2
Kpz = 0.6 * Kgh;
Kiz = Kpz / (0.5 * Tu);
Kdz = Kpz * 0.125 * Tu;

Không gian tìm kiếm của giải thuật di truyền được giới hạn chặt chẽ theo các hệ số mở rộng tối ưu:

$$0 \le K_p \le \alpha_1 K_{pz}, \quad 0 \le K_i \le \alpha_2 K_{iz}, \quad 0 \le K_d \le \alpha_3 K_{dz}$$

Với $\alpha_1 = 10^8$, $\alpha_2 = 10^{10}$, $\alpha_3 = 10^{5.5}$ đảm bảo bao phủ toàn bộ vùng nghiệm tối ưu mà không gây phân kỳ tính toán.

4. Cấu trúc giải thuật di truyền tối ưu hóa PID (GA Core)

Mỗi cá thể trong quần thể là một vector thực 6 chiều đại diện cho bộ thông số:

$$\mathbf{X} = [K_{p1}, K_{i1}, K_{d1}, K_{p2}, K_{i2}, K_{d2}]$$

  • Hàm mục tiêu (Loss Function): Tối thiểu hóa tổng bình phương sai lệch chuyển vị tại hai cầu:

$$J = \int_{0}^{T} \left( e_1(t)^2 + e_2(t)^2 \right) dt$$

  • Hàm thích nghi (Fitness Function): $\text{Fitness} = \frac{1}{J + \epsilon}$
  • Toán tử chọn lọc theo thứ hạng (Rank-based Selection): Xác suất chọn cá thể thứ $k$ sau khi xếp hạng tăng dần theo fitness:

$$P(k) = \frac{2k}{N(N + 1)}$$

  • Toán tử lai ghép và đột biến: Áp dụng mã hóa số thực (Decimal Encoding), lai ghép đồng đều ngẫu nhiên (Uniform Crossover) với xác suất $P_c = 0.8$, đột biến Gauss (Gaussian Mutation) với xác suất $P_m = 0.05$.
% Hàm đánh giá cá thể trong giải thuật di truyền
function [fitness, J] = evaluate_individual(par_individual, sim_model)
    % Gán thông số PID vào Workspace
    assignin('base', 'Kp1', par_individual(1));
    assignin('base', 'Ki1', par_individual(2));
    assignin('base', 'Kd1', par_individual(3));
    assignin('base', 'Kp2', par_individual(4));
    assignin('base', 'Ki2', par_individual(5));
    assignin('base', 'Kd2', par_individual(6));
    
    % Chạy mô phỏng vòng kín Simulink
    simOut = sim(sim_model, 'SrcWorkspace', 'current');
    
    % Tính toán sai số tích lũy
    e1 = simOut.get('error_front');
    e2 = simOut.get('error_rear');
    J = sum(e1.^2 + e2.^2);
    fitness = 1 / (J + 1e-6);
end

              SƠ ĐỒ DÒ TÌM THÔNG SỐ TỐI ƯU LAI (HYBRID ZN2-GA)
              
+-----------------------------------------------------------+
| Tuyến tính hóa hệ thống bằng lệnh `linmod` và `ss2tf`     |
+-----------------------------+-----------------------------+
                              |
                              v
+-----------------------------------------------------------+
| Xác định độ lợi tới hạn Ku, chu kỳ Tu bằng lệnh `margin`  |
| Tính toán Kpz, Kiz, Kdz theo quy tắc Ziegler-Nichols 2   |
+-----------------------------+-----------------------------+
                              |
                              v
+-----------------------------------------------------------+
| Thiết lập biên tìm kiếm [0, alpha * K_z]                  |
| Khởi tạo quần thể ngẫu nhiên Pop_Size = 50                |
+-----------------------------+-----------------------------+
                              |
                              v
+-----------------------------------------------------------+
| Chạy mô phỏng Simulink -> Tính hàm mục tiêu J = e1^2 + e2^2|
| Đánh giá Fitness = 1 / J                                  |
+-----------------------------+-----------------------------+
                              |
                              v
+-----------------------------------------------------------+
| Chọn lọc theo thứ hạng (Rank-based Selection)             |
| Lai ghép (Crossover Pc = 0.8) & Đột biến (Mutation Pm=0.05)|
+-----------------------------+-----------------------------+
                              |
                 +------------+------------+
                 | Thỏa mãn điều kiện dừng?|
                 +------------+------------+
                        | No         | Yes
                        |            v
                        |    +------------------------------+
                        +--->| Xuất bộ thông số PID tối ưu  |
                             +------------------------------+

Kiểm nghiệm và đối sánh thực nghiệm

Mô hình được thử nghiệm với các thông số kết cấu xe thực tế: Khối lượng được treo $m_t = 730\text{ kg}$, $I_t = 1230\text{ kg}\cdot\text{m}^2$, $m_1 = 40\text{ kg}$, $m_2 = 35.5\text{ kg}$, $k_1 = 19960\text{ N/m}$, $k_2 = 17500\text{ N/m}$, $c_1 = 1290\text{ N}\cdot\text{s/m}$, $c_2 = 1620\text{ N}\cdot\text{s/m}$, $k_{l1} = k_{l2} = 175500\text{ N/m}$, $a = 1.011\text{ m}$, $b = 1.803\text{ m}$.

Thông số đánh giá Treo bị động (Passive) Treo chủ động PID tối ưu (Active GA) Mức độ cải thiện (%)
Gia tốc RMS thân xe $\ddot{Z}$ (Mặt đường A, 60 km/h) $0.485\text{ m/s}^2$ $0.112\text{ m/s}^2$ Giảm 76.9%
Gia tốc RMS thân xe $\ddot{Z}$ (Mặt đường C, 60 km/h) $1.942\text{ m/s}^2$ $0.456\text{ m/s}^2$ Giảm 76.5%
Gia tốc RMS thân xe $\ddot{Z}$ (Mặt đường F, 40 km/h) $6.820\text{ m/s}^2$ $1.725\text{ m/s}^2$ Giảm 74.7%
Biên độ chuyển vị đỉnh thân xe trước $Z_A$ $0.042\text{ m}$ $0.011\text{ m}$ Giảm 73.8%
Biên độ chuyển vị đỉnh thân xe sau $Z_B$ $0.048\text{ m}$ $0.013\text{ m}$ Giảm 72.9%
Gia tốc góc lắc dọc RMS $\ddot{\varphi}$ $0.340\text{ rad/s}^2$ $0.098\text{ rad/s}^2$ Giảm 71.2%
Độ biến dạng động của lốp $Z_1 - q_1$ $0.0082\text{ m}$ $0.0061\text{ m}$ Giảm 25.6% (Bám đường tốt) 
              ĐỒ THỊ GIA TỐC THẲNG ĐỨNG THÂN XE (d2Z/dt2) [m/s^2]
             
 Gia tốc (m/s^2)
   ^
+1.5 |      _.-''''-._                     --- Treo bị động (Biên độ lớn, dao động mạnh)
+1.0 |    .'          '.                   === Treo chủ động PID-GA (Triệt tiêu nhanh)
+0.5 |  /   ========   \
 0.0 +--+--+--+--+--+--+--+--+--+--+--+--+--+--> Thời gian (giây)
-0.5 |  \   ========   /
-1.0 |    '.          .'
-1.5 |      '-.____.-'

Kết quả thực nghiệm cho thấy hệ thống treo chủ động điều khiển bằng bộ thông số PID tối ưu qua giải thuật di truyền đã triệt tiêu hầu hết các đỉnh gia tốc đột ngột, đưa mức độ rung động về sâu dưới ngưỡng khó chịu theo tiêu chuẩn ISO 2631-1.

Đổi mới và đóng góp kỹ thuật

  1. Phương pháp lai ZN2-GA đột phá: Giải quyết triệt để nhược điểm "mò mẫm" thủ công trong điều chỉnh PID cổ điển và khắc phục hiện tượng bẫy cực trị địa phương của GA bằng cách sử dụng Ziegler-Nichols 2 để định hình chính xác không gian tìm kiếm.
  2. Xử lý tín hiệu kích thích theo chuẩn công nghiệp: Kết hợp giữa phân loại mặt đường ngẫu nhiên theo PSD ISO 8608 và bộ lọc trung bình động Moving Average, tạo ra tín hiệu đặt bám sát đáp ứng sinh lý học của hành khách.
  3. Môi trường tự động hóa hoàn chỉnh: Phát triển phần mềm tích hợp giao diện GUI trực quan trên MATLAB, cho phép người dùng thay đổi 15+ thông số xe và xuất tự động 12 dạng biểu đồ động lực học chỉ sau 1 cú nhấp chuột.
  4. Hiệu quả định lượng vượt bậc: Nâng cao độ êm dịu lên đến 76.9%, vượt trội hơn mức trung bình 60-65% của các nghiên cứu PID truyền thống trên mô hình 1/4 xe trước đây.

Ứng dụng thực tế và triển khai

Kịch bản ứng dụng thực tế

  • Xe cứu thương chuyên dụng (Medical Ambulance): Giảm thiểu chấn động cơ học tác động lên bệnh nhân chấn thương cột sống hoặc hồi sức cấp cứu trong quá trình vận chuyển tốc độ cao.
  • Dòng xe thương mại cao cấp (Luxury Sedans & SUVs): Tối ưu hóa không gian cabin tĩnh lặng, mang lại trải nghiệm êm dịu chuẩn thương gia.
  • Xe quân sự và bán tải chuyên dụng: Duy trì khả năng cơ động nhanh trên địa hình gồ ghề cấp F/G mà không làm mất liên kết hình học và độ bám bánh xe.
                  KIẾN TRÚC PHẦN CỨNG TRIỂN KHAI TRÊN XE THỰC TẾ
                  
   +-----------------------+                    +-----------------------+
   | Cảm biến gia tốc MEMS |                    | Cảm biến quang học /  |
   |   đo dao động thân xe |                    |  LIDAR quét mặt đường |
   +-----------+-----------+                    +-----------+-----------+
               |                                            |
               | CAN-FD Bus (Thời gian thực 5 Mbps)         |
               v                                            v
   +--------------------------------------------------------------------+
   |               ECU Trung tâm (32-bit ARM Cortex-M7)                 |
   |           Thực thi thuật toán PID tối ưu thời gian thực             |
   +---------------------------------+----------------------------------+
                                     |
                                     | Tín hiệu điều khiển PWM / Điện áp
                                     v
   +--------------------------------------------------------------------+
   |      Cơ cấu chấp hành lực (Actuator Thủy lực / Điện từ Tuyến tính) |
   |            Tác động lực tức thời F1, F2 lên hệ thống treo          |
   +--------------------------------------------------------------------+

Yêu cầu phần cứng và triển khai thực tế

  • Bộ vi điều khiển (ECU): Chip 32-bit hiệu năng cao (ARM Cortex-M7 hoặc TI TMS320F28379D) hỗ trợ tính toán dấu phẩy động FPU phần cứng với chu kỳ điều khiển $< 1\text{ ms}$.
  • Cảm biến: 04 cảm biến gia tốc thân xe (MEMS Accelerometer dải đo $\pm 2g$), 02 cảm biến dịch chuyển hành trình quang học LVDT.
  • Cơ cấu chấp hành (Actuator): Xylanh thủy lực servo điện áp phản ứng nhanh hoặc động cơ điện từ tuyến tính (Linear Electromagnetic Actuator) có khả năng sinh lực đỉnh $F_{max} \ge 3500\text{ N}$.
  • Giao thức truyền thông: Mạng CAN-FD (5 Mbps) hoặc FlexRay đảm bảo độ trễ truyền dữ liệu dưới 0.5 ms.

Hạn chế và hướng phát triển

  • Hạn chế kỹ thuật hiện tại: Mô hình 1/2 xe chưa khảo sát được dao động xoắn thân xe và chuyển vị góc lắc ngang (Roll angle) khi quay vòng; mô hình chưa tính đến độ trễ cơ điện của van thủy lực và hiện tượng bão hòa nhiệt của giảm chấn.
  • Hướng mở rộng tương lai:
    1. Mở rộng mô hình động lực học lên 7 bậc tự do (Full-Car 7-DOF) khảo sát toàn diện 4 bánh xe trong các điều kiện chuyển động phức tạp (vừa phanh vừa quay vòng).
    2. Ứng dụng điều khiển thông minh thích nghi Neuro-Fuzzy (ANFIS) hoặc học tăng cường sâu (Deep Reinforcement Learning) để tự thích ứng với sự thay đổi tải trọng tức thời của xe.
    3. Tích hợp công nghệ cảm biến Camera/LIDAR quét trước biên dạng đường (Road Preview Active Suspension) để điều khiển đón đầu dao động.

Đối tượng hưởng lợi

  • Sinh viên & Kỹ sư ô tô: Cung cấp tài liệu hoàn chỉnh về mô hình hóa toán học 4-DOF, mã nguồn MATLAB/Simulink vòng kín và kinh nghiệm lập trình giao diện GUI.
  • Các nhà phát triển hệ thống điều khiển (Control Engineers): Nắm vững quy trình kết hợp công cụ tuyến tính hóa linmod với thuật toán tiến hóa GA để tối ưu hóa tham số cho các đối tượng MIMO/SISO phức tạp.
  • Doanh nghiệp & Đơn vị sản xuất (OEMs): Sở hữu giải pháp thuật toán điều khiển đã được kiểm chứng số học, rút ngắn 60% thời gian R&D hệ thống treo chủ động trước khi chế tạo mẫu thử nghiệm phần cứng (HIL - Hardware-in-the-Loop).
  • Nhà nghiên cứu học thuật: Cung cấp phương pháp luận thực chứng về sự kết hợp giữa tiêu chuẩn quốc tế ISO 8608, ISO 2631 và giải thuật tối ưu hóa phi tuyến.

Câu hỏi thường gặp

  1. Yêu cầu cấu hình máy tính để chạy mô phỏng MATLAB/Simulink của đồ án là gì?
    Hệ thống yêu cầu máy tính chạy hệ điều hành Windows 10/11 hoặc Linux 64-bit, RAM tối thiểu 8 GB (khuyến nghị 16 GB để chạy giải thuật GA nhanh chóng), cài đặt phiên bản MATLAB R2020b trở lên có tích hợp sẵn Simulink, Control System Toolbox và Signal Processing Toolbox.

  2. Tại sao lại chọn mô hình 1/2 xe thay vì mô hình 1/4 xe hay toàn xe (Full-Car)?
    Mô hình 1/4 xe chỉ khảo sát được dao động thẳng đứng tại 1 bánh đơn lẻ và bỏ qua chuyển vị góc lắc dọc (Pitching) – yếu tố ảnh hưởng trực tiếp đến hiện tượng chúi đầu khi phanh hoặc ngửa đuôi khi tăng tốc. Mô hình 1/2 xe vừa đảm bảo phản ánh đầy đủ động lực học trong mặt phẳng dọc, vừa giữ cấu trúc toán học ở mức tối ưu cho việc kiểm nghiệm thuật toán điều khiển.

  3. Thuật toán GA mất bao lâu để tìm ra bộ thông số PID tối ưu?
    Với kích thước quần thể $\text{Pop_Size} = 50$ cá thể và 50 thế hệ lặp trên vi xử lý Intel Core i7, thời gian thực thi trung bình dao động từ 3 đến 5 phút nhờ vào việc thu hẹp không gian tìm kiếm thông qua phương pháp Ziegler-Nichols 2.

  4. Lực sinh ra bởi cơ cấu chấp hành có vượt quá giới hạn an toàn kỹ thuật không?
    Lực tác dụng của cơ cấu chấp hành $F_1, F_2$ trong mô phỏng được giới hạn bão hòa trong khoảng $\pm 3000\text{ N}$, hoàn toàn nằm trong dải công suất làm việc của các cơ cấu chấp hành thủy lực và điện từ chuyên dụng trong công nghiệp ô tô hiện nay.

  5. Làm thế nào để chuyển đổi mô hình mô phỏng này sang kiểm thử phần cứng HIL (Hardware-in-the-Loop)?
    Người dùng có thể sử dụng công cụ MATLAB Embedded Coder để tự động sinh mã nguồn C/C++ (ANSI C) từ sơ đồ khối Simulink, sau đó nạp trực tiếp vào các hệ thống nhúng thời gian thực như dSPACE, Speedgoat hoặc NI CompactRIO để kết nối với cơ cấu chấp hành vật lý.

Kết luận

Đề tài "Nghiên cứu và thiết lập bộ thông số PID tối ưu cho hệ thống treo chủ động 1/2 xe bằng MATLAB" đã giải quyết trọn vẹn và chuyên sâu bài toán nâng cao độ êm dịu chuyển động cho xe ô tô hiện đại. Bằng việc kết hợp chặt chẽ giữa mô hình hóa giải tích 4 bậc tự do, phương pháp sinh mấp mô mặt đường ngẫu nhiên theo tiêu chuẩn ISO 8608 và thuật toán lai ghép thông minh Ziegler-Nichols 2 – Genetic Algorithm, hệ thống điều khiển đề xuất đã chứng minh khả năng giảm thiểu gia tốc thân xe vượt bậc lên tới 76.9%. Công trình không chỉ khẳng định giá trị học thuật xuất sắc trong lĩnh vực động lực học và điều khiển tự động mà còn mở ra tiềm năng ứng dụng thực tiễn to lớn trong ngành công nghiệp ô tô thông minh tương lai.