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

    網(wǎng)格變形控制技術(shù)在帶前置定子導(dǎo)管槳水動(dòng)力優(yōu)化設(shè)計(jì)中的應(yīng)用

    2011-07-11 09:39:52高同兵閆學(xué)琳陳小寧
    艦船科學(xué)技術(shù) 2011年12期
    關(guān)鍵詞:控制點(diǎn)控制技術(shù)網(wǎng)格

    劉 敏,張 寧,高同兵,閆學(xué)琳,陳小寧

    (武漢第二船舶設(shè)計(jì)研究所,湖北 武漢 430064)

    網(wǎng)格變形控制技術(shù)在帶前置定子導(dǎo)管槳水動(dòng)力優(yōu)化設(shè)計(jì)中的應(yīng)用

    劉 敏,張 寧,高同兵,閆學(xué)琳,陳小寧

    (武漢第二船舶設(shè)計(jì)研究所,湖北 武漢 430064)

    應(yīng)用網(wǎng)格變形控制技術(shù),對(duì)帶前置定子導(dǎo)管槳進(jìn)行了轉(zhuǎn)動(dòng)葉片水動(dòng)力外形的優(yōu)化設(shè)計(jì)研究,優(yōu)化方案的敞水效率較原始方案提高了2%,并能顯著減弱葉梢流動(dòng)的分離現(xiàn)象。優(yōu)化結(jié)果表明,基于網(wǎng)格變形控制技術(shù)的網(wǎng)格重構(gòu)方法能在保證網(wǎng)格相似性的條件下,解決多方案網(wǎng)格自動(dòng)生成的難題,便于進(jìn)行組合式推進(jìn)器的水動(dòng)力性能優(yōu)化設(shè)計(jì)。

    網(wǎng)格變形控制技術(shù);帶前置定子導(dǎo)管槳;水動(dòng)力優(yōu)化設(shè)計(jì)

    0 引言

    在過(guò)去的20年里,氣動(dòng)優(yōu)化設(shè)計(jì)(ADO)已經(jīng)成為計(jì)算流體力學(xué)應(yīng)用方面的一個(gè)重要研究領(lǐng)域。進(jìn)行優(yōu)化設(shè)計(jì)時(shí),需要對(duì)各方案進(jìn)行設(shè)計(jì)結(jié)果的預(yù)估評(píng)價(jià),其中CFD技術(shù)在流體性能優(yōu)化方面得到了廣泛應(yīng)用。然而,在對(duì)樣本進(jìn)行采樣分析時(shí),相應(yīng)的網(wǎng)格自動(dòng)化生成技術(shù)仍是目前的研究難點(diǎn)。網(wǎng)格生成不僅需要較大的計(jì)算資源和較長(zhǎng)的運(yùn)算時(shí)間,而且對(duì)參數(shù)變化范圍較大的樣本,極易失敗,從而導(dǎo)致樣本的缺失?;蚴巧傻木W(wǎng)格質(zhì)量較差,得到誤差較大的結(jié)果,使得優(yōu)化失敗。目前,進(jìn)行性能優(yōu)化的網(wǎng)格生成方法主要有面元法、網(wǎng)格模板和非結(jié)構(gòu)化網(wǎng)格3種。但這些方法均存在一定的局限性。如趙威[1]等采用面元法進(jìn)行螺旋槳性能優(yōu)化,避免了流域自動(dòng)網(wǎng)格生成的難題,但針對(duì)泵噴推進(jìn)器這樣的復(fù)雜結(jié)構(gòu),其適用性需要進(jìn)一步研討。Andrea[2]等采用CFD分析和優(yōu)化技術(shù)對(duì)高升力渦輪葉柵進(jìn)行了性能優(yōu)化,其中網(wǎng)格重構(gòu)采用網(wǎng)格模版的方法生成。對(duì)于幾何相似的樣本,可以采用網(wǎng)格模板的方法生成網(wǎng)格。然而對(duì)相似性較小的樣本,該方法容易生成帶負(fù)體積的計(jì)算網(wǎng)格,從而使得CFD模擬無(wú)法進(jìn)行。另外,ICEM-CFD針對(duì)結(jié)構(gòu)化網(wǎng)格提供了重復(fù)利用塊文件的網(wǎng)格生成方法,但因?yàn)樾枰付ㄍ負(fù)鋵?duì)應(yīng)關(guān)系而在自動(dòng)網(wǎng)格生成方面存在不足。目前,由于非結(jié)構(gòu)化網(wǎng)格對(duì)邊界的適應(yīng)性強(qiáng),而被廣泛用于 CFD優(yōu)化。劉艷[3]等在NACA0012翼型非結(jié)構(gòu)化網(wǎng)格的基礎(chǔ)上,采用動(dòng)彈網(wǎng)格技術(shù)生成翼型樣本網(wǎng)格。但該方法無(wú)法保證各樣本的網(wǎng)格相似性,特別是邊界層網(wǎng)格的生成。

    本文基于網(wǎng)格變形控制技術(shù)進(jìn)行了網(wǎng)格重構(gòu),由1套基本網(wǎng)格生成各方案相應(yīng)的變形網(wǎng)格,可以保證各方案網(wǎng)格的相似性,解決了多方案設(shè)計(jì)中網(wǎng)格自動(dòng)化生成的難題,便于進(jìn)行設(shè)計(jì)優(yōu)化。本文以文獻(xiàn)[4]中提供的帶定子導(dǎo)管槳為研究對(duì)象,首先對(duì)其水動(dòng)力性能進(jìn)行CFD計(jì)算,驗(yàn)證所采用網(wǎng)格及模擬方法的可行性;然后將葉片沿其積疊線進(jìn)行扭轉(zhuǎn),給定扭角范圍為-8°~8°并隨機(jī)生成樣本;再根據(jù)給定的扭角自動(dòng)生成相應(yīng)網(wǎng)格并進(jìn)行CFD性能計(jì)算;最后通過(guò)比較不同扭角下槳的水動(dòng)力性能,得到2°扭角下效率提高2%的優(yōu)化方案。該方案葉梢附近流場(chǎng)分析表明,葉梢的流動(dòng)分離得到了較好的控制。

    1 計(jì)算模型及數(shù)值模擬方法

    文獻(xiàn)[4]給出了帶定子導(dǎo)管槳外形幾何尺寸及轉(zhuǎn)子主要幾何參數(shù),并對(duì)其水動(dòng)力性能進(jìn)行了實(shí)驗(yàn)研究。其中,轉(zhuǎn)子和定子葉數(shù)分別為4葉和9葉,其三維幾何外形如圖1所示。計(jì)算域采用單流道﹑結(jié)構(gòu)化網(wǎng)格進(jìn)行劃分,葉片表面﹑葉頂及輪轂附近進(jìn)行網(wǎng)格加密處理。其中,葉片表面添加邊界層網(wǎng)格,網(wǎng)格總量約110萬(wàn),如圖2所示。采用Fluent軟件求解Navier-Stokes方程,計(jì)算區(qū)域劃分為動(dòng)、靜2個(gè)域,采用MRF方法建模,轉(zhuǎn)/靜交界面采用Mixing-Plane模型進(jìn)行定常模擬,得到轉(zhuǎn)子葉片的推力和扭矩,并計(jì)算得到其敞水效率。

    圖1 帶定子導(dǎo)管槳正視圖Fig.1 Display+Z view of ducted propeller

    圖2 導(dǎo)管槳單流道CFD網(wǎng)格Fig.2 Single flow passage mesh

    2 網(wǎng)格變形控制技術(shù)

    采用4×4×4控制點(diǎn)控制轉(zhuǎn)子區(qū)域網(wǎng)格,其中葉片流域進(jìn)口﹑出口﹑周期性邊界上控制點(diǎn)設(shè)置為固定控制點(diǎn),移動(dòng)控制點(diǎn)設(shè)置為圍繞葉片的2×2×4個(gè)控制點(diǎn),如圖3所示。采用樣條曲線連接控制點(diǎn),保證控制曲線的光滑,當(dāng)控制點(diǎn)移動(dòng)時(shí),空間樣條曲線擬和控制點(diǎn),整體網(wǎng)格依次移動(dòng),保證了整體網(wǎng)格的光滑過(guò)渡。沿輪蓋表面移動(dòng)葉頂?shù)?個(gè)控制點(diǎn),可以實(shí)現(xiàn)葉片的前彎及后彎,如圖4和圖5所示。沿周向移動(dòng)圍繞葉片的網(wǎng)格控制點(diǎn),可以實(shí)現(xiàn)葉片的扭轉(zhuǎn),如圖6所示。本文采用的網(wǎng)格變形控制技術(shù)基于空間網(wǎng)格進(jìn)行直接變形,保證了葉片表面網(wǎng)格的光順及整體網(wǎng)格質(zhì)量,并可以保證變形網(wǎng)格與原始網(wǎng)格擁有相同的網(wǎng)格數(shù)量及拓?fù)潢P(guān)系,從而確保了各方案網(wǎng)格的相似性。

    3 性能計(jì)算及分析

    本文對(duì)J=0.4﹑0.6﹑0.8﹑1.0等4個(gè)工況進(jìn)行了CFD計(jì)算,并與實(shí)驗(yàn)數(shù)據(jù)進(jìn)行了對(duì)比,如圖7所示。從圖中可以看到,推力系數(shù)與實(shí)驗(yàn)結(jié)果吻合較好,表明本文所采用模擬方法的正確性。扭矩系數(shù)在高進(jìn)速條件下誤差控制在4.5%以內(nèi),滿足工程要求。

    圖8所示為轉(zhuǎn)子葉片附近流線圖。區(qū)域1所示為葉頂分離流,區(qū)域2所示為轉(zhuǎn)子輪轂附近尾跡區(qū)域流線。在區(qū)域1位置,存在明顯的葉頂分離流動(dòng)。流線在葉片迎流邊附近即發(fā)生流動(dòng)脫離。在區(qū)域2位置,存在部分低速流體的明顯堆積,輪轂附近流線從壓力面回流到葉片吸力面,與吸力面流體發(fā)生摻混,形成1個(gè)回流漩渦區(qū)再流到下游。圖9和圖10分別為轉(zhuǎn)子葉片吸力面與壓力面表面靜壓分布等值線圖。在轉(zhuǎn)子葉片導(dǎo)邊附近,吸力面上存在1個(gè)明顯的高壓區(qū),而壓力面上存在1個(gè)低壓區(qū),這對(duì)槳推力帶來(lái)不利影響,直接影響其效率。

    圖9 轉(zhuǎn)子葉片吸力面靜壓分布等值線Fig.9 Pressure contour of blade suction

    圖10 轉(zhuǎn)子葉片壓力面靜壓分布等值線Fig.10 Pressure contour of blade pressure

    4 設(shè)計(jì)優(yōu)化及分析

    在造型設(shè)計(jì)過(guò)程中,影響三維葉片造型的參數(shù)可分為氣動(dòng)參數(shù)(如相對(duì)氣流角﹑氣流脫軌角等)和幾何參數(shù)(如葉片的彎﹑掠﹑扭﹑葉截面中弧線拱度,厚度分布等),所涉及的參數(shù)較多。在給定導(dǎo)管槳外形的條件下,針對(duì)葉片進(jìn)行性能優(yōu)化時(shí),目前研究大多針對(duì)其中比較重要的參數(shù)進(jìn)行優(yōu)化,如Ahn C S[5]等通過(guò)優(yōu)化積疊軸,調(diào)整葉片的彎掠特性來(lái)提高壓氣機(jī)性能。在對(duì)截面進(jìn)行參數(shù)化優(yōu)化時(shí),由于各截面之間的參數(shù)較多,約束關(guān)系難以給定,容易造成造型失敗,如在某個(gè)截面出現(xiàn)零厚度等非物理解。因此,部分研究使用自由曲面函數(shù)(如NURBS函數(shù)等)直接進(jìn)行三維空間曲面造型并進(jìn)行優(yōu)化。袁新[6]采用NURBS曲面對(duì)軸流壓氣機(jī)葉片進(jìn)行直接幾何造型,并集成CFD軟件進(jìn)行優(yōu)化設(shè)計(jì),完成了多級(jí)軸流壓氣機(jī)的優(yōu)化設(shè)計(jì)。

    通過(guò)對(duì)圖8~圖10中轉(zhuǎn)子葉片流線及靜壓結(jié)果的分析可知,為了提高轉(zhuǎn)子水動(dòng)力性能,并減少葉片流動(dòng)分離和堆積,可通過(guò)改變?nèi)~片安裝角(扭)及輪轂附近葉截面拱度來(lái)實(shí)現(xiàn)。這樣極大地減少了優(yōu)化參數(shù),便于進(jìn)行優(yōu)化。

    本文主要采用改變?nèi)~片扭角進(jìn)行優(yōu)化。將葉型離散為4個(gè)葉截面(如圖3所示),通過(guò)改變?nèi)~片控制點(diǎn)周向位置改變?nèi)~片扭角,相應(yīng)的網(wǎng)格采用控制網(wǎng)格變形技術(shù)實(shí)現(xiàn)。圖11給出了葉頂部截面在不同扭角下葉形坐標(biāo)對(duì)比。扭角變化范圍為-8°~8°,以30%弦長(zhǎng)處為葉片積疊線即葉截面旋轉(zhuǎn)中心,以不同扭角隨機(jī)生成樣本,采用CFD模擬得到其水動(dòng)力性能。計(jì)算所采用的進(jìn)速系數(shù)為0.6,即設(shè)計(jì)進(jìn)速點(diǎn)。本文給出了20個(gè)隨機(jī)樣本的計(jì)算結(jié)果,其水動(dòng)力性能點(diǎn)隨扭角的變化曲線如圖12所示。可見(jiàn)隨著扭角的增加推力系數(shù)降低,同時(shí)扭矩系數(shù)也降低。效率曲線隨扭角在2.2°附近達(dá)到最大后隨即降低,其值相比于原型(0°)增加2.4%。在-2°~2°扭角之間,曲線變化平緩,偏離該范圍后效率曲線梯度較大。

    圖13所示為扭角分別為-6°,2°和6°三種情況下葉頂附近流線。在扭角為-6°時(shí),流線在葉片吸力面發(fā)生明顯分離,主要是因?yàn)閬?lái)流攻角過(guò)大,導(dǎo)致葉頂附近流動(dòng)惡化;當(dāng)扭角為2°時(shí),在設(shè)計(jì)工況點(diǎn)流線未見(jiàn)分離,流動(dòng)較理想;在扭角為6°時(shí),流動(dòng)在壓力面發(fā)生分離,并生成較明顯的旋渦。根據(jù)機(jī)翼理論,當(dāng)扭角過(guò)大或過(guò)小時(shí),流動(dòng)均在葉頂葉片表面發(fā)生明顯分離,所選取的攻角不合理將導(dǎo)致性能下降。綜合圖12和圖13的分析結(jié)果,盡管在扭角為2.2°附近效率曲線取得最大值,但考慮到工況偏移后將導(dǎo)致攻角變化,而效率曲線在2.2°扭角上斜率較大,為了保證較寬廣的工況范圍,因此優(yōu)化方案的扭角取為2°。

    圖14 優(yōu)化槳與原型槳性能曲線Fig.14 Performance curves of improved propeller and original propeller

    圖14所示為采用2°扭角后優(yōu)化槳與原型槳的性能曲線對(duì)比。在設(shè)計(jì)進(jìn)速點(diǎn),效率得到約2%的提升;在高進(jìn)速情況下,由于偏離優(yōu)化槳設(shè)計(jì)工況,原型槳相反效率較高。

    5 結(jié)語(yǔ)

    本文對(duì)帶定子導(dǎo)管槳水動(dòng)力性能進(jìn)行設(shè)計(jì)優(yōu)化,采用控制網(wǎng)格變形技術(shù)對(duì)葉片扭角進(jìn)行多方案對(duì)比。結(jié)果表明,采用2°扭角的優(yōu)化方案可以明顯控制葉頂流動(dòng)分離,設(shè)計(jì)點(diǎn)敞水效率可較原始方案提高2%,并且在性能曲線上設(shè)計(jì)點(diǎn)附近的效率均可得到提高。本文的研究成果表明,基于網(wǎng)格變形控制技術(shù)的網(wǎng)格重構(gòu)方法,能在保證網(wǎng)格相似性的條件下,解決多方案網(wǎng)格自動(dòng)生成的難題,便于進(jìn)行組合式推進(jìn)器的水動(dòng)力性能優(yōu)化設(shè)計(jì)。

    [1]趙威,楊晨俊.船舶螺旋槳螺距及拱度的優(yōu)化設(shè)計(jì)研究[J].中國(guó)造船,2010,51(1):1 -8.

    ZHAO Wei,YANG Chen-jun.Reasearch on optimized design of the pitch and camber of marine propellers[J].Shipbuilding of China,2010,51(1):1 -8.

    [2]ARNONE A,BONAIUTI D,F(xiàn)OCACCI A,et al.Parametric optimization of a high-lift turbine vane[J].ASME GT -2004 -54308,2004.

    [3]劉艷,王江峰,伍貽兆.基于遺傳算法的翼型多目標(biāo)氣動(dòng)優(yōu)化設(shè)計(jì)[J].飛機(jī)設(shè)計(jì),2008,28(4):1 -5.

    LIU Yan,WANG Jiang-feng,WU Yi-zhao.Multi-objective aerodynamic optimization design of airfoil based on genetic algorithm[J].Aircraft Design,200,28(4):1 -5.

    [4]HUGHES M J,KINNAS S A.An analysis method for a ducted propellerwith pre-swirlstator blades[C].propeller/shafting 91 symposium virginia beach,virginia,september,1991.17 -18.

    [5]AHN C S,KIM K Y.Aerodynamic design optimization of an axial flow compressor roto[J].r.ASME GT - 2002 -30445,2002.

    [6]陳志鵬,袁新.全工況性能優(yōu)化在壓氣機(jī)多級(jí)環(huán)境中的應(yīng)用[J].工程熱物理學(xué)報(bào),2008,29(2):221 -224.

    CHEN Zhi-peng,YUAN Xin.Aerodynamics optimum design of compressor blades based on numerical overall performance calculation[J].JournalofEngineering Thermophysics,2008,29(2):221-224.

    Application of controled mesh deformation technology in hydrodynamic optimization design for a ducted propeller with pre-swirl stator

    LIU Min,ZHANG Ning,GAO Tong-bing,YAN Xue-lin,CHEN Xiao-ning
    (Wuhan Second Ship Design and Research Institute,Wuhan 430064,China)

    This paper applied controled mesh deformation technology to hydrodynamic profile optimization of a ducted propeller with pre-swirl stator.Open water efficiency of the optimization scheme is 2%higher than the initial scheme,and the flow separation near the blade tip is highly broken.The results of this paper indicates that mesh reconstruction method based on controled mesh deformation technology can settle the problem in the mesh automatic generation for multi-scheme design in the condition of mesh similarity and can be used in hydrodynamic optimization design for integrated propulsor.

    controled mesh deformation technology;ducted propeller with pre-swirl stator;hydrodynamic optimization design

    U664.34

    A

    1672-7649(2011)12-0019-05

    10.3404/j.issn.1672-7649.2011.12.004

    2010-12-02;

    2011-05-27

    劉敏(1980-),男,工程師,主要研究方向?yàn)榇翱傮w性能及推進(jìn)器設(shè)計(jì)。

    猜你喜歡
    控制點(diǎn)控制技術(shù)網(wǎng)格
    用全等三角形破解網(wǎng)格題
    反射的橢圓隨機(jī)偏微分方程的網(wǎng)格逼近
    對(duì)工程建設(shè)中的機(jī)械自動(dòng)化控制技術(shù)探討
    基于BTN8962TA的PVG32比例閥控制技術(shù)
    電子制作(2019年7期)2019-04-25 13:16:50
    重疊網(wǎng)格裝配中的一種改進(jìn)ADT搜索方法
    NFFD控制點(diǎn)分布對(duì)氣動(dòng)外形優(yōu)化的影響
    基于風(fēng)險(xiǎn)管理下的項(xiàng)目建設(shè)內(nèi)部控制點(diǎn)思考
    基于曲面展開(kāi)的自由曲面網(wǎng)格劃分
    景觀照明聯(lián)動(dòng)控制技術(shù)的展望
    相似材料模型中控制點(diǎn)像點(diǎn)坐標(biāo)定位研究
    中文字幕高清在线视频| 久久久久久大精品| 亚洲五月天丁香| 成年免费大片在线观看| 日日摸夜夜添夜夜添小说| 亚洲自偷自拍图片 自拍| 午夜福利高清视频| 国产av不卡久久| 国产成人精品久久二区二区91| 国产成人aa在线观看| 十八禁人妻一区二区| 精品久久蜜臀av无| 一本大道久久a久久精品| 黑人巨大精品欧美一区二区mp4| ponron亚洲| 久久精品国产清高在天天线| 又粗又爽又猛毛片免费看| 国产主播在线观看一区二区| 国产精品av久久久久免费| 久久香蕉激情| www.熟女人妻精品国产| 老熟妇乱子伦视频在线观看| 国产精品香港三级国产av潘金莲| 搡老妇女老女人老熟妇| 黄色片一级片一级黄色片| 色尼玛亚洲综合影院| 欧美性长视频在线观看| 色在线成人网| 久久久久久久午夜电影| 一级作爱视频免费观看| 成人一区二区视频在线观看| 三级男女做爰猛烈吃奶摸视频| 村上凉子中文字幕在线| 免费观看人在逋| 99在线视频只有这里精品首页| 亚洲国产精品久久男人天堂| 好看av亚洲va欧美ⅴa在| 伊人久久大香线蕉亚洲五| 1024香蕉在线观看| а√天堂www在线а√下载| 性欧美人与动物交配| 波多野结衣巨乳人妻| 国产亚洲精品久久久久久毛片| 中文字幕人成人乱码亚洲影| 老鸭窝网址在线观看| 成年女人毛片免费观看观看9| 夜夜躁狠狠躁天天躁| 婷婷亚洲欧美| 久久精品夜夜夜夜夜久久蜜豆 | 中文在线观看免费www的网站 | 91国产中文字幕| 欧美一级毛片孕妇| 一级毛片高清免费大全| 久久欧美精品欧美久久欧美| 听说在线观看完整版免费高清| 亚洲人成电影免费在线| 日日摸夜夜添夜夜添小说| 一区二区三区高清视频在线| 色综合站精品国产| xxxwww97欧美| 国产伦在线观看视频一区| 免费在线观看影片大全网站| 长腿黑丝高跟| 制服诱惑二区| 日日爽夜夜爽网站| 一区二区三区激情视频| 18禁国产床啪视频网站| 久久性视频一级片| 正在播放国产对白刺激| 婷婷精品国产亚洲av| 亚洲美女黄片视频| 亚洲av第一区精品v没综合| 久久久久久亚洲精品国产蜜桃av| 欧美乱色亚洲激情| 美女大奶头视频| 淫秽高清视频在线观看| 别揉我奶头~嗯~啊~动态视频| 午夜久久久久精精品| 欧美激情久久久久久爽电影| av天堂在线播放| 长腿黑丝高跟| 亚洲成av人片免费观看| 两个人看的免费小视频| 国语自产精品视频在线第100页| 日韩欧美国产一区二区入口| 欧美性猛交黑人性爽| 99在线视频只有这里精品首页| 男女那种视频在线观看| 国产精品99久久99久久久不卡| 波多野结衣高清无吗| 久久香蕉国产精品| 99热6这里只有精品| 国产激情欧美一区二区| 又紧又爽又黄一区二区| 亚洲 国产 在线| 日韩精品免费视频一区二区三区| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲熟妇中文字幕五十中出| 露出奶头的视频| 亚洲精华国产精华精| 国产成人av教育| 欧美成狂野欧美在线观看| 亚洲人与动物交配视频| 国产成人一区二区三区免费视频网站| 欧美黄色淫秽网站| 精品久久蜜臀av无| 欧美另类亚洲清纯唯美| 国产精品 欧美亚洲| 欧美极品一区二区三区四区| 观看免费一级毛片| 欧美大码av| 国产人伦9x9x在线观看| 高潮久久久久久久久久久不卡| 最近最新中文字幕大全电影3| 伦理电影免费视频| 亚洲在线自拍视频| 啦啦啦观看免费观看视频高清| 亚洲午夜理论影院| 在线十欧美十亚洲十日本专区| 亚洲中文字幕日韩| 99在线人妻在线中文字幕| 欧美日韩福利视频一区二区| 欧美日韩精品网址| avwww免费| av在线天堂中文字幕| 可以在线观看的亚洲视频| 精品第一国产精品| 全区人妻精品视频| 亚洲国产欧洲综合997久久,| 午夜老司机福利片| 91在线观看av| 97人妻精品一区二区三区麻豆| 香蕉丝袜av| 少妇的丰满在线观看| 亚洲成av人片在线播放无| 老司机靠b影院| 在线免费观看的www视频| 欧美成狂野欧美在线观看| 怎么达到女性高潮| 婷婷精品国产亚洲av在线| 美女高潮喷水抽搐中文字幕| 久久久久性生活片| 欧美日韩国产亚洲二区| 欧美又色又爽又黄视频| 长腿黑丝高跟| 变态另类丝袜制服| 色哟哟哟哟哟哟| www日本在线高清视频| 精品久久久久久久人妻蜜臀av| 啦啦啦观看免费观看视频高清| 欧美激情久久久久久爽电影| 12—13女人毛片做爰片一| 国产精品影院久久| 日本精品一区二区三区蜜桃| 在线观看美女被高潮喷水网站 | 听说在线观看完整版免费高清| 亚洲五月婷婷丁香| 国产麻豆成人av免费视频| 欧美国产日韩亚洲一区| 国产激情欧美一区二区| 国产精品亚洲一级av第二区| 色在线成人网| 天天添夜夜摸| 欧美乱妇无乱码| 天堂av国产一区二区熟女人妻 | 亚洲欧美日韩高清专用| 成人永久免费在线观看视频| 美女午夜性视频免费| 久久精品成人免费网站| 国产在线精品亚洲第一网站| 天天躁狠狠躁夜夜躁狠狠躁| 男插女下体视频免费在线播放| 国产精品野战在线观看| 搡老岳熟女国产| www.www免费av| 国产精品野战在线观看| 国产精品香港三级国产av潘金莲| 好男人在线观看高清免费视频| 亚洲国产欧美网| 国产成人系列免费观看| 麻豆成人午夜福利视频| 香蕉久久夜色| 一进一出好大好爽视频| 97超级碰碰碰精品色视频在线观看| 欧美精品啪啪一区二区三区| 国产精品,欧美在线| 操出白浆在线播放| 婷婷精品国产亚洲av| www.999成人在线观看| 久久久精品国产亚洲av高清涩受| 国产黄片美女视频| 久久久久久大精品| 欧美大码av| 国产三级中文精品| 夜夜爽天天搞| 精品国内亚洲2022精品成人| 国产高清激情床上av| 成人手机av| 亚洲av成人一区二区三| 欧美午夜高清在线| 国产日本99.免费观看| 国产精品一区二区精品视频观看| 欧美一区二区国产精品久久精品 | 人人妻,人人澡人人爽秒播| 最近最新中文字幕大全免费视频| xxx96com| 国产精品影院久久| 1024视频免费在线观看| 国产1区2区3区精品| 色在线成人网| a在线观看视频网站| 啦啦啦韩国在线观看视频| 国产野战对白在线观看| 丝袜人妻中文字幕| 欧美日本视频| 婷婷六月久久综合丁香| 麻豆成人午夜福利视频| 757午夜福利合集在线观看| 国产91精品成人一区二区三区| 国产精品美女特级片免费视频播放器 | 国产精品野战在线观看| 少妇被粗大的猛进出69影院| 99久久国产精品久久久| 变态另类成人亚洲欧美熟女| 麻豆国产97在线/欧美 | 亚洲欧美日韩无卡精品| 1024手机看黄色片| 国产精品1区2区在线观看.| 床上黄色一级片| 国产成人欧美在线观看| 日韩欧美 国产精品| 男女之事视频高清在线观看| 午夜日韩欧美国产| 国产精品免费一区二区三区在线| 精品久久蜜臀av无| 日韩精品青青久久久久久| 大型av网站在线播放| 看片在线看免费视频| 日韩 欧美 亚洲 中文字幕| 1024手机看黄色片| 久热爱精品视频在线9| 亚洲熟妇熟女久久| 精品无人区乱码1区二区| 日本熟妇午夜| 久久香蕉国产精品| 国产人伦9x9x在线观看| 国产亚洲精品第一综合不卡| 青草久久国产| 欧美日韩瑟瑟在线播放| 听说在线观看完整版免费高清| 亚洲av成人精品一区久久| 国产亚洲精品第一综合不卡| 日韩有码中文字幕| 成人特级黄色片久久久久久久| 嫩草影院精品99| 最新美女视频免费是黄的| 国产男靠女视频免费网站| 麻豆av在线久日| 午夜日韩欧美国产| 久久久水蜜桃国产精品网| 久久精品91蜜桃| 亚洲狠狠婷婷综合久久图片| 人妻丰满熟妇av一区二区三区| 五月伊人婷婷丁香| 午夜老司机福利片| 欧美三级亚洲精品| 久久久精品欧美日韩精品| 最新美女视频免费是黄的| 亚洲真实伦在线观看| 国产精华一区二区三区| 亚洲一区中文字幕在线| 中文字幕熟女人妻在线| 最新美女视频免费是黄的| 成人国产一区最新在线观看| 青草久久国产| 麻豆成人av在线观看| 黄色视频,在线免费观看| 成人18禁高潮啪啪吃奶动态图| 午夜福利18| 男人的好看免费观看在线视频 | 中亚洲国语对白在线视频| 亚洲av片天天在线观看| 亚洲av成人av| 成人18禁在线播放| 18禁黄网站禁片免费观看直播| 毛片女人毛片| 变态另类成人亚洲欧美熟女| 久久精品国产亚洲av高清一级| 精品国内亚洲2022精品成人| 久久精品国产清高在天天线| 91麻豆av在线| 国产精品影院久久| 国产欧美日韩一区二区精品| 久热爱精品视频在线9| 91av网站免费观看| 欧美乱码精品一区二区三区| 久久久久久久午夜电影| 亚洲 欧美一区二区三区| 美女黄网站色视频| 亚洲第一欧美日韩一区二区三区| 日韩欧美国产一区二区入口| 可以在线观看的亚洲视频| 精品一区二区三区视频在线观看免费| 久久久久国内视频| 亚洲专区国产一区二区| cao死你这个sao货| 岛国在线观看网站| 国产一区二区三区在线臀色熟女| 人人妻人人看人人澡| 国产精品一及| 欧美成人免费av一区二区三区| 男人舔奶头视频| www.www免费av| 曰老女人黄片| 男女视频在线观看网站免费 | 亚洲国产欧美一区二区综合| 精品国内亚洲2022精品成人| 亚洲熟妇中文字幕五十中出| 一级毛片女人18水好多| 五月玫瑰六月丁香| 欧美国产日韩亚洲一区| 国产高清激情床上av| 两人在一起打扑克的视频| 老汉色∧v一级毛片| 欧美成人性av电影在线观看| 99国产精品一区二区三区| 一二三四社区在线视频社区8| 欧美中文综合在线视频| 成人三级黄色视频| 久久这里只有精品19| 久久久久久免费高清国产稀缺| 国产精品 国内视频| 亚洲av日韩精品久久久久久密| 久久久久久亚洲精品国产蜜桃av| 日韩欧美国产一区二区入口| 久久这里只有精品19| 99久久久亚洲精品蜜臀av| 久久久久九九精品影院| 中文亚洲av片在线观看爽| 欧美成人性av电影在线观看| 日韩国内少妇激情av| 欧美日韩中文字幕国产精品一区二区三区| 免费一级毛片在线播放高清视频| 亚洲精品美女久久久久99蜜臀| 国产久久久一区二区三区| 美女午夜性视频免费| 性欧美人与动物交配| av天堂在线播放| 日本一区二区免费在线视频| 18禁观看日本| 午夜老司机福利片| 久久久久亚洲av毛片大全| 操出白浆在线播放| 一区二区三区高清视频在线| 亚洲人成电影免费在线| 天天躁夜夜躁狠狠躁躁| 别揉我奶头~嗯~啊~动态视频| 欧美日韩黄片免| 久久久久国内视频| 女同久久另类99精品国产91| 久久天躁狠狠躁夜夜2o2o| 国产精品一区二区三区四区免费观看 | 悠悠久久av| 一本久久中文字幕| 黑人巨大精品欧美一区二区mp4| 亚洲黑人精品在线| 亚洲一区中文字幕在线| 男女视频在线观看网站免费 | 最近在线观看免费完整版| av免费在线观看网站| 极品教师在线免费播放| 一本综合久久免费| 黄频高清免费视频| 色老头精品视频在线观看| 国产黄a三级三级三级人| 日韩欧美国产在线观看| 久久精品国产亚洲av香蕉五月| 欧美另类亚洲清纯唯美| 国产精品av视频在线免费观看| 亚洲第一欧美日韩一区二区三区| 亚洲成人中文字幕在线播放| 久久久国产成人精品二区| 2021天堂中文幕一二区在线观| 国产成人av教育| 亚洲专区中文字幕在线| bbb黄色大片| 999精品在线视频| 1024手机看黄色片| 日本免费一区二区三区高清不卡| 欧美日韩福利视频一区二区| 男人舔奶头视频| 国产成人一区二区三区免费视频网站| 国产av在哪里看| aaaaa片日本免费| 少妇粗大呻吟视频| 免费在线观看黄色视频的| 一区二区三区高清视频在线| 久久久精品国产亚洲av高清涩受| 欧美日本亚洲视频在线播放| 精品久久蜜臀av无| 久久这里只有精品中国| 成人18禁在线播放| 99热只有精品国产| 亚洲av中文字字幕乱码综合| 欧美性猛交╳xxx乱大交人| 韩国av一区二区三区四区| 亚洲人成77777在线视频| 欧美高清成人免费视频www| 亚洲av日韩精品久久久久久密| 精品高清国产在线一区| 免费高清视频大片| 国产av麻豆久久久久久久| 国产精品日韩av在线免费观看| 天堂av国产一区二区熟女人妻 | 最近最新中文字幕大全电影3| 美女黄网站色视频| 啦啦啦韩国在线观看视频| 久久精品综合一区二区三区| netflix在线观看网站| 叶爱在线成人免费视频播放| 亚洲va日本ⅴa欧美va伊人久久| 岛国在线免费视频观看| 日韩精品中文字幕看吧| 母亲3免费完整高清在线观看| 亚洲一区中文字幕在线| 中亚洲国语对白在线视频| 在线观看www视频免费| 男插女下体视频免费在线播放| 国产99白浆流出| 午夜精品久久久久久毛片777| 麻豆久久精品国产亚洲av| 真人一进一出gif抽搐免费| 午夜视频精品福利| 国产视频内射| 国产成人系列免费观看| 精品久久久久久久久久久久久| 91麻豆精品激情在线观看国产| 精品久久久久久成人av| 草草在线视频免费看| 精品第一国产精品| 亚洲无线在线观看| 亚洲最大成人中文| 桃红色精品国产亚洲av| 午夜福利高清视频| 桃色一区二区三区在线观看| 少妇人妻一区二区三区视频| 欧美极品一区二区三区四区| 91大片在线观看| 欧美午夜高清在线| 国产99久久九九免费精品| 亚洲av电影在线进入| 中文字幕av在线有码专区| 亚洲国产欧洲综合997久久,| 在线国产一区二区在线| 亚洲aⅴ乱码一区二区在线播放 | 国产麻豆成人av免费视频| 一进一出抽搐gif免费好疼| av超薄肉色丝袜交足视频| 国产免费男女视频| 欧美一区二区国产精品久久精品 | 色在线成人网| 女警被强在线播放| 国产99白浆流出| 俺也久久电影网| 国产亚洲精品av在线| 亚洲,欧美精品.| 亚洲天堂国产精品一区在线| 午夜免费观看网址| 高清在线国产一区| 男女做爰动态图高潮gif福利片| 久久精品国产99精品国产亚洲性色| 老司机靠b影院| 久9热在线精品视频| 在线观看免费视频日本深夜| 成人午夜高清在线视频| 国产成+人综合+亚洲专区| 亚洲狠狠婷婷综合久久图片| 窝窝影院91人妻| 日本精品一区二区三区蜜桃| 淫妇啪啪啪对白视频| 久久天堂一区二区三区四区| 久久人妻av系列| 五月玫瑰六月丁香| 50天的宝宝边吃奶边哭怎么回事| 99re在线观看精品视频| 成人国语在线视频| av天堂在线播放| 久久久久性生活片| 久久欧美精品欧美久久欧美| 国产麻豆成人av免费视频| 免费在线观看影片大全网站| 一级片免费观看大全| 男男h啪啪无遮挡| 亚洲18禁久久av| 黄色丝袜av网址大全| 一a级毛片在线观看| 亚洲免费av在线视频| 999久久久国产精品视频| 国产成人影院久久av| 欧美日韩精品网址| 我要搜黄色片| 啦啦啦韩国在线观看视频| 欧美黄色片欧美黄色片| √禁漫天堂资源中文www| 99国产极品粉嫩在线观看| 久久午夜亚洲精品久久| 婷婷精品国产亚洲av| 国产高清有码在线观看视频 | 午夜福利18| 非洲黑人性xxxx精品又粗又长| 国产精品久久久久久久电影 | 国产欧美日韩一区二区三| 窝窝影院91人妻| 精品免费久久久久久久清纯| 日韩中文字幕欧美一区二区| 亚洲免费av在线视频| 少妇被粗大的猛进出69影院| 亚洲av中文字字幕乱码综合| 99国产综合亚洲精品| 中文字幕av在线有码专区| 久久精品综合一区二区三区| 亚洲中文日韩欧美视频| 女生性感内裤真人,穿戴方法视频| 日本一本二区三区精品| 99在线人妻在线中文字幕| 中文字幕久久专区| 国产精品一及| 中文字幕最新亚洲高清| 国产在线观看jvid| 国产日本99.免费观看| 伊人久久大香线蕉亚洲五| 狂野欧美白嫩少妇大欣赏| 午夜福利高清视频| 桃红色精品国产亚洲av| 亚洲 欧美 日韩 在线 免费| а√天堂www在线а√下载| 99国产精品一区二区三区| 免费观看人在逋| 国产伦一二天堂av在线观看| videosex国产| 久久 成人 亚洲| 国产三级黄色录像| 中文字幕久久专区| 亚洲精品美女久久久久99蜜臀| 黑人欧美特级aaaaaa片| 少妇被粗大的猛进出69影院| 国产精品 欧美亚洲| 久久精品国产亚洲av香蕉五月| 欧美中文综合在线视频| av欧美777| 久久伊人香网站| 国产成人系列免费观看| 亚洲欧美精品综合久久99| 97人妻精品一区二区三区麻豆| 色综合欧美亚洲国产小说| 一区福利在线观看| 高清在线国产一区| 久久中文看片网| 欧美中文综合在线视频| 好看av亚洲va欧美ⅴa在| 叶爱在线成人免费视频播放| 日本一本二区三区精品| 亚洲熟妇中文字幕五十中出| 操出白浆在线播放| 亚洲中文字幕一区二区三区有码在线看 | 国产一区二区在线观看日韩 | 最近最新中文字幕大全电影3| 国产精品久久久久久精品电影| 午夜久久久久精精品| 欧美av亚洲av综合av国产av| 老司机福利观看| 2021天堂中文幕一二区在线观| 成人av在线播放网站| 搞女人的毛片| 性色av乱码一区二区三区2| 精品国产乱子伦一区二区三区| 男人舔女人下体高潮全视频| 亚洲avbb在线观看| 国产主播在线观看一区二区| 亚洲国产精品久久男人天堂| 久久精品国产亚洲av香蕉五月| 国产激情久久老熟女| 久久午夜亚洲精品久久| 久久精品人妻少妇| 一级片免费观看大全| 日韩中文字幕欧美一区二区| av福利片在线| 久久精品人妻少妇| 久久这里只有精品19| 午夜福利18| 国产激情久久老熟女| 黄色女人牲交| 午夜精品一区二区三区免费看| АⅤ资源中文在线天堂| 一区二区三区国产精品乱码| 国产精品av久久久久免费| 国产午夜精品论理片| 久久亚洲真实| 三级国产精品欧美在线观看 | 国产av一区在线观看免费| 亚洲国产欧洲综合997久久,| 国产午夜精品论理片| 这个男人来自地球电影免费观看| 伊人久久大香线蕉亚洲五| 两人在一起打扑克的视频| 久久精品夜夜夜夜夜久久蜜豆 | 国产高清视频在线播放一区| 国产区一区二久久| 午夜a级毛片| 婷婷丁香在线五月| 免费看日本二区| 悠悠久久av| 国产97色在线日韩免费| 亚洲美女黄片视频| 欧美大码av|