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

    多段翼型高精度數(shù)值模擬技術(shù)研究

    2013-08-21 11:21:32王運(yùn)濤孟德虹鄧小剛
    關(guān)鍵詞:主翼邊界層站位

    王運(yùn)濤,孟德虹,鄧小剛

    (中國(guó)空氣動(dòng)力研究與發(fā)展中心 空氣動(dòng)力學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室,四川 綿陽(yáng) 621000)

    0 引 言

    隨著計(jì)算機(jī)硬件技術(shù)和CFD技術(shù)本身的發(fā)展,采用基于雷諾平均NS方程(RANS)的數(shù)值模擬軟件已經(jīng)可以模擬真實(shí)飛行器的復(fù)雜外形及全機(jī)的復(fù)雜流場(chǎng),包括二維高升力翼型和三維帶增升裝置的全機(jī)構(gòu)型[1]?;赗ANS方程預(yù)測(cè)巡航構(gòu)型氣動(dòng)特性變化趨勢(shì)和模擬全湍流附著流動(dòng)的能力已經(jīng)逐步得到飛行器設(shè)計(jì)工程師的認(rèn)可,但高升力構(gòu)型的數(shù)值模擬可信度水平依然很低。高升力構(gòu)型的主要模擬難點(diǎn)在于:邊界層轉(zhuǎn)捩、激波/邊界層干擾、粘性尾跡區(qū)干擾、尾跡與邊界層摻混和分離流動(dòng)等等。

    為了提高CFD軟件的數(shù)值模擬精度,研究高升力構(gòu)型的流動(dòng)機(jī)理,空氣動(dòng)力學(xué)的實(shí)驗(yàn)工作者和CFD工作者付出了巨大的努力[2-4],高升力構(gòu)型的數(shù)值模擬也是許多CFD可信度專(zhuān)題會(huì)議的主題。Rumsey等人[5]對(duì)多位國(guó)外研究者的工作進(jìn)行了綜述,但這些研究工作基本上全部采用二階空間離散精度的差分格式。三階以上空間離散精度的差分格式受到幾何守恒率、對(duì)接邊界的高精度處理等技術(shù)的影響,在復(fù)雜外形的應(yīng)用上一直處于探索階段,目前高精度格式的應(yīng)用主要集中于簡(jiǎn)單構(gòu)型的復(fù)雜流動(dòng)機(jī)理研究,如氣動(dòng)噪聲機(jī)理、超聲速底部流動(dòng)等。

    鄧小剛團(tuán)隊(duì)從20世紀(jì)90年代致力于高精度計(jì)算方法研究,先后構(gòu)造了3至9階的線性耗散緊致格式(Dissipative Compact Schemes,DCS)和一類(lèi)高階緊致非線性格式(Compact Nonlinear Schemes,CNS),通過(guò)引入加權(quán)思想,又構(gòu)造了一系列加權(quán)緊致非線性格式 WCNS(Weighted Compact Nonlinear Scheme),包括隱式和顯式兩種類(lèi)型,并通過(guò)典型算例驗(yàn)證了WCNS格式的耗散特性和色散特性。近期在幾何守恒率、對(duì)接邊界高精度處理等復(fù)雜外形應(yīng)用技術(shù)方面取得重要進(jìn)展[6-7],為高精度格式的工程應(yīng)用奠定了良好的基礎(chǔ)。

    本文基于空間五階精度的 WCNS格式,綜合利用幾何守恒率和對(duì)接邊界的高精度處理技術(shù),通過(guò)求解任意坐標(biāo)系下的雷諾平均的N-S方程,采用兩種多段翼型的標(biāo)準(zhǔn)算例,考核了WCNS格式的模擬低速?gòu)?fù)雜構(gòu)型的能力。詳細(xì)研究了SA一方程湍流模型、SST兩方程湍流模型、轉(zhuǎn)捩位置對(duì)該翼型壓力分布和典型站位速度型的影響。研究結(jié)果表明,高精度WCNS格式具有模擬低速?gòu)?fù)雜構(gòu)型流場(chǎng)的能力,采用全湍流模擬方式可以較好地模擬多段翼型的壓力分布,但對(duì)邊界層速度型和阻力的模擬精度較差;模擬實(shí)驗(yàn)的轉(zhuǎn)捩位置可以顯著提高附面層內(nèi)速度型的模擬精度和阻力的預(yù)測(cè)精度。

    1 高精度計(jì)算方法

    控制方程為任意坐標(biāo)系的雷諾平均的NS方程:

    其中Q為守恒變量,E、F、G為對(duì)流通量,Ev、Fv、Gv為粘性通量。對(duì)流項(xiàng)離散格式采用5階精度加權(quán)緊致非線性格式(WCNS),粘性項(xiàng)離散格式為4階精度中心格式,對(duì)流項(xiàng)離散及邊界與近邊界的具體格式為:

    (a)內(nèi)點(diǎn)格式

    對(duì)流項(xiàng)離散采用原始變量型的 WCNS-E-5格式,設(shè)網(wǎng)格間距為h,以ξ方向?yàn)槔?/p>

    (b)邊界格式和近邊界格式

    近期的研究發(fā)現(xiàn),幾何守恒率與對(duì)接邊界的處理是制約高精度格式應(yīng)用于復(fù)雜外形的關(guān)鍵因素,這方面的最新研究進(jìn)展可參考文獻(xiàn)[6]和[7],這里不再贅述。

    表1 NLR7301兩段翼型的氣動(dòng)特性(GA P2.6%)Table 1 Aerodynamic character of NLR7301 two element airfoil

    2 NLR7301兩段翼型的高精度數(shù)值模擬

    NLR-7301兩段翼型的風(fēng)洞實(shí)驗(yàn)是在70年代末期在NLR3m×2m低速風(fēng)洞中完成的,實(shí)驗(yàn)結(jié)果包括了總體氣動(dòng)特性、壓力分布、典型站位的速度型分布等多種數(shù)據(jù)[8]。該翼型的襟翼偏角為20°主翼/襟翼重疊區(qū)均為5.3%c;實(shí)驗(yàn)構(gòu)型包含了兩種不同的縫隙寬度,一種為2.6%c,另一種為1.3%c。該翼型是典型的起飛構(gòu)型,主翼的后緣經(jīng)過(guò)仔細(xì)的修型,全流場(chǎng)沒(méi)有明顯的分離。本文采用了縫隙寬度為2.6%c的構(gòu)型,計(jì)算構(gòu)型、典型測(cè)量站位見(jiàn)圖1。

    圖1 NLR7301構(gòu)型及速度型站位Fig.1 NLR7301profile and velocity locations

    本文采用多塊對(duì)接網(wǎng)格技術(shù)模擬該兩段翼型,網(wǎng)格單元數(shù)為20萬(wàn),物面第一層距離7.5E-6,遠(yuǎn)場(chǎng)150c。,局部網(wǎng)格分布見(jiàn)圖2。

    圖2 NLR7301的計(jì)算網(wǎng)格(局部放大)Fig.2 NLR7301computational grid(local)

    ·湍流模型的比較

    采用WCNS格式、Spalart-Allmaras一方程湍流模型和Menter's k-ωSST兩方程模型,全湍流方式數(shù)值模擬了該兩段翼型的繞流流場(chǎng)。來(lái)流馬赫數(shù)0.185,攻角為6°,雷諾數(shù)2.51×106。主翼與襟翼的壓力分布與實(shí)驗(yàn)結(jié)果的比較見(jiàn)圖3,氣動(dòng)特性的數(shù)值模擬結(jié)果見(jiàn)表1。由圖中看出,采用兩種湍流模型得到壓力分布基本一致,且與實(shí)驗(yàn)結(jié)果吻合良好。但從總體氣動(dòng)特性上來(lái)看,采用兩種湍流模型得到的阻力系數(shù)明顯比實(shí)驗(yàn)值偏高,其中,采用SA模型阻力系數(shù)大43個(gè)阻力單位(1個(gè)阻力單位=0.0001),采用SST模型阻力系數(shù)大19個(gè)阻力單位。說(shuō)明采用高精度格式和全湍流方式,可以較好地模擬壓力分布,但阻力系數(shù)明顯偏大。

    圖3 攻角為6°時(shí)壓力分布比較圖Fig.3 Surface pressure coefficient(fully turbulent)

    圖4 給出了采用兩種湍流模型得到的典型站位的速度型與實(shí)驗(yàn)的比較。由圖中看出,兩種湍流模型得到速度型模擬結(jié)果基本一致,但在邊界層厚度、主翼尾跡區(qū)的模擬精度與實(shí)驗(yàn)結(jié)果差別較大,這是導(dǎo)致計(jì)算得到的阻力系數(shù)與實(shí)驗(yàn)結(jié)果差別較大的主要原因。

    圖4 湍流模型對(duì)典型站位速度型的影響Fig.4 Effect of turbulence model to velocity profiles

    ·轉(zhuǎn)捩位置的影響

    NLR-7301的邊界層測(cè)量的實(shí)驗(yàn)結(jié)果給出了典型迎角下,機(jī)翼上翼面、下翼面和襟翼上表面的轉(zhuǎn)捩位置[8],攻角6°的轉(zhuǎn)捩位置見(jiàn)表2。需要說(shuō)明的是,表2中給出的主翼上表面的轉(zhuǎn)捩位置實(shí)際上是主翼前緣的層流分離到湍流再附的起始點(diǎn)和終止點(diǎn),數(shù)值模擬的過(guò)程中將主翼上表面層流分離的起始點(diǎn)作為轉(zhuǎn)捩模擬的起始點(diǎn)。本文采用SA一方程模型,直接模擬實(shí)驗(yàn)的轉(zhuǎn)捩位置,重點(diǎn)研究了轉(zhuǎn)捩位置對(duì)阻力特性、邊界層速度型的影響。

    表2 NLR-7301兩段翼型轉(zhuǎn)捩位置(GA P2.6%)Table 2 Transition location of NLR 7301two element airfoil

    圖5給出了是否模擬轉(zhuǎn)捩位置對(duì)邊界層內(nèi)速度型的影響,包括了主翼上表面、下表面及襟翼上表面的典型站位。與采用全湍流方式相比較,模擬了實(shí)驗(yàn)的轉(zhuǎn)捩位置后,主翼上下表面的邊界層厚度,尤其是主翼下表面的速度型改善尤為明顯,襟翼上表面的邊界層厚度和主翼的尾跡區(qū)的模擬結(jié)果也有明顯的改善。模擬轉(zhuǎn)捩位置后的速度型分布與實(shí)驗(yàn)結(jié)果取得了較好的一致。邊界層內(nèi)速度型模擬精度的提高大大改善了阻力的預(yù)測(cè)精度,與實(shí)驗(yàn)結(jié)果幾乎相同(見(jiàn)表1)。

    圖5 轉(zhuǎn)捩位置對(duì)典型站位速度型的影響Fig.5 Effect of transition location to velocity profiles

    3 30P-30N三段翼型的高精度數(shù)值模擬

    麥道航空公司的30P-30N三段增升構(gòu)型是被CFD工作者廣泛采用的多段構(gòu)型之一,該翼型的前緣縫翼和后緣襟翼的偏角均為30°,前緣縫翼的縫道寬度為2.95%,外伸量為-2.5%;后緣襟翼縫道寬度1.27%,外伸量為0.25%,是典型的著陸構(gòu)型。計(jì)算構(gòu)型、速度型測(cè)量位置如圖6所示。

    本文的數(shù)值模擬采用了多塊對(duì)接結(jié)構(gòu)網(wǎng)格技術(shù),“C”型網(wǎng)格為主導(dǎo),計(jì)算區(qū)域的遠(yuǎn)場(chǎng)邊界取為120倍弦長(zhǎng)。壁面的第一排網(wǎng)格達(dá)到了1.0×10-5弦長(zhǎng),網(wǎng)格在各個(gè)剪切層附近均進(jìn)行了適當(dāng)?shù)募用埽员WC附面層內(nèi)和剪切層的數(shù)值模擬精度,網(wǎng)格規(guī)模達(dá)到了13.8萬(wàn)。圖7給出了計(jì)算網(wǎng)格的局部放大圖。

    圖6 30P-30N構(gòu)型及速度型站位Fig.6 30P-30Nprofile and velocity locations

    圖7 30P-30N的計(jì)算網(wǎng)格(局部放大)Fig.7 30P-30Ncomputational grid(local)

    ·湍流模型的比較

    采用 WCNS格式、Spalart-Allmaras一方程湍流模型和Menter's k-ωSST兩方程模型,全湍流方式數(shù)值模擬了該三段翼型的繞流流場(chǎng)。本文模擬的狀態(tài)為:馬赫數(shù)M=0.2,攻角α=19°,雷諾數(shù)Re=9×106。該攻角已接近失速攻角,這使得該狀態(tài)的流動(dòng)模擬更具挑戰(zhàn)性。計(jì)算得到的壓力分布與實(shí)驗(yàn)結(jié)果的比較見(jiàn)圖8,氣動(dòng)特性的數(shù)值模擬結(jié)果見(jiàn)表3。由圖中看出,采用兩種湍流模型得到壓力分布基本一致,且與實(shí)驗(yàn)結(jié)果吻合良好。但從氣動(dòng)特性上來(lái)看,采用兩種湍流模型得到的阻力系數(shù)相差67個(gè)阻力單位,由于沒(méi)有相應(yīng)的阻力實(shí)驗(yàn)結(jié)果,這里沒(méi)辦法進(jìn)行定量上的比較。

    圖8 攻角為19°時(shí)壓力分布比較圖Fig.8 Surface pressure coefficient(fully turbulent)

    表3 30P-30N三段翼型的氣動(dòng)特性Table 3 Aerodynamic character of 30P-30Nthree element airfoil

    圖9給出了采用兩種湍流模型得到的典型站位的速度型與實(shí)驗(yàn)的比較。由圖中看出,兩種湍流模型得到速度型模擬結(jié)果基本一致,但在縫翼尾跡區(qū)、邊界層厚度、主翼尾跡區(qū)的模擬精度與實(shí)驗(yàn)結(jié)果差別較大。

    圖9 湍流模型對(duì)速度型的影響(全湍流)Fig.9 Effect of turbulence model to velocity profiles

    ·轉(zhuǎn)捩位置的影響

    文獻(xiàn)[5]中給出了相對(duì)于收縮構(gòu)型弦長(zhǎng)c的轉(zhuǎn)捩區(qū)域(見(jiàn)表4)。其中實(shí)驗(yàn)測(cè)量值中“n/a”表示在縫翼下表面型面尖點(diǎn)以后、主翼下表面型面尖點(diǎn)以后和襟翼下表面均沒(méi)有觀察到確定的轉(zhuǎn)捩位置。為了研究轉(zhuǎn)捩位置對(duì)速度型模擬精度的影響,本文選用SA一方程模型數(shù)值模擬了實(shí)驗(yàn)的轉(zhuǎn)捩位置,對(duì)應(yīng)于實(shí)驗(yàn)值沒(méi)有明確轉(zhuǎn)捩位置的點(diǎn),本文分別選擇了縫翼下表面尖點(diǎn)前、主翼下表面尖點(diǎn)前一小段距離為轉(zhuǎn)捩區(qū)域,而在襟翼下表面則為層流區(qū)域。

    表4 計(jì)算與實(shí)驗(yàn)的轉(zhuǎn)捩位置Table 4 Transition location for calculation and test

    圖10給出了采用SA一方程模型是否模擬實(shí)驗(yàn)轉(zhuǎn)捩位置的典型站位速度型的比較??梢钥吹侥M了實(shí)驗(yàn)的轉(zhuǎn)捩位置后,典型站位的速度型得到了明顯的改善。首先是主翼上邊界層的厚度和與實(shí)驗(yàn)結(jié)果更加接近,其次是縫翼和主翼的尾跡區(qū)速度值和尾跡區(qū)的寬度均有明顯改善,但與實(shí)驗(yàn)結(jié)果相比較還存在一定差距。文獻(xiàn)[5]中指出,推遲縫翼上表面的轉(zhuǎn)捩位置可以進(jìn)一步提高縫翼尾跡區(qū)的數(shù)值模擬精度。從表3中可以看出,模擬了實(shí)驗(yàn)的轉(zhuǎn)捩位置后,阻力系數(shù)降低了78個(gè)阻力單位。

    圖10 典型站位速度型的比較(轉(zhuǎn)捩位置A)Fig.10 Velocity profile on typical station(Transition location A)

    4 結(jié) 論

    本文采用五階精度的WCNS格式和結(jié)構(gòu)對(duì)接網(wǎng)格技術(shù),綜合應(yīng)用幾何守恒率和對(duì)接邊界高精度處理技術(shù),通過(guò)求解任意坐標(biāo)系下的RANS方程,數(shù)值模擬了NLR7301二段翼型和30P-30N三段翼型的復(fù)雜流場(chǎng),主要研究了SA和SST兩種湍流模型、轉(zhuǎn)捩位置對(duì)數(shù)值模擬結(jié)果的影響。通過(guò)與相應(yīng)的實(shí)驗(yàn)結(jié)果相比較,得到以下一些基本結(jié)論:(1)采用幾何守恒率和對(duì)接邊界高精度處理技術(shù),WCNS格式可以模擬多段翼型的復(fù)雜流場(chǎng)。(2)采用全湍流模擬方式,兩種湍流模型均可以較好地模擬多段翼型的壓力分布,但速度型的計(jì)算結(jié)果與實(shí)驗(yàn)結(jié)果差別明顯。(3)模擬實(shí)驗(yàn)的轉(zhuǎn)捩位置,對(duì)氣動(dòng)力系數(shù),尤其是阻力系數(shù)影響顯著,計(jì)算結(jié)果更接近實(shí)驗(yàn)值;可以明顯提高邊界層、主翼尾跡區(qū)的數(shù)值模擬精度,但縫翼尾跡區(qū)的計(jì)算精度仍然與實(shí)驗(yàn)結(jié)果存在明顯差距,需要進(jìn)一步研究。

    致謝:本文的研究工作是基于WCNS高精度數(shù)值模擬軟件完成的。在此,對(duì)軟件的開(kāi)發(fā)者表示感謝。

    [1] TINOCO E N,BOGUE D R.Progress toward CFD for full flight envelope[J].Aeronautical Jounal,2005,109:451-460.

    [2] PAUL L JOHNSON,KENNETH M JONES,MICHAEL D MADSON.Experimental investigation of a simplified 3Dhigh lift configuration in support of CFD validation[R].AIAA 2000-4217.

    [3] STUART E ROGERS,KARLIN ROTH,STEVEN M NASH.CFD validation of high-lift flows with significant wind-tunnel effects[R].AIAA 2000-4218.

    [4] 朱自強(qiáng),陳迎春,吳宗成,等.高升力系統(tǒng)外形的數(shù)值模擬計(jì)算[J].航空學(xué)報(bào),2005,26(3):257-262.

    [5] CHRISTOPHER L.RUMSEY,SUSAN X.YING.Prediction of high lift:review of present CFD capability[J].Progress in Aerospace Sciences,2002,38:145-180.

    [6] XIAOGANG DENG,MEILIANG MAO,GUOHUA TU et al.Extending weighted compact nonlinear schemes to complex grids with character-based interface condition[J].AIAA Journal,2010,48(12):2840-2851.

    [7] XIAOGANG DENG,MEILIANG MAO,GUOHUA TU.Geometric conservation law and applications to high-order finite difference schemes with stationary grids[J].Journal of Computational Physics,2010,230(4):1100-1115.

    [8] VAN DEN BERG B.Boundary layer measurements on a two-dimensional wing with flap[R].NLR TR 79009U,1979.

    猜你喜歡
    主翼邊界層站位
    多段翼低雷諾數(shù)繞流渦-邊界層相互干擾
    三黃雞
    某型民機(jī)低速巡航構(gòu)型平尾抖振特性風(fēng)洞試驗(yàn)研究
    提高政治站位 對(duì)標(biāo)國(guó)內(nèi)一流
    建黨百年說(shuō)“站位”
    水上消防(2021年3期)2021-08-21 03:12:00
    基于HIFiRE-2超燃發(fā)動(dòng)機(jī)內(nèi)流道的激波邊界層干擾分析
    提升站位講政治 創(chuàng)新?lián)?dāng)爭(zhēng)出彩
    探究鴨式布局模型飛機(jī)
    航空模型(2017年3期)2017-07-28 21:22:01
    一類(lèi)具有邊界層性質(zhì)的二次奇攝動(dòng)邊值問(wèn)題
    非特征邊界的MHD方程的邊界層
    午夜福利高清视频| 亚洲精品美女久久久久99蜜臀| 一二三四社区在线视频社区8| 黄频高清免费视频| 国内精品久久久久久久电影| 免费观看精品视频网站| 亚洲国产高清在线一区二区三 | 国产精品电影一区二区三区| 国产乱人伦免费视频| 午夜免费成人在线视频| 亚洲无线在线观看| 日本一区二区免费在线视频| 亚洲欧美日韩高清在线视频| 亚洲 欧美 日韩 在线 免费| 男人舔女人下体高潮全视频| 最近最新中文字幕大全免费视频| 成年女人毛片免费观看观看9| 精品一区二区三区av网在线观看| 精品不卡国产一区二区三区| 久久狼人影院| 两个人视频免费观看高清| 两性夫妻黄色片| 性欧美人与动物交配| 一级毛片精品| 久久精品91蜜桃| 亚洲精品久久成人aⅴ小说| 女生性感内裤真人,穿戴方法视频| 99国产精品一区二区蜜桃av| www国产在线视频色| 成人精品一区二区免费| 久久久久久人人人人人| 久久中文看片网| 亚洲精品粉嫩美女一区| 99精品久久久久人妻精品| 高清在线国产一区| 欧美成人午夜精品| 一进一出好大好爽视频| 又大又爽又粗| 级片在线观看| 少妇 在线观看| 亚洲第一av免费看| 99国产综合亚洲精品| 一区在线观看完整版| 亚洲国产中文字幕在线视频| 91在线观看av| 无遮挡黄片免费观看| 中文字幕人妻丝袜一区二区| 午夜精品在线福利| 丰满人妻熟妇乱又伦精品不卡| 免费在线观看视频国产中文字幕亚洲| 涩涩av久久男人的天堂| 国内毛片毛片毛片毛片毛片| 国产精品亚洲av一区麻豆| 一级毛片精品| 亚洲av成人不卡在线观看播放网| 淫秽高清视频在线观看| av福利片在线| 波多野结衣巨乳人妻| 亚洲精品在线美女| 亚洲精品美女久久久久99蜜臀| 18禁裸乳无遮挡免费网站照片 | 国产成人啪精品午夜网站| 欧美日本视频| 欧美av亚洲av综合av国产av| 亚洲欧美一区二区三区黑人| 国产亚洲欧美精品永久| 亚洲色图av天堂| 亚洲国产精品sss在线观看| 国产精品爽爽va在线观看网站 | 久久久精品国产亚洲av高清涩受| 91老司机精品| 纯流量卡能插随身wifi吗| 国产精品免费一区二区三区在线| 99riav亚洲国产免费| 亚洲国产精品久久男人天堂| 日本一区二区免费在线视频| 法律面前人人平等表现在哪些方面| 欧美日本视频| 美女免费视频网站| 啪啪无遮挡十八禁网站| 欧美日韩亚洲国产一区二区在线观看| 欧美激情久久久久久爽电影 | 在线播放国产精品三级| 久久精品91蜜桃| 欧美激情 高清一区二区三区| 亚洲精华国产精华精| 午夜老司机福利片| 亚洲成国产人片在线观看| 老司机福利观看| 欧美丝袜亚洲另类 | 老鸭窝网址在线观看| 久久精品成人免费网站| 免费不卡黄色视频| 最近最新中文字幕大全电影3 | 精品国产一区二区久久| 校园春色视频在线观看| 国产激情欧美一区二区| 亚洲国产精品合色在线| 久久久久久人人人人人| 国产蜜桃级精品一区二区三区| 人人妻人人爽人人添夜夜欢视频| 久久久国产精品麻豆| 国产一区二区三区在线臀色熟女| 青草久久国产| 91麻豆av在线| 久久国产精品人妻蜜桃| 一级毛片高清免费大全| 亚洲电影在线观看av| 国产精品免费一区二区三区在线| 国产熟女午夜一区二区三区| 男男h啪啪无遮挡| 国产精品影院久久| 久久精品国产亚洲av高清一级| 欧美成人免费av一区二区三区| 国产免费男女视频| 欧美午夜高清在线| 午夜老司机福利片| 老熟妇乱子伦视频在线观看| 岛国视频午夜一区免费看| 亚洲色图 男人天堂 中文字幕| 精品久久久久久久久久免费视频| 99国产极品粉嫩在线观看| 香蕉丝袜av| 精品欧美一区二区三区在线| 动漫黄色视频在线观看| 99国产综合亚洲精品| 中文字幕色久视频| 在线观看www视频免费| 免费高清视频大片| 桃色一区二区三区在线观看| 香蕉久久夜色| 欧美性长视频在线观看| 国语自产精品视频在线第100页| 精品国产一区二区三区四区第35| videosex国产| 麻豆国产av国片精品| 中文字幕色久视频| 国产激情欧美一区二区| 精品福利观看| 久热爱精品视频在线9| 人人妻人人澡欧美一区二区 | 一卡2卡三卡四卡精品乱码亚洲| 欧洲精品卡2卡3卡4卡5卡区| 天堂影院成人在线观看| 51午夜福利影视在线观看| 午夜福利欧美成人| 欧美黑人欧美精品刺激| 19禁男女啪啪无遮挡网站| 成人免费观看视频高清| 中国美女看黄片| 国产精品永久免费网站| 久久精品亚洲精品国产色婷小说| 亚洲成人久久性| 1024香蕉在线观看| 精品人妻在线不人妻| 午夜亚洲福利在线播放| 性少妇av在线| 免费久久久久久久精品成人欧美视频| 在线观看日韩欧美| 国产成人精品久久二区二区91| 免费看a级黄色片| 悠悠久久av| 老熟妇乱子伦视频在线观看| 热re99久久国产66热| 午夜成年电影在线免费观看| 首页视频小说图片口味搜索| 久9热在线精品视频| 巨乳人妻的诱惑在线观看| 国产av在哪里看| 亚洲精品中文字幕在线视频| 黄色女人牲交| av超薄肉色丝袜交足视频| 国产麻豆69| 高清在线国产一区| 变态另类成人亚洲欧美熟女 | 国产成人精品在线电影| 两性夫妻黄色片| 精品日产1卡2卡| 中文字幕av电影在线播放| 黄网站色视频无遮挡免费观看| 变态另类丝袜制服| 两人在一起打扑克的视频| 757午夜福利合集在线观看| 久久精品国产清高在天天线| 亚洲免费av在线视频| 视频在线观看一区二区三区| 欧美国产日韩亚洲一区| 亚洲av成人不卡在线观看播放网| 中亚洲国语对白在线视频| 欧美一级a爱片免费观看看 | 欧美日韩亚洲综合一区二区三区_| 男女之事视频高清在线观看| 99国产精品一区二区蜜桃av| www.www免费av| 两人在一起打扑克的视频| 久久影院123| 狂野欧美激情性xxxx| 淫妇啪啪啪对白视频| 91九色精品人成在线观看| 免费久久久久久久精品成人欧美视频| 亚洲国产精品久久男人天堂| 看片在线看免费视频| 久久精品亚洲熟妇少妇任你| 999久久久精品免费观看国产| 国产在线观看jvid| 国产蜜桃级精品一区二区三区| 嫩草影院精品99| av福利片在线| 一级黄色大片毛片| 国产精品精品国产色婷婷| 成年女人毛片免费观看观看9| 亚洲人成网站在线播放欧美日韩| 国产精品av久久久久免费| 法律面前人人平等表现在哪些方面| 亚洲一区二区三区色噜噜| 香蕉丝袜av| 国产aⅴ精品一区二区三区波| 精品电影一区二区在线| 麻豆av在线久日| 一级毛片高清免费大全| 18禁国产床啪视频网站| 老司机午夜福利在线观看视频| 日本黄色视频三级网站网址| 久久人妻熟女aⅴ| 九色国产91popny在线| 97碰自拍视频| 久久午夜亚洲精品久久| 97人妻天天添夜夜摸| 欧美激情极品国产一区二区三区| 午夜福利成人在线免费观看| 亚洲精品在线观看二区| 午夜精品在线福利| 国产精品免费视频内射| 久久中文字幕人妻熟女| 国产激情欧美一区二区| 国产一卡二卡三卡精品| 色av中文字幕| 此物有八面人人有两片| 亚洲成av人片免费观看| 久久久久国产精品人妻aⅴ院| 国产精品二区激情视频| 国产亚洲欧美精品永久| 女性被躁到高潮视频| 女人高潮潮喷娇喘18禁视频| 精品国产一区二区三区四区第35| 国产亚洲精品第一综合不卡| 丰满的人妻完整版| 好男人在线观看高清免费视频 | 亚洲中文日韩欧美视频| 日韩免费av在线播放| 天堂动漫精品| 两人在一起打扑克的视频| 精品卡一卡二卡四卡免费| 久久国产精品男人的天堂亚洲| 久久久久久久精品吃奶| 成人18禁高潮啪啪吃奶动态图| 老熟妇乱子伦视频在线观看| 国产精品久久久人人做人人爽| 高清毛片免费观看视频网站| 啦啦啦观看免费观看视频高清 | 最近最新中文字幕大全免费视频| 可以在线观看的亚洲视频| 真人做人爱边吃奶动态| 午夜影院日韩av| av网站免费在线观看视频| 丰满人妻熟妇乱又伦精品不卡| 99久久99久久久精品蜜桃| 美女大奶头视频| 男女下面插进去视频免费观看| 校园春色视频在线观看| 一边摸一边做爽爽视频免费| 视频区欧美日本亚洲| 欧美黄色淫秽网站| 又黄又爽又免费观看的视频| 人人妻人人爽人人添夜夜欢视频| 三级毛片av免费| 好看av亚洲va欧美ⅴa在| 亚洲自拍偷在线| www.www免费av| 国产激情欧美一区二区| 国产野战对白在线观看| 日本免费a在线| 亚洲男人的天堂狠狠| 黄网站色视频无遮挡免费观看| 女生性感内裤真人,穿戴方法视频| 99久久精品国产亚洲精品| 国产精品一区二区精品视频观看| 精品乱码久久久久久99久播| 十八禁人妻一区二区| 国产99久久九九免费精品| 亚洲欧美日韩另类电影网站| 日本欧美视频一区| 99国产精品一区二区蜜桃av| 亚洲精品国产色婷婷电影| 妹子高潮喷水视频| 久久久水蜜桃国产精品网| 丝袜人妻中文字幕| 国内精品久久久久精免费| 国产精品综合久久久久久久免费 | 99国产精品一区二区蜜桃av| 欧美一级毛片孕妇| 色老头精品视频在线观看| 18禁国产床啪视频网站| 母亲3免费完整高清在线观看| 两个人看的免费小视频| ponron亚洲| 19禁男女啪啪无遮挡网站| av片东京热男人的天堂| 一本大道久久a久久精品| 波多野结衣av一区二区av| 午夜福利欧美成人| 国产亚洲精品久久久久5区| 女警被强在线播放| 夜夜躁狠狠躁天天躁| 亚洲精品一卡2卡三卡4卡5卡| 老熟妇仑乱视频hdxx| 美女高潮喷水抽搐中文字幕| 久久国产精品男人的天堂亚洲| 狠狠狠狠99中文字幕| 老司机福利观看| 午夜激情av网站| 国产高清视频在线播放一区| 国产亚洲精品第一综合不卡| 免费看美女性在线毛片视频| 伦理电影免费视频| 欧美亚洲日本最大视频资源| 亚洲人成77777在线视频| 国产精品 欧美亚洲| 久久国产乱子伦精品免费另类| bbb黄色大片| 日韩精品中文字幕看吧| 91麻豆精品激情在线观看国产| 久久久国产成人免费| 脱女人内裤的视频| 国产真人三级小视频在线观看| 亚洲黑人精品在线| 国产1区2区3区精品| 日韩成人在线观看一区二区三区| 1024视频免费在线观看| 50天的宝宝边吃奶边哭怎么回事| 亚洲欧美日韩无卡精品| 久久午夜亚洲精品久久| 免费在线观看视频国产中文字幕亚洲| 99re在线观看精品视频| 99国产综合亚洲精品| 真人一进一出gif抽搐免费| 99久久精品国产亚洲精品| 中文字幕色久视频| 丝袜美足系列| 欧美精品亚洲一区二区| 此物有八面人人有两片| 久久久国产欧美日韩av| 深夜精品福利| 色综合站精品国产| 亚洲精华国产精华精| 日韩av在线大香蕉| 欧美色视频一区免费| 久99久视频精品免费| 高潮久久久久久久久久久不卡| 大型黄色视频在线免费观看| 久久香蕉国产精品| www国产在线视频色| 在线观看免费视频日本深夜| 黄色视频不卡| 中亚洲国语对白在线视频| 免费久久久久久久精品成人欧美视频| 午夜激情av网站| 国产一区二区激情短视频| 成人18禁高潮啪啪吃奶动态图| 久久婷婷人人爽人人干人人爱 | 日韩精品免费视频一区二区三区| 国产精品 欧美亚洲| 日本 av在线| 午夜成年电影在线免费观看| 欧美午夜高清在线| 精品一区二区三区av网在线观看| 亚洲精品美女久久av网站| 怎么达到女性高潮| 黄色视频不卡| 色播在线永久视频| 免费观看人在逋| 一本大道久久a久久精品| 欧美激情高清一区二区三区| 色婷婷久久久亚洲欧美| a级毛片在线看网站| 国产精品免费视频内射| 禁无遮挡网站| ponron亚洲| 可以在线观看的亚洲视频| 亚洲全国av大片| 国产欧美日韩综合在线一区二区| av天堂在线播放| 亚洲成a人片在线一区二区| 国产精品影院久久| 中文字幕色久视频| 日本 av在线| 日本三级黄在线观看| 久热这里只有精品99| www.熟女人妻精品国产| 88av欧美| 一级毛片高清免费大全| 99国产综合亚洲精品| 又黄又粗又硬又大视频| 啪啪无遮挡十八禁网站| 一本综合久久免费| 久久久水蜜桃国产精品网| 亚洲精品国产一区二区精华液| 女同久久另类99精品国产91| 视频区欧美日本亚洲| 成人18禁在线播放| av天堂在线播放| 这个男人来自地球电影免费观看| 91精品国产国语对白视频| 国产成人系列免费观看| 禁无遮挡网站| 老司机深夜福利视频在线观看| 色综合亚洲欧美另类图片| 视频在线观看一区二区三区| 在线观看舔阴道视频| 欧美av亚洲av综合av国产av| 国产欧美日韩一区二区精品| 欧美绝顶高潮抽搐喷水| 亚洲男人天堂网一区| 女人精品久久久久毛片| 激情视频va一区二区三区| videosex国产| 欧美国产日韩亚洲一区| 国产成人欧美| 男女下面插进去视频免费观看| 欧美乱码精品一区二区三区| 色综合站精品国产| 国产精品永久免费网站| 美女午夜性视频免费| 18禁黄网站禁片午夜丰满| 在线观看免费视频日本深夜| 黄片小视频在线播放| 免费在线观看黄色视频的| 国产日韩一区二区三区精品不卡| 久热爱精品视频在线9| 岛国视频午夜一区免费看| 欧美人与性动交α欧美精品济南到| 国内精品久久久久精免费| 国产精品综合久久久久久久免费 | 久久精品影院6| 久久这里只有精品19| 咕卡用的链子| 国产又爽黄色视频| 啪啪无遮挡十八禁网站| 制服人妻中文乱码| 啦啦啦免费观看视频1| 非洲黑人性xxxx精品又粗又长| 国产99白浆流出| 亚洲五月天丁香| 成人永久免费在线观看视频| 欧美黑人欧美精品刺激| 国产av一区在线观看免费| 国产不卡一卡二| 欧美黄色淫秽网站| 国产成人精品在线电影| 久久青草综合色| 99热只有精品国产| 国产精品一区二区在线不卡| 国产精品久久久久久精品电影 | 在线播放国产精品三级| av在线播放免费不卡| 久久精品国产99精品国产亚洲性色 | 午夜福利高清视频| 日本 欧美在线| 超碰成人久久| 亚洲熟妇熟女久久| 久久国产精品人妻蜜桃| 免费在线观看完整版高清| 国产aⅴ精品一区二区三区波| 婷婷六月久久综合丁香| 国产精品亚洲av一区麻豆| 波多野结衣高清无吗| 欧美色视频一区免费| 99久久久亚洲精品蜜臀av| 两性夫妻黄色片| 悠悠久久av| 成人手机av| 欧美日韩黄片免| 久久午夜亚洲精品久久| 亚洲中文字幕一区二区三区有码在线看 | 一a级毛片在线观看| 如日韩欧美国产精品一区二区三区| 一级a爱视频在线免费观看| 欧美成人性av电影在线观看| 精品卡一卡二卡四卡免费| 变态另类丝袜制服| 我的亚洲天堂| 精品乱码久久久久久99久播| 国产av一区在线观看免费| 午夜a级毛片| 婷婷丁香在线五月| 色在线成人网| 桃色一区二区三区在线观看| 婷婷精品国产亚洲av在线| 午夜老司机福利片| 午夜福利高清视频| 久久香蕉激情| 色哟哟哟哟哟哟| 给我免费播放毛片高清在线观看| 亚洲成人免费电影在线观看| 18禁裸乳无遮挡免费网站照片 | 国产成人免费无遮挡视频| 国产午夜精品久久久久久| 国产99久久九九免费精品| 亚洲成av片中文字幕在线观看| xxx96com| 欧美黑人欧美精品刺激| 久久久久精品国产欧美久久久| 欧美在线黄色| 可以免费在线观看a视频的电影网站| 此物有八面人人有两片| 非洲黑人性xxxx精品又粗又长| 18禁观看日本| 黄色片一级片一级黄色片| 不卡一级毛片| 麻豆国产av国片精品| 欧美亚洲日本最大视频资源| 香蕉国产在线看| 亚洲一区二区三区色噜噜| 久久午夜亚洲精品久久| 亚洲精品中文字幕在线视频| 如日韩欧美国产精品一区二区三区| 动漫黄色视频在线观看| 91国产中文字幕| 夜夜看夜夜爽夜夜摸| 久久精品成人免费网站| 麻豆成人av在线观看| 日韩高清综合在线| 亚洲熟女毛片儿| 久久久久九九精品影院| 亚洲电影在线观看av| 国产欧美日韩一区二区精品| 欧美一区二区精品小视频在线| 久99久视频精品免费| 动漫黄色视频在线观看| 桃色一区二区三区在线观看| 宅男免费午夜| 久久精品亚洲精品国产色婷小说| 国产成人影院久久av| 人人妻,人人澡人人爽秒播| 国产精品永久免费网站| or卡值多少钱| 在线观看免费视频日本深夜| 欧美色欧美亚洲另类二区 | 日日摸夜夜添夜夜添小说| 高清黄色对白视频在线免费看| 黄网站色视频无遮挡免费观看| 夜夜夜夜夜久久久久| 1024香蕉在线观看| av电影中文网址| 嫁个100分男人电影在线观看| 亚洲精品一区av在线观看| av天堂久久9| 视频区欧美日本亚洲| 国语自产精品视频在线第100页| 国产极品粉嫩免费观看在线| 精品人妻在线不人妻| 成人精品一区二区免费| 久久天躁狠狠躁夜夜2o2o| 午夜两性在线视频| 欧美成狂野欧美在线观看| 好男人电影高清在线观看| 男人操女人黄网站| 色哟哟哟哟哟哟| 国产熟女午夜一区二区三区| 亚洲国产欧美日韩在线播放| 久久久久久亚洲精品国产蜜桃av| 日韩精品免费视频一区二区三区| 精品第一国产精品| 久久香蕉国产精品| 欧美久久黑人一区二区| 亚洲国产看品久久| 人妻久久中文字幕网| 国产欧美日韩一区二区三| 久久人妻av系列| 国产一区二区三区综合在线观看| 波多野结衣巨乳人妻| 黄频高清免费视频| 欧美成人午夜精品| 可以在线观看的亚洲视频| 亚洲av成人不卡在线观看播放网| 亚洲av电影不卡..在线观看| 88av欧美| 久久精品国产综合久久久| 人成视频在线观看免费观看| 亚洲人成网站在线播放欧美日韩| 亚洲一卡2卡3卡4卡5卡精品中文| 精品人妻1区二区| 亚洲avbb在线观看| 十八禁网站免费在线| 亚洲专区字幕在线| 亚洲人成电影观看| 午夜精品国产一区二区电影| 国产单亲对白刺激| 免费av毛片视频| 看片在线看免费视频| 久久久久久久精品吃奶| 国产亚洲精品久久久久久毛片| 欧美日韩一级在线毛片| 免费看a级黄色片| 天堂影院成人在线观看| 亚洲人成77777在线视频| 精品福利观看| 国产欧美日韩一区二区精品| 亚洲精品在线观看二区| 熟女少妇亚洲综合色aaa.| 免费在线观看完整版高清| 国产亚洲精品av在线| 久久天躁狠狠躁夜夜2o2o| 亚洲全国av大片| 亚洲avbb在线观看|