侯楓琳 張建橋 陳鵬 馬藝琳 楊旭 郝莉

【摘 要】采用多流體網(wǎng)格法，基于JWL狀態(tài)方程，用面向?qū)ο蟮腃++語言自行開發(fā)編制了NM—MMIC通用多物質(zhì)二維流體彈塑性程序，對有和無防護(hù)檔墻的爆炸空氣沖擊波進(jìn)行了數(shù)值模擬，并用VISC 2D可視化軟件對計算結(jié)果進(jìn)行動畫演示，從中不僅可以看到爆炸過程的整體概貌，還可以清晰地顯示細(xì)微過程。計算結(jié)果基本符合沖擊波饒過擋墻的物理現(xiàn)象和規(guī)律，這說明該物理模型和數(shù)值算法較為合理，可用于對可能發(fā)生的爆炸事故進(jìn)行風(fēng)險評估。

【關(guān)鍵詞】沖擊波；多流體網(wǎng)格法；數(shù)值模擬

0.引言

在現(xiàn)代生活中，爆炸災(zāi)害的發(fā)生與日俱增，造成建筑物和設(shè)備的破壞及人員的傷亡。大多數(shù)的工業(yè)災(zāi)害事故是由爆炸引起的。由爆炸產(chǎn)生的空氣沖擊波遇到目標(biāo)，如建筑物、軍事設(shè)施等會發(fā)生反射和繞射現(xiàn)象，研究這些現(xiàn)象對于摧毀敵方的軍事目標(biāo)和加強(qiáng)我方的防御工事，以及對危險品生產(chǎn)車間的防護(hù)都有很大的實用價值。本文采用多流體的Euler型算法，運用自行編制的NM-MMIC通用程序，模擬了在爆點周圍有無防護(hù)擋墻的二維爆炸場，并通過考慮防護(hù)擋墻的形狀對其后爆炸場的影響，給出了一些對工程設(shè)計有參考價值的結(jié)果。

1.計算方法

1.1總體方程組

爆炸沖擊波的形成過程可用二維不定常彈塑性流體力學(xué)方程組描述[1]。在沖擊波的數(shù)值模擬中，涉及到三種介質(zhì)：炸藥、土和空氣。炸藥采用TNT， 爆轟產(chǎn)物的初始條件采用瞬時爆轟模型，狀態(tài)方程分別為：對于空氣，近似采用理想氣體狀態(tài)方程，對炸藥采用JWL狀態(tài)方程，土介質(zhì)所采用的狀態(tài)方程是：

p=8.626×10ee （1）

式中，P為壓力， εpv為體積應(yīng)變。

圖1帶防護(hù)墻的初始圖

1.2數(shù)值計算方法

計算域（z×r）取150×150m，網(wǎng)格劃分（z×r）為600×600個。網(wǎng)格步長取0.25 。邊界條件為：左邊界取對稱邊界，其他三邊取連續(xù)邊界條件，計算采用顯式差分格式。應(yīng)力，密度，比內(nèi)能，速度，應(yīng)力偏量等均定義在網(wǎng)格中心，人工粘性定義在網(wǎng)格邊界中點[2]。對于含兩種物質(zhì)的混合格，運用Youngs界面處理技術(shù)，當(dāng)混合網(wǎng)格中有三種物質(zhì)時，采用簡單線界面技術(shù)（SLIC）。圖1為正梯形擋墻的計算模型。

2.數(shù)值模擬結(jié)果

本文研究了無擋墻、矩形擋墻、正梯形擋墻、左斜右直擋墻等不同形狀和位置的防護(hù)擋墻對爆炸空氣沖擊波的影響，考察了圓柱形區(qū)域的高壓氣體對地面和防護(hù)擋墻的沖擊波流場的發(fā)展規(guī)律[3]。圖2為正梯形擋墻的爆炸過程。

圖2 正梯形擋墻的爆炸過程

通過計算可以得出如下結(jié)論：

（1）靠近爆心一側(cè)的防護(hù)擋墻形狀對遠(yuǎn)場壓力的影響較大，梯形擋墻要好于矩形擋墻；而遠(yuǎn)離爆心一側(cè)的防護(hù)擋墻形狀對遠(yuǎn)場壓力的影響不明顯，梯形擋墻與左斜右直的擋墻對遠(yuǎn)場的效果差別不大。

（2）無防護(hù)擋墻的超壓要比有防護(hù)擋墻時大得多，這表明防護(hù)擋墻在減弱遠(yuǎn)場壓力方面還是起到了一定作用的。

（3）在擋墻后方高度的2倍距離處的峰值壓力比無擋墻時要大得多，這是因為在此處形成馬赫反射的緣故。這一點在設(shè)計防護(hù)擋墻時應(yīng)特別注意。

3.結(jié)論

本文針對爆炸沖擊波問題，采用多流體網(wǎng)格法，研究了防護(hù)擋墻對爆炸效應(yīng)的影響。針對不同位置、不同形狀的擋墻，給出了爆炸場的各個物理量隨時間的變化關(guān)系。在此基礎(chǔ)上重點討論了遠(yuǎn)場壓力對防護(hù)擋墻的依賴關(guān)系，為工程設(shè)計提供了十分有用的數(shù)值分析手段。

【參考文獻(xiàn)】

[1]寧建國，張文耀，鄭祖國，李偉.爆炸與沖擊數(shù)值模擬的可視化.兵工學(xué)報，第21卷增刊，2000，8：89-91.

[2]寧建國，王仲琦，趙衡陽等.沖擊波繞流的數(shù)值模擬研究[J].北京理工大學(xué)學(xué)報，1999，19（5）.

[3]郝莉.多物質(zhì)二維流體彈塑性程序的理論基礎(chǔ)及其應(yīng)用研究，北京理工大學(xué)博士學(xué)位論文，2003，8.