永磁同步電機(jī)AEKF無磁鏈直接轉(zhuǎn)矩控制

，

(福州大學(xué) 電氣工程與自動(dòng)化學(xué)院，福建 福州 350100)

1 引言

在永磁同步電動(dòng)機(jī)(PMSM)的控制方式中，直接轉(zhuǎn)矩控制(DTC)因其具有轉(zhuǎn)矩響應(yīng)速度快、控制與計(jì)算簡(jiǎn)單等特點(diǎn)被廣泛使用，傳統(tǒng)的直接轉(zhuǎn)矩控制是在定子兩相靜止坐標(biāo)系下觀測(cè)磁鏈，計(jì)算轉(zhuǎn)矩，并對(duì)兩者同時(shí)進(jìn)行滯環(huán)控制[1-3]。但考慮到控制目的僅僅是為提高電機(jī)電磁轉(zhuǎn)矩響應(yīng)的快速性，沒有理由非要將定子磁鏈幅值控制在恒定值，因此可以簡(jiǎn)化傳統(tǒng)DTC控制，引入無磁鏈DTC控制[4]，這種方法在定子磁鏈還未超過磁鏈限幅值時(shí)單獨(dú)控制電磁轉(zhuǎn)矩，由于控制方式更加直接，因此提高了電磁轉(zhuǎn)矩的響應(yīng)速度和精度，同時(shí)使轉(zhuǎn)矩脈動(dòng)有所降低。

此外，因?yàn)橹苯愚D(zhuǎn)矩控制無需旋轉(zhuǎn)坐標(biāo)變換，對(duì)轉(zhuǎn)子位置信息要求不高，所以易于實(shí)現(xiàn)無傳感器控制。目前使用較多的方法有：模型參考自適應(yīng)法[5]、滑模變結(jié)構(gòu)法[6]、擴(kuò)展卡爾曼濾波法(EKF)[7-10]等。由于EKF思路簡(jiǎn)單，應(yīng)用范圍廣泛且應(yīng)用效果尚可，因此本文采用擴(kuò)展卡爾曼濾波法，并將其引入永磁同步電機(jī)無磁鏈DTC控制中。但是，由于該方法的濾波效果受噪聲影響很大，錯(cuò)誤的噪聲協(xié)方差矩陣會(huì)使系統(tǒng)發(fā)散，而該矩陣的選取目前往往依靠試湊得到，所以本文在擴(kuò)展卡爾曼濾波的基礎(chǔ)上引入噪聲自適應(yīng)算法，目前較為常用的自適應(yīng)擴(kuò)展卡爾曼濾波器有：相關(guān)法自適應(yīng)濾波[11]、基于極大似然準(zhǔn)則的自適應(yīng)濾波[12]、Sage-Husa自適應(yīng)濾波[13-14]、協(xié)方差匹配自適應(yīng)濾波[15]等。本文將協(xié)方差匹配自適應(yīng)濾波與Sage-Husa自適應(yīng)濾波相結(jié)合，使系統(tǒng)噪聲和量測(cè)噪聲可以自適應(yīng)調(diào)節(jié)，提高系統(tǒng)的收斂性能以及觀測(cè)精度。

2 永磁同步電機(jī)數(shù)學(xué)模型

本文所有研究都是針對(duì)正弦波隱極式永磁同步電動(dòng)機(jī)，且電機(jī)定子繞組為星形連接，這類永磁同步電機(jī)應(yīng)用最為廣泛，隱極式永磁同步電機(jī)交直軸電感相同，無磁阻轉(zhuǎn)矩分量，計(jì)算較為方便。

xy坐標(biāo)系下永磁同步電動(dòng)機(jī)數(shù)學(xué)模型表示為：

(1)

(2)

由于定子磁鏈定向于x軸，所以φy=0，于是由(1)、(2)可得轉(zhuǎn)矩方程為：

(3)

φx、φy、ux、uy、ix、iy分別為定子磁鏈、電壓和電流x、y軸分量；φf為永磁體磁鏈；|φs|為定子磁鏈幅值；δ為定轉(zhuǎn)子磁鏈夾角；Rs、Ls分別表示定子電阻和電感；p為電機(jī)極對(duì)數(shù)；ωr為轉(zhuǎn)子磁鏈電角速度。

3 無磁鏈DTC原理概述

傳統(tǒng)PMSM DTC理論是根據(jù)式(3)，保持定子磁鏈幅值恒定，控制轉(zhuǎn)矩角δ的變化實(shí)現(xiàn)對(duì)電磁轉(zhuǎn)矩的實(shí)時(shí)控制。無磁鏈DTC不再保持定子磁鏈幅值恒定，只把轉(zhuǎn)矩當(dāng)做主要控制量，實(shí)現(xiàn)更快速更準(zhǔn)確的轉(zhuǎn)矩響應(yīng)。

式(3)兩邊對(duì)δ求導(dǎo)，可得：

(4)

當(dāng)開關(guān)頻率很高時(shí)，轉(zhuǎn)矩和轉(zhuǎn)矩角變化很小，這時(shí)(4)可以化簡(jiǎn)為：

(5)

圖1 定子磁鏈與電壓矢量關(guān)系圖

結(jié)合圖1，因?yàn)樵谝粋€(gè)開關(guān)周期內(nèi)磁鏈變化很小，可以認(rèn)為DA與OD垂直，因此得到轉(zhuǎn)矩角變化量：

(6)

將上式帶入(5)，可進(jìn)一步化簡(jiǎn)為：

(7)

如圖1所示，做DE垂直于AB，由于DA垂直于OD，則∠BAD=∠ODE=δ，AE可以視為Δφs在交軸上的分量，當(dāng)忽略定子電阻時(shí)，在一個(gè)開關(guān)周期內(nèi)定子磁鏈交軸分量與電壓矢量交軸分量的關(guān)系可表示為：

φq=uqT

(8)

結(jié)合式(5)電磁轉(zhuǎn)矩可表示為：

(9)

根據(jù)式(9)可以看出，為了更加突出直接轉(zhuǎn)矩直接控制轉(zhuǎn)矩這一特點(diǎn)，可以只控制電壓矢量的交軸分量，不需要考慮磁鏈，該方法控制簡(jiǎn)單快捷，理論上可以提升轉(zhuǎn)矩的動(dòng)態(tài)響應(yīng)性能。

圖2 空間電壓矢量與各坐標(biāo)系關(guān)系圖

由圖2可以得到六個(gè)非零電壓矢量交軸分量表達(dá)式為：

(10)

如圖2所示，在2π周期內(nèi)可以分為6個(gè)扇區(qū)，每個(gè)扇區(qū)選擇對(duì)電磁轉(zhuǎn)矩影響最大的一個(gè)電壓矢量，拿一扇區(qū)(0°～60°)舉例來說，由式10可知，當(dāng)需要增大轉(zhuǎn)矩時(shí)，u3q對(duì)電磁轉(zhuǎn)矩影響最大；當(dāng)需要減小轉(zhuǎn)矩時(shí)，采用u6q可以最快的降低電磁轉(zhuǎn)矩；為了減小轉(zhuǎn)矩脈動(dòng)，在一定范圍內(nèi)將零矢量作用于電機(jī)。各扇區(qū)電壓矢量的選擇如下表所示，當(dāng)τ為1時(shí)表示需要增加電磁轉(zhuǎn)矩，當(dāng)τ為0時(shí)，表示采用零矢量保持電磁轉(zhuǎn)矩，當(dāng)τ為-1時(shí)表示需要減小電磁轉(zhuǎn)矩。

表1 無磁鏈條件下電壓矢量選擇表

上述控制方法只考慮了定子磁鏈的交軸分量，但電壓矢量作用后還會(huì)使磁鏈直軸分量產(chǎn)生變化，一定情況下定子磁鏈幅值可能會(huì)超過額定值，導(dǎo)致定子繞組磁飽和，嚴(yán)重時(shí)可使電機(jī)停轉(zhuǎn)。因此需要在磁鏈超出額定值時(shí)對(duì)其限幅。XY坐標(biāo)系以定子磁鏈定向，因此在忽略定子電阻的情況下，空間電壓矢量的X軸分量可以直接決定定子磁鏈幅值，此外Y軸分量決定了電機(jī)轉(zhuǎn)矩角的變化，可以控制電磁轉(zhuǎn)矩的增減。

由圖2可以得到非零電壓矢量的XY軸分量分別為：

(11)

(12)

同樣可將2π周期分為6個(gè)扇區(qū)，拿一扇區(qū)(-30°～30°)舉例，根據(jù)式(11)、(12)此區(qū)間內(nèi)可以減小磁鏈幅值的電壓矢量有u3、u4、u5三個(gè)矢量，如果此時(shí)需要增加電磁轉(zhuǎn)矩只能選擇u3；反之需減小電磁轉(zhuǎn)矩時(shí)只能選擇u5；假如需要保持電磁轉(zhuǎn)矩需要選擇在此扇區(qū)中正負(fù)變化量相等的電壓矢量，即u4。各扇區(qū)電壓矢量選擇如表2所示。

表2 磁鏈限幅條件下電壓矢量選擇表

綜上，在磁鏈未達(dá)到額定磁鏈時(shí)參照表1進(jìn)行控制，當(dāng)磁鏈超過額定限幅，參照表2進(jìn)行控制，以上就構(gòu)成了無磁鏈直接轉(zhuǎn)矩控制的理論基礎(chǔ)。

4 AEKF觀測(cè)器設(shè)計(jì)

4.1 EKF原理概述

EKF濾波步驟可分為以下兩步：

(1)預(yù)測(cè)階段。計(jì)算狀態(tài)預(yù)測(cè)值x和預(yù)測(cè)誤差方差陣P，即：

(13)

(2)更新階段。計(jì)算EKF增益矩陣K，進(jìn)行狀態(tài)預(yù)測(cè)值和預(yù)測(cè)誤差方差陣的更新，即：

增益矩陣：

(14)

狀態(tài)預(yù)測(cè)值和預(yù)測(cè)誤差方差陣更新：

(15)

Φk,k-1≈

(16)

(17)

(18)

將以上矩陣代入式(13)～(15)經(jīng)過數(shù)次迭代后可以估計(jì)出狀態(tài)變量真實(shí)值，用于PMSM DTC中。

4.1 系統(tǒng)噪聲估計(jì)器

根據(jù)理論分析以及實(shí)驗(yàn)探究，擴(kuò)展卡爾曼濾波必須在狀態(tài)模型和噪聲統(tǒng)計(jì)特性準(zhǔn)確描述的情況下才可以達(dá)到預(yù)期效果。系統(tǒng)噪聲協(xié)方差矩陣與量測(cè)噪聲協(xié)方差矩陣的選取直接影響整個(gè)系統(tǒng)的收斂效果，雖然在電機(jī)模型下噪聲矩陣是簡(jiǎn)單白噪聲，但是目前選取噪聲協(xié)方差矩陣主要靠試湊得到。因此本文提出一種自適應(yīng)卡爾曼濾波理論，可以準(zhǔn)確的得到隨系統(tǒng)參數(shù)變化的噪聲協(xié)方差矩陣，保證系統(tǒng)快速收斂。

本文將Sage-Husa自適應(yīng)擴(kuò)展卡爾曼濾波與協(xié)方差匹配法相結(jié)合，Sage-Husa自適應(yīng)擴(kuò)展卡爾曼濾波用來估計(jì)系統(tǒng)噪聲，協(xié)方差匹配法用來估計(jì)量測(cè)噪聲。該方法既避免了Sage-Husa自適應(yīng)濾波法計(jì)算復(fù)雜的缺點(diǎn)，另一方面已有文獻(xiàn)證明，系統(tǒng)噪聲與量測(cè)噪聲不能被Sage-Husa自適應(yīng)濾波法同時(shí)觀測(cè)，所以該方法保證兩類噪聲都被觀測(cè)的前提下仍然保證了系統(tǒng)魯棒性。

Sage-Husa自適應(yīng)擴(kuò)展卡爾曼濾波時(shí)變系統(tǒng)噪聲估計(jì)器表達(dá)式如下：

(19)

(20)

根據(jù)卡爾曼迭代公式可知濾波新息為：

(21)

式(19)中b表示遺忘因子，其取值范圍為0.95～0.995，本文實(shí)驗(yàn)中取0.95。由于最后一項(xiàng)數(shù)值很小，為方便收斂，將(20)化簡(jiǎn)為系統(tǒng)噪聲協(xié)方差矩陣Qk的有偏估計(jì)：

(22)

此時(shí)將Qk參與EKF迭代過程，實(shí)現(xiàn)自適應(yīng)濾波。

4.2 量測(cè)噪聲估計(jì)器

將Yk=Hkxk+vk代入(21)，得：

(23)

由(23)式可得到新息的理論協(xié)方差為：

(24)

實(shí)際的新息協(xié)方差可表示為：

(25)

其中N為滑動(dòng)窗口大小。

將協(xié)方差實(shí)際值等效為理論值，可以得到量測(cè)噪聲協(xié)方差矩陣為：

(26)

該方法的觀測(cè)性能受滑動(dòng)窗口大小的影響，較小的滑動(dòng)窗口計(jì)算量小，但觀測(cè)結(jié)果波動(dòng)性大，較大的滑動(dòng)窗口計(jì)算量過大對(duì)硬件性能要求苛刻，本文實(shí)驗(yàn)中選取的滑動(dòng)窗口為100。因此，在前100個(gè)控制周期內(nèi)，可以先固定量測(cè)噪聲協(xié)方差矩陣，利用Sage-Husa自適應(yīng)擴(kuò)展卡爾曼濾波得到系統(tǒng)噪聲協(xié)方差矩陣的值，隨后進(jìn)行量測(cè)噪聲的估計(jì)，使得在每個(gè)控制周期內(nèi)都有至少一個(gè)噪聲協(xié)方差矩陣在修正，提高EKF觀測(cè)性能。

5 仿真與實(shí)驗(yàn)驗(yàn)證

5.1 仿真分析

為了驗(yàn)證以上理論的正確性，用Matlab/Simulink進(jìn)行仿真分析，其中AEKF算法用S-function來實(shí)現(xiàn)，電機(jī)主要參數(shù)如表1所示。

表1 永磁同步電機(jī)主要參數(shù)

仿真中轉(zhuǎn)速給定為600r/min，負(fù)載轉(zhuǎn)矩為空載啟動(dòng)，0.2s后突加負(fù)載至4N*m，分別采用傳統(tǒng)零矢量DTC、無磁鏈DTC進(jìn)行仿真分析。兩者的磁鏈、轉(zhuǎn)矩與轉(zhuǎn)速波形如圖3～6所示。

圖3 傳統(tǒng)零矢量DTC定子磁鏈波形

圖4 傳統(tǒng)零矢量DTC轉(zhuǎn)矩和轉(zhuǎn)速波形

圖5 無磁鏈DTC定子磁鏈波形

圖6 無磁鏈DTC轉(zhuǎn)矩和轉(zhuǎn)速波形

根據(jù)仿真波形可以看出，無磁鏈DTC定子磁鏈波形并非標(biāo)準(zhǔn)的圓形，在空載時(shí)定子磁鏈大部分時(shí)間小于磁鏈限幅，增大轉(zhuǎn)矩至4N*m時(shí)，磁鏈近似于圓形磁鏈但始終沒有超出限幅，此外由于在無磁鏈DTC更加直接的控制電磁轉(zhuǎn)矩，其轉(zhuǎn)矩脈動(dòng)更小，響應(yīng)速度上無磁鏈DTC更快速。

在無磁鏈DTC控制的基礎(chǔ)上加入AEKF算法，在試驗(yàn)了較多的噪聲協(xié)方差矩陣后，發(fā)現(xiàn)AEKF在許多EKF無法收斂的噪聲矩陣下都能持續(xù)收斂，且收斂速度明顯優(yōu)于EKF算法，參數(shù)矩陣選擇如下；將AEKF觀測(cè)器參與閉環(huán)，磁鏈、轉(zhuǎn)矩、轉(zhuǎn)速以及轉(zhuǎn)子位置角仿真波形如圖7～9所示。

圖7 AEKF無磁鏈DTC定子磁鏈波形

圖8 AEKF無磁鏈DTC轉(zhuǎn)矩和轉(zhuǎn)速波形

圖9 實(shí)際位置角與AEKF估計(jì)位置角對(duì)比波形

由仿真波形可以看出，AEKF無磁鏈DTC控制不僅保留了其轉(zhuǎn)矩響應(yīng)速度快的特點(diǎn)，而且磁鏈、轉(zhuǎn)矩脈動(dòng)有所減小?？梢钥吹剑珹EKF觀測(cè)器下轉(zhuǎn)子位置角波形略滯后于實(shí)際位置角，但誤差很小，因此仍能保證電機(jī)穩(wěn)定運(yùn)行。

5.2 實(shí)驗(yàn)分析

本文采用TI 公司DSP TMS320F2812作為控制核心，功率變換裝置采用富士公司7MBP50RA120智能功率模塊，永磁同步電機(jī)為隱極式正弦波伺服電機(jī)，電機(jī)主要參數(shù)與表1相同，所帶負(fù)載為一個(gè)三相同步測(cè)功機(jī)。分別對(duì)無磁鏈DTC、EKF無磁鏈DTC及AEKF無磁鏈DTC進(jìn)行實(shí)驗(yàn)分析，給定轉(zhuǎn)速為600r/min，空載啟動(dòng)待轉(zhuǎn)速穩(wěn)定后，突加2.7N*m負(fù)載轉(zhuǎn)矩。圖10～17為CCS3.3程序示波器觀測(cè)下的無磁鏈DTC、EKF無磁鏈DTC和AEKF無磁鏈DTC轉(zhuǎn)速轉(zhuǎn)矩、定子磁鏈以及位置角波形。

圖10 無磁鏈DTC轉(zhuǎn)速和轉(zhuǎn)矩波形

圖11 無磁鏈DTC定子磁鏈波形

圖12 EKF無磁鏈DTC轉(zhuǎn)速和轉(zhuǎn)矩波形

圖13 EKF無磁鏈DTC定子磁鏈波形

圖14 實(shí)際位置角與EKF估計(jì)位置角對(duì)比波形

圖15 AEKF無磁鏈DTC轉(zhuǎn)速和轉(zhuǎn)矩波形

根據(jù)圖12和圖15可以看出，EKF和AEKF無磁鏈DTC對(duì)轉(zhuǎn)速轉(zhuǎn)矩脈動(dòng)都有了明顯改善，轉(zhuǎn)矩脈動(dòng)較無磁鏈DTC減小了約0.65N*m，但EKF無磁鏈DTC在突加轉(zhuǎn)矩后收斂速度較慢，突加時(shí)刻會(huì)產(chǎn)生較大超調(diào)，而AEKF無磁鏈DTC有效的克服了這一點(diǎn)，保證了較快的轉(zhuǎn)矩響應(yīng)及收斂性，轉(zhuǎn)矩超調(diào)也明顯減?。粚?duì)比圖11、13、16可以看出加入AEKF觀測(cè)器的定子磁鏈脈動(dòng)更??；位置角波形方面EKF和AEKF觀測(cè)器都有準(zhǔn)確的觀測(cè)效果。

圖16 EKF無磁鏈DTC定子磁鏈波形

圖17 實(shí)際位置角與AEKF估計(jì)位置角對(duì)比波形

6 結(jié)論

為了優(yōu)化傳統(tǒng)DTC的控制性能，本文提出了星形連接正弦波永磁同步電動(dòng)機(jī)無磁鏈DTC控制，該方法在定子磁鏈到達(dá)限幅之前只控制電磁轉(zhuǎn)矩大小，由于控制直接作用于轉(zhuǎn)矩，因此該方法加快了轉(zhuǎn)矩響應(yīng)速度并減小了轉(zhuǎn)矩脈動(dòng)。為了進(jìn)一步優(yōu)化控制效果，在無磁鏈DTC的基礎(chǔ)上引入EKF無傳感器控制，該方法不僅提高了實(shí)驗(yàn)研究的經(jīng)濟(jì)性而且估計(jì)數(shù)據(jù)接近于真實(shí)值，使得轉(zhuǎn)矩脈動(dòng)有明顯下降，但EKF受噪聲矩陣影響較大，而噪聲矩陣目前又只能通過試湊的方法得到，這大大限制了EKF算法的收斂性，導(dǎo)致在電機(jī)負(fù)載變化時(shí)系統(tǒng)不能快速收斂。所以本文提出一種新型自適應(yīng)EKF算法，將Sage-Husa法與協(xié)方差匹配法相結(jié)合，并通過仿真和實(shí)驗(yàn)證明AEKF算法有更好的觀測(cè)效果和更快的收斂速度。

[1] Rahman M F,Zhong L,Haque M DEnamul,et al.A direct torque-controlled interior permanent-magnet synchronous motor drive without a speed sensor[J].IEEE Transactions on Energy Conversion,2003,18(1):17-22.

[2] HU Yuwen,TIAN Chun,YOU Zhiqing,et al.Direct torque control system and sensorless technique of permanent magnet synchronous motor [J].Chinese Journal of Aeronautics,2003,16(2):97-102.

[3] 陳振,劉向東,戴亞平,等.采用預(yù)期電壓矢量調(diào)制的PMSM直接轉(zhuǎn)矩控制[J].電機(jī)與控制學(xué)報(bào),2009,13(1):43-44.

[4] 胡育文,高瑾,楊建飛,等.永磁同步電動(dòng)機(jī)直接轉(zhuǎn)矩控制系統(tǒng)[M].機(jī)械工業(yè)出版社,2015:54-55+60.

[5] 佘致廷,肖岸文,孫煒,等.基于MRAS理論的無速度傳感器直接轉(zhuǎn)矩控制調(diào)速系統(tǒng)[J].電工技術(shù)學(xué)報(bào),2006,21(4):99-100.

[6] 孫振興,張興華.基于滑模觀測(cè)器的感應(yīng)電機(jī)無速度傳感器直接轉(zhuǎn)矩控制[J].電工電能新技術(shù),2012,31(4):30-31.

[7] 劉英培,萬健如,沈虹,等.基于EKF PMSM定子磁鏈和轉(zhuǎn)速觀測(cè)直接轉(zhuǎn)矩控制[J].電工技術(shù)學(xué)報(bào),2009,24(12):58-59.

[8] 劉英培,萬健如,梁鵬飛.基于擴(kuò)展卡爾曼濾波器和空間電壓矢量調(diào)制的永磁同步電機(jī)直接轉(zhuǎn)矩控制[J].中國(guó)電機(jī)工程學(xué)報(bào),2009,29(27):70-72.

[9] 陳振,劉向東,靳永強(qiáng),等.采用擴(kuò)展卡爾曼濾波磁鏈觀測(cè)器的永磁同步電機(jī)直接轉(zhuǎn)矩控制[J].中國(guó)電機(jī)工程學(xué)報(bào),2008,28(33):77-81.

[10] 刑巖,王旭,楊丹,等.改進(jìn)定子磁鏈估計(jì)器及其參數(shù)選取[J].電機(jī)與控制學(xué)報(bào),2016,20(3):32-33.

[11] 計(jì)會(huì)鳳.基于浮動(dòng)車GPS數(shù)據(jù)的動(dòng)態(tài)交通預(yù)測(cè)與誘導(dǎo)模型研究[D].遼寧:遼寧工程技術(shù)大學(xué),2009:71-73.

[12] 岳曉奎,袁建平.一種基于極大似然準(zhǔn)則的自適應(yīng)卡爾曼濾波算法[J].西北工業(yè)大學(xué)學(xué)報(bào),2005,23(4):469-472.

[13] 閔偉,周志宇,周振華.改進(jìn)型Sage-Husa卡爾曼濾波器在電壓暫降檢測(cè)中的應(yīng)用[J].電網(wǎng)技術(shù),2013,37(1):231-232.

[14] 楊學(xué)海,張繼業(yè),張晗.基于改進(jìn)Sage-Husa的自適應(yīng)無跡卡爾曼濾波的鋰離子電池SOC估計(jì)[J].電工電能新技術(shù),2016,35(1):32.

[15] 趙琳,等.非線性系統(tǒng)濾波理論[M].國(guó)防工業(yè)出版社,2012:36-41+193.