Tổng quan nghiên cứu

Các sự cố cháy nổ công nghiệp và nguy cơ mất an toàn tại các khu chứa hóa chất dễ bắt cháy đang là thách thức nghiêm trọng trên phạm vi toàn cầu, khi áp suất sóng nổ có thể đạt mức tức thời lên đến 300 kbar cùng nhiệt độ từ 3000°C đến 4000°C. Sóng xung kích phát sinh lan truyền trong không khí với vận tốc ban đầu xấp xỉ 7200 m/s, tạo ra xung lực cực đại tàn phá các công trình dân dụng liền kề. Trước thực trạng đó, luận văn thạc sĩ kỹ thuật xây dựng thực hiện năm 2020 tại Trường Đại học Kinh tế Công nghiệp Long An đã tập trung giải quyết bài toán mô phỏng sự tương tác phi tuyến giữa tải trọng nổ và cụm kết cấu khu dân cư.

Mục tiêu cốt lõi của nghiên cứu là xây dựng quy trình phân tích động lực học phi tuyến, xác định trường phân bố áp suất thoáng qua, đo lường biến thiên vận tốc và mức tiêu tán năng lượng của sóng xung kích qua thời gian thực từ 5 ms đến 6 ms. Phạm vi nghiên cứu khảo sát mô hình cụm 7 khối nhà liền kề có chiều cao từ 300 mm đến 450 mm trong không gian ba chiều chịu tác động trực tiếp từ khối thuốc nổ TNT chuẩn. Ý nghĩa khoa học và thực tiễn của đề tài thể hiện qua việc giảm thiểu hơn 90% chi phí thử nghiệm thực địa nguy hiểm, đồng thời đạt độ chuẩn xác cao với sai số kiểm chứng dưới 5% so với các công bố khoa học quốc tế, cung cấp cơ sở định lượng quan trọng cho công tác quy hoạch vùng đệm an toàn và gia cường công trình đô thị.

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 trên nền tảng cơ học môi trường liên tục phi tuyến và khí động lực học sóng xung kích cường độ cao. Hai khung lý thuyết nòng cốt bao gồm: lý thuyết đàn dẻo phi tuyến kết hợp mô hình biến cứng tốc độ biến dạng Cowper-Symonds để mô tả ứng xử kết cấu bê tông cốt thép dưới tải trọng động, và lý thuyết nhiệt động học sóng nổ Rankine-Hugoniot liên hệ các thông số gián đoạn trước và sau mặt sóng xung kích.

Phương trình trạng thái Jones-Wilkins-Lee (JWL) được ứng dụng chuẩn xác để mô hình hóa quá trình giải phóng năng lượng của khối nổ TNT với mật độ khối 1630 kg/m3, hệ số áp lực quy đổi 3,734.10^6 kPa, vận tốc kích nổ 6930 m/s và mật độ năng lượng đạt 6,0.10^6 kJ/m3. Đi cùng với đó là các định luật bảo toàn khối lượng, động lượng, năng lượng và quy luật tỷ lệ khoảng cách nổ Hopkinson-Cranz (Z = R / W^(1/3)) kết hợp cùng công thức Brode và Mills nhằm định lượng chính xác quá trình suy giảm của áp suất dư và áp lực phản xạ cực đại trên bề mặt vật cản.

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

Nghiên cứu ứng dụng phương pháp phần tử hữu hạn Arbitrary Lagrangian-Eulerian (ALE) trong môi trường phần mềm tính toán chuyên dụng ANSYS Autodyn và Explicit Dynamics. Cỡ mẫu mô hình số được phân chia thành hai miền tính toán với quy mô lưới chặt chẽ: miền lưu chất khí và khối nổ TNT (Eulerian) bao gồm 352.944 nút và 336.175 phần tử; miền công trình xây dựng (Lagrangian) sử dụng phần tử bậc cao SOLID186 20 nút gồm 17.055 nút và 12.694 phần tử.

Quy trình chọn mẫu và kiểm soát chất lượng lưới được kiểm chứng nghiêm ngặt qua 4 tiêu chuẩn hình học: độ nghiêng lệch Skewness trung bình đạt 0,20187 (tối đa 0,70483), chất lượng phần tử Element Quality trung bình đạt 0,88999, tỷ lệ cạnh Aspect Ratio trung bình 1,3596 và hệ số Jacobian trung bình đạt 1,0769. Phương pháp thu thập dữ liệu áp dụng kỹ thuật bố trí đầu đo áp lực điện tử (gauge points) tại 13 tọa độ chiến lược (từ T1 đến T13) trên các mặt đón sóng, mặt hông và vùng khuất của các khối công trình. Lý do lựa chọn phương pháp ALE là nhằm kết hợp hoàn hảo khả năng mô tả dòng khí nén vận tốc siêu âm của lưới Eulerian mà không gây méo mó phần tử, đồng thời lưới Lagrangian mô tả chính xác ứng xử biến dạng và chuyển vị của kết cấu công trình.

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 động lực học trong khoảng thời gian từ 0 đến 5 ms và 6 ms đã mang lại những phát hiện quan trọng về cơ chế lan truyền tải nổ trong đô thị:

Thứ nhất, sóng nổ hình thành xung áp suất cực đại tức thời và suy giảm phi tuyến tính cực nhanh. Tại điểm khảo sát T1 (tọa độ x = 300 mm, y = 1100 mm, z = 150 mm), áp suất sóng xung kích vọt lên mức đỉnh chỉ trong 0,8 ms đầu tiên trước khi tiêu tán dần về áp suất khí quyển xung quanh mức 100 kN/m2.

Thứ hai, cấu trúc hình học đô thị tạo ra sự che chắn và phân tán áp lực rõ rệt. Mặt đón sóng trực diện của khối nhà A và B chịu áp suất phản xạ cao hơn từ 38% đến 45% so với điểm T3 (tọa độ x = 1260 mm, y = 1000 mm, z = 400 mm) nằm tại khối nhà G ở phía sau cụm công trình.

Thứ ba, vận tốc truyền sóng nổ giảm mạnh từ mức cực đại 6930 m/s tại tâm nổ xuống dưới 1200 m/s sau khi vượt qua khoảng cách 1500 mm và va đập qua 2 dãy công trình chắn sóng đầu tiên.

Thứ tư, mật độ năng lượng của sóng nổ hấp thụ vào bề mặt kết cấu suy giảm theo quy luật bậc ba theo khoảng cách, trong đó các công trình ngoại vi hấp thụ khoảng 65% tổng động năng ban đầu của luồng khí nén.

Thảo luận kết quả

Hiện tượng gia tăng áp suất đột biến tại các mặt trước của tòa nhà được giải thích do sự phản xạ sóng trực diện, làm vận tốc hạt không khí bị triệt tiêu đột ngột và chuyển hóa toàn bộ động năng thành thế năng áp suất cục bộ. Khi sóng xung kích luồn qua các khe hẹp giữa các khối nhà cao 300 mm đến 450 mm, hiệu ứng co thắt dòng khí (hiệu ứng hẹp) đã khiến vận tốc sóng tại các ngõ ngách tăng cục bộ thêm 15% trước khi suy giảm.

Dữ liệu tính toán từ 13 điểm đo được đối chiếu trực tiếp với công bố quốc tế của các chuyên gia trong ngành, cho thấy dạng đường cong áp suất theo thời gian và giá trị áp suất đỉnh có độ trùng khớp trên 95,5%. Trong các báo cáo kỹ thuật, toàn bộ tập dữ liệu này được trực quan hóa tối ưu thông qua các biểu đồ phân bố áp suất theo thời gian (Overpressure vs Time), đồ thị so sánh vận tốc lan truyền và bảng tổng hợp xung lượng sóng phản xạ tại từng cao độ mặt đứng, giúp người đọc dễ dàng nhận diện vùng nguy hiểm trọng yếu.

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

Nhằm nâng cao năng lực kháng chấn và chống chịu tải nổ cho các công trình khu dân cư lân cận khu công nghiệp, 4 nhóm giải pháp thực thi bao gồm:

  1. Thiết lập khoảng lùi an toàn và quy hoạch vùng đệm sinh thái: Các cơ quan quản lý đô thị cần ban hành quy chuẩn khoảng cách ly tối thiểu từ 50 m đến 100 m giữa tường rào khu công nghiệp chứa hóa chất nguy hiểm và khu dân cư, giúp giảm thiểu hơn 60% áp lực sóng xung kích tác động trực tiếp lên kết cấu công trình theo lộ trình thực hiện 12 tháng.

  2. Gia cường kết cấu chịu lực mặt trước công trình: Các chủ đầu tư và đơn vị thi công cần áp dụng giải pháp bọc sợi composite polymer (FRP) hoặc tăng cường hàm lượng cốt thép đai tại các cột tầng 1 (từ cao độ 0 m đến 3,5 m), nâng cao khả năng hấp thụ năng lượng va đập thêm 30% đến 40% trong thời gian hoàn thiện 6 tháng cho mỗi dự án.

  3. Xây dựng tường chắn sóng tiêu tán năng lượng: Ban quản lý các khu công nghiệp cần lắp đặt hệ thống tường chắn bê tông xốp hoặc tường nghiêng chuyển hướng sóng nổ xung quanh các bồn chứa nhiên liệu, giúp triệt tiêu ít nhất 25% động năng luồng khí nén trong vòng 18 tháng triển khai.

  4. Ứng dụng quy trình mô phỏng số ALE trong thẩm duyệt thiết kế: Các sở xây dựng và đơn vị tư vấn thiết kế cần tích hợp bắt buộc bước phân tích động lực học phi tuyến trên các phần mềm như ANSYS Autodyn cho các công trình công cộng tập trung trên 500 người trước năm 2027.

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

Nội dung và dữ liệu thực nghiệm số trong luận văn mang lại giá trị thiết thực cho 4 nhóm đối tượng chính:

  • Kỹ sư thiết kế kết cấu công trình: Khai thác quy trình thiết lập tải trọng nổ tương đương, phương pháp tính toán ma trận độ cứng phi tuyến và các biểu thức sóng xung kích để áp dụng vào tính toán kết cấu kháng nổ cho các tòa nhà văn phòng, chung cư cao tầng.

  • Học viên cao học và nhà nghiên cứu chuyên ngành Kỹ thuật Xây dựng: Sử dụng luận văn như tài liệu tham khảo chuẩn mực về phương pháp ghép nối tương tác chất lưu - kết cấu (ALE), cách khai báo phương trình trạng thái JWL và quy trình chia lưới chất lượng cao SOLID186 trong ANSYS Workbench.

  • Cơ quan quản lý quy hoạch đô thị và Phòng cháy chữa cháy: Ứng dụng kết quả phân tích trường áp suất và vận tốc lan truyền để xây dựng bản đồ cảnh báo rủi ro, quy định bán kính an toàn phòng nổ và phê duyệt quy hoạch khu dân cư giáp ranh khu chế xuất.

  • Ban quản lý dự án và chủ đầu tư hạ tầng khu công nghiệp: Tham khảo các thông số kỹ thuật về mức độ suy giảm tải nổ để bố trí hợp lý các phân khu chức năng, nhà kho hóa chất và hệ thống công trình phụ trợ giảm chấn.

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

  1. Phương pháp Arbitrary Lagrangian-Eulerian (ALE) giải quyết hạn chế nào của phương pháp phần tử hữu hạn truyền thống? Phương pháp Lagrangian truyền thống thường gặp lỗi dừng tính toán do phần tử lưới bị méo mó nghiêm trọng khi chất lưu biến dạng lớn. Phương pháp ALE cho phép lưới Eulerian mô phỏng dòng khí nén đạt vận tốc 6930 m/s chuyển động linh hoạt qua các phần tử cố định, đồng thời lưới Lagrangian theo dõi chính xác biến dạng của tường và cột nhà.

  2. Phương trình trạng thái Jones-Wilkins-Lee (JWL) đóng vai trò gì trong bài toán mô phỏng nổ? Phương trình JWL xác định chính xác áp suất khí sinh ra trong quá trình giãn nở của sản phẩm nổ TNT từ mật độ 1630 kg/m3 và năng lượng 6,0.10^6 kJ/m3. Đây là cầu nối vật lý biến đổi năng lượng hóa học thành trường áp suất động học tác động trực tiếp lên môi trường không khí xung quanh.

  3. Tiêu chí nào đảm bảo độ tin cậy của mô hình lưới với hơn 350.000 phần tử trong nghiên cứu? Nghiên cứu kiểm soát độ tin cậy thông qua chỉ số chất lượng phần tử trung bình đạt 0,88999 và độ nghiêng lệch Skewness chỉ 0,20187. Lưới được làm mịn cục bộ quanh vị trí tâm nổ và các góc cạnh công trình, đảm bảo gradient áp suất được ghi nhận liên tục và không gây phân kỳ toán học.

  4. Tại sao áp suất phản xạ trên mặt tường công trình lại lớn hơn nhiều lần so với áp suất sóng tới trong không khí tự do? Khi mặt sóng xung kích va chạm vuông góc vào bề mặt cứng của tòa nhà, luồng khí chuyển động tốc độ cao bị chặn lại tức thời. Quá trình hãm động năng này tạo ra một sóng phản xạ ngược chiều, làm mật độ phân tử khí nén tăng đột biến và đẩy áp suất phản xạ lên gấp 2 đến 8 lần áp suất sóng tới.

  5. Kết quả mô phỏng trên mô hình thu nhỏ có thể quy đổi sang quy mô công trình thực tế như thế nào? Dữ liệu mô phỏng hoàn toàn có thể tịnh tiến sang công trình thực tế thông qua định luật đồng dạng nổ Hopkinson-Cranz với tham số khoảng cách tỷ lệ Z. Khi khối lượng thuốc nổ TNT tăng lên W lần, khoảng cách an toàn và kích thước hình học sẽ tỷ lệ thuận với căn bậc ba của W trong khi áp suất cực đại không đổi.

Kết luận

  • Luận văn đã xây dựng thành công quy trình mô phỏng 3D hiện đại đánh giá tác động của sóng nổ lên khu dân cư bằng phương pháp ALE trong ANSYS Autodyn với độ chính xác cao.
  • Kiểm chứng chặt chẽ dữ liệu trên 352.944 nút lưới miền khí và 17.055 nút lưới kết cấu, đạt sai số dưới 4,5% so với các công trình nghiên cứu quốc tế tiêu chuẩn.
  • Phân tích chi tiết quy luật tiêu tán áp suất và suy giảm vận tốc sóng nổ từ 6930 m/s xuống mức an toàn khi lan truyền qua các dãy nhà có chiều cao từ 300 mm đến 450 mm.
  • Cung cấp 4 giải pháp kỹ thuật có tính ứng dụng cao về khoảng lùi an toàn 50 m đến 100 m, gia cường mặt đứng công trình và quy chuẩn hóa thẩm duyệt thiết kế kháng chấn.
  • Hướng phát triển tiếp theo trong 24 tháng tới là mở rộng bài toán sang mô hình phá hủy rạn nứt phi tuyến của vật liệu bê tông cốt thép thực tế và tương tác lan truyền trong địa hình đô thị phức tạp.

Các kỹ sư kết cấu, nhà quy hoạch và nhà nghiên cứu hãy ứng dụng ngay mô hình tính toán số tiên tiến này vào thực tiễn dự án để tối ưu hóa an toàn công trình đô thị!