• <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輸變電設備運行維護探討
    亚洲内射少妇av| a级一级毛片免费在线观看| 老司机影院成人| 亚洲欧美一区二区三区国产| 狂野欧美激情性bbbbbb| 国产精品一区www在线观看| 日韩 亚洲 欧美在线| 天堂俺去俺来也www色官网| 国产精品久久久久久久久免| 高清在线视频一区二区三区| av专区在线播放| 午夜视频国产福利| 99久久精品一区二区三区| 久久精品国产a三级三级三级| 丝袜在线中文字幕| 99热这里只有是精品在线观看| 亚洲精品国产av蜜桃| 久久久久久人妻| 国产免费又黄又爽又色| 夜夜爽夜夜爽视频| 少妇的逼水好多| 久久精品国产亚洲av涩爱| 久久女婷五月综合色啪小说| 欧美精品一区二区大全| 国产午夜精品久久久久久一区二区三区| 26uuu在线亚洲综合色| av免费在线看不卡| 在线观看国产h片| 国产亚洲欧美精品永久| www.色视频.com| 亚洲国产精品专区欧美| 国产一区二区在线观看av| 啦啦啦在线观看免费高清www| av天堂久久9| 少妇猛男粗大的猛烈进出视频| 交换朋友夫妻互换小说| 国产永久视频网站| 久久久亚洲精品成人影院| 99视频精品全部免费 在线| 精品酒店卫生间| 永久网站在线| 午夜日本视频在线| 国内揄拍国产精品人妻在线| 一边亲一边摸免费视频| 人人澡人人妻人| 久久韩国三级中文字幕| 麻豆成人午夜福利视频| 国产无遮挡羞羞视频在线观看| 日本91视频免费播放| 80岁老熟妇乱子伦牲交| 三级经典国产精品| 亚洲精品亚洲一区二区| 草草在线视频免费看| 国产精品一区二区在线不卡| 国产亚洲精品久久久com| 亚洲精品中文字幕在线视频 | 自拍偷自拍亚洲精品老妇| 日韩强制内射视频| 麻豆精品久久久久久蜜桃| 赤兔流量卡办理| 日韩一本色道免费dvd| 91久久精品国产一区二区三区| 少妇人妻精品综合一区二区| 亚洲性久久影院| 街头女战士在线观看网站| 国产亚洲91精品色在线| 乱系列少妇在线播放| 精品一品国产午夜福利视频| 只有这里有精品99| 一级毛片久久久久久久久女| 亚洲人成网站在线播| 伊人久久精品亚洲午夜| 亚洲伊人久久精品综合| 国产在视频线精品| 秋霞伦理黄片| 在线观看免费日韩欧美大片 | 国产精品一二三区在线看| 国产伦理片在线播放av一区| 26uuu在线亚洲综合色| 精品少妇内射三级| 午夜影院在线不卡| 午夜日本视频在线| 国产高清有码在线观看视频| 亚洲国产精品国产精品| 国产高清不卡午夜福利| 日韩av免费高清视频| 曰老女人黄片| 尾随美女入室| 视频区图区小说| 丰满少妇做爰视频| 人妻 亚洲 视频| 一级,二级,三级黄色视频| 亚洲欧洲精品一区二区精品久久久 | 日产精品乱码卡一卡2卡三| 丰满迷人的少妇在线观看| 婷婷色av中文字幕| 波野结衣二区三区在线| 国产成人精品婷婷| 天堂俺去俺来也www色官网| 伊人久久精品亚洲午夜| 最近最新中文字幕免费大全7| 我要看黄色一级片免费的| 大又大粗又爽又黄少妇毛片口| 交换朋友夫妻互换小说| 午夜免费鲁丝| 十八禁高潮呻吟视频 | 国产日韩欧美视频二区| 黄色配什么色好看| 大码成人一级视频| 在线观看www视频免费| 黄色一级大片看看| xxx大片免费视频| 国产 精品1| 国产又色又爽无遮挡免| 久久精品久久久久久噜噜老黄| 卡戴珊不雅视频在线播放| 国产乱人偷精品视频| 永久网站在线| 人人澡人人妻人| 熟妇人妻不卡中文字幕| 国产免费一级a男人的天堂| 国国产精品蜜臀av免费| 国产69精品久久久久777片| 国产在线视频一区二区| 国产精品国产三级国产av玫瑰| 国产精品嫩草影院av在线观看| 人妻 亚洲 视频| 91aial.com中文字幕在线观看| 91久久精品国产一区二区成人| 亚洲精品一区蜜桃| 久久免费观看电影| 美女福利国产在线| 亚洲丝袜综合中文字幕| 欧美bdsm另类| 下体分泌物呈黄色| 中文资源天堂在线| 色吧在线观看| 婷婷色av中文字幕| 亚洲第一区二区三区不卡| 有码 亚洲区| 国产精品一二三区在线看| 少妇高潮的动态图| 国产亚洲欧美精品永久| 亚洲国产成人一精品久久久| 国产av码专区亚洲av| 国产成人精品婷婷| 狠狠精品人妻久久久久久综合| 99精国产麻豆久久婷婷| 成人影院久久| 久久精品国产a三级三级三级| 三级国产精品欧美在线观看| 观看美女的网站| 午夜精品国产一区二区电影| 国内少妇人妻偷人精品xxx网站| 综合色丁香网| 国产日韩欧美在线精品| 免费观看a级毛片全部| 精品少妇内射三级| 少妇被粗大的猛进出69影院 | 黄色日韩在线| 日韩强制内射视频| 91成人精品电影| 91久久精品电影网| 黄色毛片三级朝国网站 | 人人澡人人妻人| 免费在线观看成人毛片| 免费av不卡在线播放| 五月玫瑰六月丁香| 国产精品久久久久久久久免| 欧美激情极品国产一区二区三区 | 熟妇人妻不卡中文字幕| 色5月婷婷丁香| 在线看a的网站| 亚洲自偷自拍三级| 亚洲精品国产av成人精品| 亚洲国产欧美日韩在线播放 | 日本爱情动作片www.在线观看| 天堂俺去俺来也www色官网| 国产精品麻豆人妻色哟哟久久| av福利片在线观看| 亚洲av成人精品一二三区| 91久久精品国产一区二区成人| 亚洲av不卡在线观看| 免费大片黄手机在线观看| 成人黄色视频免费在线看| 午夜免费鲁丝| 欧美成人精品欧美一级黄| 久久久久久久久久成人| 两个人的视频大全免费| 亚洲自偷自拍三级| 久久 成人 亚洲| 国产精品99久久99久久久不卡 | 麻豆精品久久久久久蜜桃| 日韩欧美一区视频在线观看 | 男人爽女人下面视频在线观看| 亚洲精品乱码久久久v下载方式| 久热这里只有精品99| 亚洲电影在线观看av| 九九久久精品国产亚洲av麻豆| 卡戴珊不雅视频在线播放| 三上悠亚av全集在线观看 | 亚洲,一卡二卡三卡| 中文字幕制服av| 欧美激情极品国产一区二区三区 | 97超碰精品成人国产| 久久久亚洲精品成人影院| av一本久久久久| 亚洲欧美日韩另类电影网站| 久久毛片免费看一区二区三区| 国产精品国产三级专区第一集| 又爽又黄a免费视频| 久久久久久久精品精品| 久久国内精品自在自线图片| 中文在线观看免费www的网站| 国产av码专区亚洲av| 日韩制服骚丝袜av| 一区二区三区四区激情视频| 欧美日韩亚洲高清精品| 中国美白少妇内射xxxbb| 国产午夜精品久久久久久一区二区三区| 国产精品女同一区二区软件| 91精品国产九色| 欧美日本中文国产一区发布| 99热网站在线观看| 中文乱码字字幕精品一区二区三区| 久久狼人影院| 精品久久久久久电影网| 久久人人爽人人片av| 99精国产麻豆久久婷婷| 91久久精品国产一区二区三区| 建设人人有责人人尽责人人享有的| 亚洲av成人精品一区久久| 精品卡一卡二卡四卡免费| 欧美丝袜亚洲另类| 久久亚洲国产成人精品v| 伦理电影免费视频| 最近的中文字幕免费完整| 亚洲精品一二三| 中文天堂在线官网| 夫妻午夜视频| 国模一区二区三区四区视频| 91精品国产九色| 国产在线免费精品| 又大又黄又爽视频免费| 最近中文字幕高清免费大全6| 久久人人爽人人片av| 中文字幕亚洲精品专区| 熟女电影av网| 亚洲国产av新网站| 国产成人freesex在线| 最新中文字幕久久久久| a级毛片免费高清观看在线播放| 亚洲va在线va天堂va国产| 观看免费一级毛片| 国产一区二区在线观看日韩| 99久久精品国产国产毛片| 国产精品人妻久久久久久| 最新中文字幕久久久久| 亚洲精品成人av观看孕妇| 女性生殖器流出的白浆| 亚洲精品国产色婷婷电影| 国产爽快片一区二区三区| 美女脱内裤让男人舔精品视频| 亚洲精品乱久久久久久| 狠狠精品人妻久久久久久综合| 三级国产精品欧美在线观看| 男男h啪啪无遮挡| 成人午夜精彩视频在线观看| 少妇人妻久久综合中文| 国产精品伦人一区二区| 久久免费观看电影| 中文字幕精品免费在线观看视频 | 国产极品粉嫩免费观看在线 | 超碰97精品在线观看| 中文字幕人妻丝袜制服| 国产成人免费无遮挡视频| a级一级毛片免费在线观看| 亚洲美女视频黄频| 国产精品久久久久久精品电影小说| 99热6这里只有精品| 高清视频免费观看一区二区| 国产一区亚洲一区在线观看| 乱人伦中国视频| 国产精品一区二区在线观看99| 九九久久精品国产亚洲av麻豆| a 毛片基地| av在线老鸭窝| 欧美日韩视频精品一区| 国产精品久久久久久精品电影小说| 人人妻人人澡人人爽人人夜夜| 9色porny在线观看| 国产在线男女| 精品午夜福利在线看| 特大巨黑吊av在线直播| 免费高清在线观看视频在线观看| 久久国内精品自在自线图片| 国产高清不卡午夜福利| 国精品久久久久久国模美| 久久久精品免费免费高清| 少妇人妻久久综合中文| 亚洲美女搞黄在线观看| 成人黄色视频免费在线看| 男人和女人高潮做爰伦理| 国产欧美亚洲国产| 久久人人爽av亚洲精品天堂| 国产午夜精品一二区理论片| 大片电影免费在线观看免费| 丝袜喷水一区| 精品一区二区免费观看| 中文字幕制服av| 国产黄频视频在线观看| a 毛片基地| 午夜福利在线观看免费完整高清在| 日韩成人av中文字幕在线观看| 日韩电影二区| 丰满人妻一区二区三区视频av| 国产精品蜜桃在线观看| 亚洲精品,欧美精品| 精品一区二区三区视频在线| 插逼视频在线观看| 精品卡一卡二卡四卡免费| 内射极品少妇av片p| 男女边吃奶边做爰视频| 亚洲av成人精品一区久久| 少妇人妻 视频| 丰满饥渴人妻一区二区三| 中文天堂在线官网| 五月天丁香电影| 国产男女超爽视频在线观看| 在线观看免费高清a一片| 日韩成人伦理影院| 亚洲欧洲精品一区二区精品久久久 | 天天躁夜夜躁狠狠久久av| 国产在视频线精品| 大陆偷拍与自拍| 日本爱情动作片www.在线观看| 蜜桃在线观看..| 午夜激情久久久久久久| 免费少妇av软件| 色婷婷av一区二区三区视频| 亚洲不卡免费看| 日本色播在线视频| 亚洲精品中文字幕在线视频 | 春色校园在线视频观看| 欧美区成人在线视频| 免费观看无遮挡的男女| 日本黄色日本黄色录像| av天堂久久9| 久久婷婷青草| 国产精品免费大片| av专区在线播放| 在线观看三级黄色| 日本欧美视频一区| 在线观看三级黄色| 精品一品国产午夜福利视频| 日韩av免费高清视频| 一个人免费看片子| 各种免费的搞黄视频| 国产乱人偷精品视频| 日韩欧美一区视频在线观看 | 两个人免费观看高清视频 | 久久久久久久国产电影| 亚洲熟女精品中文字幕| 免费看日本二区| 国产真实伦视频高清在线观看| 日本黄色片子视频| 黄片无遮挡物在线观看| 超碰97精品在线观看| 欧美精品一区二区免费开放| 嘟嘟电影网在线观看| 国产精品福利在线免费观看| 妹子高潮喷水视频| 亚洲熟女精品中文字幕| 国产91av在线免费观看| 欧美日韩视频精品一区| 狠狠精品人妻久久久久久综合| av福利片在线| 久久久久久人妻| 少妇被粗大猛烈的视频| 亚洲av在线观看美女高潮| 国产亚洲5aaaaa淫片| 51国产日韩欧美| 成人毛片a级毛片在线播放| 日本欧美国产在线视频| 伊人亚洲综合成人网| 国产淫片久久久久久久久| 寂寞人妻少妇视频99o| 中文字幕精品免费在线观看视频 | 国产黄片美女视频| 国产免费一区二区三区四区乱码| 国产av精品麻豆| 久久午夜福利片| 久久久久久久久久久免费av| 久久午夜福利片| 99精国产麻豆久久婷婷| 一本一本综合久久| 国产免费一级a男人的天堂| 老司机亚洲免费影院| 亚洲欧美精品自产自拍| av女优亚洲男人天堂| 国产免费一级a男人的天堂| 精品一区二区三区视频在线| 亚洲av国产av综合av卡| 建设人人有责人人尽责人人享有的| 婷婷色麻豆天堂久久| 51国产日韩欧美| 汤姆久久久久久久影院中文字幕| 性色av一级| 街头女战士在线观看网站| 久久精品久久久久久久性| 久久99蜜桃精品久久| 美女视频免费永久观看网站| 日韩伦理黄色片| 99热这里只有精品一区| 亚洲伊人久久精品综合| 日韩视频在线欧美| 国国产精品蜜臀av免费| 熟女av电影| av天堂久久9| 99re6热这里在线精品视频| 亚洲色图综合在线观看| 国产精品一区二区在线观看99| 最后的刺客免费高清国语| 99视频精品全部免费 在线| 国语对白做爰xxxⅹ性视频网站| 美女脱内裤让男人舔精品视频| 久久国内精品自在自线图片| 一级,二级,三级黄色视频| 天堂俺去俺来也www色官网| 中文欧美无线码| 91精品国产九色| 免费播放大片免费观看视频在线观看| av黄色大香蕉| 亚州av有码| 日本欧美视频一区| 国产精品一区二区在线不卡| 黄色配什么色好看| 日韩欧美精品免费久久| 成人亚洲精品一区在线观看| 观看美女的网站| 人妻一区二区av| 啦啦啦视频在线资源免费观看| 三级经典国产精品| 久久国内精品自在自线图片| 久久热精品热| 热re99久久精品国产66热6| 免费观看a级毛片全部| 国产91av在线免费观看| 日韩中文字幕视频在线看片| 亚洲精华国产精华液的使用体验| 人人妻人人添人人爽欧美一区卜| 久久鲁丝午夜福利片| 嫩草影院入口| 欧美精品一区二区大全| 久久99热6这里只有精品| 国产极品天堂在线| 国产精品熟女久久久久浪| 夜夜看夜夜爽夜夜摸| 精品一区二区三区视频在线| 国产午夜精品久久久久久一区二区三区| 91精品国产九色| 国产成人91sexporn| 久久狼人影院| 天堂8中文在线网| 精品酒店卫生间| 亚洲va在线va天堂va国产| 久久久久久久久久久久大奶| 日本黄色片子视频| 色婷婷av一区二区三区视频| 成人午夜精彩视频在线观看| 国产av国产精品国产| 熟女电影av网| 久久精品国产亚洲网站| 精品国产一区二区久久| 亚洲第一av免费看| 久久人人爽av亚洲精品天堂| 欧美区成人在线视频| 熟女人妻精品中文字幕| 国产一区二区三区av在线| 久久精品国产鲁丝片午夜精品| 少妇的逼好多水| 99国产精品免费福利视频| 男男h啪啪无遮挡| 久久99一区二区三区| 丰满饥渴人妻一区二区三| 亚洲国产av新网站| a级一级毛片免费在线观看| 国产黄片视频在线免费观看| 精品国产露脸久久av麻豆| 亚洲欧美一区二区三区国产| 精华霜和精华液先用哪个| 一个人免费看片子| 欧美性感艳星| 欧美日韩综合久久久久久| 国产淫片久久久久久久久| 欧美日韩在线观看h| 亚洲图色成人| 国产精品免费大片| 国产在线免费精品| 久久久久精品性色| 欧美bdsm另类| 国产日韩欧美在线精品| 人人妻人人添人人爽欧美一区卜| 肉色欧美久久久久久久蜜桃| 久久人人爽人人片av| 搡女人真爽免费视频火全软件| 免费少妇av软件| 免费大片18禁| 精品国产一区二区三区久久久樱花| 久久精品夜色国产| 最新中文字幕久久久久| 秋霞伦理黄片| 亚洲色图综合在线观看| 两个人的视频大全免费| 久久久久久久大尺度免费视频| 亚洲一级一片aⅴ在线观看| 80岁老熟妇乱子伦牲交| 毛片一级片免费看久久久久| 久久久亚洲精品成人影院| 欧美成人精品欧美一级黄| 97超视频在线观看视频| 中文资源天堂在线| 夜夜骑夜夜射夜夜干| 亚洲一区二区三区欧美精品| 精品99又大又爽又粗少妇毛片| 国产精品人妻久久久影院| 中国美白少妇内射xxxbb| 日韩欧美精品免费久久| 自拍偷自拍亚洲精品老妇| 亚洲国产精品国产精品| 欧美另类一区| 亚洲中文av在线| 亚洲精品日本国产第一区| 五月伊人婷婷丁香| 久久婷婷青草| 亚洲成人手机| 久久免费观看电影| 日日啪夜夜爽| 狂野欧美激情性xxxx在线观看| 国产白丝娇喘喷水9色精品| av女优亚洲男人天堂| 日韩中字成人| 青春草视频在线免费观看| 亚洲欧美日韩卡通动漫| 国内少妇人妻偷人精品xxx网站| 成年av动漫网址| 国国产精品蜜臀av免费| 夫妻性生交免费视频一级片| 97精品久久久久久久久久精品| 日韩精品有码人妻一区| 少妇的逼水好多| av国产久精品久网站免费入址| 日韩欧美 国产精品| 欧美日韩亚洲高清精品| 久久人人爽人人爽人人片va| 精品久久国产蜜桃| 五月开心婷婷网| 久久久久网色| 久久国产乱子免费精品| 欧美一级a爱片免费观看看| 国产av国产精品国产| 国内少妇人妻偷人精品xxx网站| 亚洲美女搞黄在线观看| 久久久久人妻精品一区果冻| 国产高清不卡午夜福利| 国产69精品久久久久777片| 国产精品国产三级国产av玫瑰| 亚洲无线观看免费| 老司机亚洲免费影院| 日韩伦理黄色片| 国产午夜精品一二区理论片| 插阴视频在线观看视频| 中文欧美无线码| 精品一区二区三卡| 国产精品久久久久久av不卡| 亚洲精品乱码久久久久久按摩| 欧美成人午夜免费资源| 亚洲av电影在线观看一区二区三区| 免费观看av网站的网址| 免费看日本二区| √禁漫天堂资源中文www| h日本视频在线播放| 最近2019中文字幕mv第一页| 日日撸夜夜添| 国产成人一区二区在线| 一区二区三区乱码不卡18| freevideosex欧美| 午夜影院在线不卡| 少妇猛男粗大的猛烈进出视频| 久久久午夜欧美精品| 99热这里只有是精品50| 久久 成人 亚洲| 一级,二级,三级黄色视频| 在线观看人妻少妇| 啦啦啦中文免费视频观看日本| 日本黄大片高清| a级毛色黄片| 亚洲一级一片aⅴ在线观看| 国产成人午夜福利电影在线观看| 人妻人人澡人人爽人人| 婷婷色综合大香蕉| 国产欧美日韩精品一区二区| 亚洲成人手机| 国产午夜精品久久久久久一区二区三区| 精品视频人人做人人爽| 啦啦啦视频在线资源免费观看| 最新的欧美精品一区二区| av在线老鸭窝| 少妇猛男粗大的猛烈进出视频| 一区二区三区精品91| 黄色欧美视频在线观看| 久久韩国三级中文字幕| 免费高清在线观看视频在线观看| 欧美97在线视频| 午夜福利视频精品|