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

    坐底式平臺(tái)結(jié)構(gòu)疲勞強(qiáng)度分析

    2016-06-03 05:59:08李亞軍
    造船技術(shù) 2016年2期
    關(guān)鍵詞:波浪工況曲線

    陸 超, 李亞軍

    (中國(guó)艦船研究設(shè)計(jì)中心, 湖北 武漢 430064)

    ?

    坐底式平臺(tái)結(jié)構(gòu)疲勞強(qiáng)度分析

    陸超, 李亞軍

    (中國(guó)艦船研究設(shè)計(jì)中心, 湖北 武漢 430064)

    摘要采用疲勞簡(jiǎn)化分析方法和譜分析方法,使用ANSYS軟件對(duì)某坐底式平臺(tái)進(jìn)行了結(jié)構(gòu)疲勞強(qiáng)度的校核。根據(jù)作業(yè)場(chǎng)地環(huán)境條件進(jìn)行了詳細(xì)的工況劃分,運(yùn)用熱點(diǎn)應(yīng)力的譜分析法計(jì)算對(duì)關(guān)鍵節(jié)點(diǎn)進(jìn)行了細(xì)部疲勞分析。結(jié)果表明,該部分結(jié)構(gòu)設(shè)計(jì)符合疲勞強(qiáng)度要求。

    關(guān)鍵詞坐底式平臺(tái)疲勞強(qiáng)度譜分析法ANSYS

    analysisANSYS

    0引言

    坐底式海洋平臺(tái)由于其設(shè)計(jì)建造周期短、造價(jià)低廉的優(yōu)勢(shì),在河流和海灣等30 m以下的淺水域以及海床平坦的淺海區(qū)域的油氣勘探開發(fā)作業(yè)中發(fā)揮了重要的作用[1~3]。我國(guó)渤海沿岸的勝利油田、大港油田和遼河油田,由于其淺水的海況條件,在開采過程中廣泛使用了坐底式平臺(tái)進(jìn)行作業(yè)[4,5]。

    由于處在各種不利海況和復(fù)雜交變載荷作用下,平臺(tái)結(jié)構(gòu)疲勞破壞是結(jié)構(gòu)失效的重要因素。因此對(duì)平臺(tái)結(jié)構(gòu)的疲勞強(qiáng)度分析是十分必要和關(guān)鍵的。根據(jù)中國(guó)船級(jí)社頒發(fā)的《海上移動(dòng)平臺(tái)入級(jí)規(guī)范》2012版規(guī)定[6]:疲勞校核的目的是確保平臺(tái)結(jié)構(gòu)在營(yíng)運(yùn)期間具有足夠的疲勞壽命。疲勞壽命的計(jì)算結(jié)果可用作制定平臺(tái)在制造和工作期間檢驗(yàn)規(guī)程的依據(jù)。疲勞分析的范圍和方法將取決于平臺(tái)設(shè)計(jì)中所考慮的預(yù)期作業(yè)模式和區(qū)域。結(jié)構(gòu)的疲勞設(shè)計(jì)壽命應(yīng)不小于平臺(tái)的設(shè)計(jì)年限,且不小于20 年。對(duì)于可能產(chǎn)生潛在疲勞裂紋的任一焊縫和引起應(yīng)力集中的結(jié)構(gòu)形式均應(yīng)進(jìn)行抗疲勞設(shè)計(jì)。必要時(shí),應(yīng)進(jìn)行節(jié)點(diǎn)細(xì)部的疲勞分析。對(duì)于柱穩(wěn)式平臺(tái),校核部位主要包括以下內(nèi)容:撐桿;撐桿與下殼體、柱體和甲板相連處;主柱與下殼體相連處;主要結(jié)構(gòu)不連續(xù)處。

    目前海洋平臺(tái)分析校核中常用的平臺(tái)疲勞分析方法主要有疲勞簡(jiǎn)化分析方法、隨機(jī)疲勞分析方法以及確定性方法。疲勞簡(jiǎn)化分析方法,或稱簡(jiǎn)化方法,是用于抗疲勞設(shè)計(jì)的初期階段。在確定危險(xiǎn)結(jié)點(diǎn)后,需再進(jìn)行隨機(jī)疲勞分析。在簡(jiǎn)化方法中對(duì)設(shè)計(jì)參數(shù)的取用應(yīng)適當(dāng)保守。隨機(jī)疲勞分析方法又稱譜分析方法。隨機(jī)疲勞分析應(yīng)用公認(rèn)的方法和準(zhǔn)則進(jìn)行計(jì)算。疲勞分析載荷應(yīng)基于預(yù)期作業(yè)海域或北大西洋長(zhǎng)期分布的環(huán)境數(shù)據(jù)資料。載荷響應(yīng)計(jì)算考慮的浪向角應(yīng)不小于7 個(gè)。確定性方法主要基于經(jīng)驗(yàn)曲線進(jìn)行疲勞壽命估算,精確性不及隨機(jī)疲勞分析方法。本文主要依托ANSYS軟件,采用疲勞簡(jiǎn)化分析方法和隨機(jī)疲勞分析方法對(duì)某坐底式平臺(tái)進(jìn)行疲勞校核及分析。

    1環(huán)境條件

    根據(jù)平臺(tái)所處海區(qū)每月的波高浪向概率分布,可以統(tǒng)計(jì)得出平臺(tái)所處海區(qū)全年的波浪方向概率分布,如表1所示[7]。

    根據(jù)平臺(tái)所處海區(qū)每月的波高和波周期概率分布,可以統(tǒng)計(jì)得出平臺(tái)所處海區(qū)全年的波高和波周期聯(lián)合概率分布,如表2所示。

    表1 平臺(tái)所處海區(qū)全年波浪方向概率分布         單位:%

    表2 平臺(tái)所處海區(qū)全年波高和波周期聯(lián)合概率分布           單位:%

    2計(jì)算模型建立

    2.1平臺(tái)數(shù)據(jù)

    平臺(tái)在空載情況下的重量控制如表3所示。

    表3 平臺(tái)在空載情況下的重量控制表

    2.2計(jì)算模型及載荷輸入

    本項(xiàng)目中采用ANSYS軟件對(duì)坐底式平臺(tái)的結(jié)構(gòu)疲勞損傷進(jìn)行計(jì)算。根據(jù)原型數(shù)據(jù)建立有限元模型,如圖1所示。其中,船體、沉墊的外殼模型采用Shell281單元建立。Shell281適用于薄到中等厚度的殼結(jié)構(gòu),該單元有8個(gè)節(jié)點(diǎn),每個(gè)節(jié)點(diǎn)有六個(gè)自由度,分別為沿節(jié)點(diǎn)X、Y、Z方向的平動(dòng)及繞節(jié)點(diǎn)X、Y、Z軸的轉(zhuǎn)動(dòng)。Shell281單元具有應(yīng)力剛化及大變形功能,且有強(qiáng)大的非線性功能,可以對(duì)截面數(shù)據(jù)進(jìn)行定義、分析、可視化等操作。由于Shell281單元具備以上特性,因此采用該單元計(jì)算平臺(tái)立柱與下殼體連接處的熱點(diǎn)應(yīng)力是有效的。對(duì)于所有的骨材和肋板,采用Beam188單元建立。

    圖1 平臺(tái)的有限元模型

    平臺(tái)上部模塊由于并未在疲勞分析中扮演重要角色,因此在模型中并未體現(xiàn),但各部分功能模塊的具體重量則按照重量控制報(bào)告分別在上部甲板的相應(yīng)位置處予以了模擬,模擬方法采用調(diào)整等效材料密度的方式實(shí)現(xiàn)。

    根據(jù)波浪理論可知,當(dāng)結(jié)構(gòu)的直徑與波長(zhǎng)比值小于2.0時(shí)采用Morrison方法計(jì)算波浪載荷,據(jù)此對(duì)于該平臺(tái)立柱的載荷施加可采用基于Morrison理論的水動(dòng)力計(jì)算單元。為便于在計(jì)算模型中形成有效的載荷輸入,在所有樁腿模型的中心部位額外設(shè)置了pipe59單元。該類單元是ANSYS軟件中特有的水動(dòng)力計(jì)算單元,通過在單元屬性設(shè)置watertable項(xiàng)即可形成對(duì)波浪參數(shù)的輸入。pipe59單元在模型中設(shè)置為所有節(jié)點(diǎn)與相應(yīng)位置處的Shell281單元所有節(jié)點(diǎn)耦合,已形成有效的載荷傳遞操作,并保證不會(huì)對(duì)結(jié)構(gòu)整體剛度形成影響。

    2.3疲勞損傷評(píng)估位置

    為確定平臺(tái)的疲勞損傷評(píng)估位置,首先針對(duì)平臺(tái)在極端波浪工況條件下的應(yīng)力分布情況進(jìn)行了計(jì)算分析。計(jì)算結(jié)果顯示,平臺(tái)六根立柱與上部沉墊結(jié)合處的應(yīng)力水平最高,并且應(yīng)力在這些位置處表現(xiàn)出很高的應(yīng)力集中水平,如圖2、圖3所示。在對(duì)結(jié)構(gòu)主要連接處的熱點(diǎn)應(yīng)力進(jìn)行有限元分析時(shí),其主要側(cè)重點(diǎn)在于生成一個(gè)模型,并使該模型給出受焊縫影響的外側(cè)區(qū)域應(yīng)力精度。此模型應(yīng)具備將應(yīng)力外推回到焊趾的精細(xì)網(wǎng)格,以保證熱點(diǎn)應(yīng)力計(jì)算的精度。應(yīng)力集中分析的有限元模型對(duì)單元類型和網(wǎng)格尺寸通常較為敏感。在減小單元尺寸時(shí),有限元應(yīng)力在不連續(xù)點(diǎn)處會(huì)趨于無限大。因此,有必要給單元尺寸設(shè)定一個(gè)下限并使用外推到熱點(diǎn)的方法,對(duì)不同計(jì)算機(jī)程序和用戶的結(jié)果進(jìn)行對(duì)比,有一個(gè)較為統(tǒng)一的基礎(chǔ)。另一方面,為了正確地得到幾何應(yīng)力增加的情況,使位于“t/2”和“3t/2”(t為板厚)處的應(yīng)力參考點(diǎn)不處于同一單元內(nèi)至關(guān)重要。這意味著板厚數(shù)量級(jí)的單元尺寸可用于建模。在焊縫的整個(gè)幾何形狀都包括在模型描述中的情況下,如果使用實(shí)體單元模擬,位于熱點(diǎn)附近的單元尺寸必須要減少到板厚的一半[8]。

    圖2 平臺(tái)結(jié)構(gòu)整體應(yīng)力分布

    圖3 平臺(tái)結(jié)構(gòu)熱點(diǎn)應(yīng)力分析區(qū)

    單元應(yīng)力通常是高斯積分點(diǎn)推算出來的。單元類型,做一些外推以便確定代表焊趾位置的應(yīng)力是必要的。為了保存熱點(diǎn)處主應(yīng)力方向的資料,應(yīng)使用分量應(yīng)力以外推。如果用殼單元模擬,且整個(gè)焊縫的幾何形狀不包括在模型中,則必須外推到單元交界線上。如果整個(gè)焊縫幾何形狀包括在模型中(3D模型),則外推與焊趾有關(guān)。若使用8節(jié)點(diǎn)殼單元,則熱點(diǎn)可認(rèn)為位于單元交界線上。示例如圖4所示。

    圖4 單元模擬示例

    據(jù)此,在本項(xiàng)目有限元模型的建立中,對(duì)平臺(tái)的立柱與下殼體連接處的局部區(qū)域進(jìn)行了詳細(xì)的劃分,如圖5所示。由于在不同的外部載荷作用方向下,連接處焊縫周圍的應(yīng)力分布將會(huì)發(fā)生變化,因此選取平臺(tái)立柱與下殼體6個(gè)連接處每處的外側(cè)8個(gè)點(diǎn)的熱點(diǎn)應(yīng)力進(jìn)行計(jì)算。圖6為8個(gè)熱點(diǎn)的位置示意圖。

    圖5 平臺(tái)立柱與殼體連接處的局部區(qū)域網(wǎng)格劃分

    根據(jù)上述分析結(jié)果,可將平臺(tái)結(jié)構(gòu)最終的疲勞分析對(duì)象確定為:平臺(tái)六根立柱與上部沉墊結(jié)合處。

    3簡(jiǎn)化疲勞分析

    根據(jù)DNVGL-RP-C203[9]的要求,當(dāng)采用簡(jiǎn)化

    圖6 連接處8個(gè)熱點(diǎn)的位置示意圖

    疲勞分析時(shí),平臺(tái)結(jié)構(gòu)應(yīng)滿足以下疲勞損傷要求:

    (1)

    DNVGL-RP-C203規(guī)范根據(jù)不同管節(jié)點(diǎn)類型推薦了兩種S-N曲線,即WJ曲線和CJ曲線。其中,WJ曲線適用于焊接管節(jié)點(diǎn),而CJ曲線主要用于澆鑄類型的管節(jié)點(diǎn)。對(duì)于這兩種不同類型的管節(jié)點(diǎn),S-N曲線表達(dá)式中的各參數(shù)取值不同,取值結(jié)果如表4所示。管節(jié)點(diǎn)基礎(chǔ)設(shè)計(jì)S-N曲線可以寫成下式形式:

    (2)

    式中:N為在應(yīng)力范圍S下預(yù)計(jì)發(fā)生破壞時(shí)的應(yīng)力循環(huán)次數(shù);k1為一常數(shù);m為S-N曲線的反斜率。

    表4 基礎(chǔ)設(shè)計(jì)S-N曲線參數(shù)取值

    結(jié)合本項(xiàng)目平臺(tái)結(jié)構(gòu)型式,選用WJ曲線進(jìn)行疲勞分析。WJ曲線形式如圖7所示。需要注意的是,這里的S-N曲線是在標(biāo)準(zhǔn)厚度(16 mm)下得到的。因此,在疲勞分析過程中還需要考慮實(shí)際管件厚度對(duì)應(yīng)力范圍的影響,并進(jìn)行厚度修正。當(dāng)材料厚度大于標(biāo)準(zhǔn)厚度時(shí),可以按照下式進(jìn)行修正:

    圖7 標(biāo)準(zhǔn)厚度(16mm)下S-N曲線

    (3)

    式中:tref為標(biāo)準(zhǔn)厚度(16 mm);S為實(shí)際應(yīng)力范圍;S0為S-N曲線中標(biāo)準(zhǔn)厚度下的應(yīng)力范圍;t是材料厚度。依據(jù)上述S-N曲線計(jì)算得到的累積疲勞損傷還應(yīng)乘以5.0的安全系數(shù)。

    根據(jù)平臺(tái)所處海區(qū)每月的波高浪向概率分布,可以統(tǒng)計(jì)得出平臺(tái)所處海區(qū)全年的平均跨零周期為5.2 s,對(duì)應(yīng)的波高為2.0 m。對(duì)上述的平臺(tái)模型進(jìn)行該波浪條件下的瞬態(tài)動(dòng)力分析,可得六個(gè)計(jì)算位置處的應(yīng)力分布情況(見圖8)和波動(dòng)情況(18.5 MPa)。

    圖8 簡(jiǎn)化疲勞分析熱點(diǎn)應(yīng)力分布

    根據(jù)以上結(jié)果可確定WJ曲線的參數(shù)m為3,Gama函數(shù)取值根據(jù)DNVGL-RP-C203要求確定為6.0,將18.5 MPa的應(yīng)力波動(dòng)范圍帶入上述簡(jiǎn)化疲勞校核公式(式1)中,可得到簡(jiǎn)化評(píng)估的總損傷比為:D=0.66<1.0。據(jù)此,簡(jiǎn)化疲勞分析結(jié)果表明,平臺(tái)結(jié)構(gòu)的疲勞壽命滿足設(shè)計(jì)要求。

    4波浪譜疲勞分析

    4.1計(jì)算工況

    根據(jù)第3章中給出的不同波高、波周期以及波浪方向設(shè)定不同的疲勞計(jì)算工況。共計(jì)算8種不同的波浪方向(即N、EN、E、ES、S、WS、W、WN分別表示北、東北、東、東南、南、西南、西以及西北方位),每個(gè)方向下包含40種不同的波高波周期組合,總計(jì)320個(gè)計(jì)算工況。平臺(tái)的坐向及6根立柱的編號(hào)如圖9所示。

    圖9 平臺(tái)坐向及立柱編號(hào)示意圖

    4.2波浪譜疲勞分析

    在各波浪疲勞分析工況中,選取Jonswap譜描述波浪能在頻域的分布情況。

    由上一節(jié)得到的波高譜通過傳遞函數(shù)轉(zhuǎn)換,可以得到結(jié)構(gòu)各連接處的應(yīng)力波動(dòng)范圍響應(yīng)譜:

    (4)

    式中:Sηη(f)是波高功率密度函數(shù);T(f)是應(yīng)力范圍傳遞函數(shù)。

    為確定結(jié)構(gòu)各連接處應(yīng)力范圍響應(yīng)譜與波高譜之間的傳遞函數(shù)T(f),使用ANSYS有限元軟件對(duì)坐底式平臺(tái)結(jié)構(gòu)進(jìn)行瞬態(tài)動(dòng)力分析,進(jìn)而得到結(jié)構(gòu)各連接處的熱點(diǎn)應(yīng)力范圍。

    對(duì)于上一節(jié)得到的應(yīng)力范圍響應(yīng)譜,其n階譜矩:

    (5)

    應(yīng)力范圍標(biāo)準(zhǔn)差:

    (6)

    應(yīng)力響應(yīng)有效頻率:

    (7)

    結(jié)構(gòu)連接處的應(yīng)力范圍歷程可假定是一個(gè)窄帶高斯過程,其峰值服從瑞利分布。由式(6)確定的應(yīng)力范圍標(biāo)準(zhǔn)差可得瑞利分布概率密度函數(shù)p(σr):

    (8)

    該函數(shù)曲線所圍總面積為1。其中,對(duì)于大于5倍標(biāo)準(zhǔn)差的應(yīng)力范圍,其出現(xiàn)概率可以忽略。對(duì)于5倍標(biāo)準(zhǔn)差以內(nèi)的應(yīng)力范圍,可將其劃分為更小的應(yīng)力范圍區(qū)域,用各分塊應(yīng)力范圍均值代表該區(qū)域內(nèi)的應(yīng)力范圍值。各區(qū)域的出現(xiàn)概率與曲線下所圍頻帶面積相對(duì)應(yīng),第i個(gè)應(yīng)力范圍區(qū)域的出現(xiàn)概率為

    (9)

    在波浪疲勞分析工況j下,第i個(gè)應(yīng)力范圍區(qū)域的年循環(huán)次數(shù)為

    (10)

    式中:右邊第一個(gè)常數(shù)是一年中作業(yè)時(shí)間,單位為s;P%是疲勞工況j的年出現(xiàn)概率;P(σri)是第i個(gè)應(yīng)力范圍區(qū)域的出現(xiàn)概率,fe是疲勞工況j的應(yīng)力響應(yīng)有效頻率。

    將各波浪疲勞分析工況下各應(yīng)力范圍區(qū)域的熱點(diǎn)應(yīng)力范圍和年循環(huán)次數(shù)代入S-N曲線中,即可求得該應(yīng)力范圍區(qū)域的年疲勞損傷。在波浪疲勞分析中選用API RP WSD[10]規(guī)范推薦的管節(jié)點(diǎn)S-N曲線。

    由各波浪疲勞分析工況下,各應(yīng)力范圍區(qū)域的年循環(huán)次數(shù)和通過S-N曲線計(jì)算得到的臨界循環(huán)次數(shù),可求得結(jié)構(gòu)各連接處的年累積疲勞損傷Df:

    (11)

    式中:nij是疲勞分析工況j下第i個(gè)應(yīng)力范圍區(qū)域的年循環(huán)次數(shù);Nij是由S-N曲線計(jì)算得到的該應(yīng)力范圍的臨界循環(huán)次數(shù)。將兩者在不同疲勞工況和不同應(yīng)力范圍下的比值相加,就得到年累積疲勞損傷Df。

    在周期性的波浪載荷作用下,平臺(tái)立柱與下殼體各連接處的疲勞壽命Tf即為上述年累積疲勞損傷的倒數(shù):

    (12)

    4.3校核結(jié)果

    坐底式平臺(tái)結(jié)構(gòu)的波浪譜疲勞分析結(jié)果如表5所示。可以看出,平臺(tái)的最大累計(jì)疲勞損傷出現(xiàn)在3號(hào)和4號(hào)立柱與沉墊的結(jié)合位置,最小疲勞壽命分別為95年和96年。

    表5 波浪譜疲勞分析結(jié)果

    續(xù)表5 波浪譜疲勞分析結(jié)果

    5總結(jié)

    經(jīng)過簡(jiǎn)化疲勞分析的總體估計(jì),以及波浪譜疲勞分析的多工況熱點(diǎn)應(yīng)力計(jì)算,選取的6個(gè)立柱與沉墊結(jié)合處的熱點(diǎn)應(yīng)力疲勞壽命滿足平臺(tái)設(shè)計(jì)壽命要求,該部分的結(jié)構(gòu)設(shè)計(jì)對(duì)于疲勞強(qiáng)度是合理的。

    參考文獻(xiàn)

    [1]任貴永. 坐底式平臺(tái)在淺海石油勘探開發(fā)中的應(yīng)用[J]. 中國(guó)海洋平臺(tái), 1992(3):126-128.

    [2]任貴永, 孟昭瑛. 坐底式平臺(tái)特性及其在淺海開發(fā)中的應(yīng)用[J]. 中國(guó)海洋平臺(tái), 1994(3):32-35.

    [3]宋長(zhǎng)智, 張格, 唐廣銀. 老齡坐底式平臺(tái)整體強(qiáng)度分析研究[J]. 中國(guó)海洋平臺(tái), 2015(1):27-30.

    [4]田海慶. 淺海水域坐底式平臺(tái)安全定位措施與監(jiān)測(cè)技術(shù)[J]. 中國(guó)海洋平臺(tái), 2007(4):37-38.

    [5]汪張?zhí)? “中油海3號(hào)”坐底式鉆井平臺(tái)[J]. 上海造船, 2009,78(2):53-54.

    [6]中國(guó)船級(jí)社.海上移動(dòng)平臺(tái)入級(jí)規(guī)范[M]. 北京:人民交通出版社,2012.

    [7]竺艷蓉.海洋工程波浪力學(xué)[M].天津:天津大學(xué)出版社,1991.

    [8]張朝陽, 劉俊, 白艷彬,等. 基于譜分析法的深水半潛式平臺(tái)疲勞強(qiáng)度分析[J]. 海洋工程, 2012, 30(1):53-59.

    [9]DNV. Fatigue Design of Offshore Steel Structures[S].2005.

    [10]American Petroleum Institute.RP2A-WSD:Planning,Designing,and Constructing Fixed Offshore Platforms-Working Stress Design[S].Houston:American Petroleum Institute,2002.

    Analysis on Fatigue Strength of Bottom Supported Platform

    LU Chao, LI Ya-jun

    (China Ship Research and Design Center, Wuhan Hubei 430064, China)

    AbstractThe check of structural fatigue strength of bottom supported platform is accomplished based on simplified fatigue analysis method and spectral-based fatigue analysis method with the help of ANSYS. Detailed working conditions are arranged in accordance with local environmental conditions. Fatigue analyses of key joints are carried on with the use of spectral-based fatigue analysis for hot spots. The results show that the fatigue life of the structure satisfies the design strength request.

    KeywordsBottom supported platformFatigue strengthSpectral-based fatigue

    中圖分類號(hào)P75

    文獻(xiàn)標(biāo)志碼A

    作者簡(jiǎn)介:陸超(1981-),男,碩士,工程師,研究方向?yàn)楹I辖Y(jié)構(gòu)物設(shè)計(jì)評(píng)估。

    猜你喜歡
    波浪工況曲線
    未來訪談:出版的第二增長(zhǎng)曲線在哪里?
    出版人(2022年8期)2022-08-23 03:36:50
    波浪谷和波浪巖
    熱網(wǎng)異常工況的辨識(shí)
    煤氣與熱力(2022年4期)2022-05-23 12:44:44
    不同工況下噴水推進(jìn)泵內(nèi)流性能研究
    幸福曲線
    英語文摘(2020年6期)2020-09-21 09:30:40
    基于非負(fù)矩陣分解的高速列車走行部工況識(shí)別
    沿平坦凸曲線Hilbert變換的L2有界性
    波浪谷隨想
    去看神奇波浪谷
    夢(mèng)寐以求的S曲線
    Coco薇(2015年10期)2015-10-19 12:42:05
    av福利片在线| 两个人视频免费观看高清| 中国美女看黄片| 1024手机看黄色片| 国产片内射在线| 757午夜福利合集在线观看| 免费在线观看影片大全网站| 亚洲最大成人中文| 亚洲欧美精品综合一区二区三区| 一区二区三区激情视频| 免费无遮挡裸体视频| 欧美日韩亚洲综合一区二区三区_| 久久久久亚洲av毛片大全| 国产精品亚洲一级av第二区| 国产午夜精品论理片| 看黄色毛片网站| 亚洲中文日韩欧美视频| 一本综合久久免费| 日韩精品中文字幕看吧| 色哟哟哟哟哟哟| 欧美日本亚洲视频在线播放| 精华霜和精华液先用哪个| 香蕉av资源在线| 国产精品日韩av在线免费观看| 一夜夜www| 免费看十八禁软件| 精品一区二区三区四区五区乱码| 老司机福利观看| 亚洲欧洲精品一区二区精品久久久| 中文字幕av在线有码专区| 成人av一区二区三区在线看| 琪琪午夜伦伦电影理论片6080| 在线观看免费日韩欧美大片| 啦啦啦观看免费观看视频高清| 亚洲熟女毛片儿| 婷婷精品国产亚洲av在线| 日韩成人在线观看一区二区三区| 久久久久久久久中文| 欧美高清成人免费视频www| 亚洲精品在线美女| 美女午夜性视频免费| 久久香蕉国产精品| 嫁个100分男人电影在线观看| 99久久国产精品久久久| 天堂影院成人在线观看| 亚洲一区二区三区不卡视频| 老汉色∧v一级毛片| 国产一区二区在线观看日韩 | 日本a在线网址| 可以在线观看的亚洲视频| 99久久无色码亚洲精品果冻| 午夜福利高清视频| 午夜两性在线视频| 日韩大尺度精品在线看网址| 天堂av国产一区二区熟女人妻 | 欧美绝顶高潮抽搐喷水| 国产一区二区在线观看日韩 | 国产一区在线观看成人免费| 精品久久久久久久末码| 男女那种视频在线观看| 国内少妇人妻偷人精品xxx网站 | 国产精品亚洲av一区麻豆| 国产熟女xx| 最好的美女福利视频网| 嫩草影院精品99| 午夜激情福利司机影院| netflix在线观看网站| 婷婷精品国产亚洲av在线| 最近视频中文字幕2019在线8| 无人区码免费观看不卡| 国产精品亚洲一级av第二区| 男女床上黄色一级片免费看| 给我免费播放毛片高清在线观看| 国产1区2区3区精品| 日韩免费av在线播放| 久久香蕉国产精品| 国产v大片淫在线免费观看| 天天一区二区日本电影三级| 婷婷六月久久综合丁香| 亚洲成av人片免费观看| 久久 成人 亚洲| 亚洲全国av大片| 琪琪午夜伦伦电影理论片6080| 精品无人区乱码1区二区| 久久久久久久久中文| av天堂在线播放| 丝袜美腿诱惑在线| 亚洲avbb在线观看| 欧美激情久久久久久爽电影| 99热只有精品国产| 免费在线观看亚洲国产| 亚洲国产日韩欧美精品在线观看 | 国内少妇人妻偷人精品xxx网站 | 欧美另类亚洲清纯唯美| 久久香蕉精品热| 亚洲18禁久久av| 欧美日韩一级在线毛片| 天天一区二区日本电影三级| 美女高潮喷水抽搐中文字幕| 狂野欧美白嫩少妇大欣赏| 午夜日韩欧美国产| 中文字幕人成人乱码亚洲影| 丝袜人妻中文字幕| 怎么达到女性高潮| 久久国产乱子伦精品免费另类| 51午夜福利影视在线观看| 人妻久久中文字幕网| 欧美不卡视频在线免费观看 | 一区二区三区国产精品乱码| 精华霜和精华液先用哪个| 一个人观看的视频www高清免费观看 | 久久精品国产综合久久久| 日韩精品青青久久久久久| 中文字幕熟女人妻在线| 91麻豆av在线| 精品国产超薄肉色丝袜足j| 99在线人妻在线中文字幕| 一进一出抽搐动态| 两性午夜刺激爽爽歪歪视频在线观看 | 大型av网站在线播放| 国产成人啪精品午夜网站| 国产视频一区二区在线看| 亚洲激情在线av| 俺也久久电影网| 国产精品电影一区二区三区| 一级a爱片免费观看的视频| 亚洲五月天丁香| 国产又黄又爽又无遮挡在线| 每晚都被弄得嗷嗷叫到高潮| 午夜亚洲福利在线播放| 黄色丝袜av网址大全| 欧美日韩黄片免| 琪琪午夜伦伦电影理论片6080| 国产69精品久久久久777片 | 亚洲精品美女久久久久99蜜臀| 别揉我奶头~嗯~啊~动态视频| 人人妻人人看人人澡| 日韩欧美精品v在线| 欧美在线一区亚洲| 淫妇啪啪啪对白视频| 又粗又爽又猛毛片免费看| 丝袜美腿诱惑在线| 久久久久久久久免费视频了| 人人妻人人澡欧美一区二区| 18禁裸乳无遮挡免费网站照片| 高潮久久久久久久久久久不卡| 伊人久久大香线蕉亚洲五| 亚洲欧美一区二区三区黑人| 午夜久久久久精精品| 免费在线观看影片大全网站| 最近最新中文字幕大全免费视频| 中文字幕久久专区| 男女下面进入的视频免费午夜| 桃红色精品国产亚洲av| 久久久久免费精品人妻一区二区| 丰满人妻一区二区三区视频av | 嫩草影院精品99| 欧美日韩亚洲国产一区二区在线观看| 亚洲全国av大片| 国产高清有码在线观看视频 | 少妇粗大呻吟视频| 亚洲 国产 在线| 国产午夜精品论理片| 999精品在线视频| 国产主播在线观看一区二区| 床上黄色一级片| 亚洲精品中文字幕一二三四区| 日韩有码中文字幕| 伦理电影免费视频| 精品福利观看| 国产一区在线观看成人免费| 国产黄色小视频在线观看| 狂野欧美白嫩少妇大欣赏| 看免费av毛片| 精品国产亚洲在线| 亚洲激情在线av| 这个男人来自地球电影免费观看| 校园春色视频在线观看| 国产精品亚洲美女久久久| 久久精品人妻少妇| 久久香蕉国产精品| 国产熟女午夜一区二区三区| 国产探花在线观看一区二区| 亚洲熟女毛片儿| 国产精品爽爽va在线观看网站| 欧洲精品卡2卡3卡4卡5卡区| 欧美日本亚洲视频在线播放| 麻豆成人午夜福利视频| 欧美 亚洲 国产 日韩一| 他把我摸到了高潮在线观看| 99热6这里只有精品| 757午夜福利合集在线观看| 精品乱码久久久久久99久播| 亚洲avbb在线观看| 国产男靠女视频免费网站| 欧美成人一区二区免费高清观看 | 亚洲va日本ⅴa欧美va伊人久久| 少妇被粗大的猛进出69影院| 国产真实乱freesex| 国产免费男女视频| 女人被狂操c到高潮| 亚洲国产精品合色在线| 露出奶头的视频| 欧洲精品卡2卡3卡4卡5卡区| 亚洲自拍偷在线| 久久天躁狠狠躁夜夜2o2o| 男男h啪啪无遮挡| 久久亚洲精品不卡| 国产视频内射| 黄色 视频免费看| 国产黄色小视频在线观看| 夜夜躁狠狠躁天天躁| 神马国产精品三级电影在线观看 | 一进一出抽搐gif免费好疼| 黄频高清免费视频| 亚洲av中文字字幕乱码综合| 久久天堂一区二区三区四区| 两人在一起打扑克的视频| 日韩中文字幕欧美一区二区| 男女之事视频高清在线观看| 高潮久久久久久久久久久不卡| 美女免费视频网站| 国产熟女午夜一区二区三区| 亚洲va日本ⅴa欧美va伊人久久| 国产精品免费一区二区三区在线| 精品午夜福利视频在线观看一区| 18禁观看日本| tocl精华| 精品不卡国产一区二区三区| 我要搜黄色片| 久久久久久久久免费视频了| 五月伊人婷婷丁香| 麻豆成人午夜福利视频| 在线a可以看的网站| 叶爱在线成人免费视频播放| 国产探花在线观看一区二区| 国产主播在线观看一区二区| 亚洲人成电影免费在线| 日本黄色视频三级网站网址| 国产亚洲精品久久久久5区| 怎么达到女性高潮| 亚洲,欧美精品.| 亚洲第一欧美日韩一区二区三区| 午夜免费成人在线视频| 成人av在线播放网站| 欧美精品亚洲一区二区| 亚洲av第一区精品v没综合| 脱女人内裤的视频| 88av欧美| 国产精品久久电影中文字幕| 搡老岳熟女国产| 18禁黄网站禁片午夜丰满| 国产精品免费视频内射| 成在线人永久免费视频| 精品高清国产在线一区| 中文字幕人妻丝袜一区二区| 国产片内射在线| 少妇人妻一区二区三区视频| 中文字幕精品亚洲无线码一区| 久热爱精品视频在线9| 我要搜黄色片| 在线观看舔阴道视频| 亚洲熟妇中文字幕五十中出| 国产人伦9x9x在线观看| 丁香欧美五月| 欧美黄色淫秽网站| 亚洲真实伦在线观看| 国产亚洲av嫩草精品影院| 婷婷亚洲欧美| 亚洲欧美一区二区三区黑人| 久久精品国产综合久久久| 国产高清视频在线观看网站| 亚洲欧洲精品一区二区精品久久久| 少妇人妻一区二区三区视频| 久久久久久亚洲精品国产蜜桃av| 国产av在哪里看| 欧美色视频一区免费| 亚洲色图av天堂| 黄色视频,在线免费观看| 久久久久久久久免费视频了| 免费电影在线观看免费观看| 在线播放国产精品三级| 婷婷六月久久综合丁香| 蜜桃久久精品国产亚洲av| 久久人人精品亚洲av| 国产真实乱freesex| 日本成人三级电影网站| 久久亚洲真实| 欧美日韩精品网址| cao死你这个sao货| 又黄又爽又免费观看的视频| 性色av乱码一区二区三区2| 中文字幕av在线有码专区| 天堂动漫精品| 夜夜躁狠狠躁天天躁| 国产免费av片在线观看野外av| 久久草成人影院| 97碰自拍视频| 日本一区二区免费在线视频| 动漫黄色视频在线观看| 国产精华一区二区三区| 国产视频一区二区在线看| 日韩av在线大香蕉| 成年版毛片免费区| 母亲3免费完整高清在线观看| 一本大道久久a久久精品| 久久精品综合一区二区三区| 啦啦啦韩国在线观看视频| 国产精品免费一区二区三区在线| 亚洲天堂国产精品一区在线| 国产私拍福利视频在线观看| 天堂影院成人在线观看| 国产精品一及| 两个人看的免费小视频| 国产精品一及| 九色国产91popny在线| 国产免费av片在线观看野外av| 老司机福利观看| 亚洲欧美日韩无卡精品| 夜夜夜夜夜久久久久| 91字幕亚洲| 午夜精品久久久久久毛片777| 亚洲成人免费电影在线观看| 一边摸一边抽搐一进一小说| 久久久久国产精品人妻aⅴ院| 免费观看人在逋| 欧美日本视频| 听说在线观看完整版免费高清| 成人欧美大片| 女生性感内裤真人,穿戴方法视频| 精品久久久久久久久久免费视频| 精品乱码久久久久久99久播| 亚洲av成人精品一区久久| 久久99热这里只有精品18| 99久久综合精品五月天人人| 在线观看一区二区三区| 法律面前人人平等表现在哪些方面| 国产av在哪里看| 久久九九热精品免费| 久久久久免费精品人妻一区二区| 中文字幕人妻丝袜一区二区| 亚洲欧美日韩高清在线视频| 老汉色av国产亚洲站长工具| 看免费av毛片| 久久精品综合一区二区三区| 国产亚洲精品综合一区在线观看 | 亚洲欧美日韩东京热| 日韩欧美免费精品| 成人永久免费在线观看视频| 18禁美女被吸乳视频| 国产成+人综合+亚洲专区| 露出奶头的视频| 波多野结衣巨乳人妻| 国产午夜精品久久久久久| 国内精品久久久久精免费| 亚洲18禁久久av| 精品一区二区三区av网在线观看| 天堂√8在线中文| 黑人欧美特级aaaaaa片| 五月玫瑰六月丁香| 啦啦啦观看免费观看视频高清| 天堂动漫精品| 国产激情偷乱视频一区二区| 嫩草影院精品99| 日韩欧美在线乱码| 日韩欧美免费精品| 色老头精品视频在线观看| 国产高清有码在线观看视频 | 成在线人永久免费视频| 搞女人的毛片| 岛国在线免费视频观看| 老司机深夜福利视频在线观看| 欧美日韩乱码在线| 国产成人精品久久二区二区免费| 99热只有精品国产| 亚洲av成人精品一区久久| 欧美黑人欧美精品刺激| 一级毛片女人18水好多| 国产午夜精品久久久久久| 免费一级毛片在线播放高清视频| 90打野战视频偷拍视频| 国产伦一二天堂av在线观看| 欧美成人性av电影在线观看| 亚洲全国av大片| 久久亚洲真实| 日本熟妇午夜| 欧美性猛交╳xxx乱大交人| 久久伊人香网站| 99国产极品粉嫩在线观看| 老司机午夜十八禁免费视频| 国产av在哪里看| 国产精品亚洲av一区麻豆| 丁香欧美五月| 国产精品 欧美亚洲| av在线播放免费不卡| 日本在线视频免费播放| 岛国视频午夜一区免费看| 精品一区二区三区av网在线观看| 在线免费观看的www视频| 夜夜躁狠狠躁天天躁| 久久精品影院6| 免费一级毛片在线播放高清视频| av在线天堂中文字幕| 久久久久免费精品人妻一区二区| 亚洲一码二码三码区别大吗| 国产亚洲av高清不卡| www.www免费av| 小说图片视频综合网站| 亚洲色图 男人天堂 中文字幕| 999精品在线视频| 可以在线观看的亚洲视频| 欧美又色又爽又黄视频| 国产视频内射| 午夜a级毛片| 国产人伦9x9x在线观看| 国产精品一区二区三区四区久久| 中文字幕av在线有码专区| 中文字幕熟女人妻在线| 又黄又爽又免费观看的视频| a在线观看视频网站| 欧美成人午夜精品| av超薄肉色丝袜交足视频| 9191精品国产免费久久| 在线观看一区二区三区| 1024香蕉在线观看| 男人的好看免费观看在线视频 | 国产91精品成人一区二区三区| 久久久久久亚洲精品国产蜜桃av| 国产精品影院久久| 99热6这里只有精品| 午夜福利高清视频| 久久人人精品亚洲av| 亚洲黑人精品在线| 婷婷亚洲欧美| 两个人视频免费观看高清| 又黄又粗又硬又大视频| 老汉色av国产亚洲站长工具| 一本精品99久久精品77| 天天躁狠狠躁夜夜躁狠狠躁| 免费在线观看黄色视频的| 在线观看www视频免费| 久久天躁狠狠躁夜夜2o2o| 又黄又粗又硬又大视频| 婷婷六月久久综合丁香| 制服诱惑二区| 一二三四在线观看免费中文在| 琪琪午夜伦伦电影理论片6080| 制服人妻中文乱码| 99re在线观看精品视频| 无人区码免费观看不卡| 亚洲中文字幕日韩| 少妇裸体淫交视频免费看高清 | 动漫黄色视频在线观看| 国产伦人伦偷精品视频| 看片在线看免费视频| 久久精品国产亚洲av香蕉五月| 很黄的视频免费| 国产一区二区三区视频了| 久久中文字幕一级| 男女那种视频在线观看| 18禁黄网站禁片免费观看直播| 精品熟女少妇八av免费久了| 激情在线观看视频在线高清| 视频区欧美日本亚洲| 看免费av毛片| 欧美日韩乱码在线| 国产视频一区二区在线看| 欧美日韩一级在线毛片| 国产成人系列免费观看| 九色国产91popny在线| 久久精品综合一区二区三区| 搞女人的毛片| 男插女下体视频免费在线播放| 日韩免费av在线播放| 欧美中文日本在线观看视频| 免费av毛片视频| 日日夜夜操网爽| 一个人免费在线观看电影 | 国产欧美日韩一区二区三| 精品国内亚洲2022精品成人| 在线观看美女被高潮喷水网站 | 欧美成人性av电影在线观看| 全区人妻精品视频| 在线观看免费午夜福利视频| 搡老岳熟女国产| 观看免费一级毛片| 最近视频中文字幕2019在线8| 色老头精品视频在线观看| 久久九九热精品免费| 99国产极品粉嫩在线观看| 99热这里只有精品一区 | 舔av片在线| 757午夜福利合集在线观看| 禁无遮挡网站| 欧美日本亚洲视频在线播放| 中文在线观看免费www的网站 | 亚洲色图av天堂| 国产av又大| 精品人妻1区二区| 欧美中文综合在线视频| 变态另类丝袜制服| 欧美精品啪啪一区二区三区| 亚洲色图 男人天堂 中文字幕| 久久人人精品亚洲av| 在线播放国产精品三级| 露出奶头的视频| 老司机在亚洲福利影院| 国产成人aa在线观看| 男人舔女人下体高潮全视频| 国产探花在线观看一区二区| 亚洲七黄色美女视频| 亚洲中文字幕一区二区三区有码在线看 | 亚洲国产欧美一区二区综合| 夜夜看夜夜爽夜夜摸| 老司机午夜十八禁免费视频| 一进一出好大好爽视频| 欧美最黄视频在线播放免费| 国产成人精品久久二区二区91| 黄色成人免费大全| 99久久国产精品久久久| 午夜a级毛片| 精品日产1卡2卡| 美女午夜性视频免费| 国产伦人伦偷精品视频| aaaaa片日本免费| 每晚都被弄得嗷嗷叫到高潮| 一本久久中文字幕| 老司机深夜福利视频在线观看| 久久久国产精品麻豆| 久久 成人 亚洲| 久久久水蜜桃国产精品网| 99热这里只有精品一区 | 日韩欧美一区二区三区在线观看| 亚洲精品av麻豆狂野| 亚洲黑人精品在线| 成人国产一区最新在线观看| 99久久综合精品五月天人人| 亚洲精品美女久久久久99蜜臀| 老司机深夜福利视频在线观看| 香蕉国产在线看| 精华霜和精华液先用哪个| 人妻丰满熟妇av一区二区三区| 久久精品国产亚洲av高清一级| 身体一侧抽搐| 人妻久久中文字幕网| 19禁男女啪啪无遮挡网站| 激情在线观看视频在线高清| 狂野欧美白嫩少妇大欣赏| avwww免费| 国产精品野战在线观看| 国产精品久久久久久人妻精品电影| 亚洲一区二区三区色噜噜| 精品久久蜜臀av无| 国产欧美日韩精品亚洲av| 欧美精品啪啪一区二区三区| 日韩有码中文字幕| 国产蜜桃级精品一区二区三区| 亚洲av片天天在线观看| 中国美女看黄片| 日本五十路高清| 观看免费一级毛片| 亚洲精品中文字幕一二三四区| 人妻夜夜爽99麻豆av| 老司机福利观看| 欧美一区二区国产精品久久精品 | 又爽又黄无遮挡网站| 亚洲精品美女久久久久99蜜臀| 色老头精品视频在线观看| 女警被强在线播放| 他把我摸到了高潮在线观看| 欧美日韩国产亚洲二区| 好男人在线观看高清免费视频| 美女午夜性视频免费| 99国产精品99久久久久| 99久久无色码亚洲精品果冻| 亚洲全国av大片| 中文字幕最新亚洲高清| 欧美性猛交╳xxx乱大交人| 999久久久精品免费观看国产| 成熟少妇高潮喷水视频| 国产午夜精品久久久久久| 久久国产精品影院| 两人在一起打扑克的视频| 日本精品一区二区三区蜜桃| 激情在线观看视频在线高清| 老司机午夜十八禁免费视频| 免费在线观看完整版高清| 亚洲国产精品合色在线| 美女 人体艺术 gogo| 美女黄网站色视频| 日韩国内少妇激情av| 国产欧美日韩一区二区精品| 久久精品国产亚洲av香蕉五月| 欧美日韩国产亚洲二区| 国产熟女午夜一区二区三区| 久久婷婷人人爽人人干人人爱| 91在线观看av| 亚洲色图 男人天堂 中文字幕| 亚洲国产看品久久| 黑人操中国人逼视频| 香蕉国产在线看| 国产精品av视频在线免费观看| 国产午夜精品论理片| 久久久国产成人免费| 免费av毛片视频| 在线十欧美十亚洲十日本专区| 黄色片一级片一级黄色片| 91麻豆精品激情在线观看国产| 成人永久免费在线观看视频| 久99久视频精品免费| 香蕉av资源在线| 少妇熟女aⅴ在线视频| 亚洲人成伊人成综合网2020| 日韩欧美三级三区| 激情在线观看视频在线高清|