交流混合型磁軸承磁場及懸浮力特性分析

張 松，張維煜，朱熀秋

(1.南京師范大學(xué)，江蘇南京210097;2.江蘇大學(xué)，江蘇鎮(zhèn)江212013)

0 引 言

磁軸承作為一種極有前途的軸承，依據(jù)驅(qū)動電流類型不同，大致可分為兩類:①由直流電驅(qū)動的磁軸承，即稱為直流磁軸承［1］;②由交流電驅(qū)動的磁軸承，即稱為交流磁軸承［1］。功率放大器對于直流磁軸承來說，它具有成本高、體積大的缺點。驅(qū)動一個徑向磁軸承或者需要用采用4路單極性功率放大器，或采用2路雙極性功率放大器亦可;而對于交流磁軸承，由于采用三相功率逆變器驅(qū)動，因此具有應(yīng)用價格便宜、技術(shù)成熟且體積小的優(yōu)點，而且驅(qū)動一個徑向磁軸承只需要一個三相功率逆變器就可實現(xiàn)控制目的［1－4］。依據(jù)磁軸承中產(chǎn)生的懸浮力的原理不同，磁軸承大致可分為三種:偏置磁通來源于永磁體，控制磁通來源于電磁鐵的磁軸承，即為混合型磁軸承;偏置磁通與控制磁通完全來源于永磁體的磁軸承，即為被動型磁軸承;偏置磁通與控制磁通完全來源于電磁鐵，即為主動型磁軸承［3］。其中，混合型磁軸承具有明顯的優(yōu)勢，由于永磁體的存在，使得其靜態(tài)偏置磁通是永磁體負責產(chǎn)生，解決了主動磁軸承中的偏置磁通來源于電磁鐵的弊端，在很大程度上使得功率放大器的功耗有所減少［4］。又由于在混合磁軸承中，偏置磁通無需電磁鐵提供，因此使得電磁鐵的安匝數(shù)大大減少，且磁軸承的體積也有所減小。由此，交流混合型磁軸承在高速電主軸、高速飛輪儲能、無軸承電機等場合具有潛在應(yīng)用前景［1－4］，已成為磁軸承研究與發(fā)展的一個熱點方向。

本文不僅對交流混合型磁軸承的組成結(jié)構(gòu)進行了介紹，而且對該交流混合型磁軸承的工作基本原理也進行了闡述，還創(chuàng)建了該磁軸承的實體模型，利用Ansoft軟件仿真并分析了該交流磁軸承的磁通分布狀況(氣隙及磁極內(nèi)部)。并分別以磁軸承的有限元實體模型和數(shù)學(xué)模型為研究對象，利用Ansoft有限元分析軟件與MATLAB軟件對磁軸承懸浮力的非線性變化規(guī)律進行了計算分析。仿真結(jié)果表明:針對于該磁軸承內(nèi)的磁場分布情況來說，其分析結(jié)果與從理論上分析其磁場走向及分布情況吻合，該磁軸承的參數(shù)設(shè)計的是有效的;基于數(shù)學(xué)模型和有限元實體模型的交流混合型磁軸承特性分析結(jié)果一致，該磁軸承所建立的數(shù)學(xué)模型與磁軸承實際工作情況下的模型幾乎一致。

1 結(jié)構(gòu)與工作原理

圖1為交流混合型磁軸承的整體結(jié)構(gòu)視圖。從圖中可以看出，交流混合型磁軸承的定子是由一組對稱的定子片所構(gòu)成，每個定子片上均分布有3個沿圓周互呈120°的凸極，因此該磁軸承一共有6個磁極。兩個定子片對應(yīng)的凸極分布角度相同，這些凸極上分別繞制有三相控制線圈，且兩個定子片上同一角度對應(yīng)的凸極上纏繞的線圈以串聯(lián)形式連接。永磁體與隔磁鋁環(huán)位于對稱的定子片之間。

圖1 關(guān)于磁軸承整體的結(jié)構(gòu)視圖

工作原理如下:對于交流磁軸承來說，其原理是與無軸承電機的工作原理有緊密關(guān)聯(lián)的，也就是說交流磁軸承相當于其轉(zhuǎn)矩繞組極對數(shù)為0，懸浮力繞組極對數(shù)為1，且采用三相功率逆變器給懸浮力繞組提供控制電流的一種特殊的無軸承電機［5－7］。在徑向方向上如果存在擾動，則轉(zhuǎn)子因此發(fā)生偏移。此時傳感器負責采集徑向方向上的偏移x與y，此位移信號經(jīng)過控制器處理，輸出的控制信號再經(jīng)過三相功率逆變器進行驅(qū)動，輸出的三相控制信號使得作用在轉(zhuǎn)子上的由三相控制線圈和永磁體共同作用的合力方向與轉(zhuǎn)子發(fā)生偏移方向相反，其大小可以使轉(zhuǎn)子穩(wěn)定懸浮在所設(shè)定的平衡位置。

圖2為交流混合型磁軸承內(nèi)部的磁路分布情況。圖中，實線代表靜態(tài)偏置磁通的走向，是由永磁體造成的實線的起點出自于N極(永磁體)，分別流經(jīng)一側(cè)定子后進入相應(yīng)氣隙、然后穿過轉(zhuǎn)子后進入另一。側(cè)氣隙、然后進入另一側(cè)定子，之后回歸于S極;虛線代表控制磁通的走向，是由電磁造成的。

圖2 磁軸承的內(nèi)部磁路分布

2 數(shù)學(xué)模型

圖3為永磁體作用繪制的磁路圖。其中，F(xiàn)m為磁動勢的值，Φm為總磁通量，GA1、GB1、GC1、GA2、GB2、GC2分別為氣隙磁導(dǎo)值，ΦA(chǔ)1、ΦB1、ΦC1、ΦA(chǔ)2、ΦB2、ΦC2分別為氣隙磁通值。

圖3 永磁體作用情況下的磁路圖

圖4為電磁鐵作用下的磁路圖。如果轉(zhuǎn)子中心沿X、Y軸發(fā)生偏移x、y，則各氣隙處磁導(dǎo):

圖4 電磁鐵作用下的磁路圖

式中:μ0為真空磁導(dǎo)率;S為磁極面積;δ為徑向氣隙長度。

根據(jù)磁路基爾霍夫定律ΣF=0和ΣΦ=0，施加控制磁通后，各氣隙處由控制磁通和永磁偏置磁通相疊加構(gòu)成合成磁通，其中各氣隙處的偏置磁通及控制磁通:

將式(3)進行線性化處理，首先將其寫成泰勒展開的形式，由于轉(zhuǎn)子發(fā)生在平衡位置附近處的偏移量比起氣隙長度來說，微乎甚微，因此文中將二階以上無窮小量忽略，得:

根據(jù)三相電流的關(guān)系iA+iB+iC=0，可得到:

3 基于實體模型的磁場分析

本文研究的交流混合型磁軸承部分結(jié)構(gòu)及電磁參數(shù):最大承載力Fmax≥350 N，氣隙長度 δ=0.5 mm，飽和磁通密度Bsat=1.75 T，取BH曲線線性區(qū)最大磁通密度Bmax=0.9Bsat=1.575 T，軸向長度 L=38 mm，定子最大外徑D=126 mm，磁極面積S=363 mm2，控制線圈采用的漆包線，其標稱直徑為0.67 mm，線圈匝數(shù)N=366匝，永磁體采用釹鐵硼。

為防止鐵心發(fā)生磁飽和現(xiàn)象且可以充分利用永磁體的最大磁性能，直流磁軸承一般將工作點B0選定在定子鐵心磁化曲線的線性范圍內(nèi)最大磁通密度Bmax的1/2處。區(qū)別于直流磁軸承，交流混合型磁軸承的工作點B0則選在定子鐵心磁化曲線的Bmax的1/3處，為0.525 T。由于本文研究的交流混合型磁軸承有極靴的存在，雖然流過磁路中的磁通相同，但極靴面積大于磁路中其它橫截面積，所以氣隙處磁通密度小于定子鐵心的磁通密度。由于在理論計算氣隙磁通密度時，一般圍繞定子鐵心的磁化曲線展開進行，所以氣隙磁通密度值在磁軸承工作過程中應(yīng)該小于理論計算值。本文中提到的磁極處磁通密度是基于磁極中選定的一個特定橫截面上展開計算分析的。

圖5是永磁偏置情況下的氣隙磁通密度值與磁極內(nèi)磁通密度值。由圖5可知，氣隙中分布的永磁磁場是關(guān)于中心對稱的，并且磁軸承中分布的永磁磁場也是對稱的。在氣隙中的偏置磁通密度值為0.49 T，其值小于偏置磁通密度的計算值0.525 T，與理論分析相吻合，即對于氣隙磁通密度值來說，其實際工作時的磁通密度值比理論計算值小。

圖5 永磁偏置時磁極中與氣隙中對應(yīng)的磁通密度極坐標圖

圖6是永磁偏置與施加最大控制電流情況下的氣隙磁通密度值。由圖6可知，永磁偏置情況下，氣隙處磁通密度值為0.49 T，施加最大控制電流時，A相磁極下對應(yīng)氣隙處的磁通密度值達到最大，B、C相磁極下的磁通密度值幾乎抵消為零。A相磁極下氣隙處磁通密度值約為1.1 T，其值小于理論計算的最大磁通密度值Bmax=1.575 T，與理論分析一致。

圖6 兩種情況下氣隙中對應(yīng)的磁通密度極坐標圖

4 非線性分析

由式(3)可知，懸浮力Fi與磁通量(氣隙)的二次方成正比，而控制電流和位移又決定氣隙處磁通量的大小，由此Fi事實上只與x、y和ix、iy有關(guān)。圖7(a)、圖8(a)為不同位移處的力/電流關(guān)系曲線簇(基于磁軸承有限元實體模型，用Ansoft軟件對磁軸承懸浮力的非線性特性進行分析)。另外，基于磁軸承數(shù)學(xué)模型，用MATLAB軟件做出同樣條件下不同位移處的力/電流關(guān)系曲線簇，如圖7(b)、圖8(b)所示。

圖7(a)、圖7(b)按著箭頭所示的方向(轉(zhuǎn)子在x軸方向位移從0逐漸變化到0.1 mm、0.2 mm，0.3 mm、0.4 mm)時，電流 ix與轉(zhuǎn)子受力 Fx、Fy的關(guān)系曲線圖。由圖可知，受徑向二自由度之間的交叉耦合關(guān)系影響所致，徑向懸浮力Fx與控制電流ix的關(guān)系呈現(xiàn)了不對稱性。安匝數(shù)從－150安匝變化到150安匝時，曲線的線性特性很好，控制該范圍的磁軸承能達到較好的效果。可以看出隨著轉(zhuǎn)子在徑向x的方向偏移量增大，其力/電流關(guān)系曲線的線性度變差。觀察發(fā)現(xiàn)，轉(zhuǎn)子在徑向x方向0.4 mm處發(fā)生偏移時，平衡位置的線性度已經(jīng)不是很理想，但是偏移發(fā)生在0.3 mm以內(nèi)時，平衡位置的線性度仍然比較理想。說明磁軸承工作在氣隙很小的范圍內(nèi)時，為取得優(yōu)良的控制效果，磁軸承如果正常工作的氣隙比較小時，無需采用復(fù)雜的非線性控制器就可以達到目的。如果磁軸承工作在氣隙范圍比較大的情況下，那么在設(shè)計控制器時需考慮補償，由于大氣隙造成的非線性和由于三級結(jié)構(gòu)造成的非對稱性，例如可采用神經(jīng)網(wǎng)絡(luò)逆控制器等進行解耦控制或者采用參數(shù)自整定模糊PID控制算法進行分段控制。

圖8(a)、8(b)按著箭頭所示的方向(轉(zhuǎn)子在y軸方向的位移從0 逐漸變化到0.1 mm、0.2 mm、0.3 mm)，電流iy與轉(zhuǎn)子受力Fy、Fx的關(guān)系曲線圖。由圖可知，電流iy變化時，不僅產(chǎn)生徑向力Fy，也產(chǎn)生在x軸負方向上的力Fx，控制電流iy的安匝數(shù)在－150安匝到150安匝之內(nèi)變化時，F(xiàn)x很小，F(xiàn)x的上升幅度與安匝數(shù)的增加幅度成正比。

綜上，圖7顯示懸浮力Fx與控制電流ix相對于平衡位置不是對稱的，這是與該磁軸承的結(jié)構(gòu)(三極)的非對稱性所密切相關(guān)。但是圖8中顯示懸浮力Fx與控制電流ix相對于平衡位置是對稱的。因此可得出，控制該磁軸承時，在x、y方向的控制難易程度是有所區(qū)別。在x方向調(diào)節(jié)磁軸承穩(wěn)定懸浮時，其控制參數(shù)會影響y方向?qū)Υ泡S承的控制效果，而控制y方向的轉(zhuǎn)子平衡時，改變其控制參數(shù)則不會對x方向產(chǎn)生任何影響。

圖7(b)、圖8(b)相比于圖7(a)、圖8(a)是在完全理想的情況下的力/電流關(guān)系曲線，但是仍然可以反映表達此軸承的力/電流關(guān)系的趨勢，與基于磁軸承有限元實體模型進行仿真的結(jié)果作比較，相似的結(jié)果證明了該磁軸承數(shù)學(xué)模型建立準確，逼近該交流型混合磁軸承的正常運行情況。

5 結(jié) 語

針對交流混合型磁軸承，本文應(yīng)用經(jīng)典有效的建模方法對該磁軸承進行了數(shù)學(xué)模型的建立;對采用有限元Ansoft磁場分析軟件建立的磁軸承實體模型進行仿真分析，可以很清楚地觀察到該磁軸承內(nèi)的磁通分布情況;針對磁軸承的兩種模型(有限元實體模型與數(shù)學(xué)模型)，對其懸浮力非線性變化規(guī)律特性進行了計算分析。由仿真結(jié)果可知:對該磁軸承內(nèi)的磁場分布情況來說，其分析結(jié)果與從理論上分析其磁場走向及分布情況吻合、兩種模型分析下的懸浮力非線性特性分析結(jié)果一致，表明對該磁軸承的結(jié)構(gòu)參數(shù)計算結(jié)果來說，其設(shè)計是合理有效的，對數(shù)學(xué)模型的計算結(jié)果來說，其建立是正確的，且結(jié)果接近磁軸承實際運行情況下的結(jié)果。

［1］ 朱熀秋，沈玉祥，鄔清海，等.交流混合型磁軸承建模與控制系統(tǒng)［J］.中國電機工程學(xué)報，2009，29(18):100 －105.

［2］ 朱熀秋，張仲，諸德宏，等.交直流三自由度混合磁軸承結(jié)構(gòu)與有限元分析［J］.中國電機工程學(xué)報，2007，27(12):77－81.

［3］ Zhu H，Xie Z，Zhu D.Principles and parameter design for AC － DC three－ degree freedom hybrid magnetic bearing［J］.Chinese Journal of Mechanical Engineering，2006，19(4):534 －539.

［4］ Zhu H，Chen H，Xie Z，et al.Configuration and control for AC －DC three degrees of freedom hybrid magnetic bearings［C］//10th International Symposium on Magnetic Bearings.Martigny，Switzerland，2006.

［5］ 王曉琳，鄧智泉.無軸承異步電機磁場定向控制策略分析［J］.中國電機工程學(xué)報，2007，27(27):77 －82.

［6］ Zhou Y，Zhu H，Li T，et al.Decoupling control of magnetically levitated induction motor with inverse system theory［C］//CES/IEEE 5th International Power Electronics and Motion Control Conference.Shanghai，China，2006:1 －5.

［7］ Tera T，Yamauchi Y，Chiba A，et al.Performances of bearingless and sensorless induction motor drive based on mutual inductances and rotor displacements estimation［J］.IEEE Transactions on Industrial Electronics，2006，53(1):187－194.