肖國(guó)林,鄧智勇,唐國(guó)元,黃道敏,謝金輝
(1.武漢第二船舶設(shè)計(jì)研究所,湖北武漢430064;2.華科技大學(xué)船舶與海洋工程學(xué)院,湖北武漢430074;3.武漢空軍預(yù)警學(xué)院,湖北武漢430010)
在某大型武器系統(tǒng)旋轉(zhuǎn)裝置的電液位置伺服控制系統(tǒng)設(shè)計(jì)中,根據(jù)設(shè)計(jì)要求,其階躍響應(yīng)調(diào)整時(shí)間須小于2 s,且應(yīng)盡量避免系統(tǒng)出現(xiàn)超調(diào)。該旋轉(zhuǎn)裝置為低增益、大慣量的控制系統(tǒng),由于液壓系統(tǒng)的復(fù)雜性,使位置控制系統(tǒng)模型為一高階傳遞函數(shù)模型,從而對(duì)其性能優(yōu)化帶來了不便。
近年來,模型降階和降階控制器的設(shè)計(jì)一直都是控制理論中的熱門研究領(lǐng)域,并在過去的幾十年里取得了長(zhǎng)足的發(fā)展和廣泛的應(yīng)用[1]。模型降階方法已廣泛應(yīng)用于電力系統(tǒng)[2]、VLSI大規(guī)模集成電路[3]、大型空間柔性結(jié)構(gòu)[4]、飛機(jī)控制系統(tǒng)[5]等許多控制系統(tǒng)中。文獻(xiàn)[6]研究了非最小相位系統(tǒng)傳遞函數(shù)的演化模型降階方法,其分析結(jié)果表明,降階后所得的模型不僅保留了原有模型的基本特性,也使分析大大地得到了簡(jiǎn)化。
本研究針對(duì)某武器系統(tǒng)實(shí)際電液位置伺服系統(tǒng)為一高階系統(tǒng)的特點(diǎn),分析其傳遞函數(shù)的零、極點(diǎn)分布情況,為使用模型降階方法提供依據(jù),并在此基礎(chǔ)上采用降階近似的方法來設(shè)計(jì)和優(yōu)化位置伺服系統(tǒng),以期取得良好效果。
位置伺服系統(tǒng)原理如圖1所示。
圖1 位置伺服系統(tǒng)結(jié)構(gòu)原理圖
用自整角機(jī)作為比較元件,其輸入θi為自整角機(jī)發(fā)送器轉(zhuǎn)子軸的轉(zhuǎn)角,即系統(tǒng)的輸入信號(hào);θL為自整角機(jī)接收器轉(zhuǎn)子軸的轉(zhuǎn)角,即負(fù)載輸出轉(zhuǎn)角,也就是系統(tǒng)的反饋量。兩者差值經(jīng)相敏放大器作用后成為直流信號(hào),并經(jīng)功率放大后用于驅(qū)動(dòng)電液伺服閥,液壓油經(jīng)過高頻伺服閥進(jìn)入馬達(dá),驅(qū)動(dòng)馬達(dá)正、反轉(zhuǎn),從而帶動(dòng)負(fù)載轉(zhuǎn)動(dòng)。筆者用絕對(duì)式編碼器實(shí)時(shí)監(jiān)測(cè)負(fù)載角度位置,并反饋到輸入端。
根據(jù)自整角機(jī)的工作原理,自整角機(jī)在輸入量與輸出量間的誤差角θe不大的條件下,sin(θi-θL)≈θiθL=θe,因此其模型可以寫成:
自整角機(jī)輸出的交流電壓信號(hào)為us,經(jīng)過相敏放大器整流后,交流電壓信號(hào)轉(zhuǎn)換成直流電壓信號(hào)ug。ug的大小比例于us的幅值,ug的極性對(duì)應(yīng)于us的相位。直流放大器輸入為ug,輸出為差動(dòng)電流Δi。由于這些電器元件的動(dòng)態(tài)過程與液壓動(dòng)力元件相比可以看成比例環(huán)節(jié)[7],相敏放大器的增益為:
功率放大器增益為:
當(dāng)電液伺服閥的頻寬與液壓固有頻率相近時(shí),電液伺服閥的傳遞函數(shù)可用2階環(huán)節(jié)來表示[8],即:
式中:QL—電液伺服閥的輸出流量;Ksv—電液伺服閥增益;Gsv(s)—Ksv=1時(shí)電液伺服閥的傳遞函數(shù);ωsv—電液伺服閥的固有頻率;ξsv—電液伺服閥的阻尼比。
回轉(zhuǎn)裝置為慣性、黏性負(fù)載,考慮到泄漏和油的壓縮性,忽略彈性負(fù)載,則以電液伺服閥的輸出流量QL為輸入,液壓馬達(dá)軸角位移θm的輸出值為[9]:
式中:i—齒輪傳動(dòng)比;Dm—馬達(dá)每弧度排量;Kce—總的流量壓力系數(shù);Vt—馬達(dá)兩腔及連接管道總?cè)莘e;βe—液壓油體積彈性模量;TL—作用于馬達(dá)上的干擾力矩。
且:
式中:Jt—馬達(dá)及負(fù)載的轉(zhuǎn)動(dòng)慣量折算到馬達(dá)軸上的總等效馬達(dá)轉(zhuǎn)動(dòng)慣量。
另外:
式中:Bm—馬達(dá)及負(fù)載轉(zhuǎn)動(dòng)時(shí)的總黏性阻尼系數(shù),通常很小。
由式(1~7)可以畫出系統(tǒng)的方框圖如圖2所示。
圖2 位置伺服系統(tǒng)方框圖
依據(jù)圖2寫出系統(tǒng)的開環(huán)傳遞函數(shù)為:
式中:Kv—開環(huán)增益。
且:
可見:系統(tǒng)的傳遞函數(shù)模型為一高階模型,這對(duì)控制系統(tǒng)設(shè)計(jì)帶來了不便。通過分析系統(tǒng)的零、極點(diǎn)分布情況,可以得到系統(tǒng)的主導(dǎo)極點(diǎn),從而為系統(tǒng)的降階近似處理提供依據(jù)。
本研究所設(shè)計(jì)系統(tǒng)的主要參數(shù)如表1所示。
表1 系統(tǒng)參數(shù)
負(fù)載的轉(zhuǎn)動(dòng)慣量折算到馬達(dá)軸上的總等效馬達(dá)轉(zhuǎn)動(dòng)慣量Jt為:
將相應(yīng)參數(shù)值代入式(4~8),可得系統(tǒng)的開環(huán)傳遞函數(shù)為:
系統(tǒng)的零、極點(diǎn)分布圖如圖3所示。
圖3 系統(tǒng)零極點(diǎn)分布圖
由圖3可知:該高階系統(tǒng)具有兩個(gè)實(shí)數(shù)主導(dǎo)極點(diǎn)s1=-0.522,s2=-1.67。且非主導(dǎo)極點(diǎn)s3,4=-83±110j和s5=-216實(shí)部的模比主導(dǎo)極點(diǎn)實(shí)部的模大3倍以上。
經(jīng)對(duì)傳遞函數(shù)進(jìn)行分析可知:系統(tǒng)為穩(wěn)定高階系統(tǒng)。且主導(dǎo)極點(diǎn)附近無(wú)閉環(huán)零點(diǎn),因此該高階系統(tǒng)可近似成如下2階系統(tǒng):
式中:k—增益因子。
增益因子數(shù)值在高階系統(tǒng)變?yōu)?階系統(tǒng)過程中保持不變,因此有s(0)=s'(0),則可求得k=0.871 74。
故2階系統(tǒng)閉環(huán)傳遞函數(shù)為:
通過計(jì)算機(jī)計(jì)算原5階系統(tǒng)和近似3階系統(tǒng)的單位階躍響應(yīng),可知基于一對(duì)實(shí)數(shù)主導(dǎo)極點(diǎn)求取的高階系統(tǒng)單位階躍響應(yīng)與降階近似處理后的過阻尼2階系統(tǒng)的單位階躍響應(yīng)吻合度很高。
由傳遞函數(shù)可知特征方程有兩個(gè)不相等的負(fù)實(shí)根,對(duì)應(yīng)于s平面負(fù)實(shí)軸上的兩個(gè)不等實(shí)極點(diǎn)。因此該2階系統(tǒng)阻尼大于1,為過阻尼系統(tǒng),相應(yīng)的單位階躍響應(yīng)也是非周期地趨于穩(wěn)態(tài)輸出,但響應(yīng)速度較慢。系統(tǒng)需要校正,提高其響應(yīng)速度。
根據(jù)2階系統(tǒng)的閉環(huán)傳遞函數(shù)寫出其開環(huán)傳遞函數(shù)為:
一般情況下,過阻尼系統(tǒng)響應(yīng)緩慢,故通常不采用過阻尼系統(tǒng)[10]。但是,本研究的回轉(zhuǎn)裝置為低增益、大慣量的控制系統(tǒng),按設(shè)計(jì)要求應(yīng)盡量避免出現(xiàn)超調(diào),因此希望對(duì)系統(tǒng)進(jìn)行校正,使之在保持無(wú)超調(diào)的前提下提高其響應(yīng)速度。
本研究對(duì)近似2階系統(tǒng)進(jìn)行校正的思路為:在系統(tǒng)的前向通道中串入校正環(huán)節(jié),串入的校正環(huán)節(jié)為k(1+T1S)/(1+T2S),希望在校正后系統(tǒng)變?yōu)榫哂辛己脛?dòng)態(tài)特性的臨界阻尼二階系統(tǒng),因此可設(shè)置校正環(huán)節(jié)的參數(shù),使校正后系統(tǒng)的閉環(huán)傳遞函數(shù)中有一個(gè)零點(diǎn)與一個(gè)極點(diǎn)恰好相等,相互抵消。在此基礎(chǔ)上提出理想的阻尼比和調(diào)整時(shí)間,由此建立方程,得出校正環(huán)節(jié)的系數(shù)。
在系統(tǒng)中串入校正環(huán)節(jié)k(1+T1S)/(1+T2S),系統(tǒng)的開環(huán)傳遞函數(shù)變?yōu)?
則系統(tǒng)的閉環(huán)傳遞函數(shù)可寫為:
取ξ=1,調(diào)整時(shí)間為ts=1,則ωn=4.75,代入上式計(jì)算,可得到串聯(lián)校正環(huán)節(jié)為:
將校正環(huán)節(jié)加入到高階系統(tǒng)中,可得到校正后的高階系統(tǒng)開環(huán)傳遞函數(shù)為:
校正后閉環(huán)函數(shù)的階躍響應(yīng)和開環(huán)函數(shù)的Bode圖如圖4、圖5所示。
圖4 校正后高階系統(tǒng)的階躍響應(yīng)
由圖4可知,校正后的系統(tǒng)調(diào)整時(shí)間約為1 s,系統(tǒng)無(wú)超調(diào)。系統(tǒng)能在較短時(shí)間內(nèi)達(dá)到穩(wěn)定,動(dòng)態(tài)性能較好。
由圖5可知,校正后的系統(tǒng)穿越頻率為0.319 Hz,相角裕度為77°。截止頻率為3.97 Hz,幅值裕度為29.3 dB,系統(tǒng)具有良好的穩(wěn)定性。仿真所得到的調(diào)整時(shí)間約為1 s。
本研究通過對(duì)某大型旋轉(zhuǎn)裝置原理的分析,建立了系統(tǒng)傳遞函數(shù)模型,對(duì)其電液位置伺服系統(tǒng)進(jìn)行了仿真;并使用降階近似的方法,得出了其近似二階系統(tǒng)模型;提出了系統(tǒng)校正方法,并據(jù)此計(jì)算了校正環(huán)節(jié)具體參數(shù)。仿真結(jié)果表明,所設(shè)計(jì)的系統(tǒng)能較好滿足設(shè)計(jì)要求。
圖5 校正后系統(tǒng)的開環(huán)Bode圖
通過對(duì)上述系統(tǒng)采用降階近似方法進(jìn)行分析和設(shè)計(jì)的過程可知,對(duì)于類似高階位置伺服系統(tǒng),為了簡(jiǎn)化其控制系統(tǒng)設(shè)計(jì),可以分析系統(tǒng)零、極點(diǎn)分布,對(duì)于在其時(shí)間響應(yīng)中暫次要作用的非主導(dǎo)極點(diǎn)進(jìn)行簡(jiǎn)化,從而得到了近似低階系統(tǒng),并在此基礎(chǔ)上進(jìn)行了控制系統(tǒng)的設(shè)計(jì)。上述方法為同類型控制系統(tǒng)的設(shè)計(jì)提供了有益參考。
[1]杜鑫.基于LMI技術(shù)的線性系統(tǒng)模型降階與靜態(tài)輸出反饋控制器設(shè)計(jì)[D].沈陽(yáng):東北大學(xué)信息科學(xué)與工程學(xué)院,2009.
[2]王佩,劉永強(qiáng).無(wú)刷雙饋風(fēng)力發(fā)電機(jī)模型降階研究:第I部分一無(wú)刷雙饋風(fēng)力發(fā)電機(jī)多時(shí)間尺度模型[J].控制理論與應(yīng)用,2008,25(1):135-138.
[3]TAN S,HE L.Advanced Model Order Reduetion Techniques in VLSI Design[M].Cambridge University Press,2007.
[4]FREITAS F D,ROMMES J,MARTINS N.Gramian-based reduetion method applied to large sparse power system deseriptor models[J].IEEE Transactions on Power Systems,2008,23(3):1258-1270.
[5]王欣,史忠科.模型降階H∞方法在陣風(fēng)干擾著陸控制中的應(yīng)用[J].飛行力學(xué),2000,12(4):35-37.
[6]黃樟燦,陳思多.非最小相位系統(tǒng)傳遞函數(shù)的演化模型降階[J].系統(tǒng)工程與電子技術(shù),2000,22(12):60-62.
[7]高翔,朱蕓,梁棟.電液伺服系統(tǒng)的建模與仿真研究[J].海軍工程學(xué)院學(xué)報(bào),1999(4):87-90.
[8]王占林.近代電氣液壓伺服控制[M].北京:北京航空航天大學(xué)出版社,2005.
[9]易孟林,曹樹平,劉銀水.電液控制技術(shù)[M].武漢:華中科技大學(xué)出版社,2010.
[10]楊叔子,楊克沖.機(jī)械工程控制基礎(chǔ)[M].武漢:華中科技大學(xué)出版社,2005.