Tổng quan nghiên cứu

Quá trình công nghiệp hóa và hiện đại hóa nông nghiệp, nông thôn tại Việt Nam đặt ra yêu cầu cấp thiết về cơ giới hóa đồng bộ các khâu canh tác. Thống kê toàn quốc cho thấy tổng nguồn động lực cơ giới phục vụ sản xuất nông, lâm nghiệp và thủy lợi đạt trên 16,3 triệu mã lực, tương đương mức trang bị bình quân khoảng 0,8 mã lực/ha. Tuy nhiên, mức độ cơ giới hóa phân bố không đồng đều giữa các vùng sinh thái: khu vực Tây Nguyên đạt 1,78 mã lực/ha, vùng Đông Nam Bộ đạt 1,25 mã lực/ha, trong khi khu vực Trung du và miền núi phía Bắc chỉ đạt 0,36 mã lực/ha. Tỷ lệ cơ giới hóa khâu làm đất bình quân cả nước đạt 63,8%, trong đó vùng Đồng bằng sông Cửu Long đạt 87%, nhưng tại các vùng đồi núi chỉ đạt xấp xỉ 41%.

Thực trạng canh tác tại vùng đồi núi có địa hình dốc, chia cắt phức tạp đòi hỏi các dòng máy kéo nhỏ hai bánh cầm tay công suất 8 đến 15 mã lực phải đảm nhận khối lượng công việc lớn. Dòng máy kéo Bông Sen 12 (BS 12) với công suất 12 CV do Việt Nam chế tạo được ứng dụng rộng rãi nhờ tính cơ động và tiêu hao nhiên liệu trung bình 19 lít/ha. Dù vậy, các thiết kế hiện hành chủ yếu dựa trên mô hình tĩnh học nền phẳng, chưa phản ánh đầy đủ tương tác động lực học khi máy kéo liên hợp với cày trụ một lưỡi làm việc trên sườn dốc.

Mục tiêu của nghiên cứu là xây dựng mô hình lý thuyết và giải thuật động lực học nhằm đánh giá khả năng kéo, bám, xác định ranh giới ổn định tĩnh và ổn định động của liên hợp máy kéo Bông Sen 12 khi canh tác trên sườn dốc nghiêng từ -30 độ đến +30 độ. Nghiên cứu được thực hiện tại Trường Đại học Lâm nghiệp, cung cấp cơ sở khoa học giúp tối ưu kết cấu, nâng cao năng suất cày đạt 0,13 ha/giờ và giảm thiểu từ 15% đến 20% nguy cơ lật trượt mất an toàn cho người vận hành.

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 kết hợp với cơ học đất và cơ học vật rắn nhiều khối lượng. Khung lý thuyết tập trung phân tích sự cân bằng lực và mô men của liên hợp máy theo phương trình vi phân chuyển động thiết lập bằng nguyên lý D'Alembert. Các lực tác động cơ bản bao gồm lực kéo tiếp tuyến tại bánh chủ động, lực cản lăn của nền đất, lực cản dốc, lực bám tiếp xúc và tổng lực cản phá đất tác dụng lên lưỡi cày trụ.

Mô hình phá hủy đất Mohr-Coulomb và bài toán biến dạng phẳng nửa không gian được sử dụng để xác định lực cản phá đất, lực cắt đáy xới và lực ma sát giữa đất với bề mặt lưỡi cày. Hệ thống được mô hình hóa thành hệ động lực học nhiều phần tử, bao gồm khối lượng phần đầu máy kéo 500 kg liên kết qua khớp quay với bộ phận cày công tác.

Các khái niệm then chốt trong nghiên cứu gồm:

  • Lực kéo tiếp tuyến: Lực sinh ra tại vết tiếp xúc giữa bánh xe chủ động và mặt đất do mô men quay động cơ 43 Nm truyền qua hệ thống truyền động có hiệu suất 0,563.
  • Khả năng bám: Giới hạn lực truyền cực đại phụ thuộc vào tải trọng bám và hệ số bám của từng loại đất.
  • Góc ổn định tĩnh và góc ổn định động: Giới hạn góc nghiêng địa hình lớn nhất mà liên hợp máy duy trì được trạng thái cân bằng không bị trượt dọc, trượt ngang hoặc lật đổ khi đứng yên và khi chuyển động.
  • Mô men cản lăn: Mô men cản trở sự lăn của bánh xe do biến dạng đàn dẻo của lốp và nền đất.

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

Nghiên cứu kết hợp chặt chẽ giữa phương pháp giải tích lý thuyết và phương pháp mô phỏng số. Dữ liệu đầu vào thu thập từ thông số kỹ thuật thực nghiệm của máy kéo Bông Sen 12 với bán kính làm việc bánh xe 0,28 m, chiều cao trọng tâm 641 mm và chiều rộng cơ sở điều chỉnh từ 530 mm đến 810 mm.

Cỡ mẫu dữ liệu thực nghiệm được thu thập trên 12 bộ mẫu thử nghiệm lặp lại trên 4 nhóm nền đất đặc trưng trong nông lâm nghiệp: đất khô có hệ số cản lăn 0,05 và hệ số bám 0,7; đất ẩm có hệ số cản lăn 0,06 và hệ số bám 0,6; đất hoang hóa có hệ số cản lăn 0,07 và hệ số bám 0,8; đất canh tác có hệ số cản lăn 0,12 và hệ số bám 0,6. Phương pháp chọn mẫu phân tầng có chủ đích được áp dụng nhằm bao phủ toàn bộ các dải cơ lý tính đất biến thiên tại các vùng lâm sinh đồi dốc.

Lý do lựa chọn phương pháp phân tích số trên phần mềm Matlab - Simulink là khả năng giải chính xác hệ phương trình vi phân phi tuyến bậc cao mô tả quỹ đạo chuyển động và chuyển vị xoay của máy kéo trên dốc. Phương pháp này khắc phục nhược điểm của các mô hình tuyến tính hóa tĩnh trước đây, cho phép mô phỏng chính xác đáp ứng động học liên tục theo thời gian thực. Toàn bộ quy trình tính toán và mô phỏng được thực hiện liên tục trong khung thời gian 12 tháng.

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

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

Khảo sát động lực học chỉ ra lực kéo tiếp tuyến của máy kéo Bông Sen 12 biến thiên mạnh theo 6 cấp số tiến và 2 cấp số lùi. Ở số tiến I (tỷ số truyền 174,07), lực kéo đạt cực đại 15.050 N; số II đạt 8.360 N; số III đạt 5.168 N; số IV đạt 3.987 N; số V đạt 2.041 N và số VI giảm xuống 1.305 N. Ở số lùi I (tỷ số truyền 259,35), lực kéo đạt mức 22.423 N. Khả năng kéo thực tế giảm dần theo hàm bậc nhất khi độ dốc tăng từ 0 độ lên 30 độ.

Khả năng bám của máy kéo đạt giá trị lớn nhất trên đất hoang hóa với lực bám cực đại trên 3.924 N, trong khi trên đất ẩm và đất canh tác lực bám chỉ đạt xấp xỉ 2.943 N (giảm 25%). Khi độ dốc địa hình tăng, lực bám suy giảm theo hàm cosin của góc dốc.

Kết quả tính toán góc ổn định tĩnh và ổn định động xác định được các giá trị tới hạn rõ rệt giữa các điều kiện làm việc:

  • Góc ổn định tĩnh dọc theo điều kiện trượt đạt 34 độ 39 phút trên đất hoang hóa, 33 độ 59 phút trên đất khô và 30 độ 58 phút trên đất canh tác; trong khi góc ổn định theo điều kiện lật đạt cố định 35 độ 05 phút.
  • Góc ổn định tĩnh ngang theo điều kiện trượt đạt 25 độ 39 phút trên đất hoang hóa, 22 độ 47 phút trên đất khô và 19 độ 48 phút trên đất canh tác; thấp hơn nhiều so với góc lật ngang 32 độ 17 phút.
  • Góc ổn định động dọc khi máy xuống dốc không tải giảm nghiêm trọng, chỉ còn 11 độ 01 phút trên đất hoang hóa, 10 độ 02 phút trên đất khô và 9 độ 25 phút trên đất canh tác. Khi lên dốc đủ tải, góc ổn định động đạt từ 23 độ 04 phút đến 33 độ 31 phút.

Thảo luận kết quả

Hiện tượng trượt luôn xuất hiện trước hiện tượng lật khi máy kéo đứng yên trên sườn dốc, bởi góc giới hạn trượt luôn nhỏ hơn góc giới hạn lật từ 1 độ đến 12 độ tùy thuộc loại đất. Nguyên nhân là do hệ số bám của đất giới hạn lực giữ tiếp tuyến trước khi mô men trọng trường làm triệt tiêu phản lực tựa pháp tuyến. Ngược lại, khi máy kéo chuyển động xuống dốc, hiện tượng lật đổ lại xuất hiện sớm hơn hiện tượng trượt do sự dịch chuyển trọng tâm về phía trước làm phản lực tựa tại điểm tỳ bánh sau bị triệt tiêu, khiến góc ổn định động giảm mạnh xuống dưới 10 độ.

So sánh với các mô hình máy kéo lớn 4 bánh của Muller hoặc Vogel, máy kéo nhỏ 2 bánh BS 12 có khối lượng chuyển động và mô men quán tính nhỏ hơn nhiều, dẫn tới độ nhạy chấn động và dao động lắc quanh trục bánh chủ động cao hơn từ 30% đến 45%.

Dữ liệu nghiên cứu được trình bày chi tiết qua hệ thống 6 bảng số liệu tra cứu và 16 đồ thị động học. Các đồ thị mô phỏng trên Matlab - Simulink biểu diễn trực quan mối quan hệ phi tuyến giữa độ dốc địa hình, lực cản cày và hệ số ma sát trượt tại ranh giới mất ổn định, cho phép xác định chính xác vùng làm việc an toàn của liên hợp máy.

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

Dựa trên các kết quả động lực học đã thiết lập, nghiên cứu đề xuất 4 nhóm giải pháp kỹ thuật và vận hành cụ thể:

  • Khống chế chế độ số truyền và tốc độ di chuyển trên sườn dốc: Người vận hành chỉ sử dụng tay số tiến cấp I và cấp II (tương ứng vận tốc 1,6 đến 2,9 km/h) khi cày đất trên sườn dốc từ 10 độ đến 18 độ. Tuyệt đối không vận hành ở các cấp số cao hơn nhằm đảm bảo mô men xoắn 43 Nm ổn định và giảm trên 25% hệ số trượt bánh xe. Mục tiêu giảm suất tiêu hao nhiên liệu xuống dưới 19 lít/ha, triển khai ngay trong quy trình vận hành thực tế.
  • Lắp đặt bánh lồng tăng cường khả năng bám: Các cơ sở cơ khí chế tạo máy nông nghiệp (như VIKYNO, DISOCO) cần tiêu chuẩn hóa bộ bánh lồng phụ trợ đường kính tương thích lốp 6,00-12, tăng bề rộng cơ sở lên mức tối đa 810 mm. Giải pháp này giúp nâng hệ số bám trên đất ẩm và đất canh tác từ 0,6 lên 0,75, tăng lực bám thêm 20%, triển khai chế tạo thử nghiệm trong vòng 6 tháng.
  • Hiệu chỉnh tọa độ trọng tâm và góc đặt cày trụ: Nhóm kỹ sư cơ khí Đại học Lâm nghiệp cần phối hợp với nhà sản xuất để thiết kế hạ thấp trọng tâm đầu máy từ 641 mm xuống dưới 580 mm bằng cách bố trí lại giá đỡ động cơ diesel, giúp nâng góc ổn định động khi xuống dốc thêm ít nhất 3 độ (đạt trên 12 độ), hoàn thành trong khung thời gian 12 tháng.
  • Ban hành quy chuẩn an toàn lao động trên đất dốc: Cục Cơ điện Nông nghiệp và Ngành nghề Nông thôn cần xây dựng tài liệu khuyến cáo giới hạn góc nghiêng làm việc tối đa cho máy kéo cầm tay không vượt quá 10 độ 55 phút (tương đương 50% góc ổn định tĩnh), nghiêm cấm quay vòng máy ở vận tốc cực đại 4,25 km/h trên địa hình dốc trên 15 độ, hoàn thành ban hành trong vòng 9 tháng.

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

Nội dung luận văn mang giá trị ứng dụng và học thuật cao cho 4 nhóm đối tượng trọng tâm:

  • Kỹ sư thiết kế và chuyên gia chế tạo máy nông lâm nghiệp: Áp dụng các công thức tính toán lực kéo tiếp tuyến, mô men cản lăn và ma sát đất để cải tiến hệ thống truyền động (hộp số 3 tiến 1 lùi nhân đôi), tối ưu hóa phân bố tải trọng và kết cấu khung vỏ máy kéo cầm tay.
  • Giảng viên, nghiên cứu sinh và sinh viên ngành Cơ khí động lực: Khai thác hệ phương trình vi phân D'Alembert và bộ mã mô phỏng Matlab - Simulink làm học liệu tham khảo chuyên sâu cho các môn học Động lực học ô tô máy kéo, Máy nông nghiệp và Phương pháp mô phỏng số.
  • Cán bộ kỹ thuật nông trường và chủ trang trại nông lâm nghiệp: Nắm vững các bảng tra cứu góc ổn định tĩnh và động trên từng loại đất để bố trí hướng cày theo đường đồng mức, lựa chọn cấp số và tốc độ vận hành hợp lý, nâng cao năng suất cày đạt 0,13 ha/giờ và hạn chế tai nạn lao động.
  • Nhà quản lý và hoạch định chính sách cơ giới hóa nông thôn: Sử dụng dữ liệu định lượng về phân bố công suất động lực (0,36 đến 1,78 mã lực/ha) để xây dựng chiến lược đầu tư máy móc phù hợp với đặc thù địa hình đồi núi dốc tại miền Bắc và Tây Nguyên.

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

Khả năng kéo lớn nhất của máy kéo Bông Sen 12 đạt bao nhiêu Newton ở các cấp số truyền?
Ở số tiến I với tỷ số truyền 174,07, máy sinh ra lực kéo tiếp tuyến lớn nhất đạt 15.050 N, đáp ứng tốt lực cản phá đất lớn. Ở số lùi I, lực kéo đạt 22.423 N. Lực kéo tiếp tuyến giảm dần ở các cấp số cao hơn, đạt 8.360 N ở số II và chỉ còn 1.305 N ở số tiến VI.

Tại sao máy kéo BS 12 dễ bị lật khi di chuyển xuống dốc hơn khi di chuyển lên dốc?
Khi xuống dốc, trọng tâm máy kéo dịch chuyển dồn tải trọng lên bánh trước, làm phản lực tựa tại bánh sau giảm nhanh về 0. Đồng thời mô men phanh và lực cản lăn tạo xu hướng lật quanh trục trước, khiến góc ổn định động khi xuống dốc giảm mạnh xuống chỉ còn 9 độ 25 phút đến 11 độ 01 phút.

Giới hạn độ dốc an toàn tối đa cho phép máy kéo BS 12 vận hành trên ruộng đồi là bao nhiêu?
Theo tiêu chuẩn kỹ thuật an toàn động lực học, góc nghiêng cho phép làm việc trên sườn dốc lấy bằng 50% góc ổn định tĩnh nhỏ nhất, tương đương 10 độ 55 phút (khoảng 11 độ). Vận hành vượt quá ngưỡng này khiến máy kéo dễ bị trượt ngang hoặc mất lái khi quay vòng ở tốc độ 4,25 km/h.

Giải pháp lắp bánh lồng kim loại mang lại hiệu quả động lực học như thế nào?
Bánh lồng tăng diện tích tiếp xúc và tăng độ bám sâu vào tầng đất cứng qua các mấu bám, giúp nâng hệ số bám từ 0,6 lên 0,72 - 0,75 trên nền đất ẩm và đất canh tác. Điều này làm tăng lực bám thêm 18% đến 22%, giảm độ trượt quay và tiết kiệm nhiên liệu khi làm đất.

Phần mềm Matlab - Simulink đóng vai trò gì trong việc giải quyết bài toán của luận văn?
Matlab - Simulink được dùng để lập trình giải hệ phương trình vi phân chuyển động phi tuyến mô tả chuyển vị xoay và trượt của máy kéo trên sườn dốc. Phần mềm xuất ra 16 sơ đồ và đồ thị phản ánh chính xác tương tác giữa độ dốc, lực cản cày và hệ số ma sát giới hạn.

Kết luận

  • Thiết lập thành công mô hình động lực học toàn diện cho liên hợp máy kéo Bông Sen 12 với cày trụ một lưỡi trên sườn dốc bằng phương pháp giải tích và nguyên lý D'Alembert.
  • Xác định chính xác dải lực kéo tiếp tuyến từ 1.305 N đến 15.050 N trên 6 cấp số tiến và lực bám giới hạn từ 2.943 N đến 3.924 N trên 4 nhóm đất canh tác.
  • Làm sáng tỏ ranh giới góc ổn định tĩnh dọc (30 độ 58 phút đến 34 độ 39 phút), ổn định tĩnh ngang (19 độ 48 phút đến 25 độ 39 phút) và chỉ ra nguy cơ mất ổn định lật động khi xuống dốc tại góc nghiêng 9 độ 25 phút.
  • Ứng dụng thành công công cụ mô phỏng Matlab - Simulink để xây dựng biểu đồ tương quan giữa độ dốc địa hình, lực cản phá đất và chuyển vị trượt - xoay của máy.
  • Đề xuất 4 nhóm giải pháp kỹ thuật khả thi giúp nâng cao năng suất cày đạt 0,13 ha/giờ, tối ưu hóa tiêu hao nhiên liệu 19 lít/ha và đảm bảo an toàn vận hành ở độ dốc dưới 10 độ 55 phút.

Luận văn đóng góp luận cứ khoa học quan trọng cho việc hoàn thiện thiết kế chế tạo máy kéo hai bánh nội địa. Trong 12 tháng tiếp theo, nhóm nghiên cứu khuyến nghị tiếp tục chế tạo thử nghiệm cụm đối trọng hạ thấp trọng tâm và hoàn thiện quy chuẩn vận hành trên đất dốc. Độc giả và các đơn vị sản xuất cơ khí nông nghiệp quan tâm có thể khai thác toàn bộ dữ liệu phương trình và thuật toán mô phỏng trong công trình để ứng dụng vào thực tiễn sản xuất.