劉樹勇 蘇 攀 位秀雷 刁愛民

(海軍工程大學動力工程學院1) 武漢 430033) (中國人民解放軍91404部隊2) 秦皇島 066000)

0 引 言

如何有效的進行低頻隔振一直以來是振動領域的熱點問題.為了提高隔振器的隔振范圍，可減小隔振器的固有頻率，這將意味著降低系統(tǒng)剛度，從而使得系統(tǒng)的承載能力下降.為了克服這一缺點，諸多學者提出了具有高靜低動剛度性能的隔振器.高靜低動剛度隔振器由正剛度機構(gòu)和負剛度機構(gòu)并聯(lián)組成，其中正剛度機構(gòu)用于承受主要載荷，負剛度機構(gòu)用于抵消正剛度彈性元件的剛度，使得系統(tǒng)在靜態(tài)平衡位置處的剛度趨于零.Ivana等[1-6]將傾斜彈簧作為負剛度機構(gòu)與垂直彈簧并聯(lián)組成高靜低動剛度隔振器，并詳細研究了其幅頻特性及傳遞特性.Xu等[7-9]使用電磁彈簧作為負剛度機構(gòu)，設計了準零剛度隔振器，并進行了理論與實驗分析.Meng等[10]提出來一種具有等厚和變厚度蝶形彈簧準零剛度隔振器的設計，并利用平均法研究了系統(tǒng)參數(shù)對隔振器傳遞率的影響.Huang等[11]利用歐拉屈曲梁作為負剛度機構(gòu)，分析了該高靜低動剛度隔振器的靜力與動力學特性并進行了試驗研究，驗證了高靜低動剛度隔振器比線性隔振器具有更寬的隔振頻帶.

目前對含反饋的多自由度高靜低動剛度隔振系統(tǒng)動力學特征的控制規(guī)律研究較為鮮見.本文應用平均法研究反饋控制對高靜低動剛度振動系統(tǒng)動力學特性的影響，得到系統(tǒng)響應幅值和控制參數(shù)之間的關系式，并對系統(tǒng)在Van der Pol平面上的解進行分析.

1 高靜低動剛度系統(tǒng)力學特性分析

高靜低動剛度系統(tǒng)原理見圖1.兩個負剛度傾斜彈簧與一個正剛度垂直彈簧在O點連接，正剛度彈簧元件剛度為kv，兩負剛度彈簧元件的長度相同，且另一端安裝點的高度相同，安裝點距離系統(tǒng)中心點的距離為a，左側(cè)彈簧剛度為kl，右側(cè)彈簧剛度為kr，系統(tǒng)處于靜平衡位置時，正剛度彈簧元件的壓縮量為h0.

坐標x定義為點O從初始位置開始在豎直方向上的位移，向下為正，則系統(tǒng)的力與位移的關系為

(1)

進行無量綱化，可得：

(2)

圖1 高靜低動系統(tǒng)原理圖

根據(jù)剛度定義，將式(2)對位移求導得準零系統(tǒng)的無量綱剛度為

(3)

(4)

當α=1時，無量綱力-位移及剛度-位移曲線見圖2.由圖2可知，當h=hQZS時，系統(tǒng)在靜平衡位置時由兩側(cè)彈簧產(chǎn)生的負剛度與垂直彈簧產(chǎn)生的正剛度相互抵消，可獲得準零剛度特性；當hhQZS時，系統(tǒng)在靜平衡位置時剛度為負值.

圖2 量綱一的量力-位移和剛度-位移曲線圖

(5)

則剛度為

(6)

圖3 力和剛度的近似表達式與精確表達式對比曲線(α=1,h=hQZS)

高靜低動剛度隔振器無量綱力及剛度精確表達式與近似表達式對比曲線見圖3.由圖3可知，無量綱力及剛度近似表達式與精確表達式的誤差隨著位移量的增大而增大.當系統(tǒng)在平衡位置附近處位移量較小時，近似表達式與精確表達式的誤差很小，三階泰勒級數(shù)展開式能很好地模擬精確表達式.

2 高靜低動剛度系統(tǒng)反饋控制模型

兩自由度高靜低動剛度振動系統(tǒng)反饋控制模型見圖4，其中系統(tǒng)的激振力幅值為F，頻率為Ω；上層質(zhì)量塊為M1，振動位移為X1，阻尼系數(shù)為C1； 下層質(zhì)量塊為M2，振動位移為X2，下層彈簧元件剛度為K2，阻尼系數(shù)為C2.反饋控制回路中，傳感器采集上層質(zhì)量塊速度信號，作動器安裝于上下質(zhì)量塊之間，根據(jù)反饋信號控制上層質(zhì)量塊的振動，因此作動器的力可表示為：P=K1Z1.

圖4 兩自由度高靜低動剛度系統(tǒng)反饋控制模型

系統(tǒng)的微分方程為

Fcos(ΩT)+M1g-P

(7)

系統(tǒng)的量綱一的量運動微分方程為

fcos(ωt)-kz1

(8)

3 平均法分析系統(tǒng)響應

對式(8)而言，其精確解難以得到，為便于研究，采用平均法研究方程的近似解以及反饋增益對解的影響.

假設在線性條件下，式(8)基本解為

(9)

對式(9)求導，可得

(10)

在非線性條件下，由于振動的復雜性，假設可以得到與式(9)在形式上相同，但系數(shù)具有時變特征的基本解，其表達式為

(11)

同時假設系統(tǒng)的振動速度具有式(10)相同的形式：

(12)

由于方程(10)求導后可得到

(13)

聯(lián)立式(12)～(13)，得到

(14)

對式(12)求導，得到二階導數(shù)形式

(15)

將式(12)～(14)代入式(2)中，可得

(1+k-ω2)a1(t)cos(ωt-φ1(t))-ωξ1a1(t)sin(ωt-φ1(t))-fcos(ωt)+

(16)

χ+(v1+μ-ω2)a2(t)cos(ωt-φ2(t))-(ωξ2+μωξ1)a2(t)sin(ωt-φ2(t))+

3rμ(a1(t)cos(ωt-φ1(t))-a2(t)cos(ωt-φ2(t)))a2(t)a1(t)cos(ωt-φ1(t))cos(ωt-φ2(t))+

(17)

利用式(14)、式(15)和式(17)直接求出式(9)中幅值和相位參數(shù)的一階導數(shù)解析表達式.

(18a)

(18b)

(-(v1+μ-ω2)cos(ωt-φ2(t))+(μξ1+ξ2)ωsin(ωt-φ2(t)))a2(t)-

(18c)

(-(v1+μ-ω2)cos(ωt-φ2(t))+(μξ1+ξ2)ωsin(ωt-φ2(t)))a2(t)-

(18d)

為方便進行分析，將非自治系統(tǒng)轉(zhuǎn)換為自治系統(tǒng)，進行如下變換：

φ1(t)=ωt-θ1(t)，φ2(t)=ωt-θ2(t)，ωt=θ1(t)+φ1(t)

則變換后，式18a)～d)變?yōu)?/p>

(19a)

(19b)

(19c)

(19d)

假設參數(shù)a1(t)，a2(t)，φ1(t)，φ2(t)具有慢變特性，對式19a)～d)等式右邊在一個周期T( 0～2π) 范圍內(nèi)積分，并計算該區(qū)間的平均值，可得：

(20a)

[πfcosφ1(t)/ωa1(t)]-

(20b)

(20c)

(20d)

由式20b)可得，反饋參數(shù)k的存在對相位產(chǎn)生了影響，同時由于20b)中含有a1(t)和a2(t)，它們和式20a)和式20c)是非線性耦合的，因此還造成響應幅值的改變.由于相位變化與控制過程中延遲有重要關聯(lián)，因此，這一結(jié)果表明多自由度高靜低動剛度振動控制中增益對延遲具有影響.在施加控制的過程中，相位的補償不僅需要考慮系統(tǒng)本身的延遲，而且考慮控制增益參數(shù)的影響.此外，如果響應的幅值ai(t)i=1，2在t趨向于無窮大為常數(shù)時，則系統(tǒng)的解曲線收斂于一個極限環(huán)，而當ai(t)i=1,2不為常數(shù)時，則系統(tǒng)的解將非常復雜，可能出現(xiàn)擬周期甚至混沌狀態(tài).

4 系統(tǒng)動態(tài)特性分析

當系統(tǒng)無反饋時，式(8)中參數(shù)的取值為ξ1=0.2、ξ2=0.2、μ=1、f=10、ω=3、r=1、v1=1、χ=2.通過式20a)～d)可預測系統(tǒng)的長期行為，見圖5a).此時系統(tǒng)振動幅值a1(t)隨時間無限增大，趨向于常數(shù)3.758.說明系統(tǒng)處于周期運動，系統(tǒng)的運動最終將收斂到極限環(huán)上.當對系統(tǒng)進行反饋控制時，反饋增益取為k=5，其余參數(shù)不變，結(jié)果見圖5b).此時系統(tǒng)振動幅值a1(t)隨時間增大，趨向于常數(shù)2.926.表明系統(tǒng)的運動仍處于周期運動.對比可發(fā)現(xiàn)：通過反饋控制后，極限環(huán)變小.為了觀察二維平面上的收斂行為，給出了Van der pol 平面上的收斂軌線，見圖5c～d).因此，通過對解析表達式的計算分析可以了解非線性系統(tǒng)的長期特性，以及反饋控制后的響應趨向.顯然，通過反饋控制后，能夠降低振動幅值.

圖5 應用近似解析方法得到響應的特征

當反饋增益進一步增大為k=10，其余參數(shù)保持不變，見圖6.由圖 6可知，此時響應的幅值收斂到1.996 ，即隨著控制強度的增加極限環(huán)大小進一步降低，振動得到進一步抑制.

圖6 反饋增益增大時的響應幅值收斂行為

改變系統(tǒng)的參數(shù)，如阻尼減少，激勵頻率增大，觀察結(jié)果的變化.當參數(shù)的取值為ξ1=0.05，ξ2=0.01，μ=1，f=10，ω=3.8，r=1，v1=1，χ=2時進行計算，見圖7.系統(tǒng)振動幅值a1(t)，a2(t)隨時間增大始終在一定范圍內(nèi)波動，表明在不同初始條件下，系統(tǒng)的解并不是在單一的極限環(huán)上，而是在“極限環(huán)帶”上運動，此時系統(tǒng)可能處于多周期運動或混沌運動.

圖7 系統(tǒng)響應形成的極限環(huán)帶

為了說明圖7的預測結(jié)果，改變反饋控制增益大小，得到系統(tǒng)Lyapunov指數(shù)譜曲線隨反饋增益變化的規(guī)律，見圖8. 由圖8可知，反饋增益在0.02～2范圍呈現(xiàn)混沌運動，最大Lyapunov指數(shù)大于零；在2～10范圍內(nèi)，混沌運動和周期運動區(qū)域相互交錯.該結(jié)果與圖7理論預測結(jié)果相吻合.

圖8 隨反饋增益變化的Lyapunov指數(shù)譜

5 結(jié) 束 語

本文建立了含反饋控制的多自由度高靜低動剛度隔振系統(tǒng)的數(shù)學模型，應用平均法推導了多自由度高靜低動剛度隔振系統(tǒng)的響應幅值隨反饋增益變化的解析表達式，并預測了系統(tǒng)響應的久期變化趨勢，得到了Van der Pol二平面上的極限環(huán).當存在“純凈的”極限環(huán)時，系統(tǒng)的呈現(xiàn)簡單的周期振動；當極限環(huán)轉(zhuǎn)化為“環(huán)帶”時，系統(tǒng)的運動變得復雜，存在擬周期和混沌運動狀態(tài).多自由度高靜低動剛度系統(tǒng)中引入反饋控制后，既可以將系統(tǒng)的混沌響應控制到周期狀態(tài)，抑制其有害的混沌振動；也可以通過反控制進入小幅值混沌狀態(tài).