黎 森，王定奇，袁長(zhǎng)波，李 密

(中國(guó)飛行試驗(yàn)研究院發(fā)動(dòng)機(jī)所，陜西 西安 710089)

1 引 言

渦輪螺旋槳發(fā)動(dòng)機(jī)因其低速飛行時(shí)具有較高的推進(jìn)效率而廣泛應(yīng)用于執(zhí)行運(yùn)輸、巡邏、滅火等任務(wù)的飛行器上。從節(jié)約運(yùn)營(yíng)成本方面考慮，該類發(fā)動(dòng)機(jī)也受到各類小型客機(jī)、私人飛機(jī)、公務(wù)機(jī)以及農(nóng)業(yè)飛機(jī)生產(chǎn)者的青睞[1]。隨著我國(guó)航空事業(yè)的不斷發(fā)展，新型飛機(jī)對(duì)渦輪螺旋槳發(fā)動(dòng)機(jī)的性能指標(biāo)也提出了更高的要求。該類發(fā)動(dòng)機(jī)動(dòng)力來(lái)自于螺旋槳旋轉(zhuǎn)產(chǎn)生的拉力以及排氣產(chǎn)生的推力，據(jù)統(tǒng)計(jì)，螺旋槳產(chǎn)生的拉力，約占總推力的95%以上。一款新型發(fā)動(dòng)機(jī)的試飛工作貫穿于整個(gè)型號(hào)的設(shè)計(jì)、研制、鑒定、生產(chǎn)和使用全過(guò)程，如果能盡早對(duì)被試發(fā)動(dòng)機(jī)的螺旋槳拉力性能做出有效預(yù)估，就可以及時(shí)地給試飛設(shè)計(jì)提供有效參考，具有重要的工程應(yīng)用價(jià)值。

關(guān)于螺旋槳拉力性能的確定，SAE的報(bào)告給出了具體詳實(shí)的兩種方法，分別為“J”方法和“Theta”方法[2]，這兩種方法都依賴于已知的螺旋槳特性。螺旋槳的風(fēng)洞特性是獲取其拉力的有效方法之一，廣泛應(yīng)用于工程試驗(yàn)中。風(fēng)洞特性通常由風(fēng)洞試驗(yàn)或CFD仿真獲得。螺旋槳風(fēng)洞試驗(yàn)費(fèi)用大、成本高、耗時(shí)長(zhǎng)，CFD仿真則需要大量繪制網(wǎng)格，網(wǎng)格結(jié)構(gòu)復(fù)雜且計(jì)算周期長(zhǎng)。為了節(jié)約時(shí)間和經(jīng)費(fèi)成本，縮短試飛周期，本文提出根據(jù)有限的風(fēng)洞試驗(yàn)點(diǎn)以及通用螺旋槳部件特性來(lái)預(yù)估整個(gè)包線內(nèi)的螺旋槳拉力性能的方法。

2 螺旋槳通用部件特性

表征螺旋槳通用部件特性的參數(shù)通常包括表征來(lái)流相對(duì)速度的參數(shù)——前進(jìn)比J、表征發(fā)動(dòng)機(jī)用于驅(qū)動(dòng)螺旋槳旋轉(zhuǎn)所輸出的機(jī)械能大小的相對(duì)參數(shù)——功率系數(shù)Cp以及表征發(fā)動(dòng)機(jī)輸出功率被用來(lái)產(chǎn)生拉力時(shí)的能量利用率參數(shù)——效率η。具體定義如下：

(1)

(2)

(3)

其中，V、n、D、ρ、T、N分別表示來(lái)流速度、螺旋槳轉(zhuǎn)速、螺旋槳直徑、來(lái)流密度、螺旋槳拉力和螺旋槳輸出功率。通用螺旋槳部件特性會(huì)以等效率線的形式給出，如圖1所示。

(a)螺旋槳部件通用特性(二維)

(b) 螺旋槳部件通用特性(三維)

從圖中可以直觀地觀察到，螺旋槳部件特性圖在三維坐標(biāo)中形似山峰，峰的底面是前進(jìn)比和功率系數(shù)交織成的平面，等效率線即是山峰的等高線，從底向上逐漸增大。山表面上的每一個(gè)點(diǎn)可確定螺旋槳的唯一工作狀態(tài)，山峰的表面包含所有的螺旋槳工作狀態(tài)。根據(jù)部件特性計(jì)算螺旋槳拉力時(shí)，還需要拉力系數(shù)的參與，其定義如下：

(4)

其中，T為螺旋槳拉力。由式(1)～式(4)可以得出：

(5)

在飛行試驗(yàn)中，發(fā)動(dòng)機(jī)功率N可根據(jù)加裝的扭矩壓力傳感器通過(guò)式(6)間接測(cè)得：

N=K·P·n

(6)

其中，K為螺旋槳構(gòu)造常數(shù)，與螺旋槳測(cè)扭機(jī)構(gòu)的構(gòu)造相關(guān)，P為扭矩壓力。根據(jù)式(2)，可進(jìn)一步獲得功率系數(shù)Cp，而前進(jìn)比J可以根據(jù)發(fā)動(dòng)機(jī)狀態(tài)和飛行參數(shù)獲得。通過(guò)三維曲面插值算法，可以在圖1中獲得螺旋槳的工作效率η。最后結(jié)合式(4)和式(5)，可以獲得螺旋槳的拉力。

3 通用特性耦合系數(shù)的確定

由于螺旋槳部件特性通常以無(wú)量綱的形式給出，這種特點(diǎn)為特性耦合算法提供了可能。文獻(xiàn)[3]指出，在缺少部件特性或已知特性點(diǎn)的數(shù)量較少時(shí)，可以使用耦合的方法獲得適用性較好的特性。另外，通過(guò)對(duì)文獻(xiàn)[4]和文獻(xiàn)[5]中給出的螺旋槳特性進(jìn)行整理發(fā)現(xiàn)，螺旋槳的特性在形態(tài)上具有一定的相似性。據(jù)此，可以將螺旋槳的通用部件耦合到實(shí)際螺旋槳特性中，以期通過(guò)有限個(gè)螺旋槳部件的實(shí)際特性，對(duì)整體的特性進(jìn)行有效預(yù)估。

假設(shè)實(shí)際螺旋槳工作狀態(tài)在特性上表示為(J,Cp,η)，其中η是J和Cp的函數(shù)，用公式表示為：

η=f1(J,Cp)

(7)

(8)

通過(guò)試驗(yàn)獲得部分工作狀態(tài)：(Ji,Cpi,ηi),i=1,2,3,…,n，n表示已知試驗(yàn)點(diǎn)數(shù)。

為了對(duì)耦合系數(shù)進(jìn)行優(yōu)化，選擇優(yōu)化目標(biāo)函數(shù)為：

(9)

當(dāng)Of取最小值時(shí)，對(duì)應(yīng)的k1,k2,k3即為最佳耦合因子。

4 基于粒子群算法的最佳耦合因子的確定

粒子群算法起源于自然界中鳥(niǎo)群覓食行為，由美國(guó)工程師Eberhart最早提出。鳥(niǎo)群在覓食過(guò)程中，會(huì)“自發(fā)”集中在食物最密集的地區(qū)，相當(dāng)于優(yōu)化過(guò)程中的最優(yōu)解。然而這種“自發(fā)”的意識(shí)基于鳥(niǎo)的兩種能力：一是鳥(niǎo)自己的經(jīng)驗(yàn)，它會(huì)根據(jù)自己經(jīng)過(guò)的位置積累一些信息，知道哪個(gè)位置的食物最多，這個(gè)食物最多的位置相當(dāng)于優(yōu)化問(wèn)題中的局部最優(yōu)解；二是鳥(niǎo)會(huì)借鑒其它鳥(niǎo)的經(jīng)驗(yàn)，從其它鳥(niǎo)那里獲得食物最多位置信息并逐步向食物最多的地方逼近，這個(gè)食物最多的位置即是全局最優(yōu)解。在向食物最多的位置逼近過(guò)程中，鳥(niǎo)會(huì)在每個(gè)停駐點(diǎn)的位置信隨機(jī)對(duì)四周的一定范圍進(jìn)行隨機(jī)搜索，以期待這種隨機(jī)搜索可以獲得食物最多的位置信息，并對(duì)以前的全局最優(yōu)的位置進(jìn)行更正。

最佳耦合因子的確定過(guò)程是一個(gè)三維度的全局優(yōu)化問(wèn)題，第i個(gè)粒子位置可以表示為：

其中，m表示粒子數(shù)。將每個(gè)粒子代入式(10)計(jì)算出每個(gè)粒子的適應(yīng)值，記錄每個(gè)粒子經(jīng)歷的最優(yōu)位置，即局部最優(yōu)解：

整個(gè)群體所有粒子經(jīng)歷的最優(yōu)位置記錄為：

粒子的運(yùn)動(dòng)速度記為：

位置更新如下：

5 算例與仿真結(jié)果

用粒子群算法優(yōu)化耦合系數(shù)時(shí)，輸入初始條件包括各個(gè)維度的取值范圍、粒子群規(guī)模、慣性因子、速度限制值及收斂條件。求解的收斂過(guò)程如圖3所示。

圖3 目標(biāo)函數(shù)值隨迭代次數(shù)收斂過(guò)程

由收斂過(guò)程可知，經(jīng)過(guò)71次迭代，目標(biāo)函數(shù)最終收斂至

0.0513。取收斂點(diǎn)的粒子位置參數(shù)，可獲取耦合后的螺旋槳特性圖，最終按照式(1)～式(6)計(jì)算飛行過(guò)程中的拉力，并與實(shí)際螺旋槳特性進(jìn)行對(duì)比，計(jì)算結(jié)果如圖4～圖6所示。

圖4 4km穩(wěn)定平飛時(shí)不同功率狀態(tài)下發(fā)動(dòng)機(jī)拉力曲線對(duì)比

圖5 6km穩(wěn)定平飛時(shí)不同功率狀態(tài)下發(fā)動(dòng)機(jī)拉力曲線對(duì)比

圖6 8km穩(wěn)定平飛時(shí)不同功率狀態(tài)下發(fā)動(dòng)機(jī)拉力曲線對(duì)比

從仿真結(jié)果看，粒子群算法可以有效地計(jì)算通用螺旋槳部件特性的耦合系數(shù)，且通過(guò)耦合特性計(jì)算的螺旋槳拉力和根據(jù)風(fēng)洞試驗(yàn)特性計(jì)算的飛行拉力之間誤差約為2%～3%。

6 結(jié) 論

計(jì)算螺旋槳發(fā)動(dòng)機(jī)飛行中的拉力通常使用螺旋槳特性，獲取有效的螺旋槳特性需要大量的風(fēng)洞試驗(yàn)或CFD計(jì)算，周期長(zhǎng)，投入大。本文通過(guò)風(fēng)洞試驗(yàn)數(shù)據(jù)驗(yàn)證了螺旋槳通用特性耦合算法的有效性，結(jié)果表明，粒子群算法在耦合系數(shù)確定時(shí)非常有效，且使用耦合后特性計(jì)算螺旋槳發(fā)動(dòng)機(jī)飛行拉力具有一定的精度，該方法具有工程應(yīng)用價(jià)值。