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

    第二代低階面元法耦合自由尾跡計(jì)算旋翼槳葉氣動載荷

    2012-04-06 12:48:40王海鵬王華明
    關(guān)鍵詞:尾跡低階元法

    王海鵬,王華明

    (南京航空航天大學(xué) 直升機(jī)旋翼動力學(xué)重點(diǎn)實(shí)驗(yàn)室,江蘇 南京 210016)

    0 引 言

    拉普拉斯方程是分析不可壓、無粘、無旋流場的基本方程,用拉普拉斯方程計(jì)算飛行器的氣動載荷,最常用的數(shù)值方法是面元法。面元法經(jīng)過半個世紀(jì)的發(fā)展,由最初的求解流場中無升力體表面壓力分布到現(xiàn)在的全機(jī)氣動特性計(jì)算,其精度和穩(wěn)定性已得到大幅度的提升[1]。20世紀(jì)80年代末,美國NASA推出了根據(jù)低階面元法開發(fā)的計(jì)算程序“VSAERO”[2],該程序通過引入內(nèi)部速度勢邊界條件,在減少計(jì)算量的同時較大地提高了計(jì)算精度,受到世界航空界的廣泛關(guān)注,被認(rèn)為是第二代低階面元法發(fā)展的起始。

    在面元法發(fā)展的同時,旋翼自由尾跡計(jì)算方法也在不斷進(jìn)步,由最初的影響系數(shù)法[3]、時間步進(jìn)松弛迭代法[4]、到現(xiàn)在的半隱式預(yù)估修松弛迭代法[5],其計(jì)算穩(wěn)定性與精度也在逐步提高。由于研究側(cè)重點(diǎn)不同,國內(nèi)外旋翼自由尾跡計(jì)算往往集中研究尾跡形狀的求解,對于旋翼槳葉多采用簡單的二階升力線[5-6]或不考慮槳葉厚度的升力面模型[4],這種處理辦法不能完全模擬槳葉的三維幾何特性,也不能計(jì)算槳葉上下表面的壓力分布與速度分布。文獻(xiàn)[7-9]采用低階面元法耦合自由尾跡計(jì)算旋翼氣動特性,但其尾跡計(jì)算均采用啟動過程的時間步進(jìn)法,穩(wěn)定性較差,收斂速度慢。隨著計(jì)算機(jī)技術(shù)的發(fā)展,雷諾時均N-S方程(RANS)也被運(yùn)用于求解旋翼的氣動特性,但由于旋翼氣動環(huán)境的特殊性與計(jì)算方法本身的特點(diǎn),計(jì)算仍需耗費(fèi)大量的時間與資源。

    本文采用第二代低階面元法耦合自由尾跡預(yù)估修正法計(jì)算旋翼在軸流狀態(tài)下的氣動載荷。在計(jì)算中考慮槳葉的三維幾何特性,采用矩形、三角形單位平面重構(gòu)槳葉幾何外形;槳葉尾跡采用全展尾跡,其中包括三圈自由尾跡(近尾跡)和四圈遠(yuǎn)尾跡;并利用槳葉翼型風(fēng)洞試驗(yàn)數(shù)據(jù)對氣動載荷進(jìn)行粘性修正。通過和CFD方法進(jìn)行對比分析,證明了本文所采用的計(jì)算方法的優(yōu)點(diǎn)。

    1 計(jì)算方法

    1.1 第二代低階面元法

    面元法總控方程為拉普拉斯方程,其為橢圓型偏微分方程。對于三維拉普拉斯方程常利用格林定理,將其轉(zhuǎn)化為升力體表面的二維問題[10]。

    第一代低階面元法通常只在升力體表面布置一種奇點(diǎn),或在給定某種奇點(diǎn)分布下布置另一種未知強(qiáng)度奇點(diǎn),通過物面不穿透邊界條件求解未知奇點(diǎn)強(qiáng)度分布。上述計(jì)算方法往往在非控制點(diǎn)處存在法向速度“穿透”,對計(jì)算結(jié)果有較大影響。只有通過加密升力體便面的網(wǎng)格數(shù)目,才能獲得較好的計(jì)算結(jié)果。

    第二代低階面元法采用矩形、三角形單位平面重構(gòu)槳葉幾何外形,槳葉面元上布置未知強(qiáng)度的面源與偶極子,槳葉后緣脫出偶極子尾渦(全展),槳葉和尾跡模型如圖1所示。通過上述源項(xiàng)組合,空間中任意點(diǎn)P處速度勢可表示為以下形式[10]。

    其中s為物面坐標(biāo)點(diǎn)矢量;R(s,p)為物面或尾跡點(diǎn)至計(jì)算點(diǎn)的距離;n為物面點(diǎn)處垂直于物面指向流域的法向量;Φ∞p為遠(yuǎn)方自由流在計(jì)算點(diǎn)P點(diǎn)處的速度勢;S為槳葉表面;Sw為槳葉尾跡表面,δ(s)與μ(s)分別為布置于槳葉表面的源項(xiàng)與偶極子項(xiàng)。

    圖1 面元法的槳葉模型與尾跡模型Fig.1 Blade panels and its wake model

    上述方程采用混合邊界條件求解,具體即為槳葉外表面速度矢量不穿透(紐曼邊界條件)和槳葉內(nèi)表面速度勢為常值(狄利克雷邊界條件)。如下式所示:

    新邊界條件的引入使得槳葉表面的源強(qiáng)可直接顯式求出而無需方程求解,極大地減少了計(jì)算量;內(nèi)部速度勢的引入減少面元間法向速度的“泄露”,降低了面網(wǎng)格分布對計(jì)算結(jié)果的敏感性,從而提高了計(jì)算精度。

    尾跡強(qiáng)度使用后緣庫塔條件確定,其等于后緣上下面元偶極子強(qiáng)度之差。在已知尾跡形狀后,利用邊界條件并在槳葉表面、尾跡處劃分單位面元可將2式離散為線性方程組,求解此線性方程組可得未知源強(qiáng)與偶極子強(qiáng)度,進(jìn)而反推出槳葉表面速度與槳葉氣動載荷。

    1.2 自由尾跡計(jì)算

    利用尾跡不受力條件,得到尾跡總控方程[6]:

    上式為一常微分方程,其中p為尾跡點(diǎn)坐標(biāo)矢量,r為尾跡渦線,V(p)為在空間P點(diǎn)處的速度矢量合,其中包括槳葉面元,尾跡偶極子面(尾跡渦線)和遠(yuǎn)方來流的共同作用。本文采用半隱式預(yù)估修正法[6]求解尾跡節(jié)點(diǎn)位置。

    在自由尾跡計(jì)算中,渦核模型的選取尤為重要,其直接決定尾跡計(jì)算的收斂。本文使用文獻(xiàn)[6,11]中推薦使用的渦核模型,渦核半徑表達(dá)式如下式所示。

    式中rc為渦核半徑;rc0為初始渦核半徑,本文中選取0.1倍槳葉剖面弦長;a為 “Oseen”系數(shù)取1.225643;Rev為渦雷諾數(shù);v為動粘性系數(shù);t為尾跡歷程時間;a1值在0.01至0.1之間選取[5]。

    在自由尾跡計(jì)算中,評判兩次迭代的尾跡均方根差值來判斷尾跡形狀是否收斂[5]。尾跡均方根差值計(jì)算如下式所示:

    1.3 計(jì)算步驟與迭代方法

    計(jì)算時首先劃分槳葉面元,根據(jù)動量葉素理論估算平均誘導(dǎo)速度,生成初始尾跡,并劃分尾跡面元;在此基礎(chǔ)上根據(jù)第二代低階面元法求解槳葉面元奇點(diǎn)強(qiáng)度,并求出該尾跡下槳葉氣動載荷;接下來根據(jù)槳葉與尾跡奇點(diǎn)強(qiáng)度,運(yùn)用半隱式預(yù)估修正法計(jì)算旋翼新尾跡;繼而判斷尾跡形狀是否收斂,若滿足收斂標(biāo)準(zhǔn)便結(jié)束計(jì)算,未滿足則返回第二步繼續(xù)計(jì)算。詳細(xì)步驟如圖2所示。

    圖2 計(jì)算流程框圖Fig.2 Calculation flow chart

    1.4 粘性修正

    由于基于拉普拉斯方程的面元法無法考慮粘性力,在計(jì)算旋翼旋轉(zhuǎn)力矩時,會造成較大誤差。針對上述情況,通常采用“附面層修正”方法來考慮粘性作用影響,該方法引入有粘/無粘耦合迭代,用表面源模型模擬粘性帶來的附面層位移厚度效應(yīng)。該方法可較為準(zhǔn)確地計(jì)算粘性作用,但增加了面元法的計(jì)算量與計(jì)算耗時。考慮到旋翼工作時,槳葉翼型剖面攻角較小,可忽略粘性力對翼型法向載荷的影響,本文即利用面元法所計(jì)算的法向氣動載荷,根據(jù)相同雷諾數(shù)下槳葉剖面的二維翼型CFD計(jì)算數(shù)據(jù)或風(fēng)洞試驗(yàn)數(shù)據(jù),查表反推出對應(yīng)法向力系數(shù)下的軸向力系數(shù),用于修正粘性影響。此法也可快速估算槳葉不同半徑處攻角。

    2 計(jì)算結(jié)果對比

    為驗(yàn)證上述計(jì)算方法的正確性,本文選用廣泛引用的“Caradonna”旋翼作為計(jì)算算例,其翼型為NACA0012,展弦比為6,無扭轉(zhuǎn),無尖削,其他參數(shù)參見文獻(xiàn)[12]。

    在本算例計(jì)算中,對單片槳葉劃分672個面元:其中槳葉沿展向劃分10段;沿弦向上下表面劃分48段;沿槳葉厚度方向劃分4段。每片槳葉拖出3圈自由尾跡,4圈遠(yuǎn)尾跡,單周尾跡中等角度劃分24個尾跡節(jié)點(diǎn),共計(jì)1848個尾跡節(jié)點(diǎn)。

    為便于對比,本文還采用CFD方法求解雷諾時均N-S方程計(jì)算該旋翼氣動特性。在CFD計(jì)算中采用全結(jié)構(gòu)化網(wǎng)格,網(wǎng)格數(shù)量約200萬,采用剪切應(yīng)力輸運(yùn)(SST)湍流模型。

    圖3為以上兩種方法計(jì)算槳葉壓力分布與試驗(yàn)值的對比。從對比結(jié)果可以看出,以上兩種方法與試驗(yàn)結(jié)果有很好的一致性。

    圖4為以上兩種方法計(jì)算旋翼在1250r/min時,對應(yīng)5°、8°和12°總距角時旋翼拉力系數(shù)。通過對比可以發(fā)現(xiàn),采用面元法耦合自由尾跡和CFD的計(jì)算結(jié)果與試驗(yàn)結(jié)果有較好的一致性。

    圖3 槳葉剖面壓力分布對比Fig.3 The comparison chart of pressure distribution on blade sections

    圖4 旋翼拉力系數(shù)對比Fig.4 The comparison chart of rotor lift coefficient

    圖5為以上兩種方法計(jì)算旋翼在1250r/min時,對應(yīng)5°、8°和12°總距角時旋翼力矩系數(shù)。其中包括面元法修正前與修正后的計(jì)算結(jié)果。由于面元法只能考慮旋翼力矩中由槳葉誘導(dǎo)阻力造成的旋轉(zhuǎn)力矩,無法考慮由槳葉摩擦造所成的旋轉(zhuǎn)力矩,無摩擦修正的計(jì)算結(jié)果明顯小于CFD計(jì)算結(jié)果。利用槳葉翼型風(fēng)洞試驗(yàn)數(shù)據(jù)修正后,槳葉力矩計(jì)算結(jié)果與CFD計(jì)算結(jié)果有較好的一致性。

    圖5 旋翼力矩系數(shù)對比Fig.5 The comparison chart of rotor moment coefficient

    圖6為上述兩種方法計(jì)算拉力系數(shù)與試驗(yàn)值的相對誤差對比圖。從圖中可以看出,兩種方法計(jì)算精度相當(dāng)。

    圖7為上述兩種計(jì)算方法耗時對比圖,計(jì)算平臺處理器為“Xeon QC E7330”。通過對比,可見本文方所述面元法耗時僅為CFD計(jì)算耗時的1/10左右,大大節(jié)省了計(jì)算時間。

    圖6 拉力系數(shù)相對誤差對比圖Fig.6 The comparison chart of lift coefficient relative error

    圖7 計(jì)算耗時對比圖Fig.7 The comparison chart of calculation time cost

    由于本文中尾跡模型為全展自由尾跡模型,其為從槳葉后緣拖出的渦面,而非一條集中渦線,直接比較尾跡形狀有一定困難。圖8給出了1250r/min時不同總距角的尾跡計(jì)算結(jié)果,通過對比可以發(fā)現(xiàn)隨著總距角的增加,尾跡徑向收縮更為劇烈,軸向移動速度也有所增加,這符合一般計(jì)算規(guī)律。

    圖8 旋翼尾跡形狀(只顯示3圈自由尾跡)Fig.8 The figure of rotor wake geometry

    3 結(jié)論與展望

    第二代低階面元法耦合自由尾跡的計(jì)算方法能夠模擬旋翼槳葉的三維幾何特性,快速計(jì)算槳葉表面的壓力分布和旋翼的氣動載荷。

    本文提出的的方法計(jì)算旋翼在軸流狀態(tài)下的氣動特性與CFD方法的計(jì)算精度相仿,但計(jì)算耗費(fèi)的時間僅為CFD方法的十分之一。

    通過引入槳葉的揮舞運(yùn)動,本文提出的計(jì)算方法可拓展至直升機(jī)前飛時旋翼的氣動特性計(jì)算。本文提出的計(jì)算方法和旋翼氣動噪聲分析方法相結(jié)合,可以預(yù)估旋翼的氣動噪聲。

    [1]ERICKSON L L.Panel methods an introduction[R].NASA TP 2995,1990.

    [2]MASKEW B.Program VSAERO theory document[R].NASA CR 4023,1987.

    [3]BLISS D B.A new approach to the free wake problem for hovering rotors[A].Annual Forum Proceedings of the American Helicopter Society[C].1985,463-477.

    [4]徐國華.應(yīng)用自由尾跡分析的新型槳尖旋翼氣動特性研究[D].南京:南京航空航天大學(xué),1990.

    [5]BAGAI A.Contributions to the mathematical modeling of rotor flow-fields using a Pseudo-Implicit Free-Wake A-nalysis[D].College Park:University of Maryland,1995.

    [6]曾洪江,胡繼忠.一種新的自由渦尾跡計(jì)算方法[J].航空學(xué)報(bào),2004,25(6):546-550.

    [7]尹堅(jiān)平,胡章偉.旋翼表面非定常壓力脈動計(jì)算的三維自由尾跡非定常面元法[J].空氣動力學(xué)學(xué)報(bào),1997,15(2):185-191.

    [8]AHMED S R,VIDJAJA V T.Unsteady panel method calculation of pressure distribution on BO105model rotor blades[J].Journal of the American Helicopter Society,1998,43(1):47-56.

    [9]仲唯貴.直升機(jī)旋翼尾渦與尾槳干擾噪聲研究[D].南京:南京航空航天大學(xué),2006.

    [10]KATZ J,PLOTKIN A.Low-speed aerodynamics,from wing theory to panel methods[M]. McGraw-Hill,1991.

    [11]BHAGWAT M J,LEISHMAN J G.Generalized viscous vortex core models for application to free-vortex wake and aeroacoustic calculations[A].Proceedings of the 58th Annual Forum of AHS International[C].2002.

    [12]CARADONNA F X,TUNG C.Experimental and analytical studies of a model helicopter rotor in hover[R].NASA TM 81232,1981.

    猜你喜歡
    尾跡低階元法
    一種基于Radon 變換和尾跡模型的尾跡檢測算法
    換元法在解題中的運(yùn)用
    山西低階煤分布特征分析和開發(fā)利用前景
    基于離散元法的礦石對溜槽沖擊力的模擬研究
    一類具低階項(xiàng)和退化強(qiáng)制的橢圓方程的有界弱解
    Extended Fisher-Kolmogorov方程的一類低階非協(xié)調(diào)混合有限元方法
    基于EEMD-Hilbert譜的渦街流量計(jì)尾跡振蕩特性
    換元法在解題中的應(yīng)用
    “微元法”在含電容器電路中的應(yīng)用
    國內(nèi)外低階煤煤層氣開發(fā)現(xiàn)狀和我國開發(fā)潛力研究
    中國煤層氣(2015年3期)2015-08-22 03:08:23
    亚洲专区中文字幕在线| 女人被狂操c到高潮| 日韩高清综合在线| 欧美丝袜亚洲另类 | 无遮挡黄片免费观看| 午夜a级毛片| 美女高潮喷水抽搐中文字幕| 亚洲黑人精品在线| 午夜激情欧美在线| 国产亚洲精品久久久com| 日韩高清综合在线| 性欧美人与动物交配| 国产 一区 欧美 日韩| 欧美日本视频| 午夜福利在线观看免费完整高清在 | 身体一侧抽搐| 我的女老师完整版在线观看| 夜夜爽天天搞| 精品久久久久久久久久久久久| 亚洲精华国产精华液的使用体验 | 禁无遮挡网站| 99久国产av精品| 啦啦啦观看免费观看视频高清| 国产探花极品一区二区| 久久久久久九九精品二区国产| 国产一区二区在线观看日韩| 听说在线观看完整版免费高清| 俺也久久电影网| 国产亚洲欧美98| 国产aⅴ精品一区二区三区波| 色综合婷婷激情| 欧美高清成人免费视频www| 免费av不卡在线播放| 午夜福利高清视频| 国产精品国产三级国产av玫瑰| 男女之事视频高清在线观看| 少妇人妻精品综合一区二区 | 久久久久久久午夜电影| 毛片女人毛片| 成人美女网站在线观看视频| 欧美区成人在线视频| 热99在线观看视频| 日韩av在线大香蕉| 久久草成人影院| 欧美潮喷喷水| 亚洲精品在线观看二区| 亚洲av日韩精品久久久久久密| 成人无遮挡网站| 欧美3d第一页| 尾随美女入室| 日本五十路高清| 美女高潮喷水抽搐中文字幕| 亚洲黑人精品在线| xxxwww97欧美| 一本精品99久久精品77| 亚洲专区国产一区二区| 欧美+亚洲+日韩+国产| 他把我摸到了高潮在线观看| 欧美激情在线99| 国内精品久久久久久久电影| 亚洲av中文字字幕乱码综合| 日本-黄色视频高清免费观看| 成人国产麻豆网| 欧美性猛交╳xxx乱大交人| 日韩 亚洲 欧美在线| 国产 一区 欧美 日韩| 99久久九九国产精品国产免费| 18禁黄网站禁片免费观看直播| 成年女人永久免费观看视频| 国产人妻一区二区三区在| 99久久成人亚洲精品观看| 久久精品国产亚洲av涩爱 | 亚洲,欧美,日韩| 亚洲真实伦在线观看| 国产精品一区二区三区四区久久| 亚洲国产色片| 日日摸夜夜添夜夜添小说| 在线天堂最新版资源| 一a级毛片在线观看| 美女xxoo啪啪120秒动态图| 久久久色成人| 黄色欧美视频在线观看| 变态另类成人亚洲欧美熟女| 久久99热这里只有精品18| av在线天堂中文字幕| av女优亚洲男人天堂| 国产综合懂色| 免费av不卡在线播放| ponron亚洲| 91麻豆精品激情在线观看国产| 99久久精品国产国产毛片| 欧洲精品卡2卡3卡4卡5卡区| 国产乱人伦免费视频| av中文乱码字幕在线| 麻豆国产97在线/欧美| 999久久久精品免费观看国产| 久久精品国产清高在天天线| www.www免费av| x7x7x7水蜜桃| 中文字幕高清在线视频| 午夜福利在线观看免费完整高清在 | 最近在线观看免费完整版| 国产精品永久免费网站| 麻豆成人av在线观看| 3wmmmm亚洲av在线观看| 99九九线精品视频在线观看视频| 精品人妻熟女av久视频| 欧美潮喷喷水| 国产一级毛片七仙女欲春2| 日本 av在线| 国产伦人伦偷精品视频| 国产v大片淫在线免费观看| 亚洲精华国产精华液的使用体验 | 男女那种视频在线观看| 国产精品国产高清国产av| 亚洲真实伦在线观看| 色综合站精品国产| 99热6这里只有精品| 精品久久久噜噜| 中文亚洲av片在线观看爽| 一区福利在线观看| 国产高清有码在线观看视频| 一卡2卡三卡四卡精品乱码亚洲| 日韩国内少妇激情av| 人人妻,人人澡人人爽秒播| 国产成人福利小说| bbb黄色大片| 中文字幕久久专区| 国产女主播在线喷水免费视频网站 | 国产精品乱码一区二三区的特点| 神马国产精品三级电影在线观看| 成年女人永久免费观看视频| 精品无人区乱码1区二区| 国产伦精品一区二区三区四那| 一区二区三区激情视频| 在线免费观看的www视频| 亚洲va在线va天堂va国产| 亚洲图色成人| 黄色一级大片看看| 成人国产麻豆网| 欧美在线一区亚洲| 熟女人妻精品中文字幕| 日韩大尺度精品在线看网址| 日韩精品青青久久久久久| 色哟哟·www| 一本一本综合久久| 可以在线观看毛片的网站| 精品一区二区三区人妻视频| 日韩欧美三级三区| 亚洲五月天丁香| 极品教师在线视频| 黄色日韩在线| 88av欧美| 国产视频内射| 少妇高潮的动态图| 日本 av在线| 久久久精品欧美日韩精品| 久久国内精品自在自线图片| 久久久久久久久久久丰满 | 亚洲av免费在线观看| 在线观看午夜福利视频| 最近最新中文字幕大全电影3| 又紧又爽又黄一区二区| 一进一出好大好爽视频| 成年人黄色毛片网站| 成人美女网站在线观看视频| 日日啪夜夜撸| 亚洲狠狠婷婷综合久久图片| 美女高潮的动态| 亚洲,欧美,日韩| 91av网一区二区| 日日干狠狠操夜夜爽| 女的被弄到高潮叫床怎么办 | 22中文网久久字幕| 欧美成人a在线观看| 最新在线观看一区二区三区| 欧美激情久久久久久爽电影| 免费电影在线观看免费观看| 久99久视频精品免费| 国产大屁股一区二区在线视频| 变态另类成人亚洲欧美熟女| 床上黄色一级片| 久99久视频精品免费| 成人综合一区亚洲| 日韩欧美 国产精品| 九色成人免费人妻av| 国产单亲对白刺激| 又紧又爽又黄一区二区| 欧美精品国产亚洲| 99热6这里只有精品| 国产午夜精品论理片| 麻豆精品久久久久久蜜桃| 99久久精品一区二区三区| 国产三级中文精品| 亚洲七黄色美女视频| 中文字幕免费在线视频6| 亚洲经典国产精华液单| 香蕉av资源在线| 禁无遮挡网站| 在线观看舔阴道视频| 日本色播在线视频| 国产精华一区二区三区| 少妇的逼水好多| 国产三级中文精品| 91久久精品国产一区二区三区| 韩国av在线不卡| 欧美激情国产日韩精品一区| 精华霜和精华液先用哪个| 色播亚洲综合网| av黄色大香蕉| 看片在线看免费视频| 麻豆成人av在线观看| 天堂动漫精品| 亚洲成a人片在线一区二区| 国产精品不卡视频一区二区| 五月伊人婷婷丁香| 男插女下体视频免费在线播放| 中文字幕久久专区| 无遮挡黄片免费观看| av视频在线观看入口| 国产av一区在线观看免费| 99视频精品全部免费 在线| 国产精品一及| av女优亚洲男人天堂| 亚洲精华国产精华液的使用体验 | 99riav亚洲国产免费| 午夜精品一区二区三区免费看| 91午夜精品亚洲一区二区三区 | 国内精品久久久久精免费| 国产高清视频在线播放一区| 男女视频在线观看网站免费| 国产在视频线在精品| 日本黄色视频三级网站网址| 久久精品影院6| 人妻少妇偷人精品九色| 日韩,欧美,国产一区二区三区 | 热99在线观看视频| 亚洲欧美清纯卡通| 91狼人影院| 久久亚洲精品不卡| 免费观看的影片在线观看| 黄色欧美视频在线观看| 老熟妇仑乱视频hdxx| 麻豆国产97在线/欧美| 亚洲精品乱码久久久v下载方式| 亚洲精品亚洲一区二区| 亚洲四区av| 一个人看视频在线观看www免费| 国产午夜精品久久久久久一区二区三区 | 久久午夜福利片| 国产视频内射| 日韩中文字幕欧美一区二区| 久久亚洲真实| 熟妇人妻久久中文字幕3abv| 嫩草影院精品99| 久久人妻av系列| 亚洲国产精品sss在线观看| 久久九九热精品免费| 欧美人与善性xxx| 国产欧美日韩精品亚洲av| 欧美极品一区二区三区四区| 精品久久久久久,| 欧美色视频一区免费| 国产白丝娇喘喷水9色精品| 国产一区二区激情短视频| 欧洲精品卡2卡3卡4卡5卡区| 五月伊人婷婷丁香| 干丝袜人妻中文字幕| 永久网站在线| 久久亚洲精品不卡| 精品无人区乱码1区二区| 日本-黄色视频高清免费观看| 亚洲无线在线观看| 啦啦啦韩国在线观看视频| av天堂在线播放| 给我免费播放毛片高清在线观看| 亚洲va日本ⅴa欧美va伊人久久| 日韩中字成人| 99热这里只有是精品在线观看| 波野结衣二区三区在线| 综合色av麻豆| 日韩一区二区视频免费看| 亚洲国产精品sss在线观看| 亚洲欧美清纯卡通| 国产精品乱码一区二三区的特点| 88av欧美| 国产日本99.免费观看| 永久网站在线| 嫩草影院新地址| 又黄又爽又刺激的免费视频.| 1000部很黄的大片| 色精品久久人妻99蜜桃| 成人一区二区视频在线观看| 国产主播在线观看一区二区| 99热这里只有是精品在线观看| 国产女主播在线喷水免费视频网站 | 成人国产综合亚洲| 女生性感内裤真人,穿戴方法视频| 国产综合懂色| 欧美日韩精品成人综合77777| 亚洲av免费高清在线观看| 熟妇人妻久久中文字幕3abv| 男女边吃奶边做爰视频| 99热这里只有精品一区| 国产精品福利在线免费观看| 麻豆成人午夜福利视频| 婷婷精品国产亚洲av| 午夜亚洲福利在线播放| 国产成人a区在线观看| av天堂中文字幕网| 91麻豆精品激情在线观看国产| 国产成人aa在线观看| 伊人久久精品亚洲午夜| 99在线人妻在线中文字幕| 无人区码免费观看不卡| 亚洲第一区二区三区不卡| 国产一区二区在线av高清观看| 亚洲欧美精品综合久久99| 日韩欧美 国产精品| 99国产精品一区二区蜜桃av| av.在线天堂| 日日干狠狠操夜夜爽| 亚洲无线观看免费| 国产精品精品国产色婷婷| 天堂网av新在线| www.色视频.com| 欧美色欧美亚洲另类二区| 国产一区二区三区视频了| 99国产精品一区二区蜜桃av| 男人的好看免费观看在线视频| 中文在线观看免费www的网站| 搞女人的毛片| 91久久精品国产一区二区成人| 久久久久性生活片| 国内精品久久久久精免费| 国产亚洲精品久久久久久毛片| 精华霜和精华液先用哪个| 国产欧美日韩精品亚洲av| 三级毛片av免费| av在线亚洲专区| 国内精品美女久久久久久| 免费高清视频大片| 麻豆成人午夜福利视频| 久久久久久久久中文| 免费看a级黄色片| 俺也久久电影网| 中文资源天堂在线| 色综合亚洲欧美另类图片| 在线免费观看不下载黄p国产 | 国产精品一及| 99热这里只有是精品在线观看| 麻豆成人av在线观看| 亚洲最大成人av| 国产精品久久电影中文字幕| 成人国产一区最新在线观看| 搡老妇女老女人老熟妇| 校园春色视频在线观看| 国产亚洲精品久久久久久毛片| 久久精品综合一区二区三区| 男人狂女人下面高潮的视频| 国产白丝娇喘喷水9色精品| 他把我摸到了高潮在线观看| 成人性生交大片免费视频hd| 亚洲最大成人中文| 88av欧美| 极品教师在线视频| 大型黄色视频在线免费观看| 男女之事视频高清在线观看| 最近在线观看免费完整版| 色尼玛亚洲综合影院| 亚洲av中文字字幕乱码综合| 亚洲美女黄片视频| 精品免费久久久久久久清纯| 夜夜夜夜夜久久久久| 内射极品少妇av片p| 国产高清视频在线播放一区| 偷拍熟女少妇极品色| 美女高潮的动态| 国产精品国产高清国产av| 久久这里只有精品中国| 男女边吃奶边做爰视频| 亚洲人成伊人成综合网2020| 看黄色毛片网站| 久久久久国产精品人妻aⅴ院| 精品人妻视频免费看| 中国美女看黄片| 免费在线观看日本一区| 一区二区三区激情视频| 丰满的人妻完整版| 免费搜索国产男女视频| 成人亚洲精品av一区二区| 91精品国产九色| 两性午夜刺激爽爽歪歪视频在线观看| 在线观看av片永久免费下载| 国产精品爽爽va在线观看网站| 在线免费十八禁| 亚洲av中文字字幕乱码综合| 国产精品无大码| 天堂动漫精品| 日韩 亚洲 欧美在线| 久久久久久久久久成人| 午夜福利欧美成人| 一边摸一边抽搐一进一小说| 精品久久久噜噜| 婷婷精品国产亚洲av在线| 狠狠狠狠99中文字幕| 老师上课跳d突然被开到最大视频| 成人毛片a级毛片在线播放| 亚洲性夜色夜夜综合| 国产精品永久免费网站| 日韩人妻高清精品专区| 嫩草影院新地址| 国产在视频线在精品| 伊人久久精品亚洲午夜| 婷婷色综合大香蕉| 日韩av在线大香蕉| 在线播放国产精品三级| 少妇高潮的动态图| 久久精品国产自在天天线| 禁无遮挡网站| 久久久久久久亚洲中文字幕| 老熟妇乱子伦视频在线观看| 99热网站在线观看| 日本 欧美在线| 精品日产1卡2卡| 在线a可以看的网站| 欧美精品国产亚洲| av女优亚洲男人天堂| 国产伦在线观看视频一区| 亚洲自偷自拍三级| 中文在线观看免费www的网站| 婷婷六月久久综合丁香| 久久精品国产亚洲av天美| 久久久午夜欧美精品| 色视频www国产| 黄色日韩在线| 嫩草影院新地址| 国产伦精品一区二区三区四那| 亚洲无线在线观看| 精品久久久久久久人妻蜜臀av| 久久久久久久久久黄片| 免费一级毛片在线播放高清视频| 精品国产三级普通话版| 欧美成人免费av一区二区三区| 国产av在哪里看| 国产精品嫩草影院av在线观看 | 精品久久久久久,| 成年免费大片在线观看| 村上凉子中文字幕在线| 色尼玛亚洲综合影院| 国产主播在线观看一区二区| 日本三级黄在线观看| 久久人人精品亚洲av| 韩国av在线不卡| 神马国产精品三级电影在线观看| 午夜福利视频1000在线观看| 亚洲电影在线观看av| 国产亚洲av嫩草精品影院| 亚洲电影在线观看av| 成人av在线播放网站| 国产精品综合久久久久久久免费| 国产成人aa在线观看| 成人av一区二区三区在线看| 中文字幕精品亚洲无线码一区| 97碰自拍视频| 熟女人妻精品中文字幕| 精品人妻视频免费看| 国产成人一区二区在线| 国产黄片美女视频| 欧美精品国产亚洲| 久久久久性生活片| 网址你懂的国产日韩在线| 黄色女人牲交| 一进一出抽搐动态| 天天一区二区日本电影三级| 免费av毛片视频| 国产欧美日韩精品亚洲av| 他把我摸到了高潮在线观看| 18禁黄网站禁片午夜丰满| 麻豆成人av在线观看| 观看美女的网站| 成人鲁丝片一二三区免费| 欧美+亚洲+日韩+国产| 免费观看的影片在线观看| 成人精品一区二区免费| 久久这里只有精品中国| 一级黄片播放器| 成人高潮视频无遮挡免费网站| 国产高清激情床上av| 一本一本综合久久| 18+在线观看网站| 日韩欧美一区二区三区在线观看| 国内精品久久久久久久电影| h日本视频在线播放| 丰满人妻一区二区三区视频av| 麻豆精品久久久久久蜜桃| 日韩国内少妇激情av| 国国产精品蜜臀av免费| 亚洲图色成人| 免费看a级黄色片| 久久久色成人| 成年版毛片免费区| 深爱激情五月婷婷| 国内精品久久久久精免费| 精品国产三级普通话版| aaaaa片日本免费| 日韩精品中文字幕看吧| 深夜精品福利| 国产日本99.免费观看| 精品一区二区三区视频在线观看免费| 日本a在线网址| 国产亚洲精品av在线| av视频在线观看入口| h日本视频在线播放| 啪啪无遮挡十八禁网站| 精品一区二区三区视频在线观看免费| 久99久视频精品免费| av天堂在线播放| 欧美中文日本在线观看视频| 亚洲av.av天堂| 国产在线精品亚洲第一网站| 少妇的逼好多水| 成人综合一区亚洲| 欧美性猛交╳xxx乱大交人| 国产探花极品一区二区| 91久久精品国产一区二区成人| 悠悠久久av| 精品无人区乱码1区二区| 日韩欧美精品v在线| 两人在一起打扑克的视频| 日韩中字成人| 一本精品99久久精品77| 欧美黑人欧美精品刺激| 国产成人一区二区在线| 少妇猛男粗大的猛烈进出视频 | 亚洲va在线va天堂va国产| 99久久精品一区二区三区| 亚洲av成人av| 色尼玛亚洲综合影院| 免费看a级黄色片| 国产麻豆成人av免费视频| 国产极品精品免费视频能看的| 最近视频中文字幕2019在线8| 淫妇啪啪啪对白视频| 丰满的人妻完整版| 老司机福利观看| 国产精品福利在线免费观看| 国产精品美女特级片免费视频播放器| 国产精品亚洲一级av第二区| 久久6这里有精品| 国产伦在线观看视频一区| 国产精品三级大全| x7x7x7水蜜桃| 高清在线国产一区| 国产免费一级a男人的天堂| 夜夜爽天天搞| 国产探花极品一区二区| 国产精品自产拍在线观看55亚洲| 日本撒尿小便嘘嘘汇集6| 夜夜看夜夜爽夜夜摸| 免费看美女性在线毛片视频| 免费一级毛片在线播放高清视频| 久久这里只有精品中国| 蜜桃亚洲精品一区二区三区| avwww免费| 极品教师在线视频| 久久国产乱子免费精品| 在线免费观看不下载黄p国产 | 九九热线精品视视频播放| 特级一级黄色大片| 少妇猛男粗大的猛烈进出视频 | 国产精品国产高清国产av| 国产一区二区三区在线臀色熟女| 精品一区二区三区av网在线观看| 波多野结衣高清作品| 午夜免费男女啪啪视频观看 | 最后的刺客免费高清国语| 国产视频一区二区在线看| 一本一本综合久久| 午夜免费成人在线视频| 亚洲av电影不卡..在线观看| 一区福利在线观看| 日韩欧美国产在线观看| 九色成人免费人妻av| 中文字幕av在线有码专区| 国产精品爽爽va在线观看网站| 免费看光身美女| 18禁裸乳无遮挡免费网站照片| 91久久精品电影网| 久久99热6这里只有精品| 3wmmmm亚洲av在线观看| 午夜久久久久精精品| 国产一区二区在线av高清观看| 伦精品一区二区三区| 天天躁日日操中文字幕| 动漫黄色视频在线观看| 99热6这里只有精品| 3wmmmm亚洲av在线观看| 国产精品一区二区免费欧美| 露出奶头的视频| 国产精品国产高清国产av| 久久国内精品自在自线图片| 日本免费一区二区三区高清不卡| 久久久国产成人精品二区| 男女做爰动态图高潮gif福利片| 人妻丰满熟妇av一区二区三区| 黄色视频,在线免费观看| 亚洲av五月六月丁香网| 久久国内精品自在自线图片| 午夜免费成人在线视频| 97超级碰碰碰精品色视频在线观看| 午夜a级毛片| 国产精品电影一区二区三区| 午夜亚洲福利在线播放| 欧美区成人在线视频|