丁鍇

摘要：本文采用滑移網(wǎng)格方法，模擬了導(dǎo)彈飛行時舵偏做動的非定常過程。通過非定常過程的模擬計算，得到了導(dǎo)彈動態(tài)過程的氣動特性變化，相較于傳統(tǒng)方法，非定常模擬計算更加接近導(dǎo)彈飛行時的真實情況，并且能夠通過仿真手段得到導(dǎo)彈動導(dǎo)數(shù)數(shù)據(jù)。

關(guān)鍵詞：動網(wǎng)格技術(shù);導(dǎo)彈飛行;數(shù)值模擬

引言

導(dǎo)彈的飛行過程中，會通過轉(zhuǎn)動舵面，產(chǎn)生附加空氣動力，使導(dǎo)彈產(chǎn)生機(jī)動從而達(dá)到擊傷目標(biāo)的目的，因此，計算機(jī)動過程中導(dǎo)彈的瞬時流場具有重要的意義。目前在導(dǎo)彈氣動外形研制中的數(shù)值模擬多以常規(guī)的CFD 方法對特定狀態(tài)下導(dǎo)彈的氣動力進(jìn)行定常模擬分析，這種計算方法的優(yōu)點是計算時間短并且計算結(jié)果相對較為精確，但是定常模擬的方法不能模擬計算導(dǎo)彈的飛行狀態(tài)，并且也不能計算導(dǎo)彈的動導(dǎo)數(shù)信息，因此通過動網(wǎng)格方法對導(dǎo)彈的非定常氣動特性進(jìn)行計算，是對當(dāng)前計算方法的一個有益補(bǔ)充和進(jìn)步。

研究內(nèi)容

為了準(zhǔn)確便捷的對導(dǎo)彈彈道飛行過程進(jìn)行數(shù)值模擬計算，計算中采用了軸對稱導(dǎo)彈模型，采用FLUENT 對導(dǎo)彈舵偏為0度，攻角以6度 / 秒的速率進(jìn)行轉(zhuǎn)動，馬赫數(shù) Ma 為3.01的狀態(tài)進(jìn)行計算。

在計算開始時使用定常計算方法計算馬赫數(shù)3.01下0度攻角狀態(tài)的導(dǎo)彈氣動力，待計算結(jié)果穩(wěn)定后再載入控制網(wǎng)格轉(zhuǎn)動的 UDF 程序，設(shè)置網(wǎng)格轉(zhuǎn)動角速度為6度 / 秒，計算步長1000步 / 秒，在每計算1000步后保存計算結(jié)果，計算5000步后停止計算，計算導(dǎo)彈攻角從0°～30°變化的非定常過程。圖2給出了 Ma=3.01時，不同攻角下流場垂直于 z 軸的中心對稱面壓力云圖，迎風(fēng)面、背風(fēng)面壓力云圖流線圖。

隨著時間的變化，導(dǎo)彈攻角越來越大，在同一馬赫數(shù)下導(dǎo)彈彈身上的壓力分布也隨之變化，并且壓力分布變化符合氣動規(guī)律，滑移網(wǎng)格在計算導(dǎo)彈動態(tài)攻角氣動力變化上是可行的。

通過對比導(dǎo)彈計算數(shù)據(jù)和風(fēng)洞試驗數(shù)據(jù)可以大致得出滑移網(wǎng)格方法的精確度，對比分析結(jié)果如圖3所示：

從圖3可以看出，F(xiàn)LUENT 計算數(shù)據(jù)與風(fēng)洞試驗數(shù)據(jù)相似度較高，證明了此計算方法可行。3.01馬赫數(shù)下兩者符合較好，只是在大攻角時有一定差距。由于 FLUENT 計算數(shù)據(jù)在每一個時間步長也就是0.001s 就有一個數(shù)據(jù)，所以相較與風(fēng)洞試驗，F(xiàn)LUENT 計算數(shù)據(jù)很容易能夠取得動導(dǎo)數(shù)數(shù)據(jù)，并且得到的數(shù)據(jù)相較于傳統(tǒng)的工程算法得到的數(shù)據(jù)更為準(zhǔn)確。

結(jié)論

本文利用 FLUENT 對導(dǎo)彈動態(tài)飛行過程進(jìn)行了模擬計算，從結(jié)果來看，動網(wǎng)格方法可以有效的計算導(dǎo)彈的氣動力數(shù)據(jù)，對導(dǎo)彈在動態(tài)過程中的氣動力變化與響應(yīng)做出有效的模擬。