Giới thiệu dự án

Nghiên cứu động lực học và khảo sát độ êm dịu chuyển động (Ride Comfort Analysis) là một trong những bài toán trọng tâm hàng đầu trong kỹ nghệ phát triển ô tô hiện đại (Automotive NVH Engineering). Theo các báo cáo từ Hiệp hội Kỹ sư Ô tô Quốc tế (SAE International), các tiêu chí liên quan đến rung động và tiếng ồn (Noise, Vibration, and Harshness - NVH) chiếm tới hơn 30% quyết định mức độ hài lòng và trải nghiệm tiện nghi của người sử dụng xe. Tại thị trường Việt Nam, phân khúc xe bán tải (Pickup Truck) chứng kiến sự gia tăng nhanh chóng nhờ khả năng vận hành đa dụng, vừa phục vụ vận chuyển hàng hóa tải trọng lớn, vừa đáp ứng nhu cầu di chuyển gia đình trên đa dạng địa hình.

                  +----------------------------------------------+
                  |           KÍCH ĐỘNG MẶT ĐƯỜNG                |
                  |  - Sóng điều hòa (Harmonic Sine)             |
                  |  - Bán bình phương hàm sin (Half-sine Bump)  |
                  +----------------------+-----------------------+
                                         |
                                         v
                  +----------------------------------------------+
                  |         BÁNH XE & KHỐI LƯỢNG KĐT             |
                  |    M1 (Cầu trước), M2 (Cầu sau), CL1, CL2    |
                  +----------------------+-----------------------+
                                         |
                                         v
                  +----------------------------------------------+
                  |          HỆ THỐNG TREO (C, K)                |
                  |  - Trước: Độc lập tay đòn kép (C1, K1)       |
                  |  - Sau: Phụ thuộc nhíp lá (C2, K2)           |
                  +----------------------+-----------------------+
                                         |
                                         v
                  +----------------------------------------------+
                  |       THÂN XE - KHỐI LƯỢNG ĐƯỢC TREO         |
                  |    Dịch chuyển thẳng đứng (Z), Lắc dọc (φ)   |
                  +----------------------+-----------------------+
                                         |
                                         v
                  +----------------------------------------------+
                  |       ĐÁNH GIÁ CHỈ TIÊU ÊM DỊU (CAE)         |
                  |  - Tần số dao động riêng f0 (Hz)             |
                  |  - Gia tốc hiệu dụng toàn phương RMS (m/s²)  |
                  |  - Tiêu chuẩn ISO 2631-1 / TCVN 6964-1:2001  |
                  +----------------------------------------------+

Đặc thù kỹ thuật của dòng xe bán tải là kết cấu hệ thống treo sau phụ thuộc sử dụng nhíp lá (Leaf Spring) có độ cứng cao nhằm chịu tải trọng nặng, kết hợp hệ thống treo trước độc lập (Independent Front Suspension). Khi xe vận hành ở trạng thái không tải hoặc tải nhẹ qua các bề mặt đường gồ ghề, dao động kích thích từ mặt đường truyền trực tiếp qua lốp xe, cơ cấu treo và tác động tiêu cực đến thân xe. Những rung động cưỡng bức này không chỉ gây mệt mỏi, suy giảm khả năng tập trung của người lái mà còn ảnh hưởng đến hệ thần kinh, cột sống và gia tăng nguy cơ mất an toàn giao thông theo tiêu chuẩn ISO 2631-1:1997.

Mục tiêu của dự án

  1. Xây dựng mô hình toán học: Thiết lập hệ phương trình vi phân chuyển động 4 bậc tự do (4-DOF Pitch-Bounce Model) mô tả dao động ô tô trong mặt phẳng dọc dựa trên nguyên lý D'Alembert.
  2. Số hóa mô hình mô phỏng: Lập trình giải thuật tính toán số trên phần mềm MATLAB và trực quan hóa sơ đồ khối động lực học trên Simulink.
  3. Thu thập thông số thực tế: Đo đạc và trích xuất bộ thông số động học, hình học và kết cấu thực tế của xe Mitsubishi Triton từ trung tâm dịch vụ kỹ thuật ủy quyền Mitsubishi Nam Auto.
  4. Kiểm chứng đa nền tảng (Cross-Verification): Xây dựng mô hình động lực học đa vật thể (Multi-Body Dynamics - MBD) trên phần mềm thương mại CarSim để đối chuẩn và thẩm định độ tin cậy của thuật toán xây dựng trên MATLAB-Simulink.
  5. Đánh giá tiêu chuẩn êm dịu: Định lượng độ êm dịu thông qua giá trị tần số dao động riêng ($f_0$) và bình phương trung bình gia tốc dịch chuyển ($a_{RMS}$), đối chiếu với tiêu chuẩn TCVN 6964-1:2001.

Phạm vi và giới hạn nghiên cứu

  • Đối tượng khảo sát: Xe bán tải Mitsubishi Triton thế hệ trang bị hệ thống treo trước độc lập lò xo xoắn và hệ thống treo sau phụ thuộc nhíp lá.
  • Phạm vi mô hình hóa: Không gian 2D mặt phẳng dọc ($OXZ$), tập trung vào hai dạng dao động chính là dịch chuyển thẳng đứng (Bounce - $Z$) và lắc dọc thân xe (Pitch - $\varphi$). Bỏ qua dao động lắc ngang (Roll) và chuyển hướng (Yaw).
  • Kích động mặt đường: Mô hình hóa toán học với hai dạng biên dạng chuẩn: kích động điều hòa hình sin liên tục và kích động dạng va đập xung nhọn bán bình phương hàm sin (Bump road).

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

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

Trong quy trình phát triển phương tiện truyền thống, đánh giá độ êm dịu thường phụ thuộc vào các bài kiểm tra thực nghiệm tại đường thử chuyên dụng (Proving Ground). Phương pháp này đòi hỏi chi phí chế tạo mẫu thử (Physical Prototypes) lớn, thiết bị cảm biến đo lường đắt tiền và kéo dài thời gian phát triển sản phẩm.

Phương pháp đánh giá Ưu điểm Nhược điểm Chi phí & Thời gian
Thực nghiệm đường thử (Physical Proving Ground Testing) Độ chính xác thực tế cao nhất, phản ánh đúng toàn bộ tương tác vật lý. Rủi ro hỏng hóc nguyên mẫu, phụ thuộc thời tiết, khó lặp lại chuẩn xác 100% điều kiện kích động. Chi phí rất cao ($50.000 - $200.000), chu kỳ đo kiểm 3 - 6 tháng.
Mô hình 1D/2-DOF giản lược (Quarter-Car Model) Tính toán siêu nhanh, phương trình đại số đơn giản, dễ giải tích. Không thể hiện được dao động lắc dọc (Pitching), mất tương tác động lực học liên kết cầu trước - cầu sau. Chi phí thấp, độ chính xác động lực học tổng thể xe $< 60%$.
Mô phỏng phối hợp MATLAB-Simulink & CarSim (Đồ án áp dụng) Tối ưu hóa giữa độ chính xác động lực học 4-DOF và tốc độ xử lý số; kiểm chứng chéo trực tiếp; dễ dàng tinh chỉnh độ cứng ($C$) và giảm chấn ($K$). Yêu cầu chuẩn hóa chính xác bộ thông số đầu vào cơ hệ (khối lượng, mômen quán tính, ma trận độ cứng). Chi phí phần mềm thấp, thời gian mô phỏng tính bằng mili-giây, độ chính xác $> 95%$.

Phân loại yêu cầu hệ thống theo mô hình MoSCoW

  • Must-have (Bắt buộc): Hệ phương trình vi phân mô tả 4 bậc tự do của phần khối lượng được treo (Sprung Mass - $M_s$) và hai khối lượng không được treo (Unsprung Mass - $M_1, M_2$); bộ giải ODE tích phân số bậc cao; hàm kích động mặt đường dạng điều hòa và bán bình phương hàm sin; bộ lọc tính chỉ số $a_{RMS}$.
  • Should-have (Nên có): Bộ khối giao diện đồ họa Simulink phân nhánh rõ ràng (Subsystem Cầu trước, Cầu sau, Thân xe); mô hình đối chuẩn 3D đa vật thể trên CarSim; cấu hình xuất tín hiệu sang MATLAB Workspace tự động.
  • Could-have (Có thể mở rộng): Tích hợp thuật toán điều khiển mờ (Fuzzy Logic) hoặc LQR cho hệ thống treo bán tích cực (Semi-active Suspension với giảm chấn từ biến MR Damper).
  • Won't-have (Chưa thực hiện trong giai đoạn này): Mô hình hóa biến dạng phi tuyến khung gầm bằng phần tử hữu hạn (Chassis Flexible FEM) và ma sát trượt giữa các lá nhíp phi tuyến đa lớp.

Thiết kế hệ thống

Kiến trúc mô phỏng khảo sát độ êm dịu được thiết kế theo cấu trúc khối mô-đun hóa, đảm bảo tính liên kết chặt chẽ giữa toán học lý thuyết, mô hình khối chức năng và môi trường mô phỏng vật lý chuyên sâu.

Technology Stack & Thông số kỹ thuật

  • Môi trường tính toán số: MATLAB R2023b (MathWorks), sử dụng Scripting Engine và Symbolic Math Toolbox để thiết lập ma trận khối lượng $[M]$, ma trận cản $[K]$, ma trận độ cứng $[C]$.
  • Môi trường mô phỏng động học: Simulink Version 10.8 với thư viện Continuous, Math Operations, Signal Routing, và Sinks.
  • Phần mềm kiểm chứng chuyên sâu: CarSim Version 2020.0 (Mechanical Simulation Corporation) thiết lập xe cơ sở Pick-up Truck với hệ thống lái, lốp và khí động học tương ứng.
  • Bộ thông số hình học và động học xe Mitsubishi Triton:
Thông số ký hiệu Diễn giải vật lý Giá trị danh định Đơn vị tính
$M_s$ Khối lượng được treo (Sprung Mass) $1850$ $\text{kg}$
$M_1$ Khối lượng không được treo cầu trước $115$ $\text{kg}$
$M_2$ Khối lượng không được treo cầu sau $145$ $\text{kg}$
$J_y$ Mômen quán tính lắc dọc thân xe $2850$ $\text{kg}\cdot\text{m}^2$
$C_1, C_2$ Độ cứng hệ thống treo trước / sau $36000 ;/; 48000$ $\text{N/m}$
$K_1, K_2$ Hệ số cản giảm chấn treo trước / sau $2400 ;/; 3100$ $\text{N}\cdot\text{s/m}$
$C_{L1}, C_{L2}$ Độ cứng hướng kính của lốp trước / sau $220000 ;/; 220000$ $\text{N/m}$
$L, a, b$ Chiều dài cơ sở, khoảng cách từ trọng tâm tới cầu trước / sau $3.00, 1.38, 1.62$ $\text{m}$

Methodology

Quy trình nghiên cứu áp dụng mô hình chuẩn kỹ thuật V-Model kết hợp các pha kiểm chứng lặp (Iterative CAE Validation Framework):

     Xác định yêu cầu (TCVN 6964-1)              Đánh giá tiêu chuẩn Êm dịu
             \                                              /
        Khảo sát thông số xe Triton               Kiểm chứng đa nền tảng CarSim
               \                                      /
           Xây dựng Mô hình toán (4-DOF)      Tích phân tín hiệu & Tính a_RMS
                   \                              /
                   Thiết lập MATLAB & Simulink Model
  • Phase 1 (Tuần 1 - 3): Thu thập dữ liệu thực tế tại xưởng Mitsubishi Nam Auto, tổng hợp cơ sở lý thuyết dao động ô tô và tiêu chuẩn rung động toàn thân ISO 2631.
  • Phase 2 (Tuần 4 - 7): Giải phóng liên kết cơ hệ, áp dụng nguyên lý D'Alembert thiết lập hệ phương trình vi phân chuyển động 4 bậc tự do; chuẩn hóa hàm kích động mặt đường.
  • Phase 3 (Tuần 8 - 11): Lập trình giải thuật số trên MATLAB Script; xây dựng cấu trúc khối tích phân trong Simulink; tiến hành chạy mô phỏng khảo sát chuyển vị, vận tốc và gia tốc.
  • Phase 4 (Tuần 12 - 14): Cấu hình cơ hệ chi tiết trên CarSim; thiết lập đường thử mấp mô ảo tương đương; trích xuất dữ liệu dao động và kiểm tra chéo (Cross-check).
  • Phase 5 (Tuần 15 - 16): Xử lý hậu kỳ số liệu (Post-processing), tính toán các giá trị hiệu dụng gia tốc trọng số RMS và hoàn thiện đồ án.

Implementation và kết quả

Development Process & Mathematical Formulation

Dựa trên nguyên lý cân bằng động D'Alembert cho cơ hệ chuyển động trong mặt phẳng thẳng đứng:

$$\sum \mathbf{F}_z = 0, \quad \sum \mathbf{M}_y = 0$$

Hệ phương trình vi phân chuyển động 4 bậc tự do của xe ô tô hai cầu được thiết lập như sau:

  1. Phương trình chuyển động của khối lượng không được treo cầu trước ($M_1$): $$M_1 \ddot{E}1 + (K_1 + K{L1})\dot{E}1 - K_1 \dot{Z}1 + (C_1 + C{L1})E_1 - C_1 Z_1 = K{L1}\dot{q}1 + C{L1}q_1$$

  2. Phương trình chuyển động của khối lượng không được treo cầu sau ($M_2$): $$M_2 \ddot{E}2 + (K_2 + K{L2})\dot{E}2 - K_2 \dot{Z}2 + (C_2 + C{L2})E_2 - C_2 Z_2 = K{L2}\dot{q}2 + C{L2}q_2$$

  3. Phương trình dịch chuyển thẳng đứng của khối lượng được treo ($M_s$): $$M_s \ddot{Z} + (K_1 + K_2)\dot{Z} + (C_1 + C_2)Z - (a K_1 - b K_2)\dot{\varphi} - (a C_1 - b C_2)\varphi - K_1 \dot{E}_1 - C_1 E_1 - K_2 \dot{E}_2 - C_2 E_2 = 0$$

  4. Phương trình chuyển động lắc dọc thân xe quanh trục trọng tâm ($J_y$): $$J_y \ddot{\varphi} + (a^2 K_1 + b^2 K_2)\dot{\varphi} + (a^2 C_1 + b^2 C_2)\varphi - (a K_1 - b K_2)\dot{Z} - (a C_1 - b C_2)Z + a K_1 \dot{E}_1 + a C_1 E_1 - b K_2 \dot{E}_2 - b C_2 E_2 = 0$$

Dưới dạng ma trận trạng thái: $$[M]{\ddot{x}} + [K]{\dot{x}} + [C]{x} = {F(t)}$$

Trong đó: $${x} = \begin{bmatrix} E_1 & E_2 & Z & \varphi \end{bmatrix}^T$$

%==========================================================================
% Chuong trinh tinh toan dao dong va tri so RMS gia toc oto 4-DOF
% Doi tuong khao sat: Mitsubishi Triton
%==========================================================================
clc; clear; close all;

% 1. Khoi tao thong so he thong
Ms = 1850;      % Khoi luong duoc treo (kg)
M1 = 115;       % Khoi luong khong duoc treo cau truoc (kg)
M2 = 145;       % Khoi luong khong duoc treo cau sau (kg)
Jy = 2850;      % Momen quan tinh than xe (kg.m^2)
C1 = 36000;     % Do cung treo truoc (N/m)
C2 = 48000;     % Do cung treo sau (N/m)
K1 = 2400;      % He so giam chan treo truoc (N.s/m)
K2 = 3100;      % He so giam chan treo sau (N.s/m)
CL1 = 220000;   % Do cung lop truoc (N/m)
CL2 = 220000;   % Do cung lop sau (N/m)
a = 1.38;       % Khoang cach tu trong tam den cau truoc (m)
b = 1.62;       % Khoang cach tu trong tam den cau sau (m)
L = a + b;      % Chieu dai co so (m)
v_kmh = 40;     % Van toc xe (km/h)
v = v_kmh / 3.6;% Quy doi sang m/s

% 2. Thiet lap kich dong mat duong (Dang ban binh phuong hinh sin)
q0 = 0.08;      % Chieu cao go giam toc (m)
S = 2.5;        % Chieu dai go giam toc (m)
t_pass = S / v; % Thoi gian banh truoc vuot qua go
dt_delay = L / v; % Do tre thoi gian giua cau truoc va cau sau

% 3. Thiet lap he ma tran dong luc hoc
M_mat = diag([M1, M2, Ms, Jy]);

K_mat = [ (K1),       0,          -K1,         a*K1;
           0,         (K2),       -K2,        -b*K2;
          -K1,       -K2,         (K1+K2),    -(a*K1-b*K2);
           a*K1,     -b*K2,      -(a*K1-b*K2), (a^2*K1+b^2*K2) ];

C_mat = [ (C1+CL1),   0,          -C1,         a*C1;
           0,         (C2+CL2),   -C2,        -b*C2;
          -C1,       -C2,         (C1+C2),    -(a*C1-b*C2);
           a*C1,     -b*C2,      -(a*C1-b*C2), (a^2*C1+b^2*C2) ];

% 4. Tinh toan tan so dao dong rieng f0 cua than xe
A_sys = [zeros(4,4), eye(4); -inv(M_mat)*C_mat, -inv(M_mat)*K_mat];
eigenvalues = eig(A_sys);
wn = imag(eigenvalues(imag(eigenvalues) > 0));
fn = wn / (2*pi);

fprintf('--- KET QUA TAN SO DAO DONG RIENG ---\n');
fprintf('Tan so dao dong than xe f_bounce: %.2f Hz\n', fn(1));
fprintf('Tan so dao dong than xe f_pitch:  %.2f Hz\n', fn(2));
fprintf('Tan so dao dong cau truoc:        %.2f Hz\n', fn(3));
fprintf('Tan so dao dong cau sau:          %.2f Hz\n', fn(4));

% 5. Tinh toan gia tri hieu dung toan phuong a_RMS (Simulink/ODE output)
% a_RMS = sqrt( (1/T) * trapz(t, z_ddot.^2) )

Testing và validation

Mô hình được thử nghiệm trên 2 kịch bản vận hành tiêu chuẩn:

  • Scenario 1: Vận tốc $v = 40 \text{ km/h}$ di chuyển qua gờ giảm tốc dạng bán bình phương hàm sin chiều cao $q_0 = 0.08 \text{ m}$, chiều dài $S = 2.5 \text{ m}$.
  • Scenario 2: Vận tốc $v = 60 \text{ km/h}$ di chuyển trên mặt đường gợn sóng điều hòa biên độ $q_0 = 0.02 \text{ m}$, bước sóng $\lambda = 5.0 \text{ m}$.
                 KẾT QUẢ ĐỐI CHUẨN GIA TỐC THÂN XE (SCENARIO 1)
   Gia tốc z" (m/s²)
       3.0 |              /\ (Peak: 2.38 m/s²)
           |             /  \
       1.5 |            /    \   --- MATLAB/Simulink
           |  ---------/------\----- CarSim Validation
       0.0 |          /        \________
           |         /                  \________
      -1.5 |________/                            \_____________
           +---------------------------------------------------->
           0.0       0.5       1.0       1.5       2.0    Thời gian (s)

Bảng tổng hợp kết quả mô phỏng và kiểm chứng CarSim

Đại lượng khảo sát MATLAB-Simulink CarSim 2020.0 Độ lệch tương đối (%) Đánh giá tiêu chuẩn
Biên độ chuyển vị lớn nhất ($Z_{max}$) $0.0542 \text{ m}$ $0.0561 \text{ m}$ $3.38%$ Khớp hoàn toàn xu hướng dập tắt dao động
Gia tốc thân xe cực đại ($\ddot{Z}_{max}$) $2.38 \text{ m/s}^2$ $2.47 \text{ m/s}^2$ $3.64%$ Xuất hiện đỉnh tại thời điểm $t = 0.28\text{s}$
Góc lắc dọc cực đại ($\varphi_{max}$) $0.038 \text{ rad}$ $0.039 \text{ rad}$ $2.56%$ Triệt tiêu hoàn toàn sau $1.4\text{s}$
Tần số dao động riêng thân xe ($f_{01}$) $1.28 \text{ Hz}$ $1.25 \text{ Hz}$ $2.40%$ Nằm trong vùng tối ưu ($1.0 - 1.5\text{ Hz}$)
Gia tốc hiệu dụng trọng số ($a_{RMS}$) $0.485 \text{ m/s}^2$ $0.508 \text{ m/s}^2$ $4.52%$ Mức "Có cảm giác chút ít không thoải mái" theo TCVN 6964-1

Kết quả đạt được

  1. Độ chính xác mô hình: Mô hình toán học 4-DOF trên MATLAB-Simulink đạt mức tương quan cao ($R^2 > 0.96$) so với phần mềm thương mại CarSim, sai số các thông số dao động chính đều $< 5%$.
  2. Định lượng độ êm dịu:
    • Tần số dao động riêng của phần khối lượng được treo đạt $1.28 \text{ Hz}$ đối với nhún đứng và $1.46 \text{ Hz}$ đối với lắc dọc. Các giá trị này nằm hoàn hảo trong dải tần số thích ứng sinh lý học tự nhiên của cơ thể người khi đi bộ ($1.0 - 1.5 \text{ Hz}$ tương đương $60 - 90 \text{ chu kỳ/phút}$).
    • Giá trị $a_{RMS}$ tính toán đạt $0.485 \text{ m/s}^2$ trên mặt đường mấp mô xung kích, đáp ứng giới hạn an toàn sức khỏe theo tiêu chuẩn TCVN 6964-1:2001 (ngưỡng bắt đầu mệt mỏi là $0.63 \text{ m/s}^2$).
  3. Thời gian dập tắt dao động: Hệ thống giảm chấn thủy lực kép ($K_1 = 2400 \text{ Ns/m}, K_2 = 3100 \text{ Ns/m}$) giúp dập tắt hoàn toàn rung động cưỡng bức trong khoảng thời gian $t_s = 1.35 \text{ s}$, đảm bảo bánh xe không bị tách khỏi mặt đường, duy trì tối đa điều kiện bám an toàn ($D_s$).

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

  • Quy trình thẩm định mô phỏng liên hợp (Dual-Engine Validation Pipeline): Thay vì chỉ dừng lại ở các phương trình giải tích lý thuyết, đồ án đã hiện thực hóa quy trình cầu nối giữa mô hình tính toán số tự xây dựng (White-box model trên MATLAB-Simulink) và mô hình công nghiệp thương mại (Black-box model trên CarSim).
  • So sánh cải tiến hiệu năng với các giải pháp hiện hành:
Tiêu chí Đo kiểm thực nghiệm đường thử Mô hình 1/4 xe (Quarter-Car) Giải pháp Đồ án (4-DOF Dual Validation)
Khảo sát góc lắc dọc ($\varphi$) Có (Cần nhiều cảm biến IMU) Hoàn toàn không Có (Mô phỏng chính xác tương tác trục)
Thời gian thiết lập & chạy Hàng tuần Dưới $1 \text{ giây}$ Dưới $5 \text{ giây}$
Chi phí triển khai Rất tốn kém ($> 100$ triệu VNĐ) Gần như bằng 0 Tiết kiệm 90% chi phí thử nghiệm sơ bộ
Độ tin cậy kỹ thuật Chuẩn tuyệt đối Kém ($< 60%$) Cao ($> 95%$ so với CarSim)
  • Đóng góp học thuật và công nghiệp: Cung cấp bộ công cụ mã nguồn mở dạng khối (Modular Blocksets) trong Simulink được tham số hóa chuẩn cho dòng xe bán tải, tạo tiền đề cho các kỹ sư phát triển thuật toán điều khiển treo thích ứng (Adaptive/Active Suspension System).

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

+-----------------------------------------------------------------------------------+
|                        LỘ TRÌNH ỨNG DỤNG THỰC TẾ (ROADMAP)                        |
+-----------------------------------------------------------------------------------+
| Giai đoạn 1: Tiền khả thi (0 - 3 tháng)                                           |
| - Tích hợp mô hình vào chương trình giảng dạy Động lực học Ô tô.                  |
| - Xây dựng GUI MATLAB cho phép sinh viên thay đổi thông số C, K trực quan.        |
+-----------------------------------------------------------------------------------+
| Giai đoạn 2: Tối ưu hóa R&D (3 - 6 tháng)                                         |
| - Ứng dụng cho các xưởng dịch vụ nâng cấp phuộc/nhíp xe bán tải (Aftermarket).    |
| - Tự động hóa tính toán độ cứng lò xo và van tiết lưu giảm chấn tối ưu.           |
+-----------------------------------------------------------------------------------+
| Giai đoạn 3: Phát triển Điều khiển Nâng cao (6 - 12 tháng)                         |
| - Thiết kế thuật toán điều khiển bán tích cực Skyhook / Groundhook.               |
| - Nhúng thuật toán vào bộ điều khiển phần cứng ECU/dSPACE (Hardware-in-the-Loop).  |
+-----------------------------------------------------------------------------------+
  • Môi trường ứng dụng sản xuất: Các trung tâm R&D xe chuyên dụng có thể sử dụng bộ mô hình để ước lượng nhanh ảnh hưởng của việc thay đổi tải trọng thùng xe (Payload Variations từ $0 \text{ kg}$ đến $950 \text{ kg}$) lên chỉ số êm dịu mà không cần chế tạo mẫu thử nghiệm vật lý.
  • Tối ưu hóa chi phí (ROI): Việc số hóa hoàn toàn giai đoạn thử nghiệm sơ bộ giúp các đơn vị độ xe và phát triển phụ tùng giảm tới 70% số chu kỳ thử nghiệm phá hủy, tiết kiệm hàng trăm giờ vận hành trên các máy đo xóc chuyên dụng.

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

  • Hạn chế kỹ thuật hiện tại:
    • Giả thiết tuyến tính hóa: Độ cứng lốp xe ($C_{L}$) và hệ số cản ($K$) đang được coi là hằng số tuyến tính, chưa xét đến đặc tính trễ phi tuyến (Hysteresis) của cao su lốp và van giảm chấn đa tầng khi vận tốc nén/trả thay đổi đột ngột.
    • Mô hình không gian 2D bỏ qua các tương tác dao động xoắn thân xe, dao động lắc ngang (Roll) khi xe vừa xóc vừa vào cua.
  • Hướng phát triển tiếp theo:
    • Mở rộng mô hình lên Không gian 7 bậc tự do (Full-Car 7-DOF Model) kết hợp dao động nhún, lắc dọc và lắc ngang.
    • Tích hợp mặt đường ngẫu nhiên theo tiêu chuẩn quốc tế ISO 8608 (từ mặt đường loại A đến loại E) sử dụng phương pháp tạo tín hiệu lọc nhiễu trắng (White Noise Shaping Filter).
    • Nghiên cứu cơ cấu chấp hành hệ thống treo khí nén điện tử kết hợp thuật toán điều khiển tối ưu mờ $H_\infty$ / Skyhook.

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

  • Sinh viên & Học viên cao học ngành Kỹ thuật Ô tô: Tiếp cận trọn vẹn quy trình mô hình hóa toán học chuyển động cơ hệ phức tạp từ nguyên lý D'Alembert đến code hóa MATLAB và mô phỏng khối trên Simulink.
  • Kỹ sư CAE & Động lực học phương tiện (Vehicle Dynamics Engineers): Có sẵn bộ template mô phỏng 4-DOF đã được chuẩn hóa và kiểm chứng với CarSim, sẵn sàng tùy biến cho các dòng xe du lịch hoặc xe tải nhẹ.
  • Doanh nghiệp phụ tùng & Xưởng dịch vụ ô tô: Công cụ định lượng chính xác hiệu quả êm dịu khi nâng cấp các bộ kit phuộc/lò xo độ (Aftermarket Suspension Kits) cho khách hàng dựa trên số liệu đo lường khoa học.
  • Nhà nghiên cứu học thuật: Cung cấp dữ liệu cơ sở tin cậy phục vụ các công trình nghiên cứu sâu hơn về kiểm soát dao động và tối ưu hóa đa mục tiêu hệ thống treo.

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

1. Cần cấu hình phần cứng và phần mềm như thế nào để chạy mô hình mô phỏng này?

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 8GB (khuyến nghị 16GB), bộ cài đặt MATLAB/Simulink phiên bản R2020a trở lên và phần mềm CarSim phiên bản 2019.0 trở lên. Quá trình tính toán trên MATLAB hoàn tất chỉ trong $1 - 3 \text{ giây}$.

2. Tại sao mô hình 4 bậc tự do trong mặt phẳng dọc lại đủ tin cậy để đánh giá độ êm dịu?

Vì khi xe chuyển động thẳng qua các mấp mô ngang đường (như gờ giảm tốc hoặc ổ gà đối xứng), chuyển động thẳng đứng (Bounce) và quay quanh trục ngang qua trọng tâm (Pitch) là 2 dạng dao động chiếm hơn 85% tổng năng lượng kích thích truyền lên cơ thể người ngồi trong xe.

3. Làm thế nào để điều chỉnh mô hình này cho một dòng xe khác (ví dụ: SUV hoặc Sedan)?

Chỉ cần thay đổi tệp thông số đầu vào trong đoạn script MATLAB bao gồm: phân bố khối lượng ($M_s, M_1, M_2$), mômen quán tính $J_y$, tọa độ trọng tâm ($a, b$), độ cứng lò xo ($C_1, C_2$), hệ số cản ($K_1, K_2$) và độ cứng lốp ($C_{L1}, C_{L2}$) tương ứng của xe mục tiêu.

4. Tiêu chuẩn TCVN 6964-1:2001 quy định ngưỡng gia tốc RMS như thế nào đối với độ êm dịu?

Theo tiêu chuẩn, giá trị $a_{RMS} < 0.315 \text{ m/s}^2$ là trạng thái hoàn toàn thoải mái; $0.315 - 0.63 \text{ m/s}^2$ có cảm giác một chút không thoải mái; $0.5 - 1.0 \text{ m/s}^2$ gây khó chịu rõ rệt; và $> 2.0 \text{ m/s}^2$ ở mức cực kỳ khó chịu, có thể gây chấn thương nếu kéo dài.

5. Tại sao có sự chênh lệch nhỏ ($< 5%$) giữa kết quả tính toán của Simulink và CarSim?

Sự chênh lệch phát sinh do CarSim sử dụng mô hình động lực học phi tuyến 3D toàn xe (Full-vehicle Non-linear Multibody Model) có tính đến độ mềm khớp cao su liên kết (Bushings), động học góc đặt bánh xe (Suspension Kinematics elastokinematics) và biến dạng phức tạp của lốp xe, trong khi mô hình Simulink đang xét hệ tuyến tính trong không gian 2D.

Kết luận

Đồ án tốt nghiệp "Khảo sát độ êm dịu trên xe ô tô" (Nghiên cứu điển hình trên xe bán tải Mitsubishi Triton) đã giải quyết thành công bài toán cốt lõi trong kỹ thuật động lực học phương tiện: thiết lập hoàn chỉnh hệ phương trình vi phân chuyển động 4 bậc tự do, số hóa mô hình trên nền tảng tính toán MATLAB-Simulink và thẩm định chéo độ chính xác thông qua phần mềm tiêu chuẩn công nghiệp CarSim.

Kết quả nghiên cứu chứng minh hệ thống treo nguyên bản của Mitsubishi Triton đảm bảo tốt các tiêu chí về tần số dao động tự nhiên ($1.28 \text{ Hz}$) và gia tốc hiệu dụng trọng số ($a_{RMS} = 0.485 \text{ m/s}^2$), tuân thủ nghiêm ngặt các yêu cầu an toàn sinh lý học theo TCVN 6964-1:2001 và ISO 2631-1:1997. Phương pháp luận và bộ công cụ mô phỏng được xây dựng trong đồ án mở ra giải pháp số hóa hiệu quả, giúp rút ngắn thời gian và tối ưu hóa chi phí nghiên cứu phát triển cho ngành công nghiệp kỹ thuật ô tô.