艦船目標(biāo)在水下爆炸作用下毀傷效應(yīng)研究?

趙紅光

（91439部隊 大連 116041）

艦船目標(biāo)在水下爆炸作用下毀傷效應(yīng)研究?

趙紅光

（91439部隊 大連 116041）

采用數(shù)值計算以及理論分析方法對水下爆炸作用下目標(biāo)的毀傷效果進行計算分析。首先結(jié)合試驗數(shù)據(jù)驗證了數(shù)值計算方法的正確性，然后選擇計算目標(biāo)對其在不同工況下的毀傷情況進行計算，同時采用理論分析方法對數(shù)值計算結(jié)果與實船計算結(jié)果分析比對，結(jié)果表明在近場爆炸作用下，只考慮毀傷的情況下可用艙段代替目標(biāo)艦船。

爆炸；數(shù)值；計算毀傷；目標(biāo)

1 引言

由于實船水下爆炸試驗耗費大，實施困難，對環(huán)境也會造成不利影響，無法通過大規(guī)模的試驗獲得數(shù)據(jù)。而隨著計算機硬件和計算方法的不斷發(fā)展進步，水下爆炸數(shù)值研究得到了各國學(xué)術(shù)界的高度重視，已成為研究水下爆炸現(xiàn)象及機理的一項重要手段［1～7］。

而集群系統(tǒng)作為一種高性能計算平臺，其最高可對1000萬網(wǎng)格的水下爆炸仿真模型進行仿真計算，其采用的并行計算技術(shù)和搭載的128核計算軟件能實現(xiàn)對水下爆炸作用機理的更精確模擬和快速計算。

2 仿真計算方法

ABAQUS采用聲固耦合法進行水下爆炸仿真計算。將流體定義為聲學(xué)介質(zhì)，沖擊波在聲學(xué)介質(zhì)中傳播。對于同時涉及固體和聲學(xué)單元的問題，存在兩種聲學(xué)壓力場的計算公式，即總波公式和散射波公式。當(dāng)入射波載荷施加時?？偛ü较啾壬⑸洳ü礁咏诮Y(jié)構(gòu)載荷。當(dāng)聲學(xué)介質(zhì)由于流體的空化是非線性或考慮空化效應(yīng)的影響時，必須使用總波公式［8］。

假定流體為可壓縮、絕熱且運動為小位移［9］，有速度相關(guān)動量損失的流體動量方程可表示為

其中，p為流體動壓力（超過靜壓力的部分）；x為流體質(zhì)點的空間位置坐標(biāo)；v?f和 v?f為流體質(zhì)點速度和加速度；Df為流體密度；α為單位體積單位速度的力；Ki為相關(guān)的場變量（如溫度、濕度或鹽度等）。

再假定流體為無粘性、線性可壓縮，則流體本構(gòu)方程可表示為

其中，vf為流體質(zhì)點位移，Rf為流體體積模量。

為了得到直接積分瞬態(tài)分析的偏微分方程，將式（1）兩邊除以Df并對x求導(dǎo)。假定忽略α Df的微分項，將結(jié)果與式（2）的時間導(dǎo)數(shù)合并得到用壓力表示的流體微分方程：

引入任意的變化域δp，在整個流體域中對式（3）積分得到運動方程的等效弱形式：

通過ABAQUS的聲固耦合分析［10］得到流體域平衡方程：

其中，δvm為位移變分，τ為結(jié)構(gòu)阻力，σ為節(jié)點應(yīng)力，p為作用在結(jié)構(gòu)濕表面的壓力，n為指向流體的結(jié)構(gòu)表面法向量，D為結(jié)構(gòu)密度，αc為質(zhì)量比例阻尼系數(shù)，vm、v?m和 v?m分別為結(jié)構(gòu)單點位移、速度、加速度，δε為虛應(yīng)變對應(yīng)的虛位移。

這樣就得到了流體域中的結(jié)構(gòu)平衡方程。接著利用有限元法將結(jié)構(gòu)和聲學(xué)介質(zhì)進行離散化，定義壓力加載面，將Geers and Hunters模型［11］計算出的水下爆炸壓力載荷加載到離散化的表面上。最后利用顯式時間積分法通過求解離散化的式（5）和式（6）得到結(jié)構(gòu)響應(yīng)和流體域中的壓力傳播。

3 仿真方法保真度驗證

結(jié)合某實船爆炸數(shù)據(jù)，對仿真方法進行驗證。計算工況為1000kg TNT距左舷70m，水深50m沉底爆炸，龍骨沖擊因子為0.337。根據(jù)當(dāng)時爆炸時測點的實際情況選取進行比對的測點位置，分別選擇了3個應(yīng)變測點和3個加速度測點。測點具體位置分別在：01甲板55#肋位迎爆面中縱桁縱向應(yīng)變測點、主甲板25#肋位迎爆面縱桁縱向應(yīng)變測點、內(nèi)底45#肋位迎爆面第一扶強材邊中間根部板格應(yīng)變測點、01甲板27#肋位室垂向加速度測點、主甲板76#肋位垂向加速度測點、內(nèi)底47#肋位垂向加速度測點。對比結(jié)果見表1～2。

表1 應(yīng)變峰值對比

表2 加速度峰值對比

通過對峰值的相對誤差計算，得出應(yīng)變峰值平均精度77.48%，加速度峰值平均精度77.52%，平均相對誤差均在30%以內(nèi)，計算精度符合工程要求。

4 計算目標(biāo)

4.1 計算對象

計算對象以某艦為原型，截取艦中兩個艙段，仿真試驗?zāi)Ｐ腿鐖D1～4所示，圖1為艙段內(nèi)部結(jié)構(gòu)有限元模型圖，圖2為艙段加強結(jié)構(gòu)與內(nèi)部結(jié)構(gòu)圖，圖3為水域有限元模型圖，圖4為艙段仿真計算整體有限元模型圖。

4.2 仿真工況設(shè)計及測點選取

仿真試驗共設(shè)置了5個工況，工況設(shè)計主要依據(jù)水中兵器的裝藥量及使用情況，其中第一個工況采用驗證算法工況相同的沖擊因子。

模型綜合考慮迎背爆面、船體結(jié)構(gòu)、設(shè)備位置等因素，共設(shè)置了10個測點，測點位置均位于模型中橫剖面位置。

5 計算結(jié)果及分析

對數(shù)值計算結(jié)果進行分析，并與實船仿真結(jié)果進行比對。

5.1 仿真計算結(jié)果

仿真計算結(jié)果如圖5～圖7所示，圖5仿真計算水域壓力云圖，圖6仿真計算結(jié)構(gòu)應(yīng)力云圖。

選取底部中心測點處的仿真計算結(jié)果表明，雖然在結(jié)構(gòu)形式上存在一定的差異，但艙段模型和整船模型的水平方向與垂直方向的加速度、速度曲線值在同等數(shù)量級上且具有可比性，符合設(shè)計要求，可用于模擬實際戰(zhàn)場環(huán)境。

5.2 與實船仿真結(jié)果比對

為了對比分析實船與艙段之間的毀傷參數(shù)的異同，將艙段與實船方正計算結(jié)果進行比對，結(jié)果如圖7～圖10所示。

圖7～圖8為測點垂向加速度比對結(jié)果，由圖可知兩者加速度的變化趨勢基本一致，但是艙段加速度幅值大于實船。

圖9～圖10為等效靶及實船的應(yīng)力云圖，由圖可知等效靶與實船毀傷情況基本一致。由毀傷理論分析可知，在近場毀傷作用下，艙段和實船滿足材料、相當(dāng)板厚、外部形狀、內(nèi)部結(jié)構(gòu)基本與實船相同的條件下，就可保證兩者抗毀傷效果相同，比對結(jié)果驗證了該理論的正確性，為毀傷等效靶設(shè)計提供一定參考。

綜上所述，艙段模型均是采用與原始艦艇等同的艙段，沒有進行幾何的縮比改動，在沖擊加速度趨勢上基本相同，量值上不具備可比性，可作為考核設(shè)備沖擊振動的平臺使用；從兵器毀傷角度考慮，可用艙段替代實際艦船。

6 結(jié)語

通過本文的研究，得到一些有意義的結(jié)論

1）通過與試驗數(shù)據(jù)的比對，證明了聲固耦合法計算爆炸毀傷數(shù)值計算方法的正確性，為后續(xù)的仿真計算提供方法和路徑。

2）采用數(shù)值計算方法對目標(biāo)艦船的艙段及目標(biāo)的毀傷情況進行仿真計算，比對兩者的毀傷情況，發(fā)現(xiàn)兩者毀傷情況相類似，可用艙段替代實際目標(biāo)。

3）從艙段和目標(biāo)的沖擊環(huán)境來看，兩者的變化趨勢比較一致，但艙段的加速度量值高于實船，用艙段替代實船時需對沖擊環(huán)境進行修正。

4）采用數(shù)值仿真方法與同等級別整船進行了對比驗證，只考慮毀傷的情況下可用艙段代替目標(biāo)艦船。

［1］張振華，汪玉，張立軍等.船體梁在水下近距爆炸作用下反直觀動力行為的相似分析［J］.應(yīng)用數(shù)學(xué)和力學(xué)，2011（12）：78-63.

［2］ JaremaM.Didoszak.Parametric studies of DDG-81 ship shock trial simulations［D］.Naval Postgraduate School，2004.

［3］姚熊亮，楊文山，初文華等.艦船水下接觸爆炸數(shù)值研究［J］.計算物理，2010（11）：105-110.

［4］牟金磊，朱錫，黃曉明等.水下近場非接觸爆炸作用下固支方板破口計算［J］.振動與沖擊，2011：64-70.

［5］Young S.Shin.Ship shock modeling and simulation for far-field underwater explosion［J］.Computers and Struc?tures，2004，82：2211-2219.

［6］姚熊亮，李克杰，張阿漫等.艦船水下爆炸新的損傷評估體系探索［J］.中國造船，2008（09）：81-87.

［7］楊樹濤，朱永凱，焦磊.艦船舷側(cè)防護結(jié)構(gòu)水下接觸爆炸動響應(yīng)分析研究［J］.中國艦船研究，2010（08）：116-121.

［8］DassaultSystèmesSimulia Corp.ABAQUS Analysis User’s Manuals，version 6.9［Z］.RI，USA：DassaultSystèmes，2009.

［9］ DassaultSystèmesSimulia Corp.ABAQUS Theory Manu?als，version 6.9［Z］.RI，USA：DassaultSystèmes，2009.

［10］李國華、李玉節(jié)、張效慈等.浮動沖擊平臺水下爆炸沖擊譜測量與分析［J］.船舶力學(xué)，2000（4）：54-59.

［11］GEERS T L，HUNTER K S.An integrated wave-effect s model for an underwater explosion bubble［J］.Journal of the Acoustical Society of America，2002，111 （4）：1584-1601.

［12］張瑋，崔立.艦船抗沖擊試驗仿真技術(shù)［J］.科技導(dǎo)報，2009，27（14）：38-41.

Research on Warship Target Damage Effect Subject to Under Water Explosion

ZHAO Hongguang
（No.91439 Troops of PLA，Dalian 116041）

The under-water explosion damage status of warship has been calculated using numerical simulation and theory analysis method.First，validity of simulation method has been testified by comparing experiment data.Second，proper target and in?stances have been selected，and then target damage simulation has been operated.The target and warship simulation results have been analyzed by using theory method.It shows that the target can substitute warship in the damage status.

explosion，numerical compute，damage，target

U661.2

10.3969/j.issn.1672-9730.2017.11.026

Class Number U661.2

2017年5月4日，

2017年6月11日

趙紅光，男，助理工程師，研究方向：水中兵器試驗研究。