Tổng quan nghiên cứu

Thảm họa cháy rừng hàng năm gây thiệt hại hàng nghìn hecta tài nguyên và phá hủy nghiêm trọng hệ sinh thái tự nhiên. Nhằm nâng cao hiệu quả kiểm soát hỏa hoạn tại các địa hình đồi núi phức tạp, đề tài cấp Nhà nước mã số KC 07.13/06-10 do Trường Đại học Lâm nghiệp chủ trì đã thiết kế và chế tạo thành công mẫu xe chữa cháy rừng đa năng trên nền tảng xe tải việt dã Ural-4320 (6x6) có tổng khối lượng 16.615 kg và công suất động cơ 132 kW. Tuy nhiên, khi vận hành mở băng cản lửa trên bề mặt đất rừng gồ ghề với độ mấp mô biến thiên từ 0,15 m đến 0,20 m, sự kết hợp giữa kích động biên dạng địa hình và các xung lực cơ cấu làm việc tạo ra dao động phức tạp, đe dọa trực tiếp đến tính êm dịu chuyển động, độ bền kết cấu và an toàn của kíp lái gồm 5 đến 7 người.

Mục tiêu cốt lõi của nghiên cứu là thiết lập mô hình động lực học dao động trong mặt phẳng đối xứng dọc của xe chữa cháy rừng đa năng, từ đó giải và mô phỏng hệ phương trình vi phân dao động nhằm làm sáng tỏ ảnh hưởng của vận tốc di chuyển, độ gồ ghề mặt đất và xung lực công tác đến trạng thái động lực học của xe. Phạm vi nghiên cứu tập trung vào giai đoạn xe di chuyển thi công băng cản lửa thông qua tổ hợp cưa cắt cây bụi phía trước (tốc độ quay 3.000 vòng/phút) và hệ thống phay cỏ rác phía sau (tốc độ 1.000 vòng/phút). Nghiên cứu được thực nghiệm đối chứng trực tiếp tại khu vực rừng núi Luốt, Xuân Mai, Hà Nội trong giai đoạn 2010 - 2011. Kết quả của công trình mang ý nghĩa khoa học và thực tiễn sâu sắc khi cung cấp căn cứ định lượng giúp kiểm soát hệ số tải trọng động dưới ngưỡng an toàn 3,0, nâng cao tuổi thọ cụm chi tiết và tối ưu hóa chế độ vận hành cho phương tiện chữa cháy rừng chuyên dụng đầu tiên tại Việt Nam.

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 nền tảng lý thuyết dao động cơ học kỹ thuật và động lực học hệ nhiều vật (Multi-body Dynamics) ứng dụng cho ô tô máy kéo bánh hơi. Để mô tả chính xác chuyển động tương đối của các phần tử trong không gian phẳng đối xứng dọc, tác giả áp dụng phương trình vi phân Lagrange loại 2 dạng tổng quát:

$$\frac{d}{dt}\left(\frac{\partial T}{\partial \dot{q}_i}\right) - \frac{\partial T}{\partial q_i} + \frac{\partial \Pi}{\partial q_i} + \frac{\partial \Phi}{\partial \dot{q}_i} = Q_i^*$$

Hệ thống lý thuyết tích hợp ba mô hình cốt lõi: mô hình tương tác liên hợp "Đường - Xe - Người", mô hình biến dạng đàn hồi tuyến tính của lốp cao su 14.00-20 kết hợp hệ thống treo nhíp lá, và mô hình kích động điều hòa từ bề mặt địa hình lâm nghiệp.

Nghiên cứu vận dụng 5 khái niệm chuyên ngành then chốt bao gồm:

  1. Hệ 6 tọa độ suy rộng: Xác định chuyển vị thẳng đứng của trọng tâm thân xe ($z_0$), chuyển vị góc lắc dọc của thân xe ($\alpha_0$), chuyển vị thẳng đứng của chốt xoay nhíp sau ($z_C$), chuyển vị góc của nhíp sau ($\alpha$), cùng chuyển vị thẳng đứng của trục cầu trước ($z_1$) và cụm phay cỏ ($z_4$).
  2. Hệ số tải trọng động: Tỷ số giữa lực động cực đại tác dụng lên cơ cấu và tải trọng tĩnh tương ứng.
  3. Hàm kích động mặt đường ngẫu nhiên quy đổi điều hòa: Quy luật mô tả biên dạng mấp mô theo bước sóng và biên độ.
  4. Hàm hao tán Rayleigh: Đại lượng biểu thị mức độ tiêu tán năng lượng do hệ số cản giảm chấn nhớt của lốp và hệ thống treo sinh ra.
  5. Lực suy rộng không thế: Tổng hợp các ngoại lực và xung lực cắt gọt tác động tức thời lên các cụm công tác.

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

Nguồn dữ liệu của đề tài được tổng hợp từ hai kênh chính: dữ liệu trắc đạc hình học biên dạng mặt đất rừng tự nhiên và dữ liệu đo đạc thông số dao động thực nghiệm. Nhóm nghiên cứu tiến hành khảo sát trắc đạc thực địa theo phương pháp chia đoạn chuẩn hóa với cự ly 50 cm trên toàn bộ tuyến đường thử nghiệm tại rừng núi Luốt để xác định quan hệ giữa chiều dài quãng đường $s$ và chiều cao mấp mô $q$. Mối tương quan phi tuyến được xử lý và kiểm định thông qua thủ tục hồi quy Correlation Test trên phần mềm thống kê chuyên dụng XLSTAT, từ đó chuẩn hóa thành hàm kích động điều hòa có biên độ $h_0$ từ 0,15 m đến 0,20 m và bước sóng quy chuẩn $s_0$.

Về phương pháp giải tích số, hệ phương trình vi phân phi tuyến bậc hai 6 bậc tự do được lập trình giải tự động trên môi trường phần mềm Mathematica, cho phép mô phỏng toàn diện trường gia tốc và chuyển vị trọng tâm với các kịch bản vận tốc biến thiên từ 3 km/h, 5 km/h, 10 km/h đến 15 km/h. Về phương pháp thực nghiệm, tác giả áp dụng công nghệ đo lường đại lượng không điện bằng tín hiệu điện qua bộ thu thập số liệu đa kênh DMC Plus (HBM - Đức) liên kết máy vi tính qua giao diện cáp RS-232, kết hợp cảm biến gia tốc điện cảm có dải tần đáp ứng cao gắn tại trọng tâm khối lượng được treo (13.300 kg). Cỡ mẫu thực nghiệm bao gồm 4 kịch bản kiểm tra lặp lại nhiều lần trên cùng một địa hình chuẩn hóa: chạy không tải công tác ở 3 mức tốc độ (3 km/h, 5 km/h, 10 km/h) và chạy có tải phay cỏ rác ở tốc độ 5 km/h. Việc lựa chọn phương pháp này đảm bảo tính tương thích cao nhất giữa mô hình số trị và điều kiện làm việc thực tế ngoài hiện trường.

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

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

Nghiên cứu đã đưa ra các kết quả định lượng chuẩn xác về trạng thái động lực học của xe chữa cháy rừng đa năng thông qua 4 phát hiện quan trọng:

  1. Biên độ mấp mô mặt đất tỷ lệ thuận với gia tốc dao động: Khi xe chuyển động ở vận tốc cố định 5 km/h, việc tăng chiều cao mấp mô mặt đất $h_0$ từ 0,15 m lên 0,20 m làm gia tốc dao động trọng tâm tăng tuyến tính khoảng 33,3%, từ mức cực đại 1,56 m/s² lên 2,08 m/s², đòi hỏi người lái phải giảm tốc độ tức thì khi gặp chướng ngại vật lớn hơn 0,20 m.
  2. Vận tốc di chuyển tạo ra bước nhảy vọt về tải trọng động: Ở dải tốc độ thấp từ 3 km/h đến 5 km/h, gia tốc trọng tâm dao động ổn định trong khoảng 0,85 m/s² đến 1,56 m/s². Tuy nhiên, khi vận tốc tăng lên 10 km/h, gia tốc đỉnh đạt 7,62 m/s² (tăng 388% so với mức 5 km/h). Khi vận tốc đạt 15 km/h, gia tốc bùng nổ lên tới 150 m/s², khiến hệ số tải trọng động vọt lên mức 16,0, vượt quá 530% so với giới hạn an toàn cho phép của kết cấu cơ khí (ngưỡng 2,5 - 3,0).
  3. Xung lực phay cỏ có mức độ ảnh hưởng thứ cấp: Khi kích hoạt hệ thống làm sạch cỏ rác với lực tác dụng thẳng đứng cực đại $S_z = 3.178,66\text{ N}$ ở vận tốc 5 km/h, gia tốc dao động thẳng đứng tại trọng tâm xe chỉ tăng nhẹ từ 1,56 m/s² lên 1,72 m/s² (mức tăng xấp xỉ 10,25%), chứng minh cơ cấu phay va đập không gây mất ổn định nghiêm trọng đến khung gầm chính.
  4. Độ tương đồng cao giữa lý thuyết và thực nghiệm: Kết quả phân tích phổ dữ liệu xuất từ phần mềm CATMAN trên hệ thống DMC Plus cho thấy sai số biên độ gia tốc giữa mô hình giải tích Mathematica và đo đạc thực tế chỉ dao động từ 8% đến 12%, khẳng định độ tin cậy tuyệt đối của hệ phương trình động học 6 bậc tự do.

Thảo luận kết quả

Hiện tượng gia tốc dao động tăng vọt phi tuyến tính ở dải vận tốc trên 10 km/h được giải thích do tần số kích động từ mấp mô mặt đất rừng ($f = v / s_0$) tiến sát vào dải tần số dao động riêng của phần khối lượng được treo (13.300 kg) kết hợp với độ cứng nhíp trước ($c_{n1} = 210.000\text{ N/m}$) và nhíp sau ($c_{n2} = 808.000\text{ N/m}$). Sự cộng hưởng này khuếch đại lực quán tính, gây nguy cơ biến dạng dẻo khung sát-xi, nứt vỡ đường ống công nghệ và gây chấn thương cột sống cho người vận hành.

So sánh với các công trình nghiên cứu trước đây về máy kéo lâm nghiệp của các tác giả trong và ngoài nước, xe chữa cháy rừng đa năng có tính chất động học phức tạp hơn rõ rệt do tải trọng tĩnh lớn (16.615 kg) phân bố trên 3 trục cầu việt dã và ảnh hưởng của khối chất lỏng 5,5 m³ dao động trong bồn chứa. Trong các báo cáo chuyên ngành, dữ liệu này thường được biểu diễn thông qua đồ thị miền thời gian mô tả biên độ gia tốc $z_0''(t)$ và bảng so sánh ma trận độ cứng quy đổi. Kết quả nghiên cứu xác nhận rằng độ mấp mô địa hình rừng là nguồn kích động nguy hiểm nhất, trong khi lực cắt gọt cơ khí của đĩa cưa và dao phay chỉ đóng vai trò kích động cục bộ tần số cao.

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

Dựa trên các kết luận rút ra từ mô hình động lực học và dữ liệu kiểm chứng thực tế, 4 giải pháp kỹ thuật và vận hành được đề xuất:

  1. Quy chuẩn hóa dải tốc độ làm việc an toàn: Ban hành quy trình vận hành nghiêm ngặt, khống chế vận tốc di chuyển của xe khi thi công tạo băng cản lửa trong khoảng từ 3 km/h đến tối đa 5 km/h (không vượt quá 8 km/h trên địa hình có độ dốc dưới 15 độ). Giải pháp này giúp duy trì hệ số tải trọng động dưới mức 2,5 và gia tốc trọng tâm dưới 1,8 m/s². Thời gian áp dụng: Ngay lập tức. Chủ thể thực hiện: Lực lượng Kiểm lâm và kíp lái trực tiếp.
  2. Cải tiến hệ thống treo phụ trợ và ghế ngồi giảm xóc: Lắp đặt bổ sung hệ thống giảm chấn thủy lực hai chiều cho cầu sau và trang bị ghế lái chuyên dụng có cơ cấu giảm chấn khí nén với tỷ số tần số kích động trên tần số riêng đạt 0,5 - 0,6 nhằm triệt tiêu 40% rung chấn truyền lên cơ thể người lái. Thời gian triển khai: 6 - 12 tháng. Chủ thể thực hiện: Đơn vị chế tạo cơ khí và nhà máy ô tô Mê Kông.
  3. Tối ưu hóa vách ngăn chống xóc trong bồn chứa nước: Thiết kế bổ sung hệ thống vách ngăn tiêu sóng đa khoang cho bồn chứa dung tích 5,5 m³ nhằm giảm thiểu hiện tượng dồn dịch trọng tâm tức thời của khối nước khi vượt chướng ngại vật mấp mô cao 0,20 m, giúp tăng 15% tính ổn định chống lật. Thời gian hoàn thành: Trong kỳ bảo dưỡng 3 tháng. Chủ thể thực hiện: Xưởng bảo dưỡng kỹ thuật lâm nghiệp.
  4. Tích hợp cảm biến cảnh báo quá ngưỡng dao động: Lắp đặt thiết bị đo gia tốc thời gian thực tại trọng tâm cabin liên kết còi báo động để tự động nhắc nhở tài xế giảm tốc khi gia tốc rung vượt quá ngưỡng 2,0 m/s². Thời gian nghiên cứu ứng dụng: 12 - 18 tháng. Chủ thể thực hiện: Viện Nghiên cứu Cơ điện Lâm nghiệp kết hợp Cục Kiểm lâm.

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

Luận văn là tài liệu chuyên khảo mang giá trị học thuật và ứng dụng cao cho 4 nhóm đối tượng:

  1. Kỹ sư thiết kế và chế tạo máy cơ giới nông - lâm nghiệp: Nắm vững quy trình thiết lập ma trận khối lượng, độ cứng và giảm chấn 6 bậc tự do; khai thác trực tiếp bảng thông số kỹ thuật quy đổi (độ cứng nhíp, lốp xe 14.00-20) phục vụ tính toán độ bền mỏi cho các dòng xe công trình việt dã 6x6.
  2. Cán bộ quản lý phòng cháy chữa cháy rừng và lực lượng Kiểm lâm: Tiếp cận cơ sở khoa học để xây dựng cẩm nang vận hành xe chữa cháy đa năng, lựa chọn chế độ chạy máy phù hợp với từng cấp độ dốc địa hình và tổ chức thi công băng cản lửa an toàn, đạt hiệu suất cắt dọn thực bì trên 1,5 km/h.
  3. Học viên cao học và nghiên cứu sinh ngành Kỹ thuật cơ khí động lực: Sử dụng luận văn như một công trình mẫu mực về việc kết hợp giữa phương pháp giải tích vi phân Lagrange trên Mathematica và công nghệ đo đạc thực nghiệm hiện đại qua thiết bị DMC Plus/CATMAN.
  4. Các doanh nghiệp và cơ sở sản xuất ô tô chuyên dùng: Vận dụng phương pháp mô hình hóa dao động phẳng để nâng cấp khung gầm xe tải Ural-4320 thành các dòng xe cứu hộ cứu nạn, xe địa hình quân sự và máy làm đất chuyên dụng.

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

Tại sao tốc độ di chuyển của xe chữa cháy rừng khi tạo băng cản lửa không được vượt quá 10 km/h?
Kết quả tính toán động lực học chứng minh ở tốc độ 10 km/h, gia tốc dao động đạt 7,62 m/s² và tăng vọt lên 150 m/s² tại 15 km/h. Vận tốc lớn gây hiện tượng cộng hưởng giữa mấp mô địa hình và hệ thống treo, đẩy hệ số tải trọng động lên mức 16,0, vượt xa giới hạn an toàn cho phép (2,5 - 3,0), gây nguy cơ gãy nhíp và lật xe.

Xung lực sinh ra từ hệ thống phay cỏ rác phía sau ảnh hưởng như thế nào đến toàn bộ xe?
Hệ thống phay gồm 10 đĩa và 20 dao tạo ra lực cắt thẳng đứng biến thiên với chu kỳ $S_z = 3.178,66(1 - \cos 20vt)\text{ N}$. Do tần số cắt cao và phân tán lực theo chu kỳ quay, xung lực này chỉ làm tăng gia tốc dao động tổng thể của xe thêm khoảng 10,25% (từ 1,56 m/s² lên 1,72 m/s² ở 5 km/h), không gây nguy hiểm cho thân xe.

Mô hình toán học trong luận văn giải quyết bài toán dao động bằng phương pháp nào?
Tác giả thiết lập hệ 6 phương trình vi phân chuyển động thông qua phương trình Lagrange loại 2, mô tả chuyển động của 5 khối lượng quy đổi với 6 tọa độ suy rộng ($z_0, \alpha_0, z_C, \alpha, z_1, z_4$). Hệ phương trình sau đó được giải số học tự động bằng các thuật toán tích phân vi phân trên phần mềm chuyên dụng Mathematica.

Khối lượng bồn nước 5,5 m³ được xử lý như thế nào trong mô hình lý thuyết?
Do téc nước 5.500 lít được thiết kế các vách ngăn cố định bên trong nên chất lỏng chỉ dịch chuyển cục bộ trong từng khoang hẹp. Khi xe dao động nhỏ quanh vị trí cân bằng tĩnh, khối nước cùng thân xe được giả thiết quy đổi thành một vật rắn tuyệt đối duy nhất có khối lượng được treo $m_0 = 13.300\text{ kg}$ đặt tại trọng tâm xe.

Thiết bị thực nghiệm nào đã được sử dụng để xác định gia tốc dao động của xe?
Tác giả sử dụng cảm biến gia tốc nguyên lý điện cảm chính xác cao của hãng HBM (Đức) gắn cố định tại trọng tâm xe, kết nối với bộ thu thập và khuếch đại thông tin đo lường đa năng DMC Plus. Dữ liệu thời gian thực được truyền qua cổng RS-232 và ghi nhận, phân tích bằng phần mềm chuyên dụng CATMAN.

Kết luận

  • Luận văn đã xây dựng thành công mô hình động lực học dao động 6 bậc tự do cho xe chữa cháy rừng đa năng trên cơ sở gầm xe Ural-4320 trong mặt phẳng đối xứng dọc.
  • Giải pháp mô phỏng số trên phần mềm Mathematica kết hợp trắc đạc địa hình đã xác định rõ quan hệ tỷ lệ thuận giữa biên độ mấp mô mặt đất rừng và gia tốc dao động của xe.
  • Nghiên cứu chỉ ra ranh giới vận tốc an toàn khi xe thực hiện tạo băng cản lửa là từ 3 km/h đến 5 km/h; nghiêm cấm chạy vượt quá 10 km/h để tránh phá hủy kết cấu do tải trọng động tăng vọt gấp 16 lần.
  • Kết quả thực nghiệm tại rừng núi Luốt bằng hệ thống DMC Plus và phần mềm CATMAN đạt độ chuẩn xác cao với sai số so với lý thuyết chỉ từ 8% đến 12%.
  • Bước phát triển tiếp theo là hoàn thiện hệ thống giảm chấn thứ cấp cho hệ thống treo và số hóa dữ liệu địa hình mặt đất rừng trên toàn quốc.

Hãy khai thác ngay tài liệu nghiên cứu này để làm chủ công nghệ tính toán dao động và tối ưu hóa quy trình thiết kế, vận hành các dòng phương tiện cơ giới chuyên dụng phục vụ công tác phòng chống cháy rừng hiệu quả.