# A Novel Integral 5-DOFs Hybrid Magnetic Bearing with One Permanent Magnet Ring Used for Turboexpander

Authors: Bangcheng Han; Xu Liu; Shiqiang Zheng
Journal: Mathematical Problems in Engineering
 (2014)

---

## Abstract

We propose a novel combined five-degrees-of-freedom (5-DOFs) hybrid magnetic bearing (HMB) with only one permanent magnet ring (PMR) used for turboexpanders. It has two radial magnetic bearing (RMB) units; each has four poles and one thrust magnetic bearing (TMB) to control 5-DOFs. Based on one PMR, the bias flux of the two radial magnetic bearing units and the one thrust magnetic bearing unit is constructed. As a result, ultra-high-speed, lower power loss, small size, and low cost can be achieved.  Furthermore, the equivalent magnetic circuit method and 3D finite element method (FEM) are used to model and analyze the combined 5-DOFs HMB. The force-current, force-position, torque-coil currents, the torque-angle position, and the stiffness models of the combined 5-DOFs HMB are given. Moreover, its coupling problems between the RMB units and the AMB unit are also proposed in this paper. An example is given to clarify the mathematical models and the coupling problems, and the linearized models are proposed for the follow-up controller design.


---

## Body

## 1. Introduction

Magnetic bearing can provide many advantages for turbomachinery such as high-speed, long-life, higher efficiency, lower maintenance, and higher reliability. However, the greatest benefit for turboexpanders is realized by the ability of the magnetic bearings to run in the process gas, which in turn eliminates the lube oil and dry gas seals. As a result, there is no lube oil to contaminate the process. Elimination of the lubrication system also provides significant size and mass savings for off-shore applications. The turboexpanders in S2M, GE, SKF, and so forth have been equipped with the magnetic bearings [1–5]. Jumonville et al. provide useful photographs of actual magnetic bearing hardware, loading plots, and general design guidelines. A neon cryogenic turboexpander with magnetic bearing is developed by Hirai et al. [6]. The experimental equipment of active magnetic bearing on a 150 M3 turbo-oxygen gas expander is tested by Xiping et al. [7]. Active magnetic bearing applied on a gas expander is mentioned in the literature [8, 9]. A magnetically levitated 400 kW turbogenerator system for natural gas expansion is developed in the literature [10]. Schmied provides excellent information on magnetic bearing control system in a very understandable fashion [11]. The rotor of the turboexpander mentioned above is supported by active magnetic bearings in 5 degrees-of-freedom (DOFs) (Figure 1), which consists of two radial magnetic bearings with bias current, two thrust magnetic bearings with bias current, expander wheel, displacement sensors, shaft, touchdown bearings, and compressor wheel. The radial magnetic bearing and thrust magnetic bearing with bias current have better control ability and high stiffness, whereas they suffer from high power loss due to the bias current. To reduce the power consumption, hotness, and nonlinearity, the radial and axial hybrid magnetic bearings with permanent magnet providing bias flux for suspension and the electromagnet providing necessary control flux are developed [12–15]. Most of the 5-axis hybrid magnetic bearing (HMB) systems are usually composed of two radial hybrid magnetic bearings (RMBs) and one pair of axial hybrid magnetic bearings (AMBs) [16–18]. These magnetic bearing systems are easy to manufacture, but their disadvantages are high power consumption, bulky, complex process, and high cost. The reason is that each RMB has eight poles that lead to more eddy current loss. Moreover, each of RMBs and AMBs has one permanent magnet ring (PMR) that is used to provide bias flux. As many as four PMRs are required in this 5-DOFs HMB system.Figure 1
The configuration of turboexpander supported by the active magnetic bearings in 5 DOFs. 1: expander wheel; 2: displacement sensor; 3: radial magnetic bearing; 4: thrust magnetic bearing; 5: shaft; 6: touchdown bearing; 7: compressor wheel.In order to reduce the power consumption, cost, and size, it is important to reduce the number of the units by means of combining the RMB and AMB. Integrated axial magnetic bearing and radial magnetic bearing with the conical rotor are designed and analyzed [19, 20] and one downside of the integrated bearings is a strong coupling problem between the radial and axial degrees-of-freedom. In addition, the combined radial-axial hybrid magnetic bearings (CRAMBs) are proposed [21–28]. The CRAMBs mentioned above have been applied in ultra-high-speed machinery, but it may not suit the turboexpanders.In this paper, a novel combined 5-DOFs HMB with only one PMR is proposed to be used for turboexpanders in this paper. The combined 5-DOFs HMB consists of two RMB units each with four poles and one AMB unit. The bias flux of the two RMB units and AMB unit is supplied by only one PMR. The low power consumption, compact structure, and small size can be gained due to the combined 5-DOFs HMB, the permanent magnetic bias HMB, and the reduction in the number of the poles of the RMB units. The low cost is gained due to the one PMR used. Moreover, ultra-high-speed is also gained due to the small outside diameter of the rotor thrust disk. The equivalent magnetic circuit method and 3D finite element method (FEM) are used to model and analyze the combined 5-DOFs HMB. The force-current, force-position, torque-current, the torque-angular position, and the stiffness models of the integral 5-DOFs HMB are given. Moreover, its coupling problems between the RMB units and the AMB unit are also proposed. An example is given to clarify the mathematical models and the coupling problems, and the linearized models are proposed for the follow-up controller design.This paper is organized as follows. First, in Section2, the configuration of combined 5-DOFs HMB is given. Then, the equivalent magnetic circuit method is used to model and analyze the combined 5-DOFs HMB in Section 3. Section 4 is devoted to developing simulation and clarifying the linearized models. Finally, Section 5 concludes this paper. We focus on the configuration, modeling, and analysis of the combined 5-DOFs HMB; a future paper will compare with experimental results.


## 2. Analyzed Configuration

A simplified configuration of turboexpander supported by the combined 5-DOFs HMB is shown in Figure2. The assembly consists of one 5-DOFs HMB, expander wheel, displacement sensor, shaft, touchdown bearing, and compressor wheel. And the rigid body motion on the rotor-bearing system is considered in this paper. The configuration and magnetic circuit of the combined 5-DOFs HMB are shown in Figure 3. The force and torque produced by the RMB units and AMB unit are shown in Figure 4. It consists of two RMB units (RMB units A and B), one pair of AMB unit, and one axially magnetized PMR. The bias fluxes produced by the permanent magnets are shown by the solid arrow lines and the control flux is generated by y-axis (blue color) and z-axis coil currents (red color), which are shown by the dotted arrow lines in Figure 3. Each stator of the RMB unit consists of two soft magnetic pieces each of which has four equally spaced magnetic poles and four coils distributed at intervals of 90 degrees around the stator iron. The control flux produced by the radial and axial coil currents will allow adding or subtracting flux to the radial air gaps and axial air gaps, and then restoring force and torque will be produced to stabilize the rotor in 5-DOFs except the spin axis in relation to the rotor position. The bias flux of the RMB units and AMB unit is produced by only one PMR. The RMB units A and B stabilize rotor on the two axes (x-axis and y-axis) orthogonal to the spin axis (z-axis) and exert restoring torque on the rotor around x-axis or (and) y-axis when subjected to tilting disturbed torque. The AMB unit exerts restoring force on the rotor at z-axis. In a word, the low power consumption, compact structure, and small size can be gained due to the combined 5-DOFs HMB, the permanent magnetic bias HMB, and the reduction in the number of the poles of the RMB unit. The low cost is gained due to only one PMR being used. Moreover, ultra-high-speed is also gained due to the small outside diameter of the rotor thrust disk.Figure 2
The configuration of turboexpander supported by the combined 5-DOFs HMB. 1: expander wheel; 2: displacement sensor; 3: the combined 5-DOFs HMB; 4: shaft; 5: touchdown bearing; 6: compressor wheel.Figure 3
The configuration and magnetic circuit of the combined 5-DOFs HMB. (a) The section view of the RMB unitA. (b) The y-z section view of the combined 5-DOFs HMB.Figure 4
The force and moment produced by the combined 5-DOFs HMB.


## 3. Modeling and Analysis of the Combined 5-DOFs HMB

### 3.1. The Equivalent Magnetic Circuits of the Combined 5-DOFs HMB

According to construction of the combined 5-DOFs HMB, the magnetic resistance of the iron path and flux leakage is neglected, the permeability of the iron is assumed to be infinite, and their equivalent bias and control magnetic circuits can be gained and are shown in Figure5.Equivalent bias circuit and control circuits of the combined 5-DOFs HMB. (a) The equivalent bias circuit. (b) The equivalent control circuit of RMB unitA. (c) The equivalent control circuit of RMB unit B. (d) The equivalent control circuit of AMB unit.
(a)
(b)
(c)
(d)


### 3.2. The Model of Combined 5-DOFs HMB for Translation inx-Axis (or y-Axis) and z-Axis Direction

Figure3 describes the operating principle of this combined 5-DOFs HMB in y-axis and z-axis. The rotor will depart from its center position when subject disturbance force in negative y-axis direction. The control flux on the upper side of the RMB units A and B will add to the bias flux in upper radial air gap and subtract from the bias flux in the lower radial air gap. And the control flux of the AMB unit should also be calculated. A net restoring force on the rotor will be produced in positive y-axis direction.According to the structure of the combined 5-DOFs HMB (Figure4) and its equivalent magnetic circuits (Figure 5), the net force of the RMB unit A in the y-axis direction (N
r
i
A
x
=
0, where N
r is the number of winding turns and i
A
x is the control current of x-axis coil currents), f
y, is
(1)
f
A
y
=
(
ϕ
y
+
A
p
+
ϕ
y
+
A
i
-
ϕ
y
+
A
z
)
2
2
μ
0
A
r
-
(
ϕ
y
-
A
p
-
ϕ
y
-
A
i
-
ϕ
y
-
A
z
)
2
2
μ
0
A
r
,
where the constant μ
0 is the permeability of free space and is equal to 4π × 10−7 H/m; A
r is the area of the pole face of RMB unit; ϕ
y
+
A
p and ϕ
y
-
A
p are the bias flux in radial air gaps of the RMB unit A; ϕ
y
+
A
i and ϕ
y
-
A
i are the control flux in the radial air gaps of the RMB unit A produced by the y-axis coil currents; ϕ
y
+
A
z and ϕ
y
-
A
z are the control flux in the radial air gaps of the RMB unit A produced by the z-axis coil current of the AMB unit.The bias flux in Figure5 can be calculated, respectively, as follows:
(2)
ϕ
y
+
A
p
=
F
m
R
sum
R
A
R
z
1
R
A
+
R
z
1
1
R
y
+
A
,
ϕ
y
-
A
p
=
F
m
R
sum
R
A
R
z
1
R
A
+
R
z
1
1
R
y
-
A
,
ϕ
z
1
p
=
F
m
R
sum
R
A
R
z
1
R
A
+
R
z
1
1
R
z
1
,
ϕ
z
2
p
=
F
m
R
sum
R
B
R
z
2
R
B
+
R
z
2
1
R
z
2
,
where
(3)
R
sum
=
R
m
+
R
A
R
z
1
R
A
+
R
z
1
+
R
B
R
z
2
R
B
+
R
z
2
,
R
A
=
1
(
1
/
R
x
+
A
)
+
(
1
/
R
x
-
A
)
+
(
1
/
R
y
+
A
)
+
(
1
/
R
y
-
A
)
,
R
B
=
1
(
1
/
R
x
+
B
)
+
(
1
/
R
x
-
B
)
+
(
1
/
R
y
+
B
)
+
(
1
/
R
y
-
B
)
.As shown in Figures3 and 5, R
sum is the total reluctance of the bias circuit; R
x
+
A and R
x
-
A are the reluctance of the radial air gaps in x-axis for the RMB unit A; R
y
+
A and R
y
-
A are the reluctance of the radial air gaps in y-axis for the RMB unit A; R
x
+
B and R
x
-
B are the reluctance of the radial air gaps in x-axis for the RMB unit B; R
y
+
B and R
y
-
B are the reluctance of the radial air gaps in y-axis for the RMB unit B; R
z
1 and R
z
2 are the reluctance of the axial air gaps in z-axis for the AMB unit. F
m
=
H
c
h
m, where H
c is the coercive force of the permanent magnet; h
m is PMR thickness; R
m is the magnetic resistance of PMR and it can be given as
(4)
R
m
=
h
m
μ
0
μ
r
A
m
,
where μ
r is the relative permeability of permanent magnet and A
m is the middle-section area of the PMR.The reluctances of radial air gap are written by(5)
R
x
+
A
=
s
r
0
+
x
μ
0
A
r
,
R
x
-
A
=
s
r
0
-
x
μ
0
A
r
,
R
x
+
B
=
s
r
0
+
x
μ
0
A
r
,
R
x
-
B
=
s
r
0
-
x
μ
0
A
r
,
R
y
+
A
=
s
r
0
+
y
μ
0
A
r
,
R
y
-
A
=
s
r
0
-
y
μ
0
A
r
,
R
y
+
B
=
s
r
0
+
y
μ
0
A
r
,
R
y
-
B
=
s
r
0
-
y
μ
0
A
r
,
where s
r
0 is the nominal air gap of the RMB unit; y is the rotor position as deviation from the operating point in y-axis direction; x is the deviation from the operating point in x-axis direction.The reluctances of axial air gap are written by(6)
R
z
1
=
s
z
0
+
z
μ
0
A
z
,
R
z
2
=
s
z
0
-
z
μ
0
A
z
,
where s
z
0 is the nominal air gap of the AMB unit; z is the rotor position as deviation from the operating point in z-axis direction; A
z is the area of the pole face for the AMB unit.The control flux in the radial air gaps ofy-axis, ϕ
y
+
A
i and ϕ
y
-
A
i, generated by y-axis coil currents of RMB unit A (N
r
i
A
x
=
0, where N
r is the number of winding turns and i
A
x is the control current of x-axis coil currents) can be calculated as
(7)
ϕ
y
+
A
i
=
N
r
i
A
y
R
y
+
A
+
R
x
-
A
/
/
R
x
+
A
/
/
R
y
-
A
+
N
r
i
A
y
R
y
-
A
+
R
x
-
A
/
/
R
x
+
A
/
/
R
y
+
A
R
x
-
A
/
/
R
x
+
A
/
/
R
y
+
A
R
y
+
A
,
ϕ
y
-
A
i
=
N
r
i
A
y
R
y
-
A
+
R
x
-
A
/
/
R
x
+
A
/
/
R
y
+
A
+
N
r
i
A
y
R
y
+
A
+
R
x
-
A
/
/
R
x
+
A
/
/
R
y
-
A
R
x
-
A
/
/
R
x
+
A
/
/
R
y
-
A
R
y
-
A
,
where i
A
y is the control current of y-axis coil currents for the same displacement of RMB unit A.The control flux in the radial air gaps,ϕ
y
+
A
z and ϕ
y
-
A
z, which are produced by the z-axis coil currents of the AMB unit can be calculated as
(8)
ϕ
y
+
A
z
=
ϕ
z
i
[
R
m
/
/
(
R
A
+
R
B
)
]
R
A
R
A
+
R
B
1
R
y
+
A
,
ϕ
y
-
A
z
=
ϕ
z
i
[
R
m
/
/
(
R
A
+
R
B
)
]
R
A
R
A
+
R
B
1
R
y
-
A
,
where the control flux in the axial air gap is calculated by
(9)
ϕ
z
i
=
N
z
i
z
R
z
1
+
R
z
2
+
R
m
/
/
(
R
A
+
R
B
)
,
where N
z is the number of winding turns of the AMB and i
z is the control current of z-axis coil current. And ϕ
z
i
1
=
ϕ
z
i
2.Substituting (2), (7), and (8) into the right side of (1) results in the net force on y-axis direction. Using the Taylor series expansion for small values of y and i
y, we can get the following attractive force with linear terms in y-axis direction:
(10)
f
A
y
(
y
,
i
y
)
≅
f
A
y
|
x
=
y
=
z
=
0
i
A
y
=
i
z
=
0
+
k
r
i
i
y
+
k
r
s
y
.Since the net force is zero at the center position in radial and axial directions, (10) can be written by
(11)
f
A
y
(
y
,
i
A
y
)
=
k
r
i
i
y
+
k
r
s
y
,
where k
r
i is the current stiffness of the RMB unit and k
r
s is the negative position stiffness of the RMB unit caused by the increasing attractive force as the air gap is reduced. This must be overcome for stable suspension since the homogeneous solution to magnetic bearing force (MBF) is in the form of hyperbolic functions indicating that y grows with time. They can be calculated as follows:
(12)
k
r
i
=
∂
f
A
y
∂
i
A
y
|
i
A
y
=
i
z
=
0
x
=
y
=
z
=
0
=
2
μ
0
A
r
N
r
H
c
h
m
×
(
s
r
0
2
(
4
A
r
s
r
0
+
A
z
s
z
0
)
h
h
h
h
h
×
(
h
m
μ
r
A
m
+
2
s
r
0
s
z
0
4
A
r
s
z
0
+
A
z
s
r
0
)
)
-
1
,
k
r
s
=
∂
f
A
y
∂
y
|
i
A
y
=
i
z
=
0
x
=
y
=
z
=
0
=
-
2
μ
0
A
r
(
H
c
h
m
)
2
×
(
s
r
0
3
(
4
A
r
s
r
0
+
A
z
s
z
0
)
2
h
h
h
h
h
h
×
(
h
m
μ
r
A
m
+
2
s
r
0
s
z
0
4
A
r
s
z
0
+
A
z
s
r
0
)
2
)
-
1
.According to the structure of the combined 5-DOFs HMB (Figure4) and its equivalent magnetic circuits (Figure 5), the net force of the AMB unit in the z-axis direction is
(13)
f
z
=
(
ϕ
z
1
p
+
ϕ
z
i
)
2
2
μ
0
A
z
-
(
ϕ
z
2
p
-
ϕ
z
i
)
2
2
μ
0
A
z
,
where ϕ
z
1
p and ϕ
z
2
p are the bias flux in axial air gaps of the AMB unit and ϕ
z
i is the control flux in the axial air gaps of the AMB unit by the z-axis coil current.Substituting (2) and (9) for the right side of (13) results in the net force on z-axis direction. Using the Taylor series expansion for small values of z and i
z, we can get the following attractive force with linear terms in z-axis direction:
(14)
f
z
(
z
,
i
z
)
≅
f
z
|
x
=
y
=
z
=
0
i
A
y
=
i
z
=
0
+
k
z
i
i
z
+
k
z
s
z
.Since the net force is zero at the center, (14) can be written by
(15)
f
z
(
z
,
i
z
)
=
k
z
i
i
z
+
k
z
s
z
,
where k
z
i is the current stiffness of the AMB unit and k
z
s is the negative position stiffness of the AMB unit caused by the increasing attractive force as the air gap is reduced. This must be overcome for stable suspension since the homogeneous solution to MBF is in the form of hyperbolic functions indicating that z grows with time. They can be calculated as follows:
(16)
k
z
i
=
∂
f
z
∂
i
z
|
i
A
y
=
i
B
y
=
i
z
=
0
x
=
y
=
z
=
0
=
2
μ
0
N
z
H
c
h
m
×
(
(
4
A
r
s
z
0
s
r
0
+
A
z
)
h
h
h
h
h
×
(
h
m
μ
r
A
m
+
2
s
r
0
s
z
0
4
A
r
s
z
0
+
A
z
s
r
0
)
h
h
h
h
h
×
(
h
m
s
r
0
μ
r
A
m
s
r
0
+
2
A
r
h
m
+
2
s
z
0
A
z
)
)
-
1
,
k
z
s
=
∂
f
z
∂
z
|
i
A
y
=
i
B
y
=
i
z
=
0
x
=
y
=
z
=
0
=
2
μ
0
A
z
(
H
c
h
m
)
2
(
(
h
m
/
μ
r
A
m
)
+
(
2
s
r
0
s
z
0
/
(
4
A
r
s
z
0
+
A
z
s
r
0
)
)
)
2
×
A
z
s
r
0
3
-
s
r
0
2
(
4
A
r
s
z
0
+
A
z
s
r
0
)
s
z
0
(
4
A
r
s
z
0
+
A
z
s
r
0
)
3
.


### 3.3. The Model of Combined 5-DOFs HMB for Tilting

The RMB unitsA and B of the combined 5-DOFs HMB also control 2-DOFs around x-axis and y-axis which are orthogonal to the spin axis (z-axis) when the rotor is subject to disturbance torque around x-axis and y-axis. The RMB units A and B create tilting stiffness to restrain the rotor when subjected to the disturbed torques as shown in Figure 4. According to the structure of the combined 5-DOFs HMB (Figure 3) and its equivalent magnetic circuits (Figure 5), the torque, M
x, acting on the rotor can be given by (17) (N
z
i
z
=
0 and N
r
i
A
x
=
N
r
i
B
x
=
0). Consider
(17)
M
x
=
l
cos
⁡
α
[
(
ϕ
y
+
A
p
+
ϕ
y
+
A
i
)
2
-
(
ϕ
y
-
A
p
-
ϕ
y
-
A
i
)
2
2
μ
0
A
r
hhhhhhhh
+
(
ϕ
y
-
B
p
+
ϕ
y
-
B
i
)
2
-
(
ϕ
y
+
B
p
-
ϕ
y
+
B
i
)
2
2
μ
0
A
r
]
,
where l is the distance from center of the magnetic pole of the RMB units A and B to the center of rotor as shown in Figure 4. The l
*
cos
⁡
α
≅
l for small value of the tilting angle, α. Equation (17) can be simplified as
(18)
M
x
≅
l
[
(
ϕ
y
+
A
p
+
ϕ
y
+
A
i
)
2
-
(
ϕ
y
-
A
p
-
ϕ
y
-
A
i
)
2
μ
0
A
r
]
.The change in radial air gap isl
*
tan
α
≅
l
*
α for small value of the tilting angle, α. The radial air gap resistance can be written as
(19)
R
x
+
A
=
s
r
0
μ
0
A
r
,
R
x
-
A
=
s
r
0
μ
0
A
r
,
R
x
+
B
=
s
r
0
μ
0
A
r
,
R
x
-
B
=
s
r
0
μ
0
A
r
,
R
y
+
A
=
R
y
-
B
≅
s
r
0
+
l
α
μ
0
A
r
,
R
y
-
A
=
R
y
+
B
≅
s
r
0
-
l
α
μ
0
A
r
.Using the Taylor series expansion for small values ofα and i
A
y we can get the following attractive torque with linear terms:
(20)
M
x
(
i
A
y
,
α
)
≅
k
α
i
i
A
y
+
k
α
s
α
,
where k
α
i is the tilting current stiffness and k
α
s is the negative tilting displacement stiffness. They can be calculated as follows:
(21)
k
α
i
=
∂
M
x
∂
i
A
y
|
i
A
y
=
0
α
=
0
=
4
μ
0
A
r
N
r
H
c
h
m
l
×
(
4
A
r
s
z
0
s
r
0
+
A
z
s
r
0
2
s
z
0
h
h
h
h
h
h
×
(
h
m
μ
r
A
m
+
2
s
r
0
s
z
0
4
A
r
s
z
0
+
A
z
s
r
0
)
4
A
r
s
z
0
s
r
0
+
A
z
s
r
0
2
s
z
0
)
-
1
,
k
α
s
=
∂
M
x
∂
α
|
i
A
y
=
0
α
=
0
=
-
4
μ
0
A
r
(
H
c
h
m
)
2
l
2
×
(
s
r
0
(
4
A
r
s
z
0
+
A
z
s
r
0
s
z
0
)
2
h
h
h
h
h
×
(
h
m
μ
r
A
m
+
2
s
r
0
s
z
0
4
A
r
s
z
0
+
A
z
s
r
0
)
2
)
-
1
.


### 3.4. The Coupling Model of the Combined 5-DOFs HMB

#### 3.4.1. The Current Stiffness and Position of the RMB Unit Affected by the Axial Air Gap

The current stiffnessk
r
i and position stiffness k
r
s of the RMB unit in (12) do not consider the effect of the rotor in z-axis direction, namely, the axial air gap. However, the bias flux generated from the PMR flows across the radial air gaps and the axial air gaps (Figure 5(a)); the alternation of the axial air gap will affect the current stiffness and the position stiffness of the RMB and will be calculated by
(22)
k
r
i
z
=
∂
f
A
y
∂
i
A
y
|
i
A
y
=
i
z
=
0
x
=
y
=
0
=
2
μ
0
A
r
N
r
H
c
h
m
×
(
∑
(
s
z
0
-
z
)
s
r
0
4
A
r
(
s
z
0
-
z
)
+
A
z
s
r
0
)
s
r
0
2
(
4
A
r
s
r
0
+
A
z
s
z
0
+
z
)
h
h
h
h
h
h
×
(
h
m
μ
r
A
m
+
(
s
z
0
+
z
)
s
r
0
4
A
r
(
s
z
0
+
z
)
+
A
z
s
r
0
h
h
h
h
h
h
h
h
h
h
+
(
s
z
0
-
z
)
s
r
0
4
A
r
(
s
z
0
-
z
)
+
A
z
s
r
0
)
∑
)
-
1
,
k
r
s
z
=
∂
f
A
y
∂
y
|
i
A
y
=
i
z
=
0
x
=
y
=
0
=
-
2
μ
0
A
r
(
H
c
h
m
)
2
×
(
∑
+
s
r
0
(
s
z
0
-
z
)
4
A
r
(
s
z
0
-
z
)
+
A
z
s
r
0
)
2
s
r
0
3
(
4
A
r
s
r
0
+
A
z
s
z
0
+
z
)
2
h
h
h
h
h
×
(
h
m
μ
r
A
m
+
s
r
0
(
s
z
0
+
z
)
4
A
r
(
s
z
0
+
z
)
+
A
z
s
r
0
h
h
h
h
h
h
h
h
h
+
s
r
0
(
s
z
0
-
z
)
4
A
r
(
s
z
0
-
z
)
+
A
z
s
r
0
)
2
∑
)
-
1
.


#### 3.4.2. The Current Stiffness and Displacement Stiffness of the RMB Unit Affected by the Coil Current in AMB Unit

From Figure5(d), the control flux generated by the coil current in AMB unit will flow through radial air gaps of the RMB unit, and the current stiffness and displacement stiffness of the RMB unit could be affected and given by
(23)
k
r
i
i
z
=
∂
f
A
y
∂
i
A
y
|
i
A
y
=
i
B
y
=
0
x
=
y
=
z
=
0
=
2
μ
0
A
r
N
r
H
c
h
m
×
(
(
4
A
r
s
r
0
+
s
r
0
2
A
z
s
z
0
)
h
h
h
h
h
×
(
h
m
μ
r
A
m
+
2
s
r
0
s
z
0
4
A
r
s
z
0
+
A
z
s
r
0
)
(
4
A
r
s
r
0
+
s
r
0
2
A
z
s
z
0
)
)
-
1
-
μ
0
A
r
N
r
N
z
i
z
h
m
×
(
s
r
0
(
h
m
s
r
0
μ
0
A
m
s
r
0
+
2
A
r
h
m
+
2
s
z
0
A
z
)
h
h
h
h
h
×
(
μ
r
A
m
s
r
0
+
2
A
r
h
m
)
2
s
z
0
A
z
)
-
1
,
k
r
s
i
z
=
∂
f
A
y
∂
y
|
i
A
y
=
i
B
y
=
0
x
=
y
=
z
=
0
=
-
2
μ
0
A
r
s
r
0
[
∫
H
c
h
m
h
h
h
h
h
h
h
h
h
h
×
(
(
4
A
r
+
s
r
0
A
z
s
z
0
)
h
h
h
h
h
h
h
h
h
h
h
h
h
h
×
(
h
m
μ
r
A
m
+
2
s
r
0
s
z
0
4
A
r
s
z
0
+
A
z
s
r
0
)
)
-
1
h
h
h
h
h
h
h
h
h
h
h
-
N
z
i
z
(
h
m
s
r
0
)
hhhhhhhhhh
×
(
(
μ
r
A
m
s
r
0
+
2
A
r
h
m
)
2
s
z
0
A
z
)
-
1
2
s
r
0
(
h
m
s
r
0
μ
0
A
m
s
r
0
+
2
A
r
h
m
+
2
s
z
0
A
z
)
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
×
(
μ
r
A
m
s
r
0
+
2
A
r
h
m
)
2
s
z
0
A
z
)
-
1
∫
]
2
.


#### 3.4.3. The Current Stiffness and Displacement Stiffness of the AMB Unit Affected by the Radial Air Gap

The alternation of the radial air gap will affect the current stiffness and the displacement stiffness of the AMB unit and will be calculated by(24)
k
z
i
y
=
∂
f
z
∂
i
z
|
i
A
y
=
i
B
y
=
i
z
=
0
x
=
z
=
0
=
2
μ
0
N
z
H
c
h
m
(
A
r
s
r
0
-
y
+
A
r
s
r
0
+
y
+
2
A
r
s
r
0
+
A
z
s
z
0
)
×
(
s
z
0
(
∫
h
m
μ
r
A
m
hhhhhhhhh
+
2
2
A
r
s
r
0
/
(
s
r
0
2
-
y
2
)
+
2
A
r
/
s
r
0
+
A
z
/
s
z
0
∫
)
h
h
h
h
h
×
(
∫
1
h
m
/
μ
r
A
m
+
2
A
r
s
r
0
/
(
s
r
0
2
-
y
2
)
+
A
r
/
s
r
0
h
h
h
h
h
h
h
h
h
+
2
s
z
0
A
z
∫
)
)
-
1
,
k
z
s
y
=
∂
f
z
∂
z
|
i
A
y
=
i
B
y
=
i
z
=
0
x
=
z
=
0
=
2
μ
0
A
z
(
H
c
h
m
)
2
(
(
h
m
/
μ
r
A
m
)
+
(
2
s
r
0
s
z
0
/
(
4
A
r
s
z
0
+
A
z
s
r
0
)
)
)
2
×
A
z
-
s
z
0
(
(
4
A
r
/
s
r
0
)
+
(
A
z
/
s
z
0
)
)
s
z
0
4
(
(
4
A
r
/
s
r
0
)
+
(
A
z
/
s
z
0
)
)
3
.


## 3.1. The Equivalent Magnetic Circuits of the Combined 5-DOFs HMB

According to construction of the combined 5-DOFs HMB, the magnetic resistance of the iron path and flux leakage is neglected, the permeability of the iron is assumed to be infinite, and their equivalent bias and control magnetic circuits can be gained and are shown in Figure5.Equivalent bias circuit and control circuits of the combined 5-DOFs HMB. (a) The equivalent bias circuit. (b) The equivalent control circuit of RMB unitA. (c) The equivalent control circuit of RMB unit B. (d) The equivalent control circuit of AMB unit.
(a)
(b)
(c)
(d)


## 3.2. The Model of Combined 5-DOFs HMB for Translation inx-Axis (or y-Axis) and z-Axis Direction

Figure3 describes the operating principle of this combined 5-DOFs HMB in y-axis and z-axis. The rotor will depart from its center position when subject disturbance force in negative y-axis direction. The control flux on the upper side of the RMB units A and B will add to the bias flux in upper radial air gap and subtract from the bias flux in the lower radial air gap. And the control flux of the AMB unit should also be calculated. A net restoring force on the rotor will be produced in positive y-axis direction.According to the structure of the combined 5-DOFs HMB (Figure4) and its equivalent magnetic circuits (Figure 5), the net force of the RMB unit A in the y-axis direction (N
r
i
A
x
=
0, where N
r is the number of winding turns and i
A
x is the control current of x-axis coil currents), f
y, is
(1)
f
A
y
=
(
ϕ
y
+
A
p
+
ϕ
y
+
A
i
-
ϕ
y
+
A
z
)
2
2
μ
0
A
r
-
(
ϕ
y
-
A
p
-
ϕ
y
-
A
i
-
ϕ
y
-
A
z
)
2
2
μ
0
A
r
,
where the constant μ
0 is the permeability of free space and is equal to 4π × 10−7 H/m; A
r is the area of the pole face of RMB unit; ϕ
y
+
A
p and ϕ
y
-
A
p are the bias flux in radial air gaps of the RMB unit A; ϕ
y
+
A
i and ϕ
y
-
A
i are the control flux in the radial air gaps of the RMB unit A produced by the y-axis coil currents; ϕ
y
+
A
z and ϕ
y
-
A
z are the control flux in the radial air gaps of the RMB unit A produced by the z-axis coil current of the AMB unit.The bias flux in Figure5 can be calculated, respectively, as follows:
(2)
ϕ
y
+
A
p
=
F
m
R
sum
R
A
R
z
1
R
A
+
R
z
1
1
R
y
+
A
,
ϕ
y
-
A
p
=
F
m
R
sum
R
A
R
z
1
R
A
+
R
z
1
1
R
y
-
A
,
ϕ
z
1
p
=
F
m
R
sum
R
A
R
z
1
R
A
+
R
z
1
1
R
z
1
,
ϕ
z
2
p
=
F
m
R
sum
R
B
R
z
2
R
B
+
R
z
2
1
R
z
2
,
where
(3)
R
sum
=
R
m
+
R
A
R
z
1
R
A
+
R
z
1
+
R
B
R
z
2
R
B
+
R
z
2
,
R
A
=
1
(
1
/
R
x
+
A
)
+
(
1
/
R
x
-
A
)
+
(
1
/
R
y
+
A
)
+
(
1
/
R
y
-
A
)
,
R
B
=
1
(
1
/
R
x
+
B
)
+
(
1
/
R
x
-
B
)
+
(
1
/
R
y
+
B
)
+
(
1
/
R
y
-
B
)
.As shown in Figures3 and 5, R
sum is the total reluctance of the bias circuit; R
x
+
A and R
x
-
A are the reluctance of the radial air gaps in x-axis for the RMB unit A; R
y
+
A and R
y
-
A are the reluctance of the radial air gaps in y-axis for the RMB unit A; R
x
+
B and R
x
-
B are the reluctance of the radial air gaps in x-axis for the RMB unit B; R
y
+
B and R
y
-
B are the reluctance of the radial air gaps in y-axis for the RMB unit B; R
z
1 and R
z
2 are the reluctance of the axial air gaps in z-axis for the AMB unit. F
m
=
H
c
h
m, where H
c is the coercive force of the permanent magnet; h
m is PMR thickness; R
m is the magnetic resistance of PMR and it can be given as
(4)
R
m
=
h
m
μ
0
μ
r
A
m
,
where μ
r is the relative permeability of permanent magnet and A
m is the middle-section area of the PMR.The reluctances of radial air gap are written by(5)
R
x
+
A
=
s
r
0
+
x
μ
0
A
r
,
R
x
-
A
=
s
r
0
-
x
μ
0
A
r
,
R
x
+
B
=
s
r
0
+
x
μ
0
A
r
,
R
x
-
B
=
s
r
0
-
x
μ
0
A
r
,
R
y
+
A
=
s
r
0
+
y
μ
0
A
r
,
R
y
-
A
=
s
r
0
-
y
μ
0
A
r
,
R
y
+
B
=
s
r
0
+
y
μ
0
A
r
,
R
y
-
B
=
s
r
0
-
y
μ
0
A
r
,
where s
r
0 is the nominal air gap of the RMB unit; y is the rotor position as deviation from the operating point in y-axis direction; x is the deviation from the operating point in x-axis direction.The reluctances of axial air gap are written by(6)
R
z
1
=
s
z
0
+
z
μ
0
A
z
,
R
z
2
=
s
z
0
-
z
μ
0
A
z
,
where s
z
0 is the nominal air gap of the AMB unit; z is the rotor position as deviation from the operating point in z-axis direction; A
z is the area of the pole face for the AMB unit.The control flux in the radial air gaps ofy-axis, ϕ
y
+
A
i and ϕ
y
-
A
i, generated by y-axis coil currents of RMB unit A (N
r
i
A
x
=
0, where N
r is the number of winding turns and i
A
x is the control current of x-axis coil currents) can be calculated as
(7)
ϕ
y
+
A
i
=
N
r
i
A
y
R
y
+
A
+
R
x
-
A
/
/
R
x
+
A
/
/
R
y
-
A
+
N
r
i
A
y
R
y
-
A
+
R
x
-
A
/
/
R
x
+
A
/
/
R
y
+
A
R
x
-
A
/
/
R
x
+
A
/
/
R
y
+
A
R
y
+
A
,
ϕ
y
-
A
i
=
N
r
i
A
y
R
y
-
A
+
R
x
-
A
/
/
R
x
+
A
/
/
R
y
+
A
+
N
r
i
A
y
R
y
+
A
+
R
x
-
A
/
/
R
x
+
A
/
/
R
y
-
A
R
x
-
A
/
/
R
x
+
A
/
/
R
y
-
A
R
y
-
A
,
where i
A
y is the control current of y-axis coil currents for the same displacement of RMB unit A.The control flux in the radial air gaps,ϕ
y
+
A
z and ϕ
y
-
A
z, which are produced by the z-axis coil currents of the AMB unit can be calculated as
(8)
ϕ
y
+
A
z
=
ϕ
z
i
[
R
m
/
/
(
R
A
+
R
B
)
]
R
A
R
A
+
R
B
1
R
y
+
A
,
ϕ
y
-
A
z
=
ϕ
z
i
[
R
m
/
/
(
R
A
+
R
B
)
]
R
A
R
A
+
R
B
1
R
y
-
A
,
where the control flux in the axial air gap is calculated by
(9)
ϕ
z
i
=
N
z
i
z
R
z
1
+
R
z
2
+
R
m
/
/
(
R
A
+
R
B
)
,
where N
z is the number of winding turns of the AMB and i
z is the control current of z-axis coil current. And ϕ
z
i
1
=
ϕ
z
i
2.Substituting (2), (7), and (8) into the right side of (1) results in the net force on y-axis direction. Using the Taylor series expansion for small values of y and i
y, we can get the following attractive force with linear terms in y-axis direction:
(10)
f
A
y
(
y
,
i
y
)
≅
f
A
y
|
x
=
y
=
z
=
0
i
A
y
=
i
z
=
0
+
k
r
i
i
y
+
k
r
s
y
.Since the net force is zero at the center position in radial and axial directions, (10) can be written by
(11)
f
A
y
(
y
,
i
A
y
)
=
k
r
i
i
y
+
k
r
s
y
,
where k
r
i is the current stiffness of the RMB unit and k
r
s is the negative position stiffness of the RMB unit caused by the increasing attractive force as the air gap is reduced. This must be overcome for stable suspension since the homogeneous solution to magnetic bearing force (MBF) is in the form of hyperbolic functions indicating that y grows with time. They can be calculated as follows:
(12)
k
r
i
=
∂
f
A
y
∂
i
A
y
|
i
A
y
=
i
z
=
0
x
=
y
=
z
=
0
=
2
μ
0
A
r
N
r
H
c
h
m
×
(
s
r
0
2
(
4
A
r
s
r
0
+
A
z
s
z
0
)
h
h
h
h
h
×
(
h
m
μ
r
A
m
+
2
s
r
0
s
z
0
4
A
r
s
z
0
+
A
z
s
r
0
)
)
-
1
,
k
r
s
=
∂
f
A
y
∂
y
|
i
A
y
=
i
z
=
0
x
=
y
=
z
=
0
=
-
2
μ
0
A
r
(
H
c
h
m
)
2
×
(
s
r
0
3
(
4
A
r
s
r
0
+
A
z
s
z
0
)
2
h
h
h
h
h
h
×
(
h
m
μ
r
A
m
+
2
s
r
0
s
z
0
4
A
r
s
z
0
+
A
z
s
r
0
)
2
)
-
1
.According to the structure of the combined 5-DOFs HMB (Figure4) and its equivalent magnetic circuits (Figure 5), the net force of the AMB unit in the z-axis direction is
(13)
f
z
=
(
ϕ
z
1
p
+
ϕ
z
i
)
2
2
μ
0
A
z
-
(
ϕ
z
2
p
-
ϕ
z
i
)
2
2
μ
0
A
z
,
where ϕ
z
1
p and ϕ
z
2
p are the bias flux in axial air gaps of the AMB unit and ϕ
z
i is the control flux in the axial air gaps of the AMB unit by the z-axis coil current.Substituting (2) and (9) for the right side of (13) results in the net force on z-axis direction. Using the Taylor series expansion for small values of z and i
z, we can get the following attractive force with linear terms in z-axis direction:
(14)
f
z
(
z
,
i
z
)
≅
f
z
|
x
=
y
=
z
=
0
i
A
y
=
i
z
=
0
+
k
z
i
i
z
+
k
z
s
z
.Since the net force is zero at the center, (14) can be written by
(15)
f
z
(
z
,
i
z
)
=
k
z
i
i
z
+
k
z
s
z
,
where k
z
i is the current stiffness of the AMB unit and k
z
s is the negative position stiffness of the AMB unit caused by the increasing attractive force as the air gap is reduced. This must be overcome for stable suspension since the homogeneous solution to MBF is in the form of hyperbolic functions indicating that z grows with time. They can be calculated as follows:
(16)
k
z
i
=
∂
f
z
∂
i
z
|
i
A
y
=
i
B
y
=
i
z
=
0
x
=
y
=
z
=
0
=
2
μ
0
N
z
H
c
h
m
×
(
(
4
A
r
s
z
0
s
r
0
+
A
z
)
h
h
h
h
h
×
(
h
m
μ
r
A
m
+
2
s
r
0
s
z
0
4
A
r
s
z
0
+
A
z
s
r
0
)
h
h
h
h
h
×
(
h
m
s
r
0
μ
r
A
m
s
r
0
+
2
A
r
h
m
+
2
s
z
0
A
z
)
)
-
1
,
k
z
s
=
∂
f
z
∂
z
|
i
A
y
=
i
B
y
=
i
z
=
0
x
=
y
=
z
=
0
=
2
μ
0
A
z
(
H
c
h
m
)
2
(
(
h
m
/
μ
r
A
m
)
+
(
2
s
r
0
s
z
0
/
(
4
A
r
s
z
0
+
A
z
s
r
0
)
)
)
2
×
A
z
s
r
0
3
-
s
r
0
2
(
4
A
r
s
z
0
+
A
z
s
r
0
)
s
z
0
(
4
A
r
s
z
0
+
A
z
s
r
0
)
3
.


## 3.3. The Model of Combined 5-DOFs HMB for Tilting

The RMB unitsA and B of the combined 5-DOFs HMB also control 2-DOFs around x-axis and y-axis which are orthogonal to the spin axis (z-axis) when the rotor is subject to disturbance torque around x-axis and y-axis. The RMB units A and B create tilting stiffness to restrain the rotor when subjected to the disturbed torques as shown in Figure 4. According to the structure of the combined 5-DOFs HMB (Figure 3) and its equivalent magnetic circuits (Figure 5), the torque, M
x, acting on the rotor can be given by (17) (N
z
i
z
=
0 and N
r
i
A
x
=
N
r
i
B
x
=
0). Consider
(17)
M
x
=
l
cos
⁡
α
[
(
ϕ
y
+
A
p
+
ϕ
y
+
A
i
)
2
-
(
ϕ
y
-
A
p
-
ϕ
y
-
A
i
)
2
2
μ
0
A
r
hhhhhhhh
+
(
ϕ
y
-
B
p
+
ϕ
y
-
B
i
)
2
-
(
ϕ
y
+
B
p
-
ϕ
y
+
B
i
)
2
2
μ
0
A
r
]
,
where l is the distance from center of the magnetic pole of the RMB units A and B to the center of rotor as shown in Figure 4. The l
*
cos
⁡
α
≅
l for small value of the tilting angle, α. Equation (17) can be simplified as
(18)
M
x
≅
l
[
(
ϕ
y
+
A
p
+
ϕ
y
+
A
i
)
2
-
(
ϕ
y
-
A
p
-
ϕ
y
-
A
i
)
2
μ
0
A
r
]
.The change in radial air gap isl
*
tan
α
≅
l
*
α for small value of the tilting angle, α. The radial air gap resistance can be written as
(19)
R
x
+
A
=
s
r
0
μ
0
A
r
,
R
x
-
A
=
s
r
0
μ
0
A
r
,
R
x
+
B
=
s
r
0
μ
0
A
r
,
R
x
-
B
=
s
r
0
μ
0
A
r
,
R
y
+
A
=
R
y
-
B
≅
s
r
0
+
l
α
μ
0
A
r
,
R
y
-
A
=
R
y
+
B
≅
s
r
0
-
l
α
μ
0
A
r
.Using the Taylor series expansion for small values ofα and i
A
y we can get the following attractive torque with linear terms:
(20)
M
x
(
i
A
y
,
α
)
≅
k
α
i
i
A
y
+
k
α
s
α
,
where k
α
i is the tilting current stiffness and k
α
s is the negative tilting displacement stiffness. They can be calculated as follows:
(21)
k
α
i
=
∂
M
x
∂
i
A
y
|
i
A
y
=
0
α
=
0
=
4
μ
0
A
r
N
r
H
c
h
m
l
×
(
4
A
r
s
z
0
s
r
0
+
A
z
s
r
0
2
s
z
0
h
h
h
h
h
h
×
(
h
m
μ
r
A
m
+
2
s
r
0
s
z
0
4
A
r
s
z
0
+
A
z
s
r
0
)
4
A
r
s
z
0
s
r
0
+
A
z
s
r
0
2
s
z
0
)
-
1
,
k
α
s
=
∂
M
x
∂
α
|
i
A
y
=
0
α
=
0
=
-
4
μ
0
A
r
(
H
c
h
m
)
2
l
2
×
(
s
r
0
(
4
A
r
s
z
0
+
A
z
s
r
0
s
z
0
)
2
h
h
h
h
h
×
(
h
m
μ
r
A
m
+
2
s
r
0
s
z
0
4
A
r
s
z
0
+
A
z
s
r
0
)
2
)
-
1
.


## 3.4. The Coupling Model of the Combined 5-DOFs HMB

### 3.4.1. The Current Stiffness and Position of the RMB Unit Affected by the Axial Air Gap

The current stiffnessk
r
i and position stiffness k
r
s of the RMB unit in (12) do not consider the effect of the rotor in z-axis direction, namely, the axial air gap. However, the bias flux generated from the PMR flows across the radial air gaps and the axial air gaps (Figure 5(a)); the alternation of the axial air gap will affect the current stiffness and the position stiffness of the RMB and will be calculated by
(22)
k
r
i
z
=
∂
f
A
y
∂
i
A
y
|
i
A
y
=
i
z
=
0
x
=
y
=
0
=
2
μ
0
A
r
N
r
H
c
h
m
×
(
∑
(
s
z
0
-
z
)
s
r
0
4
A
r
(
s
z
0
-
z
)
+
A
z
s
r
0
)
s
r
0
2
(
4
A
r
s
r
0
+
A
z
s
z
0
+
z
)
h
h
h
h
h
h
×
(
h
m
μ
r
A
m
+
(
s
z
0
+
z
)
s
r
0
4
A
r
(
s
z
0
+
z
)
+
A
z
s
r
0
h
h
h
h
h
h
h
h
h
h
+
(
s
z
0
-
z
)
s
r
0
4
A
r
(
s
z
0
-
z
)
+
A
z
s
r
0
)
∑
)
-
1
,
k
r
s
z
=
∂
f
A
y
∂
y
|
i
A
y
=
i
z
=
0
x
=
y
=
0
=
-
2
μ
0
A
r
(
H
c
h
m
)
2
×
(
∑
+
s
r
0
(
s
z
0
-
z
)
4
A
r
(
s
z
0
-
z
)
+
A
z
s
r
0
)
2
s
r
0
3
(
4
A
r
s
r
0
+
A
z
s
z
0
+
z
)
2
h
h
h
h
h
×
(
h
m
μ
r
A
m
+
s
r
0
(
s
z
0
+
z
)
4
A
r
(
s
z
0
+
z
)
+
A
z
s
r
0
h
h
h
h
h
h
h
h
h
+
s
r
0
(
s
z
0
-
z
)
4
A
r
(
s
z
0
-
z
)
+
A
z
s
r
0
)
2
∑
)
-
1
.


### 3.4.2. The Current Stiffness and Displacement Stiffness of the RMB Unit Affected by the Coil Current in AMB Unit

From Figure5(d), the control flux generated by the coil current in AMB unit will flow through radial air gaps of the RMB unit, and the current stiffness and displacement stiffness of the RMB unit could be affected and given by
(23)
k
r
i
i
z
=
∂
f
A
y
∂
i
A
y
|
i
A
y
=
i
B
y
=
0
x
=
y
=
z
=
0
=
2
μ
0
A
r
N
r
H
c
h
m
×
(
(
4
A
r
s
r
0
+
s
r
0
2
A
z
s
z
0
)
h
h
h
h
h
×
(
h
m
μ
r
A
m
+
2
s
r
0
s
z
0
4
A
r
s
z
0
+
A
z
s
r
0
)
(
4
A
r
s
r
0
+
s
r
0
2
A
z
s
z
0
)
)
-
1
-
μ
0
A
r
N
r
N
z
i
z
h
m
×
(
s
r
0
(
h
m
s
r
0
μ
0
A
m
s
r
0
+
2
A
r
h
m
+
2
s
z
0
A
z
)
h
h
h
h
h
×
(
μ
r
A
m
s
r
0
+
2
A
r
h
m
)
2
s
z
0
A
z
)
-
1
,
k
r
s
i
z
=
∂
f
A
y
∂
y
|
i
A
y
=
i
B
y
=
0
x
=
y
=
z
=
0
=
-
2
μ
0
A
r
s
r
0
[
∫
H
c
h
m
h
h
h
h
h
h
h
h
h
h
×
(
(
4
A
r
+
s
r
0
A
z
s
z
0
)
h
h
h
h
h
h
h
h
h
h
h
h
h
h
×
(
h
m
μ
r
A
m
+
2
s
r
0
s
z
0
4
A
r
s
z
0
+
A
z
s
r
0
)
)
-
1
h
h
h
h
h
h
h
h
h
h
h
-
N
z
i
z
(
h
m
s
r
0
)
hhhhhhhhhh
×
(
(
μ
r
A
m
s
r
0
+
2
A
r
h
m
)
2
s
z
0
A
z
)
-
1
2
s
r
0
(
h
m
s
r
0
μ
0
A
m
s
r
0
+
2
A
r
h
m
+
2
s
z
0
A
z
)
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
×
(
μ
r
A
m
s
r
0
+
2
A
r
h
m
)
2
s
z
0
A
z
)
-
1
∫
]
2
.


### 3.4.3. The Current Stiffness and Displacement Stiffness of the AMB Unit Affected by the Radial Air Gap

The alternation of the radial air gap will affect the current stiffness and the displacement stiffness of the AMB unit and will be calculated by(24)
k
z
i
y
=
∂
f
z
∂
i
z
|
i
A
y
=
i
B
y
=
i
z
=
0
x
=
z
=
0
=
2
μ
0
N
z
H
c
h
m
(
A
r
s
r
0
-
y
+
A
r
s
r
0
+
y
+
2
A
r
s
r
0
+
A
z
s
z
0
)
×
(
s
z
0
(
∫
h
m
μ
r
A
m
hhhhhhhhh
+
2
2
A
r
s
r
0
/
(
s
r
0
2
-
y
2
)
+
2
A
r
/
s
r
0
+
A
z
/
s
z
0
∫
)
h
h
h
h
h
×
(
∫
1
h
m
/
μ
r
A
m
+
2
A
r
s
r
0
/
(
s
r
0
2
-
y
2
)
+
A
r
/
s
r
0
h
h
h
h
h
h
h
h
h
+
2
s
z
0
A
z
∫
)
)
-
1
,
k
z
s
y
=
∂
f
z
∂
z
|
i
A
y
=
i
B
y
=
i
z
=
0
x
=
z
=
0
=
2
μ
0
A
z
(
H
c
h
m
)
2
(
(
h
m
/
μ
r
A
m
)
+
(
2
s
r
0
s
z
0
/
(
4
A
r
s
z
0
+
A
z
s
r
0
)
)
)
2
×
A
z
-
s
z
0
(
(
4
A
r
/
s
r
0
)
+
(
A
z
/
s
z
0
)
)
s
z
0
4
(
(
4
A
r
/
s
r
0
)
+
(
A
z
/
s
z
0
)
)
3
.


## 3.4.1. The Current Stiffness and Position of the RMB Unit Affected by the Axial Air Gap

The current stiffnessk
r
i and position stiffness k
r
s of the RMB unit in (12) do not consider the effect of the rotor in z-axis direction, namely, the axial air gap. However, the bias flux generated from the PMR flows across the radial air gaps and the axial air gaps (Figure 5(a)); the alternation of the axial air gap will affect the current stiffness and the position stiffness of the RMB and will be calculated by
(22)
k
r
i
z
=
∂
f
A
y
∂
i
A
y
|
i
A
y
=
i
z
=
0
x
=
y
=
0
=
2
μ
0
A
r
N
r
H
c
h
m
×
(
∑
(
s
z
0
-
z
)
s
r
0
4
A
r
(
s
z
0
-
z
)
+
A
z
s
r
0
)
s
r
0
2
(
4
A
r
s
r
0
+
A
z
s
z
0
+
z
)
h
h
h
h
h
h
×
(
h
m
μ
r
A
m
+
(
s
z
0
+
z
)
s
r
0
4
A
r
(
s
z
0
+
z
)
+
A
z
s
r
0
h
h
h
h
h
h
h
h
h
h
+
(
s
z
0
-
z
)
s
r
0
4
A
r
(
s
z
0
-
z
)
+
A
z
s
r
0
)
∑
)
-
1
,
k
r
s
z
=
∂
f
A
y
∂
y
|
i
A
y
=
i
z
=
0
x
=
y
=
0
=
-
2
μ
0
A
r
(
H
c
h
m
)
2
×
(
∑
+
s
r
0
(
s
z
0
-
z
)
4
A
r
(
s
z
0
-
z
)
+
A
z
s
r
0
)
2
s
r
0
3
(
4
A
r
s
r
0
+
A
z
s
z
0
+
z
)
2
h
h
h
h
h
×
(
h
m
μ
r
A
m
+
s
r
0
(
s
z
0
+
z
)
4
A
r
(
s
z
0
+
z
)
+
A
z
s
r
0
h
h
h
h
h
h
h
h
h
+
s
r
0
(
s
z
0
-
z
)
4
A
r
(
s
z
0
-
z
)
+
A
z
s
r
0
)
2
∑
)
-
1
.


## 3.4.2. The Current Stiffness and Displacement Stiffness of the RMB Unit Affected by the Coil Current in AMB Unit

From Figure5(d), the control flux generated by the coil current in AMB unit will flow through radial air gaps of the RMB unit, and the current stiffness and displacement stiffness of the RMB unit could be affected and given by
(23)
k
r
i
i
z
=
∂
f
A
y
∂
i
A
y
|
i
A
y
=
i
B
y
=
0
x
=
y
=
z
=
0
=
2
μ
0
A
r
N
r
H
c
h
m
×
(
(
4
A
r
s
r
0
+
s
r
0
2
A
z
s
z
0
)
h
h
h
h
h
×
(
h
m
μ
r
A
m
+
2
s
r
0
s
z
0
4
A
r
s
z
0
+
A
z
s
r
0
)
(
4
A
r
s
r
0
+
s
r
0
2
A
z
s
z
0
)
)
-
1
-
μ
0
A
r
N
r
N
z
i
z
h
m
×
(
s
r
0
(
h
m
s
r
0
μ
0
A
m
s
r
0
+
2
A
r
h
m
+
2
s
z
0
A
z
)
h
h
h
h
h
×
(
μ
r
A
m
s
r
0
+
2
A
r
h
m
)
2
s
z
0
A
z
)
-
1
,
k
r
s
i
z
=
∂
f
A
y
∂
y
|
i
A
y
=
i
B
y
=
0
x
=
y
=
z
=
0
=
-
2
μ
0
A
r
s
r
0
[
∫
H
c
h
m
h
h
h
h
h
h
h
h
h
h
×
(
(
4
A
r
+
s
r
0
A
z
s
z
0
)
h
h
h
h
h
h
h
h
h
h
h
h
h
h
×
(
h
m
μ
r
A
m
+
2
s
r
0
s
z
0
4
A
r
s
z
0
+
A
z
s
r
0
)
)
-
1
h
h
h
h
h
h
h
h
h
h
h
-
N
z
i
z
(
h
m
s
r
0
)
hhhhhhhhhh
×
(
(
μ
r
A
m
s
r
0
+
2
A
r
h
m
)
2
s
z
0
A
z
)
-
1
2
s
r
0
(
h
m
s
r
0
μ
0
A
m
s
r
0
+
2
A
r
h
m
+
2
s
z
0
A
z
)
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
×
(
μ
r
A
m
s
r
0
+
2
A
r
h
m
)
2
s
z
0
A
z
)
-
1
∫
]
2
.


## 3.4.3. The Current Stiffness and Displacement Stiffness of the AMB Unit Affected by the Radial Air Gap

The alternation of the radial air gap will affect the current stiffness and the displacement stiffness of the AMB unit and will be calculated by(24)
k
z
i
y
=
∂
f
z
∂
i
z
|
i
A
y
=
i
B
y
=
i
z
=
0
x
=
z
=
0
=
2
μ
0
N
z
H
c
h
m
(
A
r
s
r
0
-
y
+
A
r
s
r
0
+
y
+
2
A
r
s
r
0
+
A
z
s
z
0
)
×
(
s
z
0
(
∫
h
m
μ
r
A
m
hhhhhhhhh
+
2
2
A
r
s
r
0
/
(
s
r
0
2
-
y
2
)
+
2
A
r
/
s
r
0
+
A
z
/
s
z
0
∫
)
h
h
h
h
h
×
(
∫
1
h
m
/
μ
r
A
m
+
2
A
r
s
r
0
/
(
s
r
0
2
-
y
2
)
+
A
r
/
s
r
0
h
h
h
h
h
h
h
h
h
+
2
s
z
0
A
z
∫
)
)
-
1
,
k
z
s
y
=
∂
f
z
∂
z
|
i
A
y
=
i
B
y
=
i
z
=
0
x
=
z
=
0
=
2
μ
0
A
z
(
H
c
h
m
)
2
(
(
h
m
/
μ
r
A
m
)
+
(
2
s
r
0
s
z
0
/
(
4
A
r
s
z
0
+
A
z
s
r
0
)
)
)
2
×
A
z
-
s
z
0
(
(
4
A
r
/
s
r
0
)
+
(
A
z
/
s
z
0
)
)
s
z
0
4
(
(
4
A
r
/
s
r
0
)
+
(
A
z
/
s
z
0
)
)
3
.


## 4. The Performance Analysis of the Combined 5-DOFs HMB Using Field-Circuit Method

### 4.1. The Combined 5-DOFs HMB Parameters and the 3D Finite Element Model

The combined 5-DOFs HMB parameters are given in Table1. The example aims to validate the design and analysis method based on linearized model and 3D finite element method (FEM).Table 1
Design parameters of the combined 5-DOFs HMB.
Parameters
Value
Inner diameter of the RMB unit rotor stacks (mm)
46
Outer diameter of the RMB unit rotor stacks (mm)
68
Length of the RMB unit rotor stacks (mm)
20
Inner diameter of the RMB unit stator (mm)
69
Outer diameter of the RMB unit stator (mm)
113
Length of the RMB unit stator (mm)
20
Radial pole shoe area (mm2)
734
The number of winding turns in coil of RMB unit
220
Inner diameter of the axial pole shoe (mm)
76
Outer diameter of the axial pole shoe (mm)
90
The number of winding turns in coil of AMB unit
250
The radial and axial air gap (mm)
0.5
Relative permeability of the permanent magnet ring
1.082
Outer diameter of the permanent magnet ring (mm)
126
Inner diameter of the permanent magnet ring (mm)
113
Length of the permanent magnet ring (mm)
5The linearized model is given by(25)
f
x
=
k
x
i
i
x
+
k
x
s
x
=
284.4
i
x
-
5.5195
×
1
0
5
x
,
f
y
=
k
y
i
i
y
+
k
y
s
y
=
284.4
i
y
-
5.5195
×
1
0
5
y
,
f
z
=
k
z
i
i
z
+
k
z
s
z
=
290.5
i
z
-
8.5007
×
1
0
5
z
,
M
x
=
k
α
i
i
α
+
k
α
s
α
=
39.25
i
α
-
5255.6
α
,
M
y
=
k
β
i
i
β
+
k
β
s
β
=
39.25
i
β
-
5255.6
β
,
where k
i
x and k
i
y are the current stiffness of the system in x-axis and y-axis, respectively; k
s
p is the displacement stiffness of the system in x-axis and y-axis; k
z is the displacement stiffness in z-axis; k
α and k
β are the tilting displacement stiffness around x-axis and y-axis; k
α
i and k
β
i are the tilting current stiffness around x-axis and y-axis. Moreover, the rotor can be controlled in five degrees-of-freedom.The finite element model using Ansoft Maxwell V13 is built and shown in Figure6(a) except the air. Its meshes are shown in Figures 6(b) and 6(c), and the total number of mesh is 286036.The 3D finite element model and mesh of the 5-DOFs HMB except the air. (a) The 3D finite element model. (b) The 3D finite element mesh. (c) Section view of the 3D finite element mesh.
(a)
(b)
(c)


### 4.2. The Performance Analysis for Radial and Axial Translation of the Combined 5-DOFs HMB

#### 4.2.1. The Analysis Results of the Force-Displacement and Displacement Stiffness of the RMB Unit (Using Linearized Method, Nonlinear Method, and FEM)

When the control currents in radial coil currents and axial coil current are 0 A, and the rotor position changes from −0.4 mm to 0.4 mm inx-axis (or y-axis), and hold in center position in z-axial direction, the relationship between net force and the rotor position in x-axis (or y-axis) direction is shown in Figure 7(a) using linearized method, nonlinear method, and FEM. It is obvious that the force is nonlinear. The calculated results are approximately equal for small values of the rotor in x-axis (or y-axis) and there is little difference between the nonlinear model and FEM. In addition, the displacement stiffnesses using the three calculated methods are shown in Figure 7(b). The displacement stiffness in rotor center is -
5.5195
×
10
5
N
/
m using linearized method and nonlinear method, and the displacement stiffness calculated by the FEM in rotor center is -
5.3815
×
10
5
N
/
m; the calculated error between them is 2.6%.Curves of the force-rotor position and displacement stiffness inx-axis (or y-axis) of the RMB unit (using linearized method, nonlinear method, and FEM). (a) Curves of the force-rotor position of the RMB unit. (b) Curves of the displacement stiffness of the RMB unit.
(a)
(b)


#### 4.2.2. The Analysis Results of the Force-Coil Current and Current Stiffness of the RMB Unit (Using Linearized Method, Nonlinear Method, and FEM)

When the rotor locates at radial and axial center and the coil current inx-axis ory-axis of the RMB unit changes from −1 A to 1 A (the coil current of the AMB unit is 0 A), the relationship between net force of the RMB unit and coil current is shown in Figure 8 using linearized method (see (11)), nonlinear method (see (1)), and FEM. It is obvious that the net force is linear, and the calculated results are approximately equal. In addition, the current stiffness of the RMB unit using the three calculated methods is 284.4 N/A (linearized method and nonlinear method) and 284.13 N/A (FEM); the calculated error is between linearized method and FEM 1.5%.Figure 8
Curves of the force-coil current inx-axis (or y-axis) of the RMB unit (using linearized method, nonlinear method, and FEM).


#### 4.2.3. The Analysis Results of the Force-Coil Current and Rotor Position inx-Axis (or y-Axis) of the RMB Unit (Using Nonlinear Method)

When the rotor inz-axis locates its center position and the coil current of the AMB unit is 0 A, the curves of the net force-coil current-rotor position of the RMB unit in x-axis (or y-axis) using nonlinear method are shown in Figure 9. It is obvious that it is approximately linearity for small displacement and coil current and nonlinear for large displacement and coil current. The displacement stiffness and current stiffness affected by the rotor position and the coil current in x-axis (or y-axis) of the RMB unit are shown in Figures 10 and 11, respectively. They are approximately linear for small displacement and coil current.Figure 9
Curves of the force-coil current and rotor position inx-axis (or y-axis) of the RMB unit (using nonlinear method).Figure 10
Curves of the displacement stiffness-rotor position-coil current and rotor position inx-axis (or y-axis) of the RMB unit (using nonlinear method).Figure 11
Curves of the current stiffness-coil current and rotor position inx-axis (or y-axis) of the RMB unit (using nonlinear method).


#### 4.2.4. The Analysis Results of the Force-Displacement and Displacement Stiffness of the AMB Unit (Using Linearized Method, Nonlinear Method, and FEM)

When the control currents in radial coil currents and axial coil current all are 0 A, and the rotor position changes from −0.4 mm to 0.4 mm inz-axis, and hold in center position in x-axial and y-axis direction, the relationship between net force and the rotor position in z-axis direction is shown in Figure 12(a) using linearized method (see (11)), nonlinear method (see (1)), and FEM. It is obvious that the force is also nonlinear. The calculated results are approximately equal for small values of the rotor in z-axis and there is little difference between the nonlinear model and FEM. In addition, the displacement stiffness using the three calculated methods is shown in Figure 12(b). The calculated results are approximately equal for small value of the displacement in z-axis between the linearized method and FEM. The displacement stiffness in rotor center is -
8.5007
×
10
5
N
/
m using linearized method and nonlinear method, and the displacement stiffness calculated by the FEM in rotor center is -
8.3732
×
10
5 N/m; the calculated error between them is 4.4%.The force-rotor position inz-axis of the AMB unit and its displacement stiffness (using linearized method, nonlinear method, and FEM). (a) The force-rotor position in z-axis of the AMB unit. (b) The displacement stiffness of the AMB unit.
(a)
(b)


#### 4.2.5. The Analysis Results of the Force-Coil Current and Current Stiffness of the AMB Unit (Using Linearized Method, Nonlinear Method, and FEM)

When the rotor locates at radial and axial center and the coil current inz-axis of the AMB unit changes from −1 A to 1 A (the coil currents of the RMB unit are 0 A), the relationship between net force of AMB unit and its coil current is shown in Figure 13 using linearized method (see (15)), nonlinear method (see (13)), and FEM. It is obvious that the net force is linear, and the calculated results are approximately equal. In addition, the current stiffness of the AMB unit using the three calculated methods is 290.5 N/A (linearized method and nonlinear method) and 280.3 N/A (FEM); the calculated error is 3.5% between linearized method and FEM.Figure 13
The force-coil current inz-axis of the AMB unit (using linearized method, nonlinear method, and FEM).


#### 4.2.6. The Analysis Results of the Force-Coil Current and Rotor Position inz-Axis of the AMB Unit (Using Nonlinear Method)

When the rotor locates its center position inx-axis and y-axis and the coil currents of the RMB units are 0 A, the curves of the net force-coil current-rotor position of the AMB unit in z-axis using nonlinear method are shown in Figure 14. It is obvious that the relationship of the net force-displacement in z-axis is approximately linear and the relationship of the net force-versus coil current is approximately linear for small displacement and nonlinear for large coil current. The displacement stiffness affected by the rotor position and the coil current in z-axis of the AMB unit is shown in Figure 15(a). The current stiffness affected by the rotor position and the coil current in z-axis of the AMB unit is shown in Figure 15(b). In addition, the current stiffness is only affected by the coil current of the AMB unit.Figure 14
The force-coil current and rotor position inz-axis of the AMB unit (using nonlinear method).The curves of the stiffness versus rotor position inz-axis and coil current of the AMB unit (using nonlinear method). (a) The displacement stiffness versus coil current and rotor position in z-axis of the AMB unit. (b) The current stiffness versus coil current and rotor position in z-axis of the AMB unit.
(a)
(b)


### 4.3. The Performance Analysis for Tilting of the Combined 5-DOFs HMB

#### 4.3.1. The Analysis Results of the Torque-Angular Position and Angular Displacement Stiffness of the Combined 5-DOFs HMB

The torque can be produced by the RMB unitsA and B. When the control currents in radial coil and axial coil all are 0 A and the rotor tilts within −5 mrad~5 mrad around x-axis or y-axis, the relationship between torque and angular position of rotor is shown in Figure 16(a) using linearized method, nonlinear method, and FEM. The calculated results are approximately equal for small values of the angular position of rotor and there is little difference between the nonlinear model and FEM. In addition, the tilting position stiffness is shown in Figure 16(b). The calculated results are approximately equal for small value of the angular position between the linearized method and FEM. The angular displacement stiffness in rotor center is −5255.6 Nm/rad.Curves of the torque-angular position of the combined 5-DOFs HMB and its tilting displacement stiffness (using linearized method, nonlinear method, and FEM). (a) Curves of the torque-angular position of the combined 5-DOFs HMB. (b) The tilting displacement stiffness.
(a)
(b)


#### 4.3.2.  The Analysis Results of the Torque-Coil Current and Tilting Current Stiffness of the Combined 5-DOFs HMB (Using Linearized Method, Nonlinear Method, and FEM)

When the rotor locates at radial and axial center and the coil currents of the RMB unitsA and B in x-axis or y-axis change from −1 A to 1 A (the coil current of AMB unit is 0 A), the relationship between torque and coil current of the RMB units is shown in Figure 17 using linearized method (see (15)), nonlinear method (see (13)) and FEM. It is obvious that the torque is linear, and the calculated results are approximately equal. In addition, the tilting current stiffness using the three calculated methods is 39.25 Nm/A (linearized method and nonlinear method) and 38.07 Nm/A (FEM); the calculated error between the linearized method and FEM is 3%.Figure 17
Curves of the torque-coil current in RMB unitsA and B of the combined 5-DOFs HMB (using linearized method, nonlinear method, and FEM).


#### 4.3.3. The Analysis Results of the Torque-Coil Current and Angular Position inx-Axis or y-Axis of the Combined 5-DOFs HMB (Using Linearized Method and Nonlinear Method)

When the rotor locates its center position inz-axis and the coil currents of the AMB unit are 0 A, the curves of the torque produced by the RMB units A and B versus their coil current and rotor angular position around x-axis or y-axis using nonlinear method are shown in Figure 18 using linearized method and nonlinear method. It is obvious that it is approximately linear for small values of the angular position and coil current. Moreover, the curves of the angular position stiffness versus angular position and coil currents in RMB units A and B are shown in Figure 19(a). And the curves of the tilting current stiffness versus angular position and coil currents in RMB units A and B are shown in Figure 19(b). They all are affected by the angular position and coil currents in RMB units A and B.Figure 18
Curves of the torque-coil current and angular position aroundx-axis or y-axis in the RMB units A and B of the combined 5-DOFs HMB (using linearized method and nonlinear method).Curves of the tilting displacement stiffness and the tilting current stiffness versus angular position aroundx-axis or y-axis and the coil current in the RMB units A and B of the combined 5-DOFs HMB (using linearized method and nonlinear method). (a) Tilting displacement stiffness versus angular position around x-axis or y-axis and the coil current. (b) Tilting current stiffness versus angular position around x-axis or y-axis and the coil current.
(a)
(b)


### 4.4. The Coupling Performance Analysis of the Combined 5-DOFs HMB

#### 4.4.1. The Coupling Performance Analysis for Translation of the Combined 5-DOFs HMB

Equation (22) shows that the displacement stiffness and current stiffness of the RMB unit will be affected by the rotor position in z-axis and the coil current of the AMB unit. The effect of the rotor position in z-axis on the displacement stiffness and current stiffness of the RMB units A and B is shown in Figures 20(a) and 20(b), respectively. They all are affected by the rotor position in z-axis and their change trend is opposite.The displacement stiffness and the current stiffness of the RMB unitsA and B affected by the rotor position of the AMB unit in z-axis. (a) The displacement stiffness of the RMB units A and B affected by the rotor position of the AMB unit in z-axis. (b) The current stiffness of the RMB units A and B affected by the rotor position of the AMB unit in z-axis.
(a)
(b)The effect of the coil current of the AMB unit on the displacement stiffness and current stiffness of the RMB unitsA and B is shown in Figures 21(a) and 21(b), respectively. They all are affected by the coil current of the AMB unit and their change trend is the same.The displacement and current stiffness of the RMB unitsA and B affected by the coil current of the AMB unit. (a) The displacement stiffness of the RMB units A and B affected by the coil current of the AMB unit. (b) The current stiffness of the RMB units A and B affected by the coil current of the AMB unit.
(a)
(b)The 3D curves of the displacement stiffness and current stiffness of the RMB unitsA and B affected by the rotor position in z-axis and the coil current of the AMB unit are shown in Figures 22(a) and 22(b).The displacement stiffness and the current stiffness of the RMB unitsA and B affected by the rotor position in z-axis and the coil current of the AMB unit. (a) The displacement stiffness of the RMB units A and B affected by the rotor position in z-axis and the coil current of the AMB unit. (b) The current stiffness of the RMB units A and B affected by the rotor position in z-axis and the coil current of the AMB unit.
(a)
(b)Figure23 shows that the displacement stiffness and current stiffness of the AMB unit will be affected by the rotor position in x-axis or y-axis of the RMB unit. The effect of the rotor position in x-axis or y-axis on the displacement stiffness and current stiffness of the AMB unit is shown in Figures 23(a) and 23(b), respectively. They all are affected by the rotor position in x-axis or y-axis.The displacement and the current stiffness of the AMB unit affected by the rotor position of the RMB unit inx-axis or y-axis. (a) The displacement stiffness of the AMB unit affected by the rotor position of the RMB unit. (b) The current stiffness of the AMB unit affected by the rotor position of the RMB unit.
(a)
(b)


#### 4.4.2. The Coupling Performance Analysis for Tilting of the Combined 5-DOFs HMB

The tilting displacement stiffness and tilting current stiffness of the RMB unitsA and B will be affected by the rotor position in z-axis and coil current of the AMB unit, and the influence curves of RMB units A and B are shown in Figures 24(a), 24(b), 25(a), and 25(b), respectively. It is obvious that the tilting displacement stiffness and tilting current stiffness of the RMB units A and B all are affected by the rotor position in z-axis and coil current of the AMB unit.Curves of the tilting displacement stiffness and tilting current stiffness of the RMB unitsA and B affected by the rotor position in z-axis of the AMB unit. (a) The tilting displacement stiffness of the RMB units A and B affected by the rotor position in z-axis. (b) The tilting current stiffness of the RMB units A and B affected by the rotor position in z-axis of the AMB unit.
(a)
(b)Curves of the tilting displacement stiffness and tilting current stiffness of the RMB unitsA and B affected by the rotor position in z-axis and coil current of the AMB unit. (a) Curves of the tilting displacement stiffness of the RMB units A and B affected by the rotor position in z-axis and coil current of the AMB unit. (b) Curves of the tilting current stiffness of the RMB units A and B affected by the rotor position in z-axis and coil current of the AMB unit.
(a)
(b)


## 4.1. The Combined 5-DOFs HMB Parameters and the 3D Finite Element Model

The combined 5-DOFs HMB parameters are given in Table1. The example aims to validate the design and analysis method based on linearized model and 3D finite element method (FEM).Table 1
Design parameters of the combined 5-DOFs HMB.
Parameters
Value
Inner diameter of the RMB unit rotor stacks (mm)
46
Outer diameter of the RMB unit rotor stacks (mm)
68
Length of the RMB unit rotor stacks (mm)
20
Inner diameter of the RMB unit stator (mm)
69
Outer diameter of the RMB unit stator (mm)
113
Length of the RMB unit stator (mm)
20
Radial pole shoe area (mm2)
734
The number of winding turns in coil of RMB unit
220
Inner diameter of the axial pole shoe (mm)
76
Outer diameter of the axial pole shoe (mm)
90
The number of winding turns in coil of AMB unit
250
The radial and axial air gap (mm)
0.5
Relative permeability of the permanent magnet ring
1.082
Outer diameter of the permanent magnet ring (mm)
126
Inner diameter of the permanent magnet ring (mm)
113
Length of the permanent magnet ring (mm)
5The linearized model is given by(25)
f
x
=
k
x
i
i
x
+
k
x
s
x
=
284.4
i
x
-
5.5195
×
1
0
5
x
,
f
y
=
k
y
i
i
y
+
k
y
s
y
=
284.4
i
y
-
5.5195
×
1
0
5
y
,
f
z
=
k
z
i
i
z
+
k
z
s
z
=
290.5
i
z
-
8.5007
×
1
0
5
z
,
M
x
=
k
α
i
i
α
+
k
α
s
α
=
39.25
i
α
-
5255.6
α
,
M
y
=
k
β
i
i
β
+
k
β
s
β
=
39.25
i
β
-
5255.6
β
,
where k
i
x and k
i
y are the current stiffness of the system in x-axis and y-axis, respectively; k
s
p is the displacement stiffness of the system in x-axis and y-axis; k
z is the displacement stiffness in z-axis; k
α and k
β are the tilting displacement stiffness around x-axis and y-axis; k
α
i and k
β
i are the tilting current stiffness around x-axis and y-axis. Moreover, the rotor can be controlled in five degrees-of-freedom.The finite element model using Ansoft Maxwell V13 is built and shown in Figure6(a) except the air. Its meshes are shown in Figures 6(b) and 6(c), and the total number of mesh is 286036.The 3D finite element model and mesh of the 5-DOFs HMB except the air. (a) The 3D finite element model. (b) The 3D finite element mesh. (c) Section view of the 3D finite element mesh.
(a)
(b)
(c)


## 4.2. The Performance Analysis for Radial and Axial Translation of the Combined 5-DOFs HMB

### 4.2.1. The Analysis Results of the Force-Displacement and Displacement Stiffness of the RMB Unit (Using Linearized Method, Nonlinear Method, and FEM)

When the control currents in radial coil currents and axial coil current are 0 A, and the rotor position changes from −0.4 mm to 0.4 mm inx-axis (or y-axis), and hold in center position in z-axial direction, the relationship between net force and the rotor position in x-axis (or y-axis) direction is shown in Figure 7(a) using linearized method, nonlinear method, and FEM. It is obvious that the force is nonlinear. The calculated results are approximately equal for small values of the rotor in x-axis (or y-axis) and there is little difference between the nonlinear model and FEM. In addition, the displacement stiffnesses using the three calculated methods are shown in Figure 7(b). The displacement stiffness in rotor center is -
5.5195
×
10
5
N
/
m using linearized method and nonlinear method, and the displacement stiffness calculated by the FEM in rotor center is -
5.3815
×
10
5
N
/
m; the calculated error between them is 2.6%.Curves of the force-rotor position and displacement stiffness inx-axis (or y-axis) of the RMB unit (using linearized method, nonlinear method, and FEM). (a) Curves of the force-rotor position of the RMB unit. (b) Curves of the displacement stiffness of the RMB unit.
(a)
(b)


### 4.2.2. The Analysis Results of the Force-Coil Current and Current Stiffness of the RMB Unit (Using Linearized Method, Nonlinear Method, and FEM)

When the rotor locates at radial and axial center and the coil current inx-axis ory-axis of the RMB unit changes from −1 A to 1 A (the coil current of the AMB unit is 0 A), the relationship between net force of the RMB unit and coil current is shown in Figure 8 using linearized method (see (11)), nonlinear method (see (1)), and FEM. It is obvious that the net force is linear, and the calculated results are approximately equal. In addition, the current stiffness of the RMB unit using the three calculated methods is 284.4 N/A (linearized method and nonlinear method) and 284.13 N/A (FEM); the calculated error is between linearized method and FEM 1.5%.Figure 8
Curves of the force-coil current inx-axis (or y-axis) of the RMB unit (using linearized method, nonlinear method, and FEM).


### 4.2.3. The Analysis Results of the Force-Coil Current and Rotor Position inx-Axis (or y-Axis) of the RMB Unit (Using Nonlinear Method)

When the rotor inz-axis locates its center position and the coil current of the AMB unit is 0 A, the curves of the net force-coil current-rotor position of the RMB unit in x-axis (or y-axis) using nonlinear method are shown in Figure 9. It is obvious that it is approximately linearity for small displacement and coil current and nonlinear for large displacement and coil current. The displacement stiffness and current stiffness affected by the rotor position and the coil current in x-axis (or y-axis) of the RMB unit are shown in Figures 10 and 11, respectively. They are approximately linear for small displacement and coil current.Figure 9
Curves of the force-coil current and rotor position inx-axis (or y-axis) of the RMB unit (using nonlinear method).Figure 10
Curves of the displacement stiffness-rotor position-coil current and rotor position inx-axis (or y-axis) of the RMB unit (using nonlinear method).Figure 11
Curves of the current stiffness-coil current and rotor position inx-axis (or y-axis) of the RMB unit (using nonlinear method).


### 4.2.4. The Analysis Results of the Force-Displacement and Displacement Stiffness of the AMB Unit (Using Linearized Method, Nonlinear Method, and FEM)

When the control currents in radial coil currents and axial coil current all are 0 A, and the rotor position changes from −0.4 mm to 0.4 mm inz-axis, and hold in center position in x-axial and y-axis direction, the relationship between net force and the rotor position in z-axis direction is shown in Figure 12(a) using linearized method (see (11)), nonlinear method (see (1)), and FEM. It is obvious that the force is also nonlinear. The calculated results are approximately equal for small values of the rotor in z-axis and there is little difference between the nonlinear model and FEM. In addition, the displacement stiffness using the three calculated methods is shown in Figure 12(b). The calculated results are approximately equal for small value of the displacement in z-axis between the linearized method and FEM. The displacement stiffness in rotor center is -
8.5007
×
10
5
N
/
m using linearized method and nonlinear method, and the displacement stiffness calculated by the FEM in rotor center is -
8.3732
×
10
5 N/m; the calculated error between them is 4.4%.The force-rotor position inz-axis of the AMB unit and its displacement stiffness (using linearized method, nonlinear method, and FEM). (a) The force-rotor position in z-axis of the AMB unit. (b) The displacement stiffness of the AMB unit.
(a)
(b)


### 4.2.5. The Analysis Results of the Force-Coil Current and Current Stiffness of the AMB Unit (Using Linearized Method, Nonlinear Method, and FEM)

When the rotor locates at radial and axial center and the coil current inz-axis of the AMB unit changes from −1 A to 1 A (the coil currents of the RMB unit are 0 A), the relationship between net force of AMB unit and its coil current is shown in Figure 13 using linearized method (see (15)), nonlinear method (see (13)), and FEM. It is obvious that the net force is linear, and the calculated results are approximately equal. In addition, the current stiffness of the AMB unit using the three calculated methods is 290.5 N/A (linearized method and nonlinear method) and 280.3 N/A (FEM); the calculated error is 3.5% between linearized method and FEM.Figure 13
The force-coil current inz-axis of the AMB unit (using linearized method, nonlinear method, and FEM).


### 4.2.6. The Analysis Results of the Force-Coil Current and Rotor Position inz-Axis of the AMB Unit (Using Nonlinear Method)

When the rotor locates its center position inx-axis and y-axis and the coil currents of the RMB units are 0 A, the curves of the net force-coil current-rotor position of the AMB unit in z-axis using nonlinear method are shown in Figure 14. It is obvious that the relationship of the net force-displacement in z-axis is approximately linear and the relationship of the net force-versus coil current is approximately linear for small displacement and nonlinear for large coil current. The displacement stiffness affected by the rotor position and the coil current in z-axis of the AMB unit is shown in Figure 15(a). The current stiffness affected by the rotor position and the coil current in z-axis of the AMB unit is shown in Figure 15(b). In addition, the current stiffness is only affected by the coil current of the AMB unit.Figure 14
The force-coil current and rotor position inz-axis of the AMB unit (using nonlinear method).The curves of the stiffness versus rotor position inz-axis and coil current of the AMB unit (using nonlinear method). (a) The displacement stiffness versus coil current and rotor position in z-axis of the AMB unit. (b) The current stiffness versus coil current and rotor position in z-axis of the AMB unit.
(a)
(b)


## 4.2.1. The Analysis Results of the Force-Displacement and Displacement Stiffness of the RMB Unit (Using Linearized Method, Nonlinear Method, and FEM)

When the control currents in radial coil currents and axial coil current are 0 A, and the rotor position changes from −0.4 mm to 0.4 mm inx-axis (or y-axis), and hold in center position in z-axial direction, the relationship between net force and the rotor position in x-axis (or y-axis) direction is shown in Figure 7(a) using linearized method, nonlinear method, and FEM. It is obvious that the force is nonlinear. The calculated results are approximately equal for small values of the rotor in x-axis (or y-axis) and there is little difference between the nonlinear model and FEM. In addition, the displacement stiffnesses using the three calculated methods are shown in Figure 7(b). The displacement stiffness in rotor center is -
5.5195
×
10
5
N
/
m using linearized method and nonlinear method, and the displacement stiffness calculated by the FEM in rotor center is -
5.3815
×
10
5
N
/
m; the calculated error between them is 2.6%.Curves of the force-rotor position and displacement stiffness inx-axis (or y-axis) of the RMB unit (using linearized method, nonlinear method, and FEM). (a) Curves of the force-rotor position of the RMB unit. (b) Curves of the displacement stiffness of the RMB unit.
(a)
(b)


## 4.2.2. The Analysis Results of the Force-Coil Current and Current Stiffness of the RMB Unit (Using Linearized Method, Nonlinear Method, and FEM)

When the rotor locates at radial and axial center and the coil current inx-axis ory-axis of the RMB unit changes from −1 A to 1 A (the coil current of the AMB unit is 0 A), the relationship between net force of the RMB unit and coil current is shown in Figure 8 using linearized method (see (11)), nonlinear method (see (1)), and FEM. It is obvious that the net force is linear, and the calculated results are approximately equal. In addition, the current stiffness of the RMB unit using the three calculated methods is 284.4 N/A (linearized method and nonlinear method) and 284.13 N/A (FEM); the calculated error is between linearized method and FEM 1.5%.Figure 8
Curves of the force-coil current inx-axis (or y-axis) of the RMB unit (using linearized method, nonlinear method, and FEM).


## 4.2.3. The Analysis Results of the Force-Coil Current and Rotor Position inx-Axis (or y-Axis) of the RMB Unit (Using Nonlinear Method)

When the rotor inz-axis locates its center position and the coil current of the AMB unit is 0 A, the curves of the net force-coil current-rotor position of the RMB unit in x-axis (or y-axis) using nonlinear method are shown in Figure 9. It is obvious that it is approximately linearity for small displacement and coil current and nonlinear for large displacement and coil current. The displacement stiffness and current stiffness affected by the rotor position and the coil current in x-axis (or y-axis) of the RMB unit are shown in Figures 10 and 11, respectively. They are approximately linear for small displacement and coil current.Figure 9
Curves of the force-coil current and rotor position inx-axis (or y-axis) of the RMB unit (using nonlinear method).Figure 10
Curves of the displacement stiffness-rotor position-coil current and rotor position inx-axis (or y-axis) of the RMB unit (using nonlinear method).Figure 11
Curves of the current stiffness-coil current and rotor position inx-axis (or y-axis) of the RMB unit (using nonlinear method).


## 4.2.4. The Analysis Results of the Force-Displacement and Displacement Stiffness of the AMB Unit (Using Linearized Method, Nonlinear Method, and FEM)

When the control currents in radial coil currents and axial coil current all are 0 A, and the rotor position changes from −0.4 mm to 0.4 mm inz-axis, and hold in center position in x-axial and y-axis direction, the relationship between net force and the rotor position in z-axis direction is shown in Figure 12(a) using linearized method (see (11)), nonlinear method (see (1)), and FEM. It is obvious that the force is also nonlinear. The calculated results are approximately equal for small values of the rotor in z-axis and there is little difference between the nonlinear model and FEM. In addition, the displacement stiffness using the three calculated methods is shown in Figure 12(b). The calculated results are approximately equal for small value of the displacement in z-axis between the linearized method and FEM. The displacement stiffness in rotor center is -
8.5007
×
10
5
N
/
m using linearized method and nonlinear method, and the displacement stiffness calculated by the FEM in rotor center is -
8.3732
×
10
5 N/m; the calculated error between them is 4.4%.The force-rotor position inz-axis of the AMB unit and its displacement stiffness (using linearized method, nonlinear method, and FEM). (a) The force-rotor position in z-axis of the AMB unit. (b) The displacement stiffness of the AMB unit.
(a)
(b)


## 4.2.5. The Analysis Results of the Force-Coil Current and Current Stiffness of the AMB Unit (Using Linearized Method, Nonlinear Method, and FEM)

When the rotor locates at radial and axial center and the coil current inz-axis of the AMB unit changes from −1 A to 1 A (the coil currents of the RMB unit are 0 A), the relationship between net force of AMB unit and its coil current is shown in Figure 13 using linearized method (see (15)), nonlinear method (see (13)), and FEM. It is obvious that the net force is linear, and the calculated results are approximately equal. In addition, the current stiffness of the AMB unit using the three calculated methods is 290.5 N/A (linearized method and nonlinear method) and 280.3 N/A (FEM); the calculated error is 3.5% between linearized method and FEM.Figure 13
The force-coil current inz-axis of the AMB unit (using linearized method, nonlinear method, and FEM).


## 4.2.6. The Analysis Results of the Force-Coil Current and Rotor Position inz-Axis of the AMB Unit (Using Nonlinear Method)

When the rotor locates its center position inx-axis and y-axis and the coil currents of the RMB units are 0 A, the curves of the net force-coil current-rotor position of the AMB unit in z-axis using nonlinear method are shown in Figure 14. It is obvious that the relationship of the net force-displacement in z-axis is approximately linear and the relationship of the net force-versus coil current is approximately linear for small displacement and nonlinear for large coil current. The displacement stiffness affected by the rotor position and the coil current in z-axis of the AMB unit is shown in Figure 15(a). The current stiffness affected by the rotor position and the coil current in z-axis of the AMB unit is shown in Figure 15(b). In addition, the current stiffness is only affected by the coil current of the AMB unit.Figure 14
The force-coil current and rotor position inz-axis of the AMB unit (using nonlinear method).The curves of the stiffness versus rotor position inz-axis and coil current of the AMB unit (using nonlinear method). (a) The displacement stiffness versus coil current and rotor position in z-axis of the AMB unit. (b) The current stiffness versus coil current and rotor position in z-axis of the AMB unit.
(a)
(b)


## 4.3. The Performance Analysis for Tilting of the Combined 5-DOFs HMB

### 4.3.1. The Analysis Results of the Torque-Angular Position and Angular Displacement Stiffness of the Combined 5-DOFs HMB

The torque can be produced by the RMB unitsA and B. When the control currents in radial coil and axial coil all are 0 A and the rotor tilts within −5 mrad~5 mrad around x-axis or y-axis, the relationship between torque and angular position of rotor is shown in Figure 16(a) using linearized method, nonlinear method, and FEM. The calculated results are approximately equal for small values of the angular position of rotor and there is little difference between the nonlinear model and FEM. In addition, the tilting position stiffness is shown in Figure 16(b). The calculated results are approximately equal for small value of the angular position between the linearized method and FEM. The angular displacement stiffness in rotor center is −5255.6 Nm/rad.Curves of the torque-angular position of the combined 5-DOFs HMB and its tilting displacement stiffness (using linearized method, nonlinear method, and FEM). (a) Curves of the torque-angular position of the combined 5-DOFs HMB. (b) The tilting displacement stiffness.
(a)
(b)


### 4.3.2.  The Analysis Results of the Torque-Coil Current and Tilting Current Stiffness of the Combined 5-DOFs HMB (Using Linearized Method, Nonlinear Method, and FEM)

When the rotor locates at radial and axial center and the coil currents of the RMB unitsA and B in x-axis or y-axis change from −1 A to 1 A (the coil current of AMB unit is 0 A), the relationship between torque and coil current of the RMB units is shown in Figure 17 using linearized method (see (15)), nonlinear method (see (13)) and FEM. It is obvious that the torque is linear, and the calculated results are approximately equal. In addition, the tilting current stiffness using the three calculated methods is 39.25 Nm/A (linearized method and nonlinear method) and 38.07 Nm/A (FEM); the calculated error between the linearized method and FEM is 3%.Figure 17
Curves of the torque-coil current in RMB unitsA and B of the combined 5-DOFs HMB (using linearized method, nonlinear method, and FEM).


### 4.3.3. The Analysis Results of the Torque-Coil Current and Angular Position inx-Axis or y-Axis of the Combined 5-DOFs HMB (Using Linearized Method and Nonlinear Method)

When the rotor locates its center position inz-axis and the coil currents of the AMB unit are 0 A, the curves of the torque produced by the RMB units A and B versus their coil current and rotor angular position around x-axis or y-axis using nonlinear method are shown in Figure 18 using linearized method and nonlinear method. It is obvious that it is approximately linear for small values of the angular position and coil current. Moreover, the curves of the angular position stiffness versus angular position and coil currents in RMB units A and B are shown in Figure 19(a). And the curves of the tilting current stiffness versus angular position and coil currents in RMB units A and B are shown in Figure 19(b). They all are affected by the angular position and coil currents in RMB units A and B.Figure 18
Curves of the torque-coil current and angular position aroundx-axis or y-axis in the RMB units A and B of the combined 5-DOFs HMB (using linearized method and nonlinear method).Curves of the tilting displacement stiffness and the tilting current stiffness versus angular position aroundx-axis or y-axis and the coil current in the RMB units A and B of the combined 5-DOFs HMB (using linearized method and nonlinear method). (a) Tilting displacement stiffness versus angular position around x-axis or y-axis and the coil current. (b) Tilting current stiffness versus angular position around x-axis or y-axis and the coil current.
(a)
(b)


## 4.3.1. The Analysis Results of the Torque-Angular Position and Angular Displacement Stiffness of the Combined 5-DOFs HMB

The torque can be produced by the RMB unitsA and B. When the control currents in radial coil and axial coil all are 0 A and the rotor tilts within −5 mrad~5 mrad around x-axis or y-axis, the relationship between torque and angular position of rotor is shown in Figure 16(a) using linearized method, nonlinear method, and FEM. The calculated results are approximately equal for small values of the angular position of rotor and there is little difference between the nonlinear model and FEM. In addition, the tilting position stiffness is shown in Figure 16(b). The calculated results are approximately equal for small value of the angular position between the linearized method and FEM. The angular displacement stiffness in rotor center is −5255.6 Nm/rad.Curves of the torque-angular position of the combined 5-DOFs HMB and its tilting displacement stiffness (using linearized method, nonlinear method, and FEM). (a) Curves of the torque-angular position of the combined 5-DOFs HMB. (b) The tilting displacement stiffness.
(a)
(b)


## 4.3.2.  The Analysis Results of the Torque-Coil Current and Tilting Current Stiffness of the Combined 5-DOFs HMB (Using Linearized Method, Nonlinear Method, and FEM)

When the rotor locates at radial and axial center and the coil currents of the RMB unitsA and B in x-axis or y-axis change from −1 A to 1 A (the coil current of AMB unit is 0 A), the relationship between torque and coil current of the RMB units is shown in Figure 17 using linearized method (see (15)), nonlinear method (see (13)) and FEM. It is obvious that the torque is linear, and the calculated results are approximately equal. In addition, the tilting current stiffness using the three calculated methods is 39.25 Nm/A (linearized method and nonlinear method) and 38.07 Nm/A (FEM); the calculated error between the linearized method and FEM is 3%.Figure 17
Curves of the torque-coil current in RMB unitsA and B of the combined 5-DOFs HMB (using linearized method, nonlinear method, and FEM).


## 4.3.3. The Analysis Results of the Torque-Coil Current and Angular Position inx-Axis or y-Axis of the Combined 5-DOFs HMB (Using Linearized Method and Nonlinear Method)

When the rotor locates its center position inz-axis and the coil currents of the AMB unit are 0 A, the curves of the torque produced by the RMB units A and B versus their coil current and rotor angular position around x-axis or y-axis using nonlinear method are shown in Figure 18 using linearized method and nonlinear method. It is obvious that it is approximately linear for small values of the angular position and coil current. Moreover, the curves of the angular position stiffness versus angular position and coil currents in RMB units A and B are shown in Figure 19(a). And the curves of the tilting current stiffness versus angular position and coil currents in RMB units A and B are shown in Figure 19(b). They all are affected by the angular position and coil currents in RMB units A and B.Figure 18
Curves of the torque-coil current and angular position aroundx-axis or y-axis in the RMB units A and B of the combined 5-DOFs HMB (using linearized method and nonlinear method).Curves of the tilting displacement stiffness and the tilting current stiffness versus angular position aroundx-axis or y-axis and the coil current in the RMB units A and B of the combined 5-DOFs HMB (using linearized method and nonlinear method). (a) Tilting displacement stiffness versus angular position around x-axis or y-axis and the coil current. (b) Tilting current stiffness versus angular position around x-axis or y-axis and the coil current.
(a)
(b)


## 4.4. The Coupling Performance Analysis of the Combined 5-DOFs HMB

### 4.4.1. The Coupling Performance Analysis for Translation of the Combined 5-DOFs HMB

Equation (22) shows that the displacement stiffness and current stiffness of the RMB unit will be affected by the rotor position in z-axis and the coil current of the AMB unit. The effect of the rotor position in z-axis on the displacement stiffness and current stiffness of the RMB units A and B is shown in Figures 20(a) and 20(b), respectively. They all are affected by the rotor position in z-axis and their change trend is opposite.The displacement stiffness and the current stiffness of the RMB unitsA and B affected by the rotor position of the AMB unit in z-axis. (a) The displacement stiffness of the RMB units A and B affected by the rotor position of the AMB unit in z-axis. (b) The current stiffness of the RMB units A and B affected by the rotor position of the AMB unit in z-axis.
(a)
(b)The effect of the coil current of the AMB unit on the displacement stiffness and current stiffness of the RMB unitsA and B is shown in Figures 21(a) and 21(b), respectively. They all are affected by the coil current of the AMB unit and their change trend is the same.The displacement and current stiffness of the RMB unitsA and B affected by the coil current of the AMB unit. (a) The displacement stiffness of the RMB units A and B affected by the coil current of the AMB unit. (b) The current stiffness of the RMB units A and B affected by the coil current of the AMB unit.
(a)
(b)The 3D curves of the displacement stiffness and current stiffness of the RMB unitsA and B affected by the rotor position in z-axis and the coil current of the AMB unit are shown in Figures 22(a) and 22(b).The displacement stiffness and the current stiffness of the RMB unitsA and B affected by the rotor position in z-axis and the coil current of the AMB unit. (a) The displacement stiffness of the RMB units A and B affected by the rotor position in z-axis and the coil current of the AMB unit. (b) The current stiffness of the RMB units A and B affected by the rotor position in z-axis and the coil current of the AMB unit.
(a)
(b)Figure23 shows that the displacement stiffness and current stiffness of the AMB unit will be affected by the rotor position in x-axis or y-axis of the RMB unit. The effect of the rotor position in x-axis or y-axis on the displacement stiffness and current stiffness of the AMB unit is shown in Figures 23(a) and 23(b), respectively. They all are affected by the rotor position in x-axis or y-axis.The displacement and the current stiffness of the AMB unit affected by the rotor position of the RMB unit inx-axis or y-axis. (a) The displacement stiffness of the AMB unit affected by the rotor position of the RMB unit. (b) The current stiffness of the AMB unit affected by the rotor position of the RMB unit.
(a)
(b)


### 4.4.2. The Coupling Performance Analysis for Tilting of the Combined 5-DOFs HMB

The tilting displacement stiffness and tilting current stiffness of the RMB unitsA and B will be affected by the rotor position in z-axis and coil current of the AMB unit, and the influence curves of RMB units A and B are shown in Figures 24(a), 24(b), 25(a), and 25(b), respectively. It is obvious that the tilting displacement stiffness and tilting current stiffness of the RMB units A and B all are affected by the rotor position in z-axis and coil current of the AMB unit.Curves of the tilting displacement stiffness and tilting current stiffness of the RMB unitsA and B affected by the rotor position in z-axis of the AMB unit. (a) The tilting displacement stiffness of the RMB units A and B affected by the rotor position in z-axis. (b) The tilting current stiffness of the RMB units A and B affected by the rotor position in z-axis of the AMB unit.
(a)
(b)Curves of the tilting displacement stiffness and tilting current stiffness of the RMB unitsA and B affected by the rotor position in z-axis and coil current of the AMB unit. (a) Curves of the tilting displacement stiffness of the RMB units A and B affected by the rotor position in z-axis and coil current of the AMB unit. (b) Curves of the tilting current stiffness of the RMB units A and B affected by the rotor position in z-axis and coil current of the AMB unit.
(a)
(b)


## 4.4.1. The Coupling Performance Analysis for Translation of the Combined 5-DOFs HMB

Equation (22) shows that the displacement stiffness and current stiffness of the RMB unit will be affected by the rotor position in z-axis and the coil current of the AMB unit. The effect of the rotor position in z-axis on the displacement stiffness and current stiffness of the RMB units A and B is shown in Figures 20(a) and 20(b), respectively. They all are affected by the rotor position in z-axis and their change trend is opposite.The displacement stiffness and the current stiffness of the RMB unitsA and B affected by the rotor position of the AMB unit in z-axis. (a) The displacement stiffness of the RMB units A and B affected by the rotor position of the AMB unit in z-axis. (b) The current stiffness of the RMB units A and B affected by the rotor position of the AMB unit in z-axis.
(a)
(b)The effect of the coil current of the AMB unit on the displacement stiffness and current stiffness of the RMB unitsA and B is shown in Figures 21(a) and 21(b), respectively. They all are affected by the coil current of the AMB unit and their change trend is the same.The displacement and current stiffness of the RMB unitsA and B affected by the coil current of the AMB unit. (a) The displacement stiffness of the RMB units A and B affected by the coil current of the AMB unit. (b) The current stiffness of the RMB units A and B affected by the coil current of the AMB unit.
(a)
(b)The 3D curves of the displacement stiffness and current stiffness of the RMB unitsA and B affected by the rotor position in z-axis and the coil current of the AMB unit are shown in Figures 22(a) and 22(b).The displacement stiffness and the current stiffness of the RMB unitsA and B affected by the rotor position in z-axis and the coil current of the AMB unit. (a) The displacement stiffness of the RMB units A and B affected by the rotor position in z-axis and the coil current of the AMB unit. (b) The current stiffness of the RMB units A and B affected by the rotor position in z-axis and the coil current of the AMB unit.
(a)
(b)Figure23 shows that the displacement stiffness and current stiffness of the AMB unit will be affected by the rotor position in x-axis or y-axis of the RMB unit. The effect of the rotor position in x-axis or y-axis on the displacement stiffness and current stiffness of the AMB unit is shown in Figures 23(a) and 23(b), respectively. They all are affected by the rotor position in x-axis or y-axis.The displacement and the current stiffness of the AMB unit affected by the rotor position of the RMB unit inx-axis or y-axis. (a) The displacement stiffness of the AMB unit affected by the rotor position of the RMB unit. (b) The current stiffness of the AMB unit affected by the rotor position of the RMB unit.
(a)
(b)


## 4.4.2. The Coupling Performance Analysis for Tilting of the Combined 5-DOFs HMB

The tilting displacement stiffness and tilting current stiffness of the RMB unitsA and B will be affected by the rotor position in z-axis and coil current of the AMB unit, and the influence curves of RMB units A and B are shown in Figures 24(a), 24(b), 25(a), and 25(b), respectively. It is obvious that the tilting displacement stiffness and tilting current stiffness of the RMB units A and B all are affected by the rotor position in z-axis and coil current of the AMB unit.Curves of the tilting displacement stiffness and tilting current stiffness of the RMB unitsA and B affected by the rotor position in z-axis of the AMB unit. (a) The tilting displacement stiffness of the RMB units A and B affected by the rotor position in z-axis. (b) The tilting current stiffness of the RMB units A and B affected by the rotor position in z-axis of the AMB unit.
(a)
(b)Curves of the tilting displacement stiffness and tilting current stiffness of the RMB unitsA and B affected by the rotor position in z-axis and coil current of the AMB unit. (a) Curves of the tilting displacement stiffness of the RMB units A and B affected by the rotor position in z-axis and coil current of the AMB unit. (b) Curves of the tilting current stiffness of the RMB units A and B affected by the rotor position in z-axis and coil current of the AMB unit.
(a)
(b)


## 5. Conclusion

A novel combined 5-DOFs HMB with only one PMR is proposed, used for turboexpander, in this paper. The combined 5-DOFs HMB consists of two RMB units and one AMB unit. The bias flux of the two RMB units and AMB unit is supplied by only one PMR. The low power consumption, compact structure, and small size can be gained due to the combined 5-DOFs HMB, the permanent magnetic bias HMB, and the reduction in the number of the poles of the RMB. The low cost is gained due to the one PMR used. Moreover, ultra-high-speed is also gained due to the small outside diameter of the rotor thrust disk. The equivalent magnetic circuit method and 3D finite element method (FEM) are used to model and analyze the integral 5-DOFs HMB. The force-current, force-position, torque-coil currents, torque-angular position, and the stiffness models of the integral 5-DOFs HMB are given. Moreover, its coupling problems between the RMB units and the AMB unit are also proposed. An example is given to clarify the mathematical models and the coupling problems, and the linearized models are proposed for the follow-up controller design.



---

*Source: 162561-2014-05-07.xml*