晏衛(wèi)東，馬曉明，劉延青，張 戈，魏 琳，郭 君，張阿漫

（1.中國(guó)人民解放軍某部，遼寧 葫蘆島 125000；2.哈爾濱工程大學(xué)，黑龍江 哈爾濱 150001）

0 引言

反艦導(dǎo)彈武器技術(shù)的發(fā)展，使得艦船在作戰(zhàn)環(huán)境中越來(lái)越容易受到反艦導(dǎo)彈武器的攻擊。大部分反艦導(dǎo)彈武器會(huì)在艦船結(jié)構(gòu)附近或艙室內(nèi)發(fā)生爆炸，沖擊波超壓峰值極高，會(huì)對(duì)艦船結(jié)構(gòu)產(chǎn)生巨大的破壞作用。因此，研究反艦導(dǎo)彈武器空中爆炸載荷特性及其對(duì)船體結(jié)構(gòu)毀傷機(jī)理對(duì)于提高艦船結(jié)構(gòu)在戰(zhàn)斗環(huán)境中的生命力有很重要的意義。國(guó)內(nèi)外研究人員對(duì)艦船結(jié)構(gòu)在爆炸載荷作用下的動(dòng)態(tài)響應(yīng)進(jìn)行了理論和試驗(yàn)研究，研究結(jié)果大多表明近場(chǎng)爆炸載荷會(huì)對(duì)艦船局部結(jié)構(gòu)產(chǎn)生“終點(diǎn)效應(yīng)”，即艦船結(jié)構(gòu)受到攻擊后發(fā)生最大變形以及毀滅性破壞［1，2］。此類(lèi)研究特點(diǎn)是研究成果都是建立在特定的艦船結(jié)構(gòu)以及特定載荷基礎(chǔ)之上，在應(yīng)用范圍和使用條件方面存在較大差異，在實(shí)際應(yīng)用過(guò)程中會(huì)產(chǎn)生不同程度的誤差，并且公式計(jì)算的結(jié)果并不能說(shuō)明結(jié)構(gòu)的實(shí)時(shí)動(dòng)態(tài)響應(yīng)。本文應(yīng)用應(yīng)力波及塑性動(dòng)力學(xué)理論計(jì)算了船體板架結(jié)構(gòu)在近場(chǎng)空爆載荷作用下的實(shí)時(shí)變形，得到了一些有實(shí)際工程意義的結(jié)論。

1 無(wú)限介質(zhì)近場(chǎng)爆炸載荷特性

炸藥在空中爆炸時(shí)，瞬時(shí)轉(zhuǎn)變?yōu)楦邷馗邏旱谋óa(chǎn)物，爆炸產(chǎn)物在空氣中膨脹，目標(biāo)離爆心越近，所受沖擊力越大，但是如果所受作用面積較小，盡管沖擊波壓力很高，破壞作用也只帶有局部性。

2 無(wú)限介質(zhì)中爆炸載荷驗(yàn)證及正壓區(qū)作用時(shí)間確定

2.1 載荷驗(yàn)證

由于空氣近場(chǎng)爆炸載荷的破壞性極大，因而公開(kāi)發(fā)表近空氣自由場(chǎng)爆炸載荷的數(shù)據(jù)比較少，本文采用數(shù)值計(jì)算方法對(duì)爆炸載荷計(jì)算的準(zhǔn)確性進(jìn)行驗(yàn)證。應(yīng)用多物理場(chǎng)非線(xiàn)性有限元軟件ANSYS/LS-DYNA對(duì)TNT炸藥在無(wú)限介質(zhì)場(chǎng)中爆炸進(jìn)行數(shù)值模擬計(jì)算，模擬計(jì)算的流場(chǎng)和藥包的有限元模型如圖1所示。

圖1 流場(chǎng)及藥包有限元模型圖

根據(jù)實(shí)際仿真計(jì)算結(jié)果，經(jīng)過(guò)試算、研究、分析，得到戰(zhàn)斗部炸藥和無(wú)限介質(zhì)場(chǎng)的仿真模型和計(jì)算參數(shù)。

仿真計(jì)算的有限元網(wǎng)格劃分質(zhì)量直接影響到結(jié)果精度，本文采用逐層控制劃分法對(duì)所需的無(wú)限介質(zhì)場(chǎng)進(jìn)行網(wǎng)格劃分，以提高計(jì)算結(jié)果準(zhǔn)確度。通過(guò)大量的仿真計(jì)算，得到戰(zhàn)斗部炸藥在無(wú)限介質(zhì)場(chǎng)中爆炸時(shí)的沖擊波傳播規(guī)律和距爆心不同距離的無(wú)限介質(zhì)場(chǎng)的壓力時(shí)間曲線(xiàn)，如圖2～圖4所示。

從圖2～圖4可以看出，爆炸沖擊波以球形波的方式由爆心向外在流場(chǎng)中傳播，與爆心距離越大的位置，沖擊波壓力峰值越小，正壓作用時(shí)間越大，此特點(diǎn)符合自由空氣場(chǎng)中爆炸沖擊波在空氣中的傳播規(guī)律。表1為不同距離無(wú)限介質(zhì)場(chǎng)壓力峰值經(jīng)驗(yàn)公式計(jì)算值與數(shù)值計(jì)算值比較。

圖2 不同時(shí)刻爆炸沖擊波在空氣中的傳播

由表1可知，空氣中近場(chǎng)爆炸產(chǎn)物在距離藥包較近的位置產(chǎn)生極大的壓力，初始階段隨著爆炸氣體的體積膨脹，壓力值衰減很快，隨著膨脹體積的增大，壓力衰減速度比初始階段變慢，且為逐漸衰減到零的趨勢(shì)。數(shù)值計(jì)算結(jié)果的特點(diǎn)與炸藥空中爆炸物理過(guò)程定性符合，且經(jīng)驗(yàn)公式計(jì)算值與數(shù)值模擬計(jì)算值誤差約為10%左右，說(shuō)明了計(jì)算方法的有效性。

圖3 距爆心不同單元自由場(chǎng)壓力時(shí)歷曲線(xiàn)

圖4 不同距離自由場(chǎng)壓力時(shí)歷曲線(xiàn)

表1 不同距離無(wú)限介質(zhì)場(chǎng)壓力峰值經(jīng)驗(yàn)公式計(jì)算值與數(shù)值模擬計(jì)算值比較

2.2 正壓區(qū)作用時(shí)間確定

應(yīng)用經(jīng)驗(yàn)公式確定正壓區(qū)作用時(shí)間，如表2所示。由于其適用范圍的限制，具有較大的計(jì)算誤差。

表2 不同距離正壓區(qū)作用時(shí)間經(jīng)驗(yàn)公式計(jì)算值與數(shù)值模擬計(jì)算值比較

本文對(duì)數(shù)值計(jì)算結(jié)果采用數(shù)據(jù)擬合方法，得到正壓區(qū)作用時(shí)間的回歸計(jì)算公式，具體沖擊波衰減規(guī)律計(jì)算公式參見(jiàn)文獻(xiàn)［3］。

3 爆炸載荷作用下板架結(jié)構(gòu)的變形

3.1 作用在板架結(jié)構(gòu)的爆炸載荷計(jì)算

作用在板架結(jié)構(gòu)的爆炸載荷計(jì)算與無(wú)限自由場(chǎng)爆炸載荷計(jì)算不同，當(dāng)爆炸載荷作用在板架結(jié)構(gòu)表面時(shí)，形成復(fù)雜的反射作用，具體計(jì)算方法采用文獻(xiàn)［4，5］提供的方法。

3.2 板架初始破口的計(jì)算

由于載荷強(qiáng)度高、作用時(shí)間短，加上板架結(jié)構(gòu)面積較小，可將艦船的板架結(jié)構(gòu)近似看做加筋結(jié)構(gòu)與板具有相同破壞模式，將加筋結(jié)構(gòu)均攤到面板上采用等效板厚方法來(lái)處理。

3.3 板架變形計(jì)算

為了在簡(jiǎn)化問(wèn)題的同時(shí)得到工程實(shí)際問(wèn)題可以接受的結(jié)果，板架破口計(jì)算時(shí)，將艦船板架結(jié)構(gòu)簡(jiǎn)化為軸對(duì)稱(chēng)問(wèn)題，假設(shè)所關(guān)心區(qū)域?yàn)閳A形，這樣一來(lái)問(wèn)題就簡(jiǎn)化為處理圓板問(wèn)題［6，7］。

根據(jù)文獻(xiàn)［7］研究?jī)?nèi)容，推導(dǎo)出正壓作用時(shí)間內(nèi)板架破口計(jì)算方法。板架撓度的計(jì)算分為塑性變形區(qū)的撓度計(jì)算和破口區(qū)的撓度計(jì)算。

4 算例驗(yàn)證

目前，空中近場(chǎng)接觸爆炸實(shí)驗(yàn)數(shù)據(jù)資源不足，可參考比對(duì)的數(shù)據(jù)較少，破口半徑經(jīng)驗(yàn)公式值引用吉田?。?］根據(jù)二戰(zhàn)艦船破損資料的試驗(yàn)結(jié)果。本文采用上述方法，應(yīng)用數(shù)值差分計(jì)算方法編制程序計(jì)算不同工況下船體結(jié)構(gòu)變形數(shù)據(jù)，具體數(shù)據(jù)如表3所示。

本文的仿真計(jì)算結(jié)果與經(jīng)驗(yàn)公式計(jì)算值相比數(shù)值較小，這是由于現(xiàn)代造船技術(shù)和金屬材料生產(chǎn)技術(shù)比二戰(zhàn)時(shí)期有所提升，使得艦船結(jié)構(gòu)的總體強(qiáng)度和局部強(qiáng)度都有提高，從理論角度看方法是可行的。

反艦導(dǎo)彈戰(zhàn)斗部在艦船艙室內(nèi)部爆炸會(huì)給艦船局部結(jié)構(gòu)造成嚴(yán)重破壞，2發(fā)以上的反艦導(dǎo)彈武器直接命中船體的不同部位，足以使一艘2 000t左右的艦船喪失或部分喪失戰(zhàn)斗力，如果綜合考慮火災(zāi)、戰(zhàn)斗部破片及艦船技術(shù)裝備破壞等因素，艦船的損毀情況會(huì)更嚴(yán)重，對(duì)艦船的戰(zhàn)斗力破壞會(huì)更大。

表3 艦船板架結(jié)構(gòu)等效板厚變形計(jì)算結(jié)果

從近場(chǎng)接觸爆炸板架變形實(shí)驗(yàn)來(lái)看［9，10］，艦船板架結(jié)構(gòu)的變形呈花瓣?duì)铋_(kāi)裂，花瓣?duì)铋_(kāi)裂部分會(huì)發(fā)生卷邊，本文仿真計(jì)算并未考慮這些因素，其計(jì)算結(jié)構(gòu)破口的撓度值稍偏大。

5 結(jié)論

本文采用仿真計(jì)算方法對(duì)某船體板架結(jié)構(gòu)在無(wú)限介質(zhì)中近場(chǎng)爆炸載荷作用下其動(dòng)態(tài)響應(yīng)進(jìn)行了計(jì)算，主要結(jié)論如下：①通過(guò)理論指導(dǎo)在非線(xiàn)性有限元軟件中建立了無(wú)限介質(zhì)中近場(chǎng)爆炸載荷計(jì)算的有限元模型，從仿真計(jì)算結(jié)果分析來(lái)看，此方法的計(jì)算結(jié)果對(duì)開(kāi)展艦船結(jié)構(gòu)破壞研究工程實(shí)踐具有指導(dǎo)意義；②艦船結(jié)構(gòu)破壞動(dòng)態(tài)響應(yīng)計(jì)算的初始條件可采用應(yīng)力波理論計(jì)算方法獲得，有利于提高計(jì)算精度；③通過(guò)仿真計(jì)算與經(jīng)驗(yàn)公式計(jì)算的對(duì)比分析，說(shuō)明本文所建立的無(wú)限介質(zhì)中近場(chǎng)爆炸載荷作用下艦船結(jié)構(gòu)動(dòng)態(tài)響應(yīng)的計(jì)算模型，可作為開(kāi)展艦船結(jié)構(gòu)破壞研究工程實(shí)踐中結(jié)構(gòu)變形的估算方法。

［1］姚熊亮.艦船結(jié)構(gòu)振動(dòng)沖擊與噪聲［M］.北京：國(guó)防工業(yè)出版社，2007.

［2］Calder C A，Kelly J M，Goldsmith W.Projective impact on an infinite viscoplastic plate［J］.Int J Solid Structures，1971（7）：1143-1152.

［3］Henrych J.The dynamics of explosive and its use［M］.Oxford：Elsevier Scientific Publishing Company，1979.

［4］寧建國(guó)，王成，馬天寶.爆炸與沖擊動(dòng)力學(xué)［M］.北京：國(guó)防工業(yè)出版社，2010.

［5］王禮立，應(yīng)力波基礎(chǔ)［M］.北京：國(guó)防工業(yè)出版社，2010.

［6］楊桂通，塑性動(dòng)力學(xué)［M］.北京：高等教育出版社，2000.

［7］朱錫，白雪飛，張振華.空中接觸爆炸作用下船體板架塑性動(dòng)力響應(yīng)及破口研究［J］.中國(guó)造船，2004，45（2）：43-50.

［8］吉田隆.二次世界大戰(zhàn)初期日本海軍艦船在炸彈攻擊下的損傷實(shí)例分析［J］.船の科學(xué)，1990（5）：70-81.

［9］吳有生，彭興寧，趙本立.爆炸載荷作用下艦船板架結(jié)構(gòu)的變形與破損［J］.中國(guó)造船，1995（4）：55-61.

［10］劉潤(rùn)泉，白雪飛，朱錫.艦船單元結(jié)構(gòu)模型水下接觸爆炸破口實(shí)驗(yàn)研究［J］.海軍工程大學(xué)學(xué)報(bào)，2001，13（5）：41-46.