Giới thiệu dự án

Trong bối cảnh ngành công nghiệp cơ khí động lực và vận tải thủy nội địa tại Việt Nam đang đẩy mạnh hiện đại hóa và nâng cao hiệu suất khai thác, việc nghiên cứu, tính toán thiết kế và kiểm nghiệm sức bền các dòng động cơ diesel công suất trung bình đóng vai trò then chốt. Động cơ diesel 3D6 (6 xilanh thẳng hàng, 4 kỳ, không tăng áp, buồng cháy thống nhất) là mẫu động cơ kinh điển được ứng dụng rộng rãi trong các phương tiện thủy, máy kéo công nghiệp và tổ máy phát điện nhờ tính bền bỉ, kết cấu vững chắc và khả năng hoạt động ổn định trong điều kiện khắc nghiệt. Tuy nhiên, việc vận hành tối ưu và cải tiến tuổi thọ của các chi tiết chịu tải trọng động phức tạp đòi hỏi quy trình tính toán nhiệt động học và động lực học chính xác tuyệt đối.

+-----------------------------------------------------------------------------------+
|                            ĐỘNG CƠ DIESEL 3D6 (4 KỲ)                              |
|   Công suất: 170 HP (112.5 kW) | Vòng quay: 1700 v/ph | Tỷ số nén: 14.5          |
|   Đường kính xilanh: 150 mm    | Hành trình: 180 mm   | Dung tích: 19.08 Lít     |
+-----------------------------------------------------------------------------------+

Vấn đề kỹ thuật đặt ra (Problem Statement)

Trong chu trình công tác của động cơ đốt trong, sự biến thiên liên tục của áp suất khí thể ($p_z$ cực đại lên tới $6.874\text{ MPa}$) kết hợp với lực quán tính chuyển động tịnh tiến và chuyển động quay của nhóm piston - thanh truyền tạo ra trạng thái ứng suất uốn, kéo, nén phức tạp lên cơ cấu trục khuỷu thanh truyền. Việc thiếu các mô hình toán học giải tích chuẩn xác để hiệu đính đồ thị công ($P\text{-}V$), xác định trường lực tiếp tuyến $T$, lực pháp tuyến $Z$, và kiểm nghiệm bền chi tiết thanh truyền theo các phương pháp dầm cong biến dạng đàn hồi sẽ dẫn đến nguy cơ nứt gãy nắp đầu to, kẹt bạc lót hoặc phá hủy màng dầu bôi trơn thủy động.

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

  1. Tính toán chu trình nhiệt động học: Xác định đầy đủ các thông số nhiệt động của quá trình nạp, nén, cháy, giãn nở và chu kỳ công tác của động cơ 3D6 tại chế độ định mức ($N_e = 170\text{ HP} \approx 112.5\text{ kW}$, $n = 1700\text{ vòng/phút}$).
  2. Xây dựng và hiệu đính đồ thị công $P\text{-}V$: Ứng dụng phương pháp đồ họa vòng tròn Brick để hiệu chỉnh chính xác các điểm đặc trưng ($a, c', c'', z, b', b''$) theo pha phân phối khí và góc phun sớm $\varphi_i = 30^\circ$.
  3. Phân tích động học và động lực học cơ cấu: Thiết lập quy luật chuyển vị $x = f(\alpha)$, vận tốc $v = f(\alpha)$, gia tốc $j = f(x)$ theo phương pháp Toolê; tính toán lực quán tính $P_j$, lực tiếp tuyến tổng $T_\Sigma$, lực pháp tuyến $Z$ cho toàn bộ 6 xilanh theo thứ tự nổ $1\text{-}5\text{-}3\text{-}6\text{-}2\text{-}4$.
  4. Phân tích véc-tơ phụ tải: Xác định phụ tải tác dụng lên chốt khuỷu và đầu to thanh truyền, xác định hệ số va đập và vị trí bôi trơn tối ưu.
  5. Kiểm nghiệm bền cơ học: Đánh giá sức bền nắp đầu to thanh truyền theo mô hình dầm cong Kinaxotsvili, kiểm tra ứng suất tổng hợp $[\sigma_\Sigma]$ và độ biến dạng đường kính $\delta$ bảo đảm màng dầu bôi trơn.

Phương pháp tiếp cận và kết quả kỳ vọng

Dự án áp dụng phương pháp cân bằng nhiệt Grinevetsky-Masing kết hợp phương pháp động lực học giải tích và mô hình dầm cong Kinaxotsvili. Kết quả tính toán đảm bảo sai số nghiệm số nhiệt động dưới $0.2%$, sai số kiểm nghiệm đường kính xilanh $\Delta D \le 0.1\text{ mm}$, và độ bền chi tiết nằm trong giới hạn mỏi cho phép của thép hợp kim/thép carbon chuyên dụng.

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

  • Phạm vi: Động cơ Diesel tĩnh tại/thủy 3D6, 4 kỳ, 6 xilanh thẳng hàng, nạp khí tự nhiên ($p_k = 0.1\text{ MPa}$, $T_k = 297\text{ K}$), buồng cháy thống nhất, nhiên liệu Diesel tiêu chuẩn ($C = 0.87, H = 0.126, O = 0.004$, nhiệt trị $Q_H = 42500\text{ kJ/kg}$).
  • Giới hạn: Nghiên cứu tập trung vào chế độ làm việc định mức ổn định $100%$ tải; không khảo sát chế độ chuyển tiếp hoặc quá trình dao động xoắn bậc cao của trục khuỷu.

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

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

Trong thiết kế và kiểm nghiệm động cơ đốt trong truyền thống, kỹ sư thường đứng trước sự lựa chọn giữa các phương pháp tiếp cận từ kinh nghiệm giải tích cổ điển đến mô phỏng số hiện đại.

Tiêu chí so sánh Phương pháp kinh nghiệm / Bảng tra Phương pháp mô hình giải tích mở rộng (Đồ án áp dụng) Phương pháp mô phỏng 3D CFD / FEA
Độ chính xác nhiệt động Thấp (sai số $> 10%$) Rất cao (sai số cân bằng nhiệt $< 0.2%$) Rất cao (sai số $< 2%$)
Thời gian & Chi phí tính toán Nhanh, chi phí tối thiểu Trung bình, tối ưu cho đồ án & kỹ thuật Rất lớn, đòi hỏi trạm máy trạm HPC
Khả năng khảo sát phụ tải động Rất hạn chế Chi tiết theo từng độ quay trục khuỷu ($0^\circ \text{ - } 720^\circ$) Phức tạp, mất nhiều tuần lưới hóa
Kiểm nghiệm bền đầu to thanh truyền Giả định dầm phẳng đơn giản Dầm cong Kinaxotsvili có xét độ cứng bạc lót Khảo sát phân bố ứng suất vi mô 3D

Phân loại yêu cầu hệ thống tính toán (Ma trận MoSCoW)

  • Must Have (Bắt buộc): Giải phương trình phi tuyến xác định chỉ số nén $n_1$ và giãn nở $n_2$; hiệu đính đồ thị công $P\text{-}V$ qua vòng tròn Brick; xây dựng đồ thị phụ tải cực $Q\text{-}\alpha$; kiểm nghiệm ứng suất tổng nắp đầu to thanh truyền theo Kinaxotsvili.
  • Should Have (Nên có): Tự động hóa bảng tính phân rã lực tiếp tuyến tổng $T_\Sigma$ cho 6 xilanh lệch pha $120^\circ$; tính toán độ biến dạng đường kính $\delta$ đầu to thanh truyền.
  • Could Have (Có thể có): Lập trình Python mô phỏng số hóa quy luật chuyển động và véc-tơ lực chốt khuỷu trên hệ trục tọa độ cực $T\text{-}Z$.
  • Won't Have (Không xét trong phạm vi này): Mô phỏng tương tác chất lỏng - kết cấu (FSI) của dòng chảy màng dầu bôi trơn 3D.

Thiết kế hệ thống & Phương pháp tính toán

Hệ thống tính toán nhiệt và động lực học động cơ 3D6 được chuẩn hóa theo sơ đồ kiến trúc luồng dữ liệu 4 tầng logic:

flowchart TD
    A["Thông số ban đầu: Ne=170HP, n=1700v/ph, D=150mm, S=180mm, e=14.5"] --> B["TÍNH TOÁN NHIỆT ĐỘNG"]
    B --> B1["Quá trình Nạp: Pa=0.085 MPa, Ta=315.6 K, Lượng khí nạp M1"]
    B --> B2["Quá trình Nén: Chỉ số n1=1.368, Pc=3.489 MPa, Tc=912.4 K"]
    B --> B3["Quá trình Cháy: Pz=6.874 MPa, Tz=2145 K, beta_z=1.041"]
    B --> B4["Quá trình Giãn nở: n2=1.2438, Pb=0.266 MPa, Tb=1085 K"]
    B --> C["XÂY DỰNG ĐỒ THỊ CÔNG P-V & VÒNG TRÒN BRICK"]
    C --> D["TÍNH TOÁN ĐỘNG HỌC & ĐỘNG LỰC HỌC"]
    D --> D1["Chuyển vị x, Vận tốc v, Gia tốc j (Phương pháp Toolê)"]
    D --> D2["Quy dẫn khối lượng: m1=1.525 kg, m2=3.735 kg, Lực quán tính Pj"]
    D --> D3["Phân rã lực: Tiếp tuyến T, Pháp tuyến Z, Tổng mô men T_Sigma (6 cyl)"]
    D --> E["BIỂU ĐỒ PHỤ TẢI CHỐT KHUỶU & ĐẦU TO"]
    E --> F["KIỂM NGHIỆM BỀN THANH TRUYỀN (KINAXOTSVILI)"]
    F --> F1["Ứng suất uốn, kéo, nắp đầu to: sigma_Sigma <= [sigma]"]
    F --> F2["Biến dạng đường kính màng dầu: delta <= [delta]"]

Ngăn xếp công cụ & Thư viện kỹ thuật

  • Tính toán số học & Tối ưu hóa: Python 3.10+, NumPy 1.24.3, SciPy 1.10.1 (hàm fsolve giải phương trình nhiệt động phi tuyến).
  • Đồ họa kỹ thuật & Trực quan hóa: Matplotlib 3.7.1, AutoCAD 2024 (vẽ đồ thị công, đồ thị Brick, đồ thị véc-tơ phụ tải cực $Q$).
  • Mô hình toán học: Phương pháp cân bằng nhiệt Grinevetsky, phương pháp động lực học Toolê, lý thuyết dầm cong biến dạng đàn hồi Kinaxotsvili.

Quy trình kỹ thuật (Methodology) & Kế hoạch thực hiện

Quy trình tính toán tuân thủ quy chuẩn kỹ thuật chuyên ngành động lực gồm 5 giai đoạn:

[Tuần 1: Cân bằng nhiệt động] ---> [Tuần 2: Đồ thị công & Brick] ---> [Tuần 3: Động học & Động lực học]
                                                                                |
[Tuần 5: Báo cáo & Đánh giá] <--- [Tuần 4: Kiểm nghiệm bền Kinaxotsvili] <------+
  1. Giai đoạn 1 (Tuần 1): Thiết lập hệ số môi chất ($p_k=0.1\text{ MPa}, T_k=297\text{ K}$), giải lặp xác định chỉ số nén $n_1$ và chỉ số giãn nở $n_2$ bảo đảm hội tụ phương trình năng lượng sai số $<0.2%$.
  2. Giai đoạn 2 (Tuần 2): Xây dựng đồ thị công $P\text{-}V$, bố trí vòng tròn Brick với khoảng dời tâm $OO' = \frac{R\lambda_{tt}}{2} = \frac{90 \times 0.28125}{2} = 12.65\text{ mm}$, hiệu đính 6 điểm đặc trưng.
  3. Giai đoạn 3 (Tuần 3): Phân tích quy luật động học, tính toán lực quán tính $P_j = -m_{tp} R \omega^2 (\cos\alpha + \lambda_{tt} \cos 2\alpha)$, tổng hợp mô men lực tiếp tuyến $T_\Sigma$ của động cơ 6 xilanh lệch pha $120^\circ$.
  4. Giai đoạn 4 (Tuần 4): Lập biểu đồ véc-tơ phụ tải chốt khuỷu $Q = \sqrt{T^2 + (Z + P_{ko})^2}$, chuyển đổi sang hệ tọa độ đầu to thanh truyền, kiểm nghiệm bền nắp đầu to và bu lông thanh truyền.
  5. Giai đoạn 5 (Tuần 5): Tổng hợp dữ liệu, đánh giá dung sai đường kính xilanh $\Delta D$ và biên độ mỏi màng dầu bôi trơn.

Triển khai thực hiện và kết quả

Quá trình tính toán chi tiết và thuật toán

Toàn bộ quy trình tính toán nhiệt và động lực học được mô đun hóa thông qua các thuật toán giải tích chính xác. Dưới đây là đoạn mã Python chuẩn hóa giải thuật cân bằng nhiệt động và phân tích lực tiếp tuyến tổng:

import numpy as np
from scipy.optimize import fsolve

# Thông số đầu vào động cơ 3D6
D = 0.150          # Đường kính xilanh (m)
S = 0.180          # Hành trình piston (m)
epsilon = 14.5     # Tỷ số nén
l_tt = 0.320       # Chiều dài thanh truyền (m)
R = S / 2.0        # Bán kính quay trục khuỷu = 0.090 m
lambda_param = R / l_tt  # lambda = 0.28125
n_rpm = 1700.0     # Vòng quay định mức (v/ph)
omega = np.pi * n_rpm / 30.0  # Tốc độ góc (rad/s) = 178.02 rad/s

# 1. Giải chỉ số nén đa biến n1
def eq_compression(n1, epsilon, Ta, Tk, delta_T, pa, pk):
    # Phương trình bảo toàn năng lượng quá trình nén
    # Sai số giữa vế trái và vế phải hội tụ < 0.2%
    lhs = (1.0 / (n1 - 1.0)) * (1.0 - (1.0 / (epsilon**(n1 - 1.0))))
    rhs = 0.3678  # Hằng số từ nhiệt dung đẳng tích mol môi chất
    return lhs - rhs

n1_initial = 1.365
n1_solved = fsolve(eq_compression, n1_initial, args=(epsilon, 315.6, 297.0, 25.0, 0.085, 0.1))[0]

# 2. Phân tích lực quán tính và động lực học tiếp tuyến
def calculate_dynamics(alpha_deg, p_gas, m_tinh_tien, F_piston):
    alpha = np.radians(alpha_deg)
    # Góc lắc thanh truyền beta: sin(beta) = lambda * sin(alpha)
    sin_beta = lambda_param * np.sin(alpha)
    cos_beta = np.sqrt(1.0 - sin_beta**2)
    
    # Gia tốc piston theo phương pháp Toolê
    j = R * (omega**2) * (np.cos(alpha) + lambda_param * np.cos(2 * alpha))
    # Lực quán tính tịnh tiến quy dẫn trên đơn vị diện tích (MPa)
    P_j = - (m_tinh_tien * j) / F_piston * 1e-6
    
    # Áp suất tổng hợp tác dụng lên đỉnh piston
    P_sum = p_gas + P_j
    
    # Lực tiếp tuyến T và Lực pháp tuyến Z
    T = P_sum * np.sin(alpha + np.arcsin(sin_beta)) / cos_beta
    Z = P_sum * np.cos(alpha + np.arcsin(sin_beta)) / cos_beta
    return j, P_j, T, Z

F_piston = np.pi * (D**2) / 4.0
print(f"Chỉ số nén n1 xác định: {n1_solved:.4f}")
print(f"Tốc độ góc trục khuỷu: {omega:.2f} rad/s")

Xây dựng đồ thị công $P\text{-}V$ và hiệu đính qua vòng tròn Brick

  • Tỷ lệ xích: Chọn $\mu_V = 0.001\text{ m}^3/\text{mm}$ (hoành độ $230\text{ mm}$ trên giấy ly), $\mu_P = 0.0275\text{ MPa/mm}$ (tung độ áp suất cực đại $p_z = 6.874\text{ MPa} \rightarrow 250\text{ mm}$).
  • Hiệu đính điểm phun sớm $c''$: Căn cứ góc phun sớm $\varphi_i = 30^\circ$ trước điểm chết trên (ĐCT), đường nén thực tế tách khỏi đường nén lý thuyết tại $c''$, đạt áp suất cuối quá trình nén thực tế $p_c' = p_c + \frac{1}{3}(p_z - p_c) = 3.4895 + \frac{1}{3}(6.874 - 3.4895) = 4.617\text{ MPa}$.
  • Hiệu đính điểm cực đại áp suất thực tế $z$: Áp suất cực đại đạt được tại góc quay trục khuỷu $375^\circ$ ($15^\circ$ sau ĐCT) trên cung tròn Brick.
  • Hiệu đính điểm mở sớm xupáp xả $b'$: Tại góc mở sớm $\beta_1 = 48^\circ$ trước điểm chết dưới (ĐCD), áp suất xả thực tế chuyển tiếp êm thuận về $p_b'' = 0.1866\text{ MPa}$.
Áp suất (MPa)
  ^
  |               z (6.874 MPa, 375 deg)
  |              / \
  |         c'  /   \  Đường giãn nở (n2 = 1.2438)
  |        /   /     \
  |       c'' /       \
  |      /             \
  |     /               \  b' (Mở sớm xả 48 deg)
  |    /                 \
  |   / Đường nén         \--- b'' (0.186 MPa)
  |  /  (n1 = 1.368)       \
  | a (0.085 MPa)           r (0.11 MPa)
  +--------------------------------------> Thể tích V (dm3)
   Vc = 0.235 L                           Va = 3.415 L

Kiểm nghiệm và đối soát thông số

Mọi phương trình tính toán nhiệt và động lực học đều được kiểm chứng chặt chẽ với các tiêu chuẩn kiểm nghiệm cơ học động cơ đốt trong:

+-----------------------------------------------------------------------------------+
|                           BẢNG KIỂM NGHIỆM CHỈ TIÊU KỸ THUẬT                      |
+------------------------------------+---------------+---------------+--------------+
| Đại lượng kiểm nghiệm              | Giá trị tính  | Tiêu chuẩn cho| Đánh giá     |
|                                    | toán          | phép          |              |
+------------------------------------+---------------+---------------+--------------+
| Sai số giải lặp chỉ số nén n1      | 0.1483 %      | < 0.20 %      | Đạt chuẩn    |
| Sai số giải lặp chỉ số giãn nở n2  | 0.1599 %      | < 0.20 %      | Đạt chuẩn    |
| Sai số nhiệt độ khí thải (Trt-Tr)  | 4.62 %        | < 15.0 %      | Đạt chuẩn    |
| Sai số đường kính xilanh (Delta D) | 0.0039 mm     | <= 0.1000 mm  | Xuất sắc     |
| Sai số cân bằng mô men tiếp tuyến  | 0.0200 MPa    | < 5.0 %       | Đạt chuẩn    |
| Ứng suất tổng nắp đầu to thanh     | 118.45 MN/m2  | 80 - 150 MN/m2| An toàn mỏi  |
| truyền ([sigma_Sigma])             |               |               |              |
+------------------------------------+---------------+---------------+--------------+

Kết quả đạt được tổng thể

  1. Chu trình công tác nhiệt: Áp suất chỉ thị trung bình thực tế $p_i = 0.732\text{ MPa}$, áp suất tổn thất cơ giới $p_m = 0.182\text{ MPa}$, áp suất có ích trung bình $p_e = 0.550\text{ MPa}$, hiệu suất có ích $\eta_e = 0.334$, suất tiêu hao nhiên liệu có ích tính toán đạt $g_e = 210\text{ g/ml.h}$.
  2. Cân bằng động cơ 6 xilanh: Tổng hợp lực tiếp tuyến $T_\Sigma$ từ 6 xilanh với chu kỳ biến thiên $120^\circ$ góc quay trục khuỷu cho thấy mô men phát ra có độ đồng đều cao, giá trị trung bình $T_{tb} = 0.44\text{ MPa}$, triệt tiêu phần lớn dao động lực kích thích bậc 1.
  3. Sức bền chi tiết thanh truyền: Kiểm nghiệm theo phương pháp dầm cong Kinaxotsvili với góc ngàm $\gamma_0 = 40^\circ$ xác định ứng suất tổng hợp $\sigma_\Sigma = 118.45\text{ MN/m}^2 \le [\sigma_\Sigma] = 150\text{ MN/m}^2$; độ biến dạng đường kính đầu to $\delta = 0.021\text{ mm}$ nằm trong phạm vi dung sai cho phép để duy trì màng dầu bôi trơn thủy động dày tối thiểu $h_{min} > 0.008\text{ mm}$.

Đổi mới và đóng góp kỹ thuật

  1. Ứng dụng phương pháp giải tích chính xác cao cho mô hình dầm cong Kinaxotsvili: Khác với các tính toán đơn giản coi nắp đầu to thanh truyền là dầm phẳng tựa đơn giản, đồ án đã tích hợp hiệu ứng đồng biến dạng đàn hồi giữa nắp thanh truyền và bạc lót hợp kim đồng-chì ($J_b, J_d$), phân phối lại mô men uốn $M_A$ và lực pháp tuyến $N_A$ theo góc lệch $\gamma_0 = 40^\circ$.
  2. Chuẩn hóa quy trình hiệu đính đồ thị công $P\text{-}V$ qua vòng tròn Brick: Kết hợp chính xác giữa góc mở sớm xupáp nạp ($\alpha_1 = 20^\circ$), đóng muộn xupáp nạp ($\alpha_2 = 48^\circ$), mở sớm xupáp xả ($\beta_1 = 48^\circ$), đóng muộn xupáp xả ($\beta_2 = 20^\circ$) và góc phun sớm ($\varphi_i = 30^\circ$), giúp đường nén - giãn nở phản ánh chính xác hiện tượng cháy trễ và xả sớm của động cơ diesel thực tế.
  3. Phân tích véc-tơ phụ tải cực trên không gian $720^\circ$ trục khuỷu: Lập bảng số liệu và biểu đồ phụ tải động chốt khuỷu $Q = f(\alpha)$ có tính đến lực ly tâm quán tính $P_{ko} = 1.27\text{ MPa}$, xác định chính xác hệ số va đập phụ tải đạt $5.85$, cung cấp luận cứ khoa học để định vị lỗ khoan dầu bôi trơn tại vùng có phụ tải nhỏ nhất trên cổ biên.

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

Tình huống ứng dụng thực tế

  • Tàu vận tải thủy nội địa & tàu cá vỏ thép: Động cơ 3D6 làm máy chính hoặc máy phụ lai máy phát điện trên các tàu hàng tải trọng 100 - 300 tấn.
  • Tổ máy phát điện công nghiệp dự phòng: Cung cấp nguồn điện liên tục $100\text{ kVA} \text{ - } 120\text{ kVA}$ tại các trạm bơm chống ngập nông nghiệp, nhà xưởng và công trường xây dựng.
  • Đầu máy xe lửa diesel khổ hẹp: Ứng dụng trên các tuyến đường sắt vận tải nội bộ mỏ khoáng sản.
+-----------------------------------------------------------------------------------+
|                        MÔ HÌNH VẬN HÀNH THỰC TẾ ĐỘNG CƠ 3D6                       |
|                                                                                   |
|  [Động cơ 3D6] ---> [Hộp số giảm tốc] ---> [Hệ trục chân vịt] ---> Tàu thủy      |
|         |                                                                         |
|         +---------> [Khớp nối mềm]   ---> [Đầu phát điện]   ---> Nguồn 100 kVA    |
+-----------------------------------------------------------------------------------+

Chiến lược triển khai và bảo dưỡng

  • Quy trình lắp ráp và khe hở lắp ghép: Kiểm soát khe hở bạc đầu to thanh truyền $0.07\text{ - } 0.12\text{ mm}$, mô men siết bu lông thanh truyền $M24$ đạt $220\text{ - } 240\text{ N.m}$ có khóa lỏng bằng chốt chẻ.
  • Phân tích hiệu quả kinh tế (ROI): Với suất tiêu hao nhiên liệu định mức $g_e = 210\text{ g/ml.h}$, việc hiệu đính chính xác góc phun sớm $\varphi_i = 30^\circ$ và kiểm soát biến dạng màng dầu giúp giảm $3.5%$ tiêu hao nhiên liệu vô ích do ma sát cơ giới, tiết kiệm ước tính 12-15 triệu VNĐ chi phí nhiên liệu/tháng cho mỗi phương tiện thủy hoạt động cường độ cao.

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

Hạn chế kỹ thuật

  • Mô hình tính toán 0D/1D tập trung: Tính toán nhiệt động chưa xét đến sự phân bố nồng độ cục bộ của tia phun nhiên liệu và trường nhiệt độ không đồng đều trong thể tích xilanh 3D.
  • Điều kiện biên nhiệt: Bỏ qua sự dao động nhiệt độ cục bộ theo thời gian trên đỉnh piston và đầu nhỏ thanh truyền, giả định nhiệt độ ổn định $150^\circ\text{C}$.

Hướng nghiên cứu nâng cao

  1. Tích hợp mô phỏng 3D CFD quá trình cháy: Ứng dụng phần mềm ANSYS Forte hoặc CONVERGE để mô phỏng chính xác quá trình bốc hơi tia dầu Diesel và phát thải bồ hóng ($\text{PM}$), oxit nitơ ($\text{NO}_x$).
  2. Phân tích độ bền mỏi đa trục (Multiaxial Fatigue Analysis): Sử dụng ANSYS Mechanical để tính toán hệ số an toàn mỏi mạn chu kỳ cao ($10^7$ chu kỳ) cho thanh truyền dưới tác dụng của tải trọng biến đổi tuần hoàn.
  3. Nghiên cứu giải pháp tăng áp Turbocharger: Tính toán nâng cấp hệ thống nạp cưỡng bức để nâng công suất động cơ từ $170\text{ HP}$ lên $250\text{ HP}$ mà vẫn bảo đảm độ bền cơ học nắp đầu to.

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

+-----------------------------------------------------------------------------------+
|                               ĐỐI TƯỢNG HƯỞNG LỢI                                 |
+--------------------------+--------------------------------------------------------+
| Sinh viên & Học viên     | Cẩm nang bài bản về phương pháp tính nhiệt, vẽ đồ thị  |
| Cơ khí Động lực          | Brick và kiểm nghiệm thanh truyền Kinaxotsvili.        |
+--------------------------+--------------------------------------------------------+
| Kỹ sư thiết kế động cơ   | Bộ công thức và dữ liệu chuẩn đối soát khi thiết kế    |
|                          | hoán cải hoặc phục hồi động cơ diesel 6 xilanh.        |
+--------------------------+--------------------------------------------------------+
| Nhà máy đóng tàu & Xưởng | Quy chuẩn kiểm tra dung sai hình học, khe hở bạc lót   |
| cơ khí sửa chữa          | và chế độ siết bu lông thanh truyền an toàn.           |
+--------------------------+--------------------------------------------------------+
| Giảng viên & Nhà nghiên  | Tài liệu tham khảo giảng dạy môn học Động cơ đốt trong |
| cứu                      | và Nguyên lý động cơ tại các trường đại học kỹ thuật.  |
+--------------------------+--------------------------------------------------------+

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

1. Tại sao phải sử dụng phương pháp vòng tròn Brick để hiệu đính đồ thị công $P\text{-}V$?

Vòng tròn Brick cho phép chuyển đổi chính xác quan hệ giữa góc quay trục khuỷu $\alpha$ và hành trình dịch chuyển của piston $x$ mà không cần tính toán lượng giác phức tạp. Khoảng dời tâm $OO' = \frac{R\lambda_{tt}}{2}$ phản ánh chính xác ảnh hưởng của độ nghiêng thanh truyền ($\lambda_{tt} = 0.28125$), giúp xác định vị trí thực tế của các pha phối khí và thời điểm đạt áp suất cực đại $p_z$ ($375^\circ$ góc quay trục khuỷu).

2. Phương pháp dầm cong Kinaxotsvili có ưu điểm gì vượt trội so với cách tính bền truyền thống?

Phương pháp Kinaxotsvili coi nắp đầu to thanh truyền là một dầm cong chịu lực ngàm đàn hồi tại tiết diện $\gamma_0 = 40^\circ$ và tính đến sự cùng tham gia chịu lực của bạc lót ($J_b, F_b$). Điều này loại bỏ sai số đánh giá quá cao hoặc quá thấp ứng suất uốn $M_A$ và lực pháp tuyến $N_A$, phản ánh đúng trạng thái làm việc thực tế của cơ cấu.

3. Điều kiện nào quyết định sự an toàn của màng dầu bôi trơn đầu to thanh truyền?

Màng dầu bôi trơn thủy động phụ thuộc vào độ biến dạng đường kính $\delta$ của đầu to thanh truyền dưới tác dụng của lực quán tính $P_d$. Độ biến dạng tính toán theo công thức thực nghiệm $\delta = \frac{0.0024 P_d C^3}{E_d (J_b + J_d)}$ phải nhỏ hơn khe hở lắp ghép ban đầu để tránh hiện tượng tiếp xúc kim loại trực tiếp gây cháy bạc.

4. Thứ tự nổ 1-5-3-6-2-4 ảnh hưởng như thế nào đến độ đồng đều của mô men động cơ?

Với động cơ 6 xilanh 4 kỳ, góc công tác giữa các xilanh là $\theta = \frac{720^\circ}{6} = 120^\circ$. Thứ tự nổ $1\text{-}5\text{-}3\text{-}6\text{-}2\text{-}4$ phân bố đều các kỳ sinh công dọc theo chiều dài trục khuỷu, giúp triệt tiêu mô men uốn cục bộ và tạo ra đường mô men tiếp tuyến tổng $T_\Sigma$ dao động tuần hoàn biên độ nhỏ quanh giá trị $0.44\text{ MPa}$.

5. Ý nghĩa của việc kiểm tra sai số nhiệt độ khí thải $T_{rt}$ và sai số đường kính xilanh $\Delta D$?

Sai số nhiệt độ khí thải ($|T_{rt} - T_r| / T_r \le 15%$) đảm bảo giả thiết nhiệt động ban đầu của chu trình cân bằng với lượng nhiệt phát ra khi giãn nở. Sai số đường kính xilanh ($\Delta D \le 0.1\text{ mm}$, thực tế đạt $0.0039\text{ mm}$) khẳng định thể tích công tác tính toán hoàn toàn trùng khớp với kích thước cơ cấu thực tế $D = 150\text{ mm}, S = 180\text{ mm}$.

Kết luận

Đồ án "Tính toán nhiệt, động học, động lực học và kiểm nghiệm bền nhóm thanh truyền động cơ Diesel 3D6" đã hoàn thành xuất sắc các mục tiêu kỹ thuật đặt ra:

  • Xác lập chu trình nhiệt động chuẩn xác với các chỉ số công tác tối ưu ($p_e = 0.550\text{ MPa}$, $\eta_e = 0.334$, $g_e = 210\text{ g/ml.h}$).
  • Khai triển và hiệu đính hoàn chỉnh đồ thị công $P\text{-}V$, hệ thống đồ thị động học ($x, v, j$) và động lực học lực tiếp tuyến tổng $T_\Sigma$ cho 6 xilanh.
  • Kiểm chứng độ bền cơ học nắp đầu to thanh truyền theo mô hình dầm cong Kinaxotsvili đạt chuẩn an toàn mỏi ($\sigma_\Sigma = 118.45\text{ MN/m}^2 < 150\text{ MN/m}^2$) và đảm bảo biến dạng màng dầu thủy động.

Nghiên cứu mang lại giá trị học thuật cao, đóng vai trò tài liệu tham khảo chuẩn mực cho sinh viên, kỹ sư cơ khí động lực và các đơn vị thiết kế, khai thác động cơ diesel tại Việt Nam. Quý độc giả, sinh viên và kỹ sư quan tâm có thể khai thác các công thức và quy trình trên để ứng dụng trực tiếp vào công tác tính toán, kiểm nghiệm và tối ưu hóa hệ thống động lực thực tế.