脈沖渦流多層厚度檢測與時頻特征量提取*

王志春， 王玉璽

(內(nèi)蒙古科技大學(xué) 信息工程學(xué)院，內(nèi)蒙古 包頭 014010)

0 引 言

金屬板材厚度是產(chǎn)品質(zhì)量、設(shè)備服役狀態(tài)的重要參數(shù)，如金屬板(銅、鋁、鋼板等)軋制過程中的厚度檢測、各類金屬中缺陷的檢測等[1]，因此，準確、快速測量金屬板材厚度至關(guān)重要。電磁渦流檢測是無損檢測領(lǐng)域的五大常規(guī)檢測方法之一( 射線、超聲、磁粉、滲透和渦流)，具有實現(xiàn)非接觸檢測，能穿透被測件的覆蓋層等優(yōu)點[2]。其中脈沖渦流測厚法不需要耦合劑，其可在拆除保溫層的管道壁厚的檢測得到使用，也可在高溫等惡劣環(huán)境下使用，而且操作方便、檢測速度快、易于實現(xiàn)自動化，在金屬材料的厚度檢測中應(yīng)用十分廣泛[3]。因此，在進行金屬缺陷和厚度的檢測的研究中，脈沖渦流有其重要的研究意義[4]。本文利用激勵脈沖頻率成分豐富的特點，深入開展了時域和頻域特征提取技術(shù)的研究，更好地分析被測材料的信息[5]。

1 理論基礎(chǔ)

1.1 脈沖渦流測厚理論

脈沖信號源產(chǎn)生具有一定占空比的方波，并加到激勵線圈兩端，此時在激勵線圈中就會產(chǎn)生周期的脈沖電流。激勵線圈中的電流便會感生出一個快速衰減的脈沖磁場，變化著的磁場在被測材料中感應(yīng)出瞬時渦流，向被測材料內(nèi)部傳播，并感應(yīng)出一個快速衰減的渦流磁場[6]。隨著渦流磁場的在導(dǎo)體內(nèi)的變化，再檢測線圈則會感應(yīng)出隨時間變化的電壓。由于脈沖渦流在被測材料內(nèi)部的傳播過程是逐漸衰減的，對于不同厚度的導(dǎo)體試件，得到的瞬態(tài)感應(yīng)電壓信號不同。所以，通過對檢測到的電壓信號進行分析，找到感應(yīng)電壓與被測導(dǎo)體試件厚度的相應(yīng)的規(guī)律，實現(xiàn)對導(dǎo)體材料的厚度檢測[7]。

1.2 脈沖渦流滲透深度分析

由渦流的趨膚效應(yīng)可知，渦流總是集中在接近線圈的導(dǎo)體表面。對于不同電導(dǎo)率金屬材料，在導(dǎo)體的橫截面，其電流密度分布也呈現(xiàn)出很大差異，從工件表面向其內(nèi)部衰減，并且是按負指數(shù)規(guī)律衰減

(1)

式中J0為導(dǎo)體表面的渦流密度;Jx為距導(dǎo)體表面為x處的渦流密度;f為交流電流的頻率;μ為材料的磁導(dǎo)率；σ為材料的電導(dǎo)率。

為了更確切地表明脈沖渦流在被測導(dǎo)體內(nèi)的衰減情況，引入了表示衰減特性的物理量:滲透深度δ。定義為:渦流密度衰減到導(dǎo)體表面處的1/e時，距離導(dǎo)體表面的距離

(2)

渦流的檢測深度由滲透深度δ定性反映[8]，δ變大，對應(yīng)的渦流檢測的深度也會隨之增大；反之，則會減小。但δ并不是渦流檢測所能達到的最大深度。渦流所能檢測的最大深度可達2.6倍的滲透深度，在此深度時，渦流密度已經(jīng)衰減得很小，且為表面密度的5%。通常情況下，定義2.6δ的滲透深度作為渦流的有效滲透深度。根據(jù)分析，定義渦流有效滲透深度為L，即

(3)

由式(3)可得，對不同電導(dǎo)率的金屬材料有其對應(yīng)的最大檢測厚度，進而確定在相應(yīng)的測厚范圍內(nèi)應(yīng)用脈沖渦流檢測技術(shù)[9]。

2 建模與仿真

2.1 建 模

對于被測材料，選取兩種不同電導(dǎo)率的金屬(45#鋼和鋁)疊放，進而測得不同厚度金屬對應(yīng)的感應(yīng)差壓信號。圖1為經(jīng)過網(wǎng)格剖析的差分傳感器模型。

圖1 差分傳感器模型

在不同檢測環(huán)境下，多層金屬材料厚度的測定由于其屬性不同，增加了檢測難度，為了測定金屬材料的厚度，本文根據(jù)被測材料電導(dǎo)率的不同，將其分為上下兩層的被測導(dǎo)體，并設(shè)計兩種不同的方案：1)上層被測金屬材料的電導(dǎo)率小于下層被測金屬材料，上層金屬材料選定為45#鋼，下層金屬材料選定為鋁，如圖2(a)；2)上層被測金屬材料的電導(dǎo)率大于下層被測金屬材料，上層金屬材料選定為鋁，下層金屬材料選定為45#鋼，如圖2(b)。

圖2 兩層不同材料的被測導(dǎo)體

2.2 仿真參數(shù)選取

由脈沖渦流測厚原理，將被測導(dǎo)體厚度的計算轉(zhuǎn)換為找到檢測線圈上的電壓信號。因此，只要能夠找出檢測線圈上的電壓與被測材料厚度之間的對應(yīng)關(guān)系，即可實現(xiàn)脈沖渦流對金屬導(dǎo)體的厚度測量。 在保持提離距離、激勵信號等條件不變的情況下，采用不同厚度的被測材料，得到不同厚度被測導(dǎo)體對應(yīng)的電壓，進而找出被測導(dǎo)體的厚度與電壓之間的變化關(guān)系[10]。仿真參數(shù)以及金屬材料參數(shù)如表1和表2所示 。

表1 仿真參數(shù)

表2 金屬材料的參數(shù)

3 數(shù)據(jù)處理與分析

3.1 仿真結(jié)果

由表1所選的仿真參數(shù)，并且根據(jù)仿真時所要求的約束條件，通過Ansoft仿真[11]，即可得到在不同厚度時對應(yīng)的電壓變化規(guī)律。本文僅列出了2 mm 45#鋼、2 mm及0.50 mm鋁、1 mm 45#鋼的金屬材料的波形，觀察檢測端及參考端電壓的仿真結(jié)果。

1)激勵脈沖電壓

Ansoft仿真時所選的激勵電壓的參數(shù)：脈沖電壓幅值-5～5 V，占空比50 %，激勵頻率1 kHz。

“用曲的話來說，這難道不是老公的義務(wù)嗎？但這打死也不能說。Y反復(fù)強調(diào)，S跟她說過的，跟曲的婚姻沒有愛情，也幾乎沒有性生活，曲對S沒有熱情，而S也不渴念她的身體——就是因為這句話，Y才決定跟他好的，否則，她還搞什么呢，絕不會讓S碰她。”

2)檢測端與參考端電壓

在圖3中繪出檢測端以及參考端的電壓波形。

圖3 檢測端與參考端線圈的電壓

3.2 數(shù)據(jù)分析

3.2.1 時域信號特征值提取

對脈沖渦流信號進行時域分析，通常采用差分信號進行信號的特征提取,主要有脈沖渦流差分信號的峰值、過零時間等。圖4為檢測端與參考端的感應(yīng)差分信號，其中比較兩圖可以發(fā)現(xiàn)，在同一時刻，其差分信號所對應(yīng)的幅值相反。

圖4 檢測端與參考端差分感應(yīng)信號

圖5為不同厚度對應(yīng)的峰值，可以看出，在測量不同厚度的金屬材料時，對于上層材料為45#鋼(下層材料為鋁)且厚度逐漸變化時，其對應(yīng)的峰值呈現(xiàn)一定的線性關(guān)系，并且隨著厚度的增加峰值逐漸增加；而對于上層材料為鋁(下層材料為45#鋼)，當厚度逐漸增加時，對應(yīng)的峰值呈現(xiàn)非線性的關(guān)系。圖6為不同厚度的金屬材料對應(yīng)的過零時間，可以得到當上層金屬材料為鋁(下層材料45#鋼)，且金屬鋁的厚度逐漸變換時，其對應(yīng)的過零時間與厚度呈現(xiàn)線性關(guān)系；而對于上層材料為45#鋼(下層材料鋁)，當45#鋼的厚度逐漸變化時，其過零時間幾乎保持不變。

圖5 不同厚度的輸出峰值電壓

圖6 不同厚度的過零時間

由此可見，針對不同電導(dǎo)率層金屬的測定，有其對應(yīng)不同的特征值。在上層電導(dǎo)率小于下層電導(dǎo)率的情況下，其電壓幅值作為特征值比較理想，而過零時間不明顯；在上層電導(dǎo)率大于下層電導(dǎo)率的情況下，其過零時間作為特征值比較理想，而電壓峰值呈現(xiàn)非線性關(guān)系。

3.2.2 頻域信號特征值提取

脈沖渦流擁有很寬的頻寬，當激勵信號經(jīng)過被測材料被檢測線圈捕獲后，檢測信號的頻譜中即包含了金屬材料厚度的信息[12]。由傅里葉級數(shù)展式可得，一個脈沖激勵信號可以被分解為含有基波和許多諧波成分的組合

(4)

截取一個周期的差分信號，經(jīng)過快速傅立葉變換(fast Fourier transform,FFT)后得到相頻和幅頻特性曲線，圖7(a)為不同厚度鋼的相頻特性曲線，在圖中每一條曲線均存在一個相位過零點[13]。

圖7(b)為相位過零點的放大圖，從圖中可以得到不同厚度對應(yīng)的零點逐漸左移，在鋼的厚度為1～5 mm時與對應(yīng)的厚度呈線性關(guān)系。圖7(c)為不同厚度鋁的相頻特性曲線，圖7(d)為經(jīng)過放大得到的相頻圖，在圖中也能夠得到對應(yīng)的相位過零點隨著厚度的增加左移，在鋁的厚度為0.25～1.25 mm時與對應(yīng)的厚度呈線性關(guān)系。

圖7 不同厚度的相頻特性曲線

從圖7中可以看出：當測量兩層金屬材料時，若下層電導(dǎo)率較大，基頻幅值能夠作為特征值；反之，則不能；比較圖7(b)和圖7(d)可知，對于不同電導(dǎo)率層的被測金屬，在一定范圍內(nèi)，相位也可作為其特征值。

圖8為不同厚度的材料經(jīng)過FFT之后的基頻幅值，在不斷改變45#鋼厚度(下層金屬鋁厚度不變)時，基頻幅值隨之增大，而當改變鋁的厚度(下層金屬45#鋼厚度不變)時，鋁的厚度增加時基頻幅值反而減小。

圖8 不同厚度的基頻幅值

4 結(jié) 論

1)利用脈沖渦流對電導(dǎo)率不同的多層金屬材料進行厚度測定時，當被測材料的上層電導(dǎo)率小于下層電導(dǎo)率時，檢測響應(yīng)信號的峰值、基頻幅值可作為特征，并且基頻幅值的線性度為最佳，而當被測材料的上層電導(dǎo)率大于下層電導(dǎo)率時，檢測響應(yīng)信號的過零時間可作為特征值。無論被測導(dǎo)體的電導(dǎo)率大小如何，相位均可作為其特征值，但相位作為特征值時其對應(yīng)的測量范圍較小。

2)脈沖渦流對被測材料進行測定時，對于上下兩層不同電導(dǎo)率的分布，得到的差分響應(yīng)信號在同一時刻對應(yīng)的幅值相反，如在0.5 ms時刻，被測導(dǎo)體的上層電導(dǎo)率小于下層電導(dǎo)率時幅值為負值，而對于被測導(dǎo)體的上層電導(dǎo)率大于下層電導(dǎo)率時幅值為正值。為在不同金屬電導(dǎo)率厚度的測定提供了理論依據(jù)。

[1] 張 震.脈沖渦流測厚技術(shù)研究[D].西安:西安理工大學(xué),2007.

[2] 李培源，何宇廷，杜金強，等.基于柔性渦流傳感器疲勞裂紋監(jiān)測試驗研究[J].傳感器與微系統(tǒng)，2015,34(1):24-27.

[3] Abidin I Z，Tian G Y，Wilson J，et al． Quantitative evaluation of angular defects by pulsed eddy current ther-mography[J]．NDT&E International，2010，43:537-546．

[4] 劉世杭，過玉清.脈沖渦流感應(yīng)測量方法的研究[J].東華大學(xué)學(xué)報,2001,27(6):49-52.

[5] 李 斌，孫忠林，彭芳瑜,等．葉片零件在線檢測測點規(guī)劃問題研究[J].中國機械工程，2008,19(5)：540-542.

[6] 周德強，田貴云，尤麗華,等.方向性脈沖渦流應(yīng)力檢測研究進展[J].傳感器與微系統(tǒng),2011,30(9):1-2.

[7] 周 昊.基于脈沖渦流信號主成分分析的金屬測厚研究[D].西安:西安理工大學(xué)，2008.

[8] 柯 海.脈沖渦流測厚信號斜率法研究[D].武漢.華中科技大學(xué)，2013.

[9] 康學(xué)福.基于脈沖渦流信號的金屬膜厚測量[D].天津:天津大學(xué)，2011.

[10] Ansoft在工程電磁場中的應(yīng)用[M].北京:中國水利水電出版社，2013.

[11] 徐晨曦，闕沛文，毛義梅.脈沖渦流檢測系統(tǒng)的設(shè)計[J].傳感器與微系統(tǒng)，2009,28(6):95-96.

[12] 周德強，張斌強，田貴云,等． 脈沖渦流檢測中裂紋的深度定量及分類識別[J]．儀器儀表學(xué)報，2009，30(6):1190-1194．

[13] 高軍哲.多頻渦流無損檢測的干擾抑制和缺陷檢測方法研究[D].長沙:國防科技大學(xué)，2011.