張 行 劉 彤 劉思敏 許學(xué)峰 高孟琦 閆增瑞 張仕民

(1.中國石油大學(xué)(北京)機(jī)械與儲運工程學(xué)院 2.中油國際管道公司)

0 引 言

油氣管道投產(chǎn)前，需要對新建管道內(nèi)部進(jìn)行檢查評估[1-2]。在役管道也經(jīng)常使用測徑清管器檢測其內(nèi)部的變形情況[3-5]。目前，檢測管道內(nèi)部變形情況的常用工具是管道測徑清管器。

傳統(tǒng)測徑清管器是在常規(guī)清管器上加裝測徑鋁盤[6],當(dāng)測徑清管器運行結(jié)束后，通過檢測鋁盤變形量粗略估計管道內(nèi)部發(fā)生的最大變形量。這種測徑清管器檢測精度低，不能對管道多個變形位置進(jìn)行測量，并且還存在鋁制測徑盤不可重復(fù)利用的缺點。近年來，管道內(nèi)部變形檢測方法主要有接觸式通徑檢測和非接觸式無損檢測[7-15]。王宇楠等[16]在泡沫清管器中放置磁塊，根據(jù)磁場變化對管內(nèi)殘余物質(zhì)的高度和外形等特征進(jìn)行預(yù)先測量。馬書義等[17]根據(jù)回波信號的時間歷程、幅值等信息獲取管道截面的變化。趙翰學(xué)等[18]研發(fā)了150 mm大量程、高精度渦流傳感器實現(xiàn)對管道變形的檢測。WANG H.昊等[19]為檢測直徑小于200 mm的管道而提出鋼板彈簧卡尺探測臂。接觸式通徑檢測方法不僅存在探測臂引發(fā)機(jī)械故障的風(fēng)險，而且還有檢測速度慢的缺點。非接觸式無損檢測方法檢測費用較高，設(shè)備結(jié)構(gòu)復(fù)雜，再加上管道結(jié)構(gòu)多樣以及服役環(huán)境等原因，應(yīng)用存在局限性[20-21]。

針對上述問題，本文提出一種油氣管道柔性測徑清管器。該清管器將傳統(tǒng)鋁制測徑盤改為柔性材質(zhì)，并將應(yīng)變片置于超彈性柔性盤中,其工作原理是：當(dāng)測徑清管器越過管道凹陷處時，柔性盤發(fā)生變形，應(yīng)變片也隨之發(fā)生變形，應(yīng)變儀測得測徑盤遇阻時的應(yīng)變變化，據(jù)此可以快速識別管道內(nèi)通徑變化情況，對管道變形處實現(xiàn)較準(zhǔn)確的檢測和量化。這種檢測手段既可以對管內(nèi)多個變形進(jìn)行測量，同時又能減少引發(fā)機(jī)械故障的風(fēng)險。其檢測設(shè)備在運行過程中對管壁無損傷且柔性測徑盤可重復(fù)利用。

本文以聚氨酯材質(zhì)測徑盤為研究對象，采用ABAQUS 2016有限元仿真軟件對其在凹陷管道內(nèi)運行情況進(jìn)行模擬。從放置面、縱向分布和橫向分布3方面確定應(yīng)變片放置位置，討論厚度、瓣數(shù)及聚氨酯材料硬度3個參數(shù)變化時，測徑盤應(yīng)力與應(yīng)變的分布，并進(jìn)行分析。

1 有限元模型

1.1 柔性測徑盤模型

柔性測徑盤如圖1所示。柔性測徑盤外徑Dg為管道內(nèi)徑的95%，Dg=196.64 mm；Dr表示剛體外徑，即測徑盤的內(nèi)徑,Dr=51 mm；D表示每瓣之間的距離，D=10 mm；tg為測徑盤的厚度。

圖1 柔性測徑盤Fig.1 Flexible gauging plate

實際測徑盤是用 “U”形將盤分割成若干瓣，n表示測徑盤的瓣數(shù)。討論測徑盤的厚度(tg)、瓣數(shù)(n)和聚氨酯硬度等3個參數(shù)對柔性測徑盤應(yīng)力和應(yīng)變的影響，3個參數(shù)的變化如表1所示。

表1 不同模擬情況下測徑盤參數(shù)Table 1 Gauging plate parameters under different simulation conditions

1.2 管道及凹陷模型

凹陷大小通常由深度、長度和寬度決定[22]。本文研究的凹陷如圖1b所示。其中H是凹陷深度，Dp是管道外徑，t是管道壁厚。以?203.2 mm(8 in)管道為例，Rp=109.5 mm，t=6 mm。

對于管道凹陷基于深度的定義，許多規(guī)范都提出了評價標(biāo)準(zhǔn)[23-25]，本文中凹陷深度為15 mm。由于管道和清管器都是回轉(zhuǎn)體，所以使用極坐標(biāo)系進(jìn)行分析。從管道軸向看，逆時針方向是正方向，凹陷中心位于270°。

1.3 測徑盤和帶有凹陷管道的有限元模型

在ABAQUS 2016有限元軟件中，測徑盤在凹陷管道運行時的有限元模型如圖2所示。為保證網(wǎng)格質(zhì)量，測徑盤選擇六面體單元網(wǎng)格，在與管道凹陷接觸的測徑盤邊緣進(jìn)行倒角處理并加密網(wǎng)格，共有68 248個網(wǎng)格單元。管道同樣選擇六面體單元，共有10 318個網(wǎng)格單元。在管道凹陷處，選擇四面體進(jìn)行網(wǎng)格劃分，共有5 115個網(wǎng)格單元。

圖2 測徑盤有限元模型Fig.2 Finite element model of gauging plate

管道的彈性模量是2.01 MPa，泊松比是0.3。測徑盤采用聚氨酯材料，以Mooney-Rivlin為本構(gòu)模型[26]，其邵氏硬度與材料特性常數(shù)C01、C10的關(guān)系如表2所示?？紤]管道無軸向位移，管道外壁及凹陷處受完全約束。為提高運算速率，將管道凹陷處簡化為剛體。

表2 聚氨酯材料參數(shù)Table 2 Polyurethane material parameters

2 測徑盤過管道凹陷運行分析

2.1 測徑盤過管道凹陷過程的應(yīng)力分析

第10組算例用于分析測徑盤過管道凹陷過程中應(yīng)力與應(yīng)變的變化情況。測徑盤在凹陷管道運行過程如圖3所示。4個典型位置對應(yīng)的測徑盤應(yīng)力云圖如圖4所示。以測徑盤與凹陷接觸的面為正面，反之為背面。

圖3 模擬過程及結(jié)果Fig.3 Simulation process and results

測徑盤接觸凹陷前無應(yīng)力變化。當(dāng)測徑盤繼續(xù)向前運行，底瓣接觸管道凹陷并緊貼凹陷輪廓向上“抬起”，應(yīng)力隨之增大，到達(dá)凹陷最高處時，測徑盤底瓣徑向上升到最大高度，測徑盤應(yīng)力達(dá)到最大。然后，測徑盤底瓣沿管道凹陷輪廓慢慢“落下”，應(yīng)力隨之減小，越過凹陷后測徑盤恢復(fù)原狀，底瓣所受應(yīng)力也恢復(fù)至0。

由圖4可知，測徑盤遇凹陷過程中，離測徑盤孔徑近的位置應(yīng)力更大，其應(yīng)力沿圓盤半徑至邊緣方向以“水波紋”的形式逐級遞減。在測徑盤底瓣“頸部”邊緣會出現(xiàn)應(yīng)力集中，若將應(yīng)變片放置在此處，則測徑盤經(jīng)過微小凹陷時也會產(chǎn)生較大應(yīng)力，這不能反映管道內(nèi)部真實狀況?？紤]到實際操作，在邊緣處不便加裝應(yīng)變片，所以在后續(xù)分析應(yīng)變片優(yōu)選位置時，此處不做考慮。

圖4 4個典型位置的應(yīng)力云圖Fig.4 Cloud chart of stress at 4 typical positions

2.2 確定應(yīng)變片放置面

在測徑盤遇凹陷過程中，其正面受到拉應(yīng)力，背面受到壓應(yīng)力。放置應(yīng)變片時，應(yīng)以應(yīng)變變化更大的一面，即應(yīng)變率更高的面為優(yōu)選。從測徑盤底瓣外緣沿半徑方向的中心線上等距選取10 個節(jié)點，如圖5所示。以Point2為例，比較正、背面的應(yīng)變率，其結(jié)果如圖6所示，正面的應(yīng)變變化比背面更明顯，以此可知應(yīng)在測徑盤正面放置應(yīng)變片。

圖5 測徑盤正面節(jié)點取情況Fig.5 Selection of front nodes and lines of gauging plate

圖6 Point2正、背面應(yīng)變和應(yīng)變率隨時間變化曲線Fig.6 Variation curve of front and back strain and strain rate with time at Point 2

2.3 確定應(yīng)變片縱向位置

沿圓盤半徑方向做7條線將底瓣等分，每條線的起點都為測徑盤孔徑邊緣，等距選取10個節(jié)點，如圖7所示。測徑盤每條線上不同點的應(yīng)變-位移曲線如圖8所示。首先分析第4條線上10個點的應(yīng)變情況。由圖8a可知，Point1～Point3在10個點中應(yīng)變變化顯著，其中Point2在測徑盤底端到達(dá)管道凹陷最高處時,應(yīng)變達(dá)到最大。選取Point1、Point3、Point5、Point7、Point9進(jìn)行應(yīng)變比較。從圖8b～8h可以看出：由于Point3接近應(yīng)力集中點，所以其應(yīng)變變化非常明顯；圖8c～圖8g中Point1與Point3這兩點之間應(yīng)變相差微小，后幾點之間應(yīng)變差距較大。由此可得，應(yīng)變片縱向位置介于Point1～Point3之間最合適。

圖7 測徑盤正面7條線Fig.7 7 lines on the front of the gauging plate

圖8 測徑盤每條線上不同點的應(yīng)變-位移曲線Fig.8 Strain-displacement curve of different points on each line of the gauging plate

2.4 確定應(yīng)變片橫向位置

Line1～Line7的前3個節(jié)點應(yīng)變-位移曲線如圖9所示。

圖9 不同線上的前3節(jié)點應(yīng)變-位移曲線Fig.9 Strain-displacement curve of the first three nodes on different lines

由圖9可知：第3條線上第1節(jié)點的應(yīng)變隨位移變化最明顯，第1條線變化最弱；第1條和第7條線的應(yīng)變曲線幾乎重合且變化最不明顯；由于第1和第7條線的第3節(jié)點距離測徑盤應(yīng)力集中點較近，所以不能合理有效地提供應(yīng)變數(shù)據(jù)。由此可得，應(yīng)變片的橫向位置應(yīng)處在Line2～Line6之間。綜上所述，應(yīng)變片放置在測徑盤正面每瓣根部中心位置，有利于提升檢測系統(tǒng)的靈敏度。

2.5 網(wǎng)格獨立性

將算例10用于網(wǎng)格獨立性分析。當(dāng)單元尺寸分別為1.1和1.5 mm時，計算得測徑盤底瓣位于C點時的應(yīng)力分別為2.301和2.299 MPa，二者相差0.09%，這兩種網(wǎng)格具有良好一致性。因此，構(gòu)建的網(wǎng)格具有獨立性。

3 不同參數(shù)對測徑盤過凹陷的應(yīng)力與應(yīng)變影響

由前文可以得出，Line4 的第2節(jié)點應(yīng)變變化最顯著，將此點記為M，M點是對不同參數(shù)研究最具代表性的一點。

3.1 厚度對測徑清管器測徑盤M點應(yīng)力和應(yīng)變的影響

設(shè)n=8，聚氨酯材料硬度為70 HS時，測徑盤厚度tg=12、15、18和21 mm，其他條件不變。圖10為不同厚度下測徑盤的應(yīng)力與應(yīng)變分布。由圖10可知，M點在測徑盤底瓣底端位于凹陷最高點處的應(yīng)力應(yīng)變分布，相應(yīng)的應(yīng)變云圖如圖11所示。從圖10可得，厚度從12 mm增加到21 mm，M點應(yīng)力從0.67 MPa增加到1.42 MPa，其應(yīng)變變化也有相似趨勢。除M點外，測徑盤其他位置的應(yīng)力應(yīng)變也隨著測徑盤厚度的增加而增加。

圖10 不同厚度下測徑盤的應(yīng)力與應(yīng)變分布Fig.10 Stress-strain distribution of gauging plate with different thickness

圖11 不同厚度測徑盤的應(yīng)變云圖Fig.11 Cloud chart of strain of different thickness of gauging plate

3.2 瓣數(shù)對測徑盤M點應(yīng)力和應(yīng)變的影響

當(dāng)tg=20 mm、聚氨酯材料硬度為70 HS時，測徑盤瓣數(shù)從6瓣變?yōu)?、8和9瓣，其他條件不變。圖12為M點在測徑盤底瓣底端位于凹陷最高點處的應(yīng)力與應(yīng)變分布，相應(yīng)的應(yīng)變云圖如圖13所示。從圖12可知，隨著瓣數(shù)的增加，應(yīng)力和應(yīng)變變化并不顯著。

圖12 不同瓣數(shù)下測徑盤的應(yīng)力與應(yīng)變分布Fig.12 Stress-strain distribution of gauging plate with different flap numbers

圖13 不同瓣數(shù)下測徑盤的應(yīng)變云圖Fig.13 Cloud chart of strain of gauging plate with different flap numbers

3.3 聚氨酯材料硬度對測徑盤M點應(yīng)力和應(yīng)變的影響

當(dāng)tg=20 mm、n=8時，測徑盤材料硬度從60 HS依次變?yōu)?5、70和75 HS，其他條件不變。圖14為M點在測徑盤底瓣底端位于凹陷最高點時的應(yīng)力與應(yīng)變分布，相應(yīng)的應(yīng)變云圖如圖15所示。

圖14 不同聚氨酯材料硬度下測徑盤的應(yīng)力與應(yīng)變分布Fig.14 Stress-strain distribution of gauging plate with different hardness

圖15 不同聚氨酯材料硬度測徑盤的應(yīng)變云圖Fig.15 Cloud chart of strain of gauging plate with different polyurethane hardness

從圖14可知,M點的應(yīng)力隨聚氨酯材料硬度的增加而增加，但應(yīng)變沒有實質(zhì)性的變化。這是由于隨聚氨酯材料硬度的增加，測徑盤變得不易彎曲，所以應(yīng)變變化很小。

4 結(jié) 論

(1)柔性測徑清管器既可實現(xiàn)對管內(nèi)多個變形進(jìn)行準(zhǔn)確檢測，其柔性測徑盤又可重復(fù)利用，同時可以提升檢測作業(yè)的安全性，其成本低、適用范圍廣的特點可滿足不同的作業(yè)環(huán)境要求。

(2)在柔性測徑盤越過凹陷過程中，其變形瓣的應(yīng)力沿半徑方向呈“水波紋”形式逐級遞減分布，在變形瓣“頸部”邊緣會出現(xiàn)應(yīng)力集中現(xiàn)象。

(3)應(yīng)變片放置在測徑盤正面每瓣根部中心位置時應(yīng)變變化更顯著，有利于提升檢測系統(tǒng)的靈敏度。

(4)測徑盤的厚度比瓣數(shù)和材料硬度對其應(yīng)力與應(yīng)變的影響大，瓣數(shù)對應(yīng)力和應(yīng)變的影響最小。