遭受水下爆炸的艦船縮比模型毀傷評估

程素秋，陳高杰，王樹樂

遭受水下爆炸的艦船縮比模型毀傷評估

程素秋，陳高杰，王樹樂

（中國人民解放軍 海軍91439部隊，遼寧 大連 116041）

對某型艦船縮比模型遭受水下爆炸的彈塑性動力響應進行數值仿真與實爆試驗研究，基于總縱強度毀傷和局部塑性變形的表征參數進行毀傷等級評估。結果表明，加速度響應和毀傷等級隨沖擊因子增大而增大，局部響應與船體振型相關，局部塑性變形毀傷等級最高時總縱強度毀傷等級并不是最高。研究成果可為實船爆炸試驗評估提供參考。

振動與波；艦船縮比模型；水下爆炸；動力響應；數值模擬

遭受水下爆炸的艦船結構彈塑性動力響應研究對于提高現役艦船的生命力和戰斗力，具有重要意義。因實船爆炸試驗的巨額經費、超長周期和不可重復等特點，在實際研究中常用縮比模型爆炸試驗來做先期研究。1984年，Rentz用試驗手段研究加筋平板在水下爆炸作用下的動態響應，并進行了數值計算［1］。1988年，Gifford等對具有初始裂紋焊接厚板的動態響應作了一系列試驗研究［2］。胡俊波等提出了改進型沖擊因子的計算方法，并經某目標抗水下爆炸試驗驗證［3］。陳輝等利用整船縮比模型試驗，得到沖擊環境參數與水下爆炸載荷間的關系［4］。朱錫等對水面艦艇防雷艙結構模型進行了一系列水下抗爆能力試驗［5］。張海鵬等研究對大縮尺比船模，使用畸變模型實現對原型總縱強度準確預測的方法［6］。秦健等基于π定律提出了在水下爆炸作用下加筋板模型動態響應的相似預報方法［7］。程素秋等對某典型艙段模型進行水下爆炸試驗，研究了模型結構所遭受的載荷及其動態響應［8］。高浩鵬等設計了多層隔振系統用于水下爆炸試驗時測量設備抗沖擊防護［9］。張曉陽等采用頻域法和時域法對液壓螺栓有接觸應力時抗沖擊性能進行仿真［10］。

1　艦船縮比模型設計

以某型艦為母船，按總縱強度相似原則，對母船作簡化，并忽略母船艦載設備對其質量分布的影響。縮比模型（以下簡稱船模）共分8個艙，縮比因子1：20，不另行配重。其中艏艙較長（是其它艙室的兩倍長），其余艙室長度一樣。船模共設7道橫艙壁，2層平臺和中縱艙壁，平臺2與船底外底板組成雙層底，雙層底內還有兩道縱旁艙壁加強，甲板、平臺1、平臺2的板厚4 mm，其余結構的板厚3 mm，模型結構總質量約m=429.4 kg，材料選用普通船用鋼Q235；內部結構示意圖見圖1。兩端分別設有水密間，可以保證模型在殼體破裂后仍能浮在水面。

圖1　某型艦船縮比模型結構剖面圖

2　數值仿真研究

2.1有限元模型建立

計算時考慮重力因素和殼板的應變率強化效應。有限元模型中船體網格由59 034個殼單元（S4R）組成，水域網格由3 615 054個聲學單元（AC3D8R）組成。結構與水域間用TIE約束法模擬流固耦合，采用ABAQUS/EXPLICIT求解器（聲固耦合法）計算。船模和外部水域的有限元模型分別如圖2、圖3所示。

圖2　結構有限元模型

圖3　流體模型圖

圖4　Q235鋼的動態力學性能參數

通過擬合得到：D=271.6，取ρ=7 800 kg/m3。E=2.1×1011Pa，材料極限應變取εf=0.3。

炸藥采用 Jones-Wilkins-Lee（JWL）狀態方程［12］。JWL描述爆壓P和單位體積內能e及爆炸產物體積V之間的關系，即

式中 A=373.77 GPa；B=3.747 1 GPa；R1=4.15；R2=0.9；ω=0.35；η=ρ/ρ0、ρ0=1.63 g/cm3、ρ為爆炸產物密度；TNT的e=6.0 GJ/m3。爆速為6 930 m/ s，爆壓為21 GPa。

2.2計算工況

取計算工況和實爆試驗工況完全相同。船模爆炸試驗是在某海域的近海岸完成，試驗場景如圖7所示。炸藥為1 kg RDX球形藥包和6 kg TNT球形藥包。對船模進行了七次水下爆炸，每次試驗時選用的藥包、布放距離、龍骨沖擊因子詳見表1，最后一次為毀傷試驗，沒有進行測量。

表1　船模計算工況表

2.3船模測點分布

船模上布設包括自由場壓力、加速度和應變等測點。試驗前在船模上要按照預定位置裝好傳感器。加速度、應變測點位置及標號如圖5、圖6所示，應變片貼在平臺1與平臺2間的內壁上。

圖5　平臺2加速度測點位置

圖6　底部加速度測點位置

圖7　船模爆炸試驗場景

3　船模結構彈塑性動力響應分析

3.1水下爆炸沖擊波載荷分析

Cole和Zamyshlyayev相繼提出了計算水下沖擊波載荷的經驗公式（2）與（3），來計算沖擊波壓力P（t）、峰值壓力Pm等。其中：W為藥量/kg，R為爆心距觀測點的距離/m，R0為藥包的初始半徑/m，早期的沖擊波以指數形式衰減，當t＞θ以后，沖擊波的衰減變慢，以近似于時間倒數的關系衰減。

圖8和圖9表示自由場壓力曲線，圖中是實測曲線與計算曲線的對比，二者趨勢基本相符，但實測曲線較計算曲線更復雜，所包含的載荷更豐富，反映了自由面截斷和海底反射效應（注：工況3與工況4選擇工況3；工況5與工況6選擇工況5）。自由場壓力峰值對比見表2。

圖8　仿真計算與實爆試驗的自由場壓力對比

3.2船模在水下爆炸作用下的動響應規律

船模爆炸試驗的沖擊因子是按照炸藥當量由小到大，爆距從遠及近的原則而設計的，為便于對比，前三次試驗使用相同的藥包，后三次試驗藥包一樣。另外，工況1和工況2、工況3和工況4的沖擊因子基本相同，工況5和工況6只是藥包距船模水平距離不同；工況3和工況5只是藥包TNT當量不同，這是為考察單一因素作用下的船模動響應規律而設置的。

從表3中數據可知，同一測點在相同沖擊因子作用下的響應也是有差異的，工況2、4中各測點響應值明顯高于工況1、3（個別值偏離大的是因其零漂現象引起的）；同一工況中各測點的動響應也是不同的，離爆源越近，響應值越大，由于藥包布放時偏向船模首部，因此首部測點的響應值要大于尾部測點；距離爆源位置差不多時，船模底部中間處的響應值要大于首尾兩端及1/4處；平臺2上的響應要小于底部。這是由于艦船是具有結構振型的，初始傳遞到船底的沖擊波作為船底動能的一部分，以剛體運動的形式出現，其余爆炸能量則引起艦船結構以不同結構振型作振蕩變形。當沖擊能量在船體結構中向上傳播時，高頻分量會因結構振型的作用而衰減，而低頻振型的響應就變得更加突出。

圖9　測點BOTTOM-2-L在各次工況中的加速度時間歷程比對

4　船模結構在水下爆炸作用下的毀傷評估

水面艦船在水下爆炸作用下的毀傷，可分為總體毀傷和局部毀傷?？傮w毀傷又分為：總縱強度毀傷和剩余總縱強度毀傷。前者是由中近場大當量水

下沖擊波和氣泡脈動壓力作用下艦艇結構產生鞭狀運動、及近距爆炸產生的沖擊振動與底部負壓區使艦艇產生中垂狀態的迭加作用等引起。后者是在接觸爆炸作用下造成艦艇的局部破口，引起艦艇結構間接的總體毀傷，屬于艦艇結構破損剩余強度。水面艦船局部毀傷又可分為：局部破損毀傷和局部塑性變形毀傷。前者是由近距或接觸爆炸作用下引起

的；后者是由非接觸爆炸作用下引起的。針對表1所列的工況，船模的毀傷形式有兩種：總縱強度毀傷和局部塑性變形毀傷。

表2　自由場壓力峰值對比

表3　船模爆炸試驗的加速度峰值統計

4.1總縱強度毀傷計算

式中Mu1為完好船體梁計算剖面彈性極限彎矩值；Mq1為計算狀態對應的船體梁剖面所受靜、動外力合成作用彎矩的最大設計值。即

其中Ms1為完好艦艇靜水彎矩值，Me1為氣泡載荷附加彎矩值，Md1為最大設計沖擊振動彎矩值。

4.2局部塑性變形毀傷計算

假設加筋板結構為四周固支，在爆炸載荷作用下加筋板結構發生大撓度變形，考慮到結構邊界產生塑性鉸線，最終變形的撓曲面可近似為正弦曲面。其撓曲面函數可取為

（1）板架邊界塑性鉸彎曲變形能U1

（2）板架彎曲變形能U2

（3）板架伸長變形能U3

（4）結構初始動能Ek

相對于結構響應時間，沖擊波持續時間非常短，可近似認為炸藥的沖擊波能全部被結構吸收。在非接觸水下爆炸條件下，可假設沖擊波為平面波，由沖擊波的能流密度表達式得到結構吸收的沖擊波能，即沖擊波作用后結構的初始動能

（6）能量守恒原理

根據能量守恒原理有U1+U2+U3=Ek。將式（6）—式（10）和參數代入，進行求解可以得到中心撓度w0，進而可以得到整個加筋板結構的變形場。即

4.3局部塑性變形算例分析

實爆試驗結果表明，舷側板架（上）的變形較嚴重，因此取該板架作為局部撓曲變形的分析對象。

對比表4和表5，發現同一工況針對不同參數的評估結果是不同的。評估時針對不同參數可分別作結論，然后綜合全部評估結果衡量后得出在某個沖擊因子下會達到的效果及量級。圖10是工況7的數值模擬與實爆試驗的塑性變形對比。對比圖中由變形情況可知，數值模擬與實爆試驗結果基本相符。

表4　船模總縱強度抗毀傷等級評估結果

表5　舷側板架局部撓曲變形毀傷等級評估結果

圖10　船模在工況7后的塑性變形對比

5　結語

對某型艦船縮比模型遭受水下爆炸時的彈塑性動力響應進行數值仿真與一系列實爆試驗研究。基于總縱強度毀傷和局部塑性變形對該模型進行了毀傷評估。結果表明，加速度響應和毀傷等級隨沖擊因子增大而增大，局部響應與船體振型相關，局部塑性變形毀傷等級最高時總縱強度毀傷等級并不是最高。研究成果可為實船爆炸試驗評估提供參考。

［1］RENTZ T R.An experimental investigation into the dynamicresponseofastiffenedflatplateloaded impulsively by an underwater shockwave，A123151［R］. 1984.

［2］GIFFORD L N，CARLBERG J R，WIGGS A J，et al. Explosive testing of full thickness precracked weldments ［J］.ASTM Special Technical Publication，1990（1074）：157-177.

［3］胡俊波，張志華，李慶民，等.水下爆炸對水下目標的毀傷評估研究［J］.振動與沖擊，2010，29（10）：206-210.

［4］陳輝，李玉節，潘建強，等.水下爆炸條件下不同裝藥對水面艦船沖擊環境的影響試驗研究［J］.振動與沖擊，2011，30（7）：16-20.

［5］朱錫，張振華，劉潤泉，等.水面艦艇舷側防雷艙結構模型抗爆試驗研究［J］.爆炸與沖擊，2004，24（2）：133-139.

［6］張海鵬，岳永威，蘇羅青，等.水下爆炸作用下艦船總縱強度模型實驗方案研究［J］.振動與沖擊，2012，31（6）：175-180.

［7］秦健，張振華.原型和模型不同材料時加筋板沖擊動態響應的相似預報方法［J］.爆炸與沖擊，2010，30（5）：511-516.

［8］程素秋，樊寶順，薛飛，等.水下非接觸爆炸作用下艙段模型的動態響應［J］.爆炸與沖擊，2008，28（4）：360-365.

［9］高浩鵬，黃映云，金輝，等.水下爆炸作用下多層隔振系統動態特性分析［J］.噪聲與振動控制，2015，35（3）：186-188.

［10］張曉陽，陳汝剛，解忠良，等.艦船軸系用液壓螺栓抗沖擊性能的數值仿真［J］.噪聲與振動控制，2015，35（4）：149-152.

Damage Evaluation of a Ship Scaled Model Subjected to Underwater Explosion

CHENG Su-qiu，CHEN Gao-jie，WANG Shu-le
（91439 Unit of PLA，Dalian 116041，Liaoning China）

The dynamic response of a ship scaled model subjected to underwater explosion is studied using numerical simulation and realistic explosion tests.Based on the characteristic parameters of longitudinal strength and local plastic distortion，the grade evaluation of anti-deflection capability of the model is conducted.The results show that the acceleration response and damage grade increase with the shock factors，and the local response is related to the vibration mode.When the damage grade based on the local plastic distortion reaches the highest value，the damage grade based on longitudinal strength is not necessarily the highest.The result is of significance for evaluation of ship shock tests.

vibration and wave；ship scaled model；underwater explosion；dynamic response；numerical simulation

U663.85

ADOI編碼：10.3969/j.issn.1006-1335.2016.03.015

1006-1355（2016）03-0070-06

2015-06-10

程素秋（1968-），遼寧省丹東市，女，碩士，高級工程師，從事水下爆炸試驗仿真與評估工作。E-mail：csq4028@sina.com