Tổng quan về luận án

Trong tiến trình công nghiệp hóa và hiện đại hóa, ngành công nghiệp ô tô Việt Nam đứng trước yêu cầu cấp thiết về việc nâng cao tỷ lệ nội địa hóa và tự chủ công nghệ thiết kế chế tạo theo "Chiến lược phát triển ngành công nghiệp ô tô Việt Nam đến năm 2025, tầm nhìn đến năm 2035" ban hành bởi Thủ tướng Chính phủ. Thực trạng sản xuất xe tải thương dụng trong nước kéo dài hơn hai thập kỷ qua chủ yếu dừng lại ở mức độ lắp ráp dạng CKD (Completely Knocked Down) với linh kiện ngoại nhập có chất lượng không đồng đều, thiếu hụt nền tảng tính toán thiết kế gốc. Hệ thống truyền lực (HTTL) đóng vai trò huyết mạch trong việc truyền và biến đổi mô men xoắn từ động cơ đến các bánh xe chủ động, song việc thiết kế HTTL tại Việt Nam từ trước đến nay chủ yếu dựa trên các công thức kinh nghiệm tĩnh hoặc mô men danh nghĩa tĩnh của động cơ ($M_{emax}$), bỏ qua các quá trình quá độ sinh tải trọng động cực hạn.

Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ sự thiếu vắng một hệ phương pháp luận tích hợp giữa mô hình hóa giải tích động lực học phi tuyến, kỹ thuật mô phỏng không gian năng lượng đa miền và quy trình thực nghiệm kiểm chứng trên trục quay thực tế đối với dòng xe tải thông dụng thích ứng điều kiện giao thông và tập quán vận hành đặc thù tại Việt Nam (thường xuyên quá tải và thao tác đóng ly hợp đột ngột). Nhằm giải quyết triệt để vấn đề này, luận án tiến sĩ kỹ thuật chuyên ngành Kỹ thuật Cơ khí Động lực (Mã số: 62520116) của Nghiên cứu sinh Đỗ Giao Tiến, dưới sự hướng dẫn khoa học của PGS.TS Nguyễn Trọng Hoan và TS. Nguyễn Thanh Quang tại Trường Đại học Bách Khoa Hà Nội, đã thực hiện đề tài: "Nghiên cứu tải trọng động cho thiết kế hệ thống truyền lực ô tô tải thông dụng sản xuất tại Việt Nam".

Nghiên cứu tập trung giải quyết 3 câu hỏi nghiên cứu cốt lõi:

  1. Cơ chế hình thành và quy luật lan truyền tải trọng động cực đại ($M_{max}$) trong HTTL khi khởi hành xe trong các điều kiện đóng ly hợp quá độ và quá tải là gì?
  2. Làm thế nào để xây dựng một mô hình toán học và không gian trạng thái Bond Graph đa phần tử phản ánh chính xác tương tác phi tuyến giữa các cụm chi tiết (ly hợp trượt, hộp số phân cấp, trục các đăng đàn hồi, vi sai và tương tác lốp - mặt đường)?
  3. Phổ tải trọng động thực nghiệm và hệ số xung động $K_d = M_{max}/M'_{emax}$ biến thiên ra sao theo thông số kết cấu và hành vi người lái, từ đó xác lập bộ tiêu chuẩn tải trọng phục vụ tính bền và tối ưu hóa kết cấu?

Khung lý thuyết của luận án tích hợp thuyết dao động xoắn phi tuyến cơ học máy động lực, lý thuyết truyền tải năng lượng Bond Graph của Henry M. Paynter (1959), và lý thuyết động lực học ô tô của Grishkevich (1986) cùng Xitovich (1979). Nghiên cứu thực hiện khảo sát toàn diện trên mẫu xe tải thương dụng 3 tấn nhãn hiệu LIFAN 3070G1 (tải trọng danh định 29.800 N, tự trọng 42.300 N, tổng trọng lượng 74.050 N, trang bị động cơ diesel 4102qbz công suất cực đại 81 kW tại 3.000 v/ph, mô men xoắn cực đại $M_{emax} = 290\text{ Nm}$ tại 2.000 v/ph, hộp số 5 cấp số tiến). Ý nghĩa thực tiễn của công trình là cung cấp hệ số tải trọng động chuẩn hóa, tạo bước đột phá trong việc chuyển dịch từ lắp ráp thụ động sang thiết kế chủ động các cụm truyền lực tại Việt Nam.

Literature Review và Positioning

Nghiên cứu về tải trọng động và dao động xoắn trong HTTL ô tô đã trải qua nhiều giai đoạn phát triển trên trường quốc tế với hai trường phái tiếp cận chính:

Trường phái cơ học giải tích cổ điển tập trung vào mô hình hóa tập trung khối lượng quán tính và độ đàn hồi quy dẫn. Mazziotti (1965) đặt nền móng đầu tiên khi khảo sát hiện tượng rung ồn (NVH) và dao động mô men xoắn trên trục truyền động xe tải hạng nặng. Tiếp nối hướng đi này, các nhà khoa học Xô Viết tiêu biểu như Xitovich (1979, 1985) và Grishkevich (1986) đã cụ thể hóa phương pháp xác định tải trọng cực đại khi đóng ly hợp đột ngột. Xitovich chỉ ra rằng khi nhả bàn đạp ly hợp cưỡng bức, lực quán tính tịnh tiến của đĩa ép làm gia tăng lực ép ma sát tức thời, tạo ra các đỉnh mô men xoắn cao gấp 2 đến 3 lần mô men cực đại danh nghĩa của động cơ. Grishkevich (1986) đã đề xuất hệ công thức kinh nghiệm phân kỳ theo tỷ số truyền chung ($i_{tr}$), phân lập các miền tải trọng khi $i_{tr} < 20$, $20 \le i_{tr} \le 50$, và $i_{tr} \ge 50$. Tuy nhiên, hạn chế lớn nhất của trường phái giải tích cổ điển là việc đơn giản hóa toàn bộ hệ thống về mô hình một khâu đàn hồi duy nhất, bỏ qua tính chất phân tán năng lượng, ma sát nhớt phi tuyến và hiện tượng trượt vi mô phức tạp tại các khớp nối cơ khí.

Trường phái mô hình hóa hiện đại và điều khiển tối ưu tiếp cận vấn đề thông qua các công cụ không gian trạng thái và phần tử hữu hạn. Rabeih (1997) đã chứng minh rằng mô hình rút gọn 4 khối lượng quán tính có thể phản ánh chính xác các tần số dao động riêng thấp mà không làm tổn hao độ chính xác so với mô hình bậc cao. Yuanfeng Xia và cộng sự (2005) áp dụng phương pháp phần tử hữu hạn (FEM) để giải bài toán trị riêng cho mô hình 5 bậc tự do nhằm đánh giá nguy cơ cộng hưởng xoắn. Magnus Pettersson (1997, 2003) nghiên cứu dao động xoắn phi tuyến HTTL và đề xuất thuật toán điều khiển tích cực mô men động cơ để triệt tiêu dao động giật cục (anti-jerk control), khẳng định rằng tải trọng động biên độ lớn tập trung ở miền tần số thấp (0,5 - 20 Hz từ mặt đường và quá trình đóng ly hợp), trong khi dao động tần số cao (20 - 300 Hz từ kích thích động cơ) gây mỏi và tiếng ồn. Gần đây, các công trình của Yang et al. (2007) và Choi et al. (2008) phát triển mô hình ghép bậc tự do giữa dao động xoắn HTTL và dao động thẳng đứng của hệ thống treo.

Tại Việt Nam, các nghiên cứu trước đây của Đặng Tiến Hòa (2000), Haidar Qassim (2005), Hồ Sĩ Xuân Diệu (2009), Cao Hùng Phi (2012) và Lê Văn Anh (2012) chủ yếu tập trung vào máy kéo nông nghiệp, cơ cấu vi sai, rung ồn hộp số cục bộ hoặc mài mòn bánh răng bề mặt. Công trình của Nguyễn Khắc Tuân (2011) bước đầu ứng dụng mô phỏng tải trọng động nhưng dành riêng cho cấu trúc xe hybrid.

Luận án của Đỗ Giao Tiến định vị chính xác tại giao điểm còn bỏ ngỏ: Thiết lập chuỗi phương pháp luận hoàn chỉnh từ mô hình giải tích 5 khối lượng, chuyển hóa sang mô hình mạng liên kết Bond Graph toàn diện đa cụm (xét chi tiết từng cấp số hộp số, trục các đăng, vi sai, bán trục và lốp xe), và thẩm định bằng hệ thống đo tenzo không dây trên đường thử thực tế đối với xe tải thông dụng 3 tấn. So sánh với các nghiên cứu quốc tế trên xe GAZ-53A ($K_d = 1,9 - 3,1$) và ZIL-130 ($K_d = 2,2$) của Xitovich, nghiên cứu này làm rõ tính đặc thù của động cơ diesel hiện đại kết hợp trợ lực ly hợp khí nén trong điều kiện chở quá tải 150% tại Việt Nam.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án đã mở rộng và làm sâu sắc thêm lý thuyết tải trọng động của Xitovich (1979) và Grishkevich (1986) thông qua việc phân lập và giải mã định lượng 8 trạng thái động lực học quá độ (transient dynamic states) xuất hiện khi đóng ly hợp khởi hành xe:

  1. Thiết lập mô hình giải tích 5 khối lượng quán tính ($I_1$: bánh đà và phần chủ động ly hợp; $I_2$: phần bị động ly hợp; $I_3$: bánh răng hộp số; $I_4$: bánh xe chủ động; $I_5$: khối lượng chuyển động tịnh tiến của ô tô quy dẫn) tương tác qua khâu đàn hồi trục các đăng $E$ và hai khâu liên kết ma sát phi tuyến $L_1$ (ly hợp ma sát) và $L_2$ (khâu bám lốp - mặt đường).
  2. Xây dựng ma trận điều kiện biên dựa trên quan hệ so sánh giữa vận tốc góc giới hạn của đĩa bị động $\omega_{2max} = M_c \sqrt{E/I_{23}}$ với vận tốc góc quy dẫn của bánh đà $\omega_{1max}$, kết hợp cùng ngưỡng mô men trượt giới hạn $M_{cp1}$ và $M_{cp2}$. Luận án đã chuẩn hóa hệ phương trình vi phân cho 3 nhóm trạng thái lớn:
    • Nhóm $A$ ($\omega_{2max} < \omega_{1max}$ - ly hợp $L_1$ trượt): Gồm phân nhánh $A_1$ (bánh xe bị khóa cứng, $M_{maxA1} = 2M_c$), $A_2$ (bánh xe không trượt, $M_{maxA2} \approx 2M_c$), và $A_3$ (cả $L_1$ và $L_2$ cùng trượt đồng thời, phương trình cân bằng năng lượng đàn hồi - ma sát): $$M_{maxA3} \approx M_c \frac{I_{45}}{I_{25}} + M_\varphi \frac{I_{23}}{I_{25}} + \sqrt{\left[M_c \left(2 - \frac{I_{45}}{I_{25}}\right) - M_\varphi \frac{I_{23}}{I_{25}}\right] M_c \frac{I_{45}}{I_{25}}}$$
    • Nhóm $B$ ($\omega_{2max} \ge \omega_{1max}$ - ly hợp $L_1$ dính hoàn toàn): Gồm $B_1$ (bánh xe khóa cứng, $M_{maxB1} = \omega_{1max}/\sqrt{E/I_{25}}$), $B_2$ (xe chuyển động không trượt lốp), và $B_3$ (bánh xe $L_2$ bị trượt quay trên đường).
    • Nhóm $C$ ($\omega_{2max} \ge \omega_{1max}$, chuyển trạng thái stick-slip $L_1$ từ dính sang trượt do phản lực quán tính cản vượt ngưỡng mô men ma sát định mức): Xác lập phương trình giải tích cho $C_1$ và $C_2$.

Khám phá lý thuyết quan trọng nhất là việc chứng minh mô men động cực đại thực tế không chỉ bị chặn bởi mô men ma sát tĩnh định mức $M_{c0} = \beta M_{emax}$, mà hoàn toàn có thể vượt qua ngưỡng này do tích lũy động năng xoắn và giải phóng năng lượng đàn hồi đột ngột trong khoảng thời gian siêu ngắn ($40/500\text{ s} - 250/500\text{ s}$).

Khung phân tích độc đáo

Luận án tiên phong ứng dụng phương pháp mạng liên kết (Bond Graph Theory) do Henry M. Paynter phát minh để mô hình hóa toàn diện hệ thống truyền lực ô tô tải nhiều bậc tự do. Khung phân tích được cấu trúc dựa trên việc phân rã hệ thống thành các biến trạng thái cơ bản: biến nỗ lực (Effort - $e$, tương ứng với mô men xoắn $M$) và biến dòng (Flow - $f$, tương ứng với vận tốc góc $\omega$), với tích số $P = e \cdot f$ đại diện cho dòng công suất tức thời truyền qua các liên kết cấu trúc (half-arrow bonds).

[SE: Momen Dong Co] ---> (1_truc khuyu) ---> [I: Banh Da]
                                |
                             (0_truot) ---> [R: Ma Sat Ly Hop]
                                |
                             (1_ly hop) ---> [C: Lo Xo Giam Chan]
                                |
                             [TF: Hop So 5 Cap i_hs]
                                |
                             (1_cac dang) ---> [C: Dan Hoi Truc] + [R: Can Nhot]
                                |
                             [TF: Truyen Luc Chinh i_0]
                                |
                             (0_vi sai) ---> (1_ban truc trai/phai)
                                                  |
                                               [TF: Banh Xe r_bx]
                                                  |
                                               (1_tinh tien) ---> [I: Khoi Luong Xe m] + [R: Can Lan/Gio]

Khung phân tích tích hợp 5 lớp phần tử chức năng:

  • Nguồn nỗ lực (Effort Source - $SE$): Mô tả mô men động cơ dạng chuỗi Fourier điều hòa $M(\varphi) = M_0 + \sum M_{uk} \sin(k\Omega + \varphi_k)$ đến hài bậc 4.
  • Phần tử quán tính ($I$): Biểu diễn mô men quán tính bánh đà ($I_{bd}$), cụm đĩa ma sát ($I_{lh}$), trục số ($I_{hs}$), trục các đăng ($I_{cd}$), bán trục ($I_{bt}$), bánh xe ($I_{bx}$) và khối lượng chuyển động tịnh tiến của xe ($m$).
  • Phần tử tích lũy thế năng đàn hồi ($C$): Biểu diễn độ cứng xoắn lò xo giảm chấn ly hợp ($C_{lh}$), độ cứng xoắn trục các đăng ($C_{cd}$), bán trục ($C_{bt}$) và độ đàn hồi hướng tiếp tuyến của lốp ($C_l$).
  • Phần tử tiêu tán năng lượng ($R$): Mô tả ma sát trượt Coulomb tại bề mặt đĩa ly hợp và hệ số cản nhớt ($B_{lh}, B_{hs}, B_{cd}, B_{bt}, B_{bx}$).
  • Phần tử chuyển đổi ($TF$): Biểu diễn các tỷ số truyền biến đổi động học tại hộp số ($i_{hs}$), truyền lực chính ($i_0$) và bán kính lăn bánh xe ($r_{bx}$).

Toàn bộ hệ thống phương trình trạng thái phi tuyến được thiết lập tự động từ sơ đồ Bond Graph và giải số trên môi trường Matlab-Simulink, cho phép mô phỏng chính xác các hiện tượng phi tuyến va đập và chuyển pha ma sát mà phương pháp giải tích tuyến tính truyền thống không thể thực hiện được.

Phương pháp nghiên cứu tiên tiến

Thiết kế nghiên cứu

Nghiên cứu được thiết kế theo triết lý thực chứng khoa học (positivism) kết hợp chặt chẽ giữa suy diễn lý thuyết (deductive modeling) và quy nạp thực nghiệm (empirical induction). Thiết kế nghiên cứu đa tầng (multi-level research design) được triển khai qua 3 giai đoạn đồng bộ:

[Giai doan 1: Mo hinh hoa Ly thuyet]
  - Giai tich 5 khoi luong (Dong luc hoc qua do)
  - Khong gian trang thai Bond Graph (Matlab-Simulink)
  - Phan tich tan so rieng & cong huong Fourier
             |
             v
[Giai doan 2: Thiet ke & Thuc nghiem Hien truong]
  - Xe LIFAN 3070G1 (3 muc tai: 0%, 100%, 150%)
  - Cau Tenzo khong day do Momen truc cac dang
  - Cam bien quang Omron do toc do goc & Loadcell Kubota do luc dap
             |
             v
[Giai doan 3: Kiem chung & Chuan hoa Bo Thong so]
  - So sanh doi chung Ly thuyet vs Thuc nghiem (Sai so < 10%)
  - Xac lap he so tai trong dong K_d phuc vu quy trinh thiet ke

Đối tượng nghiên cứu thực nghiệm là xe tải 3 tấn LIFAN 3070G1 nguyên bản, trong tình trạng kỹ thuật hoàn hảo, đáp ứng đầy đủ các tiêu chuẩn kiểm định an toàn kỹ thuật và bảo vệ môi trường.

Quy trình nghiên cứu rigorous

Quy trình thu thập dữ liệu thực nghiệm được thiết kế và kiểm soát vô cùng nghiêm ngặt nhằm loại trừ các sai số ngẫu nhiên và nhiễu loạn môi trường:

  • Nguyên lý và bố trí cảm biến đo mô men: Sử dụng 4 lá điện trở tenzo biến dạng (strain gauges) chuyên dụng được dán chính xác theo sơ đồ cầu Wheatstone toàn phần (full bridge) nghiêng góc $45^\circ$ so với đường sinh của trục các đăng trung gian. Sơ đồ này cho phép bù trừ triệt để ảnh hưởng của lực uốn dọc trục và sự biến thiên nhiệt độ môi trường, chỉ đo thuần túy biến dạng xoắn thuần túy.
  • Hệ thống truyền dẫn dữ liệu không dây: Để giải quyết bài toán truyền tín hiệu từ trục quay tốc độ cao (lên đến 3.000 v/ph) về máy tính trung tâm mà không bị nhiễu chổi than cơ học truyền thống, tác giả đã tích hợp module phát sóng vô tuyến vi sai gắn đồng trục cùng bộ pin nguồn độc lập, truyền dữ liệu về bộ thu cố định trên khung gầm.
  • Khuếch đại và xử lý tín hiệu: Tín hiệu sau thu nhận được đưa qua bộ khuếch đại chuyên dụng đa kênh DMC Plus kết nối máy tính tính toán, kết hợp bộ lọc thông thấp (low-pass filter) để triệt tiêu nhiễu tần số cao từ hệ thống đánh lửa và rung động vỏ xe.
  • Đo vận tốc góc và lực thao tác: Sử dụng cảm biến quang điện tử OMRON E3F-DS10C4 định vị bắt điểm phản quang trên trục để ghi nhận chính xác tức thời tốc độ quay; đồng thời gắn cảm biến đo lực chuyên dụng Kubota LBX-100L (Nhật Bản) trực tiếp trên bàn đạp ly hợp nhằm định lượng chính xác thời gian và lực tác động nhả ly hợp của lái xe ($t_c$).
  • Quy trình hiệu chuẩn (Calibration): Cầu đo tenzo trên trục các đăng được hiệu chuẩn tĩnh độc lập trước khi thử nghiệm thông qua cơ cấu tay đòn chuẩn và thiết bị gia tải lực chuẩn, xây dựng đường đặc tính chuẩn tuyến tính cao ($R^2 > 0,998$) với hệ số chuẩn tải $k_{ct} = 12,34\text{ Nm/mV}$.

Data và phân tích

Ma trận quy hoạch thực nghiệm trực giao đa yếu tố được thực hiện trên đường thử bê tông nhựa phẳng, hệ số bám cao ($\varphi = 0,75 - 0,8$), bao gồm các tổ hợp biến số vận hành:

  • Cấp số truyền hộp số: Khảo sát tại tay số 1 ($i_{hs1} = 7,66$) và tay số 2 ($i_{hs2} = 4,52$).
  • Chế độ tải trọng tĩnh: 3 mức gồm không tải (0% tải, $G_0 = 42.300\text{ N}$), đầy tải danh định (100% tải, $G_{dt} = 74.050\text{ N}$), và quá tải khắc nghiệt (150% tải, $G_{qt} = 88.950\text{ N}$).
  • Tốc độ góc ban đầu của động cơ khi đóng ly hợp: 4 mức đặc trưng ($n_e = 1.000, 1.500, 2.000, 2.300\text{ v/ph}$).
  • Thời gian đóng ly hợp ($t_c$): Biến thiên từ đóng cực nhanh/đột ngột ($t_c = 0,2 - 0,25\text{ s}$), đóng trung bình ($t_c = 0,3 - 0,35\text{ s}$), đến đóng êm dịu ($t_c > 0,5\text{ s}$).

Mỗi chế độ thực nghiệm được lặp lại tối thiểu 5 lần để lấy giá trị trung bình thống kê và loại bỏ các điểm dị biệt (outliers). Dữ liệu số học được xử lý phổ dao động và so sánh trực tiếp với kết quả chạy mô phỏng từ mô hình Bond Graph trên Matlab-Simulink.

Phát hiện đột phá và implications

Những phát hiện then chốt

Kết quả mô phỏng và đo đạc thực nghiệm đã phát hiện những quy luật động lực học mang tính đột phá:

Hệ số Tải Trọng Động (K_d) phụ thuộc Tải Trọng Xe và Cấp Số:
4.0 |                                    [Tay số 2: K_d max = 3,74]
    |                                          /
3.0 |              [Tay số 1: K_d max = 2,58] /
    |                     /                  /
2.0 |                    /                  /
    |                   /                  /
1.0 |__________________/__________________/______________
          Không tải (0%)     100% Tải         150% Tải
  1. Sự bùng nổ tải trọng động cực đại khi đóng ly hợp đột ngột: Khi khởi hành xe ở chế độ nhả bàn đạp ly hợp nhanh ($t_c = 0,25\text{ s}$), mô men xoắn tức thời xuất hiện trên trục các đăng đạt đỉnh cực đại lên tới $M_{max} = 1.683\text{ Nm}$ ở tay số 1 và $M_{max} = 1.425\text{ Nm}$ ở tay số 2 (tại chế độ quá tải 150%, $n_e = 2.000\text{ v/ph}$). So với mô men danh nghĩa lớn nhất của động cơ quy dẫn về trục các đăng, hệ số tải trọng động thực tế $K_d$ đạt giá trị từ $2,12$ đến $2,58$ ở tay số 1, và đặc biệt tăng vọt lên mức $3,15$ đến $3,74$ ở tay số 2.
  2. Quy luật nghịch đảo giữa cấp số truyền và hệ số quá tải động: Trái với nhận định trực giác cho rằng tay số thấp nhất (số 1) có tỷ số truyền lớn nhất thì chịu xung lực nguy hiểm nhất, kết quả nghiên cứu chỉ ra rằng tay số 2 mới là tay số nguy hiểm nhất về mặt hệ số tải trọng động ($K_d$). Nguyên nhân do ở tay số 1, mô men bám của bánh xe chủ động với mặt đường dễ bị vượt qua làm bánh xe bị trượt quay ($L_2$ trượt), đóng vai trò như một van an toàn tự nhiên giới hạn đỉnh mô men xoắn; ngược lại, ở tay số 2, mô men truyền xuống bánh xe nhỏ hơn giới hạn bám khiến bánh xe không trượt dính chặt với mặt đường ($L_2$ không trượt), toàn bộ động năng quán tính của bánh đà và động cơ phóng thích trực tiếp làm vặn xoắn dữ dội trục các đăng và hộp số.
  3. Ảnh hưởng phi tuyến của thời gian đóng ly hợp ($t_c$): Tải trọng động tỷ lệ nghịch phi tuyến với thời gian nhả bàn đạp ly hợp. Khi rút ngắn thời gian đóng ly hợp từ $t_c = 0,5\text{ s}$ xuống $t_c = 0,2\text{ s}$, hệ số $K_d$ tăng vọt từ $1,32$ lên $3,65$. Tuy nhiên, khi tốc độ tăng mô men ma sát đạt ngưỡng $k = 3/t_c \ge 150\text{ s}^{-1}$ (tương đương $t_c \le 0,02\text{ s}$ lý thuyết), giá trị tải trọng động tiệm cận trạng thái bão hòa, tương đương trường hợp hệ thống chịu xung lực bậc thang tức thời (step input).
  4. Phổ tần số riêng và nguy cơ cộng hưởng xoắn: Tính toán trị riêng xác định được 4 tần số dao động riêng cơ bản của HTTL: $f_1 = 3,68\text{ Hz}$, $f_2 = 12,54\text{ Hz}$, $f_3 = 45,62\text{ Hz}$, và $f_4 = 118,3\text{ Hz}$. Khi phân tích kích thích Fourier từ mô men trục khuỷu động cơ 4 kỳ 4 xi-lanh ở dải vận tốc $n_e = 2.000 - 2.300\text{ v/ph}$, tần số kích thích của hài bậc 1 ($f_{k1} = 33,3 - 38,3\text{ Hz}$) và hài bậc 2 ($f_{k2} = 66,7 - 76,7\text{ Hz}$) nằm tách biệt hoàn toàn khỏi các vùng tần số riêng nguy hiểm $f_1, f_2, f_3$, khẳng định HTTL nguyên bản không bị phá hủy bởi cộng hưởng xoắn ở chế độ làm việc định mức, song lò xo giảm chấn ly hợp cần duy trì độ cứng xoắn $C_{lh} \le 8.230\text{ Nm/rad}$ để ngăn ngừa trùng tần số khi chạy không tải (garanti).
  5. Độ tương thích cao giữa mô hình lý thuyết và thực nghiệm: Sai số tương đối giữa mô phỏng Bond Graph và số liệu thực nghiệm hiện trường tại mọi chế độ khảo sát đều nằm dưới ngưỡng $\delta \le 8,7%$, chứng minh tính vững chắc và chuẩn xác của mô hình toán.

Implications đa chiều

  • Về mặt lý thuyết cơ học: Đặt nền móng cho phương pháp phân tích liên hợp năng lượng đa miền (Multi-energy domain Bond Graph) trong khảo sát động lực học ô tô tại Việt Nam, mở rộng cơ chế giải thích hiện tượng stick-slip phi tuyến của ly hợp ma sát.
  • Về mặt thiết kế kỹ thuật: Cung cấp bộ thông số tải trọng tính toán thực tế thay thế cho các hệ số an toàn tĩnh thiếu căn cứ. Khuyến cáo các kỹ sư thiết kế cụm truyền lực xe tải nhẹ bắt buộc phải sử dụng hệ số tải trọng động thiết kế $K_d = 2,5 - 2,8$ cho các cụm trục sau hộp số khi tính toán theo tay số 1, và $K_d = 3,5 - 3,8$ khi tính toán kiểm nghiệm bền mỏi và bền dập cục bộ cho trục các đăng và bán trục ở tay số 2.
  • Về mặt chế tạo và tối ưu hóa: Đưa ra giải pháp giảm tải trọng động bằng cách tối ưu hóa đặc tính phi tuyến của lò xo giảm chấn đĩa ma sát ly hợp (giảm độ cứng xoắn kết hợp tăng góc trượt giảm chấn) giúp hấp thu đỉnh xung lực mô men khi khởi hành đột ngột.
  • Về chính sách và kiểm định: Cung cấp luận cứ khoa học để Cục Đăng kiểm Việt Nam và các cơ quan quản lý ban hành quy chuẩn thử nghiệm độ bền quá tải động lực học đối với ô tô sản xuất, lắp ráp trong nước.

Limitations và Future Research

Mặc dù đạt được những kết quả xuất sắc, luận án thẳng thắn chỉ rõ 3 giới hạn nghiên cứu:

  1. Giới hạn không gian chuyển động: Mô hình Bond Graph tập trung chủ yếu vào dao động xoắn thuần túy của hệ thống truyền lực khi xe chạy thẳng trên đường bằng phẳng, chưa tích hợp mô hình không gian 3D xét đến tương tác liên hợp với dao động lắc ngang, dao động mấp mô thẳng đứng từ hệ thống treo và lực quán tính khi xe quay vòng.
  2. Giả định ma sát lốp - mặt đường: Mô hình tương tác bánh xe - mặt đường mới chỉ áp dụng hệ số bám Coulomb đơn giản hóa với ngưỡng bám không đổi, chưa tích hợp mô hình lốp phi tuyến phức hợp Magic Formula (Pacejka) xét đến độ trượt biến thiên liên tục theo nhiệt độ và độ ẩm mặt đường.
  3. Mẫu khảo sát thực nghiệm: Nghiên cứu thực nghiệm tập trung chi tiết trên một dòng xe đại diện là LIFAN 3070G1 (3 tấn), chưa có điều kiện mở rộng đo kiểm đồng loạt trên các phân khúc xe tải trung và xe tải nặng 3 - 4 chân.

Các hướng phát triển nghiên cứu tiếp theo bao gồm:

  • Mở rộng mô hình Bond Graph liên kết không gian đa chiều giữa HTTL, hệ thống treo và hệ thống lái (Full Vehicle Powertrain-Chassis Dynamic Model).
  • Nghiên cứu ứng dụng các bộ giảm chấn xoắn thủy lực hoặc ly hợp điều khiển điện tử (Dual-Clutch / Automated Manual Transmission) với thuật toán nhả ly hợp thông minh chống giật.
  • Ứng dụng mô hình tải trọng động của luận án vào bài toán dự báo tuổi thọ mỏi tích lũy (Cumulative Fatigue Life Prediction) của hệ thống bánh răng hộp số và trục các đăng theo tiêu chuẩn ISO 6336.

Tác động và ảnh hưởng

Luận án đóng góp trực tiếp vào sự chuyển dịch công nghệ của ngành công nghiệp ô tô nội địa:

  • Ảnh hưởng học thuật: Thiết lập chuẩn mực phương pháp luận nghiên cứu thực nghiệm đo mô men xoắn động trên trục quay không dây cho các viện nghiên cứu và trường đại học kỹ thuật trên cả nước. Mô hình Bond Graph của luận án trở thành tài liệu giảng dạy chuyên sâu trong chương trình đào tạo sau đại học ngành Kỹ thuật Cơ khí Động lực.
  • Chuyển giao và sản xuất công nghiệp: Kết quả nghiên cứu đã được chia sẻ và ứng dụng trực tiếp tại Nhà máy ô tô Cổ Loa, Tổng công ty VEAM và các đơn vị sản xuất xe tải thương dụng, phục vụ trực tiếp cho các đề tài khoa học trọng điểm cấp Nhà nước (Chương trình KC.05) về thiết kế, chế tạo cụm hộp số và cầu sau xe tải đến 3 tấn.
  • Hiệu quả kinh tế - xã hội: Giúp các doanh nghiệp tối ưu hóa kích thước kết cấu chi tiết máy, giảm thiểu khối lượng kim loại dư thừa nhưng vẫn đảm bảo tuyệt đối hệ số an toàn chịu tải, giảm thiểu tỷ lệ gãy trục các đăng và vỡ bánh răng vi sai khi xe vận hành trong điều kiện quá tải tại Việt Nam, tiết kiệm hàng chục tỷ đồng chi phí bảo hành và sửa chữa hàng năm.

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

  • Nghiên cứu sinh và Giảng viên ngành Cơ khí Động lực: Kế thừa hệ phương pháp luận mô hình hóa Bond Graph và quy trình thực nghiệm đo động lực học trục quay chính xác cao.
  • Kỹ sư R&D tại các nhà máy ô tô (VEAM, Thaco, TMT, VinFast): Sử dụng trực tiếp bộ thông số hệ số tải trọng động $K_d$ và ma trận độ cứng xoắn chuẩn để tính toán, thiết kế mới hoặc cải tiến nội địa hóa HTTL.
  • Các nhà hoạch định chính sách và Quản lý kỹ thuật: Cơ sở dữ liệu khoa học tin cậy phục vụ xây dựng tiêu chuẩn kỹ thuật quốc gia (TCVN) về an toàn kết cấu và độ bền mỏi hệ thống truyền động ô tô thương mại.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và mở rộng lý thuyết nào? Trả lời: Luận án đã mở rộng xuất sắc lý thuyết tải trọng động của Xitovich (1979) bằng cách xây dựng hệ 8 trạng thái động lực học phân nhánh có xét đến tương tác chuyển trạng thái dính - trượt (stick-slip) đồng thời giữa ly hợp ma sát ($L_1$) và bánh xe chủ động ($L_2$), tích hợp trọn vẹn trong không gian trạng thái Bond Graph đa miền năng lượng.

2. Điểm cải tiến cốt lõi của phương pháp nghiên cứu so với các công trình trước đây? Trả lời: So với các nghiên cứu của Grishkevich (chỉ tính gần đúng giải tích 1 khâu đàn hồi) và Cao Hùng Phi (chỉ khảo sát rung ồn tĩnh), luận án đã đột phá bằng việc thiết lập hệ thống đo biến dạng Tenzo không dây thời gian thực trên trục các đăng quay kết hợp bộ khuếch đại DMC Plus và cảm biến lực bàn đạp Kubota, giải quyết triệt để bài toán đo tải trọng động tức thời trong điều kiện xe chạy thực tế trên đường.

3. Phát hiện nào bất ngờ nhất từ số liệu thực nghiệm và giải thích cơ học? Trả lời: Phát hiện bất ngờ nhất là hệ số tải trọng động $K_d$ ở tay số 2 ($K_{d max} = 3,74$) cao hơn đáng kể so với tay số 1 ($K_{d max} = 2,58$). Về mặt cơ học, ở tay số 1 mô men xoắn vượt quá giới hạn bám gây trượt lốp ($L_2$ trượt) giúp giải phóng năng lượng quá tải; trong khi ở tay số 2, bánh xe bám chặt mặt đường ($L_2$ dính) buộc toàn bộ năng lượng va đập quán tính phải chuyển hóa thành biến dạng xoắn đàn hồi trên trục truyền động.

4. Luận án có cung cấp quy trình tái lập thực nghiệm (Replication Protocol) không? Trả lời: Hoàn toàn có. Luận án trình bày chi tiết sơ đồ đấu nối cầu Wheatstone tenzo, thông số kỹ thuật cảm biến Omron, cảm biến lực Kubota LBX-100L, công thức lấy chuẩn tải $k_{ct} = 12,34\text{ Nm/mV}$, quy trình dán cảm biến góc $45^\circ$, và ma trận thực nghiệm trực giao đảm bảo khả năng tái lập độc lập 100%.

5. Chương trình nghiên cứu 10 năm tiếp theo được định hình ra sao? Trả lời: Định hình 3 mũi nhọn: (1) Tích hợp mô hình Bond Graph toàn xe 3D đa vật thể; (2) Tự động hóa điều khiển quá trình đóng ly hợp thông minh thích ứng tải trọng trên hộp số AMT; (3) Thiết lập bản đồ tải trọng mỏi tích lũy đa trục theo chu kỳ vận hành thực tế tại các vùng địa hình Việt Nam.

Kết luận

  1. Luận án đã giải quyết trọn vẹn và triệt để bài toán xác định tải trọng động cực hạn tác dụng lên hệ thống truyền lực ô tô tải thông dụng sản xuất tại Việt Nam, kết hợp nhuần nhuyễn giữa lý thuyết cơ học hiện đại và thực nghiệm hiện trường tiên tiến.
  2. Ứng dụng thành công và chuẩn hóa phương pháp mạng liên kết Bond Graph để mô hình hóa toàn diện các cụm chi tiết phi tuyến trong HTTL, mở ra hướng đi mới trong việc mô phỏng động lực học cơ khí phức hợp.
  3. Xác lập chính xác quy luật biến thiên của hệ số tải trọng động $K_d$, chứng minh tay số 2 là chế độ chịu quá tải động nguy hiểm nhất ($K_d$ đạt tới $3,74$ khi chở quá tải 150% và đóng ly hợp nhanh), đập tan giả định thiết kế truyền thống vốn chỉ tập trung vào tay số 1.
  4. Xây dựng bộ dữ liệu khoa học thực nghiệm có độ tin cậy cao trên xe tải LIFAN 3070G1 với sai số mô phỏng - thực nghiệm dưới 8,7%, làm chủ công nghệ đo lường tenzo không dây trên trục quay tốc độ cao.
  5. Cung cấp bộ thông số tải trọng tính toán thiết kế chuẩn mực phục vụ trực tiếp cho các nhà máy sản xuất ô tô trong nước (VEAM, Vinaxuki, Thaco), đóng góp thiết thực vào mục tiêu nâng cao hàm lượng khoa học công nghệ và tự chủ thiết kế chế tạo phương tiện giao thông mang thương hiệu Việt Nam.