齊 川 張思全

（上海海事大學(xué) 物流工程學(xué)院，中國 上海 201306）

0 引言

隨著工業(yè)技術(shù)的快速發(fā)展，許多行業(yè)對(duì)于金屬材料表面光滑度的要求越來越高，光滑度的高低直接影響到工業(yè)生產(chǎn)的安全、產(chǎn)品的質(zhì)量等，因而對(duì)金屬材料表面缺陷的檢測(cè)已經(jīng)成為工業(yè)生產(chǎn)中的一個(gè)重要環(huán)節(jié)。目前，針對(duì)金屬材料表面的缺陷，常用的無損檢測(cè)方法有射線檢測(cè)、超聲檢測(cè)、磁粉檢測(cè)、渦流檢測(cè)等。其中，渦流檢測(cè)技術(shù)是一種基于電磁感應(yīng)原理的常規(guī)無損檢測(cè)方法，它采用非接觸測(cè)量，不需要耦合劑，不需要對(duì)材料表面清洗，具有檢測(cè)速度快，靈敏度高，抗干擾能力強(qiáng)等特點(diǎn)，因此在工業(yè)生產(chǎn)和科學(xué)研究等各領(lǐng)域都得到了廣泛的應(yīng)用。[1]

金屬導(dǎo)體表面裂紋缺陷的定性、定量分析一直是渦流檢測(cè)中的重點(diǎn)問題。 本文基于有限元仿真法，利用Ansoft Maxwell 軟件建立包括線圈和導(dǎo)體平板的渦流檢測(cè)系統(tǒng)有限元模型，對(duì)不同深度的裂紋進(jìn)行有限元仿真，從電磁場分布的角度深入分析不同裂紋深度對(duì)磁感應(yīng)強(qiáng)度的影響，從而為快速、準(zhǔn)確地獲取缺陷的真實(shí)深度提供參考依據(jù)，同時(shí)也有利于渦流傳感器參數(shù)的優(yōu)化設(shè)計(jì)，并降低了渦流檢測(cè)系統(tǒng)的實(shí)驗(yàn)成本。[2]

1 基本原理

渦流檢測(cè)是建立在電磁感應(yīng)原理基礎(chǔ)之上的一種無損檢測(cè)方法，由三部分組成，即加載交變電流的檢測(cè)線圈、檢測(cè)電流的儀器和被檢測(cè)的導(dǎo)電物體。 當(dāng)通有交變電流的線圈靠近被檢測(cè)導(dǎo)電物體時(shí)，由于線圈中交變電流產(chǎn)生的一次磁場作用于導(dǎo)體，在導(dǎo)體表面感應(yīng)出渦電流，該渦電流同時(shí)產(chǎn)生一個(gè)與原磁場方向相反的二次磁場作用在檢測(cè)線圈。若導(dǎo)電物體的表面存在裂紋缺陷，由于裂紋電阻大，會(huì)切斷或降低電渦流，使其流向發(fā)生扭曲，朝著缺陷的底部和邊緣偏轉(zhuǎn)，從而影響到工件周圍感應(yīng)磁場的分布。因此，可通過測(cè)量該感應(yīng)磁場的變化，實(shí)現(xiàn)裂紋等缺陷的定性、定量分析。[4]

2 理論基礎(chǔ)

19 世紀(jì)中期，麥克斯韋在總結(jié)前人工作的基礎(chǔ)上，提出了適用于所有宏觀電磁現(xiàn)象的數(shù)學(xué)模型，稱之為麥克斯韋方程組。 麥克斯韋方程組描述了電磁場的基本特征，是渦流檢測(cè)技術(shù)的理論基礎(chǔ)。 它由四個(gè)定律組成，分別為：安培環(huán)路定律、高斯磁通定律、法拉第電磁感應(yīng)定律和高斯電通定律。[3]其微分形式如下所示：

式中E 表示電場強(qiáng)度矢量，H 表示磁場強(qiáng)度矢量，D 表示電通密度矢量，B 表示磁通密度矢量，J 表示電流密度矢量，ρ 表示電荷密度。上述各參量在均勻介質(zhì)中滿足如下關(guān)系：

其中，ε 為介電常數(shù)，μ 為磁導(dǎo)率，σ 為電導(dǎo)率。

3 有限元仿真

本文采用Ansoft Maxwell 軟件對(duì)檢測(cè)系統(tǒng)進(jìn)行電磁場有限元仿真分析， 其中Maxwell 3D 渦流場模塊用于分析導(dǎo)體中時(shí)變電流或外界交變磁場源所引起的時(shí)變磁場， 它利用自適應(yīng)分析法進(jìn)行網(wǎng)格剖分，可使求解問題的速度和精度得到提高。

3.1 模型構(gòu)建

首先建立渦流檢測(cè)模型。 模型主要由線圈、鐵芯和被測(cè)工件三部分組成。 三維模型如圖1 所示。 其中，被測(cè)工件是一塊（40*40*5）mm3的不銹鋼平板，其電導(dǎo)率為16.5MS·m-1，相對(duì)磁導(dǎo)率為1，表面正中心有一個(gè)矩形缺陷。 激勵(lì)線圈的內(nèi)徑為3.2mm，外徑為5.6mm，提離高度為1mm，匝數(shù)為800，材料屬性為系統(tǒng)默認(rèn)的銅。 鐵芯的底面半徑為4mm，高度為5mm，材料屬性為鐵氧體。 背景區(qū)域?yàn)榭諝狻?各部分材料的電導(dǎo)率、磁導(dǎo)率恒定且各向同性。

圖1 導(dǎo)電平板渦流檢測(cè)的有限元模型

3.2 參數(shù)設(shè)定與求解器設(shè)定

有限元的參數(shù)設(shè)定一般包括設(shè)置邊界條件、 激勵(lì)源和網(wǎng)格剖分等。 該模型中，由于對(duì)稱面兩側(cè)的電流大小相等，符號(hào)相反，而且磁場與對(duì)稱面相切， 所以將背景區(qū)域的左邊界和后邊界設(shè)為Zero Tangential H Field 邊界條件，即磁場強(qiáng)度切向分量恒為零邊界條件。

由于線圈為多股聚酯漆包銅線多層并繞而成，故本文采用絞線電流源strand，匝數(shù)為800。 假定線圈激勵(lì)電流在其截面上均勻分布。 激勵(lì)源電流初始相位為0。 繞組電流設(shè)為正弦激勵(lì)電流，幅值為100A，相位為0。

考慮到渦流對(duì)導(dǎo)體的集膚效應(yīng)，對(duì)被測(cè)導(dǎo)體采用基于集膚效應(yīng)滲透深度剖分設(shè)置，即導(dǎo)體表面的網(wǎng)格較為稠密，而集膚效應(yīng)層之下的網(wǎng)格則相對(duì)較為稀疏。輸入材料的相對(duì)磁導(dǎo)率、電導(dǎo)率和工作頻率，軟件便可自動(dòng)計(jì)算出集膚效應(yīng)的滲透深度。

最后一步為設(shè)置求解器。 設(shè)定激勵(lì)頻率為10kHz，選擇自適應(yīng)分析，設(shè)置計(jì)算的循環(huán)數(shù)50、每次循環(huán)改進(jìn)的百分比為30%，以及總的能量誤差5%。

4 仿真結(jié)果分析

根據(jù)Maxwell 電磁感應(yīng)理論，激勵(lì)線圈遠(yuǎn)離缺陷區(qū)域時(shí)，工件周圍的感應(yīng)磁場流動(dòng)受到的阻礙很小.磁感應(yīng)強(qiáng)度B 基本保持恒定； 遇到缺陷時(shí)，感應(yīng)磁場會(huì)發(fā)生偏轉(zhuǎn)，磁感應(yīng)強(qiáng)度會(huì)發(fā)生變化。 而不同的缺陷深度產(chǎn)生的磁感應(yīng)強(qiáng)度變化也是不同的。[5-6]為了研究裂紋深度對(duì)測(cè)量結(jié)果的影響，使裂紋的長度和寬度均為定值，將裂紋的深度依次設(shè)為1mm、2mm、3mm。在三種不同裂紋深度的情況下分別進(jìn)行仿真。仿真時(shí)，裂紋與X 軸平行，激勵(lì)線圈沿Y 軸方向水平移動(dòng)，垂直于缺陷走向；電流沿X 軸方向，平行于缺陷走向。 線圈從Y 軸-16mm 處開始掃描，每移動(dòng)2mm 作一次仿真，終點(diǎn)為Y 軸16mm，總共進(jìn)行15次仿真。 圖2 為線圈掃描到裂紋正上方時(shí)導(dǎo)電平板磁場密度分布圖。

圖2 磁場密度分布圖

圖3 磁感應(yīng)強(qiáng)度變化曲線

以激勵(lì)線圈中心正下方0.5mm的A 點(diǎn)作為參照點(diǎn)， 計(jì)算A 點(diǎn)在不同裂紋深度下的每一次仿真的磁感應(yīng)強(qiáng)度。圖5 為裂紋深度變化時(shí)A 點(diǎn)磁感應(yīng)強(qiáng)度B 的變化曲線。

從圖中可以看出， 當(dāng)線圈離裂紋較遠(yuǎn)時(shí)，A 點(diǎn)的磁感應(yīng)強(qiáng)度變化較為平穩(wěn)，基本在6.1 處作小范圍的波動(dòng)，且裂紋深度對(duì)磁感應(yīng)強(qiáng)度的影響不大； 隨著線圈向裂紋靠近，A 點(diǎn)的磁感應(yīng)強(qiáng)度的變化幅度逐漸變大， 當(dāng)A 點(diǎn)正好位于裂紋正上方時(shí)， 磁感應(yīng)強(qiáng)度的波動(dòng)達(dá)到最大，且裂紋越深，磁感應(yīng)強(qiáng)度變化越大。

5 結(jié)束語

本文基于有限元分析方法對(duì)渦流傳感器檢測(cè)系統(tǒng)進(jìn)行了數(shù)值仿真， 利用Ansoft Maxwell 軟件建立了渦流傳感器的檢測(cè)三維模型，對(duì)于不同深度的裂紋，分別使激勵(lì)線圈在導(dǎo)體上方沿裂紋的垂直方向進(jìn)行掃描并仿真， 得到了線圈磁力線分布和被測(cè)導(dǎo)體的磁場密度分布，通過計(jì)算線圈正下方磁感應(yīng)強(qiáng)度，重點(diǎn)分析了不同裂紋深度對(duì)對(duì)磁感應(yīng)強(qiáng)度的影響。 仿真結(jié)果表明：裂紋缺陷會(huì)對(duì)線圈的磁感應(yīng)強(qiáng)度發(fā)生變化，且這一變化隨裂紋深度的增加而變大。

［1］徐可北，周華俊，任吉林.渦流檢測(cè)[M].北京：機(jī)械工業(yè)出版社，2011：21-25.

［2］趙博，張洪亮.Ansoft12 在工程電磁場中的應(yīng)用[M].北京：中國水利水電出版社，2010：19-23.

［3］謝亞楠.電磁場理論[M].北京：清華大學(xué)出版社，2011：38-42.

［4］羅成剛.電渦流位移傳感器線圈電磁場仿真分析[J].傳感器與微系統(tǒng)，2008，27(1)：24-26.

［5］周德強(qiáng),張斌強(qiáng),田貴云等.脈沖渦流檢測(cè)中裂紋的深度定量及分類識(shí)別[J].儀器儀表學(xué)報(bào)，2009,30

［6］陳德志，邵可然.管道裂紋渦流檢測(cè)線圈阻抗信號(hào)的快速仿真[J].電工技術(shù)學(xué)報(bào)，2000,12.

［7］Luquire J W,Deeds W E,Dodd C V.Alternating Current Distribution Between Planar Conductors[J].Journal of Applied Physics,1970,41(10)：3983-3991.

［8］Haitao Wang，Qiufeng Luo，Xin Wang，et al.Simulation and experimental study of remote field eddy current testing on flat conductive plate[J].International Journal of Applied Electromagnetics and Mechanics，2010，33(3-4)：1261-1266.

［9］Sharma C.The Effects of Lift-Off from Wall Thinning Signal in Pulsed Eddy Current Testing[J].Journal of Magnetics，2012，17(4)：298-301.