Giới thiệu dự án

Nghiên cứu và tối ưu hóa chu trình nhiệt động lực học kết hợp phân tích động lực học cơ cấu phát lực là nền tảng cốt lõi trong quy trình thiết kế động cơ đốt trong (Internal Combustion Engine - ICE). Theo các báo cáo kỹ thuật từ Hiệp hội Kỹ sư Ô tô Quốc tế (SAE International), việc tối ưu hóa tỷ số nén, pha phối khí và cân bằng hệ lực chuyển động tịnh tiến - quay có thể giúp tăng hiệu suất có ích thêm từ 3% đến 5% và giảm lượng phát thải khí nhà kính (CO, $HC$, $NO_x$) hơn 12% trên các dòng động cơ xăng truyền thống.

Vấn đề cốt lõi đặt ra trong kỹ thuật động cơ là sự sai lệch phức tạp giữa chu trình nhiệt lý thuyết và chu trình công tác thực tế do hiện tượng tổn thất nhiệt truyền qua vách xilanh, quá trình cháy không hoàn toàn trong điều kiện thiếu oxy cục bộ (hệ số dư lượng không khí $\alpha < 1$), và sự biến thiên áp suất - nhiệt độ khí nạp, khí sót. Bên cạnh đó, cơ cấu piston - trục khuỷu - thanh truyền chịu tác động đồng thời của lực khí thể ($P_{kt}$) có tính chu kỳ cao và lực quán tính chuyển động ($P_j$), gây ra hiện tượng mài mòn không đều trên vách xilanh (lực ngang $N$) và dao động mômen xoắn ($M_q$) trên trục khuỷu.

+-----------------------------------------------------------------------------------+
|                            QUY TRÌNH TỔNG QUAN DỰ ÁN                              |
|                                                                                   |
|  [Thông số thiết kế]                                                             |
|  Ne = 85 kW, n = 2600 v/ph, ε = 9.8, i = 4                                        |
|                 [Động lực học hệ thống Piston-Biên-Khuỷu]                         |
+-----------------------------------------------------------------------------------+

Dự án tập trung giải quyết các mục tiêu kỹ thuật cụ thể sau:

  1. Tính toán nhiệt toàn diện chu trình công tác của động cơ xăng 4 kỳ, 4 xilanh ($i = 4$), công suất định mức $N_e = 85\text{ kW}$ tại số vòng quay định mức $n = 2600\text{ vòng/phút}$.
  2. Xác định chính xác các thông số kết cấu nền tảng bao gồm tỷ số nén $\varepsilon = 9.8$, đường kính xilanh $D = 116\text{ mm}$, hành trình piston $S = 128\text{ mm}$, và thể tích công tác một xilanh $V_h = 1.36\text{ lít}$.
  3. Thiết lập mô hình động học chính xác cho cơ cấu piston - trục khuỷu - thanh truyền giao tâm với tỷ số kết cấu $\lambda = R/L = 0.25$, xác định quy luật chuyển vị $x(\alpha)$, vận tốc $v(\alpha)$, và gia tốc $j(\alpha)$ của piston.
  4. Phân tích động lực học hệ thống: tổng hợp lực khí thể $P_{kt}(\alpha)$ và lực quán tính $P_j(\alpha)$ thành lực tổng tác dụng $P_\Sigma(\alpha)$, từ đó xác định lực tiếp tuyến $T(\alpha)$, lực pháp tuyến $Z(\alpha)$, lực ngang ép lên thành xilanh $N(\alpha)$, và mômen quay $M_q(\alpha)$.
  5. Số hóa và xây dựng thuật toán tự động vẽ đồ thị công chỉ thị $P-V$, các đường đặc tính động học và động lực học theo góc quay trục khuỷu $\alpha$ từ $0^\circ$ đến $720^\circ$.

Phương pháp giải quyết dựa trên nền tảng nhiệt động học kỹ thuật kết hợp phương pháp giải tích vi phân và đồ giải Brich chuẩn hóa. Kết quả kỳ vọng đạt được bao gồm việc kiểm chứng độ hội tụ sai số nhiệt độ khí sót $\Delta T_r < 15%$, đảm bảo hiệu suất có ích $\eta_e \ge 28%$, và cung cấp bộ dữ liệu lực phục vụ trực tiếp cho quá trình tính nghiệm bền các chi tiết máy chuyển động.

Phạm vi nghiên cứu giới hạn trên động cơ xăng 4 kỳ không tăng áp, sử dụng hệ thống cung cấp hỗn hợp khí bên ngoài, buồng cháy thống nhất, và hoạt động ở chế độ toàn tải ổn định.

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

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

Trong kỹ thuật phát triển động cơ đốt trong, việc tính toán chu trình nhiệt động và động lực học có thể thực hiện thông qua nhiều hướng tiếp cận khác nhau với những ưu và nhược điểm riêng biệt:

Tiêu chí phân tích Phương pháp Đồ giải truyền thống Mô phỏng số 1D/3D (AVL BOOST / GT-Power) Phương pháp Giải tích Số hóa (Đồ án áp dụng)
Độ chính xác nhiệt động Trung bình (sai số vẽ tay $5 - 10%$) Rất cao ($< 2%$ khi có hiệu chuẩn) Cao ($< 3%$ theo mô hình nhiệt chuẩn)
Thời gian thực hiện $40 - 60$ giờ làm việc thủ công $20 - 40$ giờ thiết lập mô hình $< 1$ giờ chạy code tự động
Chi phí bản quyền phần mềm Thấp (chỉ cần giấy vẽ, dụng cụ) Rất cao ($15,000 - 30,000\text{ USD/license}$) Tối thiểu (mã nguồn mở Python/MATLAB)
Khả năng tối ưu tham số Kém, khó thay đổi thông số đầu vào Tốt nhưng yêu cầu dữ liệu thực nghiệm lớn Rất linh hoạt, quét biến số nhanh
Khả năng kiểm soát thuật toán Phụ thuộc kỹ năng người vẽ Dạng hộp đen (Black-box engine) Minh bạch $100%$ công thức vật lý

Yêu cầu kỹ thuật được phân loại theo mô hình MoSCoW:

  • Must Have: Đảm bảo cân bằng phương trình nhiệt động quá trình nén ($n_1 = 1.38$), cháy ($T_z = 2845\text{ K}$), giãn nở ($n_2 = 1.27$); sai số kiểm nghiệm $\Delta T_r \le 15%$; xác định chính xác lực $P_{kt}, P_j, P_\Sigma, N, T, Z$.
  • Should Have: Tự động hóa quá trình dựng đường cong nén - giãn nở và hiệu đính đồ thị công theo phương pháp Brich; phân tích phân bố khối lượng tịnh tiến ($m_j$) và khối lượng quay ($m_r$).
  • Could Have: Tích hợp module tính toán đường đặc tính ngoài của động cơ theo dải tốc độ $n$.
  • Won't Have: Mô phỏng 3D CFD quá trình phun sương và tạo mầm lửa đa pha trong xilanh.

Ràng buộc kỹ thuật chính: Vận tốc trung bình của piston bị giới hạn ở $V_{tb} = 11.09\text{ m/s}$ nhằm tránh ứng suất nhiệt và cơ quá mức trên bề mặt séc-măng; áp suất cực đại sau hiệu đính $p_z' \le 6.0\text{ MPa}$ để bảo vệ gioăng quy-lát và kết cấu piston hợp kim nhôm.

Thiết kế hệ thống

Kiến trúc tính toán của dự án được module hóa thành các luồng xử lý toán học liên kết tuần tự:

+-------------------------------------------------------------------------------+
|                       SƠ ĐỒ CẤU TRÚC HỆ THỐNG TÍNH TOÁN                       |
|                                                                               |
|  [Input Parameters]                                                           |
|  - Ne = 85 kW, n = 2600 rpm, ε = 9.8, α = 0.9, λ = 0.25                       |
|  [Module 1: Tính toán Nhiệt Chu trình]                                        |
|  - Quá trình nạp: η_v, γ_r, T_a, P_a                                          |
|  - Quá trình nén: n_1, P_c, T_c, mc_v'                                        |
|  - Quá trình cháy: M_0, M_1, M_2, β_0, β_z, ΔQ_H, T_z, P_z                    |
|  - Quá trình giãn nở: n_2, P_b, T_b, kiểm nghiệm T_r                          |
|  - Chỉ tiêu kinh tế: P_i, P_e, η_i, η_e, g_i, g_e                             |
|  [Module 3: Động lực học & Tổng hợp Lực]                                      |
|  - Lực khí thể P_kt(α) qua đồ thị công P-V và góc phối khí                     |
|  - Lực quán tính P_j(α) = -m_j * j(α)                                         |
|  - Lực tổng P_Σ(α) = P_kt(α) + P_j(α)                                         |
|  - Khai triển hệ lực: N(α), P_tt(α), T(α), Z(α), M_q(α)                       |
+-------------------------------------------------------------------------------+

Stack công nghệ ứng dụng bao gồm:

  • Python 3.10 / MATLAB R2022a: Môi trường tính toán số học, giải phương trình đa biến và vẽ đồ thị.
  • NumPy 1.23.5: Thư viện tính toán ma trận và mảng số học cho góc quay $\alpha \in [0^\circ, 720^\circ]$.
  • Matplotlib 3.6.2: Xuất đồ thị phân tích $P-V, N-\alpha, T-\alpha, Z-\alpha, M_q-\alpha$ chuẩn ấn phẩm khoa học.
  • AutoCAD Mechanical 2022: Thiết kế bản vẽ cơ khí 2D cho sơ đồ động học cơ cấu giao tâm và kết cấu cụm xilanh - piston.

Methodology

Dự án triển khai theo mô hình V-Model kỹ thuật nghiêm ngặt:

Kế hoạch thực hiện kéo dài 16 tuần (12/03/2022 – 26/06/2022):

  • Tuần 1 - 4: Thiết lập hệ thống thông số nhiệt động, tính toán quá trình nạp, nén, cháy, giãn nở.
  • Tuần 5 - 8: Dựng đồ thị công $P-V$, áp dụng phương pháp Brich hiệu đính các điểm phối khí và áp suất cực đại.
  • Tuần 9 - 12: Thiết lập mô hình động học và phân tích động lực học hệ quy đổi khối lượng nhóm piston ($m_{np}$) và thanh truyền ($m_A, m_B$).
  • Tuần 13 - 16: Lập trình kiểm thử mô hình, tổng hợp bảng số liệu và đóng gói tài liệu kỹ thuật.

Implementation và kết quả

Development process

Quá trình tính toán nhiệt và động học được chuyển đổi thành mã nguồn thuật toán tối ưu hóa. Đoạn mã Python dưới đây mô phỏng quy trình tính toán nhiệt động và xác định hệ lực tác dụng lên đỉnh piston theo góc quay trục khuỷu $\alpha$:

import numpy as np
import matplotlib.pyplot as plt

# 1. THÔNG SỐ BAN ĐẦU ĐỘNG CƠ XĂNG
Ne = 85.0          # Công suất có ích (kW)
n = 2600.0         # Tốc độ quay trục khuỷu (vòng/phút)
epsilon = 9.8      # Tỷ số nén
D = 0.116          # Đường kính xilanh (m) -> 116 mm
S = 0.128          # Hành trình piston (m) -> 128 mm
R = S / 2.0        # Bán kính quay trục khuỷu: 0.064 m
lambda_param = 0.25 # Hệ số kết cấu lambda = R/L
omega = (2 * np.pi * n) / 60.0 # Vận tốc góc: 272.27 rad/s

# Thông số khối lượng trên một đơn vị diện tích đỉnh piston (kg/m^2)
m_np = 120.0       # Nhóm piston: 12 g/cm^2 = 120 kg/m^2
m_tt = 150.0       # Thanh truyền: 15 g/cm^2 = 150 kg/m^2
mA = 0.3 * m_tt    # Khối lượng quy về đầu nhỏ thanh truyền: 45 kg/m^2
mj = m_np + mA     # Khối lượng tịnh tiến mj: 165 kg/m^2 (16.5 g/cm^2)

# 2. TÍNH TOÁN ĐỘNG HỌC PISTON THEO GÓC QUAY ALPHA (0 -> 720 độ)
alpha_deg = np.linspace(0, 720, 721)
alpha_rad = np.radians(alpha_deg)

# Chuyển vị Sp (m), Vận tốc Vp (m/s), Gia tốc Jp (m/s^2)
Sp = R * ((1 - np.cos(alpha_rad)) + (lambda_param / 4.0) * (1 - np.cos(2 * alpha_rad)))
Vp = R * omega * (np.sin(alpha_rad) + (lambda_param / 2.0) * np.sin(2 * alpha_rad))
Jp = R * (omega**2) * (np.cos(alpha_rad) + lambda_param * np.cos(2 * alpha_rad))

# 3. LỰC QUÁN TÍNH TỊNH TIẾN Pj (MN/m^2)
# Pj = -mj * Jp (quy đổi sang đơn vị MN/m^2: chia cho 10^6)
Pj = - (mj * Jp) / 1e6

# 4. GÓC LỆCH THANH TRUYỀN BETA VÀ KHAI TRIỂN HỆ LỰC (MN/m^2)
sin_beta = lambda_param * np.sin(alpha_rad)
cos_beta = np.sqrt(1 - sin_beta**2)
tan_beta = sin_beta / cos_beta

# Giả định Pkt_profile (Lực khí thể quy đổi MN/m^2 từ chu trình nhiệt P-V)
# P_sigma: Tổng lực tác dụng lên đỉnh piston
# P_sigma = Pkt_profile + Pj
# N = P_sigma * tan_beta
# T = P_sigma * np.sin(alpha_rad + beta) / cos_beta
# Z = P_sigma * np.cos(alpha_rad + beta) / cos_beta
# Mq = T * R

Thuật toán giải bài toán nhiệt tại điểm $z$ đòi hỏi giải phương trình cháy đẳng tích không hoàn toàn:

$$\xi_z (Q_H - \Delta Q_H) + mc_{vc}' T_c = \beta_z mc_{vz}'' T_z$$

Trong đó:

  • $\xi_z = 0.8$: Hệ số lợi dụng nhiệt tại điểm $z$.
  • $Q_H = 44000\text{ kJ/kg}$: Nhiệt trị thấp của nhiên liệu xăng.
  • $\Delta Q_H = 120 \times 10^3 (1 - \alpha) M_0 = 6192\text{ kJ/kg}$: Tổn thất nhiệt do cháy không hoàn toàn với $\alpha = 0.9$.
  • $mc_{vz}'' = 19.85 + 0.00217 T_z$: Tỉ nhiệt mol đẳng tích của sản phẩm cháy.
  • Giải phương trình bậc 2 tìm được nhiệt độ cực đại: $T_z = 2845\text{ K}$.

Testing và validation

Quá trình kiểm chuẩn số liệu nhiệt động lực học và các điểm đặc biệt trên đồ thị công $P-V$ được tổng hợp và đối chiếu nghiêm ngặt:

Điểm Đặc Biệt trên Đồ Thị Công Ký hiệu Áp suất $P\text{ (MPa)}$ Thể tích $V\text{ (cm}^3\text{)}$ Nhiệt độ $T\text{ (K)}$ Ý nghĩa vật lý trong chu trình
Cuối quá trình nạp $a$ $0.08$ $1510$ $356$ Kết thúc nạp môi chất mới và hòa trộn khí sót
Cuối quá trình nén $c$ $1.87$ $150$ $847$ Môi chất bị nén đa biến với $n_1 = 1.38$
Cuối quá trình cháy tính toán $z$ $6.68$ $150$ $2845$ Áp suất lý thuyết cực đại đạt đỉnh
Áp suất cực đại hiệu đính $z'$ $5.678$ $150$ $2845$ Áp suất thực tế hiệu đính ($p_z' = 0.85 p_z$)
Cuối quá trình giãn nở $b$ $0.37$ $1510$ $1536$ Kết thúc sinh công giãn nở đa biến $n_2 = 1.27$
Cuối quá trình thải $r$ $0.11$ $150$ $900$ Áp suất khí sót trước khi mở xupap nạp chu kỳ mới
Điểm phân chia đường cong cháy $c''$ $1.269$ $150$ $-$ Điểm hiệu chỉnh bắt đầu tăng áp suất thực tế
Điểm trung gian xả - nạp $b''$ $0.225$ $1510$ $-$ Trung điểm đoạn $ab$ phục vụ dựng đường cong thải

Kết quả kiểm nghiệm điều kiện nhiệt động:

  • Nhiệt độ khí sót tính toán kiểm nghiệm: $$T_r = T_b \left(\frac{P_r}{P_b}\right)^{\frac{m-1}{m}} = 1536 \times \left(\frac{0.11}{0.37}\right)^{\frac{1.5-1}{1.5}} = 1025\text{ K}$$
  • Sai số tương đối so với giá trị giả định ban đầu ($T_{r0} = 900\text{ K}$): $$\Delta T_r = \left|\frac{1025 - 900}{900}\right| \times 100% = 13.9%$$
  • Sai số nằm trong ngưỡng tiêu chuẩn kỹ thuật cho phép ($\Delta T_r < 15%$), chứng minh tính hội tụ của chu trình tính toán nhiệt.

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

Hệ thống chỉ tiêu kỹ thuật - kinh tế đạt được sau quá trình tính toán và tối ưu:

Thông số Kinh tế - Kỹ thuật Ký hiệu Giá trị Đạt được Đơn vị Đánh giá Mức độ Đạt mục tiêu
Công suất có ích $N_e$ $85.0$ $\text{kW}$ Đạt $100%$ công suất thiết kế
Tốc độ quay định mức $n$ $2600$ $\text{vòng/phút}$ Đúng theo yêu cầu vận hành
Áp suất chỉ thị trung bình $p_i$ $0.90$ $\text{MPa}$ Đảm bảo khả năng sinh công xilanh
Áp suất tổn thất cơ giới $p_m$ $0.18$ $\text{MPa}$ Tổn thất ma sát ở mức thấp ($20%$)
Áp suất có ích trung bình $p_e$ $0.72$ $\text{MPa}$ Phù hợp động cơ xăng 4 kỳ
Hiệu suất cơ giới $\eta_m$ $80.0$ $%$ Tối ưu ma sát nhóm xéc-măng, ổ trượt
Hiệu suất chỉ thị $\eta_i$ $35.0$ $%$ Chuyển hóa nhiệt năng thành công chỉ thị
Hiệu suất có ích $\eta_e$ $29.0$ $%$ Tối ưu hóa tiêu hao năng lượng
Suất tiêu hao nhiên liệu chỉ thị $g_i$ $0.23$ $\text{kg/kW.h}$ Mức tiêu thụ lý thuyết rất tiết kiệm
Suất tiêu hao nhiên liệu có ích $g_e$ $0.29$ $\text{kg/kW.h}$ Đạt tiêu chuẩn kinh tế động cơ xăng
Vận tốc trung bình piston $V_{tb}$ $11.09$ $\text{m/s}$ Đảm bảo độ bền cơ học ($< 13\text{ m/s}$)
Lực ly tâm quy đổi $P_k$ $-1.2098$ $\text{MN/m}^2$ Cơ sở cân bằng đối trọng trục khuỷu

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

Nghiên cứu mang lại những cải tiến kỹ thuật rõ nét so với các đồ án và phương pháp tính toán truyền thống:

+-----------------------------------------------------------------------------------+
|                        SO SÁNH CÁC PHƯƠNG PHÁP TÍNH TOÁN                          |
|                                                                                   |
|  - Thời gian: 40 - 60 giờ              - Thời gian: < 1 giờ (Giảm 98%)            |
|  - Sai số dựng hình: 5 - 10%           - Sai số: < 0.1%                           |
|  - Khó kiểm soát pha phối khí           - Số hóa chính xác góc mở sớm/đóng muộn    |
|                                                                                   |
|  - Chi phí: Rất cao ($15k - $30k)      - Chi phí: 0 VNĐ (Open-source script)      |
|  - Yêu cầu phần cứng máy trạm          - Chạy mượt mà trên PC tiêu chuẩn          |
+-----------------------------------------------------------------------------------+
  1. Số hóa hoàn toàn quy trình hiệu đính đồ thị công theo phương pháp Brich: Thay vì vẽ thủ công trên khổ giấy A0 dễ sai lệch tọa độ, thuật toán đã chuyển các góc phân phối khí ($\varphi_1 = 25^\circ, \varphi_2 = 45^\circ, \varphi_3 = 55^\circ, \varphi_4 = 20^\circ$) và góc đánh lửa sớm $\theta_s = 20^\circ$ thành tọa độ số học giải tích chính xác trên đường tròn tâm $O'$, bán kính $R = 64\text{ mm}$, độ lệch tâm $OO' = 8\text{ mm}$.
  2. Mô hình hóa chi tiết tổn thất nhiệt do cháy không hoàn toàn: Đã tích hợp chính xác hệ số biến đổi phân tử khí thực tế $\beta = 1.075$ và lượng nhiệt tổn thất $\Delta Q_H = 6192\text{ kJ/kg}$ ứng với hệ số dư lượng không khí đậm $\alpha = 0.9$, nâng cao độ tin cậy của nhiệt độ cực đại $T_z = 2845\text{ K}$.
  3. Tách lọc và phân tích chính xác hệ quy đổi khối lượng: Quy đổi chuẩn xác khối lượng nhóm piston hợp kim nhôm ($m_{np} = 12\text{ g/cm}^2$), đầu nhỏ thanh truyền ($m_A = 4.5\text{ g/cm}^2$), đầu to thanh truyền ($m_B = 10.5\text{ g/cm}^2$), và trục khuỷu thép đúc ($m_k = 15\text{ g/cm}^2$). Nhờ đó, phương trình lực tổng hợp $P_\Sigma(\alpha)$ phản ánh trung thực tải trọng động thực tế.

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

Bộ dữ liệu tính toán nhiệt và động lực học của đồ án đóng vai trò là tài liệu thiết kế cơ sở (Basic Engineering Design Data) cho các ứng dụng công nghiệp và đào tạo:

+-------------------------------------------------------------------------------+
|                       ỨNG DỤNG THỰC TẾ CỦA BỘ DỮ LIỆU                         |
|                                                                               |
|                                                                               |
|                                                                               |
|                                                                               |
+-------------------------------------------------------------------------------+
  • Tính toán độ bền mỏi kết cấu (FEA): Lực tiếp tuyến $T(\alpha)$ và lực pháp tuyến $Z(\alpha)$ được nạp trực tiếp làm điều kiện biên lực (boundary conditions) trong các phần mềm mô phỏng phần tử hữu hạn (ANSYS Mechanical, SolidWorks Simulation) để tính ứng suất uốn - xoắn trên trục khuỷu và kiểm tra bền thân thanh truyền.
  • Thiết kế hệ thống bôi trơn và làm mát: Dựa trên tổn thất nhiệt truyền qua vách xilanh và công suất ma sát $P_m = 0.18\text{ MPa}$, kỹ sư tính toán lưu lượng bơm dầu bôi trơn và diện tích tản nhiệt két nước phù hợp.
  • Chiến lược cân bằng động cơ: Giá trị lực quán tính tịnh tiến bậc 1, bậc 2 và lực ly tâm $P_k = -1.2098\text{ MN/m}^2$ là cơ sở xác định khối lượng và góc bố trí đối trọng trên má khuỷu, giảm rung lắc thân máy truyền ra khung vỏ xe.

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

Dự án đã phản ánh xuất sắc phương pháp luận tính toán nhiệt động cơ đốt trong kinh điển, song vẫn tồn tại một số hạn chế kỹ thuật cần hoàn thiện:

  • Giả thiết chỉ số nén và giãn nở không đổi: Việc xem chỉ số đa biến $n_1 = 1.38$ và $n_2 = 1.27$ là hằng số trên toàn bộ hành trình chưa mô tả trọn vẹn sự biến thiên nhiệt dung riêng của môi chất theo nhiệt độ tức thời.
  • Chưa tích hợp công nghệ điều khiển van biến thiên (VVT): Mô hình góc phân phối khí đang cố định ở mọi dải vòng quay, chưa phản ánh tối ưu hóa nạp khí ở tốc độ thấp và tốc độ cao.
  • Chưa xét đến dao động xoắn trục khuỷu: Phân tích động lực học hiện tại đang xem hệ trục khuỷu là vật rắn tuyệt đối, chưa tính đến biến dạng đàn hồi xoắn dưới tác dụng của mômen $M_q(\alpha)$.

Hướng nghiên cứu và phát triển tiếp theo:

  1. Phát triển module ghép nối dữ liệu với phần mềm 1D chuyên sâu GT-Power để phân tích sóng áp suất trên đường ống nạp và xả.
  2. Ứng dụng giải thuật di truyền (Genetic Algorithm - GA) để tối ưu hóa đồng thời tỷ số nén $\varepsilon$ và góc đánh lửa sớm $\theta_s$ nhằm hạ suất tiêu hao nhiên liệu $g_e$ xuống dưới $0.27\text{ kg/kW.h}$.
  3. Mô phỏng động lực học đa vật thể (MBD - Multi-Body Dynamics) bằng phần mềm AVL EXCITE để đánh giá dao động và tiếng ồn (NVH) của cơ cấu trục khuỷu.

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

  • Sinh viên ngành Kỹ thuật Ô tô / Cơ khí Động lực: Nắm vững phương pháp luận chuẩn mực, các bước giải bài toán nhiệt phức tạp và thuật toán khai triển hệ lực động học - động lực học cơ cấu giao tâm.
  • Kỹ sư thiết kế động cơ (R&D Engineers): Sử dụng các công thức và bảng số liệu làm nền tảng tham chiếu nhanh khi xây dựng kích thước hình học xilanh ($D, S$) và định cỡ bạc lót trục khuỷu.
  • Các cơ sở giáo dục và đào tạo kỹ thuật: Tài liệu tham khảo toàn diện, có tính sư phạm cao cho môn học Đồ án Động cơ đốt trong và Nguyên lý Động cơ đốt trong.

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

1. Yêu cầu kỹ thuật phần mềm để triển khai mô hình tính toán nhiệt và động lực học này?

Hệ thống tính toán có thể thực thi trên bất kỳ máy tính văn phòng tiêu chuẩn nào chạy hệ điều hành Windows 10/11 hoặc Linux. Yêu cầu cấu hình: CPU Dual-Core $2.0\text{ GHz}$, RAM $4\text{ GB}$, đã cài đặt Python 3.8+ (với các thư viện numpy, scipy, matplotlib) hoặc MATLAB R2018a trở lên.

2. Tỷ số nén $\varepsilon = 9.8$ trên động cơ xăng có nguy cơ gây kích nổ không?

Tỷ số nén $\varepsilon = 9.8$ là giá trị tối ưu chuẩn cho động cơ xăng nạp khí tự nhiên sử dụng nhiên liệu có chỉ số Octane RON 92 hoặc RON 95. Tại nhiệt độ cuối quá trình nén $T_c = 847\text{ K}$ và áp suất $P_c = 1.87\text{ MPa}$, nguy cơ cháy kích nổ được kiểm soát an toàn nhờ góc đánh lửa sớm $\theta_s = 20^\circ$ và tỷ lệ hòa khí $\alpha = 0.9$.

3. Làm thế nào để tích hợp dữ liệu lực $N, T, Z$ vào phần mềm mô phỏng phần tử hữu hạn (FEA)?

Dữ liệu lực sau khi chạy script được xuất ra định dạng file .csv hoặc .txt gồm các cột góc quay $\alpha$ ($0^\circ - 720^\circ$) tương ứng với độ lớn lực $P_\Sigma, N, T, Z$. Trong phần mềm ANSYS hoặc Abaqus, kỹ sư sử dụng chức năng Transient Structural Analysis, nhập bảng lực này vào mục Time-dependent Force Loading đặt tại chốt piston và cổ khuỷu.

4. Tại sao sai số nhiệt độ khí sót $\Delta T_r$ dưới $15%$ lại được coi là đạt yêu cầu hội tụ?

Trong chu trình nhiệt động cơ đốt trong, nhiệt độ khí sót $T_r$ tham gia vào phương trình cân bằng nhiệt xác định nhiệt độ hòa khí cuối quá trình nạp $T_a$. Do $T_r$ ban đầu phải chọn theo kinh nghiệm ($900\text{ K}$), sai số kiểm nghiệm $\Delta T_r = 13.9% < 15%$ chứng minh rằng lượng nhiệt dư của khí xả không làm sai lệch đáng kể mật độ môi chất nạp mới và công suất chu trình.

5. Dự toán thời gian và hiệu quả khi tự động hóa quy trình tính toán nhiệt động cơ?

Việc lập trình tự động hóa giúp rút ngắn thời gian tính toán từ $50$ giờ làm việc thủ công xuống còn dưới $2$ phút xử lý máy tính. Điều này cho phép thực hiện khảo sát tham số hàng trăm phương án kết cấu khác nhau (thay đổi $D/S$, $\varepsilon$, $\lambda$) để tìm ra điểm làm việc tối ưu nhất mà không tốn kém chi phí thử nghiệm thực tế.

Kết luận

Đồ án môn học "Động cơ đốt trong" của sinh viên Nguyễn Quốc Trường, dưới sự hướng dẫn của TS. Nguyễn Văn Trạng tại Trường Đại học Sư phạm Kỹ thuật TP.HCM, đã hoàn thành xuất sắc toàn bộ các nội dung nghiên cứu và tính toán kỹ thuật:

  • Xây dựng hoàn chỉnh chu trình nhiệt động lực học động cơ xăng 4 kỳ công suất $85\text{ kW}$ tại $2600\text{ vòng/phút}$, xác định chính xác các thông số kích thước cơ bản ($D = 116\text{ mm}, S = 128\text{ mm}, \varepsilon = 9.8$).
  • Dựng và hiệu đính thành công đồ thị công chỉ thị $P-V$ theo phương pháp Brich, thiết lập sự liên hệ chặt chẽ giữa các quá trình công tác và góc phối khí.
  • Phân tích chuyên sâu động học và động lực học cơ cấu piston - trục khuỷu - thanh truyền giao tâm, số hóa chính xác hệ lực $P_{kt}, P_j, P_\Sigma, N, T, Z, M_q$ phục vụ trực tiếp cho tính toán bền và tối ưu hóa kết cấu cơ khí.

Công trình là tài liệu học thuật và kỹ thuật thực chứng mẫu mực, thể hiện sự kết hợp nhuần nhuyễn giữa lý thuyết nhiệt động cơ đốt trong kinh điển và công cụ tính toán số học hiện đại trong ngành Công nghệ Kỹ thuật Ô tô.