李春鵬，錢戰(zhàn)森，孫俠生

1.中國航空工業(yè)空氣動(dòng)力研究院 高速高雷諾數(shù)氣動(dòng)力航空科技重點(diǎn)實(shí)驗(yàn)室，沈陽 110034 2.中國航空研究院，北京 100012

遠(yuǎn)程民機(jī)具有巡航速度高、使用航程大、攜帶燃油多等特點(diǎn)，燃油消耗產(chǎn)生的重量變化導(dǎo)致巡航過程大多以非設(shè)計(jì)狀態(tài)飛行，機(jī)翼翼面流場極易出現(xiàn)激波，由此導(dǎo)致阻力激增，不但降低了飛行效率，而且還因臨近抖振邊界而影響飛行安全。當(dāng)前，通過不斷增加飛行高度的階梯巡航方式在一定程度上能夠減小巡航升力變化幅度，但隨著航空市場的快速發(fā)展，未來因空域緊張而增加的對(duì)飛行高度限制將會(huì)使上述問題更加突出?？紤]基于固定外形機(jī)翼多點(diǎn)綜合優(yōu)化設(shè)計(jì)僅僅是性能折中的一種手段，為有效提升機(jī)翼在整個(gè)巡航階段的綜合氣動(dòng)性能，研究者逐漸將重點(diǎn)集中于后緣外形隨飛行條件變化而不斷調(diào)整的變彎度機(jī)翼。應(yīng)用變彎度機(jī)翼技術(shù)不僅能夠通過后緣變形來調(diào)整機(jī)翼載荷分布實(shí)現(xiàn)巡航效率提升、抖振邊界擴(kuò)展，還能進(jìn)行陣風(fēng)減緩控制，并改善機(jī)動(dòng)性能，進(jìn)而減小機(jī)翼結(jié)構(gòu)重量、降低燃油消耗并節(jié)省成本，對(duì)未來“綠色航空”發(fā)展具有重要的促進(jìn)作用。

通過后緣外形變形來改善非設(shè)計(jì)點(diǎn)氣動(dòng)性能的民機(jī)變彎度機(jī)翼技術(shù)，最早由美國國家航空航天局（National Aeronautics Space and Administration，NASA）德萊頓（Dryden）飛行研究中心于20世紀(jì)70年代提出，針對(duì)L-1011三星客機(jī)（L-1011 TriStar）研究結(jié)果表明，變彎度機(jī)翼技術(shù)名義巡航飛行狀態(tài)可獲得1% ～ 3%的減阻收益［1-2］。波音公司和空客公司從20世紀(jì)80年代開始評(píng)估變彎度機(jī)翼概念，并結(jié)合B777-200ER、A330/A340等開展了應(yīng)用研究［3-4］，目前基于后緣剛性舵面的變彎度機(jī)翼技術(shù)已在最新的遠(yuǎn)程民機(jī)B787和A350等機(jī)型上獲得了應(yīng)用［5］。國內(nèi)針對(duì)后緣剛性舵面的變彎度機(jī)翼技術(shù)也開展了大量的研究［6-9］。

與基于后緣剛性舵面偏轉(zhuǎn)的變彎度機(jī)翼相比，基于柔性結(jié)構(gòu)/變形機(jī)構(gòu)的機(jī)翼后緣連續(xù)變彎度技術(shù)具有更大的變形空間、更強(qiáng)的翼面載荷調(diào)節(jié)能力，以及更小的飛行阻力。國內(nèi)外學(xué)者圍繞連續(xù)變彎后緣的氣動(dòng)收益問題開展了廣泛的研究。Lyu等［10］采用高精度伴隨求解器，針對(duì)CRM機(jī)翼開展了連續(xù)變彎后緣的氣動(dòng)設(shè)計(jì)研究，完成了407個(gè)飛行狀態(tài)點(diǎn)后緣變形外形優(yōu)化，系統(tǒng)分析了后緣連續(xù)變彎在設(shè)計(jì)點(diǎn)和非設(shè)計(jì)點(diǎn)的減阻收益及其對(duì)燃油消耗的影響。Burdette等［11］采用高精度求解器，針對(duì)uCRM模型開展了旨在減阻減重的后緣連續(xù)變彎優(yōu)化設(shè)計(jì)研究，明確了后緣連續(xù)變彎節(jié)省燃油消耗的2種控制機(jī)制。Chen等［12］采用低精度氣動(dòng)結(jié)構(gòu)分析方法，針對(duì)波音公司在亞音速超級(jí)綠色飛機(jī)研究（Subsonic Ultra Green Aircraft Research，SUGAR）項(xiàng)目中提出桁架支撐機(jī)翼布局，并結(jié)合連續(xù)變彎后緣襟翼開展了減阻減重優(yōu)化設(shè)計(jì)研究。Ting等［13］以NASA 通用運(yùn)輸機(jī)標(biāo)模（Generic Transport Model，GTM）模型，采用耦合跨聲速小擾動(dòng)方程和邊界層修正的渦格法，針對(duì)多個(gè)典型飛行狀態(tài)開展了連續(xù)變彎后緣襟翼減阻優(yōu)化設(shè)計(jì)研究。Bui［14］利用基于RANS（Reynolds Averaged Navier-Stokes）方程的計(jì)算流體力學(xué)（Computational Fluid Dynamics，CFD）求解器，針對(duì)“灣流”III公務(wù)機(jī)機(jī)翼開展了自適應(yīng)后緣（Adaptive Compliant Trailing Edge，ACTE）激波邊界層干擾誘導(dǎo)分離控制研究。Bartels等［15］采用降階模型從主動(dòng)氣彈控制角度針對(duì)桁架支撐機(jī)翼布局開展了連續(xù)變彎后緣襟翼應(yīng)用研究。郭同彪［16］、聶瑞［17］、梁海朝［18］等針對(duì)跨聲速翼型研究了后緣連續(xù)變彎減阻收益及影響機(jī)理。王一凡等［19］針對(duì)自然層流翼型研究了考慮整個(gè)巡航段氣動(dòng)性能的后緣連續(xù)變彎優(yōu)化設(shè)計(jì)方法。李春鵬等［20］針對(duì)帶有連續(xù)變彎后緣跨聲速翼型開展了考慮后緣自適應(yīng)變形的基本翼型多點(diǎn)優(yōu)化設(shè)計(jì)研究。王斌等［21］針對(duì)寬體客機(jī)變彎度機(jī)翼開展了后緣連續(xù)變彎氣動(dòng)優(yōu)化設(shè)計(jì)研究，對(duì)比分析了后緣變彎與前后緣同時(shí)變彎的氣動(dòng)差異。郭同彪等［22］針對(duì)民機(jī)機(jī)翼-機(jī)身-平尾構(gòu)型開展了考慮配平約束的后緣連續(xù)變彎度機(jī)翼氣動(dòng)優(yōu)化設(shè)計(jì)研究。雷銳午等［23］基于典型民機(jī)標(biāo)模開展了考慮抖振性能約束的機(jī)翼后緣變彎減阻優(yōu)化設(shè)計(jì)研究。

結(jié)合二維翼型和三維機(jī)翼開展的相關(guān)研究已充分驗(yàn)證了連續(xù)變彎后緣在改善整個(gè)巡航階段氣動(dòng)性能的技術(shù)優(yōu)勢，但在工程應(yīng)用中，連續(xù)變彎后緣氣動(dòng)外形設(shè)計(jì)還應(yīng)考慮變形結(jié)構(gòu)變形能力、飛行控制系統(tǒng)操縱邏輯等因素對(duì)弦向和展向變形外形的約束。本文以遠(yuǎn)程民機(jī)標(biāo)模為研究對(duì)象，基于機(jī)翼后緣指關(guān)節(jié)變形結(jié)構(gòu)，從保證后緣隨升力單調(diào)偏轉(zhuǎn)角度，開展連續(xù)變彎后緣外形變形矩陣氣動(dòng)優(yōu)化設(shè)計(jì)研究。首先結(jié)合機(jī)翼某展向站位翼型剖面完成連續(xù)變彎后緣的弦向變形方式篩選，然后針對(duì)全機(jī)構(gòu)型，對(duì)比分析了后緣協(xié)調(diào)變形對(duì)變彎度機(jī)翼氣動(dòng)收益的影響規(guī)律，從多點(diǎn)關(guān)聯(lián)和階梯變形角度構(gòu)建了巡航任務(wù)剖面變形矩陣，并進(jìn)一步研究了抖振抑制和阻力發(fā)散點(diǎn)減阻等非巡航任務(wù)剖面下的變形矩陣設(shè)計(jì)方法。

1 研究模型

基于遠(yuǎn)程民機(jī)標(biāo)?；練鈩?dòng)外形，參考機(jī)翼后緣原有常規(guī)舵面，在80%翼展范圍內(nèi)布置弦向和展向均與機(jī)翼光滑連接的柔性變彎后緣，設(shè)定7個(gè)變彎控制剖面以保證后緣具有足夠的變形能力，典型變彎度機(jī)翼后緣變形外形與基本外形對(duì)比見圖1。

圖1 變彎度機(jī)翼后緣變形外形對(duì)比Fig.1 Comparison of variable camber wing trailing edge shape

指關(guān)節(jié)變形結(jié)構(gòu)作為一種典型變彎后緣結(jié)構(gòu)形式，由能相互獨(dú)立偏轉(zhuǎn)的分段剛性組成，分段之間通過柔性蒙皮光滑連接，具有結(jié)構(gòu)簡單、剛度強(qiáng)度大、承載能力強(qiáng)等突出優(yōu)勢［24］。參考3段式指關(guān)節(jié)變形結(jié)構(gòu)，保留其末段具有一定弦向占比且外形不可變的特征，忽略其內(nèi)段結(jié)構(gòu)對(duì)后緣外形的限制，定義變彎控制剖面外形曲線（圖2）：根據(jù)變彎后緣剖面上下翼面確定基本外形中弧線，將中弧線按一定比例分為3段，考慮結(jié)構(gòu)變形能力限制，確定末段中弧線占比25%，且保留與之相對(duì)應(yīng)的上下翼面基本外形。結(jié)合3段中弧線的前端點(diǎn)確定后緣變彎外形控制轉(zhuǎn)軸，各段中弧線除繞其前端轉(zhuǎn)軸偏轉(zhuǎn)外，還隨其前端中弧線的轉(zhuǎn)動(dòng)而轉(zhuǎn)動(dòng)。將旋轉(zhuǎn)前的上下翼面基本外形前端與旋轉(zhuǎn)后的末段中弧線對(duì)應(yīng)上下翼面基本外形前端利用二次樣條曲線光順連接形成變形后的變彎控制剖面。

圖2 變彎控制剖面幾何外形定義Fig.2 Geometric shape definition of deformation shape control section

變彎度機(jī)翼后緣外形變形在調(diào)整機(jī)翼翼面載荷分布的同時(shí)也會(huì)產(chǎn)生一定的俯仰力矩增量，為了計(jì)及配平俯仰力矩帶來的阻力影響，需在開展氣動(dòng)特性評(píng)估之前，確定與變彎度機(jī)翼后緣變形外形的相對(duì)應(yīng)的全動(dòng)平尾配平偏度，全動(dòng)平尾及其轉(zhuǎn)軸見圖3。

圖3 全動(dòng)平尾及其轉(zhuǎn)軸Fig.3 All-movable horizontal tail and its rotation axis

2 計(jì)算方法及驗(yàn)證

2.1 計(jì)算方法

數(shù)值模擬采用三維N-S方程，在直角坐標(biāo)系中三維N-S方程守恒形式可表示為

式中：w為狀態(tài)矢量；f為無黏（對(duì)流）通矢量項(xiàng)；fv為黏性（耗散）通矢量項(xiàng)。

采用Menter’s SST湍流模型，利用有限體積法控制方程離散，通過多重網(wǎng)格加速收斂。為提高計(jì)算效率，采用半模計(jì)算。遠(yuǎn)場與計(jì)算域中心距離為20倍機(jī)身長度，網(wǎng)格為多塊結(jié)構(gòu)化網(wǎng)格，在飛行器表面附近生成棱柱層網(wǎng)格以捕捉邊界層內(nèi)的流動(dòng)特征，第1層網(wǎng)格高度為飛機(jī)參考長度的10-6，壁面網(wǎng)格法向增長率1.25，棱柱層數(shù)33層，半模網(wǎng)格點(diǎn)數(shù)800萬左右。物面采用無滑移條件，對(duì)稱面采用對(duì)稱邊界條件，遠(yuǎn)場由當(dāng)?shù)匾痪S黎曼不變量確定。

2.2 計(jì)算驗(yàn)證

選擇第5屆阻力預(yù)測會(huì)議跨聲速民機(jī)標(biāo)模CRM翼身組合體構(gòu)型來驗(yàn)證方法的可靠性。采用結(jié)構(gòu)網(wǎng)格計(jì)算，網(wǎng)格點(diǎn)為800萬左右，馬赫數(shù)Ma= 0.85，雷諾數(shù)Re= 5.0×106，升力系數(shù)CL= 0.486。計(jì)算模型及網(wǎng)格見圖4。

圖4 CRM翼身組合體模型及計(jì)算網(wǎng)格Fig.4 CRM wing-body model and computational mesh

4個(gè)展向站位機(jī)翼剖面壓力分布與試驗(yàn)［24］對(duì)比如圖5所示［25］，圖中：2yb為當(dāng)?shù)卣瓜蛘疚痪鄬?duì)稱面距離與半展長的比值。計(jì)算得到的不同站位翼面壓力系數(shù)分布形態(tài)、前緣吸力峰、激波強(qiáng)度等關(guān)鍵流場特征與試驗(yàn)基本一致，僅靠近翼梢站位剖面的激波位置略有偏差，說明本文所采用的方法具有較高精度，可用于遠(yuǎn)程民機(jī)變彎度機(jī)翼的氣動(dòng)分析和優(yōu)化研究。

圖5 CRM翼身組合體計(jì)算與試驗(yàn)翼面壓力系數(shù)對(duì)比［25］Fig.5 Pressure coefficient comparison of computational result and test data of CRM wing body［25］

3 基于代理模型的優(yōu)化方法

根據(jù)幾何外形成型方式確定變彎度機(jī)翼后緣外形優(yōu)化設(shè)計(jì)變量，利用通用優(yōu)化平臺(tái)參數(shù)化建模、網(wǎng)格自動(dòng)化生成和高精度流場計(jì)算等模塊實(shí)現(xiàn)設(shè)計(jì)變量對(duì)應(yīng)響應(yīng)值的自動(dòng)求解，選取代理優(yōu)化算法完成給定狀態(tài)的設(shè)計(jì)變量尋優(yōu)，結(jié)合后緣外形間變形關(guān)聯(lián)關(guān)系制定氣動(dòng)優(yōu)化策略，進(jìn)而完成變彎度機(jī)翼后緣外形變形矩陣氣動(dòng)優(yōu)化流程搭建。

3.1 代理優(yōu)化算法

代理優(yōu)化算法是在基于試驗(yàn)設(shè)計(jì)構(gòu)建初始代理模型的基礎(chǔ)上，按照一定的優(yōu)化加點(diǎn)準(zhǔn)則（如MSP、EI、PI、MSE、LCB等），采用傳統(tǒng)優(yōu)化算法求解相應(yīng)的子優(yōu)化問題，以較小計(jì)算代價(jià)實(shí)現(xiàn)最優(yōu)解的高效預(yù)測，并將預(yù)測最優(yōu)解及其真實(shí)響應(yīng)值添加到用于構(gòu)建代理模型的數(shù)據(jù)集，通過代理模型不斷更新來實(shí)現(xiàn)新增樣本序列向最優(yōu)解收斂［26］，代理優(yōu)化算法典型框架如圖6所示。

圖6 代理優(yōu)化算法典型框架［26］Fig.6 Framework of surrogate-based optimization algorithm［26］

代理優(yōu)化算法既可借助代理模型替換復(fù)雜、費(fèi)時(shí)的高精度數(shù)值計(jì)算，又可基于歷史數(shù)據(jù)求解優(yōu)化加點(diǎn)準(zhǔn)則定義的子優(yōu)化問題，在提高優(yōu)化效率和改善全局搜索能力等方面均具有明顯的優(yōu)勢。

本文選用拉丁超立方方法進(jìn)行試驗(yàn)設(shè)計(jì)，代理模型為Kriging模型，利用EI、MSP等加點(diǎn)準(zhǔn)則完成優(yōu)化過程中的代理模型更新，采用遺傳算法與序列二次規(guī)劃相結(jié)合的方法對(duì)優(yōu)化加點(diǎn)準(zhǔn)則定義的子優(yōu)化問題尋優(yōu)［27］。

3.2 氣動(dòng)優(yōu)化策略

區(qū)別于幾何形狀固定的機(jī)翼單目標(biāo)優(yōu)化或多目標(biāo)優(yōu)化問題，變彎度機(jī)翼后緣外形變形矩陣的氣動(dòng)優(yōu)化設(shè)計(jì)可理解為一種具有耦合關(guān)系的單目標(biāo)優(yōu)化設(shè)計(jì)問題，盡管變形矩陣內(nèi)單個(gè)設(shè)計(jì)點(diǎn)對(duì)應(yīng)的后緣外形氣動(dòng)優(yōu)化是單目標(biāo)問題，但受變形關(guān)聯(lián)關(guān)系約束，其外形優(yōu)化設(shè)計(jì)空間選取與變形矩陣內(nèi)其他設(shè)計(jì)點(diǎn)的后緣外形相關(guān)。在滿足后緣外形變形結(jié)構(gòu)約束的前提下，為了充分挖掘變彎度機(jī)翼后緣變形的氣動(dòng)收益，針對(duì)具有耦合關(guān)系的單目標(biāo)優(yōu)化問題制定的氣動(dòng)優(yōu)化策略為

1） 根據(jù)遠(yuǎn)程民機(jī)變彎度機(jī)翼后緣外形變形需求確定設(shè)計(jì)狀態(tài)矩陣。

2） 確定矩陣內(nèi)關(guān)鍵設(shè)計(jì)狀態(tài)，根據(jù)優(yōu)先程度完成全部設(shè)計(jì)狀態(tài)的分層。

3） 根據(jù)設(shè)計(jì)狀態(tài)分層關(guān)系逐步調(diào)整后緣外形優(yōu)化設(shè)計(jì)變量的取值范圍，并以此為基礎(chǔ)完成后緣外形變形矩陣的氣動(dòng)優(yōu)化。

4） 綜合對(duì)比變形矩陣內(nèi)的后緣外形，根據(jù)氣動(dòng)收益及變形關(guān)聯(lián)關(guān)系調(diào)整后緣變形外形。

此外，為提高氣動(dòng)優(yōu)化效率，以后緣基本外形阻力系數(shù)為基準(zhǔn)，只在具有減阻收益時(shí)進(jìn)行配平計(jì)算，否則將其阻力系數(shù)疊加某一特定值（明顯大于阻力系數(shù)的大量，如99）做壞點(diǎn)處理。

變彎度機(jī)翼后緣外形變形矩陣氣動(dòng)優(yōu)化流程見圖7，其中設(shè)計(jì)狀態(tài)矩陣主要是以升力系數(shù)階梯變化為特征的巡航任務(wù)剖面和以抖振、阻力發(fā)散為特征的非巡航任務(wù)剖面兩種情況確定，結(jié)合后緣外形變形需求，按變形程度篩選矩陣內(nèi)關(guān)鍵設(shè)計(jì)狀態(tài)，按優(yōu)先程度完成設(shè)計(jì)狀態(tài)分層，優(yōu)化設(shè)計(jì)變量的取值范圍隨著優(yōu)先層級(jí)的降低而逐漸收縮。對(duì)后緣變形外形完成幾何參數(shù)化建模，確定設(shè)計(jì)變量及其取值范圍，按最大的設(shè)計(jì)變量空間完成高層級(jí)設(shè)計(jì)狀態(tài)對(duì)應(yīng)的試驗(yàn)設(shè)計(jì)，獲取代理模型構(gòu)建所需的初始樣本點(diǎn)集，通過自動(dòng)化模型和網(wǎng)格生成、數(shù)值計(jì)算完成樣本對(duì)應(yīng)的氣動(dòng)性能參數(shù)求解。利用初始樣本點(diǎn)集完成代理模型構(gòu)建后，結(jié)合代理模型加點(diǎn)準(zhǔn)則定義代理子優(yōu)化問題并完成代理模型參數(shù)尋優(yōu)，將得出的待加點(diǎn)作為更新樣本在獲取其氣動(dòng)性能精確解后加入樣本點(diǎn)集并完成代理模型更新，以更新樣本的氣動(dòng)性能提升程度作為收斂條件，如果相鄰迭代樣本的氣動(dòng)性能差量小于規(guī)定值即完成優(yōu)化，輸出由更新樣本組成的后緣變形外形氣動(dòng)優(yōu)化解集，解集選取范圍為規(guī)定值范圍內(nèi)最優(yōu)氣動(dòng)性能對(duì)應(yīng)的樣本。結(jié)合高層級(jí)設(shè)計(jì)狀態(tài)優(yōu)化結(jié)果，按層級(jí)關(guān)系依次確定其他設(shè)計(jì)狀態(tài)對(duì)應(yīng)的設(shè)計(jì)變量取值范圍，并完成后緣變形外形氣動(dòng)優(yōu)化解集求解。最后圍繞設(shè)計(jì)狀態(tài)矩陣，兼顧氣動(dòng)收益和變形關(guān)聯(lián)關(guān)系，從保證各設(shè)計(jì)狀態(tài)后緣外形曲線相似的角度，對(duì)后緣變形控制剖面偏度調(diào)整，進(jìn)而確定變彎度機(jī)翼后緣外形變形矩陣。

圖7 后緣外形變形矩陣氣動(dòng)優(yōu)化設(shè)計(jì)流程Fig.7 Aerodynamic optimization design process of trailing edge deformation matrix

4 后緣變彎控制剖面變形方式影響

參考3段式指關(guān)節(jié)變形結(jié)構(gòu)確定的變彎控制剖面最多需要4個(gè)變量來描述，包括1個(gè)中弧線比例分段變量和3個(gè)分段中弧線旋轉(zhuǎn)角度變量。選取不同控制變量組合的后緣外形變形能力對(duì)比見圖8，其中4個(gè)變量的順序分別為中弧線比例分段變量、前段中弧線旋轉(zhuǎn)角度變量、中段中弧線旋轉(zhuǎn)角度變量和后段中弧線旋轉(zhuǎn)角度變量，單變量是3段中弧線共用的旋轉(zhuǎn)角度變量，而將中弧線比例分段變量固定為0.45。相對(duì)采用單變量取值為5描述的變彎控制剖面變形外形，采用4變量組合描述的變彎控制剖面變形外形在后緣點(diǎn)接近的情況下具有更大范圍的后緣變形控制能力。

圖8 后緣外形變形方案變形控制能力對(duì)比Fig.8 Deformation control capability comparison of different trailing edge schemes

選取某遠(yuǎn)程民機(jī)標(biāo)模機(jī)翼典型展向站位翼型剖面，結(jié)合不同升力系數(shù)的二維流場后緣外形變形方式的影響。按照后掠翼理論確定二維翼型設(shè)計(jì)馬赫數(shù)Ma= 0.7、雷諾數(shù)Re= 1.0×107，設(shè)計(jì)升力系數(shù)CL及設(shè)計(jì)變量取值范圍如表1所示。

二維翼型后緣變形外形氣動(dòng)優(yōu)化設(shè)計(jì)結(jié)果對(duì)比如圖9所示，優(yōu)化方案對(duì)應(yīng)二維翼型后緣外形設(shè)計(jì)變量取值如表2所示，典型升力系數(shù)下后緣外形對(duì)比如圖10所示，可以看出單變量及4變量變形方式雖然變形能力不同，但兩者獲得的氣動(dòng)收益相當(dāng)。相對(duì)升力系數(shù)CL= 0.7對(duì)應(yīng)的翼型設(shè)計(jì)點(diǎn)，升力系數(shù)偏離越遠(yuǎn)，后緣變形幅度越大，所獲得的氣動(dòng)收益也越大。對(duì)于小升力系數(shù)對(duì)應(yīng)的后緣上偏狀態(tài)，2種變形方式獲得的幾何外形略有差異，而對(duì)于大升力系數(shù)對(duì)應(yīng)的后緣上偏狀態(tài)，兩種變形方式獲得的幾何外形基本一致。

表2 二維翼型后緣外形優(yōu)化設(shè)計(jì)結(jié)果Table 2 Aerodynamic optimization results of 2D airfoil trailing edge deformation shapes

圖9 二維翼型后緣變形外形升阻比曲線Fig.9 Lift-to-drag ratio curves of different airfoil trailing edge deformation shape

圖10 二維翼型后緣優(yōu)化外形對(duì)比Fig.10 Aerodynamic optimization geometric comparison of 2D airfoil trailing edge deformation shapes

以上分析表明，對(duì)于3段中弧線旋轉(zhuǎn)角度相同的單變量控制方式，盡管變形能力相對(duì)有限，但外形變形規(guī)律更符合升力大幅變化情況下的流場控制需求，更適用于變彎度機(jī)翼后緣變形外形的氣動(dòng)優(yōu)化設(shè)計(jì)研究。

5 巡航任務(wù)剖面變形矩陣

對(duì)于飛行馬赫數(shù)和高度固定的巡航任務(wù)剖面，遠(yuǎn)程民機(jī)在定常平飛過程中的升力系數(shù)會(huì)隨著燃油消耗而逐漸減小，巡航任務(wù)剖面起始點(diǎn)燃油最多，對(duì)應(yīng)升力系數(shù)最大，巡航任務(wù)剖面終止點(diǎn)燃油最少，對(duì)應(yīng)升力系數(shù)最小。

某遠(yuǎn)程民機(jī)標(biāo)模巡航馬赫數(shù)Ma= 0.85，巡航高度H= 11 km，以升力系數(shù)CL= 0.48對(duì)應(yīng)的半油狀態(tài)作為基準(zhǔn)確定巡航任務(wù)剖面中點(diǎn)，按20%的燃油變化量確定巡航任務(wù)剖面起始及終止點(diǎn)，按10%的燃油變化量確定巡航任務(wù)剖面前段及后段中點(diǎn)，得出變彎度機(jī)翼后緣變形外形巡航任務(wù)剖面設(shè)計(jì)狀態(tài)矩陣見表3。

表3 巡航任務(wù)剖面設(shè)計(jì)狀態(tài)矩陣Table 3 Aerodynamic design matrix of cruise flight profile

5.1 單點(diǎn)最優(yōu)變形矩陣

為明確變彎度機(jī)翼后緣外形變形在不同升力系數(shù)下的最大減阻空間，分別針對(duì)巡航任務(wù)剖面中不同燃油變化對(duì)應(yīng)的設(shè)計(jì)升力系數(shù)，選取相同的設(shè)計(jì)變量空間，開展變彎度機(jī)翼后緣變形外形氣動(dòng)優(yōu)化設(shè)計(jì)，形成單點(diǎn)最優(yōu)變形矩陣。

變彎度機(jī)翼后緣外形的變彎控制剖面采用單變量控制方式描述，共計(jì)7個(gè)設(shè)計(jì)變量，每個(gè)變量的取值范圍均為-5 ～ 5，非配平及配平狀態(tài)下的優(yōu)化結(jié)果氣動(dòng)特性對(duì)比見表4，表中：ΔCD為后緣變形外形相對(duì)基本外形的減阻量，α為飛行迎角，δtail為平尾偏度。從優(yōu)化結(jié)果可知，非配平及配平狀態(tài)下的變彎度機(jī)翼后緣外形變形減阻規(guī)律基本一致，且考慮配平約束后的減阻收益略大。相對(duì)升力系數(shù)CL= 0.48對(duì)應(yīng)的半油基準(zhǔn)狀態(tài)，減阻收益隨著升力系數(shù)變化量的增加而增大，但主要的減阻收益集中在升力系數(shù)增加的高升力狀態(tài)，對(duì)于同等幅度的升力系數(shù)變化量，高升力系數(shù)時(shí)后緣外形變形所能獲得的減阻量約為低升力系數(shù)時(shí)的7 ～ 8倍。

表4 單點(diǎn)最優(yōu)變形矩陣優(yōu)化結(jié)果氣動(dòng)特性對(duì)比（Ma =0.85，H = 11 km）Table 4 Aerodynamic characteristics comparison of single-point optimal deformation matrix optimization results（Ma = 0.85，H = 11 km）

單點(diǎn)最優(yōu)變形矩陣后緣外形變彎控制剖面偏度對(duì)比見圖11，圖中：正值代表相應(yīng)變彎控制剖面下偏，負(fù)值代表相應(yīng)控制剖面上偏，下同。非配平及配平狀態(tài)下的變彎控制剖面偏轉(zhuǎn)規(guī)律基本一致，隨著設(shè)計(jì)升力系數(shù)的增加，內(nèi)側(cè)5個(gè)變彎控制剖面逐漸下偏，而外側(cè)2個(gè)變彎控制剖面則先上偏后下偏。對(duì)于非配平狀態(tài)，內(nèi)側(cè)5個(gè)變彎控制剖面在固定升力系數(shù)下的偏度相差較小，外側(cè)2個(gè)變彎控制剖面的偏度在低升力系數(shù)時(shí)與內(nèi)側(cè)相差不大，但在高升力系數(shù)狀態(tài)，特別是相對(duì)基準(zhǔn)升力增加20%時(shí)的偏度相對(duì)內(nèi)側(cè)明顯增大。對(duì)于配平狀態(tài)，以中間變彎控制剖面4為分界點(diǎn)，后緣外形變彎控制剖面偏度在固定升力系數(shù)下整體呈先減小后增加的趨勢，這種展向偏度變化也減小了外側(cè)變彎控制剖面在高升力系數(shù)時(shí)的絕對(duì)偏度。整體來看，對(duì)于完全不考慮后緣外形變形關(guān)聯(lián)關(guān)系所獲得的變形矩陣，在巡航飛行過程中，隨著燃油的不斷消耗，不同后緣外形變彎控制剖面在偏轉(zhuǎn)控制方面沒有明顯的規(guī)律性，主要表現(xiàn)為各剖面偏度增量的差別較大，在個(gè)別狀態(tài)還存在反向偏轉(zhuǎn)，后緣變形控制邏輯較為復(fù)雜，難以滿足實(shí)際工程應(yīng)用要求。

圖11 單點(diǎn)最優(yōu)變形矩陣后緣外形變彎控制剖面偏度對(duì)比（Ma = 0.85，H = 11 km）Fig.11 Deflection angles comparison of trailing edge deformation shape control sections in singlepoint optimal deformation matrix（Ma = 0.85，H = 11 km）

5.2 多點(diǎn)關(guān)聯(lián)變形矩陣

基于配平狀態(tài)下的單點(diǎn)最優(yōu)變形矩陣氣動(dòng)優(yōu)化設(shè)計(jì)結(jié)果，考慮后緣外形變形關(guān)聯(lián)關(guān)系，開展變彎度機(jī)翼后緣變形外形氣動(dòng)優(yōu)化設(shè)計(jì)，形成多點(diǎn)關(guān)聯(lián)變形矩陣。

構(gòu)建多點(diǎn)關(guān)聯(lián)變形矩陣的關(guān)鍵在不同升力系數(shù)對(duì)應(yīng)后緣變形外形氣動(dòng)優(yōu)化設(shè)計(jì)變量空間的選取。考慮到后緣外形變形在小升力系數(shù)下的減阻能力有限，額外增加設(shè)計(jì)約束難以獲得足夠的減阻收益，因此以巡航任務(wù)剖面最小設(shè)計(jì)升力系數(shù)狀態(tài)對(duì)應(yīng)的后緣變形外形為約束依次設(shè)定設(shè)計(jì)變量的取值范圍。

針對(duì)巡航任務(wù)剖面設(shè)計(jì)狀態(tài)矩陣，配平條件下的多點(diǎn)關(guān)聯(lián)變形矩陣氣動(dòng)優(yōu)化按2種思路展開（圖12）：①按單點(diǎn)最優(yōu)變形矩陣確定設(shè)計(jì)變量最大值，并將其作為不同設(shè)計(jì)點(diǎn)共用的設(shè)計(jì)空間上邊界，以鄰近較小升力系數(shù)設(shè)計(jì)點(diǎn)對(duì)應(yīng)的設(shè)計(jì)變量取值為設(shè)計(jì)空間下邊界，按照設(shè)計(jì)升力系數(shù)從小到大的順序完成多點(diǎn)關(guān)聯(lián)變形矩陣的氣動(dòng)優(yōu)化設(shè)計(jì)，優(yōu)化過程中下邊界逐漸縮?。虎诮Y(jié)合單點(diǎn)最優(yōu)變形矩陣中最大及最小升力系數(shù)設(shè)計(jì)點(diǎn)確定設(shè)計(jì)空間的初始上下邊界，按照設(shè)計(jì)升力系數(shù)從兩端到中間的順序依次完成多點(diǎn)關(guān)聯(lián)變形矩陣的氣動(dòng)優(yōu)化設(shè)計(jì)，優(yōu)化過程中上下邊界同時(shí)縮小。

圖12 多點(diǎn)關(guān)聯(lián)變形矩陣氣動(dòng)優(yōu)化設(shè)計(jì)思路Fig.12 Aerodynamic optimization design ideas of multi-point coupling deformation matrix

按照2種設(shè)計(jì)變量空間選取思路獲得的配平狀態(tài)多點(diǎn)關(guān)聯(lián)變形矩陣氣動(dòng)優(yōu)化設(shè)計(jì)結(jié)果氣動(dòng)特性對(duì)比見表5，表中：ΔCD為后緣變形外形相對(duì)基本外形的減阻量，ΔDCD為后緣變形外形相對(duì)單點(diǎn)最優(yōu)變形外形的增阻量。相比于單點(diǎn)最優(yōu)變形矩陣，下邊界逐漸縮小方式獲得的變形矩陣最大阻力系數(shù)增量僅為0.000 05，上下邊界同時(shí)縮小方式獲得的變形矩陣最大阻力系數(shù)增量僅為0.000 01。綜合優(yōu)化算法尋優(yōu)能力等因素限制，可認(rèn)為考慮后緣外形變形關(guān)聯(lián)關(guān)系的設(shè)計(jì)變量空間縮減并不會(huì)產(chǎn)生明顯的減阻收益損失。

表5 配平狀態(tài)多點(diǎn)關(guān)聯(lián)變形矩陣優(yōu)化結(jié)果氣動(dòng)特性對(duì)比（Ma = 0.85，H = 11 km）Table 5 Aerodynamic characteristics comparison of multi-point coupling deformation matrix optimization results under trimmed condition（Ma = 0.85，H = 11 km）

配平狀態(tài)多點(diǎn)關(guān)聯(lián)變形矩陣后緣外形變彎控制剖面偏度對(duì)比如圖13所示。盡管設(shè)計(jì)變量空間選取方式不同，會(huì)形成差別較大的后緣外形變形矩陣，但相比于單點(diǎn)最優(yōu)變形矩陣，多點(diǎn)關(guān)聯(lián)變形矩陣內(nèi)不同設(shè)計(jì)點(diǎn)對(duì)應(yīng)的后緣變彎控制剖面偏度具有明顯規(guī)律性，即隨著巡航過程中升力系數(shù)的減小，不同展向站位的后緣變彎控制剖面只需同向等比例的偏轉(zhuǎn)即可近似獲取最優(yōu)變形外形。

圖13 配平狀態(tài)多點(diǎn)關(guān)聯(lián)變形矩陣后緣外形變彎控制剖面偏度對(duì)比（Ma = 0.85，H = 11 km）Fig.13 Deflection angles comparison of trailing edge deformation shape control sections in multipoint coupling deformation matrix under trimmed condition（Ma = 0.85，H = 11 km）

巡航任務(wù)剖面最小和最大升力系數(shù)設(shè)計(jì)狀態(tài)對(duì)應(yīng)多點(diǎn)關(guān)聯(lián)變形矩陣后緣外形與后緣無偏轉(zhuǎn)的基本外形在配平狀態(tài)下的內(nèi)、中、外3個(gè)展向站位翼面壓力分布對(duì)比如圖14和圖15所示，其中內(nèi)側(cè)展向站位對(duì)應(yīng)后緣外形變彎控制剖面1，中間展向站位對(duì)應(yīng)變彎控制剖面4，外側(cè)展向站位對(duì)應(yīng)變彎控制剖面6。

圖14 后緣外形變形對(duì)翼面壓力分布影響（Ma = 0.85，H = 11 km，CL = 0.384配平狀態(tài)）Fig.14 Effect of trailing edge deformation on pressure coefficient distribution（Ma = 0.85，H = 11 km，CL = 0.384 trimmed condition）

圖15 后緣外形變形對(duì)翼面壓力分布影響（Ma = 0.85，H = 11 km，CL = 0.576的配平狀態(tài)）Fig.15 Effect of trailing edge deformation on pressure coefficient distribution（Ma = 0.85，H = 11 km，CL = 0.576 trimmed condition）

在小升力系數(shù)設(shè)計(jì)狀態(tài)（CL= 0.384），只有一種后緣變形外形方案，內(nèi)側(cè)展向站位后緣外形變彎控制剖面下偏0.11°，翼面壓力分布沒有明顯變化。中間展向站位上偏0.58°，外側(cè)展向站位上偏0.15°，變彎控制剖面上偏主要影響上翼面的壓力分布變化，具體表現(xiàn)為前緣吸力峰后的低壓區(qū)向后擴(kuò)展，變彎后緣前端所處的當(dāng)?shù)刎?fù)壓強(qiáng)度降低、范圍減小，翼面外側(cè)后緣附近激波強(qiáng)度降低，剖面壓心前移。在大升力系數(shù)設(shè)計(jì)狀態(tài)（CL= 0.576），按照設(shè)計(jì)變量空間選取的差異，多點(diǎn)關(guān)聯(lián)變形矩陣優(yōu)化獲得兩種后緣變形外形方案，在內(nèi)、中、外3個(gè)展向站位，下邊界逐漸縮減方案的后緣變彎控制剖面偏度分別為下偏1.6°、不偏轉(zhuǎn)（0°）和下偏0.92°，上下邊界同時(shí)縮減方案的后緣變彎控制剖面偏度分別為下偏1.1°、0.7°和0.7°，盡管2種外形方案在不同展向站位的后緣變彎控制剖面偏度差別較大，其影響僅局限于變彎后緣型面附近的激波位置變化，翼型前半部分的壓力變化基本一致，整體表現(xiàn)為負(fù)壓強(qiáng)度降低和負(fù)壓范圍的向后擴(kuò)展，剖面壓心后移。由此可知變彎度機(jī)翼減阻的關(guān)鍵在于通過后緣外形變形調(diào)整主翼部分的壓力分布，后緣當(dāng)?shù)匦兔娴膲毫Ψ植甲兓绊戄^小，同時(shí)也說明了對(duì)于同一升力系數(shù)確實(shí)存在多種后緣外形變形方案。

5.3 階梯變形矩陣

結(jié)合變彎度機(jī)翼后緣外形變形巡航任務(wù)剖面減阻收益對(duì)變彎控制剖面偏度不敏感的特點(diǎn)，基于多點(diǎn)關(guān)聯(lián)變形矩陣，從均衡巡航任務(wù)剖面減阻收益角度開展階梯變形矩陣氣動(dòng)設(shè)計(jì)。

選定設(shè)計(jì)變量空間上下邊界同時(shí)縮減的配平狀態(tài)多點(diǎn)關(guān)聯(lián)變形矩陣，綜合矩陣內(nèi)后緣變形外形變彎控制剖面偏轉(zhuǎn)規(guī)律確定偏轉(zhuǎn)基準(zhǔn)，分別以最小和最大升力系數(shù)為邊界，按巡航任務(wù)剖面內(nèi)中間升力系數(shù)與邊界升力系數(shù)的比例關(guān)系，確定與升力系數(shù)線性對(duì)應(yīng)的階梯變形矩陣。典型升力系數(shù)下的階梯變形矩陣與多點(diǎn)關(guān)聯(lián)變形矩陣對(duì)比見圖16，兩者整體偏度基本一致，僅在外側(cè)展向站位略有偏差。

圖16 階梯變形矩陣與多點(diǎn)關(guān)聯(lián)變形矩陣后緣外形變彎控制剖面偏度對(duì)比Fig.16 Deflection angles comparison of trailing edge deformation shape control sections between step deformation matrix and multi-point coupling deformation matrix

階梯變形矩陣內(nèi)后緣變形外形相對(duì)多點(diǎn)關(guān)聯(lián)變形矩陣和無偏轉(zhuǎn)基本后緣的氣動(dòng)特性對(duì)比見表6，表中ΔCD1為階梯變形矩陣相對(duì)基本后緣的阻力系數(shù)減小量，ΔCD2為階梯變形矩陣相對(duì)多點(diǎn)關(guān)聯(lián)變形矩陣的阻力系數(shù)減小量。結(jié)果表明，后緣外形變形調(diào)整后的階梯變形矩陣仍然具有理想的減阻收益，與優(yōu)化得到的多點(diǎn)關(guān)聯(lián)變形矩陣相比，阻力系數(shù)增量不超過0.000 02。

表6 配平狀態(tài)階梯變形矩陣后緣變形外形阻力系數(shù)減小量（Ma = 0.85，H = 11 km）Table 6 Drag coefficient reduction of step deformation matrix trailing edge deformation under trimmed condition（Ma = 0.85，H = 11 km）

值得注意的是，變彎度機(jī)翼后緣外形變形減阻收益隨著設(shè)計(jì)升力系數(shù)的增加而增大，根據(jù)升力系數(shù)等間距劃分設(shè)計(jì)狀態(tài)矩陣確定后緣外形階梯變形矩陣存在減阻收益失衡問題，即矩陣內(nèi)減阻收益小的后緣變形外形過多，而減阻收益大的后緣變形外形過少。為此在大升力系數(shù)區(qū)間內(nèi)針對(duì)性的加密設(shè)計(jì)點(diǎn)來構(gòu)建新的設(shè)計(jì)狀態(tài)矩陣和階梯變形矩陣，新增設(shè)計(jì)點(diǎn)下階梯變形矩陣后緣變形外形阻力系數(shù)減小量對(duì)比見表7，表中ΔCD3

表7 配平狀態(tài)新增設(shè)計(jì)點(diǎn)階梯變形矩陣后緣變形外形阻力系數(shù)減小量（Ma = 0.85，H = 11 km）Table 7 New added design points drag coefficient reduction of step deformation matrix trailing edge deformation under trimmed condition（Ma = 0.85，H = 11 km）

為階梯變形矩陣CL= 0.528對(duì)應(yīng)后緣變形外形相對(duì)基本后緣的阻力系數(shù)減小量，ΔCD4為階梯變形矩陣CL= 0.480對(duì)應(yīng)后緣變形外形相對(duì)基本后緣的阻力系數(shù)減小量，ΔCD5為階梯變形矩陣CL= 0.576對(duì)應(yīng)后緣變形外形相對(duì)基本后緣的阻力系數(shù)減小量。通過數(shù)據(jù)對(duì)比可以看出，在大升力系數(shù)區(qū)間，基于階梯變形矩陣獲得的后緣變形外形減阻量與設(shè)計(jì)升力系數(shù)基本呈線性關(guān)系。在CL= 0.480 ～ 0.576的大升力系數(shù)區(qū)間，相比于加密設(shè)計(jì)狀態(tài)矩陣，不增加設(shè)計(jì)點(diǎn)時(shí)產(chǎn)生的最大阻力系數(shù)增量不小于0.000 11（ΔCD1-ΔCD5，CL= 0.547），增加一個(gè)設(shè)計(jì)點(diǎn)（CL= 0.528）時(shí)，產(chǎn)生的最大阻力系數(shù)增量不小于0.000 04（ΔCD1-ΔCD3，CL= 0.504）。

6 非巡航任務(wù)剖面變形矩陣

考慮航空管制等實(shí)際應(yīng)用場景和相關(guān)設(shè)計(jì)規(guī)范要求，在民機(jī)機(jī)翼精細(xì)化氣動(dòng)設(shè)計(jì)中，大升力系數(shù)時(shí)的抖振特性和飛行馬赫數(shù)增加時(shí)的阻力發(fā)散特性等非巡航任務(wù)剖面性能也需要重點(diǎn)考慮。

某遠(yuǎn)程民機(jī)標(biāo)模巡航馬赫數(shù)Ma= 0.85，巡航高度H= 11 km，設(shè)計(jì)升力系數(shù)CL= 0.48，以此為基準(zhǔn)，選定1.3g過載狀態(tài)為抖振點(diǎn)（Ma=0.85，H= 11 km，CL= 0.624）開展變彎度機(jī)翼后緣外形變形抖振抑制效能分析，選定阻力發(fā)散點(diǎn)（Ma= 0.87，H= 11 km，CL= 0.48）開展減阻收益研究。

6.1 大升力系數(shù)下的抖振抑制變形

根據(jù)單點(diǎn)最優(yōu)變形矩陣和多點(diǎn)關(guān)聯(lián)變形矩陣氣動(dòng)優(yōu)化設(shè)計(jì)方法獲得抖振點(diǎn)的后緣變形外形方案變彎控制剖面偏度對(duì)比如圖17所示。考慮后緣變形外形與CL= 0.576對(duì)應(yīng)方案的相似性，將單點(diǎn)最優(yōu)抖振抑制變形方案視為下邊界逐漸縮減多點(diǎn)關(guān)聯(lián)變形方案?？梢钥闯?，盡管2種抖振抑制方案在中內(nèi)側(cè)的變彎控制剖面偏度差異明顯，但抖振抑制變形規(guī)律基本相同，都是僅在CL= 0.576對(duì)應(yīng)方案的基礎(chǔ)上上偏外側(cè)變彎控制剖面，且增加的偏度也相差不大。

圖17 配平狀態(tài)抖振抑制方案后緣外形變彎控制剖面偏度對(duì)比Fig.17 Deflection angles comparison of buffeting suppression scheme trailing edge deformation shape control sections under trimmed condition

后緣外形抖振抑制變形方案相對(duì)后緣不偏轉(zhuǎn)的基本外形減阻收益對(duì)比如表8所示，阻力系數(shù)減小量超過0.003 5，其中設(shè)計(jì)變量空間更大的單點(diǎn)最優(yōu)變形方案略大。

表8 配平狀態(tài)抖振抑制方案后緣變形外形阻力系數(shù)減小量（Ma = 0.85，H = 11 km）Table 8 Drag coefficient reduction of buffeting suppression scheme trailing edge deformation under trimmed condition（Ma = 0.85，H = 11 km）

后緣外形抖振抑制變形方案機(jī)翼表面壓力分布云圖及極限流線對(duì)比如圖18所示，典型變彎控制剖面翼面壓力分布對(duì)比如圖19所示。總體上2種后緣變形方案的流場相差不大，典型剖面的壓力分布幾乎完全重合，都是通過增加后緣載荷、降低機(jī)翼前段負(fù)壓強(qiáng)度的方法來降低翼面激波強(qiáng)度，進(jìn)而減小阻力。但受制于后緣外形變形能力限制，翼面激波強(qiáng)度仍較大，未能完全消除翼面后緣的流動(dòng)分離。

圖18 抖振點(diǎn)后緣外形變形對(duì)翼面流場分布影響（Ma =0.85，H = 11 km，CL = 0.624配平狀態(tài)）Fig.18 Effect of buffeting point trailing edge deformation on pressure coefficient distribution（Ma = 0.85，H = 11 km，CL = 0.624 trimmed condition）

圖19 抖振點(diǎn)后緣外形變形對(duì)2y/b = 0.68剖面壓力分布影響對(duì)比（Ma = 0.85，H = 11 km，CL =0.62配平狀態(tài)）Fig.19 Effect of buffeting point trailing edge deformation on pressure coefficient distribution at 2y/b=0.68（Ma = 0.85，H = 11 km，CL =0.624 trimmed condition）

6.2 阻力發(fā)散馬赫數(shù)下的減阻變形

根據(jù)單點(diǎn)最優(yōu)變形矩陣設(shè)計(jì)方法獲得非配平及配平條件下的阻力發(fā)散點(diǎn)減阻后緣變形方案變彎控制剖面偏度對(duì)比如圖20所示。在不考慮配平要求時(shí)，后緣變形外形呈內(nèi)側(cè)上偏、外側(cè)下偏的特點(diǎn)，在配平條件下，相對(duì)于Ma= 0.85的巡航狀態(tài)，后緣變形外形同樣呈內(nèi)側(cè)上偏、外側(cè)下偏的規(guī)律，但內(nèi)側(cè)偏度略小。

圖20 阻力發(fā)散減阻方案后緣外形變彎控制剖面偏度對(duì)比Fig.20 Deflection angles comparison of drag divergence scheme trailing edge deformation shape control sections

阻力發(fā)散點(diǎn)后緣變形減阻收益對(duì)比如表9所示，非配平狀態(tài)下的后緣變形相比后緣不變形的基本外形阻力系數(shù)減小0.000 25，配平狀態(tài)下的阻力系數(shù)減小量僅為0.000 11，減阻收益降低超過50%。

非配平及配平狀態(tài)下的阻力發(fā)散點(diǎn)減阻后緣變形方案與基本外形在典型變彎控制剖面的翼面壓力分布對(duì)比分別如圖21和圖22所示。此時(shí)內(nèi)側(cè)及外側(cè)變彎控制剖面上偏及下偏僅能夠影響后緣附近的激波位置，但對(duì)剖面前段的壓力幾乎沒有影響，由此也進(jìn)一步表明后緣外形變形無法有效控制主要由激波阻力決定的阻力發(fā)散。

圖21 阻力發(fā)散點(diǎn)后緣外形變形對(duì)翼面壓力分布影響（Ma = 0.87，H = 11 km，CL = 0.48非配平狀態(tài)）Fig.21 Effect of drag divergence point trailing edge deformation on pressure coefficient distribution（Ma = 0.87，H = 11 km，CL = 0.48 untrimmed condition）

7 結(jié)論

選取基于代理模型的優(yōu)化方法，以某遠(yuǎn)程民機(jī)標(biāo)模為研究對(duì)象，在選定變彎控制剖面參數(shù)化方法的基礎(chǔ)上，開展了考慮變形關(guān)聯(lián)約束的變彎度機(jī)翼后緣外形變形矩陣氣動(dòng)設(shè)計(jì)研究。針對(duì)某遠(yuǎn)程民機(jī)標(biāo)模，得出結(jié)論如下：

1） 對(duì)于參考指關(guān)節(jié)變形結(jié)構(gòu)確定的后緣外形變彎控制剖面，采用三段中弧線旋轉(zhuǎn)角度相同的單變量控制方式盡管變形能力有限，但是完全能夠滿足變彎度機(jī)翼后緣變形外形的氣動(dòng)優(yōu)化設(shè)計(jì)要求。

2） 對(duì)于飛行馬赫數(shù)和高度固定、僅升力系數(shù)隨燃油消耗而逐漸減小的巡航任務(wù)剖面，變彎度機(jī)翼后緣外形變形的減阻收益隨著升力系數(shù)相對(duì)基準(zhǔn)狀態(tài)變化量的增加而增大，且高升力狀態(tài)減阻量約為同等變化幅度對(duì)應(yīng)低升力狀態(tài)的7 ～ 8倍。

3） 在小升力系數(shù)狀態(tài)，后緣外形變形能夠擴(kuò)展前緣吸力峰后的低壓區(qū)，減小變彎后緣當(dāng)?shù)氐呢?fù)壓強(qiáng)度，壓心前移；在大升力系數(shù)狀態(tài)，后緣外形變形能夠降低主翼負(fù)壓強(qiáng)度并向后擴(kuò)展負(fù)壓范圍，壓心后移。

4） 變彎度機(jī)翼減阻的關(guān)鍵在于通過后緣外形變形調(diào)整主翼壓力分布，但對(duì)變彎控制剖面的偏度不敏感，后緣外形變形關(guān)聯(lián)關(guān)系約束不會(huì)明顯減小后緣外形變形的減阻收益。

5） 存在后緣外形變彎控制剖面偏轉(zhuǎn)規(guī)律相同、偏度與升力系數(shù)線性對(duì)應(yīng)，且減阻收益明顯的巡航任務(wù)剖面階梯變形矩陣，相比基于升力等間距變化的階梯劃分，基于減阻量等間距變化的階梯劃分能夠在整個(gè)任務(wù)剖面內(nèi)獲得更大的減阻收益。

6） 大升力區(qū)間內(nèi)的巡航任務(wù)剖面階梯變形矩陣，后緣外形變彎控制剖面偏度、后緣變形外形對(duì)應(yīng)設(shè)計(jì)升力系數(shù)、后緣外形變形減阻量3者之間均呈線性關(guān)系。

7） 對(duì)于由抖振點(diǎn)和阻力發(fā)散點(diǎn)確定的非巡航任務(wù)剖面變形矩陣，變彎度機(jī)翼后緣外形變形相對(duì)巡航任務(wù)剖面變形矩陣同樣具有明顯的規(guī)律性，抖振點(diǎn)的后緣外形變形能夠通過改變翼面載荷分布和激波強(qiáng)度而減小阻力，但不能完全翼面消除分離。阻力發(fā)散點(diǎn)的后緣外形變形對(duì)主翼部分壓力分布影響較小，無法有效控制主要由激波阻力決定的阻力發(fā)散。