薄壁小直徑管道導(dǎo)波探傷的傳感器布置方案優(yōu)化

王偉，游鵬輝，鐘萬里，謝濤，徐俊，鐘力強(qiáng)，肖曉暉(.廣東電網(wǎng)公司電力科學(xué)研究院，廣東 廣州，50080；.武漢大學(xué) 動(dòng)力與機(jī)械學(xué)院，湖北 武漢，43007)

薄壁小直徑管道導(dǎo)波探傷的傳感器布置方案優(yōu)化

王偉1，游鵬輝2，鐘萬里1，謝濤2，徐俊2，鐘力強(qiáng)1，肖曉暉2
(1.廣東電網(wǎng)公司電力科學(xué)研究院，廣東 廣州，510080；
2.武漢大學(xué) 動(dòng)力與機(jī)械學(xué)院，湖北 武漢，430072)

研究傳感器布置對(duì)導(dǎo)波沿薄壁金屬管道的頻散特性、模態(tài)轉(zhuǎn)換和缺陷定位的影響規(guī)律。其步驟為：首先，采用ANSYS軟件建立薄壁小直徑的無缺陷管道和缺陷管道有限元模型；采用瞬態(tài)動(dòng)力學(xué)分析方法，對(duì)管道端部周向各節(jié)點(diǎn)施加軸向瞬時(shí)位移載荷，模擬超聲傳感器激發(fā)L(0,2)模態(tài)波；根據(jù)脈沖-回波信號(hào)的時(shí)間歷程曲線，定義回波分辨率、反射系數(shù)。其次，分析缺陷的軸向位置精度以及4種傳感器布置方案對(duì)導(dǎo)波頻散、模態(tài)轉(zhuǎn)換、周向定位和反射系數(shù)的影響。研究結(jié)果表明：增加傳感器數(shù)量可有效抑制彎曲模態(tài)和導(dǎo)波頻散的發(fā)生，并增強(qiáng)回波信號(hào)的幅值；當(dāng)傳感器數(shù)量小于4個(gè)時(shí)，難以對(duì)裂紋缺陷進(jìn)行周向定位；反射系數(shù)隨傳感器數(shù)量的增加呈線性增大。

管道；導(dǎo)波；傳感器布置；頻散特性；反射系數(shù)

薄壁金屬管道在服役過程中由于腐蝕或老化等因素造成的缺陷常常導(dǎo)致管道發(fā)生泄漏事故，不僅造成經(jīng)濟(jì)損失，而且對(duì)生態(tài)環(huán)境造成威脅，因此，缺陷檢測(cè)對(duì)管道安全至關(guān)重要。針對(duì)金屬管道缺陷，無損檢測(cè)是一種有效的技術(shù)途徑，其中超聲導(dǎo)波因具有傳輸距離遠(yuǎn)、檢測(cè)效率高等優(yōu)點(diǎn)具有廣闊的應(yīng)用前景。近年來，超聲導(dǎo)波檢測(cè)技術(shù)發(fā)展迅速，頻散特性和信號(hào)傳播幅值是管道缺陷檢測(cè)的重要影響因素，而導(dǎo)波傳感器數(shù)量的選擇與布置方案是研究的關(guān)鍵。在導(dǎo)波理論方面，F(xiàn)LETCHER等[1]利用共同源(CSM)的方法合成集中導(dǎo)波，使反射系數(shù)增大1倍，提高了缺陷檢測(cè)的靈敏度；DEMMA等[2]發(fā)現(xiàn)雖然低頻扭轉(zhuǎn)導(dǎo)波對(duì)周向缺陷不敏感，但可以用來檢測(cè)有較大軸向尺寸的缺陷；YUEMIN等[3]推導(dǎo)了導(dǎo)波在管道中傳播方向的理論公式，并進(jìn)行了數(shù)字模擬和實(shí)驗(yàn)研究，得到傳感器布置間隙計(jì)算公式。在模擬和實(shí)驗(yàn)研究方面，ALLEYNE等[4]利用脈沖回波法激勵(lì)超聲波脈沖，并對(duì)由缺陷引起的回波信號(hào)進(jìn)行分析，通過回波幅判斷缺陷的程度；WANG等[5]使用L(0,2)模態(tài)波研究了缺陷兩側(cè)和不同尺寸參數(shù)缺陷的反射信號(hào)，更加準(zhǔn)確地提取管道缺陷尺寸，研究缺陷對(duì)導(dǎo)波的反射規(guī)律；GALVAGNI等[6]等用T(0,1)模態(tài)通過模擬和實(shí)驗(yàn)研究了考慮管道支撐下的管道導(dǎo)波傳播，繪制了頻散曲線，得出區(qū)分缺陷和簡單支撐的方法；YANG等[7]模擬不同頻率模態(tài)波檢測(cè)直管和彎管缺陷，得出管道缺陷定位誤差隨頻率的增大而增大；ALLEYNE等[8]研究指出模態(tài)T(0,1)比模態(tài)L(0,2)更易激發(fā)，但對(duì)于充液管道來說，模態(tài)L(0,2)的信噪比會(huì)降低，極易產(chǎn)生模態(tài)轉(zhuǎn)變而導(dǎo)致誤判；PETER等[9]提出利用MP處理數(shù)據(jù)，不僅提高信噪比，而且可以直接、準(zhǔn)確地判斷裂紋軸向位置。國內(nèi)近年利用超聲導(dǎo)波進(jìn)行管道檢測(cè)的理論和實(shí)驗(yàn)研究。程載斌等[10]對(duì)超聲導(dǎo)波在管道裂紋缺陷檢測(cè)進(jìn)行了數(shù)值模擬，得到了管道的周向裂紋長度和寬度對(duì)回波信號(hào)影響的關(guān)系；董為榮等[11]運(yùn)用有限元分析法，對(duì)目前現(xiàn)場(chǎng)檢測(cè)中應(yīng)用的L(0,2)及T(0,1)模態(tài)導(dǎo)波在管中傳播過程進(jìn)行數(shù)值模擬研究，給出了缺陷回波反射系數(shù)與缺陷橫截面積各影響因素之間的關(guān)系曲線，近似判定缺陷的幾何尺寸；孫廣開等[12]應(yīng)用ANSYS有限元分析程序模擬超聲導(dǎo)波的傳播及遇到缺陷的反射回波情況，較準(zhǔn)確地定位了缺陷位置；孔雙慶[13]模擬了不同模態(tài)、頻率和周期的導(dǎo)波在有缺陷管道中的回波特征，判定管道缺陷損傷程度及其軸向定位；他得安等[14]在自由管材下通過實(shí)驗(yàn)研究了導(dǎo)波幅值和波包寬度隨傳播距離的變化情況；何存富等[15]采用軸對(duì)稱激勵(lì)接收的方法，通過多組實(shí)驗(yàn)對(duì)導(dǎo)波檢測(cè)中傳感器數(shù)量和頻率特性進(jìn)行了研究，并得出反射回波幅值與傳感器的固有特性有關(guān)；馬書義等[16]計(jì)算了不同幾何尺寸管道中 L(0,2)和 L(0,1)導(dǎo)波非頻散段的限制頻率，得出限定或最小化缺陷回波分辨距離可獲得最佳的導(dǎo)波激勵(lì)信號(hào)周期及檢測(cè)頻率或頻段的結(jié)論；王悅民等[17]研究了管道中L(0,2)導(dǎo)波起始截止頻率與管道直徑和壁厚的變化關(guān)系，研究結(jié)論為在整個(gè)頻率范圍內(nèi)進(jìn)一步選擇合適的檢測(cè)頻率提供了參考。從上述研究結(jié)果可知導(dǎo)波的頻散特性和信號(hào)傳播幅值是影響管道檢測(cè)的重要因素，然而目前研究主要集中于管道尺寸對(duì)其影響方面，而對(duì)傳感器數(shù)量對(duì)導(dǎo)波的頻散特性及信號(hào)傳播幅值的影響研究較少，為此，本文作者對(duì)傳感器數(shù)量進(jìn)行研究，以便為管道實(shí)際檢測(cè)提供指導(dǎo)。為了產(chǎn)生沿管道傳播的軸對(duì)稱導(dǎo)波，需要將足夠數(shù)量的傳感器等間距地沿管道圓周方向布置；考慮到電廠實(shí)際管道直徑，傳感器數(shù)量受到限制，本文在不同傳感器數(shù)量的情況下，利用L(0,2)模態(tài)導(dǎo)波對(duì)管道裂紋檢測(cè)進(jìn)行數(shù)值模擬，研究傳感器布置對(duì)導(dǎo)波頻散特性、模態(tài)轉(zhuǎn)換及反射系數(shù)的影響，建立導(dǎo)波反射系數(shù)和傳感器數(shù)量之間的關(guān)系曲線。

1　薄壁管中的導(dǎo)波

1.1導(dǎo)波的模式及頻散現(xiàn)象

管中導(dǎo)波有3種模式：縱向?qū)РＪ?L模式)、彎曲模式(F模式)和扭轉(zhuǎn)模式(T模式)。這3種模式的導(dǎo)波分別用L(n,m)，F(xiàn)(n,m)和T(n,m)表示，其中n和m分別代表周向和徑向模態(tài)參數(shù)且均為整數(shù)。L模式和T模式為軸對(duì)稱模式，F(xiàn)模式為非軸對(duì)稱模式。

由彈性動(dòng)力學(xué)理論首先得到無限長空心圓柱體中導(dǎo)波傳播的頻散方程，通過對(duì)頻散方程求解可得到空心圓管中導(dǎo)波傳播的頻散曲線。由于導(dǎo)波的頻散和多模態(tài)特性，使得導(dǎo)波對(duì)金屬管道的缺陷定位非常困難，因此，必須根據(jù)頻散曲線的特點(diǎn)，選擇特定的導(dǎo)波模態(tài)進(jìn)行數(shù)值模擬。

本文仿真采用的是L(0,2)模態(tài)。該模態(tài)有2個(gè)顯著的優(yōu)點(diǎn)：首先，縱波L(0,2)在40~200 kHz范圍內(nèi)幾乎沒有頻散現(xiàn)象，即導(dǎo)波信號(hào)波包在傳播過程中不會(huì)發(fā)生明顯的變形而導(dǎo)致無法識(shí)別；其次，該模態(tài)傳播的速度最快，因此，最先到達(dá)的就是目標(biāo)信號(hào)，能夠提高特征識(shí)別的精度。

1.2導(dǎo)波基本檢測(cè)理論

導(dǎo)波沿管道傳播的波速只與材料的性能有關(guān)，其縱波波速為[18]

式中：cL為L(0,2)模態(tài)導(dǎo)波在管中的傳播速度，m/s；E為材料彈性模量，Pa；ρ為材料密度，kg/m3；υ為泊松比。

根據(jù)導(dǎo)波信號(hào)時(shí)域圖，可得出初次波與缺陷回波信號(hào)到達(dá)的時(shí)間間隔Δt。設(shè)缺陷離管端距離為x，則有以下關(guān)系：2x=CLΔt。經(jīng)多次測(cè)量，可以精確確定缺陷離管道激勵(lì)端的距離，并且該方法可用作導(dǎo)波波速擬合，比較擬合出的波速是否與理論波速相吻合。

2　薄壁管道中超聲導(dǎo)波的數(shù)值模擬

2.1有限元模型的建立

管道模型包括無缺陷管道和缺陷管道；管道模型幾何參數(shù)包括長度l、內(nèi)徑r和管壁厚h；材料參數(shù)包括彈性模量E、材料密度ρ和泊松比υ。建立缺陷管道模型時(shí)，直接建立帶有裂紋的三維管道模型。為了減少在劃分時(shí)網(wǎng)格單元的畸變，劃分網(wǎng)格前對(duì)管道模型進(jìn)行適當(dāng)切割，以提高模型網(wǎng)格質(zhì)量。缺陷尺寸參數(shù)包括距激發(fā)端面軸向距離l1、周向角度θ、徑向深度R和軸向長度W，如圖1所示。

圖1　管道缺陷尺寸參數(shù)Fig.1　Parameters of pipe defect

有限元模型采用solid45單元類型。根據(jù)單元和時(shí)間參數(shù)計(jì)算公式得模型單元長度LE、總計(jì)算時(shí)間T、時(shí)間步長Δ。為提高計(jì)算精度，缺陷處網(wǎng)格更小。設(shè)網(wǎng)格單元長度為LE0，缺陷管道模型如圖2所示。

實(shí)際檢測(cè)時(shí)，管道有包覆層而且被埋于水下、土壤或其他物質(zhì)中，邊界狀況復(fù)雜，不利于傳感器的安裝。本文定量地研究傳感器數(shù)量和導(dǎo)波反射系數(shù)之間的關(guān)系，因此，對(duì)管道模型的邊界條件進(jìn)行簡化。加載時(shí)，管道模型的內(nèi)外都設(shè)為真空，對(duì)管道的非激發(fā)端面施加固定約束。

圖2　管道模型示意圖Fig.2　Schematic diagram of pipe model

針對(duì)所建模型，采用ANSYS瞬態(tài)動(dòng)力學(xué)模塊中完全瞬態(tài)分析法進(jìn)行超聲導(dǎo)波檢測(cè)數(shù)值模擬。

2.2導(dǎo)波激勵(lì)與接收

在管道左端部激勵(lì)縱向模態(tài)導(dǎo)波L(0,2)，即在端部的節(jié)點(diǎn)上加載對(duì)應(yīng)的位移曲線，選擇激勵(lì)頻率為140 kHz，激勵(lì)信號(hào)采用調(diào)制的10周期正弦信號(hào)。這種激勵(lì)信號(hào)的頻域能量集中，能有效地降低能量的發(fā)散，抑制導(dǎo)波的頻散。激勵(lì)信號(hào)的表達(dá)式為[13]

式中：n為選用的單音頻數(shù)目；fc為信號(hào)中心頻率，Hz；t為時(shí)間，s。

通過對(duì)管道左端的端部周向各節(jié)點(diǎn)預(yù)加軸向瞬時(shí)位移載荷模擬縱向入射導(dǎo)波，監(jiān)測(cè)節(jié)點(diǎn)位移接收反射波，接收信號(hào)位置距導(dǎo)波信號(hào)激發(fā)端1個(gè)單元長度，如圖1與圖2所示。

2.3回波信號(hào)特征提取

裂紋管道的回波信號(hào)時(shí)程曲線如圖3所示。接收節(jié)點(diǎn)檢測(cè)到3組信號(hào)：激勵(lì)信號(hào)x(t)，從裂紋反射回來的回波信號(hào)f(t )，從管道末端反射回來的回波信號(hào)d(t)。

f(t)的波形時(shí)長相對(duì)x(t)的波形時(shí)長有所增加，即出現(xiàn)導(dǎo)波傳播過程的頻散效應(yīng)。為了描述導(dǎo)波信號(hào)在傳播過程中的頻散程度和提取反射回波信號(hào)，定義以下2個(gè)參數(shù)。

1)回波分辨率S，為端部回波信號(hào)d(t)在時(shí)域上的持續(xù)時(shí)間T1與激勵(lì)信號(hào)x(t)在時(shí)域上的持續(xù)時(shí)間T0之比，即S=T1/T0。S越大，頻散越嚴(yán)重。

2)反射系數(shù)C，為反射回波的幅值A(chǔ)1與入射波的幅值A(chǔ)0之比，即C=A1/A0。

導(dǎo)波在管道缺陷和端面處會(huì)發(fā)生反射，在缺陷處的反射系數(shù)大于同周向其他位置的反射系數(shù)。

圖3　缺陷回波時(shí)間歷程Fig.3　Time history of defect echo

3　結(jié)果與分析

3.1仿真工況

管道的幾何尺寸主要是根據(jù)實(shí)際檢測(cè)管道和仿真計(jì)算量而定。本文研究模擬的管道幾何參數(shù)為：l= 1 500 mm，r=76 mm，h=5.5 mm；缺陷尺寸參數(shù)為：l1=1 m，θ=50°，R=5.5 mm，W=1 mm。管道的材料選擇為T24鋼，其材料參數(shù)為：ρ=7.890 t/m3，E=2.0× 1011Pa，υ=0.3。

建立管道幾何模型及確定所采用的單元之后對(duì)管道模型進(jìn)行網(wǎng)格劃分。管道的圓周劃分?jǐn)?shù)由實(shí)驗(yàn)管道環(huán)向布置的傳感器個(gè)數(shù)決定[15]，周向劃分?jǐn)?shù)為加載傳感器數(shù)量的1~3倍。根據(jù)仿真計(jì)算精度和計(jì)算機(jī)配置，模型圓周方向的網(wǎng)格數(shù)為32個(gè)，厚度方向的網(wǎng)格數(shù)為2個(gè)，總網(wǎng)格數(shù)為19 200個(gè)。設(shè)置的計(jì)算參數(shù)為：T= 0.68 ms，LE0=5 mm，Δ=0.8 μs。

建立無缺陷管道和缺陷管道，分別對(duì)對(duì)稱加載4種傳感器數(shù)量的2種模型進(jìn)行數(shù)值模擬。

3.2波速擬合及軸向定位

裂紋的軸向定位是根據(jù)式(2)計(jì)算得到，通過激勵(lì)波和管道端面反射波的到達(dá)時(shí)刻計(jì)算得到Δt= 5.788 2×10-4s，導(dǎo)波的模擬速度為5 171.2 m/s。由式(1)知導(dǎo)波速度理論值為5 277.8 m/s，擬合結(jié)果較好。

根據(jù)激勵(lì)波和管道裂紋反射波的到達(dá)時(shí)刻測(cè)得Δt=4.06×10-4s，裂紋軸向距離為1.05 m。裂紋實(shí)際軸向距離為1 m，缺陷定位誤差為5%，在檢測(cè)誤差范圍內(nèi)。

3.3導(dǎo)波頻散

頻散現(xiàn)象在時(shí)域中的表現(xiàn)為：某一模態(tài)波的包絡(luò)線隨傳播距離的增加不斷拉長，嚴(yán)重降低了檢測(cè)信號(hào)的信噪比，若對(duì)應(yīng)模態(tài)導(dǎo)波頻散現(xiàn)象嚴(yán)重，則激發(fā)信號(hào)的包絡(luò)線也容易變寬。為避免不規(guī)則管道缺陷對(duì)導(dǎo)波頻散的影響，采用無缺陷管道進(jìn)行仿真，對(duì)接收信號(hào)位置處周向32個(gè)節(jié)點(diǎn)的位移時(shí)程曲線進(jìn)行簡單疊加，以消除彎曲波的影響。

布置不同傳感器數(shù)量的無缺陷管道端面反射系數(shù)見圖4。由圖4可知：當(dāng)傳感器數(shù)量由32個(gè)減小到4個(gè)時(shí)，回波分辨率S由1.12增大到1.49，L(0,2)模態(tài)端面反射回波包絡(luò)線變寬，同時(shí)信號(hào)峰值幅度衰減變快，大大影響導(dǎo)波的檢測(cè)距離，這會(huì)對(duì)實(shí)際應(yīng)用造成不利影響，如分辨率、幅度和檢測(cè)距離減小，所以，應(yīng)布置足夠的傳感器以抑制導(dǎo)波頻散。

圖4　不同傳感器數(shù)量的無缺陷管道端面回波時(shí)程曲線Fig.4　Echo time history curve of non-defective pipes end with different sensor numbers

3.4管道缺陷模態(tài)轉(zhuǎn)換

模擬對(duì)稱加載的4種傳感器數(shù)量對(duì)缺陷反射系數(shù)的影響見圖5。分析圖5可知：當(dāng)傳感器的數(shù)量增多時(shí)，激勵(lì)信號(hào)的幅值增大；缺陷反射回波大致相同，都出現(xiàn)了L(0,2)模態(tài)，且波形極清晰；當(dāng)傳感器數(shù)量為8個(gè)和4個(gè)時(shí)缺陷反射回波中還出現(xiàn)了其余模態(tài)的包絡(luò)線。通過缺陷與發(fā)射端面實(shí)際距離計(jì)算出該模態(tài)導(dǎo)波的模擬速度為3 857.9 m/s，對(duì)比頻散曲線可知該模態(tài)為彎曲模態(tài)F(3,3)，數(shù)值模擬結(jié)果與實(shí)驗(yàn)結(jié)果較吻合。分析頻散曲線可知：140 kHz處可激勵(lì)的最高模態(tài)為F(13,1)，傳感器數(shù)量大于在該頻率處激勵(lì)的最高模態(tài)對(duì)應(yīng)的數(shù)值13，由于各彎曲模態(tài)間的相互干涉，會(huì)抑制彎曲模態(tài)波，從而增強(qiáng)軸對(duì)稱模態(tài)波。因此，足夠的傳感器數(shù)量不僅可以增強(qiáng)激勵(lì)信號(hào)的幅值，而且可以有效地抑制彎曲模態(tài)。

圖5　不同傳感器數(shù)量的缺陷回波時(shí)程曲線Fig.5　Echo time history curves of defect with different sensor numbers

3.5裂紋缺陷的周向定位

缺陷的周向定位是根據(jù)管道有限元模型上各信號(hào)接收節(jié)點(diǎn)的周向位置，將各節(jié)點(diǎn)的周向弧度作為橫坐標(biāo)，節(jié)點(diǎn)對(duì)應(yīng)的反射系數(shù)為縱坐標(biāo)，繪制周向反射系數(shù)的極坐標(biāo)圖。通過分析周向反射系數(shù)圖的曲線軌跡對(duì)缺陷進(jìn)行周向定位。

圖6　不同傳感器數(shù)量的反射系數(shù)Fig.6　Reflection coefficients for different sensor numbers

傳感器數(shù)量變化時(shí)裂紋缺陷的周向定位如圖6所示。其中，圖6(b)所示是傳感器數(shù)量為4個(gè)和8個(gè)時(shí)的放大圖。分析圖6(b)可知：不同傳感器數(shù)量的缺陷反射回波的幅值變化趨勢(shì)大致相同，缺陷一側(cè)的反射系數(shù)增大速度明顯高于非缺陷一側(cè)反射系數(shù)的增大速度，因此，缺陷的導(dǎo)波反射系數(shù)在圓周方向的分布特征也越來越明顯，但當(dāng)傳感器數(shù)量為4時(shí)，已不能明顯地區(qū)分缺陷的周向位置，且反射系數(shù)幅值降低，不易識(shí)別。這說明傳感器數(shù)量的變化不會(huì)改變反射回波的幅值變化趨勢(shì)，但隨著傳感器數(shù)量增加，可以有效加強(qiáng)信號(hào)傳播幅值。

3.6傳感器數(shù)量與反射系數(shù)之間的關(guān)系

缺陷管道回波反射系數(shù)見表1。

表1　實(shí)驗(yàn)管道的反射系數(shù)Table 1　Reflection coefficients of pipe

根據(jù)表1繪制缺陷管道反射系數(shù)C和傳感器數(shù)量N的關(guān)系曲線，如圖7所示。分析圖7可知：缺陷管道的反射系數(shù)與感器數(shù)量呈現(xiàn)出良好的線性關(guān)系，這一點(diǎn)與實(shí)驗(yàn)數(shù)據(jù)所表現(xiàn)出的線性關(guān)系相吻合。缺陷管道反射系數(shù)C與傳感器數(shù)量N的線性擬合關(guān)系式為：C=0.010 8N-0.027 7。

圖7　傳感器數(shù)量與反射系數(shù)的關(guān)系Fig.7　Relationship between sensor number and reflection coefficient

4　結(jié)論

1)通過無缺陷管道的脈沖-回波信號(hào)時(shí)域振動(dòng)曲線，計(jì)算L(0,2)模態(tài)波的模擬波速，利用模擬波速進(jìn)行缺陷的軸向定位，得到其定位精度，驗(yàn)證了本文模型的正確性，包括模型的約束、導(dǎo)波的激發(fā)、計(jì)算參數(shù)等。

2)增加傳感器數(shù)量能夠抑制彎曲模態(tài)和導(dǎo)波頻散，提高特征識(shí)別的精度。

3)隨著傳感器數(shù)量的改變，各仿真工況中的缺陷回波幅值變化趨勢(shì)大致相同。當(dāng)傳感器數(shù)量小于4個(gè)時(shí)，缺陷回波反射幅值降低，難以對(duì)裂紋缺陷進(jìn)行周向定位。

4)缺陷管道反射系數(shù)C與傳感器數(shù)量N呈線性變化關(guān)系，其擬合關(guān)系式為：C=0.010 8N-0.027 7。

[1]FLETCHER S,LOWE M J S,RATASSEPP M,et al.Detection of axial cracks in pipes using foused guided waves[J].Journal of Nondestructive Evaluation,2012,31(1):56-64.

[2]DEMMA A,CAWLEY P,LOWE M,et al.The reflection of the fundamental torsional mode from cracks and notches in pipes[J]. Dynamic Method for Damage Detection in Structures,2008, 499(2):195-209.

[3]YUEMIN W,LIHUA S,CHUANJUN S,et al.Theoretical and experimental research on controlling guided waves propagation direction in pipes[C]//Electronic Measurement&Instruments (ICEMI),2011 10th International Conference.Chengdu,China, IEEE,2011:28-32.

[4]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,43(2):180-187.

[5]WANG X,PETER W T,MECHEFSKE C K,et al.Experimental investigation of reflection in guided wave-based inspection for thecharacterizationofpipelinedefects[J].Nondestructive Testing and Evaluation International,2010,43(4):365-374.

[6]GALVAGNI A,CAWLEY P.The reflection of guided waves from simple supports in pipes[J].The Journal of the Acoustical Society ofAmerica,2011,129(4):1869-1880.

[7]YANG H,WANG C.Study on simulation of non-destructive testing for pipeline defects by ultrasonic guided waves[C]// CrossStraitQuad-RegionalRadioScienceandWireless Technology Conference(CSQRWC).Harbin,China:IEEE,2011: 238-242.

[8]ALLEYNE D N,VOGT T,CAWLEY P.The choice of torsional or longitudinal excitation in guided wave pipe inspection[J]. Insight-Non-Destructive Testing and Condition Monitoring, 2009,51(7):373-377.

[9]PETER W T,WANG X.Characterization of pipeline defect in guided-waves based inspection through matching pursuit with the optimized dictionary[J].NDT&E International,2013,54: 171-182.

[10]程載斌,王志華,張立軍,等.管道超聲縱向?qū)Р鸭y檢測(cè)數(shù)值模擬[J].應(yīng)用力學(xué)學(xué)報(bào),2005,21(4):76-79. CHENGZaibin,WANGZhihua,ZHANGLijun,etal. Numerical simulation of crack detection in pipes using ultrasonic longitudinalguided-wave[J].ChineseJournalofApplied Mechanics,2005,21(4):76-79.

[11] 董為榮,帥健,許葵.管道腐蝕缺陷超聲導(dǎo)波檢測(cè)數(shù)值模擬研究[J].機(jī)械強(qiáng)度,2009,30(6):988-993. DONG Weirong,SHUAI Jian,XU Kui.Numerical simulation of corrosion detection in pipes using ultrasonic guided waves[J]. Journal of Mechanical Strength,2009,30(6):988-993.

[12]孫廣開,焦陽,李光海,等.超聲導(dǎo)波管道缺陷檢測(cè)數(shù)值模擬[J].河北工業(yè)科技,2010,27(1):18-21. SUN Guangkai,JIAO Yang,LI Guanghai,et al.Numerical simulation of defect detection in pipes using ultrasonic guided waves[J].Hebei Journal of Industrial Science and Technology, 2010,27(1):18-21.

[13]孔雙慶.管道超聲導(dǎo)波檢測(cè)的數(shù)值模擬和實(shí)驗(yàn)研究[D].大連:大連理工大學(xué)工程力學(xué)系,2011:24-41. KONG Shuangqing.Numerieal simulation and experimental investigation of defect detection in pipes using ultrasonic guided waves[D].Dalian:Dalian University of Technology.Department of Engineering Mechanics,2011:24-41.

[14]他得安,易勇,劉鎮(zhèn)清.傳播距離對(duì)管中導(dǎo)波傳播特性的影響[J].無損檢測(cè),2004,25(11):553-556. TA Dean,YI Yong,LIU Zhenqing.Influence of propagation distance on the characteristics of ultrasonic guided wave propagation in pipes[J].NDT,2004,25(11):553-556.

[15]何存富,劉增華,鄭璟瑜,等.管道導(dǎo)波檢測(cè)中傳感器數(shù)量和頻率特性研究[J].北京工業(yè)大學(xué)學(xué)報(bào),2005,30(4):393-397. HE Cunfu,LIU Zenghua,ZHENG Jingyu,et al.Study on the number and the frequency characteristic of transducers in pipe inspection using guided waves[J].Journal of Beijing University of Technology,2005,30(4):393-397.

[16]馬書義,武湛君,劉科海,等.空心圓管中導(dǎo)波頻散特性與檢測(cè)頻率選擇[J].機(jī)械工程學(xué)報(bào),2014,50(20):8-17. MAShuyi,WUZhanjun,LIUKehai,etal.Dispersion characteristics and inspection frequency selection of guided waves in hollow cylinders[J].Journal of Mechanical Engineering, 2014,50(20):8-17.

[17]王悅民,沈立華,申傳俊,等.管道導(dǎo)波無損檢測(cè)頻率選擇與管材特征關(guān)系[J].機(jī)械工程學(xué)報(bào),2009,45(8):243-248. WANGYuemin,SHENLihua,SHENChuanjun,etal. Relationshipbetweenfrequencyselectionandpipe characteristics for guided waves NDT[J].Journal of Mechanical Engineering,2009,45(8):243-248.

[18]董為榮,帥健,許葵.管道T(0,1)模態(tài)導(dǎo)波檢測(cè)數(shù)值模擬研究[J].無損檢測(cè),2008,30(3):149-152. DONG Weirong,SHUAI Jian,XU Kui.Numerical simulation study of detection in pipes using T(0,1)mode ultrasonic guided waves[J].NDT,2008,30(3):149-152.

(編輯陳燦華)

Optimization of guided wave sensors distribution along thin-walled small-diameter pipe

WANG Wei1,YOU Penghui2,ZHONG Wanli1,XIE Tao2,XU Jun2,ZHONG Liqiang1,XIAO Xiaohui2
(1.Electric Power Research Institute of Guangdong Power Grid Company,Guangzhou 510080,China;
2.College of Power and Mechanical Engineering,Wuhan University,Wuhan 430072,China)

The influence of sensor distribution scheme on the frequency dispersion characteristics,modal conversion and defect localization of guided wave along thin-walled metal pipe was investigated.The procedures were as follows.Firstly, the models of defective and non-defective finite element pipes were built using ANSYS software.Axial displacement was exerted to circumferential nodes in the end of pipe based on the transient dynamic analysis methods,and ultrasonic sensor activating L(0,2)mode was simulated.Then,according to time-history curves of the pulse-echo,characteristic parameters such as echo resolution and reflection coefficient were defined.Secondly,the axial defect location precision,as well as the influence of four sensor distribution schemes on the guided wave dispersion,modal conversion,circumferential localization and reflection coefficient were analyzed.The results show that bending modal and wave dispersion can be inhibited effectively while amplitude of echo signal can be amplified by adding sensors.It is difficult to locate crack defects circumferentially when the number of sensors is less than 4.And the reflection coefficient increases linearly with the increase of the number of sensors.

pipe;guided wave;sensor distribution;dispersion characteristics;reflection coefficient

肖曉暉，博士，教授，從事特種機(jī)器人與微操作機(jī)器人研究；E-mail:xhxiao@whu.edu.cn

TM8；TM113

A

1672-7207(2016)07-2254-06

10.11817/j.issn.1672-7207.2016.07.011

2015-09-31；

2015-11-22

國家自然科學(xué)基金資助項(xiàng)目(51175383)(Project(51175383)supported by the National Natural Science Foundation of China)