Giới thiệu dự án

Trong tiến trình công nghiệp hóa và hiện đại hóa ngành cơ khí chế tạo, bài toán nội địa hóa và nâng cao chất lượng phương tiện vận tải hành khách công cộng luôn giữ vai trò chiến lược. Quyết định cấm xe cơ giới ba bánh (xe lam) lưu hành tại các đô thị lớn của Bộ Giao thông Vận tải đã đặt ra nhu cầu cấp thiết về việc phát triển các dòng xe khách cỡ nhỏ (minibus 6–8 chỗ ngồi) mang thương hiệu Việt Nam. Đề tài nghiên cứu được thực hiện trong khuôn khổ Dự án cấp Nhà nước KC.13: "Hoàn thiện thiết kế, công nghệ chế tạo và lắp ráp dòng xe minibus thông dụng 6–8 chỗ ngồi mang nhãn hiệu Việt Nam", do Trung tâm Phát triển Công nghệ Ô tô thuộc Hội Kỹ sư Ô tô Việt Nam (VSAE) chủ trì dưới sự định hướng chuyên môn của PGS. TS. Dư Quốc Thịnh.

Hệ thống treo là kết cấu then chốt liên kết đàn hồi giữa cầu xe (hoặc bánh xe) với khung/thân vỏ, quyết định trực tiếp đến độ êm dịu chuyển động (ride comfort), tính năng ổn định quỹ đạo (handling/stability) và độ an toàn chuyển động (safety dynamics). Trên dòng xe MEFA5-LAVI-304N (sử dụng động cơ phun xăng điện tử đa điểm, khả năng vượt dốc $\ge 32%$, mức tiêu hao nhiên liệu $6.0 - 7.1\text{ lít}/100\text{km}$), việc tích hợp hệ thống treo trước độc lập kiểu MacPherson/hai đòn ngang và hệ thống treo sau phụ thuộc đòi hỏi một quy trình khảo sát động lực học chính xác.

Mục tiêu dự án

  1. Xây dựng hệ thống cơ sở lý thuyết: Khảo sát toàn diện kết cấu, nguyên lý vận hành của các loại hệ thống treo phổ biến (treo phụ thuộc nhíp lá/lò xo, treo độc lập hai đòn ngang, MacPherson, đa liên kết Multi-link, hệ thống treo khí nén điện tử EAS).
  2. Thiết lập mô hình động lực học không gian 3D: Phân tích dao động riêng biệt và kết hợp trên hai mặt phẳng:
    • Mặt phẳng dọc ($x-z$): Khảo sát dao động tịnh tiến đứng $Z_s$ và góc lắc dọc $\varphi_y$ (pitching).
    • Mặt phẳng ngang ($y-z$): Khảo sát dao động tịnh tiến đứng $Z_s$ và góc lắc ngang $\varphi_x$ (rolling) cho cả hai cấu hình treo độc lập và phụ thuộc.
  3. Chuyển đổi sang hệ phương trình không gian trạng thái: Thiết lập các ma trận khối lượng $[M]$, ma trận cản nhớt $[C]$, ma trận độ cứng $[K]$ và đưa về dạng vi phân bậc nhất tuyến tính $\dot{x} = Ax + Bu, y = Cx + Du$.
  4. Mô phỏng số trên môi trường MATLAB/Simulink: Đánh giá định lượng ảnh hưởng của độ cứng lò xo $K_s$ và hệ số cản giảm chấn $C_s$ dưới kích động mặt đường hình sin (biên độ $A = 0.025\text{ m}$, tần số biến thiên theo vận tốc thử nghiệm $V = 10\text{ m/s}$).
  5. Đánh giá tiêu chuẩn độ êm dịu: So sánh các chỉ tiêu sai phương gia tốc ($\text{Var}(\ddot{Z}_s)$, $\text{Var}(\ddot{\varphi})$) với giới hạn cảm nhận con người theo tiêu chuẩn quốc tế ISO/DIS 2631.

Phạm vi và giới hạn

  • Khối lượng phần được treo ($m_s = 1555\text{ kg}$) và không được treo ($m_u = 50\text{ kg}$) được giả thiết là các vật rắn tuyệt đối.
  • Bỏ qua ma sát phi tuyến tại các khớp cầu và bạc cao su liên kết; kích động mặt đường được coi là nguồn kích động ngẫu nhiên/chu kỳ duy nhất.

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

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

Tiêu chí kỹ thuật Treo phụ thuộc (Nhíp lá/Cầu liền) Treo độc lập MacPherson Treo Khí nén - Điện tử (EAS)
Khối lượng không treo ($m_u$) Lớn ($\ge 80\text{ kg}$/cầu) $\rightarrow$ Tải trọng động cao Nhỏ ($\le 40\text{ kg}$/cụm) $\rightarrow$ Tăng độ êm Rất nhỏ $\rightarrow$ Tối ưu hóa phản hồi mặt đường
Đặc tính đàn hồi Tuyến tính, khó điều chỉnh linh hoạt Tuyến tính/phi tuyến nhẹ qua đòn dẫn Biến thiên phi tuyến theo áp suất và tải trọng
Không gian bố trí Chiếm nhiều khoảng sáng gầm Gọn gàng, tối ưu khoang động cơ Phức tạp (cần máy nén, bình tích áp, cụm van)
Chi phí chế tạo & Bảo dưỡng Rất thấp, dễ gia công trong nước Trung bình, tối ưu hóa sản xuất quy mô lớn Rất cao, đòi hỏi ECU điều khiển thời gian thực
Độ ổn định khi vào cua Hạn chế, dễ văng đuôi Tốt, giảm thiểu góc nghiêng thân xe Tối ưu tuyệt đối (hỗ trợ chuyển chế độ Sport/Comfort)

Hệ thống treo trên xe MEFA5-LAVI-304N được tối ưu hóa theo mô hình lai: Cầu trước sử dụng treo độc lập lò xo trụ tích hợp giảm xóc thủy lực hai lớp vỏ để tăng độ êm dịu cho khoang lái; cầu sau sử dụng treo phụ thuộc kết hợp thanh ổn định ngang nhằm tối ưu hóa khả năng chịu tải của xe khách đô thị.

Thiết kế hệ thống và Cơ sở toán học

Mô hình không gian trạng thái của hệ dao động được xây dựng từ phương trình động lực học tổng quát theo định luật II Newton và nguyên lý D'Alembert:

$$[M]{\ddot{q}} + [C]{\dot{q}} + [K]{q} = {F_{stim}}$$

Chuyển đổi sang hệ phương trình trạng thái với vector trạng thái bậc 8: $x = \begin{bmatrix} q \ \dot{q} \end{bmatrix}_{8 \times 1}$:

$$\begin{cases} \dot{x} = A x + B u \ y = C x + D u \end{cases}$$

Trong đó các ma trận thành phần được tính toán theo công thức cấu trúc:

$$A = \begin{bmatrix} [0]{4\times 4} & [I]{4\times 4} \ -[M]^{-1}[K] & -[M]^{-1}[C] \end{bmatrix}{8 \times 8}, \quad B = \begin{bmatrix} [0]{4\times 4} \ [M]^{-1} \end{bmatrix}_{8 \times 4}$$

$$C = \begin{bmatrix} [I]{4\times 4} & [0]{4\times 4} \ [0]{4\times 4} & [I]{4\times 4} \ -[M]^{-1}[K] & -[M]^{-1}[C] \end{bmatrix}{12 \times 8}, \quad D = \begin{bmatrix} [0]{4\times 4} \ [0]{4\times 4} \ [M]^{-1} \end{bmatrix}{12 \times 4}$$

Công nghệ sử dụng và Thông số định lượng

  • Nền tảng mô phỏng: MATLAB R14 (v7.1) / Simulink Dynamic System Toolbox.
  • Thông số động học xe MEFA5-LAVI-304N:
    • Khối lượng phần được treo: $m_s = 1555\text{ kg}$.
    • Khối lượng không được treo (trước/sau): $m_{u1} = m_{u2} = 50\text{ kg}$.
    • Mô men quán tính lắc dọc thân xe: $J_y = 1809\text{ kg}\cdot\text{m}^2$.
    • Mô men quán tính lắc ngang thân xe: $J_x = 872\text{ kg}\cdot\text{m}^2$.
    • Mô men quán tính khối lượng không treo trục ngang: $I_u = 250\text{ kg}\cdot\text{m}^2$.
    • Chiều dài cơ sở chia theo trọng tâm: $a = 1.29\text{ m}$ (trục trước), $b = 1.06\text{ m}$ (trục sau).
    • Khoảng cách vết bánh xe: $2d_q = 2 \times 1.23\text{ m} = 2.46\text{ m}$; khoảng cách gối đỡ lò xo: $2d = 2 \times 1.71\text{ m} = 3.42\text{ m}$.
    • Độ cứng lốp xe: $K_t = 140000\text{ N/m}$; hệ số cản của lốp: $C_t = 0\text{ N}\cdot\text{s/m}$.

Implementation và kết quả

Quy trình xây dựng mã nguồn (MATLAB/Simulink)

Mã nguồn MATLAB tự động hóa quá trình thiết lập ma trận trạng thái, giải bài toán truyền dao động và tính toán độ lệch chuẩn gia tốc nhằm đánh giá độ êm dịu:

%% File: ctdoc.m - Khao sat anh huong cua do cung Ks tren mat phang doc
global tsim;
ms = 1555;      % Khoi luong duoc treo (kg)
mu1 = 50; mu2 = 50; % Khoi luong khong treo (kg)
Jy = 1809;      % Momen quan tinh quanh truc Oy (kg.m^2)
a = 1.29; b = 1.06; % Toa do trong tam den hai cau (m)
V = 10;         % Van toc chuyen dong (m/s)

% Mang khao sat do cung lo xo cau truoc
ks_range = [10000:2000:25000]; 
kt1 = 140000; kt2 = 140000; % Do cung lop (N/m)
cs1 = 1000; cs2 = 1040;     % He so giam chan thiet ke (N.s/m)
ct1 = 0; ct2 = 0;

for i = 1:length(ks_range)
    ks1 = ks_range(i);
    ks2 = ks1 + 3626; % Quan he do cung giua cau sau va cau truoc
    
    % Thiet lap ma tran khoi luong
    Md = diag([mu1, mu2, ms, Jy]);
    
    % Ma tran can nhot Cd
    Cd = [(cs1+ct1),       0,          -cs1,             a*cs1;
                 0, (cs2+ct2),         -cs2,            -b*cs2;
              -cs1,      -cs2,    (cs1+cs2),     (b*cs2-a*cs1);
             a*cs1,    -b*cs2, (b*cs2-a*cs1), (a^2*cs1 + b^2*cs2)];
         
    % Ma tran do cung Kd
    Kd = [(ks1+kt1),       0,          -ks1,             a*ks1;
                 0, (ks2+kt2),         -ks2,            -b*ks2;
              -ks1,      -ks2,    (ks1+ks2),     (b*ks2-a*ks1);
             a*ks1,    -b*ks2, (b*ks2-a*ks1), (a^2*ks1 + b^2*ks2)];
         
    % Ma tran khong gian trang thai
    Ad = [zeros(4), eye(4); -inv(Md)*Kd, -inv(Md)*Cd];
    Bd = [zeros(4); inv(Md)];
    Cd_out = [eye(4), zeros(4); zeros(4), eye(4); -inv(Md)*Kd, -inv(Md)*Cd];
    Dd_out = [zeros(4); zeros(4); inv(Md)];
    
    % Chay mo phong Simulink
    tsim = 10;
    my_opt = simset('InitialStep', 0.01, 'OutputVariables', 'txy');
    [t, x, y] = sim('doc', [0, tsim], my_opt);
    
    % Tinh toan sai phuong / do lech chuan gia toc
    var_acc_Zs(i) = std(y(:, 11), 1);   % Gia toc tinh tien dung Z''_s
    var_acc_Phi(i) = std(y(:, 12), 1);  % Gia toc goc lac doc phi''_y
end

Kết quả đo lường và Kiểm chứng theo chuẩn ISO/DIS 2631

Các giá trị kiểm thử động lực học được tổng hợp dựa trên miền giá trị tham số danh định của dòng xe MEFA5-LAVI-304N:

Mặt phẳng khảo sát Cấu hình hệ treo Giá trị tối ưu $K_s$ Giá trị tối ưu $C_s$ Sai phương gia tốc đứng $\text{Var}(\ddot{Z}_s)$ Sai phương gia tốc góc $\text{Var}(\ddot{\varphi})$ Đánh giá theo ISO 2631
Mặt phẳng dọc Độc lập trước / Phụ thuộc sau $18092\text{ N/m}$ $1000\text{ N}\cdot\text{s/m}$ $4.4 \times 10^{-4}\text{ m/s}^2$ $2.5 \times 10^{-4}\text{ rad/s}^2$ Đạt chuẩn (Êm dịu cao)
Mặt phẳng ngang Treo độc lập (Cầu trước) $18092\text{ N/m}$ $1000\text{ N}\cdot\text{s/m}$ $4.9 \times 10^{-4}\text{ m/s}^2$ $1.8 \times 10^{-4}\text{ rad/s}^2$ Đạt chuẩn (Ổn định lái)
Mặt phẳng ngang Treo phụ thuộc (Cầu sau) $21718\text{ N/m}$ $1040\text{ N}\cdot\text{s/m}$ $5.1 \times 10^{-4}\text{ m/s}^2$ $2.1 \times 10^{-4}\text{ rad/s}^2$ Đạt chuẩn (Chịu tải tốt) 
Biểu đồ phân bố sai phương gia tốc theo độ cứng Ks và giảm chấn Cs:
Gia tốc Z''s (m/s^2)
 ^
 |          * (Ks = 10000 N/m - Dao động mạnh)
 |           \
 |            \
 |             *---* (Vùng tối ưu: Ks = 18092 N/m, Cs = 1000 Ns/m -> Var = 4.4e-4 m/s^2)
 |                  \
 |                   * (Ks = 30000 N/m - Va đập cứng)
 +-------------------------------------------------------------> Độ cứng Ks (N/m)

Đổi mới và đóng góp

  1. Khái quát hóa mô hình không gian trạng thái đa bậc tự do: Thay vì giải trực tiếp các hệ phương trình vi phân phi tuyến bậc hai cồng kềnh, đề tài đã chuẩn hóa toàn bộ mô hình toán học trên cả mặt phẳng dọc và ngang về dạng không gian trạng thái cấp 8 ($\mathbb{R}^{8 \times 8}$). Phương pháp này cho phép tích hợp trực tiếp vào khối State-Space tiêu chuẩn của Simulink, rút ngắn thời gian tính toán mô phỏng đến $65%$ so với phương pháp số trị truyền thống (Runge-Kutta thuần túy trong C/C++).

  2. Đánh giá tách biệt và toàn diện cho từng cấu hình cầu xe: Phân tích độc lập hệ thống treo độc lập phía trước và phụ thuộc phía sau giúp xác định chính xác tỷ lệ phân bổ độ cứng:

    $$\frac{K_{s2}}{K_{s1}} = \frac{21718}{18092} \approx 1.20$$

    Tỷ lệ này đảm bảo tần số dao động riêng của cầu sau cao hơn cầu trước, giúp triệt tiêu hiện tượng cộng hưởng lắc dọc (pitch resonance) khi xe di chuyển qua các mấp mô liên tiếp.

  3. Bộ công cụ tham số hóa cho công nghiệp ô tô nội địa: Cung cấp trọn vẹn bộ file kịch bản MATLAB (ctdoc.m, ctdoccs.m, ngangdl.m, ngangdlcs.m, ctngangpt.m, ngangptcs.m), tạo tiền đề số liệu thực nghiệm cho dây chuyền sản xuất xe minibus thuộc Dự án KC.13.

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

Trường hợp ứng dụng thực tế

  • Kiểm định và tối ưu hóa kết cấu gầm ô tô: Áp dụng trực tiếp trong giai đoạn tiền khả thi (pre-production) để lựa chọn kích thước dây lò xo thép hợp kim và đường kính lỗ tiết lưu van giảm chấn thủy lực hai chiều.
  • Thiết kế hệ thống điều khiển bán chủ động (Semi-Active Suspension): Cung cấp mô hình toán gốc để tích hợp van điện từ thay đổi độ mở tiết lưu hoặc giảm chấn lưu chất từ tính (MR Damper) có phản hồi thời gian thực từ cảm biến gia tốc.

Lộ trình tích hợp công nghệ

  1. Giai đoạn 1 (Thiết kế & Tối ưu hóa số): Mô phỏng xác thực thông số lò xo trụ ($K_{s1} = 18092\text{ N/m}$, $K_{s2} = 21718\text{ N/m}$) và giảm chấn hai lớp vỏ ($C_{s1} = 1000\text{ N}\cdot\text{s/m}$, $C_{s2} = 1040\text{ N}\cdot\text{s/m}$).
  2. Giai đoạn 2 (Chế tạo mẫu thử & Đo đạc sàn rung): Lắp đặt cảm biến gia tốc áp điện trên khối lượng được treo ($m_s$) và không được treo ($m_u$) trên bệ thử 4 trụ thủy lực để hiệu chuẩn sai số mô hình ($\le 5%$).
  3. Giai đoạn 3 (Sản xuất loạt & Chuyển giao công nghệ): Tích hợp vào quy trình công nghệ lắp ráp chuyển giao từ tập đoàn ChangAn.

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

Hạn chế kỹ thuật

  • Mô hình lốp xe đơn giản: Coi lốp là phần tử đàn hồi tuyến tính thuần túy với độ cứng $K_t = 140000\text{ N/m}$, bỏ qua đặc tính trượt ngang phi tuyến (Mô hình lốp ma thuật Pacejka) và hệ số cản nhớt nội tại của cao su ($C_t \approx 0$).
  • Điều kiện kích động đơn kênh: Thử nghiệm tập trung trên biên dạng mấp mô hình sin đơn giản ($A = 0.025\text{ m}$), chưa bao quát hết phổ mật độ công suất mặt đường ngẫu nhiên (Power Spectral Density - PSD) theo tiêu chuẩn ISO 8608 từ loại A đến loại E.

Hướng phát triển đề xuất

  • Mở rộng mô hình động lực học không gian toàn xe 7 bậc tự do (Full-car 7-DOF) kết hợp động học hệ thống lái và phanh.
  • Ứng dụng giải thuật điều khiển tối ưu hiện đại như LQR (Linear Quadratic Regulator) hoặc H-infinity ($\mathcal{H}_\infty$) để điều khiển hệ thống treo khí nén điện tử, tự động điều chỉnh độ cao gầm xe khi vượt qua ngưỡng tốc độ $80\text{ km/h}$ và $140\text{ km/h}$.

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

  • Sinh viên & Học viên ngành Cơ khí Ô tô / Cơ điện tử: Tiếp cận phương pháp thiết lập phương trình vi phân chuyển động D'Alembert và kỹ thuật giải hệ phương trình không gian trạng thái đa biến trên Simulink.
  • Kỹ sư R&D Ô tô: Sở hữu bộ khung tính toán định lượng ($m_s, J_x, J_y, K_s, C_s$) có sẵn cấu trúc mã nguồn để tùy biến cho các dòng xe thương mại khác nhau.
  • Doanh nghiệp sản xuất, lắp ráp xe trong nước: Tiết kiệm chi phí thử nghiệm thực địa hàng trăm triệu đồng thông qua quy trình tối ưu hóa ảo trước khi gia công cơ khí chính xác.
  • Nhà nghiên cứu học thuật: Nền tảng toán học tin cậy để triển khai các đề tài chuyên sâu về hệ thống treo chủ động (Active Suspension System), hệ thống kiểm soát chống lật (Active Anti-Roll Bar).

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

1. Yêu cầu cấu hình phần cứng và phần mềm để chạy mô hình mô phỏng này là gì?

Hệ thống yêu cầu cài đặt phần mềm MATLAB R14 (v7.1) trở lên tích hợp Simulink Blockset. Do mô hình sử dụng phương pháp không gian trạng thái tuyến tính, mã nguồn có thể thực thi mượt mà trên các cấu hình máy tính cá nhân tiêu chuẩn (CPU x86/x64, RAM từ 2GB trở lên).

2. Vì sao phải phân chia nghiên cứu hệ thống treo trên hai mặt phẳng dọc và ngang độc lập?

Do kết cấu xe MEFA5-LAVI-304N có khối lượng phân bố đối xứng qua mặt phẳng dọc giữa xe, các dao động lắc dọc (pitching) và lắc ngang (rolling) có mức độ tương tác chéo rất nhỏ ở dải vận tốc khảo sát $V = 10\text{ m/s}$. Việc tách biệt thành hai mô hình phẳng giúp giảm bậc tự do, nâng cao độ chính xác số trị và làm rõ cơ chế tác động của từng thông số giảm xóc.

3. Tỷ lệ sai lệch giữa mô phỏng Simulink và thực nghiệm đường thực tế là bao nhiêu?

Với các giả thiết vật rắn tuyệt đối và kích động điều hòa, sai số giữa mô phỏng số và đo đạc thực tế trên sàn rung đối với gia tốc thân xe dao động trong khoảng $5% - 8%$, hoàn toàn nằm trong giới hạn cho phép của thiết kế kỹ thuật ô tô sơ bộ.

4. Tại sao độ cứng lò xo cầu sau ($K_{s2} = 21718\text{ N/m}$) lại được chọn lớn hơn cầu trước ($K_{s1} = 18092\text{ N/m}$)?

Cầu sau của dòng xe minibus chở khách 7–8 chỗ phải chịu tải trọng tĩnh và tải trọng động biến thiên lớn khi đủ tải (toàn bộ hành khách và khoang hành lý phía sau). Độ cứng lớn hơn giúp duy trì khoảng sáng gầm xe, hạn chế kịch hành trình giảm xóc và tối ưu hóa tính năng êm dịu khi đầy tải.

5. Có thể áp dụng mã nguồn này cho xe du lịch con (Sedan/SUV) không?

Hoàn toàn có thể. Người kỹ sư chỉ cần thay đổi các tham số đầu vào trong ma trận khối lượng $M$, khoảng cách hình học ($a, b, d, d_q$) và thông số lò xo, giảm xóc của dòng xe mục tiêu mà không cần tái cấu trúc lại hệ phương trình trạng thái.

Kết luận

Đồ án tốt nghiệp "Mô Phỏng Hệ Thống Treo Ô Tô MEFA5-LAVI-304N Bằng Simulink" của tác giả Nguyễn Xuân Toàn đã giải quyết trọn vẹn bài toán động lực học dao động phức tạp trên phương tiện vận tải hành khách nội địa. Bằng việc kết hợp chặt chẽ giữa cơ sở cơ học giải tích, phương pháp mô hình hóa không gian trạng thái hiện đại và công cụ tính toán số MATLAB/Simulink, đề tài đã xác lập thành công bộ thông số tối ưu cho hệ thống treo: độ cứng cầu trước $K_{s1} = 18092\text{ N/m}$, cầu sau $K_{s2} = 21718\text{ N/m}$, hệ số cản giảm chấn $C_{s1} = 1000\text{ N}\cdot\text{s/m}$, $C_{s2} = 1040\text{ N}\cdot\text{s/m}$. Kết quả đạt được đáp ứng nghiêm ngặt tiêu chuẩn độ êm dịu ISO/DIS 2631, đóng góp giá trị khoa học và thực tiễn trực tiếp cho Dự án cấp Nhà nước KC.13, mở ra tiền đề quan trọng cho việc phát triển các thế hệ xe thương mại thông minh tiếp theo tại Việt Nam.