李 波,肖林京,孫傳余,文藝成,肖 楠,陳忠霞
(山東科技大學(xué),青島 266590)
?
雙相勵(lì)磁8/6極開關(guān)磁阻電動(dòng)機(jī)轉(zhuǎn)矩?cái)?shù)學(xué)模型
李 波,肖林京,孫傳余,文藝成,肖 楠,陳忠霞
(山東科技大學(xué),青島 266590)
為研究雙相勵(lì)磁8/6極開關(guān)磁阻電動(dòng)機(jī)轉(zhuǎn)矩特性,文章根據(jù)電機(jī)原理和等效磁路法,給出了雙繞組間自感、互感公式;利用開關(guān)磁阻電動(dòng)機(jī)磁路分割法,求解了8/6極開關(guān)磁阻電動(dòng)機(jī)的氣隙磁導(dǎo),進(jìn)而推導(dǎo)出電機(jī)的轉(zhuǎn)矩模型公式。最后,通過有限元仿真結(jié)果,證明了轉(zhuǎn)矩?cái)?shù)學(xué)模型的正確性。研究表明:該轉(zhuǎn)矩?cái)?shù)學(xué)模型能正確描述8/6極開關(guān)磁阻電動(dòng)機(jī)無位置偏移時(shí)的轉(zhuǎn)矩輸出特性。
開關(guān)磁阻電動(dòng)機(jī);雙相勵(lì)磁;數(shù)學(xué)模型;有限元分析
開關(guān)磁阻電動(dòng)機(jī)[1-2](以下簡稱SRM),作為一種新型電機(jī),因結(jié)構(gòu)簡單、工作可靠、維護(hù)量小、啟動(dòng)性能及調(diào)速性能優(yōu)越等優(yōu)點(diǎn),受到了研究學(xué)者的廣泛關(guān)注。
8/6極SRM從勵(lì)磁方式可分為單相勵(lì)磁和雙相勵(lì)磁,文獻(xiàn)[3-4]針對(duì)上述兩種不同的勵(lì)磁方式,對(duì)其自感、互感及轉(zhuǎn)矩輸出特性進(jìn)行了理論及實(shí)驗(yàn)分析,研究表明雙相勵(lì)磁模式較單相勵(lì)磁模式,具有平均轉(zhuǎn)矩高、輸出平穩(wěn)性好的優(yōu)點(diǎn)。
轉(zhuǎn)矩性能是SRM領(lǐng)域的研究重點(diǎn),對(duì)SRM調(diào)節(jié)轉(zhuǎn)速、穩(wěn)定旋轉(zhuǎn)等具有重要意義,文獻(xiàn)[5]針對(duì)8/6極SRM,在建立樣機(jī)仿真模型,獲取并分析仿真數(shù)據(jù)的基礎(chǔ)上,給出了自互感模型、磁鏈模型以及轉(zhuǎn)矩模型,但上述模型為基于該樣機(jī)仿真數(shù)據(jù)的擬合模型,并未給出關(guān)于電機(jī)具體參數(shù)的數(shù)學(xué)公式。文獻(xiàn)[6]針對(duì)8/6極單繞組磁懸浮SRM,借鑒文獻(xiàn)[7-10]的磁路分割法,通過建立等效磁路模型、求解氣隙磁導(dǎo),推導(dǎo)出了電機(jī)轉(zhuǎn)矩模型。由于8/6極磁懸浮SRM和8/6極SRM具有相同的定轉(zhuǎn)子齒極結(jié)構(gòu),即而借鑒8/6極磁懸浮SRM的氣隙磁導(dǎo)求解方法,并應(yīng)用到8/6極SRM氣隙磁導(dǎo)的計(jì)算研究中。
文章針對(duì)雙相勵(lì)磁8/6極SRM,在Maxwell 2D仿真分析的基礎(chǔ)上,建立了等效磁路模型,并推導(dǎo)出自互感表達(dá)式;再利用磁路分割法給出了齒間氣隙磁導(dǎo)求解方法及磁導(dǎo)公式,藉此推導(dǎo)出雙相勵(lì)磁8/6極SRM的轉(zhuǎn)矩?cái)?shù)學(xué)模型;最后通過有限元仿真對(duì)比,驗(yàn)證了該數(shù)學(xué)模型的正確性。
雙相勵(lì)磁8/6極SRM,包括電機(jī)定子、電機(jī)轉(zhuǎn)子和勵(lì)磁繞組,且兩相繞組共同通電,為電機(jī)轉(zhuǎn)子提供電磁吸引力。如圖1所示,AB相勵(lì)磁繞組通電后,在定轉(zhuǎn)子齒間形成電磁通路,產(chǎn)生電磁吸引力,驅(qū)動(dòng)電機(jī)轉(zhuǎn)子逆時(shí)針轉(zhuǎn)動(dòng)。
定義逆時(shí)針旋轉(zhuǎn)方向?yàn)檎?,定轉(zhuǎn)子齒極正對(duì)時(shí),轉(zhuǎn)子位置角為0,則AB相導(dǎo)通時(shí),A相轉(zhuǎn)子位置角θa∈[-π/12,0],B相轉(zhuǎn)子位置角θb∈[-π/6,-π/12],D相轉(zhuǎn)子位置角θd∈[0,π/12],且有關(guān)系θb=θa-π/12,θd=θa+π/12,圖1中所示為θa=-π/24的位置。

圖1 雙相勵(lì)磁8/6極SRM原理圖
本文采用Maxwell 2D有限元仿真軟件,建立雙相勵(lì)磁8/6極SRM的二維模型,選擇材料,定義邊界條件,添加繞組激勵(lì),設(shè)置運(yùn)動(dòng)參量,最后完成瞬態(tài)磁場的仿真分析。其中,二維電機(jī)模型參數(shù)設(shè)置如下:定子鐵心外半徑30 mm,內(nèi)半徑19.1 mm,定子軛環(huán)厚度4.9 mm,定子齒極高6 mm,定子齒極寬7.45 mm;轉(zhuǎn)子鐵心外半徑19 mm,內(nèi)半徑5.36 mm,轉(zhuǎn)子軛環(huán)厚度5 mm,轉(zhuǎn)子齒極高8.64 mm,轉(zhuǎn)子齒極寬7.41 mm;軸向長度1 000 mm,各相繞組匝數(shù)50。
Maxwell 2D針對(duì)雙相勵(lì)磁8/6極SRM,進(jìn)行磁路仿真分析,當(dāng)AB相繞組導(dǎo)通且 時(shí),其磁路與磁密分布分別如圖2、圖3所示。

圖2 電機(jī)磁路分布圖圖3 電機(jī)磁密分布圖
由圖2知,AB相繞組磁力線主要經(jīng)由AB相相鄰定轉(zhuǎn)子齒極及AD相相鄰定轉(zhuǎn)子齒極,形成閉合磁路,少量磁力線經(jīng)由AB相相對(duì)定轉(zhuǎn)子齒極和AC相定轉(zhuǎn)子齒極,形成閉合磁路。
由圖3知,AB相磁力線主要集中于A相、B相、D相定轉(zhuǎn)子齒極處,極少量分布在C相定轉(zhuǎn)子齒極處。且仿真發(fā)現(xiàn),隨著θa減小,C相定轉(zhuǎn)子齒極處磁密會(huì)有所增加,但相較于A相、B相、D相仍可忽略不計(jì)。
為簡化雙相勵(lì)磁8/6極SRM轉(zhuǎn)矩求解,忽略漏磁通、定子交鏈轉(zhuǎn)子軛部磁通、磁鋼磁導(dǎo)和磁路飽和現(xiàn)象,并假設(shè)SRM電機(jī)轉(zhuǎn)子處于平衡位置。
以8/6極SRM的AB相繞組勵(lì)磁為例,認(rèn)為AB相通電過程中,磁力線僅經(jīng)由A、B、D三相支路導(dǎo)通,即忽略C相定轉(zhuǎn)子齒間氣隙磁導(dǎo),并繪制其等效磁路圖,如圖4所示:Na為A相繞組匝數(shù);ia為A相繞組電流;Nb為B相繞組匝數(shù);ib為B相繞組電流;Pa1,Pa2和φa1,φa2分別為A相上、下兩個(gè)氣隙處的氣隙磁導(dǎo)和氣隙磁通量;Pb1,Pb2和φb1,φb2分別為B相上、下兩個(gè)氣隙處的氣隙磁導(dǎo)和氣隙磁通量;Pd1,Pd2與φd1,φd2分別為D相上、下兩個(gè)氣隙處的氣隙磁導(dǎo)和氣隙磁通量。

圖4 AB相等效磁路圖
由電機(jī)齒極結(jié)構(gòu)的對(duì)稱性,可知Pa1=Pa2,Pb1=Pb2,Pd1=Pd2,且φa1=φa2,φb1=φb2,φd1=φd2,根據(jù)AB相等效磁路圖,可得AB相繞組自感、互感表達(dá)式:
(1)
式中:La為A相繞組自感;Lb為B相繞組自感;Mab為AB相繞組互感;P0=Pa+Pb+Pd;Pa=Pa1=Pa2;Pb=Pb1=Pb2;Pd=Pd1=Pd2。
根據(jù)文獻(xiàn)[6-10],求解氣隙磁導(dǎo)時(shí),利用磁路分割法,將齒間磁路分割為正對(duì)部分的直線磁路和非正對(duì)部分的橢圓磁路。

圖5 磁路分割圖
當(dāng)SRM轉(zhuǎn)子無位置偏移時(shí),其氣隙磁導(dǎo)[6]:
(2)

根據(jù)雙相勵(lì)磁8/6極SRM原理可知,A相轉(zhuǎn)子位置角θa∈[-π/12,0],帶入式(2),得A相齒間氣隙磁導(dǎo)Pa:
(3)
同理可知,D相轉(zhuǎn)子位置角θd∈[0,π/12],且有θd=θa+π/12,則D相齒間氣隙磁導(dǎo)Pd:
(4)
對(duì)于B相齒間氣隙磁導(dǎo) ,已知其轉(zhuǎn)子位置角θb∈[-π/6,-π/12],且θb=θa-π/12。如圖1所示,當(dāng)θb∈[-π/8,-π/12](即θa∈[-π/24,0])時(shí),B相定轉(zhuǎn)子齒極既有正對(duì)部分,又有非正對(duì)部分,與A、B相同理求解,即有:
(5)
當(dāng)θb∈[-π/6,-π/8](即θa∈[-π/12,-π/24])時(shí),B相轉(zhuǎn)子齒與B相定子齒完全偏離,僅有非正對(duì)部分氣隙磁導(dǎo),如圖6所示,此時(shí)氣隙磁導(dǎo)式(2)不再適用,Pb求解方式如下。

圖6 B相氣隙磁路示意圖
根據(jù)文獻(xiàn)[7],非正對(duì)部分微元磁導(dǎo)公式:
(6)
式中:dP2為非正對(duì)部分微元磁導(dǎo);dt為非正對(duì)部分微元長度。
首先假設(shè):(1) B相轉(zhuǎn)子齒與B相定子齒完全偏離時(shí),非正對(duì)部分微元磁導(dǎo)仍符合式(6);(2) 忽略B相轉(zhuǎn)子齒與C相定子齒間的氣隙磁導(dǎo)。根據(jù)圖6可知,非正對(duì)部分微元磁導(dǎo)dP2的積分長度t∈[(|θb|-π/8)r,π/8·r]。
將積分長度上下限,帶入式(6),并化簡得:
(7)
整合式(5)和式(7),得:
(8)
已知雙相勵(lì)磁8/6極SRM,AB相繞組儲(chǔ)能式Wab:
(9)
得雙相勵(lì)磁8/6極SRM轉(zhuǎn)矩公式T(θa):
(10)
聯(lián)立式(1)、式(3)、式(4)、式(8)、式(10),并化簡,得電機(jī)轉(zhuǎn)矩?cái)?shù)學(xué)模型:
(11)
式中:
f0=fa+fb+fd。
以雙相勵(lì)磁8/6極SRM為仿真對(duì)象,利用Maxwell2D有限元仿真軟件,得到θa∈[-15°,0°]過程中,電機(jī)轉(zhuǎn)矩與θa間的仿真輸出轉(zhuǎn)矩曲線。另,依據(jù)轉(zhuǎn)矩?cái)?shù)學(xué)模型式(11),利用MATLAB軟件,繪制轉(zhuǎn)矩T(θa)與θa間的模型計(jì)算轉(zhuǎn)矩曲線。通過對(duì)比仿真輸出轉(zhuǎn)矩與模型計(jì)算轉(zhuǎn)矩的擬合效果,驗(yàn)證轉(zhuǎn)矩?cái)?shù)學(xué)模型的正確性。
雙相勵(lì)磁8/6極SRM仿真模型,分別設(shè)定AB相勵(lì)磁繞組電流ia=ib=2A,ia=ib=3A,進(jìn)行瞬態(tài)仿真,則仿真輸出轉(zhuǎn)矩與模型計(jì)算轉(zhuǎn)矩的對(duì)比情況如圖7所示。由圖7對(duì)比可知,仿真輸出轉(zhuǎn)矩與模型計(jì)算轉(zhuǎn)矩?cái)M合效果良好,即轉(zhuǎn)矩?cái)?shù)學(xué)模型式(11)能較準(zhǔn)確地描述雙相勵(lì)磁8/6極SRM導(dǎo)通周期內(nèi)的轉(zhuǎn)矩輸出情況。

圖7 仿真輸出轉(zhuǎn)矩與模型計(jì)算轉(zhuǎn)矩對(duì)比圖
但曲線繪制過程中發(fā)現(xiàn),模型計(jì)算轉(zhuǎn)矩曲線為非連續(xù)曲線,間斷點(diǎn)為[-7.5°,T(-7.5°)],原因在于:氣隙磁導(dǎo)Pb在θa=-7.5°處,連續(xù)不可導(dǎo),而轉(zhuǎn)矩T(θa)求解時(shí),需對(duì)Pb求導(dǎo),繼而導(dǎo)致了轉(zhuǎn)矩T(θa)在該點(diǎn)處不連續(xù)。由圖7中亦可看出,由間斷采樣數(shù)據(jù)繪制出的模型計(jì)算轉(zhuǎn)矩曲線在[-7.5°,T(-7.5°)]點(diǎn)處并不圓滑。
文章通過研究雙相勵(lì)磁8/6極SRM轉(zhuǎn)矩?cái)?shù)學(xué)模型,得出如下結(jié)論:
(1) 轉(zhuǎn)矩?cái)?shù)學(xué)模型式(11)能較準(zhǔn)確地描述雙相勵(lì)磁8/6極開關(guān)磁阻電機(jī)的轉(zhuǎn)矩輸出情況;(2) 式(11)為非連續(xù)函數(shù),間斷點(diǎn)為[-7.5°,T(-7.5°)],由此導(dǎo)致,當(dāng)轉(zhuǎn)子位置角θa趨近于-7.5°時(shí),實(shí)際輸出轉(zhuǎn)矩與公式計(jì)算可能存在一定偏差。
[1] 吳建華.開關(guān)磁阻電機(jī)設(shè)計(jì)與應(yīng)用[M].北京:機(jī)械工業(yè)出版社,2000.
[2] 甘醇,吳建華,王寧,等.一種零電壓保持開通的開關(guān)磁阻電機(jī)再生制動(dòng)控制策略[J].電機(jī)與控制學(xué)報(bào),2015,19(9):8-15.
[3] 高潔,張鶴旭,米彥青.計(jì)及互感的開關(guān)磁阻電機(jī)單雙相勵(lì)磁靜態(tài)性能分析[J].電機(jī)與控制學(xué)報(bào),2012,16(11):45-51.
[4] 弓鏘,曲爾光,劉艷芳.8/6極開關(guān)磁阻電機(jī)兩相勵(lì)磁電磁轉(zhuǎn)矩分析[J].現(xiàn)代電子技術(shù),2012,35(12):176-178.
[5] 林委.計(jì)及互感的8/6極開關(guān)磁阻電機(jī)轉(zhuǎn)矩估算模型研究[D].重慶:重慶大學(xué),2014.
[6] 嵇小輔,馬滔,項(xiàng)倩雯.8/6極單繞組磁懸浮開關(guān)磁阻電機(jī)建模與分析[J].微電機(jī),2014,47(6):9-14.
[7] TAKEMOTO M,SUZUKI H,CHIBA A,et al.Improved analysis of a bearingless switched reluctance motor[J].IEEE Transactions on Industry Applications,2001,37(1):26-34.
[8] 張京軍,龍榮,張海軍,等.開關(guān)磁阻電機(jī)徑向電磁力解析建模及有限元分析[J].煤炭學(xué)報(bào),2012,37(4):700-704.
[9] 孫玉坤,吳建兵,項(xiàng)倩雯.基于有限元法的磁懸浮開關(guān)磁阻電機(jī)數(shù)學(xué)模型[J].中國電機(jī)工程學(xué)報(bào),2007,27(12):33-40.
[10] 楊艷,鄧智泉,曹鑫,等.無軸承開關(guān)磁阻電機(jī)徑向電磁力模型[J].電機(jī)與控制學(xué)報(bào),2009,13(3):377-382.
Torque Mathematic Model of 8/6 Switched Reluctance Motor with Two-Winding Excitation
LIBo,XIAOLin-jing,SUNChuan-yu,WENYi-cheng,XIAONan,CHENZhong-xia
(Shandong University of Science and Technology,Qingdao 266590,China)
In view of the torque of 8/6 switched reluctance motor (SRM) with two-winding excitation, the stator winding's self and mutual inductance were deduced based on the principle of 8/6 SRM and the method of magnetic equivalent circuit. Utilizing the method of magnetic field division, the gap permeance of 8/6 SRM was solved, and then the mathematic model of torque was derived. Finally, based on the results of finite-element analysis, the validity of the mathematic model was verified. The studies show that the mathematic model of torque can describe the output characteristics of 8/6 SRM correctly when the motor stays balanced.
switched reluctance motor; two-phase excitation; mathematic model; finite-element analysis
2016-01-13
山東省優(yōu)秀中青年科學(xué)家科研獎(jiǎng)勵(lì)基金(BS2013NJ015);青島市博士后研究人員應(yīng)用研究項(xiàng)目(01020120521);青島市應(yīng)用基礎(chǔ)研究計(jì)劃項(xiàng)目(15-9-1-66-jch)
TM352
A
1004-7018(2016)08-0019-03
李波(1991-),男,碩士研究生,研究方向?yàn)榇艖腋〖夹g(shù)。