蔣 喜 王棠昱 孫 藺 吳 瑕

1.中國石油集團(tuán)工程設(shè)計有限責(zé)任公司西南分公司，四川 成都 610041；2.西南石油大學(xué)石油工程學(xué)院，四川 成都 610500

0 前言

隨著我國管道建設(shè)的發(fā)展，長輸油氣管道不可避免的通過地形陡峭的局部山區(qū)［1］。敷設(shè)在高陡邊坡上的埋地輸氣管道，若管土間的相互摩擦力不能與管道的重量載荷相平衡，將會使邊坡上的管道沿坡道下滑。一般而言，未覆土的高陡邊坡管道，由于缺少頂端覆土的摩擦阻力，極易形成管道下滑，是高陡邊坡管道最危險的特殊邊界情況。但是目前，對高陡邊坡管道的應(yīng)力及位移研究主要集中于邊坡穩(wěn)定性或滑坡后的失效分析［2-4］。因此，對高陡邊坡上的未覆土輸氣管道進(jìn)行應(yīng)力分析研究十分必要。美國COADE公司的CAESAR II，是基于3D梁模型的應(yīng)力分析專業(yè)軟件［5］，可通過合理的工況設(shè)置準(zhǔn)確的模擬非線性載荷［6-7］。本文基于CAESAR II軟件，建立了高陡邊坡未覆土管道的應(yīng)力分析模型，研究了約束條件和邊界條件的設(shè)置方法，討論了加載工況和應(yīng)力校核標(biāo)準(zhǔn)，并通過計算實例驗證了模型的可行性。

1 管道模型的建立

長期實踐和工程應(yīng)用證明，對于管徑與壁厚之比小于100的管道，采用梁單元的精度已足夠［8］。因此管道主體采用3D梁單元模型進(jìn)行模擬，3D模型中定義各坐標(biāo)方向：

a）軸向為平行于管道軸線方向，對應(yīng)坐標(biāo)軸x；

b）縱向為垂直于管道軸線方向，當(dāng)管道水平時與重力方向平行，對應(yīng)坐標(biāo)軸y。

c）橫向為垂直于管道軸線方向，當(dāng)管道水平時垂直于重力方向，對應(yīng)坐標(biāo)軸z。

3D模型中定義各坐標(biāo)方向見圖1。

圖1 各坐標(biāo)方向的定義

管道中的彎頭同樣采用梁單元模型，但是考慮到彎頭處因扁率所導(dǎo)致的應(yīng)力集中，因此特別引入了應(yīng)力增強(qiáng)系數(shù)（Stress Intensification Factor，簡稱SIF），來考慮彎頭處的應(yīng)力集中效果。作為埋地輸氣管道，SIF的計算參照美國ASME B31.8《輸氣和配氣管道系統(tǒng)》標(biāo)準(zhǔn)的附錄E［9］。 建模過程中的管道承壓則作為均布載荷加載于管道內(nèi)壁，除此外還應(yīng)計入壓力荷載所帶來的玻爾登效應(yīng)。

由于梁模型默認(rèn)將重量均分給兩端的2個節(jié)點，因此在劃分管段網(wǎng)格時，若節(jié)點距離過大，將導(dǎo)致過多的重量錯誤的集中于節(jié)點處，使應(yīng)力校核無法通過。因此，文獻(xiàn)［6］建議在建模中：凡是由地面支承的直徑大于300 mm的管道，節(jié)點之間的距離應(yīng)小于等于20倍的管道直徑；而直徑小于等于300mm的管道，節(jié)點之間的距離則應(yīng)大于30倍管道直徑。

2 約束與邊界條件

2.1 土壤約束

未覆土的管道橫臥在管溝中，其管底與土壤連續(xù)接觸，側(cè)邊與管溝壁接觸，管道與管溝作用，其受力：

a）軸向：因管底與土壤連續(xù)接觸，因此受到土壤的摩擦力作用。

b）橫向：受到橫向制止約束，并產(chǎn)生軸向摩擦力。

c）縱向：管道受土壤向上的支承力作用，不會產(chǎn)生垂直土壤面向下的運動。

以上土壤作用可以表達(dá)為：

式中：dy為管道在縱向上的位移，m；f0為單位長度管道上因重力而產(chǎn)生的摩擦力，N；f1為單位長度管道上因橫向推力而產(chǎn)生的摩擦力，N；D為管道外徑，m；μ為土壤與管壁之間的摩擦系數(shù)［10］，按表1取值；α為坡道傾角，管道水平時取0°，管道豎直時90°，deg；ρp為鋼管密度，kg/m3；ρf為管內(nèi)流體的密度，kg/m3；T為管道所受橫向推力，N。

在CAESAR II中，這一約束設(shè)置的方法見圖2。

圖2 CAESAR II中的約束設(shè)置

表1 土壤與管壁之間的摩擦系數(shù)

2.2 邊界條件

根據(jù)未覆土高陡邊坡輸氣管道的實際情況，有三種方法可以簡化邊界條件。一是管道在坡道兩端，即坡道的底部和頂部立即結(jié)束，并設(shè)定為自由邊界；二是根據(jù)斷層邊界理論[11]，認(rèn)為坡道之外的管道將不受坡面管道移動的影響，將管道在坡道兩端處立即結(jié)束，并設(shè)定為固定邊界；三是綜合以上兩方面，將管道在坡道兩端各延長一定距離，再設(shè)定為固定邊界。三種邊界條件見圖3，在實際建模時，自然邊界和固定約束邊界在坡道兩端均只增加1個管段節(jié)點。

圖3 驗證模型的邊界條件

為了驗證不同邊界對模型的影響，采用CAESAR II建立了某管道的應(yīng)力分析模型。管道規(guī)格為Φ1219mm×26.4mm，管段坡面總長650m，坡度為55°。向管道內(nèi)加載水的自重與15MPa的水試壓壓力后，得到管道在三種邊界條件下的軸向及縱向位移見圖4～5，位移方向的規(guī)定見圖1。

三種邊界條件下的計算結(jié)果表明，整條管道都存在軸向位移，這是由非嵌固管道的熱脹效應(yīng)所產(chǎn)生的。管道的軸向位移沿坡道先增大后減小，表明管道沿坡道下滑的趨勢同時受到坡道兩端即坡底和坡頂管道的制約。管道的縱向位移與軸向位移的趨勢相同，不同之處在于管道兩側(cè)出現(xiàn)了縱向位移為0的情況，這是由于建模管道在坡底和坡頂都是水平管段，因此雖然存在軸向位移，但管道實際沿水平方向移動，并沒有縱向位移。

圖4 不同邊界條件下的軸向位移

圖5 不同邊界條件下的縱向位移

分析計算結(jié)果可知：自由邊界管道在軸向和縱向上均出現(xiàn)了遠(yuǎn)大于另兩類邊界的位移。主要原因是忽略了坡頂和坡底管道對坡道上管道的約束作用。而事實上，邊坡兩端所連接的管道可以通過管土間摩擦力來平衡坡道上管道的重力。因此，自由邊界條件忽略這一阻力作用所形成的大位移是不合理的。

另一方面，添加固定約束邊界與加長管道的自由邊界在變化趨勢上是一致的。但加長管道的位移極值較大，且出現(xiàn)極值的位置更靠后。這一現(xiàn)象是坡道管段位移帶動相鄰管道移動而出現(xiàn)的，在固定約束邊界中，無法表達(dá)這一移動。因此，加長管段的固定邊界條件是未覆土高陡邊坡輸氣管道的最佳簡化邊界條件。

R.P.Kennedy于1977年提出斷層邊界理論，指出即使在斷層錯動這樣的大變形條件下，管土之間存在較大相對位移的范圍只有十幾米到三十米左右。為了保證計算的準(zhǔn)確性，規(guī)定直徑大于300mm的管道，加長距離應(yīng)不小于200倍的管道直徑；而直徑小于等于12 in.（DN300）的管道，加長距離不應(yīng)小于300倍管道直徑。

3 加載條件與校核標(biāo)準(zhǔn)

未覆土的高陡邊坡管道在水試壓過程中所受到的基本載荷，見表2。

在管道的應(yīng)力分析中，未覆土管道應(yīng)視作非嵌固管道。根據(jù)ASMEB31.8的相關(guān)規(guī)定，一次、二次應(yīng)力均應(yīng)進(jìn)行校核，需考慮的組合工況及其校核標(biāo)準(zhǔn)見表3。

表2 管道所受基本載荷

表3 組合工況與校核標(biāo)準(zhǔn)

應(yīng)注意在計算熱應(yīng)力時采用L2-L3而非直接使用T1，其主要的原因是當(dāng)管道存在非線性約束時，單獨加載T1與采用OPE-SUS間的位移結(jié)果可能發(fā)生反向。為了保證求解結(jié)果的準(zhǔn)確性，CAESAR II一般在求解EXP工況時采用OPE-SUS的組合工況。

未覆土的管道在水重及水壓的共同作用下，可能被拉動下滑。因此，當(dāng)管道中的一次、二次應(yīng)力滿足要求后，還要對管道的軸向和縱向負(fù)位移進(jìn)行校核。

4 應(yīng)用實例

某天然氣輸送管道規(guī)格為Φ406.4mm×14.3mm，材料為L360，管道總長626m，落差248m，最大坡度達(dá)61.58°。管內(nèi)介質(zhì)密度為80 kg/m3，所在地區(qū)等級為一級，安裝溫度10℃，運行溫度20℃，運行壓力10MPa，水試壓壓力為15MPa。根據(jù)管道基本信息及安裝圖，在管道兩側(cè)添加了加長管道后的固定邊界?；贑AESAR II所建成的分析模型見圖6，計算得到的各工況條件下的應(yīng)力和位移極值見表4。

由表4可知，水試壓條件下的應(yīng)力極值遠(yuǎn)大于其它加載工況。為了進(jìn)一步分析水試壓條件下管道應(yīng)力的變化趨勢，特將該工況下的應(yīng)力值進(jìn)行排序，表5為應(yīng)力最大的5個節(jié)點數(shù)據(jù)。

圖6 建成的未覆土高陡邊坡模型

表4 各工況條件下的應(yīng)力和位移極值

表5 水試壓條件下的應(yīng)力計算結(jié)果

本實例中管道的最大軸向位移為-4.42 cm，最大縱向位移為-1.92 cm（位移方向的規(guī)定可見圖1）。管道的最大位移出現(xiàn)在節(jié)點281處，位于管道最長且最陡峭邊坡的頂端，見圖7。

綜合管道走向與應(yīng)力、位移校核結(jié)果，可以得出：

a）管道在未覆土條件下的各工況中沒有出現(xiàn)應(yīng)力超限的情況，其中應(yīng)力最大的為水試壓工況，最大應(yīng)力值占許用應(yīng)力的60.7%，最大應(yīng)力節(jié)點為139，位于坡底平臺與下一段下坡管道的彎頭處。

b）根據(jù)表5可知，水試壓條件下的最大應(yīng)力位置均為彎頭處，但其應(yīng)力值的大小并不只取決于彎頭的角度，還與前后連接的管道是否水平及坡道長短有關(guān)，因此最大角度的坡道彎頭并不一定是最大應(yīng)力點。

c）管道的最大軸向位移為-4.42 cm，最大縱向位移為-1.92 cm。管道的最大位移出現(xiàn)在節(jié)點281處，位于管道最長且最陡峭邊坡的頂端，由此可見，邊坡頂?shù)墓艿朗俏灰谱畲蟆?/p>

圖7 節(jié)點281所在位置

5 結(jié)論

針對未覆土條件下高陡邊坡上天然氣管道的實際情況，提出了該類管道的應(yīng)力分析模型。通過比選三類不同的簡化邊界條件，確定將坡道兩側(cè)有水平增長段的固定約束作為模型的邊界條件。建立了四種加載工況，分別為水試壓工況、操作工況、持續(xù)工況和熱應(yīng)力工況，并基于美國ASMEB31.8標(biāo)準(zhǔn)的要求提出了各類工況的校核方法。通過計算實例驗證了應(yīng)力計算模型和校核方法的可行性。

［1］吳 宏，李 波.西氣東輸工程隧道內(nèi)管道施工新方法［J］.油氣儲運，2003，2（6）：47-52.W u Hong，Li Bo.Latest Methods Used to Build the Pipe in Land Tunnelsin W est-EastGasProject［J］.Oil&GasStorage and Transportation，2003， 2（6）:47-52.

［2］王芝銀，袁鴻鵠，王 怡，等.管道穿越土質(zhì)邊坡斜豎井設(shè)計參數(shù)的確定方法［J］.中國石油大學(xué)學(xué)報 （自然科學(xué)版），2010，34（1）： 105-108.Wang Zhiyin，Yuan Honghu，Wang Yi，et al.New Method of Determining Design Parametersof Inclined Shaft for Oil-Gas Pipeline Crossing Along Soil Slope ［J］.Journalof China University of Petroleum（Edition ofNaturalScience），2010，34（1）:105-108.

［3］俞火明.杭-甬天然氣管道沿線某高邊坡穩(wěn)定性分析與加固方案評價［J］.中國港灣建設(shè)，2011，176（5）： 47-50.Yu Huom ing.Stability Analysisand Reinforcement Scheme of aHigh Slopealong Hang-Yong NaturalGasPipeline［J］.ChinaHarbour Engineering， 2011， 176（5）:47-50.

［4］周曉瑩，郭永華，呂曉華，等.地質(zhì)滑坡條件下不同因素對油氣管道危險性影響分析［J］.工業(yè)安全與環(huán)保，2012，38（5）：42-44.Zhou Xiaoying，Guo Yonghua，Lv Xiaohua，etal.Study on Effect of Different Factorson Pipelines Risk under Landslide［J］.Industrial Safety and Environmental Protection， 2012， 38（5）:42-44.

［5］龍忠輝.CAESARⅡ管道應(yīng)力分析軟件開發(fā)應(yīng)用［J］.化工設(shè)備與管道，2001，38（3）：50-54.Long Zhonghui.Developmentand Application of Piping Stress AnalysisSoftware CAESAR Ⅱ ［J］.ProcessEquipment&Piping， 2001， 38（3）:50-54.

［6］GeorgeM.管道應(yīng)力分析理論指南［M］.王 昕，譯.COADE Inc中國技術(shù)服務(wù)與培訓(xùn)中心，1998.22-27.George M.Pipe Stress Analysis Theory Guide［M］.COADE Inc China Technical Service and Training Center，1998.22-27.

［7］Stumpf P D，Ainsworth H A，Appleton R A，et al.Technical ReferenceManual［M］.COADE Inc Press，2003.128-221.

［8］Audibert J，Nyman K.Coefficients of Subgrade Reactions for the Design of Buried Piping ［C］.2ndASCESpecialty Conference on Structure Design of Nuclear Plant Facilities，LA，1975.235-244.

［9］Nance A D，Hayden L E，Stumpf PD，et al.B31.8-2010 Gas Transm ission and Distribution Piping Systems［S］.The American Society ofMechanicalEngineersPress，2010.21-27.

［10］帥 健，于桂杰.管道及儲罐強(qiáng)度設(shè)計［M］.北京：石油工業(yè)出版社，2006.10-11.Shuai Jian，Yu Guijie.Stress Design for Pipelines and Tanks［M］.Beijing:Petroleum Industry Press， 2006.10-11.

［11］KharaitiL，George M，Nyman D J.GuidelinesFor the Seism ic Design of O il and Gas Pipeline System ［M］.The American Society ofCivilEngineersPress，1984.70-77.