Abstract
This paper studies on the permanent magnets configuration on the transmission torque of the multi-pole bilayer magnetorheological (MR) coupling. Based on the electromagnetic field theory, the magnetic circuit models of traditional permanent magnet array (TPMA) and Halbach permanent magnet array (HPMA) are established, and the magnetic flux density within the MR fluid working gaps has been derived in order to evaluate the merits of the designed MR coupling. A 3D FE magnetic-fluid analysis has been necessary following the initial conceptual analysis, in order to study the influence of key parameters on the transmission torque. The results show that the transmission torque of the MR coupling with Halbach permanent magnet array is 33.45% higher than that of the ordinary permanent magnet array, with a same structure size. For the MR coupling with Halbach permanent magnet array, the unilateral magnetic focusing effect is better with the increase of the residual flux density of the secondary magnetic pole as well as the radial length of magnetic pole. And the single side magnetic focusing effect is the best when the main magnetic pole is 15°. The influence of the magnetic pole angle on the transmission torque has been further studied.
Keywords
Introduction
Magnetorheological (MR) fluid is a new type of smart material. Due to its fast, controllable, and reversible liquid-solid transformation under the applied magnetic field [1,2], it is widely used in brakes [3], clutches [4], dampers [5,6], MR valves [7], MR couplings [8], MR polishing machines [9], and others [10,11]. Compared to traditional mechanical transmission devices, MR transmission devices which utilizing the characteristics of MR fluid has the advantages of quick response, and adjustable torque output. Therefore, it has a mature application in various transmission scenarios.
At present, MR transmission devices mainly include disc type, cylindrical type, multi-pole type, multi-disc type, multi-cylinder type and other structures. Li et al. [12] established the coupling model of structure and magnetic circuit, and optimized the design method after simulation. By optimizing the design of the disc MR brake, the power density of the brake is improved. Yuan et al. [13] designed a disc MR brake with a circular groove, and its braking torque is 35% higher than that of the common disc MR brake. Zhao et al. [14] designed a cylindrical MR transmission controlled by shape memory alloy, and analyzed its braking performance. Wu et al. [15] designed a high-power multi-level and multi-layer MR brake based on the principle of magnetic field superposition. Compared with the traditional single-disc, multi-disc, or multi-stage single-layer MR brakes, it has a larger braking torque and power density. Wu et al. [16] proposed a multi-disc MR transmission device, which uses the electromagnetic force generated by the coil to squeeze the MR fluid, thereby improving the transmission torque. Qin et al. [17] designed a multi-cylinder MR brake, using the magnetic fluid sealing technology to reduce the unsteady friction, and deduced its current-torque model to predict the output torque according to the input current size.
Aiming at the problems of large power consumption, complicated structure, and coil heating of the coil-type MR transmission devices. Bucchi et al. [18] designed a MR clutch with a magnetic field provided by a permanent magnet, and changed the working state of the clutch by changing the position of the permanent magnet, and analyzed its transmission performance. Further, Rizzo et al. [19–21] proposed a new gap shape on this basis to reduce the torque loss in the separation operation, and the eddy current in the working process of the clutch was analyzed. Chen [22] proposed a MR coupling with a magnetic field provided by a permanent magnet, and analyzed the magnetic field distribution characteristics of the region where the MR fluid of the device is located, as well as the circumferential distribution of shear yield stress and viscosity in polar coordinates for two specified locations. Böse et al. [23] proposed a magnetorheological clutch with magnetic field provided by permanent magnet and coil. Based on this, the torque transmission performance of three magnetorheological clutches with different structures was designed, manufactured, and tested. Moghani et al. [24] proposed a hybrid magnetorheological clutch which adjusts the magnetic field generated by the permanent magnet through the excitation coil, and they carried out multi-objective optimization design for the transmission torque, clutch quality, and magnetic field strength. The conclusion is that the advantage of the permanent magnet MR transmission devices is that the coil type MR devices will not be unable to transmit torque due to power failure, but not limited to low power consumption and simple structure.
For the MR devices with permanent magnet arrays, the Halbach permanent magnet array can make the device more compact, and the power density of the device is higher than that of the traditional permanent magnet array. The unique arrangement of the Halbach permanent magnet array allows the magnetic field generated by it to concentrate on one side of the array, which is very suitable for applications where only one side of the magnetic field is needed. It can utilize most of the magnetic field, which improves the utilization rate of the magnetic field [25].
When the Halbach permanent magnet array is used, it avoids the loss and heat generated at the stator due to the eddy current effect, and saves more space than the traditional permanent magnet array.
To date, there are few researches focused on the performance of MR couplings with different permanent magnets configuration. In this paper, the effects of the permanent magnets’ configuration and the permanent magnets parameters on the performance of the multi-pole bilayer MR coupling is studied. Section 2 introduces the structure of the multi-pole bilayer permanent magnets MR coupling. Section 3 derives the transmission torque. Section 4 establishes the magnetic circuit modeling. Section 5 compares and analyzes the two configurations the multi-pole bilayer MR coupling by finite element simulation. Section 6 concludes the study.
Structural design of the multi-pole bilayer permanent magnets MR coupling
Figure 1 shows the multi-pole bilayer permanent magnets MR coupling with traditional permanent magnet array (TPMA) and Halbach permanent magnet array (HPMA). As shown in Fig. 1(a) and (b), the coupling is mainly composed of an inner stator, a MR fluid, a permanent magnet array, and a shell. An outer MR fluid gap is formed between the shell and the permanent magnet array, while an inner MR fluid gap is located between the inner stator and the permanent magnet array. The transmission torque of the coupling is generated by the shear between the MR fluid in the inner and outer gaps.

Structural diagram of the multi-pole bilayer MR coupling with: (a) TPMA; (b) HPMA.
In this design, R 1 is the inner diameter of the inner stator, R 2 is the outer diameter of the inner stator as well as the inner diameter of the inner MR fluid gap, R 3 is the outer diameter of the inner MR fluid gap, R 4 is the inner diameter of the outer MR fluid gap, and R 5 is the outer diameter of the outer MR fluid gap, R 6 is the outer diameter of the shell, z is the axial length of the inner and outer MR fluid gap, and c 1 is the main pole angle.
In terms of material selection, neodymium iron boron (NdFeB) with high magnetic energy product and coercive force is selected as the permanent magnet. MRF-132DG from Lord Company is selected as the MR fluid. The inner stator and the shell are made of electric pure iron (20#) which is magnetically conductive and not easily magnetized. Material properties are in Table 1.
Material parameters
The size of structural parameters is the same, except that the arrangement of permanent magnets is different. It is seen that the arrows indicate the direction of the residual magnetic flux of the permanent magnets. For the multi-pole bilayer MR coupling with TPMA, the permanent magnet array is composed of permanent magnets and aluminum blocks arranged at intervals along the circumferential direction, and the magnetic pole directions of adjacent permanent magnets are opposite. For the multipole bilayer MR coupling with HPMA, the arrangement of permanent magnets is a Halbach permanent magnet array in which the magnetic field direction of the main magnetic pole is along the radial direction while the magnetic field direction of the secondary magnetic pole is along the circumferential direction.
The torque transmission of multi-pole bilayer MR coupling depends on the dynamic yield stress and viscosity of MR fluids. Due to the effectiveness and simplicity, the Bingham model is commonly used in designs of the MR transmission device. According to the Bingham model, the shear stress equation of MR fluids is established as [26,27]:
The dynamic shear stress τ
y
can be obtained from the yield characteristics and magnetization characteristics of the MR fluid. For the MRF-132DG, the shear yield stress is fitted polynomials by the least squares’ method, and its expression is obtained [28]:
Shear rate
The magnetic field strength in MR fluid is evenly distributed, and the edge effect is ignored. The transmission torque generated by the inner MR fluid gap can be expressed as:
Substituting Eqs (1) and (4) into Eq. (3):
For Eq. (5), integrating it on the inner MR fluid gap:
By simplifying it, the expression of torque generated by the inner MR fluid gap is obtained [29]:
Similarly, the torque generated by the outer MR fluid gap can be expressed as:
Therefore, the total torque of the multi-pole bilayer MR coupling is:
It can be seen from Eq. (9) that the magnitude of the transmission torque is proportional to the axial width of the MR fluid gap, and is positively related to the radius of the MR fluid gap. The larger the radius of the MR fluid gap is, the bigger the transmission torque will be. In order to increase the transmission torque of the multi-pole bilayer MR coupling, the magnetic field needs to be concentrated as far as possible at the outer gap. For this reason, the influence of the arrangement of the permanent magnets on the transmission torque is analyzed.
In this section, the magnetic density in the MR gaps is firstly derived based on the magnetic circuit analysis. The magnetic circuit is modeled by the average magnetic circuit method without considering the magnetic leakage. The magnetic circuit model and simplified equivalent magnetic circuit of the multi-pole bilayer MR coupling with TPMA and HPMA are shown in Fig. 2 and Fig. 3, respectively.

Magnetic circuit node configuration and simplified magnetic circuit of TPMA: (a) node configuration of one-sixth of the bilayer MR coupling; (b) simplified magnetic circuit of one-sixth of the bilayer MR coupling.

Magnetic circuit model and simplified equivalent magnetic circuit of HPMA: (a) node configuration of one-sixth of the bilayer MR coupling; (b) simplified magnetic circuit of one-sixth of the bilayer MR coupling.
The components in the magnetic circuit of multi-pole bilayer MR coupling are divided into sections and connected at nodes according to the difference of crosssectional area and material of the magnetic circuit. R ab , R a ′ b ′ , R gh , R g ′ h ′ are the magnetic resistance of the inner gap MR fluid ab, a ′ b ′ , gh, g ′ h ′ , respectively. R bc , R b ′ c ′ , R fg , R f ′ g ′ are the magnetic resistance of the main pole bc, b ′ c ′ , fg, f ′ g ′ , respectively. R cd , R c ′ d ′ , R ef , R e ′ f ′ are the magnetic resistance of the outer MR fluid gap cd, c ′ d ′ , ef, e ′ f ′ , respectively. R de , R d ′ e ′ are the magnetic resistance of the shell de, d ′ e ′ , respectively. R ah , R a ′ h ′ are the magnetic resistance of the inner stator ah, a ′ h ′ , respectively. R i ′ j ′ is the magnetic resistance of the permanent magnets i ′ j ′ .
The magnetomotive force in the permanent magnet circuit can be expressed as:
The expression for the magnetic flux Φ obtained from Ohm’s law of the magnetic circuit:
Substituting the effective magnetic flux area S
i
and magnetoresistance of each magnetic circuit into Eq. (12) to obtain the average magnetic flux density in each magnetic circuit. Thus, the magnetoresistance can be calculated:
The reluctance of each section of the magnetic circuit is:
For the multi-pole bilayer MR coupling with TPMA, the magnetomotive force can be calculated according to the first law of magnetic circuit:
The magnetic flux Φ can be calculated by substituting the magnetomotive force and the reluctance of each section of the magnetic circuit. Then, the magnetic flux density of each section of the magnetic circuit can be obtained by substituting the result into the Eq. (12). The magnetic flux density within the inner and outer gaps can be expressed as:
According to the magnetic circuit law:
The magnetic fluxes Φ1, Φ2 and Φ3 in the magnetic circuit can be calculated from Eqs ((23))–((25)). Substituting the above results into Eq. ((12)), then the magnetic flux density of each section of the magnetic circuit for the multi-pole bilayer MR coupling with HPMA can be obtained. Thus, the magnetic flux densities
Similarly, the magnetic flux density in other sections of the magnetic circuit can be obtained, and the magnetic field strength can be calculated by using the magnetization curve of the materials. In this design, it is necessary to ensure that there is no magnetic saturation in each part of the magnetic circuit, and the magnetic field should be concentrated in the two MR fluid gaps. The magnetic circuit model of the designed MR coupling provides a reference for the design of the key dimensional parameters of the multi-pole MR coupling.
Magnetic analysis of the multi-pole bilayer MR coupling
In order to evaluate transmission torque of the multi-pole bilayer MR coupling with TPMA and HPMA, finite element analysis is employed to solve the magnetic circuit. Firstly, the values of structural parameters are carried out according to the magnetic circuit modeling which is shown in Section 4, and a 3D simulation model of the multi-pole bilayer MR coupling has been built. Then, the magnitude and direction of the residual flux density of each permanent magnet are set. After the final mesh is divided, the simulation analysis is carried out under the steady state. The values of structural parameters are shown in Table 2.
Parameters and values
Parameters and values
The total transmission torque of the multi-pole bilayer MR coupling with TPMA obtained from the simulation results is 47.02 N m, while that of the HPMA is 70.65 N m, when the residual magnetic flux density of the permanent magnets is set to 0.8 T. The simulation results show that the designed MR coupling with HPMA which has six poles can significantly improve the transmission torque compared to the same size of the MR coupling with TPMA.
Figure 4 shows the crosssectional magnetic flux density and flux direction of the two multi-pole bilayer MR couplings with TPMA and HPMA. From the simulation results, it is seen that the direction of the magnetic flux is the same as the equivalent magnetic circuit in the magnetic circuit analysis. This result verifies the correctness of the magnetic circuit design. For the multi-pole bilayer MR coupling with TPMA, as shown in Fig. 4(a), most of the magnetic flux passes through the adjacent permanent magnets because the magnetic permeability of the permanent magnets is high and the direction of the magnetic poles of two adjacent permanent magnets are opposite. Due to the distribution of the magnetic field, its magnetic permeability properties may not be considered when selecting the materials of the spaced block between the two adjacent permanent magnets. The block can be designed as a hollow structure to reduce the weight, thereby the weight of the coupling be reduced, and the torquemass ratio of the coupling could be also improved. For the multi-pole bilayer MR coupling with HPMA, as shown in Fig. 4(b), the magnetic flux density on the outer stator is higher than that of the MR coupling with TPMA, while the magnetic flux density on the inner stator is significantly lower than the MR coupling with TPMA. The results show that the multi-pole bilayer MR coupling with HPMA can well concentrate the magnetic field to one side.

Flow direction of the magnetic flux in the cross section of the multi-pole bilayer MR coupling with: (a) TPMA; (b) HPMA.
The magnetic flux density of the inner and outer MR fluid working gaps is shown in Fig. 5. It is seen that the average flux densities in the inner gap of the MR coupling with TPMA and HPMA are 0.99 T and 0.48 T, respectively, as shown in Fig. 5(a). The arc length is the product of radian and radius. The magnetic flux density within the inner gap of the multi-pole bilayer MR coupling with TPMA is larger than that MR coupling with HPMA. This is due to that the sum of the magnetic resistance of the magnetic circuit nodes a ′ b ′ , g ′ h ′ and a ′ h ′ in HPMA is larger than that of the magnetic circuit node i ′ j ′ , which results in less magnetic flux passing through the inner gap in HPMA. Figure 5(b) shows that the average flux densities in the outer gap of the MR coupling with TPMA and HPMA are 0.63 T and 0.98 T, respectively. The magnetic flux density within the outer gap of the multi-pole bilayer MR coupling with TPMA is smaller than that MR coupling with HPMA. This is due to that less magnetic flux passes through the inner gap and inner stator in the multi-pole bilayer MR coupling with HPMA, while the flux density in the outer gap of MR coupling with HPMA is greater than that of TPMA.

Magnetic flux density in the MR fluid working gaps. (a) Inner gap; (b) outer gap.
To sum up, the magnetic field strength in the outer gap of the MR coupling with HPMA is larger. It can be seen from Eq. (3) that the radius of MR fluid is positively correlated with the generated torque. Therefore, increasing the magnetic field strength of the outer gap is beneficial to increase the total torque. Its unilateral magnetic focusing effect in the outer gap is obvious, and the transmission torque of the coupling is effectively improved.
The comparison results between the calculated value and the simulation value of the gap magnetic flux density of the magnetorheological fluid of the two schemes are shown in Table 3.
Simulation and calculation values of magnetic flux density within the MR fluid gaps
The magnetic flux density on the cylindrical surfaces of the inner and outer MR fluid gaps along the axial direction is shown in Fig. 6. The distribution of magnetic flux density along the axial direction is the same, which makes the torque evenly distributed on the circumference. Due to the aluminum block and the secondary magnetic pole, lower magnetic flux density is observed in the areas where the spacer block or the secondary magnetic pole is located (dark blue region) while higher magnetic flux density is observed in the areas near the main permanent pole (red or yellow region).

Magnetic flux density distribution in the MR fluid gaps: (a) TPMA; (b) HPMA.
Furthermore, the influence of residual magnetic flux which generated by the secondary pole of the permanent magnet is studied. As shown in Fig. 7, with the increase of the residual magnetic flux density of the secondary magnetic pole, the torque generated by the inner gap decreases. This is because that the magnetic field strength in the inner gap decreases. It is also seen that the torque generated by the outer gap increases with the increase of the residual flux density of the secondary magnetic pole. This is mainly due to that the flux passing through the outer gap increases, and the magnetomotive force in the magnetic circuit increases with the increase of the residual flux density of the secondary pole. The magnetic field in the inner gap decreases gradually, which also indicates that the magnetic field gradually accumulates in the outer gap. The magnetic focusing effect can be defined by the ratio of the torque generated by the inner gap to the torque generated by the outer gap. The smaller the ratio is, the better the magnetic focusing effect is. From the simulation results, with the increase of the residual flux density of the secondary pole, the unilateral magnetic focusing effect of the multi-pole bilayer MR coupling with HPMA is more obvious.

Effect of residual flux density on transmission torque.
Figure 8 shows the influence of the inner diameter of the permanent magnet which is also the radius of the inner gap on the transmission torque of the MR coupling. In Fig. 8(a), the torque generated by the inner gap of the MR coupling with TPM increases with the increase of the inner diameter of the permanent magnet. This is due to that the working radius of the inner gap increases when the inner diameter of the inner stator increases. It can be seen from Eq. (4) that the torque increases with the increases of the radius of the inner gap. And the inner diameter of the inner stator has little effect on the transmission torque generated by the outer gap.

Effect of inner stator diameter on transmission torque: (a) TPMA; (b) HPMA.
In Fig. 8(b), the torque generated by the inner gap of the MR coupling with HPMA gradually increases with the increase of the inner diameter of the inner stator, which is also due to the increase of the working radius of MR fluid gap. However, the torque generated by the outer gap decreases with the increase of the inner diameter of the inner stator, which is also due to the decrease of the magnetomotive force of the main pole. It results in the decrease of the magnetic field strength in the outer gap. It is found that the total torque decreases with the increase of inner diameter of the inner stator. With the increase of inner diameter of the inner stator, the ratio of the torque generated by the inner gap to the outer gap gradually increases. It indicates that the unilateral magnetic focusing effect of Halbach permanent magnet array (HPMA) becomes weaker with the decrease of radial length of permanent magnet. The results show that the longer the radial length of the permanent magnet is, the better the magnetic field will be gathered on the outer gap side is.
The relationship between the transmission torque and the angle of the main pole of the permanent magnet is shown in Fig. 9. From the Fig. 9(a), it is seen that the torques generated by the inner gap, the outer gaps and the total torque of the MR coupling with TPMA increase with the increase of the angle of the main magnetic pole. According to the Eqs (11) and (13), the magnetic flux in the magnetic circuit increases with the increase of the main pole angle. Therefore, both the torques generated by the inner gap and the outer gap will increase, as well as the total torque.

Transmission torque versus the main magnetic pole angle: (a) TPMA; (b) HPMA.
As can be seen from Fig. 9(b), the torque generated by the inner gap of the MR coupling with HPMA first decreases and then increases with the angle of the main magnetic pole, and it reaching the minimum value 4.63 N m at 15°. This is due to that as the angle of the secondary magnetic pole decreases with the increases of the angle of the main magnetic pole. This makes the reluctance of the magnetic circuit node i ′ j ′ decreases gradually, and the flux passing through the inner gap decreases. However, the magnetic flux at the main magnetic pole increases with the increase of its angle, which makes the magnetic field in the inner gap begin to increase after 15°. It is also found that the torque generated by the outer gap first increases and then decreases with the increase of the main magnetic pole angle, and it reaches the maximum value 74.56 N m at 50°. This is due to that as the angle of the main magnetic pole increases, the magnetic field passing through the inner gap decreases, and the magnetic flux in the magnetic circuit increases as the angle of the main magnetic pole increases. However, at the same time, the magnetomotive force at the secondary magnetic pole decreases when the angle of the secondary magnetic pole decreases, and the magnetic flux at the secondary magnetic pole will pass through the outer gap, which makes the torque generated by the outer gap begin to decrease after 50°. The total torque increases with the increase of the main pole angle. The ratio of the torque generated by the inner to the outer gap first decreases and then increases with the increase of the angle of the main magnetic pole. When the main magnetic pole is 15°, the ratio is the smallest, and the magnetic gathering effect is the best. In this case, a single gap MR coupling design is suitable, which can further simplify its structure, and the loss of torque is smaller than that of the bilayer gap design.
The purpose of this paper is to improve the transmission torque of the multi-pole bilayer permanent magnets MR coupling by enhancing the magnetic field strength in the MR fluid gaps. Firstly, two structures of the multi-pole bilayer MR with traditional permanent magnet array (TPMA) and Halbach permanent magnet array (HPMA) are introduced. Then, the transmission torque of the MR coupling is deduced. Based on the electromagnetic field theory, the magnetic circuit models are established, and the magnetic flux density within the MR fluid working gaps has been derived in order to evaluate the merits of the designed MR coupling. Finally, a 3D FE magnetic-fluid analysis has been conducted. The results show that the torque generated by the multi-pole bilayer MR with HPMA is 33.45% higher than that of the MR coupling with TPMA when six main magnetic poles are set.
For the multi-pole bilayer MR with TPMA, the transmission torque generated by the inner and outer gaps increases with the increase of the main magnetic pole angle, and the inner diameter of the inner stator has little effect on the transmission torque generated by the outer gap.
For the multi-pole bilayer MR with HPMA, the bigger the residual magnetic flux density of the secondary magnetic pole and the radial length of the magnetic pole are, the better the unilateral magnetic concentration effect is. In this design, the transmission torque generated by the inner gap decreases first and then increases with the angle of the main magnetic pole, while the transmission torque generated by the outer gap increases first and then decreases. The results show that the torque generated by the outer gap reaches the maximum value 74.56 N m when the main magnetic pole is 50°. When the angle of the main magnetic pole is 15°, the single side magnetic gathering effect is the best.
Footnotes
Acknowledgements
This work supported by Research Program supported by the National Natural Science Foundation of China (51805444), the Young Scholars Reserve Talents Support Program of Xihua University (DC1900007176) and Talent Introduction Project of Xihua University (Z201017).
