艦船抗爆領域水下爆炸載荷研究進展

牟金磊 朱 錫

海軍工程大學 船舶與動力學院，湖北 武漢 430033

艦船抗爆領域水下爆炸載荷研究進展

牟金磊 朱 錫

海軍工程大學 船舶與動力學院，湖北 武漢 430033

海戰中水下爆炸載荷是艦船的重要威脅之一，為了給艦船抗爆研究提供準確的載荷輸入，需要對載荷進行系統研究。通過查閱大量文獻，介紹了國內外艦船抗爆領域水下爆炸載荷的研究概況，從沖擊波載荷和氣泡載荷兩個方面總結了該領域的研究進展。由于沖擊波載荷的研究比較成熟，重點對氣泡運動方程、氣泡射流等氣泡階段的研究工作進行了總結和分析。并且在分析前人工作的基礎上對該領域有待進一步解決的問題和發展趨勢進行了展望，為該領域內研究工作的開展提供一定參考。

水下爆炸；沖擊波；氣泡運動方程；氣泡射流

1 引 言

在海戰中，水雷、魚雷和深水炸彈等水中兵器在水下爆炸產生的沖擊波和氣泡脈動壓力等載荷，能夠導致艦船局部或總體結構的嚴重毀傷，使艦船喪失戰斗力。水下爆炸載荷是水面艦船和潛艇的重要威脅。為了提高艦船抗爆能力，各國海軍都非常重視水下爆炸載荷的研究。特別是各海軍強國，如美國、澳大利亞等，在裝藥的設計和性能、爆炸機理、艦船響應等方面，采用理論分析、數值模擬、模型試驗和實船試驗等手段，基本涵蓋了水下爆炸及其對艦船毀傷作用研究的主要方面，并且取得了大量成果［1－9］，大大提高了艦船的生命力。我國海軍也逐漸重視該領域的研究工作，國內眾多研究單位和學者都不同程度地投入到該領域的研究之中。

炸藥在水中爆炸，爆轟波傳到藥包表面，壓縮水介質形成具有陡峭波頭的水中沖擊波。高度壓縮的爆轟產物以氣泡的形式向外膨脹，推動周圍的水徑向流動。氣泡內的壓力隨著氣泡膨脹擴大而不斷下降，當壓力降至周圍水介質的靜壓時，由于水的慣性運動，氣泡過度膨脹，直至達到最大半徑。此時氣泡內的壓力低于周圍水的平衡壓力，周圍的水開始作反向運動，向氣泡中心聚合，使氣泡不斷地收縮，造成氣泡內部壓力不斷增加。同樣，由于聚合水流的慣性運動，氣泡被過度壓縮直至達到最小半徑。這種氣泡脈動次數可達十幾次以上，但在有限水介質中氣泡脈動的次數則要少得多。圖1顯示了水下爆炸時水中某點壓力P的時域曲線、氣泡脈動半徑R的時域曲線及與其對應的氣泡膨脹和壓縮過程。水下爆炸沖擊波載荷，峰值壓力大，但是作用時間短，可以使艦船局部變形或破裂。氣泡脈動的壓力比沖擊波小得多，然而其作用時間遠超過首次沖擊波，沖量與沖擊波相當，甚至超過沖擊波，且氣泡運動的周期與艦船整體結構頻率相當，容易使船體產生“鞭狀響應”，對艦船的整體強度有很大影響。當炸藥在離艦船較近處爆炸時，氣泡受艦船結構邊界的影響容易失穩而產生射流，在沖擊波破壞的基礎上進一步引起艦船結構的毀傷。

對水下爆炸載荷的研究始于19世紀初。第一、二次世界大戰期間，由于戰爭的需要，對于水下爆炸載荷的研究進展很快，形成了系統的理論。庫爾在其所著的《水下爆炸》一書中對爆炸沖擊波和氣泡脈動等現象進行了全面描述，揭示了各個爆炸參數之間的關系［10］，該書至今仍被世界各國公認為經典之作。另外一部水下爆炸的權威著作是前蘇聯科學家Zamyshlyayev等著的《Dynamic Loads in Underwater Explosion》，該書在庫爾研究成果的基礎上作了進一步發展，全面系統地論述了沖擊波和隨后的壓力波，采用理論和試驗相結合的方法，重點討論了自由面和底部效應、繞射效應、空化效應、沖擊波和結構之間的相互作用，并且給出了包含氣泡脈動過程的水中壓力波的解析表達式，更便于應用。

2 沖擊波載荷研究

關于沖擊波在水中的傳播，庫爾在 《水下爆炸》一書中介紹了三種沖擊波理論［10］。分別是基爾克烏特—別澤理論、基爾克烏特—布林克里理論和賓尼理論。他們都是根據流體動力學的基本方程式得到的近似解。但是三種理論的解法和所作假設的物理根據卻不同。當距離與藥包半徑之比大于10時，按基爾克烏特－布林克里理論計算出來的峰值壓力與距離的關系曲線和實驗值很一致；按基爾克烏特－別澤理論所得到的曲線則比實驗值大15%～20%；按賓尼理論所得到的曲線則比實驗值要低，但是在藥包附近則急劇增加。應該說賓尼方法較適合計算沖擊波最初的生成和傳播，基爾克烏特－別澤理論則更適合計算較遠距離上的沖擊波。相對于以上沖擊波理論而言，在工程上人們更習慣于使用由相似理論得到的經驗公式，而不同學者得到的經驗公式卻各不相同［11］。庫爾給出的指數衰減的沖擊波壓力計算經驗公式如式（1）所示，與基爾克烏特—布林克里理論計算結果和試驗結果都吻合較好，一直被各國艦船抗爆研究者引用，是世界各國公認的計算沖擊波壓力的經典公式：

式中，P（t）為時刻測點壓力；Pm為沖擊波峰值壓力，計算公式見式（2）；θ為沖擊波指數衰減的時間常數。

式中，W為炸藥重量；R為測點爆距；常數K、α為與炸藥有關的常數，通常通過試驗測定。

關于水下爆炸沖擊波的理論研究基本上都基于庫爾的《水下爆炸》一書中建立的理論框架。此后，學者們開展了一些試驗和仿真研究，在相似律、邊界對沖擊波傳播的影響等方面做了大量工作。SNAY H G詳細闡述了將爆炸相似律運用于水下爆炸時的基本前提和使用方法［12］。TEMKIN S A在小藥量試驗中發現自由表面附近水下爆炸沖擊波超壓峰值出現了非線性滯后，分析發現這是由于自由表面發射引起的［13］。NADAMITSU Y等通過理論、試驗和仿真相結合的方法對水下爆炸沖擊波的Von Neumann反射現象進行了研究，提出了與試驗結果吻合很好的理論公式［14］。SAYAPIN A等比較了幾種測量水下沖擊波壓力的方法，分析了各種測量方法的優缺點［15］。我國很多研究機構在水下爆炸試驗中歸納了沖擊波物理參數的試驗數據，并歸納了廣泛應用的計算沖擊波壓力的經驗公式系數［16］。錢勝國通過分析近自由水面水下爆炸情況下自由水面對沖擊波傳播過程中的影響和沖擊波反射后的變化規律，提出從爆深與爆炸能量溢出的關系來修正庫爾沖擊波壓力公式的觀點。張鵬翔、顧文彬等分析了淺水中水面和水底反射對沖擊波參數的影響，指出水面對沖擊波切斷影響較強，水底反射的稀疏波對沖擊波波形切斷影響較弱［17］。周方毅等對無限水介質中爆炸的沖擊波壓力計算公式在實際應用中應如何選用進行了辨析，提出了不同條件下推薦使用的公式［18］。

國內外對水下爆炸沖擊波的研究開展較早，主要依靠試驗手段研究，取得了大量研究成果。關于水下爆炸的沖擊波理論計算方法和試驗經驗公式不斷改進，近年來的主要工作集中于對沖擊波經驗公式的修正以及對適用條件的限制，為以后研究工作的開展奠定基礎，使沖擊波階段的載荷計算更加精確。但是沖擊波載荷方面仍有許多尚未解決的問題。比如藥包近場沖擊波的傳播規律尚不十分清楚，主要是由于近場的超高壓難以測量，隨著試驗技術的進步，該方面的研究也將進一步深入。另外隨著近年來非理想高能炸藥在水中兵器中的應用，對于非理想高能炸藥的沖擊波研究擺在了研究者的面前。關于非理想炸藥的沖擊波峰值、時間衰減規律等計算是否滿足相似律等一系列問題都有待進一步解決。

3 氣泡載荷研究

在20世紀80年代以前，水下爆炸研究主要集中于沖擊波方面，對氣泡載荷認識比較粗淺。但是，從80年代中期起，研究者意識到氣泡對結構的損傷可能比沖擊波更嚴重，目前水下爆炸的氣泡階段已經成為國際上研究的熱點。

爆炸沖擊波過后，爆炸產物形成的氣泡中還有相當部分能量，氣泡在與周圍水介質的相互作用下，膨脹和壓縮，產生滯后流與脈動壓力。水下爆炸產生的氣泡脈動壓力具有極強的破壞力。氣泡的脈動壓力可以引起船體的鞭狀響應［19］，導致艦船整體縱向總體屈曲破壞和大變形。當氣泡靠近船體時，氣泡收縮還可能引起射流，產生非常高的局部壓力，造成船體的局部嚴重破壞。

3.1 氣泡運動方程研究

對氣泡脈動過程的理論研究以對氣泡運動方程的研究為主線。1923年，Lamb建立了第一個氣泡膨脹運動方程，即著名的Rayleigh－Plesset方程，該方程是在假設流體不可壓縮基礎上建立的［20］。1941年，HERRING考慮了流體的可壓縮性，建立了氣泡膨脹運動方程［21］。但是他們的研究都沒有考慮氣泡脈動過程中的上升運動。1986年，HEATON等人建立了對稱軸獨立變化的橢圓體模型，研究了氣泡運動過程中的非球形效應和能量輻射損失的影響。研究結果發現，由于非球形效應使得氣泡在脈動過程中的實際上升速度低于按球形計算的理論值［22］。Temkin綜述了小藥量水中爆炸所產生的氣泡脈動壓力在水中的傳播規律，他在分析了非線性聲學的適用性之后指出在很多情形下可以忽略非線性因素。

VERNON T A［23］假設流體為不可壓、無旋理想流體，應用勢流理論，采用點源模擬爆炸氣泡，根據能量守恒原理建立了爆炸氣泡動力學運動方程，分別考慮了深水中、有浮力和自由面時的氣泡脈動問題，其建立的方程為了考慮能量的耗散，采用了形狀阻力系數，并以此為基礎引入了虛擬力的概念。VERNON模型能夠計算氣泡脈動過程中速度、壓力等的變化過程，其反映的規律與實際過程接近，尤其是第一次氣泡脈動過程，但是Vernon模型存在的主要問題是，將脈動氣泡假定為絕熱的，與外界無熱交換作用，事實上這是不可能的。這是Vernon方程不能正確預報第一次氣泡脈動以后行為的根本原因。

水下爆炸現象傳統上被視為兩個不相干的過程：沖擊波和氣泡過程，而氣泡過程初始條件的選取成為難題。 HUNTER K S（2001）［24，25］和 GEERS T L（2002）［26，27］將水下爆炸問題視為由沖擊波階段和氣泡脈動階段組成的單一過程，第一階段為第二階段提供初始條件。HUNTER K S在其博士論文中［24］，通過遠場壓力曲線的相似關系采用一個簡單的體積加速度模型來確定沖擊波階段氣泡的運動。在沖擊波階段采用體積加速度模型（Volume－Acceleration Model）以便確定隨后的氣泡脈動階段氣泡運動初始條件。氣泡脈動階段的運動將采用雙重漸進近似模型（DAA，Doubly Asymptotic Approximation）以便計算氣泡半徑a和氣泡遷移速度u。Geers－Hunter模型借鑒了水下爆炸沖擊波問題研究中的一些方法，較好地處理了水下爆炸氣泡脈動過程中的能量損耗問題，克服了不可壓流體假設下氣泡脈動模型的缺陷。并在國內外的相關研究中得到了采用，具有較好的應用前景，已被大型商業有限元程序ABAQUS所采用，用于在水下爆炸相關研究中計算氣泡脈動問題。模型中有5個控制方程：

式中，

其中，a為氣泡半徑；u為氣泡上浮位移；ρg為內部氣體密度；cg為內部氣體中聲速；ρl為外部流體密度； cl為外部流體中聲速；ζ＝ρgcg／ρlcl為氣泡內部氣體的聲阻抗和氣泡外部流體的聲阻抗比值；pI為氣泡所處位置的靜水壓力，包括大氣壓；φl0為氣泡外表面水的膨脹速度勢；φl1為平動速度勢；φg1為氣泡內表面氣體速度勢。

炸藥爆炸的絕熱狀態方程為：

式中，Kc為炸藥材料的絕熱常數；Vc為炸藥體積；V為當前的氣泡體積；γ為氣體的比熱。

該模型的不足之處主要是沒有引入邊界面對氣泡脈動的影響，且仍要采用一些借助于試驗確定的經驗參數。

國內船舶科學研究中心的張效慈等人引入余能率的概念。根據實測，各次脈動余能率在0.46～0.66之間，脈動余能呈幾何級數遞減。根據這一方法對Vernon模型進行改進，能得到與實際情況更為接近的結果［28］。梁龍河等也采用Vernon模型編制了拉氏一維不定常彈塑性流體動力學計算程序——SINL程序，對水下爆炸產生的諸如氣泡脈動規律、脈動周期、水中沖擊波壓力的變化等特性進行了研究，給出了不同裝藥水下爆炸產生的氣泡脈動的一些規律、脈動周期變化規律、氣泡和水交界面處的壓力曲線等結果［29］。 相關文獻［29，30］介紹了在Vernon模型的基礎上進行的進一步改進，考慮了水平剛性面和垂直剛性面邊界條件下的氣泡脈動運動，并建立了相應的求解方程，擴展了求解范圍。

目前，通過求解氣泡運動方程來求解載荷的方法，還僅限于簡單情況下采用，通常需要很多理想假設。在實際的工程背景下，涉及到流固耦合、氣泡動力學等學科知識，難以通過該方式求解。比如炸藥在船底中部附近近距爆炸時，若氣泡與船體結構接觸，則由于船體中部吃水被爆炸氣泡排開，艦船底部產生負壓，船體在重力和負壓作用下產生中垂，可能會折斷破壞。若氣泡與結構沒有接觸，在氣泡收縮過程中，隨著流體運動，在艦船底部會形成一個持續時間較長、范圍較大的負壓，可能會使艦船整體出現嚴重的中垂破壞。目前很少有人研究由于水下爆炸引起的艦船底部的負壓場對艦船結構的影響。特別是利用氣泡運動方程來分析負壓場的形成條件、大小、變化規律等特性的研究至今尚屬空白。

3.2 氣泡射流載荷研究

對于氣泡載荷的研究，研究人員最初只關注到了球對稱氣泡的運動。但在后來的試驗和仿真研究中發現，由于邊界的影響，氣泡會出現非球形情況：氣泡在膨脹階段被結構表面輕微地排斥開，而在坍塌階段被結構表面強烈地吸引，這時在氣泡內部將會形成一股射流，這股射流產生于遠離結構表面的一側，并且高速穿過氣泡，直到它撞擊到氣泡壁的另一側。高速射流可能會引起艦船結構的局部毀傷。

在理論研究方面，BJERKNES C A和BJERKNES V F K發現當流體中脈動的流體同相振蕩時彼此吸引，而當它們反相振蕩時又相互排斥。最初，他們試圖利用這一發現來解釋電磁機理［32］，后來這種在流體中脈動物體之間的相互作用力就被稱為Bjerknes力。這為后人研究氣泡與邊界的相互作用提供了一個非常重要的力學概念。BENJAMIN和ELLIS［33］將開爾文沖量作為一個分析工具，解釋了邊界面處氣泡的運動，為該領域的后續研究奠定了基礎。 他們以 Rm（ρΔp）1／2為特征量將開爾文沖量無量綱化，其中：Rm為氣泡的最大半徑；ρ為流體密度；ΔP＝P∞－Pc。 設 P∞為炸藥深度處無限域流場的壓力、Pc為飽和蒸汽壓（通常為常數）則開爾文沖量的表達式為：

式中，ρ為流體密度；φ為速度勢；Sb為氣泡表面；n為流體的外法線。如果認為包含氣泡的無限流體域內線性動量守恒，將外部邊界延伸到無窮大，可以得到開爾文沖量的下述關系：

式中，V為氣泡的體積；Ss＋Sf表示剛性邊界和自由面邊界；ex為浮力方向的單位矢量。因此，開爾文沖量體現了浮力和由剛性壁面或自由面引起的Bjerknes力對氣泡的作用，它決定了氣泡的運動方向。

BLAKE 等人［34，35］在前人工作的基礎上，基 于開爾文沖量推導了粗略估算邊界附近氣泡的運動方向的Blake準則。通過引入無量綱參數δ＝（ρgRm／ΔP）1／2表示浮力的作用。 從物理的角度看，它對應著最大半徑為Rm的氣泡在重力作用下運動一個半徑所花的時間與氣泡半個周期的比值。定義無量綱參數γ＝d／Rm，表示Bjerknes力的重要性。其中g為重力加速度，d為氣泡中心到邊界的距離，ΔP＝P∞－Pc，P∞為炸藥深度處無限域流場的壓力，Pc為飽和蒸汽壓（通常為常數）?？梢约俣汩_爾文沖量線將γ-δ參數域分成幾個區域，氣泡被吸引向邊界的區域以及氣泡被排斥的區域。應用該方法，可以發現分界線由γδ＝0.442給出，稱為Blake準則。Blake利用該準則判斷固壁附近氣泡運動方向，如圖 2 所示［36］。

Blake準則推導過程中，由于采用了許多假設，因此結論具有局限性。相對于理論研究，試驗研究可以更直觀的顯示氣泡運動和坍塌過程。早期的試驗主要依靠電火花、激光等產生。最初，研究人員用放電的方法來模擬水下爆炸氣泡的動態特性，典型的例子是CHAHINE等人采用水下放電的方法模擬水下爆炸氣泡的動態特性獲得了成功［37］。這種方法產生的氣泡具有安全、高效、費用低廉等優點，但是難以控制氣泡的大小、形狀以及氣泡內部的能量。利用激光脈沖產生氣泡，可以很好地控制氣泡的形狀、位置等參數。LAUTERBORN等對激光產生的氣泡進行研究，得到了氣泡坍塌產生射流的速度與氣泡中心到壁面距離的關系［38］。VOGEL等通過試驗測量了激光產生的氣泡坍塌時產生射流的周圍流場分布，清楚地顯示了氣泡膨脹、收縮、坍塌發生射流全過程的細節［39，40］。 隨著電子技術、光測技術的迅猛發展，水下爆炸試驗測試技術得到了很大的提高，特別是高速攝影技術的快速發展，使得捕捉氣泡動力學行為成為可能。近年來，直接利用藥包爆炸產生氣泡的試驗研究取得了較大的進展，BOYCE P和DEBONO S［41］在戶外水池中進行水下爆炸試驗，對剛壁附近氣泡的運動狀態進行高速攝影拍攝，捕捉了水下爆炸氣泡射流形成的過程，為水下爆炸氣泡的研究提供了寶貴的實驗數據和圖像。KLASEBOERK［42］用試驗方法研究水下爆炸氣泡與結構之間的相互作用，在水池中安裝高速攝影儀捕捉由水下爆炸產生的氣泡的膨脹、坍塌及回彈的過程。在試驗過程中，使用了幾種不同類型的結構，包括剛性和彈性的圓形板、平板，平板的變形用非接觸的激光傳感裝置測得，這些試驗研究為水下爆炸氣泡的研究提供了珍貴的資料。國內，張寒虹等利用高速陰影攝影和掃描攝影觀察了水中電爆炸的基本物理現象［43，44］。 張立等在淺水體外將CCD高速攝影機實拍和圖像處理技術引入氣泡脈動參數的測試研究，得到了0.3 g DDNP裝藥的氣泡半徑隨時間變化的運動軌跡和脈動周期，與電測結果和經驗公式計算結果都非常一致［45］。朱錫等利用高速攝影技術進行了一系列不同藥量的TNT炸藥水下爆炸試驗，對氣泡脈動周期和最大氣泡半徑等參數進行了研究［46］。不論是研究范圍還是研究深度，國內與國外先進水平都還存在很大差距。由于很多研究涉及軍事秘密，國外有價值的資料都很難獲得，這就要求國內研究者要立足于自主研究，追趕國際先進水平。

理論研究僅限于較規則的模型，試驗研究可重復性比較差，尤其是水下爆炸試驗花費巨大，受場地、試驗條件等眾多因素制約，且偶然性因素特別多。隨著計算機技術的不斷發展，數值仿真研究得到了迅速發展。國外從上世紀70年代就開始了對水下爆炸氣泡的數值仿真研究，從早期的邊界元法（BEM）［34］和有限體積法（VOF）［47－49］，發展到了近年興起的光滑粒子水動力學法 （SPH）［50－52］。WANG Q X等對自由面附近的氣泡的非線性變化進行了數值仿真，利用邊界積分法求解速度勢函數的Laplace方程，得到了氣泡坍塌形成射流的三種模式［53］。PEARSON A利用邊界積分法計算自由面下一個和下兩個氣泡的高度非線性運動，該理論可以推廣到有限個氣泡，并分析了氣泡射流問題及其引起的水面效應［54］。KALUMUCK K M利用邊界元法研究了氣泡與臨近結構之間的相互影響，提出由于結構的響應，引發了氣泡周期改變以及射流形成的結構表面的壓力特性的改變［55］。國內雖然起步較晚，但是該領域的研究越來越受到國內學者的重視，且發展迅速，取得了很多成果。戚定滿，魯傳敬等利用邊界元法對軸對稱氣泡進行研究，模擬了氣泡膨脹及潰滅的過程，對軸對稱氣泡在不同表面附近潰滅產生射流進行了較好地模擬［56－58］。方斌等利用商用軟件對不同邊界條件下的氣泡進行了仿真［59］，分析了水下爆炸氣泡的半徑、外形變化、氣泡脈動壓力、流場中的速度分布等問題。牟金磊等重點對近自由面的氣泡射流現象進行了仿真研究［60］，主要分析了射流方向與射流速度的問題。宗智等分別利用SPH法和邊界積分法，對水下爆炸二維和三維氣泡進行了仿真研究，獲得了氣泡運動規律、流場壓強變化、結構危險載荷及其分布等情況［61，62］。姚熊亮、張阿漫等以邊界元積分法為基礎，對水下爆炸氣泡進行了系統的數值仿真研究，對自由面、剛性壁面、水中結構及各種復雜邊界附近的氣泡運動特性進行了詳細研究，對影響氣泡射流的因素進行分析，得到了很多有價值的成果［63－71］。

氣泡射流的相關研究文獻很多，國內外學者通過理論、試驗和仿真各種方法對氣泡射流現象進行研究，也取得很多研究成果，對氣泡射流現象的認識也不斷深入。但是由于氣泡射流現象的復雜，涉及到多相流問題、流體動力學問題等，仍有很多方面沒有研究透徹。而且，多數文獻的研究僅僅關注氣泡本身的動態特性，對氣泡塌陷產生射流的條件、現象等進行分析，而對射流產生的載荷研究不夠充分。另外氣泡射流能否應用相似律，射流產生后對周圍的流體壓力場有何影響等問題都還有待進一步研究。其中，作為艦船結構毀傷研究輸入條件的射流載荷的相關研究無疑是非常重要的。

4 研究建議及展望

水下爆炸載荷研究是一項非常復雜的課題，其中涉及流固耦合、氣液耦合、氣液固三相耦合等難題，該領域的研究還有大量工作要開展。結合本文查閱的文獻，提出幾點研究建議。

1）近場沖擊波傳播規律研究。由于試驗條件限制，目前該類試驗開展尚不多見，隨著試驗技術的不斷提高，近場沖擊波的特性將會成為下一階段研究的重點方向。

2）利用氣泡運動方程求解艦船結構附近氣泡載荷。從理論上求解作用在艦船結構上的水下爆炸載荷，為艦船結構抗爆抗沖擊研究提供準確的載荷輸入，涉及到多個前沿學科，難度比較大，目前尚無很理想的解決方案。

3）氣泡射流產生的載荷研究。氣泡塌陷發生射流會給艦船結構造成嚴重的局部毀傷，射流形成的條件、機理以及產生的載荷大小的定量計算等研究還很不完善，尚有大量工作待開展。

隨著國內外學者對水下爆炸載荷及其作用下艦船典型的結構毀傷研究不斷深入，對于指導人們改進艦船設計水平，提高其抗爆抗沖擊能力具有重要意義，同時也能為提高水中兵器作戰能力提供參考。

［1］STETTLER J W.Damping mechanisms and their effects on the whipping response of a submerged submarine subjected to an underwater explosion［R］.ADA298743，1995.

［2］MENON S.Experimental and Numerical Studies of Underwater Explosions［R］.ADA317378， 1996.

［3］SCHNEIDER N A.Prediction of surface ship response to severe underwater explosions using a virtual underwater shock environment［R］.ADA417555， 2003.

［4］KRUEGER S R.Simulation of cylinder implosion initiated by underwater explosion［R］.ADA457530， 2006.

［5］UCAR H.Dynamic response of a catamaran－hull ship subjected to underwater explosions［R］.ADA462708， 2006.

［6］STEPHEN D B.Acceleration of a plate subject to explosive blast loading－tiral results DSTO－TN－0270［R］.Australia：Defence Science and Technology Organisation，2000.

［7］JOHN P B.The Effect of Non－spherical Collapse on Determination of Explosion Bubble Parameters，DSTO－RR－0238 ［R］.Australia：Defence Science and Technology Organisation，2002.

［8］MARIAN F，BOX P，MCLEAN A.Dynamic Impulse Calibration of 100MPa Blast Pressure Transducers ［R］.Australia：Defence Science and Technology Organisation，DSTO－TN－0525，2003.

［9］JOHN M B，MICHAEL B，TERRY T，et al.An Experimental Facility for Imaging of Medium Scale Underwater Explosions［R］.Australia：Defence Science and Technology Organisation，DSTO－TR－1432，2003.

［10］庫爾.水下爆炸［M］.北京：國防工業出版社，1960.

［11］HENRYCH C J.The dynamic of explosion and its use［M］.New York.Elsevier Scientific Publishing Company，1979.

［12］SNAY H G.The Scaling of Underwater Explosion Phenomena［R］.AD0271468， 1962.

［13］TEMKIN S A.Review of the Propagation of Pressure Pulse Produced by Small Underwater Explosion Charges［R］.NOLTR－6181，1988.

［14］NADAMITSU Y，LIU Z Y，FUJITA M，et al.Von Neumann reflection of underwater shock wave ［J］.Journal of Material Processing Technology，1999，85（1/3）：48－51.

［15］SAYAPIN A，GRINENKO A，EFIMOV S，et al.Comparison of different methods of measurement of pressure of underwater shock waves generated by electrical discharge［J］.Shock Waves， 2006，15（2）：73－80.

［16］惲壽榕，趙衡陽.爆炸力學 ［M］.北京：國防工業出版社，2005.

［17］張鵬翔，顧文彬，葉序雙.淺層水中爆炸沖擊波切斷現象淺探［J］.爆炸與沖擊，2002，22（3）：222－228.

［18］周方毅，陳曉強，張可玉，等.無限水介質中爆炸沖擊波壓力計算公式辨析［J］.爆破，2003，20（1）：7－9.

［19］李玉節，張效慈，吳有生.水下爆炸氣泡激起的船體鞭狀運動［J］.中國造船，2001，42（3）：1－7.

［20］PLESSET M S.The dynamics of cavitation bubbles ［J］.Journal of Applied Mechanics.1949，16：277－282.

［21］HERRING C.Theory of the Pulsations of the Gas Bubble Produced by an Underwater Explosions［R］.USA：National Defence Research Council Report，C4－sr20－010，1941.

［22］HEATON K C.Effects of Non－Sphericity and Radiative Energy Loss on the Migration of the Gas Bubble from Underwater Explosion［R］.ADA166823， 1986.

［23］VERNON T A.Whipping response of ship hulls from underwater explosion bubble loading［R］.ADA178096，Defence Research Establishment Atlantic DARTMOUTH，1986.

［24］KENDALL S H.Global－Shape－Function Models of an Underwater Explosion Bubble［D］.USA，Boulder：University of Colorado，2001.

［25］HUNTER K S，GEERS T L.Pressure and velocity fields produced by an underwater explosion ［J］.Journal of the Acoustical Society of America，2004，115（4）：1483－1496.

［26］GEER T L，HUNTER K S.An integrated wave－effects model for an underwater explosion bubble ［J］.Journal of the Acoustical Society of America，2002，111 （4）：1584－1601.

［27］GEER T L，Chung－Kyu Park.Optimization of the G＆H bubble model［J］.Shock and Vibration，2005，12（1）：3－8.

［28］張效慈，李玉節，趙本立.深水爆炸水動壓力場對潛艇結構的動態影響［J］.中國造船，1997，139（4）：61－68.

［29］梁龍河，曹菊珍，袁仙春.水下爆炸特性的一維球對稱數值研究［J］.高壓物理學報，2004，18（3）：203－208.

［30］方斌，朱錫，陳細弟，等.水平剛性面下方水下爆炸氣泡垂向運動的理論研究 ［J］.爆炸與沖擊，2006，26（4）：345－350.

［31］方斌，朱錫，張振華.垂直剛性面邊界條件下水下爆炸氣泡運動的理論研究［J］.海軍工程大學學報，2007，19（2）：81－85.

［32］Bjerknes.Fields of Force［M］.Columbia University Press，1966：45－47.

［33］BENJAMIN T B， ELLIS A T.The collapse of cavitation bubbles and the pressures thereby produced against solid boundaries ［J］.Philosophical Transactions of the Royal Society of London ，1966，260：221－240.

［34］BLAKE J R，GIBSON D C.Growth and collapse of a vapour cavity near a free surface ［J］.Journal of Fluid Mechanics， 1981， 111： 123－140.

［35］BLAKE J R，GIBSON D C.Cavitation bubbles near boundaries［J］.Annual Review of Fluid Mechanics，1987，19：99－123.

［36］BLAKE J R，CERONE P.A note on the impulse due to a－vapour bubble near a boundary.Journal of Math Science， 1982， 23： 383－93.

［37］CHAHINE G L，FREDERICK G S，LAMBRECHT C J，et al.Spark－Generated Bubbles Laboratory－Scale Models Of Underwater Explosionsuse For Validation Of Simulation Tools ［C］//SAVIAC Proceedings 66th Shock and Vib.Symposium，Biloxi，MS，1995，Vol.2，265－276.

［38］LAUTERBORN W，BOLLE H.Experimental investigations of cavitation－bubble in the neighbourhood of a solid boundary［J］.J.Fluid Mech.1975，72：391－399.

［39］VOGEL A， LAUTERBORN W.Time－resolved particle image velocity used in the investigation of cavitation bubble dynamics.Applied Optics，1988，27（9）：1869－1876.

［40］VOGEL A，LAUTERBORN W.Acoustic transient generation by laser－produced cavitation bubbles near solid boundaries.J Acoust Socam，1988，84（2）：719－731.

［41］BOYCE P，DEBONO S.Report of underwater explosion tests［R］.Centre Technique Des Systeme Navals，France，2003.

［42］KLASEBOER E，HUNG K C，WANG C，et al.Experimental and numerical investigation of the dynamics of an underwater explosion bubble near a resilient／rigid structure［J］.J.Fluid Mech.2005，537： 387－413.

［43］張寒虹，陳志福.水中電爆炸的高速攝影研究［J］.光子學報，2001，30（7）：889－891.

［44］張寒虹，陳志福，張弛.水中電爆炸的實驗研究［J］.爆炸與沖擊，2002，22（4）：363－367.

［45］張立，章橋龍，郭進，等.裝藥淺水爆炸的氣泡脈動參數研究［J］.煤礦爆破，2005（2）：5－8.

［46］朱錫，牟金磊，洪江波，等.水下爆炸氣泡脈動特性的試驗研究 ［J］.哈爾濱工程大學學報，2007，28 （4）：365－368.

［47］ROGERS J C W，SZYMCZAK W G，BERGER A E，et al.Numerical solution of hydrodynamic free boundary problems［J］.International Journal of Numerical Mathematics，1990，95：241－266.

［48］POPINET S，ZALESKI S.Bubble collapse near a solid boundary：A numerical study of the influence of viscosity［J］.Journal of Fluid Mechanics，2002，464：137－163.

［49］GLIMM J，GROVE J W，LI X L，et al.Three–dimensional front tracking［J］.J Sci Comput，1998，19：703～727.

［50］LIU M B，LIU G R，LAM K Y，et al.Mesh free particle simulation of the detonation process for high explosives in shaped charge unlined cavity configurations ［J］.Shock Waves，2003，12（6）：509－520.

［51］LIU M B，LIU G R，LAM K Y，et al.Smoothed particle hydrodynamics for numerical simulation of underwater explosion［J］.Computational Mechanics，2003，30（2）：106－118.

［52］LIU G R，ZONG Z，LAM K Y.Computer simulation of high explosive explosion using smoothed particle hydrodynamics methodology ［J］.Computers and Fluids，2003，32： 85－98

［53］WANG Q X，YEO K S，KHOO B C，et al.Nonlinear interaction between gas bubble and free surface ［J］.Computers and Fluids，1996，25（7）：607－628.

［54］PEARSON A，COX E，BLAKE J R，et al.Bubble interactions near a free surface ［J］.Engineering Analysis with Boundary Elements， 2004，（28）：295－313.

［55］KALUMUCK K M，DURAISWAMI R，CHAHINE G L.Bubble dynamics fluid－structure interaction simulation by coupling fluid BEM and structural FEM codes［J］.Journal of Fluid and Structures.1995，9：861－883.

［56］戚定滿， 魯傳敬.單氣泡演化及輻射噪聲 ［J］.上海交通大學學報，1998，32（12）： 50－54.

［57］戚定滿，魯傳敬，何友聲.兩氣泡運動特性研究［J］.力學季刊，2000，21（1）：16－20.

［58］LU C J，HE Y S，ZHU S Q.Transient cavity collapse in the vicinity of a flexible boundary ［J］.Hydrodynamics，1996，1305－1310.

［59］方斌，朱錫.不同邊界條件下水下爆炸氣泡的數值模擬［J］.海軍工程大學學報，2008，20（2）：85－90，112.

［60］牟金磊，朱錫，張振華.近自由面水下爆炸氣泡現象的數值仿真研究.艦船科學技術，2008，30（4）：113－116.

［61］宗智，鄒麗，劉謀斌，等.模擬二維水下爆炸問題的光滑粒子（SPH）方法［J］.水動力學研究與進展 A 輯，2007，22（1）：61－67.

［62］宗智，何亮，張恩國.水中結構物附近三維爆炸氣泡的數值模擬［J］.水動力學研究與進展 A 輯，2007，22（5）：592－602.

［63］姚熊亮，張阿漫.簡單Green函數法模擬三維水下爆炸氣泡運動［J］.力學學報，2006，38（6）：749－759.

［64］張阿漫，姚熊亮.近壁面氣泡的運動規律研究［J］.物理學報，2008，57（3）：1662－1671.

［65］張阿漫，姚熊亮.近壁面氣泡的動態特性研究［J］.中國造船，2008，49（1）：1－10.

［66］張阿漫，姚熊亮.近自由面水下爆炸氣泡運動規律研究［J］.物理學報，2008，57（1）：339－353.

［67］張阿漫，姚熊亮.單個三維氣泡的動力學特性研究［J］.應用力學學報，2008，25（1）：107－112.

［68］張阿漫，姚熊亮.復雜邊界附近氣泡的動態特性研究［J］.力學季刊，2008，29（1）：24－32.

［69］張阿漫，姚熊亮，李佳，等.不同邊界附近氣泡的三維數值模擬與實驗值對比研究［J］.中國科學G輯：物理學力學天文學，2008，38（7）：896－907.

［70］張阿漫，姚熊亮.近邊界三維水下爆炸氣泡動態特性研究［J］.爆炸與沖擊，2008，28（2）：124－130.

［71］張阿漫，姚熊亮.水深和藥量的變化對水下爆炸氣泡射流的影響研究［J］.工程力學，2008，25（3）：222－229.

Research Progress on the Underwater Explosive Load for Warship Shock Resistance

Mu Jin－lei Zhu Xi
College of Naval Architecture and Power， Naval University of Engineering，Wuhan 430033，China

The underwater explosion is one of the major threats to the warships in the sea battle field.In order to provide accurate loads input for the study of warship explosion protection，a systematic study on the loads is desired.By reviewing numerous literature both at home and abroad，researches on under water explosive load in the field of warship anti－explosion were summarized， especially the research advances on the loads by the shock wave and the air bubble.The research on the loads by shock wave is relatively mature， therefore， much more efforts must be made on the investigation and analysis on the air bubble motion equation and fluid jets from the bubble.Based on the previous work in the area， some of the problems to be solved are raised and advancements are predicted，which can provid reference for the study in this area.

underwater explosion； shock wave； motion equation of bubble； bubble fluid jets

O383.1；U661.6

A

1673－3185（2010）02－01－08

2009－09-17

國家安全重大基礎研究（國防973）項目（51335020103）

牟金磊（1980－），男，博士研究生。研究方向：艦船結構水下抗爆抗沖擊。E-mail：mjl1007＠163.com

朱 錫（1961－），男，教授，博士生導師。 研究方向：艦船結構抗爆。E-mail：zhuxi816＠ 163.com