旋翼懸停狀態(tài)槳尖渦測量方法研究

劉平安, 林永峰, 陳垚峰, 袁明川

(中國直升機設(shè)計研究所, 江西 景德鎮(zhèn) 333000)

旋翼懸停狀態(tài)槳尖渦測量方法研究

劉平安*, 林永峰, 陳垚峰, 袁明川

(中國直升機設(shè)計研究所, 江西 景德鎮(zhèn) 333000)

為了研究旋翼槳尖渦渦核結(jié)構(gòu)，采用粒子圖像測速(PIV)技術(shù)對4m直徑槳尖開孔旋翼在懸停狀態(tài)的槳尖渦流場進行了詳細測量，獲得了旋翼145°方位角處槳尖渦附近流場的速度分布，以及槳尖渦的渦核半徑、旋轉(zhuǎn)速度和渦量分布等試驗數(shù)據(jù)，并研究了PIV分析中選取不同審查窗尺寸對槳尖渦結(jié)構(gòu)測量結(jié)果的影響。為了避免槳尖渦位置不穩(wěn)定的影響，提出了基于槳尖渦流場速度特征的條件平均法，并與簡單平均法和基于渦量峰值的條件平均法進行了對比，驗證了對槳尖渦流場進行條件平均的必要性，以及所提出的條件平均法的有效性。

旋翼; 懸停; 槳尖渦; 流場; PIV

0 引 言

在直升機的飛行過程中，旋翼槳葉會脫出強烈集中的槳尖渦，槳尖渦對旋翼流場具有強烈的影響。而且這些脫出的槳尖渦與槳葉會發(fā)生相互干擾，從而在槳葉上產(chǎn)生出強烈的沖擊氣動載荷和噪聲。因此開展旋翼氣動特性和噪聲特性的研究都離不開對槳尖渦的研究。由于旋翼槳尖渦生成和流動的復(fù)雜性，需要開展精細的槳尖渦流場測量，以獲得槳尖渦的特性。

目前，作為非接觸式流場測量方法，激光多普勒測速(LDV)技術(shù)和粒子圖像測速(PIV)技術(shù)在流場測量領(lǐng)域得到了廣泛應(yīng)用。LDV擁有非常小探測體積，最大的優(yōu)勢是其具有很高的空間分辨率，使其可以準(zhǔn)確測量槳尖渦渦核內(nèi)速度梯度特別高的流場[1-2]。國外采用該方法對旋翼流場開展了大量的測量研究。

但是LDV是一種單點測量方法，測量時需要對觀測區(qū)域進行掃描，測量需要大量的時間。而PIV則提供了大區(qū)域的瞬態(tài)速度場測量手段。

而且，受模型運動、自然湍流或流場的不穩(wěn)定性等因素影響，旋翼的流場結(jié)構(gòu)通常是隨機變化的。如動態(tài)失速現(xiàn)象周期出現(xiàn)，但是每個周期的失速期都有獨立的時間歷程。而對于直升機旋翼尾跡，由于模型的揮舞、擺振等運動與旋轉(zhuǎn)頻率不完全相同，將會導(dǎo)致槳尖渦的位置和當(dāng)?shù)貥~的氣動特性在每轉(zhuǎn)都有些許不同。而PIV提供整個區(qū)域的瞬態(tài)速度場，被認為比LDV更加適合用于測量旋翼槳尖渦流場。

隨著相機分辨率的提高，PIV技術(shù)的流場測量空間分辨率也越來越高。通過對測量參數(shù)的優(yōu)化，PIV流場測量的空間分辨率已經(jīng)可以勝任槳尖渦結(jié)構(gòu)的測量。PIV技術(shù)開始被廣泛應(yīng)用于旋翼流場測量[3-7]。

國外的多項先進直升機旋翼研究項目中都采用了PIV技術(shù)開展槳尖渦流場測量(如ERATO[8]項目，ATIC-1[9]和ATIC-2[10]項目以及HART II[11]項目等)，用于分析旋翼流動機理以及噪聲機理。

國內(nèi)近年來也開展了一些采用PIV技術(shù)測量旋翼流場的研究。中國空氣動力研究與發(fā)展中心的楊永東[12]、袁紅剛[13-14]等開展了旋翼流場的測量，初步探索了槳尖渦的測量技術(shù)。黃水林[15]開展了縱列式雙旋翼尾跡特性研究。這些研究主要測量了旋翼的尾流特征以及尾跡邊界，沒有對槳尖渦的結(jié)構(gòu)測量開展詳細研究。

本文采用PIV技術(shù)，測量了4m直徑槳尖開孔旋翼145°渦齡角槳尖渦周圍的流場，并研究了審查窗尺寸，以及平均方法對槳尖渦結(jié)構(gòu)分析結(jié)果的影響規(guī)律。通過研究，掌握了利用PIV測量旋翼槳尖渦渦核結(jié)構(gòu)的方法，獲得了合理的測量結(jié)果。并得到了旋翼槳尖渦附近流場的速度分布，以及槳尖渦的渦核半徑、旋轉(zhuǎn)速度和渦量分布等試驗數(shù)據(jù)。

1 試驗設(shè)備及試驗?zāi)Ｐ?/h2>

1.1 實驗設(shè)備

槳尖流動控制對旋翼槳尖流場的影響試驗在中國直升機設(shè)計研究所懸停試驗廳進行。

使用的旋翼試驗臺為動導(dǎo)數(shù)試驗臺，試驗臺包含支架系統(tǒng)及短艙，短艙包含動力系統(tǒng)、測量系統(tǒng)、操縱控制系統(tǒng)、數(shù)據(jù)采集處理系統(tǒng)、監(jiān)視報警系統(tǒng)及振動監(jiān)視系統(tǒng)等。試驗現(xiàn)場照片如圖1所示。

PIV系統(tǒng)包括：照明激光器、同步控制器和高速數(shù)字相機。

激光系統(tǒng)光源為脈沖式雙Nd∶Yag激光器，脈沖能量為500mJ，脈沖寬度為6～8ns，片光的厚度小于1.5mm，相機的分辨率為2048pixel2048pixel。

粒子擴散裝置采用多噴口方案，噴口間距150mm，粒子播放面積為1.5m×1.5m。當(dāng)試驗?zāi)康暮驮囼灎顟B(tài)不同時，測量區(qū)域所需粒子濃度不同。為了調(diào)節(jié)測量區(qū)域粒子濃度，需要將粒子發(fā)生器噴口布置在不同的地方。為了便于發(fā)煙位置和方向的調(diào)節(jié)，將粒子發(fā)生器放到可移動支撐臺架上。

1.2 試驗?zāi)Ｐ?/p>

試驗?zāi)Ｐ蜑?m直徑級別旋翼，槳轂型式為鉸接式，槳葉片數(shù)2片，槳葉弦長c為148mm，槳轂半徑R為300mm，旋翼實度0.094，有翼型起始位置為0.26R，槳葉負扭轉(zhuǎn)為12°(從旋轉(zhuǎn)中心至槳尖)，槳尖形狀為矩形。在槳尖布置有4根金屬管，從槳尖的前緣引氣，并從槳尖端面排出。4個引氣孔和4個出氣孔分別位于相對槳尖端面前緣點0.15c、0.30c、0.45c和0.60c處。槳葉槳尖噴氣孔布置如圖2所示。

2 試驗方法

本次PIV試驗的測量平面通過槳轂中心且垂直于槳盤平面。槳葉通過時，測量平面平行于槳葉四分之一弦線。激光器位于旋翼斜下方，向旋翼槳尖方向投射。

本次試驗采用灰度互相關(guān)法獲取測量區(qū)域的速度場。在不同的測量狀態(tài)下，根據(jù)測量區(qū)域的流場速度調(diào)節(jié)跨幀時間間隔。粒子在2幀粒子圖像間的最大相對位移小于審查窗尺寸的1/4。

試驗中為了獲得所需方位角，在旋翼臺上安裝了光電碼盤，當(dāng)指定槳葉旋轉(zhuǎn)至90°方位角時提供一個脈沖信號，通過該脈沖信號觸發(fā)PIV測試系統(tǒng)開展測量工作。通過控制脈沖信號與觸發(fā)的時間間隔，可獲得指定的測量方位角。觸發(fā)延時時間t與測量方位角ψ的關(guān)系為：

式中，n為旋翼轉(zhuǎn)速。

圖3為使用長焦鏡頭獲得的槳尖渦流場粒子圖像，圖像中槳尖渦區(qū)域粒子分布比較均勻，粒子清晰可見，可用于槳尖渦流場速度的定量測量。

3 試驗結(jié)果分析處理

本次試驗的目的是使用PIV技術(shù)測量槳尖渦的速度矢量場，并用于獲取槳尖渦的特性，如最大旋轉(zhuǎn)速度，渦核半徑及渦量峰值。對圖3顯示的粒子圖像進行互相關(guān)運算得到了槳尖渦附近的速度矢量場(旋翼轉(zhuǎn)速573r/min，總距5°,拍攝方位角145°)，結(jié)果如圖4所示。由于此時的槳尖渦隨著旋翼下洗流向下運動，并且向旋翼內(nèi)側(cè)收縮，所以經(jīng)過速度疊加后，槳尖渦左上方區(qū)域的速度較高，而右下方區(qū)域速度較低。

圖5給出了圖4中橫線處的y向速度，該橫線通過槳尖渦的渦核中心，因此該曲線為槳尖渦的切向速度分布(使用槳尖速度ΩR無量綱化)。從圖中可以清晰地看到1個正峰值和1個負峰值，這是槳尖渦的渦核邊界，2個峰值之間距離一般被定義為槳尖渦的渦核直徑。

在進一步處理PIV圖像之前，先開展了PIV審查窗尺寸以及平均處理對這些渦特性數(shù)據(jù)的影響。

Fig.5 Blade tip vortex tangential velocity profile obtained from data shown in Fig.4

3.1 空間分辨率對試驗結(jié)果的影響分析

為了能夠精確地測量槳尖渦結(jié)構(gòu)，測試設(shè)備的空間分辨率必須具備足以正確解析槳尖渦渦核區(qū)域極大的速度梯度。在速度梯度測量中衡量探針空間分辨率的參數(shù)為剪切參數(shù)α。該參數(shù)為設(shè)備空間分辨率與剪切層厚度以及速度梯度的比值。具體如下：

式中，Lm為測量設(shè)備的空間分辨率，L為剪切層厚度。

對于槳尖渦，剪切層厚度為渦核半徑rc，而切向速度在渦核半徑內(nèi)從最大值降為零。因此，對于槳尖渦，上式可寫為：

在PIV測量中，測試空間分辨率就是計算速度時所使用審查窗的尺寸。Grant[16]提出，在PIV測量中，審查窗內(nèi)最大速度變化應(yīng)該小于20%。因此，對于槳尖渦渦核流場，為了使用PIV有效測量出最大切向速度，須滿足α<0.2。本文將使用不同的審查窗尺寸對試驗獲取的同一對PIV粒子圖像進行分析處理，開展空間分辨率對測量結(jié)果的影響研究。

審查窗或取樣窗是進行互相關(guān)運算時所取的計算區(qū)域的像素數(shù)，是最重要的PIV參數(shù)之一，相對于其它測試技術(shù)，其定義了PIV的探測體積。

為了獲得旋翼槳尖渦流場結(jié)構(gòu)，使用超長焦鏡頭結(jié)合移測架獲得槳尖渦附近較小范圍的視場，增大了粒子圖像的空間分辨率，測量區(qū)域尺寸為182mm×182mm。

分別選取24pixel×24pixel、32pixel×32pixel、40pixel×40pixel以及48pixel×48pixel這4種審查窗尺寸對圖3顯示的粒子圖像進行互相關(guān)運算。對應(yīng)于拍攝平面，空間分辨率分別達到2.1，2.8，3.6和4.3mm。

表1給出了使用不同審查窗尺寸獲得的渦核參數(shù)。通過對比可以看到，隨著審查窗尺寸增大，測量獲取的槳尖渦的最大旋轉(zhuǎn)速度減小，渦核半徑顯著增大。當(dāng)審查窗尺寸從24pixel×24pixel增大到48pixel×48pixel時，渦核半徑增大了約18%。可見，隨著審查窗尺寸增大，獲得的槳尖渦流場的誤差急劇增大，不再適合槳尖渦渦核結(jié)構(gòu)的測量。

表1 不同審查窗尺寸獲得的渦核參數(shù)對比Table 1 Comparison of vortex characteristics obtained with different interrogation window sizes

審查窗尺寸應(yīng)該越小越好，但是它受到粒子分布、平面內(nèi)粒子對丟失、平面外粒子對丟失、低的圖像亮度或圖像背景噪聲等因素限制。大部分參數(shù)可以根據(jù)流場情況通過調(diào)整激光脈沖間隔時間、片光厚度進行優(yōu)化。

但是還有2個因素仍然限制了窗口尺寸的減小。第一個是渦核內(nèi)部的粒子數(shù)。事實上，懸停狀態(tài)的槳尖渦具有渦核內(nèi)粒子密度低的特性，這是因為槳尖渦的強烈旋轉(zhuǎn)運動所產(chǎn)生的離心力會將粒子趕出渦核中心(圖3中槳尖渦渦核處有一個黑洞就是因為粒子密度低導(dǎo)致的)。第二個因素是粒子在進行互相關(guān)運算的2幀圖像間必須有一定量的相對運動，否則也無法測量出準(zhǔn)確的速度矢量。

通過調(diào)整相機參數(shù)、粒子濃度、激光強度，本次試驗將測量區(qū)域調(diào)整為182mm×182mm。當(dāng)審查窗尺寸取24pixel×24pixel(空間分辨率2.1mm)時，PIV測量系統(tǒng)的剪切參數(shù)α<0.2，滿足PIV測試要求。結(jié)合50%的審查窗重疊率，在槳尖渦渦核內(nèi)可以獲得22個有效的數(shù)據(jù)，可以比較細致地描繪出槳尖渦的渦核結(jié)構(gòu)。

3.2 條件平均

試驗過程中有許多不確定因素，因此，為了避免試驗中出現(xiàn)的隨機誤差，對獲得的試驗數(shù)據(jù)必須進行平均處理。

由于旋翼系統(tǒng)的振動、拍照系統(tǒng)的振動以及槳尖渦脫離的不確定性等原因，導(dǎo)致了每個測量結(jié)果中槳尖渦的位置會有較大差異。如果簡單地把獲得的瞬態(tài)數(shù)據(jù)直接進行平均，會使得槳尖渦不同部位的速度矢量疊加到一起，破壞了槳尖渦的流場結(jié)構(gòu)。本文使用2種條件平均方法處理試驗數(shù)據(jù)，并與簡單平均法進行對比。具體方法如下：

簡單平均法：將相同狀態(tài)、相同視窗下的有效瞬態(tài)處理結(jié)果相同網(wǎng)格點的速度矢量疊加起來，并除以圖像總數(shù)。

條件平均法一：條件平均的主要目的是為了避免槳尖渦位置的不穩(wěn)定性對平均結(jié)果的影響。因此，在進行平均處理前，將有效瞬態(tài)結(jié)果的坐標(biāo)原點平移至渦核中心處，然后再進行平均處理。本文以渦量的峰值所處位置作為槳尖渦的渦核中心。渦量按如下公式計算：

條件平均法二：由于PIV測量結(jié)果是速度矢量場，還可以使用槳尖渦流場中速度特征來定位槳尖渦位置。本文提出，將有效瞬態(tài)結(jié)果的坐標(biāo)原點平移至槳尖渦流場速度最小處，然后再進行平均處理。這樣可以使各個瞬態(tài)結(jié)果中速度大小和方向相似的部分疊加到一起，盡量不破壞槳尖渦結(jié)構(gòu)。

圖6～8分別為使用以上3種平均法處理的槳尖渦剖面流場結(jié)果(為了便于顯示，降低了速度矢量顯示密度，背景為渦量云圖)。對比3幅圖可以看到，使用2種條件平均法處理的試驗結(jié)果幾乎相同，而使用簡單平均法處理出的槳尖渦流場中，渦核部分的渦量明顯比2種條件平均法小。

Fig.7 Velocity vector field of blade tip vortex (conditional average: method one)

Fig.8 Velocity vector field of blade tip vortex (conditional average: method two)

圖9給出了通過3種平均法處理的槳尖渦渦量分布規(guī)律。可以看到，在槳尖渦中心區(qū)域，使用簡單平均法獲得的渦量值比條件平均法小許多，使用簡單平均法獲得的最大渦量比使用條件平均法二小了約44.5%。

圖10給出了通過3種平均法處理的槳尖渦切向速度分布規(guī)律。

圖中，2種條件平均的結(jié)果曲線幾乎重疊在一起，而簡單平均的曲線中，切向速度峰值明顯要小許多，而且2個峰值的間距(渦核直徑)也明顯變大了。在渦核內(nèi)，由簡單平均獲得的切向速度曲線的斜率、渦核內(nèi)的速度梯度變小了，這就是使用簡單平均法渦核內(nèi)渦量會顯著降低的原因。

表2給出了3種平均法得到的渦核半徑以及最大切向速度。簡單平均法獲得的槳尖渦的渦核半徑最大，約比使用條件平均法二獲得的結(jié)果大38%。而其獲得的切向速度峰值最小，約比使用條件平均法二獲得的結(jié)果小10%。2種條件平均法獲得的各項槳尖渦參數(shù)相差不多，條件平均法二結(jié)果略好一些。

表2 不同平均方法獲得的渦核參數(shù)對比Table 2 Comparison of vortex characteristics obtained with different average methods

4 結(jié) 論

通過對試驗結(jié)果的分析，可以得出以下幾點結(jié)論：

(1) 對于本文的試驗工況，當(dāng)測量區(qū)域尺寸為182mm×182mm，審查窗尺寸取24pixel×24pixel時，測量的剪切參數(shù)α<0.2，滿足PIV測試的要求。結(jié)合50%的審查窗重疊率，在槳尖渦渦核內(nèi)可以獲得22個有效的數(shù)據(jù)，可以比較細致地描繪出槳尖渦的渦核結(jié)構(gòu)。

(2) 審查窗尺寸對試驗結(jié)果影響非常大，隨著審查窗尺寸增大，測量獲取的槳尖渦的最大旋轉(zhuǎn)速度減小，渦核半徑顯著增大。當(dāng)審查窗尺寸從24pixel×24pixel增大到48pixel×48pixel時，渦核半徑增大了約18%，測量誤差顯著增大。

(3) 對比分析了3種平均法得到的渦核結(jié)構(gòu)參數(shù)。簡單平均法獲得的槳尖渦的渦核半徑最大，約比使用條件平均法二獲得的結(jié)果大38%。而其獲得的切向速度峰值最小，約比使用條件平均法二獲得的結(jié)果小10%。2種條件平均法獲得的各項槳尖渦參數(shù)相差不多，條件平均法二結(jié)果略好。

[1]Coyne A J, Mahendra J B, Leishman J D. Investigation into the rollup and diffusion of rotor tip vortices using laser Doppler velocimetry[C]. 53th Annual Forum of the American Helicopter Society, Viginia Beach, VA, 1997.

[2]Martin P B, Pugliese J G, Leishman J G. High resolution trailing vortex measurements in the wake of a hovering rotor[C]. 57th Annual Forum of the American Helicopter Society, Washington D C, 2001.

[3]Heineck J T, Yamauchi G K, Woodcock A J, et al. Application of three-component PIV to a hovering rotor wake[C]. 56th Annual Forum of the American Helicopter Society, Virginia Beach, VA, USA, 2000.

[4]Martin P B, Pugliese J G, Leishman J G, et al. Stereoscopic PIV measurements in the wake of a hovering rotor[C]. 56th Annual Forum of the American Helicopter Society, Virginia Beach, USA, 2000.

[5]Richard H, Raffel M. Rotor wake measurements: Full-scale and model tests[C]. 58th Annual Forum of the American Helicopter Society, Montreal, Canada, 2002.

[6]van der Wall B G, Richard H. Analysis methodology for 3C PIV data[C]. 31th European Rotorcraft Forum, Florence, Italy, 2005.

[7]Wadcock1 A J, Yamauchi G K, Solis E, et al. PIV measurements in the wake of a full-scale rotor in forward flight[C]. 29th AIAA Applied Aerodynamics Conference, Honolulu, Hawaii, 2011.

[8]Splettstoesser W R, van der Wall B G, Junker B, et al. The ERATO programme: wind tunnel results and proof of design for an aeroacoustically optimized rotor[C]. 25th European Rotorcraft Forum, Rome, Italy, 1999.

[9]Murashige A, Kobiki N, Tsuchihashi A, et al. ATIC aeroacoustic model rotor test at DNW[C]. 24th European Rotorcraft Forum, Marseilles, France, 1998.

[10]Murashige A, Kobiki N, Tsuchihashi A, et al. Second ATIC aeroacoustic model rotor testat DNW[C]. 26th European Rotorcraft Forum, The Hague, Netherlands, 2000.

[11]Burley C, Brooks T B, van der Wall B, et al. Rotor wake vortex definition-initial evaluation of 3-C PIV results of the hart-II study[C]. 26th European Rotorcraft Forum. Bristol, England, 2002.

[12]楊永東, 武杰. 懸停旋翼槳尖渦的試驗研究[J]. 實驗流體力學(xué), 2008, 22(3): 36-39.

Yang Y D, Wu J. Investigation of hovering rotor tip vortex[J]. Journal of Experiments in Fluid Mechanics, 2008, 22(3): 36-39.

[13]袁紅剛, 李進學(xué), 楊永東, 等. 懸停狀態(tài)下旋翼尾跡測量試驗研究[J]. 空氣動力學(xué)報, 2010, 28(3): 306-309.

Yuan H G, Li J X, Yang Y D, et al. Test investigation of measuring wake geometry for helicopter rotors in hover[J]. Acta Aerodynamica Sinica, 2010, 28(3): 306-309.

[14]袁紅剛, 李進學(xué), 楊永東, 等. 前飛狀態(tài)下旋翼尾跡測量試驗研究[J]. 實驗流體力學(xué), 2010, 24(4): 29-32.

Yuan H G, Li J X, Yang Y D, et al. Test investigation of measuring wake geometry for helicopter rotors in forward flight[J]. Journal of Experiments in Fluid Mechanics, 2010, 24(4): 29-32.

[15]黃水林, 林永峰, 黃建萍, 等. 基于PIV技術(shù)的縱列式雙旋翼尾跡特性實驗研究[J]. 空氣動力學(xué)報, 2012, 30(3): 334-339.

Huang S L, Lin Y F, Huang J P, et al. Expermental investigation on the wake characteristics of tandem twin rotors[J]. Acta Aerodynamicsa Sinica, 2012, 30(3): 334-339.

[16]Grant I. Particle image velocimetry: a review[J]. Proceedings of the Institution of Mechanical Engineers, 1997, 211: 55-76.

(編輯：李金勇)

Blade tip vortex measurements of a hovering rotor

Liu Ping’an*, Lin Yongfeng, Chen Yaofeng, Yuan Mingchuan

(China Helicopter Research and Development Institute, Jingdezhen Jiangxi 333000, China)

In order to investigate the structure of blade tip vortices, Panicle Image Velocimetry (PIV) measurements were conducted on tip vortices trailing from the rotor with slotted blade to quantify the swirl velocity, core radius and vorticity profile at a fixed wake angle of 145 degree. The influence of the PIV image analysis parameters, e.g. sampling window size, on the vortex structure parameters has been investigated. To mitigate the effects of vortex wander, a scheme to align the velocity fields based on the locations where the velocity minimum appears prior to averaging was developed. Compared to simple average, the conditional average method obtains more accurate parameters of the vortex.

rotor;hovering;tip vortex;flow field;PIV

1672-9897(2017)04-0039-06

10.11729/syltlx20160186

2016-12-03;

2017-02-13

LiuPA,LinYF,ChenYF,etal.Bladetipvortexmeasurementsofahoveringrotor.JournalofExperimentsinFluidMechanics, 2017, 31(4): 39-44. 劉平安, 林永峰, 陳垚峰, 等. 旋翼懸停狀態(tài)槳尖渦測量方法研究. 實驗流體力學(xué), 2017, 31(4): 39-44.

V211.71

A

劉平安(1985-)，男，貴州施秉人，工程師。研究方向：直升機空氣動力學(xué)。通信地址：天津空港經(jīng)濟區(qū)東三道中航直升機b座。E-mail: lpa18505097@avic.com

*通信作者 E-mail: lpa18505097@avic.com