翼身連接處幾何參數(shù)優(yōu)化匹配的流動機理仿真

葉 露，吳奕銘，袁廣田

(1.中國民用航空飛行學院，四川 廣漢 618307； 2.中國飛行試驗研究院，陜西 西安 710089)

引言

在飛機生產(chǎn)過程中，需要將飛機的各個部件與機身連接，形成一個整體，如機翼與機身連接、尾翼與機身連接等。飛機的各個部件，如機翼、機身、尾翼等，單獨放在氣流中所產(chǎn)生的阻力總和小于組成一個整體時所產(chǎn)生的阻力，這種飛機各部分之間由于氣流的相互干擾而產(chǎn)生的額外阻力，稱為干擾阻力。干擾阻力產(chǎn)生的主要原因是氣流流經(jīng)各個部件時會出現(xiàn)氣流互相干擾，改變了流場分布，使得流動過早分離，最大升力系數(shù)和臨界迎角減小，嚴重影響飛機大迎角飛行性能，尤其對起飛和著陸階段的安全性和經(jīng)濟性帶來嚴重影響[1]。

由于機翼機身連接處的氣流干擾對全機的氣動性能影響很大，本研究選取某一特定飛機的機翼-機身的連接部位著陸構(gòu)型為研究對象，采用數(shù)值模擬的方法，從連接部位的前緣、展向中段及后段3個方面開展整流參數(shù)設(shè)計研究，為連接部位幾何參數(shù)匹配方案設(shè)計提供理論基礎(chǔ)，選取的研究模型如圖1所示。

圖1 連接部位參數(shù)匹配優(yōu)化的3段

1 連接部位幾何參數(shù)優(yōu)化設(shè)計思想

導致連接部位流動分離的主要原因：首先，受機身影響，使得機翼根部前緣吸力峰增高，逆壓梯度增大，致使后緣流動分離提前；另外，機翼阻滯機身附面層發(fā)展，產(chǎn)生馬蹄渦，對下游流動產(chǎn)生不利影響；再者，機翼與機身大多采用銳角連接，造成附面層摻混堆積，增加附面層流動的不穩(wěn)定性；最后，翼身接合部后段型面局部收縮過快，外加上述3點共同作用，使得機翼后緣形成類似邊角渦的分離流動。根據(jù)流動分離機理，延緩連接部位處的氣流分離，可通過以下方法進行流動控制，如降低機翼翼根前緣吸力峰；減小由于機翼的外形突變而引起的逆壓梯度的大小與范圍，消弱或消除機翼前緣的流動分離；阻止翼根前緣馬蹄渦的形成，從而消除馬蹄渦引起的非定常流動對翼根下游流動穩(wěn)定性的影響，減小馬蹄渦尾跡引起的流動尾跡的不均勻；減小壁面相交處的逆壓梯度，以減少附面層的摻混、堆積，防止氣流過早分離[2-3]。為達到上述效果，可從機翼-機身連接部位的前緣、展向中段及后段3個方面進行參數(shù)優(yōu)化匹配設(shè)計。

2 研究方法

選取某一著陸構(gòu)型下的機身-機翼外形，采用數(shù)值模擬的方法，在原始連接參數(shù)的基礎(chǔ)上，根據(jù)上述設(shè)計思想，進行相關(guān)參數(shù)的設(shè)計，以期獲得最佳的匹配參數(shù)。

2.1 計算模型

選取已有的原始整流外形作為研究基礎(chǔ)，在重新設(shè)計連接部位參數(shù)之前，對帶原始整流的著陸構(gòu)型做氣動性能和流動機理評估，找出原始整流的優(yōu)缺點，以便在重新設(shè)計整流的過程中有所針對。帶原始整流的著陸構(gòu)型和原始整流外形及連接部位局部放大圖如圖2所示。

圖2 原始計算模型示意圖

本研究采用有限體積法數(shù)值求解雷諾平均N-S方程(RANS)，采用能夠高精度預測流動分離的k-ω兩方程湍流模型，時間推進方法為隱式方法[4-5]。

數(shù)值模擬條件為：速度為265 km/h,高度為標準海平面高度，以機翼平均氣動弦為參考長度的雷諾數(shù)Re=1.5×106，氣動力系數(shù)以各自的機翼參考面積為基準給出[1]。

2.2 網(wǎng)格生成

本研究選用六面體網(wǎng)格，基本構(gòu)型計算網(wǎng)格總數(shù)600萬，第一層網(wǎng)格尺寸為cA×10-5(cA為全機的平均氣動弦長)，附面層增長比率為1.16。機翼上表面及翼根處等流動變化劇烈的地方進行了網(wǎng)格加密，以期更好的把握流場的劇烈變化[2]，計算該模型所用網(wǎng)格如圖3所示。

圖3 原始連接參數(shù)著陸構(gòu)型的計算網(wǎng)格模型

3 連接部位幾何參數(shù)優(yōu)化設(shè)計方案

3.1 原始連接參數(shù)著陸構(gòu)型的氣動性能分析

圖4給出了原始連接參數(shù)著陸構(gòu)型下的氣動特性數(shù)據(jù)。數(shù)值模擬與風洞實驗結(jié)果相比，升力系數(shù)CL曲線沿坐標橫軸略向右平移，最大升力系數(shù)相同，臨界迎角略有增大，最大升系數(shù)為2.12，臨界迎角為12°；阻力系數(shù)CD在小迎角上吻合較好；在達到臨界迎角前，升阻比K偏小，最大升阻比同樣出現(xiàn)在迎角4°左右，為8.5，其誤差在允許的范圍內(nèi)，證明所選用的數(shù)值模擬方法可以滿足研究要求。

圖4 原始連接參數(shù)著陸構(gòu)型的氣動特性

著陸構(gòu)型下6°～16°迎角的機翼表面及翼根處流線圖，如圖5所示，當α<6°，模型表面無明顯流動分離；當α>6°，翼根流動分離隨迎角的增加而增強。從機身上的分離渦可以看出，剛開始流動分離僅出現(xiàn)在翼根表面，隨著迎角增大，分離渦在空間的尺寸越來越大，并不斷的沿機翼的展向和弦向擴大[2-3,6]。

圖5 原始連接參數(shù)著陸構(gòu)型機翼表面及翼根處流線分布

3.2 連接參數(shù)優(yōu)化設(shè)計方案1

方案1主要針對翼根前緣部分進行流動控制，減小前緣相對厚度，加大前緣后掠角，展向中段和后段整流改變較小，前緣整流與機翼前緣之間線性過渡，外形如圖6所示。

圖6 方案1外形

方案1數(shù)值模擬計算結(jié)果表明，大迎角下的氣動得到了改善，在16°迎角時，升力系數(shù)由原始整流的2.15提高到2.24。方案1翼根仿真結(jié)果如圖7所示。圖7a為截面壓力分布，橫坐標為翼根翼弦長不同站位X，縱坐標表示壓力系數(shù)Cp。由圖7a可見，其前緣負壓峰值仍然較大，導致失速特性不佳，主要原因是展向中段和后段整流較小，未能從根本上改善翼身連接部位的不光順過渡；中后段整流長度不夠，且整流不飽滿，空間擴張仍然較快；前緣整流沿機翼展向線性過渡，過渡較為劇烈，展向流動分量大[7-10]。由圖7b可知，機身、機翼上流動分離均有較大幅度的改善，但翼根區(qū)域流動仍存在明顯流動氣流分離。

圖7 方案1翼根仿真結(jié)果(α=16°)

3.3 連接參數(shù)優(yōu)化設(shè)計方案2

根據(jù)方案1存在的問題，進行改進得到連接參數(shù)優(yōu)化設(shè)計方案2：擴大展向中段的整流范圍，使翼身夾角增大，翼身過渡更加平緩；中后段整流長度增加，且使整流較飽滿；前緣整流沿機翼展向過渡采用曲線，方案2外形如圖8所示。

圖8 方案2外形

當α=16°，整流方案2全機升力系數(shù)由整流方案1的2.24增加到2.35，臨界迎角超過了16°。方案2翼根仿真結(jié)果如圖9所示，由翼根壓力分布可見，翼根處前緣吸力峰從整流方案1的11降到了9，且翼根流動分離已消弱很多，機身、機翼流動特性良好。

圖9 方案2翼根仿真結(jié)果(α=16°)

3.4 連接參數(shù)優(yōu)化設(shè)計方案3

方案3在設(shè)計方案2的外形基礎(chǔ)上，繼續(xù)加大前緣整流幅度，延伸后緣整流長度，根部相對厚度13.8%，展向?qū)挾葹?.054b/2，b為翼展，方案3外形如圖10所示。

圖10 方案3外形

圖11的計算結(jié)果表明，方案3的最大升力系數(shù)達到了2.39，相對于前2個方案，最大升力系數(shù)均有提升，臨界迎角為16°；最大升阻比為8.8，出現(xiàn)在α=4°時。

圖11 方案3的氣動性能

圖12給出了整流方案3的10°～18°迎角下翼根表面流線圖。α<14°，表面基本為附著流；α≥14°，翼根后緣處出現(xiàn)流動分離趨勢；α>16°，翼根分離發(fā)展較快，翼身接合部出現(xiàn)大范圍分離，空間旋渦特征強烈。

圖12 方案3的翼根處表面流線

綜上，通過本研究連接參數(shù)優(yōu)化匹配，相比于原始方案，對氣動性能的影響主要體現(xiàn)在改善了著陸構(gòu)型的失速特性，消除升力的非線性；其次，改善中小迎角的升阻力特性，升力系數(shù)線性段延長，升阻性能得到改善[11-12]。原始整流的最大升力系數(shù)為2.12，臨界迎角13°；整流方案3的最大升力系數(shù)為2.39，臨界迎角16°，最大升力系數(shù)提高12%以上，臨界迎角增加了3°。

4 結(jié)論

本研究通過數(shù)值模擬的方法，結(jié)合飛機連接部位流動的機理和氣動性能，在某一原始連接外形的基礎(chǔ)上，通過連接部位參數(shù)的優(yōu)化，獲得了最優(yōu)的匹配方案，并得出了從機翼-機身連接部位的前緣、展向中段及后段3個方面整流參數(shù)設(shè)計的原則，具體如下：

(1) 對于前緣整流，主要通過減小翼根處翼型的相對厚度以及使翼身型面平緩過渡來延緩附面層分離，通過增大翼根前緣后掠角來緩和前緣加速特性，降低前緣吸力峰，減弱甚至消除翼根前緣的馬蹄渦，具體參數(shù)的選取，應結(jié)合具體流態(tài)及壓力分布來確定參數(shù)，如果前緣部位的吸力峰較高且逆壓梯度較強，需要盡可能減小前緣半徑和增加前緣后掠角，并適當對前緣進行負扭轉(zhuǎn)；

(2) 對于展向中段的參數(shù)設(shè)計，主要確保與機身連接平緩過渡，具體設(shè)計的長度范圍應根據(jù)連接部位的夾角大小而定，連接部位若為銳角，此時不需要有太大的參數(shù)調(diào)整，若為鈍角，所需的展向整流幅度越大，應盡量使整流后的夾角大于90°；

(3) 對于后段整流，通過改變后段翼根翼型的曲率，減緩流管截面積的擴張，以克服逆壓梯度的作用，這樣不僅可以減小翼身之間的干擾阻力，也可以防止翼根過早出現(xiàn)流動分離。