某型裝備擊水過程仿真分析

吉禮超，夏慶升，賈 磊，李 寧

(海軍駐重慶地區軍事代表局，重慶401121)

0 引言

目前，具有搭載某型裝備的主戰平臺有很多，但具備某種作戰模式下搭載該型裝備的平臺卻很少。因此，對該型裝備進行改進，使其能滿足某型主戰裝備搭載的需要，能很好地解決某種作戰模式下搭載該型裝備的問題，提高搭載該型裝備的成功率，具有重要的軍事意義。

由于能搭載該型裝備的平臺具有超低空、大速度的能力，突防能力強，導致該型裝備擊水沖擊強度大，擊水速度快，擊水沖擊力大，容易出現損壞的嚴重后果。因此通過加強仿真分析力度，為該型裝備的設計提供有效的數據支持，具有重要意義。

1 LS-DYNA軟件及擊水仿真方案說明

作為世界上最著名的通用顯式動力分程序，LSDYNA軟件能夠模擬真實世界的各種復雜問題，特別適合求解各種二維、三維非線性結構高速碰撞、爆炸和金屬成型等非線性動力沖擊問題，同時可以求解傳熱、流體及流固耦合問題。從理論和算法而言，LS-DYNA是目前所有的顯示求解程序的鼻祖和理論基礎，在工程應用如汽車安全設計、武器系統設計、金屬成型、跌落仿真等領域被廣泛認可為是最佳的分析軟件包。

LS-DYNA軟件以ALE算法為主，兼有Lagrange和Euler算法，以顯式求解為主，兼有隱式求解功能。ALE結合拉格朗日和歐拉網格各自的優點，其時間積分采用中心差分格式，克服了各自的缺點，成為目前非線性連續介質力學中大變形分析的十分先進的有效方法。圖1是運用Ansys/LS-DYNA仿真分析的處理過程。該過程重點是網格劃分和K文件的編寫，同時由于該過程涉及到流固耦合，因此本過程計算也是仿真項目耗時最多的計算任務。

圖1 入水沖擊計算流程圖Fig.1 The shock calculation flow chat of splashing into the water

2 擊水過程仿真分析

本次仿真坐標系有大地坐標系和裝備坐標系。

大地坐標系:原點在裝備頭部，x軸過裝備軸線指向水平方向，y軸指向豎直向上，z軸垂直x-y平面。

裝備坐標系:原點在裝備質心位置，x軸與裝備軸重合指向尾部，y軸在豎直平面內垂直x軸，z軸垂直x-y平面。

選取4個典型條件下該型裝備的初始條件進行仿真分析，其中A為裝備頭部節點，B為裝備尾部節點，仿真結果如下:

圖2 擊水過程有限元模型和坐標系Pig.2 The finite element model and grid system of splashing into the water process

1)平臺速度820 km/h，高度47 m，該型裝備擊水瞬間初始條件見表1。

表1 擊水瞬間初始條件Tab.1 The instant initial condition of splashing into the water

圖3 77°(與垂直方向夾角)入水過程Fig.3 The splashing into the water process at 77°(the angle at vertical direction)

從圖3可以看出，在0.026 s時，裝備尾部撞擊到水袋邊緣，形成二次沖擊，在后面的動力特性曲線中可以體現出來。

圖4 入水瞬間加速度曲線(最大值2 302 m/s2，脈寬0.002 5 s)Fig.4 The acceleration curve at the splashing into the water moment(the most 2 302 m/s2，pulse width 0.002 5 s)

由于受到頭部水沖擊力作用，該作用力產生力矩作用，使裝備產生繞Z軸方向的旋轉，從而使頭部Y方向速度減小，尾部Y方向速度增大。0.025 s時A點的Y向速度減小到了0，隨后產生向上的速度，使裝備抬頭。

尾部加速度曲線在0.026 s時出現第2個峰值1 394 m/s2，表明裝備尾部在該時刻撞擊到水袋邊緣。從裝備姿態變化曲線可看出，裝備姿態角從75°增大到了90°，也就是裝備已經處于水平狀態。

2)平臺速度560 km/h，高度65 m，裝備擊水瞬間初始條件見表2。

表2 擊水瞬間初始條件Tab.2 The instant initial condition of splashing into the water

圖11 67°入水過程Fig.11 The splashing into the water process at 67°

從圖11可清楚地看到擊水過程產生的空化現象。淺色與深色之間的是氣、液二相狀態，不但可以看到尾部擊水，還可看到尾部擊水后卷起的水氣。

從圖12～圖15可以判斷，平臺以560 km/h的速度從72 m的高度進行投放，裝備傾角向垂直方向趨于增大，頭部徑向過載峰值為146個重力加速度(約1 460 m/s2)。

下面是2種極限狀況的擊水過程，一種是平臺以820 km/h的速度從1 000 m高度進行投放，討論該種情況主要是由于980 m投放，裝備下降過程會出現平衡速度，是傘降的一個普遍現象;另一種是以820 km/h的速度從47 m高度進行投放，裝備以47°擊水過程，該過程通過仿真可知很難出現，因為裝備是靜穩設計，以47 m投放時俯仰角在擊水時不會變動如此劇烈，但考慮到擊水過載與擊水速度有很大關系，47 m投放擊水速度最大，為防止極端情況出現，在此處還是對該情況進行了討論，主要提供過載的一種極限參考。首先是980 m擊水過程，然后是47 m，47°擊水過程的仿真分析。

3)投放速度820 km/h，高度980 m，裝備擊水瞬間初始條件見表3。

表3 擊水瞬間初始條件Tab.3 The instant initial condition of splashing into the water

4)平臺速度820 km/h，高度47 m，擊水瞬間初始條件見表4。

表4 擊水瞬間初始條件Tab.4 The instant initial condition of splashing into the water

圖21 47°入水過程Fig.21 The splashing into the water process at 47°

圖26 頭部徑向加速度曲線 (最大值1 297 m/s)Fig.26 The acceleration curve of head at radial direction(the most1 297 m/s)

3 結語

目前基于LS-DYNA做過航空類裝備項目，其入水沖擊與實驗值完全在一個量級上，與美國海軍軍械所提供的檢驗標準也基本一致，具有較高的仿真可信度。通過對以上4個典型條件下該型裝備的初始條件進行仿真分析，可以得出平臺投放該型裝備的最大擊水過載、脈寬、入水姿態角等參數，可以為該型裝備的設計提供有效的數據支持。

［1］問德蓀.工程流體力學(水力學)［M］.北京:高等教育出版社，2004.

［2］韓占忠.FLUENT-流體工程仿真計算實例與應用［M］.北京:北京理工大學出版社，2004.

［3］張國強，吳家敏.流體力學［M］.北京:機械工業出版社，2006.

［4］FINNENORE E J.工程流體力學及其工程應用工程［M］.北京:機械工業出版社，2006.

［5］高殿榮，吳曉明.工程流體力學［M］.北京:機械工業出版社，2004.

［6］沈韶華.工程力學［M］.北京:經濟科學出版社，2010.

［7］原方.工程力學［M］.北京:清華大學出版社，2006.