許飛鴻　徐趙東

(東南大學(xué)混凝土及預(yù)應(yīng)力混凝土結(jié)構(gòu)教育部重點實驗室，南京　210096)

?

多級線圈磁流變阻尼器磁路分析

許飛鴻徐趙東

(東南大學(xué)混凝土及預(yù)應(yīng)力混凝土結(jié)構(gòu)教育部重點實驗室，南京210096)

摘要:為了研究多級線圈磁流變阻尼器內(nèi)部的磁場分布，對多級線圈磁流變阻尼器的磁路進行了理論分析，根據(jù)磁路歐姆定律給出了多級線圈磁流變阻尼器有效間隙處磁感應(yīng)強度的計算方法．基于通用有限元軟件建立了多級線圈磁流變阻尼器的有限元模型，計算了不同磁場下各阻尼間隙處的磁感應(yīng)強度．通過引入高斯函數(shù)和指數(shù)函數(shù)，模擬了磁場的分布規(guī)律以及電流對磁場的影響規(guī)律，建立了多級線圈磁流變阻尼器磁場分布的計算模型，并給出了相應(yīng)計算公式．結(jié)果表明，該模型可以快速有效地計算出不同磁場下阻尼器各有效阻尼間隙處的磁感應(yīng)強度，模型計算結(jié)果與有限元仿真結(jié)果吻合良好．所提模型為多級線圈磁流變阻尼器阻尼力計算以及線圈通電優(yōu)化分析提供了基礎(chǔ)．

關(guān)鍵詞:多級線圈磁流變阻尼器;磁路;磁感應(yīng)強度;有限元仿真

引用本文:許飛鴻，徐趙東．多級線圈磁流變阻尼器磁路分析［J］．東南大學(xué)學(xué)報(自然科學(xué)版)，2016，46(1) : 100－104．DOI: 10．3969/j．issn．1001－0505．2016．01．017．

磁流變液(magnetorheological fluid，MRF)作為一種智能材料［1－2］，在磁場作用下可以瞬間由流動的液體轉(zhuǎn)變?yōu)榘牍虘B(tài)的黏塑性體，具有可控的屈服強度．基于MRF研制的磁流變阻尼器(magneto-rheological damper，MR damper)具有耗能低、出力大、響應(yīng)速度快等優(yōu)點，已廣泛應(yīng)用于土木、車輛等領(lǐng)域的振動控制中［3－6］．在對大型結(jié)構(gòu)進行振動控制時，阻尼器需要提供很大的阻尼力，通過增加勵磁線圈的匝數(shù)可以提高磁場強度，從而達到增大阻尼力的目的．然而，隨著匝數(shù)的增加，活塞長度也增加，阻尼間隙處的磁場不均勻度以及活塞運動的不穩(wěn)定性也將隨之增加．為了解決此問題，部分學(xué)者研制了多級線圈磁流變阻尼器［7－9］，即將多組線圈并聯(lián)連接，從而可以在提高阻尼力的同時提高磁場均勻度以及阻尼器的響應(yīng)時間［10－11］．

在計算磁流變阻尼器提供的阻尼力時，首先必須進行磁路分析，計算每個阻尼間隙處的磁感應(yīng)強度．多級線圈磁流變阻尼器的磁路結(jié)構(gòu)復(fù)雜，不同線圈產(chǎn)生的磁場相互疊加，與普通單線圈磁流變阻尼器的磁路大不相同．然而，現(xiàn)有文獻主要圍繞單線圈磁流變阻尼器進行磁路分析，關(guān)于多級線圈磁流變阻尼器的磁路分析則較少．鑒于此，本文以五級線圈磁流變阻尼器為例，對多級線圈磁流變阻尼器的磁路進行了理論分析與有限元仿真，并在有限元仿真結(jié)果的基礎(chǔ)上通過引入高斯函數(shù)，提出了五級線圈磁流變阻尼器的磁場分布計算模型．這種多級線圈磁流變阻尼器的磁路分析方法具有普適性，利用所提模型可以有效計算任意電流下各阻尼間隙處的磁感應(yīng)強度，并為線圈通電的優(yōu)化分析提供基礎(chǔ)．

1　理論分析

研究對象為本課題組研制的五級線圈磁流變阻尼器，其結(jié)構(gòu)如圖1所示．外缸筒材料選用45#鋼，活塞材料選用DT4電工純鐵;線圈和磁流變液之間填充環(huán)氧樹脂，防止漏磁以及線圈腐蝕．

圖1　磁流變阻尼器結(jié)構(gòu)示意圖

當(dāng)線圈通入電流時，磁導(dǎo)率大的導(dǎo)磁體將構(gòu)成磁路的通徑．該阻尼器的等效磁路如圖2所示．圖中，R1～R4分別為阻尼器兩端部分磁芯、磁軛、磁流變液以及鋼套筒的磁阻，R'1～R'4為各線圈之間公共部分的磁阻，Φ1～Φ5，Φ'1～Φ'4分別為磁路中各部分的磁通量．根據(jù)磁路歐姆定律可得［9］

圖2　磁流變阻尼器等效磁路

阻尼器各部分磁阻的大小與其材料、尺寸有關(guān)．該磁流變阻尼器為圓形，對于圓環(huán)形磁路單元，根據(jù)磁路基本理論［12］，可得環(huán)形軸向單元的磁阻為

式中，μ0為真空磁導(dǎo)率;μ為材料的相對磁導(dǎo)率; r1，r2分別為環(huán)形單元的內(nèi)徑與外徑; l為軸向單元的軸向長度; h為徑向單元的寬度．

根據(jù)式(1)～(3)，可計算得到區(qū)域k的磁通量，進而計算每個部分的磁感應(yīng)強度，即環(huán)形徑向單元的磁阻為

式中，Bk，Sk為區(qū)域k的磁感應(yīng)強度及面積．

在磁流變阻尼器的實際工作環(huán)境下，根據(jù)結(jié)構(gòu)每一時刻的動力響應(yīng)，采用最優(yōu)控制算法計算出每一時刻阻尼器所需提供的最優(yōu)控制力，而阻尼器輸出控制力的大小是由每一個阻尼間隙處的磁感應(yīng)強度決定的，因此，如何快速有效地確定不同通電情況下每個阻尼間隙處的磁感應(yīng)強度是非常重要的．磁流變阻尼器中各材料的磁化特性都是非線性的，相對磁導(dǎo)率并不是一個定值，故通過求解式(1)～(4)得到各部分的磁感應(yīng)強度較為繁瑣．鑒于此，本文對磁流變阻尼器的磁路進行了有限元仿真分析，并在仿真結(jié)果的基礎(chǔ)上，通過引入高斯函數(shù)來模擬磁場在不同阻尼間隙處的分布情況，最終提出了磁流變阻尼器在不同電流下的磁場分布計算模型．

2　有限元仿真

磁流變阻尼器內(nèi)部結(jié)構(gòu)尺寸如表1所示，采用通用有限元軟件對磁流變阻尼器進行仿真分析，根據(jù)結(jié)構(gòu)對稱性采用二維軸對稱模型，單元選用二維八節(jié)點四邊形實體單元PLANE53．忽略漏磁效應(yīng)，即磁導(dǎo)體內(nèi)部靠近邊界磁感應(yīng)強度方向平行于界面．圖3為磁流變阻尼器的幾何模型及有限元模型．

表1　磁流變阻尼器設(shè)計參數(shù)

圖3　磁流變阻尼器模型

將磁流變阻尼器內(nèi)部5個線圈和6個阻尼間隙分別進行編號(見圖1)．根據(jù)對稱性，分別對1#，2#，3#線圈施加電流，得到3個線圈分別通電時內(nèi)部磁場的分布情況，施加電流I = 0. 6，1. 2，1. 8，2. 1，2. 4 A．由于阻尼力的大小是由每個阻尼間隙處的平均磁感應(yīng)強度決定的，因此在以下分析中均采用每個阻尼間隙處的平均磁感應(yīng)強度．當(dāng)I = 0. 6 A時，1#，2#，3#線圈分別通電時磁場在每個阻尼間隙處的分布情況見圖4．

圖4　單個線圈分別通電時的磁場分布

3　磁場分布計算模型

式中，B0. 6表示I =0. 6 A時阻尼器阻尼間隙處的磁感應(yīng)強度; x為位置坐標(biāo)，以圖3中左端點(即阻尼器模型左部端點)為坐標(biāo)原點，向右為正方向; a，b，c為待定參數(shù)．

當(dāng)I =0. 6 A時，采用高斯函數(shù)對1#，2#，3#線圈分別通電時的磁場分布進行模擬，結(jié)果見圖5．由圖可知，利用高斯函數(shù)可以較好地模擬不同線圈通電時磁場在各阻尼間隙處的分布情況．

為了反映施加電流大小對阻尼間隙處磁感應(yīng)強度的影響，對I =0. 6，1. 2，1. 8，2. 1，2. 4 A時1#，2#，3#線圈分別通電時阻尼間隙處的磁感應(yīng)強度進行計算，發(fā)現(xiàn)對每個線圈單獨通電時，磁感應(yīng)強度隨電流的變化具有相同的規(guī)律，可以采用指數(shù)函數(shù)

由圖4可以看出，單個線圈通電時，磁場在阻尼器各阻尼間隙處的分布服從高斯分布，故可采用高斯函數(shù)來模擬磁場分布，計算公式為對電流的影響進行模擬，計算公式為

式中，f為電流對磁感應(yīng)強度的影響系數(shù); m，n，p，q分別為待定參數(shù)．

根據(jù)計算結(jié)果，采用指數(shù)函數(shù)對1#，2#，3#線圈分別通電時的f值進行擬合，結(jié)果見圖6．由圖可知，f(0. 6) =1;且3條曲線幾乎重合，表明電流的影響相同．

圖5　高斯函數(shù)擬合結(jié)果

圖6　電流對磁感應(yīng)強度的影響

結(jié)合式(5)和(6)，即可得到多級線圈磁流變阻尼器磁場分布的計算模型為

式中，B(I)為施加不同電流時的磁感應(yīng)強度．

利用式(7)便可計算得到任意線圈單獨施加電流時磁流變阻尼器每個阻尼間隙處的磁感應(yīng)強度值．然后，根據(jù)不同線圈產(chǎn)生磁場疊加的原理［9］，可以計算任意線圈組合通電時每個阻尼間隙處的磁感應(yīng)強度值，計算公式如下:

式中，Br為第r個阻尼間隙處總的磁感應(yīng)強度; Btr為第t個線圈在第r個阻尼間隙處產(chǎn)生的磁感應(yīng)強度，設(shè)磁場方向向上為正; It為第t個線圈所施加電流的大?。?/p>

圖7　3個線圈間隔通電時磁感應(yīng)強度分布圖

為了驗證本文所提計算模型的有效性，以3個線圈間隔通電情況(即磁流變阻尼器1#，3#，5#線圈分別施加相同電流)為例，將所提模型的計算值與有限元仿真值進行對比，結(jié)果如圖7所示．由圖可知出，計算結(jié)果與有限元仿真結(jié)果吻合良好，表明本文所提出的多級線圈磁流變阻尼器磁場分布計算模型是有效的．

4　結(jié)語

本文以五級線圈磁流變阻尼器為例，基于磁路歐姆定律，對多級線圈磁流變阻尼器的磁路進行了理論分析，給出了多級線圈磁流變阻尼器的阻尼間隙處磁感應(yīng)強度計算方法．對多級線圈磁流變阻尼器的磁路進行了有限元仿真，并在仿真結(jié)果的基礎(chǔ)上，通過引入高斯函數(shù)和指數(shù)函數(shù)，建立了多級線圈磁流變阻尼器內(nèi)部磁場分布的計算模型．模型計算結(jié)果與有限元仿真結(jié)果吻合良好，說明該模型可以快速有效地計算任意通電工況下各阻尼間隙處的磁感應(yīng)強度，為磁流變阻尼器通電優(yōu)化分析提供了基礎(chǔ)．

參考文獻(References)

［1］Ashour O，Rogers C A，Kordonsky W．Magnetorheological fluids: materials，characterization，and devices ［J］．Journal of Intelligent Material Systems and Structures，1996，7(2) : 123－130．

［2］Jolly M R，Bender J W，Carlson J D．Properties and applications of commercial magnetorheological fluids ［J］．Journal of Intelligent Material Systems and Structures，1999，10(1) : 5－13．

［3］Yao G Z，Yap F F，Chen G，et al．MR damper and its application for semi-active control of vehicle suspension system［J］．Mechatronics，2002，12(7) : 963－973．

［4］管友海，黃維平．MR阻尼器在海洋平臺半主動振動控制中的應(yīng)用［J］．中國海洋平臺，2002，17(3) : 25－28．DOI: 10. 3969/j．issn．1001－4500. 2002. 03. 007．Guan Youhai，Huang Weiping．The application of mr damper in the semi-active vibration control for offshore platform［J］．China Offshore Platform，2002，17(3) : 25－28．DOI: 10. 3969/j．issn．1001－4500. 2002. 03. 007．(in Chinese)

［5］Xu Z D，Shen Y P，Guo Y Q．Semi-active control of structures incorporated with magnetorheological dampers using neural networks［J］．Smart Materials and Structures，2003，12(1) : 80－87．

［6］Xu Z D，Sha L F，Zhang X C，et al．Design，performance test and analysis on magnetorheological damper for earthquake mitigation［J］．Structural Control and Health Monitoring，2013，20(6) : 956－970．

［7］Liao C R，Zhao D X，Xie L，et al．A design methodology for a magnetorheological fluid damper based on a multi-stage radial flow mode［J］．Smart Materials and Structures，2012，21(8) : 085005．

［8］Hu G，Long M，Huang M，et al．Design，analysis，prototyping，and experimental evaluation of an efficient double coil magnetorheological valve［J］．Advances in Mechanical Engineering，2014，6: 403410．

［9］Zheng J J，Li Z C，Koo J H，et al．Magnetic circuit design and multiphysics analysis of a novel MR damper for applications under high velocity［J］．Advances in Mechanical Engineering，2014，6: 402501．

［10］楊廣強，Spencer B F，Carlson Jr J D，et al．足尺磁流變阻尼器的建模及動態(tài)特性［J］．地震工程與工程振動，2001，21(4) : 8－23．Yang Guangqiang，Spencer B F，Carlson Jr J D，et al．Large-scale MR fluid dampers: modeling，and dynamic performance considerations［J］．Earthquake Engineering and Engineering Vibration，2001，21(4) : 8－23．(in Chinese)

［11］關(guān)新春，歐進萍．磁流變減振驅(qū)動器的響應(yīng)時間試驗與分析［J］．地震工程與工程振動，2002，22(6) : 96－102．DOI: 10. 3969/j．issn．1000－1301．2002．06．016．Guan Xinchun，Ou Jinping．Response time experiment and analysis of magnetorheological fluid actuator［J］．Earthquake Engineering and Engineering Vibration，2002，22(6) : 96－102．DOI: 10. 3969/j．issn．1000－1301. 2002. 06. 016．(in Chinese)

［12］周瑾，吳友海．支承磁懸浮軸承的磁流變阻尼器磁場分析［J］．磁性材料及器件，2010，41(6) : 17－20．DOI: 10. 3969/j．issn．1001－3830. 2010. 06. 005．Zhou Jin，Wu Youhai．The magnetic field analysis of magneto-rheological damper support for magnetic bearing［J］．Journal of Magnetic Materials and Devices，2010，41(6) : 17－20．DOI: 10. 3969/j．issn．1001－3830. 2010. 06. 005．(in Chinese)

Magnetic circuit analysis on multi-coil magnetorheological damper

Xu Feihong Xu Zhaodong
(Key Laboratory of Concrete and Prestressed Concrete Structures of Ministry of Education，Southeast University，Nanjing 210096，China)

Abstract:In order to study the magnetic field distribution of a multi-coil magnetorheological damper，the magnetic circuit of the multi-coil magnetorheological damper is theoretically analyzed，and the computational method for the magnetic induction intensity in the effective damping gaps of the multi-coil magnetorheological damper is proposed based on the magnetic circuit Ohm's law．The finite element model for the multi-coil magnetorheological damper is established by using the general finite element software，and the magnetic induction intensities of the effective damping gaps in different magnetic fields are calculated．The computational model for the magnetic field distribution of the multi-coil magnetorheological damper is established by using the Gaussian function and the exponential function to simulate the magnetic field distribution law and the influence of the current on the magnetic field，and the corresponding computational formulas are also given．The results show that the magnetic induction intensities of the damping gaps can be calculated effectively and quickly by using the proposed model in different magnetic fields，and the calculation results by the proposed model agree well with the simulation results by the finite element model．The proposed model lays a foundation for the damping force calculation and the optimal analysis of current inputting．

Key words:multi-coil magnetorheological damper; magnetic circuit; magnetic induction intensity; finite element simulation

基金項目:國家中青年科技創(chuàng)新領(lǐng)軍人才支持計劃資助項目、江蘇省杰出青年基金資助項目(BK20140025)、江蘇高校優(yōu)勢學(xué)科建設(shè)工程資助項目(CE02-2-33)、江蘇省普通高校研究生科研創(chuàng)新計劃資助項目(KYLX15(－)0088)．

收稿日期:2015-06-29．

作者簡介:許飛鴻(1989—)，男，博士生;徐趙東(聯(lián)系人)，男，博士，教授，博士生導(dǎo)師，xuzhdgyq@ seu．edu．cn．

DOI:10．3969/j．issn．1001－0505．2016．01．017

中圖分類號:TU352. 1

文獻標(biāo)志碼:A

文章編號:1001－0505(2016) 01-0100-05