大口徑曲線管道在運(yùn)營(yíng)階段的溫度應(yīng)力分析

秋水

1 引言
市政工程、石油工程中經(jīng)常采用管道來(lái)進(jìn)行水、
天然氣、石油等的輸送。這類管道在運(yùn)營(yíng)過程中，由
于氣溫變化或其他原因( 如原油預(yù)熱) 等造成管內(nèi)
流體的溫度變化通常使管壁內(nèi)產(chǎn)生不可忽視的溫度
應(yīng)力。對(duì)于管道的溫度應(yīng)力，國(guó)外的學(xué)者Iimura 和
Nishio［1］進(jìn)行了試驗(yàn)研究，結(jié)果表明在高溫流體作
用下管壁將產(chǎn)生壓應(yīng)力，而當(dāng)管道冷卻后，管壁殘留
可觀的拉應(yīng)力。Bakic［2］等提出了考慮變溫歷史的
管壁溫度應(yīng)力計(jì)算公式，結(jié)果表明變溫歷史造成的
溫度應(yīng)力影響管道使用壽命。Taler 和Lubecki［3］通
過有限差分法獲得了高溫蒸汽管道內(nèi)流體和管道的
溫度場(chǎng)分布。國(guó)內(nèi)在較早時(shí)候也對(duì)該領(lǐng)域進(jìn)行了研
究，如華東石油學(xué)院［4］對(duì)埋地重油熱輸管道的溫度
應(yīng)力進(jìn)行了現(xiàn)場(chǎng)測(cè)試，獲得軸向應(yīng)力和軸向位移的
分布規(guī)律，并提出軸向應(yīng)力和環(huán)向應(yīng)力的計(jì)算公式。
蔡強(qiáng)康，呂英民［5］采用數(shù)值方法對(duì)埋地?zé)彷敼芫�的
彎管溫度應(yīng)力進(jìn)行了影響參數(shù)分析，結(jié)果表明彎管
曲率半徑與彎管外徑之比，彎管角度，側(cè)向土壓力等
是重要的參數(shù)。郭瑞平等［6］對(duì)L 形直埋供熱管道
進(jìn)行數(shù)值模擬，得到內(nèi)力和位移沿管長(zhǎng)的分布規(guī)律。
已有研究主要集中于管徑為200 ～ 800mm 的小口徑
管道，大口徑管道較少涉及。實(shí)際上，大口徑管道在
運(yùn)營(yíng)時(shí)所承受的溫度應(yīng)力要比小口徑的大，因此更
有必要對(duì)工程中常用的大口徑管道進(jìn)行溫度應(yīng)力分
析。另外，實(shí)際工程中為避讓障礙物，經(jīng)常采用曲線
線形進(jìn)行管道敷設(shè)，由于曲率影響，這類曲線形管道
的溫度應(yīng)力要比直線管道的大。已有研究主要對(duì)直
建筑科學(xué)第28 卷
線形管道［1-3］或者管道的彎頭段分析［4-6］，而對(duì)曲線
形管道甚少涉及。由上可見，目前對(duì)于大口徑曲線
形管道的溫度應(yīng)力分析非常缺乏，工程上幾乎沒有
可供借鑒的例子和建議，有必要進(jìn)行深入的分析研
究。
本文以青草沙水源地原水工程嚴(yán)橋支線工程曲
線段頂管工程為背景，采用數(shù)值方法研究了大口徑
曲線管道在運(yùn)營(yíng)階段下的溫度應(yīng)力，并進(jìn)行了影響
參數(shù)分析，分析結(jié)果可為類似工程提供建議。
2 工程概況
青草沙水源地原水工程管道采用頂管法施工，
設(shè)計(jì)輸水規(guī)模為719 萬(wàn)m3 /d。嚴(yán)橋支線工程是其
中最大的輸水支線工程，設(shè)計(jì)規(guī)模為440 萬(wàn)m3 /d，
采用二根內(nèi)徑為3600mm 的鋼頂管進(jìn)行頂管施工。
J42B 至J45 井管段為避開障礙物，采用曲線形的管
道，其最小曲率半徑約為880m，弧長(zhǎng)約為170m［7］，
該段曲線管道為本文的計(jì)算分析對(duì)象。
3 有限元分析模型及參數(shù)
按照實(shí)際尺寸進(jìn)行建模，管道的中心軸線長(zhǎng)度
為170m，曲率半徑為880m，內(nèi)徑3600mm，壁厚
34mm。采用20 節(jié)點(diǎn)的六面體二次實(shí)體單元，整個(gè)
三維有限元模型包括24480 個(gè)單元，122664 個(gè)節(jié)
點(diǎn)。沿管壁厚度方向布置三個(gè)單元�？紤]到管道長(zhǎng)
度較大，所以其兩端可近似認(rèn)為沿軸向有固定約束。
有限元模型的總體視圖和局部視圖如圖1，2
所示。
圖1 有限元模型總體視圖
圖2 有限元模型的局部視圖
管道材料是鋼材，采用線彈性模型，密度為
7800kg /m3，線膨脹系數(shù)為1. 3 × 10 － 5，彈性模量為
200GPa，泊松比為0. 25。
4 管道溫度應(yīng)力的參數(shù)影響分析
在數(shù)值計(jì)算中，對(duì)如下幾組參數(shù)進(jìn)行參數(shù)影響
分析: ( 1) 管道內(nèi)外壁溫差值5℃，10℃，15℃，20℃
四種情況( 內(nèi)壁溫度比外壁高) ; ( 2) 管道整體升溫
5℃， 10℃， 15℃， 20℃四種情況; ( 3) 整體升溫20℃
時(shí)，管道內(nèi)徑3600mm 而壁厚采用48、34、24mm 三
種情況，即徑厚比分別為75、105、150; ( 4) 整體升溫
20℃時(shí)，管道軸線曲率半徑R 采用800m，880m，
3000m， 5000m 和無(wú)窮大( 即直線) 五種情況。
分析的結(jié)果包括: ( 1) 彎曲內(nèi)側(cè)的外壁中部應(yīng)
力σA
( 圖3) ; ( 2) 彎曲內(nèi)側(cè)的內(nèi)壁中部應(yīng)力σB
; ( 3)
彎曲外側(cè)的內(nèi)壁中部應(yīng)力σC
; ( 4) 彎曲外側(cè)的外壁
中部應(yīng)力σD
; ( 5) 彎曲內(nèi)側(cè)的內(nèi)外壁中部應(yīng)力差σE
= σB － σA
; ( 6) 彎曲外側(cè)的內(nèi)外壁中部應(yīng)力差σF =
σC － σD。
圖3 管道示意圖
4. 1 管道內(nèi)外壁溫差
由圖4 可見，由于曲率存在，彎曲內(nèi)側(cè)的σA
和
σB
均呈中間受壓，兩端部受拉的模式。σA
以受拉
為主，隨溫差的增大，σA
在管道兩端部隨之增大，在
中部基本不變; σB
由于內(nèi)壁溫度高而以受壓為主，
隨溫差的增大，σB
在管中部增加較快，在兩端增加
較慢。彎曲外側(cè)的σC
和σD
均呈中間受拉、兩端部
受壓的模式，都隨著溫差的增大而增大，且由于內(nèi)壁
溫度較高，σC
以受壓為主，而σD
以受拉為主。由表
1 可見，曲線管道內(nèi)側(cè)和外側(cè)的內(nèi)外壁應(yīng)力差均隨
著溫差的增大而呈線性增大，且外側(cè)的應(yīng)力差比內(nèi)
側(cè)的應(yīng)力差略大。
4. 2 管道整體升溫值
由圖5 可見，在整體升溫時(shí)，管壁的應(yīng)力模式與
溫差情況一樣，但應(yīng)力幅值較大。彎曲內(nèi)側(cè)的σA