劉春芳 榮 剛
沈陽工業(yè)大學,沈陽,110870
加工中心雙懸浮系統(tǒng)支持向量機α階逆解耦時變滑??刂?/p>
劉春芳榮剛
沈陽工業(yè)大學,沈陽,110870
龍門數(shù)控加工中心移動機械橫梁可采用雙電磁懸浮系統(tǒng)共同懸浮來消除摩擦的影響,而這兩個電磁懸浮系統(tǒng)存在著一定的耦合關系。當數(shù)控機床工作時,由于受力不平衡或者擾動等因素影響,這種耦合的存在會使兩個電磁懸浮系統(tǒng)的懸浮氣隙受到影響,并降低加工精度。分析得出了橫梁發(fā)生繞質心旋轉和上下平移時兩個電磁懸浮系統(tǒng)的耦合定量關系,為了消除耦合的影響,采用支持向量機逼近雙電磁懸浮系統(tǒng)的α階逆系統(tǒng),將所得到的α階逆系統(tǒng)串聯(lián)在原系統(tǒng)前構成偽線性復合系統(tǒng),從而將原系統(tǒng)解耦成兩個獨立的SISO偽線性系統(tǒng)。針對解耦后的偽線性SISO系統(tǒng)設計了時變滑模變結構控制器,可使系統(tǒng)在任意初始狀態(tài)下系統(tǒng)狀態(tài)變量都能直接到達系統(tǒng)的滑模面上,消除了狀態(tài)變量到達滑模面的過程,以最短的時間實現(xiàn)了滑模變結構控制,實現(xiàn)了對參數(shù)攝動和外部干擾的全局魯棒性。仿真實驗結果表明,該方案不需要被控對象精確的數(shù)學模型即可采用支持向量機逼近被控系統(tǒng)的α階逆系統(tǒng),可有效地對耦合的懸浮系統(tǒng)實現(xiàn)解耦。解耦后的SISO偽線性系統(tǒng)采用時變滑模變結構控制具有響應速度快和魯棒性強等特點。
α階逆系統(tǒng);解耦;支持向量機;時變滑模變結構
隨著科學技術的發(fā)展,人們對零件加工精度的要求也越來越高。傳統(tǒng)的龍門移動數(shù)控機床移動橫梁與導軌完全接觸存在著摩擦,即使加入潤滑劑也不能完全消除摩擦影響。為了消除摩擦的影響,利用兩個電磁懸浮系統(tǒng)將移動橫梁懸浮起來,可提高加工精度。由于兩個電磁懸浮系統(tǒng)控制同一個移動機械橫梁,所以它們之間存在著一定的耦合關系,即當橫梁傾斜或旋轉時,兩個電磁懸浮系統(tǒng)的懸浮氣隙會同時發(fā)生變化,從而兩個電磁懸浮系統(tǒng)的電磁力和其他參數(shù)也會受到影響。以往的文獻忽略了兩個電磁懸浮系統(tǒng)的耦合關系,只是通過設計良好的同步控制器來減小兩個電磁懸浮系統(tǒng)的同步誤差[1]。耦合的存在會降低工件的加工精度,從而降低系統(tǒng)的穩(wěn)定性。
針對兩個系統(tǒng)存在的耦合關系,本文采用α階逆系統(tǒng)解耦控制方法消除耦合對系統(tǒng)穩(wěn)定懸浮的影響,其原理為:將原系統(tǒng)的α階逆系統(tǒng)與原系統(tǒng)串聯(lián)就可將雙電磁耦合懸浮系統(tǒng)解耦成兩個獨立的SISO系統(tǒng),并可將獨立的系統(tǒng)簡化為偽線性系統(tǒng),降低SISO系統(tǒng)控制器的設計難度[2]。逆系統(tǒng)解耦方法需要被控對象數(shù)學模型精確可知,但在工程實踐中很難做到,而且非線性被控對象復雜多變,所以逆系統(tǒng)模型就更加難以建立了。
支持向量機(SVM)的出現(xiàn)為解決非線性系統(tǒng)的逆系統(tǒng)建模難的問題提供了一種有效的方法。支持向量機可以逼近任意非線性函數(shù),這為其用于對系統(tǒng)的辨識提供了理論依據(jù)[3-4]。
本文在分析了龍門數(shù)控加工中心雙電磁懸浮系統(tǒng)耦合情況的基礎上,證明了系統(tǒng)的可逆性。由于被控系統(tǒng)的精確數(shù)學模型難以建立,所以利用支持向量機可以逼近任意非線性函數(shù)的特點,可辨識出被控系統(tǒng)的α階逆系統(tǒng)。得到的逆系統(tǒng)與被控對象串聯(lián)構成偽線性系統(tǒng),從而將MIMO系統(tǒng)解耦成多個SISO。解耦后的單系統(tǒng)采用專家PID閉環(huán)控制,增加了懸浮系統(tǒng)的快速性和魯棒性。仿真實驗表明,支持向量機α階逆系統(tǒng)可以很好地解決系統(tǒng)的耦合問題。
1.1單電磁懸浮系統(tǒng)數(shù)學模型
圖1為龍門移動數(shù)控加工中心機械結構圖。
1.主軸單元 2.懸浮氣隙 3.箱式橫梁 4.X方向直線電機 5.導向單元 6.伺服單元 7.懸浮電磁鐵 8.切削刀具圖1 雙懸浮系統(tǒng)數(shù)控機床部分結構圖
從圖1可以看出,移動橫梁由兩個電磁懸浮系統(tǒng)共同懸浮,為了分析得出兩個電懸浮系統(tǒng)之間的耦合關系,首先需要從建立單電磁懸浮系統(tǒng)數(shù)學模型開始。圖2為單電磁懸浮系統(tǒng)結構圖,圖中,Φ1為漏磁通;Φm為氣隙磁通。
圖2 單電磁懸浮系統(tǒng)結構圖
由文獻[1]得單電磁懸浮系統(tǒng)數(shù)學表達式為
(1)
式中,c(t)為氣隙高度;F(i,c)為電磁力;Fd為外部干擾;m為橫梁的質量;i(t)、u(t)分別為控制電流和電壓;μ0為真空磁導率,μ0=4π×10-7;A為鐵芯面積;N為線圈匝數(shù);R為有效的氣隙磁阻。
由式(1)可以看出,電磁懸浮系統(tǒng)由運動方程、電磁力方程和電壓方程組成,其中電磁力與電流的平方成正比,與懸浮氣隙的平方成反比,因此電磁懸浮系統(tǒng)是典型的非線性系統(tǒng)。
1.2橫梁雙懸浮系統(tǒng)耦合分析
雙電磁懸浮系統(tǒng)耦合示意圖見圖3,圖3中,c1、c2、f1、f2、l、θ分別為兩個電磁鐵的傳感器位移、兩個電磁鐵的電磁力、橫梁的一半長度和橫梁旋轉的角度。
圖3 雙電磁懸浮系統(tǒng)耦合示意圖
其中,轉動角:
(2)
豎直方向上的位移:
(3)
旋轉方向上位移:
(4)
由式(3)、式(4)得
(5)
豎直方向上的合力:
fg=f1+f2
(6)
質心運動動力學方程和繞質心轉動動力學方程分別為
(7)
(8)
(9)
經過力的坐標與傳感器中坐標之間的變換及式(7)、式(9)得加速度與力的關系為
(10)
式(10)表明兩個電磁懸浮系統(tǒng)的加速度存在著耦合關系。
根據(jù)式(1)可以得出電磁力與電壓之間的關系為
(11)
式中,l0為電磁力系數(shù);i1、i2為雙電磁鐵線圈的控制電流;u0為真空磁導率。
通過式(1)和式(10)可以得出系統(tǒng)狀態(tài)方程如下:
(12)
計算輸出變量對時間的導數(shù)可以得到:
(13)
將式(11)代入式(13),式(13)會直接變?yōu)檩斎胱兞繛閡1、u2的方程組。
Jacobi矩陣
(14)
(15)
對于耦合的雙電磁懸浮系統(tǒng),可采用α階逆系統(tǒng)方法進行解耦,其基本思想為:將建立的MIMO被控系統(tǒng)的α階逆系統(tǒng)串聯(lián)在原系統(tǒng)前構成偽線性復合系統(tǒng),原系統(tǒng)會被解耦成多個獨立的SISO系統(tǒng),并且獨立系統(tǒng)具有線性傳遞性質,即偽線性系統(tǒng)。解耦后的獨立子系統(tǒng)可以采用線性控制方法,簡化了系統(tǒng)的結構,適合于工程實踐。
則稱系統(tǒng)Πα為Σ系統(tǒng)的逆系統(tǒng)。
逆系統(tǒng)方法的實現(xiàn)必須滿足兩個條件:第一,被控對象數(shù)學模型精確可知,第二,非線性模型的逆必須求解出來。非線性模型復雜多變,逆系統(tǒng)模型不易求出。由于支持向量機具有逼近任意非線性函數(shù)的功能且具有風險函數(shù)最小化的特點,因此本文采用支持向量機來逼近雙電磁懸浮系統(tǒng)的α階逆系統(tǒng)來達到對耦合系統(tǒng)進行解耦的目的。其基本原理為:通過非線性內積核函數(shù)將數(shù)據(jù)從低維空間Rn映射到Hilbert高維特征空間,然后在高維空間建立一個線性回歸函數(shù),其表達形式為
f(x)=wTφ(x)+b
(16)
式中,w為權值;φ(x)為非線性映射;b為閾值。
當輸入樣本訓練集(xi,ui)時(i=1,2,…,j),其中,xi為懸浮氣隙采集點,ui為依據(jù)結構風險最小化原理,則支持向量機回歸學習最優(yōu)逼近應使得風險函數(shù)最?。?/p>
(17)
通過引用對偶定理、拉格朗日函數(shù)及核函數(shù),將式(17)等價為如下凸二次規(guī)劃問題:
(18)
(19)
最后可以構造非線性回歸函數(shù)[5-6]如下:
(20)
支持向量機擬合α階逆系統(tǒng)流程如圖4所示。
圖4 支持向量機擬合α階逆系統(tǒng)流程
本文設計了時變滑模變結構控制器對解耦后的單懸浮系統(tǒng)進行控制。時變滑模變結構改進了滑模變結構控制或其改進方法在系統(tǒng)的狀態(tài)變量未到達設定的滑模面之前無法實現(xiàn)的不足,使系統(tǒng)在任意初始狀態(tài)下的狀態(tài)變量都能直接到達系統(tǒng)的滑模面上,取消了狀態(tài)變量到達滑模面的過程,以最短的時間實現(xiàn)滑模變結構控制,實現(xiàn)對參數(shù)攝動和外部干擾的全局魯棒性[7]。時變滑模變結構控制器設計如下:
由式(1)可知,該磁懸浮系統(tǒng)是一個典型的非線性系統(tǒng),為設計方便,可以通過相應的坐標變換,將其變換為仿射型非線性系統(tǒng)。
(21)
設位置指令為r,選取位置誤差z1、速度誤差z2、加速度誤差z3為狀態(tài)變量,則誤差表達式可表示為
(22)
因此可得到新坐標系下的磁懸浮系統(tǒng)仿射非線性模型表達式:
(23)
為了討論問題方便,這里作如下假設。
(1)φ(z,t)和d(z,t)滿足:
φmin≤φi(z,t)≤φmax
(24)
dmin≤d(z,t)≤dmax
(25)
(2)b(z,t)滿足:
bmin≤b(z,t)≤bmaxb(z,t)>0
(26)
時變滑模面可設計為
(27)
f(t)=f(0)e-λ tλ>0
(28)
根據(jù)上述時變滑模面選取原則,懸浮系統(tǒng)的時變滑模面設計為
σ(z)=σ1(z)-σ0(z)e-λ t
(29)
σ1(z)=a1z1+a2z2+a3z3
時變滑模變結構控制律由線性控制和切換控制兩部分組成,即ud=uc+uvss。其中,線性控制律為
(30)
(31)
式中,kc為正的常數(shù)。
切換控制律為
uvss=(εs+|λσ0e-λ t|/bmin)sgnσ
(32)
其中,λ為嚴格正常數(shù),σ0為約束函數(shù),式(32)中的增益項εs滿足如下條件:
(33)
雙電磁懸浮系統(tǒng)參數(shù)如下:懸浮橫梁質量m=568 kg;等效質量mθ=0.79 kg;鐵芯面積A=0.0125 m2;電磁鐵繞組匝數(shù)N=340;電阻R=0.65 Ω;期望懸浮氣隙c=0.002 m。本支持向量機核函數(shù)選擇高斯函數(shù),高斯核函數(shù)第一個參數(shù)v=0.05,第二個參數(shù)w=0.01。懲罰因子C=0.6,不敏感損失函數(shù)參數(shù)ε=0.02,時變滑??刂破鞯膮?shù)a1=300,a2=30,a3=1,k=3000,λ=10。仿真系統(tǒng)見圖5。其中封裝的雙電磁懸浮系統(tǒng)見圖6。
圖5 支持向量機α階逆系統(tǒng)解耦系統(tǒng)框圖
圖6 雙電磁懸浮系統(tǒng)仿真框圖
對單懸浮系統(tǒng)施加1000N干擾雙電磁懸浮系統(tǒng)PID控制仿真結果見圖7。
1.受到干擾的電磁懸浮系統(tǒng) 2.未受到干擾的電磁懸浮系統(tǒng)圖7 雙電磁懸浮系統(tǒng)未解耦PID控制仿真曲線
由圖7可得,單電磁懸浮系統(tǒng)受到干擾時另一個系統(tǒng)也會受到影響,這說明兩個電磁懸浮系統(tǒng)存在明顯的耦合關系。
圖8表明支持向量機α階逆系統(tǒng)解耦控制可以很好地使兩個耦合的電磁懸浮系統(tǒng)解耦成兩個獨立的系統(tǒng)互不干擾。圖9表明時變滑模變結構控制較PID控制具有響應速度快、超調小、魯棒性強的特點。
1.受到干擾的電磁懸浮系統(tǒng) 2.未受到干擾的電磁懸浮系統(tǒng)圖8 支持向量機α階逆解耦PID控制仿真曲線
1.受到干擾的電磁懸浮系統(tǒng) 2.未受到干擾的電磁懸浮系統(tǒng)圖9 支持向量機α階逆解耦時變滑模變結構仿真曲線
本文對龍門雙電磁懸浮系統(tǒng)耦合情況進行了分析,得出了兩個電磁懸浮系統(tǒng)的耦合關系,通過可逆性分析,證明了耦合雙電磁懸浮系統(tǒng)是可逆的。為了使耦合的雙電磁懸浮系統(tǒng)解耦出兩個獨立的單系統(tǒng),本文采用α階逆系統(tǒng)對耦合系統(tǒng)進行解耦。由于非線電磁懸浮系統(tǒng)的α階逆系統(tǒng)的精確模型不易得出,因此本文采用支持向量機逼近雙電磁懸浮系統(tǒng)的α階逆系統(tǒng)。解耦后的獨立偽線性系統(tǒng)采用時變滑模變結構控制。仿真結果表明,支持向量機α階逆系統(tǒng)具有很好的解耦控制效果,使兩個懸浮系統(tǒng)互不干擾;時變滑模變結構控制的獨立偽線性系統(tǒng)具有響應速度快、魯棒性強的特點。
[1]劉春芳,邢銀龍,王通.加工中心懸浮系統(tǒng)基于優(yōu)化函數(shù)二級模糊控制[J].組合機床與自動化加工技術,2012,6(6):60-63.
LiuChunfang,XingYinlong,WangTong.TheTwoLevelFuzzyControlofMachiningCenterSuspensionSystemBaseonOptimizationFunction[J].ModularMachineTool&AutomaticManufacturingTechnique,2012,6(6):60-63.
[2]孫曉東,朱火晃秋.基于神經網(wǎng)絡逆系統(tǒng)理論無軸承異步電動機解耦控制[J].電工技術學報,2010,25(1):43-49.
SunXiaodong,ZhuHuangqiu.DecouplingControlofBearinglessInductionMotorsBasedonNeuralNetworkInverseSystemMethod[J].TransactionsofChinaElectrotechnicalSociety,2010,25(1):43-49.
[3]陳杰,朱琳.基于混合最小二乘支持向量機網(wǎng)絡模型的非線性系統(tǒng)辨識[J].控制理論與應用,2010,27(3):303-309.
ChenJie,ZhuLin.NewIdentificationApproachforNonlinearSystemsBasedontheCombinationNetworkModelofLeastSquaresandSupportVectorMachines[J].ControlTheory&Applications,2010,27(3):303-309.
[4]SunYukun,WangBo,JiXiaofu,etal.MultivariableDecouplingControlBasedonFuzzy-neuralNetWorkαth-orderInverseSysteminFermentationProcess[J].ControlTheory&Applications,2010,27(2):188-192.
[5]劉國海,張錦,趙文祥,等.兩電機變頻系統(tǒng)的支持向量機廣義逆內模解耦控制[J].中國電機工程學報,2011,31(6):85-91.
LiuGuohai,ZhangJin,ZhaoWenxiang,etal.InternalModelControlBasedonSupportVectorMachinesGeneralizedInverseforTwo-motorVariableFrequencySystemApplications[J].ProceedingsoftheCSEE,2011,31(6):85-91.
[6]WenXiangjun,ZhangYunong,YanWeiwu,etal.NonlinearDecouplingControllerDesignBasedonLeastSquaresSupportVectorRegression[J].JournalofZhejiangUniversityScienceA,2006, 7(2):275-284.
[7]胡慶雷,馬廣富,姜野,等. 三軸穩(wěn)穩(wěn)定撓性衛(wèi)星姿態(tài)機動時變滑模變結構和主動振動制[J].控制理論與應用,2009,26(2):122-126.
HuQinglei,MaGuangfu,JingYe,etal.VariableStructureControlwithTime-varyingSlidingModeandVibrationControlforFlexibleSatellite[J].ControlTheory&Applications,2009,26(2):122-126.
(編輯王艷麗)
SVMαOrder Inverse System Decoupling Time-varying Sliding Mode Control of Double Suspension Systems of Machining Center
Liu ChunfangRong Gang
Shenyang University of Technology,Shenyang,110870
In order to eliminate the friction between moving beam and guide rail in gantry CNC machine, two suspension systems were adopted to suspend the moving beam jointly, so there were coupling relation between two suspension systems. When the beam suffering unbalance force or external disturbance as the CNC machine working, the existing coupling relationship would influence the two suspension system’s air gap, so the machining accuracy of compent was reduced. This paper analysed and concluded the coupling quantitative relationship when beam rolling around the center of mass and moving up and down.In order to eliminate the influence of coupling this paper adopted SVM to approach the two coupling suspension systems’sαorder inverse system. The obtainedαorder inverse system in series to coupling systems will be decoupled to two SISO pseudo linear systems.The decoupled SISO system adopted time-varying sliding mode control method. The method can make the system reach to the sliding mode surface directly at any original state, eliminating the process of original state reach to the sliding mode surface. So the sliding control can be achieved using shortest time and can realize global robustness of parameter perturbation and external disturbance. The simulation results show that the proposed method can decouple the coupling system effectively and can approach to the controlled object’sαorder inverse system without precise mathematical model. The time-varying sliding mode control method has fast response speed and strong robustness ability.
αorder inverse system; decoupling; support vector machine(SVM); time-varying sliding mode
2013-07-16
國家自然科學基金資助項目(50805098)
TP273< class="emphasis_italic">DOI
:10.3969/j.issn.1004-132X.2015.05.018
劉春芳,女,1975年生。沈陽工業(yè)大學電氣工程學院教授。主要研究方向為伺服系統(tǒng)先進控制策略及應用。榮剛,男,1988年生。沈陽工業(yè)大學電氣工程學院碩士研究生。