Tổng quan nghiên cứu

Trong ngành lâm nghiệp hiện đại, khâu bốc dỡ gỗ rừng trồng là một trong những công đoạn nặng nhọc, nguy hiểm nhất, chiếm tới khoảng 40% chi phí nhân công và thời gian toàn bộ chu kỳ khai thác. Nhằm cơ giới hóa khâu sản xuất này, đề tài khoa học công nghệ cấp Nhà nước KC-07-26-05 đã nghiên cứu, chế tạo thử nghiệm thành công cơ cấu tay thủy lực gắn sau máy kéo bánh hơi Shibaura SD 2843 với tầm vươn tối đa 3,8 m, chiều cao nâng 4,3 m, tải trọng nâng danh định 2000 N và tổng trọng lượng bản thân kết cấu đạt 2950 N. Tuy nhiên, các tính toán thiết kế ban đầu chủ yếu dựa trên phương pháp sức bền vật liệu truyền thống và sử dụng các hệ số tải trọng động theo kinh nghiệm tĩnh, dẫn đến tình trạng kết cấu kim loại bị thừa bền ở các đoạn thân dầm nhưng lại tiềm ẩn nguy cơ mất an toàn tại các điểm tập trung ứng suất khi làm việc ở giai đoạn quá độ.

Vấn đề cốt lõi của nghiên cứu là xác định chính xác trạng thái ứng suất và biến dạng thực tế của các bộ phận cơ bản trên tay thủy lực trong các thời điểm chuyển tiếp tải trọng phức tạp, bao gồm: thời điểm bắt đầu nhấc tải khỏi mặt đất, thời điểm hạ tải kèm phanh hãm đột ngột, và thời điểm khởi động quay cần bốc. Mục tiêu cụ thể của luận văn thạc sĩ kỹ thuật do học viên Đinh Thị Hợi thực hiện dưới sự hướng dẫn của Phó Giáo sư, Tiến sĩ Nguyễn Nhật Chiêu tại Trường Đại học Lâm nghiệp (năm 2009) là xây dựng mô hình không gian ba chiều (3D), mô phỏng động lực học và ứng dụng phương pháp phần tử hữu hạn để đánh giá định lượng phân bố ứng suất, độ võng biến dạng của cẳng tay, cánh tay và trụ quay. Nghiên cứu cung cấp luận cứ khoa học chuẩn xác để tối ưu hóa kết cấu cơ khí theo hướng giảm 10% đến 15% khối lượng kim loại bản thân, nâng cao tải trọng nâng hữu ích vượt mức 2000 N và thiết lập chế độ làm việc tối ưu cho máy kéo trên địa hình sườn dốc lâm nghiệp 10 đến 20 độ.

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 được xây dựng dựa trên sự tích hợp của hai nền tảng lý thuyết cơ học kinh điển kết hợp với công nghệ cơ tin học hiện đại:

  • Lý thuyết phương pháp phần tử hữu hạn (FEM): Luận văn áp dụng mô hình phần tử hữu hạn lấy chuyển vị nút làm ẩn số gốc. Dưới tác dụng của ngoại lực, phương trình vi phân cân bằng của hệ kết cấu được thiết lập dựa trên định luật đàn hồi Hooke và tính tương thích biến dạng liên tục tại các điểm nút, từ đó suy ra ten-xơ biến dạng và trường ứng suất tương đương von Mises trên toàn bộ miền thể tích kết cấu.
  • Lý thuyết động lực học máy khai thác lâm nghiệp: Khung lý thuyết dựa trên các công trình nghiên cứu của các nhà khoa học lâm nghiệp thế giới, đặc biệt là trường phái cơ học máy Liên Xô với các mô hình phương trình vi phân dao động của hệ nhiều vật rắn liên hợp (máy kéo - tay thủy lực - tải trọng gỗ) mô tả bởi phương trình vi phân Lagrange loại II.

Các khái niệm then chốt được chuẩn hóa trong nghiên cứu bao gồm:

  • Giai đoạn quá độ (transient phase): Khoảng thời gian chuyển tiếp rất ngắn (từ 0,1 đến 0,5 giây) khi cơ cấu thay đổi đột ngột vận tốc và gia tốc, phát sinh lực quán tính lớn gấp nhiều lần tải trọng tĩnh.
  • Hệ số động lực học ($K$): Tỷ số giữa tải trọng động cực đại đo được trong quá trình chuyển động quá độ và tải trọng tĩnh danh định của hệ thống.
  • Phần tử khối bậc cao Solid 186: Kiểu phần tử đẳng tham số dạng khối 20 nút, mỗi nút có 3 bậc tự do tịnh tiến ($U_X, U_Y, U_Z$), cho phép mô tả chính xác ứng xử uốn, xoắn và trường ứng suất bậc hai của vật liệu thép hợp kim CT5 có mô đun đàn hồi $E = 215 \cdot 10^9 \text{ Pa}$, hệ số Poisson $\nu = 0,3$ và khối lượng riêng $\rho = 7850 \text{ kg/m}^3$.

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

Nghiên cứu sử dụng phương pháp kết hợp chặt chẽ giữa mô phỏng số giải tích và đo đạc thực nghiệm trên mô hình máy thực:

  • Nguồn dữ liệu và mô hình hóa: Toàn bộ kích thước hình học thực tế của 3 cụm chi tiết cơ bản gồm cẳng tay (chiều dài 2,5 m), cánh tay (chiều dài 2,4 m) và trụ quay (chiều cao 2,0 m) được mô hình hóa không gian 3D chi tiết trên phần mềm SolidWorks, sau đó chuyển đổi qua định dạng chuẩn IGES để nạp trực tiếp vào môi trường phần mềm phân tích phần tử hữu hạn ANSYS.
  • Cỡ mẫu và phương pháp chọn mẫu phân tích: Mô hình phần tử hữu hạn thực hiện phân chia mạng lưới tự do (free meshing) có kiểm soát toàn cục cấp độ 4 (Basis Global Control), tạo ra cỡ mẫu tính toán gồm hàng chục nghìn phần tử khối Solid 186 và nút tọa độ. Phương pháp chọn mẫu tập trung cục bộ được áp dụng tại các vùng tiếp giáp nguy hiểm như lỗ xỏ chốt trục bản lề, gân tăng cứng đế trụ đỡ và các tai gá xilanh thủy lực, nơi có gradient biến thiên ứng suất rất cao.
  • Phương pháp đo đạc thực nghiệm: Sử dụng phương pháp đo các đại lượng không điện bằng điện thông qua hệ thống cảm biến đo lực đầu cần, cảm biến ten-xơ đo biến dạng điện trở (strain gauge) gắn tại các điểm xung yếu trên cẳng tay và trụ quay. Dữ liệu thực nghiệm được thu thập, khuếch đại và xử lý thông qua hệ thu thập dữ liệu chuyên dụng DMC Plus / DMC Laplus kết hợp phần mềm Kaleidagraph và Microsoft Excel qua 15 chu kỳ thử nghiệm độc lập ở các chế độ tải trọng khác nhau.
  • Lý do lựa chọn phương pháp: Phương pháp phần tử hữu hạn trên ANSYS cho phép phân tích chính xác các chi tiết có hình dạng hình học không gian phức tạp chịu tải trọng đa trục mà các công thức giải tích sức bền vật liệu đơn giản không thể giải được; đồng thời, việc đối chứng với số liệu đo thực nghiệm trực tiếp trên máy kéo Shibaura SD 2843 đảm bảo độ tin cậy khoa học tuyệt đối cho các kết quả mô phỏng số.

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 trên phần mềm ANSYS kết hợp đối chứng thực nghiệm đo đạc lực đầu cần đã mang lại những phát hiện định lượng then chốt về trạng thái chịu tải của tay thủy lực ở các giai đoạn làm việc khác nhau:

  • Phát hiện 1 - Tải trọng động ở giai đoạn hạ tải và phanh đạt giá trị cao nhất: Thực nghiệm đo lực tại đầu cần ngoạm gỗ ở tầm vươn xa nhất 3,8 m và tải nâng tối đa 2000 N xác định được hệ số động lực học cực đại $K = 1,897$. Giá trị này cao hơn đáng kể so với giai đoạn bắt đầu nâng tải ($K = 1,667$) và giai đoạn bắt đầu quay cần ($K = 1,658$). Điều này chứng minh quá trình hãm phanh khi hạ tải tạo ra xung lực quán tính nguy hiểm nhất cho toàn bộ hệ thống, làm tăng tải trọng tác dụng thực tế lên tới 89,7% so với tính toán tĩnh.
  • Phát hiện 2 - Phản lực khớp liên kết xilanh tăng vọt trong giai đoạn quá độ: Tại khớp nối giữa trụ xilanh nâng hạ và cánh tay, lực phản lực tổng hợp theo phương thẳng đứng ($R_Z$) tăng từ mức tĩnh $42.601,6\text{ N}$ lên tới $71.034,72\text{ N}$ khi bắt đầu nâng tải (tăng 66,74%) và đạt đỉnh điểm $80.795,76\text{ N}$ trong quá trình hạ tải và phanh (tăng 89,65%). Tương tự, lực tác dụng lên gối ổ đỡ của trụ quay tăng từ $22.954\text{ N}$ (trạng thái tĩnh) lên mức cực đại $43.533,25\text{ N}$ ở giai đoạn phanh hãm.
  • Phát hiện 3 - Phân bố ứng suất tập trung không đều trên các dầm kim loại: Kết quả phân tích trường ứng suất tương đương von Mises chỉ ra rằng ứng suất cực đại xuất hiện cục bộ tại mép các lỗ lắp chốt xilanh nâng hạ cánh tay và các đường hàn gân gia cường tại chân đế trụ xoay, vượt ngưỡng 160 MPa ở giai đoạn hạ tải phanh. Trong khi đó, phần lớn diện tích thân hộp của cánh tay và cẳng tay chỉ chịu mức ứng suất trung bình dưới 45 MPa đến 65 MPa, thấp hơn rất nhiều so với giới hạn bền mỏi cho phép của thép CT5 (giới hạn chảy khoảng 235 MPa), chứng tỏ kết cấu thân dầm đang bị dư thừa kim loại nghiêm trọng.
  • Phát hiện 4 - Ảnh hưởng của mô men xoay khởi động: Khi tay thủy lực bắt đầu quay góc 60 độ, mô men chủ động từ động cơ thủy lực truyền tới trụ quay đạt giá trị $M = 956\text{ Nm}$, sinh ra lực tiếp tuyến tác động từ bộ truyền xích đôi lên trụ quay đạt $34.649,03\text{ N}$ và lực phản lực ngang $R_X = 49.483,55\text{ N}$ tại khớp nối cánh tay, tạo ra trạng thái uốn xoắn đồng thời phức tạp cho trục đứng.

Thảo luận kết quả

Nguyên nhân dẫn đến hiện tượng tải trọng động tăng vọt tại giai đoạn hạ tải và phanh là do vận tốc dòng chất lỏng công tác qua van phân phối bị chặn đột ngột khi đóng cửa van, kết hợp với quán tính khối lượng của bó gỗ 2000 N ở đầu cần có cánh tay đòn dài 3,8 m tạo ra dao động uốn đàn hồi với biên độ lớn.

Dữ liệu nghiên cứu được biểu diễn trực quan thông qua mô đun xử lý kết quả ANSYS POST1 dưới dạng các biểu đồ đường đồng mức màu phân bố ứng suất von Mises và trường véc-tơ chuyển vị tổng thể. Đồng thời, mô đun POST26 được ứng dụng để trích xuất các đồ thị hàm biến dạng theo thời gian tại các nút phần tử xung yếu. Dữ liệu này tương thích hoàn toàn với các đồ thị dao động lực thu được từ phần mềm Kaleidagraph trong quá trình xử lý tín hiệu cảm biến điện trở.

So sánh với các nghiên cứu trước đây (như các tính toán thiết kế của đề tài KC-07-26 vốn sử dụng hệ số động chung theo kinh nghiệm từ 1,25 đến 1,35), kết quả của luận văn đã chứng minh hệ số động thực tế lên tới 1,897. Sự chênh lệch khoảng 40% này giải thích rõ ràng lý do vì sao một số chi tiết tai gối chốt bị biến dạng cục bộ sau thời gian dài vận hành, trong khi khối lượng toàn bộ cụm tay bốc lại quá nặng (2950 N). Mức sai số giữa mô phỏng phần tử hữu hạn và thực nghiệm đo ten-xơ chỉ dao động trong khoảng 6,5% đến 9,2%, khẳng định mô hình số SolidWorks - ANSYS có độ tin cậy rất cao và phản ánh trung thực trạng thái làm việc thực tế của máy kéo.

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

Dựa trên các kết quả khảo sát ứng suất và biến dạng ở giai đoạn quá độ, luận văn đưa ra 4 nhóm giải pháp kỹ thuật cụ thể nhằm hoàn thiện kết cấu và tối ưu hóa vận hành tay thủy lực:

  • Tối ưu hóa hình học dầm chịu lực: Tiến hành giảm chiều dày thành thép tấm tổ hợp của dầm cánh tay (dài 2,4 m) và dầm cẳng tay (dài 2,5 m) tại các vùng chịu mô men uốn thấp từ 6 mm xuống 4,5 mm. Giải pháp này giúp cắt giảm từ 10% đến 14% tổng trọng lượng bản thân của tay thủy lực (giảm khoảng 300 N đến 410 N trên tổng 2950 N), thực hiện bởi nhóm kỹ sư phòng thiết kế cơ khí trong thời gian 3 tháng.
  • Gia cường cục bộ các vị trí tập trung ứng suất: Bổ sung các tấm đệm lót tăng cứng cục bộ dạng vành khuyên quanh các lỗ xỏ chốt trục xilanh nâng hạ cánh tay và mở rộng bán kính góc lượn gân đế trụ xoay từ $R = 5\text{ mm}$ lên $R = 12\text{ mm}$. Mục tiêu giảm đỉnh ứng suất cục bộ xuống dưới mức 140 MPa, hoàn thành chế tạo bổ sung trong vòng 1 tháng trước khi lắp ráp thử nghiệm.
  • Tích hợp van giảm chấn chống sốc thủy lực: Thiết kế lắp đặt bổ sung cụm van tiết lưu giảm chấn một chiều và van chống sốc thủy lực điều chỉnh áp suất tại các đường dầu dẫn vào xilanh nâng hạ và xilanh co duỗi. Mục tiêu kỹ thuật là triệt tiêu hiện tượng dội áp, hạ thấp hệ số động lực học khi phanh $K$ từ 1,897 xuống dưới ngưỡng 1,35, do xưởng cơ điện chế tạo thực hiện trong 2 tháng.
  • Ban hành quy chuẩn quy trình vận hành an toàn: Quy định người điều khiển máy kéo Shibaura SD 2843 không thực hiện đồng thời hai thao tác xoay cần và hạ tải đột ngột khi mang tải cực đại 2000 N ở tầm vươn xa 3,8 m; giới hạn tốc độ góc quay cần không vượt quá 0,15 rad/s. Hướng dẫn kỹ thuật này cần được các đơn vị lâm trường và xí nghiệp khai thác gỗ áp dụng ngay trong thực tế sản xuất.

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

Luận văn thạc sĩ kỹ thuật này là tài liệu chuyên khảo có giá trị ứng dụng cao cho 4 nhóm đối tượng chính:

  • Kỹ sư thiết kế và chuyên viên R&D máy công trình: Ứng dụng quy trình mô hình hóa tích hợp giữa SolidWorks, CosmosMotion và ANSYS để tính toán, tối ưu hóa kết cấu các loại cần trục, tay máy bốc dỡ thủy lực lắp trên máy kéo bánh hơi hoặc bánh xích.
  • Các doanh nghiệp cơ khí nông lâm nghiệp: Khai thác toàn bộ thông số kết cấu, kích thước phân mảnh và bản vẽ chi tiết để chế tạo nội địa hóa các dòng tay thủy lực cỡ nhỏ lắp sau máy kéo nông nghiệp, thay thế thiết bị nhập khẩu đắt tiền với giá thành giảm từ 45% đến 60%.
  • Giảng viên, nghiên cứu sinh và học viên cao học ngành Cơ khí: Sử dụng làm tài liệu tham khảo chuyên sâu trong giảng dạy và nghiên cứu về động lực học máy lâm nghiệp, lý thuyết tính toán phần tử hữu hạn và phương pháp đo đạc thực nghiệm các đại lượng động cơ học.
  • Cán bộ quản lý kỹ thuật tại các lâm trường khai thác gỗ: Nắm vững các giới hạn chịu tải động và các chế độ làm việc quá độ nguy hiểm để xây dựng quy trình bảo dưỡng định kỳ, kiểm tra vết nứt mỏi tại các khớp liên kết chốt và nâng cao tuổi thọ vận hành cho đội máy bốc dỡ.

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

1. Tại sao việc nghiên cứu giai đoạn quá độ lại mang tính quyết định đối với độ bền tay thủy lực? Trong quá trình vận hành bốc dỡ gỗ, các thời điểm bắt đầu nâng tải, hạ tải kèm phanh hoặc bắt đầu quay sinh ra gia tốc đột ngột, làm phát sinh lực quán tính lớn. Kết quả nghiên cứu chứng minh hệ số tải trọng động tại giai đoạn hạ tải phanh đạt tới $K = 1,897$, khiến phản lực xilanh tăng tới 89,65% so với tính toán tĩnh thông thường.

2. Phần mềm SolidWorks và ANSYS đóng vai trò gì trong quy trình tính toán của đề tài? SolidWorks được sử dụng để xây dựng mô hình hình học 3D trung thực của các chi tiết phức tạp như cẳng tay, cánh tay và trụ quay. Mô hình sau đó được xuất qua định dạng IGES sang ANSYS để chia lưới phần tử khối bậc cao Solid 186 (20 nút), gán tải trọng động và giải bài toán trường ứng suất - biến dạng chính xác.

3. Vì sao kết cấu ban đầu của tay thủy lực lại bị thừa bền ở thân dầm nhưng lại dễ hỏng ở các khớp? Do phương pháp thiết kế truyền thống áp dụng hệ số an toàn tĩnh đồng đều cho toàn bộ kết cấu. Mô phỏng ANSYS cho thấy ứng suất thân dầm chỉ đạt 45 đến 65 MPa (dư bền lớn), trong khi mép các lỗ chốt xilanh nâng hạ chịu ứng suất tập trung vượt 160 MPa do phải chịu lực nén đỉnh lên tới $80.795,76\text{ N}$.

4. Giải pháp giảm trọng lượng bản thân mang lại hiệu quả kinh tế - kỹ thuật ra sao? Việc giảm độ dày thép ở các vùng ít chịu lực giúp giảm khoảng 350 N trọng lượng kết cấu (từ 2950 N ban đầu). Điều này cho phép tăng trực tiếp tải trọng gỗ nâng hữu ích của mỗi chuyến bốc vượt mức 2000 N, đồng thời cải thiện độ ổn định chống lật dọc của máy kéo Shibaura SD 2843 khi làm việc trên sườn dốc 10 đến 20 độ.

5. Kết quả nghiên cứu của luận văn có thể áp dụng cho các dòng máy kéo nông nghiệp khác không? Hoàn toàn có thể áp dụng. Toàn bộ quy trình thuật toán mô hình hóa trên ANSYS, phương pháp chia lưới phần tử hữu hạn và các công thức xác định tải trọng động lực học trong luận văn đều có thể chuyển giao trực tiếp để thiết kế tay thủy lực cho các dòng máy kéo có công suất tương đương từ 25 đến 50 mã lực.

Kết luận

  • Xây dựng hoàn chỉnh mô hình số 3D: Luận văn đã thiết lập thành công mô hình hình học 3D trên SolidWorks và chuyển đổi chính xác sang mô hình phần tử hữu hạn trên ANSYS sử dụng phần tử khối 20 nút Solid 186 cho toàn bộ các cụm chi tiết chính của tay thủy lực.
  • Xác định chính xác hệ số động lực học: Thông qua đo đạc thực nghiệm, đề tài đã chỉ rõ hệ số tải trọng động cực đại xuất hiện ở giai đoạn hạ tải và phanh đạt $K = 1,897$, giai đoạn bắt đầu nâng tải đạt $K = 1,667$ và giai đoạn bắt đầu quay cần đạt $K = 1,658$.
  • Định vị các vùng tập trung ứng suất nguy hiểm: Phân tích định lượng đã xác định được phản lực xilanh nâng hạ đạt đỉnh $80.795,76\text{ N}$ và phản lực ổ đỡ trụ xoay đạt $43.533,25\text{ N}$, chỉ ra các vị trí then chốt cần gia cường góc lượn và vành đệm chốt.
  • Khẳng định tiềm năng tối ưu hóa giảm trọng lượng: Chứng minh cơ sở khoa học cho phép giảm từ 10% đến 15% khối lượng kim loại thân dầm cẳng tay và cánh tay, giúp giảm tải trọng bản thân và tăng năng suất bốc dỡ gỗ.
  • Đóng góp thực tiễn cho cơ giới hóa lâm nghiệp: Đề xuất giải pháp bổ sung van thủy lực giảm chấn và quy trình vận hành an toàn cho máy kéo Shibaura SD 2843 khi bốc xếp gỗ rừng trồng cự ly ngắn.

Kế hoạch phát triển tiếp theo dự kiến triển khai chế tạo mẫu tay thủy lực cải tiến trong vòng 6 đến 12 tháng, kết hợp thử nghiệm kiểm định độ bền mỏi thực tế tại các lâm trường vùng núi phía Bắc. Quý bạn đọc, nhà nghiên cứu và kỹ sư quan tâm hãy tải toàn văn tài liệu luận văn thạc sĩ kỹ thuật của tác giả Đinh Thị Hợi để tra cứu toàn bộ ma trận số liệu tính toán và ứng dụng ngay vào các dự án cơ giới hóa máy lâm nghiệp.