張思進+杜偉霞+殷珊

摘 要：振動篩系統(tǒng)是一類非光滑度很高的多參數(shù)非線性動力系統(tǒng)，傳統(tǒng)的分岔準則無法直接適用，這里采用了新的不依賴于特征值計算的顯式分岔臨界準則，以實現(xiàn)振動篩系統(tǒng)雙Hopf分岔的反控制.首先，根據(jù)系統(tǒng)的運動方程得到Poincaré映射在不動點處的線性化矩陣；然后，對系統(tǒng)施加線性反饋控制器，得到受控后的Poincaré映射，根據(jù)分岔臨界準則求得雙Hopf分岔的顯式臨界條件；最后，通過模態(tài)疊加法對兩類系統(tǒng)分別進行了數(shù)值模擬.結(jié)果顯示，在相同的系統(tǒng)參數(shù)下線性反饋控制器通過調(diào)整控制參數(shù)，可以有效地實現(xiàn)雙Hopf分岔的反控制.雙Hopf分岔可以提高一些振動機械的工作效率，具有一定的實際意義.

關鍵詞：雙Hopf分岔反控制；非光滑；顯式分岔臨界準則；線性反饋控制器

中圖分類號：O322； TH113 文獻標志碼：A

Anti-control of Double Hopf Bifurcation of VibrationRating Griddle System

ZHANG Sijin，DU Weixia， YIN Shan

（College of Mechanical and Vehicle Engineering，Hunan University，Changsha 410082，China）

Abstract：The vibrating griddle system is a kind of nonlinear multi-parameter dynamical system with high non-smoothness，and traditional bifurcation criterion cant be applied directly. Therefore，a new explicit bifurcation criterion independent on eigenvalue calculation was adopted to achieve the anti-control of double Hopf bifurcation for vibrating griddle system. Firstly，the Poincaré mapping at the fixed point was obtained according to the motion equation of the system. Then，a linear feedback controller was applied to the system to obtain the controlled Poincaré mapping，and the explicit critical condition of double Hopf bifurcation was obtained according to the bifurcation critical criterion. Finally，numerical simulation of the two types of systems was carried out by the modal superposition method. The results show that the linear feedback controller can effectively achieve the anti-control of double Hopf bifurcation by adjusting the control parameters under the same system parameters. Double Hopf bifurcation can improve the working efficiency of some vibration machinery，so it has certain practical significance.

Key words：double Hopf bifurcation anti-control；non-smooth；explicit critical criterion of bifurcation；linear feedback controller

近幾十年來，由于振動篩結(jié)構簡單、處理能力強、工作可靠，得到了很快的發(fā)展，并廣泛應用在煤炭、采礦、建筑、冶金等多種行業(yè).對于振動篩的研究基本上是以生產(chǎn)實際為導向，力求提高篩分效率，提高產(chǎn)品質(zhì)量，降低能源消耗.文獻[1-2]采用數(shù)值分析與試驗相結(jié)合的方法對其動態(tài)特性和內(nèi)部應力進行了分析，結(jié)果表明振動篩固有頻率和工作頻率保持在合理范圍內(nèi)可避免共振.文獻[3-4]運用多軸疲勞分析的臨界平面法對大型直線振動篩進行了疲勞壽命數(shù)值分析，運用應力-強度干涉模型對其可靠性進行了預測.文獻[5]記錄了振動篩的研究進展情況，提出了盡管研究已經(jīng)進展得很好，但是還有很多未解決的問題，像物料與篩面碰撞中的分岔、混沌等非線性現(xiàn)象.

由于振動篩的強非線性和非光滑性，使得振動篩的動力學行為相當復雜，因此振動篩的分岔和混沌行為成了近年來非線性動力學領域的熱門研究課題.白亮亮等[6]對振動篩系統(tǒng)的非線性響應進行了分析，在進行仿真后，發(fā)現(xiàn)了其存在Flip分岔、Neimark-Sacker分岔、環(huán)面倍化分岔以及Hopf-Hopf分岔，并向混沌演化，提出了當振動篩與物料有相同的振動頻率時，物料的分離效果最佳.張永祥等[7]基于Poincaré 映射研究了系統(tǒng)的余維三分岔和非常規(guī)的混沌演化過程，在分岔點附近發(fā)現(xiàn)了三角形吸引子、“五角星型”、“輪胎型”概周期吸引子.雖然近年來不少

研究者開始研究分岔的反控制[8-13]，但是有關振動篩的分岔反控制研究還是很少，文獻[14-15]研究了二自由度振動落砂機周期倍化分岔的反控制.endprint

本文通過振動篩的力學模型，對物料與篩面碰撞中的分岔進行了探究，并把文獻[11]提出的反控制理論應用在振動篩上，通過數(shù)值計算，得出原系統(tǒng)在線性反饋控制器的作用下通過調(diào)節(jié)控制參數(shù)得到雙Hopf分岔解，實現(xiàn)了雙Hopf分岔的反控制.對于一些振動系統(tǒng)像落砂機、篩分機、振動棒、攪拌機等，雙Hopf分岔產(chǎn)生的環(huán)面概周期運動可以提高該類系統(tǒng)的工作效率，比如可以使攪拌機攪拌得更均勻，使篩分機篩分得更快，因此，實現(xiàn)雙Hopf分岔的反控制具有一定的實際意義.

1 運動方程以及Poincaré映射

振動篩系統(tǒng)的力學簡化模型如圖1所示，質(zhì)量為M1和M2的振子分別由剛度為K1和K2的兩個線性彈簧以及阻尼系數(shù)為C1和C2的線性阻尼器相連接，振子只在垂直方向上運動，對M1施加簡諧激振力Psin（Ωt+τ）.質(zhì)量塊M3在重力作用下從上方落下并與M1產(chǎn)生碰撞，改變速度后以新的初值運動，如此往復.設該模型里的阻尼是Rayleigh型比例阻尼，碰撞恢復系數(shù)為R.

振動篩系統(tǒng)的運動微分方程是：

M11+（C1+C2）1-C22+（K1+K2）q1-K2q2=Psin（ΩT+τ），

M22-C21+C22-K2q1+K2q2=0，3=-g.（1）

式中：q1，q2，q3分別是振子M1，M2，M3的位移；M1和M3之間的間隙是e.

將以上微分運動方程無量綱化為：

100m12+2ζ1+c-c-kk12

+1+k-k-kkx1x2=10sin（ωt+τ），3=-l（2）

其中無量綱量為：

m=M2/M1，k=K2/K1，c=C2/C1，

ζ=C1/（2M1ωn1），ωn1=K1/M1，

xi=qiK1/P，ω=ΩM1K1，t=TK1M1，

l=M1gP，δ=eK1/P，γ=M3/M1，

（i=1，2，3）

當q1-q3=e時，即x1-x3=δ，系統(tǒng)將發(fā)生碰撞.由碰撞前后系統(tǒng)動量守恒和恢復系數(shù)的定義，振子M1與M3碰撞前后的速度關系為：

1-+γ3-=1++γ3+，

1+-3+=-R（1--3-），（3）

式中：1-，1+和3-，3+分別表示振子M1，M3碰撞前后的瞬時速度；R表示碰撞恢復系數(shù).

設振子碰撞前的瞬時為0時刻，則下一次振子碰撞前的瞬時為2nπ/ω，n∈Z.系統(tǒng)周期運動的初始條件是：

x1（0）=x3（0）=x1（2nπ/ω）

=x3（2nπ/ω）=x10，

x2（0）=x2（2nπ/ω）=x20，

2（0）=2（2nπ/ω）=20，

1（0）=1-γR1+γ1（2nπ/ω）+

γ（1+R）1+γ3（2nπ/ω）=10，

3（0）=1+R1+γ1（2nπ/ω）+

γ-R1+γ3（2nπ/ω）=30（4）

下面對系統(tǒng)進行解耦，得到其周期解的表達式.方程（2）可以表示成如下形式：

M+C+Kx=Fsin（ωt+τ），

3=l.（5）

式中：M=[1 0；0 m]，C=2ζK，F(xiàn)=[1；0]，

K=[1+k -k；-k k].

令Ψ表示方程（5）的正則模態(tài)矩陣，ω1，ω2表示系統(tǒng)的固有頻率，經(jīng)坐標變換得：

x=Ψξ，ξ=[ξ1 ξ2]T，（6）

將（6）式代入（5）式，解耦得：

[I]+2ζ[Λ]+[Λ]ξ=sin（ωt+τ），

3=-l.（7）

式中：2ζΛ=diag[2ζω21，2ζω22]，=ΨTF，

Λ=diag[ω21，ω22].

由模態(tài)疊加法可得方程（7）的通解：

xi=∑2j=1Ψij（e-ηjt（ajcosωdjt+bjsinωdjt）+

Ajsin（ωt+τ）+Bjcos（ωt+τ）），

i=∑2j=1Ψij（e-ηjt（（bjωdjt-ηjaj）cosωdjt-

（ηjbj+ajωdj）sinωdjt）+ωAjcos（ωt+τ）-

Bjω（sin（ωt+τ）），（i=1，2），

x3=C2+C1t-l2t2，

3=C1-lt.（8）

式中：Ψij是正則模態(tài)矩陣Ψ的元素；

ηj=ζω2j，ωdj=ω2j-η2j；aj，bj是積分常數(shù)，由系統(tǒng)的初始條件和模態(tài)參數(shù)確定；Aj，Bj是振幅常數(shù)：

Aj=ω2j-ω2（ω2j-ω2）2+（2ηjω）2j，

Bj=-2ηjω（ω2j-ω2）2+（2ηjω）2j.（9）

這里j=∑2k=1Ψkjfk0，f10=F1，f20=0.

選擇以下截面：

∑={（x1，1，x2，2，x3，3，θ）∈R6×

S1|x1=x3，1=1+，3=3+}，

作為系統(tǒng)的Poincaré截面，其中θ=ωtmod（2π），則X*=（x10，10，x20，20，30，τ）T

表示系統(tǒng)在Poincaré 截面上的周期不動點.參考文獻[7]可建立Poincaré 映射：

Xk+1=f（μ，X*）（10）

式中：μ=（ζ，ω）為分岔參數(shù).

Poincaré 映射（10）在不動點X*處的線性化矩陣為：

Df（μ，X*）=f（σμ，ΔX）ΔX|（μ，X*），（11）endprint

式中：ΔX為X的擾動矢量.

2 系統(tǒng)的雙Hopf分岔反控制

2.1 反饋控制系統(tǒng)

對系統(tǒng)（2）施加線性反饋控制器：

M+C+Kx+β（-p）+

α（x-xp）=Fsin（ωt+τ），

3=-l（12）

1-+γ3-=1++γ3+，

1+-3+=-R（1--3-）（13）

式中：α和β是線性反饋控制增益矩陣.由于該控制方法并不改變原系統(tǒng)的周期解，所以xp=（x1p，x2p）T和p=（1p，2p）T仍然是原系統(tǒng)的周期解位移和速度.令y=x-xp，=-p，則（12）式可化為：

M*+C*+K*y=

F（sin（ωt+τ）-sin（ωt+τ0））（14）

式中：M*=M，C*=C+β，K*=K+α通過適當?shù)淖儞Q可以得到（12）式的通解：

xi=∑2j=1Φij{（e-jt（jcosdjt+jsindjt）+

j[sin（ωt+τ）-sin（ωt+τ0）]+

j[cos（ωt+τ）-cos（ωt+τ0）]}+xip，

i=∑2j=1Φij{（e-jt[（jdjt-jj）cosdjt-

（jj+jdj）sindjt）]+

jω[cos（ωt+τ）-cos（ωt+τ0）]-

jω[sin（ωt+τ）-sin（ωt+τ0）]}+ip，

x3=C4+C3t-l2t2，

3=C3-lt.（15）

式中：Φij是方程（14）的正則模態(tài)矩陣；j=ζ2j，dj=2j-2j，j是受控系統(tǒng)的固有頻率；j，j，C3，C4，τ是積分常數(shù)，由系統(tǒng)的初始條件和周期性條件確定；j，j是振幅常數(shù)，由下式求出：

j=2j-ω2（2j-ω2）2+（2jω）2j，

j=-2jω（2j-ω2）2+（2jω）2j，

j=∑2k=1Φkjfk0.（16）

2.2 受控系統(tǒng)的Poincaré 映射

因為該反饋控制方法并不改變原系統(tǒng)的周期解，故選擇Poincaré 截面：

σ′={（x1，1，x2，2，x3，3，θ）∈R6×

S|x1=x3，1=1+，3=3+}，

式中：θ=ωtmod（2π）；S=Rmod（2π）；R表示實數(shù).基于原周期解可以建立Poincaré 映射為：

X′k+1=H（μ，ε，Xk），（17）

式中：Xk=（x1k，1k，x2k，2k，3k，τk）T，μ和ε=（k1，k2）分別是分岔參數(shù)和控制參數(shù).

2.3 受控系統(tǒng)產(chǎn)生雙Hopf分岔的顯式臨界條件

傳統(tǒng)分岔的臨界準則，需要在參數(shù)空間內(nèi)逐點取值以求出特征值并驗證該特征值是否滿足分岔的臨界準則，以此來確定控制參數(shù)，這樣既盲目、具有不確定性又需要計算很長時間.而通過極點配置方法獲得的控制參數(shù)，其對于控制參數(shù)的物理意義、控制魯棒性和橫截條件存在機理不明確問題.所以，本文參考文獻[12]推導出不直接依賴于特征值計算的發(fā)生雙Hopf分岔的顯式臨界條件，以獲得符合條件的控制參數(shù).

設映射在不動點X*=（x10，10，x20，20，30，τ0）T處的線性化矩陣為：

DXH（μ，ε，X*）=H（μ，ε，X*）X*=Q（μ，ε），（18）

其特征多項式有以下形式：

Pμ（λ）=λ6+m1λ5+m2λ4+m3λ3+m4λ2+m5λ+m6，（19）

式中：mi=mi（μ，ε）是與分岔參數(shù)μ和控制參數(shù)ε有關的實數(shù).

若映射（17）的雅可比矩陣的特征多項式（19）在分岔點處滿足如下條件：

（T1） 特征值條件

Pμ0（-1）>0，Pu0（1）>0，

Δ±n-1（μ0）=0，Δ±n-3（μ0）=0，

Δ±j（μ0）>0，

（j=n-2，n-4，…，1（or 2））

（T2）橫截條件

2Δ（μ）dμ2μ=μ0≠0，

（T3） 非共振條件

cos（2π/m）≠ψ，ψ=1-Pμ0（1）Δ-32Δ+4，

m=3，4，5，…，

式中：Pμ0（λ），Δ±j（μ0）的具體含義和表達式可參考文獻[12].

則映射在分岔點處發(fā)生非共振的雙Hopf分岔.則可得雙Hopf分岔的顯式臨界條件（具體條件可參考文獻[11]）：

1-m6-m4+m6m4+m1m5+m2m6-

m25-m26-m6+m21-m2m26+m6m1m5=0，（20a）

1-m2m1-m3m2-m4m3-m5m4-m6

-m31-m4m1-m5m2-m6m3

-m4-m51-m6m1m2

-m5-m601m1

-m60001=0.（20b）

3 振動篩系統(tǒng)的數(shù)值模擬

3.1 原系統(tǒng)的數(shù)值模擬

選取振動篩系統(tǒng)的一組參數(shù)為：

m=0.18，c=k=0.182，ζ=0.01，

R=064，l=0784，γ=0.16.

使用模態(tài)疊加法解耦得到系統(tǒng)的解析解以后直接編程計算，可得質(zhì)塊M1的速度1隨著ω變化的分岔圖如圖2所示：

由圖2可以看到系統(tǒng)隨著ω的減小發(fā)生了分岔，通過MATLAB求得精確分岔點是ω=4.911 1，這里探究分岔點ω=4.911 1附近的動力學行為.endprint

通過計算可以得到該系統(tǒng)線性化矩陣的六個特征值分別為：

λ1，2=0.530 0±0.840 9i，λ3，4=0.031 9±

0.099 94i，

λ5，6=-0.465 8±0.418 3i，

其中|λ1，2|=0.994 0，|λ3，4|=1.000 0，|λ5，6|=0626 0，

可以看出有一對復共軛特征值位于單位圓上，系統(tǒng)發(fā)生Hopf分岔，形成擬周期碰振運動.在分岔點ω=4.911 1處原系統(tǒng)的相平面圖如圖3所示，對應的映射圖如圖4所示.

當ω=4.9111時，映射圖為一吸引不變?nèi)?，通過觀察這兩個圖可以看出，質(zhì)塊M1最終趨于穩(wěn)定的擬周期運動.

3.2 受控系統(tǒng)雙Hopf分岔的反控制

選取線性反饋增益矩陣為

α=[k1 -k2；-k2 k1]，β=2ζα，

式中：k1和k2是控制參數(shù).則由以上顯式臨界條件可以得到控制參數(shù)分岔圖，如圖5所示.

圖5中被黑色的粗實線圍成的白色區(qū)域里的點可以滿足條件中的每個不等式，白色區(qū)域里的兩條細黑實線是由（20a）和（20b）兩個等式得到的，通過MATLAB計算出兩條直線的交點，取其中一個交點為（0.835 9，0.668 67），在交點附近取一對控制參數(shù)（0.937，0.758 76），在該控制增益參數(shù)下可得列圖6.

可以看出在該控制參數(shù)下，映射圖6呈現(xiàn)出了“輪胎型”概周期吸引子，得到了原系統(tǒng)的雙Hopf分岔解.

以上對分岔點處的動力學進行了分析和分岔反控制，接下來看非分岔點處的分岔反控制效果，在其它系統(tǒng)參數(shù)不變的情況下，取ω=4.96同樣可以得到6個特征值：

λ1，2=0.538 4±0.835 5i，λ3，4=0.047 5±

0.998 04i，

λ5，6=-0.466 0±0.419 1i，

其中|λ1，2|=0.994 0，|λ3，4|=0.999 1，|λ5，6|=0626 8，可以看出特征值都在單位圓內(nèi)，原系統(tǒng)處于穩(wěn)定的周期運動，如圖7和圖8所示.

同樣可得控制參數(shù)分岔圖，如圖9所示.

通過MATLAB計算出兩個等式（20a）和（20b）的交點，在其中一個交點附近取一對控制參數(shù)（0876 3，0.698 8），在該控制增益參數(shù)下可得到圖10.

由圖可以看出，無論是在分岔點處還是非分岔點處，都可以得到原系統(tǒng)的雙Hopf分岔解，同時也體現(xiàn)了顯式分岔臨界準則的優(yōu)越性和普適性.

4 結(jié) 論

本文建立了振動篩的力學模型，得到分段系統(tǒng)的周期解表達式.通過建立Poincaré 映射，得到映射在不動點處的線性化矩陣，通過分析可知該系統(tǒng)發(fā)生了Hopf分岔.然后應用分岔的顯式臨界準則對原系統(tǒng)進行了雙Hopf分岔反控制的分析，通過調(diào)節(jié)控制增益得到了雙Hopf分岔解.該方法的優(yōu)點就是可以不需要直接計算系統(tǒng)雅克比矩陣的特征值.這一理論結(jié)果有助于今后對多自由度振動篩系統(tǒng)的優(yōu)化設計，也有利于提高振動篩系統(tǒng)的篩分效率，具有一定的實際意義.

參考文獻

[1] 呂高常.大型高頻振動篩動態(tài)特性分析[J].礦山機械，2016，44（6）：58-61.

L Gaochang.Dynamic characteristics analysis on large high-frequency vibrating screen[J].Mining & Processing Equi-pment，2016，44（6）：58-61.（In Chinese）

[2] 趙環(huán)帥.大型高頻振動篩動態(tài)特性及內(nèi)部應力分布特點分析[J].礦山研究與開發(fā)，2013，33（5）：85-89.

ZHAO Huanshuai.Analysis on dynamic behavior and internal stress distribution of large high frequency vibrating screen[J].Mining Research & Development，2013，33（5）：85-89.（In Chinese）

[3] 張則榮，樊智敏，王永巖.大型直線振動篩的疲勞壽命及可靠性分析[J].煤炭學報，2014，39（S1）：257-261.

ZHANG Zerong，F(xiàn)AN Zhimin，WANG Yongyan.Fatigue life and reliability analysis of large linear vibrating screen[J].Journal of China Coal Society，2014，39（S1）：257-261.（In Chinese）

[4] 孫旖.大型振動篩的可靠性研究與應用[J] .礦山機械，2016，44（11）：54-58.

SUN Yi.Reliability study and application of large-scale vibrating screen[J].Mining & Processing Equipment， 2016，44（11）：54-58.（In Chinese）

[5] 陳會征，陳予恕.振動篩分機械相關問題的研究進展[C]//中國力學學會.第十三屆全國非線性振動暨第十屆全國非線性動力學和運動穩(wěn)定性學術會議摘要集.中國天津，2011：870-876.

CHEN Huizheng，CHEN Yushu.Research progress of vibration screening machinery related problems[C]// Chinese Society of Mechanics.Thirteenth National Symposium on Nonlinear Vibration and the 10th National Conference on Nonlinear Dynamics and Motion Stability.Tianjing，China，2011：870-876.（In Chinese）endprint

[6] 白亮亮，李萬祥，沙世民.直線篩分機系統(tǒng)的非線性響應分析[J] .科學技術與工程，2015，15（35）：36-41.

BAI Liangliang，LI Wanxiang，SHA Shimin.Nonlinear response analysis of linear screening system[J].Science Technology and Engineering，2015，15（35）：36-41.（In Chinese）

[7] 張永祥，孔貴芹，俞建寧.振動篩系統(tǒng)的兩類余維三分岔與非常規(guī)混沌演化[J].物理學報，2008，57（10）：6182-6187.

ZHANG Yongxiang，KONG Guiqin，YU Jianning.Two codimension -3 bifurcation and non-typical routes to chaos of a shaker system[J].Acta Physica Sinica，2008，57（10）：6182-6187.（In Chinese）

[8] CHEN Z，YU P.Controlling and anti-controlling Hopf bifurcations in discrete maps using polynomial functions[J].Chaos Soltions & Fractals，2005，26（4）：1231-1248.

[9] 劉素華，唐駕時.四維Qi系統(tǒng)零平衡點的Hopf 分岔反控制[J].物理學報，2008，57（10）：6162-6168.

LIU Suhua，TANG Jiashi.Anti-control of Hopf bifurcation at zero equilibrium of 4D Qi system[J].Acta Physica Sinica，2008，57（10）：6162-6168.（In Chinese）

[10]CHENG Z S.Anti-control of Hopf bifurcation for Chens system through washout filters[J].Neurocomputing，2010，73（16）：3139-3146.

[11]伍新，文桂林，徐慧東，等.慣性式?jīng)_擊振動落砂機周期倍化分岔的反控制[J].固體力學學報，2015，36（1）：28-34.

WU Xin，WEN Guilin，XU Huidong，et al.Anticontrolling period-doubling bifurcation of an inertial impact shaker system[J].Chinese Journal of Solid Mechanics，2015，36（1）.：28-34.（In Chinese）

[12]伍新，文桂林，何麗萍，等.振動落砂機系統(tǒng)的擬周期碰撞設計[J].湖南大學學報：自然科學版，2016，43（4）：38-43.

WU Xin，WEN Guilin，HE Liping，et al.Design of quasi-periodic impact motion of an impact shaker system[J].Journal of Hunan University： Natural Sciences，2016，43（4）：38-43.（In Chinese）

[13]徐慧東，文桂林，伍新，等.三自由度含間隙碰撞振動系統(tǒng)Poincaré 映射 Hopf-Hopf交互分岔的反控制[J].振動工程學報，2015，28（6）：952-959.

XU Huidong，WEN Guilin，WU Xin，et al.Anti-controlling Hopf-Hopf interaction bifurcation on Poincaré map of a three-degree-of-fredom vibro-impact system with clearance[J].Journal of Vibration Enginering，2015，28（6）：952-959.（In Chinese）

[14]WEN G L，CHEN S J，JIM Q T.A new criterion of period-doubling bifurcation in maps and its application to an inertial impact shaker[J].Journal of Sound and Vibration，2008，311（1）：212-223.

[15]伍新，文桂林，徐慧東，等.三自由度含間隙碰撞振動系統(tǒng)Neimark-Sacker分岔的反控制[J].物理學報，2015，64（20）：93-100.

WU Xin，WEN Guilin，XU Huidong，et al.Anti-controlling Neimark-Sacker bifurcation of a three-degree-of-fredom vibration system with clearan-ce[J].Acta Physica Sinica，2015，64（20）：93-100.（In Chinese）endprint