Tổng quan về luận án

Cơ giới hóa nông nghiệp và lâm nghiệp đóng vai trò then chốt trong việc tối ưu hóa chuỗi giá trị sản xuất, giảm thiểu tổn thất sau thu hoạch và nâng cao năng suất lao động. Theo số liệu của Tổng cục Thống kê (2021), cả nước hiện có hơn 400.000 máy kéo nông nghiệp đa năng với tổng công suất xấp xỉ 5 triệu mã lực, ghi nhận mức tăng trưởng 60% trong giai đoạn 2011–2021. Tuy nhiên, tính chất mùa vụ đặc thù khiến thời gian vận hành hữu ích của máy kéo đa năng chỉ đạt khoảng 2–3 tháng mỗi năm, chủ yếu phục vụ khâu làm đất, dẫn đến hiệu suất khai thác tài sản và thu hồi vốn đầu tư rất thấp. Việc chuyển đổi máy kéo bánh lốp đa năng thành liên hợp máy (LHM) kéo – xúc lật thông qua việc lắp đặt bộ công tác xúc lật dạng module phía trước (điển hình trên các dòng máy công suất trung bình 30–40 HP như Yanmar 3000 lắp bộ xúc Y3510FLH) là giải pháp đột phá nhằm mở rộng dải công tác sang các nhiệm vụ bốc xếp vật liệu rời (đất, cát, sỏi, ngũ cốc, phân bón), san lấp mặt bằng và dọn dẹp trang trại.

Vấn đề cốt lõi đặt ra là việc bổ sung cụm xúc lật phía trước và đối trọng phía sau làm thay đổi căn bản đặc trưng khối lượng, dịch chuyển trọng tâm tổng thể lên cao và về phía trước, biến đổi mô men quán tính khối lượng theo không gian 3 chiều ($J_x, J_y, J_z$) phụ thuộc trực tiếp vào góc nâng cần $\alpha$ và mức tải $k_Q$. Trong điều kiện di chuyển mang tải trên nền đường mấp mô không bằng phẳng, sự phân bố lại tải trọng làm tăng tải đột biến lên cầu trước (cầu dẫn hướng) và giảm tải trọng bám của cầu sau (cầu chủ động), dễ dẫn đến mất ổn định hướng, trượt quay, tách bánh khỏi mặt đất hoặc lật xe nguy hiểm. Các nghiên cứu trước đây tại Việt Nam và quốc tế hầu như chỉ xem xét máy kéo chạy không tải, máy xúc chuyên dụng có kết cấu cầu trước cứng vững (không có khớp cân bằng), hoặc khảo sát ổn định của gầu xúc ở trạng thái đứng yên mà chưa giải quyết trọn vẹn bài toán động lực học không gian của LHM máy kéo – xúc lật khi di chuyển mang tải trên nền mấp mô.

Nghiên cứu tập trung giải quyết hệ thống câu hỏi và giả thuyết khoa học:

  • Câu hỏi nghiên cứu 1 (RQ1): Các yếu tố kết cấu (vị trí lắp đặt, khối lượng đối trọng $Q_{DT}$, góc nâng cần $\alpha$) và điều kiện vận hành (vận tốc $V$, tải trọng $Q_{VL}$, biên độ mấp mô mặt đường $h_0$) tương tác như thế nào đối với các chỉ tiêu ổn định chuyển động và khả năng điều khiển của LHM?
  • Câu hỏi nghiên cứu 2 (RQ2): Quy luật phân bố tải trọng pháp tuyến động và phản lực tiếp xúc bánh xe – mặt đường diễn biến ra sao khi LHM vượt qua mấp mô đơn và mấp mô tuần hoàn có độ lệch pha không gian?
  • Giả thuyết nghiên cứu 1 (H1): Việc xác lập cấu hình tải trọng hợp lý với đối trọng $Q_{DT} \ge 450\text{ kg}$, góc nâng cần $\alpha \in [-40^\circ, -30^\circ]$ ở dải tốc độ $V < 3\text{ km/h}$ sẽ duy trì phản lực pháp tuyến tiếp xúc bánh xe liên tục ($P_z > 0$), ngăn chặn triệt để hiện tượng mất ổn định lật ngang và trượt hướng.
  • Giả thuyết nghiên cứu 2 (H2): Mô hình động lực học không gian tích hợp khớp liên kết thăng bằng cầu trước và hàm kích thích mấp mô tuần hoàn lệch pha giải tích cho phép mô phỏng chính xác đáp ứng dao động không gian với hệ số tương quan thực nghiệm $R^2 > 0.95$.

Khung lý thuyết của luận án được xây dựng dựa trên cơ học hệ nhiều vật (Multi-body Dynamics), nguyên lý biến phân D'Alembert – Lagrange loại hai và lý thuyết tương tác bánh lốp – mặt đường (Terramechanics). Nghiên cứu được triển khai thực nghiệm tại Viện Phát triển Công nghệ Cơ – Điện và Bộ môn Động lực, Khoa Cơ – Điện, Học viện Nông nghiệp Việt Nam trên đối tượng thực là máy kéo Yanmar 3000 (động cơ 3T84, 30 HP, khối lượng bản thân 1.420 kg, chiều dài cơ sở 1.845 mm) tích hợp bộ xúc lật Y3510FLH (sức nâng định mức 450 kg, dung tích gầu $0,2\text{ m}^3$, chiều cao nâng cực đại 2.700 mm).

Literature Review và Positioning

Nghiên cứu động lực học và ổn định của máy di động trên nền mấp mô được hình thành từ nhiều nhánh học thuật chuyên sâu. Nhánh thứ nhất tập trung vào máy xúc lật chuyên dụng: Ikonen (2006) thiết lập mô hình phẳng 2 bậc tự do (2-DOF) nửa xe theo phương dọc để khảo sát dao động khi nâng hạ tải. Debeleac (2007) mở rộng thành mô hình không gian 5 bậc tự do (5-DOF) mô tả rung động kết cấu máy xúc chuyên dụng, và tiếp tục nâng cấp lên mô hình phẳng 3 bậc tự do (Debeleac, 2014) khảo sát hệ số tải trọng động của cần và gầu. Popescu và cộng sự (2004, 2009) mô phỏng động học dẫn động thủy lực IF-65 lắp trên máy kéo U-650 và đánh giá ổn định dọc khi phanh trên dốc. Nhánh thứ hai nghiên cứu động lực học máy kéo nông nghiệp: Zhen Li và cộng sự (2015) thiết lập mô hình không gian 9 bậc tự do phân tích chuyển động trên dốc ngang; Ahmadi (2014) khảo sát dao động toàn xe có tính đến khớp lắc thăng bằng cầu trước. Tại Việt Nam, Võ Văn Hường (2014), Nguyễn Ngọc Tú (2016), Tạ Tuấn Hưng (2017) và Trần Văn Tùng (2017) đã đi sâu vào động lực học đoàn xe kéo moóc và ô tô tải nhiều cầu trên địa hình lâm nghiệp.

                      ┌─────────────────────────────────────────────────────────┐
                      │    HỆ THỐNG MÁY CHUYÊN DỤNG & ĐOÀN XE NÔNG - LÂM NGHIỆP  │
                      └────────────────────────────┬────────────────────────────┘
                                                   │
         ┌─────────────────────────────────────────┴─────────────────────────────────────────┐
         ▼                                                                                   ▼
┌───────────────────────────────────────────────────┐               ┌───────────────────────────────────────────────────┐
│     MÁY XÚC LẬT CHUYÊN DỤNG & HỆ CÔNG TÁC         │               │         MÁY KÉO BÁNH LỐP & ĐOÀN XE KÉO            │
│  - Ikonen (2006): Mô hình dọc 2-DOF               │               │  - Zhen Li et al. (2015): Mô hình không gian 9-DOF│
│  - Debeleac (2007, 2014): Mô hình 5-DOF & 3-DOF   │               │  - Ahmadi (2014): Khớp lắc thăng bằng cầu trước   │
│  - Popescu et al. (2004, 2009): Nâng hạ & phanh   │               │  - Võ Văn Hường, Nguyễn Ngọc Tú, Tạ Tuấn Hưng:    │
│  * Hạn chế: Khung vỏ cứng, thiếu khớp lắc thăng   │               │    Động lực học đoàn xe, ổn định quay vòng        │
│    bằng, gầu cố định hoặc khảo sát tĩnh.          │               │  * Hạn chế: Không tích hợp bộ công tác xúc tải    │
│                                                   │               │    phía trước, thông số khối lượng cố định.       │
└─────────────────────────┬─────────────────────────┘               └─────────────────────────┬─────────────────────────┘
                          │                                                                   │
                          └─────────────────────────────────┬─────────────────────────────────┘
                                                            │
                                                            ▼
                                ┌───────────────────────────────────────────────────────┐
                                │          KHOẢNG TRỐNG KHOA HỌC (RESEARCH GAP)         │
                                │   Thiếu mô hình động lực học không gian tích hợp      │
                                │   đồng thời: Khớp thăng bằng cầu trước + Khối lượng   │
                                │   & quán tính biến thiên đa cấp (m_k, J_i) + Kích thích│
                                │   mấp mô không gian 3 chiều có lệch pha tiếp xúc điểm. │
                                └───────────────────────────┬───────────────────────────┘
                                                            │
                                                            ▼
                                ┌───────────────────────────────────────────────────────┐
                                │     ĐÓNG GÓP ĐỘT PHÁ CỦA LUẬN ÁN (PHẠM TRỌNG PHƯỚC)    │
                                │  1. Mô hình không gian tích hợp tham số khối lượng    │
                                │     biến thiên theo góc nâng α và mức tải k_Q.        │
                                │  2. Kích thích mấp mô sin lệch pha (S_P = S_T, d).     │
                                │  3. Trắc diện ảo cung tròn thực nghiệm tiếp xúc điểm. │
                                │  4. Bộ thông số vận hành an toàn chuẩn hóa.           │
                                └───────────────────────────────────────────────────────┘

Trong bức tranh tổng quan đó tồn tại một cuộc tranh luận học thuật sâu sắc giữa hai trường phái tiếp cận:

  1. Trường phái mô hình hóa tổng quát nguyên khối: Giả định hệ thống máy là một khối rắn duy nhất có tọa độ trọng tâm cố định nhằm đơn giản hóa việc giải ma trận phương trình vi phân dao động (điển hình trong các phân tích của Popescu et al., 2009 và Debeleac, 2014). Hạn chế lớn của trường phái này là triệt tiêu hoàn toàn sự biến thiên của tọa độ trọng tâm cục bộ và ten-xơ quán tính khi nâng hạ cần, dẫn đến sai số tích lũy nghiêm trọng khi đánh giá ổn định lật động.
  2. Trường phái cơ hệ nhiều vật ghép nối động lực: Đòi hỏi mô hình hóa chi tiết từng bậc tự do của cơ cấu nâng, xi lanh thủy lực và biến dạng khung (như tiếp cận của Zhen Li et al., 2015). Tuy nhiên, cách tiếp cận này lại đưa vào quá nhiều tham số không xác định (nhiễu tham số), làm phức tạp hóa quá trình kiểm chứng thực nghiệm trên thực địa.

Luận án đã định vị chính xác khoảng trống nghiên cứu bằng cách dung hòa hai cách tiếp cận: xây dựng mô hình hệ nhiều vật có cấu trúc tham số hóa theo từng trạng thái công tác cố định (State-dependent Parameter Lumped Model). Theo đó, LHM được xem xét như một máy kéo đặc biệt với các đặc trưng quán tính ($m_k, a, b, h, J_x, Jy, J_z$) biến thiên liên tục theo từng chế độ tải ($k_Q$) và vị trí góc nâng cần ($\alpha$), đồng thời tích hợp đầy đủ khớp lắc thăng bằng cầu trước quanh trục dọc $O_1x_1$. Đây là bước tiến vượt bậc so với mô hình của Debeleac (2007) vốn áp dụng cho khung xe cứng của máy xúc chuyên dụng, và hoàn thiện hơn mô hình của Zhen Li (2015) nhờ bổ sung đầy đủ cụm công tác xúc lật và đối trọng.

Đó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 lý thuyết động lực học máy kéo nông nghiệp thông qua việc tái cấu trúc hệ phương trình vi phân chuyển động phi tuyến Lagrange loại hai cho cơ hệ có liên kết động học phi đối xứng:

$$\frac{d}{dt}\left(\frac{\partial L}{\partial \dot{q}_i}\right) - \frac{\partial L}{\partial q_i} = Q_i$$

Hệ tọa độ suy rộng được xác lập đầy đủ bao gồm các chuyển vị thẳng đứng ($z$), góc lắc dọc ($\varphi$), góc lắc ngang thân xe ($\theta$) và góc lắc ngang độc lập của cầu trước quanh chốt thăng bằng ($\theta_1$). Ma trận quán tính khối lượng được xây dựng lại hoàn toàn bằng cách tích hợp ten-xơ quán tính của 4 cụm phần tử chính: thân máy kéo cơ sở ($m_K$), cụm cần nâng ($m_C$), cụm gầu xúc chứa vật liệu ($m_G = m_{G0} + m_{VL}$) và cụm đối trọng đuôi ($m_{DT}$):

$$J_x = J_{xk} + m_k r_{xk}^2 + J_{xC} + m_C r_{xC}^2 + J_{xG} + m_G r_{xG}^2 + J_{xDT} + m_{DT} r_{xDT}^2$$

$$J_y = J_{yk} + m_k r_{yk}^2 + J_{yC} + m_C r_{yC}^2 + J_{yG} + m_G r_{yG}^2 + J_{yDT} + m_{DT} r_{yDT}^2$$

Mô hình lý thuyết này tạo ra một bước chuyển đổi mô hình (paradigm shift) trong tính toán ổn định: thay vì xem xét các chỉ tiêu ổn định tĩnh tĩnh bức thông qua góc nghiêng giới hạn $\alpha_{\lim}$, luận án chuyển sang đánh giá ổn định động thông qua tiêu chí duy trì phản lực pháp tuyến dương liên tục tại 4 điểm tiếp xúc bánh xe ($P_{z1}, P_{z2}, P_{z3}, P_{z4} > 0$) dưới tác động của kích thích động học không gian ngẫu nhiên và tuần hoàn.

                          ┌─────────────────────────────────────────────────────────────┐
                          │     HỆ THỐNG PHƯƠNG TRÌNH ĐỘNG LỰC HỌC KHÔNG GIAN TỔNG QUÁT   │
                          └──────────────────────────────┬──────────────────────────────┘
                                                         │
               ┌─────────────────────────────────────────┼─────────────────────────────────────────┐
               ▼                                         ▼                                         ▼
┌───────────────────────────────┐         ┌───────────────────────────────┐         ┌───────────────────────────────┐
│     ĐẶC TRƯNG HÌNH HỌC        │         │      ĐẶC TRƯNG ĐỘNG LỰC       │         │     TƯƠNG TÁC TIẾP XÚC        │
│        VÀ QUÁN TÍNH           │         │       KHỚP CẦU TRƯỚC          │         │        BÁNH LỐP - ĐẤT         │
│  - Trọng tâm tổng thể:        │         │  - Chuyển vị góc lắc độc lập  │         │  - Lực đàn hồi và giảm chấn:  │
│    a(α, k_Q), b, h            │         │    cầu trước: θ_1 quanh O_1x_1│         │    F_zi = C_i·Δz_i + k_i·Δż_i │
│  - Ten-xơ quán tính:          │         │  - Khử mô men lật ngang truyền│         │  - Hàm kích thích mặt đường:  │
│    J_x(α), J_y(α), J_z(α)     │         │    lên thân xe tại cầu trước  │         │    q(t) = h_0·sin(ωt + φ_0)   │
└──────────────┬────────────────┘         └──────────────┬────────────────┘         └──────────────┬────────────────┘
               │                                         │                                         │
               └─────────────────────────────────────────┼─────────────────────────────────────────┘
                                                         │
                                                         ▼
                                ┌─────────────────────────────────────────────────┐
                                │     HỆ PHƯƠNG TRÌNH VI PHÂN CHUYỂN ĐỘNG (4.1)   │
                                │   [M(α, k_Q)]·{q̈} + [C]·{q̇} + [K]·{q} = {F(t)}  │
                                │   Tọa độ suy rộng: q = [z, φ, θ, θ_1]^T         │
                                └────────────────────────┬────────────────────────┘
                                                         │
                                                         ▼
                                ┌─────────────────────────────────────────────────┐
                                │        CHỈ TIÊU ĐÁNH GIÁ ỔN ĐỊNH ĐỘNG HỌC       │
                                │  - Ổn định tiếp xúc: P_zi > 0 (chống tách bánh) │
                                │  - Ổn định lật ngang: θ < θ_lim, α_lim          │
                                │  - Khả năng điều khiển lái: P_zn ≥ P_zn_cp      │
                                └─────────────────────────────────────────────────┘

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

Khung phân tích của luận án tích hợp đồng thời 3 trụ cột lý thuyết:

  1. Lý thuyết động lực học cơ hệ đa vật liên kết đàn hồi: Mô tả các chuyển vị tương đối của thân xe, cầu xe và bánh xe thông qua hệ số độ cứng lốp ($C_i$) và hệ số cản giảm chấn ($k_i$).
  2. Lý thuyết kích thích động học mặt đường không gian hai vệt bánh: Luận án đề xuất phương pháp kích thích mấp mô hình sin có cùng bước sóng ($S_P = S_T$) nhưng lệch pha nhau một khoảng $d$ ($0 \le d \le 0.5 S_P$):

$$q_P(t) = h_{0P} \sin\left(\frac{2\pi V}{S_P} t\right)$$

$$q_T(t) = h_{0T} \sin\left(\frac{2\pi V}{S_T} t + \frac{2\pi d}{S_P}\right)$$

Giải pháp này tạo ra đồng thời dao động uốn dọc ($\varphi$) và xoắn ngang ($\theta$) tuần hoàn, cho phép khảo sát tách biệt các tần số cộng hưởng nguy hiểm của hệ thống. 3. Lý thuyết bám – trượt Terramechanics: Xác lập quan hệ giữa phản lực tiếp tuyến ($P_k$), lực cản lăn ($P_f = f \cdot P_z$), mô men cản quay vòng ($M_{cq}$) và điều kiện bám ($P_k \le \varphi \cdot P_z$) trong điều kiện hệ số bám thấp của nền đất nông nghiệp.

Điều kiện biên của mô hình được xác định chặt chẽ: góc lắc của chốt thăng bằng cầu trước dao động trong giới hạn hình học cho phép; lốp xe duy trì tiếp xúc điểm tại đáy vết biến dạng đàn hồi; vật liệu trong gầu xúc được coi là khối đồng nhất gắn cứng với vỏ gầu ở từng chu kỳ khảo sát.

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

Thiết kế nghiên cứu

Luận án tuân thủ thế giới quan thực chứng hậu kỳ (Post-positivism) kết hợp chặt chẽ giữa suy diễn toán học (mô hình hóa giải tích) và quy nạp thực nghiệm (kiểm chứng hiện trường). Thiết kế nghiên cứu được phân tầng 3 cấp độ:

  • Cấp độ 1 (Tương tác vi mô lốp – nền): Xác định thực nghiệm hệ số lực bám $\varphi$, hệ số cản lăn $f$ trên hai bề mặt thực tế: nền đất nông nghiệp và đường bê tông.
  • Cấp độ 2 (Đặc trưng khối lượng & Quán tính vĩ mô): Đo đạc chu kỳ dao động tự do $T_x, T_y, T_z$ của máy kéo cơ sở thông qua hệ thống giá treo dao động chuyên dụng để suy biến ten-xơ quán tính thực tế ($J_{xk}, J_{yk}, J_{zk}$).
  • Cấp độ 3 (Động lực học toàn hệ thống): Mô phỏng số trên môi trường Matlab-Simulink kết hợp chạy thực nghiệm LHM vượt mấp mô đơn nhân tạo trên bãi thử chuyên dụng.

Quy trình nghiên cứu rigorous

Quy trình thu thập và xử lý dữ liệu được thiết lập với độ chuẩn xác cơ học cao:

  • Xác định ten-xơ quán tính máy kéo: Sử dụng phương pháp dao động con lắc vật lý, ghi nhận tín hiệu dao động góc quanh trục $Y$ cho chu kỳ $T_y = 0.7\text{ s}$, quanh trục $X$ cho chu kỳ $T_x = 2.72\text{ s}$, và quanh trục $Z$ cho chu kỳ $T_z = 6.7\text{ s}$.
  • Hệ thống cảm biến đo đạc không điện:
    • Cảm biến khoảng cách laser quang học gắn tại các góc dầm cầu và khung xe đo chuyển vị thẳng đứng $z_i$ và góc lắc $\theta, \varphi$.
    • Cảm biến áp suất dầu thủy lực và lưu lượng kế giám sát trạng thái xi lanh nâng hạ.
    • Cảm biến tốc độ vòng quay bánh xe và góc quay trục lái.
  • Thu thập và xử lý tín hiệu: Toàn bộ tín hiệu điện áp từ cảm biến được đưa qua bộ thu thập dữ liệu đa kênh (Data Acquisition Unit - DAQ), kết nối máy tính qua giao thức tốc độ cao và xử lý thời gian thực bằng phần mềm chuyên dụng DasyLab 10.
  • Phương pháp trắc diện ảo cung tròn đột phá: Khi bánh xe trước ($r_{k1}$) và bánh xe sau ($r_{k2}$) có kích thước chênh lệch lớn lăn qua cùng một chướng ngại vật mấp mô đơn, tâm bánh xe sẽ dịch chuyển theo quỹ đạo khác nhau. Luận án đề xuất mô hình mấp mô thực nghiệm dạng mặt cong cung tròn, sử dụng một trắc diện ảo tương đương trong mô hình toán để đồng nhất hóa điều kiện tiếp xúc điểm độc lập với đường kính bánh xe, loại bỏ hoàn toàn sai số hình học tiếp xúc của các phương pháp thử nghiệm truyền thống.
┌──────────────────────────────────────────────────────────────────────────────────────────────────┐
│                             QUY TRÌNH THỰC NGHIỆM & KIỂM CHỨNG MÔ HÌNH                            │
└────────────────────────────────────────────────┬─────────────────────────────────────────────────┘
                                                 │
            ┌────────────────────────────────────┴────────────────────────────────────┐
            ▼                                                                         ▼
┌───────────────────────────────────────┐                 ┌───────────────────────────────────────┐
│     KHẢO SÁT THÔNG SỐ CƠ SỞ & TẢI     │                 │       ĐO ĐẠC ĐẶC TRƯNG QUÁN TÍNH      │
│  - Yanmar 3000: m = 1420 kg, L = 1845 │                 │  - Dao động con lắc vật lý đo T_i:    │
│  - Bộ xúc Y3510FLH: m_nâng = 450 kg   │                 │    * Trục Y: T_y = 0,70 s             │
│  - Nền thử nghiệm: Đất ruộng & bê tông│                 │    * Trục X: T_x = 2,72 s             │
│  - Đo lực cản lăn f & lực bám phi     │                 │    * Trục Z: T_z = 6,70 s             │
└───────────────────┬───────────────────┘                 └───────────────────┬───────────────────┘
                    │                                                         │
                    └────────────────────────────┬────────────────────────────┘
                                                 │
                                                 ▼
                    ┌─────────────────────────────────────────────────────────┐
                    │     BÀI TOÁN MÔ PHỎNG SỐ TRÊN MATLAB - SIMULINK         │
                    │  - Giải hệ vi phân chuyển động 4 bậc tự do              │
                    │  - Kích thích: Mấp mô đơn & Mấp mô sin lệch pha         │
                    │  - Xuất đáp ứng: Gia tốc z̈, φ̈, θ̈, phản lực P_zi         │
                    └────────────────────────────┬────────────────────────────┘
                                                 │
                                                 ▼
                    ┌─────────────────────────────────────────────────────────┐
                    │       THỰC NGHIỆM ĐỐI CHỨNG VỚI HỆ THỐNG DASYLAB 10     │
                    │  - Cảm biến laser đo chuyển vị & góc nghiêng khung      │
                    │  - Trắc diện thực nghiệm cung tròn (tiếp xúc điểm)      │
                    │  - Chạy thử ở cấp số truyền 1 và cấp số truyền 2        │
                    └────────────────────────────┬────────────────────────────┘
                                                 │
                                                 ▼
                    ┌─────────────────────────────────────────────────────────┐
                    │       ĐÁNH GIÁ ĐỘ TIN CẬY & TỐI ƯU HÓA THÔNG SỐ         │
                    │  - Hệ số tương quan lý thuyết - thực nghiệm:            │
                    │    R^2 = 0,97 ÷ 0,98 (Độ tin cậy rất cao)               │
                    │  - Xác lập biên thông số: Q_DT ≥ 450 kg, V < 3 km/h     │
                    └─────────────────────────────────────────────────────────┘

Data và phân tích

Dữ liệu mô phỏng giải hệ phương trình vi phân chuyển động được thực hiện trên phần mềm Matlab-Simulink với thuật toán tích phân Runge-Kutta bậc 4/5 (ode45) thích ứng bước thời gian $\Delta t = 0.001\text{ s}$.

Kết quả kiểm chứng giữa đường cong chuyển vị lý thuyết và đồ thị đo đạc thực nghiệm ở các cấp số truyền (Số truyền 1: $V = 1.6\text{ km/h}$; Số truyền 2: $V = 2.8\text{ km/h}$) đạt độ khớp nối rất cao. Phân tích thống kê hồi quy cho thấy hệ số tương quan đạt $R^2 = 0.97 \div 0.98$, khẳng định độ tin cậy và tính vững chắc của mô hình toán học đã thiết lập.

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

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

  1. Xác lập ngưỡng đối trọng tối thiểu triệt tiêu nguy cơ lật và bốc đầu: Khi LHM di chuyển mang đầy tải định mức ($Q_{VL} = 450\text{ kg}$), nếu không bố trí đối trọng ($Q_{DT} = 0$), phản lực pháp tuyến cầu sau giảm xuống xấp xỉ bằng không khi xe leo mấp mô, gây mất hoàn toàn lực bám đẩy và mất ổn định hướng chuyển động. Nghiên cứu chứng minh rằng bắt buộc phải trang bị khối lượng đối trọng $Q_{DT} \ge 450\text{ kg}$ phía sau đuôi máy kéo để phân bố lại tỷ lệ tải trọng bám cầu sau đạt trên 40% trọng lượng toàn bộ LHM.
Cấu hình đối trọng ($Q_{DT}$) Mức tải gầu ($Q_{VL}$) Trạng thái phản lực cầu sau ($P_{zk}$) Nguy cơ động lực học Khả năng điều khiển
$Q_{DT} = 0\text{ kg}$ $450\text{ kg}$ (Đầy tải) $P_{zk} \to 0\text{ kN}$ khi gặp mấp mô Bốc đầu, trượt quay bánh chủ động Mất lái hoàn toàn
$Q_{DT} = 250\text{ kg}$ $450\text{ kg}$ (Đầy tải) $P_{zk}$ dao động mạnh, biên độ lớn Tách bánh tức thời khi dao động Lái rất nặng, mất ổn định hướng
$Q_{DT} = 450\text{ kg}$ (Chuẩn) $450\text{ kg}$ (Đầy tải) $P_{zk} > 0.4 \cdot G_{total}$ (Ổn định) Duy trì tiếp xúc liên tục, chống lật Ổn định hướng tốt, đảm bảo lực bám
  1. Xác định vùng góc nâng an toàn động lực học: Nghiên cứu phát hiện góc nâng cần hợp lý khi mang tải di chuyển là $\alpha \in [-40^\circ, -30^\circ]$ (tương ứng với độ cao đáy gầu cách mặt đất $200 \div 400\text{ mm}$). Việc nâng cần lên vị trí cao hơn ($\alpha > 0^\circ$) làm tọa độ trọng tâm $h$ tăng vọt, gây khuếch đại gia tốc lắc ngang $\ddot{\theta}$ lên gấp $2.3$ lần khi đi qua mấp mô đơn và làm giảm góc nghiêng tới hạn chống lật ngang $\alpha_{\lim}$ từ $38^\circ$ xuống dưới $18^\circ$.

  2. Giới hạn vận tốc an toàn theo biên độ mấp mô: Khi di chuyển đầy tải ($Q_{VL} = 450\text{ kg}$), vận tốc di chuyển an toàn bắt buộc phải duy trì ở ngưỡng $V < 3\text{ km/h}$. Nếu vận tốc vượt qua ngưỡng này, độ chênh lệch biên độ sóng mấp mô hai vệt bánh $\Delta h_0 = |h_{0P} - h_{0T}| > 0.07\text{ m}$ sẽ kích phát hiện tượng cộng hưởng dao động góc lắc ngang thân xe, làm bánh xe bên cao tách hoàn toàn khỏi mặt đất ($P_{zT} = 0$ hoặc $P_{zP} = 0$).

  3. Quá tải hệ thống lái cơ khí và yêu cầu trợ lực thủy lực: Kết quả đo đạc và tính toán chỉ ra rằng, do tải trọng dồn lên cầu trước khi đầy tải tăng hơn 65% so với nguyên bản máy kéo Yanmar 3000, mô men cản quay vòng tại cầu trước $M_{cq}$ tăng mạnh, dẫn đến lực tác động lên vô lăng $F_{VL}$ vượt xa giới hạn cho phép theo Tiêu chuẩn Việt Nam (TCVN) về công thái học vận hành máy nông nghiệp. Do đó, việc bổ sung cơ cấu trợ lực lái thủy lực là yêu cầu kỹ thuật bắt buộc để đảm bảo an toàn lao động.

Implications đa chiều

  • Về mặt lý thuyết: Bổ sung phương pháp luận hoàn chỉnh trong việc mô hình hóa cơ hệ máy nông nghiệp đa năng có liên kết phụ tải biến thiên, làm phong phú giáo trình động lực học ô tô – máy kéo chuyên sâu.
  • Về thiết kế chế tạo cơ khí: Cung cấp bộ thông số kích thước, trọng lượng lắp ghép và sơ đồ bố trí tối ưu của bộ xúc lật Y3510FLH trên khung gầm máy kéo cỡ 30 HP, mở đường cho các doanh nghiệp cơ khí nội địa chế tạo module xúc lật lắp lẫn tiêu chuẩn hóa với chi phí thấp.
  • Về quy chuẩn vận hành và chính sách: Làm căn cứ khoa học để xây dựng sổ tay hướng dẫn an toàn kỹ thuật quốc gia cho máy nông nghiệp đa năng, ngăn ngừa tai nạn lật xe chết người trong canh tác và xây dựng nông thôn mới.

Limitations và Future Research

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

  1. Giả thiết khung xe tuyệt đối cứng: Mô hình toán chưa tính đến độ xoắn mềm cục bộ của sắt-xi khung gầm máy kéo khi chịu tải trọng xoắn lệch cực đại giữa cầu trước và cầu sau.
  2. Mô hình lốp tuyến tính hóa: Đặc trưng lốp xe ($C_i, k_i$) được mô hình hóa theo dạng tuyến tính đơn giản hóa, chưa phản ánh hết tính chất phi tuyến và hiện tượng trễ đàn hồi (hysteresis) của lốp áp suất thấp trên nền đất bùn lầy thụt sâu.
  3. Tính chất vật liệu tải cố định: Nghiên cứu giả định khối lượng vật liệu trong gầu xúc ($m_{VL}$) phân bố tĩnh tuyệt đối, chưa xét đến hiệu ứng dịch chuyển trọng tâm cục bộ của chất lỏng hoặc bùn thải lỏng sánh khi xe lắc mạnh.

Chương trình nghiên cứu tiếp theo (Future Research Agenda) cần mở rộng vào các hướng:

  • Nghiên cứu hệ thống điều khiển chủ động góc nâng cần thủy lực (Active Boom Suspension) nhằm tự động dập tắt dao động lắc dọc khi xe chạy tốc độ cao.
  • Thiết lập mô hình tương tác đất – lốp phi tuyến 3D bằng phương pháp phần tử hữu hạn phiếm định (DEM-FEM) trên nền đất yếu lầy thụt.
  • Tích hợp hệ thống lái tự động thông minh (Autonomous Navigation) dựa trên thuật toán kiểm soát ổn định chống lật theo thời gian thực.

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

Nghiên cứu mở ra tiềm năng trích dẫn học thuật cao trong các tạp chí chuyên ngành Cơ khí Động lực, Kỹ thuật Nông nghiệp (Agricultural Engineering) và Cơ học Hệ nhiều vật. Về mặt kinh tế – xã hội, việc tận dụng hơn 400.000 máy kéo hiện có để biến thành máy xúc lật đa năng giúp tiết kiệm hàng nghìn tỷ đồng chi phí đầu tư mua sắm máy chuyên dùng nhập khẩu từ nước ngoài.

Giải pháp nâng cao hiệu suất sử dụng máy kéo từ 2–3 tháng lên 8–10 tháng/năm, trực tiếp gia tăng thu nhập cho nông dân và các hợp tác xã nông nghiệp, đồng thời thúc đẩy cơ giới hóa khâu thu gom chất thải chăn nuôi, phụ phẩm sinh khối và xây dựng hạ tầng nông thôn bền vững.

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

  • Nghiên cứu sinh và Giảng viên ngành Cơ khí: Tiếp cận một mô hình toán học giải tích chặt chẽ, có khả năng kế thừa cho các loại liên hợp máy nông nghiệp khác (máy kéo – máy khoan lỗ, máy băm phụ phẩm, máy đào rãnh).
  • Kỹ sư R&D tại các nhà máy chế tạo máy nông nghiệp: Sử dụng trực tiếp bộ bản vẽ tính toán, sơ đồ phân bố tải trọng và kết cấu gá lắp module xúc lật Y3510FLH để thương mại hóa sản phẩm.
  • Hợp tác xã nông nghiệp và Chủ trang trại: Được trang bị quy trình vận hành an toàn chuẩn mực (tải trọng, tốc độ, góc nâng), loại trừ rủi ro lật xe và hư hỏng hệ thống truyền lực.
  • Cơ quan quản lý và Đăng kiểm máy nông nghiệp: Có cơ sở dữ liệu thực nghiệm chuẩn xác để ban hành các tiêu chuẩn an toàn kỹ thuật đối với máy kéo cải tiến, chuyển đổi công năng.

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ì? Đó là việc mở rộng hệ phương trình vi phân chuyển động phi tuyến Lagrange loại hai cho cơ hệ máy kéo liên hợp có ten-xơ quán tính ($J_x, J_y, J_z$) và tọa độ trọng tâm thay đổi liên tục theo hàm của góc nâng cần $\alpha$ và hệ số đầy tải $k_Q$, đồng thời tích hợp chính xác động học của khớp lắc thăng bằng cầu trước ($O_1x_1$).

2. Điểm cải tiến phương pháp luận so với các nghiên cứu quốc tế là gì? So với mô hình của Debeleac (2007) (vốn áp dụng cho khung cứng máy xúc chuyên dụng) và Zhen Li (2015) (mô hình máy kéo không có tải phía trước), luận án đã đề xuất phương pháp kích thích mấp mô tuần hoàn lệch pha giải tích ($S_P = S_T, d$) và kỹ thuật thực nghiệm mấp mô cung tròn ảo tiếp xúc điểm độc lập với sự chênh lệch đường kính bánh trước – sau.

3. Phát hiện nào gây bất ngờ nhất về mặt số liệu? Mặc dù vận tốc di chuyển rất thấp ($V < 3\text{ km/h}$), nhưng nếu độ lệch biên độ mấp mô hai bên bánh chỉ cần vượt quá $0.07\text{ m}$ khi xe mang đầy tải không có đối trọng, phản lực pháp tuyến cầu sau lập tức triệt tiêu hoàn toàn, biến máy kéo thành cơ hệ mất kiểm soát hướng nghiêm trọng.

4. Giao thức lặp lại thực nghiệm (Replication Protocol) có được cung cấp rõ ràng không? Luận án mô tả chi tiết toàn bộ cấu hình thử nghiệm: máy kéo Yanmar 3000 (động cơ 3T84), bộ xúc Y3510FLH, vị trí lắp đặt cảm biến laser, sơ đồ khối phần mềm DasyLab 10 và các bước lấy mẫu tín hiệu với độ trễ thấp, đảm bảo khả năng tái lập thí nghiệm tuyệt đối.

5. Định hướng nghiên cứu 10 năm tới là gì? Phát triển hệ thống cân bằng điện tử tích cực (Electronic Stability Program - ESP) cho máy nông nghiệp đa năng và ứng dụng công nghệ điều khiển thủy lực thích nghi (Adaptive Electro-Hydraulic System) chống rung lắc gầu xúc tự động.

Kết luận

Luận án của NCS. Phạm Trọng Phước là một công trình nghiên cứu công phu, giải quyết trọn vẹn bài toán khoa học và kỹ thuật thực tiễn về ổn định chuyển động của liên hợp máy kéo – xúc lật:

  1. Xây dựng thành công mô hình động lực học không gian giải tích hoàn chỉnh có tính đến sự biến thiên liên tục của các đặc trưng khối lượng và khớp lắc thăng bằng cầu trước.
  2. Đề xuất phương pháp kích thích mấp mô sin không gian lệch pha và mô hình thực nghiệm mấp mô đơn cung tròn ảo tiếp xúc điểm có giá trị học thuật cao.
  3. Xác lập chính xác miền giá trị thông số vận hành an toàn tối ưu khi đầy tải: Đối trọng $Q_{DT} \ge 450\text{ kg}$; Góc nâng cần $\alpha \in [-40^\circ, -30^\circ]$; Vận tốc $V < 3\text{ km/h}$; Độ chênh mấp mô $\Delta h_0 \le 0.07\text{ m}$.
  4. Chứng minh thực nghiệm yêu cầu bắt buộc phải lắp đặt hệ thống trợ lực lái thủy lực nhằm đảm bảo điều kiện công thái học và an toàn điều khiển.
  5. Kiểm chứng mô hình với độ tin cậy cơ học và thống kê vượt trội ($R^2 = 0.97 \div 0.98$).
  6. Mở ra hướng chuyển giao công nghệ module hóa bộ xúc lật Y3510FLH trên dải máy kéo 30–40 HP, mang lại hiệu quả kinh tế – kỹ thuật và an sinh xã hội sâu rộng cho nền nông nghiệp Việt Nam.