(蘇州大學(xué) 機電工程學(xué)院,蘇州 215021)
永磁同步電機(PMSM)因具有效率高、環(huán)保節(jié)能、高效能、結(jié)構(gòu)類型多樣等特點,在航空航天、船舶電力推進、風(fēng)機水泵、電梯、空調(diào)等領(lǐng)域行業(yè)都有著極大的應(yīng)用[1,2]。由于PI控制器的結(jié)構(gòu)簡單、實現(xiàn)方便、性能優(yōu)良、且對控制對象參數(shù)變化不敏感,因此被大量應(yīng)用于PMSM驅(qū)動系統(tǒng)[3]。而PI控制器參數(shù)的設(shè)置將直接影響控制器的控制效果,如果參數(shù)設(shè)置不當(dāng),系統(tǒng)的性能將大打折扣。
傳統(tǒng)的參數(shù)整定方法,如Z-N法,臨界比例法等,容易陷入局部最優(yōu),且調(diào)參時需要人工進行大量試驗,效率較低,且在試驗過程中,試驗人員往往以經(jīng)驗為主,根據(jù)過往調(diào)參時所使用過的調(diào)參經(jīng)驗進行試湊,因此得到的試驗結(jié)果會帶有較強的主觀性和隨機性,整定出來的控制器性能通常也會較差,無法滿足性能的要求[4]。隨著控制理論的不斷發(fā)展,出現(xiàn)了一些采用遺傳算法、比例增益法、卡爾曼濾波和Elman神經(jīng)網(wǎng)絡(luò)結(jié)合法來整定速度環(huán)控制器參數(shù)的方法,但在線的高階矩陣求逆計算會占用大量的系統(tǒng)資源,對硬件要求較高,大大限制了這些方法的應(yīng)用[6]。為此,本文基于速度環(huán)控制器的數(shù)學(xué)模型,以相角裕度作為衡量系統(tǒng)穩(wěn)定性的指標(biāo),通過系統(tǒng)的開環(huán)傳遞函數(shù)來計算出控制器的參數(shù),并在保證系統(tǒng)穩(wěn)定的前提下,進一步加快系統(tǒng)的響應(yīng),對于速度環(huán)控制器參數(shù)整定方法的研究具有十分重要的意義。
由于永磁同步電機伺服系統(tǒng)是一個耦合性較強的非線性系統(tǒng),為了簡化永磁同步電機數(shù)學(xué)模型的推導(dǎo)過程,須作以下五點假設(shè):
1)電機定子繞組的三相在空間上對稱分布;
2)空間磁場呈正弦分布;
3)不計磁滯和渦流損耗;
4)忽略鐵心飽和;
5)忽略齒槽轉(zhuǎn)矩及其他擾動影響。
根據(jù)上述五點假設(shè),采用id=0的磁場定向控制策略,則電機的電壓方程如下:
式中,Ud和Uq分別d軸和q軸定子電壓,id和iq分別d軸和q軸定子電流,Ld和Lq分別d軸和q軸定子繞組電感,Rs為定子繞組電阻,ωr為電機電角速度,Ψf為電機永磁體磁鏈。
電機的轉(zhuǎn)矩方程為:
式中,Te為電機的電磁轉(zhuǎn)矩,TL為負(fù)載轉(zhuǎn)矩,J為負(fù)載轉(zhuǎn)動慣量,ω為電機機械角速度,P為電機的磁極對數(shù),B為粘滯阻尼。
常見的速度環(huán)控制器結(jié)構(gòu)為雙閉環(huán)結(jié)構(gòu),內(nèi)環(huán)為電流環(huán),速度環(huán)控制器的系統(tǒng)方塊圖如圖1所示。其中,速度控制器為PI控制器,傳遞函數(shù)為:
式中,Kp為PI控制器的比例增益系數(shù),Ki為PI控制器的積分時間常數(shù)。
圖1 速度環(huán)控制器的系統(tǒng)方塊圖
電流環(huán)作為速度環(huán)控制器的內(nèi)環(huán)部分,由于其帶寬遠(yuǎn)遠(yuǎn)高于速度環(huán)的期望截止頻率,因此可以將電流環(huán)的閉環(huán)傳遞函數(shù)等效為單位增益的慣性環(huán)節(jié)[7],其傳遞函數(shù)為:
式中,ωc為電流環(huán)的帶寬。
由式(4)可知,電機電流到電磁轉(zhuǎn)矩的傳遞函數(shù)為:
忽略電機軸上粘滯阻尼對系統(tǒng)的影響,由式(3)可知,轉(zhuǎn)矩到轉(zhuǎn)速的傳遞函數(shù)為:
為了抑制速度信號中所混入的高頻噪聲對系統(tǒng)的干擾,需要對速度信號做濾波處理,通常選用數(shù)字低通濾波來實現(xiàn),其傳遞函數(shù)為:
式中,ωf為數(shù)字濾波器的帶寬。
因此,整個速度環(huán)控制器的開環(huán)傳遞函數(shù)為:
其系統(tǒng)結(jié)構(gòu)圖如圖2所示。
圖2 速度環(huán)控制器的系統(tǒng)結(jié)構(gòu)圖
通過觀察速度環(huán)控制器的開環(huán)傳遞函數(shù)表達式,發(fā)現(xiàn)在ω=Ki處存在一個零點,記為Z1;在ω=0處存在兩個極點,分別記為P1、P2;在ω=ωf處存在一個極點,記為P3,在ω=ωc處存在一個極點,記為P4,該系統(tǒng)為最小相位系統(tǒng)。根據(jù)零點和極點位置分布情況的不同,系統(tǒng)的工作狀態(tài)及Bode圖圖像也會不同,這里假設(shè)電流環(huán)的帶寬總是大于速度濾波器的帶寬,即極點P4處的頻率大于極點P3處的頻率,總共會有以下三種可能情況:
1)零點Z1處頻率低于極點P3處頻率;
2)零點Z1處頻率介于極點P3以及極點P4處頻率中間;
3)零點Z1處頻率大于極點P4處頻率。
對于情況一,其Bode圖趨勢大致如圖3所示,通過分析可知,當(dāng)零點Z1處頻率低于極點P3處頻率時,由于交接頻率ω=Ki處微分環(huán)節(jié)的作用,相頻特性曲線從-180°開始表現(xiàn)為上升的趨勢,之后又由于兩個慣性環(huán)節(jié)的作用,曲線開始慢慢下降,因此存在一個極大值點,只要增益設(shè)置得當(dāng),就可以得到充足的相角裕度來保證系統(tǒng)的穩(wěn)定性,一旦增益設(shè)置的太大或太小,系統(tǒng)工作都將會不穩(wěn)定。
圖3 零點Z1處頻率低于極點P3處頻率Bode圖
對于情況二,其Bode圖趨勢大致如圖4所示,通過分析可知,當(dāng)零點Z1處頻率介于極點P3以及極點P4處頻率中間時,由于交接頻率ω=ωf處慣性環(huán)節(jié)的作用,相頻特性曲線從-180°開始便呈現(xiàn)出下降的趨勢,即便在之后會因微分環(huán)節(jié)而帶來曲線走勢的部分上升,但相角裕度始終小于零,系統(tǒng)呈現(xiàn)出一種不穩(wěn)定的工作狀態(tài)。
圖4 零點Z1處頻率介于極點P3以及極點P4處頻率中間Bode圖
對于情況三,其Bode圖趨勢大致如圖5所示,通過分析可知,當(dāng)零點Z1處頻率大于極點P4處頻率中間時,和情況二類似,由于交接頻率ω=ωf處慣性環(huán)節(jié)的作用,相頻特性曲線從-180°開始便呈現(xiàn)出一種下降的趨勢,最后經(jīng)過微分環(huán)節(jié)作用曲線開始上升,并最終趨于-270°,因此該曲線存在一個極小值點,由于相角裕度始終小于零,系統(tǒng)呈現(xiàn)出一種不穩(wěn)定的工作狀態(tài)。
三種情況中,只有情況一系統(tǒng)可以穩(wěn)定工作。
圖5 零點Z1處頻率大于極點P4處頻Bode圖
相角裕度值的大小決定了系統(tǒng)的穩(wěn)定性,對于PI控制器的調(diào)試,通常要求增益穿越頻率處于開環(huán)相位的峰值附近[10],由于電流環(huán)帶寬一般遠(yuǎn)大于濾波器的帶寬,因此在進行穩(wěn)定性分析時,可考慮忽略極點P4所帶來的影響,因此速度環(huán)控制器的開環(huán)傳遞函數(shù)可簡化為:
為了滿足使增益穿越頻率處于開環(huán)相位的峰值附近這一設(shè)計條件,現(xiàn)假設(shè)零點Z1處頻率和極點P3處頻率關(guān)系為:
則有:
由于微分環(huán)節(jié)和慣性環(huán)節(jié)的相頻特性分別為arctanTω和-arctanTω,對于速度環(huán)控制器,其簡化過后的開環(huán)傳遞函數(shù)的相頻特性曲線可表示為:
解得:
故速度環(huán)控制器的增益穿越頻率的位置應(yīng)處于Bode圖中零點Z1處頻率所在位置和極點P3處頻率所在位置的中間,因為增益穿越頻率處的增益為單位增益,故有:
通過對式(17)求解,可得:
將式(16)所解得的ω代入式(14),可得相角裕度為:
上述過程中求得的Kp,Ki即為整定參數(shù),這兩個參數(shù)的適用條件為電流環(huán)帶寬遠(yuǎn)大于濾波器的帶寬,一般為5倍或者5倍以上頻率。將整定得到的參數(shù)代入速度環(huán)控制器的傳遞函數(shù),得到整定后控制器Bode圖如圖6所示,此時相角裕度為極大值。
通過對式(20)分析可知,隨著τ取值的不斷增大,相角裕度的極大值也會隨之增大,即系統(tǒng)的穩(wěn)定性也會隨之加強。
圖6 整定后的速度環(huán)控制器Bode圖
采用Maxon的EC-i 40系列449464型號電機作為被控對象,該款電機的特性如表1所示。
表1 電機特性參數(shù)
在Simulink中搭建速度環(huán)控制器系統(tǒng),模型系統(tǒng)按照圖2所示的速度環(huán)控制器來搭建,如圖7所示。搭建的如圖8所示的實驗平臺,按圖中所示編號順序依次為1號:DSP控制板,2號:開關(guān)電源,3號:電機驅(qū)動板,4號:Maxon電機。
實驗中選擇電流環(huán)的帶寬和濾波器的帶寬分別為3000Hz和500Hz,為了展示控制效果,分別進行在相同負(fù)載情況下選取不同τ值以及在相同τ值情況下選取不同負(fù)載一共兩組參數(shù)整定仿真和實驗。
選用轉(zhuǎn)動慣量為100gcm2的負(fù)載,并選用τ值為5,10,20進行實驗,整定得到的參數(shù)如表2所示。
表2 100gcm2負(fù)載下不同τ取值下的整定參數(shù)
圖7 Simulink中搭建的速度環(huán)控制器
圖8 速度環(huán)控制器實驗平臺
通過Simulink仿真和實驗所得到的速度階躍波形分別如圖9(a)和圖9(b)所示,其中,藍(lán)色曲線為τ=5時的響應(yīng)曲線,紅色曲線為τ=10時的響應(yīng)曲線,綠色曲線為τ=20時的響應(yīng)曲線。通過對圖像的分析可知,無論是用仿真還是實驗,所得到的階躍響應(yīng)波形均具有如下特性:
1)隨著τ值的增加,系統(tǒng)的超調(diào)量會隨之減小;
2)隨著τ值的增加,系統(tǒng)的上升時間隨之增大。
從這兩個特性中可以看出,即隨著τ取值的不斷增加,系統(tǒng)的相角裕度會增大,因此系統(tǒng)的穩(wěn)定性就得到了提高,但相應(yīng)的系統(tǒng)的響應(yīng)速度會有一些下降,這與最初的分析相符。
圖9 速度階躍波形
選用τ值為20,轉(zhuǎn)動慣量為100gcm2和500gcm2的負(fù)載進行實驗,整定得到的參數(shù)如表3所示。
表3 τ值為20下不同負(fù)載的整定參數(shù)
通過Simulink仿真和實驗所得到的速度階躍波形分別如圖10(a)和圖10(b)所示,其中,紅色曲線為轉(zhuǎn)動慣量為100gcm2的響應(yīng)曲線,藍(lán)色曲線為轉(zhuǎn)動慣量為500gcm2的響應(yīng)曲線,通過觀察曲線可知,當(dāng)負(fù)載加大時,系統(tǒng)的上升時間變長,即系統(tǒng)的響應(yīng)速度變慢,如果想要提高系統(tǒng)的響應(yīng)速度,可以適當(dāng)?shù)臏p小τ的取值以此來提高系統(tǒng)的響應(yīng)速度,不過這樣會犧牲部分系統(tǒng)的穩(wěn)定性。
圖10 速度階躍波形
本文提出了一種基于速度環(huán)控制器結(jié)構(gòu)模型,采用頻率響應(yīng)的設(shè)計方法來對PMSM伺服系統(tǒng)的速度環(huán)控制器參數(shù)進行整定的方法,通過對速度環(huán)控制器傳遞函數(shù)的零點和極點分布情況的不同,討論了可使系統(tǒng)工作穩(wěn)定的零點和極點的分布情況,并依據(jù)使系統(tǒng)相角裕度最大的設(shè)計原則對速度環(huán)控制器的參數(shù)進行了設(shè)計和整定。仿真和實驗表明,采用這套方法所整定的控制器參數(shù)具有可靠性高、響應(yīng)快的特點,對于電機轉(zhuǎn)速具有良好的控制效果,在實際工程應(yīng)用中具有一定的應(yīng)用價值。