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

    LNG船靠離泊操縱運(yùn)動(dòng)數(shù)值模擬分析

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

    ,,,,

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

    LNG船靠離泊操縱運(yùn)動(dòng)數(shù)值模擬分析

    鄭坤1,時(shí)光志1,張雨新2,周毅1,陸晟2

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

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

    靠離泊;操縱運(yùn)動(dòng);數(shù)值模擬;CFD

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

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

    1 工作水域概況

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

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

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

    2 船型及計(jì)算工況

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

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

    表1 船型主尺度 m

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

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

    表2 計(jì)算工況

    3 數(shù)值模擬方法

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

    3.1 控制方程

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

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

    (1)

    (2)

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

    SSTk-ω模型湍動(dòng)能和湍動(dòng)耗散率的求解方程為

    (3)

    (4)

    VOF控制方程為

    (5)

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

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

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

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

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

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

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

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

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

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

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

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

    船體的運(yùn)動(dòng)軌跡如圖5所示。

    圖5 靠泊時(shí)船體運(yùn)動(dòng)軌跡(時(shí)間間隔為30 s,虛線為示意軌跡)

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

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

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

    圖6 離泊時(shí)船體運(yùn)動(dòng)軌跡(時(shí)間間隔為15 s,虛線為示意軌跡)

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

    5 結(jié)論

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

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

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

    [3] 史宏達(dá),范天會(huì),韓桂萍.東營(yíng)港擴(kuò)建工程碼頭船舶靠離泊數(shù)值模擬研究[J].中國(guó)海洋大學(xué)學(xué)報(bào),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] 時(shí)光志,周志峰.某液化天然氣運(yùn)輸船推進(jìn)器選型研究[J].船海工程,2014,43(5):132-137.

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

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

    [10] 吳召華,陳作鋼,代燚.基于體積力法的船體自航性能數(shù)值預(yù)報(bào)[J].上海交通大學(xué)學(xué)報(bào),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è)計(jì)制造

    猜你喜歡
    船舶
    船舶避碰路徑模糊控制系統(tǒng)
    計(jì)算流體力學(xué)在船舶操縱運(yùn)動(dòng)仿真中的應(yīng)用
    CM節(jié)點(diǎn)控制在船舶上的應(yīng)用
    基于改進(jìn)譜分析法的船舶疲勞強(qiáng)度直接計(jì)算
    《船舶》2022 年度征訂啟事
    船舶(2021年4期)2021-09-07 17:32:22
    船舶!請(qǐng)加速
    BOG壓縮機(jī)在小型LNG船舶上的應(yīng)用
    船舶 揚(yáng)帆奮起
    軍工文化(2017年12期)2017-07-17 06:08:06
    船舶壓載水管理系統(tǒng)
    小型船舶艉軸架設(shè)計(jì)
    船海工程(2015年4期)2016-01-05 15:53:30
    久久久成人免费电影| 听说在线观看完整版免费高清| 亚洲美女视频黄频| 亚洲精品色激情综合| 国产精品一及| 国产私拍福利视频在线观看| 97人妻精品一区二区三区麻豆| 深夜精品福利| 久久人人爽人人爽人人片va| 美女黄网站色视频| 熟妇人妻久久中文字幕3abv| 亚洲专区国产一区二区| 亚洲人成网站在线播放欧美日韩| 欧美三级亚洲精品| 一进一出抽搐gif免费好疼| 男人舔奶头视频| 日本色播在线视频| 色综合亚洲欧美另类图片| 日本精品一区二区三区蜜桃| 狂野欧美激情性xxxx在线观看| 国产成人a∨麻豆精品| 久久午夜亚洲精品久久| 能在线免费观看的黄片| 成年女人毛片免费观看观看9| 老师上课跳d突然被开到最大视频| 成人精品一区二区免费| 欧美人与善性xxx| 又爽又黄无遮挡网站| 乱人视频在线观看| avwww免费| 午夜福利在线在线| 午夜影院日韩av| 亚洲成人中文字幕在线播放| 久久精品国产亚洲av香蕉五月| 色吧在线观看| 黑人高潮一二区| 一级毛片电影观看 | 精品欧美国产一区二区三| 51国产日韩欧美| 免费看a级黄色片| 嫩草影视91久久| 久久欧美精品欧美久久欧美| 少妇被粗大猛烈的视频| 亚洲综合色惰| 久久99热这里只有精品18| 国产三级中文精品| 国产精品一区二区性色av| 国产一级毛片七仙女欲春2| 国产伦在线观看视频一区| 国产在视频线在精品| 成人综合一区亚洲| 搡老熟女国产l中国老女人| 免费av不卡在线播放| ponron亚洲| 五月伊人婷婷丁香| 国产一区亚洲一区在线观看| 免费大片18禁| av视频在线观看入口| 欧美在线一区亚洲| 国产aⅴ精品一区二区三区波| 变态另类成人亚洲欧美熟女| 久久中文看片网| 男人的好看免费观看在线视频| 中文在线观看免费www的网站| 在线a可以看的网站| 久久精品国产亚洲av香蕉五月| 国产成人a∨麻豆精品| 亚洲不卡免费看| 国产黄片美女视频| 深夜精品福利| 1024手机看黄色片| 尾随美女入室| 亚洲美女视频黄频| 好男人在线观看高清免费视频| 亚洲在线自拍视频| 日韩欧美 国产精品| 99在线人妻在线中文字幕| 亚洲欧美日韩高清专用| 久久人妻av系列| 久久亚洲国产成人精品v| 亚洲国产日韩欧美精品在线观看| 国产亚洲精品久久久久久毛片| 免费看a级黄色片| 国产成年人精品一区二区| 亚洲七黄色美女视频| 国产欧美日韩一区二区精品| 久久久国产成人精品二区| 国产成人a区在线观看| 又黄又爽又免费观看的视频| 身体一侧抽搐| 欧美日韩在线观看h| 我的女老师完整版在线观看| 嫩草影院精品99| 最后的刺客免费高清国语| 国产单亲对白刺激| av在线老鸭窝| 欧美日韩国产亚洲二区| 成年女人看的毛片在线观看| 国产精品无大码| 国产精品一区二区三区四区久久| 男女之事视频高清在线观看| 最近在线观看免费完整版| 一卡2卡三卡四卡精品乱码亚洲| 少妇的逼水好多| 插逼视频在线观看| 日韩欧美精品v在线| 丰满人妻一区二区三区视频av| 免费观看人在逋| 老司机影院成人| 中文字幕免费在线视频6| 亚洲人成网站在线观看播放| 国产视频一区二区在线看| 亚洲精品日韩在线中文字幕 | 国产精品福利在线免费观看| 国产精品人妻久久久影院| 最近中文字幕高清免费大全6| 国产淫片久久久久久久久| 精品久久久久久成人av| 日日摸夜夜添夜夜添小说| 久久鲁丝午夜福利片| 我的老师免费观看完整版| www日本黄色视频网| 色综合站精品国产| 亚洲精品国产av成人精品 | 自拍偷自拍亚洲精品老妇| 深夜精品福利| 国产成人a∨麻豆精品| 韩国av在线不卡| 亚洲av成人精品一区久久| 国产毛片a区久久久久| av在线观看视频网站免费| 国产精品三级大全| 亚洲精品国产av成人精品 | 久久精品91蜜桃| 嫩草影视91久久| 自拍偷自拍亚洲精品老妇| 一个人看的www免费观看视频| 22中文网久久字幕| 精品一区二区免费观看| 国产一区二区三区在线臀色熟女| 亚洲五月天丁香| 在线看三级毛片| 久久久精品大字幕| 欧美+日韩+精品| 一进一出抽搐gif免费好疼| 精品人妻偷拍中文字幕| 成人国产麻豆网| 精品午夜福利在线看| 色综合站精品国产| 日本a在线网址| 99国产极品粉嫩在线观看| 自拍偷自拍亚洲精品老妇| 色综合色国产| 午夜福利高清视频| 国产高清激情床上av| 国产成人a∨麻豆精品| 亚洲欧美精品综合久久99| 中文字幕av在线有码专区| 悠悠久久av| 亚洲国产欧美人成| 国产麻豆成人av免费视频| 最近最新中文字幕大全电影3| 久久人妻av系列| 又黄又爽又刺激的免费视频.| 国产精品国产高清国产av| 91午夜精品亚洲一区二区三区| 亚洲一级一片aⅴ在线观看| 国产高清有码在线观看视频| 好男人在线观看高清免费视频| 欧洲精品卡2卡3卡4卡5卡区| 国产精品99久久久久久久久| 日产精品乱码卡一卡2卡三| 国产国拍精品亚洲av在线观看| 国产日本99.免费观看| 国产毛片a区久久久久| 国产一区亚洲一区在线观看| 亚洲av电影不卡..在线观看| 亚洲激情五月婷婷啪啪| 亚洲国产精品成人综合色| 国产精品一区二区免费欧美| 色噜噜av男人的天堂激情| 日本一二三区视频观看| 国产成人一区二区在线| 天堂动漫精品| 欧美高清性xxxxhd video| 1000部很黄的大片| 亚洲欧美日韩卡通动漫| 18+在线观看网站| 国产高清激情床上av| 可以在线观看的亚洲视频| 亚洲欧美日韩东京热| 国产精品国产三级国产av玫瑰| 三级毛片av免费| 老熟妇乱子伦视频在线观看| 亚洲性夜色夜夜综合| 亚洲不卡免费看| 22中文网久久字幕| 俺也久久电影网| 日韩欧美在线乱码| 丰满的人妻完整版| 久久草成人影院| 亚洲一区二区三区色噜噜| 综合色丁香网| 国产伦精品一区二区三区四那| 精品久久久久久久久亚洲| 99久国产av精品国产电影| 亚洲,欧美,日韩| 在线免费观看的www视频| 日本撒尿小便嘘嘘汇集6| 欧美最新免费一区二区三区| 日韩欧美免费精品| 成人综合一区亚洲| 成人国产麻豆网| 国产黄色视频一区二区在线观看 | 亚洲精品色激情综合| 亚洲精品久久国产高清桃花| 特级一级黄色大片| 免费不卡的大黄色大毛片视频在线观看 | 国产精品电影一区二区三区| 香蕉av资源在线| 国产日本99.免费观看| 99视频精品全部免费 在线| 中文字幕熟女人妻在线| 欧美一级a爱片免费观看看| 97人妻精品一区二区三区麻豆| 久久九九热精品免费| 免费看光身美女| 亚洲内射少妇av| 欧美极品一区二区三区四区| 最好的美女福利视频网| 精品人妻熟女av久视频| 永久网站在线| 久久久久国内视频| 成年版毛片免费区| 国产一区二区三区在线臀色熟女| 联通29元200g的流量卡| 99久久中文字幕三级久久日本| 美女高潮的动态| 亚洲一区二区三区色噜噜| 99国产精品一区二区蜜桃av| 国产午夜精品论理片| 国产一区亚洲一区在线观看| 精品久久久久久久久久久久久| av在线观看视频网站免费| 久久国内精品自在自线图片| 三级经典国产精品| 国产精品精品国产色婷婷| 久久久a久久爽久久v久久| 深夜a级毛片| 午夜爱爱视频在线播放| 不卡一级毛片| 午夜免费激情av| 日韩,欧美,国产一区二区三区 | 丰满人妻一区二区三区视频av| 欧美中文日本在线观看视频| 国产精品久久久久久亚洲av鲁大| 亚洲国产精品成人久久小说 | 免费av毛片视频| 最新在线观看一区二区三区| 日日啪夜夜撸| 国产精品美女特级片免费视频播放器| 午夜福利在线观看免费完整高清在 | 亚洲高清免费不卡视频| 亚洲欧美日韩东京热| 97在线视频观看| 国产大屁股一区二区在线视频| 午夜福利18| 亚洲最大成人手机在线| 三级男女做爰猛烈吃奶摸视频| 免费大片18禁| 久久亚洲国产成人精品v| 日本精品一区二区三区蜜桃| 内地一区二区视频在线| 亚洲中文字幕一区二区三区有码在线看| 国产男人的电影天堂91| 欧美成人a在线观看| 亚洲电影在线观看av| 中文字幕熟女人妻在线| 九九爱精品视频在线观看| 亚洲中文字幕日韩| 亚洲无线观看免费| 一a级毛片在线观看| 午夜免费激情av| 亚洲精华国产精华液的使用体验 | 国产精品亚洲一级av第二区| 国产精品一区二区三区四区久久| 国产美女午夜福利| 日日摸夜夜添夜夜添av毛片| 少妇被粗大猛烈的视频| 国产亚洲精品久久久com| 欧美潮喷喷水| 少妇人妻一区二区三区视频| 日韩av在线大香蕉| 一本精品99久久精品77| 99久久久亚洲精品蜜臀av| 深爱激情五月婷婷| 欧洲精品卡2卡3卡4卡5卡区| 国产午夜福利久久久久久| 欧美性猛交黑人性爽| 人妻久久中文字幕网| 精品熟女少妇av免费看| 国产精品精品国产色婷婷| 久久人人爽人人爽人人片va| 亚洲第一电影网av| 久久99热这里只有精品18| 特大巨黑吊av在线直播| 欧美日韩在线观看h| 中文字幕精品亚洲无线码一区| av.在线天堂| 亚洲婷婷狠狠爱综合网| av天堂中文字幕网| 日韩欧美精品免费久久| 99视频精品全部免费 在线| 久久韩国三级中文字幕| 欧美不卡视频在线免费观看| 欧美zozozo另类| 特级一级黄色大片| 校园人妻丝袜中文字幕| 久久久久久久久久成人| 国产黄色小视频在线观看| 久久久欧美国产精品| 亚洲精品粉嫩美女一区| 欧美最黄视频在线播放免费| 国产av不卡久久| 国产在线精品亚洲第一网站| 99久久精品热视频| 亚洲第一电影网av| 香蕉av资源在线| 色尼玛亚洲综合影院| 欧美bdsm另类| 亚洲不卡免费看| 少妇被粗大猛烈的视频| a级毛色黄片| 网址你懂的国产日韩在线| 国产黄a三级三级三级人| 国产日本99.免费观看| 欧美最新免费一区二区三区| 九色成人免费人妻av| 男人狂女人下面高潮的视频| 在线a可以看的网站| 少妇的逼好多水| 国产v大片淫在线免费观看| 国产成人精品久久久久久| 欧美性感艳星| 国产大屁股一区二区在线视频| 午夜福利在线在线| 一边摸一边抽搐一进一小说| 午夜免费男女啪啪视频观看 | 免费看日本二区| 九九在线视频观看精品| a级毛片a级免费在线| 亚洲久久久久久中文字幕| 亚洲一区二区三区色噜噜| а√天堂www在线а√下载| 国产在线精品亚洲第一网站| 亚洲成人精品中文字幕电影| 春色校园在线视频观看| 尾随美女入室| 免费不卡的大黄色大毛片视频在线观看 | 国产淫片久久久久久久久| 色综合色国产| 日韩强制内射视频| 精品人妻一区二区三区麻豆 | 99九九线精品视频在线观看视频| 久久久欧美国产精品| 一级毛片电影观看 | 午夜激情福利司机影院| 亚洲精品一区av在线观看| 韩国av在线不卡| 日产精品乱码卡一卡2卡三| 久久人人爽人人片av| 亚洲精品日韩在线中文字幕 | 日日摸夜夜添夜夜添av毛片| 少妇裸体淫交视频免费看高清| 国产午夜精品论理片| 久久综合国产亚洲精品| 久久久久国产网址| 久久天躁狠狠躁夜夜2o2o| .国产精品久久| 久久午夜福利片| 一级黄色大片毛片| 毛片女人毛片| 99久久久亚洲精品蜜臀av| 国内精品宾馆在线| 大香蕉久久网| 最近视频中文字幕2019在线8| 亚洲美女搞黄在线观看 | 1024手机看黄色片| 成年版毛片免费区| 美女xxoo啪啪120秒动态图| 国内久久婷婷六月综合欲色啪| 亚洲精品在线观看二区| 99久久精品国产国产毛片| 国产精品乱码一区二三区的特点| 男女边吃奶边做爰视频| 日本在线视频免费播放| 日本与韩国留学比较| 好男人在线观看高清免费视频| h日本视频在线播放| 听说在线观看完整版免费高清| 99热这里只有是精品50| 3wmmmm亚洲av在线观看| 狂野欧美白嫩少妇大欣赏| 女同久久另类99精品国产91| 综合色丁香网| 日韩一区二区视频免费看| 久久精品国产清高在天天线| 成人二区视频| 天天躁夜夜躁狠狠久久av| 午夜老司机福利剧场| 久久久久国产网址| 听说在线观看完整版免费高清| 久久久精品大字幕| 久久这里只有精品中国| 一级黄片播放器| 精品熟女少妇av免费看| 日韩av在线大香蕉| 18+在线观看网站| 免费高清视频大片| 亚洲不卡免费看| 精华霜和精华液先用哪个| 夜夜夜夜夜久久久久| 日韩av不卡免费在线播放| 亚洲最大成人手机在线| 三级男女做爰猛烈吃奶摸视频| 在线观看美女被高潮喷水网站| 国内精品久久久久精免费| 在线观看一区二区三区| 成年免费大片在线观看| 高清午夜精品一区二区三区 | 日韩高清综合在线| 日本熟妇午夜| 黄色视频,在线免费观看| 国产精品电影一区二区三区| 国产精品日韩av在线免费观看| 中文资源天堂在线| 卡戴珊不雅视频在线播放| 中文亚洲av片在线观看爽| 老司机影院成人| 男人舔奶头视频| 干丝袜人妻中文字幕| 男女下面进入的视频免费午夜| 欧美高清性xxxxhd video| 亚洲一级一片aⅴ在线观看| 国产精品电影一区二区三区| 在线观看美女被高潮喷水网站| 国产精品一区二区三区四区免费观看 | 国产欧美日韩精品一区二区| 俄罗斯特黄特色一大片| 亚洲av二区三区四区| 给我免费播放毛片高清在线观看| 大型黄色视频在线免费观看| 免费看日本二区| 久久婷婷人人爽人人干人人爱| 观看免费一级毛片| 久久久精品94久久精品| 晚上一个人看的免费电影| 国产一级毛片七仙女欲春2| 久久天躁狠狠躁夜夜2o2o| 国内精品一区二区在线观看| 嫩草影院精品99| 麻豆一二三区av精品| 国产白丝娇喘喷水9色精品| 一级毛片久久久久久久久女| 免费一级毛片在线播放高清视频| 国产男人的电影天堂91| 久久人妻av系列| www.色视频.com| 成人鲁丝片一二三区免费| 日日啪夜夜撸| 狂野欧美白嫩少妇大欣赏| 1024手机看黄色片| 日韩国内少妇激情av| 一区二区三区四区激情视频 | 日韩在线高清观看一区二区三区| 亚洲内射少妇av| 99热精品在线国产| 国内揄拍国产精品人妻在线| 免费黄网站久久成人精品| 久久精品夜夜夜夜夜久久蜜豆| 久久99热6这里只有精品| 色在线成人网| 久久99热6这里只有精品| 最好的美女福利视频网| 亚洲精品456在线播放app| avwww免费| 亚洲av电影不卡..在线观看| 亚洲成人中文字幕在线播放| 91精品国产九色| 一级a爱片免费观看的视频| 免费看av在线观看网站| 村上凉子中文字幕在线| 欧美色欧美亚洲另类二区| 日日啪夜夜撸| aaaaa片日本免费| 亚洲精品亚洲一区二区| 成人av在线播放网站| 成人毛片a级毛片在线播放| 乱系列少妇在线播放| 欧美一区二区精品小视频在线| 插阴视频在线观看视频| 亚洲经典国产精华液单| 欧美在线一区亚洲| 黄色一级大片看看| 成人午夜高清在线视频| 日韩,欧美,国产一区二区三区 | 熟女电影av网| 高清毛片免费观看视频网站| 在线a可以看的网站| 精品久久久久久久末码| 欧美最新免费一区二区三区| 少妇猛男粗大的猛烈进出视频 | 色av中文字幕| 国内揄拍国产精品人妻在线| 免费黄网站久久成人精品| 精品国产三级普通话版| 亚洲不卡免费看| 色综合亚洲欧美另类图片| 国产乱人偷精品视频| 精品久久久久久久久av| 99久久九九国产精品国产免费| 国产欧美日韩一区二区精品| www日本黄色视频网| 成人综合一区亚洲| а√天堂www在线а√下载| 不卡一级毛片| 精品午夜福利视频在线观看一区| 欧美精品国产亚洲| 搞女人的毛片| 久久午夜亚洲精品久久| 日本色播在线视频| 国产探花在线观看一区二区| 日日摸夜夜添夜夜爱| 精品福利观看| 俺也久久电影网| 国产一级毛片七仙女欲春2| 亚洲欧美日韩高清专用| 悠悠久久av| 乱系列少妇在线播放| 精品人妻偷拍中文字幕| 国产精品美女特级片免费视频播放器| 午夜福利高清视频| 中文字幕精品亚洲无线码一区| 久久欧美精品欧美久久欧美| 国产视频内射| 国产 一区精品| 蜜桃亚洲精品一区二区三区| 久久亚洲精品不卡| 国产 一区 欧美 日韩| 亚洲av成人av| 午夜日韩欧美国产| 免费大片18禁| 伦理电影大哥的女人| 六月丁香七月| 在线a可以看的网站| 99久国产av精品国产电影| 大型黄色视频在线免费观看| 女人被狂操c到高潮| 精品久久久久久久久亚洲| 久久久久久久久久久丰满| 国产高清视频在线观看网站| 成人精品一区二区免费| 欧美一区二区亚洲| 一个人看视频在线观看www免费| 性插视频无遮挡在线免费观看| 一级黄色大片毛片| 免费观看的影片在线观看| 国产女主播在线喷水免费视频网站 | 我要搜黄色片| 亚洲欧美日韩东京热| 99久久无色码亚洲精品果冻| 女同久久另类99精品国产91| 国内精品宾馆在线| 天天躁夜夜躁狠狠久久av| 日本撒尿小便嘘嘘汇集6| 一级a爱片免费观看的视频| 中文字幕免费在线视频6| 深夜精品福利| 久久久久国产网址| 一进一出抽搐动态| 久久这里只有精品中国| 亚洲精华国产精华液的使用体验 | 日韩欧美国产在线观看| 一级毛片久久久久久久久女| 日本黄大片高清| 精品99又大又爽又粗少妇毛片| 一级毛片电影观看 | 十八禁国产超污无遮挡网站| 我要搜黄色片| 一区二区三区高清视频在线| 亚洲精品乱码久久久v下载方式| 国产精品爽爽va在线观看网站| 久久精品91蜜桃| 精品人妻一区二区三区麻豆 | 免费人成在线观看视频色| 嫩草影院入口| 亚洲自拍偷在线| 亚洲丝袜综合中文字幕| 97热精品久久久久久| 村上凉子中文字幕在线| 免费人成在线观看视频色| 亚洲av第一区精品v没综合| 成人亚洲欧美一区二区av| 日本一二三区视频观看| 人妻少妇偷人精品九色| 欧美+亚洲+日韩+国产| 又爽又黄无遮挡网站| 久久久久久久久大av| a级一级毛片免费在线观看| 我要看日韩黄色一级片| 亚洲国产日韩欧美精品在线观看| 国产精品久久久久久av不卡| 久久精品91蜜桃| 非洲黑人性xxxx精品又粗又长| 国语自产精品视频在线第100页|