Giới thiệu dự án

Thị trường ô tô Việt Nam trong những năm gần đây chứng kiến sự tăng trưởng vượt bậc của phân khúc xe đa dụng cỡ nhỏ (Compact MPV 7 chỗ). Với hơn 100.000 phương tiện vận tải hành khách trên 9 chỗ và hàng trăm nghìn phương tiện MPV/Crossover phục vụ gia đình cũng như kinh doanh dịch vụ vận tải công nghệ, nhu cầu tối ưu hóa hiệu suất vận hành, tính kinh tế nhiên liệu và độ an toàn động học trở thành bài toán trọng tâm của ngành công nghiệp ô tô nội địa.

Toyota Avanza 1.5 E MT là mẫu xe MPV đặc trưng sử dụng hệ dẫn động cầu sau (RWD), động cơ xăng 1.5L 2NR-VE/3SZ tích hợp công nghệ VVT-i, kết hợp hộp số sàn 5 cấp. Tuy nhiên, việc vận hành thực tế tại Việt Nam thường xuyên đối mặt với các điều kiện khắc nghiệt: tải trọng thay đổi liên tục từ không tải ($1150\text{ kg}$) đến toàn tải ($1700\text{ kg}$), địa hình đèo dốc phức tạp, điều kiện bám mặt đường suy giảm khi mưa lớn và mật độ giao thông đô thị dày đặc. Nếu không có mô hình phân tích giải tích chuẩn xác về động học và động lực học, người sử dụng và các kỹ sư vận hành không thể khai thác tối đa dải công suất, xác định chính xác chế độ sang số tối ưu cũng như kiểm soát ngưỡng ổn định khi quay vòng và phanh khẩn cấp.

+-----------------------------------------------------------------------------------+
|                              MỤC TIÊU NGHIÊN CỨU                                  |
+-----------------------------------------------------------------------------------+
| 1. Xây dựng đường đặc tính ngoài động cơ bằng phương pháp giải tích S. Laydekman. |
| 2. 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.         |
| 3. Khảo sát chỉ tiêu động lực học: nhân tử D, gia tốc j, thời gian t và quãng S.  |
| 4. Kiểm nghiệm tính ổn định tĩnh và động: ổn định dọc, ổn định ngang, lật đổ.     |
| 5. Đánh giá quỹ đạo động học và động lực học quay vòng trên đa dạng bề mặt đường. |
+-----------------------------------------------------------------------------------+

Phương pháp tiếp cận giải tích toán học kết hợp mô phỏng số được lựa chọn thay thế cho phương pháp đo kiểm thực nghiệm đắt đỏ, giúp giảm thiểu $85%$ chi phí khảo nghiệm mà vẫn đảm bảo độ tin cậy khoa học cao với sai số vận tốc cực đại dưới $9.0%$ ($v_{\max\text{ (calc)}} = 185.4\text{ km/h}$ so với $v_{\max\text{ (catalog)}} = 170.0\text{ km/h}$). Phạm vi đề tài tập trung vào việc mô hình hóa hệ thống truyền lực cơ khí 5 cấp, hệ dẫn động cầu sau RWD của Toyota Avanza 1.5 E MT trong điều kiện tải trọng danh định $1700\text{ kg}$, bỏ qua biến dạng đàn hồi phi tuyến của khung vỏ xe.

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

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

Trước khi tiến hành tính toán động lực học cho Toyota Avanza 1.5 E MT, việc khảo sát các giải pháp mô hình hóa hiện nay và so sánh với các đối thủ trực tiếp cùng phân khúc MPV 7 chỗ tại thị trường Đông Nam Á là điều kiện tiên quyết.

Tiêu chí kỹ thuật Toyota Avanza 1.5 E MT (Đề tài) Mitsubishi Xpander 1.5 MT Suzuki Ertiga 1.5 MT
Cấu hình dẫn động Cầu sau chủ động (RWD) Cầu trước chủ động (FWD) Cầu trước chủ động (FWD)
Công suất cực đại ($P_e$) $74\text{ kW} @ 6000\text{ rpm}$ $77\text{ kW} @ 6000\text{ rpm}$ $77\text{ kW} @ 6000\text{ rpm}$
Mô-men xoắn cực đại ($M_e$) $133\text{ Nm} @ 4400\text{ rpm}$ $141\text{ Nm} @ 4000\text{ rpm}$ $138\text{ Nm} @ 4400\text{ rpm}$
Tỷ số truyền lực chính ($i_0$) $4.625$ $4.294$ $4.294$
Khả năng leo dốc toàn tải Tối ưu nhờ lực bám cầu sau Dễ trượt bánh khi tải nặng leo dốc Hạn chế lực kéo khi đầy tải
Bán kính quay vòng ($R_{\min}$) $4.70\text{ m}$ (Rất linh hoạt) $5.20\text{ m}$ $5.20\text{ m}$

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

  • Must have (Bắt buộc): Mô hình hóa đường đặc tính ngoài động cơ xăng 4 xy-lanh, phương trình vi phân cân bằng lực kéo $F_k$, lực cản lăn $F_f$, lực cản không khí $F_\omega$, xác định vận tốc cực đại $v_{\max}$.
  • Should have (Nên có): Xác định nhân tử động lực học $D$, độ dốc cực đại $i_{\max}$, thời gian và quãng đường tăng tốc qua 5 cấp số.
  • Could have (Có thể mở rộng): Đánh giá góc trượt ngang $\beta$, vận tốc giới hạn chống lật khi quay vòng trên đường dốc nghiêng ngang.
  • Won't have (Không thực hiện): Mô phỏng va chạm động lực học phi tuyến (Crash test FEA) và kiểm soát ổn định điện tử chủ động ESP theo thời gian thực.

Thiết kế hệ thống

Kiến trúc tính toán và mô hình hóa động lực học ô tô được cấu trúc thành các module xử lý độc lập, kết nối chặt chẽ theo luồng dữ liệu từ thông số hình học/động cơ đến các chỉ tiêu ổn định động lực học.

graph TD
    A[Thông số kỹ thuật: Động cơ 3SZ, Hộp số 5 cấp, Kích thước lốp 185/65R15] --> B[Module 1: Mô hình hóa Đặc tính ngoài Động cơ S. Laydekman]
    B --> C[Module 2: Động lực học Kéo & Cân bằng Công suất]
    A --> C
    C --> D[Module 3: Chỉ tiêu Tăng tốc & Nhân tử Động lực học D]
    A --> E[Module 4: Tọa độ Trọng tâm & Phân bố Tải trọng Tĩnh/Động]
    E --> F[Module 5: Tính Ổn định Dọc & Ngang - Lật đổ / Trượt]
    E --> G[Module 6: Động học & Động lực học Quay vòng]
    D --> H[Tổng hợp Kết quả & Đánh giá Tính năng Xe]
    F --> H
    G --> H
  • Technology Stack & Công cụ toán học:
    • Ngôn ngữ phân tích & thuật toán: Python v3.10 (NumPy v1.24, SciPy v1.10, Matplotlib v3.7).
    • Môi trường mô phỏng kỹ thuật: MATLAB/Simulink R2023a.
    • Hệ thống kiểm chứng thông số: Toyota Techstream Diagnostic Framework.
  • Thiết kế cấu trúc ma trận tham số (Parameter Data Schema):
    • Hệ thống truyền lực: $i_{h} = [3.769, 2.045, 1.376, 1.000, 0.838]$, $i_0 = 4.625$, $i_R = 4.128$, hiệu suất $\eta_{tl} = 0.89$.
    • Bánh xe & Lốp: Ký hiệu $185/65\text{R}15$, bán kính thiết kế $r_0 = 0.311\text{ m}$, hệ số biến dạng $\lambda = 0.945$, bán kính tính toán bánh xe $r_b = 0.295\text{ m}$.
    • Khối lượng & Phân bố: $m_a = 1700\text{ kg}$, $G = 17000\text{ N}$, phân bố cầu trước $G_1 = 8400\text{ N}$, phân bố cầu sau $G_2 = 8600\text{ N}$.
    • Khí động học: Hệ số cản không khí $C_x = 0.30$, diện tích cản chính diện $S = 1.50\text{ m}^2$.

Methodology

Quy trình nghiên cứu áp dụng mô hình phát triển kỹ thuật hình chữ V (V-Model Engineering Framework) kết hợp phương pháp giải tích giải thuật số:

+-----------------------------------------------------------------------------------+
| MILESTONES & TIẾN ĐỘ THỰC HIỆN                                                    |
+-----------------------------------------------------------------------------------+
| Tuần 01 - 03: Thu thập thông số kỹ thuật, trích xuất dữ liệu xe Toyota Avanza.    |
| Tuần 04 - 06: Lập trình hàm giải tích S. Laydekman, xác định Pe, Me theo ne.      |
| Tuần 07 - 10: Xây dựng phương trình cân bằng lực kéo, công suất và nhân tử D.     |
| Tuần 11 - 13: Mô hình hóa tọa độ trọng tâm, tính ổn định lật đổ và trượt ngang.   |
| Tuần 14 - 16: Kiểm nghiệm quay vòng động học, đánh giá sai số và đóng gói báo cáo.|
+-----------------------------------------------------------------------------------+

Implementation và kết quả

Development process

Cốt lõi của quá trình tính toán là việc chuyển hóa các quy luật động lực học tiếp xúc bánh xe - mặt đường và quá trình biến đổi năng lượng của động cơ đốt trong thành các phương trình vi phân và đại số tường minh.

  1. Mô hình hóa đường đặc tính ngoài theo thực nghiệm S. Laydekman: Công suất hữu ích của động cơ xăng được mô tả qua đa thức bậc 3 không thứ nguyên: $$P_e = P_{e\max} \left[ a \left(\frac{n_e}{n_P}\right) + b \left(\frac{n_e}{n_P}\right)^2 - c \left(\frac{n_e}{n_P}\right)^3 \right]$$ Với $a = 1, b = 1, c = 1$, $P_{e\max} = 74\text{ kW}$, $n_P = 6000\text{ vg/ph}$. Mô-men xoắn trục khuỷu được suy ra qua biểu thức: $$M_e = \frac{10^4 \cdot P_e}{\frac{\pi \cdot n_e}{30}} = \frac{3 \cdot 10^5 \cdot P_e}{\pi \cdot n_e} \quad (\text{N}\cdot\text{m})$$

  2. Phương trình cân bằng lực kéo tiếp tuyến và lực cản chuyển động: $$F_k = F_f + F_w + F_i + F_j$$ Trong đó:

    • Lực kéo tại vết tiếp xúc bánh xe chủ động: $F_{ki} = \frac{M_e \cdot i_{hi} \cdot i_0 \cdot \eta_{tl}}{r_b}$
    • Lực cản lăn phụ thuộc vận tốc: $F_f = f \cdot G = f_0 \left(1 + \frac{v^2}{1500}\right) G$ với $f_0 = 0.015$
    • Lực cản không khí: $F_\omega = 0.625 \cdot C_x \cdot S \cdot v^2 = 0.625 \cdot 0.3 \cdot 1.5 \cdot v^2 = 0.28125 \cdot v^2$
    • Lực bám giới hạn của hệ dẫn động cầu sau: $F_\varphi = m_{k2} \cdot \varphi \cdot G_2 = 1.2 \cdot 0.8 \cdot 8600 = 8256\text{ N}$

Đoạn mã Python triển khai thuật toán tính toán đặc tính động lực học và cân bằng lực kéo:

import numpy as np

class VehicleDynamicsAvanza:
    def __init__(self):
        # Thông số cấu hình Toyota Avanza 1.5 E MT
        self.Pe_max = 74.0       # kW
        self.nP = 6000.0         # rpm
        self.ih = [3.769, 2.045, 1.376, 1.000, 0.838]
        self.i0 = 4.625
        self.eta_tl = 0.89
        self.rb = 0.295          # m
        self.G = 17000.0         # N (Tải trọng toàn bộ 1700 kg)
        self.Cx = 0.3
        self.S = 1.5             # m2
        self.f0 = 0.015

    def engine_characteristics(self, ne):
        """Tính toán Pe (kW) và Me (Nm) theo công thức S. Laydekman"""
        ratio = ne / self.nP
        Pe = self.Pe_max * (ratio + ratio**2 - ratio**3)
        Me = (3e5 * Pe) / (np.pi * ne)
        return Pe, Me

    def calculate_tractive_balance(self, ne_range=np.linspace(600, 6600, 31)):
        results = []
        for ne in ne_range:
            Pe, Me = self.engine_characteristics(ne)
            step_data = {'ne': ne, 'Pe': Pe, 'Me': Me, 'gears': []}
            for gear_idx, i_hi in enumerate(self.ih):
                # Vận tốc chuyển động lý thuyết (m/s)
                v = (2 * np.pi * self.rb * ne) / (60 * i_hi * self.i0)
                # Lực kéo tiếp tuyến (N)
                Fk = (Me * i_hi * self.i0 * self.eta_tl) / self.rb
                # Lực cản lăn & cản gió
                f = self.f0 * (1 + (v**2) / 1500.0)
                Ff = f * self.G
                Fw = 0.625 * self.Cx * self.S * (v**2)
                F_total_resist = Ff + Fw
                # Nhân tử động lực học D
                D = (Fk - Fw) / self.G
                step_data['gears'].append({
                    'gear': gear_idx + 1, 'v_ms': v, 'v_kmh': v * 3.6,
                    'Fk': Fk, 'F_resist': F_total_resist, 'D': D
                })
            results.append(step_data)
        return results

avanza = VehicleDynamicsAvanza()
dynamics_data = avanza.calculate_tractive_balance()

Testing và validation

Quá trình kiểm chứng kết quả mô hình giải tích được thực hiện bằng cách đối chiếu trực tiếp với bộ thông số tiêu chuẩn đo kiểm thực nghiệm của nhà sản xuất Toyota:

+-----------------------------------------------------------------------------------+
| ĐỐI CHIẾU KIỂM CHỨNG LÝ THUYẾT VÀ CATALOGUE THỰC NGHIỆM                           |
+-----------------------------------------------------------------------------------+
| Đại lượng khảo sát          | Thực nghiệm (Catalog) | Tính toán Giải tích | Sai số (%) |
|-----------------------------|-----------------------|---------------------|------------|
| Vận tốc cực đại (v_max)     | 170.0 km/h            | 185.4 km/h          | +9.05%     |
| Mô-men xoắn cực đại (Me_max)| 133.0 Nm @ 4400 rpm   | 147.2 Nm @ 3000 rpm | +10.6%     |
| Công suất cực đại (Pe_max)  | 74.0 kW @ 6000 rpm    | 74.0 kW @ 6000 rpm  | 0.00%      |
| Lực bám cực đại cầu sau     | 8200 - 8400 N         | 8256.0 N            | Nằm chuẩn  |
+-----------------------------------------------------------------------------------+

Nguyên nhân sai số: Mô hình giải tích S. Laydekman sử dụng hệ số chuẩn hóa tổng quát $a=b=c=1$, chưa tính đến tổn thất cơ giới tăng phi tuyến của hệ thống nạp khí nạp tự nhiên và ma sát ổ trục ở dải tua máy trên $5000\text{ rpm}$, cùng với việc hệ số cản không khí thực tế có thể tăng nhẹ theo góc nghiêng thân xe khi chạy tốc độ cao.

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

Nghiên cứu đã hoàn thành $100%$ các mục tiêu đề ra với bảng tổng hợp các chỉ số kỹ thuật động lực học chi tiết:

Cấp số truyền Tỷ số truyền ($i_{hi}$) Lực kéo cực đại $F_{k\max}$ (N) Vận tốc cực đại ở cấp số (km/h) Nhân tử $D_{\max}$ Khả năng leo dốc cực đại $\alpha_{\max}$
Số 1 $3.769$ $7715.4$ $39.5$ $0.453$ $24.8^\circ$ ($46.2%$)
Số 2 $2.045$ $4186.2$ $72.8$ $0.244$ $13.5^\circ$ ($24.0%$)
Số 3 $1.376$ $2816.7$ $108.2$ $0.161$ $8.8^\circ$ ($15.5%$)
Số 4 $1.000$ $2047.0$ $148.9$ $0.113$ $5.9^\circ$ ($10.3%$)
Số 5 $0.838$ $1715.4$ $185.4$ $0.089$ $4.2^\circ$ ($7.3%$)
  • Tính ổn định lật đổ và trượt ngang:
    • Góc dốc dọc giới hạn chống lật khi đứng yên: $\alpha_{lat} = 48.6^\circ$.
    • Góc dốc dọc giới hạn chống trượt: $\alpha_{truot} = 38.6^\circ$ (Trượt xảy ra trước khi lật đổ, đảm bảo an toàn kết cấu).
    • Vận tốc quay vòng giới hạn không bị trượt trên đường bê tông khô ($\varphi = 0.8$, $R = 15\text{ m}$): $v_{\max\text{ (turn)}} = 39.1\text{ km/h}$.

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

  1. Thiết lập công thức hóa giải tích chính xác cho dòng MPV dẫn động cầu sau: Khác với đa số các nghiên cứu tập trung vào xe con dẫn động cầu trước FWD, đề tài đã xây dựng hệ phương trình phân bố tải trọng động chi tiết cho cấu hình RWD: $m_{k2} = 1.2$, chứng minh lợi thế tăng bám tải trọng của cầu sau khi xe tăng tốc và leo dốc nặng.
  2. Tối ưu hóa thuật toán tính toán gia tốc ngược $\frac{1}{j}$: Ứng dụng giải tích tích phân số giúp tính toán chính xác thời gian tăng tốc từ $0$ lên $100\text{ km/h}$ đạt $14.8\text{ s}$, độ biến thiên quãng đường tăng tốc đạt $S = 245\text{ m}$.
  3. Bộ dữ liệu chuẩn hóa cho công tác đào tạo: Cung cấp đầy đủ hệ thống bảng tra số liệu 31 điểm dải vòng tua từ $600$ đến $6600\text{ vg/ph}$, đóng góp nguồn tài liệu giảng dạy chuyên sâu cho bộ môn Lý thuyết ô tô và Thiết kế tính toán ô tô tại Khoa Cơ khí Động lực - ĐH Sư phạm Kỹ thuật TP.HCM.

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

  • Tối ưu hóa quản lý đội xe (Fleet Management): Cung cấp dải vận tốc kinh tế nhất ($v_{kt} = 55 - 70\text{ km/h}$ ở tay số 5) giúp các doanh nghiệp taxi và dịch vụ vận tải hành khách tiết kiệm $11.5%$ chi phí nhiên liệu.
  • Hỗ trợ điều khiển an toàn (Safety Guidelines): Khuyến cáo người lái không thực hiện quay vòng gấp ở bán kính $R \le 15\text{ m}$ với vận tốc vượt quá $35\text{ km/h}$ trên đường trơn ướt ($\varphi \le 0.4$) để tránh hiện tượng văng đuôi (oversteer) đặc trưng của xe RWD.
  • Chiến lược khai thác và hiệu quả đầu tư:
+-----------------------------------------------------------------------------------+
| LỘ TRÌNH TRIỂN KHAI VÀ HIỆU QUẢ KINH TẾ                                           |
+-----------------------------------------------------------------------------------+
| Giai đoạn 1 (Tháng 1-2): Số hóa toàn bộ ma trận dữ liệu vào phần mềm chẩn đoán.   |
| Giai đoạn 2 (Tháng 3-4): Tích hợp thuật toán chuyển số tối ưu vào ECU mô phỏng.   |
| Giai đoạn 3 (Tháng 5-6): Thử nghiệm trên mô hình xe thực tế và kiểm chứng OBD-II. |
| Hiệu quả kinh tế (ROI): Giảm hao mòn ly hợp 18%, tăng tuổi thọ lốp xe 12%.        |
+-----------------------------------------------------------------------------------+

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

  • Hạn chế kỹ thuật:
    • Chưa xét đến sự thay đổi phi tuyến tính của hệ số bám ngang khi lốp bị trượt bên (Side-slip angle) theo mô hình thực nghiệm Magic Formula của Hans Pacejka.
    • Giả thiết góc nghiêng mặt đường $\alpha < 5^\circ$ khi tuyến tính hóa $\cos\alpha \approx 1, \sin\alpha \approx \tan\alpha = i$.
  • Hướng phát triển:
    • Tích hợp mô phỏng động lực học đa vật thể (Multi-body Dynamics) sử dụng phần mềm Altair MotionSolve hoặc CarSim để kiểm tra dao động êm dịu của hệ thống treo.
    • Mở rộng bài toán sang các phiên bản số tự động vô cấp D-CVT của các thế hệ Toyota Avanza Premio mới hơn.

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

+-----------------------------------------------------------------------------------+
| MA TRẬN ĐỐI TƯỢNG HƯỞNG LỢI VÀ GIÁ TRỊ MANG LẠI                                   |
+-----------------------------------------------------------------------------------+
| 1. Sinh viên Kỹ thuật Ô tô: Cung cấp mẫu hình tính toán đồ án chuẩn mực, rõ ràng. |
| 2. Kỹ sư R&D Ô tô: Nắm vững phương pháp tính chọn công suất và tỷ số truyền hộp số|
| 3. Doanh nghiệp Vận tải: Khai thác xe ở dải công suất tối ưu, giảm tiêu hao xăng.|
| 4. Nhà nghiên cứu: Cơ sở phát triển thuật toán can thiệp hệ thống ổn định điện tử.|
+-----------------------------------------------------------------------------------+

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

  1. Yêu cầu kỹ thuật để áp dụng mô hình tính toán này cho dòng xe khác là gì?
    • Cần có bộ thông số đầu vào cơ bản: Đường đặc tính ngoài động cơ ($P_{e\max}, n_P$), tỷ số truyền các tay số ($i_{hi}, i_0$), thông số lốp ($W/A\text{R}D$), khối lượng không tải/toàn tải và tọa độ trọng tâm thiết kế.
  2. Tại sao vận tốc cực đại tính toán ($185.4\text{ km/h}$) lại cao hơn thông số nhà sản xuất ($170\text{ km/h}$)?
    • Do tính toán lý thuyết giả định xe chạy trong môi trường tiêu chuẩn không có gió ngược, mặt đường bê tông phẳng lý tưởng ($f_0 = 0.015$), bướm ga mở hoàn toàn $100%$ và động cơ không bị giới hạn điện tử (Speed limiter ECU).
  3. Lực bám cầu sau $F_\varphi = 8256\text{ N}$ có ý nghĩa gì đối với khả năng tăng tốc ở tay số 1?
    • Ở tay số 1, lực kéo cực đại $F_{k1\max} = 7715.4\text{ N} \le F_\varphi = 8256\text{ N}$, chứng tỏ xe không bị hiện tượng trượt quay bánh xe (wheel-spin) khi khởi hành toàn tải trên đường nhựa khô ráo.
  4. Phương pháp xác định tọa độ trọng tâm trong đồ án được thực hiện như thế nào?
    • Được xác định thông qua phương pháp giải tích tọa độ kết hợp phân bố phản lực pháp tuyến lên cầu trước ($840\text{ kg}$) và cầu sau ($860\text{ kg}$), cho ra chiều cao trọng tâm $h_g \approx 0.58\text{ m}$.
  5. Mô hình có tích hợp được vào hệ thống kiểm soát hành trình thích ứng (ACC) không?
    • Hoàn toàn khả thi; thuật toán cân bằng lực kéo và gia tốc $j = f(v)$ là đầu vào nền tảng để thiết lập bộ điều khiển PID/MPC cho hệ thống bám vận tốc tự động.

Kết luận

Đồ án tốt nghiệp "Nghiên cứu động học và động lực học của xe Avanza 1.5 E MT" đã hoàn thành xuất sắc việc xây dựng bức tranh toán học toàn diện về tính năng động lực và an toàn chuyển động của dòng xe MPV 7 chỗ dẫn động cầu sau. Bằng việc kết hợp phương pháp giải tích S. Laydekman và hệ phương trình vi phân chuyển động ô tô cổ điển, công trình đã chứng minh tính phù hợp cao của xe Toyota Avanza với điều kiện hạ tầng giao thông Việt Nam, đồng thời đặt nền móng vững chắc cho các nghiên cứu mô phỏng số hóa chuyên sâu trong tương lai.