馬翠紅 王俊琦 孟凡偉
摘 要: 針對(duì)傳統(tǒng)PID控制算法不能滿足純滯后、大慣性的蓄熱式加熱爐溫度控制系統(tǒng)工藝需求的問(wèn)題,提出基于動(dòng)態(tài)矩陣控制的爐溫控制算法。首先介紹動(dòng)態(tài)矩陣控制的原理,然后將其運(yùn)用在蓄熱式加熱爐的溫度控制系統(tǒng)中,通過(guò)Matlab軟件對(duì)爐膛溫度模型進(jìn)行仿真,最后將本算法與傳統(tǒng)PID控制算法進(jìn)行比較分析。仿真結(jié)果表明,應(yīng)用動(dòng)態(tài)矩陣控制算法對(duì)蓄熱式加熱爐溫度進(jìn)行控制時(shí)效果更佳,動(dòng)態(tài)性能更優(yōu),是一種擁有廣泛應(yīng)用前景的控制算法。
關(guān)鍵詞: 動(dòng)態(tài)矩陣控制; PID控制; 控制算法; 蓄熱式加熱爐; 溫度控制系統(tǒng); 模型仿真
中圖分類號(hào): TN876?34; TP391.9 文獻(xiàn)標(biāo)識(shí)碼: A 文章編號(hào): 1004?373X(2018)19?0104?04
Abstract: Since the traditional PID control algorithm can′t satisfy the technical requirements of regenerative heating furnace temperature control system for pure lag and big inertia, a furnace temperature control algorithm based on dynamic matrix control is proposed. The control principle of dynamic matrix is proposed, and applied to the temperature control system of regenerative heating furnace. The Matlab software is used to simulate the furnace box temperature model. The proposed algorithm is compared with the traditional PID control algorithm. The simulation results show that the dynamic matrix control algorithm has better control effect and higher dynamic performance for regenerative heating furnace temperature control, and is a perfect control algorithm with widespread application prospect.
Keywords: dynamic matrix control; PID control; control algorithm; regenerative heating furnace; temperature control system; model simulation
蓄熱式加熱爐是基于蓄熱式燃燒技術(shù)的大型熱工設(shè)備,目前我國(guó)普遍使用的爐型為三段式加熱爐,其結(jié)構(gòu)主要包含預(yù)熱段、加熱段以及均熱段[1]。預(yù)熱段在進(jìn)料端,此處溫度相對(duì)較低,預(yù)熱方式主要是通過(guò)煙氣余熱進(jìn)行熱傳遞,提高設(shè)備對(duì)熱能的利用率。加熱段一般位于爐體的中間,其主要任務(wù)是為爐子提供足夠的熱量對(duì)爐內(nèi)的產(chǎn)品進(jìn)行加熱,加熱段爐體溫度相對(duì)較高,能夠使產(chǎn)品迅速達(dá)到所需的溫度。均熱段的位置處于爐子的出料端,在此段爐內(nèi)溫度與產(chǎn)品加熱的預(yù)定溫度相差較小,均熱段能夠確保鋼坯出爐時(shí)各斷面溫度盡可能分布均勻。蓄熱式加熱爐設(shè)置了若干組成對(duì)排列的蓄熱室,工作時(shí)可預(yù)熱待燃燒的空氣或煤氣,具有較高的余熱回收利用率及含氮污染物排放低等優(yōu)點(diǎn),目前普遍應(yīng)用于冶金工業(yè)的加工生產(chǎn)工序中[2]。
隨著冶金鋼鐵工業(yè)的發(fā)展,國(guó)內(nèi)外企業(yè)對(duì)加熱爐溫度控制系統(tǒng)的控制品質(zhì)的要求越來(lái)越高,蓄熱式加熱爐溫度控制系統(tǒng)具有大滯后、大慣性、干擾因素多等特點(diǎn),運(yùn)用傳統(tǒng)PID控制算法很難達(dá)到較為理想的效果,因此需要提出一種新型的控制算法來(lái)滿足加熱工藝的需求[3?4]。動(dòng)態(tài)矩陣控制是21世紀(jì)以來(lái)在冶金鋼鐵領(lǐng)域受到越來(lái)越多應(yīng)用的一種增量式預(yù)測(cè)控制算法。本文分析動(dòng)態(tài)矩陣控制原理,應(yīng)用動(dòng)態(tài)矩陣控制算法對(duì)蓄熱式加熱爐爐溫控制系統(tǒng)進(jìn)行調(diào)控,借助Matlab仿真軟件分別對(duì)動(dòng)態(tài)矩陣控制算法和傳統(tǒng)PID控制算法控制下的爐溫模型進(jìn)行仿真,并對(duì)二者的性能指標(biāo)進(jìn)行對(duì)比分析。
動(dòng)態(tài)矩陣控制算法與其他預(yù)測(cè)控制算法一樣,是一種基于模型的控制。它采用的模型是被控系統(tǒng)的階躍響應(yīng)這一非參數(shù)模型,該算法對(duì)系統(tǒng)模型精度要求低,動(dòng)態(tài)性能好,能夠適用于大滯后、大慣性系統(tǒng),屬于最優(yōu)控制算法[5]。動(dòng)態(tài)矩陣控制算法利用過(guò)程模型預(yù)測(cè)系統(tǒng)在人為控制下將要發(fā)生的動(dòng)態(tài)行為,并且根據(jù)已知的約束條件及要求通過(guò)滾動(dòng)優(yōu)化和反饋校正優(yōu)化當(dāng)前與未來(lái)時(shí)刻的控制量,使系統(tǒng)輸出曲線盡可能的與設(shè)定的參考軌跡相符。因此,動(dòng)態(tài)矩陣控制一般分為預(yù)測(cè)模型、滾動(dòng)優(yōu)化以及反饋校正三個(gè)主要部分[6?7],其工作原理如圖1所示。
1.1 預(yù)測(cè)模型
動(dòng)態(tài)矩陣控制采用被控系統(tǒng)的階躍響應(yīng)模型,給定階躍響應(yīng)[a(t)]后,根據(jù)周期[T]得到各采樣點(diǎn)的值[ai=a(iT),i=1,2,…,N]。[N]為模型時(shí)域長(zhǎng)度,它的選取能反映被控對(duì)象的全部動(dòng)態(tài)特性,應(yīng)保證[aN]約等于穩(wěn)態(tài)值,即[aN≈a(∞)]。根據(jù)線性系統(tǒng)具有比例性和疊加性,可得預(yù)測(cè)模型方程[8?9]:
式中:[A]為動(dòng)態(tài)系數(shù)矩陣;[P]為預(yù)測(cè)時(shí)域長(zhǎng)度,[M]為控制時(shí)域長(zhǎng)度,在進(jìn)行參數(shù)選擇時(shí),一般情況下遵循[M≤P≤N]的原則。
1.2 滾動(dòng)優(yōu)化
動(dòng)態(tài)矩陣控制算法與其他控制算法一樣,是一種按一定優(yōu)化規(guī)則確定控制量的預(yù)測(cè)控制算法。滾動(dòng)優(yōu)化一般可分為調(diào)節(jié)與跟蹤兩個(gè)大問(wèn)題。調(diào)節(jié)問(wèn)題,首先確定從某一時(shí)刻開(kāi)始的未來(lái)[M]個(gè)控制增量,然后在控制增量的作用下,盡最大可能地使未來(lái)[P]個(gè)時(shí)刻的輸出預(yù)測(cè)值近似于設(shè)定的期望值。跟蹤問(wèn)題,確定從某時(shí)刻開(kāi)始的[M]個(gè)控制增量,盡可能地使輸出預(yù)測(cè)值能夠跟蹤期望軌跡。為了使控制達(dá)到更加滿意的效果,將輸出量沿著預(yù)定曲線達(dá)到設(shè)定值,該曲線即為參考軌跡:
滾動(dòng)優(yōu)化并不是一次離線完成的,而是多次在線優(yōu)化的結(jié)果。事實(shí)上,動(dòng)態(tài)矩陣控制算法在每一時(shí)刻都會(huì)按相應(yīng)的優(yōu)化算法計(jì)算出從這一時(shí)刻起的控制增量向量。從原理上講,每次給出的增量向量均不相同,對(duì)每個(gè)向量分量而言,都是滾動(dòng)優(yōu)化計(jì)算的結(jié)果。
1.3 反饋校正
在建立控制系統(tǒng)的預(yù)測(cè)模型時(shí),通常采用的模型參數(shù)[ai]只是控制系統(tǒng)的粗略模型。由于實(shí)際工作時(shí)系統(tǒng)模型參數(shù)變化、模型失配、外界干擾等影響,使根據(jù)預(yù)測(cè)模型得出的預(yù)測(cè)值會(huì)偏離實(shí)際輸出值,如果不加校正,這種偏差會(huì)越來(lái)越大,必要的作法是引入動(dòng)態(tài)校正環(huán)節(jié),利用差值實(shí)時(shí)對(duì)預(yù)測(cè)值進(jìn)行校正[10],即:
蓄熱式加熱爐溫度控制系統(tǒng)是一個(gè)復(fù)雜的大滯后系統(tǒng),控制設(shè)計(jì)時(shí)一般采用一階(或者二階)延遲系統(tǒng)近似表示蓄熱式加熱爐溫度控制系統(tǒng)模型。
2.1 動(dòng)態(tài)矩陣控制算法仿真
采用動(dòng)態(tài)矩陣控制算法對(duì)被控對(duì)象[G(s)]在階躍響應(yīng)下進(jìn)行多次仿真,取得最優(yōu)控制效果時(shí)的輸出曲線和控制作用曲線如圖2所示,此時(shí)的控制參數(shù)為采樣周期[Ts=10],模型時(shí)域長(zhǎng)度[N=100],預(yù)測(cè)時(shí)域長(zhǎng)度[P=25],控制時(shí)域長(zhǎng)度[M=3]。
2.2 PID控制算法仿真
采用PID控制算法對(duì)被控對(duì)象[G(s)]在階躍響應(yīng)下進(jìn)行仿真,選用目前較為常用的齊格勒尼科爾斯(Ziegler?Nichols)整定方法對(duì)PID參數(shù)進(jìn)行整定,這種方法是利用Simulink軟件仿真等方法得到單位階躍響應(yīng)曲線,觀察曲線形狀為S形,過(guò)轉(zhuǎn)折點(diǎn)作切線,經(jīng)計(jì)算得到切線方程為:
2.3 仿真結(jié)果對(duì)比分析
將上述方法得到的動(dòng)態(tài)矩陣控制與PID控制的單位階躍響應(yīng)曲線仿真結(jié)果進(jìn)行對(duì)比,如圖4所示。通過(guò)對(duì)比分析,得出在PID控制算法下系統(tǒng)達(dá)到穩(wěn)定的時(shí)間較長(zhǎng),且超調(diào)量相對(duì)較大,控制效果并不理想,而在動(dòng)態(tài)矩陣控制算法下系統(tǒng)達(dá)到穩(wěn)定的時(shí)間較短,超調(diào)量相對(duì)較小,系統(tǒng)震蕩較小,控制效果更優(yōu)。
由上述兩種控制算法下蓄熱式加熱爐溫度控制系統(tǒng)仿真結(jié)果對(duì)比得知,目前PID控制算法相對(duì)簡(jiǎn)單且技術(shù)成熟,但在處理諸如蓄熱式加熱爐溫度控制系統(tǒng)之類的純滯后、大慣性系統(tǒng)時(shí)整定參數(shù)相對(duì)繁瑣,常常出現(xiàn)小幅震蕩以及相對(duì)較大的超調(diào)量,使溫度控制系統(tǒng)需要較長(zhǎng)的時(shí)間才能重新達(dá)到穩(wěn)定狀態(tài),很難得到理想的控制效果。動(dòng)態(tài)矩陣控制算法超調(diào)量小,系統(tǒng)達(dá)到穩(wěn)定的時(shí)間短、震蕩較小,具有較強(qiáng)的穩(wěn)定性和魯棒性,控制效果更加理想。與此同時(shí),PID控制在生產(chǎn)過(guò)程中得到的傳遞函數(shù)模型很難與實(shí)際模型完全吻合,而動(dòng)態(tài)矩陣控制具有對(duì)系統(tǒng)模型精度要求低、控制精度高,無(wú)模型誤差等優(yōu)點(diǎn),能夠克服蓄熱式加熱爐在加熱工作時(shí)的不確定性、滯后性與非線性,并且能夠有效地處理過(guò)程被控量和操縱變量中的各種元素。
動(dòng)態(tài)矩陣控制是一種基于計(jì)算機(jī)控制技術(shù),具有預(yù)測(cè)功能的增量式算法,該算法能夠?qū)⒛P皖A(yù)測(cè)出的信息進(jìn)行綜合分析與運(yùn)用。動(dòng)態(tài)矩陣控制還利用反饋校正和滾動(dòng)優(yōu)化提高算法的穩(wěn)、準(zhǔn)、快特性。隨著計(jì)算機(jī)控制技術(shù)的發(fā)展和人們對(duì)矩陣?yán)碚摰纳钊胙芯?,更加?yōu)化的動(dòng)態(tài)矩陣控制算法將會(huì)在蓄熱式加熱爐溫度控制系統(tǒng)這種純滯后的工業(yè)控制系統(tǒng)中得到廣泛應(yīng)用。
注:本文通訊作者為王俊琦。
參考文獻(xiàn)
[1] 陳迪安,李春生,苗為人,等.大型在線液態(tài)出渣蓄熱式高溫硅鋼加熱爐設(shè)計(jì)研究和應(yīng)用[J].礦冶,2013,22(z1):157?160.
CHEN Dian, LI Chunsheng, MIAO Weiren, et al. Design and application of large online liquid slag storage high temperature silicon steel heating furnace [J]. Mining and metallurgy, 2013, 22(S1): 157?160.
[2] 布煥存.蓄熱式高溫空氣燃燒技術(shù)的應(yīng)用[J].鋼鐵研究學(xué)報(bào),2005,17(5):1?6.
BU Huancun. Application of regenerative high temperature air combustion technology [J]. Journal of steel research, 2005, 17(5): 1?6.
[3] 劉承林.基于蓄熱式燃燒技術(shù)的加熱爐控制系統(tǒng)[D].天津:天津理工大學(xué),2014.
LIU Chenglin. Heating furnace control system based on regene?rative combustion technology [D]. Tianjin: Tianjin University of Technology, 2014.
[4] 石學(xué)勇.蓄熱式加熱爐的模糊控制器[J].冶金能源,2011,30(5):61?62.
SHI Xueyong. Fuzzy controller of regenerative reheating furnace [J]. Metallurgical energy, 2011, 30(5): 61?62.
[5] 董洪霞.動(dòng)態(tài)矩陣控制算法的改進(jìn)及其仿真研究[D].鞍山:遼寧科技大學(xué),2012.
DONG Hongxia. Improvement of dynamic matrix control algorithm and its simulation research [D]. Anshan: Liaoning University of Science and Technology, 2012.
[6] 高強(qiáng),李航.基于單變量與多變量系統(tǒng)的模型預(yù)測(cè)控制研究[J].計(jì)算機(jī)工程與設(shè)計(jì),2013,34(9):3266?3272.
GAO Qiang, LI Hang. Research on model predictive control based on single?variable and multivariable systems [J]. Computer engineering and design, 2013, 34(9): 3266?3272.
[7] 甄子洋.預(yù)見(jiàn)控制理論及應(yīng)用研究進(jìn)展[J].自動(dòng)化學(xué)報(bào),2016,42(2):172?188.
ZHEN Ziyang. Advances in predictive control theory and application [J]. Acta automatica Sinica, 2016, 42(2): 172?188.
[8] 楊麗華,趙文杰.基于Matlab動(dòng)態(tài)矩陣控制算法設(shè)計(jì)仿真[J].儀器儀表用戶,2012,19(4):64?66.
YANG Lihua, ZHAO Wenjie. Matlab dynamic matrix control algorithm design simulation [J]. Instrument user, 2012, 19(4): 64?66.
[9] 侯立剛,徐利軍,蘇成利,等.基于緊格式線性化的無(wú)模型動(dòng)態(tài)矩陣控制[J].控制工程,2013,20(2):246?249.
HOU Ligang, XU Lijun, SU Chengli, et al. Model?free dynamic matrix control based on compact format linearization [J]. Control engineering, 2013, 20(2): 246?249.
[10] 孫建平,齊園園.動(dòng)態(tài)矩陣控制在電加熱爐溫度控制中的應(yīng)用[J].計(jì)算機(jī)仿真,2013,30(6):386?388.
SUN Jianping, QI Yuanyuan. Application of dynamic matrix control in temperature control of electric heating furnace [J]. Computer simulation, 2013, 30(6): 386?388.
[11] 趙春鋒,孫興朋,王東興.基于Matlab/Simulink下加熱爐傳遞函數(shù)的PID仿真試驗(yàn)[J].節(jié)能,2010,29(7):20?22.
ZHAO Chunfeng, SUN Xingpeng, WANG Dongxing. PID si?mulation test of transfer function of heating furnace based on Matlab/Simulink [J]. Energy saving, 2010, 29(7): 20?22.
[12] 王正林,王勝開(kāi),陳國(guó)順,等.Matlab/Simulink與控制系統(tǒng)仿真[M].2版.北京:電子工業(yè)出版社,2008:223?228.
WANG Zhenglin, WANG Shengkai, CHEN Guoshun, et al. Matlab/Simulink and control system simulation [M]. 2nd ed. Beijing: Publishing House of Electronics Industry, 2008: 223?228.