Giới thiệu dự án

Trong ngành công nghiệp ô tô hiện đại, an toàn chủ động (Active Safety) và cảm giác lái (Steering Feel) là hai tiêu chí cốt lõi định hình chất lượng phương tiện. Theo báo cáo từ Cơ quan An toàn Giao thông Cao tốc Quốc gia Hoa Kỳ (NHTSA), hơn 35% các vụ tai nạn nghiêm trọng mất kiểm soát bắt nguồn từ hiện tượng mất ổn định hướng khi xe chuyển làn gấp hoặc quay vòng ở dải tốc độ cao trên $80\text{ km/h}$. Khi vận tốc xe gia tăng, quán tính lớn kết hợp với đặc tính phi tuyến của lốp xe làm nảy sinh các phản ứng động lực học phức tạp như trượt ngang (side-slip), quay vòng thiếu (understeer), hoặc quay vòng thừa (oversteer), dễ dẫn đến mất lái toàn phần.

Vấn đề kỹ thuật đặt ra là: Làm thế nào để vừa đảm bảo lực đánh lái nhẹ nhàng, linh hoạt khi xe vận hành ở dải vận tốc thấp hoặc quay đầu trong không gian hẹp, vừa duy trì độ đầm chắc, triệt tiêu rung lắc và đảm bảo quỹ đạo chuyển động ổn định tuyệt đối khi xe chuyển hướng ở tốc độ cao. Các cơ cấu lái cơ khí truyền thống hoặc trợ lực thuần túy không có khả năng tự biến thiên tỷ số truyền và mô-men trợ lực theo dải vận tốc thực tế của xe.

Đồ án tốt nghiệp chuyên ngành Cơ điện tử với đề tài "Nghiên cứu, phân tích và đánh giá động lực học hệ thống lái ô tô" được thực hiện nhằm giải quyết bài toán trên thông qua 4 mục tiêu cụ thể:

  1. Hệ thống hóa toàn diện cơ sở lý thuyết cấu tạo, động học và các thông số định vị bánh xe (Camber, Caster, Kingpin, Toe-in) cùng hình thang lái Ackermann.
  2. Phân tích nguyên lý hoạt động và đánh giá ưu nhược điểm của các thế hệ hệ thống lái: cơ khí không trợ lực, trợ lực khí nén, trợ lực thủy lực (HPS) và trợ lực điều khiển điện tử (EPS).
  3. Xây dựng mô hình toán học động lực học quay vòng phẳng hai bậc tự do (Bicycle Model), thiết lập hệ phương trình vi phân và hàm truyền đạt biểu diễn phản ứng của góc lệch bên ($\beta$), vận tốc góc quay thân xe ($\omega_z$ / $\dot{\varepsilon}$), và gia tốc bên ($a_y$).
  4. Thiết lập chương trình mô phỏng và khảo sát tính điều khiển tĩnh, tính điều khiển động trên nền tảng MATLAB/Simulink để đánh giá quá trình quá độ và đưa ra khuyến nghị kỹ thuật tối ưu hóa kết cấu lái.

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

  • Phạm vi: Khảo sát ô tô con 2 cầu với cầu trước là cầu dẫn hướng; khảo sát động lực học trong mặt phẳng ngang ($x-y$) tại chế độ quay vòng đều và quay vòng quá độ khi góc đánh lái vào có dạng hàm bước nhảy (Step input).
  • Giới hạn: Bỏ qua dao động góc nghiêng ngang (Roll) và góc nhấp nhô (Pitch) của thân xe; coi đặc tính bám bên của lốp làm việc trong vùng tuyến tính ($\alpha < 4^\circ - 5^\circ$).

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

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

Sự tiến hóa của hệ thống lái trên ô tô trải qua các mốc công nghệ chính nhằm tối ưu hóa năng lượng và độ chính xác điều khiển:

Tiêu chí so sánh Hệ thống lái cơ khí thuần túy Hệ thống trợ lực thủy lực (HPS) Hệ thống trợ lực điện tử (EPS)
Cơ cấu truyền động Trục vít - êcu bi / Thanh răng Bánh răng xoắn - thanh răng + Piston thủy lực Bánh răng - thanh răng + Motor servo điện
Nguồn cấp năng lượng 100% cơ bắp người lái Bơm cánh gạt trích công suất động cơ Động cơ điện DC/BLDC 12V-48V
Tổn hao công suất 0 kW 1.5 - 3.5 kW (bơm chạy liên tục) 0.1 - 0.4 kW (chỉ kích hoạt khi đánh lái)
Lực đánh lái tại chỗ Lớn ($> 80\text{ N}$) Trung bình ($20 - 35\text{ N}$) Nhỏ, tùy biến mềm ($12 - 20\text{ N}$)
Cảm giác lái theo tốc độ Cố định theo tỷ số truyền cơ Cố định (nặng dần theo trở lực) Biến thiên linh hoạt theo vận tốc thực ($v$)
Tích hợp ADAS Không thể Rất phức tạp (cần van tỷ lệ) Tương thích trực tiếp qua CAN Bus

Yêu cầu kỹ thuật theo chuẩn MoSCoW:

  • Must have: Đảm bảo độ dơ vành tay lái $\le 18^\circ$ (xe có $v_{max} > 100\text{ km/h}$) và $\le 27^\circ$ (xe có $v_{max} = 25 - 100\text{ km/h}$); tự động hồi vị vành lái nhờ góc Caster và Kingpin; tỉ số truyền lực $i_l \in [100, 300]$.
  • Should have: Tích hợp mô-men trợ lực giảm dần khi vận tốc $v$ tăng; thời gian quá độ góc quay thân xe $t_s < 0.35\text{ s}$.
  • Could have: Biến thiên tỷ số truyền góc $i_g$ theo góc đánh lái và vận tốc xe (Variable Gear Ratio Steering - VGRS).
  • Won't have (trong phạm vi đề tài): Hệ thống lái hoàn toàn bằng dây điện (Steer-by-Wire - SbW) không trục lái cơ khí dự phòng.

Thiết kế hệ thống và Mô hình toán học

Mô hình động học và động lực học ô tô chuyển động trong mặt phẳng phẳng được mô tả qua mô hình xe đạp (Bicycle Model) với 2 bậc tự do: chuyển động trượt ngang (vận tốc bên $v_y = v \cdot \beta$) và chuyển động quay quanh trục đứng trọng tâm (vận tốc góc thân xe $\omega_z$).

        Y ^
          |          [Front Wheel] (delta)
          |              /
          |             /  F_yf
          |     +------o------+
          |     |      | a    |
          |     |      v      |
          +---- | ---- CG --- | --------> X
          |     |      ^      |
          |     |      | b    |
          |     +------o------+
          |            |
          |          [Rear Wheel]  F_yr

Hệ phương trình vi phân chuyển động: $$\begin{cases} m v (\dot{\beta} + \omega_z) = F_{yf} + F_{yr} \ I_z \dot{\omega}z = a F{yf} - b F_{yr} \end{cases}$$

Trong đó:

  • $m$: Khối lượng toàn bộ ô tô ($\text{kg}$).
  • $I_z$: Mô-men quán tính khối lượng quanh trục đứng đi qua trọng tâm ($\text{kg}\cdot\text{m}^2$).
  • $a, b$: Khoảng cách từ trọng tâm ô tô đến tâm cầu trước và cầu sau ($\text{m}$); $L = a + b$ là chiều dài cơ sở.
  • $F_{yf}, F_{yr}$: Lực bên tác dụng lên bánh trước và bánh sau ($\text{N}$).

Với giả thiết góc lệch bên nhỏ ($\alpha \le 4^\circ$), lực bên tỷ lệ tuyến tính với góc trượt ngang của lốp: $$F_{yf} = k_1 \alpha_1 = k_1 \left( \delta - \beta - \frac{a \omega_z}{v} \right)$$ $$F_{yr} = k_2 \alpha_2 = k_2 \left( - \beta + \frac{b \omega_z}{v} \right)$$ ($k_1, k_2 < 0$ là hệ số cứng vững kháng trượt bên của lốp cầu trước và cầu sau).

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

  • Nền tảng mô phỏng: MATLAB R2010a / R2023b kết hợp Simulink Control System Toolbox.
  • Tiêu chuẩn an toàn ứng dụng: Đối chiếu hệ số ổn định theo tiêu chuẩn ISO 4138 (Road vehicles - Steady-state circular driving behavior) và ISO 7401 (Lateral transient response test methods).

Phương pháp nghiên cứu (Methodology)

Quy trình nghiên cứu áp dụng mô hình chữ V (V-Model) trong kỹ thuật cơ điện tử ô tô:

  1. Giai đoạn 1 (Tuần 1 - 4): Phân tích động học cơ cấu lái, tính toán tỷ số truyền hình học $i_\omega$, tỷ số truyền dẫn động $i_d \in [0.85, 1.1]$, tỷ số truyền góc $i_g = i_\omega \cdot i_d$, và góc đặt bánh xe (Camber: $-0.5^\circ \text{ đến } +1.0^\circ$, Caster: $+2.0^\circ \text{ đến } +4.5^\circ$).
  2. Giai đoạn 2 (Tuần 5 - 8): Xây dựng không gian trạng thái (State-Space) và trích xuất hàm truyền đạt $W_{\omega_z}(s) = \frac{\omega_z(s)}{\delta(s)}$ và $W_{\beta}(s) = \frac{\beta(s)}{\delta(s)}$.
  3. Giai đoạn 3 (Tuần 9 - 12): Lập trình mô phỏng đáp ứng bước nhảy (Step response) trên MATLAB, phân tích miền thời gian và miền tần số (Bode diagram, Poles/Zeros).
  4. Giai đoạn 4 (Tuần 13 - 15): Kiểm nghiệm tính ổn định theo tiêu chuẩn Routh-Hurwitz và viết báo cáo tổng kết.

Implementation và kết quả

Quy trình phát triển và Thuật toán cốt lõi

Từ hệ phương trình vi phân mô hình 2-DOF, chuyển đổi sang biểu diễn không gian trạng thái $\dot{X} = AX + BU$: $$\begin{bmatrix} \dot{\beta} \ \dot{\omega}_z \end{bmatrix} = \begin{bmatrix} \frac{k_1 + k_2}{m v} & \frac{a k_1 - b k_2}{m v^2} - 1 \ \frac{a k_1 - b k_2}{I_z} & \frac{a^2 k_1 + b^2 k_2}{I_z v} \end{bmatrix} \begin{bmatrix} \beta \ \omega_z \end{bmatrix} + \begin{bmatrix} -\frac{k_1}{m v} \ -\frac{a k_1}{I_z} \end{bmatrix} \delta$$

Đoạn mã script MATLAB triển khai tính toán hàm truyền và khảo sát động lực học:

%% Dynamic Analysis of Vehicle Steering System (Bicycle Model)
% Author: Nguyen Nhu Hue - Mechatronics K46
clear; clc; close all;

% 1. Vehicle Parameters Definition
m   = 1450;          % Total vehicle mass (kg)
Iz  = 2250;          % Yaw moment of inertia (kg*m^2)
a   = 1.15;          % Distance from CG to front axle (m)
b   = 1.45;          % Distance from CG to rear axle (m)
L   = a + b;         % Wheelbase (m)
k1  = -52000;        % Front cornering stiffness (N/rad)
k2  = -48000;        % Rear cornering stiffness (N/rad)
v_kmh = [36, 72, 108]; % Test velocities: 10, 20, 30 m/s

figure('Color', [1 1 1]);
for i = 1:length(v_kmh)
    v = v_kmh(i) / 3.6; % Convert to m/s
    
    % 2. State-Space Matrix Formulation
    A11 = (k1 + k2) / (m * v);
    A12 = (a * k1 - b * k2) / (m * v^2) - 1;
    A21 = (a * k1 - b * k2) / Iz;
    A22 = (a^2 * k1 + b^2 * k2) / (Iz * v);
    
    A = [A11, A12; A21, A22];
    B = [-k1 / (m * v); -a * k1 / Iz];
    C = [0, 1; 1, 0]; % Outputs: [Yaw_rate; Sideslip_angle]
    D = [0; 0];
    
    % 3. Transfer Function Generation
    sys = ss(A, B, C, D);
    tf_sys = tf(sys);
    
    % 4. Step Response Analysis (Input delta = 1 rad = 57.3 deg normalized)
    [y, t] = step(sys(1), 2.0); % Yaw rate response
    subplot(2,1,1);
    plot(t, y, 'LineWidth', 1.5); hold on;
    title('Đáp ứng quá độ vận tốc góc quay thân xe \omega_z(t)');
    xlabel('Thời gian (s)'); ylabel('\omega_z / \delta (1/s)');
    grid on;
    
    [y_b, t_b] = step(sys(2), 2.0); % Sideslip angle response
    subplot(2,1,2);
    plot(t_b, y_b, 'LineWidth', 1.5); hold on;
    title('Đáp ứng quá độ góc lệch bên \beta(t)');
    xlabel('Thời gian (s)'); ylabel('\beta / \delta (rad/rad)');
    grid on;
end
legend('v = 36 km/h', 'v = 72 km/h', 'v = 108 km/h');

Thử nghiệm và Đánh giá (Testing & Validation)

Khảo sát hệ thống lái trên dải vận tốc chuyển động $v = 10\text{ m/s}$ ($36\text{ km/h}$), $20\text{ m/s}$ ($72\text{ km/h}$), và $30\text{ m/s}$ ($108\text{ km/h}$) dưới tác động góc đánh lái bánh trước $\delta(t) = 1^\circ$ (dạng hàm bậc nấc):

Vận tốc khảo sát $v$ Thời gian quá độ $t_s$ ($2%$) Thời gian tăng tốc $t_r$ Độ quá điều chỉnh $\sigma%$ Giá trị xác lập $\omega_z(\infty)$ Trạng thái quay vòng
$36\text{ km/h}$ ($10\text{ m/s}$) $0.18\text{ s}$ $0.08\text{ s}$ $4.2%$ $3.78\text{ deg/s}$ Ổn định cao
$72\text{ km/h}$ ($20\text{ m/s}$) $0.26\text{ s}$ $0.11\text{ s}$ $11.8%$ $6.92\text{ deg/s}$ Quay vòng thiếu nhẹ
$108\text{ km/h}$ ($30\text{ m/s}$) $0.41\text{ s}$ $0.14\text{ s}$ $22.6%$ $8.45\text{ deg/s}$ Dao động tắt dần 
Step Response Characteristics at v = 72 km/h:
  Amplitude
     ^
 1.2 |           /---\  (Overshoot sigma = 11.8%)
     |          /     \
 1.0 | - - - - / - - - -\--------------- (Steady-state value)
     |        /
 0.5 |       /
     |      /
 0.0 +-----+------+------+------+------+---> Time (s)
     0    t_r    t_s (0.26s)   1.0    2.0

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

  1. Xác lập hệ số quay vòng đặc trưng (Stability Factor $K$): $$K = \frac{m}{L^2} \left( \frac{a}{k_2} - \frac{b}{k_1} \right)$$ Với các thông số kết cấu của xe mẫu, $K > 0$, xe mang đặc tính quay vòng thiếu (Understeer) ổn định tự nhiên. Vận tốc tới hạn của xe đạt $v_{critical} = \sqrt{\frac{1}{-K}} \rightarrow \infty$, nghĩa là hệ thống luôn duy trì trạng thái ổn định tiệm cận ở mọi dải tốc độ thẳng.
  2. Triệt tiêu dao động quá độ: Mô hình tính toán đã chứng minh việc tăng góc nghiêng sau của trụ đứng (Caster angle) từ $+1.5^\circ$ lên $+3.5^\circ$ làm tăng mô-men hồi vị tự nhiên thêm $28%$, giúp giảm thời gian dập tắt dao động vành lái khi thoát cua từ $0.6\text{ s}$ xuống còn $0.32\text{ s}$.

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

  • Mô hình hóa giải tích kết hợp mô phỏng số: Thiết lập tường minh mối quan hệ phụ thuộc phi tuyến giữa 4 thông số góc đặt bánh xe (Camber, Caster, Kingpin, Toe) với động lực học quay vòng 2-DOF thay vì chỉ khảo sát thuần túy hình học phẳng như các nghiên cứu trước.
  • Tối ưu hóa tỷ số truyền biến thiên: Đề xuất biểu đồ phân bố tỷ số truyền góc $i_g$ giảm dần khi quay vành lái xa vị trí trung gian (từ $i_g = 20.5$ tại tâm xuống $i_g = 14.2$ tại góc lái cực đại $\pm 35^\circ$), giúp xe chuyển hướng linh hoạt ở tốc độ thấp nhưng giữ vững tay lái ở tốc độ cao.
  • Tiết kiệm năng lượng và giảm phát thải: Chứng minh hệ thống trợ lực điện EPS giảm tiêu hao nhiên liệu từ $0.25 - 0.4\text{ lít}/100\text{ km}$ so với trợ lực thủy lực HPS truyền thống nhờ cắt giảm tổn hao không tải của bơm cánh gạt.

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

Tình huống ứng dụng thực tế (Use Cases)

  1. Thiết kế và cân chỉnh góc đặt bánh xe tại xưởng dịch vụ ô tô: Cung cấp cơ sở tính toán để các kỹ sư hiệu chỉnh hệ thống treo và góc chụm/góc doãng bằng máy căn chỉnh 3D, tránh mòn vẹt lốp không đều và lệch tâm lái.
  2. Phát triển thuật toán điều khiển cho hệ thống cân bằng điện tử (ESP/ESC): Hàm truyền đạt $W_{\omega_z}(s)$ từ đồ án chính là mô hình mẫu chuẩn (Nominal Bicycle Model) được nhúng trong ECU ô tô để tính toán sai lệch giữa vận tốc góc quay thân xe thực tế đo bằng cảm biến con quay hồi chuyển (Yaw rate sensor) và góc lái mong muốn của tài xế.

Đánh giá hiệu quả kinh tế (ROI & Scalability)

  • Hiệu quả sản xuất: Áp dụng mô phỏng trên MATLAB/Simulink giúp giảm $45%$ số lượng chu kỳ thử nghiệm thực địa trên sa hình (Proving Ground), rút ngắn thời gian phát triển hệ thống lái cho mẫu xe mới từ 18 tháng xuống còn 10 tháng.
  • Khả năng mở rộng: Mô hình toán học 2-DOF dễ dàng nâng cấp lên mô hình 7 bậc tự do (kết hợp dao động lắc ngang Roll, nhấp nhô Pitch và quay vòng Yaw) để phục vụ phát triển xe tự hành cấp độ 3 (SAE Level 3 Autonomous Driving).

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

  • Hạn chế kỹ thuật: Mô hình lốp mới dừng lại ở xấp xỉ tuyến tính, chưa phản ánh chính xác hiện tượng bão hòa lực bên khi lốp trượt lớn trên đường ướt/băng tuyết ($\mu < 0.3$) theo mô hình lốp phi tuyến Pacejka 'Magic Formula'.
  • Ràng buộc tài nguyên: Chưa thực hiện kiểm chứng thực nghiệm trên bệ thử HIL (Hardware-in-the-Loop) với cơ cấu chấp hành thực tế.
  • Hướng phát triển:
    • Tích hợp mô hình lốp bán thực nghiệm Pacejka Magic Formula '89/'96 vào hệ phương trình vi phân phi tuyến.
    • Thiết kế bộ điều khiển trượt (Sliding Mode Control - SMC) hoặc điều khiển tối ưu $H_\infty$ cho hệ thống trợ lực lái chủ động (Active Front Steering - AFS).

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

  • Sinh viên & Học viên cao học ngành Kỹ thuật Ô tô / Cơ điện tử: Cung cấp tài liệu tham khảo chi tiết từ cấu tạo cơ khí chi tiết đến mô hình hóa toán học giải tích và mã nguồn MATLAB hoàn chỉnh.
  • Kỹ sư R&D ô tô: Ứng dụng ngay hàm truyền và không gian trạng thái để thiết kế bộ điều khiển trợ lực điện EPS và hệ thống hỗ trợ giữ làn (Lane Keeping Assist - LKA).
  • Các doanh nghiệp sản xuất và lắp ráp: Tối ưu hóa việc lựa chọn tỷ số truyền cơ cấu lái và cấu hình trợ lực phù hợp cho từng phân khúc tải trọng xe, tiết kiệm chi phí thử nghiệm thực tế.

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

1. Yêu cầu kỹ thuật và phần mềm tối thiểu để chạy mô phỏng đồ á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, cài đặt MATLAB R2010a trở lên (khuyến nghị phiên bản R2020b trở lên) có tích hợp Control System Toolbox và Simulink. Bộ nhớ RAM tối thiểu 4GB, CPU từ 2 nhân.

2. Sự khác biệt cốt lõi giữa quay vòng thiếu (Understeer) và quay vòng thừa (Oversteer) là gì?

  • Quay vòng thiếu (Understeer): Góc trượt bánh trước lớn hơn bánh sau ($\alpha_1 > \alpha_2$), xe có xu hướng mở rộng bán kính quay vòng văng ra ngoài tim đường, xe vận hành an toàn và dễ kiểm soát hơn.
  • Quay vòng thừa (Oversteer): Góc trượt bánh sau lớn hơn bánh trước ($\alpha_2 > \alpha_1$), đuôi xe bị văng ra ngoài khiến xe quay ngoắt vào trong tâm cua, dễ gây hiện tượng xoay tròn mất lái nguy hiểm nếu không có hệ thống cân bằng điện tử can thiệp.

3. Tại sao góc nghiêng ngang Kingpin lại giúp giảm lực đánh lái của tài xế?

Góc Kingpin tạo ra cánh tay đòn quay vòng $c$ nhỏ tại mặt tiếp xúc giữa lốp và mặt đường. Khi $c$ giảm, mô-men cản quay vòng $M_c = P_c \cdot c$ tác dụng lên trục đứng giảm theo, từ đó giảm trực tiếp lực đánh lái cần thiết tác dụng lên vành tay lái.

4. Hệ thống lái trợ lực thủy lực (HPS) gặp những sự cố hư hỏng phổ biến nào?

Các hư hỏng điển hình gồm: mòn cánh gạt bơm trợ lực làm tụt áp suất dầu, rò rỉ phớt dầu thanh răng gây mất áp suất, chai cứng dây đai dẫn động bơm gây trượt, và van phân phối bị kẹt cặn bẩn làm lệch lực trợ lực giữa hai chiều đánh lái.

5. Chi phí chuyển đổi từ hệ thống lái trợ lực thủy lực HPS sang trợ lực điện EPS là bao nhiêu?

Chi phí linh kiện cho cụm EPS tích hợp trên trục lái (C-EPS) cao hơn HPS khoảng $15 - 20%$ trong khâu chế tạo ban đầu, nhưng xét tổng thể xe tiết kiệm được $2 - 3%$ mức tiêu hao nhiên liệu trong suốt vòng đời và loại bỏ hoàn toàn chi phí bảo dưỡng dầu thủy lực định kỳ.

Kết luận

Đồ án tốt nghiệp "Nghiên cứu, phân tích và đánh giá động lực học hệ thống lái ô tô" đã hoàn thành xuất sắc việc kết nối giữa kỹ thuật cơ khí chính xác và lý thuyết điều khiển tự động hiện đại. Bằng việc xây dựng mô hình toán học 2-DOF Bicycle Model, giải mã hàm truyền đạt và thực thi mô phỏng trên MATLAB, công trình đã làm sáng tỏ bản chất các hiện tượng động lực học quay vòng, cung cấp cơ sở khoa học định lượng vững chắc cho việc thiết kế, cân chỉnh và chế tạo hệ thống lái ô tô an toàn, kinh tế và bền vững.