Abstract

This paper analyzes the effects of time delay on the stability of the rotation modes for the magnetically suspended flywheel (MSFW) with strong gyroscopic effects. A multi-input multioutput system is converted into a single-input single-output control system with complex coefficient by variable reconstruction, and the stability equivalence of the systems before and after variable reconstruction is proven. For the rotation modes, the stability limits and corresponding vibration frequencies are found as a function of nondimensional magnetic stiffness and damping and nondimensional parameters of rotor speed and time delay. Additionally, the relationship between cross feedback control system stability and time delay is investigated. And an effective phase compensation method based on cross-channel is further presented. Simulation and experimental results are presented to demonstrate the correctness of the stability analysis method and the superiority of the phase compensation strategy.

1. Introduction

With the development of magnetic bearing technology, a magnetically suspended flywheel (MSFW) represents an exciting alternative to the traditional mechanical flywheel due to its inherent superior features such as contact free operation, zero friction, adjustable damping, and stiffness as well as the potential for low vibration and super high rotational speeds [1, 2]. As a result, it has been widely employed in power industry [3, 4], transportation, spaceflight [58], and so on [9, 10].

To maximize the energy storage capacity, a MSFW usually chooses a radially thick rotor so that the principle mass moment of inertial is larger than the transversal counterpart, that is, . In this case, the gyroscopic effects are especially strong, which puts a strain on the stability analysis of the system.

Particularly, to improve the transient response and to reduce power consumption of a magnetic bearing, the permanent-magnetic-biased hybrid magnetic bearings (PMHMBs) are applied [1]. This bias field causes the rotor to behave as it was supported on a spring with negative stiffness. There is no time delay associated with the bias field since it is caused by a permanent magnet. However, the time delay of the electromagnetic bearing always exists since each of the components involved in it and its control system have a time delay associated with them. For an analog device, the time delay can be approximated by a linear time delay low-pass filter. The amplitude of the transfer function is assumed to be a constant and the phase is linear. However, for the digital devices, this time delay is the loop time for each measurement [11]. The total time delay is the sum of the individual time delay [12]. All of these make the stability analysis of the MSFW more complex and difficult.

In recent years, considerable researches have been carried out for the stability analysis of the rotor supported by magnetic bearings with time delay. Ji [13, 14] analyzed the effect of time delays on the nonlinear dynamic behavior of the magnetic bearing system. Wang and Liu [15], Li et al. [16], and Ji [17] studied the Hopf bifurcation and its stability in a magnetic bearing system with time delays. However, these works only considered a simple single-degree-of-freedom (DOF) magnetic bearing without gyroscopic effects. As for the multi-DOF cases with high rotor speed and strong gyroscopic effects, Kascak et al. [18] employed nondimensional magnetic stiffness and damping as well as nondimensional system parameters to describe the stability limits of a PD controller. However, they did not consider the cross feedback controller. Kascak et al. [19] presented the stability boundaries of a spinning rotor with parametrically excited gyroscopic system. Bauomy [20] performed the stability analysis of a rotor active magnetic bearing system with time varying stiffness. Inayat-Hussain [21] studied the nonlinear dynamics of a magnetically supported rigid rotor in auxiliary bearings. Andersen et al. [22] and Sugai et al. [23] analyzed the rotation modes stability of the magnetically suspended rotor with passive cylinder-magnet bearings. Additionally, Pei [24] presented the stability boundaries of a spinning rotor with parametrically excited gyroscopic system. Fan and Pan [25] described the stability of rotor-shaft system including electromagnetic exciters. All in all, few literatures studied the multi-DOF magnetic bearings with time delays and gyroscopic effects.

To resolve the gyroscopic effects, considerable research has been conducted and various approaches have been proposed [2632]. Among these methods, decentralized PD plus filtered cross feedback control is the most well-known and easy-to-implement approach [29]. To further improve the system nutation stability, different phase compensation methods are developed, such as unsymmetrical current sampling resistance method [11, 12] and inverse system method with phase compensation [33, 34], where the phase compensations are performed with both decentralized and cross-channels. For simplicity, it is called the combinational compensation method in this paper. In fact, there are other phase compensation strategies in addition to this, such as phase compensation with decentralized channels and phase compensation with cross-channels. Thus, the question naturally raised is what is the most effective way for the system stability among the three phase compensation methods? This has not been effectively resolved.

To resolve the problems above, this paper describes a modeling method which is used to theoretically map the stability boundaries of the MSFW controller with cross feedback and time delay. As for the rotation modes, the relationships between the time delays and the critical nutation frequency are further analytically analyzed, and then an effective phase compensation strategy is proposed.

This paper is organized as follows. First, in Section 2, the system is modeled with complex coefficient. Then, its rotation modes stability is further studied in Section 3. Simulation and experiments are developed in Section 4. Finally, Section 5 concludes this paper.

2. Modeling of the MSFW System and Its Variable Reconstruction

2.1. Modeling of the MSFW System

Here, we review the model of a rigid rotor supported by magnetic bearings, which has been established in [7]. Figure 1 shows the coordinate system definition of a MSFW, where is the geometric center of the magnetic bearing stator; refers to the distance between the central point of radial magnetic bearing and point ; is the rotor speed, -, -, and -axis form the generalized coordinate system of the rotor position; and are the torques in the and directions, respectively; and are the rotor angular displacements about the - and -axis. , , , , , and are the magnetic forces along the magnetic bearing coordinate system -, -, -, -, -, and -axis.

Then, the dynamic model of a rigid rotor supported by PMHMBs can be described aswhere is the mass of the flywheel rotor; , , and are the moments of inertia of the rotor about the -, -, and -axis, respectively, and ; and , , and are the linear displacements of the mass center of the rotor from point in the -, -, and -axis, respectively.

Figure 2 shows the schematic of the cross feedback control system, where , , and are the transfer functions of the basic controller (such as a PID controller), the cross feedback controller, and the antialias filter, respectively.

As for an axisymmetric rotor supported by PMHMBs, the magnetic force can be linearized as [7]where and denote, respectively, the current stiffness and displacement stiffness of the channel. Note that the four radial channels of the MSFW are often symmetric; for simplicity, and denote the current and displacement stiffness of the radial magnetic bearings, respectively.

For PMHMBs, there is no time delay associated with the bias field since it is produced by a permanent magnet [24], and the delay only exists in the electromagnetic field. Accordingly, the closed-loop model of the control system with the cross feedback controller can be given by [33, 35]where denotes the power amplifier of the HMBs, is the proportional gain of the displacement sensor, is the distance from the geometric center of the rotor to the sensor, and is the control system delay.

It can be seen from (3) that the cross feedback control only influences the system rotation modes stability and does not endanger the system translation modes stability. As for the translation modes of the MSFW, their stability analysis can be equivalent to that of the single-DOF magnetic bearing control system, which has been clearly analyzed. Accordingly, this paper focuses on the rotation modes stability analysis of the MSFW system with cross feedback control system.

The rotation modes model can be described as

2.2. Variable Reconstruction and Its Stability Equivalence Analysis

As for the rotation modes model (4), it is a two-input and two-output control system. To simplify the stability analysis of the MSFW rotation modes, variable reconstruction is employed in this paper.

From the coordinate definition of the MSFW shown in Figure 1, leads 90 degrees. Accordingly, we can define , where is the imaginary number unit and . Note that the dissymmetry of (4) and multiplying the second equation by and then adding the result to the fourth equation yield the differential equations:where .

According to the control system theory, the necessary and sufficient conditions of the rotation modes stability are as follows:

Accordingly,

As for complex-coefficient system (5), the necessary and sufficient condition of the stable control system is

Note that ; that is, and .

Therefore, the stability of system (4) is equivalent to that of (5). Then, a two-input two-output system is converted into a single-input single-output control system with complex coefficient. That is, we can get the stability of system (4) by analyzing the stability of system (5).

3. Stability Analysis of Rotation Modes of the MSFW

3.1. Stability Analysis of the Rotation Modes

According to the classical small signal stability analysis method, assuming an eigenvalue solution of the form , then

Substituting (9) into (5) yields the characteristic equation:

Define the nondimensional parameters as follows:Then, the eigenvalue equation can be simplified as follows:

The eigenvalue equation has the similar solution to that of the translation mode under the condition the rotor speed of the MSFW is zero. If the eigenvalue is complex, then the imaginary part of the eigenvalue is the circular frequency of the vibration that multiplies the time delay; that is

Substituting (13) into (12) yields

Solving the real and imaginary portions of the eigenvalue equation, the nondimension damping obtainsAnd the nondimensional stiffness achieves

If the real part of the eigenvalue is positive, the vibrations grow in time and the system is unstable. Contrarily, if is negative, the vibrations decay in time and the system is stable. The case that equals zero defines the stability boundary as follows:

Note that this paper mainly focuses on the nutation modes stability analysis and phase compensation for MSFW with cross feedback controller and time delay, the possible hidden oscillations [36], Bifurcation and global stability analysis being the next topics in the near future.

3.2. Relationship between Cross Feedback Control System Stability and Time Delay

To further analyze how the control system delay influences the nutation stability, the following analysis has been developed, and an effective phase compensation strategy is proposed.

The eigenvalue equation of (5) can be resolved as

Suppose its eigenvalue is ; then, substituting this into (18) yields(1)If , then we have . The system has a simple steady state bifurcation.(2)If , the system is critically stable; (19) can be simplified as

Moving the terms of trigonometric functions into the left-hand side and then squaring and adding these two equations yield a quadratic equation for :

When , the above equation can be further simplified as

Sumptuously, multiplying and to the two equations of (19) and then adding them together yield

Note that as, for a MSFW control system without phase compensation, (23) can be rewritten as

Accordingly, the second equation of (20) achieves

Note that, for the MSFW with high rotor speed, . As for positive , we haveThat is, .

Combining (24) and (26) yields

Combing (25) and (26) obtainsDeriving the above equation with yields

Note that and in practice; therefore, . That is, the larger the delay time is, the smaller the critical nutation frequency is.

Theorem 1. If the rotor speed is limited to zero, , and if the rotor speed is high enough, .

Proof. Multiplying the first and the second item of (25) by and , respectively, and then adding them together yieldsThen,whereTherefore, Hence,When , we haveNoting that yieldsSubstituting (36) into (34) yieldsObviously, if the rotor speed is limited to zero, , compared with the cross-channel, the decentralized channel is the main factor which endangers the system nutation stability. If the rotor speed is high enough, , the delay time of the cross-channel is the main factor which influences the system nutation stability. This completes the proof.

According to Theorem 1, we can further draw the following corollary.

Corollary 2. As for the control system (23), if , then the stability of the rotation modes of the MSFW has nothing to do with the rotor speed.

Proof. Since , (24) can be simplified as The stability analysis is similar to that of the translation modes of the MSFW. This completes the proof.

Accordingly, as for the MSFW system with high rotor speed, it is the delay time of the cross-channel that influences seriously the system nutation stability instead of the decentralized channel. As a result, cross-channel phase compensation is more effective than that of the decentralized channel at high speeds and the nutation stability of the MSFW can be further improved.

Similarly, it can be proven that the cross-channel compensation method is more effective than that of the combinational compensation method. Accordingly, at high rotor speeds, cross-channel compensation is the most effective phase compensation method of the three compensation methods.

4. Simulation and Experimental Study

In order to demonstrate the effectiveness and validity of the proposed stability analysis and phase compensation strategy, and to reveal how closely the theory represents the physical system, both simulation and experiments have been developed.

4.1. Simulation and Experimental Setup

The large magnetically suspended control moment gyroscope (MSCMG), as shown in Figure 3, is employed as the simulation and experimental setup since its MSFW has strong gyroscopic effects.

The well-known decentralized PD plus filter cross feedback controller is employed. The main parameters of the MSCMG and the controllers used in simulation and experiments have been listed in Tables 1 and 2. In Table 1, is the air gap of the radial magnetic bearings; and are the current stiffness and displacement stiffness of the radial magnetic bearings under the conditions and ; and and are the coil resistance and inductance of the radial magnetic bearings. In Table 2, and are the proportional and feedback coefficients of the current controller for the MSFW system, which are the same as those of the proposed controller; and are the proportional, integral, and differential coefficients of the PD plus cross feedback controller for the MSFW system; and are the cross coefficients of the high-pass and low-pass filters, respectively; is the total cross coefficient; and and are the cut-off frequencies of the two-order high-pass and low-pass filters.

Considering the bandwidths of interest of the current loop and the position loop of the MSFW and the simplicity of the realization, second-order antialias filters are introduced into both the displacement and current AD sampling of the MSFW system, whose transfer function is given as

As far as the realization and implementation in hardware are concerned, a TMS320C32 digital signal processor (DSP) is employed. Both the sampling time and servo time are set to 150 μs.

4.2. Simulation Results

Figure 4 shows a plot of the nondimensional damping versus nondimensional stiffness for different growth rate. Gyroscopic effects result in either nutation (forward whirl) or precession (backward whirl). It also shows the growth rate, (real part of the eigenvalue), and the frequency of oscillation, (imaginary part of the eigenvalue). The nondimensional stiffness is given by (16) and the nondimensional damping is given by (15). If the real part of the eigenvalue is defined as positive, the vibrations grow in time and the system is unstable. If is defined as negative, the vibrations decay in time and the system is stable. The condition that equals zero defines the stability boundary. The imaginary part of the eigenvalue, , is a parameter which varies from 0 to . The nondimensional speed of the flywheel is 0.1 for nutation, 0 for zero rotor speed, and −0.1 for precession. Nondimensional negative stiffness of 0.2 is assumed. From Figure 4, it can be drawn that a stable region is similar to that of the single-DOF magnetic bearing except it is now limited by nutation.

Figure 5 shows the dynamic stability map for a MSFW with time delay for various values of rotor speed, . Note that the stiffness must be greater than the negative stiffness for the system to be stable. The area from the trajectory to the left line () represents the stability limits resulting from the forward whirl (the real part of the eigenvalue, , equals zero). The stable region becomes smaller for higher nondimensional rotor speeds. As for the traditional decentralized PD control, there is no stable region for nondimensional rotor speeds higher than 0.9. Contrarily, as for the decentralized PD plus cross feedback control, the nondimensional rotor speeds can reach about 1.7, and there is larger stable region than that of the case without cross feedback control at every rotor speed. This demonstrates that the cross feedback control algorithm can effectively improve the system stability, which is in accordance with the analysis made above.

To further prove the correctness of Theorem 1, that is, the phase compensation method, the positive and negative frequency characteristics analysis method [33] is employed. At Hz, three different phase compensation methods have been compared, that is, the decentralized channels compensation, the cross-channel compensation, and the combinational compensation. To be fair, the three methods employ the same compensation filters:

The control channel transfer function of the MSFW system without phase compensation can be described as

According to the positive and negative frequency characteristics stability criterion [33], the control channel transfer function with the decentralized channel can be described as

Similarly, the control channel transfer function with the cross-channel compensation can be described asAnd the control channel transfer function with the combinational phase compensation can be described as

Figure 5 shows the comparative simulation results among different compensation methods, where “without compensation” demonstrates that the control system does not introduce any phase compensation. According to the nutation and precession stability criterion of the magnetically suspended with strong gyroscopic effects [33], at rotor speed , the rotor speed stability margin can be described aswhere is the gain of the control channel at the frequency .

From Figure 6, the critical whirling frequencies and their corresponding gain can be obtained: , , , , , , , and . According to the nutation stability criterion [33], it can be resolved that Hz, Hz, Hz, and Hz. It is obvious that, among the three phase compensation methods, the cross-channel phase compensation improves the system nutation stability most greatly. It is necessary to note that the traditional combinational compensation method does have some effect in improving the system nutation stability. However, it is not the most effective way compared with the cross-channel compensation. Additionally, the decentralized channel compensation method endangers the system nutation stability margin instead of improving it. All in all, among the three compensation methods, the cross-channel compensation is the most effective way, which is in accordance with the analysis made above.

4.3. Experimental Results

To further prove the effectiveness and superiority of the phase compensation method, a comparison of experiments between the cross-channel compensation and the decentralized channels compensation has been performed. At , these two different compensation methods are performed. Figures 7 and 8 show the relationship between the nutation amplitude and the rotor speed with the two different compensation methods, respectively.

As shown in Figure 7, the nutation amplitude increases heavily with the rotor speed with the decentralized compensation method. And it increases to about −49 dB at 45 Hz, which demonstrates that the system will lose its nutation stability if the rotor speed continues to increase. On the contrary, the nutation amplitude keeps at about −54 dB in the rotor speed of 40 to 70 Hz with the cross-channel compensation method as shown in Figure 8. This shows that the system nutation stability has been improved greatly by the cross-channel compensation instead of the decentralized channel compensation.

5. Conclusion

This paper studied the effects of time delay on the stability of the rotation modes and presented an effective phase compensation method based on the cross-channels for the MSFW with strong gyroscopic effects. Analysis, simulation, and experimental results yield the following conclusion:(1)The larger the rotor speed is, the smaller the stability region is.(2)The forward whirl limited the region of stability.(3)As for high speed MSFW, cross feedback control can effectively reject the nutation and precession instability.(4)The delay of the cross-channel is the main factor influencing the nutation stability of the high-speed MSFW system, and the cross-channel compensation is the most effective way to improve the nutation stability.

To sum up, the presented local stability analysis method can effectively demonstrate the system stability region of the MSFW with strong gyroscopic effects and time delay, and the presented phase compensation method can greatly improve the system nutation stability. The global stability and semistable oscillations are the next topics in the near future.

Competing Interests

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

Acknowledgments

This work was supported by the National Natural Science Foundation of China under Grants 51475472 and 61403396.