全入射角度下平板沖擊波的壁壓載荷及 局部空化特性*

李海濤，朱石堅(jiān)，陳志堅(jiān)，牟金磊

(1.海軍工程大學(xué)科研部，湖北 武漢 430033； 2.海軍工程大學(xué)艦船工程系，湖北 武漢 430033)

當(dāng)水下爆炸沖擊波傳播到結(jié)構(gòu)邊界時(shí)，會(huì)產(chǎn)生反射、透射等現(xiàn)象，結(jié)構(gòu)物的存在造成沖擊波載荷特性發(fā)生變化。結(jié)構(gòu)邊界處沖擊波壁壓載荷的確定是研究水下爆炸作用下結(jié)構(gòu)早期響應(yīng)的前提條件，可為研究后續(xù)爆炸氣泡運(yùn)動(dòng)引起的結(jié)構(gòu)響應(yīng)提供初始運(yùn)動(dòng)參數(shù)；同時(shí)，壁壓載荷的確定也可為數(shù)值模擬及理論計(jì)算提供相對(duì)準(zhǔn)確的載荷輸入條件。

對(duì)平板結(jié)構(gòu)的壁壓載荷及結(jié)構(gòu)初期響應(yīng)，主要利用平面波近似法對(duì)簡(jiǎn)單結(jié)構(gòu)進(jìn)行理論求解[1-3]，而對(duì)有限尺度平板結(jié)構(gòu)的準(zhǔn)確響應(yīng)解析模型尚未充分研究。對(duì)于有限尺度結(jié)構(gòu)處的沖擊波壁壓載荷，金輝等[4]、程素秋等[5]、樊寶順等[6]對(duì)實(shí)船及艙段模型的結(jié)構(gòu)沖擊波壁壓進(jìn)行了初步的實(shí)驗(yàn)研究，比較了壁壓載荷曲線與自由場(chǎng)壓力曲線的變化規(guī)律，但未就壁壓載荷衰減變化特性以及與入射角度之間的關(guān)系等問題進(jìn)行深入的探討。

本文中以平面波理論為基礎(chǔ)，推導(dǎo)無(wú)限平板在全入射角度下的沖擊波壁壓載荷理論計(jì)算公式，利用前期沖擊波壁壓實(shí)驗(yàn)數(shù)據(jù)對(duì)該理論公式進(jìn)行修正，提出一種壁壓修正方法和經(jīng)驗(yàn)公式，深入分析斜入射角度下壁壓載荷的變化特性，并根據(jù)壁壓的負(fù)壓特性對(duì)平板局部空化進(jìn)行初步研究，為后續(xù)研究局部空化對(duì)結(jié)構(gòu)響應(yīng)的影響奠定基礎(chǔ)。

1 平板壁壓計(jì)算公式

圖1 任意角度入射沖擊波作用下的無(wú)限平板Fig.1 Infinite plate subjected to an inclined plane shockwave

將艦船結(jié)構(gòu)在沖擊波作用下的早期響應(yīng)過程簡(jiǎn)化為無(wú)限平板在平面波作用下的響應(yīng)過程，平板兩側(cè)分別為水和空氣介質(zhì)，整個(gè)簡(jiǎn)化物理模型如圖1所示。對(duì)于無(wú)限平板結(jié)構(gòu)，當(dāng)任意角度入射沖擊波pi作用到平板邊界時(shí)，會(huì)在平板邊界面上形成反射波pr和透射波pt，入射波和反射波疊加即為結(jié)構(gòu)流固耦合面處的壁壓載荷。在不同介質(zhì)中傳播的沖擊波滿足平面波反射和折射定理。

理論計(jì)算過程中，進(jìn)一步引入如下假設(shè)：

(1)水為不可壓縮、無(wú)旋、無(wú)黏的理想流體；

(2)與平板相鄰的水質(zhì)點(diǎn)和平板具有相同的運(yùn)動(dòng)速度；

(3)忽略壓縮波對(duì)平板變形的影響，平板兩側(cè)邊界質(zhì)點(diǎn)速度相同；

(4)假設(shè)不同入射角度沖擊波同時(shí)到達(dá)平板，忽略因入射沖擊波作用時(shí)間差異而引起的平板運(yùn)動(dòng)參數(shù)變化。

圖1中α為入射角和反射角，β為透射角，ρa(bǔ)、ρw和ca、cw為空氣、水的密度和聲速。根據(jù)連續(xù)性邊界條件，在斜入射條件下，平板兩側(cè)質(zhì)點(diǎn)的法向速度必須相等。入射面一側(cè)質(zhì)點(diǎn)速度為(pi-pr)cosα/(ρwcw)，透射面一側(cè)質(zhì)點(diǎn)速度為ptcosβ/(ρa(bǔ)ca)，因此

vp=(pi-pr)cosα/(ρwcw)=ptcosβ/(ρa(bǔ)ca)

(1)

(2)

根據(jù)平面波折射定理，入射角度α和透射角β之間滿足如下方程：

(3)

方程(1)～(2)中，pr、pt和β為未知量，其他均為已知量，聯(lián)立求解得斜入射條件下平板邊界理論壁壓

(4)

2 壁壓公式的實(shí)驗(yàn)修正

有學(xué)者對(duì)某特定實(shí)船及結(jié)構(gòu)模型的沖擊波壁壓峰值進(jìn)行了測(cè)量及數(shù)據(jù)分析，獲得了一些斜入射條件下沖擊波壁壓實(shí)驗(yàn)值。以該實(shí)驗(yàn)數(shù)據(jù)為基礎(chǔ)，同時(shí)結(jié)合小尺度平板模型壁壓實(shí)驗(yàn)數(shù)據(jù)，可以對(duì)已建立的斜入射條件下無(wú)限平板沖擊波壁壓計(jì)算公式進(jìn)行初步修正。

2.1 壁壓峰值修正

當(dāng)有限尺度結(jié)構(gòu)受水下爆炸沖擊波作用時(shí)，實(shí)際壁壓值比理論計(jì)算值偏小，需引入壁壓修正系數(shù)Cp(Cp為壁壓實(shí)驗(yàn)峰值pwe和理論峰值pwt之比)。本文中作者前期對(duì)垂直入射條件下沖擊波壁壓峰值修正工作開展了初步研究，在系列小尺寸鋼板爆炸沖擊實(shí)驗(yàn)基礎(chǔ)上，獲得了部分工況條件下壁壓峰值修正系數(shù)Cp的實(shí)驗(yàn)值，初步總結(jié)出經(jīng)驗(yàn)公式如下[7]：

(5)

式中：me為炸藥當(dāng)量，R為炸藥與測(cè)點(diǎn)的距離，系數(shù)K=0.614 2、φ1=-0.072 2、φ2=0.007 8。

針對(duì)斜入射條件，擬在垂直入射壁壓修正系數(shù)公式基礎(chǔ)上，進(jìn)一步考慮入射角度對(duì)修正系數(shù)的影響，將全入射角度沖擊波壁壓峰值修正系數(shù)Cp表示為如下形式：

(6)

式中：α為入射角度。當(dāng)α= 0時(shí)，f(α)=1，此時(shí)公式(6)即為垂直入射時(shí)的平板壁壓峰值修正公式。

利用前期壁壓實(shí)驗(yàn)數(shù)據(jù)對(duì)函數(shù)f(α)進(jìn)行擬合分析，確定其函數(shù)形式，具體實(shí)驗(yàn)工況及數(shù)據(jù)見表1。其中，工況1～3來自于本文中作者新完成的1m尺度方板水下爆炸實(shí)驗(yàn)；工況4～5來源于文獻(xiàn)[4]；工況6～7來源于文獻(xiàn)[5]；工況8～15來源于文獻(xiàn)[6]。

通過對(duì)表1中不同工況下壁壓實(shí)驗(yàn)數(shù)據(jù)的分析，確定f(α)的函數(shù)形式為：

f(α)=2.289 4cos3α-4.062 2cos2α+2.900 4cosα-0.090 718

(7)

表1 不同爆炸工況下沖擊波壁壓峰值實(shí)驗(yàn)值及理論值對(duì)比Table 1 Comparison of theoretical and experimental results of wall pressure in test cases

注：工況1～3、5～8、13～15為函數(shù)f(α)擬合樣本點(diǎn)工況；工況4、9～12為考查點(diǎn)工況。

圖2 函數(shù)f(α)計(jì)算值與實(shí)驗(yàn)值的比較Fig.2 Comparison of theoretical and experimental results of f(α)

圖2中對(duì)函數(shù)f(α)擬合值及實(shí)驗(yàn)值進(jìn)行了比較。從圖中可以看出，通過實(shí)驗(yàn)數(shù)據(jù)擬合得到的不同斜入射角度下的f(α)函數(shù)曲線，對(duì)于考查實(shí)驗(yàn)工況數(shù)據(jù)也有較好的預(yù)報(bào)精度；隨著入射角度的增大，壁壓峰值修正系數(shù)值非線性減小，這與實(shí)驗(yàn)數(shù)據(jù)反映的規(guī)律基本一致；當(dāng)沖擊波入射角度較小時(shí)，該擬合曲線預(yù)報(bào)精度較高，而對(duì)于入射角度接近于90°的情況，本函數(shù)提供的計(jì)算數(shù)據(jù)會(huì)存在一定的預(yù)報(bào)誤差。

2.2 衰減常數(shù)的選取

實(shí)驗(yàn)研究發(fā)現(xiàn)，沖擊波壁壓曲線在衰減過程中也基本滿足指數(shù)衰減規(guī)律，但與自由場(chǎng)壓力相比，壁壓曲線的衰減速度更快。表2中列出了不同工況下理論壁壓衰減常數(shù)θpt與實(shí)驗(yàn)壁壓衰減常數(shù)θpe的對(duì)比情況，表中ε表示誤差。從表中數(shù)據(jù)可以看出，對(duì)于壁壓衰減常數(shù)而言，理論值和實(shí)驗(yàn)值相差較小，均在10%以內(nèi)，可基本認(rèn)為實(shí)驗(yàn)壁壓曲線的早期衰減過程與理論壁壓曲線一致。在目前實(shí)驗(yàn)數(shù)據(jù)偏少的情況下，有限尺度平板處壁壓衰減常數(shù)可以參考無(wú)限平板壁壓衰減常數(shù)進(jìn)行計(jì)算。

表2 不同爆炸工況下沖擊波衰減常數(shù)實(shí)驗(yàn)值及理論值對(duì)比Table 2 Comparison of theoretical and experimental results of θ in test cases

2.3 壁壓修正公式

(8)

式中:θpt為由式(4)計(jì)算得到的沖擊波壁壓衰減常數(shù)。

3 分析與討論

3.1 壁壓載荷特性

為分析有限尺度平板沖擊波壁壓載荷實(shí)際變化特性，針對(duì)尺寸為1 m×1 m、厚1 mm的鋼質(zhì)平板開展了水下爆炸實(shí)驗(yàn)研究。實(shí)驗(yàn)時(shí)，所取沖擊波入射角度較小，基本為垂直入射狀態(tài)。圖3給出了表1中工況2、3條件下沖擊波壁壓曲線理論值和實(shí)驗(yàn)值的比較情況。

圖3 沖擊波壁壓理論值和實(shí)驗(yàn)值的比較Fig.3 Comparison of theoretical and experimental results of wall pressure

圖4 相同沖擊因子下入射角度對(duì) 壁壓曲線的影響Fig.4 Effect of incidence angle on wall-pressure value

圖5 入射角度對(duì)同一點(diǎn)處沖擊波 壓力衰減常數(shù)的影響Fig.5 Effect of incidence angle on θ

從圖3可以看出，沖擊波壁壓衰減的前期過程基本滿足指數(shù)衰減規(guī)律，且衰減速度很快；衰減后期階段，理論和實(shí)驗(yàn)壁壓曲線均出現(xiàn)負(fù)壓值，這是反射稀疏波與入射波疊加的效果；在沖擊波前期衰減過程中，經(jīng)指數(shù)函數(shù)修正后的壁壓值與實(shí)驗(yàn)值吻合較好，但因指數(shù)函數(shù)本身特點(diǎn)，修正壁壓值不會(huì)出現(xiàn)負(fù)壓情況，這與實(shí)驗(yàn)壁壓曲線不一致，而無(wú)限平板理論壁壓公式能夠較好地反映沖擊波壁壓衰減后期的變化特性。因此，指數(shù)函數(shù)壁壓修正曲線僅適用于前期衰減過程，對(duì)于后期階段，仍需要通過理論壁壓公式分析其后續(xù)變化規(guī)律。

為了分析入射角度對(duì)沖擊波壁壓曲線的影響特性，選定50 g TNT炸藥在3 mm厚平板下方爆炸為研究工況，且保持沖擊因子為0.112 kg0.5/m不變。從圖4可以看出，針對(duì)平板某一入射點(diǎn)而言，隨著入射角度的增大，沖擊波壁壓峰值雖然保持不變，但其壁壓曲線衰減過程加快，且衰減速度與入射角度變化過程存在非線性關(guān)系；入射角度越小，沖擊波壁壓曲線出現(xiàn)最低負(fù)壓值的時(shí)間越晚，且最低負(fù)壓絕對(duì)值越大，該處因負(fù)壓存在而形成局部空化的能力增強(qiáng)。

為進(jìn)一步分析全入射角度下沖擊波壁壓曲線衰減常數(shù)的變化特點(diǎn)，選定50 g TNT炸藥在3 mm厚平板下方爆炸為研究工況，且保持沖擊因子為0.224 kg0.5/m不變。從圖5可以看出，對(duì)于同一爆源，垂直入射處的沖擊波壁壓曲線衰減時(shí)間最長(zhǎng)；隨著入射角度的增大，其壁壓衰減常數(shù)減小；而對(duì)同一入射點(diǎn)處的自由場(chǎng)壓力曲線而言，入射角度越大(即考查點(diǎn)爆距越大)，壓力衰減過程越長(zhǎng)；對(duì)于所選計(jì)算工況而言，垂直入射時(shí)的壁壓衰減常數(shù)約為自由場(chǎng)壓力曲線的1/5，隨著入射角度的增大，兩種壓力曲線的衰減常數(shù)值相差越大。

3.2 負(fù)壓及空化特性

從壁壓載荷特性研究可以看出，沖擊波壁壓衰減后期會(huì)出現(xiàn)負(fù)壓值，隨著爆炸工況的變化，最低負(fù)壓值的大小及出現(xiàn)時(shí)間均會(huì)發(fā)生變化。選定厚度τ=15mm平板下方垂直爆距H=10 m的工況為研究對(duì)象，進(jìn)一步分析藥量me及入射角度α變化對(duì)最低負(fù)壓值pwn的影響特性。從圖6可以看出，隨著入射角度的增大，沖擊波最低負(fù)壓的絕對(duì)值非線性減小，當(dāng)α<20°時(shí)，最低負(fù)壓值變化緩慢，當(dāng)20°<α<70°時(shí)，最低負(fù)壓值存在快速變化過程，當(dāng)α>70°時(shí)，最低負(fù)壓值的變化又趨于平緩；同一入射角度條件下，隨著藥量的增加，最低負(fù)壓的絕對(duì)值是增大的，說明大藥量更容易造成平板處發(fā)生空化現(xiàn)象，但藥量的增加對(duì)造成最低負(fù)壓絕對(duì)值增大的效果是逐步減小的，最低負(fù)壓存在一個(gè)極限值。

當(dāng)平板沖擊波壁壓出現(xiàn)負(fù)壓時(shí)，在平板邊界面上可能會(huì)出現(xiàn)局部空化現(xiàn)象。文獻(xiàn)[8]在分析過程中，取水的空化壓力為-3.0MPa，本文中也選取該值作為空化臨界壓力，對(duì)典型爆炸工況下無(wú)限平板局部空化的形成范圍進(jìn)行初步研究。

由于幾何對(duì)稱性，炸藥爆炸時(shí)會(huì)在平板上形成圓形的局部空化區(qū)域。選定500kgTNT炸藥在無(wú)限平板下方爆炸作為分析工況，研究不同平板厚度時(shí)空化區(qū)域半徑Rc與垂直爆距H之間的變化規(guī)律。從圖7可以看出，對(duì)于同一厚度平板，隨著垂直爆距H的增大，平板上形成局部空化的面積會(huì)出現(xiàn)先增大后減小的過程，且當(dāng)H超過某一值時(shí)，整個(gè)平板會(huì)因?yàn)樽畹拓?fù)壓值達(dá)不到空化條件而無(wú)法出現(xiàn)局部空化；隨著平板厚度τ的增加，同一爆距條件下，爆炸形成局部空化的能力增強(qiáng)，空化區(qū)域范圍也將增大，這說明板厚對(duì)空化的形成也有重要影響：板厚越大，壁壓的最低負(fù)壓絕對(duì)值越大，平板在更大爆距范圍內(nèi)出現(xiàn)空化現(xiàn)象的可能性增大。

圖6 不同藥量條件下入射角度對(duì) 最小壁壓值的影響Fig.6 Effect of incidence angle on minimum wall-pressure value

圖7 不同板厚條件下垂直爆距對(duì) 局部空化區(qū)域半徑的影響Fig.7 Effect of minimum stand-off distance on radius of local cavitation zone

4 結(jié) 論

(1)全入射角度下有限尺度平板結(jié)構(gòu)沖擊波壁壓公式僅需要在無(wú)限平板理論壁壓公式基礎(chǔ)上進(jìn)行峰壓修正，修正后的壁壓曲線與實(shí)際壁壓曲線吻合較好；隨著入射角度的增加，壁壓曲線修正系數(shù)非線性減小。

(2)相同沖擊因子條件下，隨著入射角度的增大，沖擊波壁壓形成最低負(fù)壓的時(shí)間越早，其負(fù)壓絕對(duì)值也相對(duì)減小，形成局部空化的能力降低。

(3)相同條件下，隨著藥量或者平板厚度的增大，壁壓的最低負(fù)壓絕對(duì)值也增大，形成局部空化的能力增強(qiáng)；炸藥在某一垂直爆距范圍內(nèi)變化時(shí)，局部空化范圍會(huì)出現(xiàn)先增大后減小的過程，超過該爆距范圍，整個(gè)平板會(huì)因?yàn)樽畹拓?fù)壓值達(dá)不到空化條件而無(wú)法出現(xiàn)局部空化。

[1] Jiang J, Olson M D.Rigid-plastic analysis of underwater blast loaded stiffened plates[J].International Journal of Mechanic and Science, 1995,37(8):843-859.

[2] Makinen K.Cavitation models for structures excited by a plane shock wave[J].Journal of Fluids and Structures, 1998,12:85-101.

[3] Santiago L D.Fluid-interaction and cavitation effects on a surface ship model due to an underwater explosion[R].California: Naval Postgraduate School, 1997.

[4] 金輝,周學(xué)濱,周華,等.水下爆炸中自由場(chǎng)壓力和船體壁壓的測(cè)量與分析[J].海軍工程大學(xué)學(xué)報(bào)，2009,21(5):82-87.Jin Hui, Zhou Xue-bin, Zhou Hua, et al.Measurement and analysis of free-field pressure and ship hull pressure under water explosion[J].Journal of Naval University of Engineering, 2009,21(5):82-87.

[5] 程素秋,樊寶順,薛飛,等.水下非接觸爆炸作用下艙段模型的動(dòng)態(tài)響應(yīng)[J].爆炸與沖擊，2008,28(4):360-366.Cheng Su-qiu, Fan Bao-shun, Xue Fei, et al.Dynamical response measurement of a cabin model subjected to noncontact underwater explosion[J].Explosion and Shock Waves, 2008,28(4):360-366.

[6] 樊寶順,程素秋,韓峰.艦船艙段模型在水下爆炸作用下的壁壓分析[J].中國(guó)艦船研究,2009,4(5):20-22.Fan Bao-shun, Cheng Su-qiu, Han Feng.Wall pressure analysis of cabin model subjected to underwater explosion[J].Chinese Journal of Ship Research, 2009,4(5):20-22.

[7] 李海濤,朱錫,張振華.水下爆炸球面沖擊波作用下船體梁的剛塑性動(dòng)響應(yīng)特性[J].工程力學(xué)，2010,27(10):202-207.Li Hai-tao, Zhu Xi, Zhang Zhen-hua.Dynamic rigid-plasitc response of ship-like beam subjected to underwater spherical shockwaves[J].Engineering Mechanics, 2010,27(10):202-207.

[8] 顧文彬,蘇青笠,劉建青,等.水下平面爆炸沖擊波作用下空化區(qū)域形成及其特征[J].爆破,2004，21(4)：8-11.Gu Wen-bin, Su Qing-li, Liu Jian-qiang, et al.Formation process of cavitation region affected by underwater plane blast shock wave and its characteristics[J].Blasting, 2004,21(4):8-11.