彭 科，胡 凡，張為華，周 張

(1.國(guó)防科學(xué)技術(shù)大學(xué) 航天科學(xué)與工程學(xué)院,長(zhǎng)沙 410073;2.江南工業(yè)集團(tuán)有限公司,長(zhǎng)沙 410200)

?

柵格翼氣動(dòng)特性及其應(yīng)用研究綜述

彭 科1，胡 凡1，張為華1，周 張2

(1.國(guó)防科學(xué)技術(shù)大學(xué) 航天科學(xué)與工程學(xué)院,長(zhǎng)沙 410073;2.江南工業(yè)集團(tuán)有限公司,長(zhǎng)沙 410200)

從實(shí)驗(yàn)研究、工程估算、數(shù)值模擬三方面綜述了柵格翼氣動(dòng)特性研究方法及研究現(xiàn)狀；概述了柵格翼主要?dú)鈩?dòng)特性；總結(jié)了弧形、整體后(前)掠、局部后掠等格翼新穎幾何構(gòu)型研究進(jìn)展；全面概述了柵格翼在國(guó)內(nèi)外各類型飛行器上應(yīng)用情況及主要用途。針對(duì)格翼設(shè)計(jì)應(yīng)用過(guò)程的瓶頸問題，分析展望了柵格翼相關(guān)領(lǐng)域需要關(guān)注和解決的關(guān)鍵問題。

柵格翼；氣動(dòng)特性；幾何構(gòu)型；應(yīng)用動(dòng)態(tài)

0 引言

柵格翼是由外部框架和內(nèi)部柵格布置形成的空間多升力面系統(tǒng)[1]，其基本構(gòu)型分為框架式與蜂窩式，蜂窩式又分為正置和斜置。柵格翼起源于飛行器設(shè)計(jì)初期所采用的由多翼面組成的升力系統(tǒng)，但由于多翼面升力系統(tǒng)在當(dāng)時(shí)的技術(shù)條件下存在結(jié)構(gòu)、重量、阻力、工藝等方面的缺陷，長(zhǎng)期以來(lái)飛行器設(shè)計(jì)中單翼升力面占主導(dǎo)地位。隨著技術(shù)的進(jìn)步，人們對(duì)起源于多翼面升力系統(tǒng)的柵格翼的性能特點(diǎn)認(rèn)識(shí)不斷深入，用其作為各類飛行器氣動(dòng)穩(wěn)定面與控制面的優(yōu)勢(shì)不斷凸顯。俄羅斯、美國(guó)、德國(guó)和中國(guó)等都加強(qiáng)了柵格翼相關(guān)研究工作，取得了一系列理論與應(yīng)用成果。

柵格翼幾何結(jié)構(gòu)復(fù)雜，氣動(dòng)問題是其設(shè)計(jì)過(guò)程中的核心關(guān)鍵問題，對(duì)柵格翼氣動(dòng)特性國(guó)內(nèi)外研究進(jìn)展進(jìn)行綜述，把握其研究方法、性能特點(diǎn)、應(yīng)用前景，具有十分重要的意義。本文綜述柵格翼氣動(dòng)特性研究方法、研究成果與國(guó)內(nèi)外應(yīng)用動(dòng)態(tài)，展望柵格翼相關(guān)領(lǐng)域需要關(guān)注和解決的關(guān)鍵問題，為相關(guān)領(lǐng)域研究人員提供參考。

1 柵格翼氣動(dòng)特性研究方法進(jìn)展

柵格翼氣動(dòng)特性研究方法包含實(shí)驗(yàn)、工程估算以及數(shù)值計(jì)算三方面。

1.1 實(shí)驗(yàn)研究

柵格翼外形復(fù)雜，其外流場(chǎng)結(jié)構(gòu)、氣動(dòng)特征、氣動(dòng)分析方法與平面翼差異較大，為掌握柵格翼外流場(chǎng)特征及基本氣動(dòng)特性，建立其氣動(dòng)特性分析方法。國(guó)內(nèi)外從風(fēng)洞實(shí)驗(yàn)、水洞實(shí)驗(yàn)、彈道靶實(shí)驗(yàn)與飛行實(shí)驗(yàn)等方面廣泛開展了柵格翼氣動(dòng)特性研究。

美國(guó)從1985年至上世紀(jì)末密集開展了大量不同構(gòu)形柵格翼亞、跨、超聲速風(fēng)洞實(shí)驗(yàn)，研究了柵格翼基本構(gòu)型氣動(dòng)特性、弧形柵格翼氣動(dòng)特性、跨聲速壅塞現(xiàn)象、減阻構(gòu)型、幾何參數(shù)對(duì)氣動(dòng)特性影響、后掠效應(yīng)、平面翼與柵格翼性能異同等[2]；1996年，美國(guó)進(jìn)行了柵格翼彈頭的火箭助推飛行試驗(yàn)與空投試驗(yàn)，驗(yàn)證了柵格翼作為氣動(dòng)穩(wěn)定部件的良好性能[2]；近年來(lái)，美國(guó)進(jìn)行了“獵戶座”逃逸飛行器縮比模型亞、跨、超聲速氣動(dòng)特性風(fēng)洞實(shí)驗(yàn)研究，驗(yàn)證了柵格翼布局逃逸飛行器靜穩(wěn)定性，對(duì)比了常規(guī)與整體后掠?xùn)鸥褚須鈩?dòng)性能[3]。英國(guó)采用風(fēng)洞實(shí)驗(yàn)研究對(duì)比了柵格翼與平面翼縱橫向靜穩(wěn)定特性、控制特性、鉸鏈力矩特性與阻力特性等方面的異同。國(guó)外還開展了柵格翼彈頭彈道靶實(shí)驗(yàn)研究，測(cè)得了柵格翼彈頭滾轉(zhuǎn)力矩系數(shù)與滾轉(zhuǎn)阻尼系數(shù)[4-6]，得到了流場(chǎng)結(jié)構(gòu)陰影圖像。

國(guó)內(nèi)，南京理工大學(xué)經(jīng)過(guò)多次實(shí)驗(yàn)，研究了亞、跨、超聲速條件下格寬翼弦比對(duì)柵格翼氣動(dòng)特性的影響；探索了減小柵格翼超聲速阻力的方法和途徑；分別研究了采用無(wú)尾翼、“T”型尾翼、柵格尾翼時(shí)鴨式布局遠(yuǎn)程彈在跨/超聲速段的氣動(dòng)特性，得到了不同馬赫數(shù)及攻角狀態(tài)下不同尾翼構(gòu)形對(duì)遠(yuǎn)程彈氣動(dòng)特性尤其是鴨舵滾轉(zhuǎn)控制特性的影響[7-11]。北京航空航天大學(xué)陸中榮教授等采用氫氣泡法和絲線法在水洞和低速風(fēng)洞中進(jìn)行了柵格翼繞流特性實(shí)驗(yàn)[12]，根據(jù)實(shí)驗(yàn)結(jié)果分析了柵格翼自由渦系與各網(wǎng)孔內(nèi)流態(tài)，分析了流動(dòng)機(jī)理。中國(guó)空氣動(dòng)力研究與發(fā)展中心建立了國(guó)內(nèi)首個(gè)全尺寸柵格舵試驗(yàn)平臺(tái)[13]，實(shí)現(xiàn)了全尺寸柵格舵模型氣動(dòng)外形、舵控系統(tǒng)、結(jié)構(gòu)的一體化試驗(yàn)驗(yàn)證；此外，該中心還成功研制出了動(dòng)態(tài)響應(yīng)能力良好、滿足小量鉸鏈力矩測(cè)量要求的天平[14]。中國(guó)航天空氣動(dòng)力技術(shù)研究院研究了柵格翼展開過(guò)程氣動(dòng)力測(cè)量試驗(yàn)方法[15]，實(shí)現(xiàn)了全尺寸柵格翼超聲速條件下自動(dòng)連續(xù)展開。

實(shí)驗(yàn)手段加深了人們對(duì)柵格翼繞流特性與氣動(dòng)特性的理解，研究成果促進(jìn)了工程估算與數(shù)值計(jì)算的不斷發(fā)展，為柵格翼設(shè)計(jì)提供有力支撐。

1.2 工程估算

以大量實(shí)驗(yàn)現(xiàn)象、實(shí)驗(yàn)數(shù)據(jù)為基礎(chǔ)，國(guó)內(nèi)外學(xué)者結(jié)合不同馬赫數(shù)狀態(tài)柵格翼流場(chǎng)結(jié)構(gòu)，根據(jù)部件組拆法，考慮各部件間干擾建立了柵格翼翼身組合體氣動(dòng)特性工程估算模型。

1.2.1 柵格翼流場(chǎng)結(jié)構(gòu)

柵格翼幾何結(jié)構(gòu)復(fù)雜，其繞流結(jié)構(gòu)隨馬赫數(shù)變化而發(fā)生演變，對(duì)其氣動(dòng)特性影響顯著，需根據(jù)不同馬赫數(shù)條件下流場(chǎng)結(jié)構(gòu)特點(diǎn)確定相應(yīng)氣動(dòng)特性工程估算方法[2]。引入臨界馬赫數(shù)表征柵格翼流場(chǎng)結(jié)構(gòu)，不同馬赫數(shù)狀態(tài)柵格翼流場(chǎng)結(jié)構(gòu)與流動(dòng)特征如表1所示。

表1 柵格翼流場(chǎng)結(jié)構(gòu)與流動(dòng)特征

Table 1 Grid fin flow field Sstructure and characteristics

各臨界馬赫數(shù)的計(jì)算是確定特定馬赫數(shù)下柵格翼流場(chǎng)結(jié)構(gòu)的關(guān)鍵。文獻(xiàn)[2]將柵格翼內(nèi)部流動(dòng)等效為二維噴管流動(dòng)建立了柵格翼第一臨界馬赫數(shù)計(jì)算方法。第二臨界馬赫數(shù)Macr2與格片前緣半徑等參數(shù)相關(guān)，計(jì)算較為復(fù)雜，目前多依靠實(shí)驗(yàn)數(shù)據(jù)比擬得到。第三臨界馬赫數(shù)可根據(jù)格片鍥角采用理論公式計(jì)算，當(dāng)格片為鈍前緣時(shí)采用等效鍥方法計(jì)算，精度要求較高時(shí)可采用數(shù)值方法計(jì)算[16]。

1.2.2 單獨(dú)柵格翼氣動(dòng)特性工程估算模型

對(duì)于單獨(dú)柵格翼，其氣動(dòng)力包含表面法向壓力與切向摩擦力兩部分貢獻(xiàn)。摩擦力貢獻(xiàn)目前多采用經(jīng)驗(yàn)公式計(jì)算[17-18]；表面法向壓力貢獻(xiàn)采用基于無(wú)粘假設(shè)的工程模型計(jì)算，本節(jié)對(duì)該領(lǐng)域國(guó)內(nèi)外研究方法進(jìn)行綜述。

亞聲速時(shí)，根據(jù)實(shí)驗(yàn)結(jié)果，由柵格翼格片脫出的自由渦系形成互相平行的旋渦壁[12,19]，該類流動(dòng)的氣動(dòng)力計(jì)算適于采用基于升力面理論的渦格法。文獻(xiàn)[18,20-23]采用該方法計(jì)算了不同類型柵格翼亞聲速氣動(dòng)特性，文獻(xiàn)[18]對(duì)馬赫數(shù)影響作了修正。此外，臧勇基于第二代低階面元法VSAERO，建立了亞聲速、小攻角狀態(tài)飛行器流場(chǎng)模型與氣動(dòng)載荷計(jì)算模型，進(jìn)而計(jì)算了柵格翼亞聲速氣動(dòng)特性[24]。

柵格翼跨聲速段格孔內(nèi)部流動(dòng)復(fù)雜，氣動(dòng)特性工程估算難度較大，目前國(guó)內(nèi)外相關(guān)研究較少，計(jì)算理論尚不成熟。Macr1≤Ma≤1情況下多采用考慮邊界層修正的一維流理論計(jì)算[25]，在1≤Ma≤Macr2根據(jù)理論公式計(jì)算脫體激波前后氣流參數(shù)，近似認(rèn)為波后流動(dòng)方向不變，根據(jù)波后參數(shù)采用亞聲速狀態(tài)方法計(jì)算柵格舵所受氣動(dòng)力。

關(guān)于柵格翼超聲速氣動(dòng)特性工程估算，沈遐齡等不考慮格片間激波干擾，根據(jù)各格片理論升力線斜率求得正置柵格翼超聲速氣動(dòng)特性[20]。文獻(xiàn)[26-27]給出了斜置蜂窩式柵格翼對(duì)正置蜂窩式柵格翼的等效條件。文獻(xiàn)[25]根據(jù)平面翼超聲速氣動(dòng)特性計(jì)算模型，考慮端板效應(yīng)，求得柵格翼各格片表面壓力分布，進(jìn)而得到其超聲速氣動(dòng)特性。Theerthamalai P等考慮格片間激波干擾，基于激波-膨脹波理論發(fā)展了較為完整的柵格翼超聲速氣動(dòng)特性工程估算方法[17]。文獻(xiàn)[28]基于面元法建立柵格翼翼身組合體高超聲速氣動(dòng)特性計(jì)算模型。

1.2.3 柵格翼翼身組合體部件干擾模型

對(duì)于柵格翼翼身組合體，需考慮彈身對(duì)柵格翼干擾、柵格翼對(duì)彈身干擾、柵格翼各格片間干擾。彈身對(duì)柵格翼的氣動(dòng)干擾計(jì)算方法常用的方法有：偶極子模擬彈身法、上洗流場(chǎng)模型法、平面渦格法，其中偶極子模擬彈身法所做的簡(jiǎn)化最少，精度最高，是用來(lái)計(jì)算柵格翼與彈身間相互氣動(dòng)干擾的最佳方法；上洗流場(chǎng)模型法、平面渦格法計(jì)算較簡(jiǎn)單，在對(duì)計(jì)算精度要求不高的情況下可采用[18]。對(duì)大攻角狀態(tài)彈體脫體渦影響，目前國(guó)內(nèi)外多采用經(jīng)驗(yàn)公式估算[18,22]。柵格翼對(duì)彈身干擾主要基于鏡像法與經(jīng)驗(yàn)公式計(jì)算[17]。柵格翼各格片間干擾在亞聲速狀態(tài)需做考慮，文獻(xiàn)[25]采用1.2.2節(jié)方法將各格片同時(shí)劃分面元聯(lián)合求解，所得結(jié)果即反應(yīng)了柵格舵舵片間干擾，跨、超聲速段小型運(yùn)載火箭各舵片間干擾可忽略不計(jì)。

現(xiàn)有工程估算方法效率較高，但精度不能滿足柵格翼外形設(shè)計(jì)與氣動(dòng)特性仿真等工程需要，相關(guān)理論方法尚待進(jìn)一步研究。

1.3 數(shù)值模擬

基于大量假設(shè)與簡(jiǎn)化的工程估算方法存在精度瓶頸，數(shù)值模擬方法被引入到柵格翼氣動(dòng)特性研究。該方法能用較少假設(shè)給出詳盡、準(zhǔn)確、廣泛的信息，在柵格翼氣動(dòng)問題研究過(guò)程中發(fā)揮了不可替代的作用，拓展了柵格翼氣動(dòng)特性研究方式。

關(guān)于柵格翼氣動(dòng)特性數(shù)值模擬，國(guó)內(nèi)外公開發(fā)表的文獻(xiàn)所取數(shù)學(xué)模型、求解方法多為常規(guī)CFD方法。其幾何結(jié)構(gòu)復(fù)雜，數(shù)值模擬難度主要集中于網(wǎng)格劃分。目前國(guó)內(nèi)外柵格翼翼身組合體氣動(dòng)特性數(shù)值模擬網(wǎng)格劃分方法主要包括：多塊對(duì)接結(jié)構(gòu)網(wǎng)格[29]、四面體非結(jié)構(gòu)網(wǎng)格[30]、直角非結(jié)構(gòu)網(wǎng)格[31]、結(jié)構(gòu)/非結(jié)構(gòu)混合網(wǎng)格[32-33]重疊結(jié)構(gòu)網(wǎng)格[34]等，后3種方法最為常用。直角非結(jié)構(gòu)網(wǎng)格劃分方法簡(jiǎn)單、省時(shí)、自動(dòng)化程度較高，便于網(wǎng)格自適應(yīng)處理；結(jié)構(gòu)/非結(jié)構(gòu)混合網(wǎng)格方法結(jié)合了結(jié)構(gòu)網(wǎng)格和非結(jié)構(gòu)網(wǎng)格各自的優(yōu)點(diǎn)，減小網(wǎng)格生成難度，提高網(wǎng)格質(zhì)量，計(jì)算精度較高；重疊結(jié)構(gòu)網(wǎng)格劃分方法可有效提高計(jì)算精度。

數(shù)值模擬方法對(duì)于柵格翼的氣動(dòng)特性預(yù)示精度可滿足工程需求，但計(jì)算效率較低。

2 柵格翼氣動(dòng)特性研究成果綜述

2.1 柵格翼主要?dú)鈩?dòng)特性

柵格翼作為新型氣動(dòng)部件，在升力特性、鉸鏈力矩特性、折疊特性、強(qiáng)度特性等方面較傳統(tǒng)平面翼具有不可替代的優(yōu)點(diǎn)，主要性能特點(diǎn)為[38]：

(1)升力特性好，大攻角狀態(tài)仍能保持良好線性，不易失速；

(2)壓心漂移小，可有效減小鉸鏈力矩，降低對(duì)舵機(jī)的要求，減輕舵機(jī)質(zhì)量；

(3)弦向尺寸小，可緊貼彈體折疊安裝，減小尺寸，利于存儲(chǔ)，折疊狀態(tài)可在自身空氣動(dòng)力矩作用下自動(dòng)打開；

(4)強(qiáng)度質(zhì)量比大，有利于降低結(jié)構(gòu)質(zhì)量；

(5)阻力較大，可用作減速部件，合理設(shè)計(jì)降低阻力后可拓寬使用范圍。

此外，柵格翼還有以下特點(diǎn)：

(1)作為空間多升力面系統(tǒng)，柵格尾翼縱向格片對(duì)升力仍有貢獻(xiàn)，比平板尾翼利于提高全箭靜穩(wěn)定性，可提供更大法向過(guò)載，提高機(jī)動(dòng)性能。值得關(guān)注的是，背風(fēng)區(qū)柵格翼大攻角時(shí)其根部出現(xiàn)流動(dòng)分離，導(dǎo)致部分升力面當(dāng)?shù)毓ソ菫樨?fù)值，對(duì)升力的貢獻(xiàn)低于其他翼面[30];

(2)對(duì)于采用鴨舵控制滾轉(zhuǎn)的火箭/導(dǎo)彈，在攻角非零情況下，受脫體渦影響，尾翼為平面翼時(shí)易出現(xiàn)滾轉(zhuǎn)控制耦合問題，甚至產(chǎn)生控制反效，采用柵格翼時(shí)滾轉(zhuǎn)控制更為有效，為解決鴨式布局滾轉(zhuǎn)耦合問題提供了思路[39];

(3)在大攻角狀態(tài)柵格尾翼仍可保持偏航方向穩(wěn)定性，也具備大攻角狀態(tài)橫向機(jī)動(dòng)的能力；

(4)柵格翼跨聲速階段產(chǎn)生氣流壅塞現(xiàn)象[40-41]，對(duì)其氣動(dòng)特性影響顯著，控制效率降低[38]，阻力增加，且第一臨界馬赫數(shù)附近法向力系數(shù)、俯仰力矩系數(shù)、壓心系數(shù)變化皆較大，對(duì)飛行器靜穩(wěn)定性產(chǎn)生較大影響[42-43]。

2.2 柵格翼幾何構(gòu)型研究進(jìn)展

對(duì)于柵格翼基本構(gòu)形，影響其氣動(dòng)性能的主要因素有邊框厚度、格片厚度、柵格長(zhǎng)寬比、格寬翼弦比與翼型等。陳少松等采用風(fēng)洞實(shí)驗(yàn)方法對(duì)不同形狀柵格翼的減阻特性進(jìn)行了分析，得到了柵格翼邊框尺寸、邊框截面形狀、格柵形式、格寬翼弦比等對(duì)柵格翼阻力的影響，結(jié)果表明，(1)柵格翼的邊框剖面形狀對(duì)柵格翼的阻力影響較大，遠(yuǎn)大于格片厚度帶來(lái)的影響；(2)減少柵格數(shù)目的情況下，柵格翼阻力降低，同時(shí)帶來(lái)升力的損失；(3)柵格翼格寬翼弦比通過(guò)影響內(nèi)部流態(tài)影響其升阻力特性，在一定馬赫數(shù)、攻角狀態(tài)下，存在最佳格寬翼弦比使得柵格翼升阻比最大[8-10]。吳曉軍等運(yùn)用數(shù)值模擬方法研究了柵格長(zhǎng)寬比對(duì)柵格氣動(dòng)特性的影響，結(jié)果表明，在所選擇的長(zhǎng)寬比范圍內(nèi)，長(zhǎng)寬比對(duì)阻力的影響較小，而對(duì)法向力影響明顯，可根據(jù)需要優(yōu)化長(zhǎng)寬比[44]。

隨著研究的深入，出現(xiàn)了弧形柵格翼、整體后掠?xùn)鸥褚砼c前緣局部后掠?xùn)鸥褚淼确浅Ｒ?guī)構(gòu)形柵格翼，提升了性能。

2.2.1 弧形柵格翼

文獻(xiàn)[45]介紹了弧形柵格翼概念，弧形柵格翼折疊安裝可與彈體完全貼合，最大限度減小尺寸?；⌒螙鸥褚砜稍O(shè)計(jì)為凸面迎風(fēng)與凹面迎風(fēng)兩種形式。兩種形式的弧形柵格翼在跨聲速段升阻比與平面柵格翼差異不大，亞、超聲速段皆高于平面柵格翼[46]；迎風(fēng)面為凸面時(shí)其升力系數(shù)與升阻比皆略大于迎風(fēng)面為凹面的情況；迎風(fēng)面為凹面時(shí)弧形柵格翼展開所需機(jī)構(gòu)力與其所受氣動(dòng)力方向一致，便于展開。目前弧形柵格翼設(shè)計(jì)多以凹面為迎風(fēng)面(圖1)。

圖1 弧形柵格翼Fig.1 Curved grid fin

2.2.2 整體后(前)掠?xùn)鸥褚?/p>

文獻(xiàn)[45]引出了整體后掠與整體前掠?xùn)鸥褚砀拍?，各自又分為“斜”與“平”2種類型(圖2、圖3)。

“斜型”柵格翼格片平面與彈體軸線夾角δ等于翼面后掠角Λ，來(lái)流攻角為零時(shí)格片當(dāng)?shù)毓ソ遣粸榱?。?duì)于“斜型”柵格翼，當(dāng)后(前)掠角一致時(shí)，后掠與前掠方式升、阻力特性基本一致，在前/后掠角在約±30°以內(nèi)時(shí)，法向力基本不隨后(前)掠角變化，30°后(前)掠角狀態(tài)阻力約為零后(前)掠角狀態(tài)的2～3倍，“斜型”后掠?xùn)鸥褚砜稍诒３诛w行器靜穩(wěn)定性的同時(shí)用作減速部件[45]?！捌叫汀焙?前)掠?xùn)鸥褚砜捎行p少阻力，后掠狀態(tài)較前掠狀態(tài)減阻效果更好?！矮C戶座”逃逸飛行器采用“平型”整體后掠?xùn)鸥褚硖岣哽o穩(wěn)定性能[3]。

(a)斜型 (b)平型

(a)斜型 (b)平型

2.2.3 前緣局部后掠?xùn)鸥褚?/p>

受平面翼后掠減阻效應(yīng)啟發(fā)，Guyot等提出了局部后掠?xùn)鸥褚砀拍頪47](圖4)，其零升阻力較常規(guī)構(gòu)型柵格翼約減少30%～40%，比整體后掠?xùn)鸥褚頊p阻效果更好，并可有效減緩跨聲速段壅塞現(xiàn)象，減小阻力[48-49]。局部后掠?xùn)鸥褚矸譃闁鸥窠唤狱c(diǎn)后掠(模型GP)與格片中心點(diǎn)后掠兩類(模型GV)，前者減阻更為有效[11]。

前緣局部后掠?xùn)鸥褚硪话阋硐逸^長(zhǎng)，一定程度上影響了其折疊性能，另外加工更為復(fù)雜，目前尚未見前緣局部后掠?xùn)鸥褚淼南嚓P(guān)工程應(yīng)用。

(a)常規(guī)構(gòu)型

(b)后掠?xùn)鸥褚?/p>

3 柵格翼國(guó)內(nèi)外應(yīng)用動(dòng)態(tài)

柵格翼因其獨(dú)特優(yōu)勢(shì)引起了世界各國(guó)的重視，前蘇聯(lián)/俄羅斯、美國(guó)、德國(guó)以及中國(guó)等都加強(qiáng)了柵格翼的理論與試驗(yàn)研究工作，將其成功用于各類飛行器[50]。

3.1 前蘇聯(lián)/俄羅斯

前蘇聯(lián)于20世紀(jì)50年代初就對(duì)柵格翼進(jìn)行了研究，進(jìn)行了大量氣動(dòng)計(jì)算和實(shí)驗(yàn)，開展了柵格翼結(jié)構(gòu)、強(qiáng)度、生產(chǎn)工藝的研究，形成了柵格翼設(shè)計(jì)方法，并成功將柵格翼用作聯(lián)盟號(hào)宇宙飛船救生逃逸系統(tǒng)的穩(wěn)定、減速部件，1963年實(shí)現(xiàn)了首次發(fā)射。N-1探月火箭一級(jí)采用柵格翼作為穩(wěn)定面，在轉(zhuǎn)運(yùn)過(guò)程中柵格翼折疊，1969年進(jìn)行了首次發(fā)射。

70年代前后，前蘇聯(lián)科學(xué)家開始在各類導(dǎo)彈設(shè)計(jì)中廣泛應(yīng)用柵格翼。戰(zhàn)略彈道導(dǎo)彈中，野人(Savage，北約代號(hào)SS-13)I型與II型一級(jí)使用4個(gè)柵格翼作為穩(wěn)定面；罪人(Sinner，北約代號(hào)SS-16)、佩刀(Saber，北約代號(hào)SS-20)I型與II型、鐮刀(Sickle，北約代號(hào)SS-25)、鐮刀 B(Sickle B，北約代號(hào)SS-27)一級(jí)尾段有8個(gè)柵格翼，其中4個(gè)為穩(wěn)定面，另外4個(gè)為控制面。

戰(zhàn)術(shù)彈道導(dǎo)彈中，薄板(Scaleboard，北約代號(hào)SS-12，兩級(jí)固體導(dǎo)彈)二級(jí)尾部安裝有4個(gè)弧形柵格翼，一級(jí)飛行時(shí)柵格翼折疊，以減小一級(jí)飛行靜不穩(wěn)定度，改善一級(jí)操縱性，二級(jí)飛行時(shí)柵格翼展開，并作為控制面與燃?xì)舛媛?lián)動(dòng)；圣甲蟲(Scarab，北約代號(hào)SS-21)、蜘蛛(Spider，北約代號(hào)SS-23)采用可折疊柵格翼，作為控制面與燃?xì)舛媛?lián)動(dòng)。

空空導(dǎo)彈中蝰蛇(Adder，北約代號(hào)AA-12，俄羅斯稱R-77)及其各改進(jìn)型采用可折疊柵格翼作為控制面，R-77導(dǎo)彈是世界上首型采用柵格翼作為控制面的空空導(dǎo)彈，充分利用了柵格翼各方面優(yōu)點(diǎn)，是目前柵格翼最成功的應(yīng)用案例[50]。對(duì)于空射導(dǎo)彈，隨著現(xiàn)代戰(zhàn)爭(zhēng)對(duì)隱身性能的需求不斷增強(qiáng)，機(jī)載武器、傳感器等掛載實(shí)現(xiàn)內(nèi)埋是該類武器系統(tǒng)發(fā)展的必然趨勢(shì)，使用可折疊的柵格翼是減小空射導(dǎo)彈體積、實(shí)現(xiàn)內(nèi)埋的有效手段[51]。

俄羅斯最新的潛射巡航導(dǎo)彈熾熱(Sizzler，北約代號(hào)SS-N-27，俄羅斯稱3M-54 Klub)助推級(jí)采用折疊柵格翼，待發(fā)狀態(tài)緊貼彈體折疊，發(fā)射后柵格翼打開，可有效解決導(dǎo)彈水中潛行靜穩(wěn)定問題。

3.2 美國(guó)

美國(guó)柵格翼相關(guān)研究始于20世紀(jì)80年代中后期，1991年產(chǎn)生了第一個(gè)關(guān)于柵格翼技術(shù)的專利[51]。美國(guó)在MOAB(Massive Ordnance Air Blast bomb)炸彈上實(shí)現(xiàn)了柵格翼的首次應(yīng)用，將柵格翼用作穩(wěn)定面與控制面，MOAB炸彈于2003年3月11日開展了首次試驗(yàn)。此外，美國(guó)在MOP(Massive Ordnance Penetrator)侵徹彈中同樣使用了柵格翼。美國(guó)空軍干擾彈MALD(Miniature Air-Launched Decoy)采用復(fù)雜構(gòu)形可折疊柵格翼，可從B-52轟炸機(jī)、F-16戰(zhàn)斗機(jī)等平臺(tái)發(fā)射以干擾對(duì)方防空力量[52]。為提高反潛作戰(zhàn)能力，解決聲吶浮標(biāo)的高空準(zhǔn)確投放問題，美國(guó)于2005年開展了相關(guān)技術(shù)研究，設(shè)計(jì)了空投飛行器(ADV，Air Delivery Vehicle)，該飛行器在高空從飛機(jī)上投放后，可在柵格翼控制下機(jī)動(dòng)到指定位置并釋放聲吶浮標(biāo)，該技術(shù)也可用于其他作戰(zhàn)環(huán)境下的載荷空投。此外，美國(guó)在“獵戶座”逃逸飛行器[3]及無(wú)逃逸塔的逃逸飛行器方案[53]中同樣使用了柵格翼。

3.3 其他

德國(guó)也開展了柵格翼的理論研究和試驗(yàn)研制工作。2002年2月，德國(guó)試飛了其首款使用柵格翼進(jìn)行姿態(tài)控制的高超聲速導(dǎo)彈HFK-E0，其飛行馬赫數(shù)達(dá)6.66。

國(guó)內(nèi)自20世紀(jì)90年代初開始了對(duì)柵格翼的研究，“長(zhǎng)征二號(hào)F”型運(yùn)載火箭逃逸飛行器使用柵格翼作為氣動(dòng)穩(wěn)定部件，快舟小型固體運(yùn)載火箭在國(guó)內(nèi)實(shí)現(xiàn)了柵格舵作為氣動(dòng)控制部件的首次應(yīng)用，標(biāo)志著我國(guó)柵格翼設(shè)計(jì)應(yīng)用水平邁上了新的臺(tái)階。

4 關(guān)鍵技術(shù)問題分析與展望

在柵格翼獨(dú)特性能優(yōu)勢(shì)的驅(qū)動(dòng)下，國(guó)內(nèi)外關(guān)于其氣動(dòng)特性研究與應(yīng)用實(shí)踐方面都取得一系列的進(jìn)展與突破，體現(xiàn)在以下方面：

(1)氣動(dòng)特性研究手段不斷發(fā)展。實(shí)驗(yàn)研究投入不斷增大，為準(zhǔn)確把握其流場(chǎng)結(jié)構(gòu)與性能特點(diǎn)發(fā)揮了不可替代的作用；初步建立了柵格翼氣動(dòng)特性工程估算方法；采用數(shù)值求解方法在其性能與新穎構(gòu)型研究方面取得諸多成果。

(2)對(duì)于其氣動(dòng)特性的認(rèn)識(shí)不斷深入，研究提出了弧形、全局后掠、局部后掠等新穎柵格翼構(gòu)型，為進(jìn)一步提高柵格翼性能提供了新的思路。

(3)在多個(gè)國(guó)家、不同類型飛行器、多種型號(hào)上得以成功應(yīng)用。

為進(jìn)一步發(fā)揮柵格翼作為氣動(dòng)穩(wěn)定面與控制面的性能優(yōu)勢(shì)，克服其缺點(diǎn)，以下方面值得進(jìn)一步關(guān)注與研究：

(1)柵格翼氣動(dòng)特性工程計(jì)算方法研究。工程計(jì)算方法效率較高，是柵格翼初步設(shè)計(jì)過(guò)程不可或缺的仿真手段?；诖罅繉?shí)驗(yàn)與數(shù)值仿真結(jié)果修正現(xiàn)有工程計(jì)算模型，建立其氣動(dòng)特性計(jì)算經(jīng)驗(yàn)?zāi)Ｐ?，是亟待開展的工作。

(2)柵格翼外形設(shè)計(jì)準(zhǔn)則研究。目前公開發(fā)表的文獻(xiàn)鮮有直接面向柵格翼設(shè)計(jì)的成果，亟需建立柵格翼格數(shù)、格柵形式、格寬翼弦比等基本參數(shù)設(shè)計(jì)準(zhǔn)則，加強(qiáng)柵格翼設(shè)計(jì)模型、流程、方法研究。

(3)跨聲速段動(dòng)態(tài)不穩(wěn)定性研究。柵格翼在跨聲速段特定馬赫數(shù)條件下俯仰力矩距離變化，動(dòng)態(tài)不穩(wěn)定問題突出，采用非定常數(shù)值模擬等手段對(duì)其產(chǎn)生過(guò)程、產(chǎn)生機(jī)理進(jìn)行研究，對(duì)降低柵格翼應(yīng)用成本與應(yīng)用風(fēng)險(xiǎn)具有現(xiàn)實(shí)意義。

(4)隱身性能研究。柵格翼隱身性能較差，開展其隱身特性定量研究對(duì)其應(yīng)用具有指導(dǎo)意義。

(5)高超聲速狀態(tài)氣動(dòng)熱特性研究、氣動(dòng)力/熱耦合計(jì)算研究、防熱結(jié)構(gòu)設(shè)計(jì)。柵格翼優(yōu)良的高超聲速性能必將帶來(lái)其高超聲速領(lǐng)域應(yīng)用的不斷拓展，但高超聲速條件下其復(fù)雜繞流結(jié)構(gòu)導(dǎo)致的氣動(dòng)熱問題同樣突出，開展高超聲速狀態(tài)氣動(dòng)熱特性研究、氣動(dòng)力/熱耦合計(jì)算研究與防熱結(jié)構(gòu)設(shè)計(jì)研究重要理論與實(shí)踐意義。

(6)結(jié)構(gòu)與加工工藝研究。柵格翼加工工藝較為復(fù)雜，高超聲速條件下防熱結(jié)構(gòu)的引入使得該問題更為突出，在當(dāng)前最新工藝成果的基礎(chǔ)上開展其加工工藝研究對(duì)推動(dòng)?xùn)鸥褚響?yīng)用領(lǐng)域的拓展具有重要意義。

5 結(jié)束語(yǔ)

隨著相關(guān)領(lǐng)域技術(shù)的不斷發(fā)展，柵格翼作為氣動(dòng)穩(wěn)定面與控制面的獨(dú)特性能優(yōu)勢(shì)必將使其應(yīng)用領(lǐng)域不斷拓展。結(jié)合運(yùn)載火箭、彈道導(dǎo)彈、巡航導(dǎo)彈、智能彈藥、再入飛行器、水下航行體等多種不同對(duì)象，開展柵格翼應(yīng)用研究必將產(chǎn)生巨大社會(huì)與軍事效益。

[1] 吳甲生，雷娟棉.制導(dǎo)兵器氣動(dòng)布局與氣動(dòng)特性[M].北京:國(guó)防工業(yè)出版社,2008.

[2] Wm David Washington,Mark S Miller.Experimental investigations of grid fin aerodynamics:a synopsis of nine wind tunnel and three flight tests[C]//RTO AVT Symposium on Missile Aerodynamics.Sorrento,Italy:1998.

[3] Daniel A Pruzan,Michael R Mendenha,William C Rose,David M Schuster.Grid fin stabilization of the orion launch abort vehicle[C]//29th AIAA Applied Aerodynamics Conference.Honolulu,Hawaii:2011.

[4] Dupuis A,Berner C.Aerodynamic aspects of a grid finned projectile at subsonic and supersonic velocities[C]//19th International Symposium of Ballistics.Interlaken,Switzerland:2001.

[5] Berner C,Dupuis A.Wind tunnel tests of a grid finned projectile configuration[C]//39th Aerospace Sciences Meeting & Exhibit.Reno,Nevada:2001.

[6] Fournier E Y.Wind tunnel investigation of grid fin and conventional planar control surfaces[C]//39th Aerospace Sciences Meeting & Exhibit.Reno,Nevada:2001.

[7] 鄧帆，陳少松，譚獻(xiàn)忠，等.不同尾翼鴨式布局遠(yuǎn)程彈在跨/超聲速的氣動(dòng)特性實(shí)驗(yàn)研究[J].實(shí)驗(yàn)流體力學(xué)，2010,24(2):46-50.

[8] 陳少松，張建葉.格寬翼弦比對(duì)柵格翼氣動(dòng)特性的影響[J].南京理工大學(xué)學(xué)報(bào),2005,29(6):693-696.

[9] 陳少松，徐琴，王福華，等.格柵翼減阻特性研究[J].流體力學(xué)實(shí)驗(yàn)與測(cè)量， 2001,15(4):7-11.

[10] 陳少松，徐琴，王福華，等.格柵翼組合體的超音速氣動(dòng)特性研究[J].彈道學(xué)報(bào)， 2000,12(2):50-54.

[11] 鄧帆，陳少松.柵格翼外形特征對(duì)減阻影響的研究[J].實(shí)驗(yàn)流體力學(xué)，2011,25(3):10-15.

[12] 陸中榮，王海文.柵格翼繞流特性的實(shí)驗(yàn)研究[J].流體力學(xué)實(shí)驗(yàn)與測(cè)量，2002,16(1):21-26.

[13] 陳建中，趙忠良，涂正光，等.柵格舵氣動(dòng)與操縱特性高速風(fēng)洞試驗(yàn)技術(shù)研究[J].實(shí)驗(yàn)流體力學(xué)，2011,25(6):82-87.

[14] 陳建中，易國(guó)慶，彭超，等.全尺寸柵格舵鉸鏈力矩天平研制[J].實(shí)驗(yàn)流體力學(xué)， 2012,26(1):79-82.

[15] 吳凱，張尚彬，蔣坤，等.柵格翼展開過(guò)程測(cè)力試驗(yàn)方法[C]//中國(guó)空氣動(dòng)力學(xué)會(huì)測(cè)控技術(shù)專委會(huì)第六屆四次學(xué)術(shù)交流會(huì)文集.三亞:2013.

[16] 張亮，王淑華，姜貴慶.鈍化前緣對(duì)柵格翼激波干擾與熱流分布的影響[J].宇航學(xué)報(bào),2010,31(2):324-328.

[17] Theerthamalai P,Nagarathinam M.Aerodynamic analysis of grid-fin configurations at supersonic speeds[J].Journal of Spacecraft and Rockets,2006,43(4):750-756.

[18] Theerthamalai P.Aerodynamic characterization of grid fins at subsonic speeds[J].Journal of Aircraft,2007,44(2):694-697.

[19] LU Zhong-rong,SU Wen-han.Hairpin vortices in the wake of grid fin[J].Chinese Journal of Aeronautics,2001,14(4):200-204.

[20] 沈遐齡，王麗麗.柵格翼和機(jī)身組合體的氣動(dòng)特性計(jì)算分析[J].北京航空航天大學(xué)學(xué)報(bào),2001,27(1):54-56.

[21] 余永剛.亞跨音速柵格翼氣動(dòng)特性研究[D].南京:南京理工大學(xué),2007.

[22] John E Burkhalter,Roy J Hartfield,Todd M Leleuxt.Nonlinear aerodynamic analysis of grid fin configurations[J].Journal of Aircraft,1995,32(3):547-554.

[23] Ross A Brooks,John E Burkhaltert.Experimental and analytical analysis of grid fin configurations[J].Journal of Aircraft,1989,26(9):885-887.

[24] 臧勇.柵格舵/彈身組合體亞音速氣動(dòng)特性研究[D].北京:北京理工大學(xué),2000.

[25] Richard W Kretzschmar,John E Burkhalter.Aerodynamic prediction methodology for grid fins[C]//RTO AVT Symposium on Missile Aerodynamics.Sorrento,Italy:1998.

[26] 任懷宇，蔡昱，孫洪森.高超聲速柵格翼氣動(dòng)特性研究[J].宇航學(xué)報(bào),2010,31(1):44-49.

[27] 陸中榮，沈遐齡，童自力,等.柵格翼空氣動(dòng)力特性的計(jì)算與分析[J].北京航空航天大學(xué)學(xué)報(bào),1996,22(5):575-580.

[28] 周張，胡凡，于勇.基于面元法的柵格翼翼身組合體高超聲速氣動(dòng)特性研究[J].固體火箭技術(shù),2012,35(4):442-446.

[29] Sun Y,Khalid M.A Cfd investigation of grid fin missiles[C]//34th AIAAlASME/SAE/ASEE Joint Propulsion Conference & Exhibit.Cleveland,OH:1998.

[30] James Despirito,Jr Milton E Vaughn,W David Washington.Cfd investigation of canard-controlled missile with planar and grid fins in supersonic flow[C]//AIAA Atmospheric Flight Mechanics Conference and Exhibit.Monterey,California:2002.

[31] James E Kless,Michael J Aftosmis.Analysis of grid fins for launch abort vehicle using a cartesian euler solver[C]//29th AIAA Applied Aerodynamics Conference.Honolulu,Hawaii:2011.

[32] Deng You-qi,Ma Ming-sheng,Zheng Ming.Navier-stokes computation of grid fin missile using hybrid structured-unstructured grids[J].Chinese Journal of Aeronautics,2006,19(4):304-308.

[33] 劉剛，肖中云，江雄,等.混合網(wǎng)格方法在柵格翼數(shù)值模擬中的應(yīng)用研究[J].空氣動(dòng)力學(xué)學(xué)報(bào),2007,25(2):230-277.

[34] 王學(xué)占，王立強(qiáng).基于結(jié)構(gòu)重疊網(wǎng)格的柵格翼導(dǎo)彈流場(chǎng)數(shù)值仿真[J].航空兵器,2012(1):44-47.

[35] Montgomery C Hughson,Eric L Blades,Edward A Luke,Gregg L Abate.Analysis of lattice grid tailfin missiles in high-speed flow[C]//25th AIAA Applied Aerodynamics Conference.Miami,FL:2007.

[36] Philippe Reynier,Uwe Reisch,José-Maria Longo,Rolf Radespiel.Flow predictions around a missile with lattice wings using the actuator disc concept[J].Aerospace Science and Technology,2004(8):377-389.

[37] 李錦，江定武，毛枚良,等.BGK-NS格式在復(fù)雜流動(dòng)中的應(yīng)用研究[J].空氣動(dòng)力學(xué)報(bào),2013,31(4):449-453.

[38] Simpson G M,Sadler A J.Controls:A comparison with conventional,planar fins[C]//RTO AVT Symposium.Sorrento,Italy:1998.

[39] James Despirito,Jr Milton E Vaughn,W David Washington.Numerical investigation of canard-controlled missile with planar and grid fins[J].Journal of Spacecraft and Rockets,2003,40(3):363-370.

[40] Montgomery C Hughson,Eric L Blades,Gregg L Abate.Transonic aerodynamic analysis of lattice grid tail fin missiles[C]//24th Applied Aerodynamics Conference.San Francisco,California:2006.

[41] Karl S Orthner.Aerodynamic analysis of lattice grid fins in transonic flow[D].Ohio:Air Force Institute Of Technology,2004.

[42] Gregg Abate,Gerald Winchenbach,Wayne Hathaway.Transonic aerodynamic and scaling issues for lattice fin projectiles tested in a ballistics range[C]//19th International Symposium of Ballistics.Interlaken,Switzerland:2001.

[43] Gregg Abate,Ralf Duckerschein,Wayne Hathaway.Subsonic/transonic free-flight tests of a generic missile with grid fins[C]//38th AIAA Aerospace Sciences Meeting and Exhibit.Reno,NV:2000.

[44] 吳曉軍，陳紅全，鄧有奇,等.柵格繞流數(shù)值模擬研究[J].空氣動(dòng)力學(xué)學(xué)報(bào),2009,27(1):78-82.

[45] Wm David Washington,Pamela F Booth,Mark S Miller.Curvature and leading edge sweep back effects on grid fin aerodynamic characteristics[J].AIAA-93-3480-Cp,1993:659-667.

[46] 譚獻(xiàn)忠，鄧帆，陳少松.翼面氣動(dòng)外形對(duì)柵格翼減阻的影響[J].實(shí)驗(yàn)力學(xué),2013,28(2):255-260.

[47] Daniel Guyot,Erich Schulein.Novel locally swept lattice wings for missile control at high speeds[C]//45th AIAA Aerospace Sciences Meeting and Exhibit.Reno,Nevada:2007.

[48] Marco Debiasi,Zeng Yan,Chng Tat Loon.Swept-back grid fins for transonic drag reduction[C]//28th AIAA Applied Aerodynamics Conference.Chicago,Illinois:2010.

[49] Jinsheng Cai.Numerical study on choked flow over grid-fin configurations[J].Journal of Spacecraft and Rockets,2009,46(5):949-956.

[50] 黎漢華，石玉紅.柵格翼國(guó)內(nèi)外研究現(xiàn)狀及發(fā)展趨勢(shì)[J].導(dǎo)彈與航天運(yùn)載技術(shù),2008(6):27-30.

[51] Jerome P Fanucci,Michael J King,James J Gorman.Method for producing lattice fin for missiles or other fluid-born bodies[P].US006928715B2.Woburn,MA,US:2005.08.16.

[52] Stephen R Perillo,Donald J Atkins.Challenges and emerging trends in store separation engineering-an air force seek eagle office perspective[C]//47th AIAA Aerospace Sciences Meeting Including The New Horizons Forum and Aerospace Exposition.Orlando,Florida:2009.

[53] Paul V Tartabini,Michael G Gilbert,James R Beaty.Flight performance feasibility studies for the max launch abort system[C]//AIAA Atmospheric Flight Mechanics (AFM)Conference.Boston,MA:2013.

(編輯：呂耀輝)

Review of aerodynamic characteristics and application of grid fin

PENG Ke1,HU Fan1,ZHANG Wei-hua1,ZHOU Zhang2

(1.College of Aerospace Science and Engineering, National University of Defense Technology, Changsha 410073, China;2.Jiangnan Industries Group CO.,LTD,Changsha 410200,China)

Researd methods on aerodynamic characteristics of grid fin,were summarized in fields of experimental research,engineering prediction,and numerical simulation.Major aerodynamic characteristics of grid fin,latest geometrical configurations research progresses about curved grid fin and global&local swept grid fin,as well as application trends at home and abroad were discussed.Critical concerns about grid fin were analyzed in view of the bottleneck problems of design and application.

grid fin;aerodynamic characteristics;geometrical configurations;application trends

2014-06-20；

：2014-10-20。

國(guó)家自然科學(xué)基金(51105368)；國(guó)防科技大學(xué)優(yōu)秀研究生創(chuàng)新資助項(xiàng)目(S130105)。

彭科(1989—)，男，博士生，研究方向?yàn)轱w行器總體設(shè)計(jì)與氣動(dòng)外形設(shè)計(jì)優(yōu)化。E-mail：pengke_pk@163.com

V412

A

1006-2793(2015)04-0458-07

10.7673/j.issn.1006-2793.2015.04.002