劉春芳 榮 剛
沈陽工業大學,沈陽,110870
加工中心雙懸浮系統支持向量機α階逆解耦時變滑模控制
劉春芳榮剛
沈陽工業大學,沈陽,110870
龍門數控加工中心移動機械橫梁可采用雙電磁懸浮系統共同懸浮來消除摩擦的影響,而這兩個電磁懸浮系統存在著一定的耦合關系。當數控機床工作時,由于受力不平衡或者擾動等因素影響,這種耦合的存在會使兩個電磁懸浮系統的懸浮氣隙受到影響,并降低加工精度。分析得出了橫梁發生繞質心旋轉和上下平移時兩個電磁懸浮系統的耦合定量關系,為了消除耦合的影響,采用支持向量機逼近雙電磁懸浮系統的α階逆系統,將所得到的α階逆系統串聯在原系統前構成偽線性復合系統,從而將原系統解耦成兩個獨立的SISO偽線性系統。針對解耦后的偽線性SISO系統設計了時變滑模變結構控制器,可使系統在任意初始狀態下系統狀態變量都能直接到達系統的滑模面上,消除了狀態變量到達滑模面的過程,以最短的時間實現了滑模變結構控制,實現了對參數攝動和外部干擾的全局魯棒性。仿真實驗結果表明,該方案不需要被控對象精確的數學模型即可采用支持向量機逼近被控系統的α階逆系統,可有效地對耦合的懸浮系統實現解耦。解耦后的SISO偽線性系統采用時變滑模變結構控制具有響應速度快和魯棒性強等特點。
α階逆系統;解耦;支持向量機;時變滑模變結構
隨著科學技術的發展,人們對零件加工精度的要求也越來越高。傳統的龍門移動數控機床移動橫梁與導軌完全接觸存在著摩擦,即使加入潤滑劑也不能完全消除摩擦影響。為了消除摩擦的影響,利用兩個電磁懸浮系統將移動橫梁懸浮起來,可提高加工精度。由于兩個電磁懸浮系統控制同一個移動機械橫梁,所以它們之間存在著一定的耦合關系,即當橫梁傾斜或旋轉時,兩個電磁懸浮系統的懸浮氣隙會同時發生變化,從而兩個電磁懸浮系統的電磁力和其他參數也會受到影響。以往的文獻忽略了兩個電磁懸浮系統的耦合關系,只是通過設計良好的同步控制器來減小兩個電磁懸浮系統的同步誤差[1]。耦合的存在會降低工件的加工精度,從而降低系統的穩定性。
針對兩個系統存在的耦合關系,本文采用α階逆系統解耦控制方法消除耦合對系統穩定懸浮的影響,其原理為:將原系統的α階逆系統與原系統串聯就可將雙電磁耦合懸浮系統解耦成兩個獨立的SISO系統,并可將獨立的系統簡化為偽線性系統,降低SISO系統控制器的設計難度[2]。逆系統解耦方法需要被控對象數學模型精確可知,但在工程實踐中很難做到,而且非線性被控對象復雜多變,所以逆系統模型就更加難以建立了。
支持向量機(SVM)的出現為解決非線性系統的逆系統建模難的問題提供了一種有效的方法。支持向量機可以逼近任意非線性函數,這為其用于對系統的辨識提供了理論依據[3-4]。
本文在分析了龍門數控加工中心雙電磁懸浮系統耦合情況的基礎上,證明了系統的可逆性。由于被控系統的精確數學模型難以建立,所以利用支持向量機可以逼近任意非線性函數的特點,可辨識出被控系統的α階逆系統。得到的逆系統與被控對象串聯構成偽線性系統,從而將MIMO系統解耦成多個SISO。解耦后的單系統采用專家PID閉環控制,增加了懸浮系統的快速性和魯棒性。仿真實驗表明,支持向量機α階逆系統可以很好地解決系統的耦合問題。
1.1單電磁懸浮系統數學模型
圖1為龍門移動數控加工中心機械結構圖。

1.主軸單元 2.懸浮氣隙 3.箱式橫梁 4.X方向直線電機 5.導向單元 6.伺服單元 7.懸浮電磁鐵 8.切削刀具圖1 雙懸浮系統數控機床部分結構圖
從圖1可以看出,移動橫梁由兩個電磁懸浮系統共同懸浮,為了分析得出兩個電懸浮系統之間的耦合關系,首先需要從建立單電磁懸浮系統數學模型開始。圖2為單電磁懸浮系統結構圖,圖中,Φ1為漏磁通;Φm為氣隙磁通。

圖2 單電磁懸浮系統結構圖
由文獻[1]得單電磁懸浮系統數學表達式為
(1)
式中,c(t)為氣隙高度;F(i,c)為電磁力;Fd為外部干擾;m為橫梁的質量;i(t)、u(t)分別為控制電流和電壓;μ0為真空磁導率,μ0=4π×10-7;A為鐵芯面積;N為線圈匝數;R為有效的氣隙磁阻。
由式(1)可以看出,電磁懸浮系統由運動方程、電磁力方程和電壓方程組成,其中電磁力與電流的平方成正比,與懸浮氣隙的平方成反比,因此電磁懸浮系統是典型的非線性系統。
1.2橫梁雙懸浮系統耦合分析
雙電磁懸浮系統耦合示意圖見圖3,圖3中,c1、c2、f1、f2、l、θ分別為兩個電磁鐵的傳感器位移、兩個電磁鐵的電磁力、橫梁的一半長度和橫梁旋轉的角度。

圖3 雙電磁懸浮系統耦合示意圖
其中,轉動角:
(2)
豎直方向上的位移:
(3)
旋轉方向上位移:
(4)
由式(3)、式(4)得
(5)
豎直方向上的合力:
fg=f1+f2
(6)
質心運動動力學方程和繞質心轉動動力學方程分別為
(7)
(8)

(9)
經過力的坐標與傳感器中坐標之間的變換及式(7)、式(9)得加速度與力的關系為
(10)
式(10)表明兩個電磁懸浮系統的加速度存在著耦合關系。

根據式(1)可以得出電磁力與電壓之間的關系為
(11)
式中,l0為電磁力系數;i1、i2為雙電磁鐵線圈的控制電流;u0為真空磁導率。
通過式(1)和式(10)可以得出系統狀態方程如下:
(12)
計算輸出變量對時間的導數可以得到:
(13)
將式(11)代入式(13),式(13)會直接變為輸入變量為u1、u2的方程組。
Jacobi矩陣
(14)


(15)
對于耦合的雙電磁懸浮系統,可采用α階逆系統方法進行解耦,其基本思想為:將建立的MIMO被控系統的α階逆系統串聯在原系統前構成偽線性復合系統,原系統會被解耦成多個獨立的SISO系統,并且獨立系統具有線性傳遞性質,即偽線性系統。解耦后的獨立子系統可以采用線性控制方法,簡化了系統的結構,適合于工程實踐。


則稱系統Πα為Σ系統的逆系統。
逆系統方法的實現必須滿足兩個條件:第一,被控對象數學模型精確可知,第二,非線性模型的逆必須求解出來。非線性模型復雜多變,逆系統模型不易求出。由于支持向量機具有逼近任意非線性函數的功能且具有風險函數最小化的特點,因此本文采用支持向量機來逼近雙電磁懸浮系統的α階逆系統來達到對耦合系統進行解耦的目的。其基本原理為:通過非線性內積核函數將數據從低維空間Rn映射到Hilbert高維特征空間,然后在高維空間建立一個線性回歸函數,其表達形式為
f(x)=wTφ(x)+b
(16)
式中,w為權值;φ(x)為非線性映射;b為閾值。
當輸入樣本訓練集(xi,ui)時(i=1,2,…,j),其中,xi為懸浮氣隙采集點,ui為依據結構風險最小化原理,則支持向量機回歸學習最優逼近應使得風險函數最小:
(17)

通過引用對偶定理、拉格朗日函數及核函數,將式(17)等價為如下凸二次規劃問題:
(18)


(19)
最后可以構造非線性回歸函數[5-6]如下:

(20)
支持向量機擬合α階逆系統流程如圖4所示。

圖4 支持向量機擬合α階逆系統流程
本文設計了時變滑模變結構控制器對解耦后的單懸浮系統進行控制。時變滑模變結構改進了滑模變結構控制或其改進方法在系統的狀態變量未到達設定的滑模面之前無法實現的不足,使系統在任意初始狀態下的狀態變量都能直接到達系統的滑模面上,取消了狀態變量到達滑模面的過程,以最短的時間實現滑模變結構控制,實現對參數攝動和外部干擾的全局魯棒性[7]。時變滑模變結構控制器設計如下:
由式(1)可知,該磁懸浮系統是一個典型的非線性系統,為設計方便,可以通過相應的坐標變換,將其變換為仿射型非線性系統。

(21)
設位置指令為r,選取位置誤差z1、速度誤差z2、加速度誤差z3為狀態變量,則誤差表達式可表示為
(22)
因此可得到新坐標系下的磁懸浮系統仿射非線性模型表達式:
(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)
根據上述時變滑模面選取原則,懸浮系統的時變滑模面設計為
σ(z)=σ1(z)-σ0(z)e-λ t
(29)
σ1(z)=a1z1+a2z2+a3z3

時變滑模變結構控制律由線性控制和切換控制兩部分組成,即ud=uc+uvss。其中,線性控制律為
(30)
(31)
式中,kc為正的常數。
切換控制律為
uvss=(εs+|λσ0e-λ t|/bmin)sgnσ
(32)
其中,λ為嚴格正常數,σ0為約束函數,式(32)中的增益項εs滿足如下條件:

(33)
雙電磁懸浮系統參數如下:懸浮橫梁質量m=568 kg;等效質量mθ=0.79 kg;鐵芯面積A=0.0125 m2;電磁鐵繞組匝數N=340;電阻R=0.65 Ω;期望懸浮氣隙c=0.002 m。本支持向量機核函數選擇高斯函數,高斯核函數第一個參數v=0.05,第二個參數w=0.01。懲罰因子C=0.6,不敏感損失函數參數ε=0.02,時變滑模控制器的參數a1=300,a2=30,a3=1,k=3000,λ=10。仿真系統見圖5。其中封裝的雙電磁懸浮系統見圖6。

圖5 支持向量機α階逆系統解耦系統框圖

圖6 雙電磁懸浮系統仿真框圖
對單懸浮系統施加1000N干擾雙電磁懸浮系統PID控制仿真結果見圖7。

1.受到干擾的電磁懸浮系統 2.未受到干擾的電磁懸浮系統圖7 雙電磁懸浮系統未解耦PID控制仿真曲線
由圖7可得,單電磁懸浮系統受到干擾時另一個系統也會受到影響,這說明兩個電磁懸浮系統存在明顯的耦合關系。
圖8表明支持向量機α階逆系統解耦控制可以很好地使兩個耦合的電磁懸浮系統解耦成兩個獨立的系統互不干擾。圖9表明時變滑模變結構控制較PID控制具有響應速度快、超調小、魯棒性強的特點。

1.受到干擾的電磁懸浮系統 2.未受到干擾的電磁懸浮系統圖8 支持向量機α階逆解耦PID控制仿真曲線

1.受到干擾的電磁懸浮系統 2.未受到干擾的電磁懸浮系統圖9 支持向量機α階逆解耦時變滑模變結構仿真曲線
本文對龍門雙電磁懸浮系統耦合情況進行了分析,得出了兩個電磁懸浮系統的耦合關系,通過可逆性分析,證明了耦合雙電磁懸浮系統是可逆的。為了使耦合的雙電磁懸浮系統解耦出兩個獨立的單系統,本文采用α階逆系統對耦合系統進行解耦。由于非線電磁懸浮系統的α階逆系統的精確模型不易得出,因此本文采用支持向量機逼近雙電磁懸浮系統的α階逆系統。解耦后的獨立偽線性系統采用時變滑模變結構控制。仿真結果表明,支持向量機α階逆系統具有很好的解耦控制效果,使兩個懸浮系統互不干擾;時變滑模變結構控制的獨立偽線性系統具有響應速度快、魯棒性強的特點。
[1]劉春芳,邢銀龍,王通.加工中心懸浮系統基于優化函數二級模糊控制[J].組合機床與自動化加工技術,2012,6(6):60-63.
LiuChunfang,XingYinlong,WangTong.TheTwoLevelFuzzyControlofMachiningCenterSuspensionSystemBaseonOptimizationFunction[J].ModularMachineTool&AutomaticManufacturingTechnique,2012,6(6):60-63.
[2]孫曉東,朱火晃秋.基于神經網絡逆系統理論無軸承異步電動機解耦控制[J].電工技術學報,2010,25(1):43-49.
SunXiaodong,ZhuHuangqiu.DecouplingControlofBearinglessInductionMotorsBasedonNeuralNetworkInverseSystemMethod[J].TransactionsofChinaElectrotechnicalSociety,2010,25(1):43-49.
[3]陳杰,朱琳.基于混合最小二乘支持向量機網絡模型的非線性系統辨識[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]劉國海,張錦,趙文祥,等.兩電機變頻系統的支持向量機廣義逆內模解耦控制[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]胡慶雷,馬廣富,姜野,等. 三軸穩穩定撓性衛星姿態機動時變滑模變結構和主動振動制[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年生。沈陽工業大學電氣工程學院教授。主要研究方向為伺服系統先進控制策略及應用。榮剛,男,1988年生。沈陽工業大學電氣工程學院碩士研究生。