參數慢變磁流變非線性懸架系統主共振研究

董小閔，王小龍，古曉科

（重慶大學機械傳動國家重點實驗室，重慶 400044）

參數慢變磁流變非線性懸架系統主共振研究

董小閔，王小龍，古曉科

（重慶大學機械傳動國家重點實驗室，重慶 400044）

針對磁流變阻尼器隨工作溫度變化導致阻尼力發生改變，進而降低磁流變懸架系統性能的問題，基于傳熱學方程和參數慢變非線性振動理論研究溫度變化對磁流變非線性懸架系統的主共振動力學行為的影響。根據傳熱學理論建立磁流變阻尼器工作區域能量微分方程，用數值方法求出溫度變化規律，用改進的Bingham模型計算磁流變阻尼力；建立慢變參數磁流變非線性懸架系統動力學方程，利用平均法求解主共振激勵下的懸架動行程的穩態幅值響應。研究表明：隨著運行時間增加，磁流變阻尼器的溫度逐漸升高，導致磁流變阻尼力大幅降低，懸架動行程穩態振幅及其變化量越大，從而使得懸架系統的性能偏離設計目標。

慢變參數；磁流變阻尼器；溫度；非線性振動

磁流變阻尼器作為一種智能阻尼器件，具有結構簡單、響應快、能耗低、動態范圍寬等優點，在結構振動控制領域受到了極大的關注。但近年來隨著研究的深入，研究人員發現，磁流變阻尼器件會由于環境或自身原因導致溫度的升高，引起輸出阻尼力大幅降低［1］，導致系統動力學行為偏離設計目標，影響系統性能［1－2］。因此，有必要理解溫度對系統動力學行為的影響規律，從而為不確定性懸架系統的優化和控制策略的調整提供研究基礎。

溫度引起磁流變阻尼器輸出阻尼力的變化過程與控制周期相比，是一個緩慢變化的過程，磁流變振動系統屬于參數慢變的非線性動力學系統，對其動力學行為的分析在理論上有一定的難度。目前國內外在溫度對磁流變阻尼器的影響研究方面主要集中于磁流變阻尼器溫升理論建模和控制補償兩個方面，而溫度對于磁流變減振系統動力學行為影響方面的研究還很少。如Gordaninejad等［3］利用集總參數法建立了磁流變阻尼器溫度的理論模型，實驗驗證該理論模型可以較好地預測溫升。Liu等［4］研究了考慮溫度補償的1／4車輛懸架的天棚阻尼控制。Wilson等［5］針對一種遠征戰斗機（Expeditionary Fighting Vehicle）座椅，通過建立考慮溫度變化的非線性Bingham模型（Nonlinear Bingham Plastic Model）提高了半主動天棚阻尼控制的魯棒性。

目前慢變參數的非線性振動研究主要關注于系統質量和剛度發生緩慢變化時的解析方法，對于系統阻尼（如磁流變系統）緩慢發生變化的研究還比較少。如Kevorkian［6］的研究表明，如果不是選擇恰當的快、慢尺度，傳統的多時間尺度法失效。聞邦椿等［7］應用KBM法研究含有物料慢變參數的振動壓實系統；陳予?。?］應用KBM法研究具有任意個準循環坐標的多自由度系統的非定常解。李以農等［9］針對一類具有參數慢變的非線性振動系統，討論了慢變系統的分析方法，結合轉子系統的慢變剛度對其振動特性進行了研究。

為此，本文針對汽車磁流變懸架系統在運行過程中動力學行為隨溫度發生變化的問題，提出利用傳熱學和參數慢變非線性振動理論分析溫度變化對汽車懸架系統動力學行為的影響。

1 溫度影響下磁流變阻尼器輸出阻尼力計算

磁流變阻尼器作為一種能量耗散裝置，能將外界輸入的振動能轉換為自身的熱力學能并耗散，該熱力學能則表現為磁流變阻尼器在服役條件下的內部溫升。Gordaninejad等［3］假設系統各組成部分在任意時刻具有相同的溫度，利用集總參數法研究了磁流變阻尼器的溫度隨時間變化的規律。集總參數法是非穩態導熱過程中物體的溫度場的一種求解方法。在實際應用中，要求系統各點過余溫度的最大偏差不超過5%，對應的特征畢奧數≤0．1。

假設液壓筒內壁的溫度與磁流變液的溫度相同，利用集總參數法，建立磁流變阻尼器的能量微分方程

導出的熱量以對流傳熱的方式傳熱到周圍空氣中，忽略熱輻射傳導方式，對流傳熱過程可表示為

式中：h為液壓筒壁和周圍空氣間的自然對流傳熱系數，As為阻尼器的外表面積，Θ0為阻尼器外圍溫度。

磁流變阻尼器中磁流變液的熱力學能，由磁流變阻尼器做功和勵磁線圈的電阻熱組成。

式中：F（t）為與溫度有關的磁流變阻尼力，X（t）為隨時間變化的位移激勵，R、I分別為勵磁線圈的電阻和加載的電流。

設激勵為正弦激勵，表示為

目前國內外學者在計算磁流變阻尼力時，考慮了溫度變化因素的有非線性動態逆向模型（Nonlinear Dynamic Inverse Model，NLDIM）［2］、非線性加權雙曲正切模型（Nonlinear Weighted Hyper-Tangent Model）［5］、溫度唯象模型［10］、指數形式模型［4］和改進Herschel-Bulkley［11－12］等模型，但除了Herschel-Bulkley之外，其中絕大多數模型都不便于動力學方程的計算和分析。因此，本文利用Herschel-Bulkley模型計算考慮溫度變化時的磁流變阻尼力。

Herschel-Bulkley本構關系模型為

為簡化計算，取n＝1，Herschel-Bulkley模型即為改進的Bingham模型，從而得到磁流變阻尼力

式中，C1和C2分別為與粘性阻尼和可控阻尼力相關的系數，由磁流變阻尼器的關鍵結構尺寸參數而定。

式（7）和式（8）可由數值方法聯立求解，將式（9）代入式（7），即可求得阻尼力關于時間的表達式。

在得出磁流變阻尼力隨溫度的變化規律后，研究溫度變化引起的阻尼力變化對主共振下汽車懸架系統動力學行為的影響。

2 慢變參數非線性汽車懸架系統主共振分析

結合第2節計算得到的磁流變阻尼力的表達式，以下研究如圖1所示的主共振汽車磁流變懸架系統，建立動力學方程

式中：m為系統的質量，x和x分別為系統的垂直速度和加速度，x0為位移激勵，k為剛度，F（τ）為溫度慢變的磁流變阻尼力，由式（6）～式（9）確定，x（0）＝x （0）＝0

圖1 單自由度汽車磁流變懸架系統Fig．1 SDOF vehicle MR suspension system

目前求解慢變參數動力學問題有效的方法主要有漸進法、多尺度法和平均法［7－9，13］。但漸進法只能用于求解激勵力較小的動力學問題。多尺度法雖然可以用來求解較大激勵力的動力學問題，但當動力學方程中有如可控Bingham力形式的靜載荷時，不方便求解永年項。因此，本文選用平均法來求解。

在求得穩態響應幅值后，通過數值方法分析穩態響應幅值隨時間的變化規律。

通過穩態振幅a對溫度T求導數，可得到穩態振幅隨溫度的慢變規律，其表達式如下：

3 數值仿真

［3，15－16］，并結合實際的汽車磁流變懸架系統，選取仿真參數：As＝0．039 m2，R＝6 Ω，h＝2．5 W／m－2℃，λ＝0．000 85，m＝240 kg，k＝16 000 N／m，a＝0．003 8，c（0）＝1 100 Ns·m－1，b＝0．002 5。

圖3 磁流變阻尼力變化規律Fig．3 MR damping force-time history

圖4為溫度變化時懸架動行程的穩態振幅，加載電流分別為0．5 A、1．0 A、1．5 A和2．0 A。結合圖2給出的磁流變阻尼器內部溫度的變化規律，從圖4中可以看出，隨著運行時間增長，溫度的升高，汽車磁流變懸架系統動行程穩態振幅變大，即溫度的升高導致汽車磁流變懸架系統的隔振效果變差。此外，隨著加載電流的增大，穩態振幅的最大變化量依次為1．1%、1．9%、3．2% 和4．6%，呈逐漸增大的趨勢。

圖5為不同加載電流下懸架動行程穩態振幅隨慢變參數溫度的靈敏度分析。結合圖2和圖4中給出的磁流變阻尼器內部溫度變化規律和汽車磁流變懸架動行程穩態振幅隨電流的變化規律，可以看出，隨著電流的逐漸增大，磁流變阻尼器內部的溫度升高的越快，汽車懸架系統動行程穩態振幅也逐漸增大，同時，穩態振幅對溫度變化的靈敏程度也逐漸增大。

圖4 汽車磁流變懸架動行程穩態振幅隨電流的變化Fig．4 Steady-state amplitude of suspension stroke as a function of input electric current for MR damper

4 結 論

本文針對汽車磁流變懸架系統在運行過程中阻尼力隨溫度發生變化的問題，利用傳熱學和參數慢變非線性振動理論研究了非穩態溫度慢變汽車磁流變懸架系統的主共振動力學行為。隨著磁流變阻尼器運行時間增長，溫度的逐漸升高，引起磁流變阻尼力大幅降低，導致懸架系統隔振性能偏離設計目標。此外，隨著加載電流逐漸增大，磁流變阻尼器內部溫度的升高越快，因此也導致懸架系統動行程穩態振幅的變化量也越大，穩態振幅對溫度的靈敏度也越高，隔振性能逐漸變差。研究表明，為了充分發揮磁流變阻尼器的性能，在非穩態實時控制時需要進行合理的補償，控制補償的研究和實驗驗證將在后續工作中進一步進行。

圖5 不同加載電流懸架動行程穩態振幅靈敏度分析Fig．5 Steady-state amplitude sensitivity analysis of suspension stroke as a function of input electric current for MR damper

參考文獻

［1］Batterbee D，Sims N D．Temperature sensitive controller performance of MR Dampers［J］．Journal of Intelligent Material Systems and Structures，2008：1－13．

［2］Jumani S S．An invertible open-loop nonlinear dynamic temperature dependent MR damper model［D］．Master of Science In Mechanical Engineering，2010．

［3］Gordaninejad F，Breese D G．Heating of Magneto-rheological fluid dampers［J］．Journal of Intelligent Material Systems and Structures，1999，10（8）：634－645．

［4］Liu Y，Gordaninejad F，Evrensel C A，et al．Temperature dependent skyhook control of HMMWV suspension using a failsafe magneto-rheological damper［C］．Pro．SPIE 5054，SmartStructuresandMaterials2003：Industrialand Commercial Applications of Smart Structures Technologies，332，San Dieo，United States，2003．

［5］Wilson N L，Wereley N M，Choi Y T，et al．Performance robustnessofamagneto-rheologicalseatsuspensionto temperature variations using skyhook control［C］．Proc．SPIE 7288，Active and Passive Smart Structures and Intelligent Systems，San Dieo，United States，2009．

［6］Kevorkian J，Cole J D，Perturbation methods in applied Mathematics［M］．Springer New York，1981．

［7］辛麗麗，梁繼輝，聞邦椿．考慮物料慢變參數的振動壓實系統分析［J］．東北大學學報，2009，30（4）：585－588．

XIN Li-li，LIANG Ji-hui，WEN Bang-chun．Analysis of vibratingcompactorsysteminvolvingslowlyvarying parameters of compacted material［J］．Journal of Northeastern University：Natural Science，2009，30（4）：585－588．

［8］陳予?。蔷€性振動［M］．天津：科學技術出版社，1983．

［9］李以農，鄭玲，聞邦椿．一類具有參數慢變的非線性振動系統［J］．重慶大學學報，2000，23（6）：16－18．

LI Yi-nong，ZHENG Ling，WEN Bang-chun．Study on the nonlinear vibration system with slow-changing parameters［J］．Journal of Chongqing University：Natural Science Edition，2000，23（6）：16－18．

［10］徐趙東，李愛群，程文瀼，等．磁流變阻尼器帶質量元素的溫度唯象模型［J］．工程力學，2005，22（2）：144－148．

XU Zhao-dong，LI Ai-qun，CHENG Wen-rang，et al．A temperature phenomenological model with mass element of magneto-rheological damper［J］．Engineering Mechanics，2005，22（2）：144－148．

［11］Sahin H，Wang X J，Gordaninejad F．A new model for yield stress of magneto-rheological greases／gels under combined effects of magnetic field and temperature［J］．Proc．SPIE 7288，Active and Passive Smart Structures and Intelligent Systems，San Dieo，United States，2009．

［12］SahinH，WangXJ，GordaninejadF．Temperature Dependence of Magneto-rheological Materials［J］．Journal of Intelligent Material Systems and Structures，2009：1－8．

［13］聞邦椿，李以農，韓清凱．非線性振動理論中的解析解法及工程應用［M］．沈陽：東北大學出版社，2001．

［14］李韶華，楊紹普．采用改進Bingham模型的非線性汽車懸架的主共振［J］．振動與沖擊，2006，25（4）：109－111．

LI Shao-hua，YANG Shao-pu．Primary resonance of a nonlinear vehicle suspension system using a modified Bingham model［J］．Journal of Vibration and Shock，2006，25（4）：109－111．

［15］郭樹起．磁流變阻尼器動態特征分析及其在車輛懸架中的應用［D］．北京：北京交通大學，2005．

［16］王四棋，余淼，浮潔，等．正弦激勵作用下磁流變阻尼器溫升理論與試驗研究［J］．機械工程學報，2013，49（8）：123 －128．

WANGSi-qi，YUMiao，FUJie，etal．Heatingof magnetorheological dampers under sine motion：theoretical andexperimentalstudy［J］．JournalofMechanical Engineering，2013，49（8）：123－128．

Main resonance analysis for a MR suspension system with slowly varying parameters

DONG Xiao-min，WANG Xiao-long，GU Xiao-ke
（State Key Laboratory of Mechanical Transmission，Chongqing University，Chongqing 400044，China）

The output damping force of a magneto-rheological（MR）damper is reduced with increase in working temperature．As a result，the performance of a MR suspension system is degraded．On the basis of the theory of heat transfer and nonlinear vibration，the main resonance of a MR suspension dynamic system with slowly varying parameters was analyzed here．The energy differential equations at the working area of the MR damper were established and solved with the numerical method．The modified Bingham model was applied to calculate the output damping force of the MR damper．The nonlinear dynamic equation of the MR suspension system with slowly varying parameters was formulated．The steady amplitude response of the MR suspension system under the excitation of main resonance was solved with the averaging method．The calculated results showed that the working temperature rises with increase in running time of the MR damper，consequently，the output damping force of the MR damper decreases significantly and the vibration amplitude variation of the suspension system increases，the performance of the MR suspension system is degraded．

slowly varying parameter；MR damper；temperature；nonlinear vibration

TH113．1

A

10．13465／j．cnki．jvs．2014．23．022

國家自然科學基金（51275539，60804018）；中央高?；穑–DJZR12110058，CDJZR13135553）

2013－09－12 修改稿收到日期：2013－12－19

董小閔男，博士，教授，1975年生