王飛宇, 田井呈, 卓克瓊, 趙朝會(huì)
(上海電機(jī)學(xué)院 電氣學(xué)院,上?!?01306)
?
基于改進(jìn)模型參考自適應(yīng)算法的永磁同步電機(jī)轉(zhuǎn)動(dòng)慣量辨識(shí)
王飛宇,田井呈,卓克瓊,趙朝會(huì)
(上海電機(jī)學(xué)院 電氣學(xué)院,上海201306)
針對使用固定增益系數(shù)的模型參考自適應(yīng)(MRAS)算法在穩(wěn)定性與收斂響應(yīng)速度不可同時(shí)兼得的問題,介紹了基于MRAS算法的永磁同步電機(jī)轉(zhuǎn)動(dòng)慣量辨識(shí)的原理,研究了該算法中自適應(yīng)增益系數(shù)對電機(jī)轉(zhuǎn)動(dòng)慣量辨識(shí)結(jié)果的影響。探討了一種依據(jù)辨識(shí)參數(shù)偏差而實(shí)時(shí)改變的變增益系數(shù)的方法,并通過MATLAB/Simulink對參數(shù)辨識(shí)系統(tǒng)進(jìn)行了仿真驗(yàn)證。仿真結(jié)果表明,該方法收斂速度快,參數(shù)辨識(shí)結(jié)果振蕩幅值小,穩(wěn)定性好。
永磁同步電機(jī); 模型參考自適應(yīng); 轉(zhuǎn)動(dòng)慣量辨識(shí)
交流永磁同步電機(jī)伺服控制系統(tǒng)被廣泛應(yīng)用于高精度、高性能的場合。然而,伺服系統(tǒng)的動(dòng)態(tài)性能會(huì)受電機(jī)內(nèi)部參數(shù)、外部擾動(dòng)轉(zhuǎn)矩和一些非線性因素的影響[1]。由于轉(zhuǎn)動(dòng)慣量對伺服系統(tǒng)的控制性能影響較大,是控制系統(tǒng)主要的性能指標(biāo)[2],因此辨識(shí)出整個(gè)伺服系統(tǒng)的轉(zhuǎn)動(dòng)慣量變得非常重要。
國內(nèi)外研究人員對轉(zhuǎn)動(dòng)慣量辨識(shí)的方法做了大量的研究?,F(xiàn)有的轉(zhuǎn)動(dòng)慣量辨識(shí)方案可歸納為兩類: 離線轉(zhuǎn)動(dòng)慣量辨識(shí)方案和在線轉(zhuǎn)動(dòng)慣量辨識(shí)方案。離線轉(zhuǎn)動(dòng)慣量辨識(shí)方案是在系統(tǒng)控制程序運(yùn)行前通過直接測試法、加減速法和人工軌跡法等[3- 6]對慣量進(jìn)行辨識(shí),將測得的參數(shù)提供給控制程序使用。離線式辨識(shí)方法需要對大量的數(shù)據(jù)進(jìn)行存儲(chǔ)和運(yùn)算,不具有實(shí)時(shí)性且辨識(shí)精度差,這些缺點(diǎn)限制了其應(yīng)用于高性能的控制系統(tǒng)的范圍[7]。隨著大量的現(xiàn)代控制理論應(yīng)用于電機(jī)控制領(lǐng)域,逐漸形成一套在線式辨識(shí)理論,可對運(yùn)行中的系統(tǒng)的轉(zhuǎn)動(dòng)慣量等參數(shù)進(jìn)行實(shí)時(shí)辨識(shí),并將所辨識(shí)的參數(shù)實(shí)時(shí)提供給系統(tǒng)的控制程序使用,實(shí)時(shí)性較強(qiáng),控制精度高。常用的在線式轉(zhuǎn)動(dòng)慣量辨識(shí)方法有: 最小二乘法、卡爾曼濾波器、觀測器法和模型參考自適應(yīng)法(Model Reference Adaptive System, MRAS)等。文獻(xiàn)[8]使用遞推最小二乘法對直流電機(jī)的轉(zhuǎn)動(dòng)慣量進(jìn)行辨識(shí),但算法需要加高階濾波器,適用于控制簡單的系統(tǒng)。文獻(xiàn)[9]對具有遺忘因子的最小二乘法進(jìn)行了改進(jìn),克服了辨識(shí)器的波動(dòng)問題。文獻(xiàn)[10]提出了一種基于擴(kuò)展卡爾曼濾波的辨識(shí)方法,可對系統(tǒng)的轉(zhuǎn)動(dòng)慣量和負(fù)載轉(zhuǎn)矩等實(shí)時(shí)辨識(shí),但需要進(jìn)行復(fù)雜的計(jì)算和非線性濾波的處理。使用狀態(tài)觀測器法[11]對轉(zhuǎn)動(dòng)慣量進(jìn)行辨識(shí)需要通過對模型的非線性、多變量和多狀態(tài)量等進(jìn)行處理,存在多種誤差[12]。MRAS法因具有算法簡單、便于使用在數(shù)字控制系統(tǒng)等優(yōu)點(diǎn),研究人員對其做了大量研究。文獻(xiàn)[13]結(jié)合MRAS理論和閉環(huán)磁鏈觀測器構(gòu)造出異步電機(jī)的速度自適應(yīng)觀測器,成功地對電機(jī)的轉(zhuǎn)速進(jìn)行了辨識(shí)。文獻(xiàn)[14]將最小二乘法加入了MRAS辨識(shí)系統(tǒng),提出了一種改進(jìn)的MRAS參數(shù)辨識(shí)方法,改善了辨識(shí)速度和振蕩幅度。文獻(xiàn)[15]將變結(jié)構(gòu)控制與模型自適應(yīng)系統(tǒng)結(jié)合,提出了一種變結(jié)構(gòu)MRAS方法,對參數(shù)辨識(shí)的動(dòng)靜態(tài)性能較好。文獻(xiàn)[16]對MRAS做了改進(jìn),通過設(shè)計(jì)兩個(gè)自適應(yīng)增益系數(shù)大小不同的系統(tǒng),在一定程度上改善了辨識(shí)模型對增益參數(shù)的依賴。
本文通過分析MRAS系統(tǒng)中自適應(yīng)增益系數(shù)對永磁同步電機(jī)轉(zhuǎn)動(dòng)慣量辨識(shí)的收斂速度與參數(shù)穩(wěn)定性的影響,給出了一種依據(jù)辨識(shí)誤差提供增益系數(shù)的控制方案,通過仿真得出該方案具有抑制辨識(shí)結(jié)果波動(dòng)和收斂速度快等優(yōu)勢。
1.1電機(jī)數(shù)學(xué)模型
永磁同步電機(jī)的數(shù)學(xué)模型建立在不計(jì)渦流損耗、磁滯損耗和忽略鐵心飽和等條件下,在id=0的矢量控制系統(tǒng)下對永磁同步電機(jī)轉(zhuǎn)動(dòng)慣量辨識(shí),其中運(yùn)功方程(1)與轉(zhuǎn)矩方程(2)是對轉(zhuǎn)動(dòng)慣量辨識(shí)的基本方程。
(1)
(2)
式中:J——轉(zhuǎn)動(dòng)慣量;
ωr——轉(zhuǎn)子角速度;
Te——電磁轉(zhuǎn)矩;
TL——負(fù)載轉(zhuǎn)矩;
p——電機(jī)極對數(shù);
ψf——永磁鐵的勵(lì)磁磁鏈;
iq——兩相旋轉(zhuǎn)q軸電流分量。
1.2MRAS轉(zhuǎn)動(dòng)慣量辨識(shí)
MRAS方法由參考模型和可調(diào)模型組成,且兩個(gè)模型具有相同物理意義的輸入量和輸出量,依據(jù)兩個(gè)模型輸出量的誤差構(gòu)成合適的自適應(yīng)辨識(shí)算法。當(dāng)通過算法調(diào)節(jié)兩個(gè)模型的輸出量無限逼近時(shí),便可辨識(shí)出實(shí)際物理模型中的待辨識(shí)參數(shù)。永磁同步電機(jī)轉(zhuǎn)動(dòng)慣量辨識(shí)方法中參考模型與可調(diào)模型的確立過程如下。
將式(1)離散和簡化可得
ωm(k)=2ωm(k-1)-ωm(k-2)+
(3)
式中:T——系統(tǒng)的采樣周期。
在高性能的永磁同步電機(jī)控制系統(tǒng)中,在一個(gè)采樣周期內(nèi)負(fù)載轉(zhuǎn)矩的變化可近似保持不變,則有
ωm(k)=2ωm(k-1)-ωm(k-2)+
bΔTe(k-1)
(4)
其中:b=T/J;
ΔTe(k-1)——一個(gè)采樣周期的電磁轉(zhuǎn)矩的差值。
將式(4)作為參考模型,則可得出可調(diào)模型的方程:
ωg(k)=2ωm(k-1)-ωm(k-2)+
bg(k-1)ΔTe(k-1)
(5)
式中:ωg——估計(jì)的速度信號;
ωm——速度信號。
由文獻(xiàn)[17]可得出永磁同步電機(jī)轉(zhuǎn)動(dòng)慣量MRAS規(guī)律:
(6)
式中:β——自適應(yīng)增益系數(shù);
Δω(k)——速度信號與估計(jì)速度信號的差值。
依據(jù)MRAS算法的原理以及以上推導(dǎo)過程,可得到轉(zhuǎn)動(dòng)慣量辨識(shí)的原理圖,如圖1所示。
圖1 轉(zhuǎn)動(dòng)慣量辨識(shí)原理圖
在MRAS算法中,自適應(yīng)增益系數(shù)β的取值大小對轉(zhuǎn)動(dòng)慣量辨識(shí)的結(jié)果影響較大。一般在辨識(shí)系統(tǒng)中,自適應(yīng)增益系數(shù)一旦選定,在辨識(shí)的過程中就不會(huì)變化,這就給參數(shù)辨識(shí)帶來了問題: 自適應(yīng)增益系數(shù)越大,參數(shù)辨識(shí)的收斂速度越快,但辨識(shí)結(jié)果的波動(dòng)性越大;自適應(yīng)系數(shù)越小,參數(shù)辨識(shí)波動(dòng)性越小,但參數(shù)辨識(shí)結(jié)果的收斂速度越慢。實(shí)際過程中綜合考慮收斂速度與波動(dòng)性,選取一個(gè)合適的折中值,但這樣會(huì)犧牲掉大增益系數(shù)時(shí)的收斂速度快和小增益系數(shù)時(shí)波動(dòng)性小的優(yōu)勢。
為了充分利用自適應(yīng)增益系數(shù)在最大值和最小值時(shí)的優(yōu)勢,設(shè)計(jì)了一個(gè)依據(jù)辨識(shí)參數(shù)偏差的變化而變化的可變增益系數(shù)的方法,使在輸入偏差較大的階段選擇較大的自適應(yīng)增益系數(shù),在輸入偏差較小的階段選擇較小的自適應(yīng)增益系數(shù),獲得收斂速度快和辨識(shí)結(jié)果波動(dòng)小的參數(shù)辨識(shí)結(jié)果。變增益系數(shù)β值的推導(dǎo)公式為
β=βmin+[e(J)]2βmax
(7)
式中:βmin——穩(wěn)定性最好的自適應(yīng)增益系數(shù);
e(J)——轉(zhuǎn)動(dòng)慣量實(shí)際值與辨識(shí)值偏差的標(biāo)幺化取值,范圍[0,1];
βmax——最大收斂自適應(yīng)增益系數(shù)。
可變增益系數(shù)β值的推導(dǎo)流程圖與變增益系數(shù)參數(shù)辨識(shí)系統(tǒng)原理圖如圖2所示。
圖2 變增益系數(shù)法流程圖和變增益系數(shù)參數(shù)辨識(shí)系統(tǒng)原理圖
通過MATLAB/Simulink搭建仿真模型,采用id=0的矢量控制。待辨識(shí)的電機(jī)參數(shù)如下: 極對數(shù)p=4;永磁體磁鏈0.175Wb;電機(jī)相繞組電阻2.875Ω;直軸、交軸電感均為0.85mH;轉(zhuǎn)動(dòng)慣量J=0.00085kg·m2。轉(zhuǎn)動(dòng)慣量辨識(shí)系統(tǒng)仿真圖如圖3所示。
依據(jù)參考模型式(4)、可調(diào)模型式(5)和自適應(yīng)率算法式(6),在MATLAB/Simulink中搭建轉(zhuǎn)動(dòng)慣量辨識(shí)子模塊如圖4所示。
3.1自適應(yīng)增益系數(shù)對轉(zhuǎn)動(dòng)慣量辨識(shí)的影響
在式(6)中,自適應(yīng)增益系數(shù)β的取值會(huì)影響轉(zhuǎn)動(dòng)慣量辨識(shí)的結(jié)果。取自適應(yīng)增益系數(shù)為β=30,β=20,β=10,β=5對永磁同步電機(jī)的轉(zhuǎn)動(dòng)慣量辨識(shí),仿真結(jié)果如圖5所示。
由圖5可知,自適應(yīng)增益系數(shù)β=30時(shí),收斂時(shí)間為0.05s,收斂速度最快;當(dāng)β=20時(shí),收斂時(shí)間為0.07s,收斂速度較快;當(dāng)β=10時(shí),收斂時(shí)間為0.14s,收斂速度稍慢;當(dāng)β=5時(shí),收斂時(shí)間為0.25s,收斂速度最慢。所以,自適應(yīng)增益系數(shù)越大,轉(zhuǎn)動(dòng)慣量的辨識(shí)波形收斂速度越快。另外,隨著自適應(yīng)增益系數(shù)逐漸增大,收斂速度變化的趨勢降低,當(dāng)達(dá)到一定收斂速度以后,增大自適應(yīng)系數(shù)對收斂速度的改變作用變小。
圖3 永磁同步電機(jī)轉(zhuǎn)動(dòng)慣量辨識(shí)系統(tǒng)仿真框圖
圖4 MRAS轉(zhuǎn)動(dòng)慣量辨識(shí)算法仿真圖
圖5 不同增益系數(shù)時(shí)轉(zhuǎn)動(dòng)慣量辨識(shí)波形圖
在t[0.3-0.38],J[0.0007-0.00095]區(qū)間取值,可得到在不同增益系數(shù)下的局部放大圖,如圖6所示。
由圖6(a)可知,β=30時(shí)波形振蕩幅度最大;由圖6(d)可知,β=5時(shí)波形振蕩幅值最小。綜合分析得出: 自適應(yīng)增益系數(shù)越大,轉(zhuǎn)動(dòng)慣量辨識(shí)的結(jié)果振蕩幅度越大,波動(dòng)性越強(qiáng);自適應(yīng)增益系數(shù)越小,辨識(shí)結(jié)果振蕩幅度越小,辨識(shí)精度越高。
3.2變自適應(yīng)增益系數(shù)轉(zhuǎn)動(dòng)慣量辨識(shí)的仿真
在MATLAB中按圖2流程圖對自適應(yīng)增益系數(shù)β進(jìn)行仿真,采取s函數(shù)編寫變增益系數(shù),取βmax=30、βmin=5,初始轉(zhuǎn)動(dòng)慣量取值為J=0.00085kg·m2。仿真采用變增益系數(shù)法和固定增益值為(β=5、β=30)作對比分析,轉(zhuǎn)動(dòng)慣量辨識(shí)的仿真波形如圖7所示。
為了對比分析采用變增益系數(shù)法和β=5、β=30對參數(shù)辨識(shí)結(jié)果的振蕩幅值和穩(wěn)定性的影響,在t[0.3-0.33],J[0.0007-0.00095]區(qū)間取值,可得出局部放大圖如圖7(b)所示。
圖6 不同增益系數(shù)時(shí)辨識(shí)結(jié)果局部放大圖
圖7 變增益系數(shù)法和β=5、β=30的轉(zhuǎn)動(dòng)慣量辨識(shí)
由圖7(a)對比波形可知: 變增益系數(shù)法與固定增益β=30時(shí)轉(zhuǎn)動(dòng)慣量辨識(shí)結(jié)果的收斂時(shí)間約為0.04s,收斂速度較快;當(dāng)β=5時(shí),收斂時(shí)間為0.25s,收斂時(shí)間長,收斂速度較慢。由圖7(b) 振蕩幅度比較波形可知,變增益系數(shù)法與固定增益β=5時(shí)辨識(shí)結(jié)果的振蕩幅度較小,辨識(shí)誤差為2.3%;當(dāng)β=30時(shí),振蕩幅度較大,辨識(shí)參數(shù)誤差為10.5%。所以,采用變增益系數(shù)法的永磁同步電機(jī)轉(zhuǎn)動(dòng)慣量辨識(shí),具有收斂速度快,辨識(shí)精度高,能有效抑制辨識(shí)結(jié)果振蕩的優(yōu)點(diǎn)。
本文在分析MRAS轉(zhuǎn)動(dòng)慣量辨識(shí)的基礎(chǔ)上,給出了一種變增益系數(shù)算法的研究和實(shí)現(xiàn)方法,并通過軟件仿真驗(yàn)證了理論的正確性。
可以得出,采用變增益系數(shù)法的轉(zhuǎn)動(dòng)慣量辨識(shí)系統(tǒng)在收斂速度方面比采用固定增益系數(shù)的轉(zhuǎn)動(dòng)慣量辨識(shí)收斂速度快;在穩(wěn)定性方面,轉(zhuǎn)動(dòng)慣量辨識(shí)精度可提高8%,穩(wěn)定性較好。因此,變增益系數(shù)法可同時(shí)獲得穩(wěn)定性與收斂速度兩方面的優(yōu)勢,具有收斂時(shí)間短、抑制辨識(shí)參數(shù)振蕩和參數(shù)辨識(shí)精度高等特點(diǎn)。
[1]丁信忠,張承瑞,李虎修,等.永磁同步電機(jī)的轉(zhuǎn)動(dòng)慣量辨識(shí)及狀態(tài)估計(jì)[J].山東大學(xué)學(xué)報(bào)(工學(xué)版),2012,42(2): 70-77.
[2]DU S, ZHAO S H, CHEN Y S. Inertia identification for speed control of PMSM servo motor[J]. College of Electrical Engineering, 2010,32(3): 1- 4.
[3]徐東,王田苗,魏洪興.一種基于簡化模型的永磁同步電機(jī)轉(zhuǎn)動(dòng)慣量辨識(shí)和誤差補(bǔ)償[J].電工技術(shù)學(xué)報(bào),2013,28(2): 126-131.
[4]盧少武,唐小琦,宋寶.伺服系統(tǒng)轉(zhuǎn)動(dòng)慣量辨識(shí)及其應(yīng)用[J].微電機(jī),2011,44(10): 41- 43.
[5]劉旭,阮毅,張朝藝.一種異步電機(jī)轉(zhuǎn)動(dòng)慣量辨識(shí)方法[J].電機(jī)與控制應(yīng)用,2013,40(5): 29-35.
[6]韓亞榮,邱鑫,朱德明.永磁交流伺服系統(tǒng)轉(zhuǎn)動(dòng)慣量辨識(shí)方法[J].電工電能新技術(shù),2013,32(3): 36- 40.
[7]梁驕雁.永磁同步電機(jī)轉(zhuǎn)動(dòng)慣量辨識(shí)的研究[D].南京: 南京航空航天大學(xué),2011.
[8]GUO Y J, HUANG L P, MURAMATSU M. Reserch on inertia identification and auto-tuning of speed controller for AC servo system[C]∥Proceedings of the Power Conversion Conference-Osaka, 2002: 896-901.
[9]毛丁輝,邱建琪,史涔溦.基于轉(zhuǎn)動(dòng)慣量的異步電機(jī)參數(shù)自整定系統(tǒng)研究[J].機(jī)電工程,2015,32(6): 830-835.
[10]夏加寬,盛麗君.轉(zhuǎn)臺(tái)系統(tǒng)轉(zhuǎn)動(dòng)慣量辨識(shí)與控制器參數(shù)自調(diào)整[J].沈陽工業(yè)大學(xué)學(xué)報(bào),2008,30(6): 605-609.
[11]NIU L, XU D G, YANG M, et al. On-line inertia identification algorithm for PI parameters optimization in speed loop[J]. IEEE Transactions on Power Electronics, 2015,30(2): 849-859.
[12]姚磊,邱鑫,王慧貞,等.永磁交流伺服系統(tǒng)轉(zhuǎn)動(dòng)慣量辨識(shí)方法[J].電機(jī)與控制應(yīng)用,2013,40(5): 29-35.
[13]姬宣德.基于模型參考自適應(yīng)的異步電動(dòng)機(jī)無速度傳感器DTC[J].微特電機(jī),2010,38(5): 63-67.
[14]張華強(qiáng),嚴(yán)洪峰,冷艷禮,等.基于模型參考自適應(yīng)的PMSM在線辨識(shí)研究[J].電機(jī)與控制應(yīng)用,2015,42(12): 3-6.
[15]王慶龍,張崇巍,張興.基于變結(jié)構(gòu)模型參考自適應(yīng)系統(tǒng)的永磁同步電機(jī)轉(zhuǎn)速辨識(shí)[J].中國電機(jī)工程學(xué)報(bào),2008,28(9): 71-75.
[16]沈艷霞,劉永欽.基于改進(jìn)型模型參考自適應(yīng)的PMSM參數(shù)辨識(shí)[J].電氣傳動(dòng),2009,39(5): 47-50.
[17]王偉華,肖曦.基于電感辨識(shí)的PMSM電流自適應(yīng)增量預(yù)測控制[J].電機(jī)與控制學(xué)報(bào),2014,18(2): 75-82.
Inertia Identification of Permanent Magnet Synchronous Motor Based on Improved Model Reference Adaptive System Algorithm
WANGFeiyu,TIANJingcheng,ZHUOKeqiong,ZHAOChaohui
(College of Electrical Engineering, Shanghai Dianji University, Shanghai 201306, China)
According to the algorithm’s problem which could not acquire performance both in the stability and response speed, introduced the principle of the inertia identification of permanent magnet synchronous motor using model reference adaptive, and studied the function of adaptive gain coefficient in PMSM’s rotational inertia identification algorithm. Present a novel method to identify the inertia which could change the gain coefficient in real time was presented, and the simulation in MATLAB/Simulink was given. The simulation results showed that the novel method had the performance of fast convergent speed, small oscillation amplitude results and good stability.
permanent magnet synchronous motor(PMSM); model reference adaptive system(MRAS); inertia identification
王飛宇(1990—),男,碩士,研究方向?yàn)殡姍C(jī)驅(qū)動(dòng)和控制。
趙朝會(huì)(1963—),男,博士,教授,研究方向?yàn)殡娏﹄娮蛹半娏鲃?dòng)。
TM 351
A
1673-6540(2016)08- 0063- 05
2016-03-23