• <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ì)
    色综合亚洲欧美另类图片| 狂野欧美激情性xxxx在线观看| 国产一区二区亚洲精品在线观看| 麻豆国产97在线/欧美| 乱人视频在线观看| 国产精品爽爽va在线观看网站| 中文字幕熟女人妻在线| 久久精品影院6| 亚洲人成网站在线播| av专区在线播放| 国产探花极品一区二区| 特级一级黄色大片| 九色国产91popny在线| 赤兔流量卡办理| 免费av不卡在线播放| 久久久久国产精品人妻aⅴ院| 啦啦啦观看免费观看视频高清| 日本黄色片子视频| 国产av一区在线观看免费| 亚洲在线观看片| 亚洲午夜理论影院| 免费无遮挡裸体视频| 精品一区二区三区视频在线观看免费| 久久久久久九九精品二区国产| 看黄色毛片网站| 国产免费男女视频| 欧美不卡视频在线免费观看| 久久久久久久久久久丰满 | 特大巨黑吊av在线直播| 免费看av在线观看网站| 老师上课跳d突然被开到最大视频| 亚洲在线观看片| 91在线观看av| 变态另类成人亚洲欧美熟女| 亚洲专区国产一区二区| 日本黄大片高清| 久久热精品热| 亚洲真实伦在线观看| 神马国产精品三级电影在线观看| 久久国产乱子免费精品| 99热只有精品国产| 国产精品无大码| 久久精品影院6| 亚洲内射少妇av| 欧美日韩乱码在线| 免费大片18禁| 亚洲七黄色美女视频| 午夜福利在线观看吧| 在线免费十八禁| 中文在线观看免费www的网站| 成人欧美大片| 两性午夜刺激爽爽歪歪视频在线观看| 一级黄色大片毛片| 亚洲第一区二区三区不卡| 男女啪啪激烈高潮av片| 欧美bdsm另类| 亚洲国产欧洲综合997久久,| 搡老妇女老女人老熟妇| 亚洲欧美日韩高清专用| 别揉我奶头 嗯啊视频| 国产伦精品一区二区三区四那| 伊人久久精品亚洲午夜| 91在线观看av| 亚洲va日本ⅴa欧美va伊人久久| 国产精品人妻久久久久久| 日日摸夜夜添夜夜添小说| 色综合色国产| 国内精品美女久久久久久| 最好的美女福利视频网| 国产精品一区二区三区四区久久| 校园人妻丝袜中文字幕| 欧美一区二区国产精品久久精品| 欧美激情在线99| 免费在线观看成人毛片| 欧美丝袜亚洲另类 | 老熟妇仑乱视频hdxx| 成人国产麻豆网| 久久久久久久久大av| 亚洲av日韩精品久久久久久密| 一个人免费在线观看电影| 最近最新免费中文字幕在线| 国产精品美女特级片免费视频播放器| 91麻豆av在线| 免费人成在线观看视频色| 亚洲,欧美,日韩| 伦理电影大哥的女人| 亚洲人成网站在线播放欧美日韩| 永久网站在线| 久久热精品热| 一a级毛片在线观看| 美女高潮喷水抽搐中文字幕| 亚洲精品国产成人久久av| 狠狠狠狠99中文字幕| 久久九九热精品免费| 男人狂女人下面高潮的视频| 国产精品一区二区免费欧美| 99热精品在线国产| 欧美不卡视频在线免费观看| 黄色日韩在线| 国产爱豆传媒在线观看| 午夜a级毛片| 亚洲av日韩精品久久久久久密| 日韩欧美精品免费久久| 国产日本99.免费观看| 中文在线观看免费www的网站| 熟妇人妻久久中文字幕3abv| 身体一侧抽搐| 草草在线视频免费看| a级毛片a级免费在线| 18禁裸乳无遮挡免费网站照片| 午夜激情欧美在线| 99热这里只有是精品50| 日本五十路高清| 天美传媒精品一区二区| 真实男女啪啪啪动态图| 亚洲av中文av极速乱 | 不卡视频在线观看欧美| 国产精品,欧美在线| 亚洲最大成人手机在线| 91在线观看av| 丰满的人妻完整版| 在现免费观看毛片| 国产伦在线观看视频一区| 亚洲成人久久性| 日本五十路高清| 成人国产一区最新在线观看| 色噜噜av男人的天堂激情| 他把我摸到了高潮在线观看| 91久久精品国产一区二区三区| 在线观看午夜福利视频| 亚洲一区高清亚洲精品| 成人综合一区亚洲| 国产精品久久久久久久电影| 舔av片在线| 国产v大片淫在线免费观看| 亚洲经典国产精华液单| 丰满人妻一区二区三区视频av| 色精品久久人妻99蜜桃| 2021天堂中文幕一二区在线观| 老司机深夜福利视频在线观看| 国产熟女欧美一区二区| 无遮挡黄片免费观看| 日韩欧美三级三区| 91久久精品国产一区二区成人| 国产中年淑女户外野战色| 亚洲av成人精品一区久久| 日韩国内少妇激情av| 亚洲国产精品合色在线| 亚洲美女黄片视频| 亚州av有码| 观看美女的网站| 欧美日韩亚洲国产一区二区在线观看| 麻豆久久精品国产亚洲av| 精品久久国产蜜桃| 亚洲久久久久久中文字幕| 99视频精品全部免费 在线| 男人舔奶头视频| 最近视频中文字幕2019在线8| 国产日本99.免费观看| 日韩精品中文字幕看吧| 精品一区二区免费观看| 久久6这里有精品| 欧美国产日韩亚洲一区| 久久香蕉精品热| 国产精品无大码| 欧美xxxx黑人xx丫x性爽| 欧美一级a爱片免费观看看| 国产乱人视频| 色精品久久人妻99蜜桃| 看片在线看免费视频| 精品人妻视频免费看| 亚洲欧美日韩卡通动漫| 看片在线看免费视频| 国产午夜精品久久久久久一区二区三区 | 国产大屁股一区二区在线视频| 欧美另类亚洲清纯唯美| 看片在线看免费视频| 动漫黄色视频在线观看| 国产精品av视频在线免费观看| 校园人妻丝袜中文字幕| 亚洲欧美清纯卡通| 色精品久久人妻99蜜桃| aaaaa片日本免费| 亚洲男人的天堂狠狠| 国内精品久久久久精免费| 特级一级黄色大片| 成人av一区二区三区在线看| 亚洲自偷自拍三级| 一边摸一边抽搐一进一小说| 亚洲avbb在线观看| 久久精品91蜜桃| 丰满乱子伦码专区| 精品久久久久久久久久久久久| 一区二区三区高清视频在线| 免费搜索国产男女视频| 日韩一区二区视频免费看| 亚洲精品日韩av片在线观看| 精品人妻熟女av久视频| 精品不卡国产一区二区三区| 偷拍熟女少妇极品色| 日韩欧美国产在线观看| 日本色播在线视频| 国产精品美女特级片免费视频播放器| 国产黄片美女视频| 亚洲狠狠婷婷综合久久图片| 国产成人a区在线观看| 亚洲精品一区av在线观看| 一区二区三区免费毛片| 国产午夜精品论理片| 他把我摸到了高潮在线观看| 久久精品国产99精品国产亚洲性色| 女同久久另类99精品国产91| 亚洲性夜色夜夜综合| 高清在线国产一区| 村上凉子中文字幕在线| 嫁个100分男人电影在线观看| 精品久久久久久久人妻蜜臀av| 午夜a级毛片| 人妻少妇偷人精品九色| 久久久久国产精品人妻aⅴ院| 久久人人爽人人爽人人片va| 国产精品久久电影中文字幕| 亚洲乱码一区二区免费版| 99久久无色码亚洲精品果冻| 在线观看舔阴道视频| 91麻豆精品激情在线观看国产| 亚洲欧美激情综合另类| 99久久久亚洲精品蜜臀av| 日本a在线网址| 国产高清三级在线| 国产一区二区三区av在线 | 欧美黑人欧美精品刺激| 人妻久久中文字幕网| 亚洲一区二区三区色噜噜| 最近中文字幕高清免费大全6 | 波多野结衣高清无吗| 午夜福利在线观看吧| 成人欧美大片| 99热只有精品国产| 特大巨黑吊av在线直播| 自拍偷自拍亚洲精品老妇| 国模一区二区三区四区视频| 国产探花在线观看一区二区| www.色视频.com| 中文字幕精品亚洲无线码一区| 精品福利观看| 亚洲av第一区精品v没综合| 久久久久久久亚洲中文字幕| 久久婷婷人人爽人人干人人爱| 亚洲四区av| 欧美潮喷喷水| 国产高清不卡午夜福利| 一进一出抽搐gif免费好疼| 97人妻精品一区二区三区麻豆| 亚洲av成人av| 日本欧美国产在线视频| 精品一区二区三区av网在线观看| 久久精品91蜜桃| 精品久久国产蜜桃| 99热6这里只有精品| 亚洲性久久影院| 亚洲av免费在线观看| 欧美中文日本在线观看视频| 久久久午夜欧美精品| 欧美xxxx黑人xx丫x性爽| 亚洲无线观看免费| 床上黄色一级片| 欧美国产日韩亚洲一区| 色噜噜av男人的天堂激情| 美女 人体艺术 gogo| 中文亚洲av片在线观看爽| 午夜亚洲福利在线播放| 韩国av一区二区三区四区| 深夜精品福利| 高清在线国产一区| 日本 av在线| 在线免费观看不下载黄p国产 | 能在线免费观看的黄片| 在线观看66精品国产| 欧美日韩乱码在线| 久久久久免费精品人妻一区二区| 国产黄a三级三级三级人| 亚洲图色成人| 3wmmmm亚洲av在线观看| 桃色一区二区三区在线观看| 国产精品综合久久久久久久免费| 99九九线精品视频在线观看视频| 给我免费播放毛片高清在线观看| 成年免费大片在线观看| 欧美高清性xxxxhd video| 高清日韩中文字幕在线| 成人鲁丝片一二三区免费| 精品一区二区三区人妻视频| 最近最新免费中文字幕在线| 亚洲av五月六月丁香网| 欧美性猛交╳xxx乱大交人| 欧美另类亚洲清纯唯美| av天堂中文字幕网| 日本色播在线视频| 一个人免费在线观看电影| 欧美色欧美亚洲另类二区| 特级一级黄色大片| www.色视频.com| 99在线视频只有这里精品首页| 亚洲av中文字字幕乱码综合| 少妇人妻一区二区三区视频| 超碰av人人做人人爽久久| 亚洲在线自拍视频| 99久久精品热视频| 亚洲人成网站在线播| 国产午夜精品久久久久久一区二区三区 | 琪琪午夜伦伦电影理论片6080| 日本一本二区三区精品| 亚洲一级一片aⅴ在线观看| 99热这里只有精品一区| 色噜噜av男人的天堂激情| av女优亚洲男人天堂| 精品久久久久久久人妻蜜臀av| 欧美最新免费一区二区三区| 色噜噜av男人的天堂激情| 国产伦人伦偷精品视频| 精品日产1卡2卡| 日本-黄色视频高清免费观看| 综合色av麻豆| 午夜福利成人在线免费观看| 亚洲人成网站在线播放欧美日韩| 精品午夜福利在线看| 亚洲欧美清纯卡通| 久久人妻av系列| av天堂中文字幕网| 热99在线观看视频| 少妇裸体淫交视频免费看高清| 色av中文字幕| 国产高清不卡午夜福利| 国产中年淑女户外野战色| 熟女电影av网| 一边摸一边抽搐一进一小说| 九九在线视频观看精品| 一区二区三区免费毛片| 国产老妇女一区| 男女视频在线观看网站免费| .国产精品久久| 久久久久久久精品吃奶| 免费无遮挡裸体视频| 日本爱情动作片www.在线观看 | 免费无遮挡裸体视频| а√天堂www在线а√下载| 国产在线男女| 亚洲av成人av| 欧美潮喷喷水| 亚洲av日韩精品久久久久久密| 窝窝影院91人妻| 亚洲18禁久久av| 日本一本二区三区精品| 亚洲国产欧洲综合997久久,| 大又大粗又爽又黄少妇毛片口| 国产精品一及| 亚洲在线观看片| 99热只有精品国产| 亚洲三级黄色毛片| 麻豆国产97在线/欧美| 大又大粗又爽又黄少妇毛片口| 美女高潮喷水抽搐中文字幕| 亚洲在线自拍视频| 天堂动漫精品| 日日夜夜操网爽| 亚洲内射少妇av| 99热精品在线国产| 九九热线精品视视频播放| 亚洲在线自拍视频| 国产精品福利在线免费观看| 男女之事视频高清在线观看| 久99久视频精品免费| 99精品在免费线老司机午夜| 日本免费一区二区三区高清不卡| 少妇裸体淫交视频免费看高清| 精品久久久久久久久久久久久| 免费在线观看影片大全网站| 一夜夜www| 国产精品一区二区免费欧美| 2021天堂中文幕一二区在线观| 精品久久久久久久末码| a级一级毛片免费在线观看| 国产精品福利在线免费观看| 1024手机看黄色片| 精品人妻偷拍中文字幕| 久99久视频精品免费| 22中文网久久字幕| 美女 人体艺术 gogo| 人人妻,人人澡人人爽秒播| 国产日本99.免费观看| 自拍偷自拍亚洲精品老妇| 日本黄色片子视频| av国产免费在线观看| 日本撒尿小便嘘嘘汇集6| 2021天堂中文幕一二区在线观| 最新在线观看一区二区三区| 国产一区二区在线观看日韩| 久久久久精品国产欧美久久久| 日日撸夜夜添| 国产综合懂色| 九色成人免费人妻av| 国产欧美日韩精品一区二区| 国产一区二区三区在线臀色熟女| 久久精品国产亚洲av涩爱 | 黄色女人牲交| 精品久久久久久久久久免费视频| 精品久久国产蜜桃| 中文资源天堂在线| 免费大片18禁| 亚洲色图av天堂| 日本 欧美在线| 亚洲人成网站在线播放欧美日韩| 1024手机看黄色片| 日日啪夜夜撸| 又爽又黄无遮挡网站| 亚洲成人精品中文字幕电影| 99热这里只有精品一区| 国产精品免费一区二区三区在线| 日日夜夜操网爽| 亚洲无线在线观看| 色5月婷婷丁香| 一级毛片久久久久久久久女| 国产高清激情床上av| 久久精品国产亚洲av天美| 日日撸夜夜添| 窝窝影院91人妻| 淫秽高清视频在线观看| 特级一级黄色大片| 日本 欧美在线| 亚洲国产高清在线一区二区三| 最近在线观看免费完整版| 大型黄色视频在线免费观看| 欧美丝袜亚洲另类 | 午夜久久久久精精品| 国产精品人妻久久久影院| 日日夜夜操网爽| 亚洲狠狠婷婷综合久久图片| 久久久久国产精品人妻aⅴ院| 亚洲精华国产精华精| 一区二区三区激情视频| 成人性生交大片免费视频hd| 亚洲乱码一区二区免费版| 女的被弄到高潮叫床怎么办 | 窝窝影院91人妻| 嫩草影院精品99| 国产探花极品一区二区| 在线观看午夜福利视频| 久久久久久九九精品二区国产| 亚洲欧美激情综合另类| 一a级毛片在线观看| 99riav亚洲国产免费| 国产探花极品一区二区| 国国产精品蜜臀av免费| 久久人妻av系列| 午夜影院日韩av| 男插女下体视频免费在线播放| 亚洲av免费高清在线观看| 欧美黑人欧美精品刺激| 午夜老司机福利剧场| 两性午夜刺激爽爽歪歪视频在线观看| 国产一区二区在线av高清观看| 色尼玛亚洲综合影院| 深夜a级毛片| 午夜激情欧美在线| 免费在线观看影片大全网站| 18禁裸乳无遮挡免费网站照片| 久久久久国产精品人妻aⅴ院| 狂野欧美激情性xxxx在线观看| 国产精品亚洲美女久久久| 99久久无色码亚洲精品果冻| 三级国产精品欧美在线观看| 在线看三级毛片| 在线观看一区二区三区| 天堂网av新在线| 嫩草影院精品99| 国产精品国产高清国产av| 一本久久中文字幕| 嫩草影视91久久| 欧美中文日本在线观看视频| 搡女人真爽免费视频火全软件 | 精品欧美国产一区二区三| 久久久午夜欧美精品| 色av中文字幕| 国产男人的电影天堂91| 日韩在线高清观看一区二区三区 | 国产亚洲精品综合一区在线观看| 91在线精品国自产拍蜜月| 波野结衣二区三区在线| 午夜爱爱视频在线播放| 婷婷色综合大香蕉| 日本一本二区三区精品| 18禁在线播放成人免费| 日本成人三级电影网站| 午夜激情欧美在线| 国产高清激情床上av| 人妻丰满熟妇av一区二区三区| av在线天堂中文字幕| 欧美最新免费一区二区三区| 可以在线观看的亚洲视频| 精品人妻视频免费看| 老师上课跳d突然被开到最大视频| 琪琪午夜伦伦电影理论片6080| 国产精品乱码一区二三区的特点| a在线观看视频网站| 国产三级中文精品| 日本三级黄在线观看| 国产黄a三级三级三级人| 欧美绝顶高潮抽搐喷水| 99久久精品国产国产毛片| 国产av在哪里看| 九九爱精品视频在线观看| 天堂√8在线中文| 久久精品国产亚洲av涩爱 | 美女xxoo啪啪120秒动态图| 国产精品1区2区在线观看.| 久久九九热精品免费| 村上凉子中文字幕在线| 久久国产精品人妻蜜桃| 美女xxoo啪啪120秒动态图| 三级男女做爰猛烈吃奶摸视频| 国产美女午夜福利| avwww免费| 亚洲三级黄色毛片| 一区二区三区四区激情视频 | 国产乱人伦免费视频| 99精品久久久久人妻精品| 五月伊人婷婷丁香| 精品国内亚洲2022精品成人| 成人午夜高清在线视频| 国产高潮美女av| 干丝袜人妻中文字幕| 亚洲一级一片aⅴ在线观看| 中国美白少妇内射xxxbb| 久久久久国内视频| av在线观看视频网站免费| 级片在线观看| av在线天堂中文字幕| 亚洲人成网站在线播放欧美日韩| aaaaa片日本免费| 久久久久久久久大av| 国产精品人妻久久久久久| 午夜福利在线观看吧| 此物有八面人人有两片| 精品人妻一区二区三区麻豆 | av视频在线观看入口| 人妻夜夜爽99麻豆av| 亚洲欧美精品综合久久99| 国产亚洲91精品色在线| 日韩精品中文字幕看吧| 乱系列少妇在线播放| 精品久久久久久久久亚洲 | 校园人妻丝袜中文字幕| 欧美日韩黄片免| 性欧美人与动物交配| 国产精品免费一区二区三区在线| 国产成人一区二区在线| 国内精品宾馆在线| 黄片wwwwww| 国产伦精品一区二区三区四那| 99久久无色码亚洲精品果冻| 国产三级中文精品| 中文字幕人妻熟人妻熟丝袜美| av在线天堂中文字幕| 欧美三级亚洲精品| 婷婷丁香在线五月| 国产一区二区亚洲精品在线观看| 黄色日韩在线| 99热这里只有是精品在线观看| 亚州av有码| 人妻制服诱惑在线中文字幕| 久久久色成人| 久久精品久久久久久噜噜老黄 | 成人午夜高清在线视频| 日本免费一区二区三区高清不卡| 久久久久久久久久久丰满 | 欧美最黄视频在线播放免费| 99久久无色码亚洲精品果冻| 久久久午夜欧美精品| 老女人水多毛片| 搡老妇女老女人老熟妇| 国内少妇人妻偷人精品xxx网站| 国产久久久一区二区三区| 久久久精品欧美日韩精品| 国产精品久久久久久久电影| 中出人妻视频一区二区| 国产一区二区激情短视频| 国产白丝娇喘喷水9色精品| 一级黄片播放器| 亚州av有码| 日韩精品中文字幕看吧| 中文字幕精品亚洲无线码一区| 国产精品,欧美在线| 免费看光身美女| 国产美女午夜福利| 国产私拍福利视频在线观看| 不卡视频在线观看欧美| 国产美女午夜福利| 国产91精品成人一区二区三区| 国产探花极品一区二区| 熟女人妻精品中文字幕| 91麻豆精品激情在线观看国产| 熟妇人妻久久中文字幕3abv| 成人av在线播放网站| 他把我摸到了高潮在线观看| 欧美激情国产日韩精品一区| 午夜爱爱视频在线播放| 香蕉av资源在线| 亚洲欧美日韩东京热| 好男人在线观看高清免费视频| 最新中文字幕久久久久| 国产亚洲精品综合一区在线观看| 亚洲四区av| 亚洲美女视频黄频| 欧美+日韩+精品|