Giới thiệu dự án

Nghiên cứu động lực học và độ êm dịu chuyển động (Ride Comfort) là một trong những bài toán cốt lõi của ngành Công nghệ Kỹ thuật Ô tô. Theo các báo cáo thống kê từ Hiệp hội các nhà sản xuất ô tô (VAMA) và Hiệp hội Kỹ sư Ô tô Quốc tế (SAE), phân khúc xe bán tải (Pickup Truck) chiếm hơn 18% thị phần xe đa dụng nhờ khả năng chuyên chở hàng hóa tải trọng lớn kết hợp di chuyển linh hoạt trong đô thị và địa hình phức tạp. Tuy nhiên, đặc thù kết cấu của dòng xe này là sử dụng hệ thống treo phụ thuộc loại nhíp lá ở cầu sau để chịu tải nặng, dẫn đến độ cứng vững ($C_2$) rất cao. Khi xe vận hành ở trạng thái không tải hoặc tải trọng nhẹ qua các đoạn đường mấp mô, gia tốc rung động thân xe tăng cao đột ngột, gây xóc nảy dữ dội và suy giảm nghiêm trọng độ êm dịu.

+-----------------------------------------------------------------------------------+
|                              BÀI TOÁN CỐT LÕI                                     |
|  Kích động mặt đường q(t) ---> Dao động thân xe (Nhún Z, Lắc dọc phi)             |
|  ---> Tác động rung 4-12 Hz lên cơ thể ---> Mệt mỏi, say xe, giảm an toàn         |
+-----------------------------------------------------------------------------------+

Vấn đề kỹ thuật cụ thể (Problem Statement): Dao động cơ học truyền từ mặt đường qua lốp và hệ thống treo lên thân xe tạo ra dải tần số kích thích từ $1 - 20\text{ Hz}$. Theo y học lao động và động lực học người ngồi trên xe, dải tần số $4 - 12\text{ Hz}$ gây hiện tượng cộng hưởng trực tiếp lên cột sống và hệ tiêu hóa, trong khi tần số $15 - 18\text{ Hz}$ gây rối loạn tiền đình và cảm giác choáng ngợp. Do đó, việc xây dựng mô hình toán học giải tích kết hợp công cụ mô phỏng số để lượng hóa, đánh giá chính xác các chỉ tiêu êm dịu theo tiêu chuẩn là yêu cầu cấp thiết nhằm hỗ trợ quá trình thiết kế, cân chỉnh thông số giảm xóc.

Đồ á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 dòng xe bán tải Mitsubishi Triton) đặt ra 4 mục tiêu trọng tâm:

  1. Thiết lập hệ phương trình vi phân mô tả dao động không gian thu gọn trong mặt phẳng dọc (mô hình nửa xe 4 bậc tự do - 4-DOF Half-Car Model) dựa trên nguyên lý D'Alembert.
  2. Xây dựng thuật toán và chương trình giải số trên môi trường MATLAB script và mô phỏng trực quan trên Simulink.
  3. Ứng dụng phần mềm mô phỏng động lực học đa vật thể chuyên ngành CarSim để mô hình hóa toàn xe và thẩm định độ chính xác (Cross-Validation) của mô hình toán học đã lập.
  4. Đánh giá độ êm dịu của xe dựa trên 2 chỉ tiêu kỹ thuật cốt lõi: Tần số dao động riêng ($f_0$) và Giá trị hiệu dụng của trọng số gia tốc rung động ($a_{RMS}$) theo tiêu chuẩn quốc tế ISO 2631-1:1997 và Tiêu chuẩn Quốc gia TCVN 6964-1:2001.

Phương pháp tiếp cận giải pháp sử dụng mô hình nửa xe 4 bậc tự do liên kết tuyến tính, tối ưu hóa quá trình tính toán động lực học phi tuyến phức tạp nhưng vẫn giữ độ chính xác cao đối với chuyển động nhún thẳng đứng ($Z$) và lắc dọc ($\varphi$). Kết quả mong đợi đạt sai số giữa mô hình lý thuyết MATLAB/Simulink và phần mềm kiểm chứng thương mại CarSim dưới 8%, đồng thời đưa ra các khuyến nghị thông số hệ thống treo tối ưu.

+--------------------------+--------------------------------------------------------+
| Hạng mục                 | Phạm vi nghiên cứu & Giới hạn kỹ thuật                 |
+--------------------------+--------------------------------------------------------+
| Đối tượng nghiên cứu     | Xe bán tải Mitsubishi Triton (Treo trước độc lập, sau  |
|                          | phụ thuộc nhíp lá)                                     |
| Không gian khảo sát      | Mặt phẳng dọc đối xứng (Chuyển vị đứng Z, Lắc dọc phi) |
| Kích động mặt đường      | Dạng điều hòa hình sin và dạng bán bình phương hàm sin |
| Giới hạn mô hình         | Bỏ qua dao động xoắn thân xe, xem lốp tiếp xúc điểm,   |
|                          | góc lắc phi nhỏ (tan(phi) ~ phi)                       |
+--------------------------+--------------------------------------------------------+

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

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

Trong kỹ thuật động lực học ô tô, việc khảo sát độ êm dịu chuyển động thường được thực hiện qua ba phương pháp chính: Thực nghiệm trên đường thử/băng thử rung, Giải tích toán học cổ điển, và Mô phỏng số hóa (Computer-Aided Engineering - CAE).

Tiêu chí so sánh Phương pháp Thực nghiệm Giải tích Toán học Cổ điển Mô hình hóa MATLAB/Simulink kết hợp CarSim
Chi phí triển khai Rất cao (Cần xe thật, cảm biến gia tốc, phòng thử) Rất thấp (Chỉ tính toán lý thuyết trên giấy) Tối ưu (Chi phí bản quyền phần mềm và máy trạm)
Thời gian chu kỳ 3 - 6 tháng cho mỗi đợt thử nghiệm Nhanh nhưng giới hạn phương trình bậc thấp Cực nhanh (Mô phỏng tức thời trong vài giây)
Độ linh hoạt thông số Kém (Phải thay thế cơ cấu lò xo, thanh xoắn, nhíp thật) Kém (Khó khảo sát nhiều kịch bản phi tuyến) Rất cao (Thay đổi tham số $M_s, C, K$ trực tiếp)
Độ chính xác Tuyệt đối theo điều kiện môi trường đo Thấp do giản lược quá nhiều biến số Cao ($>92%$ tương thích với dữ liệu thực tế)

Yêu cầu kỹ thuật theo thang đo ưu tiên MoSCoW:

  • Must have (Bắt buộc): Hệ phương trình vi phân mô tả đầy đủ 4 tọa độ suy rộng ($E_1, E_2, Z, \varphi$); bộ giải ODE bậc 4 tích hợp tính $a_{RMS}$; mô hình mặt đường chuẩn gờ giảm tốc và sóng sin.
  • Should have (Nên có): Bộ giao diện tương thích đồng bộ dữ liệu hai chiều giữa MATLAB Workspace và CarSim S-Function.
  • Could have (Có thể có): Mở rộng khảo sát đặc tính phi tuyến của giảm chấn thủy lực hai chiều (nén - trả).
  • Won't have (Chưa thực hiện đợt này): Khảo sát dao động lắc ngang (Roll) và dao động xoắn khung gầm 3D.

Thiết kế hệ thống

Mô hình dao động nửa xe 4 bậc tự do trong mặt phẳng dọc được mô tả qua sơ đồ cấu trúc động lực học:

graph TD
    subgraph Road_Excitation ["Nguồn kích động mặt đường"]
        q1["Kích động bánh trước: q1(t)"]
        q2["Kích động bánh sau: q2(t) = q1(t - L/v)"]
    end

    subgraph Unsprung_Masses ["Khối lượng không được treo (KLKĐT)"]
        M1["Cầu trước M1 (Chuyển vị E1)"]
        M2["Cầu sau M2 (Chuyển vị E2)"]
    end

    subgraph Suspension_System ["Hệ thống treo & Lốp"]
        Tire1["Lốp trước: CL1, KL1"]
        Tire2["Lốp sau: CL2, KL2"]
        Susp1["Treo trước: C1, K1"]
        Susp2["Treo sau: C2, K2"]
    end

    subgraph Sprung_Mass ["Khối lượng được treo (KLĐT)"]
        Body["Thân xe (Ms, Jy) - Chuyển vị đứng Z & Góc lắc dọc phi"]
    end

    q1 --> Tire1 --> M1 --> Susp1 --> Body
    q2 --> Tire2 --> M2 --> Susp2 --> Body

Hệ thống phương trình vi phân dao động viết dưới dạng ma trận trạng thái:

$$[M]{\ddot{x}} + [K]{\dot{x}} + [C]{x} = {F(t)}$$

Trong đó:

  • Vector tọa độ suy rộng: ${x} = \begin{bmatrix} E_1 & E_2 & Z & \varphi \end{bmatrix}^T$
  • Ma trận khối lượng: $$[M] = \begin{bmatrix} M_1 & 0 & 0 & 0 \ 0 & M_2 & 0 & 0 \ 0 & 0 & M_s & 0 \ 0 & 0 & 0 & J_y \end{bmatrix}$$
  • Ma trận cản giảm chấn: $$[K] = \begin{bmatrix} K_1+K_{L1} & 0 & -K_1 & aK_1 \ 0 & K_2+K_{L2} & -K_2 & -bK_2 \ -K_1 & -K_2 & K_1+K_2 & -(aK_1-bK_2) \ aK_1 & -bK_2 & -(aK_1-bK_2) & a^2K_1+b^2K_2 \end{bmatrix}$$
  • Ma trận độ cứng: $$[C] = \begin{bmatrix} C_1+C_{L1} & 0 & -C_1 & aC_1 \ 0 & C_2+C_{L2} & -C_2 & -bC_2 \ -C_1 & -C_2 & C_1+C_2 & -(aC_1-bC_2) \ aC_1 & -bC_2 & -(aC_1-bC_2) & a^2C_1+b^2C_2 \end{bmatrix}$$
  • Vector lực kích động từ mặt đường: $${F(t)} = \begin{bmatrix} C_{L1}q_1 + K_{L1}\dot{q}1 \ C{L2}q_2 + K_{L2}\dot{q}_2 \ 0 \ 0 \end{bmatrix}$$
+-----------------------------------------------------------------------------------+
| THÔNG SỐ ĐẦU VÀO XE MITSUBISHI TRITON KHẢO SÁT                                   |
+-----------------------------------------------------------------------------------+
| Khối lượng được treo (Ms): 1850 kg          | Mômen quán tính dọc (Jy): 2850 kg.m2|
| Khối lượng không treo trước (M1): 115 kg   | Khối lượng không treo sau (M2): 140 kg|
| Độ cứng treo trước (C1): 38500 N/m          | Độ cứng treo sau (C2): 52000 N/m    |
| Hệ số cản treo trước (K1): 2400 N.s/m       | Hệ số cản treo sau (K2): 3100 N.s/m |
| Độ cứng lốp trước/sau (CL1=CL2): 210000 N/m| Chiều dài cơ sở (L = a + b): 3.00 m  |
| Khoảng cách trọng tâm - cầu trước (a): 1.38m| Trọng tâm - cầu sau (b): 1.62 m     |
+-----------------------------------------------------------------------------------+

Công nghệ và phiên bản phần mềm sử dụng:

  • MathWorks MATLAB R2023b (Xử lý ma trận, thuật toán giải tích và tính toán thống kê).
  • MathWorks Simulink v10.8 (Mô hình hóa sơ đồ khối dynamic continuous-time).
  • Mechanical Simulation CarSim 2020.0 (Mô phỏng động lực học đa vật thể 3D).

Methodology

Nghiên cứu áp dụng quy trình tính toán kỹ thuật tuần tự theo chu trình V-Model:

[Xác định thông số hình học/động lực học]
                  │
                  ▼
[Xây dựng hệ phương trình vi phân D'Alembert]
                  │
                  ▼
[Lập trình mã nguồn MATLAB & Mô hình Simulink] <───> [Thiết lập mô hình CarSim]
                  │                                           │
                  └───────────────┬───────────────────────────┘
                                  ▼
                     [So sánh & Kiểm nghiệm kết quả]
                                  │
                                  ▼
                 [Đánh giá theo TCVN 6964-1:2001]
  • Đánh giá rủi ro kỹ thuật: Vòng lặp đại số (Algebraic Loop) phát sinh trong Simulink khi giải hệ vi phân bậc cao được triệt tiêu bằng cách sử dụng các khối Unit Delay (độ trễ lấy mẫu $T_s = 0.001\text{ s}$) và thuật toán giải số ode45 (Dormand-Prince) với bước tích phân biến đổi (Variable-step).

Implementation và kết quả

Development process

Quá trình triển khai được cấu trúc thành 4 giai đoạn rõ ràng:

  • Phase 1: Thu thập dữ liệu hình học, thông số khối lượng và độ cứng thực tế của xe Mitsubishi Triton tại xưởng dịch vụ ủy quyền Mitsubishi Nam Auto.
  • Phase 2: Xây dựng script giải bài toán ma trận trên MATLAB để tính toán tần số dao động riêng và dạng dao động (Eigenvalues & Eigenvectors).
  • Phase 3: Xây dựng sơ đồ khối Simulink cho từng cụm cầu trước, cầu sau và thân xe; cấu hình các hàm kích động $q_1(t)$ và $q_2(t) = q_1(t - \tau)$ với thời gian trễ $\tau = L/v$.
  • Phase 4: Thiết lập địa hình thử nghiệm trên CarSim (3D Road Surface) tương ứng với hàm kích động điều hòa và gờ đơn hình sin bán bình phương.

Đoạn mã nguồn MATLAB giải hệ phương trình vi phân và tính toán gia tốc hiệu dụng RMS:

% Script tinh toan dao dong va chi so RMS do em diu
clc; clear; close all;

% 1. Khoi tao thong so dong luc hoc Mitsubishi Triton
Ms = 1850; Jy = 2850; M1 = 115; M2 = 140;
C1 = 38500; C2 = 52000; K1 = 2400; K2 = 3100;
CL1 = 210000; CL2 = 210000; KL1 = 0; KL2 = 0;
a = 1.38; b = 1.62; L = a + b; v_kmh = 40; v = v_kmh / 3.6;

% 2. Thiet lap ma tran M, K, C
M = diag([M1, M2, Ms, Jy]);
K = [K1+KL1, 0, -K1, a*K1;
     0, K2+KL2, -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 = [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];

% 3. Tinh toan tan so dao dong rieng (Natural Frequencies)
[V, D] = eig(inv(M)*C);
omega_n = sqrt(diag(D));
f_n = omega_n / (2 * pi); % f1 (nhun than xe), f2 (lac doc), f3, f4 (cau xe)

% 4. Mo phong phan ung voi kich dong mat duong ban binh phuong hinh sin
q0 = 0.08; S = 1.2; % Chieu cao go 80mm, chieu dai 1.2m
t_bump = S / v; t_sim = 3.0; dt = 0.001;
t = 0:dt:t_sim;
tau = L / v;

q1 = zeros(size(t)); q2 = zeros(size(t));
for i = 1:length(t)
    if t(i) <= t_bump
        q1(i) = q0 * (sin(pi * v * t(i) / S))^2;
    end
    if t(i) >= tau && t(i) <= (tau + t_bump)
        q2(i) = q0 * (sin(pi * v * (t(i) - tau) / S))^2;
    end
end

% 5. Tinh toan gia tri RMS cua gia toc than xe theo TCVN 6964-1:2001
% Gia su vector gia toc tinh tien Z_ddot da duoc giai tu ODE45/Simulink
% a_RMS = sqrt( (1/T) * trapz(t, Z_ddot.^2) )

Testing và validation

Nghiên cứu tiến hành mô phỏng và kiểm nghiệm trên 2 kịch bản mấp mô mặt đường tiêu chuẩn:

  1. Kịch bản 1 - Kích động điều hòa: Biên độ mấp mô $q_0 = 0.05\text{ m}$, bước sóng mấp mô $S = 5.0\text{ m}$, vận tốc xe thay đổi từ $20 - 80\text{ km/h}$.
  2. Kịch bản 2 - Kích động gờ giảm tốc (Bán bình phương hình sin): Chiều cao gờ $q_0 = 0.08\text{ m}$, chiều dài $S = 1.2\text{ m}$, vận tốc xe $v = 30\text{ km/h}$ và $v = 50\text{ km/h}$.
+-----------------------------------------------------------------------------------+
| KẾT QUẢ SO SÁNH BIÊN ĐỘ GIA TỐC VÀ DỊCH CHUYỂN (TẠI VẬN TỐC 40 KM/H)              |
+------------------------------------+-----------------+---------------+------------+
| Thông số khảo sát                  | MATLAB/Simulink | Phần mềm CarSim| Sai số (%) |
+------------------------------------+-----------------+---------------+------------+
| Dịch chuyển đứng cực đại Zmax (m)  | 0.0462          | 0.0485        | 4.74%      |
| Gia tốc đứng cực đại Zddot_max (m/s2)| 2.842         | 2.975         | 4.47%      |
| Góc lắc dọc cực đại phi_max (rad)  | 0.0315          | 0.0331        | 4.83%      |
| Gia tốc góc lắc dọc max (rad/s2)   | 1.954           | 2.062         | 5.24%      |
| Trọng số gia tốc a_RMS (m/s2)      | 0.482           | 0.511         | 5.67%      |
+------------------------------------+-----------------+---------------+------------+

+-----------------------------------------------------------------------------------+
| BIỂU ĐỒ ĐÁP ỨNG GIA TỐC THÂN XE THEO THỜI GIAN KHI QUA GỜ GIẢM TỐC               |
|                                                                                   |
|  Gia tốc Z'' (m/s2)                                                               |
|   3.0 ┼               /---\ (CarSim)                                              |
|   2.0 ┼              /     \                                                      |
|   1.0 ┼    /---\    /       \   /---\ (MATLAB)                                    |
|   0.0 ┼───/─────\──/─────────\─/─────\───────────────> Thời gian t (s)            |
|  -1.0 ┼  /       \/           V       \                                           |
|  -2.0 ┼ /                              \                                          |
|  -3.0 ┼/                                \                                         |
|       0.0      0.5      1.0     1.5     2.0     2.5                               |
+-----------------------------------------------------------------------------------+

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

  1. Về mặt tần số riêng ($f_0$):

    • Tần số dao động nhún thân xe: $f_1 = 1.18\text{ Hz}$ (Tương ứng $70.8\text{ lần/phút}$).
    • Tần số dao động lắc dọc thân xe: $f_2 = 1.35\text{ Hz}$ (Tương ứng $81.0\text{ lần/phút}$).
    • Kết quả này hoàn toàn nằm trong dải tần số êm dịu tối ưu cho xe chở người và xe bán tải ($60 - 85\text{ lần/phút}$ tương đương $1.0 - 1.4\text{ Hz}$).

  2. Về chỉ tiêu gia tốc hiệu dụng RMS ($a_{RMS}$):

    • Trên đường mấp mô điều hòa ở dải tốc độ $30 - 50\text{ km/h}$, giá trị $a_{RMS}$ dao động từ $0.38 - 0.49\text{ m/s}^2$.
    • Đối chiếu với Bảng 2.1 (TCVN 6964-1:2001), mức gia tốc này thuộc dải $0.315 - 0.63\text{ m/s}^2$, tạo cảm giác "Chỉ có một chút ít không thoải mái", chứng minh kết cấu hệ thống treo nguyên bản của Mitsubishi Triton đảm bảo tốt điều kiện tiện nghi cho hành khách.

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

Nghiên cứu mang lại các đóng góp khoa học và kỹ thuật thực tiễn:

  • Phương pháp luận tích hợp: Thay vì chỉ giải tích thuần túy hoặc mô phỏng đơn lẻ, đồ án xây dựng quy trình kết nối đối sánh song song giữa mô hình trạng thái không gian ma trận trên MATLAB, mô hình luồng tín hiệu Simulink và mô hình động lực học đa vật thể CarSim.
  • Tính chính xác và hiệu quả thời gian: Mô hình toán 4-DOF đạt độ chính xác cao với sai số dưới $6%$ so với phần mềm CAE công nghiệp CarSim, nhưng giảm hơn $85%$ thời gian thiết lập và tính toán.
  • Bằng chứng số liệu thực tế: Toàn bộ dữ liệu đầu vào ($M_s, M_1, M_2, C_1, C_2, a, b$) được thu thập và đo đạc trực tiếp trên xe thật tại xưởng dịch vụ chính hãng Mitsubishi Nam Auto, đảm bảo tính thực chứng cao thay vì sử dụng số liệu lý thuyết giả định.

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

Tình huống ứng dụng thực tế

+-----------------------------------------------------------------------------------+
| CÁC TÌNH HUỐNG ỨNG DỤNG THỰC TẾ CỦA MÔ HÌNH                                       |
+-----------------------------------------------------------------------------------+
| 1. Tối ưu hóa hệ số cản K1, K2 của giảm chấn độ (Aftermarket Shock Absorbers)     |
| 2. Cân chỉnh độ cứng nhíp lá C2 khi hoán cải xe bán tải thành xe cứu thương/chở tiền|
| 3. Tích hợp thuật toán điều khiển cho hệ thống treo bán chủ động (Semi-Active)   |
| 4. Xây dựng bài thí nghiệm số cho sinh viên ngành Kỹ thuật Ô tô                   |
+-----------------------------------------------------------------------------------+

Yêu cầu triển khai hệ thống (System Requirements)

  • Phần cứng (Hardware):
    • CPU: Intel Core i5 / AMD Ryzen 5 thế hệ 8 trở lên (Tối thiểu 4 cores, 2.5 GHz).
    • RAM: Tối thiểu 8 GB (Khuyến nghị 16 GB).
    • Ổ cứng: Trống tối thiểu 10 GB SSD.
  • Môi trường phần mềm (Software Stack):
    • Hệ điều hành: Windows 10/11 64-bit hoặc Linux Ubuntu 20.04/22.04 LTS.
    • MATLAB & Simulink: Phiên bản R2018b trở lên (Toolbox: Control System Toolbox, Signal Processing Toolbox).
    • CarSim: Phiên bản 2019.0 trở lên.

Hướng dẫn cài đặt và chạy mô hình

  1. Khởi động MATLAB, điều hướng thư mục làm việc đến thư mục chứa mã nguồn: 
    cd('/path/to/Triton_Ride_Comfort_Project');
  2. Nạp file thông số xe vào không gian làm việc (Workspace): 
    run('parameters_triton.m');
  3. Khởi chạy mô phỏng Simulink: 
    sim('Half_Car_4DOF_Model.slx');
  4. Xuất đồ thị và tính toán chỉ số RMS: 
    run('plot_results_and_calc_rms.m');

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

Hạn chế kỹ thuật hiện tại

  • Mô hình phẳng 4 bậc tự do chỉ khảo sát dao động trong mặt phẳng dọc đối xứng, chưa đánh giá được dao động lắc ngang (Roll) khi xe quay vòng hoặc đi vào mặt đường có độ gồ ghề lệch pha giữa vệt bánh trái và vệt bánh phải.
  • Xem đặc tính lốp và giảm chấn là tuyến tính thuần túy, chưa xét đến tính phi tuyến của van giảm chấn thủy lực và tính trễ đàn hồi (Hysteresis) của lốp xe cao su.

Hướng nghiên cứu mở rộng

  • Nâng cấp mô hình lên Mô hình không gian 7 bậc tự do (7-DOF Full-Car Model) bao gồm 4 chuyển vị bánh xe, 1 chuyển vị nhún thân xe, 1 góc lắc dọc (Pitch) và 1 góc lắc ngang (Roll).
  • Nghiên cứu ứng dụng giải thuật điều khiển thông minh (Fuzzy Logic Controller / LQR - Linear Quadratic Regulator) trên hệ thống treo bán chủ động sử dụng giảm chấn lưu chất từ tính (Magnetorheological Damper - MR Damper) nhằm giảm thêm $25 - 35%$ giá trị gia tốc RMS thân xe.

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

  • Sinh viên & Học viên cao học Kỹ thuật Ô tô: Nắm vững phương pháp lập phương trình động lực học vi phân bằng nguyên lý D'Alembert và kỹ năng mô phỏng chuyên sâu trên MATLAB/Simulink và CarSim.
  • Kỹ sư R&D và Thiết kế ô tô: Sở hữu công cụ tính toán nhanh để đánh giá sơ bộ chất lượng êm dịu trước khi sản xuất thử nghiệm hoặc hoán cải hệ thống treo.
  • Các xưởng nâng cấp (Tuning/Aftermarket Workshops): Có cơ sở lý thuyết định lượng để tư vấn lựa chọn bộ phuộc giảm xóc và nhíp độ phù hợp cho từng mục đích sử dụng (off-road, chở hàng nặng, chở gia đình).
  • Nhà nghiên cứu học thuật: Cung cấp bộ tham số thực nghiệm chi tiết và mã nguồn mở chuẩn xác phục vụ các nghiên cứu mở rộng về hệ thống điều khiển khung gầm chủ động.

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ô phỏng?

Hệ thống yêu cầu máy tính chạy Windows 10/11 hoặc Linux, RAM tối thiểu 8GB, cài đặt MATLAB R2018b trở lên kèm Simulink. Nếu muốn kiểm chứng trực quan 3D đa vật thể, cần cài đặt thêm phần mềm CarSim 2019 trở lên.

2. Mô hình nửa xe 4-DOF có khả năng mở rộng lên mô hình toàn xe 7-DOF không?

Hoàn toàn khả thi. Bằng cách bổ sung thêm 3 phương trình vi phân cho 2 bánh xe phía đối diện và góc lắc ngang quanh trục X (Roll), ma trận $[M], [K], [C]$ sẽ được mở rộng từ kích thước $4 \times 4$ lên $7 \times 7$ mà không làm thay đổi cấu trúc thuật toán giải.

3. Làm thế nào để tích hợp mô hình này với dữ liệu đo đạc thực tế từ cảm biến?

Dữ liệu gia tốc từ cảm biến IMU (Inertial Measurement Unit) gắn trên xe thật có thể được lưu dưới dạng file .csv hoặc .mat, sau đó đưa vào MATLAB để lọc tín hiệu (Butterworth Filter) và tính toán trực tiếp $a_{RMS}$ để đối soát với kết quả mô phỏng.

4. Tại sao giá trị gia tốc RMS lại quan trọng hơn giá trị gia tốc cực đại?

Gia tốc cực đại ($Z_{max}$) chỉ phản ánh biên độ tức thời tại một điểm va chạm, trong khi gia tốc hiệu dụng RMS ($a_{RMS}$) tích phân toàn bộ năng lượng rung động tác động lên cơ thể người trong suốt một khoảng thời gian, phản ánh chính xác mức độ mệt mỏi và tổn hại sức khỏe theo ISO 2631.

5. Chi phí và thời gian triển khai mô hình mô phỏng này so với thử nghiệm thực địa?

Triển khai mô hình số giúp tiết kiệm hơn $90%$ chi phí so với việc thuê đường thử chuyên dụng và gắn hàng chục cảm biến đắt tiền, đồng thời rút ngắn chu kỳ đánh giá từ hàng tuần xuống chỉ còn vài giây xử lý tính toán.

Kết luận

Đồ án "Khảo sát độ êm dịu trên xe ô tô" đã hoàn thành xuất sắc các mục tiêu nghiên cứu đề ra: Thiết lập thành công mô hình toán học nửa xe 4 bậc tự do, lập trình giải bài toán dao động trên MATLAB/Simulink và kiểm nghiệm tương thích với phần mềm công nghiệp CarSim với sai số dưới $6%$. Kết quả khảo sát trên xe Mitsubishi Triton khẳng định tần số dao động riêng ($1.18\text{ Hz}$ và $1.35\text{ Hz}$) cùng chỉ số gia tốc hiệu dụng $a_{RMS}$ ($0.38 - 0.49\text{ m/s}^2$) hoàn toàn đáp ứng tiêu chuẩn êm dịu TCVN 6964-1:2001 và ISO 2631-1:1997. Công trình mở ra hướng ứng dụng thực tiễn mạnh mẽ trong công tác nghiên cứu, giảng dạy và tối ưu hóa hệ thống treo trên các phương tiện giao thông vận tải hiện đại.