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

    30P-30N多段翼型復雜流場數(shù)值模擬技術(shù)研究

    2010-11-08 06:18:40王運濤王光學張玉倫
    空氣動力學學報 2010年1期
    關(guān)鍵詞:主翼站位構(gòu)型

    王運濤,王光學,張玉倫

    (1.空氣動力學國家重點實驗室,四川綿陽 621000;2.中國空氣動力研究與發(fā)展中心,四川綿陽 621000)

    0 引 言

    隨著計算機硬件技術(shù)和CFD技術(shù)本身的發(fā)展,采用基于雷諾平均NS方程(RANS)的數(shù)值模擬軟件已經(jīng)可以模擬真實飛行器的復雜外形及全機的復雜流場,包括二維高升力翼型和三維帶增升裝置的全機構(gòu)型[1]?;赗ANS方程預(yù)測巡航構(gòu)型氣動特性變化趨勢和模擬全湍流附著流動的能力已經(jīng)逐步得到飛行器設(shè)計工程師的認可,但高升力構(gòu)型的數(shù)值模擬可信度水平依然很低。采用CFD手段尚很難準確模擬高升力構(gòu)型的最大升力系數(shù)及失速攻角,特別是對于存在明顯分離區(qū)的復雜流動,準確預(yù)測分離流動的開始和發(fā)展,以及雷諾數(shù)效應(yīng)依然是CFD的難點之一。

    為了研究高升力構(gòu)型的流動機理,提高CFD軟件的數(shù)值模擬精度,空氣動力學的試驗工作者和CFD工作者付出了巨大的努力[2-4],高升力構(gòu)型的數(shù)值模擬也是許多CFD可信度專題會議的主題。高升力構(gòu)型的主要模擬難點在于:邊界層轉(zhuǎn)捩、激波/邊界層干擾、粘性尾跡區(qū)干擾、尾跡與邊界層摻混和分離流動等等。盡管采用RANS方程模擬高升力構(gòu)型存在諸多困難,尤其是工程湍流模型的適用范圍眾說紛紜,但采用RANS方程和工程湍流模型依然是模擬飛行器復雜構(gòu)型的主要手段。Rumsey等人[5]對多位國外研究者的工作進行了綜述,研究表明,考慮氣動彈性影響和轉(zhuǎn)捩位置有助于提高多段翼型的數(shù)值模擬精度。隨著網(wǎng)格生成技術(shù)、大規(guī)模并行計算技術(shù)和工程湍流模型研究工作的不斷深入,采用RANS方程模擬高升力構(gòu)型的可信度水平有望達到一個新的高度。計算結(jié)果與試驗結(jié)果進行對比不應(yīng)該局限于升力、阻力等總體氣動特性,還應(yīng)該包含壓力分布和粘性剪切層數(shù)據(jù)的精細比較。

    本文通過求解任意坐標系下的雷諾平均的N-S方程,采用多塊對接結(jié)構(gòu)網(wǎng)格技術(shù),在與相應(yīng)試驗結(jié)果對比的基礎(chǔ)上,詳細研究了SA一方程湍流模型、SST兩方程湍流模型、不同的轉(zhuǎn)捩位置對該翼型壓力分布和典型站位速度型的影響。本文的研究結(jié)果表明,采用全湍流模擬方式可以較好地模擬該多段翼型的壓力分布,但對速度型的模擬精度較差;模擬試驗的轉(zhuǎn)捩位置可以改善主翼附面層與前緣縫翼邊界層發(fā)展的模擬精度;采用微吸氣技術(shù)推遲前緣縫翼的轉(zhuǎn)捩位置,可以進一步提高主翼上縫翼尾跡區(qū)的數(shù)值模擬精度。

    1 計算構(gòu)型和計算方法

    麥道航空公司的30P-30N三段增升構(gòu)型是被CFD工作者廣泛采用的多段構(gòu)型之一,該翼型的前緣縫翼和后緣襟翼的偏角均為30°,前緣縫翼的縫道寬度為2.95%,外伸量為-2.5%;后緣襟翼縫道寬度1.27%,外伸量為0.25%,是典型的著陸構(gòu)型。30P-30N三段翼型曾作為1993年NASA Langley舉行的CFD Challenge Workshop算例,其風洞試驗主要是90年代在NASA Langley的低湍流度增壓風洞中完成的,主要包括測壓試驗和邊界層測量等,試驗雷諾數(shù)5×106~16×106。試驗中采用了側(cè)壁吸氣方法來保證近似的二維流動特性,但試驗結(jié)果表明,攻角大于16°后,流動的三維效應(yīng)較強,這會對試驗與計算結(jié)果的比較產(chǎn)生一定的影響。計算構(gòu)型、速度型測量位置如圖1所示。本文模擬的狀態(tài)為:馬赫數(shù)M=0.2,攻角 α=19°,雷諾數(shù)Re=9×106。該攻角已接近失速攻角,這使得該狀態(tài)的流動模擬更具挑戰(zhàn)性。

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

    圖1 30P-30N構(gòu)型及速度型站位Fig.1 30P-30N profile locations

    圖2 30P-30N的計算網(wǎng)格(局部放大)Fig.2 30P-30N computational grid(local)

    本文的數(shù)值模擬軟件采用了中國空氣動力研究與發(fā)展中心自行研發(fā)的CFD軟件TRIP2.0(TRIsonic Platform Version2.0)。該軟件采用結(jié)構(gòu)網(wǎng)格技術(shù)和有限體積方法,通過數(shù)值求解三維任意坐標系下的Reynolds-Averaged Navier-Stokes(RANS)方程 ,獲得繞流復雜飛行器的繞流流場。在本文的數(shù)值模擬中,離散方程組的求解采用LU-SGS方法,無粘通量的離散選擇了三階精度的ROE格式,粘性通量的離散采用中心格式,湍流模型采用了SA一方程和SST兩方程模型,計算中采用了多重網(wǎng)格技術(shù)加速收斂。

    2 全湍流模擬結(jié)果與分析

    圖3給出了采用SA一方程模型和SST兩方程湍流模型壁面壓力分布的計算結(jié)果,同時給出了相應(yīng)的試驗結(jié)果。沒有考慮試驗的轉(zhuǎn)捩位置,湍流模型方程全流場求解。必須指出的是,采用全湍流模擬方式并不意味著流場駐點以后立刻產(chǎn)生湍流邊界層,一般而言采用全湍流模擬方式得到的數(shù)值解本身是按照湍流模型自身定義的轉(zhuǎn)捩位置而開始轉(zhuǎn)捩的,只不過與試驗的轉(zhuǎn)捩位置相比較大大提前[5]。

    圖3 表面壓力系數(shù)的比較(全湍流)Fig.3 Surface pressure coefficient(fully turbulent)

    由圖3中可以看出,采用兩種湍流模型得到的縫翼、主翼和襟翼上的壓力分布幾乎相同,并均與試驗結(jié)果吻合。兩種湍流模型的數(shù)值模擬結(jié)果僅在襟翼的上表面存在細微的差別,SA模型的計算結(jié)果略低于SST模型的計算結(jié)果。采用兩種湍流模型得到的總體氣動特性存在比較明顯的差別,兩種湍流模型得到的升力系數(shù)相差0.08,阻力系數(shù)相差0.0016,相比較而言,SST模型得到的升力系數(shù)更接近試驗結(jié)果。

    圖4給出了采用兩種湍流模型得到的六個典型站位上的速度型及相應(yīng)的試驗結(jié)果(具體站位見圖1)。由圖中可見,采用全湍流方式得到的典型站位上的速度型與試驗結(jié)果在定性上是接近的,定量上差別明顯。在主翼的三個站位上數(shù)值模擬結(jié)果與試驗結(jié)果的差別主要在于主翼的邊界層偏厚和縫翼尾跡區(qū)速度偏低、寬度偏大三個方面,但兩種湍流模型的結(jié)果彼此基本吻合。在襟翼的三個站位上,兩種湍流模型的計算結(jié)果差別明顯,在主翼尾跡區(qū)的模擬上,SST模型過低地預(yù)測了主翼尾跡區(qū)的速度,SA模型則與試驗結(jié)果吻合較好;兩種湍流模型均過低地預(yù)測了縫翼尾跡區(qū)的速度值、過高地預(yù)測了縫翼尾跡區(qū)的寬度,在最后一個站位上,SST模型的縫翼尾跡區(qū)的數(shù)值結(jié)果有所改善??傊?采用全湍流方式模擬該三段構(gòu)型,主翼邊界層厚度、縫翼的尾跡區(qū)和主翼的尾跡區(qū)的計算結(jié)果均與試驗結(jié)果差別較大,SA模型在模擬主翼尾跡區(qū)的結(jié)果略優(yōu),SST模型則在模擬縫翼尾跡區(qū)上結(jié)果略優(yōu)(襟翼站位上)。

    3 模擬試驗轉(zhuǎn)捩位置的結(jié)果與分析

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

    表1 計算與試驗的轉(zhuǎn)捩位置(α=19°)Table1 Transition location for calculation and test(α=19°)

    表面壓力分布的計算結(jié)果(本文沒有給出)表明,模擬試驗的轉(zhuǎn)捩位置得到的壓力分布與采用全湍流模擬方式得到的壓力分布基本一致,且均與試驗結(jié)果吻合良好,只是在襟翼的上表面模擬試驗轉(zhuǎn)捩位置后得到的壓力系數(shù)略高于全湍流計算的結(jié)果和試驗值,從本文以下的討論中可以看出,這主要是由于襟翼上方邊界層的摻混模擬的差別造成的。

    圖5給出了采用SST兩方程模型是否模擬試驗轉(zhuǎn)捩位置的典型站位速度型的比較。模擬了試驗的轉(zhuǎn)捩位置后,典型站位的速度型得到了明顯的改善。首先是主翼上邊界層的厚度和與試驗結(jié)果更加接近,其次是縫翼和主翼的尾跡區(qū)速度值和尾跡區(qū)的寬度均有明顯改善,但與試驗結(jié)果相比較還存在一定差距。

    表2給出了采用兩種湍流模型和是否模擬試驗轉(zhuǎn)捩位置得到的升力系數(shù)、阻力系數(shù)以及相應(yīng)的試驗值??梢钥吹?模擬了試驗的轉(zhuǎn)捩位置后,升力系數(shù)略有增加,阻力系數(shù)比全湍流方式降低了100個阻力單位(1阻力單位=0.0001)。

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

    表2 氣動特性的比較(α=19°)Table2 Aerodynamic character(α=19°)

    4 延遲轉(zhuǎn)捩位置的計算結(jié)果與分析

    從第3節(jié)的討論中可以看到,模擬了試驗的轉(zhuǎn)捩位置后,雖然縫翼的尾跡區(qū)模擬精度有了比較明顯的改善,但與試驗結(jié)果之間依然存在較大的差距,本文作者認為縫翼上表面的邊界層的產(chǎn)生和發(fā)展是決定其尾跡區(qū)模擬精度的主要原因之一。本節(jié)將討論縫翼上表面轉(zhuǎn)捩位置變化對其尾跡區(qū)發(fā)展的影響。

    在數(shù)值模擬過程發(fā)現(xiàn),將試驗測定的縫翼上表面轉(zhuǎn)捩位置進一步向下游方向移動而保持轉(zhuǎn)捩區(qū)長度不變,縫翼上表面的逆力梯度會導致縫翼上表面頭部出現(xiàn)較大范圍的分離區(qū)。為了避免縫翼頭部的分離,本文在數(shù)值模擬中采用了微吸氣技術(shù),微吸氣的起始位置在縫翼頭部頂點,終止位置在確定的轉(zhuǎn)捩區(qū)之前。采用上述技術(shù),本文將縫翼上表面的起始轉(zhuǎn)捩位置向下游推遲到-0.073的位置(轉(zhuǎn)捩位置B),而保持轉(zhuǎn)捩區(qū)長度不變,其它的轉(zhuǎn)捩參數(shù)與轉(zhuǎn)捩位置A相同。

    圖6給出了采用SST兩方程模型和轉(zhuǎn)捩位置B的數(shù)值模擬結(jié)果,圖中同時還給出了全湍流和轉(zhuǎn)捩位置A的數(shù)值模擬結(jié)果。由圖中可以看出,將轉(zhuǎn)捩位置推遲到B點后,縫翼尾跡區(qū)的模擬精度有了進一步的提高,尤其是在圖1中所示的前四個站位上(在x/c=0.1075站位上,試驗結(jié)果可能存在誤差[6])。本節(jié)和第三節(jié)的計算結(jié)果表明,SST模型具有較好的尾跡區(qū)模擬精度,導致計算與試驗縫翼尾跡區(qū)速度型差異的主要原因不是湍流模型本身,而是縫翼上表面湍流邊界層起始位置。

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

    5 結(jié) 論

    本文采用TRIP2.0軟件和結(jié)構(gòu)對接網(wǎng)格技術(shù),通過求解任意坐標系下的RANS方程,數(shù)值模擬了30P-30N三段翼型的復雜流場,主要研究了SA和SST兩種湍流模型、轉(zhuǎn)捩位置對數(shù)值模擬結(jié)果的影響。通過與相應(yīng)的試驗結(jié)果相比較,得到以下一些基本結(jié)論:

    (1)采用全湍流模擬方式,兩種湍流模型均可以較好地模擬該多段翼型的壓力分布,但速度型的計算結(jié)果與試驗結(jié)果差別明顯;

    (2)模擬試驗的轉(zhuǎn)捩位置,對壓力分布的計算結(jié)果影響較小,對氣動力系數(shù),尤其是阻力系數(shù)影響顯著;可以明顯提高速度型的數(shù)值模擬精度,但縫翼尾跡區(qū)的計算精度仍然與試驗結(jié)果存在明顯差距;

    (3)采用微吸氣技術(shù)推遲縫翼的轉(zhuǎn)捩位置,可以進一步提高縫翼尾跡區(qū)的數(shù)值模擬精度;

    (4)本文采用的湍流模型本身具有較強的尾跡模擬能力,轉(zhuǎn)捩位置是導致計算與試驗在速度型上差異的主要因素之一。

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

    [2]PAUL L JOHNSON,KENNETH M JONES,MICHAEL DMADSON.Experimental investigation of a simplified 3D high 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]朱自強,陳迎春,吳宗成等.高升力系統(tǒng)外形的數(shù)值模擬計算[J].航空學報,2005,26(3):257-262.

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

    [6]CHRISTOPHER L RUMSEY,THOMAS BGATSKI,SUSAN X YING,et al.Prediction of high-lift f lows using turbulent closure models[R].AIAA-97-2260.

    猜你喜歡
    主翼站位構(gòu)型
    多段翼低雷諾數(shù)繞流渦-邊界層相互干擾
    航空學報(2023年12期)2023-07-28 10:44:12
    三黃雞
    某型民機低速巡航構(gòu)型平尾抖振特性風洞試驗研究
    提高政治站位 對標國內(nèi)一流
    建黨百年說“站位”
    水上消防(2021年3期)2021-08-21 03:12:00
    分子和離子立體構(gòu)型的判定
    人大建設(shè)(2019年5期)2019-10-08 08:55:10
    探究鴨式布局模型飛機
    航空模型(2017年3期)2017-07-28 21:22:01
    航天器受迫繞飛構(gòu)型設(shè)計與控制
    遙感衛(wèi)星平臺與載荷一體化構(gòu)型
    亚洲 国产 在线| 免费久久久久久久精品成人欧美视频| 这个男人来自地球电影免费观看| 亚洲一码二码三码区别大吗| 欧美日韩中文字幕国产精品一区二区三区 | 国产av又大| 成人手机av| 国产成人影院久久av| 午夜福利视频1000在线观看 | 曰老女人黄片| 欧美成人午夜精品| 欧美一区二区精品小视频在线| 999久久久国产精品视频| 久久人人精品亚洲av| 国产av在哪里看| 国产蜜桃级精品一区二区三区| 亚洲av美国av| 黑丝袜美女国产一区| 十八禁人妻一区二区| 精品国产一区二区久久| www.精华液| 精品第一国产精品| 久久久久久久久中文| 一边摸一边抽搐一进一出视频| 国产视频一区二区在线看| 12—13女人毛片做爰片一| 免费看a级黄色片| 国产av一区二区精品久久| 亚洲精品美女久久久久99蜜臀| 色综合站精品国产| 狠狠狠狠99中文字幕| 在线观看免费日韩欧美大片| 久99久视频精品免费| 99久久综合精品五月天人人| 国产日韩一区二区三区精品不卡| 黄色毛片三级朝国网站| 亚洲自偷自拍图片 自拍| 9191精品国产免费久久| 亚洲国产精品合色在线| 亚洲中文字幕日韩| 波多野结衣高清无吗| 一级a爱片免费观看的视频| 1024香蕉在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 久久久水蜜桃国产精品网| 中文字幕另类日韩欧美亚洲嫩草| 日韩大码丰满熟妇| 国产欧美日韩综合在线一区二区| 久久精品国产99精品国产亚洲性色 | av天堂久久9| 真人做人爱边吃奶动态| 超碰成人久久| 成人特级黄色片久久久久久久| 伊人久久大香线蕉亚洲五| 国产又爽黄色视频| 在线免费观看的www视频| 国产精品爽爽va在线观看网站 | 免费在线观看完整版高清| 91成年电影在线观看| 满18在线观看网站| 精品欧美国产一区二区三| 精品一区二区三区四区五区乱码| 少妇的丰满在线观看| 在线播放国产精品三级| 丁香六月欧美| 国产乱人伦免费视频| 日韩精品青青久久久久久| 国产主播在线观看一区二区| 欧美黑人欧美精品刺激| 午夜亚洲福利在线播放| 中亚洲国语对白在线视频| 午夜精品久久久久久毛片777| 国产亚洲欧美精品永久| 看免费av毛片| 日韩精品中文字幕看吧| 亚洲精品中文字幕在线视频| 国产真人三级小视频在线观看| 天堂√8在线中文| 99国产精品一区二区蜜桃av| 日本一区二区免费在线视频| 亚洲精品国产一区二区精华液| 人妻丰满熟妇av一区二区三区| 搞女人的毛片| 99久久久亚洲精品蜜臀av| 国产亚洲精品第一综合不卡| 麻豆成人av在线观看| 久久国产精品影院| 乱人伦中国视频| 国产私拍福利视频在线观看| 极品教师在线免费播放| 久久性视频一级片| 欧美丝袜亚洲另类 | 欧美日韩亚洲国产一区二区在线观看| 亚洲情色 制服丝袜| 久久精品国产综合久久久| 国产亚洲av高清不卡| 精品久久久久久,| 日本精品一区二区三区蜜桃| 久久中文字幕人妻熟女| 男人的好看免费观看在线视频 | 亚洲成av片中文字幕在线观看| 少妇裸体淫交视频免费看高清 | 在线观看免费视频日本深夜| 久久久久久人人人人人| 一进一出抽搐gif免费好疼| 国产精品香港三级国产av潘金莲| 国产成年人精品一区二区| 亚洲av第一区精品v没综合| 欧美人与性动交α欧美精品济南到| 精品久久久久久成人av| 丝袜美足系列| 成人特级黄色片久久久久久久| 国产男靠女视频免费网站| 国产精品免费一区二区三区在线| 免费观看人在逋| 欧美色视频一区免费| 亚洲avbb在线观看| 亚洲av片天天在线观看| 亚洲精品在线美女| 两个人免费观看高清视频| 亚洲九九香蕉| 好男人在线观看高清免费视频 | 日韩欧美国产在线观看| 日日爽夜夜爽网站| 午夜激情av网站| 纯流量卡能插随身wifi吗| 亚洲片人在线观看| 黄色视频不卡| 高清黄色对白视频在线免费看| 精品国产一区二区久久| 欧美国产日韩亚洲一区| 黑人操中国人逼视频| 亚洲人成网站在线播放欧美日韩| 亚洲自拍偷在线| 欧美日本中文国产一区发布| 性色av乱码一区二区三区2| 一进一出抽搐动态| 国产精品 欧美亚洲| 9191精品国产免费久久| 中文亚洲av片在线观看爽| 欧美在线黄色| 亚洲av日韩精品久久久久久密| 久久久久国内视频| 久久久国产成人精品二区| 人人澡人人妻人| 亚洲va日本ⅴa欧美va伊人久久| 97碰自拍视频| 国产高清激情床上av| 成人亚洲精品一区在线观看| 十八禁人妻一区二区| 成人亚洲精品av一区二区| 午夜日韩欧美国产| 久久精品91蜜桃| 色在线成人网| АⅤ资源中文在线天堂| 69av精品久久久久久| 国产成人av教育| 国产精品一区二区免费欧美| 一边摸一边抽搐一进一出视频| 精品人妻1区二区| 日韩欧美一区视频在线观看| 黄色视频,在线免费观看| 亚洲黑人精品在线| 成人国语在线视频| 黄色a级毛片大全视频| 九色亚洲精品在线播放| 美女午夜性视频免费| 久久中文字幕人妻熟女| 麻豆国产av国片精品| 在线播放国产精品三级| 色婷婷久久久亚洲欧美| 大型黄色视频在线免费观看| 欧美黑人欧美精品刺激| 丝袜人妻中文字幕| 大码成人一级视频| 一级毛片精品| 黄色 视频免费看| 天堂动漫精品| 成人手机av| 亚洲精品国产色婷婷电影| 国产乱人伦免费视频| 巨乳人妻的诱惑在线观看| 日韩欧美国产一区二区入口| 亚洲 国产 在线| 一二三四社区在线视频社区8| av超薄肉色丝袜交足视频| 99久久精品国产亚洲精品| 50天的宝宝边吃奶边哭怎么回事| 亚洲色图av天堂| 亚洲国产毛片av蜜桃av| 91精品三级在线观看| 亚洲五月婷婷丁香| 黄色丝袜av网址大全| 久热爱精品视频在线9| 91成人精品电影| 日日干狠狠操夜夜爽| 大码成人一级视频| 手机成人av网站| 啦啦啦免费观看视频1| 最近最新中文字幕大全免费视频| 成人国语在线视频| 亚洲国产日韩欧美精品在线观看 | 我的亚洲天堂| 亚洲欧美精品综合一区二区三区| 国产真人三级小视频在线观看| 国产精品av久久久久免费| 国产欧美日韩综合在线一区二区| 免费观看精品视频网站| 精品国产超薄肉色丝袜足j| 欧美老熟妇乱子伦牲交| 午夜福利在线观看吧| 久久久久精品国产欧美久久久| 欧美精品亚洲一区二区| 美女午夜性视频免费| 人妻丰满熟妇av一区二区三区| 999久久久精品免费观看国产| 午夜视频精品福利| 国产主播在线观看一区二区| 大型av网站在线播放| 多毛熟女@视频| 久久 成人 亚洲| 国产乱人伦免费视频| av视频免费观看在线观看| 99在线人妻在线中文字幕| 久久久久久久精品吃奶| 免费久久久久久久精品成人欧美视频| 色老头精品视频在线观看| 一区福利在线观看| 色精品久久人妻99蜜桃| 91大片在线观看| 成人av一区二区三区在线看| 黄色视频不卡| 99久久国产精品久久久| 看片在线看免费视频| 老司机靠b影院| 国产一区二区激情短视频| 最新在线观看一区二区三区| 久久精品亚洲精品国产色婷小说| 在线观看一区二区三区| 一级,二级,三级黄色视频| 国内久久婷婷六月综合欲色啪| 亚洲精品一卡2卡三卡4卡5卡| 欧美午夜高清在线| www.www免费av| 黄色a级毛片大全视频| 悠悠久久av| 午夜福利免费观看在线| 亚洲人成伊人成综合网2020| 丝袜美足系列| 最好的美女福利视频网| 欧美成人性av电影在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 一级毛片高清免费大全| 波多野结衣av一区二区av| 免费人成视频x8x8入口观看| 午夜久久久在线观看| 色老头精品视频在线观看| 男人舔女人下体高潮全视频| 日韩av在线大香蕉| www.熟女人妻精品国产| 午夜两性在线视频| 在线免费观看的www视频| 亚洲色图 男人天堂 中文字幕| 中文字幕久久专区| 久久婷婷成人综合色麻豆| 琪琪午夜伦伦电影理论片6080| 国产麻豆69| 亚洲国产欧美网| 国产真人三级小视频在线观看| a级毛片在线看网站| 国产日韩一区二区三区精品不卡| 丝袜在线中文字幕| 免费高清在线观看日韩| av片东京热男人的天堂| 免费在线观看视频国产中文字幕亚洲| 日韩成人在线观看一区二区三区| 日本a在线网址| 欧美中文综合在线视频| 精品不卡国产一区二区三区| 99精品欧美一区二区三区四区| 操出白浆在线播放| 久久精品国产综合久久久| 国产xxxxx性猛交| 免费看十八禁软件| 精品一区二区三区av网在线观看| 老汉色∧v一级毛片| 91成年电影在线观看| 久久草成人影院| 亚洲精品一卡2卡三卡4卡5卡| 亚洲男人天堂网一区| 国产av又大| 在线国产一区二区在线| 成人三级黄色视频| 电影成人av| 黑人欧美特级aaaaaa片| 一区二区三区高清视频在线| 90打野战视频偷拍视频| 欧美成人免费av一区二区三区| 国产视频一区二区在线看| 亚洲熟妇熟女久久| 99久久精品国产亚洲精品| 色老头精品视频在线观看| 国产私拍福利视频在线观看| 一级作爱视频免费观看| 女人爽到高潮嗷嗷叫在线视频| 国产不卡一卡二| 欧美成狂野欧美在线观看| 亚洲国产欧美网| 国产一区在线观看成人免费| 一边摸一边抽搐一进一小说| 窝窝影院91人妻| 国产精品久久久av美女十八| 丝袜人妻中文字幕| 中文字幕高清在线视频| 久久久国产欧美日韩av| 亚洲av熟女| 在线播放国产精品三级| 亚洲成人久久性| 国产伦人伦偷精品视频| 国产精品永久免费网站| 老汉色∧v一级毛片| 天天添夜夜摸| svipshipincom国产片| 一级a爱视频在线免费观看| 97人妻精品一区二区三区麻豆 | 欧美老熟妇乱子伦牲交| 可以在线观看的亚洲视频| 国产免费男女视频| 一进一出好大好爽视频| 真人一进一出gif抽搐免费| 亚洲国产日韩欧美精品在线观看 | 免费无遮挡裸体视频| 国产欧美日韩一区二区三| 国产成人av教育| 亚洲av五月六月丁香网| 99国产精品免费福利视频| 久久精品国产99精品国产亚洲性色 | 老司机在亚洲福利影院| АⅤ资源中文在线天堂| 夜夜夜夜夜久久久久| 少妇 在线观看| 国产精品爽爽va在线观看网站 | 香蕉国产在线看| netflix在线观看网站| 欧美一区二区精品小视频在线| 国产一区二区在线av高清观看| 电影成人av| 正在播放国产对白刺激| 伊人久久大香线蕉亚洲五| а√天堂www在线а√下载| 午夜免费鲁丝| 亚洲美女黄片视频| 亚洲精品久久国产高清桃花| 看片在线看免费视频| 亚洲片人在线观看| 亚洲成av人片免费观看| 看黄色毛片网站| 久久精品国产亚洲av香蕉五月| 热99re8久久精品国产| 欧美不卡视频在线免费观看 | 无限看片的www在线观看| 一进一出好大好爽视频| 久久久精品国产亚洲av高清涩受| 欧美不卡视频在线免费观看 | 久久国产精品影院| 露出奶头的视频| 日韩欧美在线二视频| 999久久久国产精品视频| 亚洲黑人精品在线| 丁香六月欧美| 国产主播在线观看一区二区| 久99久视频精品免费| 成人国产综合亚洲| 亚洲精品美女久久av网站| 老司机午夜福利在线观看视频| 亚洲一区二区三区色噜噜| 狠狠狠狠99中文字幕| 精品国内亚洲2022精品成人| 亚洲第一av免费看| 午夜精品久久久久久毛片777| 一边摸一边抽搐一进一出视频| 十八禁网站免费在线| 精品日产1卡2卡| 女人被躁到高潮嗷嗷叫费观| 亚洲aⅴ乱码一区二区在线播放 | 一个人观看的视频www高清免费观看 | 成年人黄色毛片网站| 熟妇人妻久久中文字幕3abv| 午夜免费观看网址| 免费不卡黄色视频| 亚洲成人免费电影在线观看| 18禁美女被吸乳视频| 中文字幕久久专区| 国产高清激情床上av| 日本vs欧美在线观看视频| 国产精品1区2区在线观看.| or卡值多少钱| 久久人妻av系列| 美女大奶头视频| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲午夜精品一区,二区,三区| 国产精品日韩av在线免费观看 | 免费在线观看完整版高清| 性欧美人与动物交配| 亚洲一区二区三区色噜噜| 久久亚洲真实| 国产成人精品久久二区二区免费| 国产av一区二区精品久久| x7x7x7水蜜桃| 欧美日韩黄片免| 免费人成视频x8x8入口观看| 看黄色毛片网站| 国产私拍福利视频在线观看| 老司机午夜福利在线观看视频| 日本免费一区二区三区高清不卡 | 久久国产精品影院| 国内毛片毛片毛片毛片毛片| 精品久久久久久久久久免费视频| 欧美在线黄色| 日本一区二区免费在线视频| 久久国产精品影院| 亚洲av电影在线进入| 夜夜躁狠狠躁天天躁| 久热这里只有精品99| 国产成人欧美在线观看| 国产高清视频在线播放一区| 亚洲五月天丁香| 欧美 亚洲 国产 日韩一| 亚洲精品粉嫩美女一区| 黄色丝袜av网址大全| 美女免费视频网站| 可以在线观看的亚洲视频| 久热爱精品视频在线9| 日韩欧美免费精品| 成人欧美大片| 亚洲片人在线观看| 十分钟在线观看高清视频www| 老司机靠b影院| a在线观看视频网站| 色老头精品视频在线观看| 日韩视频一区二区在线观看| 夜夜躁狠狠躁天天躁| 久热这里只有精品99| 一卡2卡三卡四卡精品乱码亚洲| 久久人人97超碰香蕉20202| 国产成人影院久久av| 亚洲自偷自拍图片 自拍| videosex国产| 国产99久久九九免费精品| 色综合亚洲欧美另类图片| 成熟少妇高潮喷水视频| 成人永久免费在线观看视频| 日韩欧美国产在线观看| 午夜福利18| 日韩有码中文字幕| 精品国产超薄肉色丝袜足j| 一级片免费观看大全| 好看av亚洲va欧美ⅴa在| 女生性感内裤真人,穿戴方法视频| 黑人欧美特级aaaaaa片| 日韩欧美一区二区三区在线观看| 久久久国产成人精品二区| 亚洲伊人色综图| 免费在线观看黄色视频的| 日韩中文字幕欧美一区二区| 人成视频在线观看免费观看| 欧美绝顶高潮抽搐喷水| 亚洲性夜色夜夜综合| 亚洲 欧美一区二区三区| 日韩 欧美 亚洲 中文字幕| 两人在一起打扑克的视频| 亚洲 国产 在线| 别揉我奶头~嗯~啊~动态视频| 少妇 在线观看| 色播亚洲综合网| 亚洲最大成人中文| 少妇裸体淫交视频免费看高清 | 12—13女人毛片做爰片一| 国产亚洲av嫩草精品影院| 麻豆成人av在线观看| 色播亚洲综合网| 欧美乱色亚洲激情| 欧美日韩亚洲综合一区二区三区_| 午夜日韩欧美国产| 亚洲精品美女久久久久99蜜臀| 人人澡人人妻人| 宅男免费午夜| 久久午夜综合久久蜜桃| 黄色女人牲交| 多毛熟女@视频| 在线观看免费视频网站a站| 久久草成人影院| 亚洲欧美日韩无卡精品| 亚洲人成伊人成综合网2020| 纯流量卡能插随身wifi吗| 别揉我奶头~嗯~啊~动态视频| 国产99白浆流出| 午夜精品在线福利| 中文亚洲av片在线观看爽| 又大又爽又粗| 国产视频一区二区在线看| 久久久久国内视频| 国产乱人伦免费视频| 亚洲国产欧美一区二区综合| 国语自产精品视频在线第100页| 久久亚洲真实| 欧美不卡视频在线免费观看 | 人人妻,人人澡人人爽秒播| 久9热在线精品视频| 身体一侧抽搐| 一区二区三区精品91| 黑人巨大精品欧美一区二区mp4| 日本 欧美在线| 一二三四在线观看免费中文在| 日本 欧美在线| 精品福利观看| 色精品久久人妻99蜜桃| 亚洲欧美精品综合一区二区三区| 电影成人av| 侵犯人妻中文字幕一二三四区| 中亚洲国语对白在线视频| 变态另类成人亚洲欧美熟女 | 一级黄色大片毛片| 18禁裸乳无遮挡免费网站照片 | 在线观看舔阴道视频| 露出奶头的视频| 国产亚洲精品av在线| 男女之事视频高清在线观看| 少妇熟女aⅴ在线视频| 99国产精品一区二区蜜桃av| 美女国产高潮福利片在线看| 欧美成人性av电影在线观看| 中文字幕精品免费在线观看视频| 少妇裸体淫交视频免费看高清 | 一夜夜www| 少妇熟女aⅴ在线视频| 日日干狠狠操夜夜爽| 老司机在亚洲福利影院| 国产精品久久久久久人妻精品电影| 手机成人av网站| 最近最新免费中文字幕在线| 久热爱精品视频在线9| 久久久久九九精品影院| 午夜两性在线视频| tocl精华| 亚洲成人久久性| 日韩有码中文字幕| 日韩欧美国产一区二区入口| 操美女的视频在线观看| 无遮挡黄片免费观看| 狂野欧美激情性xxxx| 一区二区日韩欧美中文字幕| 91大片在线观看| 国产精品九九99| 在线国产一区二区在线| 免费观看人在逋| 美女国产高潮福利片在线看| 亚洲精华国产精华精| 国产99白浆流出| 国产欧美日韩综合在线一区二区| 乱人伦中国视频| 久久草成人影院| 日韩欧美一区视频在线观看| 级片在线观看| 欧美激情极品国产一区二区三区| 丝袜美足系列| 日日干狠狠操夜夜爽| 国产精品98久久久久久宅男小说| 久久久久九九精品影院| 一区二区三区高清视频在线| aaaaa片日本免费| 丁香欧美五月| 日本在线视频免费播放| 精品乱码久久久久久99久播| 看免费av毛片| 国产成人精品在线电影| 国产成+人综合+亚洲专区| 久久久久久久精品吃奶| 夜夜看夜夜爽夜夜摸| 久99久视频精品免费| 国产野战对白在线观看| 校园春色视频在线观看| 国产一卡二卡三卡精品| 免费在线观看日本一区| 国产av精品麻豆| 亚洲av五月六月丁香网| 桃红色精品国产亚洲av| 亚洲成人久久性| 免费在线观看黄色视频的| 人人妻人人澡欧美一区二区 | 精品久久久久久,| 日韩欧美在线二视频| 美女高潮喷水抽搐中文字幕| 99国产精品99久久久久| 黄频高清免费视频| 一级a爱视频在线免费观看| 欧美成狂野欧美在线观看| 亚洲全国av大片| 首页视频小说图片口味搜索| 日韩欧美一区视频在线观看| 人人澡人人妻人| 亚洲av电影在线进入| 亚洲国产精品合色在线| 禁无遮挡网站| 亚洲精品在线观看二区| 一卡2卡三卡四卡精品乱码亚洲| 美女高潮喷水抽搐中文字幕| 成人精品一区二区免费| 久久久久久久精品吃奶| 精品一区二区三区视频在线观看免费| 三级毛片av免费| 欧美在线黄色| 我的亚洲天堂| 日韩欧美免费精品| 一级黄色大片毛片|