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

    航行體沾濕區(qū)域?qū)张菸擦鹘Y(jié)構(gòu)的影響

    2019-11-06 01:26:02王威王聰李聰慧宋武超
    兵工學(xué)報 2019年10期
    關(guān)鍵詞:渦量尾流空泡

    王威,王聰,李聰慧,宋武超

    (哈爾濱工業(yè)大學(xué) 航天學(xué)院,黑龍江 哈爾濱 150001)

    0 引言

    水下武器借助通氣超空泡的包裹,可以大幅度地減小水下阻力、提高航行速度。超空泡多相流動特性研究始終是各國學(xué)者關(guān)注的熱點問題,其中空泡尾流特性是通氣超空泡多相流動研究的重要方面,空泡尾部的雙渦管泄氣現(xiàn)象多在低弗勞德數(shù)Fr條件下產(chǎn)生,水洞實驗也易于觀察,因此針對雙渦管尾流特性的研究備受各國學(xué)者青睞。

    從工程應(yīng)用的角度看,帶有后體的超空泡航行體多相流特性也是需要考慮的問題,袁緒龍等[7]和陳偉政等[8]通過水洞實驗分析了沾濕區(qū)域及重力對空泡形態(tài)的非對稱性影響。航行體后體與流動介質(zhì)相接觸[9],沾濕區(qū)域是維持航行體受力平衡的關(guān)鍵。沾濕區(qū)域從空泡下部穿出,必然會影響空泡多相尾流結(jié)構(gòu),這與無后體的自由閉合空泡尾流存在較大差異。Yu等[10]對航行體后體沾濕區(qū)域的升力進行了研究,分析了升力突變的原因。Wang等[11]研究了超空泡航行體的尾流結(jié)構(gòu),采用定常計算方法定性分析了空泡尾部流場壓力和渦量分布特性,指出了沾濕區(qū)域?qū)u量方向的影響,并給出了渦管內(nèi)部流動參數(shù)的變化規(guī)律。

    從相關(guān)文獻報道可以看出,對圓盤空化器模型的研究多側(cè)重于空泡尾流結(jié)構(gòu),對帶有后體的航行體模型研究多側(cè)重于流體動力,但研究后體沾濕區(qū)域?qū)叫畜w空泡尾流結(jié)構(gòu)的影響并不多見。數(shù)值模擬方法是研究超空泡尾流結(jié)構(gòu)的重要方法,其優(yōu)點是易于捕捉實驗方法較難獲得的流場參數(shù)變化特性,分析流場參數(shù)非定常變化。

    本文利用數(shù)值模擬方法,在文獻[11]的基礎(chǔ)上,采用非定常計算方法,進一步定量分析低Fr條件下航行體攻角逐漸變化時沾濕區(qū)域形成過程中空泡尾流場壓力及渦量分布的演化特性,揭示沾濕區(qū)域?qū)张菸擦鹘Y(jié)構(gòu)的影響機理,給出沾濕區(qū)域比例與空泡尾流結(jié)構(gòu)的轉(zhuǎn)變關(guān)系。

    1 數(shù)學(xué)模型及驗證

    1.1 航行體模型

    根據(jù)低Fr條件(5

    航行體模型如圖1所示,該模型由圓盤空化器、圓錐段、圓柱段、圓柱形尾噴管組成。航行體模型圓柱段直徑為D,航行體全長L=12D,錐段長度為L/2,圓盤空化器直徑Dn=0.3D,航行體質(zhì)心位置在xc=6.5D處,坐標(biāo)系原點建立在航行體質(zhì)心。

    圖1 計算模型Fig.1 Computation model

    1.2 計算流域邊界條件

    計算流域的邊界如圖2所示(圖中p∞為遠場邊界壓力),計算域為圓柱形,直徑為20D,長度為48D,重力方向沿y軸負向,計算流域采用六面體結(jié)構(gòu)化網(wǎng)格劃分。航行體表面采用壁面邊界,在空化器后部設(shè)有通氣孔,采用質(zhì)量入口條件;流域入口和流域圓柱面邊界都采用速度入口條件,流域出口采用水深5 m壓力出口條件。

    圖2 計算域及邊界條件Fig.2 Computational domain and boundary condition

    Zhou等[12]研究了無限流域尺度對數(shù)值模擬結(jié)果的影響,指出p∞的等值線應(yīng)在計算流域內(nèi)部。本文利用p∞的等值面判斷三維無限流域數(shù)值模擬對流域的最小尺度要求。從圖2可以發(fā)現(xiàn),p∞的等值面并未與計算流域邊界交叉,且在計算流域內(nèi)部,表明本文的流域邊界尺度用于模擬無限流域情況下的通氣空泡流動是合理的。

    比較不同網(wǎng)格密度條件下航行體全濕狀態(tài)時的流體動力系數(shù)隨網(wǎng)格數(shù)量Ne的變化,結(jié)果如圖3所示。圖3中,CD為阻力系數(shù),CL為升力系數(shù)。

    圖3 不同網(wǎng)格密度下航行體的流體動力系數(shù)Fig.3 Hydrodynamic coefficients of vehicle under different grid densities

    從圖3中可以看出:在0.6×106~0.9×106網(wǎng)格密度條件下,航行體的阻力系數(shù)CD和升力系數(shù)CL均隨著網(wǎng)格數(shù)量而變化;在0.9×106~1.2×106的網(wǎng)格密度條件下,航行體的阻力和升力系數(shù)不再隨網(wǎng)格數(shù)量變化。采用0.9×106以上網(wǎng)格數(shù)量進行數(shù)值模擬時,計算結(jié)果不會因為網(wǎng)格疏密程度而造成計算誤差,綜合考慮數(shù)值計算速度等因素,本文選取網(wǎng)格數(shù)量Ne=1.0×106的條件進行數(shù)值模擬研究。

    1.3 模型驗證

    本文的通氣率系數(shù)為

    Cq=Q/(v∞D(zhuǎn)2),

    (1)

    式中:Q為氣體體積流量。

    在相同的Cq和v∞條件下,將文獻[10]的水洞實驗結(jié)果和本文數(shù)值模擬結(jié)果對比,如圖4和圖5所示。圖4分別給出了數(shù)值模擬和實驗空泡形態(tài)發(fā)展過程,上方是數(shù)值模擬結(jié)果(本文以空氣體積分數(shù)為0.5作為空泡界面),下方是實驗結(jié)果,可見空泡尾部都出現(xiàn)了渦管泄氣的現(xiàn)象。

    圖4 數(shù)值模擬結(jié)果和實驗結(jié)果對比Fig.4 Comparison of simulated and experimental results

    圖5 數(shù)值模擬和實驗的空泡尺度對比Fig.5 Comparison of simulated and experimental cavity dimensions

    空泡尺度的對比結(jié)果如圖5所示。由圖5可以看出,空泡長度Lc和空泡直徑Dc的變化曲線具有很好的一致性,并且數(shù)值模擬結(jié)果和實驗結(jié)果的空泡尺度相差在6%以內(nèi)。圖4和圖5表明,數(shù)值模擬結(jié)果與實驗結(jié)果吻合較好,選用的數(shù)學(xué)模型可以很好地預(yù)測通氣空泡多相流動現(xiàn)象。

    2 模擬結(jié)果分析與討論

    2.1 空泡尾流結(jié)構(gòu)

    航行體在不同工況下的空泡形態(tài)如表1所示。

    表1 不同工況下航行體的空泡形態(tài)

    Tab.1 Supercavity shapes under different working conditions

    比較表1中工況1和工況2,發(fā)現(xiàn):在無沾濕區(qū)域條件下,通氣超空泡的尾流在重力環(huán)境中呈現(xiàn)雙渦管的泄氣形式;工況1僅起對比作用,表明在給定流速條件下,不計重力環(huán)境的壓差因素時,空泡尾部為單一渦管的流動結(jié)構(gòu)。

    航行體后體沾濕區(qū)域主要為航行體提供縱平面內(nèi)平衡的升力作用[13-14]。比較表1中工況2、工況3、工況4可以發(fā)現(xiàn),隨著攻角α的增加,航行體沾濕區(qū)域面積逐步增大,航行體超空泡尾流結(jié)構(gòu)從雙渦發(fā)展到三渦,表明沾濕區(qū)域?qū)张菸擦鹘Y(jié)構(gòu)存在影響。

    表1中工況5僅起與工況4的對比作用,給出了不計環(huán)境壓差因素時航行體的超空泡尾部形態(tài)。比較工況4、工況5可以發(fā)現(xiàn),在相同攻角(α=1.8°)條件下,兩種工況中航行體的超空泡形狀和沾濕區(qū)域都非常相似。二者環(huán)境壓差不同,進一步表明航行體沾濕區(qū)域是影響超空泡尾部流場結(jié)構(gòu)的重要因素。

    航行體直線運動時空泡形態(tài)演化過程如圖6所示。由圖6可見,0.50 s以后,通氣超空泡已經(jīng)包裹整個航行體,形態(tài)尺度趨于穩(wěn)定,空化數(shù)σ≈0.028,此時壓力曲線平直,整個航行體表面壓力分布趨于平穩(wěn),空泡內(nèi)部壓力差別較小[11]。

    圖6 航行體空泡演化Fig.6 Evolution of cavity of vehicle

    本文對壓力系數(shù)Cp及渦量Ω分別定義如下:

    (2)

    (3)

    式中:p為流域壓力;Ωx、Ωy、Ωz分別表示x軸、y軸、z軸方向的渦量,渦量的符號按照右手螺旋定則取值;vx、vy、vz分別表示x軸、y軸、z軸方向的流動速度。

    空泡在0.30 s形成自由閉合,航行體從0.30 s開始勻速改變攻角至0.34 s結(jié)束,航行體的攻角α從0°增加到1.8°. 航行體尾流場的壓力及渦量變化如圖7所示。

    圖7 空泡尾部壓力及渦量云圖Fig.7 Contours of pressure and vorticity at the tail of supercavity

    統(tǒng)計圖7中空泡尾部流場最右側(cè)截面在航行體攻角改變過程中空泡尾渦的非定常變化,結(jié)果如表2所示。表2中S0為沾濕區(qū)域面積,S為航行體表面積。

    表2 空泡尾流橫截面渦方向的非定常變化

    Tab.2 Unsteady variation of vortex direction in cross section of cavity wake

    從表2中可以看出:t取值為0.300~0.310 s時,空泡尾部的雙渦管逐漸形成,流域環(huán)境壓差使空泡尾部主要以雙渦管閉合為主,雖然有小部分沾濕區(qū)域影響,但S0/S只有0.87%左右,空泡尾部的渦量方向在此期間保持一致,為逆時針和順時針,如圖7(b)所示;t取值為0.310~0.319 s時,空泡尾部雙渦逐漸消失,隨著S0/S的不斷增大,空泡尾部形成四渦,該情況在圣安東尼瀑布實驗室的研究中也有出現(xiàn)[5],實驗指出四渦情況發(fā)生在超空泡尾流從回射流到雙渦管之間的轉(zhuǎn)變過程中,由于數(shù)值模擬研究中選擇空氣體積分數(shù)不同以及數(shù)值方法在處理細微氣泡的旋渦結(jié)構(gòu)時還與實驗結(jié)果存在一定差異,故數(shù)值模擬暫時還不能精細地反映實驗中渦管細節(jié);t取值為0.319~0.321 s時,隨著沾濕區(qū)域?qū)翰钣绊懙脑龃?,四渦中的上方雙渦逐漸減小,下方雙渦逐漸增大,渦量云圖如圖7(c)所示;t取值為0.321~0.322 s時,隨著S0/S的進一步擴大,空泡尾流上方雙渦逐漸被下方雙渦取代;t取值為0.322~0.340 s時,空泡尾部上方的渦管尺度逐漸增大,形成三渦管尾流,上方渦管流域的渦量很小,下方雙渦的渦量保持順時針和逆時針方向,如圖7(d)所示;t取值為0.340~0.440 s時,航行體攻角保持固定,空泡形態(tài)趨于穩(wěn)定,尾渦形態(tài)如圖7(e)所示。

    2.2 空泡尾部壓力

    航行體尾部渦管附近的流線方向如圖8所示,從圖8(a)和圖8(b)中的流線方向?qū)Ρ瓤梢园l(fā)現(xiàn),沾濕區(qū)域致使空泡尾部橫截面中渦管附近的流線方向發(fā)生了改變。

    圖8 空泡尾部渦管附近的壓力和流線Fig.8 Pressure and streamline near vortex tube at cavity tail

    圖9 攻角、沾濕區(qū)域和空泡尾流高壓區(qū)無量綱壓力差系數(shù)隨時間變化曲線Fig.9 Changes of angle-of-attack, wetted area dimensionless pressure difference in the high pressure region of wake flow of cavity with time

    由圖9可以看出:t取值為0.300~0.318 s時,受流域壓差影響,上方高壓區(qū)的壓力小于下方,雙渦管之間的流動方向向上,在空泡尾部形成一對中心上卷的雙渦;t取值為0.318~0.322 s時,航行體的沾濕區(qū)域開始逐漸影響尾流壓力分布,使上下方高壓區(qū)逐漸融合,壓力差減??;t取值為0.322~0.340 s時,攻角不斷增加,沾濕區(qū)域繞流的誘導(dǎo)作用對尾流壓力分布的影響進一步擴大,高壓區(qū)再次分成上下兩個部分,上方高壓區(qū)的壓力大于下方,壓力差值逐漸增加,雙渦管之間的介質(zhì)轉(zhuǎn)為向下流動,在空泡尾部形成一對中心下卷的雙渦;t取值為0.340 s后,航行體攻角不在變化,沾濕區(qū)域趨于穩(wěn)定。

    2.3 空泡尾部渦量

    為進一步分析四渦出現(xiàn)期間沾濕區(qū)域?qū)ξ矞u方向的影響規(guī)律,考慮到流域環(huán)境壓差和沾濕區(qū)域?qū)张菸矞u存在相反方向的影響,可以通過空泡尾部上下流線的速度分布進行深入研究[7-8],如圖10所示??张菘v向平面上圍繞整個空泡閉合輪廓線的環(huán)量Γ[15]近似表示為

    (4)

    式中:v為空泡輪廓線速度矢量;s為空泡輪廓線長度矢量;l為空泡輪廓線周長;vx1和vx2分別為空泡尾部上、下表面流線的x軸方向速度。

    圖10 空泡尾部上下流線的速度vx1和vx2Fig.10 Velocities vx1 and vx2 of upper and lower streamlines at cavity tail

    由于空泡上下表面的速度差異在空泡尾部處較大,用vx1曲線和vx2曲線所圍的面積可以近似地表征空泡尾部的渦量變化,如圖11所示。圖11(a)中空泡剛好閉合在航行體尾部,受流域環(huán)境壓差作用,空泡下表面的速度略大[11],形成雙渦尾流。圖11(b)中空泡尾流既受環(huán)境壓差影響,又受到不斷擴大的沾濕區(qū)域影響,空泡下表面速度vx2受沾濕區(qū)域阻擋逐漸小于vx1,渦量方向(順時針區(qū)域和逆時針區(qū)域的面積之差)變化較快,空泡尾部形成了四渦流動結(jié)構(gòu)。圖11(c)中,空泡尾流隨著沾濕區(qū)域的進一步擴大,空泡上表面的速度逐漸大于下表面,渦量方向發(fā)生變化,再次形成方向相反的雙渦尾流。

    圖11 空泡尾部上下流線的速度比較(g=9.8 m/s2)Fig.11 Comparison of velocities of upper and lower streamlines at cavity tail(g=9.8 m/s2)

    進一步研究沾濕區(qū)域和流域環(huán)境壓差對空泡尾渦影響的主次關(guān)系,對比分析了不考慮環(huán)境壓差條件下(如工況1,假設(shè)g=0)航行體攻角變化時,沾濕區(qū)域?qū)张菸膊繙u量的影響。如圖12所示,0.300 s時航行體攻角α=0°,空泡尾部不受沾濕區(qū)域影響,航行體空泡尾部上下速度差異較小,空泡橫截面上有4個尾渦。0.340 s時航行體攻角α=1.8°,沾濕面積比例為36%左右,沾濕區(qū)域影響了航行體尾部的流速分布,空泡尾部出現(xiàn)對稱分布的渦量。

    圖12 空泡尾部上下流線的速度比較(g=0 m/s2)Fig.12 Comparison of velocities of upper and lower streamlines at cavity tail(g=0 m/s2)

    將vx1和vx2沿x軸方向積分,兩條曲線所圍的無量綱面積A為

    (5)

    用A的統(tǒng)計結(jié)果分析比較考慮沾濕區(qū)域引起的流速差異和環(huán)境壓差引起的流速差異,比較二者對空泡尾部渦量方向的影響關(guān)系,如圖13所示。

    圖13 沾濕區(qū)域?qū)张菸矞u的影響規(guī)律Fig.13 Influence law of wetted area on cavity wake vortex

    從圖13中可以看出,沾濕區(qū)域影響的曲線和流域環(huán)境壓差影響的曲線相交在S0/S處于7~10%之間:當(dāng)S0/S<7%時,環(huán)境壓差對空泡尾部渦量方向起到主導(dǎo)作用,呈現(xiàn)一對中心上卷的雙渦;當(dāng)7%10%時,沾濕區(qū)域繞流的誘導(dǎo)作用對空泡尾部渦量方向起到主導(dǎo)作用,呈現(xiàn)一對中心下卷的雙渦。另外,由于流動過程中存在多相流之間復(fù)雜的擾動影響,兩條曲線交點處S0/S的數(shù)值是一個范圍。

    3 結(jié)論

    本文研究了超空泡航行體沾濕區(qū)域?qū)张菸擦鹘Y(jié)構(gòu)的影響,對空泡尾流結(jié)構(gòu)的壓力及渦量的非定常變化進行了深入分析,得到了沾濕區(qū)域?qū)张菸擦鹘Y(jié)構(gòu)的影響規(guī)律。主要得出以下結(jié)論:

    1)揭示了低弗勞德數(shù)條件下通氣超空泡尾流結(jié)構(gòu)特性。航行體無攻角時,空泡尾部形成一對中心上卷的雙渦;航行體帶有攻角的條件下,空泡尾部形成一對中心下卷的雙渦。

    2)給出了沾濕區(qū)域?qū)ξ矞u方向的影響規(guī)律。在考慮重力環(huán)境影響的條件下,沾濕區(qū)域比例小于7%時,環(huán)境壓差是空泡尾渦方向的主導(dǎo)因素;沾濕區(qū)域比例大于10%時,沾濕區(qū)域是空泡尾渦方向的主導(dǎo)因素。

    3)指出了空泡的四渦尾流結(jié)構(gòu)產(chǎn)生的原因。沾濕區(qū)域的比例處于7%~10%之間時,在環(huán)境壓差和沾濕區(qū)域的雙重作用下,空泡尾部流場中出現(xiàn)了四渦流動結(jié)構(gòu)。

    致謝本文撰寫過程中得到了哈爾濱工業(yè)大學(xué)航天學(xué)院曹登慶教授和徐敏強教授的大力幫助與支持,在此對兩位老師表示衷心的感謝!

    猜你喜歡
    渦量尾流空泡
    水下航行體雙空泡相互作用數(shù)值模擬研究
    含沙空化對軸流泵內(nèi)渦量分布的影響
    自由表面渦流動現(xiàn)象的數(shù)值模擬
    飛機尾流的散射特性與探測技術(shù)綜述
    基于LPV的超空泡航行體H∞抗飽和控制
    基于CFD的對轉(zhuǎn)槳無空泡噪聲的仿真預(yù)報
    船海工程(2015年4期)2016-01-05 15:53:28
    錐形流量計尾流流場分析
    水面艦船風(fēng)尾流效應(yīng)減弱的模擬研究
    航態(tài)對大型船舶甲板氣流場的影響
    The application of numerical simulation of delta wing with blunt leading edge using RANS/LES hybrid method
    国产精品久久久久久久久免| 亚洲无线观看免费| 久久久久国产网址| a 毛片基地| 又粗又硬又长又爽又黄的视频| 亚洲欧美精品自产自拍| 丁香六月天网| 男女国产视频网站| 丝袜喷水一区| www.色视频.com| 少妇熟女欧美另类| 观看美女的网站| 九九爱精品视频在线观看| 国产精品一区二区性色av| 日韩成人av中文字幕在线观看| 亚洲精品乱码久久久久久按摩| 99热这里只有是精品50| 美女内射精品一级片tv| av在线观看视频网站免费| 午夜视频国产福利| 国产日韩欧美视频二区| 日本av手机在线免费观看| 国产成人91sexporn| 五月开心婷婷网| 天堂俺去俺来也www色官网| 内地一区二区视频在线| 91精品伊人久久大香线蕉| 日本vs欧美在线观看视频 | 少妇人妻一区二区三区视频| 乱人伦中国视频| 成人无遮挡网站| 亚洲国产欧美在线一区| av播播在线观看一区| 两个人免费观看高清视频 | 中文天堂在线官网| 在线观看人妻少妇| 好男人视频免费观看在线| 桃花免费在线播放| 内射极品少妇av片p| 老熟女久久久| 日本av手机在线免费观看| 一本色道久久久久久精品综合| 免费观看在线日韩| 精品视频人人做人人爽| 久热久热在线精品观看| av国产精品久久久久影院| 欧美+日韩+精品| 国产欧美日韩一区二区三区在线 | 爱豆传媒免费全集在线观看| 你懂的网址亚洲精品在线观看| 99热这里只有精品一区| 久久午夜福利片| 2018国产大陆天天弄谢| 91aial.com中文字幕在线观看| 日本黄色日本黄色录像| 国产日韩欧美亚洲二区| 亚洲欧美日韩另类电影网站| 国产男女内射视频| 自拍偷自拍亚洲精品老妇| 亚洲国产av新网站| 亚洲情色 制服丝袜| 寂寞人妻少妇视频99o| 国产在线男女| 欧美成人精品欧美一级黄| 亚洲精品日本国产第一区| 九色成人免费人妻av| 日本午夜av视频| 永久网站在线| 精品视频人人做人人爽| 黄色视频在线播放观看不卡| av一本久久久久| 男女啪啪激烈高潮av片| 十八禁网站网址无遮挡 | 亚洲真实伦在线观看| 国产欧美日韩精品一区二区| 女的被弄到高潮叫床怎么办| 久久精品国产亚洲av天美| 少妇精品久久久久久久| 日日啪夜夜爽| 亚洲精品一区蜜桃| 精品亚洲成国产av| 国产免费又黄又爽又色| 欧美成人精品欧美一级黄| 亚洲国产成人一精品久久久| 久久久久精品性色| 99热6这里只有精品| 久久久久精品性色| 在线免费观看不下载黄p国产| 国产黄色视频一区二区在线观看| 免费在线观看成人毛片| 成人二区视频| 超碰97精品在线观看| 蜜桃久久精品国产亚洲av| 国产精品.久久久| 亚洲人成网站在线观看播放| 在线观看免费视频网站a站| 国产精品国产av在线观看| 亚洲精品一二三| 成人午夜精彩视频在线观看| 伊人久久精品亚洲午夜| 亚洲中文av在线| 亚洲av二区三区四区| 99视频精品全部免费 在线| 久久人妻熟女aⅴ| 人妻夜夜爽99麻豆av| 国产精品国产三级专区第一集| 国产一区二区在线观看日韩| 亚洲一级一片aⅴ在线观看| 国产日韩一区二区三区精品不卡 | 亚洲精品,欧美精品| 久久久久久久久大av| 国产精品国产三级国产av玫瑰| 亚洲欧洲精品一区二区精品久久久 | 你懂的网址亚洲精品在线观看| av.在线天堂| 亚洲激情五月婷婷啪啪| 国产午夜精品久久久久久一区二区三区| 自拍偷自拍亚洲精品老妇| 精品久久久久久电影网| 中文精品一卡2卡3卡4更新| 大话2 男鬼变身卡| 精品国产一区二区三区久久久樱花| 一个人看视频在线观看www免费| av视频免费观看在线观看| 久久精品夜色国产| 国产精品一区www在线观看| 日日摸夜夜添夜夜爱| 亚洲精品久久久久久婷婷小说| 99精国产麻豆久久婷婷| 少妇高潮的动态图| 亚洲欧美成人精品一区二区| 欧美日韩av久久| 大香蕉97超碰在线| 9色porny在线观看| 亚洲av福利一区| 婷婷色麻豆天堂久久| 亚洲欧美精品自产自拍| 国产精品一区二区三区四区免费观看| 国产深夜福利视频在线观看| 日本欧美国产在线视频| 青青草视频在线视频观看| a级毛片免费高清观看在线播放| 汤姆久久久久久久影院中文字幕| 成人亚洲欧美一区二区av| 高清视频免费观看一区二区| 蜜桃在线观看..| 国产精品蜜桃在线观看| 免费观看的影片在线观看| 男男h啪啪无遮挡| 在线观看美女被高潮喷水网站| 我要看日韩黄色一级片| 人妻一区二区av| 国产精品一区二区性色av| 欧美少妇被猛烈插入视频| 美女福利国产在线| 精品久久久久久电影网| av天堂中文字幕网| 91精品一卡2卡3卡4卡| 国产日韩一区二区三区精品不卡 | 国产精品人妻久久久影院| 日韩伦理黄色片| 亚洲av中文av极速乱| 人妻夜夜爽99麻豆av| 亚洲欧美一区二区三区黑人 | 美女福利国产在线| 久久ye,这里只有精品| 99久久中文字幕三级久久日本| 少妇被粗大猛烈的视频| 人人妻人人澡人人看| 精品人妻偷拍中文字幕| 两个人的视频大全免费| h日本视频在线播放| 看免费成人av毛片| 久热这里只有精品99| 高清毛片免费看| 中文字幕制服av| 性色av一级| 午夜激情福利司机影院| 曰老女人黄片| 最近手机中文字幕大全| 国产淫片久久久久久久久| 丰满少妇做爰视频| 能在线免费看毛片的网站| 久久青草综合色| 国产高清有码在线观看视频| 久久午夜综合久久蜜桃| 久久ye,这里只有精品| 大片电影免费在线观看免费| 久久久久久久久大av| 又大又黄又爽视频免费| 最新中文字幕久久久久| 亚洲熟女精品中文字幕| 2022亚洲国产成人精品| 69精品国产乱码久久久| 大码成人一级视频| 日韩不卡一区二区三区视频在线| 午夜av观看不卡| 亚洲av中文av极速乱| 国产一区二区三区av在线| 自拍偷自拍亚洲精品老妇| 日韩中字成人| 如日韩欧美国产精品一区二区三区 | 3wmmmm亚洲av在线观看| 国产av精品麻豆| 尾随美女入室| 午夜免费观看性视频| a级片在线免费高清观看视频| 十八禁网站网址无遮挡 | 美女xxoo啪啪120秒动态图| 国产黄频视频在线观看| 日韩制服骚丝袜av| 亚洲欧洲精品一区二区精品久久久 | 一二三四中文在线观看免费高清| 国产爽快片一区二区三区| 亚洲欧洲精品一区二区精品久久久 | 亚洲天堂av无毛| 少妇被粗大猛烈的视频| 丰满乱子伦码专区| 久久 成人 亚洲| 色视频www国产| 成人美女网站在线观看视频| 色94色欧美一区二区| 日本午夜av视频| 国产日韩一区二区三区精品不卡 | 日韩大片免费观看网站| 日本猛色少妇xxxxx猛交久久| 毛片一级片免费看久久久久| 啦啦啦啦在线视频资源| 国产欧美日韩精品一区二区| 欧美日韩综合久久久久久| 能在线免费看毛片的网站| 26uuu在线亚洲综合色| 高清av免费在线| 久久精品久久精品一区二区三区| 最后的刺客免费高清国语| 香蕉精品网在线| 精品视频人人做人人爽| 黑丝袜美女国产一区| 在线观看免费日韩欧美大片 | 久久久国产一区二区| av卡一久久| 亚洲欧美日韩另类电影网站| 丝瓜视频免费看黄片| 欧美日本中文国产一区发布| 亚洲精品aⅴ在线观看| 日韩不卡一区二区三区视频在线| 91久久精品国产一区二区三区| 黄色毛片三级朝国网站 | 久久久午夜欧美精品| 一级片'在线观看视频| 欧美另类一区| 久久99热这里只频精品6学生| 插阴视频在线观看视频| 亚洲精品日本国产第一区| 亚洲精品中文字幕在线视频 | 一级毛片黄色毛片免费观看视频| 国产精品成人在线| 五月玫瑰六月丁香| 91久久精品国产一区二区成人| 欧美日韩av久久| 新久久久久国产一级毛片| 久久鲁丝午夜福利片| 日韩人妻高清精品专区| 国产伦精品一区二区三区视频9| 一本—道久久a久久精品蜜桃钙片| 免费大片18禁| 国产 精品1| av线在线观看网站| a级片在线免费高清观看视频| 国产有黄有色有爽视频| 免费播放大片免费观看视频在线观看| 国产高清国产精品国产三级| 曰老女人黄片| 国产日韩欧美亚洲二区| 久久久久久久久久久免费av| 国产精品福利在线免费观看| 国产精品人妻久久久影院| 精品国产国语对白av| 欧美丝袜亚洲另类| 蜜桃在线观看..| 国产精品一区二区性色av| 九九久久精品国产亚洲av麻豆| 黄色怎么调成土黄色| 欧美日韩av久久| 久久久国产欧美日韩av| 国产在线免费精品| 最新的欧美精品一区二区| 十分钟在线观看高清视频www | 欧美丝袜亚洲另类| 两个人免费观看高清视频 | 综合色丁香网| 国产日韩欧美视频二区| 国产欧美日韩精品一区二区| 夫妻性生交免费视频一级片| 精品国产乱码久久久久久小说| 人人澡人人妻人| 国产美女午夜福利| 亚洲av成人精品一区久久| h日本视频在线播放| 欧美bdsm另类| 国产精品女同一区二区软件| 国产精品熟女久久久久浪| 老熟女久久久| 看十八女毛片水多多多| 国产精品成人在线| 久久ye,这里只有精品| 狂野欧美白嫩少妇大欣赏| 欧美97在线视频| 肉色欧美久久久久久久蜜桃| 日日撸夜夜添| 精品少妇内射三级| 女的被弄到高潮叫床怎么办| 色5月婷婷丁香| 91aial.com中文字幕在线观看| 18禁在线播放成人免费| 国产无遮挡羞羞视频在线观看| 国产伦在线观看视频一区| 在线观看av片永久免费下载| 99久久中文字幕三级久久日本| 国产精品久久久久久av不卡| 国产精品久久久久久精品古装| 久久毛片免费看一区二区三区| 男的添女的下面高潮视频| 久久久久视频综合| 美女cb高潮喷水在线观看| 亚洲av在线观看美女高潮| 国产黄片美女视频| 少妇人妻一区二区三区视频| 9色porny在线观看| 亚洲欧美中文字幕日韩二区| 国产精品麻豆人妻色哟哟久久| 久久久久久久久久人人人人人人| 久久久久精品久久久久真实原创| 美女xxoo啪啪120秒动态图| 国产黄色视频一区二区在线观看| 精品一区二区三卡| 午夜福利在线观看免费完整高清在| 亚洲av.av天堂| 久久午夜福利片| 一级爰片在线观看| 久久人妻熟女aⅴ| 亚洲欧美成人综合另类久久久| 国产日韩一区二区三区精品不卡 | 亚洲中文av在线| 寂寞人妻少妇视频99o| 男人和女人高潮做爰伦理| 亚洲不卡免费看| 人人妻人人澡人人爽人人夜夜| 亚洲精品,欧美精品| 国产精品女同一区二区软件| 另类精品久久| 久久久久久久久久成人| 欧美日韩视频精品一区| 午夜福利视频精品| 国产深夜福利视频在线观看| 国产有黄有色有爽视频| 一边亲一边摸免费视频| 麻豆乱淫一区二区| 在线播放无遮挡| 欧美精品一区二区大全| 草草在线视频免费看| 2021少妇久久久久久久久久久| 亚洲精品一二三| 免费播放大片免费观看视频在线观看| 婷婷色综合www| 男女无遮挡免费网站观看| 国产一区二区三区av在线| 五月天丁香电影| 久久久欧美国产精品| 黑丝袜美女国产一区| 五月玫瑰六月丁香| 在线 av 中文字幕| 国产毛片在线视频| 五月伊人婷婷丁香| 能在线免费看毛片的网站| 观看免费一级毛片| 国产永久视频网站| 久久人人爽av亚洲精品天堂| 三上悠亚av全集在线观看 | 视频中文字幕在线观看| 久久久久久久大尺度免费视频| 国产精品.久久久| 国产一区二区在线观看av| 亚洲国产色片| 国产精品不卡视频一区二区| 亚洲精品日韩av片在线观看| 免费观看的影片在线观看| 人人妻人人看人人澡| 视频区图区小说| 国产 精品1| 国产精品一区二区三区四区免费观看| 日韩中文字幕视频在线看片| 日韩一区二区三区影片| av天堂中文字幕网| 黑人高潮一二区| 日韩一区二区三区影片| 高清视频免费观看一区二区| www.色视频.com| 最近2019中文字幕mv第一页| 少妇高潮的动态图| 亚洲国产欧美在线一区| 91精品国产九色| 天堂中文最新版在线下载| 久久人人爽av亚洲精品天堂| 精品亚洲乱码少妇综合久久| 久久久久久久久久久久大奶| 视频中文字幕在线观看| 精品亚洲成国产av| 日韩一本色道免费dvd| 91精品伊人久久大香线蕉| 日韩精品免费视频一区二区三区 | 精品久久国产蜜桃| 色视频www国产| 久久久久久久久久久免费av| 亚洲欧洲日产国产| 国产欧美日韩一区二区三区在线 | 国产av码专区亚洲av| 99九九线精品视频在线观看视频| av黄色大香蕉| 最新的欧美精品一区二区| 91久久精品国产一区二区成人| 国产欧美日韩精品一区二区| 精品国产乱码久久久久久小说| 亚洲精品国产av成人精品| 日韩成人伦理影院| 亚洲伊人久久精品综合| 国产日韩欧美视频二区| 免费观看在线日韩| 日韩成人av中文字幕在线观看| 夜夜爽夜夜爽视频| 亚洲欧洲国产日韩| 少妇人妻精品综合一区二区| 亚洲av二区三区四区| 亚洲欧洲日产国产| 蜜桃久久精品国产亚洲av| 边亲边吃奶的免费视频| 免费大片18禁| 日韩一区二区视频免费看| 最黄视频免费看| a 毛片基地| 久久国产精品大桥未久av | 久久狼人影院| 久久精品国产自在天天线| 精品人妻熟女av久视频| .国产精品久久| 久久久久久久精品精品| 久久久欧美国产精品| av免费在线看不卡| 啦啦啦啦在线视频资源| 成人毛片a级毛片在线播放| 国产免费视频播放在线视频| 亚洲精华国产精华液的使用体验| 18禁裸乳无遮挡动漫免费视频| 午夜久久久在线观看| 国产69精品久久久久777片| 两个人的视频大全免费| 丰满少妇做爰视频| 久久精品国产亚洲av涩爱| 亚洲国产精品成人久久小说| 久久久久久久久久久免费av| av一本久久久久| 国产一区有黄有色的免费视频| 91aial.com中文字幕在线观看| 熟女人妻精品中文字幕| 色94色欧美一区二区| 成人漫画全彩无遮挡| 国产伦精品一区二区三区四那| 18禁裸乳无遮挡动漫免费视频| 香蕉精品网在线| 天天躁夜夜躁狠狠久久av| 最近手机中文字幕大全| 成年人免费黄色播放视频 | 久久国内精品自在自线图片| 成年人免费黄色播放视频 | 妹子高潮喷水视频| 国产成人精品无人区| 最后的刺客免费高清国语| 精品人妻熟女毛片av久久网站| 国产亚洲最大av| av在线老鸭窝| 亚洲综合色惰| 国产精品久久久久久久久免| 日韩av不卡免费在线播放| 亚洲四区av| 老司机影院毛片| 亚洲欧美日韩卡通动漫| 国产成人午夜福利电影在线观看| 日本wwww免费看| 成人二区视频| 亚洲一区二区三区欧美精品| 亚洲四区av| 中文字幕人妻丝袜制服| 美女福利国产在线| 少妇被粗大的猛进出69影院 | 韩国av在线不卡| 边亲边吃奶的免费视频| 我要看黄色一级片免费的| 免费av不卡在线播放| 久久99热这里只频精品6学生| kizo精华| 大码成人一级视频| 十八禁高潮呻吟视频 | 另类精品久久| 国产极品天堂在线| 狂野欧美激情性bbbbbb| 亚洲欧美一区二区三区黑人 | 99久久精品国产国产毛片| 国产精品一区二区三区四区免费观看| 国产日韩一区二区三区精品不卡 | 免费在线观看成人毛片| 中文字幕人妻丝袜制服| 亚洲色图综合在线观看| 赤兔流量卡办理| 国产亚洲欧美精品永久| 精品久久久噜噜| 最近最新中文字幕免费大全7| 久久久久久久久久久丰满| 国产一区二区在线观看日韩| 亚洲av综合色区一区| 蜜臀久久99精品久久宅男| 狂野欧美白嫩少妇大欣赏| 中文天堂在线官网| 久久久久精品性色| 天天操日日干夜夜撸| 国内精品宾馆在线| 插逼视频在线观看| 99久久综合免费| 美女脱内裤让男人舔精品视频| 精品少妇内射三级| 久久鲁丝午夜福利片| 精品少妇久久久久久888优播| 人妻一区二区av| 精品一区在线观看国产| 国产又色又爽无遮挡免| 大片电影免费在线观看免费| 精品亚洲成a人片在线观看| 婷婷色综合大香蕉| 免费看日本二区| 69精品国产乱码久久久| 亚洲国产成人一精品久久久| 最黄视频免费看| 亚洲国产欧美日韩在线播放 | 最近的中文字幕免费完整| 亚洲精品自拍成人| av女优亚洲男人天堂| 男的添女的下面高潮视频| 欧美xxxx性猛交bbbb| 在线观看免费日韩欧美大片 | 亚洲人与动物交配视频| 成人毛片a级毛片在线播放| 欧美精品亚洲一区二区| 国产永久视频网站| 国产美女午夜福利| 街头女战士在线观看网站| 婷婷色av中文字幕| 久久99精品国语久久久| 亚洲成人一二三区av| 天天躁夜夜躁狠狠久久av| 午夜福利视频精品| 人妻一区二区av| 亚州av有码| 日本黄色片子视频| 亚洲成色77777| 久久久国产一区二区| 免费在线观看成人毛片| 成年女人在线观看亚洲视频| 人体艺术视频欧美日本| 这个男人来自地球电影免费观看 | 久久人人爽人人片av| 男女免费视频国产| 午夜91福利影院| 亚洲欧美日韩东京热| 亚洲,一卡二卡三卡| 精品酒店卫生间| 最后的刺客免费高清国语| 国产成人精品无人区| 狂野欧美白嫩少妇大欣赏| 国产黄频视频在线观看| 18禁动态无遮挡网站| 深夜a级毛片| 国产精品久久久久久久久免| 日韩精品有码人妻一区| 丁香六月天网| 爱豆传媒免费全集在线观看| 亚洲av成人精品一区久久| 亚洲内射少妇av| 内射极品少妇av片p| av卡一久久| 免费观看av网站的网址| av女优亚洲男人天堂| 日日啪夜夜撸| 国产69精品久久久久777片| 欧美bdsm另类| 国产亚洲av片在线观看秒播厂| 搡女人真爽免费视频火全软件| 亚洲av成人精品一区久久| 国产精品三级大全| 日韩精品免费视频一区二区三区 | 亚洲人成网站在线播| 亚洲欧美成人综合另类久久久| av线在线观看网站| 高清av免费在线| 亚洲,欧美,日韩| 免费大片18禁| 久久国产精品大桥未久av | 91午夜精品亚洲一区二区三区| 在线亚洲精品国产二区图片欧美 | 青青草视频在线视频观看| 欧美区成人在线视频| 一级黄片播放器| 亚洲欧美成人综合另类久久久| 欧美97在线视频| 9色porny在线观看| 人人妻人人澡人人爽人人夜夜| 色婷婷久久久亚洲欧美|