湯 超,趙雪松
(安徽工程大學機械與汽車工程學院,安徽 蕪湖 241000)
泵是一種輸送液體或增加液體能量的流體機械,在各類泵中以單級單吸離心泵的應用最廣,砂漿泵屬于其中一種.據(jù)調查,泵的耗電量約占世界總發(fā)電量的1/4,而光單級泵的耗電量就占到總電量的1/12.面對當今工業(yè)生產節(jié)能降耗的要求,在泵行業(yè)開發(fā)出節(jié)能高效的單級單吸離心泵產品顯得尤為重要.
目前人們在泵的可靠性方面研究多一些,對其節(jié)能調速控制的研究還比較少.國內尚未見到將自適應控制運用于砂漿泵的研究,也未對泵的高效轉速控制范圍進行分析,劃分的控制范圍過大,控制效果往往不能令人滿意[1].市場上有一些泵的電機轉速是恒定的,當遇到某些特殊工況時,需手動調節(jié)閥門,這種調節(jié)方式不夠精確且勞動強度大.最主要的缺陷是這樣操作即便達到工況要求,電機仍保持恒定轉速.如果此時泵在低于設計流量的情況下工作,就會加速損耗過流部件并造成很大的能源浪費.
綜上所述,展開對泵的節(jié)能調速控制研究具有明顯的經濟和環(huán)保意義.目前常用的控制方法有PID控制、模糊控制、自適應控制和魯棒控制[2].本文采用傳統(tǒng)PID控制和模型參考自適應模糊控制對砂漿泵系統(tǒng)進行仿真對比研究,比較兩種控制方案的結果,選出實現(xiàn)泵節(jié)能調速的最好方案.
節(jié)能調速控制系統(tǒng)的框圖如圖1所示.
圖1 砂漿泵節(jié)能調速控制系統(tǒng)框圖
如圖1所示,砂漿泵節(jié)能調速控制系統(tǒng)由砂漿泵、模糊控制器、變頻器、轉速傳感器和電機組成,控制量為給定轉速,被控制量為實際轉速.控制過程為:比較輸入給定轉速與實際轉速,用模糊控制器控制差值信號,將所得到的輸出變量去控制變頻器,變頻器控制電機轉速,進而達到砂漿泵的調速控制目的[3].
砂漿泵的數(shù)學模型可簡單等效為純滯后慣性環(huán)節(jié)和帶死區(qū)慣性環(huán)節(jié)的串聯(lián),表示為
式中,K為總增益,τ為純滯后時間常數(shù),T1為慣性時間常數(shù),T2為變頻器時間常數(shù).
參考市場上砂漿泵實際參數(shù),取K=20,T1=10,T2=4,τ=4.代入(1)式可得砂漿泵系統(tǒng)的數(shù)學模型為
根據(jù)本文控制方法要求,需建立一個合適的與實際系統(tǒng)相配套的參考模型,它會影響到后面模糊控制器的性能.在這里,砂漿泵系統(tǒng)采用典型的二階模型.傳遞函數(shù)為
根據(jù)相關參數(shù)選取原則,得到砂漿泵系統(tǒng)參考模型為
PID是當前最常用的一種控制方式,因其控制器結構簡單、控制效果不錯,已有很普遍的應用.本文先用傳統(tǒng)PID控制進行仿真研究.通過分析計算,選出一個控制效果最佳的PID控制器,在SIMULINK中搭建模型,如圖2所示.仿真結果如圖3所示.
由圖3可以看出:系統(tǒng)在第1 s內無響應,運行到20 s后才達到穩(wěn)定,且超調量明顯較大,難以滿足系統(tǒng)節(jié)能調速控制的要求.對于非線性、時變和參數(shù)不定的系統(tǒng),傳統(tǒng)PID無法做到精準控制,有明顯缺陷,因此針對本文砂漿泵調速控制系統(tǒng)其控制效果并不理想.
圖2 傳統(tǒng)PID控制的SIMULINK模型
圖3 傳統(tǒng)PID控制的仿真結果
模糊控制是當前控制領域的一大熱門,它結構簡單,魯棒性好,與自適應控制相結合,對于那些非線性、大時滯的復雜系統(tǒng)有著更好的控制效果.圖4為模型參考自適應模糊控制的工作原理框圖[4],在整個系統(tǒng)中輸入信號r是一單位階躍函數(shù),輸出信號ym為一單調上升曲線趨近于1.當對象的結構或參數(shù)發(fā)生變化時,系統(tǒng)的自適應機構發(fā)生作用,不斷修正兩個輸入的增益信號E和˙E.Yp要能跟上ym,當t→∞時,有
模型參考自適應模糊控制系統(tǒng)原理圖如圖4所示.
圖4 模型參考自適應模糊控制系統(tǒng)原理圖
本文采用控制更精準的電流控制,即輸入輸出的信號都為電流信號.
3.2.1 建立隸屬度函數(shù)
建立e和e˙的隸屬度函數(shù):
(6)式中ym和yp分別表示參考模型和實際系統(tǒng)在單位階躍輸入下的輸出量.
(8)式中K1是輸入變誤差e的量化因子,(9)式中K2是誤差變化˙e的量化因子,(10)式中K3是輸出校正量u的量化因子.設定模糊控制器的輸入輸出論值域為[-6,6].在對稱的隸屬函數(shù)中,以鐘型、梯形和三角形隸屬度函數(shù)的應用最廣,鐘型隸屬度函數(shù)相較另外兩種函數(shù),其曲線上的點平滑且都為非零點,因此本文全部選用鐘型隸屬度函數(shù)[5-6],則E和˙E的隸屬度函數(shù)如圖5所示.
圖5 E和˙E的隸屬度函數(shù)
輸入變量E的語言變量取值:負大(NB),負中(NM),負小(NS),零(ZO),正小(PS),正中(PM),正大(PB).E的隸屬度賦值表如表1所示.
表1 E的賦值表
同理,輸入變量˙E取值也是:負大(NB),負中(NM),負小(NS),零(ZO),正小(PS),正中(PM),正大(PB).˙E隸屬度賦值表如表2所示.
表2 ˙E的賦值表
模糊控制器的輸出模糊論域通過量化因子轉化,取值為[-6,6],即 e=[-emax,emax]=[-6,6],Δe= [- Δemax,Δemax]= [-6,6].由經驗得到E和˙E的基本論域分別為:
對它們量化因子,得到:
同理,根據(jù)經驗取基本論域 u=[-umax,umax]=[-1,1],而模糊論域 U={c}= [-nc,-nc+1,…,-1,0,1,…,nc-1,nc] = [-1,1],這樣就可求出輸出變量U的量化因子K3:
輸出變量U的隸屬度函數(shù)如圖6所示.
圖6 輸出變量U的隸屬度函數(shù)
輸出變量 U取值:負大(NB),負中(NM),負小(NS),零(ZO),正小(PS),正中(PM),正大(PB).其賦值表如表3所示.
表3 U的隸屬度賦值表
3.2.2 建立控制規(guī)則
通過以上的計算分析得出E和˙E的控制規(guī)則,如表4所示.
表4 E和˙E的控制規(guī)則
3.2.3 計算模糊關系
根據(jù)模型參考自適應模糊原理,輸入變量,誤差E和誤差變化˙E要不斷修正,因此它們是平行關系用“and”表示,而輸出變量U用“then”表示,邏輯關系為:If E=NB and˙E=NB then U=PB.
總的模糊關系為
即:
根據(jù)以上分析,在MATLAB中建立起模糊控制器,接下來進行仿真.
按圖4所示系統(tǒng)結構搭建控制系統(tǒng)的模型,如圖7所示.
圖7 自適應模糊控制系統(tǒng)的Simulink模型
將時間設為0~10 s,進行仿真后的結果如圖8所示.橫坐標為系統(tǒng)響應時間;縱坐標為控制量轉速.
圖8 模型參考自適應模糊控制仿真結果
由圖8可以看出:因系統(tǒng)存在純滯后環(huán)節(jié),第1 s內無響應,到第4 s時就達到穩(wěn)定,并且系統(tǒng)超調量較小,對比圖3可以看出模型參考自適應模糊控制比PID控制的效果更理想[7-10].
1)傳統(tǒng)PID控制對于簡單的線性系統(tǒng)控制效果較好,但當運用到砂漿泵這種復雜的非線性系統(tǒng)中時,它易受干擾而產生超調,延長系統(tǒng)達到穩(wěn)定狀態(tài)時間,超調量大,不能滿足砂漿泵節(jié)能調速控制的要求.
2)模型參考自適應模糊控制將泵實際系統(tǒng)與參考模型的輸出信號作比較,通過模糊控制器對誤差劃分論域,最后給出相應的控制量,最終系統(tǒng)在誤差趨近于0時達到穩(wěn)定.采用模型參考自適應模糊控制,砂漿泵調速控制系統(tǒng)穩(wěn)定時間大大縮短,超調量小,可作為砂漿泵節(jié)能調速的控制方案.
[1]劉釗,張珊珊.變量泵控制方式及其應用[J].中國工程機械學報,2004(3):304-307.
[2]安寧,邱瑋煒,戚烜.智能控制綜述[J].技術與市場,2010(5):10-11.
[3]馮浩.變頻調速技術在離心泵節(jié)能中的應用[J].石油化工設計,2009(3):45-47,21.
[4]擺玉龍,楊利君,柴乾隆.基于系統(tǒng)辨識的模型參考自適應控制[J].自動化與儀器儀表,2011(3):23-26.
[5]邢清華,劉付顯.直覺模糊集隸屬度與非隸屬度函數(shù)的確定方法[J].控制與決策,2009(3):393-397.
[6]趙永娟,孫華東.基于Matlab的模糊PID控制器的設計和仿真[J].微計算機信息,2009(1):48-49,40.
[7]郝朝會,孫傳祝,蘇夏侃.自適應模糊PID控制在茶葉殺青機中的應用[J].農機化研究,2013(2):201-204.
[8]劉修廷,潘宏俠.自適應模糊神經控制系統(tǒng)及其MATLAB實現(xiàn)[J].機械工程與自動化,2010(3):162-164.
[9]徐彤,崔高健,曲永利.模糊控制集中供熱控制器的MATLAB 仿真[J].區(qū)域供熱,2013(1):39-42.
[10]Wu C Y,Zhao J,Zheng B C.State tracking of MRAC systems in the presence of controller temporary failure based on a switching method[J].Mathematical Problems in Engineering,2013,213:1-9.