蔡海潮 尚振東 晁紅軍

1.河南科技大學(xué)機(jī)電工程學(xué)院，洛陽，471003

2.機(jī)械裝備先進(jìn)制造河南省協(xié)同創(chuàng)新中心，洛陽，471003

3.國(guó)家電網(wǎng)河南電力公司，洛陽，471000

0 引言

作為五大運(yùn)輸工具之一，管道在石油石化等行業(yè)中得到廣泛使用，其目的是輸送石油、天然氣等能源物質(zhì)。管道在長(zhǎng)期使用過程中，管內(nèi)運(yùn)輸介質(zhì)的腐蝕性會(huì)導(dǎo)致管內(nèi)壁減薄現(xiàn)象的發(fā)生，與此同時(shí)，埋置于地下或暴露在空氣中的管外壁又極易受到土壤或潮濕空氣的腐蝕影響，而且極易受到外界損傷，繼而造成管路泄漏和爆管事故。為保護(hù)管道免受腐蝕或外力損傷，往往在其內(nèi)外表面附著一層包覆層，但包覆層的存在又使得常用的超聲、渦流、磁粉及目視等無損檢測(cè)技術(shù)無法對(duì)管道中的缺陷逐點(diǎn)進(jìn)行檢測(cè)［1］。超聲導(dǎo)波檢測(cè)技術(shù)是一種新興的無損檢測(cè)方法［2］，特別適合對(duì)長(zhǎng)輸管道進(jìn)行快速大面積檢測(cè)，檢測(cè)時(shí)只需剝開管道某一位置的包覆層，并在該位置放置超聲換能器激勵(lì)出所需的超聲導(dǎo)波，則導(dǎo)波沿著管道可傳播較遠(yuǎn)的距離，此時(shí)在距接收端較遠(yuǎn)處設(shè)置換能器來接收傳播的導(dǎo)波信號(hào)，即可獲得關(guān)于整個(gè)管道的狀態(tài)信息［3］。

包覆層的存在會(huì)使管道中傳播的超聲導(dǎo)波受到影響。KWUN等［4］通過實(shí)驗(yàn)初步研究了在有煤焦油磁漆防腐層的地埋管道中扭轉(zhuǎn)模態(tài)的衰減特性，研究結(jié)果表明，扭轉(zhuǎn)模態(tài)的衰減會(huì)隨著激勵(lì)頻率的增加而增大。ALLEYNE等［5］利用低頻L（0，2）模態(tài)對(duì)地埋管道中的腐蝕缺陷進(jìn)行了檢測(cè)，結(jié)果表明，超聲導(dǎo)波在外表面包覆一層黏性介質(zhì)的管道中傳播時(shí)會(huì)產(chǎn)生一定程度的衰減，其衰減與介質(zhì)的厚度以及超聲導(dǎo)波的波結(jié)構(gòu)和頻率有關(guān)。SHIN等［6］采用超聲導(dǎo)波縱向模態(tài)在帶聚乙烯包覆層管道中進(jìn)行了試驗(yàn)，通過可變角斜探頭來調(diào)節(jié)相速度及激勵(lì)頻率，產(chǎn)生最佳模態(tài)以檢測(cè)該類管道中的缺陷。BARSHINGER 等［7?9］對(duì)帶黏彈性包覆層管道中縱向模態(tài)的傳播特性進(jìn)行了研究，并在帶煤焦油黏彈性包覆層管道中進(jìn)行了檢測(cè)試驗(yàn)，研究結(jié)果表明，頻率為190 kHz的L（0，2）模態(tài)和頻率為630 kHz的L（0，3）模態(tài)在包覆層管道中傳播時(shí)衰減較小，而頻率為740 kHz的L（0，3）模態(tài)則由于衰減較大，無法用于缺陷檢測(cè)，但沒有對(duì)適合管道檢測(cè)的低頻L（0，2）模態(tài)的傳播特性作進(jìn)一步研究。HUA等［10］采用ABAQUS軟件建立了不同包覆層材料的有限元模型，結(jié)果表明，包覆層的厚度對(duì)導(dǎo)波傳播有較大的影響。KIRBY等［11］研究了縱向模態(tài)在帶包覆層管道中的傳播，結(jié)果表明，包覆層材料的切變性對(duì)縱向模態(tài)的傳播有較大的影響，但沒有進(jìn)一步研究包覆層材料密度對(duì)導(dǎo)波傳播的影響。

本文根據(jù)理論分析結(jié)果，采用數(shù)值模擬和實(shí)驗(yàn)研究的手段分析了瀝青包覆層性能對(duì)導(dǎo)波縱向模態(tài)傳播的影響，然后通過實(shí)驗(yàn)在帶黏彈性包覆層管道中激勵(lì)和接收縱向模態(tài)，并將實(shí)驗(yàn)結(jié)果與理論進(jìn)行了比較和分析。

1 帶包覆層管道超聲導(dǎo)波傳播模型

依據(jù)空心管道的導(dǎo)波傳播模型，可建立帶包覆層管道的超聲導(dǎo)波傳播模型。帶包覆層管道的鋼管外層和包覆層滿足鋼及包覆層的控制方程。研究帶包覆層管道的超聲導(dǎo)波傳播特性，需要考慮管道中的導(dǎo)波傳播與包覆層之間的相互作用，因此，建立帶包覆層管道的頻散方程要同時(shí)考慮管道內(nèi)外表面的邊界條件，管道外表面和包覆層內(nèi)表面之間的徑向位移和徑向應(yīng)力均是連續(xù)的。對(duì)于帶包覆層管道，包覆層外表面的應(yīng)力為零，它與管徑的外表面作用時(shí)應(yīng)滿足位移和應(yīng)力連續(xù)的邊界條件；而管道內(nèi)表面處界面為自由表面，相應(yīng)的應(yīng)力分量為零。聯(lián)立管道內(nèi)外表面的邊界條件可以建立一組特征方程，該方程即為帶包覆層管道的頻散方程，由于該方程不存在解析解的形式，故只能通過數(shù)值方法求解，求解后可獲得帶包覆層管道的頻散關(guān)系。對(duì)帶包覆層管道進(jìn)行數(shù)值求解所得的根波數(shù)為復(fù)數(shù)，波數(shù)復(fù)數(shù)根的實(shí)部為相速度頻散曲線，而波數(shù)復(fù)數(shù)根的虛部為衰減頻散曲線。

計(jì)算所用的管道密度為7 850 kg m3，縱波波速為5 900 m s，橫波波速為5 900 m s，鋼管外徑為273 mm，壁厚為7 mm；包覆層為黏彈性材料瀝青，其密度為1 500 kg m3，縱波波速為1 860 m s，橫波波速為750 m s，縱波衰減系數(shù)為0.023。

圖1 含包覆層管道頻散曲線Fig.1 dispersion curve of coating pipe

圖1表明，在0～400 kHz范圍內(nèi)存在3個(gè)縱向模態(tài)，而且在低頻段，相速度和群速度頻散曲線與空管中的基本相同。但隨著頻率的增加，各模態(tài)的頻散特性與空管時(shí)有較大區(qū)別。圖1b中，L（0,2）模態(tài)在30～200 kHz范圍內(nèi)頻散較小，而且由于包覆層的存在，各模態(tài)導(dǎo)波傳播時(shí)會(huì)引起一定程度的衰減；圖1c中，L（0,2）模態(tài)在低頻段范圍內(nèi)衰減較小，說明受包覆層的影響較小，因而在對(duì)管道缺陷進(jìn)行檢測(cè)時(shí)，為使管道中傳播的能量泄漏較少，以保證檢測(cè)距離和缺陷靈敏度，可以選擇較低頻段的L（0,2）模態(tài)進(jìn)行檢測(cè)。

2 帶包覆層管道導(dǎo)波傳播有限元模擬

管道外有包覆層存在時(shí)，管外壁的受激波動(dòng)問題可歸結(jié)為管壁彈性體和包覆層之間的非定常耦聯(lián)波動(dòng)問題，管外壁的變形能與包覆層變形能之間將發(fā)生相互作用，因此，模擬研究帶包覆層管道中導(dǎo)波傳播特性時(shí)需要考慮導(dǎo)波對(duì)包覆層的作用以及包覆層對(duì)導(dǎo)波傳播的影響。求解帶包覆層管道耦合波動(dòng)問題時(shí)采用有限元模擬的方法［12］。

采用有限元模擬，管道材料為X70管道鋼，材料的力學(xué)參數(shù)屬性見表1，采用Comsol有限元軟件模擬時(shí)，管道為彈性體，其單元模型物理場(chǎng)設(shè)定為固體力學(xué)模塊。

表1 X70鋼的材料屬性Tab.1 Material properties of X70 steel

管道外包覆層一般為黏彈性材料，在有限元模擬時(shí)假設(shè)該部分為彈性體，其單元模型也用固體力學(xué)模塊來表示，但需要引入Kelvin-Voigt模型來描述包覆層引起的導(dǎo)波能量衰減。有限元模型中包覆層的材料屬性見表2。

表2 包覆層的材料屬性Tab.2 Material properties of coating

包覆層管道有限元模擬時(shí)，管道兩端面邊界條件設(shè)定為自由邊界，包覆層厚度設(shè)置為2 mm，包覆層的外壁邊界條件設(shè)定為自由邊界。根據(jù)管道中導(dǎo)波波長(zhǎng)和包覆層中波長(zhǎng)［13］，管道部分網(wǎng)格單元長(zhǎng)度設(shè)為8 mm，包覆層部分網(wǎng)格單元長(zhǎng)度設(shè)為2 mm。根據(jù)管道幾何尺寸和單元長(zhǎng)度，計(jì)算獲得有限元模擬時(shí)的時(shí)間參數(shù)，綜合分析最終設(shè)定時(shí)間步長(zhǎng)為0.3μs，總計(jì)算時(shí)間為1.1 ms。

有限元模擬中沿管道軸向施加位移載荷作為激勵(lì)，模擬產(chǎn)生縱向模態(tài)超聲導(dǎo)波。有限元模擬完成后，分別提取管道激勵(lì)端和另一端8個(gè)點(diǎn)的時(shí)間位移信號(hào)并做數(shù)據(jù)平均處理，獲得導(dǎo)波在管道中傳播的反射回波信號(hào)和透射信號(hào)。

3 包覆層對(duì)超聲導(dǎo)波傳播特性的影響

本文在帶包覆層管道三維有限元模型的基礎(chǔ)上，研究了瀝青包覆層性能對(duì)管道中超聲導(dǎo)波傳播衰減特性的影響，管道上瀝青包覆層的覆蓋長(zhǎng)度為0.6 m，選擇5周期、中心頻率分別為40 kHz、70 kHz、100 kHz、160 kHz的脈沖信號(hào)作為激勵(lì)在管道中激勵(lì)縱向模態(tài)超聲導(dǎo)波。分別研究激勵(lì)頻率、包覆層密度和包覆層縱波衰減系數(shù)對(duì)導(dǎo)波衰減特性的影響。不同激勵(lì)頻率下的導(dǎo)波衰減結(jié)果見圖2。結(jié)果表明，隨著激勵(lì)頻率的增大，管道中的超聲導(dǎo)波衰減逐漸加大，這與理論分析結(jié)果一致，說明激勵(lì)頻率越大，導(dǎo)波在管道中傳播時(shí)泄漏進(jìn)包覆層的能量越多。圖2還表明，在同一激勵(lì)頻率下，隨著包覆層長(zhǎng)度的增大，導(dǎo)波的衰減也會(huì)進(jìn)一步增大，說明超聲導(dǎo)波在包覆層管道中傳播時(shí)，其檢測(cè)能力會(huì)大幅下降。

圖2不同頻率下的導(dǎo)波衰減Fig.2 Guided wave attenuation of different frequencies

為研究包覆層密度對(duì)導(dǎo)波衰減特性的影響，利用建立的帶包覆層管道三維有限元模型，采用不同的激勵(lì)頻率分別研究不同密度的瀝青包覆層對(duì)管中超聲導(dǎo)波傳播特性的影響。不同激勵(lì)頻率下包覆層密度改變時(shí)的導(dǎo)波衰減結(jié)果見圖3。

圖3 不同包覆層密度下的導(dǎo)波衰減Fig.3 Guided wave attenuation of different coating density

圖3結(jié)果表明，同一密度下隨激勵(lì)頻率的增大，導(dǎo)波衰減增大。同一激勵(lì)頻率下隨包覆層密度的逐漸增大，導(dǎo)波的衰減也逐漸增大，而且呈線性變化趨勢(shì)。從不同頻率的衰減曲線還可看出，不同激勵(lì)頻率下的導(dǎo)波衰減特性曲線斜率不同，頻率越大，斜率越大，說明在激勵(lì)頻率較大時(shí)，導(dǎo)波能量的衰減對(duì)包覆層的密度更敏感，而頻率越低，斜率越小，導(dǎo)波能量的衰減對(duì)包覆層密度的變化不敏感。

圖4 不同縱波衰減系數(shù)下的導(dǎo)波衰減Fig.4 Guided wave attenuation of different longitudinal wave attenuation coefficient

本文在有限元模型的基礎(chǔ)上，采用不同的激勵(lì)頻率分別研究了不同縱波衰減系數(shù)的瀝青包覆層對(duì)管中超聲導(dǎo)波傳播特性的影響。不同激勵(lì)頻率下包覆層縱波衰減系數(shù)改變時(shí)的導(dǎo)波衰減結(jié)果見圖4。圖4結(jié)果表明，同一縱波衰減系數(shù)下隨激勵(lì)頻率的增大，導(dǎo)波衰減逐漸增大。當(dāng)激勵(lì)頻率增大到160 kHz時(shí)，導(dǎo)波衰減明顯增大。而在同一激勵(lì)頻率下，隨包覆層縱波衰減系數(shù)的逐漸增大，導(dǎo)波的衰減呈線性變化趨勢(shì)，略微增大。比較不同頻率的衰減曲線可看出，任一激勵(lì)頻率下的導(dǎo)波衰減特性曲線斜率均較小，表明導(dǎo)波能量的衰減對(duì)包覆層縱波衰減系數(shù)的變化敏感程度較低。

有限元分析結(jié)果表明，激勵(lì)頻率對(duì)帶包覆層管道的檢測(cè)影響較大，在用低頻縱向模態(tài)超聲導(dǎo)波對(duì)帶低密度包覆層管道檢測(cè)時(shí)，其靈敏度較高，檢測(cè)效果較好。

4 實(shí)驗(yàn)驗(yàn)證

為了驗(yàn)證有限元模擬的結(jié)果，本文對(duì)帶包覆層管道進(jìn)行了實(shí)驗(yàn)。實(shí)驗(yàn)用管道幾何尺寸與有限元仿真管道幾何尺寸相同。實(shí)驗(yàn)用檢測(cè)裝置和換能器布置方式見圖5。檢測(cè)系統(tǒng)由壓電換能器、信號(hào)激勵(lì)接收裝置和計(jì)算機(jī)組成，壓電換能器采用長(zhǎng)度伸縮型壓電陶瓷片，其尺寸為20mm×3 mm×0.5 mm，信號(hào)激勵(lì)接收裝置為美國(guó)物理聲學(xué)公司的多通道激勵(lì)接收卡，為使得激勵(lì)的信號(hào)在傳播過程中頻散較小，實(shí)驗(yàn)選取激勵(lì)產(chǎn)生單音頻信號(hào)用來激勵(lì)超聲換能器，激勵(lì)信號(hào)在管道中傳播并被另一端的接收換能器接收，接收到的導(dǎo)波信號(hào)經(jīng)放大后送入計(jì)算機(jī)進(jìn)行分析處理。

圖5 實(shí)驗(yàn)裝置示意圖Fig.5 Schematic diagram of experiment device

為開展包覆層對(duì)導(dǎo)波傳播特性的影響實(shí)驗(yàn)，在SBS防水層上刷萬能膠并粘貼在管道表面用來模擬帶包覆層管道工況，SBS防水層材料由瀝青制成，其材料屬性由實(shí)驗(yàn)測(cè)得，參數(shù)與表2參數(shù)基本相同。

實(shí)驗(yàn)中，導(dǎo)波信號(hào)激勵(lì)頻率為40～160 kHz，分別獲得激勵(lì)頻率為 40k Hz、70 kHz、100 kHz、160 kHz時(shí)的導(dǎo)波L（0,2）模態(tài)接收信號(hào)。有無包覆層時(shí)的端面回波波形見圖6。

圖6 激勵(lì)頻率為100 kHz時(shí)端面回波波形Fig.6 Guided wave echo waveform at 100 kHz excitation frequency

對(duì)圖6中波形進(jìn)行波速分析可知，激勵(lì)的模態(tài)為L(zhǎng)（0,2）模態(tài)，該模態(tài)頻散小，而且在有無包覆層管道中均可以接收到兩個(gè)端面回波，但有包覆層時(shí)端面回波信號(hào)幅值明顯減小。按照衰減特性計(jì)算公式αi=-20lg（AiA0）（Ai表示第i個(gè)波包的反射幅值，A0表示第一個(gè)波包的反射幅值），獲得不同頻率下包覆層對(duì)導(dǎo)波傳播的衰減值，不同頻率下導(dǎo)波傳播的衰減結(jié)果見圖7。

圖7 不同頻率下導(dǎo)波衰減結(jié)果Fig.7 Guided wave attenuation results of different frequencies

由圖7可看出，帶包覆層管道中導(dǎo)波的衰減隨頻率的增大而增大，而且隨著包覆層長(zhǎng)度的增大，導(dǎo)波傳播能量逐漸減小，實(shí)驗(yàn)結(jié)果與仿真結(jié)果的衰減規(guī)律基本一致，證明了有限元模擬結(jié)果的正確性。

為觀察包覆層密度對(duì)導(dǎo)波傳播特性的影響，實(shí)驗(yàn)過程中在管道表面粘貼不同密度的包覆層進(jìn)行導(dǎo)波傳播檢測(cè)實(shí)驗(yàn)，分別獲得不同激勵(lì)頻率時(shí)不同包覆層密度的L（0,2）模態(tài)衰減情況，見表3。

表3 不同密度包覆層下L（0,2）模態(tài)的衰減情況Tab.3 Attenuation of L（0,2）mode under different coating density dB

表3表明，在同一頻率下，不同包覆層密度其衰減值不同，密度較小的包覆層其衰減較小，且隨著激勵(lì)頻率的增大其衰減變化趨勢(shì)與數(shù)值模擬結(jié)果一致。

5 結(jié)論

采用L（0,2）模態(tài)導(dǎo)波對(duì)帶瀝青包覆層管道進(jìn)行檢測(cè)是有效的，在激勵(lì)頻率較大的情況下，包覆層密度和縱波衰減系數(shù)會(huì)對(duì)導(dǎo)波能量衰減產(chǎn)生較大的影響，研究結(jié)果為工程管道的大范圍在役檢測(cè)提供了理論指導(dǎo)。

［1］ ROSE J L.GUIDED WAVES Successes and Chal?lenges for Ultrasonic Testing in NDT and SHM［J］.Materials Evaluation,2010，68（5）：495?500.

［2］ 何存富,劉增華,孫雅欣,等.基于超聲導(dǎo)波技術(shù)對(duì)彎管中缺陷檢測(cè)的實(shí)驗(yàn)研究［J］.中國(guó)機(jī)械工程,2005,16（18）：1662?1665.HE Cunfu，LIU Zenghua，SUN Yaxin，et al.Experi?mental Study on Defect Detection in Curved Pipes Using Ultrasonic Guided Waves Technique［J］.Chi?na Mechanical Engineering，2005,16（18）：1662?1665.

［3］ CARANDENTE R，CAWLEY P.A Method to Es?timate the Size of Corrosion Patches with Guided Waves in Pipes［J］.Review of Progress in Quantita?tive Nondestructive Evaluation,2012，1430（1）：1929?1936.

［4］ KWUN H，KIM S Y，CHOI M S，et al.Torsional Guided?wave Attenuation in Coal?tar?enamel?coat?ed，Buried Piping［J］.Ndt&E International,2004，37（8）：663?665.

［5］ ALLEYNE D N，PAVLAKOVIC B，LOWE M J S，et al.Rapid Long?range Inspection of Chemical Plant Pipework Using Guided Waves［J］.Key Engineering Materials,2001,270/273（1）：434?441.

［6］ SHIN H J，SONG S J，PARK J S.Long Range In?spection of Polyethylene Coated Steel Pipes by Us?ing Guided Waves［C］//Review of Progress in Quantitative Nondestructive Evaluation.Brunswick，2001:1194?1201.

［7］ BARSHINGER J，ROSE J L，AVIOLI M J.Guided Wave Resonance Tuning for Pipe Inspection［J］.Journal of Pressure Vessel Technology?Transac?tions of the ASME,2002，124（3）：303?310.

［8］ BARSHINGER J N，ROSE J L.Ultrasonic Guided Wave Propagation in Pipes with Viscoelastic Coat?ings［M］.Review of Progress in Quantitative Nonde?structive Evaluation,2002,615:239?246.

［9］ BARSHINGER J N，ROSE J L.Guided Wave Prop?agation in an Elastic Hollow Cylinder Coated with a Viscoelastic Material［J］.IEEE Transactions on Ul?trasonics Ferroelectrics and Frequency Control,2004，51（11）：1547?1556.

［10］ HUA J，ROSE J L.Guided Wave Inspection Pene?tration Power in Viscoelastic Coated Pipes［J］.Or Insight,2010，52（4）：195?205.

［11］ KIRBY R，ZLATEV Z，MUDGE P.On the Scatter?ing of Longitudinal Elastic Waves from Axisym?metric Defects in Coated Pipes［J］.Journal of Sound and Vibration,2013，332（20）：5040?5058.

［12］ 沈立華，王悅民，鄭明雪，等.基于磁致伸縮效應(yīng)導(dǎo)波的數(shù)值模擬和實(shí)驗(yàn)研究［J］.中國(guó)機(jī)械工程，2007,18（5）：532?535.SHEN Lihua，WANG Yuemin，ZHENG Mingxue，et al.Numerical Simulation and Experimental Re?search on Guided Waves Based on Mangentostric?tiveEffect［J］.ChinaMechanicalEngineering，2007,18（5）：532?535.

［13］ 何存富，李偉，吳斌.扭轉(zhuǎn)模態(tài)導(dǎo)波檢測(cè)管道縱向缺陷的數(shù)值模擬［J］.北京工業(yè)大學(xué)學(xué)報(bào)，2007,33（10）:1009?1013.HE Cunfu，LI Wei，WU Bin.Numerical Simulation of Longitudinal Crack Detection in Pipes Using Ul?trasonic Guided wave T（0,1）Mode［J］.Journal of Beijing University of Technology，2007，33（10）:1009?1013.*