武柏安 回學(xué)文 龍海洋 李耀剛

摘要：磁流變阻尼器由于自身獨(dú)特的優(yōu)點(diǎn)，在工程上受到了廣泛的關(guān)注。為了描述磁流變阻尼器的非線性和遲滯性，首先在現(xiàn)有的實(shí)驗(yàn)設(shè)備上進(jìn)行了阻尼器的標(biāo)定試驗(yàn)搜集數(shù)據(jù)，其次，在Simulink環(huán)境下建立了磁流變阻尼器的Dahl模型，對模型中的參數(shù)進(jìn)行了辨識(shí)，最終進(jìn)行了仿真對比。

關(guān)鍵詞：磁流變阻尼器;參數(shù)辨識(shí);Dahl模型

0 ?引言

磁流變阻尼器是一種內(nèi)部包含磁流變液體的阻尼器，施加電流的變化會(huì)改變電磁體中磁性粒子的密度跟強(qiáng)度，從而改變磁流變阻尼器的阻尼系數(shù)，作為一種新興的減振器，這種半主動(dòng)裝置可以應(yīng)用在車身與底盤之間的懸架上，能夠很好地緩解車身內(nèi)部發(fā)動(dòng)機(jī)與其他部件的產(chǎn)生的振動(dòng)和由于地面不平產(chǎn)生的顛簸。這種半主動(dòng)裝置能夠提高車輛行駛時(shí)的平順性和安全性。

磁流變阻尼器的位移與阻尼力之間存在著高強(qiáng)度的非線性關(guān)系，與固定阻尼不同的是，磁流變阻尼器的平均阻尼是可變的，它的遲滯曲線和飽和水平也與固定阻尼不同，給磁流變阻尼器的建模帶來了很大的難度。這些經(jīng)典的模型比如Bingham[1]模型、Bouc-Wen[2]模型和多項(xiàng)式[3]模型等。2015年Soltane等[3]提出了Bingham的正則化模型，在Bingham模型的基礎(chǔ)上引入正則化參數(shù)的方法，將不連續(xù)的Bingham方程替換為連續(xù)方程，可以非常逼近原模型的值。2019年祝世興等[4]在Bingham的模型的基礎(chǔ)上，將阻尼器的速度分段，建立了線性模型和非線性，結(jié)果線性模型在低速區(qū)不具有很高的精度，非線性模型能夠很好地模擬出阻尼器的特性。2005年H.H. Tsang等[5]針提出了一種簡化的逆模型，通過任意的控制方法可以達(dá)到控制磁流變阻尼器的作用。2019年張麗霞等[6]建立了磁流變阻尼器的魔術(shù)公式，其模型簡單且物理意義明確并獲得了誤差控制在8%以內(nèi)的模型。本文對Lord的一款磁流變阻尼器進(jìn)行參數(shù)化建模，在不同工況下進(jìn)行了實(shí)驗(yàn)研究，采用得到的數(shù)據(jù)在Simulink中建立了Dahl模型并對模型中的參數(shù)進(jìn)行了辨識(shí)，并進(jìn)行了仿真對比。

1 ?阻尼器標(biāo)定試驗(yàn)

如圖1所示，本次試驗(yàn)選用Lord公司生產(chǎn)的磁流變阻尼器，它的工作參數(shù)為：最大行程248mm，外徑42.1mm，活塞桿直徑10mm，最大輸出阻尼力2447N，最大輸入電流1A響應(yīng)時(shí)間小于15ms，所示實(shí)驗(yàn)采用的測試裝置為杭州億恒科技生產(chǎn)的阻尼器試驗(yàn)系統(tǒng)，該系統(tǒng)有PC機(jī)，電液伺服控制器、阻尼試驗(yàn)臺(tái)、液壓油源系統(tǒng)和溫度控制系統(tǒng)組成。

試驗(yàn)位移加載方式為正弦波加載，實(shí)驗(yàn)條件如表1所示。

圖2、圖3所示為阻尼器在位移15mm、頻率2Hz不同電流下的特性曲線，從中可以看出強(qiáng)烈的滯回特性和非線性。不同工況下的曲線相似，得到多組實(shí)驗(yàn)數(shù)據(jù)，為后面的建模進(jìn)行了鋪墊。

2 ?Dahl模型的建立

Dahl阻尼器數(shù)學(xué)模型，由以下兩個(gè)方程式組成：

其中：F為磁流變阻尼器的輸出力，N;x為阻尼器活塞與缸體的相對位移，mm;A為滯回曲線的形狀參數(shù)，s/mm;Fd為可調(diào)庫倫摩擦力，N;f0為阻尼器內(nèi)部初始力，N;K0為剛度系數(shù)，N/mm;C0為粘滯阻尼系數(shù)，Ns/mm;

將Dahl模型在Simulink中進(jìn)行搭建，模型有5個(gè)未知參數(shù)，本文采用Simulink中的Simulink Design Optimization進(jìn)行識(shí)別。輸入數(shù)據(jù)為阻尼器的實(shí)驗(yàn)測得的位移和速度，輸出數(shù)據(jù)為實(shí)驗(yàn)測得的阻尼力。由于實(shí)驗(yàn)數(shù)據(jù)較多，計(jì)算量較大，所以將各個(gè)工況下得到的辨識(shí)結(jié)果進(jìn)行平均計(jì)算。

參數(shù)最終的擬合結(jié)果用公式（2）表示：

通過導(dǎo)出的擬合模型來進(jìn)行驗(yàn)證，將模擬結(jié)果與圖中所示的MR阻尼器實(shí)驗(yàn)條件下的進(jìn)行比較。

圖5中相同顏色的圖線代表同一電流下的阻尼力位移圖線，實(shí)線代表仿真值，虛線代表實(shí)際值，從圖中可以看出，仿真值和實(shí)際值在位移上存在偏差，那是因?yàn)樵谧鲎枘崃μ匦栽囼?yàn)時(shí)，主要是因?yàn)樵囼?yàn)時(shí)設(shè)備誤差造成。經(jīng)數(shù)值比較以后，存在2mm左右的誤差，屬于正常誤差而且在最大位移時(shí)，磁流變阻尼器的速度方向發(fā)生變化，這時(shí)候產(chǎn)生的沖擊振動(dòng)使得仿真值和實(shí)驗(yàn)值之間的誤差較大，使模型存在一定的誤差。

3 ?結(jié)論

本文采用億恒公司的試驗(yàn)設(shè)備進(jìn)行阻尼器標(biāo)定試驗(yàn)，在Simulink環(huán)境下建立的磁流變阻尼器的Dahl模型，對模型的參數(shù)進(jìn)行了辨識(shí)，并進(jìn)行了仿真對比。為磁流變阻尼器的應(yīng)用做出了理論準(zhǔn)備。

參考文獻(xiàn)：

[1]SPENCER Jr B F， DYKE S J， SAIN M K， et al. PHENOMENOLOGICAL MODEL FOR MAGNETORHEOLOGICAL DAMPERS[J]. Journal of Engineering Mechanics， 1997， 123（3）：230-238.

[2]GUO S， YANG S， PAN C. Dynamic modeling of magnetorheological damper behaviors[J]. Journal of Intelligent material systems and structures， 2006， 17（1）：3-14.

[3]SOLTANE S， MONTASSAR S， MEKKI O B， et al. A hysteretic Bingham model for MR dampers to control cable vibrations[J]. Journal of Mechanics of Materials and Structures， 2015， 10（2）：195-206.

[4]祝世興，郝新琛.基于Bingham模型的磁流變阻尼器模型改進(jìn)研究[J].機(jī)床與液壓，2019，47（17）：98-102，128.

[5]TSANG H H， SU R K L， CHANDLER A M. Simplified inverse dynamics models for MR fluid dampers[J]. Engineering Structures， 2006， 28（3）：327-341.

[6]張麗霞，龐齊齊，潘福全，等.磁流變減振器魔術(shù)公式模型在懸架控制中的應(yīng)用[J].中國機(jī)械工程，2019：1-8.

[7]PENG G R， LI W H， DU H， et al. Modelling and identifying the parameters of a magneto-rheological damper with a force-lag phenomenon[J]. Applied Mathematical Modelling， 2014， 38（15-16）：3763-3773.