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

    水下設備懸垂法安裝動力定位耦合數(shù)值模擬

    2017-12-06 06:45:28孫麗萍朱梟猛艾尚茂李彬
    哈爾濱工程大學學報 2017年11期
    關鍵詞:船舶設備系統(tǒng)

    孫麗萍,朱梟猛,艾尚茂,李彬

    (1.哈爾濱工程大學 船舶工程學院,黑龍江 哈爾濱 150001; 2.哈爾濱工程大學 水下機器人技術(shù)重點實驗室,黑龍江 哈爾濱 150001)

    水下設備懸垂法安裝動力定位耦合數(shù)值模擬

    孫麗萍1,朱梟猛1,艾尚茂1,李彬2

    (1.哈爾濱工程大學 船舶工程學院,黑龍江 哈爾濱 150001; 2.哈爾濱工程大學 水下機器人技術(shù)重點實驗室,黑龍江 哈爾濱 150001)

    懸垂法是目前較為先進的深水、超深水大型水下裝備安裝方法,針對考慮安裝船動力定位的懸垂法安裝過程數(shù)值模擬問題,本文結(jié)合三維勢流理論、集中質(zhì)量法以及Morison方程建立了安裝系統(tǒng)分析模型。利用PID控制器模擬動力定位系統(tǒng),設計了推力分配單元,并考慮動力定位對安裝系統(tǒng)的作用,最終實現(xiàn)安裝船動力定位模式下的水下設備懸垂法安裝時域耦合數(shù)值模擬。模擬結(jié)果顯示:懸垂下放過程中安裝纜張力逐漸增加,避免了突變載荷,且下放過程平緩;船的波頻運動對安裝系統(tǒng)影響較大;水下設備在下放初始階段發(fā)生縱搖,且對橫向載荷較為敏感;設備質(zhì)量影響運動速度和縱搖振蕩角度。

    懸垂法; 動力定位; 水下安裝; 時域耦合; 數(shù)值模擬; 推進器; 波頻運動

    海洋油氣開發(fā)依賴于水下生產(chǎn)系統(tǒng),其深水安裝技術(shù)是建設生產(chǎn)系統(tǒng)的基礎。懸垂法(pendulous installation method,PIM)作為一種新型水下設備安裝方法,具有成本低、適用深水、避免共振,且對安裝船舶要求較低,施工方便等優(yōu)點,適合對水下增壓系統(tǒng)、管匯等大型水下設備進行深水安裝[1],然而該方法實際案例較少,仍需進行深入研究。

    Petrobras公司利用數(shù)值模擬、水池試驗研究了管匯的水動力性能以及不同形狀物體下放運動特征及纜索受力等,并進行全尺寸模型試驗,最終利用PIM法成功將重達280 t和200 t的水下管匯安裝至1 845 m及1 900 m水深[2-4]。國內(nèi)學者利用Fluent、Moses等軟件對PIM法安裝過程中船-纜-體耦合系統(tǒng)進行模擬[5-6];對水下分離器PIM法下放過程進行風險分析,建立了風險評估體系[7]。

    懸垂法安裝水深較大,錨泊定位將不適用,且錨泊線易與安裝纜干涉,因此安裝船舶需利用動力定位(dynamic positioning,DP)系統(tǒng)保持位置及艏向,在PIM法安裝中對動力定位過程的模擬具有重要意義。

    本文建立安裝系統(tǒng)數(shù)學模型,然后建立DP系統(tǒng)簡化模型。在起重船動力定位模式下,進行吊放作業(yè)時域耦合數(shù)值模擬。在三用工作船動力定位模式下,進行水下設備懸垂運動下放過程時域耦合數(shù)值模擬。解決了PIM法安裝中吊放作業(yè)及懸垂運動與水面船舶動力定位的耦合問題。得到運動、推力和纜張力等計算結(jié)果,分析了運動相互影響以及安裝纜張力的影響,并進行敏感性分析。

    1 安裝系統(tǒng)數(shù)學模型

    水下設備PIM法安裝作業(yè)涉及相關船舶在吊放以及懸垂下放過程中進行六自由度運動,其運動方程為[8]

    式中:M為包含附加質(zhì)量的船體慣性矩陣,C(ν)為科里奧利向心力矩陣,D(ν)為阻尼矩陣,g(η)為回復力,τ為船舶所受載荷,η為地球固定坐標系下船舶位置,v為隨體坐標系下船舶運動速度,J為坐標轉(zhuǎn)換矩陣。

    通常船舶在海上作業(yè)將受到波浪、風和流的影響,一階波浪力影響船舶波頻運動,二階波浪力、風力、流力和推進器推力影響船舶低頻運動[9]。DP系統(tǒng)主要對船舶低頻運動進行控制。

    一階波浪力采用三維勢流理論進行求解。認為不規(guī)則波是若干規(guī)則波的疊加,設規(guī)則波數(shù)為N,ωj、aj分別為波譜中各規(guī)則波頻率與波幅,fj為與之對應的Froude-Krylov力和繞射力,εj為隨機相位角,則一階波浪力的計算公式為

    采用二次傳遞函數(shù)計算二階波浪力,忽略和頻力,根據(jù)Newman近似有二階波浪力計算公式[10]:

    式中Tjj為傳遞函數(shù)。

    根據(jù)OCIMF規(guī)范推薦方法求解風和流載荷[11]:

    式中:Cx、Cy、Cm為風/流力系數(shù),ρ為流體密度,V為相對速度,Ax、Ay、Am為受力面積。

    安裝纜采用集中質(zhì)量法進行動態(tài)分析,假設纜線由若干集中質(zhì)量點和無質(zhì)量直線段纜元組成[12]。纜元僅模擬纜線軸向和扭轉(zhuǎn)特性,質(zhì)量、重力、浮力、流體作用力等均集中作用于纜元兩端節(jié)點處。

    纜元中軸向彈簧阻尼單元的有效張力為

    纜元間彎矩為

    纜元扭矩為

    式中:Kr為扭轉(zhuǎn)剛度,τ為扭轉(zhuǎn)角,L0為單元初始長度,Dr為扭轉(zhuǎn)阻尼。

    作用于纜線的阻力采用Morison方程求解。

    附加質(zhì)量的影響通常在局部x、y、z方向上分別計算,對于每個方向纜元受纜加速度引起的額外慣性力和周圍流體加速運動的作用力。各方向附加質(zhì)量作用為

    Fa=CaMFAL+(1+Ca)MFAF

    式中:Ca是該方向附加質(zhì)量系數(shù),MF流體質(zhì)量,AL為該方向纜加速度分量,AF為該方向流體加速度分量。

    水下安裝設備進行剛體六自由度運動,流體力采用Morison方程求解。

    水面工作船舶、安裝纜、水下設備組成水下安裝系統(tǒng),該系統(tǒng)的運動方程為

    M(p,a)+C(p,v)+K(p)=F(p,v,t)

    式中:M為系統(tǒng)慣性力,C為系統(tǒng)阻尼力,K為系統(tǒng)剛度力,p、v、a分別為位置、速度和加速度,t為時間。

    2 動力定位系統(tǒng)

    本文采用PID控制器模擬DP控制系統(tǒng):

    DP船舶的運動系統(tǒng)通常是過驅(qū)動的,將推力分配轉(zhuǎn)化為最優(yōu)化問題來求解。目標是在各推進器發(fā)出所需總推力和力矩的同時達到最低的能源消耗,同時減小機械磨損。推力分配目標函數(shù)為[14]

    約束條件為推力及變化率上下限、角度及變化率上下限。其中s=τ-B(α)u,τ為總推力,B為方位角矩陣,α為各推進器方位角,u為各推進器推力。三用工作船主推進器為槳舵組合,并配備三臺槽道推進器。起重船采用七臺全回轉(zhuǎn)推進器,如圖1所示。

    圖1 起重船推進器位置Fig.1 Position of thrusters on crane vessel

    3 懸垂法安裝時域模擬

    PIM法安裝初始狀態(tài)如圖2所示,水下設備被起重船吊放至水深50 m處,三用工作船利用DP系統(tǒng)穩(wěn)定在預定位置。切斷起重船吊纜,安裝設備即開始擺錘下放運動。

    圖2 初始狀態(tài)Fig.2 Initial condition

    安裝設備基本參數(shù)如表1所示。起重船總長205 m,垂線間長185 m,型寬39.2 m,吃水8.1 m, 排水量48 430 t。與三用工作船的主尺度以及推進器配置參數(shù)分別如表2所示。

    表1 安裝設備基本參數(shù)

    表2 三用工作船主要參數(shù)

    安裝纜由三種成分組成,連接安裝船的部分為鏈,中間部分為聚酯纜,下端為鋼纜。這樣既可以保證安裝纜強度,又能降低重量并避免共振。安裝纜的各部分參數(shù)如表3所示。

    表3 安裝纜主要參數(shù)

    該設備目標安裝海域水深1 500 m,風速10 m/s,采用JONSWAP波能譜模擬海浪不規(guī)則波,有義波高1.5 m,譜峰周期6.7 s,各深度流速如表4所示。假設風、浪、流載荷沿同一方向作用于安裝系統(tǒng)。為保證水下懸垂下放安全施工,設環(huán)境載荷方向為180°,即安裝船處于迎浪狀態(tài)。

    表4 海流速度

    4 計算結(jié)果與討論

    4.1吊放過程數(shù)值模擬

    設備通過起重船吊機以0.1 m/s速度下放至水下50 m,之后保持10 min以校核吊纜張力。吊放過程總模擬時間設為1 200 s。

    圖3展示了起重船縱向運動和縱向總推力。起重船在DP系統(tǒng)控制下,縱向運動具有1 m以內(nèi)的定位精度。設備吊放過程對起重船升沉及縱搖平衡位置影響很小,主要由于其重量遠小于起重船排水量。推進器在該海況下處于低功率狀態(tài),最大總推力132 kN,有很大余量保證吊放過程的定位精度。

    圖4展示了水下設備吊放過程的垂向位移、垂向速度、縱向位移和縱向速度。設備垂向運動受到起重船升沉及縱搖影響,而縱向運動則受起重船縱向運動影響,其運動趨勢與起重船相吻合。

    圖3 起重船計算結(jié)果Fig.3 Simulation results of crane vessel

    圖4 水下設備計算結(jié)果Fig.4 Simulation results of subsea equipment

    圖5為吊纜張力曲線,張力隨起重船及水下設備運動而產(chǎn)生動態(tài)響應,作業(yè)過程纜張力處于安全范圍內(nèi),且水下設備入水后浮力增加導致張力降低。

    圖5 吊纜張力Fig.5 Sling tension

    4.2懸垂下放過程數(shù)值模擬

    水下設備懸垂下放過程中,三用工作船處于DP模式。圖6(a)~(c)展示了工作船水平三自由度運動曲線。設備對工作船縱向運動影響較大。下放初始階段由于纜未完全張緊,安裝物對工作船影響較小。安裝纜逐漸張緊時工作船受到縱向外載荷而產(chǎn)生較大位移,達到5.2 m,DP系統(tǒng)具有足夠定位能力,使其位置穩(wěn)定。由圖6(d)、(e)知設備及纜對工作船縱搖及垂蕩平衡位置發(fā)生變化。濾波后得到縱搖平衡位置變?yōu)?0.4°,垂蕩平衡位置變?yōu)?0.1 m。圖7為工作船水平方向總控制力。工作船受環(huán)境載荷及水下設備作用發(fā)生縱向位移,縱向推力隨之大幅增加??v向最大總推力達到751 kN。

    圖8(a)、(b)展示了水下設備懸垂法安裝下擺運動過程中垂向速度與縱向速度時歷。在接近海底時,水下設備開始出現(xiàn)高頻垂向震蕩,這種震蕩主要由于工作船運動導致。下放初期,由于安裝纜并未完全拉緊,工作船的運動對水下設備影響很小。此外,在下放初始階段,受流體及安裝纜共同作用產(chǎn)生力矩導致水下設備出現(xiàn)縱搖振蕩,隨下放過程的進行,縱搖角度逐漸恢復為零,如圖8(c)所示。圖8(d)展示了水下設備的運動軌跡,并對初始階段前33 s進行了逐幀示意,圖中可以發(fā)現(xiàn)水下設備在懸垂下放的初始階段產(chǎn)生了類似羽毛在空氣中下落的飄蕩運動。

    圖6 三用工作船運動曲線Fig.6 Motion curves of AHTS

    圖7 三用工作船推進器總控制力Fig.7 Resulting control forces of AHTS

    PIM法安裝過程中,安裝纜結(jié)構(gòu)的可靠性關系著水下設備安裝過程中的安全性。圖9(a)為安裝纜張力時歷曲線。圖9(b)安裝纜最大張力校核,展示了安裝纜各個長度位置下整個模擬過程中的最大張力,其中0 m位置表示安裝纜上端。分析圖9(a)可知工作船波頻運動對安裝纜張力影響較大,使其在安裝設備穩(wěn)定在水下后出現(xiàn)高頻振蕩。圖9(c)為安裝纜張力經(jīng)傅里葉變換后得到安裝纜張力譜曲線。安裝纜在0.67~1.21 rad/s頻率范圍內(nèi)響應較大,而該頻率對應于工作船高頻運動頻率,印證了懸垂法中工作船的波頻運動對安裝纜張力影響較大的結(jié)論。

    為驗證計算模型的通用性,進行敏感性分析,主要考慮不同設備質(zhì)量、波浪周期以及環(huán)境載荷方向?qū)Y(jié)果的影響。表5為設備質(zhì)量影響,結(jié)果顯示隨著質(zhì)量增加,最大垂向速度增加,最大縱搖角減小。圖10展示了不同波浪周期下安裝纜張力,波浪周期影響張力動態(tài)范圍。圖11展示了環(huán)境載荷方向?qū)υO備水平運動軌跡的影響,設備水下運動對橫向載荷敏感。

    表5 設備質(zhì)量影響

    圖8 水下設備計算結(jié)果Fig.8 Simulation results of subsea equipment

    圖9 安裝纜張力計算結(jié)果Fig.9 Simulation results of line tension

    圖10 波浪周期對纜張力影響Fig.10 Influence of wave period on line tension

    圖11 環(huán)境載荷方向影響Fig.11 Influence of environment load direction

    5 結(jié)論

    1)水下設備縱向運動趨勢與起重船基本一致。

    2)在懸垂下放過程中,水下設備及安裝纜主要影響工作船縱蕩、縱搖及垂蕩運動,安裝纜逐漸張緊使工作船產(chǎn)生較大縱向位移,推進器縱向控制力明顯增大,DP系統(tǒng)將工作船穩(wěn)定在指定位置。

    3)水下設備在下放初始階段產(chǎn)生縱搖振蕩;安裝纜、水下設備受流體阻力作用,其擺錘下落運動過程較為平緩,下放需要一定時間;下放階段對橫向載荷較為敏感;安裝纜張力在下放過程中逐漸增加,避免了突變載荷。

    4)工作船波頻運動導致安裝纜張力產(chǎn)生動態(tài)響應,本文范圍內(nèi)波浪周期增大時纜張力動態(tài)范圍增大。

    5)隨著設備質(zhì)量增加,最大垂向速度增加,最大縱搖角反而減小。

    [1] WANG A M, ZHU S H, ZHU X H, et al. Pendulous installation method and its installation analysis for a deepwater manifold in South China Sea[C]// Proceedings of the 23rd International Offshore and Polar Engineering. Alaska, USA, 2013: 774-784.

    [2] FERNANDE A C, SALES J S Jr. The concomitant model testing approach for the development of the pendulous installation method of heavy devices in deep water[C]// Proceedings of the 26th International Conference on Offshore Mechanics and Arctic Engineering. California, USA, 2007: 489-497.

    [3] FERNANDE A C, MINEIRO F P, RIBEIRO M. Estimation of the added mass and drag coefficients relevant to the pendulous installation method by several experimental and numerical methods[C]// Proceedings of 25th International Conference on Offshore Mechanics and Arctic Engineering. Hamburg, Germany, 2006.

    [4] FERNANDE A C, NEVES C R, SALES J S Jr, et al. Installation of deepwater manifolds by the pendulous method under the light of deep ocean basin model testing and numerical[C]// Proceedings of 25th International Conference on Offshore Mechanics and Arctic Engineering. Hamburg, Germany, 2006.

    [5] 姚麗琳. 深水水下安裝下放運動分析技術(shù)研究[D]. 上海:上海交通大學, 2011.

    YAO Lilin. Research on hydrodynamic analysis of subsea deepwater installation[D]. Shanghai: Shanghai JiaoTong University, 2011.

    [6] 張萬里. 基于多體非線性耦合理論的下垂式安裝方法數(shù)值模擬[D].天津:天津大學, 2012.

    ZHANG Wanli. Numerical anaysis of pendulous installation method based on nonlinear multi-body coupling theory[D]. Tianjin: Tianjin University, 2012.

    [7] 趙夢一, 葛玖浩, 李偉, 等. 超深水水下分離器擺錘法下放風險分析[J]. 中國安全生產(chǎn)科學技術(shù), 2015, 11(9): 151-157.

    ZHAO Mengyi, GE Jiuhao, LI Wei, et al. Risk assessment on pendulous installation of subsea separator in ultra deepwater[J]. Journal of safety science and technology, 2015, 11(9): 151-157.

    [8] FOSSEN T I. Guidance and control of ocean vehicles[M]. UK: John Wiley and Sons, 1994.

    [9] FALTINSEN. Sea loads on ships and offshore structures[M]. UK: Cambridge University Press, 1990.

    [10] S?RENSEN A J. A survey of dynamic positioning control systems[J]. Annual reviews in control, 2011, 35(1): 123-136.

    [11] OCIMF. Prediction of wind and current loads on VLCCs[S]. London: Witherby & Co., 1994.

    [12] 宋林峰, 孫麗萍, 王德軍. 深水S型鋪管托管架-船體-管線耦合分析[J]. 哈爾濱工程大學學報, 2013, 34(4): 451-420.

    SONG Linfeng, SUN Liping, WANG Dejun. Coupling analysis of stinger-lay barge-pipeline of S-lay installation in deep water[J]. Journal of Harbin Engineering University, 2013, 34(4): 451-420.

    [13] WU C L, SUN L P, AI S M, et al. Coupling motion analysis on a dynamic-positioning S-laying[J]. Journal of ship mechanics, 2015, 19(3): 249-263.

    [14] FOSSEN T I. Handbook of marine craft hydrodynamics and motion control[M].West Sussex:John Wiley & Sons, 2011.

    本文引用格式:

    孫麗萍,朱梟猛,艾尚茂,等. 水下設備懸垂法安裝動力定位耦合數(shù)值模擬[J]. 哈爾濱工程大學學報, 2017, 38(11): 1689-1695.

    SUN Liping, ZHU Xiaomeng, AI Shangmao, et al. Coupled numerical simulation of pendulous installation method for subsea equipment with dynamic positioning[J]. Journal of Harbin Engineering University, 2017, 38(11): 1689-1695.

    Couplednumericalsimulationofpendulousinstallationmethodforsubseaequipmentwithdynamicpositioning

    SUN Liping1, ZHU Xiaomeng1, AI Shangmao1, LI Bin2

    (1.College of Shipbuilding Engineering, Harbin Engineering University, Harbin 150001, China; 2.Science and Technology on Underwater Vehicle Laboratory, Harbin Engineering University, Harbin 150001, China)

    The pendulous installation method is advanced for lowering large subsea equipment in deepwater and ultra-deepwater, which has a few successful cases. In order to implement the numerical simulation of pendulous installation considering the dynamic positioning process of vessel, the analysis model was established based on three-dimensional potential theory, lumped mass method and Morison equations. PID controller was employed and the thrust allocation unit was designed in the dynamic positioning system of the vessel. The effect of the dynamic positioning on the installation system was taken into account comprehensively. The coupled time domain simulation of pendulous installation of the equipment was accomplished with the vessels operating in dynamic positioning mode. The results show that the tension increases gradually and sudden load is avoided. The lowering progress is smooth. The wave frequency motion has great influence on the installation system. Pitch occurs at the initial period of the pendulous motion. The motion is sensitive to lateral loads. The velocity and pitch fluttering angle are influenced by the change of mass of equipment.

    pendulous installation method; dynamic positioning; subsea installation; coupled time domain; numerical simulation; thrusters; wave frequency motion

    10.11990/jheu.201608052

    http://www.cnki.net/kcms/detail/23.1390.u.20170427.1321.010.html

    U661

    A

    1006-7043(2017)11-1689-07

    2016-08-25.

    網(wǎng)絡出版日期:2017-04-27.

    國家高技術(shù)船舶科研計劃項目(KY10100150023).

    孫麗萍(1962-), 女, 教授,博士生導師;

    朱梟猛(1988-), 男,博士研究生.

    朱梟猛,E-mail: gavinzhu@126.com.

    猜你喜歡
    船舶設備系統(tǒng)
    計算流體力學在船舶操縱運動仿真中的應用
    諧響應分析在設備減振中的應用
    Smartflower POP 一體式光伏系統(tǒng)
    《船舶》2022 年度征訂啟事
    船舶(2021年4期)2021-09-07 17:32:22
    WJ-700無人機系統(tǒng)
    ZC系列無人機遙感系統(tǒng)
    北京測繪(2020年12期)2020-12-29 01:33:58
    船舶!請加速
    基于MPU6050簡單控制設備
    電子制作(2018年11期)2018-08-04 03:26:08
    連通與提升系統(tǒng)的最后一塊拼圖 Audiolab 傲立 M-DAC mini
    500kV輸變電設備運行維護探討
    美女内射精品一级片tv| 久久久欧美国产精品| 日本黄大片高清| 国内揄拍国产精品人妻在线| 一个人免费在线观看电影| 成人性生交大片免费视频hd| 国产免费又黄又爽又色| 边亲边吃奶的免费视频| 午夜激情久久久久久久| 男人爽女人下面视频在线观看| 国产91av在线免费观看| 大香蕉久久网| 亚洲欧美精品专区久久| 国产精品日韩av在线免费观看| 午夜福利视频1000在线观看| 国产日韩欧美在线精品| 插逼视频在线观看| 黑人高潮一二区| av在线播放精品| 色视频www国产| 午夜免费男女啪啪视频观看| 中文字幕av成人在线电影| 午夜激情久久久久久久| 最近的中文字幕免费完整| 日韩电影二区| 亚洲精品日韩在线中文字幕| 一级二级三级毛片免费看| 少妇人妻精品综合一区二区| 免费看不卡的av| 欧美精品一区二区大全| 精品久久国产蜜桃| 在线免费观看不下载黄p国产| 精品国内亚洲2022精品成人| 国产免费福利视频在线观看| 国产一区二区亚洲精品在线观看| 非洲黑人性xxxx精品又粗又长| 国产欧美日韩精品一区二区| freevideosex欧美| 久久精品夜夜夜夜夜久久蜜豆| 2018国产大陆天天弄谢| 2022亚洲国产成人精品| 精品熟女少妇av免费看| 亚洲欧美精品专区久久| 国产成人91sexporn| 国产亚洲最大av| 久久午夜福利片| 国产麻豆成人av免费视频| 亚洲av不卡在线观看| 国产成人a区在线观看| 免费看不卡的av| 真实男女啪啪啪动态图| 天堂网av新在线| 久久久久免费精品人妻一区二区| 你懂的网址亚洲精品在线观看| 亚洲成人精品中文字幕电影| 免费av毛片视频| 久久久国产一区二区| 成人漫画全彩无遮挡| 亚洲国产精品成人综合色| 国产黄色小视频在线观看| 丰满少妇做爰视频| 能在线免费观看的黄片| 国产精品蜜桃在线观看| 天堂影院成人在线观看| 国产午夜精品一二区理论片| 免费看光身美女| 亚洲成色77777| 日韩一本色道免费dvd| 九九久久精品国产亚洲av麻豆| 亚洲内射少妇av| 国产黄片视频在线免费观看| 简卡轻食公司| 亚洲精品国产av蜜桃| av女优亚洲男人天堂| 日本猛色少妇xxxxx猛交久久| 成人毛片a级毛片在线播放| 日本-黄色视频高清免费观看| 午夜福利在线在线| 一级毛片 在线播放| 午夜免费男女啪啪视频观看| 国产麻豆成人av免费视频| 狠狠精品人妻久久久久久综合| 国产亚洲av片在线观看秒播厂 | 国内精品宾馆在线| 久久久久九九精品影院| 国产亚洲精品久久久com| 如何舔出高潮| 国产淫语在线视频| 天堂av国产一区二区熟女人妻| 小蜜桃在线观看免费完整版高清| 波野结衣二区三区在线| 国产精品无大码| 日韩大片免费观看网站| 最近视频中文字幕2019在线8| 69人妻影院| 性色avwww在线观看| 国产中年淑女户外野战色| 亚洲图色成人| 亚洲av中文字字幕乱码综合| 国产毛片a区久久久久| 午夜福利在线观看免费完整高清在| 国产亚洲5aaaaa淫片| 色吧在线观看| 91久久精品国产一区二区三区| 一级毛片久久久久久久久女| 欧美日韩国产mv在线观看视频 | 亚洲欧美日韩无卡精品| 搡女人真爽免费视频火全软件| 淫秽高清视频在线观看| 免费观看无遮挡的男女| 麻豆成人午夜福利视频| 亚洲国产成人一精品久久久| 精品国内亚洲2022精品成人| 久久鲁丝午夜福利片| 日本-黄色视频高清免费观看| 亚洲精品,欧美精品| 久久久久久久久久黄片| 日本免费a在线| 亚洲va在线va天堂va国产| 日韩一区二区三区影片| 免费黄网站久久成人精品| 免费少妇av软件| 特级一级黄色大片| .国产精品久久| 亚洲精品日韩在线中文字幕| 久久精品国产自在天天线| 久久久久久久午夜电影| 性插视频无遮挡在线免费观看| 老司机影院成人| 午夜激情福利司机影院| 亚洲精品国产av蜜桃| 久久这里有精品视频免费| 亚洲色图av天堂| 男的添女的下面高潮视频| 精品久久久久久久久亚洲| 国产黄a三级三级三级人| 午夜福利在线观看吧| 亚洲人与动物交配视频| 午夜久久久久精精品| 十八禁网站网址无遮挡 | 欧美精品一区二区大全| 菩萨蛮人人尽说江南好唐韦庄| 97在线视频观看| 毛片一级片免费看久久久久| 一级黄片播放器| 精品久久久久久久久久久久久| 日韩欧美一区视频在线观看 | 一区二区三区高清视频在线| 热99在线观看视频| a级毛片免费高清观看在线播放| 一级毛片我不卡| 亚洲成人一二三区av| 免费观看精品视频网站| 欧美人与善性xxx| 午夜福利高清视频| 国产 亚洲一区二区三区 | 欧美日韩视频高清一区二区三区二| 国产亚洲精品久久久com| 国产单亲对白刺激| 欧美一级a爱片免费观看看| 91精品一卡2卡3卡4卡| 美女脱内裤让男人舔精品视频| 免费高清在线观看视频在线观看| 久久久久久伊人网av| 欧美另类一区| 午夜老司机福利剧场| 婷婷色av中文字幕| 国产日韩欧美在线精品| 国产精品一区二区三区四区免费观看| 国产女主播在线喷水免费视频网站 | 99热这里只有是精品在线观看| 高清毛片免费看| 可以在线观看毛片的网站| 精品酒店卫生间| 日韩一本色道免费dvd| 五月玫瑰六月丁香| 午夜老司机福利剧场| 99热这里只有精品一区| 国产美女午夜福利| 国产在线一区二区三区精| 中国美白少妇内射xxxbb| 亚洲精品aⅴ在线观看| 99热网站在线观看| 丝袜美腿在线中文| 青春草国产在线视频| 欧美一级a爱片免费观看看| h日本视频在线播放| 国产午夜精品久久久久久一区二区三区| 免费电影在线观看免费观看| 欧美激情在线99| 日本黄色片子视频| 国产激情偷乱视频一区二区| 伦理电影大哥的女人| 国产精品无大码| 有码 亚洲区| 卡戴珊不雅视频在线播放| 国产一级毛片在线| 中文字幕久久专区| 亚洲av二区三区四区| 亚洲精品日韩av片在线观看| 综合色av麻豆| 精品欧美国产一区二区三| 成年av动漫网址| 伊人久久精品亚洲午夜| 国产黄频视频在线观看| 精品久久久久久久久av| 男人和女人高潮做爰伦理| 欧美极品一区二区三区四区| 国产精品不卡视频一区二区| 男女啪啪激烈高潮av片| 精品一区二区免费观看| 亚洲精品久久午夜乱码| 亚洲精品成人av观看孕妇| 精品久久久久久久久久久久久| 伦理电影大哥的女人| 亚洲国产精品sss在线观看| 国产精品日韩av在线免费观看| 午夜福利视频1000在线观看| 精品久久久久久成人av| 国产成人午夜福利电影在线观看| 大陆偷拍与自拍| 91av网一区二区| 一级毛片电影观看| 2021少妇久久久久久久久久久| 中文字幕制服av| 免费看不卡的av| 麻豆久久精品国产亚洲av| 亚洲欧美成人综合另类久久久| 狂野欧美激情性xxxx在线观看| 国产精品嫩草影院av在线观看| 色尼玛亚洲综合影院| 九草在线视频观看| 日日啪夜夜撸| 精品久久国产蜜桃| av又黄又爽大尺度在线免费看| 日日干狠狠操夜夜爽| 午夜亚洲福利在线播放| 国产片特级美女逼逼视频| 日日啪夜夜爽| 99re6热这里在线精品视频| 国内精品美女久久久久久| 国产一区二区三区av在线| 高清av免费在线| 久久午夜福利片| 韩国高清视频一区二区三区| 亚洲,欧美,日韩| 亚洲欧洲国产日韩| 男插女下体视频免费在线播放| 久久精品熟女亚洲av麻豆精品 | 男人舔女人下体高潮全视频| 亚洲精品久久午夜乱码| 免费在线观看成人毛片| 国产综合懂色| 肉色欧美久久久久久久蜜桃 | 亚洲自拍偷在线| or卡值多少钱| 久久草成人影院| 色5月婷婷丁香| 日韩一本色道免费dvd| 国产精品一区二区三区四区久久| 嫩草影院入口| 啦啦啦啦在线视频资源| 日韩电影二区| 老司机影院毛片| 久久99精品国语久久久| 国产亚洲午夜精品一区二区久久 | 久久6这里有精品| 国产精品美女特级片免费视频播放器| 日本一本二区三区精品| 麻豆成人av视频| 街头女战士在线观看网站| 偷拍熟女少妇极品色| a级毛色黄片| 精品少妇黑人巨大在线播放| 日韩精品青青久久久久久| 亚洲伊人久久精品综合| 女人十人毛片免费观看3o分钟| 久久精品国产亚洲网站| 亚洲激情五月婷婷啪啪| 精品人妻视频免费看| 日本免费a在线| 男人狂女人下面高潮的视频| 午夜视频国产福利| 国产午夜福利久久久久久| 自拍偷自拍亚洲精品老妇| 一级a做视频免费观看| 男女视频在线观看网站免费| 嘟嘟电影网在线观看| 国产激情偷乱视频一区二区| 在线a可以看的网站| 日本色播在线视频| 亚洲精品国产成人久久av| 午夜老司机福利剧场| 中文字幕亚洲精品专区| 97超碰精品成人国产| 欧美日韩一区二区视频在线观看视频在线 | 噜噜噜噜噜久久久久久91| 春色校园在线视频观看| videos熟女内射| 三级经典国产精品| 欧美日韩精品成人综合77777| 国产精品女同一区二区软件| 久久人人爽人人片av| 看黄色毛片网站| 国产黄色免费在线视频| 日韩强制内射视频| 亚洲成人久久爱视频| 小蜜桃在线观看免费完整版高清| 欧美精品一区二区大全| 久久久久久久久大av| 亚洲婷婷狠狠爱综合网| 亚洲国产精品国产精品| 成人亚洲欧美一区二区av| 亚洲av.av天堂| 性色avwww在线观看| 国产大屁股一区二区在线视频| 一个人看的www免费观看视频| 国产一区二区在线观看日韩| 久久精品久久精品一区二区三区| 免费在线观看成人毛片| 内射极品少妇av片p| 日本与韩国留学比较| 国产精品爽爽va在线观看网站| 日本黄色片子视频| 色5月婷婷丁香| 91aial.com中文字幕在线观看| 一级毛片久久久久久久久女| 欧美3d第一页| 午夜视频国产福利| 午夜福利高清视频| 久久久久久久午夜电影| 人妻系列 视频| 内射极品少妇av片p| 亚洲精品成人久久久久久| 日本av手机在线免费观看| 成人鲁丝片一二三区免费| kizo精华| 亚洲熟女精品中文字幕| av在线天堂中文字幕| 午夜福利成人在线免费观看| 成人高潮视频无遮挡免费网站| 久久久久精品久久久久真实原创| 中文乱码字字幕精品一区二区三区 | .国产精品久久| 亚洲精品乱久久久久久| 亚洲最大成人手机在线| 夜夜爽夜夜爽视频| 国内少妇人妻偷人精品xxx网站| 青青草视频在线视频观看| 蜜桃亚洲精品一区二区三区| 欧美人与善性xxx| 国产不卡一卡二| videossex国产| 久久久精品免费免费高清| 日产精品乱码卡一卡2卡三| 日韩强制内射视频| 99热这里只有是精品在线观看| 亚洲在久久综合| 亚洲欧洲国产日韩| 大片免费播放器 马上看| 蜜桃久久精品国产亚洲av| 三级经典国产精品| 最近最新中文字幕免费大全7| 久久午夜福利片| 国产免费又黄又爽又色| 亚洲激情五月婷婷啪啪| 色综合站精品国产| 精品久久久久久成人av| 国产免费福利视频在线观看| 日韩三级伦理在线观看| 亚洲欧美一区二区三区国产| 久久久欧美国产精品| 亚洲欧美成人综合另类久久久| 国模一区二区三区四区视频| 国产亚洲av片在线观看秒播厂 | 精品久久久久久久久亚洲| 国精品久久久久久国模美| 超碰97精品在线观看| h日本视频在线播放| 国产精品一区二区三区四区久久| 麻豆成人午夜福利视频| 麻豆成人av视频| 菩萨蛮人人尽说江南好唐韦庄| 2021少妇久久久久久久久久久| 久久这里有精品视频免费| 亚洲欧美精品专区久久| 男女视频在线观看网站免费| 又黄又爽又刺激的免费视频.| 22中文网久久字幕| 亚洲av电影在线观看一区二区三区 | 高清av免费在线| 99热6这里只有精品| 一个人看的www免费观看视频| 亚洲欧洲国产日韩| 午夜精品在线福利| 天天一区二区日本电影三级| 亚州av有码| 亚洲精品第二区| 亚洲av二区三区四区| 一级a做视频免费观看| 国产一级毛片在线| 免费观看av网站的网址| 最近视频中文字幕2019在线8| 亚洲久久久久久中文字幕| 白带黄色成豆腐渣| 最近最新中文字幕免费大全7| 麻豆乱淫一区二区| .国产精品久久| 亚洲18禁久久av| 干丝袜人妻中文字幕| 大香蕉久久网| 精品午夜福利在线看| 国产单亲对白刺激| 成人毛片60女人毛片免费| 成年免费大片在线观看| 久久亚洲国产成人精品v| 亚洲精品日韩av片在线观看| 亚洲精品国产av成人精品| 国产伦一二天堂av在线观看| 一级片'在线观看视频| 久久精品国产亚洲网站| 黄色日韩在线| 一边亲一边摸免费视频| 婷婷色av中文字幕| 在线天堂最新版资源| 日本免费a在线| 日韩欧美国产在线观看| 日本-黄色视频高清免费观看| 久久久久久国产a免费观看| 亚洲美女搞黄在线观看| 我的老师免费观看完整版| 青春草亚洲视频在线观看| 中文字幕免费在线视频6| 永久网站在线| 狂野欧美激情性xxxx在线观看| 欧美日韩亚洲高清精品| 国产有黄有色有爽视频| 久久久久久久午夜电影| 欧美 日韩 精品 国产| 免费观看av网站的网址| or卡值多少钱| 久久久久久久久久久丰满| 国产永久视频网站| 色播亚洲综合网| 可以在线观看毛片的网站| 国产不卡一卡二| 国产精品久久久久久精品电影小说 | 一本一本综合久久| 三级国产精品片| 欧美日韩综合久久久久久| 国产高清国产精品国产三级 | 国产精品一及| 男人狂女人下面高潮的视频| 国产免费福利视频在线观看| 一个人免费在线观看电影| 联通29元200g的流量卡| 男女下面进入的视频免费午夜| 国产精品久久久久久精品电影小说 | 日韩av不卡免费在线播放| 丰满少妇做爰视频| 亚洲无线观看免费| 美女黄网站色视频| 超碰97精品在线观看| 久久久久久九九精品二区国产| 一区二区三区四区激情视频| 赤兔流量卡办理| 九九久久精品国产亚洲av麻豆| 久久人人爽人人片av| 三级国产精品欧美在线观看| 一区二区三区四区激情视频| 成人午夜精彩视频在线观看| 精品一区二区免费观看| 久久99热6这里只有精品| 日本黄大片高清| 亚洲国产最新在线播放| 久久久久久久久中文| 国产 亚洲一区二区三区 | 国产美女午夜福利| .国产精品久久| 久久久精品欧美日韩精品| 在线观看人妻少妇| 国产精品福利在线免费观看| 亚洲va在线va天堂va国产| 亚洲欧美日韩东京热| 国产av不卡久久| 欧美高清性xxxxhd video| 亚洲精品自拍成人| 最后的刺客免费高清国语| 久久久久免费精品人妻一区二区| 免费不卡的大黄色大毛片视频在线观看 | 精品久久久久久久久久久久久| 久久精品人妻少妇| 午夜福利在线观看吧| 亚洲,欧美,日韩| 天堂√8在线中文| 国产91av在线免费观看| 久久久久久九九精品二区国产| 毛片一级片免费看久久久久| av又黄又爽大尺度在线免费看| 亚洲精品第二区| 国产熟女欧美一区二区| 国产综合精华液| 可以在线观看毛片的网站| 精品国产三级普通话版| 午夜激情欧美在线| 在线天堂最新版资源| 久久99热这里只频精品6学生| 一级毛片aaaaaa免费看小| 久久99热这里只有精品18| 三级国产精品片| 亚洲精品视频女| 成年人午夜在线观看视频 | 中文欧美无线码| 午夜精品在线福利| www.av在线官网国产| 美女xxoo啪啪120秒动态图| 中文字幕av在线有码专区| 一边亲一边摸免费视频| 亚洲人成网站高清观看| 日日撸夜夜添| 成人亚洲精品av一区二区| 自拍偷自拍亚洲精品老妇| 国产精品一二三区在线看| 国内精品美女久久久久久| 久久久久九九精品影院| 岛国毛片在线播放| 99九九线精品视频在线观看视频| av专区在线播放| 国产一区亚洲一区在线观看| 十八禁国产超污无遮挡网站| xxx大片免费视频| 国产大屁股一区二区在线视频| 99热网站在线观看| 亚洲不卡免费看| av国产免费在线观看| 天堂俺去俺来也www色官网 | 国产精品久久久久久精品电影| 亚洲人成网站在线播| 哪个播放器可以免费观看大片| 麻豆乱淫一区二区| 91午夜精品亚洲一区二区三区| 中文乱码字字幕精品一区二区三区 | 亚洲精品,欧美精品| 亚洲欧美一区二区三区国产| 午夜福利网站1000一区二区三区| 亚洲av免费高清在线观看| 大陆偷拍与自拍| 丝袜美腿在线中文| 国产白丝娇喘喷水9色精品| 91在线精品国自产拍蜜月| 国产探花在线观看一区二区| 精品久久久久久电影网| 蜜桃久久精品国产亚洲av| 哪个播放器可以免费观看大片| 色视频www国产| 街头女战士在线观看网站| 久久精品久久久久久久性| 蜜臀久久99精品久久宅男| 国产中年淑女户外野战色| 在线观看av片永久免费下载| 熟女人妻精品中文字幕| 两个人的视频大全免费| 好男人在线观看高清免费视频| 亚洲电影在线观看av| 亚洲最大成人手机在线| 最近中文字幕2019免费版| 午夜福利在线观看吧| 一个人看的www免费观看视频| 丝袜美腿在线中文| 黄色一级大片看看| 男女视频在线观看网站免费| 三级男女做爰猛烈吃奶摸视频| 久久久久九九精品影院| 全区人妻精品视频| 免费观看无遮挡的男女| 大香蕉久久网| 日日摸夜夜添夜夜爱| 精品酒店卫生间| 亚洲国产精品成人综合色| www.av在线官网国产| 免费观看无遮挡的男女| 国产毛片a区久久久久| 女人十人毛片免费观看3o分钟| 99re6热这里在线精品视频| 校园人妻丝袜中文字幕| 看免费成人av毛片| 91在线精品国自产拍蜜月| 国产成人精品福利久久| 男女视频在线观看网站免费| 国产成人精品福利久久| 在线观看免费高清a一片| 精品不卡国产一区二区三区| 性色avwww在线观看| 一级毛片黄色毛片免费观看视频| 2021少妇久久久久久久久久久| 精品一区二区三区人妻视频| 国内精品宾馆在线| 国产免费福利视频在线观看| 国产精品熟女久久久久浪| 亚洲精品aⅴ在线观看| 成人av在线播放网站| 日韩一本色道免费dvd| av女优亚洲男人天堂| 日韩av不卡免费在线播放| 亚洲精品成人av观看孕妇| 日本与韩国留学比较| 男人狂女人下面高潮的视频| 免费播放大片免费观看视频在线观看| 精品人妻一区二区三区麻豆| 亚洲av免费在线观看| 成人一区二区视频在线观看| 国产一区二区亚洲精品在线观看| 欧美高清成人免费视频www| 国产综合精华液| 欧美激情久久久久久爽电影|