王 偉,王淑紅,梁力波,張一博
(太原理工大學(xué) 電氣與動(dòng)力工程學(xué)院,太原 030024)
隨著國(guó)家節(jié)能政策的推廣執(zhí)行,永磁同步電動(dòng)機(jī)由于體積小,質(zhì)量輕,效率高,運(yùn)行安全可靠,其研究和推廣不斷升級(jí),在航空航天、工業(yè)生產(chǎn)、軍用機(jī)械、新能源汽車(chē)及風(fēng)電等領(lǐng)域得到了普遍應(yīng)用[1-2]。為了提高內(nèi)置式永磁同步電機(jī)(Interior Permanent Magnet Synchronous Motor, IPMSM)的帶載能力和調(diào)速范圍,國(guó)內(nèi)外學(xué)者做了一些研究。文獻(xiàn)[3]將迭代法應(yīng)用到控制系統(tǒng)中,得到需要的d、q軸電流給定,具有一定的工程實(shí)用性,但其過(guò)程比較復(fù)雜;文獻(xiàn)[4]先離散測(cè)試電機(jī)各個(gè)狀態(tài)下的參數(shù),然后制作成表格,可以實(shí)時(shí)根據(jù)工況查詢所需電流給定值,但是該方法需要大量存儲(chǔ)空間;文獻(xiàn)[5]在轉(zhuǎn)矩給定的情況下,利用定步長(zhǎng)改變電流矢量角從而改變電流軌跡,但是在變負(fù)載時(shí),系統(tǒng)穩(wěn)定性很難保證。
本文從實(shí)驗(yàn)室一臺(tái)內(nèi)置式永磁同步電機(jī)自身的參數(shù)出發(fā),針對(duì)MTPA控制策略中d、q軸電流給定值無(wú)法獲得解析解的問(wèn)題,提出了一種多項(xiàng)式曲線擬合方法求解高次方程組,得到了不同轉(zhuǎn)速要求時(shí)的交直軸電流給定值,同時(shí)采用梯度下降法,利用轉(zhuǎn)矩梯度和電壓差對(duì)MTPA指令進(jìn)行修正,解決了電機(jī)id=0矢量控制下無(wú)法額定運(yùn)行的問(wèn)題,實(shí)現(xiàn)了電機(jī)全范圍調(diào)速。仿真和實(shí)驗(yàn)表明,在本文控制算法下,電機(jī)控制精度高,具有較寬的調(diào)速范圍和較強(qiáng)的帶載能力。
永磁同步電機(jī)是一個(gè)多變量、強(qiáng)耦合的非線性高階系統(tǒng)。本文基于轉(zhuǎn)子磁場(chǎng)定向理論,將非線性的電磁關(guān)系線性化,建立了PMSM的數(shù)學(xué)模型。首先將一些比較小的參數(shù)和因素省去,做出以下假設(shè)[6]:
(1)永磁體產(chǎn)生主磁場(chǎng),轉(zhuǎn)子上無(wú)阻尼繞組。
(2)電機(jī)定轉(zhuǎn)子表面光滑,氣隙均勻,反電勢(shì)是正弦波。
(3)電機(jī)三相繞組對(duì)稱(chēng),空間互差120°,轉(zhuǎn)子結(jié)構(gòu)中心對(duì)稱(chēng)。
(4)不計(jì)磁路飽和、磁滯、渦流對(duì)電機(jī)的影響,將鐵磁部分磁路視為線性。
(5)不考慮電機(jī)參數(shù)的變化。
依據(jù)以上假設(shè),可得d-q軸旋轉(zhuǎn)坐標(biāo)系下PMSM的數(shù)學(xué)模型為[7]
定子電壓方程:
(1)
(2)
定子磁鏈方程:
ψd=Ldid+ψf
(3)
ψq=Lqiq
(4)
電磁轉(zhuǎn)矩方程:
(5)
機(jī)械運(yùn)動(dòng)方程:
(6)
式中,ud、uq分別為定子電壓d軸、q軸分量;id、iq分別為定子電流d軸、q軸分量;ψd、ψq分別為定子磁鏈d軸、q軸分量;Ld、Lq分別為定子繞組自感d軸、q軸分量;Te為電磁轉(zhuǎn)矩,TL為負(fù)載轉(zhuǎn)矩;ωm為轉(zhuǎn)子機(jī)械角速度;ωe=np*ωm為轉(zhuǎn)子電角速度;ψf為永磁體磁鏈;Rs為定子電阻;J為電機(jī)轉(zhuǎn)動(dòng)慣量;B為粘質(zhì)摩擦系數(shù);np為電機(jī)極對(duì)數(shù)。
逆變器作為電機(jī)供電設(shè)備,其輸出能力影響著電機(jī)端電壓[8-9]。電機(jī)在本體設(shè)計(jì)時(shí),其結(jié)構(gòu)和絕緣決定了電機(jī)的電流承受能力。弱磁控制時(shí),需要考慮到電壓和電流兩方面的約束。永磁同步電動(dòng)機(jī)電流、電壓、轉(zhuǎn)矩軌跡約束曲線如圖1所示。
圖1 弱磁控制下電壓、電流、轉(zhuǎn)矩軌跡
(1)電壓極限方程;
(Lqiq)2+(Ldid+ψf)≤(Ulim/ωe)2
(7)
(2)電流極限方程:
IPMSM的電流矢量幅值受到電機(jī)額定電流限制,其表達(dá)式為
在dq坐標(biāo)平面內(nèi),式(8)的函數(shù)關(guān)系是以原點(diǎn)為圓心、以ilim為半徑的圓形曲線,一般稱(chēng)之為電流極限圓,電流矢量同樣被約束在這個(gè)圓形區(qū)域內(nèi)。
綜上,考慮到去磁電流上限,IPMSM在調(diào)速過(guò)程中,電壓極限橢圓和電流極限圓的交集OABC,就是電機(jī)穩(wěn)定運(yùn)行區(qū)域。
忽略飽和因素帶來(lái)的影響,實(shí)驗(yàn)室用IPMSM樣機(jī)參數(shù)如表1所示。
表1 樣機(jī)參數(shù)
采用SVPWM控制策略,逆變器為電壓型,直流電壓為537 V,忽略逆變器側(cè)IGBT管壓降,逆變器輸出相電壓基波峰值最大為311 V,控制系統(tǒng)中坐標(biāo)變換采用恒幅值變換方式,則合成電壓矢量幅值和相電壓幅值相等。利用電機(jī)參數(shù)和隱函數(shù)畫(huà)圖指令可以在dq軸平面研究電機(jī)的運(yùn)行范圍。
圖2 矢量控制分析
圖3 MTPA控制分析
圖4 弱磁控制分析
圖2~圖4分別為不同控制方式下電機(jī)運(yùn)行狀態(tài)分析圖,其中橢圓曲線為電壓極限橢圓,正圓曲線為電流極限圓,DF曲線和OH曲線為MTPA軌跡曲線,當(dāng)轉(zhuǎn)矩一定時(shí),電機(jī)電流運(yùn)行點(diǎn)必須在紅色圓和黑色橢圓內(nèi)。圖2為id=0矢量控制下的分析圖。因?yàn)椴捎胕d=0矢量控制,電流被限制在過(guò)AC的直線上。A點(diǎn)表明電機(jī)在額定轉(zhuǎn)速下,帶載最大為3.99 N·m,C點(diǎn)表明電機(jī)在額定轉(zhuǎn)矩下最大轉(zhuǎn)速為983 r/min,說(shuō)明id=0的矢量控制對(duì)這種電機(jī)不再適用。圖3為MTPA控制下的分析圖。電機(jī)電流運(yùn)行軌跡被限制在過(guò)EFG的曲線上,E點(diǎn)表明電機(jī)在額定轉(zhuǎn)速下,帶載最大為5.962 N·m,E點(diǎn)表明電機(jī)在額定轉(zhuǎn)矩下最大轉(zhuǎn)速為1235 r/min。比較圖2和圖3,MTPA控制策略明顯提高了電機(jī)的帶載能力和調(diào)速范圍。盡管采用了MTPA的控制方式,電機(jī)依舊無(wú)法到達(dá)額定運(yùn)行狀態(tài),為了進(jìn)一步提高電機(jī)的轉(zhuǎn)速和轉(zhuǎn)矩,電機(jī)的電流穩(wěn)態(tài)運(yùn)行點(diǎn)必須向d軸負(fù)向偏移,如圖4所示。本文采用梯度下降法進(jìn)行弱磁,對(duì)MTPA曲線上的電流指令進(jìn)行修正,使得電機(jī)運(yùn)行范圍增大,額定負(fù)載時(shí)電機(jī)的轉(zhuǎn)速可達(dá)1600 r/min(I點(diǎn)),空載時(shí)轉(zhuǎn)速可達(dá)1900 r/min(G點(diǎn))。
MTPA曲線上的交直軸電流給定值由Te與id、iq關(guān)系式?jīng)Q定:
由式(9)可見(jiàn),該關(guān)系式為二元四階高次方程組,無(wú)法得到d、q軸電流對(duì)應(yīng)的解析解。本文對(duì)兩個(gè)高次方程所對(duì)應(yīng)的平面曲線進(jìn)行擬合,在Maltlab環(huán)境下,畫(huà)出約束關(guān)系曲線,然后在曲線上取得足夠多的點(diǎn),利用cftool工具箱,得到擬合曲線如圖5所示。
圖5 MTPA擬合曲線
利用電機(jī)參數(shù),可得多項(xiàng)式擬合函數(shù)為
(10)
(11)
該擬合曲線的SSE(和方差)為0.001344,RMSE(均方差)為0.006929,R-squre(確定系數(shù))為1,表明多項(xiàng)式對(duì)高次方程的擬合具有較高的精度。
弱磁控制下的電流計(jì)算公式很復(fù)雜,為了避免巨大運(yùn)算量,用電壓差形成閉環(huán),并結(jié)合轉(zhuǎn)矩的梯度對(duì)MTPA曲線上電流值做出修正,使電機(jī)工作在弱磁區(qū)內(nèi),響應(yīng)迅速,魯棒性好。IPMSM恒轉(zhuǎn)矩方向表達(dá)示為
令
電流沿著恒轉(zhuǎn)矩梯度修正,α為修正系數(shù),修正值為
圖6 永磁同步電動(dòng)機(jī)弱磁系統(tǒng)結(jié)構(gòu)框圖
(1)額定負(fù)載起動(dòng)
圖7 額定負(fù)載運(yùn)行
(2)額定轉(zhuǎn)速運(yùn)行
圖8 額定轉(zhuǎn)速運(yùn)行
圖7~圖8表示采用id=0矢量控制下電機(jī)的轉(zhuǎn)速、dq軸電流和轉(zhuǎn)矩波形。圖7可以看出,電機(jī)帶額定負(fù)載起動(dòng),所能達(dá)到的最大速度為983 r/min,此時(shí)dq軸電流分別為0 A和5.04 A,對(duì)應(yīng)圖2中C點(diǎn) (0,5.019) ;圖10可以看出,電機(jī)轉(zhuǎn)速額定轉(zhuǎn)速下所能帶最大負(fù)載為4 Nm,此時(shí)dq軸電流分別為0 A和1.44 A,對(duì)應(yīng)圖2中A點(diǎn)(0,1.43)。理論和仿真基本吻合。
(1)額定負(fù)載起動(dòng)
圖9 額定負(fù)載運(yùn)行
(2)額定轉(zhuǎn)速運(yùn)行
圖10 額定負(fù)載運(yùn)行
圖9~圖10表示采用MTPA控制下電機(jī)轉(zhuǎn)速、dq軸電流和轉(zhuǎn)矩波形。圖9可以看出,電機(jī)帶額定負(fù)載起動(dòng),電機(jī)所能達(dá)到的最大速度為1233 r/min,此時(shí)電機(jī)的dq軸電流分別為-1.8 A和4.02 A,對(duì)應(yīng)圖3F中點(diǎn) (-1.7935,4.0184)。圖10可以看出,電機(jī)轉(zhuǎn)速為1500 r/min,所能帶最大負(fù)載為5.967 N,此時(shí)dq軸電流分別為-0.5 A和2 A,對(duì)應(yīng)圖3中D點(diǎn)(-0.510,1.996),理論和仿真基本吻合。MTPA比矢量控制更合適該種電機(jī)。
圖11 基于MTPA的弱磁控制仿真圖形
為了更好地驗(yàn)證本文所提控制方案的正確性與實(shí)用性,設(shè)計(jì)了主回路和控制回路,搭建了2.2 kW永磁同步電機(jī)實(shí)驗(yàn)平臺(tái)。系統(tǒng)采用TI公司的TMS320F2812為控制芯片,采用CAN通信實(shí)現(xiàn)下位機(jī)與上位機(jī)的數(shù)據(jù)傳輸,采樣頻率設(shè)定為4 kHz。實(shí)驗(yàn)平臺(tái)如圖12所示。
圖12 永磁同步電機(jī)實(shí)驗(yàn)平臺(tái)
直流發(fā)電機(jī)作為永磁同步電機(jī)負(fù)載,調(diào)節(jié)發(fā)電機(jī)勵(lì)磁和電樞回路電阻可以調(diào)節(jié)負(fù)載大小。
圖13 基于MTPA的弱磁控制流程圖
基于MTPA的弱磁控制流程如圖13所示,通過(guò)光電編碼器實(shí)時(shí)反饋電機(jī)轉(zhuǎn)角和轉(zhuǎn)速信息,轉(zhuǎn)角用于系統(tǒng)坐標(biāo)變換,轉(zhuǎn)速參與轉(zhuǎn)速環(huán)控制,得到輸出給定轉(zhuǎn)矩。根據(jù)電機(jī)實(shí)際參數(shù),設(shè)定曲線擬合軌跡,生成MTPA算法下電流運(yùn)行指令,電機(jī)在該指令下運(yùn)行,同時(shí)合成電壓矢量參與電壓閉環(huán)弱磁,合成矢量電壓幅值一旦達(dá)到弱磁開(kāi)啟條件,系統(tǒng)生成電流修正指令,使得電機(jī)電流軌跡向弱磁區(qū)移動(dòng)。
本文給出了矢量控制和基于MTPA弱磁控制下的實(shí)驗(yàn)結(jié)果。
圖14 矢量控制下990 r/min額定帶載運(yùn)行
圖14為矢量控制下,電機(jī)990 r/min額定帶載運(yùn)行時(shí)轉(zhuǎn)速和dq軸電流波形。在矢量控制下,帶額定負(fù)載運(yùn)行,電機(jī)起動(dòng)電流較大,2.5 s轉(zhuǎn)速達(dá)到穩(wěn)定。電機(jī)的最大轉(zhuǎn)速為990 r/min,此時(shí)dq軸電流分別為0 A和5 A,對(duì)應(yīng)圖2中的C點(diǎn)(0,5.019)。
圖15 基于MTPA弱磁算法下1600 r/min額定帶載運(yùn)行
圖16 基于MTPA弱磁算法下1900 r/min空載運(yùn)行
圖15為基于MTPA弱磁算法下電機(jī)1600 r/min額定負(fù)載運(yùn)行時(shí)轉(zhuǎn)速和dq軸電流波形。在MTPA算法控制下,電機(jī)在前1.3 s內(nèi)保持輸出轉(zhuǎn)矩限幅值,此時(shí)由MTPA算法得到id、iq分別-2 A和5.2 A,隨轉(zhuǎn)速的上升,MTPA上的點(diǎn)已經(jīng)無(wú)法滿足運(yùn)行要求,1.3 s電機(jī)達(dá)到弱磁開(kāi)啟條件,電流在梯度下降法的作用下,沿著轉(zhuǎn)矩梯度修正。為了保持高轉(zhuǎn)速,d軸去磁電流幅值變大,同時(shí)為滿足轉(zhuǎn)矩平衡關(guān)系,q軸電流有所下降,最終達(dá)到穩(wěn)定運(yùn)行狀態(tài)。電機(jī)最大轉(zhuǎn)速達(dá)到1600 r/min,比矢量控制提高了610 r/min,dq軸電流分別為-4 A和3 A,對(duì)應(yīng)圖4中Ⅰ點(diǎn)(-3.8698,3.2653)。圖16為基于MTPA弱磁算法下電機(jī)1900 r/min空載運(yùn)行時(shí)的轉(zhuǎn)速、dq軸電流波形??梢园l(fā)現(xiàn),電機(jī)穩(wěn)定運(yùn)行在1600 r/min時(shí),dq軸電流為-3 A和0.3 A,對(duì)應(yīng)圖4中G點(diǎn)(-3.0774,0.30498)。實(shí)驗(yàn)結(jié)果與理論分析和仿真較吻合,驗(yàn)證了本文所提MTPA的弱磁控制算法在內(nèi)置式永磁同步電機(jī)控制系統(tǒng)中的可行性和正確性。
本文從電機(jī)本身參數(shù)的角度出發(fā),提出了一種基于MTPA的弱磁控制算法。利用曲線擬合,解決了MTPA控制算法中d、q軸電流給定值無(wú)法獲得解析解的問(wèn)題,得到了不同轉(zhuǎn)速要求時(shí)的交直軸電流給定值,并利用轉(zhuǎn)矩的梯度和電壓閉環(huán)對(duì)MTPA指令進(jìn)行修正,擴(kuò)大了電機(jī)的運(yùn)行范圍。實(shí)驗(yàn)結(jié)果表明,在特定電機(jī)參數(shù)下,id=0的矢量控制無(wú)法使電機(jī)達(dá)到額定運(yùn)行狀態(tài),而本文所提出的基于MTPA的弱磁算法,使得電機(jī)在額定負(fù)載下,到了1600 r/min,空載條件下,電機(jī)轉(zhuǎn)速達(dá)到1900 r/min,控制精度較高,運(yùn)行性能好,提高了電機(jī)的帶載能力,拓寬了電機(jī)的調(diào)速范圍。