深海管道非線性屈曲理論計(jì)算方法

余建星，李智博，杜尊峰，傅明煬，卞雪航，楊 源

(天津大學(xué)建筑工程學(xué)院，天津 300072)

在深水環(huán)境條件下，管道一旦出現(xiàn)結(jié)構(gòu)失效，將對工程建設(shè)和維修造成重大損失。由于管道在生產(chǎn)制造過程中，其幾何尺寸和材料性質(zhì)的分布存在一定程度的不均勻性，形成管道結(jié)構(gòu)的局部缺陷，使管道在一定壓力條件下發(fā)生屈曲失穩(wěn)。在管道的鋪設(shè)和使用過程中，如果發(fā)生碰撞或腐蝕等現(xiàn)象，將進(jìn)一步增加管道局部缺陷的程度，降低管道屈曲失穩(wěn)所需的壓潰壓力。對于大徑厚比情況，管道處于小變形的彈性階段發(fā)生屈曲失穩(wěn)，此時(shí)的幾何變形和材料本構(gòu)皆為線性關(guān)系。然而，對于小徑厚比情況而言，發(fā)生屈曲失穩(wěn)時(shí)幾何變形和材料本構(gòu)皆可能處于非線性階段。對于存在初始缺陷的情況，其與徑厚比的耦合影響，更使得管道的屈曲失穩(wěn)表現(xiàn)為明顯的非線性特征。

周承倜［1］對海底管道的屈曲進(jìn)行了研究，介紹了8種不同控制屈曲發(fā)生的屈曲壓力。Fabian［2］對彈性范圍內(nèi)管道在彎矩和軸向拉力聯(lián)合作用下的非線性屈曲響應(yīng)進(jìn)行過分析。Wilhoit［3］等研究了管道在彎矩和軸向拉力聯(lián)合作用下的彈塑性響應(yīng)，但未考慮加載路徑的影響。Palmer［4］于1975年首次揭示深水管道鋪設(shè)過程中，管道的局部屈曲失穩(wěn)會(huì)導(dǎo)致屈曲沿管道軸線傳播的現(xiàn)象。Bai［5］應(yīng)用有限元軟件分析了管道在外壓、軸向拉力和彎矩共同作用下的屈曲性能，并進(jìn)行了初始橢圓度、殘余應(yīng)力、應(yīng)變硬化和加載路徑的敏感性分析。Dama［6］通過試驗(yàn)和數(shù)值模擬研究了受損管道在內(nèi)壓和彎矩作用下的承載力，以及在交變荷載作用下的疲勞剩余壽命預(yù)測，為局部屈曲受損管道的繼續(xù)使用提供了安全評估基礎(chǔ)。Toscano［7］等人建立了三維有限元模型，分析了深水海底管道整體式止屈器的穿越壓力和機(jī)理，并與試驗(yàn)結(jié)果進(jìn)行對比。Xue［8］應(yīng)用ABAQUS軟件分析了非均勻和受腐蝕海底管道在外壓作用下的屈曲和屈曲傳播特性，發(fā)現(xiàn)了管道對稱和反對稱屈曲模態(tài)主要與管道截面不均勻部分的相對厚度和角度有關(guān)。Estefen［9］通過小比例試驗(yàn)結(jié)合薄殼理論分析評估了完好管道在外部靜水壓力和彎矩作用下的極限強(qiáng)度，以及受損管道在外部靜水壓力作用下的剩余強(qiáng)度。Kyriakide.S［10］等對管道的屈曲問題進(jìn)行了研究，分析了徑厚比、材料特性、初始橢圓度和加載路徑對管道局部屈曲的影響。綜上所述，在深海管道屈曲方面，國內(nèi)研究較少，且較國外尚有較大差距。針對海底管道屈曲過程中的幾何非線性和材料非線性因素，建立屈曲準(zhǔn)則，利用虛功原理求解含橢圓度缺陷的管道屈曲變形及其臨界壓力，分析其影響規(guī)律，并進(jìn)行全比例試驗(yàn)和用有限元模擬進(jìn)行驗(yàn)證，取得很好結(jié)果。此分析方法為我國深海油氣開發(fā)的安全保障提供理論支持，為進(jìn)一步深入研究屈曲傳播和止屈技術(shù)奠定基礎(chǔ)。

1 理論模型

管道在外壓作用下發(fā)生橢圓度變形。對于沿軸向具有缺陷程度近乎相同的較長管道而言，在其前屈曲階段，這種橢圓度變形沿軸向幾乎是一致的，直至缺陷程度略微大于其它截面的部位首先發(fā)生屈曲失穩(wěn)為止。管道壓潰壓力的確定，與失穩(wěn)之后的后屈曲階段無關(guān)，僅取決于軸向變形一致的前屈曲階段，可以以某一截面的變形代替管道整體的變形。故此，建立二維環(huán)模型，如圖1所示，并基于以下假定:

1)管壁厚度沿環(huán)向均勻分布;

2)圓環(huán)的彎曲變形滿足平斷面假定;

3)圓環(huán)徑向應(yīng)力與其環(huán)向應(yīng)力相比較小，可以忽略不計(jì);

4)考慮材料的彈塑性本構(gòu)關(guān)系;

5)認(rèn)為在變形過程中，材料保持各向同性;

6)管道長度足夠長，端部約束條件的影響可以忽略;

7)考慮截面上存在沿軸向一致的軸向應(yīng)力和軸向應(yīng)變的影響;

8)管道壓潰時(shí)，截面發(fā)生雙凹形變形。

坐標(biāo)系的建立包括整體笛卡爾坐標(biāo)系和局部曲線坐標(biāo)系。曲線坐標(biāo)系取圓環(huán)變形前的構(gòu)形作為參考，v，w分別表示環(huán)向和徑向的位移分量，R為變形前的中性面的半徑，管壁的中性面上各點(diǎn)由坐標(biāo)R，θ表示，如圖2所示。

圖1 管道變形截面示意Fig.1 Cross-section pipe deformation

圖2 中性面位移的定義Fig.2 Definition of the neutral plane displacement

2 基本方程

2.1 幾何方程

將環(huán)向應(yīng)變?chǔ)?分為薄膜應(yīng)變和彎曲應(yīng)變 κθθ兩部分的疊加:

其中，z為中性層法向坐標(biāo)(m)。薄膜應(yīng)變和彎曲應(yīng)變分別為:

2.2 本構(gòu)方程

管道的本構(gòu)關(guān)系采用基于J2塑性流動(dòng)理論的彈塑性模型。

各向異性的屈服函數(shù):

2.3 能量方程

根據(jù)虛功原理，管道的平衡方程可以表達(dá)成如下形式［10］:

外力功增量的變分 為外壓在管壁變形方向作功沿環(huán)向的積分

3 數(shù)值離散

3.1 定解條件

1)邊界條件

在變形過程中，管壁的位移函數(shù)ν、ω需滿足對稱的邊界條件:

2)初始條件

屈曲失穩(wěn)因結(jié)構(gòu)本身或施加載荷的非均勻性而導(dǎo)致。在本模型中，載荷的施加和材料特性具有對稱性和各向同性，因此，需要在初始條件中賦予結(jié)構(gòu)一定的幾何缺陷，即初始橢圓度。當(dāng)初始橢圓度為極小值時(shí)，代表實(shí)際中截面為圓形的管道屈曲情況。

假設(shè)管道具有軸向統(tǒng)一的初始橢圓度Δ0(Δ0=a/R)，此時(shí)中性面上沿環(huán)向各點(diǎn)的徑向初始位移ˉw和環(huán)向初始位移:

由此，在變形過程中任意質(zhì)點(diǎn)處的環(huán)向應(yīng)變表達(dá)式改寫為:

3.2 位移函數(shù)

假設(shè)管道的屈曲為橢圓化屈曲，根據(jù)截面變形的邊界條件，則位移函數(shù)可寫成三角級數(shù)的展開式形式:

位移函數(shù)展開式形式和階數(shù)，在很大程度上影響著分析結(jié)果的精度。

3.3 非線性方程組

將位移離散函數(shù)代入幾何方程、本構(gòu)方程，可得到離散形式的應(yīng)力、應(yīng)變表達(dá)式。對于每個(gè)載荷步，利用虛功原理，將離散形式的能量方程對各個(gè)位移函數(shù)參數(shù)求變分，可得一組非線性方程組，如式(14)所示。每個(gè)載荷步下，此方程組的平衡對應(yīng)著結(jié)構(gòu)的靜力平衡狀態(tài)。

4 求解和分析

4.1 計(jì)算求解

根據(jù)初始缺陷所對應(yīng)的管道結(jié)構(gòu)參數(shù)，進(jìn)行準(zhǔn)靜態(tài)的增量加載，直至管道最大徑向位移隨壓力的增加而開始降低為止?；玖鞒倘鐖D3所示。

4.2 理論分析與有限元分析、壓潰試驗(yàn)及DNV規(guī)范的對比

管道的ABAQUS有限元模型采用C3D8I單元，這樣相比于殼單元可以沿壁厚方向劃分更多的網(wǎng)格，以提高分析精度?？紤]對稱性建立管道模型，在管道兩端施加全固端約束，在圓柱殼的外表面施加均布壓力。根據(jù)理論分析得到壓潰壓力值，設(shè)定一個(gè)高于此值的壓力，運(yùn)用RIKS方法進(jìn)行分析計(jì)算。圖4所示為節(jié)點(diǎn)位移云圖;圖5所示為管道模型的加載壓力與管道被壓縮體積之間的關(guān)系曲線，圖中壓力的峰值即為管道的壓潰壓力。

圖3 分析流程示意Fig.3 Analysis process

圖4 節(jié)點(diǎn)位移云圖Fig.4 Nodal displacement cloud

圖5 加載壓力隨管道壓縮體積變化的曲線過程Fig.5 Curve of the compressed pipe volume with the load

在進(jìn)行理論分析和有限元分析的同時(shí)，進(jìn)行了管道屈曲壓潰試驗(yàn)，試驗(yàn)結(jié)果可與前者進(jìn)行對比驗(yàn)證。試驗(yàn)采用等比例模型，壓力艙體內(nèi)部可容納8 m長的管道試件。試件兩端約束條件為固支，試件外直徑為16英寸(約406 mm)和12英寸(325 mm)兩種。通常，海底管道的屈曲行為與以下因素有關(guān):軸向載荷、環(huán)境溫度、邊界條件、初始缺陷、靜水壓力等。對于管道這種薄壁結(jié)構(gòu)而言，當(dāng)其處于中空狀態(tài)時(shí)最有可能在外部載荷作用下發(fā)生屈曲，這是最危險(xiǎn)的情況。由于實(shí)際的海底管道表面要敷設(shè)保溫材料，可以認(rèn)為這種狀態(tài)下管道內(nèi)部的溫度是恒定的，而且海水環(huán)境溫度也可以視為恒定值。因此，本試驗(yàn)中的工況為管道內(nèi)部中空，端部固支，外部逐漸加載靜水壓，未加載軸力和彎矩，不考慮溫度變化對管道屈曲行為的影響。

圖6為試驗(yàn)管件壓潰前后的形狀對比，可見壓潰后管件變形非常嚴(yán)重，管件變形中間最扁，兩端過度逐漸變形變小。

圖6 試驗(yàn)管件壓潰前后對比Fig.6 Comparison of test tube before and after collapse

通過壓潰試驗(yàn)，深水海底管道屈曲試驗(yàn)測試系統(tǒng)得到了加壓試驗(yàn)的壓力增加過程，從中可以看到壓力從增加到臨界載荷瞬時(shí)降至接近零的全過程，其最高點(diǎn)即為管件的臨界載荷，如圖7所示。

將管道尺寸為325×10、406×10、325×6(mm)的三種管道的壓潰壓力理論結(jié)果與有限元模擬結(jié)果、試驗(yàn)結(jié)果、規(guī)范計(jì)算結(jié)果進(jìn)行比較。管材規(guī)格為API-5LX65，彈性模量 E=2.06×1011Pa，屈服極限 σ0=448 MPa。結(jié)果如表1所示。表中的初始橢圓度值和試驗(yàn)結(jié)果均為對三種尺寸管道進(jìn)行多次管件試驗(yàn)的平均統(tǒng)計(jì)結(jié)果。

圖7 試驗(yàn)中壓力隨時(shí)間加載曲線Fig.7 Curve of test pressure load with time

表1 分析結(jié)果對比Tab.1 Comparison of the results of different motheds

從表1中可以看出，理論計(jì)算方法的壓潰壓力與試驗(yàn)結(jié)果、有限元模擬及規(guī)范值都很接近，證明理論計(jì)算方法的計(jì)算精度的可靠性。

4.3 橢圓度對壓潰壓力的影響

采用理論分析方法，對徑厚比在10～55之間的五種管道在不同初始橢圓度缺陷情況下的壓潰壓力進(jìn)行計(jì)算。計(jì)算分析結(jié)果表明:管道壓潰壓力隨徑厚比的增加和橢圓度缺陷程度的增大而下降，不同的徑厚比受橢圓度影響程度不同，計(jì)算結(jié)果如圖8所示。對于徑厚比為54.17、40.64、32.5、20、10五種情況，橢圓度為0.5%時(shí)，壓潰壓力相對于無初始橢圓度分別下降13.3%、20%、26.7%、16.2%、5.5%。

圖8 不同徑厚比的管道的壓潰壓力-橢圓度曲線關(guān)系Fig.8 Curve of collapse pressure with initial ovility

5 結(jié)語

在考慮軸向均勻一致變形的情況下，建立基于圓環(huán)的管道力學(xué)模型，利用虛功原理建立平衡方程，選用合適的位移離散函數(shù)進(jìn)行方程離散并數(shù)值求解，得到含橢圓度缺陷的管道均勻壓潰的前屈曲階段變形響應(yīng)及壓潰壓力。通過試驗(yàn)、有限元分析結(jié)果、規(guī)范對其進(jìn)行對比驗(yàn)證，計(jì)算了不同橢圓度對壓潰壓力的影響。分析結(jié)果表明:1)該方法的計(jì)算結(jié)果精度可靠，適用于對深水管道的屈曲壓潰強(qiáng)度評估;2)在該方法中，位移離散函數(shù)形式、階數(shù)的選取和管道材料本構(gòu)關(guān)系影響最終結(jié)果的準(zhǔn)確性，對這兩者的選取要盡量精確;3)沿軸向分布一致的橢圓度缺陷對管道的壓潰壓力有明顯影響，減弱程度因不同徑厚比而異，上述計(jì)算曲線為工程應(yīng)用提供參考。

［1］ 周承倜，馬 良.海底管道屈曲及其傳播現(xiàn)象［J］.中國海上油氣，1994，6(6):1-8.

［2］ Fabian O.Collapse of cylindrical elastic tubes under combined bending pressure and axial loads［J］.International Journal of Solids and Structures，1977，13:1257-1270.

［3］ Wilhoit J C，Merwin J E.The effect of axial tension on moment carrying capacity of line pipe stressed beyond the elastic limit［C］∥Proceedings of the Offshore Technology Conference.1971:OTC1335.

［4］ Palmer A C，Martin J H.Buckle propagation in submarine pipelines［J］.Nature(London)，1975，254:46-48.

［5］ Bai Yong，Bai Qiang.Subsea Pipelines and Risers［M］.Oxford:Elsevier Science Ltd，2005.

［6］ Dama E，Karamanos S A.Failure of locally buckled pipelines［J］.ASME，Journal of Pressure Vessel Technology，2007，129:272-279.

［7］ Toscano R G，Mantovano L O，Amenta P M，et al.Collapse arrestors for deepwater pipelines［J］.Cross-over Mechanisms，Computer＆ Structures，2008，86:728-743.

［8］ Xue Jianghong.A non-linear finite-element analysis of buckle propagation in subsea corroded pipelines［J］.Finite Elements in Analysis and Design，2006，42:1211-1219.

［9］ Estefen S F.Collapse behavior of intact and damaged deepwater pipelines and the influence of the reeling method of installation［J］.Journal of Constructional Steel Research，1999，50:99-114.

［10］Kyriakide S，Corona E.Mechanics of Submarine Pipelines［M］.Oxford:Elsevier Science Ltd，2007.