周德強(qiáng)，閆向陽(yáng)，尤麗華，張秋菊

（1．江南大學(xué)機(jī)械工程學(xué)院，江蘇無(wú)錫 214122;2．無(wú)錫國(guó)盛精密模具有限公司，江蘇無(wú)錫 214024）

0 引言

脈沖渦流無(wú)損檢測(cè)方法是目前渦流檢測(cè)領(lǐng)域的一個(gè)研究熱點(diǎn)，激勵(lì)脈沖包含的頻譜寬，可提取的缺陷信息豐富［1］，對(duì)深層缺陷檢測(cè)能力強(qiáng)，應(yīng)用前景廣闊。傳統(tǒng)的圓柱形脈沖渦流傳感器可以對(duì)表面和近表面缺陷進(jìn)行定量檢測(cè)［2］，但激勵(lì)磁場(chǎng)直接穿過(guò)檢測(cè)線圈，使得激勵(lì)信號(hào)對(duì)檢測(cè)信號(hào)產(chǎn)生很大的干擾。文獻(xiàn)［3］指出，在矩形載流線圈激勵(lì)下，感應(yīng)電磁場(chǎng)在金屬鋁板中的衰減率較傳統(tǒng)的渦流激勵(lì)線圈要小，矩形載流線圈激勵(lì)的渦流滲透深度更深。此外，脈沖激勵(lì)矩形線圈具有如下優(yōu)點(diǎn)［4～7］:在裂紋檢測(cè)中由于提離效應(yīng)影響小而對(duì)微小裂紋具有高靈敏度;能有效消除邊緣效應(yīng)的影響;由于具有方向特性，矩形探頭能夠產(chǎn)生均勻的渦流分布等。我國(guó)學(xué)者何赟澤等人也對(duì)脈沖渦流矩形傳感器開(kāi)展了大量研究，提出了蝶形圖新型特征量等對(duì)缺陷進(jìn)行定量與分類識(shí)別［8～10］。因此，開(kāi)展脈沖渦流矩形傳感器的優(yōu)化設(shè)計(jì)十分必要。

脈沖渦流矩形傳感器可采用2種方式進(jìn)行檢測(cè):脈沖渦流矩形傳感器平行與垂直被檢測(cè)對(duì)象。本文以脈沖渦流矩形傳感器平行于被檢測(cè)對(duì)象為模型，采用Comsol Multiphysics有限元仿真軟件建立了脈沖渦流矩形傳感器檢測(cè)模型，通過(guò)電導(dǎo)率的改變，提取磁場(chǎng)的瞬態(tài)信號(hào)。通過(guò)改變矩形傳感器的尺寸參數(shù)即改變長(zhǎng)、寬、高的比值，分析傳感器尺寸參數(shù)的變化對(duì)傳感器靈敏度和線性度的影響。然后根據(jù)仿真結(jié)果，設(shè)計(jì)制作了3種尺寸的矩形探頭，并使用這些探頭對(duì)鋁板的表面缺陷進(jìn)行了檢測(cè)，通過(guò)對(duì)比仿真和實(shí)驗(yàn)的結(jié)果，得出了矩形傳感器的最佳尺寸參數(shù)。

1 仿真模型的建立和參數(shù)設(shè)置

1．1 仿真模型的建立

本文采用Comsol Multiphysics建立了脈沖渦流矩形探頭傳感器三維模型。在AC/DC模塊下選擇物理場(chǎng)為磁場(chǎng)，并且在瞬態(tài)求解下進(jìn)行仿真分析。整體建模結(jié)構(gòu)如圖1所示。

圖1 脈沖渦流矩形探頭有限元仿真模型圖Fig 1 Finite element simulation model of pulsed eddy current rectangular probe

1．2 電流設(shè)定與其他參數(shù)說(shuō)明

圖2為線圈的俯視圖。激勵(lì)電流密度大小定義變量為curd，curd=N×0．5 ×sgnt/A，N=T×H/d2。其中，N為線圈的匝數(shù)，0．5為設(shè)定的電流大小，t為時(shí)間，A為線圈的橫截面積，T表示線圈的厚度，實(shí)驗(yàn)設(shè)定為3 mm，H為線圈的高度，d為漆包線直徑，擬定為0．3 mm。按照逆時(shí)針?lè)较蛟O(shè)定電流密度方向，磁場(chǎng)方向z軸為正向。按照?qǐng)D示①～⑧的順序?yàn)楦鱾€(gè)線圈添加外部電流密度，此模型的電流只需要添加x，y2個(gè)方向，z方向始終為0。如模塊1，只在x方向有電流，故x方向外部電流密度設(shè)為curdA/m2，y方向?yàn)?;模塊4其電流方向?yàn)?35°角，則x方向外部電流密度為－curd/sqrt（2）A/m2，y方向?yàn)閏urd/sqrt（2）A/m2，依次完成對(duì)矩形線圈的電流設(shè)定。

圖2 矩形線圈的俯視圖Fig 2 Top view of the rectangular coil

被測(cè)試樣為鋁板，其尺寸（長(zhǎng)×寬×高）為100 mm×100 mm×10 mm，相對(duì)磁導(dǎo)率為1，電導(dǎo)率為 σ=3．774×107S/m。激勵(lì)電流幅值大小為0．5 A，頻率為100 Hz。傳感器放置在試樣上方0．5 mm處，即提離高度為0．5 mm。

2 仿真實(shí)驗(yàn)與結(jié)果分析

本文采用單因素輪換法進(jìn)行仿真實(shí)驗(yàn)，即控制矩形線圈的2個(gè)變量不變，改變另外一個(gè)變量來(lái)探究傳感器的靈敏度。首次實(shí)驗(yàn)保持長(zhǎng)寬相等，改變高度來(lái)仿真，即按比例變化為1∶1∶h，h表示變化的高度，取值從1變化到2，每次遞增0．25。比例基值為20 mm，線圈厚度為3 mm。仿真安排如表1所示。

表1 高度變化的激勵(lì)線圈尺寸參數(shù)表Tab 1 Excitation coil size parameters of height change

每次仿真求取0．005 s內(nèi)500個(gè)點(diǎn)的磁場(chǎng)變化情況，提取試樣上方0．5 mm處的磁通密度z軸分量作為研究對(duì)象。電導(dǎo)率沒(méi)有變化時(shí)的數(shù)據(jù)，稱之為參考信號(hào)，電導(dǎo)率逐次增加10%求取信號(hào)。采用各個(gè)電導(dǎo)率變化仿真得到的信號(hào)減去參考信號(hào)得到差分信號(hào)。差分信號(hào)有明顯的峰值信息，可以有效地提高傳感器的靈敏度。如圖3所示，①為參考信號(hào)，即電導(dǎo)率沒(méi)有變化的信號(hào);②為電導(dǎo)率變化50%的信號(hào);③為參考信號(hào)減去電導(dǎo)率變化信號(hào)得到的差分信號(hào)。

圖3 參考信號(hào)與差分信號(hào)等關(guān)系圖Fig 3 Reference signal and differential signal etc

圖4為提取不同電導(dǎo)率變化的差分信號(hào)峰值，該曲線反應(yīng)的是當(dāng)長(zhǎng)、寬、高比值為1∶1∶2的時(shí)差分信號(hào)結(jié)果。其斜率反映矩形線圈傳感器的靈敏度與線性度。通過(guò)對(duì)比不同尺寸參數(shù)下的差分信號(hào)峰值特征，得出不同尺寸傳感器的優(yōu)劣，進(jìn)而優(yōu)化選擇出最佳的尺寸參數(shù)。

圖4 差分信號(hào)峰值隨電導(dǎo)率變化Fig 4 Differential signal peak value vs conductivity changes

圖5為長(zhǎng)寬比例不變，高度變化的綜合圖。通過(guò)對(duì)比可知在長(zhǎng)、寬、高比例為1∶1∶1．5時(shí)靈敏度最高，線性度也最好，即為最優(yōu)比值。

同理對(duì)長(zhǎng)和寬的變化實(shí)施第二組實(shí)驗(yàn)，此次保持寬和高的比為1∶1，變化長(zhǎng)度進(jìn)行優(yōu)化選擇。實(shí)驗(yàn)安排如表2。

同樣在不同尺寸下提取差分信號(hào)峰值進(jìn)行分析，其結(jié)果如圖6所示。由于長(zhǎng)、寬因素在線圈中無(wú)特定量，即長(zhǎng)、寬為等價(jià)的，通過(guò)對(duì)比長(zhǎng)、寬時(shí)的曲線數(shù)據(jù)可知，當(dāng)長(zhǎng)、寬比例為2∶1時(shí)為最優(yōu)比值。

圖5 不同高度下的差分信號(hào)峰值隨電導(dǎo)率的變化Fig 5 Peak value of differential signal with height changes and conductivity changes

表2 長(zhǎng)度變化激勵(lì)線圈尺寸參數(shù)表Tab 2 Excitation coil size parameters of height change

因此，當(dāng)長(zhǎng)、寬、高綜合考慮時(shí)，由于高是相對(duì)于長(zhǎng)、寬的1．5 倍，而長(zhǎng)、寬優(yōu)化比為 2∶1，因此，對(duì)比長(zhǎng)、寬、高比值為 2∶1∶1．5，2∶1∶3，2∶1∶2．25 時(shí)的數(shù)據(jù)，即長(zhǎng)、寬、高尺寸分別為（40，20，30 mm），（40，20，60 mm），（40，20，45 mm），其結(jié)果如圖7所示。從圖7可知，當(dāng)長(zhǎng)、寬、高比值為2∶1∶1．5時(shí)，脈沖渦流矩形傳感器的靈敏度、線性度最好。

圖6 長(zhǎng)寬比值變化時(shí)差分信號(hào)峰值圖Fig 6 Differential signal peak value diagram of length and width ratio change

圖7 三組比值的差分信號(hào)峰值圖Fig 7 Differential signal peak value diagram of three groups of ratio

3 實(shí)驗(yàn)驗(yàn)證與結(jié)果分析

3．1 實(shí)驗(yàn)驗(yàn)證系統(tǒng)設(shè)計(jì)

根據(jù)仿真結(jié)果，實(shí)驗(yàn)制作了3個(gè)不同尺寸的矩形線圈骨架。按照仿真的最佳三組比例，以20 mm為基值，3個(gè)線圈骨架的尺寸分別為40 mm×20 mm×30 mm，40 mm×20 mm×45mm，40mm×20mm×60mm。使用直徑為0．36mm的漆包線繞制的線圈厚度與仿真使用的厚度相同，為3 mm。驗(yàn)證實(shí)驗(yàn)采用試樣的缺陷深度變化來(lái)代替電導(dǎo)率的變化。被測(cè)試樣為鋁板，尺寸（長(zhǎng)×寬×厚）為150 mm×50 mm×10 mm。缺陷寬度為 2 mm，深度依次為 2，4，6，8 mm，實(shí)驗(yàn)試塊如圖8所示。板材厚度為10 mm，缺陷可認(rèn)為電導(dǎo)率為0，缺陷檢測(cè)實(shí)驗(yàn)近似模擬被測(cè)試件電導(dǎo)率的變化。因此，可以通過(guò)此方法來(lái)驗(yàn)證仿真結(jié)果。實(shí)驗(yàn)使用3個(gè)探頭分別在沒(méi)有缺陷和不同深度的缺陷處進(jìn)行測(cè)量和數(shù)據(jù)采集。

圖8 試樣缺陷示意圖Fig 8 Schematic diagram of specimen defects

3．2 實(shí)驗(yàn)結(jié)果分析

本實(shí)驗(yàn)使用DAQ2010數(shù)據(jù)采集卡進(jìn)行數(shù)據(jù)采集。采樣頻率為10 kHz，脈沖波的上升沿觸發(fā)，采樣點(diǎn)數(shù)為1 000點(diǎn)，共采集10個(gè)周期的數(shù)據(jù)。對(duì)信號(hào)進(jìn)行平均處理，避免單個(gè)周期的偶然性，提高實(shí)驗(yàn)的精確度。采集信號(hào)后的處理步驟如圖9。

圖9 數(shù)據(jù)處理流程圖Fig 9 Flow chart of data processing

如圖10所示，原始信號(hào)中由于直流分量的存在，導(dǎo)致信號(hào)電壓整體均在零以上。

圖10 采集到的原始信號(hào)Fig 10 Collected original signal

圖11為原始信號(hào)減去各自的均值的結(jié)果，可以看出:對(duì)電路信號(hào)中的直流分量進(jìn)行了去除。由于信號(hào)中存在一定的噪聲，所以，必須經(jīng)過(guò)濾波處理。經(jīng)過(guò)軟件濾波然后獲取差分信號(hào)如圖12所示，信號(hào)的噪聲明顯減小。將圖12的10個(gè)周期的信號(hào)平均到一個(gè)周期，這樣可以進(jìn)一步消除噪聲。其周期平均濾波結(jié)果如圖13所示。

圖11 原始信號(hào)減去均值得到的信號(hào)Fig 11 Signal obtained from original signal subtracted by the mean value

圖12 差分信號(hào)Fig 12 Differential signal

從圖13可以看出:隨著缺陷深度的增加，差分信號(hào)峰值增大。為此，本文提取差分信號(hào)峰值進(jìn)行靈敏度與線性度分析。其3個(gè)矩形脈沖渦流探頭的實(shí)驗(yàn)結(jié)果如圖14所示。從圖14中可以看出:尺寸為40 mm×20 mm×30 mm的傳感器線性度和靈敏度最好。實(shí)驗(yàn)結(jié)果驗(yàn)證了仿真結(jié)果的有效性，證實(shí)當(dāng)長(zhǎng)、寬、高比值為2∶1∶1．5時(shí)，脈沖渦流矩形傳感器的靈敏度、線性度最好。

圖13 平均到一個(gè)周期的信號(hào)Fig 13 Signal of averaged one period

圖14 不同尺寸傳感器的靈敏度分析Fig 14 Sensitivity analysis of different size sensors

4 結(jié)論

本文通過(guò)仿真和實(shí)驗(yàn)探究了脈沖渦流矩形傳感器尺寸參數(shù)對(duì)對(duì)鋁板檢測(cè)的靈敏度影響規(guī)律。研究發(fā)現(xiàn):脈沖渦流矩形傳感器靈敏度、線性度與傳感器長(zhǎng)寬、高度有關(guān)。在長(zhǎng)寬不變的情況下，當(dāng)長(zhǎng)、寬、高比為1∶1∶1．5時(shí)，傳感器的線性度和靈敏度最高;在高不變時(shí)，長(zhǎng)、寬比為2∶1時(shí)傳感器最優(yōu)。通過(guò)對(duì)比 2∶1∶1．5，2∶1∶3，2∶1∶2．25 得到長(zhǎng)、寬、高最佳比為2∶1∶1．5。本文通過(guò)對(duì)矩形傳感器尺寸的仿真和實(shí)驗(yàn)，找出了最優(yōu)化比值，有利于進(jìn)一步提高脈沖渦流矩形傳感器的檢測(cè)性能。

［1］ 楊賓峰，張 輝，趙玉豐，等．基于新型脈沖渦流傳感器的裂紋缺陷定量檢測(cè)技術(shù)［J］．空軍工程大學(xué)學(xué)報(bào):自然科學(xué)版，2011，12（1）:73 －77．

［2］ 徐志遠(yuǎn)，武新軍，黃 琛，等．激勵(lì)參數(shù)和試樣電磁參數(shù)對(duì)脈沖渦流檢測(cè)的影響的仿真分析［J］．無(wú)損檢測(cè)，2011，30（6）:1－4．

［3］ Hoshikawa H，Koyama K．Eddy current distribution using parameters normalized by standard penetration depth［J］．Materials Evaluation，1999，9:587 －593．

［4］ Morozov M，Tian Guiyun，Withers P J．The pulsed eddy current response to applied loading of various aluminium alloys［J］．NDT＆E International，2010，43:493 －500．

［5］ Fava J，Ruch M．Design，construction and characterization of ECT sensors with tectangular planar coils［J］．Insight，2004，46（5）:268－274．

［6］ Greenhouse H M．Design of planar rectangular microelectronic inductors［J］．IEEE T Parts Hyb Pac（PHP），1974，10（2）:101 －109．

［7］ Theodoulids T P，Kriezis E E．Impedance evaluation of rectangular coils for eddy current testing of planar media［J］．NDT＆E International，2002，35:407 － 414．

［8］ He Yunze，Luo Feilu，Pan Mengchun．Defect edge identification with rectangular pulsed eddy current sensor based on transient response signals［J］．NDT＆E International，2010，43:409 －415．

［9］ He Yunze，Luo Feilu，Pan Mengchun．Defect classification based on rectangular pulsed eddy current sensor in different directions［J］．Sensors and Actuators A:Physical，2010，157:26 － 31．

［10］ He Yunze，Luo Feilu，Pan Mengchun．Defect characterization based on pulsed eddy current imaging technique［J］．Sensors and Actuators A:Physical，2010，164:1 － 7．