1. Introduction
Magnetic bearings have been investigated extensively for the last decades. Unlike heteropolar bearings [1-2], Homopolar magnetic bearings have a unique biasing scheme that directs the bias flux flow into the active pole plane where it energizes the working air gaps, and then returns through the dead pole plane and the shaft sleeve. Use of rare earth permanent magnets such as Samarium-Cobalt (Sm-Co) and Neodymium-Iron-Boron (Nd-Fe-B) yields a very high efficiency when the permanent magnets are used as the source of bias flux to energize the air gaps and electromagnets are used to supply control fluxes in the active plane. Unique flux paths of mixed axial/radial orientation in homopolar design help to lower eddy current power losses significantly [3].
Maslen and Meeker [4] introduced a fault tolerant 8-pole heteropolar magnetic bearing actuator with independently controlled currents. Flux coupling in a heteropolar magnetic bearing allows the remaining coils to produce force resultants identical to the unfailed bearing, if the remaining coil currents are properly redistributed. Na and Palazzolo [5] also investigated the optimized realization of fault-tolerant magnetic bearing actuators and experimentally showed it on a flexible rotor such that rotor displacements after failure can be maintained close to the displacements before failure for up to all combinations of 4 coils failed and certain combinations of 5 coils failed out of 8 coils.
The present work describes the theory and following numerical analysis for the novel fault-tolerant homopolar magnetic bearing. Energy efficient homopolar magnetic bearings with fault tolerant capability may find great use in some applications such as flywheel energy storage systems and momentum wheels.
2. Bearing Model
The schematic drawing of a 6-active pole, double plane, permanent magnet biased homopolar magnetic bearing is shown in Fig. 1.
Fig. 1 Homopolar magnetic bearing
Fig. 2 Equivalent magnetic circuit for 6 pole homopolar magnetic bearing
Assuming that eddy current effects and material path reluctances are neglected, Maxwell’s equations are reduced to the equivalent magnetic circuit for the homopolar magnetic bearing as shown in Fig. 2.
The reluctance in air gap j of the active pole plane is;
\(R _ { j } = \frac { g _ { j } } { \mu _ { 0 } a _ { 0 } }\)
\(g _ { j } = g _ { 0 } - x \operatorname { cos } \theta _ { j } - y \operatorname { sin } \theta _ { j }\)
The parameters , , , , and represent the permeability of air, the pole face area of the active pole, the rotor positions, and nominal air gap, respectively. The permanent magnets are modeled as a source, , and the return path reluctance . Applying Ampere’s loop law to the magnetic circuit results in a matrix equation.
\(\)
\(R \Phi = H + N I\)
The magnetic circuit equation (4) for a homopolar magnetic bearing has a permanent magnet bias forcing term that energize the working air gaps. The coil current vector then provides only AC control currents. The fault tolerant, 6-active pole, homopolar magnetic bearing utilizes 6 independent coils each driven by its power amplifiers. The currents distributed to the 6-active-pole bearing are generally expressed as a 6×2 distribution matrix and control voltage vector . The current vector is;
\(\left. \begin{array} { l } { I = T v _ { c } } \\ { T = \left[ \begin{array} { l l } { T _ { x } } & { T _ { y } } \end{array} \right] , v _ { c } = \left[ \begin{array} { l } { v _ { c x } } \\ { v _ { c y } } \end{array} \right] } \end{array} \right. \)
The parameters and represent the x and y control voltages. A typical current distribution scheme (coil winding scheme) for a homopolar magnetic bearing is;
\(\tilde { T } = \left[ \begin{array} { c c } { \sqrt { 3 } / 2 } & { 1 / 2 } \\ { 0 } & { 1 } \\ { - \sqrt { 3 } / 2 } & { 1 / 2 } \\ { - \sqrt { 3 } / 2 } & { - 1 / 2 } \\ { 0 } & { - 1 } \\ { \sqrt { 3 } / 2 } & { - 1 / 2 } \end{array} \right]\)
The feedback control voltages and , determined with any type of control law and measured rotor motions, are distributed to each pole via in normal operation, and create effective stiffness and damping of the bearing to suspend the rotor around the bearing center position. If some coils fails,
the full (6×1) current vector is related to the reduced current vector by introducing a failure map matrix .
\(I = W I\)
For example the matrix for the 4th-6th coil failed bearing is described as;
\(W = \left[ \begin{array} { l l l l } { 1 } & { 0 } & { 0 } & { 0 } \\ { 0 } & { 1 } & { 0 } & { 0 } \\ { 0 } & { 0 } & { 1 } & { 0 } \\ { 0 } & { 0 } & { 0 } & { 0 } \\ { 0 } & { 0 } & { 0 } & { 1 } \\ { 0 } & { 0 } & { 0 } & { 0 } \end{array} \right. \)
The reduced distribution matrix is defined as;
\(\hat { I } = \hat { T } v _ { c }\)
\(\left. \begin{array} { l } { \hat { T } = [ \hat { T } _ { x } \hat { T } _ { y } ] } \\ { \hat { T _ { x } } = [ t _ { 1 } t _ { 2 } t _ { 3 } ] ^ { T } , \hat { T } _ { y } = [ t _ { 4 } t _ { 5 } t _ { 6 } ] ^ { T } } \end{array} \right. \)
The reduced distribution matrix is to be determined such that the magnetic forces should remain very much invariant before and after failure. The flux density vector can be determined from the magnetic circuit equation in Eq. (3). Leakage and fringing effects can be empirically determined and are simply derated by the fringing factor . The flux density vector in the air gaps is described as;
\(B = \zeta A ^ { - 1 } R ^ { - 1 } ( H + N I )\)
where the pole face area matrix is . The flux density vector is then reformulated as;
\(B = G v\)
\(\left. \begin{array} { l } { G = [ G _ { b } H G _ { c } \hat { T } _ { x } G _ { c } \hat { T } _ { y } ] \cdot v = [ 1 v _ { c x } v _ { c y } ] ] ^ { T } } \\ { G _ { b } = \zeta A ^ { - 1 } R ^ { - 1 } , G _ { c } = \zeta A ^ { - 1 } R ^ { - 1 } N W } \end{array} \right. \)
Magnetic forces developed in the active pole plane are then described as;
\(f _ { x } = v ^ { T } M _ { x } v\)
\(f _ { y } = v ^ { T } M _ { y } v\)
\(M _ { x } ( \hat { T } ) = - G ^ { T } \frac { \partial D } { \partial x } G , M _ { y } ( \hat { T } ) = - G ^ { T } \frac { \partial D } { \partial y } G\)
and where the air gap energy matrix is;
\(D = \operatorname { diag } ( [ g _ { j } a _ { 0 } / ( 2 \mu _ { 0 } ) ] )\)
Employing an optimal current distribution matrix may decouple the linearized forces of the failed bearing, and even maintain the same decoupled magnetic forces as those of an unfailed magnetic bearing. The necessary conditions to yield the same decoupled magnetic control forces are;
\(M _ { x } = k _ { v } \left[ \begin{array} { c c c } { 0 } & { \frac { 1 } { 2 } } & { 0 } \\ { \frac { 1 } { 2 } } & { 0 } & { 0 } \\ { 0 } & { 0 } & { 0 } \end{array} \right] , M _ { y } = k _ { v } \left[ \begin{array} { c c c } { 0 } & { 0 } & { \frac { 1 } { 2 } } \\ { 0 } & { 0 } & { 0 } \\ { \frac { 1 } { 2 } } & { 0 } & { 0 } \end{array} \right]\)
If the distribution matrix is determined such that Eq. (15) satisfies Eq. (17), the nonlinear magnetic forces in Eqs. (13) and (14) become linearized at bearing center positions.
\(f _ { x } = k _ { v } v _ { c x } , f _ { y } = k _ { v } v _ { c y }\)
Equations (15) and (17) can be written in 18 scalar forms, and then boils down to 10 algebraic equations if redundant terms are eliminated. The equality constraints to yield the same control forces before and after failure are;
\(\left. \begin{array} { l } { h _ { 1 } ( \hat { T } ) = \hat { T } _ { x } ^ { T } Q _ { x 0 } \hat { T } _ { x } = 0 } \\ { h _ { 2 } ( \hat { T } ) = \hat { T } _ { y } ^ { T } Q _ { x 0 } \hat { T } _ { y } = 0 } \\ { h _ { 3 } ( \hat { T } ) = H ^ { T } Q _ { b x 0 } \hat { T } _ { y } = 0 } \\ { h _ { 4 } ( \hat { T } ) = \hat { T _ { x } ^ { T } } Q _ { x 0 } \hat { T } _ { y } = 0 } \\ { h _ { 5 } ( \hat { T } ) = H ^ { T } Q _ { b x 0 } \hat { T } _ { y } = k _ { v } / 2 } \\ { h _ { 6 } ( \hat { T } ) = \hat { T _ { x } ^ { T } } Q _ { y 0 } \hat { T } _ { x } = 0 } \\ { h _ { 7 } ( \hat { T } ) = \hat { T _ { y } ^ { T } } Q _ { y 0 } \hat { T } _ { y } = 0 } \\ { h _ { 8 } ( \hat { T } ) = H ^ { T } Q _ { b y 0 } \hat { T } _ { x } = 0 } \\ { h _ { 9 } ( \hat { T } ) = \hat { T } _ { x } ^ { T } Q _ { y 0 } \hat { T } _ { y } = 0 } \\ { h _ { 10 } ( \hat { T } ) = H ^ { T } Q _ { b y 0 } \hat { T } _ { y } = k _ { v } / 2 } \end{array} \right. \)
where
\(\)
The criterion for choosing the best candidate is the one that will yield the maximum load capacity prior to any saturation. To accomplish this a distribution matrix can be determined
by using the Lagrange Multiplier method to minimize the Euclidean norm of the flux density vector . The cost function is defined as;
\(J = B ( \hat { T } ) ^ { T } P B ( \hat { T } )\)
where the diagonal weighting matrix is also selected to maximize the load capacity. The Lagrange Multiplier method is then used to solve for that satisfies Eq. (19). Define:
\(L ( \hat { T } ) = B ( \hat { T } ) ^ { T } P B ( \hat { T } ) + \sum _ { j = 1 } ^ { 10 } \lambda _ { j } h _ { j } ( \hat { T } )\)
Partial differentiation of Eq. (21) with respect to and leads to 16 nonlinear algebraic equations to solve for and .
\(\Psi = \left[ \begin{array} { c } { \Psi _ { 1 } ( t , \lambda ) } \\ { \Psi _ { 2 } ( t , \lambda ) } \\ { \cdot } \\ { \cdot } \\ { \Psi _ { 15 } ( t , \lambda ) } \\ { \Psi _ { 16 } ( t , \lambda ) } \end{array} \right] = 0\)
\(\Psi _ { i } = \frac { \partial L } { \partial t _ { i } } = 0 , i = 1,2 , \ldots , 6\)
\(\Psi _ { 6 + j } = h _ { j } ( \hat { T } ) = 0 , j = 1,2 , \ldots , 10\)
3. Numerical Analysis
Various initial guesses are tested to find the solution of Eq.(22). Some examples of distribution matrices are calculated for the 6-pole homopolar magnetic bearing with the nominal air gap (0.508mm), pole face area (602mm2), number of coil turns (50turns). It is assumed that permanent magnets are selected to produce bias flux density of 0.8 Tesla in the air gaps of the active pole plane. A distribution matrix for an 6-pole homopolar bearing with the 5th-6th coils failed operation is calculated as;
\(T _ { 56 } = \left[ \begin{array} { c c } { 0.9669 } & { 1.6434 } \\ { - 1.1149 } & { 1.3394 } \\ { - 0.6164 } & { 1.6591 } \\ { - 1.8801 } & { - 0.0161 } \\ { 0 } & { 0 } \\ { 0 } & { 0 } \end{array} \right]\)
A distribution matrix with the 4th-5th-6th coils failed operation is calculated as;
\(T _ { 456 } = \left[ \begin{array} { c c } { 1.6776 } & { 1.0501 } \\ { 0.0003 } & { 1.9262 } \\ { - 1.6781 } & { 1.0501 } \\ { 0 } & { 0 } \\ { 0 } & { 0 } \\ { 0 } & { 0 } \end{array} \right]\)
Similarly, the distribution matrices can be calculated for a failed homopolar bearing up to any combinations of 3 coils failed out of 6 coils.
4. Control Simulations
The schematic of a fault-tolerant control scheme utilizing distribution matrices developed in the previous section is shown in Fig. 3.
The controller consists of two independent parts, which are a feedback voltage control law and an adaptive current distribution mechanism. Though any control algorithm for magnetic bearing systems appearing in the literature can be utilized with the fault tolerant scheme, for sake of illustration, a simple PD feedback control law is used to stabilize the system. While the feedback control law remains unaltered during the failure the appropriate current distribution matrix can be continuously updated using an adaptive current distribution mechanism. The fault-tolerant control scheme can be easily implemented in a physical controller (DSP). The distribution matrix of T~ is switched to T56 and T456 when 2 adjacent coils failed at 0.05 seconds and then 3 adjacent coils failed at 0.1 seconds. The rotor speed is held constant at 10,000 RPM.
Fig. 3 Schematic of fault-tolerant control system
Fig. 4 Current plot for a series of failures
Fig. 5 Flux density plot for a series of failures
Fig. 6 Displacement plot for a series of failures
Fig. 7 Orbit plot for a series of failures
Fig. 4 shows transient response of the current inputs to the outboard bearing from the normal unfailed operation through the 5-6th coils and 4-5-6th coils of the outboard bearing failed at 0.05 seconds and 0.1 seconds, respectively. Fig. 5 shows the corresponding transient response of the flux densities in the outboard bearing for the 5-6th coils and 4-5-6th coils failed operation.
Fig. 6 and Fig. 7 show transient rotor displacements and orbits at the outboard bearing for the 5-6th coils and 4-5-6th coils failed operation respectively. Figs. 4-7 indicates that very much the same rotordynamic responses are maintained throughout the series of failure events, while currents and fluxes in the homopolar magnetic bearing change significantly.
5. Conclusion
A fault-tolerant control scheme is developed for the permanent magnet biased, 6-active- pole, homopolar magnetic bearings.
The homopolar bearing actuator using the fault tolerant control algorithm can preserve the same linearized magnetic forces by redistributing the remaining currents even if some components such as coils or power amplifiers suddenly fail. The distribution matrix T of control voltages is determined by using the Lagrange Multiplier optimization with equality constraints for a failed bearing in a manner that the load capacity should be maximized. Simulations show that very much the same vibrations (orbits or displacements) are maintained throughout failure events while currents and fluxes change significantly with different distribution scheme. The distribution matrices can be calculated for a failed homopolar bearing up to all combinations of 3 coils failed out of 6 coils. The fault tolerance of homopolar magnetic bearings is achieved at the expense of additional hardware requirements such as independent coils (power amplifiers), fault detection system, additional DSP controller channels.
Acknowledgements
This paper was supported by Kyungnam University Research Fund.
참고문헌
- P. E. Allaire, D. W. Lewis, J. D. Knight, "Active Vibration Control of a Single Mass Rotor on Flexible Supports," Journal of the Franklin Institute vol. 315 pp. 211-222, (1983). https://doi.org/10.1016/0016-0032(83)90025-X
- J. Salm, G. Schweitzer, "Modeling and Control of a Flexible Rotor with Magnetic Bearing," Proceedings of the Third International Conference on Vibrations in Rotating Machinery pp. 553-561, (1984).
- Y. Fan, A. Lee, F. Hsiao, "Design of a Permanent/Electromagnetic Magnetic Bearing-Controlled Rotor System." Journal of the Franklin Institute, vol. 334B, pp. 337-356, (1997). https://doi.org/10.1016/S0016-0032(96)00101-9
- E. H. Maslen, D. C. Meeker, "Fault To lerance of Magnetic Bearings by Generalized Bias Current Linearization." IEEE Transactions on Magnetics, vol. 31, pp. 2304-2314, (1995). https://doi.org/10.1109/20.376229
- U. J., Na, A. B., Palazzolo, and A., Provenza, "Test and Theory Correlation Study for a Flexible Rotor on Fault-Tolerant Magnetic Bearings." ASME Journal of Vibration and Acoustics, vol. 124, pp. 359-366, (2002). https://doi.org/10.1115/1.1467652