馮長(zhǎng)先， 單小磊， 王志謙

(1．蘭州交通大學(xué) 機(jī)電工程學(xué)院，甘肅 蘭州 730070;2．蘭州交通大學(xué) 數(shù)理與軟件工程學(xué)院，甘肅 蘭州 730070)

齒輪高速轉(zhuǎn)動(dòng)時(shí)，由于齒側(cè)間隙、時(shí)變嚙合剛度、內(nèi)外部激勵(lì)以及安裝誤差等因素的存在，使輪齒間產(chǎn)生嚙合、脫和、再嚙合重復(fù)沖擊的強(qiáng)非線性振動(dòng)，并在一定的參數(shù)條件下產(chǎn)生混沌運(yùn)動(dòng)，對(duì)齒輪系統(tǒng)的動(dòng)態(tài)特性有著不良影響［1］。改善齒輪系統(tǒng)的動(dòng)態(tài)特性的有效方法之一是抑制混沌運(yùn)動(dòng)，使系統(tǒng)運(yùn)動(dòng)回到周期的有序狀態(tài)。目前混沌控制已成為學(xué)者研究的一個(gè)熱點(diǎn)［2-7］。文獻(xiàn)［2］～文獻(xiàn)［4］分別將彈簧剛度反饋法、時(shí)滯控制法和耦合反饋控制法應(yīng)用到一類單自由度振動(dòng)碰撞系統(tǒng)，完成了將系統(tǒng)的混沌運(yùn)動(dòng)抑制到周期軌道的控制目標(biāo)。文獻(xiàn)［5］分析了一類兩自由度碰撞振動(dòng)系統(tǒng)，通過(guò)調(diào)節(jié)外加正弦驅(qū)動(dòng)力，將系統(tǒng)混沌運(yùn)動(dòng)穩(wěn)定到周期軌道。文獻(xiàn)［6］運(yùn)用OGY 控制原理，以混沌吸引子內(nèi)部不穩(wěn)定周期軌道為控制目標(biāo)，通過(guò)對(duì)系統(tǒng)外部激勵(lì)參數(shù)連續(xù)施加小擾動(dòng)信號(hào)，實(shí)現(xiàn)了系統(tǒng)混沌運(yùn)動(dòng)的穩(wěn)定化，由于該方法需要精確的計(jì)算出不動(dòng)點(diǎn)的位置以及目標(biāo)周期軌道在Poincar˙e 截面上的穿越點(diǎn)，故很難應(yīng)用于高維非線性系統(tǒng)。

因此，在文獻(xiàn)［7］自適應(yīng)脈沖微擾法的基礎(chǔ)上，提出了根據(jù)系統(tǒng)運(yùn)動(dòng)與自適應(yīng)速度k 的變化規(guī)律圖，有目的的選取k 值，從而將系統(tǒng)的混沌運(yùn)動(dòng)控制到不同的周期軌道。將該方法運(yùn)用到間隙非線性齒輪系統(tǒng)，并且對(duì)加入高斯白噪聲后的系統(tǒng)進(jìn)行了控制，數(shù)值模擬結(jié)果表明該方法設(shè)計(jì)的控制器結(jié)構(gòu)簡(jiǎn)單、反應(yīng)靈敏、具有較強(qiáng)的魯棒性。

1 力學(xué)模型及其動(dòng)力學(xué)方程分析

圖1 中T1，T2分別表示系統(tǒng)的輸入轉(zhuǎn)矩和負(fù)載轉(zhuǎn)矩;θ1，θ2分別表示主動(dòng)輪和從動(dòng)輪的扭轉(zhuǎn)角位移;I1，I2分別為主動(dòng)輪和從動(dòng)輪的轉(zhuǎn)動(dòng)慣量; rb1，rb2分別為主、從動(dòng)輪的基圓半徑;c 為輪齒嚙合的阻尼系數(shù);2δ 為沿嚙合線上的齒側(cè)間隙;e( τ) 為沿齒輪基圓切向度量的綜合誤差，e( τ) = eacos( wt +φc) ，ea為誤差激勵(lì)幅值，φc為誤差初相位角。

圖1 單級(jí)齒輪系統(tǒng)間隙非線性動(dòng)力學(xué)模型

考慮時(shí)變嚙合剛度得到該模型量綱一化后的微分方程

將其化為狀態(tài)空間的微分方程組如下

其中

式中，x 為量綱一化位移;ζ 為量綱一化阻尼;k0為量綱一化時(shí)變嚙合剛度幅值;Fm為量綱一化內(nèi)部平均載荷;Fah為量綱一化外部平均載荷;Ω 為量綱一化激勵(lì)頻率。

2 脈沖自適應(yīng)控制

設(shè)n 維非線性系統(tǒng)為

式中，x = ( x1，x2，x3，…，xn) 為系統(tǒng)變量;f 為關(guān)于x 的非線性函數(shù);a = ( a1，a2，…，am) 為系統(tǒng)參變量。采用系統(tǒng)參變量脈沖控制混沌系統(tǒng)時(shí)，其控制算法可概述如下:

選擇系統(tǒng)某一參變量aλ(1 ≤λ ≤m) 作為微擾脈沖對(duì)象。從t0時(shí)刻開(kāi)始，每隔Δt 時(shí)刻將a′加入到系統(tǒng)中，a′ 正比于參變量aλ的形式為

式中，h = 0，1，2，…，式中具有兩個(gè)可控參量:脈沖作用時(shí)間間隔Δt 和脈沖強(qiáng)度r。當(dāng)Δt 和r 取得適當(dāng)值時(shí)，就可以把混沌運(yùn)動(dòng)鎮(zhèn)定到周期軌道或不動(dòng)點(diǎn)。其中，脈沖強(qiáng)度r 的選擇是實(shí)現(xiàn)控制的主要環(huán)節(jié)。

為此，設(shè)計(jì)自適應(yīng)控制器˙r = kxij( i ∈［1，n］，j = 0，1，2，…) 來(lái)選取適當(dāng)?shù)拿}沖強(qiáng)度r。系統(tǒng)參變量xi為自適應(yīng)控制器的輸入變量，xi可以是x 的任意分量，k 為自適應(yīng)速度。自適應(yīng)控制器的表達(dá)式最終可表示為

式中，E 可為任意常數(shù)，j +1 為自適應(yīng)控制器的階。j 在理論上可為無(wú)限大，但在實(shí)際應(yīng)用時(shí)為了實(shí)施方便，既不失一般性，又避免選擇的參數(shù)過(guò)多，可取E = 0，且j 為低階。

3 數(shù)值模擬

取參數(shù)齒寬b =25 mm，間隙δ =0．01 mm，齒數(shù)z1=25，z2=25，模數(shù)m =3 mm。計(jì)算參數(shù)k =0．1，Ω =1．0，F(xiàn)m=0．1，F(xiàn)ah=0．2，初值矢量為(0，0)［8］。

系統(tǒng)隨阻尼比ζ 從0．05 至0．09 變化的分岔圖和最大lyapunov 指數(shù)圖如圖2 和圖3 所示。

當(dāng)ζ =0．055 時(shí)，最大lyapunov 指數(shù)λmax=0．035 3，在此參數(shù)條件下系統(tǒng)為混動(dòng)運(yùn)動(dòng)狀態(tài)，圖4、圖5分別為x1的相圖和時(shí)間歷程圖。

選取系統(tǒng)(2) 式中的ζ 作為微擾對(duì)象，以x2作為控制器的輸入變量。則自適應(yīng)控制器的形式為:r =kx23/3( j = 2，E = 0) 。從t0時(shí)刻開(kāi)始，在t0+ hΔt 時(shí)刻系統(tǒng)方程為

其余時(shí)間仍為(2) 式。

在使用自適應(yīng)脈沖控制方法對(duì)混沌運(yùn)動(dòng)進(jìn)行控制時(shí)，控制的另一個(gè)關(guān)鍵環(huán)節(jié)是如何選擇自適應(yīng)速度k，不同的k 將得到不同的周期軌道控制目標(biāo)。但若隨意選取k 后，其控制后的目標(biāo)信息無(wú)法預(yù)測(cè)，更有可能無(wú)法達(dá)到控制的目的。因此，為了確切的對(duì)目標(biāo)進(jìn)行控制，避免對(duì)k 值的毫無(wú)目的的試探，通過(guò)對(duì)受控系統(tǒng)在Poincar˙e 截面上的信息來(lái)確定控制后的目標(biāo)結(jié)果。圖6 是通過(guò)Poincar˙e 截面法獲得的受控系統(tǒng)在ζ = 0．055，其它數(shù)據(jù)同上的變化規(guī)律圖。

圖2 系統(tǒng)隨阻尼比ζ 的分岔圖

圖3 最大lyapunov 指數(shù)圖

圖4 ζ=0．055 相圖

圖5 ζ=0．055 時(shí)間歷程圖

采用四階Runge-kutta 法求解，步長(zhǎng)0．01，開(kāi)始控制時(shí)間n =3 000，為檢驗(yàn)控制的穩(wěn)定性，共積分5 000 步，圖7 為h =2 的數(shù)值仿真結(jié)果，圖中數(shù)據(jù)從1 000 步開(kāi)始截取。

通過(guò)對(duì)圖6 的分析，可以很方便的選取k 值來(lái)獲得目標(biāo)軌道的周期數(shù)，由控制效果圖和相圖可以看出系統(tǒng)的混沌運(yùn)動(dòng)迅速的轉(zhuǎn)化為穩(wěn)定的周期運(yùn)動(dòng)。圖7 中( a)( b) ，( c) ( d) ，( e) ( f) 分別將混沌運(yùn)動(dòng)鎮(zhèn)定到周期二軌道、周期四軌道和周期八軌道的控制效果圖和相圖。經(jīng)過(guò)試驗(yàn)還發(fā)現(xiàn)改變j 仍可以獲得不同的周期軌道，而且當(dāng)j 逐漸增大，其他參數(shù)不變時(shí)，在一定程度上可以加快混沌運(yùn)動(dòng)穩(wěn)定到不動(dòng)點(diǎn)或周期軌道的速度。

對(duì)于實(shí)際系統(tǒng)還應(yīng)考慮安裝誤差，齒形誤差，及其它外界影響，假設(shè)系統(tǒng)(2) 的變量x2受到幅值為σ 的正態(tài)分布的高斯白噪聲，即

式中，N 為標(biāo)準(zhǔn)正態(tài)分布的高斯白噪聲序列，其均值為0，方差為1。取σ =0．01，k =0．01，系統(tǒng)在n =3 000時(shí)進(jìn)行控制，圖中數(shù)據(jù)從n =1 000 開(kāi)始截取，仿真結(jié)果如圖8 所示。

圖6 系統(tǒng)隨自適應(yīng)速度k 的變化規(guī)律圖

圖7 選取不同的k 值將混沌運(yùn)動(dòng)控制到不同的周期軌道

將圖8( a) ( b) 和圖7( a) ( b) 比較可知，在系統(tǒng)受到噪聲干擾與未受干擾時(shí)的控制效果基本一致，這說(shuō)明使用該方法控制混沌時(shí)能夠有效抑制噪聲干擾，具有一定的魯棒性。

另分別選取系統(tǒng)(2) 式中的pm和pah作為微擾對(duì)象，仍以x2作為控制器的輸入變量，其控制效果圖為圖9。

圖8 加入噪聲背景后的控制效果圖和相圖

圖9 控制效果圖

圖9( a) 中，k = -0． 1，σ = 0．01，共疊代9 000 步，數(shù)值從第1 000 步開(kāi)始截取，從第3 000 步開(kāi)始加入控制，5 000 ～7 000 為控制關(guān)閉，第7 000 步開(kāi)始實(shí)施控制并加入高斯白噪聲;圖9( b) 中，k = -0． 07，σ = 0．01，共疊代10 000 步，數(shù)值從第5 000 步開(kāi)始截取，從第6 000 步開(kāi)始加入控制，6 000 ～7 000 為控制關(guān)閉，第8 500 步加入高斯白噪聲并重新實(shí)施控制。從圖中可以看出以pm和pah作為微擾對(duì)象，選取不同的k 值時(shí)，仍可以將系統(tǒng)的混沌運(yùn)動(dòng)鎮(zhèn)定到周期軌道，并在加入噪聲背景后具有較好的控制效果。

4 結(jié)論

基于間隙齒輪系統(tǒng)非線性動(dòng)力學(xué)模型，采用自適應(yīng)脈沖控制法，成功的將其混沌運(yùn)動(dòng)穩(wěn)定到周期軌道，對(duì)實(shí)際的齒輪傳動(dòng)設(shè)計(jì)有一定的理論指導(dǎo)意義。該控制方法的主要原理是每隔Δt 對(duì)控制參量施加一次脈沖信號(hào)，數(shù)值模擬表明了該方法的有效性和較好的魯棒性，而且該方法還具有如下優(yōu)點(diǎn):

(1) 該方法不必預(yù)先求出系統(tǒng)的不動(dòng)點(diǎn)和周期軌道的詳細(xì)信息，也不必等到系統(tǒng)運(yùn)行到控制目標(biāo)附近時(shí)才施加控制，增加了可操控性。

(2) 該方法可以根據(jù)系統(tǒng)運(yùn)動(dòng)隨自適應(yīng)速度k 的變化規(guī)律，來(lái)選擇預(yù)期的控制目標(biāo)軌道，克服了只能用試錯(cuò)法來(lái)不斷調(diào)整控制參量的缺陷。

(3) 可以通過(guò)改變控制器的參數(shù)來(lái)適當(dāng)?shù)母淖兛刂祈憫?yīng)的速度，適用于實(shí)時(shí)在線遞推處理，從而在工程中有一定的應(yīng)用價(jià)值。

［1］李潤(rùn)方，王建軍． 齒輪系統(tǒng)動(dòng)力學(xué)——振動(dòng)、沖擊、噪聲［M］．北京: 科學(xué)出版社，1997．

［2］Silvio，Ibere，Ricardo ． Damping control law for a chaotic impact oscillator［J］． Chaos，Solition and Fractals，2007，32:745-750．

［3］Edward Ott，Celso Grebogi，James A．Yorke controlling chaos［J］．Physical Review Letters，1990，64(11) :1196-1199．

［4］茍向峰，羅冠煒，呂小紅． 含雙側(cè)約束碰撞振動(dòng)系統(tǒng)的混沌控制［J］．機(jī)械科學(xué)與技術(shù)，2011，30(8) :1262-1266．

［5］王艷，張曉娟，王帥．兩自由度碰撞振動(dòng)系統(tǒng)的混沌運(yùn)動(dòng)及控制［J］．石家莊鐵道學(xué)院學(xué)報(bào)，2008，21(4) :63-66．

［6］陳學(xué)森，董海軍，劉曉寧．含時(shí)變嚙合剛度的間隙非線性齒輪系統(tǒng)的混沌控制［J］．機(jī)械科學(xué)與技術(shù)，2006，25(9) :1035-1037．

［7］龔禮華． 基于自適應(yīng)脈沖微擾實(shí)現(xiàn)混沌控制的研究［J］．物理學(xué)報(bào)，2005，54(8) :3502-3506．

［8］郜志英，沈允文，董海軍． 齒輪系統(tǒng)倍周期分岔和混沌層次結(jié)構(gòu)的研究［J］．機(jī)械工程學(xué)報(bào)，2005，41(4) :44-48．