Research on the Flow and Transmission Performance of Magnetorheological Fluid between Two Discs

: The viscoplastic flow of magnetorheological fluid in a disc was analyzed based on the Navier–Stokes Momentum Equation, and the yielded and unyielded decomposition surfaces were obtained. For the shear ‐ thinning phenomenon of magnetorheological fluid, the magnetorheological properties of the magnetorheological fluid were described based on the Herschel–Bulkley model. Then, the relationship between torque, magnetic field, material, size and motion was established, and the magnetic field and working gap were optimized and analyzed. The results show that in the unyielding region, the magnetorheological fluid flows rigidly. In the yielding region, it flows as a viscous fluid. The degree of error of the proposed torque equation decreased gradually with the increase to current, as observed by experimental comparison.


Introduction
Magnetorheological fluid (MRF) is a solid-liquid two-phase smart material whose rheological properties are controlled by an applied magnetic field. It consists of a base fluid, additives and micron-sized magnetic particles uniformly distributed in the base fluid [1]. MRF can complete the reversible transition from a liquid to a solid-like state in milliseconds under a magnetic field. The shear yield stress and apparent viscosity of the material can be continuously controlled with the applied magnetic field [2]. As the magnetorheological effect transformation is reversible and easy to control, rapid response and other characteristics give MRF in the field of transmission an extensive range of application prospects [3]. Magnetorheological (MR) clutches, MR dampers, MR brakes, etc. are a new class of power transfer devices based on the rheological effect that use MRF as the power transfer medium and regulate torque through an applied magnetic field to achieve the purpose of output control [4,5].
Theoretical research and the structural designs of MR actuators are the focus of researchers currently exploring the practical engineering applications of MRF. Choi et al. [6] introduced the torque characteristics of electrorheological and magnetorheological clutches and conducted experiments on the response speed and durability of disc-type MR clutches to verify the feasibility of an MR transmission by comparing the experimental data of two types of clutches. Huang et al. [7] established a relationship expression between the magnetorheological effect generated by a cylindrical MR brake, the required volume of MRF and the minimum gap thickness, to provide a theoretical basis for the design of a cylindrical MR drive. Farjoud et al. [8,9] deduced the torque equation of a cylindrical MR brake at a high shear rate based on the Herschel-Bulkley model and determined the plug flow region in the working gap of MRF through modeling to obtain the relationship between the plug flow region and the viscous torque. Fernandez et al. [10] optimized the structural design of an MR clutch which used permanent magnets by means of a finite element simulation and adopted a multifaceted magnetic field blocking mechanism which significantly improved the torque transfer capability of the transmission. Rizzo [11] proposed a magnetic design method for MR clutches by applying nonlinear materials to amend the magnetic axial force direction of a permanent magnet system. Thakur et al. [12] performed magnetic field analysis for radial MR clutches with different slot types, summarized the relationship between torque performance, disc radius and magnetic field based on their experimental results, and optimized the design of the working gap. Zhang et al. [13] proposed a torque estimation method for disc-type MR devices that considered the centrifugal effects and introduced radial stresses into magnetorheological constitution analysis to improve the accuracy of torque evaluation. Wang et al. [14] presented a computational fluid dynamics simulation model for simulating the behavior of a wet MR clutch and validated the model with experimental data to complement the hardware-based MR clutch design process.
The above research was focused on discussing the influence of a single factor on the torque of MR device, but rarely analyzed the flowing nature of MRF, and rarely considered the shear-thinning phenomenon of MRF in its working process. In terms of device optimization, the analysis of the influence of multiple factors on the transfer torque is rarely considered. Especially in the question of the optimal arrangement of the magnetic field of MR devices, the magnetic field distribution in the working gap of an MRF is highly uneven, resulting in a low MR transfer torque. In order to solve the above-mentioned problems, the viscoplastic flow of the MRF in a disk-type gap was analyzed, and the torque equation was established considering the shear-thinning phenomenon of MRF. Based on electromagnetic field theory and the finite element method, the distribution of the magnetic field in the disc-type working gap of MRF was analyzed. The optimal magnetic field design scheme was obtained through structural adjustment. This article provides a valuable reference for the design and manufacture of disc-type MR transmission devices.

Design of an MR Transmission Device
The working principle of the disc-type MR transmission device is shown in Figure 1. The support shell is fixed with the frame by screws to ensure that the coil does not rotate relatively, and the thickness of the air gap between the support shell and the output disc needs to be taken into account. Copper gaskets and sealing packing in the transmission device are magnetic field-blocking parts that are used to improve the magnetic flux distribution. The MR transmission device work process is as follows: (1) In the initial state, the drive shaft starts to rotate under the action of the input torque from the power source. At this time, the excitation coil is not energized, and the torque transmitted by the zero-field viscosity of the MRF alone is not enough to drive the output disk to rotate. (2) After the excitation coil is energized, the magnetic flux generated by the coil passes through the gap of the MRF perpendicularly, the magnetic particles in the MRF are arranged into a chain structure along the magnetic flux direction, and the shear yield stress gradually increases, driving the output disk to start rotating. When the current increases, until the MRF reaches magnetic saturation, the device transmission torque reaches the extreme value and the driven-disc rotation remains stable. (3) In the event that the system needs to stop the power input, the coil is disconnected, and the shear yield stress of MRF decreases rapidly. Finally, the input disk and the output disk are disconnected, and the output disk stops rotating. Figure 2 shows the shear model between disks and MRF. Consider two discs with radius r and distance h (h<<r), where the working gap between the discs is filled with MRF. The right disc rotates along the z-axis with angular velocity ω 1 and the MRF between the discs is sheared, which drives the output disc to rotate with angular velocity ω 2 . As the working disk in the sudden change of speed, the impact of the working disk rotation on the fluid with time changes shows a certain hysteresis. The inertial force and vibration of the MRF during the drive process have some influence on the flow of the MRF, which increases the uncertainty of the MRF transient analysis [15]. Therefore, to negate the effect of response time on MRF flow, the steady-state assumption is used to analyze the parallel flow of the MRF in the working gap between the two disks. The following conditions are assumed: (1) the MRF is considered incompressible; (2) the flow of the MRF is steady-state; (3) there is no flow of the MRF in the radial direction; (4) the flow of the MRF is only a function of the radius; (5) ignore the volume force; (6) the magnetic field strength is uniformly distributed in the working gap; (7) the pressure in the MRF is along the thickness direction invariant; and (8) ignore the inertial force.

Analysis of MRF Flow
According to the above assumptions, in the cylindrical coordinate system (r, θ, z), the distribution of the flow velocity is assumed to be where v r , v θ and v z are the velocities along the radius direction, rotation direction and thickness direction, respectively, and ω(z) is the angular velocity of the MRF in the thickness direction, which is a function of z.
Assume that the pressure gradient along the rotation direction is dp/dθ. For constant flow, ∂v θ /∂t = 0, there is no change in the r direction. The flow equation along the θ direction can be expressed as [16] where η is the apparent viscosity of the MRF, which is related to the magnetic field strength H, and can be expressed as: Since the two flat plates are parallel, h = h 0 , and neither viscosity η nor (ω 1 + ω 2 ) is equal to zero, it can only be the case that dp/dθ = 0. When the magnetic field strength H varies, the magnetic shear yield stress of the MRF is τy(H). The interface of σ zθ = τy(H) is defined as the decomposition surface of the yielded and unyielded regions. Assuming that zy is the dimension at the decomposition plane between the yielded and unyielded regions, the Navier-Stokes Momentum Equations for the MRF flowing between the two slit flat plates as circumferential shear are [17] (pp. 115-118).
where σ θθ is the tangential pressure along the θ direction and τ rθ and τ zθ are the shear stresses of the MRF.
Since the angular velocity is constant along the flow direction, τ rθ = 0. Since the two flat plates are parallel to each other, the pressure σ θθ is constant, and the above equation can be simplified as From Equation (11), it can be determined that τ zθ is a constant, independent of z. Suppose the load torque on the driven shaft is T L , and the shear stress of the MRF required to make the output disk rotate is τ L , then By the Integral Median Theorem, the integration of Equation (12) is where τ* is the value of the shear stress at point r 0 in the interval (r 3 , r 2 ), and when (r 2r 3 )/r 3 <<1, τ* is approximately equal to the value at r 0 = (r 2 + r 3 )/2. Based on the Bingham model, the dynamic shear yield stress after the yielding of the MRF is τ* = τy(H) + η 0 γ , where η 0 is the zero-field viscosity at H = 0, and the shear strain rate γ ̇ is denoted as [18] In the yielding region, the angular velocity is obtained from Equations (13) and (14) as Integrating Equation (15) and taking the boundary conditions when In the unyielding region, the angular velocity is Taking the boundary conditions when z = zy, ω(z) = ω 2 , the dimensions at the decomposition surface of the yielded and unyielded regions can be obtained as Meanwhile, the output angular velocity of the disc-type MR drive can be expressed as As shown in Equation (19), when the structure, material and load of the disc-type MR drive are fixed, the output angular velocity of the device increases linearly with increasing shear yield stress.

Analysis of Torque
Due to the fact that the Bingham model can only effectively describe the rheological properties of the post-shear yield region and cannot describe the shear thinning phenomenon, in order to describe the phenomenon of the shear thinning or thickening of the MRF viscosity during shear rheology more accurately, the 1-D flow model of the MRF can be described by the Herschel-Bulkley model [8] where m and k are constants greater than zero, which are determined experimentally, and m is the flow behavior index of the MRF. When m = 1, the Herschel-Bulkley model degenerates to the Bingham model; when m < 1 (or m > 1), the shear thinning (or thickening) of the reacting liquid is observed. The equivalent zero-field viscosity in the Herschel-Bulkey model is For the flow of the MRF between two discs, assume that the working area of the disc is a ring of radius from r 3 to r 2 , and take a micro-circular area dA at radius r; since dA = 2πrdr, the torque M transferred by the whole disc can be expressed as From the Herschel-Bulkley model shown in Equation (20), the torque M generated by the magnetorheological fluid after yielding can be composed of two parts, which are the magnetic torque M H generated by the MRF magnetic shear-yielding stress and the viscous torque M η generated by the MRF viscosity.
Combining Equations (20) and (22), the torque M H transferred due to the shear yield stress generated by the chain-formation mechanism of the MRF is Combining Equations (20)(21)(22), the viscous torque Mη transferred due to the shear yield stress generated by the viscosity of the MRF is Therefore, the torque M transferred by the MRF between the two discs can be expressed as

Magnetorheological Fluid Properties
The MRF used in this article was MRF-1 [18] and the material parameters are shown in Table 1. The relationship between the shear yield stress τy(H) and the magnetic field intensity H is derived from the experimental data, as shown in Figure 3.  From Figure 3a, it can be seen that the shear yield stress grows approximate linearly along with the increase of magnetic field intensity when the magnetic field intensity is lower than 100 kA/m. MRF starts to express the magnetic saturation phenomenon when the magnetic field exceeds 100 kA/m, at which point the shear yield stress changes more flatly with the magnetic field. The relationship between the magnetic field intensity H and the magnetic shear yield stress τy(H) is established from Figure 3 and fitted to the corresponding equation as According to Figure 3b, it can be concluded that the viscosity of MRF decreases with an increase to the shear rate, and the phenomenon of shear thinning occurs. M = 0.75 and k = 0.085Pa•s m were derived from the experimental data, as the speed error ∆ω was 10 rad/s.

Simplified Model of the Magnetic Circuit
According to the magnetic Ohm's Law [19], we performed a magnetic field analysis of the working gap of the transmission device. The simplified equivalent model of the transmission device and the equivalent model of the magnetic circuit of radii r 1~r7 , the widths L 1~L4 , magnetoresistances R 1~R9 , etc., are shown in Figure 4 and Table 2.   The input disk, output disk and support shell of the MR drive are made of Steel-1020, the material of the drive shaft and driven shaft are Q235, the surrounding is set as an air boundary. and the thickness of the working gap is h. The thickness of the air gap h' at R 3 and R 8 was 0.35 mm. According to the magnetic Ohm's Law, the magnetomotive force F m generated by the excitation coil is [20] where N is the number of turns of the coil (t) and I is the current passed into the excitation coil (A). The ratio of the magnetomotive force F m to the magnetic flux Φ is defined as the magnetic reluctance R.
From Figure 4, the magnetic circuit can be divided into a total of eight parts, labelled 1~9 [19], and the magnetoresistance of each part is where μ 1 and μ 2 are the absolute permeability of Steel-1020 and MRF, respectively, and μ = μ 0 μ r .
The total magnetoresistance of the magnetic circuit R M is obtained by summing up the magnetoresistance of each segment as follows Therefore, the total magnetic flux of the whole magnetic circuit is According to Kirchhoff's Law, the magnetic flux density in the working gap is Substituting Equation (30) into Equation (31), the relationship between magnetic induction and current can be obtained as   where N is the number of turns of the excitation coil and μ 0 is the vacuum permeability.
Taking μ 0 = 4π × 10 −7 N/A 2 , μ r is the relative permeability of the MRF, which is related to the volume fraction, diameter size and temperature of the magnetic powder particles in the MRF, and which can be measured experimentally [21,22]. Combined with Equations (19) and (33), the change to the output angular speed is positively correlated with the increase to the current when the load torque is constant, and the MR drive output speed is at its maximum when the current increases to the point of magnetic saturation of the MRF [23]. Define the distance between the yielded and unyielded partition interface of the magnetorheological fluid to the input disk as the yielded thickness hτ. The working gap thickness of MRF h is set to 1.0 mm, and the maximum transferable torque is 130 N•m. From Equation (13), it can be known that when the transmission is at full load and the speed difference ∆ω is 0 rad/s, the required minimum MRF shear yield stress τ min can be expressed as   when the MRF is in a non-yielding state, zy = 0, τmin = 37.54 kPa. Therefore, when the magnetic shear yield stress is 38 kPa, the load torque T L and the yield thickness of MRF h τ are shown in Figure 5.  As can be seen from Figure 5, when the load torque is 136.1~138.2 N•m, the MRF transitions from a rigid-plastic flow state to shear yield flow state, and the yield thickness of the MRF multiplies from 0.18 mm to 1 mm; the MRF is wholly in shear yield at T L = 138.2 N•m. The partial yielding flow of the MRF occurs in a narrow range as a transition between the unyielding and fully yielding states.

Analysis of Magnetic Field
From Equations (26) and (33), it can be seen that the magnetic field intensity is related to the current, the working gap size, and the permeability of materials. The calculations were performed for the above variables to determine the magnetic field strength in the working gap.
The minimum magnetic field intensity H min = 67.82 kA/m within the MRF gap under a rated load was obtained by bringing τ = 38 kPa into Equation (26). Assuming that no magnetic saturation occurs in all parts of the MR drive, and the materials of the input disc, output disc and support shell are Steel-1020, the relative permeability of the housing and MRF are 1200 and 5, respectively. From Equation (29), the magnetoresistance is 1.06 × 10 5 H −1 . The design value of the coil parameter NI is 1500 A•t for a magnetic induction intensity H greater than 67.82 kA/m in the working gap of the MRF.
The relative magnetic permeability μ r of the material changes with the magnetic field intensity H when the current of the excitation coil changes. In order to improve the accuracy of the magnetic field analysis, the Finite Element Method (FEM) was used to simulate the magnetic field in the whole working gap. Set the number of turns of the coil N = 500 t, the coil was placed through the currents from 1~4 A, respectively. The 2D axisymmetric model in Figure 4 was divided into 8909 cells using a 2D triangular mesh. We analyzed the influence of the copper gaskets on the transmission using the FEM, obtaining the results shown in Figures 6 and 7. As presented in Figure 6 and Figure 7, the overall magnetic flux distribution within the working gap of the MR transmission is uniform, and the working gap can be divided into four segments according to the variations in the magnetic field intensity. When the working gap radial distance r is from 0 to 5 mm, the magnetic field intensity increases rapidly to extremum with increases to the radial length. As the working gap radial distance is located in the two segments of r = 5~20 mm and r = 30~45 mm, the magnetic field intensity curve changes relatively smoothly and the magnetic field is more stable. With the radial distance from 20 mm to 30 mm, the magnetic field intensity curve decreases rapidly to 10 kA/m in the packing seal barrier interval. While r ranges from 45 to 50 mm, the magnetic field declines exponentially from the stable value.
In the absence of copper gaskets, the magnetic field in the working gap of the MRF on both sides is asymmetric. In summary, for optimal torque characteristics, the working gap thickness of MRF h is 1.0 mm, the number of turns of the coil N is 400 t, and the current I is 0~3 A. These measurements were used as as the final optimization scheme of the transmission device.

Flow Analysis of Magnetorheological Fluid
According to Equations (3) and (20), the apparent viscosity η of the MRF after yielding is From Equation (35), it can be observed that the apparent viscosity of MRF is inversely proportional to the shear strain rate. Along with the increased magnetic field intensity, the effect of zero-field viscosity η 0 on the apparent viscosity becomes smaller. When the magnetic field intensity is 10 kA/m, 30 kA/m, 50 kA/m and 70 kA/m, respectively, the relationships between apparent viscosity and the shear strain rate of the MRF were as shown in Figure 8.
As can be seen from Figure 8, along with the increase of magnetic field intensity, the effect of the shear strain rate on the apparent viscosity increases. The value of the apparent viscosity tends to stabilize when the shear rate is greater than 800 s −1 . To facilitate the analysis of the plug flow region (the unyielding region of the MRF), the ratio of the minimum shear yield stress τ min required by the load to the MRF magnetic shear yield stress τy(H) is defined as the plug flow coefficient λ. Then, λ can be expressed as The magnetic field intensity data at 3A in Figure 7a was brought into Equations (26) and (36) to obtain the plug flow coefficient in the working gap, as shown in Figure 9, where hz is the distance across the gap, which indicates the distance between the plane where the MRF is located and the input disc.
As can be seen from Figure 9, when λ ≥ 1, the magnetorheological fluid is in the unyielding state, and it flows rigidly in the working gap. Since λ << 1 near the wall of the working gap of the transmission, the magnetorheological fluid near the wall of the working gap is in the yielding flow state. According to the plug flow coefficient in Figure 9a, the flow velocity distribution of the MRF in the working gap was plotted as Figure 9b. From Figure 9b, it can be seen that the rotational speed of the disc-type magnetorheological transmission with the coil placed at one side is more related to the thickness of the gap, and the flow velocity distribution of the MRF in the unyielding region (regions A and B) is closer to the rotational speed of the input disc. In regions E, C and D, the magnetic field intensity is low, which leads to the yielding of the MRF, and the MRF is in a viscoplastic flow state. The flow velocity distribution in the cross section where the radius of disc r is 104 and 81 mm is intercepted in regions E and C, respectively, as shown in Figure 10. The flow velocity of MRF varies widely along the thickness direction in the E, B, and C regions, the flow velocity does not decrease according to the linear law because the viscosity of the MRF undergoes shear thinning with the increase to the shear strain rate. The MRF in this region is not fully cured, which is one of the reasons for the unstable transmission performance.
According to the analysis of Figure 5, when the load torque is further increased, the magnetorheological fluid in the working gap changes from a rigid to a viscoplastic flow state. The load torque of the transmission was set to 137.7 N•m, the speed difference ∆ω = 22.74 rad/s was known from Equations (18) and (19), and the yield thickness was 0.52 mm. The flow of MRF in the working gap is shown in Figure 11. As can be seen from Figure 11, the magnetic field intensity near the input disc was slightly greater than near the output disc when the shear yield stress required by the load exceeded the magnetic shear yield stress of the MRF. The MRF thus yields to flow near the output disc, and a rotational speed difference will appear. Figure 11b shows the flow velocity of the MRF in the cross section where the radius of the disc r was 89 mm, and the MRF in the yielding region was in a laminar flow state.
When hz > 0.52 mm, the MRF in the A region yields everywhere, and the flow rate of the MRF at this time is When hz ≤ 0.52 mm, the magnetorheological fluid in the A region does not yield, and the magnetorheological fluid in in a solid-like flow state. At this time, when the flow velocity of the magnetorheological fluid is ω(z) = ω 1 , the input disk and the output disk rotate synchronously.

Magnetorheological Torque Test
An MR transmission performance test bench was built to verify the correctness of the torque equation, as shown below ( Figure 12 The primary devices involved in the test system included a coupling, a ZH07 speedtorque sensor (range 500 N•m, indexing value 0.01 N•m, maximum speed 4000 rpm), a T16-1-L bevel gear reducer (transmission ratio i = 3), an MCK-S dynamic torque digital display meter, a magnetic powder brake, an RXN-3010D digital DC power supply, a YQ-3000-V7 inverter speed regulator, a 1LE0001-1DA4 three-phase asynchronous motor (power 18.5 kW, rated speed 2930 rpm), and a ZL3A-5S torque controller. The type of magnetorheological fluid was MRF-1 [18]. The internal structure of the MR drive is shown below ( Figure 13):

Analysis of Experimental Results
During the experiment, different currents were passed into the excitation coil of the MR drive and the braking torque of the magnetic powder brake was continuously adjusted. The relationship between the current and the torque transmitted by the MRF was derived from the value measured by the torque sensor.
The Finite Element Analysis (FEA) results in Figure 7a were substituted into Equation (25), and the values obtained were r 2 = 0.129 mm and r 3 = 0.079 mm. The input disk speed ω 1 was 100 rad/s, the output disk speed ω 2 was 90 rad/s, the Herschel-Bulkley model coefficient m is 0.75, k was 0.085 Pa•s m and the thickness of the working gap h was 0.001 m.
The relationship between the current passed into the excitation coil and the torque of the MR drive are shown in Figure 14. As can be observed from Figure 14, the MR torque can be divided into roughly two stages according to the rate of change of the transmitted torque.
The transmission torque increased significantly with the current when the current was between 0.5 A and 1.5 A. When the current reached 0.5 A, the experimentally derived torque was 26.57 N•m, and the theoretically calculated torque was 48.99 N•m. The tested torque was 82.39 N•m when the current was 1.0 A, and the theoretically calculated torque was 93.10 N•m, with a transmission torque error of 11.50%.
The current I was in the range of 1.5 A to 3.0 A in the second stage. At this time, the torque of the transmission device increased slowly with the current. When the current was 3.0A, the experimentally derived torque was 129.84 N•m, and the theoretically calculated torque was 133.19 N•m with a torque error of 2.52%.
Due to problems such as the settling of the magnetic particles in the MRF [24], especially when the current is low, the magnetic force between magnetic particles decreased, the zero-field viscosity of the magnetorheological fluid and the shear yield stress under the action of the magnetic field decreased, the magnetorheological fluid performance was unstable, and the magnetic field distribution in the working gap was not uniform, resulting in a lower torque performance of the transmission than the calculated value.
The torque sensor was installed between the output disc and the magnetic powder brake, and the braking torque of the magnetic powder brake was set to 125 N•m. The input disc speed ω 1 was controlled by the inverter and stabilized at 100 rad/s. The relationship between the current of the coil I and the output disc speed ω 2 was determined according to the value taken by the speed-torque sensor and compared with the theoretical result of Equation (19).
As can be observed from Figure 15, when the current increased to 2.3 A, the output disk started to rotate under the shear yield stress of the MRF, and it grew linearly with the increase of current, with the yield flow area of the MRF decreasing. Because of the fluctuation of the motor speed and additional friction during the experiment, there were some errors between the experimental and theoretical values, but the overall trend was consistent. The results indicate that the magnitude of the current significantly affects the angular velocity distribution of the flow of MRF in the working gap.

Conclusions
In this research, magnetic field gap optimization and torque equation derivation were carried out for a disc-type MR drive. The equivalent reluctance was solved based on the magnetic Ohm's Law, and the relationship between the magnetic field strength, coil current, structure size and material was derived. The magnetic flux distribution was analyzed by the FEA, and the maximum magnetic field strength in the working gap of the MRF was improved by the reasonable arrangement of magnetoresistive elements.
The working principle of the drive was elaborated, the flow velocity distribution of MRF in the working gap was analyzed, the output disc speed equation was derived based on the momentum equation, the torque equation of the disc MR drive was derived based on the Herschel-Bulkley model, and the shear torque generated by MRF at different currents was calculated. This research shows that: (1) When the coil is energized with DC current, the magnetic field strength under the steady-state magnetic field changes with the variations of the current, working gap size, and material magnetization, satisfying a cubic polynomial relationship with the shear yield stress of the MRF. (2) Magnetic field optimization through the reasonable use of magnetoresistive elements can effectively increase the torque performance of the MR drive and reduce the influence of the working gap thickness on the magnetic field. Under the same current, the magnetic field strength in the gap of the MRF can be enhanced by up to 45.71% through magnetic field optimization. (3) When the MRF has a continuous circumferential flow between two slit plates, the output angular velocity of the disc-type MR drive is positively correlated with the change in current, and the yield region of the MRF expands with the increase of the load torque, but does not satisfy the linear relationship. (4) When the torque of the drive is greater than the rated design value, the MRF rapidly changes from rigid flow to full-yield flow, and the partial yield flow state of the MRF exists only in a narrow torque zone. The flow mode of MRF in the yield flow region is laminar flow, and the flow velocity varies greatly with the thickness direction. (5) The trend of the experimentally tested curve was generally consistent with that of the theoretically calculated curve. The experimental comparison shows that the degree of error of the torque equation used in this paper decreases gradually with the increase of current.