秦 嵐，羅應(yīng)立，趙海森

(華北電力大學(xué) 電氣與電子工程學(xué)院，北京 102206)

0 引言

大功率低速直驅(qū)永磁同步電機(jī)與傳統(tǒng)帶有減速機(jī)構(gòu)的驅(qū)動(dòng)系統(tǒng)相比，具有結(jié)構(gòu)簡(jiǎn)單、體積小以及效率高等優(yōu)點(diǎn);而永磁體的存在使得其功率因數(shù)較高，同時(shí)還可以降低激磁電流產(chǎn)生的銅耗。此類電機(jī)可直接驅(qū)動(dòng)機(jī)械負(fù)載，淘汰了冗重的減速機(jī)構(gòu)，在需要大轉(zhuǎn)矩低轉(zhuǎn)速的場(chǎng)合具有獨(dú)特優(yōu)勢(shì)。然而，受極數(shù)多、每極每相槽數(shù)小等因素的影響，其諧波含量較大[1～5]且對(duì)轉(zhuǎn)矩特性產(chǎn)生較大影響。

針對(duì)低速永磁電機(jī)轉(zhuǎn)矩方面的研究，文獻(xiàn)[6]從改進(jìn)控制策略的方向出發(fā)，采用一種重復(fù)控制和PI控制相結(jié)合的方案來降低轉(zhuǎn)矩/轉(zhuǎn)速波動(dòng);文獻(xiàn)[7]研究了不同轉(zhuǎn)子槽楔的形狀和寬度對(duì)于削弱永磁電機(jī)轉(zhuǎn)矩波動(dòng)的效果;文獻(xiàn)[8]提出一系列改善低速永磁電機(jī)轉(zhuǎn)矩波動(dòng)的方法，如磁極斜極、調(diào)整極弧系數(shù)、調(diào)整磁鋼尺寸、定子斜槽等。從收集到的文獻(xiàn)看，對(duì)于轉(zhuǎn)子極靴上的槽數(shù)以及氣隙大小對(duì)轉(zhuǎn)矩波動(dòng)的影響尚缺乏研究。本文在利用時(shí)步有限元法對(duì)低速直驅(qū)永磁同步電機(jī)轉(zhuǎn)矩特性進(jìn)行研究過程中發(fā)現(xiàn)，在眾多的設(shè)計(jì)參數(shù)中，氣隙長(zhǎng)度、極靴開槽個(gè)數(shù)這兩個(gè)設(shè)計(jì)參數(shù)會(huì)對(duì)電機(jī)動(dòng)態(tài)轉(zhuǎn)矩產(chǎn)生較明顯影響，因此，本文重點(diǎn)介紹這兩個(gè)參數(shù)對(duì)動(dòng)態(tài)轉(zhuǎn)矩以及穩(wěn)態(tài)轉(zhuǎn)矩特性的影響。

1 低速直驅(qū)永磁電機(jī)特點(diǎn)

1.1 極槽配合

多極少槽型的電機(jī)為了能夠形成對(duì)稱3相繞組，其槽數(shù)和極對(duì)數(shù)，需滿足如下前提條件[5]:

式中:Q為定子槽數(shù);p為極對(duì)數(shù);k為正整數(shù);GCD()為求取最大公約數(shù) (Greatest Common Divisor，GCD)。文獻(xiàn)[5]中提出2～24極的槽極配合方式及其每極每相槽數(shù)。本文以額定轉(zhuǎn)速200 r/min的電機(jī)作為研究實(shí)例，其定子槽數(shù)為36，極數(shù)為30，每極每相槽數(shù)僅為q=2/5。

1.2 結(jié)構(gòu)特點(diǎn)

電機(jī)基本尺寸如表1所示。定子結(jié)構(gòu)與普通三相感應(yīng)電動(dòng)機(jī)相同，永磁體以切向式結(jié)構(gòu)嵌入實(shí)心轉(zhuǎn)子內(nèi)部[9]，其外側(cè)采用導(dǎo)電槽楔。表1中最后兩項(xiàng)分別為“極靴小槽數(shù)”與“氣隙長(zhǎng)度”。極靴小槽數(shù)為0～6，表示極靴表面從不開槽至開6個(gè)小槽等多個(gè)不同結(jié)構(gòu)方案; “氣隙長(zhǎng)度”為2.5～6 mm也對(duì)應(yīng)于多個(gè)不同的結(jié)構(gòu)方案。

表1 電機(jī)模型基本尺寸Tab.1 Basic size of the motor model

電機(jī)計(jì)算模型的橫截面局部結(jié)構(gòu)如圖1所示。其中1為永磁體;2是用鋁合金制成的轉(zhuǎn)子導(dǎo)電槽楔;3為低碳鋼制成的轉(zhuǎn)子鐵心;4為極靴表面小槽;5為定子槽。為了進(jìn)一步研究冷卻方案，文中定子側(cè)采用矩形開口槽方案。

圖1 模型電機(jī)局部區(qū)域結(jié)構(gòu)簡(jiǎn)圖Fig.1 Local structure diagram of the motor model

1.3 繞組排列方式

考慮到電機(jī)為多極少槽的形式，每極每相槽數(shù)小于1/2，采用槽電勢(shì)星形圖方法[10]來確定繞組排列。該電機(jī)可以劃分為3個(gè)單元電機(jī) (10極12個(gè)槽)。參考文獻(xiàn)[2]和[4]，對(duì)比不同繞組節(jié)距后確定繞組節(jié)距為1，最終繞組排列方式如圖2所示。

2 低速同步電機(jī)動(dòng)態(tài)過程的計(jì)算

2.1 時(shí)步有限元計(jì)算模型

對(duì)電機(jī)起動(dòng)等瞬變過程的分析，目前比較準(zhǔn)確的方法是時(shí)步有限元法，該方法將電機(jī)橫截面離散成有限個(gè)單元的集合，并將磁場(chǎng)、電路、運(yùn)動(dòng)方程耦合到一起，能夠充分考慮鐵磁材料的飽和、起動(dòng)過程中轉(zhuǎn)子導(dǎo)條集膚效應(yīng)等因素，通過給定電源電壓，可直接計(jì)算出磁密、定子電流、電磁轉(zhuǎn)矩、轉(zhuǎn)速等物理量隨時(shí)間變化的規(guī)律[11]。在建立電機(jī)的有限元模型時(shí)需作如下假設(shè):(1)定子鐵心外表面和轉(zhuǎn)子鐵心內(nèi)表面無漏磁;(2)忽略鐵心渦流效應(yīng);(3)端部參數(shù)均采用集中參數(shù)處理。根據(jù)上述假設(shè)，電機(jī)電磁場(chǎng)的邊值問題可以寫成如下的統(tǒng)一形式[11]:

式中:Ω為求解區(qū)域;Γ1為定子鐵心外圓邊界;Γ2為永磁體的邊界;A為矢量磁位;μ為相對(duì)磁導(dǎo)率;Js為電流密度;σ為電導(dǎo)率;v1和v1為不同材料的相對(duì)磁阻率;δs為永磁體的等效面電流密度。

電機(jī)模型的主要求解區(qū)域如圖1中的編號(hào)所示，不同求解區(qū)域的材料屬性如表2所示。

圖2 繞組連接展開圖Fig.2 Expanded windings connection diagram

表2 求解區(qū)域材料屬性Tab.2 Material properties of the solving area

永磁電動(dòng)機(jī)采用星接方式，定子側(cè)方程為[11]:

式中:帶下標(biāo)的e與u代表各對(duì)應(yīng)相的相電勢(shì)及相電壓;r與 Ls為每相電阻及漏電感。運(yùn)動(dòng)方程為[12]:

式中:J為電機(jī)及負(fù)載的轉(zhuǎn)動(dòng)慣量;Te為電磁轉(zhuǎn)矩，根據(jù)時(shí)步有限元計(jì)算所得氣隙磁場(chǎng)采用麥克斯韋應(yīng)力法[13，14]計(jì)算得到;Tm為負(fù)載轉(zhuǎn)矩;Ω 為轉(zhuǎn)子的機(jī)械角速度;θ為轉(zhuǎn)子機(jī)械角。

2.2 典型狀態(tài)的磁場(chǎng)分布

通過時(shí)步有限元計(jì)算，可以得到起動(dòng)過程及穩(wěn)態(tài)任一時(shí)刻的磁場(chǎng)分布，進(jìn)而求出轉(zhuǎn)矩。堵轉(zhuǎn)情況下的磁力線分布及穩(wěn)態(tài)時(shí)齒頂中心一個(gè)點(diǎn)的磁密隨時(shí)間變化的曲線，如圖3所示。

單元電機(jī)中空載氣隙磁密的空間分布曲線如圖4(a)所示，其諧波分析結(jié)果如圖4(b)所示?？梢钥闯?，單元電機(jī)中永磁體磁極所產(chǎn)生的空載氣隙磁密中的基波幅值最大，是所需要的工作波。

3 主要設(shè)計(jì)參數(shù)對(duì)動(dòng)態(tài)轉(zhuǎn)矩的影響

3.1 氣隙長(zhǎng)度對(duì)動(dòng)態(tài)轉(zhuǎn)矩的影響

進(jìn)行對(duì)比分析計(jì)算的氣隙分別為2.5 mm，4 mm與6 mm，其極靴表面槽數(shù)均為6個(gè)，槽內(nèi)的導(dǎo)條為紫銅。不同氣隙下起動(dòng)過程對(duì)比如圖5所示。

定轉(zhuǎn)子磁場(chǎng)極對(duì)數(shù)相同，轉(zhuǎn)速相同，磁場(chǎng)相互作用產(chǎn)生恒定力矩;定轉(zhuǎn)子磁場(chǎng)極對(duì)數(shù)相同，轉(zhuǎn)速不同，相互作用產(chǎn)生波動(dòng)力矩[13]。由于該電機(jī)中諧波磁場(chǎng)較大，導(dǎo)致穩(wěn)態(tài)轉(zhuǎn)矩存在不同周期和幅值的波動(dòng)。由此可以看出氣隙長(zhǎng)度為2.5 mm，4 mm以及6 mm時(shí)，電機(jī)穩(wěn)態(tài)轉(zhuǎn)矩波動(dòng)幅值分別約為161 N·m，120 N·m以及57 N·m。

圖5 不同氣隙長(zhǎng)度的轉(zhuǎn)速、轉(zhuǎn)矩對(duì)比Fig.5 Speed and torque comparison of different air gap length

3.2 轉(zhuǎn)子極靴表面開槽對(duì)動(dòng)態(tài)轉(zhuǎn)矩的影響

轉(zhuǎn)子槽的多少?zèng)Q定了齒諧波磁勢(shì)的大小。齒諧波的次數(shù)為Q/p±1，使齒諧波次數(shù)增高，其影響就小。針對(duì)這個(gè)問題，考慮增加轉(zhuǎn)子槽數(shù)Q2，方法是在轉(zhuǎn)子極靴上均勻開小槽，槽中嵌入紫銅導(dǎo)條。固定氣隙為4 mm，對(duì)于轉(zhuǎn)子極靴表面不開槽、開1個(gè)小槽、均勻開6個(gè)小槽三種情況，對(duì)其動(dòng)態(tài)轉(zhuǎn)矩進(jìn)行對(duì)比。極靴表面不同個(gè)數(shù)小槽計(jì)算區(qū)域如圖6所示，所開小槽寬度為2 mm，深度為10 mm。起動(dòng)過程的轉(zhuǎn)矩以及穩(wěn)態(tài)轉(zhuǎn)矩曲線對(duì)比如圖7所示。

由圖7可見，轉(zhuǎn)子極靴表面如果不開小槽，則穩(wěn)態(tài)轉(zhuǎn)矩波動(dòng)幅值約為341 N·m，開1個(gè)小槽，則穩(wěn)態(tài)轉(zhuǎn)矩波動(dòng)幅值約為159 N·m，一旦轉(zhuǎn)子極靴表面小槽數(shù)達(dá)到6個(gè)，則穩(wěn)態(tài)轉(zhuǎn)矩波動(dòng)幅值僅為120 N·m。

3.3 時(shí)步有限元計(jì)算結(jié)果可信度分析

為了考核文中所用時(shí)步有限元程序的正確性，課題組前期已對(duì)現(xiàn)有22 kW自起動(dòng)永磁電機(jī)起動(dòng)及運(yùn)行性能進(jìn)行了實(shí)測(cè)和計(jì)算對(duì)比，兩者基本吻合[12～14]。間接驗(yàn)證了本文時(shí)步有限元計(jì)算結(jié)果的可信性。

4 結(jié)論

(1)穩(wěn)態(tài)轉(zhuǎn)矩波動(dòng)量隨氣隙增加而顯著減小，當(dāng)固定極靴表面開槽個(gè)數(shù)為6個(gè)，氣隙為2.5 mm，4 mm及6 mm時(shí)，空載穩(wěn)態(tài)轉(zhuǎn)矩波動(dòng)量分別約為161 N·m，120 N·m以及57 N·m。

(2)極靴表面是否開槽及開槽個(gè)數(shù)對(duì)穩(wěn)態(tài)轉(zhuǎn)矩波動(dòng)量影響很大。對(duì)于氣隙4 mm的情況，如果極靴表面不開槽，空載轉(zhuǎn)矩波動(dòng)量約為341 N·m，開1個(gè)槽，則該波動(dòng)量約為159 N·m，如果槽數(shù)達(dá)到6個(gè)，則空載轉(zhuǎn)矩波動(dòng)量?jī)H為120 N·m。

此外，經(jīng)計(jì)算發(fā)現(xiàn)轉(zhuǎn)子槽楔導(dǎo)條材料、極靴槽內(nèi)導(dǎo)條材料以及定子槽深同樣會(huì)對(duì)轉(zhuǎn)矩產(chǎn)生影響，上述因素對(duì)轉(zhuǎn)矩的影響有待進(jìn)一步研究。

[1]Freddy Magnussen，Heinz Lendenmann.Parasitic effects in PM machines with concentrated windings[J].IEEE Transon Industry Applications， 2007， 43(5):1223-1232.

[2]Nicola Bianchi，Silverio Bolognani，Michele Dai Pré，et al.Design considerations for fractional-slot winding configurations of synchronous machines[J].IEEE Trans on Industry Applications，2006，42(4):997-1006.

[3]Ayman M.EL-Refaie.Fractional-Slot Concentrated-Windings Synchronous Permanent Magnet Machines:Opportunities and Challenges[J].IEEE Trans on Industrial Electronics，2010，57(1):107-121.

[4]Mihail V Cistelecan1，F(xiàn)ernando J T E Ferreira，Mihail Popescu1.Three phase tooth-concentrated multiple-layer fractional windings with low space harmonic content[C].IEEE Energy Conversion Congress and Exposition(ECCE)，2010.1399-1405.

[5]Jér?me Cros，Philippe Viarouge，Synthesis of High Performance PM Motors With Concentrated Windings[J].IEEE Trans on Energy Conversion，2002，17(2):248-253.

[6]儲(chǔ)劍波，胡育文，黃文新，等.一種抑制永磁同步電機(jī)轉(zhuǎn)速脈動(dòng)的方法[J].電工技術(shù)學(xué)報(bào)，2009，24(12):43-49.Chu Jianbo，Hu Yuwen，Huang Wenxin，et al.Suppressing speed ripples of permanent magnetic synchronous motor based on a Method[J].Transactions of China Electrotechnical Society，2009，24(12):43-49.

[7]楊玉波，王秀和，丁婷婷.一種削弱永磁同步電動(dòng)機(jī)齒槽轉(zhuǎn)矩的方法[J].電機(jī)與控制學(xué)報(bào)，2008，12(5):520-523.Yang Yubo，Wang Xiuhe，Ding Tingting.Method of reducing cogging torque of solid-rotor permanent magnet synchronous motors[J].Electric Machines and Control，2008，12(5):520-523.

[8]肖杭，金敏捷，譚娃，等.低速永磁同步電動(dòng)機(jī)的設(shè)計(jì)研究[J]，微特電機(jī)，2001，(6):26-28.Xiao Hang，Jin Minjie，Tan Wa，et al.Design study of permanent magnet synchronous motor with low speed[J].Small＆ Special Machines，2001，(6):26-28.

[9]唐任遠(yuǎn).現(xiàn)代永磁電機(jī)理論與設(shè)計(jì)[M].北京:機(jī)械工業(yè)出版社，1997.

[10]湯蘊(yùn)璆，史乃，沈文豹.電機(jī)理論與運(yùn)行[M].北京:水利電力出版社，1984.

[11]李志強(qiáng).自起動(dòng)永磁電機(jī)動(dòng)態(tài)過程的有限元分析與實(shí)驗(yàn)研究[D].北京:華北電力大學(xué)，2009.

[12]羅應(yīng)立，李志強(qiáng)，趙海森，等.自起動(dòng)稀土永磁電機(jī)斷相運(yùn)行時(shí)步有限元分析[J].電機(jī)與控制學(xué)報(bào)，2009，13(3):383-388.Luo Yingli，Li Zhiqiang，zhao Haisen，et al.Openphase operation analysis of line-start rare-earth permanent magnet motor with time-stepping FEM[J].Electric Machines and Control，2009，13(3):383-388.

[13]付媛，羅應(yīng)立，陳偉華，等.自起動(dòng)永磁電動(dòng)機(jī)的穩(wěn)態(tài)轉(zhuǎn)矩波動(dòng)研究[J].電機(jī)與控制學(xué)報(bào)，2010，14(11):1-6.Fu Yuan;Luo Yingli;Chen Weihua，et al.Research on torque fluctuation of line-start PM motor based on time-step finite element method[J].Electric Machines and Control，2010，14(11):1-6.

[14]姚有光，湯蘊(yùn)璆.用麥克斯韋應(yīng)力法計(jì)算凸極同步電機(jī)的電磁轉(zhuǎn)矩[J].電機(jī)與控制學(xué)報(bào)，1989，12(1):20-25.Yao Youguang，Tang Yunqiu.Numerical calculation of tangential forces and electromagnetic torque of salient pole synchronous machine by maxwell stress method[J].Electric Machines and Control，1989，12(1):20-25.