Research Article  Open Access
Shaolin Ran, Yefa Hu, Huachun Wu, Xin Cheng, "Resonance Vibration Control for AMB Flexible Rotor System Based on µSynthesis Controller", Mathematical Problems in Engineering, vol. 2018, Article ID 4362101, 16 pages, 2018. https://doi.org/10.1155/2018/4362101
Resonance Vibration Control for AMB Flexible Rotor System Based on µSynthesis Controller
Abstract
The resonance vibration control of flexible rotor supported on active magnetic bearings (AMB) is a challenging issue in the industrial applications. This work addresses the application of robust control method to the resonance vibration control for AMB flexible rotor while passing through the critical speed. This modelbased method shows great superiority to handling flexible mode vibration, which can guarantee robust stability and performance when encountering modal perturbation. First, the designed flexible rotorAMB test rig is briefly introduced. Then the system modeling is described in detail including flexible rotor, power amplifier, displacement sensors and magnetic actuator and rotordynamics are analyzed. Model validation is carried out by sine sweeping test. Finally, the μsynthesis controller is designed. The simulation and experimental results indicate that the designed μsynthesis controller, which shows great robustness to modal perturbation, can effectively suppress the resonance vibration of flexible rotor and achieve supercritical operation.
1. Introduction
Active magnetic bearing (AMB) is a typically mechatronical bearing which can make the rotor suspend stably through active feedback control. It has many characters such as contactless, nofriction, and nolubrication [1]. Benefiting from these characters, AMBs enable the machines with a very high rotational speed and has been widely used in compressors [2], motors [3], and flywheel system [4]. High rotational speed brings about many benefits such as high energy efficiency and compactness for the energysaving and emissionreduction [5]. In these cases, because of structure limitation, the high speed machines have to operate above the bending critical speeds [6]. However, the flexible rotor encounters resonance before reaching its rated speed. Resonance vibration can cause serious damage to the device. Hence, the resonance vibration suppression of flexible rotor becomes a challenging issue.
Compared with general control of flexible space structure, such as flight control [7] and robotic [8] and hard disk control [9], control of rotor’s flexible mode in AMB system can be a more challenging issue, because the AMB is inherently openloop unstable system, which has fundamental limitations for closedloop stability. Subject to gyroscopic effect and rotor internal damping, rotor dynamic is inherently speed dependent, which complicates the controller design [10]. Meanwhile, due to the nonlinear material properties in laminations, shrink fits, assembling and manufacture error, the rotor modal perturbation should be considered during controller design [11]. Therefore, the resonance vibration control of flexible rotor has been the focus of much research over the past decades.
PID is the most widely used control method in the industry for its simple structure. So, researchers resort to PID controller integrated with different general second order filters to achieve good resonance vibration control effect for AMB flexible rotor [12–14]. In these literatures, the optimal damping and optimal phase angle are derived in detail based on the plant phase information nearby the rotor bending mode, then phase lead filter and notch filter are connected in series after the PID controller based on the optimal phase angle. The added filters extend the application of PID for flexile mode control, which also complicate the control system nevertheless. This standard control method does not allow for taking into account all factors that can influence the static and dynamic quality of controlled process [15].
The robust control method, such as H∞ and μsynthesis, can take direct account of system uncertainties in the controller design procedure for actively reducing vibration of rotor systems [16]. The early research focuses on the constructions and solving procedure of the robust control problem [17–19] and mainly focuses on the rigid rotor [20–22]. Then, many scholars turn to flexible rotor robust control based on finite element method modeling. H∞ control method is successfully applied to the AMB flexible rotor system and the analytical and experimental results show the effectiveness of the H∞ control [23, 24]. However, the H∞ control is based on unstructured uncertainties. When encountered with parametric uncertainties, the synthesized controller may be conservative. Compared to H∞ control, the μsynthesis can address explicitly the parametric uncertainty in system and specify the performance specification to reduce the controller conservatism. Jinxiang Zhou et al. applied μsynthesis control method to the AMB system for dealing with parameters uncertainty of linear electromagnetic force model and weighting functions are introduced based on the objectives of AMBs performance. The experimental results showed desired performance of the AMB with μ controller [25]. In [26], rotor instability phenomenon affected by aerodynamic crosscoupled stiffness forces is studied by designing a rotor test rig with magnetic bearings and a μsynthesis controller was designed. The research results indicate that μsynthesis controller can achieve better control effect on the stability of rotor 1st bending mode under crosscoupled stiffness excitation. Although μsynthesis controller can promote the performance of AMB system on the rotor resonance vibration, some degree of system modeling accuracy must be guaranteed.
In this paper, a transparent procedure of control oriented rotorAMB system modeling and synthesis controller design is presented. Both the simulation and experimental results demonstrate that the designed controller has better rotor resonance vibration suppression performance and high robustness to system uncertainty to make the flexible rotor exceed the first bending critical speed smoothly.
2. The Flexible RotorAMB Test Rig
The purpose of the flexible rotorAMB rest rig is to study rotor vibration control and supercritical operation as shown in Figure 1.
Two radial magnetic bearings are placed at both ends of the rotor. An exciter AMB in quarter spans of the rotor is used to exert different excitation on the rotor. The rotor is connected with a 4 kW motor through a flexible coupling. Four eddy current displacement sensors are installed besides two bearings respectively to measure the displacement of the rotor in and directions. When the control is off, the rotor is down on the auxiliary bearings.
The rotor is 0.93 m long and weighs around 18.14 kg with three laminated steel journals. Two radial support AMBs are at the nondriven end (NDE) and driven end (DE). Only the NDE and DE support AMBs are utilized for control. The sketch of the flexible rotor in this work is shown in Figure 2.
3. Modeling of AMB System
μSynthesis is fundamentally a modelbased control method requiring a detailed theoretical model of the system. Therefore, the modeling of the AMB system is very important. The AMB system includes controller, power amplifier, rotor and stator, and displacement sensors. A typical control diagram of the magnetic bearing system is shown in Figure 3, which shows the signal exchange of each electromechanical component. In order to get a modelbased controller with better performance, a relatively accurate system model should be obtained.
3.1. Modeling of the AMB
The support AMB stators are 8pole heteropolar bearings with laminated siliconiron. The orientation of the quarter magnet is offset 45° from vertical, allowing the rotor weight to be distributed equally between the two control channels. This ensures the and dynamic will be nearly identical. The stator is shown in Figure 4.
(a)
(b)
Affected by current and gap, the magnetic force is inherently nonlinear. In order to simplify the control design, an approximate linear model can be as follows:where is current stiffness and is negative displacement stiffness. and only depend on the physical characteristics.
There will always be discrepancy between the actual bearing performance and analytical model due to the errors in the manufacture and linearization process. Therefore, in order to obtain the exact actuator performance, the current stiffness and displacement stiffness are determined by test data which contribute to the system modeling. When the rotor is suspended by PID controller, in direction of the rotor, the control current will be adjusted to keep the rotor at the equilibrium position under a force exerted by the pull meter. In this way, a series of forces and currents data will be recorded, and the relationship between forces and currents is obtained as shown in Figure 5.
Then, when the pull meter is removed and the rotor is suspended at different positions in the direction, a series of displacements and currents data is obtained, and their relationship is shown in Figure 6.
The slope of the forcecurrent curve is the current stiffness , which can be obtained by linear fitting. The displacement stiffness can be derived as [27]where is the slope of the currentdisplacement curve, which can also be obtained linear fitting.
Finally, = 213.5 N/A, = 1110.2 N/mm.
3.2. Modeling and Analysis of Flexible Rotor
In this work, the finite element (FE) method is employed to model the flexible rotor. As shown in Figure 7, the rotor is divided into 56 units with 57 nodes. Based on the Timoshenko beam theory, the motion equation of the rotor is as follows:where , , , and denote the mass matrix, internal shaft damping matrix, gyroscopic matrix, and stiffness matrix, respectively. They are all symmetric and positive definite except for that the gyroscopic matrix is skewsymmetric; represents translations along the x and y axes and angular displacement about the and axes. The 4 × 1 vector F_{mag} denotes the electromagnetic force. B_{1} is the distribution matrix of magnetic force; is external force acting on the rotor (mainly are unbalanced forces). Details regarding the FE modeling procedures for rotors are available in [28].
Based on the developed FE model, the lumpedmass method is employed to obtain the modal frequencies of the rotor in the freefree state, which are 152.56 Hz, 423.59 Hz, 866.65 Hz, 1312.34 Hz, and 1752.92 Hz, respectively. However, the shrinkfitting silicon steel sheet and sleeves on the rotor will undoubtedly promote the bending stiffness, which has an important influence on the rotor modal frequency. Therefore, it is necessary to update the theoretical FE model of the rotor according to the experimental modal data of the rotor.
Experimental modal analysis is an effective method for updating rotor model. In order to get the real modal parameters of the rotor, a LMS vibration and noise testing system is employed to conduct modal test. The rotor is hung with flexible string to emulate the freefree state, and finally the first five bending modal frequencies of the rotor are obtained, as shown in Figure 8.
(a) Modal impact test
(b) Modal test results
Generally, the physical parameters of the rotor, such as material density, mass, and geometry, are constant. However, in the FE modeling process, the elastic modulus of elements becomes an engineering assumption due to influence of the shrinkfitting components. Thus, elastic modulus of the rotor with shrinkfitting components can be taken as a manual parameter, and the FE model can be made approximate to the real rotor characteristics by adjusting the parameters. After lots of trial and error, the updated model modal frequencies are approximate to the modal test data as listed in Table 1.

An important graphic tool for rotordynamic analysis is the undamped freefree mode shape plot, which is shown in Figure 9.
The undamped modal shape can reflect the behavior of the rotor at the sensors and actuators position which can make it clear to us whether the distribution of sensors and actuators is reasonable. Figure 9 shows that the first bending mode of the rotor has good controllability and observability from the control perspective.
Analyzing the rotor critical speed under different support stiffness is very important for predicting the resonance interval. Compared with ball bearing, the support stiffness of AMB is quite low. Thus, rotor critical speed is sensitive to support stiffness changes. This phenomenon is very important in the controller design because the controller should be robust to the rotor frequency perturbation. The rotor critical speed map is provided in Figure 10. It is indicated that there is a slight critical speed perturbation. This uncertainty must be considered in the controller design.
The order of updated rotor model obtained through modal test is very large, which will make controller design difficult. Therefore, a reduced order rotor model is needed. Generally, modal truncation is used to reduce the rotor model order. In this work, the first two bending modals of the rotor are preserved.
The gyroscopic effect can be ignored, because the axial/radial inertia ratio of the flexible rotor is always small. An axisymmetric rotor can be decoupled on the two vertical planes. If only the transverse vibration characteristics of a single plane ( plane) of the rotor are considered, its motion equations can be simplified as follows:
The equation can be transferred to modal space through the modal transfer matrix . Usually, the internal damping of the rotor is small. Thus, to better reflect the behavior of physical rotor, low levels of modal damping (0.3%) were added. The mass normalized motion equation can be obtained.where u = Φη; ζ = 0.003.
The rotor model above can be expressed as statespace form.The magnetic force can be assembled to the statespace equation of rotor. Finally, the orderreduced rotor model can be obtained.where , .
3.3. Modeling of Other Electrical Components
Power amplifier is a key component of the closedloop system which receives control voltage and generates enough current for AMB coils. The power amplifier dynamics is affected by the inductance of AMB coils and eddy current effects due to rotor motion. The dashed line in Figure 11 indicates the frequency response of power amplifier while the rotor is centered using copper foil. The transfer function fitted to the experimental frequency response is as follows:where K_{a} = 1, T_{a} = 0.00012, ω_{a} = 530 Hz, and ξ = 0.22. Since all amplifiers are alike, one transfer function is representative of each of their individual performance characteristics.
The transfer function can be converted into a statespace model as shown:
Four eddy current displacement sensors are employed to measure the rotor lateral vibration in the x and y direction. The sensors are equipped with a lowpass filter to remove the carrier signal. According to the parameters provided by the sensor manufacturer, the integrated transfer function of the displacement sensors and the antialiasing filter can be modeled aswhere = 5000, = 0.000056.
3.4. Model Validation
The AMB is openloop unstable, so system identification must be conducted in closedloop condition. Figure 12 shows the scheme of the sweeping measurement. When the rotor is suspended at the equilibrium position under PID controller, the sine sweeping excitation signal is superimposed into the input end of the power amplifier (PA). The control signal superimposed with sine sweep and sensor signal is sent to the frequency response analyzer. In this way, the frequency response characteristic of control plant can be obtained as shown in Figure 13.
4. Uncertainty Analysis Based on Gap Metric
RotorAMB system has many uncertainties, such as complicated geometry, nonlinear material properties in magnetic bearing laminations, shrink fits, assembling, and manufacture error. These uncertainties must be considered during controller design. However, taking all the uncertainties into consideration in the controller design will result in large conservatism, which will decrease the controller performance and increase the controller order. Thus, the type of uncertainty should be classified according to its effect on closedloop system robustness, determining which uncertainty can be easily suppressed by any feedback controller and which uncertainty would have serious impact on closedloop system stability. The uncertainty which has serious impact on closedloop system stability must be explicitly addressed to reduce conservatism in the controller design.
In this work, gap metric is employed for uncertainty analysis because it can evaluate the effect of openloop plant uncertainty from the closedloop perspective without assuming a specific feedback controller [29].
The vgap metric is defined in [30] aswhere and are the corresponding normalized right and left coprime factorization of the plant P_{1} and P_{2}; wno(g) denotes the winding number about the origin of g(s), as s follows the standard Nyquist Dcontour.
Nonconservative representation of uncertainty is the key to robust controller design. Based on the developed control plant model and physical system characteristics, some obvious uncertainties are extracted including loop gain, current stiffness , displacement stiffness , and rotor first bending mode frequency . Uncertainty of and can reach 15%~20%; especially have large changes due to the rotor displacement variation. Uncertainty of is usually set to 3%~5% in the robust controller design. The loop gain uncertainty can reach 10%. Table 2 shows the vgap metric value of each uncertainty when they have maximum variation.

It is indicated from Table 2 that the gap metric of loop gain, and , is small. These uncertainties can be well suppressed by any feedback controller. However, the gap metric value of rotor first bending frequency is very large even if it has a small variation, which must be addressed explicitly [29].
Model error and changing operating condition always result in modal frequency uncertainty. Gyroscopic effect and shrink fit property due to speed and temperature variation will induce modal frequency variation. Therefore, it is very important to design AMB feedback systems that are robust to modal frequency uncertainty.
5. μSynthesis Controller Design
For complex linear systems, synthesis is a multivariable controller design technique, which is highly suited for complex plants and MIMO control problems because it shows great superiority to handling flexible mode vibration. Synthesis takes the robust stability and performance as the optimized specifications explicitly and takes system uncertainties into account at the same time. These features are suitable for the flexible mode control of rotorAMB system. Although synthesis is currently still rarely used, it shows high technical potential for industrial applications with high robustness to plant and other uncertainties. A detailed discussion of the advantage of synthesis for the control design of rotorAMB system is given in [16].
We have made great modeling efforts for rotorAMB system above and a relatively accurate system model was obtained. Based on vgap, system uncertainties are analyzed. Special attention is paid to the design procedure of the modelbased synthesis controller in the following section. The block diagram of μsynthesis feedback loop design is shown in Figure 14, where e is the tracking error, u is the feedback control input, d represents the external disturbance, r is the reference signal, and z_{1} and z_{2} are weighted output. Δ is the parametric uncertainty, representing the rotor modal perturbation. G is the physical system and K is feedback controller which needs to be synthesized. W_{u}, W_{e}, W_{d}, and W_{r} are weighting functions to capture the performance specifications.
The system gain from disturbance inputs to performance outputs iswhere is the system sensitivity function and is the complementary sensitivity function. The synthesis procedure is to find a controller making structured singular value
5.1. Rotor Modal Frequency Uncertainty Representation
Appropriately accounting for model uncertainties is a crucial issue in the scope of robust control. The rotor modal frequency perturbation is parameter uncertainty with a very large gap metric value. System closedloop robust stability is very sensitive to its variation. Hence, it should be addressed explicitly. In this paper, the real uncertainty representation is employed. The uncertain modal state matrix for the bending mode in a single plane can be expressed:where is the bending mode frequency. is the damping ratio of the mode. When the uncertainty is small, the quadratic term of can be omitted and the above result holds. The singular value plot of the uncertain plant dynamics is shown in Figure 15. The plot close to the first and second bending mode frequency has dense lines indicating the uncertainties we characterize.
5.2. Weighting Function Selection
The weighting function associated with the synthesis controller design can be difficult to calculate accurately and few formal techniques exist for the proper selection. Therefore, the selection of weighting functions can be directed based on the prior experience [10].
Generally, the weight is used to bound the system sensitivity peak. Thus, is usually chosen as a constantwhere represents the sensitivity peak according to ISO148393.
Weight is usually taken as a tune button to bound GS and is selected between 1/5 and 1. In this paper, weight is as
is the performance specification weight, which can strongly influence the steady error, overshoot and settling time of the rotorAMB system. A constant weighting on sensitivity often results in a controller similar to a PD controller. In order to balance the rotor sag, a small integration is needed. Therefore, where is the desired bandwidth and A = 0.0733 is the integration constant; gain a = 2.6 is used as additional parameter to help and bound the S and GS.
In general, which bounds the control sensitivity function is often selected as second order highpass filter, enforcing a sharper controller gain rolloff at high frequency. After much trial and error, the final is chosen as where 6 represents the rolloff frequency and add high frequency poles so that the inverse of the transfer function is proper.
Bode magnitude plots of the above weighting functions are shown in Figures 16 and 17. In general, these weighting functions are proper transfer functions.
The weighting functions defined above are SISO transfer function or scalar constants. However, in the controller synthesized procedure, the weights should be MIMO transfer function matrices. The performance requirement and dynamics of the four channels are similar. Therefore, the same weights are assigned to each channel.
5.3. Synthesis Results and Evaluation of Controller
The Matlab function dksyn is used to perform the DK iteration steps of the μ controller synthesis. After several iterations, a 30th order controller is obtained. However, such a high order controller needs heavy calculation burden and is difficult to implement with hardware limitations. Hence, an orderreduced controller with 14 states is retained based on balanced stochastic model truncation method. The Hankel singular values of the fullorder controller are shown in Figure 18.
The singular value of the controller is shown in Figure 19. In the low frequency range, the controller has integral effect which is enforced by weighting function . This can guarantee the rotor better steady error and static stiffness to balance rotor weight. In the middle frequency range, the controller is similar to the PD controller with stiffness and damping to stabilize the rigid body modes. In the high frequency, the controller places a series of poles and zeros to stabilize the flexible rotor bending modes and enforce sharp gain rolloff.
Sensitivity function is an important mean to evaluate the stability margin of the system [31]. Although the sensitivity peak is high in the rigid body mode, the low frequency uncertainty does not show a serious problem of robustness. On the contrary, the sensitivity in the flexible mode becomes more important. As can be seen from the vgap value in the previous section, a slight perturbation at the flexible mode leads to serious closedloop robust stability problems. Hence, in order to verify the robustness of μsynthesis controller to rotor flexible mode perturbation, the output sensitivity functions of the DEDE channel with ±5% natural frequency perturbation are shown in Figures 20 and 21. It is indicated that system sensitivity at the flexible mode still has low bound with ±5% natural frequency perturbation. The designed controller has high robustness to modal frequency uncertainty.
The system compliance transfer function is shown in Figure 22. In contrast to the sensitivity function measurement with the sole purpose of system robustness assessment, the compliance transfer function measurement constitutes a means for assessing the controlled system’s performance, most notably in terms of damping of resonance [27]. From the compliance transfer function plot, the damping ratio of rotor first bending mode can be obtained with 17.35%, which demonstrates the μsynthesis controller enforces much active damping for the flexible mode.
6. Experimental Validation
The controller designed in Section 5 is implemented on the flexible rotorAMB test rig. The experimental setup is shown in Figure 23, including dSPACE, monitor computer, amplifier, and flexible rotor test rig. The amplifier is a current feedback PWM switching amplifier. The PWM switching frequency is 26 kHz and the amplifier bandwidth is 870 Hz. The eddy current displacement sensor we used is MICROEPSILON made in Germany and its static resolution reaches 0.05 μm. A lowpass filter is connected in series after the displacement sensor to remove the high frequency noise. The μsynthesis controller is discretized at 20 kHz and executed via the dSPACE DS1103 rapid control prototype. Detailed data of the experiment platform are listed in Table 3.

First, the initial levitation test is conducted with μsynthesis controller. Figure 24 shows the rotor positon from the rest state to the equilibrium position under μsynthesis controller within 0.1 s, which indicates good dynamic performance.
Sensitivity is an important specification of the robust stability margin. The sensitivity function can be measured directly by using online frequency sweep test [27]. Figure 25 shows the output sensitivity of the system under μsynthesis controller at the DE and NDE. It is indicated that the test results show good agreement with model prediction. The measured sensitivity peak is 9.405 dB located in Zone A [31]. The sensitivity nearby the rotor first bending mode under μsynthesis controller is much lower, which indicates that the μsynthesis controller has better robustness to the modal frequency uncertainty.
The compliance transfer functions at standstill are obtained by a swept sine measurement under μsynthesis controller. The nominal model was compared with the test data as shown in Figure 26. Based on the compliance measurement data, the damping ratios at the flexible mode can be estimated, which are 10.93% and 12.89%, respectively, at the DE and NDE. The μ controller enforced much active damping on the first flexible mode.
The rotor runup test is conducted. Figure 27 shows the rotor orbit changes with speed under μ controller. The flexible rotor successfully passes the first bending critical speed and reaches the rated speed 14400 rpm. The rotor orbit at the NDE shows obvious rigid mode critical speed peak. However, the orbit at the DE does not show obvious rigid mode critical speed peak due to the influence of flexible coupling. The rotor vibration level of steadystate operation at 14400 rpm is located in Zone A [32].
7. Conclusions
In this paper, the resonance vibration control for AMB flexible rotor based on μsynthesis controller is investigated. A flexible rotor test rig is designed and the system model is obtained through the combination of system identification and theoretical modeling method. With the help of gap metric tool, rotor modal frequency perturbation has significant effect on the system stability, which must be addressed explicitly. The vgap value reaches 0.9835 under ±5% modal perturbation. Then, the synthesis controller is designed based on the uncertain system model. The simulation and experimental results show that the μ controller is robust to the modal perturbation and enforces much active damping on the first flexible mode (damping ratio reaches 12.89%). In the runup test, the rotor amplitude is within 45 μm nearby the first bending critical speed, which indicates that the designed controller has good resonance vibration control performance for AMB flexible rotor system.
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This research is financially supported by the National Natural Science Foundation of China (No. 51575411).
References
 H. Yefa, Z. Zude, and J. Zhengfeng, Basic Theory and Application of Magnetic Bearings, China Machine Press, 2006.
 B. Han, Q. Xu, and Q. Yuan, “Multiobjective optimization of a combined radialaxial magnetic bearing for magnetically suspended compressor,” IEEE Transactions on Industrial Electronics, vol. 63, no. 4, pp. 2284–2293, 2016. View at: Google Scholar
 Y. Le and K. Wang, “Design and optimization method of magnetic bearing for highspeed motor considering eddy current effects,” IEEE/ASME Transactions on Mechatronics, vol. 21, no. 4, pp. 2061–2072, 2016. View at: Publisher Site  Google Scholar
 C. Wang and J. Tang, “Design and mathematical analysis of a novel reluctance forcetype hybrid magnetic bearing for flywheel with gimballing capability,” Mathematical Problems in Engineering, vol. 2013, Article ID 836058, 17 pages, 2013. View at: Publisher Site  Google Scholar
 Q. Yuan, B. Han, and B. Dong, “The multiple objective optimization of highspeed rotor supported by magnetic bearing in BLDCM,” International Journal of Applied Electromagnetics and Mechanics, vol. 46, no. 3, pp. 663–673, 2014. View at: Google Scholar
 M. A. Rahman, A. Chiba, and T. Fukao, “Super high speed electrical machines  Summary,” in Proceedings of the 2004 IEEE Power Engineering Society General Meeting, pp. 1272–1275, USA, June 2004. View at: Google Scholar
 S. Islam, P. X. Liu, and A. El Saddik, “Robust control of fourrotor unmanned aerial vehicle with disturbance uncertainty,” IEEE Transactions on Industrial Electronics, vol. 62, no. 3, pp. 1563–1571, 2015. View at: Publisher Site  Google Scholar
 C. Zhang, G. Zhao, and K. J. Tseng, “Robust decoupling control of a threelink robot arm directly driven by permanent magnet synchronous motors,” International Journal of Applied Electromagnetics and Mechanics, vol. 36, no. 4, pp. 279–293, 2011. View at: Google Scholar
 R. Conway, J. Choi, R. Nagamune, and R. Horowitz, “Robust trackfollowing controller design in hard disk drives based on parameter dependent lyapunov functions,” IEEE Transactions on Magnetics, vol. 46, no. 4, pp. 1060–1068, 2010. View at: Publisher Site  Google Scholar
 L. Guoxin, Robust stabilization of rotoractive magnetic bearing systems [Ph.D. thesis], University of Virginia, 2007.
 A. C. Wroblewski, Model identification, updating, and validation of an active magnetic bearing highspeed machining spindle for precision machining operation [Ph.D. Dissertation], Cleveland State University, 2011.
 T. Enqiong, F. Jiancheng, and H. Bangcheng, “Active vibration control of the flexible rotor in high energy density magnetically suspended motor with mode separation method,” Journal of Engineering for Gas Turbines and Power, vol. 137, Article ID 082503, 2015. View at: Google Scholar
 S. Zheng, H. Li, C. Peng, and Y. Wang, “Experimental investigations of resonance vibration control for noncollocated AMB flexible rotor systems,” IEEE Transactions on Industrial Electronics, vol. 64, no. 3, pp. 2226–2235, 2017. View at: Publisher Site  Google Scholar
 E. Tang, B. Han, and Y. Zhang, “Optimum compensator design for the flexible rotor in magnetically suspended motor to pass the first bending critical speed,” IEEE Transactions on Industrial Electronics, vol. 63, no. 1, pp. 343–354, 2016. View at: Publisher Site  Google Scholar
 Z. Gosiewski and A. Mystkowski, “The robust control of magnetic bearings for rotating machinery,” Solid State Phenomena, vol. 113, pp. 125–130, 2006. View at: Publisher Site  Google Scholar
 E. H. Maslen and J. T. Sawicki, “Musynthesis for magnetic bearings: Why use such a complicated tool,” ASME International Mechanical Engineering Congress and Exposition, pp. 1103–1112, 2007. View at: Google Scholar
 A. Lanzon and P. Tsiotras, “Robust control of energy momentum wheels supported on active magnetic bearings using H∞ loopshaping and μsynthesis,” in Proceedings of the 15th World Congress of the International Federation of Automatic Control, 2002, pp. 121–126, Spain, July 2002. View at: Google Scholar
 H. Seto and T. Namerikawa, “An H∞ control system design of a magnetic bearing considering periodic disturbance,” in Proceedings of the IEEE Conference on Control Applications, pp. 227–232, 2005. View at: Google Scholar
 Z. Gosiewski and A. Mystkowski, “Robust control of active magnetic suspension: analytical and experimental results,” Mechanical Systems and Signal Processing, vol. 22, no. 6, pp. 1297–1303, 2008. View at: Publisher Site  Google Scholar
 N. S. Gibson, Hinfinity Control of Active Magnetic Bearings: An Intelligent Uncertainty Modeling Approach [Ph.D. Dissertation], North Carolina State University, 2004.
 S. J. M. Steyn, Multivariable H∞ control for an active magnetic bearing flywheel system [M.S. Dissertation], North West University, 2010.
 A. Noshadi, J. Shi, W. S. Lee, P. Shi, and A. Kalam, “System identification and robust control of multiinput multioutput active magnetic bearing systems,” IEEE Transactions on Control Systems Technology, vol. 24, no. 4, pp. 1227–1239, 2016. View at: Publisher Site  Google Scholar
 R. P. Jastrzebski, K. M. Hynynen, and A. Smirnov, “H∞ control of active magnetic suspension,” Mechanical Systems and Signal Processing, vol. 24, no. 4, pp. 995–1006, 2010. View at: Publisher Site  Google Scholar
 Y.W. Long, Z.Y. Xie, and X. Xu, “H∞ robust control strategy for an active magnetic bearing,” Journal of Vibration and Shock, vol. 32, no. 23, pp. 115–120, 2013. View at: Google Scholar
 J. Zhou, S. Zheng, and B. Han, “Application of μ synthesis to the control of active magnetic bearings with parametric uncertainty,” in Proceedings of the 20th International Conference on Electrical Machines and Systems (ICEMS), pp. 1–5, August 2017. View at: Google Scholar
 S. E. Mushi, Z. Lin, and P. E. Allaire, “Stability analysis for a flexible rotor on active magnetic bearings subject to aerodynamic loads,” in Proceedings of the 12th international symposium on magnetic bearings (ISMB '12), pp. 22–25, 2010. View at: Google Scholar
 G. Schweitzer and E. H. Maslen, Magnetic Bearing: Theory, Design, and Application to Rotating Machinery, Springer, Berlin, Germany, 2009.
 M. I. Friswell, Dynamics of Rotating Machines, Cambridge University Press, Cambridge, UK, 2010.
 G. Li, Z. Lin, and P. Allaire, “Uncertainty classification for rotorAMB systems,” in Proceedings of the International Symposium on Magnetic Bearings, 2006. View at: Google Scholar
 G. Vinnicombe, “Uncertainty and feedback—H∞ loopshaping and the vGap metric,” Uncertainty and feedback, pp. 117–119, 2001. View at: Google Scholar
 ISO 148393, Mechanical vibration  Vibration of rotating machinery equipped with active magnetic bearings  Part 3: Evaluation of stability margin, International Organization for Standardization ISO, 2006.
 ISO 148392, Mechanical vibrationVibration of rotating machinery equipped with active magnetic bearings  Part 2: Evaluation of vibration, International Organization for Standardization ISO, 2006.
Copyright
Copyright © 2018 Shaolin Ran et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.