牛寶良

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

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

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

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

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

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

QZS隔振器的垂直彈簧設(shè)計(jì)與普通隔振器相同，根據(jù)已知負(fù)載質(zhì)量、隔振器的彈簧剛度，允許的最大變形量確定垂直彈簧的參數(shù)，然后設(shè)計(jì)水平彈簧的尺寸、剛度、預(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ù)剛度、組合后的零剛度曲線見(jiàn)圖2。在Y向±10 mm范圍內(nèi)，達(dá)到準(zhǔn)零剛度。

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

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

2.1 隔力模型

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

2.2 隔幅模型

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

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

圖4 QZS彈簧仿真模型

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

2.3 沖擊隔離模型

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

圖7 沖擊激勵(lì)波形

3 仿真分析

3.1 隔力仿真

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

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

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

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

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

3.2 隔幅仿真

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

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

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

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

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

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

3.3 隔沖擊仿真

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

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

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

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

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

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

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

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

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

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

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

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

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

4 結(jié)束語(yǔ)

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

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

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

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

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

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

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

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

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

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

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

[7] 徐道臨，張?jiān)掠ⅲ芗酉?，?一種準(zhǔn)零剛度隔振器的特性分析與實(shí)驗(yàn)研究[J].振動(dòng)與沖擊，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].長(zhǎng)沙:湖南大學(xué)，2015.

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

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

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

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

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

[17] 柴凱，楊慶超，朱石堅(jiān)，等.矩形永磁鐵的磁力計(jì)算及其應(yīng)用[J].噪聲與振動(dòng)控制，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.