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

    基于網(wǎng)格框架的非結(jié)構(gòu)附面層網(wǎng)格生成技術(shù)

    2022-01-06 02:47:18龐宇飛肖素梅何雨陽謝冬香
    計算力學(xué)學(xué)報 2021年6期
    關(guān)鍵詞:附面層邊界點棱柱

    王 碩, 龐宇飛, 肖素梅, 劉 楊, 齊 龍, 何雨陽, 謝冬香

    (1.西南科技大學(xué) 制造科學(xué)與工程學(xué)院,綿陽 621010; 2.中國空氣動力研究與發(fā)展中心 空氣動力學(xué)國家重點實驗室,綿陽 621000)

    1 引 言

    隨著計算機(jī)科學(xué)的發(fā)展,計算力大幅度提高,使得計算流體力學(xué)(computational fluid danamics)在各個行業(yè)中變得越來越重要,而網(wǎng)格生成在CFD中扮演著重要的角色,且是數(shù)值計算的重要步驟。網(wǎng)格中附面層網(wǎng)格是為了更好地捕捉近物面薄層的流動物理特性,影響著CFD計算的精度。目前生成附面層網(wǎng)格的自動性及可靠性較差,而且需要大量的人工交互來完成,使得附面層網(wǎng)格的生成費時費力。近年來,很多學(xué)者采用混合網(wǎng)格的策略生成半結(jié)構(gòu)化的棱柱網(wǎng)格,以完成附面層網(wǎng)格的生成[1-5]。

    前沿層進(jìn)法[6-10]為經(jīng)典附面層網(wǎng)格生成方法。該方法先選取物面網(wǎng)格為前沿層,然后基于該前沿層逐層沿法向推進(jìn),生成的網(wǎng)格作為前沿層,繼續(xù)重復(fù)上述操作。文獻(xiàn)[11]基于前沿層進(jìn)法,發(fā)展了一套附面層網(wǎng)格生成方法,其中包括前沿法向和附面層生成高度的調(diào)整方法。

    張來平等[12]針對附面層三棱柱網(wǎng)格進(jìn)行了研究,提出了一種聚合的方法。首先生成各向異性的四面體網(wǎng)格,然后再將四面體網(wǎng)格進(jìn)行聚合,從而生成三棱柱網(wǎng)格。但是聚合后會有少數(shù)質(zhì)量較差的網(wǎng)格單元。孫巖[13]提出一種交互式棱柱網(wǎng)格生成方法,通過人工交互生成附面層邊界點的空間推進(jìn)網(wǎng)格,然后利用徑向基函數(shù)插值方法,得到附面層每一層內(nèi)部的網(wǎng)格點。該方法利用人工交互生成邊界點作為約束條件,有效避免了附面層網(wǎng)格交叉,并可以對局部網(wǎng)格進(jìn)行調(diào)整,最終生成高質(zhì)量的棱柱網(wǎng)格。但該方法需要人工交互生成附面層網(wǎng)格的邊界,且附面層推進(jìn)過程中會存在一定的誤差,需要后期進(jìn)行優(yōu)化。孫巖等[14]提出了一種基于約束框架的棱柱網(wǎng)格生成方法。首先,利用TFI(Transfinite Interpolation)生成粗背景結(jié)構(gòu)框架網(wǎng)格,再通過背景網(wǎng)格和表面網(wǎng)格的對應(yīng)關(guān)系插值得到附面層網(wǎng)格。該方法在附面層網(wǎng)格質(zhì)量控制和局部修改方面有較高的效率和顯著的優(yōu)勢,但背景框架網(wǎng)格生成比較困難。Lu等[15]提出一種結(jié)構(gòu)多塊附面層網(wǎng)格自動生成技術(shù),是一種由外而內(nèi)的附面層網(wǎng)格生成技術(shù)。首先,從表面網(wǎng)格中提取幾何特征,再基于幾何特征構(gòu)造附面層網(wǎng)格框架,然后在框架約束內(nèi)通過TFI生成結(jié)構(gòu)附面層網(wǎng)格。該方法能夠快速生成結(jié)構(gòu)附面層網(wǎng)格,但結(jié)構(gòu)網(wǎng)格的生成由于需要進(jìn)行網(wǎng)格拓?fù)湓O(shè)計,不僅費時,且對軟件使用者的經(jīng)驗有一定要求。

    本文在前文的研究基礎(chǔ)上,針對非結(jié)構(gòu)附面層網(wǎng)格生成的自動性和局部網(wǎng)格質(zhì)量控制能力進(jìn)行研究,提出了一種棱柱附面層網(wǎng)格生成方法。

    2 附面層棱柱網(wǎng)格生成方法

    本文發(fā)展了文獻(xiàn)[15]的以分而治之的附面層網(wǎng)格生成思想,用構(gòu)造的方法生成非結(jié)構(gòu)附面層網(wǎng)格。該方法生成的網(wǎng)格精度和全局一致性較高,能夠有效避免復(fù)雜外形附面層網(wǎng)格局部及全局交叉現(xiàn)象。

    2.1 方法思路

    引進(jìn)自動化的表面網(wǎng)格分片算法[16],自動將表面網(wǎng)格分成區(qū)域可控的多個網(wǎng)格面片,并以此生成多塊網(wǎng)格框架;利用結(jié)構(gòu)網(wǎng)格附面層自動生成技術(shù)[15]的框架線算法,并加以改進(jìn),利用徑向基函數(shù)插值[17]方法生成框架線中的輪廓線;在框架線的基礎(chǔ)上,利用徑向基函數(shù)插值方法構(gòu)造附面層頂面網(wǎng)格;通過線性插值生成附面層頂層網(wǎng)格點與表面網(wǎng)格點之間的附面層內(nèi)部網(wǎng)格點,進(jìn)而生成附面層網(wǎng)格。避免了文獻(xiàn)[13]在外形變化劇烈的復(fù)雜條件下,精細(xì)附面層網(wǎng)格生成的插值精度擾動抑制難題。圖1為附面層網(wǎng)格生成算法的流程。

    圖1 附面層網(wǎng)格生成示意圖

    2.2 徑向基函數(shù)插值

    徑向基函數(shù)插值在動網(wǎng)格中已經(jīng)有較好應(yīng)用。該方法通過選取基點之間的距離構(gòu)建插值基,能夠適應(yīng)復(fù)雜邊界。徑向基函數(shù)插值的一般表現(xiàn)形式為[18]

    (1)

    根據(jù)式(1)得到關(guān)于插值系數(shù)αi的線性方程組,表示成矩陣的形式為

    (2)

    求解矩陣方程(2)可以得到插值系數(shù)αi。然后通過式(1)能得到任一點r的函數(shù)值及位移。

    2.3 生成流程

    算法的詳細(xì)說明如下。

    (1) 輸入表面網(wǎng)格R。

    (2) 支撐線的生成。按照一定的幾何特性選取部分邊界點作為支撐點,如圖1(a)選取A,B,C和D四個支撐點。按照網(wǎng)格點的法向方向和輸入的附面層高度h,算出支撐線AE,BF,CG和DH??梢酝ㄟ^平均支撐點周圍網(wǎng)格單元的法向得到每個支撐點的法向方向,

    (3)

    式中np為支撐點的法向方向,Nf為支撐點p相鄰網(wǎng)格單元總數(shù),ni為第i個網(wǎng)格單元的法向。

    該方法在個別點法向生成存在一定問題,這是附面層生成求法向所面臨的共同難題。

    (3) 輪廓線的生成。定義徑向基函數(shù)插值的被插值函數(shù)f(r)為選取基點的位移,r為表面網(wǎng)格邊界點。點A和點B的位移變形量為h,將點A和點B作為徑向基函數(shù)插值的基點,首先計算出插值系數(shù)αi,然后將線AB上的非支撐網(wǎng)格點進(jìn)行插值計算,得到線AB上非支撐網(wǎng)格點的位移變形量,從而得到EF上對應(yīng)的網(wǎng)格點。同理得到FG,F(xiàn)H和HE上對應(yīng)的網(wǎng)格點,生成輪廓線上的網(wǎng)格點Ro。

    對比文獻(xiàn)[15]輪廓線的生成是線性連接兩個相鄰的支撐線端點,這需要較好的網(wǎng)格邊界。而本文方法可以處理比較復(fù)雜的網(wǎng)格邊界的輪廓線生成,對網(wǎng)格拓?fù)湟筝^低。

    (4) 附面層頂層網(wǎng)格點的生成。求解插值系數(shù)αi時,徑向基函數(shù)插值中f(r)為表面邊界網(wǎng)格點的位移ΔR,r為表面網(wǎng)格邊界點,已知N個邊界點Rb為基點,邊界點的位移變形量為

    ΔR=Ro-Rb

    (4)

    利用徑向基函數(shù)插值可以得到求解插值系數(shù)的線性方程組

    ΔRx=Mαx, ΔRy=Mαy, ΔRz=Mαz

    (5~7)

    在求解表面網(wǎng)格內(nèi)部點的位移時,被插值函數(shù)f(r)為待求值,內(nèi)部網(wǎng)格的位移表示為ΔV。

    ΔV=Dα=DM-1ΔR=HΔR

    (8)

    式(8)寫成矩陣的形式為

    (9)

    求解矩陣方程(9)計算出表面網(wǎng)格內(nèi)部點的位移值ΔV,然后得到附面層頂層網(wǎng)格點Rt。把表面網(wǎng)格點的拓?fù)潢P(guān)系映射到附面層頂層網(wǎng)格點上,得到附面層頂層網(wǎng)格的拓?fù)潢P(guān)系,生成附面層頂層網(wǎng)格。

    (5 )附面層內(nèi)部網(wǎng)格點的生成。連接RtR,將其按照一定的增長率分成k份,得到內(nèi)部網(wǎng)格點Rp,通過表面網(wǎng)格拓?fù)潢P(guān)系和附面層頂層網(wǎng)格的鏈接關(guān)系,生成三棱柱附面層網(wǎng)格。

    2.4 算法測試

    圖2所示為法向推進(jìn)法和構(gòu)造法的二維算法。圖2(a)表明法向推進(jìn)法在局部因推進(jìn)過高會出現(xiàn)網(wǎng)格交叉的現(xiàn)象,一般通過降低附面層高度來避免,而圖2(b)通過框架線的約束,能有效解決此問題。

    圖2 法向推進(jìn)法和構(gòu)造法

    3 算例驗證

    通過4種典型幾何外形驗證本文提出的附面層網(wǎng)格生成方法,分別是V型槽、圓球、吊艙和F6機(jī)翼。測試在3.5 GHz主頻、16 GB內(nèi)存的計算機(jī)上進(jìn)行。

    3.1 V型槽

    V型槽是各種幾何形體中的一種常見結(jié)構(gòu)。如飛行器的機(jī)翼與機(jī)身連接處等。圖3為V型槽的附面層網(wǎng)格,圖中V型槽槽底附近的附面層網(wǎng)格不是沿著法向生長,而是在框架線和附面層頂層網(wǎng)格的約束下,沿著與法向有一定角度的方向生長,保證了附面層網(wǎng)格生成時不會出現(xiàn)交叉現(xiàn)象。

    V型槽表面網(wǎng)格共有1074個網(wǎng)格點。設(shè)置20層附面層層數(shù)、0.001 mm的首層高度和1.2的增長率。附面層的網(wǎng)格單元最小角為22.5°,最大角為135.3°,生成時間為0.4 s。

    圖3 V型槽的附面層網(wǎng)格

    3.2 圓球

    球面是最基本的幾何曲面。圖4所示為圓球的附面層網(wǎng)格,可以看出附面層網(wǎng)格的正交性較好,整體附面層網(wǎng)格質(zhì)量較好。

    圓球表面網(wǎng)格共有5024個網(wǎng)格點。設(shè)置20層附面層生長層數(shù)、0.001 mm的首層高度和1.2的增長率。附面層的網(wǎng)格單元最小角為25.4°,最大角為124.3°,生成時間為1.2 s。

    圖4 圓球的附面層網(wǎng)格

    3.3 吊艙

    吊艙是外凸面和內(nèi)凹面都存在的幾何形體。前緣是環(huán)形的凸曲面,后緣是個環(huán)形薄片。圖5(a)為吊艙表面網(wǎng)格和輪廓線。圖5(b,c)分別為吊艙的切片圖,圖5(d)為局部切片圖。從圖5(b~d)可以看出,采用本文方法生成的附面層網(wǎng)格,在吊艙特殊位置處,有效地避免了附面層網(wǎng)格生成過程中網(wǎng)格點法向量交叉的問題,提高了附面層網(wǎng)格質(zhì)量。

    圖5 吊艙的附面層網(wǎng)格

    吊艙表面網(wǎng)格共有19097個網(wǎng)格點。設(shè)置20層附面層生長層數(shù)、0.001 mm的首層高度和1.2的增長率。附面層的網(wǎng)格單元最小角為6.5°,最大角為171°,生成時間為21 s。

    3.4 F6機(jī)翼

    圖6(a)所示為F6機(jī)翼的表面網(wǎng)格和輪廓線。圖6(b)為F6機(jī)翼附面層的切片圖。圖6(c)為F6局部切片圖。從圖6(b,c)可以看出,通過本方法進(jìn)行附面層網(wǎng)格生成,在F6機(jī)翼的前后緣以及機(jī)翼邊界處,可以有效避免利用法向推進(jìn)法帶來的法向交叉問題。因為本方法選取少數(shù)網(wǎng)格點作為基點進(jìn)行法向生成,所以避免了交叉現(xiàn)象,較好地提高了網(wǎng)格質(zhì)量。

    圖6 F6機(jī)翼的附面層網(wǎng)格

    機(jī)翼表面網(wǎng)格共有12970個網(wǎng)格點。設(shè)置20層附面層生長層數(shù)、0.001 mm的首層高度和1.2的增長率。附面層的網(wǎng)格單元最小角為5.2°,最大角為161°,生成時間為31.4 s。

    4 結(jié) 論

    本文提出了一種基于框架的附面層棱柱網(wǎng)格生成方法,并通過幾種典型外形對該方法進(jìn)行了測試,成功生成了附面層網(wǎng)格,得到以下結(jié)論。通過構(gòu)建的框架線和附面層頂層網(wǎng)格點對附面層內(nèi)部網(wǎng)格點的約束,有效避免了附面層網(wǎng)格的交叉現(xiàn)象,保證了附面層網(wǎng)格的生成;表面網(wǎng)格發(fā)生變化時,僅需要對發(fā)生變化的局部區(qū)域重新生成附面層網(wǎng)格,即可快速重新生成整體的附面層網(wǎng)格,大幅降低了修改附面層網(wǎng)格的時間;大量減少了人工交互的工作,為進(jìn)一步附面層網(wǎng)格生成自動化提供了參考。

    本文在生成附面層內(nèi)部的網(wǎng)格點時,運用了簡單線性插值。表面網(wǎng)格形狀比較復(fù)雜的時候,會造成附面層網(wǎng)格正交性較差的現(xiàn)象。將復(fù)雜的表面網(wǎng)格進(jìn)行合理網(wǎng)格分區(qū),達(dá)到復(fù)雜表面網(wǎng)格簡單化的目的,或者選擇一種合適插值方法能有效地解決這個問題,下一步工作將進(jìn)一步嘗試和完善。

    猜你喜歡
    附面層邊界點棱柱
    道路空間特征與測量距離相結(jié)合的LiDAR道路邊界點提取算法
    層次化點云邊界快速精確提取方法研究
    基于數(shù)值模擬的流場附面層邊緣識別方法
    純位移線彈性方程Locking-Free非協(xié)調(diào)三棱柱單元的構(gòu)造分析
    立足概念,注重推理——以棱柱為例
    空間垂直關(guān)系錯解剖析
    超聲壓氣機(jī)葉柵流場的數(shù)值模擬與試驗驗證
    基于AT89C52單片機(jī)的三棱柱旋轉(zhuǎn)黑板的研究
    一種去除掛網(wǎng)圖像鋸齒的方法及裝置
    電腦與電信(2014年6期)2014-03-22 13:21:06
    空氣流量組合測量耙的研制及試驗
    日韩欧美 国产精品| 男女啪啪激烈高潮av片| 看黄色毛片网站| 99久久无色码亚洲精品果冻| 色综合亚洲欧美另类图片| 亚洲国产精品sss在线观看| 色哟哟哟哟哟哟| 国产白丝娇喘喷水9色精品| 淫妇啪啪啪对白视频| 成人鲁丝片一二三区免费| 一本久久中文字幕| 女人十人毛片免费观看3o分钟| 亚洲精品在线观看二区| 在线天堂最新版资源| 国产在线男女| 久久国产精品人妻蜜桃| 麻豆av噜噜一区二区三区| 天美传媒精品一区二区| 制服丝袜大香蕉在线| 久久久久久久午夜电影| 丰满人妻一区二区三区视频av| 国产一区二区三区av在线 | 身体一侧抽搐| 亚洲经典国产精华液单| 亚洲成av人片在线播放无| 日韩大尺度精品在线看网址| 成人亚洲精品av一区二区| aaaaa片日本免费| 国产亚洲精品久久久com| 欧美成人a在线观看| 高清日韩中文字幕在线| 久久精品久久久久久噜噜老黄 | 欧美激情久久久久久爽电影| 精品久久久久久久人妻蜜臀av| 最近最新中文字幕大全电影3| 99riav亚洲国产免费| 久久久久久久久久黄片| 高清日韩中文字幕在线| 色哟哟哟哟哟哟| 在线观看66精品国产| 日本撒尿小便嘘嘘汇集6| 99久久精品国产国产毛片| 国模一区二区三区四区视频| 噜噜噜噜噜久久久久久91| 在线免费观看的www视频| 一进一出好大好爽视频| 国产伦精品一区二区三区视频9| 国产午夜精品久久久久久一区二区三区 | 不卡视频在线观看欧美| 99久久精品一区二区三区| 精品久久久久久久久久久久久| 午夜精品久久久久久毛片777| 一区二区三区免费毛片| 九九热线精品视视频播放| 1000部很黄的大片| 小蜜桃在线观看免费完整版高清| 欧美高清性xxxxhd video| 国产精品久久久久久久久免| 国产成人av教育| 久久午夜亚洲精品久久| 性色avwww在线观看| 18禁在线播放成人免费| 亚洲av成人av| 美女高潮喷水抽搐中文字幕| 在线观看一区二区三区| 美女大奶头视频| 一级a爱片免费观看的视频| 日本免费一区二区三区高清不卡| 精品不卡国产一区二区三区| 亚洲av二区三区四区| 成年版毛片免费区| 亚洲专区国产一区二区| 国产毛片a区久久久久| 久久久久久久久久成人| 中文资源天堂在线| 亚洲午夜理论影院| 午夜福利高清视频| 联通29元200g的流量卡| 联通29元200g的流量卡| 免费观看精品视频网站| 国内精品久久久久精免费| 十八禁国产超污无遮挡网站| 一级av片app| 91麻豆av在线| 日本黄色片子视频| 欧美3d第一页| 亚洲一区高清亚洲精品| 在线看三级毛片| xxxwww97欧美| 免费看av在线观看网站| 久久草成人影院| 男人狂女人下面高潮的视频| 乱系列少妇在线播放| av视频在线观看入口| 国产精品精品国产色婷婷| 嫩草影院精品99| 国产一区二区三区在线臀色熟女| 精品久久久久久久久av| 深夜精品福利| 色哟哟哟哟哟哟| 又黄又爽又刺激的免费视频.| 高清毛片免费观看视频网站| 国产主播在线观看一区二区| 日韩国内少妇激情av| 欧美丝袜亚洲另类 | 最近最新免费中文字幕在线| 精品人妻视频免费看| 国产精品亚洲美女久久久| 欧美日韩亚洲国产一区二区在线观看| 亚洲熟妇中文字幕五十中出| h日本视频在线播放| 熟女电影av网| 午夜福利高清视频| 男女啪啪激烈高潮av片| 搡女人真爽免费视频火全软件 | 免费av观看视频| 国产精品99久久久久久久久| 国产精品一区二区三区四区免费观看 | 人妻少妇偷人精品九色| 非洲黑人性xxxx精品又粗又长| 日韩 亚洲 欧美在线| 久久精品国产亚洲网站| 中文字幕精品亚洲无线码一区| 97超级碰碰碰精品色视频在线观看| 国产成人av教育| 啦啦啦啦在线视频资源| 美女cb高潮喷水在线观看| 亚洲av免费高清在线观看| 麻豆国产av国片精品| 亚洲一区高清亚洲精品| 精品午夜福利视频在线观看一区| 国产精品电影一区二区三区| 国产在视频线在精品| 欧美潮喷喷水| 91狼人影院| 黄色日韩在线| 国内精品美女久久久久久| 午夜福利在线在线| 18禁黄网站禁片免费观看直播| 日韩中字成人| 88av欧美| 欧美日韩中文字幕国产精品一区二区三区| 偷拍熟女少妇极品色| 亚洲精品久久国产高清桃花| 国产精品嫩草影院av在线观看 | 国产又黄又爽又无遮挡在线| 九色国产91popny在线| 成人精品一区二区免费| 一个人看的www免费观看视频| 国产女主播在线喷水免费视频网站 | av中文乱码字幕在线| 男女下面进入的视频免费午夜| 99热只有精品国产| 欧美日韩中文字幕国产精品一区二区三区| 超碰av人人做人人爽久久| 日韩av在线大香蕉| 最后的刺客免费高清国语| 无人区码免费观看不卡| 校园春色视频在线观看| 久久久久久久亚洲中文字幕| 欧美高清成人免费视频www| 中文资源天堂在线| 国产毛片a区久久久久| 九九在线视频观看精品| 窝窝影院91人妻| 免费av观看视频| 又黄又爽又免费观看的视频| 精品一区二区三区视频在线观看免费| 亚洲第一区二区三区不卡| 2021天堂中文幕一二区在线观| 久久久久久伊人网av| 日韩强制内射视频| 成人特级av手机在线观看| 丝袜美腿在线中文| 成年人黄色毛片网站| 国产麻豆成人av免费视频| 夜夜夜夜夜久久久久| 中文字幕人妻熟人妻熟丝袜美| 18+在线观看网站| av福利片在线观看| 国产三级中文精品| 精品国产三级普通话版| 亚洲18禁久久av| 亚洲国产精品成人综合色| 久99久视频精品免费| 精品午夜福利视频在线观看一区| 联通29元200g的流量卡| 亚洲va日本ⅴa欧美va伊人久久| 一级黄色大片毛片| 国产精品久久久久久亚洲av鲁大| 97人妻精品一区二区三区麻豆| 亚洲av免费在线观看| av在线天堂中文字幕| www.www免费av| 亚洲欧美清纯卡通| aaaaa片日本免费| 别揉我奶头 嗯啊视频| 99精品在免费线老司机午夜| 欧美+日韩+精品| 麻豆久久精品国产亚洲av| 午夜福利在线在线| 精品一区二区免费观看| 日韩强制内射视频| 日韩av在线大香蕉| 热99re8久久精品国产| 欧洲精品卡2卡3卡4卡5卡区| 日日摸夜夜添夜夜添小说| 成人国产综合亚洲| 人人妻,人人澡人人爽秒播| 亚洲av.av天堂| 噜噜噜噜噜久久久久久91| 午夜福利在线观看吧| bbb黄色大片| 黄色丝袜av网址大全| 日韩,欧美,国产一区二区三区 | 久久久久久久久中文| 麻豆成人av在线观看| 亚洲欧美日韩无卡精品| 国产国拍精品亚洲av在线观看| 女人十人毛片免费观看3o分钟| 能在线免费观看的黄片| 欧美国产日韩亚洲一区| 波多野结衣巨乳人妻| 欧美三级亚洲精品| 午夜福利18| 变态另类丝袜制服| 色精品久久人妻99蜜桃| 免费看日本二区| 如何舔出高潮| 日韩一区二区视频免费看| 欧美最新免费一区二区三区| 精品人妻1区二区| 88av欧美| 日本免费a在线| 久久久色成人| 亚洲av不卡在线观看| 十八禁国产超污无遮挡网站| 白带黄色成豆腐渣| 日日啪夜夜撸| 麻豆国产97在线/欧美| 久久久精品欧美日韩精品| 春色校园在线视频观看| 女的被弄到高潮叫床怎么办 | 热99在线观看视频| 国产一区二区三区在线臀色熟女| 久久九九热精品免费| 免费看av在线观看网站| 天堂av国产一区二区熟女人妻| 久久精品91蜜桃| 男女下面进入的视频免费午夜| 嫩草影院精品99| 免费观看的影片在线观看| 日本三级黄在线观看| 国产精品不卡视频一区二区| 亚洲成人久久性| 丰满人妻一区二区三区视频av| 最近中文字幕高清免费大全6 | 亚洲天堂国产精品一区在线| 午夜日韩欧美国产| av.在线天堂| 免费av毛片视频| 欧美激情在线99| 久久精品国产亚洲网站| 午夜视频国产福利| 中文亚洲av片在线观看爽| 国内揄拍国产精品人妻在线| 18禁在线播放成人免费| 国产真实乱freesex| 一本精品99久久精品77| 老熟妇乱子伦视频在线观看| 99国产极品粉嫩在线观看| 在线免费观看的www视频| 日本免费a在线| 夜夜夜夜夜久久久久| 人妻丰满熟妇av一区二区三区| 99在线人妻在线中文字幕| 色综合站精品国产| 毛片女人毛片| 亚洲久久久久久中文字幕| 熟妇人妻久久中文字幕3abv| 成人国产一区最新在线观看| 变态另类丝袜制服| 国产亚洲精品久久久com| 婷婷丁香在线五月| 国产伦精品一区二区三区四那| 一区二区三区激情视频| 男人舔女人下体高潮全视频| 亚洲狠狠婷婷综合久久图片| 国产成人aa在线观看| 精品久久久噜噜| 精品一区二区三区人妻视频| 国产精品一及| 亚洲欧美清纯卡通| 国产久久久一区二区三区| 黄色配什么色好看| 无遮挡黄片免费观看| 亚洲av第一区精品v没综合| 特大巨黑吊av在线直播| 欧美高清成人免费视频www| 91久久精品电影网| 亚洲av成人av| 99久久精品国产国产毛片| 日韩欧美在线二视频| 99久久九九国产精品国产免费| www日本黄色视频网| 国内毛片毛片毛片毛片毛片| 别揉我奶头~嗯~啊~动态视频| 99热这里只有是精品50| 国产精品98久久久久久宅男小说| 婷婷色综合大香蕉| 最近视频中文字幕2019在线8| 99riav亚洲国产免费| 精品久久久久久成人av| 天天一区二区日本电影三级| 欧美日韩精品成人综合77777| 亚洲av熟女| 精品99又大又爽又粗少妇毛片 | 中文字幕熟女人妻在线| 很黄的视频免费| 亚洲成av人片在线播放无| 波多野结衣高清无吗| a级毛片a级免费在线| 亚洲欧美激情综合另类| 在线国产一区二区在线| 悠悠久久av| 亚洲国产欧洲综合997久久,| 欧美成人a在线观看| 精品日产1卡2卡| 国产亚洲精品av在线| 热99在线观看视频| 欧美xxxx性猛交bbbb| 麻豆成人av在线观看| 中文字幕熟女人妻在线| 亚洲在线自拍视频| 村上凉子中文字幕在线| 亚洲四区av| 看黄色毛片网站| 日韩精品青青久久久久久| 美女高潮喷水抽搐中文字幕| 欧美3d第一页| 亚洲中文字幕日韩| 国内毛片毛片毛片毛片毛片| 久久久久性生活片| 成人国产综合亚洲| 久久久久久伊人网av| 亚洲七黄色美女视频| 日日夜夜操网爽| 久久婷婷人人爽人人干人人爱| 亚洲人与动物交配视频| 一级黄色大片毛片| 日本免费a在线| 欧美潮喷喷水| 国产私拍福利视频在线观看| 午夜激情欧美在线| av福利片在线观看| 听说在线观看完整版免费高清| 精品人妻一区二区三区麻豆 | 九色国产91popny在线| 亚洲真实伦在线观看| 国产精华一区二区三区| 看免费成人av毛片| 丰满乱子伦码专区| 中文在线观看免费www的网站| 成人三级黄色视频| 国产欧美日韩一区二区精品| av在线亚洲专区| 亚州av有码| 国产v大片淫在线免费观看| 国产午夜福利久久久久久| 69人妻影院| 国产 一区 欧美 日韩| 日本与韩国留学比较| 亚洲国产色片| 婷婷精品国产亚洲av| 国国产精品蜜臀av免费| 亚洲天堂国产精品一区在线| 亚洲成人中文字幕在线播放| 精品一区二区三区视频在线| 网址你懂的国产日韩在线| 我要看日韩黄色一级片| 亚洲在线自拍视频| 国产 一区精品| 成年人黄色毛片网站| 国产伦在线观看视频一区| 国产精品久久久久久精品电影| 人妻少妇偷人精品九色| 人妻丰满熟妇av一区二区三区| 一级a爱片免费观看的视频| 成年人黄色毛片网站| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲自偷自拍三级| 成年女人毛片免费观看观看9| av在线蜜桃| 亚洲aⅴ乱码一区二区在线播放| 国产欧美日韩精品一区二区| 久久久成人免费电影| 啦啦啦啦在线视频资源| 永久网站在线| 高清毛片免费观看视频网站| 亚洲色图av天堂| 毛片一级片免费看久久久久 | 国产一区二区激情短视频| 成人欧美大片| 国产一区二区在线av高清观看| 可以在线观看毛片的网站| 国产一区二区激情短视频| 国产在线男女| 噜噜噜噜噜久久久久久91| 国产精品不卡视频一区二区| 国产精品女同一区二区软件 | 最后的刺客免费高清国语| 国产亚洲精品久久久久久毛片| 亚洲精品久久国产高清桃花| 全区人妻精品视频| 亚洲最大成人手机在线| 麻豆久久精品国产亚洲av| 国产真实乱freesex| 精品人妻偷拍中文字幕| 99热6这里只有精品| 最好的美女福利视频网| 免费看日本二区| 免费观看精品视频网站| 免费不卡的大黄色大毛片视频在线观看 | 一级黄片播放器| 亚洲精品亚洲一区二区| xxxwww97欧美| 久久国产乱子免费精品| 精品99又大又爽又粗少妇毛片 | 免费人成在线观看视频色| 有码 亚洲区| 人人妻人人看人人澡| 听说在线观看完整版免费高清| 伦精品一区二区三区| av视频在线观看入口| 波野结衣二区三区在线| 俺也久久电影网| 精品一区二区三区人妻视频| 日韩 亚洲 欧美在线| 国产男靠女视频免费网站| 中出人妻视频一区二区| 久久这里只有精品中国| 亚洲精品乱码久久久v下载方式| 国产精品久久久久久精品电影| 国产色爽女视频免费观看| 十八禁网站免费在线| 看片在线看免费视频| 噜噜噜噜噜久久久久久91| 精品久久久噜噜| 91av网一区二区| 欧洲精品卡2卡3卡4卡5卡区| 欧美黑人欧美精品刺激| 欧美高清成人免费视频www| 我要搜黄色片| 亚洲av二区三区四区| 一本久久中文字幕| 男女啪啪激烈高潮av片| 中出人妻视频一区二区| 亚洲av.av天堂| 联通29元200g的流量卡| 免费大片18禁| 色综合色国产| 国内精品宾馆在线| 日韩欧美国产一区二区入口| 色av中文字幕| 亚洲精品色激情综合| 精品久久久久久久人妻蜜臀av| 日韩欧美免费精品| 欧美最新免费一区二区三区| 在线免费十八禁| 人人妻人人澡欧美一区二区| 亚洲不卡免费看| 五月玫瑰六月丁香| 中国美白少妇内射xxxbb| 亚洲欧美日韩卡通动漫| 免费观看在线日韩| aaaaa片日本免费| 欧美bdsm另类| 成人美女网站在线观看视频| 噜噜噜噜噜久久久久久91| 欧美日韩瑟瑟在线播放| 亚洲欧美清纯卡通| 男人舔奶头视频| 欧美成人免费av一区二区三区| 亚洲av.av天堂| 免费一级毛片在线播放高清视频| 少妇裸体淫交视频免费看高清| 婷婷精品国产亚洲av| 久久久色成人| 无遮挡黄片免费观看| 亚洲男人的天堂狠狠| 美女高潮喷水抽搐中文字幕| 在线播放无遮挡| 亚洲黑人精品在线| 悠悠久久av| 亚洲一区高清亚洲精品| 日日撸夜夜添| 人人妻人人看人人澡| 欧美zozozo另类| 永久网站在线| 啦啦啦啦在线视频资源| 国产91精品成人一区二区三区| a级一级毛片免费在线观看| 亚洲av一区综合| 欧美日韩综合久久久久久 | 欧美区成人在线视频| 日韩 亚洲 欧美在线| 亚洲精品成人久久久久久| 久久国产精品人妻蜜桃| 听说在线观看完整版免费高清| 国产乱人伦免费视频| 久久久国产成人免费| 午夜福利18| 天堂影院成人在线观看| or卡值多少钱| av天堂中文字幕网| 热99re8久久精品国产| 少妇丰满av| 国产色婷婷99| 日本爱情动作片www.在线观看 | 免费电影在线观看免费观看| 亚洲午夜理论影院| 国模一区二区三区四区视频| 此物有八面人人有两片| 欧美国产日韩亚洲一区| 成人特级av手机在线观看| 人妻制服诱惑在线中文字幕| 日本欧美国产在线视频| 精品人妻视频免费看| 最近最新中文字幕大全电影3| 亚洲最大成人av| 国产精品1区2区在线观看.| 国产精品一区www在线观看 | 天堂网av新在线| 亚洲av一区综合| 日韩中文字幕欧美一区二区| 麻豆av噜噜一区二区三区| 高清日韩中文字幕在线| 中国美白少妇内射xxxbb| 日韩欧美精品免费久久| av中文乱码字幕在线| av黄色大香蕉| 久久久精品欧美日韩精品| 三级男女做爰猛烈吃奶摸视频| 热99re8久久精品国产| 草草在线视频免费看| 美女被艹到高潮喷水动态| 两个人视频免费观看高清| 国产免费一级a男人的天堂| 久久草成人影院| 十八禁网站免费在线| 国产精品乱码一区二三区的特点| 99在线视频只有这里精品首页| 婷婷色综合大香蕉| 窝窝影院91人妻| 少妇猛男粗大的猛烈进出视频 | 日韩精品中文字幕看吧| 成人性生交大片免费视频hd| 精品国产三级普通话版| 赤兔流量卡办理| 免费搜索国产男女视频| 日韩人妻高清精品专区| 亚洲国产精品久久男人天堂| 露出奶头的视频| 国产精品亚洲一级av第二区| xxxwww97欧美| 免费不卡的大黄色大毛片视频在线观看 | 久久久久免费精品人妻一区二区| 久久久久国内视频| 村上凉子中文字幕在线| 淫妇啪啪啪对白视频| 欧美成人性av电影在线观看| 欧美一区二区精品小视频在线| 免费黄网站久久成人精品| 丰满乱子伦码专区| 女同久久另类99精品国产91| 两人在一起打扑克的视频| 免费不卡的大黄色大毛片视频在线观看 | 搞女人的毛片| 日韩亚洲欧美综合| 久久精品久久久久久噜噜老黄 | 天美传媒精品一区二区| 两个人的视频大全免费| 亚洲久久久久久中文字幕| 淫秽高清视频在线观看| av福利片在线观看| 夜夜夜夜夜久久久久| 欧美黑人欧美精品刺激| 国内精品一区二区在线观看| 午夜爱爱视频在线播放| 国产免费av片在线观看野外av| 韩国av一区二区三区四区| 毛片女人毛片| 欧美bdsm另类| 久久久久久久久中文| 一本久久中文字幕| 18禁黄网站禁片免费观看直播| 国产激情偷乱视频一区二区| 日韩欧美精品免费久久| 91午夜精品亚洲一区二区三区 | 91av网一区二区| av黄色大香蕉| 亚洲国产欧美人成| 不卡视频在线观看欧美| 性色avwww在线观看| 中文字幕av成人在线电影| 中出人妻视频一区二区| 国产私拍福利视频在线观看| 国产激情偷乱视频一区二区| 日本爱情动作片www.在线观看 | 成人毛片a级毛片在线播放| 香蕉av资源在线| 免费一级毛片在线播放高清视频| 校园人妻丝袜中文字幕| 99视频精品全部免费 在线| 一进一出抽搐动态|