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

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

0 引言

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

由于能搭載該型裝備的平臺(tái)具有超低空、大速度的能力，突防能力強(qiáng)，導(dǎo)致該型裝備擊水沖擊強(qiáng)度大，擊水速度快，擊水沖擊力大，容易出現(xiàn)損壞的嚴(yán)重后果。因此通過(guò)加強(qiáng)仿真分析力度，為該型裝備的設(shè)計(jì)提供有效的數(shù)據(jù)支持，具有重要意義。

1 LS-DYNA軟件及擊水仿真方案說(shuō)明

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

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

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

2 擊水過(guò)程仿真分析

本次仿真坐標(biāo)系有大地坐標(biāo)系和裝備坐標(biāo)系。

大地坐標(biāo)系:原點(diǎn)在裝備頭部，x軸過(guò)裝備軸線指向水平方向，y軸指向豎直向上，z軸垂直x-y平面。

裝備坐標(biāo)系:原點(diǎn)在裝備質(zhì)心位置，x軸與裝備軸重合指向尾部，y軸在豎直平面內(nèi)垂直x軸，z軸垂直x-y平面。

選取4個(gè)典型條件下該型裝備的初始條件進(jìn)行仿真分析，其中A為裝備頭部節(jié)點(diǎn)，B為裝備尾部節(jié)點(diǎn)，仿真結(jié)果如下:

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

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

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

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

從圖3可以看出，在0.026 s時(shí)，裝備尾部撞擊到水袋邊緣，形成二次沖擊，在后面的動(dòng)力特性曲線中可以體現(xiàn)出來(lái)。

圖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)

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

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

2)平臺(tái)速度560 km/h，高度65 m，裝備擊水瞬間初始條件見(jiàn)表2。

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

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

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

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

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

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

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

4)平臺(tái)速度820 km/h，高度47 m，擊水瞬間初始條件見(jiàn)表4。

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

圖21 47°入水過(guò)程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 結(jié)語(yǔ)

目前基于LS-DYNA做過(guò)航空類裝備項(xiàng)目，其入水沖擊與實(shí)驗(yàn)值完全在一個(gè)量級(jí)上，與美國(guó)海軍軍械所提供的檢驗(yàn)標(biāo)準(zhǔn)也基本一致，具有較高的仿真可信度。通過(guò)對(duì)以上4個(gè)典型條件下該型裝備的初始條件進(jìn)行仿真分析，可以得出平臺(tái)投放該型裝備的最大擊水過(guò)載、脈寬、入水姿態(tài)角等參數(shù)，可以為該型裝備的設(shè)計(jì)提供有效的數(shù)據(jù)支持。

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

［2］韓占忠.FLUENT-流體工程仿真計(jì)算實(shí)例與應(yīng)用［M］.北京:北京理工大學(xué)出版社，2004.

［3］張國(guó)強(qiáng)，吳家敏.流體力學(xué)［M］.北京:機(jī)械工業(yè)出版社，2006.

［4］FINNENORE E J.工程流體力學(xué)及其工程應(yīng)用工程［M］.北京:機(jī)械工業(yè)出版社，2006.

［5］高殿榮，吳曉明.工程流體力學(xué)［M］.北京:機(jī)械工業(yè)出版社，2004.

［6］沈韶華.工程力學(xué)［M］.北京:經(jīng)濟(jì)科學(xué)出版社，2010.

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