水下爆炸作用下艦船總縱強度模型實驗方案研究

張海鵬，岳永威，蘇羅青，張阿漫，姚熊亮

(哈爾濱工程大學 船舶工程學院，哈爾濱 150001)

提高艦船的抗爆抗沖擊性能和生命力是保證艦船發揮其作戰能力［1］的一項重要前提，各國海軍十分重視該方面的研究。對艦船進行實船爆炸實驗是評估艦船總縱強度最準確有效的方式，但需要耗費大量的資金，且實施較為困難。因此，許多發展中國家在評估艦船水下爆炸作用下的總縱強度［2-3］方面，基本上采用艦船縮比模型開展水下爆炸實驗來評估，此時會產生這樣的矛盾:若縮尺比較大，則實驗費用居高不下，開展縮比模型實驗的意義大為降低;若縮尺比較小，縮比模型的尺寸會非常小，模型在制造工藝上存在較大困難，且換算到實船上結果的準確性也得不到保證，我國雖然開展了一些縮比模型的水下爆炸實驗，但其縮尺比均較小，模型與實船既未能達到完全幾何相似，也不存在畸變相似關系，因此其實驗結果難以評估與考核實船的抗爆抗沖能力。

在評估艦船總縱強度水下爆炸實驗中，若采用合適的縮尺比，再以此艦船縮比模型為基準模型，根據彈性結構畸變準則，選擇合適的畸變參數對基準模型中幾何尺寸較小的板梁進行畸變，使用畸變模型進行實驗。這樣，不僅能解決模型在制造工藝中的困難，還可以在降低實驗成本的前提下實現準確預測艦船原型總縱強度的目的。本文的加載方式考慮了水下爆炸沖擊波和氣泡載荷的一次相似律，并利用彈性結構畸變準則對基準模型進行畸變，通過數值計算手段得到滿足工程精度、制造工藝、經濟性等要求的艦船總縱強度模型實驗方案。

1 縮比模型水下爆炸相似規律

多年來，國內外許多學者一直致力于水下爆炸“相似律”的研究和探索，尋求通過縮比模型實驗預測原型結構動響應特性的方法，現今也取得了許多具有學術價值［4］的成果。水下爆炸是一個很復雜的物理過程，分為沖擊波和氣泡脈動兩個階段，流體的特性也會在不同的階段發生變化，如果簡單的將模型加載方式用同一個沖擊波、氣泡脈動相似準則處理，得到的模型實驗數據結論將與原型發生偏差，甚至背離。因此需要對水下爆炸產生的沖擊波和氣泡的相似律分別討論。另外在工程實際中，某些時候縮比模型的相似條件不能滿足，如模型中的板或梁的幾何尺寸按縮尺比進行縮尺，縮比后的模型在制作中會產生諸多問題，如焊接難度大、所需規格市場上不存在，因此就必須利用彈性結構的完全相似理論及畸變準則對基準模型的結構進行改造得到畸變模型，最終實現對原型總縱強度的準確預測。

1.1 沖擊波載荷相似律

自由場中水下爆炸沖擊波現象區別于一般水動力學現象，它的作用力是瞬態的并且具有很強的沖量，水動力學中一些相似常數(如Fr數)可以不必討論。本文僅考慮自由場中的水下爆炸問題，要使自由場水下爆炸中沖擊波階段的過程相似，必須使各基本相似量間滿足一定的相似條件，即:

若用裝藥量W表示，即r0=，相似條件可利用經驗公式［7］由已知的方程導出:

同理可求得:λE=λi=λθ=λl，即當模型的幾何尺寸、加載形式與原型相似時，兩者的壓力峰值是相同的，而脈沖持續時間、沖量、比沖量、能量密度將按特征長度相似常數(λl)換算。

1.2 氣泡脈動載荷相似律

對于自由場中水下爆炸氣泡脈動問題，本文假定氣泡膨脹時周圍的水產生徑向運動，而沒有垂向位移，即不考慮與氣泡上浮時流體靜壓變化有關的流體靜浮力的影響［8］。在這個假定的條件下，要使自由場水下爆炸中氣泡脈動階段的過程相似，必須使各基本相似量間滿足一定的相似條件，即:

式中:Pd為爆深D(m)處靜水壓(Pa)，等于D+D0(m)水柱壓強，D0為大氣壓的水柱壓強;為氣泡脈動壓力峰值;T為一次氣泡脈動周期;θ為衰減系數;ρ0為裝藥密度;r0為裝藥半徑;Ep為單位質量氣泡能;為氣泡最大半徑;ρw為水密度。選定Pd、ρw和r0為基本物理量，同樣由π定理得到各物理量函數關系式為:

以一次氣泡脈動周期項T項作為因變項，并用裝藥質量及水柱表示，經轉換得到Ep、和各物理量與周期的函數關系式分別用G3、G4和G5表示［6］:

由經驗公式［6］:

類比后求得關系式［6］:

綜上所述，在不考慮重力效應影響的情況下，自由場中氣泡脈動相似的的充分必要條件是:(D+D0)是個常量。要實現這一充要條件，必須使用將實驗水池做封閉增壓處理，在(D+D0)一定時，當模型的幾何尺寸、加載形式與原型相似時，在一次氣泡脈動階段，氣泡最大半徑、脈動周期以及脈動頻率按特征長度相似常數(λl)換算，氣泡脈動峰值不變。事實上，以目前國內的實驗條件，通過對水池進行封閉增壓來實現水下爆炸一次氣泡脈動的相似難以實現，大部分模型試驗主要在露天水池中進行，此時(D+D0)就不可能是一個常量，對氣泡脈動的模擬就存在誤差，然而，要得到在不考慮對實驗水池進行封閉增壓處理時，大縮尺比船模水下爆炸實驗進行原型總縱強度預報的誤差范圍，還需要嚴格的理論推導及大量實驗，本文主要介紹水下爆炸實驗模型的選取方法，這里不做詳細討論。

1.3 艦船模型畸變相似準則

船體板格通常為矩形板，按邊界不同大致可分為剛性固定、彈性固定、彈性支持和簡支。下面以均布荷重P(t)作用下四邊剛性固定的矩形板為例，討論其板厚不能按特征長度相似常數λ縮小時的畸變修正方法。此時將均布荷重P(t)作用下四邊剛性固定的矩形板作為一個系統，取表征系統這一現象的物理量:均布荷重P(t)、時間t、板厚d，t時刻的位移U(t)、t時刻的應變ε(t)、t時刻的應力σ(t)、材料泊松比μ、結構固有頻率f、系統總能量Q、材料密度ρ、t時刻的結構尺寸L(t)和彈性模量E，共12個物理量。當縮比模型與原型使用同種材料時有:

當板厚不能按特征長度相似常數λL(t)縮小時，定義完全幾何相似模型為基準模型(上標“m”代表基準模型的∏項和物理量)，以基準模型為基礎僅對板厚進行畸變的模型定義為畸變模型(上標“b”代表畸變模型的∏項和物理量)。對于基準模型和畸變模型，僅有相似常數λd≠λL(t)，其它各相似常數均為1時，縮比模型發生畸變，定義“板厚畸變系數”為 φ，根據文獻［9-11］有:

設畸變模型和基準模型時間的因變∏項的比值為δt，稱之為“時間預測系數”。即:

時間預測系數δt僅是畸變系數φ的函數，通過大量實驗測量或數值計算的方法可得到δt隨φ變化的規律，進而確定其函數形式和最終結構畸變準則。相關文獻表明［11］，彈性體結構基于該方法確定的畸變準則進行畸變，得到的模型計算結果與基準模型有很高的一致性，并且在制作工藝上的可控性好。

2 計算結果分析

為選定艦船總縱強度模型實驗的最優方案，本文對四種艦船縮尺實驗方案進行水下爆炸數值計算分析，并對實際模型橫剖面彎矩預測偏差進行總結，得到既能使實驗模型在工程精度范圍內實現對原型橫剖面彎矩的準確預測，同時又滿足制造工藝的可控性和實驗成本的經濟性的最終實驗方案。

2.1 模型實驗計算工況設定

本文以XXX艦有限元模型為原型，如圖1所示，選取攻擊武器為MDM-X型水雷，爆心位于船舯正下方，見圖2，裝藥質量相當于m=Wkg的TNT，調整爆距R=Z0m，使其氣泡周期［7］與船體濕模態一階固有頻率相同，使氣泡脈動載荷激起船體共振，探究XXX艦此典型工況下的總縱強度。

為了評估縮比模型與原模型的數值計算結果，本文選取了延船長方向上三處典型剖面部位進行橫剖面彎矩時歷曲線輸出，分別為船舯橫剖面、X1橫剖面、X2橫剖面，具體示意圖如圖3，其中X1=X2。

2.2 完全幾何相似模型實驗方案

為了選擇合適的縮尺比，本研究提出四種完全幾何縮比模型實驗方案與原型進行對比分析，采用自由場中水下爆炸一次沖擊波和氣泡脈動相似理論對模型進行相似加載，并考慮實驗模型的經濟性及可操作性等因素，確定合適的縮比模型，為下一步畸變模型實驗方案做準備。對上述艦船完全幾何縮小，艦船、藥包原型與縮比模型相似關系見表1，其中Dp為爆深D(m)處的水柱壓強，(D0)p=10.33 m為大氣壓的水柱壓強。

表1 XXX型艦原型與模型相似關系Tab.1 Similarity relation between the model and the archetype of XXX ship

下面分別給出XXX艦及其縮比模型在靜水中遭受水下爆炸載荷作用時各個考核剖面的無量綱橫剖面彎矩時歷曲線。此處的原型與縮比模型彎矩時歷曲線由相應剖面的輸出點擬合而得，將各模型時歷曲線的輸出點橫坐標t按時間縮尺比放大，同時將各輸出點縱坐標M按相應的彎矩縮尺比放大，即得到相似處理后的模擬時歷曲線。無量綱橫剖面彎矩時歷曲線的原型與縮比模型相比，縮比模型的橫剖面彎矩縮尺比按特征長度相似常數三次方即換算，縮比模型的時間縮尺比按特征長度相似常數λl換算，相似處理后的時歷曲線結果如圖4～圖9(數據已經過無量綱處理):

從上圖可以看出，應用自由場中水下爆炸一次沖擊波和氣泡脈動相似理論處理后的縮比模型與原型的典型考核橫剖面彎矩時歷曲線基本吻合。應用完全幾何相似模型，可以通過模型實驗非常準確地預測原型結構遭受水下爆炸載荷時的結構響應特性及其總縱強度。

2.3 畸變模型實驗方案

完全幾何相似模型中的板或梁的幾何尺寸(如厚度)按縮尺比進行縮尺時，會出現如焊接難度大、鋼材規格難以實現等問題，必須使用畸變模型。模型畸變后，圖3所考慮的考核剖面抗彎剛度必定發生變化，保證畸變模型預測精度的關鍵在于合理選取畸變系數φ，由前文推導可知，預測系數δt是畸變系數φ的函數，但二者的函數形式并不確定，而且不同模型的畸變系數函數形式也不盡相同，實際上，預測系數δi與φ的函數關系可以通過一系列模型試驗來確定，但這樣做需要耗費大量的人力物力，顯然工程上是不可取的。因此，本文采用有限元預測系數法尋找預測系數δi隨畸變系數φ的變化規律。具體方案如下:首先定義完全幾何相似模型為基準模型，以基準模型為基礎，僅對板厚進行畸變的模型定義為畸變模型，建立系列畸變模型和基準模型之間的預測系數δi=F(φ)函數關系，結合完全幾何相似理論，實現對原型動響應特性的預報。

數值實驗結果表明，當畸變系數φ=3(板厚為3h的模型Model 4)時，畸變模型對原型的預測誤差最小。下面分別給出當畸變系數φ=3時，XXX艦的基準模型及各畸變模型的在靜水中遭受水下爆炸載荷作用無量綱橫剖面彎矩時歷曲線，及經過相似處理的無量綱橫剖面彎矩時歷曲線對比圖。

選取以上橫剖面的動響應峰值進行預測偏差分析，見表2:

表2 典型部位彎矩峰值預測偏差分析Tab.2 Deviation analysis of moment peak of typical position

表3 畸變實驗模型評估表Tab.3 Evaluation table of distortion experimental model

3 結論

本文基于自由場中水下爆炸一次沖擊波和氣泡脈動相似律，通過數值實驗手段研究了在大縮尺比船模水下爆炸實驗中，使用畸變模型實現對原型總縱強度準確預測的方法。通過對實驗結果的分析可得到的結論及建議如下:

(2)應用完全幾何相似模型，可以通過模型實驗準確地預測原型結構遭受水下爆炸氣泡脈動載荷時的結構響應特性及其總縱強度。

［1］汪 玉.艦船現代沖擊理論及應用［M］.北京:科學出版社，2005:21-29.

［2］徐培譽.水下爆震對結構之響應［D］.臺北:國立臺灣大學，2004.

［3］ Keil K A.The response of ship to underwater explosion［J］.SNAME，1961，69(4):366-410.

［4］程素秋，寧永成，張 臣，等.相似理論在水下爆炸模型實驗中的應用［J］.船舶科學技術，2008，30(3):96-99.

［5］張 亮，李云波.流體力學［M］.哈爾濱:哈爾濱工程大學出版社，2007:179-187.

［6］張效慈.水下爆炸實驗相似準則［J］.船舶力學，2007，2(1):108-118.

［7］庫爾P，羅耀杰，等譯.水下爆炸［M］.北京:國防工業出版社，1960.

［8］Klaseboer E，Khoo B C，Hung K C.Dynamics of an oscillating bubble near a floating structure［J］.Journal of Fluid and Structure，2005，21:395-412.

［9］ Gel’fand B E，et al.Similarity criteria for underwater explosion［J］.Combustion，Explosion，and Shock Waves，2004，40(2):214-218.

［10］徐 挺，著.相似方法及其應用［M］.北京:機械工程出版社，1995.

［11］謝建林.瞬態載荷作用下彈塑性結構動響應相似性研究［D］.哈爾濱:哈爾濱工程大學，2008:24-55.