劉籌資，查恩堯，歐陽(yáng)和平，程永舟，郭錦澤，鄭長(zhǎng)青

(1.長(zhǎng)沙理工大學(xué)水利工程學(xué)院，長(zhǎng)沙 410114；2.中交第四航務(wù)工程勘察設(shè)計(jì)院有限公司，廣州 510000；3.湖南水利水電職業(yè)技術(shù)學(xué)院，長(zhǎng)沙 410100；4.珠海爆破新技術(shù)開發(fā)公司，廣東 珠海 519099)

水下爆炸技術(shù)作為一種快速有效的水下施工手段，在航道疏浚[1]、水下巖石開挖[2]、水下結(jié)構(gòu)拆除[3-4]等工程中應(yīng)用廣泛。水下爆炸在不同條件下會(huì)形成多階段的物理現(xiàn)象，當(dāng)水下爆炸發(fā)生在水面附近時(shí)，除沖擊波傳播、氣泡脈動(dòng)現(xiàn)象外，氣泡在脈動(dòng)過(guò)程中還會(huì)與水面相互作用在水面附近形成形態(tài)不同的水冢[5]。在重力作用下，水冢沖擊水面形成的波浪，會(huì)對(duì)周圍結(jié)構(gòu)特別是艦船等浮體結(jié)構(gòu)產(chǎn)生較大影響。因此開展水下爆炸氣泡運(yùn)動(dòng)對(duì)水面形態(tài)變化影響的研究，對(duì)深入了解近水面水下爆炸現(xiàn)象，指導(dǎo)水下爆炸實(shí)施時(shí)對(duì)周邊結(jié)構(gòu)物的防護(hù)具有重要意義。

在爆炸興波方面，沈國(guó)光等[6-7]基于MAC算法將爆炸氣泡與自由液面的相互作用過(guò)程簡(jiǎn)化，這一方法雖然能夠簡(jiǎn)化水下爆炸的過(guò)程，但賦予氣泡的內(nèi)壓是以與試驗(yàn)水面現(xiàn)象的異同作為標(biāo)準(zhǔn)的，因此只能作為特定裝藥量的計(jì)算，此外由于模擬中忽略了氣泡脈動(dòng)、射流的過(guò)程，對(duì)于興波機(jī)理的討論是不充分的。Safiyari等[8-9]基于SPH方法，將興波過(guò)程簡(jiǎn)化為空腔的擴(kuò)張、坍縮，進(jìn)而引起水面波動(dòng)的過(guò)程，這一研究同樣忽略了氣泡脈動(dòng)、射流現(xiàn)象在水面形態(tài)形成中的關(guān)鍵作用。王志凱[10]基于ALE算法研究了水下爆炸氣泡破碎興波對(duì)浮動(dòng)沖擊平臺(tái)的影響，該研究中雖然考慮了氣泡脈動(dòng)，但其計(jì)算工況中，水深距離參數(shù)達(dá)到了2.7，依據(jù)文獻(xiàn)[11-13]，該條件下氣泡與自由液面相互作用形成的水冢對(duì)自由液面的擾動(dòng)很小，能夠形成的興波波浪也較小，并且由于模型設(shè)置了浮體結(jié)構(gòu)，水冢的發(fā)展以及興波波浪的形成也會(huì)被抑制。劉云龍等[12]基于BEM方法也討論了氣泡破碎興波對(duì)浮體結(jié)構(gòu)的影響，但僅限于氣泡脈動(dòng)成簡(jiǎn)化的錐形水冢沖擊水面形成的波浪與浮體結(jié)構(gòu)的相互作用。李汪諱等[14]、Wang等[15]基于BEM方法對(duì)水下爆炸興波現(xiàn)象進(jìn)行了模擬，但文獻(xiàn)[14]未討論各種水冢形態(tài)對(duì)水面興波的影響，文獻(xiàn)[15]簡(jiǎn)化了水冢形態(tài)，忽略不連續(xù)水體，將氣泡運(yùn)動(dòng)在水面附近形成的水冢分別簡(jiǎn)化為水下空腔和水面水柱兩種形態(tài)。前人在興波方面的研究由于采用的BEM方法對(duì)不連續(xù)邊界的模擬存在困難，基本都建立在一定的水冢形態(tài)簡(jiǎn)化條件下進(jìn)行，這與實(shí)際的水面現(xiàn)象不相符，對(duì)于水面波動(dòng)發(fā)展過(guò)程的討論也存在一定的誤差。

基于此，本數(shù)值模擬研究采用LS-DYNA非線性有限元軟件，基于ALE算法，結(jié)合水冢形成過(guò)程及其特征幾何參數(shù)、氣泡運(yùn)動(dòng)特性與水面波動(dòng)過(guò)程，就不同水深距離參數(shù)下水下爆炸氣泡運(yùn)動(dòng)對(duì)水面形態(tài)的影響機(jī)理開展研究。

1 氣泡運(yùn)動(dòng)數(shù)值模型的建立

1.1 算法選擇與幾何模型的建立

由于炸藥爆炸的氣體產(chǎn)物在水下形成的氣泡存在著膨脹收縮過(guò)程，在氣泡運(yùn)動(dòng)的影響下水面會(huì)存在復(fù)雜的變形，同時(shí)為了捕捉氣泡的運(yùn)動(dòng)特性，對(duì)于氣泡邊界的捕捉也具有一定的要求。因此，本研究采用ALE算法進(jìn)行三維數(shù)值模擬。模型中涉及到的炸藥、水和空氣均采用ALE算法。

為了兼顧計(jì)算效率，采用三維1/4模型，在對(duì)稱面上施加對(duì)稱邊界條件，在空氣域和水域外邊界處設(shè)置pressure outflow邊界單元層，以模擬無(wú)限邊界條件，模型的x、y方向和裝藥下方的水域深度應(yīng)至少為氣泡脈動(dòng)最大半徑的4倍，空氣層的高度則需滿足水冢中心水柱向上發(fā)展的要求。為了提高計(jì)算效率同時(shí)保證模擬結(jié)果的正確性，在氣泡最大半徑1.5倍范圍內(nèi)采用與裝藥半徑尺寸相同的均勻網(wǎng)格，其他區(qū)域采用由均勻網(wǎng)格區(qū)域至流域外邊界逐漸漸變的網(wǎng)格，漸變網(wǎng)格最小處與均勻網(wǎng)格相同，最大處網(wǎng)格尺寸為均勻網(wǎng)格的10倍。幾何模型如圖1所示。

圖1 幾何模型Fig.1 Geometric model

1.2 材料模型的設(shè)置

模型主要涉及炸藥、水和空氣3種物質(zhì)，其中炸藥采用JWL狀態(tài)方程進(jìn)行描述，水采用Grüneisen狀態(tài)方程進(jìn)行描述，空氣采用線性多項(xiàng)式狀態(tài)方程進(jìn)行描述。各材料的具體參數(shù)如表1～表3所示[16-17]。

表1 TNT炸藥的材料方程和狀態(tài)方程參數(shù)取值

表2 水的材料方程和狀態(tài)方程參數(shù)取值

表3 空氣的材料方程和狀態(tài)方程參數(shù)取值

1.3 氣泡運(yùn)動(dòng)模型有效性驗(yàn)證

為了驗(yàn)證氣泡運(yùn)動(dòng)模型的有效性，設(shè)置與文獻(xiàn)[18]實(shí)驗(yàn)相同的工況，使用5.2 g TNT的球形裝藥(與實(shí)驗(yàn)中4 g PENT炸藥等效)，裝藥半徑為0.009 132 m，模型計(jì)算尺寸為1 m×1 m×4 m，漸變網(wǎng)格劃分尺寸最小為0.009 m，最大為0.1 m，均勻網(wǎng)格處網(wǎng)格尺寸為0.009 m。

由自由場(chǎng)條件下文獻(xiàn)[18]中實(shí)驗(yàn)結(jié)果和數(shù)值模擬結(jié)果的對(duì)比(見圖2)可以發(fā)現(xiàn)，在自由場(chǎng)條件下數(shù)值模擬的氣泡運(yùn)動(dòng)趨勢(shì)與文獻(xiàn)中的各個(gè)階段保持一致；由在近自由液面條件下文獻(xiàn)[18]中實(shí)驗(yàn)結(jié)果和數(shù)值模擬結(jié)果進(jìn)行的對(duì)比(見圖3)可知，兩者氣泡運(yùn)動(dòng)趨勢(shì)大致一致。

圖2 自由場(chǎng)條件下氣泡運(yùn)動(dòng)趨勢(shì)對(duì)比Fig.2 Comparison of bubble movement trends under free field conditions

圖3 近自由液面條件下氣泡運(yùn)動(dòng)趨勢(shì)對(duì)比Fig.3 Comparison of bubble movement trends under the condition of near free liquid surface

對(duì)于水下爆炸氣泡脈動(dòng)的第一周期及氣泡半徑計(jì)算，已經(jīng)形成了能夠應(yīng)用于工程實(shí)際的經(jīng)驗(yàn)公式[19]：

(1)

(2)

式中：R為氣泡第一次脈動(dòng)的氣泡最大半徑，m；T為脈動(dòng)周期，s；W為裝藥當(dāng)量，kg；D為裝藥位置處水深，m；10.3為大氣壓強(qiáng)對(duì)應(yīng)液柱高度，m；K4、K5為與炸藥相關(guān)的計(jì)算參數(shù)。

常見氣泡脈動(dòng)參數(shù)取值如表4所示[19]。

表4 氣泡脈動(dòng)經(jīng)驗(yàn)公式計(jì)算參數(shù)

另外也對(duì)自由場(chǎng)和近水面條件下的氣泡半徑時(shí)程進(jìn)行了分析，數(shù)值模擬結(jié)果數(shù)據(jù)與文獻(xiàn)[18]中的工況數(shù)據(jù)進(jìn)行對(duì)比分析發(fā)現(xiàn)兩者之間的誤差在10%以內(nèi)(見表5)，為了體現(xiàn)研究的嚴(yán)謹(jǐn)性又對(duì)經(jīng)驗(yàn)公式[19]的結(jié)果進(jìn)行了計(jì)算對(duì)比，分析發(fā)現(xiàn)在兩種條件下氣泡的半徑和周期變化誤差均在10%以內(nèi)(見表5)，說(shuō)明氣泡運(yùn)動(dòng)數(shù)值模型與實(shí)驗(yàn)結(jié)果一致性良好。

表5 氣泡半徑與周期誤差計(jì)算結(jié)果

1.4 網(wǎng)格尺寸對(duì)計(jì)算精度的影響

為了討論網(wǎng)格尺寸大小對(duì)計(jì)算精度的影響，以5 kg TNT裝藥，在水面下1 m處爆炸為例，在理論氣泡脈動(dòng)第一周期最大半徑的1.5倍范圍內(nèi)分別劃分0.9、1.0、1.1、1.2倍裝藥半徑(0.09 m)的均勻網(wǎng)格，其余水域網(wǎng)格劃分由內(nèi)至外逐漸增大，最大處為內(nèi)部均勻網(wǎng)格尺寸的10倍。模型大小為20 m×20 m×30 m，水域高15 m，空氣域高15 m。

將計(jì)算所得的氣泡體積等效為規(guī)則圓的半徑并與經(jīng)驗(yàn)式(1)結(jié)果進(jìn)行對(duì)比，得到了不同網(wǎng)格條件下的氣泡第一周期脈動(dòng)最大半徑(見表6)，可以看出在幾種工況下，不同的網(wǎng)格尺寸的相對(duì)誤差均不超過(guò)5%。其中，當(dāng)網(wǎng)格尺寸與裝藥半徑為0.9倍和1.0倍時(shí)，氣泡第一周期脈動(dòng)最大半徑分別為2.629 m和2.633 m，兩者與經(jīng)驗(yàn)公式的理論半徑都較為接近，但采用0.9倍裝藥半徑與網(wǎng)格尺寸的計(jì)算時(shí)長(zhǎng)相較于1.0倍裝藥半徑與網(wǎng)格尺寸約增加一半。在綜合考慮計(jì)算精度和計(jì)算效率的條件下，本文計(jì)算模型中心區(qū)域的均勻網(wǎng)格尺寸選取為1.0倍裝藥半徑大小。

表6 不同網(wǎng)格條件下氣泡參數(shù)計(jì)算結(jié)果

2 氣泡運(yùn)動(dòng)對(duì)水面形態(tài)的影響

2.1 模型設(shè)置

由于驗(yàn)證模型采用的是5.2 g的球形TNT炸藥，不便對(duì)氣泡運(yùn)動(dòng)對(duì)水面形態(tài)變化的影響進(jìn)行研究，因此建立尺寸為10 m×10 m×30 m的模型，水域、空氣域各高15 m，采用5 kg球形TNT裝藥。

2.2 裝藥深度對(duì)氣泡運(yùn)動(dòng)特性的影響

水下爆炸氣泡在不同裝藥水深條件下將會(huì)出現(xiàn)不同的運(yùn)動(dòng)特性，為了研究裝藥深度對(duì)其運(yùn)動(dòng)特性的影響，首先定義裝藥水深H為水下裝藥位置距水面的距離，炸藥起爆位置在水面以下1～5 m的模型中心處，設(shè)置9種不同工況，各個(gè)工況之間間隔0.5 m。

由不同裝藥深度條件下氣泡垂直方向位移的變化情況(見圖4)可以看出，在H為1.0～3.0 m時(shí)，氣泡垂直方向位移呈現(xiàn)先上后下的變化趨勢(shì)，且移動(dòng)大小隨裝藥水深的增大而減小，當(dāng)H>3.0 m時(shí)，氣泡垂直方向位移逐漸向上移動(dòng)。

注：a～i分別代表H為1.0、1.5、2.0、2.5、3.0、3.5、4.0、4.5、5.0 m。圖4 水深對(duì)氣泡垂直方向位移的影響Fig.4 Influence of water depth on the vertical displacement of bubble

由不同裝藥深度下氣泡半徑的變化情況(見圖5)可以發(fā)現(xiàn)，氣泡脈動(dòng)第一周期時(shí)長(zhǎng)隨著裝藥深度的增大逐漸延長(zhǎng)，從以上分析可知在2.0～3.0 m裝藥水深范圍內(nèi)存在一個(gè)臨界水深。當(dāng)裝藥水深小于臨界水深時(shí)，氣泡脈動(dòng)第一周期最大半徑與脈動(dòng)第一周期時(shí)長(zhǎng)隨裝藥水深的增大而增大，氣泡中心點(diǎn)的位移先向上后向下，且移動(dòng)大小隨裝藥水深增大而減??；當(dāng)裝藥水深大于臨界水深時(shí)，氣泡脈動(dòng)第一周期最大半徑與脈動(dòng)第一周期時(shí)長(zhǎng)隨水深的增大而減小，氣泡中心點(diǎn)向上位移大小隨水深增大而增大。

注：a～i分別代表H為1.0、1.5、2.0、2.5、3.0、3.5、4.0、4.5、5.0 m。圖5 不同裝藥深度氣泡半徑變化關(guān)系Fig.5 Change relationship of bubble radius with different charge depth

2.3 水冢形成過(guò)程

氣泡運(yùn)動(dòng)的變化也會(huì)導(dǎo)致水面水冢形態(tài)發(fā)生變化，因此對(duì)水冢的形成過(guò)程進(jìn)行研究非常有必要。為了研究水冢的形成過(guò)程，對(duì)近水面條件下5 kg球形TNT炸藥在改變不同水深距離參數(shù)時(shí)產(chǎn)生的不同類型的水冢形成過(guò)程進(jìn)行研究。同時(shí)為了與水面中心水柱相區(qū)分，定義在水面附近形成的水體形態(tài)稱之為水冢。經(jīng)前人研究發(fā)現(xiàn)，水?，F(xiàn)象的形成與水深距離參數(shù)γf=H/Rmax有關(guān)[5]，H為炸藥裝藥位置距水面的距離，m；Rmax為氣泡脈動(dòng)過(guò)程的最大半徑，m。依據(jù)水下爆炸氣泡最大半徑公式，在水深距離參數(shù)γf范圍在0.000 1～1.1內(nèi)，討論不同水深距離參數(shù)γf下水冢形態(tài)的特征及形成機(jī)理(空腔半徑為氣泡脈動(dòng)第一周期所形成的實(shí)際最大半徑)。

當(dāng)起爆位置距水面相近時(shí)，會(huì)產(chǎn)生零碎型水冢(見圖6)，此時(shí)水深距離參數(shù)γf為0.000 195，裝藥水深0.000 4 m，氣泡膨脹形成的水下空腔半徑為2.053 5 m。炸藥起爆后，形成的氣泡突出水面與空氣相連通而破碎(見圖6a、圖6b)；在t=354 ms氣泡膨脹至最大半徑，形成了水下空腔，水面水體形成了零碎型水冢；在重力作用下空腔四周水體向內(nèi)匯聚并形成匯聚水柱(見圖6d、圖6e)；t=2 094 ms時(shí)刻匯聚水柱高度到達(dá)最大值；隨后受重力作用沖擊水面形成波浪傳播。

圖6 零碎型水冢發(fā)展過(guò)程Fig.6 Development process of fragmented water mound

當(dāng)水深距離參數(shù)γf增大到0.082 9，此時(shí)裝藥水深下降到0.215 m，氣泡膨脹形成的水下空腔半徑為2.589 m時(shí)會(huì)形成飛濺型水冢(見圖7)。炸藥起爆后，在t=6、54 ms時(shí)刻氣泡突出水面發(fā)生破碎，并在t=174 ms時(shí)刻達(dá)到最大半徑，形成空腔；t=294、594 ms時(shí)刻，由于氣泡的帶動(dòng)，使得水體在空腔上部聚集，產(chǎn)生背向射流，形成飛濺型水冢。接著便在t=774、1 194 ms時(shí)刻產(chǎn)生類矩形的水面水柱，并在t=1 614 ms時(shí)刻達(dá)到最大值。

圖7 飛濺型水冢發(fā)展過(guò)程Fig.7 Development process of splashing water mound

接著當(dāng)水深距離參數(shù)γf增加到0.549時(shí)，裝藥水深為1.5 m，氣泡脈動(dòng)第一周期實(shí)際最大半徑為2.732 m會(huì)形成酒杯型水冢(見圖8)。在t=6 ms時(shí)刻，氣泡開始膨脹，并在t=118 ms時(shí)刻達(dá)到最大半徑；t=258 ms時(shí)刻，氣泡在收縮過(guò)程中，由于受到自由液面Bjerknes力的排斥作用，會(huì)形成向下的氣泡射流，進(jìn)而帶動(dòng)水流形成向上的射流，并逐漸發(fā)展為水面中心水柱；t=318 ms時(shí)刻，氣泡在向下射流的帶動(dòng)下，分裂為兩部分，并在t=338 ms時(shí)刻達(dá)到第一周期的最小半徑；t=478 ms時(shí)刻氣泡再次膨脹，并且中部和上部氣泡相融合，同時(shí)在氣泡的推動(dòng)和中心水柱的阻礙下，向上運(yùn)動(dòng)的水體在中心水柱周圍形成了環(huán)狀水體，環(huán)狀水體與中心水柱一起構(gòu)成了酒杯型水冢。

圖8 酒杯型水冢發(fā)展過(guò)程Fig.8 Development process of wine glass water mound

最后當(dāng)水深距離參數(shù)γf增大到0.998時(shí)，裝藥水深為2.692 m，此時(shí)氣泡脈動(dòng)第一周期實(shí)際最大半徑為2.697 m，形成的水冢形態(tài)為皇冠型水冢(見圖9)。在t=6 ms時(shí)刻，氣泡開始膨脹，并在t=118 ms時(shí)刻達(dá)到最大半徑；t=258 ms時(shí)刻，氣泡在收縮過(guò)程中，由于受到自由液面Bjerknes力的排斥作用，會(huì)形成向下的氣泡射流，進(jìn)而帶動(dòng)水流形成向上的射流，并逐漸發(fā)展為水面中心水柱；t=318 ms時(shí)刻氣泡收縮至最小半徑，并在t=338 ms時(shí)刻開始二次脈動(dòng)，同時(shí)向上的水射流也在不斷發(fā)展，并促進(jìn)了中心水柱的發(fā)展，環(huán)狀水體也在氣泡膨脹的推動(dòng)下開始突出水面向上發(fā)展，并與中心水柱一起構(gòu)成了皇冠型水冢。

圖9 皇冠型水冢發(fā)展過(guò)程Fig.9 Development process of crown-shaped water mound

2.4 水冢氣泡運(yùn)動(dòng)特性

對(duì)上述4種不同類型的水冢氣泡運(yùn)動(dòng)特性進(jìn)行分析，由于零碎型水冢和飛濺型水冢工況下的水下爆炸氣泡，炸藥裝藥水深較淺，不能形成穩(wěn)定的氣泡脈動(dòng)周期和射流，因此不對(duì)其運(yùn)動(dòng)特性進(jìn)行討論。

對(duì)于酒杯型水冢和皇冠型水冢工況下的氣泡，裝藥水深能夠滿足氣泡脈動(dòng)所需的條件，可以形成穩(wěn)定的氣泡脈動(dòng)與射流。通過(guò)兩者氣泡脈動(dòng)半徑、周期的對(duì)比(見圖10a)可以看出，酒杯型水冢氣泡脈動(dòng)第一周期的最大半徑為2.732 m，周期為336 ms；皇冠型水冢氣泡脈動(dòng)第一周期的最大半徑為2.697 m，周期為360 ms。這表明氣泡最大半徑隨著水深的增大而減小，氣泡脈動(dòng)周期則相反，這一現(xiàn)象與文獻(xiàn)[20]得到的結(jié)論相同。

從氣泡中心點(diǎn)在垂直方向的位移變化(見圖10b)可以看出，在膨脹階段，酒杯型水冢的氣泡中心點(diǎn)在垂直方向上的位移明顯高于皇冠型水冢，這是由于酒杯型水冢氣泡所處位置水深較淺，氣泡受到的外部靜水壓強(qiáng)也較小，因此膨脹階段的阻力也是偏小的。在收縮階段，由于受到Bjerknes力的作用形成向下射流，因此在射流和氣泡收縮趨勢(shì)的帶動(dòng)下，氣泡中心位置不斷下移。

圖10 酒杯型水冢、皇冠型水冢氣泡運(yùn)動(dòng)特征對(duì)比Fig.10 Comparison of bubble movement characteristics between wine glass and crown-shaped water mound

2.5 不同類型水冢的特征參數(shù)對(duì)比

對(duì)比上述4種不同水冢的特征尺寸(見圖11)可知，幾種類型的水冢寬度基本處于4～5 m，與氣泡最大直徑相近，呈現(xiàn)出一定的相似性，結(jié)合上文對(duì)水冢發(fā)展過(guò)程進(jìn)行詳細(xì)的研究，可以知道水冢的形成過(guò)程與氣泡的運(yùn)動(dòng)息息相關(guān)，氣泡的運(yùn)動(dòng)特性決定了水冢形態(tài)。

圖11 不同類型水冢特征尺寸對(duì)比Fig.11 Comparison of the characteristic sizes of different types of water mounds

零碎型和飛濺型水冢的最大寬度是氣泡在破碎前形成的空腔直徑。飛濺型水冢寬度大于零碎型水冢，這是由于零碎型水冢的起爆水深更淺，氣泡會(huì)更快速的膨脹到空氣中發(fā)生破碎。

酒杯型水冢和皇冠型水冢的水冢寬度，主要與氣泡進(jìn)行第二周期脈動(dòng)推動(dòng)水面形成的射流位置相關(guān)，皇冠型水冢中心水柱的向上速度更低，形成的中心水柱底部寬度更大，使得氣泡第二周期脈動(dòng)形成的環(huán)狀水體在更靠外側(cè)的位置形成，水冢寬度也因此更大。

對(duì)于零碎型水冢和飛濺型水冢的中心水柱高度，形成的主要原因是空腔四周水體在重力作用下向空腔底部匯聚形成的匯聚水柱。由圖6c及圖7c可知，飛濺型水?？涨恢車w范圍更大，且空腔高度更高，從能量角度分析，飛濺型水冢形成的匯聚水柱應(yīng)該會(huì)比零碎型水冢匯聚水柱更高，但由于飛濺型水冢在空腔頂部形成的交匯向下射流與匯聚水流的運(yùn)動(dòng)方向相反，阻礙了匯聚水柱的形成，使得飛濺型水冢形成的中心水柱高度更低，這一現(xiàn)象也體現(xiàn)在飛濺型水冢的匯聚水柱的上下斷面寬度基本相同，而零碎型水冢的匯聚水柱為錐形。

對(duì)比酒杯型水冢和皇冠型水冢的水柱高度，由圖11可知，酒杯型水冢形成的中心水柱高度為35.213 m，皇冠型水冢形成的中心水柱高度僅為4.437 m，兩者相差較大。這是由于皇冠型水冢所處位置水深較大，向下氣泡射流的劇烈程度相較于酒杯型水冢是偏低的，因此形成的向上射流速度也較低。在向上射流形成的時(shí)刻，垂直方向距氣泡頂部1 m位置處，皇冠型水冢向上射流速度為9.27 m/s，而酒杯型水冢的向上射流速度為29.8 m/s，同時(shí)皇冠型水冢氣泡上部水體厚度更大，向上射流要穿透水體進(jìn)入空氣中形成中心水柱所要消耗的能量也更大。綜上原因，在皇冠型水冢與酒杯型水冢的工況下雖然都能形成穩(wěn)定的氣泡脈動(dòng)、射流，但形成的中心水柱高度存在極大的差距。

2.6 不同類型水冢的水面波動(dòng)特性

水下爆炸引起的水面波動(dòng)在一定程度上會(huì)對(duì)水上結(jié)構(gòu)物施加波浪壓力，對(duì)于類似核彈等超大裝藥量情況下的水下爆炸，形成的水面波動(dòng)甚至可以等效為海嘯波。根據(jù)前文所形成的的水冢，分析不同類型水冢形成的水面波動(dòng)特性。

距炸藥氣泡中心5 m處不同形態(tài)水冢形成的水面波動(dòng)如圖12所示。對(duì)于零碎型水冢的水面波動(dòng)，水冢形態(tài)形成的800 ms以內(nèi)，由于氣泡破碎在100 ms左右已經(jīng)破碎，自由液面只是在空腔膨脹過(guò)程中，水面出現(xiàn)了較低程度的上升。在1 000～2 200 ms左右的液面下降是空腔匯聚水柱開始形成，到2 650 ms匯聚水柱達(dá)到最高時(shí)刻之間，中心水柱不斷升高，四周水體向中心區(qū)域流動(dòng)，在一定的時(shí)間延遲之后帶動(dòng)5 m處自由液面下降。自2 650 ms水柱達(dá)到最高點(diǎn)之后，在重力作用下匯聚水柱開始沖擊水面形成水面波動(dòng)，并持續(xù)向外傳播，2 200～4 000 ms水面上升，此時(shí)水面最大的上升高度為31.570 cm。

圖12 不同形態(tài)水冢形成的水面波動(dòng)時(shí)程變化Fig.12 Time-history of water surface fluctuation formed by different forms of water mounds

對(duì)于飛濺型水冢的水面波動(dòng)，在水冢形成階段，相較于零碎型水冢，飛濺型水冢對(duì)水面的影響更加明顯，這主要是由于其形成的空腔范圍更大，氣泡的膨脹帶動(dòng)水體向上運(yùn)動(dòng)但并未破碎飛濺，形成了連續(xù)的水面隆起，帶動(dòng)了監(jiān)測(cè)點(diǎn)處的自由液面位置上升。在400～1 900 ms之間，監(jiān)測(cè)點(diǎn)處自由液面下降，一方面是由于前一階段氣泡脈動(dòng)破碎之后水面的隆起消失，更主要的是與零碎型水冢相同，空腔四周水體開始匯聚形成匯聚水柱，帶動(dòng)周圍水體向中心匯集，使水面降低。在2 800 ms時(shí)刻匯聚水柱達(dá)到最大高度，此后開始在重力作用下沖擊水面，形成波浪向外傳播。監(jiān)測(cè)點(diǎn)處的水面最大上升高度為64.060 cm，大致為零碎型水冢的2倍，雖然飛濺型水柱高度比零碎型水柱高度大0.325 m，兩者的體積分別為38.780 m3和15.450 m3，但是飛濺型水冢所具有的重力勢(shì)能更大，因此沖擊水面形成的興波波浪高度也更大。

對(duì)于酒杯型水冢的水面波動(dòng)，在水冢形成初期，氣泡迅速膨脹頂起的水體在氣泡收縮階段的下落形成了初期的水面波動(dòng)。在400 ms之后，隨著中心水柱和四周環(huán)形水體的上升，四周水體向中心流動(dòng)，監(jiān)測(cè)點(diǎn)的水面略微下降。隨后四周的環(huán)形水體和中心水柱根部附近的大量水體在重力作用下持續(xù)沖擊水面，形成了較大的水面波動(dòng)，監(jiān)測(cè)點(diǎn)最高的水面波動(dòng)位移為83.085 cm。在最高水面波動(dòng)通過(guò)監(jiān)測(cè)點(diǎn)后，可以看到后續(xù)出現(xiàn)一個(gè)小的水面波動(dòng)，以及一個(gè)持續(xù)的水面上升過(guò)程，則是殘余的部分環(huán)狀水體和部分下落的中心水柱沖擊水面形成的殘余水面波動(dòng)。

對(duì)于皇冠型水冢的水面波動(dòng)，與酒杯型水冢類似，在水冢形成初期，氣泡膨脹頂起的水面下墜形成了一定程度的水面波動(dòng)。伴隨著水冢的形成，四周水流向中心流動(dòng)，使得監(jiān)測(cè)點(diǎn)處自由液面出現(xiàn)一定程度的下降。隨后在中心水柱和四周環(huán)狀水體形成的水冢下墜沖擊下，自由液面形成波動(dòng)?；使谛退Ｔ诒O(jiān)測(cè)點(diǎn)處形成的最大水面波動(dòng)位移為53.640 cm，相對(duì)酒杯型水冢較低，形成該現(xiàn)象的原因一方面是由于中心水柱的高度更低，僅為4.337 m。另一方面是由于氣泡多次脈動(dòng)形成的自由液面破碎情況更加劇烈，后續(xù)的環(huán)狀水體與之前形成的環(huán)狀水體相互沖擊消耗了能量。從圖3中可以觀察到在最高水面波動(dòng)通過(guò)監(jiān)測(cè)點(diǎn)之后，還形成了兩個(gè)小的水面波動(dòng)，前者是部分殘余的環(huán)狀水體沖擊水面形成的，后者是后續(xù)中心水柱上部水體沖擊水面形成的。

3 結(jié)論

1)隨著水深距離參數(shù)的增大，水下爆炸氣泡逐漸由不具備完整周期的破碎狀態(tài)變?yōu)榫哂忻}動(dòng)射流的完整周期狀態(tài)，同時(shí)由于Bjerknes力減小，氣泡向上運(yùn)動(dòng)趨勢(shì)逐漸增大。

2)隨水深距離參數(shù)的增大，水冢形態(tài)依次為零碎型、飛濺型、酒杯型、皇冠型；前兩者形成水冢形態(tài)與氣泡破碎時(shí)帶動(dòng)的水體運(yùn)動(dòng)有關(guān)，水冢的特征尺寸受氣泡半徑大小和水冢形成機(jī)理影響；后兩者形成的水冢形態(tài)是氣泡脈動(dòng)和射流現(xiàn)象疊加造成的，水冢的特征尺寸與氣泡大小及形成的射流速度相關(guān)。

3) 隨著水深距離參數(shù)的增大，零碎型水冢和飛濺型水冢的興波機(jī)理為：在重力作用下，氣泡破碎后空腔四周水體在空腔底部匯聚形成的匯聚水柱沖擊自由液面形成水面波動(dòng)；酒杯型水冢和皇冠型水冢在重力作用下，由組成水冢的中心水柱和周圍環(huán)狀水體，沖擊自由液面形成水面波動(dòng)。