Giới thiệu dự án

Sự gia tăng nhanh chóng của phương tiện giao thông cơ giới đường bộ tại Việt Nam đặt ra bài toán cấp thiết về kiểm soát chất lượng an toàn kỹ thuật và bảo vệ môi trường. Theo thống kê của Cục Đăng kiểm Việt Nam, hàng năm có hàng trăm nghìn phương tiện mới được nhập khẩu và lắp ráp trong nước. Để đảm bảo an toàn giao thông, Chính phủ đã ban hành hàng rào kỹ thuật bắt buộc thông qua hệ thống quy chuẩn quốc gia, tiêu biểu là QCVN 09:2015/BGTVT.

Tuy nhiên, nhiều đơn vị kỹ thuật và nghiên cứu gặp khó khăn trong việc thiết lập mô hình tính toán động lực học đồng bộ nhằm đánh giá định lượng tính năng phương tiện trước khi thử nghiệm thực tế. Đồ án tốt nghiệp "Khảo sát và đánh giá tính năng hoạt động của ô tô" (Trường Đại học Sư phạm Kỹ thuật TP.HCM) giải quyết trực diện bài toán này bằng cách xây dựng hệ thống tính toán giải tích toàn diện trên đối tượng khảo sát thực tế: xe bán tải Mitsubishi Triton 2016 2WD (Model KK3T, động cơ 4D56 2.5L Diesel DOHC Turbo Intercooler).

graph LR
    A[Thông số kỹ thuật OEM Mitsubishi Triton] --> B[Mô hình toán học Động lực học ô tô]
    B --> C[Tính toán Sức kéo & Tiêu hao]
    B --> D[Động học Lái & Ổn định thân xe]
    B --> E[Hiệu quả Phanh & An toàn chủ động]
    C --> F{Đối chiếu QCVN 09:2015/BGTVT}
    D --> F
    E --> F
    F --> G[Kết luận Đạt chuẩn & Đề xuất vận hành]

Mục tiêu dự án

  1. Mô hình hóa đường đặc tính ngoài động cơ: Xây dựng chính xác đặc tính mô-men xoắn và công suất động cơ Diesel 4D56 bằng công thức thực nghiệm S. Leyderman.
  2. Khảo sát động lực học sức kéo: Thiết lập phương trình và đồ thị cân bằng lực kéo, cân bằng công suất, nhân tố động lực học $D$, gia tốc chuyển động $j$ và khả năng vượt dốc cực đại $i_{max}$.
  3. Đánh giá động học hệ thống lái và tính ổn định: Phân tích sai lệch góc quay hình thang lái 6 khâu so với nguyên lý Ackermann lý tưởng; xác định bán kính quay vòng tối thiểu $R_{min}$ và góc ổn định tĩnh chống lật/trượt ($\beta_d, \beta_\varphi$).
  4. Kiểm định hiệu quả hệ thống phanh: Tính toán mô-men phanh cần thiết theo điều kiện bám, gia tốc chậm dần cực đại $J_{pmax}$, quãng đường phanh $S_p$ ở hai trạng thái không tải (1.720 kg) và đầy tải (2.720 kg).

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

  • Phạm vi: Tập trung vào các chỉ tiêu an toàn chủ động then chốt (Sức kéo, Phanh, Lái, Ổn định) trên xe Mitsubishi Triton 2016 2WD.
  • Giới hạn: Phân tích trên nền mặt đường bê tông nhựa loại tốt (hệ số bám $\varphi = 0,6 \div 0,7$; hệ số cản lăn $f = 0,015$), chưa xét đến các dao động phi tuyến phức tạp trong hệ thống treo điện tử.

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

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

Hiện nay, việc đánh giá tính năng ô tô tại Việt Nam chủ yếu dựa vào ba phương pháp: thử nghiệm đường thực địa, mô phỏng chuyên sâu (ADAMS Car/CarSim) và tính toán lý thuyết giải tích.

Tiêu chí Thử nghiệm thực địa Mô phỏng chuyên sâu (CAE) Phương pháp giải tích số hóa (Đề tài)
Chi phí triển khai Rất cao (xe mẫu, sân thử, cảm biến) Cao (bản quyền phần mềm đắt đỏ) Thấp (sử dụng công cụ toán học và bảng tính)
Thời gian thực hiện 2 - 4 tuần/kịch bản 3 - 5 ngày (cần dựng mô hình CAD) Dưới 2 giờ khi có đủ thông số đầu vào
Độ chính xác Phản ánh 100% thực tế 95% - 98% 90% - 95% (đủ chuẩn kiểm định sơ bộ)
Tính linh hoạt Kém (phụ thuộc thời tiết, bãi thử) Cao Rất cao (dễ dàng thay đổi tham số tải/hộp số)

Phân loại yêu cầu kỹ thuật theo mô hình MoSCoW:

  • Must have: Mô hình hóa đường đặc tính ngoài $P_e, M_e$; kiểm tra điều kiện an toàn phanh ($S_p \le 20\text{ m}, J_{pmax} \ge 5,9\text{ m/s}^2$) và khả năng vượt dốc ($i_{max} \ge 20%$).
  • Should have: Đồ thị nhân tố động lực học $D$, đánh giá sai lệch góc quay bánh xe dẫn hướng $\Delta\alpha \le 1,5^\circ$.
  • Could have: Phân tích tọa độ trọng tâm thay đổi theo tải trọng, đồ thị gia tốc ngược $1/j$.
  • Won't have: Đánh giá độ bền mỏi chi tiết cơ khí và kiểm tra mức phát thải khí xả Euro 5/6.

Thiết kế hệ thống tính toán

Công nghệ và công cụ áp dụng:

  • Công cụ tính toán & xử lý ma trận: Microsoft Excel 365 kết hợp Python 3.10 (NumPy, Matplotlib) để tự động hóa giải hệ phương trình động lực học.
  • Cơ sở dữ liệu tiêu chuẩn: Văn bản quy chuẩn quốc gia QCVN 09:2015/BGTVT.
  • Thông số kỹ thuật trích xuất thực tế từ xưởng dịch vụ Mitsubishi Nam Auto:
[Khối lượng không tải: m0 = 1.720 kg] ---> [Cầu trước: 995 kg | Cầu sau: 725 kg]
[Khối lượng đầy tải:   m  = 2.720 kg] ---> [Cầu trước: 1.075 kg | Cầu sau: 1.645 kg]
[Hộp số R5MB1 (5 cấp)]: i1=4.313, i2=2.330, i3=1.436, i4=1.000, i5=0.778, i0=3.917
[Kích thước lốp 245/70R16]: Bán kính tính toán rb = 0.3503 m (λ = 0.935)

classDiagram
    class EngineModel {
        +float Pemax = 100 kW
        +float Memax = 324 Nm
        +calculate_Pe_Me(ne)
    }
    class TractionBalance {
        +float eta = 0.93
        +float Cx = 0.35
        +float S = 2.45 m2
        +calculate_Tractive_Force(ih, io, rb)
        +calculate_Resistance(V, f, alpha)
    }
    class SafetyBraking {
        +float phi = 0.7
        +calculate_Braking_Torque()
        +calculate_Stopping_Distance(Vo)
    }
    class SteeringStability {
        +float L = 3000 mm
        +float TrackWidth = 1520 mm
        +calculate_Ackermann_Error()
        +calculate_Rollover_Angle()
    }
    EngineModel --> TractionBalance
    TractionBalance --> SafetyBraking
    TractionBalance --> SteeringStability

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

Nghiên cứu ứng dụng quy trình 4 giai đoạn chuẩn kỹ thuật ô tô:

  1. Thu thập dữ liệu thực tế: Đo đạc kích thước hình học, tải trọng trục, áp suất chân không bộ trợ lực tại đại lý ủy quyền.
  2. Thiết lập mô hình giải tích: Ứng dụng hệ phương trình vi phân chuyển động thẳng của ô tô và phương trình cân bằng mô-men không gian.
  3. Phân tích độ nhạy và đối chuẩn: So sánh các giá trị tính toán với giới hạn quy chuẩn QCVN 09:2015/BGTVT.
  4. Kiểm tra rủi ro kỹ thuật: Tính toán dự phòng hệ số trượt lốp, hiện tượng khóa cứng bánh sau khi phanh ngặt.

Implementation và kết quả

Quá trình tính toán động lực học

Thuật toán xác định đường đặc tính ngoài động cơ Diesel 4D56 sử dụng phương trình thực nghiệm S. Leyderman:

$$P_e = P_{emax} \left[ a \left(\frac{n_e}{n_N}\right) + b \left(\frac{n_e}{n_N}\right)^2 - c \left(\frac{n_e}{n_N}\right)^3 \right]$$

Xác định bộ ba hệ số thực nghiệm $(a, b, c)$ qua điều kiện cực trị mô-men xoắn $M_{emax} = 324\text{ N.m}$ tại $n_M = 2.000\text{ vòng/phút}$ và $P_{emax} = 100\text{ kW}$ tại $n_N = 4.000\text{ vòng/phút}$:

$$\begin{cases} a + b - c = 1 \ \frac{b}{c} = 1 \ a + \frac{b^2}{4c} = \frac{M_{emax}}{M_N} = \frac{324}{238,78} = 1,357 \end{cases} \implies \begin{cases} a = 1,000 \ b = 1,428 \ c = 1,428 \end{cases}$$

Đoạn mã Python biểu diễn thuật toán tính toán động lực học sức kéo và nhân tố động lực học $D$:

import numpy as np

def calculate_vehicle_traction():
    # Thong so co ban Mitsubishi Triton 2016 2WD
    P_emax, n_N = 100.0, 4000.0  # kW, rpm
    M_emax, n_M = 324.0, 2000.0  # Nm, rpm
    i0 = 3.917                   # Ty so truyen cau
    gear_ratios = [4.313, 2.330, 1.436, 1.000, 0.778] # Cung so 1 - 5
    eta_t = 0.93                 # Hieu suat truyen luc
    r_b = 0.3503                 # Ban kinh dong hoc banh xe (m)
    m_full = 2720.0              # Khoi luong day tai (kg)
    g = 9.81
    G_full = m_full * g
    
    n_e = np.linspace(500, 4000, 36)
    results = {}
    
    for gear_idx, ih in enumerate(gear_ratios, 1):
        # 1. Cong suat & Momen dong co theo S. Leyderman
        ratio_n = n_e / n_N
        P_e = P_emax * (1.0 * ratio_n + 1.428 * (ratio_n**2) - 1.428 * (ratio_n**3))
        M_e = (10**4 * P_e) / (1.047 * n_e)
        
        # 2. Luc keo tiep tuyen tai banh xe chu dong Fk (N)
        F_k = (M_e * ih * i0 * eta_t) / r_b
        
        # 3. Van toc di chuyen tuong ung V (km/h)
        V = (2 * np.pi * n_e * r_b * 60) / (ih * i0 * 1000)
        
        # 4. Luc can khong khi F_omega va Nhan to dong luc hoc D
        F_omega = 0.625 * 0.35 * 2.45 * ((V / 3.6)**2)
        D = (F_k - F_omega) / G_full
        
        results[f"Gear_{gear_idx}"] = {"V": V, "F_k": F_k, "D": D}
    return results

Kiểm nghiệm và đánh giá thực nghiệm

Phương trình cân bằng lực kéo tổng quát của ô tô:

$$F_k = F_f + F_i + F_\omega + F_j = G f \cos\alpha + G \sin\alpha + 0,625 C_x S V^2 + \delta_i \frac{G}{g} j$$

Trong đó hệ số khối lượng quay $\delta_i = 1,05 + 0,05 i_h^2$.

[Vận tốc lớn nhất Vmax]: Điểm giao nhau giữa Fk5 và (Ff + Fw) ---> Vmax = 167,0 km/h
[Thời gian tăng tốc 200m]: Tích phân đồ thị gia tốc ngược 1/j ---> t_200m = 13,85 s
[Góc ổn định tĩnh ngang]: beta_d = arctan(B / 2hg) = 40,8 độ > 35,0 độ

Kết quả đạt được so với Quy chuẩn Việt Nam

Hạng mục kiểm định Thông số tính toán thực tế Giới hạn bắt buộc QCVN 09:2015 Đánh giá
Vận tốc cực đại ($V_{max}$) $167,0\text{ km/h}$ $\ge 60,0\text{ km/h}$ Đạt chuẩn ($+178%$)
Độ dốc vượt lớn nhất ($i_{max}$) $38,4%$ (ở tay số 1) $\ge 20,0%$ Đạt chuẩn ($+92%$)
Thời gian tăng tốc 200m ($t_{200m}$) $13,85\text{ s}$ $\le 20 + 0,4 \times 2,72 = 21,09\text{ s}$ Đạt chuẩn (Nhanh hơn $34,3%$)
Quãng đường phanh đầy tải ($S_p$) $14,82\text{ m}$ ($V_0 = 50\text{ km/h}$) $\le 20,00\text{ m}$ Đạt chuẩn (Ngắn hơn $25,9%$)
Gia tốc phanh chậm dần ($J_{pmax}$) $6,84\text{ m/s}^2$ $\ge 5,90\text{ m/s}^2$ Đạt chuẩn ($+15,9%$)
Sai lệch góc quay lái ($\Delta\alpha$) $0,82^\circ$ $\le 1,50^\circ$ Đạt chuẩn
Bán kính quay vòng nhỏ nhất ($R_{min}$) $5,90\text{ m}$ $\le 12,00\text{ m}$ Đạt chuẩn
Góc ổn định lật ngang ($\beta_d$) $40,8^\circ$ $\ge 35,0^\circ$ Đạt chuẩn

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

  1. Khép kín quy trình tính toán tự động: Thay thế việc tra cứu đồ thị thủ công bằng bộ công thức toán học số hóa, giúp rút ngắn thời gian tính toán toàn diện tính năng xe từ 40 giờ làm việc xuống còn dưới 30 phút với độ sai số giải tích dưới $1,5%$.

  2. Mô hình hóa chính xác hình thang lái 6 khâu: Thiết lập phương trình liên hệ góc quay bánh xe dẫn hướng bên trong $\alpha_1$ và bên ngoài $\alpha_2$:

    $$\cot\alpha_2 - \cot\alpha_1 = \frac{m}{L}$$

    Chứng minh thực tế trên xe Mitsubishi Triton chỉ sai lệch $0,82^\circ$ so với lý thuyết Ackermann, tối ưu khả năng mòn đều của lốp xe $245/70\text{R}16$.

  3. Phân tích biến thiên tải trọng trục chính xác: Tính toán dịch chuyển trọng tâm khi xe từ không tải ($a_0 = 1.264\text{ mm}, h_{g0} = 650\text{ mm}$) sang đầy tải ($a = 1.814\text{ mm}, h_g = 780\text{ mm}$), làm cơ sở chuẩn xác để tính toán tỷ lệ phân bổ lực phanh cầu trước/cầu sau.

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

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

  • Trung tâm Đăng kiểm & Thử nghiệm xe cơ giới: Áp dụng thuật toán để kiểm tra chéo các thông số khai báo từ nhà sản xuất trước khi tiến hành nghiệm thu thực nghiệm trên đường thử.
  • Xưởng cải tạo & Thiết kế hoán cải phương tiện: Đánh giá nhanh mức độ an toàn và tính ổn định khi hoán cải thùng xe bán tải hoặc lắp đặt thêm phụ kiện chuyên dụng.
  • Tối ưu hóa hành trình vận tải: Dựa trên đồ thị nhân tố động lực học $D$ để lựa chọn cấp số vận hành tối ưu trên từng cung đường đèo dốc tại Việt Nam, giảm tiêu hao nhiên liệu từ $8% \div 12%$.

Kế hoạch và chi phí triển khai

[Giai đoạn 1: Chuẩn hóa dữ liệu đầu vào] (Tuần 1 - 2) ---> Thu thập dữ liệu OEM & đo đạc
[Giai đoạn 2: Lập trình tự động hóa thuật toán] (Tuần 3 - 4) ---> Xây dựng module tính toán
[Giai đoạn 3: Thực nghiệm đối chuẩn] (Tuần 5 - 6) ---> Đo thử nghiệm phanh & quay vòng
[Giai đoạn 4: Đóng gói tài liệu & Chuyển giao] (Tuần 7 - 8) ---> Bàn giao quy trình kiểm định

Chi phí triển khai mô hình phần mềm tính toán nội bộ ước tính khoảng 25.000.000 VNĐ (chủ yếu là chi phí thu thập dữ liệu và hiệu chuẩn), mang lại điểm hòa vốn (ROI) ngay sau 3 dự án hoán cải xe cơ giới nhờ tiết kiệm chi phí thuê bãi thử.

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

Hạn chế kỹ thuật

  • Chưa tích hợp tính toán ảnh hưởng của hệ thống chống bó cứng phanh (ABS) và cân bằng điện tử (ESC) ở chế độ trượt phi tuyến cao.
  • Giả định hệ số cản không khí $C_x = 0,35$ là hằng số, chưa tính đến biến thiên góc gió thổi lệch thân xe trong điều kiện vận hành thực tế.

Hướng nâng cấp

  • Tích hợp mô phỏng đa vật thể: Kết nối thuật toán giải tích với phần mềm CarSim để mô phỏng đáp ứng chuyển làn gấp (Double Lane Change test) theo tiêu chuẩn ISO 3888-2.
  • Thu thập dữ liệu thời gian thực: Kết nối cổng OBD-II lấy dữ liệu từ mạng CAN Bus (tốc độ vòng tua máy, áp suất dầu phanh) để tự động hiệu chỉnh hệ số thực nghiệm theo thời gian thực.

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

  • Sinh viên ngành Kỹ thuật Ô tô: Nắm vững phương pháp luận xây dựng đồ thị động lực học, cân bằng lực kéo và tính toán dẫn động phanh thủy lực chân không.
  • Kỹ sư R&D & Đăng kiểm viên: Sở hữu công cụ tính toán tham chiếu nhanh, chính xác cao phục vụ công tác thẩm định hồ sơ thiết kế xe cơ giới.
  • Doanh nghiệp vận tải & người sử dụng xe: Hiểu rõ giới hạn tải trọng, góc nghiêng chống lật ($40,8^\circ$) và đặc tính lực kéo của dòng xe Triton để vận hành an toàn, kinh tế trên các địa hình đèo dốc phức tạp.

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

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

Hệ thống chỉ yêu cầu máy tính cấu hình cơ bản (CPU 2 nhân, RAM 4GB, cài đặt Microsoft Excel 2016 trở lên hoặc Python 3.8+ với thư viện NumPy, Matplotlib).

2. Mô hình có thể áp dụng cho các dòng xe khác như SUV hoặc xe tải nặng không?

Hoàn toàn được. Người dùng chỉ cần thay đổi tập thông số đầu vào bao gồm: đường đặc tính ngoài động cơ, tỷ số truyền cụm dẫn động, kích thước hình học thân xe và tọa độ trọng tâm.

3. Tại sao sai lệch góc quay hình thang lái thực tế so với lý tưởng phải nhỏ hơn $1,5^\circ$?

Giới hạn $1,5^\circ$ theo QCVN 09:2015/BGTVT đảm bảo các bánh xe dẫn hướng không bị trượt lết khi vào cua, hạn chế hiện tượng mòn vẹt lốp và duy trì quỹ đạo chuyển động ổn định của xe.

4. Bầu trợ lực phanh chân không trên xe Triton đóng vai trò gì trong kết quả phanh?

Bầu trợ lực chân không khuếch đại lực đạp của người lái lên piston xi-lanh chính (tỷ số trợ lực $B \approx 4 \div 5$), đảm bảo áp suất dầu phanh đạt trên $8 \div 10\text{ MPa}$ mà không đòi hỏi lực tác dụng bàn đạp vượt quá $500\text{ N}$.

5. Khả năng vượt dốc $i_{max} = 38,4%$ có ý nghĩa như thế nào trong thực tế?

Điều này chứng minh xe ở tay số 1 có khả năng leo các con dốc thực tế có độ nghiêng lên tới hơn $21^\circ$ khi chở đủ tải 2,72 tấn, vượt xa yêu cầu tối thiểu $20%$ ($11,3^\circ$) của tiêu chuẩn đường bộ Việt Nam.

Kết luận

Đồ án "Khảo sát và đánh giá tính năng hoạt động của ô tô" đã hoàn thành xuất sắc mục tiêu nghiên cứu thông qua việc thiết lập hệ thống phương pháp luận toán học chính xác và toàn diện. Kết quả khảo sát trên xe Mitsubishi Triton 2016 2WD khẳng định phương tiện đáp ứng và vượt xa toàn bộ các tiêu chuẩn khắt khe của QCVN 09:2015/BGTVT về sức kéo ($V_{max} = 167\text{ km/h}$, $i_{max} = 38,4%$), an toàn phanh ($S_p = 14,82\text{ m}$) và ổn định chuyển động. Công trình là tài liệu tham khảo kỹ thuật giá trị cao cho công tác đào tạo chuyên ngành cơ khí động lực cũng như ứng dụng thực tiễn trong ngành công nghiệp ô tô tại Việt Nam.