丁 華，何宇廷，杜金強，張海威，伍黎明

(空軍工程大學工程學院，西安710038)

對傳統渦流傳感器的陣列化構造技術作為一個新興分支已逐漸成為渦流無損檢測技術和結構健康監(jiān)測技術中的研究重點。陣列化構造技術不僅能夠實現對傳統渦流檢測技術在檢測空間上的微觀拓展，提高其檢測靈敏度，而且也能夠實現在檢測空間上的宏觀拓展，提高其檢測速度。目前，渦流陣列傳感器已逐漸應用于蒸汽管道廣布微裂紋、核反應堆殼體裂紋、高溫部件表面裂紋、管棒條型材、大面積金屬曲面裂紋的無損檢測和在線監(jiān)測等領域[1－8]。

在渦流陣列傳感器技術研究中，傳感器正向模型和逆向模型的建立是非常重要的部分，其是傳感器結構設計和參數優(yōu)化的基礎。同國外相比，國內目前對渦流陣列傳感器技術中的模型研究并不多，同時建立的陣列物理模型還只是將傳統單線圈物理模型進行簡單的合并，正向模型和逆向模型還只是側重于通過有限元方法進行純數值分析[9－11]。而國外研究者在陣列傳感器結構物理模型和正向、逆向解析或數值模型研究等方面進行了很多的理論和技術創(chuàng)新[12－13]，其中又以 MWM(meandering winding magnetometer)為代表[14－15]。在結構設計上，MWM是一種柔性平面渦流陣列傳感器，由單個激勵線圈和多個感應單元構成，感應單元可以是傳統的線圈也可以是GMR(Giant magnetoresistive)傳感元;在正向、逆向模型研究上，MWM的平面化設計使研究的電磁場量空間從三維空間降低到二維平面空間，解析模型構建的難度大大降低，而同數值法相比，解析法在逆向模型的構建上具有相當大的優(yōu)勢。

在Yanko Sheiretov[16]的研究基礎之上，本文針對MWM技術體系的一類分布式激勵源渦流陣列傳感器的半解析模型展開研究，提出了該傳感器半解析模型的直接FFT(Fast Fourier Transform)構建方法。該構建方法簡化了該分布式激勵源渦流陣列傳感器半解析模型的構建過程，同時也提高了建模精度，為后續(xù)傳感器的結構改進和參數優(yōu)化設計奠定基礎。

1 分布式激勵源渦流陣列傳感器

分布式激勵源渦流陣列傳感器的概念首先由James R.Melcher提出并經 N.J.Goldfine[17]等人得到發(fā)展。圖1和圖2為分布式激勵源渦流陣列傳感器的幾何結構示意圖和激勵源電流大小分布圖[16]。

圖1 分布式激勵源渦流陣列傳感器的幾何結構示意圖

圖2 激勵源大小和方向分布圖

圖1中，傳感器由在y向上等間距隔開的多直導線組成，相鄰導線間距為h，在y向總寬度為L。在傳感器導線中通激勵電流，激勵電流的大小和方向如圖2所示，即激勵電流滿足正弦包線分布，Iz(y)表示電流沿z向流動，大小和方向隨y向變化。圖1中，由于傳感器導線厚度較小，將傳感器簡化成平面，傳感器下方為空氣，上方為被測金屬結構，其磁導率為μ，電導率為σ。激勵源正弦包線分布設計的目的在于提高激勵源所產生的激勵場的主波長，從而減低激勵場向外擴散的衰減速率，提高傳感器的檢測深度。由Maxwell微分方程組同時不考慮位移電流因素(激勵頻率小于10 MHz)可推得磁矢A的空間擴散方程[18]

在圖1中，導線電流沿z向流動，則空間中磁矢只存在z分量，通過分離變量法，可得在各向同性層狀介質空間內磁矢A的解析通解[19]

其中 ky=2nπ/λ，n=1，2，…，Az0，n和 AzΔ，n為層狀空間上下界面處磁矢A的傅里葉變換系數，λ為激勵場空間分布主波長，Δ為介質層的厚度。由式(2)可知，主波長λ越大則激勵場在場量空間上的衰減速率就越慢。

圖3是寬度為L的分布式激勵源渦流陣列傳感器在x=0界面處的激勵電流正弦變換的頻譜圖，q是波數。由圖3可知，該類型的傳感器具有0.837L的主波長，即采取圖2所示的分布式激勵源可在空間得到與激勵源幾何分布尺寸相近的激勵場空間分布主波長。同其它類型的陣列傳感器相比，該類型的傳感器在測量深度裂紋上具有較大優(yōu)勢[16]。

圖3 寬度為L的分布式激勵源正弦變換頻譜圖

2 傳感器半解析模型的直接FFT構建

針對分布式激勵源渦流陣列傳感器的半解析模型構建問題，Yanko Sheiretov首先通過卷積運算得到傳感器在x=0界面處的線電流密度正弦變換解析解，通過離散化得到線電流密度的傅里葉變換系數離散值;然后根據層狀介質空間上下界面處場量的傅里葉變換系數傳遞關系矩陣，并使用洛倫茲因子來等效激勵源導線寬度對傅里葉變換系數的影響得到x=0界面處場量的傅里葉變換系數;最后通過反FFT變換得到x=0界面處場量的空間分布圖。在整個半解析建模過程中，不僅建模過程較為復雜，而且通過分析模型結果，發(fā)現其使用洛倫茲因子來等效激勵源導線寬度對模型影響的做法導致激勵源導線內部的線電流密度分布與實際情況不符。本文通過直接對傳感器在x=0界面處線電流密度進行虛擬采樣，而后通過FFT和反FFT得到x=0界面處場量的空間分布，不僅簡化了半解析模型的建模過程而且也將導線寬度直接考慮在模型之中。

2.1 傳感器x=0界面處線電流密度虛擬采樣及其FFT

設分布式激勵源渦流陣列傳感器的寬度為L，第m個激勵源導線中心坐標為ym，h=L/16，激勵源所滿足的正弦包線幅值為Id，激勵源導線寬度為W。第m個激勵源導線中的線電流密度分布滿足

即在傳感器寬度范圍內，線電流密度滿足式(4)，在傳感器寬度范圍之外和激勵源間隔內的線電流密度為0。根據線電流密度的分布規(guī)律通過離散化和數值方法可以得到傳感器在x=0界面處線電流密度的傅里葉變換系數 Kz，n(x=0)。

2.2 層狀介質空間中場量傅里葉變換系數傳遞關系矩陣

考慮如圖4所示的層狀介質空間，其中AzΔ，n和Az0，n分別是層狀介質空間上下界面處磁矢Az的傅里葉變換系數，HzΔ，n和 Hz0，n分別是上下界面處 y 向磁場強度 Hy的傅里葉變換系數，BxΔ，n和分別 Bx0，n是上下界面處x向磁感應強度Bx的傅里葉變換系數，Δ為層狀介質空間厚度。

圖4 層狀介質空間示意圖

由式(2)并根據B=▽×A，可得

則層狀介質空間中上下界面處磁矢Az和y向磁場強度Hy、x向磁感應強度Bx的傅里葉變換系數傳遞關系矩陣為

2.3 磁感應強度Bx的反FFT求解

考慮一簡化情況，即將分布式激勵源渦流陣列傳感器置于空氣中(假定空氣中的電磁參數與真空相同)，針對傳感器上方與下方的空氣層，其介質層厚度Δ為無窮大，式(6)簡化為

式中，μa是真空磁導率。在x=0界面處的電磁場量滿足邊界條件

綜合上述式(8)、式(9)和式(10)，x=0界面處的磁感應強度的傅里葉變換系數Bx，n(x=0)與該界面處線電流密度的傅里葉變換系數Kz，n(x=0)存在如下關系

則由 2.1 節(jié)中 Kz，n(x=0)可求得 Bx，n(x=0)，通過反FFT變換即可得到x=0界面處磁感應強度Bx(x=0)的空間分布，同時根據磁場強度和磁感應強度之間的本構關系也可得到磁場強度Hx(x=0)的空間分布。

2.4 模型結果分析

設分布式激勵源渦流陣列傳感器的寬度L=0.05 m，相鄰激勵源導線間距h=L/16，正弦包線幅值Id=1 A，激勵源寬度W=L/256，采樣點數N=1 024×32，將傳感器置于空氣中。圖5為經反FFT后得到的磁場強度Hx(x=0)空間分布圖。

圖5 Hx(x=0)分布圖

為驗證建立的半解析模型準確性，同時考慮將分布式激勵源渦流陣列傳感器放置于空氣中的簡化情況，本文使用經典電磁理論計算傳感器中心處磁感應強度大小并與半解析模型分析結果相比較。

傳感器中激勵源導線z向長度相對于傳感器總寬度L在數量級上要大很多，可以將激勵源導線簡化為無限長直導線。根據安培定律，無限長直導線電流I在真空中所產生的磁感應強度滿足[20]

式中，r是空間場強點到直導線的垂直距離。對于分布式激勵源渦流陣列傳感器，激勵源寬度W同傳感器寬度L相比可忽略不計，則其傳感器中心處的磁感應強度Bx(x=0，y=0)為

在圖5中，傳感器中心處的磁感應強度為168.1μa，兩者幾乎相等，驗證了本文提出的半解析模型的準確性。

在文獻[17]中，Yanko Sheiretov在半解析模型推導過程中使用洛倫茲因子來等效激勵源導線寬度。圖6為使用洛倫茲因子后重新反FFT變換后得到的某個激勵源在其寬度范圍之內的線電流密度分布。而實際上，在低頻激勵下，導線內部線電流密度應該是均勻分布，即使在高頻激勵所導致的集服效應的影響下，導線內部的線電流密度也應該是中間小，兩邊大。

圖6 激勵源線電流密度分布

3 模型仿真應用

在渦流陣列傳感器應用領域，傳感器與被測試件之間提離距離和被測試件電導率的變化可以用來表征被測試件結構微觀完整性的變化情況，如被測試件表面點蝕的有無、深淺或是其表面微裂紋的萌生、擴展等。針對本文所研究的分布式激勵源渦流陣列傳感器，圖7為其一種簡單的應用物理模型，以Bx(x=0，y=0)為特征量，通過本文建立的半解析模型，建立提離距離、被測試件電導率同Bx(x=0，y=0)特征量的相互關系。

圖7中，Δ1為提離距離，Δ2為被測試件厚度，Δ3為空氣層厚度。

圖7 分布式激勵源渦流陣列傳感器應用物理模型

設置傳感器外形參數同2.2節(jié)，被測試件厚度Δ2為0.002 m，空氣層厚度Δ3為無窮大，激勵頻率為1 MHz。根據建立的半解析模型，圖8為Hx(x=0，y=0)隨Δ1和σ變化的三維圖。

圖8 Hx(x=0，y=0)隨提離距離和電導率變化的三維圖

分析圖8結果可知，相比于將傳感器置于空氣中，被測試件上的渦流場對空間磁場的擾動作用很大，不僅大大降低了Hx(x=0，y=0)的幅值而且還會改變其相位。本節(jié)所述的應用模型只考慮Hx(x=0，y=0)處的單特征量，對應與單傳感元的設計，通過在z向上放置多個磁場傳感元(GMR傳感元等)可以很容易實現陣列化設計。

4 結論

本文針對MWM技術體系的一類分布式激勵源渦流陣列傳感器的半解析模型展開研究，提出了傳感器半解析模型的直接FFT構建方法。

(1)分布式激勵源渦流陣列傳感器能夠在空間中激勵出主波長同其幾何分布尺寸相近的激勵場。同其它傳感器陣列化構造化技術相比，其主波長要大的多，從而降低激勵場在場量空間中的衰減速率，提高了傳感器檢測深度裂紋的能力。

(2)在Yanko Sheiretov基礎之上，提出了分布式激勵源渦流陣列傳感器半解析模型的直接FFT構建方法。建模過程和數值計算結果表明運用此方法不僅簡化了半解析模型的建立過程，而且直接將激勵源導線寬度考慮在建模過程之中，提高了建模精度。

(3)分布式激勵源渦流陣列傳感器對其周圍介質空間的電磁屬性，如電導率等較為敏感，結合該傳感器可以檢測深度裂紋的能力以及對周圍介質空間電磁屬性的敏感性，可望將該類型的傳感器運用于金屬結構的在線監(jiān)測之中。

[1]Renato Gracin.Advanced Eddy Currents System for Inspection of VVER Steam Generator Tubing[C]//7th international Conference on Nuclear Option in Countries with Small and Medium Electricity Grids，Dubrovnik，Croatia，2008.

[2]Mitsuru Uesaka，Kazumi Hakuta.Eddy-Current Testing by Flexible Microloop Magnetic Sensor Array[J].IEEE Trans on Magnetics，1998，34(4):2287－2297.

[3]Yanko Sheiretov，Dave Grundy.MWM-Array Sensors for In Situ Monitoring of High-Temperature Components in Power Plants[J].IEEE Sensors Journal，2009，9(11):1527－1536.

[4]Joubert P Y，Le Bihan Y.Multi Sensor Probe and Defect Classificationin Eddy Current tubing Inspection[J].Sensors and Actuators A:Physical，2006，129(1－2):10－14.

[5]丁天懷，陳祥林.電渦流傳感器陣列測試技術[J].測試技術學報，2006，20(1):1－5.

[6]楊賓峰，羅飛路，曹雄恒，等.脈動渦流腐蝕成像陣列傳感器應用研究[J].傳感技術學報，2005，18(1):112－115.

[7]劉波，羅飛路，侯良潔.平板表層缺陷檢測渦流陣列傳感器的設計[J].傳感技術學報，2011，24(5):679－683.

[8]樓敏珠，郭韻.基于渦流陣列傳感器的高強度螺栓掃描檢測成像系統設計[J].無損檢測，2009，31(6):485－487.

[9]劉波，羅飛路.多激勵渦流陣列探頭阻抗分析[J].中國電機工程學報，2010，30(18):122－126.

[10]何永勃，邵雨果.陣列渦流傳感器互感信息三圍有限元仿真[J].中國民航大學學報，2010，28(1):29－32.

[11]趙冬明.陣列傳感器下金屬板中缺陷擾動的磁場量分布研究[D].長沙:國防科學技術大學，2004.

[12]Zaoui A，Feliachi M，Doirat V，et al.A Fast 3D Modeling of Arrayed Eddy Current Sensor[C]//12th Biennial IEEE Conference on Electromagnetic Fields Computation，Miami，USA，2006.

[13]Z Tae-Ok Kima，Gil-Seung Leea.Modeling of Eddy Current Sensor Using Geometric and Electromagnetic Data[C]//Journal of Mechanical Science and Technology，2007，21:465－475.

[14]Vladimir Zilberstein，DarrellSchlicker.MWM EddyCurrent Sensors for Monitoring of Crack Initiation and Growth During Fatigue Tests and in Service[J].International Journal of Fatigue，2001，23:S477－S485.

[15]Vladimir Zilberstein，Karen Walrath.MWM Eddy-Current Arrays for Crack Initiation and Growth Monitoring[J].International Journal of Fatigue，2003，25:S1147－S1155.

[16]Yanko Sheiretov.Deep Penetration Magnetoquasistatic Sensors[D].Cambridge，MA，USA:Massachusetts Institute of Technology，2001.

[17]Goldfine N J，Washabaugh A，Schlicker D E.High Resolution Inductive Array Imaging of Buried Objects[C]//SPIE AeroSense 2000，Orlando，Fl，2000.

[18]范孟豹.多層導電結構電渦流檢測的解析建模研究[D].杭州:浙江大學，2009.

[19]James R Melcher.Continuum Electromechanics[M].Cambridge，MA:MIT Press，1981:2.32－2.33

[20]毛均杰，何建國.電磁場理論[M].長沙:國防科學技術大學出版社，1998:61.