董瑞佳,董嬪,謝強
(1.唐山工業(yè)職業(yè)技術(shù)學(xué)院自動化工程系,河北唐山 063299;2.河南工業(yè)職業(yè)技術(shù)學(xué)院機械工程學(xué)院,河南南陽 473000;3.武漢理工大學(xué)機電工程學(xué)院,湖北武漢 430070)
電液作動器(Electro-Hydraulic Actuator,EHA)集成了液壓伺服系統(tǒng)中的動力元件、控制元件和執(zhí)行元件,因其功率密度大、精度高、體積小等優(yōu)點,被廣泛應(yīng)用于航空、船舶和冶金等領(lǐng)域。但受系統(tǒng)中存在的動靜態(tài)摩擦力特性、機構(gòu)飽和以及環(huán)境變化等非線性因素的影響,電液作動器中的伺服系統(tǒng)表現(xiàn)為較強非線性[1-2]。尤其在高精度的運動控制中,摩擦成為影響高精度伺服機構(gòu)誤差的主要因素之一,主要表現(xiàn)為電液作動器的低速爬行和極限環(huán)震蕩,從而影響伺服系統(tǒng)的動靜態(tài)特性,產(chǎn)生較大的跟蹤誤差[3-4]。
為了減低摩擦對作動器伺服系統(tǒng)帶來的不利影響,目前學(xué)者和工程人員主要采用降低摩擦因數(shù)和摩擦補償兩種方式。通過潤滑和表面處理降低動靜摩擦因數(shù)來降低摩擦對系統(tǒng)影響的改善措施只能起到降低爬行臨界速度的效果,不能完全消除爬行,而且受工藝水平的影響,摩擦因數(shù)只能在一定范圍內(nèi)變動[5-6]。因此從控制的角度進行摩擦補償成為了國內(nèi)外學(xué)者研究重點,并取得了豐厚的研究成果。TRAN等[7]分別采用LuGre摩擦模型和修正的LuGre摩擦模型對液壓缸進行仿真分析,結(jié)果表明修正后的摩擦模型能更準確地預(yù)測液壓缸的行為。PRABEL 和ASCHEMANN[8]提出一種自適應(yīng)反饋控制的摩擦補償策略,通過非線性觀測器對前饋摩擦模型中的參數(shù)進行估計,改善了系統(tǒng)的位置跟蹤性能。國內(nèi)魏瓊等人[9]通過在低速和高速工況下分別進行試驗驗證,并與 PID 控制方法進行比較,結(jié)果表明非線性自適應(yīng)反步摩擦補償控制方法有效改善了起步階段的動態(tài)滯后現(xiàn)象,減少了低速工況下爬行及高速工況下的黏滑振蕩現(xiàn)象,提高了系統(tǒng)的響應(yīng)速度及跟蹤精度。
凌明祥和朱長春[10]基于Stribeck效應(yīng)建立了相對完善的電液伺服振動試驗臺非線性摩擦力理論模型和待辨識參數(shù)的目標函數(shù)。黨選舉等提出了一種基于不完全微分的改進自適應(yīng)反步積分滑??刂?ID-BIABISMC),實現(xiàn)了基于摩擦補償?shù)囊簤核欧到y(tǒng)位置控制,將反演控制、積分滑模面滑模控制與摩擦補償相結(jié)合,解決了液壓系統(tǒng)的非結(jié)構(gòu)不確定性問題[11]。
為解決摩擦力不確定性與非線性影響電液作動器低速工況位置精度不足問題,提出一種前饋補償+ESO(Extended State Observer)的控制策略,基于前饋補償和擴張狀態(tài)觀測器和LuGre摩擦力模型,建立基于LuGre動態(tài)摩擦模型的系統(tǒng)模型和前饋補償+ESO的控制器,結(jié)合仿真平臺驗證所提出摩擦補償策略的合理性。
早在16世紀中期,人們就開始了對摩擦力的研究,隨著科技和測量技術(shù)的提升,不同的摩擦力特性也逐漸地被研究者揭示出來,如黏性摩擦、Stribeck效應(yīng)、預(yù)滑移、摩擦記憶等。雖然至今國內(nèi)外學(xué)者關(guān)于摩擦機制并沒有十分完善的科學(xué)解釋,但國內(nèi)外學(xué)者基于不同的摩擦力特性提出了不同的摩擦模型。按照摩擦模型是否由微分方程描述可大致分為靜態(tài)摩擦模型和動態(tài)摩擦模型。靜態(tài)摩擦模型主要包括庫侖摩擦模型、黏性摩擦模型、指數(shù)摩擦模型等。動態(tài)摩擦模型主要包括Dahl模型、鬃毛模型、LuGre模型等。法國學(xué)者DE-WIT結(jié)合Dahl和鬃毛模型思想提出了LuGre摩擦模型[12-13]。
LuGre摩擦模型是一個較為完善的動態(tài)摩擦模型,用兩個接觸面間彈性鬃毛的平均偏移來表征摩擦力的動態(tài)行為。接觸面間通過彈性鬃毛接觸,受外力作用,預(yù)滑動階段通過鬃毛彈性形變產(chǎn)生摩擦力,形變量超過一定范圍,產(chǎn)生滑動,穩(wěn)態(tài)運動時鬃毛的變形量和滑動速度成正比。其中鬃毛的平均變形量表示為
(1)
由式(1)可知,穩(wěn)態(tài)的鬃毛平均變形量為
(2)
式中:sign(?)為符號函數(shù)。
(3)
(4)
式中:λ為與靜摩擦相關(guān)的無因次量,通常取λ=2;fc為庫侖摩擦力;fs為最大靜摩擦力;vs為Stribeck速度。
則LuGre模型描述的摩擦力為
(5)
式中:δ0為剛度系數(shù);δ1為阻尼系數(shù);δ2為黏滯系數(shù)。
LuGre模型可以真實地描述目前所觀測到的Stribeck效應(yīng)、預(yù)滑移等大多數(shù)靜態(tài)特性以及瞬態(tài)特性,且因其簡潔的表達形式,受到了控制界的高度重視。式(1)、式(2)、式(3)描述了LuGre動態(tài)摩擦模型,其中參數(shù)fc、fs、vs、δ0、δ1、δ2由實驗數(shù)據(jù)辨識得到。
由式(3)可以看出,LuGre模型中存在不連續(xù)函數(shù)sign(·)。不連續(xù)函數(shù)的存在使得系統(tǒng)在零速附近LuGre模型不可微,給后續(xù)控制帶來困難。因此,國內(nèi)外許多學(xué)者通過設(shè)計連續(xù)的擬合函數(shù)描述摩擦力在零速附近的情況。基于雙曲正切函數(shù)將摩擦力模型描述為
(6)
則
(7)
根據(jù)使用要求和場合的不同,電液作動器的執(zhí)行器有對稱液壓缸和非對稱液壓缸,因非對稱缸在結(jié)構(gòu)尺寸方面的獨特優(yōu)勢,被廣泛使用。但由于兩腔的作用面積不同,非對稱液壓缸的正反方向動態(tài)特性不同,且兩腔流量不平衡,很難建立系統(tǒng)精確的傳遞函數(shù),以伸出單向過程的傳遞函數(shù)來近似系統(tǒng)傳遞函數(shù)的方式容易導(dǎo)致系統(tǒng)的建模誤差增加。因此采用狀態(tài)空間方程來描述閥控非對稱缸系統(tǒng),其結(jié)構(gòu)簡圖如圖1所示。
圖1 閥控非對稱缸的系統(tǒng)
液壓缸活塞動力平衡方程為
(8)
式中:x為液壓缸的伸出位移;m為活塞及負載折算的總質(zhì)量;B為黏性摩擦系數(shù);F為負載;K為彈性負載剛度;A1為液壓缸無桿腔的作用面積;A2為液壓缸有桿腔的作用面積;p1為液壓缸無桿腔的壓力;p2為液壓缸有桿腔的壓力。
非對稱缸兩端的流量連續(xù)方程為
(9)
式中:x為液壓缸的伸出位移;V1為進油腔的初始容積;V2為回油腔的初始容積;βe為油液的彈性模量;Q1為流入液壓缸無桿腔的流量;Q2為流出回油腔的流量;Ct為液壓缸泄漏系數(shù)。
Q1和Q2根據(jù)伺服閥負載流量方程表示為
(10)
式中:Kq1、Kq2為伺服閥閥芯左右兩端流量增益;ps為泵源壓力;pr為回油側(cè)壓力;xv為伺服閥閥芯位移。
其中定義函數(shù)s(x)為
(11)
為簡化系統(tǒng)模型,假設(shè)伺服閥為對稱且匹配的,則有Kq1=Kq2=Kq。且有
(12)
式中:Cd為流量系數(shù);w為伺服閥的面積梯度;ρ為液壓油密度。
定義系統(tǒng)狀態(tài)變量為
(13)
將伺服閥簡化為比例環(huán)節(jié)則xv=kiu,且kqi=kikq,則由式(8)(9)(10)可得如下系統(tǒng)狀態(tài)方程
(14)
由于LuGre摩擦模型可以描述大多數(shù)的摩擦行為,因此為提高系統(tǒng)的伺服性能,結(jié)合式(1)(4)(5)(7)描述的LuGre摩擦模型,則閥控非對稱缸的伺服系統(tǒng)方程可寫為
(15)
式中:fm為LuGre摩擦力。
(16)
已知在伺服系統(tǒng)中,隨著壓力、溫度系統(tǒng)內(nèi)部條件以及外部干擾等的變化,液壓油彈性模量βe、負載質(zhì)量m等參數(shù)會隨之變化。為方便控制器設(shè)計,對于所提出的基于LuGre摩擦模型的伺服系統(tǒng)狀態(tài)空間方程,假設(shè)系統(tǒng)參數(shù)A1、A2、Kq、Ct、m、βe、ps、pr和摩擦力參數(shù)δ0、δ1、δ2為已知固定參數(shù)。
伺服系統(tǒng)因具有一定的相位滯后,輸出量的反饋值經(jīng)過一定滯后疊加在控制量上,而摩擦前饋補償控制可以通過建立精確的摩擦模型通過輸入信號或擾動按照已知模型直接反饋到控制量上,可以在被控制量發(fā)生變化前消除擾動的影響,并不受系統(tǒng)滯后的影響,因此前饋控制可以有效提升高精度伺服系統(tǒng)跟蹤性能。根據(jù)式(5)所建立的基于雙曲正切近似的連續(xù)可微摩擦力模型可將系統(tǒng)的前饋控制描述為圖2。
圖2 前饋控制結(jié)構(gòu)框圖
摩擦前饋補償控制雖能提高伺服系統(tǒng)的跟蹤性能,但是需要了解系統(tǒng)的精確模型,實際的動態(tài)摩擦過程受速度和平均鬃毛變形量的影響,因此前饋摩擦補償僅能補償靜態(tài)摩擦力,并且靜態(tài)摩擦力的參數(shù)辨識和擬合均會影響補償精度,導(dǎo)致補償效果下降。
摩擦前饋補償控制雖能提高伺服系統(tǒng)的跟蹤性能,但是需要了解系統(tǒng)的精確模型,對模型的實際參數(shù)依賴性較高,且系統(tǒng)存在許多外部擾動,前饋摩擦很難完全消除這些誤差。實際的動態(tài)摩擦過程受速度和平均鬃毛變形量z的影響,因此本文作者設(shè)計摩擦前饋補償器和擴張狀態(tài)觀測器(ESO)結(jié)合的控制方式,控制原理圖如圖3所示。
圖3 前饋控制+ESO控制結(jié)構(gòu)框圖
考慮系統(tǒng)
(17)
將不確定干擾設(shè)定為狀態(tài)x4=fn(x1,x2,t),則觀測器可設(shè)計為
(18)
誤差的非線性加權(quán)函數(shù)形式表示為
(19)
式中:α3<α2<α1,取α3=0.1、α2=0.5、α1=1.5;β1、β2、β3、β4為控制器可調(diào)增益系數(shù);δ>0。
擴張狀態(tài)觀測器的控制量可設(shè)計為
u=β1fal(e1,α1,δ)+β2fal(e2,α2,δ)+β3fal(e3,α3,δ)+g(z1,z2,z3,t)/b
(20)
為驗證摩擦前饋補償+ESO控制方式的有效性,使用提出的控制方式和系統(tǒng)空間方程搭建了Simulink仿真平臺,對電液伺服系統(tǒng)進行了位置跟蹤仿真,給出的系統(tǒng)仿真參數(shù)如表1所示。
表1 電液促動器伺服系統(tǒng)仿真參數(shù)
使用MATLAB中的Simulink模塊搭建仿真平臺,并將文中所提出的控制方法的控制效果進行對比分析。仿真給定位置指令為xd=0.1sin(0.1t),所設(shè)計控制器增益參數(shù)設(shè)置為β1=2、β2=0.2、β3=0.1、β4=0.05。取系統(tǒng)的初始狀態(tài)為x1(0)=0,x2(0)=0.01,x3(0)=0,x4(0)=0,x5(0)=0。系統(tǒng)位置跟蹤仿真結(jié)果如圖4—圖6所示,跟蹤誤差曲線如圖7所示。
圖4 PID控制位置跟蹤曲線 圖5 摩擦前饋控制位置跟蹤曲線
圖6 前饋控制+ESO控制位置跟蹤曲線 圖7 不同控制方式跟蹤誤差曲線
對比幾種控制效果,由于干擾和系統(tǒng)參數(shù)的變化,PID控制在系統(tǒng)的伸出和縮回過程均存在較大的跟蹤誤差;摩擦前饋補償控制方式下,由于非對稱缸伸出和縮回系統(tǒng)結(jié)構(gòu)參數(shù)的變化,系統(tǒng)縮回過程的誤差明顯大于伸出過程;而所設(shè)計的摩擦前饋補償+ESO控制方式具有更小的跟蹤誤差,系統(tǒng)伸出和縮回過程能夠很好地跟蹤指令信號,跟蹤誤差僅為常規(guī)PID控制的1/4,跟蹤精度達到0.2 mm。
(1)通過連續(xù)的雙曲正切函數(shù)擬合LuGre摩擦模型中的分段函數(shù)g(x),并在常規(guī)閥控非對稱缸的系統(tǒng)狀態(tài)空間方程中引入新的狀態(tài)變量,即LuGre動態(tài)摩擦力中的鬃毛平均變形量z,以此來建立伺服系統(tǒng)精確的狀態(tài)空間方程。
(2)為提高對動態(tài)摩擦力以及外部干擾的補償效果,提高伺服系統(tǒng)跟蹤精度,結(jié)合非線性加權(quán)函數(shù),將系統(tǒng)的未建模部分和外部干擾擴張為系統(tǒng)狀態(tài)空間變量;依據(jù)系統(tǒng)結(jié)構(gòu)和參數(shù),設(shè)計系統(tǒng)的擴張狀態(tài)觀測器,生成系統(tǒng)的控制量u;綜合建立的LuGre摩擦模型,設(shè)計前饋補償+ESO的控制方式。仿真結(jié)果表明:與PID和前饋補償控制相比,系統(tǒng)具有更好的跟蹤精度。