任 鵬， 張 偉， 葉 楠， 黃 威

（1哈爾濱工業(yè)大學 高速撞擊研究中心，哈爾濱150080；2江蘇科技大學 船海學院，江蘇 鎮(zhèn)江212003）

1 引 言

大型艦艇作為海軍的重要裝備在維護國家海權(quán)，保護海上活動安全等方面起著舉足輕重的作用。然而隨著現(xiàn)代精確制導武器和高速水中兵器的快速發(fā)展，艦艇的安全受到了日趨嚴重的威脅，如何提高艦艇結(jié)構(gòu)的抗爆抗沖擊能力成為了衡量艦艇生命力最重要的技術(shù)指標之一。但由于炸藥爆炸實驗存在危險程度高，花費大，數(shù)據(jù)采集精度和廣度無法全部保證等缺點，導致沖擊波載荷作用下艦艇結(jié)構(gòu)的抗爆抗沖擊性能實驗無法大規(guī)模展開，嚴重影響了艦艇防護結(jié)構(gòu)的發(fā)展?；诖?，劍橋大學的Deshpande和Fleck等[1]提出了一種柱形非藥式水下爆炸裝置，并利用該裝置對具有吸能核心材料的夾層板結(jié)構(gòu)進行了沖擊波加載試驗，基于Taylor在1941年提出的水下一維沖擊波動量傳遞理論，比較了耦合算法和解耦算法對于夾層吸能材料在沖擊波載荷作用下變形問題的適用性；美國西北大學的Espinosa、Mori和Latourte等[2-4]利用類錐型非藥式水下爆炸模擬裝置結(jié)合高速攝影技術(shù)對多種結(jié)構(gòu)板進行了水下大當量爆炸沖擊波加載實驗，得到了爆炸沖擊載荷作用下合金板、復(fù)合材料板和夾芯板的毀傷變形情況；Huson等[5]基于反力墻技術(shù)自制了非藥式水下爆炸沖擊波模擬裝置，該裝置能夠?qū)^大的結(jié)構(gòu)進行水下沖擊波加載實驗；國內(nèi)的李伯松等[6]也對非藥式爆炸模擬技術(shù)進行了探索，成功設(shè)計出了爆炸模擬裝置，但該裝置無法對水下爆炸問題進行模擬；哈爾濱工業(yè)大學的任鵬和張偉等[7-9]對非藥式水下爆炸技術(shù)進行了深入研究，利用該技術(shù)對層合結(jié)構(gòu)板進行了水下爆炸沖擊波加載試驗，得到了層合板的動態(tài)響應(yīng)規(guī)律。這些研究在一定程度上緩解了實驗室進行大當量水下爆炸試驗所面臨的問題。本研究利用實驗與仿真相結(jié)合的方法，對飛片及活塞對水下爆炸沖擊波初始峰值及其衰減規(guī)律的影響進行了分析，建立了非藥式水下爆炸沖擊波特性與飛片及活塞相關(guān)參數(shù)的表達式，為后續(xù)艦艇結(jié)構(gòu)的水下抗沖擊實驗研究奠定了基礎(chǔ)。

2 實驗研究

2.1 實驗裝置

非藥式水下爆炸沖擊波模擬實驗裝置如圖1（a）所示，該裝置包括一級輕氣炮、激光測速系統(tǒng)、水容器及沖擊波測試系統(tǒng)，利用輕氣炮發(fā)射飛片正撞擊活塞可以在水容器中產(chǎn)生水下爆炸沖擊波[1-2]。

圖1 非藥式水下爆炸沖擊波模擬實驗裝置示意圖Fig.1 Experimental set-up sketches of non-explosive underwater shock simulation

該裝置中水容器長L=406 mm，其中水容器內(nèi)部柱形部分長L1=76 mm，壓力傳感器位置L2=89 mm，L3=140 mm，活塞直徑D1=66 mm，靶板受沖擊部分直徑D2=152 mm，靶板整體直徑D3=256 mm，水艙內(nèi)部散射角 γ=7°，如圖1（b）所示。

2.2 實驗結(jié)果

圖2為實驗測得的典型壓力歷程曲線。其中圖2（a）和圖2（b）分別為質(zhì)量是0.282 kg、速度是60.27 m/s和質(zhì)量為0.572 kg、速度為70.55 m/s的飛片正撞擊活塞，由A傳感器測得的水中壓力歷程曲線。從這兩個實驗中可以發(fā)現(xiàn)，當沖擊波傳播到A傳感器位置時，該處的壓力急劇增加到峰值載荷，然后以指數(shù)形式衰減。且由圖3可見，在A處測得的質(zhì)量為0.572 kg的飛片正撞擊活塞對應(yīng)的壓力峰值98 MPa，壓力從峰值載荷衰減到零持續(xù)時間約0.13 ms，均高于質(zhì)量為0.282 kg的飛片正撞擊活塞產(chǎn)生的壓力峰值69.2 MPa和壓力峰值衰減時間0.1 ms。

圖2 數(shù)值模擬結(jié)果和實驗結(jié)果對比Fig.2 Comparison of the results from calculations and test

3 數(shù)值模擬方法

為了更加深入地探究非藥式水下爆炸沖擊波峰值特性與飛片及活塞間的關(guān)系，利用AUTODYN-2D程序，對非藥式水下爆炸沖擊波加載裝置進行了數(shù)值模擬研究。

3.1 數(shù)學模型

非藥式水下爆炸問題涉及流固耦合，單一的Lagrange或Euler算法均無法準確地描述整個物理過程。本研究采用Euler-Lagrange耦合算法，對流體介質(zhì)的波動與水艙的相互作用進行了耦合計算。其中，水介質(zhì)及附近區(qū)域定義為Euler網(wǎng)格，艙壁結(jié)構(gòu)定義為Lagrange網(wǎng)格，水介質(zhì)采用均布網(wǎng)格劃分。并在非藥式水下爆炸沖擊波模擬裝置的水平中軸線上定義高斯觀測點，用以分析水中爆炸沖擊波的傳播規(guī)律。

3.2 材料模型

水的Shock狀態(tài)方程能夠較好地描述水中爆炸沖擊波的傳播和衰減[10]，具體參數(shù)見表1。

表1 水的材料性能參數(shù)Tab.1 Material parameters of water

非藥式水下爆炸實驗裝置中的飛片和活塞材料均為S-7模具鋼，艙體材料和靶板材料為4340鋼，具體參數(shù)分別見表2和表3。

表2 4340鋼材料性能參數(shù)Tab.2 Material parameters of 4340 steel

表3 S-7鋼材料性能參數(shù)Tab.3 Material parameters of S-7 steel

表2，3中：ρ0為水的初始密度，c0為水中聲速，s1為狀態(tài)方程常數(shù)，Γ0為格林愛森常數(shù)，E為楊氏模量，v為泊松比，σ0為屈服強度，K為硬化常數(shù)，n為硬化指數(shù)，ε0為應(yīng)變率常數(shù)。

3.3 數(shù)值模擬結(jié)果有效性驗證

為了驗證有限元方法對該類問題研究的有效性，首先對上述實驗進行了數(shù)值模擬，模擬中所用靶板厚度為25 mm。數(shù)值模擬結(jié)果與實驗結(jié)果的對比如圖2所示。從圖2可見，數(shù)值模擬得到的壓力歷程曲線和實驗測量得到的曲線基本符合良好。重要特征如：第一個壓力峰值及其指數(shù)衰減過程均與實驗結(jié)果保持了良好的一致性，從而驗證了計算結(jié)果的可靠性，由此表明，本文所建立的數(shù)值模型能夠較好地模擬非藥式水下爆炸沖擊問題。圖3給出了飛片質(zhì)量為0.282 kg、速度為60.27 m/s時，撞擊活塞得到的壓力波波陣面仿真進程，其中時間零點為飛片撞擊活塞時，由圖3可見，沖擊波波陣面在傳播過程中保持了良好的平面性。

圖3 沖擊波波陣面歷程圖Fig.3 Time histories of the shock wave front

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

本研究對飛片質(zhì)量為0.265 kg、0.795 kg、1.325 kg、1.855 kg 和 2.385 kg，飛片速度分別為 50 m/s、100 m/s、150 m/s、200 m/s和250 m/s的工況進行數(shù)值計算，為了便于比較分析，活塞厚度分別取12 mm、20 mm和30 mm。

圖4 活塞厚度為12 mm時水下爆炸沖擊波初始峰值與飛片速度的關(guān)系Fig.4 Relationship between shock wave strength and flyer velocity when the piston length is 12 mm

4.1 飛片對水下爆炸沖擊波強度的影響

圖4為五種質(zhì)量飛片在不同速度條件下正撞擊活塞得到的非藥式水下爆炸沖擊波初始峰值P變化曲線。從圖4可以看出，當活塞厚度為12 mm時，不同質(zhì)量的飛片在同一撞擊速度條件下，產(chǎn)生的水下爆炸初始沖擊波峰值基本相同，且沖擊波初始峰值隨飛片撞擊速度的增加而增加，成線性趨勢。這說明，當活塞一定時，非藥式水下爆炸沖擊波初始峰值的大小僅與飛片的初始撞擊速度有關(guān)，飛片的質(zhì)量對其并無顯著影響。

4.2 活塞對水下爆炸沖擊波強度的影響

圖5分別為質(zhì)量為0.265 kg和1.325 kg的飛片在不同速度條件下分別撞擊12 mm、20 mm和30 mm活塞產(chǎn)生的水下沖擊波初始峰值。從圖5（a）中可見，當飛片質(zhì)量為0.265 kg時，隨著活塞厚度的增加，飛片速度為50 m/s，100 m/s和200 m/s時，20 mm厚活塞和30 mm厚活塞所對應(yīng)的沖擊波初始峰值與相同工況下12 mm厚活塞對應(yīng)的初始峰值一致性良好，但飛片速度為150 m/s時，30 mm厚活塞對應(yīng)的沖擊波初始峰值要明顯小于12 mm和20 mm厚活塞所對應(yīng)的初始峰值，但該現(xiàn)象在飛片速度為200 m/s時明顯減弱，沖擊波初始峰值的分布規(guī)律沒有因此改變。類似的情況在飛片質(zhì)量為1.325 kg的工況條件下同樣出現(xiàn)。綜合圖5（a）和圖5（b）可以發(fā)現(xiàn)，當飛片質(zhì)量確定時，活塞厚度的改變不會改變沖擊波初始峰值隨飛片速度變化的整體趨勢。

圖5 不同質(zhì)量飛片撞擊活塞得到的沖擊波峰值載荷Fig.5 Relationship between shock wave strength and piston length

4.3 飛片對水下爆炸沖擊波衰減的影響

飛片質(zhì)量分別為0.265 kg、0.795 kg、1.325 kg和1.855 kg時，飛片在0～250 m/s速度范圍內(nèi)正撞擊12 mm厚活塞產(chǎn)生的水下爆炸初始沖擊波沿非藥式水下爆炸裝置中軸線的衰減趨勢如圖6所示，其中P0分別為對應(yīng)各工況的初始沖擊波峰值，將各高斯觀測點測得的初始沖擊波進行無量綱化。從圖中可見，各初始沖擊波隨時間呈指數(shù)型衰減。當一定質(zhì)量的飛片以50 m/s撞擊活塞時，水下爆炸初始沖擊波傳播0.23 ms后，無量綱化的沖擊波強度趨近于零，這與其它速度條件下的沖擊波衰減趨勢存在差異，但隨著飛片撞擊速度的增加，初始沖擊波的衰減趨勢趨于穩(wěn)定。這說明該質(zhì)量范圍內(nèi)的飛片在初始撞擊速度為0～250 m/s范圍內(nèi)時，飛片的初始撞擊速度并不影響水下非藥式爆炸初始沖擊波隨時間的衰減趨勢。同時從圖6中可見，飛片的質(zhì)量影響了初始沖擊波的衰減速度。隨著飛片質(zhì)量的增加，初始沖擊波隨時間的衰減趨勢逐漸減緩。

4.4 活塞對水下爆炸沖擊波衰減的影響

圖7所示為不同質(zhì)量飛片在0～250 m/s范圍內(nèi)正撞擊不同厚度活塞產(chǎn)生的初始沖擊波衰減時間常數(shù)變化趨勢。由圖7可見，在該活塞厚度變化范圍內(nèi)，不同質(zhì)量飛片所對應(yīng)的初始沖擊波衰減時間常數(shù)隨活塞質(zhì)量的變化規(guī)律基本相同，初始沖擊波的衰減時間常數(shù)隨著活塞厚度的增加而逐漸減小，但并不改變初始沖擊波衰減時間常數(shù)隨飛片質(zhì)量增加而增加的整體趨勢。

圖6 不同質(zhì)量飛片正撞擊12 mm厚活塞產(chǎn)生的初始沖擊波隨時間衰減趨勢Fig.6 Variation of the dimensionless shock loading with time

4.5 計算結(jié)果分析

基于上述數(shù)值模擬結(jié)果，可以發(fā)現(xiàn)，非藥式水下爆炸初始沖擊波峰值只與飛片的撞擊速度有關(guān)。根據(jù)聲學近似理論，非藥式水下爆炸初始沖擊波強度可表示為：

其中：ρw為水的密度，vw為水中聲速。

設(shè)活塞與水的接觸位置為x=0處，如圖1所示，當活塞與水作用瞬間，在x=0位置有：

其中：vp為飛片撞擊活塞后，活塞獲得的初始速度，vw為水在x=0位置處的粒子速度。圖8為通過設(shè)置在活塞上的高斯觀測點得到的活塞初始速度隨飛片撞擊速度變化曲線。由該圖可見，飛片速度vf與活塞速度vp成線性關(guān)系，擬合得到：

圖7 初始沖擊波衰減時間常數(shù)與活塞質(zhì)量的關(guān)系Fig.7 Relation between decay constant and piston mass

圖8 活塞運動速度隨飛片初始撞擊速度變化曲線Fig.8 Relations between piston velocity and flyer velocity

圖9 f1 （mf ）與飛片質(zhì)量的關(guān)系Fig.9 Relations between f1 （mf ）and flyer mass

其中：k=1.01。

將（2）、（3）式代入（1）式得到非藥式水下爆炸沖擊波初始峰值與飛片速度的關(guān)系為：

如圖4、5所示，理論預(yù)測值與仿真實驗得到的非藥式水下爆炸沖擊波初始峰值吻合良好，當活塞厚度為12 mm時，理論值略小于仿真實驗值，但隨著活塞厚度的增加，該差異逐漸縮小。

基于上述數(shù)值模擬結(jié)果，發(fā)現(xiàn)非藥式水下爆炸沖擊波衰減時間常數(shù)與飛片質(zhì)量、活塞質(zhì)量、水中聲速、水的密度以及活塞直徑有關(guān)，推導出非藥式水下爆炸沖擊波初始峰值的衰減時間常數(shù)與飛片和活塞之間的關(guān)系，非藥式水下爆炸沖擊波衰減時間常數(shù)可表述為：

由圖7可知，在本文所做實驗工況范圍內(nèi)，沖擊波衰減時間常數(shù)與活塞質(zhì)量成線性關(guān)系，即：

其中：f1（mf）為飛片質(zhì)量mf的函數(shù)，a為待定系數(shù)，mp為活塞質(zhì)量。由圖7可見，當飛片質(zhì)量分別為0.265 kg、0.795 kg和 1.855 kg時，擬合得到的線性關(guān)系斜率近似相同，而飛片質(zhì)量為1.325 kg時擬合直線的斜率與其他三種工況有一定差異，故將擬合值平均得到a的值為311.89。代入（6）式得：

圖9為由（7）式得到的f1（mf）隨飛片質(zhì)量mf變化而產(chǎn)生的變化規(guī)律。 由圖9可見，f1（mf）的一般形式可以表示為：

其中：b，c為待定系數(shù)，擬合得到：

將（9）式代入（7）式得到預(yù)測非藥式水下爆炸沖擊波初始峰值衰減時間常數(shù)的經(jīng)驗公式：

5 結(jié) 論

針對非藥式水下爆炸沖擊波模擬實驗這一特定問題，對具有不同速度、不同質(zhì)量的飛片及不同厚度的活塞進行了數(shù)值模擬，結(jié)果表明：非藥式水下爆炸沖擊波初始峰值與飛片的質(zhì)量和活塞的厚度無關(guān)，與飛片初始速度成線性關(guān)系；飛片的初始撞擊速度對初始沖擊波的衰減時間常數(shù)無顯著影響，其僅隨飛片和活塞的質(zhì)量改變而改變。通過仿真實驗數(shù)據(jù)進行回歸分析，得到了預(yù)測非藥式水下爆炸沖擊波初始峰值和初始沖擊波衰減時間常數(shù)的經(jīng)驗公式。

[1]Despande V S,Heaver A,Fleck N A.An underwater shock simulator[J].Proc.R.Soc.A,2006,462:1021-1041.

[2]Espinosa H D,Lee S,Moldovan N.A novel fluid structure interaction experiment to investigate deformation of structural elements subjected to impulsive loading[J].Experimental Mechanics,2006,46:805-824.

[3]Mori L F,Queheillalt D T,Wadley H N G,Espinosa H D.Deformation and failure modes of I-core sandwich structures subjected to underwater impulsive loads[J].Experimental Mechanics,2009,49:257-275.

[4]Latourte F,Gregoire D,Zenkert D,et al.Failure mechanisms in composite panels subjected to underwater impulsive loads[J].Journal of the Mechanics and Physics of Solids,2011,59:1623-1646.

[5]Huson P,Asaro R J,Stewart L,Hegemier G A.Non-explosive methods for simulating blast loading of structures with complex geometries[J].International Journal of Impact Engineering,2010,38:546-557.

[6]Li Bosong,He Yongsheng,Han Nairen,Du Kui.Brief introduction of shock testing machine for simulating blast vibration[J].Explosion and Shock Waves,2002,22(1):79-82.(in Chinese)

[7]Ren Peng,Zhang Wei,Guo Zitao,Wei Gang.Numerical simulation for deformation of multi-layer steel plates under underwater impulsive loading[J].Journal of Harbin Institute of Technology,2012,19(2):99-103.

[8]Ren Peng,Zhang Wei,Huang Wei,et al.Deformation mode and strain field analysis of clamped air-back circular plate subjected to underwater explosive loading[J].Journal of Ship Mechanics,2013,17(11):1339-1344.

[9]Ren Peng,Zhang Wei,Huang Wei,et al.Research on non-explosive underwater shock loading device[J].Explosion and Shock Wave,2014,34(3):334-339.(in Chinese)

[10]Xu Yuxin,Wang Shushan,Li Yuan.Study on numerical simulation of the underwater explosive[J].Journal of Projectiles,Rockets,Missiles and Guidance,2009,29(6):95-97.(in Chinese)