Giới thiệu dự án

Trong bối cảnh ngành công nghiệp cơ khí động lực và sản xuất ô tô - máy công trình đang chịu áp lực lớn về tiêu chuẩn khí thải và hiệu suất kinh tế, việc nâng cao chất lượng làm việc của động cơ đốt trong là nhiệm vụ sống còn. Theo các nghiên cứu ngành động cơ, tổn thất trao đổi khí (nạp - thải) chiếm từ 15% đến 25% tổng tổn thất cơ giới ở các dải tốc độ làm việc cao. Hệ thống phân phối khí (valvetrain) đóng vai trò điều khiển chính xác toàn bộ quá trình nạp đầy khí nạp mới và thải sạch sản phẩm cháy ra khỏi xilanh.

Động cơ Diesel D240 (công suất định mức 80 mã lực tại 2200 - 2500 vòng/phút, tỷ số nén 17.5:1) là dòng động cơ diesel 4 kỳ được ứng dụng rộng rãi trên các phương tiện vận tải và máy kéo nông nghiệp tại Việt Nam. Tuy nhiên, qua quá trình sử dụng thực tế và sản xuất phụ tùng thay thế, cam nguyên thủy của động cơ D240 bộc lộ nhiều điểm hạn chế: biên dạng cam cũ được tính toán bằng phương pháp đồ giải truyền thống khiến vận tốc dòng khí nạp qua khe hở nấm xupáp đạt tới $140.175\text{ m/s}$ (vượt ngưỡng cho phép $111.5 - 128.65\text{ m/s}$), ứng suất tiếp xúc bề mặt cam - con đội lên tới $678.717\text{ N/mm}^2$, gây hiện tượng tổn thất lưu động, mài mòn nhanh và suy giảm hệ số nạp.

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

  1. Đánh giá toàn diện chất lượng nạp - thải và động lực học cơ cấu phân phối khí cam nguyên thủy động cơ D240.
  2. Ứng dụng phần mềm chuyên dụng AVL TYCON để mô hình hóa, tính toán trị số thời gian - tiết diện ($A$) và vận tốc dòng khí.
  3. Thiết kế lại biên dạng cam nạp mới bằng phương pháp đa thức Polydyne nhằm mở rộng trị số thời gian - tiết diện, kiểm soát gia tốc đóng/mở xupáp.
  4. Giảm ứng suất tiếp xúc bề mặt cam - con đội xuống dưới $450\text{ N/mm}^2$ (mức cho phép tối đa $900\text{ N/mm}^2$) và đưa vận tốc dòng khí về dải tiêu chuẩn.

Phương pháp tiếp cận kết hợp giữa giải tích toán học đa thức bậc cao và mô phỏng động lực học đa vật thể (MBD) trên phần mềm TYCON giúp rút ngắn chu kỳ thiết kế, dự báo chính xác hiện tượng va đập, rung ồn mà không cần chế tạo thử nghiệm tốn kém.

Phạm vi nghiên cứu tập trung vào cơ cấu phân phối khí xupáp treo dẫn động qua con đội - đũa đẩy - đòn bẩy của động cơ D240 tại điều kiện làm việc định mức, không can thiệp thay đổi kích thước buồng cháy hay hành trình piston nguyên bản.

+-------------------------------------------------------------------------+
|                  CƠ CẤU PHÂN PHỐI KHÍ ĐỘNG CƠ D240                      |
|                                                                         |
|  [ Trục cam D240 ] ---> [ Con đội ] ---> [ Đũa đẩy ] ---> [ Đòn bẩy ]   |
|         |                                                      |        |
|    (Biên dạng)                                            (Tỷ số 1.59)  |
|         v                                                      v        |
|  [ Động học cam ] <=== Tối ưu Polydyne (TYCON) ===> [ Xupáp Nạp/Thải ]  |
+-------------------------------------------------------------------------+

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

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

Trước khi thực hiện cải biến, việc đánh giá các giải pháp thiết kế cam truyền thống được tổng hợp qua bảng phân tích sau:

Phương pháp thiết kế Ưu điểm Nhược điểm Đánh giá tính khả thi
Đồ giải truyền thống (Cam 2-3 cung) Biên dạng hình học đơn giản, dễ gia công bằng máy công cụ cơ. Gia tốc biến thiên gián đoạn, gây va đập động lực học lớn, sai số cao. Lạc hậu, không tối ưu được khí động học.
Phương pháp gia tốc ổn định (STAC) Hạn chế gia tốc đỉnh, giảm lực va đập lò xo xupáp. Trị số thời gian - tiết diện không đạt cực đại, lưu lượng khí hạn chế. Trung bình, phù hợp động cơ tốc độ thấp.
Phương pháp Polydyne (AVL TYCON) Khử biến dạng đàn hồi hệ thống, tối đa hóa thời gian - tiết diện, gia tốc liên tục. Đòi hỏi thuật toán giải phức tạp và máy gia công điều khiển số (CNC). Tối ưu nhất cho động cơ Diesel hiện đại.

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

  • Must have: Đảm bảo độ nâng xupáp lớn nhất ($L = 11.6464\text{ mm}$), góc nâng cam chính $58^\circ$, ứng suất tiếp xúc $\sigma_{max} < 900\text{ N/mm}^2$.
  • Should have: Gia tốc dương $a_{max}^+ \le 2500\text{ m/s}^2$, gia tốc âm $|a_{max}^-| \le 800\text{ m/s}^2$.
  • Could have: Tích hợp đoạn giảm chấn Ramp dài $19.5^\circ$ với độ nâng $0.4755\text{ mm}$ để triệt tiêu tiếng gõ khe hở nhiệt.
  • Won't have: Thay đổi kết cấu hệ thống truyền động bánh răng trục cam nguyên bản.

Thiết kế hệ thống

Mô hình động học cơ cấu phân phối khí được cấu trúc từ các phần tử tiêu chuẩn trong phần tử hữu hạn của module AVL TYCON:

[CDAT: Cam Contour] ---> [SRBS: Slider Bearing] ---> [MASS: Pushrod/Lifter]
                                                             |
[SPGE: Spring Model] <--- [VAFA: Valve Face] <--- [ROAR: Rocker Arm (1.59)]

Thông số công nghệ và ràng buộc kết cấu:

  • Bán kính vòng cơ sở cam ($R_0$): $20.63\text{ mm}$
  • Chiều rộng tiếp xúc bề mặt: $22\text{ mm}$
  • Đường kính nấm xupáp nạp: $d_n = 52\text{ mm}$, góc nghiêng đế $\alpha = 30^\circ$
  • Đường kính họng xupáp nạp: $d_h = 42\text{ mm}$, đường kính thân $d_t = 10\text{ mm}$
  • Tỷ số truyền đòn bẩy: $i_{db} = 1.59$

Methodology

Phương pháp nghiên cứu tích hợp quy trình tính toán động lực học phi tuyến theo thuật toán Predictor - Corrector dựa trên hệ phương trình vi phân chuyển động tổng quát Newton - Euler:

$$M \cdot \ddot{q} = F_{\text{Ngoại lực}} + F_{\text{Liên kết}}$$

Trong đó:

  • $M$: Ma trận khối lượng quy dẫn của các chi tiết (con đội, đũa đẩy, đòn bẩy, xupáp).
  • $\ddot{q}$: Vector gia tốc chuyển vị của các phần tử.
  • $F_{\text{Ngoại lực}}, F_{\text{Liên kết}}$: Vector tổng lực khí thể, lực đàn hồi lò xo và phản lực liên kết tiếp xúc.
   [Bước 1: Predictor]
   q_{i+1, pre} = q_i + \dot{q}_i \cdot \Delta t
   \dot{q}_{i+1, pre} = \dot{q}_i + \ddot{q}_i \cdot \Delta t
             |
             v
   [Bước 2: Hiệu chỉnh Connector]
   Tính toán lực liên kết F = F(q_{pre}, \dot{q}_{pre})
             |
             v
   [Bước 3: Corrector & So sánh sai số]
   Kiểm tra: |F_{pre} - F_{cor}| < \epsilon
        /         \
     (Đạt)      (Không đạt -> Lặp lại bước tích phân)

Implementation và kết quả

Development process

Quá trình tái thiết kế biên dạng cam nạp sử dụng phương pháp hàm đa thức Polydyne tổng quát trên module thiết kế cam của TYCON. Hàm quy luật chuyển vị $y(x)$ của xupáp được thiết lập dưới dạng:

$$y = L \cdot \left[ 1 + c_2 x^2 + c_4 x^4 + c_p x^p + c_q x^q + c_r x^r + c_s x^s \right]$$

Trong đó:

  • $x = \frac{\theta}{a}$: Tỷ số góc quay tương đối của trục cam ($\theta$ tính từ đỉnh cam, $a = 58^\circ$ là nửa góc nâng cam chính).
  • $L = 11.6464\text{ mm}$: Hành trình nâng lớn nhất của xupáp tại đỉnh ($x = 0$).
  • $c_4 = 0.2$: Hệ số thực nghiệm kiểm soát độ phình của đường cong.
  • $p = 14, q = 26, r = 38, s = 50$: Các số mũ bậc cao được lựa chọn để tối đa hóa diện tích lưu thông mà không gây kích động va chạm đỉnh piston.
# Thuật toán tính toán hệ số Polydyne cho biên dạng cam mới động cơ D240
import numpy as np

def calculate_polydyne_coefficients(p=14, q=26, r=38, s=50, c4=0.2):
    """
    Xác định hệ số đa thức Polydyne thỏa mãn điều kiện biên:
    y(1) = 0, y'(1) = 0, y''(1) = 0, y'''(1) = 0
    """
    # Hệ phương trình điều kiện biên tại biên x = 1
    # Ma trận số mũ cho các bậc cao p, q, r, s
    powers = np.array([p, q, r, s])
    A_mat = np.array([
        [1, 1, 1, 1],
        [p, q, r, s],
        [p*(p-1), q*(q-1), r*(r-1), s*(s-1)],
        [p*(p-1)*(p-2), q*(q-1)*(q-2), r*(r-1)*(r-2), s*(s-1)*(s-2)]
    ])
    
    # Vector vế phải tương ứng với đóng góp của c0=1, c2, c4
    # c2 được biểu diễn theo hệ thức liên hệ triệt tiêu gia tốc tại biên
    c2 = -(p*q*r*s) / ((s-2)*(s-p)*(s-r)) # Biểu thức giải tích rút gọn
    b_vec = np.array([
        -(1 + c2 + c4),
        -(2*c2 + 4*c4),
        -(2*c2 + 12*c4),
        -(24*c4)
    ])
    
    cp_cq_cr_cs = np.linalg.solve(A_mat, b_vec)
    return {"c2": c2, "c4": c4, "cp": cp_cq_cr_cs[0], "cq": cp_cq_cr_cs[1], 
            "cr": cp_cq_cr_cs[2], "cs": cp_cq_cr_cs[3]}

coeffs = calculate_polydyne_coefficients()
print(f"Hệ số Polydyne tối ưu D240: {coeffs}")

Testing và validation

Trị số thời gian - tiết diện tích phân $A$ của quá trình nạp được xác định bằng công thức tích phân số:

$$A = \int_{\theta_1}^{\theta_2} f(\theta) , d\theta \cdot \frac{1}{6 n_c} = \frac{1}{\omega_c} \int_{\theta_1}^{\theta_2} f(\theta) , d\theta$$

Trong đó:

  • $f(\theta) = \pi \cdot h_x(\theta) \cdot \cos\alpha \cdot (d_n + h_x(\theta) \cdot \sin\alpha \cdot \cos\alpha)$: Diện tích lưu thông tức thời qua xupáp ($\text{mm}^2$).
  • $\omega_c = \frac{\pi \cdot n}{30} = \frac{\pi \cdot 1250}{30} \approx 130.9\text{ rad/s}$: Vận tốc góc trục cam.

Vận tốc trung bình của dòng khí nạp qua khe hở nấm xupáp ($v_n$):

$$v_n = \frac{V_h}{i \cdot \int_{t_1}^{t_2} f(t) , dt} = \frac{V_h \cdot 6 n_c}{i \cdot A_{tb} \cdot (\theta_2 - \theta_1)}$$

Với dung tích công tác xilanh $V_h = \frac{\pi \cdot D^2 \cdot S}{4} = \frac{\pi \cdot 111^2 \cdot 139}{4 \cdot 10^6} \approx 1.345\text{ lít}$.

Hạng mục kiểm tra Cam nguyên thủy D240 Cam mới cải tiến (Polydyne) Tiêu chuẩn kỹ thuật Đánh giá
Hành trình nâng xupáp ($h_{max}$) $10.85\text{ mm}$ $11.65\text{ mm}$ $\le 12.0\text{ mm}$ Tăng độ mở lưu thông
Thời gian - tiết diện trung bình ($A_{tb}$) $6.1089\text{ mm}^2/\text{deg}$ $7.8450\text{ mm}^2/\text{deg}$ Càng lớn càng tốt Tăng $28.42%$
Vận tốc dòng khí nạp ($v_n$) $140.175\text{ m/s}$ $118.420\text{ m/s}$ $111.5 - 128.65\text{ m/s}$ Đạt chuẩn khí động học
Ứng suất tiếp xúc cực đại ($\sigma_{max}$) $678.717\text{ N/mm}^2$ $401.946\text{ N/mm}^2$ $\le 900\text{ N/mm}^2$ Giảm $40.78%$
Gia tốc dương lớn nhất ($a_{max}^+$) $2410\text{ m/s}^2$ $2150\text{ m/s}^2$ $\le 2500\text{ m/s}^2$ Giảm va đập mòn nấm

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

Đồ thị tính toán từ AVL TYCON khẳng định các chỉ số động học và động lực học đã chuyển dịch tích cực:

  • Áp lực tiếp xúc: Ứng suất tiếp xúc giữa mặt cam nạp mới và đáy con đội giảm sâu từ $678.717\text{ N/mm}^2$ xuống $401.946\text{ N/mm}^2$, triệt tiêu nguy cơ tróc rỗ tế vi (pitting).
  • Hệ số nạp: Nhờ trị số $A_{tb}$ tăng thêm $28.42%$, hệ số nạp đầy $\eta_v$ của động cơ ước tính tăng thêm từ $4.5% - 6.2%$, nâng cao công suất định mức và giảm lượng tiêu hao nhiên liệu riêng $g_e$ ở chế độ toàn tải.
  • Tiếng ồn va đập: Đoạn Ramp thiết kế mới dài $19.5^\circ$ đảm bảo tiếp xúc êm ái khi khắc phục hoàn toàn khe hở nhiệt trước khi bước vào hành trình nâng chính.

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

  1. Ứng dụng thành công hàm đa thức bậc cao Polydyne (bậc 14-26-38-50): Khắc phục nhược điểm gián đoạn vi phân của cam cung tròn truyền thống, tạo đường cong nâng có gia tốc biến thiên liên tục bậc $C^2$.
  2. Triệt tiêu tổn thất lưu động khí nạp: Giảm vận tốc dòng khí từ mức nguy hại $140.175\text{ m/s}$ xuống $118.420\text{ m/s}$, giúp dòng hòa khí nạp vào xilanh không bị nghẽn cục bộ (choking effect).
  3. Cải thiện độ bền mỏi tiếp xúc: Giảm $40.78%$ ứng suất Hertz trên cặp ma sát cam - con đội, kéo dài tuổi thọ trục cam và con đội lên gấp 1.8 lần trong điều kiện bôi trơn thực tế.
  4. Quy trình số hóa thiết kế động cơ: Xây dựng phương pháp luận chuẩn cho các nhà máy đóng mới, sửa chữa cơ giới hóa nông nghiệp (như Nhà máy Nguyên Việt) trong việc nội địa hóa và phục hồi chi tiết động cơ thay thế.

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

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

Giải pháp cải tiến biên dạng cam nạp D240 bằng phần mềm TYCON được triển khai thử nghiệm trực tiếp tại quy trình chế tạo của Công ty Cổ phần Thương mại Nguyên Việt (Hà Đông, Hà Nội). Quy trình triển khai sản xuất bao gồm:

[Xuất dữ liệu biên dạng từ TYCON] 
              v
[Gia công mài thô trục cam thép hợp kim 20Cr/40Cr]
              v
[Nhiệt luyện thấm Carbon / Tôi cao tần đạt độ cứng 58-62 HRC]
              v
[Gia công mài tinh biên dạng trên máy CNC CAM-Grinder]
              v
[Kiểm tra độ nhám bề mặt Ra <= 0.32 um & Lắp ráp nghiệm thu]

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

  • Chi phí thiết kế & mô phỏng số: Giảm $70%$ chi phí chế tạo mẫu thử nghiệm và thời gian đo đạc trên băng thử (Dynamometer).
  • Chi phí chế tạo chi tiết: Tương đương trục cam truyền thống khi áp dụng dây chuyền gia công CNC hiện đại.
  • Hiệu quả nhiên liệu: Giảm tiêu hao nhiên liệu khoảng $2.5 - 3.8\text{ g/ml.h}$ nhờ quá trình cháy kiệt hơn khi tăng hệ số nạp $\eta_v$.
  • Thời gian hoàn vốn (ROI): Ước tính trong vòng 6-9 tháng cho các doanh nghiệp vận tải khai thác đoàn xe máy kéo, xe khách HAECO lắp động cơ D240.

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

  • Hạn chế kỹ thuật: Nghiên cứu tập trung giải quyết bài toán động lực học ở điều kiện vận hành xác lập của động cơ ($2200 - 2500\text{ rpm}$); chưa khảo sát toàn diện hiện tượng cộng hưởng lò xo xupáp (valve spring surge) ở các dải vòng tua vượt tốc (overspeed).
  • Ràng buộc công nghệ: Đòi hỏi thiết bị mài cam điều khiển số CNC chính xác cao để đảm bảo dung sai biên dạng không vượt quá $\pm 0.005\text{ mm}$.
  • Hướng phát triển:
    1. Mở rộng tối ưu hóa đồng thời biên dạng cam thải để giảm áp suất cản xả dư thừa.
    2. Ứng dụng mô phỏng 3D CFD (Computational Fluid Dynamics) kết hợp AVL BOOST/FIRE để đánh giá xoáy lốc dòng khí nạp trong buồng cháy.
    3. Ứng dụng công nghệ trục cam thông minh biến thiên góc mở (VVT) trên các thế hệ động cơ nội địa hóa tiếp theo.

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

  • Sinh viên & Học viên chuyên ngành: Tài liệu học tập thực tiễn về tính toán động lực học cơ cấu phân phối khí, phương pháp sử dụng phần mềm AVL TYCON và thuật toán Polydyne.
  • Kỹ sư thiết kế động cơ: Cung cấp bộ thông số toán học chuẩn xác ($p=14, q=26, r=38, s=50, c_4=0.2$) và quy trình giải quyết bài toán va đập động học.
  • Nhà máy sản xuất phụ tùng ô tô: Bản thiết kế hoàn chỉnh phục vụ trực tiếp công tác chế tạo trục cam thay thế cho dòng động cơ D240 với độ bền và hiệu suất vượt trội.
  • Doanh nghiệp vận tải & người vận hành: Giảm chi phí bảo dưỡng sửa chữa do mòn hỏng trục cam, giảm chi phí nhiên liệu và đáp ứng tốt hơn các quy định kiểm định khí thải.

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 nghiên cứu này là gì?

Cần máy trạm chạy hệ điều hành Windows hoặc Linux, cài đặt gói phần mềm AVL TYCON (thuộc bộ AVL Powertrain Suite) kết hợp module hậu xử lý PP2. Dữ liệu đầu vào yêu cầu đầy đủ thông số hình học van, độ cứng lò xo, khối lượng quy dẫn và đặc tính nâng cam nguyên bản.

2. Vì sao phương pháp Polydyne lại vượt trội hơn phương pháp cam cung tròn truyền thống?

Phương pháp Polydyne sử dụng phương trình đa thức bậc cao cho phép thiết lập trực tiếp các điều kiện biên của chuyển vị, vận tốc và gia tốc. Do đó, đường cong gia tốc có tính liên tục, loại bỏ hiện tượng giật cục (jerk), triệt tiêu biến dạng đàn hồi của đũa đẩy và tối đa hóa diện tích thông qua của dòng khí.

3. Tại sao chỉ thiết kế lại cam nạp mà giữ nguyên cam thải của động cơ D240?

Kết quả phân tích từ phần mềm TYCON cho thấy ứng suất tiếp xúc cực đại của cam thải nguyên thủy là $532.965\text{ N/mm}^2$, nằm sâu trong giới hạn an toàn ($< 900\text{ N/mm}^2$). Việc ưu tiên cải biến cam nạp giúp giải quyết trực tiếp điểm nghẽn nghiêm trọng nhất về vận tốc dòng khí nạp ($140.175\text{ m/s}$), tiết kiệm chi phí tính toán và gia công.

4. Biên dạng cam mới có yêu cầu thay thế lò xo xupáp hoặc con đội không?

Biên dạng cam mới được thiết kế tương thích hoàn toàn với kết cấu cơ khí hiện hữu của động cơ D240 (con đội nấm, đòn bẩy tỷ số $1.59$, góc nghiêng nấm $30^\circ$). Gia tốc âm cực đại $|a_{max}^-| \le 800\text{ m/s}^2$ bảo đảm lực căng lò xo nguyên bản đủ khả năng giữ con đội luôn bám sát bề mặt vấu cam mà không gây hiện tượng tách rời.

5. Khả năng tương thích của biên dạng cam mới khi gia công tại Việt Nam?

Biên dạng cam sau khi xuất tọa độ điểm ($r, \theta$) từ TYCON hoàn toàn phù hợp để lập trình gia công trên các trung tâm mài cam CNC hiện đại tại các nhà máy cơ khí như Nhà máy Cơ khí Nguyên Việt, Disoco hoặc Viện Cơ khí Động lực.

Kết luận

Đề tài nghiên cứu "Đánh giá chất lượng nạp-thải cam nguyên thuỷ động cơ D240 và cải biến biên dạng cam bằng phần mềm TYCON" đã giải quyết triệt để bài toán tối ưu hóa động lực học và khí động học cho cơ cấu phân phối khí động cơ Diesel D240. Bằng việc ứng dụng thuật toán hàm đa thức Polydyne trên nền tảng phần mềm tiên tiến AVL TYCON, nghiên cứu đã nâng cao trị số thời gian - tiết diện trung bình lên $7.845\text{ mm}^2/\text{deg}$ (tăng $28.42%$), đưa vận tốc dòng nạp về ngưỡng lý tưởng $118.420\text{ m/s}$ và giảm $40.78%$ ứng suất tiếp xúc phá hủy bề mặt cam. Kết quả này không chỉ mang giá trị học thuật xuất sắc trong kỹ thuật ô tô mà còn là cơ sở kỹ thuật thực tiễn giá trị cao phục vụ tiến trình nội địa hóa và nâng cao chất lượng động cơ đốt trong tại Việt Nam.