基于FLUENT的磁流變阻尼器阻尼特性分析

吳俊岐，侯保林

(南京理工大學 機械工程學院， 南京 210094)

【機械制造與檢測技術】

基于FLUENT的磁流變阻尼器阻尼特性分析

吳俊岐，侯保林

(南京理工大學 機械工程學院， 南京 210094)

以某新型磁流變阻尼器為研究對象，利用前處理軟件ICEM-CFD和FLUENT建立阻尼器內部空腔的流場模型，針對FLUENT無法直接模擬流動過程中材料性質變化，提出利用UDF來對仿真計算中使用的Herschel-Bulkley黏塑性流體模型進行定義的計算方法，實現牛頓流與黏塑性流的混合仿真計算。通過對比不同磁通密度以及不同活塞桿運動速度條件下阻尼器的阻尼力輸出，分析磁流變阻尼器的阻尼特性，得到磁流變液對阻尼器阻尼力輸出的影響規律。

磁流變阻尼器；FLUENT；阻尼特性分析；建模仿真

磁流變液(Magnetorheological Fluids,MRF)是一種新型的功能型材料，是由微米級的磁化顆粒混合于某種低磁導率的液體中制備而成。主要成分包括：軟磁性顆粒、基載液以及為了防止顆粒沉降而包覆在顆粒表面的添加劑，磁流變液可以在一定的磁場條件下產生磁流變效應。磁流變效應是在無磁場作用下，磁流變液呈牛頓流體的特性，而在強磁場的作用下，其流變特性發生急劇的變化，表現為類似固體的性質，撤去磁場時又恢復其流動特性的現象[1]。整個磁流變效應的響應過程非常迅速，在毫秒級。整個變化過程中液體的黏度連續，變化過程可逆且易于控制。基于上述優點，磁流變液在土木建筑，航空航天，軍事裝備等領域的應用前景十分廣闊。目前國內國外對磁流變液都進行了大量的研究，因此磁流體的生產也取得了較大的進展。

磁流變阻尼器是一種以磁流變液為基液的新型阻尼器，可以便于實現阻尼力的實時連續控制。傳統的阻尼器基液為普通液體，因此阻尼器所能承受的阻尼力有限且不可控，如傳統的火炮后坐阻尼器，當阻尼器的各主要結構尺寸確定后，阻尼器的阻尼力輸出也將確定，此時火炮的后坐行程與后坐速度都是近似固定值，磁流變阻尼器的出現改變了這一現狀，并且可以基于現有的阻尼器進行改進，無需完全重新設計。本文以某新型磁流變阻尼器為研究對象，提出了以利用FLUENT為仿真軟件，結合使用UDF(User-Defined Function)編寫的第三方程序為控制程序來模擬磁流變阻尼器的內部牛頓流與黏塑性流的計算方法。

1 磁流變阻尼器的工作原理

磁流變阻尼器的內部磁場由線圈環繞軟磁材料而成的電磁鐵提供。線圈環繞的區域外側，無磁流變效應，這部分區域成為非激活區；而在阻尼間隙的其他區域，磁場方向垂直于液體流動方向，會產生磁流變效應，可以通過控制線圈中電流的大小來改變磁流變效應程度，稱這部分區域為激活區。為了有效地改變阻尼力的大小，在阻尼器的中間軟磁體和內軟磁體上從左到右各布置了兩股線圈。

圖1是磁流變阻尼器的簡圖。阻尼器在受到外部作用力的狀態下，活塞桿帶動磁流變閥向左移動，內部的磁流變液通過阻尼間隙向右移動。流經阻尼間隙的流變區域的過程中，磁流變液受到磁場的作用，表觀上體現出的變化是黏度增大，內部剪切力增大，會瞬間轉變為接近固體狀態的黏塑性體[2]，隨著流體壓力的持續增大，磁流變體會克服材料的剪切屈服應力點繼續流動，流出磁流變閥后，阻尼間隙中的磁流變液恢復牛頓流體狀態。

1.阻尼器外筒; 2.外軟磁體; 3.外阻尼間隙; 4.中間軟磁體; 5.內阻尼間隙; 6.線圈; 7.內軟磁體

圖1 阻尼器簡圖

2 磁流變液動力學模型研究

關于磁流變阻尼器的分析與仿真計算，最根本的還是基于磁流變液的黏度、屈服極限、臨界應變率等材料特性的計算。磁流變液在磁場的作用下，內部的軟磁性材料被磁化，形成磁偶極子，彼此間相互吸引，從而形成抗剪力。在受到外磁場作用的前提下，磁流變液的剪切應力會隨著剪切應變的變化而變化，整個過程分為漸變區、屈服區、穩定區3個階段[3-4]。磁流變液的材料特性在每個區段都有各自的特點。

常用磁流變阻尼器動力學模型有Bingham模型、Herschel-Bulkley模型、非線性雙黏度模型、修正的Bingham模型等等。Bingham模型簡單易懂，現有的許多文獻采用的多是Bingham模型描述磁流變液在磁場作用下的行為，但實際情況下，磁流變液在剪切率變化的過程中會出現剪切致稠或稀化現象，Bingham無法描述這一現象。Herschel-Bulkley模型則可以模擬上述現象，因此使用 Herschel-Bulkley模型進行計算(如圖2所示)。在模型中,引入一個表征磁流變體剪切稀化的系數因子,在不同的外加載荷下具有不同的系數。

Herschel-Bulkley模型的本構方程為[5-6]

(1)

(2)

由此可知材料黏度是隨著剪切應變率冪次方的變化而變化的，此時Herschel-Bulkley模型的本構方程可表示為

(3)

當n=1時，Herschel-Bulkley模型就退化成了Bingham模型，nlt;1表示剪切稀化，ngt;1表示剪切致稠[7-8]。

為便于計算，流體黏度直接使用ANSYS/FLUENT的Herschel-Bulkley黏度表達式，查閱ANSYS/FLUENT的幫助文件，得到Herschel-Bulkley模型黏度變化規律定義如下：

(4)

（1）PoW。尋找隨機數n，使得：SHA(SHA(v||hp||t||n||hm))

τ0(B)=-86.51B3+122.56B2+17.42B-0.877

(5)

由式(5)可以看出，剪切屈服應力大小與磁場強度B有關，相同的磁流體前提下，通過控制磁場強度，即線圈的電流強度就可以控制黏度和剪切力的變化，進而控制阻尼器的阻尼力。

3 阻尼器阻尼力計算

3.1 流場的建模分析

首先利用三維建模軟件對流場區域進行建模，之后用ICEM-CFD(Integrated Computer Engineering and Manufacturing code for Computational Fluid Dynamics)對流場區域進行網格劃分，采用三維建模方法，網格全部采用六面體的非結構化網格。由圖1可以看出，整個流場區域的激活區與非激活區是連接在一起的，磁流變液在流經激活區與非激活區的過程當中液體的黏度以及屈服應力是變化的，由于無法通過FLUENT直接模擬這一現象，因此考慮使用UDF控制流體的黏度和剪切屈服應力的變化，核心思想是在ICEM中將激活區和非激活區定義為不同的part，之后通過UDF讀取part的ID進而對不同區域賦予不同的材料屬性，如圖3所示是UDF的程序流程框圖。由于UDF是基于C語言的一種開發環境，因此程序的格式與C語言程序格式完全相同，程序中用到的指令都是UDF的庫文件命令，將程序掛載到FLUENT中之后編譯，在設置材料黏度的時候直接使用。

圖3 UDF程序流程

本文用流動模式模擬剪切模式下流體的力學特征[9]，由于仿真計算中用到的是不可壓縮黏性流體，因此求解域的入口速度即為阻尼器活塞桿的運動速度；出口處的相對壓力設置為0，壁面為固定壁面。

3.2 阻尼力計算分析

關于磁流變阻尼器的分析與仿真計算，最根本的還是基于磁流變液的黏度、屈服極限、臨界剪切應變率等材料特性的計算。磁流變液在磁形將網格導入FLUENT中進行計算，圖4所示的是計算所用材料的屈服應力—磁通密度關系曲線，根據對磁流變閥處的分析計算(見圖5)，可以得到激活區中的最大磁感應強度約為0.471T，利用式(5)計算出對應的入口速度以及不同的磁通密度(不大于0.5T)來計算對應的流場的壓力，進而得到磁流變阻尼器的阻尼特性。

圖4 MRF-132LD磁流變液屈服應力與磁通密度關系

圖5 磁流變閥內電流與磁通密度關系

3.3 計算結果

可以看出，流體黏度的變化是按照程序的設定而分布的，在非激活區中黏度為0.3 Pa·s，在激活區中則是按照Herschel-Bulkley模型規律變化。再由計算的壓力數據，最終得到的不同磁場情況下磁流變阻尼器在不同活塞速度下所產生的阻尼力，如圖7所示。

圖6 流場黏度分布

圖7 阻尼器阻尼力與活塞速度關系

從計算結果可以看出，在速度不大于2 m/s的情況下，磁流變效應對阻尼器的阻尼力影響較大，電流從0 A增大到2 A可以使阻尼力翻倍增長，活塞桿速度為2 m/s時，0.5T時的阻尼力比0T時的阻尼力增加了124.3%，這是因為在速度較小的情況下，流體的黏塑性力大于慣性力，因此黏度的變化可以極大地影響阻尼力的變化[10]。當速度超過2 m/s時，流體的慣性力增大，因此調節黏度的大小對阻尼力的影響有限。當活塞桿的速度為10 m/s時，0.5T時的阻尼力比0T時的阻尼力增加9.8%，阻尼器的可控性依舊很明顯。因此，在受到相同外力條件下，磁流變阻尼器活塞桿的移動距離是可調的。

用計算所用阻尼器來承受火炮后坐力，火炮的后坐運動方程如下[11-12]：

(6)

其中：mh為后坐部分質量；x為后坐行程；t為后坐時間；Fpt為炮膛合力，同一彈丸有相同的Fpt；FMR為阻尼器阻力；Ff為復進機力。由式(6)可以看出，對于同一彈丸，使用相同火炮發射時，改變磁流變阻尼器中電流的大小就可以改變火炮的后坐行程與后坐速度。

4 結論

仿真結果表明，在一定范圍內(磁通密度飽和前)，通過改變線圈內電流的大小改變磁流變閥內部的磁場強度，可以有效改變磁流變液的剪切應力，達到控制磁流變阻尼器阻尼力的目的。

UDF程序不只針對文中使用的模型而專門編寫的，它也適用于其他阻尼器模型。針對不同的磁流變液以及阻尼器的結構，通過修改程序的參數就能夠進行仿真。

[1] 張掙鑫,劉黔會,黃方林.一種在ANSYS中模擬磁流變液的方法[J].貴州大學學報:自然科學版，2011,28(4)：110-113.

[2] 張紅輝,廖昌榮,余淼,等.考慮磁場的磁流變效應Bingham建模[J].機械工程學報,2006,42(2):33-36.

[3] 姚金光,晏華.高性能磁流變液研究的進展[J].材料開發與應用,2009,24(2):62-67.

[4] 李鴿鴿,李鴿鴿,劉保國,等.磁流變體在機械工程中的應用綜述[J].機械強度,2015(2):219-225.

[5] 李占衛,李治軍.磁流變阻尼器動力學模型的研究現狀[J].機械制造與自動化,2012,41(1):142-145.

[6] GUO N Q,DU H,LI W H.Finite element analysis and simulation evaluation of a magnetorheological valve[J].International Journal of Advanced Manufacturing Technology,2003,21(6):438-445.

[7] JOHN S,WERELEY N M.Nondimensional quasi-steady analysis of magnetorheological dampers utilizing a Herschel-Bulkley model with preyield viscosity[C]//Proceedings of SPIE-The International Society for Optical Engineering,2003:53-65.

[8] LEE D Y,CHOI Y T,WERELEY N M.Performance analysis of ER/MR impact damper systems using herschel-bulkley model[J].Journal of Intelligent Material Systems amp; Structures,2002,13(7/8):525-531.

[9] 侯保林,MEHDI,Ahmadian.沖擊載荷作用下磁流變阻尼器的建模與分析[J].機械工程學報,2006,42(4):173-178.

[10] 安德森，姚朝暉,周強.計算流體力學入門:Computational fluid dynamics:the basics with applications[M].北京：清華大學出版社,2010.

[11] 高躍飛.火炮反后坐裝置設計[M].北京：國防工業出版社,2010.

[12] 李延成.沖擊載荷下磁流變緩沖器半主動控制研究[D].南京：南京理工大學,2007.

(責任編輯唐定國)

AnalysisofDampingCharacteristicsofMagnetorheological(MR)DamperBasedonFLUENT

WU Junqi, HOU Baolin

(School of Mechanical Engineering, Nanjing University of Science and Technology, Nanjing 210094, China)

A new type of MR damper is taken as the research object, by using preprocessing software ICEM-CFD and FLUENT, it establishes the flow field model of the internal cavity of the damper. For the fact that FLUENT simulation can not directly change the properties of the material during flow process, it proposed the use of UDF to define the Herschel-Bulkley viscoplastic fluid model used in simulation calculations. And then the hybrid simulation calculation of Newton flow and viscoplastic flow is realized. By comparing the different magnetic flux densities and the damping force output of the dampers under different piston rod motion conditions. We analyzed the damping characteristics of MR damper, and the effect rule of MR fluid on the damping force output of the damper was obtained.

magnetorheological damper; FLUENT; damping characteristics analysis; modeling simulation

2017-07-02；

2017-07-31

國防973項目子課題(6132490102)

吳俊岐(1992—)，男，碩士研究生，主要從事機械設計與仿真研究；侯保林(1965—)，男，教授，博士生導師，主要從事火炮、自動武器、彈藥工程研究。

10.11809/scbgxb2017.11.031

本文引用格式：吳俊岐，侯保林.基于FLUENT的磁流變阻尼器阻尼特性分析[J].兵器裝備工程學報，2017(11):142-145.

formatWU Junqi, HOU Baolin.Analysis of Damping Characteristics of Magnetorheological(MR) Damper Based on FLUENT[J].Journal of Ordnance Equipment Engineering,2017(11):142-145.

TH113

A

2096-2304(2017)11-0142-04