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

    全冷式超大型液化氣船波浪載荷預(yù)報(bào)與分析

    2016-10-12 05:38:50周清華梁小軍李小靈
    中國(guó)航海 2016年3期
    關(guān)鍵詞:液艙船體波浪

    周清華, 梁小軍, 李小靈, 陳 兵

    (江南造船(集團(tuán))有限責(zé)任公司, 上海 201913)

    全冷式超大型液化氣船波浪載荷預(yù)報(bào)與分析

    周清華, 梁小軍, 李小靈, 陳 兵

    (江南造船(集團(tuán))有限責(zé)任公司, 上海 201913)

    全冷式超大型液化氣船(Very Large Gas Carrier,VLGC)的運(yùn)動(dòng)及波浪載荷預(yù)報(bào)對(duì)船體、圍護(hù)系統(tǒng)和液艙支座的結(jié)構(gòu)安全性評(píng)估至關(guān)重要。對(duì)此,以83 000 m3VLGC為例,采用WALCS程序系統(tǒng)作為主要分析工具,應(yīng)用三維線性勢(shì)流理論,在頻域內(nèi)預(yù)報(bào)該VLGC船體運(yùn)動(dòng)與波浪載荷的傳遞函數(shù)和長(zhǎng)期預(yù)報(bào)極值。以第二液貨艙重心處的加速度和舯剖面的垂向波浪彎矩為研究對(duì)象,比較分析WALCS預(yù)報(bào)值與IGC設(shè)計(jì)值、UR-S11規(guī)范值及SESAM預(yù)報(bào)值的差異。分析結(jié)果可為VLGC的結(jié)構(gòu)安全性評(píng)估提供載荷輸入和參考。

    全冷式超大型液化氣船(VLGC);加速度;垂向波浪彎矩;傳遞函數(shù);長(zhǎng)期預(yù)報(bào)

    Abstract: Ship motion and wave load prediction of Very Large Gas Carrier(VLGC) is essential for the safety assessment of hull structure, containment system and tank support. An 83 000 m3VLGC is taken as the object of study for illustration. The hydrodynamic analysis software WALCS based on three dimensional linear potential flow theory is used as the primary analytical tool to predict the hydrodynamic characteristics of ship motion and wave load in frequency domain. The gravity acceleration of second cargo tank and the vertical wave bending moment of midship section are calculated. The WALCS forecast values are compared with those from IGC rule,UR-S11 standard and SESAM. The calculated results can be a reference of the excitation input for the structural safety assessment of VLGC.

    Keywords: VLGC; acceleration; vertical wave bending moment; transfer function; long term prediction

    83 000 m3全冷式超大型液化氣船(Very Large Gas Carrier,VLGC)[1]是江南造船(集團(tuán))有限責(zé)任公司研制的高技術(shù)、高附加值船舶,打破了國(guó)外船廠的技術(shù)壟斷。相比常規(guī)船型,VLGC對(duì)船體波浪載荷預(yù)報(bào)的要求更高,主要原因有:

    (1) 船長(zhǎng)超過(guò)200 m,質(zhì)量分布特殊,規(guī)范計(jì)算法不足以反映其在惡劣海況下的載荷響應(yīng)特征;

    (2) 液貨加速度產(chǎn)生內(nèi)部液體壓力,應(yīng)考慮船體運(yùn)動(dòng)的長(zhǎng)期預(yù)報(bào);

    (3) 圍護(hù)系統(tǒng)的安全性要求高,該船設(shè)置有4個(gè)A型獨(dú)立棱形液艙,并通過(guò)特殊的支座支撐在貨艙內(nèi),承受液艙和貨物的質(zhì)量及各類(lèi)靜、動(dòng)載荷,由于支座只承受壓力而不傳遞拉力,因此在進(jìn)行VLGC結(jié)構(gòu)分析時(shí)需進(jìn)行接觸分析或大量迭代運(yùn)算。[2]

    由此可見(jiàn),為保證船體、圍護(hù)系統(tǒng)和液艙支座的結(jié)構(gòu)安全性,準(zhǔn)確計(jì)算出船體運(yùn)動(dòng)加速度和波浪載荷尤為重要。目前船舶在波浪中的運(yùn)動(dòng)和載荷響應(yīng)預(yù)報(bào)方法[3]已日臻成熟,主要有切片理論、三維理論、非線性時(shí)域理論及考慮彈振的水彈性理論等。這些方法中,基于勢(shì)流理論的三維頻域線性計(jì)算方法[4]在工程中應(yīng)用較為廣泛。船級(jí)社和相關(guān)科研院校已開(kāi)發(fā)出相應(yīng)的軟件,如DNV/SESAM,BV/HYDROSTAR和MIT/WAMIT等。這里選取WALCS程序系統(tǒng)對(duì)83 000 m3VLGC的運(yùn)動(dòng)及波浪載荷進(jìn)行預(yù)報(bào)與分析。該船的技術(shù)參數(shù)為:總長(zhǎng)226 m;垂線間長(zhǎng)215 m;型寬36.6 m;型深22.2 m;設(shè)計(jì)吃水11.4 m;載重量54 000 t;設(shè)計(jì)航速16.8 kn。為驗(yàn)證WALCS數(shù)值計(jì)算的有效性,以加速度和垂向波浪彎矩為研究對(duì)象,比較分析WALCS預(yù)報(bào)值與IGC設(shè)計(jì)值、UR-S11規(guī)范值及SESAM預(yù)報(bào)值的差異。

    1 WALCS波浪載荷計(jì)算原理

    WALCS 軟件可用來(lái)計(jì)算低速航行于有限和無(wú)限水深規(guī)則波中的船舶或海洋平臺(tái)的運(yùn)動(dòng)、脈動(dòng)壓力及剖面載荷的穩(wěn)態(tài)響應(yīng)。

    基于線性勢(shì)流理論,流場(chǎng)總速度勢(shì)φ(x,y,z,t)[5]可表示為

    (1)

    式(1)中:φT(x,y,z)=φI(x,y,z)+φD(x,y,z)+φR(x,y,z);φS為定常興波勢(shì);U為航速;ω為遭遇頻率;φT為非定常速度勢(shì);φI為入射勢(shì),已知;φD和φR分別為繞射勢(shì)及輻射勢(shì),均為非定常擾動(dòng)勢(shì)。

    頻域內(nèi)低航速船舶運(yùn)動(dòng)非定常擾動(dòng)勢(shì)需滿足以下定解條件。

    1) 域內(nèi)條件[L]

    (2)

    2) 線性自由面條件[F]

    (3)

    3) 物面條件[S]

    (4)

    4) 底部條件[B]

    (5)

    5) 遠(yuǎn)方輻射條件[R]:相應(yīng)的遠(yuǎn)方輻射條件。

    式(2)~式(5)中:φj(j=1,2,…,6)為單位復(fù)速度的j模式運(yùn)動(dòng)輻射勢(shì);φ7為單位波幅的繞射勢(shì);v=ω2/g;n為船舶表面指向內(nèi)部的單位法線向量。

    根據(jù)剛體動(dòng)力學(xué)原理導(dǎo)出規(guī)則波中船舶運(yùn)動(dòng)微分方程為

    [K]{η(t)}={f(t)}

    (6)

    式(6)中:[M]為船舶質(zhì)量矩陣;[Ma]為附加質(zhì)量矩陣;[C]為阻尼矩陣;[K]為靜恢復(fù)力矩陣;{η(t)}為六自由度剛體運(yùn)動(dòng);{f(t)}為波浪干擾力。

    采用高斯消去法對(duì)式(6)進(jìn)行求解,得到船舶運(yùn)動(dòng)響應(yīng)的穩(wěn)態(tài)解。利用線性化的伯努利方程,并計(jì)入靜水壓力變化部分的貢獻(xiàn),得到作用于船體濕表面上總的脈動(dòng)壓力為

    P(x,y,z,t)=Re{P(x,y,z)eiωt}

    (7)

    求解船體運(yùn)動(dòng)響應(yīng)和脈動(dòng)壓力之后,應(yīng)用達(dá)朗貝爾原理計(jì)算船體剖面內(nèi)的波浪誘導(dǎo)力和力矩。

    2 實(shí)船計(jì)算

    2.1計(jì)算模型

    WALCS采用基于自由面格林函數(shù)的三維頻域方法,通過(guò)在船體濕表面分布奇點(diǎn)來(lái)求解流場(chǎng)速度勢(shì),要求創(chuàng)建水動(dòng)力面元模型。模型網(wǎng)格的劃分以四邊形單元為主,其大小取為4倍的縱骨間距,網(wǎng)格總計(jì)1 968個(gè)。定義船殼網(wǎng)格為濕表面,指定波浪水動(dòng)壓力的矢量法向?yàn)橛闪黧w指向船殼表面。圖1為VLGC水動(dòng)力面元模型。

    圖1 VLGC水動(dòng)力面元模型

    計(jì)算波浪載荷時(shí)需提供準(zhǔn)確的質(zhì)量模型。為計(jì)算剖面載荷,采用WALCS提供的分段質(zhì)量模型方式創(chuàng)建整船質(zhì)量模型。

    2.2傳遞函數(shù)

    為充分觀察不同浪向角的波浪對(duì)船體運(yùn)動(dòng)和波浪載荷的影響規(guī)律,共選取13個(gè)浪向角(0°~180°,間隔15°,各浪向的概率分布認(rèn)為是均勻的)、29個(gè)波浪頻率(0.10~1.50 rad/s,步長(zhǎng)取0.05),水深為無(wú)限水深;采用臨界阻尼法對(duì)橫搖阻尼進(jìn)行修正,臨界阻尼系數(shù)取0.05;沿船長(zhǎng)從艉垂線至艏垂線選取21個(gè)剖面作為計(jì)算剖面。選取液艙裝載率為100%時(shí)的典型工況作為計(jì)算工況,不存在液艙晃蕩問(wèn)題。

    船體運(yùn)動(dòng)和主要載荷控制參數(shù)的傳遞函數(shù)(RAO)見(jiàn)圖2~圖4。

    圖2中:不同浪向下的垂蕩和橫搖幅頻響應(yīng)曲線吻合程度均較低,表明浪向?qū)Υw運(yùn)動(dòng)響應(yīng)的影響較大;橫浪時(shí),垂蕩和橫搖響應(yīng)達(dá)到最大。圖3中,垂向波浪彎矩幅頻響應(yīng)曲線為單峰曲線,峰值出現(xiàn)在舯剖面附近。圖4中,垂向波浪剪力幅頻響應(yīng)曲線為雙峰曲線,峰值出現(xiàn)在距艉垂線1/4船長(zhǎng)和3/4船長(zhǎng)的剖面附近。載荷峰值均出現(xiàn)在浪向角為0°,頻率為0.50~0.55 rad/s時(shí),即船舶在迎浪航行、遭遇波長(zhǎng)/船長(zhǎng)為0.95~1.15的波浪時(shí),垂向波浪彎矩和垂向波浪剪力響應(yīng)最大??梢?jiàn),WALCS預(yù)報(bào)值與實(shí)際規(guī)律一致。

    a) 垂蕩

    b) 橫搖

    a) 不同浪向(遭遇頻率為0.5 rad/s)

    b) 不同頻率(迎浪)

    a) 不同浪向

    b) 不同頻率

    2.3長(zhǎng)期預(yù)報(bào)

    長(zhǎng)期預(yù)報(bào)中,選取北大西洋波浪散布圖作為海浪統(tǒng)計(jì)資料,采用ISSC雙參數(shù)譜模擬散布圖中的海況,運(yùn)用威布爾分布擬合長(zhǎng)期分布,分別考慮長(zhǎng)峰波和短峰波2種波浪形式;短峰波的擴(kuò)散函數(shù)采用f(x)=2cos2θ/π,其中θ為波系相對(duì)于主浪向的夾角。

    對(duì)第二液艙重心處的垂向加速度和舯剖面的垂向波浪彎矩進(jìn)行長(zhǎng)期預(yù)報(bào),計(jì)算結(jié)果見(jiàn)圖5。

    圖5中,在北大西洋海況下,垂向加速度和垂向波浪彎矩的長(zhǎng)期預(yù)報(bào)極值與超越概率的對(duì)數(shù)的負(fù)值呈正比關(guān)系。

    3 結(jié)果比較與分析

    3.1加速度

    a) 垂向加速度(組合浪向)

    b) 垂向波浪彎矩(組合浪向)

    《國(guó)際散裝運(yùn)輸液化氣體船舶構(gòu)造和設(shè)備規(guī)則》(IGC規(guī)則)給出加速度分量的指導(dǎo)公式是船舶在北大西洋中以概率水平為10-8運(yùn)動(dòng)而產(chǎn)生的加速度,適用于船長(zhǎng)>50 m的船舶。[6]指導(dǎo)公式為

    (8)

    (9)

    (10)

    表1 加速度的IGC設(shè)計(jì)值和WALCS預(yù)報(bào)值 g

    由表1可知:對(duì)比3個(gè)方向的加速度計(jì)算結(jié)果,縱向最小,垂向其次,橫向最大;WALCS的縱向加速度預(yù)報(bào)值遠(yuǎn)小于IGC設(shè)計(jì)值,最大偏差為50.0%;橫向加速度的IGC設(shè)計(jì)值位于WALCS在10-7~10-8概率水平下的預(yù)報(bào)值區(qū)間內(nèi);WALCS在 10-8概率水平下的垂向加速度預(yù)報(bào)值略低于IGC設(shè)計(jì)值,偏差為5.2%。經(jīng)分析,IGC指導(dǎo)公式中的加速度與航速呈正比,選取設(shè)計(jì)航速計(jì)算,而WALCS選取零航速計(jì)算,兩者的差異對(duì)縱向加速度計(jì)算的影響較大?;诰€性頻域方法的WALCS在高航速下的縱向加速度預(yù)報(bào)精度還有待通過(guò)波浪載荷模型試驗(yàn)進(jìn)一步驗(yàn)證,橫向加速度和垂向加速度的計(jì)算結(jié)果在可接受的范圍內(nèi)。

    3.2垂向波浪彎矩

    國(guó)際船級(jí)社協(xié)會(huì)(International Association of Classification Societies,IACS)于1991年提出UR-S11,即對(duì)常規(guī)船舶的總縱強(qiáng)度提出統(tǒng)一要求。[7]該要求規(guī)定舯部垂向波浪彎矩設(shè)計(jì)值按以下公式求出。

    中垂Ms=0.11CL2B(Cb+0.7)

    (11)

    中拱Mh=0.19CL2BCb

    (12)

    式(11)和式(12)中:L,B和Cb分別為船長(zhǎng)、船寬及方形系數(shù);系數(shù)C由式(13)確定。

    (13)

    由式(13)可知,IACS設(shè)計(jì)值僅與L,B和Cb有關(guān),未計(jì)入船體線型、質(zhì)量分布和航行條件等對(duì)船體載荷有重要影響的因素,相當(dāng)于大量船舶的統(tǒng)計(jì)平均值。[8]

    表2為垂向波浪彎矩的規(guī)范計(jì)算值和直接預(yù)報(bào)值;圖6為WADAM計(jì)算的10-8概率水平下舯部垂向波浪彎矩的長(zhǎng)期預(yù)報(bào)極值曲線;圖7為WASIM計(jì)算的規(guī)則波激勵(lì)下舯部垂向合成彎矩時(shí)歷曲線。

    表2 垂向波浪彎矩的規(guī)范計(jì)算值和長(zhǎng)期預(yù)報(bào)值 MN·m

    圖6 10-8概率水平下舯部垂向波浪彎矩長(zhǎng)期預(yù)報(bào)極值曲線

    圖7 規(guī)則波激勵(lì)下的舯部垂向合成彎矩時(shí)歷曲線

    由表2可知:采用線性方法預(yù)報(bào),3種計(jì)算程序的長(zhǎng)峰波預(yù)報(bào)值比較吻合,偏差不超過(guò)3%;與長(zhǎng)峰波相比,彎矩的短峰波預(yù)報(bào)值在10-6概率水平和10-8概率水平下均要小約11%,該結(jié)論與文獻(xiàn)[8]中的計(jì)算結(jié)果一致;UR-S11規(guī)范值均位于線性程序在零航速10-6~10-8概率水平下的長(zhǎng)峰波預(yù)報(bào)值區(qū)間內(nèi),UR-S11規(guī)范計(jì)算的均值與10-8概率水平下的短峰波預(yù)報(bào)值更為接近;考慮自由液面、瞬時(shí)效應(yīng)和航速等非線性因素之后,中垂?fàn)顟B(tài)下的垂向波浪彎矩值大于中拱狀態(tài)下的垂向波浪彎矩值;WASIM非線性預(yù)報(bào)值最大,其均值與10-8概率水平下的長(zhǎng)峰波預(yù)報(bào)值較為接近。

    4 結(jié)束語(yǔ)

    1) 采用WALCS對(duì)VLGC的主要載荷控制參數(shù)進(jìn)行長(zhǎng)期預(yù)報(bào),預(yù)報(bào)值與UR-S11規(guī)范值和SESAM預(yù)報(bào)值具有可比性,計(jì)算結(jié)果適用于VLGC結(jié)構(gòu)安全性評(píng)估。

    2) IGC規(guī)則要求縱向限位支座應(yīng)承受液艙整體向前的沖擊力,其加速度為0.5g,遠(yuǎn)大于船體縱向加速度。因此,選取0.5g的縱向加速度作為結(jié)構(gòu)安全性評(píng)估的主要載荷,在實(shí)際工程應(yīng)用中是偏于安全的。

    3) 實(shí)際海況多為具有明顯主方向的短峰波,長(zhǎng)峰波中的載荷預(yù)報(bào)值往往偏大,且其大小與考慮非線性因素后的預(yù)報(bào)值相當(dāng)。建議VLGC船體載荷用長(zhǎng)峰波計(jì)算,這樣處理是趨于保守的。

    4) 考慮非線性因素后的VLGC中拱和中垂波浪彎矩并不相等,頻域線性計(jì)算方法無(wú)法考慮該特點(diǎn)。因此,有必要對(duì)其進(jìn)行修正或采取時(shí)域非線性計(jì)算方法。

    [1] 李海洲.江南造船升級(jí)版VLGC“破繭成蝶”[N].中國(guó)船舶報(bào),2014-04-04(8).

    [2] 胡可一,李小靈.中國(guó)首艘VLGC 的設(shè)計(jì)和建造[J].中國(guó)造船,2014,55(4): 132-139.

    [3] 張海彬,趙耕賢.水動(dòng)力分析在海洋結(jié)構(gòu)物設(shè)計(jì)中的應(yīng)用[J].中國(guó)海洋平臺(tái),2008,23(1):1-6.

    [4] 李輝,張藝瀚,任慧龍.船舶波浪載荷長(zhǎng)短期預(yù)報(bào)方法及其特性研究[J].華中科技大學(xué)學(xué)報(bào)(自然科學(xué)版),2013,41(12): 112-116.

    [5] 戴仰山,沈進(jìn)威,宋競(jìng)正.船舶波浪載荷[M].北京:國(guó)防工業(yè)出版社,2007:26-45.

    [6] IMO. International Code for the Construction and Equipment of Ships Carrying Liquefied Gases in Bulk(IGC Code)[S]. 2006.

    [7] 顧學(xué)康,陳瑞章,沈進(jìn)威,等.海浪統(tǒng)計(jì)資料對(duì)船舶波浪彎矩設(shè)計(jì)值的影響及其意義[J].船舶力學(xué),1998, 2(5):50-62.

    [8] 陳瑞章,沈進(jìn)威.波浪彎矩設(shè)計(jì)值與長(zhǎng)期預(yù)報(bào)理論計(jì)算值的比較[J].中國(guó)造船,1996(4):41-44.

    WaveLoadPredictionandAnalysisofVeryLargeGasCarrier

    ZHOUQinghua,LIANGXiaojun,LIXiaoling,CHENBing

    (Jiangnan Shipyard (Group) Co., Ltd., Shanghai 201913, China)

    U661.1

    A

    2016-04-21

    工業(yè)和信息化部高技術(shù)船舶科研項(xiàng)目(工信部聯(lián)裝[2012]534號(hào))

    周清華(1984—),男,湖北監(jiān)利人,工程師,碩士,從事船舶與海洋結(jié)構(gòu)物設(shè)計(jì)制造研究。E-mail:supershiper@163.com

    1000-4653(2016)03-0068-05

    猜你喜歡
    液艙船體波浪
    B型LNG液艙支座縱骨趾端處表面裂紋擴(kuò)展計(jì)算
    船體行駛過(guò)程中的壓力監(jiān)測(cè)方法
    波浪谷和波浪巖
    基于CFD的大型船舶液艙晃蕩研究
    波浪谷隨想
    去看神奇波浪谷
    考慮晃蕩效應(yīng)的獨(dú)立B型LNG液艙結(jié)構(gòu)多目標(biāo)優(yōu)化
    海洋工程(2016年2期)2016-10-12 05:08:07
    FPSO與運(yùn)輸船旁靠時(shí)液艙晃蕩與船舶運(yùn)動(dòng)耦合效應(yīng)分析
    焊接殘余應(yīng)力對(duì)船體結(jié)構(gòu)疲勞強(qiáng)度的影響分析
    焊接(2015年9期)2015-07-18 11:03:51
    赴美軍“仁慈”號(hào)醫(yī)院船駐船體會(huì)
    岛国在线观看网站| 国产亚洲欧美在线一区二区| 天堂俺去俺来也www色官网| 中国国产av一级| 欧美xxⅹ黑人| 天天添夜夜摸| 久久中文字幕一级| 国产在线视频一区二区| 99精国产麻豆久久婷婷| 久久精品成人免费网站| 老司机午夜十八禁免费视频| 国产片内射在线| 国产高清国产精品国产三级| 丰满迷人的少妇在线观看| www.精华液| 黑人欧美特级aaaaaa片| 一本综合久久免费| 麻豆乱淫一区二区| 国产精品久久久久久精品古装| 秋霞在线观看毛片| 亚洲精品中文字幕在线视频| 午夜福利视频在线观看免费| 午夜福利乱码中文字幕| 久热这里只有精品99| 国产精品国产av在线观看| 久热这里只有精品99| 国产免费现黄频在线看| 国产精品av久久久久免费| 亚洲精品日韩在线中文字幕| 国产精品av久久久久免费| 中文字幕另类日韩欧美亚洲嫩草| 成人黄色视频免费在线看| 亚洲av国产av综合av卡| 国产黄色免费在线视频| 免费观看a级毛片全部| 欧美日韩福利视频一区二区| 一区二区av电影网| 热99re8久久精品国产| 久久久久国产一级毛片高清牌| 成人三级做爰电影| 黑人猛操日本美女一级片| 搡老乐熟女国产| 搡老乐熟女国产| 亚洲 欧美一区二区三区| 久久 成人 亚洲| 欧美乱码精品一区二区三区| 老司机在亚洲福利影院| 精品人妻在线不人妻| 欧美xxⅹ黑人| 高清av免费在线| 我要看黄色一级片免费的| av天堂在线播放| www.av在线官网国产| 中文欧美无线码| 国产人伦9x9x在线观看| 亚洲精品国产区一区二| 国产欧美亚洲国产| 精品一区二区三区四区五区乱码| 9191精品国产免费久久| tube8黄色片| 91成年电影在线观看| 中文字幕人妻熟女乱码| 亚洲精品av麻豆狂野| 老司机靠b影院| 日本精品一区二区三区蜜桃| 嫁个100分男人电影在线观看| 免费久久久久久久精品成人欧美视频| 一级黄色大片毛片| 18禁国产床啪视频网站| 国产一区二区三区av在线| 搡老乐熟女国产| 一边摸一边做爽爽视频免费| 久久久水蜜桃国产精品网| 91国产中文字幕| 亚洲国产精品一区三区| 妹子高潮喷水视频| 性高湖久久久久久久久免费观看| 在线av久久热| 亚洲一区中文字幕在线| 在线永久观看黄色视频| 成年人黄色毛片网站| 男女国产视频网站| 亚洲精品中文字幕在线视频| 91成年电影在线观看| 国产精品秋霞免费鲁丝片| 日本av手机在线免费观看| 另类精品久久| 国产三级黄色录像| 亚洲男人天堂网一区| 久久久久久久国产电影| 亚洲精品久久午夜乱码| 啦啦啦中文免费视频观看日本| 久久久国产一区二区| 丰满人妻熟妇乱又伦精品不卡| 日本av手机在线免费观看| 我要看黄色一级片免费的| 五月天丁香电影| 亚洲一卡2卡3卡4卡5卡精品中文| 日本vs欧美在线观看视频| 一级毛片电影观看| 真人做人爱边吃奶动态| 国产老妇伦熟女老妇高清| 不卡av一区二区三区| 男女高潮啪啪啪动态图| 丰满少妇做爰视频| 亚洲精品国产av蜜桃| 久久毛片免费看一区二区三区| 精品国产乱码久久久久久男人| 别揉我奶头~嗯~啊~动态视频 | 中文字幕精品免费在线观看视频| 亚洲 国产 在线| 欧美黑人精品巨大| 免费观看av网站的网址| bbb黄色大片| 午夜精品久久久久久毛片777| 国产成人一区二区三区免费视频网站| 脱女人内裤的视频| 国产视频一区二区在线看| 国产淫语在线视频| 美女脱内裤让男人舔精品视频| 久久精品久久久久久噜噜老黄| 久久热在线av| 午夜激情av网站| 国产xxxxx性猛交| 五月开心婷婷网| 免费人妻精品一区二区三区视频| 欧美另类亚洲清纯唯美| 99九九在线精品视频| 性色av一级| 亚洲精品一区蜜桃| 一二三四社区在线视频社区8| 别揉我奶头~嗯~啊~动态视频 | av不卡在线播放| 国产黄频视频在线观看| 首页视频小说图片口味搜索| 免费久久久久久久精品成人欧美视频| 亚洲黑人精品在线| 亚洲情色 制服丝袜| 精品亚洲成a人片在线观看| 亚洲成人手机| 精品第一国产精品| 国产精品欧美亚洲77777| 亚洲激情五月婷婷啪啪| av有码第一页| 黄色视频,在线免费观看| 男女床上黄色一级片免费看| 老汉色av国产亚洲站长工具| 亚洲av美国av| avwww免费| 巨乳人妻的诱惑在线观看| 国产精品一区二区精品视频观看| 少妇裸体淫交视频免费看高清 | 国产精品欧美亚洲77777| 欧美精品人与动牲交sv欧美| 久久天堂一区二区三区四区| 蜜桃国产av成人99| 91精品伊人久久大香线蕉| 国产亚洲欧美精品永久| 菩萨蛮人人尽说江南好唐韦庄| 美女大奶头黄色视频| 亚洲精品国产精品久久久不卡| 国产精品麻豆人妻色哟哟久久| 国产精品一二三区在线看| 一级黄色大片毛片| 又黄又粗又硬又大视频| 精品久久久精品久久久| 一本色道久久久久久精品综合| bbb黄色大片| 在线观看免费日韩欧美大片| 国产黄色免费在线视频| 手机成人av网站| 精品少妇久久久久久888优播| 久久影院123| 欧美大码av| 天堂8中文在线网| 丰满人妻熟妇乱又伦精品不卡| 久久香蕉激情| 男女免费视频国产| 亚洲精品一二三| 精品视频人人做人人爽| 国产精品国产av在线观看| 免费在线观看完整版高清| 日本黄色日本黄色录像| 欧美老熟妇乱子伦牲交| www.999成人在线观看| 十八禁网站免费在线| 日韩视频一区二区在线观看| 电影成人av| 亚洲熟女精品中文字幕| 性少妇av在线| 丰满人妻熟妇乱又伦精品不卡| 日日夜夜操网爽| 国产精品熟女久久久久浪| 欧美日韩黄片免| 新久久久久国产一级毛片| 亚洲成av片中文字幕在线观看| 免费观看人在逋| tocl精华| 精品久久久久久电影网| 亚洲黑人精品在线| 欧美精品人与动牲交sv欧美| 女性生殖器流出的白浆| 久久毛片免费看一区二区三区| 一区福利在线观看| 老熟女久久久| 啦啦啦啦在线视频资源| 免费少妇av软件| 美女午夜性视频免费| 80岁老熟妇乱子伦牲交| 丝袜美腿诱惑在线| 伊人亚洲综合成人网| 操出白浆在线播放| 中文精品一卡2卡3卡4更新| 日日爽夜夜爽网站| 国产精品免费视频内射| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品二区激情视频| 最新的欧美精品一区二区| 亚洲综合色网址| 黄色a级毛片大全视频| 老司机影院成人| 亚洲 国产 在线| 国产精品久久久人人做人人爽| 叶爱在线成人免费视频播放| 亚洲七黄色美女视频| 青草久久国产| 极品人妻少妇av视频| 日韩一卡2卡3卡4卡2021年| 亚洲精品乱久久久久久| 少妇的丰满在线观看| 欧美一级毛片孕妇| 亚洲成国产人片在线观看| 少妇人妻久久综合中文| e午夜精品久久久久久久| 啪啪无遮挡十八禁网站| 69av精品久久久久久 | 国产成人系列免费观看| 精品视频人人做人人爽| 深夜精品福利| 黑丝袜美女国产一区| 色精品久久人妻99蜜桃| 成年人黄色毛片网站| 少妇的丰满在线观看| 男女高潮啪啪啪动态图| 99国产精品99久久久久| 少妇 在线观看| 999久久久国产精品视频| 波多野结衣av一区二区av| 宅男免费午夜| 在线看a的网站| 日本撒尿小便嘘嘘汇集6| 亚洲精品久久午夜乱码| 亚洲精品成人av观看孕妇| 国产又爽黄色视频| 90打野战视频偷拍视频| 美女主播在线视频| 国产成人啪精品午夜网站| 久久精品国产综合久久久| 水蜜桃什么品种好| 国产高清视频在线播放一区 | 满18在线观看网站| 国产av又大| 欧美国产精品一级二级三级| 下体分泌物呈黄色| 国产日韩欧美在线精品| 亚洲精品美女久久久久99蜜臀| 视频区欧美日本亚洲| 亚洲av欧美aⅴ国产| 国产欧美日韩综合在线一区二区| av又黄又爽大尺度在线免费看| 我的亚洲天堂| 久久人人爽人人片av| 亚洲av成人不卡在线观看播放网 | 成人国产av品久久久| 亚洲av欧美aⅴ国产| 国产精品亚洲av一区麻豆| 免费在线观看日本一区| 激情视频va一区二区三区| 久久国产精品影院| av免费在线观看网站| xxxhd国产人妻xxx| 免费在线观看日本一区| 国产男女内射视频| 久久性视频一级片| 日韩大码丰满熟妇| 久久人妻熟女aⅴ| 国产黄色免费在线视频| 亚洲欧美成人综合另类久久久| 久久影院123| 午夜福利视频精品| 黄色 视频免费看| 亚洲,欧美精品.| 中亚洲国语对白在线视频| 国产精品亚洲av一区麻豆| 色播在线永久视频| 久久久久久久国产电影| 美女高潮喷水抽搐中文字幕| 色精品久久人妻99蜜桃| 欧美中文综合在线视频| 亚洲欧美精品综合一区二区三区| 国产精品香港三级国产av潘金莲| 天天操日日干夜夜撸| 18禁裸乳无遮挡动漫免费视频| 亚洲av片天天在线观看| 老汉色∧v一级毛片| 日本91视频免费播放| h视频一区二区三区| 黄色视频不卡| 满18在线观看网站| 热99re8久久精品国产| 日韩大片免费观看网站| 淫妇啪啪啪对白视频 | 欧美亚洲日本最大视频资源| 亚洲第一av免费看| 五月开心婷婷网| 大香蕉久久成人网| 国产男人的电影天堂91| 精品少妇一区二区三区视频日本电影| avwww免费| 男女之事视频高清在线观看| 精品一区二区三区四区五区乱码| 黑人欧美特级aaaaaa片| 老熟妇乱子伦视频在线观看 | 国产伦理片在线播放av一区| 欧美在线一区亚洲| 国产成人欧美在线观看 | 男女无遮挡免费网站观看| 亚洲av电影在线进入| 日韩一卡2卡3卡4卡2021年| netflix在线观看网站| 极品少妇高潮喷水抽搐| videos熟女内射| 精品国产一区二区三区久久久樱花| 久久这里只有精品19| 9热在线视频观看99| www.精华液| 成年人午夜在线观看视频| 国产不卡av网站在线观看| 免费高清在线观看视频在线观看| 久久国产亚洲av麻豆专区| 欧美精品高潮呻吟av久久| 啪啪无遮挡十八禁网站| 午夜免费鲁丝| 亚洲少妇的诱惑av| 亚洲av欧美aⅴ国产| 精品国产乱子伦一区二区三区 | 国产亚洲精品久久久久5区| √禁漫天堂资源中文www| 性色av一级| 大香蕉久久网| 在线av久久热| 一区二区三区激情视频| 国产男女超爽视频在线观看| 亚洲成国产人片在线观看| 久久国产精品影院| 国产亚洲av片在线观看秒播厂| 久久亚洲精品不卡| 色婷婷av一区二区三区视频| 激情视频va一区二区三区| 国产亚洲精品久久久久5区| 国产成人免费无遮挡视频| 少妇被粗大的猛进出69影院| 夜夜夜夜夜久久久久| 一边摸一边做爽爽视频免费| 日韩视频一区二区在线观看| 多毛熟女@视频| 超碰97精品在线观看| 丝袜美腿诱惑在线| 黄色a级毛片大全视频| 国产一区有黄有色的免费视频| 在线观看免费日韩欧美大片| 又黄又粗又硬又大视频| 国产欧美亚洲国产| av一本久久久久| 久久精品国产a三级三级三级| 国产成人精品久久二区二区91| 黄色毛片三级朝国网站| 久久天堂一区二区三区四区| 最新在线观看一区二区三区| 午夜激情av网站| 日韩,欧美,国产一区二区三区| 性少妇av在线| 日韩 亚洲 欧美在线| 国产亚洲欧美在线一区二区| √禁漫天堂资源中文www| 韩国精品一区二区三区| 欧美 亚洲 国产 日韩一| 午夜福利一区二区在线看| 一进一出抽搐动态| 九色亚洲精品在线播放| 免费在线观看完整版高清| 中文字幕人妻熟女乱码| 一级片免费观看大全| 亚洲性夜色夜夜综合| 精品欧美一区二区三区在线| 久久久久久人人人人人| 汤姆久久久久久久影院中文字幕| www日本在线高清视频| 日韩熟女老妇一区二区性免费视频| 中国国产av一级| 男人添女人高潮全过程视频| 一区二区av电影网| 亚洲国产成人一精品久久久| 国产一区二区三区av在线| 欧美精品啪啪一区二区三区 | 一区二区三区精品91| 精品福利观看| 黑丝袜美女国产一区| 久久久久久久精品精品| 男人添女人高潮全过程视频| 黄色片一级片一级黄色片| 麻豆乱淫一区二区| av超薄肉色丝袜交足视频| 精品少妇内射三级| 日韩视频在线欧美| 美女主播在线视频| 免费日韩欧美在线观看| 亚洲成av片中文字幕在线观看| 婷婷丁香在线五月| 美女视频免费永久观看网站| 最近最新中文字幕大全免费视频| 淫妇啪啪啪对白视频 | 午夜福利在线观看吧| 一本久久精品| 国产一区二区激情短视频 | 美女高潮喷水抽搐中文字幕| 日本猛色少妇xxxxx猛交久久| 三上悠亚av全集在线观看| 永久免费av网站大全| 免费一级毛片在线播放高清视频 | 啦啦啦啦在线视频资源| 天堂中文最新版在线下载| 一区二区三区四区激情视频| 欧美午夜高清在线| 亚洲第一欧美日韩一区二区三区 | 欧美精品亚洲一区二区| 日韩欧美免费精品| 久久 成人 亚洲| 人妻一区二区av| www.精华液| 国产精品麻豆人妻色哟哟久久| 精品国产一区二区久久| 国产精品一二三区在线看| 啪啪无遮挡十八禁网站| 国产成人欧美在线观看 | 日韩视频一区二区在线观看| 啦啦啦视频在线资源免费观看| 99re6热这里在线精品视频| 亚洲欧洲精品一区二区精品久久久| 一区福利在线观看| av又黄又爽大尺度在线免费看| 久久精品国产a三级三级三级| 最黄视频免费看| 亚洲熟女毛片儿| 亚洲av欧美aⅴ国产| 五月天丁香电影| 黄频高清免费视频| 国产成+人综合+亚洲专区| 欧美激情极品国产一区二区三区| 在线av久久热| 十分钟在线观看高清视频www| 欧美日韩亚洲高清精品| 国产一区二区三区av在线| 国产视频一区二区在线看| 最近最新免费中文字幕在线| 两性午夜刺激爽爽歪歪视频在线观看 | 日韩大片免费观看网站| 黑人巨大精品欧美一区二区蜜桃| 欧美日韩视频精品一区| 久久国产精品影院| 天堂中文最新版在线下载| 中国国产av一级| 男女无遮挡免费网站观看| 国产成人影院久久av| 手机成人av网站| 欧美日韩亚洲综合一区二区三区_| 中国国产av一级| 国产精品久久久人人做人人爽| 老司机午夜福利在线观看视频 | www日本在线高清视频| 女性生殖器流出的白浆| 午夜免费鲁丝| 欧美黑人精品巨大| 精品欧美一区二区三区在线| 欧美+亚洲+日韩+国产| 免费看十八禁软件| 国产精品熟女久久久久浪| 日本a在线网址| 免费一级毛片在线播放高清视频 | 90打野战视频偷拍视频| 岛国在线观看网站| 成年人黄色毛片网站| 人人妻,人人澡人人爽秒播| 日韩中文字幕视频在线看片| 欧美激情 高清一区二区三区| 精品高清国产在线一区| 国产精品久久久久久人妻精品电影 | 欧美日韩亚洲高清精品| 可以免费在线观看a视频的电影网站| 久久精品aⅴ一区二区三区四区| 黄色怎么调成土黄色| 精品国产一区二区三区久久久樱花| 精品福利观看| 欧美日韩一级在线毛片| 少妇精品久久久久久久| 久久综合国产亚洲精品| 90打野战视频偷拍视频| 欧美精品亚洲一区二区| bbb黄色大片| 美女中出高潮动态图| 色视频在线一区二区三区| 亚洲九九香蕉| 亚洲成人免费电影在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 国产日韩欧美在线精品| 久久久久久久久免费视频了| 香蕉丝袜av| 亚洲一区中文字幕在线| 久9热在线精品视频| 涩涩av久久男人的天堂| 久久午夜综合久久蜜桃| 大陆偷拍与自拍| 精品福利观看| 国产极品粉嫩免费观看在线| 精品一品国产午夜福利视频| 日日夜夜操网爽| 满18在线观看网站| 飞空精品影院首页| 色综合欧美亚洲国产小说| 亚洲久久久国产精品| 一级毛片精品| 日韩熟女老妇一区二区性免费视频| bbb黄色大片| 免费一级毛片在线播放高清视频 | 99香蕉大伊视频| 丝袜在线中文字幕| 日本av免费视频播放| 老汉色∧v一级毛片| 麻豆国产av国片精品| 亚洲国产成人一精品久久久| 99热国产这里只有精品6| 自线自在国产av| 99久久精品国产亚洲精品| 精品亚洲成a人片在线观看| 日本vs欧美在线观看视频| 午夜福利视频在线观看免费| 亚洲精品乱久久久久久| 搡老岳熟女国产| 美女扒开内裤让男人捅视频| 日韩电影二区| 欧美中文综合在线视频| 久久久欧美国产精品| 成人18禁高潮啪啪吃奶动态图| 中文字幕人妻丝袜一区二区| 亚洲色图综合在线观看| 中文字幕精品免费在线观看视频| 久久99一区二区三区| 国产精品99久久99久久久不卡| 汤姆久久久久久久影院中文字幕| 国产无遮挡羞羞视频在线观看| 日本wwww免费看| 国产亚洲av片在线观看秒播厂| 国产亚洲午夜精品一区二区久久| 99re6热这里在线精品视频| 国产成人免费无遮挡视频| 女性生殖器流出的白浆| 老熟妇乱子伦视频在线观看 | 亚洲欧美激情在线| 成人国语在线视频| 久久久久精品人妻al黑| 女性被躁到高潮视频| 69精品国产乱码久久久| 日韩欧美一区视频在线观看| 99久久99久久久精品蜜桃| 桃红色精品国产亚洲av| 午夜影院在线不卡| 日日摸夜夜添夜夜添小说| 久久天躁狠狠躁夜夜2o2o| 国产深夜福利视频在线观看| 中国国产av一级| 51午夜福利影视在线观看| 悠悠久久av| svipshipincom国产片| 精品久久久久久久毛片微露脸 | 交换朋友夫妻互换小说| 日韩制服丝袜自拍偷拍| 青春草亚洲视频在线观看| 国产男女超爽视频在线观看| 亚洲精品在线美女| 久久久国产一区二区| 91成年电影在线观看| 精品少妇一区二区三区视频日本电影| 欧美黄色淫秽网站| 亚洲精品美女久久av网站| 女性被躁到高潮视频| 国产成人av教育| 少妇被粗大的猛进出69影院| 在线观看舔阴道视频| 久久性视频一级片| 青春草视频在线免费观看| 亚洲国产欧美网| 亚洲美女黄色视频免费看| 国产xxxxx性猛交| 国产成人精品久久二区二区免费| 免费人妻精品一区二区三区视频| 欧美在线黄色| 精品卡一卡二卡四卡免费| 欧美成狂野欧美在线观看| 搡老岳熟女国产| 热re99久久国产66热| 国产色视频综合| 久久亚洲精品不卡| 9191精品国产免费久久| 18在线观看网站| 欧美av亚洲av综合av国产av|