任 帥，馬貴陽，王海燕

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

基于ANSYS的隧道內(nèi)輸氣管道應(yīng)力分析

任 帥，馬貴陽，王海燕

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

隧道內(nèi)的彎管處通常是應(yīng)力集中的地方，使得管道處于危險甚至失效狀態(tài)。為了研究隧道內(nèi)彎管在不同情況下的受力，利用有限元分析軟件ANSYS建立彎管的有限元模型，明確了隧道內(nèi)垂直L型和Z型彎管的應(yīng)力分布情況和應(yīng)力集中點。研究結(jié)果表明：在試壓和運行狀態(tài)下，彎頭處所受的應(yīng)力最大，管道內(nèi)壓是引起管道應(yīng)力的主要因素。隧道內(nèi)管道最大應(yīng)力隨管道半徑的增大而逐漸增大，呈近似線性關(guān)系。不同形狀的彎管結(jié)構(gòu)分布會對彎管的應(yīng)力分布產(chǎn)生影響。L型彎管的最大應(yīng)力隨兩臂長的增加先增大后減小，Z型彎管的最大應(yīng)力隨相互平行的兩彎管的增大而逐漸增大。

隧道；輸氣管道；試壓工況；有限元分析；半徑；臂長

隨著我國輸氣管道建設(shè)高峰時期的到來，在山區(qū)敷設(shè)管道時，為了減少植被破壞、防止水土流失等自然災(zāi)害的發(fā)生，隧道穿越形式被廣泛采納[1]。隧道內(nèi)的油氣管道的設(shè)計目前沒有相對可靠的理論依據(jù)，大部分設(shè)計必須依靠經(jīng)驗參數(shù)，此外，管道在隧道內(nèi)的應(yīng)力分布情況模糊，使得管道的運行存在很大的隱患[2]。導(dǎo)致管道失效的因素有管道腐蝕，管道疲勞，設(shè)計不合理等，除此之外還有一個經(jīng)常被忽略的重要原因，就是在隧道內(nèi)存在一些特殊管段，諸如彎頭等常常會產(chǎn)生應(yīng)力集中，這使得管道不能滿足強度要求而失效[3-6]。國外的管道在“L”和“Z”段可能產(chǎn)生橫向位移處增加軟泡材料，如廢聚氨脂，泡沫塑料。而目前我國并沒有采取相應(yīng)措施，在小管徑和溫差較低時彎頭部位應(yīng)力未達到材料的強度極限，仍可正常運行，但當(dāng)運行管道管徑，溫差較大時會造成管道破裂[7,8]。

管道施工中一個非常重要的程序是管道試壓，不僅是對管道工程在生產(chǎn)中能否在設(shè)計壓力下安全運行的一次檢驗，也是對管道施工質(zhì)量，材料性能，管道整體性能的一次綜合檢驗[9,10]?！对秃吞烊粴廨斔凸艿来┛缭皆O(shè)計規(guī)范》6.2.1中規(guī)定，大中型水域，鐵路及Ⅱ級以上的公路穿越管段，必須獨立進行強度試壓和嚴密性試驗，合格后再與相鄰管段連接。本文利用有限元軟件ANASYS對隧道內(nèi)彎管建立有限元模型，分析由于半徑,結(jié)構(gòu)參數(shù)的不同而使彎管受到不同的應(yīng)力應(yīng)變，找出危險點，為設(shè)計施工提供一定理論指導(dǎo)。

1 管道的應(yīng)力分類及校核[11,12]

因為管道上有載荷作用，就會產(chǎn)生各種各樣的應(yīng)力。由于外部和內(nèi)部作用的影響，諸如壓力、重力、溫度都可能引起管道相對應(yīng)應(yīng)力，這些應(yīng)力又可分為一次應(yīng)力，二次應(yīng)力和峰值應(yīng)力。

由管道的內(nèi)壓、自重和其他外載產(chǎn)生的應(yīng)力稱為一次應(yīng)力。于浦義等編寫的石油化工壓力管道設(shè)計手冊中規(guī)定：內(nèi)壓、自重和外載產(chǎn)生的環(huán)向應(yīng)力、軸向應(yīng)力和剪應(yīng)力等應(yīng)滿足外部力和力矩的平衡條件。如果應(yīng)力過大，管道產(chǎn)生的變形將無限制的增大（即沒有自限性），最終導(dǎo)致破壞。一次應(yīng)力的大小是衡量管道能否安全運行的標(biāo)準(zhǔn)之一，為了確保一次應(yīng)力不變，安裝條件、內(nèi)壓、溫度和外載必須保持不變。一次應(yīng)力值小干或等于操作溫度下材料的許用應(yīng)力， 則管道可以安全運行。

由于管道變形受阻而引起的正應(yīng)力和剪應(yīng)力稱為二次應(yīng)力，其不是直接與外力相平衡，而是通過管道各部分變形來適應(yīng)的。同時，二次應(yīng)力作用的范圍限制在局部范圍內(nèi)。熱脹推力的作用可以引起管道局部屈服，從而產(chǎn)生少量塑性變形，此時，椎力就不會再增加，塑性變形也不再發(fā)展，這種現(xiàn)象稱為有自限性。

根據(jù)ASME．B31.8《氣體運輸和分配系統(tǒng)的管道標(biāo)準(zhǔn)》，因為隧道內(nèi)的管道處于比較封閉的環(huán)境，并建有固定墩，支墩等輔助設(shè)施，所以靜力分析主要校核自重W，內(nèi)壓P，溫度T。管道總應(yīng)力疊加計算式（1）：

式中L1為計算管道系統(tǒng)對邊界的推力。

一次應(yīng)力疊加計算式：L2=W+P

式中L2為管道系統(tǒng)中各點的一次應(yīng)力值。

二次應(yīng)力疊加計算式：L3= L1-L2

式中L3為管道系統(tǒng)中各點的二次應(yīng)力值。

輸氣管道穿越工程在正常投產(chǎn)運行前需要進行單獨試壓，根據(jù)SY/T0015.1-1998《原油和天然氣輸送管道穿越工程設(shè)計規(guī)范—穿跨工程》中規(guī)定，油氣輸送管道的穿越管段試壓壓力，試壓時間和試壓介質(zhì)分別按國家標(biāo)準(zhǔn)GB50251和GB50253規(guī)定執(zhí)行，強度驗算壓力按下（2）式計算，

式中

:Pt—試壓壓力，MPa；

δ—管道設(shè)計壁厚，mm；

σs—管道規(guī)定的屈服極限，MPa；

D—管道外徑，mm。

2 管道模型的建立及應(yīng)力分析

2.1 模型的建立

隧道內(nèi)管道走向的平面分布如圖1所示。隧道內(nèi)彎管可以分為垂直型L型和Z型。管道直徑813 MPa，壁厚22 mm，運行壓力為6.4 MPa，根據(jù)式1，試壓壓力24 MPa，試壓介質(zhì)為水。管內(nèi)溫度為50 ℃，管外平均溫度為10 ℃在ANSYS中選用pipe16，pipe18單元建立彎管有限元模型。定義第一類第二類材料的彈性模量為206 GPa，密度為7 850 kg/m3,泊松比為0.3。用命令“RUN”定義管道分別沿Y軸的延伸量為4.6、0.4 m，沿X軸的延伸量為0.4、4.6 m。用命令“BEND”定義彎管處半徑為0.3 m，建立垂直L型彎管。用命令“RUN”定義管道分別沿X軸的延伸量為4.6、0.4 m，沿Y軸的延伸量為0.4、3.2、0.4 m，沿X軸的延伸量為0.4，4.6 m。用命令“BEND”定義兩個彎管處半徑為0.3 m，建立垂直Z型管道。通過以上步驟生成的管道彎管模型如圖2所示。

圖1 隧道穿越管道走向平面分布圖Fig.1 The Distribution map of thepipeline thruogh tunnel

圖2 彎管有限元模型Fig.2 The finite element model of bending pipeline

2.2 加載求解

定義Y方向的重力加速度為9.8 m/s2，定義兩端點所有節(jié)點自由度。加載完成后進入求解器求解。通過通用后處理器顯示應(yīng)力分布情況。垂直L型Z型彎管在試壓工況下的應(yīng)力分布如圖3所示。

圖3 彎管在試壓工況下的應(yīng)力分布Fig.3 The stress distribution of bending pipeline under pressure testing condition

3 數(shù)值模擬結(jié)果分析

3.1 運行工況不同對最大應(yīng)力的影響

通過上述分析，兩類彎管在不同運行工況下X方向的最大應(yīng)力值如下表1所示。從表1可知，試壓工況下“L”，“Z”型管道在彎管處的應(yīng)力最大，“Z”型彎管處的最大應(yīng)力大于“L”型的，且最大應(yīng)力發(fā)生在“Z”型管靠下邊的彎管處。隧道內(nèi)彎管在設(shè)計安裝時必須考慮試壓工況，使得管道能夠滿足試壓時的強度要求。

表1 不同運行工況彎管應(yīng)力分布情況表Table 1 The stress distribution of bending pipeline under different conditions

3.2 管道半徑對最大應(yīng)力的影響

由于在試壓工況下管道需承受的應(yīng)力最大，所以下面只考慮在試壓工況下半徑對最大應(yīng)力的影響。圖4為管道半徑變化對最大應(yīng)力的影響，由圖可知，最大應(yīng)力隨半徑呈近似線性關(guān)系。當(dāng)隧道內(nèi)管道為大半徑時需考慮在彎管處增加保護措施以防管道破裂。

圖4 管道半徑變化對最大應(yīng)力的影響Fig.4 The effect of pipe radius on the maximum stress

3.3 管道結(jié)構(gòu)尺寸對最大應(yīng)力的影響

改變“L”，“Z”型彎管的臂長，但保持兩臂相等，“Z”型垂直臂保持4 m不變，相互平行的兩臂變換。臂長變化隨最大應(yīng)力值的關(guān)系如下表2，由表可知，“L”型彎管，隨著臂長變長先增大后減小?！癦”型彎管隨臂長增加最大應(yīng)力也增加。當(dāng)“Z”型彎管兩平行臂長保持10 m不變，垂直臂長由4米逐漸增加到10 m，最大應(yīng)力逐漸減小，見圖5。

表2 水平臂長變換對最大應(yīng)力的影響Table 2 The parallel arm length influence on the maximum stress distribution of bending pipeline

圖5 垂直臂長對“Z”型彎管的應(yīng)力的影響情況Fig.5 Influence of vertical arm length on the stress distribution of “Z” bending pipeline

4 結(jié) 論

（1）在試壓工況下，垂直L型與Z型彎管所受應(yīng)力最大且最大應(yīng)力產(chǎn)生的位置都在彎管處，加之彎管處一般有焊縫，所以彎頭處為危險面。

（2）隨著隧道內(nèi)管道半徑增加，彎管處的最大應(yīng)力值增加。

（3）對于“L”型彎管，最大應(yīng)力隨兩臂長的增大先增大后減小。對于“Z”型彎管，最大應(yīng)力隨相互平行臂長的增加而增大。在實際設(shè)計安裝時，根據(jù)隧道實際深度、走向結(jié)合以上結(jié)論，調(diào)整特殊管段處的結(jié)構(gòu)尺寸，必要時在彎頭處添加軟泡材料，使隧道內(nèi)管道可以安全運行。

［1］ 陸悅.川氣東送管道工程山嶺隧道內(nèi)管道安裝設(shè)計[J].科技與企業(yè)，2011(4):60-61.

［2］ 岳進才.壓力管道技術(shù)[M].北京：中國石化出版社，2006：13-16.

［3］ 吳曉楠，鮮燕，劉源海，等.清管過程中隧道內(nèi)輸氣管道應(yīng)力分析[J]天然氣與石油，2012,30（2）:1-3.

［4］向勇.西氣東輸中衛(wèi)黃河隧道設(shè)計[J].石油工程建設(shè)，2011，37(3):51-57.

［5］彭星煜,梁光川,張鵬,等.人工爆破地震作用下輸氣管道動力響應(yīng)分析[J].天然氣工業(yè)，2012,32（11）：81-84.

［6］唐永進. 壓力管道應(yīng)力分析[M]. 2版.北京：中國石化出版社，2009：51-66.

［7］郭瑞平，李廣信，陳輪.“L”型直埋供熱管道受力規(guī)律數(shù)值模擬[J].計算結(jié)構(gòu)力學(xué)及其應(yīng)運，1996,13（3）：327-332.

［8］郭瑞平，李廣信.“Z”型直埋供熱管道受力特性的研究[J].煤炭學(xué)報，1995,20（6）:658-663.

［9］吳曉南，舒浩紋，昝林峰，等.試壓工況下盾構(gòu)隧道內(nèi)輸氣管道應(yīng)力分析[J] 天然氣工業(yè)，2013,33（3）:73-77.

［10］余運復(fù)，張繼霞，吳俊松，等.國內(nèi)外輸油氣站工藝管道試壓標(biāo)準(zhǔn)差異分析[J].全面腐蝕控制，2014(6)：74-78.

［11］于浦義,張德姜,唐永進,等.石油化工壓力管道設(shè)計手冊[M].北京:化學(xué)工業(yè)出版社,2010.

［12］ 隋之鋒. 輸氣管道特殊管段應(yīng)力分析及安全性研究[D]. 中國石油大學(xué)(華東), 2009. DOI:10.7666/d.y1543071.

Stress Analysis of Gas Pipeline in Tunnel Based on ANSYS

REN Shuai，MA Gui-yang，WANG Hai-yan
（Liaoning Shihua University, Liaoning Fushun 113001，China）

The elbows in the tunnel are usually the stress concentration point. They make the pipeline be in dangerous or even failure states. In order to study the stress condition in tunnel under different conditions, the finite element analysis software ANSYS was used to build the finite element model to carry out the numerical simulation in order to find out the stress distribution and the stress concentration points in type “L” and type “Z” pipelines. The results show that, under pressure testing and operating condition, the pressure of the elbow is the highest. The inner pressure in the pipeline is the main factor to cause the stress. The maximum stress of pipeline in the tunnel gradually increases along with the pipe radius increasing, and there is an approximate linear relationship between them. The different structure of elbow can affect the stress distribution, the maximum stress on “L” pipe first increases and then decreases with the increase of arm length, but it always increases on “Z” pipe.

Tunnel; Gas pipeline; Pressure testing condition; Finite element analysis；Radius; Arm length

TE 832

： A

： 1671-0460（2015）11-2652-03

2015-04-28

任帥（1989-），男，內(nèi)蒙古自治區(qū)鄂爾多斯市人，在讀碩士，現(xiàn)就讀于遼寧石油化工大學(xué)，研究方向：管道泄漏分析。E-mail：943178963@qq.com。