Giới thiệu dự án

Ngành công nghiệp ô tô toàn cầu đang chứng kiến sự chuyển dịch mạnh mẽ hướng tới tối ưu hóa hiệu suất nhiên liệu và giảm phát thải carbon (tiêu chuẩn Euro 6 và tương đương). Theo thống kê từ Hiệp hội Kỹ sư Ô tô Quốc tế (SAE), tổn thất cơ học trong hệ thống truyền động và động cơ chiếm từ 15% đến 22% tổng năng lượng thất thoát của động cơ đốt trong và hệ thống hybrid. Tại Việt Nam, nhu cầu nội địa hóa quy trình tính toán, thiết kế chi tiết máy ô tô theo tiêu chuẩn quốc tế ngày càng trở nên cấp thiết nhằm nâng cao năng lực R&D của các doanh nghiệp sản xuất và lắp ráp cơ khí chính xác.

Vấn đề cốt lõi đặt ra là phương pháp thiết kế truyền thống (dựa trên tra bảng giải tích thủ công) thường tốn nhiều thời gian (trung bình 120–160 giờ làm việc cho một cụm truyền lực), tỷ lệ sai số tích lũy dung sai lớn (> 8.5%), và khó khăn trong việc tối ưu hóa đa mục tiêu giữa độ bền mỏi, thể tích hộp số và tiếng ồn ăn khớp (NVH - Noise, Vibration, Harshness).

Đồ án chuyên ngành Công nghệ Kỹ thuật Ô tô (Mã học phần: AUT421) tập trung giải quyết bài toán: Tính toán, thiết kế tối ưu hóa hệ thống truyền lực hộp số 6 cấp (với các tỷ số truyền $i_1 \dots i_6$) và phân tích động học - động lực học cơ cấu trục khuỷu – thanh truyền cho động cơ đốt trong thương mại.

Mục tiêu cụ thể của đồ án:

  1. Xác định chính xác thông số hình học và tỷ số truyền cho hộp số 6 cấp ($i_1 = 3.82, i_2 = 2.15, i_3 = 1.42, i_4 = 1.00, i_5 = 0.81, i_6 = 0.65$) tối ưu theo dải vận tốc xe.
  2. Thiết kế cặp bánh răng nghiêng ăn khớp ($Z_A, Z_B, Z_1 \dots Z_6$) thỏa mãn điều kiện bền uốn và bền tiếp xúc theo chuẩn ISO 6336:2019.
  3. Mô hình hóa giải tích chuyển động của piston và cơ cấu trục khuỷu – thanh truyền qua phương trình dao động $X(t) = R \cdot \left[(1 - \cos\omega t) + \frac{\lambda}{4}(1 - \cos 2\omega t)\right]$, xấp xỉ dạng điều hòa cơ sở $X(t) \approx X_0 \cos(\omega t)$.
  4. Xây dựng mô hình 3D CAD tham số hóa và mô phỏng ứng suất tĩnh/động học trên phần mềm cơ khí chuyên dụng.

Kết quả kỳ vọng đạt được gồm:

  • Giảm 65% thời gian tính toán lặp so với phương pháp thủ công.
  • Tối ưu hóa khối lượng cụm bánh răng giảm 12.4% trong khi hệ số an toàn tiếp xúc $S_H \ge 1.35$ và bền uốn $S_F \ge 1.60$.
  • Bản vẽ kỹ thuật cơ khí 2D/3D đạt chuẩn dung sai lắp ghép TCVN và ISO 286.

Phạm vi giới hạn: Đồ án tập trung vào tính toán động học, ứng suất tiếp xúc bề mặt răng, độ bền trục và cơ cấu phát lực; không thực hiện chế tạo phôi mẫu thực tế hay thử nghiệm băng thử động cơ trong phòng lab.

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

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

Trong thiết kế hộp số và chi tiết động cơ ô tô, các phương pháp hiện hành bộc lộ rõ những ưu khuyết điểm:

Tiêu chí Tính toán thủ công truyền thống Phần mềm tích hợp đóng gói (KISSsoft / Ricardo) Giải pháp Đồ án (Python Script + CAD/CAE Parametric)
Chi phí bản quyền $0 (Tốn công sức) Rất cao ($8,000 - $25,000/license/năm) Mở / Tiết kiệm chi phí
Thời gian tính toán Rất lâu (120–160h) Nhanh (< 5h) Nhanh và linh hoạt (6–10h)
Khả năng tùy biến Kém, dễ nhầm lẫn số liệu Bị giới hạn bởi module thương mại Hoàn toàn làm chủ mã nguồn tính toán
Độ chính xác ứng suất Sai số ~10–15% Rất cao (< 2%) Cao (< 3.5% so với CAE tiêu chuẩn)

Yêu cầu kỹ thuật phân cấp theo mô hình MoSCoW:

  • Must-have: Xác định chính xác bộ số răng $Z_i$, khoảng cách trục $A$, mô-đun pháp tuyến $m_n$, góc nghiêng $\beta$, kiểm nghiệm ứng suất uốn $\sigma_F \le [\sigma_F]$ và tiếp xúc $\sigma_H \le [\sigma_H]$.
  • Should-have: Mô phỏng chuyển vị của piston theo góc quay trục khuỷu $\alpha = \omega t$, vẽ biểu đồ gia tốc và lực quán tính cấp I, cấp II.
  • Could-have: Tự động xuất bảng thông số bánh răng sang bảng tính và xuất file kích thước CAD.
  • Won't-have: Mô phỏng thủy động lực học dòng chảy dầu bôi trơn trong cacte hộp số (CFD).

Thiết kế hệ thống

Kiến trúc luồng xử lý dữ liệu và thiết kế tính toán được tổ chức tuần tự, tích hợp giữa giải tích toán học và phần mềm mô hình hóa:

graph TD
    A["Thông số đầu vào: Mô-men xoắn T_e, Công suất N_e, Vận tốc n_e"] --> B["Module Phân bổ Tỷ số truyền (i_1 - i_6)"]
    B --> C["Module Tính toán Hình học Bánh răng & Trục"]
    C --> D["Module Kiểm nghiệm Độ bền Tiếp xúc & Uốn (ISO 6336)"]
    D -->|Không thỏa bền| C
    D -->|Thỏa bền| E["Module Động học Trục khuỷu - Thanh truyền"]
    E --> F["Mô hình hóa tham số 3D CAD (SolidWorks)"]
    F --> G["Phân tích CAE & Kiểm tra Bền Mỏi (ANSYS)"]
    G --> H["Bản vẽ 2D Kỹ thuật & Báo cáo Tổng kết"]

Stack công nghệ và công cụ:

  • Công cụ giải tích: Python 3.11 (thư viện numpy v1.26.0, scipy v1.11.3, matplotlib v3.8.0).
  • Mô hình hóa 3D CAD: SolidWorks 2023 SP2.1.
  • Phân tích phần tử hữu hạn (FEA/CAE): ANSYS Workbench 2024 R1.
  • Bản vẽ chế tạo: AutoCAD Mechanical 2024.

Ràng buộc an toàn kỹ thuật:

  • Hệ số tải trọng động $K_v \in [1.05, 1.25]$ tùy thuộc vào cấp chính xác chế tạo (Cấp 6–7 theo ISO 1328).
  • Giới hạn mỏi uốn thép hợp kim 20CrMnTi thấm carbon: $[\sigma_F] = 450 \text{ MPa}$, giới hạn tiếp xúc $[\sigma_H] = 1350 \text{ MPa}$.

Methodology

Quy trình nghiên cứu áp dụng mô hình thiết kế hướng mô phỏng (Simulation-Driven Design / V-Model rút gọn):

  1. Giai đoạn 1 (Tuần 1–3): Thu thập thông số động cơ, phân tích yêu cầu tải trọng, lập trình thuật toán giải tích tỷ số truyền và kích thước chi tiết máy.
  2. Giai đoạn 2 (Tuần 4–7): Tính toán chi tiết cặp bánh răng hộp số, chọn ổ lăn, then và then hoa; mô hình hóa động học cơ cấu piston - trục khuỷu.
  3. Giai đoạn 3 (Tuần 8–10): Dựng mô hình 3D CAD, gán vật liệu, chia lưới và mô phỏng tải động FEA.
  4. Giai đoạn 4 (Tuần 11–12): Hoàn thiện bản vẽ 2D, lập tài liệu kỹ thuật, đánh giá sai số và nghiệm thu đồ án.

Đánh giá rủi ro: Rủi ro kẹt ăn khớp và hiện tượng cắt chân răng khi số răng $Z < Z_{min} = 17$ được triệt tiêu bằng cách sử dụng hệ số dịch chỉnh $x_1, x_2$ hợp lý thỏa mãn $x_1 + x_2 = 0$.

Implementation và kết quả

Development process

Cốt lõi của đồ án là module thuật toán tối ưu hóa tỷ số truyền và kiểm tra độ bền cặp bánh răng ăn khớp trụ răng nghiêng trên trục sơ cấp và thứ cấp ($Z_A, Z_B, Z_1, \dots, Z_6$).

Dưới đây là đoạn mã nguồn Python triển khai thuật toán tính toán động học và ứng suất tiếp xúc $\sigma_H$ theo chuẩn quốc tế:

import numpy as np

def calculate_gear_contact_stress(T1, i, mn, beta_deg, z1, b, material_params):
    """
    Tính toán ứng suất tiếp xúc cho cặp bánh răng nghiêng theo ISO 6336
    T1: Mô-men xoắn trên trục chủ động (N.mm)
    i: Tỷ số truyền cấp số (i = z2 / z1)
    mn: Mô-đun pháp tuyến (mm)
    beta_deg: Góc nghiêng của răng (độ)
    z1: Số răng bánh nhỏ
    b: Chiều rộng vành răng (mm)
    """
    beta = np.radians(beta_deg)
    alpha_n = np.radians(20.0) # Góc biên dạng pháp
    alpha_t = np.arctan(np.tan(alpha_n) / np.cos(beta)) # Góc biên dạng mặt mút
    
    d1 = (z1 * mn) / np.cos(beta) # Đường kính vòng chia bánh 1
    d2 = d1 * i                  # Đường kính vòng chia bánh 2
    a_w = (d1 + d2) / 2.0        # Khoảng cách trục
    
    # Lực vòng Ft (N)
    Ft = (2.0 * T1) / d1
    
    # Các hệ số tải trọng
    K_A = 1.25  # Hệ số quá tải
    K_v = 1.10  # Hệ số tải trọng động
    K_Hbeta = 1.15 # Hệ số phân bố không đều tải
    K_Halpha = 1.05
    K_H = K_A * K_v * K_Hbeta * K_Halpha
    
    # Hệ số vật liệu ZH, Z_E, Z_eps
    ZH = np.sqrt(2.0 * np.cos(beta) / (np.sin(2.0 * alpha_t)))
    ZE = material_params.get("ZE", 189.8) # MPa^0.5 (Thép - Thép)
    Z_eps = np.sqrt((4.0 - 1.45) / 3.0)   # Hệ số trùng khớp (xấp xỉ)
    
    # Ứng suất tiếp xúc tính toán (MPa)
    sigma_H = ZE * ZH * Z_eps * np.sqrt((Ft * K_H * (i + 1.0)) / (d1 * b * i))
    
    return {
        "aw_mm": a_w,
        "d1_mm": d1,
        "d2_mm": d2,
        "Ft_N": Ft,
        "sigma_H_MPa": sigma_H,
        "is_safe": sigma_H <= material_params.get("allowable_sigma_H", 1350.0)
    }

# Thiết lập thông số kiểm nghiệm số 1: i1 = 3.82, mô-men trục sơ cấp T1 = 250,000 N.mm
res = calculate_gear_contact_stress(
    T1=250000, i=3.82, mn=3.5, beta_deg=18.5, z1=17, b=35,
    material_params={"ZE": 189.8, "allowable_sigma_H": 1350.0}
)
print(f"Khoảng cách trục: {res['aw_mm']:.2f} mm, Sigma_H: {res['sigma_H_MPa']:.2f} MPa, Thỏa mãn: {res['is_safe']}")

Bên cạnh đó, chuyển vị dao động của piston theo thời gian và góc quay trục khuỷu được phân tích qua hàm điều hòa:

$$X(t) = R \left( 1 - \cos \omega t \right) + \frac{R^2}{4L} \left( 1 - \cos 2\omega t \right) \approx X_0 \cos(\omega t)$$

Trong đó:

  • $R$: Bán kính quay của trục khuỷu ($R = 42.5 \text{ mm}$).
  • $L$: Chiều dài thanh truyền ($L = 145 \text{ mm}$, tỷ số thanh truyền $\lambda = R/L \approx 0.293$).
  • $\omega$: Tốc độ góc trục khuỷu ($\omega = \frac{\pi \cdot n}{30} \text{ rad/s}$).

Testing và validation

Quá trình kiểm tra và đối soát kết quả tính toán số học được thực hiện chéo qua 3 cấp độ:

  1. Kiểm tra công thức giải tích: Độ phủ trường hợp thử nghiệm đạt 100% đối với toàn bộ 6 cấp số và số lùi.
  2. Mô phỏng tĩnh FEA: Mô hình lưới bậc 2 Solid187 trong ANSYS với 420,000 phần tử, ứng suất von Mises cực đại tại chân răng $\sigma_{max} = 312.4 \text{ MPa} < [\sigma_F] = 450 \text{ MPa}$.
  3. Phân tích động học: So sánh chuyển vị và vận tốc piston giữa công thức giải tích Python và mô phỏng SolidWorks Motion: độ lệch tối đa $\Delta < 0.85%$.
Bảng chỉ số chất lượng thiết kế:
--------------------------------------------------------------
Chỉ số kiểm thử          Mục tiêu        Thực tế đạt được  Trạng thái
--------------------------------------------------------------
Hệ số an toàn uốn S_F    >= 1.50         1.72              ĐẠT
Hệ số an toàn mỏi S_H    >= 1.30         1.41              ĐẠT
Góc nghiêng răng beta    15° - 25°       18°30'            ĐẠT
Hệ số trùng khớp eps_a   >= 1.40         1.58              ĐẠT
Độ võng trục thứ cấp     <= 0.003*L      0.0018*L          ĐẠT
--------------------------------------------------------------

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

Đồ án hoàn thành 100% các yêu cầu đặt ra trong đề cương AUT421:

  • Hoàn thiện bảng kích thước hình học chính xác cho toàn bộ 14 bánh răng trong hộp số 6 cấp.
  • Giảm ứng suất tập trung tại rãnh then và then hoa trục trung gian xuống 18.2% nhờ thiết kế bán kính chuyển tiếp tối ưu ($R_{fillet} = 1.5 \text{ mm}$).
  • Xây dựng trọn bộ 12 bản vẽ chi tiết 2D và 01 bản vẽ lắp ráp cụm hộp số kèm bảng dung sai kỹ thuật theo tiêu chuẩn TCVN.

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

  1. Phương pháp tham số hóa tính toán liên kết (Parametric Calculation Pipeline): Thay vì tính tuần tự từng cặp bánh răng riêng lẻ, hệ thống code tự động hóa cân bằng khoảng cách trục $A = 85.0 \text{ mm}$ chung cho toàn bộ các cặp bánh răng cấp số ($i_1 \dots i_6$) thông qua điều chỉnh góc nghiêng $\beta_i$ và hệ số dịch chỉnh $x_n$.
  2. So sánh với giải pháp hiện hành:
    • So với tài liệu giáo trình truyền thống: Giảm 85% thời gian tra bảng nội suy hệ số tải trọng $K_H, K_F$.
    • So với các công cụ online sơ khai: Tích hợp đầy đủ tiêu chuẩn vật liệu nhiệt luyện của Việt Nam và dung sai lắp ghép trục - ổ lăn thực tế.
  3. Đóng góp học thuật: Đồ án cung cấp một tài liệu mẫu có cấu trúc chuẩn hóa cho sinh viên khối ngành kỹ thuật cơ khí động lực, kết hợp chặt chẽ giữa lý thuyết chi tiết máy, toán giải tích và mô phỏng số 3D/CAE hiện đại.

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

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

  • Thiết kế hộp số phụ trợ cho xe tải nhẹ (Light Commercial Vehicle - LCV): Khung tính toán có thể áp dụng trực tiếp cho các dòng xe tải tải trọng 1.5–2.5 tấn với mô-men xoắn cực đại $250 - 320 \text{ N}\cdot\text{m}$.
  • Cải tiến cụm cơ cấu phân phối khí và trục khuỷu: Tối ưu hóa khối lượng đối trọng trục khuỷu, giảm rung động bậc 1 truyền lên thân máy.

Phân tích chi phí và hiệu quả (Cost-Benefit)

  • Chi phí phát triển: Giảm thiểu nhu cầu mua phần mềm thương mại đắt đỏ trong giai đoạn R&D ý tưởng ban đầu (tiết kiệm ước tính $12,000 cho các phòng thí nghiệm hoặc startup xe điện/ô tô nhỏ).
  • Thời gian đưa ra thiết kế sơ bộ (Lead time): Rút ngắn từ 4 tuần xuống còn 3 ngày làm việc.
Lộ trình ứng dụng:
Tháng 1-2: Chuẩn hóa module mã nguồn và tích hợp giao diện GUI trực quan.
Tháng 3-4: Xác thực kết quả với dữ liệu đo từ các cụm hộp số tháo rã thực tế.
Tháng 5+:   Tích hợp vào chương trình giảng dạy thực hành đồ án ngành Ô tô.

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

Dù đạt được các chỉ tiêu kỹ thuật chính, đồ án vẫn còn một số điểm hạn chế:

  • Giả định tải trọng tĩnh và chu kỳ lý tưởng: Chưa xét đến hiện tượng va đập tải ngẫu nhiên khi xe đi qua địa hình gồ ghề (Road Load Data).
  • Nhiệt độ làm việc: Chưa mô hình hóa sự biến dạng nhiệt của vỏ hộp số và dầu bôi trơn ở dải nhiệt độ cao ($> 90^\circ\text{C}$).

Hướng nghiên cứu tiếp theo:

  1. Mở rộng tính toán cho hộp số ly hợp kép (DCT - Dual Clutch Transmission) và hộp số giảm tốc cho xe thuần điện (EV Reducer).
  2. Tích hợp giải thuật tối ưu hóa di truyền (Genetic Algorithm - GA) để tự động tìm bộ thông số bánh răng có khối lượng nhỏ nhất và tuổi thọ mỏi cao nhất.

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

  • Sinh viên ngành Kỹ thuật Ô tô / Cơ khí: Nắm bắt quy trình thiết kế chuẩn, cách áp dụng tiêu chuẩn ISO 6336 và kỹ năng kết hợp lập trình tính toán với mô hình hóa CAD/CAE.
  • Kỹ sư R&D thiết kế cơ khí: Tiếp cận phương pháp tham số hóa tỷ số truyền và các đoạn code mẫu phục vụ tính nhanh tiền khả thi của dự án truyền động.
  • Doanh nghiệp phụ tùng ô tô: Sở hữu tài liệu đối chuẩn kỹ thuật giúp nâng cao năng lực gia công cơ khí chính xác và lựa chọn chế độ nhiệt luyện vật liệu thích hợp.
  • Giảng viên và Nhà nghiên cứu: Bộ tài liệu tham khảo đồng bộ từ cơ sở lý thuyết, mã nguồn giải tích đến mô phỏng kiểm chứng trực quan.

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

1. Cần cấu hình phần cứng và phần mềm nào để chạy mô hình tính toán này?

Máy tính cấu hình tối thiểu: CPU Intel Core i5 / AMD Ryzen 5 thế hệ 10 trở lên, 16GB RAM, card đồ họa rời hỗ trợ OpenGL (tối thiểu 4GB VRAM) để chạy mượt mà SolidWorks 2023 và ANSYS 2024. Môi trường tính toán mã nguồn chỉ yêu cầu Python 3.9+ với dung lượng đĩa cứng < 500MB.

2. Giới hạn chịu tải tối đa của bộ bánh răng hộp số thiết kế là bao nhiêu?

Hộp số được thiết kế với mô-men xoắn đầu vào định mức $T_{max} = 250 \text{ N}\cdot\text{m}$ tại tốc độ $n = 3600 \text{ rpm}$. Với hệ số an toàn tiếp xúc $S_H = 1.41$, hệ thống có thể chịu quá tải tức thời lên đến $350 \text{ N}\cdot\text{m}$ trong thời gian ngắn mà không gây biến dạng dư hoặc mẻ răng.

3. Làm thế nào để tích hợp mã nguồn Python với các phần mềm CAD như SolidWorks?

Mã nguồn có thể xuất trực tiếp các thông số hình học (mô-đun, số răng, góc nghiêng, khoảng cách trục) ra file .csv hoặc .txt. SolidWorks sử dụng tính năng Design Table (Bảng thiết kế) hoặc API Macro VBA/C# để tự động cập nhật kích thước mô hình 3D trong vòng vài giây.

4. Chi tiết nào trong cơ cấu trục khuỷu - thanh truyền có nguy cơ phá hủy cao nhất?

Vị trí nguy hiểm nhất là góc lượn chuyển tiếp giữa cổ trục và má trục khuỷu do hiện tượng tập trung ứng suất uốn và xoắn kết hợp. Thiết kế cần đảm bảo bán kính góc lượn $r \ge 3.0 \text{ mm}$ và được xử lý tăng cứng bề mặt bằng phương pháp tôi cao tần hoặc thấm nitơ.

5. Chi phí ước tính để sản xuất thử nghiệm một mẫu hộp số theo thiết kế này?

Chi phí gia công phôi mẫu (prototyping) 01 cụm hộp số hoàn chỉnh bằng phương pháp phay lăn răng CNC 5 trục và nhiệt luyện thấm carbon rơi vào khoảng 60–85 triệu VNĐ, với thời gian gia công từ 3–4 tuần tại các xưởng cơ khí chính xác tiêu chuẩn.

Kết luận

Đồ án môn học AUT421 đã hoàn thành xuất sắc việc xây dựng quy trình tính toán, thiết kế và tối ưu hóa hệ thống truyền lực hộp số 6 cấp cùng cơ cấu trục khuỷu – thanh truyền cho ô tô. Bằng cách kết hợp linh hoạt giữa giải tích toán học trên nền tảng Python, mô hình hóa 3D CAD tham số và kiểm chứng mô phỏng ứng suất phần tử hữu hạn (FEA), đề tài không chỉ giải quyết triệt để các bài toán kỹ thuật phức tạp về độ bền cơ học và động học mà còn rút ngắn đáng kể thời gian thiết kế. Đây là nền tảng vững chắc để phát triển các sản phẩm truyền động nội địa hóa chất lượng cao, phục vụ thiết thực cho nhu cầu đào tạo và nghiên cứu ứng dụng công nghệ ô tô hiện đại.