# Theoretical research of mechanical behaviour of magneto-rheological fluid/Teorinis magnetoreologiniu skysciu tyrimas.

1. Introduction

Design of new features in machines and devices, reducing dimensions of equipment and wide use of computer control requires implementation of materials with special properties. In order to create smart computer-controlled drives with desired dynamic characteristics, smart materials became desired in many applications. One of successful computer controlled motion application is an implementation of magneto-rheological materials in such drives. Magneto-rheological fluids (MRF) are a class of materials whose rheological properties may be rapidly altered by applying a magnetic field [1].These materials usually are liquids with magnetically polarized ferrous particles suspended in a carrier liquid. The particle size in MRFs has on the order of micrometers and is 1-3 orders of magnitude larger than the same particles in ferromagnetic fluids [2]. Main interest to MRF is caused due to their feature to change its viscosity from magnitude of applied external magnetic field as well as liquid shape change. This feature makes MRF attractive to various commercial applications, such as dampers [3, 4] and brakes [4, 5]. This paper is intended to evaluate possibility to utilize MRF for use in computer controlled actuators.

2. Formulation of research using finite element method

Device designing with MRF requires estimation of behaviour of MRF in static and dynamic magnetic field and possibility to reach desired characteristics. Theoretical modelling of MRF is hardly possible due inability directly to solve non-linear differential equations. The way to do that is an application of finite element analysis. From the perspective of FEA methodology, there are necessary to perform such steps for FEA in case of MRF [6] (Fig. 1):

1. Building the geometry of object;

2. Building the finite element model;

3. Solving the finite element model;

4. Analysing the results.

Mathematical modelling starts with geometry of object (Fig. 1). 2D or 3D CAD programs are mostly used for geometrical modelling. CAD model should approximately, with some boundary simplification, correspond to modelled object.

CAD model is often sophisticated geometrically, therefore, in order to perform calculations using finite element method, it is idealized which helps to simplify its geometry and eliminate elements that do not affect calculations and to exchange thin panes of model to surfaces, etc. Idealization of a model takes place in order to simplify the task and make the model able to mesh as well as to shorten calculation time.

[FIGURE 1 OMITTED]

From magneto-rheological fluids documentation that provides Yield Stress--Magnetic Field Strength dependencies, we can calculate yield stress values in accordance to calculated magnetic field strength. Achieved values can be used for calculation of total pressure drop in a system.

3. The goal of the research

The main goal of an executed theoretical research is to determine pressure fluctuation change to magneto-rheological fluids in the closed volume under influence of external magnetic field, thus causing fluid chamber to change their shape and deliver displacement to the system. In order to reach the goal, the following tasks have been set: to determine how external magnetic field influences general magneto-rheological fluid's pressure fluctuation; to determine parameters that define magneto-rheological fluid's property to change its state when exposed to a magnetic field.

4. Research methodology

The pressure drop developed in a device based on pressure driven flow mode is commonly assumed to result from the sum of a viscous component [DELTA][P.sub.[eta]] and a field dependent induced yield stress component [DELTA][P.sub.[tau]], and can by approximated by [7]:

[DELTA]P = [DELTA][P.sub.[eta]] + [DELTA][P.sub.[tau]]H = [12[eta]QL/[g.sup.3]w] + [c[[tau].sub.y](H)L/g], (1)

where L is the length, g is the gap between fluid walls, w is the width of the flow channel between the fixed poles, Q is the volumetric flow rate, [eta] is the viscosity with no applied field, [[tau].sub.y] is the yield stress developed in response to an applied field H. The parameter c is a function of the flow velocity profile and has a value ranging from 2 (for [DELTA][P.sub.[eta]] / [DELTA]P[tau] less than ~1) to a maximum value of 3 (for [DELTA][P.sub.[eta]] / [DELTA]P[tau] less than ~100) [8].

Behavior of controllable fluids is often represented as a Bingham plastic having variable yield strength [7]:

[[tau].sub.total] = [[tau].sub.y](H) + [[eta].sub.p][??], (2)

where [[tau].sub.total] is the total yield stress, [tau](H) is the strength caused by the applied magnetic field H, [[eta].sub.p] is the magnetic field--independent plastic viscosity defined as the slope of the shear stress versus shear strain rate relationship, [??] is the shear rate.

Below the yield stress, material behaves viscoelastically [7]:

[[tau].sub.totoal] = G[gamma] ; [[tau].sub.totoal] < [[tau].sub.y] (3)

where G is the complex material modulus.

From the formula (1) it is evident, that in order to determine how external magnetic field influences magneto-rheological fluid's pressure loss, it is needed to determine yield stress [[tau].sub.y]. In majority of magneto-rheological fluid's paperwork there is graphs presenting yield stress vs. magnetic field strength, therefore, once magnetorheological fluid's magnetic field strength is calculated, yield stress figures are available.

To determine magneto-rheological fluid's and Magnetic field Strength (H) a Femm 4.2 program was used. The program enables solving low frequency electromagnetic problems on two-dimensional planar and axisymmetric domains [9].

Femm 4.2 program has three modules:

1. Femm.exe program is designed for creation of geometry, defining of physical properties of materials and boundary conditions of model;

2. Triangle.exe program splits the analysed mathematical model into triangle elements;

3. Fkern.exe is a solver that solves differential equations in order to obtain a solution.

In magnetostatic problems time is invariant, therefore [9]:

[??] x H = J; (4)

[??] x B = 0, (5)

where H is field intensity, J is divergence, B is flux density.

Connection between B and H in each material is expressed as follows [9]:

B = [mu]H, (6)

where [mu] is permeability.

For nonlinear materials, the permeability is function of B [9]:

[mu] = B/H(B). (7)

Femm goes about finding a field that satisfies (4) and (6) via magnetic vector potential approach. Flux density is written in terms of the vector potential, A, as [9]:

B = [??] x A, (8)

where A is vector potential.

This definition of B always satisfies (5). The, (4) can be rewritten as:

[??] x ([1/[mu](B)] [??] x A) = J. (9)

For linear isotropic materials (assuming the Coulomb gauge, [??] x A = 0), (9) reduces to:

-[1/[mu]][[??].sup.2] A = J. (10)

Femm retains the form (9), so that magnetostatic problems with a nonlinear B - H relationship can be solved.

To execute the research with Femm.exe program there was a two-dimensional solenoid model created and above the solenoid there were MRF fluids situated in an aluminium box.

[FIGURE 2 OMITTED]

Internal height of the box is 1 mm. Magneto-rheological fluid is isolated from external environment by 1 mm aluminium walls. In the research solenoid, MRF and electromagnetic forces were used that were arranged symmetrically and vertically to the axis, therefore the research is axisymmetric meaning that there was a half of the model used for the research. Accurate model dimensions are provided in a Fig. 2. To strengthen magnetic field of solenoid there was an iron core used which had 2000 copper coils of 0.63 mm in diameter winded around iron core.

The analysed model with Triangle.exe program was split into 11852 elements. In total 6133 nodes were generated. Elements number of finite elements mesh was enlarged in those areas, where measures of the model were smaller (Fig. 3).

[FIGURE 3 OMITTED]

In the research there were four Lord Corporation magneto-reological fluids of different characteristics used that are displayed in documentation [10-13].

The research was executed under changing strength of solenoid's flow from 0 to 2 A. After each experiment an average magnetic field strength H and magnetic flux density B, crossing notional horizontal fluid's symmetry axis was calculated.

After comparing magnetic lines arrangement in Figs. 4 and 5 that are formed under influence of electric streams flowing through solenoid's strands, we can see that magneto-rheological fluids distorts magnetic field's lines of force (Fig. 6).

[FIGURE 4 OMITTED]

[FIGURE 5 OMITTED]

[FIGURE 6 OMITTED]

Generated calculations of magnetic field strength and magnetic flux density figures dependency from strength of the stream flowing through strands of solenoid are presented in Tables 1 and 2 as well as Figs. 7 and 8.

[FIGURE 7 OMITTED]

The Fig. 8 shows, that in range of current strength from 0 to 2 A, magnetic induction changes almost linearly.

From researched magneto-rheological fluids documentation that provides Yield Stress - Magnetic Field Strength dependencies, we can take the figures of yield stress in accordance to the figures of magnetic field strength.

[FIGURE 8 OMITTED]

The results are presented in Table 3 and Fig. 9. Additional pressure generated by magnetic field strength causes change of magneto-rheological fluid surface, therefore the system can be treated as mechanical.

[FIGURE 9 OMITTED]

4. Conclusions

After theoretical research of MRF using FEM and modified algorithm of internal pressure definition, some conclusions are made:

1. There is created new methodology, which enables to find mechanical behavior of MRF in magnetic field.

2. Internal pressure from magnetic field in MRF can be obtained for static and dynamic magnetic load case, movement of liquid pressurized membrane requires special CFD analysis.

Accepted September 05, 2013

References

[1.] Schwartz, M. 2001. Encyclopedia of Smart Materials Set, 1st ed. Wiley-Interscience, 2259 p.

[2.] Li, D. 2008. Encyclopedia of Microfluidics and Nanofluidics, 1st ed. Springer, 2226 p. http://dx.doi.org/10.1007/978-0-387-48998-8.

[3.] Odenbach, S. 2009. Colloidal Magnetic Fluids: Basics, Development and Application of Ferrofluids, Softcover reprint of hardcover 1st ed. 2009 ed. 440 p.

[4.] Odenbach, S. 2002. Magnetoviscous Effects in Ferrofluids, 1st ed. Springer, 151 p.

[5.] Varadan, V.K.; Vinoy, K.J.; Gopalakrishnan, S. 2006. Smart Material Systems and MEMS: Design and Development Methodologies, 1st ed. Wiley, 418 p. http://dx.doi.org/10.1002/0470093633.

[6.] Dassault Systemes. 2008. S.A. COSMOSWorks Designer, Concord--Massachusetts, 394 p.

[7.] Phillips, R.W. 1969. Engineering applications of fluids with a variable yield stress, Ph.D. Thesis, Univeristy of California, Berkeley.

[8.] Jolly, M.; Bender, J.; Carlson, J. 1999. Properties and applications of commercial magnetorheological fluids, Journal of Intelligent Material Systems and Structures, 10(1): 5-13.

[9.] Meeker, D. 2010. Finite Element Method Magnetics, Femm 4.2 User's Manual, 158 p.

[10.] MRF-336AG Magneto-Rheological Fluid <http:/ www.lord.com/products-and-solutions/magneto-rheo logical-%28mr%29/product.xml/1640/2> Accessed 2012 Feb 3.

[11.] MRF-241ES Magneto-Rheological Fluid <http:// www.lord.com/products-and-solutions/magneto-rheo logical%28mr%29/product.xml/1640/2> Accessed 2012 Feb 3.

[12.] MRF-132EG Magneto-Rheological Fluid <http:// www.lord.com/products-and-solutions/magneto-rheo logical-%28mr%29/product.xml/1645/2> Accessed 2012 Feb 3.

[13.] MRF-140CG Magneto-Rheological Fluid <http:// www.lord.com/products-and-solutions/magnetorheological-%28mr%29/product.xml/1646/2> Accessed 2012 Feb.

A. Klevinskis, Vilnius Gediminas Technical University, Basanaviciaus 28, 03224 Vilnius, Lithuania, E-mail: andrius.klevinskis@dok.vgtu.lt

V. Bucinskas, Vilnius Gediminas Technical University, Basanaviciaus 28, 03224 Vilnius, Lithuania, E-mail: vytautas.bucinskas@vgtu.lt

D. Udris, Vilnius Gediminas Technical University, Naugarduko 41, 03227 Vilnius, Lithuania, E-mail: dainius.udris@vgtu.lt

http://dx.doi.org/ 10.5755/j01.mech.19.5.5527
```Table 1
Calculated magnetic field strength values
of different magneto-rheological fluids

Current,   MRF-132AD,   MRF-140CG,   MRF-241ES,   MRF-336AG,   Air,
A          kA/m         kA/m         kA/m         kA/m         kA/m

0.2        2.0          1.4          1.8          1.9          3.0
0.5        5.2          3.6          4.5          4.7          7.7
1          10.4         7.4          9.4          9.6          15.5
1.5        15.6         11.3         14.8         14.9         23.2
2          20.9         15.7         20.8         20.6         30.9

Table 2
Calculated magnetic induction of different
magneto-rheological fluids

Curent,   MRF-132AD,   MRF-140CG,   MRF-241ES,   MRF-336AG,   Air,
A         T            T            T            T            T

0.2       0.028        0.030        0.029        0.029        0.023
0.5       0.071        0.074        0.073        0.072        0.055
1         0.142        0.148        0.144        0.144        0.11
1.5       0.213        0.221        0.214        0.214        0.172
2         0.283        0.292        0.282        0.283        0.230

Table 3
Current vs. yield stress

A          Pa           Pa           Pa           Pa

0.2        0            0            0            0
0.5        2.4          2.7          3            1
1          3.6          4            6            1.5
1.5        4.5          6            8            2.5
2          6            8            9            4.7
```