• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

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

    2017-03-09 07:56:29劉士明楊衛(wèi)東虞志浩
    振動(dòng)與沖擊 2017年3期
    關(guān)鍵詞:小翼后緣氣動(dòng)力

    劉士明, 楊衛(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

    猜你喜歡
    小翼后緣氣動(dòng)力
    漢字獵人(一)
    我家養(yǎng)了一只紙精靈(二)
    我家養(yǎng)了一只紙精靈(四)
    飛行載荷外部氣動(dòng)力的二次規(guī)劃等效映射方法
    我是霸王龍
    機(jī)翼后緣連續(xù)變彎度對(duì)客機(jī)氣動(dòng)特性影響
    側(cè)風(fēng)對(duì)拍動(dòng)翅氣動(dòng)力的影響
    柔性后緣可變形機(jī)翼氣動(dòng)特性分析
    TNF-α和PGP9.5在椎體后緣離斷癥軟骨終板的表達(dá)及意義
    高速鐵路接觸線覆冰后氣動(dòng)力特性的風(fēng)洞試驗(yàn)研究
    欧美色视频一区免费| 国产真实乱freesex| 成人三级黄色视频| 三级男女做爰猛烈吃奶摸视频| 婷婷丁香在线五月| 久久久久九九精品影院| 91麻豆精品激情在线观看国产| 午夜福利在线观看免费完整高清在 | 亚洲乱码一区二区免费版| 久久久久免费精品人妻一区二区| 亚洲国产中文字幕在线视频| 国产综合懂色| 制服人妻中文乱码| 成人三级做爰电影| 1024手机看黄色片| 一本久久中文字幕| 免费在线观看亚洲国产| 香蕉国产在线看| bbb黄色大片| 中出人妻视频一区二区| 国内久久婷婷六月综合欲色啪| 久久热在线av| 黄色片一级片一级黄色片| 亚洲欧美一区二区三区黑人| 一进一出抽搐gif免费好疼| 久久国产精品人妻蜜桃| 欧美色视频一区免费| 综合色av麻豆| www.自偷自拍.com| 黑人操中国人逼视频| 欧美日韩一级在线毛片| 国产精华一区二区三区| 亚洲熟妇中文字幕五十中出| 最好的美女福利视频网| 亚洲国产高清在线一区二区三| 手机成人av网站| 亚洲欧美日韩高清在线视频| 欧美成人免费av一区二区三区| 成人国产综合亚洲| 最新美女视频免费是黄的| 精品久久蜜臀av无| 成在线人永久免费视频| www.自偷自拍.com| 国产亚洲精品综合一区在线观看| 97超级碰碰碰精品色视频在线观看| 搡老妇女老女人老熟妇| 一进一出好大好爽视频| 日韩精品青青久久久久久| 成人高潮视频无遮挡免费网站| 成人三级做爰电影| 午夜福利视频1000在线观看| or卡值多少钱| 久久人人精品亚洲av| 成人永久免费在线观看视频| 欧美极品一区二区三区四区| 人妻久久中文字幕网| 中出人妻视频一区二区| 校园春色视频在线观看| 欧美黄色淫秽网站| 中文字幕精品亚洲无线码一区| 久久这里只有精品19| 欧美性猛交黑人性爽| 男人舔女人的私密视频| 国产精品野战在线观看| av片东京热男人的天堂| 一夜夜www| 亚洲国产看品久久| 亚洲专区国产一区二区| 亚洲在线观看片| 在线观看免费午夜福利视频| 亚洲熟妇熟女久久| 麻豆久久精品国产亚洲av| 国产成人av激情在线播放| 成年人黄色毛片网站| 久久香蕉国产精品| 少妇人妻一区二区三区视频| 亚洲色图av天堂| 午夜久久久久精精品| 亚洲国产色片| 色吧在线观看| 又粗又爽又猛毛片免费看| 欧美乱码精品一区二区三区| 99久久综合精品五月天人人| 99久久无色码亚洲精品果冻| 欧美午夜高清在线| 亚洲成人免费电影在线观看| 日日夜夜操网爽| 欧美黄色淫秽网站| 一本久久中文字幕| 国产成年人精品一区二区| 国产欧美日韩精品亚洲av| 日本免费a在线| 一个人免费在线观看电影 | 国产单亲对白刺激| 青草久久国产| 欧美性猛交╳xxx乱大交人| 色吧在线观看| 国产男靠女视频免费网站| 亚洲国产高清在线一区二区三| 亚洲va日本ⅴa欧美va伊人久久| 日韩大尺度精品在线看网址| 欧美另类亚洲清纯唯美| 国产成人一区二区三区免费视频网站| 午夜福利在线观看吧| 少妇裸体淫交视频免费看高清| 长腿黑丝高跟| 99久久精品热视频| 亚洲国产精品sss在线观看| 一进一出抽搐gif免费好疼| 两个人的视频大全免费| av福利片在线观看| 亚洲色图av天堂| 成人高潮视频无遮挡免费网站| 特级一级黄色大片| 久久精品亚洲精品国产色婷小说| 波多野结衣巨乳人妻| 啦啦啦韩国在线观看视频| 男女那种视频在线观看| 日本三级黄在线观看| 特大巨黑吊av在线直播| 国产成人一区二区三区免费视频网站| 国产亚洲av嫩草精品影院| 嫩草影院入口| 嫩草影院精品99| 成年女人永久免费观看视频| 久久精品影院6| 男女那种视频在线观看| 女人被狂操c到高潮| 久99久视频精品免费| a在线观看视频网站| 一夜夜www| 欧美日韩乱码在线| 国产精品爽爽va在线观看网站| 亚洲精品色激情综合| 欧美黑人欧美精品刺激| 91麻豆精品激情在线观看国产| 黄色女人牲交| 日韩欧美一区二区三区在线观看| 巨乳人妻的诱惑在线观看| 午夜影院日韩av| 人人妻人人澡欧美一区二区| 亚洲精品在线观看二区| 日韩欧美 国产精品| 亚洲熟女毛片儿| 精品不卡国产一区二区三区| 亚洲午夜精品一区,二区,三区| 国产免费男女视频| 91在线观看av| 午夜日韩欧美国产| 波多野结衣巨乳人妻| 99国产精品99久久久久| 美女被艹到高潮喷水动态| 免费观看精品视频网站| 亚洲av成人精品一区久久| 最近最新中文字幕大全电影3| 亚洲男人的天堂狠狠| 欧美极品一区二区三区四区| 99久久综合精品五月天人人| 国产单亲对白刺激| 国产精品亚洲av一区麻豆| 99久久精品国产亚洲精品| x7x7x7水蜜桃| 成年免费大片在线观看| 非洲黑人性xxxx精品又粗又长| 熟女人妻精品中文字幕| 亚洲人成电影免费在线| 别揉我奶头~嗯~啊~动态视频| 亚洲色图 男人天堂 中文字幕| 禁无遮挡网站| 1024手机看黄色片| 一级作爱视频免费观看| 一进一出抽搐动态| 午夜激情欧美在线| 国产亚洲欧美98| 亚洲欧美精品综合久久99| 午夜福利成人在线免费观看| 精品久久久久久成人av| 欧美中文综合在线视频| 日本三级黄在线观看| 高清在线国产一区| 午夜成年电影在线免费观看| 亚洲精品久久国产高清桃花| 国产精品 国内视频| 色精品久久人妻99蜜桃| 麻豆av在线久日| 久久香蕉国产精品| 中文字幕熟女人妻在线| 亚洲九九香蕉| 午夜亚洲福利在线播放| 亚洲av熟女| 中国美女看黄片| 国产精华一区二区三区| 别揉我奶头~嗯~啊~动态视频| 精品一区二区三区四区五区乱码| 国产亚洲欧美98| 熟女少妇亚洲综合色aaa.| 99久久成人亚洲精品观看| 国产综合懂色| 成人欧美大片| 欧美乱码精品一区二区三区| 99国产精品一区二区三区| 成年版毛片免费区| 亚洲电影在线观看av| 国产精品综合久久久久久久免费| 成人国产一区最新在线观看| 午夜福利18| 最新中文字幕久久久久 | 成在线人永久免费视频| 51午夜福利影视在线观看| 国产精品久久视频播放| 久久午夜亚洲精品久久| 国产成人一区二区三区免费视频网站| 国产单亲对白刺激| 国产av一区在线观看免费| 亚洲 欧美一区二区三区| 久久精品国产综合久久久| 国产精品国产高清国产av| 变态另类成人亚洲欧美熟女| 黑人巨大精品欧美一区二区mp4| 我的老师免费观看完整版| 国产av在哪里看| 国产不卡一卡二| 亚洲美女视频黄频| 国产精品,欧美在线| 日韩精品青青久久久久久| 成人国产综合亚洲| 日本在线视频免费播放| 香蕉国产在线看| 最近视频中文字幕2019在线8| 黑人操中国人逼视频| 99久久国产精品久久久| а√天堂www在线а√下载| 久久99热这里只有精品18| 黄色 视频免费看| 老司机深夜福利视频在线观看| 国产av一区在线观看免费| 久久这里只有精品19| 免费看十八禁软件| 噜噜噜噜噜久久久久久91| 日韩欧美一区二区三区在线观看| 波多野结衣高清作品| 国产黄色小视频在线观看| 国产亚洲精品久久久com| 久99久视频精品免费| 色av中文字幕| 一级毛片精品| 欧美日韩精品网址| 老熟妇仑乱视频hdxx| 国产精品,欧美在线| 操出白浆在线播放| 亚洲无线观看免费| av天堂在线播放| 久久中文字幕人妻熟女| 色视频www国产| 搡老熟女国产l中国老女人| 熟女电影av网| 男人舔女人下体高潮全视频| 亚洲熟妇熟女久久| 国内毛片毛片毛片毛片毛片| 久久久久久久精品吃奶| 久久九九热精品免费| 久99久视频精品免费| 久久草成人影院| 久久人妻av系列| 看黄色毛片网站| 成年版毛片免费区| 麻豆久久精品国产亚洲av| 国产三级在线视频| 亚洲自拍偷在线| 亚洲自拍偷在线| 婷婷精品国产亚洲av在线| 欧美+亚洲+日韩+国产| 午夜免费成人在线视频| 91字幕亚洲| 亚洲av成人不卡在线观看播放网| 国产av在哪里看| 午夜福利欧美成人| 在线免费观看的www视频| 男女床上黄色一级片免费看| 国产欧美日韩精品亚洲av| 国产精品九九99| 嫩草影院精品99| 91久久精品国产一区二区成人 | 亚洲人成网站在线播放欧美日韩| 日韩免费av在线播放| 男女做爰动态图高潮gif福利片| 日韩欧美三级三区| 亚洲中文av在线| 国产高清有码在线观看视频| 给我免费播放毛片高清在线观看| 久99久视频精品免费| 国产精品亚洲av一区麻豆| 成人18禁在线播放| 中文字幕久久专区| 一区二区三区国产精品乱码| 婷婷精品国产亚洲av| av在线蜜桃| 午夜影院日韩av| 热99在线观看视频| 成年女人永久免费观看视频| 偷拍熟女少妇极品色| 久久午夜亚洲精品久久| 狂野欧美白嫩少妇大欣赏| 又爽又黄无遮挡网站| 午夜久久久久精精品| 国产欧美日韩精品一区二区| 亚洲国产欧美网| 日韩精品中文字幕看吧| 免费无遮挡裸体视频| 亚洲国产高清在线一区二区三| 成人永久免费在线观看视频| 国产在线精品亚洲第一网站| 久久草成人影院| 免费观看的影片在线观看| 亚洲九九香蕉| 日本精品一区二区三区蜜桃| 久久国产乱子伦精品免费另类| 亚洲专区国产一区二区| 日日夜夜操网爽| 校园春色视频在线观看| 日韩有码中文字幕| 精品久久久久久久末码| av在线天堂中文字幕| 变态另类丝袜制服| 国产成人精品无人区| 亚洲成人久久爱视频| 欧洲精品卡2卡3卡4卡5卡区| 一级毛片高清免费大全| 又黄又爽又免费观看的视频| 国产精品综合久久久久久久免费| 少妇丰满av| 露出奶头的视频| 亚洲精华国产精华精| 久久久久性生活片| 黑人欧美特级aaaaaa片| 又黄又粗又硬又大视频| 麻豆久久精品国产亚洲av| 美女大奶头视频| 久久久久久久久久黄片| 国产黄色小视频在线观看| 最近在线观看免费完整版| 精品日产1卡2卡| 少妇的逼水好多| 亚洲中文字幕一区二区三区有码在线看 | 亚洲一区二区三区色噜噜| 久久天堂一区二区三区四区| 99久久无色码亚洲精品果冻| 两人在一起打扑克的视频| 一级毛片精品| 曰老女人黄片| 97碰自拍视频| 日韩欧美三级三区| 日韩av在线大香蕉| 女警被强在线播放| 一级毛片女人18水好多| 欧美精品啪啪一区二区三区| 一级a爱片免费观看的视频| 亚洲欧美一区二区三区黑人| 久久精品国产99精品国产亚洲性色| 91九色精品人成在线观看| 午夜福利在线在线| 黄色视频,在线免费观看| xxxwww97欧美| 欧美日本视频| 精品人妻1区二区| a级毛片a级免费在线| 午夜激情福利司机影院| 色综合婷婷激情| 亚洲精品一区av在线观看| 亚洲九九香蕉| 狂野欧美激情性xxxx| 国产精品乱码一区二三区的特点| 18禁黄网站禁片免费观看直播| 国产高清有码在线观看视频| 国产高清videossex| 日韩中文字幕欧美一区二区| 黄片小视频在线播放| 国产亚洲精品久久久com| 日韩大尺度精品在线看网址| 人妻夜夜爽99麻豆av| 久久久久久久午夜电影| 非洲黑人性xxxx精品又粗又长| 久久精品国产综合久久久| 国产野战对白在线观看| 51午夜福利影视在线观看| 亚洲av电影不卡..在线观看| 亚洲精品一卡2卡三卡4卡5卡| 香蕉av资源在线| 亚洲色图 男人天堂 中文字幕| 99热精品在线国产| 嫩草影视91久久| 免费看美女性在线毛片视频| av片东京热男人的天堂| 国产精品乱码一区二三区的特点| 少妇丰满av| 麻豆av在线久日| 国产成人av激情在线播放| 国产高潮美女av| 听说在线观看完整版免费高清| 少妇的丰满在线观看| 国语自产精品视频在线第100页| 欧美午夜高清在线| tocl精华| 成人精品一区二区免费| 日本与韩国留学比较| 国产成人精品无人区| 黄色丝袜av网址大全| 日本黄大片高清| 成人特级av手机在线观看| 女人被狂操c到高潮| 成人国产综合亚洲| 亚洲激情在线av| 欧美成狂野欧美在线观看| 看片在线看免费视频| www.自偷自拍.com| 搡老妇女老女人老熟妇| 精品久久久久久成人av| 亚洲熟女毛片儿| 麻豆国产97在线/欧美| 最新在线观看一区二区三区| 激情在线观看视频在线高清| 1024香蕉在线观看| 国产成人aa在线观看| 亚洲专区字幕在线| 男女做爰动态图高潮gif福利片| 国产精华一区二区三区| 精品久久久久久久人妻蜜臀av| 精品久久久久久,| 18美女黄网站色大片免费观看| 国产探花在线观看一区二区| 亚洲真实伦在线观看| 免费搜索国产男女视频| 91老司机精品| 午夜免费成人在线视频| 国产精品久久久久久人妻精品电影| 男女那种视频在线观看| 亚洲真实伦在线观看| 日韩av在线大香蕉| 亚洲在线观看片| 免费看a级黄色片| 麻豆av在线久日| 丁香六月欧美| 美女免费视频网站| 久久久久九九精品影院| 成人精品一区二区免费| 无遮挡黄片免费观看| www.www免费av| 久久国产精品人妻蜜桃| 一个人看的www免费观看视频| 老汉色av国产亚洲站长工具| 国产单亲对白刺激| 全区人妻精品视频| 欧美日韩中文字幕国产精品一区二区三区| 欧美黄色淫秽网站| 国产精品久久久久久久电影 | АⅤ资源中文在线天堂| 欧美黑人巨大hd| 中文字幕av在线有码专区| 美女大奶头视频| 一本一本综合久久| 欧美日韩国产亚洲二区| 国产91精品成人一区二区三区| 非洲黑人性xxxx精品又粗又长| 手机成人av网站| 成人三级黄色视频| 成人无遮挡网站| 亚洲午夜理论影院| 一进一出好大好爽视频| 国内精品久久久久久久电影| 又紧又爽又黄一区二区| 一个人免费在线观看的高清视频| 蜜桃久久精品国产亚洲av| 久久久久久九九精品二区国产| 欧美激情在线99| 国产激情偷乱视频一区二区| 热99在线观看视频| 99久久国产精品久久久| 又紧又爽又黄一区二区| 成在线人永久免费视频| 久久午夜综合久久蜜桃| 亚洲美女黄片视频| 亚洲成av人片免费观看| 不卡av一区二区三区| 九色成人免费人妻av| 精品99又大又爽又粗少妇毛片 | 日本黄大片高清| av视频在线观看入口| 免费一级毛片在线播放高清视频| 久久这里只有精品中国| 在线看三级毛片| 国产极品精品免费视频能看的| 51午夜福利影视在线观看| 精品国产乱子伦一区二区三区| 色噜噜av男人的天堂激情| 亚洲人与动物交配视频| 国产精品98久久久久久宅男小说| 噜噜噜噜噜久久久久久91| 观看免费一级毛片| av国产免费在线观看| bbb黄色大片| 亚洲熟女毛片儿| 欧美zozozo另类| netflix在线观看网站| 岛国视频午夜一区免费看| 亚洲av电影不卡..在线观看| 日韩欧美国产一区二区入口| 桃红色精品国产亚洲av| 亚洲国产精品合色在线| 淫秽高清视频在线观看| 最新美女视频免费是黄的| 国产精品久久久人人做人人爽| 一边摸一边抽搐一进一小说| 欧美绝顶高潮抽搐喷水| 亚洲国产欧美人成| 久久人妻av系列| 男女下面进入的视频免费午夜| 国产日本99.免费观看| 老汉色∧v一级毛片| 麻豆av在线久日| 搞女人的毛片| 久久久国产精品麻豆| 麻豆一二三区av精品| 黄片小视频在线播放| 国产成人啪精品午夜网站| 色在线成人网| 国产成人av教育| 色哟哟哟哟哟哟| 十八禁人妻一区二区| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲精品国产精品久久久不卡| 国产 一区 欧美 日韩| 国内揄拍国产精品人妻在线| 欧美丝袜亚洲另类 | 婷婷精品国产亚洲av| 日本免费a在线| 国产伦精品一区二区三区四那| 18禁国产床啪视频网站| 校园春色视频在线观看| 午夜福利18| 亚洲成人久久爱视频| 亚洲人成网站在线播放欧美日韩| 最近最新中文字幕大全免费视频| 9191精品国产免费久久| 国产一区二区在线av高清观看| 国产精品一区二区免费欧美| 久久久成人免费电影| 成年女人毛片免费观看观看9| av天堂在线播放| 国产成人一区二区三区免费视频网站| av在线蜜桃| 成年人黄色毛片网站| 久久香蕉精品热| 午夜福利在线观看吧| 欧美日韩瑟瑟在线播放| 欧美日韩乱码在线| 好看av亚洲va欧美ⅴa在| 90打野战视频偷拍视频| 啦啦啦观看免费观看视频高清| 国产91精品成人一区二区三区| 在线免费观看不下载黄p国产 | 婷婷亚洲欧美| а√天堂www在线а√下载| 亚洲专区字幕在线| av在线天堂中文字幕| 在线观看舔阴道视频| 亚洲狠狠婷婷综合久久图片| 中文字幕高清在线视频| 国产亚洲欧美在线一区二区| 丰满的人妻完整版| 国产成人啪精品午夜网站| 久久这里只有精品中国| 熟女少妇亚洲综合色aaa.| 亚洲乱码一区二区免费版| 亚洲人成伊人成综合网2020| 97碰自拍视频| 精品无人区乱码1区二区| 97碰自拍视频| 久久久久性生活片| av中文乱码字幕在线| 99久久精品国产亚洲精品| 人人妻,人人澡人人爽秒播| 国产成人一区二区三区免费视频网站| 大型黄色视频在线免费观看| 亚洲成人久久爱视频| 18禁观看日本| 亚洲人成伊人成综合网2020| 此物有八面人人有两片| 国内精品久久久久久久电影| 女人高潮潮喷娇喘18禁视频| av片东京热男人的天堂| 久久国产乱子伦精品免费另类| 少妇的丰满在线观看| 国产真人三级小视频在线观看| 成人三级做爰电影| 热99在线观看视频| 99久久成人亚洲精品观看| 国产91精品成人一区二区三区| 日本一本二区三区精品| 国产精品免费一区二区三区在线| 亚洲av成人精品一区久久| 久久精品影院6| 精品久久久久久久人妻蜜臀av| 露出奶头的视频| 又黄又粗又硬又大视频| 亚洲av熟女| 黑人巨大精品欧美一区二区mp4| 禁无遮挡网站| 国产爱豆传媒在线观看| 99久久精品国产亚洲精品| 精华霜和精华液先用哪个| 国产久久久一区二区三区| 97超视频在线观看视频| 一级毛片精品| 国内毛片毛片毛片毛片毛片| 激情在线观看视频在线高清|