李強，操小龍

(中國航天科工飛航技術研究院 北京機電工程研究所， 北京 100074)

0 引 言

對于在大氣層內(nèi)的高速飛行器，沖壓發(fā)動機具有比沖高、推重比大、可長時間工作的特點。其進氣道通過激波系的壓縮使來流減速增壓，以滿足發(fā)動機的使用要求，是沖壓發(fā)動機的關鍵部件[1]。低空高速飛行時，該類型飛行器的進氣道壓力可達來流壓力的數(shù)十倍，對結構安全產(chǎn)生較大威脅。因此進氣道壓力載荷的預示和限制，成為飛行器研制過程中必須解決的關鍵問題之一。國外亞燃沖壓發(fā)動機PTV的研制中，采用了進氣道壓力限制技術，通過控制燃料供應量以降低由于進氣道超壓導致的結構失效風險[2]。

進氣道壓力載荷預示一般使用數(shù)值仿真和風洞試驗方法。白曉征等[3]對馬赫數(shù)3.5飛行條件下整流罩開啟過程中的進氣道非定常流動進行了數(shù)值模擬，壓力振蕩峰值達到了來流壓力的51.6倍；熊志平等[4]利用有限體積法研究了超聲速軸對稱進氣道流場，其中馬赫數(shù)3.0、進氣道超臨界狀態(tài)的擴張段壓力峰值為來流的17.5倍；張曉旻等[5]開展了進氣道堵蓋打開過程非穩(wěn)態(tài)流場研究，其中馬赫數(shù)2.5情況下出口堵蓋處的壓力峰值為來流壓力的20倍以上；衛(wèi)峰等[6]研究了曲外錐乘波前體進氣道在超聲速風洞中的壓力特性，其中馬赫數(shù)4.0下的進氣道抗反壓能力為來流壓力的33倍。與進氣道壓力經(jīng)驗公式研究相關的文獻鮮有發(fā)表，僅有少數(shù)工程設計人員從事該領域研究。蒲曉航等[7]針對彈用ATR 進氣道，建立了進氣道超聲速性能快速預估模型，絕大多數(shù)來流及背壓組合工況下的總壓恢復系數(shù)計算誤差不超過10%。

本文以典型頜下進氣高速飛行器為研究對象，使用Fluent軟件對其進氣道流場進行數(shù)值仿真，研究不同馬赫數(shù)、攻角、側滑角及余氣系數(shù)條件下的進氣道壓力特性；針對進氣道載荷快速估算及限制設計需要，以馬赫數(shù)和余氣系數(shù)為主要變量，采用最小二乘法擬合進氣道壓力載荷經(jīng)驗公式；使用經(jīng)驗公式反算真實飛行狀態(tài)下的進氣道壓力，并與飛行試驗測壓結果進行對比，對經(jīng)驗公式的精度進行驗證。

1 進氣道壓力特性計算

1.1 計算模型

本文以美國ASALM飛行器進氣道為研究對象，如圖1所示。

圖1 ASALM飛行器進氣道構型[8]

該進氣道為混壓式超聲速進氣道，采用兩級錐壓縮，進氣道迎風面捕獲面為扇形，內(nèi)通道為扇環(huán)形通道，在尾部轉(zhuǎn)彎導流到飛行器中部變?yōu)閳A形。該進氣道具有良好的攻角特性，能夠適應馬赫數(shù)2.5～4.5工作范圍要求[9-10]。

1.2 數(shù)值模擬方法

采用Fluent軟件，利用有限體積法求解定常三維Navier-Stokes方程，湍流模型為k-ωSST模型，采用隱式離散方法進行時間推進。計算網(wǎng)格為六面體結構化網(wǎng)格，在進氣道壁面附近適當加密。進氣道出口設置為壓力出口邊界。通過風洞試驗和發(fā)動機試車射流試驗測量得到。壁面為絕熱、無滑移、固體邊界。流場計算時先選擇一階精度，待一階精度計算收斂，正激波在進氣道內(nèi)穩(wěn)定下來后，再選擇二階精度進行計算[11]。

計算參數(shù)為：來流馬赫數(shù)2.5～3.5，攻角-3°～+3°，側滑角0°～3°，余氣系數(shù)(A)1.2～5.0(對應進氣道出口反壓/靜壓比范圍為4.8～33.0)。需要說明的是，余氣系數(shù)是實際空氣流量與供入燃料完全燃燒所需理論空氣流量的比值，一般為大于1.0的正數(shù)，其值越小表明燃料供應越多，燃燒室內(nèi)壓力越大[12]。本文中，變余氣系數(shù)是通過改變出口壓力邊界條件實現(xiàn)的，具體方法為：首先在風洞試驗中測量進氣道捕獲流量，然后在發(fā)動機試車射流試驗中根據(jù)已知流量，通過改變供油量實現(xiàn)不同余氣系數(shù)條件，同時測量發(fā)動機熱態(tài)工作下的進氣道出口壓力，作為數(shù)值模擬的輸入。

1.3 計算結果

余氣系數(shù)1.2條件下不同馬赫數(shù)的進氣道壁面壓比如圖2所示，截取進氣道壓縮錐到S型擴張段的上母線部分，Pdl為進氣道壓力與來流壓力之比?？梢钥闯觯航?jīng)過兩級錐面壓縮，進氣道內(nèi)壓力逐漸升高，到擴張段前基本趨于平穩(wěn)；同等余氣系數(shù)條件下，隨著來流馬赫數(shù)增大，前體錐面激波對氣流的壓縮能力不斷提高，進氣道捕獲流量逐漸增大，進氣道壓力也隨之升高；馬赫數(shù)由2.5增大到3.5，進氣道擴張段壓比由12增大到33，增大約175%。

圖2 不同馬赫數(shù)來流的進氣道壁面壓比

馬赫數(shù)2.5條件下、不同余氣系數(shù)的進氣道壁面壓比如圖3所示。

圖3 進氣道壁面壓比與余氣系數(shù)的關系

從圖3可以看出：在同等來流馬赫數(shù)條件下，進氣道捕獲流量不變；但隨著余氣系數(shù)增大，發(fā)動機燃料供應量減小，燃燒室壓力降低，進氣道壓力也隨之降低，結尾激波不斷后移(如圖4所示)，這與文獻[13]中的激波結構變化趨勢類似。余氣系數(shù)由1.2增大到3.5，進氣道擴張段壓比從12減小到5左右，降低約58%。在發(fā)動機正常包線內(nèi)，進氣道工作在超臨界狀態(tài)，燃燒室壓力僅影響結尾激波后的流場，進氣道喉道前的壓力不受余氣系數(shù)影響。這對簡化進氣道載荷的分析是有利的。

圖4 不同余氣系數(shù)狀態(tài)的進氣道激波結構

馬赫數(shù)分別為2.5和3.0，余氣系數(shù)1.2，不同攻角和側滑角條件下的進氣道壁面壓比如圖5所示。

(a) Ma=2.5

(b) Ma=3.0

從圖5可以看出：對于頜下進氣道，攻角由負變正，進氣道捕獲空氣流量增大，因此在同等余氣系數(shù)條件下，進氣道壓力略有升高，其中攻角+3°狀態(tài)的壓比較攻角-3°狀態(tài)大約5%；側滑角對進氣道壓力影響較小，側滑3°狀態(tài)的進氣道壓力較側滑角0°狀態(tài)小約1%?？傮w來看，正常工作包線內(nèi)，較小的攻角、側滑角對進氣道壓力影響不明顯。

2 進氣道壓力經(jīng)驗公式研究

在飛行器進氣道結構設計時，考慮各種拉偏情況的彈道點往往數(shù)以萬計，都開展數(shù)值計算難以實現(xiàn)，只能從中篩選嚴酷情況進行載荷分析。通過數(shù)值仿真和風洞試驗數(shù)據(jù)的歸納總結，將進氣道壓力擬合為若干變量的函數(shù)，不僅方便載荷嚴酷情況的篩選，也能滿足工程上快速估算的需要。另外，在進氣道載荷較大、需要進行限制時，函數(shù)形式的進氣道壓力經(jīng)驗公式也更方便彈道的設計。因此進氣道壓力經(jīng)驗公式研究在型號研制過程中有較迫切的需求。

基于數(shù)值仿真掌握的進氣道壓力載荷規(guī)律特性，開展經(jīng)驗公式的擬合研究。研究時主要采用以下思路：

(1) 進氣道結構主要薄弱點在亞聲速擴張段，該處載荷不僅與飛行馬赫數(shù)、高度、燃料供應量有關，而且受任務需求限制，彈道的馬赫數(shù)-高度剖面難以作大的調(diào)整；可供控制的主要是燃料供應量，可使用余氣系數(shù)進行定量描述。

(2) 與零攻角、零側滑角狀態(tài)相比，攻角±3°狀態(tài)的進氣道擴張段壓力變化小于3%，側滑角3°的變化小于1%。在小姿態(tài)角范圍內(nèi)，可忽略攻角、側滑角對進氣道壓力的影響。

(3) 將有量綱的進氣道擴張段壓力分解為來流壓力和無量綱進氣道壓比的乘積。前者主要與飛行高度有關，后者與飛行馬赫數(shù)和余氣系數(shù)有關。進氣道壓力可描述為

Pinlet(Ma,H,A)=P∞(H)×Pdl(Ma,A)

(1)

式中：Pinlet、P∞、Pdl和A分別為進氣道壓力、來流壓力、無量綱進氣道壓比和余氣系數(shù)。P∞可根據(jù)國軍標[14]或文獻[15]中的公式確定。

(2)

臨界壓比Pdl(A=1.2)與馬赫數(shù)有關，使用最小二乘法擬合為馬赫數(shù)的線性函數(shù)[16]，如式(3)所示。式(3)與仿真數(shù)據(jù)的對比如圖6所示。

Pdl(A=1.2)=-15.25+9.943Ma

(3)

圖6 式(3)與仿真數(shù)據(jù)的對比

(4)

圖7 式(4)與仿真數(shù)據(jù)的對比

使用式(2)～式(4)計算馬赫數(shù)2.5～3.0、余氣系數(shù)1.2～5.0狀態(tài)的進氣道擴張段壓比，并與數(shù)值仿真的對比如圖8所示。

圖8 進氣道壓比經(jīng)驗公式與仿真值的對比

從圖8可以看出：對于圖8中的11個狀態(tài)點，除1個點擬合誤差較大外，其他點擬合誤差均小于5%，平均擬合誤差為2.46%?？紤]到經(jīng)驗公式主要用于進氣道壓力快速評估和彈道控制律設計，因此該精度能夠滿足型號工程研制階段的載荷預示要求。

3 飛行試驗驗證

該飛行器兩次飛行試驗中，在發(fā)動機工作窗口對進氣道擴張段壁面壓力進行測量，采集頻率為50 Hz，采集時長分別為59.8和5.2 s。試驗后利用飛行試驗遙測的靜壓、馬赫數(shù)、余氣系數(shù)，使用式(1)～式(4)對進氣道擴張段壓力進行反算，其與飛行試驗測量值的對比如圖9所示。

(a) 飛行試驗1

(b) 飛行試驗2

從圖9可以看出：在測試窗口內(nèi)，計算的進氣道載荷變化趨勢和遙測值符合較好，有效測試時間內(nèi)大部分計算值與飛行試驗測量值相差小于5%，載荷峰值最大相差13.6%。經(jīng)分析，誤差一方面與經(jīng)驗公式本身的擬合誤差有關，另一方面也與飛行試驗遙測誤差及飛行器姿態(tài)的快速變化有關?？紤]到經(jīng)驗公式主要用于設計階段的載荷情況篩選和快速估算，該精度已經(jīng)能夠滿足工程應用需要。

4 結 論

(1) 超聲速進氣道壓比隨馬赫數(shù)增大而增大，隨余氣系數(shù)增大而減?。惠^小的攻角、側滑角對進氣道壓力的影響不明顯。

(2) 進氣道壓力主要受馬赫數(shù)、高度和余氣系數(shù)影響，基于無量綱化、解耦方法擬合的經(jīng)驗公式具有較高的精度，可用于載荷快速評估和彈道控制律設計。