Resonance Vibration Control for AMB Flexible Rotor System Based on μ-Synthesis Controller

The resonance vibration control of flexible rotor supportedon activemagnetic bearings (AMB) is a challenging issue in the industrial applications.This work addresses the application of robust controlmethod to the resonance vibration control for AMBflexible rotor while passing through the critical speed. This model-based 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.


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, no-friction, and no-lubrication [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 energy-saving and emission-reduction [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 open-loop unstable system, which has fundamental limitations for closed-loop 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][13][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][18][19] and mainly focuses on the rigid rotor [20][21][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 cross-coupled 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 cross-coupled 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 rotor-AMB 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.

The Flexible Rotor-AMB Test Rig
The purpose of the flexible rotor-AMB 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.

Modeling of AMB System
-Synthesis is fundamentally a model-based 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 model-based 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 silicon-iron.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.
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 force-current 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 current-displacement curve, which can also be obtained linear fitting.Finally,   = 213.5 N/A,   = -1110.2N/mm.

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 skew-symmetric;  represents translations along the x and y axes and angular displacement about the  and  axes.The 4 × 1 vector F  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 lumped-mass method is employed to obtain the modal frequencies of the rotor in the free-free state, which are 152.56Hz, 423.59 Hz, 866.65 Hz, 1312.34Hz, and 1752.92Hz, respectively.However, the shrink-fitting 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 free-free state, and finally the first five bending modal frequencies of the rotor are obtained, as shown in Figure 8.
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 shrink-fitting components.Thus, elastic modulus of the rotor with shrink-fitting 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 free-free 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.
The rotor model above can be expressed as state-space form.
The magnetic force can be assembled to the state-space equation of rotor.Finally, the order-reduced rotor model can be obtained. where

Modeling of Other Electrical Components.
Power amplifier is a key component of the closed-loop 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  = 1, T  = 0.00012,   = 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  a () 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 low-pass 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  s () can be modeled as where  s = 5000,  s = 0.000056.

Model Validation.
The AMB is open-loop unstable, so system identification must be conducted in closed-loop 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.

Uncertainty Analysis Based on V-Gap Metric
Rotor-AMB 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 closed-loop system robustness, determining which uncertainty can be easily suppressed by any feedback controller and Frequency (Hz) Frequency (Hz) which uncertainty would have serious impact on closed-loop system stability.The uncertainty which has serious impact on closed-loop system stability must be explicitly addressed to reduce conservatism in the controller design.
In this work, V-gap metric is employed for uncertainty analysis because it can evaluate the effect of open-loop plant uncertainty from the closed-loop perspective without assuming a specific feedback controller [29].
The v-gap metric is defined in [30] as where   and G ( = 1, 2) 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 D-contour.
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  1 .Uncertainty of   and  x can reach 15%∼20%; especially   have large changes due to the rotor displacement variation.Uncertainty of  1 is usually set to 3%∼5% in the robust controller design.The loop gain uncertainty can reach 10%.Table 2 shows the v-gap metric value of each uncertainty when they have maximum variation.
It is indicated from Table 2 that the V-gap metric of loop gain,   and   , is small.These uncertainties can be well suppressed by any feedback controller.However, the V-gap metric value of rotor first bending frequency is very large even Frequency (Hz)

Mathematical Problems in Engineering
Frequency (Hz) Frequency (Hz) Frequency (Hz) Frequency (Hz) Frequency (Hz)   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.

𝜇-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 rotor-AMB 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 rotor-AMB system is given in [16].
We have made great modeling efforts for rotor-AMB system above and a relatively accurate system model was obtained.Based on v-gap, 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  , W  , W  , and W  are weighting functions to capture the performance specifications.
The system gain from disturbance inputs to performance outputs is where  is the system sensitivity function and  is the complementary sensitivity function.The synthesis procedure is to find a controller  making structured singular value

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 V-gap metric value.System closed-loop 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  th 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.

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 constant where  S = 3 represents the sensitivity peak according to ISO14839-3.
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 rotor-AMB 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  0 = 30 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 high-pass filter, enforcing a sharper controller gain roll-off at high frequency.After much trial and error, the final   is chosen as where  2 = 486 represents the roll-off frequency and  1 = 8 2 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.

Synthesis Results and Evaluation of Controller.
The Matlab function dksyn is used to perform the D-K 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 balanced stochastic model truncation method.The Hankel singular values of the full-order 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 roll-off.
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 v-gap value in the previous section, a slight perturbation at the flexible mode leads to serious closed-loop robust stability problems.Hence, in order to verify the robustness of -synthesis controller to rotor flexible mode perturbation, the output sensitivity functions of the DE-DE 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.

Experimental Validation
The controller designed in Section 5 is implemented on the flexible rotor-AMB 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 MICRO-EPSILON made in Germany and its static resolution reaches 0.05 m.A low-pass 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  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 run-up 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 steady-state operation at 14400 rpm is located in Zone A [32].

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 V-gap metric tool, rotor modal frequency perturbation has significant effect on the system stability, which must be addressed explicitly.The v-gap 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 run-up 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.

Figure 3 :
Figure 3: Typical control schematic diagram of magnetic bearing system.

Figure 5 :
Figure 5: Relationship of forces and currents.

Figure 6 :
Figure 6: Relationship of currents and displacements.

Figure 7 :
Figure 7: The FE model of the rotor.

Figure 10 :
Figure 10: Relationship between critical speed and support stiffness.

Figure 15 :
Figure 15: Singular value plot of uncertain system.

Figure 24 :
Figure 24: Levitation response under  controller at the DE.

Table 2 :
V-gap metric value of uncertainties.

Table 3 :
Parameter of experiment platform.