李大興 夏革非 張華東 李文龍

（國網(wǎng)冀北電力有限公司承德供電公司 承德 067000）

1 引言

儲能系統(tǒng)在微電網(wǎng)中有非常大的市場前景，對電網(wǎng)的電能質(zhì)量、電網(wǎng)穩(wěn)定性以及供電可靠性都有很大的提升。飛輪儲能系統(tǒng)是一種能量轉(zhuǎn)換的儲能裝置，突破化學(xué)電池的局限，用物理的方法實現(xiàn)儲能。和其他基于化學(xué)性質(zhì)的儲能技術(shù)相比，飛輪儲能在微電網(wǎng)中具有廣闊的應(yīng)用前景。飛輪儲能系統(tǒng)中，支承高速飛輪的軸承技術(shù)是制約飛輪儲能效率和壽命的關(guān)鍵因素之一。磁軸承由于無機(jī)械摩擦、無噪聲而被廣泛應(yīng)用。無軸承電機(jī)是一種磁電機(jī)，沒有任何傳統(tǒng)的軸承或磁軸承[1-4]。它的新穎之處在于轉(zhuǎn)子懸浮力所需的電磁力是由主磁場自身產(chǎn)生的，而不是有獨立的磁軸承產(chǎn)生。為了在主磁場中同時產(chǎn)生轉(zhuǎn)矩力和洛倫茲力，需要專門設(shè)計所需要的定子和轉(zhuǎn)子。目前已經(jīng)出現(xiàn)的各種定轉(zhuǎn)子結(jié)構(gòu)中，還沒有一個轉(zhuǎn)子擁有一個大的轉(zhuǎn)矩承載力和很好的可控洛倫磁力。

通過有限元對各種不同轉(zhuǎn)子結(jié)構(gòu)（籠式轉(zhuǎn)子、磁阻轉(zhuǎn)子、外嵌永久磁鐵轉(zhuǎn)子和內(nèi)嵌永久磁鐵轉(zhuǎn)子）的無軸承電機(jī)轉(zhuǎn)矩和磁懸浮力進(jìn)行了比較性研究[5-6]，鑒于此，本文提出了一種新型的混合轉(zhuǎn)子結(jié)構(gòu)，目的是為了改進(jìn)無軸承電機(jī)的轉(zhuǎn)矩和懸浮力性能。根據(jù)得出的懸浮力模型，提出了一種新的懸浮力矢量控制策略。

2 不同轉(zhuǎn)子結(jié)構(gòu)的轉(zhuǎn)矩和懸浮力產(chǎn)生機(jī)理

為了產(chǎn)生轉(zhuǎn)矩和懸浮力，無軸承電機(jī)的定子有兩套不同極對數(shù)的定子繞組。理論上，所有的交流無刷轉(zhuǎn)子都能用于無軸承電機(jī)[7,8]。但是，為了得出哪一種轉(zhuǎn)子結(jié)構(gòu)能產(chǎn)生較大的轉(zhuǎn)矩和懸浮力，因此有必要對不同轉(zhuǎn)子結(jié)構(gòu)的無軸承電機(jī)做一個比較。如圖1 所示。

圖1 無軸承電機(jī)轉(zhuǎn)子結(jié)構(gòu)Fig.1 Rotor structures of bearingless motor

通過有限元（FEM）方法研究了懸浮力大小，相對于四個參考轉(zhuǎn)子的懸浮繞組電流。工作條件假定為：電機(jī)的氣隙是相等的，都為0.6mm，轉(zhuǎn)矩繞組的固定電流為0.45A。結(jié)果如圖2 所示，鼠籠轉(zhuǎn)子和磁阻轉(zhuǎn)子分別有最大和最小懸浮力。對于PM轉(zhuǎn)子懸浮力的承載力大致是一樣的。

圖2 不同轉(zhuǎn)子結(jié)構(gòu)懸浮力大小的比較Fig.2 Comparison of levitation forces for different rotor structures

在同樣工作條件下，不同轉(zhuǎn)子結(jié)構(gòu)計算出的電磁轉(zhuǎn)矩如圖3 所示。結(jié)果顯示無軸承電機(jī)PM 轉(zhuǎn)子產(chǎn)生的電磁轉(zhuǎn)矩比籠型和磁阻式轉(zhuǎn)子電機(jī)的大。同時也可以得出懸浮繞組電流參數(shù)對轉(zhuǎn)矩的影響是較小的，甚至可以忽略，對懸浮力控制來說是一個好的特點。

圖3 不同轉(zhuǎn)子結(jié)構(gòu)的轉(zhuǎn)矩大小比較Fig.3 Comparison of torque for different rotor structures

3 新型的混合轉(zhuǎn)子結(jié)構(gòu)

圖2 和圖3 的結(jié)果表明：對無軸承電機(jī)來說，外嵌式PM 轉(zhuǎn)子是一個較好地選擇，和其他的轉(zhuǎn)子結(jié)構(gòu)相比，其產(chǎn)生的轉(zhuǎn)矩是最大的，懸浮力也是較大的。對于在轉(zhuǎn)子和定子之間沒有偏移的電機(jī)來說上面的結(jié)論是正確。事實上，由于各種原因轉(zhuǎn)子和定子之間的偏移或多或少還是存在的。例如，如果空氣氣隙不等，在氣隙小的一邊將產(chǎn)生電磁吸力。為了研究偏移原因，用有限元計算出的轉(zhuǎn)子不同偏移量的軸向力如圖4 所示，從圖4 中可以看出，徑向力隨著轉(zhuǎn)子位移偏移量的增大而線性增加。為了使轉(zhuǎn)子回到中心位置，懸浮繞組就需要一個大的電流。但是，永磁電機(jī)通過懸浮繞組電流產(chǎn)生的氣隙磁通密度與通過永久磁鐵產(chǎn)生的相比是較小的，主要因為懸浮繞組電流產(chǎn)生的磁通還要通過一個較大的氣隙，這與空氣氣隙的數(shù)量和永久磁鐵的高度有關(guān)。因此，僅僅通過增大懸浮繞組電流而使轉(zhuǎn)子回到中心位置是很困難的。

圖4 轉(zhuǎn)子不同偏移量所受到的軸向力Fig.4 Radial force due to eccentricity

本文提出的新型無軸承電機(jī)的轉(zhuǎn)子結(jié)構(gòu)如圖5所示。兩個輔助轉(zhuǎn)子鐵心加到永磁電機(jī)的轉(zhuǎn)子上，可以增大懸浮力。混合的轉(zhuǎn)子結(jié)構(gòu)擁有永磁電機(jī)和籠型電機(jī)的優(yōu)點。

圖5 新型無軸承電機(jī)轉(zhuǎn)子結(jié)構(gòu)Fig.5 A new rotor structure of bearingless motor

為了驗證新型轉(zhuǎn)子結(jié)構(gòu)的有效性，圖6 給出了不同轉(zhuǎn)子結(jié)構(gòu)當(dāng)轉(zhuǎn)子偏離中心0.05mm 時，徑向力相對于懸浮繞組電流圖。圖6 中的四條曲線對應(yīng)于下列四種不同情況：

1.無附加鐵心的永磁電機(jī)，g0＝1.0mm；

2.有附加鐵心且氣隙長度相等的永磁電機(jī)（2lc＝15mm，gc=g0＝1.0mm）；

3.有附加鐵心且氣隙長度較小的永磁電機(jī)（2lc＝15mm，gc=0.8mm，g0＝1.0mm）；

4.有附加鐵心且氣隙小，氣隙較小的永磁電機(jī)（2lc＝10mm，gc=0.5mm，g0＝1.0mm）；

從曲線1 可以看出，對于永磁電機(jī)僅僅通過增加懸浮繞組電流不能補(bǔ)償徑向力。由曲線2 和4 可以看出用附加的轉(zhuǎn)子鐵心或減小附加轉(zhuǎn)子鐵心的空氣氣隙來補(bǔ)償徑向力是有效的。

圖6 不同轉(zhuǎn)子結(jié)構(gòu)徑向懸浮力對比Fig.6 Comparison of radial forces for different rotor structures

4 無軸承電機(jī)懸浮力控制方法

新型結(jié)構(gòu)的無軸承電機(jī)有四極三相轉(zhuǎn)矩繞組和兩極兩相的懸浮力繞組，在x-y靜止坐標(biāo)系中，三相轉(zhuǎn)矩繞組可以變換成等效的兩相繞組。如果給四極和兩極繞組提供頻率相同的、對稱的兩相交流電流，則在x-y靜止坐標(biāo)系中，相電流可以表示成：

四極和兩極繞組電流將分別產(chǎn)生兩個氣隙磁場。通過Maxwell張量法，轉(zhuǎn)子表面的徑向力可以計算出來，將該力沿轉(zhuǎn)子表面合成，x與y方向力合成結(jié)果可以表示為：

系數(shù)k的大小由下面的方程決定：

式中，N4和N2分別是四極和兩極繞組每相的等效匝數(shù)；r、1l 及 g0分別為轉(zhuǎn)子內(nèi)半徑，附加槽長度和氣隙長度。由式（1）～式（4）可以得出x和y方向力的表達(dá)式為

由式（8）、（9）可以看出，通過瞬時定子電流可以計算出x、y方向的懸浮力。但是x、y方向力的大小取決于繞組電流。為了簡化控制系統(tǒng)，最好減弱x、y方向繞組電流的相互作用。在x-y坐標(biāo)系中相電流可以變換成d-q旋轉(zhuǎn)坐標(biāo)系中新變量，x-y參考坐標(biāo)系與d-q參考坐標(biāo)系轉(zhuǎn)換電路如圖7所示。兩參考坐標(biāo)中變量的關(guān)系式如下

如果i4放在d軸上，那么角度θ由θ=atan(i4y/ i4x)來決定，如圖7所示。規(guī)定好i4的方向之后，由于i4和i2有相同的角頻率，則d-q坐標(biāo)系中的電流就變成直流成分，方程式（5）可以化簡成：

圖7 xy 坐標(biāo)與dq 坐標(biāo)的轉(zhuǎn)換關(guān)系圖Fig.7 Transformation diagram of xy coordinate and dq coordinate

利用新變量，懸浮力的表達(dá)式可以簡化為

由式（7）可以得出，在同步旋轉(zhuǎn)坐標(biāo)系中懸浮力控制轉(zhuǎn)換為對i2的控制。

5 實驗結(jié)果

提出了基于DSP TMS320C32 和Xinlinx CPLD混和轉(zhuǎn)子結(jié)構(gòu)的無軸承電機(jī)控制系統(tǒng)，上述的懸浮力矢量控制也被應(yīng)用。圖8 和圖9 顯示了在速度達(dá)到750r/min 時，懸浮力矢量控制和無懸浮力矢量控制時轉(zhuǎn)子的偏移量。每個圖形的最頂端是由渦流傳感器檢測的x 及y 方向的轉(zhuǎn)子偏移量。圖8 顯示的是無矢量控制，轉(zhuǎn)子偏移量大，由于無懸浮力控制降低了軸承的利用率。在懸浮力的矢量控制下，轉(zhuǎn)子偏移量較小，轉(zhuǎn)子的懸浮情況如圖9 所示。

圖8 無轉(zhuǎn)子矢量控制時轉(zhuǎn)子偏移軌跡Fig.8 Rotor offset trajectory without vector control

圖9 轉(zhuǎn)子矢量控制時轉(zhuǎn)子偏移軌跡Fig.9 Rotor offset trajectory with vector control

6 結(jié)論

提出的混和轉(zhuǎn)子結(jié)構(gòu)的無軸承電機(jī)不僅能產(chǎn)生大的電磁力矩和懸浮力，而且也實現(xiàn)了一種基于電流矢量的懸浮力控制方法。樣機(jī)測試結(jié)果表明，新型無軸承電機(jī)的設(shè)計和控制策略是可行的。

[1]Okada Y,Miyamoto S,Ohishi T.Levitation and torque control of internal permanent magnet type bearingless motor[J].IEEE Transactions on Control Systems Technology,1996,4(5):565-571.

[2]Suzuki T,Chiba A,Rahman M A and Fukao T.An air-gap-fluxorientation vector controller for stable operation of bearingless induction motors[J].IEEE Transactions on Industry Applications,2000,36(4):1069-1076.

[3]Fang J C,Zheng S Q,Han B C.AMB vibration control for structural resonance of double-gimbal control moment gyro with high-speed magnetically suspended rotor[J].IEEE/ASME Transactions on Mechatronics,2013,18(1):32-43.

[4]Ren Y,Fang J C.High-stability and fast-response twisting motion control for the magnetically suspended rotor system in a control moment gyro[J].IEEE/ASME Transactions on Mechatronics,2013,18(5):1625-1634.

[5]Zhang W Y,Zhu H Q.Improved model and experiment for AC-DC three-degree of freedom hybrid magnetic bearing[J].IEEE Transactions on Magnetics[J].2013,49(11):5554-5565.

[6]Zhang W Y,Zhu H Q.Precision modeling method specifically for AC magnetic bearings[J].IEEE Transactions on Magnetics,2013,49(11):5543-5553.

[7]Darbandi S M,Behzad M,Salarieh H,Mehdigholi H.Linear output feedback control of a three-pole magnetic bearing[J].IEEE/ASME Transactions on Mechatronics,2014,19(4):1323-1330.

[8]李德廣，劉淑琴.磁懸浮軸承數(shù)字控制的穩(wěn)定性分析及預(yù)補(bǔ)償算法[J].電工技術(shù)學(xué)報，2011，26(6):108-112.Li Deguang,Liu Shuqin.Stability analysis of the AMB digital control system and its predictioncompensatory algorithm[J].Transactions of China Electrotechnical Society,2011，26(6):108-112.