成旭華,姚建勇,樂貴高

(南京理工大學(xué)機(jī)械工程學(xué)院,290014,南京)

電液伺服系統(tǒng)具有動(dòng)態(tài)響應(yīng)快、體積小、輸出力/轉(zhuǎn)矩大及控制精確度高等突出優(yōu)點(diǎn),在國防、航空航天、民用等領(lǐng)域得到了廣泛應(yīng)用。電液伺服系統(tǒng)是一個(gè)典型的非線性系統(tǒng),包括各種非線性特性和建模不確定性,如摩擦非線性、電液伺服閥流量壓力非線性等非線性特性,電液伺服閥流量增益、泄漏系數(shù)、負(fù)載質(zhì)量、執(zhí)行機(jī)構(gòu)的黏性摩擦系數(shù)、液壓油彈性模量等未知參數(shù)以及外部干擾、未建模動(dòng)態(tài)、死區(qū)等建模不確定性[1-2],這些非線性限制了電液伺服系統(tǒng)高性能控制器的發(fā)展。為了提高非線性系統(tǒng)的控制性能,越來越多的研究人員設(shè)計(jì)了各種非線性控制器[3-6],滑模[7]、自適應(yīng)[8]、自適應(yīng)魯棒[9-10]等方法相繼運(yùn)用到液壓非線性運(yùn)動(dòng)控制中。

液壓系統(tǒng)中部分參數(shù)的不確定性比較大,如液壓油彈性模量可從100 MPa到1 400 MPa,泄漏系數(shù)可從1×10-13m5/(N·s)到1×10-10m5/(N·s)等。傳統(tǒng)自適應(yīng)控制對(duì)參數(shù)自適應(yīng)初值比較敏感,系統(tǒng)存在強(qiáng)參數(shù)不確定性時(shí),初值的選取將影響系統(tǒng)的瞬態(tài)響應(yīng)性能,甚至閉環(huán)穩(wěn)定性。針對(duì)此類問題,多模自適應(yīng)控制被提出,且受到廣泛關(guān)注。文獻(xiàn)[11]針對(duì)線性系統(tǒng)建立含有多個(gè)固定模型和兩個(gè)自適應(yīng)模型,提出一種有停留時(shí)間的切換策略,并分析了系統(tǒng)的穩(wěn)定性。李曉理等設(shè)計(jì)了一種基于有界擾動(dòng)分區(qū)的多模型自適應(yīng)控制,解決了含有界擾動(dòng)的離散時(shí)間系統(tǒng)輸出跟蹤問題,提高了系統(tǒng)控制性能[12]。邵彭真等設(shè)計(jì)了一種多模型切換自適應(yīng)控制器,解決了軋機(jī)電液伺服系統(tǒng)在空載和負(fù)載兩種不同工況的切換控制。針對(duì)線性系統(tǒng),Han等提出一種二階多模型自適應(yīng)控制[14],運(yùn)用凸集的概念,減小了辨識(shí)模型的個(gè)數(shù),且具有一定的魯棒性能[15]。隨后,陳杰等提出非線性系統(tǒng)的二階多模型自適應(yīng)控制[16-17],但沒有考慮存在外部干擾等不確定性的情況。自適應(yīng)控制在存在外部干擾、未建模動(dòng)態(tài)等不確定性時(shí),系統(tǒng)跟蹤性能將會(huì)受到影響。

本文針對(duì)電液伺服系統(tǒng)的特點(diǎn),建立了系統(tǒng)的數(shù)學(xué)模型以及多個(gè)辨識(shí)模型,在辨識(shí)模型中設(shè)計(jì)非線性魯棒項(xiàng),以減小外部干擾等不確定性的影響,并基于反步法和李雅普諾夫分析法設(shè)計(jì)了一種考慮強(qiáng)參數(shù)不確定性和不確定非線性的多模型魯棒自適應(yīng)控制器(MMRAC),使系統(tǒng)得到了漸進(jìn)跟蹤的性能,提高了系統(tǒng)的瞬態(tài)響應(yīng)性能。實(shí)驗(yàn)結(jié)果證明了所設(shè)計(jì)控制器的有效性。

1 系統(tǒng)建模與問題描述

本文所研究的電液伺服系統(tǒng)如圖1所示,圖中右側(cè)的負(fù)載被伺服液壓缸驅(qū)動(dòng)。控制器設(shè)計(jì)的目標(biāo)是在電液伺服系統(tǒng)的輸出下使得負(fù)載位移盡可能地跟蹤任意指定的位移軌跡。電液伺服系統(tǒng)的動(dòng)力學(xué)方程可以表示為

(1)

圖1 電液伺服系統(tǒng)示意圖

忽略伺服液壓缸的外泄漏,則伺服液壓缸的壓力動(dòng)態(tài)方程為[1]

(2)

式中:Vt為伺服液壓缸的總?cè)莘e;βe為有效油液彈性模量;Ct為內(nèi)泄漏系數(shù);QL為負(fù)載流量;q(t)為建模誤差。

采用高響應(yīng)的電液伺服閥,則閥芯位移xv與控制輸入u近似為比例環(huán)節(jié),即xv=kiu,因此得到電液伺服閥流量方程

(3)

(4)

(5)

對(duì)于確定的液壓系統(tǒng),有以下合理的假設(shè)。

假設(shè)1系統(tǒng)的指令信號(hào)yd(t)=x1d(t)∈C3且有界;實(shí)際液壓系統(tǒng)在一般工況下滿足0

假設(shè)2系統(tǒng)參數(shù)不確定性θ的范圍已知,即

θp∈Ωpθ{θp:θpmin≤θp≤θpmax}

(6)

式中:θpmin=[θp1min,θp2min,θp3min,θp4min]T、θpmax=[θp1max,θp2max,θp3max,θp4max]T是向量θp已知的上下界,θp3min>0。

(7)

在控制器設(shè)計(jì)分析過程中,將用到以下引理。

引理1?x∈R,a>0,有

0≤xtanh(x/a)≤|x|

(8)

引理2?x≥0,y>0或x>0,y≥0有

x/x+y≤1

(9)

2 控制器的設(shè)計(jì)

2.1 辨識(shí)模型及參數(shù)自適應(yīng)的設(shè)計(jì)

對(duì)系統(tǒng)進(jìn)行在線辨識(shí),根據(jù)系統(tǒng)狀態(tài)方程(5)設(shè)計(jì)多個(gè)辨識(shí)模型如下

(10)

如上所示,各辨識(shí)模型的形式是一樣的,只是參數(shù)估計(jì)的初值不同。

(11)

式中:j=1,2,3,4表示為第幾個(gè)參數(shù);τij為參數(shù)自適應(yīng)函數(shù),在后續(xù)控制器設(shè)計(jì)中決定其具體形式。

給定如下參數(shù)自適應(yīng)率

(12)

式中Γ>0為正定對(duì)角矩陣。

對(duì)于任意的參數(shù)自適應(yīng)函數(shù)τi,不連續(xù)映射式(12)滿足以下性質(zhì)

P1θi∈Ωθi{θi:θpmin≤θi≤θpmax}

(13)

(14)

2.2 控制器的設(shè)計(jì)

基于各個(gè)辨識(shí)模型,用反步設(shè)計(jì)方法設(shè)計(jì)相應(yīng)的控制器ui。

定義emi1=xi1-x1d,則

(15)

令αi1為xi2的虛擬控制輸入,定義αi1與xi2的誤差為emi2=xi2-αi1,則

(16)

設(shè)計(jì)虛擬控制輸入

(17)

式中k1>0為可調(diào)增益。

將式(17)代入式(16)得

(18)

令αi2為xi3的虛擬控制輸入,定義αi2與xi3的誤差為emi3=xi3-αi2,則

(19)

設(shè)計(jì)虛擬控制輸入

(20)

式中k2>0為可調(diào)增益。

將式(20)代入式(19)得

(21)

對(duì)emi3求導(dǎo),得到其動(dòng)態(tài)方程如下

(22)

基于式(22)設(shè)計(jì)控制器

ui=uia+uis

(23)

式中k3>0為可調(diào)增益。

將控制器式(23)代入式(22)得

(24)

2.3 切換策略的設(shè)計(jì)

對(duì)系統(tǒng)進(jìn)行在線辨識(shí)并基于辨識(shí)模型設(shè)計(jì)控制器,其動(dòng)機(jī)是通過辨識(shí)模型的穩(wěn)定使得被控對(duì)象穩(wěn)定。因此,如果使用多個(gè)辨識(shí)模型,具有最佳性能的模型也將產(chǎn)生最佳的控制性能。為了提高系統(tǒng)控制性能,需要在任何時(shí)刻有可靠的辨識(shí)性能指標(biāo)。本文采用的切換方案是對(duì)基于系統(tǒng)辨識(shí)誤差的指標(biāo)函數(shù)進(jìn)行監(jiān)視,每次切換時(shí)選擇對(duì)應(yīng)于指標(biāo)函數(shù)值最小的辨識(shí)模型的控制器作為系統(tǒng)的當(dāng)前控制輸入。

本文選擇的性能指標(biāo)函數(shù)為

(25)

式中:α>0、β>0為對(duì)應(yīng)項(xiàng)的權(quán)重,分別決定實(shí)時(shí)和歷史辨識(shí)誤差在性能指標(biāo)中的權(quán)重。

切換函數(shù)

S(t)=sub(Ji(t)=Jmin(t),i=1,2,…,N)

u=us(t)

(26)

式中:sub(·)為取下標(biāo)函數(shù);Jmin(t)為t時(shí)刻各個(gè)辨識(shí)模型的性能指標(biāo)函數(shù)中的最小值。將基于性能指標(biāo)函數(shù)最小的辨識(shí)模型設(shè)計(jì)的控制器切換為當(dāng)前控制器,即u=us(t),且切換間隙Tmin>ε,ε為任意小的數(shù),以防止無限切換。

2.4 控制性能及穩(wěn)定性分析

定義系統(tǒng)真實(shí)位移和系統(tǒng)給定指令的誤差z=xp1-x1d。

定理對(duì)于電液伺服系統(tǒng)式(5),取參數(shù)自適應(yīng)函數(shù)

τi=φei3

(27)

式中:φ=[-xp2,-xp3,-gu,-1]T。

取M的自適應(yīng)率

(28)

式中γ>0為參數(shù)自適應(yīng)增益。

選取λ1,λ2,λ3,k1,k2,k3使得如下定義的矩陣為正定矩陣

則控制器式(23)及切換函數(shù)式(26)保證閉環(huán)系統(tǒng)中所有信號(hào)均有界,可獲得漸近穩(wěn)定跟蹤,即當(dāng)t→∞時(shí),z→0,并且提高了系統(tǒng)的瞬態(tài)響應(yīng)性能。

證明假設(shè)某一時(shí)刻控制器為ui,取李雅普諾夫函數(shù)為

(29)

對(duì)其求導(dǎo)得

(30)

將參數(shù)自適應(yīng)率式(11)(28)代入式(30)得

(31)

式中Ei=[emi1,emi2,emi3ei1,ei2,ei3]T。

(32)

對(duì)式(32)兩邊進(jìn)行積分,得

(33)

z=xp1-x1d=(xi1-x1d)-

(xi1-xp1)=emi1-ei1

(34)

由上述證明過程可知,每一個(gè)控制器都可以使系統(tǒng)穩(wěn)定。因此,控制器式(23)和切換函數(shù)式(26)可以使得系統(tǒng)任意時(shí)刻都是穩(wěn)定的,即t→∞時(shí),z→0,且把離真實(shí)系統(tǒng)最近的辨識(shí)模型所設(shè)計(jì)的控制器切換為系統(tǒng)當(dāng)前控制器,提高了系統(tǒng)的瞬態(tài)響應(yīng)性能。

3 實(shí)驗(yàn)結(jié)果及分析

建立如圖2所示的閥控電液位置伺服系統(tǒng),其系統(tǒng)參數(shù)與文獻(xiàn)[3]相同。系統(tǒng)的供油壓力Ps=10 MPa,Pr=0.08 MPa,采樣時(shí)間為0.5 ms。

圖2 電液位置伺服系統(tǒng)實(shí)驗(yàn)平臺(tái)

對(duì)比以下5種控制器以驗(yàn)證所提出控制器的有效性。

(1)MMRAC。本文所提出的多模型魯棒自適應(yīng)控制器,控制器參數(shù)k1=1 000,k2=200,k3=50,λ1=1 000,λ2=400,λ3=100,γ=0.05,Γ=diag[20 000,0.2,1×10-5,50]。建立了4個(gè)辨識(shí)模型,θ1(0)=[-1×105,-100,15,0]T,θ2(0)=[-5×105,-500,30,0]T,θ3(0)=[-1×106,-1 000,50,0]T,θ4(0)=[-2×106,-2 000,90,0]T,Mi(0)=10,δ(t)=10 000/(t2+1)。取α=1,β=1,初始控制器選為u2。

(2)魯棒自適應(yīng)控制器(RAC)。傳統(tǒng)的RAC,參數(shù)自適應(yīng)初值θ(0)=[-5×105,-500,30,0]T,其余參數(shù)同MMRAC。

(3)多模型自適應(yīng)控制器(MMAC)。與MMRAC的區(qū)別在于辨識(shí)模型中沒有非線性魯棒項(xiàng),控制器參數(shù)同MMRAC。

(4)PI控制器?？刂破鲄?shù)kp=8 000,ki=2 000。

(5)VFPI控制器。該控制器是速度前饋PI控制器,速度增益kv選為0.028 1 V·s/mm,這種控制結(jié)構(gòu)在工業(yè)中被廣泛使用,其增益與PI控制器相同。

采用最大跟蹤誤差Me、平均跟蹤誤差μ、跟蹤誤差的標(biāo)準(zhǔn)差σ等指標(biāo)[6]來評(píng)估以上5種控制器的性能。

為驗(yàn)證所提出算法的有效性,選取兩種不同工況進(jìn)行實(shí)驗(yàn)驗(yàn)證。

工況1:取系統(tǒng)指令x1d=10arctan(sin(πt))·[1-e-t]/0.785 4 (mm),該指令保證了假設(shè)1的條件。MMRAC作用下系統(tǒng)位置跟蹤、控制器輸入、伺服液壓缸兩腔的壓力值、切換函數(shù)分別如圖3～圖6所示,5種控制器作用下系統(tǒng)位置跟蹤誤差如圖7所示。5種控制器作用下最后兩周期的性能指標(biāo)如表1所示。

圖3 MMRAC作用下系統(tǒng)的位置跟蹤

圖4 MMRAC的控制輸入

圖5 伺服液壓缸兩腔的壓力值

圖6 MMRAC控制器切換圖

由位置跟蹤誤差圖7及表1中后兩周期的性能指標(biāo)可以看出,MMRAC的各項(xiàng)性能指標(biāo)均最優(yōu)。由于RAC參數(shù)自適應(yīng)初值選取得不合適,導(dǎo)致其跟蹤性能較差,而MMRAC采用了多個(gè)辨識(shí)模型,其瞬態(tài)響應(yīng)較RAC有了明顯的提高。與MMAC相比,MMRAC通過在辨識(shí)模型中加入非線性魯棒項(xiàng),抑制了外部干擾、未建模動(dòng)態(tài)等不確定性的影響,提高了系統(tǒng)辨識(shí)的精度。由于MMRAC基于更為準(zhǔn)確的辨識(shí)模型設(shè)計(jì)相應(yīng)控制器,從而提高了系統(tǒng)的控制性能。從圖6可以看出,MMRAC的初始控制輸入為u2,在經(jīng)過一個(gè)采樣周期后,系統(tǒng)的控制輸入切換到u1,而后面S(t)的值不再變化,說明已經(jīng)切換到最佳的控制器。PI控制器對(duì)模型不確定性有一些魯棒性,其跟蹤誤差比其他4種控制器都要差。VFPI控制器中加入了速度前饋,利用速度前饋消除主要的模型動(dòng)態(tài),同時(shí)通過大的反饋增益抑制其他未考慮的因素,但其性能指標(biāo)仍不如MMRAC。

圖7 工況1下各控制器位置跟蹤誤差比較圖

控制器MeμσPI0.265 10.136 60.079 1VFPI0.062 50.021 20.013 3RAC0.195 30.071 60.046 9MMAC0.054 90.018 50.011 6MMRAC0.038 60.017 90.007 6

工況2:為了進(jìn)一步驗(yàn)證算法的有效性,選取系統(tǒng)指令x1d=30arctan(sin(0.4πt))[1-e-t]/0.785 4 (mm)。5種控制器的位置跟蹤誤差如圖8所示,最后兩周期的性能指標(biāo)如表2所示。

同樣,從實(shí)驗(yàn)結(jié)果可以看出,MMRAC的各項(xiàng)性能指標(biāo)均優(yōu)于其他4種控制器,從而進(jìn)一步驗(yàn)證了所提出算法的有效性。

圖8 工況2下各控制器位置跟蹤誤差比較圖

控制器MeμσPI0.230 30.134 10.062 1VFPI0.055 10.016 00.012 2RAC0.190 90.069 60.042 1MMAC0.045 70.008 50.007 1MMRAC0.032 90.007 80.005 7

4 結(jié) 論

本文針對(duì)電液伺服系統(tǒng)存在的強(qiáng)參數(shù)不確定性和不確定非線性,設(shè)計(jì)了一種多模型魯棒自適應(yīng)控制器。通過設(shè)計(jì)多個(gè)含有非線性魯棒項(xiàng)的辨識(shí)模型對(duì)系統(tǒng)進(jìn)行在線辨識(shí),并基于辨識(shí)模型設(shè)計(jì)相應(yīng)的控制器。采用基于辨識(shí)誤差的性能指標(biāo)函數(shù)作為切換依據(jù),選取最佳辨識(shí)模型及控制器,并把該控制器作為當(dāng)前控制器,解決了傳統(tǒng)自適應(yīng)控制對(duì)參數(shù)自適應(yīng)初值敏感的問題,提高了系統(tǒng)的瞬態(tài)響應(yīng)性能。通過辨識(shí)模型中的非線性魯棒項(xiàng),減小了干擾、未建模動(dòng)態(tài)等不確定非線性的影響。運(yùn)用李雅普諾夫穩(wěn)定性理論,證明了控制器的穩(wěn)定性。實(shí)驗(yàn)結(jié)果表明,通過設(shè)計(jì)多個(gè)含有非線性魯棒項(xiàng)的辨識(shí)模型,可以有效抑制建模不確定性的影響,提高系統(tǒng)跟蹤精度,與MMAC相比,穩(wěn)態(tài)精度提高了28%左右,驗(yàn)證了本文所提算法的有效性。