崔 凱徐應(yīng)洲 肖 堯 李廣利

(中國科學(xué)院力學(xué)研究所高溫氣體動力學(xué)國家重點實驗室，北京100190)

(中國科學(xué)院大學(xué)工程科學(xué)學(xué)院,北京100049)

乘波體壓縮面變化對其氣動性能影響分析1)

崔 凱2)徐應(yīng)洲 肖 堯 李廣利

(中國科學(xué)院力學(xué)研究所高溫氣體動力學(xué)國家重點實驗室，北京100190)

(中國科學(xué)院大學(xué)工程科學(xué)學(xué)院,北京100049)

乘波體是一種利用激波包裹特性獲得高升阻比的高速飛行器構(gòu)型.已有研究中，乘波體氣動性能的改善主要依賴于給定源流場條件下的前緣型線優(yōu)化.本文采用數(shù)值優(yōu)化和計算流體力學(xué)模擬為主要手段分析了乘波體壓縮面變化對其氣動性能的影響，以期有效拓展乘波體的設(shè)計空間.主要內(nèi)容如下：首先給出了一種基于表面局部變形的乘波體設(shè)計方法.其次結(jié)合運用增量修正參數(shù)化方法、計算流體力學(xué)分析和微分演化算法構(gòu)造了乘波體壓縮面外形氣動優(yōu)化設(shè)計流程，以一種橢圓錐形流場生成的乘波體作為基準構(gòu)型開展了無黏優(yōu)化.之后從優(yōu)化結(jié)果中選擇升阻比遞增的6個典型構(gòu)型進行前緣鈍化處理后，基于N-S方程對其氣動性能進行了評估.最后綜合依據(jù)無黏/黏性計算結(jié)果分析了乘波體壓縮面變化對其氣動性能的影響.結(jié)果表明該部分形狀的改變對乘波體氣動性能影響十分明顯，在升力面積不變的條件下，乘波體壓縮面形狀變化可導(dǎo)致其升阻比出現(xiàn)成倍變化，即使在升力不減條件下，升阻比較基準構(gòu)型也可獲得超過14%的提升.此外，還可導(dǎo)致乘波體相對壓心系數(shù)出現(xiàn)明顯偏移.

乘波體，優(yōu)化設(shè)計，增量參數(shù)化方法，微分演化算法，計算流體力學(xué)

引言

目前高超聲速飛行器的氣動布局主要包括升力體、軸對稱旋成體、翼身融合體和乘波體等幾類[1]，其中乘波體[2-4]是一種較為獨特的構(gòu)型概念.該構(gòu)型基于給定流場，通過激波面切割和流線追蹤方式設(shè)計.在設(shè)計飛行條件下(給定馬赫數(shù)、攻角等)，高速飛行時產(chǎn)生的弓形激波完全附著于飛行器的邊緣，上下表面沒有流動泄露，激波后高壓區(qū)可完全包裹于飛行器的下半部分，使飛行器升阻比有效提升.因飛行時好像乘在激波之上，故稱為“乘波體”.近年來，隨著高超聲速飛行器相關(guān)研究的不斷升溫，有關(guān)乘波體的各類研究也逐漸成為熱點之一.

由于乘波體可基于任意超聲速流場生成，因而源自不同流場的乘波體設(shè)計方法和性能分析是主要研究內(nèi)容之一.早期工作主要圍繞一些相對簡單的源流場，如正圓錐流場[5]、傾斜圓錐、橢圓錐及帶有縱向曲率的旋成體流場[6-8]、楔錐流場等[9].近年來，隨著數(shù)值模擬技術(shù)的不斷發(fā)展，一些更加復(fù)雜的流場，如任意截面錐形體流場[10-11]、卡門曲線為母線的旋成體流場[12]和多級壓縮錐體流場[13]等亦被用于生成乘波體.在密切錐類方法基礎(chǔ)上，也進一步發(fā)展了密切流場法[14]和密切內(nèi)錐設(shè)計方法等[15-16].此外，近期還有直接對錐體及其流場進行切割的被動乘波體設(shè)計方法[17-18].上述研究均在不同程度上改善了乘波構(gòu)型的綜合性能.在此基礎(chǔ)上，以乘波體為基礎(chǔ)的創(chuàng)新構(gòu)型和高超聲速飛行器新概念也不斷涌現(xiàn)[19-25]，這意味著乘波體已逐漸開始從單純的概念研究逐漸向?qū)嵱没A段推進.

伴隨乘波體相關(guān)研究的深入，氣動外形優(yōu)化也逐漸成為此領(lǐng)域重要研究內(nèi)容.這一領(lǐng)域的研究始于20世紀60年代.最具代表性的工作為 Bowcutt等[26]和Corda等[27]于1987年和1988年分別針對圓錐流場和冪次率母線軸對稱體流場生成的乘波體，采用參考溫度法計算摩擦阻力，并以非線性單純形算法作為優(yōu)化驅(qū)動器對前緣線進行優(yōu)化，獲得了一系列“黏性優(yōu)化乘波體”.基于類似的優(yōu)化流程，后續(xù)工作通過結(jié)合使用不同的乘波體設(shè)計方法及不同優(yōu)化方法，特別是現(xiàn)代啟發(fā)式算法，將這一研究由單目標優(yōu)化逐漸拓展至多目標優(yōu)化設(shè)計[19-20,28-32].

依據(jù)其設(shè)計流程，乘波體由其源流場和位于激波面上的前緣線共同確定.從目前公開發(fā)表文獻看，乘波體的外形優(yōu)化均在給定源流場條件下，通過改變前緣線形狀進行優(yōu)化設(shè)計，對乘波體下表面(即壓縮面)形狀的研究相對較少.由于乘波體壓縮面為迎流面，是乘波體升力和阻力的主要來源，因此其外形變化必然會導(dǎo)致乘波體氣動性能的大幅改變.文獻[33-34]以黏性優(yōu)化乘波體為基礎(chǔ)，分別采用風(fēng)洞實驗和計算流體力學(xué)方法對兩個前緣線俯視形狀相同，但壓縮面具有明顯差別的乘波體進行了氣動性能分析，結(jié)果表明兩種構(gòu)型的升阻比差別超過5%，該文結(jié)果表明乘波體壓縮面的形狀變化對其氣動性能的影響不可忽略.因此，對壓縮面形狀的合理設(shè)計以及優(yōu)化問題進行分析和研究，對拓展乘波體及乘波飛行器的設(shè)計空間，獲得氣動性能更為優(yōu)良的乘波體無疑具有積極意義.但從研究現(xiàn)狀看，這一方向的工作尚未看到公開文獻發(fā)表.

本文主要基于上述背景，從應(yīng)用基礎(chǔ)研究角度出發(fā)，以橢圓錐導(dǎo)乘波體為實例，開展了乘波體壓縮面變化對其氣動性能的影響研究.具體工作內(nèi)容簡述如下：首先對給定前緣線條件下(等價于升力面積不變)乘波體壓縮面合理設(shè)計問題進行了分析，在此基礎(chǔ)上基于超聲速流場擾動僅向下游傳播的特點，給出了一種基于表面局部變形的乘波體構(gòu)型設(shè)計方法.之后以無黏優(yōu)化設(shè)計和黏性分析作為手段，分析了乘波體壓縮面形狀變化對其氣動性能的影響.結(jié)果表明，即使在考慮黏性及鈍化條件下，乘波體壓縮面的變化也可導(dǎo)致其升力、阻力、升阻比和壓心等參數(shù)出現(xiàn)大幅改變.

1 乘波體壓縮面設(shè)計問題分析

1.1 源自給定流場的乘波體方法簡介

源于給定流場的乘波體設(shè)計方法在相關(guān)文獻中已有較多介紹[8-13]，此處僅進行簡單闡述.乘波體生成流程如圖1所示.首先，給定一個任意三維幾何體作為生成體，在設(shè)計條件下，應(yīng)用數(shù)值求解歐拉方程可獲得其無黏流場，同時捕捉到激波面；其次，在激波面上定義一條三維曲線，即乘波體的前緣線；隨后，從前緣線上的每一點出發(fā)，沿來流方向進行流線追蹤，所有流線構(gòu)成的曲面做為乘波體的下表面(壓縮面)；最后，采用自由流面方法設(shè)計乘波體上表面并完成設(shè)計.由于流面上沒有法向流動，因而在設(shè)計條件下，乘波體下表面的繞流流場與源無黏流場一致，同時由于乘波體的前緣線精確位于激波面上，下表面經(jīng)壓縮激波后產(chǎn)生的高壓區(qū)不會泄露到上部，因此可保證其良好的氣動性能.

圖1 乘波體設(shè)計流程圖示Fig.1 Sketch of waverider construction

1.2 給定前緣線條件下乘波體壓縮面設(shè)計問題分析

從前述乘波體的設(shè)計流程中不難看出，乘波體的外形主要由兩個因素確定：(1)生成乘波體的源流場，(2)位于激波面上的前緣線.源流場給定后，激波面可唯一確定，如何在激波面上合理選擇前緣線，即前緣線優(yōu)化問題已有較多研究.然而，源流場與激波面并非一一對應(yīng)，即同一激波面可能對應(yīng)多個源流場，圖2基于二維圖示對此進行說明.當給定某一圓錐及來流條件時，可獲得激波面形狀.假定將圓錐母線從中間某一位置開始進行有限幅度的上彎或下彎，如圖2中虛線所示，由于前半部分仍保持錐體外形，因此激波面仍將與初始圓錐相同，但波后流場不同.此時，當在此激波面上給定相同的前緣線時，基于不同的流場中也會必然獲得下表面(壓縮面)形狀迥異的乘波體構(gòu)型.由于乘波體壓縮面為迎流面，其外形的變化必然導(dǎo)致乘波體的氣動性能產(chǎn)生明顯差異.

圖2 具有相同激波面的源幾何體Fig.2 Generating bodies with the same shock layer

由上可知，乘波體壓縮面的變化等價于改變生成乘波體的源流場.依據(jù)這一思路，雖然可以基于前述乘波體設(shè)計流程來逐次修改其壓縮面外形，但將大幅增加流程的復(fù)雜性，同時，如果考慮源幾何體的三維局部變形，源流場的求解只能依賴于三維計算流體力學(xué)分析，這也將帶來計算消耗的大幅增加.此外，采用這一思路是將乘波體的參數(shù)化間接采用源幾何體的參數(shù)化替代，這也增加了參數(shù)化方法選擇和設(shè)計空間界定的難度.針對上述問題，本文基于超聲速流場的特點，采用表面局部變形方法生成乘波體，并基于增量修正方法進行參數(shù)化.

2 基于表面局部變形的乘波體修正方法和增量參數(shù)化設(shè)計

乘波體的設(shè)計馬赫數(shù)一般在高超聲速區(qū)域，在此條件下，流場中的擾動僅向下游傳播，且擾動范圍限于激波錐之內(nèi).乘波體一般為扁平體，當壓縮面的中后部分區(qū)域產(chǎn)生有限幅度的外形變化時，若能合理限制變形區(qū)域和幅度，其對流場所產(chǎn)生的擾動將不會干擾乘波體的邊緣區(qū)域，因此其乘波特性，即激波附著邊緣及高壓區(qū)包裹特性將不會受到破壞，可生成新的乘波體外形.顯然，采用這樣的乘波體生成方式，可大幅簡化前述分析和優(yōu)化流程.

由于乘波體壓縮面為流面，無法解析表達，若采用樣條曲面對其進行擬合，設(shè)計參數(shù)個數(shù)及各設(shè)計參數(shù)的控制區(qū)域均難以控制.基于上述情況，本文采用增量修正方法[25]進行參數(shù)化設(shè)計.簡而言之，該方法是將外形的參數(shù)化轉(zhuǎn)化為其修正量的參數(shù)化，并在網(wǎng)格坐標系統(tǒng)下將修正量疊加至原有外形，數(shù)學(xué)表達式如下

式中，i和j分別為控制域內(nèi)的面網(wǎng)格點編號為控制域內(nèi)基準外形網(wǎng)格點的三維空間坐標向量，Δri,j為該點坐標修正量向量，ri,j為修正后的該點外形坐標值向量.Ax,Ay和Az分別為3個空間坐標方向的最大增量值，Cx,Cy和Cz分別對應(yīng)3個方向的控制函數(shù)，p為3個方向控制函數(shù)的參變量集，實際中作為設(shè)計變量使用.圖3給出了一個實例，圖中，左側(cè)上下兩圖分別為初始外形和修正量形狀，右側(cè)為疊加后的外形.

圖3 外形參數(shù)化實例Fig.3 An example of shape parameterization

基于本文問題的特點并從縮減計算量的角度考慮，參數(shù)化主要針對乘波體壓縮面的對稱(軸向)和尾緣(展向)剖面的兩條邊界曲線，乘波體前緣線固定不變.沿軸向的控制函數(shù)曲線采用廣義拋物線進行參數(shù)化，其歸一化的表達形式為

其中，P1為直線段與曲線段的交接點，在優(yōu)化設(shè)計中作為設(shè)計變量.沿展向的控制函數(shù)曲線采用廣義余弦函數(shù)進行參數(shù)化，其歸一化的表達形式為

其中，P2和P3為余弦曲線兩端直線段與曲線段的交接點，在優(yōu)化設(shè)計中作為設(shè)計變量.此外，修正量的外形基于邊界線的變化插值獲得.采用上述參數(shù)化方法，通過參數(shù)P1和P3的合理選擇，可保證變形區(qū)域距離邊界保留足夠距離，進而保證乘波體壓縮面邊緣區(qū)域不發(fā)生變形，保持其乘波特性.

3 乘波體壓縮面無黏優(yōu)化和結(jié)果分析

基于優(yōu)化設(shè)計分析乘波體壓縮面形狀變化對其氣動性能的影響，可在設(shè)計空間內(nèi)獲得充足的研究樣本.優(yōu)化設(shè)計基于無黏分析進行，其主要原因如下：(1)依據(jù)設(shè)計原則，乘波體僅在尖前緣及無黏條件下方可完全保證其乘波特性，因而基于無黏計算可嚴格測試前述基于局部表面變形的乘波體生成和參數(shù)化設(shè)計方法；(2)基于無黏分析可大幅縮減計算量及提高優(yōu)化流程的魯棒性.后續(xù)計算中，自由來流條件均給定為飛行馬赫數(shù)6和高度25km.

3.1 基準乘波體構(gòu)型及性能參數(shù)

基準乘波體構(gòu)型采用基于任意錐體流場的乘波體外形設(shè)計軟件包 GCBWRG[10]，源于截面長寬比為1.618的橢圓錐體流場生成.長度給定為2m，其上表面為水平自由來流面.圖4給出了飛行攻角0°條件下乘波體沿軸向不同截面壓力分布云圖.可以看出，高壓區(qū)均完全包裹于乘波體的下方，乘波性能良好.圖5繪出了該乘波體升阻比隨飛行攻角的變化曲線，考慮底部阻力時，0°攻角下的升阻比為5.65，最大升阻比出現(xiàn)于2°攻角狀態(tài)，數(shù)值為6.55，對應(yīng)升力系數(shù)Cl=0.06235(以俯視投影為參考面積，其數(shù)值為1.4492m2).

圖4 基準乘波體(半模)及其截面壓力分布Fig.4 Baseline waverider(half model)and its pressure contours at dif f erent cross section

圖5 基準乘波體升阻比隨攻角變化曲線Fig.5 Variation of the lift-to-drag ratio values with fligh angle of attack for baseline waverider

3.2 優(yōu)化目標和流程

優(yōu)化目標給定為升阻比最大化，同時將升力系數(shù)不減作為不等式約束.優(yōu)化模型的數(shù)學(xué)表達式為

式(4)中，p為設(shè)計變量構(gòu)成的矢量，K為升阻比，Cl為升力系數(shù).將約束條件采用采用懲罰函數(shù)法處理并進行歸一化后，式(4)的約束優(yōu)化問題可轉(zhuǎn)化為如下的無約束優(yōu)化問題

式 (5)中，K0為基準乘波體的最大升阻比 (K0=6.5)，λ為懲罰因子，這里固定取為100.采用歸一化表達的目的為避免升阻比和升力系數(shù)間的數(shù)值不匹配.

優(yōu)化設(shè)計的設(shè)計變量共5個，包括外形最大增量值(控制外形變化量)、式(2)和式 (3)中的參數(shù)P1～P3(控制外形變化區(qū)域和形狀),以及飛行攻角.上述5個參數(shù)及其給定的上下邊界如表1所示.表中，Hmax為最大增量值，α為飛行攻角(度).

表1 設(shè)計變量的上下邊界Table 1 Boundary values of design space

優(yōu)化算法采用微分演化算法[35]，該方法屬于進化類算法的一種，主要思路為從一個初始種群出發(fā)，經(jīng)過變異、交叉、遺傳等操作來獲得目標函數(shù)的最優(yōu)值.該方法與遺傳算法的思路相似，具有不需要計算目標函數(shù)的梯度，魯棒性強，全局最優(yōu)值搜索能力較強等優(yōu)點.同時，由于該方法直接采用實數(shù)編碼，且在染色體的交叉選擇及種群更新等方面的處理優(yōu)于遺傳算法，因此更加高效.優(yōu)化算法的流程圖如圖6所示.

圖6 微分演化算法流程圖Fig.6 Flowchart of dif f erential evolution algorithm

3.3 優(yōu)化結(jié)果及分析

圖7 升阻比(a)和升力系數(shù)(b)收斂圖Fig.7 Convergence history of theL/D(a)and the lift coefficient(b)

圖7分別給出了優(yōu)化目標(升阻比)和約束(升力系數(shù))隨迭代進行的歷程.經(jīng)優(yōu)化后，乘波體的最大升阻比為8.39，較基準構(gòu)型提升29.1%，升力系數(shù)為0.06248，略高于基準構(gòu)型的升力系數(shù)(0.2%).上述結(jié)果表明通過壓縮面優(yōu)化，可在滿足升力系數(shù)不減條件下有效改善乘波體的升阻比.此外，整個優(yōu)化歷程中的升阻比最大值為9.69，較基準構(gòu)型提升近50%，但其升力系數(shù)僅為0.04042，不滿足約束條件，這表明升力約束為緊約束.

為評估本文乘波體生成方法的有效性，從所有優(yōu)化算例中遴選出變形量達到上下邊界值，即向上(圖8(a))和向下(圖8(b))變形幅度最大的兩個構(gòu)型進行驗證.圖8分別給出了上述兩個構(gòu)型在設(shè)計條件下(0°飛行攻角)的不同截面壓力分布云圖.從圖中可以看出，兩個乘波體構(gòu)型的壓縮面形狀具有較大差別，但在設(shè)計狀態(tài)下，高壓區(qū)仍完全包裹于飛行器下部，即乘波特性明顯保持.圖9給出了上述兩個構(gòu)型乘波體尾緣處(X=2m)的截面壓力分布與基準乘波體的等值線圖比較，圖中左半部分均為基準乘波體，右半側(cè)分別為構(gòu)型A和構(gòu)型B，為清晰表達，僅給出激波面附近的壓力等值線(2200Pa～3200Pa).從圖中可見，經(jīng)壓縮面變形后，兩個乘波體的激波面均與基準外形激波面保持一致.上述結(jié)果充分表明了前述乘波體生成方法的有效性.

圖8 兩種乘波體的不同截面壓力云圖Fig.8 Pressure contours at dif f erent cross-sections of the two waveriders

圖9 乘波體尾緣壓力等值線比較Fig.9 Pressure contours comparison at the ending edge

兩個構(gòu)型的設(shè)計變量數(shù)值及設(shè)計條件下的主要氣動參數(shù)結(jié)果如表2所示.從表中可以看出，兩個構(gòu)型的最大變形量(Hmax)分別達到上下邊界.從氣動參數(shù)比較看，二者的升力系數(shù)差別約為7倍，阻力系數(shù)差別達20倍以上，升阻比亦有接近4倍的差距.相對壓心變化也超過一倍.上述結(jié)果表明壓縮面的形狀改變可導(dǎo)致乘波體的氣動性能出現(xiàn)大幅變化，設(shè)計空間獲得極大拓展.此外，結(jié)果還表明乘波體與高升阻比并不等價，若壓縮面選擇不當，盡管構(gòu)型具有乘波特性，但亦可能導(dǎo)致其升阻比大幅下降.

表2 兩種邊界構(gòu)型的設(shè)計變量和氣動參數(shù)值Table 2 Values of design variables and aerodynamic parameters of the two configuration

4 考慮黏性和邊緣鈍化的典型構(gòu)型評估

基于前述的無黏優(yōu)化設(shè)計結(jié)果，對于本文算例而言，乘波體壓縮面的設(shè)計空間可得到基本界定.為進一步分析，從前述優(yōu)化結(jié)果中選擇具有不同升阻比的 6個典型構(gòu)型 (無黏升阻比變化范圍為2.87～7.86)，采用加材料方法對其邊緣進行厚度為10mm的鈍化處理，所得外形如圖10所示.之后基于N-S方程對上述構(gòu)型的氣動性能進行了數(shù)值評估.

圖10 典型構(gòu)型外形圖Fig.10 Geometries of typical waverider configuration

計算網(wǎng)格采用多塊結(jié)構(gòu)網(wǎng)格，見圖11.為保證可比性，所有算例均基于相同的網(wǎng)格拓撲結(jié)構(gòu)和網(wǎng)格參數(shù)，總網(wǎng)格數(shù)為232萬，貼近壁面的網(wǎng)格進行了加密處理，近壁首層網(wǎng)格尺度為0.2mm，以保證邊界層的捕捉.計算采用總變差減小格式，時間推進選用隱式格式，采用k-ε湍流模型.

圖11 鈍化前緣乘波體黏性分析計算網(wǎng)格示意圖Fig.11 Grid structure for viscous analysis of blunt-edge waveriders

在設(shè)計狀態(tài)(飛行攻角0°)下，上述乘波體尾緣截面壓力分布云圖與基準構(gòu)型的比較如圖12所示.圖中左側(cè)部分均對應(yīng)基準構(gòu)型，右側(cè)部分對應(yīng)前述6種典型構(gòu)型.從圖12可以明顯看出，盡管由于鈍化和黏性的影響，前緣激波出現(xiàn)脫體，下表面高壓區(qū)有部分泄漏，但所有構(gòu)型的激波位置、激波面附近壓力分布及上表面附近的壓力分布均與基準構(gòu)型保持一致，即經(jīng)表面局部變形后的構(gòu)型所具有的乘波特性與基準構(gòu)型相同.該結(jié)果再一次驗證了本文提出的乘波體生成方法的有效性.

表3列出了上述6種典型構(gòu)型及基準乘波體構(gòu)型的黏性/無黏氣動性能評估結(jié)果.表中Cl為升力系數(shù)，Cd為阻力系數(shù)，L/D為升阻比，Xcp為縱向相對壓心系數(shù).C0表示基準乘波體構(gòu)型(黏性評估時亦對邊緣進行了鈍化)，C1～C6與圖10中給出的構(gòu)型對應(yīng)，α為飛行攻角.由于無黏優(yōu)化過程中該參數(shù)為可變量，因此在黏性評估中亦采用無黏計算時的攻角，以保證可比性.從表中的計算結(jié)果可以看出，首先，在考慮前緣鈍化及黏性效應(yīng)條件下，乘波體的升力與無黏分析結(jié)果差別不大，但阻力明顯增大.其次，乘波體的俯視投影面積不變，濕潤面積差別較小，同時邊緣鈍化的形狀也完全一致，因此較無黏結(jié)果而言，盡管阻力系數(shù)明顯增大，但對所有構(gòu)型而言，阻力系數(shù)的增量相對變化不大，均在0.05左右.

圖12 不同乘波體尾緣壓力分布云圖比較Fig.12 Pressure contours comparison at trailing edge plane of dif f erent waveriders

表3 典型構(gòu)型黏性/無黏氣動參數(shù)比較Table 3 Comparison of aerodynamic parameters based on dif f erent numerical models for typical waveriders

典型構(gòu)型的選擇是依據(jù)無黏升阻比遞增原則，因此升阻比基本呈線性增加.考慮邊緣鈍化和黏性后，由于阻力增加，升阻比明顯下降，但按構(gòu)型順序仍呈增加趨勢.當乘波體厚度較小時，鈍化邊緣的脫體激波阻力及摩擦阻力所占總阻力的比重明顯增加，因此無黏和有黏條件下升阻比差別明顯增加，這一點在構(gòu)型C5和C6的結(jié)果中體現(xiàn)尤為明顯.由此可知，盡管考慮邊緣鈍化和黏性條件的升阻比極值與無黏條件下的升阻比極值不同，但在相當大的范圍內(nèi)，無黏和有黏升阻比變化規(guī)律相同.因此實際的優(yōu)化設(shè)計中可以考慮先開展無黏優(yōu)化，之后在無黏優(yōu)化結(jié)果附近縮減所搜空間，再進行黏性優(yōu)化或分析，可以有效減小計算量，提高計算效率.此外，在考慮鈍化和黏性條件下，構(gòu)型C1和C6的升阻比相差超過一倍.相對于初始C0構(gòu)型而言，在升力不減條件下，升阻比也有超過14%的改善(構(gòu)型C6).這表明通過壓縮面的改變提升飛行器氣動性能具有極大的優(yōu)化空間.

圖13給出了上述乘波體壓縮面壓力分布云圖比較，從圖中可見，隨壓縮面外形變化，物面壓力分布差異明顯，這是導(dǎo)致相應(yīng)構(gòu)型氣動參數(shù)發(fā)生變化的原因.此外，壓縮面壓力分布的變化還導(dǎo)致俯仰力矩出現(xiàn)大幅變化，使構(gòu)型的縱向壓心產(chǎn)生偏移，如表3所示.

圖13 不同壓縮面乘波體壓力分布比較Fig.13 Pressure contours on lower surfaces of waveriders with dif f erent compression surfaces

5 結(jié)論和討論

本文給出了一種基于表面局部變形的乘波體設(shè)計方法，并在此基礎(chǔ)上結(jié)合數(shù)值優(yōu)化和計算流體力學(xué)方法評估了壓縮面變化對乘波體氣動特性的影響.結(jié)果表明，所給出的乘波體設(shè)計方法簡便易行，與增量參數(shù)化方法相結(jié)合，可在較大范圍內(nèi)修改乘波體壓縮面，且充分保持構(gòu)型的乘波特性.數(shù)值優(yōu)化結(jié)果表明，在前緣線(升力面積)不變的條件下，乘波體的升阻力系數(shù)及壓心系數(shù)均具有較大的優(yōu)化設(shè)計空間.在實際應(yīng)用中，可結(jié)合實際需求，同時考慮升阻比、穩(wěn)定性和配平特性等開展優(yōu)化設(shè)計，應(yīng)可取得明顯效果.

對于本文所給出的乘波體設(shè)計方法而言，表面變形幅度和變形區(qū)域的選擇對于構(gòu)型能否保持乘波特性十分關(guān)鍵.增量參數(shù)化設(shè)計方法基于網(wǎng)格坐標系統(tǒng)改變外形，便于變形區(qū)域的界定，因此十分適合于壓縮面優(yōu)化.但目前的方法在變形區(qū)域的選擇中主要依賴于經(jīng)驗，在后續(xù)工作中，基于空氣動力學(xué)基本理論，在不同設(shè)計條件下對變形區(qū)域及變形最大幅度的合理選擇開展研究，是一個十分值得探索的問題.

在本文工作中，乘波體的上表面采用自由來流面設(shè)計，并在優(yōu)化和分析中保持不變.但從構(gòu)型可以直觀看出，此類乘波體的容積相對較小，通過優(yōu)化，壓縮面趨于平緩，隨升阻比增加，容積減小更為明顯.而在實際中，如果考慮容積需求，乘波體的上表面只能采用壓縮面設(shè)計，這也將產(chǎn)生額外的阻力.此時應(yīng)同時考慮上下表面形狀變化開展優(yōu)化，對于這類問題，本文的思路和方法可直接推廣使用.此外，當在實際設(shè)計中考慮橫航向穩(wěn)定性時，采用本文的方法對乘波體壓縮面形狀進行修改和優(yōu)化也將十分便利.

1 葉友達.近空間高速飛行器氣動特性研究與布局設(shè)計優(yōu)化.力學(xué)進展，2009，39(6)：683-694(Ye Youda.Study on aerodynamic characteristic and design of optimization for high speed near space vehicles.Advances in Mechanics，2009，39(6)：683-694(in Chinese))

2 Ferguson F，Dasque N，Dhanasar M,et al.Waverider design and analysis.AIAA Paper 2015-3508,2015

3 Lunan DA.Waverider，a revised chronology.AIAA Paper 2015-3529,2015

4 Ferguson F，Dasque N，Dhanasar M,et al.The design,analysis and performance evaluation of waverider configuration for hypersonic vehicle applications.AIAA Paper 2015-1008,2015

5 Jones JG，Moore KC，Pike J,et al.A method for designing lifting configuration for high supersonic speeds,using axisymmetric fl w fieldsIngenieur-Archiv，1968，37：56-72

6 Rasmussen ML.Waverider configuration derived from inclined circular and elliptic cones.Journal of Spacecraft and Rockets，1980，17(5)：537-545

7 Doty RT，Rasmussen ML.Approximation for hypersonic fl w past an inclined cone.AIAA Journal，1973，11：1310-1315

8 Rasmussen ML，Clement LW.Cone-derived waveriders with longitudinal curvature.AIAA Paper 84-2100,1984

9 Takashima N，Lewis MJ.Waverider configuration based on nonaxisymmetric fl w field for engine-airframe integration，AIAA Pa-per 94-0380,1994

10 Cui K，Yang GW.Waverider configuration derived from general conical FlowfieldsActa Mechanica Sinica，2007，23(3)：247-255

11 Cui K，Yang GW.The ef f ect of conical fl wfield on the performance of waveriders at mach 6.Chinese Science Bulletin，2007，52(1)：51-64

12 Ding F，Liu J,Chen CB,et al.Novel approach for design of a waverider vehicle generated from axisymmetric supersonic fl ws past a pointed von Karman ogive.Aerospace Science and Technology，2015，42：297-308

13 呂偵軍，王江峰，伍貽兆等.多級壓縮錐導(dǎo)乘波體設(shè)計與分析.宇航學(xué)報，2015，36(5)：518-523(LZhenjun,Wang Jiangfeng,Wu Yizhao,et al.Design and analysis of multistage compression conederived waverider configurationJournal of Astronautics,2015, 36(5)：518-523(in Chinese))

14 Rodi PE.The osculating fl wfielmethod of waverider geometry generation.AIAA Paper 2005-0511,2005

15 賀旭照，倪鴻禮.密切內(nèi)錐乘波體設(shè)計方法和性能分析.力學(xué)學(xué)報，2011，43(5)：803-808(He Xuzhao,Ni Hongli.Osculating inward turning cone(OIC)wave rider-design methods and performace analysis.Chinese Journal of Theoretical and Applied Mechanics, 2011,43(5)：803-808(in Chinese))

16 賀旭照，倪鴻禮.密切曲面錐乘波體--設(shè)計方法和性能分析.力學(xué)學(xué)報，2011，43(6)：1077-1082(He Xuzhao,Ni Hongli.Osculating curved cone(occ)waverider：design methods and performance analysis.Chinese Journal of Theoretical and Applied Mechanics,2011, 43(6)：1077-1082(in Chinese))

17 Lyu YC,Jiang CW,Gao ZX,et al.Passive waverider method and its validation.AIAA Paper 2014-4346,2014

18 Hu SY,Jiang CW,Gao ZX,et al.Design of periodic cruise vehicle based on the passive waverider method.AIAA Paper 2015-4546, 2015

19 Steelant J,Langener T,Matteo FD,et al.Conceptual design of the high-speed propelled experimental fligh test vehicle HEXAFLY. AIAA Paper 2015-3539,2015

20 Pezzella G,Marini M,Reimann B,et al.Aerodynamic design analysis of the HEXAFLY-INT hypersonic glider.AIAA Paper 2015-3644,2015

21 Ding F，Liu J，Shen CB,et al.Novel inlet-airframe integration methodology for hypersonic waverider vehicles.Acta Astronautica，2015，111：178-197

22 Liu J,Ding F,Huang W,et al.Novel approach for designing a hypersonic gliding-cruising dual waverider vehicle.Acta Astronautica, 2014,102：81-88

23 Li YQ，An P，Pan CJ,et al.Integration methodology for waveriderderived hypersonic inlet and vehicle forebody.AIAA Paper 2014-3229,2014

24 Cui K，Hu SC，Li GL,et al.Conceptual design and aerodynamic evaluation of hypersonic airplane with double flankinair inlets.Science China Technology Science，2013，56(8)：1980-1988

25 胡守超,崔凱,李廣利等.基于實驗設(shè)計方法的高超聲速飛機前緣型線優(yōu)化分析.力學(xué)學(xué)報,2016,48(2)：290-299(Hu Shouchao, Cui Kai,Li Guangli,et al.Optimization and analysis of the leading edge shape for hypersonic airplanes based on DOE methods.Chinese Journal of Theoretical and Applied Mechanics,2016,48(2)：290-299(in Chinese))

26 Bowcutt KG,Anderson JD,Capriotti D.Viscous optimized hypersonic waveriders.AIAA Paper 87-0272,1987

27 Corda S,Anderson JD.Viscous optimized hypersonic waveriders designed from axisymmertirc fl w fieldsAIAA Paper 88-0396, 1988

28 張鋒濤,崔凱,楊國偉等.基于神經(jīng)網(wǎng)絡(luò)技術(shù)的乘波體優(yōu)化設(shè)計.力學(xué)學(xué)報,2009,41(3)：418-424(Zhang Fengtao,Cui Kai,Yang Guowei,et al.Optimization design of waverider based on the artificia neural networks.Chinese Journal of Theoretical and Applied Mechanics,2009,41(3)：418-424(in Chinese))

29 Lobbia MA,Suzuki K.Experimental investigation of a Mach 3.5 waverider designed using computational flui dynamics.AIAA Journal,2015,53(6)：1590-1601

30 Ryan KM,Lewis MJ,Yu KH.Comparison of robust optimization methods applied to hypersonic vehicle design.Journal of Aircraft, 2015,52(5)：1510-1523

31 Lobbia MA.Optimization of waverider-derived crew reentry vehicles using a rapid aerodynamics analysis approach.AIAA Paper 2015-0757,2015

32 Lobbia MA.Multidisciplinary design optimization of hypersonic transport configuration using waveriders.AIAA Paper 2014-2359, 2014

33 Bauer SXS.Analysis of Two Viscous Optimized Waveriders//Proceeding of the First International Hypersonic Waverider Symposium,1990,Univ.of Maryland,College Park,MD,USA

34 Charles EC.Interpretation of waverider performance data using computational fluidynamics.Journal of Aircraft,1994,31(5)：1095-1100

35 Stron R,Price K.Dif f erential evolution—a simple and efficient heuristic for global optimization over continuous spaces.Journal of Global Optimization,1997,11：341-359

EFFECT OF COMPRESSION SURFACE DEFORMATION ON AERODYNAMIC PERFORMANCES OF WAVERIDERS1)

Cui Kai2)Xu Yingzhou Xiao Yao Li Guangli

(State Key Lab of High-Temperature Gas Dynamics，Institute of Mechanics，Chinese Academy of Sciences，Beijing100190，China)

(School of Engineering Science，University of Chinese Academy of Sciences，Beijing100049，China)

A waverider is a type of hypersonic lifting body that has the entire bow shock underneath the body as well as attached to the leading edge when flyin at its design Mach number.Present research for improving the aerodynamic performance of waveriders mainly focused on searching an optimal profil of the leading edge on the condition of given a specifi generatingfl wfield Inordertofurtherextendthedesignspaceofwaveriders,anoveldesignmethodthatisbased on a local shape deformation technique is presented in this paper.Moreover,an inviscid analysis-based optimization study was carried out to research the ef f ect of compression surface deformation on aerodynamic performances of waveriders by integrating the increment-based parameterization method,the computational flui dynamic analysis,and the dif f erential evolutionalgorithm.Afterwards,sixselectedwaveriderconfiguration werepolishedtobluntleadingedges,andthentheir aerodynamic performances were evaluated by solving the Navier-Stokes equations.The results show that both theL/Dand the relative pressure center coefficient of the waveriders produce significan changes with the variation of compressionsurface shape.Among all waveriders,the maximal dif f erence of theL/Dis more than double.Even by considering the lift constraint,the increment of theL/Dis more than 14 percent in comparison with the baseline configuration In addition, the value of relative migration of the relative pressure center coefficients is remarkable.

waverider,aerodynamic shape optimization,increment-based parameterization method,dif f erential evolution algorithm,computational flui dynamics

O354.4,V211.24

A doi：10.6052/0459-1879-16-041

2016-02-01收稿，2016-10-11錄用,2016-10-13網(wǎng)絡(luò)版發(fā)表.

1)國家自然科學(xué)基金資助項目(11372324，11572333).

2)崔凱，副研究員，主要研究方向：飛行器構(gòu)型設(shè)計和優(yōu)化.E-mail：kcui@imech.ac.cn

崔凱，徐應(yīng)洲，肖堯，李廣利.乘波體壓縮面變化對其氣動性能影響分析.力學(xué)學(xué)報,2017,49(1)：75-83

Cui Kai,Xu Yingzhou,Xiao Yao,Li Guangli.E ff ect of compression surface deformation on aerodynamic performances of waveriders.Chinese Journal of Theoretical and Applied Mechanics,2017,49(1)：75-83