內(nèi)壓作用下含體積型缺陷彎管極限載荷研究*

王佳音，帥 健，劉道乾，孫 偉，許 葵，呂志陽

(1.中國石油大學(xué)(北京) 安全與海洋工程學(xué)院，北京 102249；2.中國石化銷售股份有限公司華中分公司，湖北 武漢 430000)

0 引言

油氣管道中彎管特殊的幾何特征導(dǎo)致其受力不同，一旦在彎頭處產(chǎn)生腐蝕或沖蝕等缺陷，在承受內(nèi)壓或其他外載荷作用時(shí)，彎管處更容易出現(xiàn)應(yīng)力集中，使其成為管路中容易失效的部分[1-2]。國內(nèi)外因油氣彎管失效引起的事故眾多，因此，研究彎管處腐蝕缺陷對(duì)管道失效的影響具有現(xiàn)實(shí)意義。

目前，許多學(xué)者在該領(lǐng)域開展了相關(guān)研究：陳鋼等[3]利用數(shù)值分析研究了內(nèi)壓和彎矩作用下局部減薄彎頭極限載荷的變化規(guī)律并討論了彎頭的典型失效模式；Oh等[4]在假定局部減薄形狀為矩形的情況下定量研究了面內(nèi)彎矩作用時(shí)彎頭的塑性行為；楊秀娟等[5]分析了單一內(nèi)壓作用下含局部減薄缺陷彎管的極限載荷變化情況，得出缺陷尺寸、位置和彎曲半徑等對(duì)極限載荷的影響規(guī)律；Kim等[6]在單一內(nèi)壓作用下，對(duì)含局部減薄彎頭進(jìn)行了一系列爆破試驗(yàn)，分析了不同減薄位置和缺陷幾何形狀(如減薄深度、長度和圓周角)等對(duì)彎頭失效壓力的影響；Lee等[7]使用ANSYS非線性有限元分析軟件，預(yù)測(cè)海底管道腐蝕缺陷彎頭的爆破壓力，發(fā)現(xiàn)彎頭的爆破壓力不僅與缺陷尺寸相關(guān)，還與缺陷位置的變化相關(guān)。

綜上所述，已有文獻(xiàn)主要集中在對(duì)局部減薄缺陷極限載荷的研究上，忽略了體積型缺陷處的應(yīng)力集中，本文通過有限元分析主要研究內(nèi)壓作用下點(diǎn)蝕及溝槽型腐蝕缺陷彎管的塑性極限載荷。

1 有限元模型

1.1 有限元模型構(gòu)建

選用高階三維20節(jié)點(diǎn)的固體結(jié)構(gòu)單元，構(gòu)建含體積型缺陷90°彎管的有限元模型，該單元具有塑性、超彈性和大變形等能力，應(yīng)用該單元能夠較為準(zhǔn)確地得到含缺陷彎管的應(yīng)力分布[8]。有限元網(wǎng)格劃分如圖1所示。為清楚地描述缺陷位置，采用環(huán)形坐標(biāo)系，如圖2所示，其中，θ為缺陷時(shí)鐘角度，(°)；ψ為軸向角度，(°)；r為徑向，mm。

圖1 有限元模型圖及網(wǎng)格劃分Fig.1 Finite element model and mesh generation

圖2 環(huán)形坐標(biāo)系Fig.2 Annular coordinate system

體積型缺陷類型如圖3所示。缺陷長度和寬度較小時(shí)，將缺陷形狀簡化為球形，如圖3(a)所示；缺陷為細(xì)長型時(shí)，為避免缺陷端部出現(xiàn)應(yīng)力集中，端部采用橢球形、槽身采用柱面進(jìn)行建模，如圖3(b)所示。

1.2 材料性能及管道規(guī)格

管材選用X80鋼材料，其目前在中國得到廣泛應(yīng)用[9-10]。X80鋼的基本力學(xué)性能見表1。為能夠真實(shí)反映彎管的受力變形和極限承載能力，同時(shí)考慮材料的幾何非線性和材料非線性。根據(jù)試驗(yàn)結(jié)果，選擇屈服應(yīng)力σs=641 MPa，抗拉強(qiáng)度σb=740 MPa。材料模型選用Ramberg-Osgood模型[11]，表達(dá)式如式(1)所示：

(1)

表1 X80鋼基本力學(xué)性能Table 1 Basic mechanical properties of X80 steel

1.3 約束和載荷

只考慮彎管受到內(nèi)壓作用，在管道內(nèi)表面施加均勻分布的內(nèi)壓載荷P0。對(duì)彎管一端進(jìn)行3個(gè)平動(dòng)自由度的固定，對(duì)另一端施加等效軸向應(yīng)力以平衡內(nèi)壓，施加的等效軸向應(yīng)力σa為：

(2)

式中：D和Di分別為管子的外徑和內(nèi)徑，mm；σa為等效軸向應(yīng)力，MPa；P0為內(nèi)壓載荷，MPa。

1.4 失效準(zhǔn)則的確定

把抗拉強(qiáng)度σb作為管道失效的判定標(biāo)準(zhǔn)，認(rèn)為當(dāng)缺陷底部中心點(diǎn)的Von Mises等效應(yīng)力達(dá)到抗拉強(qiáng)度時(shí)，判定管線失效。為進(jìn)一步驗(yàn)證該準(zhǔn)則的適用性，建立無缺陷彎管的有限元模型，復(fù)現(xiàn)文獻(xiàn)[12]中爆破試驗(yàn)的過程，并將本文模擬結(jié)果與試驗(yàn)數(shù)據(jù)進(jìn)行對(duì)比，有限元計(jì)算結(jié)果Von Mises等效應(yīng)力云圖如圖4。

圖4 Von Mises等效應(yīng)力(MPa)Fig.4 Von Mises equivalent stress(MPa)

爆破試驗(yàn)彎管為φ406.4 mm×9.5 mm，材料為X60，彎曲半徑為5D，試驗(yàn)測(cè)得爆破壓力值為29.1 MPa。利用有限元計(jì)算分析，等效應(yīng)力達(dá)到抗拉強(qiáng)度σb=595 MPa(X60的抗拉強(qiáng)度極限)時(shí)，得到的極限載荷為29.236 MPa，誤差為0.4%，結(jié)果基本一致。經(jīng)過與文獻(xiàn)中爆破試驗(yàn)數(shù)據(jù)對(duì)比，驗(yàn)證了選擇抗拉強(qiáng)度σb作為判定彎管失效的失效準(zhǔn)則是合理的。

1.5 驗(yàn)證模型的準(zhǔn)確性

為驗(yàn)證模型加載方式及邊界條件等施加是否準(zhǔn)確，選取文獻(xiàn)[13]中含體積型缺陷的算例進(jìn)行有限元結(jié)果復(fù)現(xiàn)。文獻(xiàn)[13]中模型為φ1 219 mm×19.89 mm，材料為X80，彎曲半徑為6D的彎管，缺陷長、寬、深分別為：165，50，9.945 mm。文獻(xiàn)中管道失效判定標(biāo)準(zhǔn)選用2倍彈性斜率法，故計(jì)算時(shí)亦選用2倍彈性斜率法，求解極限載荷值如圖5所示。

圖5 極限載荷求解Fig.5 Solution of ultimate load

本文計(jì)算結(jié)果為14.100 MPa，與文獻(xiàn)中的結(jié)果13.930 MPa基本一致，誤差為1.2%。因此，可以驗(yàn)證本模型載荷及邊界條件施加方式合理。

2 含體積型缺陷彎管極限載荷影響因素分析

本文主要研究缺陷幾何尺寸(缺陷的長度和深度)，缺陷位置(時(shí)鐘位置)和彎曲半徑對(duì)彎管承載能力的影響。

令Θ為多邊形,頂點(diǎn)表示為Q1,Q2,…,Qn。此外,引用|·|表示歐式距離,例如|AB|表示點(diǎn)A與點(diǎn)B間的歐式距離,||表示線的歐式長度。pΘ表示Θ的邊。令A(yù)、B是Θ邊上的點(diǎn),而分別表示從A至B的逆時(shí)針、順時(shí)針邊。

2.1 缺陷長度、深度對(duì)塑性極限載荷的影響

Choi[14]通過大量模擬實(shí)驗(yàn)得出缺陷寬度對(duì)極限載荷的影響可以忽略不計(jì)，故本文不考慮缺陷寬度的影響，只考慮缺陷長度和深度的影響。

計(jì)算模型幾何尺寸為φ1 219 mm×19.89 mm，保持彎曲半徑為4D和缺陷寬度不變，只改變?nèi)毕莸拈L度和深度，分別研究缺陷位于內(nèi)拱θ=0°，幾何中心線θ=90°，外拱θ=180°時(shí)極限載荷的變化趨勢(shì)，如圖6～11所示(d/t為相對(duì)缺陷深度，即缺陷深度與壁厚的比值)。缺陷長度選取范圍為50～2 000 mm；缺陷深度選取壁厚的20%，30%，40%，50%，60%，70% 6個(gè)等級(jí)。

圖6 內(nèi)拱處極限載荷隨缺陷相對(duì)深度的變化Fig.6 Variation of ultimate load with relative depth of defect at intrados

圖7 內(nèi)拱處極限載荷隨缺陷長度的變化Fig.7 Variation of ultimate load with length of defect at intrados

圖8 幾何中心線處極限載荷隨缺陷相對(duì)深度的變化Fig.8 Variation of ultimate load with relative depth of defect at geometry centerline

圖9 幾何中心線處極限載荷隨缺陷長度的變化Fig.9 Variation of ultimate load with length of defect at geometry centerline

圖10 外拱處極限載荷隨缺陷相對(duì)深度的變化Fig.10 Variation of ultimate load with relative depth of defect at extrados

圖11 外拱處極限載荷隨缺陷長度的變化Fig.11 The variation of failure pressure with the lengthsat extrados

由計(jì)算結(jié)果可知：

1)圖6，8，10表明不論是在內(nèi)拱、外拱還是幾何中心線處，彎管極限載荷隨缺陷深度增加而不斷減小，隨著缺陷長度增大，極限載荷隨深度增加而減小的速度會(huì)逐漸加快。

2)圖7，9，11表明不論在內(nèi)拱、外拱還是中心線處，彎管極限載荷隨缺陷長度增加而減小。當(dāng)缺陷深度不太深時(shí)，極限載荷減小的速度緩慢，隨著缺陷深度變深，極限載荷減小的速度加快。同一缺陷深度，隨著缺陷長度增加到一定程度，彎管極限載荷趨于一個(gè)定值，這個(gè)臨界缺陷長度約等于管子外徑。

2.2 缺陷位置對(duì)于塑性極限載荷的影響

為研究各種時(shí)鐘角度對(duì)極限載荷的影響規(guī)律，改變?nèi)毕莸臅r(shí)鐘角度θ，分別計(jì)算了θ為0°，30°，45°，60°，90°，120°，135°，180°時(shí)的極限載荷。計(jì)算模型為φ1 219 mm×22 mm，R=4D，缺陷長度為500，1 000 mm，缺陷深度均為50%t，70%t，結(jié)果如圖12所示。

圖12 不同位置極限載荷Fig.12 Ultimate loads at different positions

由圖12，13可以看出：

1)缺陷位于內(nèi)拱，即θ=0°時(shí)，彎管的極限載荷最小；缺陷位于外拱，即θ=180°時(shí)，極限載荷最大；缺陷位于幾何中心線時(shí)，極限載荷值處于兩者之間。

2.3 彎曲半徑對(duì)于塑性極限載荷的影響

保持缺陷尺寸相同，改變彎管彎曲半徑，研究彎管極限載荷變化規(guī)律，分析計(jì)算結(jié)果如圖13所示，其中相對(duì)彎曲半徑為R/D。由圖13可知，隨彎曲半徑增大，彎管內(nèi)拱處極限載荷逐漸增加，外拱處極限載荷逐漸減小，中心線處極限載荷不隨彎曲半徑變化而改變，內(nèi)拱和外拱的極限載荷最終趨于定值。

圖13 不同彎曲半徑的極限載荷Fig.13 Ultimate loads under different bending radiuses

3 極限內(nèi)壓預(yù)測(cè)公式

3.1 擬合公式

文獻(xiàn)[15]中提出含體積型缺陷直管的極限內(nèi)壓預(yù)測(cè)公式：

(3)

(4)

式中：A，B，C為參數(shù)；Plimit為含體積型缺陷直管的極限內(nèi)壓，MPa；t為管道壁厚，mm；D為管道外徑，mm；σb為抗拉強(qiáng)度，MPa；d為缺陷深度，mm；L為缺陷長度，mm。根據(jù)影響因素的分析，并對(duì)不同管徑、壁厚等多組管材進(jìn)行有限元計(jì)算，重新擬合參數(shù)并附加1個(gè)新的函數(shù)作為彎管系數(shù)，如式(5)所示。

(5)

根據(jù)討論出的含缺陷彎管各因素對(duì)管道失效影響的規(guī)律，提出彎管系數(shù)的函數(shù)形式為：

(6)

式中：a，b，c為參數(shù)；f2為彎管系數(shù)；θ為缺陷時(shí)鐘角度，(°)；R為彎曲半徑，mm。根據(jù)大量模擬計(jì)算，應(yīng)用Matlab擬合得到：

(7)

最終得到的含體積型缺陷彎管極限內(nèi)壓預(yù)測(cè)公式(8)～(11)。

(8)

(9)

(10)

(11)

式中：Pelbow為含體積型缺陷彎管的極限內(nèi)壓，MPa；P1為無缺陷直管極限載荷，MPa。

3.2 公式驗(yàn)證

將有限元模擬結(jié)果與公式計(jì)算結(jié)果進(jìn)行對(duì)比，并作誤差分析，如表2所示。由表2可見，預(yù)測(cè)公式結(jié)果與有限元結(jié)果基本一致，最大誤差不超過5%，公式基本可以滿足應(yīng)用的需要。

4 結(jié)論

1)彎管極限載荷隨缺陷長度和深度增加而減小。且同一缺陷深度下，缺陷長度達(dá)到一個(gè)臨界值時(shí)，極限載荷值將不再下降，這個(gè)臨界值約等于彎管直徑。

2)缺陷的時(shí)鐘位置不同，極限載荷值不同。缺陷位于外拱時(shí)極限載荷值最大，位于內(nèi)拱時(shí)極限載荷值最小，位于幾何中心線時(shí)的極限載荷值處于二者之間；且時(shí)鐘角度與極限載荷之間的變化規(guī)律呈現(xiàn)余弦函數(shù)的關(guān)系。

表2 極限內(nèi)壓公式誤差分析Table 2 Error analysis on formula of ultimate internal pressure

3)彎管極限載荷隨彎曲半徑增大，內(nèi)拱處極限載荷增大，外拱處極限載荷降低。

4)提出含缺陷彎管的極限內(nèi)壓預(yù)測(cè)公式，該公式可基本滿足工程應(yīng)用的需要。