數值模擬在港口安全評價中的應用

顧國強 黃文棟 楊偉明 孫華文 彭修乾

摘要：本文利用數值模擬手段對油罐車爆炸引起的儲油罐破裂隱患進行研究，通過有限元軟件LS-DYNA，采用目前比較通用的炸藥模型參數和計算方法，對空氣中TNT爆炸沖擊波傳播進行了數值模擬，分析典型場景的爆炸影響范圍和油罐車集中停放數量限值，為油罐車停車場的布置和停放提供可量化的建議，提高了油罐車停車場的設計科學性。研究表明在實際行車布置中應重點避免超過十倍特別是二十倍以上單體油罐車的集中爆炸，控制車輛的行駛間距，減少多輛隱患車同時存在一個區域的可能性，并盡量避免具有重大爆炸隱患的車輛停放在儲油罐臨近位置。

Abstract： This paper uses numerical simulation to study the hidden danger of oil tank rupture caused by oil tanker explosion. Through the finite element software LS-DYNA， the current common explosive model parameters and calculation methods are used to carry out TNT explosion shock wave propagation in air. The numerical simulation is carried out to analyze the explosion impact range of typical scenes and the limit number of tank trucks parked， which provides quantifiable suggestions for the layout and parking of tanker parking lots， and improves the scientific design of tanker parking lot. Studies have shown that in the actual driving arrangement， it is important to avoid more than ten times， especially more than twenty times， of concentrated explosion of single tanker trucks， control the driving distance of vehicles， reduce the possibility of multiple vehicles in one area at the same time， and try to avoid Vehicles with major explosion hazards are parked in close proximity to the storage tank.

關鍵詞：爆炸;仿真分析;數值模擬

Key words： explosion;simulation analysis;numerical simulation

中圖分類號：O643.2+21? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?文獻標識碼：A? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 文章編號：1006-4311（2019）25-0285-04

0? 引言

由于爆炸破壞本身的特點，一旦引發爆炸將會造成人員和物體的重大傷害。國內外對于爆炸造成的人員傷害、工程破壞和安全防護方面都進行了相關的探索[1]。

在人員傷害方面，主要是兩個方面：一方面，統計分析爆炸事故中的人員傷亡數據;另一方面，利用動物實驗來等效確定人員傷害。大致研究方向為人員傷害的沖擊波壓強閾值和人員安全距離分級[3]等方面。

在結構破壞方面，一方面利用爆炸實驗進行有限數據的測試，另一方面，針對爆炸實驗條件控制復雜、影響參數眾多等缺點，通過理論分析、有限元、光滑粒子等計算方法來分析爆炸對工程結構的影響。

我國在爆炸對人員傷害和工程結構影響方面，主要以人防結構為研究對象，分析遭受武器打擊時的人員傷害和結構安全抗爆能力，但少有普通民用建筑方面的研究，因此進行民用建筑內的人員安全和結構防爆性能分析具有重大的意義。近年來隨著爆炸事故的增多和安全生產意識的增強，國內相關機構增加了此方面的技術研究投入，但是距離實際應用仍用不足。

目前的不足主要是以下方面：一是分析過于簡化，當量較小，參數變量較少，無法準確分析復雜因素下的爆炸影響;另一方面，受限于分析計算處理設備，大范圍的數值模擬工具無法應用，因而無法開展城市或者大區域級別的爆炸分析。

本文通過理論經驗和數值模擬相結合的手段，針對油罐車爆炸對既有儲油罐體的影響距離進行分析，為油罐車行車路徑和停車位置的爆炸影響進行量化分析，以期能夠減小油罐車爆炸對周邊目標物的影響。

1? 港口安全問題研究

港口的安全問題一直是港口企業管理的重要的內容，隨著安全形勢日益嚴峻，有必要在原來人工經驗管理的基礎上采用更加科學的手段量化分析各種危險源的影響，為科學決策提供參考。

本文利用數值模擬手段對油罐車爆炸引起的儲油罐破裂隱患進行研究，分析典型場景的爆炸影響范圍和油罐車集中停放數量限值，為油罐車停車場的布置和停放位置提供可量化的建議。

主要的分析工況有：

①單輛槽車不同位置爆炸影響范圍確定，計算停車場區域內1、2、3、4四個點位的“單體”爆炸影響范圍。

②同一爆炸點多輛槽車影響范圍計算，計算停車場區域5點位，分別具有十、二十、三十個“單體”的爆炸影響范圍。

③監測點達到相同破壞時不同爆炸點的槽車數量計算，假設爆炸影響到6點位（見附件1），測算1、2、3、4、5點位（見附件1）分別需要的單體數量。

2? 數值計算模型

目前爆炸方面的理論分析，涉及眾多學科，因而只能針對一些簡單爆炸工況或簡化后的復雜工況獲得解析解。

LS-DYNA作為重要的有限元程序，可以處理各種復雜的近場遠場爆炸計算工況，處理沖擊波反射、透射和應力波傳播等問題。

2.1 超壓影響閾值

爆炸沖擊波是爆炸中最主要的破壞形式，衡量沖擊波殺傷力的參數有峰值超壓、沖量和動壓等。在試驗測試基礎上，統計得到的人體和建筑的爆炸傷害破壞超壓值范圍。

2.2 理論公式

理論分析作為爆炸影響分析研究方法的一種，主要針對于規則結構計算物體相應和破壞效果，后續工況計算本文采用經典的Henych空中爆炸沖擊波經驗公式與數值模擬手段進行模型驗證，公式具體形式為：

2.3 計算模型

根據停車場和相關結構實際尺寸以及分析目的，充分考慮爆炸涉及的爆炸點、門衛室及儲油罐等結構，采用殼體和梁來模擬混凝土墻體、柱子、門和窗戶等構件等結構，不同部件相交區域采用共節點方式連接，采用歐拉網格與停車場、門衛室、樓房等結構網格通過流固耦合的處理方式耦合，進行爆炸壓力波對結構的效應分析。

本次計算模型空氣域采用400*400*400m的規模，炸藥當量體積約為14.3m3，采用六面體單元對塊體部件進行網格劃分，采用梁單元對梁柱結構進行網格劃分，采用殼單元對薄壁結構進行網格劃分，并進行網格穩定和敏感性測試，最終選擇0.5m作為炸藥和空氣單元網格基本尺寸，劃分的模型網格數量為147.8萬，節點數量151.6萬。

2.4 載荷和邊界設置

2.4.1 載荷的處理

爆炸載荷通過將換算的當量體積單元轉換成TNT材料屬性，并將爆炸點設置在TNT單元的中心處。

2.4.2 約束的處理

爆炸計算模型的約束主要有以下五個方面：

①空氣域邊界設置無反射邊界;②建筑結構柱子底部設置約束;③炸藥單元與空氣單元共節點;④空氣與結構單元交界處設置耦合;⑤樓房結構的樓板和柱子對應節點進行節點合并。

2.5 材料參數

其中炸藥TNT和玻璃采用的材料參數如表4和表5所示。

3? 計算結果和分析

利用大型通用有限元軟件LS-DYNA，模擬包括TNT炸藥起爆、爆轟波傳播、沖擊波與物體相互作用及結構相應等過程。

3.1 單個油罐車爆炸范圍計算

本工況為單個車體當量TNT爆炸影響范圍計算，影響范圍分為對人和對車輛建筑兩個大方向，圖4給出了單個儲油罐當量超壓峰值隨距離的變化。

根據2.1節以及相關文獻[2][4]，人員傷害的死亡半徑、重傷半徑和輕傷半徑對應的超壓沖擊波數值接近但不完全相同，本文為方便計算初步采用爆炸沖擊波死亡半徑超壓90kPa、重傷超壓44kPa、輕傷超壓17kPa來進行影響范圍計算。通過對圖4中距離和超壓的關系插值，得到上述三個影響半徑值，如表6所示。

對物體影響方面，表7列舉不同破壞程度對應的超壓范圍值，為方便計算對比，取影響范圍下限值進行影響半徑的計算，根據爆炸點距離儲油罐的距離，點位①和③沒有影響，但是②和④的破壞范圍會覆蓋儲油罐。

3.2 不同數量油罐車影響

表8和表9給出了⑤點位不同數量油罐車爆炸產生的沖擊波超壓峰值，根據停車場和周圍建筑物的空間關系，單體油罐車爆炸主要造成二層樓房、停車場爆炸側車輛和此區域內容人員傷亡;十倍單體油罐車爆炸主要造成南北門衛室以東停車場區域內車輛破壞、近停車場一側的儲油罐毀壞和此區域內人員死亡;二十倍單體油罐車爆炸主要造成過半停車場內車輛破壞和人員死亡、所有儲油罐破壞;三十倍單體油罐車爆炸主要造成絕大部分停車場內車輛破壞和人員死亡、所有儲油罐毀壞。

根據上述計算結果，應重點避免超過十倍特別是二十倍以上單體油罐車的集中爆炸，實際行車中需控制車輛的行駛間距，減少多輛車同時存在一個區域的可能性。

3.3 儲油罐破壞所需爆炸車輛

表10給出了引起儲油罐破裂所需的不同點位爆炸車輛數量，①點位距離儲油罐268m，需要4個單體罐車同時爆炸才能引起儲油罐破壞;②與④點位距離儲油罐距離分別為99m和53m，處于單體罐車破壞范圍之內;③號點位距離最近儲油罐230m，需要2個單體罐車同時爆炸才能引起儲油罐破壞;⑤號點位距離儲油罐256m，需要3個單體罐車同時爆炸才能引起儲油罐破壞。

根據上述結果，實際停放過程中盡量避免具有相應數量的重大爆炸隱患車輛停放在儲油罐臨近位置。

4? 結束語

綜上所述，在實際行車布置中應重點避免超過十倍特別是二十倍以上單體油罐車的集中爆炸，控制車輛的行駛間距，減少多輛隱患車同時存在一個區域的可能性，并盡量避免具有重大爆炸隱患的車輛停放在儲油罐臨近位置。

本文利用數值模擬手段對油罐車爆炸引起的儲油罐破裂隱患進行研究，通過有限元軟件LS-DYNA并采用目前比較通用的炸藥模型參數和計算方法，對空氣中TNT爆炸沖擊波傳播進行了數值模擬，分析典型場景的爆炸影響范圍和油罐車集中停放數量限值，為油罐車停車場的布置和停放提供可量化的建議，提高了油罐車停車場的設計科學性。

參考文獻：

[1]亨利奇.爆炸動力學及其應用[M].熊建國，譯. 北京：科學出版社，1987.

[2]孫艷馥.爆炸沖擊波對人體損傷和防護分析[J].火炸藥學報，31（4）.

[3]李錚.空氣沖擊波作用下的安全距離[J].爆炸與沖擊，10（2）.

[4]王文和.天然氣管道泄漏火球事故后果模擬評價[J].中國安全生產科學技術，8（1）.

[5]葉曉華.軍事爆破工程[M].北京：解放軍出版社，1999，8.

[6]韋灼彬.鋼筋混凝土樁基梁板碼頭爆炸毀傷及搶修技術研究[D].天津：天津大學，2005，6.

[7]章冠人.凝聚炸藥起爆動力學[M].北京：國防工業出版社， 1991.