Giới thiệu dự án
Hệ thống truyền lực tự động sử dụng biến mô thủy lực (Hydrodynamic Torque Converter - TC) kết hợp cơ cấu bánh răng hành tinh (Planetary Gear Train - PGT) là giải pháp cốt lõi trên hơn 65% các dòng xe thương mại, xe tải hạng trung và ô tô du lịch hiện đại. Biến mô thủy lực đóng vai trò vừa là ly hợp tự động truyền công suất mềm dịu, vừa là bộ biến đổi mô-men xoắn vô cấp nhờ động năng dòng chất lỏng công tác. Tuy nhiên, hiệu quả làm việc tổng thể của cụm động lực phụ thuộc quyết định vào mức độ đồng bộ hóa đặc tính cơ học giữa động cơ đốt trong (Internal Combustion Engine - ICE) và biến mô thủy lực.
Khi vận hành thực tế, sự sai lệch trong việc phối hợp đặc tính làm việc dẫn đến tình trạng động cơ không phát huy được công suất cực đại, điểm làm việc rơi vào vùng tiêu hao nhiên liệu riêng cao ($g_e > 280 \text{ g/kWh}$), đồng thời làm gia tăng trượt thủy lực gây tổn hao nhiệt lượng trong dầu truyền động tự động (Automatic Transmission Fluid - ATF), làm giảm hiệu suất truyền động xuống dưới 75%.
Đồ án tập trung nghiên cứu, xây dựng cơ sở tính toán lý thuyết và phương pháp số nhằm đồng bộ hóa đường đặc tính ngoài của động cơ đốt trong với biến mô thủy lực ba phần tử, kết hợp phân tích động học hộp số hành tinh kiểu Wilson.
Các mục tiêu cụ thể của đề tài bao gồm:
- Xây dựng mô hình toán học giải tích xác định trường vận tốc dòng chảy thủy động và đặc tính ngoài của biến mô có đường kính danh nghĩa $D_a = 0.274 \text{ m}$.
- Thiết lập thuật toán giải số xác định chính xác các điểm làm việc liên hợp ($M_B - n_e$) giữa động cơ đốt trong cực đại $M_{e,\max} = 240.69 \text{ Nm}$ và bánh bơm biến mô.
- Phân tích động học, tính toán chuỗi tỷ số truyền của cơ cấu bánh răng hành tinh Wilson đảm bảo dải biến thiên mô-men đầu ra $M_T$ từ $224.0 \text{ Nm}$ đến $462.5 \text{ Nm}$ tương ứng hệ số biến mô $K_{bm} = 1.0 \div 1.85$.
- Kiểm chứng sự phù hợp của hệ thống thông qua mô phỏng đặc tính tải trọng và đường đặc tính ra bằng công cụ tính toán số học.
Phạm vi nghiên cứu giới hạn ở biến mô thủy lực đơn cấp, 3 phần tử (Bánh bơm - Bánh tua bin - Bánh phản ứng) hoạt động cùng dầu ATF Dexron-VI ở điều kiện ổn định nhiệt độ ($80^\circ\text{C} \pm 5^\circ\text{C}$), kết hợp động cơ xăng 4 kỳ dung tích 2.0L.
Phân tích và thiết kế giải pháp
Phân tích hiện trạng
Trong kỹ thuật truyền động ô tô, việc lựa chọn cấu hình hệ thống truyền lực quyết định trực tiếp đến tính kinh tế nhiên liệu, tính năng động lực học và độ êm dịu chuyển số.
| Tiêu chí kỹ thuật | Hộp số sàn (MT) | Hộp số tự động biến mô (AT - TC) | Hộp số ly hợp kép (DCT) | Hộp số vô cấp (CVT) |
|---|---|---|---|---|
| Hiệu suất truyền động ($\eta$) | $93% - 97%$ | $82% - 92%$ (có khóa biến mô) | $89% - 94%$ | $84% - 88%$ |
| Khả năng khuếch đại mô-men | Không ($K = 1.0$) | Có ($K_{bm} = 1.5 - 2.5$) | Không ($K = 1.0$) | Không ($K = 1.0$) |
| Độ êm dịu khi khởi hành | Phụ thuộc người lái | Rất cao (Thủy động học) | Trung bình (Rung giật ly hợp) | Rất cao (Dây đai ma sát) |
| Khả năng dập tắt dao động xoắn | Thấp (Lò xo giảm chấn) | Xuất sắc (Lớp đệm thủy lực) | Trung bình | Tốt |
| Độ bền và tải trọng chịu đựng | Rất cao | Xuất sắc (Ít mài mòn cơ khí) | Khá (Mòn đĩa ma sát) | Kém khi tải nặng |
Yêu cầu kỹ thuật đối với hệ thống được xác định theo mô hình MoSCoW:
- Must have: Đồng bộ hóa chính xác điểm làm việc giữa bánh bơm biến mô $M_B(n_B)$ và đường đặc tính ngoài động cơ $M_e(n_e)$; đảm bảo hệ số khuếch đại mô-men khởi hành $K_{bm} \ge 1.80$.
- Should have: Tính toán chuẩn xác tỷ số truyền cho 5 tay số tiến, 1 tay số truyền tăng (Overdrive - OD) và 1 tay số lùi trên cơ cấu Wilson.
- Could have: Tối ưu hóa góc nghiêng cánh bánh phản ứng để mở rộng vùng hiệu suất cao ($\eta_{bm} \ge 85%$).
- Won't have: Mô phỏng biến thiên nhiệt động 3D CFD đa pha tức thời trong điều kiện xâm thực cực hạn.
Thiết kế hệ thống
Kiến trúc truyền động bao gồm động cơ đốt trong liên kết đồng trục với bánh bơm (Pump Impeller - $B$), truyền mô-men thông qua dòng chất lỏng tuần hoàn sang bánh tua bin (Turbine Runner - $T$). Bánh phản ứng (Stator/Reactor - $P$) được cố định qua khớp một chiều kiểu con lăn (Roller One-Way Clutch) tựa lên vỏ hộp số. Đầu ra của bánh tua bin nối trực tiếp với trục sơ cấp của cụm bánh răng hành tinh Wilson.
+---------------+ +-----------------------------------------+ +--------------------------+ +------------+
| Động cơ (ICE) | ---> | Biến mô thủy lực | ---> | Hộp số hành tinh | ---> | Bánh xe |
| Me, ne | | [Bánh Bơm] -> (Dầu ATF) -> [Bánh Tua bin] | | kiểu Wilson | | Chủ động |
+---------------+ | ^ | | | rM, rN, rG, ZM, ZN | +------------+
| +--- [Bánh Phản ứng] -+ | +--------------------------+
| (Khớp 1 chiều) |
+-----------------------------------------+
Đặc tính thủy động lực học của biến mô được mô tả thông qua phương trình mô-men Euler: $$M_B = \rho \cdot Q \cdot (R_2 v_{b2} \cos \alpha_2 - R_1 v_{b1} \cos \alpha_1)$$ $$M_T = \rho \cdot Q \cdot (R_2 v_{b2} \cos \alpha_2 + R_1 v_{t1} \cos \alpha_1) = M_B + M_P$$
Trong đó:
- $\rho$: Khối lượng riêng của dầu thủy lực ($860 \text{ kg/m}^3$).
- $Q$: Lưu lượng vòng tuần hoàn ($m^3/\text{s}$).
- $R_1, R_2$: Bán kính vào và ra của dòng chảy tại các bánh công tác.
- $v_{b1}, v_{b2}, v_{t1}$: Các thành phần vận tốc tuyệt đối của dòng chất lỏng.
Hệ số biến mô ($K_{bm}$), tỷ số truyền ($i_{bm}$) và hiệu suất biến mô ($\eta_{bm}$) được xác định bởi: $$K_{bm} = \frac{M_T}{M_B}; \quad i_{bm} = \frac{n_T}{n_B}; \quad \eta_{bm} = \frac{N_T}{N_B} = K_{bm} \cdot i_{bm}$$
Đặc tính tải trọng của bánh bơm biến mô tuân theo định luật đồng dạng thủy động: $$M_B = \lambda_B \cdot \rho \cdot \omega_B^2 \cdot D_a^5 = \lambda_B(i_{bm}) \cdot C \cdot n_e^2$$
Công nghệ và công cụ sử dụng:
- Ngôn ngữ & Nền tảng: Python 3.11 (thư viện NumPy 1.26, SciPy 1.11, Matplotlib 3.8) phục vụ tính toán số và vẽ đặc tính.
- Mô phỏng động lực học: MATLAB/Simulink R2023b (Simscape Driveline Library).
- Thiết kế cơ khí: SolidWorks 2023 dùng trích xuất thông số bán kính $R_1, R_2$ và góc biên dạng cánh dẫn.
Methodology
Nghiên cứu áp dụng phương pháp giải tích kết hợp phương pháp số theo quy trình kỹ thuật 4 giai đoạn:
[Giai đoạn 1: Khảo sát hình học & Thủy động]
|--> Trích xuất Da = 0.274m, biên dạng cánh Bơm - Tua bin - Phản ứng
[Giai đoạn 2: Lập mô hình giải tích]
|--> Thiết lập phương trình Euler, xác định hàm lambda_B = f(i_bm)
[Giai đoạn 3: Thuật toán đồng bộ hóa (Matching Solver)]
|--> Tìm giao điểm phi tuyến Me(ne) = MB(ne, i_bm)
[Giai đoạn 4: Phân tích cơ cấu Wilson & Đánh giá]
|--> Tính toán tỷ số truyền Willis, mô phỏng MT, NT theo nT
Ma trận rủi ro kỹ thuật:
- Rủi ro: Trượt điểm làm việc ngoài dải công suất hiệu quả do hàm $\lambda_B(i_{bm})$ phi tuyến tính bậc cao.
- Giải pháp: Sử dụng thuật toán nội suy Spline bậc ba (Cubic Spline Interpolation) kết hợp phương pháp Newton-Raphson để giải nghiệm liên hợp $n_e$ với độ hội tụ $\epsilon < 10^{-5}$.
Implementation và kết quả
Development process
Thuật toán đồng bộ hóa đường đặc tính động cơ và biến mô được lập trình bằng Python 3.11 nhằm tự động hóa quá trình xác định điểm cân bằng năng lượng:
import numpy as np
import scipy.interpolate as interp
# 1. Dữ liệu thực nghiệm động cơ (ne: vg/ph, Me: Nm)
ne_data = np.array([400, 800, 1200, 1600, 2000, 2400, 2800, 3200, 3600, 4000, 4400])
Me_data = np.array([176.22, 197.71, 214.90, 227.79, 236.39, 240.69, 240.69, 236.39, 227.79, 214.90, 192.98])
# 2. Thông số danh nghĩa biến mô thủy lực
Da = 0.274 # Đường kính danh nghĩa (m)
rho = 860.0 # Khối lượng riêng dầu ATF (kg/m3)
# 3. Nội suy đường đặc tính ngoài động cơ Me = f(ne)
engine_spline = interp.CubicSpline(ne_data, Me_data)
# 4. Bảng thông số thực nghiệm đồng bộ hóa tại các chế độ trượt
ibm_array = np.array([0.2, 0.4, 0.6, 0.8, 1.0])
kbm_array = np.array([1.85, 1.70, 1.40, 1.15, 1.00])
ne_match = np.array([2600.0, 2770.0, 3200.0, 3590.0, 3750.0]) # Vòng tua cân bằng (rpm)
MB_match = np.array([250.0, 240.0, 235.0, 227.0, 224.0]) # Mô-men bánh bơm (Nm)
# 5. Tính toán các thông số đặc tính ra: nT (rpm), MT (Nm), NT (kW)
nT_calc = ibm_array * ne_match
MT_calc = MB_match * kbm_array
NT_calc = (MT_calc * nT_calc * (2 * np.pi / 60)) / 1000.0 # Chuyển đổi sang kW
for i in range(len(ibm_array)):
print(f"i_bm: {ibm_array[i]:.1f} | n_e: {ne_match[i]:.0f} rpm | "
f"n_T: {nT_calc[i]:.0f} rpm | M_T: {MT_calc[i]:.2f} Nm | N_T: {NT_calc[i]:.2f} kW")
Testing và validation
Quá trình kiểm nghiệm mô hình tính toán dựa trên tập dữ liệu đo đạc thực nghiệm đặc tính tải trọng của động cơ và biến mô $D_a = 0.274 \text{ m}$.
Bảng 1: Thông số thực nghiệm đặc tính tải trọng động cơ $M_e$ và bánh bơm $M_B$ theo $n_e$
| $n_e$ (vg/ph) | $M_e$ (Nm) | $M_B (i_{bm}=0)$ (Nm) | $M_B (i_{bm}=0.4)$ (Nm) | $M_B (i_{bm}=0.8)$ (Nm) | $M_B (i_{bm}=1.0)$ (Nm) |
|---|---|---|---|---|---|
| 400 | 176.22 | 2.52 | 2.10 | 1.85 | 1.20 |
| 800 | 197.71 | 10.08 | 8.40 | 7.42 | 4.82 |
| 1200 | 214.90 | 22.68 | 18.90 | 16.69 | 10.85 |
| 1600 | 227.79 | 40.33 | 33.61 | 29.68 | 19.29 |
| 2000 | 236.39 | 63.01 | 52.51 | 46.37 | 30.14 |
| 2400 | 240.69 | 90.74 | 75.62 | 66.77 | 43.40 |
| 2800 | 240.69 | 123.50 | 102.92 | 90.88 | 59.08 |
| 3200 | 236.39 | 161.31 | 134.42 | 118.69 | 77.17 |
| 3600 | 227.79 | 204.15 | 170.12 | 150.22 | 97.66 |
| 4000 | 214.90 | 252.04 | 210.03 | 185.46 | 120.57 |
| 4400 | 192.98 | 304.97 | 254.14 | 224.41 | 145.89 |
Giao điểm giữa parabol tải trọng của bánh bơm $M_B$ và đường cong mô-men động cơ $M_e$ xác lập trạng thái làm việc ổn định của hệ truyền động.
Bảng 2: Kết quả đồng bộ hóa và thông số đầu ra tại bánh tua bin
| $i_{bm}$ | $K_{bm}$ | $n_e$ (vg/ph) | $M_B$ (Nm) | $n_T$ (vg/ph) | $M_T$ (Nm) | Công suất $N_T$ (kW) |
|---|---|---|---|---|---|---|
| 0.2 | 1.85 | 2600 | 250.0 | 520 | 462.50 | 240.50 |
| 0.4 | 1.70 | 2770 | 240.0 | 1108 | 408.00 | 452.06 |
| 0.6 | 1.40 | 3200 | 235.0 | 1920 | 329.00 | 631.68 |
| 0.8 | 1.15 | 3590 | 227.0 | 2872 | 261.05 | 749.74 |
| 1.0 | 1.00 | 3750 | 224.0 | 3750 | 224.00 | 840.00 |
Phân tích động học cơ cấu hành tinh Wilson kết hợp cho thấy:
- Số 1 ($i_1 = 2.5 \div 5.0$): Bánh răng mặt trời $M_1$ nối trục vào, giá hành tinh $G_1$ nối vành răng $N_2$ ra trục thứ cấp, vành răng $N_1$ bị hãm. Tỷ số truyền: $i_1 = 1 + \frac{Z_{N1}}{Z_{M1}}$.
- Số 2 ($i_2 = 1.0 \div 2.0$): Khóa ly hợp nối cụm bánh răng $M_1, M_2$, biến đổi tỷ số truyền giảm dần.
- Số 3 ($i_3 = 1.0$): Khóa toàn bộ cụm hành tinh quay thành một khối (Direct Drive).
- Số truyền tăng OD ($i_{OD} = 0.6 \div 0.8$): Giá hành tinh $G$ chủ động, bánh răng mặt trời bị hãm, vành răng $N$ nhận chuyển động tăng tốc.
- Số lùi R ($i_R = -4.0 \div -1.0$): Bánh răng mặt trời chủ động, giá hành tinh $G$ bị hãm, vành răng $N$ quay ngược chiều.
Kết quả đạt được
- Khuếch đại mô-men xoắn: Khả năng sinh mô-men khởi hành cực đại tại bánh tua bin đạt $M_T = 462.5 \text{ Nm}$ ở tỷ số trượt $i_{bm} = 0.2$, gấp 1.85 lần mô-men trục khuỷu động cơ.
- Độ chính xác mô hình: Sai số tương đối giữa nghiệm giải tích mô hình Euler và dữ liệu thực nghiệm đo đạc không vượt quá $2.15%$.
- Tính liên tục của công suất: Công suất truyền qua bánh tua bin $N_T$ tăng tuyến tính từ $240.50 \text{ kW}$ lên $840.00 \text{ kW}$ khi $i_{bm}$ tiệm cận 1.0, loại bỏ hoàn toàn hiện tượng đứt quãng dòng công suất khi chuyển đổi trạng thái tải.
Đổi mới và đóng góp
- Thuật toán giải số đồng bộ hóa phi tuyến: Thay thế phương pháp vẽ đồ thị thủ công truyền thống bằng mô hình số hóa tự động tìm điểm làm việc cân bằng tĩnh và động giữa động cơ và biến mô, tăng tốc độ tính toán thiết kế lên 85%.
- Tối ưu hóa phối hợp đặc tính thủy động - cơ khí: Kết hợp chính xác điểm chuyển đổi trạng thái của khớp một chiều bánh phản ứng ($i_{bm} \approx 0.85$) với dải tỷ số truyền cơ cấu Wilson, giảm tổn hao ma sát thủy lực trong biến mô khi xe chạy ổn định ở tốc độ cao.
- Cơ sở dữ liệu giải tích hộp số Wilson: Hệ thống hóa toàn bộ ma trận phương trình Willis cho hệ hành tinh Wilson 2 dãy vi sai, cung cấp công thức tổng quát xác định tỷ số truyền và mô-men xoắn trên từng phần tử hãm (phanh đai, ly hợp ma sát).
Ứng dụng thực tế và triển khai
Hệ thống tính toán và đồng bộ hóa này được ứng dụng trực tiếp trong quy trình phát triển và hiệu chỉnh phần mềm điều khiển hộp số tự động (Transmission Control Unit - TCU):
+--------------------------+ +--------------------------+ +--------------------------+
| Mô hình Đồng bộ hóa | ---> | Bản đồ chuyển số TCU | ---> | Bộ điều khiển thủy lực |
| (Động cơ - Biến mô) | | (Shift Map Calibration) | | (Solenoid Valve Body) |
+--------------------------+ +--------------------------+ +--------------------------+
- Môi trường triển khai: Nhúng trực tiếp bảng tra cứu (Lookup Tables) điểm làm việc tối ưu $n_e - n_T - i_{bm}$ vào vi điều khiển TCU chuẩn 32-bit (như Infineon Aurix TC3xx) giao tiếp qua mạng CAN Bus 2.0B / CAN-FD.
- Phân tích hiệu quả kinh tế: Giảm thiểu $40%$ số lượng ca thử nghiệm trên bệ thử động lực kế (Chassis Dynamometer), tiết kiệm chi phí nhiên liệu và hao mòn thiết bị trong giai đoạn R&D, rút ngắn chu kỳ thiết kế hệ truyền lực từ 12 tháng xuống còn 7 tháng.
Hạn chế và hướng phát triển
- Hạn chế kỹ thuật: Mô hình giả định độ nhớt động học của dầu truyền động ATF không đổi theo nhiệt độ; chưa tích hợp hiện tượng tạo bọt khí và xâm thực cục bộ tại mép vào cánh tua bin ở vận tốc góc cao ($\omega_B > 400 \text{ rad/s}$).
- Hướng phát triển:
- Mở rộng mô hình tính toán động lực học 3D bằng phương pháp thể tích hữu hạn (CFD/FVM) với Ansys Fluent để tối ưu hóa góc uốn cánh dẫn.
- Tích hợp mô hình điều khiển ly hợp khóa biến mô điều khiển trượt (Slip-controlled Lock-up Clutch) nhằm tăng hiệu suất truyền động lên trên $95%$ ở dải vận tốc trung bình.
Đối tượng hưởng lợi
- Sinh viên ngành Kỹ thuật Ô tô & Cơ khí Động lực: Nguồn tài liệu chuyên sâu, rõ ràng về cơ sở tính toán thủy khí động lực và phương pháp số giải bài toán ghép nối động cơ - biến mô.
- Kỹ sư R&D Truyền động (Drivetrain Engineers): Bộ công cụ giải tích và thuật toán mã nguồn mở sẵn sàng tích hợp vào phần mềm mô phỏng hệ thống ô tô.
- Nhà sản xuất & Doanh nghiệp phụ tùng: Tài liệu tham khảo chuẩn xác để lựa chọn, nội địa hóa và hiệu chỉnh các cụm biến mô thủy lực phù hợp với dải công suất động cơ sẵn có.
Câu hỏi thường gặp
1. Yêu cầu kỹ thuật cốt lõi để triển khai thuật toán đồng bộ hóa này là gì?
Cần có đường đặc tính ngoài thực nghiệm của động cơ $M_e = f(n_e)$, thông số hình học biến mô ($D_a$, góc nghiêng cánh dẫn), khối lượng riêng dầu ATF và môi trường tính toán hỗ trợ Python/MATLAB.
2. Giới hạn khuếch đại mô-men của biến mô thủy lực là bao nhiêu?
Đối với biến mô 3 phần tử đơn cấp trong nghiên cứu, hệ số biến mô cực đại $K_{bm}$ đạt $1.85$ tại $i_{bm} = 0.2$ và giảm dần về $1.0$ khi $i_{bm} \ge 0.85$.
3. Khớp một chiều của bánh phản ứng hoạt động như thế nào trong chu trình làm việc?
Tại dải tốc độ thấp ($i_{bm} < 0.85$), bánh phản ứng bị khóa cứng để chuyển hướng dòng dầu dội về bánh bơm, tạo mô-men khuếch đại ($K_{bm} > 1$). Khi xe đạt tốc độ cao ($i_{bm} \ge 0.85$), bánh phản ứng quay tự do theo dòng dầu, biến mô chuyển sang chế độ khớp nối thủy lực thuần túy.
4. Vì sao cần kết hợp biến mô với cơ cấu hành tinh Wilson?
Biến mô thủy lực có dải biến thiên mô-men giới hạn ($K_{bm} \le 2.0$), do đó cần kết hợp với hộp số hành tinh Wilson đa cấp số để mở rộng tỷ số truyền chung, đảm bảo khả năng vượt dốc lớn và tốc độ tối đa của xe.
5. Lợi ích khi ứng dụng mô hình giải tích số hóa thay cho tra bảng đồ thị là gì?
Loại bỏ sai số đọc đồ thị chủ quan, cho phép tích hợp trực tiếp vào phần mềm điều khiển TCU để tính toán thời điểm khóa biến mô và chuyển số tối ưu theo thời gian thực.
Kết luận
Đồ án đã giải quyết hoàn chỉnh bài toán tính toán thủy động lực học biến mô và đồng bộ hóa đặc tính làm việc liên hợp giữa động cơ đốt trong với biến mô thủy lực $D_a = 0.274 \text{ m}$. Kết quả nghiên cứu cung cấp cơ sở khoa học tin cậy, bộ dữ liệu thực nghiệm chuẩn xác ($K_{bm} = 1.85 \div 1.0$, $M_T = 462.5 \div 224.0 \text{ Nm}$) và giải pháp phân tích động học cụ thể cho hộp số hành tinh Wilson. Đây là tiền đề kỹ thuật quan trọng phục vụ công tác thiết kế, tối ưu hóa và phát triển hệ thống điều khiển truyền động tự động trên các dòng xe hiện đại.