Tổng quan nghiên cứu

Trong ngành vận tải hàng hóa chuyên dùng, tai nạn giao thông liên quan đến phương tiện chở chất lỏng nguy hiểm luôn tiềm ẩn rủi ro thảm khốc về người, tài sản và môi trường. Theo các báo cáo thống kê an toàn giao thông tại Bắc Mỹ, tai nạn liên quan đến xe tải hạng nặng chiếm khoảng 28% tổng số vụ tai nạn nghiêm trọng, trong đó các dòng xe xi-téc sơmi rơ-moóc chiếm tới 62,67% số vụ lật xe và mất ổn định hướng khi phanh hoặc chuyển làn. Hiện tượng sóng sánh của chất lỏng bên trong bồn chứa khi chưa điền đầy là nguyên nhân cốt lõi tạo ra các lực và mômen thủy động lực học phi tuyến, làm dịch chuyển tức thời tọa độ trọng tâm và gây đảo lộn phân bố tải trọng giữa các trục bánh xe.

Luận văn thạc sĩ chuyên ngành Kỹ thuật Cơ khí Động lực của tác giả Lê Trung Hiếu tại Trường Đại học Bách Khoa – Đại học Quốc gia Thành phố Hồ Chí Minh tập trung giải quyết bài toán: Nghiên cứu ảnh hưởng chuyển động của xăng bên trong bồn nhiên liệu đến phân bố tải trọng của sơmi rơ-moóc chở xăng khi phanh. Mục tiêu cụ thể của nghiên cứu là xây dựng mô hình mô phỏng số ba chiều chất lượng cao nhằm khảo sát sự dao động sóng sánh của 40 mét khối xăng (tương đương khối lượng hàng hóa 29.600 kg) trong bồn chứa hình trụ elip có kích thước tiết diện 2.500 x 1.960 mm và chiều dài 10.600 mm.

Nghiên cứu được giới hạn khảo sát trong điều kiện phanh giảm tốc thẳng đều với gia tốc 3 m/s2 và 5 m/s2 từ các dải vận tốc ban đầu 40 km/h, 60 km/h và 80 km/h, dưới các mức nạp chất lỏng 50%, 70% và 90%. Kết quả của công trình mang ý nghĩa khoa học và thực tiễn sâu sắc, cung cấp cơ sở định lượng chuẩn xác giúp các nhà sản xuất xe chuyên dùng tối ưu hóa kết cấu vách ngăn, qua đó giảm thiểu từ 18% đến hơn 30% xung lực va đập tác dụng lên chốt kéo và cụm cầu sau, nâng cao rõ rệt hệ số an toàn khi vận hành sơmi rơ-moóc trên mạng lưới giao thông đường bộ.

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 cơ học lưu chất đa pha và động lực học đoàn xe sơmi rơ-moóc chuyên dụng:

  • Lý thuyết dòng chảy đa pha và bề mặt tự do: Sử dụng hệ phương trình vi phân đạo hàm riêng Navier-Stokes kết hợp phương trình liên tục cho dòng lưu chất ba chiều không nén được và phụ thuộc thời gian. Để theo dõi sự biến dạng phức tạp của bề mặt phân chia giữa hai pha xăng và không khí trong quá trình sóng sánh biên độ lớn, mô hình thể tích chất lỏng Volume of Fluid được thiết lập với điều kiện tổng tỷ lệ thể tích trong từng ô lưới tính toán luôn bằng 1.
  • Lý thuyết mô hình dòng chảy rối: Áp dụng mô hình rối k - epsilon tiêu chuẩn và mô hình Shear Stress Transport hai phương trình để giải quyết các xoáy cuộn năng lượng và độ nhớt rối sinh ra khi dòng nhiên liệu va chạm mạnh vào các bề mặt vách ngăn.
  • Động lực học phân bố tải trọng đoàn xe: Hệ thống các công thức cân bằng lực tĩnh và động xác định trọng tâm sơmi rơ-moóc G43-BX40-02 có chiều dài cơ sở 7.930 mm. Chiều cao trọng tâm đoàn xe thay đổi từ 1.630 mm ở trạng thái không tải (tổng khối lượng 9.415 kg) lên 2.240 mm ở trạng thái đầy tải (tổng khối lượng 39.015 kg), với khoảng cách trọng tâm đến tâm chốt kéo là 4.880 mm và đến tâm trục sau là 3.050 mm.

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

Nghiên cứu sử dụng phương pháp tính toán động lực học chất lỏng số hóa Computational Fluid Dynamics dựa trên nền tảng phương pháp thể tích hữu hạn Finite Volume Method trên phần mềm Ansys Fluent kết hợp môi trường lập trình mở rộng:

  • Nguồn dữ liệu và cỡ mẫu mô phỏng: Cỡ mẫu nghiên cứu gồm 18 kịch bản mô phỏng động học toàn diện được xây dựng từ ma trận kết hợp giữa 3 dải vận tốc (40, 60, 80 km/h), 2 cấp gia tốc phanh (3, 5 m/s2) và 3 mức nạp nhiên liệu (50%, 70%, 90% thể tích bồn). Các cấu hình vách ngăn gồm: mô hình cơ bản không vách chắn sóng, mô hình 5 vách ngăn kết hợp 1 vách chắn sóng, và mô hình 5 vách ngăn kết hợp 2 vách chắn sóng với các chiều cao vách chắn sóng được khảo sát ở mức 30%, 40% và 50% chiều cao bồn.
  • Kỹ thuật chia lưới và chọn mẫu: Sử dụng mạng lưới thể tích không cấu trúc với mật độ làm mịn cao tại các vùng tiếp giáp thành bồn và mép vách ngăn elip nhằm triệt tiêu hiện tượng khuếch tán số và đảm bảo sự hội tụ nghiệm chính xác.
  • Rationale lựa chọn phương pháp: Phương pháp thể tích hữu hạn kết hợp mô hình VOF được lựa chọn vì khả năng giải quyết chính xác các bài toán cơ học chất lỏng phi tuyến có biên tự do chuyển động mạnh, giúp tiết kiệm hơn 80% chi phí và thời gian so với việc thực hiện các thử nghiệm thực địa nguy hiểm với chất lỏng dễ cháy nổ.

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

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

Phân tích dữ liệu từ 18 kịch bản mô phỏng đã làm sáng tỏ các quy luật biến thiên tải trọng cơ bản của sơmi rơ-moóc chở xăng trong quá trình phanh:

  • Mức nạp nhiên liệu quyết định mức độ sóng sánh: Mức điền đầy 50% thể tích bồn tạo ra biên độ sóng sánh dữ dội nhất và gây biến thiên tải trọng lớn nhất lên hệ thống liên kết. Khi mức nạp tăng lên 70%, hiện tượng dao động giảm dần và đạt trạng thái ổn định cơ học cao nhất ở mức nạp 90%, do không gian tự do cho chất lỏng dồn dịch bị thu hẹp đáng kể.
  • Tác động cộng hưởng của vận tốc và gia tốc phanh: Khi gia tốc phanh tăng từ 3 m/s2 lên 5 m/s2 ở cùng vận tốc 80 km/h, lực xô đẩy dọc trục của 29.600 kg xăng dồn về phía trước làm tải trọng tức thời tác dụng lên chốt kéo tăng vọt hơn 25% đến 35% so với giá trị tĩnh 15.015 kg. Đồng thời, tải trọng phân bố lên cụm 3 trục sau bị sụt giảm nghiêm trọng so với tải tĩnh 24.000 kg, làm suy giảm lực bám mặt đường và tăng nguy cơ trượt lết bánh xe.
  • Hiệu quả triệt tiêu xung lực của vách chắn sóng: Bồn chứa có cấu hình 5 vách ngăn kết hợp 1 hoặc 2 vách chắn sóng phụ giúp giảm biên độ dao động tải trọng lên chốt kéo từ 18% đến 28% so với bồn không có vách chắn sóng. Đáng chú ý, việc gia tăng chiều cao vách chắn sóng từ 30% lên 50% chiều cao bồn đem lại hiệu quả ngăn chặn sóng sánh vượt trội hơn hẳn so với việc chỉ gia tăng số lượng vách chắn sóng.

Thảo luận kết quả

Khi xe thực hiện quá trình phanh khẩn cấp, khối chất lỏng dịch chuyển tự do về phía trước tạo thành sóng va đập thủy động cực mạnh tác dụng lên nắp đáy trước và các vách ngăn. Khi sử dụng vách chắn sóng có chiều cao đạt 50% (tương đương 1.250 mm), dòng xăng bị cưỡng bức đổi hướng và tạo thành các vùng xoáy cục bộ, giúp tiêu tán động năng sóng sánh nhanh hơn gấp 2 lần so với vách cao 30%.

Kết quả này hoàn toàn tương thích với các kết luận trong công trình nghiên cứu của Ranganathan và Kang về sự dịch chuyển trọng tâm chất lỏng khi phanh, đồng thời phát triển chi tiết hơn nhờ mô hình hóa 3D trực quan. Trên các đồ thị biểu diễn biến thiên tải trọng theo thời gian, đường đặc tính phản lực tại chốt kéo xuất hiện dạng sóng hình sin tắt dần; cấu hình vách chắn sóng cao 50% giúp rút ngắn thời gian dập tắt dao động hoàn toàn xuống chỉ còn khoảng 1,5 đến 2,0 giây sau khi xe dừng hẳn.

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

  • Chuẩn hóa kết cấu vách ngăn xi-téc: Các doanh nghiệp sản xuất và lắp ráp rơ-moóc chuyên dùng cần áp dụng bắt buộc cấu hình tối thiểu 5 vách ngăn chính kết hợp hệ thống vách chắn sóng có chiều cao đạt từ 40% đến 50% chiều cao tiết diện elip (khoảng 1.000 đến 1.250 mm). Mục tiêu giảm tối thiểu 25% biến thiên tải trọng động lên chốt kéo, hoàn thành áp dụng vào quy trình thiết kế kỹ thuật trong vòng 6 đến 12 tháng.
  • Quy định dải tải trọng vận hành an toàn: Doanh nghiệp vận tải xăng dầu cần ban hành quy chế vận hành nội bộ, nghiêm cấm việc lưu thông bồn chở xăng ở mức nạp lơ lửng nguy hiểm từ 40% đến 60% dung tích thiết kế; ưu tiên phương án vận chuyển ở mức nạp trên 85% thể tích hoặc điều chuyển rỗng dưới 10% thể tích, thực hiện ngay lập tức trên toàn bộ đội xe.
  • Tích hợp hệ thống phân bổ lực phanh điện tử EBS thông minh: Các đơn vị nghiên cứu và chế tạo cần phối hợp thử nghiệm lắp đặt hệ thống van điều khiển lực phanh thích ứng theo tải trọng động tức thời, bù trừ cho sự dịch chuyển chiều cao trọng tâm từ 1.630 mm lên 2.240 mm khi phanh, lộ trình thử nghiệm trong 18 tháng.
  • Kiểm soát thông số vận hành và đào tạo lái xe: Đơn vị quản lý vận tải cần đào tạo kỹ năng phanh giảm tốc êm dịu cho tài xế, duy trì gia tốc phanh dưới 3 m/s2 và kiểm soát vận tốc xe bồn dưới 60 km/h trên các cung đường đèo dốc hoặc trơn trượt để hạn chế dao động thủy động.

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

  • Kỹ sư thiết kế cơ khí và phương tiện ô tô chuyên dùng: Tiếp cận phương pháp tính toán tọa độ trọng tâm, phân bố tải tĩnh (chốt kéo 15.015 kg, trục sau 24.000 kg) và ứng dụng phần mềm mô phỏng số để tối ưu hóa kết cấu thép SS400 dày 4 mm trên xi-téc chở xăng.
  • Doanh nghiệp kinh doanh và logistics vận tải xăng dầu: Nắm bắt các vùng rủi ro động lực học theo mức nạp 50%, 70%, 90% để xây dựng phương án xếp dỡ hàng hóa và quy trình chạy xe an toàn trên các tuyến quốc lộ.
  • Giảng viên, học viên cao học và sinh viên ngành Cơ khí Động lực: Khai thác quy trình ứng dụng công cụ CFD Fluent, phương pháp thể tích hữu hạn và mô hình VOF làm tài liệu học thuật mẫu mực cho các đề tài tương tác chất lỏng - kết cấu.
  • Cán bộ đăng kiểm và cơ quan quản lý tiêu chuẩn giao thông: Sử dụng các số liệu thực nghiệm mô phỏng về tải trọng trục và gia tốc phanh 3 đến 5 m/s2 để hoàn thiện các quy chuẩn kỹ thuật quốc gia về an toàn cho phương tiện chở hàng nguy hiểm.

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

  • Tại sao mức điền đầy 50% lại gây mất ổn định lớn nhất cho xe bồn khi phanh?
    Ở mức nạp 50%, bề mặt tự do của chất lỏng đạt diện tích tiếp xúc lớn nhất trong bồn elip 2.500 x 1.960 mm, tạo không gian chuyển dịch tối đa cho 14.800 kg xăng dồn về phía trước, sinh ra mômen quán tính và lực va đập thủy động cực đại làm đảo lộn phân bố tải trọng.

  • Tăng chiều cao vách chắn sóng có ưu thế gì so với việc tăng số lượng vách?
    Tăng chiều cao vách lên 40% hoặc 50% giúp chặn đứng dòng chảy tràn bề mặt của lớp chất lỏng phía trên, triệt tiêu động năng sóng sánh hiệu quả hơn và giảm đáng kể tổng khối lượng bản thân bồn so với việc bổ sung quá nhiều vách ngăn mỏng (tiết kiệm hơn 15% khối lượng vật liệu thép).

  • Trọng tâm của sơmi rơ-moóc G43-BX40-02 thay đổi như thế nào giữa hai trạng thái tải?
    Khi không tải (9.415 kg), chiều cao trọng tâm là 1.630 mm và cách tâm chốt kéo 6.100 mm. Khi chở đầy 40 mét khối xăng (tổng tải 39.015 kg), chiều cao trọng tâm tăng lên 2.240 mm và tọa độ dọc dịch chuyển về phía trước, cách chốt kéo 4.880 mm.

  • Mô hình thể tích chất lỏng VOF giải quyết bài toán sóng sánh xăng dầu như thế nào?
    Mô hình VOF chia miền tính toán thành lưới hữu hạn, theo dõi tỷ lệ thể tích từng pha lỏng và khí trong mỗi ô kiểm tra, giúp mô tả chính xác chuyển động tách sóng, xoáy cuộn và áp suất phân bố lên thành bồn mà không cần dựng lại lưới chuyển động.

  • Gia tốc phanh ảnh hưởng như thế nào đến tải trọng tác dụng lên chốt kéo?
    Gia tốc phanh tăng từ 3 m/s2 lên 5 m/s2 làm lực quán tính dọc trục tăng tỷ lệ thuận, khiến phản lực đè lên chốt kéo tăng thêm từ 2.000 kg đến hơn 4.500 kg so với trạng thái tĩnh, làm nhẹ cầu sau và giảm khả năng dẫn hướng của đầu kéo.

Kết luận

  • Luận văn đã mô hình hóa thành công hiện tượng sóng sánh ba chiều của 29.600 kg xăng trong sơmi rơ-moóc G43-BX40-02 khi phanh bằng phương pháp thể tích hữu hạn VOF.
  • Nghiên cứu chỉ ra rằng mức nạp 50% thể tích là trạng thái nguy hiểm nhất, trong khi mức nạp 90% đạt độ ổn định dọc tối ưu nhất khi phanh.
  • Cấu hình 5 vách ngăn kết hợp vách chắn sóng cao 40% đến 50% chiều cao bồn được chứng minh là giải pháp tối ưu giúp giảm hơn 25% biến thiên tải trọng lên chốt kéo và trục sau.
  • Kết quả mô phỏng định lượng chính xác sự gia tăng áp lực động khi vận tốc phanh tăng từ 40 km/h lên 80 km/h và gia tốc phanh tăng từ 3 m/s2 lên 5 m/s2.
  • Công trình đặt nền tảng khoa học vững chắc cho việc thiết kế, cải tiến thùng xi-téc chuyên dụng và xây dựng quy chuẩn an toàn giao thông đường bộ tại Việt Nam.

Hướng phát triển tiếp theo của đề tài trong giai đoạn 12 đến 24 tháng tới là mở rộng mô phỏng bài toán động học phi tuyến trong trường hợp xe vừa phanh vừa chuyển hướng quay vòng trên đường nghiêng. Các kỹ sư thiết kế, nhà khoa học và doanh nghiệp vận tải hãy ứng dụng ngay các thông số kết cấu và khuyến nghị vận hành trong luận văn để nâng cao chất lượng phương tiện và bảo đảm an toàn giao thông tối đa.