李東風(fēng), 李立州, 肖豐樂(lè), 李彩霞

(中北大學(xué) 機(jī)電工程學(xué)院, 山西 太原 030051)

0 引言

隨著高速直升機(jī)的發(fā)展,如何降低直升機(jī)氣動(dòng)阻力功率消耗越來(lái)越受到人們的重視。直升機(jī)氣動(dòng)阻力主要由機(jī)身、主旋翼槳轂和起落架產(chǎn)生[1]。研究表明,現(xiàn)代直升機(jī)主旋翼槳轂氣動(dòng)阻力約占總阻力的30%左右[2],減小槳轂的氣動(dòng)阻力是降低直升機(jī)功率消耗的一個(gè)切實(shí)可行的措施。

20世紀(jì)70年代,人們通過(guò)槳轂的迎風(fēng)面積根據(jù)經(jīng)驗(yàn)預(yù)測(cè)槳轂阻力[3];90年代初,NASA采用風(fēng)洞試驗(yàn)預(yù)測(cè)槳轂阻力[4-5]。隨著風(fēng)洞試驗(yàn)成本的不斷增加和數(shù)值分析技術(shù)的發(fā)展,大量研究對(duì)采用計(jì)算流體力學(xué)(Computational Fluid Dynamics,CFD)技術(shù)預(yù)測(cè)槳轂氣動(dòng)阻力的準(zhǔn)確性進(jìn)行了評(píng)估:美國(guó)的貝爾直升機(jī)公司在NRTC/CR項(xiàng)目中采用CFD技術(shù)準(zhǔn)確地預(yù)測(cè)了槳轂的氣動(dòng)阻力[6];德國(guó)和法國(guó)的SHANEL項(xiàng)目檢驗(yàn)了CFD技術(shù)預(yù)測(cè)槳轂氣動(dòng)阻力和偏航力矩的準(zhǔn)確性[7]。目前,國(guó)內(nèi)還沒(méi)有針對(duì)直升機(jī)主旋翼槳轂氣動(dòng)阻力的相關(guān)研究。CFD技術(shù)在槳轂阻力預(yù)測(cè)方面應(yīng)用較晚,其主要原因是:(1)槳轂結(jié)構(gòu)復(fù)雜,難以對(duì)流域劃分高質(zhì)量的結(jié)構(gòu)網(wǎng)格;(2)槳轂各零件之間的距離較近,存在較強(qiáng)的氣流干擾;(3)槳轂為旋轉(zhuǎn)組件且轉(zhuǎn)速較高,所處氣流環(huán)境十分復(fù)雜。

本文采用CFD技術(shù)對(duì)旋翼槳轂氣動(dòng)性能進(jìn)行分析,為設(shè)計(jì)階段提供更為準(zhǔn)確的預(yù)測(cè)槳轂阻力分析方法。

1 槳轂的主要?dú)鈩?dòng)特征

槳轂包含多個(gè)部件且多為鈍體,鈍體繞流必然伴隨著氣流分離。由于空氣具有粘性,因此鈍體繞流存在流動(dòng)分離、流態(tài)轉(zhuǎn)換、旋渦脫落、渦致振動(dòng)等復(fù)雜特性。鈍體按其形狀可分為兩類(lèi),一類(lèi)表面呈光滑曲線,如圓柱、球等,繞流邊界層在物面的分離點(diǎn)不固定,分離位置和繞流特性隨著雷諾數(shù)的變化有較明顯的變化;另一類(lèi)表面有角點(diǎn),如矩形、立方體等,角點(diǎn)成為固定的分離點(diǎn),并可能在物面再附和重新分離,這類(lèi)物體的繞流在較大的雷諾數(shù)范圍內(nèi)變化不明顯[8]。

槳轂的氣動(dòng)阻力分為與氣流分離有關(guān)的壓差阻力和與氣流粘性有關(guān)的摩擦阻力兩部分,其中壓差阻力約占總阻力的90%～95%。壓差阻力是由分離渦在鈍體上脫落引起的[9],其大小與流動(dòng)分離點(diǎn)的位置有關(guān)。若氣流分離時(shí)邊界層內(nèi)為層流流態(tài),則分離點(diǎn)靠前、壓差阻力高,使總阻力大;若氣流分離時(shí)為紊流流態(tài),則對(duì)應(yīng)的流動(dòng)分離點(diǎn)延后,壓差阻力降低,使總阻力減小[10]。

2 槳轂?zāi)Ｐ?/h2>

槳轂?zāi)Ｐ图爸饕考鐖D1所示,計(jì)算域如圖2所示。進(jìn)行槳轂旋轉(zhuǎn)狀態(tài)CFD分析時(shí),所建立的流場(chǎng)域分為靠近槳轂的旋轉(zhuǎn)區(qū)域和外部的靜止區(qū)域,旋轉(zhuǎn)區(qū)域把槳轂全部包裹起來(lái)。由于槳轂結(jié)構(gòu)復(fù)雜,采用非結(jié)構(gòu)網(wǎng)格生成槳轂流域網(wǎng)格[11]。由于靠近槳轂的氣流變化劇烈,因此靠近槳轂的網(wǎng)格劃分較密,網(wǎng)格密度自槳轂表面向計(jì)算域外邊界由密到疏過(guò)渡,解決了網(wǎng)格過(guò)多、計(jì)算時(shí)間過(guò)長(zhǎng)的問(wèn)題,槳轂附近網(wǎng)格分布如圖3所示。

直升機(jī)前飛時(shí)槳轂向前傾斜一定的角度,槳轂旋轉(zhuǎn)角速度為15.13 rad/s,旋轉(zhuǎn)軸為槳轂的中心軸。流場(chǎng)網(wǎng)格模型的主要邊界類(lèi)型如圖2所示,來(lái)流速度為81.39 m/s,出口壓力為標(biāo)準(zhǔn)大氣壓,采用Spalart-Allmaras湍流模型[12]。

在CFD分析時(shí)需要考慮四種工況下槳轂的氣動(dòng)特性:槳轂在非旋轉(zhuǎn)和旋轉(zhuǎn)狀態(tài)下的迎風(fēng)方位,即0°方位和45°方位,迎風(fēng)方位如圖4所示。

圖1 槳轂主要部件Fig.1 Main parts of hub

圖2 計(jì)算域和邊界條件Fig.2 Computational domain and boundary conditions

圖3 槳轂附近網(wǎng)格分布Fig.3 Mesh distribution around hub

圖4 槳轂迎風(fēng)方位Fig.4 Headwind bearing of hub

3 仿真結(jié)果及分析

本文采用Fluent軟件對(duì)槳轂的外流場(chǎng)網(wǎng)格模型進(jìn)行計(jì)算,槳轂表面壓力分布情況如圖5所示。

圖5 槳轂?zāi)Ｐ捅砻鎵毫?chǎng)云圖Fig.5 Surface pressure distributions of hub

由圖可以看出,在槳葉柄和中央件迎風(fēng)前緣發(fā)生較大的氣流分離。在實(shí)際情況中槳葉柄與槳葉相連,槳葉柄前緣不會(huì)出現(xiàn)大規(guī)模的氣流分離,所以氣流分離所產(chǎn)生的壓差阻力忽略不計(jì)。

各工況下槳轂的阻力和相對(duì)于槳轂質(zhì)心的力矩如表1所示。由表1可以看出,與氣流分離有關(guān)的壓差阻力占總阻力的絕大部分。在槳轂處于非旋轉(zhuǎn)狀態(tài)45°方位時(shí)的阻力最大,這與迎風(fēng)面積有關(guān)。槳轂處于45°方位時(shí),槳轂在非旋轉(zhuǎn)和旋轉(zhuǎn)狀態(tài)下的阻力幾乎相等,而在0°方位時(shí)有差別。這是由于槳轂在45°方位角時(shí)的迎風(fēng)面積等于旋轉(zhuǎn)時(shí)的平均迎風(fēng)面積,與文獻(xiàn)[13]的結(jié)論相符;另外,非旋轉(zhuǎn)與旋轉(zhuǎn)狀態(tài)的阻力相差不大,與文獻(xiàn)[7]的風(fēng)洞試驗(yàn)結(jié)果一致。

非旋轉(zhuǎn)和旋轉(zhuǎn)狀態(tài)的速度云圖如圖6所示,表明槳轂的氣動(dòng)阻力與槳轂是否旋轉(zhuǎn)無(wú)關(guān)。

表1 槳轂的阻力和力矩Table 1 Drag and moment of hub

圖6 槳轂附近速度分布云圖Fig.6 Velocity distribution around hub

非旋轉(zhuǎn)和旋轉(zhuǎn)狀態(tài)的槳轂表面壓力分布如圖7所示?？梢钥闯?旋轉(zhuǎn)狀態(tài)的槳轂迎風(fēng)面的中央件、槳葉柄和螺距控制桿壓力一側(cè)升高,另一側(cè)降低。由庫(kù)塔-儒科夫斯基升力定理可知,處于旋轉(zhuǎn)狀態(tài)的槳轂繞流中存在環(huán)量與來(lái)流疊加,槳轂一側(cè)表面流速加快,另一側(cè)表面流速減緩。由伯努利定理可知流速加快的一側(cè)壓力比另一側(cè)小,由此產(chǎn)生壓差,導(dǎo)致繞y軸的側(cè)偏力矩明顯增大,如表1所示。槳轂的側(cè)偏力矩對(duì)直升機(jī)機(jī)身產(chǎn)生反扭矩,一般通過(guò)尾槳產(chǎn)生推力、或拉力通過(guò)直升機(jī)尾梁形成偏轉(zhuǎn)力矩來(lái)抵消槳轂和旋翼的反扭矩。

圖7 槳轂表面壓力分布云圖Fig.7 Surface pressure distributions of hub

4 結(jié)論

本文采用CFD技術(shù)對(duì)直升機(jī)前飛狀態(tài)的槳轂進(jìn)行氣動(dòng)分析,得出以下結(jié)論:

(1)對(duì)流場(chǎng)域劃分非結(jié)構(gòu)網(wǎng)格可以提高網(wǎng)格生成的效率,采用CFD技術(shù)完全可以預(yù)測(cè)槳轂的阻力和力矩。

(2)槳轂氣動(dòng)阻力的主要來(lái)源是由氣流分離引起的壓差阻力,槳轂處于45°方位角阻力大于0°方位角阻力,說(shuō)明槳轂阻力與與槳轂的迎風(fēng)面積有關(guān)。

(3)處于非旋轉(zhuǎn)和旋轉(zhuǎn)狀態(tài)的槳轂阻力幾乎一樣,說(shuō)明阻力與槳轂是否旋轉(zhuǎn)無(wú)關(guān)。

(4)旋轉(zhuǎn)狀態(tài)的槳轂繞y軸轉(zhuǎn)矩比非旋轉(zhuǎn)狀態(tài)明顯增大,說(shuō)明旋轉(zhuǎn)是側(cè)偏力矩的來(lái)源。

參考文獻(xiàn)：

[1] 路錄祥,王新洲,王遇波.直升機(jī)結(jié)構(gòu)與設(shè)計(jì)[M].北京:航空工業(yè)出版社,2009:36-41.

[2] Van Dam C P.Recent experience with different methods of drag prediction[J].Progress in Aerospace Sciences,1999,35(8):751-798.

[3] Sheehy T W,Clark D R.A method for predicting helicopter hub drag[R].ADA021201,1976.

[4] Young L A,Graham D R,Stroub R H.Experimental investigation of rotorcraft hub and shaft fairing drag reduction[J].Journal of Aircraft,1987,24(12):861-867.

[5] Martin D M,Mort R W,Squires P K,et al.Hub and pylon fairing integration for helicopter drag reduction [C]//AHS,Annual Forum 47th.Phoenix,AZ,1991:897-912.

[6] Hill M J,Louis M E.Rotating hub drag prediction methodology[R].American Helicopter Society Specialists’ Conference on Future Vertical Lift Aircraft Design,American Helicopter Society International,Inc,2012.

[7] Costes M,Raddatz J,Borie S,et al.Advanced rotorcraft aeromechanics studies in the French-German SHANEL project[R].35th European Rotorcraft Forum,2009.

[8] 曹廣龍,李旭東,楊德軍.二維不可壓縮粘性流體繞鈍體流動(dòng)的數(shù)值模擬[J].甘肅科技,2012,28(8):57-59.

[9] Gregory J W,Porter C O,McLaughlin T E.Circular cylinder wake control using spatially distributed plasma forcing[R].AIAA-2008-198,2008.

[10] 何川,辛明道.圓柱面上微小后向臺(tái)階層流繞流場(chǎng)的數(shù)值研究[J].工程熱物理學(xué)報(bào),1997,18(5):604-607.

[11] Lee J K, Kwon O J. Predicting aerodynamic rotor-fuselage interactions by using unstructured meshes [J]. Transactions of the Japan Society for Aeronautical and Space Sciences, 2002, 44(146):208-216.

[12] Spalart P R,Allmaras S R.A one-equation turbulence model for aerodynamic flows[R].AIAA-1992-0439,1992.

[13] Sheehy T W.A general review of helicopter rotor hub drag data [J].Journal of the American Helicopter Society,1977,22(2):2-10.