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

    極端波浪載荷下導(dǎo)管架平臺(tái)抗倒塌性能評估

    2014-07-02 01:40:05朱本瑞陳國明
    關(guān)鍵詞:結(jié)構(gòu)方法

    朱本瑞,陳國明,林 紅,黃 超,劉 康

    (中國石油大學(xué)海洋油氣裝備與安全技術(shù)研究中心,山東青島266580)

    極端波浪載荷下導(dǎo)管架平臺(tái)抗倒塌性能評估

    朱本瑞,陳國明,林 紅,黃 超,劉 康

    (中國石油大學(xué)海洋油氣裝備與安全技術(shù)研究中心,山東青島266580)

    提出載荷序列增量分析(LSIA)方法,對極端波浪載荷下導(dǎo)管架平臺(tái)結(jié)構(gòu)極限狀態(tài)進(jìn)行分析。LSIA能夠考慮不同波高及可能的甲板上浪載荷對平臺(tái)側(cè)向載荷分布的影響,“追蹤”平臺(tái)在不同載荷重現(xiàn)期作用下的真實(shí)響應(yīng),從而克服當(dāng)前Pushover的缺點(diǎn)?;贚SIA等效承載能力曲線,提出倒塌重現(xiàn)期NCP,使得平臺(tái)抗倒塌能力評估指標(biāo)具體化。以K型和X型兩種不同結(jié)構(gòu)形式的導(dǎo)管架平臺(tái)為例,采用LSIA方法對其進(jìn)行抗倒塌性能評估,并與Pushover進(jìn)行分析對比。結(jié)果表明,LSIA具有更高的評估精度,能夠確定平臺(tái)結(jié)構(gòu)真實(shí)的倒塌狀態(tài)與失效模式。

    極端波浪;導(dǎo)管架平臺(tái);LSIA;評估指標(biāo);極限承載能力;失效模式

    臺(tái)風(fēng)引起的極端波浪載荷是海洋平臺(tái)倒塌破壞的主要原因,2002—2005年臺(tái)風(fēng)Lili、Ivan、Katrina和Rita相繼摧毀的122座導(dǎo)管架平臺(tái)中,大部分是由極端波浪引起的過載造成的[1]。GEORGE[2]認(rèn)為北海平臺(tái)倒塌失效與上浪載荷直接相關(guān)。目前,導(dǎo)管架平臺(tái)設(shè)計(jì)規(guī)范通常不考慮上浪導(dǎo)致的極大作用力,而是通過規(guī)定最小甲板高程避免上浪現(xiàn)象的發(fā)生,這種采用“最小甲板氣隙”的作法,對于某些沉降嚴(yán)重油田(如北海的Ekofisk油田沉降達(dá)到8 m,并以每年7~14 cm的速度下沉[3]),顯然不能夠?yàn)槠淙珘勖芷诘陌踩厶峁┯行ПU?。靜力推覆分析(Pushover)是確定結(jié)構(gòu)極限狀態(tài)及倒塌行為的有效方法[4-5]。針對導(dǎo)管架平臺(tái)而言,該方法是對結(jié)構(gòu)施加100年一遇的環(huán)境載荷,并將其逐級進(jìn)行放大,直至結(jié)構(gòu)發(fā)生倒塌,然后采用儲(chǔ)備強(qiáng)度系數(shù)(RSR)來衡量平臺(tái)結(jié)構(gòu)的抗力水平。目前,API、ISO以及DNV規(guī)范普遍推薦使用該方法[6-8],然而該方法存在以下不足:(1)造成平臺(tái)倒塌的環(huán)境載荷較大,其概率往往遠(yuǎn)小于1/100,因此基于100年一遇環(huán)境載荷進(jìn)行放大并不能夠真實(shí)反映平臺(tái)倒塌時(shí)受到的側(cè)向載荷分布;(2)平臺(tái)結(jié)構(gòu)倒塌失效時(shí)的波冠高度有可能大于甲板設(shè)計(jì)高程,而Pushover沒有考慮上浪載荷因素;(3)Pushover得到的極限承載能力曲線無法反映不同波高載荷作用下平臺(tái)結(jié)構(gòu)的響應(yīng),因而對平臺(tái)結(jié)構(gòu)的安全評估提供的信息不足[9]。為此,筆者提出載荷序列增量分析(LSIA)方法,以克服當(dāng)前Pushover方法的缺點(diǎn)。LSIA能夠充分考慮平臺(tái)側(cè)向載荷分布的變化及甲板上浪載荷因素,且其確定的評估指標(biāo)更加直接化、具體化。

    1 LSIA方法分析原理

    1.1 分析步驟

    與Pushover方法不同,LSIA以不同重現(xiàn)期的環(huán)境載荷作為漸進(jìn)加載序列來進(jìn)行極限承載能力分析,其原理如圖1所示。記N={N1,N2,…,Nn},Ni為第i個(gè)載荷序列對應(yīng)的重現(xiàn)期,其載荷參數(shù)可表示為{Vw,H,T,D,Vuc}i,Vw、H、T、D、Vuc分別為重現(xiàn)期Ni時(shí)的風(fēng)速、波高、周期、水深和流速。由此,對平臺(tái)在每種重現(xiàn)期Ni對應(yīng)的載荷作用下分別進(jìn)行一次靜力非線性分析,從而得到平臺(tái)結(jié)構(gòu)基底剪力與位移(或傾覆力矩與轉(zhuǎn)角)的關(guān)系。顯然,當(dāng)載荷重現(xiàn)期N較小時(shí),結(jié)構(gòu)響應(yīng)為線彈性,如圖1(a)、 (b)所示;當(dāng)N足夠大時(shí),結(jié)構(gòu)進(jìn)入非線性響應(yīng)階段,如圖1(c)、(d)所示;這樣,逐漸增大環(huán)境載荷重現(xiàn)期N,最終將得到一種載荷使得平臺(tái)結(jié)構(gòu)處于臨界倒塌狀態(tài),記其載荷重現(xiàn)期為NCP,微量增加NCP將使得結(jié)構(gòu)位移瞬間增大,平臺(tái)發(fā)生倒塌破壞,記此時(shí)結(jié)構(gòu)響應(yīng)曲線為倒塌承載能力曲線(collapse bearing capacity curve,CBCC),如圖1(e)所示;進(jìn)一步將平臺(tái)倒塌前每個(gè)Ni計(jì)算得到的基底剪力與位移繪于同一坐標(biāo)系中,即得到LSIA方法的等效承載能力曲線(equivalent bearing capacity curve,EBCC),如圖1(f)所示。

    顯然,相對于Pushover,LSIA方法不僅能夠確定平臺(tái)結(jié)構(gòu)倒塌極限狀態(tài),而且可以綜合反映平臺(tái)在不同重現(xiàn)期載荷作用下的結(jié)構(gòu)響應(yīng)信息,從而為服役海洋平臺(tái)全尺度載荷強(qiáng)度下的安全評估提供數(shù)據(jù)支持。由于LSIA方法是基于特定場址載荷要素的概率分布來確定加載序列,所以能夠考慮風(fēng)、浪、流各種載荷因素的相關(guān)性,使得平臺(tái)受到的側(cè)向載荷分布與真實(shí)海洋環(huán)境相一致,從而避免了當(dāng)前Pushover方法采用“單一”強(qiáng)度載荷進(jìn)行“盲目”放大的弊端。

    圖1 LSIA原理示意圖Fig.1 Sketch map of load sequence incremental analysis principle

    1.2 評估指標(biāo)

    基于Pushover方法確定的極限承載能力曲線(UBCC)可以細(xì)分為“線彈性響應(yīng)”、“彈塑性響應(yīng)”和“倒塌響應(yīng)”3個(gè)階段,如圖2所示。由圖2可知, 3個(gè)階段分別由設(shè)計(jì)環(huán)境載荷Rd、構(gòu)件開始屈服載荷Ry以及極限載荷Ru三個(gè)評估點(diǎn)確定。由此,定義儲(chǔ)備強(qiáng)度系數(shù)RSR=Ru/Rd以評估結(jié)構(gòu)整體抗力水平,延性系數(shù)γ=su/sy以表征結(jié)構(gòu)的塑性變形能力,以及“倒塌響應(yīng)”階段曲線變化趨勢表征結(jié)構(gòu)的倒塌性態(tài),文獻(xiàn)[5]認(rèn)為平臺(tái)達(dá)到極限狀態(tài)后其承載能力迅速下降時(shí)為脆性倒塌,緩慢降低則為延性倒塌,如圖2中曲線①和②所示。

    圖2 極限承載能力曲線Fig.2 Ultimate bearing capacity curve

    LSIA繼承了當(dāng)前Pushover方法的優(yōu)點(diǎn),并賦予了承載能力曲線新的物理意義,即載荷重現(xiàn)期。參照FEMA350[10]抗震性能評估思想,采用LSIA方法時(shí),將上述3個(gè)評估點(diǎn)重新定義為正常使用(normal operation)、可以使用(immediate operation)和預(yù)防倒塌(collapse prevention),對應(yīng)的載荷重現(xiàn)期分別記為NNO、NIO和NCP。由此,當(dāng)環(huán)境載荷重現(xiàn)期N<NNO時(shí),平臺(tái)結(jié)構(gòu)能夠正常使用;當(dāng)NNO<N<NIO時(shí),平臺(tái)結(jié)構(gòu)運(yùn)行風(fēng)險(xiǎn)較低,可以使用;當(dāng)平臺(tái)遭遇的環(huán)境載荷重現(xiàn)期NIO<N<NCP時(shí),則應(yīng)發(fā)布預(yù)警信息,擬定應(yīng)急方案,以防止結(jié)構(gòu)發(fā)生倒塌帶來的經(jīng)濟(jì)損失或人員傷亡;而當(dāng)監(jiān)測的環(huán)境載荷重現(xiàn)期N>NCP時(shí),平臺(tái)結(jié)構(gòu)將不足以抵抗外界載荷而發(fā)生倒塌破壞,此時(shí)應(yīng)對平臺(tái)采取必要的風(fēng)險(xiǎn)控制措施,以降低其失效的后果,如減少平臺(tái)質(zhì)量、關(guān)閉井口、人員撤離等。由此,將承載能力曲線劃分為藍(lán)、黃、橙、紅4級“服役狀態(tài)”監(jiān)測窗口,可為臺(tái)風(fēng)多發(fā)海域的平臺(tái)結(jié)構(gòu)完整性管理、風(fēng)險(xiǎn)預(yù)警以及安全決策提供一定參考。

    此外,倒塌重現(xiàn)期NCP與載荷參數(shù)直接相關(guān),即可以確定對應(yīng)的倒塌波高HCP或倒塌風(fēng)速VWCP等,從而使得衡量平臺(tái)結(jié)構(gòu)極限承載能力的參數(shù)更加具體化、形象化和直觀化,更利于結(jié)構(gòu)工程師的理解與應(yīng)用。另一方面,基于該方法確定的NCP可用于平臺(tái)結(jié)構(gòu)可靠性分析中,近似估算平臺(tái)結(jié)構(gòu)倒塌失效的概率,即Pf=P{N>NCP}。

    2 環(huán)境載荷計(jì)算

    采用LSIA進(jìn)行極端波浪載荷下平臺(tái)抗倒塌性能評估時(shí),需要額外考慮兩方面的計(jì)算:一是計(jì)算不同重現(xiàn)期對應(yīng)的環(huán)境參數(shù),以確定加載序列環(huán)境載荷值;二是當(dāng)波高大于平臺(tái)甲板高程時(shí),應(yīng)考慮甲板上浪載荷引起的抨擊載荷,即須確定上浪載荷計(jì)算方法。

    2.1 環(huán)境參數(shù)極值

    極端波浪載荷往往是由臺(tái)風(fēng)或超強(qiáng)臺(tái)風(fēng)引起的一種發(fā)生頻率相對較高的極值載荷,可采用三參數(shù)Weibull極值分布對極端環(huán)境載荷要素進(jìn)行推算[11],于是有

    由此,計(jì)算得到環(huán)境載荷參數(shù)極值xp為

    式中,參數(shù)α、μ和ξ分別為尺度參數(shù)、位置參數(shù)和形狀參數(shù),可根據(jù)極大似然估計(jì)法利用MS EXCEL“規(guī)劃求解”功能計(jì)算得到[12]。

    2.2 甲板上浪載荷

    迄今為止,甲板上浪載荷的計(jì)算尚沒有被普遍認(rèn)可的方法。現(xiàn)有的方法大致可分為側(cè)面輪廓法和細(xì)節(jié)構(gòu)件法兩大類。側(cè)面輪廓法是根據(jù)有效甲板濕水面積及水質(zhì)點(diǎn)壓力,采用拖曳力公式或動(dòng)量公式進(jìn)行計(jì)算;細(xì)節(jié)構(gòu)件方法則分別計(jì)算作用于甲板每根桿件上的波浪載荷,通過屏蔽系數(shù)考慮結(jié)構(gòu)構(gòu)件之間的交互作用,然后進(jìn)行疊加得到總載荷。鑒于后者需要建立所有桿件與設(shè)備模型,對計(jì)算機(jī)及軟件要求高,故本文選用側(cè)面輪廓法,其主要模型包括:API模型、DNV拍擊力模型、SHELL模型和MSL模型等,各種模型的詳細(xì)討論可參考文獻(xiàn)[13-14],采用API模型[6]計(jì)算甲板上浪時(shí),有

    式中,Fwd為甲板上浪載荷,N;ρ為海水密度,kg/m3; Cd為拖曳力系數(shù);Vx為甲板淹濕頂部處水質(zhì)點(diǎn)水平速度,m/s;Vuc為與波浪運(yùn)動(dòng)同向的流速,m/s;awkf為波浪運(yùn)動(dòng)系數(shù),臺(tái)風(fēng)時(shí)取0.88;acbf為導(dǎo)管架的海流阻擋系數(shù);A為甲板淹濕輪廓在波浪方向的投影面積,m2;Hwd為甲板上浪高度,m;B為甲板寬度,m。

    計(jì)算甲板上浪關(guān)鍵在于確定甲板淹濕位置處的波面高度及水質(zhì)點(diǎn)水平速度,可根據(jù)Stokes 5th波理論建立的色散關(guān)系方程組采用牛頓法編程求解出波長L和系數(shù)λ,然后代入波面公式和速度公式計(jì)算得到[15]。

    3 算例分析

    3.1 有限元模型及環(huán)境載荷

    以兩種不同結(jié)構(gòu)形式的導(dǎo)管架平臺(tái)作為算例,應(yīng)用LSIA方法對其抗倒塌性能進(jìn)行評估,兩平臺(tái)分別記為“K平臺(tái)”和“X平臺(tái)”,如圖3所示。K平臺(tái)的導(dǎo)管架由7層構(gòu)成,4根樁管成雙斜對稱結(jié)構(gòu), ROWA和ROWB斜撐為K型布置;上部組塊為3層,寬25 m,長45 m,質(zhì)量為6.033 kt(包括設(shè)備質(zhì)量);底層甲板高程為15 m;設(shè)計(jì)水深107 m;X平臺(tái)如圖3(b)所示,ROWA和ROWB斜撐為倒K型與K型交叉布置組成的X型結(jié)構(gòu),其余參數(shù)與K平臺(tái)相同。

    圖3 海洋平臺(tái)有限元模型Fig.3 Finite element model of platforms

    考慮到分析涉及結(jié)構(gòu)非線性倒塌階段力學(xué)行為研究,建模時(shí)應(yīng)充分考慮模擬單元的力學(xué)特性及材料非線性。選用3D彈塑性單元PIPE288建立平臺(tái)導(dǎo)管架部分,材料選取雙線性隨動(dòng)強(qiáng)化模型。PIPE288能夠通過OCEAN模塊自動(dòng)模擬計(jì)算波浪載荷,從而大大降低了波浪載荷下導(dǎo)管架平臺(tái)結(jié)構(gòu)力學(xué)分析的計(jì)算工作量。

    將上部組塊結(jié)構(gòu)進(jìn)行簡化,并采用線彈性單元PIPE16進(jìn)行模擬;組塊質(zhì)量平均分布在頂部4個(gè)節(jié)點(diǎn)上,采用MASS21單元模擬;忽略樁-土非線性相互作用,導(dǎo)管架底部采用固支約束,即在泥面處約束4條樁腿的所有自由度。

    對中國科學(xué)院南海海洋研究所1948—2008年間的環(huán)境載荷資料進(jìn)行統(tǒng)計(jì)分析,由極大似然估計(jì)法得到不同環(huán)境要素的Weibull分布參數(shù)估計(jì)值,見表1。

    表1 環(huán)境載荷Weibull參數(shù)估計(jì)Table 1 Weibull parameters estimate of environmental loads

    由表1中估計(jì)參數(shù),帶入式(2),可計(jì)算得到該海域不同重現(xiàn)期對應(yīng)的風(fēng)速、波高、周期及流速的極值,見表2,其中Hmax=1.721Hs,Tmax=1.27Tz,其推導(dǎo)過程詳見文獻(xiàn)[16]。表2中的上浪載荷由式(3)計(jì)算得到,其中系數(shù)Cd取2.5,acbf取0.8,B為45 m。

    表2 環(huán)境載荷參數(shù)及上浪載荷部分?jǐn)?shù)據(jù)Table 2 Partial environmental load parameters and wave-in-deck loads

    由表2可知,平臺(tái)100年一遇設(shè)計(jì)環(huán)境參數(shù)為:波高22.4 m,周期14.1 s、表層流速2.355 m/s,最大風(fēng)速44.2 m/s,從而計(jì)算得平臺(tái)正面(0°)最大波流載荷為15.2518 MN(對應(yīng)最大相位角為32°),風(fēng)載為1.8083 MN。此外,應(yīng)注意到當(dāng)載荷重現(xiàn)期為2000 a時(shí),平臺(tái)開始出現(xiàn)甲板上浪載荷;當(dāng)重現(xiàn)期為5.8×105a時(shí),上浪高度為4.6 m,上浪載荷為34.716 MN,達(dá)到了設(shè)計(jì)環(huán)境載荷的2倍左右,可見上浪載荷作用力之大。

    3.2 極限承載能力與性能評估

    采用當(dāng)前Pushover和LSIA方法對K平臺(tái)與X平臺(tái)進(jìn)行推覆分析,分別得到兩座平臺(tái)的極限承載能力曲線UBCC和倒塌承載能力曲線CBCC如圖4所示。

    圖4 當(dāng)前Pushover與LSIA方法分析結(jié)果對比Fig.4 Result comparison of current Pushover and LSIA

    由圖4可知:采用當(dāng)前Pushover方法時(shí),K平臺(tái)與X平臺(tái)均具有較高的強(qiáng)度儲(chǔ)備,RSR分別為4.87和5.50;但其延性性能差別較大,延性系數(shù)γ分別為1.60和3.88;從倒塌階段的響應(yīng)分析可知,平臺(tái)X極限承載能力下降緩慢,倒塌模式表現(xiàn)為延性倒塌,而K平臺(tái)則表現(xiàn)脆性倒塌。由此可知,X平臺(tái)在極限承載能力、延性能力以及倒塌模式方面均表現(xiàn)為優(yōu)良的力學(xué)特性,說明平臺(tái)空間拓?fù)浣Y(jié)構(gòu)合理,具有較好的載荷傳遞路徑及冗余性能。相對而言, K平臺(tái)在抵抗外界環(huán)境載荷,特別是載荷超過其極限承載能力時(shí)表現(xiàn)為脆弱性,不利于結(jié)構(gòu)在強(qiáng)臺(tái)風(fēng)或超強(qiáng)臺(tái)風(fēng)下“生存”。

    對比Pushover極限承載能力曲線UBCC與LSIA倒塌承載能力曲線CBCC可知,用兩種方法得到的K平臺(tái)承載能力曲線變化規(guī)律不同,顯然, Pushover過高地評估了K平臺(tái)的極限承載能力;而X平臺(tái)兩曲線變化趨勢一致,平臺(tái)在達(dá)到極限承載能力后,承載能力均緩慢下降,表現(xiàn)為良好的延性性能。

    根據(jù)LISA原理,作出兩座平臺(tái)的等效承載能力曲線EBCC如圖5所示。由圖5可知,平臺(tái)結(jié)構(gòu)在彈性響應(yīng)階段,當(dāng)發(fā)生上浪載荷時(shí),等效承載能力曲線出現(xiàn)拐點(diǎn),說明上浪載荷對平臺(tái)載荷分布及結(jié)構(gòu)響應(yīng)影響較大,進(jìn)一步確定兩座平臺(tái)的性能評估參數(shù)見表3。

    由表3可知,K平臺(tái)與X平臺(tái)倒塌重現(xiàn)期NCP分別為3.9×105a和5.8×105a,此時(shí)對應(yīng)的環(huán)境載荷參數(shù):倒塌波高HCP分別為32.9和33.3 m,風(fēng)速VWCP分別為64.5和65.2 m/s,因此若平臺(tái)遭遇的超強(qiáng)臺(tái)風(fēng)強(qiáng)度大于該倒塌環(huán)境載荷時(shí),平臺(tái)極有可能發(fā)生倒塌破壞。相對當(dāng)前Pushover分析結(jié)果,LSIA得到的K平臺(tái)和X平臺(tái)極限承載能力分別為75.163 7和78.802 4 MN,RSR分別降低9.5%和16.0%??梢?若不考慮平臺(tái)側(cè)向載荷分布模式以及甲板上浪等因素的影響,僅采用基于100年一遇的環(huán)境載荷進(jìn)行推覆分析,將使得評估結(jié)果偏于保守,容易造成較大的誤差。因此,對服役海洋平臺(tái)特別是抗力衰減嚴(yán)重的老齡平臺(tái)進(jìn)行極限承載能力評估時(shí),建議采用LSIA方法以獲得更加精確的結(jié)果。

    圖5 LSIA等效承載能力曲線Fig.5 Load sequence incremental analysis EBCC

    表3 基于LSIA的平臺(tái)性能評估參數(shù)Table 3 Performance evaluation parameters based on LSIA

    3.3 倒塌失效模式分析

    鑒于平臺(tái)承載能力曲線與結(jié)構(gòu)力學(xué)特性相關(guān),為進(jìn)一步闡明兩種方法得到的抗力曲線變化趨勢的迥異,揭示真實(shí)平臺(tái)結(jié)構(gòu)倒塌機(jī)制,提取平臺(tái)結(jié)構(gòu)的塑性應(yīng)變分布如圖6所示。

    由圖6(a)可知:兩種分析方法得到的K平臺(tái)失效模式發(fā)生顯著變化,從而導(dǎo)致其承載能力曲線變化規(guī)律不同。Pushover分析時(shí),K平臺(tái)主要表現(xiàn)為中層斜撐失效后,導(dǎo)致導(dǎo)管架ROW1和ROW2兩面所有樁腿過載而發(fā)生整體倒塌;LSIA則表現(xiàn)為上層斜撐失效使得平臺(tái)ROW2面樁腿發(fā)生彎曲失效,且平臺(tái)上部組塊與導(dǎo)管架連接構(gòu)件均發(fā)生塑性變形,顯然,這是由于甲板上浪載荷使得平臺(tái)側(cè)向載荷分布更加集中于頂部,從而引起塑性分布區(qū)上移,這進(jìn)一步說明考慮甲板上浪載荷的必要性。

    圖6 K平臺(tái)及X平臺(tái)倒塌失效模式對比Fig.6 Collapse failure modes comparison of platform K and platform X

    由圖6(b)可知,與K平臺(tái)結(jié)果相似,LSIA方法亦使得X平臺(tái)導(dǎo)管架上部塑性區(qū)擴(kuò)大,但兩種方法得到的平臺(tái)失效模式并沒有發(fā)生顯著的變化,均表現(xiàn)為中層斜撐與水平撐失效后,樁腿因無支撐長度增加而發(fā)生失穩(wěn),最終使得平臺(tái)形成機(jī)構(gòu)而發(fā)生倒塌。顯然,不同載荷分布模式并沒有改變結(jié)構(gòu)最終的倒塌模式,這說明X斜撐具有良好的結(jié)構(gòu)連續(xù)性與載荷傳遞功能,從而使得平臺(tái)上部載荷有效的傳遞給導(dǎo)管架結(jié)構(gòu),在抵抗外界載荷時(shí)表現(xiàn)為良好的抗力水平和延性性能。

    綜合K平臺(tái)與X平臺(tái)分析結(jié)果可知,當(dāng)前Pushover方法無法考慮載荷分布變化的影響,在進(jìn)行平臺(tái)結(jié)構(gòu)性能評估時(shí)結(jié)果偏于保守,特別針對發(fā)生上浪載荷的平臺(tái),甚至可能得到錯(cuò)誤的結(jié)論。此外,上浪載荷使得導(dǎo)管架與上部組塊連接構(gòu)件發(fā)生塑性變形,這些連接構(gòu)件處于飛濺區(qū),腐蝕嚴(yán)重,易成為平臺(tái)結(jié)構(gòu)的弱點(diǎn),因此若連接構(gòu)件抗剪強(qiáng)度不足或無法將載荷有效傳遞給導(dǎo)管架結(jié)構(gòu),一旦上浪載荷較大,則有可能導(dǎo)致平臺(tái)上部組塊側(cè)翻,從而發(fā)生臺(tái)風(fēng)過后觀察到的倒塌失效模式。

    4 結(jié) 論

    (1)LSIA方法能夠更好地“追蹤”平臺(tái)結(jié)構(gòu)真實(shí)的響應(yīng)以及對應(yīng)的載荷信息,在進(jìn)行平臺(tái)結(jié)構(gòu)抗倒塌性能評估時(shí),能夠有效地考慮甲板上浪載荷對平臺(tái)側(cè)向載荷分布的影響,具有更高的評估精度;基于該方法確定的評估指標(biāo)更加直觀、具體,是當(dāng)前Pushover方法的補(bǔ)充。

    (2)極端波浪載荷引起的上浪載荷對平臺(tái)側(cè)向載荷分布以及倒塌失效模式影響較大,巨大的甲板上浪載荷可能引發(fā)平臺(tái)上部組塊的側(cè)翻,對服役平臺(tái)安全評估時(shí)(尤其發(fā)生嚴(yán)重沉降的平臺(tái)),有必要考慮可能的甲板上浪載荷。

    (3)不同結(jié)構(gòu)形式的導(dǎo)管架平臺(tái)抗倒塌能力不同,X型構(gòu)造增加了平臺(tái)結(jié)構(gòu)的連續(xù)性,在抵抗外界水平載荷時(shí),具有良好的載荷傳遞路徑與力學(xué)特性,從而能夠有效提高平臺(tái)結(jié)構(gòu)的抗力與延性,降低結(jié)構(gòu)倒塌的風(fēng)險(xiǎn)。

    [1] PUSKAR F J,SPONG R E,KU A,et al.Assessment of fixed offshore platform performance in hurricane Ivan [R].OTC 18325,2006.

    [2] GEORGE Z.Wave crest heights and deck damage in hurricanes Ivan,Katrina and Rita[R].OTC 18620,2007.

    [3] MARLEY M,ETTERDAL B,GRIGORIAN H.Structural reliability assessment of Ekofisk jackets under extreme loading[R].OTC 13190,2001.

    [4] MOAN T.Development of accidental collapse limit state criteria for offshore structures[J].Structural Safety, 2009,31:124-135.

    [5] 朱本瑞,陳國明,劉康,等.深水導(dǎo)管架平臺(tái)抗臺(tái)風(fēng)極限承載能力[J].石油學(xué)報(bào),2014,35(1):152-158. ZHU Ben-rui,CHEN Guo-ming,LIU Kang,et al.Typhoon-resistant ultimate bearing capacity of deep-water jacket platforms[J].Acta Petrolei Sinica,2014,35(1): 152-158.

    [6] API RP 2A-WSD.Recommended practice for planning, design and constructing fixed offshore platforms—working stress design[S].American Petroleum Institute, 2007.

    [7] ISO 19902.Petroleum and natural gas industries-fixedsteel offshore structures[S].International Organization for Standardization,2007.

    [8] DNV.ULTIGUIDE-best practice guideline for use of nonlinear analysis methods in documentation of ultimate limit states for jacket type offshore structures[S].Norway:Det Norske Veritas,1999.

    [9] GOLAFSHANI A A,BAGHERI V,EBRAHIMIAN H,et al.Incremental wave analysis and its application to performance-based assessment of jacket platforms[J].Journal of Constructional Steel Research,2011,67(3):1649-1657.

    [10] FEMA-350.Recommended seismic design criteria for new steel moment-frame buildings[S].Washington D.C.:SAC Joint Venture,Federal Emergency Management Agency, 2000.

    [11] 陳團(tuán)海,陳國明,許亮斌.基于臺(tái)風(fēng)驗(yàn)證載荷的平臺(tái)時(shí)變可靠性分析與更新[J].中國石油大學(xué)學(xué)報(bào):自然科學(xué)版,2011,35(3):129-134. CHEN Tuan-hai,CHEN Guo-ming,XU Liang-bin. Time-dependent reliability analysis and updating of offshore platforms based on typhoon proof load[J].Journal of China University of Petroleum(Edition of Natural Science),2011,35(3):129-134.

    [12] 史景釗,任學(xué)軍,陳新昌,等.一種三參數(shù)Weibull分布極大似然估計(jì)的求解方法[J].河南科學(xué),2009, 27(7):832-834. SHI Jing-zhao,REN Xue-jun,CHEN Xin-chang,et al. A solution method for 3-parameter Weibull distribution of maximum likelihood estimation[J].Henan Science, 2009,27(7):832-834.

    [13] RAAIJ K,GUDMESTAD O T.Wave-in-deck loading on fixed steel jackets decks[J].Marine Structures,2007, 20(3):164-184.

    [14] RAAIJ K.Dynamic behaviour of jackets exposed to wavein-deck forces[D].Norway:University of Stavanger, 2005.

    [15] 董勝,石湘.海洋工程數(shù)值計(jì)算方法[M].青島:中國海洋大學(xué)出版社,2007:237-239.

    [16] 文圣常,余宙文.海洋理論與計(jì)算原理[M].北京:科學(xué)出版社,1984.

    (編輯 沈玉英)

    Anti-collapse performance assessment of jacket offshore platforms in extreme storm waves

    ZHU Ben-rui,CHEN Guo-ming,LIN Hong,HUANG Chao,LIU Kang
    (Centre for Offshore Engineering and Safety Technology in China University of Petroleum,Qingdao 266580,China)

    A novel approach called load sequence incremental analysis(LSIA)was established to estimate the ultimate limit states and accurate behaviors of jacket platforms against extreme environmental loading.This approach could take into account the effects of variation in wave height and wave-in-deck loading,and track the actual behavior of platform under different load return periods.Thereby,the shortages of current Pushover analysis(CPA)were overcome.Based on the equivalent bearing capacity curve obtained by LSIA,a parameter called collapse return period(NCP)was introduced,which made the evaluation indicator of platform performance more specifical.On this basis,LSIA was used for the anti-collapse performance assessment of K platform and X platform.The results show that LSIA can achieve better assessment precision than CPA,and be able to determine the real collapse state and failure mode of platforms.

    extreme wave;jacket platform;LSIA;evaluation index;ultimate bearing capacity;failure mode

    TE 951

    :A

    1673-5005(2014)03-0135-07

    10.3969/j.issn.1673-5005.2014.03.022

    2013-12-21

    國家自然科學(xué)基金項(xiàng)目(51079159,51209218);中央高?;究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)(13CX06079A)

    朱本瑞(1986-),男,博士研究生,主要研究方向?yàn)楹Q笫脱b備強(qiáng)度與可靠性等。E-mail:zhubenrui@163.com。

    猜你喜歡
    結(jié)構(gòu)方法
    《形而上學(xué)》△卷的結(jié)構(gòu)和位置
    論結(jié)構(gòu)
    中華詩詞(2019年7期)2019-11-25 01:43:04
    新型平衡塊結(jié)構(gòu)的應(yīng)用
    模具制造(2019年3期)2019-06-06 02:10:54
    學(xué)習(xí)方法
    可能是方法不對
    論《日出》的結(jié)構(gòu)
    用對方法才能瘦
    Coco薇(2016年2期)2016-03-22 02:42:52
    四大方法 教你不再“坐以待病”!
    Coco薇(2015年1期)2015-08-13 02:47:34
    賺錢方法
    捕魚
    最近视频中文字幕2019在线8| 一级毛片精品| 亚洲中文字幕日韩| 国产欧美日韩一区二区三| 视频区欧美日本亚洲| 天堂√8在线中文| 亚洲欧美精品综合久久99| 午夜精品在线福利| 两性夫妻黄色片| 亚洲精品久久国产高清桃花| 亚洲精品久久国产高清桃花| 免费看十八禁软件| 久久久久国内视频| 1024手机看黄色片| 亚洲国产色片| 狂野欧美白嫩少妇大欣赏| 桃色一区二区三区在线观看| 免费在线观看成人毛片| 欧美3d第一页| 亚洲欧洲精品一区二区精品久久久| 在线十欧美十亚洲十日本专区| 亚洲精品在线观看二区| 免费看美女性在线毛片视频| 97超级碰碰碰精品色视频在线观看| 在线观看日韩欧美| www日本黄色视频网| 久久久成人免费电影| 成年女人永久免费观看视频| 精品久久久久久久末码| 久久草成人影院| 看黄色毛片网站| 我的老师免费观看完整版| 黄色女人牲交| 搡老妇女老女人老熟妇| 亚洲五月婷婷丁香| 午夜精品在线福利| 午夜福利成人在线免费观看| 国产伦在线观看视频一区| 久久精品亚洲精品国产色婷小说| 99热精品在线国产| svipshipincom国产片| 欧美中文综合在线视频| 人妻丰满熟妇av一区二区三区| 99国产精品一区二区蜜桃av| 变态另类丝袜制服| 国产精品久久久久久人妻精品电影| 久久九九热精品免费| 色av中文字幕| 国产伦一二天堂av在线观看| 99久久久亚洲精品蜜臀av| 午夜福利欧美成人| 中文字幕熟女人妻在线| 男人和女人高潮做爰伦理| 国产v大片淫在线免费观看| 制服人妻中文乱码| 欧美日本视频| 手机成人av网站| 亚洲精品在线美女| 日本五十路高清| 免费电影在线观看免费观看| 在线观看美女被高潮喷水网站 | 制服丝袜大香蕉在线| 香蕉国产在线看| 老司机午夜福利在线观看视频| 久久人妻av系列| 欧美日本视频| 精品国产超薄肉色丝袜足j| 国产精品久久久久久亚洲av鲁大| 久久99热这里只有精品18| 成年女人永久免费观看视频| 中亚洲国语对白在线视频| 成人无遮挡网站| 国产视频内射| 午夜成年电影在线免费观看| 脱女人内裤的视频| 色老头精品视频在线观看| 午夜视频精品福利| 日韩 欧美 亚洲 中文字幕| 亚洲熟妇中文字幕五十中出| 亚洲九九香蕉| 久久久久国内视频| 欧美黄色淫秽网站| 在线播放国产精品三级| 久久精品国产综合久久久| 国产精品99久久99久久久不卡| 日本黄大片高清| 欧美av亚洲av综合av国产av| 免费高清视频大片| 日韩av在线大香蕉| 精品一区二区三区四区五区乱码| 亚洲第一电影网av| 亚洲午夜理论影院| 亚洲熟女毛片儿| 1000部很黄的大片| 久久亚洲精品不卡| 久久精品综合一区二区三区| 久久久久九九精品影院| 狂野欧美激情性xxxx| 亚洲第一电影网av| 国产毛片a区久久久久| 日韩欧美 国产精品| 男女下面进入的视频免费午夜| 欧美黑人欧美精品刺激| 最近最新中文字幕大全免费视频| 精品一区二区三区视频在线 | 女警被强在线播放| 中出人妻视频一区二区| 亚洲欧美日韩高清在线视频| 两性午夜刺激爽爽歪歪视频在线观看| 国产精品国产高清国产av| 99久久久亚洲精品蜜臀av| 国产精品久久久久久亚洲av鲁大| 日本五十路高清| 国产精品av视频在线免费观看| 变态另类丝袜制服| 69av精品久久久久久| 18禁裸乳无遮挡免费网站照片| 99久久国产精品久久久| 法律面前人人平等表现在哪些方面| 啦啦啦观看免费观看视频高清| 日本成人三级电影网站| 此物有八面人人有两片| 精品国产美女av久久久久小说| 天天添夜夜摸| 亚洲精华国产精华精| 国产成人一区二区三区免费视频网站| 国产69精品久久久久777片 | 宅男免费午夜| 免费大片18禁| 无人区码免费观看不卡| 可以在线观看毛片的网站| 午夜两性在线视频| 一a级毛片在线观看| 国产精品av久久久久免费| 国产欧美日韩一区二区精品| 欧美极品一区二区三区四区| 国产成人av激情在线播放| 天天一区二区日本电影三级| 999久久久精品免费观看国产| 男女午夜视频在线观看| 国产精品一区二区三区四区久久| 亚洲精品美女久久av网站| 夜夜爽天天搞| 国产69精品久久久久777片 | avwww免费| 精品乱码久久久久久99久播| 欧美日韩中文字幕国产精品一区二区三区| 一本久久中文字幕| 欧美另类亚洲清纯唯美| 国产熟女xx| 我的老师免费观看完整版| av欧美777| 亚洲性夜色夜夜综合| 1024手机看黄色片| 日本撒尿小便嘘嘘汇集6| av片东京热男人的天堂| 欧美成人免费av一区二区三区| 国产av麻豆久久久久久久| 999久久久精品免费观看国产| 午夜免费观看网址| 精品久久久久久久毛片微露脸| 成年人黄色毛片网站| 特大巨黑吊av在线直播| 超碰成人久久| 神马国产精品三级电影在线观看| 久久精品国产亚洲av香蕉五月| 他把我摸到了高潮在线观看| 免费人成视频x8x8入口观看| 国产v大片淫在线免费观看| 亚洲av片天天在线观看| 黄色 视频免费看| 亚洲精品美女久久久久99蜜臀| 一夜夜www| 性欧美人与动物交配| 波多野结衣巨乳人妻| 18禁观看日本| 亚洲国产欧美人成| 日韩高清综合在线| 天堂影院成人在线观看| 欧美日韩精品网址| 巨乳人妻的诱惑在线观看| 欧美性猛交╳xxx乱大交人| 两人在一起打扑克的视频| 中文字幕av在线有码专区| 日本五十路高清| 久久久精品欧美日韩精品| 国产成人一区二区三区免费视频网站| 亚洲天堂国产精品一区在线| 无限看片的www在线观看| 免费在线观看影片大全网站| 亚洲激情在线av| 久久天躁狠狠躁夜夜2o2o| 亚洲国产中文字幕在线视频| 国产乱人伦免费视频| 亚洲激情在线av| 欧美zozozo另类| 国产激情久久老熟女| 亚洲欧美日韩东京热| 国产野战对白在线观看| 亚洲精品美女久久av网站| 非洲黑人性xxxx精品又粗又长| 小蜜桃在线观看免费完整版高清| 亚洲 国产 在线| 国产成人福利小说| 国产精品 国内视频| 婷婷丁香在线五月| 久久精品国产综合久久久| 日本五十路高清| 噜噜噜噜噜久久久久久91| 成人午夜高清在线视频| 亚洲国产色片| 高清在线国产一区| 免费无遮挡裸体视频| 国产精品久久久av美女十八| 嫩草影院精品99| 午夜激情福利司机影院| 99国产精品99久久久久| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲真实伦在线观看| 久久精品国产亚洲av香蕉五月| 一个人免费在线观看的高清视频| 免费观看人在逋| 十八禁人妻一区二区| 99热6这里只有精品| 看片在线看免费视频| 午夜视频精品福利| 视频区欧美日本亚洲| 一进一出抽搐动态| 可以在线观看毛片的网站| 日本黄色视频三级网站网址| 久久精品国产清高在天天线| av视频在线观看入口| 成人av在线播放网站| 亚洲avbb在线观看| 99国产精品一区二区三区| www.自偷自拍.com| 淫秽高清视频在线观看| aaaaa片日本免费| 中文字幕av在线有码专区| 男人舔奶头视频| 久久香蕉国产精品| 在线看三级毛片| 啦啦啦韩国在线观看视频| 91字幕亚洲| 观看美女的网站| 真人一进一出gif抽搐免费| 久久久久精品国产欧美久久久| 婷婷亚洲欧美| 亚洲欧美日韩无卡精品| 我要搜黄色片| 国产单亲对白刺激| 激情在线观看视频在线高清| 日本 av在线| 夜夜爽天天搞| 日本三级黄在线观看| 一个人看视频在线观看www免费 | 欧美一区二区精品小视频在线| 亚洲熟妇熟女久久| 免费在线观看成人毛片| 免费看光身美女| 亚洲精品粉嫩美女一区| 中文在线观看免费www的网站| 亚洲精品一区av在线观看| 好男人在线观看高清免费视频| 757午夜福利合集在线观看| 国产亚洲欧美在线一区二区| 中文字幕久久专区| 亚洲av免费在线观看| 欧美一级毛片孕妇| 亚洲专区中文字幕在线| 国产真实乱freesex| 国产精品99久久99久久久不卡| av在线天堂中文字幕| 夜夜看夜夜爽夜夜摸| 黄色视频,在线免费观看| 精品国产乱子伦一区二区三区| 国产又色又爽无遮挡免费看| 国产v大片淫在线免费观看| 成人国产一区最新在线观看| 在线观看免费午夜福利视频| 亚洲成人久久性| 中亚洲国语对白在线视频| 国产伦一二天堂av在线观看| 久久久久久国产a免费观看| 九色国产91popny在线| avwww免费| 国产97色在线日韩免费| 免费看a级黄色片| 欧美色视频一区免费| 国产精品1区2区在线观看.| 久久中文字幕一级| 91九色精品人成在线观看| 88av欧美| 日本一本二区三区精品| 日韩欧美三级三区| 久久久久精品国产欧美久久久| 精华霜和精华液先用哪个| 无人区码免费观看不卡| 好男人电影高清在线观看| 免费人成视频x8x8入口观看| 日韩精品中文字幕看吧| 国产高清有码在线观看视频| 欧美乱色亚洲激情| 少妇裸体淫交视频免费看高清| 国产激情欧美一区二区| 一二三四社区在线视频社区8| 日韩欧美国产在线观看| 久久天堂一区二区三区四区| 国产高清videossex| 五月伊人婷婷丁香| 久久草成人影院| 精品人妻1区二区| 精品一区二区三区视频在线 | 18禁国产床啪视频网站| 日韩欧美国产在线观看| 18禁国产床啪视频网站| 观看免费一级毛片| avwww免费| 99久久综合精品五月天人人| 天天躁狠狠躁夜夜躁狠狠躁| 久久久国产成人精品二区| 日本免费一区二区三区高清不卡| 欧美日韩一级在线毛片| 999久久久国产精品视频| 久久久色成人| av在线蜜桃| a在线观看视频网站| 曰老女人黄片| 日韩av在线大香蕉| 国产精品98久久久久久宅男小说| 精品久久久久久久末码| 国产亚洲av高清不卡| 91久久精品国产一区二区成人 | 日本黄色视频三级网站网址| 给我免费播放毛片高清在线观看| 在线a可以看的网站| 男女视频在线观看网站免费| 一二三四在线观看免费中文在| 91久久精品国产一区二区成人 | 亚洲中文字幕一区二区三区有码在线看 | 国产精品电影一区二区三区| 亚洲精品粉嫩美女一区| 一区二区三区国产精品乱码| or卡值多少钱| 亚洲成av人片在线播放无| 亚洲av日韩精品久久久久久密| 国产激情欧美一区二区| 亚洲av日韩精品久久久久久密| 国产激情欧美一区二区| 久久久国产精品麻豆| 成人一区二区视频在线观看| 欧美黑人欧美精品刺激| 脱女人内裤的视频| 欧美性猛交黑人性爽| 日日夜夜操网爽| 哪里可以看免费的av片| 国产精品野战在线观看| 亚洲狠狠婷婷综合久久图片| 韩国av一区二区三区四区| 最近视频中文字幕2019在线8| 小蜜桃在线观看免费完整版高清| 亚洲狠狠婷婷综合久久图片| 欧美乱色亚洲激情| 99国产精品一区二区蜜桃av| 国产精品一区二区精品视频观看| 久久天躁狠狠躁夜夜2o2o| 亚洲精品一卡2卡三卡4卡5卡| 色播亚洲综合网| 老司机在亚洲福利影院| 18禁裸乳无遮挡免费网站照片| 在线观看免费视频日本深夜| 亚洲成av人片免费观看| 两个人的视频大全免费| 国产成人一区二区三区免费视频网站| 人人妻人人澡欧美一区二区| 精品久久蜜臀av无| 国产精品av视频在线免费观看| 法律面前人人平等表现在哪些方面| 国产主播在线观看一区二区| 亚洲欧洲精品一区二区精品久久久| 免费av毛片视频| 国产精品永久免费网站| 色吧在线观看| 夜夜躁狠狠躁天天躁| 亚洲18禁久久av| 免费在线观看日本一区| 淫秽高清视频在线观看| 超碰成人久久| 色精品久久人妻99蜜桃| 国产一级毛片七仙女欲春2| 好男人在线观看高清免费视频| 我的老师免费观看完整版| 麻豆成人午夜福利视频| 亚洲 国产 在线| 亚洲人成伊人成综合网2020| 美女大奶头视频| 国产v大片淫在线免费观看| 色噜噜av男人的天堂激情| 国产精品 欧美亚洲| 国产欧美日韩一区二区三| 真实男女啪啪啪动态图| 国产欧美日韩精品亚洲av| 成人欧美大片| 久久久精品欧美日韩精品| 久久精品亚洲精品国产色婷小说| 五月玫瑰六月丁香| 欧美一级毛片孕妇| 国产精品99久久99久久久不卡| 国产91精品成人一区二区三区| 精品熟女少妇八av免费久了| 免费av不卡在线播放| 在线观看美女被高潮喷水网站 | 国产91精品成人一区二区三区| 成人鲁丝片一二三区免费| 久久久久久九九精品二区国产| 欧美成人一区二区免费高清观看 | 香蕉丝袜av| 中文字幕久久专区| 91麻豆精品激情在线观看国产| 成人亚洲精品av一区二区| 最好的美女福利视频网| 亚洲成人久久爱视频| 久久香蕉精品热| 国产成人精品无人区| 中文资源天堂在线| 久久精品综合一区二区三区| 高潮久久久久久久久久久不卡| 国产蜜桃级精品一区二区三区| 99久久精品国产亚洲精品| 亚洲精品一区av在线观看| 中出人妻视频一区二区| 1024香蕉在线观看| 黑人欧美特级aaaaaa片| 成人特级av手机在线观看| 一本综合久久免费| 免费搜索国产男女视频| 成人特级黄色片久久久久久久| 久9热在线精品视频| 欧美+亚洲+日韩+国产| 国产精品亚洲av一区麻豆| 美女午夜性视频免费| 国产高清videossex| 亚洲av成人av| 一本综合久久免费| 欧美xxxx黑人xx丫x性爽| 欧美日韩乱码在线| 亚洲国产精品合色在线| 午夜免费激情av| 免费在线观看影片大全网站| 日本黄色片子视频| 久久99热这里只有精品18| 夜夜躁狠狠躁天天躁| 真人做人爱边吃奶动态| 日本免费a在线| 亚洲午夜精品一区,二区,三区| 一个人观看的视频www高清免费观看 | 久久精品国产亚洲av香蕉五月| 国产精品一区二区免费欧美| or卡值多少钱| 老司机在亚洲福利影院| 好男人在线观看高清免费视频| 国产精品久久久久久亚洲av鲁大| 午夜福利成人在线免费观看| 深夜精品福利| 亚洲av成人一区二区三| 一个人观看的视频www高清免费观看 | 亚洲av美国av| 国产成人精品久久二区二区免费| 欧美在线一区亚洲| 国产极品精品免费视频能看的| 美女 人体艺术 gogo| 国产精品久久久av美女十八| 国产又黄又爽又无遮挡在线| 色综合站精品国产| 色综合婷婷激情| 亚洲av熟女| 亚洲精品在线美女| 日本黄色片子视频| 1024香蕉在线观看| 9191精品国产免费久久| 嫩草影院入口| 999久久久国产精品视频| 久久草成人影院| 18禁国产床啪视频网站| 亚洲国产欧美网| 免费搜索国产男女视频| 成人av一区二区三区在线看| 老司机午夜福利在线观看视频| 久久精品国产清高在天天线| 人人妻人人看人人澡| 国产高清有码在线观看视频| 不卡一级毛片| 欧美黄色淫秽网站| 丰满人妻一区二区三区视频av | 在线观看舔阴道视频| 69av精品久久久久久| 日本黄大片高清| 最新在线观看一区二区三区| 免费电影在线观看免费观看| netflix在线观看网站| 母亲3免费完整高清在线观看| 麻豆久久精品国产亚洲av| 成人鲁丝片一二三区免费| 国产精品一及| 91九色精品人成在线观看| 麻豆av在线久日| 又粗又爽又猛毛片免费看| 美女扒开内裤让男人捅视频| 国产成人福利小说| 国产精品国产高清国产av| 1024手机看黄色片| 五月伊人婷婷丁香| 国产探花在线观看一区二区| 国产伦人伦偷精品视频| 欧美日韩黄片免| 欧美中文日本在线观看视频| 色吧在线观看| 一本综合久久免费| 久久性视频一级片| 2021天堂中文幕一二区在线观| svipshipincom国产片| 久久精品综合一区二区三区| 网址你懂的国产日韩在线| 麻豆成人午夜福利视频| 国产又色又爽无遮挡免费看| 欧美最黄视频在线播放免费| 精品免费久久久久久久清纯| 999精品在线视频| 狂野欧美激情性xxxx| 日本精品一区二区三区蜜桃| 两性午夜刺激爽爽歪歪视频在线观看| 香蕉国产在线看| 91av网一区二区| 人妻丰满熟妇av一区二区三区| 日本a在线网址| 三级国产精品欧美在线观看 | 五月玫瑰六月丁香| 天天一区二区日本电影三级| 久久久国产成人免费| 亚洲国产看品久久| 国产av麻豆久久久久久久| 免费搜索国产男女视频| www.自偷自拍.com| 日韩欧美国产在线观看| 国产精品av久久久久免费| 日韩中文字幕欧美一区二区| 国产亚洲欧美98| 亚洲精品色激情综合| 一个人看视频在线观看www免费 | 天堂av国产一区二区熟女人妻| 日韩av在线大香蕉| 男人舔奶头视频| 午夜久久久久精精品| 嫩草影院精品99| 国产精品久久久久久久电影 | 日本熟妇午夜| 成人性生交大片免费视频hd| 久久精品aⅴ一区二区三区四区| 一区二区三区国产精品乱码| 这个男人来自地球电影免费观看| 亚洲乱码一区二区免费版| 亚洲av成人精品一区久久| 91在线精品国自产拍蜜月 | 国产一区二区在线av高清观看| 免费无遮挡裸体视频| 在线免费观看不下载黄p国产 | 脱女人内裤的视频| 啪啪无遮挡十八禁网站| 波多野结衣高清无吗| 最近最新中文字幕大全电影3| 日本免费a在线| 国产高潮美女av| 日本三级黄在线观看| 国产精品自产拍在线观看55亚洲| 国产精品久久久久久久电影 | 天天一区二区日本电影三级| 精品熟女少妇八av免费久了| 日韩欧美一区二区三区在线观看| 91九色精品人成在线观看| 亚洲色图 男人天堂 中文字幕| 美女午夜性视频免费| 精品久久久久久久人妻蜜臀av| 日本 欧美在线| 国产野战对白在线观看| 精品久久久久久久久久免费视频| 亚洲av日韩精品久久久久久密| h日本视频在线播放| 熟女人妻精品中文字幕| 精华霜和精华液先用哪个| 天堂av国产一区二区熟女人妻| 99国产精品一区二区蜜桃av| 999久久久国产精品视频| 欧美另类亚洲清纯唯美| 欧美一区二区国产精品久久精品| 亚洲人成网站在线播放欧美日韩| 老司机午夜福利在线观看视频| 亚洲国产欧美人成| 中文亚洲av片在线观看爽| 亚洲欧美精品综合久久99| 久久久久久大精品| 亚洲熟女毛片儿| 色综合欧美亚洲国产小说| 国产欧美日韩精品亚洲av| 一二三四在线观看免费中文在| 欧美黄色淫秽网站| 国产一区二区在线av高清观看| 在线观看免费视频日本深夜| 99久久综合精品五月天人人| 午夜福利免费观看在线| 国产精品久久久久久久电影 | 欧美不卡视频在线免费观看| 欧美日韩福利视频一区二区| 91麻豆精品激情在线观看国产| 757午夜福利合集在线观看| 99在线视频只有这里精品首页|