倪有源，葛木明，黃 亞，何 強

(1.合肥工業(yè)大學(xué)，合肥 230009;2.工業(yè)節(jié)電與電能質(zhì)量控制省級協(xié)同創(chuàng)新中心，合肥 230601)

無刷直流永磁屏蔽電機參數(shù)計算

倪有源1, 2，葛木明1，黃 亞1，何 強1

(1.合肥工業(yè)大學(xué)，合肥 230009;2.工業(yè)節(jié)電與電能質(zhì)量控制省級協(xié)同創(chuàng)新中心，合肥 230601)

屏蔽電機主要用于各類液體循環(huán)系統(tǒng)中。由于無刷直流永磁屏蔽電機的端部漏感較大，二維方法無法精確計算電機性能，需要用三維數(shù)值方法計算。采用三維瞬態(tài)有限元方法分析了一臺無刷直流永磁屏蔽電機的電磁場，并運用能量攝動法計算了不同轉(zhuǎn)子位置角時電機定子繞組的自感和互感值，且對其進(jìn)行了傅里葉分析。在此基礎(chǔ)上，采用場路耦合的方法，計算了該電機在額定工況下的各種損耗、電磁轉(zhuǎn)矩和效率等性能參數(shù)。

無刷直流；永磁屏蔽電機；三維瞬態(tài)有限元法；場路耦合法

0 引 言

管道屏蔽電泵廣泛用于各類液體循環(huán)系統(tǒng)中。管道屏蔽電泵中使用的電機是屏蔽電機[1]，為了隔離電機中的電氣部分和傳輸液體，在定子內(nèi)徑和轉(zhuǎn)子外徑處分別安裝了一層不銹鋼屏蔽套。

傳統(tǒng)的屏蔽電機采用感應(yīng)式結(jié)構(gòu)，但隨著節(jié)能標(biāo)準(zhǔn)的不斷提高，感應(yīng)屏蔽電機將逐漸被淘汰。永磁屏蔽電機可以減小屏蔽電機轉(zhuǎn)子的鐵耗，提高屏蔽泵系統(tǒng)的效率。因此對永磁屏蔽電機的研究具有重要的理論價值和工程價值。國內(nèi)外學(xué)者對永磁屏蔽電機展開了廣泛的研究。江蘇大學(xué)分析研究了永磁屏蔽電機的電磁設(shè)計和效率[2-3]。哈爾濱理工大學(xué)提出了雙屏蔽復(fù)合轉(zhuǎn)子永磁屏蔽電機，并分析研究其損耗和性能[4]。同時，一些學(xué)者研究了永磁屏蔽電機各損耗[5-6]，分析了永磁屏蔽電機的軸向電磁力[7]，研究了自起動永磁屏蔽電機[8]，計算了永磁同步屏蔽電機的三維磁場[9]等。總之，永磁屏蔽電機由于效率較高，是未來的發(fā)展方向。

無刷直流永磁屏蔽電機的精確計算直接影響到電機性能。由于無刷直流永磁屏蔽電機的端部漏感較大，二維方法無法精確計算。本文采用三維瞬態(tài)有限元法計算了一臺無刷直流永磁屏蔽電機的三維電磁場和繞組電感，并且采用場路耦合法計算了該電機額定狀態(tài)下的各項性能參數(shù)。

1 無刷直流永磁屏蔽電機數(shù)學(xué)模型

無刷直流永磁屏蔽電機采用兩相導(dǎo)通星形三相六狀態(tài)的工作方式。在這種工作方式下共有六個工作狀態(tài)，每個工作狀態(tài)下電機兩相導(dǎo)通[10]。

無刷直流永磁屏蔽電機的數(shù)學(xué)模型[11]主要由電壓方程、反電勢方程、電磁轉(zhuǎn)矩方程和機械運動方程等組成。

1.1 電壓方程

無刷直流永磁屏蔽電機的定子電壓方程：

(1)

式中：Ua，Ub和Uc分別為定子繞組三相電壓；ea，eb和ec分別為定子繞組的三相反電動勢；ia，ib和ic分別為定子繞組三相電流；R為定子繞組每相電阻；L和M分別為定子繞組的自感和互感。

1.2 反電動勢方程

無刷直流永磁屏蔽電機每相感應(yīng)電動勢：

(2)

式中：E為每相感應(yīng)電動勢；Bδ為氣隙磁密；l為導(dǎo)體軸向有效長度；D為電樞內(nèi)徑；W為每相繞組串聯(lián)匝數(shù)；n為電機轉(zhuǎn)速。

1.3 電磁轉(zhuǎn)矩方程

無刷直流永磁屏蔽電機的電磁轉(zhuǎn)矩方程：

(3)

式中：ωs為轉(zhuǎn)子機械角速度。

1.4 機械運動方程

無刷直流永磁屏蔽電機的機械運動方程：

(4)

式中：J為電機的轉(zhuǎn)動慣量；TL為電機的負(fù)載轉(zhuǎn)矩；B為摩擦系數(shù)。

2 無刷直流永磁屏蔽電機三維磁場分析

2.1 電機結(jié)構(gòu)參數(shù)

分析的無刷直流永磁屏蔽電機用于屏蔽電泵系統(tǒng)中。電機的額定電壓為24V，額定轉(zhuǎn)速為3 000r/min，額定轉(zhuǎn)矩為110mN·m。電機的主要結(jié)構(gòu)參數(shù)如表1所示，其中，定子內(nèi)徑和轉(zhuǎn)子外徑都是不包括屏蔽套的尺寸。

表1 無刷直流永磁屏蔽電機的結(jié)構(gòu)參數(shù)

無刷直流永磁屏蔽電機的繞組直接影響電機性能，因此優(yōu)化設(shè)計出一組適合的繞組參數(shù)尤為重要，表2為經(jīng)過優(yōu)化設(shè)計后的永磁屏蔽電機的繞組參數(shù)。

表2 無刷直流永磁屏蔽電機的繞組參數(shù)

2.2 電機三維磁場分析

根據(jù)無刷直流永磁屏蔽電機的結(jié)構(gòu)及材料參數(shù)，利用有限元軟件構(gòu)建電機的三維有限元模型，建立的電機三維模型如圖1所示。該電機采用6槽4極結(jié)構(gòu)，定子采用牌號為50W470的硅鋼片。轉(zhuǎn)子無鐵心，完全由永磁體構(gòu)成，永磁材料為鐵氧體，極弧系數(shù)為1。不僅結(jié)構(gòu)簡單，便于生產(chǎn)制造，而且成本低。

(a)定子(b)轉(zhuǎn)子

圖1 無刷直流永磁屏蔽電機的三維模型

電機端部的漏磁場是三維分布的，并決定了端部漏感。由于該電機的定子結(jié)構(gòu)導(dǎo)致端部漏電感較大，用二維有限元法無法精確計算，因此需要采用三維有限元方法計算。運用三維瞬態(tài)有限元法對電磁場進(jìn)行分析，得到任意時刻電機的三維磁場分布。某時刻定轉(zhuǎn)子的磁密分布如圖2所示。

(a)定子磁密(b)轉(zhuǎn)子磁密

圖2 無刷直流永磁屏蔽電機內(nèi)的磁密分布

屏蔽套在電機運行中會產(chǎn)生渦流，形成渦流損耗，對電機性能造成一定的影響，因此需要精確計算。定子屏蔽套的渦流損耗分布如圖3所示。通過后處理功能，即可計算獲得定子屏蔽套的三維渦流損耗。

圖3 定子屏蔽套的三維渦流損耗分布

2.3 電樞繞組的電感計算

無刷直流永磁屏蔽電機三相電樞繞組分別相差120°電角度。定子與轉(zhuǎn)子的相對位置影響電機的電磁場分布，從而影響電樞繞組的自感和互感。于是可以通過計算電機內(nèi)部不同位置的三維電磁場。如果電流有十分微小的攝動，磁場能量也相應(yīng)改變，于是可計算繞組電感，這種方法稱為能量攝動法。運用該方法，可獲得三相電樞繞組的電感參數(shù)。經(jīng)過計算，A相繞組的自感與定轉(zhuǎn)子的相對位置角θ之間的關(guān)系曲線如圖4所示。

圖4 A相繞組的自感波形

自感Laa的傅里葉級數(shù)形式：

(5)

式中：n為諧波次數(shù)；A0為自感直流分量；φn為n次諧波的相位角。

對Laa進(jìn)行FFT諧波分析，得到的各次諧波幅值如表3所示。

表3 自感Laa各次諧波成分

以定轉(zhuǎn)子的相對位置角θ為橫坐標(biāo)，A相繞組與B相繞組間的互感波形如圖5所示。

圖5 A相繞組與B相繞組的互感波形

同樣對Mab進(jìn)行FFT分析，得到的各次諧波幅值如表4所示。

表4 互感Mab各次諧波成分

3 無刷直流永磁屏蔽電機性能參數(shù)計算

本文采用場路耦合的方法對電機性能參數(shù)進(jìn)行計算。圖6為使用Maxwell軟件外電路搭建的無刷直流永磁屏蔽電機的電路部分，通過控制該部分電路，可以實現(xiàn)兩相導(dǎo)通星形三相六狀態(tài)。

建立無刷直流永磁屏蔽電機有限元模型和控制電路模型后，就可以利用場路耦合法計算該電機的參數(shù)性能。

無刷直流永磁屏蔽電機的損耗主要包括定子鐵心損耗、定子繞組銅耗、屏蔽套渦流損耗以及轉(zhuǎn)子旋轉(zhuǎn)的水摩擦損耗等。

圖6 無刷直流永磁屏蔽電機的控制電路圖

額定轉(zhuǎn)速下，無刷直流永磁屏蔽電機定子鐵心損耗隨時間的變化曲線如圖7所示。于是得到穩(wěn)定運行時的定子鐵心損耗為0.568W。

圖7 無刷直流永磁屏蔽電機的鐵耗

額定轉(zhuǎn)速下，無刷直流永磁屏蔽電機定子繞組銅耗隨時間的變化曲線如圖8所示。于是得到穩(wěn)定運行時定子繞組銅耗為9.396W。

圖8 無刷直流永磁屏蔽電機的銅耗

通過三維有限元法計算屏蔽套渦流損耗，得到屏蔽套渦流損耗值為0.804W。另外計算水摩擦損耗為0.479W。將各項損耗相加，得到無刷直流永磁屏蔽電機的總損耗∑p=11.25W。

額定轉(zhuǎn)速下，無刷直流永磁屏蔽電機三相電流隨時間的變化曲線如圖9所示?？梢缘玫椒€(wěn)定時三相電流的有效值為4.04A。

圖9 無刷直流永磁屏蔽電機的三相電流

額定轉(zhuǎn)速下，無刷直流永磁屏蔽電機感應(yīng)電動勢隨時間的變化曲線如圖10所示，進(jìn)行FFT諧波分析，得到感應(yīng)電動勢的基波幅值為8.06V。

圖11為額定轉(zhuǎn)速下無刷直流永磁屏蔽電機的電磁轉(zhuǎn)矩波形。對數(shù)據(jù)進(jìn)行分析，可以得到穩(wěn)定時電磁轉(zhuǎn)矩的平均值為117.3mN·m，峰峰值為58.71mN·m，顯然轉(zhuǎn)矩脈動很大。除了齒槽轉(zhuǎn)矩以及電樞反應(yīng)引起的轉(zhuǎn)矩脈動，很大一部分是由電流換相引起的。

圖10 無刷直流永磁屏蔽電機的三相感應(yīng)電動勢

圖11 無刷直流永磁屏蔽電機的電磁轉(zhuǎn)矩

由于無刷直流永磁屏蔽電機的電磁轉(zhuǎn)矩為117.3 mN·m，可計算出輸出功率為36.48 W。

因此，無刷直流永磁屏蔽電機的效率：

(6)

4 結(jié) 語

無刷直流永磁屏蔽電機的性能參數(shù)需要精確計算。本文采用三維瞬態(tài)有限元方法分析了一臺24 V無刷直流永磁屏蔽電機的三維磁場分布，并用能量攝動法計算獲得電樞繞組的自感和互感值。然后利用場路耦合法，計算獲得屏蔽套的渦流損耗以及定子的鐵耗和銅耗，從而得到電機的電磁轉(zhuǎn)矩和效率等參數(shù)。本文對于無刷直流永磁屏蔽電機的設(shè)計及性能分析提供了一定的理論參考。

[1] 國珍，丁新平，田艷兵．復(fù)合轉(zhuǎn)子屏蔽電動機運行特性及電磁參數(shù)研究[J]．微特電機，2012，40(3)：22-26．

[2] 高翠蘭．永磁屏蔽泵的設(shè)計方法研究[D]．鎮(zhèn)江：江蘇大學(xué)，2009．

[3] 王利偉．高效屏蔽泵的性能研究[D]．鎮(zhèn)江：江蘇大學(xué)，2010．

[4] 梁艷萍，張建濤，索文旭．雙屏蔽復(fù)合轉(zhuǎn)子電機渦流損耗分析[J]．中國電機工程學(xué)報，2009，29(24)：78-83．

[5] BURKHARDT Y，HUTH G，URSCHEL SY．Eddy current losses in PM canned motors[C]//The XIX International Conference on Electrical Machines，Rome，Italy.2010：1-7．

[6] 張曉晨，李偉力，曹君慈．屏蔽電機屏蔽損耗與電機性能的計算與分析[J]．哈爾濱工業(yè)大學(xué)學(xué)報，2007，39(9)：1422-1426．

[7] AN Y J，LIU G M，MENG Z J，et al．Analysis and experiment of magnetic force in permanent magnet axial thrust balance structure for canned motor pump[C]//The 5th International Conference on Intelligent Control and Information Processing，Dalian，China.2010：732-734．

[8] PERALTA-SANCHEZ E，SMITH AC，RODRIGUEZ-RIVAS J J．Steady-state analysis of a canned line-start PM motor[J]．IEEE Transactions on Magnetics，2011，47(10)：4080-4083．

[9] 倪有源，黃亞，趙亮．新型結(jié)構(gòu)永磁屏蔽電機三維磁場分析和電感計算[J]．電工技術(shù)學(xué)報，2015，30(1)：98-104．

[10] 汪俊杰．新型無刷直流電機驅(qū)動系統(tǒng)研究[D]．南京：南京航空航天大學(xué)，2008．

[11] 黃亞. 新型永磁屏蔽電機本體與控制系統(tǒng)研究[D]. 合肥：合肥工業(yè)大學(xué), 2014.

Parameters Computation for Brushless DC PM Canned Motor

NIYou-yuan1, 2，GEMu-ming1，HUANGYa1，HEQiang1

(1.Hefei University of Technology，Hefei 230009，China;2.Collaborative Innovation Center of Industrial Energy-Saving and Power Quality Control，Hefei 230601，China)

Canned motors are widely used in all kinds of liquid circulation systems. Due to the large end leakage inductance of the brushless DC permanent magnet (PM) canned motor, the traditional two dimensional (2D) method can't be used to compute its performance accurately, so a three dimensional (3D) numerical method is required. The electromagnetic fields of the motor were analyzed using a 3D transient finite element method (FEM). Besides, using a energy perturbation method, the self-inductance and mutual-inductance values of the stator windings were computed in terms of different rotor positions. On this basis, the performance parameters of the motor under the rated condition were calculated using a field-circuit coupled method.

brushless DC； PM canned motor； 3D transient FEM； field-circuit coupled method

2015-11-11

安徽省自然科學(xué)基金項目(1508085ME89)

TM33

A

1004-7018(2016)05-0006-04

倪有源(1976-), 男,博士,副教授,研究方向為電機設(shè)計及其控制。