戰(zhàn)斗機大迎角/過失速機動下的進氣道氣動特性

向歡，楊應(yīng)凱，謝錦睿，吳永勝

中國航空工業(yè)成都飛機設(shè)計研究所，成都 610091

超機動性是第4代戰(zhàn)斗機的基本特征之一，而過失速機動是戰(zhàn)斗機具有超強機動能力的重要標志。過失速機動能力的具備可使戰(zhàn)斗機機頭指向速度大幅提高，近距格斗能力顯著改善，同時可突破傳統(tǒng)飛行禁區(qū)，最大限度地挖掘戰(zhàn)機潛能。要實現(xiàn)過失速機動和全包線無憂慮飛行能力，優(yōu)異的進發(fā)匹配特性至關(guān)重要，因為大迎角和過失速機動狀態(tài)下進氣畸變特性急劇惡化，帶來巨大的發(fā)動機穩(wěn)定裕度損失，極易引起發(fā)動機喘振停車，給飛行安全帶來嚴重威脅。因此戰(zhàn)斗機要想獲得過失速機動能力，必須摸清大迎角及過失速機動飛行狀態(tài)下進氣道的氣動特性。

常見的過失速機動包括眼鏡蛇機動、赫伯斯特機動、榔頭機動、輪盤機動等多種，具有飛行迎角大、速度低、三軸角速率大以及速度/方向/高度大幅變化等特點，同時進氣道內(nèi)部流動存在強烈的非定常和非線性特征，因此給理論研究、計算仿真、試驗技術(shù)、驗證手段等方面都提出了巨大挑戰(zhàn)。

20世紀90年代，美國先后通過F-15 S/MTD、F-16 MATV、F-15 ACTIVE、F/A-18A HARV等計劃開展了大量的大迎角進氣道特性研究工作。其中在F/A-18A HARV項目中開展了大迎角進發(fā)相容性飛行試驗，主要研究了大迎角穩(wěn)定飛行時進氣道氣動特性[1]、快速機動等動態(tài)飛行過程與穩(wěn)定飛行條件下進發(fā)匹配特性的差異[2-4]、大迎角大偏航速率快速脫離機動飛行過程中進氣道流動特性及發(fā)動機失速原因[5]、CFD計算仿真對大迎角進氣道特性的預(yù)測能力[6-11]、極端復(fù)雜工況下發(fā)動機換算流量的測量方法[12-13]等問題，取得了豐碩的成果。研究發(fā)現(xiàn)，穩(wěn)定飛行時隨著迎角和側(cè)滑角增大，進氣道總壓恢復(fù)系數(shù)降低、峰值周向畸變和紊流度均增大，且迎角30°以上時此影響顯著加劇[1]?？焖俑┭鰴C動飛行時，進氣道總壓恢復(fù)系數(shù)與對應(yīng)穩(wěn)定飛行狀態(tài)差異不大，而周向畸變指數(shù)則普遍低于穩(wěn)定飛行狀態(tài)的[3-4]。而大迎角下大偏航速率快速脫離機動飛行時出現(xiàn)了12次發(fā)動機失速現(xiàn)象[14]，進一步研究表明除了進氣畸變指數(shù)、壓氣機穩(wěn)定裕度損失以外，飛機三軸運動角速率是發(fā)動機失速的重要影響因素[15]。

而中國在大迎角/過失速機動下氣動特性研究方面起步較晚。近年來，中國空氣動力研究與發(fā)展中心巫朝君等[16]采用雙力矩電機同步驅(qū)動研制了快速俯仰機動進氣道試驗裝置，建立了戰(zhàn)斗機進氣道非定常性能試驗方法。成都飛機設(shè)計研究所楊應(yīng)凱[17]對某型飛機進氣道在快速俯仰機動條件下的動態(tài)特性進行了試驗研究，獲得了一定迎角變化范圍內(nèi)進氣道氣動特性，發(fā)現(xiàn)進氣道總壓恢復(fù)系數(shù)和畸變指數(shù)都出現(xiàn)了很強的非定常遲滯效應(yīng)，與相同迎角下的穩(wěn)態(tài)結(jié)果存在明顯差異。在計算仿真方面，南京航空航天大學胡鈴心等[18]采用基于非結(jié)構(gòu)嵌套網(wǎng)格的算法對一種鴨式布局戰(zhàn)斗機眼鏡蛇機動下的外部氣動力特性進行了研究，分析了非定常流動機理；而成都飛機設(shè)計研究所向歡[19-21]則應(yīng)用動態(tài)嵌套網(wǎng)格對戰(zhàn)斗機在快速俯仰機動下的進氣道動態(tài)氣動特性進行了研究，同樣捕捉了非定常遲滯效應(yīng)，并分析了其影響因素與機理。但目前國內(nèi)開展的研究都僅限于地面仿真與試驗，而大迎角/過失速機動飛行試驗尚處空白，天地相關(guān)性問題的研究也無從談起。

因此，有必要通過大迎角風洞試驗和計算仿真相結(jié)合的地面驗證以及飛行試驗等手段，系統(tǒng)完備地研究進氣道在大迎角穩(wěn)態(tài)/過失速機動下的氣動特性及其天地相關(guān)性，解決過失速機動下的非定常非線性研究難題。

1 大迎角穩(wěn)態(tài)特性

1.1 計算仿真和風洞試驗

以某戰(zhàn)斗機為研究對象，首先通過計算仿真和風洞試驗相結(jié)合的地面驗證手段研究進氣道在大迎角穩(wěn)定飛行狀態(tài)下的氣動特性。

計算仿真采用基于有限體積法求解三維可壓縮Navier-Stokes方程的InHouse計算程序，其中無黏對流通量的空間離散采用具有較高間斷分辨率的二階迎風Roe格式，黏性擴散通量的離散則采用二階中心差分格式，湍流模型采用Realizablek-ε模型，近壁面采用標準壁面函數(shù)，時間推進采用定常隱式算法。通過ICEM商業(yè)軟件生成全結(jié)構(gòu)化網(wǎng)格(Structural Gird)以提高計算精度，流場計算域網(wǎng)格結(jié)點數(shù)約2 000萬，壁面第1層 網(wǎng)格法向高度取1.0×10-5m。

采用全金屬縮尺模型開展大迎角穩(wěn)態(tài)風洞試驗，試驗方法和測試手段不再贅述，試驗狀態(tài)面向過失速機動任務(wù)。如圖1所示，根據(jù)大迎角滾轉(zhuǎn)、赫伯斯特、榔頭、眼鏡蛇、直升機等5個過失速機動動作在不同進入條件下的控制律仿真結(jié)果，確保馬赫數(shù)Ma、迎角α、側(cè)滑角β等進氣道風洞試驗條件范圍對飛行試驗的全覆蓋。

運用上述計算仿真和風洞試驗方法，建立了包含高度、馬赫數(shù)、迎角、側(cè)滑角、發(fā)動機狀態(tài)等影響因素在內(nèi)的全包線進氣道氣動特性數(shù)據(jù)庫。圖2為不同迎角下進氣道穩(wěn)態(tài)氣動性能，圖中φ為流量系數(shù)，σ為總壓恢復(fù)系數(shù)?？梢钥吹?，迎角為0°時進氣道總壓恢復(fù)性能最高，隨著迎角的不斷增大，特別是大于30°以后，進氣道唇口分離愈發(fā)嚴重，導(dǎo)致進氣道性能顯著降低，同時不同迎角下計算仿真與風洞試驗結(jié)果吻合良好。

圖1 典型過失速機動仿真曲線

圖2 進氣道大迎角穩(wěn)態(tài)氣動性能

為了更好地揭示進氣道內(nèi)部流動形成發(fā)展機理，圖3列舉了不同迎角下進氣道內(nèi)部主要流動特征?？梢钥吹剑∮菚r進氣道出口圖譜中存在3個低壓區(qū)，分別位于左上、右上和正下方3個方位，其中左上、右上2個低壓區(qū)是由氣流繞過機頭兩側(cè)進入進氣口發(fā)展成S彎管道對轉(zhuǎn)旋渦形成的，而唇口繞流進入內(nèi)管道后一直沿著壁面發(fā)展最終附面層越積越厚形成正下方低壓區(qū)。隨著迎角增大，氣流不再繞機頭流動而是直接沖入進氣口發(fā)展成S彎管道對轉(zhuǎn)旋渦，而下唇口繞流的一部分分支沿側(cè)壁面向上與對渦匯合，另一分支則繼續(xù)沿下壁面發(fā)展，因此左上、右上2個低壓區(qū)勢力范圍加大，并有向下運動的趨勢。當迎角繼續(xù)增大，對渦氣流在管道后段下甩，形成與唇口繞流進一步合攏的發(fā)展態(tài)勢，因此出口截面上對渦低壓區(qū)也與唇口繞流低壓區(qū)匯合。此后，迎角進一步增大，對渦氣流與唇口繞流的匯合作用更加強烈，進氣道出口最終呈現(xiàn)出高壓區(qū)在上、低壓區(qū)在下的圖譜形態(tài)。

圖3 不同迎角下進氣道流動特性

1.2 飛行試驗

通過在發(fā)動機進口截面處加裝進氣道測量耙開展大迎角穩(wěn)態(tài)進發(fā)匹配特性的飛行試驗研究，測量系統(tǒng)由穩(wěn)態(tài)和動態(tài)2部分組成，與風洞試驗一致。此處研究的大迎角穩(wěn)態(tài)，是指飛行高度、馬赫數(shù)、迎角、側(cè)滑角、發(fā)動機油門桿等主要影響參數(shù)均穩(wěn)定保持1 s時間以上的飛行狀態(tài)。

圖4為不同迎角下飛行試驗所得進氣道出口總壓恢復(fù)系數(shù)、綜合畸變指數(shù)(W)與風洞試驗結(jié)果的對比。可以看到，與風洞試驗趨勢一致，隨著迎角的增大進氣道氣動性能持續(xù)下降，表現(xiàn)為總壓恢復(fù)的降低與綜合畸變指數(shù)的增大?？傮w來講，風洞試驗和飛行試驗的各參數(shù)天地一致性良好?？倝夯謴?fù)系數(shù)在中等迎角處存在最大約0.02的差異，而小迎角和大迎角下差異均小于0.01，風洞試驗結(jié)果線性度高于飛行試驗；綜合畸變指數(shù)在小迎角范圍內(nèi)存在小于1%的差異，而大迎角時吻合度高。

圖5為某大迎角下飛行試驗測得進氣道出口截面總壓分布云圖與風洞試驗結(jié)果的對比，可見二者高低壓區(qū)分布一致，均呈現(xiàn)出高壓區(qū)在上、低壓區(qū)在下的總體形態(tài)，與圖3進氣道流動發(fā)展機理分析結(jié)論相吻合。

因此，不管是總體氣動特征參數(shù)還是內(nèi)管道流動形態(tài)，天地相關(guān)性均吻合良好。表明本文采用的研究手段和方法是合理正確的，可有效支撐大迎角進氣道氣動特性的研究。

圖4 進氣道大迎角穩(wěn)態(tài)氣動性能飛行試驗與風洞試驗的結(jié)果對比

圖5 進氣道大迎角總壓分布云圖飛行試驗與風洞試驗的結(jié)果對比

2 過失速機動非定常非線性特性

2.1 計算仿真

為掌握過失速機動下進氣道非定常非線性氣動特性與規(guī)律，揭示進氣道動態(tài)特性機理，采用基于動態(tài)嵌套網(wǎng)格的URANS算法進行仿真。計算程序中時間推進采用雙時間步非定常隱式算法，物理時間步長Δt取5.0×10-6s。

建立圖6所示包含2個同心球體的計算域，分別為飛行空域的外部流場和包含飛機在內(nèi)的內(nèi)部流場，采用ICEM軟件生成全結(jié)構(gòu)嵌套網(wǎng)格，內(nèi)外流場邊界相互交叉。在模擬過失速機動的計算過程中外部流場域靜止不動，而包含飛機機體在內(nèi)的內(nèi)部流場網(wǎng)格隨著飛機一起繞重心作三軸旋轉(zhuǎn)運動，飛機每運動一步都需要搜尋內(nèi)、外流場間的嵌套邊界和宿主單元，并交換流場信息。

為簡化問題只考慮戰(zhàn)斗機快速俯仰機動過程，即只通過迎角變化近似模擬眼鏡蛇機動。給定內(nèi)部流場網(wǎng)格隨著飛機一起繞重心作單軸旋轉(zhuǎn)運動的規(guī)律為正弦函數(shù)，同時研究俯仰速率對進氣道動態(tài)特性的影響，給定圖7所示的4種運動規(guī)律。需要指出的是，因正弦函數(shù)的特點，運動過程中不同迎角下的角速度(即俯仰速率)不同，表現(xiàn)為中間高、兩端低，因此圖中F1、F2、F3、F4均為平均俯仰速率。

計算得到的進氣道非定常非線性特性如圖8所示。從圖8(a)中可以看到，進氣道總壓恢復(fù)動態(tài)特性總體上隨迎角增大而降低，與穩(wěn)態(tài)特性一致，但在俯仰機動周期內(nèi)形成一個“∞”形的交叉閉合環(huán)路。正迎角部分，上仰過程性能高于下俯過程，呈現(xiàn)出明顯的遲滯現(xiàn)象，且20°～80°迎角范圍內(nèi)差異最大，而最大迎角附近二者差異不大，表明瞬時俯仰速率越大遲滯現(xiàn)象越嚴重。這是由于飛機姿態(tài)改變所引起的進氣道內(nèi)部流場重新建立需要一定時間，進氣道出口流場的瞬時反應(yīng)始終滯后于飛機姿態(tài)的實時變化。而飛機俯仰速率越大，進氣道出口流場反應(yīng)相對于飛機姿態(tài)變化的滯后作用也就越大。此現(xiàn)象在圖8(b)中得到印證，平均俯仰速率越大，“∞”形交叉閉合環(huán)路呈現(xiàn)出向外擴張一環(huán)套一環(huán)、愈加飽滿的趨勢。此外，整個機動歷程中進氣道性能最差狀態(tài)點出現(xiàn)在迎角最大處。

圖8 過失速機動計算非定常非線性特性

2.2 飛行試驗

同樣通過在發(fā)動機進口截面處加裝進氣道測量耙的方法，開展了多個過失速機動動作下進發(fā)匹配動態(tài)特性的飛行試驗。

選取典型的眼鏡蛇機動進行分析，其迎角飛行參數(shù)如圖9所示?？梢钥吹?，進入機動動作后飛機快速拉桿增大迎角直至到達最大迎角，然后低頭改出。總體上呈現(xiàn)出非線性特征，與計算仿真給定理想的正弦函數(shù)有所區(qū)別。

圖10為眼鏡蛇機動下進氣道氣動特性的動態(tài)響應(yīng)。可以看到，進入機動后隨著飛機快速拉桿到最大迎角，進氣道總壓恢復(fù)系數(shù)顯著降低、綜合畸變指數(shù)急劇增大；從最大迎角到改出機動一段，總壓恢復(fù)系數(shù)逐漸增大、綜合畸變指數(shù)則明顯回落；整個機動歷程的進氣道特征參數(shù)曲線仍然是非線性的。

圖9 眼鏡蛇機動飛行參數(shù)

圖10 眼鏡蛇機動下進氣道特性動態(tài)響應(yīng)

從圖11進氣道動態(tài)特性隨迎角的變化規(guī)律可以看到，眼鏡蛇機動飛行周期內(nèi)總壓恢復(fù)系數(shù)和綜合畸變指數(shù)特性均表現(xiàn)為一個閉合的遲滯環(huán)路，上仰過程的特征參數(shù)均高于下俯過程，呈現(xiàn)出明顯的遲滯現(xiàn)象，而進氣道性能最差狀態(tài)點同樣出現(xiàn)在最大迎角處，與計算仿真結(jié)論一致。

因此，過失速機動下進氣道動態(tài)特性的天地相關(guān)性同樣吻合良好，表明本文采用的研究方法可有效解決過失速機動下進氣道流動非定常非線性的研究難題。

圖11 眼鏡蛇機動下進氣道非定常非線性特性

3 結(jié) 論

本文通過計算仿真和風洞試驗相結(jié)合的地面驗證以及飛行試驗等手段，系統(tǒng)研究了戰(zhàn)斗機在大迎角和過失速機動下的進氣道氣動特性。

1) 大迎角穩(wěn)定飛行下進氣道性能的天地相關(guān)性吻合良好，眼鏡蛇機動下進氣道非定常氣動特性呈現(xiàn)出明顯的遲滯效應(yīng)。

2) 獲得了戰(zhàn)斗機在大迎角和過失速機動下的進氣道氣動特性，建立了過失速機動下進氣道非定常非線性特性問題的研究方法。