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

    基于高精度Boussinesq方程的三維淺水晃蕩數(shù)值研究

    2021-06-01 03:41:12袁心怡劉祖源
    關(guān)鍵詞:液艙波峰邊界條件

    袁心怡, 蘇 焱, 劉祖源

    (武漢理工大學(xué) 交通學(xué)院, 武漢 430063)

    晃蕩可以描述為部分填充的容器中液體的運(yùn)動(dòng).液化天然氣(LNG)船、浮式液化天然氣裝置(FLNG)等采用液艙形式的液貨船在部分裝載的情況下,即使受到小幅激勵(lì),也有可能發(fā)生劇烈的晃蕩現(xiàn)象,對(duì)液艙結(jié)構(gòu)產(chǎn)生不利影響,威脅船舶航行安全.

    研究晃蕩問題多采用數(shù)值方法,基于黏性流理論在流域內(nèi)求解Navier-Stokes方程或基于勢流理論在流域內(nèi)求解Laplace方程實(shí)現(xiàn)數(shù)值模擬.Wang等[1]基于比例邊界有限元法開發(fā)了一種研究擋板對(duì)艙內(nèi)晃蕩影響的新方法,經(jīng)過收斂性和研究,發(fā)現(xiàn)該方法具有計(jì)算簡單、精度高、速度快等優(yōu)點(diǎn).Zhao等[2]基于勢流理論建立了LNG液艙在強(qiáng)制激勵(lì)下的非線性晃蕩數(shù)值模型,并加入了人工阻尼模型考慮耗能,研究了部分充液的矩形容器在平移和旋轉(zhuǎn)激勵(lì)下的晃蕩特性.Hu等[3]利用邊界元法求解控制方程,采用線性自由表面條件研究了多種形狀帶擋板的二維液艙在小幅激勵(lì)下晃蕩的固有頻率和振型,并提出了帶擋板的矩形液艙內(nèi)晃蕩的固有頻率計(jì)算經(jīng)驗(yàn)公式.ünal等[4]研究了帶T形擋板的液艙內(nèi)的晃蕩現(xiàn)象,考察了擋板高度、充液深度等參數(shù)對(duì)于艙壁水動(dòng)力載荷和自由表面高程的影響,計(jì)算采用了兩種不同的方法,分別使用黏性的層流和紊流求解器.李金龍等[5]對(duì)一種新的幾何有限體積法isoAdvector進(jìn)行修正,使之可應(yīng)用于液艙晃蕩的模擬中,并基于不同的有限體積法對(duì)各種工況進(jìn)行了模擬,發(fā)現(xiàn)修正后的方法所獲得的自由面位置和整體水動(dòng)力載荷載荷精度更高,且能較好地模擬波面的翻卷和破碎.Fu等[6]提出了一種基于局部徑向基函數(shù)和二階顯式龍格庫塔法的半拉格朗日無網(wǎng)格模型來分析二維艙內(nèi)的晃蕩現(xiàn)象,通過文獻(xiàn)對(duì)比證明了該方法的有效性和準(zhǔn)確性.

    對(duì)于強(qiáng)非線性的晃蕩問題,粒子法可以很好地捕捉自由面的大變形.但對(duì)于低充液率的容器,穩(wěn)定和精確地模擬容器內(nèi)長時(shí)間的晃蕩具有難度.Green等[7]提出了一種有效的光滑粒子流體動(dòng)力學(xué)(SPH)方法來克服這些困難,使邊界條件、耗散項(xiàng)和算法的誤差最小化,與試驗(yàn)數(shù)據(jù)比較發(fā)現(xiàn),數(shù)值模擬的結(jié)果非常精確.

    布西內(nèi)斯克(Boussinesq)方程是一類描述波浪傳播變形的數(shù)學(xué)模型.Bingham等[8]建立了數(shù)值離散方法,主要用于近岸變化底部處波浪場精確模擬.Antuono等[9]使用Boussinesq型平均水深方程建立了一個(gè)二維模型用于分析淺水中液體晃蕩現(xiàn)象.Su等[10]采用高精度Boussinesq型速度勢方程,對(duì)二維矩形艙內(nèi)的晃蕩現(xiàn)象進(jìn)行了數(shù)值模擬,結(jié)果表明基于Boussinesq型方程的數(shù)值模型在預(yù)測小激勵(lì)幅值晃蕩方面具有較好的性能.

    本研究旨在探討淺水、微幅激勵(lì)條件下三維液艙內(nèi)的晃蕩現(xiàn)象.假定液艙內(nèi)為無黏無旋不可壓流體,流域內(nèi)存在速度勢.將總速度勢分解為兩部分,一部分為滿足拉普拉斯方程和壁面邊界條件的特解,另一部分采用基于高精度速度勢型Boussinesq方程的數(shù)值模型求解.首先給出三維液艙中流體受強(qiáng)迫運(yùn)動(dòng)的控制方程和邊界條件,然后基于有限差分法對(duì)自由面進(jìn)行空間離散,建立數(shù)值模型.通過與文獻(xiàn)中二維晃蕩的結(jié)果進(jìn)行對(duì)比驗(yàn)證了數(shù)值模型的有效性,并對(duì)三維工況的結(jié)果進(jìn)行分析.

    1 數(shù)學(xué)模型

    考慮一個(gè)長為l、寬為b、靜止水深為h的液艙,如圖1所示,固定坐標(biāo)系O-xyz位于靜止水面上,原點(diǎn)位于靜水面的幾何中心,z軸垂直向上.液艙底部記為z=-h(x,y),流體自由表面記為z=η(x,y).

    圖1 晃蕩運(yùn)動(dòng)模型Fig.1 Sloshing motion model

    假定艙內(nèi)為無黏無旋不可壓流體,基于勢流理論,流域內(nèi)存在一個(gè)速度勢Φ.記水平速度為u,垂向速度為w,則速度勢與速度的關(guān)系為

    (1)

    ηt+Φxηx+Φyηy-Φz=0,z=η

    (2)

    (3)

    (4)

    ?Φ/?n=v·n

    (5)

    其中,式(2)為自由面運(yùn)動(dòng)學(xué)邊界條件;式(3)為自由面動(dòng)力學(xué)邊界條件;式(4)為拉普拉斯方程,Xij(i,j=x,y,z)表示物理量X對(duì)i方向求一階導(dǎo)數(shù),對(duì)j方向求二階導(dǎo)數(shù);式(5)為物面不可穿透條件;g為重力加速度;v=[vxvyvz]為液艙的運(yùn)動(dòng)速度;n為物面的法向量,由液艙內(nèi)指向外.對(duì)于慣性系中定義的自由面邊界條件,利用如下關(guān)系將其轉(zhuǎn)換到隨液艙運(yùn)動(dòng)坐標(biāo)系中:

    則可以給出隨艙運(yùn)動(dòng)坐標(biāo)系下,自由面滿足運(yùn)動(dòng)學(xué)和動(dòng)力學(xué)邊界條件:

    (8)

    (9)

    由于液艙的運(yùn)動(dòng)始終是已知量,在數(shù)值計(jì)算時(shí),將總速度勢分為兩部分:Φ=φ+φ,其中φ是一個(gè)滿足拉普拉斯方程和物面不可穿透條件的特解.僅考慮橫蕩、縱蕩、垂蕩運(yùn)動(dòng)時(shí),可以直接寫出φ的表達(dá)式如下:

    φ=xvx+yvy+zvz

    (10)

    然后由疊加原理可導(dǎo)出待求解的分速度勢φ所滿足的物面邊界條件和控制方程:

    (11)

    將Φ=φ+φ以及式(10)引入到式(8)、(9)中,可以得到φ滿足的自由面邊界條件如下:

    ηt+ηxφx+ηyφy-φz=0

    (12)

    (13)

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

    (14)

    則可用新變量表示自由面條件并整理如下:

    (15)

    (16)

    (17)

    式中:參數(shù)δ?1,表示展開面是緩慢變化的.由上述遞推關(guān)系式可推導(dǎo)速度勢和垂向速度的表達(dá)式如下:

    (18)

    (19)

    將表達(dá)垂向速度的無窮級(jí)數(shù)截?cái)嘣?階時(shí),以上兩式中:

    由此可得到速度勢和垂向速度的表達(dá)式,將其代入到底部運(yùn)動(dòng)學(xué)邊界條件a1、a3的表達(dá)式中得:

    (25)

    整理后可得底部方程的表達(dá)式為

    (26)

    數(shù)值計(jì)算過程中,空間導(dǎo)數(shù)計(jì)算采用有限差分法,時(shí)間步進(jìn)采用四階龍格庫塔方法.計(jì)算程序基于MATLAB語言編寫,求解線性方程組使用MATLAB中的直接求解器,該求解器可根據(jù)系數(shù)矩陣的特征選擇高效的求解方法.

    3 二維工況驗(yàn)證

    為驗(yàn)證數(shù)值模型的有效性,將三維液艙的寬度設(shè)置為小量并在相同工況下與文獻(xiàn)中的二維結(jié)果[10]進(jìn)行對(duì)比.三維液艙的寬度b=0.1 m,長度和充液高度都與二維工況相同,分別為l=1.175 m,h=0.06 m.液艙受長度方向的正弦激勵(lì)做縱蕩運(yùn)動(dòng),激勵(lì)幅值A(chǔ)=0.39 cm,記線性理論預(yù)測的共振頻率ωr=2.042 5 s-1,激勵(lì)頻率ω′分別為0.98ωr、1.02ωr及1.08ωr.

    圖2依次展示了激勵(lì)頻率為0.98ωr、1.02ωr及1.08ωr時(shí)左壁面處的液面高度變化對(duì)比圖,圖中:t′為時(shí)間與激勵(lì)周期之比,η′為自由面高度與充液高度之比.3個(gè)工況下的波高變化都呈現(xiàn)出強(qiáng)非線性現(xiàn)象,波峰陡峭但波谷平坦.當(dāng)外部激勵(lì)頻率處于共振頻率附近時(shí),矩形液艙內(nèi)的晃蕩運(yùn)動(dòng)形式對(duì)外部激勵(lì)頻率的變化非常敏感,在3種工況下觀察到了完全不同的運(yùn)動(dòng)形式,穩(wěn)定狀態(tài)下的一個(gè)變化周期內(nèi)分別包含3個(gè)波峰、2個(gè)波峰和單個(gè)波峰.相應(yīng)地,在液艙內(nèi)可以觀察到對(duì)應(yīng)數(shù)量的行進(jìn)波在艙內(nèi)來回移動(dòng).

    圖2 液面高度時(shí)歷對(duì)比圖Fig.2 Comparisons of time history of surface height

    可以看出,三維數(shù)值計(jì)算程序能很好地捕捉自由面高度的變化和晃蕩運(yùn)動(dòng)形式對(duì)于外部激勵(lì)頻率變化的敏感性.

    4 三維工況

    計(jì)算選擇長寬相等的方型液艙,將二維液艙的長度縮小一半,三維液艙模型參數(shù)為l=b=0.587 5 m,h=0.03 m.僅考慮縱蕩、橫蕩或二者混合的正弦形式運(yùn)動(dòng),A=0.1 cm,線性理論預(yù)測的艙內(nèi)晃蕩運(yùn)動(dòng)共振頻率為ωr=2.888 1 s-1.考慮到二維情況下,外部激勵(lì)頻率在共振頻率附近變化時(shí),晃蕩運(yùn)動(dòng)會(huì)出現(xiàn)特定的波高變化形式,三維工況計(jì)算也在共振頻率附近進(jìn)行.

    數(shù)值計(jì)算過程中,將靜止液面劃分成了正方形網(wǎng)格,x和y方向都有25個(gè)節(jié)點(diǎn).在二維液艙中,波高觀測點(diǎn)為左壁面處第1號(hào)節(jié)點(diǎn),而在三維液艙中,波高觀測點(diǎn)為混合運(yùn)動(dòng)合成矢量方向上的第1號(hào)節(jié)點(diǎn).

    如圖3所示,首先選取了ω′=0.98ωr,將二維工況和三維單方向運(yùn)動(dòng)工況進(jìn)行對(duì)比.可以看出,兩種情況下自由面的高度的變化幅值和周期都是一致的,液艙的寬度并沒有對(duì)艙內(nèi)流體的運(yùn)動(dòng)產(chǎn)生影響.接著將三維單方向運(yùn)動(dòng)工況和混合運(yùn)動(dòng)工況進(jìn)行了對(duì)比,從圖中可以看出,兩種工況下波高變化的周期是一致的,將單方向運(yùn)動(dòng)的波高加倍以后,幅值的大小也有很好的一致性.

    如圖4所示,計(jì)算了其他激勵(lì)頻率下,混合運(yùn)動(dòng)工況中波高觀測點(diǎn)處的波高變化.為了保證艙內(nèi)的流體運(yùn)動(dòng)達(dá)到穩(wěn)定狀態(tài),計(jì)算時(shí)間超過200個(gè)激勵(lì)周期.在1.02ωr激勵(lì)頻率下,觀察到了和二維工況中類似的現(xiàn)象:在一個(gè)變化周期內(nèi)出現(xiàn)兩個(gè)波峰.但是單波峰的現(xiàn)象卻出現(xiàn)在1.06ωr激勵(lì)頻率下,并且在1.08ωr激勵(lì)頻率下,觀測點(diǎn)波高變化的非線性減弱,說明三維情況下出現(xiàn)強(qiáng)非線性行進(jìn)波的激勵(lì)頻率域有所縮短.

    圖3 0.98ωr下波高對(duì)比圖Fig.3 Comparisons of wave height at 0.98ωr

    圖4 各頻率下波高變化圖Fig.4 Change of wave height at different frequencies

    數(shù)值計(jì)算程序可以實(shí)時(shí)顯示自由面變化,以1.06ωr激勵(lì)頻率為例,圖5顯示了數(shù)值模擬過程中某一時(shí)刻整個(gè)自由面形狀.在各個(gè)運(yùn)動(dòng)方向上,各自形成了行進(jìn)波,兩條波浪相遇的區(qū)域形成了一個(gè)波峰點(diǎn).類似地,在1.02ωr激勵(lì)頻率下,會(huì)在各個(gè)方向上觀察到兩條行進(jìn)波,且自由面上觀察到2×2個(gè)波峰點(diǎn).而在0.98ωr激勵(lì)頻率下,會(huì)在自由面上觀察到9條行進(jìn)波和3×3個(gè)波峰點(diǎn).

    圖5 自由面形狀圖Fig.5 Shape of free surface

    5 結(jié)論

    本研究基于高精度Boussinesq型速度勢方程對(duì)二維和三維液艙中在小幅激勵(lì)下的淺水晃蕩現(xiàn)象進(jìn)行了數(shù)值模擬.主要結(jié)論有:

    (1) 首次采用Boussinesq方程實(shí)現(xiàn)了三維淺水晃蕩模擬,建立的數(shù)值模型能很好地捕捉三維液艙內(nèi)的強(qiáng)非線性晃蕩運(yùn)動(dòng)以及晃蕩運(yùn)動(dòng)形式對(duì)于外部激勵(lì)頻率變化的敏感性,且具有較高的計(jì)算效率.

    (2) 三維數(shù)值模擬中,觀察到了4種不同的晃蕩運(yùn)動(dòng)形式.在各個(gè)工況下,會(huì)觀察到不同數(shù)量的行進(jìn)波,行進(jìn)波數(shù)量取決于外部激勵(lì)的頻率.在行進(jìn)波十字交匯的區(qū)域,會(huì)疊加產(chǎn)生波峰.

    (3) 在靠近共振頻率的工況下,即使是微幅激勵(lì)的耦合,都在艙內(nèi)產(chǎn)生了大幅的波浪.對(duì)于遠(yuǎn)離共振頻率的工況,相同外部激勵(lì)頻率下,與二維相比,三維情況下的非線性有所減弱.

    (4) 目前數(shù)值模型僅考慮單方向的平動(dòng)或兩方向耦合運(yùn)動(dòng),未來將對(duì)其他運(yùn)動(dòng)形式(單方向轉(zhuǎn)動(dòng)或耦合)作進(jìn)一步的研究.

    猜你喜歡
    液艙波峰邊界條件
    B型LNG液艙支座縱骨趾端處表面裂紋擴(kuò)展計(jì)算
    基于CFD的大型船舶液艙晃蕩研究
    一類帶有Stieltjes積分邊界條件的分?jǐn)?shù)階微分方程邊值問題正解
    帶有積分邊界條件的奇異攝動(dòng)邊值問題的漸近解
    作用于直立堤墻與樁柱的波峰高度分析計(jì)算
    考慮晃蕩效應(yīng)的獨(dú)立B型LNG液艙結(jié)構(gòu)多目標(biāo)優(yōu)化
    海洋工程(2016年2期)2016-10-12 05:08:07
    兒童標(biāo)準(zhǔn)12導(dǎo)聯(lián)T波峰末間期的分析
    FPSO與運(yùn)輸船旁靠時(shí)液艙晃蕩與船舶運(yùn)動(dòng)耦合效應(yīng)分析
    Dynamic Loads and Wake Prediction for Large Wind Turbines Based on Free Wake Method
    帶Robin邊界條件的2維隨機(jī)Ginzburg-Landau方程的吸引子
    午夜精品国产一区二区电影 | 久久久色成人| 久久韩国三级中文字幕| 亚洲av熟女| 最近中文字幕高清免费大全6| 亚洲av一区综合| 日韩欧美免费精品| 精品久久久久久久末码| 狂野欧美白嫩少妇大欣赏| 国产色婷婷99| 不卡一级毛片| 亚洲综合色惰| 春色校园在线视频观看| 欧美一级a爱片免费观看看| av天堂中文字幕网| 久久精品国产亚洲av香蕉五月| 天堂√8在线中文| 春色校园在线视频观看| 免费观看精品视频网站| 午夜福利高清视频| 露出奶头的视频| 成人av在线播放网站| 中文字幕人妻熟人妻熟丝袜美| 又黄又爽又免费观看的视频| avwww免费| 观看美女的网站| 少妇熟女欧美另类| 亚洲无线观看免费| 精华霜和精华液先用哪个| 国产精品精品国产色婷婷| 最近2019中文字幕mv第一页| 91在线观看av| 国产三级中文精品| 免费大片18禁| 精品午夜福利视频在线观看一区| 欧美人与善性xxx| 国产三级在线视频| 天美传媒精品一区二区| 一卡2卡三卡四卡精品乱码亚洲| 一个人看视频在线观看www免费| 12—13女人毛片做爰片一| 日本精品一区二区三区蜜桃| 99久久精品国产国产毛片| 日韩欧美国产在线观看| 寂寞人妻少妇视频99o| 国产精品国产三级国产av玫瑰| 一区二区三区免费毛片| 亚洲欧美精品综合久久99| 少妇的逼水好多| 日韩一区二区视频免费看| 床上黄色一级片| 日本免费一区二区三区高清不卡| 国产亚洲精品av在线| 网址你懂的国产日韩在线| 亚洲成人久久性| 久久久午夜欧美精品| 亚洲电影在线观看av| av在线蜜桃| 村上凉子中文字幕在线| 亚洲综合色惰| 久久综合国产亚洲精品| 国产精品久久久久久久电影| 毛片女人毛片| 欧美成人精品欧美一级黄| 日本黄色片子视频| aaaaa片日本免费| 美女高潮的动态| 亚洲人与动物交配视频| 精品一区二区免费观看| 女生性感内裤真人,穿戴方法视频| 亚洲最大成人av| 国产一区二区在线av高清观看| 国产黄a三级三级三级人| 日韩成人av中文字幕在线观看 | av在线老鸭窝| 真人做人爱边吃奶动态| 国产精品不卡视频一区二区| 国产精品国产高清国产av| 亚洲专区国产一区二区| 亚洲精品在线观看二区| 中国国产av一级| aaaaa片日本免费| 伦精品一区二区三区| 在线a可以看的网站| 午夜福利18| 欧美一区二区亚洲| 欧美国产日韩亚洲一区| 国产精品一区二区性色av| АⅤ资源中文在线天堂| 97超级碰碰碰精品色视频在线观看| 91狼人影院| 日本色播在线视频| 一本久久中文字幕| 人人妻人人看人人澡| 久久欧美精品欧美久久欧美| 欧美性感艳星| 国产在视频线在精品| 日韩,欧美,国产一区二区三区 | 日日摸夜夜添夜夜添小说| 精品一区二区三区av网在线观看| 国产精品99久久久久久久久| 久久精品人妻少妇| 99热网站在线观看| 丰满的人妻完整版| 国产精品国产高清国产av| 97碰自拍视频| 日韩欧美三级三区| 菩萨蛮人人尽说江南好唐韦庄 | 中出人妻视频一区二区| 免费人成在线观看视频色| 女同久久另类99精品国产91| 精品人妻视频免费看| 淫秽高清视频在线观看| 久99久视频精品免费| 欧美成人精品欧美一级黄| 免费观看精品视频网站| 久久久久国产精品人妻aⅴ院| 1000部很黄的大片| 自拍偷自拍亚洲精品老妇| 日本精品一区二区三区蜜桃| 久久韩国三级中文字幕| 亚洲国产精品成人久久小说 | 亚洲人成网站高清观看| 成人高潮视频无遮挡免费网站| 国产精品99久久久久久久久| 日本-黄色视频高清免费观看| 99久久精品热视频| 一区二区三区免费毛片| 观看免费一级毛片| 亚洲国产高清在线一区二区三| 国产精品av视频在线免费观看| 九色成人免费人妻av| 午夜老司机福利剧场| 国产又黄又爽又无遮挡在线| 亚洲经典国产精华液单| 欧美日韩在线观看h| 哪里可以看免费的av片| 亚洲五月天丁香| 午夜福利成人在线免费观看| 啦啦啦韩国在线观看视频| 一级黄色大片毛片| 国产午夜精品久久久久久一区二区三区 | 天堂影院成人在线观看| 免费av不卡在线播放| 又粗又爽又猛毛片免费看| 日韩中字成人| 99国产精品一区二区蜜桃av| 午夜激情福利司机影院| 国产精品人妻久久久久久| h日本视频在线播放| 一区二区三区免费毛片| 男人狂女人下面高潮的视频| 久久精品国产鲁丝片午夜精品| 亚洲国产日韩欧美精品在线观看| 大香蕉久久网| 久久热精品热| 午夜福利在线在线| 欧美性感艳星| 久久精品国产亚洲网站| 九九爱精品视频在线观看| 国产亚洲精品av在线| 麻豆久久精品国产亚洲av| 久久久久久久久大av| 少妇的逼好多水| 亚洲国产精品合色在线| 不卡视频在线观看欧美| 神马国产精品三级电影在线观看| 亚洲精品粉嫩美女一区| 嫩草影院新地址| 99热6这里只有精品| 成人欧美大片| 黄色一级大片看看| 美女cb高潮喷水在线观看| 国产亚洲av嫩草精品影院| 日本欧美国产在线视频| 精品久久久久久久久久久久久| 成人高潮视频无遮挡免费网站| 人人妻人人澡欧美一区二区| 乱人视频在线观看| 99热全是精品| 久久久久免费精品人妻一区二区| 欧美+亚洲+日韩+国产| 亚洲精品粉嫩美女一区| 日韩 亚洲 欧美在线| 亚洲av电影不卡..在线观看| 国国产精品蜜臀av免费| 69人妻影院| 一本一本综合久久| 国产69精品久久久久777片| 赤兔流量卡办理| 免费黄网站久久成人精品| av专区在线播放| 亚洲色图av天堂| 国产真实伦视频高清在线观看| www.色视频.com| 18禁在线播放成人免费| 美女 人体艺术 gogo| 成人欧美大片| 亚洲一区二区三区色噜噜| 丰满人妻一区二区三区视频av| 夜夜爽天天搞| 午夜福利成人在线免费观看| 亚洲成av人片在线播放无| 不卡一级毛片| 草草在线视频免费看| 国产一区二区亚洲精品在线观看| 免费观看的影片在线观看| 男人狂女人下面高潮的视频| 国产成人a区在线观看| 国产极品精品免费视频能看的| 国产成人91sexporn| 免费观看人在逋| 51国产日韩欧美| 成人漫画全彩无遮挡| 男女做爰动态图高潮gif福利片| 成人国产麻豆网| 午夜精品在线福利| 亚洲欧美精品综合久久99| 免费av观看视频| 波野结衣二区三区在线| 男人舔女人下体高潮全视频| 97在线视频观看| 男女啪啪激烈高潮av片| 欧美+亚洲+日韩+国产| 亚洲av成人av| 美女高潮的动态| 国产精品野战在线观看| 午夜福利在线观看免费完整高清在 | 在线天堂最新版资源| 亚洲精品亚洲一区二区| 国产精品日韩av在线免费观看| 亚洲高清免费不卡视频| 日本在线视频免费播放| 欧美日韩综合久久久久久| 国产大屁股一区二区在线视频| 啦啦啦啦在线视频资源| 天堂影院成人在线观看| 搡老熟女国产l中国老女人| 国产黄片美女视频| 最新中文字幕久久久久| 亚洲高清免费不卡视频| 直男gayav资源| 免费av不卡在线播放| 国产一区二区在线观看日韩| 嫩草影院新地址| 国产精品久久电影中文字幕| 女同久久另类99精品国产91| 不卡视频在线观看欧美| 亚洲av中文字字幕乱码综合| 99久久九九国产精品国产免费| 九九久久精品国产亚洲av麻豆| 综合色av麻豆| 欧美精品国产亚洲| 午夜福利视频1000在线观看| 日本一本二区三区精品| 日本-黄色视频高清免费观看| 亚洲av成人av| 男女视频在线观看网站免费| 一本精品99久久精品77| 看十八女毛片水多多多| 亚洲av成人av| 色视频www国产| 免费看光身美女| 美女大奶头视频| 女人十人毛片免费观看3o分钟| 老司机午夜福利在线观看视频| 国产aⅴ精品一区二区三区波| 精品无人区乱码1区二区| 国产在线男女| 久久国产乱子免费精品| 欧美激情在线99| 一本精品99久久精品77| 日韩欧美 国产精品| 18禁在线播放成人免费| 午夜福利成人在线免费观看| 免费av不卡在线播放| 日韩三级伦理在线观看| 秋霞在线观看毛片| 国产黄片美女视频| 久久中文看片网| 在现免费观看毛片| 18+在线观看网站| 永久网站在线| 国产欧美日韩精品亚洲av| 午夜激情欧美在线| 国产精品伦人一区二区| 免费观看人在逋| 别揉我奶头~嗯~啊~动态视频| 亚洲人成网站在线观看播放| 中文字幕av成人在线电影| 深爱激情五月婷婷| 欧美成人免费av一区二区三区| 搞女人的毛片| 成人三级黄色视频| 国产高清有码在线观看视频| 在线播放国产精品三级| 亚洲一区二区三区色噜噜| 淫秽高清视频在线观看| 晚上一个人看的免费电影| 一区二区三区高清视频在线| 麻豆成人午夜福利视频| 1000部很黄的大片| 好男人在线观看高清免费视频| 一级黄色大片毛片| 中国美女看黄片| 两性午夜刺激爽爽歪歪视频在线观看| 日韩欧美免费精品| 欧美+日韩+精品| 日韩强制内射视频| 美女被艹到高潮喷水动态| 在线免费观看的www视频| 校园春色视频在线观看| 久久精品久久久久久噜噜老黄 | 免费观看精品视频网站| 十八禁网站免费在线| 国内久久婷婷六月综合欲色啪| 欧美日韩在线观看h| 精品无人区乱码1区二区| 免费在线观看影片大全网站| 麻豆成人午夜福利视频| 欧美一区二区精品小视频在线| 毛片女人毛片| 国产精品久久久久久久电影| 国产精华一区二区三区| 亚洲国产高清在线一区二区三| 美女内射精品一级片tv| 亚洲av二区三区四区| 偷拍熟女少妇极品色| 久久久精品大字幕| 午夜影院日韩av| 男女边吃奶边做爰视频| 国产高潮美女av| 日本 av在线| 成人午夜高清在线视频| 国产精品一二三区在线看| 搞女人的毛片| 亚洲电影在线观看av| 亚洲欧美精品综合久久99| 一本精品99久久精品77| 国产精品不卡视频一区二区| 97人妻精品一区二区三区麻豆| 91在线精品国自产拍蜜月| 99热这里只有精品一区| 在线观看66精品国产| 寂寞人妻少妇视频99o| 亚洲三级黄色毛片| 91av网一区二区| 能在线免费观看的黄片| 天天一区二区日本电影三级| 免费高清视频大片| 97超级碰碰碰精品色视频在线观看| 国产白丝娇喘喷水9色精品| 亚洲国产欧美人成| 亚洲不卡免费看| 国产爱豆传媒在线观看| 99久久无色码亚洲精品果冻| 日韩欧美在线乱码| 男女做爰动态图高潮gif福利片| 免费在线观看成人毛片| 亚洲专区国产一区二区| 观看美女的网站| 久久精品国产鲁丝片午夜精品| 99久久久亚洲精品蜜臀av| 美女内射精品一级片tv| 亚洲,欧美,日韩| 国产一区二区三区在线臀色熟女| 啦啦啦韩国在线观看视频| 最新中文字幕久久久久| 免费在线观看成人毛片| 永久网站在线| 男女视频在线观看网站免费| 国产成人一区二区在线| 天堂√8在线中文| 精品久久久久久成人av| 99国产极品粉嫩在线观看| 久久午夜亚洲精品久久| 久久鲁丝午夜福利片| 免费不卡的大黄色大毛片视频在线观看 | 男女边吃奶边做爰视频| 亚洲人成网站在线播| 91久久精品电影网| 男女视频在线观看网站免费| 亚洲国产日韩欧美精品在线观看| 蜜臀久久99精品久久宅男| 午夜a级毛片| 国产精品一区www在线观看| 成人高潮视频无遮挡免费网站| 欧美性猛交黑人性爽| 亚洲欧美日韩高清在线视频| 亚洲综合色惰| 欧洲精品卡2卡3卡4卡5卡区| 免费看a级黄色片| 午夜激情福利司机影院| 亚洲中文字幕一区二区三区有码在线看| 亚洲七黄色美女视频| 国产亚洲精品av在线| 草草在线视频免费看| 女的被弄到高潮叫床怎么办| 男人舔奶头视频| 国产美女午夜福利| 三级国产精品欧美在线观看| 18禁裸乳无遮挡免费网站照片| 久久精品国产亚洲av涩爱 | 看非洲黑人一级黄片| 99热这里只有是精品在线观看| 悠悠久久av| 国模一区二区三区四区视频| 国产精品一区二区三区四区免费观看 | 熟女人妻精品中文字幕| 高清午夜精品一区二区三区 | 国产人妻一区二区三区在| 亚洲成人av在线免费| 91在线精品国自产拍蜜月| 国产精品99久久久久久久久| 免费高清视频大片| 亚洲av成人精品一区久久| 国产精品一及| 色尼玛亚洲综合影院| 搡老岳熟女国产| 亚洲熟妇中文字幕五十中出| 亚洲精品国产成人久久av| 成人av在线播放网站| 欧美激情国产日韩精品一区| 日韩欧美三级三区| 国产视频内射| 一级黄色大片毛片| 69av精品久久久久久| 观看美女的网站| 久久久久性生活片| 国产极品精品免费视频能看的| 性色avwww在线观看| 婷婷精品国产亚洲av| 日韩欧美免费精品| 波多野结衣巨乳人妻| 精品久久久久久久久av| 亚洲天堂国产精品一区在线| 成人三级黄色视频| 桃色一区二区三区在线观看| 国产视频一区二区在线看| 毛片一级片免费看久久久久| 中国美白少妇内射xxxbb| АⅤ资源中文在线天堂| 亚洲人成网站高清观看| 久久婷婷人人爽人人干人人爱| 免费看日本二区| 我的老师免费观看完整版| 97热精品久久久久久| 亚洲色图av天堂| а√天堂www在线а√下载| 久久久久免费精品人妻一区二区| 日韩强制内射视频| 人妻丰满熟妇av一区二区三区| 国产精品免费一区二区三区在线| 久久精品国产99精品国产亚洲性色| 国产探花极品一区二区| 极品教师在线视频| 日本五十路高清| 国产美女午夜福利| 欧美一区二区精品小视频在线| 成人性生交大片免费视频hd| 亚洲欧美精品自产自拍| 在现免费观看毛片| 你懂的网址亚洲精品在线观看 | 黄色一级大片看看| 亚洲av熟女| 久久久久国产网址| 免费观看人在逋| 99精品在免费线老司机午夜| 免费av毛片视频| 99热这里只有精品一区| 精品久久久久久久久亚洲| 啦啦啦啦在线视频资源| 99国产极品粉嫩在线观看| 在线免费十八禁| 精品一区二区三区视频在线观看免费| 在线观看午夜福利视频| 少妇熟女aⅴ在线视频| 深爱激情五月婷婷| 熟妇人妻久久中文字幕3abv| 国产黄色小视频在线观看| 看免费成人av毛片| 观看美女的网站| 又黄又爽又免费观看的视频| 两个人视频免费观看高清| 国内精品一区二区在线观看| а√天堂www在线а√下载| 日韩,欧美,国产一区二区三区 | 国产高潮美女av| 日本欧美国产在线视频| 女同久久另类99精品国产91| 国产精品久久久久久av不卡| 看免费成人av毛片| 一个人免费在线观看电影| 成人欧美大片| 97人妻精品一区二区三区麻豆| 一个人观看的视频www高清免费观看| 人人妻人人看人人澡| 十八禁网站免费在线| 午夜精品一区二区三区免费看| avwww免费| 少妇熟女欧美另类| 精品一区二区三区人妻视频| 免费观看的影片在线观看| 亚洲av五月六月丁香网| 国产一区二区三区在线臀色熟女| 日韩一区二区视频免费看| 亚洲熟妇中文字幕五十中出| av在线老鸭窝| 我的女老师完整版在线观看| 国产精品久久电影中文字幕| 精华霜和精华液先用哪个| 欧美色欧美亚洲另类二区| 91午夜精品亚洲一区二区三区| 好男人在线观看高清免费视频| 插逼视频在线观看| 日韩高清综合在线| 变态另类成人亚洲欧美熟女| 日韩av在线大香蕉| 欧美高清成人免费视频www| 国产一区二区亚洲精品在线观看| 国内精品久久久久精免费| 我的女老师完整版在线观看| 国内久久婷婷六月综合欲色啪| 香蕉av资源在线| 国产精品一区二区免费欧美| 精品午夜福利在线看| 国产探花在线观看一区二区| 国产成人福利小说| 日韩欧美在线乱码| 日本a在线网址| 老熟妇乱子伦视频在线观看| 一区二区三区免费毛片| 亚洲欧美精品自产自拍| 国产精品美女特级片免费视频播放器| 人人妻人人澡欧美一区二区| 欧美成人精品欧美一级黄| 日本 av在线| 亚洲图色成人| 免费电影在线观看免费观看| 亚洲在线观看片| 天堂动漫精品| 精品久久久噜噜| 欧美一区二区国产精品久久精品| 亚洲成av人片在线播放无| 久久精品影院6| 国产欧美日韩精品一区二区| 91久久精品国产一区二区三区| 日日摸夜夜添夜夜添av毛片| 午夜亚洲福利在线播放| 亚洲欧美日韩东京热| 夜夜爽天天搞| 男人狂女人下面高潮的视频| 在线a可以看的网站| 欧美成人一区二区免费高清观看| 欧美区成人在线视频| 99久国产av精品国产电影| 老熟妇仑乱视频hdxx| 99热6这里只有精品| 91av网一区二区| 国产精品精品国产色婷婷| 午夜视频国产福利| 欧美精品国产亚洲| 精品一区二区三区av网在线观看| 国产一区二区亚洲精品在线观看| 欧美日本视频| 天美传媒精品一区二区| 内射极品少妇av片p| 在线观看免费视频日本深夜| 91久久精品国产一区二区成人| 婷婷精品国产亚洲av| 美女内射精品一级片tv| 99视频精品全部免费 在线| 久久久久九九精品影院| 成人av一区二区三区在线看| 最近2019中文字幕mv第一页| 国产午夜精品久久久久久一区二区三区 | 色在线成人网| 国产欧美日韩精品亚洲av| 午夜精品一区二区三区免费看| 国产精品无大码| 欧美中文日本在线观看视频| 国产高清激情床上av| 乱码一卡2卡4卡精品| 日韩精品青青久久久久久| a级毛片免费高清观看在线播放| 97超碰精品成人国产| 观看美女的网站| 亚洲欧美日韩高清在线视频| 国产午夜精品论理片| 国产亚洲91精品色在线| 亚洲人成网站在线观看播放| 久久婷婷人人爽人人干人人爱| 美女xxoo啪啪120秒动态图| 中文资源天堂在线| 日本免费a在线| 人妻制服诱惑在线中文字幕| 91久久精品国产一区二区成人| 免费av不卡在线播放| 毛片女人毛片| 久久热精品热| 97超级碰碰碰精品色视频在线观看| 国产午夜精品久久久久久一区二区三区 | 日本a在线网址| 免费av不卡在线播放| 亚洲最大成人手机在线| 在线观看av片永久免费下载| 97超级碰碰碰精品色视频在线观看| 日本色播在线视频| 九九爱精品视频在线观看| 丰满的人妻完整版| 免费人成在线观看视频色| 岛国在线免费视频观看| 夜夜看夜夜爽夜夜摸| 狂野欧美白嫩少妇大欣赏| 亚洲欧美日韩无卡精品|