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

    基于Kriging 模型和網(wǎng)格變形的翼型優(yōu)化

    2022-04-20 08:20:40胡勇軍高一星
    科學(xué)技術(shù)創(chuàng)新 2022年11期
    關(guān)鍵詞:控制點(diǎn)氣動(dòng)網(wǎng)格

    胡勇軍 高一星

    (貴州貴飛飛機(jī)設(shè)計(jì)研究院有限公司,貴州 安順 561000)

    機(jī)翼通常是飛行器產(chǎn)生升力和阻力的主要部件,而在機(jī)翼設(shè)計(jì)中作為基本要素之一的翼型對飛行器氣動(dòng)特性和飛行性能都有很大的影響,因此選用擁有突出的氣動(dòng)特性的翼型對飛行器氣動(dòng)設(shè)計(jì)非常重要。早期,傳統(tǒng)的翼型設(shè)計(jì)方式都是首先取得已有翼型的試驗(yàn)資料,然后再結(jié)合翼型理論和一定的設(shè)計(jì)經(jīng)驗(yàn)等,加以剪裁和修形,最后才通過風(fēng)洞試驗(yàn)進(jìn)行驗(yàn)證。該設(shè)計(jì)過程往往要求工程設(shè)計(jì)人員必須具備較高水平的專業(yè)技術(shù)知識(shí)和較豐富的工程設(shè)計(jì)經(jīng)驗(yàn),隨著現(xiàn)代飛機(jī)設(shè)計(jì)中設(shè)計(jì)周期的急劇降低,這種方法使用的越來越少。二十世紀(jì)七十年代起,由于計(jì)算機(jī)技術(shù)的快速發(fā)展,計(jì)算流體力學(xué)(CFD)同時(shí)獲得了飛速發(fā)展,并成為空氣動(dòng)力學(xué)中發(fā)展最快的方向之一,也成為了飛行器設(shè)計(jì)過程中必不可少的技術(shù)手段。CFD 為飛機(jī)翼型設(shè)計(jì)帶來了全新的技術(shù)手段,將原由風(fēng)洞完成的大量試驗(yàn)變?yōu)榱嘶贑FD 的大量仿真計(jì)算再加少量風(fēng)洞試驗(yàn)驗(yàn)證,從而大大縮短了周期,同時(shí)提高了獲得具有優(yōu)秀氣動(dòng)特性的翼型的可能性。

    在基于CFD 的翼型設(shè)計(jì)方法中,數(shù)值優(yōu)化設(shè)計(jì)方法是常用的方法之一[1-2]。該方法可以選取目標(biāo)翼型的指定迎角的升力系數(shù)、升阻比、阻力、俯仰力矩等氣動(dòng)特性作為設(shè)計(jì)目標(biāo),疊加諸如相對厚度分布等設(shè)計(jì)約束,通過遺傳算法等方法,快速獲得滿足要求的翼型。國內(nèi)外均有許多學(xué)者進(jìn)行了該方向的研究。氣動(dòng)外形優(yōu)化設(shè)計(jì)的基本過程就是不斷調(diào)用參數(shù)計(jì)算模型進(jìn)行計(jì)算的過程,如果該過程使用高精度模型用于計(jì)算,計(jì)算過程將會(huì)變得非常復(fù)雜,為降低優(yōu)化過程計(jì)算的復(fù)雜程度,使用的模型最好使用近似方法建立,即用一個(gè)簡單近似函數(shù)代替高精度的復(fù)雜的計(jì)算模型,比如采用局部近似、全局近似等。在全局近似法中,Kriging 模型因其所具有的非線性擬合能力強(qiáng)、近似模型參數(shù)選取靈活性高等優(yōu)勢,在翼型優(yōu)化等氣動(dòng)優(yōu)化設(shè)計(jì)實(shí)踐中得到了廣泛應(yīng)用,如Meunier 就通過該方法建立了機(jī)翼氣動(dòng)力模型[3]等。

    本文在基準(zhǔn)翼型上,采用網(wǎng)格變形和CFD 方法獲得試驗(yàn)樣本點(diǎn)數(shù)據(jù),通過Kriging 模型優(yōu)化設(shè)計(jì)[4-5]出具有要求的目標(biāo)特性的翼型,采用該翼型設(shè)計(jì)的機(jī)翼隨全機(jī)一起進(jìn)行風(fēng)洞試驗(yàn),飛機(jī)氣動(dòng)特性滿足全機(jī)設(shè)計(jì)要求,成功驗(yàn)證了翼型設(shè)計(jì)。

    1 Kriging 模型

    響應(yīng)面法的基本思想是假設(shè)一個(gè)包含一些未知參數(shù)的相對簡單的顯式函數(shù)表達(dá)式,通過試驗(yàn)設(shè)計(jì)對給定的設(shè)計(jì)集合進(jìn)行試驗(yàn),然后用插值等方式確定表達(dá)式中的未知參數(shù),得到顯式的響應(yīng)面模型以逼近目標(biāo)函數(shù),用來估計(jì)非試驗(yàn)點(diǎn)的目標(biāo)響應(yīng)值。Kriging 模型方法屬于一種響應(yīng)面法,在1951 年由南非的Krige 首次提出,該方法由于有對試驗(yàn)點(diǎn)無偏估計(jì)、非線性逼近能力較好等優(yōu)點(diǎn),在高度非線性的優(yōu)化問題中得以大量應(yīng)用。該方法主要利用對未知點(diǎn)附近的信息的加權(quán)線性組合來估計(jì)該點(diǎn)的值,其本質(zhì)就是通過將估計(jì)值的誤差方差最小化來確定加權(quán)系數(shù)。

    Kriging 模型[6]由線性回歸部分(多項(xiàng)式部分)和非參數(shù)部分(多項(xiàng)式的偏離部分)組成。假設(shè)在一個(gè)d 維設(shè)計(jì)空間中,有n 個(gè)樣本點(diǎn)的集合X=[x1,x2,…,xn]T,單個(gè)樣本點(diǎn)為xi=[x,x…,x],樣本點(diǎn)處的響應(yīng)值為Y=[y1,y2,…,yn]T,則采用kriging 函數(shù)構(gòu)造如下形式近似模型。

    其中,f (X,β) 為用線性組合的多項(xiàng)式表示的回歸函數(shù),是線性回歸部分,β 為回歸模型系數(shù),整個(gè)部分也可理解為在設(shè)計(jì)空間中建立的全局近似,該多項(xiàng)式的階次通常可取零階、一階、二階等,一般不會(huì)取高于二階,而據(jù)尹大偉等研究[7],通?;貧w模型的階次對精度的影響并不是很明顯;z(X)是均值為零的非參數(shù)模型,屬于隨機(jī)函數(shù)部分,它也可理解為平均數(shù)為零的隨機(jī)的服從正態(tài)分布但協(xié)方差不為0 的局部逼近。隨機(jī)函數(shù)z(X)代表了真實(shí)響應(yīng)與假設(shè)的全局逼近的局部偏差,其統(tǒng)計(jì)學(xué)特性及樣本點(diǎn)之間協(xié)方差為:

    其中,σ2為隨機(jī)過程方差;R 是相關(guān)函數(shù)矩陣;r 則是與二個(gè)樣本點(diǎn)xi,xj之間的距離相關(guān)的空間相關(guān)函數(shù);θ為空間相關(guān)函數(shù)的比例參數(shù)向量,代表了兩個(gè)樣本點(diǎn)之間的相關(guān)性隨兩點(diǎn)距離增大而衰減的程度??臻g相關(guān)函數(shù)可以是高斯函數(shù)、指數(shù)函數(shù)、樣條函數(shù)等。在使用高斯函數(shù)時(shí),輸入向量xi、xj的第k 個(gè)分量的空間相關(guān)函數(shù)為:

    其中,θk是組成向量θ 中的第k 個(gè)量,dk為:

    相關(guān)函數(shù)的選取既具體決定了代理模型中如何擬合數(shù)據(jù),也反映了與全局逼近的局部偏差,一般選取相對比較簡單的相關(guān)函數(shù)。在本文中,選擇的是ANSYS WORKBENCH 自帶的高斯函數(shù)。

    2 網(wǎng)格變形

    在翼型優(yōu)化中,對翼型的處理和描述有很多種方法,比如采用CST 參數(shù)化描述翼型[8],采用網(wǎng)格變形獲得翼型等。采用CST 參數(shù)化方法描述翼型時(shí),參數(shù)化表示方法直接影響翼型設(shè)計(jì)方法的合理性,進(jìn)而使得翼型表示方法成為采用該方法設(shè)計(jì)翼型的重要領(lǐng)域之一。而網(wǎng)格變形是通過對網(wǎng)格模型的形狀進(jìn)行修改的幾何處理技術(shù),網(wǎng)格變形技術(shù)避免了反復(fù)進(jìn)行CAD 設(shè)計(jì)更改和再生成網(wǎng)格的環(huán)節(jié),是近年來幾何處理領(lǐng)域的熱點(diǎn)研究問題之一。用于網(wǎng)格變形的方法通常有自由變形技術(shù)、骨架驅(qū)動(dòng)變形技術(shù)、基于網(wǎng)格曲面變形技術(shù)、基于徑向基函數(shù)的網(wǎng)格變形技術(shù)等。

    本文在基準(zhǔn)翼型網(wǎng)格基礎(chǔ)上,使用ANSYS 自帶的網(wǎng)格快速變形技術(shù),建立局部網(wǎng)格變形區(qū)域的拓?fù)淇刂茀?shù),通過參數(shù)值變化控制網(wǎng)格的隨移變形,改變網(wǎng)格節(jié)點(diǎn)位置,實(shí)現(xiàn)基于純網(wǎng)格的設(shè)計(jì)變形,獲得不同的氣動(dòng)外形網(wǎng)格。

    3 流動(dòng)控制方程及邊界條件

    本文CFD 計(jì)算使用的控制方程組為基于Reynolds時(shí)均法的方程組,湍流模型為SST K-ω 模型,邊界條件為壓力遠(yuǎn)場邊界條件。方程組的求解方法則選擇了基于壓力求解的方法,通過耦合隱式方法進(jìn)行有限體積離散,而梯度離散方法采用基于網(wǎng)格點(diǎn)的Green-Gauss 方法。本文采用NASA-LANGLEY LS(1)-0413 作為初始翼型,進(jìn)行多目標(biāo)和約束的翼型優(yōu)化設(shè)計(jì)。

    4 翼型氣動(dòng)優(yōu)化設(shè)計(jì)

    4.1 設(shè)計(jì)變量選取

    初始翼型網(wǎng)格見圖1。圖中同時(shí)給出了翼型上部網(wǎng)格的變形區(qū)域,變形區(qū)域從翼型前緣到后緣,可用于上部網(wǎng)格變形的7 個(gè)控制點(diǎn)均勻分布在弦線上且不含前后緣點(diǎn),變形方向?yàn)樯舷伦冃?,從前往后分別為控制點(diǎn)1至控制點(diǎn)7。同理,用于翼型下部網(wǎng)格變形的7 個(gè)控制點(diǎn)也均勻分布在弦線上且不含前后緣點(diǎn),變形方向也為上下方向,從前往后分別為控制點(diǎn)8 至控制點(diǎn)14。控制點(diǎn)合計(jì)14 個(gè),控制點(diǎn)變形參數(shù)作為輸入變量。為了分別控制翼型上表面和下表面變形,雖然控制點(diǎn)分布相同,但區(qū)分用于上下表面控制的控制點(diǎn)。網(wǎng)格變形過程中前后緣點(diǎn)位置保持不變。為避免極端情況下出現(xiàn)負(fù)體積和上下表面交錯(cuò)的情況,各允許變形控制點(diǎn)的變形參數(shù)范圍需在優(yōu)化前通過網(wǎng)格變形測試獲得,控制點(diǎn)參數(shù)設(shè)置時(shí)作為限制條件??刂泣c(diǎn)1 至控制點(diǎn)12 的輸入?yún)?shù)均為0.05、控制點(diǎn)13 輸入?yún)?shù)為0.03、控制點(diǎn)14 輸入?yún)?shù)為0 時(shí),網(wǎng)格變形結(jié)果見圖2。在輸出變量中,包括設(shè)計(jì)目標(biāo)所需的2 個(gè)迎角的升力系數(shù)、阻力系數(shù)。此外,為了控制厚度分布,控制點(diǎn)2 與控制點(diǎn)9 的差量,控制點(diǎn)3 與控制點(diǎn)10 的差量,控制點(diǎn)4 與控制點(diǎn)11 的差量,這3 個(gè)量也作為輸出變量。輸出變量合計(jì)7 個(gè)。

    圖1 基準(zhǔn)翼型網(wǎng)格

    圖2 某輸入?yún)?shù)下網(wǎng)格變形結(jié)果

    4.2 目標(biāo)和設(shè)計(jì)約束

    目標(biāo)翼型在Ma=0.4 巡航狀態(tài)需具有較高的升阻比、較小的阻力以具有較長的航時(shí),在Ma=0.15 大迎角狀態(tài)還需具有較高的升力系數(shù)。由于飛機(jī)目標(biāo)為按照機(jī)翼安裝角-3 度設(shè)計(jì)并在0 度迎角附近巡航,翼型設(shè)計(jì)目標(biāo)為:(1)Ma=0.4,在迎角-3 度升力系數(shù)不低于0.15 且盡量大,阻力系數(shù)不高于0.0053;(2)Ma=0.15,在迎角10 度時(shí)升力系數(shù)不低于1.7 且盡量大。設(shè)計(jì)約束條件主要是厚度條件,相對厚度可以減小但不能低于12%。

    4.3 樣本分布生成

    在建立Kriging 模型前,需要完成試驗(yàn)設(shè)計(jì)(DOE)以獲得樣本點(diǎn)。試驗(yàn)設(shè)計(jì)中,通過取樣策略確定樣本點(diǎn)的數(shù)量和樣本空間分布。常用的樣本設(shè)計(jì)方法一般有正交設(shè)計(jì)、中心復(fù)合設(shè)計(jì)、均勻設(shè)計(jì)、隨機(jī)取點(diǎn)設(shè)計(jì)、拉丁超立方設(shè)計(jì)等方法,也可以自己創(chuàng)建樣本設(shè)計(jì)方法。抽取的樣本點(diǎn)希望覆蓋比較均勻,又希望顯著減少樣本數(shù)量以降低計(jì)算規(guī)模和計(jì)算周期。本文中,采用的是ANSYS WORKBENCH 中的Custom+Sampling 方法先建立總數(shù)為300 的初始試驗(yàn)樣本,在必要時(shí)將后續(xù)優(yōu)化驗(yàn)證結(jié)果也添加到試驗(yàn)樣本中以有效填充設(shè)計(jì)空間。

    4.4 優(yōu)化設(shè)計(jì)結(jié)果與分析

    在ANSYS WORKBENCH 中采用上述方法建立優(yōu)化流程,并采用CFD 方法獲得試驗(yàn)樣本作為訓(xùn)練樣本。在獲取了各樣本點(diǎn)的具體數(shù)據(jù)之后,采用Kriging 模型構(gòu)建通過翼型設(shè)計(jì)參數(shù)計(jì)算氣動(dòng)特性的近似模型,在上述約束條件下采用多目標(biāo)優(yōu)化方法進(jìn)行翼型優(yōu)化,不進(jìn)行目標(biāo)函數(shù)歸一化處理以便直接獲得設(shè)計(jì)點(diǎn)處的設(shè)計(jì)值。優(yōu)化設(shè)計(jì)前后翼型的氣動(dòng)特性結(jié)果對比見表1,優(yōu)化設(shè)計(jì)前后的翼型形狀對比見圖3。優(yōu)化后,翼型相對厚度從13%減小到12.5%。從表1 可見,經(jīng)過優(yōu)化設(shè)計(jì)后,在Ma=0.4 時(shí),阻力系數(shù)增加較小,但升力系數(shù)顯著提高,升阻比提高較多;低速大迎角時(shí)的升力系數(shù)也有所增加。

    表1 優(yōu)化設(shè)計(jì)前后對比

    圖3 優(yōu)化前后翼型

    5 全機(jī)試驗(yàn)驗(yàn)證

    采用上述優(yōu)化后的翼型生成機(jī)翼,裝配于機(jī)身,按照1:7.5 的比例縮比設(shè)計(jì)和制造了風(fēng)洞試驗(yàn)?zāi)P停贔L-8風(fēng)洞中完成了全機(jī)風(fēng)洞試驗(yàn)。全機(jī)試驗(yàn)?zāi)P鸵妶D4,全機(jī)典型狀態(tài)試驗(yàn)結(jié)果見圖5、圖6。采用優(yōu)化后的翼型設(shè)計(jì)生成的機(jī)翼及全機(jī),在典型狀態(tài)的氣動(dòng)力試驗(yàn)結(jié)果滿足了設(shè)計(jì)要求,達(dá)到了預(yù)期目的。

    圖4 試驗(yàn)?zāi)P蛨D

    圖5 風(fēng)洞試驗(yàn)結(jié)果(升力系數(shù)~阻力系數(shù))

    圖6 風(fēng)洞試驗(yàn)結(jié)果(升力系數(shù)~迎角)

    6 結(jié)論

    本文通過使用網(wǎng)格變形技術(shù)和CFD 方法獲得翼型優(yōu)化的試驗(yàn)樣本,采用Kriging 模型進(jìn)行了多目標(biāo)多約束條件下的低速翼型優(yōu)化設(shè)計(jì)研究,獲得了滿足設(shè)計(jì)目標(biāo)的翼型。將優(yōu)化獲得的翼型隨全機(jī)進(jìn)行了風(fēng)洞試驗(yàn)驗(yàn)證,全機(jī)氣動(dòng)特性滿足設(shè)計(jì)目標(biāo)要求。上述方法設(shè)計(jì)的翼型通過全機(jī)試驗(yàn)驗(yàn)證,表明該翼型優(yōu)化設(shè)計(jì)方法可以快速有效地用于飛機(jī)翼型設(shè)計(jì),極大地推動(dòng)了目標(biāo)飛機(jī)進(jìn)一步的氣動(dòng)設(shè)計(jì)工作。

    猜你喜歡
    控制點(diǎn)氣動(dòng)網(wǎng)格
    用全等三角形破解網(wǎng)格題
    中寰氣動(dòng)執(zhí)行機(jī)構(gòu)
    基于NACA0030的波紋狀翼型氣動(dòng)特性探索
    反射的橢圓隨機(jī)偏微分方程的網(wǎng)格逼近
    基于反饋線性化的RLV氣動(dòng)控制一體化設(shè)計(jì)
    重疊網(wǎng)格裝配中的一種改進(jìn)ADT搜索方法
    NFFD控制點(diǎn)分布對氣動(dòng)外形優(yōu)化的影響
    基于風(fēng)險(xiǎn)管理下的項(xiàng)目建設(shè)內(nèi)部控制點(diǎn)思考
    基于曲面展開的自由曲面網(wǎng)格劃分
    相似材料模型中控制點(diǎn)像點(diǎn)坐標(biāo)定位研究
    99国产精品免费福利视频| 国内毛片毛片毛片毛片毛片| www.自偷自拍.com| 激情视频va一区二区三区| 久久婷婷成人综合色麻豆| 18禁观看日本| 涩涩av久久男人的天堂| 在线视频色国产色| 精品国产一区二区三区四区第35| 国产在线观看jvid| 一a级毛片在线观看| 久久精品国产亚洲av香蕉五月| 99国产综合亚洲精品| 麻豆国产av国片精品| 丰满人妻熟妇乱又伦精品不卡| 亚洲自偷自拍图片 自拍| 好男人电影高清在线观看| 亚洲成av片中文字幕在线观看| 色精品久久人妻99蜜桃| 成人18禁高潮啪啪吃奶动态图| 手机成人av网站| 日韩成人在线观看一区二区三区| 天堂中文最新版在线下载| 精品久久久久久电影网| 久久精品国产综合久久久| 国产精品综合久久久久久久免费 | 三上悠亚av全集在线观看| 宅男免费午夜| 一边摸一边抽搐一进一小说| 亚洲国产中文字幕在线视频| 又紧又爽又黄一区二区| 欧美激情 高清一区二区三区| 午夜福利在线观看吧| 校园春色视频在线观看| 亚洲专区国产一区二区| 18禁国产床啪视频网站| 两性夫妻黄色片| 中文字幕最新亚洲高清| 一区二区三区激情视频| 长腿黑丝高跟| 人人妻人人澡人人看| xxxhd国产人妻xxx| 97碰自拍视频| 99国产综合亚洲精品| 老司机午夜十八禁免费视频| 国产免费男女视频| 国产xxxxx性猛交| 国内毛片毛片毛片毛片毛片| 高清av免费在线| 欧美激情久久久久久爽电影 | 国产成人精品在线电影| 欧美精品一区二区免费开放| 日日摸夜夜添夜夜添小说| 妹子高潮喷水视频| 女性生殖器流出的白浆| e午夜精品久久久久久久| netflix在线观看网站| 好男人电影高清在线观看| 国产一区二区在线av高清观看| 欧美色视频一区免费| 淫妇啪啪啪对白视频| 成人特级黄色片久久久久久久| 精品国产美女av久久久久小说| 最近最新中文字幕大全免费视频| 国产精品 欧美亚洲| 国产精品久久久久成人av| 91成人精品电影| 少妇裸体淫交视频免费看高清 | 69精品国产乱码久久久| av国产精品久久久久影院| 国产亚洲精品久久久久久毛片| 午夜a级毛片| 91av网站免费观看| 狂野欧美激情性xxxx| 色播在线永久视频| 午夜免费成人在线视频| 美女高潮到喷水免费观看| 成人永久免费在线观看视频| 欧美乱码精品一区二区三区| 国产亚洲av高清不卡| 精品国产超薄肉色丝袜足j| 999精品在线视频| 亚洲国产欧美日韩在线播放| 黄网站色视频无遮挡免费观看| 国产无遮挡羞羞视频在线观看| www.自偷自拍.com| 天天影视国产精品| 80岁老熟妇乱子伦牲交| 老汉色av国产亚洲站长工具| 嫩草影院精品99| 亚洲精品中文字幕一二三四区| 久久天躁狠狠躁夜夜2o2o| 亚洲黑人精品在线| 亚洲成国产人片在线观看| av国产精品久久久久影院| 麻豆久久精品国产亚洲av | 亚洲欧美激情综合另类| 一本大道久久a久久精品| 涩涩av久久男人的天堂| 大型av网站在线播放| 成年女人毛片免费观看观看9| 国产成+人综合+亚洲专区| 亚洲成人精品中文字幕电影 | 在线观看免费高清a一片| 伊人久久大香线蕉亚洲五| 天天添夜夜摸| 男女高潮啪啪啪动态图| 动漫黄色视频在线观看| 久久久久久大精品| 亚洲国产中文字幕在线视频| 每晚都被弄得嗷嗷叫到高潮| 久久人妻福利社区极品人妻图片| 日本三级黄在线观看| 亚洲欧美激情在线| 淫秽高清视频在线观看| 高清av免费在线| 在线观看免费视频日本深夜| 欧美日韩中文字幕国产精品一区二区三区 | 男女床上黄色一级片免费看| 亚洲一卡2卡3卡4卡5卡精品中文| netflix在线观看网站| 波多野结衣av一区二区av| 亚洲一码二码三码区别大吗| 国产精品一区二区免费欧美| 亚洲成人免费av在线播放| 满18在线观看网站| 我的亚洲天堂| 无遮挡黄片免费观看| 在线看a的网站| 久久伊人香网站| 中文亚洲av片在线观看爽| 十八禁网站免费在线| 无人区码免费观看不卡| 老熟妇乱子伦视频在线观看| 欧美老熟妇乱子伦牲交| 女人精品久久久久毛片| 亚洲人成网站在线播放欧美日韩| 午夜福利在线免费观看网站| 99在线视频只有这里精品首页| 国产片内射在线| 老司机深夜福利视频在线观看| 日本黄色视频三级网站网址| 亚洲性夜色夜夜综合| 免费在线观看亚洲国产| 亚洲五月色婷婷综合| 亚洲av成人av| 国产亚洲欧美精品永久| 三上悠亚av全集在线观看| 另类亚洲欧美激情| 国产亚洲欧美在线一区二区| 国产深夜福利视频在线观看| 亚洲精品久久成人aⅴ小说| 日韩人妻精品一区2区三区| 久久久国产一区二区| 好看av亚洲va欧美ⅴa在| 亚洲全国av大片| 日韩视频一区二区在线观看| 在线观看免费日韩欧美大片| 亚洲一区高清亚洲精品| 亚洲色图综合在线观看| 亚洲三区欧美一区| 精品一区二区三卡| 亚洲av第一区精品v没综合| 色婷婷av一区二区三区视频| 成人精品一区二区免费| 老司机靠b影院| 日韩欧美三级三区| 免费看十八禁软件| √禁漫天堂资源中文www| 中文字幕精品免费在线观看视频| 99国产极品粉嫩在线观看| 亚洲精品美女久久av网站| 操出白浆在线播放| 国产精品免费视频内射| 久久人妻福利社区极品人妻图片| 成人影院久久| 天天影视国产精品| 99久久国产精品久久久| 侵犯人妻中文字幕一二三四区| 黄色视频不卡| 老司机福利观看| 亚洲精品中文字幕在线视频| cao死你这个sao货| 韩国av一区二区三区四区| 国产免费av片在线观看野外av| 国产欧美日韩综合在线一区二区| 亚洲人成网站在线播放欧美日韩| 日本a在线网址| 久久午夜综合久久蜜桃| 亚洲国产精品999在线| 久久天躁狠狠躁夜夜2o2o| 麻豆一二三区av精品| 久久亚洲精品不卡| 午夜老司机福利片| 国产99白浆流出| 久久中文字幕人妻熟女| 脱女人内裤的视频| 一进一出好大好爽视频| 十分钟在线观看高清视频www| 亚洲三区欧美一区| 欧美在线一区亚洲| 精品久久久精品久久久| 国产黄a三级三级三级人| 可以免费在线观看a视频的电影网站| 欧美成人午夜精品| 国产成人精品无人区| 国产99久久九九免费精品| 无遮挡黄片免费观看| 国产一区二区三区在线臀色熟女 | 波多野结衣av一区二区av| 大香蕉久久成人网| 国产乱人伦免费视频| 色综合站精品国产| 午夜久久久在线观看| 男人舔女人的私密视频| 18禁裸乳无遮挡免费网站照片 | 国产区一区二久久| 一边摸一边做爽爽视频免费| 黄频高清免费视频| 看黄色毛片网站| 国产精品国产av在线观看| 久久这里只有精品19| 亚洲精品成人av观看孕妇| 免费高清在线观看日韩| 黄片播放在线免费| 99re在线观看精品视频| 无人区码免费观看不卡| 麻豆久久精品国产亚洲av | 国产精品亚洲一级av第二区| 久久久国产精品麻豆| 久久国产亚洲av麻豆专区| 黄色成人免费大全| 久久人人精品亚洲av| 90打野战视频偷拍视频| 99香蕉大伊视频| 夜夜躁狠狠躁天天躁| 久久久久久免费高清国产稀缺| 深夜精品福利| 欧美日韩福利视频一区二区| 啦啦啦在线免费观看视频4| 精品电影一区二区在线| videosex国产| 三级毛片av免费| 亚洲欧美一区二区三区久久| 久久人人精品亚洲av| 在线看a的网站| 日韩欧美在线二视频| 一级黄色大片毛片| 长腿黑丝高跟| 国产日韩一区二区三区精品不卡| 日韩欧美一区视频在线观看| 99久久综合精品五月天人人| 日韩精品青青久久久久久| 在线观看午夜福利视频| 亚洲va日本ⅴa欧美va伊人久久| 日本黄色日本黄色录像| 丰满迷人的少妇在线观看| 丁香欧美五月| 精品一区二区三区四区五区乱码| 99久久久亚洲精品蜜臀av| 亚洲黑人精品在线| 757午夜福利合集在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 老鸭窝网址在线观看| 韩国av一区二区三区四区| 亚洲精品美女久久久久99蜜臀| 国产亚洲欧美精品永久| 波多野结衣高清无吗| a在线观看视频网站| 超碰成人久久| 丝袜美腿诱惑在线| 丝袜美足系列| 久久精品人人爽人人爽视色| 日韩大尺度精品在线看网址 | www.www免费av| 欧美乱色亚洲激情| 露出奶头的视频| xxx96com| 精品乱码久久久久久99久播| 亚洲欧洲精品一区二区精品久久久| 波多野结衣av一区二区av| 免费在线观看黄色视频的| 在线观看免费视频日本深夜| 日韩av在线大香蕉| 色婷婷久久久亚洲欧美| 国产伦人伦偷精品视频| 欧美日韩国产mv在线观看视频| 人人妻人人爽人人添夜夜欢视频| 成人18禁高潮啪啪吃奶动态图| 国产精品香港三级国产av潘金莲| 中文字幕色久视频| 波多野结衣av一区二区av| 热99re8久久精品国产| 欧美日韩黄片免| 1024视频免费在线观看| 亚洲精品一卡2卡三卡4卡5卡| 久久中文字幕一级| 国产熟女xx| 国产精品国产av在线观看| 国产国语露脸激情在线看| 日日爽夜夜爽网站| 怎么达到女性高潮| 久久久久久久久免费视频了| 波多野结衣高清无吗| 久久精品91无色码中文字幕| 精品高清国产在线一区| 国产91精品成人一区二区三区| 欧美日韩福利视频一区二区| 国产精品免费视频内射| 成人免费观看视频高清| 成年女人毛片免费观看观看9| 久久精品91蜜桃| av在线天堂中文字幕 | 手机成人av网站| 午夜福利在线免费观看网站| 我的亚洲天堂| 男女做爰动态图高潮gif福利片 | 精品无人区乱码1区二区| 一二三四在线观看免费中文在| 亚洲人成77777在线视频| 久久午夜综合久久蜜桃| 欧美av亚洲av综合av国产av| 香蕉国产在线看| 又黄又粗又硬又大视频| 亚洲男人天堂网一区| 午夜老司机福利片| 777久久人妻少妇嫩草av网站| 国产一区二区三区综合在线观看| 久久草成人影院| 久久久久久人人人人人| 老鸭窝网址在线观看| 国产成人av激情在线播放| 午夜影院日韩av| 久久久久国内视频| 午夜日韩欧美国产| 男人舔女人下体高潮全视频| 久久九九热精品免费| 亚洲一区二区三区不卡视频| 两个人免费观看高清视频| 国产精品美女特级片免费视频播放器 | 中国美女看黄片| 亚洲国产欧美网| 丰满的人妻完整版| 国产单亲对白刺激| 国产又爽黄色视频| 精品久久久久久成人av| 亚洲精品久久成人aⅴ小说| 老司机福利观看| 巨乳人妻的诱惑在线观看| 亚洲欧美激情综合另类| 亚洲美女黄片视频| 老司机午夜福利在线观看视频| 村上凉子中文字幕在线| 少妇的丰满在线观看| 亚洲精华国产精华精| 欧美精品亚洲一区二区| 18禁裸乳无遮挡免费网站照片 | 亚洲情色 制服丝袜| 男女做爰动态图高潮gif福利片 | 黑人操中国人逼视频| 国产精品爽爽va在线观看网站 | 桃色一区二区三区在线观看| 村上凉子中文字幕在线| 婷婷丁香在线五月| 日韩欧美一区视频在线观看| 搡老乐熟女国产| 国产欧美日韩一区二区三| 黄色丝袜av网址大全| 精品国产美女av久久久久小说| 黄色丝袜av网址大全| 巨乳人妻的诱惑在线观看| 99久久综合精品五月天人人| 国产高清激情床上av| 免费av中文字幕在线| 在线永久观看黄色视频| xxx96com| 国产午夜精品久久久久久| 日日摸夜夜添夜夜添小说| 美女高潮到喷水免费观看| 麻豆一二三区av精品| 一个人免费在线观看的高清视频| 欧美乱色亚洲激情| 69av精品久久久久久| 欧美日韩av久久| 日韩人妻精品一区2区三区| 国产精品免费视频内射| 久久性视频一级片| 国产欧美日韩一区二区三区在线| 亚洲欧美激情综合另类| 亚洲第一av免费看| 多毛熟女@视频| 亚洲激情在线av| 黑人欧美特级aaaaaa片| 亚洲精品国产精品久久久不卡| 丁香六月欧美| 香蕉丝袜av| 欧美国产精品va在线观看不卡| 高清黄色对白视频在线免费看| 亚洲色图av天堂| av超薄肉色丝袜交足视频| 国产精品98久久久久久宅男小说| 久久性视频一级片| 大型黄色视频在线免费观看| 日韩一卡2卡3卡4卡2021年| 色在线成人网| 欧美中文综合在线视频| 国产精品免费一区二区三区在线| 国产精品电影一区二区三区| 久久久国产精品麻豆| 精品久久久精品久久久| 欧美在线一区亚洲| 日韩中文字幕欧美一区二区| 日韩有码中文字幕| 亚洲少妇的诱惑av| 欧美另类亚洲清纯唯美| 精品乱码久久久久久99久播| 看片在线看免费视频| 大陆偷拍与自拍| 久久精品亚洲熟妇少妇任你| 午夜福利在线观看吧| 亚洲 欧美 日韩 在线 免费| 欧美日韩一级在线毛片| 久久人人精品亚洲av| 成人三级做爰电影| a级片在线免费高清观看视频| 精品久久久久久久久久免费视频 | 久久国产精品人妻蜜桃| 五月开心婷婷网| 好看av亚洲va欧美ⅴa在| 满18在线观看网站| 一区二区三区激情视频| 人人妻人人澡人人看| 久久久国产欧美日韩av| 国产精品 欧美亚洲| 国产成年人精品一区二区 | 亚洲av熟女| 97人妻天天添夜夜摸| 国产精华一区二区三区| 999久久久国产精品视频| 桃色一区二区三区在线观看| 久久亚洲精品不卡| 黄色怎么调成土黄色| 淫秽高清视频在线观看| 美女国产高潮福利片在线看| 两性夫妻黄色片| 成人国产一区最新在线观看| www.999成人在线观看| 人妻丰满熟妇av一区二区三区| 欧美日韩中文字幕国产精品一区二区三区 | 9色porny在线观看| 亚洲av日韩精品久久久久久密| 91九色精品人成在线观看| 51午夜福利影视在线观看| 精品福利观看| 99在线视频只有这里精品首页| 99久久人妻综合| 人妻久久中文字幕网| 91麻豆精品激情在线观看国产 | 18禁美女被吸乳视频| 变态另类成人亚洲欧美熟女 | 欧美最黄视频在线播放免费 | 一级毛片精品| 18禁美女被吸乳视频| 亚洲精品成人av观看孕妇| 99精品久久久久人妻精品| 村上凉子中文字幕在线| 久久久久精品国产欧美久久久| 深夜精品福利| 一本综合久久免费| a级毛片在线看网站| 免费在线观看影片大全网站| 首页视频小说图片口味搜索| 亚洲激情在线av| 免费在线观看日本一区| 99国产综合亚洲精品| 国产精品99久久99久久久不卡| 如日韩欧美国产精品一区二区三区| 亚洲伊人色综图| 国产激情久久老熟女| 99热只有精品国产| 精品卡一卡二卡四卡免费| 交换朋友夫妻互换小说| 老司机靠b影院| 首页视频小说图片口味搜索| 1024香蕉在线观看| 国产成人精品久久二区二区免费| 啪啪无遮挡十八禁网站| 黄色 视频免费看| 桃色一区二区三区在线观看| 免费av毛片视频| 看片在线看免费视频| 亚洲七黄色美女视频| 97人妻天天添夜夜摸| 亚洲一卡2卡3卡4卡5卡精品中文| 搡老岳熟女国产| 黑人欧美特级aaaaaa片| 欧美日韩亚洲综合一区二区三区_| 国产精品免费视频内射| 国产成人一区二区三区免费视频网站| av欧美777| 高潮久久久久久久久久久不卡| 欧美日韩亚洲综合一区二区三区_| 一二三四在线观看免费中文在| 亚洲人成电影免费在线| 一a级毛片在线观看| 黄色怎么调成土黄色| 欧美激情高清一区二区三区| 一区二区三区激情视频| 欧美丝袜亚洲另类 | 男人的好看免费观看在线视频 | 免费av中文字幕在线| 欧美激情极品国产一区二区三区| 欧美日韩瑟瑟在线播放| 国产精品成人在线| 色在线成人网| 国产在线精品亚洲第一网站| 中文字幕另类日韩欧美亚洲嫩草| 动漫黄色视频在线观看| 国产成人欧美在线观看| a级片在线免费高清观看视频| 精品国产国语对白av| 中国美女看黄片| 波多野结衣高清无吗| 成人18禁在线播放| 一个人免费在线观看的高清视频| 变态另类成人亚洲欧美熟女 | 精品久久久精品久久久| 99riav亚洲国产免费| 人人妻,人人澡人人爽秒播| 黑人猛操日本美女一级片| 欧美日韩中文字幕国产精品一区二区三区 | 看黄色毛片网站| 亚洲精华国产精华精| 欧美+亚洲+日韩+国产| 国产精品久久电影中文字幕| 国产黄a三级三级三级人| 日韩精品中文字幕看吧| 操出白浆在线播放| 国产av一区二区精品久久| 亚洲中文日韩欧美视频| 一级,二级,三级黄色视频| www国产在线视频色| 深夜精品福利| 在线十欧美十亚洲十日本专区| 国产色视频综合| 丁香欧美五月| 可以在线观看毛片的网站| 免费搜索国产男女视频| 亚洲av成人不卡在线观看播放网| 久久久久精品国产欧美久久久| 国产一区二区三区在线臀色熟女 | 国产精品秋霞免费鲁丝片| 国产蜜桃级精品一区二区三区| 精品卡一卡二卡四卡免费| 又大又爽又粗| 在线视频色国产色| 看免费av毛片| av网站免费在线观看视频| www.自偷自拍.com| 国产精品一区二区在线不卡| 中文字幕另类日韩欧美亚洲嫩草| 亚洲精品中文字幕在线视频| 校园春色视频在线观看| 久久精品人人爽人人爽视色| 亚洲人成网站在线播放欧美日韩| 亚洲成国产人片在线观看| 十八禁网站免费在线| 亚洲一区中文字幕在线| 一进一出抽搐动态| 一进一出好大好爽视频| 色播在线永久视频| 大香蕉久久成人网| 久久国产精品男人的天堂亚洲| 一级毛片精品| 亚洲一区二区三区欧美精品| 久久人人精品亚洲av| 怎么达到女性高潮| 一个人观看的视频www高清免费观看 | 久久精品91蜜桃| 99久久综合精品五月天人人| 亚洲中文av在线| 黑人巨大精品欧美一区二区mp4| 首页视频小说图片口味搜索| 日本黄色日本黄色录像| 亚洲精品一卡2卡三卡4卡5卡| e午夜精品久久久久久久| 91精品三级在线观看| 老司机在亚洲福利影院| 日本vs欧美在线观看视频| 亚洲久久久国产精品| 欧美日韩精品网址| 亚洲成a人片在线一区二区| www.999成人在线观看| 777久久人妻少妇嫩草av网站| 操美女的视频在线观看| 国产欧美日韩精品亚洲av| av免费在线观看网站| 十分钟在线观看高清视频www| 亚洲熟妇熟女久久| 亚洲欧美精品综合一区二区三区| 三级毛片av免费| 88av欧美| 精品福利观看| 亚洲免费av在线视频| 国产黄a三级三级三级人| 97人妻天天添夜夜摸| 亚洲一区高清亚洲精品| 91精品三级在线观看| 手机成人av网站| 很黄的视频免费| 国产欧美日韩精品亚洲av| 在线观看66精品国产| 老司机午夜福利在线观看视频| videosex国产| 日韩中文字幕欧美一区二区|