冉龍飛，高文浩，吳 棟，馬天驕

（遼寧石油化工大學(xué)，遼寧 撫順113001）

0 前 言

管道運(yùn)輸是油氣資源配送的主要方式，在國(guó)民經(jīng)濟(jì)中占有重要地位。近年來(lái)，隨著西氣東輸?shù)纫慌L(zhǎng)距離油氣輸送管線的建成和投產(chǎn)，國(guó)內(nèi)管線已全網(wǎng)貫通，并與國(guó)外的油氣資源管線相接。由于管道輸送介質(zhì)易燃易爆等特性，一旦失效，就必然造成災(zāi)難性人員傷亡和經(jīng)濟(jì)損失，并污染管道周?chē)h(huán)境。因此，保證管道安全可靠運(yùn)行對(duì)實(shí)際生產(chǎn)生活具有重要意義。

國(guó)內(nèi)管道失效形式很多，管道懸空是管道失效的主要形式之一。造成管道懸空的主要原因有：①由于地面塌陷、地震等災(zāi)害造成的地殼內(nèi)部構(gòu)造變化導(dǎo)致的管道懸空；②由于滑坡、泥石流等災(zāi)害造成的地殼外部環(huán)境突變導(dǎo)致的管道懸空；③由于凍土、鹽漬土等特殊土體導(dǎo)致的管道懸空。

在長(zhǎng)輸管線運(yùn)行過(guò)程中，受到外界各種載荷的作用，導(dǎo)致管道處于懸空狀態(tài)的案例時(shí)有發(fā)生。如2010年，蘭—成—渝輸油管道德陽(yáng)段管道受洪水沖刷，管道大面積懸空，管道多處被拉斷，導(dǎo)致部分管線停輸，造成嚴(yán)重?fù)p失。通過(guò)力學(xué)計(jì)算，其允許的極限懸空長(zhǎng)度為44.8 m，但實(shí)際懸空380 m管道仍未發(fā)生斷裂。說(shuō)明對(duì)處于懸空工況下的管道應(yīng)力分析已不能滿足實(shí)際的分析要求。

本研究通過(guò)理論計(jì)算得出懸空管道的極限應(yīng)變和容許應(yīng)變，并應(yīng)用ABAQUS有限元分析軟件模擬在管道不同懸空長(zhǎng)度下管道的應(yīng)變?nèi)≈?。?dāng)應(yīng)變?nèi)≈到咏碚撚?jì)算得到的容許應(yīng)變值時(shí)，對(duì)應(yīng)的管道懸空長(zhǎng)度為該管材管道的極限懸空長(zhǎng)度，以此作為管道工程評(píng)定的標(biāo)準(zhǔn)，保證管道的安全運(yùn)行[1-2]。

1 埋地懸空管道應(yīng)力和應(yīng)變的理論計(jì)算

1.1 力學(xué)模型的選擇

在實(shí)際工程模擬計(jì)算中，懸空管道力學(xué)模型可簡(jiǎn)化為以下三種類(lèi)型：

（1）簡(jiǎn)支梁模型。該模型僅考慮懸空段管道的受力情況，不考慮管－土之間相互作用對(duì)未懸空段管道的影響。

（2）彈性地基梁模型。該模型將管道未懸空段看做半無(wú)限長(zhǎng)符合Winkler假設(shè)的彈性地基梁，即認(rèn)為兩端未懸空段的土壤物性、管道及管道變形關(guān)于中心對(duì)稱(chēng)。

（3）彈塑性地基模型。該模型認(rèn)為管道處于懸空工況時(shí)，受自重載荷、輸送介質(zhì)重力載荷等多種均布載荷的共同作用，使管道產(chǎn)生撓曲變形，而且埋地段管道由于受到管道軸向收拉伸力及土壤阻力作用，是管道產(chǎn)生彎曲變形[3]。

上述三種方法，簡(jiǎn)支梁模型未考慮懸空段與未懸空段之間的相互影響，全彈性地基梁模型僅考慮了管道的彈性變形，未考慮管道的塑性變形，均與實(shí)際情況不符。因此采用彈性地基梁模型，既考慮了懸空管道的彈性變形和塑性變形，又對(duì)埋地懸空管道的實(shí)際受力情況作了簡(jiǎn)化。

1.2 應(yīng)力計(jì)算

埋地懸空管道的受力簡(jiǎn)化模型如圖1所示。由圖可見(jiàn)，埋地管道懸空段受到均布荷載q作用。均布載荷包括管道自重、輸送介質(zhì)質(zhì)量和懸空管道埋入端上方土對(duì)管道的壓力等。管道兩端受到當(dāng)量軸向力S0的作用，當(dāng)埋地管道懸空長(zhǎng)度為L(zhǎng)時(shí)，管道埋入端受到彎矩M0的作用。在實(shí)際分析計(jì)算時(shí)，將懸空管道兩端非懸空段看作半無(wú)限長(zhǎng)符合Winkler假設(shè)的彈性地基梁，近似認(rèn)為土壤物性、埋地管道的受力和變形關(guān)于c-c軸對(duì)稱(chēng)[4]。

圖1 埋地懸空管道的受力簡(jiǎn)化模型

懸空管道的撓度微分方程為

式中：E—彈性模量；

I—慣性矩；

M0—管道在x=0處的彎矩；

v0—管道在x=0處的撓度；

q—管道所受載荷；

L—管道懸空長(zhǎng)度；

S0—懸空管道的當(dāng)量軸向力，S0=N0-Np，

N0—管道軸力，拉力為正值，壓力為負(fù)值；

p—管道承受的內(nèi)壓；

d—管道內(nèi)徑。

根據(jù)埋地懸空管道土壤物性、受力的對(duì)稱(chēng)性以及邊界條件，可得管道在x=L/2處的撓度為

1.3 應(yīng)變計(jì)算

基于應(yīng)力的管道強(qiáng)度設(shè)計(jì)方法一般適用于以應(yīng)力為控制參量的載荷作用下的彈性設(shè)計(jì)。而基于應(yīng)變的設(shè)計(jì)方法一般適用于在發(fā)生地震、塌陷、洪水、滑坡等重大地質(zhì)災(zāi)害時(shí)，管道會(huì)發(fā)生變形，此時(shí)作用在管道上的載荷是以位移為控制參量來(lái)衡量的。采用基于應(yīng)力的設(shè)計(jì)方法難以使管道滿足強(qiáng)度設(shè)計(jì)的要求，因此，將基于應(yīng)變?cè)O(shè)計(jì)方法的作為基于應(yīng)力設(shè)計(jì)的補(bǔ)充。兩種設(shè)計(jì)方法的工作區(qū)域和控制參量不同，導(dǎo)致計(jì)算過(guò)程也不相同。

基于應(yīng)變的設(shè)計(jì)方法主要是確定管道處于懸空狀態(tài)時(shí)將要承受的容許應(yīng)變和管道本身的極限應(yīng)變。為防止在懸空狀態(tài)下管道過(guò)量變形使管道出現(xiàn)裂紋，甚至拉斷，導(dǎo)致管道失效，在實(shí)際理論計(jì)算過(guò)程中，存在拉伸應(yīng)變、壓縮應(yīng)變和橢圓化變形三種類(lèi)型的應(yīng)變[5]。

1.3.1 拉伸應(yīng)變

管道處于懸空工況時(shí)，可能存在缺陷，拉伸應(yīng)變?nèi)≈颠_(dá)到或超過(guò)極限應(yīng)變?nèi)≈禃r(shí)，可能導(dǎo)致管道斷裂，拉伸應(yīng)變?nèi)≈祽?yīng)滿足（3）式要求

式中：εtf—縱向或周向的因子化拉伸應(yīng)變；

φεt—拉伸應(yīng)變阻力因子；

εtcrit—管壁或焊接部位的極限拉伸應(yīng)變。

拉伸極限應(yīng)變依據(jù)CSA Z662—2007[6]計(jì)算得到，容許拉伸應(yīng)變的阻力因子取值為0.7。

1.3.2 壓縮應(yīng)變

管道懸空時(shí)未懸空段管道承受壓縮載荷或者彎矩，管道會(huì)產(chǎn)生壓應(yīng)力和應(yīng)變，當(dāng)管壁最大壓應(yīng)變達(dá)到或者超過(guò)臨界水平，管壁會(huì)出現(xiàn)局部屈曲或褶皺。因此，壓縮應(yīng)變?nèi)≈祽?yīng)滿足（4）式要求

式中：εcf—縱向或周向的因子化壓縮應(yīng)變；

φεc—壓縮應(yīng)變阻力因子；

εccrit—軸向或周向壓縮極限應(yīng)變。

依據(jù)CSA Z662—2007，容許壓縮應(yīng)變的阻力因子取值為0.8。

當(dāng)考慮內(nèi)壓時(shí)，壓縮應(yīng)變應(yīng)按照（5）式計(jì)算

式中：t—管道壁厚；

D—管道外徑；

pi—最大設(shè)計(jì)內(nèi)壓；

1.3.3 橢圓化變形

依據(jù)美國(guó)船級(jí)社（ABS）設(shè)計(jì)規(guī)范，為防止管道嚴(yán)重橢圓化導(dǎo)致局部屈曲，應(yīng)變極限滿足（7）式要求

pe—最小外部靜水壓力；

Es—管道的彈性模量；

Fy—最小屈服強(qiáng)度。

當(dāng)不考慮內(nèi)壓時(shí)，依據(jù)薄殼褶皺理論，許可壓縮應(yīng)變?yōu)?/p>

式中：ε—彎曲應(yīng)變；

εb—管道純彎曲時(shí)的屈曲應(yīng)變；

f0—外圓率，運(yùn)輸過(guò)程中要

求f0不超過(guò)0.75%。

2 有限元分析

埋地管道隨土層運(yùn)動(dòng)而發(fā)生變形，基于埋地懸空管道的受力特點(diǎn)，應(yīng)用ABAQUS有限元分析軟件進(jìn)行模擬計(jì)算。ABAQUS是一套功能齊全的工程模擬有限元分析軟件，可以解決模擬件結(jié)構(gòu)的各種問(wèn)題，范圍包括從簡(jiǎn)單的線性分析問(wèn)題到復(fù)雜的非線性分析問(wèn)題。ABAQUS針對(duì)埋地管道開(kāi)發(fā)了專(zhuān)門(mén)模擬管道與周?chē)寥乐g相互作用的虛擬單元，即管－土相互作用單元（PSI）。ABAQUS提供的管－土相互作用單元包括：二維單元PSI24和PSI26，三維單元PSI34和PSI36。埋地管道的模擬可以用梁、管或彎接頭單元來(lái)進(jìn)行，地基受力和管道與周?chē)寥乐g的相互作用通過(guò)PSI單元來(lái)模擬。管－土相互作用單元如圖2所示，PSI單元與下面模擬管道的梁、管或彎接頭單元共用節(jié)點(diǎn)，地面的實(shí)際運(yùn)動(dòng)情況在模擬過(guò)程中通過(guò)施加邊界條件來(lái)加以描述[7-8]。

圖2 管－土相互作用單元簡(jiǎn)圖

在不考慮風(fēng)載荷、溫差、外界震動(dòng)、初始裝配應(yīng)力等其他因素影響，對(duì)懸空管道受力情況進(jìn)行模擬分析，確定采取何種單元來(lái)模擬懸空管道的實(shí)際受力情況。ABAQUS可供選擇的模擬單元有梁、管或彎接頭單元。本研究采用梁、管單元的實(shí)際受力情況。ABAQUS可供選擇的模擬單元有梁、管或彎接頭單元。本研究采用梁、管單元來(lái)模擬埋地懸空管道的受力情況，根據(jù)懸空管道的受力特點(diǎn)采用彈塑性地基模型，對(duì)懸空管道的實(shí)際工況進(jìn)行簡(jiǎn)化來(lái)模擬。在計(jì)算過(guò)程中管－土之間的相互作用通過(guò)在管單元節(jié)點(diǎn)上的軸向土彈簧、豎向土彈簧和水平橫向土彈簧來(lái)近似代替。土彈簧模型如圖3所示，該模型將管體周?chē)翆雍?jiǎn)化為等效彈塑性彈簧，通過(guò)彈簧剛度和自由度來(lái)定義管－土之間的作用。土彈簧剛度依據(jù) 《油氣輸送管道線路工程抗震技術(shù)規(guī)范》（GB 50470—2008）附錄E要求計(jì)算得到[9]。管－土相互作用屬于非線性問(wèn)題，應(yīng)用ABAQUS有限元分析軟件能模擬非線性計(jì)算非線性問(wèn)題。

圖3 土彈簧模型簡(jiǎn)圖

3 算 例

懸空埋地管道管材為X70，管道外徑1219mm，壁厚為18.4 mm，管材密度為7 800 kg/m3，彈性模量為206 GPa，泊松比為0.3，設(shè)計(jì)系數(shù)為0.72，管道的最小屈服強(qiáng)度為485 MPa，管道的工作壓力為10 MPa。管道安裝時(shí)的溫度為27℃，發(fā)生懸空時(shí)的溫度為23℃，管道埋深為1.5 m，土體密度為1 700 kg/m3，土的內(nèi)摩擦角為30°，地基系數(shù)k0=4×104kN/m。計(jì)算得到管道自重為55.1 kN/m，管道上方土壤的重力為48.9 kN/m，單位長(zhǎng)度橫向土彈簧剛度為7.15×105N/m2，垂直向上土彈簧剛度為2.68×105N/m2。在不考慮軸向載荷的情況下，應(yīng)用ABAQUS軟件計(jì)算得到在不同懸空長(zhǎng)度下管道的應(yīng)力、應(yīng)變和撓度取值見(jiàn)表1[10]。

表1 不同懸空長(zhǎng)度下管道ABAQUS有限元分析的數(shù)值計(jì)算結(jié)果

圖4 不同懸空長(zhǎng)度下管道的應(yīng)力變化曲線

應(yīng)用ABAQUS有限元分析軟件得到該工況下管道的數(shù)值解，可見(jiàn)隨著管道懸空長(zhǎng)度的增加，最大應(yīng)力，最大位移以及應(yīng)變值都增大。該管材管道在不同懸空長(zhǎng)度下管道的應(yīng)力變化曲線如圖4所示。

通過(guò)理論計(jì)算得到管道的極限壓縮應(yīng)變?nèi)≈禐?.77%，容許壓縮應(yīng)變?nèi)≈禐?.62%。軟件模擬得到管道的應(yīng)變?nèi)≈到咏艿赖娜菰S壓縮應(yīng)變時(shí)對(duì)應(yīng)的懸空長(zhǎng)度取值為350 m。管道懸空長(zhǎng)度超過(guò)該值時(shí)，管道可能被拉斷，存在安全隱患，應(yīng)采取防護(hù)和監(jiān)測(cè)措施。

4 結(jié) 語(yǔ)

對(duì)于處于懸空工況下的埋地管道，應(yīng)力分析不能滿足實(shí)際要求，采用基于應(yīng)變的分析方法能對(duì)處于變形的管道作出安全評(píng)定。對(duì)于算例中的管材而言，當(dāng)管道的極限懸空長(zhǎng)度達(dá)到350 m時(shí)，管道的應(yīng)變值為0.61%，接近管道的最大允許應(yīng)變，一旦超過(guò)該應(yīng)變值，該管材的管道可能被拉斷，存在安全隱患。因此在實(shí)際工程中，可將極限懸空長(zhǎng)度作為管線運(yùn)行的安全評(píng)定準(zhǔn)則。

［1］由小川，莊茁.高壓天然氣管線在地質(zhì)災(zāi)害下的失效分析［J］.天然氣工業(yè)，1999，19（04）：77-81.

［2］王滬毅.輸氣管線在地質(zhì)災(zāi)害中的力學(xué)行為研究［D］.西安：西北工業(yè)大學(xué)，2003.

［3］馬廷霞，吳錦強(qiáng)，唐愚，等.成品油管道的極限懸空長(zhǎng)度研究［J］.西南石油大學(xué)學(xué)報(bào)（自然科學(xué)版），2012 （04）：165-173.

［4］羅金恒，趙新偉，王峰會(huì)，等.地質(zhì)災(zāi)害下懸空管道的應(yīng)力分析及計(jì)算［J］.壓力容器，2006，23（04）：23-26.

［5］王峰會(huì)，趙新偉.高壓鋼管黃土塌陷情況下的力學(xué)分析與計(jì)算［J］.油氣儲(chǔ)運(yùn)，2004，23 （04）：6-8.

［6］ CSA Z662—2007，Oil and Gas Pipeline System［S］.

［7］符圣聰.用于埋管抗震計(jì)算的地面位移和地基反力系數(shù)［R］.北京：中國(guó)建筑科學(xué)研究院，1996.

［8］邢靜忠，柳春圖.線彈性土壤中埋設(shè)懸跨管道的屈曲分析［J］.工程力學(xué)，2008，25（10）：72-75.

［9］王小龍，姚安林.埋地鋼管局部懸空的撓度和內(nèi)力分析［J］.工程力學(xué)，2008，25（08）：218-222.

［10］唐永進(jìn).壓力管道應(yīng)力分析［M］.北京：中國(guó)石化出版社，2003：96-97.