• <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方程的邊界層
    免费看不卡的av| 女人高潮潮喷娇喘18禁视频| 美女福利国产在线| 亚洲欧洲国产日韩| 大陆偷拍与自拍| 伊人久久国产一区二区| 嫩草影院入口| 免费不卡的大黄色大毛片视频在线观看| 免费观看在线日韩| 美女午夜性视频免费| av女优亚洲男人天堂| 91久久精品国产一区二区三区| 91国产中文字幕| 日本黄色日本黄色录像| 亚洲精品一区蜜桃| 一级毛片我不卡| 国产精品免费大片| 欧美日韩成人在线一区二区| 久久精品国产亚洲av涩爱| 男人舔女人的私密视频| 人妻人人澡人人爽人人| 多毛熟女@视频| 香蕉国产在线看| 国产精品一区二区在线观看99| 午夜福利视频在线观看免费| 午夜激情av网站| 国产成人av激情在线播放| 国产国语露脸激情在线看| 欧美日韩一区二区视频在线观看视频在线| 日韩欧美一区视频在线观看| 少妇熟女欧美另类| 中文字幕av电影在线播放| 99久久中文字幕三级久久日本| av网站免费在线观看视频| 91精品三级在线观看| 国产乱人偷精品视频| 蜜桃国产av成人99| 看免费av毛片| 亚洲成国产人片在线观看| 男女边摸边吃奶| 国产在线免费精品| 天天影视国产精品| 丰满饥渴人妻一区二区三| 99久久综合免费| 国产成人精品久久二区二区91 | 国产精品秋霞免费鲁丝片| 丰满迷人的少妇在线观看| 午夜福利视频精品| 曰老女人黄片| 超碰97精品在线观看| 国产午夜精品一二区理论片| 免费高清在线观看日韩| av网站免费在线观看视频| 亚洲av国产av综合av卡| 久久久久精品久久久久真实原创| 青春草亚洲视频在线观看| 高清欧美精品videossex| 亚洲欧美精品综合一区二区三区 | 人成视频在线观看免费观看| 日产精品乱码卡一卡2卡三| 99热国产这里只有精品6| 色哟哟·www| 国产片内射在线| 国产1区2区3区精品| videossex国产| 丝袜脚勾引网站| 国产一级毛片在线| 91成人精品电影| 80岁老熟妇乱子伦牲交| 国产一区二区三区综合在线观看| 中文字幕人妻丝袜制服| 国产极品天堂在线| 国产成人91sexporn| 国产片内射在线| 美国免费a级毛片| 精品少妇一区二区三区视频日本电影 | 天天操日日干夜夜撸| 国产熟女午夜一区二区三区| 色播在线永久视频| 亚洲伊人色综图| 欧美日韩综合久久久久久| 日韩一卡2卡3卡4卡2021年| 日韩欧美精品免费久久| 秋霞伦理黄片| 亚洲熟女精品中文字幕| 一级毛片我不卡| 国产成人精品一,二区| 国产成人精品久久二区二区91 | 啦啦啦中文免费视频观看日本| 久久久久久人人人人人| 国产熟女欧美一区二区| 国产精品 欧美亚洲| 秋霞在线观看毛片| 久久精品人人爽人人爽视色| 自拍欧美九色日韩亚洲蝌蚪91| av不卡在线播放| 国产xxxxx性猛交| 老鸭窝网址在线观看| 男女高潮啪啪啪动态图| 建设人人有责人人尽责人人享有的| 欧美黄色片欧美黄色片| 亚洲av成人精品一二三区| 午夜日韩欧美国产| 男女无遮挡免费网站观看| 日韩一本色道免费dvd| 国产高清国产精品国产三级| 久久久久久人人人人人| 久久99一区二区三区| 日本猛色少妇xxxxx猛交久久| 欧美精品亚洲一区二区| 黄色一级大片看看| 国产成人91sexporn| 老司机亚洲免费影院| 日韩电影二区| 欧美日韩精品成人综合77777| 在线观看免费高清a一片| 日日啪夜夜爽| 国产精品国产三级专区第一集| 日韩成人av中文字幕在线观看| 日韩制服丝袜自拍偷拍| 成人毛片a级毛片在线播放| 丰满乱子伦码专区| 伊人久久大香线蕉亚洲五| 香蕉丝袜av| 欧美精品国产亚洲| 最近2019中文字幕mv第一页| 99re6热这里在线精品视频| 精品国产乱码久久久久久小说| 亚洲欧美一区二区三区久久| 久久这里有精品视频免费| 这个男人来自地球电影免费观看 | 国产一级毛片在线| 亚洲一码二码三码区别大吗| 制服人妻中文乱码| 国产成人欧美| 亚洲欧美精品自产自拍| 校园人妻丝袜中文字幕| 国产免费又黄又爽又色| 婷婷色综合大香蕉| 亚洲精品日本国产第一区| 午夜福利乱码中文字幕| 亚洲美女视频黄频| 久久精品夜色国产| 精品少妇黑人巨大在线播放| 欧美日韩亚洲国产一区二区在线观看 | 国产精品偷伦视频观看了| 啦啦啦中文免费视频观看日本| 伊人久久国产一区二区| 国产精品偷伦视频观看了| 热re99久久国产66热| 少妇精品久久久久久久| 制服诱惑二区| 一级毛片我不卡| 精品一品国产午夜福利视频| 欧美bdsm另类| 九色亚洲精品在线播放| 国产男女超爽视频在线观看| 亚洲精品国产色婷婷电影| 飞空精品影院首页| 天天操日日干夜夜撸| 色视频在线一区二区三区| 亚洲精品一区蜜桃| 青草久久国产| 亚洲国产欧美网| 男人添女人高潮全过程视频| 最新中文字幕久久久久| 侵犯人妻中文字幕一二三四区| 国产精品二区激情视频| 日韩中文字幕视频在线看片| 国产欧美日韩一区二区三区在线| 天天操日日干夜夜撸| 黄色怎么调成土黄色| 日日啪夜夜爽| 午夜av观看不卡| 午夜91福利影院| 日韩一区二区视频免费看| 久久99蜜桃精品久久| av卡一久久| 青春草视频在线免费观看| 制服诱惑二区| 免费在线观看视频国产中文字幕亚洲 | 一本色道久久久久久精品综合| 在线免费观看不下载黄p国产| 在线天堂中文资源库| 亚洲精品乱久久久久久| 亚洲国产成人一精品久久久| 日韩一本色道免费dvd| 高清不卡的av网站| 亚洲视频免费观看视频| 亚洲国产欧美在线一区| 欧美成人午夜免费资源| 麻豆av在线久日| 久久韩国三级中文字幕| 国产精品不卡视频一区二区| 免费在线观看视频国产中文字幕亚洲 | 菩萨蛮人人尽说江南好唐韦庄| 成年人午夜在线观看视频| 99re6热这里在线精品视频| 一本色道久久久久久精品综合| 亚洲三区欧美一区| 永久免费av网站大全| 人人澡人人妻人| 黑人巨大精品欧美一区二区蜜桃| freevideosex欧美| 久久国产亚洲av麻豆专区| 久久精品熟女亚洲av麻豆精品| 男人操女人黄网站| 最近2019中文字幕mv第一页| 狠狠婷婷综合久久久久久88av| 国产一区二区 视频在线| 新久久久久国产一级毛片| 男女啪啪激烈高潮av片| 夜夜骑夜夜射夜夜干| 18禁国产床啪视频网站| 爱豆传媒免费全集在线观看| 狠狠精品人妻久久久久久综合| 男女免费视频国产| 大片免费播放器 马上看| 男人添女人高潮全过程视频| 久久久精品国产亚洲av高清涩受| 久久久久久免费高清国产稀缺| 国产av码专区亚洲av| 中文字幕最新亚洲高清| 亚洲第一av免费看| 亚洲欧洲国产日韩| 青春草国产在线视频| 国产av精品麻豆| 国产日韩一区二区三区精品不卡| 欧美av亚洲av综合av国产av | 成人18禁高潮啪啪吃奶动态图| 美女视频免费永久观看网站| www.av在线官网国产| 久久99蜜桃精品久久| 久久人人爽av亚洲精品天堂| 青春草国产在线视频| 欧美人与善性xxx| 欧美精品一区二区免费开放| 免费观看性生交大片5| 如何舔出高潮| av有码第一页| 久久这里只有精品19| 春色校园在线视频观看| 男女无遮挡免费网站观看| 精品人妻偷拍中文字幕| 一级毛片电影观看| 国产乱来视频区| 国产在线视频一区二区| 国产精品一区二区在线不卡| 亚洲国产看品久久| 国产精品秋霞免费鲁丝片| 人妻系列 视频| 一级毛片黄色毛片免费观看视频| 久久 成人 亚洲| 国产精品香港三级国产av潘金莲 | 成人手机av| 中文乱码字字幕精品一区二区三区| 免费大片黄手机在线观看| 久久久精品94久久精品| 精品视频人人做人人爽| 可以免费在线观看a视频的电影网站 | 国产xxxxx性猛交| 日韩精品有码人妻一区| 人妻系列 视频| 天堂中文最新版在线下载| 免费观看性生交大片5| 国产av一区二区精品久久| av片东京热男人的天堂| 午夜精品国产一区二区电影| 两性夫妻黄色片| 午夜激情av网站| 在线天堂最新版资源| 啦啦啦视频在线资源免费观看| 免费看不卡的av| 亚洲国产欧美在线一区| 亚洲伊人色综图| 亚洲熟女精品中文字幕| 久久久久久免费高清国产稀缺| 老鸭窝网址在线观看| 一区二区三区乱码不卡18| 国产在线视频一区二区| 人人妻人人添人人爽欧美一区卜| kizo精华| av一本久久久久| 国产xxxxx性猛交| 极品人妻少妇av视频| 亚洲精品在线美女| 2022亚洲国产成人精品| 免费黄色在线免费观看| 交换朋友夫妻互换小说| 男人爽女人下面视频在线观看| 国产精品一区二区在线观看99| 人人妻人人添人人爽欧美一区卜| 欧美黄色片欧美黄色片| 亚洲av电影在线观看一区二区三区| 青春草视频在线免费观看| 亚洲综合色惰| 久久精品国产a三级三级三级| 波多野结衣一区麻豆| 韩国av在线不卡| 国产女主播在线喷水免费视频网站| 天天躁日日躁夜夜躁夜夜| 国产成人精品福利久久| 亚洲情色 制服丝袜| 考比视频在线观看| 欧美少妇被猛烈插入视频| 丝袜喷水一区| 一本—道久久a久久精品蜜桃钙片| 久久久久久人妻| 街头女战士在线观看网站| 欧美日韩精品网址| 精品福利永久在线观看| 老司机影院成人| 国产男女超爽视频在线观看| 黑人欧美特级aaaaaa片| 日日啪夜夜爽| 免费观看av网站的网址| av福利片在线| 黑人猛操日本美女一级片| 精品一区在线观看国产| www.熟女人妻精品国产| 亚洲欧美清纯卡通| 一级爰片在线观看| 亚洲精品美女久久久久99蜜臀 | 成年人午夜在线观看视频| 亚洲国产看品久久| 母亲3免费完整高清在线观看 | av线在线观看网站| 精品福利永久在线观看| 久久 成人 亚洲| 一区二区日韩欧美中文字幕| 又大又黄又爽视频免费| 国产 精品1| 成人亚洲精品一区在线观看| www.熟女人妻精品国产| 国产在线一区二区三区精| 色94色欧美一区二区| 亚洲综合色网址| 王馨瑶露胸无遮挡在线观看| 中文字幕亚洲精品专区| 国产伦理片在线播放av一区| 亚洲国产精品一区二区三区在线| 日韩制服骚丝袜av| 欧美日韩视频高清一区二区三区二| 在现免费观看毛片| 丝袜脚勾引网站| 不卡视频在线观看欧美| 欧美最新免费一区二区三区| 18在线观看网站| 国产亚洲午夜精品一区二区久久| 久久午夜福利片| 成年动漫av网址| 老司机影院毛片| av天堂久久9| 久久国产精品男人的天堂亚洲| 久久精品久久久久久噜噜老黄| 日韩 亚洲 欧美在线| 波野结衣二区三区在线| 午夜激情久久久久久久| 啦啦啦啦在线视频资源| 国产福利在线免费观看视频| 国产免费一区二区三区四区乱码| 欧美黄色片欧美黄色片| 一边亲一边摸免费视频| 97在线人人人人妻| 欧美日韩成人在线一区二区| 欧美国产精品va在线观看不卡| 大香蕉久久网| 久久久久人妻精品一区果冻| 亚洲欧美色中文字幕在线| 免费在线观看黄色视频的| 侵犯人妻中文字幕一二三四区| 岛国毛片在线播放| 欧美中文综合在线视频| 亚洲伊人色综图| 成人手机av| 日本欧美国产在线视频| 精品人妻熟女毛片av久久网站| 成人国产麻豆网| 免费观看无遮挡的男女| 如何舔出高潮| 考比视频在线观看| 国产一区二区三区av在线| 亚洲精品日本国产第一区| 国产成人精品在线电影| 国产精品一区二区在线不卡| 最近2019中文字幕mv第一页| 日韩熟女老妇一区二区性免费视频| 成人国产麻豆网| 一级黄片播放器| 黄色一级大片看看| 日韩中文字幕视频在线看片| 欧美成人精品欧美一级黄| 男女下面插进去视频免费观看| 亚洲内射少妇av| 亚洲精品国产av蜜桃| 午夜福利网站1000一区二区三区| 国产高清国产精品国产三级| 99热全是精品| h视频一区二区三区| 日日撸夜夜添| 久久久久久久大尺度免费视频| 日韩av在线免费看完整版不卡| 高清欧美精品videossex| 国产精品.久久久| 中国国产av一级| 日韩精品免费视频一区二区三区| 色婷婷av一区二区三区视频| 婷婷色av中文字幕| 天天躁夜夜躁狠狠久久av| 中文字幕精品免费在线观看视频| 久久鲁丝午夜福利片| 青春草亚洲视频在线观看| 日韩电影二区| 各种免费的搞黄视频| 丝瓜视频免费看黄片| 亚洲激情五月婷婷啪啪| 亚洲欧美一区二区三区久久| 赤兔流量卡办理| 十分钟在线观看高清视频www| 婷婷色av中文字幕| 日韩av不卡免费在线播放| 亚洲国产最新在线播放| 99国产精品免费福利视频| 欧美日韩亚洲国产一区二区在线观看 | 在线 av 中文字幕| 人人妻人人澡人人爽人人夜夜| 欧美日韩视频高清一区二区三区二| 亚洲精品久久成人aⅴ小说| 菩萨蛮人人尽说江南好唐韦庄| 美女视频免费永久观看网站| 日韩av不卡免费在线播放| www日本在线高清视频| 午夜福利在线免费观看网站| 国产成人av激情在线播放| 1024香蕉在线观看| 国产精品三级大全| 免费少妇av软件| 搡女人真爽免费视频火全软件| 亚洲国产色片| 色网站视频免费| 国产一区亚洲一区在线观看| 97精品久久久久久久久久精品| 老汉色av国产亚洲站长工具| 久久这里只有精品19| 高清欧美精品videossex| 秋霞伦理黄片| 在线观看免费日韩欧美大片| 久久青草综合色| 国产精品香港三级国产av潘金莲 | 母亲3免费完整高清在线观看 | 女的被弄到高潮叫床怎么办| 丝袜美腿诱惑在线| 国产综合精华液| 国产av一区二区精品久久| 搡女人真爽免费视频火全软件| 中文字幕精品免费在线观看视频| 日本黄色日本黄色录像| 有码 亚洲区| 亚洲精品一区蜜桃| 亚洲一区二区三区欧美精品| 精品久久久精品久久久| 国产不卡av网站在线观看| 哪个播放器可以免费观看大片| 亚洲精品国产av成人精品| 国产午夜精品一二区理论片| 日本色播在线视频| 咕卡用的链子| av视频免费观看在线观看| 日韩中字成人| 青春草国产在线视频| 男女啪啪激烈高潮av片| 水蜜桃什么品种好| 一级片免费观看大全| 亚洲国产精品成人久久小说| 男女下面插进去视频免费观看| 亚洲四区av| 亚洲欧美精品自产自拍| av不卡在线播放| 亚洲欧美一区二区三区黑人 | 免费观看在线日韩| 国产综合精华液| 免费少妇av软件| 亚洲av国产av综合av卡| 欧美日韩国产mv在线观看视频| 免费在线观看黄色视频的| 久久久久精品性色| 精品少妇黑人巨大在线播放| 亚洲av男天堂| 丰满饥渴人妻一区二区三| 午夜福利在线观看免费完整高清在| 国产精品一区二区在线不卡| 伊人久久国产一区二区| 电影成人av| 大片免费播放器 马上看| 亚洲欧美成人精品一区二区| 少妇人妻久久综合中文| 一边摸一边做爽爽视频免费| 亚洲一码二码三码区别大吗| 久久精品久久精品一区二区三区| 国产精品不卡视频一区二区| 午夜久久久在线观看| 国产av一区二区精品久久| 欧美另类一区| 考比视频在线观看| 欧美精品国产亚洲| 黑人巨大精品欧美一区二区蜜桃| 成人黄色视频免费在线看| 国产日韩欧美视频二区| 在线看a的网站| 一级毛片我不卡| 日韩熟女老妇一区二区性免费视频| 一级毛片电影观看| 午夜影院在线不卡| 亚洲av男天堂| av线在线观看网站| 免费人妻精品一区二区三区视频| 亚洲av免费高清在线观看| 看十八女毛片水多多多| 久久影院123| 亚洲国产看品久久| 国产精品秋霞免费鲁丝片| 人妻系列 视频| 成人毛片a级毛片在线播放| 午夜91福利影院| 欧美日韩一区二区视频在线观看视频在线| av不卡在线播放| 伊人亚洲综合成人网| 免费黄色在线免费观看| 亚洲精品,欧美精品| 黄色 视频免费看| 日本wwww免费看| 99久久人妻综合| 日本色播在线视频| 免费看av在线观看网站| 久久久亚洲精品成人影院| 免费在线观看视频国产中文字幕亚洲 | 久久人人97超碰香蕉20202| 久久精品亚洲av国产电影网| 亚洲四区av| 最近中文字幕2019免费版| 看免费av毛片| 久久av网站| 国产高清国产精品国产三级| 男女无遮挡免费网站观看| 成人国产av品久久久| 亚洲国产精品一区三区| 国产精品二区激情视频| 久久97久久精品| 夫妻性生交免费视频一级片| 黄色毛片三级朝国网站| 久久久国产精品麻豆| 热re99久久国产66热| 国产免费现黄频在线看| 国产精品久久久久久av不卡| 久久久久精品人妻al黑| 综合色丁香网| 深夜精品福利| 青青草视频在线视频观看| 久久亚洲国产成人精品v| 亚洲国产精品一区三区| 国产成人精品福利久久| 一本久久精品| 波野结衣二区三区在线| 国产免费视频播放在线视频| 国产毛片在线视频| 亚洲av电影在线观看一区二区三区| 久久精品aⅴ一区二区三区四区 | 免费在线观看视频国产中文字幕亚洲 | 一区二区日韩欧美中文字幕| 啦啦啦在线观看免费高清www| 国产又爽黄色视频| 少妇的丰满在线观看| 午夜福利乱码中文字幕| 国产综合精华液| 夫妻性生交免费视频一级片| 国产老妇伦熟女老妇高清| 亚洲精品视频女| 18禁动态无遮挡网站| 日韩大片免费观看网站| 80岁老熟妇乱子伦牲交| 久久99热这里只频精品6学生| 黑丝袜美女国产一区| 国产精品久久久久久精品电影小说| 国产精品偷伦视频观看了| www.自偷自拍.com| 国产精品亚洲av一区麻豆 | 99国产精品免费福利视频| xxxhd国产人妻xxx| 亚洲国产日韩一区二区| 午夜福利视频精品| 色网站视频免费| 99久久精品国产国产毛片| 满18在线观看网站| 欧美日韩精品网址| 亚洲激情五月婷婷啪啪| 韩国高清视频一区二区三区| 制服诱惑二区| 久热久热在线精品观看| 日韩中字成人| 久久人人97超碰香蕉20202| 日本vs欧美在线观看视频| 青春草视频在线免费观看| www.熟女人妻精品国产| 午夜福利,免费看| 午夜福利视频精品| 精品福利永久在线观看| 亚洲国产欧美日韩在线播放| 精品久久蜜臀av无| 蜜桃在线观看..| 亚洲国产欧美网| 91午夜精品亚洲一区二区三区| 亚洲精品在线美女| 观看av在线不卡| 涩涩av久久男人的天堂| 中文精品一卡2卡3卡4更新|