基于ANSYS的管道變形應(yīng)力分析

趙福興　周鵬飛　錢英豪

（江蘇省特種設(shè)備安全監(jiān)督檢驗研究院蘇州分院）

?

基于ANSYS的管道變形應(yīng)力分析

趙福興*周鵬飛錢英豪

（江蘇省特種設(shè)備安全監(jiān)督檢驗研究院蘇州分院）

摘要對某公司埋地管道受外部載荷作用發(fā)生變形這一問題，借助于ANSYS軟件的應(yīng)力分析功能，通過施加位移載荷求出管道變形后的應(yīng)力分布，為管道的安全評價提供依據(jù)。

關(guān)鍵詞管道應(yīng)力分析變形位移有限元

*趙福興，男，1967年生，工程師。蘇州市，215031。

0　引言

隨著工業(yè)技術(shù)的發(fā)展，管道在供熱、供煤氣、給水、排水、農(nóng)業(yè)灌溉、水力工程、各種工業(yè)裝置以及長距離輸送石油和天然氣中的應(yīng)用越來越廣。管道的安全與人民群眾的生命和財產(chǎn)安全息息相關(guān)。為了保障管道的安全，我國從管道的設(shè)計、制造、安裝到使用，制訂了一系列的標(biāo)準(zhǔn)。但管道在安裝及使用的過程中，往往會出現(xiàn)一些設(shè)計時未知的外部載荷的作用，使其發(fā)生變形。如果管道變形相對較大，這時候就需要對管道的安全性進行評價，并以此來決定管道是否可以繼續(xù)使用。由于外部的載荷往往比較復(fù)雜，導(dǎo)致管道變形也很復(fù)雜。因此，在管道的使用過程中如何對其發(fā)生的變形進行應(yīng)力及安全性評價，成為亟待解決的問題。

本文對某公司埋地管道受外部載荷（堆土）作用發(fā)生變形這一問題，采用ANSYS軟件對管道的變形進行數(shù)值模擬，以便了解管道的真實的應(yīng)力分布情況，為管道的安全評估及預(yù)測提供依據(jù)。

1 ANSYS管道計算介紹

ANSYS是融結(jié)構(gòu)、熱、流體、電磁、聲學(xué)于一體的大型通用有限元分析軟件，其公司也是目前世界上CAE行業(yè)中最大的公司。該軟件具有完備的前處理功能、強大的求解器、方便的后處理功能以及多種實用的二次開發(fā)功能等，被廣泛應(yīng)用于航空航天、機械制造、土木工程等眾多工業(yè)領(lǐng)域及科學(xué)研究中[1]。通過ANSYS建立所求問題的模型，劃分單元，求解有限個數(shù)值，就可以近似模擬真實環(huán)境的未知量。借助于ANSYS進行結(jié)構(gòu)強度和剛度分析，可以清楚地了解整個結(jié)構(gòu)的應(yīng)力分布以及變形情況。

目前用ANSYS對管道進行分析大都是基于實體單元，這主要是因為采用SOLID實體單元，可以很好地模擬局部應(yīng)力分布狀況，從而對管道的局部進行分析。但對于一根很長的管道進行整體分析，采用SOLID實體單元無疑會使計算的量非常龐大，而且難以保證計算的準(zhǔn)確性。為此，ANSYS針對管道中各種常見的情況，開發(fā)了六種管單元，即三維彈性直管單元PIPE16、三維彈性T形管單元PIPE17、彈性彎管單元PIPE18、塑形直管單元PIPE20、塑形彎管（彎管頭）單元PIPE60以及能夠模擬海洋波浪、水流和纜索的單元PIPE59。管單元的出現(xiàn)，極大地方便了ANSYS在管道分析中的應(yīng)用。隋之鋒、郝點、陳海峰等人借助于ANSYS的管單元對兩端固定的水平直管、垂直L形彎管和Z形管三種典型厚壁管道在重力下的變形及應(yīng)力分布進行了模擬[2]；劉純、胡波濤等人借助于ANSYS PIPING MODELS對拔伯克動力有限公司362.5 MW機組汽缸至除氧器抽氣管道進行了建模和應(yīng)力計算[3]。這些研究結(jié)果均表明，管單元可以準(zhǔn)確、方便地解決管道的分析問題。但目前基于管單元的分析還不多，許多問題還有待進一步探索和檢驗。

2　實例分析

2.1問題描述

某段天然氣高壓埋地管線（PN4.0，DN600）安裝完畢后，由于堆土導(dǎo)致管道隨土體一起發(fā)生了水平偏移，現(xiàn)需對其投用的可行性進行分析評定。管道發(fā)生的位移（或變形）如圖1所示。

圖1　管道變形示意圖

天然氣高壓管線計算參數(shù)如表1所示。材料在計算溫度（工況）下的力學(xué)性能參數(shù)如表2所示。

表1　管線計算參數(shù)

表2　材料在計算工況下的力學(xué)性能參數(shù)

2.2問題分析

由于堆土對管道的作用力比較復(fù)雜，無法采用施加外力載荷的方法來求解。由文獻(xiàn)[4]可以知道：對于彈性連續(xù)體的有限元，其等效節(jié)點力與位移，以及應(yīng)力與位移之間存在如下的關(guān)系：

式中F——作用在單元上的靜力等效力；

k——單元剛度矩陣；

δ——單元的節(jié)點位移；

F0——由初應(yīng)變引起的力。

式中σ——作用在節(jié)點上的應(yīng)力；

D——彈性矩陣；

B——單元應(yīng)變矩陣；ξ0——初始應(yīng)變。

由以上兩式可以看出，等效節(jié)點力與位移存在著對應(yīng)關(guān)系，通過施加位移也可以找到對應(yīng)的應(yīng)力。因此可采用施加位移載荷的方法來求解應(yīng)力強度。位移載荷可以通過現(xiàn)場測量獲得。為保證計算結(jié)果更加逼近實際情況，在保證測量每個變形極大值點的同時，應(yīng)盡量多地測量位移點。位移測試數(shù)據(jù)取點及位移量如圖1及表3所示。

表3　管道測量位移偏移量

2.3有限元分析

2.3.1模型的建立

考慮到所分析的管道為一直管，且變形較小，因此采用PIPE16彈性直管單元。如果應(yīng)力較大，超過屈服強度，可再考慮采用PIPE20塑性直管單元。根據(jù)表2輸入材料的彈性模量和泊松比，根據(jù)表1輸入管道的基本參數(shù)，于是即可建立幾何模型。選擇的分析標(biāo)準(zhǔn)為ASME B31.1《動力管道》。

2.3.2施加載荷

管道兩端處施加X、Y、Z約束；管道內(nèi)表面施加壓力載荷4.0 MPa；在圖1所示的測量點施加位移載荷，位移量根據(jù)表3選取。

管道模型及約束載荷如圖2所示。

圖2　管道模型及約束載荷

2.3.3求解及查看結(jié)果

載荷施加完后進入求解器進行求解。求解完之后進入后處理器顯示應(yīng)力分布和位移分布云圖，分別如圖3、圖4、圖5、圖6、圖7所示，其中變形圖為放大4.7倍之后的效果圖。

圖3　整體變形及應(yīng)力云圖

圖4　1段應(yīng)力及變形云圖

圖5　2段應(yīng)力及變形云圖

圖6　3段應(yīng)力及變形云圖

圖7　4段應(yīng)力及變形云圖

由以上應(yīng)力云圖可以看出：（1）管道最大應(yīng)力為343.968 MPa，小于屈服強度415 MPa，發(fā)生彈性變形，所以采用PIPE16彈性直管單元是合適的。（2）在管道發(fā)生突變的地方，局部應(yīng)力較大，而在管道變形比較平滑的區(qū)域，應(yīng)力較小。（3）施加位移載荷之后的變形圖基本保持與實際變形位移一致，其應(yīng)力云圖可以近似真實地反映實際應(yīng)力分布及大小。

3　結(jié)論

（1）參照GB 50251—2003《輸氣管道工程設(shè)計規(guī)范》4.2節(jié)中的地區(qū)等級劃分及強度設(shè)計系數(shù)，該段管道所屬地區(qū)等級為二級地區(qū)，強度設(shè)計系數(shù)取0.6，常溫下材料的許用應(yīng)力為249 MPa，管道發(fā)生位移之后的最大Misses應(yīng)力為343.968 MPa，

大于249 MPa，所以管道強度要求不合格。

（2）根據(jù)以上分析，管道在突變地區(qū)局部應(yīng)力較大，而大部分管道比較平滑的區(qū)域應(yīng)力均遠(yuǎn)小于許用應(yīng)力，所以建議對突變地區(qū)管道進行更換和重新布管，改變受力狀況，以確保管道的安全運行。

（3）通過ANSYS模擬的變形與堆土之后的變形情況基本一致。實例表明：借助于ANSYS強大的管道分析功能，可以近似真實地反映管道的變形及受力狀況，為管道的安全使用和報廢提供參考。

參考文獻(xiàn)

[1]余偉煒，高炳軍.ANSYS在機械與化工裝備中的應(yīng)用[M].北京：中國水利水電出版社，2007.

[2]隋之鋒，郝點，陳海峰.ANSYS在兩端固定厚壁管道應(yīng)力分析中的應(yīng)用[J].貴州化工,2008（12）：34-37.

[3]劉純，胡波濤.應(yīng)用ANSYS計算汽水管道應(yīng)力[J].鍋爐技術(shù),2005（4）：31-33

[4]周昌玉，賀小華.有限元分析的基本方法及工程應(yīng)用[M].北京：化學(xué)工業(yè)出版社，2006.

Stress Analysis of Pipeline Deformation Based on ANSYS

Zhao Fuxing Zhou Pengfei Qian Yinghao

Abstract：On account of the buried pipeline deformation due to the external load,the stress analysis function of the software ANSYS is applied so that the stress distribution of the pipeline after the deformation is obtained by exerting displacement load,which then becomes the reference basis for the pipeline safety evaluation.

Key words：Pipeline; Stress analysis; Deformation; Displacement; Finite element

收稿日期：（2015-07-11）

中圖分類號TQ 050.2

DOI：10.16759/j.cnki.issn.1007-7251.2016.04.012