水下舷外發射水密隔膜運動特性數值仿真研究

吳明晨，蘇 杭，趙 錚

(1.中船重工713研究所六室， 河南 鄭州 450052； 2.南京理工大學 能源與動力工程學院， 江蘇 南京 210094)

水下發射裝置的研究及其設計，對于提高魚雷、潛射導彈等武器的出水品質具有十分重要的意義。對于有關彈丸的水下發射過程的研究，李亞男[1]通過建立熱力學模型，對小型運載器的彈丸發射過程進行了仿真；李志華[2]通過建立魚雷的運動方程，作出了魚雷的彈道曲線，并確定魚雷不會產生很大的袋深現象；程廣濤[3]對水下發射系統進行了數學建模，確定了發射過程中發射系統參數和武器運動參數。

水下發射的主要困難在于整流罩和耐壓殼之間的空間狹小，對發射裝置水密隔膜的結構設計要求較高，水密隔膜在張開或散布過程中不能與鄰筒發生干涉。目前水下發射裝置水密隔膜多采用不系留方式，即水密隔膜破裂后與發射筒分離。該方式雖然結構簡單，但無法控制水密隔膜破片的散布范圍，在多筒同時發射或小間隔發射時隔膜破片容易碰撞鄰筒彈體。本文設計了一種水下弦外發射水密隔膜系留結構，通過轉軸連接水密隔膜和發射筒，使破裂后的水密隔膜系留在發射筒上，控制隔膜翻轉角度避免與鄰筒發生干涉。

由于發射筒為多單元聯裝結構，單元間距離較小，水密隔膜在張開過程中所占據的空間逐漸增大，因此，需要計算在不同發射條件下，水密隔膜的最大張開角度。本文主要研究水密隔膜在發射過程中的運動特點，不考慮發射裝置的動力部分及相關影響。

1 計算模型的建立

本文所研究的水下舷外發射裝置由發射筒、水密隔膜、滑塊滑軌、動力系統和彈體組成?；瑝K位于彈體尾部，動力系統推動滑塊使彈體沿滑軌在發射筒內加速運動，彈頭沖破水密隔膜后彈體發射出筒。計算模型由水域、魚雷、魚雷發射筒、水密隔膜和相對應安裝在發射筒上的鉸鏈轉軸所組成。

計算域后部為魚雷，魚雷前部為球冠，后部為空心的圓柱體，材料為鋼，模型質量與實際魚雷質量相同；魚雷前部為四瓣水密隔膜，材料為有機玻璃，扁球殼狀；水域為給定水深的海水；魚雷發射筒和水密隔膜將魚雷與水域分隔開；自行設計鉸鏈轉軸機構，在Solid Edge軟件環境下進行建模，其結構模型圖如圖1，計算域模型如圖2。

圖1 鉸鏈轉軸模型三視圖及三維圖

圖2 計算域示意圖

2 有限元分析模型的建立

流固耦合計算需要同時考慮流場沖擊和彈丸運動的相互作用，包括流場計算、固體力學計算、耦合計算等。在保證計算精度的條件下，需要對計算模型簡化，以減少計算時間，提高計算速度。

采用cm-g-μs單位制對計算域進行建模，本系統的有限元模型由水域、魚雷發射筒、水密隔膜、彈丸和鉸鏈轉軸組成。經過ICEM軟件進行網格劃分后，得到鉸鏈轉軸的有限元網格如圖3，整個計算域的有限元網格如圖4所示。

圖3 鉸鏈轉軸網格三視圖及三維圖

圖4 計算域網格

其中水介質采用歐拉網格建模， GRUNEISEN狀態方程[4]；彈丸、發射筒和鉸鏈轉軸均為鋼材料，采用拉格朗日網格建模，MAT_JOHNSON_COOK材料模型[5]；水密隔膜采用有機玻璃材料，采用MAT_ELASTIC彈性材料模型[6]。計算中采用的水介質和發射裝置的材料及狀態方程參數如表1、表2所示，未給出的參數取 LS-DYNA軟件和工程實踐中的默認值。

表1 不同模型的材料參數

表2 水介質的材料參數

約束設置與接觸設置為：在水-結構邊界，通過Euler/Lagrange罰函數耦合算法建立約束方程，將發射裝置結構與水域流體進行耦合計算[7]；采用轉軸實現水密隔膜與發射筒的連接，鉸鏈轉軸的固連面采用面面固連接觸，接觸面采用自動面面接觸。為節約求解時間和計算量，將發射裝置放在流場的中心對稱位置。對于外流場設置，通過定義不同的水深-密度曲線間接模擬不同水深時的環境壓力[8]；在水域邊界面施加透射邊界約束，模擬無限水域。

對于初始條件設置，設定彈丸以給定速度曲線運動；設定發射裝置以恒定速度運動，即以一定的艇速運動。

對于控制設置，為了避免沙漏現象的產生，采用沙漏控制語句*CONTROL_HOURGLASS，在語句中的默認沙漏系數為0.03；*CONTROL_TIMESTEP控制語句的時間步長比例因子為0.9，有效的防止了負體積的產生；基本計算時間為0.6 s，每隔0.6 ms輸出一個d3plot結果文件。

3 水密隔膜破裂過程分析

通過基于 LS-DYNA程序的多物質ALE 算法[9]，對魚雷發射過程進行了三維有限元模擬，并采用 LS-PREPOST 軟件對 LS-DYNA 求解后得到的d3plot文件進行后處理[10]。在仿真過程中，為了減少計算量，外流場只選取了離彈丸較近的部分。

圖5為在發射過程中，發射裝置在不同時刻下(t=0 s，0.10 s，0.15 s，0.20 s)的狀態變化圖。由圖5可以看出，隨著時間的進行，水密隔膜繞著鉸鏈轉軸轉動，從而水密隔膜由封閉狀態緩慢張開，使彈丸有空間出筒。在發射過程中水密隔膜與彈丸沒有直接接觸。

4 不同發射條件下的水密隔膜轉角變化

4.1 不同彈速條件下的影響

計算了不同彈速條件、水深150 m、無艇速條件下，水密隔膜的轉動角度隨時間的曲線，如圖6。

由圖6可以看出，水密隔膜轉角先在t=0.2 s左右快速增大，后保持一定角度不變，當彈丸出筒后，由于彈丸尾部產生的低壓區造成較大的壓差阻力[11]，因此在發射裝置外部有海水補入，導致發射裝置產生振動[12]，水密隔膜轉角達到最大值，并且在狀態穩定以后轉角保持不變。在不同的彈速條件下，水密隔膜的轉角曲線趨勢相同，且最大轉角基本保持不變，大致為48°。

圖5 發射裝置及彈丸在不同時刻下的狀態圖

圖6 彈速不同時的水密隔膜轉角-時間曲線

4.2 不同水深條件下的影響

計算了不同水深條件、彈速12 m/s、無艇速條件下，水密隔膜轉角隨時間的曲線如圖7。

圖7 不同水深條件下的水密隔膜轉角曲線

由圖7可以看出，不同水深條件下，水密隔膜轉角曲線的趨勢及峰值均相同，因此可以確定，不同水深條件對水密隔膜轉角基本沒有影響。

4.3 不同艇速條件下的影響

圖8是在水深150 m、彈速12 m/s、不同艇速條件下的水密隔膜轉角隨時間的曲線圖。

圖8 艇速不同時的水密隔膜轉角-時間曲線

將其與圖7所示的無艇速條件下水密隔膜轉角曲線相對比后可以看出：

在t=0-0.15 s左右的運動起始階段，由于發射裝置有艇速，彈丸實際速度比無艇速時大，因此在相同時間下，彈丸所推動的水的體積也相對增大，水密隔膜轉角速度相比無艇速時明顯增加；在t=0.15-0.45 s，為彈身出筒階段，由于彈丸的實際速度為艇速加上推力所造成的速度，彈丸給水密隔膜的力也會相對增大，因此其在48°時，水密隔膜仍會繼續轉動。當其轉動到61°左右，水密隔膜在軸向保持相對平衡狀態，轉角在61°左右基本不變；在t=0.45～0.6 s，彈丸已出筒，由于發射筒內需要海水的補入，因此水密隔膜轉角會有一定程度的振蕩，振蕩過后基本保持穩定。整個發射過程中，水密隔膜轉角的最大值為61°。

5 結論

1) 在給定不同彈速條件下，水密隔膜轉角先迅速增加，后基本不變，隨后由與彈丸出筒時產生的振蕩達到峰值并有所回落，水密隔膜在發射過程中的最大轉角為48°；

2) 在給定不同水深條件下，水密隔膜轉角曲線基本不變，最大轉角為48°；

3) 在給定艇速條件下，水密隔膜轉角速度相比無艇速時有所增大，在轉動61°左右時，水密隔膜在軸向達到相對平衡狀態，并且在達到最大轉角之后，轉角曲線基本保持穩定，整個發射過程中水密隔膜的最大轉角為61°。

由于在建模的過程中考慮到了水密隔膜的安裝方式，和海水水壓的作用，可以為水密隔膜結構模型設計提供一定指導，也可以為魚雷發射裝置設計提供參考。下一步設計應考慮水密隔膜實際結構對于發射裝置的影響。

[1] 李亞男,葉慧娟,張西勇,等.水壓平衡式魚雷發射裝置發射小型運載器內彈道仿真[J].兵器裝備工程學報,2016,37(2):41-44.

[2] 李志華,閻勇,吳吉偉.魚雷發射裝置安全性仿真分析[J].四川兵工學報,2011,32(4):7-8.

[3] 程廣濤,孔巖峰,張振山.液壓平衡式水下武器發射系統仿真分析[J].兵工學報,2009,30(7):915-919.

[4] KUBOTA S,SABURI T,OGATA Y,et al.Numerical Modeling of Shock Initiation in PETN by Unified EOS[J].Science & Technology of Energetic Materials,2010,71(3):92-97.

[5] GREGORY.Failure Modeling of Titanium 6Al-4V and Aluminum 2024-T3 With the Johnson-Cook Material Model,September 2003.

[6] GURTIN M E,MURDOCH A I.A Continuum Theory of Elastic Material Surfaces[J].Archive for Rational Mechanics and Analysis,1975,57(4):291-323.

[7] 朱洪來,白象忠.流固耦合問題的描述方法及分類簡化準則[J].工程力學,2007,24(10):92-99.

[8] AIAA.Dynamics Impact Tolerance of Shuttle Rcc Leading Edge Panels Using Ls-Dyna[C].Aiaa/Asme/Sae/Asee Joint Propulsion Conference and Exhibit.2005.

[9] ZHANG Z,SANG J,ZHANG Y.Inflation Analyses of a Circular Elastomeric Membrane With Hydrostatic Pressure[J].Theoretical and Applied Mechanics Letters,2012,2(6):15-18.

[10] 凡金金.沙粒碰撞平板模型的數值仿真研究[D].南京：南京航空航天大學,2016.

[11] 劉傳龍,張宇文,王亞東,等.提拉活塞式導彈水下發射過程數值模擬[J].上海交通大學學報,2015(5):651-656.

[12] 傅德彬,于殿君,張志勇.潛射導彈離筒后海水倒灌效應數值分析[J].固體火箭技術,2012,35(2):157-160.