Tổng quan nghiên cứu

Nguồn tài nguyên rừng Việt Nam đóng vai trò then chốt trong việc bảo vệ môi trường sinh thái và cung cấp nguyên liệu cho các ngành kinh tế. Theo thống kê điều tra, diện tích đất lâm nghiệp cả nước đạt khoảng 12,699 triệu ha, chiếm tỷ lệ 38,2% tổng diện tích tự nhiên. Trong cơ cấu rừng, rừng tự nhiên chiếm 81,5% tương đương 10,2 triệu ha, còn rừng trồng chiếm 18,5% tương ứng 2,3 triệu ha. Nhằm định hướng phát triển lâm nghiệp bền vững, Chiến lược phát triển lâm nghiệp Việt Nam giai đoạn 2006 – 2020 đã đặt mục tiêu trồng mới 1,0 triệu ha rừng đến năm 2010 và bổ sung 1,5 triệu ha trong giai đoạn kế tiếp, đáp ứng nhu cầu nguyên liệu cho ngành công nghiệp giấy, ván dăm và chế biến gỗ xuất khẩu.

Trước thập niên 1990, hoạt động khai thác rừng tự nhiên chủ yếu sử dụng các thiết bị chuyên dùng công suất lớn như máy kéo bánh xích TDT-75, TT-4 hoặc máy kéo bánh hơi LKT-80. Tuy nhiên, đặc thù rừng trồng hiện nay có trữ lượng cây nhỏ, vị trí phân tán và địa hình đồi dốc phức tạp với độ dốc dọc đường vận xuất cho phép lên đến 25%. Việc đầu tư máy kéo chuyên dùng nhập ngoại đắt tiền không còn phù hợp với khả năng tài chính của các nông lâm trường. Hướng đi tối ưu là cải tiến các dòng máy kéo nông nghiệp 4 bánh có công suất từ 35 mã lực trở lên phục vụ vận xuất gỗ theo phương pháp kéo nửa lết.

Khi vận hành trên đường vận xuất lâm nghiệp mấp mô, máy kéo phải chịu tác động kép từ kích động mặt đường và lực vẫy giật từ đuôi bó gỗ. Dao động ngang xuất hiện làm suy giảm nghiêm trọng khả năng bám, độ êm dịu và tính ổn định chống lật của phương tiện. Luận văn tập trung nghiên cứu ảnh hưởng của dao động ngang đến khả năng ổn định của máy kéo Shibaura-3000A, thiết lập hệ phương trình vi phân động lực học và xác lập chế độ vận hành tối ưu nhằm đảm bảo an toàn kỹ thuật cho liên hợp máy.

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

Khung lý thuyết áp dụng

Luận văn ứng dụng lý thuyết dao động ô tô - máy kéo kết hợp phương trình vi phân cơ học Lagrange loại 2 để mô tả chuyển động của cơ hệ trong không gian 3 chiều. Hệ động lực học liên hợp bao gồm 4 khối lượng quán tính cơ bản: thân máy kéo kèm thiết bị tời cáp, cặp bánh trước, cặp bánh sau và khối gỗ vận xuất. Khung lý thuyết tích hợp ba khái niệm nền tảng: quá trình kích động ngẫu nhiên dừng và ecgôđich của mấp mô mặt đường, ma trận độ cứng đàn hồi của hệ thống lốp và ma trận hệ số cản hao tán năng lượng.

Mô hình khảo sát 4 bậc tự do suy rộng gồm: chuyển vị thẳng đứng của trọng tâm máy kéo Z, góc lắc dọc thân máy alpha quanh trục ngang, góc lắc ngang của cầu trước beta1 và góc lắc ngang của cầu sau beta2 quanh trục dọc chuyển động. Mặt đường vận xuất được mô hình hóa dưới dạng hàm điều hòa có bước sóng S0 = 1,0 m và biên độ mấp mô h0 biến thiên từ 0 m đến 1,0 m. Khung lý thuyết đã bổ sung thành phần dịch chuyển ngang của trọng tâm y0 = (h0 + z0)*sin(beta) với chiều cao trọng tâm 552 mm, giúp mô phỏng chính xác hơn hiện tượng mất ổn định ngang so với các mô hình dao động phẳng truyền thống.

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

Nghiên cứu lựa chọn đối tượng điển hình là máy kéo nông nghiệp Shibaura-3000A trang bị tời cáp kéo gỗ, động cơ diesel 4 kỳ công suất 35 mã lực (tương đương 26 kW). Thông số hình học và tải trọng của máy gồm: tổng khối lượng bản thân 1477 kg, phân bố tải tĩnh lên cầu trước là 650 kg (phản lực Z1 = 6380 N) và cầu sau là 827 kg (phản lực Z2 = 8112 N), chiều dài cơ sở 1815 mm, khoảng cách vệt bánh trước và sau đều bằng 1200 mm.

Dữ liệu mô men quán tính được xác định thông qua phương pháp thực nghiệm con lắc vật lý và treo đàn hồi, ghi nhận mô men quán tính của thân máy với trục ngang Jx = 305 kgm2 và trục dọc Jy = 911 kgm2, mô men quán tính cầu trước Jx1 = 42 kgm2. Thông số đàn hồi lốp xe được đo đạc thực nghiệm: cặp lốp trước 8-16 có độ cứng C1 = 358.808 N/m, hệ số cản K1 = 4.715 Ns/m; cặp lốp sau 11-28 có độ cứng C2 = 437.086 N/m, hệ số cản K2 = 4.290 Ns/m. Hệ phương trình vi phân ma trận phi thuần nhất bậc hai được giải bằng phương pháp giải tích trên chương trình máy tính ngôn ngữ Pascal, kết hợp thuật toán phân tích ma trận trực giao SVD để trích xuất tần số dao động riêng và xử lý biểu đồ trên Excel. Phương pháp giải tích kết hợp SVD được lựa chọn vì đảm bảo độ hội tụ cao và triệt tiêu sai số số học khi tính toán cơ hệ nhiều bậc tự do.

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

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

Thứ nhất, việc bổ sung thành phần dao động lắc ngang của trọng tâm máy kéo làm thay đổi căn bản phản lực động học so với mô hình phẳng truyền thống. Khi biên độ mấp mô mặt đường tăng từ h0 = 0 m lên h0 = 1,0 m, biên độ chuyển vị đứng Z và góc lắc dọc alpha tăng từ 28% đến 45%, đồng thời góc lắc ngang cầu trước beta1 và cầu sau beta2 tăng mạnh hơn 60%, làm gia tăng đáng kể nguy cơ mất cân bằng ngang khi xe chở nặng.

Thứ hai, kết quả khảo sát trên 7 mức tải trọng gỗ gồm 160 kg, 200 kg, 240 kg, 280 kg, 320 kg, 360 kg và 400 kg cho thấy khối lượng gỗ tỷ lệ thuận với mức độ kích động dao động. Ở vận tốc di chuyển v = 2,0 m/s (7,2 km/h), khi tải trọng gỗ tăng từ 160 kg lên 400 kg, chuyển vị đứng Z tăng từ 0,05 m lên đỉnh điểm 0,35 m, biên độ góc lắc cầu sau beta2 tăng thêm 32%, làm giảm khả năng bám đường của bánh sau.

Thứ ba, tốc độ di chuyển là tham số nhạy cảm nhất đối với độ ổn định của liên hợp máy. Ở chế độ vận tốc v = 1,2 m/s (4,32 km/h), các giá trị biên độ dao động góc duy trì dưới ngưỡng an toàn 0,08 rad trên toàn bộ dải tải trọng từ 160 kg đến 400 kg. Khi vận tốc tăng lên v = 2,0 m/s trên mặt đường có mấp mô h0 = 1,0 m, phản lực tiếp đất của bánh trước có thời điểm giảm về 0 N, gây hiện tượng mất khả năng dẫn hướng cục bộ.

Thứ tư, nghiên cứu xác định góc dốc giới hạn lật tĩnh của máy kéo bánh hơi khi đầy tải đạt từ 35 độ đến 40 độ khi lên dốc và trên 60 độ khi xuống dốc. Tuy nhiên, trong điều kiện kéo động học trên đường vận xuất có độ dốc 25%, tác động dao động ngang từ đuôi cây gỗ làm giảm góc dốc an toàn làm việc thực tế xuống còn khoảng 18 độ đến 22 độ.

Thảo luận kết quả

Các kết quả động lực học được mô tả trực quan qua các đồ thị hàm truyền dao động theo khối lượng gỗ và vận tốc di chuyển. Đồ thị phản ánh rõ nét sự phân hóa giữa hai chế độ vận tốc v = 1,2 m/s và v = 2,0 m/s: đường đặc tính biên độ dao động Z, alpha, beta1, beta2 dốc đứng hơn rõ rệt khi tốc độ vượt qua ngưỡng 1,8 m/s, cho thấy sự cộng hưởng năng lượng từ mấp mô mặt đường vào lốp sau có độ cứng cao 437.086 N/m.

So với các công trình trước đây chỉ khảo sát dao động thẳng đứng hoặc bỏ qua độ cao trọng tâm 552 mm, mô hình không gian trong luận văn đã giải thích thỏa đáng hiện tượng phân bố lại tải trọng không đều giữa hai vệt bánh xe trái và phải. Cầu trước máy kéo có khớp bản lề xoay giúp giảm biên độ lắc truyền lên khung máy, nhưng cầu sau liên kết cứng phải hấp thụ toàn bộ xung lực xoắn ngang. Hiện tượng lật máy kéo trong thực tế lâm nghiệp thường xảy ra ở cụm cầu sau khi xe vào cua hoặc đi qua mấp mô đơn bên ở tốc độ trên 5 km/h. Do đó, việc khống chế tốc độ kéo gỗ dưới 1,2 m/s trên đường vận xuất lâm nghiệp cấp 3 là hoàn toàn có cơ sở khoa học xác thực.

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

Dựa trên kết quả tính toán động lực học, luận văn đưa ra 4 khuyến nghị kỹ thuật và giải pháp tổ chức sản xuất cụ thể:

Thứ nhất, kiểm soát chặt chẽ chế độ vận tốc và tải trọng vận xuất. Người vận hành máy kéo Shibaura-3000A cần duy trì tốc độ làm việc không vượt quá 1,2 m/s (khoảng 4,3 km/h) khi kéo khối lượng gỗ từ 320 kg đến 400 kg. Trên các cung đường có độ dốc từ 15% đến 25%, bắt buộc phải giảm vận tốc xuống dưới 1,0 m/s và vận hành ở tầng số chậm để triệt tiêu dao động lắc ngang.

Thứ hai, chuẩn hóa chất lượng thi công và bảo dưỡng đường vận xuất gỗ. Các lâm trường cần tiến hành san ủi mặt đường theo đúng tiêu chuẩn đường cấp 3 lâm nghiệp, loại bỏ toàn bộ gốc cây, đá tảng và mô đất cao trên 0,2 m. Mặt đường vận xuất phụ cần giữ biên độ mấp mô h0 dưới 0,3 m, duy trì bề rộng nền đường tối thiểu từ 2,5 m đến 3,0 m và bố trí rãnh thoát nước mặt hợp lý.

Thứ ba, tối ưu hóa áp suất khí nén và bảo trì hệ thống lốp máy kéo. Định kỳ kiểm tra và duy trì áp suất lốp trước 8-16 ở mức 1,8 đến 2,2 kG/cm2 và lốp sau 11-28 ở mức 2,0 đến 2,5 kG/cm2. Áp suất tiêu chuẩn giúp ổn định độ cứng lốp trước ở mức 358.808 N/m và lốp sau ở mức 437.086 N/m, tối đa hóa hệ số dập tắt dao động K1 = 4.715 Ns/m và K2 = 4.290 Ns/m.

Thứ tư, cải tiến cơ cấu treo giữ gỗ và tời cáp phía sau máy kéo. Đơn vị cơ khí chế tạo nên lắp đặt bổ sung cụm khớp xoay giảm chấn thủy lực tại điểm liên kết cáp kéo nhằm hấp thụ phản lực vẫy ngang từ đuôi cây gỗ. Đồng thời, cần bố trí thêm đối trọng phụ từ 50 kg đến 80 kg tại đầu máy kéo để giữ phản lực tiếp đất của cầu trước luôn đạt trên 2000 N, đảm bảo khả năng dẫn hướng an toàn.

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

Luận văn là nguồn tư liệu chuyên môn hữu ích cho 4 nhóm đối tượng chính:

Thứ nhất, kỹ sư thiết kế và chế tạo máy nông lâm nghiệp: Tiếp cận mô hình động lực học không gian 4 bậc tự do và bộ thông số quán tính (Jx = 305 kgm2, Jy = 911 kgm2) để phục vụ tính toán cải tiến, chuyển đổi công năng máy kéo nông nghiệp thành máy vận xuất gỗ chuyên dùng.

Thứ hai, cán bộ quản lý kỹ thuật tại các công ty lâm nghiệp và lâm trường: Nắm vững các định mức tải trọng vận xuất an toàn từ 160 kg đến 400 kg, thiết lập quy trình vận hành phương tiện chuẩn hóa giúp giảm thiểu 20% nguy cơ lật đổ và tiết kiệm 15% nhiên liệu tiêu hao.

Thứ ba, giảng viên và nghiên cứu sinh chuyên ngành Cơ khí động lực: Sử dụng tài liệu như một công trình mẫu mực về phương pháp kết hợp giữa đo đạc thực nghiệm dao động đàn hồi và giải thuật toán ma trận SVD trên ngôn ngữ lập trình kỹ thuật.

Thứ tư, học viên cao học và sinh viên ngành Cơ giới hóa Nông – Lâm nghiệp: Bổ sung kiến thức chuyên sâu về lý thuyết dao động ô tô - máy kéo, phương pháp lập phương trình Lagrange loại 2 và kỹ năng phân tích tương tác động học giữa bánh lốp với nền đường tự nhiên.

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

Mô hình dao động có tính đến lắc ngang mang lại lợi thế gì so với mô hình phẳng truyền thống? Mô hình dao động không gian bổ sung chuyển vị lắc ngang quanh trục dọc với chiều cao trọng tâm 552 mm. Điều này cho phép xác định chính xác góc lắc cầu trước beta1 và cầu sau beta2 khi máy kéo chịu kích động mấp mô không đối xứng, phản ánh trung thực nguy cơ mất ổn định lật ngang mà mô hình phẳng 2 bậc tự do bỏ qua.

Máy kéo nông nghiệp Shibaura-3000A công suất 35 mã lực có phù hợp để vận xuất gỗ rừng trồng không? Máy kéo Shibaura-3000A có công suất 35 mã lực, trọng lượng 1477 kg cùng 4 cấp tốc độ trục thu công suất hoàn toàn đáp ứng tốt việc vận xuất gỗ rừng trồng theo phương pháp nửa lết. Máy vận hành hiệu quả với tải trọng gỗ từ 160 kg đến 400 kg trên các tuyến đường dốc dưới 25%.

Vì sao vận tốc kéo gỗ được khuyến nghị duy trì ở mức dưới 1,2 m/s? Kết quả mô phỏng cho thấy khi tốc độ tăng từ 1,2 m/s lên 2,0 m/s trên mấp mô h0 = 1,0 m, biên độ chuyển vị đứng Z tăng từ 0,05 m lên 0,35 m và góc lắc cầu sau tăng trên 30%. Tốc độ dưới 1,2 m/s (4,32 km/h) đảm bảo bánh xe luôn bám đường và triệt tiêu nguy cơ trượt lật.

Đặc tính đàn hồi của bánh lốp đóng vai trò như thế nào trong bài toán ổn định dao động? Bánh lốp là bộ phận đàn hồi duy nhất ở cầu sau máy kéo. Với độ cứng lốp trước C1 = 358.808 N/m và lốp sau C2 = 437.086 N/m cùng hệ số cản K1 = 4.715 Ns/m, K2 = 4.290 Ns/m, việc duy trì đúng áp suất lốp giúp dập tắt nhanh các dao động kích động từ mặt đường truyền lên thân máy.

Đường vận xuất lâm nghiệp cấp 3 cần đáp ứng những tiêu chuẩn kỹ thuật cốt lõi nào? Theo quy phạm kỹ thuật lâm nghiệp, đường vận xuất cấp 3 cho máy kéo bánh hơi phải có nền đường cứng, độ dốc dọc không vượt quá 25%, mặt đường được san ủi phẳng với biên độ mấp mô dưới 0,3 m và không có độ dốc ngang lớn để đảm bảo an toàn tuyệt đối khi kéo gỗ.

Kết luận

Công trình nghiên cứu đã giải quyết trọn vẹn bài toán động lực học về dao động ngang và độ ổn định của máy kéo nông nghiệp khi vận xuất gỗ rừng trồng với các đóng góp cốt lõi:

  • Xây dựng thành công mô hình động lực học không gian 4 bậc tự do có tính đến chuyển vị lắc ngang của trọng tâm máy kéo khi kéo gỗ nửa lết.
  • Xác định hệ thống thông số thực nghiệm chuẩn xác của máy kéo Shibaura-3000A với độ cứng lốp trước C1 = 358.808 N/m, lốp sau C2 = 437.086 N/m và mô men quán tính Jx = 305 kgm2, Jy = 911 kgm2.
  • Làm rõ quy luật ảnh hưởng của biên độ mấp mô mặt đường (0 m đến 1,0 m) và tốc độ di chuyển (1,2 m/s đến 2,0 m/s) đến chuyển vị đứng Z cùng các góc lắc dao động alpha, beta1, beta2.
  • Khẳng định tính khả thi và an toàn của máy kéo Shibaura-3000A khi vận chuyển tải trọng gỗ từ 160 kg đến 400 kg ở dải vận tốc dưới 1,2 m/s trên đường dốc dưới 25%.
  • Đề xuất đồng bộ 4 nhóm giải pháp kỹ thuật về vận hành, duy tu đường lâm nghiệp, hiệu chỉnh áp suất lốp và cải tiến cơ cấu tời cáp giảm chấn.

Hướng nghiên cứu tiếp theo cần tập trung phát triển mô hình dao động ngẫu nhiên phi tuyến trong điều kiện mặt đường trơn trượt mùa mưa giai đoạn 2026 – 2030. Các đơn vị khai thác lâm nghiệp và xưởng cơ khí chế tạo hãy áp dụng ngay các thông số tối ưu từ luận văn để nâng cao năng suất khai thác và bảo đảm an toàn lao động trong sản xuất lâm nghiệp.