雙定子杯形轉(zhuǎn)子感應(yīng)電機(jī)磁路計算

姚江帆，歐陽斌，翟小飛

（海軍工程大學(xué)艦船綜合電力技術(shù)國防科技重點實驗室，武漢 430033）

雙定子杯形轉(zhuǎn)子感應(yīng)電機(jī)磁路計算

姚江帆，歐陽斌，翟小飛

（海軍工程大學(xué)艦船綜合電力技術(shù)國防科技重點實驗室，武漢 430033）

為分析雙定子杯形轉(zhuǎn)子電機(jī)性能，分別采用傳統(tǒng)磁路計算法和分布磁路法對其進(jìn)行磁路計算，研究內(nèi)定子和外定子各部分的磁通密度和磁壓降。迭代計算給定相電壓或勵磁電流下的空載特性和激磁電感并與實驗進(jìn)行對比，結(jié)果與實測值相吻合，證明了分布磁路法的準(zhǔn)確性和普適性。

雙定子 感應(yīng)電機(jī) 磁路計算 分布磁路法 空載特性

0 引言

雙定子杯形轉(zhuǎn)子感應(yīng)電機(jī)是一種特殊結(jié)構(gòu)的多相電機(jī)，其結(jié)構(gòu)如圖1所示，其內(nèi)外兩套獨立定子繞組以及空心杯形薄壁轉(zhuǎn)子的結(jié)構(gòu)使其具有轉(zhuǎn)動慣量低、響應(yīng)速度快、轉(zhuǎn)矩密度大、轉(zhuǎn)矩波動小等優(yōu)點[1]。小型雙定子杯形轉(zhuǎn)子電機(jī)在多個領(lǐng)域都有應(yīng)用，文獻(xiàn)[2]介紹了一種應(yīng)用于機(jī)器人直驅(qū)的高扭矩電機(jī)，文獻(xiàn)[3]詳細(xì)設(shè)計了一種應(yīng)用于能量回收領(lǐng)域的小型雙定子杯形轉(zhuǎn)子電機(jī)，文獻(xiàn)[4]介紹了杯形測速發(fā)電機(jī)的計算。這些電機(jī)雖然在外觀上類似本文所研究電機(jī)的縮小版，但它們的內(nèi)定子均未開槽，沒有嵌放繞組，所以從本質(zhì)上看仍是三相感應(yīng)電機(jī)，與本文電機(jī)在具體結(jié)構(gòu)、功率/轉(zhuǎn)矩密度上都有較大差異，因而表現(xiàn)出不同的性能。

為研究該雙定子杯形轉(zhuǎn)子電機(jī)的性能，首先應(yīng)進(jìn)行電機(jī)的磁路計算?，F(xiàn)有的主磁路計算方法主要有：解析磁路法，磁網(wǎng)格法和有限元法[5]。傳統(tǒng)的解析磁路法把電機(jī)分成每段包含一對極的若干扇形段，其中每一個扇形段的磁場簡化為閉合磁回路，根據(jù)全電流定律，通過計算磁回路的磁場強(qiáng)度得到勵磁磁勢，再分別計算氣隙磁密、定子齒部磁密和軛部磁密，經(jīng)過飽和系數(shù)迭代最后得到電機(jī)的空載曲線。傳統(tǒng)磁路法依賴極弧系數(shù)、氣隙磁場波形系數(shù)和軛部磁壓降矯正系數(shù)等需要查表的近似處理方式，這種方法針對正弦供電的電機(jī)雖然能快速求出最大磁密，但對于越來越普遍的多相電機(jī)適用性不強(qiáng)。近年來，有限元法發(fā)展迅速，得到廣泛應(yīng)用，但由于其計算周期長，靈活度較低，不適合電機(jī)前期設(shè)計和后期優(yōu)化。分布磁路法[6]可以有效解決上述問題，這種方法一方面對磁路進(jìn)行了網(wǎng)格剖分，不使用過多修正系數(shù)，能較細(xì)致地刻畫出電機(jī)氣隙磁場波形，在適用性和靈活性上優(yōu)于傳統(tǒng)磁路法；另一方面它通過MATLAB編程即可實現(xiàn)，計算快速方便，優(yōu)于有限元法。

圖1 雙定子杯形轉(zhuǎn)子電機(jī)結(jié)構(gòu)示意圖

本文分別用傳統(tǒng)的解析磁路法和分布磁路法對雙定子電機(jī)進(jìn)行磁路計算，求出電機(jī)空載特性曲線和激磁電感參數(shù)并與實驗結(jié)果作對比，驗證了新方法在雙定子杯形轉(zhuǎn)子電機(jī)磁路計算中的準(zhǔn)確性和有效性。

2 傳統(tǒng)磁路法簡要回顧

根據(jù)電機(jī)設(shè)計中介紹的傳統(tǒng)磁路法[7]，構(gòu)建雙定子杯形轉(zhuǎn)子電機(jī)磁路模型，以外定子為例，設(shè) 1～2段為氣隙段，2～3段為外定子齒段，3～4段為外定子軛段，內(nèi)定子方法一致，如圖2所示。不同于傳統(tǒng)三相異步電機(jī)，本文研究的雙定子電機(jī)轉(zhuǎn)子是薄壁光滑、沒有繞組的，因此轉(zhuǎn)子不是空載磁路的一部分。內(nèi)定子和外定子繞組共同建立的旋轉(zhuǎn)磁場作用于轉(zhuǎn)子，使其感應(yīng)出電磁轉(zhuǎn)矩。

根據(jù)全電流定律，沿磁力線取正向積分，有

式中H為磁場強(qiáng)度，l為磁路積分路徑，右邊表示對應(yīng)每極的勵磁磁勢。

由圖2中對閉合磁回路的分段，式（1）變?yōu)?/p>

將勵磁磁勢的計算轉(zhuǎn)化為磁路各段磁壓降的計算。其中，氣隙磁壓降為

式中，δ為氣隙徑向長度， Kδ為計及齒槽效應(yīng)引入的Carter系數(shù)，氣隙磁密最大值Bδ由每極磁通Φ、計算極弧系數(shù)α’p、電機(jī)極距τ和磁路軸向有效長度lef確定，如下式：

式中，E為相電勢，KNm為氣隙磁場波形系數(shù)，氣隙磁場呈理想正弦分布時KNm等于1.11，為2/π，電機(jī)齒部越飽和，越大。Kdp為繞組系數(shù)，f為電流頻率，Nφ為相繞組串聯(lián)匝數(shù)。KNm和α通過預(yù)取的飽和系數(shù)Ks查曲線KNm= f(Ks)和= f(Ks)得到。

圖2 雙定子電機(jī)傳統(tǒng)磁路模型

根據(jù)磁通連續(xù)原理，以內(nèi)定子為例，齒部磁密Bti、軛部磁密Bji分別為

式中，t為內(nèi)、外定子齒距，lFe為計及疊壓系數(shù)的定子鐵心軸向長度，bti近似為離齒最狹1/3齒高處的內(nèi)定子齒寬，hji為軛高。

當(dāng)齒磁密超過1.8T時，進(jìn)行如下修正：其中，bs為槽寬，Hti是查 1.8T以上磁化材料曲線的磁場強(qiáng)度。

計算出準(zhǔn)確的齒部磁密后查表獲取磁場強(qiáng)度即可算出齒部磁壓降。再將其代入飽和系數(shù)定義式

計算軛部磁壓降時，應(yīng)注意軛部磁壓降校正系數(shù)Cji的引用

式中，Cji與軛尺寸、極對數(shù)和 Bji有關(guān)，Lji為每極齒聯(lián)軛磁路計算長度。

得到各段磁壓降后，可以求出閉合回路的磁壓降和各電壓下產(chǎn)生最大磁密的勵磁電流需求，如下式

式中，Im為勵磁電流，m為相數(shù)，p為極對數(shù)。

傳統(tǒng)磁路法的核心是把氣隙段、齒部、軛部分別等效為整體考慮，計算其各自的最大磁密。雖然程序編寫較簡單，迭代收斂也比較迅速，但在實際情況中，齒部和軛部的磁勢、磁密都不是均勻的，因此引入下面的分布磁路法。

3 分布磁路法計算過程

分布磁路法是基于多個磁回路進(jìn)行迭代處理的數(shù)值計算方法，依據(jù)磁勢滿足線性疊加關(guān)系，對定子進(jìn)行等間隔周向分塊，通過氣隙磁密迭代計算出圓周氣隙中心線上各節(jié)點的磁密、磁勢及合成磁勢，進(jìn)而計算出激磁電感等。此計算方法所做的假設(shè)有：齒部磁力線為徑向、軛部磁力線為周向、引入Carter系數(shù)計及內(nèi)外定子齒槽效應(yīng)等，與傳統(tǒng)磁路法的假設(shè)一致。

3.1 磁路模型建立與劃分

在雙定子杯形轉(zhuǎn)子電機(jī)半個極距范圍內(nèi)沿周向等間隔分塊，對應(yīng)電角度從 0到 π/2，沿徑向進(jìn)行磁路計算分區(qū)：氣隙區(qū)（Ⅰ）、內(nèi)定子齒部區(qū)（Ⅱ）、內(nèi)定子軛部區(qū)（Ⅲ）、外定子齒部區(qū)（Ⅳ）、外定子軛部區(qū)（Ⅴ）。在磁路模型中，沿周向等角度均勻分為N塊，對應(yīng)氣隙中心線上N+1個等間隔節(jié)點，如圖3所示。

圖3 磁路模型徑向分區(qū)、周向分塊

3.2 階梯波磁勢表達(dá)

雙定子杯形轉(zhuǎn)子電機(jī)的勵磁磁勢為階梯波形式，在半個極距下，對 N+1個等間隔節(jié)點，設(shè)N1為內(nèi)定子、外定子齒部區(qū)域包含的節(jié)點數(shù)，N2為槽部區(qū)域包含的節(jié)點數(shù)。

設(shè)Im0為空載勵磁電流，q為每極每相線圈數(shù)，Nc為線圈匝數(shù)，q2為半個極距內(nèi)包含的內(nèi)/外定子齒數(shù)，k為計數(shù)值，取值范圍從1～q2，interp1為MATLAB線性插值函數(shù)。根據(jù)繞組布置，有半個極距下的勵磁磁動勢F(i)表達(dá)如下：

i=1～3N1/2+N2，即對應(yīng)第一個半個齒部、第一個槽部和第二個齒部時，

i=3N1/2+N2+1～3N1/2+2N2，即對應(yīng)第二個槽部時，

i=3N1/2+2N2+1～5N1/2+2N2，即對應(yīng)第三個齒部

時，

以此類推，在對應(yīng)第(k+1)個齒部，即

i=(2k-1)N1/2+kN2+1～(2k+1)N1/2+kN2時，

i=(2k+1)N1/2+kN2+1～(2k+1)N1/2+(k+1)N2時，

特別地，對應(yīng)第(q2-1)個齒部時，

在第(q2+1)個齒部，即最后半個齒部區(qū)域節(jié)點上，i=(2q2-1)N1/2+q2N2+1～N+1時，有= 0 。

3.3 各段回路磁密和磁壓降計算

給定勵磁電流Im0一個初值，同時預(yù)設(shè)一個飽和系數(shù)值kst，第i節(jié)點處氣隙磁密Bδ(i)為

在齒不太飽和時，認(rèn)為磁場全由齒中通過。即i=(2k-1)N1/2+kN2+1～(2k+1)N1/2+kN2時，第i節(jié)點處內(nèi)定子、外定子齒部磁密Bti(i)和Bto(i)為

式中，Bδav為齒部對應(yīng)節(jié)點氣隙磁密的平均值，bti和 bto分別為內(nèi)定子、外定子齒寬。其余節(jié)點的齒部磁密為0。

當(dāng)齒磁密大于1.8T時，認(rèn)為因磁路飽和磁通并未全部從齒中通過，有少部分流進(jìn)槽中，故引入類似傳統(tǒng)磁路法中的修正關(guān)系：

對于軛部磁密，同樣根據(jù)磁通連續(xù)性原理，第1節(jié)點與第i節(jié)點間氣隙中心面上的徑向磁通等于第i節(jié)點處軛部截面上的周向磁通，因此第i節(jié)點處內(nèi)定子、外定子軛部磁密Bji(i)和Bjo(i)分別為

式中，hji、hjo分別為內(nèi)定子軛高和外定子軛高。

根據(jù)第i節(jié)點處內(nèi)定子、外定子齒部和軛部磁密，查材料磁化曲線得到對應(yīng)的磁場強(qiáng)度Hti(i)、Hji(i)和 Hto(i)、Hjo(i)。因此閉合磁回路各段的磁壓降分別為

式中，hti、hto分別為內(nèi)定子、外定子齒高，Lj1、Lj2分別為內(nèi)定子、定子每極齒聯(lián)軛磁路計算長度。

經(jīng)過第i節(jié)點的閉合回路總磁壓降為

3.4 氣隙磁密迭代

由安培環(huán)路定律，閉合磁路中的總磁壓降等于磁勢。對1～N范圍內(nèi)的每一個節(jié)點對應(yīng)的閉合磁路，磁壓降應(yīng)等于該節(jié)點的勵磁磁勢。對兩個磁勢向量，認(rèn)為它們之間的相對距離d小于給定精度ε就表示磁勢相等，即

若式(32)不滿足給定精度，則對氣隙磁密Bδ(i)進(jìn)行如下修正，代回計算程序重新計算，直到滿足精度。第i節(jié)點處的修正氣隙磁密B′δ(i)為

式中，kB為修正系數(shù)，取值范圍為 0.05～0.5，一般飽和程度較低時取大值，飽和程度高時取小值，取大值可能發(fā)散，取小值迭代次數(shù)增多。

根據(jù)磁勢的奇偶對稱性將氣隙磁密值由半個極距延拓至一對極范圍，即從1/4個周期延拓到1個周期。用MATLAB中傅里葉函數(shù)分解得到氣隙磁密基波幅值，以及各次諧波分布。

4 計算結(jié)果與實驗結(jié)果比較驗證

4.1 氣隙磁密分布分析

根據(jù)前述兩種雙定子杯形轉(zhuǎn)子電機(jī)磁路計算方法，利用MATLAB軟件編寫計算程序，分別用兩種方法計算電機(jī)的氣隙磁密分布。其中，分布磁路法計算值延拓后一個周期下的氣隙磁密分布如圖4所示。

可以看出，分布磁路法計算出的氣隙磁密分布體現(xiàn)了階梯波磁勢的特征。而傳統(tǒng)磁路法利用計算極弧系數(shù)和極弧計算長度描繪氣隙磁密分布，只能粗略地描述氣隙磁密變化趨勢同計算極弧系數(shù)的關(guān)系（齒部越飽和，氣隙磁場波形越平，計算極弧系數(shù)越大），不能像分布磁路法這樣從階梯波磁勢的角度表現(xiàn)更貼近真實的磁密變化。這一點說明分布磁路法相對傳統(tǒng)磁路法計算結(jié)果更加豐富，便于后續(xù)作其他分析。

4.2 空載特性比較驗證

空載實驗由外部三相電源拖動雙定子杯形轉(zhuǎn)子電機(jī)加速到空載轉(zhuǎn)速，在空載運行穩(wěn)定后，通過Synergy采集電機(jī)轉(zhuǎn)速、內(nèi)定子和外定子三相電流、兩相線電壓數(shù)據(jù)以及電機(jī)輸入功率。取95Hz這一工作點進(jìn)行數(shù)據(jù)分析，將得到的空載特性曲線同計算結(jié)果進(jìn)行對比。由于空載磁路計算時忽略端部漏電感和槽部漏電感，而實驗所測得的端部電壓包含這兩項，因此將計算值計及定子總漏感的影響，得到傳統(tǒng)方法與分布磁路法計算結(jié)果與實驗結(jié)果對比如圖5所示。

由圖5可知，分布磁路法和傳統(tǒng)磁路法的理論計算值與實驗結(jié)果總體誤差在 4%以內(nèi)，兩種方法計算結(jié)果相差不大，說明利用這兩種方法對雙定子杯形轉(zhuǎn)子電機(jī)進(jìn)行磁路計算較為準(zhǔn)確，電機(jī)空載特性與預(yù)期結(jié)果也相差不大。

進(jìn)一步分析，可知當(dāng)飽和程度較低時，計算值與實測值吻合度較高；隨著氣隙磁密飽和加深，傳統(tǒng)磁路法計算得到的相同電壓下的勵磁電流需求較實測偏小。這一現(xiàn)象是磁路模型簡化考慮、引入修正系數(shù)進(jìn)行等效共同造成的，用分布磁路法也產(chǎn)生同樣的結(jié)果，其主要原因在于Carter系數(shù)計及定子齒槽效應(yīng)[8]。卡特系數(shù)是在不考慮飽和并假設(shè)槽型為無限深矩形槽的理想情況下利用許克變換所得的解析式，所以在飽和以及槽深有限的情況下，實際等效氣隙長度并不等于計算方法中用卡特系數(shù)計算出的等效氣隙長度。

由于分布磁路法中磁路模型按網(wǎng)格劃分，磁場描述更加細(xì)致，減少了對參數(shù)查表的依賴，所以該方法較傳統(tǒng)方法更加靈活，通用性和一般性較好。

圖5 兩種方法計算值與實驗空載特性對比

4.3 激磁電感求取與對比

由空載磁路計算得到氣隙基波磁密分布和最大值后，可由式(34)和式(35)分別算出每極基波磁通最大值Φ1m和感應(yīng)電勢基波有效值E1。

式中，Bδ1m為氣隙基波磁密最大值，Kdp1為基波繞組系數(shù)。

在給定基波相電壓有效值計算磁路的情況下，可由空載等值電路推導(dǎo)出激磁電抗，由式(36)給出。雙定子杯形轉(zhuǎn)子感應(yīng)電機(jī)理想空載運行時轉(zhuǎn)差率為 0，忽略激磁電阻和定子電阻，其空載等值電路為激磁電抗和定子漏抗串聯(lián)的一條支路。

式中，Im1為基波激磁電流有效值，U1為相電壓有效值，Xm1為激磁電抗，Xl為定子漏抗。

進(jìn)而由式(37)計算出激磁電感參數(shù)，即

表 1給出此電機(jī)在額定工況下傳統(tǒng)磁路法（方法1）和分布式磁路法（方法2）對激磁電感的計算結(jié)果與實測值及相對誤差?？梢钥闯?，兩種方法計算出的激磁電感參數(shù)相差不大，與實測值相對誤差也較小。

表1 激磁電感計算值與實測值比較

5 結(jié)語

本文分別用傳統(tǒng)磁路法和分布磁路法對雙定子杯形轉(zhuǎn)子電機(jī)展開了磁路計算，得出該電機(jī)各部分的磁密分布和空載特性曲線，并據(jù)此計算出激磁電感參數(shù)，將理論計算值與實驗值進(jìn)行了對比，結(jié)果均較為準(zhǔn)確。比較分析兩種方法，分布磁路法普遍適用性更強(qiáng)且計算結(jié)果更加豐富，在計算雙定子杯形轉(zhuǎn)子電機(jī)的磁路中具有推廣意義。

[1] 柴鳳, 崔淑梅, 程樹康, 雙定子電機(jī)發(fā)展?fàn)顩r及展望[J]. 微特電機(jī), 2001, 05: 23-26.

[2] Ross Welburn, Lira L. Ultra high torque motor System for direct drive robotics[C]. Robots 8 Conference Proceeding Applications, Detroit, Michigan, 1984: 1963-1971.

[3] Mladen V. Terzic, Dragan S. Mihic, Slobodan N. Vukosavi?. Design of high-speed, low-inertia induction machines with drag-cup rotor[J]. IEEE Transactions on Energy Conversion, 2014, 29(1) : 169-177.

[4] 黃璞, 王雙壽, 嚴(yán)關(guān)森, 空心杯轉(zhuǎn)子異步測速發(fā)電機(jī)的計算[J]. 微特電機(jī), 1976, 04: 51-62.

[5] 管少華, 程思為, 王東, 等. 用于感應(yīng)電機(jī)設(shè)計的分布磁路法快速迭代方法[J]. 電機(jī)與控制應(yīng)用, 2016, 43(6): 17-22.

[6] WANG D, WU XZ, CHEN JQ, et al. A distributed magnetic circuit approach to analysis of multiphase induction machines with non-sinusoidal supply[J]. IEEE Transactions on Energy Conversion, 2015, 30(2): 522-532.

[7] 陳世坤. 電機(jī)設(shè)計[M]. 北京: 機(jī)械工業(yè)出版社, 2002: 23-42.

[8] 王東, 吳新振, 馬偉明, 等. 非正弦供電十五相感應(yīng)電機(jī)磁路計算方法[J]. 中國電機(jī)工程學(xué)報, 2009, 29(12): 58-64.

Magnetic Circuit Calculation of Dual-stator Induction Machine with Drag-cup Rotor

Yao Jiangfan, Ouyang Bin, Zhai Xiaofei
(National Key Laboratory of Science and Technology on Vessel Integrated Power System, Naval University of Engineering，Wuhan 430033, China)

In order to analyze the performance of dual-stator induction machine, magnetic calculation is made by conventional magnetic approach and distributed magnetic circuit approach(DMCA) respectively. Magnetic flux density and magnetic potential drop in each sections of inner stator and outer stator are studied. No-load characteristic and magnetizing inductance are calculated by iterative method under given phase voltages or exciting currents and compared with experiments. The agreement between calculated results and measured ones demonstrates the accuracy and universality of DMCA.

dual-stator；induction machine；magnetic calculation；distributed magnetic circuit approach；no-load characteristic

TM346

A

1003-4862（2017）02-0048-06

2016-10-12

國家自然科學(xué)基金資助項目（51407190）。

姚江帆(1992-)，男，碩士研究生。研究方向：集成化發(fā)電技術(shù)。E-mail:cnjfyao@163.com

歐陽斌(1981-)，男，博士，副教授。研究方向：電機(jī)設(shè)計分析及控制技術(shù)。