賈一君,張 楠,安兆喆,吳世玲

(北京建筑大學 城市軌道交通車輛服役性能保障北京市重點實驗室,北京 102616)

0 引 言

同步是自然界中一種廣泛存在的現(xiàn)象。機械振動系統(tǒng)的同步是指振動系統(tǒng)中兩個或者多個運動部件以相同或者相近的速度或者角速度運行[1]。20世紀60年代,蘇聯(lián)科學家BLEKHMAN I I[2,3]對雙電機的振動同步問題進行了研究,并提出了振動系統(tǒng)自同步理論。

我國學者聞邦椿院士等人[4-7]在自同步理論研究中引入了小參數(shù)平均法,對同一平面內(nèi)雙電機振動同步問題進行了研究,并得出了系統(tǒng)同步性條件與穩(wěn)定性條件。張學良等人[8-11]設計了平面單質體三機及多機驅動自同步振動系統(tǒng)。田曉沖等人[12]針對四電機驅動的振動篩進行了同步性研究及力學分析。陳曉哲等人[13]將機體排布方式從平面引至空間,研究了雙機同軸排布的自同步振動理論。方潘等人[14]對空間三機振動同步系統(tǒng)進行了研究,得到了兩同軸電機與第三電機排布位置及轉子質量對同步性、穩(wěn)定性的影響。

隨著現(xiàn)代同步理論的發(fā)展以及對振動系統(tǒng)同步狀態(tài)精度要求的提高,諸多新型算法被引入振動系統(tǒng)同步控制之中。DONG C F等人[15]將模糊算法與PID參數(shù)整定相結合,設計了基于主從控制結構的系統(tǒng)同步控制器。姜嬌等人[16]采用矢量控制策略及PI控制算法,實現(xiàn)了反向轉動雙機驅動振動篩電機的快速反應。黃志龍等人[17]提出了一種相鄰交叉耦合控制策略,實現(xiàn)了四激振器激勵下,振動機械-物料系統(tǒng)的精確控制。

上述同步控制的實現(xiàn)依賴于大量人工經(jīng)驗和復雜的控制結構,且存在控制精度不高、達到穩(wěn)態(tài)所需時間較長等問題。因此,利用新型智能優(yōu)化算法簡化控制結構、代替人工經(jīng)驗、實現(xiàn)高精度控制等具有一定的實際意義。

麻雀搜索算法(sparrow search algorithm,SSA)主要是根據(jù)自然界中麻雀特有的覓食以及反捕食行為規(guī)律,加以改進得到的一種新型智能算法[18]14-15。

筆者將麻雀算法與直流電機滑?？刂葡嘟Y合,設計出基于主電機轉速和主從電機相位差的同步控制器,對雙機驅動機械振動系統(tǒng)主從電機進行轉速與相位差的高精度控制。

1 系統(tǒng)動力學模型

1.1 雙機驅動振動系統(tǒng)動力學模型

雙機驅動振動系統(tǒng)在結構上由剛性振動殼體、激振器、支撐彈簧組成,激振器1、2呈平面對稱排布,支撐彈簧對稱安裝在固定架上。

在系統(tǒng)穩(wěn)態(tài)工況下,兩個激振器進行速度相等、方向相反的圓周運動。

雙機驅動振動系統(tǒng)動力學模型如圖1所示。

圖1 雙機驅動振動系統(tǒng)動力學模型m1,m2—偏心轉子質量;O1,O2—偏心轉子回轉中心;r1,r2—偏心轉子回轉半徑;φ1,φ2—偏心轉子相位

為簡化分析,筆者在圖1中用偏心轉子表示激振電機。

系統(tǒng)的動能為:

(1)

式中:m—機體質量;Jm—機體轉動慣量;Ji—激振器偏心轉子繞各自回轉中心的轉動慣量;xi,yi—激振器偏心塊在Oxy坐標系中的橫縱坐標。

yi可具體表示為:

xi=x+(licosβi+rcosφi)cosψ+(lisinβi+rsinφi)sinψ
yi=y-(licosβi+rcosφi)sinψ+(lisinβi+rsinφi)cosψ

(2)

式中:βi—振動系統(tǒng)偏心轉子回轉中心到機體質心的連線與x軸夾角。

系統(tǒng)的勢能為:

(3)

式中:ky,ky—x向、y向彈簧剛度;Lx,Ly—機體質心到x向、y方向彈簧連接點的距離。

系統(tǒng)能量散逸函數(shù)為:

(4)

式中:fx,fy—x向、y向彈簧阻尼系數(shù);fi—激振電機各轉軸的阻尼。

對于整個系統(tǒng),其拉格朗日方程為:

(5)

系統(tǒng)的廣義力為:

Qi={Qx,Qy,Qz,Qφ1,Qφ2,Qφ3,Qφ4}T

(6)

其中:Qx=Qy=Qz=0;Qφi=Tei，Tei—激振電機的電磁轉矩。

筆者將求得的系統(tǒng)動能方程、勢能方程以及能量散逸函數(shù)表達式代入拉格朗日方程,可以得到系統(tǒng)x,y,ψ向的運動微分方程和兩個偏心轉子的回轉方程:

(7)

1.2 直流電機數(shù)學模型

在電機選擇方面,考慮到直流電機有結構簡單、制造容易、維護方便、運行可靠等特點,并且相對于交流電機,直流電機在啟動、調(diào)速、制動、協(xié)調(diào)運轉方面的性能更優(yōu),筆者采用直流電機驅動方式,以滿足實際應用的復雜性與合理性。

為了主要研究振動系統(tǒng)穩(wěn)定時的同步情況,筆者將直流電機的電樞電阻設置為常值函數(shù),直流電機機械特性方程如下:

(8)

式中:Tm—電機電磁轉矩;Tf—電機負載轉矩;J—電樞轉動慣量;ω—電機轉速。

電磁轉矩與電樞電流有如下關系:

Tm=Ktφia

(9)

式中:φ—電機磁極的磁通量;Kt—電機結構常數(shù);ia—電樞電流。

2 麻雀優(yōu)化變結構滑模控制器

雙機驅動振動系統(tǒng)的控制系統(tǒng)主要包括雙機驅動機械振動系統(tǒng)動力學模塊、主從控制器、直流電機-偏心轉子模塊。

雙機驅動振動系統(tǒng)控制系統(tǒng)框圖如圖2所示。

圖2 雙機驅動振動系統(tǒng)控制系統(tǒng)框圖

首先,筆者設定目標轉速ω,直流電機1(即主電機)轉速ω1與目標轉速ω相比較后,經(jīng)主電機角速度控制器處理得到主電機電流信號i1,電流信號經(jīng)主直流電機-偏心轉子模塊得到電磁轉矩Te1、轉子相位φ1等相關參數(shù)信號，進入雙機驅動振動系統(tǒng)力學模塊;

同時,直流電機2(即從電機)轉速積分后得到相位φ2,與主電機相位φ1相比較后,經(jīng)從電機相位控制器處理得到從電機電流信號i2,電流信號經(jīng)從直流電機-偏心轉子模塊得到電磁轉矩Te2、轉子相位φ2等相關參數(shù)信號，進入雙機驅動振動系統(tǒng)力學模塊。力學系統(tǒng)解算得到系統(tǒng)x、y、ψ運動信號返回直流電機-偏心轉子模塊。

2.1 指數(shù)趨近率滑模

滑模變結構控制是一種針對復雜控制系統(tǒng)設計的變結構控制策略,可以實現(xiàn)被控系統(tǒng)在切換面沿設定軌跡滑動直至趨于穩(wěn)定的效果?；＿\動包含兩個階段,即趨近運動和滑動模態(tài)。系統(tǒng)從任意初始狀態(tài)向切換面趨近,最終運動至切換面的過程叫做趨近運動,即s→0。

根據(jù)滑模變結構原理,滑?？蛇_性條件只能確保任意一位置的運動點在空間中的有限時間內(nèi)到達切換面,但并不規(guī)定其具體的運動軌跡。因此,可以嘗試利用趨近率的方法確保趨近運動過程中的動態(tài)品質。

指數(shù)趨近率的作用可以看作:使趨近速度逐漸由一個相對較大值向0趨近,這個過程所耗時相對較少,且運動點可以以一個相對較小的速度到達切換面,采用如下形式的指數(shù)趨近率:

(10)

式中:E1，k1—均為正數(shù)。

2.2 主電機轉速控制器設計

取主電機狀態(tài)變量為:

(11)

式中:ω*—給定轉速;ω1—主電機轉速。

將式(8，9)代入式(11),并對其進行求導,可得:

(12)

(13)

系統(tǒng)滑模面為:

s=c1x1+x2

(14)

c1為正數(shù),滿足Hurwitz條件。

由式(10，14)可得控制量i1表達式:

(15)

2.3 從電機相位控制器設計

取從電機狀態(tài)變量為:

(16)

系統(tǒng)滑模面為:

(17)

同樣,采用指數(shù)趨近率,有:

(18)

式中:E2，k2—均為正數(shù)。

由以上式子可得控制量表達式:

(19)

2.4 主電機麻雀優(yōu)化變結構滑?？刂?/h3>

麻雀算法具有尋優(yōu)能力強、收斂速度快等優(yōu)點,是一種新穎的智能算法。利用麻雀算法對式(15)中E1、k1進行參數(shù)尋優(yōu),可以取得較好的抑抖效果和趨近速度。

設計適應度函數(shù)為:

F=0.5|s|

(20)

在麻雀群體中有發(fā)現(xiàn)者與加入者兩種類型的麻雀,發(fā)現(xiàn)者負責尋找食物并報警給群體提供覓食方向,加入者根據(jù)發(fā)現(xiàn)者提供的信息不斷更新自身位置直至獲取食物[18]16-18。

假設群體中共有n只麻雀,待優(yōu)化變量數(shù)量為m,則群體初始位置X表示為如下形式:

(21)

發(fā)現(xiàn)者根據(jù)如下規(guī)則不斷向最優(yōu)位置Xi,j搜索迫近:

(22)

式中:t—迭代次數(shù);itermax—最大迭代次數(shù);R—預警值;ST—安全值;Q—服從正態(tài)分布的隨機數(shù);L—元素全部為1的1×m的矩陣。

其中:i=1,2,3,…,n;j=1,2,3,…,m。

加入者根據(jù)如下規(guī)則不斷更新位置Yi,j:

(23)

式中:Yw—全局最差位置;Yp—當前搜索條件下全局最優(yōu)位置。

其優(yōu)化流程如圖3所示。

圖3 麻雀優(yōu)化變結構滑?？刂葡到y(tǒng)框圖

經(jīng)過一定次數(shù)的迭代后,麻雀算法給出當前時刻下最優(yōu)的切換項系數(shù)E1、K1。

3 同步控制系統(tǒng)仿真分析

為了驗證該系統(tǒng)控制同步設計的合理性,筆者利用Simulink搭建電機-振動系統(tǒng)仿真模型,如圖4所示。

圖4 雙機驅動振動系統(tǒng)同步控制仿真模型

采用龍格-庫塔算法和變步長控制策略,系統(tǒng)仿真參數(shù)如表1所示。

表1 振動系統(tǒng)仿真參數(shù)

當主、從電機角速度與相位差均趨于穩(wěn)定,則振動系統(tǒng)進入同步狀態(tài)[19]。

同步控制仿真實驗結果如圖5所示。

由圖5(a)可以看出,主電機角速度(圖5(a)中用實線表示)在經(jīng)歷0 s～20 s的調(diào)整階段后,平滑到達預設定的角速度20 rad/s;從電機角速度曲線(圖中用虛線表示)在滑?？刂破髯饔孟赂S主電機角速度平滑達到預設定的角速度,且實現(xiàn)了較好的跟蹤效果;

圖5 同步控制仿真實驗結果

圖5(b)顯示兩電機相位差最初為-2 rad,在滑?？刂破鞯淖饔孟?兩電機相位差在經(jīng)歷一個較大波動后也實現(xiàn)穩(wěn)定,最終穩(wěn)定在0 rad左右;

圖5(c,d)為機體x、y方向位移。圖5(c)表明:在經(jīng)歷了約40 s的調(diào)整階段后,機體x方向位移最終趨于0,振動系統(tǒng)x方向趨于靜止;

圖5(d)表明:在經(jīng)歷了約30 s的調(diào)整階段后,機體y方向位移最終趨于-0.02 m～0.02 m的弦類往復波動;

圖5(e)為機體整體位移,從中可以看出,在經(jīng)歷一段時間的無規(guī)律波動后,振動系統(tǒng)逐漸過渡為x方向位移為0、y方向位移為-0.02 m～0.02 m的弦類往復波動(圖中線條密集部分);

圖5(f)為振動機體Ψ方向擺動,最終為小范圍內(nèi)的穩(wěn)態(tài)弦類波動。雙機驅動振動系統(tǒng)主要應用于直線式往復振動工作場合。

圖5結果表明:系統(tǒng)設計滿足生產(chǎn)實際要求,且主從電機角速度、相位差以及振動系統(tǒng)整體運動均處于穩(wěn)定狀態(tài),可見雙機驅動振動系統(tǒng)的麻雀滑模同步控制是有效的、穩(wěn)定的。

4 結束語

筆者應用Lagrange方程建立了雙直流電機驅動型振動系統(tǒng)的機電耦合模型,設計了基于主電機角速度差與主從電機相位差的主從變結構滑?？刂破?并利用麻雀優(yōu)化算法對主電機滑模切換項系數(shù)進行了優(yōu)化。

研究結果表明:

(1)采用麻雀優(yōu)化滑?？刂扑惴梢詫⒅麟姍C轉速與預設轉速之差控制在0.5 rad/s之內(nèi);從電機轉速跟蹤效果較好,主從電機相位差可穩(wěn)定在0.1 rad之內(nèi);

(2)在主從電機控制器作用下,振動機體水平方向基本維持不動,豎直方向表現(xiàn)為-0.02 m～0.002 m的穩(wěn)態(tài)弦類波動,滿足振動類機械實際工作的要求。

實際研究發(fā)現(xiàn),在振動機械的生產(chǎn)作業(yè)中,為滿足差異化工藝需求,往往需要3臺及以上振動電機同步工作。

在后續(xù)的研究中,筆者計劃將麻雀滑?？刂扑惴☉糜趯Χ嗯_電機的同步控制之中。