冉龍飛，高文浩，吳 棟

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

0 前 言

由于長輸油氣管道服役環(huán)境復(fù)雜，管道敷設(shè)不可避免地會經(jīng)過各種地質(zhì)災(zāi)害頻發(fā)區(qū)，管道運行存在潛在風(fēng)險。例如1994年，長慶油田元-悅段輸油管道受洪水沖刷，導(dǎo)致管道大面積懸空，另外上游的泥石流將管道沖毀26處，管道多處被拉斷，導(dǎo)致管道停輸，造成嚴(yán)重?fù)p失，說明埋地管道懸空是造成管道失效的主要原因之一。因為在地質(zhì)災(zāi)害頻發(fā)區(qū)可能存在大跨度的埋地管道懸空，在管道自重、輸送介質(zhì)質(zhì)量等載荷共同作用下，管道出現(xiàn)下垂，較大的下垂高度導(dǎo)致管道變形，甚至拉斷管道，造成嚴(yán)重的事故[1]。因此，研究埋地懸空管道的應(yīng)力分布，對管道的安全運行具有實際意義。

1 懸空管道的受力模型

埋地懸空管道兩端受到土體的彈性支撐，可將管道看作兩端受彈性支撐的梁結(jié)構(gòu)，因此懸空管道的受力模型可簡化為變形梁模型[2-3]，如圖1所示。

圖1 懸空管道力學(xué)模型

懸空段管道受到管道自重、輸送介質(zhì)質(zhì)量、懸空管道埋入端上方土對管道的壓力等均布載荷q的共同作用，并受到集中載荷作用。基于Winkler線性理論將懸空管段兩端非懸空段管道看作半無限長的彈性地基梁，近似認(rèn)為土壤物性、埋地管道的受力和變形關(guān)于軸對稱。

2 載荷計算

埋地管道在運行過程中受到各種載荷作用，靜載荷主要包括管道自重、管道輸送介質(zhì)質(zhì)量、防腐材料等保護(hù)層重、內(nèi)壓和溫度變化對管道的作用、管道上覆土層對管道的作用等，動載荷主要包括地層塌陷、沖溝、滑坡、地震等地質(zhì)災(zāi)害所產(chǎn)生的載荷[4]。

2.1 重力荷載

（1）管道自重載荷

式中：Gρ—管子的質(zhì)量；

ρ—管材密度；

D—管材外徑；

d—管材內(nèi)徑；

L—管道長度。

（2）輸送介質(zhì)載荷

式中：Gj—管子所含介質(zhì)的質(zhì)量；

ρi—輸送介質(zhì)的密度。

（3）防腐材料載荷

式中：Gb—管子上防腐材料質(zhì)量；

ρb—防腐材料密度；

dc—防腐層的平均直徑；

tb—防腐層平均厚度。

（4） 土壤載荷

2.2 應(yīng) 力

實際運行狀態(tài)下，懸空管道受到內(nèi)外壓力和溫度變化的共同作用。所受應(yīng)力中對懸空管道強度和穩(wěn)定性造成影響的應(yīng)力分別為軸向應(yīng)力、環(huán)向應(yīng)力和剪應(yīng)力[5]。

（1）軸向應(yīng)力。地下管道產(chǎn)生軸向應(yīng)力的原因是溫度的變化和環(huán)向應(yīng)力的泊松效應(yīng)，依照（4）式計算：

式中：σ—管道的軸向應(yīng)力；

σt—管道溫差產(chǎn)生的應(yīng)力；

σa—管道環(huán)向應(yīng)力；

E—管材彈性模量；

ΔT—管道溫度變化量；

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

忘記是從哪天開始，她每天在網(wǎng)上沖浪的時候，他都會在。他們會談很多心事，包括小時候的頑劣，但是，卻避免著現(xiàn)實。

α—管道材料的線膨脹系數(shù)；

μ—泊松比；

t—管道壁厚。

（2）剪應(yīng)力。局部懸空的埋地管道剪力很大，不能忽略。在埋地管道圓形橫截面上下邊緣處為零，在橫截面中性軸處絕對值最大。局部懸空的圓形埋地管道橫截面中性軸處的剪應(yīng)力為：

式中：τ—剪應(yīng)力；

Fz—管道豎向的總壓力。

（3）環(huán)向應(yīng)力。環(huán)向應(yīng)力由管道輸送介質(zhì)內(nèi)壓產(chǎn)生，內(nèi)壓是確定管道應(yīng)力狀態(tài)的主要作用力之一，內(nèi)壓產(chǎn)生的環(huán)向應(yīng)力σh為：

3 管-土相互作用關(guān)系

管-土之間相互作用，將土對管道的支撐反力用彈簧近似代替，根據(jù)場地土的性質(zhì)和管道參數(shù)，采用三個方向上土彈簧表示[6-7]。

軸向彈簧剛度

式中：k0—土靜壓力系數(shù)；

D—管道外徑；

Cs—回填土內(nèi)聚系數(shù)；

k—粘合系數(shù)；

δ—管土截面摩擦角。

橫向土彈簧剛度

式中：C0～C4為計算Nqh的系數(shù)，按表1要求取值。

表1 計算Nqh時C0～C4的取值

豎向土彈簧剛度

式中：準(zhǔn)—內(nèi)磨擦角。

相鄰?fù)翉椈芍g距離，經(jīng)過反復(fù)驗證，間距滿足式（15）和式（16）要求時，誤差可以控制在較小的范圍內(nèi)[5]。

式中：β—特征系數(shù)；

E—管材彈性模量；

I—慣性矩。

4 有限元計算

4.1 計算模型

管-土相互作用屬于復(fù)雜的非線性問題，模擬此類問題采用ABAQUS有限元分析軟件。該軟件對復(fù)雜的非線性問題能夠進(jìn)行較為準(zhǔn)確的計算，而且該軟件開發(fā)了專門針對埋地管道的模擬單元，即管-土相互作用單元（PSI）。對懸空管道受力進(jìn)行分析，計算中不考慮風(fēng)載、溫差、外界震動、初始裝配應(yīng)力等其他影響因素。ABAQUS可供選擇的模擬單元有梁、管或彎接頭單元。埋地管道與土壤之間的相互作用通過管-土相互作用單元來模擬[8]。本研究采用梁、管單元來模擬埋地懸空管道的受力特點。工程中采用簡化模型來分析管土間的相互作用關(guān)系。目前常用的管土模型主要包括彈性地基梁、土彈簧和非線性接觸三種。彈性地基梁模型基于靜力分析，假設(shè)管線周圍土體均勻分布，將管線視為彈性梁模型，主要考慮土體最終位移對管線的作用。彈性地基梁模型理論簡單明了，易于手算，因而被工程界廣泛接受，并被美國輸油（氣）抗震規(guī)范所采用。懸空管道兩端埋入土壤中，管道向下變形過程中受到軸向力的作用，并且承受彎矩，因此采用彈性地基梁模型，模擬實際懸空管道工況。埋地管道土層運動情況通過管-土相互作用單元（PSI）來模擬，如圖2所示[9-10]。

圖2 管-土相互作用單元簡圖

該單元只有位移自由度，該單元的一面或者一邊與下面的模擬管道的梁、管或彎接頭單元共用節(jié)點，而另一側(cè)的節(jié)點代表代表地面，并用于施加邊界條件來描述遠(yuǎn)場地面的運動。懸空管道埋入端土對管道的地基抗力用具有一定剛度的彈簧近似代替。

4.2 參數(shù)選取及載荷計算

參數(shù)：管材鋼級X70，管徑1 219 mm，壁厚18.4 mm的埋地懸空管道，設(shè)計壓力12 MPa，彈性模量210 GPa，泊松比0.3，慣性矩7.8×10-4m4。管材密度7 800 kg/m3，最小屈服強度483 MPa，埋深1.5 m。土壤密度1 900 kg/m3，土壤地基系數(shù)為40 000 kN/m3，土內(nèi)摩擦角為30°。

計算得出管道所承受載荷：管道自重Gp=5.3×103N/m；埋地段管道上方土層的重力q=3.4×104N/m。

表2 不同懸空長度下管道的最大應(yīng)力和最大位移

4.3 計算結(jié)果

運用ABAQUS有限元軟件計算得出埋地管道在不同懸空長度下管道的最大應(yīng)力，見表2。

根據(jù)表2中的應(yīng)力取值，可見隨著管道懸空長度的增加，管道應(yīng)力也在增加，如圖3所示。該管材管道規(guī)定的最小屈服強度為483 MPa，分析得出管道在懸空長度L=90m時，管道的屈服強度為482.1 MPa，管道材料接近于屈服階段。

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

5 結(jié) 論

懸空管道埋入端采用彈性地基梁模型，通過理論計算得到埋入段上覆土的重力和何種重力載荷。管-土相互作用通過彈簧近似代替，采用GB 50470—2008《油氣輸送管道線路工程抗震技術(shù)規(guī)范》，計算得到懸空管道埋入段在橫向、軸向和豎直方向上土彈簧剛度。運用ABAQUS有限元軟件計算得到該材質(zhì)管道在懸空長度L=90 m時，近似接近于屈服階段。說明采用基于應(yīng)力的管道設(shè)計方法存在局限性，因為管道在屈服階段，管道應(yīng)力變化不顯著，在該階段對管道作出安全評定時應(yīng)采取基于應(yīng)變的設(shè)計方法。

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

[2]王峰會，趙新偉.高壓鋼管黃土塌陷情況下的力學(xué)分析與計算[J].油氣儲運，2004，23（04）：6-8.

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

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

[5]杜景水，馬廷霞，王維斌.基于小撓度理論的懸空管道力學(xué)分析[J].油氣儲運，2009，21（07）：16-18.

[6]童華，祝效華，練章華，等.坍塌和沖溝作用下埋地管道大變形分析[J].石油機械，2007，35（11）：29-32.

[7]博雷西 A P，賽德博 O M，西利 F B，等.高等材料力學(xué)[M].北京：科學(xué)出版社，1987：39-67.

[8]唐永進(jìn).壓力管道應(yīng)力分析[M].北京：中國石化出版社，2003.

[9]馬廷霞，吳錦強，唐愚，等.成品油管道的極限懸空長度研究[J].西南石油大學(xué)學(xué)報，2012 （04）：165-173.

[10]趙林，馮啟民.埋地管線有限元建模方法研究[J].地震工程與工程振動，2001，21（02）：53-57.