# Unbalanced Magnetic Pull Effect on Stiffness Models of Active Magnetic Bearing due to Rotor Eccentricity in Brushless DC Motor Using Finite Element Method

Authors: Bangcheng Han; Shiqiang Zheng; Xiquan Liu
Journal: Mathematical Problems in Engineering
 (2013)

---

## Abstract

We firstly report on an investigation into the unbalanced magnetic pull (UMP) effect on the static stiffness models of radial active magnetic bearing (RAMB) in brushless DC motor (BDCM)  in no-loaded and loaded conditions using the finite element method (FEM). The influences of the UMP on the force-control current, force-position, current stiffness, and position stiffness of RAMB are clarified in BDCM with 100 kW rated power. We found the position stiffness to be more susceptible to UMP. The primary source of UMP is the permanent magnets of BDCM. In addition, the performance of RAMB is affected by the UMP ripples during motor commutation and also periodically affected by the angular position of rotor. The characteristic curves of RAMB force versus control current (or rotor position) and angular position of rotor affected by the UMP are given. The method is useful in design and optimization of RAMB in magnetically suspended BDCMs.


---

## Body

## 1. Introduction

The number of commercial/industrial applications for magnetic bearings is large and still growing because of the unique advantages: oil-free operation, vibration control, and on-line monitoring. Nowadays, magnetic bearings have been widely used in refrigeration compressors [1], turbo expanders, blower industry (specifically in the wastewater treatment area), pumps, and high-speed machines [2–5]. These systems are composed with magnetic bearings (MBs) and motor, and the different components are assembled on the rotor shaft which is  supported by the MBs. The rotor supported by the MBs would be naturally combined by the mechanical model.  Due to manufacturing imperfectness of mechanical parts and/or the use of magnetic bearings in high-speed motors, the rotor eccentricity is inevitable. The magnetic field in motors will create tangential electromagnetic forces on the rotor which will produce torque. The radial forces named as unbalanced magnetic pull (UMP) [6] are generated at the same time in which case the UMP is developed and in interaction with the rotor supported by magnetic bearings and they may cause significant effects on the static mathematical model of magnetic bearings. An analytical model for analyzing permanent magnet (PM) motors is presented for predicting of magnetic field and magnetic forces including UMP [7]. A simple analytical model which is able to account for the damping effects of a cage rotor on the UMP when the rotor is dynamically eccentric is presented [8]. The UMPs which act upon the rotor of a salient-pole synchronous generator due to the eccentric motion of the rotor shaft in the presence of magnetic field in no-loaded and loaded conditions are calculated using the finite element method [9]. Dynamic eccentricity up to 80% is put into the model of a large brushless rare-earth permanent magnet motors model and a variety of simulations carried out to investigate the UMP [10]. A general analytical model is developed to predict the UMP, which results in permanent-magnet brushless ac and dc machines having a diametrically asymmetric disposition of slots and phase windings [11]. An investigation on the UMP in ferrite-magnet fractional-slot brushless permanent-magnet motors due to either magnetic asymmetry or static rotor eccentricity is reported [12]. An analytic model for analyzing the UMP induced by electromagnetic structure and drive current of the BDCM is presented, and the model can simplify effectively the analysis procedure [13]. The analytical justifications for the use of a numerical impulse method to calculate the effects of the slot harmonics on the UMP are presented [14].  Maxwell stress tensor (MST) method is used to compute the UMP based on the 2D finite element (FE) calculation [15]. The UMP in a three-phase generator, caused by relative eccentricity, is calculated theoretically, and its effect on vibration of an eccentric rotor is studied [16]. The vibration sources of rotor in permanent magnet (PM) machine include mainly the centrifugal force generated by eccentricity and the UMP are reported. The finite element analysis combined with the Newmark method was used to calculate the nonlinear unbalance response due to eccentricity and the unbalanced magnetic pull [17]. Influences of rotor mass eccentric on suspension performance of bearingless switched reluctance motors were analyzed [18]. An on-line identification method of suspension force and magnetic unbalance pull force parameters in buried-type IPM bearingless motors has been presented [19]. A sort of beat appears in monitored radial displacements, caused by the interaction between the unbalance response of the rotor, and the UMP in induction motor is reported [20].The emphases of the literatures mentioned previously are the analytical model of UMP and its effect on vibration, but few literatures focus on the UMP’s effects on the stiffness models of AMB. The magnetic forces have been computed from the results of static magnetic field analysis. Thus, the potential circulating currents in the windings and eddy currents in conducting parts caused by eccentric rotor motion do not affect the magnetic forces, but they will make the core hotter. The main scope of this paper is to calculate and analyze the static UMP which can have an effect on the performance of radial active magnetic bearing (RAMB) due to static rotor eccentricity in BLDCM which is supported by AMBs. The UMP in 100 kW BLDCM is calculated owning to the static rotor eccentricity in no-loaded and loaded conditions by using FEM, and the performance of the RAMB is also analyzed by the characteristic curves of force-control current and force-position. Then, the influences of the UMP on the force-control current, the force-position, the current stiffness, and the position stiffness of RAMB are analyzed. The method is useful in design and optimization of RAMB in the magnetically suspended PM motor’s design stage.


## 2. Analysis of the UMP Effect on Characteristics of the RAMB

The configuration of a BDCM which is supported by two RAMBs and a pair of axial AMBs is shown in Figure1, and the configurations of BDCM and RAMB are shown in Figure 2. In BDCMs, the motor torque T and the UMP  fUMP are produced due to the rotor eccentricity, and the UMP fUMP will affect the characteristics of RAMBs.Figure 1
The configuration of BDCM which is supported by AMBs. (1) Rotor, (2) radial position sensors, (3) RAMB, (4) BDCM, (5) axial AMB, and (6) axial position sensor.The configurations of the BDCM and the RAMB. (a) The configuration of BDCM. (b) The configuration of RAMB.
(a)
(b)From Figure2(b), the non-linear net magnetic force [21] that is exerted on the rotor in differential driving modes can be given for x-axis as
(1)fx=14μ0n2A[(i0+ix)2(s0-x)2-(i0-ix)2(s0+x)2]cosβ,
with  k=(1/4)μ0n2Acosβ, where μ0 is permeability of free space (μ0=4π×10-7Vs/Am), n is the number of turns in the electromagnet’s coil, A is the cross-section of the iron, β is the angle under which the magnetic forces affect the rotor for each of the poles (in the case of RAMB with 8 poles, β=22.5°), the bias current i0 is passed through both coils, the control current ix is added to the coil exerting forces in the positive direction and subtracted from the opposite coil, s0 is the nominal air gap, and the displacement x is the distance of rotor depart from its center position in x-axis direction.Consideringx≪s0 and ix≪i0, (1) can be linearized by using Taylor’s series. It yields the following typical AMB relation:
(2)fx=4ki0s02ix+4ki02s03x=kiix+ksx,
where ki=4ki0/s02 is the current stiffness and ks=4ki02/s03 is the displacement stiffness. The two kinds stiffness are constants depending on the magnetic bearing geometry (air gap, dimensions, material, and the coil turns) and the bias current without considering the UMP.The UMP will have an effect on the net magnetic force of RAMB in BDCM, and the net magnetic force discussed in (1) is rewritten as
(3)fx′=14μ0n2A[(i0+ix)2(s0-x)2-(i0-ix)2(s0+x)2]cosβ+fUMPcosα2,
where α is the angle between the vector of the UMP fUMP  and the x-axis. From (3), the net magnetic force, the current stiffness, and the displacement stiffness are affected by the UMP. There are two RAMBs in the system shown in Figure 1, and thereby the UMP is divided by 2.


## 3. Finite Element Model for Calculation of the UMP Effect on RAMB

The eccentric motion of the rotor supported by AMBs in no-loaded or loaded operation will result in the UMP which acts upon the rotor and consequently on the magnetic bearings of the BDCM. In this paper, the goal is to use the FEM to simulate the UMP and calculate the results of the UMP effect on the characteristics of RAMB in BDCM.An example of magnetically suspended BDCM is given, and the design parameters of the BDCM and the RAMB are shown in Tables1 and 2, respectively. The example is used to clarify the UMP effect on the characteristics of RAMB in magnetically suspended BDCM for no-loaded and loaded conditions, and so forth.Table 1
Designed parameters of the BDCM.
Parameters
Value
Weight of the rotor, (kg)
23.6
Number of pole pairs
2
Rated power, (kW)
100
Rated voltage, (V)
480
Maximum speed, (r/min)
32000
Number of stator slots
24
Connection
Star
Winding arrangement (series)
Double layer
Parallel branches of the stator winding
4
Rotor outer diameter, (mm)
89
Air gap length, (mm)
4
Stator inner diameter, (mm)
97
Axial length, (mm)
150Table 2
Designed parameters of the RAMB.
Parameters
Value
Stator outer diameter, (mm)
228
Stator inner diameter, (mm)
89.8
Air gap length, (mm)
0.4
Rotor outer diameter, (mm)
89
Rotor inner diameter, (mm)
60
Number of turns
160
Bias current, (A)
1
Axial length, (mm)
16


## 4. Simulation Results and Discussion

### 4.1. The Load Currents in BDCM and Finite Element Model

The FEM simulations used to calculate the UMP’s effects on characteristics of RAMB resulting from the eccentric of rotor for no-loaded and loaded conditions have been carried out in this section. As shown in Table3, the UMP calculations at different loads are conducted in the following three cases.(a)
The UMP calculation without stator current which is used to clarify the UMP produced by the permanent magnet in BDCM effect on the characteristics.(b)
The UMP calculation for no-loaded condition which is used to clarify the UMP produced by the permanent magnet and no-loaded condition in BDCM effect on the characteristics.(c)
The UMP calculation for loaded condition which is used to clarify the UMP produced by the permanent magnet and rated loaded condition in BDCM effect on the characteristics.Table 3
UMP calculation at different loads.
Item
Stator current (A)
Without stator current in BDCM
0
No load in BDCM
3.2
Rated load in BDCM
208The 2-dimensional (2D) finite element (FE) models of BDCM and RAMB are established and are shown in Figure2 except the air. The B-H curve of ferromagnetic material which is used for motor stator, RAMB’s stator, and RAMB’s rotor is shown in Figure 3 when the nonlinear analysis is considered. The software of Ansoft Maxwell is used to calculate electromagnetic field. The 2D FE mesh model of BDCM is shown in Figure 4(a), the air is not shown to give a better view of model, and there are 42748 elements totally. The FE mesh model of RAMB is shown in Figure 4(b), the air is also not shown, and there are 25942 elements in total. The vector potential boundary is available for use in Maxwell 2D.Figure 3
B-H curve for the ferromagnetic material.The 2-dimentional FE mesh model of BDCM and RAMB except the air. (a) The 2-dimentional FE mesh model of BDCM except the air. (b) The FE mesh model of RAMB except the air.
(a)
(b)The lines of magnetic flux for BDCM in loaded condition are shown in Figure5(a). The lines of magnetic flux for RAMB with bias current are shown in Figure 5(b).The lines of magnetic flux for the BDCM and the RAMB. (a) The lines of magnetic flux for BDCM in loaded condition. (b) The lines of magnetic flux for RAMB with bias current.
(a)
(b)


### 4.2. The Static Linearized Model of the RAMB and the UMP Produced by Permanent Magnet in BDCM

In the rotor center, the linearized model of RAMB without considering UMP through the 2D FEM is given in (4). The current stiffness ki=297.2 N/A, and the position stiffness ks=-772020 N/m:
(4)fx=kiix+ksx=297.2ix-772020x.The UMP is produced by the permanent magnet in BDCM without stator current. The following analysis considers the eccentricity and angular position of rotor. The eccentricity changes from −200μm to 200 μm and the angular position of motor rotor changes from 0° to 180°. The characteristic curves of UMP-eccentricity-angular position of rotor are shown in Figure 6. The UMP is only calculated in x-axis as shown in Figure 2(a). It is obvious that the characteristic curves are linear; the UMP is determined by the eccentricity and is independent of the angular of motor rotor. The UMP is 38.7 N with the eccentricity of −200 μm. The relationship between the UMP and the eccentricity can be given as follows:
(5)fUMP=kepex=-193220ex,
where kep is the slope of  fUMP versus the eccentricity and the ex is the eccentricity. This is a disadvantage because the additional magnetic bearing force (MBF) is necessary to overcome the negative slope coefficient, limiting the achievable control performance.Figure 6
The UMP produced by permanent magnet in BDCM versus eccentricity and angular position of motor rotor.


### 4.3. The Characteristic Curves of the MBF Affected by the UMP of BDCM with No-Loaded Condition

The relationship between the MBF of RAMB and the control current is calculated by FEM and is shown in Figure7. The characteristic curves of MBF versus the control current which is affected by the UMP produced by the permanent magnet in BDCM for the different values of eccentricity are also shown in Figure 7. The current stiffness curves of RAMB are shown in Figure 8, and these 5 curves are coincident. It is obvious that the UMP only changes the amplitude of MBF without changing its current stiffness and the linear relationship between the MBF and the control current. The additional MBF produced by control current is necessary to overcome the UMP. For example, the 38.7 N provided by RAMB is necessary to overcome the UMP with eccentricity of −200 μm. The mathematic model of the MBF versus control current affected by the UMP can be given as follows:
(6)fi=kiix+kepex≈297.2ix-193220ex.Figure 7
The characteristic curves of the MBF versus control current affected by the UMP which is produced by the PM in BDCM.Figure 8
The characteristic curves of the current stiffness affected by the UMP which is produced by the PM in BDCM.It is different from the characteristic curves of MBF versus control current affected by the UMP; the UMP affects the amplitude and linear relationship between the MBF versus rotor position of the RAMB, as the characteristic curves are shown in Figure9. The UMP changes the equilibrium position of the RAMB; the additional control current of RAMB is necessary to maintain system stability. The position stiffness of RAMB is also varied by the UMP; and is shown in Figure 10. The absolute value of position stiffness in rotor center position increases from 772020 N/m (No UMP) to 960346 N/m (affected by UMP), namely, increases by 24.4%. The position stiffness when UMP is not generated is different from those when the rotor is eccentric. The reason is there is no UMP when the rotor located its center, and only the torque is generated by the BDCM, but the UMP will be generated by the PM and (or) the control current in BDCM due to the rotor eccentric, then the characteristic curves of the MBF will be modified by the UMP. Thus, the position stiffness model is modified. The rate of the change of the UMP produced by BDCM versus the rotor eccentricity is the same, so it is equal between the eccentricities of 100 μm and 200 μm. The position stiffness affected by the different values of eccentricity is the same, and their curves are coincident. The mathematic model of the MBF versus rotor position affected by the UMP can be given as follows:
(7)fs=kspx+kepex≈-960346x-193220ex.Figure 9
The characteristic curves of the MBF versus rotor position affected by the UMP of BDCM with no-loaded condition.Figure 10
The characteristic curves of the position stiffness affected by the UMP of BDCM with no-loaded condition.Equations (6) and (7) are substituted into (3) considering the effect of the UMP produced by permanent magnet in BDCM in x-axis, which gives
(8)fx′=kiix+kspx+2kepexcos0°=297.2ix-960346x-386440ex.The UMP is produced by the permanent magnet and stator current (3.2 A) in BDCM for no-loaded operation. The following analysis considers the eccentricity and angular position of rotor. The eccentricity changes from −200μm to 200 μm, and the angular position of motor rotor changes from 0° to 180°. The characteristic curves of UMP-eccentricity-angular position of rotor are shown in Figure 11. The UMP is also calculated in x-axis as shown in Figure 2(a). It is obvious that the characteristic curves in no-loaded condition are the same as the UMP produced by the permanent magnet in BDCM.  The UMP produced by the stator current in no-loaded condition is very small and can be neglected for the large air gap.Figure 11
The UMP in no-loaded condition versus rotor eccentricity and angular position of motor rotor.


### 4.4. The Characteristic Curves of the MBF Affected by the UMP of BDCM in Loaded Condition

The previous analysis shows that the UMP for no-loaded condition or without stator current in BDCM is determined by the value of eccentricity and is independent of the angular of motor rotor, but the UMP for loaded (stator current: 208 A) condition is not only affected by the eccentricity but also affected by the angular position of motor rotor with the larger value of eccentricity. The characteristic curves of the UMP in loaded condition versus eccentricity and angular position of motor rotor are shown in Figure12. For example, the value of UMP will change from 37.1 N to 40.8 N with the eccentricity with −200 μm when the motor rotates from 0–180 degree, and the maximum error of UMP is 3.7 N.The UMP in loaded condition (stator current: 208 A) versus eccentricity and angular position of motor rotor. (a) 2D. (b) 3D.
(a)
(b)The UMP is produced by the permanent magnet and stator current in BDCM for loaded condition (stator current: 208 A). The characteristic curves of MBF versus the control current affected by the UMP in loaded condition for the different values of eccentricity are shown in Figure13, and the current stiffness curves of RAMB shown in Figure 14 are coincident. It is obvious that the UMP for loaded operation only changes the amplitude of MBF rather than changing its current stiffness and the linear relationship between the MBF and the control current. The additional MBF is necessary to overcome the UMP, such as 40.8 N provided by the RAMB is necessary to overcome the UMP with eccentricity of −200 μm. The UMP increases from 38.7 N to 40.8 N compared to the UMP produced by permanent magnet in BDCM, and the difference is produced by the stator current in BDCM. The mathematic model of the MBF versus control current affected by the UMP in loaded condition can be given as follows:
(9)fi=kiix+kepex≈297.2ix-203041ex,
where kep increases by 5% compared with (6); the calculation results show clearly that the UMP is mainly produced by the permanent magnet in BDCM.Figure 13
The characteristic curves of the MBF versus control current affected by the UMP for loaded operation in BDCM.Figure 14
The characteristic curves of the current stiffness affected by the UMP for loaded operation in BDCM.It is also different from the characteristic curves of MBF versus control current affected by the UMP in BDCM for loaded operation; both the amplitude and the linear relationship between the MBF versus rotor position of the RAMB are affected by the UMP; the characteristic curves are shown in Figure15. The UMP changes the equilibrium position of RAMB; the additional control current of RAMB is necessary to maintain system stability. The position stiffness of RAMB is also varied by the UMP; and is shown in Figure 16. The absolute value of position stiffness in rotor center position increases from 772020 N/m (No UMP) to 968599 N/m (affected by UMP), namely, increasing by 25.5%. The position stiffness affected by the different value of eccentricity is also the same, and these curves are coincident. The mathematic model of the MBF versus rotor position affected by the UMP can be given as follows:
(10)fs=kspx+kepex≈-968599x-203041ex.Figure 15
The characteristic curves of the MBF versus rotor position affected by the UMP for loaded operation in BDCM.Figure 16
The characteristic curves of the position stiffness affected by the UMP for loaded operation in BDCM.Equations (9) and (10) are substituted into (3) considering the effect of the UMP in BDCM for loaded operation in x-axis, which gives
(11)fx′=kiix+kspx+2kepexcos0°=297.2ix-968599x-406082ex.The commutation behavior of BDCM also affects the UMP. The control currents in BDCM will change according to the rotor position, and the PM of the BDMC is not axisymmetrical. Thus, the value of UMP will change at different commutation position due to the rotor eccentricity. So, the UMP during commutation will affect the MBF. The characteristic curves of the MBF versus control current and the MBF versus rotor position are shown in Figures17 and 18 with different eccentricities. It can be calculated that the difference during commutation is 3.44 N for the MBF versus control current and 3.42 N for the MBF versus rotor position affected by the UMP in BDCM for loaded operation with eccentricity −200 μm.Figure 17
The characteristic curves of the MBF versus control current affected by the UMP during commutation in BDCM for loaded operation with eccentricity±200 μm.Figure 18
The characteristic curves of the MBF versus rotor position affected by the UMP during commutation in BDCM for loaded operation with eccentricity±200 μm.The characteristic curves of the UMP in loaded condition versus angular position of rotor for different eccentricities are shown in Figure19. The value of UMP will change from 37.1 N to 40.8 N with the eccentricity of −200 μm when the motor rotates from 0–180 degree, the maximum difference of UMP is 3.7 N, and the period of the UMP ripples is 30°. The reason is that the BDCM with 24 stator slots and 4 magnetic poles is calculated in this paper.Figure 19
The characteristic curves of the UMP in loaded condition versus angular position of rotor for different eccentricity.


## 4.1. The Load Currents in BDCM and Finite Element Model

The FEM simulations used to calculate the UMP’s effects on characteristics of RAMB resulting from the eccentric of rotor for no-loaded and loaded conditions have been carried out in this section. As shown in Table3, the UMP calculations at different loads are conducted in the following three cases.(a)
The UMP calculation without stator current which is used to clarify the UMP produced by the permanent magnet in BDCM effect on the characteristics.(b)
The UMP calculation for no-loaded condition which is used to clarify the UMP produced by the permanent magnet and no-loaded condition in BDCM effect on the characteristics.(c)
The UMP calculation for loaded condition which is used to clarify the UMP produced by the permanent magnet and rated loaded condition in BDCM effect on the characteristics.Table 3
UMP calculation at different loads.
Item
Stator current (A)
Without stator current in BDCM
0
No load in BDCM
3.2
Rated load in BDCM
208The 2-dimensional (2D) finite element (FE) models of BDCM and RAMB are established and are shown in Figure2 except the air. The B-H curve of ferromagnetic material which is used for motor stator, RAMB’s stator, and RAMB’s rotor is shown in Figure 3 when the nonlinear analysis is considered. The software of Ansoft Maxwell is used to calculate electromagnetic field. The 2D FE mesh model of BDCM is shown in Figure 4(a), the air is not shown to give a better view of model, and there are 42748 elements totally. The FE mesh model of RAMB is shown in Figure 4(b), the air is also not shown, and there are 25942 elements in total. The vector potential boundary is available for use in Maxwell 2D.Figure 3
B-H curve for the ferromagnetic material.The 2-dimentional FE mesh model of BDCM and RAMB except the air. (a) The 2-dimentional FE mesh model of BDCM except the air. (b) The FE mesh model of RAMB except the air.
(a)
(b)The lines of magnetic flux for BDCM in loaded condition are shown in Figure5(a). The lines of magnetic flux for RAMB with bias current are shown in Figure 5(b).The lines of magnetic flux for the BDCM and the RAMB. (a) The lines of magnetic flux for BDCM in loaded condition. (b) The lines of magnetic flux for RAMB with bias current.
(a)
(b)


## 4.2. The Static Linearized Model of the RAMB and the UMP Produced by Permanent Magnet in BDCM

In the rotor center, the linearized model of RAMB without considering UMP through the 2D FEM is given in (4). The current stiffness ki=297.2 N/A, and the position stiffness ks=-772020 N/m:
(4)fx=kiix+ksx=297.2ix-772020x.The UMP is produced by the permanent magnet in BDCM without stator current. The following analysis considers the eccentricity and angular position of rotor. The eccentricity changes from −200μm to 200 μm and the angular position of motor rotor changes from 0° to 180°. The characteristic curves of UMP-eccentricity-angular position of rotor are shown in Figure 6. The UMP is only calculated in x-axis as shown in Figure 2(a). It is obvious that the characteristic curves are linear; the UMP is determined by the eccentricity and is independent of the angular of motor rotor. The UMP is 38.7 N with the eccentricity of −200 μm. The relationship between the UMP and the eccentricity can be given as follows:
(5)fUMP=kepex=-193220ex,
where kep is the slope of  fUMP versus the eccentricity and the ex is the eccentricity. This is a disadvantage because the additional magnetic bearing force (MBF) is necessary to overcome the negative slope coefficient, limiting the achievable control performance.Figure 6
The UMP produced by permanent magnet in BDCM versus eccentricity and angular position of motor rotor.


## 4.3. The Characteristic Curves of the MBF Affected by the UMP of BDCM with No-Loaded Condition

The relationship between the MBF of RAMB and the control current is calculated by FEM and is shown in Figure7. The characteristic curves of MBF versus the control current which is affected by the UMP produced by the permanent magnet in BDCM for the different values of eccentricity are also shown in Figure 7. The current stiffness curves of RAMB are shown in Figure 8, and these 5 curves are coincident. It is obvious that the UMP only changes the amplitude of MBF without changing its current stiffness and the linear relationship between the MBF and the control current. The additional MBF produced by control current is necessary to overcome the UMP. For example, the 38.7 N provided by RAMB is necessary to overcome the UMP with eccentricity of −200 μm. The mathematic model of the MBF versus control current affected by the UMP can be given as follows:
(6)fi=kiix+kepex≈297.2ix-193220ex.Figure 7
The characteristic curves of the MBF versus control current affected by the UMP which is produced by the PM in BDCM.Figure 8
The characteristic curves of the current stiffness affected by the UMP which is produced by the PM in BDCM.It is different from the characteristic curves of MBF versus control current affected by the UMP; the UMP affects the amplitude and linear relationship between the MBF versus rotor position of the RAMB, as the characteristic curves are shown in Figure9. The UMP changes the equilibrium position of the RAMB; the additional control current of RAMB is necessary to maintain system stability. The position stiffness of RAMB is also varied by the UMP; and is shown in Figure 10. The absolute value of position stiffness in rotor center position increases from 772020 N/m (No UMP) to 960346 N/m (affected by UMP), namely, increases by 24.4%. The position stiffness when UMP is not generated is different from those when the rotor is eccentric. The reason is there is no UMP when the rotor located its center, and only the torque is generated by the BDCM, but the UMP will be generated by the PM and (or) the control current in BDCM due to the rotor eccentric, then the characteristic curves of the MBF will be modified by the UMP. Thus, the position stiffness model is modified. The rate of the change of the UMP produced by BDCM versus the rotor eccentricity is the same, so it is equal between the eccentricities of 100 μm and 200 μm. The position stiffness affected by the different values of eccentricity is the same, and their curves are coincident. The mathematic model of the MBF versus rotor position affected by the UMP can be given as follows:
(7)fs=kspx+kepex≈-960346x-193220ex.Figure 9
The characteristic curves of the MBF versus rotor position affected by the UMP of BDCM with no-loaded condition.Figure 10
The characteristic curves of the position stiffness affected by the UMP of BDCM with no-loaded condition.Equations (6) and (7) are substituted into (3) considering the effect of the UMP produced by permanent magnet in BDCM in x-axis, which gives
(8)fx′=kiix+kspx+2kepexcos0°=297.2ix-960346x-386440ex.The UMP is produced by the permanent magnet and stator current (3.2 A) in BDCM for no-loaded operation. The following analysis considers the eccentricity and angular position of rotor. The eccentricity changes from −200μm to 200 μm, and the angular position of motor rotor changes from 0° to 180°. The characteristic curves of UMP-eccentricity-angular position of rotor are shown in Figure 11. The UMP is also calculated in x-axis as shown in Figure 2(a). It is obvious that the characteristic curves in no-loaded condition are the same as the UMP produced by the permanent magnet in BDCM.  The UMP produced by the stator current in no-loaded condition is very small and can be neglected for the large air gap.Figure 11
The UMP in no-loaded condition versus rotor eccentricity and angular position of motor rotor.


## 4.4. The Characteristic Curves of the MBF Affected by the UMP of BDCM in Loaded Condition

The previous analysis shows that the UMP for no-loaded condition or without stator current in BDCM is determined by the value of eccentricity and is independent of the angular of motor rotor, but the UMP for loaded (stator current: 208 A) condition is not only affected by the eccentricity but also affected by the angular position of motor rotor with the larger value of eccentricity. The characteristic curves of the UMP in loaded condition versus eccentricity and angular position of motor rotor are shown in Figure12. For example, the value of UMP will change from 37.1 N to 40.8 N with the eccentricity with −200 μm when the motor rotates from 0–180 degree, and the maximum error of UMP is 3.7 N.The UMP in loaded condition (stator current: 208 A) versus eccentricity and angular position of motor rotor. (a) 2D. (b) 3D.
(a)
(b)The UMP is produced by the permanent magnet and stator current in BDCM for loaded condition (stator current: 208 A). The characteristic curves of MBF versus the control current affected by the UMP in loaded condition for the different values of eccentricity are shown in Figure13, and the current stiffness curves of RAMB shown in Figure 14 are coincident. It is obvious that the UMP for loaded operation only changes the amplitude of MBF rather than changing its current stiffness and the linear relationship between the MBF and the control current. The additional MBF is necessary to overcome the UMP, such as 40.8 N provided by the RAMB is necessary to overcome the UMP with eccentricity of −200 μm. The UMP increases from 38.7 N to 40.8 N compared to the UMP produced by permanent magnet in BDCM, and the difference is produced by the stator current in BDCM. The mathematic model of the MBF versus control current affected by the UMP in loaded condition can be given as follows:
(9)fi=kiix+kepex≈297.2ix-203041ex,
where kep increases by 5% compared with (6); the calculation results show clearly that the UMP is mainly produced by the permanent magnet in BDCM.Figure 13
The characteristic curves of the MBF versus control current affected by the UMP for loaded operation in BDCM.Figure 14
The characteristic curves of the current stiffness affected by the UMP for loaded operation in BDCM.It is also different from the characteristic curves of MBF versus control current affected by the UMP in BDCM for loaded operation; both the amplitude and the linear relationship between the MBF versus rotor position of the RAMB are affected by the UMP; the characteristic curves are shown in Figure15. The UMP changes the equilibrium position of RAMB; the additional control current of RAMB is necessary to maintain system stability. The position stiffness of RAMB is also varied by the UMP; and is shown in Figure 16. The absolute value of position stiffness in rotor center position increases from 772020 N/m (No UMP) to 968599 N/m (affected by UMP), namely, increasing by 25.5%. The position stiffness affected by the different value of eccentricity is also the same, and these curves are coincident. The mathematic model of the MBF versus rotor position affected by the UMP can be given as follows:
(10)fs=kspx+kepex≈-968599x-203041ex.Figure 15
The characteristic curves of the MBF versus rotor position affected by the UMP for loaded operation in BDCM.Figure 16
The characteristic curves of the position stiffness affected by the UMP for loaded operation in BDCM.Equations (9) and (10) are substituted into (3) considering the effect of the UMP in BDCM for loaded operation in x-axis, which gives
(11)fx′=kiix+kspx+2kepexcos0°=297.2ix-968599x-406082ex.The commutation behavior of BDCM also affects the UMP. The control currents in BDCM will change according to the rotor position, and the PM of the BDMC is not axisymmetrical. Thus, the value of UMP will change at different commutation position due to the rotor eccentricity. So, the UMP during commutation will affect the MBF. The characteristic curves of the MBF versus control current and the MBF versus rotor position are shown in Figures17 and 18 with different eccentricities. It can be calculated that the difference during commutation is 3.44 N for the MBF versus control current and 3.42 N for the MBF versus rotor position affected by the UMP in BDCM for loaded operation with eccentricity −200 μm.Figure 17
The characteristic curves of the MBF versus control current affected by the UMP during commutation in BDCM for loaded operation with eccentricity±200 μm.Figure 18
The characteristic curves of the MBF versus rotor position affected by the UMP during commutation in BDCM for loaded operation with eccentricity±200 μm.The characteristic curves of the UMP in loaded condition versus angular position of rotor for different eccentricities are shown in Figure19. The value of UMP will change from 37.1 N to 40.8 N with the eccentricity of −200 μm when the motor rotates from 0–180 degree, the maximum difference of UMP is 3.7 N, and the period of the UMP ripples is 30°. The reason is that the BDCM with 24 stator slots and 4 magnetic poles is calculated in this paper.Figure 19
The characteristic curves of the UMP in loaded condition versus angular position of rotor for different eccentricity.


## 5. Conclusions

In this paper, the investigation of the UMP effect on the stiffness models of AMB due to the rotor static eccentricity in BDCM which is supported by AMBs is reported. The nonlinear FEM is used to calculate the UMP effect on the characteristics of RAMB in BDCM with 100 kW rated power in no-loaded and loaded conditions, such as the force-current, the current stiffness, the force-displacement, and the displacement stiffness.The analytic results verified that the force-current, the force-displacement, and the displacement stiffness of RAMB are more susceptible to UMP, but the UMP has no effect on the current stiffness of RAMB. The primary source of UMP is the permanent magnets of BDCM, while the control current in BDCM has little influence on the UMP. For example, when the BDCM operates at full load, more than 95% of the UMP is produced by the permanent magnet of BDCM.It will generate UMP ripples during motor commutation; the UMP after commutation is different from that before. The difference will affect the characteristics of the RAMB. The performance of RAMB is also periodically affected by the angular position of rotor and eccentricity. The amplitude of vibration is determined by the value of eccentricity, and the dynamic models would be considered in the future.The method is useful in design and optimization of AMB in magnetically suspended BDCMs and is of concern in design stage.



---

*Source: 745912-2013-07-14.xml*