• <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
    给我免费播放毛片高清在线观看| 亚洲欧美日韩无卡精品| 国产一区亚洲一区在线观看| 亚洲av第一区精品v没综合| 午夜福利视频1000在线观看| 欧美高清性xxxxhd video| 国产伦精品一区二区三区四那| 精品不卡国产一区二区三区| 欧美bdsm另类| 亚洲aⅴ乱码一区二区在线播放| 一进一出抽搐动态| 又黄又爽又刺激的免费视频.| 中文资源天堂在线| 午夜福利高清视频| 日韩三级伦理在线观看| 日韩亚洲欧美综合| 99久久人妻综合| 亚洲无线观看免费| 精品国内亚洲2022精品成人| 男人舔女人下体高潮全视频| 日本黄色视频三级网站网址| 女同久久另类99精品国产91| 国内精品一区二区在线观看| 成人高潮视频无遮挡免费网站| 婷婷亚洲欧美| 国产伦在线观看视频一区| 一个人看视频在线观看www免费| 岛国毛片在线播放| 亚洲人成网站在线播放欧美日韩| 在线观看66精品国产| 在线天堂最新版资源| 色哟哟哟哟哟哟| 久久99精品国语久久久| 久久精品夜夜夜夜夜久久蜜豆| 国产精品久久久久久精品电影| 国产成年人精品一区二区| 亚洲va在线va天堂va国产| 亚洲精品乱码久久久v下载方式| 日韩国内少妇激情av| 亚洲精品日韩在线中文字幕 | 日韩高清综合在线| 久久久久久久久久久丰满| 国产毛片a区久久久久| 亚洲欧美精品专区久久| 搡老妇女老女人老熟妇| 99热这里只有精品一区| 亚洲成人av在线免费| 在线观看一区二区三区| 午夜老司机福利剧场| 色尼玛亚洲综合影院| 99久久九九国产精品国产免费| 国产精品综合久久久久久久免费| 中文亚洲av片在线观看爽| 99国产极品粉嫩在线观看| 日日撸夜夜添| 国产精品久久久久久精品电影| 禁无遮挡网站| 日本爱情动作片www.在线观看| 超碰av人人做人人爽久久| 狂野欧美激情性xxxx在线观看| 一个人免费在线观看电影| 欧美在线一区亚洲| 久久午夜亚洲精品久久| 久久精品国产亚洲网站| 老司机福利观看| 观看免费一级毛片| 97超碰精品成人国产| 国产精品国产高清国产av| 人妻夜夜爽99麻豆av| 此物有八面人人有两片| 国产不卡一卡二| 亚洲精品粉嫩美女一区| 亚洲欧美日韩卡通动漫| 亚洲四区av| 高清在线视频一区二区三区 | 欧美变态另类bdsm刘玥| 国产精品福利在线免费观看| 天堂√8在线中文| 一夜夜www| 国产精品99久久久久久久久| 久久精品夜色国产| 99riav亚洲国产免费| 日本欧美国产在线视频| 非洲黑人性xxxx精品又粗又长| 国产成人aa在线观看| 青春草国产在线视频 | 免费黄网站久久成人精品| 日韩视频在线欧美| 欧美日韩综合久久久久久| 夜夜看夜夜爽夜夜摸| 精品一区二区三区人妻视频| 2022亚洲国产成人精品| av国产免费在线观看| 国产伦精品一区二区三区视频9| 爱豆传媒免费全集在线观看| 免费大片18禁| 精品一区二区三区视频在线| 久久人人爽人人片av| 老师上课跳d突然被开到最大视频| 高清毛片免费看| 国产日本99.免费观看| 一本久久中文字幕| 国产精品久久久久久久久免| 国产伦精品一区二区三区四那| av在线亚洲专区| 日本黄大片高清| 热99在线观看视频| 青春草亚洲视频在线观看| 国产亚洲av嫩草精品影院| 久久久久久久久久久丰满| 看黄色毛片网站| 尤物成人国产欧美一区二区三区| 精品99又大又爽又粗少妇毛片| 精品国内亚洲2022精品成人| 欧美最新免费一区二区三区| 日日啪夜夜撸| 中文亚洲av片在线观看爽| 亚洲精品日韩av片在线观看| 国产精品女同一区二区软件| 日本免费a在线| 亚洲国产精品成人综合色| 亚洲无线在线观看| 亚洲精品粉嫩美女一区| 色综合站精品国产| 夜夜爽天天搞| 久久久国产成人免费| 女的被弄到高潮叫床怎么办| 99精品在免费线老司机午夜| 白带黄色成豆腐渣| 欧美日韩综合久久久久久| 亚洲欧美中文字幕日韩二区| 六月丁香七月| 久久午夜福利片| 少妇人妻精品综合一区二区 | 国产探花在线观看一区二区| 国产午夜精品一二区理论片| 天天躁夜夜躁狠狠久久av| 深爱激情五月婷婷| 国内久久婷婷六月综合欲色啪| 亚洲第一区二区三区不卡| 欧美精品国产亚洲| 欧美成人免费av一区二区三区| 婷婷精品国产亚洲av| 黑人高潮一二区| 国产美女午夜福利| 国产黄片美女视频| 国产一区二区在线av高清观看| 在线免费十八禁| 国产午夜精品一二区理论片| 久久99热这里只有精品18| 午夜精品在线福利| 久久久久久久久久成人| 亚洲一级一片aⅴ在线观看| 国产人妻一区二区三区在| 国产午夜精品论理片| 亚洲美女视频黄频| 美女 人体艺术 gogo| 麻豆成人av视频| www日本黄色视频网| 亚洲高清免费不卡视频| 日韩欧美在线乱码| 草草在线视频免费看| 国产精华一区二区三区| 日韩成人伦理影院| 我的女老师完整版在线观看| 国产片特级美女逼逼视频| 亚洲国产精品国产精品| 热99re8久久精品国产| 欧美成人一区二区免费高清观看| videossex国产| 国产黄色视频一区二区在线观看 | 久久人人爽人人片av| 精品国内亚洲2022精品成人| 波野结衣二区三区在线| 亚洲自拍偷在线| 91在线精品国自产拍蜜月| 尤物成人国产欧美一区二区三区| 国产极品天堂在线| 欧美成人一区二区免费高清观看| 97超碰精品成人国产| 亚洲欧洲国产日韩| 热99re8久久精品国产| 少妇熟女aⅴ在线视频| 美女国产视频在线观看| 青春草国产在线视频 | 一级二级三级毛片免费看| 国产久久久一区二区三区| 国内精品美女久久久久久| 最近视频中文字幕2019在线8| 成熟少妇高潮喷水视频| 国内少妇人妻偷人精品xxx网站| 日韩成人av中文字幕在线观看| 色尼玛亚洲综合影院| 国产爱豆传媒在线观看| 性插视频无遮挡在线免费观看| 亚洲在线自拍视频| 婷婷色综合大香蕉| 简卡轻食公司| 亚洲国产精品合色在线| 一级毛片电影观看 | 国产精品麻豆人妻色哟哟久久 | 网址你懂的国产日韩在线| 欧美精品国产亚洲| 久久精品国产亚洲av天美| 日本黄大片高清| 亚洲av一区综合| 国产又黄又爽又无遮挡在线| 欧美bdsm另类| 欧美色欧美亚洲另类二区| 人人妻人人澡人人爽人人夜夜 | 精品一区二区三区人妻视频| 亚洲精品影视一区二区三区av| 国产av麻豆久久久久久久| 亚洲成人精品中文字幕电影| 国产亚洲精品久久久com| 亚洲性久久影院| 久久中文看片网| 国产视频内射| 成人毛片a级毛片在线播放| 日韩成人av中文字幕在线观看| 国产高清三级在线| 亚洲成人中文字幕在线播放| 国产黄色视频一区二区在线观看 | 亚洲国产高清在线一区二区三| 啦啦啦观看免费观看视频高清| 午夜精品一区二区三区免费看| 国产亚洲精品av在线| 美女脱内裤让男人舔精品视频 | 国产一区亚洲一区在线观看| 国产精华一区二区三区| 十八禁国产超污无遮挡网站| 我的老师免费观看完整版| 国产精品三级大全| 最近最新中文字幕大全电影3| 亚洲精品久久久久久婷婷小说 | 成年免费大片在线观看| 一本精品99久久精品77| 国产高清不卡午夜福利| 深夜a级毛片| 欧美性感艳星| 亚洲人与动物交配视频| 久久久久久久久中文| 中文字幕av成人在线电影| 国产精品,欧美在线| 综合色丁香网| 日韩一区二区三区影片| 日韩成人伦理影院| 村上凉子中文字幕在线| 日本免费a在线| 午夜福利成人在线免费观看| 看十八女毛片水多多多| 久久精品国产清高在天天线| 亚洲va在线va天堂va国产| 久久久精品欧美日韩精品| 欧美色视频一区免费| 男女那种视频在线观看| 国产精品.久久久| 色哟哟哟哟哟哟| 国产伦精品一区二区三区四那| .国产精品久久| 成人国产麻豆网| 久久久久久大精品| 午夜免费激情av| 亚洲国产精品国产精品| 中文亚洲av片在线观看爽| 亚洲精品色激情综合| 99久国产av精品国产电影| 99久久成人亚洲精品观看| av在线观看视频网站免费| 春色校园在线视频观看| 亚洲欧美日韩东京热| 亚洲无线在线观看| 99久久精品一区二区三区| 婷婷六月久久综合丁香| 免费看日本二区| 国内少妇人妻偷人精品xxx网站| 成人无遮挡网站| 久久久久免费精品人妻一区二区| 国国产精品蜜臀av免费| 国产精品电影一区二区三区| www.色视频.com| 女人被狂操c到高潮| 91久久精品电影网| 99久久精品国产国产毛片| 久久鲁丝午夜福利片| 国产一区二区三区在线臀色熟女| 可以在线观看毛片的网站| 久久久成人免费电影| 嘟嘟电影网在线观看| 国产乱人视频| 久久久国产成人免费| 黄色日韩在线| 久久韩国三级中文字幕| 国产精品乱码一区二三区的特点| 久久久精品欧美日韩精品| 99热精品在线国产| 18禁黄网站禁片免费观看直播| 99久久无色码亚洲精品果冻| 亚洲成人精品中文字幕电影| 亚洲成人久久性| 亚洲熟妇中文字幕五十中出| 少妇人妻精品综合一区二区 | 你懂的网址亚洲精品在线观看 | av免费在线看不卡| 精品人妻熟女av久视频| 欧美激情国产日韩精品一区| 日日撸夜夜添| 秋霞在线观看毛片| 欧美日韩综合久久久久久| 久久久精品94久久精品| 亚洲精品乱码久久久v下载方式| 久久久久久九九精品二区国产| 欧美zozozo另类| 国产 一区 欧美 日韩| 国产免费一级a男人的天堂| a级毛片a级免费在线| 免费在线观看成人毛片| 亚洲av二区三区四区| 欧美一区二区国产精品久久精品| 一卡2卡三卡四卡精品乱码亚洲| 成年女人看的毛片在线观看| 男女边吃奶边做爰视频| 极品教师在线视频| 婷婷亚洲欧美| 精品久久久久久久久av| 午夜激情欧美在线| 国产成人freesex在线| 久久精品国产自在天天线| 卡戴珊不雅视频在线播放| 国产91av在线免费观看| 熟妇人妻久久中文字幕3abv| 久久这里只有精品中国| 1024手机看黄色片| 欧美日韩国产亚洲二区| 久久这里只有精品中国| 国产精品嫩草影院av在线观看| 国语自产精品视频在线第100页| 精品免费久久久久久久清纯| 精品少妇黑人巨大在线播放 | 国产亚洲精品久久久com| 三级男女做爰猛烈吃奶摸视频| 69人妻影院| 一个人看视频在线观看www免费| 青春草亚洲视频在线观看| 欧美变态另类bdsm刘玥| 日韩强制内射视频| 免费av观看视频| 色哟哟·www| 久久久精品大字幕| 国产黄a三级三级三级人| 黄色欧美视频在线观看| 九色成人免费人妻av| 亚洲av成人精品一区久久| av免费在线看不卡| 久久99热6这里只有精品| 国产黄片美女视频| 中文字幕久久专区| 一进一出抽搐gif免费好疼| 变态另类丝袜制服| 国产麻豆成人av免费视频| 亚洲欧美日韩东京热| 亚洲熟妇中文字幕五十中出| 啦啦啦啦在线视频资源| 全区人妻精品视频| av在线天堂中文字幕| 久久久久久久久中文| 久久久久久久午夜电影| 久久久欧美国产精品| 床上黄色一级片| 欧美潮喷喷水| 亚洲国产欧洲综合997久久,| 免费观看精品视频网站| 中文精品一卡2卡3卡4更新| 国产av一区在线观看免费| 国产精品av视频在线免费观看| 欧美成人免费av一区二区三区| 亚洲成人中文字幕在线播放| 精品久久久久久成人av| 国产一区二区在线观看日韩| 亚洲精品乱码久久久v下载方式| 日韩精品有码人妻一区| 天堂影院成人在线观看| 免费电影在线观看免费观看| 精品久久久久久久久久免费视频| 深夜精品福利| av专区在线播放| 国产黄a三级三级三级人| 美女被艹到高潮喷水动态| 亚洲国产欧美人成| 国产精品不卡视频一区二区| 一级毛片久久久久久久久女| 少妇熟女aⅴ在线视频| 女人十人毛片免费观看3o分钟| 久久精品国产亚洲av涩爱 | 日日摸夜夜添夜夜添av毛片| 校园人妻丝袜中文字幕| 国产白丝娇喘喷水9色精品| 男女视频在线观看网站免费| 在线天堂最新版资源| 久久九九热精品免费| 亚洲成人久久性| 亚洲五月天丁香| 国产在线男女| 国产精品麻豆人妻色哟哟久久 | 哪里可以看免费的av片| 久久精品夜色国产| 免费无遮挡裸体视频| 亚洲熟妇中文字幕五十中出| 99热这里只有精品一区| 少妇人妻一区二区三区视频| 国产91av在线免费观看| 国产精品国产三级国产av玫瑰| 精品久久久噜噜| 有码 亚洲区| 国产在线男女| 国产黄色小视频在线观看| 亚洲欧美日韩东京热| 一区二区三区四区激情视频 | 一本一本综合久久| 亚洲精品亚洲一区二区| 国产69精品久久久久777片| 亚洲成人中文字幕在线播放| 国产精品久久电影中文字幕| 免费不卡的大黄色大毛片视频在线观看 | 久久久久久久久久黄片| 国产精品电影一区二区三区| 99九九线精品视频在线观看视频| 国产精品一区二区三区四区免费观看| 嫩草影院入口| 国产成人精品久久久久久| 精品人妻熟女av久视频| 亚洲欧美精品专区久久| 亚洲久久久久久中文字幕| 美女国产视频在线观看| 国产精品电影一区二区三区| 18禁在线无遮挡免费观看视频| 国产白丝娇喘喷水9色精品| 97热精品久久久久久| 亚洲国产高清在线一区二区三| 久久久色成人| 久久99精品国语久久久| 欧美成人精品欧美一级黄| 日韩精品有码人妻一区| 国产高清三级在线| 一级毛片电影观看 | 波野结衣二区三区在线| 日韩亚洲欧美综合| 亚洲av中文av极速乱| 性欧美人与动物交配| 欧美bdsm另类| 人妻久久中文字幕网| 日本五十路高清| 亚洲av熟女| 亚洲欧美日韩无卡精品| 春色校园在线视频观看| 少妇高潮的动态图| 国产精品日韩av在线免费观看| 午夜爱爱视频在线播放| 中文亚洲av片在线观看爽| 69人妻影院| 免费观看精品视频网站| 亚洲av.av天堂| 成人亚洲精品av一区二区| 久久99蜜桃精品久久| 国产精品美女特级片免费视频播放器| 网址你懂的国产日韩在线| 一级二级三级毛片免费看| 看黄色毛片网站| 一边亲一边摸免费视频| 人人妻人人看人人澡| 观看美女的网站| av天堂中文字幕网| 99久久精品国产国产毛片| 大型黄色视频在线免费观看| 夜夜夜夜夜久久久久| 日韩精品有码人妻一区| 午夜精品一区二区三区免费看| 淫秽高清视频在线观看| 岛国毛片在线播放| 亚洲欧美日韩高清专用| 国产精品福利在线免费观看| 国内久久婷婷六月综合欲色啪| 中文精品一卡2卡3卡4更新| 18+在线观看网站| av视频在线观看入口| 中文亚洲av片在线观看爽| 日韩 亚洲 欧美在线| 国产精品无大码| 久久久久久久久久久丰满| 久久人人爽人人爽人人片va| 中国美白少妇内射xxxbb| 免费搜索国产男女视频| 亚洲国产高清在线一区二区三| 国产麻豆成人av免费视频| 五月伊人婷婷丁香| 观看免费一级毛片| 久久久精品欧美日韩精品| 欧美不卡视频在线免费观看| 亚洲av免费高清在线观看| 日产精品乱码卡一卡2卡三| 精品久久久久久成人av| 一边亲一边摸免费视频| 熟妇人妻久久中文字幕3abv| 国产成人影院久久av| 波多野结衣高清无吗| 老司机福利观看| 亚洲18禁久久av| 亚洲无线观看免费| 能在线免费看毛片的网站| 此物有八面人人有两片| 亚洲精品乱码久久久v下载方式| 又粗又硬又长又爽又黄的视频 | 91aial.com中文字幕在线观看| 亚洲精品乱码久久久v下载方式| 长腿黑丝高跟| 久久精品国产亚洲av涩爱 | 日韩一本色道免费dvd| 晚上一个人看的免费电影| 色吧在线观看| 中文字幕av在线有码专区| av在线天堂中文字幕| 波野结衣二区三区在线| 日本免费a在线| 欧美区成人在线视频| 日韩 亚洲 欧美在线| 日韩精品有码人妻一区| 精品国内亚洲2022精品成人| 精品免费久久久久久久清纯| 精品一区二区三区人妻视频| 亚洲精品456在线播放app| 97超视频在线观看视频| 国产伦理片在线播放av一区 | 国产欧美日韩精品一区二区| 国产av不卡久久| 可以在线观看毛片的网站| 少妇人妻一区二区三区视频| 波多野结衣高清作品| 别揉我奶头 嗯啊视频| 少妇丰满av| 国产精品伦人一区二区| 国产三级在线视频| 特大巨黑吊av在线直播| 欧美精品国产亚洲| 成人高潮视频无遮挡免费网站| 国产亚洲av嫩草精品影院| 久久99蜜桃精品久久| 国产又黄又爽又无遮挡在线| 深爱激情五月婷婷| 亚洲电影在线观看av| 日本黄色片子视频| 久久久久久久久久黄片| 伊人久久精品亚洲午夜| 1024手机看黄色片| 日韩中字成人| 日韩一区二区视频免费看| 一区二区三区四区激情视频 | 97热精品久久久久久| 日日啪夜夜撸| 国产亚洲av片在线观看秒播厂 | 精品久久久噜噜| av在线天堂中文字幕| 国产精品一区二区在线观看99 | 一进一出抽搐gif免费好疼| 欧美极品一区二区三区四区| 床上黄色一级片| 白带黄色成豆腐渣| 最近2019中文字幕mv第一页| 国产精品久久久久久精品电影| kizo精华| 日本一二三区视频观看| 国产精品无大码| 精品久久久噜噜| 最近中文字幕高清免费大全6| 我的老师免费观看完整版| 特大巨黑吊av在线直播| 色哟哟·www| 变态另类成人亚洲欧美熟女| 国产视频内射| 国产一区亚洲一区在线观看| 看免费成人av毛片| 人妻少妇偷人精品九色| 一区福利在线观看| 国产亚洲精品av在线| 久久久久九九精品影院| 国产 一区精品| 国产精品野战在线观看| 欧美成人一区二区免费高清观看| 国产三级在线视频| 国产日韩欧美在线精品| 嫩草影院新地址| 欧美精品一区二区大全| 亚洲人成网站在线播| 丰满乱子伦码专区| 亚洲人成网站在线播放欧美日韩| 成人二区视频| 亚洲精品粉嫩美女一区| 日本成人三级电影网站| 国产69精品久久久久777片| 中文精品一卡2卡3卡4更新| 国产毛片a区久久久久| 午夜激情欧美在线| 国产精品国产三级国产av玫瑰| 日本熟妇午夜| 最近2019中文字幕mv第一页| 免费av毛片视频| 两性午夜刺激爽爽歪歪视频在线观看| 成人性生交大片免费视频hd| www.色视频.com| 熟妇人妻久久中文字幕3abv| 亚洲欧美中文字幕日韩二区| avwww免费| av在线亚洲专区| av女优亚洲男人天堂| 菩萨蛮人人尽说江南好唐韦庄 | 欧美不卡视频在线免费观看|