水下和空中爆炸沖擊波傳播特性對比分析

張社榮，孔 源，王高輝,2

(1.天津大學(xué) 水利工程仿真與安全國家重點實驗室，天津 300072; 2.武漢大學(xué) 水資源與水電工程科學(xué)國家重點實驗室，武漢 430072)

隨著精確制導(dǎo)武器的快速發(fā)展以及國內(nèi)外恐怖襲擊和意外爆炸事件不斷發(fā)生，對重大建筑物結(jié)構(gòu)的安全構(gòu)成了巨大的威脅。建筑物結(jié)構(gòu)遭受爆炸沖擊荷載的主要來源有空中爆炸和水下爆炸。如船艦結(jié)構(gòu)和水工大壩結(jié)構(gòu)可能遭受來自制導(dǎo)炸彈的空中爆炸、魚雷的水下爆炸等沖擊荷載作用。然而由于水和空氣兩種介質(zhì)的物理屬性存在較大差異，且爆炸產(chǎn)物與炸彈周圍介質(zhì)的相互作用效應(yīng)不同，沖擊波在水和空氣中的傳播特性存在較大的差異，導(dǎo)致其對結(jié)構(gòu)的損傷機制、破壞程度及防護設(shè)計均有所不同。當(dāng)炸彈在近自由面水下或空中爆炸時，由于近自由面反射產(chǎn)生的稀疏波或沖擊波與入射沖擊波相互作用，自由面界面效應(yīng)對沖擊波的傳播特性產(chǎn)生了較大影響，將導(dǎo)致近自由面爆炸的荷載特征與自由場爆炸不同。因此研究沖擊波在不同介質(zhì)中的傳播特性及界面效應(yīng)對其傳播規(guī)律的影響，對結(jié)構(gòu)的抗爆防護設(shè)計具有重要意義。

爆炸是一種多學(xué)科交叉耦合的物理現(xiàn)象，由于它的復(fù)雜性，早期研究主要以理論和試驗研究為主[1-2]。近年來，隨著計算機技術(shù)及實驗手段的不斷進步，與爆炸相關(guān)的理論、實驗和數(shù)值模擬研究快速發(fā)展，使得數(shù)值模擬爆炸成為可能。在沖擊波對比分析方面，如寧心等[3]采用試驗的方法對比研究了水下和空中爆炸沖擊波傳播速度和物理參數(shù)；Rajendran等[4]分別對空中和水下爆炸沖擊波傳播的特性進行了理論分析；黃建松等[5]從理論上對比分析了水下和空中爆炸對船員的沖擊損傷效應(yīng)；李順波等[6]分析了爆炸沖擊波在水、土、混凝土中的衰減規(guī)律。在近自由面爆炸方面，Zamyshlyaev[7]采用試驗方法研究了深水爆炸情況下的自由面對沖擊波的影響；張鵬翔等[8]探討了淺層水中爆炸水底、水面對沖擊波的切斷現(xiàn)象，崔杰等[9]基于無網(wǎng)格SPH方法研究了近自由面水下爆炸切斷現(xiàn)象的特性及產(chǎn)生機理，Xie等[10]采用MGFM方法研究了近自由水面水下爆炸產(chǎn)生氣穴現(xiàn)象的特性和機理。以上研究主要側(cè)重從試驗和理論上分析自由場水下和空中沖擊波傳播的物理屬性，而關(guān)于界面效應(yīng)對沖擊波傳播特性的影響則主要集中在近自由面水下爆炸的水面切斷效應(yīng)方面。

本文以顯式動力分析程序AUTODYN為平臺，通過建立自由場水下和空中爆炸耦合數(shù)值模型，采用Euler方法，對比分析爆炸沖擊波在水下和空中的傳播特性。同時考慮沖擊波與自由面反射的稀疏波相互作用過程，模擬了近自由面水下和空中爆炸過程，研究了近自由面對水下和空中爆炸沖擊波傳播特性的影響。研究成果可為自由場和近自由面水下和空中爆炸荷載的確定提供參考，并為結(jié)構(gòu)的抗爆防護設(shè)計提供基礎(chǔ)。

1 材料模型及狀態(tài)方程

1.1 炸藥

在爆轟性能的考察和計算中，炸藥爆轟產(chǎn)物的狀態(tài)方程必不可少，在爆轟高溫高壓條件下，知道某時刻產(chǎn)物的組成和各組分的熱力學(xué)參數(shù)，即可按照某種混合法則建立爆轟產(chǎn)物的狀態(tài)方程。高能炸藥材料采用JWL狀態(tài)方程[11]描述了爆轟壓力P和每單位體積內(nèi)能E及相對體積V的關(guān)系：

(1)

式中:P為爆轟壓力；V為爆轟產(chǎn)物體積和炸藥初始體積之比；e為炸藥的初始內(nèi)能；A、B、R1、R2、ω為特征參數(shù)，各參數(shù)取值如下：A=3.738×1011Pa，B=3.747×1011Pa，R1=4.15，R2=0.9，ω=0.35，D=6.93×103m/s，e=5.999×109J/m3，PCJ=2.10×1010Pa，ρ0=1.63×103kg/m3。

1.2 水

炸藥在水介質(zhì)中爆炸時，在裝藥內(nèi)形成高溫高壓的爆轟產(chǎn)物，其壓力遠遠大于周圍水介質(zhì)的靜壓力，在水介質(zhì)中會產(chǎn)生水中沖擊波和氣泡脈動現(xiàn)象。水介質(zhì)采用多項式狀態(tài)方程[11]進行描述，其形式由根據(jù)壓縮狀態(tài)的不同而定。

當(dāng)水為壓縮狀態(tài)(μ>0)時，其狀態(tài)方程為：

P=A1μ+A2μ2+A3μ3+(B0+B1μ)ρ0e

(2)

當(dāng)水為膨脹狀態(tài)(μ<0)時，其狀態(tài)方程為：

P=T1μ+T2μ2+B0ρ0e

(3)

式中:μ=ρ/ρ0-1，ρ0是初始密度；A1、A2、A3、B0、B1、T1和T2是由AUTODYN材料庫直接賦值的常數(shù)；e是比內(nèi)能，其定義式如下：

e=(ρgh+p0)/ρB0

(4)

其中:ρ和h分別是水的密度和深度，g和p0分別是重力加速度和大氣壓強。各參數(shù)取值如下：ρ0=1×103kg/m3，A1=2.2×109Pa，A2=9.54×109Pa，A3=1.457×1010Pa，B0=0.28，B1=0.28，T1=2.2×109Pa，T2=0 Pa。

1.3 空氣

炸藥在空氣中爆炸時，瞬時轉(zhuǎn)變?yōu)楦邷馗邏旱谋óa(chǎn)物。爆炸產(chǎn)物在空氣中進行膨脹，其結(jié)果是在爆炸產(chǎn)物內(nèi)形成稀疏波。同時，爆炸產(chǎn)物強烈壓縮空氣，在空氣中形成爆炸空氣沖擊波。空氣采用Mat-Null材料模型和理想氣體狀態(tài)方程[11]：

(5)

式中:E是比內(nèi)能，γ是絕熱指數(shù)(取1.4)，ρ0是空氣初始密度(取1.225 kg/m3)，ρ是當(dāng)前密度。

2 爆炸沖擊波傳播特性對比分析

2.1 沖擊波經(jīng)驗公式

水下爆炸中，峰值超壓往往比靜水壓力大了幾個數(shù)量級，因此對于淺水爆炸，靜水壓力可忽略不計，峰值超壓也被簡化稱作峰值壓力。沖擊波的壓力時程曲線在某一點由瞬時壓力迅速攀升到峰值壓力，接著又呈指數(shù)形式衰減，其沖擊波壓力時程曲線、峰值壓力及沖量根據(jù)Cole經(jīng)驗公式[12]有：

P(t)=Pme-t/θ

(6)

(7)

(8)

(9)

式中:Pm為沖擊波峰值壓力，Pa；θ為沖擊波指數(shù)衰減時間常數(shù)；I為壓力沖量，N·s/m2；W為TNT炸藥質(zhì)量，kg；R為測點到爆心的距離，m；k，α，l，β是與炸藥性能有關(guān)的經(jīng)驗參數(shù)，對于TNT炸藥取值如下：k=5.33×106Pa，α=1.13，l=5 880，β=0.89。

空中爆炸時，峰值超壓在爆心距較大處比較小，因此大氣壓不可忽略。本文采用經(jīng)典的Henrych·J空中爆炸沖擊波經(jīng)驗公式[13]與數(shù)值模擬進行比較。下面給出具體形式：

ΔP(t)=ΔPf(1-t/tp)e-α0t/tp

(10)

tp=0.001(0.107+0.444Z+

0.05≤Z≤3

(12)

(13)

式中:ΔPf為沖擊波峰值超壓，MPa；Z是比例距離，m/kg1/3；R為爆心到測點的距離，m；W為TNT炸藥質(zhì)量，kg；tp為沖擊波超壓持續(xù)時間，ms；t為正壓作用時間，ms；α0為由試驗決定的衰減系數(shù)。

沖擊波沖量由Kinney的經(jīng)驗公式計算得到[14]，其具體形式如下，單位為kN·s/m2：

(14)

2.2 不同爆炸方式下沖擊波傳播特性對比分析

為了研究不同爆炸方式下沖擊波傳播特性，本文通過建立二維軸對稱自由場水下和空中爆炸有限元數(shù)值模型，對自由場水下和空中爆炸沖擊波傳播過程進行模擬。其計算模型如圖1所示，計算區(qū)域為15 m×20 m，計算網(wǎng)格尺寸為25 mm。球形TNT裝藥量為900 kg，起爆點位于炸藥的中心。水、空氣和炸藥均采用Euler算法(Euler算法適合于描述液體和氣體的行為，自由邊界面和材料的交界面通過固定的Euler網(wǎng)格來表達，故大變形或者有流動的情形不會導(dǎo)致網(wǎng)格畸變)。為了模擬自由場水下和空中爆炸，在模型截斷邊界處施加無反射邊界條件(Flow_out邊界條件)，使得沖擊波在人工截斷邊界處無反射。炸藥在水下爆炸時，將主要產(chǎn)生沖擊波壓力和氣泡脈動壓力，而文中主要對比分析沖擊波在水下和空中兩種介質(zhì)中的傳播特性，故在水下爆炸時暫不考慮氣泡脈動壓力的影響。

圖2給出了爆心距6 m處的水下和空中爆炸沖擊波典型壓力時程曲線。由圖2(a)可知，當(dāng)水下爆炸沖擊波傳播到目標(biāo)處時，沖擊波壓力迅速到達峰值，約為96.7 MPa，隨后沖擊波壓力呈指數(shù)衰減。由圖2(b)可知，當(dāng)炸藥在自由場空中爆炸時，沖擊波壓力時程曲線正壓階段與水下爆炸沖擊波傳播規(guī)律相同，超壓峰值為2.1 MPa，比水下爆炸該點處的峰值壓力小很多；沖擊波壓力由峰值衰減到標(biāo)準(zhǔn)大氣壓后，繼續(xù)衰減，即出現(xiàn)負(fù)壓區(qū)。相對正壓階段，負(fù)壓階段壓力較小，在簡化空中爆炸沖擊波荷載時，往往可忽略負(fù)壓區(qū)，只考慮正壓超壓作用。

圖1 自由場水下/空中爆炸計算模型

圖2 爆心距6 m處的壓力時程曲線

圖3給出了水下和空中爆炸沖擊波傳播的峰值壓力和沖量對比。在近爆區(qū)域，由于試驗測量誤差較大，數(shù)值模擬得到的水下和空中爆炸沖擊波峰值壓力和沖量與經(jīng)驗公式得到的峰值壓力和沖量存在一定差異。而隨著爆心距的增大，數(shù)值模擬得到數(shù)值與采用經(jīng)驗公式得到的數(shù)值較接近，說明了數(shù)值模擬技術(shù)的可靠性。由圖3(a)可知，在距炸藥中心相同距離處，水下爆炸沖擊波峰值壓力遠大于空中爆炸，在爆心距1 m處，水下爆炸沖擊波峰值壓力達1 661 MPa，而空中爆炸時，沖擊波峰值壓力衰減至88 MPa，沖擊波壓力在空氣中衰減較快。水下爆炸沖擊波峰值平均壓力是空中爆炸下的55.7倍(經(jīng)驗平均比值為60.4)，主要由于水的密度大，壓縮性較小(通常認(rèn)為是不可壓縮的流體)，而空氣密度小，可壓縮性大，沖擊波在空氣中傳播時能量大量且快速地耗散在空氣中。由圖3(b)可知，在距炸藥中心相同距離處，水下爆炸沖擊波沖量遠大于空中爆炸，在爆心距1 m處，水下爆炸沖擊波沖量達324 kPa·s，而空中爆炸時，沖擊波壓力沖量已衰減至7 kPa·s，沖擊波壓力在空氣中衰減較快。水下爆炸沖擊波平均沖量是空中爆炸下的66.5倍(經(jīng)驗平均比值為62.3)。

圖3 水下和空中爆炸沖擊波傳播的峰值壓力和沖量對比

圖4 不同爆炸方式下爆炸沖擊波到達時間對比圖

圖4給出了水下和空中爆炸沖擊波傳播到不同爆心距處的時間對比。由圖4可知，在近爆區(qū)域，空中沖擊波傳播速度較水下快，主要由于空氣可壓縮性較大，沖擊波在空中傳播的初始階段幾乎是以炸藥的爆轟速度傳播，而水的可壓縮性較小，沖擊波傳播速度在水下衰減較快，使得近爆區(qū)域沖擊波在空中傳播速度比在水下傳播快。而隨著爆心距的增加，由于空中爆炸能量耗散較快，沖擊波傳播速度很快衰減至空中聲速，故水下沖擊波傳播速度較空中快。

圖5 近自由面水下/空中爆炸計算模型

3 沖擊波傳播界面效應(yīng)分析

為了研究近自由面對水下和空中爆炸沖擊波傳播特性的影響，通過建立二維軸對稱近自由面爆炸有限元數(shù)值耦合模型，對近自由面水下和空中爆炸沖擊波傳播過程進行模擬。其計算模型如圖5所示，計算區(qū)域為20 m×40 m，計算網(wǎng)格尺寸為25 mm。球形TNT裝藥量為900 kg，起爆點位于炸藥的中心。水、空氣和炸藥均采用Euler算法。在模型截斷邊界處施加無反射邊界條件(Flow_out邊界條件)，以消除邊界的沖擊波反射。研究近自由面水下爆炸時，A區(qū)域為空氣，B區(qū)域為水，炸藥起爆深度為10 m；研究近自由面空中爆炸時，A區(qū)域為水，B區(qū)域為空氣。

圖6 近自由面水下爆炸模型的壓力云圖

圖7 近自由面空中爆炸模型的壓力云圖

圖6和圖7分別給出了近自由面水下和空中爆炸沖擊波壓力傳播過程，圖8給出了自由面附近區(qū)域的沖擊波峰值壓力變化過程曲線。由圖6可知，當(dāng)炸藥在水下起爆5.2ms后，沖擊波傳播到自由水面處(圖6(a))，并與自由水面相互作用，將在自由水面處產(chǎn)生向水下反射的強烈稀疏波以及向空中傳播的透射沖擊波。由于空氣的聲阻抗較水要小很多，傳播到空氣中的透射沖擊波強度比反射稀疏波弱很多，沖擊波峰值壓力由水下的47 MPa降低到空中的0.03 MPa(圖8(a))，故在圖6(b)～(d)中幾乎觀察不到空氣中透射波的存在。而反射的稀疏波強度較大，與入射沖擊波相互作用使得靠近自由水面的水壓力迅速降低，將在自由水面下方產(chǎn)生氣穴現(xiàn)象，如圖6(b)所示(t=6 ms)。隨著稀疏波的繼續(xù)傳播，自由水面下的氣穴將越來越大，如圖6(c)～(d)所示。當(dāng)炸藥在近自由面空中爆炸時，沖擊波在5.5 ms傳播到自由面處(圖7(a))，并與自由面相互作用，將在自由面處產(chǎn)生向空中反射的沖擊波以及向水下傳播的透射沖擊波。反射的沖擊波與入射沖擊波相互作用，使界面處出現(xiàn)沖擊波增強效應(yīng)。傳播到水下的透射沖擊波強度比反射沖擊波強很多，峰值壓力由空中的0.8 MPa上升到水下的5.7 MPa(圖8(b))。

圖9和圖10分別給出了近自由面水下和空中爆炸下距自由面不同距離(該距離記為d)處的壓力時程變化曲線。由圖9可知，當(dāng)水下爆炸沖擊波在自由面反射形成的稀疏波到達后，測點的壓力迅速下降，不能按照自由場爆炸的壓力衰減規(guī)律繼續(xù)衰減，波形似被稀疏波切斷，稱為自由水面沖擊波切斷效應(yīng)。它使測點處的壓力受到削弱，因而減弱了沖擊波對外做功的能力。距自由水面距離越近，沖擊波削弱越多。而當(dāng)炸藥在近自由面空中爆炸時，自由面反射的沖擊波到達后，測點的壓力迅速上升，出現(xiàn)沖擊波增強效應(yīng)，亦使得測點處的壓力未按自由場爆炸的壓力衰減規(guī)律繼續(xù)衰減，如圖10所示。

圖8 近自由面沖擊波峰值壓力變化過程曲線

4 結(jié) 論

本文基于顯式動力分析程序AUTODYN，建立了自由場和近自由面水下和空中爆炸的耦合數(shù)值模型，對比分析自由場水下和空中爆炸沖擊波傳播特性，并考慮沖擊波與自由面反射的稀疏波相互作用過程，研究了近自由面對水下和空中爆炸沖擊波傳播特性的影響主要得到以下結(jié)論：

圖9 近自由面和自由場水下爆炸壓力時程對比曲線

圖10 近自由面和自由場空中爆炸壓力時程對比曲線

(1) 爆炸沖擊波在水和空氣兩種介質(zhì)中的傳播特性存在較大差異。水下爆炸時，沖擊波壓力瞬時上升到峰值隨后指數(shù)衰減，與空中爆炸沖擊波的正相超壓區(qū)間特征相似；水下爆炸沖擊波峰值壓力和沖量均比空中爆炸大很多，且衰減慢，對結(jié)構(gòu)的潛在破壞能力較強，對可能遭受水下和空中兩種爆炸方式作用的結(jié)構(gòu)，應(yīng)重點關(guān)注水下爆炸下的抗爆設(shè)防。

(2) 自由界面對沖擊波傳播特性存在較大的影響。對于近自由面水下爆炸，自由面反射的稀疏波與入射沖擊波相互作用使得靠近自由水面的水壓力迅速降低，波形似被稀疏波切斷，并在自由水面下方發(fā)生氣穴現(xiàn)象。而在近自由面空中爆炸出現(xiàn)沖擊波增強效應(yīng)。

[1] Arons A B. Underwater explosion shock wave parameters at large distances from the charge [J]. The Journal of the Acoustical Society of America, 1954, 26: 343.

[2] 高勇軍, 王偉策，陳曉波. 淺層水中爆炸沖擊波壓力的測試與分析[J]. 爆破, 1999, 16(1): 9-13.

GAO Yong-jun, WANG Wei-ce, CHEN Xiao-bo. Testing and anslysis on shockwave pressure generated by explosion under shallow water [J].Blasting, 1999,16(1):9-13.

[3] 寧心, 李曉炎, 楊志煥,等. 水下沖擊波和空氣沖擊波傳播速度及物理參數(shù)的對比研究[J]. 解放軍醫(yī)學(xué)雜志, 2004, 29(2): 97-99.

NING Xin, LI Xiao-yan,YANG Zhi-huan, et al. A comparative study on the propagation speed and physical parameters of underwater blast wave and air blast wave [J]. Med J Chin PLA, 2004, 29(2): 97-99.

[4] Rajendran R, Lee J M. Blast loaded plates [J]. Marine Structures, 2009, 22(2): 99-127.

[5] 黃建松, 汪玉. 水下和空中爆炸艦員沖擊損傷效應(yīng)分析[J]. 噪聲與振動控制, 2012, 32(6): 13-16.

HUANG Jian-song, WANG Yu. Analysis of shock injury of ship personnel due to underwater and air explosion [J]. Noise and Vibration Control, 2012, 32(6): 13-16.

[6] 李順波, 東兆星, 齊燕軍,等. 爆炸沖擊波在不同介質(zhì)中傳播衰減規(guī)律的數(shù)值模擬[J]. 振動與沖擊, 2009, 28(7): 115-117.

LI Shun-bo, DONG Zhao-xing, QI Yan-jun, et al. Numerical simulation for spread decay of blasting shock wave in different media [J].Journal of Vibration and Shock, 2009, 28(7): 115-117.

[7] Zamyshlyayev B V. Dynamic loads in underwater explosion,AD=757183 [R]. 1973.

[8] 張鵬翔, 顧文彬, 葉序雙. 淺層水中爆炸沖擊波切斷現(xiàn)象淺探[J]. 爆炸與沖擊. 2002,22(3):221-228.

ZHANG Peng-xiang, GU Wen-bin, YE Xu-shuang. Discussions of blasting shock waves cutoff in shallow-layer water [J]. Explosion and Shock Waves, 2002,22(3):221-228.

[9] 崔杰, 楊文山, 李世銘,等. 近自由面水下爆炸沖擊波切斷效應(yīng)研究[J]. 船舶力學(xué), 2012, 16(5): 465-471.

CUI Jie, YANG Wen-shan, LI Shi-ming, et al. Research on the cutoff effect of shock Wave induced by underwater explosion near free surface [J]. Journal of Ship Mechanics. 2012, 16(5): 465-471.

[10] Xie W F, Liu T G, Khoo B C. The simulation of cavitating flows induced by underwater shock and free surface interaction [J]. Applied Numerical Mathematics, 2007, 57(5): 734-745.

[11] ANSYS Inc. AUTODYN user manual version 13; 2010.

[12] Cole R H. Underwater explosions [M]. New York: Dover Publications, 1965.

[13] Henrych J, Major R. The dynamics of explosion and its use [M]. New York: Elsevier Scientific Publishing company, 1979.

[14] Kinney G F, Graham K J. Explosive shocks in air[M]. Berlin and New York, Springer-Verlag, 1985.