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

    水下航行體超空泡外形控制的數(shù)值研究

    2011-09-03 11:57:28胡曉郜冶熊永亮
    關(guān)鍵詞:錐角空泡空化

    胡曉,郜冶,熊永亮

    (哈爾濱工程大學(xué)航天與建筑工程學(xué)院,黑龍江哈爾濱150001)

    航行體在高速入水和水下高速運(yùn)動過程中,會產(chǎn)生空泡(或超空泡),空泡的外形和尺寸會對航行體的水動力特性產(chǎn)生重大影響,從而影響實(shí)際的水彈道以及對航行體的控制[1-2].

    要維持理想的流動體制,對超空泡外形的控制就非常必要,目前國內(nèi)公開的文獻(xiàn)中主要是通過改變空化器線性和外形參數(shù)來控制空泡形態(tài)[3-5],Savchenko Yu.N對錐形空化器和Σ形空化器開展了研究,給出了它們的阻力特性與形狀參數(shù)和空化數(shù)的函數(shù)關(guān)系式,同時還對Σ形空化器的穩(wěn)定性進(jìn)行了深入的研究分析[6];或者通過采用非零攻角的空化器來實(shí)現(xiàn)非對稱超空泡流動模式,從而達(dá)到控制水下高速航行體運(yùn)動姿態(tài)和受力的目的,Waid[7]和Semenneko[8]對帶攻角的空化器形成的空泡形態(tài)進(jìn)行了研究,黃海龍等[9]通過數(shù)值模擬的方法詳細(xì)分析了不同攻角的圓盤空化器生成自然超空泡的形態(tài)特征及升阻力特性.

    當(dāng)航行體運(yùn)動速度增加并接近蒸汽空泡狀態(tài)時,可以通過改變空化器阻力系數(shù)cx0來控制超空泡參數(shù),隨著空化器阻力系數(shù)的變化,超空泡尺寸也會相應(yīng)改變.

    本文對不同中心錐角的變阻力空化器在空泡外形的控制能力進(jìn)行了數(shù)值研究,比較了不同中心錐角空化器的阻力系數(shù)和空泡尺寸響應(yīng)速度;并提出了一種新的空泡控制方案,數(shù)值模擬結(jié)果也證明了該方案在阻力控制方面的靈敏性和可操作性.

    1 控制方程與數(shù)值模型

    1.1 控制方程

    利用 FLUENT6.2中的 Mixture多相流模型[10],對混合相的連續(xù)性方程和動量方程進(jìn)行求解.

    混合相連續(xù)性方程:

    式中:xi為笛卡爾坐標(biāo)向量,ρm為混合介質(zhì)密度:

    式中:αk、ρk分別表示第k相的體積分?jǐn)?shù)和密度,ui

    為混合物的質(zhì)量平均速度:

    式中:uk為第k相的速度矢量.

    混合相動量方程:

    式中:p為壓力,μt為湍流粘性系數(shù),μm為混合粘性系數(shù),其表達(dá)式如下:

    式中:μk為第k相的粘度.

    湍流模型采用標(biāo)準(zhǔn)k-ε兩方程模型,近壁區(qū)域的流動采用標(biāo)準(zhǔn)壁面函數(shù)法處理.

    1.2 自然空化模型

    FLUENT6.2中的自然空化模型[11]模擬的空泡流包括水和水蒸汽介質(zhì),空泡流建模采用混合均質(zhì)流理論,整個混合物允許對流,蒸汽相傳輸方程如下:

    式中:f為蒸汽相質(zhì)量分?jǐn)?shù),Re和Rc是蒸汽產(chǎn)生和潰滅項(xiàng):

    式中:下標(biāo)l和v分別表示液相和汽相,σ為液體的表面張力,Ce和Cc為經(jīng)驗(yàn)常數(shù),其值分別為0.02和 0.01.

    1.3 動網(wǎng)格技術(shù)

    采用動態(tài)分層法對變形區(qū)域網(wǎng)格進(jìn)行更新,與動邊界相鄰的網(wǎng)格層(圖1中的層j)依據(jù)其網(wǎng)格高度h被切割或與下一網(wǎng)格層(圖1中的層i)合并,具體法則如下:

    hmin>(1+αs)hideal,j層網(wǎng)格被分割;hmin<αchideal,j層網(wǎng)格被合并(潰滅).

    其中,hmin為網(wǎng)格層j的最小網(wǎng)格高度,hideal為理想網(wǎng)格高度,αs、αc分別為網(wǎng)格層分割和潰滅因子.

    圖1 動態(tài)分層法Fig.1 Dynamic layering method

    2 數(shù)值計(jì)算

    2.1 傳統(tǒng)變阻力空化器

    圖2是文獻(xiàn)[12-13]中提到的變阻力空化器.空化器由帶尖銳前緣的外腔體1,以及安裝于腔體內(nèi)部,與1有著相同直徑Dn的中心錐形元件2組成.通過改變外體1和中心元件2的相對位置,可以使空化器在不改變直徑的情況下阻力發(fā)生變化.

    圖2 傳統(tǒng)變阻力空化器示意圖Fig.2 Scheme of the traditional variab le-drag cavitator

    試驗(yàn)[6]研究表明,空化器阻力系數(shù)cx和工作行程x的變化幅度均隨中心元件錐角2β的減小而增大,可以看出,中心錐角2β=60°時,在工作行程內(nèi)可得到阻力系數(shù)幾乎3倍的變化.

    由于以上試驗(yàn)是基于穩(wěn)態(tài)過程進(jìn)行的,并不能反映超空泡外形隨阻力變化的瞬態(tài)響應(yīng),本文展開了基于動網(wǎng)格技術(shù)的空泡流瞬態(tài)數(shù)值計(jì)算.圖3為中心元件錐角2β=90°時變阻力空化器附近的網(wǎng)格劃分.

    計(jì)算中設(shè)置入口速度為30 m/s,空化數(shù)為0.174.初始時刻,中心元件最左端與外體尖緣部分位于同一豎直平面(x=h),首先進(jìn)行自然空泡流的穩(wěn)態(tài)計(jì)算,并以此作為動態(tài)計(jì)算的初始流場,之后中心元件在12 s內(nèi)以uc=0.1m/s的速度往前移動,其中空泡面定義為蒸汽相體積分?jǐn)?shù)為0.5的等值面.

    圖3 90°錐角變阻力空化器附近網(wǎng)格Fig.3 Grid around the variable-drag cavitator,2β =90°

    2.2 變阻力空化器瞬態(tài)空泡流模擬結(jié)果分析

    2.2.1 不同中心錐角空化器的空泡流計(jì)算結(jié)果

    圖4給出了2β=60°、90°的變阻力空化器在不同工作行程x處的空泡外形.可以看出:1)雖然空化器一直向前運(yùn)動,但在中心元件完全露出來之前(0≤x≤h),空泡長度幾乎沒有變化;2)隨著中心元件的繼續(xù)運(yùn)動,其尖緣部分開始空化,空泡長度也逐漸減小,達(dá)到穩(wěn)定后不再變化;3)中心錐角越大,穩(wěn)定狀態(tài)所對應(yīng)的工作行程越大:2β=90°時,空泡長度在x=-0.8Dn之后便不再變化,而且由于空化器阻力在x=-0.3Dn處的突然降低,在原來較大的空泡與新生成的較小空泡之間形成了明顯的階梯式邊界(見圖4(b)),這與文獻(xiàn)[6]中觀察到的空泡形態(tài)相符.

    圖4 不同工作行程處的空泡外形Fig.4 Cavity shapes at different working stokes

    圖5給出了60°、90°和120°中心錐角空化器的阻力系數(shù)計(jì)算結(jié)果(其中阻力系數(shù)Cd=F/0.5ρv2S,F(xiàn)為阻力,S為空化器的最大橫截面積).圖5(a)為90°中心錐角空化器的阻力系數(shù)計(jì)算值與試驗(yàn)數(shù)據(jù)的對比,由曲線可以看出,模擬值與試驗(yàn)值變化趨勢大體一致,與模擬結(jié)果相比,穩(wěn)態(tài)試驗(yàn)中的阻力隨工作行程的增加,其增長速度更為緩慢.由圖5(b)可以看出,中心元件的運(yùn)動使空化器的阻力系數(shù)發(fā)生了明顯變化,其變化的幅度隨中心錐角的減小而增大,其中60°錐角的變阻力空化器在工作行程內(nèi)可得到阻力系數(shù)約2.4倍的變化,這與試驗(yàn)數(shù)據(jù)是相吻合的;阻力系數(shù)在初始時刻最大,其最大值與圓盤空化器的阻力系數(shù)相接近;隨著中心錐角的減小,空化器對應(yīng)的阻力系數(shù)下限值也有所降低.

    圖5 阻力系數(shù)隨工作行程的變化曲線Fig.5 Curves of drag coefficients of the cavitator on working stokes

    圖6為無量綱空泡尺寸的變化曲線,由圖6知,空泡長度Lc和最大直徑Dc的變化趨勢基本一致:即在x=0.0Dn處空泡尺寸由最大值開始減小,x=-0.80Dn處,空泡長度及最大直徑降至最低,之后不再變化;Lc、Dc的變化幅度隨錐角的減小而增大:2β=60°時,工作行程內(nèi)可得到空泡長度1.15倍和最大直徑1.29倍的變化.不同中心錐角的變阻力空化器可獲得的最大空泡尺寸與中心元件形狀無關(guān),當(dāng)工作行程值x接近中心元件高度h時(見圖2(a)),空泡尺寸達(dá)到最大.

    結(jié)合阻力系數(shù)和空泡尺寸的變化曲線可以看出,兩者變化并不同步,空泡長度及最大直徑變化要滯后于阻力的變化,初始時刻(x=h),阻力系數(shù)和空泡尺寸取得最大值,之后隨著工作行程值x的逐漸減小,阻力系數(shù)開始減小,在x=-0.50Dn達(dá)到最小,而空泡長度及最大直徑在x=0.0Dn處才開始變化,在x=-0.80Dn處降至最低.

    圖6 不同錐角變阻力空化器的空泡尺寸變化曲線Fig.6 Curves of variance of cavity size for the variable-drag cavitatorswith different cone angles

    2.2.2 中心元件速度對模擬結(jié)果的影響

    為研究中心元件的運(yùn)動速度對空泡流場結(jié)構(gòu)的影響,本文選取中心錐角2β=60°的變阻力空化器,對其中心元件速度uc=1 m/s時的瞬態(tài)空泡流進(jìn)行了模擬計(jì)算.

    圖7為不同uc下的阻力系數(shù)Cd曲線,由圖知2種速度下的阻力系數(shù)變化趨勢基本一致:即隨著工作行程值x的減小,阻力系數(shù)逐漸下降,降至最低值后不再變化,中心元件速度為1 m/s時,計(jì)算的阻力系數(shù)值要略大些.

    圖7 阻力系數(shù)隨工作行程變化曲線,2β=60°Fig.7 Curves of drag coefficient on working stoke,2β =60°

    空泡尺寸隨工作行程變化曲線如圖8,隨著工作行程值x的減小,空泡長度Lc和最大直徑Dc并沒有馬上減小,而是經(jīng)過一段時間的遲滯才開始變化;由于與速度為0.1 m/s的計(jì)算結(jié)果相比,中心元件速度為1 m/s時得到的阻力系數(shù)值要稍微大些,故該速度下計(jì)算的空泡尺寸也有小幅增加,由于遲滯效應(yīng),1 m/s速度下計(jì)算的空泡長度在x=0.466Dn時開始增加,直到工作行程值x=-0.80Dn,2種中心元件運(yùn)動速度下的空泡長度不再變化;而空泡最大直徑的增加與阻力系數(shù)同步.

    圖8 不同中心元件速度的空泡尺寸隨工作行程變化曲線,2β=60°Fig.8 Curves of cavity size on working stoke at different velocity magnitudes of central element,2β =60°

    結(jié)合圖7、8可以看出,中心元件運(yùn)動速度地大小對變阻力空化器的空泡流模擬結(jié)果有一定影響,隨著速度的增大,空泡尺寸及阻力系數(shù)值均有所增加,其中空泡長度的增加幅度最大達(dá)到了10%,而空泡最大直徑和阻力系數(shù)受速度影響較小.

    2.3 改進(jìn)的可變側(cè)向力空化器

    2.3.1 控制方案與內(nèi)部結(jié)構(gòu)

    上述變阻力空化器雖然可以比較快速的改變航行體阻力和空泡尺寸,然而要想獲得良好的控制效果,需要其控制機(jī)構(gòu)能夠輕易提供數(shù)十萬牛頓的力,這一控制方法很難實(shí)施.考慮到空化器側(cè)面的壓力較低,而且通過增加空化器有效面積同樣可以改變航行體的受力,本文提出了如圖9所示的可變側(cè)向力空化器的控制方案:在錐形空化器內(nèi)部安裝一外徑小于Dn的可動元件3,通過可動元件3的橫向運(yùn)動來控制元件2的位置,元件2、3均設(shè)計(jì)成4瓣,控制元件2內(nèi)部共有8對弧形鋼片4,其中每瓣的兩端均有一對弧形鋼片,可動元件3由錐形可動部分和圓柱操縱桿連接而成,4瓣圓柱操縱桿分別控制4瓣圓錐部分的橫向位移,從而改變控制元件的周向位移hn值,使前端空化器1的有效工作面積變化,進(jìn)而改變空泡尺寸;當(dāng)可動元件任意兩側(cè)位移不一致時,控制元件具有不對稱性,前端錐形空化器產(chǎn)生不對稱受力,使空泡形態(tài)發(fā)生變化,可以有效的控制航行體運(yùn)動姿態(tài),有關(guān)可變側(cè)向力空化器的具體內(nèi)容可參考文獻(xiàn)[14].

    圖9 可變側(cè)向力空化器的控制方案Fig.9 Control scheme of the variable-lateral force cavitator

    2.3.2 空泡流模擬結(jié)果

    圖10、11給出了45°錐角空化器的控制元件2上半周部分突出1/10Dn時的空泡流計(jì)算結(jié)果,由于元件2的不對稱性,在空化器上產(chǎn)生豎直向下的升力,導(dǎo)致空泡上下截面不對稱(如圖10所示).3種空化數(shù)下得到的空化器阻力系數(shù)分別為0.620、0.553和0.544,相比于圓錐空化器,可變側(cè)向力空化器在相同空化數(shù)下的阻力系數(shù)增加了近60%,證實(shí)了其在阻力控制方面的靈敏性.

    圖10 可變側(cè)向力空化器的空泡形態(tài)Fig.10 Cavity shapes of the variable-lateral force cavitator

    空化數(shù)為0.190時的空化器前端縱剖面的壓力系數(shù)分布如圖11所示,可以看出,凸出的元件2導(dǎo)致了不對稱的壓力分布,從而在空化器上產(chǎn)生了相當(dāng)于30%阻力的升力,升力系數(shù)達(dá)到0.177;此外壓力系數(shù)等值線顯示,控制元件2外表面的壓力非常低,這保證了元件的可操作性.

    圖11 可變側(cè)向力空化器前端壓力系數(shù)等值線,σ=0.190Fig.11 Contours of pressure coefficient in front of the variable-lateral force cavitator,σ =0.190

    3 結(jié)論

    本文利用FLUENT6.2的自然空化模型和動網(wǎng)格技術(shù)對水下超空泡航行體的瞬態(tài)空泡流進(jìn)行了數(shù)值模擬,得到了以下幾點(diǎn)結(jié)論:

    1)傳統(tǒng)變阻力空化器的阻力響應(yīng)速度較快,空泡外形的變化則要滯后于阻力變化;隨著工作行程值的減小,不同錐角的變阻力空化器的阻力系數(shù)和空泡尺寸也逐漸減小;

    2)中心錐角越小,對應(yīng)的空化器工作行程及阻力系數(shù)(空泡尺寸)的變化幅度也就越大,2β=60°時在工作行程內(nèi)可得到阻力系數(shù)幾乎3倍的變化;

    3)中心元件運(yùn)動速度對模擬結(jié)果也有一定影響,速度值從0.1m/s增至1m/s時,空泡尺寸有小幅增加;

    4)數(shù)值模擬結(jié)果顯示可變側(cè)向力空化器對阻力和升力的控制非常靈敏.

    [1]MAY A.Water entry and cavity-running behavior ofmissles[R].AD-A020259,1975.

    [2]WOLFE wP.Drag of revolution in supercavitating flow[J].Journal of Fluids Engineering,1989,112:23-30.

    [3]賈力平,張嘉鐘,于開平,等.空化器線形與超空泡減阻效果關(guān)系研究[J].船舶工程,2006,28(2):20-23.JIA Lingping,ZHANG Jiazhong,YU Kaiping,et al.Research on the relation of cavitator profile and supercavity dragreduction[J].ShipEngineering,2006,28(2):20-23.

    [4]顧建農(nóng),張志宏,高永琪,等.充氣頭型對超空泡軸對稱體阻力特性影響的試驗(yàn)研究[J].兵工學(xué)報(bào),2004,25(6):766-769.GU Jiangnong,ZHANG Zhihong,GAO Yongqi,et al.An experimental study about the influence of ventilated headforms on the drag characters of a super cavitation axisymmetric body[J].Acta Armamentarll,2004,25(6):766-769.

    [5]賈力平,王聰,于開平,等.空化器參數(shù)對通氣超空泡形態(tài)影響的實(shí)驗(yàn)研究[J].工程力學(xué),2007,24(3):159-164.JIA Liping,WANG Cong,YU Kaiping,etal.Experimental investigation of cavitator parameters effecting on ventilated supercavity shape[J].Engineering Mechanics,2007,24(3):159-164.

    [6]SAVCHENKO Y N.Control of supercavitation flowand stability of supercavitating motion of bodies[C]//The RTO AVT Lectures Series on"Supercavitating Flows",RTO Lecture Series 005.Ottawa:St.Joseph,2002:313-341.

    [7]WAID R L.Cavity shapes for circular disks at angles of attack[R].California Institute of Technology.Report No.E-73.4,1975.

    [8]SEMENENKO V N.Artificial supercavitation physics and calculation[C]//RTO/AVT lecture series on"supercavitating flows",RTO Lecture Series 005,Ottawa:St.joseph,2002:195-227.

    [9]黃海龍,王聰,黃文虎,等.變攻角圓盤空化器生成自然超空泡數(shù)值模擬[J].北京理工大學(xué)學(xué)報(bào),2008,28(2):100-103.HUANG Hailong,WANG Cong,HUANG Wenhu,et al.Numerical simulation on natural supercavity based on diskcavitator with variable attack angle[J].Transactions of Beijing Institute of Technology,2008,28(2):100-103.

    [10]MANNINEN M,TAIVASSALO V,KALLIO S.On the mixturemodel for multiphase flow[M].Espoo:Technical Research Center of Finland,VTT Publications,1996:3-67.

    [11]SINGHALA K,ATHAVALEmM,LIH Y,etal.Mathematical basis and validation of the full cavitation model[J].Journal of Fluids Engineering,2002,124:617-624.

    [12]FU Huiping,LU Chuanjing,LIJie.Numercial research on drag reduction characteristics of supercavitating body of revolution[J].Journal of ShipMechanics,2004,8(3):1-8.

    [13]SAVCHENKO Y N,VLASENKO Y D,SEMENENKO V N.Experimental study of high-speed cavitation flows[J].Jof Fluid Mechanics Research,1998,26(3):103-111(in Russian).

    [14]哈爾濱工程大學(xué).可變側(cè)向力空化器[P].中國:20091007,2007.3,2011,4.21.

    猜你喜歡
    錐角空泡空化
    錐角比對雙錐藥型罩射流成型影響的數(shù)值模擬
    功率超聲作用下鋼液中空化泡尺寸的演變特性
    鋼鐵釩鈦(2023年5期)2023-11-17 08:48:34
    水下航行體雙空泡相互作用數(shù)值模擬研究
    高鐵箱梁預(yù)應(yīng)力夾片式錨具錐角的數(shù)值分析
    錐形避雷針避雷效果最優(yōu)錐角研究
    三維扭曲水翼空化現(xiàn)象CFD模擬
    不同運(yùn)動形式下水物相互作用空化數(shù)值模擬
    基于LPV的超空泡航行體H∞抗飽和控制
    基于CFD的對轉(zhuǎn)槳無空泡噪聲的仿真預(yù)報(bào)
    船海工程(2015年4期)2016-01-05 15:53:28
    SPH在水下高速物體空泡發(fā)展模擬中的應(yīng)用
    国产精品熟女久久久久浪| 日韩成人av中文字幕在线观看| 99久国产av精品国产电影| 午夜福利在线观看免费完整高清在| 国产亚洲最大av| 大陆偷拍与自拍| 能在线免费观看的黄片| 777米奇影视久久| 三级经典国产精品| 天天躁日日操中文字幕| 欧美三级亚洲精品| 久热久热在线精品观看| 最近2019中文字幕mv第一页| 汤姆久久久久久久影院中文字幕 | 69人妻影院| 在线观看av片永久免费下载| 成人美女网站在线观看视频| 免费看美女性在线毛片视频| av在线观看视频网站免费| 蜜桃久久精品国产亚洲av| 国产av国产精品国产| 日韩一本色道免费dvd| 中文字幕人妻熟人妻熟丝袜美| 高清av免费在线| 天天一区二区日本电影三级| 国产精品一区www在线观看| 免费电影在线观看免费观看| 中文资源天堂在线| 亚洲高清免费不卡视频| 日韩一本色道免费dvd| 精品久久久久久久久久久久久| 观看美女的网站| .国产精品久久| 国产视频内射| 亚洲欧洲日产国产| 亚洲av日韩在线播放| 极品少妇高潮喷水抽搐| 特大巨黑吊av在线直播| 嫩草影院精品99| 极品少妇高潮喷水抽搐| 国产免费视频播放在线视频 | 久久久精品94久久精品| 在线a可以看的网站| 日韩欧美国产在线观看| 亚洲最大成人av| av在线蜜桃| 波多野结衣巨乳人妻| 亚洲精品国产成人久久av| 精品熟女少妇av免费看| 亚洲av不卡在线观看| 69人妻影院| 汤姆久久久久久久影院中文字幕 | 婷婷色麻豆天堂久久| av在线观看视频网站免费| 七月丁香在线播放| 国产 亚洲一区二区三区 | 熟妇人妻久久中文字幕3abv| 日韩一区二区三区影片| 国产 亚洲一区二区三区 | 伦精品一区二区三区| 国产亚洲一区二区精品| 免费大片18禁| 亚洲激情五月婷婷啪啪| 美女被艹到高潮喷水动态| 最后的刺客免费高清国语| 日韩精品有码人妻一区| 国产成人精品久久久久久| 国产成人精品一,二区| 午夜福利视频1000在线观看| 久久精品久久久久久噜噜老黄| 欧美激情国产日韩精品一区| 亚洲精品中文字幕在线视频 | 亚洲精品成人av观看孕妇| 日本猛色少妇xxxxx猛交久久| 久久精品国产亚洲av涩爱| 国产一级毛片在线| 欧美bdsm另类| 国产精品麻豆人妻色哟哟久久 | 淫秽高清视频在线观看| 欧美+日韩+精品| 成人二区视频| 舔av片在线| 久久99蜜桃精品久久| 国产精品99久久久久久久久| 日韩 亚洲 欧美在线| 国产亚洲91精品色在线| 熟妇人妻久久中文字幕3abv| 免费黄色在线免费观看| 国产视频首页在线观看| 天天躁夜夜躁狠狠久久av| 午夜久久久久精精品| ponron亚洲| 精品欧美国产一区二区三| 一个人观看的视频www高清免费观看| 精品人妻熟女av久视频| 日日撸夜夜添| 大香蕉久久网| 亚洲在线观看片| 男人和女人高潮做爰伦理| 亚洲精品色激情综合| 日韩一本色道免费dvd| 不卡视频在线观看欧美| 熟妇人妻不卡中文字幕| 听说在线观看完整版免费高清| 在线观看av片永久免费下载| 国产成人一区二区在线| 人妻系列 视频| 一级毛片我不卡| 久久韩国三级中文字幕| 日韩大片免费观看网站| 国产毛片a区久久久久| 欧美xxxx性猛交bbbb| 日韩在线高清观看一区二区三区| 久久精品夜夜夜夜夜久久蜜豆| 国产色爽女视频免费观看| 午夜福利视频1000在线观看| 成人特级av手机在线观看| 中文字幕av成人在线电影| 国产成人freesex在线| 18禁动态无遮挡网站| 卡戴珊不雅视频在线播放| 日韩一本色道免费dvd| 中文资源天堂在线| 亚洲精品乱码久久久v下载方式| 久久97久久精品| 国产成人福利小说| av免费观看日本| 大香蕉久久网| 国产在线男女| 国产精品蜜桃在线观看| 午夜福利在线观看吧| 男女国产视频网站| 色综合亚洲欧美另类图片| 久久99精品国语久久久| 国产男女超爽视频在线观看| 直男gayav资源| 午夜免费激情av| av在线天堂中文字幕| 免费看a级黄色片| 看免费成人av毛片| videos熟女内射| 男人爽女人下面视频在线观看| 国产亚洲午夜精品一区二区久久 | 欧美日韩亚洲高清精品| 国精品久久久久久国模美| 夫妻性生交免费视频一级片| 国产在线一区二区三区精| 亚洲精品国产成人久久av| 国产女主播在线喷水免费视频网站 | 久久久久久久久久人人人人人人| 免费看日本二区| 亚洲av免费高清在线观看| or卡值多少钱| 久久久久久久久大av| 亚洲欧美日韩无卡精品| 看十八女毛片水多多多| 人人妻人人看人人澡| 成人二区视频| 亚洲av一区综合| 国产一区亚洲一区在线观看| 婷婷色麻豆天堂久久| 少妇高潮的动态图| 国产伦理片在线播放av一区| 联通29元200g的流量卡| 一级av片app| 日本欧美国产在线视频| 卡戴珊不雅视频在线播放| 美女黄网站色视频| 街头女战士在线观看网站| 国产精品久久久久久久电影| 免费看光身美女| 人妻一区二区av| 日韩在线高清观看一区二区三区| av专区在线播放| 欧美日韩亚洲高清精品| 国产一级毛片在线| 麻豆乱淫一区二区| 亚洲av免费高清在线观看| 国产成人a∨麻豆精品| 亚洲国产欧美人成| 亚洲av国产av综合av卡| 插逼视频在线观看| 亚洲成色77777| 久久精品人妻少妇| 国产精品一区二区三区四区免费观看| 精品人妻一区二区三区麻豆| 蜜桃久久精品国产亚洲av| 熟妇人妻不卡中文字幕| 久99久视频精品免费| 国国产精品蜜臀av免费| 亚洲怡红院男人天堂| 亚洲真实伦在线观看| 大陆偷拍与自拍| 国产成人午夜福利电影在线观看| 日本黄色片子视频| 插阴视频在线观看视频| 狂野欧美激情性xxxx在线观看| 我的女老师完整版在线观看| 欧美日韩精品成人综合77777| 色网站视频免费| 国产成人一区二区在线| 一级爰片在线观看| 亚洲精品一区蜜桃| 免费黄色在线免费观看| 日本一二三区视频观看| 啦啦啦中文免费视频观看日本| 有码 亚洲区| 亚洲欧洲国产日韩| 少妇人妻一区二区三区视频| 午夜福利高清视频| 寂寞人妻少妇视频99o| 欧美日韩综合久久久久久| 久久久精品免费免费高清| 欧美高清性xxxxhd video| 国产亚洲一区二区精品| a级毛色黄片| 五月天丁香电影| 久久韩国三级中文字幕| 日日干狠狠操夜夜爽| 嫩草影院新地址| 久久久精品94久久精品| 少妇的逼好多水| 欧美精品国产亚洲| 久热久热在线精品观看| 久久精品夜夜夜夜夜久久蜜豆| 亚洲美女视频黄频| 久久热精品热| 免费观看av网站的网址| 特级一级黄色大片| 真实男女啪啪啪动态图| 成年av动漫网址| 女的被弄到高潮叫床怎么办| 亚洲人成网站在线播| 中文字幕免费在线视频6| 嘟嘟电影网在线观看| 国产v大片淫在线免费观看| 韩国高清视频一区二区三区| 麻豆av噜噜一区二区三区| 韩国av在线不卡| 简卡轻食公司| 久久国内精品自在自线图片| 日本一二三区视频观看| 高清欧美精品videossex| 春色校园在线视频观看| 男人和女人高潮做爰伦理| 免费观看的影片在线观看| 亚洲av免费高清在线观看| 简卡轻食公司| 久久久久性生活片| 建设人人有责人人尽责人人享有的 | av免费观看日本| 日本与韩国留学比较| 免费观看a级毛片全部| 亚洲怡红院男人天堂| 菩萨蛮人人尽说江南好唐韦庄| 国产伦精品一区二区三区视频9| 亚洲成人中文字幕在线播放| 97人妻精品一区二区三区麻豆| 国产黄频视频在线观看| 日日啪夜夜撸| eeuss影院久久| 99久久精品热视频| 国产成人91sexporn| 久久久久久伊人网av| 亚洲熟妇中文字幕五十中出| 美女黄网站色视频| 国内揄拍国产精品人妻在线| 久久精品夜夜夜夜夜久久蜜豆| 免费av毛片视频| 99热网站在线观看| 亚洲av.av天堂| 中文精品一卡2卡3卡4更新| 嫩草影院新地址| 最近2019中文字幕mv第一页| 性插视频无遮挡在线免费观看| 国产高清国产精品国产三级 | 午夜爱爱视频在线播放| 亚洲精品国产av成人精品| 人妻少妇偷人精品九色| 色视频www国产| 国产91av在线免费观看| 赤兔流量卡办理| 岛国毛片在线播放| 国产免费福利视频在线观看| 久久久午夜欧美精品| 五月伊人婷婷丁香| av又黄又爽大尺度在线免费看| 国产精品国产三级国产av玫瑰| 国产男人的电影天堂91| 久久久久久久午夜电影| 超碰av人人做人人爽久久| 丰满少妇做爰视频| 免费黄色在线免费观看| 成人午夜高清在线视频| 国产精品一区二区三区四区免费观看| 中文字幕av成人在线电影| 一级毛片我不卡| 国产大屁股一区二区在线视频| 免费观看无遮挡的男女| 日韩中字成人| 男人舔奶头视频| 热99在线观看视频| 我的女老师完整版在线观看| 国产成人91sexporn| 草草在线视频免费看| 我的女老师完整版在线观看| 我要看日韩黄色一级片| 亚洲婷婷狠狠爱综合网| 亚洲精品日本国产第一区| 国产69精品久久久久777片| 国产在视频线在精品| 日韩在线高清观看一区二区三区| 久久久久久久久久成人| 亚洲第一区二区三区不卡| 2021天堂中文幕一二区在线观| av免费在线看不卡| 国产精品日韩av在线免费观看| 日韩av在线大香蕉| 在现免费观看毛片| 国产69精品久久久久777片| 亚洲不卡免费看| 又黄又爽又刺激的免费视频.| 人人妻人人看人人澡| 插逼视频在线观看| 久久久午夜欧美精品| xxx大片免费视频| 国产一区二区在线观看日韩| 国产欧美另类精品又又久久亚洲欧美| 97在线视频观看| 69人妻影院| 国产成人福利小说| 欧美性猛交╳xxx乱大交人| 哪个播放器可以免费观看大片| 亚洲伊人久久精品综合| 国产亚洲午夜精品一区二区久久 | 99re6热这里在线精品视频| 欧美+日韩+精品| 免费观看在线日韩| 国产亚洲5aaaaa淫片| 一级毛片 在线播放| or卡值多少钱| 欧美激情在线99| 国产精品熟女久久久久浪| 国产中年淑女户外野战色| 亚洲成人一二三区av| 精品久久久噜噜| 一个人看的www免费观看视频| 大话2 男鬼变身卡| 国产亚洲5aaaaa淫片| 精品国产三级普通话版| 国产亚洲最大av| 国产精品久久视频播放| 久久精品夜色国产| 久久久久久久国产电影| 波多野结衣巨乳人妻| 在线天堂最新版资源| 亚洲av中文av极速乱| 亚洲熟妇中文字幕五十中出| 91久久精品国产一区二区成人| 一级黄片播放器| 精品99又大又爽又粗少妇毛片| 午夜久久久久精精品| 亚洲成色77777| 精品人妻视频免费看| 一边亲一边摸免费视频| 禁无遮挡网站| 亚洲av电影不卡..在线观看| 搞女人的毛片| 22中文网久久字幕| 黑人高潮一二区| av国产久精品久网站免费入址| 狠狠精品人妻久久久久久综合| 亚洲国产日韩欧美精品在线观看| 国产精品.久久久| 亚洲熟女精品中文字幕| 卡戴珊不雅视频在线播放| 国产乱来视频区| 国产精品麻豆人妻色哟哟久久 | 亚洲欧美成人综合另类久久久| 国产亚洲午夜精品一区二区久久 | 国产精品99久久久久久久久| 91午夜精品亚洲一区二区三区| 色综合亚洲欧美另类图片| 日本免费a在线| 国产黄色视频一区二区在线观看| 欧美高清性xxxxhd video| 97超碰精品成人国产| 亚洲av日韩在线播放| 在线免费观看的www视频| 国语对白做爰xxxⅹ性视频网站| 国产麻豆成人av免费视频| 伊人久久国产一区二区| 一个人看的www免费观看视频| 草草在线视频免费看| 在线观看免费高清a一片| 中文字幕免费在线视频6| 国产精品一区二区三区四区免费观看| 日韩av在线免费看完整版不卡| 久久人人爽人人爽人人片va| 婷婷色综合www| 99热这里只有是精品50| 免费观看av网站的网址| 国产精品国产三级国产av玫瑰| 美女cb高潮喷水在线观看| 中文在线观看免费www的网站| 白带黄色成豆腐渣| 亚洲伊人久久精品综合| 成人无遮挡网站| 极品少妇高潮喷水抽搐| 你懂的网址亚洲精品在线观看| 日韩大片免费观看网站| 日本一二三区视频观看| 国产色婷婷99| 在线观看av片永久免费下载| 激情五月婷婷亚洲| 麻豆av噜噜一区二区三区| 国产探花极品一区二区| 精品国产三级普通话版| 久久97久久精品| 色尼玛亚洲综合影院| 国产熟女欧美一区二区| 亚洲精品色激情综合| 男插女下体视频免费在线播放| 日韩欧美国产在线观看| 大陆偷拍与自拍| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产一级毛片七仙女欲春2| 国产精品一区二区三区四区久久| 老司机影院成人| 联通29元200g的流量卡| ponron亚洲| 淫秽高清视频在线观看| 国产69精品久久久久777片| 免费看光身美女| 国产成人免费观看mmmm| 成人无遮挡网站| 国产成人精品一,二区| 国产精品综合久久久久久久免费| 亚洲欧洲日产国产| 麻豆精品久久久久久蜜桃| 亚洲久久久久久中文字幕| 特级一级黄色大片| av在线蜜桃| 国产伦精品一区二区三区视频9| 国产精品久久视频播放| 99热这里只有是精品50| 欧美精品一区二区大全| 激情 狠狠 欧美| 超碰97精品在线观看| 黄色一级大片看看| 黄片wwwwww| 18禁裸乳无遮挡免费网站照片| 黄色欧美视频在线观看| 国产免费福利视频在线观看| 久久久久久久午夜电影| 久久精品国产亚洲网站| 不卡视频在线观看欧美| 免费播放大片免费观看视频在线观看| 插阴视频在线观看视频| 亚洲国产最新在线播放| 国产成人精品婷婷| 51国产日韩欧美| 中文资源天堂在线| 26uuu在线亚洲综合色| 久久久久九九精品影院| 伊人久久国产一区二区| 永久免费av网站大全| 日韩强制内射视频| 蜜桃久久精品国产亚洲av| 国产一级毛片七仙女欲春2| 国产高清国产精品国产三级 | 精品亚洲乱码少妇综合久久| 国产午夜福利久久久久久| 淫秽高清视频在线观看| 欧美高清性xxxxhd video| 国产女主播在线喷水免费视频网站 | 嫩草影院精品99| 晚上一个人看的免费电影| 国产熟女欧美一区二区| 青春草亚洲视频在线观看| 色尼玛亚洲综合影院| 免费观看在线日韩| 一夜夜www| 免费av观看视频| 国产黄色免费在线视频| 欧美激情在线99| 中文字幕久久专区| 亚洲激情五月婷婷啪啪| 综合色丁香网| 真实男女啪啪啪动态图| 国内少妇人妻偷人精品xxx网站| 欧美精品一区二区大全| 舔av片在线| 伊人久久国产一区二区| 亚洲精品亚洲一区二区| 国精品久久久久久国模美| 老师上课跳d突然被开到最大视频| 91久久精品电影网| 免费观看性生交大片5| 免费少妇av软件| 国产三级在线视频| 久久午夜福利片| 国产免费又黄又爽又色| 中文字幕av成人在线电影| 日韩成人av中文字幕在线观看| 日韩欧美国产在线观看| 亚洲av国产av综合av卡| 男女视频在线观看网站免费| 青春草视频在线免费观看| 久久国产乱子免费精品| 亚洲四区av| 大香蕉久久网| 国产精品三级大全| 特大巨黑吊av在线直播| 天堂中文最新版在线下载 | 大香蕉久久网| 国产成人精品婷婷| 欧美 日韩 精品 国产| 麻豆久久精品国产亚洲av| 免费看美女性在线毛片视频| videos熟女内射| 亚洲欧美日韩无卡精品| 亚洲精品影视一区二区三区av| 中文字幕久久专区| 亚洲熟女精品中文字幕| 一本一本综合久久| 精品少妇黑人巨大在线播放| 国产精品人妻久久久影院| 又爽又黄无遮挡网站| 三级国产精品欧美在线观看| 国产免费视频播放在线视频 | 一级毛片 在线播放| 三级国产精品片| 身体一侧抽搐| 两个人视频免费观看高清| 国产色婷婷99| 久久99热这里只有精品18| 日韩一区二区三区影片| 高清午夜精品一区二区三区| 欧美激情在线99| 中国美白少妇内射xxxbb| 免费不卡的大黄色大毛片视频在线观看 | 国产久久久一区二区三区| 日韩人妻高清精品专区| 久久久久九九精品影院| 欧美人与善性xxx| 久久久久久国产a免费观看| 亚洲美女搞黄在线观看| 国产黄频视频在线观看| 欧美xxxx性猛交bbbb| 亚洲欧美一区二区三区黑人 | 真实男女啪啪啪动态图| 日韩中字成人| 久久久久性生活片| 国产免费一级a男人的天堂| 成人无遮挡网站| 一个人看视频在线观看www免费| 丰满乱子伦码专区| 免费少妇av软件| 亚洲国产av新网站| 一级av片app| 国产三级在线视频| 国产黄色视频一区二区在线观看| 免费观看av网站的网址| 嫩草影院入口| www.av在线官网国产| 亚洲av一区综合| 80岁老熟妇乱子伦牲交| 亚洲av成人精品一二三区| 国产精品伦人一区二区| 亚洲乱码一区二区免费版| 国产精品久久久久久精品电影小说 | 别揉我奶头 嗯啊视频| 成人毛片a级毛片在线播放| 日本猛色少妇xxxxx猛交久久| 少妇人妻精品综合一区二区| 欧美日本视频| 日产精品乱码卡一卡2卡三| 搡老妇女老女人老熟妇| 少妇的逼水好多| av免费观看日本| 人妻少妇偷人精品九色| av在线蜜桃| 久久久色成人| 身体一侧抽搐| 日韩亚洲欧美综合| 国产精品一区二区性色av| 日本欧美国产在线视频| 肉色欧美久久久久久久蜜桃 | 久久精品综合一区二区三区| 国产色婷婷99| 亚洲国产欧美在线一区| 亚洲人成网站高清观看| 最近最新中文字幕免费大全7| 日日啪夜夜爽| 菩萨蛮人人尽说江南好唐韦庄| 久久久久性生活片| 国产av在哪里看| 国产精品一二三区在线看| 成人毛片a级毛片在线播放| 成人性生交大片免费视频hd| 成人毛片a级毛片在线播放| 国产精品熟女久久久久浪| 国产av码专区亚洲av| 国产精品久久视频播放| 亚洲色图av天堂| 国产精品久久久久久av不卡| 十八禁网站网址无遮挡 | 又粗又硬又长又爽又黄的视频| 中文字幕人妻熟人妻熟丝袜美| 99热网站在线观看| 伊人久久国产一区二区| 国产成人免费观看mmmm| 乱码一卡2卡4卡精品| 我的女老师完整版在线观看|