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

    An Analysis of Mean Wave Drift Forces of FPSO Systems by Higher-order Boundary Element Method

    2014-01-19 05:49:42WANGKeZHANGZhiqiangXUWang
    船舶力學 2014年9期
    關鍵詞:張志強大連理工大學工程力學

    WANG Ke,ZHANG Zhi-qiang,XU Wang

    (State Key Laboratory of Structural Analysis for Industrial Equipment,Department of Engineering Mechanics,Dalian University of Technology,Dalian 116023,China)

    An Analysis of Mean Wave Drift Forces of FPSO Systems by Higher-order Boundary Element Method

    WANG Ke,ZHANG Zhi-qiang,XU Wang

    (State Key Laboratory of Structural Analysis for Industrial Equipment,Department of Engineering Mechanics,Dalian University of Technology,Dalian 116023,China)

    Mean wave drift forces and moments are critical for FPSO systems in mooring,riser installation and offloading.In this paper,higher-order boundary element method(HOBEM)is applied and symmetry is considered to calculate the second order hydrodynamic interactions of Schiehallion type FPSO in NorthSea of England.The surge,sway and yaw mean drift forces on FPSO are computed by near-field method based on direct pressure integration.Comparing to conventional method,current algorithm has high accuracy and can converge steadily even for high frequency waves.It is shown that in order to control FPSO safely under severe sea condition,the FPSO systems should change its weathervaning angle to head sea all the time.

    FPSO systems; mean wave drift forces and moments; higher-order boundary element method;near field method

    1 Introduction

    FPSO Systems(Floating Production,Storage and Offloading System)have become a very important equipment of the exploration and extraction of hydrocarbons from the sea-bottom,and are considered the only feasible solution for safe and cost efficient oil and gas production at large depth,especially in areas with severe wave environments.FPSO systems normally are located hundred miles from onshore zone and can not return to harbor like traditional transfer barge.This kind of production units is very prone to the effects of the weather,in particular as a result of their requirement to remain in the same position for several months or permanent in specific area.The designing of FPSO systems should consider two kinds of environmental loads.One is hydrodynamic forces and motions and the other is mean wave drift forces and moments.Hydrodynamic forces and motions are important to assure FPSO systems safely used in severe conditions.The FPSO systems are like a weathervane ship which is turret free to rotate around a vertical riser systems,this head wave working mechanism allows the ship to be aligned with the resultant of the environmental action,minimizing the loads on the mooringlines and also reducing rolling angle of the ship.On the other hand,as one kind of floating platforms,this ship shaped oil production systems will undergo slow motions from their mean position.The slow motions are excited by the slowly varying steady forces from the environment due to the weak restoring forces of the mooring system and the highly tuned low-frequency resonance of the system.It is therefore vital for the design of moorings and riser systems that the magnitudes of the motions are predicted well and that they are kept within acceptable bounds.The slowly varying steady forces are called mean wave drift forces and can be computed from quadratic products of the first-order quantities.It is well known that two different methods are available for computing the wave drift forces:the near-field method based on the direct pressure integration and the far-field method based on the momentumconservation principle.The near-field method(e.g.Ogilvie,1983)gives six degrees of forces but computations are rather complicated because various components must be evaluated.The far-field method(Maruo,1960;Newman,1967)is relatively easy to compute but gives only the total force on ship in horizontal plane with three degree of motions:surge,sway and yaw.

    FPSO systems is a very important naval architecture and should be safely maintained at any severe environmental condition,the wave loads prediction of this kind of platform must be as accurate as possible.Based on linear radiation and diffraction theory,hydrodynamic responses of FPSO can be calculated and used to engineering application during initial stage investigation(Liu Yingzhong et al,1986;Faltinsen,1990).In order to improve the numerical precision of second order mean wave drift force,in present paper the first order radiation and diffraction problem are solved by higher-order boundary element method and the symmetry of body geometry is considered,the mean wave drift forces are obtained by near field method for full six individual components.Basic theory and formulations for mean wave drift forces are expressed in section 2;section 3 discusses the result of mean wave drift forces about Schiehallion type FPSO in NorthSea of England.The conclusion of advantage of higher-order boundary element and understanding of mean wave drift forces are drawn in last section.

    2 Theory of boundary value problem

    The Cartesian coordinate system (x,y, )z defined in Fig.1 is stationary relative to the undisturbed position of the free surface and body.The origin is located on the free surface with z axis positive upward.The assumption of a potential flow permits the definition of the flow velocity as the gradient of the velocity potential Φ satisfying the Laplace equation

    in the fluid domain.The harmonic time dependence allows the definition of a complex velocity potential φ,related to Φ by Fig.1 FPSO sketch

    where Re denotes the real part,ω is the frequency of the incident wave and t is time.The ensuing boundary-value problem will be expressed in terms of the complex velocity potential φ,with the understanding that the product of all complex quantities with the factor e-iωtapplies.The linearized form of the free-surface condition is

    where K=ω2/g is deep water wave number,and g is the acceleration of gravity.The velocity potential of the incident wave is defined by

    The linearization of the problem permits the decomposition of the velocity potential φ into the radiation (φR)and diffraction (φD)components;while diffraction potential (φD)includes incident potential (φ0)and scattering potential (φ7),the velocity potential φ can thus be expressed in following:

    The constant ξjdenotes the complex amplitudes of the body oscillatory motion in its six rigid-body degrees of freedom,and φjdenotes the corresponding unit-amplitude radiation potentials.The velocity potential φ7represents the scattered disturbance of the incident wave by the body fixed at its undisturbed position.On the undisturbed position of the body boundary,the radiation and diffraction potentials are subject to the conditions:

    3 Boundary integral equation

    The boundary value problem(Eq.1 to Eq.7)of radiation and scattering potential above can be solved in terms of forming integral equation on body surface through Green second theorem.The radiation potential should satisfy following equation on body surface:

    here r=(x,y,z )denotes field point,r0= (ξ,η,ζ )denotes source points,Sbdenotes body wetted surface in calm water.From the computational point of view,equation(9)has some advantages over equation(8)in terms of CPU time and the requirement of storage space,because of the relative simplicity of the right-hand side.The Green function G( r; r0)is referred to as the wave source potential.It is the velocity potential at the point r due to a point source of strength-4π located at the point r0.It satisfies the free-surface and radiation conditions,and in finite water depth is defined by

    here J0()is the Bessel function of zero order.In equation(10)the Fourier k-integration is indented above the pole on the real axis in order to enforce the radiation condition.The calculation of Green function and wave body radiation and diffraction are referred to Newman(1985,1986,1988).

    4 Formulation discretization

    The integral equations are approximated by discretizing the surfaces of integration by boundary elements.The surfaces cover the body and interior water plane area.The values inside an element are represented in terms of nodal values and shape functions.On the discretized body surface and on the inner water plane area,both triangular and quadrilateral elements are used.By introducing shape function hi(ξ, )η in local intrinsic coordinates in an element,a variable and its derivatives can be represented by the values at the nodes as follows:

    where s is the number of nodes in the element.Using the shape functions above and the transformation relationship between local intrinsic coordinates and the global coordinates,we can discretise the integral equation for the first order problem as:

    here J( ξ,η )is the Jacobian for the coordinate transformation.NBand NF′are the numbers of elements on the body surface and on the interior water plane area respectively.

    5 Formulations of mean wave drift force

    The mean wave drift force is crucial to mooring system of FPSO.Normally there are two methods to calculate mean wave drift force:far field method and near field method(Kim,Yue,1990;Kashiwagi,Endob,Yamaguchi,2005).Far field method can obtain three components in horizontal plane:surge,sway and yaw;while near field method can obtain the whole six components:surge,sway,heave,roll,pitch,yaw.Far field method can get surge and sway mean wave drift forces with very high precision but yaw moment converges slowly.Yaw mean drift force can be approximated by two dimensional slip-method with fairy good result for slender body,but to FPSO and TLP type platform with complicated body geometry,yaw moment is not easily converged in the range of high wave frequency.The other disadvantage of far field method is that it can only obtain a final total force whether the single or multiple bodies are considered.Near field method calculates wave force through pressure integration upon the whole body surface,it can get all six components of wave force with converged results with enough element number and higher-order method.The drawback of near field method is that it needs very complicated mathematical manipulation.Here we use near field method and the formu-

    lations to calculate the mean wave drift forcesand momentsare listed as follows:

    6 Analysis of mean wave drift force of FPSO

    From viewpoint of engineering application,the horizontal plane mean wave drift force surge,sway and yaw are thought to be the crucial criterion for designing of mooring systems of FPSO.This section discusses mean wave drift force of one kind of Schiehallion type FPSO currently working in NorthSea of England.Considering regular incident wave with unit amplitude,the wave direction is listed in Fig.1.The longitude of FPSO is 238m,width is 45m,draft is 11.5m,water depth is 240m,the gravity center of vehicle is 5.96 under water line.Fig.2 is the mesh of FPSO,since symmetry method is utilized only one plane mesh plot is given.In each symmetry plane the total element number is 420 and the number of control node is 915.

    Fig.3 Surge drift force Fig.4 Sway drift force

    Figs.3-5 show mean wave drift forces and moments surge,sway and yaw of FPSO under different incident waves.Wave drift force means continuous average wave resistance acting on body.In these figures,x axis means shallow water wave frequency,y axis means wave force or moment,they are both expressed with dimensional form.The surge and sway drift forces are both minus under current wave direction,but in order to make the description easily we discuss drift force based on formulation FX=-X,FY=-Y.In each wave direction,the frequency interval is 0.05 and totally 21 frequencies in the range between 0 to 1.05 are calculated.

    From Figs.3-5,it can find that when ω<0.4,mean wave drift forces are very small and FPSO is always in steady motion under any wave direction;generally for ω>0.4,surge,sway and yaw drift force or moment will increase with frequency but converge to different values for ω>1.0.Around frequency ω>1.0,the sway drift force will decreasingly converge to different constants with increasing of wave directions.Under this situation for ω>0.4,if FPSO is not efficiently constrained,it will move slowly from its initial position.While in the range 0.4<ω<0.8,drift forces and moment vary rapidly and 0.05 frequency interval is still large for the calculation of second order wave loads.Although the wave directions are different,the maximum values of surge,sway and yaw drift forces or moment are occurring in frequency ranges 0.55<ω<0.7 which corresponding to wave period 8.97s<T<11.4s.

    From Fig.3,in the frequency range 0.4<ω<0.8,the wave directions causing larger surge force are around 120°≤β≤150°,even though the maximum value of surge does not exist in this field but the whole surge force within this range varies slightly and from the viewpoint of ship manipulation,it is easy to control its motions;for ω>0.8,the surge force will decreasingly converge to constant values under wave directions 105°≤β≤180°,the force value is between 8.0×104N and 1.0×105N.To beam wave direction β=90°,the surge force is very small but not equal to zero,this is mainly because the geometry of FPSO is symmetric to x axis not y axis.This finding reveals the precision of model meshing and application of higher order boundary element method guarantees the drift forces converge efficiently at high frequency when ω>0.8.

    To sway drift force,the varying tendency is similar to surge drift force,the difference is that when ω>0.8 for different wave directions,the sway drift force will increase with frequency and converge to different values around frequency ω>1.0.The final steady values of sway drift force increase from head wave to beam wave and the beam wave β=90°will produce maximal force.Comparing with Fig.3 and Fig.4,the maximal surge drift force in beam wave is around 1.0×104N while maximal sway drift force is around 1.0×106N.The latter is nearly 100 times of the former.In head wave direction β=180°,the maximal surge drift force is around 2.2×105N while sway drift force is zero because of x axis symmetry of FPSO;at other wave directions,the relations of maximal sway drift force and surge drift force are:10 times(β=105°);4 times(β=120°);3 times(β=135°);2 times(β=150°);0.77 times(β=165°).From the respect of real working conditions of FPSO,the maximal surge drift force is around 2.2×105N(β=180°)and maximal sway drift force is around 1.0×106N (β=90°),the sway drift force is 5 times of surge drift force.

    Comparing to Fig.3 and Fig.4,Fig.5 is the result of yaw drift moment and has complicated varying curves.Generally it can find that at frequency ω<0.4,MZ=0 exists under any wave directions which means FPSO at this time is in absolutely safe motion condition.But it should be pointed out that at frequency ω=0.3,an irregular frequency occurs which will cause resonance of body motion at low wave frequency.At frequency ω>0.8,only head wave can be thought to be in steady condition since at this case the yaw drift moment is zero because of symmetry of FPSO.At beam wave the yaw drift moment is small but can not be ignored,the yaw drift moment at other wave directions are large.It also can be found that in high wave frequency,the surge and sway drift will converge to constant values,but yaw drift moment still varies with the changing of wave direction from head wave to beam wave,its tendency is following the changing of formulation ksin β×kcos β [Newman,1974],where k is wave number.

    6 Conclusions

    The precision of mean wave drift force of FPSO can be improved by higher-order boundary element method and considering geometrical symmetry.The surge,sway and yaw drift force are three important factors in design of mooring systems of FPSO.The most dangerous mean wave drift force occurs when wave incident direction is in the range 90°≤β≤135°,while wave frequency is in the range 0.4<ω<0.8.In order to control FPSO safely under severe sea condition,the FPSO systems should change its weathervaning angle to head sea all the time,since sway and yaw drift force are both zero under that time.

    [1]Liu Yingzhong,Miao Guoping.Ship motion in waves[M].Shanghai:University of Shanghai Jiaotong Press,1986.

    [2]Faltinsen O M.Sea loads on ships and offshore structures[M].Cambridge University Press,Cambridge,UK,1990.

    [3]Ogilvie T F.Second order hydrodynamic effects on ocean platforms[C]//Inernational Workshop on Ship and Platform Motion.Berkeley,USA,1983:205-265.

    [4]Maruo H.The drift of a body floating on waves[J].Journal of Ship Research,1960,4:1-10.

    [5]Newman J N.The drift force and moment on ships in waves[J].Journal of Ship Research,1967,11(1):51-60.

    [6]Newman J N.Second-order slowly varying forces on vessels in irregular waves[C]//Proc Int Symp Dynm Marine Vehicles Struct Waves.IMechE,London,1974.

    [7]Newman J N,Sclavounos P D.The computation of wave loads on large offshore structures[C]//BOSS’88 Conference.Trondheim,Norway,1988.

    [8]Newman J N.Algorithms for the free-surface green function[J].Journal of Engineering Mathematics,1985,19:57-67.

    [9]Newman J N.Distribution of sources and dipoles over a quadrilateral[J].Journal of Engineering Mathematics,1986,20:113-126.

    [10]Kim M H,Yue D K P.The complete second-order diffraction solution for an axis-symmetric body.Part 2:Bi-chromatic incident waves and body motions[J].Journal of Fluid Mechanics,1990,211:557-593.

    [11]Kashiwagi M,Endob K,Yamaguchi H.Wave drift forces and moments on two ships arranged side by side in waves[J].Ocean Engineering,2005,32:529-555.

    應用高階邊界元法求解FPSO型采油平臺波浪漂移力的研究

    王 科,張志強,許 旺

    (大連理工大學工業(yè)裝備結構分析國家重點實驗室,工程力學系,遼寧 大連116024)

    波浪平均漂移力是影響FPSO系統(tǒng)系泊、立管安裝和卸載的重要外在荷載。文中針對英國北海正投入使用的薩哈林型式FPSO,應用高階邊界單元方法對其二階水動力特性進行了計算分析,計算中考慮了物體的幾何對稱性,F(xiàn)PSO系統(tǒng)所受到的縱蕩,橫蕩以及首搖方向的波浪平均漂移力采用近場解法在濕表面上直接進行壓力積分獲得。與常數(shù)單元方法相比,該方法具有很高的精度,特別是對于高頻波浪能得到可靠的收斂解。研究表明為了使得FPSO系統(tǒng)在惡劣海況下安全使用,F(xiàn)PSO系統(tǒng)必須始終保持迎浪狀態(tài)工作。

    FPSO;波浪平均漂移力;高階邊界元方法;近場解法

    U661.1 U674.38+1

    A

    王 科(1970-),男,大連理工大學工程力學系講師;

    許 旺(1983-),男,大連理工大學碩士生。

    U661.1 U674.38+1

    A

    1007-7294(2010)06-0587-09

    date:2009-08-21

    Supported by the Scientific Research Foundation for the Returned Overseas Chinese Scholars,State Education Ministry,2007[24]

    Biography:WANG Ke(1970-),male,lecturer of Dalian University of Technology,E-mail:kwang@dlut.edu.cn,Tel:0411-84707330.

    張志強(1984-),男,大連理工大學博士生;

    猜你喜歡
    張志強大連理工大學工程力學
    Origin of tradeoff between movement velocity and attachment duration of kinesin motor on a microtubule
    《工程力學》微課教學模式探索
    Research on the Globalization of English in the Internet era
    大東方(2019年1期)2019-09-10 20:30:40
    這才叫大藏家——張志強先生和他的精絕藏石
    寶藏(2018年12期)2019-01-29 01:50:50
    偽隨機碼掩蔽的擴頻信息隱藏
    工程力學
    志強的石
    中華奇石(2014年12期)2014-07-09 18:30:22
    虛擬仿真技術在復雜工程力學問題實驗教學中的應用
    中泰化學與大連理工大學簽署戰(zhàn)略合作框架協(xié)議
    中國氯堿(2014年11期)2014-02-28 01:05:06
    大學物理和工程力學在機構運動分析中的應用
    物理與工程(2014年5期)2014-02-27 11:23:15
    国产91av在线免费观看| 亚洲av二区三区四区| 国国产精品蜜臀av免费| 中文精品一卡2卡3卡4更新| 极品人妻少妇av视频| 免费播放大片免费观看视频在线观看| 少妇丰满av| 色婷婷av一区二区三区视频| 我的女老师完整版在线观看| 我的女老师完整版在线观看| 一区二区三区精品91| 久久99一区二区三区| 日本黄大片高清| 国产在线免费精品| 蜜桃在线观看..| 亚洲在久久综合| 日日摸夜夜添夜夜添av毛片| 亚洲丝袜综合中文字幕| 97超碰精品成人国产| 国产成人freesex在线| 新久久久久国产一级毛片| 久久久久久久久久久丰满| 久久精品国产a三级三级三级| 久久精品国产a三级三级三级| 免费少妇av软件| 两个人免费观看高清视频 | 我的老师免费观看完整版| 成年美女黄网站色视频大全免费 | 自线自在国产av| 国产精品嫩草影院av在线观看| 国产免费视频播放在线视频| 亚洲婷婷狠狠爱综合网| 99九九在线精品视频 | 岛国毛片在线播放| 国产精品一区二区性色av| 亚洲内射少妇av| 日本av免费视频播放| 视频区图区小说| 大片电影免费在线观看免费| 久久久久人妻精品一区果冻| 97超视频在线观看视频| 国国产精品蜜臀av免费| av在线播放精品| 高清欧美精品videossex| 女性被躁到高潮视频| 国产精品嫩草影院av在线观看| 男人和女人高潮做爰伦理| 欧美精品亚洲一区二区| 国产无遮挡羞羞视频在线观看| 免费在线观看成人毛片| 狂野欧美白嫩少妇大欣赏| 韩国av在线不卡| 欧美成人精品欧美一级黄| 老熟女久久久| 人人妻人人看人人澡| 看免费成人av毛片| 国产熟女欧美一区二区| 高清黄色对白视频在线免费看 | 久久久久久久精品精品| 18禁在线播放成人免费| 婷婷色综合www| 男人和女人高潮做爰伦理| 国产精品99久久99久久久不卡 | 91aial.com中文字幕在线观看| 日韩伦理黄色片| 精品午夜福利在线看| 久久人人爽人人爽人人片va| 99九九线精品视频在线观看视频| 一级毛片久久久久久久久女| 99久国产av精品国产电影| 18禁在线无遮挡免费观看视频| 大片免费播放器 马上看| 18禁裸乳无遮挡动漫免费视频| 精品久久国产蜜桃| 国产视频首页在线观看| 国产精品三级大全| xxx大片免费视频| 18禁裸乳无遮挡动漫免费视频| 亚洲国产精品一区三区| 久久久久久久精品精品| 啦啦啦视频在线资源免费观看| av在线播放精品| 国产亚洲精品久久久com| 91久久精品国产一区二区成人| 最后的刺客免费高清国语| 在现免费观看毛片| 最近的中文字幕免费完整| 99热这里只有是精品在线观看| 熟女人妻精品中文字幕| 丝袜脚勾引网站| 亚洲国产色片| 少妇 在线观看| 欧美成人精品欧美一级黄| 观看免费一级毛片| 国产日韩欧美在线精品| 欧美激情国产日韩精品一区| 精品熟女少妇av免费看| 边亲边吃奶的免费视频| 黄片无遮挡物在线观看| 日韩不卡一区二区三区视频在线| 亚洲av男天堂| 国产一区亚洲一区在线观看| 国产精品国产三级国产av玫瑰| 国产欧美另类精品又又久久亚洲欧美| 色哟哟·www| 久久人人爽av亚洲精品天堂| 精品少妇久久久久久888优播| 日韩成人av中文字幕在线观看| 人人妻人人添人人爽欧美一区卜| 国产成人a∨麻豆精品| 18禁在线无遮挡免费观看视频| 国产亚洲精品久久久com| 欧美丝袜亚洲另类| 高清在线视频一区二区三区| 日韩视频在线欧美| 免费大片黄手机在线观看| 欧美日韩一区二区视频在线观看视频在线| 亚洲精品亚洲一区二区| 搡女人真爽免费视频火全软件| 男女边吃奶边做爰视频| 一级毛片久久久久久久久女| 少妇猛男粗大的猛烈进出视频| 国内精品宾馆在线| 妹子高潮喷水视频| 新久久久久国产一级毛片| 国产一区有黄有色的免费视频| 久久久久国产网址| 全区人妻精品视频| 亚洲经典国产精华液单| 女人久久www免费人成看片| 99re6热这里在线精品视频| 性高湖久久久久久久久免费观看| 91久久精品电影网| 汤姆久久久久久久影院中文字幕| 五月伊人婷婷丁香| 噜噜噜噜噜久久久久久91| 久久 成人 亚洲| 国产精品三级大全| av国产久精品久网站免费入址| 亚洲精品成人av观看孕妇| 久久久久久人妻| 蜜桃在线观看..| 午夜福利在线观看免费完整高清在| 亚洲精品亚洲一区二区| 国产欧美另类精品又又久久亚洲欧美| 日日摸夜夜添夜夜爱| 乱系列少妇在线播放| 在线亚洲精品国产二区图片欧美 | 亚洲av免费高清在线观看| 国产精品伦人一区二区| 久久精品国产亚洲av天美| av视频免费观看在线观看| 女人久久www免费人成看片| 亚洲天堂av无毛| 国产亚洲最大av| 国产黄片视频在线免费观看| 99久国产av精品国产电影| 99热这里只有是精品在线观看| 中文字幕人妻熟人妻熟丝袜美| 久久精品熟女亚洲av麻豆精品| 国产黄片美女视频| 欧美bdsm另类| 中国国产av一级| 日韩 亚洲 欧美在线| 国产一级毛片在线| 欧美+日韩+精品| av网站免费在线观看视频| 深夜a级毛片| 色哟哟·www| 老司机影院成人| 丝瓜视频免费看黄片| 国产在线视频一区二区| 亚洲欧美精品专区久久| 一本色道久久久久久精品综合| 亚洲欧洲精品一区二区精品久久久 | 最近中文字幕高清免费大全6| 青春草国产在线视频| 一级毛片aaaaaa免费看小| 精品国产露脸久久av麻豆| 亚洲av成人精品一二三区| 久久久久人妻精品一区果冻| 国产免费一区二区三区四区乱码| 性色av一级| 久久精品国产鲁丝片午夜精品| 亚洲在久久综合| 亚洲精品亚洲一区二区| 亚洲一级一片aⅴ在线观看| 黄色毛片三级朝国网站 | 午夜日本视频在线| 免费不卡的大黄色大毛片视频在线观看| 卡戴珊不雅视频在线播放| 狂野欧美白嫩少妇大欣赏| 日本欧美视频一区| 国产精品一区二区在线观看99| 成人国产麻豆网| 一区二区三区四区激情视频| 午夜福利影视在线免费观看| 91精品国产九色| 成年女人在线观看亚洲视频| 人体艺术视频欧美日本| 亚洲一区二区三区欧美精品| 午夜久久久在线观看| 性色avwww在线观看| 欧美 日韩 精品 国产| 日韩在线高清观看一区二区三区| 久久国产精品大桥未久av | 日韩亚洲欧美综合| 大香蕉97超碰在线| 大又大粗又爽又黄少妇毛片口| 亚洲国产精品一区二区三区在线| 亚洲精品乱久久久久久| 亚洲国产最新在线播放| 亚洲av不卡在线观看| 国产在线视频一区二区| 日韩伦理黄色片| 一二三四中文在线观看免费高清| 啦啦啦啦在线视频资源| 国产精品99久久久久久久久| 日本av免费视频播放| 9色porny在线观看| av不卡在线播放| 国产乱人偷精品视频| 人人妻人人澡人人爽人人夜夜| 交换朋友夫妻互换小说| 一边亲一边摸免费视频| 最后的刺客免费高清国语| 欧美精品高潮呻吟av久久| 国产熟女欧美一区二区| 国内揄拍国产精品人妻在线| 久久久久久久精品精品| 国产精品久久久久久久久免| 国产伦精品一区二区三区四那| 久久鲁丝午夜福利片| 日韩av在线免费看完整版不卡| 中国三级夫妇交换| 亚洲国产精品专区欧美| 国产av国产精品国产| 欧美成人精品欧美一级黄| 看免费成人av毛片| 亚洲精品日韩av片在线观看| 亚洲欧洲精品一区二区精品久久久 | 亚洲欧美一区二区三区黑人 | 国产成人a∨麻豆精品| 免费看不卡的av| 欧美老熟妇乱子伦牲交| 日本av免费视频播放| 人人妻人人添人人爽欧美一区卜| 日韩一区二区三区影片| 99国产精品免费福利视频| 五月伊人婷婷丁香| 国模一区二区三区四区视频| 国产日韩欧美亚洲二区| 色网站视频免费| 亚洲综合色惰| av国产久精品久网站免费入址| 午夜精品国产一区二区电影| 久久久久久久久久成人| 69精品国产乱码久久久| av卡一久久| 亚洲不卡免费看| 免费av不卡在线播放| 99热网站在线观看| 国产av国产精品国产| 热re99久久精品国产66热6| 久久婷婷青草| 国产精品成人在线| 在线观看美女被高潮喷水网站| 免费看av在线观看网站| 人人妻人人看人人澡| 亚洲精品日韩av片在线观看| 国精品久久久久久国模美| av在线app专区| 三级国产精品欧美在线观看| 777米奇影视久久| 中文乱码字字幕精品一区二区三区| 在线观看一区二区三区激情| 少妇高潮的动态图| 精品久久久久久久久亚洲| 亚洲国产精品一区三区| 天堂俺去俺来也www色官网| 夜夜爽夜夜爽视频| 国产精品一区二区性色av| 各种免费的搞黄视频| 日日爽夜夜爽网站| 内地一区二区视频在线| 日本免费在线观看一区| 有码 亚洲区| 麻豆精品久久久久久蜜桃| 人妻制服诱惑在线中文字幕| 美女国产视频在线观看| 久久久久网色| 中文资源天堂在线| 午夜av观看不卡| 一二三四中文在线观看免费高清| 丁香六月天网| 天天躁夜夜躁狠狠久久av| 观看美女的网站| 国产精品蜜桃在线观看| 五月开心婷婷网| 亚洲精品乱久久久久久| 精品人妻熟女av久视频| 国产伦在线观看视频一区| 99热国产这里只有精品6| 欧美日韩国产mv在线观看视频| 精品一区二区三卡| av在线观看视频网站免费| 欧美另类一区| 秋霞在线观看毛片| 十八禁网站网址无遮挡 | 精品人妻熟女av久视频| 色5月婷婷丁香| 精品久久久久久久久亚洲| 国语对白做爰xxxⅹ性视频网站| 精品一区二区三卡| 街头女战士在线观看网站| 黄色毛片三级朝国网站 | 日本av免费视频播放| 国产免费一区二区三区四区乱码| 丁香六月天网| 国产黄色免费在线视频| 免费观看无遮挡的男女| 中文在线观看免费www的网站| 91久久精品电影网| 日本爱情动作片www.在线观看| 欧美成人精品欧美一级黄| 亚洲电影在线观看av| 亚洲欧洲精品一区二区精品久久久 | 久久久久久久久久成人| 国产精品秋霞免费鲁丝片| 午夜91福利影院| 午夜视频国产福利| 少妇人妻一区二区三区视频| 一级毛片我不卡| 在线播放无遮挡| 免费播放大片免费观看视频在线观看| 男人和女人高潮做爰伦理| 色5月婷婷丁香| 成人国产麻豆网| 在现免费观看毛片| 日韩成人伦理影院| 国产成人91sexporn| a 毛片基地| 少妇被粗大的猛进出69影院 | 日日摸夜夜添夜夜爱| 亚洲三级黄色毛片| 街头女战士在线观看网站| av国产久精品久网站免费入址| 亚洲人成网站在线观看播放| 成人黄色视频免费在线看| 亚洲国产最新在线播放| h视频一区二区三区| 69精品国产乱码久久久| 人妻人人澡人人爽人人| 国产午夜精品一二区理论片| 街头女战士在线观看网站| 亚洲国产最新在线播放| 日日摸夜夜添夜夜添av毛片| 18禁裸乳无遮挡动漫免费视频| 国产乱来视频区| 欧美三级亚洲精品| 免费观看性生交大片5| 青春草视频在线免费观看| 日本色播在线视频| 高清欧美精品videossex| 亚洲,欧美,日韩| 国产成人91sexporn| 免费在线观看成人毛片| 国产有黄有色有爽视频| 久久ye,这里只有精品| 亚洲精品乱码久久久v下载方式| 纵有疾风起免费观看全集完整版| 国产av国产精品国产| 精品一品国产午夜福利视频| 日韩伦理黄色片| 简卡轻食公司| 国产淫语在线视频| 一区在线观看完整版| 国产午夜精品久久久久久一区二区三区| 久久久精品94久久精品| 亚洲不卡免费看| 日韩欧美 国产精品| 偷拍熟女少妇极品色| 国产一区有黄有色的免费视频| 一区二区三区精品91| 欧美日韩国产mv在线观看视频| 熟女av电影| 亚州av有码| 中文字幕免费在线视频6| 久久人妻熟女aⅴ| 久久青草综合色| 国产在线男女| 国产精品国产三级国产av玫瑰| 久久99热6这里只有精品| 中文字幕人妻丝袜制服| 国产视频首页在线观看| 欧美精品高潮呻吟av久久| 一区二区三区精品91| 婷婷色综合www| 成年人免费黄色播放视频 | 免费观看的影片在线观看| 熟女av电影| 777米奇影视久久| 99热这里只有精品一区| 午夜日本视频在线| 一级毛片我不卡| 免费大片黄手机在线观看| 91久久精品国产一区二区三区| 在线观看人妻少妇| 亚洲天堂av无毛| 人妻少妇偷人精品九色| h日本视频在线播放| 青青草视频在线视频观看| 一级爰片在线观看| 婷婷色综合大香蕉| 欧美精品亚洲一区二区| 成人亚洲精品一区在线观看| 成人特级av手机在线观看| 精品亚洲成国产av| 哪个播放器可以免费观看大片| 免费观看的影片在线观看| 国产成人freesex在线| 国产精品免费大片| 中文字幕亚洲精品专区| 美女大奶头黄色视频| 欧美日韩国产mv在线观看视频| 国产探花极品一区二区| 97在线视频观看| 精品人妻一区二区三区麻豆| 这个男人来自地球电影免费观看 | 久久热精品热| 午夜福利视频精品| 伦精品一区二区三区| 在线观看三级黄色| 国产深夜福利视频在线观看| 男女边摸边吃奶| 春色校园在线视频观看| 欧美激情极品国产一区二区三区 | 久久鲁丝午夜福利片| 免费av中文字幕在线| 国产日韩欧美亚洲二区| 亚洲精品日韩av片在线观看| 狂野欧美激情性bbbbbb| 欧美xxⅹ黑人| av卡一久久| 18+在线观看网站| 亚洲国产欧美在线一区| 成年美女黄网站色视频大全免费 | 久久午夜综合久久蜜桃| 国产免费又黄又爽又色| av免费在线看不卡| 九九爱精品视频在线观看| 男人舔奶头视频| 午夜福利影视在线免费观看| 性高湖久久久久久久久免费观看| 9色porny在线观看| 日韩成人伦理影院| 我要看日韩黄色一级片| 国产色婷婷99| 精品一区二区免费观看| 蜜桃久久精品国产亚洲av| 国产精品无大码| 久久99热6这里只有精品| 男的添女的下面高潮视频| 2021少妇久久久久久久久久久| 国产在线免费精品| 成人漫画全彩无遮挡| 国产伦精品一区二区三区四那| 五月玫瑰六月丁香| 精品人妻一区二区三区麻豆| 日本av免费视频播放| 水蜜桃什么品种好| 亚洲第一区二区三区不卡| 久久午夜福利片| www.av在线官网国产| 午夜免费鲁丝| 亚洲国产精品999| 国产亚洲最大av| 男人添女人高潮全过程视频| av免费在线看不卡| 新久久久久国产一级毛片| 看非洲黑人一级黄片| 91精品国产九色| 久久99蜜桃精品久久| 亚洲欧美日韩另类电影网站| 国产永久视频网站| av不卡在线播放| 亚洲电影在线观看av| videossex国产| .国产精品久久| 午夜福利视频精品| 国产成人freesex在线| 色视频在线一区二区三区| 91精品一卡2卡3卡4卡| a级毛片免费高清观看在线播放| 观看免费一级毛片| 美女国产视频在线观看| 亚洲av在线观看美女高潮| 久久精品国产自在天天线| 亚洲av综合色区一区| 国产精品偷伦视频观看了| 国产日韩欧美视频二区| 你懂的网址亚洲精品在线观看| 卡戴珊不雅视频在线播放| 亚洲中文av在线| 国产免费视频播放在线视频| 精品酒店卫生间| 黄色毛片三级朝国网站 | 免费av中文字幕在线| 国产高清三级在线| 丰满人妻一区二区三区视频av| 在线观看人妻少妇| 少妇被粗大猛烈的视频| 如何舔出高潮| av福利片在线| 纯流量卡能插随身wifi吗| 99热国产这里只有精品6| 午夜精品国产一区二区电影| 国产极品粉嫩免费观看在线 | 国产精品久久久久久av不卡| 亚洲国产日韩一区二区| 久久久久久久久大av| 亚洲精品中文字幕在线视频 | 欧美日韩视频精品一区| 97精品久久久久久久久久精品| 欧美区成人在线视频| 91久久精品国产一区二区成人| www.av在线官网国产| 少妇猛男粗大的猛烈进出视频| 亚洲国产日韩一区二区| 日本爱情动作片www.在线观看| 人妻制服诱惑在线中文字幕| 国精品久久久久久国模美| 99re6热这里在线精品视频| 日韩成人av中文字幕在线观看| 亚洲久久久国产精品| 精品少妇久久久久久888优播| 99re6热这里在线精品视频| 黄色怎么调成土黄色| www.av在线官网国产| 亚洲成人手机| 一级毛片黄色毛片免费观看视频| 久久国内精品自在自线图片| 午夜福利在线观看免费完整高清在| 亚洲四区av| 成人毛片60女人毛片免费| 婷婷色综合www| 寂寞人妻少妇视频99o| 一本大道久久a久久精品| 午夜福利网站1000一区二区三区| 女性被躁到高潮视频| 91午夜精品亚洲一区二区三区| 亚洲精品第二区| 欧美亚洲 丝袜 人妻 在线| 国产亚洲一区二区精品| 日本午夜av视频| 亚洲不卡免费看| 97超视频在线观看视频| 久久久国产欧美日韩av| 亚洲自偷自拍三级| h视频一区二区三区| 亚洲av电影在线观看一区二区三区| 在线精品无人区一区二区三| 欧美日韩视频高清一区二区三区二| 麻豆成人av视频| 久久97久久精品| 国产白丝娇喘喷水9色精品| 好男人视频免费观看在线| 永久网站在线| 少妇 在线观看| 国产一区二区三区av在线| 成人二区视频| 又大又黄又爽视频免费| 老熟女久久久| 亚洲av福利一区| 国产又色又爽无遮挡免| 久久久久久久精品精品| 成人亚洲精品一区在线观看| 夜夜骑夜夜射夜夜干| 欧美精品一区二区大全| 建设人人有责人人尽责人人享有的| 少妇人妻一区二区三区视频| 内地一区二区视频在线| 一级毛片 在线播放| 久久这里有精品视频免费| 在线观看国产h片| 国产精品99久久久久久久久| 高清不卡的av网站| 在线天堂最新版资源| 亚洲欧美一区二区三区国产| 亚洲欧美精品专区久久| 欧美成人午夜免费资源| 久久久久久久久大av| 久久精品熟女亚洲av麻豆精品| 99九九线精品视频在线观看视频| 精品人妻熟女毛片av久久网站| 热re99久久国产66热| 欧美xxⅹ黑人| 国产日韩欧美视频二区| 26uuu在线亚洲综合色| 免费黄频网站在线观看国产| 97在线视频观看| 人妻制服诱惑在线中文字幕| 人妻人人澡人人爽人人| 边亲边吃奶的免费视频| 建设人人有责人人尽责人人享有的| 精华霜和精华液先用哪个| 中文在线观看免费www的网站| 日本猛色少妇xxxxx猛交久久| 欧美+日韩+精品| 六月丁香七月| 中文天堂在线官网| 国产一级毛片在线| 亚洲伊人久久精品综合| 精品久久久噜噜| 欧美日韩视频高清一区二区三区二|