Tổng quan nghiên cứu

Cháy rừng là thảm họa môi trường nghiêm trọng, đe dọa trực tiếp đến hệ sinh thái và gây thiệt hại lớn về kinh tế với hàng nghìn hecta rừng bị thiêu rụi mỗi năm. Tại Việt Nam, địa hình đồi núi hiểm trở và mạng lưới đường lâm nghiệp gồ ghề khiến vận tốc di chuyển trung bình của các phương tiện cơ giới chữa cháy suy giảm từ 40% đến 50% so với khi hoạt động trên đường giao thông thông thường. Nhằm khắc phục hạn chế thiếu hụt phương tiện chữa cháy cơ giới chuyên dụng trong nước, Đề 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 cơ sở khung gầm xe việt dã ba cầu chủ động Ural-4320.

Tuy nhiên, do mang theo bồn chứa dung tích 5,5 mét khối nước cùng hệ thống thiết bị công tác thủy lực nặng nề (tổng khối lượng toàn bộ xe khi chở đầy nước lên tới 16.615 kg), các kích động mấp mô từ mặt đường đất rừng tạo ra tải trọng động biến đổi chu kỳ rất lớn. Dao động mạnh không chỉ làm giảm độ êm dịu chuyển động, gây mệt mỏi cho tổ vận hành từ 5 đến 7 người mà còn đẩy nhanh hiện tượng mỏi kim loại, giảm tuổi thọ cụm chi tiết và tiềm ẩn nguy cơ mất ổn định chống lật.

Luận văn thạc sĩ kỹ thuật của tác giả Lê Thị Minh Vượng, thực hiện năm 2010 tại Trường Đại học Lâm nghiệp dưới sự hướng dẫn khoa học của PGS. Nguyễn Nhật Chiêu, tập trung giải quyết bài toán cốt lõi: Xây dựng mô hình toán học mô tả dao động, lập và giải hệ phương trình vi phân dao động của xe chữa cháy rừng đa năng trong mặt phẳng đối xứng dọc khi di chuyển trên đường lâm nghiệp, kết hợp đo kiểm thực nghiệm gia tốc dao động thẳng đứng. Nghiên cứu mang ý nghĩa thực tiễn then chốt trong việc cung cấp cơ sở dữ liệu định lượng để hoàn thiện hệ thống treo, tối ưu hóa phân bố tải trọng và xác định dải vận tốc vận hành an toàn cho mẫu xe chữa cháy rừng đầu tiên được chế tạo 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

Luận văn kế thừa và phát triển lý thuyết dao động ô tô máy kéo hiện đại dựa trên mô hình tổng thể tương tác "Đường - Xe - Người" được khởi xướng bởi Mitschke (1970) và Schiehlen (1980). Để phân tích các chuyển động phức hợp của xe nhiều cầu trên địa hình lâm nghiệp, tác giả áp dụng nguyên lý cơ học giải tích với phương trình vi phân Lagrange loại II. Mô hình động lực học được xây dựng trong hệ trục tọa độ không gian Descartes OXYZ, tập trung khảo sát dao động trong mặt phẳng thẳng đứng dọc OXZ đi qua trọng tâm xe với 5 tọa độ suy rộng tương ứng 5 bậc tự do độc lập.

Các khái niệm chuyên ngành nền tảng được định nghĩa và tính toán chuẩn xác bao gồm:

  1. Độ êm dịu chuyển động (Ride comfort): Chỉ tiêu đánh giá thông qua tần số dao động, gia tốc dao động thẳng đứng và thời gian tác động lên cơ thể người vận hành theo tiêu chuẩn ngành 22 TCN 336-06.
  2. Khối lượng được treo ($m = 13.300\text{ kg}$) và khối lượng không được treo ($m_s = 3.315\text{ kg}$): Tỷ số phân bố khối lượng $\gamma \approx 4,01$ đóng vai trò quyết định đến việc giảm xung lực truyền từ mặt đường lên khung vỏ.
  3. Mô men quán tính ($J = 24.212\text{ kg}\cdot\text{m}^2$): Đại lượng đặc trưng cho khả năng chống chuyển vị góc lắc dọc của thân xe quanh trục ngang OY.
  4. Hàm kích động mấp mô mặt đường: Mô hình hóa biên dạng trắc diện dọc đường đất rừng dưới dạng hàm tuần hoàn sóng điều hòa với biên độ dao động từ 1,8 cm đến 4,0 cm và bước sóng từ 0,6 m đến 1,0 m.

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 mô hình hóa lý thuyết giải tích và phương pháp đo đạc thực nghiệm các đại lượng không điện bằng tín hiệu điện:

  • Nguồn dữ liệu và thông số đầu vào: Toàn bộ thông số kết cấu được trích xuất từ mẫu xe thực tế thuộc Đề tài KC 07.13/06-10 và kiểm chứng qua các nghiên cứu của Học viện Kỹ thuật Quân sự. Xe trang bị động cơ công suất 132 kW, mô-men xoắn cực đại 667 N.m, 3 cầu chủ động 6x6, khối lượng cầu trước 655 kg, hai cầu sau mỗi cầu 590 kg, khoảng cách từ trọng tâm đến cầu trước $l_1 = 2,79\text{ m}$, đến điểm xoay nhíp sau $l_2 = 1,46\text{ m}$.
  • Phương pháp chọn mẫu và khảo sát mặt đường: Cỡ mẫu khảo sát thực địa bao gồm 4 dạng mặt đường lâm nghiệp đặc trưng có trắc diện xác định với độ dốc di chuyển dưới 15 độ. Dữ liệu kích động đầu vào gồm dạng xung đơn chiều cao gờ 0,2 m và dạng hàm sin liên tục chiều cao mấp mô 0,04 m bước sóng 1,0 m.
  • Lý do lựa chọn phương pháp phân tích: Phương pháp phương trình Lagrange loại II được lựa chọn vì tính ưu việt tuyệt đối trong việc thiết lập hệ phương trình vi phân cho các cơ hệ holonom dừng có nhiều liên kết đàn hồi và cản nhớt phi tuyến. Hệ phương trình vi phân dạng ma trận $M\ddot{q} + B\dot{q} + Cq = P\dot{h} + Qh$ được giải và mô phỏng số thông qua phần mềm MATLAB 7.0 kết hợp công cụ đồ họa khối Simulink, cho phép phân tích tức thời đáp ứng động lực học của hệ thống treo nhíp lá và nhíp đòn gánh cân bằng cầu sau.
  • Phương pháp thực nghiệm kiểm chứng: Thực hiện đo gia tốc dao động thẳng đứng của thân xe bằng hệ thống chuyển đổi đo lường số gồm cảm biến gia tốc áp điện gắn tại trọng tâm, khối khuếch đại tín hiệu kỹ thuật số DMC Plus và phần mềm điều khiển chuyên dụng DMC Laplus. Quá trình đo thực nghiệm lặp lại 5 lần trên cùng một cung đường thử nghiệm tại nhà máy để triệt tiêu sai số ngẫu nhiên, sau đó xuất dữ liệu dưới định dạng ASCII để đối chiếu trực tiếp với 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

Phân tích mô phỏng trên MATLAB/Simulink kết hợp dữ liệu thực nghiệm đã xác lập 4 phát hiện kỹ thuật mang tính đột phá về động lực học của dòng xe chữa cháy rừng chuyên dùng:

  1. Khả năng dập tắt dao động thẳng đứng của trọng tâm xe: Khi xe gặp kích thích xung bước nhảy đơn với chiều cao mấp mô $h = 0,2\text{ m}$, chuyển dịch thẳng đứng $z$ của trọng tâm đạt giá trị cực đại ban đầu là 0,015 m, sau đó suy giảm nhanh chóng và tắt hẳn hoàn toàn sau thời gian khoảng 2,8 giây. Biên độ dao động giảm hơn 65% ngay sau chu kỳ đầu tiên nhờ sự phối hợp tối ưu giữa hệ thống giảm chấn thủy lực cầu trước ($k_{n1} = 2.400\text{ N}\cdot\text{s/m}$) và độ cứng nhíp sau ($c_{n2} = 808.000\text{ N/m}$).
  2. Kiểm soát tuyệt vời chuyển vị góc lắc dọc của thân xe: Chuyển vị góc lắc dọc thân xe quanh trục ngang OY tại vị trí trọng tâm chỉ dao động trong biên độ cực nhỏ từ -0,004 rad đến +0,010 rad và triệt tiêu sau 3,0 giây. Điều này khẳng định kết cấu cụm treo nhíp đòn gánh cân bằng liên kết giữa cầu giữa và cầu sau (với mô men quán tính nhíp sau $J_c = 289\text{ kg}\cdot\text{m}^2$) có khả năng triệt tiêu góc lắc dọc cực kỳ hiệu quả, hạn chế tối đa nguy cơ lật xe khi chạy trên dốc nghiêng lâm nghiệp.
  3. Phản ứng dao động của cầu trước và cầu sau: Tâm cầu trước có chuyển dịch thẳng đứng cực đại $z_{11} \approx 0,005\text{ m}$ dưới kích động đơn, trong khi điểm xoay nhíp sau $z_c$ ghi nhận chuyển dịch cực đại $0,016\text{ m}$ ($1,6\times 10^{-2}\text{ m}$). Mặc dù chịu tải trọng lớn hơn (phản lực tĩnh cầu sau $G_2 = 11.005\text{ kg}$, chiếm 66,2% tổng trọng lượng xe), cụm cầu sau vẫn đảm bảo sự cân bằng đàn hồi ổn định nhờ cấu trúc hai nhíp đòn gánh đối xứng.
  4. Mô tả tương quan số liệu qua đồ thị và bảng ma trận: Dữ liệu mô phỏng từ các khối Scope trong Simulink hiển thị đường cong dao động tắt dần dạng hàm mũ rõ nét. Bảng ma trận độ cứng $C$ và ma trận cản $B$ với 5 ẩn số nút phản ánh tính đối xứng cao, minh chứng các tham số quy đổi lốp xe ($c_1 = c_2 = c_3 = 1.060.000\text{ N/m}$) và giảm chấn lốp ($k_1 = k_2 = k_3 = 1.000\text{ N}\cdot\text{s/m}$) hoàn toàn tương thích với điều kiện tải trọng danh định của khung gầm Ural-4320.

Thảo luận kết quả

Nguyên nhân chính giúp xe chữa cháy rừng đa năng duy trì tính ổn định vượt trội là nhờ cơ cấu treo đòn cân bằng tại hai cầu sau. Khi một bánh xe sau rơi vào hố mấp mô hoặc leo qua mô đất cao $0,2\text{ m}$, đòn gánh sẽ tự xoay quanh trục liên kết A một góc $\theta$, làm giảm biên độ dịch chuyển truyền thẳng đứng lên khung xe xuống còn một nửa so với hệ thống treo độc lập thông thường.

Khi so sánh với các nghiên cứu dao động của ô tô tải ba cầu thông thường trong các công trình của Đỗ Tiến Vũ hay Võ Văn Hường, xe chữa cháy rừng đa năng có điểm khác biệt cốt lõi: Téc nước $5,5\text{ m}^3$ làm nâng cao trọng tâm xe từ $h_c = 1,038\text{ m}$ (khi xe không có nước) lên $h_c = 1,18\text{ m}$ (khi chở đầy nước). Việc gia tăng khối lượng được treo lên $13.300\text{ kg}$ làm tăng mô men quán tính toàn xe lên $24.212\text{ kg}\cdot\text{m}^2$, từ đó kéo dài nhẹ chu kỳ dao động riêng nhưng lại mang lại hiệu ứng tích cực là làm giảm tần số dao động xuống vùng tần số êm dịu thích hợp cho cơ thể con người (dưới 120 dao động/phút), giúp người lái duy trì khả năng kiểm soát tay lái chính xác khi len lỏi qua các gốc cây và rãnh sâu.

Kết quả đo thực nghiệm trên hệ thống DMC Plus có sự tương thích cao với mô hình lý thuyết trên Simulink, sai số biên độ gia tốc dao động chỉ dao động trong khoảng 8% đến 12%, khẳng định độ tin cậy vững chắc của mô hình toán 5 bậc tự do đã thiết lập.

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

Dựa trên các phân tích động lực học chuyên sâu, luận văn đề xuất 4 nhóm giải pháp kỹ thuật và tổ chức vận hành nhằm nâng cao độ êm dịu và tính an toàn của phương tiện:

  1. Hiệu chỉnh áp suất lốp theo điều kiện địa hình: Áp dụng hệ thống tự động điều chỉnh áp suất lốp 14.00-20 HC14, giảm áp suất từ mức danh định $3,5\text{ bar}$ xuống dải $2,0 - 2,5\text{ bar}$ khi xe tiến vào đường đất lâm nghiệp gồ ghề. Giải pháp này giúp giảm độ cứng lốp quy đổi từ $1.060.000\text{ N/m}$ xuống xấp xỉ $850.000\text{ N/m}$, tăng diện tích tiếp xúc mặt đất và giảm thêm 18% xung lực va đập truyền lên cầu xe. Thời gian áp dụng: Ngay trong quý 1 giai đoạn bàn giao xe. Chủ thể thực hiện: Kỹ thuật viên và lái xe thuộc lực lượng Kiểm lâm.
  2. Nâng cấp cụm giảm chấn thủy lực cầu trước: Tiến hành thay thế bộ giảm chấn thủy lực nguyên bản của cầu trước bằng loại giảm chấn hai lớp có tích hợp van điều tiết lưu lượng dầu hai chiều, nâng hệ số cản quy đổi từ $k_{n1} = 2.400\text{ N}\cdot\text{s/m}$ lên $2.800\text{ N}\cdot\text{s/m}$. Mục tiêu giảm 15% gia tốc dao động thẳng đứng tại cabin lái trong vòng 6 tháng tới. Chủ thể thực hiện: Nhà máy ô tô Mê Kông phối hợp với nhóm nghiên cứu Trường Đại học Lâm nghiệp.
  3. Thiết lập quy chuẩn dải tốc độ vận hành an toàn: Ban hành quy trình vận hành giới hạn tốc độ tối đa cho xe khi chở đầy 5,5 tấn nước: Duy trì tốc độ $15 - 20\text{ km/h}$ trên các tuyến đường trục lâm nghiệp cấp phối đá, giảm xuống $8 - 12\text{ km/h}$ trên đường đất tự nhiên mấp mô gợn sóng, và giữ dưới $1,5\text{ km/h}$ khi vừa di chuyển vừa kích hoạt đĩa cưa cắt thực bì mở băng cản lửa. Thời gian áp dụng: Quý 2 trước mùa cao điểm khô hanh. Chủ thể thực hiện: Ban Quản lý Rừng phòng hộ và Chi cục Kiểm lâm các tỉnh.
  4. Lắp đặt hệ thống vách ngăn chống sóng chất lỏng trong bồn nước: Bổ sung hệ thống vách ngăn tiêu âm chia bồn chứa $5,5\text{ m}^3$ thành 4 khoang liên thông có van một chiều nhằm triệt tiêu hiện tượng dồn sóng nước tự do (dịch chuyển khối lượng tức thời) khi phanh gấp hoặc leo dốc nghiêng 15 độ, hạn chế phát sinh mô men lắc phụ lên khung xe. Thời gian hoàn thiện: 3 tháng trước khi đưa vào sản xuất hàng loạt. Chủ thể thực hiện: Đơn vị chế tạo cơ khí thuộc Đề tài KC 07.13/06-10.

Đố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 có giá trị học thuật và ứng dụng thực tiễn cao cho 4 nhóm đối tượng sau:

  1. Kỹ sư thiết kế và chế tạo cơ giới hóa nông - lâm nghiệp: Tiếp cận phương pháp lập mô hình động lực học 5 bậc tự do và quy trình tính toán quy đổi tham số đàn hồi, cản nhớt để thiết kế hệ thống treo cho các dòng máy kéo, máy vận xuất gỗ và xe chuyên dụng đặc chủng.
  2. Học viên cao học và nghiên cứu sinh ngành Kỹ thuật ô tô - Cơ khí động lực: Tham khảo chi tiết cách ứng dụng phương trình Lagrange loại II, kỹ thuật lập trình M-file và mô phỏng hệ phương trình vi phân nhiều vật trên công cụ MATLAB/Simulink phục vụ đề tài nghiên cứu dao động phi tuyến.
  3. Cán bộ quản lý kỹ thuật và chỉ huy lực lượng Kiểm lâm, PCCC rừng: Nắm vững các chỉ số kỹ thuật cơ bản (dung tích bồn nước $5,5\text{ m}^3$, chiều dài ống phun $100\text{ m}$, áp lực bơm $150\text{ N/cm}^2$, khả năng vượt dốc $\le 15^\circ$) để tổ chức phương án tác chiến, phân bổ lộ trình cơ động và bảo dưỡng phương tiện đúng quy trình kỹ thuật.
  4. Các doanh nghiệp sản xuất, hoán cải ô tô chuyên dụng: Khai thác các kết quả thực nghiệm đo dao động bằng thiết bị DMC Plus để tối ưu hóa quy trình lắp ráp bồn téc và trang thiết bị chuyên dụng trên các dòng xe cơ sở việt dã 6x6, tiết kiệm từ 40% đến 50% chi phí nghiên cứu phát triển so với nhập khẩu trọn gói từ nước ngoài.

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

1. Tại sao nghiên cứu dao động của xe chữa cháy rừng lại quan trọng hơn các dòng xe tải thông thường?
Xe chữa cháy rừng đa năng phải thường xuyên di chuyển trên địa hình đồi núi không có mặt đường bằng phẳng với tải trọng nước rất lớn ($5,5\text{ m}^3$). Dao động quá mức sẽ làm giảm 40 - 50% tốc độ cơ động, gây mệt mỏi cho tổ lái và dễ phá hủy hệ thống bơm cao áp cùng các chi tiết cơ điện tử tinh vi trên xe.

2. Mô hình động lực học trong luận văn sử dụng bao nhiêu bậc tự do và những tọa độ nào?
Nghiên cứu sử dụng mô hình 5 bậc tự do trong mặt phẳng đối xứng dọc, tương ứng với 5 tọa độ suy rộng: Chuyển dịch thẳng đứng trọng tâm xe ($z$), chuyển vị góc lắc dọc thân xe ($\phi$), chuyển dịch thẳng đứng tâm cầu trước ($z_{11}$), chuyển dịch tâm xoay nhíp sau ($z_c$), và góc xoay đòn gánh nhíp sau ($\theta$).

3. Khối lượng bồn nước 5,5 tấn ảnh hưởng như thế nào đến độ êm dịu của xe?
Khối lượng nước làm tăng tổng khối lượng được treo lên $13.300\text{ kg}$ (đạt tỷ số phân bố khối lượng $\gamma \approx 4,01$ nằm trong khoảng tối ưu từ 4 đến 5). Điều này làm tăng mô men quán tính $J = 24.212\text{ kg}\cdot\text{m}^2$, giúp kéo dài chu kỳ dao động riêng và dập tắt các dao động tần số cao gây hại cho con người.

4. Thiết bị thực nghiệm nào được ứng dụng để kiểm chứng mô hình lý thuyết?
Tác giả sử dụng phương pháp đo các đại lượng không điện bằng điện với thiết bị thu nhận tín hiệu kỹ thuật số đa kênh DMC Plus, phối hợp phần mềm DMC Laplus và các cảm biến gia tốc áp điện gắn trực tiếp lên khung xe theo tiêu chuẩn ngành 22 TCN 336-06.

5. Kết quả nghiên cứu có thể ứng dụng trực tiếp cho các loại xe chữa cháy khác không?
Phương pháp luận xây dựng ma trận khối lượng $M$, ma trận cản $B$, ma trận độ cứng $C$ và thuật toán mô phỏng Simulink trong luận văn hoàn toàn có thể áp dụng mở rộng cho mọi dòng xe tải việt dã 3 cầu (6x6 hoặc 6x4) hoán cải có lắp đặt bồn téc chất lỏng công tác.

Kết luận

  • Luận văn đã xây dựng thành công mô hình toán học dao động 5 bậc tự do trong mặt phẳng thẳng đứng dọc cho mẫu xe chữa cháy rừng đa năng đầu tiên chế tạo tại Việt Nam trên cơ sở khung gầm xe Ural-4320.
  • Thiết lập hệ phương trình vi phân chuyển động bằng phương pháp Lagrange loại II và giải số hoàn chỉnh trên phần mềm MATLAB/Simulink 7.0 dưới các dạng kích động mấp mô mặt đường thực tế.
  • Xác định chính xác các thông số động lực học cốt lõi bằng thực nghiệm và lý thuyết: Khối lượng toàn bộ 16.615 kg, chiều cao trọng tâm chở đầy tải $h_c = 1,18\text{ m}$, mô men quán tính $J = 24.212\text{ kg}\cdot\text{m}^2$, tỷ số phân bố khối lượng đạt ngưỡng chuẩn $\gamma = 4,01$.
  • Kiểm chứng độ tin cậy của mô hình qua thực nghiệm đo gia tốc bằng hệ thống DMC Plus, chứng minh hệ thống treo nhíp trước và đòn gánh cân bằng cầu sau có khả năng dập tắt dao động kích thích chỉ sau chưa đầy 3 giây.
  • Đề xuất gói 4 giải pháp đồng bộ về áp suất lốp, nâng cấp giảm chấn, vách ngăn chống xóc téc nước và giới hạn tốc độ di chuyển từ $15 - 20\text{ km/h}$, đóng góp quan trọng vào việc hoàn thiện sản phẩm Đề tài cấp Nhà nước KC 07.13/06-10.

Bước tiếp theo của hướng nghiên cứu này là mở rộng bài toán dao động không gian 3 chiều có xét đến tác động lắc ngang khi vào cua trên sườn đồi dốc lớn hơn 15 độ và xây dựng mô hình dao động liên hợp máy kéo phụ trợ. Hãy liên hệ ngay với các cơ sở đào tạo chuyên ngành Cơ khí Động lực và Lâm sinh để nhận trọn bộ tài liệu mô phỏng M-file và khai thác ứng dụng thực tiễn cho đơn vị của bạn.