(1.武漢理工大學(xué)船舶動(dòng)力工程技術(shù)交通行業(yè)重點(diǎn)實(shí)驗(yàn)室 湖北武漢430063；2.國(guó)家水運(yùn)安全工程技術(shù)研究中心可靠性工程研究所 湖北武漢430063)

潤(rùn)滑油中的金屬磨粒主要為摩擦副摩擦磨損過(guò)程中的副產(chǎn)物，包含了豐富的關(guān)于機(jī)械設(shè)備磨損部位、磨損程度以及磨損類(lèi)型的信息，對(duì)其進(jìn)行監(jiān)測(cè)對(duì)機(jī)械設(shè)備狀態(tài)的動(dòng)態(tài)監(jiān)測(cè)、設(shè)備的故障預(yù)知及延長(zhǎng)機(jī)械設(shè)備的使用壽命都具有重要意義[1]。

電感式磨粒監(jiān)測(cè)技術(shù)作為磨粒監(jiān)測(cè)的重要技術(shù)手段，在金屬磨粒監(jiān)測(cè)領(lǐng)域展現(xiàn)出了廣闊的應(yīng)用前景。圖1示出了電感式磨粒傳感器的工作原理，感應(yīng)線(xiàn)圈的電感變化主要由磨粒的磁化作用和渦流作用決定。當(dāng)感應(yīng)線(xiàn)圈中存在非鐵磁性磨粒時(shí)(如圖1(a)所示)，磨粒產(chǎn)生的渦流作用會(huì)對(duì)原磁場(chǎng)產(chǎn)生削弱作用，導(dǎo)致線(xiàn)圈的等效電感減小。當(dāng)感應(yīng)線(xiàn)圈中存在鐵磁性金屬磨粒時(shí)(如圖1(b)所示)，磁化作用和渦流作用共同作用，由于感應(yīng)線(xiàn)圈的激勵(lì)頻率較低，此時(shí)渦流作用很小，因此總磁通量主要由磁化作用所決定，感應(yīng)線(xiàn)圈的等效電感增加[2]。根據(jù)電感原理開(kāi)發(fā)的傳感器中最具代表性的是GasTOPS公司研制的MetalSCAN傳感器，其已經(jīng)在民航客機(jī)和風(fēng)力發(fā)電機(jī)中得到廣泛應(yīng)用，并且可以全流量識(shí)別等效直徑為70 μm的鐵磁性磨粒以及等效直徑為270 μm的非鐵磁性磨粒[3]。但該傳感器在實(shí)際使用過(guò)程中仍有兩點(diǎn)不足：①該傳感器采用的三螺管線(xiàn)圈結(jié)構(gòu)，在實(shí)際應(yīng)用過(guò)程中不能有效識(shí)別大量連續(xù)性流動(dòng)的磨粒；②該傳感器利用渦流作用和磁化作用，可實(shí)現(xiàn)金屬磨粒磁性的識(shí)別，但不能識(shí)別金屬磨粒的材質(zhì)。

圖1 電感式磨粒傳感器工作原理Fig 1 Working principle of inductive abrasive sensor

為了彌補(bǔ)MetalSCAN傳感器的不足，范紅波等[4]、LI和JIANG[5]在三螺管式磨粒傳感器的基礎(chǔ)上，通過(guò)縮短感應(yīng)線(xiàn)圈的長(zhǎng)度，分別開(kāi)發(fā)了雙螺管式和單螺管式磨粒傳感器。這2種傳感器雖大大降低了磨粒連續(xù)性對(duì)傳感器的影響，但仍不能彌補(bǔ)電感式傳感器無(wú)法識(shí)別磨粒材質(zhì)的不足。而通過(guò)增大感應(yīng)線(xiàn)圈的激勵(lì)頻率，提升渦流作用在磨粒識(shí)別中的地位，不僅可以縮短感應(yīng)線(xiàn)圈的長(zhǎng)度，而且還可以利用不同材質(zhì)磨粒表現(xiàn)出的不同渦流特性實(shí)現(xiàn)磨粒材質(zhì)識(shí)別的目的[6]?；诖?，本文作者在電感式磨粒傳感器的基礎(chǔ)上，通過(guò)放大渦流作用在磨粒識(shí)別中的地位，提出利用電渦流原理監(jiān)測(cè)潤(rùn)滑油中金屬磨粒的方法；同時(shí)結(jié)合ANSYS Maxwell仿真軟件從理論上驗(yàn)證該方法的可行性，并對(duì)電渦流磨粒傳感器的主要參數(shù)進(jìn)行了仿真分析，為今后傳感器的設(shè)計(jì)和優(yōu)化提供了理論指導(dǎo)。

1 電渦流磨粒傳感器的基本原理

電渦流磨粒傳感器的工作原理如圖2所示，當(dāng)感應(yīng)線(xiàn)圈通以高頻電流I1時(shí)，根據(jù)法拉第電磁感應(yīng)定律可知，感應(yīng)線(xiàn)圈周?chē)鷷?huì)產(chǎn)生交變磁場(chǎng)H1，若感應(yīng)線(xiàn)圈檢測(cè)區(qū)域內(nèi)無(wú)金屬磨?？拷鼤r(shí)，線(xiàn)圈產(chǎn)生的磁場(chǎng)能量全部損失；當(dāng)感應(yīng)線(xiàn)圈檢測(cè)區(qū)域內(nèi)存在金屬磨粒時(shí)，由于磨粒在圓周方向可以等效為一圈圈的閉合電路，閉合電路中磁通量不斷變化，則在磨粒表面產(chǎn)生電渦流I2，磨粒的渦流場(chǎng)也會(huì)產(chǎn)生一個(gè)交變磁場(chǎng)H2，H2與H1方向相反，從而導(dǎo)致感應(yīng)線(xiàn)圈的磁場(chǎng)發(fā)生變化。從能量損耗的角度來(lái)看，金屬磨粒流經(jīng)感應(yīng)線(xiàn)圈后，磨粒內(nèi)存在著渦流損耗與磁損耗，能量損耗會(huì)使感應(yīng)線(xiàn)圈的電動(dòng)勢(shì)降低[7]。

圖2 渦流檢測(cè)原理Fig 2 Principle of eddy current test

從圖2所示的渦流檢測(cè)原理可知，影響電渦流磨粒傳感器輸出的因素主要來(lái)源于金屬磨粒及感應(yīng)線(xiàn)圈2個(gè)方面。結(jié)合電感式磨粒傳感器的研究可知，金屬磨粒方面主要包括磨粒材質(zhì)、磨粒尺寸及磨粒形貌等；感應(yīng)線(xiàn)圈方面主要包括線(xiàn)圈匝數(shù)、線(xiàn)圈內(nèi)徑、線(xiàn)圈激勵(lì)頻率及線(xiàn)圈纏繞方式等[8]。為驗(yàn)證電渦流原理應(yīng)用于磨粒監(jiān)測(cè)的可行性，需要對(duì)不同材質(zhì)和尺寸的磨粒進(jìn)行仿真分析；另外，線(xiàn)圈激勵(lì)頻率和線(xiàn)圈內(nèi)徑是傳感器最重要的2個(gè)參數(shù)，對(duì)其進(jìn)行分析可為傳感器設(shè)計(jì)提供理論指導(dǎo)。下面利用ANSYS Maxwell有限元仿真軟件分別進(jìn)行仿真分析。

2 傳感器磁場(chǎng)仿真

2.1 仿真過(guò)程

ANSYS Maxwell是一款專(zhuān)門(mén)用于求解電磁場(chǎng)問(wèn)題的有限元仿真軟件，下面根據(jù)ANSYS Maxwell的仿真流程，簡(jiǎn)要介紹一下文中所使用的仿真模型及求解方法：

(1)選擇求解器：文中在電感式磨粒傳感器的基礎(chǔ)上，通過(guò)放大渦流作用在磨粒識(shí)別中的地位，提出了利用電渦流原理監(jiān)測(cè)潤(rùn)滑油中金屬磨粒的方法，這會(huì)用到磁場(chǎng)模塊中的渦流場(chǎng)分析，因此選擇渦流場(chǎng)求解器。

(2)建模：參考文獻(xiàn)[8]中電感式磨粒傳感器的建模方法，將電渦流金屬磨粒傳感器簡(jiǎn)化為金屬磨粒、感應(yīng)線(xiàn)圈和真空求解域，建立如圖3所示的三維仿真模型。

(3)材料屬性設(shè)置：計(jì)算過(guò)程中要用到的材料主要有：銅(感應(yīng)線(xiàn)圈)，真空(求解域)，硅、銅、鐵、鎳(磨粒)，其材料屬性均采用材料庫(kù)中的默認(rèn)值。

(4)網(wǎng)格劃分：結(jié)合模型的實(shí)際情況，采用自適應(yīng)網(wǎng)格劃分。

(5)激勵(lì)源設(shè)置：給感應(yīng)線(xiàn)圈添加一個(gè)大小為1 A，相位為0°，類(lèi)型為Stranded的電流源激勵(lì)。

(6)設(shè)定求解參數(shù)：添加渦流場(chǎng)求解設(shè)置，設(shè)置最大收斂步數(shù)為20，其余參數(shù)采用軟件的默認(rèn)值。

(7)后處理：借鑒文獻(xiàn)[10]中用渦流密度的值表征渦流作用大小的方法，將磨粒的渦流密度云圖輸出，通過(guò)渦流密度的分析研究各參數(shù)對(duì)傳感器輸出的影響。

圖3 仿真模型圖Fig 3 Simulation model

2.2 模型可靠性分析

根據(jù)畢奧-薩伐爾定律可知，圓形線(xiàn)圈中軸線(xiàn)的磁場(chǎng)分布規(guī)律[11]為

式中：B為磁場(chǎng)強(qiáng)度；μ0為真空磁導(dǎo)率；I為線(xiàn)圈中的電流；X為距離線(xiàn)圈中心的距離；R為圓形線(xiàn)圈內(nèi)徑。

由上式可知，圓形線(xiàn)圈中軸線(xiàn)處磁場(chǎng)強(qiáng)度沿線(xiàn)圈中心呈對(duì)稱(chēng)分布，最大磁場(chǎng)強(qiáng)度為線(xiàn)圈中心處，而離線(xiàn)圈中心越遠(yuǎn)的地方，其磁場(chǎng)強(qiáng)度越小。對(duì)比分析圖4，可知磨粒產(chǎn)生的渦流密度沿線(xiàn)圈軸線(xiàn)方向上出現(xiàn)分層分布現(xiàn)象，即渦流密度在線(xiàn)圈中心最大，向線(xiàn)圈兩邊逐漸減小，這與圓形線(xiàn)圈中軸線(xiàn)上的磁場(chǎng)分布規(guī)律相吻合。

將圖4中金屬磨粒沿A-A截面進(jìn)行剖分，將剖面的電渦流密度云圖輸出，如圖5所示?？芍?，渦流密度在磨粒內(nèi)部分布不均勻，在磨粒的外表面一側(cè)最密集，隨著深度的增加，渦流密度逐漸衰減。由電流的集膚效應(yīng)[12]可知，當(dāng)導(dǎo)體通以交變電流時(shí)，電流將會(huì)聚集于導(dǎo)體表層，而非平均分布于整個(gè)導(dǎo)體的截面積中，這與圖5所示的規(guī)律相吻合。

圖4 磨粒外表面渦流密度云圖Fig 4 Eddy density cloud map on the outer surface of debris

圖5 磨粒剖面渦流密度云圖Fig 5 Eddy density cloud map in the profile of debris

綜合磁場(chǎng)分布規(guī)律及集膚效應(yīng)的分析可知，該模型的仿真結(jié)果符合客觀物理規(guī)律，下面采用上述模型和計(jì)算方法分別對(duì)不同尺寸和材質(zhì)的磨粒進(jìn)行分析，并探討不同激勵(lì)頻率及不同線(xiàn)圈內(nèi)徑時(shí)對(duì)分析結(jié)果的影響。

2.3 仿真結(jié)果及分析

由于渦流密度云圖難以實(shí)現(xiàn)數(shù)值化分析，所以文中利用ANSYS Maxwell自帶的“Marker”功能，在磨粒外表面靠近線(xiàn)圈中心處(即圖5中磨粒的最外圈)拾取100個(gè)點(diǎn)，將這100個(gè)點(diǎn)的渦流密度平均值作為表征磨粒電渦流作用強(qiáng)弱的特征參數(shù)，利用該參數(shù)得到不同磨粒尺寸、不同磨粒材質(zhì)以及不同激勵(lì)頻率及不同線(xiàn)圈內(nèi)徑下的仿真結(jié)果，如圖6所示，各仿真參數(shù)設(shè)置如表1所示。

圖6 不同參數(shù)下仿真結(jié)果Fig 6 Simulation results at different pqrameters(a) at different debris radius;(b)at different debris materials;(c)at different excitation frequency;(d)at different diameter of inner coils表1 仿真參數(shù)表Table 1 Simulation parameter Table

仿真參數(shù)abcd磨粒半徑r/μm變量200200200磨粒材質(zhì)銅變量銅銅激勵(lì)頻率f/MHz280280變量280線(xiàn)圈內(nèi)徑d/μm1 6001 6001 600變量

2.3.1 不同磨粒尺寸的仿真結(jié)果分析

馮炳華和杜永平[13]、吳超等人[14]分別從理論推導(dǎo)和仿真分析的角度，得到了電感式傳感器的電感變化量與油液磨粒尺寸成3次方關(guān)系的結(jié)論。結(jié)合圖6(a)可知，磨粒尺寸越大，磨粒產(chǎn)生的渦流密度越大，曲線(xiàn)的斜率也逐漸增加，通過(guò)對(duì)數(shù)據(jù)擬合發(fā)現(xiàn)磨粒尺寸與渦流密度成3次方的關(guān)系。因此與電感式磨粒傳感器相同，電渦流磨粒傳感器的輸出與磨粒尺寸也是立方的關(guān)系。

2.3.2 不同磨粒材質(zhì)的仿真結(jié)果分析

結(jié)合圖6(b)可知，相同的磁場(chǎng)環(huán)境中不同材質(zhì)的磨粒產(chǎn)生的渦流密度不同：非金屬磨粒(硅)不產(chǎn)生電渦流作用，鐵磁性磨粒鐵、鎳產(chǎn)生的電渦流作用大于非鐵磁性磨粒銅，鐵磁性磨粒中鎳的電渦流作用大于鐵的電渦流作用。

產(chǎn)生上述結(jié)果的原因是在非導(dǎo)體內(nèi)，不會(huì)出現(xiàn)電與磁之間的感應(yīng)與轉(zhuǎn)化作用，因此非導(dǎo)體不產(chǎn)生電渦流作用，而金屬導(dǎo)體由于其具有不同的電導(dǎo)率及磁導(dǎo)率，因此產(chǎn)生的電渦流作用也不相同，所以可以根據(jù)磨粒產(chǎn)生的電渦流大小識(shí)別判斷磨粒的材質(zhì)。

2.3.3 不同激勵(lì)頻率的仿真結(jié)果分析

線(xiàn)圈的激勵(lì)頻率是決定鐵磁性磨粒識(shí)別時(shí)磁化作用和渦流作用貢獻(xiàn)率的決定性參數(shù)，文獻(xiàn)[13]表明，電感式磨粒傳感器的激勵(lì)頻率一般為1～2.5 MHz，當(dāng)激勵(lì)頻率足夠高時(shí)，渦流作用將遠(yuǎn)遠(yuǎn)大于磨粒的磁化作用。結(jié)合圖6(c)可知，低頻激勵(lì)時(shí)，磨粒的渦流作用較小，隨著激勵(lì)頻率的提高，磨粒的渦流作用越大，曲線(xiàn)的斜率越小。

為減小磨粒磁化作用對(duì)鐵磁性金屬磨粒識(shí)別的影響，應(yīng)盡量選用較高的線(xiàn)圈激勵(lì)頻率。但線(xiàn)圈的激勵(lì)頻率過(guò)高，磨粒信號(hào)中將引入大量的高頻噪聲，并且由于電氣元件的限制，激勵(lì)頻率也不可能無(wú)限大。因此，激勵(lì)頻率的選取應(yīng)綜合考慮以上兩方面的因素。

2.3.4 不同線(xiàn)圈內(nèi)徑的仿真結(jié)果分析

線(xiàn)圈內(nèi)徑是決定電渦流傳感器流量和可靠性的關(guān)鍵因素，結(jié)合圖6(d)可知，采用電渦流原理監(jiān)測(cè)磨粒時(shí)，線(xiàn)圈的內(nèi)徑越大，磨粒產(chǎn)生的渦流密度越小，曲線(xiàn)的斜率越小。因此在設(shè)計(jì)傳感器時(shí)，應(yīng)盡量減小線(xiàn)圈內(nèi)徑，從而提高傳感器的靈敏度。同時(shí)還應(yīng)該考慮傳感器的流量和所使用的環(huán)境，因?yàn)殡S著線(xiàn)圈內(nèi)徑的減小，傳感器的流量急劇減小，同時(shí)流道被大磨粒堵塞的風(fēng)險(xiǎn)也增加。

3 結(jié)論

(1)電渦流作用可以識(shí)別磨粒尺寸，磨粒產(chǎn)生的渦流作用與磨粒尺寸成3次方關(guān)系。

(2)電渦流作用可以識(shí)別磨粒材質(zhì)，不同材質(zhì)的磨粒在相同的磁場(chǎng)環(huán)境中產(chǎn)生的渦流作用不同：非金屬磨粒(硅)不產(chǎn)生電渦流作用，不同材質(zhì)的金屬磨粒具有不同的電導(dǎo)率及磁導(dǎo)率，所產(chǎn)生的渦流作用也不同，在文中仿真模型中渦流作用由大到小依次為鎳磨粒、鐵磨粒、銅磨粒。

(3)線(xiàn)圈的激勵(lì)頻率越高，磨粒的渦流作用越大。在設(shè)計(jì)傳感器時(shí)，應(yīng)綜合考慮磁化作用和高頻噪聲的影響，選取合適的激勵(lì)頻率，從而保證傳感器的靈敏度和信號(hào)的信噪比。

(4)線(xiàn)圈的內(nèi)徑越大，磨粒的渦流作用越小。因此，在保證傳感器流量正常和流道不堵塞的前提下，應(yīng)盡量減小感應(yīng)線(xiàn)圈的內(nèi)徑，從而提高傳感器檢測(cè)靈敏度。