準(zhǔn)零剛度隔振器的時域仿真

牛寶良

(中國工程物理研究院 總體工程研究所， 四川 綿陽 621999)

普通隔振器總存在一個共振頻率，在共振頻率及其附近起不到隔振作用，反而有放大作用，這個問題一直困擾著振動工程界。1980年初Alabuzhev等[1]基于正負(fù)剛度并聯(lián)原理開發(fā)了一些準(zhǔn)零剛度(QZC)隔振器，并簡單驗(yàn)證了機(jī)構(gòu)的隔振特性。與未使用負(fù)剛度的機(jī)構(gòu)相比，準(zhǔn)零剛度隔振機(jī)構(gòu)具有顯著的低頻隔振效果。之后，國內(nèi)有學(xué)者一直在研究準(zhǔn)零剛度隔振器。湖南大學(xué)彭獻(xiàn)等[2-6]從20世紀(jì)80年代后期開始了相關(guān)研究；徐道臨團(tuán)隊[7-14]近年來做了大量的研究，提出了幾種線振動的QZS隔振器構(gòu)型及相應(yīng)的分析方法，并進(jìn)行了試驗(yàn)研究，也提出了扭振動的QZS構(gòu)型和分析；南京航空航天大學(xué)王勇等[15]在QZS隔振器的應(yīng)用方面進(jìn)行了探索，包括車-座椅-人耦合模型的響應(yīng)分析[15]、含有QZS的有源無源混合隔振控制[16]。海軍工程大學(xué)柴凱等[17]提出了基于磁性負(fù)剛度的QZS隔振器及其設(shè)計方法。國外也有相關(guān)的研究報道，Kyoung-Rock Kim[18]研究了寬負(fù)載范圍的QZS設(shè)計；Pham Van Trung[19]研究了有源無源QZS隔振器，消除了無源QZS遇到的放大現(xiàn)象；Ali Abolfathi[20]研究了QZS參數(shù)失調(diào)情況下的隔振效果。總之，作為一種低頻隔振性能優(yōu)越的隔振器，準(zhǔn)零剛度隔振器在構(gòu)型設(shè)計、性能分析等方面正處于持續(xù)研究中，相關(guān)的研究成果多是基于頻域的。準(zhǔn)零剛度隔振器含有顯著的非線性特性，時域仿真能準(zhǔn)確地描述其非線性特性，其仿真結(jié)果也更加直觀明了。因此，本文對準(zhǔn)零剛度隔振器進(jìn)行建模和時域仿真。

1 準(zhǔn)零剛度隔振器

圖1 QZS隔振器結(jié)構(gòu)示意圖

圖2 正剛度、負(fù)剛度、QZS剛度力-位移曲線

準(zhǔn)零剛度隔振器一般的構(gòu)成辦法是正剛度彈簧、負(fù)剛度彈簧并聯(lián)，合理搭配正負(fù)剛度參數(shù)，使得在零位附近一個區(qū)間剛度為0或者接近0。最典型的代表是Carrella等[7]將3個線性或者非線性彈簧進(jìn)行組合，通過參數(shù)優(yōu)化得到的一種QZS系統(tǒng)。本文所要仿真的也是這種準(zhǔn)零剛度隔振器。一個垂直向的彈簧支承負(fù)載質(zhì)量，是正彈簧剛度，兩個水平向?qū)數(shù)膹椈蓸?gòu)成豎直向的負(fù)剛度彈簧(在零位附近的一個區(qū)間)。QZC隔振器結(jié)構(gòu)示意圖如圖1所示。

QZS隔振器的垂直彈簧設(shè)計與普通隔振器相同，根據(jù)已知負(fù)載質(zhì)量、隔振器的彈簧剛度，允許的最大變形量確定垂直彈簧的參數(shù)，然后設(shè)計水平彈簧的尺寸、剛度、預(yù)壓縮量，使得在一定范圍內(nèi)綜合彈簧剛度基本為0。不失一般性，給定一組參數(shù)，質(zhì)量m=10 kg，豎向彈簧剛度ky=19 613 N/m，阻尼系數(shù)為44.287 N/(m/s)(阻尼比為0.05)，橫向彈簧剛度kx1=kx2=15ky。它的正剛度、負(fù)剛度、組合后的零剛度曲線見圖2。在Y向±10 mm范圍內(nèi)，達(dá)到準(zhǔn)零剛度。

2 準(zhǔn)零剛度隔振器建模

為通過仿真了解QZS隔振器的隔振性能，建立了基于Matlab/Simulink的仿真模型。

2.1 隔力模型

這類隔振系統(tǒng)希望傳遞到基礎(chǔ)的力更小，建立的仿真模型如圖3所示，其中QZS彈簧仿真模型如圖4所示。

2.2 隔幅模型

這類隔振系統(tǒng)希望基礎(chǔ)振動傳遞到質(zhì)量m上的運(yùn)動更小，激勵是基礎(chǔ)的運(yùn)動。建立的仿真模型如圖5所示，其中QZS彈簧仿真模型如圖6所示。

圖3 QZS隔振器(隔力)仿真模型

圖4 QZS彈簧仿真模型

圖5 QZS隔振器(隔幅)仿真模型

2.3 沖擊隔離模型

這類隔振系統(tǒng)希望基礎(chǔ)沖擊傳遞到質(zhì)量m上的運(yùn)動更小。激勵是基礎(chǔ)的沖擊運(yùn)動，建立的仿真模型與圖5類似，區(qū)別在于激勵波形。這里給出的激勵是一個由半正弦3次方構(gòu)造的一個沖擊位移波形，它具有速度零起始零結(jié)束、加速度零起始零結(jié)束的特性，可以模擬路面凸起。它的加速度、位移、速度波形如圖7所示。仿真模型與圖5、6基本一致，區(qū)別只是基礎(chǔ)位移的生成。它的參數(shù)只有2個：位移幅值、脈沖寬度。

圖7 沖擊激勵波形

3 仿真分析

3.1 隔力仿真

低于這個臨界狀態(tài)，比如位移為3 mm時，則呈現(xiàn)非常好的衰減，力比最大值為0.65，而且隨著頻率增加呈現(xiàn)單調(diào)下降的趨勢，如圖9所示。

高于這個臨界狀態(tài)，比如預(yù)期位移為5 mm時，它低頻段呈現(xiàn)嚴(yán)重的惡化，2.6 Hz以下呈現(xiàn)放大狀態(tài)，最大放大5倍，大于2.6 Hz以后，進(jìn)入良好衰減狀態(tài)，如圖10所示。

局部的時域波形如圖11所示。在放大階段，位移達(dá)到±15 mm時，加速度呈現(xiàn)非常明顯的非正弦特征。究其原因，就是位移太大，進(jìn)入了彈簧的非零區(qū)，而且這個非零區(qū)由于負(fù)剛度彈簧的非線性使得綜合剛度呈現(xiàn)非線性，因此響應(yīng)加速度呈現(xiàn)非線性，從而導(dǎo)致波形失真。

圖8 位移達(dá)到±4 mm時的力比曲線

圖10 位移達(dá)到±15 mm時的力比曲線

3.2 隔幅仿真

隔幅，就是希望基礎(chǔ)激勵時質(zhì)量m的運(yùn)動幅值小。這里用m的振動位移/基礎(chǔ)的激勵位移(以下簡稱位移比)來衡量。

通過仿真發(fā)現(xiàn)：對于所給定的這組參數(shù)，當(dāng)振動位移達(dá)到4 mm時，傳遞到基礎(chǔ)的力與激勵力之比略超1，如圖12所示。這是一個臨界狀態(tài)。

低于這個臨界狀態(tài)，比如位移為3 mm時，則呈現(xiàn)非常好的衰減，力比最大值為0.65，而且隨著頻率增加呈單調(diào)下降的趨勢，如圖13所示。

高于這個臨界狀態(tài)，比如基礎(chǔ)振動位移為5 mm時，其低頻段呈現(xiàn)嚴(yán)重的惡化。質(zhì)量m的振動位移達(dá)到17.526 mm，位移比曲線如圖14所示，2.6 Hz以下呈現(xiàn)放大狀態(tài)，最大放大3.8倍。從放大轉(zhuǎn)衰減的過渡過程的時域曲線見圖15。加速度也嚴(yán)重失真，但是與隔力時的加速度波形失真有所不同。

圖12 基礎(chǔ)位移達(dá)到±4 mm時的位移比曲線

圖14 基礎(chǔ)位移達(dá)到±5 mm時的位移比曲線

3.3 隔沖擊仿真

隔沖擊，就是希望在基礎(chǔ)激勵時質(zhì)量m的運(yùn)動幅值小。這里用m的振動位移/基礎(chǔ)的激勵位移(以下簡稱位移比)來衡量。同時也給出速度比、加速度比，以方便比較。

后來我一直沒去景花廠。阿花邀請過我，我就說忙。忙不過是個借口，實(shí)際上是我不想去。我是個老實(shí)人，不喜歡趨炎附勢那一套。雖然我對美女從不敬而遠(yuǎn)之，但當(dāng)美女和老板合二為一時，我會牢牢抓住心猿意馬的韁繩，敬而遠(yuǎn)之了。還有，大家都是做拋光的，總往哪兒跑怕別人會有想法，所以，不去為好。

如前所述，基礎(chǔ)沖擊激勵采用沖擊位移，它的參數(shù)是幅值、脈寬。以隔振器的共振周期作為參考。

當(dāng)激勵位移取10 mm、激勵脈寬取隔振器共振周期(不考慮負(fù)剛度時的普通隔振器)時，普通隔振器質(zhì)量m的振動響應(yīng)如圖16所示。

QZS隔振器的響應(yīng)如圖17所示。響應(yīng)加速度峰值比為0.089 839，速度峰值比為0.226 85，位移峰值比為0.428 6。由圖17可見呈現(xiàn)良好的隔振效果。

當(dāng)激勵脈寬取隔振器共振周期的2倍時，響應(yīng)加速度峰值比為0.197 21，速度峰值比為0.454 96，位移峰值比為0.656 44。

當(dāng)激勵脈寬取隔振器共振周期的4倍時，響應(yīng)加速度峰值比為0.633 11，速度峰值比為0.774 13，位移峰值比為0.917 08。

當(dāng)激勵脈寬取隔振器共振周期的8倍時，響應(yīng)加速度峰值比為1.159 8，速度峰值比為1.279 5，位移峰值比為1.050 8。

當(dāng)激勵脈寬取隔振器共振周期的40倍時，響應(yīng)加速度峰值比為0.977 83，速度峰值比為1.089 9，位移峰值比為1.026 1。

當(dāng)激勵脈寬取隔振器共振周期的0.5倍時，響應(yīng)加速度峰值比為0.041 171，速度峰值比為0.103 88，位移峰值比為0.298 74。

圖16 普通隔振器的沖擊響應(yīng)

圖18 QZS隔振器的沖擊響應(yīng)(40×1/fn)

從以上數(shù)據(jù)可以得出：基礎(chǔ)振動為10 mm時，在脈寬4倍共振周期以下，位移、速度、加速度都是衰減的，脈寬越小，衰減越好；在脈寬4倍共振周期以上至8倍共振周期，基本在幅值比為1附近，有輕微放大；在8倍以上，位移比、速度比、加速度比都基本上為1，因?yàn)檫@時頻率很低，質(zhì)量m基本上與基礎(chǔ)同步運(yùn)動。

4 結(jié)束語

本文建立了基于正負(fù)彈簧剛度合成的QZS隔振器的時域仿真模型，開展了隔力、隔幅、隔沖擊仿真。

仿真模型運(yùn)行正確快速，為QZS隔振器的參數(shù)設(shè)計提供了方便。

盡管只對一組QZS隔振器數(shù)據(jù)進(jìn)行了仿真，但仿真結(jié)果揭示了客觀規(guī)律，即：在激勵幅值比較小的情況下，QZS隔振器無論是隔力、隔幅都非常好，全頻段不會出現(xiàn)共振放大現(xiàn)象；激勵幅值超過臨界值以后，則會在低頻段出現(xiàn)放大現(xiàn)象，這提示設(shè)計者在應(yīng)用QZS隔振器時需要對激勵的幅值、頻率范圍加以分析預(yù)判。

對于沖擊基礎(chǔ)激勵，QZS隔振器性能非常好，在激勵位移達(dá)到10 mm(對于本文QZS隔振器，10 mm對應(yīng)零剛度區(qū)域)的情況下，在小脈寬(小于4倍共振周期)時，衰減效果很好，脈寬越小衰減越顯著。40倍共振周期以上，基本不衰減不放大。總之，對于沖擊來說，QZS隔振器是適合的。

本文僅作了參數(shù)設(shè)計與仿真，具體的效果還與QZS隔振器的工程實(shí)現(xiàn)結(jié)構(gòu)相關(guān)，還需進(jìn)一步研究。

[1] 憂高法.新巧準(zhǔn)零剛麼隔振平臺減振機(jī)理及動力學(xué)性態(tài)研究[D].蕪湖:安徽工程大學(xué)，2016.

[2] 陳樹年，黃天擇，彭獻(xiàn).汽車座椅彈性元件剛度的探討[J].湖南大學(xué)學(xué)報(自然科學(xué)版)，1989，16(3):14-20.

[3] 彭獻(xiàn)，陳樹年，宋福磐.負(fù)剛度的工作原理及應(yīng)用初探[J].湖南大學(xué)學(xué)報(自然科學(xué)版)，1992，19(4):89-94.

[4] 彭獻(xiàn)，黎大志，陳樹年.準(zhǔn)零剛度隔振器及其彈性特性設(shè)計[J].振動、測試與診斷，1997，17(4):44-46.

[5] 彭獻(xiàn)，張施詳.一種準(zhǔn)零剛度被動隔振系統(tǒng)的非線性共振響應(yīng)分析[J].湖南大學(xué)學(xué)報(自然科學(xué)版):2011，38(8):34-39.

[6] 彭獻(xiàn)，張施詳.一種準(zhǔn)零剛度被動隔振系統(tǒng)的靜力與線性動力特性分析[J].力學(xué)季刊，2012，33(3):492- 498.

[7] 徐道臨，張月英，周加喜，等.一種準(zhǔn)零剛度隔振器的特性分析與實(shí)驗(yàn)研究[J].振動與沖擊，2014，33(11):208-213.

[8] XU D，ZHANG Y，ZHOU J，et al.On the analytical and experimental assessment of the performance of a quasi-zero-stiffness isolator[J].Journal of Vibration and Control，2014，20(15):2314-2325.

[9] ZHOU J，XU D，BISHOP S.A torsion quasi-zero stiffness vibration isolator[J].Journal of Sound and Vibration，2015，338:121-133.

[10] ZHOU J，WANG X，XU D，et al.Nonlinear dynamic characteristics of a quasi-zero stiffness vibration isolator with cam-roller-spring mechanisms[J].Journal of Sound and Vibration，2015，346:53-69.

[11] 王心龍.CRSM準(zhǔn)零剛度隔振器的非線性隔振特性及實(shí)驗(yàn)研究[D].長沙:湖南大學(xué)，2015.

[12] 周杰.碟形橡膠準(zhǔn)零隔振器的設(shè)計與特性分析[D].長沙:湖南大學(xué)，2015.

[13] 周加喜，徐道臨.單層準(zhǔn)零剛度浮筏隔振特性分析及混沌化[J].中國科技論文，2014，8(8):911-915.

[14] 徐道臨，周杰，周加喜，等.碟形橡膠準(zhǔn)零剛度隔振器的設(shè)計和特性分析[J].湖南大學(xué)學(xué)報(自然科學(xué)版)， 2015，42(8):22-28.

[15] 王勇，李舜酩，程春.基于準(zhǔn)零剛度隔振器的車-座椅-人耦合模型動態(tài)特性研究[J].振動與沖擊，2016，35(15):190-196.

[16] 王勇，李舜酩，程春，等.立方速度反饋控制的準(zhǔn)零剛度隔振器動力學(xué)特性分析[J].振動工程學(xué)報，2016，29(2):305-312.

[17] 柴凱，楊慶超，朱石堅，等.矩形永磁鐵的磁力計算及其應(yīng)用[J].噪聲與振動控制，2016，36(6):51-55.

[18] KIM K R，YOU Y H，AHN H J.Optimal design of a QZS isolator using flexures for a wide range of payload[J].Int J Precis Eng Manuf，2013，14(6):911-917.

[19] VAN T P，KIM K R，AHN H J.A nonlinear control of an QZS isolator with flexures based on a lyapunov function[J].International Journal of Precision Engineering and Manufacturing，2013，14(6):919-924.

[20] ABOLFATHI A，BRENNAN M J，WATERS T P，et al.On the effects of mistuning a force-excited system containing a quasi-zero-stiffness vibration isolator[J].Journal of Vibration and Acoustics，2015，137(4):044502-0044507.