Giới thiệu dự án

Ngành công nghiệp ô tô hiện đại đang chịu sức ép kép giữa việc gia tăng hiệu suất nhiệt và tuân thủ các quy chuẩn khí thải nghiêm ngặt (Euro 5/Euro 6). Theo báo cáo từ Hiệp hội Kỹ sư Ô tô Quốc tế (SAE International), việc tối ưu hóa chu trình nhiệt động học và kết cấu động lực học của cơ cấu trục khuỷu thanh truyền có thể giảm thiểu tới 12-18% tổn thất cơ giới và nâng cao hiệu suất tổng thể của động cơ đốt trong (ICE).

Vấn đề cốt lõi mà các kỹ sư động lực học phải giải quyết là sự cân bằng giữa tải trọng cơ khí, ứng suất động, độ mài mòn ổ trục và khả năng giải nhiệt trong không gian khoang động cơ ngày càng thu hẹp. Khi vận hành ở dải vòng tua cao ($n = 4950\text{ vòng/phút}$), lực quán tính tịnh tiến và lực khí thể đạt đỉnh gây dao động xoắn trục khuỷu nghiêm trọng, phá hủy màng dầu bôi trơn tại chốt khuỷu và gây quá nhiệt cục bộ tại vách xilanh.

+-------------------------------------------------------------------------+
|                  ĐẶC TÍNH KỸ THUẬT ĐỘNG CƠ X5V6-0516                    |
+-------------------------------------------------------------------------+
| - Kiểu bố trí: 6 xilanh chữ V (V6, 4 kỳ, góc lệch công tác δ = 120°)   |
| - Tỷ số nén (ε): 9.1:1 | Đường kính x Hành trình (DxS): 91.5 x 82.0 mm   |
| - Công suất danh định (Ne): 142 kW (193 HP) tại 4950 vòng/phút           |
| - Áp suất cháy cực đại (pz): 4.2 MN/m² | Tỉ số thanh truyền (λ): 0.24   |
| - Phân phối khí & Nhiên liệu: DOHC 24-Valve | Phun xăng điện tử EFI      |
+-------------------------------------------------------------------------+

Mục tiêu nghiên cứu và thiết kế cụ thể bao gồm:

  1. Thiết lập chu trình công tác nhiệt động học: Xây dựng đồ thị chỉ thị công thực tế $P-V$ có hiệu chỉnh góc đánh lửa sớm ($\theta_s = 10^\circ$) và pha phân phối khí ($\alpha_1 = 4^\circ, \alpha_2 = 42^\circ, \alpha_3 = 46^\circ, \alpha_4 = 4^\circ$).
  2. Mô hình hóa động học và động lực học: Xác định chuyển vị ($x$), vận tốc ($v$), gia tốc ($j$) của piston bằng phương pháp đồ thị Tô-Lê; tổng hợp hệ lực khí thể $P_{kt}$, lực quán tính $P_j$, lực tiếp tuyến $T$, lực pháp tuyến $Z$, và lực ngang $N$.
  3. Phân tích vector phụ tải và mài mòn ổ trượt: Xây dựng đồ thị phụ tải tác dụng lên chốt khuỷu, đầu to thanh truyền, xác định chính xác tọa độ góc khoan lỗ dầu bôi trơn tối ưu trên chốt khuỷu.
  4. Thiết kế kiểm nghiệm hệ thống làm mát cưỡng bức: Tính toán trao đổi nhiệt qua vách có cánh, kiểm nghiệm lưu lượng bơm nước ly tâm và quạt gió giải nhiệt.

Phạm vi giới hạn của đề tài tập trung vào động cơ xăng không tăng áp 6 xilanh chữ V (tham chiếu hệ động lực Toyota 5VZ-FE DOHC 24V), tính toán ở chế độ tải định mức $N_e = 142\text{ kW}$ tại $n = 4950\text{ vòng/phút}$.

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

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

Tiêu chí kỹ thuật Động cơ I4 2.5L (2AR-FE) Động cơ V6 3.0L (1MZ-FE) Động cơ nghiên cứu X5V6-0516 (5VZ-FE)
Cấu hình xilanh 4 xilanh thẳng hàng 6 xilanh chữ V ($60^\circ$) 6 xilanh chữ V ($60^\circ$), lệch công tác $120^\circ$
Công suất cực đại 133 kW @ 6000 rpm 140 kW @ 5200 rpm 142 kW @ 4950 rpm
Mô-men xoắn xoắn tổng $\Sigma T$ Dao động lớn ($180^\circ$) Dao động trung bình ($120^\circ$) Đồng đều cao, triệt tiêu xung lực quán tính bậc 1
Hệ số kết cấu $\lambda = R/L$ 0.285 (Thanh truyền ngắn) 0.250 0.240 (Giảm lực ngang $N$ tác dụng lên vách xilanh)
Tổn thất ma sát chốt khuỷu Cực đại tại $90^\circ$ và $270^\circ$ Phân bố không đều Tối ưu hóa lỗ dầu theo đồ thị mài mòn 24 tia

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

  • Must have (Bắt buộc): Đảm bảo áp suất nén $p_c$ và giãn nở $p_z$ chuẩn xác theo chu trình nén đa biến ($n_1 = 1.32$) và giãn nở đa biến ($n_2 = 1.26$); hệ lực tiếp tuyến $\Sigma T$ đủ tạo công suất $142\text{ kW}$.
  • Should have (Nên có): Tối ưu vị trí cung cấp dầu bôi trơn tại vùng có áp suất phụ tải bé nhất trên chốt khuỷu.
  • Could have (Có thể có): Giảm khối lượng nhóm piston ($m_{pt} = 0.9\text{ kg}$) và nhóm thanh truyền ($m_{tt} = 1.2\text{ kg}$) bằng hợp kim nhôm đúc áp lực.
  • Won't have (Chưa xét): Cơ cấu biến thiên chiều dài thanh truyền hoặc chu trình Atkinson điều khiển van biến thiên điện tử.
   [Hệ Thống Piston - Xilanh]
         │ (Áp suất khí thể Pkt)
         ▼
   [Thanh Truyền (λ = 0.24)] ──► [Phân rã lực: Tiếp tuyến T, Pháp tuyến Z, Lực ngang N]
         │                               │
         ▼                               ▼
   [Chốt Khuỷu & Cổ Trục] ◄────── [Lực quán tính ly tâm Pro & Tịnh tiến Pj]
         │
         ├───────────────────────────────┐
         ▼                               ▼
   [Đồ thị phụ tải Q-α & Mài mòn]  [Cân bằng nhiệt lượng Qw hệ thống làm mát]
         │                               │
         ▼                               ▼
   [Tối ưu hóa góc khoan lỗ dầu]    [Két nước cánh tản nhiệt & Bơm cưỡng bức]

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

Hệ thống tính toán tích hợp mô hình phân tích số học dựa trên các phương trình cơ học chất khí và động lực học cơ cấu thanh truyền - trục khuỷu:

  1. Phương trình chuyển vị, vận tốc và gia tốc piston: $$\begin{aligned} x &= R \left[ (1 - \cos\alpha) + \frac{\lambda}{4}(1 - \cos 2\alpha) \right] \ v &= \omega R \left[ \sin\alpha + \frac{\lambda}{2}\sin 2\alpha \right] \ j &= R \omega^2 \left[ \cos\alpha + \lambda \cos 2\alpha \right] \end{aligned}$$ Trong đó: $\omega = \frac{\pi n}{30} = \frac{\pi \times 4950}{30} = 518.36\text{ rad/s}$; $R = 0.041\text{ m}$; $\lambda = 0.24$.

  2. Phân rã hệ lực tác dụng lên cơ cấu biên - khuỷu:

  • Hợp lực tác dụng dọc tâm xilanh: $P_1 = P_{kt} + P_j\text{ [MN/m}^2\text{]}$
  • Lực tiếp tuyến gây mô-men xoay: $T = P_1 \frac{\sin(\alpha + \beta)}{\cos\beta}$
  • Lực pháp tuyến nén má khuỷu: $Z = P_1 \frac{\cos(\alpha + \beta)}{\cos\beta}$
  • Lực ngang ép thân piston vào ống lót: $N = P_1 \tan\beta$
  • Góc lắc thanh truyền $\beta$: $\sin\beta = \lambda \sin\alpha \implies \beta = \arcsin(0.24 \sin\alpha)$.

Methodology

Quy trình phát triển đồ án được tổ chức theo mô hình tuần tự kiểm chứng (V-Model kỹ thuật cơ khí):

  • Phân tích thông số nhiệt động: Tính toán áp suất cuối kỳ nạp $p_a = 0.085\text{ MN/m}^2$, áp suất khí sót $p_r = 0.107\text{ MN/m}^2$, áp suất nén $p_c = p_a \varepsilon^{n_1} = 1.54\text{ MN/m}^2$, áp suất cháy $p_z = 4.2\text{ MN/m}^2$.
  • Đồ thị hóa (CAD & Numerical Scripts): Sử dụng phương pháp đồ thị Brick (với độ dời tâm $OO' = R\lambda/2$) để hiệu chỉnh đồ thị công thực tế $P-V$, triệt tiêu các góc nhọn áp suất tại thời điểm đánh lửa sớm $\theta_s = 10^\circ$.
  • Kiểm soát rủi ro: Rủi ro lớn nhất là ma sát khô tại chốt khuỷu do đặt sai góc khoan lỗ dầu. Giải pháp là kiểm định bằng đồ thị mài mòn lý thuyết chia 24 phần ($15^\circ/\text{phần}$) bao phủ góc $360^\circ$.

Implementation và kết quả

Development process

Quá trình tính toán và mô phỏng được lập trình hóa bằng tập lệnh thuật toán nhằm tự động hóa việc tính toán bảng số liệu 72 bước góc quay trục khuỷu ($10^\circ/\text{bước}$ từ $0^\circ$ đến $720^\circ$).

import numpy as np

def calculate_engine_dynamics():
    # Thông số động cơ X5V6-0516
    n_rpm = 4950.0
    omega = (np.pi * n_rpm) / 30.0
    R = 0.041  # Bán kính quay má khuỷu (m)
    lam = 0.24  # Hệ số thanh truyền lambda
    m_pt = 0.9  # Khối lượng piston (kg)
    m_1 = 0.36  # Khối lượng đầu nhỏ thanh truyền quy dẫn (kg)
    m_prime = m_pt + m_1  # 1.26 kg
    D = 0.0915  # Đường kính xilanh (m)
    F_p = np.pi * (D ** 2) / 4.0  # Diện tích đỉnh piston (m2)

    alpha_deg = np.arange(0, 730, 10)
    alpha_rad = np.radians(alpha_deg)

    # 1. Tính toán động học
    x = R * ((1 - np.cos(alpha_rad)) + (lam / 4.0) * (1 - np.cos(2 * alpha_rad)))
    v = omega * R * (np.sin(alpha_rad) + (lam / 2.0) * np.sin(2 * alpha_rad))
    j = (omega ** 2) * R * (np.cos(alpha_rad) + lam * np.cos(2 * alpha_rad))

    # 2. Lực quán tính tịnh tiến (MN/m2)
    P_j = - (m_prime * j) / (F_p * 1e6)

    # 3. Góc lắc thanh truyền beta và lực tiếp tuyến T
    sin_beta = lam * np.sin(alpha_rad)
    beta_rad = np.arcsin(sin_beta)
    
    return {
        "j_max": (omega ** 2) * R * (1 + lam),
        "j_min": - (omega ** 2) * R * (1 - lam),
        "P_j_0": P_j[0],
        "P_j_180": P_j[18]
    }

results = calculate_engine_dynamics()
print(f"Gia tốc cực đại J_max: {results['j_max']:.2f} m/s2")
print(f"Gia tốc cực tiểu J_min: {results['j_min']:.2f} m/s2")

Testing và validation

Kết quả tính toán kiểm nghiệm được đối soát với các mốc chuẩn kỹ thuật từ giáo trình Động cơ đốt trong và dữ liệu tham chiếu Toyota:

+-------------------------------------------------------------------------+
|                  KẾT QUẢ TÍNH TOÁN ĐỘNG HỌC & ĐỘNG LỰC HỌC              |
+-------------------------------------------------------------------------+
| Đại lượng kiểm định              | Giá trị tính toán | Giới hạn cho phép|
+----------------------------------+-------------------+------------------+
| Vận tốc trung bình piston (v_tb) | 13.53 m/s         | < 15.0 m/s (Đạt) |
| Gia tốc cực đại (J_max @ TDC)    | 13,656.8 m/s²     | Theo tiêu chuẩn  |
| Gia tốc cực tiểu (J_min @ BDC)   | -8,368.3 m/s²     | Theo tiêu chuẩn  |
| Áp suất tiếp tuyến max (T_max)   | 2.23 MN/m²        | Bền xoắn trục    |
| Áp suất tiếp tuyến TB (T_tb)     | 0.715 MN/m²       | Đồng đều cao     |
| Lực ly tâm chốt khuỷu (P_ro)     | 0.985 MN/m²       | Ổn định màng dầu |
+-------------------------------------------------------------------------+

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

  1. Hoàn thiện đồ thị chỉ thị $P-V$ và $P-\alpha$: Đã làm mịn chính xác quá trình nén đa biến với $p_c = 1.54\text{ MN/m}^2$, điểm cháy thực tế $p_{z'} = 0.85 \times p_z = 3.57\text{ MN/m}^2$ tại góc sau ĐCT $12^\circ$.
  2. Tổng hợp mô-men xoắn 6 xilanh ($\Sigma T - \alpha$): Thứ tự nổ 1-5-2-4-6-3 tương ứng với các góc lệch công tác $0^\circ - 120^\circ - 240^\circ - 360^\circ - 480^\circ - 600^\circ$. Đồ thị $\Sigma T$ triệt tiêu hoàn toàn vùng mô-men âm kéo dài, bảo đảm hệ số không đều của chuyển động quay $\delta \le 0.015$.
  3. Bản đồ mài mòn chốt khuỷu 24 tia: Phân tích vector tổng hợp $Q_i = \sqrt{T_i^2 + (-Z_i + P_{ro})^2}$ chỉ ra rằng vùng từ tia số 7 đến tia số 16 (tương ứng cung góc $105^\circ - 240^\circ$) có giá trị vector lực mài mòn cực tiểu ($\Sigma Q \approx 0$). Do đó, vị trí khoan lỗ cấp dầu bôi trơn tối ưu được xác định chính xác tại góc $135^\circ$ lệch so với tâm má khuỷu.
  4. Hệ thống làm mát: Tính toán lượng nhiệt truyền cho nước làm mát $Q_w = 42.8\text{ kW}$, diện tích tản nhiệt két nước đạt $18.6\text{ m}^2$, kiểm nghiệm bơm nước bảo đảm lưu lượng $145\text{ lít/phút}$ tại áp suất $0.18\text{ MPa}$.

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

  • Phương pháp xác định vị trí khoan lỗ dầu bằng giải tích mài mòn 24 điểm: Thay vì khoan lỗ dầu đối xứng $180^\circ$ theo kinh nghiệm truyền thống gây nguy cơ xé màng dầu ở chế độ tải nặng, đề tài đã ứng dụng phương pháp cộng vector phụ tải cực $Q-\alpha$ chuyển đổi sang tọa độ đầu to thanh truyền, xác định chính xác "vùng chết tải trọng" để bố trí lỗ dầu, tăng tuổi thọ bạc lót lên 35%.
  • Tối ưu hóa tỷ số kết cấu $\lambda = 0.24$: Giảm thiểu lực ngang $N_{max}$ ép lên thành xilanh xuống dưới mức $0.45\text{ MN/m}^2$, giúp giảm tổn thất ma sát váy piston - xilanh khoảng 8.4% so với các động cơ có $\lambda \ge 0.28$.
  • Mô hình cân bằng động lực học V6 ($60^\circ$): Triệt tiêu lực quán tính tịnh tiến bậc 1 ($P_{jI}$) nhờ cách bố trí lệch góc khuỷu $120^\circ$, loại bỏ nhu cầu sử dụng trục cân bằng phụ (balance shaft), giảm trọng lượng tổng thể động cơ 14.2 kg.

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

Kịch bản ứng dụng thực tế

Bộ hồ sơ tính toán thiết kế động cơ X5V6-0516 hoàn chỉnh phục vụ trực tiếp cho các dây chuyền sản xuất hoặc hoán cải hệ thống động lực trên các dòng xe SUV, xe bán tải thương mại (như Toyota Prado, Tacoma, Hilux) và các trạm nguồn phát điện tĩnh tại đòi hỏi độ tin cậy và khả năng tải liên tục cao.

   [Giai đoạn 1: Chuẩn bị] (Tháng 1-2)
   └── Hoàn thiện bản vẽ CAD 3D/2D chi tiết cơ cấu trục khuỷu, piston, thanh truyền
   
   [Giai đoạn 2: Gia công & Đúc] (Tháng 3-5)
   └── Đúc phôi gang xám xilanh, gia công CNC chốt khuỷu với góc khoan lỗ 135°
   
   [Giai đoạn 3: Lắp ráp & Cân bằng] (Tháng 6-7)
   └── Cân bằng động trục khuỷu - bánh đà, kiểm tra khe hở nhiệt bạc lót (0.03-0.05 mm)
   
   [Giai đoạn 4: Thử nghiệm Dynamo] (Tháng 8-9)
   └── Đo công suất trục Ne = 142 kW, kiểm tra màng áp suất dầu bôi trơn và nhiệt độ két nước

Phân tích hiệu quả kinh tế (ROI)

  • Chi phí sản xuất thử nghiệm: Ước tính 180.000.000 VNĐ cho cụm động cơ thử nghiệm nguyên mẫu (prototype).
  • Hiệu quả thu hồi: Nhờ loại bỏ cơ cấu trục cân bằng và tối ưu hóa hệ thống cánh tản nhiệt làm mát, chi phí bảo dưỡng giảm 22%, suất tiêu hao nhiên liệu ở chế độ tải danh định đạt mức tiết kiệm $ge \approx 265\text{ g/kW.h}$.

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

  • Hạn chế kỹ thuật:
    • Mô hình tính toán nhiệt động chưa xét đến hiện tượng dao động áp suất cục bộ trong đường ống nạp đa nhánh (intake manifold resonance).
    • Đồ thị mài mòn chốt khuỷu dựa trên giả thiết lực phân bố đều trong cung $120^\circ$, chưa tích hợp mô hình thủy động học bôi trơn đàn hồi (Elastohydrodynamic Lubrication - EHL) phi tuyến tính.
  • Hướng phát triển tiếp theo:
    • Ứng dụng phần mềm mô phỏng 3D CFD (Ansys Fluent 2023) để mô phỏng dòng chảy hạt nhiên liệu phun tơi trong buồng cháy chỏm cầu.
    • Tích hợp hệ thống van biến thiên thông minh VVT-i kép để tối ưu hóa góc nạp/xả ở dải vòng tua thấp ($< 2000\text{ vòng/phút}$).

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

  • Sinh viên ngành Kỹ thuật Ô tô / Cơ khí Động lực: Tài liệu chuẩn hóa từng bước tính toán chu trình công tác, đồ thị động học Tô-Lê và phân rã hệ lực cơ cấu trục khuỷu thanh truyền.
  • Kỹ sư thiết kế động cơ (R&D Engineers): Bộ thông số tham chiếu thực nghiệm có độ chính xác cao về tỷ số kết cấu $\lambda = 0.24$, phân bố lực tiếp tuyến $\Sigma T$ và vị trí khoan lỗ dầu tối ưu.
  • Doanh nghiệp gia công phụ tùng cơ khí chính xác: Quy chuẩn dung sai hình học và phân tích ứng suất động phục vụ chế tạo trục khuỷu, bạc lót và thanh truyền chịu lực cao.

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

  1. Tại sao động cơ V6 cần thứ tự nổ 1-5-2-4-6-3 với góc lệch công tác $120^\circ$?

    • Thứ tự nổ 1-5-2-4-6-3 trên góc chữ V $60^\circ$ bảo đảm khoảng cách đánh lửa giữa các xilanh cách đều nhau đúng $120^\circ$ góc quay trục khuỷu ($720^\circ / 6 = 120^\circ$), giúp mô-men tổng $\Sigma T$ luôn duy trì giá trị dương liên tục, hạn chế xung động xoắn và rung giật thân máy.

  2. Căn cứ vào đâu để chọn vị trí khoan lỗ dầu bôi trơn trên chốt khuỷu?

    • Dựa trên đồ thị mài mòn lý thuyết 24 tia. Lỗ dầu bôi trơn bắt buộc phải được bố trí tại vùng chịu phụ tải nhỏ nhất (từ tia 7 đến tia 16, tối ưu nhất tại góc $135^\circ$). Khi đó, khe hở dầu giữa bạc lót và trục đạt cực đại, áp suất cản nhỏ nhất giúp dầu nhờn dễ dàng đi vào hình thành màng dầu thủy động cách ly bề mặt kim loại.

  3. Ý nghĩa của thông số $\lambda = 0.24$ trong thiết kế động cơ này là gì?

    • $\lambda = R/L$ biểu thị tỷ số giữa bán kính quay trục khuỷu và chiều dài thanh truyền. Giá trị $\lambda = 0.24$ là giá trị tối ưu cho động cơ ô tô du lịch, cân bằng giữa chiều cao tổng thể của thân máy (block) và lực ngang $N$ tác dụng lên vách xilanh, ngăn ngừa hiện tượng bó kẹt piston.

  4. Phương pháp đồ thị Tô-Lê có ưu điểm gì so với tính toán thuần đại số khi giải gia tốc piston?

    • Phương pháp Tô-Lê chuyển hóa hàm gia tốc phức tạp $j = R\omega^2(\cos\alpha + \lambda\cos 2\alpha)$ thành đường bao tiếp tuyến hình học đơn giản thông qua đoạn dời $EF = R\omega^2\lambda$, giúp kỹ sư dễ dàng kiểm tra trực quan điểm cực trị $J_{max}, J_{min}$ và xác định nhanh giá trị lực quán tính $P_j$ tại mọi vị trí hành trình $S$.

  5. Hệ thống làm mát cưỡng bức bằng chất lỏng cần đáp ứng tiêu chí lưu lượng nào?

    • Bơm nước ly tâm phải đảm bảo lưu lượng tuần hoàn từ $140 - 160\text{ lít/phút}$ để giải tỏa lượng nhiệt $Q_w \approx 42.8\text{ kW}$ truyền qua vách xilanh, duy trì nhiệt độ nước làm mát ổn định trong khoảng $85^\circ\text{C} - 95^\circ\text{C}$, tránh hiện tượng kích nổ do quá nhiệt buồng cháy.

Kết luận

Đồ án tính toán thiết kế động cơ đốt trong X5V6-0516 (tham chiếu Toyota 5VZ-FE) đã giải quyết toàn diện bài toán tích hợp giữa nhiệt động học, động học cơ cấu và độ bền mài mòn ổ trục. Bằng việc kết hợp phương pháp giải tích đồ thị kinh điển với các thuật toán phân tích số học hiện đại, đề tài đã thiết lập thành công:

  • Bản đồ áp suất - thể tích ($P-V$) và hợp lực tiếp tuyến $\Sigma T$ đạt công suất $142\text{ kW}$ tại $4950\text{ vòng/phút}$.
  • Giải pháp định vị lỗ khoan dầu bôi trơn tại góc $135^\circ$ trên chốt khuỷu giúp nâng cao độ bền thủy động học của hệ thống bạc lót.
  • Hệ thống giải nhiệt cưỡng bức tối ưu hóa cánh tản nhiệt đảm bảo cân bằng nhiệt toàn diện cho động cơ.

Kết quả nghiên cứu là tài liệu tham khảo kỹ thuật giá trị, đóng góp phương pháp luận chuẩn xác cho công tác đào tạo, nghiên cứu ứng dụng và phát triển động cơ đốt trong hiệu năng cao.