# Simulation Study on the Motion of Magnetic Particles in Silicone Rubber-Based Magnetorheological Elastomers

Authors: Zhiqiang Xu; Heng Wu; Qiuliang Wang; Liyin Yi; Jun Wang
Journal: Mathematical Problems in Engineering
 (2019)

---

## Abstract

Magnetorheological elastomer (MRE) is an intelligent composite material and has been widely used in various fields such as vibration reduction and sensing. MRE has an excellent magnetorheological effect through the chaining of its internal magnetic particles. Current studies on MREs mainly focus on the preparation of materials and characterization of mechanical properties. However, very few studies have been conducted on the mechanism of magnetic particle motion during MRE curing. Based on the silicone rubber-based MRE, the motion mechanism of magnetic particles during curing was explored through numerical simulation. First, we analyzed the magnetic force and viscous force of magnetic particles in MRE and discussed the equations of motion of magnetic particles under applied magnetic field. Further, we established a uniform magnetic field model through the finite element method and simulated the motion of two magnetic particles under the magnetic field. Finally, we discussed the effects of particle distribution angles, particle radii, applied magnetic field strength, and distance between particles on particle velocity and displacement. The results show that the distance between particles has the greatest influence on the motion of magnetic particles, and the size of the distance between particles will affect the contact time of the particles, thus affecting the chain formation of magnetic particles in the MRE.


---

## Body

## 1. Introduction

Magnetorheological elastomer (MRE) is a smart composite material that consists of a polymeric matrix with embedded micron-sized magnetic particles, then solidifying to obtain an elastic material with variable stiffness under an applied magnetic field [1, 2]. MRE overcomes the shortcomings of magnetorheological fluid (MRF), such as leakage, easy sedimentation, and instability [3–5]. MRE has a wide range of engineering applications because its variable stiffness characteristics, such as tuning dampers, dampers, sensors, and magnetorheological elastic polishing bodies [6–9]. The excellent mechanical properties of MRE are determined by the internal magnetic particles, which will form a magnetic particle chain parallel to the direction of the magnetic field under applied magnetic field [10]. The anisotropic MRE with different mechanical properties was finally obtained by controlling the strength of the applied magnetic field, and the better the arrangement of particles, the better the mechanical properties of MRE [11]. Therefore, in order to improve the mechanical properties of MRE and make it more widely used in engineering applications, further study on its motion of magnetic particles is essential.At present, the research on magnetic particles in MRE mainly focuses on the types of magnetic particles, particle size, concentration of particles, and the influence of the applied magnetic field on the magnetic particle chain [12–14]. However, there are few studies on the chain formation of magnetic particles in MRE under magnetic field, the existing research generally analyzes the formation of chain by magnetic dipole theory [15, 16] or study on the final arrangement of magnetic particle chains by changing the influence of different external conditions [17], and there is less quantitative analysis of magnetic particle motion during MRE curing. In addition, the magnetic particle motion under the applied magnetic field was difficult to observe because of the limitation of particle size, so the motion of the particle was often described by numerical simulation [18–21].Based on the existing research, the equation of motion of magnetic particles in silicone rubber has been established. The motion of particles under magnetic field was simulated by the finite element method. The formation of magnetic particle chains was explained by the change of velocity and displacement during particle motion, and the effects of two magnetic particles at different distribution angles, particle radius, applied magnetic field, and different distance between the particles on motion process were discussed. Meanwhile, the contact time of magnetic particles under different conditions was analyzed, and the results show that changing the distance between two magnetic particles under the same ratio significantly increases the efficiency of particle chain formation. Moreover, the larger the distance between the particles, the longer the contact time, which provides an effective theoretical basis for preparing a silicone rubber-based MRE with excellent performance.


## 2. Theoretical Model

### 2.1. Movement Equation

In the preparation process of MRE, anisotropic MRE and isotropic MRE are obtained by controlling the presence or absence of a magnetic field in the curing stage [22]. The magnetic particles in the isotropic MRE are uniformly distributed in the matrix, and the magnetic particles inside the anisotropic MRE form a special chain structure under the magnetic field. Figure 1 is a schematic diagram showing the distribution of magnetic particles of an anisotropic MRE during the curing stage.Figure 1
Schematic diagram of chain formation under magnetic field in the curing stage (a) before curing and (b) after curing under magnetic field.
(a) (b)Analysis of the force of the magnetic particles in the MRE during the curing stage mainly includes the magnetic force by the applied magnetic field, the viscous force by the silicone rubber matrix, buoyancy, gravity, the interaction force between the particles, and the Brownian force [23]. The magnetic particles produce motion under the combined action of these forces, so the following equation is obtained according to Newton's formula:(1)midv→idt=F→mi+F→vi+F→ri+F→gi+F→bi+F→Biwheremi is the mass of magnetic particle i,v→i denotes the particle velocity, F→mi is the applied magnetic force, F→vi is the viscous force, F→ri is the interaction force between particles, and F→gi,F→bi, and F→Bi represent the gravity, buoyancy, and Brownian forces, respectively.The magnetic particles used in this study are micron-sized particles; thus, gravity and buoyancy can be considered to be balanced; thus, they are not considered. Further, the magnetic and viscous forces of micron-sized particles in incompressible Newtonian liquids are 200 times that of other forces[18]. Therefore, to facilitate calculation for this, the last four terms in (1) are omitted, and only the magnetic and viscous forces of particles moving in the silicone rubber matrix are considered, and the equation of movement of particles can be simplified to (2)midv→idt=F→mi+F→viIn this case, the ordinary differential equation of displacement can be expressed as follows:(3)u→¨=F→mi-F→vimiThe ordinary differential equation of velocity can also be expressed as follows:(4)v→i˙=F→mi-F→vimi(5)u→i˙=v→iwhere u→i and v→i represent the velocity and displacement of particle movement, respectively.


### 2.2. Magnetic Force

The magnetic force on the magnetic particles can be divided into two aspects: one is the force produced by the magnetization of particles under the applied magnetic field and the other is the force between the magnetic dipoles. Assuming that the magnetic particles in the magnetic field are spherical and uniform in size, the magnetic moment of the magnetic particle can be expressed as follows:(6)m→=VM→=43πR3χH→whereV is the volume of a single particle,V=(4pi∗R∧3)/3, R is the radius, M→ represents magnetization, M→=χH→, χ is the magnetic susceptibility, and H→ is the applied magnetic field strength.The magnetic force between the magnetic particles i and j is as follows:(7)F→mij=m→·∇B→(8)=∑j≠iμ0m→im→j4πrij3d→i·d→j-3d→i·r^ijd→j·r^ijThe magnetic force acting on particle i in the applied magnetic field is as follows:(9)F→mi=μ0m→iH→d→i·h^Based on the above equation, it is concluded that the comprehensive magnetic force of particle i is as follows:(10)F→mi=F→mij+F→mi(11)F→mi=3μ0m→im→j4πrij4∑j≠id→i·d→jr^ij+d→i·r^ijd→j+d→j·r^ijd→i-5d→i·r^ijd→j·r^ijr^ij+μ0m→iH→d→i·h^When only two magnetic particles act, the motion of the particles is a two-dimensional plane motion, so (11) is rewritten:(12)F→mi=4πμ0χ2H2R63rij41-5cos2⁡θijr^ij+2cos⁡θijh^Then the magnetic field force F→mi is decomposed in the x and y directions and obtained its matrix expression:(13)F→mxiF→myi=4πμ0χ2H2R63rij41-5cos2⁡θijsin⁡θij3-5cos2⁡θijcos⁡θijSince the direction of the magnetic force of the particle is related to the angle in the equation (13), so let 1-5cos2⁡θijsin⁡θij=0 and get θij=π/2 or θij=arccos⁡(15/5). Therefore, if the particles are to form a particle chain parallel to the direction of the magnetic field, the angle between the two particles θij<arccos⁡5/5 is required to be satisfied, at this time, and particles appear as close to each. If θij>arccos⁡5/5, the magnetic particles repel each other.As shown in Figure2, μ0 is the vacuum permeability; m→i and m→j are the particle magnetic moments; d→i and d→j are the direction vectors of the magnetic moment of particles i and j; r^ij is the direction vectors from particle i to particle j; H→ is the direction vector of the applied magnetic field intensity; h^ is the direction vector of the applied magnetic field intensity; rij denotes the distance between the centers of two particles i and j and satisfies the following requirement:(14)rij=r^ijFigure 2
Force diagram of particles i and j in a 2D plane under applied magnetic field.


### 2.3. Viscous Force

In a static incompressible fluid, the magnetic particles are subjected to the viscous forceF→v of the matrix when they move in the polymer matrix. The particle size is less than 3 mm; thus, according to the Stokes [24], (15)F→v=-6πRηdrdtR is the radius of the particle and η is the dynamic viscosity. The selected silicone rubber is purchased from Dow Corning, USA, and the viscosity is 3.5 Pa·s at room temperature. dr/dt is the velocity of the particle motion.


## 2.1. Movement Equation

In the preparation process of MRE, anisotropic MRE and isotropic MRE are obtained by controlling the presence or absence of a magnetic field in the curing stage [22]. The magnetic particles in the isotropic MRE are uniformly distributed in the matrix, and the magnetic particles inside the anisotropic MRE form a special chain structure under the magnetic field. Figure 1 is a schematic diagram showing the distribution of magnetic particles of an anisotropic MRE during the curing stage.Figure 1
Schematic diagram of chain formation under magnetic field in the curing stage (a) before curing and (b) after curing under magnetic field.
(a) (b)Analysis of the force of the magnetic particles in the MRE during the curing stage mainly includes the magnetic force by the applied magnetic field, the viscous force by the silicone rubber matrix, buoyancy, gravity, the interaction force between the particles, and the Brownian force [23]. The magnetic particles produce motion under the combined action of these forces, so the following equation is obtained according to Newton's formula:(1)midv→idt=F→mi+F→vi+F→ri+F→gi+F→bi+F→Biwheremi is the mass of magnetic particle i,v→i denotes the particle velocity, F→mi is the applied magnetic force, F→vi is the viscous force, F→ri is the interaction force between particles, and F→gi,F→bi, and F→Bi represent the gravity, buoyancy, and Brownian forces, respectively.The magnetic particles used in this study are micron-sized particles; thus, gravity and buoyancy can be considered to be balanced; thus, they are not considered. Further, the magnetic and viscous forces of micron-sized particles in incompressible Newtonian liquids are 200 times that of other forces[18]. Therefore, to facilitate calculation for this, the last four terms in (1) are omitted, and only the magnetic and viscous forces of particles moving in the silicone rubber matrix are considered, and the equation of movement of particles can be simplified to (2)midv→idt=F→mi+F→viIn this case, the ordinary differential equation of displacement can be expressed as follows:(3)u→¨=F→mi-F→vimiThe ordinary differential equation of velocity can also be expressed as follows:(4)v→i˙=F→mi-F→vimi(5)u→i˙=v→iwhere u→i and v→i represent the velocity and displacement of particle movement, respectively.


## 2.2. Magnetic Force

The magnetic force on the magnetic particles can be divided into two aspects: one is the force produced by the magnetization of particles under the applied magnetic field and the other is the force between the magnetic dipoles. Assuming that the magnetic particles in the magnetic field are spherical and uniform in size, the magnetic moment of the magnetic particle can be expressed as follows:(6)m→=VM→=43πR3χH→whereV is the volume of a single particle,V=(4pi∗R∧3)/3, R is the radius, M→ represents magnetization, M→=χH→, χ is the magnetic susceptibility, and H→ is the applied magnetic field strength.The magnetic force between the magnetic particles i and j is as follows:(7)F→mij=m→·∇B→(8)=∑j≠iμ0m→im→j4πrij3d→i·d→j-3d→i·r^ijd→j·r^ijThe magnetic force acting on particle i in the applied magnetic field is as follows:(9)F→mi=μ0m→iH→d→i·h^Based on the above equation, it is concluded that the comprehensive magnetic force of particle i is as follows:(10)F→mi=F→mij+F→mi(11)F→mi=3μ0m→im→j4πrij4∑j≠id→i·d→jr^ij+d→i·r^ijd→j+d→j·r^ijd→i-5d→i·r^ijd→j·r^ijr^ij+μ0m→iH→d→i·h^When only two magnetic particles act, the motion of the particles is a two-dimensional plane motion, so (11) is rewritten:(12)F→mi=4πμ0χ2H2R63rij41-5cos2⁡θijr^ij+2cos⁡θijh^Then the magnetic field force F→mi is decomposed in the x and y directions and obtained its matrix expression:(13)F→mxiF→myi=4πμ0χ2H2R63rij41-5cos2⁡θijsin⁡θij3-5cos2⁡θijcos⁡θijSince the direction of the magnetic force of the particle is related to the angle in the equation (13), so let 1-5cos2⁡θijsin⁡θij=0 and get θij=π/2 or θij=arccos⁡(15/5). Therefore, if the particles are to form a particle chain parallel to the direction of the magnetic field, the angle between the two particles θij<arccos⁡5/5 is required to be satisfied, at this time, and particles appear as close to each. If θij>arccos⁡5/5, the magnetic particles repel each other.As shown in Figure2, μ0 is the vacuum permeability; m→i and m→j are the particle magnetic moments; d→i and d→j are the direction vectors of the magnetic moment of particles i and j; r^ij is the direction vectors from particle i to particle j; H→ is the direction vector of the applied magnetic field intensity; h^ is the direction vector of the applied magnetic field intensity; rij denotes the distance between the centers of two particles i and j and satisfies the following requirement:(14)rij=r^ijFigure 2
Force diagram of particles i and j in a 2D plane under applied magnetic field.


## 2.3. Viscous Force

In a static incompressible fluid, the magnetic particles are subjected to the viscous forceF→v of the matrix when they move in the polymer matrix. The particle size is less than 3 mm; thus, according to the Stokes [24], (15)F→v=-6πRηdrdtR is the radius of the particle and η is the dynamic viscosity. The selected silicone rubber is purchased from Dow Corning, USA, and the viscosity is 3.5 Pa·s at room temperature. dr/dt is the velocity of the particle motion.


## 3. Numerical Simulation

Based on the force analysis of the magnetic particles in the second section, it is shown that the magnetic particles move under magnetic force and viscous force. In this section, we bring the derived motion equation into the finite element software to simulate the two magnetic particles under the magnetic field, which links the magnetic particle motion theory with the numerical simulation.First, a two-dimensional plane is selected to establish an empty model, and then the physical field interface is added, including magnetic field no current, moving mesh, and global ordinary differential and differential algebraic equations. Finally, transient research is added to three physics. The magnetic particles under the action of the field are simulated. The magnetic field no current module generates the applied magnetic field required for the movement of the magnetic particles. The moving mesh realizes the movement of the magnetic particles at different time steps. The global ordinary differential and differential algebra equations provide calculations for the velocity and displacement of the particles. The specific flow chart of the simulation is shown in Figure3.Figure 3
Simulation flow chart.


### 3.1. Establishment of Permanent Magnetic Field

The applied magnetic field adopts a uniform magnetic field to ensure the accuracy of the calculation result. As shown in Figure4, the main magnetic field is generated by a pair of horizontally placed permanent magnets, and the magnetic particle moving region is located in the middle of the permanent magnet. Aiming to form a uniform magnetic field having the same magnetic induction intensity in the horizontal and vertical directions, we used a yoke at the end of the magnet, as shown in the red part of Figure 4. The magnetic field design parameters are shown in Table 1.Table 1
Simulation parameters of uniform magnetic field.
Parameter Value Unit Movement field (d1∗d2) 40∗40 um Magnetic size (L1∗L2) 80∗50 um Width of magnetic yoke (s) 10 umFigure 4
Magnetic field generating device (the uniform magnetic field is generated by a pair of rectangular permanent magnets and yokes, and the middle square area is the moving area of the magnetic particles).In order to verify the uniformity of the two-dimensional plane magnetic field, the cut line 1 and the cut line 2 were, respectively, plotted in the horizontal and vertical midpoint positions of the moving region of the particle, and the magnetic flux density on the cut line was calculated. As shown in Figure5, it can be found that, after the application of the yoke, the variation of the magnetic induction intensity in the moving direction of the particles in the horizontal direction and the vertical direction is remarkably reduced with respect to the nonapplied yoke and thus can be regarded as a uniform magnetic field.Figure 5
The magnitude of the induced intensity on the horizontal and vertical lines (the solid line indicates the distribution of the magnetic field without the yoke, and the dotted line indicates the distribution of the magnetic field with the yoke).


### 3.2. Simulation Results and Discussion

Since the movement of magnetic particles in MRE directly affects its mechanical properties, the single factor analysis method was used to analyze the influence of four different factors on motion of particles in chain formation, including the distribution angle between two magnetic particles, particle radius, the strength of applied magnetic field, and distance between the particles. In the progress of simulation, the relative magnetic permeability of the particles is 3000, and the permanent magnetic field described in Section3.1 is used as the applied magnetic field of the particles.


#### 3.2.1. Influence of Distribution Angles on the Movement of Magnetic Particles

The different types of distribution angles of magnetic particles in the matrix result in different movements of the particles under the magnetic field. In this section, the relationship between the different distribution angles (θ) of two magnetic particles and their motion parameters was discussed. The angle of distribution was selected to be 0°, 30°, 60°, and 90°, the applied magnetic field strength is 50 mT, the radius of the magnetic particles is 10 μm, and the initial distribution distance of the two particles is 60 μm. Figure 6 shows the movement of the particles at different times, the black circle in the figure represents the initial distribution of the particles, and the red arrow represents the direction of movement of the particles. From the motion state of the particles, it can be found that the two particles have two effects of attraction and repulsion. Whenθ =0°, 30°, and 60°, the particles become mutually attracted, and whenθ =90°, the particles repel each other; this is because when the angle between the particles is 90°; it is known from (13) that θ=90°>arccos⁡5/5, so the particles are mutually exclusive.Figure 6
Variation of particle movement position with time in different particle distribution angles ((a), (b), (c), and (d) corresponding to the distribution angle between particlesθ = 0°, 30°, 60°, and 90°).The three distribution angles of the two particles attracted to each other in Figure6 were selected, and the changes of velocity and displacement during the movement of the particles were calculated, as shown in (a), (b), and (c) of Figure 7. As shown in Figure 7, the contact time to attract each other is different, where in the corresponding particle contact time isθ= 0°, 30°, and 60° and t = 0.0032s, t= 0.0042s, and t= 0.0136s, and asθ increases, the contact time is longer. As shown in Figures 7(a), 7(b), and 7(c) enlarged view, comparing the velocities of the particles at t= 0.002s, the average velocity v=(v1+v2)/2 corresponding toθ= 0°, 30°, and 60° gradually decreases, indicating that the smallerθ is, the larger moving speed of the particles is in the same time, thereby the contact time of the two particles is shorten. At the same time, Figures 7(a), 7(b), and 7(c) show the displacement change of the different distribution angles; asθ increases, u1_x and u2_x also gradually increases. This is because the particles first move in the x direction to reduce the angleθ and then attract each other until the particles contact. Also, it can be seen from the displacements in the y direction that the u1_y and u2_y increase during the movement, indicating that the particles have a vertical movement tendency from the beginning of the movement, and this trend was more obvious when the particles are in contact. Finally, due to the particularity of the grid calculation, when the particle motion stops, the calculation still uses the contact velocity as the final velocity. In fact, the velocity of the particle is 0 at this time, that is, when the velocity of the particle tends to be constant, and it can be regarded as the particles come into contact and the movement stops.Figure 7
The velocity and displacement of two particles with time at different distribution angles (the two particles corresponding to (a), (b), and (c) are 0°, 30°, and 60°, respectively).
(a) (b) (c)


#### 3.2.2. Effect of Different Particle Radius on the Movement of Magnetic Particles

MRE usually use micron-sized magnetic particles, so the radius of particles also is a significant effect on the magnetorheological effect of MRE. To analyze the movement of magnetic particles with different radii under magnetic field, 5um, 10um, 15um, and 20um of magnetic particles radius were simulated, and the changes of velocity and displacement parameters during the motion were obtained. In order to consider the influence of the magnetic field in the horizontal direction and the vertical direction, the angleθ between the particle connection and the applied magnetic field strength was 45°. The distance between the two particles was 60 um.As shown in Figure8, the contact time and displacement of the particles become smaller with the larger particles radius; it is because the distance between the particles is constant, with the particle radius increasing, and the corresponding particle volume and mass increase relatively, and then the magnetic force of the particle increases, so that the contact takes less time, which demonstrated that controlling the size of the magnetic particles can effectively control the particle contact time.Figure 8
Variations of velocity and displacement during the mutual attraction of four different magnetic particle radii.


#### 3.2.3. Effect of Different Strength of the Applied Magnetic Field on Particle Motion

The greater the applied magnetic field strength, the greater the magnetic force of the particles. In order to reflect the influence of different applied magnetic field strength on the movement of magnetic particles, four kinds of uniform magnetic fields were selected in this study. The magnetic field strengths were 50mT, 100mT, 150mT, and 200mT, respectively, the corresponding magnetic particles have a radius of 5 um, the distance between the particles is 30 um, and the angle between the particles is 45°.It can be seen from Figure9 that the velocity and displacement of the two particles change significantly with the strength of the applied magnetic field increases, the velocity of the particles increases slowly especially at the beginning of the motion, but the particles rapidly increase upon contact; at the same time, the displacement of the two particles in the x and y directions gradually increases, indicating that the particles move in parallel and perpendicular to the direction of the magnetic field. When the particles are in contacting, the displacement in the y direction increases rapidly, and the contact of the particles forms a particle chain parallel to the direction of the applied magnetic field. From the contact time of the two particles, the greater the magnetic field strength, the shorter the contact time of the particles. The particle contact time is 0.003s in 50 mT, while the time is 0.0002s in 200mT. That is, the time for the two particles to be chained is reduced by 93%. Because the stronger the magnetic field strength, the stronger the magnetic force of the particles under the same magnetic permeability. It shows that increasing the strength of the applied magnetic field can significantly improve the efficiency of particle contact into the chain.Figure 9
Variations in the velocity and displacement of magnetic particles with different applied magnetic field strengths.


#### 3.2.4. Effect of Different Distance between Particles on Particle Motion

Due to the random distribution of magnetic particles, there are different distribution distances between the two particles. In this section, four different particle distribution distances were selected, which are 4∗R, 6∗R, 8∗R, and 10∗R. The magnetic field is 50mT, the radius of the magnetic particle R is 5um, and the angle between the particle connection and the magnetic field is 45°. As shown in Figure 10, the longer the distance between the particles, the longer the contact time of the particles. It can be seen from the displacement curve of the particles that the larger the distance between the particles that the larger the displacement in the x and y directions, the longer the contact time.Figure 10
Variation of velocity and displacement of particles with time at different particle distribution distances.


#### 3.2.5. Effect of Different Conditions on Contact Time of Magnetic Particles

Based on the analysis of the previous sections, the contact time t of the two particles at different angles (θ), particle radius (R), applied magnetic field strength (B), and the distance (L) between the particles was considered, and the influence of the factor on the motion state of the particles was judged. Select initial conditions θ0,R0,B0,L0 for 30°, 10um, 100mT, and 50um and the corresponding scale factor λ=0.5, 1.0, 1.5, and 2.0 under a single factor. The contact time of each factor under the four scale factors was discussed. The corresponding proportional conditions are shown in Table 2 and the calculation results are shown in Table 3.Table 2
Proportional condition setting.
Factor λ = 0.5 λ = 1.0 λ = 1.5 λ = 2.0 θ (°) 0.5 ∗ θ 0 , R 0 , B 0 , L 0 θ 0 , R 0 , B 0 , L 0 1.5 ∗ θ 0 , R 0 , B 0 , L 0 2 ∗ θ 0 , R 0 , B 0 , L 0 R(um) θ 0 , 0.5 ∗ R 0 , B 0 , L 0 θ 0 , R 0 , B 0 , L 0 θ 0 , 1.5 ∗ R 0 , B 0 , L 0 θ 0 , 2 ∗ R 0 , B 0 , L 0 B(mT) θ 0 , R 0 , 0.5 ∗ B 0 , L 0 θ 0 , R 0 , B 0 , L 0 θ 0 , R 0 , 1.5 B 0 , L 0 θ 0 , R 0 , 2 ∗ B 0 , L 0 L(um) θ 0 , R 0 , B 0 , 0.5 ∗ L 0 θ 0 , R 0 , B 0 , L 0 θ 0 , R 0 , B 0 , 1.5 ∗ L 0 θ 0 , R 0 , B 0 , 2 ∗ L 0Table 3
Calculated contact time (t).
Factor λ = 0.5 λ = 1.0 λ = 1.5 λ = 2.0 θ (°) 4.32E-4 5.15E-4 7.60E-4 1.52E-3 R(um) 2.90E-4 5.15E-4 1.03E-4 2.65E-5 B(mT) 2.00E-3 5.15E-4 2.38E-4 1.40E-4 L(um) 1.59E-5 5.15E-4 2.38E-3 6.42E-3According to Table3, the change of particle contact time under different conditions was obtained, as shown in Figure 11. The height of the histogram reflects the contact time t of the two particles, and the red arrow reflects the change of the contact time of the particles under different ratios. As can be seen from Figure 11, changing the ratio of the four conditions has a significant effect on the contact time of the two particles, comparing the rate of change of particle contact time and the change rate of contact time corresponding to the distanceL between the two particles is the largest, reaching 1143%. Secondly, the particle radiusR, the magnetic field strengthB and the distribution angleθ, and the corresponding contact time change rate are 345%, 288%, and 191%, which means that if the particle chain formation efficiency is improved, the particle contact time is shortened; this is achieved by increasing the concentration of the particles to reduce the distance between the particles.Figure 11
Change in contact time of two magnetic particles under different conditions.For the magnetic particles in the silicone rubber-based MRE, this section discusses the effects of different distribution angles, particle radii, applied magnetic field strength, and the distance between the particles on the velocity and displacement during mutual attraction, and the results show that these three different conditions have different effects on the movement of the particles. Secondly, the contact time of the particles under different ratios was calculated, and the distance between the particles has the greatest influence on the contact time, and the influence of the particle distribution angle on the contact time of the particles is relatively small.


## 3.1. Establishment of Permanent Magnetic Field

The applied magnetic field adopts a uniform magnetic field to ensure the accuracy of the calculation result. As shown in Figure4, the main magnetic field is generated by a pair of horizontally placed permanent magnets, and the magnetic particle moving region is located in the middle of the permanent magnet. Aiming to form a uniform magnetic field having the same magnetic induction intensity in the horizontal and vertical directions, we used a yoke at the end of the magnet, as shown in the red part of Figure 4. The magnetic field design parameters are shown in Table 1.Table 1
Simulation parameters of uniform magnetic field.
Parameter Value Unit Movement field (d1∗d2) 40∗40 um Magnetic size (L1∗L2) 80∗50 um Width of magnetic yoke (s) 10 umFigure 4
Magnetic field generating device (the uniform magnetic field is generated by a pair of rectangular permanent magnets and yokes, and the middle square area is the moving area of the magnetic particles).In order to verify the uniformity of the two-dimensional plane magnetic field, the cut line 1 and the cut line 2 were, respectively, plotted in the horizontal and vertical midpoint positions of the moving region of the particle, and the magnetic flux density on the cut line was calculated. As shown in Figure5, it can be found that, after the application of the yoke, the variation of the magnetic induction intensity in the moving direction of the particles in the horizontal direction and the vertical direction is remarkably reduced with respect to the nonapplied yoke and thus can be regarded as a uniform magnetic field.Figure 5
The magnitude of the induced intensity on the horizontal and vertical lines (the solid line indicates the distribution of the magnetic field without the yoke, and the dotted line indicates the distribution of the magnetic field with the yoke).


## 3.2. Simulation Results and Discussion

Since the movement of magnetic particles in MRE directly affects its mechanical properties, the single factor analysis method was used to analyze the influence of four different factors on motion of particles in chain formation, including the distribution angle between two magnetic particles, particle radius, the strength of applied magnetic field, and distance between the particles. In the progress of simulation, the relative magnetic permeability of the particles is 3000, and the permanent magnetic field described in Section3.1 is used as the applied magnetic field of the particles.


### 3.2.1. Influence of Distribution Angles on the Movement of Magnetic Particles

The different types of distribution angles of magnetic particles in the matrix result in different movements of the particles under the magnetic field. In this section, the relationship between the different distribution angles (θ) of two magnetic particles and their motion parameters was discussed. The angle of distribution was selected to be 0°, 30°, 60°, and 90°, the applied magnetic field strength is 50 mT, the radius of the magnetic particles is 10 μm, and the initial distribution distance of the two particles is 60 μm. Figure 6 shows the movement of the particles at different times, the black circle in the figure represents the initial distribution of the particles, and the red arrow represents the direction of movement of the particles. From the motion state of the particles, it can be found that the two particles have two effects of attraction and repulsion. Whenθ =0°, 30°, and 60°, the particles become mutually attracted, and whenθ =90°, the particles repel each other; this is because when the angle between the particles is 90°; it is known from (13) that θ=90°>arccos⁡5/5, so the particles are mutually exclusive.Figure 6
Variation of particle movement position with time in different particle distribution angles ((a), (b), (c), and (d) corresponding to the distribution angle between particlesθ = 0°, 30°, 60°, and 90°).The three distribution angles of the two particles attracted to each other in Figure6 were selected, and the changes of velocity and displacement during the movement of the particles were calculated, as shown in (a), (b), and (c) of Figure 7. As shown in Figure 7, the contact time to attract each other is different, where in the corresponding particle contact time isθ= 0°, 30°, and 60° and t = 0.0032s, t= 0.0042s, and t= 0.0136s, and asθ increases, the contact time is longer. As shown in Figures 7(a), 7(b), and 7(c) enlarged view, comparing the velocities of the particles at t= 0.002s, the average velocity v=(v1+v2)/2 corresponding toθ= 0°, 30°, and 60° gradually decreases, indicating that the smallerθ is, the larger moving speed of the particles is in the same time, thereby the contact time of the two particles is shorten. At the same time, Figures 7(a), 7(b), and 7(c) show the displacement change of the different distribution angles; asθ increases, u1_x and u2_x also gradually increases. This is because the particles first move in the x direction to reduce the angleθ and then attract each other until the particles contact. Also, it can be seen from the displacements in the y direction that the u1_y and u2_y increase during the movement, indicating that the particles have a vertical movement tendency from the beginning of the movement, and this trend was more obvious when the particles are in contact. Finally, due to the particularity of the grid calculation, when the particle motion stops, the calculation still uses the contact velocity as the final velocity. In fact, the velocity of the particle is 0 at this time, that is, when the velocity of the particle tends to be constant, and it can be regarded as the particles come into contact and the movement stops.Figure 7
The velocity and displacement of two particles with time at different distribution angles (the two particles corresponding to (a), (b), and (c) are 0°, 30°, and 60°, respectively).
(a) (b) (c)


### 3.2.2. Effect of Different Particle Radius on the Movement of Magnetic Particles

MRE usually use micron-sized magnetic particles, so the radius of particles also is a significant effect on the magnetorheological effect of MRE. To analyze the movement of magnetic particles with different radii under magnetic field, 5um, 10um, 15um, and 20um of magnetic particles radius were simulated, and the changes of velocity and displacement parameters during the motion were obtained. In order to consider the influence of the magnetic field in the horizontal direction and the vertical direction, the angleθ between the particle connection and the applied magnetic field strength was 45°. The distance between the two particles was 60 um.As shown in Figure8, the contact time and displacement of the particles become smaller with the larger particles radius; it is because the distance between the particles is constant, with the particle radius increasing, and the corresponding particle volume and mass increase relatively, and then the magnetic force of the particle increases, so that the contact takes less time, which demonstrated that controlling the size of the magnetic particles can effectively control the particle contact time.Figure 8
Variations of velocity and displacement during the mutual attraction of four different magnetic particle radii.


### 3.2.3. Effect of Different Strength of the Applied Magnetic Field on Particle Motion

The greater the applied magnetic field strength, the greater the magnetic force of the particles. In order to reflect the influence of different applied magnetic field strength on the movement of magnetic particles, four kinds of uniform magnetic fields were selected in this study. The magnetic field strengths were 50mT, 100mT, 150mT, and 200mT, respectively, the corresponding magnetic particles have a radius of 5 um, the distance between the particles is 30 um, and the angle between the particles is 45°.It can be seen from Figure9 that the velocity and displacement of the two particles change significantly with the strength of the applied magnetic field increases, the velocity of the particles increases slowly especially at the beginning of the motion, but the particles rapidly increase upon contact; at the same time, the displacement of the two particles in the x and y directions gradually increases, indicating that the particles move in parallel and perpendicular to the direction of the magnetic field. When the particles are in contacting, the displacement in the y direction increases rapidly, and the contact of the particles forms a particle chain parallel to the direction of the applied magnetic field. From the contact time of the two particles, the greater the magnetic field strength, the shorter the contact time of the particles. The particle contact time is 0.003s in 50 mT, while the time is 0.0002s in 200mT. That is, the time for the two particles to be chained is reduced by 93%. Because the stronger the magnetic field strength, the stronger the magnetic force of the particles under the same magnetic permeability. It shows that increasing the strength of the applied magnetic field can significantly improve the efficiency of particle contact into the chain.Figure 9
Variations in the velocity and displacement of magnetic particles with different applied magnetic field strengths.


### 3.2.4. Effect of Different Distance between Particles on Particle Motion

Due to the random distribution of magnetic particles, there are different distribution distances between the two particles. In this section, four different particle distribution distances were selected, which are 4∗R, 6∗R, 8∗R, and 10∗R. The magnetic field is 50mT, the radius of the magnetic particle R is 5um, and the angle between the particle connection and the magnetic field is 45°. As shown in Figure 10, the longer the distance between the particles, the longer the contact time of the particles. It can be seen from the displacement curve of the particles that the larger the distance between the particles that the larger the displacement in the x and y directions, the longer the contact time.Figure 10
Variation of velocity and displacement of particles with time at different particle distribution distances.


### 3.2.5. Effect of Different Conditions on Contact Time of Magnetic Particles

Based on the analysis of the previous sections, the contact time t of the two particles at different angles (θ), particle radius (R), applied magnetic field strength (B), and the distance (L) between the particles was considered, and the influence of the factor on the motion state of the particles was judged. Select initial conditions θ0,R0,B0,L0 for 30°, 10um, 100mT, and 50um and the corresponding scale factor λ=0.5, 1.0, 1.5, and 2.0 under a single factor. The contact time of each factor under the four scale factors was discussed. The corresponding proportional conditions are shown in Table 2 and the calculation results are shown in Table 3.Table 2
Proportional condition setting.
Factor λ = 0.5 λ = 1.0 λ = 1.5 λ = 2.0 θ (°) 0.5 ∗ θ 0 , R 0 , B 0 , L 0 θ 0 , R 0 , B 0 , L 0 1.5 ∗ θ 0 , R 0 , B 0 , L 0 2 ∗ θ 0 , R 0 , B 0 , L 0 R(um) θ 0 , 0.5 ∗ R 0 , B 0 , L 0 θ 0 , R 0 , B 0 , L 0 θ 0 , 1.5 ∗ R 0 , B 0 , L 0 θ 0 , 2 ∗ R 0 , B 0 , L 0 B(mT) θ 0 , R 0 , 0.5 ∗ B 0 , L 0 θ 0 , R 0 , B 0 , L 0 θ 0 , R 0 , 1.5 B 0 , L 0 θ 0 , R 0 , 2 ∗ B 0 , L 0 L(um) θ 0 , R 0 , B 0 , 0.5 ∗ L 0 θ 0 , R 0 , B 0 , L 0 θ 0 , R 0 , B 0 , 1.5 ∗ L 0 θ 0 , R 0 , B 0 , 2 ∗ L 0Table 3
Calculated contact time (t).
Factor λ = 0.5 λ = 1.0 λ = 1.5 λ = 2.0 θ (°) 4.32E-4 5.15E-4 7.60E-4 1.52E-3 R(um) 2.90E-4 5.15E-4 1.03E-4 2.65E-5 B(mT) 2.00E-3 5.15E-4 2.38E-4 1.40E-4 L(um) 1.59E-5 5.15E-4 2.38E-3 6.42E-3According to Table3, the change of particle contact time under different conditions was obtained, as shown in Figure 11. The height of the histogram reflects the contact time t of the two particles, and the red arrow reflects the change of the contact time of the particles under different ratios. As can be seen from Figure 11, changing the ratio of the four conditions has a significant effect on the contact time of the two particles, comparing the rate of change of particle contact time and the change rate of contact time corresponding to the distanceL between the two particles is the largest, reaching 1143%. Secondly, the particle radiusR, the magnetic field strengthB and the distribution angleθ, and the corresponding contact time change rate are 345%, 288%, and 191%, which means that if the particle chain formation efficiency is improved, the particle contact time is shortened; this is achieved by increasing the concentration of the particles to reduce the distance between the particles.Figure 11
Change in contact time of two magnetic particles under different conditions.For the magnetic particles in the silicone rubber-based MRE, this section discusses the effects of different distribution angles, particle radii, applied magnetic field strength, and the distance between the particles on the velocity and displacement during mutual attraction, and the results show that these three different conditions have different effects on the movement of the particles. Secondly, the contact time of the particles under different ratios was calculated, and the distance between the particles has the greatest influence on the contact time, and the influence of the particle distribution angle on the contact time of the particles is relatively small.


## 3.2.1. Influence of Distribution Angles on the Movement of Magnetic Particles

The different types of distribution angles of magnetic particles in the matrix result in different movements of the particles under the magnetic field. In this section, the relationship between the different distribution angles (θ) of two magnetic particles and their motion parameters was discussed. The angle of distribution was selected to be 0°, 30°, 60°, and 90°, the applied magnetic field strength is 50 mT, the radius of the magnetic particles is 10 μm, and the initial distribution distance of the two particles is 60 μm. Figure 6 shows the movement of the particles at different times, the black circle in the figure represents the initial distribution of the particles, and the red arrow represents the direction of movement of the particles. From the motion state of the particles, it can be found that the two particles have two effects of attraction and repulsion. Whenθ =0°, 30°, and 60°, the particles become mutually attracted, and whenθ =90°, the particles repel each other; this is because when the angle between the particles is 90°; it is known from (13) that θ=90°>arccos⁡5/5, so the particles are mutually exclusive.Figure 6
Variation of particle movement position with time in different particle distribution angles ((a), (b), (c), and (d) corresponding to the distribution angle between particlesθ = 0°, 30°, 60°, and 90°).The three distribution angles of the two particles attracted to each other in Figure6 were selected, and the changes of velocity and displacement during the movement of the particles were calculated, as shown in (a), (b), and (c) of Figure 7. As shown in Figure 7, the contact time to attract each other is different, where in the corresponding particle contact time isθ= 0°, 30°, and 60° and t = 0.0032s, t= 0.0042s, and t= 0.0136s, and asθ increases, the contact time is longer. As shown in Figures 7(a), 7(b), and 7(c) enlarged view, comparing the velocities of the particles at t= 0.002s, the average velocity v=(v1+v2)/2 corresponding toθ= 0°, 30°, and 60° gradually decreases, indicating that the smallerθ is, the larger moving speed of the particles is in the same time, thereby the contact time of the two particles is shorten. At the same time, Figures 7(a), 7(b), and 7(c) show the displacement change of the different distribution angles; asθ increases, u1_x and u2_x also gradually increases. This is because the particles first move in the x direction to reduce the angleθ and then attract each other until the particles contact. Also, it can be seen from the displacements in the y direction that the u1_y and u2_y increase during the movement, indicating that the particles have a vertical movement tendency from the beginning of the movement, and this trend was more obvious when the particles are in contact. Finally, due to the particularity of the grid calculation, when the particle motion stops, the calculation still uses the contact velocity as the final velocity. In fact, the velocity of the particle is 0 at this time, that is, when the velocity of the particle tends to be constant, and it can be regarded as the particles come into contact and the movement stops.Figure 7
The velocity and displacement of two particles with time at different distribution angles (the two particles corresponding to (a), (b), and (c) are 0°, 30°, and 60°, respectively).
(a) (b) (c)


## 3.2.2. Effect of Different Particle Radius on the Movement of Magnetic Particles

MRE usually use micron-sized magnetic particles, so the radius of particles also is a significant effect on the magnetorheological effect of MRE. To analyze the movement of magnetic particles with different radii under magnetic field, 5um, 10um, 15um, and 20um of magnetic particles radius were simulated, and the changes of velocity and displacement parameters during the motion were obtained. In order to consider the influence of the magnetic field in the horizontal direction and the vertical direction, the angleθ between the particle connection and the applied magnetic field strength was 45°. The distance between the two particles was 60 um.As shown in Figure8, the contact time and displacement of the particles become smaller with the larger particles radius; it is because the distance between the particles is constant, with the particle radius increasing, and the corresponding particle volume and mass increase relatively, and then the magnetic force of the particle increases, so that the contact takes less time, which demonstrated that controlling the size of the magnetic particles can effectively control the particle contact time.Figure 8
Variations of velocity and displacement during the mutual attraction of four different magnetic particle radii.


## 3.2.3. Effect of Different Strength of the Applied Magnetic Field on Particle Motion

The greater the applied magnetic field strength, the greater the magnetic force of the particles. In order to reflect the influence of different applied magnetic field strength on the movement of magnetic particles, four kinds of uniform magnetic fields were selected in this study. The magnetic field strengths were 50mT, 100mT, 150mT, and 200mT, respectively, the corresponding magnetic particles have a radius of 5 um, the distance between the particles is 30 um, and the angle between the particles is 45°.It can be seen from Figure9 that the velocity and displacement of the two particles change significantly with the strength of the applied magnetic field increases, the velocity of the particles increases slowly especially at the beginning of the motion, but the particles rapidly increase upon contact; at the same time, the displacement of the two particles in the x and y directions gradually increases, indicating that the particles move in parallel and perpendicular to the direction of the magnetic field. When the particles are in contacting, the displacement in the y direction increases rapidly, and the contact of the particles forms a particle chain parallel to the direction of the applied magnetic field. From the contact time of the two particles, the greater the magnetic field strength, the shorter the contact time of the particles. The particle contact time is 0.003s in 50 mT, while the time is 0.0002s in 200mT. That is, the time for the two particles to be chained is reduced by 93%. Because the stronger the magnetic field strength, the stronger the magnetic force of the particles under the same magnetic permeability. It shows that increasing the strength of the applied magnetic field can significantly improve the efficiency of particle contact into the chain.Figure 9
Variations in the velocity and displacement of magnetic particles with different applied magnetic field strengths.


## 3.2.4. Effect of Different Distance between Particles on Particle Motion

Due to the random distribution of magnetic particles, there are different distribution distances between the two particles. In this section, four different particle distribution distances were selected, which are 4∗R, 6∗R, 8∗R, and 10∗R. The magnetic field is 50mT, the radius of the magnetic particle R is 5um, and the angle between the particle connection and the magnetic field is 45°. As shown in Figure 10, the longer the distance between the particles, the longer the contact time of the particles. It can be seen from the displacement curve of the particles that the larger the distance between the particles that the larger the displacement in the x and y directions, the longer the contact time.Figure 10
Variation of velocity and displacement of particles with time at different particle distribution distances.


## 3.2.5. Effect of Different Conditions on Contact Time of Magnetic Particles

Based on the analysis of the previous sections, the contact time t of the two particles at different angles (θ), particle radius (R), applied magnetic field strength (B), and the distance (L) between the particles was considered, and the influence of the factor on the motion state of the particles was judged. Select initial conditions θ0,R0,B0,L0 for 30°, 10um, 100mT, and 50um and the corresponding scale factor λ=0.5, 1.0, 1.5, and 2.0 under a single factor. The contact time of each factor under the four scale factors was discussed. The corresponding proportional conditions are shown in Table 2 and the calculation results are shown in Table 3.Table 2
Proportional condition setting.
Factor λ = 0.5 λ = 1.0 λ = 1.5 λ = 2.0 θ (°) 0.5 ∗ θ 0 , R 0 , B 0 , L 0 θ 0 , R 0 , B 0 , L 0 1.5 ∗ θ 0 , R 0 , B 0 , L 0 2 ∗ θ 0 , R 0 , B 0 , L 0 R(um) θ 0 , 0.5 ∗ R 0 , B 0 , L 0 θ 0 , R 0 , B 0 , L 0 θ 0 , 1.5 ∗ R 0 , B 0 , L 0 θ 0 , 2 ∗ R 0 , B 0 , L 0 B(mT) θ 0 , R 0 , 0.5 ∗ B 0 , L 0 θ 0 , R 0 , B 0 , L 0 θ 0 , R 0 , 1.5 B 0 , L 0 θ 0 , R 0 , 2 ∗ B 0 , L 0 L(um) θ 0 , R 0 , B 0 , 0.5 ∗ L 0 θ 0 , R 0 , B 0 , L 0 θ 0 , R 0 , B 0 , 1.5 ∗ L 0 θ 0 , R 0 , B 0 , 2 ∗ L 0Table 3
Calculated contact time (t).
Factor λ = 0.5 λ = 1.0 λ = 1.5 λ = 2.0 θ (°) 4.32E-4 5.15E-4 7.60E-4 1.52E-3 R(um) 2.90E-4 5.15E-4 1.03E-4 2.65E-5 B(mT) 2.00E-3 5.15E-4 2.38E-4 1.40E-4 L(um) 1.59E-5 5.15E-4 2.38E-3 6.42E-3According to Table3, the change of particle contact time under different conditions was obtained, as shown in Figure 11. The height of the histogram reflects the contact time t of the two particles, and the red arrow reflects the change of the contact time of the particles under different ratios. As can be seen from Figure 11, changing the ratio of the four conditions has a significant effect on the contact time of the two particles, comparing the rate of change of particle contact time and the change rate of contact time corresponding to the distanceL between the two particles is the largest, reaching 1143%. Secondly, the particle radiusR, the magnetic field strengthB and the distribution angleθ, and the corresponding contact time change rate are 345%, 288%, and 191%, which means that if the particle chain formation efficiency is improved, the particle contact time is shortened; this is achieved by increasing the concentration of the particles to reduce the distance between the particles.Figure 11
Change in contact time of two magnetic particles under different conditions.For the magnetic particles in the silicone rubber-based MRE, this section discusses the effects of different distribution angles, particle radii, applied magnetic field strength, and the distance between the particles on the velocity and displacement during mutual attraction, and the results show that these three different conditions have different effects on the movement of the particles. Secondly, the contact time of the particles under different ratios was calculated, and the distance between the particles has the greatest influence on the contact time, and the influence of the particle distribution angle on the contact time of the particles is relatively small.


## 4. Conclusions

In this paper, through the theoretical analysis and simulation calculation of magnetic particles in MRE, the following conclusions are obtained:(1)
The force analysis of the magnetic particles in the silicone rubber-based MRE was carried out, and the motion model of the magnetic particles under the applied magnetic field was established.(2)
The motion of magnetic particles in the MRE were simulated under the uniform magnetic field by the finite element simulation method, and the effects of different particle distribution angles, different particle sizes, different applied magnetic fields, and different particle distances on the velocity and displacement during particle motion were discussed; the simulated result reveal that the distance between particles has the greatest influence on contact time.(3)
The simulation on motion of two magnetic particles provided guiding for optimizing the chaining formation conditions of MRE in the curing process.



---

*Source: 8182651-2019-07-18.xml*