王 蕓
(山西交通職業(yè)技術(shù)學院,山西太原 030006)
軋機作為鋼鐵企業(yè)重要生產(chǎn)設備,其輥縫控制系統(tǒng)控制性能直接影響軋帶品質(zhì)[1]。只有提高輥縫控制系統(tǒng)控制精度和抗干擾性能,才能獲得精度高、符合生產(chǎn)需求的軋帶。
樊立萍等[2]將模糊邏輯控制和PID控制結(jié)合,構(gòu)成模糊自適應PID應用于軋機輥縫控制系統(tǒng),仿真結(jié)果表明該控制方法比傳統(tǒng)的PID控制有更好的控制效果;崔佳梅等[3]通過PID、神經(jīng)網(wǎng)絡、無模型自適應控制3種算法對軋機輥縫控制系統(tǒng)進行仿真比較,結(jié)果表明無模型自適應控制收斂快,穩(wěn)態(tài)誤差小,且能提高對擾動和模型變化的魯棒性;方一鳴等[4]針對軋機輥縫控制系統(tǒng)存在非線性特性、參數(shù)不確定性以及控制輸入前具有不確定系數(shù)的問題,提出了一種自適應反步滑??刂品椒?,仿真結(jié)果表明所設計自適應反步滑??刂破髂軐崿F(xiàn)快速準確跟蹤,并且對參數(shù)變化具有較強魯棒性;王憲等[5]針對軋機輥縫控制系統(tǒng)內(nèi)部參數(shù)不同和外界干擾不一致造成軋輥壓下位置不同步的問題,設計了基于模糊輸出反饋的輥縫伺服同步控制系統(tǒng),仿真結(jié)果表明基于模糊輸出反饋的控制系統(tǒng)不僅具有較短的同步時間并且跟隨誤差更?。焕罱ㄐ鄣萚6]提出一種基于Anti-windup的抗飽和魯棒動態(tài)輸出反饋控制算法,將所提出的算法應用于某650 mm軋機輥縫控制系統(tǒng),其結(jié)果驗證了所提出算法的有效性;付興建等[7]針對軋機輥縫控制系統(tǒng)采用了魯棒H∞控制器對參數(shù)攝動等不確定性進行了仿真,結(jié)果表明所設計的H∞控制器具有良好的跟蹤特性和對參數(shù)攝動的魯棒性;彭珍瑞等[8]在分析數(shù)控機床進給伺服系統(tǒng)數(shù)學模型的基礎上,采用一種基于人工魚群算法優(yōu)化PID控制器參數(shù),并與傳統(tǒng)的Ziegler-Nichols法進行比較.仿真實驗結(jié)果表明:人工魚群PID控制方法的效果明顯優(yōu)于傳統(tǒng)的PID控制,具有良好的動態(tài)和穩(wěn)態(tài)性能;徐秀芬[9]針對電液比例系統(tǒng)PID控制參數(shù)整定問題,提出了基于蟻群算法的PID參數(shù)優(yōu)化方案,并給出了具體的實現(xiàn)步驟,結(jié)果表明:VACA-PID控制器具有良好的靜、動態(tài)性能,完全能達到電液比例控制系統(tǒng)的要求。
以上研究成果可以看出,目前針對軋機輥縫控制系統(tǒng)的研究較多,但主要集中于系統(tǒng)階躍響應性能的研究,研究工況單一,控制器無法滿足軋機復雜工況對系統(tǒng)控制性能的要求。因此,本研究介紹了軋機輥縫控制系統(tǒng),建立了系統(tǒng)開環(huán)傳遞函數(shù),采用人工魚群算法和蟻群算法對PID控制器進行參數(shù)優(yōu)化,基于Simulink比較了優(yōu)化后的系統(tǒng)響應性能和抗干擾性能。
控制軋機工作輥運動的執(zhí)行機構(gòu)為2個液壓缸,如圖1所示,可以控制工作輥在垂直方向運動。液壓缸位移由測厚儀反饋給計算機,其運動由伺服閥閥芯運動進行控制,從而實現(xiàn)對軋機輥縫的精確控制。軋機輥縫控制系統(tǒng)如圖2所示。
分析圖2可以得出,軋機輥縫控制系統(tǒng)為液壓缸位置控制,其數(shù)學模型建立過程目前比較完善,因此本研究將系統(tǒng)關(guān)鍵部分數(shù)學模型列出[10]。
伺服放大器數(shù)學模型為:
I=K1Δu
(1)
式中,Δu—— 輸入電壓,V
I—— 輸出電流,A
K1—— 放大系數(shù),A/V
伺服閥數(shù)學模型為:
(2)
式中,xv—— 閥芯位移,m
Kb—— 閥芯位移與電流增益系數(shù),m/A
圖1 軋機輥縫液壓原理圖
圖2 軋機輥縫控制系統(tǒng)
測厚儀數(shù)學模型為:
Uf=Kfxp
(3)
式中,Kf—— 厚度反饋增益,V/m
Uf—— 反饋電壓,V
xp—— 活塞輸出位移,m
活塞位移xp對閥芯位移xv的傳遞函數(shù)為:
(4)
式中,Kq為伺服閥流量增益,m2/s。
活塞位移xp對負載干擾FL的傳遞函數(shù)為:
(5)
式中,Vt為液壓彈簧剛度最小時的總工作容積,m3。
下面確定系統(tǒng)中主要參數(shù)。液壓缸選用原有液壓缸,半徑為39.4 mm。流量系數(shù)取值范圍0.6~0.62,本研究取0.62。放大器放大系數(shù)及伺服閥固有頻率、節(jié)流口面積梯度、油液密度、厚度反饋增益等可查閱相關(guān)獲得,系統(tǒng)主要參數(shù)見表1。
表1 系統(tǒng)主要參數(shù)
將各參數(shù)代入式(4)和式(5),活塞位移xp對閥芯位移xv的傳遞函數(shù)為:
(6)
活塞位移xp對負載干擾FL的傳遞函數(shù)為:
(7)
軋機輥縫控制系統(tǒng)為典型的液壓缸位置控制,工業(yè)上常采用PID控制器,常規(guī)PID控制器如下式所示:
(8)
本研究選用更適用于計算機控制的位置式PID控制器[11]。假設采樣周期為T,采樣時間為t,則離散PID可表示為:
(9)
將式(9)代入式(8)中即可提到位置式PID控制器的表達式,表達式如式(10)。
式(10)中,β為積分項開關(guān)系數(shù),β取值遵循下式:當ek≤ε時,β=1;當ek>ε時,β=0。其中ε為閾值。ε取值將直接影響系統(tǒng)性能,根據(jù)軋機實際工況,這里選取ε為0.01。
PID控制器中最重要的是其3個參數(shù)的優(yōu)化。下面分別采用蟻群算法和人工魚群算法對PID控制器的3個參數(shù)進行優(yōu)化。
首先采用蟻群算法優(yōu)化PID 3個參數(shù)[12]。蟻群算法可模仿真實的螞蟻行為,螞蟻之間通過一種外激素物質(zhì)進行信息傳遞,螞蟻在運動路徑上留下這種信息素,蟻群集體行為表現(xiàn)出一種信息正反饋現(xiàn)象。蟻群算法優(yōu)化PID結(jié)構(gòu)框圖如圖3所示,其優(yōu)化具體流程如圖4所示。
圖3 蟻群算法優(yōu)化PID參數(shù)結(jié)構(gòu)框圖
圖4 蟻群算法優(yōu)化PID參數(shù)流程圖
圖4中信息素更新變化按式(11)所示:
τ(xn,ynm,t)=ρτ(xn,y,t)+Δτ(xn,ynm,t)
(11)
(12)
Δτk(xn,ynm,t)=Q/Fk
(13)
式中,ρ為揮發(fā)因子,0<ρ<1;Fk為第k只螞蟻的目標函數(shù)值;Q為螞蟻完成一次探索釋放信息素總量,為一個常數(shù)。
根據(jù)蟻群算法優(yōu)化PID參數(shù)流程在MATLAB軟件中進行編程。部分編程程序如圖5所示。
圖5 蟻群算法優(yōu)化PID參數(shù)程序
具體過程為:設置螞蟻數(shù)量為40,揮發(fā)系數(shù)ρ取0.8,信息啟發(fā)因子α取0.4,最大迭代次數(shù)NC取100。限定截割頭位姿調(diào)節(jié)器3個參數(shù)Kp,Ti和Td范圍設置為(0,30)。采用ITAE作為目標函數(shù),如式(14)所示:
(14)
經(jīng)過100代迭代,可獲得螞蟻算法優(yōu)化后的3個參數(shù),如表2所示。
表2 螞蟻算法優(yōu)化參數(shù)結(jié)果
再次采用人工魚群算法優(yōu)化PID 3個參數(shù)[13]。人工魚群算法通過模仿魚群易在富集營養(yǎng)較多的地方進行聚集的特點對魚群的覓食行為進行最優(yōu)求解。人工魚群算法具有全局快速收斂的優(yōu)點,相比蟻群算法具有更強的跳出局部最優(yōu)的能力。如圖6所示為其結(jié)構(gòu)框圖。人工魚群算法優(yōu)化PID 3個參數(shù)流程,如圖7所示。
圖6 人工魚群算法優(yōu)化PID控制器
圖7 人工魚群算法優(yōu)化PID參數(shù)流程圖
根據(jù)人工魚群算法優(yōu)化PID參數(shù)流程在MATLAB軟件中進行編程。部分編程程序如圖8所示。
圖8 人工魚群算法優(yōu)化PID控制器程序
具體過程如下:設置人工魚的初始種群大小為50,步長為0.8,視野visual=40,嘗試次數(shù)try-umber=40,擁擠度因子σ=10,最大迭代次數(shù)為200,搜索范圍設置為(0,30)。采用ITAE的倒數(shù)作為目標函數(shù)[14],如式(15)所示:
(15)
經(jīng)過200代迭代,魚群適應度提高,可獲得人工魚群算法優(yōu)化后的3個參數(shù),如表3所示。
圖9 截割頭位姿控制系統(tǒng)仿真模型
表3 人工魚群算法優(yōu)化參數(shù)結(jié)果
由上述所建立的軋機輥縫控制系統(tǒng)傳遞函數(shù),在Simulink中建立系統(tǒng)仿真模型,如圖9所示。本研究對系統(tǒng)施加階躍信號以及不同頻率隨機信號,比較經(jīng)2種算法優(yōu)化的系統(tǒng)響應性能和抗干擾性能。
首先對模型加入經(jīng)2種算法優(yōu)化的PID 3個參數(shù),對系統(tǒng)施加階躍信號,仿真得到基于2種算法優(yōu)化的系統(tǒng)階躍響應曲線,如圖10所示。
圖10 階躍響應曲線對比
本研究采用超調(diào)量、調(diào)整時間和穩(wěn)態(tài)誤差3個指標對階躍信號響應曲線進行評價,評價結(jié)果對比如表4所示。
從表4可得,人工魚群算法相比于蟻群算法,階躍響應曲線超調(diào)量下降了13.725%,調(diào)整時間縮短了20.287%,穩(wěn)態(tài)誤差降低了29.167%,因此經(jīng)人工魚群算法優(yōu)化的系統(tǒng)響應性能優(yōu)于蟻群算法。
表4 階躍響應結(jié)果對比
下面對系統(tǒng)施加階躍信號的同時,施加干擾力信號,比較2種算法優(yōu)化的系統(tǒng)抗干擾性能。
仿真得到基于2種算法優(yōu)化的系統(tǒng)在干擾力作用下的階躍響應曲線,如圖11所示。
圖11 階躍響應曲線對比(施加干擾力)
表5為用2種算法優(yōu)化的系統(tǒng)在干擾力作用下響應曲線超調(diào)量、調(diào)整時間和穩(wěn)態(tài)誤差的結(jié)果對比。
表5 干擾力作用下仿真結(jié)果對比
對比數(shù)據(jù)發(fā)現(xiàn):在干擾力作用下,采用人工魚群算法優(yōu)化的系統(tǒng),響應曲線超調(diào)量下降了12.576%,調(diào)整時間縮短了14.583%,穩(wěn)態(tài)誤差降低了25%,因此人工魚群算法優(yōu)化的系統(tǒng)抗干擾性能優(yōu)于蟻群算法。
最后采用Random Number模塊對系統(tǒng)施加1,2,4,8 Hz隨機信號。仿真得到系統(tǒng)隨機信號響應曲線如圖12~圖15所示。
圖12 1 Hz隨機信號系統(tǒng)響應曲線對比
圖13 2 Hz隨機信號系統(tǒng)響應曲線對比
圖14 4 Hz隨機信號系統(tǒng)響應曲線對比
圖15 8 Hz隨機信號系統(tǒng)響應曲線對比
不同頻率隨機信號下的系統(tǒng)響應波動范圍對比如表6所示。
表6 不同頻率隨機信號的響應波動范圍對比
從表6可以看出:隨機信號頻率增加,2種算法下的系統(tǒng)響應曲線波動范圍均逐漸減?。徊煌l率下,相比蟻群算法,采用人工魚群算法的隨機信號響應曲線波動范圍較小,因此采用人工魚群算法優(yōu)化的系統(tǒng)對隨機信號的響應性能優(yōu)于蟻群算法。
為提高軋機輥縫控制精度,建立了軋機輥縫控制系統(tǒng)開環(huán)傳遞函數(shù),分別采用蟻群算法和人工魚群算法對PID 3個參數(shù)進行優(yōu)化,并對優(yōu)化后的系統(tǒng)控制性能進行了仿真分析。主要得出以下結(jié)論:
針對軋機輥縫控制系統(tǒng),人工魚群算法相比蟻群算法,在加入干擾力和不加干擾力兩種情況下,系統(tǒng)階躍信號響應超調(diào)量下降了12%以上,調(diào)整時間縮短了14%以上,穩(wěn)態(tài)誤差降低了25%以上;對于不同頻率隨機信號來說,采用人工魚群算法的系統(tǒng)隨機信號響應曲線波動范圍均小于采用蟻群算法的系統(tǒng)隨機信號響應曲線波動范圍。因此經(jīng)人工魚群算法優(yōu)化的軋機輥縫控制系統(tǒng)對階躍信號和隨機信號的響應性能以及對干擾力的抗干擾性能均優(yōu)于蟻群算法優(yōu)化的系統(tǒng)。