999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

島礁地形下半潛平臺水動力性能研究

2015-04-25 10:11:48程小明劉小龍中國船舶科學(xué)研究中心江蘇無錫214082
船舶力學(xué) 2015年7期
關(guān)鍵詞:方向分析

吳 波,程小明,田 超,丁 軍,劉小龍(中國船舶科學(xué)研究中心,江蘇 無錫 214082)

島礁地形下半潛平臺水動力性能研究

吳 波,程小明,田 超,丁 軍,劉小龍
(中國船舶科學(xué)研究中心,江蘇 無錫 214082)

文章基于一島礁實(shí)地的地形數(shù)據(jù),針對該地形下的半潛式平臺進(jìn)行了水動力及運(yùn)動響應(yīng)分析,重點(diǎn)研究了島礁地形對平臺運(yùn)動響應(yīng)的影響,并結(jié)合模型試驗(yàn)結(jié)果進(jìn)行了對比分析。結(jié)果表明應(yīng)用目前常用的以線性波浪理論為基礎(chǔ)的數(shù)值方法,通過引入對地形與平臺之間的繞射與輻射水動力的干擾計(jì)算,可以對島礁地形下的半潛平臺進(jìn)行水動力和運(yùn)動分析。研究所得到的島礁條件下的平臺運(yùn)動響應(yīng)結(jié)果與模型試驗(yàn)值對比較為吻合。

島礁地形;半潛平臺;水動力性能;數(shù)值分析

0 引 言

海洋中蘊(yùn)含有豐富的油氣資源和漁業(yè)資源,開發(fā)和利用海洋包括島礁的資源已經(jīng)成為經(jīng)濟(jì)與技術(shù)發(fā)展的一個(gè)重要趨勢。島礁周邊環(huán)繞著寬度數(shù)百米至數(shù)千米水深較淺的礁盤,這些礁盤不僅具有豐富的漁業(yè)和旅游資源,而且發(fā)揮著消波護(hù)島、防止海水侵蝕島體本身的重要作用。應(yīng)用浮式平臺作為近島礁的浮動基地,在不影響島礁生態(tài)環(huán)境、礁盤和海島陸地形態(tài)的同時(shí),可以為海洋及島礁的建設(shè)與資源開發(fā)提供后勤基地。島礁半潛式平臺布放于近島礁海域,通過棧橋與島礁相連,作為島礁的生產(chǎn)和建設(shè)基地,平臺具有人員作業(yè)生活、船舶停靠、物資裝卸等功能;采用浮式平臺結(jié)構(gòu),具有部署快速、使用靈活的特點(diǎn)。圖1所示為島礁與島礁半潛式平臺的示意圖。

圖1 近島礁半潛平臺Fig.1 Semi-submersible platform near a reef island

對于海洋平臺水動力性能的分析,國內(nèi)外在該領(lǐng)域開展了大量的工作。王世圣[1]針對深水半潛平臺,建立半潛式鉆井平臺的三維水動力模型,采用勢流理論和Morison公式相結(jié)合的方法,應(yīng)用SEASAM進(jìn)行計(jì)算分析,得到半潛平臺運(yùn)動響應(yīng)的傳遞函數(shù)和固有周期。朱航[2]針對風(fēng)浪作用下的HYSY-981半潛式平臺,開展時(shí)域的數(shù)值模擬分析,得到平臺的運(yùn)動響應(yīng)結(jié)果。張威[3]對1 500 m的深海半潛式平臺進(jìn)行研究,通過頻域計(jì)算得到平臺的運(yùn)動響應(yīng),并基于模型試驗(yàn)結(jié)果進(jìn)行對比分析。肖宇維[4]應(yīng)用AQWA研究了三種型式的張力腿平臺在規(guī)則波下六自由度運(yùn)動的響應(yīng)幅值算子,并對比分析了水動力性能的特點(diǎn)。姚彥龍[5]針對三角形張力腿平臺開展頻域水動力分析,得到該平臺的運(yùn)動響應(yīng)、波浪激勵(lì)力,并在此基礎(chǔ)上對平臺的性能進(jìn)行評估,對布局進(jìn)行優(yōu)化。蒙占彬[6]基于1 500 m水深的環(huán)境條件,對傳統(tǒng)式、延伸式和海星式三種型式的張力腿平臺的運(yùn)動響應(yīng)及波浪長期預(yù)報(bào)響應(yīng)幅值進(jìn)行了對比,分析各自的水動力性能及優(yōu)缺點(diǎn)。龍?bào)銜蟍7]采用時(shí)域耦合分析方法,研究了Truss Spar平臺的耦合動力特性問題,分析平臺的運(yùn)動響應(yīng)。Mansour[8]對新型半潛平臺開展數(shù)值計(jì)算分析,對比常規(guī)型的半潛平臺的RAO值進(jìn)行研究。Clauss[9]針對GVA4000半潛平臺進(jìn)行數(shù)值計(jì)算,所得到的波浪誘導(dǎo)運(yùn)動數(shù)值結(jié)果與試驗(yàn)結(jié)果對比吻合良好。Lowa[10]在時(shí)域和頻域中對浮式平臺進(jìn)行分析,研究平臺的一階和二階運(yùn)動。國內(nèi)外的研究對象大多針對深水條件下的海洋平臺,在淺水及近島礁條件下的平臺水動力分析方面則極少見到相關(guān)的文章。

本文針對島礁地形下半潛平臺的水動力性能進(jìn)行研究分析,應(yīng)用AQWA軟件對有島礁及無島礁兩種情況開展數(shù)值計(jì)算,島礁地形對平臺水動力及運(yùn)動響應(yīng)的影響通過將島礁模擬成固定物體的濕表面加以考慮。通過對比數(shù)值計(jì)算結(jié)果與模型試驗(yàn)數(shù)據(jù),研究了近島礁條件對平臺水動力及運(yùn)動性能的影響,可為近島礁條件下平臺的性能研究和系泊系統(tǒng)的設(shè)計(jì)提供分析基礎(chǔ)。

1 基本理論

數(shù)值分析的基本理論基于線性理論,在頻域內(nèi)進(jìn)行計(jì)算,考慮島礁地形下的半潛平臺運(yùn)動響應(yīng)分析,總的速度勢可以表示為[11]:

1.1 入射問題

在數(shù)值計(jì)算中,對于設(shè)定水深的波浪入射勢可以表示為:

式中:d為水深,k為波數(shù),θ為浪向,應(yīng)用AQWA的頻域數(shù)值分析基于波幅A為1 m的規(guī)則波。基于線性化后的伯努利方程,可以得到平臺表面的壓力分布,再通過積分得到平臺的入射波浪力(Froude-Krylov力):

1.2 繞射問題

在數(shù)值計(jì)算中,對于島礁地形下平臺的繞射問題,繞射勢的控制方程及邊界條件表示為:

△φD=0 在流域內(nèi)

式中:物面的邊界條件如(6)式所示,入射勢φI已知,繞射勢φD可以表示為如(7)式所示:

式中:G為格林函數(shù):

1.3 輻射問題

島礁地形下平臺的輻射問題的數(shù)值計(jì)算,輻射勢的控制方程及邊界條件表示如下:

△φR=0 在流域內(nèi)

式中:物面的邊界條件在固定的地形表面法向速度為0,在平臺表面則等于平臺的運(yùn)動速度,輻射勢φR可以表示為:

平臺所受到的輻射波浪力可以表示為:

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

2.1 半潛平臺參數(shù)

本文研究的近島礁平臺為半潛形式,由上部甲板、中間立柱、下部浮筒與橫撐組成,平臺示意圖見圖2,平臺的主尺度參數(shù)見表1。

表1 半潛平臺主尺度Tab.1 Main dimensions of the semi-submersible platform

圖2 半潛平臺Fig.2 Semi-submersible platform

2.2 島礁地形模型

地形模型基于島礁實(shí)地地形的測量數(shù)據(jù),考慮到模型單元數(shù)及計(jì)算時(shí)間的局限,在水動力計(jì)算中地形的形狀為200 m(沿平臺寬度方向)×150 m(沿平臺長度方向)的長方形,完全浸沒于水面以下。地形周邊的海底深度為35 m,而地形高度則隨不同的水平位置而變化,在平臺中心位置處的水深為10 m。

圖3 島礁地形下半潛平臺AQWA數(shù)值分析模型Fig.3 AQWA model of the semi-submersible platform near a reef island

在近島礁海域,由于波浪傳播速度隨水深變淺而減小,來自不同方向的波浪沿著礁盤浪花區(qū)傳播后的方向會逐漸趨于垂直于島礁的坡岸,因此研究沿此方向傳播的波浪對平臺運(yùn)動的影響具有重要的意義。圖4為島礁地形下平臺的坐標(biāo)以及運(yùn)動方向定義,與模型試驗(yàn)研究的入射波浪和平臺運(yùn)動的方向定義相一致,波浪沿著X方向入射,平臺的縱蕩沿著Y方向,橫蕩沿著X方向,由此按照右手系法確定平臺相應(yīng)的六自由度運(yùn)動方向。

圖4 島礁地形下半潛平臺坐標(biāo)與運(yùn)動方向定義Fig.4 Coordinates and motion directions definition of the semi-submersible platform near a reef island

數(shù)值計(jì)算在頻域中模擬1 m波幅的規(guī)則波,浪向?yàn)?°沿著X方向。作為對比,無島礁地形下的平臺計(jì)算水深為10 m,而考慮島礁地形的入射波水深為35 m,在平臺布放處的水深為10 m,相關(guān)的數(shù)值計(jì)算條件參見表2。

表2 數(shù)值計(jì)算條件Tab.2 Numerical calculation condition

3 水動力性能分析

3.1 島礁地形對平臺響應(yīng)影響分析

圖5~10所示為有島礁地形和無島礁地形時(shí)半潛平臺運(yùn)動響應(yīng)的對比結(jié)果。圖中基于圓點(diǎn)的擬合曲線所示為無島礁地形條件下單個(gè)平臺的運(yùn)動響應(yīng)RAO數(shù)值規(guī)律[12],基于三角形點(diǎn)的擬合曲線所示為島礁地形下平臺的運(yùn)動響應(yīng)RAO數(shù)值規(guī)律。橫坐標(biāo)為波浪周期,縱蕩(Surge)、橫蕩(Sway)和垂蕩(Heave)的RAO分別為1 m波幅下平臺在該方向的位移,而橫搖(Roll)、縱搖(Pitch)和艏搖(Yaw)的RAO分別表示1 m波幅下平臺在該方向的轉(zhuǎn)動角度。

圖5 平臺縱蕩運(yùn)動響應(yīng)RAOFig.5 Surge RAO of the platform

圖6 平臺橫蕩運(yùn)動響應(yīng)RAOFig.6 Sway RAO of the platform

圖7 平臺垂蕩運(yùn)動響應(yīng)RAOFig.7 Heave RAO of the platform

圖8 平臺橫搖運(yùn)動響應(yīng)RAOFig.8 Roll RAO of the platform

圖9 平臺縱搖運(yùn)動響應(yīng)RAOFig.9 Pitch RAO of the platfor

圖10 平臺艏搖運(yùn)動響應(yīng)RAOFig.10 Yaw RAO of the platform

針對無島礁地形下的半潛平臺,在零度浪向的入射波浪作用下,平臺的縱蕩、縱搖與艏搖運(yùn)動由于平臺關(guān)于入射波的對稱性而接近于0;而由于島礁的存在,波浪受到非均勻地形的影響其方向與波高發(fā)生變化,其關(guān)于平臺的對稱性被破壞,因此位于設(shè)定布放點(diǎn)的半潛平臺在縱蕩、縱搖與艏搖方向均會發(fā)生運(yùn)動,縱蕩在波周期14 s時(shí)產(chǎn)生1.09 m位移,而縱搖與艏搖在波周期18 s時(shí)出現(xiàn)峰值,分別為2.99°與2.59°。

對于受沿X方向的浪向影響最為主要的橫蕩、橫搖與垂蕩運(yùn)動,島礁地形下平臺的運(yùn)動響應(yīng)隨著波周期的變化規(guī)律整體上與無島礁地形的規(guī)律相近,但由于固定的地形對波高的變化產(chǎn)生影響,因此島礁地形下半潛平臺的運(yùn)動響應(yīng)數(shù)值總體偏大,平臺的橫蕩在波周期25 s時(shí)最大為8.05 m,平臺的垂蕩及橫搖運(yùn)動在波周期接近自然周期時(shí)出現(xiàn)響應(yīng)峰值,同時(shí)受到島礁地形對波高的影響,平臺的垂蕩與橫搖運(yùn)動在有地形條件下的峰值數(shù)值約為無地形條件下平臺的2倍,垂蕩運(yùn)動在17 s時(shí)出現(xiàn)峰值為3.47 m,而橫搖在波周期20 s時(shí)運(yùn)動最大為16.25°。

3.2 平臺響應(yīng)與試驗(yàn)結(jié)果對比分析

平臺模型試驗(yàn)在江蘇科技大學(xué)波浪水池進(jìn)行,水池底部按島礁實(shí)測地形數(shù)據(jù)及模型縮尺比建造了假底以模擬地形的影響。圖11為模擬地形的假底的建造示意圖,圖12所示為地形條件下半潛平臺模型試驗(yàn)圖。實(shí)驗(yàn)中進(jìn)行了規(guī)則波作用下平臺運(yùn)動響應(yīng)的試驗(yàn)研究,由此得到平臺六自由度RAO數(shù)值。本文中應(yīng)用AQWA開展數(shù)值分析的島礁地形及平臺布放位置與試驗(yàn)相一致。圖13~18所示為數(shù)值計(jì)算與模型試驗(yàn)所得到的島礁地形下平臺運(yùn)動響應(yīng)的RAO對比圖。

圖11 模型試驗(yàn)的地形搭建Fig.11 Construction of the reef model for experiment

圖12 半潛平臺模型試驗(yàn)Fig.12 Model experiment of the platform

模型試驗(yàn)?zāi)M的波周期范圍為T=4.38~16.26 s。在此范圍內(nèi)數(shù)值計(jì)算與模型試驗(yàn)的運(yùn)動響應(yīng)變化趨勢總體相近。對于受沿X方向的浪向影響最為主要的橫蕩、橫搖與垂蕩運(yùn)動,數(shù)值計(jì)算結(jié)果與模型試驗(yàn)值對比較為吻合;而對于平臺的縱蕩、縱搖和艏搖運(yùn)動,由于地形處沿Y方向斜坡的存在在一定程度上影響波浪的方向及波高變化,對比結(jié)果在波周期10~16 s存在一定的差異,對比試驗(yàn)數(shù)據(jù),數(shù)值計(jì)算的縱蕩值偏大而縱搖值偏小。這意味著數(shù)值模型中波浪沿與來波垂直方向的演化與實(shí)驗(yàn)結(jié)果有一定差異,而造成這種差異的原因,部分可以歸結(jié)為數(shù)值模型中模擬的地形范圍較實(shí)驗(yàn)中構(gòu)造的地形偏小,當(dāng)然也可能是計(jì)算所依據(jù)的線性波理論的局限性。整體而言,應(yīng)用AQWA計(jì)算得出的島礁地形下平臺的運(yùn)動響應(yīng)結(jié)果與模型試驗(yàn)值對比吻合較好。

圖13 平臺縱蕩運(yùn)動響應(yīng)RAO對比Fig.13 Comparison of the Surge RAO

圖14 平臺橫蕩運(yùn)動響應(yīng)RAO對比Fig.14 Comparison of the Sway RAO

圖15 平臺垂蕩運(yùn)動響應(yīng)RAO對比Fig.15 Comparison of the Heave RAO

圖16 平臺橫搖運(yùn)動響應(yīng)RAO對比Fig.16 Comparison of the Roll RAO

圖17 平臺縱搖運(yùn)動響應(yīng)RAO對比Fig.17 Comparison of the Pitch RAO

圖18 平臺艏搖運(yùn)動響應(yīng)RAO對比Fig.18 Comparison of the Yaw RAO

3.3 島礁地形下平臺短期預(yù)報(bào)分析

島礁地形下平臺在不規(guī)則波中的運(yùn)動響應(yīng)通過三小時(shí)短期預(yù)報(bào)得出。入射波應(yīng)用JONSWAP譜,有義波高Hs為3 m,波譜峰周期Tp分別為6 s、8 s、10 s、12 s和14 s,譜峰提升因子γ選取3.0,如表3所示,其中數(shù)值為島礁地形下平臺短期預(yù)報(bào)所得到的運(yùn)動響應(yīng)最大值。

短期預(yù)報(bào)的結(jié)果表明,對于考慮地形條件的平臺,近島礁波浪方向下3 m有義波高的平臺平動運(yùn)動主要體現(xiàn)在橫蕩與垂蕩,短期預(yù)報(bào)得到的運(yùn)動最大值隨著Tp的增大而增大;而對于平臺的轉(zhuǎn)動運(yùn)動,受到來浪的作用,平臺的橫搖短期預(yù)報(bào)的最大值約在4.5°至6.5°范圍,同時(shí)由于島礁地形的特點(diǎn),來浪方向下平臺縱搖的短期預(yù)報(bào)最大值約為2.34°,艏搖運(yùn)動的最大值可達(dá)4.68°。

表3 平臺響應(yīng)短期預(yù)報(bào)Tab.3 Short-term prediction of platform responses

4 結(jié) 論

論文基于一島礁實(shí)地的地形數(shù)據(jù),應(yīng)用AQWA對島礁地形下的半潛式平臺開展數(shù)值計(jì)算和分析,重點(diǎn)研究島礁地形下平臺的運(yùn)動響應(yīng),結(jié)合模型試驗(yàn)數(shù)據(jù)進(jìn)行對比分析,得出以下主要結(jié)論:

(1)應(yīng)用AQWA可以對島礁地形下的半潛式平臺進(jìn)行數(shù)值計(jì)算模擬,所得到的地形條件下平臺的運(yùn)動響應(yīng)結(jié)果與模型試驗(yàn)值對比吻合良好;

(2)對比無島礁地形下平臺的運(yùn)動響應(yīng),島礁地形下平臺的運(yùn)動響應(yīng)數(shù)值及峰值大小總體上偏大,主要原因與入射波受島礁影響波高總體增大及方向變化有關(guān);

(3)與島礁地形下入射浪向相關(guān)的平臺橫蕩、垂蕩與橫搖運(yùn)動響應(yīng)的計(jì)算結(jié)果與試驗(yàn)值對比的吻合度相對較高,而縱蕩、縱搖與艏搖運(yùn)動的計(jì)算值與試驗(yàn)數(shù)值對比存在一定的差異。

本文應(yīng)用數(shù)值分析方法,通過有無島礁地形以及與模型試驗(yàn)結(jié)果的對比,為近島礁條件下平臺的性能研究和系泊系統(tǒng)設(shè)計(jì)提供了一種分析方法。

[1]王世圣,謝 彬,曾恒一,馮 瑋,李曉平,張海濱.3000米深水半潛式鉆井平臺運(yùn)動性能研究[J].中國海上油氣, 2007,19(4):277-284. Wang Shisheng,Xie Bin,Zeng Hengyi,Feng Wei,Zhang Haibin.Study on motion performance of 3 000 meters deepwater semi-submersible drilling platform[J].China Offshore Oil and Gas,2007,19(4):277-284.

[2]朱 航,馬 哲,翟剛軍,歐進(jìn)萍.風(fēng)浪作用下HYSY-981半潛式平臺動力響應(yīng)的數(shù)值模擬[J].振動與沖擊,2010,29 (9):113-118. Zhu Hang,Ma Zhe,Zhai Gangjun,Ou Jinping.Numerical analysis for the HYSY-981 semisubmersible platform motion responses in wind and wave[J].Journal of Vibration and Shock,2010,29(9):113-118.

[3]張 威.深海半潛式鉆井平臺水動力性能分析[D].上海:上海交通大學(xué),2006. Zhang Wei.Hydrodynamic analysis for a deep-water semi-submersible drilling platform[D].Shanghai:Shanghai Jiaotong University,2006.

[4]肖宇維,孫樹民.基于AQWA的輕型張力腿平臺型式研究[J].科學(xué)技術(shù)與工程,2012,12(33):8987-8993. Xiao Yuwei,Sun Shumin.A study of light tension leg platform types based on AQWA[J].Science Technology and Engineering,2012,12(33):8987-8993.

[5]姚彥龍.三角形張力腿平臺運(yùn)動性能研究[D].哈爾濱:哈爾濱工程大學(xué),2013. Yao Yanlong.Triangular tension leg platform motion performance research[D].Harbin:Harbin Engineering University, 2013.

[6]蒙占彬,李 牧,楊樹耕.張力腿平臺運(yùn)動性能分析研究[J].中國海洋平臺,2013,28(2):28-34. Meng Zhanbin,Li Mu,Yang Shugeng.The analysis and research of TLP’s motion performance[J].China Offshore Platform,2013,28(2):28-34.

[7]龍?bào)銜?李 巍,尤云祥.風(fēng)浪流中Truss Spar平臺耦合動力響應(yīng)分析[J].中國海洋平臺,2013,28(1):29-35. Long Xiaoye,Li Wei,You Yunxiang.Coupled dynamic analysis for a deep-water truss spar platform in wind,wave and current[J].China Offshore Platform,2013,28(1):29-35.

[8]Mansour A M,Huang E W.H-shaped pontoon deepwater floating production semisubmersible[C]//Proceedings of the 26th International Conference on Offshore Mechanics and Arctic Engineering.California,USA.OMAE2007-29385.

[9]Clauss G F,Schmittner C E,Stutz K.Freak wave impact on semisubmersibles time-domain analysis of motions and forces [C]//Proceedings of the 13th International Offshore and Polar Engineering Conference,2003,JSC-371.

[10]Lowa Y M,Langleyb R S.Time and frequency domain coupled analysis of deepwater floating production systems[J].Applied Ocean Research,2006,28:371-385.

[11]MOLIN.海洋工程水動力學(xué)[M].北京:國防工業(yè)出版社,2012.

[12]程小明,楊 鵬,吳 波.島礁半潛式生產(chǎn)生活平臺水動力性能分析[R].無錫:中國船舶科學(xué)研究中心科技報(bào)告, 2013. Cheng Xiaoming,Yang Peng,Wu Bo.Investigation on hydrodynamic performance of a semi-submersible platform for production and living[R].Wuxi:China Ship Scientific Research Center Report,2013.

Investigation on hydrodynamic performance of a semi-submersible platform near a reef island

WU Bo,CHENG Xiao-ming,TIAN Chao,DING Jun,LIU Xiao-long
(China Ship Scientific Research Center,Wuxi 214082,China)

The hydrodynamic performance of a semi-submersible platform positioned near a reef island is investigated both numerically and experimentally.The numerical results of the platform motions are compared with experimental data,which indicates that by considering the hydrodynamic interaction between the platform and the reef island in a conventional linear wave theory based on numerical method,the hydrodynamics and motions of the semi-submersible platform near the reef island can be numerically computed with a reasonable agreement with experimental data.

reef island;semi-submersible platform;hydrodynamic performanc;numerical analysis

U656.6

A

10.3969/j.issn.1007-7294.2015.07.003

1007-7294(2015)07-0773-09

2015-04-16

國家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃資助(2013CB036102);工信部高技術(shù)船舶科研項(xiàng)目資助工信部聯(lián)裝[2012]533

吳 波(1988-),男,工程師,E-mail:paulwu_china@163.com;

程小明(1957-),男,研究員,國家千人計(jì)劃學(xué)者。

猜你喜歡
方向分析
2022年組稿方向
2022年組稿方向
2021年組稿方向
2021年組稿方向
隱蔽失效適航要求符合性驗(yàn)證分析
2021年組稿方向
電力系統(tǒng)不平衡分析
電子制作(2018年18期)2018-11-14 01:48:24
電力系統(tǒng)及其自動化發(fā)展趨勢分析
位置與方向
中西醫(yī)結(jié)合治療抑郁癥100例分析
主站蜘蛛池模板: 亚洲欧洲日产无码AV| 国产青榴视频在线观看网站| 999福利激情视频| 日韩天堂在线观看| 亚洲精品午夜无码电影网| 69视频国产| 色偷偷男人的天堂亚洲av| 国产日本一线在线观看免费| 91在线免费公开视频| 粉嫩国产白浆在线观看| 亚洲侵犯无码网址在线观看| 91成人在线观看视频| 思思热在线视频精品| 伊人婷婷色香五月综合缴缴情| 欧美日韩资源| 中文国产成人精品久久一| 国产人人乐人人爱| 国产流白浆视频| 亚洲欧美自拍中文| 97se亚洲综合在线天天| av一区二区三区高清久久| 中文字幕永久视频| 一本大道香蕉中文日本不卡高清二区| 国产一在线| 亚洲国产成人综合精品2020 | 国产欧美日韩va另类在线播放| 亚洲无码高清视频在线观看| 国产欧美日韩va另类在线播放| 草逼视频国产| 无码内射中文字幕岛国片 | 亚洲精品手机在线| 亚洲香蕉伊综合在人在线| 秘书高跟黑色丝袜国产91在线 | 国产十八禁在线观看免费| 久久国产毛片| 欧美无专区| 国产哺乳奶水91在线播放| 成年午夜精品久久精品| 谁有在线观看日韩亚洲最新视频| 国产一区成人| 国产真实乱子伦视频播放| 精品人妻AV区| 亚洲福利片无码最新在线播放| 亚洲日韩精品无码专区97| 国内精品一区二区在线观看| 国产96在线 | 亚洲精品人成网线在线 | 久久青草精品一区二区三区| a级毛片网| 制服丝袜一区二区三区在线| 999国产精品| 91精选国产大片| 亚洲中文制服丝袜欧美精品| 久久国产精品夜色| 国产69精品久久久久孕妇大杂乱 | 精品久久国产综合精麻豆| 国产成人高清精品免费| 亚洲一级毛片免费看| 久久综合色视频| 亚洲色婷婷一区二区| 国产精品视频系列专区| 制服丝袜亚洲| av在线5g无码天天| 国精品91人妻无码一区二区三区| 黄色免费在线网址| 免费无码AV片在线观看中文| 亚洲一区网站| 国产乱人激情H在线观看| 国产免费久久精品99re不卡 | 日韩欧美国产三级| 亚洲欧美不卡视频| 欧美a级在线| 精品国产电影久久九九| 亚洲中文字幕av无码区| 日本欧美午夜| 亚洲综合精品第一页| 91www在线观看| 国产91高跟丝袜| 国产成人精品亚洲日本对白优播| 国产综合网站| 1769国产精品免费视频| 亚洲一区二区三区麻豆|