# Modeling and Analysis of Coupling Performance of Dynamic Stiffness Models for a Novel Combined Radial-Axial Hybrid Magnetic Bearing

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

---

## Abstract

The combined radial-axial magnetic bearing (CRAMB) with permanent magnet creating bias flux can reduce the size, cost, and mass and save energy of the magnetic bearing. The CRAMB have three-degree-of-freedom control ability, so its structure and magnetic circuits are more complicated compared to those of the axial magnetic bearing (AMB) or radial magnetic bearing (RMB). And the eddy currents have a fundamental impact on the dynamic performance of the CRAMB. The dynamic stiffness model and its cross coupling problems between different degrees of freedom affected for the CRAMB are proposed in this paper. The dynamic current stiffness and the dynamic displacement stiffness models of the CRAMB are deduced by using the method of equivalent magnetic circuit including eddy current effect, but the dynamic current stiffness of the RMB unit is approximately equal to its static current stiffness. The analytical results of an example show that the bandwidth of the dynamic current stiffness of the AMB unit and the dynamic displacement stiffness of the CRAMB is affected by the time-varying control currents or air gap, respectively. And the dynamic current stiffness and the dynamic displacement stiffness between the AMB unit and the RMB unit are decoupled due to few coupling coefficients.


---

## Body

## 1. Introduction

Most of 5-axis active magnetic bearing systems (MBSs) are usually composed of two radial magnetic bearing (RMB) units and one axial magnetic bearing (AMB) unit [1–5]. These magnetic bearing systems are the easiest way to be produced, but they are also tending to bulky, high power and high cost. In order to reduce the size, cost, and save energy and increase the high-power density of the MBS, which is important to reduce the number of the units by means of furthermore combination the AMB and RMB. A combined radial-axial magnetic bearing (CRAMB) which is named as 3-axis magnetic bearing (MB) is designed for use in an ultra-high-speed machine [6]. The magnetic forces and coupling problems of a combined AMB and RMB are analyzed [7]. The integrated AMB and RMB with conical rotor are designed and analyzed [8, 9], and one downside of the integrated bearings is a strong coupling problem between the radial and axial degrees of freedom. An AC-DC 3-DOF hybrid magnetic bearing is proposed and designed [10]. Structure and control method of an AC-DC 3-DOF hybrid magnetic bearing is introduced in the literature [11–14]. A 3-DOF axial hybrid magnetic bearing [15] is proposed, but the structure and processing technic are very complicated, and its rotational power loss will be large at high speed. The structure of a 3-DOF magnetic bearing without large thrust disk rotor is introduced and designed [16]. An integrated radial-axial magnetic bearing without axial disk is presented in the literature, and its static performance is analyzed by using 3-D FEM. The static mathematical model of a radial-axial magnetic bearing is introduced in the literatures [17, 18].Most literatures mentioned above described the structure, design, analysis, and control methods of the CRAMBs, which are mainly consider the static mathematical model without including the eddy-current effects. But few literatures on the dynamic stiffness modeling and coupling problems of CRAMB with permanent magnet bias are described when considering the dynamic characteristics affected by the eddy currents caused by time-varying air gap and control currents. Some parts of the CRAMB are very difficult to be split into laminated which can reduce the eddy-current loss, then the eddy currents will result in magnitude reduction and phase lag of radial and axial magnetic force and make the core hotter. That is to say, the eddy current has a serious impact on the dynamic performances of the CRAMB.Due to the coupled magnetic, electrical, and mechanical domains, design and analysis of CRAMB are very difficult. In particular, its dynamic characteristics are affected by the coupled domains mentioned above. As a result, to predict the dynamic operating characteristics is necessary to avoid costly trials. A few simple dynamic models including eddy currents are described for magnetic actuators [19–25]. The finite element method (FEM) is also used to analyze the eddy-current effects on the active thrust MBs. The electromechanical performance of magnetic actuators is analyzed using a coupled structural and electromagnetic FEM [26, 27]. The dynamic characteristics of MB based on field-circuit coupled method are proposed [28]. The FEM is used to analyze the operating characteristics of the electromagnet with permanent magnets [29]. The dynamic stiffness models of active thrust magnetic bearing are also analyzed by FEM and analytical method [30, 31]. The eddy current loss in surface layer of a laminated core [32] and a solid core [33] for radial MB is analyzed. And the experimental methods are used to investigate the eddy-current effects. The force/current relationship of a solid thrust MB is tested by Allaire et al. [34]. DeWeese et al. investigated different rotor and stator configurations to find the one with minimal eddy currents [35]. The limitations on the closed-loop performance of active thrust MB affected by the eddy currents are investigated [36].The emphases of dynamic stiffness mentioned above are the active thrust MB which has one degree of freedom (DOF), and the structure and magnetic flux path are easy to design and analyze. The structure and magnetic circuits of the CRAMB are more complicated compared to those of the thrust MB or radial MB, and the cross-coupling problems between different DOFs may be produced. But few literatures on dynamic characteristics and the coupling problems of the CRAMB are published. So, it would have higher research value to investigate the accurate dynamic factor model and the coupling problems including eddy currents rather than the static mathematic model. The dynamic stiffness models and the dynamic coupling models including eddy-current effects for CRAMB with permanent magnet bias are presented and derived analytically by using equivalent magnetic circuit method in this paper. The dynamic current stiffness models and the dynamic displacement stiffness models of the RMB unit and AMB unit are given. And the dynamic coupling models between axial and radial bearing units are also given. An example is also given in this paper.


## 2. Analysis Model and Method

### 2.1. Analysis Model

Figure1 shows the configuration of the CRAMB which combines the radial and axial magnetic bearing unit. The axially magnetized permanent magnet ring (PMR) is used to produce bias flux for both the RMB and the AMB unit. The laminations are used only for radial stator and rotor core, and the other parts are made of solid material. Smaller outside diameter of the rotor thrust disk is used. Therefor the lower rotational drag, the lower stresses, and the compact structure are attained. The magnetic ring is used to eliminate the saturation of the magnetic circuit for the stator core of the RMB unit. The bias flux and the control flux paths of the RMB and AMB unit are depicted in Figure 2. The solid line indicates the bias flux paths generated by the PMR and the dotted line denotes the control flux paths generated by the control current of the AMB and the RMB unit. As can be seen, their control flux produced by the control current in coils will pass through air gap of the RAM and the AMB rather than the PMR, and therefore small control currents are enough since they do not have to compensate for the high magnetic reluctance of the PMR. The control flux has to be superimposed by the bias flux in air gap for producing the axial and radial magnetic force.Figure 1
Configuration of the CRAMB with permanent magnet ring.Figure 2
The magnetic circuits of the combined radial-axial magnetic bearing with permanent magnet ring.The radial and axial magnetic bearing force and flux can be derived conveniently with the static magnetic circuit model without consideration of the eddy-current effect. However, the eddy currents caused by the time-varying air gap and the control currents will affect the dynamic characteristics of the CRAMB system. Based on an approximate analysis method [22] considering the eddy-current effect inside a nonlaminated electromagnetic actuator, the CRAMB geometry is divided into thirteen elements as shown in Figure 3. The four elements of the air gap between rotor and stator of the RMB are named as element 13. The eddy-current loss can be reduced by using laminated stator core and rotor core of the RMB unit in CRAMB. According to the equivalent magnetic circuit models of the CRAMB, the effective reluctances are separated into two parts: one is the static magnetic reluctance of element, and the other one is the dynamic magnetic reluctance which is a half-order term related to frequency.Figure 3
Division of the geometry of the CRAMB: 13 elements simply divided.The elements 1–5, 7, and 9–12 can be divided in the same way as in Sun et al.’s models [30] and Zhu et al.’s models [23]. The PMR is used for element 2 which is used to provide bias flux and can also be simplified in the same way. The elements 1, 2, 11, and 12 are located in the bias flux path. Therefore, the eddy-current effect should be considered in mathematical model. The simplified effective reluctance of each element can be used to design, analyze, and control purpose. The simplified calculation of the effective magnetic reluctances for these elements is given as follow:
(1)
R
1
=
R
1
0
+
R
1
i
=
h
1
π
μ
0
μ
r
(
r
6
2
-
r
5
2
)
+
h
1
2
π
r
5
j
ω
σ
1
μ
0
μ
r
,
where R
1
0
=
h
1
/
π
μ
0
μ
r
(
r
6
2
-
r
5
2
)
,
R
1
i
=
(
h
1
/
2
π
r
5
)
j
ω
σ
1
/
μ
0
μ
r
(2)
R
2
=
R
2
0
+
R
2
i
=
h
2
π
μ
0
μ
r
(
r
6
2
-
r
5
2
)
+
h
2
2
π
r
5
j
ω
σ
2
μ
0
μ
r
,
where R
2
0
=
h
2
/
π
μ
0
μ
r
(
r
6
2
-
r
5
2
)
,
R
2
i
=
(
h
2
/
2
π
r
5
)
j
ω
σ
2
/
μ
0
μ
r
(3)
R
3
=
R
3
0
+
R
3
i
=
h
4
π
μ
0
μ
r
(
r
6
2
-
r
5
2
)
+
h
4
2
π
r
5
j
ω
σ
1
μ
0
μ
r
,
where R
3
0
=
h
4
/
π
μ
0
μ
r
(
r
6
2
-
r
5
2
),  R
3
i
=
(
h
4
/
2
π
r
5
)
j
ω
σ
1
/
μ
0
μ
r
(4)
R
4
=
R
4
0
+
R
4
i
=
ln
(
r
6
/
r
4
)
2
π
μ
0
μ
r
h
5
+
ln
(
r
6
/
r
4
)
2
π
j
ω
σ
1
μ
0
μ
r
,
where R
4
0
=
ln
(
r
6
/
r
4
)
/
2
π
μ
0
μ
r
h
5,  R
4
i
=
(
ln
(
r
6
/
r
4
)
/
2
π
)
j
ω
σ
1
/
μ
0
μ
r
(5)
R
5
=
R
5
0
+
R
5
i
=
h
6
π
μ
0
μ
r
(
r
4
2
-
r
3
2
)
+
h
6
2
π
r
3
j
ω
σ
1
μ
0
μ
r
,
where R
5
0
=
h
6
/
π
μ
0
μ
r
(
r
4
2
-
r
3
2
), R
5
i
=
(
h
6
/
2
π
r
3
)
j
ω
σ
1
/
μ
0
μ
r
(6)
R
z
1
=
R
6
=
s
0
+
z
π
μ
0
(
r
4
2
-
r
3
2
)
,
R
z
10
=
R
6
0
=
s
0
π
μ
0
(
r
4
2
-
r
3
2
)
,
z
=
0
,
R
7
=
R
7
0
+
R
7
i
=
ln
(
r
4
/
r
3
)
2
π
μ
0
μ
r
h
7
+
ln
(
r
4
/
r
3
)
2
π
j
ω
σ
1
μ
0
μ
r
,
where R
7
0
=
ln
(
r
4
/
r
3
)
/
2
π
μ
0
μ
r
h
7,  R
7
i
=
(
ln
(
r
4
/
r
3
)
/
2
π
)
j
ω
σ
1
/
μ
0
μ
r
(7)
R
z
2
=
R
8
=
s
0
-
z
π
μ
0
(
r
4
2
-
r
3
2
)
,
R
z
0
=
R
z
20
=
R
8
0
=
R
z
10
=
s
0
π
μ
0
(
r
4
2
-
r
3
2
)
,
z
=
0
,
R
9
=
R
9
0
+
R
9
i
=
R
5
=
R
5
0
+
R
5
i
,
where R
9
0
=
h
6
/
π
μ
0
μ
r
(
r
4
2
-
r
3
2
), R
9
i
=
(
h
6
/
2
π
r
3
)
j
ω
σ
1
/
μ
0
μ
r
(8)
R
10
=
R
10
0
+
R
10
i
=
R
4
=
R
4
0
+
R
4
i
,
where R
10
0
=
ln
(
r
6
/
r
4
)
/
2
π
μ
0
μ
r
h
5, R
10
i
=
(
ln
(
r
6
/
r
4
)
/
2
π
)
j
ω
σ
1
/
μ
0
μ
r
(9)
R
11
=
R
11
0
+
R
11
i
=
ln
(
r
3
/
r
1
)
2
π
μ
0
μ
r
h
7
+
ln
(
r
3
/
r
1
)
2
π
j
ω
σ
1
μ
0
μ
r
,
where R
11
0
=
ln
(
r
3
/
r
1
)
/
2
π
μ
0
μ
r
h
7, R
11
i
=
(
ln
(
r
3
/
r
1
)
/
2
π
)
j
ω
σ
1
/
μ
0
μ
r
(10)
R
12
=
R
12
0
+
R
12
i
=
h
9
π
μ
0
μ
r
(
r
1
2
-
r
8
2
)
+
h
9
2
π
r
8
j
ω
σ
1
μ
0
μ
r
,
where R
12
0
=
h
9
/
π
μ
0
μ
r
(
r
1
2
-
r
8
2
), R
12
i
=
(
h
9
/
2
π
r
8
)
j
ω
σ
1
/
μ
0
μ
r
(11)
R
13
=
R
13
0
=
360
(
r
3
-
r
2
)
π
μ
0
h
10
(
r
3
+
r
2
)
,
where R
j
0 is the static magnetic reluctance (the first term in R
j for element j) and R
j
i is the eddy-current magnetic reluctance (the second term in R
j for element j) with j
=
1,2
,
…
,
1
3, and the parameters and their dimensions are given in Figure 4 and Table 1. R
13 is the magnetic reluctance of the radial air gap. R
6 and R
8 are the magnetic reluctances of the axial air gap. The equivalent magnetic circuits including the bias magnetic circuit and the control magnetic circuits of the CRAMB without the leakage flux are shown in Figure 5. We focus on the investigation of the eddy-current effects, and the equivalent magnetic reluctances for the leakage flux paths would be included in the future model.Table 1
Parameters of The CRAMB.
Parameters
Description
Value
r
1
Inner radius of the RMB rotor, (mm)
25
r
2
Outer radius of the RMB rotor, (mm)
30
r
7
Inner radius of the RMB stator, (mm)
30.5
A
r
Radial pole shoe area, (mm2)
1036
N
r
The number of winding turns in coil of RMB
220
r
9
Outer radius of the RMB stator, (mm)
60
r
6
Outer radius of the yoke, (mm)
70
r
3
Inner radius of the axial pole shoe, (mm)
31
r
4
Outer radius of the axial pole shoe, (mm)
34
r
5
Inner radius of the yoke, (mm)
62
r
8
Inner radius of the RMB rotor, (mm)
10
h
1
Length of the stator return ring of RMB, (mm)
31
h
2
Thickness of the permanent ring, (mm)
10
h
3
Thickness of the element 10, (mm)
8
h
4
Thickness of the element 3, (mm)
24
h
5
Thickness of the element 4, (mm)
8
h
6
Thickness of the element 5, (mm)
8
h
7
Thickness of the thrust plane, (mm)
7
h
8
Thickness of the element 9, (mm)
8
h
9
Thickness of the element 12, (mm)
8
s
0
radial and axial air gap, (mm)
0.5
N
z
The number of winding turns in coil of AMB
120
m
The mass of the rotor, (kg)
7.2
σ
1
The conductivity of yoke and thrust disk, (S/m)
8 × 106
μ
0
Permeability of free space, (H/m)
4
π
×
10
-
7
μ
r
Relative permeability of the permanent magnet
1.082
σ
2
Resistivity of the magnet ring, (Ω·m)
0.9
E
-
6Figure 4
The dimension of the CRAMB.The equivalent magnetic circuits of the CRAMB. (a) The equivalent magnetic circuit of the bias flux paths. (b) The equivalent magnetic circuit of control flux path of the RMB unit. (c) The equivalent magnetic circuit of control flux path of the AMB unit.
(a)
(b)
(c)


### 2.2. Dynamic Current Stiffness Model for the AMB Unit of CRAMB

It is assumed that the CRAMB system is linear, and the AMB unit of the CRAMB excited a sinusoidal varying ac currenti
z
=
I
z
sin
(
ω
t
) (where the amplitude, I
z, is very small). When the thrust disk located its center position, the length of air gap (the element 6, 8) is s
0.The control fluxϕ
z
i produced by a sinusoidal varying ac current i
z in the axial air gap of AMB unit and the bias flux ϕ
z
m
0 in the axial air gap of AMB unit produced by the PMR, the control flux in the left air gap (element 8) will add to the bias flux, and will subtract from the bias flux in the right air gap (element 6) as shown in Figure 3. The total flux, ϕ
z
sum
i
(
i
=
1
,
2
), in the left and the right air gap of AMB unit is given as
(12)
ϕ
z
sum
1
=
ϕ
z
m
0
+
ϕ
z
i
,
ϕ
z
sum
2
=
ϕ
z
m
0
-
ϕ
z
i
,
where the static bias flux in the axial air gap of AMB unit is calculated as
(13)
ϕ
z
m
=
F
m
(
1
1
/
∑
k
=
3
6
R
k
0
+
1
/
∑
j
=
8
10
R
j
0
R
1
0
+
R
2
0
+
R
7
0
+
R
11
0
+
R
12
0
+
R
13
4
+
1
1
/
∑
k
=
3
6
R
k
0
+
1
/
∑
j
=
8
10
R
j
0
)
-
1
.The permanent magnet material of the PMR is Nd-Fe-B, its magnetomotive force isF
m
=
H
c
h
2, and H
c is the coercive force of the PMR; h
2 is thickness of the PMR. The permeability of the magnetic ring and the thrust disk is compared to that of the air and the PMB, so (13) can be simplified as
(14)
ϕ
z
m
≈
F
m
R
2
0
+
R
13
/
4
+
R
z
0
/
2
.And the static bias flux in the air gap of the AMB can be calculated as(15)
ϕ
z
m
0
=
ϕ
z
1
0
=
ϕ
z
2
0
=
ϕ
z
m
2
=
F
m
2
R
2
0
+
R
13
/
2
+
R
z
0
.The dynamic flux in the air gap of AMB unit produced by the ac current is calculated as(16)
ϕ
z
i
=
N
z
I
z
sin
(
ω
t
)
∑
k
=
1
10
R
k
=
N
z
I
z
R
z
i
0
sin
(
ω
t
)
1
+
R
e
(
ω
)
/
R
z
i
0
=
N
z
I
z
sin
(
ω
t
)
R
z
i
0
(
1
+
λ
)
,
where R
z
i
0 is the total static magnetic reluctance of the equivalent magnetic circuit of control flux path for AMB unit, R
e
(
ω
) is the total eddy-current magnetic reluctance of the equivalent magnetic circuit of control flux path for AMB unit, and λ is the ratio of the total eddy-current magnetic reluctance to the total static magnetic reluctance R
e
(
ω
)/R
z
i
0:
(17)
R
z
i
0
=
R
6
0
+
R
8
0
,
R
e
(
ω
)
=
R
3
i
+
R
4
i
+
R
5
i
+
R
7
i
+
R
9
i
+
R
10
i
.R
e
(
ω
) is the magnetic reluctance including the eddy current effect, and it will be increased with the addition of frequency of the time-varying air gap and the control currents. The bandwidth will be increased with the addition of the air gap.The AMB unit of the CRAMB will produce a net restoring forcef
z
i on the thrust disk in z-axis direction. Since the net force, f
z
i, is calculated by
(18)
f
z
i
=
ϕ
z
1
2
2
μ
0
A
z
-
ϕ
z
2
2
2
μ
0
A
z
,
where the constant μ
0 is the permeability of free space, A
z is the area of pole face of the AMB unit; substituting (12) into (18) gives
(19)
f
z
i
=
(
ϕ
z
m
0
+
ϕ
z
i
)
2
2
μ
0
A
z
-
(
ϕ
z
m
0
-
ϕ
z
i
)
2
2
μ
0
A
z
=
2
ϕ
z
m
0
ϕ
z
i
μ
0
A
z
.The dynamic current stiffness of the AMB unit can be derived as(20)
k
z
i
i
=
d
f
z
i
d
i
z
=
2
ϕ
z
m
0
N
z
μ
0
A
z
R
z
i
0
(
1
+
λ
)
=
k
z
i
0
1
+
λ
,
where
(21)
k
z
i
0
=
2
ϕ
z
m
0
N
z
μ
0
A
z
R
z
i
0
is the static current stiffness of the AMB unit.


### 2.3. Dynamic Displacement Stiffness Model for the AMB Unit of CRAMB

When the coil currents in the RMB and the AMB unit are zero, the air gap in the RMB remains stationary, but the air gap in the AMB unit varies sinusoidally about the nominal values
0; that is, z
1
=
s
0
-
δ
sin
(
ω
t
) in left air gap, and z
2
=
s
0
+
δ
sin
(
ω
t
) in right air gap; then the magnetic reluctance of the air gap will also vary sinusoidally about the nominal value. Thus, the reluctance of the two air gaps in AMB unit includes two parts, one is the static reluctance, R
z
i
0
(
i
=
1
,
2
), and the other one is the dynamic reluctance, R
z
i
d
(
i
=
1,2
). Consider the following:
(22)
R
z
1
=
R
z
1
0
+
R
z
1
d
=
s
0
μ
0
A
z
+
δ
sin
(
ω
t
)
μ
0
A
z
,
R
z
2
=
R
z
2
0
+
R
z
2
d
=
s
0
μ
0
A
z
-
δ
sin
(
ω
t
)
μ
0
A
z
,
where
(23)
R
z
1
0
=
R
z
2
0
=
s
0
μ
0
A
z
,
R
z
1
d
=
-
R
z
2
d
=
δ
sin
(
ω
t
)
μ
0
A
z
.If all the eddy-current reluctances and leakage reluctances are ignored, the total static magnetic reluctanceR
met in the bias magnetic circuit is
(24)
R
met
=
R
1
0
+
R
2
0
+
R
7
0
+
R
11
0
+
R
12
0
+
R
13
4
+
1
1
/
∑
k
=
3
6
R
k
0
+
1
/
∑
j
=
8
10
R
j
0
.Since the relative permeability of the stator and rotor core of the AMB and the RMB unit is larger, their static magnetic reluctance can be ignored in the equivalent magnetic circuit of the bias flux path for calculating the dynamic displacement stiffness. The total static magnetic reluctance in (24) can be simplified as
(25)
R
met
=
R
met
0
+
R
met
e
≈
R
2
0
+
R
13
4
+
1
1
/
R
z
1
+
1
/
R
z
2
,
where,
(26)
R
met
0
=
R
2
0
+
R
13
4
+
R
z
1
0
R
z
2
0
R
z
1
0
+
R
z
2
0
,
R
met
e
=
(
R
z
2
0
-
R
z
1
0
)
R
z
1
d
-
(
R
z
1
d
)
2
R
z
1
0
+
R
z
2
0
.Sinceδ is smaller compared to s
0, then R
z
i
d
(
i
=
1
,
2
) is smaller compared to R
z
i
0
(
i
=
1
,
2
), and R
met
e in (26) can be expanded into the power series and approximated by the first two terms as the higher order terms can be ignored. Consider the following:
(27)
R
met
0
=
R
2
0
+
R
13
4
+
R
z
1
0
R
z
2
0
R
z
1
0
+
R
z
2
0
,
R
met
e
≈
(
R
z
2
0
-
R
z
1
0
)
R
z
1
d
R
z
1
0
+
R
z
2
0
.The flux in the air gap of AMB unit can be calculated as(28)
ϕ
z
1
=
F
m
R
met
0
+
R
met
e
1
1
/
(
R
z
1
0
+
R
z
1
d
)
+
1
/
(
R
z
2
0
-
R
z
1
d
)
1
R
z
1
0
+
R
z
1
d
,
ϕ
z
2
=
F
m
R
met
0
+
R
met
e
1
1
/
(
R
z
1
0
+
R
z
1
d
)
+
1
/
(
R
z
2
0
-
R
z
1
d
)
1
R
z
2
0
-
R
z
1
d
.Substituting (27) for the right side of (28), since R
z
i
d is smaller compared to R
z
i
0, (28) can be expanded into the power series and approximated by the first two terms as the higher order terms can be ignored. Consider the following:
(29)
ϕ
z
1
≈
F
m
R
met
0
R
z
1
0
R
z
1
0
+
R
z
2
0
[
1
-
R
z
1
d
R
met
0
(
R
z
2
0
-
R
z
1
0
R
z
1
0
+
R
z
2
0
)
-
R
z
1
d
R
z
2
0
]
,
ϕ
z
2
≈
F
m
R
met
0
R
z
2
0
R
z
1
0
+
R
z
2
0
[
1
-
R
z
1
d
R
met
0
(
R
z
2
0
-
R
z
1
0
R
z
1
0
+
R
z
2
0
)
+
R
z
1
d
R
z
1
0
]
.The Equation (29) can be modified due to R
z
1
0
=
R
z
2
0.                              (30)
ϕ
z
1
≈
F
m
2
R
met
0
(
1
-
A
R
z
1
d
R
met
0
-
R
z
1
d
R
z
2
0
)
=
ϕ
z
0
(
1
-
A
R
z
1
d
R
met
0
-
R
z
1
d
R
z
2
0
)
,
ϕ
z
2
≈
F
m
2
R
met
0
(
1
-
A
R
z
1
d
R
met
0
+
R
z
1
d
R
z
1
0
)
=
ϕ
z
0
(
1
-
A
R
z
1
d
R
met
0
+
R
z
1
d
R
z
1
0
)
,
where A
=
(
R
z
2
e
t
-
R
z
1
e
t
)
/
(
R
z
1
0
+
R
z
2
0
).The AMB unit of the CRAMB will produce a net restoring forcef
z
i on the thrust disk in z-axis direction, since the force is given by
(31)
f
z
s
=
(
ϕ
z
1
)
2
2
μ
0
A
z
-
(
ϕ
z
2
)
2
2
μ
0
A
z
=
ϕ
z
0
2
2
μ
0
A
z
×
[
(
1
-
A
R
z
1
d
R
met
0
-
R
z
1
d
R
z
2
0
)
2
-
(
1
-
A
R
z
1
d
R
met
0
+
R
z
1
d
R
z
1
0
)
2
]
,
where ϕ
z
0
=
ϕ
z
1
0
=
ϕ
z
2
0
=
F
m
/
2
R
met
0 is the static bias flux in the air gap of the AMB unit.SinceR
z
i
d
(
i
=
1
,
2
) is smaller compared to R
z
i
0
(
i
=
1
,
2
), (31) can be expanded into the power series and approximated by the first two terms as the higher order terms can be ignored. Consider the following:
(32)
f
z
s
≈
-
2
ϕ
z
0
2
(
μ
0
A
z
)
2
δ
sin
(
ω
t
)
R
z
1
0
.Static displacement stiffness of the AMB unit is given as(33)
k
z
s
0
=
d
f
z
s
d
(
δ
sin
(
ω
t
)
)
=
-
2
ϕ
z
0
2
μ
0
A
z
s
0
=
-
1
2
μ
0
A
z
s
0
(
F
m
R
2
0
+
R
13
/
4
+
R
z
1
0
/
2
)
2
.In the static analysis model mentioned above, the eddy-current magnetic reluctance is not considered. However, the dynamic field should also cause eddy current, and the static model could be considered to include eddy-current magnetic reluctance as in the dynamic model. Then the dynamic flux including eddy-current effect is given as(34)
ϕ
z
1
d
=
ϕ
z
1
0
(
1
-
B
R
z
1
d
R
sum
e
t
-
R
z
1
d
R
z
2
e
t
)
,
ϕ
z
2
d
=
ϕ
z
2
0
(
1
-
B
R
z
1
d
R
sum
e
t
+
R
z
1
d
R
z
1
e
t
)
,
where,
(35)
R
sum
e
t
=
R
sum
0
+
R
sum
ω
,
R
sum
0
=
R
1
0
+
R
2
0
+
R
7
0
+
R
11
0
+
R
12
0
+
R
13
4
+
R
z
1
0
2
,
R
sum
ω
=
∑
m
=
1
2
R
m
i
+
R
7
i
+
∑
n
=
11
12
R
n
i
+
1
1
/
∑
j
=
3
6
R
j
+
1
/
∑
k
=
8
10
R
k
-
R
z
1
0
2
,
B
=
R
z
2
e
t
-
R
z
1
e
t
R
z
1
e
t
+
R
z
2
e
t
,
R
z
1
e
t
=
∑
j
=
3
6
R
j
,
R
z
2
e
t
=
∑
k
=
8
10
R
k
.The permeability of the magnetic ring, stator core, rotor stator, and the thrust disk is compared to that of the air and the PMR, so (35) can be modified as
(36)
R
sum
0
≈
R
2
0
+
R
13
4
+
R
z
1
0
2
,
R
sum
ω
≈
∑
m
=
1
2
R
m
i
+
R
7
i
+
∑
n
=
11
12
R
n
i
+
1
1
/
∑
j
=
3
6
R
j
i
+
1
/
∑
k
=
8
10
R
k
i
-
R
z
1
0
2
,
R
z
1
e
t
=
R
3
i
+
R
4
i
+
R
5
i
+
R
6
0
,
R
z
2
e
t
=
R
8
0
+
R
9
i
+
R
10
i
.The net force of the AMB unit considering the eddy-current magnetic reluctance is given by(37)
f
z
s
=
(
ϕ
z
1
d
)
2
2
μ
0
A
z
-
(
ϕ
z
2
d
)
2
2
μ
0
A
z
=
ϕ
z
0
2
2
μ
0
A
z
[
(
1
-
B
R
z
1
d
R
sum
e
t
-
R
z
1
d
R
z
2
e
t
)
2

-
(
1
-
B
R
z
1
d
R
sum
e
t
+
R
z
1
d
R
z
1
e
t
)
2
]

.SinceR
z
i
d is smaller compared to R
z
i
0, (37) can be expanded into the power series and approximated by the first two terms as the higher order terms can be ignored. Consider the following:
(38)
f
z
s
d
≈
-
2
ϕ
z
0
2
(
μ
0
A
z
)
2
(
1
2
R
z
2
e
t
+
1
2
R
z
1
e
t
)
δ
sin
(
ω
t
)
.The dynamic displacement stiffness considering the eddy-current magnetic reluctance is calculated as(39)
k
z
s
d
=
d
f
z
s
d
d
(
δ
sin
(
ω
t
)
)
=
-
2
ϕ
z
0
2
μ
0
A
z
s
0
(
R
z
2
0
2
R
z
2
e
t
+
R
z
1
0
2
R
z
1
e
t
)
=
k
z
s
0
(
R
z
2
0
2
R
z
2
e
t
+
R
z
1
0
2
R
z
1
e
t
)
.


### 2.4. Dynamic Current Stiffness Model for the RMB Unit of CRAMB

The RMB unit controls the rotor movement alongx- or (/and) y-axis orthogonal to the spin axis (z-axis). When the rotor deviates radially from its “suspended” position at x-axis (the magnetomotive force at y-axis,  N
r
i
y
=
0), the net restoring force of the RMB unit at x-axis, f
x
i, can be given as
(40)
f
x
i
=
ϕ
x
1
2
2
μ
0
A
r
-
ϕ
x
2
2
2
μ
0
A
r
,
where,
(41)
ϕ
x
1
=
ϕ
r
m
+
ϕ
x
1
i
,
ϕ
x
2
=
ϕ
r
m
-
ϕ
x
2
i
,
ϕ
r
m
=
ϕ
m
0
4
.The flux in the air gap of the RMB unit produced by thex-axis coil current, i
x, is calculated as
(42)
ϕ
x
1
i
=
N
r
i
x
R
x
1
+
1
/
(
1
/
R
x
2
+
1
/
R
y
1
+
1
/
R
y
2
)
+
N
r
i
x
R
x
2
+
1
/
(
1
/
R
x
1
+
1
/
R
y
1
+
1
/
R
y
2
)
×
1
1
/
R
x
1
+
1
/
R
y
1
+
1
/
R
y
2
1
R
x
1
,
ϕ
x
2
i
=
N
r
i
x
R
x
2
+
1
/
(
1
/
R
x
1
+
1
/
R
y
1
+
1
/
R
y
2
)
+
N
r
i
x
R
x
1
+
1
/
(
1
/
R
x
2
+
1
/
R
y
1
+
1
/
R
y
2
)
×
1
1
/
R
x
2
+
1
/
R
y
1
+
1
/
R
y
2
1
R
x
2
.The nominal magnetic reluctances of the air gap in RMB unit are equal to each other when the rotor locates its center position,R
x
=
R
x
1
=
R
x
2
=
R
y
=
R
y
1
=
R
y
2. Then (42) can be modified as
(43)
ϕ
x
i
=
ϕ
x
2
i
=
ϕ
x
2
i
=
N
r
i
x
R
x
.Substituting (43) for the right side of (40), (40) can be modified as
(44)
f
x
i
=
2
ϕ
r
m
μ
0
A
r
N
r
i
x
R
x
.The eddy-current reluctance can be ignored due to laminated rotor and stator core, and the controls current do not pass through the other nonlaminated materials. So the dynamic current stiffness of the RMB unit will be equal to its static stiffness. Consider the following:(45)
k
r
i
i
=
k
r
i
0
=
d
f
x
i
d
i
x
=
2
ϕ
r
m
N
r
μ
0
A
r
R
x
=
N
r
2
μ
0
A
r
R
x
F
m
R
2
0
+
R
13
/
4
+
R
z
1
0
/
2
.


### 2.5. Dynamic Displacement Stiffness Model for the RMB Unit of CRAMB

When the coil currents in the RMB and the AMB unit are zero, the air gap in the AMB unit remains stationary, but the air gap in the RMB unit varies sinusoidally about the nominal values
0 in x-axis; that is, x
1
=
s
0
+
δ
sin
(
ω
t
) in upper air gap; and x
2
=
s
0
-
δ
sin
(
ω
t
) in lower air gap, then the magnetic reluctance of the air gap will also vary sinusoidally about the nominal value. Thus, the reluctance of the two air gaps in the RMB unit includes two parts, one is the static reluctance, R
x
0, and the other one is the dynamic reluctance, R
x
d. Consider the following:
(46)
R
x
1
=
R
x
0
+
R
x
d
=
s
0
μ
0
A
r
+
δ
sin
(
ω
t
)
μ
0
A
r
,
R
x
2
=
R
x
0
-
R
x
d
=
s
0
μ
0
A
r
-
δ
sin
(
ω
t
)
μ
0
A
r
,
where
(47)
R
x
0
=
s
0
μ
0
A
r
,
R
x
d
=
δ
sin
(
ω
t
)
μ
0
A
r
.If all the leakage reluctances and eddy-current reluctance are ignored, the total reluctanceR
met in the bias magnetic circuit is
(48)
R
met
=
∑
j
=
1
2
R
j
0
+
R
7
0
+
∑
k
=
11
12
R
k
0
+
1
2
/
R
y
+
1
/
R
x
1
+
1
/
R
x
2
+
1
1
/
∑
k
=
3
6
R
k
0
+
1
/
∑
j
=
8
10
R
j
0
.The relative permeability of the stator and rotor core is lager compared to the air gap, then their static reluctance can be ignored in the equivalent magnetic circuit of the bias flux path for calculated the dynamic displacement stiffness. The total reluctance in (48) can be simplified as
(49)
R
met
=
R
met
0
+
R
met
e
≈
R
2
0
+
1
2
/
R
y
+
1
/
R
x
1
+
1
/
R
x
2
+
R
z
1
0
2
,
where,
(50)
R
met
0
=
R
2
0
+
R
x
0
4
+
R
z
1
0
2
,
R
met
e
=
1
2
/
R
x
0
+
1
/
R
x
1
+
1
/
R
x
2
-
R
x
0
4
.Since dynamic magnetic reluctance is smaller compared to static magnetic reluctance of the air gap,R
met
e can be expanded into the power series and approximated by the first two terms as the higher order terms can be ignored, and the R
met
e will approximately equal to zero, R
met
e
≈
0.The flux in the upper and lower air gap of RMB inx-axis is calculated, respectively. Consider the following:
(51)
ϕ
x
1
=
F
m
R
met
1
2
/
R
x
0
+
1
/
R
x
1
+
1
/
R
x
2
1
R
x
1
,
ϕ
x
2
=
F
m
R
met
1
2
/
R
x
0
+
1
/
R
x
1
+
1
/
R
x
2
1
R
x
2
.Substituting (46) for the right side of (51), since R
x
d is smaller compared to R
x
0, (51) can be expanded into the power series and approximated by the first two terms as the higher order terms can be ignored. Consider the following:
(52)
ϕ
x
1
=
F
m
4
R
met
(
1
-
R
x
d
R
x
0
)
=
ϕ
x
0
(
1
-
R
x
d
R
x
0
)
ϕ
x
2
=
F
m
4
R
met
(
1
+
R
x
d
R
x
0
)
=
ϕ
x
0
(
1
+
R
x
d
R
x
0
)
,
where ϕ
x
0
=
F
m
/
4
R
met.The net force produced by RMB unit inx-axis is given by
(53)
f
x
s
=
(
ϕ
x
1
)
2
2
μ
0
A
r
-
(
ϕ
x
2
)
2
2
μ
0
A
r
=
ϕ
x
0
2
2
μ
0
A
r
[
(
1
-
R
x
d
R
x
0
)
2
-
(
1
+
R
x
d
R
x
0
)
2
]
=
-
2
ϕ
x
0
2
(
μ
0
A
r
)
2
δ
sin
(
ω
t
)
R
x
0
.The static displacement stiffness is calculated by(54)
k
r
s
0
=
d
f
x
s
d
(
δ
sin
(
ω
t
)
)
=
-
2
ϕ
x
0
2
μ
0
A
r
s
0
.In the static analysis model mentioned above, the eddy-current magnetic reluctance is not considered. However, the dynamic field should also cause eddy current, and the static model could be considered to include eddy-current magnetic reluctance as in the dynamic model. Then the dynamic flux in air gap of RMB unit inx-axis including eddy-current effect is given as
(55)
ϕ
x
1
=
F
m
R
sum
e
t
1
2
/
R
x
0
+
1
/
R
x
1
+
1
/
R
x
2
1
R
x
1
,
ϕ
x
2
=
F
m
R
sum
e
t
1
2
/
R
x
0
+
1
/
R
x
1
+
1
/
R
x
2
1
R
x
2
,
where R
sum
e
t
=
R
sum
0
+
R
sum
ω, R
sum
0, and R
sum
ω are given as
(56)
R
sum
0
=
∑
m
=
1
2
R
m
0
+
R
7
0
+
∑
n
=
11
12
R
n
0
+
1
2
/
R
x
0
+
1
/
R
x
1
+
1
/
R
x
2
+
R
z
1
0
2
,
R
sum
ω
=
∑
m
=
1
2
R
m
i
+
R
7
i
+
∑
n
=
11
12
R
n
i
+
1
1
/
∑
j
=
3
6
R
j
+
1
/
∑
k
=
8
10
R
k
-
R
z
1
0
2
.The permeability of the magnetic ring, stator core, rotor stator, and the thrust disk is compared to that of the air and the PMB, so (56) can be modified as
(57)
R
sum
0
=
R
2
0
+
1
2
/
R
x
0
+
1
/
R
x
1
+
1
/
R
x
2
+
R
z
1
0
2
,
R
sum
ω
=
∑
m
=
1
2
R
m
i
+
R
7
i
+
∑
n
=
11
12
R
n
i
+
1
1
/
∑
j
=
3
6
R
j
+
1
/
∑
k
=
8
10
R
k
-
R
z
1
0
2
.And (55) can be expanded into the power series and approximated by the first two terms as the higher order terms can be ignored.                                (58)
ϕ
x
1
=
F
m
R
sum
e
t
(
1
-
R
x
d
R
x
0
)
=
ϕ
x
0
R
met
R
sum
e
t
(
1
-
R
x
d
R
x
0
)
,
ϕ
x
2
=
F
m
R
sum
e
t
(
1
+
R
x
d
R
x
0
)
=
ϕ
x
0
R
met
R
sum
e
t
(
1
+
R
x
d
R
x
0
)
.The net force produced by the RMB unit considering the eddy-current magnetic reluctance is modified by(59)
f
x
s
d
=
(
ϕ
x
1
)
2
2
μ
0
A
r
-
(
ϕ
x
2
)
2
2
μ
0
A
r
.Substituting (58) for the right side of (59), (59) is modified as
(60)
f
x
s
d
=
(
ϕ
x
1
)
2
2
μ
0
A
r
-
(
ϕ
x
2
)
2
2
μ
0
A
r
=
-
2
ϕ
x
0
2
(
μ
0
A
r
)
2
(
R
met
R
sum
e
t
)
2
δ
sin
(
ω
t
)
R
x
0
.The dynamic displacement stiffness of RMB unit considering the eddy-current effect is given by(61)
k
r
s
d
=
d
f
x
s
d
(
δ
sin
(
ω
t
)
)
=
-
2
ϕ
x
0
2
μ
0
A
z
s
0
(
R
met
R
sum
e
t
)
2
=
k
r
s
0
(
R
met
R
sum
e
t
)
2
.


## 2.1. Analysis Model

Figure1 shows the configuration of the CRAMB which combines the radial and axial magnetic bearing unit. The axially magnetized permanent magnet ring (PMR) is used to produce bias flux for both the RMB and the AMB unit. The laminations are used only for radial stator and rotor core, and the other parts are made of solid material. Smaller outside diameter of the rotor thrust disk is used. Therefor the lower rotational drag, the lower stresses, and the compact structure are attained. The magnetic ring is used to eliminate the saturation of the magnetic circuit for the stator core of the RMB unit. The bias flux and the control flux paths of the RMB and AMB unit are depicted in Figure 2. The solid line indicates the bias flux paths generated by the PMR and the dotted line denotes the control flux paths generated by the control current of the AMB and the RMB unit. As can be seen, their control flux produced by the control current in coils will pass through air gap of the RAM and the AMB rather than the PMR, and therefore small control currents are enough since they do not have to compensate for the high magnetic reluctance of the PMR. The control flux has to be superimposed by the bias flux in air gap for producing the axial and radial magnetic force.Figure 1
Configuration of the CRAMB with permanent magnet ring.Figure 2
The magnetic circuits of the combined radial-axial magnetic bearing with permanent magnet ring.The radial and axial magnetic bearing force and flux can be derived conveniently with the static magnetic circuit model without consideration of the eddy-current effect. However, the eddy currents caused by the time-varying air gap and the control currents will affect the dynamic characteristics of the CRAMB system. Based on an approximate analysis method [22] considering the eddy-current effect inside a nonlaminated electromagnetic actuator, the CRAMB geometry is divided into thirteen elements as shown in Figure 3. The four elements of the air gap between rotor and stator of the RMB are named as element 13. The eddy-current loss can be reduced by using laminated stator core and rotor core of the RMB unit in CRAMB. According to the equivalent magnetic circuit models of the CRAMB, the effective reluctances are separated into two parts: one is the static magnetic reluctance of element, and the other one is the dynamic magnetic reluctance which is a half-order term related to frequency.Figure 3
Division of the geometry of the CRAMB: 13 elements simply divided.The elements 1–5, 7, and 9–12 can be divided in the same way as in Sun et al.’s models [30] and Zhu et al.’s models [23]. The PMR is used for element 2 which is used to provide bias flux and can also be simplified in the same way. The elements 1, 2, 11, and 12 are located in the bias flux path. Therefore, the eddy-current effect should be considered in mathematical model. The simplified effective reluctance of each element can be used to design, analyze, and control purpose. The simplified calculation of the effective magnetic reluctances for these elements is given as follow:
(1)
R
1
=
R
1
0
+
R
1
i
=
h
1
π
μ
0
μ
r
(
r
6
2
-
r
5
2
)
+
h
1
2
π
r
5
j
ω
σ
1
μ
0
μ
r
,
where R
1
0
=
h
1
/
π
μ
0
μ
r
(
r
6
2
-
r
5
2
)
,
R
1
i
=
(
h
1
/
2
π
r
5
)
j
ω
σ
1
/
μ
0
μ
r
(2)
R
2
=
R
2
0
+
R
2
i
=
h
2
π
μ
0
μ
r
(
r
6
2
-
r
5
2
)
+
h
2
2
π
r
5
j
ω
σ
2
μ
0
μ
r
,
where R
2
0
=
h
2
/
π
μ
0
μ
r
(
r
6
2
-
r
5
2
)
,
R
2
i
=
(
h
2
/
2
π
r
5
)
j
ω
σ
2
/
μ
0
μ
r
(3)
R
3
=
R
3
0
+
R
3
i
=
h
4
π
μ
0
μ
r
(
r
6
2
-
r
5
2
)
+
h
4
2
π
r
5
j
ω
σ
1
μ
0
μ
r
,
where R
3
0
=
h
4
/
π
μ
0
μ
r
(
r
6
2
-
r
5
2
),  R
3
i
=
(
h
4
/
2
π
r
5
)
j
ω
σ
1
/
μ
0
μ
r
(4)
R
4
=
R
4
0
+
R
4
i
=
ln
(
r
6
/
r
4
)
2
π
μ
0
μ
r
h
5
+
ln
(
r
6
/
r
4
)
2
π
j
ω
σ
1
μ
0
μ
r
,
where R
4
0
=
ln
(
r
6
/
r
4
)
/
2
π
μ
0
μ
r
h
5,  R
4
i
=
(
ln
(
r
6
/
r
4
)
/
2
π
)
j
ω
σ
1
/
μ
0
μ
r
(5)
R
5
=
R
5
0
+
R
5
i
=
h
6
π
μ
0
μ
r
(
r
4
2
-
r
3
2
)
+
h
6
2
π
r
3
j
ω
σ
1
μ
0
μ
r
,
where R
5
0
=
h
6
/
π
μ
0
μ
r
(
r
4
2
-
r
3
2
), R
5
i
=
(
h
6
/
2
π
r
3
)
j
ω
σ
1
/
μ
0
μ
r
(6)
R
z
1
=
R
6
=
s
0
+
z
π
μ
0
(
r
4
2
-
r
3
2
)
,
R
z
10
=
R
6
0
=
s
0
π
μ
0
(
r
4
2
-
r
3
2
)
,
z
=
0
,
R
7
=
R
7
0
+
R
7
i
=
ln
(
r
4
/
r
3
)
2
π
μ
0
μ
r
h
7
+
ln
(
r
4
/
r
3
)
2
π
j
ω
σ
1
μ
0
μ
r
,
where R
7
0
=
ln
(
r
4
/
r
3
)
/
2
π
μ
0
μ
r
h
7,  R
7
i
=
(
ln
(
r
4
/
r
3
)
/
2
π
)
j
ω
σ
1
/
μ
0
μ
r
(7)
R
z
2
=
R
8
=
s
0
-
z
π
μ
0
(
r
4
2
-
r
3
2
)
,
R
z
0
=
R
z
20
=
R
8
0
=
R
z
10
=
s
0
π
μ
0
(
r
4
2
-
r
3
2
)
,
z
=
0
,
R
9
=
R
9
0
+
R
9
i
=
R
5
=
R
5
0
+
R
5
i
,
where R
9
0
=
h
6
/
π
μ
0
μ
r
(
r
4
2
-
r
3
2
), R
9
i
=
(
h
6
/
2
π
r
3
)
j
ω
σ
1
/
μ
0
μ
r
(8)
R
10
=
R
10
0
+
R
10
i
=
R
4
=
R
4
0
+
R
4
i
,
where R
10
0
=
ln
(
r
6
/
r
4
)
/
2
π
μ
0
μ
r
h
5, R
10
i
=
(
ln
(
r
6
/
r
4
)
/
2
π
)
j
ω
σ
1
/
μ
0
μ
r
(9)
R
11
=
R
11
0
+
R
11
i
=
ln
(
r
3
/
r
1
)
2
π
μ
0
μ
r
h
7
+
ln
(
r
3
/
r
1
)
2
π
j
ω
σ
1
μ
0
μ
r
,
where R
11
0
=
ln
(
r
3
/
r
1
)
/
2
π
μ
0
μ
r
h
7, R
11
i
=
(
ln
(
r
3
/
r
1
)
/
2
π
)
j
ω
σ
1
/
μ
0
μ
r
(10)
R
12
=
R
12
0
+
R
12
i
=
h
9
π
μ
0
μ
r
(
r
1
2
-
r
8
2
)
+
h
9
2
π
r
8
j
ω
σ
1
μ
0
μ
r
,
where R
12
0
=
h
9
/
π
μ
0
μ
r
(
r
1
2
-
r
8
2
), R
12
i
=
(
h
9
/
2
π
r
8
)
j
ω
σ
1
/
μ
0
μ
r
(11)
R
13
=
R
13
0
=
360
(
r
3
-
r
2
)
π
μ
0
h
10
(
r
3
+
r
2
)
,
where R
j
0 is the static magnetic reluctance (the first term in R
j for element j) and R
j
i is the eddy-current magnetic reluctance (the second term in R
j for element j) with j
=
1,2
,
…
,
1
3, and the parameters and their dimensions are given in Figure 4 and Table 1. R
13 is the magnetic reluctance of the radial air gap. R
6 and R
8 are the magnetic reluctances of the axial air gap. The equivalent magnetic circuits including the bias magnetic circuit and the control magnetic circuits of the CRAMB without the leakage flux are shown in Figure 5. We focus on the investigation of the eddy-current effects, and the equivalent magnetic reluctances for the leakage flux paths would be included in the future model.Table 1
Parameters of The CRAMB.
Parameters
Description
Value
r
1
Inner radius of the RMB rotor, (mm)
25
r
2
Outer radius of the RMB rotor, (mm)
30
r
7
Inner radius of the RMB stator, (mm)
30.5
A
r
Radial pole shoe area, (mm2)
1036
N
r
The number of winding turns in coil of RMB
220
r
9
Outer radius of the RMB stator, (mm)
60
r
6
Outer radius of the yoke, (mm)
70
r
3
Inner radius of the axial pole shoe, (mm)
31
r
4
Outer radius of the axial pole shoe, (mm)
34
r
5
Inner radius of the yoke, (mm)
62
r
8
Inner radius of the RMB rotor, (mm)
10
h
1
Length of the stator return ring of RMB, (mm)
31
h
2
Thickness of the permanent ring, (mm)
10
h
3
Thickness of the element 10, (mm)
8
h
4
Thickness of the element 3, (mm)
24
h
5
Thickness of the element 4, (mm)
8
h
6
Thickness of the element 5, (mm)
8
h
7
Thickness of the thrust plane, (mm)
7
h
8
Thickness of the element 9, (mm)
8
h
9
Thickness of the element 12, (mm)
8
s
0
radial and axial air gap, (mm)
0.5
N
z
The number of winding turns in coil of AMB
120
m
The mass of the rotor, (kg)
7.2
σ
1
The conductivity of yoke and thrust disk, (S/m)
8 × 106
μ
0
Permeability of free space, (H/m)
4
π
×
10
-
7
μ
r
Relative permeability of the permanent magnet
1.082
σ
2
Resistivity of the magnet ring, (Ω·m)
0.9
E
-
6Figure 4
The dimension of the CRAMB.The equivalent magnetic circuits of the CRAMB. (a) The equivalent magnetic circuit of the bias flux paths. (b) The equivalent magnetic circuit of control flux path of the RMB unit. (c) The equivalent magnetic circuit of control flux path of the AMB unit.
(a)
(b)
(c)


## 2.2. Dynamic Current Stiffness Model for the AMB Unit of CRAMB

It is assumed that the CRAMB system is linear, and the AMB unit of the CRAMB excited a sinusoidal varying ac currenti
z
=
I
z
sin
(
ω
t
) (where the amplitude, I
z, is very small). When the thrust disk located its center position, the length of air gap (the element 6, 8) is s
0.The control fluxϕ
z
i produced by a sinusoidal varying ac current i
z in the axial air gap of AMB unit and the bias flux ϕ
z
m
0 in the axial air gap of AMB unit produced by the PMR, the control flux in the left air gap (element 8) will add to the bias flux, and will subtract from the bias flux in the right air gap (element 6) as shown in Figure 3. The total flux, ϕ
z
sum
i
(
i
=
1
,
2
), in the left and the right air gap of AMB unit is given as
(12)
ϕ
z
sum
1
=
ϕ
z
m
0
+
ϕ
z
i
,
ϕ
z
sum
2
=
ϕ
z
m
0
-
ϕ
z
i
,
where the static bias flux in the axial air gap of AMB unit is calculated as
(13)
ϕ
z
m
=
F
m
(
1
1
/
∑
k
=
3
6
R
k
0
+
1
/
∑
j
=
8
10
R
j
0
R
1
0
+
R
2
0
+
R
7
0
+
R
11
0
+
R
12
0
+
R
13
4
+
1
1
/
∑
k
=
3
6
R
k
0
+
1
/
∑
j
=
8
10
R
j
0
)
-
1
.The permanent magnet material of the PMR is Nd-Fe-B, its magnetomotive force isF
m
=
H
c
h
2, and H
c is the coercive force of the PMR; h
2 is thickness of the PMR. The permeability of the magnetic ring and the thrust disk is compared to that of the air and the PMB, so (13) can be simplified as
(14)
ϕ
z
m
≈
F
m
R
2
0
+
R
13
/
4
+
R
z
0
/
2
.And the static bias flux in the air gap of the AMB can be calculated as(15)
ϕ
z
m
0
=
ϕ
z
1
0
=
ϕ
z
2
0
=
ϕ
z
m
2
=
F
m
2
R
2
0
+
R
13
/
2
+
R
z
0
.The dynamic flux in the air gap of AMB unit produced by the ac current is calculated as(16)
ϕ
z
i
=
N
z
I
z
sin
(
ω
t
)
∑
k
=
1
10
R
k
=
N
z
I
z
R
z
i
0
sin
(
ω
t
)
1
+
R
e
(
ω
)
/
R
z
i
0
=
N
z
I
z
sin
(
ω
t
)
R
z
i
0
(
1
+
λ
)
,
where R
z
i
0 is the total static magnetic reluctance of the equivalent magnetic circuit of control flux path for AMB unit, R
e
(
ω
) is the total eddy-current magnetic reluctance of the equivalent magnetic circuit of control flux path for AMB unit, and λ is the ratio of the total eddy-current magnetic reluctance to the total static magnetic reluctance R
e
(
ω
)/R
z
i
0:
(17)
R
z
i
0
=
R
6
0
+
R
8
0
,
R
e
(
ω
)
=
R
3
i
+
R
4
i
+
R
5
i
+
R
7
i
+
R
9
i
+
R
10
i
.R
e
(
ω
) is the magnetic reluctance including the eddy current effect, and it will be increased with the addition of frequency of the time-varying air gap and the control currents. The bandwidth will be increased with the addition of the air gap.The AMB unit of the CRAMB will produce a net restoring forcef
z
i on the thrust disk in z-axis direction. Since the net force, f
z
i, is calculated by
(18)
f
z
i
=
ϕ
z
1
2
2
μ
0
A
z
-
ϕ
z
2
2
2
μ
0
A
z
,
where the constant μ
0 is the permeability of free space, A
z is the area of pole face of the AMB unit; substituting (12) into (18) gives
(19)
f
z
i
=
(
ϕ
z
m
0
+
ϕ
z
i
)
2
2
μ
0
A
z
-
(
ϕ
z
m
0
-
ϕ
z
i
)
2
2
μ
0
A
z
=
2
ϕ
z
m
0
ϕ
z
i
μ
0
A
z
.The dynamic current stiffness of the AMB unit can be derived as(20)
k
z
i
i
=
d
f
z
i
d
i
z
=
2
ϕ
z
m
0
N
z
μ
0
A
z
R
z
i
0
(
1
+
λ
)
=
k
z
i
0
1
+
λ
,
where
(21)
k
z
i
0
=
2
ϕ
z
m
0
N
z
μ
0
A
z
R
z
i
0
is the static current stiffness of the AMB unit.


## 2.3. Dynamic Displacement Stiffness Model for the AMB Unit of CRAMB

When the coil currents in the RMB and the AMB unit are zero, the air gap in the RMB remains stationary, but the air gap in the AMB unit varies sinusoidally about the nominal values
0; that is, z
1
=
s
0
-
δ
sin
(
ω
t
) in left air gap, and z
2
=
s
0
+
δ
sin
(
ω
t
) in right air gap; then the magnetic reluctance of the air gap will also vary sinusoidally about the nominal value. Thus, the reluctance of the two air gaps in AMB unit includes two parts, one is the static reluctance, R
z
i
0
(
i
=
1
,
2
), and the other one is the dynamic reluctance, R
z
i
d
(
i
=
1,2
). Consider the following:
(22)
R
z
1
=
R
z
1
0
+
R
z
1
d
=
s
0
μ
0
A
z
+
δ
sin
(
ω
t
)
μ
0
A
z
,
R
z
2
=
R
z
2
0
+
R
z
2
d
=
s
0
μ
0
A
z
-
δ
sin
(
ω
t
)
μ
0
A
z
,
where
(23)
R
z
1
0
=
R
z
2
0
=
s
0
μ
0
A
z
,
R
z
1
d
=
-
R
z
2
d
=
δ
sin
(
ω
t
)
μ
0
A
z
.If all the eddy-current reluctances and leakage reluctances are ignored, the total static magnetic reluctanceR
met in the bias magnetic circuit is
(24)
R
met
=
R
1
0
+
R
2
0
+
R
7
0
+
R
11
0
+
R
12
0
+
R
13
4
+
1
1
/
∑
k
=
3
6
R
k
0
+
1
/
∑
j
=
8
10
R
j
0
.Since the relative permeability of the stator and rotor core of the AMB and the RMB unit is larger, their static magnetic reluctance can be ignored in the equivalent magnetic circuit of the bias flux path for calculating the dynamic displacement stiffness. The total static magnetic reluctance in (24) can be simplified as
(25)
R
met
=
R
met
0
+
R
met
e
≈
R
2
0
+
R
13
4
+
1
1
/
R
z
1
+
1
/
R
z
2
,
where,
(26)
R
met
0
=
R
2
0
+
R
13
4
+
R
z
1
0
R
z
2
0
R
z
1
0
+
R
z
2
0
,
R
met
e
=
(
R
z
2
0
-
R
z
1
0
)
R
z
1
d
-
(
R
z
1
d
)
2
R
z
1
0
+
R
z
2
0
.Sinceδ is smaller compared to s
0, then R
z
i
d
(
i
=
1
,
2
) is smaller compared to R
z
i
0
(
i
=
1
,
2
), and R
met
e in (26) can be expanded into the power series and approximated by the first two terms as the higher order terms can be ignored. Consider the following:
(27)
R
met
0
=
R
2
0
+
R
13
4
+
R
z
1
0
R
z
2
0
R
z
1
0
+
R
z
2
0
,
R
met
e
≈
(
R
z
2
0
-
R
z
1
0
)
R
z
1
d
R
z
1
0
+
R
z
2
0
.The flux in the air gap of AMB unit can be calculated as(28)
ϕ
z
1
=
F
m
R
met
0
+
R
met
e
1
1
/
(
R
z
1
0
+
R
z
1
d
)
+
1
/
(
R
z
2
0
-
R
z
1
d
)
1
R
z
1
0
+
R
z
1
d
,
ϕ
z
2
=
F
m
R
met
0
+
R
met
e
1
1
/
(
R
z
1
0
+
R
z
1
d
)
+
1
/
(
R
z
2
0
-
R
z
1
d
)
1
R
z
2
0
-
R
z
1
d
.Substituting (27) for the right side of (28), since R
z
i
d is smaller compared to R
z
i
0, (28) can be expanded into the power series and approximated by the first two terms as the higher order terms can be ignored. Consider the following:
(29)
ϕ
z
1
≈
F
m
R
met
0
R
z
1
0
R
z
1
0
+
R
z
2
0
[
1
-
R
z
1
d
R
met
0
(
R
z
2
0
-
R
z
1
0
R
z
1
0
+
R
z
2
0
)
-
R
z
1
d
R
z
2
0
]
,
ϕ
z
2
≈
F
m
R
met
0
R
z
2
0
R
z
1
0
+
R
z
2
0
[
1
-
R
z
1
d
R
met
0
(
R
z
2
0
-
R
z
1
0
R
z
1
0
+
R
z
2
0
)
+
R
z
1
d
R
z
1
0
]
.The Equation (29) can be modified due to R
z
1
0
=
R
z
2
0.                              (30)
ϕ
z
1
≈
F
m
2
R
met
0
(
1
-
A
R
z
1
d
R
met
0
-
R
z
1
d
R
z
2
0
)
=
ϕ
z
0
(
1
-
A
R
z
1
d
R
met
0
-
R
z
1
d
R
z
2
0
)
,
ϕ
z
2
≈
F
m
2
R
met
0
(
1
-
A
R
z
1
d
R
met
0
+
R
z
1
d
R
z
1
0
)
=
ϕ
z
0
(
1
-
A
R
z
1
d
R
met
0
+
R
z
1
d
R
z
1
0
)
,
where A
=
(
R
z
2
e
t
-
R
z
1
e
t
)
/
(
R
z
1
0
+
R
z
2
0
).The AMB unit of the CRAMB will produce a net restoring forcef
z
i on the thrust disk in z-axis direction, since the force is given by
(31)
f
z
s
=
(
ϕ
z
1
)
2
2
μ
0
A
z
-
(
ϕ
z
2
)
2
2
μ
0
A
z
=
ϕ
z
0
2
2
μ
0
A
z
×
[
(
1
-
A
R
z
1
d
R
met
0
-
R
z
1
d
R
z
2
0
)
2
-
(
1
-
A
R
z
1
d
R
met
0
+
R
z
1
d
R
z
1
0
)
2
]
,
where ϕ
z
0
=
ϕ
z
1
0
=
ϕ
z
2
0
=
F
m
/
2
R
met
0 is the static bias flux in the air gap of the AMB unit.SinceR
z
i
d
(
i
=
1
,
2
) is smaller compared to R
z
i
0
(
i
=
1
,
2
), (31) can be expanded into the power series and approximated by the first two terms as the higher order terms can be ignored. Consider the following:
(32)
f
z
s
≈
-
2
ϕ
z
0
2
(
μ
0
A
z
)
2
δ
sin
(
ω
t
)
R
z
1
0
.Static displacement stiffness of the AMB unit is given as(33)
k
z
s
0
=
d
f
z
s
d
(
δ
sin
(
ω
t
)
)
=
-
2
ϕ
z
0
2
μ
0
A
z
s
0
=
-
1
2
μ
0
A
z
s
0
(
F
m
R
2
0
+
R
13
/
4
+
R
z
1
0
/
2
)
2
.In the static analysis model mentioned above, the eddy-current magnetic reluctance is not considered. However, the dynamic field should also cause eddy current, and the static model could be considered to include eddy-current magnetic reluctance as in the dynamic model. Then the dynamic flux including eddy-current effect is given as(34)
ϕ
z
1
d
=
ϕ
z
1
0
(
1
-
B
R
z
1
d
R
sum
e
t
-
R
z
1
d
R
z
2
e
t
)
,
ϕ
z
2
d
=
ϕ
z
2
0
(
1
-
B
R
z
1
d
R
sum
e
t
+
R
z
1
d
R
z
1
e
t
)
,
where,
(35)
R
sum
e
t
=
R
sum
0
+
R
sum
ω
,
R
sum
0
=
R
1
0
+
R
2
0
+
R
7
0
+
R
11
0
+
R
12
0
+
R
13
4
+
R
z
1
0
2
,
R
sum
ω
=
∑
m
=
1
2
R
m
i
+
R
7
i
+
∑
n
=
11
12
R
n
i
+
1
1
/
∑
j
=
3
6
R
j
+
1
/
∑
k
=
8
10
R
k
-
R
z
1
0
2
,
B
=
R
z
2
e
t
-
R
z
1
e
t
R
z
1
e
t
+
R
z
2
e
t
,
R
z
1
e
t
=
∑
j
=
3
6
R
j
,
R
z
2
e
t
=
∑
k
=
8
10
R
k
.The permeability of the magnetic ring, stator core, rotor stator, and the thrust disk is compared to that of the air and the PMR, so (35) can be modified as
(36)
R
sum
0
≈
R
2
0
+
R
13
4
+
R
z
1
0
2
,
R
sum
ω
≈
∑
m
=
1
2
R
m
i
+
R
7
i
+
∑
n
=
11
12
R
n
i
+
1
1
/
∑
j
=
3
6
R
j
i
+
1
/
∑
k
=
8
10
R
k
i
-
R
z
1
0
2
,
R
z
1
e
t
=
R
3
i
+
R
4
i
+
R
5
i
+
R
6
0
,
R
z
2
e
t
=
R
8
0
+
R
9
i
+
R
10
i
.The net force of the AMB unit considering the eddy-current magnetic reluctance is given by(37)
f
z
s
=
(
ϕ
z
1
d
)
2
2
μ
0
A
z
-
(
ϕ
z
2
d
)
2
2
μ
0
A
z
=
ϕ
z
0
2
2
μ
0
A
z
[
(
1
-
B
R
z
1
d
R
sum
e
t
-
R
z
1
d
R
z
2
e
t
)
2

-
(
1
-
B
R
z
1
d
R
sum
e
t
+
R
z
1
d
R
z
1
e
t
)
2
]

.SinceR
z
i
d is smaller compared to R
z
i
0, (37) can be expanded into the power series and approximated by the first two terms as the higher order terms can be ignored. Consider the following:
(38)
f
z
s
d
≈
-
2
ϕ
z
0
2
(
μ
0
A
z
)
2
(
1
2
R
z
2
e
t
+
1
2
R
z
1
e
t
)
δ
sin
(
ω
t
)
.The dynamic displacement stiffness considering the eddy-current magnetic reluctance is calculated as(39)
k
z
s
d
=
d
f
z
s
d
d
(
δ
sin
(
ω
t
)
)
=
-
2
ϕ
z
0
2
μ
0
A
z
s
0
(
R
z
2
0
2
R
z
2
e
t
+
R
z
1
0
2
R
z
1
e
t
)
=
k
z
s
0
(
R
z
2
0
2
R
z
2
e
t
+
R
z
1
0
2
R
z
1
e
t
)
.


## 2.4. Dynamic Current Stiffness Model for the RMB Unit of CRAMB

The RMB unit controls the rotor movement alongx- or (/and) y-axis orthogonal to the spin axis (z-axis). When the rotor deviates radially from its “suspended” position at x-axis (the magnetomotive force at y-axis,  N
r
i
y
=
0), the net restoring force of the RMB unit at x-axis, f
x
i, can be given as
(40)
f
x
i
=
ϕ
x
1
2
2
μ
0
A
r
-
ϕ
x
2
2
2
μ
0
A
r
,
where,
(41)
ϕ
x
1
=
ϕ
r
m
+
ϕ
x
1
i
,
ϕ
x
2
=
ϕ
r
m
-
ϕ
x
2
i
,
ϕ
r
m
=
ϕ
m
0
4
.The flux in the air gap of the RMB unit produced by thex-axis coil current, i
x, is calculated as
(42)
ϕ
x
1
i
=
N
r
i
x
R
x
1
+
1
/
(
1
/
R
x
2
+
1
/
R
y
1
+
1
/
R
y
2
)
+
N
r
i
x
R
x
2
+
1
/
(
1
/
R
x
1
+
1
/
R
y
1
+
1
/
R
y
2
)
×
1
1
/
R
x
1
+
1
/
R
y
1
+
1
/
R
y
2
1
R
x
1
,
ϕ
x
2
i
=
N
r
i
x
R
x
2
+
1
/
(
1
/
R
x
1
+
1
/
R
y
1
+
1
/
R
y
2
)
+
N
r
i
x
R
x
1
+
1
/
(
1
/
R
x
2
+
1
/
R
y
1
+
1
/
R
y
2
)
×
1
1
/
R
x
2
+
1
/
R
y
1
+
1
/
R
y
2
1
R
x
2
.The nominal magnetic reluctances of the air gap in RMB unit are equal to each other when the rotor locates its center position,R
x
=
R
x
1
=
R
x
2
=
R
y
=
R
y
1
=
R
y
2. Then (42) can be modified as
(43)
ϕ
x
i
=
ϕ
x
2
i
=
ϕ
x
2
i
=
N
r
i
x
R
x
.Substituting (43) for the right side of (40), (40) can be modified as
(44)
f
x
i
=
2
ϕ
r
m
μ
0
A
r
N
r
i
x
R
x
.The eddy-current reluctance can be ignored due to laminated rotor and stator core, and the controls current do not pass through the other nonlaminated materials. So the dynamic current stiffness of the RMB unit will be equal to its static stiffness. Consider the following:(45)
k
r
i
i
=
k
r
i
0
=
d
f
x
i
d
i
x
=
2
ϕ
r
m
N
r
μ
0
A
r
R
x
=
N
r
2
μ
0
A
r
R
x
F
m
R
2
0
+
R
13
/
4
+
R
z
1
0
/
2
.


## 2.5. Dynamic Displacement Stiffness Model for the RMB Unit of CRAMB

When the coil currents in the RMB and the AMB unit are zero, the air gap in the AMB unit remains stationary, but the air gap in the RMB unit varies sinusoidally about the nominal values
0 in x-axis; that is, x
1
=
s
0
+
δ
sin
(
ω
t
) in upper air gap; and x
2
=
s
0
-
δ
sin
(
ω
t
) in lower air gap, then the magnetic reluctance of the air gap will also vary sinusoidally about the nominal value. Thus, the reluctance of the two air gaps in the RMB unit includes two parts, one is the static reluctance, R
x
0, and the other one is the dynamic reluctance, R
x
d. Consider the following:
(46)
R
x
1
=
R
x
0
+
R
x
d
=
s
0
μ
0
A
r
+
δ
sin
(
ω
t
)
μ
0
A
r
,
R
x
2
=
R
x
0
-
R
x
d
=
s
0
μ
0
A
r
-
δ
sin
(
ω
t
)
μ
0
A
r
,
where
(47)
R
x
0
=
s
0
μ
0
A
r
,
R
x
d
=
δ
sin
(
ω
t
)
μ
0
A
r
.If all the leakage reluctances and eddy-current reluctance are ignored, the total reluctanceR
met in the bias magnetic circuit is
(48)
R
met
=
∑
j
=
1
2
R
j
0
+
R
7
0
+
∑
k
=
11
12
R
k
0
+
1
2
/
R
y
+
1
/
R
x
1
+
1
/
R
x
2
+
1
1
/
∑
k
=
3
6
R
k
0
+
1
/
∑
j
=
8
10
R
j
0
.The relative permeability of the stator and rotor core is lager compared to the air gap, then their static reluctance can be ignored in the equivalent magnetic circuit of the bias flux path for calculated the dynamic displacement stiffness. The total reluctance in (48) can be simplified as
(49)
R
met
=
R
met
0
+
R
met
e
≈
R
2
0
+
1
2
/
R
y
+
1
/
R
x
1
+
1
/
R
x
2
+
R
z
1
0
2
,
where,
(50)
R
met
0
=
R
2
0
+
R
x
0
4
+
R
z
1
0
2
,
R
met
e
=
1
2
/
R
x
0
+
1
/
R
x
1
+
1
/
R
x
2
-
R
x
0
4
.Since dynamic magnetic reluctance is smaller compared to static magnetic reluctance of the air gap,R
met
e can be expanded into the power series and approximated by the first two terms as the higher order terms can be ignored, and the R
met
e will approximately equal to zero, R
met
e
≈
0.The flux in the upper and lower air gap of RMB inx-axis is calculated, respectively. Consider the following:
(51)
ϕ
x
1
=
F
m
R
met
1
2
/
R
x
0
+
1
/
R
x
1
+
1
/
R
x
2
1
R
x
1
,
ϕ
x
2
=
F
m
R
met
1
2
/
R
x
0
+
1
/
R
x
1
+
1
/
R
x
2
1
R
x
2
.Substituting (46) for the right side of (51), since R
x
d is smaller compared to R
x
0, (51) can be expanded into the power series and approximated by the first two terms as the higher order terms can be ignored. Consider the following:
(52)
ϕ
x
1
=
F
m
4
R
met
(
1
-
R
x
d
R
x
0
)
=
ϕ
x
0
(
1
-
R
x
d
R
x
0
)
ϕ
x
2
=
F
m
4
R
met
(
1
+
R
x
d
R
x
0
)
=
ϕ
x
0
(
1
+
R
x
d
R
x
0
)
,
where ϕ
x
0
=
F
m
/
4
R
met.The net force produced by RMB unit inx-axis is given by
(53)
f
x
s
=
(
ϕ
x
1
)
2
2
μ
0
A
r
-
(
ϕ
x
2
)
2
2
μ
0
A
r
=
ϕ
x
0
2
2
μ
0
A
r
[
(
1
-
R
x
d
R
x
0
)
2
-
(
1
+
R
x
d
R
x
0
)
2
]
=
-
2
ϕ
x
0
2
(
μ
0
A
r
)
2
δ
sin
(
ω
t
)
R
x
0
.The static displacement stiffness is calculated by(54)
k
r
s
0
=
d
f
x
s
d
(
δ
sin
(
ω
t
)
)
=
-
2
ϕ
x
0
2
μ
0
A
r
s
0
.In the static analysis model mentioned above, the eddy-current magnetic reluctance is not considered. However, the dynamic field should also cause eddy current, and the static model could be considered to include eddy-current magnetic reluctance as in the dynamic model. Then the dynamic flux in air gap of RMB unit inx-axis including eddy-current effect is given as
(55)
ϕ
x
1
=
F
m
R
sum
e
t
1
2
/
R
x
0
+
1
/
R
x
1
+
1
/
R
x
2
1
R
x
1
,
ϕ
x
2
=
F
m
R
sum
e
t
1
2
/
R
x
0
+
1
/
R
x
1
+
1
/
R
x
2
1
R
x
2
,
where R
sum
e
t
=
R
sum
0
+
R
sum
ω, R
sum
0, and R
sum
ω are given as
(56)
R
sum
0
=
∑
m
=
1
2
R
m
0
+
R
7
0
+
∑
n
=
11
12
R
n
0
+
1
2
/
R
x
0
+
1
/
R
x
1
+
1
/
R
x
2
+
R
z
1
0
2
,
R
sum
ω
=
∑
m
=
1
2
R
m
i
+
R
7
i
+
∑
n
=
11
12
R
n
i
+
1
1
/
∑
j
=
3
6
R
j
+
1
/
∑
k
=
8
10
R
k
-
R
z
1
0
2
.The permeability of the magnetic ring, stator core, rotor stator, and the thrust disk is compared to that of the air and the PMB, so (56) can be modified as
(57)
R
sum
0
=
R
2
0
+
1
2
/
R
x
0
+
1
/
R
x
1
+
1
/
R
x
2
+
R
z
1
0
2
,
R
sum
ω
=
∑
m
=
1
2
R
m
i
+
R
7
i
+
∑
n
=
11
12
R
n
i
+
1
1
/
∑
j
=
3
6
R
j
+
1
/
∑
k
=
8
10
R
k
-
R
z
1
0
2
.And (55) can be expanded into the power series and approximated by the first two terms as the higher order terms can be ignored.                                (58)
ϕ
x
1
=
F
m
R
sum
e
t
(
1
-
R
x
d
R
x
0
)
=
ϕ
x
0
R
met
R
sum
e
t
(
1
-
R
x
d
R
x
0
)
,
ϕ
x
2
=
F
m
R
sum
e
t
(
1
+
R
x
d
R
x
0
)
=
ϕ
x
0
R
met
R
sum
e
t
(
1
+
R
x
d
R
x
0
)
.The net force produced by the RMB unit considering the eddy-current magnetic reluctance is modified by(59)
f
x
s
d
=
(
ϕ
x
1
)
2
2
μ
0
A
r
-
(
ϕ
x
2
)
2
2
μ
0
A
r
.Substituting (58) for the right side of (59), (59) is modified as
(60)
f
x
s
d
=
(
ϕ
x
1
)
2
2
μ
0
A
r
-
(
ϕ
x
2
)
2
2
μ
0
A
r
=
-
2
ϕ
x
0
2
(
μ
0
A
r
)
2
(
R
met
R
sum
e
t
)
2
δ
sin
(
ω
t
)
R
x
0
.The dynamic displacement stiffness of RMB unit considering the eddy-current effect is given by(61)
k
r
s
d
=
d
f
x
s
d
(
δ
sin
(
ω
t
)
)
=
-
2
ϕ
x
0
2
μ
0
A
z
s
0
(
R
met
R
sum
e
t
)
2
=
k
r
s
0
(
R
met
R
sum
e
t
)
2
.


## 3. Coupling Performance Analysis of the Dynamic Stiffness Models of CRAMB

Since the equivalent magnetic circuits of control flux path of the AMB and RMB unit do not influence each other, the current stiffness of AMB (or RMB) unit will not be affected by the control current of the RMB (or AMB) unit. But the equivalent bias circuits of the AMB and the RMB pass through the PMB, and the dynamic models of the AMB and the RMB unit may influence each other due the air gap varying of the AMB or RMB unit.


### 3.1. The Dynamic Current Stiffness of the AMB Unit Affected by the Rotor Position of RMB Unit inx-Axis

Based on the small variation of the rotor position of the RMB unit inx-axis, x
1
=
s
0
+
x, x
2
=
s
0
-
x, and the remains stationary in y-axis (y
1
=
y
2
=
s
0).The static bias flux in the air gap of the AMB unit is modified as(62)
ϕ
z
m
0
r
=
ϕ
z
1
r
0
=
ϕ
z
2
r
0
=
ϕ
m
0
r
2
=
F
m
R
2
0
+
1
/
(
2
/
R
x
0
+
1
/
R
x
1
+
1
/
R
x
2
)
+
R
z
1
0
/
2
=
F
m
R
2
0
+
1
/
(
2
μ
0
A
r
/
s
0
+
2
μ
0
A
r
s
0
/
(
s
0
2
-
x
2
)
)
+
R
z
1
0
/
2
.And the dynamic current stiffness of the AMB unit is modified due to rotor position varying inx-axis. Consider the following:
(63)
k
z
i
r
i
=
2
N
z
μ
0
A
z
R
z
i
0
(
1
+
λ
)
×
F
m
R
2
0
+
1
/
(
2
μ
0
A
r
/
s
0
+
2
μ
0
A
r
s
0
/
(
s
0
2
-
x
2
)
)
+
R
z
1
0
/
2
=
k
z
i
i
R
2
0
+
R
13
/
4
+
R
z
1
0
/
2
R
2
0
+
1
/
(
2
μ
0
A
r
/
s
0
+
2
μ
0
A
r
s
0
/
(
s
0
2
-
x
2
)
)
+
R
z
1
0
/
2
=
k
z
i
i
λ
r
z
,
where λ
r
z
=
(
R
2
0
+
R
13
/
4
+
R
z
1
0
/
2
)
/
(
R
2
0
+
1
/
(
2
μ
0
A
r
/
s
0
+
2
μ
0
A
r
s
0
/
(
s
0
2
-
x
2
)
)
+
R
z
1
0
/
2
) is the coefficient which reflects the dynamic current stiffness dependence of the AMB unit on the rotor position in x-axis (or y-axis). In other words, the coefficient, λ
r
z, reflects the coupling extent of the dynamic current stiffness of the AMB unit affected by the rotor position in x-axis or y-axis.From (63), the x term is included in the dynamic current stiffness model which will be affected by the rotor position of the RMB unit in x-axis or y-axis.


### 3.2. The Dynamic Displacement Stiffness of the AMB Unit Affected by the Rotor Position of RMB Unit inx-Axis

TheR
sum
0 in (36) will be modified when considering the small variation of the rotor position of the RMB unit in x-axis. Consider the following:
(64)
R
sum
0
=
R
2
0
+
1
2
μ
0
A
r
/
s
0
+
2
μ
0
A
r
s
0
/
(
s
0
2
-
x
2
)
+
R
z
1
0
R
z
2
0
R
z
1
0
+
R
z
2
0
.And the dynamic displacement stiffness is modified due to rotor position varying inx-axis. Consider the following:
(65)
k
z
s
r
d
=
-
1
2
μ
0
A
z
s
0
×
(
F
m
R
2
0
+
1
/
(
2
μ
0
A
r
/
s
0
+
2
μ
0
A
r
s
0
/
(
s
0
2
-
x
2
)
)
+
R
z
1
0
/
2
)
2
×
(
R
z
1
0
2
R
z
2
e
t
+
R
z
2
0
2
R
z
1
e
t
)
=
k
z
s
d
(
R
2
0
+
R
13
/
4
+
R
z
1
0
/
2
)
2
(
R
2
0
+
1
/
(
2
μ
0
A
r
/
s
0
+
2
μ
0
A
r
s
0
/
(
s
0
2
-
x
2
)
)
+
R
z
1
0
/
2
)
2
=
k
z
s
d
λ
r
z
2
.From (65), the x term is included in the dynamic displacement stiffness model which will be affected by the rotor position of the RMB unit in x-axis or y-axis. And λ
r
z
2 is the coefficient which reflects the dynamic displacement stiffness of the AMB unit dependence on the rotor position in x-axis (or y-axis). In other words, the coefficient, λ
r
z
2, reflects the coupling extent of the dynamic displacement stiffness of the AMB unit affected by the rotor position in x-axis or y-axis.


### 3.3. The Dynamic Current Stiffness of the RMB Unit Affected by the Rotor Position of AMB Unit inz-Axis

Based on the small variation of the rotor position of the AMB unit inz-axis, z
1
=
s
0
+
z, z
2
=
s
0
-
z, and the remains stationary in y-axis (R
y
1
=
R
y
2
≈
R
x
0). The magnetic reluctance of the air gap for the AMB unit can be modified as
(66)
R
z
1
=
s
0
+
z
μ
0
A
z
,
R
z
2
=
s
0
-
z
μ
0
A
z
.The static bias flux in the air gap of the RMB unit affected by the rotor position inz-axis varying is modified as
(67)
ϕ
r
m
=
ϕ
m
0
4
=
F
m
/
4
R
2
0
+
R
13
/
4
+
1
/
(
1
/
R
z
1
+
1
/
R
z
2
)
=
F
m
/
4
R
2
0
+
R
13
/
4
+
(
s
0
2
-
z
2
)
/
2
μ
0
A
z
s
0
.And the dynamic current stiffness of the RMB unit affected by the rotor position varying inz-axis is modified as
(68)
k
r
i
z
i
=
d
f
x
i
d
i
x
=
2
ϕ
r
m
N
r
μ
0
A
r
R
x
=
N
r
2
μ
0
A
r
R
x
F
m
R
2
0
+
R
x
/
4
+
(
s
0
2
-
z
2
)
/
2
μ
0
A
z
s
0
=
k
r
i
i
R
2
0
+
R
x
/
4
+
R
z
1
0
/
2
R
2
0
+
R
x
/
4
+
(
s
0
2
-
z
2
)
/
2
μ
0
A
z
s
0
=
k
r
i
i
λ
z
r
,
where λ
z
r
=
(
R
2
0
+
R
x
/
4
+
R
z
1
0
/
2
)
/
(
R
2
0
+
R
x
/
4
+
(
s
0
2
-
z
2
)
/
2
μ
0
A
z
s
0
).Base on (65), the z term is included in the dynamic current stiffness model which will be affected by the rotor position varying in z-axis.


### 3.4. The Dynamic Displacement Stiffness of the RMB Unit Affected by the Rotor Position of AMB Unit inz-Axis

The total static magnetic reluctance will be modified when considering the small variation of the rotor position of the AMB unit inz-axis. Consider the following:
(69)
R
met
=
R
2
0
+
R
x
0
4
+
1
1
/
R
z
1
0
+
1
/
R
z
2
0
=
R
2
0
+
R
x
0
4
+
s
0
2
-
z
2
2
μ
0
A
z
s
0
.The static bias flux in the air gap of the RMB unit is modified when considering the rotor position varying of the AMB unit inz-axis. Consider the following:
(70)
ϕ
x
0
=
F
m
4
R
met
=
F
m
4
(
R
2
0
+
R
x
0
/
4
+
(
s
0
2
-
z
2
)
/
2
μ
0
A
z
s
0
)
.The total dynamic magnetic reluctance affected by the eddy-current effect is modified as(71)
R
sum
e
t
=
∑
m
=
1
2
R
m
+
R
7
+
∑
n
=
11
12
R
n
+
1
2
/
R
x
0
+
1
/
R
x
1
+
1
/
R
x
2
+
1
1
/
∑
j
=
3
6
R
j
+
1
/
∑
k
=
8
10
R
k
+
s
0
2
-
z
2
2
μ
0
A
z
s
0
.The relative permeability of the stator and rotor core is larger compared to that of the air gap and the permanent magnet. Then (71) can be simplified.Due to rotor position varying inz-axis, the dynamic displacement stiffness of the RMB unit is modified as
(72)
k
r
s
z
d
=
d
f
x
s
d
(
δ
sin
(
ω
t
)
)
=
-
F
m
2
/
(
R
2
0
+
R
x
0
/
4
+
(
s
0
2
-
z
2
)
/
2
μ
0
A
z
s
0
)
2
8
μ
0
A
z
s
0
(
R
met
R
sum
e
t
)
2
=
k
r
s
d
λ
z
r
2
.Based on (72), the z term is included in the dynamic displacement stiffness model which will be affected by the rotor position varying in z-axis.


## 3.1. The Dynamic Current Stiffness of the AMB Unit Affected by the Rotor Position of RMB Unit inx-Axis

Based on the small variation of the rotor position of the RMB unit inx-axis, x
1
=
s
0
+
x, x
2
=
s
0
-
x, and the remains stationary in y-axis (y
1
=
y
2
=
s
0).The static bias flux in the air gap of the AMB unit is modified as(62)
ϕ
z
m
0
r
=
ϕ
z
1
r
0
=
ϕ
z
2
r
0
=
ϕ
m
0
r
2
=
F
m
R
2
0
+
1
/
(
2
/
R
x
0
+
1
/
R
x
1
+
1
/
R
x
2
)
+
R
z
1
0
/
2
=
F
m
R
2
0
+
1
/
(
2
μ
0
A
r
/
s
0
+
2
μ
0
A
r
s
0
/
(
s
0
2
-
x
2
)
)
+
R
z
1
0
/
2
.And the dynamic current stiffness of the AMB unit is modified due to rotor position varying inx-axis. Consider the following:
(63)
k
z
i
r
i
=
2
N
z
μ
0
A
z
R
z
i
0
(
1
+
λ
)
×
F
m
R
2
0
+
1
/
(
2
μ
0
A
r
/
s
0
+
2
μ
0
A
r
s
0
/
(
s
0
2
-
x
2
)
)
+
R
z
1
0
/
2
=
k
z
i
i
R
2
0
+
R
13
/
4
+
R
z
1
0
/
2
R
2
0
+
1
/
(
2
μ
0
A
r
/
s
0
+
2
μ
0
A
r
s
0
/
(
s
0
2
-
x
2
)
)
+
R
z
1
0
/
2
=
k
z
i
i
λ
r
z
,
where λ
r
z
=
(
R
2
0
+
R
13
/
4
+
R
z
1
0
/
2
)
/
(
R
2
0
+
1
/
(
2
μ
0
A
r
/
s
0
+
2
μ
0
A
r
s
0
/
(
s
0
2
-
x
2
)
)
+
R
z
1
0
/
2
) is the coefficient which reflects the dynamic current stiffness dependence of the AMB unit on the rotor position in x-axis (or y-axis). In other words, the coefficient, λ
r
z, reflects the coupling extent of the dynamic current stiffness of the AMB unit affected by the rotor position in x-axis or y-axis.From (63), the x term is included in the dynamic current stiffness model which will be affected by the rotor position of the RMB unit in x-axis or y-axis.


## 3.2. The Dynamic Displacement Stiffness of the AMB Unit Affected by the Rotor Position of RMB Unit inx-Axis

TheR
sum
0 in (36) will be modified when considering the small variation of the rotor position of the RMB unit in x-axis. Consider the following:
(64)
R
sum
0
=
R
2
0
+
1
2
μ
0
A
r
/
s
0
+
2
μ
0
A
r
s
0
/
(
s
0
2
-
x
2
)
+
R
z
1
0
R
z
2
0
R
z
1
0
+
R
z
2
0
.And the dynamic displacement stiffness is modified due to rotor position varying inx-axis. Consider the following:
(65)
k
z
s
r
d
=
-
1
2
μ
0
A
z
s
0
×
(
F
m
R
2
0
+
1
/
(
2
μ
0
A
r
/
s
0
+
2
μ
0
A
r
s
0
/
(
s
0
2
-
x
2
)
)
+
R
z
1
0
/
2
)
2
×
(
R
z
1
0
2
R
z
2
e
t
+
R
z
2
0
2
R
z
1
e
t
)
=
k
z
s
d
(
R
2
0
+
R
13
/
4
+
R
z
1
0
/
2
)
2
(
R
2
0
+
1
/
(
2
μ
0
A
r
/
s
0
+
2
μ
0
A
r
s
0
/
(
s
0
2
-
x
2
)
)
+
R
z
1
0
/
2
)
2
=
k
z
s
d
λ
r
z
2
.From (65), the x term is included in the dynamic displacement stiffness model which will be affected by the rotor position of the RMB unit in x-axis or y-axis. And λ
r
z
2 is the coefficient which reflects the dynamic displacement stiffness of the AMB unit dependence on the rotor position in x-axis (or y-axis). In other words, the coefficient, λ
r
z
2, reflects the coupling extent of the dynamic displacement stiffness of the AMB unit affected by the rotor position in x-axis or y-axis.


## 3.3. The Dynamic Current Stiffness of the RMB Unit Affected by the Rotor Position of AMB Unit inz-Axis

Based on the small variation of the rotor position of the AMB unit inz-axis, z
1
=
s
0
+
z, z
2
=
s
0
-
z, and the remains stationary in y-axis (R
y
1
=
R
y
2
≈
R
x
0). The magnetic reluctance of the air gap for the AMB unit can be modified as
(66)
R
z
1
=
s
0
+
z
μ
0
A
z
,
R
z
2
=
s
0
-
z
μ
0
A
z
.The static bias flux in the air gap of the RMB unit affected by the rotor position inz-axis varying is modified as
(67)
ϕ
r
m
=
ϕ
m
0
4
=
F
m
/
4
R
2
0
+
R
13
/
4
+
1
/
(
1
/
R
z
1
+
1
/
R
z
2
)
=
F
m
/
4
R
2
0
+
R
13
/
4
+
(
s
0
2
-
z
2
)
/
2
μ
0
A
z
s
0
.And the dynamic current stiffness of the RMB unit affected by the rotor position varying inz-axis is modified as
(68)
k
r
i
z
i
=
d
f
x
i
d
i
x
=
2
ϕ
r
m
N
r
μ
0
A
r
R
x
=
N
r
2
μ
0
A
r
R
x
F
m
R
2
0
+
R
x
/
4
+
(
s
0
2
-
z
2
)
/
2
μ
0
A
z
s
0
=
k
r
i
i
R
2
0
+
R
x
/
4
+
R
z
1
0
/
2
R
2
0
+
R
x
/
4
+
(
s
0
2
-
z
2
)
/
2
μ
0
A
z
s
0
=
k
r
i
i
λ
z
r
,
where λ
z
r
=
(
R
2
0
+
R
x
/
4
+
R
z
1
0
/
2
)
/
(
R
2
0
+
R
x
/
4
+
(
s
0
2
-
z
2
)
/
2
μ
0
A
z
s
0
).Base on (65), the z term is included in the dynamic current stiffness model which will be affected by the rotor position varying in z-axis.


## 3.4. The Dynamic Displacement Stiffness of the RMB Unit Affected by the Rotor Position of AMB Unit inz-Axis

The total static magnetic reluctance will be modified when considering the small variation of the rotor position of the AMB unit inz-axis. Consider the following:
(69)
R
met
=
R
2
0
+
R
x
0
4
+
1
1
/
R
z
1
0
+
1
/
R
z
2
0
=
R
2
0
+
R
x
0
4
+
s
0
2
-
z
2
2
μ
0
A
z
s
0
.The static bias flux in the air gap of the RMB unit is modified when considering the rotor position varying of the AMB unit inz-axis. Consider the following:
(70)
ϕ
x
0
=
F
m
4
R
met
=
F
m
4
(
R
2
0
+
R
x
0
/
4
+
(
s
0
2
-
z
2
)
/
2
μ
0
A
z
s
0
)
.The total dynamic magnetic reluctance affected by the eddy-current effect is modified as(71)
R
sum
e
t
=
∑
m
=
1
2
R
m
+
R
7
+
∑
n
=
11
12
R
n
+
1
2
/
R
x
0
+
1
/
R
x
1
+
1
/
R
x
2
+
1
1
/
∑
j
=
3
6
R
j
+
1
/
∑
k
=
8
10
R
k
+
s
0
2
-
z
2
2
μ
0
A
z
s
0
.The relative permeability of the stator and rotor core is larger compared to that of the air gap and the permanent magnet. Then (71) can be simplified.Due to rotor position varying inz-axis, the dynamic displacement stiffness of the RMB unit is modified as
(72)
k
r
s
z
d
=
d
f
x
s
d
(
δ
sin
(
ω
t
)
)
=
-
F
m
2
/
(
R
2
0
+
R
x
0
/
4
+
(
s
0
2
-
z
2
)
/
2
μ
0
A
z
s
0
)
2
8
μ
0
A
z
s
0
(
R
met
R
sum
e
t
)
2
=
k
r
s
d
λ
z
r
2
.Based on (72), the z term is included in the dynamic displacement stiffness model which will be affected by the rotor position varying in z-axis.


## 4. Example and Results

Based on the methods mentioned above, the CRAMB with permanent magnet is designed. The configuration of the CRAMB studied is shown in Figure1. The related parameters are shown in Table 1.The static current stiffness and displacement stiffness of the CRAMB could be calculated by linearized model with a small perturbation of the current or air gap. The static current stiffness,k
z
i
0, and displacement stiffness, k
z
s
0, of the AMB unit are 825.6 N/A and −8.0067 × 105 N/m, respectively. The static current stiffness, k
r
i
0, and displacement stiffness, k
r
s
0, of the RMB unit are 756.8 N/A and −1.5131 × 106 N/m, respectively.The magnitude and phase plots of the dynamic current stiffness of the AMB unit are shown in Figure6. It is obvious that the dynamic current stiffness is affected by the varying frequency of the control current of the AMB unit. And this can be shown in Figure 7. Figure 7 gives the magnitude of the ratio of the total dynamic reluctance to the total static reluctance in equivalent control flux path of the AMB unit, 20
lo
g
10
(
R
e
(
ω
)
/
R
z
i
0
). The dynamic reluctance of the control flux path of the AMB unit will increase with the varying frequency of the control current in the AMB coil current.Frequency response of the dynamic current stiffness of AMB unit. (a) The magnitude of the dynamic current stiffness of AMB unit. (b) The magnitude of the ratio of the dynamic current stiffness to the static current stiffness of the AMB unit. (c) The phase of the ratio,k
z
i
i
/
k
z
i
0, of the dynamic current stiffness to the static current stiffness of the AMB unit.
(a)
(b)
(c)Figure 7
The magnitude of the ratio of the total dynamic reluctance to the total static reluctance in equivalent control flux path of the AMB unit,20
lo
g
10
(
R
e
(
ω
)
/
R
z
i
0
).The magnitude and phase plots of the dynamic displacement stiffness of the AMB unit are shown in Figure8. It is obvious that the dynamic displacement stiffness is also affected by the varying frequency of the rotor position of the AMB unit in z-axis. From Figures 6 and 8, we can see that the response bandwidth of the dynamic displacement is higher compared to that of the dynamic current stiffness. And this can be explained by Figure 9, which shows the magnitude of the ratio of the total dynamic reluctance to the total static reluctance in equivalent bias flux path of the CRAMB, 20
lo
g
10
(
R
sum
e
t
/
R
sum
0
). The dynamic reluctance of the bias flux path will increase with the varying frequency.Frequency response of the dynamic displacement stiffness of the AMB unit. (a) The magnitude of the dynamic displacement stiffness of the AMB unit. (b) The magnitude of the ratio of the dynamic displacement stiffness to the static displacement stiffness of the AMB unit. (c) The phase of the ratio,k
z
s
d
/
k
z
s
0, of the dynamic displacement stiffness to the static displacement stiffness of the AMB unit.
(a)
(b)
(c)Figure 9
The magnitude of the ratio of the total dynamic reluctance to the total static reluctance in equivalent bias flux path of the CRAMB,20
lo
g
10
(
R
sum
e
t
/
R
sum
0
).The magnitude and phase plots of the dynamic displacement stiffness of the RMB unit are shown in Figure10. It is obvious that the dynamic displacement stiffness is also affected by the varying frequency of the rotor position of the RMB unit in x-axis. From Figures 8 and 10, we can see that the response bandwidth of the dynamic displacement of the AMB unit is higher compared to that of the dynamic displacement stiffness of the RMB unit.Frequency response of the dynamic displacement stiffness of the RMB unit. (a) The dynamic displacement stiffness of the RMB unit. (b) The ratio of the dynamic displacement stiffness to the static displacement stiffness of the RMB unit (c) Phase ofk
r
s
d
/
k
r
s
0.
(a)
(b)
(c)The magnitude of the dynamic current stiffness of the AMB unit versus varying frequency of the control current in AMB unit and the rotor position of RMB unit inx-axis (or y-axis) is shown in Figure 11. And the magnitude of the dynamic displacement stiffness of the AMB unit versus varying frequency of the rotor position of the AMB unit in z-axis and the rotor position of RMB unit in x-axis is shown in Figure 12. It is obvious that the dynamic current stiffness and the dynamic displacement stiffness of the AMB unit are determined by the varying frequency of the control current in the AMB unit or the rotor position in z-axis rather than the rotor position of the RMB unit in x-axis (or y-axis). The reason is the coefficients (shown in Figure 13) which are little affected by the rotor position in x-axis (or y-axis) direction.Figure 11
The magnitude of the dynamic current stiffness versus varying frequency of the control current in AMB unit and the rotor position of RMB unit inx-axis.Figure 12
The magnitude of the dynamic displacement stiffness versus varying frequency of the rotor position of the AMB unit inz-axis and the rotor position of RMB unit in x-axis.Figure 13
The coefficientsλ
r
z and λ
r
z
2.The magnitude of the dynamic current stiffness of the RMB unit versus varying frequency of the control current in RMB unit and the rotor position inz-axis is shown in Figure 14. And the magnitude of the dynamic displacement stiffness of the RMB unit versus varying frequency of the rotor position in x-axis (or y-axis) and the rotor position in z-axis is shown in Figure 15. It is obvious that the dynamic current stiffness and the dynamic displacement stiffness of the RMB unit are determined by the varying frequency of the control current in the RMB unit or the rotor position in x-axis (or y-axis) rather than the rotor position of the AMB unit in z-axis. The reason is the coefficients (shown in Figure 16) which are little variation affected by the rotor position in z-axis.Figure 14
The magnitude of the dynamic current stiffness of the RMB unit versus varying frequency of the control current of the RMB unit and the rotor position inz-axis.Figure 15
The magnitude of the dynamic displacement stiffness of the RMB versus varying frequency of the rotor position inx-axis and the rotor position of the AMB unit in z-axis.Figure 16
The coefficientsλ
z
r and λ
z
r
2.


## 5. Conclusion

The CRAMB with permanent magnet creating bias flux can reduce the size, cost, mass, and save energy of the magnetic bearing system [6, 7, 12, 16, 17]. The structure and its simplified dynamic stiffness models of a novel CRAMB with permanent magnet bias are presented when considering the eddy-current effect in this paper. And the dynamic cross-coupling models between axial and radial bearing units are also given using the method of equivalent magnetic circuit including eddy-current effect. An example is given and the analytical results show that the eddy-current effects have a fundamental impact on the dynamic stiffness of the CRAMB, but the dynamic current stiffness of the RMB unit is approximately equal to its static current stiffness. The bandwidths of the dynamic current stiffness of the AMB unit and the dynamic displacement stiffness of the CRAMB are affected by the time-varying control currents or air gap respectively. The rotor position varying in x-axis (or y-axis) has little influence on the dynamic current stiffness and the dynamic displacement stiffness of the AMB unit due to the few coupling coefficients. And the rotor position varying in z-axis also has little influence on the dynamic current stiffness and the dynamic displacement stiffness of the RMB unit due to the little coupling coefficients. The described method can simulate the dynamic behavior of the CRAMB with permanent magnet.


## Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of the paper.



---

*Source: 340140-2014-01-12.xml*