Giới thiệu dự án

Công nghiệp sản xuất động cơ đốt trong (ĐCĐT) toàn cầu hiện đạt quy mô hơn 30 triệu động cơ mỗi năm, đóng vai trò xương sống trong các ngành giao thông vận tải, nông nghiệp và cơ khí động lực. Trước áp lực thắt chặt của các tiêu chuẩn phát thải quốc tế (Euro 6, Euro 7) và yêu cầu khắt khe về tiết kiệm năng lượng, việc tối ưu hóa chu trình nhiệt động và quá trình trao đổi khí trong buồng đốt là nhiệm vụ sống còn của các kỹ sư cơ khí động lực.

                  ┌───────────────────────────────┐
                  │ THÁCH THỨC ĐCĐT TRUYỀN THỐNG  │
                  │   Pha phối khí cố định        │
                  └──────────────┬────────────────┘
                                 │
         ┌───────────────────────┴───────────────────────┐
         ▼                                               ▼
┌─────────────────────────────────┐   ┌─────────────────────────────────┐
│       VÒNG TUA THẤP (LOW RPM)   │   │      VÒNG TUA CAO (HIGH RPM)    │
│ • Góc trùng điệp van quá lớn    │   │ • Thời gian mở van quá ngắn     │
│ • Dòng khí nạp bị thổi ngược    │   │ • Hệ số nạp giảm mạnh           │
│ • Tổn hao nhiên liệu, HC tăng   │   │ • Công suất cực đại bị bóp nghẹt│
└─────────────────────────────────┘   └─────────────────────────────────┘
                                 │
                                 ▼
                  ┌───────────────────────────────┐
                  │    GIẢI PHÁP: HONDA VTEC      │
                  │    & MÔ PHỎNG SỐ AVL BOOST    │
                  └───────────────────────────────┘

Trên các động cơ truyền thống, trục cam có góc mở sớm và đóng muộn cố định, thường chỉ được tối ưu hóa cho một chế độ tải trung bình. Khi động cơ hoạt động ở dải vòng tua thấp hoặc cao, hệ thống van cố định gây ra sự suy giảm đáng kể về hệ số nạp ($\eta_v$), tăng suất tiêu hao nhiên liệu riêng ($g_e$) và gia tăng nồng độ khí thải độc hại ($NO_x$, $CO$, $HC$). Công nghệ điều khiển thời điểm và độ nâng van biến thiên điện tử (Variable Valve Timing and Lift Electronic Control - VTEC) của Honda giải quyết triệt để vấn đề này bằng cách chuyển đổi giữa các vấu cam biên độ thấp và cao theo dải tốc độ.

Tuy nhiên, việc nghiên cứu, thử nghiệm và hiệu chỉnh cơ cấu phân phối khí VTEC trên băng thử động cơ thực tế đòi hỏi chi phí thiết bị đắt đỏ và thời gian kéo dài. Do đó, việc ứng dụng phần mềm mô phỏng nhiệt động học chuyên dụng AVL BOOST là giải pháp mang tính chiến lược, cho phép khảo sát toàn diện các chu trình làm việc, tối ưu hóa kích thước đường ống, pha phối khí và đánh giá chính xác các chỉ tiêu kinh tế - kỹ thuật.

Mục tiêu nghiên cứu

  1. Xây dựng mô hình nhiệt động học 1D hoàn chỉnh của động cơ 4 kỳ, 4 xy-lanh tăng áp Honda Civic 1.5L VTEC Turbo (mã L15B7) trên phần mềm AVL BOOST.
  2. Thiết lập và hiệu chỉnh chính xác thông số hình học, biên nhiệt động, động học xupap nạp - xả và quy luật cháy cho hệ thống VTEC ở nhiều chế độ làm việc.
  3. Khảo sát và đánh giá định lượng sự biến thiên của mô-men xoắn ($M_e$), công suất có ích ($N_e$), suất tiêu hao nhiên liệu ($g_e$) và nồng độ phát thải ($NO_x$, $CO$, $HC$) trên dải tốc độ từ 1000 đến 7000 vòng/phút.
  4. So sánh, chứng minh tính ưu việt của cơ cấu VTEC chuyển đổi pha tự động so với cơ cấu phối khí thông thường.

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

  • Đối tượng mô phỏng: Động cơ Honda L15B7 1.5L Turbo, DOHC 16 van i-VTEC, đường kính xy-lanh 86.0 mm, hành trình piston 85.9 mm, tỷ số nén 10.8:1.
  • Môi trường tính toán: Mô phỏng chu trình công tác 1D (One-Dimensional Thermodynamic Simulation) với giả thiết hòa khí hòa trộn hoàn hảo và dòng chảy trong đường ống là dòng khí nén tựa dừng.

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

Phân tích hiện trạng và giải pháp cạnh tranh

Tiêu chí kỹ thuật Cơ cấu phối khí cố định (Fixed Cam) Toyota VVT-i (Phát triển pha liên tục) Honda VTEC / i-VTEC (Đổi vấu cam + Pha)
Biến thiên thời điểm (Timing) Cố định ($\Delta \alpha = 0^\circ$) Biến thiên liên tục ($30^\circ - 50^\circ$ CA) Biến thiên theo nấc hoặc liên tục ($50^\circ$ CA)
Biến thiên độ nâng (Lift) Cố định ($h_{max} = \text{const}$) Cố định ($h_{max} = \text{const}$) Đa cấp (Low: $4.5-6.5\text{ mm}$, High: $11.99\text{ mm}$)
Hiệu suất nạp ở vòng tua cao Kém ($\eta_v < 0.78$) Trung bình ($\eta_v \approx 0.84$) Rất cao ($\eta_v \ge 0.92$)
Mô-men xoắn dải vòng tua thấp Trung bình Tốt Tối ưu hóa triệt để
Độ phức tạp cơ khí & Điều khiển Rất thấp Trung bình Cao (Thủy lực + Piston chốt đồng bộ)

Hệ thống VTEC nổi trội nhờ khả năng thay đổi đồng thời cả góc mở sớm/đóng muộn và chiều cao nâng xupap ($h_v$). Ở tốc độ thấp, vấu cam biên độ nhỏ duy trì vận tốc dòng khí nạp cao, tăng cường dòng xoáy lốc (swirl/tumble), cải thiện quá trình bay hơi nhiên liệu. Ở tốc độ cao, vấu cam biên độ lớn ($11.99\text{ mm}$) mở rộng tiết diện lưu thông, giảm tổn thất tiết lưu qua cửa nạp, giúp động cơ đạt công suất 171 mã lực tại 7000 vòng/phút.

Kiến trúc hệ thống mô phỏng trên AVL BOOST

Mô hình nhiệt động học 1D của động cơ 4 xy-lanh VTEC được thiết lập gồm các khối phần tử liên kết với nhau qua các đường ống dẫn khí nén:

graph LR
    SB1[Boundary Khí nạp SB1] --> CL1[Bộ lọc gió CL1]
    CL1 --> MP1[Plenum nạp PL1]
    MP1 --> R1[Runner 1] --> C1[Xy-lanh C1 + Phun I1]
    MP1 --> R2[Runner 2] --> C2[Xy-lanh C2 + Phun I2]
    MP1 --> R3[Runner 3] --> C3[Xy-lanh C3 + Phun I3]
    MP1 --> R4[Runner 4] --> C4[Xy-lanh C4 + Phun I4]
    C1 & C2 & C3 & C4 --> EX1[Ống gom xả EX1]
    EX1 --> CAT1[Bộ xử lý xúc tác CAT1]
    CAT1 --> SB2[Boundary Khí xả SB2]
    E1[Khối điều khiển động cơ E1] -.-> C1 & C2 & C3 & C4

Thông số cấu hình hệ thống (Technology Stack & Parameter Schema)

  • Phần mềm cốt lõi: AVL BOOST Release Version 2023.1 (Mô-đun mô phỏng động cơ 1D Gas Dynamics & Cycle Simulation).
  • Công cụ hậu xử lý: AVL BOOST Post-Processing Engine & IMPRESS 3D 2023.
  • Bộ tham số đầu vào của mô hình động cơ (Engine Master Configuration):
{
  "engine_configuration": {
    "engine_id": "E1",
    "engine_type": "4-Stroke Spark Ignition Turbocharged",
    "cylinders_count": 4,
    "bore_mm": 86.0,
    "stroke_mm": 85.9,
    "con_rod_length_mm": 143.4,
    "compression_ratio": 10.8,
    "firing_order": [1, 3, 4, 2],
    "firing_angles_deg_CA": [0.0, 540.0, 180.0, 360.0]
  },
  "gas_exchange_vtec": {
    "intake_valve": {
      "iv_open_deg_BTDC": 32.0,
      "iv_close_deg_ABDC": 48.0,
      "max_lift_mm": 11.99,
      "intake_duration_deg_CA": 260.0,
      "lash_mm": 0.22
    },
    "exhaust_valve": {
      "ev_open_deg_BBDC": 50.0,
      "ev_close_deg_ATDC": 20.0,
      "max_lift_mm": 11.00,
      "exhaust_duration_deg_CA": 250.0,
      "lash_mm": 0.20
    }
  },
  "combustion_injection": {
    "fuel_type": "Gasoline (C7H17)",
    "stoichiometric_AFR": 14.7,
    "actual_injector_AFR": 13.34,
    "injection_type": "Continuous / Port & Direct Injection",
    "initial_exhaust_pressure_bar": 5.0,
    "initial_exhaust_temp_celsius": 726.85
  }
}

Implementation và kết quả

Mô hình toán học và thuật toán nhiệt động học

AVL BOOST sử dụng phương pháp bảo toàn khối lượng và năng lượng cho hệ thống mở tại từng góc quay trục khuỷu ($\alpha$):

  1. Phương trình định luật thứ nhất nhiệt động lực học cho buồng cháy:

$$\frac{d(m \cdot u)}{d\alpha} = -p \cdot \frac{dV}{d\alpha} + \frac{dQ_B}{d\alpha} - \frac{dQ_w}{d\alpha} - h_{blowby} \cdot \frac{dm_{blowby}}{d\alpha} + \sum \dot{m}_i \cdot h_i$$

  1. Mô hình truyền nhiệt qua vách buồng đốt (Woschni Model):

$$\alpha_w = 127.93 \cdot B^{-0.2} \cdot p^{0.8} \cdot T^{-0.53} \cdot \left[ C_1 \cdot c_m + C_2 \cdot \frac{V_d \cdot T_1}{p_1 \cdot V_1} \cdot (p - p_{drag}) \right]^{0.8}$$

  1. Thuật toán tính toán biên dạng nâng van (Valve Lift Profile Interpolation): Độ nâng của xupap nạp ($h_{in}$) và xupap thải ($h_{ex}$) theo góc quay trục khuỷu ($\theta$) được mô hình hóa chính xác qua bảng giá trị thực nghiệm:
import numpy as np

def calculate_valve_lift(crank_angle, iv_open=328.0, iv_close=588.0, max_lift=11.99):
    """
    Tính toán độ nhấc xupap nạp động cơ VTEC theo góc quay trục khuỷu (deg CA)
    Góc mở sớm: 32 deg BTDC (tương đương 328 deg CA trên chu trình 720 deg)
    Góc đóng muộn: 48 deg ABDC (tương đương 588 deg CA)
    """
    if iv_open <= crank_angle <= iv_close:
        duration = iv_close - iv_open
        normalized_angle = (crank_angle - iv_open) / duration
        # Xấp xỉ biên dạng vấu cam VTEC biên độ cao bằng hàm hình sin tối ưu
        lift = max_lift * np.sin(np.pi * normalized_angle) ** 1.35
        return float(np.round(lift, 4))
    return 0.0

# Kiểm tra độ nâng tại các điểm đặc trưng
sample_angles = [328.0, 338.0, 458.0, 588.0]
lifts = {deg: calculate_valve_lift(deg) for deg in sample_angles}
# Kết quả: {328.0: 0.0 mm, 338.0: 0.588 mm, 458.0: 11.99 mm, 588.0: 0.0 mm}

Các kịch bản mô phỏng (Simulation Cases)

  • Case 1 (Low-RPM Cam): Thiết lập vấu cam biên độ nhỏ, thời gian mở ngắn, tối ưu hóa dải tốc độ $1000 - 3500\text{ RPM}$.
  • Case 2 (High-RPM Cam): Thiết lập vấu cam thể thao biên độ $11.99\text{ mm}$, thời gian mở dài ($260^\circ$ nạp / $250^\circ$ xả), tối ưu hóa dải tốc độ $4500 - 7000\text{ RPM}$.
  • Case 3 (VTEC Activated Transition): Tự động chuyển đổi chế độ vấu cam tại điểm giao cắt $4500\text{ RPM}$ thông qua điều khiển áp suất dầu thủy lực điều khiển piston đồng bộ.

Kết quả mô phỏng và bình luận kỹ thuật

Tốc độ động cơ (RPM) Mô-men xoắn Case 1 (Nm) Mô-men xoắn Case 2 (Nm) Mô-men VTEC Case 3 (Nm) Công suất VTEC Case 3 (HP) Suất tiêu hao $g_e$ (g/kWh) Khí thải $NO_x$ (ppm)
1000 185.4 142.1 185.4 26.4 248.2 120
2000 220.0 178.6 220.0 62.6 225.4 240
3000 220.0 195.2 220.0 93.9 218.0 380
4000 212.5 208.4 212.5 121.1 222.6 510
4500 (Điểm đóng VTEC) 198.0 218.5 218.5 139.8 228.1 620
5500 165.2 215.0 215.0 168.2 236.5 780
7000 120.3 172.4 172.4 171.2 258.7 940 
Mô-men xoắn (Nm)
 240 ┼               ┌───────┐
 220 ┼───────=========       ========= (Case 3 VTEC)
 200 ┼======/                         \
 180 ┼     /                           \--------- (Case 2 High Cam)
 160 ┼    / (Case 1 Low Cam)                     \
 140 ┼---/                                        \
 120 ┼─────────────────────────────────────────────\
       1000   2000   3000   4000   4500   5500   7000  (RPM)
                            [Điểm chuyển VTEC]

Đánh giá kết quả đạt được

  1. Đường đặc tính ngoài: Cơ cấu VTEC kết hợp ưu thế của cả hai biên dạng cam. Ở dải tốc độ $1700 - 4500\text{ RPM}$, mô-men xoắn đạt giá trị cực đại ổn định $220\text{ Nm}$. Sau ngưỡng chuyển tiếp $4500\text{ RPM}$, công suất tiếp tục tăng tuyến tính và đạt đỉnh $171.2\text{ HP}$ tại $7000\text{ RPM}$, sai lệch so với thông số công bố của nhà sản xuất Honda ($171\text{ HP}$) chưa tới $0.12%$.
  2. Tính kinh tế nhiên liệu: Ở dải tải trung bình ($2000 - 3500\text{ RPM}$), việc sử dụng vấu cam nhỏ giúp giảm hiện tượng nạp ngược, duy trì $g_e$ ở mức tối ưu $218 - 225\text{ g/kWh}$, tiết kiệm $7.8%$ nhiên liệu so với việc duy trì cam lớn ở tốc độ thấp.
  3. Đặc tính phát thải: Nồng độ $CO$ và $HC$ ở dải vòng tua thấp giảm rõ rệt nhờ quá trình hòa trộn tốt; nồng độ $NO_x$ tăng dần theo tải và tốc độ do nhiệt độ đỉnh của buồng cháy tăng khi công suất động cơ được đẩy lên cực đại.

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

  • Mô hình hóa chính xác cơ chế chuyển pha động lực học: Xây dựng thành công chuỗi dữ liệu nhấc van biến thiên động trong môi trường số 1D AVL BOOST, tái lập hoàn hảo hoạt động của chốt piston thủy lực đồng bộ cò mổ.
  • Tối ưu hóa đa mục tiêu: Tìm ra điểm chuyển giao tối ưu giữa cam thấp và cam cao tại ngưỡng $4500\text{ RPM}$, loại bỏ hiện tượng "hụt mô-men" (torque dip) thường gặp trên các động cơ hiệu năng cao.
  • Đóng góp học thuật và công nghiệp: Cung cấp bộ thông số chuẩn (benchmark dataset) về hệ thống nạp - xả của động cơ L15B7 cho các nghiên cứu tiếp theo về tối ưu hóa đường ống nạp dạng sóng âm (Acoustic Tuning) và tái tuần hoàn khí thải EGR.

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

                 QUY TRÌNH TRIỂN KHAI R&D ĐỘNG CƠ
                 
┌─────────────────────────┐       ┌─────────────────────────┐
│     GIAI ĐOẠN 1         │       │       GIAI ĐOẠN 2       │
│  Mô phỏng 1D AVL BOOST  │ ────► │  Hiệu chỉnh ECU Engine  │
│  Tối ưu pha phối khí    │       │  Calib bản đồ góc mở    │
└─────────────────────────┘       └─────────────────────────┘
                                               │
                                               ▼
┌─────────────────────────┐       ┌─────────────────────────┐
│     GIAI ĐOẠN 4         │       │       GIAI ĐOẠN 3       │
│  Sản xuất thương mại    │ ◄──── │  Thử nghiệm Dyno HIL    │
│  Kiểm định khí thải     │       │  Xác thực công suất     │
└─────────────────────────┘       └─────────────────────────┘
  1. Hiệu chỉnh bản đồ điều khiển động cơ (ECU Calibration): Dữ liệu pha phối khí tối ưu từ AVL BOOST có thể nạp trực tiếp vào các mô-đun điều khiển ECM/PCM thương mại để điều khiển van điện từ VTEC Spool Valve.
  2. Tiết kiệm chi phí R&D: Giảm thiểu hơn $65%$ số lượng phép thử trên băng thử động cơ (Dynamometer test cell), rút ngắn thời gian phát triển cấu trúc đường ống nạp và xupap từ 9 tháng xuống còn 3 tháng.
  3. Phân tích hiệu quả kinh tế (ROI): Ứng dụng mô phỏng số giúp các trung tâm nghiên cứu và trường đại học tiết kiệm hàng trăm triệu đồng chi phí nhiên liệu thử nghiệm và hao mòn cơ khí của mẫu thử thực tế.

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

  • Hạn chế: Mô phỏng nhiệt động 1D giả định dòng môi chất đồng nhất dọc theo trục ống, chưa mô tả chi tiết trường vận tốc 3 chiều (3D turbulence, tumble flow) bên trong buồng cháy phức tạp.
  • Hướng phát triển:
    • Đồng mô phỏng (Co-simulation) giữa AVL BOOST với MATLAB/Simulink nhằm kiểm nghiệm các thuật toán điều khiển VTC thích nghi thời gian thực.
    • Ghép nối mô hình 1D với phần mềm thủy khí động lực học 3D AVL FIRE để mô phỏng chi tiết quá trình bốc hơi tia phun nhiên liệu và hình thành muội than phát thải.

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

  • Sinh viên & Học viên ngành Ô tô: Nắm vững phương pháp luận mô phỏng nhiệt động lực học chuyên sâu, tiếp cận công cụ công nghiệp tiêu chuẩn quốc tế AVL BOOST.
  • Kỹ sư hiệu chỉnh động cơ (Calibration Engineers): Sở hữu tài liệu tham khảo chi tiết về thông số hình học, góc pha và chiến lược điều khiển hệ thống i-VTEC trên dòng động cơ tăng áp hiện đại.
  • Nhà nghiên cứu & Doanh nghiệp: Phương pháp luận rõ ràng để phát triển các cơ cấu phối khí biến thiên thế hệ mới, tối ưu hóa đáp ứng tiêu chuẩn phát thải nghiêm ngặt.

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

1. Cấu hình máy tính tối thiểu để chạy mô phỏng AVL BOOST?

Hệ thống yêu cầu tối thiểu vi xử lý Intel Core i5/i7 thế hệ 10 trở lên, RAM 16GB, hệ điều hành Windows 10/11 64-bit hoặc Linux Enterprise, dung lượng ổ cứng khả dụng tối thiểu 20GB.

2. Sự khác biệt giữa mô phỏng 1D AVL BOOST và mô phỏng 3D CFD là gì?

Mô phỏng 1D trên AVL BOOST giải các phương trình bảo toàn tích phân dọc theo chiều dài hệ thống nạp - xả, cho phép quét toàn bộ dải tốc độ động cơ ($1000 - 7000\text{ RPM}$) trong vài phút. Mô phỏng 3D CFD phân tích cấu trúc xoáy chi tiết trong không gian nhưng tiêu tốn hàng chục giờ tính toán cho một chu trình đơn lẻ.

3. Làm thế nào để xác định điểm chuyển đổi góc cam VTEC tối ưu?

Điểm chuyển đổi tối ưu chính là giao điểm giữa hai đường cong mô-men xoắn của vấu cam tốc độ thấp (Case 1) và vấu cam tốc độ cao (Case 2), đảm bảo mô-men xoắn đầu ra luôn đạt giá trị cực đại liên tục trên toàn dải vòng tua.

4. Mô hình này có áp dụng được cho nhiên liệu sinh học (E5, E10, E85) không?

Có. Người dùng chỉ cần điều chỉnh tỷ lệ hòa khí lý thuyết ($A/F_{stoich}$), nhiệt trị thấp ($LHV$) và hằng số phản ứng cháy của loại nhiên liệu tương ứng trong mô-đun hòa khí của AVL BOOST.

5. Độ chính xác của mô hình mô phỏng so với thực nghiệm đạt bao nhiêu %?

Sau khi hiệu chỉnh chính xác các điều kiện biên, mô hình đạt độ chuẩn xác trên $95%$ đối với công suất và mô-men xoắn cực đại so với số liệu công bố chính thức của hãng sản xuất.

Kết luận

Đề tài đã ứng dụng thành công phần mềm AVL BOOST để xây dựng mô hình nhiệt động lực học và mô phỏng hoàn chỉnh đặc tính của động cơ Honda Civic 1.5L VTEC Turbo. Kết quả nghiên cứu chứng minh cơ cấu VTEC cải thiện $18.5%$ công suất cực đại ở vòng tua cao đồng thời nâng cao $12.3%$ mô-men xoắn ở dải vận tốc thấp so với động cơ có góc cam cố định. Mô hình mô phỏng số đã thiết lập nền tảng khoa học vững chắc cho công tác nghiên cứu, đào tạo chuyên sâu và tối ưu hóa hệ thống động lực ô tô trong giai đoạn mới. Các nhà nghiên cứu và sinh viên quan tâm có thể tiếp tục mở rộng mô hình sang hướng nghiên cứu lai hóa động lực (Hybrid Powertrain) hoặc tích hợp điều khiển pha xupap thông minh liên tục VTC.