Tổng quan nghiên cứu

Trong quá trình công nghiệp hóa và hiện đại hóa ngành lâm nghiệp Việt Nam, việc cơ giới hóa các khâu thu hoạch và vận chuyển gỗ đóng vai trò then chốt nhằm giải phóng sức lao động thủ công và nâng cao năng suất. Số liệu thống kê cho thấy thị trường nội địa tiêu thụ tới 98% sản lượng gỗ tròn, 92% gỗ xẻ và khoảng 80% nguyên liệu giấy. Tuy nhiên, khâu vận xuất và vận chuyển gỗ trung gian trong cự ly từ 10 km đến 15 km trên địa hình đồi núi phức tạp vẫn gặp nhiều hạn chế về mặt kỹ thuật, làm gia tăng tổn thất và ảnh hưởng sức khỏe công nhân vận hành.

Vấn đề cốt lõi đặt ra là sự rung xóc mạnh của các tổ hợp máy kéo nông lâm nghiệp khi làm việc trên các tuyến đường rừng gồ ghề. Đề tài tập trung nghiên cứu đối tượng liên hợp máy kéo MTZ-82 công suất 80 mã lực kết hợp rơ moóc một trục chở gỗ có tải trọng thiết kế 40 kN (tương đương 4 tấn). Mục tiêu cụ thể của nghiên cứu là thiết lập mô hình toán học động lực học không gian trong mặt phẳng đối xứng dọc, phân tích các chỉ tiêu dao động thẳng đứng và dao động góc, từ đó xác định nguyên nhân gây mất êm dịu và đề xuất các giải pháp kỹ thuật cải tiến hệ thống treo và phân bố tải trọng.

Phạm vi nghiên cứu được thực hiện tại các vùng rừng trồng trọng điểm phía Bắc trong giai đoạn hoàn thiện sản phẩm cơ khí nông lâm nghiệp năm 2010. Kết quả nghiên cứu mang ý nghĩa thực tiễn to lớn khi giúp giảm biên độ gia tốc dao động của ghế lái xuống dưới ngưỡng 1,8 m/s², cải thiện chỉ số độ êm dịu K trong dải an toàn từ 10 đến 25 đơn vị, nâng cao tuổi thọ kết cấu khung gầm và đảm bảo an toàn lao động cho người vận hành máy kéo lâm nghiệp.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu ứng dụng lý thuyết động lực học ô tô máy kéo hiện đại và cơ học hệ nhiều vật để khảo sát tương tác động lực giữa xe, rơ moóc và mặt đường. Ba lý thuyết và mô hình nền tảng được tích hợp bao gồm: lý thuyết dao động cơ học phi tuyến áp dụng phương trình vi phân Lagrange loại II, lý thuyết tiếp xúc đàn hồi giữa bánh lốp và nền đường đất mấp mô, cùng tiêu chuẩn đánh giá độ êm dịu chuyển động dựa trên hệ số hấp thụ năng lượng dao động K của cơ thể người.

Các khái niệm then chốt được lượng hóa bao gồm:

  • Tần số dao động riêng: Số chu kỳ dao động tự do của liên hợp máy trong một đơn vị thời gian, phụ thuộc vào độ cứng quy đổi 5.770 N/m của cụm ghế ngồi và độ cứng lốp sau 432.000 N/m.
  • Hệ số độ êm dịu chuyển động K: Chỉ tiêu thể hiện cảm giác rung xóc của người lái dựa trên gia tốc hiệu dụng bình phương trung bình trong dải tần số 1 Hz đến 5 Hz.
  • Hàm kích động mặt đường xác định: Quy luật biên dạng mấp mô sóng đường tuần hoàn với biên độ sóng từ 1,8 cm đến 4,0 cm và bước sóng biến thiên từ 0,6 m đến 1,0 m.

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

Nghiên cứu sử dụng kết hợp phương pháp mô hình hóa giải tích, mô phỏng số trên máy tính và đo kiểm thực nghiệm trực tiếp. Bộ dữ liệu thông số đầu vào bao gồm 16 đại lượng hình học và động học của máy kéo MTZ-82 (khối lượng bản thân 3.450 kg, mô men quán tính 1.973 kg.m²) và rơ moóc chở gỗ (khối lượng quy đổi 800 kg tại điểm nối moóc, mô men quán tính 1.567 kg.m²).

Cỡ mẫu khảo sát thực nghiệm mặt đường được chọn ngẫu nhiên qua 4 trắc diện dọc đại diện cho các tuyến đường lâm nghiệp chuyên dụng tại vùng đồi núi phía Bắc. Phương pháp lấy mẫu chọn điểm đo tập trung tại các vị trí mấp mô gợn sóng đặc thù với chu kỳ lặp lại liên tục.

Lý do lựa chọn phương pháp lập phương trình Lagrange loại II và mô phỏng trên nền tảng phần mềm Matlab Simulink là khả năng giải quyết chính xác hệ 4 phương trình vi phân bậc hai liên kết động tuyến tính hóa với 4 tọa độ suy rộng: chuyển vị thẳng đứng trọng tâm máy kéo, góc lắc dọc thân xe, chuyển vị thẳng đứng của ghế lái và góc lắc dọc thân rơ moóc. Phương pháp thực nghiệm đo độ cứng và hệ số cản của lò xo ghế lái sử dụng tải trọng nén 20 kg (gây biến dạng 0,034 m) cùng hệ thống đếm chu kỳ dao động tắt dần đạt chu kỳ trung bình 1,886 giây, bảo đảm sai số đo lường dưới 5% so với chuẩn kỹ thuật.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Quá trình mô phỏng động lực học trên miền thời gian 20 giây với vận tốc di chuyển 8 km/h đến 12 km/h trên biên dạng đường gợn sóng biên độ 0,04 m đã mang lại các phát hiện kỹ thuật quan trọng:

  • Biên độ gia tốc thẳng đứng: Trọng tâm máy kéo dao động tuần hoàn duy trì biên độ gia tốc trong khoảng từ -1,8 m/s² đến +1,8 m/s². Khi chưa có giảm chấn phụ, gia tốc đỉnh tại sàn buồng lái có thể vượt quá 2,5 m/s², gây mệt mỏi nhanh cho người điều khiển.
  • Tần số dao động cưỡng bức: Đồ thị thời gian ghi nhận 4 chu kỳ dao động hoàn chỉnh trong khoảng thời gian 2,0 giây, tương ứng tần số kích động là 2,0 Hz. Tần số này nằm sát vùng nhạy cảm cộng hưởng sinh học của nội tạng cơ thể người (dao động từ 1,5 Hz đến 3,0 Hz).
  • Chuyển vị góc thân xe và rơ moóc: Góc lắc dọc của máy kéo biến thiên trong phạm vi nhỏ từ -0,008 rad đến +0,008 rad, trong khi góc xoay của rơ moóc quanh điểm nối khớp đạt tới biên độ 0,1 rad (tương đương lệch góc khoảng 5,7 độ). Chuyển vị góc lớn này làm tăng 25% lực uốn động tác dụng lên khớp nối moóc.
  • Hiệu quả giảm xóc ghế ngồi: Độ cứng lò xo ghế ở mức 5.770 N/m giúp giảm độ biến thiên chuyển dịch thẳng đứng của người lái xuống còn 0,012 m, hạn chế 40% biên độ xung kích so với dao động truyền trực tiếp từ sàn máy kéo.

Thảo luận kết quả

Các kết quả mô phỏng động lực học được tổng hợp và hiển thị trực quan thông qua đồ thị chuyển vị dạng sóng tuần hoàn theo thời gian (Hình mô phỏng từ 0 đến 20 giây) và bảng ma trận các hệ số khối lượng, độ cứng, lực cản cản dịu. Dữ liệu biểu đồ gia tốc khẳng định dao động của hệ thống là dao động duy trì không tắt khi nguồn kích thích mấp mô mặt đường còn tiếp diễn.

Nguyên nhân chính dẫn đến biên độ lắc dọc lớn của rơ moóc là do liên hợp nguyên bản không có hệ thống treo trung gian; toàn bộ biến dạng đàn hồi chỉ phụ thuộc vào lốp xe với độ cứng cao lên tới 679.000 N/m. So sánh với các nghiên cứu dao động của xe tải lâm nghiệp chuyên dụng ở Bắc Âu, việc thiếu cụm nhíp đàn hồi khiến phản lực động từ mặt đường truyền thẳng lên trục rơ moóc tăng thêm 35%.

Việc phân tích dữ liệu trên đồ thị phổ tần số và bảng so sánh các chế độ vận tốc cho thấy khi máy kéo di chuyển ở dải tốc độ trên 15 km/h qua đoạn mấp mô bước sóng 1,0 m, hiện tượng va đập cục bộ tăng vọt 50%, làm hệ số êm dịu K vượt mức 30 đơn vị. Điều này khẳng định sự cần thiết phải can thiệp đồng bộ vào kết cấu nhíp giảm chấn và quy trình vận hành tốc độ.

Đề xuất và khuyến nghị

Dựa trên kết quả mô phỏng và thực nghiệm động lực học, 4 nhóm giải pháp kỹ thuật và vận hành được đề xuất:

  • Lắp đặt hệ thống treo phụ thuộc loại nhíp lá cho rơ moóc: Trang bị bộ nhíp lá bán elip kết hợp giảm chấn ống thủy lực tác dụng hai chiều giữa cầu xe và khung rơ moóc. Mục tiêu giảm 35% gia tốc dao động truyền lên khung rơ moóc, thực hiện trong vòng 6 tháng bởi các xưởng cơ khí lâm nghiệp địa phương.
  • Tối ưu hóa thông số bộ giảm chấn ghế lái: Căn chỉnh độ cứng lò xo ghế ở mức 5.500 N/m đến 5.800 N/m và hệ số cản giảm chấn đạt 170 Ns/m đến 180 Ns/m nhằm giữ tần số dao động riêng của ghế trong dải tối ưu 1,0 Hz đến 1,2 Hz. Trách nhiệm thuộc về đội ngũ kỹ thuật viên bảo dưỡng máy trước mỗi chu kỳ khai thác gỗ 3 tháng.
  • Quy chuẩn hóa dải vận tốc vận chuyển trên đường lâm nghiệp: Ban hành quy định giới hạn dải tốc độ làm việc an toàn từ 8 km/h đến 12 km/h khi kéo tải đầy 40 kN trên các cung đường đất đồi gợn sóng có chiều cao mấp mô trên 2,0 cm, giúp giảm 45% nguy cơ lật và trượt lết moóc. Giải pháp áp dụng ngay cho công nhân vận hành lâm trường.
  • Cải tiến vị trí tâm quay khớp nối móc kéo: Hạ thấp trọng tâm điểm nối moóc cách mặt đất 0,45 m và tăng chiều dài đòn nối moóc lên 1,3 m để giảm 20% mô men kích động lắc dọc truyền ngược lại cầu sau máy kéo MTZ-82. Thực hiện thiết kế chuẩn hóa bởi Viện Nghiên cứu Cơ điện Nông nghiệp trong lộ trình 12 tháng.

Đối tượng nên tham khảo luận văn

  • Kỹ sư thiết kế máy nông lâm nghiệp: Nắm vững phương pháp thiết lập hệ phương trình vi phân dao động nhiều bậc tự do, ma trận độ cứng và ma trận cản để ứng dụng tính toán thiết kế các tổ hợp máy kéo rơ moóc mới.
  • Cán bộ quản lý kỹ thuật tại các lâm trường và công ty chế biến gỗ: Ứng dụng quy trình đánh giá độ êm dịu và điều chỉnh áp suất lốp, tải trọng chuyên chở thực tế nhằm tăng 30% độ bền phương tiện và giảm thiểu chi phí sửa chữa khung gầm.
  • Học viên cao học và sinh viên chuyên ngành Cơ giới hóa, Kỹ thuật ô tô máy kéo: Sử dụng làm tài liệu tham khảo chuẩn mực về mô hình hóa động lực học trên Matlab Simulink và phương pháp thực nghiệm đo lường đại lượng không điện bằng thiết bị số.
  • Chuyên gia an toàn lao động và công thái học lâm nghiệp: Khai thác dữ liệu về phổ tần số rung động và hệ số độ êm dịu K để xây dựng định mức thời gian lái máy an toàn, ngăn ngừa các bệnh nghề nghiệp về cột sống cho công nhân lái máy kéo.

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

Tại sao dao động của máy kéo MTZ-82 khi kéo rơ moóc lại nguy hiểm hơn khi chạy đơn lẻ?
Khi kéo rơ moóc chở đầy 40 kN gỗ, trọng lượng bổ sung tác dụng lên điểm nối moóc tạo ra mô men kích động lắc dọc biến thiên liên tục. Phản lực này làm thay đổi tải trọng bám của cầu trước máy kéo tới 30%, gây mất ổn định hướng lái và gia tăng biên độ lắc dọc thân xe.

Phương pháp Lagrange loại II có ưu điểm gì trong bài toán dao động máy kéo?
Phương pháp Lagrange loại II cho phép biểu diễn tường minh động năng, thế năng và hàm hao tán thông qua 4 tọa độ suy rộng độc lập. Cách tiếp cận này loại bỏ hoàn toàn các phản lực liên kết phức tạp bên trong cơ hệ, giúp việc thành lập hệ phương trình vi phân ma trận trở nên chính xác và thuận lợi cho giải số.

Hệ số độ êm dịu K bao nhiêu là đạt chuẩn đối với người vận hành máy kéo lâm nghiệp?
Theo tiêu chuẩn công thái học ô tô máy kéo, hệ số K trong khoảng 10 đến 25 đơn vị là ngưỡng an toàn cho người vận hành làm việc liên tục trong ca 8 giờ. Khi hệ số K vượt quá 30, cơ thể người lái sẽ nhanh chóng xuất hiện cảm giác mệt mỏi, giảm phản xạ điều khiển và tăng nguy cơ chấn thương cột sống.

Độ cứng của lốp xe ảnh hưởng như thế nào đến độ êm dịu của liên hợp máy?
Lốp xe đóng vai trò là phần tử đàn hồi sơ cấp nối trực tiếp với mặt đường. Lốp sau máy kéo có độ cứng 432.000 N/m và lốp rơ moóc có độ cứng 679.000 N/m; nếu bơm lốp quá căng, tần số dao động riêng tăng vọt, làm tăng 40% gia tốc chấn động truyền trực tiếp lên khung xe.

Làm cách nào để ứng dụng kết quả nghiên cứu vào thực tế sản xuất tại lâm trường?
Lâm trường có thể ứng dụng ngay việc giới hạn tốc độ di chuyển từ 8 km/h đến 12 km/h khi chở nặng, đồng thời lắp bổ sung hệ thống nhíp lá chuyên dụng cho các moóc chở gỗ một trục hiện có và định kỳ cân chỉnh lò xo ghế lái theo thông số độ cứng 5.770 N/m.

Kết luận

  • Xây dựng thành công mô hình dao động không gian 4 bậc tự do trong mặt phẳng đối xứng dọc cho liên hợp máy kéo MTZ-82 kết hợp rơ moóc chở gỗ 40 kN.
  • Thiết lập hệ phương trình vi phân động lực học dạng ma trận bằng phương pháp Lagrange loại II và giải tự động trên nền tảng phần mềm Matlab Simulink.
  • Xác định chính xác các thông số dao động then chốt với tần số kích động 2,0 Hz, biên độ gia tốc trọng tâm máy kéo dao động trong dải an toàn -1,8 m/s² đến +1,8 m/s².
  • Đóng góp giải pháp lắp đặt hệ thống treo nhíp lá đàn hồi cho rơ moóc và tối ưu hóa hệ số cản lò xo ghế lái ở mức 173 Ns/m.
  • Đề ra lộ trình triển khai thử nghiệm diện rộng trong 12 tháng tới nhằm hoàn thiện quy chuẩn kỹ thuật khai thác gỗ rừng trồng bền vững.

Quý độc giả, kỹ sư và nhà nghiên cứu quan tâm đến mô hình tính toán động lực học chi tiết cùng tệp mã nguồn mô phỏng Matlab Simulink có thể liên hệ trực tiếp bộ môn Cơ giới hóa để cùng hợp tác ứng dụng và phát triển đề tài vào thực tiễn sản xuất.