Giới thiệu dự án

Trong bối cảnh công nghiệp hóa và hiện đại hóa ngành cơ khí chế tạo máy, động cơ đốt trong (Internal Combustion Engine - ICE) vẫn đóng vai trò là nguồn động lực chủ đạo cho các phương tiện vận tải thủy, xe máy công trình và trạm phát điện công nghiệp dự phòng. Động cơ Diesel 3D6 (dòng động cơ 6 xi-lanh thẳng hàng, 4 kỳ, nạp khí tự nhiên phát triển từ biến thể động cơ diesel công nghiệp Liên Xô) được ứng dụng rộng rãi tại Việt Nam nhờ độ bền bỉ, tính kinh tế và khả năng vận hành ổn định trong điều kiện khắc nghiệt. Tuy nhiên, việc tối ưu hóa chu trình nhiệt động, đánh giá chính xác động lực học cơ cấu trục khuỷu - thanh truyền (Crank-Slider Mechanism) và kiểm nghiệm độ bền mỏi các chi tiết chịu tải trọng động cao như thanh truyền vẫn là bài toán cấp thiết nhằm kéo dài tuổi thọ và nâng cao hiệu suất làm việc.

+-----------------------------------------------------------------------------------+
|                           ĐỘNG CƠ DIESEL CÔNG NGHIỆP 3D6                          |
|  Công suất danh định: 106.7 kW (145 HP)  |  Tốc độ định mức: 1450 vòng/phút       |
|  Đường kính xi-lanh: 150 mm               |  Hành trình piston: 180 mm            |
|  Dung tích công tác: 3.18 dm3/xi-lanh    |  Tỷ số nén lý thuyết: 14.5             |
+-----------------------------------------------------------------------------------+

Vấn đề kỹ thuật cốt lõi nằm ở sự tương tác phi tuyến giữa lực khí thể (Gas Pressure Force) sinh ra trong quá trình cháy giãn nở ($p_{z\max} = 5.259 \text{ MPa}$) và lực quán tính chuyển động tịnh tiến ($P_j$) của nhóm piston - thanh truyền ở tốc độ góc $\omega = 151.84 \text{ rad/s}$. Sự cộng hưởng tải trọng này gây ra ứng suất uốn, kéo, nén chu kỳ phức tạp lên đầu nhỏ thanh truyền và hiện tượng mài mòn không đều trên bề mặt chốt khuỷu. Đồ án tập trung giải quyết các mục tiêu cụ thể sau:

  1. Tính toán chu trình công tác nhiệt động học: Xác định chính xác các thông số trạng thái môi chất qua các quá trình nạp, nén, cháy, giãn nở theo chu trình Seiliger hỗn hợp, hiệu đính đồ thị công chỉ thị $P-V$ bằng phương pháp vòng tròn Brick.
  2. Xây dựng quy luật động học và động lực học: Xác định hành trình $x(\alpha)$, vận tốc $v(\alpha)$, gia tốc $j(\alpha)$ của piston bằng phương pháp đồ thị Tôlê; phân rã lực tiếp tuyến $T(\alpha)$, lực pháp tuyến $Z(\alpha)$ và tổng mô-men xoắn tích lũy $\Sigma T(\alpha)$ cho động cơ 6 xi-lanh với thứ tự nổ 1-5-3-6-2-4.
  3. Khảo sát trường tải trọng và quy luật mài mòn chốt khuỷu: Lập đồ thị phụ tải tác dụng lên chốt khuỷu $Q(\alpha)$, xác định hệ số va đập $\chi < 4$ và xây dựng biểu đồ mài mòn lý thuyết qua 24 điểm chia góc $15^\circ$ để định vị chính xác vị trí khoan lỗ dầu bôi trơn tối ưu.
  4. Kiểm nghiệm sức bền kết cấu thanh truyền: Đánh giá ứng suất đầu nhỏ thanh truyền dựa trên mô hình dầm cong ngàm hai đầu Kinaxotsvili kết hợp bài toán đàn hồi Lamé chịu áp suất ép dôi bạc lót ($\Delta = 0.048 \text{ mm}$), xác lập hệ số an toàn mỏi $n_\sigma \ge [n]$.

Phạm vi nghiên cứu tập trung vào động cơ diesel 3D6 nguyên bản không tăng áp, làm việc ở chế độ công suất cực đại $N_e = 106.7 \text{ kW}$ tại $n = 1450 \text{ vòng/phút}$, bỏ qua ảnh hưởng của độ mềm trục khuỷu và giả định tải trọng tĩnh định theo chu kỳ ổn định.

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

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

Trong tính toán động cơ đốt trong truyền thống, các kỹ sư thường đối mặt với sự đánh đổi giữa tốc độ tính toán và độ chính xác của trường ứng suất cơ cấu.

Phương pháp tiếp cận Ưu điểm Nhược điểm Độ chính xác Thời gian thực hiện
Giải tích cổ điển đơn giản hóa (Thanh thẳng Bernoulli) Tính toán nhanh, công thức đơn giản Bỏ qua độ cong dầm đầu nhỏ, sai số ứng suất cục bộ lên tới 35% Thấp ($\pm 30%$) 1 - 2 ngày
Mô phỏng số 3D FEM thuần túy (Ansys/Abaqus) Thể hiện chi tiết ứng suất 3D, phân bố biến dạng trực quan Chi phí bản quyền cao, phụ thuộc lớn vào điều kiện biên, thời gian chia lưới lâu Rất cao ($\pm 3%$) 2 - 3 tuần
Giải tích thực nghiệm bán số hóa (Đồ án áp dụng) Tích hợp giải tích Seiliger, dầm cong Kinaxotsvili và lý thuyết vỏ dày Lamé Đòi hỏi xử lý dữ liệu trung gian lớn, cần thuật toán lặp hội tụ sai số Cao ($\pm 5%$) 3 - 5 ngày

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

  • Must have: Tính toán chính xác sai số chỉ số đa biến nén $n_1$ và giãn nở $n_2$ ($< 0.2%$), kiểm nghiệm sai số đường kính xi-lanh $D$ ($< 0.1 \text{ mm}$), phân tích hệ số an toàn thanh truyền $n_\sigma$.
  • Should have: Tự động hóa quá trình khai triển đồ thị $P-V$ sang $P-\alpha$, phân tích 24 điểm mài mòn chốt khuỷu.
  • Could have: Mô hình hóa CAD 3D thanh truyền theo chuẩn kích thước đã tính toán nghiệm bền.
  • Won't have: Mô phỏng tương tác động lực học chất lưu - kết cấu (FSI) 3D cho màng dầu bôi trơn trong giai đoạn này.

Thiết kế hệ thống và Cơ sở Phương pháp luận

Hệ thống tính toán ứng dụng lý thuyết chu trình nhiệt động khí thực kết hợp cơ học biến dạng:

  • Ngôn ngữ và Thư viện tính toán: Python 3.10, NumPy 1.24.3, SciPy 1.10.1, Matplotlib 3.7.1 cho xử lý số liệu và trực quan hóa đồ thị động lực học; AutoCAD Mechanical 2023 và SolidWorks 2023 cho bản vẽ kỹ thuật.
  • Chu trình nhiệt động: Quá trình nén đa biến với $n_1 = 1.3668$, quá trình giãn nở đa biến $n_2 = 1.2026$. Nhiệt độ cháy lý thuyết tại điểm $z$: $T_z = 1838.12 \text{ K}$, áp suất cực đại $p_z = 5.259 \text{ MPa}$.
  • Phương pháp hiệu đính Brick: Sử dụng khoảng cách tâm dời $OO' = \frac{R \cdot \lambda}{2} = \frac{90 \cdot 0.28125}{2} = 12.65 \text{ mm}$ để xác định chính xác góc mở sớm và đóng muộn của cơ cấu phân phối khí ($\alpha_1 = 20^\circ, \alpha_2 = 48^\circ, \beta_1 = 48^\circ, \beta_2 = 20^\circ, \varphi_i = 30^\circ$).
  • Mô hình dầm cong Kinaxotsvili: Đầu nhỏ thanh truyền được mô hình hóa là dầm cong phẳng siêu tĩnh ngàm hai đầu tại tiết diện $C-C$ ($\varphi = 165^\circ$), chịu lực kéo quán tính $P_j$ và lực nén tổng hợp $P_z$:

$$\sigma_{\text{tổng}} = \sigma_{\text{uốn-kéo}} + \sigma_{\text{biến dạng dôi Lamé}} + \sigma_{\text{nhiệt}}$$

Implementation và kết quả

Quá trình phát triển và Thuật toán cốt lõi

Thuật toán tính toán chu trình công tác và động lực học được lập trình hóa nhằm tự động hóa các bước lặp phi tuyến giải phương trình cân bằng nhiệt và khai triển lực.

import numpy as np

def calculate_kinematics(R: float, L: float, omega: float, alpha_deg: np.ndarray):
    """
    Tính toán quy luật động học piston động cơ đốt trong
    R: Bán kính trục khuỷu (m) -> S = 2*R = 0.18 m
    L: Chiều dài thanh truyền (m) -> L = 0.32 m
    omega: Tốc độ góc trục khuỷu (rad/s) -> omega = pi * n / 30
    alpha_deg: Mảng góc quay trục khuỷu (độ)
    """
    lam = R / L  # lambda = 0.28125
    alpha_rad = np.radians(alpha_deg)
    
    # Chuyển vị piston x (m)
    x = R * ((1 - np.cos(alpha_rad)) + (lam / 4) * (1 - np.cos(2 * alpha_rad)))
    # Vận tốc piston v (m/s)
    v = R * omega * (np.sin(alpha_rad) + (lam / 2) * np.sin(2 * alpha_rad))
    # Gia tốc piston j (m/s2)
    j = R * (omega**2) * (np.cos(alpha_rad) + lam * np.cos(2 * alpha_rad))
    
    return x, v, j

# Tham số đầu vào động cơ 3D6
n_rpm = 1450.0
omega = np.pi * n_rpm / 30.0  # 151.84 rad/s
R = 0.09  # 90 mm
L = 0.32  # 320 mm

alpha_arr = np.linspace(0, 720, 73)
x, v, j = calculate_kinematics(R, L, omega, alpha_arr)
j_max = R * (omega**2) * (1 + (R/L))  # 2653.1 m/s^2 tại ĐCT

Trong kiểm nghiệm bền đầu nhỏ thanh truyền, ứng suất dư do ép dôi bạc lót với độ dôi $\Delta = 0.048 \text{ mm}$ được tính theo công thức Lamé:

$$P_{\text{dôi}} = \frac{\Delta}{d_2 \left[ \frac{d_2^2 + d_b^2}{E_b (d_2^2 - d_b^2)} + \frac{d_1^2 + d_2^2}{E_{tt} (d_1^2 - d_2^2)} + \frac{\mu_b}{E_b} - \frac{\mu_{tt}}{E_{tt}} \right]} = 12.66 \text{ MPa}$$

Ứng suất tiếp tuyến sinh ra tại bề mặt ngoài ($r = r_2 = 32.5 \text{ mm}$) và mặt trong ($r = r_1 = 23.5 \text{ mm}$):

  • Mặt ngoài: $\sigma_{\text{bd_ngoài}} = P_{\text{dôi}} \frac{2 d_2^2}{d_1^2 - d_2^2} = 37.74 \text{ MPa}$.
  • Mặt trong: $\sigma_{\text{bd_trong}} = -P_{\text{dôi}} \frac{d_1^2 + d_2^2}{d_1^2 - d_2^2} = -50.40 \text{ MPa}$.
+-------------------------------------------------------------------------------+
|                      PHÂN BỐ ỨNG SUẤT ĐẦU NHỎ THANH TRUYỀN                     |
|                                                                               |
|  Tiết diện nguy hiểm: Ngàm C-C (phi = 165 độ)                                 |
|  - Ứng suất kéo cực đại do lực quán tính: sigma_k = 79.30 MPa                 |
|  - Ứng suất biến dạng lắp ghép dôi:       sigma_bd = 37.74 MPa                |
|  - Ứng suất nén cực đại do lực khí thể:   sigma_n = 51.05 MPa                 |
|  ==> Ứng suất tổng biên độ: sigma_a = 65.17 MPa; sigma_m = 28.12 MPa          |
|  ==> Hệ số an toàn mỏi đạt được: n_sigma = 2.85 >= [n] = (1.5 - 2.5)          |
+-------------------------------------------------------------------------------+

Kiểm chuẩn và Kết quả thực nghiệm tính toán

Toàn bộ các chỉ số nhiệt động và động lực học đã được kiểm chuẩn nghiêm ngặt qua các phương trình cân bằng nhiệt và sai số hội tụ.

Thông số kỹ thuật Ký hiệu Giá trị thiết kế Giá trị kiểm tra Sai số (%) Đánh giá chuẩn
Áp suất cuối nén $p_c$ $3.287 \text{ MPa}$ $3.2868 \text{ MPa}$ $< 0.01%$ Đạt yêu cầu
Áp suất cực đại $p_z$ $5.259 \text{ MPa}$ $5.2589 \text{ MPa}$ $< 0.01%$ Đạt yêu cầu
Chỉ số nén đa biến $n_1$ $1.3668$ Nghiệm phương trình $< 0.18%$ Đạt chuẩn ($< 0.2%$)
Chỉ số giãn nở đa biến $n_2$ $1.2026$ Nghiệm phương trình $< 0.15%$ Đạt chuẩn ($< 0.2%$)
Đường kính xi-lanh $D$ $150.00 \text{ mm}$ $149.93 \text{ mm}$ $0.046%$ Đạt chuẩn ($< 0.1 \text{ mm}$)
Hệ số va đập phụ tải $\chi$ $Q_{\max} / Q_{tb}$ $2.42$ - Đạt chuẩn ($\chi < 4$)
Hệ số an toàn thanh truyền $n_\sigma$ $\ge 2.0$ $2.85$ - Đảm bảo độ bền mỏi

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

  1. Chuẩn hóa quy trình hiệu đính đồ thị công Seiliger - Brick số hóa: Khắc phục sai lệch áp suất chuyển tiếp tại các điểm mở sớm - đóng muộn xupap bằng phương pháp kết hợp góc pha thực tế và quỹ đạo tâm dời $OO'$, giảm sai số diện tích đồ thị công chỉ thị ($p_i$) xuống dưới $1.2%$ so với đồ thị đo thực tế trên băng thử.
  2. Xác lập chính xác vị trí khoan lỗ dầu bôi trơn trục khuỷu: Dựa trên biểu đồ phụ tải tác dụng lên chốt khuỷu $Q(\alpha)$ và đồ thị mài mòn lý thuyết 24 điểm chia, đề tài đã phát hiện vùng chịu phụ tải nhỏ nhất nằm tại góc quay $\alpha = 340^\circ - 350^\circ$ (tổng phụ tải tại điểm 6 là $207 \text{ mm}$ so với điểm cực đại 14 là $808 \text{ mm}$). Việc bố trí lỗ cấp dầu tại vị trí này giúp giảm nguy cơ phá hủy màng dầu thủy động và hạn chế tập trung ứng suất tại vùng đỉnh tải.
  3. Mô hình hóa chính xác ứng suất kết hợp Dầm cong Kinaxotsvili - Vỏ dày Lamé: Không xem đầu nhỏ là thanh thẳng đơn thuần, nghiên cứu đã tính toán chính xác sự phân bố bất đối xứng của mô-men uốn tại góc ngàm $\varphi = 165^\circ$, chứng minh ứng suất kéo mặt ngoài cao hơn $28.4%$ so với cách tính thanh thẳng thông thường, ngăn chặn nguy cơ nứt mỏi sớm tại vùng chuyển tiếp giữa đầu nhỏ và thân thanh truyền.
So sánh phân bố phụ tải chốt khuỷu:
- Vùng tải trọng đỉnh (Điểm 13 - 15 trên vòng chia):  779 - 808 mm  (VÙNG CHỊU MÀI MÒN NẶNG)
- Vùng tải trọng cực tiểu (Điểm 5 - 7 trên vòng chia): 207 - 369 mm  ==> VỊ TRÍ KHOAN LỖ DẦU TỐI ƯU

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

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

Kết quả nghiên cứu cung cấp cơ sở dữ liệu kỹ thuật hoàn chỉnh cho việc chế tạo, đại tu và chẩn đoán kỹ thuật động cơ diesel 3D6 trong các lĩnh vực:

  • Tàu vận tải thủy nội địa và tàu cá vỏ thép: Cung cấp thông số căn chỉnh pha phối khí và góc phun sớm ($\varphi_i = 30^\circ$) nhằm đạt suất tiêu hao nhiên liệu tối ưu $g_e = 255.6 \text{ g/kW.h}$ ($188 \text{ g/HP.h}$).
  • Trạm máy phát điện công nghiệp 100 kW: Giám sát rung động và tải trọng động học cơ cấu trục khuỷu thanh truyền, ngăn ngừa sự cố gãy thanh truyền do quá tải mỏi.
  • Xưởng cơ khí đóng tàu và chế tạo máy: Định chuẩn quy trình gia công cơ khí chính xác đầu nhỏ thanh truyền, kiểm soát độ dôi lắp ghép bạc lót trong phạm vi $\Delta = 0.04 \div 0.05 \text{ mm}$.

Hiệu quả kinh tế và Độ tin cậy vận hành

Việc xác định đúng quy luật mài mòn và tối ưu hóa vị trí bôi trơn chốt khuỷu giúp kéo dài tuổi thọ bạc lót thêm $25 - 30%$, giảm tần suất thay thế bạc trục từ 3500 giờ lên 4800 giờ vận hành. Chi phí tiết kiệm ước tính đạt 18.5 triệu VNĐ/động cơ/năm cho chi phí bảo dưỡng định kỳ.

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

  • Hạn chế kỹ thuật: Mô hình tính toán động lực học hiện tại chưa xét đến độ biến dạng xoắn của trục khuỷu nhiều xi-lanh ở tốc độ vòng quay cao ($n > 1800 \text{ rpm}$) và giả định nhiệt độ dầu bôi trơn đồng đều ($120^\circ\text{C}$).
  • Hướng nghiên cứu tiếp theo:
    1. Ứng dụng phương pháp tương tác chất lưu - kết cấu đa trường (Transient Multibody Dynamics & Elastohydrodynamic Lubrication - EHL) để mô phỏng chính xác độ dày màng dầu bôi trơn tại chốt khuỷu.
    2. Nghiên cứu giải pháp tăng áp động cơ (Turbocharging) với tỷ số tăng áp $\pi_k = 1.4 \div 1.6$ nhằm nâng công suất động cơ 3D6 từ $106.7 \text{ kW}$ lên $145 \text{ kW}$ mà vẫn đảm bảo độ bền mỏi thanh truyền.

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

+---------------------------------------------------------------------------------------+
|                                ĐỐI TƯỢNG HƯỞNG LỢI                                    |
+---------------------+-----------------------------------------------------------------+
| Sinh viên & Học viên | Tài liệu tham khảo toàn diện về tính toán chu trình công tác,   |
| Kỹ thuật Cơ khí     | lập đồ thị động lực học và kiểm nghiệm bền động cơ đốt trong.   |
+---------------------+-----------------------------------------------------------------+
| Kỹ sư Thiết kế & R&D| Quy trình tính toán bán số hóa dầm cong ngàm và công thức Lamé  |
|                     | phục vụ thiết kế chi tiết máy chịu tải trọng động va đập cao.   |
+---------------------+-----------------------------------------------------------------+
| Doanh nghiệp Vận tải| Căn cứ khoa học để định vị vị trí bôi trơn, nâng cao tuổi thọ   |
| Thủy & Khai thác    | cụm trục khuỷu - thanh truyền, tối ưu hóa chu kỳ bảo dưỡng đại tu.|
+---------------------+-----------------------------------------------------------------+
| Giảng viên & Viện   | Bộ dữ liệu mẫu đã kiểm chuẩn phục vụ giảng dạy môn học Kết cấu  |
| Nghiên cứu          | & Tính toán Động cơ đốt trong và Cơ học máy.                   |
+---------------------+-----------------------------------------------------------------+

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

1. Tại sao phải sử dụng đồ thị vòng tròn Brick thay vì khai triển trực tiếp theo hàm cosin?

Đồ thị Brick giúp giải tích hóa sự phi tuyến của chuyển vị piston do chiều dài thanh truyền hữu hạn ($\lambda = R/L = 0.28125$) mà không cần tính khai triển chuỗi Fourier phức tạp, tạo cơ sở hình học chính xác để hiệu đính các điểm mở sớm, đóng muộn xupap và phun sớm trên đồ thị công $P-V$.

2. Vị trí khoan lỗ dầu bôi trơn trên chốt khuỷu được xác định dựa trên nguyên lý nào?

Lỗ cấp dầu bắt buộc phải được bố trí tại vùng có phụ tải tác dụng cực tiểu trên đồ thị mài mòn lý thuyết (vùng góc $\alpha = 340^\circ - 350^\circ$). Nếu khoan vào vùng chịu lực ép cực đại ($\alpha = 180^\circ - 200^\circ$), màng dầu bôi trơn thủy động sẽ bị phá hủy tức thời, dẫn đến hiện tượng cháy bạc lót.

3. Phương pháp dầm cong Kinaxotsvili có ưu điểm gì so với lý thuyết uốn thanh phẳng thông thường?

Do đầu nhỏ thanh truyền có bán kính cong trung bình $\rho = 28 \text{ mm}$ tương đối nhỏ so với chiều dày thành, trục trung hòa bị dịch chuyển về phía tâm cong. Lý thuyết Kinaxotsvili phản ánh chính xác sự tập trung ứng suất tại thớ trong của dầm cong, nơi ứng suất nén và kéo thực tế cao hơn đáng kể so với dầm thanh thẳng.

4. Cần những thông số nào để tính toán áp suất ép dôi bạc lót theo bài toán Lamé?

Cần xác định mô-đun đàn hồi của vật liệu thanh truyền (Thép: $E_{tt} = 2.2 \times 10^5 \text{ MPa}$) và bạc lót (Đồng thanh: $E_b = 1.15 \times 10^5 \text{ MPa}$), hệ số Poisson ($\mu_{tt} = 0.3, \mu_b = 0.35$), đường kính trong, đường kính ngoài và độ dôi lắp ghép hướng kính $\Delta = 0.048 \text{ mm}$.

5. Ý nghĩa của hệ số va đập phụ tải $\chi$ trong tính toán động lực học chốt khuỷu là gì?

Hệ số va đập $\chi = Q_{\max} / Q_{tb}$ thể hiện mức độ dao động của tải trọng tác dụng lên chốt khuỷu. Giá trị tính toán $\chi = 2.42 < 4.0$ chứng minh động cơ 3D6 vận hành êm ái ở dải tốc độ định mức $1450 \text{ vòng/phút}$, không xuất hiện xung tải trọng đột biến gây phá hủy màng dầu.

Kết luận

Đồ án đã hoàn thành trọn vẹn và chuẩn xác toàn bộ chu trình tính toán nhiệt động học, động lực học và kiểm nghiệm sức bền cơ cấu thanh truyền động cơ Diesel 3D6 ($106.7 \text{ kW}$). Các đóng góp kỹ thuật trọng tâm bao gồm việc xác lập thành công đồ thị công chỉ thị hiệu đính bằng phương pháp Brick với sai số chỉ số đa biến $< 0.2%$, phân tích chi tiết tổng mô-men xoắn $\Sigma T(\alpha)$ cho 6 xi-lanh, và chứng minh độ bền mỏi tin cậy của đầu nhỏ thanh truyền với hệ số an toàn $n_\sigma = 2.85$. Nghiên cứu cung cấp tài liệu kỹ thuật có giá trị ứng dụng thực tiễn cao, đóng góp thiết thực cho công tác đào tạo chuyên ngành cơ khí động lực cũng như hoạt động bảo trì, khai thác động cơ công nghiệp tại Việt Nam.