一種新型可視鐵譜儀的靜磁場設(shè)計及應(yīng)用

黃文杰， 左洪福(. 淮陰工學(xué)院 交通工程學(xué)院，江蘇 淮安 3003； . 南京航空航天大學(xué) 民航學(xué)院，江蘇 南京 006)

?

一種新型可視鐵譜儀的靜磁場設(shè)計及應(yīng)用

黃文杰1， 左洪福2
(1. 淮陰工學(xué)院 交通工程學(xué)院，江蘇 淮安 223003； 2. 南京航空航天大學(xué) 民航學(xué)院，江蘇 南京 210016)

以在線式全流量可視鐵譜儀為背景，對靜磁場的磁路進行了設(shè)計。采用軸向充磁，雙環(huán)異極嵌套安裝方式，兩環(huán)的間隙采用不導(dǎo)磁的材料做成的圓環(huán)作為襯套，并從磁荷模型出發(fā)，利用標量磁位法，引入廣義二項式定理對同心雙環(huán)嵌套磁體產(chǎn)生的磁感應(yīng)強度進行了解析計算。通過對比磁感應(yīng)強度的理論值和實際測量值發(fā)現(xiàn)解析表達式很好地反映了靜磁場的磁感應(yīng)強度分布規(guī)律。將所設(shè)計的靜磁場應(yīng)用到某新型渦噴航空發(fā)動機性能監(jiān)測中,對監(jiān)測到的磨粒進行了分析。分析結(jié)果間接驗證了所設(shè)計的靜磁場在潤滑油高流速條件下吸附鐵磁顆粒的有效性。

在線式全流量鐵譜儀； 靜磁場； 磁路設(shè)計； 解析計算； 性能監(jiān)測

0 引 言

由于航空發(fā)動機關(guān)鍵部件的摩擦副工作在潤滑油路中且在潤滑油路中利用磁場便于收集磨粒，因此，磁場在航空發(fā)動機的故障診斷和狀態(tài)監(jiān)測中有著廣泛的應(yīng)用。以滑油中的磨粒監(jiān)測為例，早期的監(jiān)測設(shè)備以離線為主，如DMAS智能化鐵譜分析儀[1]、OLF-4鐵譜儀[2]等。隨著監(jiān)測技術(shù)的發(fā)展，一批具有代表性在線全流量監(jiān)測傳感器應(yīng)運而生，其中應(yīng)用最廣泛的是加拿大GOSTOP公司生產(chǎn)的MetalSCAN油液磨損金屬監(jiān)測傳感器，它能夠全液流在線監(jiān)測鐵磁顆粒[3]。由于稀土永磁材料的加工水平與工藝不斷提高，合理設(shè)計的靜磁體可以產(chǎn)生各種強度和分布的磁場，并且具有無能耗、體積小、成本低、控制方便、工作可靠等優(yōu)點，尤其是環(huán)形磁體，因其便于在油路中安裝，易于產(chǎn)生均勻分布的磁場，在滑油磨粒監(jiān)測方面具有廣泛的應(yīng)用前景，而靜磁體應(yīng)用的關(guān)鍵是對靜磁體磁場分布的設(shè)計。在這方面，國內(nèi)外對磁場分布的研究不斷深入。在國內(nèi)，運用磁荷模型[4]、分子電流模型[5]給出了矩形靜磁場分布的解析表達式；在國外，在采用橢圓積分的基礎(chǔ)上[6-7]，或運用磁荷模型對弧形靜磁場進行研究[8-9]，或采用半解析方法計算環(huán)形靜磁場分布[10-11]，并采用該方法對使用永磁體的裝置進行優(yōu)化[12]。通過對以往研究成果的分析，不難發(fā)現(xiàn)，國內(nèi)外對矩形靜磁場分布的研究比較全面，雖然采用半解析方法計算了環(huán)形靜磁場分布，但是，缺少真實意義上的解析公式，嚴重制約了環(huán)形靜磁場的設(shè)計和分析。而環(huán)形靜磁場作為在線式全流量鐵譜儀的核心部件，其磁場的分布決定著儀器在實際工作中能否有效吸附鐵磁顆粒，從而能否準確對航空發(fā)動機的性能作出準確預(yù)測。因此，真實意義上的解析公式的缺失阻礙了環(huán)形靜磁場在航空發(fā)動機的故障診斷和狀態(tài)監(jiān)測中使用和推廣。

本文以在線式全流量可視鐵譜儀為背景，引入廣義二項式定理對靜磁場的磁路進行設(shè)計，并對同心雙環(huán)嵌套磁體的磁場分布進行解析計算，通過對比磁感應(yīng)強度的理論值和實際測量值，發(fā)現(xiàn)解析表達式很好地反映了靜磁場的磁感應(yīng)強度分布規(guī)律。在此基礎(chǔ)上，將所設(shè)計的靜磁場應(yīng)用到某新型渦噴航空發(fā)動機性能監(jiān)測中，對監(jiān)測到的磨粒進行了分析。

1 靜磁場磁路設(shè)計

磁體采用牌號為N35燒結(jié)釹鐵硼永磁材料，其性能參數(shù)，剩磁Br=1.170～1.220 T，矯頑力HCB≥868 kA/m，內(nèi)稟矯頑力HCJ≥955 kA/m，磁能積BHmax=263～287 kJ/m3。采用軸向充磁、雙環(huán)異極嵌套安裝的方式，兩環(huán)間隙用不導(dǎo)磁材料做成的圓環(huán)作為襯套。

由矢量疊加原理，該磁體在空間任一點的磁感應(yīng)強度可由內(nèi)圓環(huán)在該點的磁感應(yīng)強度與外圓環(huán)在同一點的磁感應(yīng)強度進行矢量疊加。因此，可先計算任一圓環(huán)產(chǎn)生的磁感應(yīng)強度，然后再進行矢量疊加。同心雙環(huán)嵌套磁體的長度為a，內(nèi)徑為d′，外徑為D′，取z軸為磁體的對稱軸，磁體中心為坐標原點建立坐標系Oxyz，如圖1所示。

根據(jù)磁荷模型[13]與高斯定律[14]，設(shè)φm為磁標量勢；ρm為磁荷分布密度；μ0為真空磁導(dǎo)率，得：

(1)

令原點磁標量勢φm(0,0,0)=0；d為點P到原點的距離；φm(x0,y0,z0)為點P的磁標量勢，則由式(1)得出點P(x0,y0,z0)的標量磁位F(x0,y0,z0)為

(a)裝配實物圖(b)示意圖

圖1 同心雙環(huán)嵌套磁體

F(x0,y0,z0)=φm(x0,y0,z0)-φm(0,0,0)=

(2)

(3)

(4)

積分區(qū)域為封閉的磁體表面S。采用圓柱坐標系，令x=rcosφ,y=rsinφ,z=±a/2，并把場點取在φ=0的平面上，即x0=r0,y0=0,z0=z0，則式(4)可用圓柱坐標系表示。對于強各向異性場的燒結(jié)釹鐵硼永磁材料的磁體，退磁場徑向分量Hmr的作用可以忽略，磁化強度M可認為始終是平行于軸線且均勻分布的。因此，用圓柱坐標系表示的式(4)可用來描述磁場分布。但是，如果需要計算空間確定點的磁感應(yīng)強度，通常需要采用復(fù)雜的有限元分析，而本文在求解過程中利用廣義二項式定理[15]，將不可積因子轉(zhuǎn)變?yōu)榭煞e因子，從而簡化了計算。

根據(jù)矢量疊加原理，同心雙環(huán)嵌套磁場所產(chǎn)生的徑向磁感應(yīng)強度Htmr(r0,z0)與軸向磁感應(yīng)強度Htmz(r0,z0)分別為內(nèi)外環(huán)形磁場產(chǎn)生的徑向磁感應(yīng)強度Himr(r0,z0)、Homr(r0,z0)與軸向磁感應(yīng)強度Himz(r0,z0)、Homz(r0,z0)的疊加，即：

(5)

2 計算結(jié)果分析

為了驗證建立的靜磁場磁感應(yīng)強度解析表達式的正確性，使用Matlab[16]并利用3維高精度數(shù)字式特斯拉計對三種磁體沿徑向和軸向的磁場分量進行了測量。通過對比理論值、測量值并結(jié)合場強分布圖，對解析表達式進行分析。

同心雙環(huán)嵌套磁體尺寸規(guī)格為：里圓環(huán)的內(nèi)徑和外徑分別為18 mm和24 mm，外圓環(huán)的內(nèi)徑和外徑分別為36 mm和58 mm，軸向厚度7.5 mm。利用同心雙環(huán)嵌套磁體的解析表達式計算距離永磁環(huán)上(或下)表面一定距離處的磁場分布，選取高度為2 mm進行計算。相應(yīng)地，在此高度上用三維特斯拉計對磁感應(yīng)強度進行多點測試。徑向和軸向磁場分量的測量結(jié)果和計算結(jié)果如圖2所示。

圖2 同心雙環(huán)嵌套磁體磁感分布

圖2中：曲線1表示磁體的磁感應(yīng)強度理論值；曲線2表示磁體的磁感應(yīng)強度測量值；曲線3表示磁體的徑向磁感應(yīng)強度測量值；曲線4表示磁體的徑向磁感應(yīng)強度理論值；曲線5表示磁體的軸向磁感應(yīng)強度理論值；曲線6表示磁體的軸向磁感應(yīng)強度測量值。由圖2可知，在原點處軸向磁感應(yīng)強度最大，然后逐漸減弱，在內(nèi)圓環(huán)外邊緣處減少到最小，即約21 mm處，而徑向磁感應(yīng)強度則由于嵌套的原因在17 mm處由0變至最大，而后逐漸減小，在24 mm處達到最小。

3 監(jiān)測試驗

將所設(shè)計的磁體放入在線式可視鐵譜儀的鐵磁磨粒收集器中,并把在線式可視鐵譜儀放置在圖3所示的試驗系統(tǒng)中。由圖3可知，潤滑油從油箱中被泵出后進入航空發(fā)動機，然后再進入在線式全流量可視鐵譜儀的鐵磁磨粒收集器。在鐵磁磨粒收集器中有觀測區(qū)，環(huán)形嵌套磁體所產(chǎn)生的磁場覆蓋該區(qū)域，當一次監(jiān)測完成時, 在線式全流量可視鐵譜儀的收集器清洗機構(gòu)開始工作，清除觀測區(qū)沉積的鐵磁磨粒，以待下一次監(jiān)測。

試驗時,采用長城4010號合成航空潤滑油，潤滑油流速1.376 m/s，被監(jiān)測航空發(fā)動機為某新型渦噴發(fā)動機，監(jiān)測時間為16 h為一時間段。采用南京航空航天大學(xué)航空安全與保障技術(shù)研究所研制的改進型DMAS II對第5時間段采集的圖像進行分析,分析結(jié)果如表1所示。統(tǒng)計結(jié)果表明，在第5時間段中共收集到22顆磨粒,絕大多數(shù)磨粒尺寸小于10～15 μm，屬于正常磨損。雖然有嚴重滑動產(chǎn)生的磨粒和切削磨粒，但是數(shù)量很少,這說明該航空發(fā)動機處于磨合期，產(chǎn)生少量的此類磨粒是正常的,而且這與該新型發(fā)動機的實際情況是相符合的。

圖3 油液在線監(jiān)測試驗系統(tǒng)

表1 磨粒統(tǒng)計結(jié)果 %

4 結(jié) 語

本研究主要有以下結(jié)論：對新型的在線式可視鐵譜儀的靜磁場磁路進行了設(shè)計。試驗證明，該靜磁場能夠在高流速下吸附鐵磁性磨粒。采用廣義二項式定理對靜磁場的磁感應(yīng)強度進行了解析計算。仿真結(jié)果表明，解析表達式具有較高的計算精度。進行了某新型渦噴發(fā)動機的試驗監(jiān)測研究，監(jiān)測結(jié)果與發(fā)動機的實際情況相符合。

[1] 左洪福，吳振鋒，陸海琴，等.DMAS智能化鐵譜分析系統(tǒng)及其應(yīng)用[J].江蘇航空，1999， 11(Z1)：60-64.

[2] 趙 方，劉 巖，謝友柏. OLF-4在線鐵譜儀的采樣數(shù)據(jù)及其特征提取的研究[J].機械科學(xué)與技術(shù)，2000，19(2)：310-312.

[3] 任國全,張培林,楊 寧,等. 基于SFS的磨粒圖像三維重構(gòu)[J]. 軍械工程學(xué)院學(xué)報，2010(2)：49-52.

[4] 張 鋼，張 堅，張海龍,等.基于等效磁荷法用蒙特卡洛法計算永磁軸承磁力[J].軸承，2013(10)：1-4.

[5] Gou Xiao-fan，Yang Yong，Zheng Xiao-jing. Analytic expression of magnetic field distribution of rectangular permanent magnets[J]. Applied Mathematics and Mechanics (English Edition)，2004，25(3)：297-306.

[6] Selvaggi J，Salo S，Kwon O，etal. Computation of the external magnetic field, near-field or far-field from a circular cylindrical magnetic source using toroidal functions[J]. IEEE Trans. Magn.，2007，43(4)：1153-1156.

[7] Harris C T，Haw D W，Handler W B，etal.Application and experimental validation of an integral method for simulation of gradient-induced eddy currents on conducting surfaces during magnetic resonance imaging[J]. Physics in Medicine and Biology，2013，58(12)：4367-4379.

[8] Ravaud R, Lemarquand G，Lemarquand V，etal. Discussion about the analytical calculation of the magnetic field created by permanent magnets[J]. Progress in Electromagnetics Research B，2009，11：281-297.

[9] Ravaud R，Lemarquand G，Lemarquand V，etal. Discussion about the analytical calculation of the magnetic field created by permanent magnets[J]. Progress in Electromagnetics Research B，2009，11：281-297.

[10] Azzerboni B，Saraceno G. Three-dimensional calculation of the magnetic field created by current-carrying massive disks[J]. IEEE Trans Magn，1998，4(5)：2601-2604.

[11] S?ren B A, S?ren E,Ilmar F S. Dynamics and stability of rigid rotors levitated by passive cylinder-magnet bearings and driven/supported axially by pointwise contact clutch[J]. Journal of Sound and Vibration，2013，332(25)：6637-6658.

[12] Furlani E P. Permanent magnet and electromechanical devices：materials，analysis and applications[M]. San Diego：Academic Press，2001.

[13] 易 中. 場論基礎(chǔ)[M]. 北京：冶金工業(yè)出版社，2013.

[14] 歐維義.場的數(shù)學(xué)描寫方法[M].長春：吉林人民出版社，1983.

[15] Goodstein R L. Fundamental concepts of mathematics [M]. New York：Pergamon Press，1979.

[16] Xiao Rujing，Yang Xiangyu，Zhu Wenjuan，etal. Modeling and simulation of dual-rotor permanent-magnet synchronous motor based on MATLAB[C]//2009 12thInternational Conference on Electrical Machines and Systems，2009.11，Tower Hall Funabori，Tokyo，Japan.Piscataway：IEEE，2009：15-18.

Design and Application of Magnetic Field of a New Online Ferrograph

HUANGWen-jie1,ZUOHong-fu2
(1．Transportation College，Huaiyin Institute of Technology，Huai’an 223003，China; 2．Civil Aviation College，Nanjing University of Aeronautics and Astronautics，Nanjing 210016，China)

With the background of a new online full flow visual ferrograph, a magnetic circuit of the static magnetic field is designed, it is axial magnetized, dual-rings nested assembly with different poles, and is filled with non-magnetic materials as the bush between the rings. Starting with the magnetic charge model, combining with the scalar magnetic potential method and using the generalized binomial theorem, the magnetic inductions of the concentric nested dual-rings are respectively analytically calculated. Comparing the theoretical value of magnetic induction with the actual measured value shows that the analytical expression well reflects the rules of magnetic induction in the static magnetic field. On this basis, the designed static magnetic field is used to monitor the performance of a new turbojet aero-engine, and then the collected debris is analyzed. The analytical result indirectly verifies that the static magnetic field can effectively adsorb ferromagnetic particles under the condition of the high volume of lubricant oil flow.

online full flow ferrograph; static magnetic field; magnetic circuit design; analytical calculation; monitoring performance

2014-04-24

國家自然科學(xué)基金項目(60939003)

黃文杰(1977-)，男，江蘇淮安人，講師，主要從事設(shè)備狀態(tài)監(jiān)測研究。Tel.:0517-83559165；E-mail:huangwj@hyit.edu.cn

V 263.6

A

1006-7167(2015)02-0054-03