(陜西科技大學(xué) 電氣與信息工程學(xué)院,西安 710021)
在電力電子設(shè)備實(shí)驗(yàn)與研發(fā)過程中,控制器是最為核心也最關(guān)鍵的一部分,目前高性能的電力電子設(shè)備大多采用數(shù)字控制,在微控制芯片上利用軟件代碼實(shí)現(xiàn)控制算法。隨著對(duì)系統(tǒng)性能要求不斷提高,導(dǎo)致控制算法愈趨復(fù)雜,傳統(tǒng)的人工編寫代碼模式漸漸顯示出不足,如:開發(fā)周期長,代碼測試與驗(yàn)證費(fèi)時(shí)費(fèi)力等。為了解決傳統(tǒng)研發(fā)過程中的這些問題,提高開發(fā)效率,MathWorks公司提出了基于模型的設(shè)計(jì)方法[1],在MATLAB中設(shè)計(jì)完成控制器后可以自動(dòng)生成實(shí)現(xiàn)控制器的軟件代碼,與底層驅(qū)動(dòng)代碼集成,在微處理器上實(shí)現(xiàn)控制器,省略了人工編寫控制器代碼的步驟,并且可以很方便地對(duì)代碼進(jìn)行測試和驗(yàn)證??梢允寡芯咳藛T集中精力在控制算法的研究和控制器的設(shè)計(jì)上,而不需要投入過多的精力去關(guān)心控制器的實(shí)現(xiàn)過程,采用基于模型的設(shè)計(jì)方法可以提高代碼的可靠性、可移植性、易測試及易讀性,也可以將整個(gè)開發(fā)周期縮短至傳統(tǒng)開發(fā)模式的20%~50%[2]。
針對(duì)目前研究熱點(diǎn)光伏并網(wǎng),利用基于模型設(shè)計(jì)的方法,以單相光伏并網(wǎng)逆變器為例,從其基本原理出發(fā)直到最終的實(shí)際實(shí)現(xiàn),展示了基于模型的設(shè)計(jì)方法的整個(gè)過程及其優(yōu)點(diǎn)。
如圖1所示,單相光伏并網(wǎng)逆變系統(tǒng)主要包含:電網(wǎng)、光伏電池板、單相逆變橋、濾波器及控制器5部分。單相并網(wǎng)逆變系統(tǒng)的目的是將太陽能轉(zhuǎn)化為電流的形式注入電網(wǎng),且并網(wǎng)電流與電網(wǎng)電壓同頻同相,即實(shí)現(xiàn)單位功率因數(shù),向電網(wǎng)輸送有功電流。
圖1 系統(tǒng)框圖Fig.1 System diagram
其中:Udc為直流側(cè)電壓為MPPT(最大功率點(diǎn)跟蹤)的輸出值,也是直流側(cè)電壓給定值為有功電流幅值給定值;eac為電網(wǎng)電壓為并網(wǎng)電流瞬時(shí)給定值;il為實(shí)際并網(wǎng)電流。虛線框中的部分為控制器部分,在控制芯片上由軟件實(shí)現(xiàn)。
根據(jù)檢測直流側(cè)電壓電流,用電壓擾動(dòng)法,改變直流側(cè)給定電壓與上次輸出功率比較,從而決定電壓的改變方向,最終實(shí)現(xiàn)MPPT;MPPT輸出的給定電壓與直流側(cè)實(shí)際電壓做差,再通過PI控制器,實(shí)現(xiàn)直流側(cè)電壓的控制,應(yīng)當(dāng)注意的正負(fù),其正負(fù)號(hào)與輸出的并網(wǎng)電流方向有關(guān),總之應(yīng)保證當(dāng)增大時(shí)應(yīng)減小,即根據(jù)功率平衡來判斷;對(duì)電網(wǎng)電壓進(jìn)行鎖相可以得到與eac同頻同相的單位正弦波,與相乘可以得到瞬時(shí)電流給定值,最后通過無差拍控制器,使實(shí)際并網(wǎng)電流跟蹤給定電流,實(shí)現(xiàn)單位功率因數(shù)并網(wǎng)。
單相光伏并網(wǎng)逆變器常用的電流控制器有PI、電流滯環(huán)控制、準(zhǔn)比例諧振控制等。PI控制器對(duì)交流信號(hào)不能實(shí)現(xiàn)無靜差跟蹤;電流滯環(huán)控制器雖具有動(dòng)態(tài)快、誤差小等優(yōu)點(diǎn),但其開關(guān)頻率不固定而且不利于數(shù)字化實(shí)現(xiàn)及濾波設(shè)計(jì)[3];準(zhǔn)比例諧振控制器可以實(shí)現(xiàn)無靜差跟蹤,但動(dòng)態(tài)響應(yīng)慢。針對(duì)本系統(tǒng)的特點(diǎn),設(shè)計(jì)了無差拍控制器,其具有算法簡單,動(dòng)態(tài)響應(yīng)快、精度高等優(yōu)點(diǎn)[4-5]。
圖1中系統(tǒng)電壓電流關(guān)系為
其中:Uab為逆變橋兩橋臂之間的輸出電壓;L為濾波電感;R為濾波電感內(nèi)阻 (包括導(dǎo)線電阻);eac為電網(wǎng)電壓;il為電感電流也是并網(wǎng)電流。
采用雙極性調(diào)制,有:
其中:d為占空比;Udc為直流側(cè)電壓。將式(2)帶入式(1)并離散化可得:
其中:Ts為采樣時(shí)間;il(k+1)為系統(tǒng)下一時(shí)刻的給定電流值;il(k)為當(dāng)前電流值,因此可由式(3)計(jì)算出占空比信號(hào)從而實(shí)現(xiàn)對(duì)并網(wǎng)電流的控制。
傳統(tǒng)開發(fā)模式大多是先在MATLAB環(huán)境下搭建系統(tǒng)仿真模型,設(shè)計(jì)控制器并仿真驗(yàn)證,驗(yàn)證通過后人工編寫代碼,然后對(duì)代碼進(jìn)行測試并修復(fù)bug,最后在實(shí)際系統(tǒng)中測試。這樣開發(fā)效率就比較低下,因?yàn)槿斯懘a往往會(huì)引入bug,而且代碼測試不方便。
而基于模型的設(shè)計(jì)方法是建立好仿真模型并測試通過后,可以對(duì)控制器自動(dòng)生成代碼,而且可以通過SIL(軟件在環(huán))測試生成的代碼與MATLAB中設(shè)計(jì)的控制器是否等效,可以通過PIL(處理器在環(huán))測試代碼的執(zhí)行時(shí)間及效率,還可以通過MAT-LAB提供的工具測試覆蓋度及代碼復(fù)雜度等,而這些對(duì)于人工代碼是很難完成的。
在搭建系統(tǒng)模型的時(shí)候應(yīng)根據(jù)實(shí)際系統(tǒng)選取仿真中的參數(shù),如圖2所示,系統(tǒng)硬件部分由電網(wǎng)、逆變橋、濾波電感及直流電源組成,為了實(shí)現(xiàn)簡單采用直流電壓源代替了光伏電池板,因此并網(wǎng)有功電流幅值為外部直接給定;軟件部分由控制器(Controller)和脈寬調(diào)制(PWM)模塊組成,最終均要在微控制器上實(shí)現(xiàn),其中控制器部分就是將來要自動(dòng)生成代碼的部分,脈寬調(diào)制模塊由底層程序?qū)崿F(xiàn)。
圖2 系統(tǒng)仿真模型Fig.2 System simulation model
其中,采樣時(shí)間Ts=0.00005 s,濾波電感內(nèi)阻R=0.26 Ω,濾波電感L=3 mH,電網(wǎng)電壓有效值eac= 220 V,逆變器直流側(cè)電壓Udc=350 V,并網(wǎng)逆變器的額定功率為1.5 kW,額定并網(wǎng)電流幅值為10 A。
圖3為控制器部分的內(nèi)部結(jié)構(gòu)主要由鎖相環(huán)和無差拍控制器組成,其中無差拍控制器的表達(dá)式如式(3)所示。
圖3 控制器內(nèi)部結(jié)構(gòu)Fig.3 Inner structure of controller
圖4為并網(wǎng)電流及電網(wǎng)電壓波形,為了使并網(wǎng)電流和電網(wǎng)電壓能更好地一起顯示,故將并網(wǎng)電流放大了20倍與電網(wǎng)電壓同時(shí)顯示,圖中橫坐標(biāo)代表時(shí)間軸,縱坐標(biāo)代表電壓和電流軸,電壓100 V/格,電流5 A/格;其中t<0.04 s時(shí)并網(wǎng)電流幅值給定值可以看出并網(wǎng)電流與電網(wǎng)電壓同頻同相,且在t=0.04 s時(shí)并網(wǎng)電流能快速響應(yīng)幅值變?yōu)?0 A的正弦波電流。
圖4 并網(wǎng)電流及電網(wǎng)電壓波形Fig.4 Current and voltage waveform
圖5是5 A和10 A并網(wǎng)電流的諧波分析結(jié)果,圖中橫坐標(biāo)代表頻率,單位是Hz;縱坐標(biāo)為諧波電流相對(duì)于基波電流的相對(duì)值。圖5(a)為時(shí)并網(wǎng)電流的諧波分析結(jié)果,其電流THD(總諧波畸變率)=3.94%,圖5(b)為時(shí)并網(wǎng)電流的諧波分析結(jié)果,其并網(wǎng)電流THD=1.95%,均滿足國標(biāo)的THD<5%的要求。
圖5 并網(wǎng)電流諧波分析Fig.5 Harmonic analysis of current
搭建完系統(tǒng)模型并仿真驗(yàn)證通過后即可對(duì)系統(tǒng)進(jìn)行配置自動(dòng)生成代碼,對(duì)代碼生成的相關(guān)關(guān)鍵配置選項(xiàng)如表1所示。
表1 代碼生成的相關(guān)配置Tab.1 Configuration of the code generation
分別代表了定步長離散求解器步長為0.00005 s;嵌入式目標(biāo)文件;模塊化代碼格式;TI公司c2000系列處理器;代碼以高執(zhí)行效率為優(yōu)先目標(biāo)。配置完以上選項(xiàng)后還應(yīng)對(duì)代碼原型進(jìn)行控制,以更好地與底層程序接口。
可以分別配置代碼中初始化函數(shù)名,函數(shù)名及輸入輸出參數(shù)如圖6所示。
圖6 函數(shù)原型控制界面Fig.6 Interface of function prototype control
在配置完以上選項(xiàng)后就可以生成代碼,在生成的代碼中只需要看其中的.h文件,找到函數(shù)名,通過函數(shù)名,調(diào)用可與底層代碼集成。最終生成的代碼中函數(shù)原型如圖7所示,與期望的函數(shù)原型一致。
圖7 函數(shù)原型Fig.7 Function prototype
基于模型設(shè)計(jì)的最大優(yōu)點(diǎn)就是容易實(shí)現(xiàn)代碼的驗(yàn)證與測試,在與底層代碼集成之前需要先驗(yàn)證MATLAB自動(dòng)生成的代碼與搭建的控制器是等效的。在配置選項(xiàng)中選中SIL,即可產(chǎn)生一個(gè)S-function模塊,是自動(dòng)生成的 C代碼經(jīng)過了封裝,可以在MATLAB中調(diào)用。因此可以對(duì)比自動(dòng)生成的代碼和設(shè)計(jì)的控制器是否等效,給其和控制器相同的輸入判斷兩者的輸出之差是否為零,即可判斷兩者是否等效。如圖8所示,SIL模塊為經(jīng)過封裝自動(dòng)生成的C代碼,經(jīng)測試其與搭建的控制在相同輸入的條件下,兩者的輸出差值始終為零,即兩者等效。之后即可與底層代碼集成,在實(shí)際系統(tǒng)中實(shí)驗(yàn)測試。
圖8 代碼等效性測試Fig.8 Equivalence test of code
按照以上仿真參數(shù)搭建實(shí)際系統(tǒng),將MATLAB生成的代碼與底層代碼集成,在TI公司的TMS320F28335微處理器上實(shí)現(xiàn)控制器。實(shí)驗(yàn)結(jié)果如圖9(a)、(b)所示,圖中橫坐標(biāo)代表時(shí)間軸,縱坐標(biāo)代表電壓和電流軸,電壓100 V/格,電流7.5 A/格;分別為實(shí)測電網(wǎng)電壓與幅值為5 A和10 A的并網(wǎng)電流波形,可以看出電壓電流同頻同相,實(shí)現(xiàn)了單位功率因數(shù)并網(wǎng),與仿真結(jié)果基本相同。
圖9 實(shí)測電壓電流波形Fig.9 Measured voltage and current waveform
采用基于模型的設(shè)計(jì)方法設(shè)計(jì)了單相光伏并網(wǎng)逆變器,利用MATLAB自動(dòng)生成了代碼,然后再對(duì)生成的代碼做SIL測試其與控制器的等效性,最后與底層集成并搭建硬件實(shí)現(xiàn)了并網(wǎng)功能,實(shí)驗(yàn)結(jié)果與仿真結(jié)果相同。證明了該方法的優(yōu)點(diǎn):加快開發(fā)速度,容易實(shí)現(xiàn)代碼的驗(yàn)證、測試等工作,系統(tǒng)代碼升級(jí)方便等。對(duì)電力電子裝置的研發(fā)具有重要意義。
[1] 鄧焱弢,駱光照,陳哲,等.基于模型設(shè)計(jì)的處理器在回路聯(lián)合仿真系統(tǒng)[J].測控技術(shù),2011,30(3):45-48.
[2] 馮璨,李剛,趙剛.基于模型設(shè)計(jì)思想的功率測量技術(shù)研究[J].電測與儀表,2014,51(19):1-7.
[3] 楊露露,時(shí)斌,聞楓.基于改進(jìn)的無差拍控制方法在單相并網(wǎng)逆變器中的應(yīng)用[J].水電能源科學(xué),2014,32(5):172-175.
[4] 黃文俊,余曉建,沈永林.單相并網(wǎng)逆變器的Deadbeat控制[J].電力電子技術(shù),2007,41(2):10-12.
[5] 單竹杰,林明耀,顧娟,等.單相光伏并網(wǎng)逆變器的無差拍控制研究[J].電工電氣,2009,28(1):5-7.