,, , ,
(上海工程技術(shù)大學(xué) 電子電氣工程學(xué)院,上海 201620)
感應(yīng)耦合式電能傳輸(inductively coupled power transmission,ICPT)系統(tǒng)利用交變磁場實現(xiàn)電能的無線傳輸,其主要由交變電源、電磁耦合線圈、諧振補償網(wǎng)絡(luò)三部分構(gòu)成。和傳統(tǒng)的有線充電方式相比,ICPT系統(tǒng)采用非接觸方式實現(xiàn)電能的無線傳輸并具備安全、便捷、可靠、低維修成本等優(yōu)勢。隨著國內(nèi)外學(xué)者對ICPT系統(tǒng)中諧振補償網(wǎng)絡(luò)、磁耦合線圈、系統(tǒng)控制策略和輻射安全等方面的理論研究不斷深入,該系統(tǒng)在電動汽車、智能家居、器官移植等設(shè)備的充電系統(tǒng)中得到廣泛應(yīng)用[1]。
ICPT系統(tǒng)為高階非線性系統(tǒng),所以系統(tǒng)參數(shù)攝動和外部擾動等不確定因素都會對該系統(tǒng)的魯棒性產(chǎn)生不可避免的影響。為提高ICPT系統(tǒng)在不確定因素影響下的魯棒性,ICPT系統(tǒng)需通過一定的控制方法實現(xiàn)系統(tǒng)的閉環(huán)控制。例如,比例、積分、微分(proportion,integral and differential,PID)控制具有良好的信號跟蹤能力,然而在高階非線性系統(tǒng)中PID的控制參數(shù)很難確定,其控制結(jié)果不是很理想[2]。李硯玲等人采用μ控制方法提高π型ICPT系統(tǒng)對諧振頻率攝動的魯棒性[3],然而通過該方法設(shè)計的控制器階數(shù)較高,因此無法得到實際應(yīng)用。H∞優(yōu)化控制方法通過跟蹤外部擾動信號、不確定因素和實際系統(tǒng)與近似模型的誤差來提高系統(tǒng)的魯棒性,然而利用該方法分析閉環(huán)攝動系統(tǒng)的魯棒性較為保守。為避免系統(tǒng)魯棒性的分析結(jié)果過于保守,本文采用μ理論對系統(tǒng)的魯棒性進行分析。
ICPT系統(tǒng)中原邊和副邊線圈之間的空氣間隙較大,故采用諧振補償網(wǎng)絡(luò)來提高系統(tǒng)能量傳輸?shù)男?。諧振補償網(wǎng)絡(luò)主要有串/串聯(lián)(series series,SS)、串/并聯(lián)(series parallel,SP)、并/并聯(lián)(parallel parallel ,PP)、并/串聯(lián)(parallel series ,PS)四種形式[4]。其中,SS型諧振補償網(wǎng)絡(luò)的原邊和副邊的諧振補償電容均與電磁耦合線圈串聯(lián);SP型諧振補償網(wǎng)絡(luò)的原邊諧振補償電容與電磁耦合線圈串聯(lián),而副邊諧振補償電容與電磁耦合線圈并聯(lián);PP型諧振補償網(wǎng)絡(luò)的原邊和副邊諧振補償電容均與電磁耦合線圈并聯(lián);PS型諧振補償網(wǎng)絡(luò)的原邊諧振補償電容與電磁耦合線圈并聯(lián),而副邊諧振補償電容與電磁耦合線圈串聯(lián)。
本文選用SP型諧振拓撲作為ICPT系統(tǒng)的諧振補償網(wǎng)絡(luò)并利用廣義狀態(tài)空間平均(generalized state space averaging,GSSA)方程將該系統(tǒng)轉(zhuǎn)為線性系統(tǒng),采用線性分式變換(linear fractional transformation,LFT)方法將ICPT系統(tǒng)中的標稱和不確定因素分離并建立結(jié)構(gòu)化不確定性模型。在結(jié)構(gòu)化不確定性模型的基礎(chǔ)上利用H∞優(yōu)化控制方法設(shè)計魯棒控制器,仿真并分析ICPT閉環(huán)攝動系統(tǒng)的魯棒性。
圖1為SP型ICPT系統(tǒng)的電路拓撲,由圖可知交變電源主要由直流電源Edc和全橋逆變電路構(gòu)成,該電源產(chǎn)生的交變信號通過電磁耦合線圈產(chǎn)生交變磁場,從而將原邊產(chǎn)生的電能通過耦合線圈傳輸?shù)礁边?。原邊與副邊電磁耦合線圈的電感量分別為Lp和Ls,由于兩者存在一定的空氣間隙,互感M偏小使得電能傳輸過程中會有大量的能量損耗,因此系統(tǒng)原邊和副邊均采用諧振補償網(wǎng)絡(luò)來減少ICPT系統(tǒng)的無功功率,從而提高系統(tǒng)的能量傳輸效率,原邊諧振補償電容和副邊諧振補償電容的電容量分別為Cp和Cs。副邊諧振補償網(wǎng)絡(luò)輸出的電壓通過全橋整流電路、LC濾波電路給負載供電,該系統(tǒng)的負載可等效為電阻Rd。
圖1 SP型ICPT系統(tǒng)的電路拓撲
為實現(xiàn)ICPT系統(tǒng)原邊與副邊諧振補償網(wǎng)絡(luò)的相互獨立,可對該系統(tǒng)原邊與副邊耦合電路進行解耦。其中副邊到原邊的反射阻抗Zr見公式(1),公式中Zs為副邊電路阻抗,ω為系統(tǒng)工作頻率。由全橋整流電路、濾波電路和系統(tǒng)負載構(gòu)成的整體可等效為電阻Rl,該電阻的阻值見公式(2)。
(1)
(2)
根據(jù)上述理論得到的系統(tǒng)解耦等效電路如圖2所示。交變電源uac可等效為隨時間周期性變化的方波信號。原邊解耦等效電路中Rr和Lr構(gòu)成反射阻抗Zr。副邊等效電路通過 Thevinin定理實現(xiàn)電路拓撲的轉(zhuǎn)換,副邊解耦補償電路的電感Ls和電容Cs構(gòu)成并聯(lián)諧振補償電路。
ICPT系統(tǒng)中諧振補償網(wǎng)絡(luò)產(chǎn)生的諧振頻率對整個系統(tǒng)的工作狀態(tài)產(chǎn)生重要影響,只有當(dāng)該系統(tǒng)的諧振頻率與交變電源產(chǎn)生的交變頻率基本一致時,系統(tǒng)處于最佳工作狀態(tài)。此時ICPT系統(tǒng)的工作頻率即為系統(tǒng)諧振頻率ω0,故其原邊和副邊諧振補償網(wǎng)絡(luò)的諧振頻率見公式(3)。
(3)
當(dāng)該系統(tǒng)處于諧振工作狀態(tài)時,原邊諧振補償網(wǎng)絡(luò)可視為短路,副邊諧振補償網(wǎng)絡(luò)可視為斷路,副邊補償電容Cs兩端的電壓us的表達式見公式(4)。
(4)
ICPT系統(tǒng)對應(yīng)的等效電路如圖3所示。圖中EdcSi(t)為交變電源產(chǎn)生的方波信號,其中Si(t)的初相角記為0,其表達式如公式(5)所示。方波信號的交變頻率與全橋逆變電路中MOS管的開關(guān)頻率相同,其對應(yīng)的角頻率記為ω;原邊電磁耦合線圈的電感Lp和原邊諧振補償電容Ls串聯(lián)構(gòu)成原邊諧振補償網(wǎng)絡(luò),副邊電磁耦合線圈的電感Ls和副邊諧振補償電容Cs并聯(lián)構(gòu)成副邊諧振補償網(wǎng)絡(luò);電壓受控源jωMis和jωMip共同構(gòu)成電磁耦合線圈的互感等效模型;電流受控源idSr(t)和電壓受控源udSr(t)共同構(gòu)成全橋整流電路等效模型,其中Sr(t)的表達式見公式(6)。
(5)
(6)
將ICPT系統(tǒng)中流經(jīng)電感的電流和電容兩端的電壓作為狀態(tài)變量,得到的狀態(tài)空間方程如公式(7)的方程組所示。因為方程中存在Si(t)和Sr(t)這兩個線性時變量,故需建立ICPT系統(tǒng)的GSSA方程[5]來使其線性化。
(7)
為建立GSSA方程,首先需對系統(tǒng)中的狀態(tài)變量和線性時變量進行傅里葉變換。當(dāng)系統(tǒng)處于諧振工作狀態(tài)時,流經(jīng)原副邊諧振網(wǎng)絡(luò)的交變信號可視為正弦波,因此狀態(tài)變量的基波分量和一次諧波可近似為直流和交流分量,其中交流分量可通過實部和虛部來體現(xiàn)。該模型的狀態(tài)變量x見公式(8)。
x=[Re
lm
Re
(8)
由第一小節(jié)的公式(4)可知us和uac的初相角均為0,故這兩個線性時變量經(jīng)傅里葉變換后得到的一次傅里葉系數(shù)相同,即:
(9)
根據(jù)傅里葉變換的相關(guān)性質(zhì)和公式(7)建立ICPT系統(tǒng)的GSSA方程,其對應(yīng)的表達形式如公式(10)的方程組所示。其中該系統(tǒng)的輸入量u為直流電源Edc,輸出量y為負載電阻Rd的兩端電壓ud,廣義矩陣A、B、C、D構(gòu)成系統(tǒng)的廣義傳遞矩陣G為線性常數(shù)矩陣[6],故GSSA方程將該非線性系統(tǒng)轉(zhuǎn)為線性系統(tǒng)。
(10)
ICPT系統(tǒng)會受到系統(tǒng)參數(shù)攝動和外部擾動等不確定因素的影響,考慮到該系統(tǒng)的實際參數(shù)無法準確確定,故利用線性分式變換(LFT)的方法將系統(tǒng)的不確定因素從標稱因素分離出來,建立ICPT系統(tǒng)結(jié)構(gòu)化不確定性模型[7],其結(jié)構(gòu)化不確定性模型方程為:
Γ(Δ,M)=M22+M12(I-M11Δ-1)M21
(11)
其中:分塊矩陣M11、M12、M21、M22構(gòu)成ICPT系統(tǒng)的標稱矩陣M,該矩陣為系統(tǒng)標稱因素的集合,Δ為不確定因素構(gòu)成的系統(tǒng)攝動矩陣。
取負載電阻Rd的阻值變化作為ICPT系統(tǒng)的參數(shù)攝動,其攝動形式為:
Rd=Rdnom(1+pφ)
(12)
式中,Rdnom為負載電阻Rd的恒定值;p為負載電阻Rd的變化范圍,在本文中p值取1;φ為負載電阻Rd的攝動量,其變化范圍在-1和1之間。
由公式(7)可知負載電阻Rd的倒數(shù)為系統(tǒng)的狀態(tài)空間系數(shù),因此根據(jù)公式(11)對Rd的倒數(shù)進行LFT變換,其表達形式為:
(13)
將上述Rd攝動情況下的理論分析應(yīng)用在ICPT系統(tǒng)中便得到該系統(tǒng)在參數(shù)攝動影響下的上LFT的結(jié)構(gòu)化不確定模型,該模型如圖4所示。由圖可知,當(dāng)ICPT系統(tǒng)存在參數(shù)攝動這一不確定因素時,在原系統(tǒng)的輸入信號w和輸出信號z的基礎(chǔ)上新增了y和u作為擾動輸入和輸出信號。
圖4 ICPT系統(tǒng)的結(jié)構(gòu)化不確定性模型
圖中Δ為系統(tǒng)攝動矩陣,該擾動系統(tǒng)中Δ為Rd的攝動量φ并以對角矩陣的形式存在。考慮到Δ的具體形式未知,故需對ICPT系統(tǒng)標稱因素構(gòu)成的廣義標稱傳遞矩陣GSP建立如公式(14)所示的GSSA方程,其下標SP表示ICPT系統(tǒng)采用SP型諧振補償網(wǎng)絡(luò)。
(14)
根據(jù)分塊矩陣的特性可將GSP分成:
(15)
將公式(15)代入公式(11)便得到該擾動系統(tǒng)的參考模型,即:
Γ(Δ,GSP)=GSP-22+GSP-12(I-GSP-11Δ-1)GSP-21
(16)
通過上述的理論分析得到的系統(tǒng)參數(shù)如表1所示。
表1 ICPT系統(tǒng)電氣參數(shù)
外部擾動對ICPT系統(tǒng)的影響不容忽視,故對ICPT系統(tǒng)引入系統(tǒng)參數(shù)攝動和外部擾動等不確定因素并進行閉環(huán)控制[8],其設(shè)計框圖如圖5所示。由圖可知,PSP為擾動系統(tǒng)的廣義標稱對象,其結(jié)果可根據(jù)第二小節(jié)的GSP理論分析計算得出;Δ是系統(tǒng)攝動矩陣,其大小可以通過矩陣范數(shù)來衡量;ref為系統(tǒng)的參考量,其與實際輸出量的偏差通過控制器K進行控制來保證系統(tǒng)的穩(wěn)定輸出。
圖5 ICPT閉環(huán)攝動系統(tǒng)的設(shè)計框圖
d為系統(tǒng)的外部擾動且其具體形式未知,若僅用常數(shù)來表示擾動引起輸出信號的不確定性,則系統(tǒng)全頻帶不確定性的幅值變化相同,這樣擴大了系統(tǒng)的擾動范圍從而無法得到高性能的控制效果。因此ICPT系統(tǒng)進行閉環(huán)控制時,需要采用加權(quán)函數(shù)[9]Wp和Wu對擾動輸出信號的范圍進行限定,得到:
(17a)
(17b)
由于系統(tǒng)在不確定因素的作用下,在PSP變化的同時,也對閉環(huán)傳遞函數(shù)TSP造成影響,故將PSP的變化量對TSP的影響程度視為系統(tǒng)敏感度[6]S,其定義為:
(18)
由于敏感度和加權(quán)函數(shù)對系統(tǒng)誤差信號跟蹤及干擾抑制等指標具有重要影響,而這些指標對系統(tǒng)的穩(wěn)定性和控制性能的評估起到?jīng)Q定性作用,故由此可得到系統(tǒng)標稱穩(wěn)定和性能指標判定條件為:
(19)
考慮到系統(tǒng)不確定因素的大小是由攝動矩陣的無窮范數(shù)來衡量的,因此本文采用H∞優(yōu)化控制方法設(shè)計閉環(huán)控制器,用該方法設(shè)計的控制器不僅滿足如公式(19)所示的判定條件,還需保證ICPT系統(tǒng)中外部干擾輸入d到評估輸出e_p、e_u的傳遞矩陣無窮范數(shù)小于給定值γ??刂破鱇的設(shè)計流程如圖6所示。
圖6 ICPT系統(tǒng)的閉環(huán)控制器設(shè)計流程圖
為驗證閉環(huán)控制器的控制性能,需對ICPT閉環(huán)控制系統(tǒng)的魯棒性進行分析。為方便實際系統(tǒng)的魯棒性分析,可將圖6所示的閉環(huán)控制結(jié)構(gòu)框圖簡化為系統(tǒng)攝動矩陣Δ和標稱閉環(huán)傳遞矩陣Gclp兩部分并構(gòu)成ICPT閉環(huán)系統(tǒng)的結(jié)構(gòu)化不確定性模型,故該系統(tǒng)魯棒性的判定條件見公式(19)。根據(jù)系統(tǒng)攝動矩陣Δ不確定因素的不同,可通過魯棒穩(wěn)定性和魯棒性能[6]兩方面來綜合衡量系統(tǒng)的魯棒性。為避免結(jié)果的保守性和不必要性,采用μ理論計算Gclp的結(jié)構(gòu)奇異值大小作為該系統(tǒng)魯棒性的判定依據(jù)并記為μΔ(Gclp)[10]。
(20)
當(dāng)系統(tǒng)攝動矩陣Δ只有系統(tǒng)參數(shù)攝動引起的不確定因素時,得到的μΔ(Gclp)大小可作為系統(tǒng)魯棒穩(wěn)定性的衡量標準。圖7為閉環(huán)系統(tǒng)的魯棒穩(wěn)定性,其中的虛線和實線分別代表μ值隨頻率變化的上限和下限變化趨勢。
圖7 ICPT閉環(huán)系統(tǒng)的魯棒穩(wěn)定性
由圖可知,當(dāng)角頻帶在102~105rad/s范圍內(nèi),因系統(tǒng)參數(shù)攝動和外部擾動的作用下μ值在10-3~1之間呈現(xiàn)不規(guī)則變化趨勢。在該范圍內(nèi)μ的最大值為0.91,即系統(tǒng)自身攝動矩陣的無窮范數(shù)滿足小于1/0.91這一條件時,μ值小于0.91。由第二小節(jié)可知系統(tǒng)自身攝動矩陣的無窮范數(shù)小于1,故ICPT系統(tǒng)中表征系統(tǒng)魯棒穩(wěn)定性的μ值均小于1,證明控制器K能夠保證系統(tǒng)的魯棒穩(wěn)定性。
ICPT系統(tǒng)的魯棒性能可通過系統(tǒng)擾動對象集合的干擾抑制、誤差信號跟蹤等方面來體現(xiàn)。當(dāng)系統(tǒng)攝動矩陣Δ在系統(tǒng)自身不確定因素的基礎(chǔ)上加入性能不確定因素構(gòu)成增廣攝動矩陣時,得到的μΔ(Gclp)大小可作為系統(tǒng)魯棒性能的衡量標準。
圖8表征閉環(huán)系統(tǒng)的魯棒性能,其中的虛線和實線分別代表表征魯棒性能的上限和下限變化趨勢。由圖可知,當(dāng)角頻帶在102~104rad/s范圍內(nèi),因系統(tǒng)自身攝動的作用下μ值呈現(xiàn)不規(guī)則變化趨勢,而當(dāng)角頻帶在104~105rad/s范圍內(nèi)μ值為恒定值,由此說明在加權(quán)函數(shù)和控制器的共同作用下,可抑制外部擾動d對系統(tǒng)的干擾。在102~105rad/s角頻帶范圍內(nèi),μ的最大值為0.92,即系統(tǒng)增廣攝動矩陣的無窮范數(shù)小于1/0.92這一條件時μ值小于0.92,故該系統(tǒng)有效地實現(xiàn)誤差信號的跟蹤,并達到魯棒性能的要求。
圖8 ICPT閉環(huán)系統(tǒng)的魯棒性能
為更好地驗證控制器K對系統(tǒng)魯棒性的影響,本文仿真得到該系統(tǒng)在不確定因素影響下的暫態(tài)響應(yīng)并進行分析。圖9為ICPT系統(tǒng)的開環(huán)和閉環(huán)暫態(tài)響應(yīng)。從圖中可以看出系統(tǒng)在開環(huán)狀態(tài)下的響應(yīng)時間為0.03 s且輸出電壓的超調(diào)量變大。加入控制器K后,輸出電壓的響應(yīng)時間縮短為0.01 s左右且沒有超調(diào)量,同時當(dāng)負載電壓達到穩(wěn)定狀態(tài)時的電壓紋波減小,由此可見控制器K實現(xiàn)系統(tǒng)暫態(tài)響應(yīng)的穩(wěn)定輸出。
圖9 ICPT開環(huán)和閉環(huán)系統(tǒng)的暫態(tài)響應(yīng)圖
本文通過H∞優(yōu)化控制方法設(shè)計的魯棒控制器并建立ICPT閉環(huán)控制系統(tǒng)。首先采用GSSA方程將基于SP諧振補償網(wǎng)絡(luò)的ICPT系統(tǒng)線性化,然后通過LFT方法將ICPT系統(tǒng)中不確定因素從標稱因素中分離并建立該系統(tǒng)的結(jié)構(gòu)化不確定性模型??紤]到系統(tǒng)不確定因素的大小可通過系統(tǒng)攝動矩陣的無窮范數(shù)來衡量,故采用H∞優(yōu)化控制方法設(shè)計ICPT的魯棒控制器,該控制器可消除系統(tǒng)參數(shù)攝動和外部擾動對ICPT系統(tǒng)的影響,與此同時,在該控制器基礎(chǔ)上建立的ICPT閉環(huán)控制系統(tǒng)即使在不確定因素的影響下也能保證標稱閉環(huán)傳遞矩陣的結(jié)構(gòu)奇異值小于1,具有很好的魯棒性能。基于H∞魯棒控制的ICPT閉環(huán)控制系統(tǒng)可廣泛應(yīng)用于電動汽車、智能家居等設(shè)備的電源充電領(lǐng)域,在提高電能傳輸效率的同時,也能提高系統(tǒng)輸出的電能質(zhì)量。