磁流變阻尼器擬靜力力學特性及力學模型

張香成 毋貴斌 趙軍

摘要：為研究磁流變阻尼器（MRD）在擬靜態作用下的力學性能，制作了一個四線圈剪切閥式MRD，在不同電流和不同位移幅值下采用三角位移加載制度對MRD進行了擬靜態試驗，研究了MRD的力學性能隨電流和位移幅值的變化規律，分析了阻尼力波動現象產生的原因;根據試驗結果提出了一種擬靜態作用下修正的力學模型，對模型參數進行了識別，對力學模型進行了驗證。結果表明：所用阻尼器在擬靜態試驗過程中阻尼力波動現象顯著，無電流時阻尼力隨著位移幅值增加而略微減小，有電流時阻尼力隨著位移幅值增加而稍微增大;阻尼力隨電流的增大明顯增大，當電流超過1.5 A時，MRD達到磁飽和狀態，阻尼力增幅明顯降低;提出的力學模型能夠較好地描述擬靜力下阻尼力隨電流和位移的改變而產生的變化特性。

關鍵詞：磁流變阻尼器;擬靜態試驗;力學模型;參數識別

中圖分類號：TU352.11 文獻標志碼：A 文章編號：20966717（2020）03004608

Abstract： In order to study the mechanical properties of Magnetorheological Damper （MRD） under quasistatic action， a fourcoil shear valve MRD was fabricated. A quasistatic test was carried on the MRD by using the triangular displacement loading system under different currents and displacements. In the test， the variation of the mechanical properties of MRD with the changes of current and displacement was studied， and the causes of the damping force fluctuation phenomenon were analyzed. Based on the test results， a modified mechanical model under quasistatic loading was proposed， the parameters of the model were identified， and the mechanical model was verified. The results show that the damping force of MRD fluctuates remarkably during the quasistatic test process. The damping force decreases slightly with the increase of displacement amplitude when there is no current， and increases slightly with the increase of displacement amplitude when there is current. The damping force increases obviously with the increase of current， and when the current exceeds 1.5A， MRD reaches magnetic saturation and the increase of damping force decreases obviously. The mechanical model proposed in this paper can well describe the variation of damping force with the change of current and displacement under quasistatic force.

Keywords：magnetorheological damper; quasistatic test; mechanical model; parameter identification

磁流變阻尼器（Magnetorheological damper， MRD）是一種可調節阻尼的減振控制裝置，具有結構簡單、出力大、阻尼力連續可調且控制效果好的優點，在土木工程減振防災方面展現出了良好的應用前景，已成為目前結構振動控制的重要裝置[17]。近年來，學者們對MRD做了大量研究，梅真等[8]對MRD進行了動力學試驗，建立了雙曲正切滯回模型和BP神經網絡正向、逆向力學模型，并對這兩種模型進行驗證;杜成斌等[9]對MRD在低頻下力學性能進行研究，并對非線性滯回雙粘性模型進行了修正;張香成等[10]基于米氏方程提出了米氏模型，并通過試驗加以驗證。

對設置MRD的結構（或結構構件）進行擬靜力試驗是研究其力學性能、耗能性能和二者相互作用機理的經濟、可靠手段，然而，關于MRD擬靜力性能的研究較少，因此，在研究設置MRD結構（或結構構件）的力學性能之前，應當先對MRD的擬靜力性能進行研究。Zekeriya等[11]采用準靜態分析的方法研究了MRD中的非牛頓流動，并對MRD進行了計算流體動力學分析;Mohammad等[12]對MRD中的磁流變流體流動進行了準靜態分析，說明了磁流變流體慣性產生的影響，但這些研究并未描述MRD的擬靜力力學性能。

筆者采用位移加載方式，利用三角波加載函數在低速下對自制的MRD進行擬靜態加載，得到不同工況下MRD的阻尼力位移曲線。通過對結果進行對比分析，研究MRD在不同位移幅值、不同電流下的力學性能，分析試驗中MRD阻尼力波動等現象出現的理論原因?；谠囼灲Y果提出一種擬靜力作用下修正的MRD阻尼力計算模型，并對力學模型中的參數進行了識別。

1阻尼器擬靜態試驗

1.1阻尼器結構

采用四線圈剪切閥式阻尼器，阻尼器多級活塞材料為電工純鐵，其余部分為45號鋼，其結構如圖1所示。基座連接板利用螺栓將阻尼器與下方基座連接，線圈導線從MRD多級活塞和活塞桿中心的圓柱形孔道引出。阻尼器各部分參數見表1。

阻尼器內部填充的磁流變液采用自制方法獲得，其組分為羰基鐵粉、硫酸鈉、液體石蠟、皂土、石墨以及甲基硅油等，其中，羰基鐵粉占60%、甲基硅油占28%（質量分數）。磁流變液的相關性能在中國科學技術大學智能材料與振動控制實驗室測定，在磁場強度為0.44 T時，當剪切速率分別為0.01、0.34/s-1時，磁流變液的剪切應力分別為9 456.2、9 492.2 Pa，二者相差0.38%;當剪切速率為100/s-1時，磁流變液在不同磁場下的剪切應力見表2。

1.2試驗加載

試驗通過作動器進行位移控制，作動器的最大出力為100 kN，采用擬靜態下的三角波對阻尼器進行位移加載;基座連接板和基座中心均開有圓孔，線圈導線穿過圓孔后從基座側面引出，并與直流電源連接;試驗過程中，MRD由穩壓直流電源供電，且MRD 4個線圈中的電流等級保持一致，電流等級從小到大依次為0、0.5、1.0、1.5、2.0、2.3、2.6 A，在每個電流等級下，均采用圖2所示的加載制度對MRD進行加載，共進行了7次加載試驗，未進行重復性試驗，試驗裝置見圖3。

2試驗結果分析

試驗通過改變加載電流、阻尼器位移和加載時間，測得阻尼器在不同工況下的試驗數據，通過MATLAB軟件繪制，得到相關的阻尼力位移曲線。由于圖2中MRD試驗的最小、最大加載速度分別為0.03、0.33 mm/s，與之對應的阻尼間隙處磁流變液的剪切速率在該剪切速率范圍內，磁流變液的剪切應力相差極小，因此，在MRD試驗結果分析中忽略速率的影響。

2.1力學性能隨位移幅值的變化規律

MRD在電流為0、1.5、2.6 A時，不同位移幅值下的試驗結果分別如圖4、圖5和圖6所示。

從圖4可以看出，當MRD通電電流為0 A時，在MRD位移幅值較小時，阻尼器所提供的阻尼力值約為0.8 kN。隨著MRD位移幅值的增加，MRD所提供的阻尼力逐漸降低。當MRD位移幅值為30 mm時，MRD所提供的阻尼力值約為0.5 kN。此時，MRD在較小位移幅值下阻尼力較大，而在較大位移幅值時阻尼力較小。

從圖5、圖6可以看出，隨著MRD通電電流增加，阻尼器所提供的阻尼力也隨之增加。在通電的情況下，當位移幅值較小時，MRD的阻尼力也較小;隨著MRD位移幅值的增加，阻尼力也會增加。此外，由圖5、圖6可以看出，阻尼力迅速上升至最大值后，會在一定區間內不斷波動。

由于MRD試驗為靜力加載，其速度極慢，速度所產生的粘滯力很微小，可忽略不計;在低電流時，磁流變液本身剪切強度低，剪切力與擠壓力也較小。此時，MRD內部摩擦力占主導地位。

摩擦力主要由活塞桿與密封圈（MRD采用O型圈密封）間相對滑動提供，摩擦力大小的變化可分為兩個階段：啟動階段活塞桿開始運動，摩擦力逐漸增大，直到達到O型圈啟動摩擦力，由于在達到相對滑動前，O型圈會產生預位移，密封圈自身產生變形，所提供摩擦力較最大靜摩擦力更大;之后為滑動階段，此時活塞桿開始穩定運動，O型圈處于滑動狀態，預位移得到釋放，O型圈與摩擦桿之間的摩擦力趨于穩定，此時為滑動摩擦，摩擦力比最大靜摩擦力要小。在MRD位移幅值較小時，阻尼器活塞的運動速度極慢，密封圈與活塞桿之間一直處于靜摩擦滑動摩擦靜摩擦不斷交替的狀態[14]，此時的摩擦力以靜摩擦為主，摩擦力較大;在MRD位移幅值增加后，阻尼器活塞運動速度變快，密封圈與活塞桿之間穩定處于動摩擦狀態，摩擦力較小。因此，在低電流時，摩擦力在阻尼力中處于主導地位，出現了小位移幅值下阻尼力較大而在大位移幅值時阻尼力較小的現象。

由式（1）可知，阻尼力部分包括擠壓力Fn，而Fn與活塞移動位移成正比[15]。當阻尼器線圈通電流后，磁流變液抗壓屈服強度迅速增加，此時，隨著活塞移動位移的增加，擠壓力也隨之增加，進而導致阻尼力隨著位移的增加而增加。

2.2力學性能隨電流的變化規律

不同電流下阻尼器的工作狀況如圖7和圖8所示。當電流較小時，MRD所提供的阻尼力也較小;隨著電流的增加，阻尼力也在不斷增加;電流達到1.5 A后，阻尼力增加不再明顯，說明MRD已經達到磁飽和狀態。

2.3阻尼力波動現象成因分析

試驗采用剪切閥式阻尼器，由式（1）可知，在線圈通上一定電流后，磁流變液剪切力將成為阻尼力的主要部分。

由阻尼力的試驗曲線可以看出，在MRD活塞桿的運動過程中，阻尼力沒有穩定地保持在某一數值上，而是在某一區間內上下波動。

這種現象出現的原因是：如圖9（a）所示，磁流變液中的磁性粒子在磁場的作用下整齊排列成鏈[16]，在活塞與缸體產生相對滑動時，阻尼通道中的磁流變液層間也會產生相對滑動，且隨著剪切速率增加，顆粒鏈的傾斜角也隨之增加，如圖9（b）所示。在外界磁場和剪切作用的共同作用下，顆粒鏈將不斷斷裂重組，直至達到動態平衡狀態[17]。由于磁流變液并非完全勻質，且擬靜力加載下，剪切速率低，加載時間長，磁性粒子間的動態平衡無法保持完全穩定，隨著活塞運動，打亂、破壞的磁性粒子鏈數量也在變化。因此，在磁性粒子鏈不斷被打亂、重組的過程中，阻尼力也在小范圍內不斷波動。

3力學模型及參數識別

3.1阻尼器計算模型

由于在擬靜態加載方式下，阻尼器活塞相對運動速度低，Cu·可忽略不計;MRD的擠壓阻尼力Fn與活塞移動位移u成正比。

3.2計算模型參數識別

首先確定Ff。在電流為0 A時，無外加磁場，磁流變液呈現低粘度的牛頓流體特性，此時，磁流變液剪切力Fτ和擠壓力Fn均為0，MRD的阻尼力由摩擦力Ff提供。由此可得出，在各循環時，MRD摩擦力Ff的大小如表3所示。

在擬靜態加載下，加載速度極低，MRD剪切阻尼力Fτ的大小只與磁流變液剪切屈服強度有關，因此，在相同電流下，MRD剪切阻尼力Fτ不變。

將求得的摩擦力及不同工況下試驗結果帶入式（3）中，通過非線性擬合，可得到不同電流下的系數k，以及不同加載電流下MRD的剪切阻尼力Fτ，如表4、表5所示。

由于Fτ是MRD在不同磁感應強度下的屈服應力，與MRD的位移幅值和激勵頻率無關，只與電流強度有關。Fn是MRD在不同磁感應強度下的擠壓應力，只與磁場強度以及活塞與筒壁距離有關，式（2）中，u為活塞移動位移，而磁場強度與電流強度正相關，故k為與電流強度有關的系數。

采用非線性最小二乘法，根據式（5）、式（6）對Fτ和k進行非線性擬合，可得式中各參數值分別為：Fτ0=23.03、a=1.464、k0=561.4、b=0.304 1。圖10、圖11分別為參數Fτ和k的試驗值和擬合值之間的對比。以圖10為例，在電流分別為0、0.5、1.0、1.5、2.0、2.3、2.6 A時，試驗測得Fτ分別為0、11.58、18.33、20.24、21.33、22.73、22.39 kN，相應電流下的擬合值為0、11.95、17.70、20.47、21.79、22.23、22.51 kN，試驗值和擬合值非常接近。

3.3模型對比驗證

為了驗證模型精度，將得到的模型結果與試驗所得滯回曲線進行對比，如圖12、圖13所示。同時，通過試驗數據中實測位移計算各數據點對應的模型值，與實測的阻尼力進行比較計算，取平均值得到了相應電流下試驗測得阻尼力與模型計算阻尼力的平均誤差，見表6、表7。

從圖12和圖13可以看出，所得力學模型與試驗結果較吻合，不僅能很好地描述出MRD在擬靜態加載制度下阻尼力隨電流的增加而增加的特征，也能夠模擬阻尼力在不同電流下隨位移幅值變化的特性，即：在無電流時，阻尼力隨著位移幅值的增加而減小;在通電流時，阻尼力隨著位移幅值的增加而增加。但是，該模型不能描述MRD靜態試驗時阻尼力隨著活塞的運動而產生的上下波動情況。

從表6和表7可以看出，24組數據中共有21組誤差小于等于5%;誤差較大的3組數據：電流0 A時，位移幅值2 mm，電流2.0 A時，位移幅值6、8 mm，均為較小位移幅值時的工況。由于MRD在制作時無法保證磁流變液完全填滿缸筒內部，會在磁流變液中存在氣泡，這種情況在位移幅值較小時會導致模型誤差相對較大。而在位移幅值較大的工況中，這一模型誤差所造成的影響將會降低，因此，試驗數據與力學模型間的誤差也更小。

4結論

1）擬靜態試驗下，無電流時，小位移值下阻尼力較大，而在大位移時阻尼力較小;有電流時，小位移值下阻尼力較小，而在大位移時阻尼力較大。

2）MRD所提供的阻尼力隨著電流的增加而增加;當電流達到1.5 A以后，阻尼力增加不再明顯，說明MRD已經達到磁飽和狀態。

3）在擬靜態試驗過程中，MRD的阻尼力會產生明顯的波動現象。

4）利用試驗數據對提出的擬靜態加載下修正的力學模型進行參數識別，結果表明，提出的力學模型能夠較好地描述擬靜力下MRD阻尼力隨電流和位移的改變而產生的變化特性。

由于課題組試驗樣本不夠豐富，力學模型還不能與MRD的基本參數（如表1中的線圈匝數、活塞桿直徑等，表2中的剪切應力）進行關聯，并且該模型不具備普適性，因此，不能直接預測其他MRD的擬靜力力學性能，在以后的研究中需要以此為目標進行改進。參考文獻：

[1] ?FERDAUS M M， RASHID M M， HASAN M H， et al. Optimal design of magnetorheological damper comparing different configurations by finite element analysis [J]. Journal of Mechanical Science and Technology， 2014， 28（9）： 36673677.

[2] 胡國良， 張佳偉. 磁流變阻尼器結構優化設計研究現狀[J]. 機床與液壓， 2019， 47（1）： 145150.

HU G L， ZHANG J W. State of the art review on structural optimization design of magnetorheological damper [J]. Machine Tool & Hydraulics， 2019， 47（1）： 145150. （in Chinese）

[3] ?CHOOI W W， OYADIJI S O. Design， modelling and testing of magnetorheological （MR） dampers using analytical flow solutions [J]. Computers & Structures， 2008， 86（3/4/5）： 473482.

[4] ?DAVIS L C. Model of magnetorheological elastomers [J]. Journal of Applied Physics， 1999， 85（6）： 33483351.

[5] ?SHOU M J. Dynamic behavior of magnetorheological energy absorber under impact loading [J]. Journal of Mechanical Engineering， 2019， 55（1）： 72.

[6] ?LI R， ZHOU M J， WU M J， et al. Semiactive predictive control of isolated bridge based on magnetorheological elastomer bearing [J]. Journal of Shanghai Jiaotong University （Science）， 2019， 24（1）： 6470.

[7] ?GURUBASAVARAJU T M， KUMAR H， MAHALINGAM A. An approach for characterizing twintube shearmode magnetorheological damper through coupled FE and CFD analysis [J]. Journal of the Brazilian Society of Mechanical Sciences and Engineering， 2018， 40（3）： 139.

[8] 梅真， 高毅超， 郭子雄. 磁流變阻尼器動力性能測試與建模[J]. 振動.測試與診斷， 2017， 37（3）： 553559， 632.

MEI Z， GAO Y C， GUO Z X.Dynamic testing and modeling of a magnetorheological damper [J]. Journal of Vibration， Measurement & Diagnosis， 2017， 37（3）： 553559， 632. （in Chinese）

[9] 于國軍， 杜成斌. 一種復合磁流變阻尼器的設計及性能試驗[J]. 力學季刊， 2014， 35（1）： 131138.

YU G J， DU C B. Design and experimental studies on composite magnetorheological damper [J]. Chinese Quarterly of Mechanics， 2014，35（1）： 131138. （in Chinese）

[10] 張香成， 徐趙東， 王紹安， 等. 磁流變阻尼器的米氏模型及試驗驗證[J]. 工程力學， 2013， 30（3）： 251255.

ZHANG X C， XU Z D， WANG S A， et al. Michaelismenten model of magnetorheological damper and test verification [J]. Engineering Mechanics， 2013， 30（3）： 251255. （in Chinese）

[11] ?PARLAK Z， ENGIN T. Timedependent CFD and quasistatic analysis of magnetorheological fluid dampers with experimental validation [J]. International Journal of Mechanical Sciences， 2012， 64（1）： 2231.

[12] ?ZOLFAGHARIAN M M， KAYHANI M H， NOROUZI M， et al. Parametric investigation of twin tube magnetorheological dampers using a new unsteady theoretical analysis [J]. Journal of Intelligent Material Systems and Structures， 2019， 30（6）： 878895.

[13] ?SPENCER B F J， DYKE S J， SAIN M K， et al. Phenomenological model for magnetorheological dampers [J]. Journal of Engineering Mechanics， 1997， 123（3）： 230238.

[14] 馬春紅， 白少先， 康盼. 氟橡膠O型圈低壓氣體密封黏滯摩擦特性實驗[J]. 摩擦學學報， 2014， 34（2）： 160164.

MA C H， BAI S X， KANG P. Experiment of viscous friction characteristics of fluorous rubber Orings at low gas seal pressure [J]. Tribology， 2014， 34（2）： 160164. （in Chinese）

[15] 姚軍， 張進秋， 賈進峰， 等. 擠壓式磁流變減振器力學模型研究[J]. 中國工程機械學報， 2015， 13（6）： 497503.

YAO J， ZHANG J Q， JIA J F， et al. Study on mechanics model of squeeze flow type of magnetoreological fluid damper [J]. Chinese Journal of Construction Machinery， 2015， 13（6）： 497503. （in Chinese）

[16] 朱俊濤， 徐趙東， 張香成. 磁流變彈性體主鏈吸附模型[J]. 東南大學學報（自然科學版）， 2011， 41（2）： 342346.

ZHU J T， XU Z D， ZHANG X C. Main chain adsorption model of magnetorheological elastomer [J]. Journal of Southeast University（Natural Science Edition）， 2011， 41（2）： 342346. （in Chinese）

[17] 阮曉輝. 磁流變液力學性能及其應用研究[D]. 合肥： 中國科學技術大學， 2017.

RUAN X H. Study on the mechanical properties and applications of magnetorheological fluid [D]. Hefei： University of Science and Technology of China， 2017.（in Chinese）

[18] 張香成， 何尚文， 李倩， 等. 鉛磁流變阻尼器的試驗及計算模型[J]. 工程力學， 2016， 33（10）： 123128.

ZHANG X C， HE S W， LI Q， et al. Test and mathematic model of lead magnetorheological damper [J]. Engineering Mechanics， 2016， 33（10）： 123128. （in Chinese）

（編輯王秀玲）