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

    微型共軸雙旋翼氣動(dòng)性能數(shù)值模擬與試驗(yàn)分析

    2017-11-07 10:06:35紀(jì)玉霞汪長(zhǎng)煒
    實(shí)驗(yàn)流體力學(xué) 2017年5期
    關(guān)鍵詞:共軸旋翼拉力

    雷 瑤,紀(jì)玉霞,汪長(zhǎng)煒

    (福州大學(xué) 機(jī)械工程及自動(dòng)化學(xué)院,福州 350116)

    微型共軸雙旋翼氣動(dòng)性能數(shù)值模擬與試驗(yàn)分析

    雷 瑤*,紀(jì)玉霞,汪長(zhǎng)煒

    (福州大學(xué) 機(jī)械工程及自動(dòng)化學(xué)院,福州 350116)

    為研究懸停狀態(tài)下旋翼的間距對(duì)微型共軸雙旋翼氣動(dòng)性能的影響,文中通過(guò)搭建試驗(yàn)平臺(tái)對(duì)間距比h/r分別為0.32、0.38、0.45、0.51、0.58、0.65和0.75下的共軸雙旋翼進(jìn)行氣動(dòng)性能測(cè)試,以測(cè)量不同旋翼轉(zhuǎn)速下所得共軸雙旋翼的拉力和功耗對(duì)共軸雙旋翼氣動(dòng)布局進(jìn)行優(yōu)化,試圖找出具有最佳氣動(dòng)特性的共軸旋翼布局。另外,通過(guò)試驗(yàn)誤差分析確定了相應(yīng)的拉力系數(shù)、功率系數(shù)和功率載荷,且試驗(yàn)誤差均小于2%。同時(shí),為更直觀得到不同間距下氣流干擾對(duì)旋翼系統(tǒng)氣動(dòng)性能的影響,文中采用數(shù)值模擬得到了不同間距比下旋翼的流線分布和壓力分布。最后,對(duì)比試驗(yàn)結(jié)果,綜合分析旋翼間氣動(dòng)干擾的影響,最終得到間距比h/r為0.38時(shí)的共軸雙旋翼具有最佳的氣動(dòng)布局。研究結(jié)果表明,懸停狀態(tài)的共軸雙旋翼可以通過(guò)改變間距大大提高氣動(dòng)性能,且同一間距下轉(zhuǎn)速越大雖然旋翼間干擾越強(qiáng)烈,但此時(shí)開(kāi)始出現(xiàn)耦合,使得系統(tǒng)的氣動(dòng)性能可能更好,同時(shí),由于上下旋翼間的相互誘導(dǎo),雖然轉(zhuǎn)速較低時(shí)升力較小,但是功耗明顯低于高轉(zhuǎn)速,使得系統(tǒng)具有更大的功率載荷。

    共軸雙旋翼;懸停狀態(tài);氣動(dòng)干擾;數(shù)值模擬;氣動(dòng)布局

    0 引 言

    旋翼作為整機(jī)動(dòng)力部件,旋翼數(shù)量、旋翼拉力和功耗決定了系統(tǒng)載荷和續(xù)航時(shí)間[1]。共軸雙旋翼,由于上下旋翼反轉(zhuǎn),所有功率均用于提供垂直方向拉力,具有較好的懸停性能。此外,共軸雙旋翼的尺寸小、重量輕、結(jié)構(gòu)緊湊、氣動(dòng)效率高,并且可以縮小整體尺寸至所需重量[2-4],對(duì)軍用和民用領(lǐng)域具有重要意義。

    然而旋翼間存在的氣動(dòng)干擾和復(fù)雜的流場(chǎng)環(huán)境會(huì)影響系統(tǒng)的性能和效率[5-6],對(duì)于微型共軸雙旋翼,其上下旋翼重疊區(qū)域內(nèi)的氣動(dòng)干擾較為嚴(yán)重[7-8],而其重疊區(qū)域的大小直接由旋翼的間距決定。此外,懸停是共軸雙旋翼消耗較高功率[9]的工作狀態(tài),在此狀態(tài)下,上下旋翼槳尖渦碰撞使得懸停狀態(tài)下的旋翼系統(tǒng)內(nèi)部氣動(dòng)干擾嚴(yán)重且耦合復(fù)雜。

    國(guó)內(nèi)外針對(duì)共軸雙旋翼間的氣動(dòng)干擾均開(kāi)展了相關(guān)研究。Bohorquez[10]對(duì)半徑為172mm微型共軸飛行器進(jìn)行試驗(yàn),發(fā)現(xiàn)只有當(dāng)軸向間距超過(guò)35%旋翼半徑時(shí),上旋翼才受下旋翼少量的影響。Syal[11]使用動(dòng)量理論研究半徑為3810mm的Harrington旋翼,得到性能最優(yōu)的共軸雙旋翼系統(tǒng),并發(fā)現(xiàn)旋翼間的相互干擾使上下旋翼的性能下降。雷瑤[12-13]等設(shè)計(jì)試驗(yàn)平臺(tái)分析不同間距下的共軸雙旋翼氣動(dòng)性能,并通過(guò)風(fēng)洞進(jìn)行水平和垂直的來(lái)流試驗(yàn),發(fā)現(xiàn)旋翼間干擾小、性能較佳的軸向間距及水平來(lái)流對(duì)旋翼性能影響較大。在試驗(yàn)研究方面,Lakshminarayan[14]等發(fā)展高精度的計(jì)算平臺(tái)來(lái)研究具有不同氣動(dòng)布局的微型共軸雙旋翼中流場(chǎng)的氣動(dòng)干擾與間距的關(guān)系。Coleman[15]通過(guò)大量的試驗(yàn)數(shù)據(jù)分析了旋翼間距和載荷分配的問(wèn)題。馬楊超[16]等采用水洞實(shí)驗(yàn)測(cè)量技術(shù)對(duì)懸停狀態(tài)下的共軸雙旋翼流場(chǎng)特性進(jìn)行了研究并得到流場(chǎng)矢量圖和渦量云圖數(shù)據(jù)。姬樂(lè)強(qiáng)[17]等設(shè)計(jì)了上、下旋翼不同安裝形式的試驗(yàn)裝置并通過(guò)風(fēng)洞試驗(yàn)研究了共軸雙旋翼的氣動(dòng)特性。

    目前,專門(mén)針對(duì)旋翼間距對(duì)氣動(dòng)性能影響的研究相對(duì)較少且試驗(yàn)設(shè)計(jì)復(fù)雜,特別是對(duì)于懸停狀態(tài)下的微型共軸旋翼系統(tǒng)最佳間距的確定。另外,大多數(shù)研究還僅僅只是對(duì)模型旋翼進(jìn)行試驗(yàn),在旋翼系統(tǒng)性能測(cè)試的過(guò)程中,外界擾動(dòng)和測(cè)量誤差常被忽略,且綜合試驗(yàn)與模擬結(jié)果的分析比較少。

    綜上,本文通過(guò)搭建試驗(yàn)平臺(tái)對(duì)懸停狀態(tài)下多個(gè)軸向間距的共軸雙旋翼的氣動(dòng)特性進(jìn)行了測(cè)試,并確定了試驗(yàn)誤差。另外,結(jié)合數(shù)值模擬通過(guò)旋翼間流線的變化分析了不同間距上下旋翼間的氣動(dòng)干擾,進(jìn)而得出了最佳共軸間距,優(yōu)化了共軸雙旋翼的布局。

    1 理論分析

    本文所采用的旋翼半徑r為200mm,弦長(zhǎng)為35mm,槳葉2個(gè),轉(zhuǎn)速范圍1600~2300r/min,槳尖馬赫數(shù)變化范圍0.1~0.14,槳尖雷諾數(shù)范圍0.8×105~1.16×105。定義旋翼間距為h,為忽略旋翼尺寸帶來(lái)的影響,定義無(wú)量綱旋翼間距比為h/r。

    共軸雙旋翼最佳氣動(dòng)布局要求能產(chǎn)生較大拉力的同時(shí)降低功耗,在實(shí)測(cè)數(shù)據(jù)基礎(chǔ)上結(jié)合數(shù)值模擬理論分析旋翼間流場(chǎng)干擾耦合,其優(yōu)化設(shè)計(jì)過(guò)程如下:

    (1) 確定共軸間距范圍。為避免旋翼間碰撞并考慮到整機(jī)空間尺寸,間距范圍在0.2r~0.8r之間。為細(xì)化間距變化程度,根據(jù)初步試驗(yàn)結(jié)果,分別取間距比h/r為0.32、0.38、0.45、0.51、0.58、0.65和0.75;

    (2) 搭建旋翼試驗(yàn)平臺(tái)測(cè)試工作轉(zhuǎn)速下各間距的系統(tǒng)拉力和功耗;

    (3) 計(jì)算功率載荷,根據(jù)氣動(dòng)性能確定共軸間距的優(yōu)劣;

    (4) 結(jié)合數(shù)值模擬,根據(jù)得到的各間距流線分布規(guī)律分析氣動(dòng)干擾對(duì)氣動(dòng)性能的影響,對(duì)試驗(yàn)結(jié)果進(jìn)行分析解釋作為理論分析支撐;

    (5) 對(duì)比試驗(yàn)和數(shù)值計(jì)算,確定最終最佳間距作為共軸雙旋翼氣動(dòng)布局。

    2 試驗(yàn)研究

    通過(guò)試驗(yàn)平臺(tái)的搭建來(lái)對(duì)共軸雙旋翼氣動(dòng)特性進(jìn)行測(cè)量可以反映旋翼真實(shí)的工作狀態(tài),提供可信的數(shù)據(jù)。

    2.1搭建試驗(yàn)平臺(tái)

    旋翼氣動(dòng)特性試驗(yàn)平臺(tái)主要包括:動(dòng)力系統(tǒng),由直流電源(型號(hào):格氏ace鋰聚合物電池)、直流無(wú)刷電機(jī)(型號(hào):EM2835)和PWM調(diào)速系統(tǒng)組成;傳動(dòng)系統(tǒng),包括齒輪傳動(dòng)系統(tǒng)和內(nèi)外旋翼軸,旋翼軸由2個(gè)電機(jī)分別帶動(dòng)使上、下旋翼反向旋轉(zhuǎn),每個(gè)旋翼?yè)碛歇?dú)立傳動(dòng)機(jī)構(gòu),減小了引入誤差。另外,在旋翼軸上加工不同位置的定位孔可以改變旋翼間距;測(cè)量系統(tǒng),主要測(cè)量旋翼的轉(zhuǎn)速、拉力以及功率,轉(zhuǎn)速由轉(zhuǎn)速計(jì)讀出(型號(hào):TM-5010K,精度:±0.01%±1d),拉力通過(guò)微型壓力傳感器獲得(型號(hào):PLD204D-19,精度:0.5%%F.S),功率則根據(jù)記錄的電流和電壓值后處理得到;支座,起支撐旋翼傳動(dòng)系統(tǒng)的作用,試驗(yàn)中將旋翼反裝來(lái)避免產(chǎn)生地面效應(yīng)。為保證共軸雙旋翼系統(tǒng)實(shí)現(xiàn)扭矩平衡,固定上旋翼轉(zhuǎn)速調(diào)整下旋翼轉(zhuǎn)速使系統(tǒng)扭矩為0,上下旋翼轉(zhuǎn)速誤差在2%以內(nèi)。

    試驗(yàn)裝置示意如圖1所示,傳感器、轉(zhuǎn)速計(jì)以及電表測(cè)量值經(jīng)數(shù)據(jù)采集卡傳輸?shù)诫娔X中,以便計(jì)算和整理出不同間距下共軸雙旋翼的拉力、轉(zhuǎn)速和功率,這3個(gè)物理量可用來(lái)衡量旋翼系統(tǒng)的氣動(dòng)特性。

    圖1 試驗(yàn)裝置示意圖Fig.1 Sketch of experimental setup

    2.2試驗(yàn)結(jié)果分析

    拉力是決定旋翼系統(tǒng)載荷最直接的因素,不同間距下共軸雙旋翼拉力隨轉(zhuǎn)速變化如圖2所示,隨轉(zhuǎn)速增加,各間距下共軸系統(tǒng)的拉力和功率均增大。間距比h/r為0.32和0.38時(shí),旋翼拉力值具有明顯優(yōu)勢(shì),且當(dāng)轉(zhuǎn)速大于2050r/min時(shí),拉力增幅擴(kuò)大。這是由于在低轉(zhuǎn)速下,槳葉對(duì)低雷諾數(shù)下產(chǎn)生的粘性效應(yīng)更敏感,旋翼受到的阻力較大,影響整體的拉力水平。當(dāng)間距比大于0.45時(shí),共軸拉力減小趨于穩(wěn)定,并逐漸與旋翼的距離無(wú)關(guān)。由此得出小間距下的旋翼整體拉力較大,而隨著間距增大,上下洗效應(yīng)減弱,旋翼間的有利干擾減小,整體性能下降。

    圖2 不同間距下共軸雙旋翼拉力Fig.2 Thrust of a coaxial rotor with different spacing

    在保持較大拉力的同時(shí),功率損耗是衡量系統(tǒng)氣動(dòng)性能的另外一個(gè)重要指標(biāo),各間距功耗隨轉(zhuǎn)速變化如圖3所示,對(duì)比轉(zhuǎn)速1600~2300r/min內(nèi)旋翼消耗的功率,間距比0.38均小于間距比0.32,其他間距比下的功率消耗大體呈現(xiàn)出一致的趨勢(shì)。結(jié)合圖2~3,僅當(dāng)間距比為0.38時(shí),該旋翼系統(tǒng)具有最大拉力的同時(shí)擁有較低功耗。

    選取典型轉(zhuǎn)速1600r/min(低速)、1980r/min(中速)和2280r/min(高速)3種情況,共軸雙旋翼拉力隨間距變化如圖4所示,當(dāng)間距比為0.3~0.45時(shí),旋翼獲得較好的拉力水平,明顯觀察到當(dāng)間距比為0.38時(shí),旋翼產(chǎn)生的拉力最大,該間距比下的拉力優(yōu)勢(shì)隨轉(zhuǎn)速擴(kuò)大。當(dāng)間距比為0.45~0.75時(shí),隨間距增大,拉力值波動(dòng)較小并趨于穩(wěn)定。綜合圖2~4分析可得,間距比為0.38的共軸雙旋翼表現(xiàn)出較為優(yōu)異的氣動(dòng)性能,并且在高轉(zhuǎn)速下的旋翼性能得到進(jìn)一步提高。

    圖3 不同間距下共軸雙旋翼功率Fig.3 Power consumption of a coaxial rotor with different spacing

    圖4 典型轉(zhuǎn)速下共軸拉力隨間距比變化Fig.4 Coaxial thrust vs.rotor spacing at typical rotational speed

    為研究旋翼間的干擾對(duì)系統(tǒng)拉力的影響,3種情況下的旋翼拉力隨轉(zhuǎn)速變化如圖5所示,對(duì)比單旋翼的拉力、共軸雙旋翼的拉力以及單獨(dú)上下旋翼的拉力和,發(fā)現(xiàn)伴隨轉(zhuǎn)速增大,共軸雙旋翼的拉力增長(zhǎng)速度大于單旋翼。這是由于懸停狀態(tài)下,上下旋翼相鄰區(qū)域的相互誘導(dǎo)使共軸雙旋翼的拉力快速上升。此時(shí)對(duì)比共軸雙旋翼與單獨(dú)上下旋翼的拉力之和曲線,顯然共軸雙旋翼的拉力略小于其單獨(dú)上下旋翼產(chǎn)生的拉力總和。

    圖5 旋翼拉力隨轉(zhuǎn)速變化Fig.5 Thrust vs.rotational speed

    這是由于下旋翼處于上旋翼的尾跡中,上旋翼受下旋翼上洗的影響,其中下旋翼的工作環(huán)境更差,損失更嚴(yán)重,共軸旋翼間的相互干擾使上下旋翼的拉力相對(duì)同工作狀態(tài)下的單獨(dú)旋翼均有一定的損失。在槳盤(pán)載荷相同的情況下,功率載荷PL(Power Loading)是衡量旋翼懸停效率最重要的指標(biāo)之一,定義如下:

    式中:CT為拉力系數(shù);CP為功率系數(shù);Ω為轉(zhuǎn)速。其中,拉力系數(shù)和功率系數(shù)分別定義如下:

    式中:R為旋翼半徑;T為拉力;Q為扭矩;ρ為空氣密度;A為旋翼面積;P為功率。

    因此,功率載荷簡(jiǎn)化為:

    因此,得到間距比為0.38的共軸雙旋翼功率載荷隨轉(zhuǎn)速變化如圖6所示。低轉(zhuǎn)速下的功率載荷PL明顯比較大,這是由于低轉(zhuǎn)速?zèng)Q定低功率,雖然低轉(zhuǎn)速同時(shí)會(huì)導(dǎo)致拉力水平降低,但功率的下降速度小于拉力,故低轉(zhuǎn)速下的功率載荷具有優(yōu)勢(shì),這與Bohorquez[10]的觀點(diǎn)一致。

    圖6 功率載荷隨轉(zhuǎn)速變化Fig.6 Power Loading (PL) vs.Rotational speed

    2.3試驗(yàn)誤差分析

    試驗(yàn)平臺(tái)的測(cè)量誤差主要有2個(gè):轉(zhuǎn)速測(cè)量誤差和壓力傳感器誤差。其中轉(zhuǎn)速的誤差值與轉(zhuǎn)子磁鐵數(shù)量有關(guān),試驗(yàn)中采用24個(gè)磁鐵,因此任意轉(zhuǎn)速下精度為1/24。每一秒內(nèi),磁鐵經(jīng)過(guò)霍爾傳感器的次數(shù),存在誤差1/24×60=2.5r/min。另外,壓力傳感器精度為0.02%F.S.。

    另外,拉力計(jì)算誤差與轉(zhuǎn)速成比例增長(zhǎng)。應(yīng)用Kline-McClintock方程式得到的拉力系數(shù)誤差的計(jì)算如下:

    則:

    同理可得功率系數(shù)的計(jì)算誤差:

    代入試驗(yàn)測(cè)量的旋翼拉力、功率和轉(zhuǎn)速值,由Kline-McClintock方程計(jì)算得拉力系數(shù)、功率系數(shù)和功率載荷的平均誤差分別為1.2%、1.1%和1.5%。

    3 數(shù)值模擬

    3.1計(jì)算模型及邊界條件設(shè)置

    本文以ANSYS為仿真平臺(tái)建立微型共軸雙旋翼的有限元模型,主要對(duì)旋翼的旋轉(zhuǎn)域和外部空氣域進(jìn)行幾何建模、網(wǎng)格劃分以及湍流模型和邊界條件的設(shè)置。為提高流場(chǎng)模擬精度,在流場(chǎng)物理量梯度較大的區(qū)域進(jìn)行網(wǎng)格細(xì)化,控制網(wǎng)格總體數(shù)目并保證旋翼模型周圍的網(wǎng)格密度,整個(gè)計(jì)算域的網(wǎng)格總數(shù)為371585。旋翼及流場(chǎng)的網(wǎng)格劃分如圖7所示。

    圖7 共軸雙旋翼及流場(chǎng)的網(wǎng)格劃分Fig.7 Mesh generation of coaxial rotor and flow field

    由于旋翼尺寸較小,雷諾數(shù)范圍接近105,因此流體按不可壓流處理。設(shè)遠(yuǎn)場(chǎng)的流體速度為0,遠(yuǎn)場(chǎng)邊界速度為旋轉(zhuǎn)速度,另外,N-S 方程滿足無(wú)滑移條件,物面邊界速度同樣等于旋翼旋轉(zhuǎn)速度。初始條件為靜止流場(chǎng),采用二階迎風(fēng)格式計(jì)算無(wú)粘通量,采用一階迎風(fēng)格式計(jì)算粘性通量,速度場(chǎng)和壓力場(chǎng)的耦合采取 SIMPLE法。選擇S-A湍流模型,由于壁面用于限制fluid和solid區(qū)域[18],所以取邊界條件為壁面條件。壓力插值選擇Standard,梯度插值選擇Least Squares Cell Based。

    為了驗(yàn)證數(shù)值模擬方法的有效性,本文選取標(biāo)準(zhǔn)翼型NACA0012對(duì)比文獻(xiàn)[15]中的仿真結(jié)果,得到拉力系數(shù)在低雷諾數(shù)環(huán)境下的變化規(guī)律,計(jì)算結(jié)果對(duì)比如圖8所示。從圖中可以看出,在槳尖雷諾數(shù)范圍為0.8×105~1.0×105內(nèi),本文的計(jì)算結(jié)果雖然較文獻(xiàn)略高,但是2種計(jì)算結(jié)果整體吻合較好,證實(shí)了本文的數(shù)值模擬方法是可行的。

    圖8 計(jì)算結(jié)果對(duì)比Fig.8 Comparison of computational results

    3.2流線分析

    共軸雙旋翼在懸停狀態(tài)時(shí),翼間存在強(qiáng)烈的氣動(dòng)干擾。隨著槳葉轉(zhuǎn)動(dòng),后緣產(chǎn)生尾渦相互誘導(dǎo)并漂移,開(kāi)始對(duì)流場(chǎng)氣動(dòng)特性產(chǎn)生影響。為能直觀地從流場(chǎng)結(jié)構(gòu)等細(xì)節(jié)來(lái)觀察共軸雙旋翼周圍及重疊區(qū)域的流場(chǎng)特點(diǎn),各間距下的徑向流線分布如圖9所示。從圖中可明顯看出由于旋翼間的氣動(dòng)干擾,旋翼附近產(chǎn)生的渦流使得流線周向發(fā)散,間距較小時(shí)渦流相互作用強(qiáng)烈,當(dāng)渦流干擾發(fā)生在旋翼附近時(shí)會(huì)對(duì)共軸升力和功耗產(chǎn)生影響。當(dāng)間距增加,流線開(kāi)始耦合在一起時(shí),這種影響開(kāi)始變得特別明顯??傮w來(lái)說(shuō),共軸雙旋翼重疊區(qū)域的流場(chǎng)存在耦合干擾,渦流周向流動(dòng)明顯。隨著間距增大,上下旋翼間的氣動(dòng)干擾減弱。對(duì)比試驗(yàn)結(jié)果,小間距內(nèi)強(qiáng)烈的氣流耦合一定程度上提高了旋翼系統(tǒng)的性能。這是由于當(dāng)旋翼間流場(chǎng)相互作用到一定程度時(shí)這種氣動(dòng)干擾可能會(huì)抵消一部分功耗,這部分減小的功耗使得整體雙旋翼的功率載荷得到提高。此外,渦流撞擊到下旋翼可減小共軸旋翼系統(tǒng)的不穩(wěn)定性。

    3.3槳尖壓強(qiáng)分析

    旋翼槳尖壓強(qiáng)分布如圖10所示,不同間距比下的旋翼槳尖處均存在負(fù)壓區(qū)域,當(dāng)轉(zhuǎn)速較高時(shí)槳尖附近會(huì)產(chǎn)生低壓,從而產(chǎn)生吸附作用,使得上下旋翼槳尖處吸附的槳尖渦脫落并相互誘導(dǎo),最終影響流場(chǎng)分布。此外,對(duì)比旋翼間壓力分布,發(fā)現(xiàn)間距比為0.38的旋翼間壓力較小,與旋翼外部空氣域的壓力差值較大,表明小間距下的共軸雙旋翼與周圍空氣的相互作用更強(qiáng)烈。

    3.4試驗(yàn)與模擬結(jié)果對(duì)比

    為了驗(yàn)證試驗(yàn)測(cè)量的結(jié)果,對(duì)旋翼在轉(zhuǎn)速范圍為1600~2300r/min內(nèi)進(jìn)行數(shù)值模擬。拉力系數(shù)、功率系數(shù)的試驗(yàn)值與計(jì)算值如圖11所示,在轉(zhuǎn)速范圍為1600~1950r/min時(shí),試驗(yàn)值較數(shù)值模擬計(jì)算的拉力系數(shù)略高,功率系數(shù)較計(jì)算值略低,但整體曲線走勢(shì)趨于一致且與文獻(xiàn)[10]的結(jié)果吻合。

    圖11 試驗(yàn)與計(jì)算結(jié)果對(duì)比圖Fig.11 Comparison of experimental and computational results

    4 結(jié) 論

    通過(guò)搭建旋翼氣動(dòng)性能測(cè)試平臺(tái)對(duì)共軸雙旋翼在不同間距下的拉力和功耗進(jìn)行了測(cè)量,分析了其氣動(dòng)特性并計(jì)算了試驗(yàn)誤差,結(jié)合數(shù)值模擬對(duì)比試驗(yàn)結(jié)果分析了旋翼間不同間距下的干擾耦合對(duì)氣動(dòng)特性的影響,并確定了最佳的共軸雙旋翼氣動(dòng)布局,主要結(jié)論如下:

    (1) 共軸拉力隨間距增加,當(dāng)間距達(dá)到一定值時(shí),拉力趨于穩(wěn)定。綜合間距比0.32~0.75下的試驗(yàn)結(jié)果,發(fā)現(xiàn)間距比為0.38的微型共軸雙旋翼系統(tǒng)具有最佳的氣動(dòng)性能,并且轉(zhuǎn)速的提高可以使拉力增幅擴(kuò)大。

    (2) 對(duì)于間距相同的共軸雙旋翼,由于功耗損失小,當(dāng)槳盤(pán)載荷相同時(shí),低轉(zhuǎn)速下的功率載荷較大。

    (3) 計(jì)算所得拉力系數(shù)、功率系數(shù)和功率載荷的誤差均小于2%,系統(tǒng)的測(cè)量誤差隨轉(zhuǎn)速降低而減小,且試驗(yàn)與數(shù)值計(jì)算結(jié)果趨于一致,驗(yàn)證了結(jié)果的有效性。

    (4) 間距比為0.38時(shí),渦流清晰、旋翼間干擾強(qiáng)烈、周向流線槳渦明顯,上下旋翼強(qiáng)烈的氣流耦合提高了旋翼性能。小間距下的旋翼系統(tǒng)穩(wěn)定性更好,旋翼與周圍空氣的作用更強(qiáng)烈。

    [1]唐正飛,王暢,高卓飛.微型旋翼懸停狀態(tài)氣動(dòng)性能分析方法[J].南京航空航天大學(xué)學(xué)報(bào),2011,43(3):357-362.Tang Z F,Wang C,Gao Z F.Analysis method for prediction of micro-rotor performance in hover[J].Journal of Nanjing University of Aeronautics &Astronautics,2011,43(3):357-362.

    [2]朱正,招啟軍,李鵬.懸停狀態(tài)共軸剛性雙旋翼非定常流動(dòng)干擾機(jī)理[J].航空學(xué)報(bào),2016,37 (2):568-578.Zhu Z,Zhao Q J,Li P.Unsteady flow interaction mechanism of coaxial rigid rotors in hover[J].Acta Aeronautica et Astronautica Sinica,2016,37(2):568-578.

    [3]高卓飛.微型共軸雙旋翼氣動(dòng)特性分析方法與試驗(yàn)研究[D].南京:南京航空航天大學(xué),2011.Gao Z F.Analytical methods and experimental research on micro coaxial rotor’s aerodynamic performance[D].Nanjing:Nanjing University of Aeronautics and Astronautics,2011.

    [4]Yoon S,Lee H C,Pulliam T P.Computational study of flow interactions in coaxial rotors[C].The AHS Technical Meeting on Aeromechanics Design for Vertical Lift,San Francisco:NASA Ames Research Center,2016:20-22.

    [5]王強(qiáng),陳銘,王保兵,等.旋翼幾何參數(shù)對(duì)共軸雙旋翼懸停性能的影響分析[J].航空動(dòng)力學(xué)報(bào),2014,29(6):1434-1443.Wang Q,Chen M,Wang B B,et al.Effect of geometric parameters of rotor on hovering performance of coaxial rotor[J].Journal of Aerospace Power,2014,29(6):1434-1443.

    [6]雷瑤.Hex-rotor無(wú)人機(jī)多旋翼流場(chǎng)數(shù)值模擬與試驗(yàn)研究[D].長(zhǎng)春:中國(guó)科學(xué)院長(zhǎng)春光學(xué)精密機(jī)械與物理研究所,2013.Lei Y.Aerodynamics of a hex-rotor SUAV:numerical simulation and experimental study[D].Changchun:Changchun Institute of Optics,Fine Mechanics and Physics Chinese Academy of Sciences,2013.

    [7]王暢.微型共軸雙旋翼氣動(dòng)特性分析方法與實(shí)驗(yàn)研究[D].南京:南京航空航天大學(xué),2010.Wang C.Analytical methods and experimental research of micro-rotor aerodynamic characteristics[D].Nanjing:Nanjing University of Aeronautics and Astronautics,2010.

    [8]葉靚,徐國(guó)華.共軸式雙旋翼懸停流場(chǎng)和氣動(dòng)力的CFD計(jì)算[J].空氣動(dòng)力學(xué)學(xué)報(bào),2012,30(4):437-442.Ye L,Xu G H.Calculation on flow field and aerodynamic force of coaxial rotors in hover with CFD method[J].Acta Aerodynamica Sinica,2012,30(4):437-442.

    [9]Xu H Y,Ye Z Y.Numerical simulation of unsteady flow around forward flight helicopter with coaxial rotors[J].Chinese Journal of Aeronautics,2011,24(1):1-7.

    [10]Bohorquez F.Rotor hover performance and system design of an efficient coaxial rotary wing micro air vehicle[D].Maryland:University of Maryland,2007.

    [11]Syal M.Contributions to the aerodynamic optimization of a coaxial rotor system[D].Maryland:University of Maryland,2008.

    [12]Lei Y,Bai Y,Xu Z J,et al.An experimental investigation on aerodynamic performance of a coaxial rotor system with different rotor spacing and wind speed[J].Experimental Thermal and Fluid Science,2013:779-785.

    [13]Lei Y,Bai Y,Xu Z.Wind effect on aerodynamic optimization for non-planar rotor pairs using full-scale measurements[J].Journal of Intelligent &Robotic Systems (2017):doi:10.1007/s10846-017-0476-y.

    [14]Lakshminarayan V K,Baeder J D.Computational investigation of micro-scale coaxial rotor aerodynamics in hover[J].Journal of Aircraft,2010,47(3):940-955.

    [15]Coleman C P.A survey of theoretical and experimental coaxial rotor aerodynamic research:NASA TP 3675[R].Washington D C:NASA,1997.

    [16]馬楊超,于世美,鄧彥敏.共軸式雙旋翼懸停誘導(dǎo)速度場(chǎng)的PIV實(shí)驗(yàn)研究[J].實(shí)驗(yàn)流體力學(xué),2012,26(1):16-20.Ma Y C,Yu S M,Deng Y M.PIV experimental investigation of coaxial rotors’ induced velocity field in hover[J].Journal of Experiments in Fluid Mechanics,2012,26(1):16-20.

    [17]姬樂(lè)強(qiáng),朱清華,李建波.共軸雙旋翼及孤立旋翼自轉(zhuǎn)氣動(dòng)特性試驗(yàn)研究[J].實(shí)驗(yàn)流體力學(xué),2013,27(5):7-10.Ji L Q,Zhu Q H,Li J B.Experimental investigation on aerodynamic characteristics of auto-rotating coaxial twin rotor and single rotor[J].Journal of Experiments in Fluid Mechanics,2013,27(5):7-10.

    [18]白越,曹萍,高慶嘉,等.六轉(zhuǎn)子微型飛行器及其低雷諾數(shù)下的旋翼氣動(dòng)性能仿真[J].空氣動(dòng)力學(xué)學(xué)報(bào),2011,29(3):325-329.Bai Y,Cao P,Gao Q J,et al.Characterization of a multi rotor micro air vehicle and simulation of airfoil aerodynamics at low Reynolds[J].Acta Aerodynamica Sinica,2011,29(3):325-329.

    雷瑤(1985-),女,重慶人,博士,講師。研究方向:多旋翼無(wú)人機(jī)氣動(dòng)布局優(yōu)化及控制。通信地址:福建省福州市福州地區(qū)大學(xué)新區(qū)學(xué)園路2號(hào)(350116)。E-mail :yaolei@ fzu.edu.cn

    (編輯:楊 娟)

    Numericalsimulationandexperimentalstudyonaerodynamicsofthemicrocoaxialrotors

    Lei Yao*,Ji Yuxia,Wang Changwei

    (School of Mechanical Engineering and Automation,Fuzhou University,Fuzhou 350116,China)

    In order to study the effect of rotor spacing on the aerodynamic performance of small-scale coaxial rotor in hover,a platform is designed to test the aerodynamic characteristics of the coaxial rotor system with different rotor spacing ratiosh/rof 0.32,0.38,0.45,0.51,0.58,0.65 and 0.75.Rotor thrust and power obtained from measurement with different rotary speeds are used to estimate the aerodynamic characteristics of the rotor system.In the meantime,the uncertainty of the thrust coefficient,power coefficient and power loading (PL),and the calculated errors are all less than 2%.Distributions of streamline and pressure with different rotor spacing ratios are obtained from numerical simulation to analyze the interference between top and bottom rotors.Finally,in comparison with test results,the aerodynamic configuration is optimized with the rotor spacing ratio of 0.38.The results prove that the performance of coaxial rotor in hover can be improved with proper rotor spacing.With certain rotor spacing,aerodynamic performance at higher rotational speed is more prominent when the interference between the rotors is beneficial for the aerodynamics,and greater power loading is obtained at lower rotational speed with smaller interference.

    coaxial rotor;in hover;aerodynamic interference;numerical simulation;aerodynamic configuration

    V211.7

    A

    1672-9897(2017)05-0067-07

    10.11729/syltlx20160193

    2016-12-13;

    2017-04-07

    國(guó)家自然科學(xué)基金資助項(xiàng)目(51505087);福建省教育廳資助項(xiàng)目(JA15054)

    *通信作者 E-mail:yaolei@fzu.edu.cn

    LeiY,JiYX,WangCW.Numericalsimulationandexperimentalstudyonaerodynamicsofthemicrocoaxialrotors.JournalofExperimentsinFluidMechanics,2017,31(5):67-73.雷 瑤,紀(jì)玉霞,汪長(zhǎng)煒.微型共軸雙旋翼氣動(dòng)性能數(shù)值模擬與試驗(yàn)分析.實(shí)驗(yàn)流體力學(xué),2017,31(5):67-73.

    猜你喜歡
    共軸旋翼拉力
    改進(jìn)型自抗擾四旋翼無(wú)人機(jī)控制系統(tǒng)設(shè)計(jì)與實(shí)現(xiàn)
    大載重長(zhǎng)航時(shí)油動(dòng)多旋翼無(wú)人機(jī)
    基于CFD的螺旋槳拉力確定方法
    基于STM32的四旋翼飛行器的設(shè)計(jì)
    電子制作(2019年9期)2019-05-30 09:41:48
    共軸剛性雙旋翼非定常氣動(dòng)干擾載荷分析
    自不量力
    共軸共聚焦干涉式表面等離子體顯微成像技術(shù)
    跟蹤導(dǎo)練(三)(3)
    共軸剛性旋翼直升機(jī)旋翼控制相位角問(wèn)題分析
    等動(dòng)拉力和重力拉力蝶泳劃臂動(dòng)作生物力學(xué)分析
    最好的美女福利视频网| 日韩一卡2卡3卡4卡2021年| 中文字幕最新亚洲高清| 变态另类成人亚洲欧美熟女 | 久久精品人人爽人人爽视色| 真人一进一出gif抽搐免费| 国产激情久久老熟女| 又黄又粗又硬又大视频| 人人妻,人人澡人人爽秒播| 亚洲熟妇中文字幕五十中出 | 亚洲成人久久性| 欧美激情 高清一区二区三区| av在线天堂中文字幕 | 国产精品免费视频内射| 亚洲片人在线观看| 18美女黄网站色大片免费观看| 亚洲成国产人片在线观看| 在线国产一区二区在线| 欧美日本亚洲视频在线播放| 一夜夜www| 97超级碰碰碰精品色视频在线观看| 免费一级毛片在线播放高清视频 | 久久久水蜜桃国产精品网| 18禁裸乳无遮挡免费网站照片 | 两性午夜刺激爽爽歪歪视频在线观看 | 欧美日韩中文字幕国产精品一区二区三区 | 日日干狠狠操夜夜爽| 国产成人一区二区三区免费视频网站| 91精品国产国语对白视频| 亚洲人成网站在线播放欧美日韩| 欧美老熟妇乱子伦牲交| 丝袜人妻中文字幕| 黄片播放在线免费| 悠悠久久av| 国产精品一区二区三区四区久久 | 人人妻人人爽人人添夜夜欢视频| 一区在线观看完整版| 最新美女视频免费是黄的| 国产乱人伦免费视频| 精品一区二区三卡| 男女下面插进去视频免费观看| 女人精品久久久久毛片| 又黄又粗又硬又大视频| 日韩欧美一区视频在线观看| 久久久久国内视频| 午夜免费鲁丝| 不卡一级毛片| 成年人免费黄色播放视频| 成年人黄色毛片网站| 9色porny在线观看| 国产高清视频在线播放一区| 久久精品人人爽人人爽视色| 成年女人毛片免费观看观看9| 一本综合久久免费| av在线播放免费不卡| а√天堂www在线а√下载| 久久久国产成人精品二区 | 午夜福利在线观看吧| 欧美中文综合在线视频| 亚洲中文字幕日韩| 亚洲精品国产精品久久久不卡| 真人一进一出gif抽搐免费| 人人妻人人添人人爽欧美一区卜| 大型黄色视频在线免费观看| 欧美激情高清一区二区三区| 日韩欧美在线二视频| 国产精品国产高清国产av| 欧美在线黄色| 久久中文字幕一级| 成人国语在线视频| 日日夜夜操网爽| 亚洲欧美激情在线| 亚洲精品粉嫩美女一区| 首页视频小说图片口味搜索| 超碰成人久久| 一a级毛片在线观看| cao死你这个sao货| 法律面前人人平等表现在哪些方面| 欧美不卡视频在线免费观看 | 久久影院123| 亚洲成人免费电影在线观看| 色综合站精品国产| 久久草成人影院| 亚洲国产精品合色在线| 国产精品免费一区二区三区在线| 99热国产这里只有精品6| 麻豆一二三区av精品| 老熟妇仑乱视频hdxx| 女性生殖器流出的白浆| 老司机午夜福利在线观看视频| 欧美国产精品va在线观看不卡| 国产一区在线观看成人免费| a在线观看视频网站| 国产精品一区二区免费欧美| 国产一区二区三区视频了| 一二三四社区在线视频社区8| 国产欧美日韩一区二区精品| 国产亚洲精品一区二区www| 久久国产乱子伦精品免费另类| 日本a在线网址| 成熟少妇高潮喷水视频| 亚洲一区二区三区色噜噜 | 久久午夜综合久久蜜桃| 天天躁夜夜躁狠狠躁躁| 午夜福利在线免费观看网站| 亚洲国产中文字幕在线视频| 亚洲一区二区三区色噜噜 | 久久热在线av| 天天躁夜夜躁狠狠躁躁| 国内毛片毛片毛片毛片毛片| 欧美一级毛片孕妇| 妹子高潮喷水视频| 亚洲九九香蕉| 男女高潮啪啪啪动态图| 999久久久精品免费观看国产| 国产欧美日韩精品亚洲av| 日本一区二区免费在线视频| 久久久久国产精品人妻aⅴ院| av视频免费观看在线观看| 超色免费av| 免费观看精品视频网站| 天天躁夜夜躁狠狠躁躁| 日本精品一区二区三区蜜桃| 国产麻豆69| 国产aⅴ精品一区二区三区波| 国产精品久久久久成人av| 在线观看一区二区三区| 午夜精品久久久久久毛片777| 麻豆一二三区av精品| 无遮挡黄片免费观看| 亚洲aⅴ乱码一区二区在线播放 | 中文字幕人妻丝袜制服| www.熟女人妻精品国产| 999久久久国产精品视频| 国产成人精品在线电影| 国内久久婷婷六月综合欲色啪| 日韩中文字幕欧美一区二区| 大陆偷拍与自拍| 99热国产这里只有精品6| 久久午夜综合久久蜜桃| 一级a爱视频在线免费观看| 精品乱码久久久久久99久播| 亚洲成a人片在线一区二区| 在线十欧美十亚洲十日本专区| 日韩欧美一区视频在线观看| 久久这里只有精品19| 国产在线精品亚洲第一网站| 国产亚洲精品久久久久久毛片| 久久午夜综合久久蜜桃| 涩涩av久久男人的天堂| 一级片免费观看大全| 91大片在线观看| ponron亚洲| 国内毛片毛片毛片毛片毛片| 精品一区二区三区av网在线观看| 波多野结衣av一区二区av| 中文字幕另类日韩欧美亚洲嫩草| 两个人看的免费小视频| 亚洲成人精品中文字幕电影 | 欧美乱色亚洲激情| www.999成人在线观看| 国产精品秋霞免费鲁丝片| 精品国产国语对白av| 日本欧美视频一区| 成年人黄色毛片网站| 日日摸夜夜添夜夜添小说| 欧美av亚洲av综合av国产av| 欧美日韩视频精品一区| 在线十欧美十亚洲十日本专区| 宅男免费午夜| 人人妻人人澡人人看| 热99国产精品久久久久久7| 欧美午夜高清在线| 国产成人欧美在线观看| 久久精品国产清高在天天线| 少妇裸体淫交视频免费看高清 | 亚洲精品在线观看二区| 中出人妻视频一区二区| 久久香蕉精品热| 国产成人精品在线电影| 中文字幕色久视频| 女人爽到高潮嗷嗷叫在线视频| 成人特级黄色片久久久久久久| 欧美成人免费av一区二区三区| 久久久久久久久中文| 热re99久久精品国产66热6| 少妇被粗大的猛进出69影院| 免费观看人在逋| 精品日产1卡2卡| 精品久久久精品久久久| 波多野结衣av一区二区av| 岛国在线观看网站| 日韩成人在线观看一区二区三区| 在线永久观看黄色视频| 黄频高清免费视频| 亚洲欧美激情综合另类| 亚洲国产欧美网| 每晚都被弄得嗷嗷叫到高潮| 久久精品国产清高在天天线| 亚洲五月天丁香| xxx96com| 日韩欧美国产一区二区入口| 99在线人妻在线中文字幕| 国产欧美日韩一区二区精品| 日韩大尺度精品在线看网址 | 99国产精品一区二区三区| 成人三级黄色视频| 一夜夜www| 婷婷丁香在线五月| 在线十欧美十亚洲十日本专区| 在线视频色国产色| 大型av网站在线播放| svipshipincom国产片| 9热在线视频观看99| 男女下面进入的视频免费午夜 | 最好的美女福利视频网| 国产熟女午夜一区二区三区| 老司机在亚洲福利影院| 国产成人精品在线电影| 免费观看精品视频网站| 十八禁网站免费在线| 婷婷丁香在线五月| 免费看a级黄色片| 久久久久九九精品影院| 欧美日韩精品网址| 亚洲性夜色夜夜综合| 最新美女视频免费是黄的| 无遮挡黄片免费观看| 两性夫妻黄色片| 中文字幕精品免费在线观看视频| 成人免费观看视频高清| 成人18禁高潮啪啪吃奶动态图| 婷婷精品国产亚洲av在线| 麻豆av在线久日| 亚洲av成人一区二区三| 中文字幕人妻丝袜一区二区| 国内久久婷婷六月综合欲色啪| 午夜免费成人在线视频| 69精品国产乱码久久久| 亚洲成人免费电影在线观看| 夜夜爽天天搞| 无限看片的www在线观看| 长腿黑丝高跟| 国产人伦9x9x在线观看| 一区福利在线观看| 久久草成人影院| 国产精品综合久久久久久久免费 | 久99久视频精品免费| 九色亚洲精品在线播放| 高清黄色对白视频在线免费看| 国产精品野战在线观看 | 在线观看免费高清a一片| 国产免费av片在线观看野外av| 性少妇av在线| 三上悠亚av全集在线观看| 夜夜爽天天搞| 在线十欧美十亚洲十日本专区| 亚洲中文av在线| 高清在线国产一区| 亚洲精品久久成人aⅴ小说| 色婷婷av一区二区三区视频| 亚洲全国av大片| 久久久久九九精品影院| 久热爱精品视频在线9| av在线天堂中文字幕 | 他把我摸到了高潮在线观看| 大码成人一级视频| 村上凉子中文字幕在线| 99久久国产精品久久久| 一进一出抽搐gif免费好疼 | 久久久国产欧美日韩av| 国产精品美女特级片免费视频播放器 | 亚洲 欧美一区二区三区| 亚洲精品一区av在线观看| 久久性视频一级片| 亚洲成人免费电影在线观看| 久久精品成人免费网站| 成年版毛片免费区| 日日干狠狠操夜夜爽| 色综合婷婷激情| 国产成年人精品一区二区 | 亚洲熟妇中文字幕五十中出 | 久久久久久久午夜电影 | 亚洲性夜色夜夜综合| 久久青草综合色| 亚洲激情在线av| 亚洲第一欧美日韩一区二区三区| 亚洲成人久久性| 国产有黄有色有爽视频| 黑人操中国人逼视频| 88av欧美| 9色porny在线观看| 男女下面进入的视频免费午夜 | 久久人人97超碰香蕉20202| 久久人妻福利社区极品人妻图片| 亚洲国产欧美网| 日本黄色日本黄色录像| 精品一区二区三区av网在线观看| 一进一出好大好爽视频| 国产高清videossex| 国产一区在线观看成人免费| 在线观看午夜福利视频| 久久精品国产综合久久久| 欧美午夜高清在线| 久久久国产成人免费| 欧美国产精品va在线观看不卡| 免费av毛片视频| 国产国语露脸激情在线看| 两性午夜刺激爽爽歪歪视频在线观看 | 999精品在线视频| 五月开心婷婷网| 老司机午夜福利在线观看视频| 亚洲va日本ⅴa欧美va伊人久久| 婷婷六月久久综合丁香| 欧美精品一区二区免费开放| 天堂影院成人在线观看| 女人被躁到高潮嗷嗷叫费观| 免费高清在线观看日韩| 久久伊人香网站| 身体一侧抽搐| 久久人妻福利社区极品人妻图片| 搡老乐熟女国产| 亚洲专区中文字幕在线| 满18在线观看网站| 人人妻人人添人人爽欧美一区卜| 日韩欧美一区视频在线观看| 久久久久国产精品人妻aⅴ院| 俄罗斯特黄特色一大片| 一区二区三区国产精品乱码| 亚洲av日韩精品久久久久久密| 脱女人内裤的视频| 欧美激情极品国产一区二区三区| 精品国内亚洲2022精品成人| 日本黄色日本黄色录像| 一级毛片精品| 国产成人精品在线电影| 在线观看舔阴道视频| 国产亚洲精品第一综合不卡| 国产亚洲欧美98| 97碰自拍视频| 精品无人区乱码1区二区| 亚洲 国产 在线| 一级作爱视频免费观看| 女人被躁到高潮嗷嗷叫费观| 久99久视频精品免费| 亚洲国产欧美日韩在线播放| 在线观看免费午夜福利视频| 日韩免费av在线播放| 69精品国产乱码久久久| 大型黄色视频在线免费观看| 亚洲人成伊人成综合网2020| 天天躁夜夜躁狠狠躁躁| 亚洲人成电影免费在线| 成人国产一区最新在线观看| av福利片在线| 久久久久久人人人人人| 99精品久久久久人妻精品| 免费在线观看影片大全网站| tocl精华| 欧美日韩亚洲高清精品| 99热国产这里只有精品6| 老熟妇仑乱视频hdxx| av福利片在线| 亚洲精品在线观看二区| aaaaa片日本免费| www.自偷自拍.com| 久久久久久亚洲精品国产蜜桃av| 老鸭窝网址在线观看| 精品福利观看| 亚洲一区二区三区欧美精品| 国产亚洲av高清不卡| 免费看十八禁软件| 亚洲精品成人av观看孕妇| 久久久久国内视频| 99国产综合亚洲精品| 亚洲欧美一区二区三区久久| 黄频高清免费视频| 亚洲色图av天堂| 国产精品 国内视频| 97人妻天天添夜夜摸| 亚洲黑人精品在线| 激情在线观看视频在线高清| 久久人妻av系列| 美女午夜性视频免费| 亚洲国产看品久久| 高清毛片免费观看视频网站 | 90打野战视频偷拍视频| 欧美乱码精品一区二区三区| 91成年电影在线观看| 日日夜夜操网爽| 国产亚洲欧美在线一区二区| 亚洲av成人一区二区三| 亚洲欧美一区二区三区黑人| 亚洲一区中文字幕在线| 久久精品91无色码中文字幕| 一区二区三区激情视频| 久9热在线精品视频| 999久久久国产精品视频| 无人区码免费观看不卡| 无遮挡黄片免费观看| 精品欧美一区二区三区在线| 在线永久观看黄色视频| 日韩大尺度精品在线看网址 | 久久久精品国产亚洲av高清涩受| 午夜免费成人在线视频| 精品日产1卡2卡| 日韩欧美一区视频在线观看| 9色porny在线观看| 五月开心婷婷网| 91在线观看av| 狠狠狠狠99中文字幕| 欧美日韩精品网址| 国产黄色免费在线视频| 超碰97精品在线观看| 久久精品亚洲av国产电影网| 每晚都被弄得嗷嗷叫到高潮| 色婷婷av一区二区三区视频| 99re在线观看精品视频| 亚洲精品国产色婷婷电影| 亚洲av美国av| 91精品国产国语对白视频| 久久久精品欧美日韩精品| 日韩免费高清中文字幕av| 日日干狠狠操夜夜爽| av免费在线观看网站| 一区二区三区精品91| 欧美日韩亚洲高清精品| 久久久久久亚洲精品国产蜜桃av| 久久国产乱子伦精品免费另类| 久久久国产一区二区| 少妇被粗大的猛进出69影院| 看免费av毛片| 丝袜美腿诱惑在线| 老熟妇乱子伦视频在线观看| 18禁美女被吸乳视频| 中文字幕色久视频| 一区福利在线观看| 50天的宝宝边吃奶边哭怎么回事| 亚洲av日韩精品久久久久久密| 亚洲一区二区三区不卡视频| av超薄肉色丝袜交足视频| 亚洲成人免费电影在线观看| 人人妻人人澡人人看| 波多野结衣av一区二区av| 女人精品久久久久毛片| 欧美一区二区精品小视频在线| 18禁裸乳无遮挡免费网站照片 | 亚洲avbb在线观看| 国产91精品成人一区二区三区| 亚洲在线自拍视频| 久久久精品欧美日韩精品| 免费日韩欧美在线观看| 国产黄色免费在线视频| 亚洲专区中文字幕在线| 99久久综合精品五月天人人| 亚洲视频免费观看视频| 久久国产亚洲av麻豆专区| 精品国产乱码久久久久久男人| 亚洲午夜精品一区,二区,三区| 久久精品成人免费网站| 亚洲专区中文字幕在线| 婷婷丁香在线五月| 欧美黄色片欧美黄色片| 国产亚洲精品综合一区在线观看 | 黄色 视频免费看| 性少妇av在线| 80岁老熟妇乱子伦牲交| 一级毛片精品| 日日干狠狠操夜夜爽| 久久香蕉精品热| 国产欧美日韩综合在线一区二区| 国产又爽黄色视频| 人人澡人人妻人| 淫秽高清视频在线观看| 在线观看一区二区三区激情| 成人18禁高潮啪啪吃奶动态图| 国产精品九九99| 老司机深夜福利视频在线观看| 久久精品91无色码中文字幕| 精品电影一区二区在线| 亚洲第一欧美日韩一区二区三区| 亚洲第一青青草原| 丝袜美足系列| 啪啪无遮挡十八禁网站| 老司机靠b影院| 久久久久久亚洲精品国产蜜桃av| 国产欧美日韩一区二区三| 国产精品电影一区二区三区| 黄频高清免费视频| 激情视频va一区二区三区| 女警被强在线播放| 午夜a级毛片| 日韩大尺度精品在线看网址 | 婷婷丁香在线五月| 亚洲久久久国产精品| av电影中文网址| 国产精品久久视频播放| 成人国产一区最新在线观看| 欧美中文综合在线视频| 国产精品九九99| 亚洲国产精品合色在线| 麻豆国产av国片精品| 亚洲精品一区av在线观看| 国内久久婷婷六月综合欲色啪| 国产午夜精品久久久久久| 欧美精品啪啪一区二区三区| 真人一进一出gif抽搐免费| 色在线成人网| 人人澡人人妻人| 久久人妻福利社区极品人妻图片| 亚洲三区欧美一区| av片东京热男人的天堂| 国产欧美日韩一区二区三| 看黄色毛片网站| 欧美在线一区亚洲| 精品电影一区二区在线| 欧洲精品卡2卡3卡4卡5卡区| 成年人免费黄色播放视频| 欧美性长视频在线观看| 成年人黄色毛片网站| 亚洲成人精品中文字幕电影 | 日本 av在线| 人人澡人人妻人| 男女床上黄色一级片免费看| 国产成人av激情在线播放| 99精品欧美一区二区三区四区| 久久人人爽av亚洲精品天堂| 久久精品影院6| 日日摸夜夜添夜夜添小说| 琪琪午夜伦伦电影理论片6080| 久久 成人 亚洲| 99国产精品一区二区蜜桃av| 老熟妇仑乱视频hdxx| 美女扒开内裤让男人捅视频| 久久香蕉国产精品| 水蜜桃什么品种好| 三级毛片av免费| 国产精品久久久人人做人人爽| 久久国产亚洲av麻豆专区| 母亲3免费完整高清在线观看| 好男人电影高清在线观看| svipshipincom国产片| 久久久久久大精品| 侵犯人妻中文字幕一二三四区| 在线观看免费午夜福利视频| 精品人妻在线不人妻| 久久 成人 亚洲| 母亲3免费完整高清在线观看| 亚洲 国产 在线| 亚洲第一欧美日韩一区二区三区| av欧美777| 91大片在线观看| 亚洲精品国产精品久久久不卡| 亚洲国产欧美网| 亚洲精品中文字幕一二三四区| 老熟妇乱子伦视频在线观看| 日韩大码丰满熟妇| 欧美日韩亚洲高清精品| 日本vs欧美在线观看视频| 日本五十路高清| 国产97色在线日韩免费| 色在线成人网| 精品欧美一区二区三区在线| 在线观看免费视频网站a站| 久久久久精品国产欧美久久久| 国产精品免费一区二区三区在线| 50天的宝宝边吃奶边哭怎么回事| 亚洲精品粉嫩美女一区| 9191精品国产免费久久| 日本免费a在线| 色老头精品视频在线观看| 丝袜人妻中文字幕| 欧美日韩亚洲高清精品| 操美女的视频在线观看| 99re在线观看精品视频| 亚洲欧洲精品一区二区精品久久久| 男女床上黄色一级片免费看| 真人做人爱边吃奶动态| 少妇裸体淫交视频免费看高清 | a级毛片黄视频| 久久久久国内视频| 97碰自拍视频| 久久热在线av| 大型av网站在线播放| 黑人巨大精品欧美一区二区蜜桃| 亚洲精品国产色婷婷电影| 中文字幕另类日韩欧美亚洲嫩草| 少妇粗大呻吟视频| 夫妻午夜视频| 成在线人永久免费视频| 身体一侧抽搐| 法律面前人人平等表现在哪些方面| 天天添夜夜摸| 动漫黄色视频在线观看| 日日夜夜操网爽| 18美女黄网站色大片免费观看| 色在线成人网| 中文字幕av电影在线播放| 久久久久久久久久久久大奶| 免费在线观看影片大全网站| 日本vs欧美在线观看视频| 九色亚洲精品在线播放| 国产精品免费一区二区三区在线| 午夜成年电影在线免费观看| 18禁国产床啪视频网站| 日日夜夜操网爽| av在线天堂中文字幕 | 欧美日本亚洲视频在线播放| 怎么达到女性高潮| 亚洲人成伊人成综合网2020| 午夜福利,免费看| 欧美日韩精品网址| 18禁裸乳无遮挡免费网站照片 | 久久久久国产精品人妻aⅴ院| 久久久精品国产亚洲av高清涩受|