Open access peer-reviewed chapter

Electromagnetic Levitation System for Active Magnetic Bearing Wheels

Written By

Yonmook Park

Submitted: 28 May 2016 Reviewed: 13 December 2016 Published: 31 May 2017

DOI: 10.5772/67227

From the Edited Volume

Bearing Technology

Edited by Pranav H. Darji

Chapter metrics overview

1,763 Chapter Downloads

View Full Metrics

Abstract

In this chapter, the author presents an electromagnetic levitation system for active magnetic bearing wheels. This system consists of a rotor, a shaft, a cover, and a base. The author derives a meaningful electromagnetic force by using the singular value decomposition. The author develops a control system using the proportional‐integral‐derivative controller to control the position of the rotor and regulate the two gimbal angles of the rotor. The author gives the numerical simulation and experimental results on the control of the electromagnetic levitation system.

Keywords

  • active magnetic bearing
  • electromagnetic levitation system
  • motion control

1. Introduction

As a reaction wheel in spacecraft, a ball bearing wheel, a magnetic bearing wheel, and an active magnetic bearing wheel have mainly been used. First, a ball bearing wheel uses a ball bearing to maintain the separation between the bearing races. Ball bearings reduce rotational friction and support radial and axial loads by using at least two races to contain balls and transmit the loads through balls. Ball bearings tend to have a lower load capacity than other kinds of rolling element bearings mainly due to the small contact area between balls and races. Also, ball bearings should be lubricated periodically with a lubricant such as oil and grease for ball bearings to operate properly [1]. Next, a magnetic bearing is used in a magnetic bearing wheel. A magnetic bearing supports a load by the magnetic levitation principle. In magnetic bearing wheels, permanent magnets are used to carry a wheel, a control system is used to hold a wheel stable, and power is used when a levitated wheel deviates from its target position. A magnetic bearing wheel also requires a back‐up bearing in case of control system or power failure and during initial start‐up conditions. A magnetic bearing has two kinds of instabilities. One is that attractive magnets provide an unstable static force that decreases at distant distances and increases at close distances. The other is that magnetism gives rise to oscillations that may cause loss of suspension if driving forces are present [1]. Finally, in an active magnetic bearing wheel, a rotating shaft is levitated by the principle of electromagnetic suspension. A wheel is supported in an active magnetic bearing wheel without physical contact. The contactless operation of the active magnetic bearing wheels eliminates the need of lubrication of the bearing components, which allows them to operate cleanly. Moreover, it can accommodate irregularities in the mass distribution automatically, which allows it to spin around its center of mass with very low vibration, and can suppress the nutation and precession of the rotor effectively. The components of an active magnetic bearing wheel are an active magnetic bearing, a wheel, a control system, an electromagnet assembly, power amplifiers, and gap sensors. This bias current is mediated by a control system that offsets the bias current by equal but opposite perturbations of current as the rotor deviates by a small amount from its center position [1].

The active magnetic bearing wheel exhibits very lower vibration than ball bearing wheels and magnetic bearing wheels. Thus, it is a desirable reaction wheel for the spacecraft attitude control since vibration is the critical factor for the high precision spacecraft attitude control. The active magnetic bearing is the very important component among components of active magnetic bearing wheels. Due to this importance, various kinds of active magnetic bearings have been developed and their control methods have been studied (e.g., [210]).

In this chapter, the author presents an electromagnetic levitation system for active magnetic bearing wheels. This system consists of a rotor, a shaft, a cover, and a base. Also, this system does not include a mechanism for spinning the rotor around its rotating axis. The author derives a meaningful electromagnetic force by using the singular value decomposition [11]. The proportional‐integral‐derivative (PID) controller is used to control the position of the rotor and regulate the two gimbal angles of the rotor. The author gives the numerical simulation and experimental results on the control of the electromagnetic levitation system.

Advertisement

2. Electromagnetic levitation system

In Figure 1, the schematic of the electromagnetic levitation system developed in this chapter is given. The cover protects the rotor, and the base supports the rotor, shaft, and cover. This system can levitate the rotor up to 0.8 mm from the ground in the z‐axis, rotate the rotor, and gimbal the rotor within a small angle of ±0.2°.

Figure 1.

Schematic representation of the electromagnetic levitation system.

The dynamic equations of motion of the electromagnetic levitation system are given as follows:

mz ¨= FzE1
Iɸ¨ = TxE2
Iθ¨= TyE3

In Eqs. (1)(3), φ and θ are the gimbal angles of the rotor in the x‐ and y‐axes, respectively, z is the displacement of the rotor in the z‐axis, Fz is the control force in the z‐axis, Tx and Ty are the control torques applied to the rotor in the x‐ and y‐axes, respectively, m = 0.72 kg is the mass of the rotor, and I = 877.367 × 10-6 kg m2 is the inertia of the rotor for the x‐ and y‐axes.

Let us consider the four pairs of electromagnets shown in Figure 1. Then, the control inputs Fz in Eq. (1), Tx in Eq. (2), and Ty in Eq. (3) can be represented as follows:

i=14Fei = Fz + mgE4
De(Fe2  Fe4) = TxE5
De(Fe3  Fe1) = TyE6

In Eqs. (4)(6), Fei is the electromagnetic force generated by the ith pair of electromagnets, and g = 9.8 m/s2 is the acceleration of gravity. Then, Eqs. (4)(6) can be written as follows:

AFe = u + ugE7

where

A[10De 1De0 10De 1De0 ]E8

Fe[Fe1 Fe2 Fe3 Fe4]T, u [Fz Tx Ty]T, and ug [mg 0 0]T. After designing the control inputs Fz, Tx, and Ty, the four electromagnetic forces Fei, i = 1,…,4 have to be determined by Eq. (7). Among solutions for Eq. (7), the minimal norm solution is derived by using the singular value decomposition [11]. Let the singular value decomposition of the matrix AR3×4 in Eq. (8) be UΣVT and define

A++UTE9

where UR3×3 and VR4×4 are orthogonal matrices, Σ ∈ R3×4, and A+R4×3 and Σ+R4×3 denote the pseudoinverse matrices of the matrices A and Σ, respectively. With some calculations, we obtain the following by the singular value decomposition of the matrix A in Eq. (8)

V=12[1111 0202 2020 1111 ]E10
Σ+=12[1000 02De00 002De0 ]  [ Σ1101×3 ]E11
U=[100010001 ]E12

In Eq. (11), 01×3 implies the 1 × 3 zero matrix. Then, with A+ of Eq. (9), we obtain

Fe=A+(u+ug)=[14(Fz + mg)  12TyDe 14(Fz + mg) + 12TxDe14(Fz + mg) + 12TyDe14(Fz + mg)  12TxDe]E13

The following condition holds for any other solution F^e  to Eq. (7) [12]

Fe2 < F^e2E14

where ‘2 denotes the Euclidean norm of a vector (i.e., for a vector xRn, x2  i=1nxi2).

Since the two gimbal angles, φ and θ, are very small, we can approximate sin(φ) φ and sin(θ) θ by the small‐angle approximation. Thus, the displacement from the bottom surface of the ith gap sensor to the top surface of the rotor in the z‐axis can be calculated as follows:

lg1=Lg1+ δLg1=Lg1z+Dgsin(θ)Dgsin(φ)Lg1z+Dgθ DgφE15
lg2=Lg2+ δLg2=Lg2zDgsin(θ)Dgsin(φ)Lg2z Dgθ DgφE16
lg3=Lg3+ δLg3=Lg3zDgsin(θ)+Dgsin(φ)Lg2z Dgθ + DgφE17
lg4=Lg4+ δLg4=Lg4z +Dgsin(θ)+Dgsin(φ)Lg4z +Dgθ + DgφE18

Then, from Eqs. (15) to (18), the system state z, φ, and θ can be calculated as follows:

z=14(i=14Lgi  i=14lgi)E19
φ =14Dg[(lg3 + lg4  lg1  lg2)  (Lg3 + Lg4  Lg1  Lg2)]E20
θ= 14Dg[(lg1 + lg4  lg2  lg3)  (Lg1 + Lg4  Lg2  Lg3)]E21

Similarly, the displacement from the bottom surface of the ith pair of electromagnets to the top surface of the rotor in the z‐axis can be calculated as follows:

le1=Le1+ δLe1=Le1z +Desin(θ)Le1z +DeθE22
le2=Le2+ δLe2=Le2z Desin(φ)Le2z DeφE23
le3=Le3+ δLe3=Le3z Desin(θ)Le3z DeθE24
le4=Le4+ δLe4=Le4z +Desin(φ)Le4z +DeφE25

By the Maxwell's equation [13], the following equation is obtained for the control currents supplied to the coils of the four pairs of electromagnets

ii=2leinFeiµ0G , i=1,,4E26

where μ0 = 4π × 10-7 N/A2 is the permeability constant of free space, n = 240 is the number of coil turn, ii is the control current of the ith pair of electromagnets, and G = 50.265 × 10-6 m2 is the cross‐sectional area of a pair of electromagnets. The author limits each control current in Eq. (26) by 1 A.

Advertisement

3. Numerical simulation and experimental results

The author gives the numerical simulation and experimental results on the control of the electromagnetic levitation system in this section.

The author uses the following discretized PID controller to control the displacement of the rotor in Eq. (19) and the two gimbal angles of the rotor in Eqs. (20) and (21).

u(t) = Kpxe(t) + KiTsi=otxe(i) + Kd(xe(t)  xe(t1))TsE27

where Ts is the sampling time and given by Ts = 1 ms, xextx is the error between the system state x ≜ [z φ θ]T and the target system state xt ≜ [zt φt θt]T and Kp ≜ diag[Kpz, K, K], Ki ≜ diag[Kiz, K, K], and Kd ≜ diag[Kdz, K, K] denote the 3 × 3 diagonal positive definite matrices.

The target position and target gimbal angles of the rotor are set to be zt = 0.3 mm and φt = 0° and θt = 0°, respectively. The author initially decides the feedback gains of the PID controller in Eq. (27) that can achieve the control objective by adopting the well‐known Ziegler‐Nichols method [14] and then finely tunes the feedback gains of the PID controller in Eq. (27) by an experiment. Also, the antiwindup compensator is used to make the overshot as small as possible. According to the Ziegler‐Nichols method [14], first, Ki = Kd = diag[0, 0, 0] are set, and the proportional gain Kp is then increased until the system just oscillates. The proportional gain is then multiplied by 0.6, and the integral and derivative gains are calculated as Kp = 0.6 Km, Ki = Kp(ωm/π), and Kd = Kp(0.25 π/ωm) where Km ≜ diag[Kmz, K, K] denotes the 3 × 3 diagonal positive definite matrix with the gain elements at which the proportional system oscillates, and ωm diag[ωmz, ω, ω] denotes the 3 × 3 diagonal positive definite matrix with the oscillation frequency elements. As a result, the feedback gains of the PID controller in Eq. (27) are chosen as follows: Kpz = 1000, K = K = 20, Kiz = 5000, K = K = 8.33, Kdz = 100, and K = K = 0.013. Consequently, we see that the dominant feedback gains in this PID controller are Kp and Ki, and thus, one can obtain Kmz ≅ 1666.67, K = K 33.33, ωm ≅ 15.71 rad/s and ω = ω ≅ 1.31 rad/s.

The control flow diagram of the system is shown in Figure 2. After we measure the displacements from the bottom surfaces of the four gap sensors to the top surface of the rotor, we calculate the displacement of the rotor and the two gimbal angles of the rotor by Eqs. (19)(21), respectively. The control input u ≜ [Fz Tx Ty]T is made by the PID controller in Eq. (27). Then, we calculate the four electromagnetic forces and the displacements from the bottom surfaces of the four pairs of electromagnets to the top surface of the rotor by Eq. (13) and the Eq. (22)(25), respectively. After we calculate the control currents by Eq. (26), they pass through the current limiters and are supplied to the coils of the four pairs of electromagnets by the power electronics.

Figure 2.

Control flow diagram of the electromagnetic levitation system.

In the numerical simulation, it will be demonstrated that the electromagnetic force Fe of Eq. (13) satisfies the condition of Eq. (14) with respect to the following electromagnetic force F^e, which is another solution to Eq. (7)

F^e=[mg 12Fz + 12(Tx  Ty)De  12mgTyDe + mg12Fz  12(Tx  + Ty)De  12mg][F^e1F^e2F^e3F^e4]TE28

With the initial system state given by zini = 0 mm, z˙ini = 0 mm/s, the values are φini = θini = 0° and φ˙ini = θ˙ini = 0°/s. In Figures 35, the numerical simulation results on the control of the electromagnetic levitation system using Fe of Eq. (13) and F^eof Eq. (28) are shown. As shown in Figure 3, the rotor reaches the target position zt = 0.3 mm, and the control forces using Fe of Eq. (13) and F^e of Eq. (28) show the same behaviors, respectively. In Figure 4, we see that each control current using Fe of Eq. (13) reaches the same value, the control currents of i1 and i3 using F^e1 and F^e3 of Eq. (28), respectively, reach the same value, and the control currents of i2 and i4 using F^e2  and F^e4  of Eq. (28), respectively, reach the same value. In Figure 5, we see that the Euclidean norm of Fe of Eq. (13) is about 3.162 times smaller than that of F^e of Eq. (28) in the steady‐state region. Therefore, the numerical simulation results shown in Figure 5 illustrate that Fe of Eq. (13) satisfies the condition of Eq. (14) with respect to F^eof Eq. (28).

Figure 3.

Time histories of the position of the rotor and the control force, which are obtained by the numerical simulation.

Figure 4.

Time histories of the control currents, which are obtained by the numerical simulation.

Figure 5.

Time histories of the Euclidean norm of the electromagnetic force, which are obtained by the numerical simulation.

The experimental results on the control of the system using Fei, i = 1,…,4 of Eq. (13) are shown in Figures 69. As shown in Figures 6 and 7, the PID controller successfully levitates the rotor at the target position zt = 0.3 mm with well regulating the two gimbal angles. The trajectories of control currents are shown in Figure 8. And the Euclidean norm of Fe of Eq. (13) is shown in Figure 9. As shown in Figures 4, 5, 8, and 9, the trajectories of control currents and Euclidean norm of Fe of Eq. (13) obtained by the experiment move around the values obtained by the numerical simulation.

Figure 6.

Time histories of the system state of the electromagnetic levitation system using Fe of Eq. (13), which are obtained by the experiment.

Figure 7.

Time histories of the control inputs of the electromagnetic levitation system using Fe of Eq. (13), which are obtained by the experiment.

Figure 8.

Time histories of the control currents of the electromagnetic levitation system using Fe of Eq. (13), which are obtained by the experiment.

Figure 9.

Time histories of the Euclidean norm of Fe of Eq. (13), which are obtained by the experiment.

In order to simulate an external disturbance, a human hand presses down hard on the rotor to the ground at about 20 s, and it is removed from the rotor momentarily. In Figure 6, we see that, after we remove the external disturbance from the rotor, the system state becomes a steady state within 2 s. Also, in Figure 8, we see that, as the external disturbance applied to the rotor increases, each control current increases to resist the external disturbance. It should be remarked that the operating parameters like an applied load and a speed of the rotor may influence on the design of electromagnetic levitation system because these parameters make an impact on the system dynamics.

Advertisement

4. Conclusion

In this chapter, the electromagnetic levitation system was developed as a prototype for developing active magnetic bearing wheels. A control system was developed to control the position and two gimbal angles of the rotor. The experimental results demonstrated that the control system can control the position of the rotor and regulate the two gimbal angles. The refinement of the electromagnetic levitation system for the development of active magnetic bearing wheels is the further research topic.

Advertisement

Acknowledgments

The author would like to thank the Satellite Technology Research Center at the Korea Advanced Institute of Science and Technology for its support to develop the electromagnetic levitation system.

References

  1. 1. Schweitzer G., Maslen E.H., editors. Magnetic bearings. Berlin, Germany: Springer; 2009. pp. 1–535.
  2. 2. Trumper D.L., Olson S.M., Subrahmanyan P.K. Linearizing control of magnetic suspension systems. IEEE Trans. Control Syst. Technol. 1997;5(4):427–438.
  3. 3. Hong S.‐K., Langari R. Robust fuzzy control of a magnetic bearing system subject to harmonic disturbances. IEEE Trans. Control Syst. Technol. 2000;8(2):366–371.
  4. 4. Lindlau J.D., Knospe C.R. Feedback linearization of an active magnetic bearing with voltage control. IEEE Trans. Control Syst. Technol. 2002;10(1):21–31.
  5. 5. Mazenc F., de Queiroz M.S., Malisoff M., Gao F. Further results on active magnetic bearing control with input saturation. IEEE Trans. Control Syst. Technol. 2006;14(5):914–919.
  6. 6. Sivrioglu S. Adaptive backstepping for switching control active magnetic bearing system with vibrating base. IET Control Theory Appl. 2007;1(4):1054–1059.
  7. 7. Jugo J., Lizarraga I., Arredondo I. Nonlinear modelling and analysis of active magnetic bearing systems in the harmonic domain: a case study. IET Control Theory Appl. 2008;2(1):61–71.
  8. 8. Tang J., Fang J., Wen W. Superconducting magnetic bearings and active magnetic bearings in attitude control and energy storage flywhee for spacecraft. IEEE Trans. Appl. Supercond. 2012;22(6):5702109.
  9. 9. Schuhmann T., Hofmann W., Werner R. Improving operational performance of active magnetic bearings using kalman filter and state feedback control. IEEE Trans. Ind. Electron. 2012;59(2):821–829.
  10. 10. Fang J., Zheng S., Han B. Attitude sensing and dynamic decoupling based on active magnetic bearing of MSDGCMG. IEEE Trans. Instrum. Meas.. 2012;61(2):338–348.
  11. 11. Strang G. Linear Algebra and Its Applications. 3rd ed. Florida, USA: Harcourt Brace Jovanovich Publishers; 1988. pp. 1–505.
  12. 12. Golub G.H., Van Loan C.F. Matrix computations. 3rd ed. Maryland, USA: The Johns Hopkins University Press; 1996. pp. 1–728.
  13. 13. Maxwell J.C. A dynamical theory of the electromagnetic field. Phil. Trans. R. Soc. Lond. 1865;155:459–512.
  14. 14. Ogata K. Modern control engineering. 5th ed. New Jersey, USA: Prentice‐Hall; 1996. pp. 1–894.

Written By

Yonmook Park

Submitted: 28 May 2016 Reviewed: 13 December 2016 Published: 31 May 2017