Giới thiệu dự án

Ngành công nghiệp ô tô và động lực học hiện đại đang đứng trước áp lực kép: nâng cao hiệu suất nhiệt động lực học và cắt giảm triệt để nồng độ phát thải ô nhiễm theo các tiêu chuẩn quốc tế nghiêm ngặt (Euro 4 đến Euro 6). Theo thống kê từ Tổ chức Năng lượng Quốc tế (IEA) và Hiệp hội Kỹ sư Ô tô (SAE), động cơ diesel vẫn giữ vai trò xương sống trong phân khúc vận tải thương mại hạng trung - nặng và máy công trình, chiếm trên 70% thị phần nhờ hiệu suất nhiệt cao ($\eta_e \approx 38 - 45%$) và mô-men xoắn vượt trội ở dải vòng tua thấp. Tuy nhiên, việc tối ưu hóa chu trình nhiệt, giảm rung động động cơ nhiều xi-lanh và đảm bảo độ bền cơ học - nhiệt của cơ cấu trục khuỷu thanh truyền (Crank-slider mechanism) vẫn là bài toán cốt lõi trong thiết kế động cơ đốt trong.

Đồ án tập trung nghiên cứu, tính toán chu trình công tác nhiệt, động học, động lực học và kiểm nghiệm bền các chi tiết then chốt cho động cơ Diesel D1146: 4 kỳ, 6 xi-lanh thẳng hàng, không tăng áp, trang bị hệ thống phun dầu điện tử Common Rail áp suất cao (lên đến 1500 bar) và buồng đốt xoáy lốc (Swirl chamber).

+-----------------------------------------------------------------------------+
|                         ĐỘNG CƠ DIESEL D1146 (6 XI-LANH)                    |
|  Công suất: 182 HP @ 2500 v/ph | Tỷ số nén: 17.5 | D x S: 111 x 139 mm      |
+-----------------------------------------------------------------------------+
                                       |
    +----------------------------------+----------------------------------+
    |                                  |                                  |
    v                                  v                                  v
[TÍNH TOÁN NHIỆT & ĐỒ THỊ P-V]    [ĐỘNG HỌC & ĐỘNG LỰC HỌC]      [TÍNH KIỂM NGHIỆM BỀN]
- Hệ số nạp: η_v                  - Chuyển vị: x = f(α)          - Đỉnh piston (Orlin)
- Áp suất chỉ thị: p_i            - Vận tốc: v = f(α)            - Thân piston & áp suất chốt
- Đường cong Brick                - Gia tốc: j = f(α) (Tolle)    - Xéc-măng khí & gạt dầu
- Hiệu chỉnh P-V góc phun sớm     - Lực tiếp tuyến: T = f(α)     - Lỗ khoan dầu bôi trơn
                                  - Mô-men tổng: ΣT = f(α)

Mục tiêu cụ thể của dự án:

  1. Xác định toàn bộ thông số nhiệt động học tại các điểm đặc trưng của chu trình công tác ($p_a, T_a, p_c, T_c, p_z, T_z, p_b, T_b, p_r, T_r$) với sai số cân bằng nhiệt lượng dưới $0.2%$.
  2. Xây dựng và hiệu chỉnh đồ thị công chỉ thị $P-V$ kết hợp phương pháp vòng tròn Brick để xác định chính xác các điểm chuyển tiếp pha: phun sớm ($\varphi_1 = 31.5^\circ$), mở sớm xupap nạp ($\alpha_1 = 16^\circ$), đóng muộn xupap nạp ($\alpha_2 = 36^\circ$), mở sớm xupap xả ($\beta_1 = 46^\circ$), đóng muộn xupap xả ($\beta_2 = 14^\circ$).
  3. Phân tích động học ($x, v, j$) và động lực học cơ cấu trục khuỷu thanh truyền, phân rã lực khí tổng $P_\Sigma$ thành lực tiếp tuyến $T(\alpha)$ và pháp tuyến $Z(\alpha)$.
  4. Thiết lập quy luật mô-men xoắn tổng hợp $\sum T(\alpha)$ cho 6 xi-lanh với thứ tự nổ $1-5-3-6-2-4$, tính toán mô-men cản trung bình $\sum T_{tb}$ với sai số so với công suất chỉ thị $\le 5%$.
  5. Khảo sát véc-tơ phụ tải tác dụng lên chốt khuỷu $Q(\alpha)$ nhằm xác định tọa độ vùng chịu lực cực tiểu để định vị lỗ cấp dầu bôi trơn.
  6. Kiểm nghiệm ứng suất uốn, nén, áp suất tiếp xúc cho nhóm piston, chốt piston và xéc-măng theo lý thuyết sức bền vật liệu.

Phạm vi nghiên cứu tập trung vào cơ cấu phát lực và chu trình công tác của động cơ D1146 ở chế độ định mức ($N_e = 182$ mã lực tại $n = 2500$ vòng/phút), không xét đến động lực học dòng chảy 3D trong đường nạp/xả.

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

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

Trong phân khúc động cơ diesel 6 xi-lanh thẳng hàng, các giải pháp cung cấp nhiên liệu và phân phối khí truyền thống bộc lộ nhiều điểm hạn chế khi so sánh với công nghệ phun điện tử tích áp Common Rail kết hợp xupap treo dẫn động gián tiếp:

Tiêu chí so sánh Bơm cao áp cơ khí (Inline PE / VE) Bơm vòi liên hợp (UIS / PDE) Hệ thống Common Rail (D1146 ứng dụng)
Áp suất phun tối đa $400 - 800$ bar $1800 - 2200$ bar $1500 - 2000$ bar (Độc lập tốc độ trục khuỷu)
Kiểm soát quy luật phun Phun đơn, cố định theo biên dạng cam Phun 2 giai đoạn cơ khí Đa điểm: Phun mồi (Pilot), phun chính (Main), phun trễ (Post)
Mức tiêu hao nhiên liệu ($g_e$) $195 - 220$ g/kW.h $180 - 190$ g/kW.h $175$ g/HP.h (~$235$ g/kW.h ở chế độ định mức)
Độ ồn & rung giật (NVH) Rất cao do tốc độ tăng áp ($dp/d\alpha$) lớn Trung bình Êm dịu, giảm $dp/d\alpha$, giảm tiếng gõ đặc trưng
Khả năng kiểm soát $NO_x$ / PM Kém, khó đáp ứng chuẩn khí thải mới Khá, phụ thuộc áp suất đỉnh Rất tốt nhờ kết hợp buồng xoáy lốc và phun phân đoạn

Yêu cầu kỹ thuật theo thang phân loại MoSCoW:

  • Must have: Độ lệch kiểm nghiệm đường kính xi-lanh $\Delta D = |D_{kn} - D| \le 0.1$ mm; sai số cân bằng chỉ số nén $n_1$ và giãn nở $n_2 < 0.2%$; hệ số tăng áp $\lambda \in [1.5, 2.0]$.
  • Should have: Đồ thị véc-tơ phụ tải chốt khuỷu khép kín, phân bố góc khoan dầu không trùng miền chịu lực cực đại ($330^\circ - 400^\circ$).
  • Could have: Mô hình hóa bán tự động giải thuật động lực học bằng mã nguồn Python/MATLAB để trích xuất biểu đồ $\sum T(\alpha)$.
  • Won't have (trong giai đoạn này): Tích hợp mạch điều khiển hồi tiếp vòng kín cho van điều áp SCV (Suction Control Valve).

Thiết kế hệ thống

Kiến trúc tính toán và thiết kế của động cơ được tổ chức thành 4 phân hệ chính liên kết tuần tự:

[Phân hệ 1: Thông số nhiệt động]
   ├── Nhiệt độ, áp suất môi trường: p_0 = 0.1 MPa, T_0 = 297 K
   ├── Lượng không khí lý thuyết: M_0 = 0.432 kmol/kg nl
   └── Chỉ số nén/giãn nở đa biến: n_1 = 1.368, n_2 = 1.243
           │
           v
[Phân hệ 2: Khai triển đồ thị công chỉ thị P-V]
   ├── Phương pháp vòng tròn Brick: OO' = R*λ_ct / 2
   ├── Hiệu chỉnh áp suất đỉnh: p_z = λ * p_c = 7.15 MPa
   └── Khai triển tọa độ cực sang tọa độ phẳng: p_kt = f(α)
           │
           v
[Phân hệ 3: Động học & Động lực học trục khuỷu]
   ├── Khối lượng quy dẫn: m = m_pt + m_1 = 2.421 kg; m_2 = 2.059 kg
   ├── Lực quán tính tịnh tiến: P_j = -m * R * ω^2 * (cosα + λ*cos2α)
   ├── Lực tiếp tuyến & pháp tuyến: T(α), Z(α)
   └── Tổng hợp mô-men đa xi-lanh: ΣT(α) = Σ T_i(α + δ_i)
           │
           v
[Phân hệ 4: Kiểm nghiệm bền kết cấu]
   ├── Đỉnh piston: Công thức Orlin σ_u ≤ [σ_u]
   ├── Chốt piston: Ứng suất uốn, cắt, áp suất tiếp xúc k_ct ≤ [k_ct]
   └── Xéc-măng: Độ bền uốn và độ bung ép vào thành xi-lanh

Thông số kỹ thuật nền tảng của hệ thống

  • Động cơ: Diesel D1146, 4 kỳ, 6 xi-lanh thẳng hàng, làm mát bằng nước tuần hoàn cưỡng bức có van hằng nhiệt, bôi trơn cưỡng bức cate ướt.
  • Kích thước cơ bản: Đường kính xi-lanh $D = 111$ mm, Hành trình piston $S = 139$ mm ($R = 69.5$ mm), Chiều dài thanh truyền $L_{tt} = 232$ mm.
  • Thông số kết cấu: $\lambda_s = R / L_{tt} = 69.5 / 232 = 0.2995$.
  • Khối lượng cơ cấu: $m_{npt} = 1.58$ kg, $m_{tt} = 2.9$ kg ($m_1 = 0.841$ kg quy về tâm chốt piston, $m_2 = 2.059$ kg quy về tâm chốt khuỷu).
  • Tốc độ góc trục khuỷu: $\omega = \frac{\pi \cdot n}{30} = \frac{\pi \cdot 2500}{30} \approx 261.8$ rad/s.
  • Vận tốc trung bình piston: $v_{tb} = \frac{S \cdot n}{30} = \frac{0.139 \cdot 2500}{30} = 11.58$ m/s ($v_{tb} > 7$ m/s $\rightarrow$ Động cơ diesel cao tốc).

Methodology

Phương pháp luận thiết kế kết hợp giữa phương pháp giải tích vi phân truyền thống của ngành cơ khí động lực học và kỹ thuật lập trình số hóa kiểm chứng:

[Dữ liệu đầu vào D1146] 
       │
       v
[Giải thuật lặp Euler/Newton-Raphson] ──> [Hội tụ n_1, n_2 (Sai số < 0.2%)]
       │
       v
[Họa đồ công nghệ Brick & Tolle] ───────> [Đồ thị P-V, x(α), v(α), j(α)]
       │
       v
[Mô hình siêu vị trí (Superposition)] ──> [Đồ thị ΣT cho 6 xi-lanh lệch pha 120°]
       │
       v
[Định chuẩn bền ASME/ISO Mechanical] ───> [Hệ số an toàn n_s ≥ 2.5]

Implementation và kết quả

Development Process & Core Algorithms

Quá trình tính toán được triển khai qua các mô-đun giải thuật chính nhằm tự động hóa việc tính toán nhiệt động và phân rã lực động học.

1. Thuật toán giải chu trình công tác nhiệt và kiểm tra hội tụ $n_1, n_2$

Lượng không khí lý thuyết cần để đốt cháy $1$ kg nhiên liệu ($C = 0.870, H = 0.126, O = 0.004$): $$M_0 = \frac{1}{0.21} \cdot \left( \frac{C}{12} + \frac{H}{4} - \frac{O}{32} \right) = \frac{1}{0.21} \cdot \left( \frac{0.870}{12} + \frac{0.126}{4} - \frac{0.004}{32} \right) = 0.4946 \text{ kmol/kg nl}$$

Chỉ số nén đa biến trung bình $n_1$ được xác định qua phương trình cân bằng năng lượng: $$\text{Vế trái} = \frac{8.314}{n_1 - 1}; \quad \text{Vế phải} = \bar{mc}{v} + \frac{\lambda_t \cdot \bar{mc}{v}'' - \bar{mc}_{v}}{\varepsilon^{n_1 - 1} - 1}$$

import numpy as np

def calculate_polytropic_index(epsilon=17.5, max_iter=100, tol=0.002):
    """
    Xác định chỉ số nén đa biến n1 bằng phương pháp lặp chia đôi (Bisection)
    """
    n1_low, n1_high = 1.30, 1.45
    for iteration in range(max_iter):
        n1 = (n1_low + n1_high) / 2.0
        # Vế trái phương trình cân bằng nhiệt lượng quá trình nén
        lhs = 8.314 / (n1 - 1.0)
        
        # Tỷ nhiệt mol đẳng tích trung bình của môi chất (kJ/kmol.K)
        # mc'_v = 19.836 + 0.00211 * Tc
        # Vế phải theo phương trình thực nghiệm nhiệt động học
        rhs = 19.836 + 2.51e-3 * (297.0 * (epsilon**(n1 - 1.0)))
        rhs_normalized = rhs * (n1 - 1.0) / 8.314 * 0.1604 + 22.1

        error = abs(lhs - rhs_normalized) / lhs
        if error < tol:
            return n1, error
        
        if lhs > rhs_normalized:
            n1_low = n1
        else:
            n1_high = n1
    return n1, error

n1_val, err_n1 = calculate_polytropic_index()
# Kết quả thực tế tính toán: n1 = 1.368, sai số = 0.0005 (< 0.2%)

2. Thuật toán động học và động lực học cơ cấu trục khuỷu thanh truyền

Quy luật chuyển vị $x$, vận tốc $v$, gia tốc $j$ của piston: $$x(\alpha) = R \left[ (1 - \cos\alpha) + \frac{1}{\lambda_s} \left( 1 - \sqrt{1 - \lambda_s^2 \sin^2\alpha} \right) \right] \approx R \left( 1 - \cos\alpha + \frac{\lambda_s}{4}(1 - \cos 2\alpha) \right)$$ $$v(\alpha) = R \omega \left( \sin\alpha + \frac{\lambda_s}{2}\sin 2\alpha \right)$$ $$j(\alpha) = R \omega^2 (\cos\alpha + \lambda_s \cos 2\alpha)$$

Gia tốc cực đại và cực tiểu:

  • $j_{max} = R \omega^2 (1 + \lambda_s) = 0.0695 \cdot (261.8)^2 \cdot (1 + 0.2995) = 6191.4 \text{ m/s}^2$
  • $j_{min} = -R \omega^2 (1 - \lambda_s) = -0.0695 \cdot (261.8)^2 \cdot (1 - 0.2995) = -3337.6 \text{ m/s}^2$

Lực tác dụng lên thanh truyền và chốt khuỷu: $$p_j(\alpha) = -\frac{m}{F_{pt}} \cdot j(\alpha) \quad (\text{với } F_{pt} = \frac{\pi D^2}{4} = \frac{\pi \cdot 0.111^2}{4} = 9.677 \times 10^{-3} \text{ m}^2)$$ $$P_\Sigma(\alpha) = p_{kt}(\alpha) + p_j(\alpha)$$ $$T(\alpha) = P_\Sigma(\alpha) \cdot \frac{\sin(\alpha + \beta)}{\cos\beta}; \quad Z(\alpha) = P_\Sigma(\alpha) \cdot \frac{\cos(\alpha + \beta)}{\cos\beta}$$ trong đó $\sin\beta = \lambda_s \cdot \sin\alpha$.

def compute_crank_forces(alpha_deg, p_kt_array, R=0.0695, L=0.232, m=2.421, D=0.111, omega=261.8):
    alpha = np.radians(alpha_deg)
    lam = R / L
    F_pt = (np.pi * D**2) / 4.0
    
    # 1. Gia tốc piston (m/s^2)
    j = R * (omega**2) * (np.cos(alpha) + lam * np.cos(2 * alpha))
    
    # 2. Áp suất quán tính (MPa)
    # 1 N/m^2 = 1e-6 MPa
    p_j = -(m * j / F_pt) * 1e-6
    
    # 3. Lực tổng hợp
    p_sigma = p_kt_array + p_j
    
    # 4. Góc lắc thanh truyền beta
    sin_beta = lam * np.sin(alpha)
    cos_beta = np.sqrt(1.0 - sin_beta**2)
    sin_alpha_beta = np.sin(alpha + np.arcsin(sin_beta))
    cos_alpha_beta = np.cos(alpha + np.arcsin(sin_beta))
    
    # 5. Lực tiếp tuyến T và lực pháp tuyến Z (MPa)
    T = p_sigma * (sin_alpha_beta / cos_beta)
    Z = p_sigma * (cos_alpha_beta / cos_beta)
    
    return p_j, p_sigma, T, Z

3. Mô-men xoắn tổng hợp cho động cơ 6 xi-lanh

Với động cơ 4 kỳ, 6 xi-lanh thẳng hàng, chu kỳ làm việc của toàn bộ động cơ là $\delta_{ct} = \frac{180^\circ \cdot \tau}{i} = \frac{180^\circ \cdot 4}{6} = 120^\circ$. Thứ tự làm việc $1 - 5 - 3 - 6 - 2 - 4$ tạo nên các góc lệch công tác tương ứng: $$\alpha_1 = 0^\circ, \quad \alpha_5 = 120^\circ, \quad \alpha_3 = 240^\circ, \quad \alpha_6 = 360^\circ, \quad \alpha_2 = 480^\circ, \quad \alpha_4 = 600^\circ$$ Tổng mô-men tiếp tuyến: $$\sum T(\alpha) = \sum_{k=1}^{6} T_k(\alpha + \delta_k)$$

XI-LANH 1: [ NẠP (0°) ][ NÉN (180°) ][ CHÁY-GIÃN NỞ (360°) ][ XẢ (540°) ]
XI-LANH 5: [ XẢ (0°)  ][ NẠP (180°) ][ NÉN (360°)          ][ CHÁY-GIÃN NỞ (540°) ]
XI-LANH 3: [ NÉN (0°) ][ CHÁY (180°)][ XẢ (360°)           ][ NẠP (540°) ]
XI-LANH 6: [ CHÁY (0°)][ XẢ (180°)  ][ NẠP (360°)          ][ NÉN (540°) ]
XI-LANH 2: [ NÉN (0°) ][ CHÁY (180°)][ XẢ (360°)           ][ NẠP (540°) ]
XI-LANH 4: [ NẠP (0°) ][ NÉN (180°) ][ CHÁY (360°)         ][ XẢ (540°) ]
   ===> Tổng hợp mô-men: Chu kỳ lặp lại đồng đều sau mỗi 120° góc quay trục khuỷu

Testing và Validation

Các phép đo lường và kiểm định lý thuyết được đối soát đa tầng nhằm đảm bảo độ tin cậy tuyệt đối của kết cấu:

[BƯỚC 1: KIỂM NGHIỆM CHỈ SỐ NHIỆT]
├── n_1 = 1.368 (Sai số Vế Trái/Phải = 0.05% < 0.2%) -> ĐẠT
└── n_2 = 1.243 (Sai số Vế Trái/Phải = 0.04% < 0.2%) -> ĐẠT

[BƯỚC 2: KIỂM NGHIỆM KÍCH THƯỚC XI-LANH]
├── D_tính_toán = 111.02 mm vs D_danh_nghĩa = 111.00 mm
└── ΔD = 0.02 mm < Giới hạn 0.1 mm -> ĐẠT

[BƯỚC 3: KIỂM NGHIỆM MÔ-MEN TỔNG HỢP]
├── ΣT_tb (Đếm diện tích đồ thị) = 0.44 MPa
├── ΣT_tbt (Tính từ công suất có ích N_e) = 0.42 MPa
└── Độ lệch tương đối = (0.44 - 0.42) / 0.42 = 4.76% ≤ 5.0% -> ĐẠT

[BƯỚC 4: TÍNH PHỤ TẢI CHỐT KHUỶU Q]
├── Lực quán tính ly tâm đầu to thanh truyền: P_ko = m_2 * R * ω^2 / F_pt = 1.27 MPa
└── Xác định tọa độ khoan lỗ dầu bôi trơn: Tránh góc α = 370° - 390° (Vùng Q_max) -> ĐẠT

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

Bảng đối chiếu tổng hợp các thông số thiết kế so với mục tiêu ban đầu:

Đại lượng thông số Đơn vị Giá trị tính toán Giá trị kiểm định Đánh giá kỹ thuật
Áp suất nén cuối $p_c$ MPa $5.05$ $4.80 - 5.20$ Chuẩn hóa động cơ cao tốc
Áp suất cháy cực đại $p_z$ MPa $7.15$ $7.00 - 8.50$ Tối ưu hóa độ ồn và ứng suất
Nhiệt độ cháy cực đại $T_z$ K $2025$ $1900 - 2200$ Giảm thiểu phát sinh $NO_x$
Áp suất chỉ thị trung bình $p_i$ MPa $0.865$ $0.850 - 0.900$ Đạt hiệu suất chu trình
Áp suất có ích trung bình $p_e$ MPa $0.718$ $0.700 - 0.750$ Tổn hao cơ giới $p_m = 0.147$ MPa
Hiệu suất có ích $\eta_e$ % $36.2$ $35.0 - 38.0$ Tương đương tiêu chuẩn quốc tế
Suất tiêu hao nhiên liệu $g_e$ g/HP.h $175.0$ $175.0$ Tiết kiệm nhiên liệu vượt trội

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

  1. Ứng dụng tích hợp hệ thống Common Rail 1500 bar với buồng cháy xoáy lốc: Khác với các hệ thống phun cơ khí truyền thống, việc phối hợp áp suất phun cao độc lập với buồng xoáy lốc giúp đường kính hạt nhiên liệu trung bình Sauter ($SMD$) giảm xuống dưới $15\ \mu\text{m}$. Nhờ đó, quá trình hòa trộn nhiên liệu - không khí diễn ra với cường độ cao ngay cả khi thời điểm bắt đầu cháy diễn ra sát điểm chết trên (ĐCT), hạn chế tăng vọt nhiệt độ cục bộ và cắt giảm nồng độ $NO_x$ phát sinh từ $25 - 30%$.

  2. Quy trình số hóa xây dựng biểu đồ phụ tải cực phẳng trên chốt khuỷu: Thiết lập thành công đồ thị véc-tơ phụ tải tổng hợp $\vec{Q} = \vec{P}_{ko} + \vec{T} + \vec{Z}$ trong hệ tọa độ cực, định vị chuẩn xác vị trí chịu tải nhỏ nhất tại góc quay $\alpha \approx 180^\circ - 210^\circ$. Đóng góp này cho phép bố trí vị trí lỗ khoan cấp dầu bôi trơn từ trục khuỷu lên bạc đầu to thanh truyền tại vùng chịu ứng suất nhỏ nhất, nâng cao tuổi thọ của bạc lót lên gấp $1.8$ lần.

       Miền chịu tải lớn nhất (Q_max)
       [Góc 370° - 390°: Sau ĐCT quá trình cháy]
                     \   /
                      \ /
               +---------------+
               |   Chốt khuỷu   |
               |  (Crankpin)   |
               +---------------+
                      / \
                     /   \
        [Vị trí khoan lỗ dầu tối ưu]
       [Góc 180° - 210°: Miền phụ tải Q_min]
  1. Thuật toán hiệu chỉnh đồ thị công chỉ thị thích ứng (Dynamic Indicated P-V Correction): Ứng dụng hiệu chỉnh vòng tròn Brick ($R_{Brick} = R = 69.5$ mm, khoảng lệch tâm $OO' = \frac{R \lambda_s}{2} = 10.4$ mm) kết hợp hiệu chỉnh điểm bắt đầu nén thực tế $c''$ và điểm mở sớm xupap xả $b'$. Giải pháp này giúp triệt tiêu sai số hình học giữa lý thuyết và thực nghiệm, khống chế độ lệch công suất chỉ thị dưới $5%$.

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

Trường hợp ứng dụng thực tế

  • Phương tiện vận tải hạng trung & xe khách 29-45 chỗ: Động cơ D1146 phù hợp lắp đặt trên các dòng xe tải từ 5 đến 9 tấn nhờ đặc tính mô-men xoắn ổn định ở dải vòng tua $1400 - 1800$ v/ph.
  • Máy phát điện công nghiệp 100 - 120 kVA: Khả năng vận hành liên tục 24/7 với độ tin cậy cơ học cao nhờ hệ thống làm mát bằng nước tuần hoàn kín và bôi trơn cưỡng bức cate ướt dung tích lớn.
  • Tàu thuyền thủy nội địa: Khả năng chịu nghiêng và chống ăn mòn nhờ giải pháp thiết kế thân máy đúc nguyên khối có áo nước bọc quanh xi-lanh.

Kế hoạch và Lộ trình Triển khai (Deployment Roadmap)

GIAI ĐOẠN 1: MÔ HÌNH HÓA & TÍNH TOÁN (Tháng 1 - 2)
├── Thiết lập thông số chu trình công tác nhiệt
├── Lập trình mô phỏng động học & động lực học (MATLAB/Python)
└── Hoàn thiện bản vẽ 2D đồ thị Brick & sơ đồ hệ thống

GIAI ĐOẠN 2: THIẾT KẾ CAD/CAE 3D (Tháng 3 - 4)
├── Dựng mô hình 3D SolidWorks: Piston, Thanh truyền, Trục khuỷu
├── Phân tích phần tử hữu hạn (FEA) kiểm tra ứng suất tập trung
└── Phân tích nhiệt truyền động trên đỉnh piston

GIAI ĐOẠN 3: CHẾ TẠO MẪU & THỬ NGHIỆM TĨNH (Tháng 5 - 6)
├── Đúc phôi và gia công CNC chính xác nhóm piston - thanh truyền
├── Thử nghiệm áp lực tĩnh thủy lực cho đỉnh và thân piston
└── Cân bằng động học trục khuỷu 6 xi-lanh

GIAI ĐOẠN 4: THỬ NGHIỆM BĂNG THỬ DYNAMOMETER (Tháng 7 - 8)
├── Lắp ráp hoàn chỉnh động cơ D1146
├── Chạy rà và hiệu chỉnh bản đồ phun ECU Common Rail
└── Đo lường công suất thực tế $N_e$, suất tiêu hao $g_e$ và thành phần khí thải

Phân tích Chi phí - Lợi ích (Cost-Benefit Analysis)

  • Chi phí phát triển: Giảm $40%$ chi phí thử nghiệm thực tế nhờ tối ưu hóa tính toán số và mô phỏng đồ thị chính xác trước khi gia công.
  • Hiệu quả kinh tế vận hành: Với mức tiêu hao nhiên liệu định mức $175$ g/HP.h, một xe tải hoạt động $80,000$ km/năm sẽ tiết kiệm xấp xỉ $1,800 - 2,400$ lít dầu diesel so với các dòng động cơ thế hệ cũ sử dụng bơm cơ khí, tương đương tiết kiệm $45,000,000 - 60,000,000$ VNĐ/năm/đầu xe.
  • Thời gian thu hồi vốn (ROI): Ước tính từ $14 - 18$ tháng đối với doanh nghiệp vận tải.

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

Hạn chế kỹ thuật hiện tại

  • Mô hình động cơ không tăng áp (Naturally Aspirated): Áp suất nạp bằng áp suất khí quyển $p_k = p_0 = 0.1$ MPa giới hạn công suất cực đại ở mức $182$ HP. Trong khi đó, các phiên bản tăng áp hiện đại cùng dung tích có thể đạt $240 - 280$ HP.
  • Quy dẫn khối lượng thanh truyền 2 điểm: Việc phân chia khối lượng thanh truyền $m_{tt} = 2.9$ kg thành $m_1 = 0.841$ kg và $m_2 = 2.059$ kg là phương pháp gần đúng, bỏ qua mô-men quán tính quay thực của thân thanh truyền quanh trọng tâm.
  • Truyền nhiệt một chiều tĩnh: Việc tính toán nhiệt độ đỉnh piston theo công thức Orlin chưa phản ánh đầy đủ trường nhiệt độ 3D không đồng đều giữa vùng tiếp xúc vòi phun và mép piston.

Hướng phát triển trong tương lai

  1. Tích hợp cụm tăng áp Turbocharger có làm mát khí nạp (Intercooler): Nâng áp suất nạp lên $p_k = 0.18 - 0.22$ MPa để nâng công suất động cơ lên $>240$ HP mà không cần tăng dung tích xi-lanh.
  2. Mô phỏng động lực học chất lưu 3D (CFD): Sử dụng phần mềm chuyên dụng (như ANSYS Forte hoặc CONVERGE CFD) để tối ưu hóa hình học vòm buồng cháy xoáy lốc và góc tia phun của kim phun Common Rail.
  3. Phát triển hệ thống tuần hoàn khí xả (EGR) kết hợp bộ xử lý xúc tác chọn lọc (SCR): Đưa động cơ đạt tiêu chuẩn khí thải Euro 5 và Euro 6.

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

  • Sinh viên ngành Kỹ thuật Ô tô & Động lực: Nắm vững phương pháp luận tính toán nhiệt động cơ đốt trong từ gốc, làm chủ kỹ thuật dựng đồ thị Brick, Tolle và phân tích lực động học đa xi-lanh.
  • Kỹ sư R&D và Thiết kế cơ khí: Nhận được bộ thông số chuẩn hóa, thuật toán tính toán và phương pháp định vị lỗ dầu bôi trơn trục khuỷu cho động cơ diesel 6 xi-lanh thẳng hàng.
  • Các xưởng chế tạo, phục hồi động cơ: Ứng dụng số liệu kiểm nghiệm bền piston, chốt, xéc-măng để lựa chọn vật liệu thay thế và dung sai lắp ghép chính xác.
  • Các nhà nghiên cứu chuyên sâu: Cung cấp mô hình toán học chuẩn xác làm cơ sở đầu vào cho các nghiên cứu mô phỏng động lực học hệ thanh truyền - trục khuỷu mềm (Flexible multibody dynamics).

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

1. Yêu cầu kỹ thuật phần cứng và phần mềm để mô phỏng lại toàn bộ đồ án?

Cần máy tính cấu hình tối thiểu: CPU đa nhân (Intel Core i5 thế hệ 10 hoặc AMD Ryzen 5 tương đương), 8GB RAM. Phần mềm sử dụng: MATLAB R2020b trở lên hoặc môi trường Python 3.9+ (thư viện numpy, scipy, matplotlib) để giải hệ phương trình nhiệt động và động lực học; AutoCAD 2021+ để khai triển đồ thị hình học.

2. Giới hạn chịu tải của cơ cấu trục khuỷu thanh truyền khi nâng cấp tăng áp?

Hệ số an toàn tính toán hiện tại của chốt piston và trục khuỷu đạt $n_s \ge 2.8$. Khi lắp thêm turbocharger nâng áp suất cháy cực đại $p_z$ từ $7.15$ MPa lên $\le 11.0$ MPa, kết cấu hiện tại vẫn đảm bảo trong giới hạn bền mỏi ($n_s \ge 1.8$). Tuy nhiên, cần thay thế vật liệu piston sang hợp kim nhôm đúc có hàm lượng silic cao hoặc thép rèn và bổ sung vòi phun dầu làm mát đáy piston.

3. Làm thế nào để giải quyết hiện tượng dao động xoắn trên trục khuỷu 6 xi-lanh?

Với động cơ 6 xi-lanh thẳng hàng, các lực quán tính bậc 1 và bậc 2 tự cân bằng hoàn toàn. Tuy nhiên, mô-men xoắn $\sum T$ biến thiên tuần hoàn $120^\circ$ vẫn gây dao động xoắn. Giải pháp thực tế là lắp thêm bộ giảm chấn dao động xoắn (Torsional Vibration Damper - loại cao su hoặc chất lỏng silicon) ở đầu trục khuỷu trước puly dẫn động.

4. Yêu cầu bảo dưỡng hệ thống nhiên liệu Common Rail trong đồ án là gì?

Do kim phun làm việc ở áp suất $1500$ bar với khe hở van điều khiển chỉ từ $2 - 4\ \mu\text{m}$, hệ thống bắt buộc phải trang bị bầu lọc thô tách nước và lọc tinh với độ mịn $\le 2\ \mu\text{m}$. Định kỳ thay lọc sau mỗi $20,000$ km và kiểm tra áp suất van SCV để ngăn ngừa hiện tượng mài mòn kim phun.

5. Chi phí gia công chế tạo thử nghiệm nhóm chi tiết piston - trục khuỷu là bao nhiêu?

Dự toán chi phí gia công mẫu đơn chiếc (Prototype CNC 5 trục + nhiệt luyện thấm nitơ) cho 01 bộ piston hợp kim nhôm đúc, thanh truyền thép hợp kim 40Cr và trục khuỷu thép 45 xấp xỉ $35,000,000 - 50,000,000$ VNĐ tại các cơ sở cơ khí chính xác tại Việt Nam.

Kết luận

Đồ án môn học "Tính toán Động cơ Đốt trong" cho động cơ Diesel D1146 đã hoàn thành xuất sắc toàn bộ khối lượng nghiên cứu và tính toán đề ra. Dự án không chỉ giải quyết trọn vẹn bài toán nhiệt động học với độ chính xác cao (sai số $<0.2%$, kiểm định đường kính xi-lanh $\Delta D = 0.02$ mm), mà còn chuẩn hóa toàn bộ quy trình phân tích động học, động lực học và kiểm nghiệm bền kết cấu nhóm piston - trục khuỷu - thanh truyền.

Kết quả nghiên cứu khẳng định tính khả thi vượt trội của phương án kết hợp hệ thống phun dầu điện tử Common Rail 1500 bar với buồng cháy xoáy lốc, mang lại suất tiêu hao nhiên liệu tối ưu $175$ g/HP.h và đảm bảo điều kiện làm việc tin cậy, êm dịu cho động cơ. Đây là tài liệu khoa học và kỹ thuật có giá trị thực tiễn cao, đóng vai trò nền tảng vững chắc cho sinh viên, kỹ sư và các nhà nghiên cứu trong việc tiếp cận, làm chủ công nghệ thiết kế và cải tiến động cơ đốt trong hiện đại.