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

    LNG船靠離泊操縱運動數(shù)值模擬分析

    2017-11-01 12:46:38,,,,
    船海工程 2017年5期
    關(guān)鍵詞:船舶

    ,,,,

    (1.中海油能源發(fā)展股份有限公司采油服務(wù)分公司,天津 300452;2.上海船舶研究設(shè)計院,上海 201203)

    LNG船靠離泊操縱運動數(shù)值模擬分析

    鄭坤1,時光志1,張雨新2,周毅1,陸晟2

    (1.中海油能源發(fā)展股份有限公司采油服務(wù)分公司,天津 300452;2.上海船舶研究設(shè)計院,上海 201203)

    基于CFD方法對廣西防城港內(nèi)小型LNG船自力靠離泊操縱運動進(jìn)行數(shù)值模擬,分析風(fēng)、浪、流作用下靠離泊過程中船體的運動響應(yīng)。數(shù)值計算結(jié)果表明:該型LNG船能夠通過自力實現(xiàn)靠離泊操縱。靠泊運動中,當(dāng)受到離岸方向風(fēng)、浪、流作用時,可以在轉(zhuǎn)艏過程中通過右槳產(chǎn)生推力,減小船體漂離碼頭的距離,當(dāng)風(fēng)、浪、流為攏岸方向時,可以將右舵轉(zhuǎn)至180°,使右槳產(chǎn)生拉力,以防止船體向碼頭移動時速度過大。在離泊時,為了安全考慮,在攏岸方向風(fēng)、浪、流作用時,船舶在轉(zhuǎn)艏前需要與岸邊保持較大距離。

    靠離泊;操縱運動;數(shù)值模擬;CFD

    船舶操縱運動的數(shù)值模擬方法大多基于上世紀(jì)70年代日本船舶操縱運動數(shù)學(xué)模型研究小組(mathematical modeling group,MMG)開發(fā)的分離式模型,該計算模型將螺旋槳推力、舵力、錨鏈力通過水動力系數(shù)加入到船舶運動方程中,同時將環(huán)境載荷,如風(fēng)、浪、流對船體作用力,通過一系列系數(shù)添加到運動方程中。MMG模型比較簡單實用,因此在船舶操縱運動問題中已有較多的應(yīng)用[1-3]。但如何比較準(zhǔn)確地獲得MMG模型中各個水動力系數(shù)一直是限制MMG模型應(yīng)用的難點。近年來,計算流體力學(xué)(computational fluid dynamics, CFD)為操縱性數(shù)值模擬提供了另一種更加精細(xì)地描述船舶操縱運動的預(yù)報工具。CFD直接對船體運動流場進(jìn)行求解,不需要做過多的簡化,也不需要通過經(jīng)驗公式或者試驗來確定水動力系數(shù),能夠比較真實地模擬船、槳、舵之間流動的相互干擾過程,因此,在一些復(fù)雜的操縱性問題中得到較廣泛的應(yīng)用,如多自由度回轉(zhuǎn)[4]、限制水域的靠離泊[5]、波浪中橫甩問題[6]等。

    擬對一小型LNG船在廣西防城港漁澫港區(qū)內(nèi)的靠離泊問題進(jìn)行數(shù)值模擬研究,該船為3萬m3LNG船,為國內(nèi)首艘搭載全回轉(zhuǎn)電力推進(jìn)系統(tǒng)的小型LNG船,由于該船頻繁進(jìn)出防城港,其運營模式將對國內(nèi)中小型LNG的水上運輸有重要的參考價值[7-8]??紤]到該船操縱性良好的特點,為降低航運成本,對該船在港區(qū)內(nèi)的依靠自力實現(xiàn)靠離泊問題進(jìn)行分析。船體運動過程基于CFD進(jìn)行計算,通過對不同環(huán)境條件下靠泊和離泊過程進(jìn)行模擬,分析3萬m3LNG船在防城港內(nèi)的靠離泊能力,為防城港駕引人員的實際操作提供參考。

    1 工作水域概況

    圖1 廣西防城港漁澫港區(qū)相關(guān)布置示意

    工作水域為廣西防城港漁澫港區(qū),見圖1。目前,港區(qū)已建成泊位40個,其中萬t級以上泊位25個,5萬t級以上泊位15個,最大設(shè)計靠泊能力20萬t級,碼頭岸線總長7.745 km,年通過能力3 701萬t。分析選擇的泊位為第5作業(yè)區(qū)的501號泊位,位于漁澫港區(qū)第五作業(yè)區(qū)的南部,北側(cè)為規(guī)劃的漁澫港區(qū)第五作業(yè)區(qū)502號泊位,南側(cè)為已建成的廣西東灣5 萬t級液體化工碼頭,與中電防城港電廠隔江相望。碼頭按8萬m3LNG船和5萬m3LPG船設(shè)計,兼靠3萬m3LNG船。營運貨種為LNG和LPG等液體散貨[9]。

    回旋水域布置在碼頭前方,寬度為80 m?;匦虬?萬m3LNG 船和5萬m3LPG 船設(shè)計,由于碼頭前方有一條規(guī)劃航道和一條防城港電廠煤碼頭航道,為了減少回旋水域?qū)降赖挠绊?,回旋水域采用橢圓形布置,長軸取為600 m,短軸取為480 m??紤]到進(jìn)港航道水深現(xiàn)狀,回旋水域底標(biāo)高按照現(xiàn)有航道底標(biāo)高取值為-10 m。回旋水域尺度可滿足3萬m3LNG船不乘潮進(jìn)港,同時滿足5萬m3LPG船乘潮不滿載進(jìn)港。

    2 船型及計算工況

    3萬m3LNG船的船型主尺度見表1。船體幾何形狀見圖2,艏部裝有1個艏側(cè)推進(jìn)器,側(cè)推孔直徑為2.25 m,艏側(cè)推最大推力為196 kN,功率為1 500 kW。為簡化計算,在數(shù)值模擬中沒有建立艏側(cè)推的螺旋槳模型,而是將側(cè)推力直接施加到船體上。艉部為雙吊艙全回轉(zhuǎn)推進(jìn)器,螺旋

    槳為5葉槳定距槳,直徑為3.7 m,無導(dǎo)流管。舵的最大回轉(zhuǎn)速度為2.0 r/min。

    表1 船型主尺度 m

    圖2 船體幾何形狀(鍘視)

    根據(jù)港區(qū)自然條件,確定4個模擬工況,風(fēng)速均為13.8 m/s,對應(yīng)蒲氏6級風(fēng),具體環(huán)境參數(shù)見表2。分別計算每種工況滿載靠泊和壓載離泊2種情況。

    表2 計算工況

    3 數(shù)值模擬方法

    船體運動和螺旋槳計算采用軟件商業(yè)CFD(computational fluid dynamics)軟件STAR-CCM+進(jìn)行。由于操縱過程中船體運動范圍較大,因此,采用重疊網(wǎng)格技術(shù)(over set)處理船體運動。為了減小計算量,計算中將螺旋槳簡化為虛擬盤。虛擬盤法根據(jù)給定的螺旋槳敞水?dāng)?shù)據(jù)對螺旋槳盤附近的水流施加體積力,從而產(chǎn)生推力。該方法在船、槳、舵相互影響的研究方面已有較多應(yīng)用[10-11],并表現(xiàn)出較高計算效率和令人滿意的計算精度,是研究船、槳相互干擾問題中一種常用的處理方法。此外,由于港區(qū)內(nèi)波浪較小,難以在CFD模擬中精確模擬波浪的演化,因此,波浪力是通過勢流軟件Hydrostar進(jìn)行計算得到的,然后帶入到CFD求解過程中,在CFD計算過程中根據(jù)船體與波浪的角度,插值得到每個時間步船體受到的波浪力。這種處理方法的主要優(yōu)點是能夠較為合理地簡化計算問題,在基于MMG模型的數(shù)值模擬中應(yīng)用較多。Skejic和MFaltinsen[12]通過研究認(rèn)為船舶在波浪中的慢漂運動主要是由二階波浪力引起的,一階波浪力對其操縱運動影響并不顯著,因此,在靠離泊計算時只考慮二階波浪力。由于船舶靠離泊過程的主要特點是低速、大漂角,因此計算時船舶運動主要考慮縱向和橫向平移運動以及艏搖3個自由度。

    3.1 控制方程

    計算中控制方程為Reynolds平均Navier-Stokes方程,即RANS方程,湍流模型采用SSTk-ω兩方程湍流模型,自由面基于VOF(volume of fluid)方法求解。對流項離散采用二階精度的迎風(fēng)格式。

    控制方程包括連續(xù)性方程和動量守恒方程,對于不可壓縮流體為

    (1)

    (2)

    式中:i,j為變量在坐標(biāo)系中第i和j方向的分量;p為壓力,u為速度矢量,ρ為流體密度,f為質(zhì)量力。

    SSTk-ω模型湍動能和湍動耗散率的求解方程為

    (3)

    (4)

    VOF控制方程為

    (5)

    式中:φ為流體相體積分?jǐn)?shù)。

    3.2 計算域及網(wǎng)格劃分

    計算域取為長方體,大小為:船前向前延伸1.0Lpp,2側(cè)和底部向外延伸1.5Lpp。底部為壁面邊界條件,4周各面為來流和壓力出口,具體設(shè)置根據(jù)各個工況實際水流方向確定。計算網(wǎng)格類型為非結(jié)構(gòu)6面體網(wǎng)格,采用STAR-CCM+中自帶網(wǎng)格生成器Trimmer生成,計算網(wǎng)格見圖3。

    圖3 計算域網(wǎng)格

    為更好地計算船體附近的流動特征,在船體表面生成了邊界層網(wǎng)格,同時對船體表面以及船底與水底之間流場進(jìn)行了網(wǎng)格加密,見圖4。

    圖4 船體表面網(wǎng)格

    為了更好地捕捉水面興波,對靜水面附近網(wǎng)格在高度方向進(jìn)行加密。為解決船體大幅度運動帶來的動邊界問題,采用重疊網(wǎng)格進(jìn)行處理,計算中生成2套網(wǎng)格:船體附近區(qū)域和背景網(wǎng)格。船體附近網(wǎng)格隨船體做浮體運動,背景網(wǎng)格僅作平移運動,平移速度等于船體的移動速度。這樣既保證了船體附近的網(wǎng)格質(zhì)量,同時避免使用過大的計算域。此外,為保證船體旋轉(zhuǎn)后周圍流場的計算精度,在船體旋轉(zhuǎn)區(qū)域范圍內(nèi)對網(wǎng)格進(jìn)行加密。計算網(wǎng)格總數(shù)為330萬,其中船體部分網(wǎng)格為185萬,背景區(qū)域網(wǎng)格數(shù)為145萬。

    4 計算結(jié)果及分析

    4.1 靠泊過程數(shù)值模擬

    由圖1可見,船舶靠泊操作過程大體程序是先將船駛?cè)牖剞D(zhuǎn)區(qū)中心附近,然后開始向左舷側(cè)做轉(zhuǎn)艏操作,當(dāng)轉(zhuǎn)至船體舷側(cè)與碼頭平行時,再通過動力定位能力逐漸平移至碼頭,實現(xiàn)靠泊。數(shù)值模擬主要關(guān)注回轉(zhuǎn)區(qū)調(diào)頭過程。計算初始時刻假定船體已駛?cè)牖剞D(zhuǎn)區(qū),航速為0,左舵已轉(zhuǎn)至90°位置,左側(cè)螺旋槳開始工作,轉(zhuǎn)速為90 r/min。右槳和右舵的操作需要根據(jù)港內(nèi)風(fēng)、浪、流的方向和大小而定,在實際靠泊過程中駕駛員依靠經(jīng)驗結(jié)合港區(qū)環(huán)境和船體運動速度進(jìn)行實時調(diào)整。這種交互式的操縱策略難以數(shù)值模擬,因此采用一種簡化的操縱策略??紤]到工況1和工況3中,風(fēng)、浪、流為離岸方向,在回轉(zhuǎn)過程中需要右槳提供推力,因此右舵為0°,右槳在船頭指向碼頭方向時開始工作,將船體推向碼頭,具體工作參數(shù)見表3。

    表3 靠泊過程中右側(cè)螺旋槳工作規(guī)律

    相反,在工況2和工況4中,風(fēng)、浪、流為攏岸方向,這時需要右槳提供拉力幫助船體定位,因此右舵處于180°,螺旋槳產(chǎn)生向后的拉力防止船體向碼頭方向產(chǎn)生過大的運動速度。需要說明的是,這種操縱控制策略不是最優(yōu)的,但通過這種簡單的操縱控制得到的船體靠泊運動響應(yīng)可以為實際駕引人員提供參考。

    船體的運動軌跡如圖5所示。

    圖5 靠泊時船體運動軌跡(時間間隔為30 s,虛線為示意軌跡)

    由圖5可見,船舶能夠在回轉(zhuǎn)區(qū)通過自力完成回轉(zhuǎn)運動,船體移動距離較小。風(fēng)、浪、流對船體艏搖角速度影響不大,對船體的位移有一定的影響。在工況1和工況2中,通過右槳的配合,船體轉(zhuǎn)艏后所處位置較為合理。在工況4中,由于右槳工作時間較短,船體轉(zhuǎn)艏后向回轉(zhuǎn)區(qū)西側(cè)偏移一定距離,給后面的動力定位平移至碼頭增加了難度,因此實際操作中,應(yīng)當(dāng)較早啟動右槳。

    4.2 離泊過程數(shù)值模擬

    由于碼頭前面水域較小,航道中心線離碼頭橫向距離較大,因此船舶離港時需要先通過艏側(cè)推和尾推將船體平移至離碼頭一段距離,然后開始轉(zhuǎn)艏操作。這里主要模擬轉(zhuǎn)艏過程中船體的運動情況。由于港口航道與碼頭的角度約為72°,所需回轉(zhuǎn)角度較靠泊時小,因此,在計算離泊時沒有開啟右槳,僅通過左槳和艏側(cè)推實現(xiàn)離泊操作。計算得到的船體運動軌跡見圖6。

    圖6 離泊時船體運動軌跡(時間間隔為15 s,虛線為示意軌跡)

    由圖6中可見,水流對船體操縱運動影響較大,不同工況下船體偏移距離有較大區(qū)別,在工況1和工況3中水流向離岸方向,船體不需要離泊較遠(yuǎn)即可進(jìn)行轉(zhuǎn)艏操作,但對于攏岸風(fēng)、浪、流作用時,如工況4,風(fēng)浪流向船施加碼頭方向力,因此船體需要平行離開碼頭較遠(yuǎn)距離再開始轉(zhuǎn)艏操作。整體上看,船舶能夠依靠自力安全離泊。

    5 結(jié)論

    從數(shù)值計算結(jié)果看,在右槳和右舵的輔助配合下,船舶能夠比較好地完成靠離泊運動。相對來講,靠泊過程所需時間較長,受環(huán)境影響較大。為更好地控制船體與碼頭之間的橫向距離,當(dāng)風(fēng)、浪、流為離岸方向時,可以通過右槳產(chǎn)生推力,減小船體漂離碼頭的距離;當(dāng)風(fēng)、浪、流為攏岸方向時,可以將右舵轉(zhuǎn)至180°,使右槳產(chǎn)生拉力,以防止船體向碼頭移動時速度過大。在離泊操縱過程中,船體轉(zhuǎn)艏前需要與碼頭保持一定的安全距離,在攏岸方向風(fēng)、浪、流作用時,所需的安全距離較大。

    [1] 楊鹽生.船舶靠離泊操縱數(shù)學(xué)模型的研究[J].大連海事大學(xué)學(xué)報,1996,22(4):11-15.

    [2] 朱冬健,馬寧,顧解忡,等.波浪中船舶操縱性數(shù)值預(yù)報及自航模驗證[J].中國艦船,2015,10(1):76-82.

    [3] 史宏達(dá),范天會,韓桂萍.東營港擴建工程碼頭船舶靠離泊數(shù)值模擬研究[J].中國海洋大學(xué)學(xué)報,2007,37(5):841-846.

    [4] WANG J H, ZhAO W W, WAN D C. Free maneuvering simulation of ONR tumblehome using overset grid method in naoe-FOAM-SJTU solver[C]. 31stSymposium on Naval Hydrodynamics, California, 2016:11-16.

    [5] WANG H M, ZOU Z J, TIAN X M. Computation of theviscous hydrodynamicforces on a KVLCC2 model moving obliquely in shallow water[J]. Journal of shanghai jiaotong university (science), 2009,14(2):241-244.

    [6] CARRICA P M, SADAT-HISSEINI H, STERM F. CFD analysis of broaching for a model surface combatant with explicit simulation of moving rudders and rotating propellers[J]. Computers & fluids, 2012,53:117-132.

    [7] 時光志,周志峰.某液化天然氣運輸船推進(jìn)器選型研究[J].船海工程,2014,43(5):132-137.

    [8] 何金平.中小型LNG運輸船推進(jìn)系統(tǒng)選型分析[J].船海工程,2014,43(5):142-144.

    [9] 趙斌.大型LNG船舶安全進(jìn)港靠泊防城港東灣液體化工碼頭的思考[J].港口經(jīng)濟,2015(7):59-61.

    [10] 吳召華,陳作鋼,代燚.基于體積力法的船體自航性能數(shù)值預(yù)報[J].上海交通大學(xué)學(xué)報,2013,47(6):943-949.

    [11] SIMONSEN C D, STERM F. RANS maneuvering simulation of essoosaka with rudder and a body-force propeller[J]. Journal of ship research, 2005,49(2):98-120.

    [12] SKEJJC R, FALTINSEN O M. A unified seakeeping and maneuvering analysis of a monohull in Regular IncidentWaves[C]. Proc. of the 7th International Conference on Hyd., Ischia, Italy, 2006:97-104.

    Numerical Investigation of Maneuvering Behavior of a LNG Carrier during Berthing and Departing

    ZHENGKun1,SHIGuang-zhi1,ZHANGYu-xin2,ZHOUYi1,LUSheng2

    (1.CNOOC Energy Technology & Services-Oil Production Services Company, Tianjin 300452, China;2.Shanghai Merchant Ship Design & Research Institute, Shanghai 201203, China)

    Numerical simulations of berthing and departure of a small LNG ship were carried out based on CFD method. The maneuvering behaviors of ship interacting with wind, wave and tide, during berthing and departing were analyzed. Numerical results showed that the computed ship is capable of berthing and departing by itself. When direction of wind, wave and tide is offshore, right side propeller should work, imposing thrust on the ship. The right side propeller needs pull the ship to prevent a large velocity of ship when the wind, wave and tide come alongside shore. During departing, ship should move offshore a large distance before turning when the wind, wave and tide direct shore.

    berthing and departing; maneuvering; numerical simulation; CFD

    U661.33

    A

    1671-7953(2017)05-0048-05

    10.3963/j.issn.1671-7953.2017.05.013

    2017-07-12

    修回日期:2017-08-31

    鄭坤(1987—),男,工學(xué)碩士,助理工程師

    研究方向:船舶與海洋結(jié)構(gòu)物設(shè)計制造

    猜你喜歡
    船舶
    船舶避碰路徑模糊控制系統(tǒng)
    計算流體力學(xué)在船舶操縱運動仿真中的應(yīng)用
    CM節(jié)點控制在船舶上的應(yīng)用
    基于改進(jìn)譜分析法的船舶疲勞強度直接計算
    《船舶》2022 年度征訂啟事
    船舶(2021年4期)2021-09-07 17:32:22
    船舶!請加速
    BOG壓縮機在小型LNG船舶上的應(yīng)用
    船舶 揚帆奮起
    軍工文化(2017年12期)2017-07-17 06:08:06
    船舶壓載水管理系統(tǒng)
    中國船檢(2017年3期)2017-05-18 11:33:09
    小型船舶艉軸架設(shè)計
    船海工程(2015年4期)2016-01-05 15:53:30
    在线精品无人区一区二区三| 亚洲精华国产精华液的使用体验| 国产乱人偷精品视频| 高清毛片免费看| 久久久久精品久久久久真实原创| 嘟嘟电影网在线观看| 成人国产av品久久久| 在线观看免费日韩欧美大片 | 老司机影院成人| 美女xxoo啪啪120秒动态图| 免费观看无遮挡的男女| 精品久久久久久久久av| 美女福利国产在线| 日韩成人av中文字幕在线观看| 一级毛片我不卡| 日韩亚洲欧美综合| 性色av一级| 免费不卡的大黄色大毛片视频在线观看| 80岁老熟妇乱子伦牲交| 丝袜在线中文字幕| 人妻系列 视频| 久久久久久久国产电影| av黄色大香蕉| 大又大粗又爽又黄少妇毛片口| 精品人妻熟女毛片av久久网站| 极品少妇高潮喷水抽搐| 国产综合精华液| 自拍欧美九色日韩亚洲蝌蚪91 | 亚洲av在线观看美女高潮| 建设人人有责人人尽责人人享有的| 婷婷色av中文字幕| 80岁老熟妇乱子伦牲交| 一个人免费看片子| 少妇熟女欧美另类| 国产精品99久久久久久久久| 国产伦理片在线播放av一区| 午夜福利网站1000一区二区三区| 国产一区二区在线观看日韩| 亚洲精品视频女| 97精品久久久久久久久久精品| 国产亚洲91精品色在线| 久久久国产一区二区| 26uuu在线亚洲综合色| 亚洲av中文av极速乱| 伦理电影大哥的女人| 久久久国产欧美日韩av| 51国产日韩欧美| 日日摸夜夜添夜夜爱| 校园人妻丝袜中文字幕| 久久人人爽av亚洲精品天堂| 极品人妻少妇av视频| 久久久久视频综合| 日韩成人av中文字幕在线观看| 水蜜桃什么品种好| 狠狠精品人妻久久久久久综合| 中文在线观看免费www的网站| 色5月婷婷丁香| 日本黄大片高清| 国产 一区精品| av免费在线看不卡| 啦啦啦在线观看免费高清www| 多毛熟女@视频| 美女xxoo啪啪120秒动态图| 日日爽夜夜爽网站| 99久久综合免费| 亚洲av免费高清在线观看| 97精品久久久久久久久久精品| 亚洲国产成人一精品久久久| 亚洲国产精品999| 免费看光身美女| 国产成人精品久久久久久| 观看免费一级毛片| 欧美日韩视频高清一区二区三区二| 你懂的网址亚洲精品在线观看| 亚洲av成人精品一区久久| 99热全是精品| 国产成人精品婷婷| 久久精品国产自在天天线| 久久婷婷青草| 观看免费一级毛片| 免费少妇av软件| 嘟嘟电影网在线观看| 国产亚洲精品久久久com| 国产精品人妻久久久影院| 欧美日韩亚洲高清精品| 中国三级夫妇交换| 插逼视频在线观看| 亚洲中文av在线| 蜜桃在线观看..| 夫妻午夜视频| 国产在线男女| av在线app专区| 国语对白做爰xxxⅹ性视频网站| 91久久精品电影网| 久久狼人影院| 欧美 亚洲 国产 日韩一| 日韩强制内射视频| 黑人猛操日本美女一级片| 亚洲伊人久久精品综合| 国产亚洲5aaaaa淫片| 日本av免费视频播放| 欧美高清成人免费视频www| 国产精品一区二区三区四区免费观看| 欧美精品一区二区免费开放| 各种免费的搞黄视频| 国产亚洲午夜精品一区二区久久| 人人妻人人添人人爽欧美一区卜| 亚洲精品日本国产第一区| 亚洲av不卡在线观看| 亚洲人成网站在线观看播放| 亚洲欧美一区二区三区黑人 | 哪个播放器可以免费观看大片| 在线观看免费日韩欧美大片 | 2018国产大陆天天弄谢| 黄色毛片三级朝国网站 | 内射极品少妇av片p| 中文天堂在线官网| 国产一区亚洲一区在线观看| 热99国产精品久久久久久7| 免费黄色在线免费观看| 大又大粗又爽又黄少妇毛片口| 国产有黄有色有爽视频| 欧美老熟妇乱子伦牲交| 欧美日本中文国产一区发布| tube8黄色片| 丝袜在线中文字幕| 久久婷婷青草| 国产在线免费精品| 精品一品国产午夜福利视频| 国产av精品麻豆| 汤姆久久久久久久影院中文字幕| 多毛熟女@视频| 亚洲成人手机| 午夜精品国产一区二区电影| 欧美激情国产日韩精品一区| 观看av在线不卡| 我的老师免费观看完整版| 国产精品秋霞免费鲁丝片| 亚洲自偷自拍三级| 又粗又硬又长又爽又黄的视频| 91久久精品国产一区二区三区| 曰老女人黄片| 三上悠亚av全集在线观看 | 免费黄频网站在线观看国产| 午夜福利,免费看| 一个人看视频在线观看www免费| 亚洲精品日韩在线中文字幕| 国产一区二区三区av在线| 亚洲av.av天堂| 亚洲高清免费不卡视频| 国产深夜福利视频在线观看| 中国美白少妇内射xxxbb| 男女啪啪激烈高潮av片| 欧美日韩视频精品一区| 国产精品三级大全| 亚洲精品成人av观看孕妇| 人人妻人人添人人爽欧美一区卜| 各种免费的搞黄视频| 欧美日韩精品成人综合77777| 日韩中文字幕视频在线看片| videossex国产| 国产一级毛片在线| 91精品一卡2卡3卡4卡| 最近中文字幕2019免费版| 国产午夜精品久久久久久一区二区三区| 有码 亚洲区| 中文字幕av电影在线播放| 高清不卡的av网站| av在线播放精品| a 毛片基地| 亚洲国产欧美日韩在线播放 | 中文欧美无线码| 91久久精品国产一区二区三区| 美女福利国产在线| 欧美 亚洲 国产 日韩一| 亚洲国产最新在线播放| 久久人人爽av亚洲精品天堂| 99久久中文字幕三级久久日本| 亚洲国产日韩一区二区| 午夜福利影视在线免费观看| 汤姆久久久久久久影院中文字幕| 国产成人aa在线观看| 国产老妇伦熟女老妇高清| 欧美xxxx性猛交bbbb| 曰老女人黄片| 国语对白做爰xxxⅹ性视频网站| 黑人猛操日本美女一级片| 97超视频在线观看视频| 亚洲内射少妇av| 久久久久视频综合| 少妇被粗大猛烈的视频| 最新中文字幕久久久久| 日本vs欧美在线观看视频 | 在线观看免费日韩欧美大片 | 亚洲无线观看免费| 在线精品无人区一区二区三| 午夜免费男女啪啪视频观看| 亚洲在久久综合| 亚洲激情五月婷婷啪啪| 国产日韩一区二区三区精品不卡 | 精品久久国产蜜桃| 91精品国产九色| 欧美老熟妇乱子伦牲交| 久久久久久久久久人人人人人人| 2022亚洲国产成人精品| 亚洲精品国产成人久久av| 日韩精品有码人妻一区| 精品久久久久久电影网| 在线天堂最新版资源| 国产成人免费观看mmmm| 亚洲欧洲精品一区二区精品久久久 | 成人免费观看视频高清| 国产伦精品一区二区三区四那| 免费看光身美女| 国产成人午夜福利电影在线观看| 国产精品.久久久| 伊人亚洲综合成人网| 久久久久久久久久久免费av| 97在线人人人人妻| 国产无遮挡羞羞视频在线观看| 亚洲精品一区蜜桃| 如何舔出高潮| 久久青草综合色| 国产精品嫩草影院av在线观看| av线在线观看网站| 亚洲欧美日韩另类电影网站| 国产亚洲午夜精品一区二区久久| 熟女电影av网| 国产在线男女| 国产视频内射| 色视频www国产| 精品久久久噜噜| 22中文网久久字幕| 亚洲伊人久久精品综合| 精品一区二区免费观看| 一区二区av电影网| 少妇被粗大猛烈的视频| 少妇人妻一区二区三区视频| 卡戴珊不雅视频在线播放| 熟妇人妻不卡中文字幕| 九草在线视频观看| 亚州av有码| 国产精品99久久99久久久不卡 | 国产精品不卡视频一区二区| 久久久久国产网址| 亚洲激情五月婷婷啪啪| 欧美bdsm另类| 美女国产视频在线观看| 99九九线精品视频在线观看视频| 91午夜精品亚洲一区二区三区| 99热网站在线观看| 91久久精品国产一区二区三区| 嫩草影院入口| 青春草视频在线免费观看| 久久久午夜欧美精品| av国产精品久久久久影院| 老司机亚洲免费影院| 国产在线男女| 看非洲黑人一级黄片| 成人无遮挡网站| 老女人水多毛片| 亚洲不卡免费看| 99九九在线精品视频 | 亚洲四区av| a级一级毛片免费在线观看| 国产在线视频一区二区| 精品久久久久久久久亚洲| 精品久久久精品久久久| 国产精品国产三级国产专区5o| 两个人免费观看高清视频 | 欧美日韩综合久久久久久| 日本色播在线视频| 成年人午夜在线观看视频| 能在线免费看毛片的网站| 欧美精品亚洲一区二区| 美女cb高潮喷水在线观看| 国产精品99久久99久久久不卡 | 欧美人与善性xxx| 日韩不卡一区二区三区视频在线| 国产精品欧美亚洲77777| 亚洲av日韩在线播放| 蜜桃久久精品国产亚洲av| av免费观看日本| 男的添女的下面高潮视频| 国产视频首页在线观看| 蜜桃久久精品国产亚洲av| 全区人妻精品视频| 国产精品女同一区二区软件| 久久韩国三级中文字幕| 日韩 亚洲 欧美在线| 精品一区二区三区视频在线| 国产在视频线精品| 久久 成人 亚洲| 亚洲婷婷狠狠爱综合网| 欧美 亚洲 国产 日韩一| 99久国产av精品国产电影| 少妇熟女欧美另类| 亚洲美女搞黄在线观看| 天天躁夜夜躁狠狠久久av| 人人澡人人妻人| 草草在线视频免费看| 看非洲黑人一级黄片| 男男h啪啪无遮挡| 久久 成人 亚洲| 亚洲成人手机| 精品熟女少妇av免费看| 观看av在线不卡| 国产av国产精品国产| 妹子高潮喷水视频| 精品少妇久久久久久888优播| 嫩草影院新地址| 亚洲国产精品一区三区| 午夜免费鲁丝| 高清黄色对白视频在线免费看 | 日韩制服骚丝袜av| 婷婷色综合www| 国产精品伦人一区二区| 少妇高潮的动态图| 国产成人精品婷婷| 亚洲av在线观看美女高潮| 免费看光身美女| 搡女人真爽免费视频火全软件| 最近最新中文字幕免费大全7| 日韩,欧美,国产一区二区三区| 97精品久久久久久久久久精品| 久久久久久久大尺度免费视频| 免费少妇av软件| 久久人妻熟女aⅴ| 精品久久久久久久久av| 中文字幕制服av| 亚洲欧美精品自产自拍| 九色成人免费人妻av| 日韩欧美精品免费久久| 丰满饥渴人妻一区二区三| 曰老女人黄片| 91久久精品国产一区二区成人| 99久国产av精品国产电影| 人妻制服诱惑在线中文字幕| 人人澡人人妻人| 黄色欧美视频在线观看| 日韩制服骚丝袜av| 在线 av 中文字幕| 久久久a久久爽久久v久久| 一级毛片黄色毛片免费观看视频| 欧美bdsm另类| 亚洲欧美清纯卡通| 80岁老熟妇乱子伦牲交| 最新的欧美精品一区二区| 六月丁香七月| 亚洲av二区三区四区| 久久 成人 亚洲| 国产毛片在线视频| 精品一区在线观看国产| 日韩精品免费视频一区二区三区 | 国产美女午夜福利| 欧美变态另类bdsm刘玥| 久久热精品热| 成人18禁高潮啪啪吃奶动态图 | 91精品国产国语对白视频| 99久久精品一区二区三区| 国产亚洲欧美精品永久| 日本欧美视频一区| 日韩制服骚丝袜av| 国产乱人偷精品视频| 99久久人妻综合| 少妇被粗大猛烈的视频| 亚洲欧美一区二区三区国产| 亚洲久久久国产精品| 婷婷色麻豆天堂久久| 久久久精品94久久精品| 国产亚洲一区二区精品| 国产精品99久久久久久久久| 日韩强制内射视频| 精品少妇黑人巨大在线播放| 欧美另类一区| 国产精品一区二区在线不卡| 亚洲人成网站在线播| 国产永久视频网站| 青青草视频在线视频观看| 久久精品国产鲁丝片午夜精品| 精品熟女少妇av免费看| 成人影院久久| 天堂俺去俺来也www色官网| 国产精品人妻久久久久久| 欧美国产精品一级二级三级 | 日韩制服骚丝袜av| 国产高清有码在线观看视频| 人人妻人人澡人人爽人人夜夜| 久久6这里有精品| 性高湖久久久久久久久免费观看| 26uuu在线亚洲综合色| 一本大道久久a久久精品| 中文字幕免费在线视频6| 精品久久久噜噜| 男的添女的下面高潮视频| 国精品久久久久久国模美| 亚洲欧美成人精品一区二区| 欧美另类一区| 久热久热在线精品观看| 免费人妻精品一区二区三区视频| 久久国产亚洲av麻豆专区| 乱码一卡2卡4卡精品| 午夜av观看不卡| 午夜福利影视在线免费观看| 久久久久久人妻| 亚洲国产精品999| 91久久精品电影网| 女人精品久久久久毛片| 精品久久久久久电影网| 三级经典国产精品| 我的女老师完整版在线观看| 中文字幕精品免费在线观看视频 | 最近中文字幕高清免费大全6| 久久这里有精品视频免费| 观看av在线不卡| av卡一久久| 国产伦精品一区二区三区四那| 日产精品乱码卡一卡2卡三| 伦理电影免费视频| 晚上一个人看的免费电影| 久久久国产精品麻豆| 久久久久视频综合| 香蕉精品网在线| 黄色日韩在线| 亚洲内射少妇av| 国产在线一区二区三区精| 十八禁网站网址无遮挡 | av不卡在线播放| 国产精品蜜桃在线观看| 一区二区三区精品91| 亚洲精品国产成人久久av| 久久人人爽人人片av| 久久久久久久久久久久大奶| av卡一久久| 插逼视频在线观看| 午夜老司机福利剧场| 男人添女人高潮全过程视频| 亚洲综合色惰| 老女人水多毛片| 美女国产视频在线观看| 免费黄色在线免费观看| 亚洲av中文av极速乱| 嘟嘟电影网在线观看| 欧美日韩综合久久久久久| 在线 av 中文字幕| 国产高清不卡午夜福利| 亚洲高清免费不卡视频| av视频免费观看在线观看| 男女边吃奶边做爰视频| 精品视频人人做人人爽| 91久久精品国产一区二区三区| 两个人免费观看高清视频 | 两个人免费观看高清视频 | 久久狼人影院| 国产精品秋霞免费鲁丝片| 免费看av在线观看网站| 少妇人妻久久综合中文| 老司机亚洲免费影院| 看十八女毛片水多多多| 久久精品久久精品一区二区三区| 中文资源天堂在线| 国产av国产精品国产| 人人妻人人看人人澡| 久久久国产一区二区| 欧美 亚洲 国产 日韩一| 亚洲精品第二区| 在线看a的网站| 亚洲成人手机| 国产色婷婷99| 国产一区二区三区av在线| 另类精品久久| 大香蕉久久网| av在线播放精品| 黄色一级大片看看| 免费观看性生交大片5| 国产成人一区二区在线| 视频中文字幕在线观看| h日本视频在线播放| 欧美精品高潮呻吟av久久| 国产毛片在线视频| 日韩强制内射视频| 午夜免费观看性视频| 天天操日日干夜夜撸| videossex国产| 亚洲精华国产精华液的使用体验| 18禁裸乳无遮挡动漫免费视频| 男女国产视频网站| 97在线视频观看| 国产免费一级a男人的天堂| 成人国产麻豆网| 亚洲欧美一区二区三区国产| 九草在线视频观看| 99久久中文字幕三级久久日本| 三级国产精品片| 天天躁夜夜躁狠狠久久av| 国产精品久久久久成人av| 国产精品99久久久久久久久| 18禁在线播放成人免费| 亚洲国产精品999| 色哟哟·www| 亚洲av成人精品一区久久| 久久99精品国语久久久| 久久99热这里只频精品6学生| 熟妇人妻不卡中文字幕| 精品少妇黑人巨大在线播放| 国产女主播在线喷水免费视频网站| 在线观看三级黄色| 午夜福利影视在线免费观看| 久久久a久久爽久久v久久| 成人毛片60女人毛片免费| 欧美少妇被猛烈插入视频| 免费观看的影片在线观看| 一个人免费看片子| 国产熟女欧美一区二区| av国产久精品久网站免费入址| 99精国产麻豆久久婷婷| 成年女人在线观看亚洲视频| 日日摸夜夜添夜夜添av毛片| 亚洲欧美精品专区久久| 人人妻人人添人人爽欧美一区卜| 三级国产精品片| 老司机影院毛片| 一区二区av电影网| 嘟嘟电影网在线观看| 一级毛片 在线播放| 在线亚洲精品国产二区图片欧美 | 日韩av不卡免费在线播放| 在线 av 中文字幕| 国产精品国产三级专区第一集| 国产乱来视频区| 制服丝袜香蕉在线| 亚州av有码| 51国产日韩欧美| 少妇高潮的动态图| 少妇的逼好多水| 精品亚洲乱码少妇综合久久| 只有这里有精品99| 乱系列少妇在线播放| 久久久午夜欧美精品| 嘟嘟电影网在线观看| 最新中文字幕久久久久| 乱人伦中国视频| 国产精品国产三级国产av玫瑰| 午夜影院在线不卡| 最新的欧美精品一区二区| 成年av动漫网址| 女人精品久久久久毛片| 十八禁高潮呻吟视频 | 亚洲欧洲日产国产| 热re99久久精品国产66热6| 国产精品人妻久久久久久| 久久久久视频综合| 一个人免费看片子| 十分钟在线观看高清视频www | 九九在线视频观看精品| 亚洲综合精品二区| 欧美日韩综合久久久久久| 免费在线观看成人毛片| 一区在线观看完整版| 少妇熟女欧美另类| 免费大片黄手机在线观看| 日韩中字成人| 久久 成人 亚洲| 欧美3d第一页| 欧美xxⅹ黑人| 国产精品免费大片| 青春草视频在线免费观看| 多毛熟女@视频| 青春草视频在线免费观看| 亚洲精品国产色婷婷电影| 水蜜桃什么品种好| 亚洲av欧美aⅴ国产| 亚洲精品日本国产第一区| 人人妻人人澡人人爽人人夜夜| 国产成人a∨麻豆精品| 久久av网站| 3wmmmm亚洲av在线观看| 最后的刺客免费高清国语| 久热这里只有精品99| 哪个播放器可以免费观看大片| 亚洲国产色片| 欧美精品人与动牲交sv欧美| 在线 av 中文字幕| 交换朋友夫妻互换小说| 少妇 在线观看| 欧美日韩综合久久久久久| 国产片特级美女逼逼视频| 午夜日本视频在线| 在线亚洲精品国产二区图片欧美 | 人体艺术视频欧美日本| 一级毛片久久久久久久久女| 亚洲图色成人| 街头女战士在线观看网站| 亚洲色图综合在线观看| 国产色婷婷99| 黑人高潮一二区| 亚洲人成网站在线观看播放| 大片免费播放器 马上看| 久久国产精品大桥未久av | 久久久国产精品麻豆| 好男人视频免费观看在线| 精品国产露脸久久av麻豆| 国产 精品1| 九九久久精品国产亚洲av麻豆| 少妇的逼水好多| 亚洲综合精品二区| 亚洲美女黄色视频免费看| 啦啦啦啦在线视频资源| 最黄视频免费看| 超碰97精品在线观看| 少妇被粗大猛烈的视频| 久久精品国产亚洲av天美| 女人精品久久久久毛片| 国产伦精品一区二区三区四那| 亚洲综合精品二区|