王彤，周瑾，孟凡許，籍文韜

(南京航空航天大學(xué) 機電學(xué)院，江蘇 南京 210016)

0 引言

隨著列車運行速度的不斷提高，列車輪軌之間垂向振動增加，被動減振器由于其性能不可調(diào)已無法滿足減振需求，而主動減振器又具有所需能量大、裝置復(fù)雜的缺點，半主動減振裝置成為近些年列車減振研究的熱點[1]。其基本思想是根據(jù)振動環(huán)境與系統(tǒng)的狀態(tài)調(diào)節(jié)減振器的阻尼特性[2]。對于列車的半主動減振，目前主要有改變阻尼介質(zhì)黏度與改變液壓節(jié)流面積兩種方式。常見的改變阻尼介質(zhì)的減振器有磁流變阻尼器、電流變液阻尼器[3]。對于列車磁流變半主動減振器的研究，目前大多是動力學(xué)建模與仿真的研究[4]，根據(jù)列車減振參數(shù)設(shè)計磁流變阻尼器的結(jié)構(gòu)參數(shù)的相關(guān)研究較少。本文設(shè)計用于某型列車模型垂向減振的剪切閥式磁流變減振器，為根據(jù)列車參數(shù)設(shè)計用于列車垂向磁流變減振器結(jié)構(gòu)參數(shù)的問題提供了一種解決方法與思路，具有很強的工程實用價值。

1 列車磁流變阻尼器設(shè)計方法

磁流變阻尼器的工作模式分為閥式、剪切式、擠壓式與剪切閥式4種，其中剪切閥式磁流變阻尼器因其結(jié)構(gòu)形式簡單、出力大而得到廣泛應(yīng)用。對于磁流變阻尼器，出力大小和出力范圍是衡量磁流變阻尼器性能的重要條件，一個設(shè)計良好的磁流變阻尼器需要滿足出力范圍且具有廣泛的調(diào)節(jié)范圍。

磁流變阻尼器的阻尼力直接影響列車減振效果。過小的阻尼力不能完全控制住列車的振動，阻尼器力過大的阻尼器可調(diào)范圍減小、尺寸大、成本高[5]。如何確定列車磁流變阻尼器阻尼力，并根據(jù)阻尼力快速方便地計算結(jié)構(gòu)參數(shù)是設(shè)計過程中必須考慮的問題。針對此問題，提出一種基于列車整車阻尼比與磁流變阻尼器出力模型的優(yōu)化設(shè)計方法，首先將列車模型簡化為單自由度系統(tǒng)，使用阻尼比公式計算列車模型垂向減振磁流變阻尼器的最大阻尼力。然后以最大阻尼力為優(yōu)化性能約束條件，以使磁流變阻尼器的可調(diào)范圍最大為目標(biāo)值，使用Matlab優(yōu)化工具箱得到磁流變阻尼器的結(jié)構(gòu)參數(shù)。以某型列車模型垂向減振磁流變阻尼器為例，設(shè)計方法原理及過程如下所述。

2 列車垂向減振阻尼器阻尼力計算

所設(shè)計磁流變阻尼器的應(yīng)用對象為某型號列車的按比例縮小的列車模型。列車模型由1節(jié)車廂以及2個轉(zhuǎn)向架構(gòu)成，如圖1。車廂采用鋼制骨架蒙鐵皮，長約6 m。車廂與轉(zhuǎn)向架之間使用筒形空氣彈簧連接，轉(zhuǎn)向架與輪對之間由于彈簧剛度較大看作固連。

圖1 縮比列車模型車體

列車模型具體參數(shù)如表1。

表1 列車模型參數(shù)

將列車模型簡化為單自由度模型， 如圖2。左側(cè)為列車模型的示意圖，由1個車廂與2個轉(zhuǎn)向架組成，車廂質(zhì)量為m2，轉(zhuǎn)向架質(zhì)量為m1，在車廂與1個轉(zhuǎn)向架之間有2個空氣彈簧與二系阻尼器，其剛度系數(shù)與阻尼系數(shù)為k2與c2。

圖2 列車單自由度整車模型

所設(shè)計的磁流變阻尼器用于代替列車二系阻尼器。將整個系統(tǒng)簡化為單自由度系統(tǒng)并作如下假設(shè)：系統(tǒng)沿縱向與橫向方向?qū)ΨQ，僅考慮車身的垂向運動；車身質(zhì)量平均分配給每個空氣彈簧與二系阻尼器，且使用相同型號的空氣彈簧與二系阻尼器；將車廂簡化為質(zhì)量塊，車廂質(zhì)量為md=m2，4個空氣彈簧的總剛度系數(shù)kd=4·k2，4個二系垂向阻尼器的總阻尼系數(shù)cd=4·c2。則列車整車阻尼比為：

(1)

根據(jù)相關(guān)文獻[6-7]，列車的垂向振動的整車阻尼比一般在0.2～0.45之間。使用式(1)計算阻尼比在0.2與0.45時阻尼器的阻尼系數(shù)，結(jié)果如表2第2行。磁流變阻尼器的阻尼力具有很強的非線性以及滯回性，此處將其簡化為線性阻尼，認為阻尼力與速度成正比，且車廂垂向的振動速度范圍一般為0～0.3 m/s，此處按0.3m/s取。那么按照整車共安裝4個磁流變阻尼器，則單個磁流變阻尼器阻尼力大小如表2第3行。

表2 不同阻尼比時單個磁流變阻尼器參數(shù)

由于阻尼比在0.2～0.45間是一般被動阻尼器達到的阻尼比優(yōu)選值，考慮到磁流變阻尼器的阻尼力可由控制器控制其阻尼力，不同的控制策略原理所需最大、最小阻尼力均不同；為了保證可以為控制策略留出足夠的控制余量，將整車阻尼比的范圍擴大到0.65，如表2最右一列所示，故所設(shè)計的磁流變阻尼器的最大阻尼系數(shù)為7 666 N·s/m，最大阻尼力為2 300 N。

3 二系阻尼系數(shù)對列車模型振動傳遞率的影響

二系的阻尼系數(shù)對兩自由度列車模型振動傳遞率的影響：可改變二系阻尼系數(shù)為0 N·s/m、1 000 N·s/m、3 000 N·s/m、5 000 N·s/m、7 000 N·s/m、10 000 N·s/m，一系阻尼系數(shù)保持為0 N·s/m，系統(tǒng)振動傳遞率曲線如圖3所示。

圖3 不同二系阻尼系數(shù)下二自由度系統(tǒng)振動傳遞率曲線

從圖3可以看出，振動傳遞率曲線主要分為5個區(qū)域，第1個區(qū)域在0.1Hz～2Hz，隨著頻率的升高，振動傳遞率的值變化較小，保持在0db，即基礎(chǔ)振動與二層質(zhì)量mv的振動幅值相同，振動沒有被放大，在這一區(qū)域內(nèi)阻尼系數(shù)的變化對振動傳遞率影響很?。坏?個區(qū)域在2Hz～5Hz，為一階共振區(qū)域，共振點約在3.5Hz，在一階共振區(qū)，隨著二系阻尼系數(shù)的增大，振動傳遞率顯著下降，即增大二系阻尼系數(shù)可以減小基礎(chǔ)振動的傳遞；第3個區(qū)域為5Hz～10Hz，隨著頻率的增大，振動傳遞率逐漸減小到0db以下，表明二層質(zhì)量mv的振動小于基礎(chǔ)振動，隨著阻尼系數(shù)的增大，振動傳遞率減小，在此區(qū)間使用越小阻尼系數(shù)的減振效果越好；第4個區(qū)域為10Hz～30Hz，此區(qū)域?qū)儆诙A共振區(qū)域，整體趨勢是隨著阻尼系數(shù)的增大，振動傳遞率減小；但是在共振點，阻尼系數(shù)增大，振動傳遞率降低；第5個區(qū)域為30Hz～100Hz，在此區(qū)域，隨著頻率的增大，振動傳遞率顯著降低，增大阻尼系數(shù)后振動傳遞率增大，越小阻尼系數(shù)的減振效果越好。

整體來說，要獲得較好的減振效果，在一、二階共振點的位置需要有大阻尼(3 000～10 000 N·s/m)；在其他的頻率范圍間，小阻尼(0～1 000 N·s/m)可以獲得較好的減振效果；改變二系阻尼系數(shù)可以有效控制系統(tǒng)的減振效果。

4 剪切閥式磁流變阻尼器阻尼力出力模型及最小阻尼系數(shù)

剪切式閥式磁流變阻尼器的出力模型主要基于磁流變液賓漢本構(gòu)方程與液壓公式，將磁流變阻尼器的出力簡化為與速度相關(guān)的粘滯力以及與磁流變液剪切強度相關(guān)的庫倫力。這種模型將剪切閥式磁流變阻尼器的出力與結(jié)構(gòu)參數(shù)相結(jié)合，將磁流變阻尼器的出力簡化為線性，可以較為方便地估算在任意結(jié)構(gòu)參數(shù)、運動速度、線圈電流下的磁流變阻尼器的出力范圍。

通過對4種不同規(guī)格的磁流變阻尼器[8-10]的理論出力以及實驗出力的比較，得到適應(yīng)2 300 N出力大小的剪切閥式磁流變阻尼器出力模型表達式如下：

(2)

其中：Fτ是阻尼力隨磁場而改變的庫侖力，即可控阻尼力；Fη是與速度有關(guān)的粘滯阻尼力；Ff為密封件與加工精度的問題會產(chǎn)生的摩擦力。公式中Ap表示受到剪切力作用的截面積，活塞直徑為D，活塞桿直徑為d，活塞與缸筒間隙為h，活塞的有效長度為l，阻尼間隙間的受剪面積為πDl，以上結(jié)構(gòu)參數(shù)如圖4所示。τ是磁流變液的剪切應(yīng)力，τy是受磁場影響的磁流變液的屈服應(yīng)力，η是磁流變液的塑性粘度，sgn為符號函數(shù)，v為活塞的運動速度。

圖4 剪切閥式磁流變阻尼器結(jié)構(gòu)參數(shù)示意圖

根據(jù)出力模型公式，磁流變阻尼器總的出力由與磁場相關(guān)的庫侖力和磁場無關(guān)的粘滯阻尼力組成，當(dāng)阻尼器的電流為0時，磁流變阻尼器為粘滯阻尼力，符合F=c·v。此時電流相關(guān)項均為0，阻尼力Fz為：

(3)

其中：Ff是由于密封圈或加工造成的初始摩擦力，占阻尼器的總力較小，為了方便求解阻尼系數(shù)，將其簡化。此時的阻尼力Fz為：

(4)

阻尼系數(shù)cz為:

(5)

即在電流為0 A的時候，阻尼系數(shù)由式(5)計算，此時的阻尼系數(shù)為定值，只與磁流變液的黏度系數(shù)和結(jié)構(gòu)參數(shù)相關(guān)，與速度無關(guān)；阻尼系數(shù)cz為磁流變阻尼器的最小阻尼系數(shù)。

5 磁流變阻尼器結(jié)構(gòu)參數(shù)優(yōu)化設(shè)計

設(shè)計一個固定出力的磁流變阻尼器時，需根據(jù)出力模型預(yù)先估算阻尼器的參數(shù)，比如活塞桿直徑、有效長度、外筒的內(nèi)徑；各參數(shù)的取值并不唯一，不同的結(jié)構(gòu)參數(shù)有可能對應(yīng)相同的最大阻尼力，但是各參數(shù)下的最小阻尼力會不同，即阻尼器阻尼力可調(diào)范圍不同。所以首先研究各參數(shù)對出力的影響，再對其進行優(yōu)化。

使用Matlab編寫程序研究基本參數(shù)(最大出力為2 300 N的磁流變阻尼器參數(shù)：外徑50 mm，工作長度30 mm，活塞桿直徑12 mm，工作間隙1 mm，最大電流3 A時磁流變液的屈服強度為30 kPa)下工作間隙h，工作長度l，活塞桿直徑d，活塞直徑D改變時對最大阻尼力的影響。

圖5 各參數(shù)對磁流變阻尼器出力的影響

然后使用Matlab優(yōu)化工具箱進行優(yōu)化，優(yōu)化的目標(biāo)是使磁流變阻尼器的最小阻尼力(即最小阻尼系數(shù))達到最小值，優(yōu)化的性能約束是使磁流變阻尼器的最大阻尼力達到固定值，結(jié)構(gòu)約束是活塞桿直徑、缸筒內(nèi)徑、有效長度在規(guī)定的范圍內(nèi)。

優(yōu)化目標(biāo)是最小阻尼系數(shù)，屬于求解約束條件下非線性極小值問題，使用fminimax函數(shù)；優(yōu)化的性能約束是使磁流變阻尼器最大阻尼力達到2 300 N。磁流變阻尼器參數(shù)優(yōu)化化過程如下：

1) 在Matlab中建立目標(biāo)函數(shù)，即式(5)。

目標(biāo)函數(shù)因變量為阻尼系數(shù)，自變量為結(jié)構(gòu)參數(shù)，包括有效長度、活塞直徑、活塞桿直徑。

2) 建立性能約束函數(shù)，即式(2)。

約束函數(shù)因變量為阻尼力，自變量為結(jié)構(gòu)參數(shù)，函數(shù)的實質(zhì)為阻尼器出力模型。

3) 建立結(jié)構(gòu)約束

有效長度、缸筒內(nèi)徑、活塞桿直徑的初始參數(shù)與變化范圍為：

x0=[0.03 0.05 0.01] ；

xl=[0.025 0.04 0.005] ；

xu=[0.035 0.06 0.015] ；

其中：x0是有效長度、活塞直徑、活塞桿直徑的初始參數(shù)，分別為30mm、50mm、10mm；xl為結(jié)構(gòu)參數(shù)的變化最小尺寸，xu為結(jié)構(gòu)參數(shù)的變化最大尺寸。

4) 調(diào)用優(yōu)化函數(shù)fminimax得到最終的優(yōu)化結(jié)果為：Xopt=0.031 0，0.047 8，0.013 0，即優(yōu)化后的有效長度，活塞直徑，活塞桿直徑的尺寸分別為31 mm，47.8 mm，13 mm，磁流變阻尼器的理論計算最大阻尼力為2 288 N；fopt=449.5，即優(yōu)化后的最小阻尼系數(shù)為449.5 N·s/m。

根據(jù)上面對磁流變阻尼器阻尼力的估算與優(yōu)化，并考慮加工的因素將尺寸進行圓整，確定最大目標(biāo)出力為2 300 N的磁流變阻尼器的主要結(jié)構(gòu)參數(shù)如表3，三維模型如圖6。

表3 優(yōu)化后的剪切閥式磁流變阻尼器結(jié)構(gòu)參數(shù)

圖6 剪切閥式磁流變阻尼器三維模型

6 結(jié)語

1) 列車磁流變阻尼器的設(shè)計方法基于列車整車阻尼比與磁流變阻尼器出力模型。該方法首先將列車模型簡化為單自由度系統(tǒng)，使用阻尼比公式計算列車模型垂向減振磁流變阻尼器的最大阻尼力。然后以最大阻尼力為優(yōu)化性能約束條件，以使磁流變阻尼器的可調(diào)范圍最大為目標(biāo)值，使用Matlab優(yōu)化工具箱得到磁流變阻尼器的結(jié)構(gòu)參數(shù)。

2) 以某型列車模型為例，表明利用最大阻尼力與出力模型對磁流變阻尼器進行結(jié)構(gòu)性能參數(shù)的設(shè)計方法是可行的，此方法可以較為方便快速地估算任意參數(shù)的列車垂向減振磁流變阻尼器的結(jié)構(gòu)參數(shù)。

3) 優(yōu)化設(shè)計的剪切閥式阻尼器達到了預(yù)期的設(shè)計目標(biāo)。預(yù)期設(shè)計目標(biāo)為2 300 N，最終設(shè)計的磁流變阻尼器的最大阻尼力為2 288 N，并且此時最小阻尼系數(shù)最小，阻尼器可調(diào)范圍最大。完成磁流變阻尼器的加工以及性能測試，驗證所設(shè)計的磁流變阻尼器是否與理論最大出力一致是下一步的研究工作。

[1] 唐啟志. 基于磁流變阻尼器的高速列車橫向半主動減振控制方法的研究[D]. 長沙:中南大學(xué),2012.

[2] Lee S H,Park J H,et al.Performance valuation of passive and semi-active magnetorheological dampers[J]. Journal of Mechanical Science & Technology,2007, 21(6):913-918.

[3] 劉宏友,曾京,李莉. 國內(nèi)外半主動懸掛系統(tǒng)的研發(fā)及應(yīng)用綜述[J]. 國外鐵道車輛,2012,49(6):12-18.

[4] 宗路航. 磁流變阻尼器的動力學(xué)模型及其在車輛懸架中的應(yīng)用研究[D]. 合肥:中國科學(xué)技術(shù)大學(xué),2013.

[5] 宗路航,陳現(xiàn)敏,郭朝陽,等. 磁流變阻尼器半主動控制方法的實驗研究[J]. 實驗力學(xué),2010,25(2):143-150.

[6] 顧永輝. 軌道車輛垂向振動性能研究及減振器阻尼參數(shù)選擇[D]. 上海:上海交通大學(xué),2008.

[7] 李異. 鐵路液壓減振器的應(yīng)用研究[D]. 成都:西南交通大學(xué),2007.

[8] 劉建軍. 新型磁流變阻尼器及其智能控制方法研究[D]. 天津: 天津大學(xué),2009.

[9] 關(guān)新春,歐進萍. 磁流變耗能器的阻尼力模型及其參數(shù)確定[J]. 振動與沖擊, 2001(1):5-8.

[10] 歐進萍,關(guān)新春. 磁流變耗能器及其性能 [J]. 地震工程與工程振動,1998,18(3):75-82.