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

    超深水水下分離器承壓性能數(shù)值分析與試驗(yàn)研究

    2017-01-17 08:36:08李秀美陳國(guó)明張慎顏葛玖浩劉秀全
    關(guān)鍵詞:樣機(jī)筒體分離器

    李秀美, 陳國(guó)明, 李 偉, 張慎顏, 葛玖浩, 劉秀全

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

    超深水水下分離器承壓性能數(shù)值分析與試驗(yàn)研究

    李秀美, 陳國(guó)明, 李 偉, 張慎顏, 葛玖浩, 劉秀全

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

    以中國(guó)自主設(shè)計(jì)制造的2 000 m超深水水下分離器為研究對(duì)象,基于圖算法和一致缺陷模態(tài)法分析分離器樣機(jī)極限承載能力;建立包含嵌入式接管子模型的水下分離器精細(xì)數(shù)值模型,研究馬鞍形焊縫應(yīng)力分布規(guī)律;設(shè)計(jì)高壓艙試驗(yàn)驗(yàn)證水下分離器樣機(jī)的物理承壓能力,并對(duì)比分析高壓艙試驗(yàn)數(shù)據(jù)和數(shù)值計(jì)算結(jié)果。研究表明:水下分離器樣機(jī)失穩(wěn)模式為塑性失穩(wěn),下放工況應(yīng)該主要評(píng)估其結(jié)構(gòu)強(qiáng)度;接管馬鞍形焊縫的鞍點(diǎn)和冠點(diǎn)以及N1接管焊縫的角度φ為45°和135°位置為應(yīng)力特征點(diǎn);試驗(yàn)樣機(jī)在高壓艙試驗(yàn)過(guò)程中整體應(yīng)力水平處于線彈性變形階段,能夠滿足2 000 m水深的承壓強(qiáng)度要求,模擬結(jié)果和試驗(yàn)結(jié)果吻合較好。

    超深水水下分離器; 嵌入式接管; 子模型; 馬鞍形焊縫; 高壓艙試驗(yàn)

    深水水下分離器可在海底完成氣液分離和油水分離,有效解決深海油氣生產(chǎn)帶來(lái)的流動(dòng)安全保障與經(jīng)濟(jì)效益的平衡問(wèn)題。深水水下分離器在下放到目標(biāo)水深海域時(shí),承受超高外壓載荷,須評(píng)估承壓結(jié)構(gòu)強(qiáng)度和穩(wěn)定性。由于深水水下分離器這種厚壁短圓筒結(jié)構(gòu)發(fā)生塑性失穩(wěn)所需的外壓極大,承重結(jié)構(gòu)強(qiáng)度是否滿足要求成為評(píng)估重點(diǎn)。目前,深水水下承壓殼體結(jié)構(gòu)強(qiáng)度的研究大都不考慮殼體接管的影響,將模型簡(jiǎn)化成筒體結(jié)構(gòu)進(jìn)行分析,最常見(jiàn)的簡(jiǎn)化形式為海底承壓管道。Urthaler等[1]分析了深水水下承壓結(jié)構(gòu)的薄膜應(yīng)力、軸向應(yīng)力和徑向應(yīng)力分布特點(diǎn)。An等[2]建立了包含一種新型夾層材料的深水海底夾層管道的有限元模型,分析該承壓結(jié)構(gòu)在外壓和縱向彎曲載荷下的極限強(qiáng)度。Li[3]、Netto[4]、Lee[5]等采用數(shù)值計(jì)算和試驗(yàn)方法分析了含腐蝕缺陷的海底管道承壓結(jié)構(gòu)的剩余強(qiáng)度。Zheng等[6]研究了帶有凹陷的單壁管和管中管結(jié)構(gòu)在外壓載荷下強(qiáng)度問(wèn)題;深水水下分離器在海底完成油氣水三相分離,接管是不可忽略的功能性結(jié)構(gòu),為典型的嵌入式接管。將水下分離器承壓殼體簡(jiǎn)化成筒體難以獲得典型多接管承壓結(jié)構(gòu)的應(yīng)力特征。嵌入式接管環(huán)焊縫為復(fù)雜的、大直徑開(kāi)孔的馬鞍形焊縫,這種焊縫造成管內(nèi)外變形增大從而加劇了焊縫處的應(yīng)力集中,是整個(gè)超深水水下分離器殼體結(jié)構(gòu)的薄弱環(huán)節(jié)。Choo等[7]采用極限分析方法研究了外壓載荷下復(fù)合夾層結(jié)構(gòu)接管和容器相交區(qū)域的強(qiáng)度。葛玖浩等[8]建立多接管水下分離器數(shù)值模型,分析了接管參數(shù)對(duì)整體承壓結(jié)構(gòu)應(yīng)力分布的影響。李志剛等[9]采用有限元分析方法校核了1 500 m水下分離器的總體塑性變形,但未對(duì)測(cè)點(diǎn)的選取進(jìn)行科學(xué)論證。上述文獻(xiàn)從不同角度開(kāi)展多接管承壓結(jié)構(gòu)整體強(qiáng)度分析研究,但考慮包含馬鞍形焊縫的嵌入式接管結(jié)構(gòu),并從“數(shù)?!焙汀拔锬!苯嵌认到y(tǒng)研究2 000 m超深水水下分離器結(jié)構(gòu)應(yīng)力響應(yīng)方面的文獻(xiàn)未見(jiàn)報(bào)道。鑒于此,筆者以中國(guó)自主設(shè)計(jì)的2 000 m超深水水下分離器為研究對(duì)象,在分析其極限承載能力的基礎(chǔ)上,確定下放工況應(yīng)該主要考慮承壓結(jié)構(gòu)的強(qiáng)度問(wèn)題,設(shè)計(jì)高壓艙試驗(yàn)驗(yàn)證分離器樣機(jī)的承壓性能。

    1 水下分離器極限承載能力分析

    1.1 圖算法

    GB 150《壓力容器》規(guī)定:外壓厚壁容器必須首先評(píng)估其穩(wěn)定性,然后再進(jìn)行承壓結(jié)構(gòu)強(qiáng)度的校核。作為GB 150外壓容器設(shè)計(jì)的推薦做法,圖算法可以避開(kāi)Mises公式求解臨界壓力的繁瑣計(jì)算,不受彈性失穩(wěn)和彈塑性失穩(wěn)條件的限制,通過(guò)查圖即可獲得外壓容器的失穩(wěn)臨界壓力。圖算法首先通過(guò)承壓結(jié)構(gòu)的幾何參數(shù)在外壓應(yīng)變系數(shù)A曲線中確定系數(shù)A,然后通過(guò)系數(shù)A在外壓應(yīng)力系數(shù)B曲線中確定系數(shù)B;由于目標(biāo)水下分離器樣機(jī)屬于典型的厚壁容器,最后按照下式計(jì)算臨界失穩(wěn)壓力pcr[10]:

    其中

    (1)

    目標(biāo)水下分離器筒體長(zhǎng)度為6 000 mm,內(nèi)徑為1 200 mm,有效厚度為85 mm,兩端采用半球形封頭封閉。根據(jù)GB 150圖4-1得到計(jì)算長(zhǎng)度L為筒體的總長(zhǎng)度加上每個(gè)凸形封頭曲面深度的1/3,其值為6 457 mm,由GB 150圖4-2查得系數(shù)A=4.5×10-2,由GB 150圖4-6查得系數(shù)B=216 MPa,最終根據(jù)式(1)得到目標(biāo)水下分離器的臨界失穩(wěn)壓力pcr=49.9 MPa。

    1.2 一致缺陷模態(tài)法

    圖算法是以不含接管的外壓筒體為分析對(duì)象,計(jì)算得到的臨界失穩(wěn)壓力未考慮目標(biāo)水下分離器接管的影響。一致缺陷模態(tài)分析方法以數(shù)值分析軟件建立的實(shí)體有限元模型為基礎(chǔ),在外壓容器穩(wěn)定性評(píng)估領(lǐng)域得到廣泛運(yùn)用。該方法認(rèn)為一階屈曲模態(tài)代表結(jié)構(gòu)最不利的幾何形狀偏差,體現(xiàn)了結(jié)構(gòu)失穩(wěn)時(shí)的形態(tài),分析時(shí)將特征值屈曲分析得到一階模態(tài)乘以一個(gè)系數(shù)模擬承壓結(jié)構(gòu)的初始幾何缺陷[11],并進(jìn)行非線性屈曲分析求得臨界失穩(wěn)壓力。

    參照文獻(xiàn)[8]建立水下分離器樣機(jī)整體有限元模型,進(jìn)行特征值屈曲分析得到的一階屈曲模態(tài)如圖1所示,一階特征屈曲壓力p1=160.98 MPa?;诓牧侠硐霃椝苄岳碚摵蛶缀未笞冃卫碚?進(jìn)行非線性屈曲分析得到筒體典型位置處的載荷-位移曲線如圖2所示。由圖1可知,分離器樣機(jī)失穩(wěn)時(shí),整體承壓結(jié)構(gòu)呈現(xiàn)接管兩側(cè)筒體向內(nèi)凹陷、失穩(wěn)波數(shù)n=2的變形。圖2中的載荷度α表示外壓載荷和一階特征屈曲壓力p1的比值,由圖2可知,當(dāng)載荷度為0.343時(shí),殼體結(jié)構(gòu)位移出現(xiàn)“無(wú)窮大”變化趨勢(shì),承壓結(jié)構(gòu)失穩(wěn),失穩(wěn)極限載荷pcr=55.2 MPa,大于圖算法分析得到的臨界載荷,GB 150圖算法分析結(jié)果相對(duì)保守。

    由經(jīng)典厚壁圓筒應(yīng)力分析理論可知,目標(biāo)水下分離器在臨界失穩(wěn)壓力作用下發(fā)生塑性失穩(wěn),在外壓達(dá)到該臨界壓力之前,水下分離器首先出現(xiàn)塑性變形。目標(biāo)水下分離器下放工況所承受的最大外壓為20 MPa,遠(yuǎn)小于圖算法和一致缺陷模態(tài)法分析得到的臨界失穩(wěn)壓力,那么在該壓力下分離器殼體結(jié)構(gòu)的強(qiáng)度是否滿足要求是亟需解決的問(wèn)題。水下分離器樣機(jī)筒體和嵌入式接管之間的馬鞍形焊縫是整個(gè)承壓結(jié)構(gòu)的薄弱環(huán)節(jié),須建立包含嵌入式接管的精細(xì)數(shù)值模型進(jìn)行分析。

    圖1 水下分離器一階屈曲模態(tài)Fig.1 First order buckling mode of subsea separator prototype

    圖2 筒體典型位置載荷-位移曲線Fig.2 Load-displacement curve of typical position in shell

    2 嵌入式接管子模型分析

    子模型技術(shù)又稱“切割邊界位移”技術(shù),基于圣維南原理,在獲得局部精確解的前提下,同時(shí)不增加整體模型復(fù)雜性和計(jì)算量,是有效的工程分析方法[12-13]。采用子模型技術(shù)分析水下分離器結(jié)構(gòu)強(qiáng)度的基本做法是:建立并分析整體模型(較為粗糙的模型,忽略倒角、焊縫等一些細(xì)節(jié)特征),采用切割邊界的方法,把模型局部區(qū)域從整體模型中分割出來(lái),整體模型在切割邊界上的計(jì)算位移作為子模型的位移邊界條件,然后對(duì)該局部區(qū)域細(xì)化網(wǎng)格計(jì)算,最后驗(yàn)證切割邊界的合理性。

    2.1 水下分離器整體模型分析

    圖3 水下分離器整體有限元模型分析結(jié)果Fig.3 Analysis result of subsea separator overall finite model

    2.2 嵌入式接管子模型分析

    以N1接管為例,嵌入式接管的子模型如圖4所示。馬鞍形焊縫在PRO/E中通過(guò)等厚掃描切割后經(jīng)過(guò)變截面掃描生成[14],將生成的接管模型導(dǎo)入ANSYS中進(jìn)行網(wǎng)格劃分。整個(gè)子模型單元均為六面體SOLID95單元,焊縫處網(wǎng)格進(jìn)行了細(xì)化,筒體厚度方向有8層單元。子模型單元數(shù)量接近50 000,節(jié)點(diǎn)數(shù)量超過(guò)210 000,具有較高精度。將整體模型插值提供的邊界條件施加到子模型中,并施加和整體模型一樣的載荷約束條件。焊縫經(jīng)過(guò)消氫處理和消除熱應(yīng)力熱處理,不考慮熱影響區(qū)材料性能的改變。N1接管子模型Mises應(yīng)力如圖5所示,接管應(yīng)力最大值仍然處在內(nèi)壁與筒體相交的位置,倒角的存在明顯減弱了接管和筒體處的應(yīng)力集中。焊縫應(yīng)力沿焊縫方向分布不均勻,在鞍點(diǎn)處出現(xiàn)應(yīng)力極小值。

    圖4 N1接管子模型Fig.4 Submodel of nozzle N1

    圖5 N1接管子模型Mises等效應(yīng)力云圖Fig.5 Mises stress nephogram of nozzle N1 submodel

    采用路徑映射技術(shù)[15]得到馬鞍形焊縫焊趾處應(yīng)力隨角度φ變化曲線如圖6所示,角度φ轉(zhuǎn)向定義為冠點(diǎn)逆時(shí)針至鞍點(diǎn)方向。由圖6可知,應(yīng)力沿焊縫方向分布不均勻,存在明顯的波峰和波谷。在冠點(diǎn)和鞍點(diǎn)出現(xiàn)應(yīng)力極小值,在φ為45°和135°對(duì)稱位置附近出現(xiàn)應(yīng)力極大值,這些位置是N1接管馬鞍形焊縫的應(yīng)力特征點(diǎn)。左波峰應(yīng)力極值大于右波峰,這是由于施加了軸向的加速度導(dǎo)致的。

    為驗(yàn)證切割邊界的準(zhǔn)確性,定義如圖4所示的三條切割路徑。在整體模型中定義相同的切割路徑,進(jìn)行路徑插值比對(duì),結(jié)果如圖7所示。整體模型和子模型切割邊界上應(yīng)力分布基本一致,可以認(rèn)為切割邊界的選取基本合適。

    圖6 N1接管馬鞍形焊縫焊趾應(yīng)力分布Fig.6 Stress distribution at saddle-shape weld toe of nozzle N1

    圖7 N1接管切割邊界路徑應(yīng)力分布對(duì)比Fig.7 Stress distribution comparison of cutting paths of nozzle N1

    圖8 其他接管馬鞍形焊縫焊趾應(yīng)力分布Fig.8 Stress distribution at saddle-shape weld toe of other nozzles

    由圖8可知,其他接管的應(yīng)力分布規(guī)律大致相同,接管位置對(duì)馬鞍形焊縫應(yīng)力分布基本沒(méi)有影響,兩個(gè)應(yīng)力波峰向鞍點(diǎn)漂移,最終在鞍點(diǎn)處合并成一個(gè)應(yīng)力波峰,應(yīng)力波谷消失,應(yīng)力特征點(diǎn)只有兩個(gè)冠點(diǎn)和一個(gè)鞍點(diǎn)。

    3 水下分離器高壓艙試驗(yàn)研究

    2 000m超深水水下分離器樣機(jī)接管采用法蘭密封,全部焊縫經(jīng)過(guò)X射線探傷檢測(cè)后未發(fā)現(xiàn)缺陷。為驗(yàn)證水下分離器樣機(jī)下放工況下物理承壓結(jié)構(gòu)強(qiáng)度和本文建立的精細(xì)數(shù)值模型的精度,在中船重工第七O二研究所水下工程結(jié)構(gòu)試驗(yàn)室進(jìn)行高壓艙試驗(yàn)。

    3.1 試驗(yàn)設(shè)計(jì)

    (1) 試驗(yàn)裝置。

    水下分離器高壓艙試驗(yàn)主要結(jié)構(gòu)如圖9所示。

    圖9 水下分離器高壓艙試驗(yàn)裝置Fig.9 Hyperbaric chamber experiment device of subsea separator prototype

    由于試驗(yàn)樣機(jī)總長(zhǎng)度超出壓力筒內(nèi)徑,不能將試件水平放入壓力筒內(nèi),為此定做吊耳和裙座以便于模型的吊裝及保護(hù)。水泵打壓形成試驗(yàn)壓力,應(yīng)變及壓力數(shù)據(jù)由UCAM-70A數(shù)字測(cè)量系統(tǒng)采集,現(xiàn)場(chǎng)試驗(yàn)情況如圖10所示。

    (2)測(cè)點(diǎn)布置方案。

    水下分離器承壓結(jié)構(gòu)的基本構(gòu)件為筒體和封頭,測(cè)點(diǎn)須涵蓋封頭、筒體的典型位置。封頭和筒體連接處存在較大邊緣應(yīng)力,需要設(shè)置測(cè)點(diǎn)捕捉該位置的應(yīng)力特征。由整體有限元模型分析結(jié)果可知,接管外壁和支座應(yīng)力水平較低,分離器外壁最大應(yīng)力出現(xiàn)在墊板位置,墊板位置為應(yīng)力特征點(diǎn),須布置測(cè)點(diǎn)。接管子模型分析結(jié)果顯示,N1接管馬鞍形焊縫的應(yīng)力特征點(diǎn)為冠點(diǎn)、鞍點(diǎn)以及φ=45°和135°的位置,其他接管馬鞍形焊縫應(yīng)力特征點(diǎn)為冠點(diǎn)和鞍點(diǎn)。為驗(yàn)證子模型切割邊界適用性,應(yīng)該在切割邊界設(shè)置測(cè)點(diǎn)。

    試驗(yàn)前對(duì)貼片位置進(jìn)行打磨,打磨之后使用瞬干496膠水固定應(yīng)變片,最后使用南大703膠完成測(cè)點(diǎn)的密封。測(cè)量所用雙向直角應(yīng)變片編號(hào)以筒體軸向?yàn)閱翁?hào),周向?yàn)殡p號(hào),型號(hào)為BA120-2BB;所用三向直角應(yīng)變花編號(hào)以逆時(shí)針?lè)较驗(yàn)槠鹗?型號(hào)為BA120-2CA。水下分離器試驗(yàn)樣機(jī)測(cè)點(diǎn)布置方案如圖11所示,所有應(yīng)變片貼在外表面,總共設(shè)置115片,測(cè)點(diǎn)分布覆蓋數(shù)值分析得到的應(yīng)力特征點(diǎn):應(yīng)變片E1、E2為溫度補(bǔ)償片;E3-E9為封頭典型位置測(cè)點(diǎn);E11-E14用于測(cè)量筒體封頭過(guò)渡區(qū)的邊緣應(yīng)力;E80-E87為遠(yuǎn)離結(jié)構(gòu)不連續(xù)處的筒體典型位置;E21/E22、E29/E30等測(cè)點(diǎn)位于子模型切割邊界;E108-E115分布在分離器墊板中心處;其余測(cè)點(diǎn)基本分布在嵌入式接管馬鞍形焊縫應(yīng)力特征位置;考慮到N1接管焊縫φ=45°和135°位置的對(duì)稱性,只布置一側(cè)的測(cè)點(diǎn)。

    圖10 分離器樣機(jī)現(xiàn)場(chǎng)試驗(yàn)Fig.10 Field test of subsea separator prototype

    圖11 水下分離器樣機(jī)貼片圖Fig.11 Patch figure of subsea separator prototype

    (3)加卸載程序。

    高壓艙試驗(yàn)載荷設(shè)計(jì)基于水下分離器下放到目標(biāo)水深這種最危險(xiǎn)工況:僅受20MPa外壓,內(nèi)壓為零。模型試驗(yàn)分3次加載,加載速率不大于0.5MPa/min。前兩次為預(yù)加載試驗(yàn),最高壓力為20MPa,加載程序同為:0→4→8→12→16→20→16→12→8→4→0 (MPa);每個(gè)階段穩(wěn)壓2min并進(jìn)行3次應(yīng)變測(cè)量,當(dāng)壓力升至20MPa時(shí),保壓20min并每5min進(jìn)行一次應(yīng)變測(cè)量。第3次為正式加載試驗(yàn),加壓過(guò)程不設(shè)置穩(wěn)壓臺(tái)階,緩慢加壓到設(shè)計(jì)壓力20MPa,每0.5MPa記錄數(shù)據(jù)1次;當(dāng)壓力升至20MPa時(shí),保壓20min并每5min進(jìn)行一次應(yīng)變測(cè)量;最后緩慢卸載外壓到0MPa,每0.5MPa記錄數(shù)據(jù)1次。

    3.2 試驗(yàn)結(jié)果

    試驗(yàn)過(guò)程中壓力穩(wěn)定良好,未出現(xiàn)異常響聲,試驗(yàn)結(jié)束后試件結(jié)構(gòu)未見(jiàn)明顯變形和局部屈曲失穩(wěn),法蘭管口未見(jiàn)泄漏。按要求在加卸載過(guò)程中對(duì)結(jié)構(gòu)進(jìn)行應(yīng)變測(cè)量并記錄相應(yīng)的壓力值,部分測(cè)點(diǎn)在加卸載過(guò)程中的應(yīng)變-壓力曲線如圖12所示(應(yīng)變單位為με,1με=10-6)。試驗(yàn)中水下分離器除測(cè)點(diǎn)E48以外的所有測(cè)點(diǎn)應(yīng)變-壓力關(guān)系曲線基本呈線性且卸載時(shí)回復(fù)良好,證明試驗(yàn)壓力內(nèi)各測(cè)點(diǎn)處的結(jié)構(gòu)變形處于材料彈性范圍內(nèi)。

    測(cè)點(diǎn)在20 MPa下應(yīng)變測(cè)量值分布如圖13所示。所有測(cè)點(diǎn)均為壓應(yīng)變,封頭典型位置兩向應(yīng)變基本一致,筒體典型位置測(cè)點(diǎn)的環(huán)向應(yīng)變明顯大于軸向應(yīng)變,基本上成倍數(shù)關(guān)系,這些現(xiàn)象與理論相符。應(yīng)變最大值出現(xiàn)在分離器墊板位置(對(duì)應(yīng)測(cè)點(diǎn)為E112-E113處,沿圓周方向),這與數(shù)值分析得出的結(jié)論一致。測(cè)點(diǎn)軸向應(yīng)變值在E69和E105處發(fā)生突變,應(yīng)變值急劇減小,這是因?yàn)樵撎幫搀w內(nèi)部存在整流板,對(duì)筒體結(jié)構(gòu)有一定支撐能力,軸向應(yīng)變較小。

    圖12 部分測(cè)點(diǎn)應(yīng)變-壓力曲線Fig.12 Strain-pressure curves of some measure points

    圖13 20 MPa下試驗(yàn)樣機(jī)應(yīng)變測(cè)量結(jié)果Fig.13 Strain measurement result of subsea separator prototype under 20 MPa

    4 試驗(yàn)結(jié)果與模擬結(jié)果對(duì)比分析

    根據(jù)原型樣機(jī)貼片方案,考慮貼片位置誤差(環(huán)向±1°,軸向±5mm),在數(shù)值模型中提取與試驗(yàn)樣機(jī)相對(duì)應(yīng)的56個(gè)測(cè)點(diǎn)應(yīng)變數(shù)據(jù)(焊縫處應(yīng)變結(jié)果從接管子模型中獲得)與試驗(yàn)數(shù)據(jù)進(jìn)行對(duì)比分析。

    4.1 應(yīng)力特征點(diǎn)結(jié)果對(duì)比

    筒體典型位置、N1接管以及其他接管馬鞍形焊縫應(yīng)力特征點(diǎn)結(jié)果對(duì)比如圖14所示。由圖14可知,由Lamé公式計(jì)算得到的筒體典型位置應(yīng)力理論解和模擬結(jié)果、試驗(yàn)結(jié)果基本一致,這說(shuō)明接管對(duì)遠(yuǎn)離結(jié)構(gòu)不連續(xù)處的筒體應(yīng)力場(chǎng)影響很小,只會(huì)造成周圍應(yīng)力場(chǎng)的重新分布。接管馬鞍形焊縫模擬結(jié)果和試驗(yàn)結(jié)果有較高一致性,高壓艙試驗(yàn)?zāi)軌虿蹲今R鞍形焊縫的應(yīng)力特征,測(cè)點(diǎn)布置方案合理。

    4.2 總體誤差分析

    數(shù)值模型總體誤差分析結(jié)果如圖15所示。由圖15可知,封頭、筒體典型位置以及大多數(shù)焊縫處模擬結(jié)果和試驗(yàn)結(jié)果吻合較好,一小部分焊縫測(cè)點(diǎn)誤差在20%以上,測(cè)點(diǎn)E48和E9誤差在40%以上。模擬結(jié)果誤差在10%以內(nèi)的數(shù)據(jù)占75.22%,在20%以內(nèi)的數(shù)據(jù)占93.8%,焊縫區(qū)域平均誤差為7.27%,總體平均誤差為7.93%。這說(shuō)明本文中建立的數(shù)值模型具有較高的計(jì)算精度,屬于精確建模,其分析結(jié)果能夠準(zhǔn)確模擬典型焊接承壓結(jié)構(gòu)在外壓載荷下的應(yīng)力響應(yīng),可以指導(dǎo)水下分離器結(jié)構(gòu)設(shè)計(jì)和強(qiáng)度校核。

    4.3 誤差原因分析

    數(shù)值模型采用的是對(duì)稱結(jié)構(gòu),對(duì)稱位置的應(yīng)變結(jié)果相同,E108-E115測(cè)點(diǎn)處對(duì)稱位置試驗(yàn)結(jié)果不一致,這正說(shuō)明試驗(yàn)樣機(jī)在加工制造過(guò)程中引入誤差是不可避免的。從圖12中可知測(cè)點(diǎn)E48應(yīng)變-壓力曲線不成線性,且卸載回復(fù)效果較差,說(shuō)明該測(cè)點(diǎn)處應(yīng)變片測(cè)得數(shù)據(jù)可信度不高,造成模擬結(jié)果相對(duì)誤差較大。測(cè)點(diǎn)E9靠近封頭直線邊,封頭在沖壓成型過(guò)程中會(huì)造成局部減薄和厚度增加,在靠近直線邊封頭厚度增加量最大[16],在現(xiàn)場(chǎng)測(cè)量中該處的厚度達(dá)到87mm以上,造成該測(cè)點(diǎn)處應(yīng)變減小。

    5 結(jié) 論

    (1) 首次提出的2 000m超深水水下分離器下放工況承壓性能評(píng)估方法在分析結(jié)構(gòu)極限承載能力的基礎(chǔ)上,確定了下放工況應(yīng)首先考慮結(jié)構(gòu)強(qiáng)度問(wèn)題。建立的精細(xì)數(shù)值模型及設(shè)計(jì)的高壓艙試驗(yàn)檢驗(yàn)了數(shù)模-物模一致性,數(shù)模結(jié)果指導(dǎo)物模試驗(yàn)設(shè)計(jì),物模試驗(yàn)驗(yàn)證了數(shù)模精度。

    (2)N1接管馬鞍形焊縫在φ=45°和φ=135°位置存在兩個(gè)應(yīng)力波峰,鞍點(diǎn)和冠點(diǎn)為應(yīng)力極小值位置;樣機(jī)接管位置不影響焊縫應(yīng)力分布,其他接管焊縫應(yīng)力分布相同,只在鞍點(diǎn)位置存在一個(gè)應(yīng)力波峰,冠點(diǎn)為應(yīng)力極小值位置。

    (3)自主設(shè)計(jì)的超深水水下分離器整體應(yīng)力水平處于線彈性變形階段,能夠滿足下放危險(xiǎn)工況下的強(qiáng)度要求。模擬結(jié)果總體平均相對(duì)誤差為7.93%,數(shù)值模型具有較高精度。

    [1] URTHALER Y, MCFARLANE R, BROOKES D. Effect of seawater pressure on the static and fatigue strength of subsea production equipment used in deepwater oil and gas fields: proceedings of the ASME Pressure Vessels and Piping Conference 2009, Czech republic[C]. Prague: American Society of Mechanical Engineers, 2010.

    [2] AN C, CASTELLO X, OLIVEIRA A M, et al. Limit strength of new Sandwich Pipes with strain hardening cementitious composites (SHCC) core: finite element modelling: ASME 2012: 31st International Conference on Ocean, Offshore and Arctic Engineering, Rio de Janeiro[C]. Brazil: American Society of Mechanical Engineers, 2012.

    [3] LI H, PAN X, REN H, et al. Study on the residual strength assessment method on corroded subsea pipeline based on reference stress method: proceedings of the 28th International Conference on Ocean, Offshore and Arctic Engineering 2009, Honolulu[C].Hawaii: American Society of Mechanical Engineers, 2009.

    [4] NETTO T A, FERRAZ U S, BOTTO A. Residual strength of corroded pipelines under external pressure: a simple assessment: 6th International Pipeline Conference, Calgary[C]. Alberta: American Society of Mechanical Engineers, 2007.[5] LEE G H, POURARIA H, SEO J K, et al. Burst strength behaviour of an aging subsea gas pipeline elbow in different external and internal corrosion-damaged positions[J]. International Journal of Naval Architecture and Ocean Engineering, 2015,7(3):435-451.[6] ZHENG J, PALMER A, BRUNNING P, et al. Indentation and external pressure on subsea single wall pipe and pipe-in-pipe[J]. Ocean Engineering, 2014,83:125-132.

    [7] CHOO Y S, HORNE M R, MONTAGUE P. Strength of intersections in composite sandwich vessels under external pressure: proceedings of the 1st International Offshore and Polar Engineering Conference[C]. International Society of Offshore and Polar Engineers (ISOPE), 1991.

    [8] 葛玖浩,李偉,陳國(guó)明,等.2 000 m超深水水下分離器承壓結(jié)構(gòu)強(qiáng)度分析[J].石油機(jī)械,2015,43(2):60-64. GE Jiuhao, LI Wei, CHEN Guoming, et al. Strength analysis on bearing structure of subsea separator for 2 000 m ultra deepwater[J]. China Petroleum Machinery, 2015,43(2):60-64.

    [9] 李志剛,劉培林,高杰,等.基于分析設(shè)計(jì)直接法的水下分離器結(jié)構(gòu)總體塑性變形校核[J].中國(guó)石油大學(xué)學(xué)報(bào)(自然科學(xué)版),2013,37(3):137-140. LI Zhigang, LIU Peilin, GAO Jie, et al. Gross plastic deformation design check for subsea separator based on direct route method for design by analysis[J]. Journal of China University of Petroleum (Edition of Natural Science), 2013,37(3):137-140.

    [10] 壓力容器: GB 150-2011[S].北京:中國(guó)標(biāo)準(zhǔn)出版社,2011.

    [11] 鄧志軍,陳冰冰,鄭浣琪,等.外壓模擬計(jì)算中圓筒初始幾何偏差描述方法的研究[J].機(jī)械工程學(xué)報(bào),2015,51(6):66-72. DENG Zhijun, CHEN Bingbing, ZHENG Huanqi, et al. Research on initial geometric deviation description for numerical simulation of cylindrical shells under external pressure[J]. Journal of Mechanical Engineering, 2015,51(6):66-72.

    [12] ARAI K, YODO K, OKADA H, et al. Ultra-large scale fracture mechanics analysis using a parallel finite element method with submodel technique[J]. Finite Elements in Analysis and Design, 2015,105:44-55.

    [13] WANG H, LI A, HU R, et al. Accurate stress analysis on steel box girder of long span suspension bridges based on multi-scale submodeling method[J]. Advances in Structural Engineering, 2010,13(4):727-740.

    [14] 陳永秋,白雪.帶鈍邊馬鞍形焊接坡口數(shù)控切割數(shù)學(xué)模型[J].焊接學(xué)報(bào),2010,31(7):91-94. CHEN Yongqiu, BAI Xue. Mathematical model for NC cutting saddle type of welding groove with edge[J]. Transactions of the China Welding Institution, 2010,31(7):91-94.

    [15] 陳團(tuán)海,陳國(guó)明.T型焊接管節(jié)點(diǎn)應(yīng)力集中系數(shù)數(shù)值分析[J].焊接學(xué)報(bào),2010,31(11):45-48. CHEN Tuanhai, CHEN Guoming. Numerical analysis on stress concentration factor of joints in welded T-tube[J]. Transactions of the China Welding Institution, 2010,31(11):45-48.

    [16] 夏春申,徐貴珍.厚壁半球封頭及大型中厚橢圓封頭沖壓工藝探討[J].壓力容器,1997(4):72-75. XIA Chunshen, XU Guizhen. A study on press of process for thick-walled hemispherical head and large-scale medium thick-walled elliptical head[J]. Pressure Vessel Technology, 1997(4):72-75.

    (編輯 沈玉英)

    Numerical analysis and experimental study on pressure bearing performance of ultra-deepwater subsea separator

    LI Xiumei, CHEN Guoming, LI Wei, ZHANG Shenyan, GE Jiuhao, LIU Xiuquan

    (CentreforOffshoreEngineeringandSafetyTechnologyinChinaUniversityofPetroleum,Qingdao266580,China)

    With an independently designed 2 000 m subsea separator as the research object, the ultimate bearing capacity of the separator prototype was analyzed based on the nomography method and the conformable imperfection mode method. A fine subsea separator numerical model containing embedded nozzle submodel was established to explore the stress distribution of saddle-shape weld. A hyperbaric chamber experiment was designed to validate the pressure-bearing performance of the subsea separator, and the experiment data were compared with the simulation data. Research results indicate that the instability mode of the separator is plastic instability and the structural strength should be assessed in the installation condition. The stress-characteristic locations are at the coronal and saddle points for saddle-shape weld in all nozzles and in the position ofφ=45° andφ=135° for saddle-shape weld in nozzle N1. In the experiment, the overall stress level of the subsea separator is in an elastic deformation stage and the prototype can satisfy the pressure-bearing strength requirements at 2 000 m water depth. The numerical analysis results are highly consistent with the experimental results.

    ultra-deepwater subsea separator; embedded nozzle; submodel; saddle-shape weld; hyperbaric chamber experiment

    2016-04-22

    國(guó)家高技術(shù)研究發(fā)展計(jì)劃項(xiàng)目(2013AA09A213);中央高?;究蒲袠I(yè)務(wù)費(fèi)專項(xiàng) (14CX06123A);山東省自然科學(xué)基金聯(lián)合專項(xiàng)(ZR2014EL018)

    李秀美(1989-),男,博士研究生,研究方向?yàn)楹Q蟪袎航Y(jié)構(gòu)完整性管理與評(píng)價(jià)。E-mail:lixiumei2006 @126.com。

    陳國(guó)明(1962-),男,教授,博士,博士生導(dǎo)師,研究方向?yàn)楹Q笥蜌夤こ碳把b備、油氣安全工程。E-mail: offshore@126.com。

    1673-5005(2016)06-0126-09

    10.3969/j.issn.1673-5005.2016.06.016

    TE 58

    A

    李秀美,陳國(guó)明,李偉,等. 超深水水下分離器承壓性能數(shù)值分析與試驗(yàn)研究[J]. 中國(guó)石油大學(xué)學(xué)報(bào)(自然科學(xué)版), 2016,40(6):126-134.

    LI Xiumei, CHEN Guoming, LI Wei, et al. Numerical analysis and experimental study on pressure bearing performance of ultra-deepwater subsea separator[J]. Journal of China University of Petroleum (Edition of Natural Science), 2016,40(6):126-134.

    猜你喜歡
    樣機(jī)筒體分離器
    基于ADAMS虛擬樣機(jī)的門座起重機(jī)動(dòng)力學(xué)仿真研究
    b型管板與筒體溫差應(yīng)力的分析計(jì)算和評(píng)定
    化工管理(2021年7期)2021-05-13 00:46:04
    回轉(zhuǎn)窯筒體對(duì)接操作方法
    一種臥式筒體糞污發(fā)酵裝置的筒體設(shè)計(jì)與分析
    基于三維數(shù)字樣機(jī)的運(yùn)載器裝配工藝規(guī)劃技術(shù)
    強(qiáng)制循環(huán)流化床物料分離效率的研究及應(yīng)用
    旋風(fēng)分離器切向進(jìn)料接管焊縫修理與改進(jìn)建議
    苜蓿蠕變特性的虛擬樣機(jī)仿真研究
    HAZOP理論分析在分離器作業(yè)中的應(yīng)用探討
    安全(2015年1期)2016-01-19 06:19:07
    便攜式正滲透凈水袋原理樣機(jī)的耐用性能研究
    大型av网站在线播放| 后天国语完整版免费观看| 国产国语露脸激情在线看| 午夜精品在线福利| 一个人免费在线观看的高清视频| 少妇裸体淫交视频免费看高清 | 视频区图区小说| 中文字幕色久视频| 极品人妻少妇av视频| 99热国产这里只有精品6| 欧美日韩福利视频一区二区| 超色免费av| 在线天堂中文资源库| 欧美黑人欧美精品刺激| 人妻一区二区av| 黄片小视频在线播放| 巨乳人妻的诱惑在线观看| 精品国产亚洲在线| 777久久人妻少妇嫩草av网站| 亚洲自偷自拍图片 自拍| 亚洲avbb在线观看| 在线天堂中文资源库| 在线视频色国产色| 在线观看免费高清a一片| 国产精品亚洲av一区麻豆| 国产熟女午夜一区二区三区| 搡老乐熟女国产| 九色亚洲精品在线播放| 欧美黄色片欧美黄色片| av超薄肉色丝袜交足视频| 国产成人欧美在线观看 | 19禁男女啪啪无遮挡网站| 亚洲五月婷婷丁香| 精品国内亚洲2022精品成人 | 男女下面插进去视频免费观看| 女人被狂操c到高潮| 在线国产一区二区在线| 色在线成人网| 99国产综合亚洲精品| 午夜老司机福利片| 又黄又粗又硬又大视频| 在线观看免费视频网站a站| 久久久久久久精品吃奶| 国产男女内射视频| 麻豆成人av在线观看| 午夜福利乱码中文字幕| 黄片大片在线免费观看| 精品国内亚洲2022精品成人 | 18在线观看网站| 夫妻午夜视频| 亚洲专区国产一区二区| 亚洲成av片中文字幕在线观看| 久久精品亚洲av国产电影网| 中文字幕av电影在线播放| 久久久久久久国产电影| 国产一卡二卡三卡精品| a级片在线免费高清观看视频| 国产高清激情床上av| 日本五十路高清| 黄色毛片三级朝国网站| 久久 成人 亚洲| 国产单亲对白刺激| 操出白浆在线播放| bbb黄色大片| 三级毛片av免费| 成人黄色视频免费在线看| 性色av乱码一区二区三区2| 精品国产乱码久久久久久男人| 国产成人系列免费观看| 国产在线一区二区三区精| 欧美久久黑人一区二区| 国产成人av激情在线播放| 看片在线看免费视频| 国内久久婷婷六月综合欲色啪| av网站在线播放免费| 国产av又大| 欧美成狂野欧美在线观看| 在线观看免费午夜福利视频| 国内久久婷婷六月综合欲色啪| 久久国产亚洲av麻豆专区| 50天的宝宝边吃奶边哭怎么回事| 色尼玛亚洲综合影院| 久久人妻av系列| 国产成人精品久久二区二区免费| 麻豆成人av在线观看| 少妇粗大呻吟视频| 少妇粗大呻吟视频| 又紧又爽又黄一区二区| 久久精品亚洲熟妇少妇任你| 久久草成人影院| 免费看十八禁软件| 美女高潮喷水抽搐中文字幕| 999久久久国产精品视频| 水蜜桃什么品种好| 在线观看免费视频日本深夜| 久久久久久久国产电影| 美女午夜性视频免费| 天堂√8在线中文| 精品福利永久在线观看| 免费在线观看黄色视频的| www.精华液| 在线观看舔阴道视频| 极品少妇高潮喷水抽搐| 夜夜爽天天搞| 国产激情欧美一区二区| 免费黄频网站在线观看国产| 久久亚洲真实| av国产精品久久久久影院| 国产高清视频在线播放一区| 久久人妻熟女aⅴ| e午夜精品久久久久久久| 我的亚洲天堂| 黑人操中国人逼视频| 美国免费a级毛片| 一级毛片女人18水好多| 最近最新中文字幕大全电影3 | 一级a爱片免费观看的视频| 日本黄色视频三级网站网址 | 精品国产美女av久久久久小说| 久久婷婷成人综合色麻豆| 美国免费a级毛片| 一本大道久久a久久精品| 成熟少妇高潮喷水视频| 亚洲avbb在线观看| 热re99久久国产66热| 女人高潮潮喷娇喘18禁视频| 高清黄色对白视频在线免费看| 欧洲精品卡2卡3卡4卡5卡区| 悠悠久久av| 午夜福利免费观看在线| 久久久久久亚洲精品国产蜜桃av| 久热爱精品视频在线9| 亚洲男人天堂网一区| 黄色片一级片一级黄色片| 9191精品国产免费久久| av视频免费观看在线观看| 久久性视频一级片| 美国免费a级毛片| 久久精品aⅴ一区二区三区四区| 黄片大片在线免费观看| 欧美乱色亚洲激情| 欧美黑人精品巨大| 久久人妻熟女aⅴ| 国产色视频综合| 男人的好看免费观看在线视频 | 777米奇影视久久| 中文字幕人妻熟女乱码| 亚洲一码二码三码区别大吗| 99国产精品一区二区蜜桃av | 亚洲一区高清亚洲精品| 国产精品偷伦视频观看了| 两性夫妻黄色片| 亚洲精华国产精华精| 久久中文字幕人妻熟女| 亚洲人成77777在线视频| 黄色视频,在线免费观看| 日韩视频一区二区在线观看| 欧美黑人欧美精品刺激| 欧美激情极品国产一区二区三区| 久久香蕉国产精品| 久久香蕉国产精品| 亚洲精品一二三| 国产xxxxx性猛交| tube8黄色片| 人人妻,人人澡人人爽秒播| av欧美777| 啦啦啦 在线观看视频| 亚洲精品自拍成人| 99久久国产精品久久久| 在线观看66精品国产| 91精品三级在线观看| 国产亚洲欧美98| 午夜激情av网站| 国产aⅴ精品一区二区三区波| 色老头精品视频在线观看| a级片在线免费高清观看视频| 18禁黄网站禁片午夜丰满| 亚洲一区中文字幕在线| 欧美午夜高清在线| e午夜精品久久久久久久| 两性午夜刺激爽爽歪歪视频在线观看 | 黑人欧美特级aaaaaa片| 建设人人有责人人尽责人人享有的| 美女午夜性视频免费| 人人妻人人澡人人爽人人夜夜| 中亚洲国语对白在线视频| 欧美日韩一级在线毛片| 久久草成人影院| 国产男女内射视频| 91大片在线观看| 国产精品乱码一区二三区的特点 | 日本黄色视频三级网站网址 | 成在线人永久免费视频| 久久精品国产亚洲av高清一级| 91av网站免费观看| 国产亚洲欧美在线一区二区| 日韩精品免费视频一区二区三区| 国产av一区二区精品久久| 黄色女人牲交| 欧美午夜高清在线| 丰满的人妻完整版| 人妻 亚洲 视频| 一区福利在线观看| 亚洲精品中文字幕在线视频| 岛国在线观看网站| 成人永久免费在线观看视频| 一边摸一边抽搐一进一出视频| 国产精品国产高清国产av | av中文乱码字幕在线| 超碰97精品在线观看| 老司机亚洲免费影院| 成年人午夜在线观看视频| 久久精品亚洲精品国产色婷小说| 9191精品国产免费久久| 淫妇啪啪啪对白视频| 国产高清国产精品国产三级| 色尼玛亚洲综合影院| 日日夜夜操网爽| 欧美精品啪啪一区二区三区| 国产深夜福利视频在线观看| avwww免费| 男女之事视频高清在线观看| 国产欧美亚洲国产| 日韩熟女老妇一区二区性免费视频| 啦啦啦 在线观看视频| 俄罗斯特黄特色一大片| 成人黄色视频免费在线看| 欧美成狂野欧美在线观看| 欧美精品一区二区免费开放| tocl精华| 国产亚洲欧美98| 久久久国产成人精品二区 | 亚洲av欧美aⅴ国产| 天天操日日干夜夜撸| 久久人妻熟女aⅴ| 少妇裸体淫交视频免费看高清 | 成在线人永久免费视频| 久久久久精品人妻al黑| 国产精品av久久久久免费| 日本黄色日本黄色录像| 日本精品一区二区三区蜜桃| 久久香蕉国产精品| 中文字幕另类日韩欧美亚洲嫩草| 国产有黄有色有爽视频| 亚洲aⅴ乱码一区二区在线播放 | 少妇猛男粗大的猛烈进出视频| 窝窝影院91人妻| 久久人人爽av亚洲精品天堂| 91在线观看av| 午夜激情av网站| 老司机在亚洲福利影院| 老司机在亚洲福利影院| 久久婷婷成人综合色麻豆| 人人妻,人人澡人人爽秒播| 亚洲国产毛片av蜜桃av| 在线免费观看的www视频| 99在线人妻在线中文字幕 | 国产成人一区二区三区免费视频网站| 久久久久久亚洲精品国产蜜桃av| 精品久久久久久久久久免费视频 | 国产精品av久久久久免费| 黄色视频,在线免费观看| 中文字幕人妻丝袜制服| 狠狠婷婷综合久久久久久88av| 久热这里只有精品99| 久热这里只有精品99| 99热只有精品国产| 国产av精品麻豆| 69av精品久久久久久| 久久午夜亚洲精品久久| 一进一出抽搐动态| 成人国语在线视频| 中文字幕最新亚洲高清| 99精品在免费线老司机午夜| 久久精品亚洲熟妇少妇任你| 久久国产精品影院| 国产精品综合久久久久久久免费 | av中文乱码字幕在线| 美女午夜性视频免费| 在线观看免费视频网站a站| 80岁老熟妇乱子伦牲交| 亚洲av成人av| 免费女性裸体啪啪无遮挡网站| 99国产综合亚洲精品| 51午夜福利影视在线观看| 国产伦人伦偷精品视频| 99国产精品一区二区蜜桃av | 免费一级毛片在线播放高清视频 | 日韩制服丝袜自拍偷拍| 日韩制服丝袜自拍偷拍| 国产真人三级小视频在线观看| 亚洲欧美色中文字幕在线| 国产一卡二卡三卡精品| 日韩视频一区二区在线观看| 久久久久国内视频| 99香蕉大伊视频| 色综合欧美亚洲国产小说| 9色porny在线观看| 韩国av一区二区三区四区| 色婷婷久久久亚洲欧美| 老熟妇乱子伦视频在线观看| 一边摸一边做爽爽视频免费| 黄色成人免费大全| 好男人电影高清在线观看| 黄色视频不卡| 国产精品亚洲av一区麻豆| 啦啦啦视频在线资源免费观看| 成年人午夜在线观看视频| 一区二区三区激情视频| 天天躁夜夜躁狠狠躁躁| 国产真人三级小视频在线观看| 亚洲精品久久成人aⅴ小说| 91av网站免费观看| 日本撒尿小便嘘嘘汇集6| 脱女人内裤的视频| 国产乱人伦免费视频| 国产亚洲一区二区精品| 国产深夜福利视频在线观看| 精品人妻在线不人妻| 久久人妻福利社区极品人妻图片| 免费观看a级毛片全部| 国产激情久久老熟女| 一级毛片精品| 久久久精品区二区三区| 精品乱码久久久久久99久播| 制服诱惑二区| 国产人伦9x9x在线观看| 欧美日韩国产mv在线观看视频| 国产无遮挡羞羞视频在线观看| 满18在线观看网站| 性少妇av在线| 成人精品一区二区免费| 国产亚洲精品久久久久久毛片 | а√天堂www在线а√下载 | 久久中文看片网| 国产精品国产高清国产av | 午夜免费鲁丝| 侵犯人妻中文字幕一二三四区| 国产精品美女特级片免费视频播放器 | 欧美久久黑人一区二区| 美国免费a级毛片| 男女午夜视频在线观看| 国产男女超爽视频在线观看| 国产精品亚洲av一区麻豆| 久久国产亚洲av麻豆专区| 亚洲专区中文字幕在线| 99精品欧美一区二区三区四区| 黄片大片在线免费观看| 国产aⅴ精品一区二区三区波| 亚洲精品在线观看二区| 午夜老司机福利片| 黑人巨大精品欧美一区二区蜜桃| 国产极品粉嫩免费观看在线| 国产精品综合久久久久久久免费 | 欧美黑人欧美精品刺激| 看免费av毛片| 大香蕉久久网| 国产成人精品久久二区二区免费| videosex国产| 色播在线永久视频| 天天添夜夜摸| 美国免费a级毛片| 99国产综合亚洲精品| 真人做人爱边吃奶动态| 久久久久久久午夜电影 | 亚洲全国av大片| 麻豆国产av国片精品| 精品国产超薄肉色丝袜足j| 日韩大码丰满熟妇| 动漫黄色视频在线观看| 中文欧美无线码| 9热在线视频观看99| 男女免费视频国产| 老司机亚洲免费影院| 亚洲视频免费观看视频| 亚洲精品乱久久久久久| 日韩成人在线观看一区二区三区| 热99re8久久精品国产| 在线观看66精品国产| 亚洲av第一区精品v没综合| 91老司机精品| 国产日韩欧美亚洲二区| 亚洲精品中文字幕一二三四区| 天天添夜夜摸| 国产精品 欧美亚洲| 久久久久久免费高清国产稀缺| 欧美午夜高清在线| 99热国产这里只有精品6| 亚洲中文日韩欧美视频| 他把我摸到了高潮在线观看| 国产97色在线日韩免费| 黑人巨大精品欧美一区二区mp4| 欧美激情极品国产一区二区三区| 老司机福利观看| 欧美亚洲日本最大视频资源| 日本a在线网址| 少妇猛男粗大的猛烈进出视频| 日韩精品免费视频一区二区三区| 日韩有码中文字幕| 欧美日韩精品网址| 岛国毛片在线播放| av网站在线播放免费| 女人被躁到高潮嗷嗷叫费观| 精品免费久久久久久久清纯 | 女人高潮潮喷娇喘18禁视频| 99riav亚洲国产免费| 咕卡用的链子| 一本一本久久a久久精品综合妖精| 久久久久国产一级毛片高清牌| 妹子高潮喷水视频| 亚洲精品中文字幕一二三四区| 欧美最黄视频在线播放免费 | 日本黄色视频三级网站网址 | netflix在线观看网站| 久久天堂一区二区三区四区| 欧美日韩一级在线毛片| 女性被躁到高潮视频| 91精品三级在线观看| 精品午夜福利视频在线观看一区| 动漫黄色视频在线观看| 国产有黄有色有爽视频| 久久人人爽av亚洲精品天堂| 亚洲精品在线观看二区| 国产精品偷伦视频观看了| 亚洲专区字幕在线| 深夜精品福利| 夜夜夜夜夜久久久久| 亚洲精品中文字幕在线视频| 亚洲av日韩在线播放| 久久人人爽av亚洲精品天堂| 精品亚洲成国产av| 青草久久国产| 岛国在线观看网站| 五月开心婷婷网| 99re在线观看精品视频| 男人的好看免费观看在线视频 | 黄片大片在线免费观看| 黑人操中国人逼视频| av欧美777| 黄色片一级片一级黄色片| 欧美激情久久久久久爽电影 | 亚洲av片天天在线观看| 午夜福利乱码中文字幕| 国产欧美日韩一区二区精品| 99国产精品一区二区三区| 日韩欧美国产一区二区入口| 黄频高清免费视频| 国产高清videossex| 成年人午夜在线观看视频| 欧美激情极品国产一区二区三区| 午夜成年电影在线免费观看| 亚洲精品在线观看二区| 成年动漫av网址| 男女免费视频国产| 在线观看免费高清a一片| 久久精品国产亚洲av香蕉五月 | 亚洲三区欧美一区| 高清黄色对白视频在线免费看| 亚洲 欧美一区二区三区| 少妇被粗大的猛进出69影院| 18禁黄网站禁片午夜丰满| 亚洲熟女毛片儿| 久久人妻福利社区极品人妻图片| 精品国内亚洲2022精品成人 | av网站免费在线观看视频| 又黄又粗又硬又大视频| 久久久国产精品麻豆| 我的亚洲天堂| 久久精品aⅴ一区二区三区四区| 精品久久久久久电影网| 麻豆乱淫一区二区| 欧美日韩乱码在线| 国产99久久九九免费精品| 日韩欧美一区视频在线观看| 黄色片一级片一级黄色片| 在线观看免费视频日本深夜| 久久中文字幕一级| 欧美最黄视频在线播放免费 | 麻豆成人av在线观看| 久久人妻熟女aⅴ| 日日摸夜夜添夜夜添小说| 亚洲,欧美精品.| 免费久久久久久久精品成人欧美视频| 久久精品国产a三级三级三级| 久久精品亚洲av国产电影网| 激情在线观看视频在线高清 | 91在线观看av| 亚洲精品国产精品久久久不卡| 多毛熟女@视频| 国产在线观看jvid| 一级作爱视频免费观看| 亚洲精品成人av观看孕妇| 天堂中文最新版在线下载| 成人精品一区二区免费| 亚洲第一青青草原| 久久久国产一区二区| 女人高潮潮喷娇喘18禁视频| 99香蕉大伊视频| 黄网站色视频无遮挡免费观看| 啪啪无遮挡十八禁网站| 在线观看舔阴道视频| 99精品欧美一区二区三区四区| 黑人猛操日本美女一级片| 村上凉子中文字幕在线| 在线观看日韩欧美| 国产精品久久久久久精品古装| 国产精品久久久久久人妻精品电影| 真人做人爱边吃奶动态| 捣出白浆h1v1| 中文字幕最新亚洲高清| 国产高清激情床上av| 中出人妻视频一区二区| 精品亚洲成a人片在线观看| 亚洲人成伊人成综合网2020| 纯流量卡能插随身wifi吗| 久久性视频一级片| 国产乱人伦免费视频| 性色av乱码一区二区三区2| 亚洲av片天天在线观看| 国产欧美亚洲国产| 高潮久久久久久久久久久不卡| 国产av一区二区精品久久| 色综合欧美亚洲国产小说| av中文乱码字幕在线| 日韩有码中文字幕| 久久久久久人人人人人| 青草久久国产| 天堂俺去俺来也www色官网| a级毛片在线看网站| 国产精品98久久久久久宅男小说| 国产亚洲精品久久久久5区| 韩国av一区二区三区四区| 午夜福利在线免费观看网站| 黄色怎么调成土黄色| 高清视频免费观看一区二区| 成人18禁高潮啪啪吃奶动态图| 热99久久久久精品小说推荐| 欧美精品av麻豆av| 另类亚洲欧美激情| 亚洲avbb在线观看| 国产xxxxx性猛交| 久久久久国产一级毛片高清牌| 亚洲国产毛片av蜜桃av| 日本撒尿小便嘘嘘汇集6| av天堂在线播放| 精品无人区乱码1区二区| 老熟女久久久| 国产欧美日韩综合在线一区二区| 女人被躁到高潮嗷嗷叫费观| 中文字幕人妻熟女乱码| 亚洲精品中文字幕在线视频| 欧美成狂野欧美在线观看| 不卡一级毛片| 国产色视频综合| 不卡一级毛片| 欧美日韩一级在线毛片| 一进一出抽搐gif免费好疼 | 久9热在线精品视频| 欧美日韩黄片免| 男人操女人黄网站| 久久婷婷成人综合色麻豆| 一区在线观看完整版| 久久久久久免费高清国产稀缺| 香蕉国产在线看| 亚洲精品国产区一区二| 成人手机av| 超碰97精品在线观看| 免费观看a级毛片全部| 日韩免费av在线播放| 99久久国产精品久久久| 亚洲伊人色综图| 少妇的丰满在线观看| 久久国产精品男人的天堂亚洲| 国产亚洲欧美精品永久| 天天影视国产精品| 大码成人一级视频| 亚洲国产看品久久| 国产男女内射视频| 久久亚洲精品不卡| av天堂久久9| 人人妻人人爽人人添夜夜欢视频| 黑人欧美特级aaaaaa片| 日韩三级视频一区二区三区| 99精品在免费线老司机午夜| 成人免费观看视频高清| 亚洲午夜理论影院| а√天堂www在线а√下载 | 男女下面插进去视频免费观看| 岛国在线观看网站| 亚洲一区二区三区欧美精品| 动漫黄色视频在线观看| 中文字幕人妻熟女乱码| 乱人伦中国视频| av片东京热男人的天堂| 免费少妇av软件| www.熟女人妻精品国产| 亚洲精品一二三| 制服人妻中文乱码| 久久久久精品人妻al黑| 成人精品一区二区免费| 99精国产麻豆久久婷婷| 欧美乱色亚洲激情| 亚洲一码二码三码区别大吗| 水蜜桃什么品种好| 日韩大码丰满熟妇| 精品电影一区二区在线| 两人在一起打扑克的视频| 欧美在线黄色| svipshipincom国产片| 亚洲第一青青草原| 久久国产亚洲av麻豆专区| 国产精品av久久久久免费| 国产免费男女视频| 久久99一区二区三区| 亚洲三区欧美一区|