后緣小翼智能旋翼減振效果影響因素分析

劉士明, 楊衛(wèi)東, 虞志浩, 吳 杰

(1.南京航空航天大學(xué) 航空宇航學(xué)院 直升機(jī)旋翼動(dòng)力學(xué)國(guó)家級(jí)重點(diǎn)實(shí)驗(yàn)室，南京 210000;2.江蘇科技大學(xué) 船舶與海洋工程學(xué)院，鎮(zhèn)江 212000)

后緣小翼智能旋翼減振效果影響因素分析

劉士明1, 楊衛(wèi)東1, 虞志浩1, 吳 杰2

(1.南京航空航天大學(xué) 航空宇航學(xué)院 直升機(jī)旋翼動(dòng)力學(xué)國(guó)家級(jí)重點(diǎn)實(shí)驗(yàn)室，南京 210000;2.江蘇科技大學(xué) 船舶與海洋工程學(xué)院，鎮(zhèn)江 212000)

建立了帶后緣小翼智能旋翼氣動(dòng)彈性載荷計(jì)算模型及減振優(yōu)化分析方法。模型考慮剛體后緣小翼的氣動(dòng)力與慣性力對(duì)彈性槳葉系統(tǒng)的影響，使用粘性渦粒子法結(jié)合翼型查表法計(jì)算旋翼氣動(dòng)載荷，采用力積分法計(jì)算槳葉與槳轂載荷，構(gòu)造了包含槳葉根部扭轉(zhuǎn)及槳轂振動(dòng)載荷為目標(biāo)函數(shù)的優(yōu)化問(wèn)題，基于最速下降-黃金分割組合優(yōu)化算法尋找最佳小翼偏轉(zhuǎn)規(guī)律。研究發(fā)現(xiàn)，建立的后緣小翼載荷控制方法有效，可降低振動(dòng)目標(biāo)函數(shù)70%。槳葉的彈性扭轉(zhuǎn)使后緣小翼能有效實(shí)施減振，但彈性扭轉(zhuǎn)對(duì)小翼氣動(dòng)力矩的放大作用使減振時(shí)通常伴隨著槳葉扭轉(zhuǎn)載荷增大的現(xiàn)象。

直升機(jī)；旋翼；振動(dòng)；后緣小翼；優(yōu)化方法

旋翼是直升機(jī)最重要的部件之一，既是直升機(jī)的升力面又是操縱面。前飛時(shí)旋翼工作環(huán)境復(fù)雜，由于氣流不對(duì)稱(chēng)、周期變距、槳渦干擾等原因，槳葉上的周期交變載荷傳到機(jī)身上給直升機(jī)帶來(lái)嚴(yán)重的振動(dòng)問(wèn)題，槳轂交變載荷是直升機(jī)機(jī)身振動(dòng)的源頭。主動(dòng)控制后緣小翼(Actively Controlled Flap, ACF)智能旋翼通過(guò)控制安裝在槳葉外端后緣處的小翼按合理規(guī)律偏轉(zhuǎn)，改變氣動(dòng)力分布及槳葉彈性響應(yīng)，達(dá)到降低直升機(jī)振動(dòng)水平的目的，是一種重要的直升機(jī)振動(dòng)主動(dòng)控制方法。90年代初壓電材料等智能材料迅速發(fā)展，使這種直升機(jī)振動(dòng)主動(dòng)控制方法成為可能。

目前，國(guó)際上對(duì)后緣小翼智能旋翼的振動(dòng)與噪聲控制研究已全面展開(kāi)[1]，以理論、風(fēng)洞試驗(yàn)、飛行測(cè)試等方式進(jìn)行，其中理論研究是設(shè)計(jì)新構(gòu)型和試驗(yàn)研究的前提。CHOPRA等[2-3]分別建立了后緣小翼智能旋翼氣動(dòng)彈性分析模型，論證了后緣小翼在振動(dòng)控制方面的作用。后緣小翼減振研究多圍繞小翼對(duì)4片槳葉旋翼的4/rev槳轂載荷影響展開(kāi)，而對(duì)于某些狀態(tài)的特定方向，8/rev諧波載荷也不容忽視。從DATTA等[4-5]進(jìn)行旋翼載荷分析的計(jì)算結(jié)果中可以發(fā)現(xiàn)，改變0～2/rev的氣動(dòng)力輸入會(huì)影響3～10/rev的氣動(dòng)力計(jì)算結(jié)果，說(shuō)明低階氣動(dòng)力確實(shí)會(huì)對(duì)旋翼系統(tǒng)的高階響應(yīng)產(chǎn)生影響，使用低階小翼輸入控制高階載荷值得嘗試。在研究時(shí)還發(fā)現(xiàn)，使用后緣小翼進(jìn)行減振時(shí)可能增加低階槳葉載荷，槳葉載荷的低階諧波成分是旋翼疲勞問(wèn)題的主要根源。同時(shí)，后緣小翼的偏轉(zhuǎn)對(duì)槳尖響應(yīng)與根部扭轉(zhuǎn)力矩有很大影響[6]，而根部扭轉(zhuǎn)力矩與操縱拉桿載荷密切相關(guān)，使用后緣小翼進(jìn)行旋翼減振時(shí)可能帶來(lái)的額外根部扭轉(zhuǎn)力矩不容忽視。研究小翼減振機(jī)理、尋找根部扭矩增加的根源對(duì)指導(dǎo)工程設(shè)計(jì)很有價(jià)值。

本文推導(dǎo)旋翼動(dòng)力學(xué)方程，建立后緣小翼智能旋翼氣動(dòng)彈性綜合分析模型，借助最速下降-黃金分割組合優(yōu)化算法研究后緣小翼對(duì)智能旋翼的載荷控制效果，并以構(gòu)造加權(quán)目標(biāo)函數(shù)、假設(shè)無(wú)小翼氣動(dòng)力矩、假設(shè)結(jié)構(gòu)扭轉(zhuǎn)剛硬等方法，討論分析槳轂載荷控制和槳葉扭轉(zhuǎn)特性、根部扭轉(zhuǎn)力矩間的關(guān)系。

1 智能旋翼氣彈模型

1.1 旋翼氣動(dòng)彈性方程

考慮剛性后緣小翼慣性力和氣動(dòng)力對(duì)彈性槳葉旋翼系統(tǒng)的影響，基于Hamilton原理構(gòu)建后緣小翼智能旋翼系統(tǒng)動(dòng)力學(xué)方程

(1)

式中:δTf、δTb、δUb和δWa分別為小翼動(dòng)能、槳葉動(dòng)能、槳葉應(yīng)變能和旋翼氣動(dòng)力做功的變分。

為描述槳葉的彈性運(yùn)動(dòng)及小翼的剛體運(yùn)動(dòng)，建立慣性坐標(biāo)系I，槳轂旋轉(zhuǎn)坐標(biāo)系H，槳葉變形前坐標(biāo)系B，變形后剖面坐標(biāo)系E，以及小翼坐標(biāo)系C，如圖1所示。C系原點(diǎn)位于小翼轉(zhuǎn)軸處，C2正方向朝小翼前緣，C1方向同E1，C3方向由右手定則確定。

槳葉上質(zhì)點(diǎn)A經(jīng)剛體運(yùn)動(dòng)和彈性變形后，在槳轂旋轉(zhuǎn)坐標(biāo)系H中的坐標(biāo)為

(2)

式中：Re是變形后A所在參考剖面坐標(biāo)原點(diǎn)在槳轂旋轉(zhuǎn)坐標(biāo)系中的坐標(biāo)，CBE是E坐標(biāo)系到B坐標(biāo)系的坐標(biāo)轉(zhuǎn)換矩陣，(0,η,ζ)T是A點(diǎn)在變形后坐標(biāo)系中的位置。

小翼的運(yùn)動(dòng)預(yù)先給定，不作為自由度，與槳葉上A點(diǎn)處于同一參考剖面上的小翼上點(diǎn)B在槳轂旋轉(zhuǎn)坐標(biāo)系H中的變形后坐標(biāo)為

(3)

式中：(0,ηc,ζc)T是B點(diǎn)在小翼坐標(biāo)系中的坐標(biāo)，XH是小翼轉(zhuǎn)軸到槳葉變距軸線的距離。

圖1 旋翼系統(tǒng)坐標(biāo)系介紹Fig.1 Coordinate of rotor

點(diǎn)B相對(duì)于慣性系在槳轂系中的加速度為

(4)

式中：Ω為槳轂旋轉(zhuǎn)角速度矢量。

根據(jù)d’Alembert原理，以慣性力虛功的形式表示后緣小翼的動(dòng)能變分，并從中分離出質(zhì)量矩陣Mf和非線性廣義力向量FMf：

(5)

用相同方法得到槳葉的動(dòng)能變分

(6)

在方程建立時(shí)將帶小翼的有限元單元中小翼的廣義質(zhì)量矩陣和非線性廣義力向量加到對(duì)應(yīng)位置。

基于幾何精確大變形梁模型推導(dǎo)槳葉應(yīng)變能[7]：

δUb=?(δεTQε)dAdx

(7)

氣動(dòng)力對(duì)槳葉做功由虛位移和外載荷點(diǎn)乘得到

(8)

式中：Gb是廣義力系數(shù)矩陣,Fa和Ma是變形后坐標(biāo)系中的氣動(dòng)力與氣動(dòng)力矩，qb是廣義自由度，F(xiàn)A是氣動(dòng)力虛功引起的廣義力。

使用CFD方法制作翼型數(shù)據(jù)表，氣彈迭代時(shí)由馬赫數(shù)、翼型迎角、小翼偏角三個(gè)量插值計(jì)算帶小翼翼型氣動(dòng)力的環(huán)量部分，包含翼型迎角與小翼偏角的速度與加速度影響的氣動(dòng)力非環(huán)量部分由Theodorsen模型[8]計(jì)算。槳盤(pán)平面的誘導(dǎo)速度分布[9]由黏性渦粒子法計(jì)算，位置x處誘導(dǎo)速度vi為

(9)

式中：N是粒子數(shù)目，k是光滑參數(shù)，x是空間位置，K(ρ)是Biot-Savart核函數(shù)，α是粒子的渦矢量。

1.2 方程求解和載荷計(jì)算

將式(5)～(8)代入式(1)中得到旋翼動(dòng)力學(xué)微分方程，由于高頻彈性扭轉(zhuǎn)及離心力的剛化作用，此方程通常是剛性的，本文采用隱式梯形公式與牛頓迭代算法在位形空間中求解方程[10。求解氣彈耦合動(dòng)力學(xué)方程得到槳葉的彈性響應(yīng)，使用力積分法計(jì)算槳葉結(jié)構(gòu)載荷，從待求載荷的徑向位置到槳尖的剖面載荷積分得到參考點(diǎn)處的結(jié)構(gòu)載荷。結(jié)構(gòu)載荷、慣性力、外力形成平衡力系，因此得到剖面慣性力和氣動(dòng)力后，剖面結(jié)構(gòu)載荷可寫(xiě)成

(10)

面力在剖面坐標(biāo)系中的坐標(biāo)。

(11)

式中：Nb為槳葉片數(shù)，F(xiàn)為載荷向量，帶預(yù)錐角βp的旋翼，t為不同槳葉方位角ψ處的坐標(biāo)變換矩陣。

2 優(yōu)化問(wèn)題建立與求解

為獲取槳轂振動(dòng)載荷控制所需的最佳小翼偏角，建立如下優(yōu)化問(wèn)題

設(shè)計(jì)變量：以2～5/rev小翼操縱的傅里葉級(jí)數(shù)作為設(shè)計(jì)變量，記為δf

δf=[δ2c,δ2s,δ3c,δ3s,δ4c,δ4s,δ5c,δ5s]T

(12)

在一個(gè)槳葉旋轉(zhuǎn)周期內(nèi)，小翼偏角可表示為

(13)

目標(biāo)函數(shù)：為滿(mǎn)足研究需求，本文研究中共用到以下三種目標(biāo)函數(shù)：

(1)J4：4/rev槳轂振動(dòng)控制

以計(jì)算穩(wěn)態(tài)后的無(wú)量綱槳轂振動(dòng)載荷諧波幅值建立振動(dòng)目標(biāo)函數(shù)。對(duì)四片槳葉旋翼而言，4/rev槳轂載荷是機(jī)身振動(dòng)的主要來(lái)源，僅考慮4/rev槳轂載荷時(shí)目標(biāo)函數(shù)J4可表示為

(14)

(2)J48：4/rev與8/rev槳轂振動(dòng)控制

對(duì)于某些狀態(tài)，8/rev載荷幅值也比較明顯。可在優(yōu)化目標(biāo)中同時(shí)考慮4/rev與8/rev槳轂載荷，構(gòu)建新的振動(dòng)目標(biāo)函數(shù)J48，表達(dá)式為

(15)

(3)J48T：槳轂振動(dòng)與根部扭矩加權(quán)控制

研究減振與根部扭矩的關(guān)系，構(gòu)建扭矩目標(biāo)函數(shù)

(16)

構(gòu)建的槳轂振動(dòng)與根部扭矩加權(quán)目標(biāo)函數(shù)為

(17)

優(yōu)化問(wèn)題求解時(shí)，使用基于梯度的最速下降法確定尋優(yōu)方向，使用黃金分割法進(jìn)行一維尋優(yōu)。參照?qǐng)D2，最速下降-黃金分割組合優(yōu)化算法的計(jì)算流程為：

①初始化n=0。

(18)

置k=1。

③根據(jù)梯度的物理意義，沿梯度的反方向標(biāo)量場(chǎng)下降最快，確定小翼操縱的更新方向

(19)

(20)

⑧令k=k+1，轉(zhuǎn)⑤。

圖2 優(yōu)化問(wèn)題求解流程圖Fig.2 Flowchart of optimization problem

3 模型與驗(yàn)證

3.1 計(jì)算模型與假設(shè)

以全尺寸BO105旋翼為理論計(jì)算的基準(zhǔn)旋翼，基本參數(shù)如下表，槳葉結(jié)構(gòu)參數(shù)見(jiàn)文獻(xiàn)[13]。單片小翼重0.2 kg，小翼質(zhì)心位于小翼轉(zhuǎn)軸后方0.011 m處，小翼位置與尺寸如圖3。為簡(jiǎn)化分析模型，提出以下假設(shè)：模仿壓電堆、金屬框等對(duì)槳葉質(zhì)量特性的影響，帶小翼的槳葉段在原質(zhì)心處增加10%槳葉質(zhì)量的額外質(zhì)量；忽略小翼對(duì)槳葉剖面剛度、慣量、重心位置等的影響；忽略高階小翼操縱對(duì)旋翼配平的影響。

表1 模型旋翼結(jié)構(gòu)參數(shù)

圖3 后緣小翼位置與尺寸Fig.3 Position and dimension of trailing edge flap

3.2 分析模型驗(yàn)證

以帶0.2c小翼的NACA23012翼型壓力系數(shù)分布驗(yàn)證本文計(jì)算翼型氣動(dòng)力的能力，計(jì)算狀態(tài)為：α=8°，δ=15°,Ma=0.105,Re=2.19×106。從下圖可見(jiàn)，計(jì)算值與試驗(yàn)值[14]吻合很好，本文使用的翼型氣動(dòng)力模型足夠精確，可用于帶后緣小翼的旋翼氣動(dòng)彈性分析。

圖4 帶小翼翼型表面壓力系數(shù)分布驗(yàn)證Fig.4 Cp distribution of airfoil with trailing edge flap

以BO105旋翼風(fēng)洞實(shí)驗(yàn)結(jié)果和CAMRAD/JA計(jì)算結(jié)果[15]驗(yàn)證本旋翼模型計(jì)算旋翼載荷的能力，計(jì)算狀態(tài)為μ=0.197，αs=-4.8°，CT/σ=0.071。從圖5可以看出，本文建立的計(jì)算模型能夠很好地預(yù)測(cè)槳葉揮舞彎矩，捕捉到了0.57R處旋翼后行邊揮舞彎矩的變化情況?？傮w上計(jì)算精度優(yōu)于使用自由尾跡/準(zhǔn)定常氣動(dòng)模型計(jì)算的CAMRAD/JA。

圖5 旋翼載荷驗(yàn)證, r/R=0.57Fig.5 Validation of rotor load, r/R=0.57

4 結(jié)果與討論

4.1 雙頻控制結(jié)果

目標(biāo)函數(shù)僅考慮4/rev槳轂載荷時(shí)，有效降低了4/rev載荷幅值，但對(duì)于某些計(jì)算狀態(tài)，槳轂垂向8/rev載荷幅值較大，在進(jìn)行振動(dòng)控制時(shí)會(huì)出現(xiàn)8/rev槳轂載荷幅值增大的情況，如圖6(b)。在目標(biāo)函數(shù)中考慮8/rev槳轂載荷的影響，可在幾乎不增加小翼操縱幅值的情況下有效降低目標(biāo)載荷，小翼偏轉(zhuǎn)規(guī)律如圖7。由于“葉珊效應(yīng)”[16]的存在，低階的氣流擾動(dòng)會(huì)引起對(duì)應(yīng)階次的槳葉振動(dòng)，而槳葉振動(dòng)會(huì)引起更高階的氣動(dòng)載荷，由于這種氣動(dòng)力與槳葉彈性間復(fù)雜耦合的存在，2～5/rev小翼操縱不僅對(duì)4/rev槳轂載荷控制有效，對(duì)8/rev槳轂載荷控制同樣有效，能實(shí)現(xiàn)低階小翼輸入對(duì)高階載荷的控制，在4/rev與8/rev兩個(gè)頻率同時(shí)減振時(shí)，槳轂振動(dòng)目標(biāo)函數(shù)降低了70%左右。

圖6 無(wú)量綱槳轂載荷控制比較Fig．6Dimensionlesshubloadscontroleffect圖7小翼偏轉(zhuǎn)規(guī)律Fig．7Trailingedgeflapdeflectionangle

偏角為正時(shí)小翼后緣向下偏，如90°方位角附近，小翼下偏后改變了剖面等效弦長(zhǎng)與來(lái)流的夾角，有效迎角增大，氣動(dòng)力系數(shù)增大，參考圖8(c)。同時(shí)，小翼下偏后在翼型下表面后緣處形成較強(qiáng)的高壓區(qū)，如圖4，此處壓差到變距軸線位置的力臂較大，形成明顯的低頭力矩，如圖8(d)，俯仰力矩的改變直接影響槳葉的彈性扭轉(zhuǎn)分布。受控小翼偏轉(zhuǎn)時(shí)改變翼型附近流場(chǎng)分布，帶小翼段的槳葉直接發(fā)生氣動(dòng)力和氣動(dòng)力矩的變化，小翼產(chǎn)生的氣動(dòng)力使槳葉的整體揮舞和扭轉(zhuǎn)響應(yīng)發(fā)生變化，從而影響旋翼整體的氣動(dòng)力分布，如圖8(a)、(b)、(e)、(f)。

小翼按最優(yōu)偏轉(zhuǎn)規(guī)律運(yùn)動(dòng)時(shí)氣動(dòng)力作用下的槳轂載荷幅值最小，但此時(shí)通常伴隨著根部扭轉(zhuǎn)力矩增大的現(xiàn)象。如下圖，小翼受2～5/rev控制時(shí)對(duì)應(yīng)階次的槳葉根部扭轉(zhuǎn)力矩幅值均高于基準(zhǔn)槳葉，載荷控制時(shí)帶來(lái)了額外的低頻振動(dòng)。

4.2 扭轉(zhuǎn)問(wèn)題討論

以J48T作為此處研究的目標(biāo)函數(shù)，改變扭矩權(quán)重系數(shù)wT，比較槳轂載荷控制與槳葉根部扭轉(zhuǎn)力矩的關(guān)系。從下圖可以看出，wT越小，槳轂振動(dòng)控制的效果越好，最佳可減振70%左右(圖中A點(diǎn)處)，但同時(shí)根部扭轉(zhuǎn)力矩會(huì)增大。在保證根部扭轉(zhuǎn)力矩不增加的情況下進(jìn)行載荷控制，只能減振不到40%(圖中B點(diǎn)處)，僅為最佳減振效果的57%。為使后緣小翼充分發(fā)揮減振效果，有必要在根部扭轉(zhuǎn)力矩方面做一些讓步。

使用小翼進(jìn)行槳轂載荷控制時(shí)槳葉根部扭轉(zhuǎn)載荷增大，部分源于小翼偏轉(zhuǎn)產(chǎn)生的附加俯仰力矩，而小翼引起的槳葉彈性扭轉(zhuǎn)響應(yīng)改變后，槳葉整體的氣動(dòng)力分布變化又對(duì)這一現(xiàn)象起了一定的放大作用。為進(jìn)一步分析小翼減振與槳葉扭轉(zhuǎn)間的關(guān)系，對(duì)兩種假設(shè)的極限情況進(jìn)行分析：

(1)增大扭轉(zhuǎn)剛度，分析小翼對(duì)扭轉(zhuǎn)剛硬槳葉的減振效果；

(2)計(jì)算時(shí)人工消除掉小翼偏轉(zhuǎn)引起的翼型俯仰力矩變化，分析一種理想的“無(wú)力矩后緣小翼”對(duì)旋翼的減振效果。

如圖11，實(shí)線對(duì)應(yīng)基準(zhǔn)槳葉的槳尖扭轉(zhuǎn)響應(yīng)，小翼受控時(shí)扭轉(zhuǎn)響應(yīng)發(fā)生變化，扭轉(zhuǎn)幅值增大?！盁o(wú)力矩小翼”受控時(shí)僅依靠小翼引起的升力系數(shù)和阻力系數(shù)變化影響旋翼系統(tǒng)，使扭轉(zhuǎn)響應(yīng)發(fā)生的變化很小。當(dāng)槳葉扭轉(zhuǎn)剛硬時(shí)，槳葉不發(fā)生扭轉(zhuǎn)。

扭轉(zhuǎn)剛硬槳葉和“無(wú)力矩小翼”都難以通過(guò)槳葉的扭轉(zhuǎn)改變迎角分布，因而很難對(duì)小翼的效果起到放大的作用，所以這兩種情況的載荷控制效果都不理想，如圖12?？梢?jiàn)，后緣小翼智能旋翼在載荷控制時(shí)槳葉的柔性和小翼偏轉(zhuǎn)引起的氣動(dòng)俯仰力矩缺一不可。

(a)0．68R升力分布(c)0．75R升力分布(e)0．82R升力分布

(b)0．68R俯仰力矩(d)0．75R俯仰力矩(f)0．82R俯仰力矩圖8 振蕩氣動(dòng)力分布Fig．8Vibratoryaerodynamicforcedistribution

圖9 控制前后槳葉根部扭轉(zhuǎn)力矩比較Fig．9Comparisonoftorsionalmomentatbladeroot圖10 減振與扭轉(zhuǎn)Pareto曲線Fig．10Paretocurveofvibrationcontrolandtorsionmoment圖11 槳尖扭轉(zhuǎn)響應(yīng)Fig．11Bladetiptorsionresponse

圖12 無(wú)量綱4/rev槳轂載荷Fig.12 Dimensionless 4/rev hub loads

5 結(jié) 論

(1)由2～5/rev后緣小翼操縱不僅能對(duì)4/rev槳轂載荷實(shí)行控制，對(duì)8/rev的高階載荷控制同樣有效。在4、8/rev雙頻率同時(shí)控制時(shí)減振效果可達(dá)70%。

(2)后緣小翼智能旋翼槳轂振動(dòng)控制的同時(shí)通常伴隨著槳葉根部扭矩增大的現(xiàn)象，帶來(lái)附加的低頻振動(dòng)，進(jìn)行槳葉設(shè)計(jì)時(shí)需考慮可能存在的疲勞問(wèn)題。

(3)翼型后緣小翼對(duì)扭轉(zhuǎn)剛硬槳葉的載荷抑制效果較差。若忽略小翼偏轉(zhuǎn)引起的俯仰力矩，則小翼無(wú)法有效實(shí)現(xiàn)減振。

[ 1 ] FRIEDMANN P P. On-blade control of rotor vibration, noise, and performance: Just around the corner? The 33rd Alexander Nikolsky Honorary Lecture[J]. Journal of the American Helicopter Society, 2014, 59(4): 1-37.

[ 2 ] FALLS J, DATTA A, CHOPRA I. Integrated trailing-edge flaps and servotabs for helicopter primary control[J]. Journal of the American Helicopter Society, 2010, 55(3): 1-15.

[ 3 ] PATT D, LIU L, FRIEDMANN P P. Rotorcraft vibration reduction and noise prediction using a unified aeroelastic response simulation[J]. Journal of the American Helicopter Society, 2005, 50(1):95-106.

[ 4 ] DATTA A. Fundamental understanding, prediction and validation of rotor vibratory loads in steady-level flight[D]. College Park: University of Maryland, 2004.

[ 5 ] DATTA A, CHOPRA I. Validation and understanding of UH-60A vibratory loads in steady level flight[J]. Journal of the American Helicopter Society, 2004, 49(3): 271-287.

[ 6 ] KIM J S, SMITH E C, WANG K W. Helicopter blade loads control via multiple trailing-edge flaps[C]//AHS 62nd Annual Forum. Phoenix AZ: AHS International, 2006: 1000-1010.

[ 7 ] 虞志浩, 楊衛(wèi)東, 張呈林. 基于Broyden法的旋翼多體系統(tǒng)氣動(dòng)彈性分析[J]. 航空學(xué)報(bào), 2012, 33(12): 2171-2182. YU Zhihao, YANG Weidong, ZHANG Chenglin. Aeroelasticity analysis of rotor multibody system based on Broyden method[J]. Acta Aeronautica et Astronautica Sinica, 2012, 33(12): 2171-2182.

[ 8 ] THEODORSEN T, GARRICK I E. Nonstationary flow about a wing-aileron-tab combination including aerodynamic balance[R]. NACA Report 736, 1942.

[ 9 ] 譚劍鋒, 王浩文, 吳超, 等. 基于非定常面元/黏性渦粒子混合法的旋翼/平尾非定常氣動(dòng)干擾[J]. 航空學(xué)報(bào)，2014，35(3): 643-656. TAN Jianfeng, WANG Haowen, WU Chao, et al. Rotor/empennage un-steady aerodynamic interaction with unsteady pan-el/viscous vortex particle hybrid method[J]. Acta Aero-nautica et Astronautica Sinica, 2014, 35(3): 643-656.

[10] 吳杰, 楊衛(wèi)東, 虞志浩. 旋翼槳葉結(jié)構(gòu)載荷計(jì)算方法比較研究[J]. 振動(dòng)與沖擊, 2014, 33(7): 210-214. WU Jie, YANG Weidong, YU Zhihao. Comparison among rotor blade structural load calculation methods[J]. Journal of Vibration and Shock, 2014, 33(7): 210-214.

[11] KORATKAR N A, CHOPRA I. Analysis and testing of mach-scaled rotor with trailing edge flaps[J]. AIAA Journal, 2000, 38(7): 1113-1124.

[12] 張柱，黃文俊，楊衛(wèi)東. 后緣小翼型智能旋翼槳葉模型設(shè)計(jì)分析與試驗(yàn)研究[J]. 南京航空航天大學(xué)學(xué)報(bào), 2011, 43(3): 296-301. ZHANG Zhu, HUANG Wenjun, YANG Weidong. Design analysis and test of smart rotor blades model with trailing edge flaps[J]. Journal of Nanjing University of Aeronautics & Astronautics, 2011, 43(3): 296-301.

[13] STALEY J A. Validation of rotorcraft flight simulation program through correlation with flight data for soft-in-plane hingeless rotors[R]. AMRDL-TR-75-50, 1976.

[14] HARRIS T A, LOWRY J G. Pressure distribution over an NACA 23012 airfoil with a fixed slot and a slotted flap[R]. NACA-TR-633, 1942.

[15] PETERSON R L, MAIER T. Correlation of wind tunnel and flight test results of a full-scale hingeless rotor[C]//AHS Aeromechanics Specialists Conference. Washington DC: AHS International, 1994.

[16] 張曉谷. 直升機(jī)動(dòng)力學(xué)設(shè)計(jì)[M]. 北京: 航空工業(yè)出版社, 1995.

Influence factors analysis for smart vibration control of a rotor wing with trailing edge flaps

LIU Shiming1, YANG Weidong1, YU Zhihao1, WU Jie2

(1. National Key Laboratory of Rotorcraft Aeromechanics, College of Aerospace Engineering,Nanjing University of Aeronautics and Astronautics, Nanjing 210000, China;2. School of Naval Architecture & Ocean Engineering,Jiangsu University of Science and Technology, Zhenjiang 212000, China)

An aeroelastic load computation model and the optimization method of vibration control for a smart rotor wing with trailing edge flags were developed. The influences of aerodynamic forces and inertial forces of rigid body trailing edge flaps on its elastic blades system were considered. The viscous vortex particle method with the airfoil look-up table was used to compute the aerodynamic load of the rotor wing. The vibratory blade and hub loads were predicted with the force integration method. An optimal algorithm combining the steepest decent method and the golden section algorithm was derived for the defined objective function including the blade torsion and hub vibration loads. It was demonstrated that the proposed model can effectively control the vibratory loads, and the objective function can be reduced by about 70%; the torsion of elastic blades ensures the trailing edge flaps can effectively reduce vibration, while it amplifies the aerodynamic pitch moment of the trailing edge flaps, the latter brings additional blade torsional moment.

helicopter; rotor wing; vibration; trailing edge flap; optimization method

國(guó)家自然科學(xué)基金(11272148)；江蘇高校優(yōu)勢(shì)學(xué)科建設(shè)工程；旋翼動(dòng)力學(xué)國(guó)防科技重點(diǎn)實(shí)驗(yàn)室基金(9140C400401140C40183)

2015-10-30 修改稿收到日期：2016-01-18

劉士明 男，博士生，1991年1月生

楊衛(wèi)東 男，博士，教授，1967年5月生

V211.47

A

10.13465/j.cnki.jvs.2017.03.022