李欣

摘 要：該文對(duì)導(dǎo)彈氣動(dòng)特性快速估算軟件Datcom的計(jì)算缺陷進(jìn)行修正。結(jié)合某無(wú)翼式導(dǎo)彈的風(fēng)洞試驗(yàn)數(shù)據(jù)，利用國(guó)內(nèi)經(jīng)典的氣動(dòng)特性快速估算資料，以無(wú)翼式布局為例，針對(duì)Datcom軟件法向力和壓心的計(jì)算缺陷進(jìn)行修正，有效提高其計(jì)算精度，提出針對(duì)該類(lèi)氣動(dòng)布局的修正方法。計(jì)算結(jié)果表明，該文提出的修正方法行之有效，修正后的計(jì)算結(jié)果與風(fēng)洞試驗(yàn)結(jié)果較為吻合。

關(guān)鍵詞：氣動(dòng)特性；快速估算；Datcom

中圖分類(lèi)號(hào)：TJ760 文獻(xiàn)標(biāo)志碼：A

0 引言

工程算法能夠借助解析表達(dá)式，快速估算氣動(dòng)特性。如果能夠借助或者參考風(fēng)洞試驗(yàn)和CFD數(shù)據(jù)，對(duì)其中的重要參數(shù)進(jìn)行修正，提高其計(jì)算精度，使之達(dá)到或接近工程應(yīng)用的要求，便可大量節(jié)省試驗(yàn)成本，縮短研制周期。同時(shí)與CFD計(jì)算相比，工程算法簡(jiǎn)易便捷，可作為導(dǎo)彈總體多學(xué)科優(yōu)化的氣動(dòng)學(xué)科代理模型。目前，對(duì)導(dǎo)彈氣動(dòng)特性進(jìn)行計(jì)算，最常用的方法是使用快速估算軟件Missle Datcom。Datcom軟件是美國(guó)NASA通過(guò)大量的風(fēng)洞試驗(yàn)，對(duì)導(dǎo)彈頭部、彈身、彈翼、舵面等各個(gè)部件的氣動(dòng)特性，以及部件之間的耦合因子進(jìn)行數(shù)據(jù)擬合，進(jìn)而得到的導(dǎo)彈氣動(dòng)特性快速計(jì)算程序，氣動(dòng)特性快速計(jì)算的主要參數(shù)主要通過(guò)插值擬合來(lái)獲取。Datcom軟件是目前通用的氣動(dòng)特性快速計(jì)算軟件，可在氣動(dòng)外形方案論證階段對(duì)外形進(jìn)行快速的優(yōu)化設(shè)計(jì)，但大量實(shí)踐經(jīng)驗(yàn)表明，Datcom的計(jì)算仍存在缺陷，而且針對(duì)不同的布局還不一樣，因此對(duì)該軟件的計(jì)算數(shù)據(jù)進(jìn)行“去粗取精”，獲得準(zhǔn)確的結(jié)果以此來(lái)指導(dǎo)實(shí)踐十分重要。

該文以無(wú)翼式布局導(dǎo)彈全彈法向力和壓心計(jì)算為例，分析氣動(dòng)快速估算軟件Datcom對(duì)該類(lèi)氣動(dòng)布局的計(jì)算缺陷，歸納其產(chǎn)生該種缺陷的原因，借助經(jīng)典氣動(dòng)特性工程估算《箭彈空氣動(dòng)力特性分析與計(jì)算》（以下簡(jiǎn)稱(chēng)《箭彈》）中的氣動(dòng)參數(shù)數(shù)據(jù)，建立針對(duì)該類(lèi)布局的計(jì)算修正方法，對(duì)Datcom中的氣動(dòng)參數(shù)進(jìn)行修正，提高其計(jì)算精度，使之與風(fēng)洞試驗(yàn)數(shù)據(jù)吻合。

1 Datcom計(jì)算缺陷

1.1 Datcom計(jì)算缺陷原因分析

通過(guò)與《箭彈》相關(guān)數(shù)據(jù)對(duì)比，可分析出針對(duì)無(wú)翼式布局，Datcom的4條計(jì)算缺陷。

（1）在大攻角、低馬赫區(qū)域，Datcom計(jì)算的壓心系數(shù)明顯偏大。與《箭彈》相比。在低馬赫區(qū)，Datcom的頭部壓心系數(shù)很大，而《箭彈》中的系數(shù)很小，在大攻角、低馬赫段，如果用《箭彈》的頭部壓心數(shù)據(jù)將其替換，就可以使Datcom在低馬赫壓心系數(shù)降低。

（2）在跨聲速段，Datcom計(jì)算的全彈壓心、法向力都偏大。在跨聲速段，與《箭彈》相比，Datcom的舵面法向力偏大，可用《箭彈》的數(shù)據(jù)將其替換，就可以同時(shí)降低法向力和壓心。

（3）在小攻角、高馬赫區(qū)，Datcom計(jì)算的壓心系數(shù)偏大，法向力偏小。（如圖1所示，圖中α為攻角，Ma為導(dǎo)彈飛行馬赫數(shù)）在高馬赫區(qū)，《箭彈》中的頭部法向力系數(shù)明顯比要Datcom中要大，使用《箭彈》中的數(shù)據(jù)與其進(jìn)行互補(bǔ)，可以有效在小攻角、高馬赫段解決法向力偏小，壓心系數(shù)偏大的問(wèn)題。

（4）在大攻角、高馬赫區(qū)，Datcom計(jì)算的壓心系數(shù)偏小，法向力偏小。Datcom的翼身干擾因子數(shù)據(jù)，在高馬赫大攻角段比《箭彈》中偏小，尤其翼身干擾因子KW（B）隨馬赫數(shù)增加下降很明顯。如果用《箭彈》中的數(shù)據(jù)將其替換，就可以同時(shí)增大法向力和壓心系數(shù)。

2.2 Datcom計(jì)算缺陷修正方法

上文在闡述計(jì)算缺陷的同時(shí)，也闡述了修正方法的依據(jù)，下面要介紹具體實(shí)施的過(guò)程。

（1）編寫(xiě)文件讀取程序?qū)atcom后處理文件中的所需參數(shù)讀取出來(lái)。

（2）利用《箭彈》中的數(shù)據(jù)，對(duì)Datcom中需要修正的數(shù)據(jù)進(jìn)行修改，按照《箭彈》中的公式重新計(jì)算法向力和壓心。

依據(jù)之前分析的無(wú)翼式布局的計(jì)算缺陷，經(jīng)過(guò)試探，大概制定了修正方案。

（1）在Ma<1.3并且α>22°，Ma≥4并且10°<α<20°這2個(gè)區(qū)域，用《箭彈》的頭部壓心數(shù)據(jù)；

（2）在Ma>2 并且 α<5°這個(gè)區(qū)域，采用《箭彈》中的頭部法向力數(shù)據(jù)，用下式計(jì)算；

CBN1=CBN+0.8（CBNT-CBN） （1）

其中，CBN1為修正后的頭部法向力數(shù)據(jù)，CBN為Datcom的頭部法向力數(shù)據(jù)，CBNT為《箭彈》中的頭部法向力數(shù)據(jù)；

（3）在1

（4）在Ma*sin（α）>1這個(gè)區(qū)域，采用《箭彈》中的翼身干擾因子數(shù)據(jù)。

3 修正結(jié)果

經(jīng)上述方法修正后，Datcom的計(jì)算精度可得到明顯提高，如圖2所示。圖中“—*—”為修正前的Datcom的計(jì)算值，“— —”是基于以上所提方案修正后的值。

從圖2可以看出，針對(duì)該無(wú)翼式導(dǎo)彈氣動(dòng)外形，采取提出修正方法后，壓心精度有了很明顯的提高。

4 結(jié)論

該文提出的修正方法，可有效修正Datcom對(duì)于無(wú)翼式布局導(dǎo)彈法向力和壓心的計(jì)算缺陷，修正方法的關(guān)鍵，在于歸納出Datcom對(duì)同一類(lèi)氣動(dòng)布局的計(jì)算缺陷，以及產(chǎn)生缺陷的原因，再利用《箭彈》的數(shù)據(jù)對(duì)關(guān)鍵參數(shù)進(jìn)行修正。筆者也歸納過(guò)正常布局的Datcom計(jì)算缺陷，同樣存在一致性，可利用《箭彈》的數(shù)據(jù)來(lái)修正。

參考文獻(xiàn)

[1]本刊編寫(xiě)組.箭彈空氣動(dòng)力特性分析與計(jì)算[M].北京：國(guó)防工業(yè)出版社，1979.

[2]徐敏，安效民.飛行器空氣動(dòng)力特性分析與計(jì)算方法[M].西北工業(yè)大學(xué)出版社，2011.

[3]Willian.Blake， Stanley.Lash. MISSILE DATCOM USERS GUIDE[M].Ohio：Wright Patterson Air Force Base.1997.