The control systems applied on active magnetic bearing are several. A perfect levitation is characterized by maintaining the operating point condition that is characterized by the center of stator coincident with the geometric center of shaft. The first controller implemented for this purpose is PID controller that is characterized by an algorithm that leads the amplifier to produce control current until the operating point condition is not reached, this is obtained by an integration operator. The effect of an integrator is essential but not necessary for a centered levitation for example in the robust control characterized by a dynamic model depended on plant of system so that it depends on angular speed as LQR controller does. In LQR there is not integrator so there is not a perfectly centered section of shaft with center of stator. On contrary PID controller does not depend on angular speed and it can be easily implemented according some simple rules. Predictive control is another interesting controller characterized by a multiple controller operating in different condition in order to get the minimum of cost function, but also in this case the angular speed is introduce for the same reason discussed before.

1. Introduction

Active Magnetic Bearings (AMBs) use electromagnets to attract the ferromagnetic cape winding the rotor which is free to rotate with no physical contact with the bearing. This operation, called active magnetic levitation, is unstable unless of a certain control’s algorithm performed respecting the imposed constraining, [1]. In order to achieve a stable levitation, an active feedback control of the current in the magnetic coils is necessary. As it might be expected, a variety of control schemes are used and a variety of studies have been done for AMBs control. The dynamic system, however, depends above all on the rotor’s angular speed, as the case for any rotor dynamic system because of the gyroscopic effect [2]. The gyroscopic effect leads the system into instability phenomena which must be considered to achieve stable levitation. The rotor motion is characterized by translation along 𝑥-𝑦 directions and rotation of rotor around those axes [1, 3, 4]. Moreover, the transformation coordinates allow relating the displacements captured by the sensors and the displacements of the section located on the middle plane of the bearings location. Without a suitable control system, no magnetic levitation is possible [1]. During last twenty years, a fast number of control systems were applied on magnetic levitation in order to provide enough acknowledge about the capability of rotor to maintain the contactless between the rotor and stator. Obviously while there are many control systems which are not able to maintain the operating point position of the controlled sections without a further algorithm such as the integrator of a PID controller [5], some other control systems need the entire state vector to create the feedback such as the optimal control [3] characterized by a matrix whose number of column is equal to the dimension of state vector. The cutting edge of control systems is represented by 𝜇-synthesis, loop-shaping design procedure, and H-infinity robust control, the latter with its variant sub(𝐻)∞. The reason is not only the recovering of operating point position without integrator but also the possibility to delete the use of some sensors to capture further components belonging to the state vector, a problem that is usually solved by the introduction of an observer. There is something more convenient than the previous simplification which is more related to the performances of the machine in using robust control. The advantages derived by using robust control are its possibility in case of presence of dynamic perturbation such as, no modeled dynamic, neglected nonlinearities, effects of reduced-order models, system-parameter variation due to environmental changes, hysteresis, and torn and worn factors. Moreover, it is used also in case there is a presence of sensor and actuator noise. An application of 𝜇-synthesis controller was for the Army's weapon system. The test fixture is patterned after the Apache helicopter's 30 mm gun and has tunable nonlinearities which may be representative not only of the nonlinearities of the gun, but also to other mechanical systems as well. The models of the test fixture which were available at the time of the work are also described. The goal of pointing the gun is to reduce dispersions of fired gun rounds on targets. The resulting 𝜇-synthesis design, when connected with a nonlinear simulation, exhibited limit-cycle behavior of unacceptable amplitude [6]. Due to high surface speed and active control capabilities, active magnetic bearings hold great promise for high speed machining spindles. The control problem posed by this application is examined, and the development of an advanced prototype is reviewed. A 𝜇-synthesis framework is proposed for this problem, and it is shown that the minimization of the susceptibility to machining chatter may be easily put into this framework. In addition to handling uncertainties in sensor and actuator components, this formulation may also include an uncertainty representing the range of cutting tools for the spindle, [7].

The proposed control algorithms are developed using 𝜇-analysis to obtain robust stability and robust performance in simulation investigation. In simulations work, three different active vibration control algorithms are used. A similar approach was applied in [8] where a comparison between three different controllers was performed in order to analyze the differences on the dynamic behavior. Many other applications of robust control which were performed through loop-shaping design procedure such as in [9] where an ğ»âˆž controller was performed by evolution optimization to control a robot arm.

The loop-shaping method is commonly used also to obtain tradeoffs of robust stability and robust performance. This technique is a particular optimization problem to guarantee closed loop stability at all frequencies [10].

2. Mathematical Model

The particular configuration shown in this work considers a rotor with four degrees of freedom with eight poles for each active magnetic bearing, having a slope of 45° with regard to horizontal direction so that the force’s resultant supports the rotor along the 𝑥 and 𝑦 directions crossing the center of mass and rotation around it as Figure 1 shows.

The system is subjected to a state of uncertainty about its mass, cross, and polar moment of inertia dictated by the parameters 𝛿𝑚,𝛿𝐼𝑃,and𝛿𝐼𝑇 in the ranges, respectively, 𝑃𝑚,𝑃𝐼𝑇,andP𝐼𝑃. The equation of motion is referred to the center of gravity, and it has the following expression:𝐌𝑑2𝐪𝑔(𝑡)𝑑𝑡2+Ω⋅𝐆𝑑𝐪𝑔(𝑡)𝑑𝑡=𝐁Θmag𝐁mag𝐟𝐢𝑐(𝑡),𝐪𝑏,(𝑡)(1) where âŽ¡âŽ¢âŽ¢âŽ¢âŽ¢âŽ¢âŽ£ğŒ=𝑚0000𝐼𝑇0000𝑚0000ğ¼ğ‘‡âŽ¤âŽ¥âŽ¥âŽ¥âŽ¥âŽ¥âŽ¦âŽ¡âŽ¢âŽ¢âŽ¢âŽ¢âŽ¢âŽ£,𝐆=0000000−𝐼𝑃00000ğ¼ğ‘ƒâŽ¤âŽ¥âŽ¥âŽ¥âŽ¥âŽ¥âŽ¦,𝐁00Θmag=⎡⎢⎢⎢⎢⎢⎣⎤⎥⎥⎥⎥⎥⎦,𝐁cos(𝛼)0−sin(𝛼)00cos(𝛼)0sin(𝛼)sin(𝛼)0cos(𝛼)00−sin(𝛼)0cos(𝛼)mag=⎡⎢⎢⎢⎢⎢⎣1010−𝑙𝑏𝐴0𝑙𝑏𝐵001010𝑙𝑏𝐴0âˆ’ğ‘™ğ‘ğµâŽ¤âŽ¥âŽ¥âŽ¥âŽ¥âŽ¥âŽ¦,𝐟𝐢𝑐(𝑡),𝐪𝑏(=âŽ¡âŽ¢âŽ¢âŽ¢âŽ¢âŽ¢âŽ£ğ¹ğ‘¡)𝑥𝐴𝑖𝑐𝑥𝐴(𝑡),𝑥𝑏𝐴𝐹(𝑡)𝑦𝐴𝑖𝑐𝑦𝐴(𝑡),𝑦𝑏𝐴𝐹(𝑡)𝑥𝐵𝑖𝑐𝑥𝐵(𝑡),𝑥𝑏𝐵𝐹(𝑡)𝑦𝐵𝑖𝑐𝑦𝐵(𝑡),ğ‘¦ğ‘ğµî€¸âŽ¤âŽ¥âŽ¥âŽ¥âŽ¥âŽ¥âŽ¦,𝐪(𝑡)𝑔𝑥(𝑡)=𝑔(𝑡)𝛾(𝑡)𝑦𝑔(𝑡)𝜗(𝑡)𝐓,𝑚=𝑚1+𝑃𝑚𝛿𝑚=𝑚+𝑚⋅𝑃𝑚𝛿𝑚=𝑚+Δ𝐼𝑚,𝑇=𝐼𝑇1+𝑃𝐼𝑇𝛿𝐼𝑇=𝐼𝑇+𝐼𝑇⋅𝑃𝐼𝑇𝛿𝐼𝑇𝑇=𝐼𝑇+Δ𝐼𝑇,𝐼𝑃=𝐼𝑃1+𝑃𝐼𝑃𝛿𝐼𝑃=𝐼𝑃+𝐼𝑃⋅𝑃𝐼𝑃𝛿𝐼𝑃=𝐼𝑃+Δ𝐼𝑃.(2)

The last expression leads to structured uncertainties matrix such as𝐌=𝐌+𝐌𝐏𝐌Δ𝐌∈ℜ𝟒×𝟒,𝐆=𝐆+𝐆𝐏𝐆Δ𝐆∈ℜ𝟒×𝟒,(3) whereâŽ¡âŽ¢âŽ¢âŽ¢âŽ¢âŽ¢âŽ£ğŒ=0𝑚000𝐼𝑇0000𝑚0000ğ¼ğ‘‡âŽ¤âŽ¥âŽ¥âŽ¥âŽ¥âŽ¥âŽ¦,âŽ¡âŽ¢âŽ¢âŽ¢âŽ¢âŽ¢âŽ£ğ†=0000000−𝐼𝑃00000ğ¼ğ‘ƒâŽ¤âŽ¥âŽ¥âŽ¥âŽ¥âŽ¥âŽ¦,00𝐌𝐏𝑀Δ𝑀=ğŒâŽ¡âŽ¢âŽ¢âŽ¢âŽ¢âŽ¢âŽ£ğ‘ƒğ‘š0000𝑃𝐼𝑇0000𝑃𝑚0000ğ‘ƒğ¼ğ‘‡âŽ¤âŽ¥âŽ¥âŽ¥âŽ¥âŽ¥âŽ¦âŽ¡âŽ¢âŽ¢âŽ¢âŽ¢âŽ¢âŽ£ğ›¿ğ‘š0000𝛿𝐼𝑇0000𝛿𝑚0000ğ›¿ğ¼ğ‘‡âŽ¤âŽ¥âŽ¥âŽ¥âŽ¥âŽ¥âŽ¦,𝐆𝐏𝐺Δ𝐺=ğ†âŽ¡âŽ¢âŽ¢âŽ¢âŽ¢âŽ¢âŽ£00000𝑃𝐼𝑃000000000ğ‘ƒğ¼ğ‘ƒâŽ¤âŽ¥âŽ¥âŽ¥âŽ¥âŽ¥âŽ¦âŽ¡âŽ¢âŽ¢âŽ¢âŽ¢âŽ¢âŽ£00000𝛿𝐼𝑃000000000ğ›¿ğ¼ğ‘ƒâŽ¤âŽ¥âŽ¥âŽ¥âŽ¥âŽ¥âŽ¦,(4) respectively, the mass and gyroscopic nominal matrix and structured uncertainties matrices relative to mass and gyroscopic effect.

By introducing a transformation of coordinates (5), which has a great consideration in the relation of the displacements between the section relative to bearing location and the sensors, the system is analyzed according the equation of motion (6): 𝐪𝑏(𝑡)=𝐁Θdisp𝐁𝐓mag𝐪𝑔𝐪(𝑡),sensor(𝐁𝑡)=Θdisp𝐁sensor𝐪𝑔(𝐪𝑡),sensor𝐁(𝑡)=Θdisp𝐁sensor𝐁Θdisp𝐁𝐓mag−1𝐪𝑏𝐟(𝑡),mag_to_𝑔(𝑡)=𝐁Θmag𝐁mag𝐟𝐢𝑐(𝑡),𝐪𝑏,(𝑡)(5)𝐌𝑑2𝐪𝑔𝑑𝑡2+Ω⋅𝐆𝑑𝐪𝑔𝑑𝑡=𝐁Θmag𝐁mag𝐟𝐮(𝑡),𝐪𝑏,(𝑡)(6) where 𝐁Θdisp=⎡⎢⎢⎢⎢⎢⎣⎤⎥⎥⎥⎥⎥⎦,𝐁cos(𝛼)sin(𝛼)00−sin(𝛼)cos(𝛼)0000cos(𝛼)sin(𝛼)00−sin(𝛼)cos(𝛼)sensor=⎡⎢⎢⎢⎢⎢⎣1−𝑙𝑠𝐴00001𝑙𝑠𝐴1𝑙𝑠𝐵00001âˆ’ğ‘™ğ‘ ğµâŽ¤âŽ¥âŽ¥âŽ¥âŽ¥âŽ¥âŽ¦,(7)𝐪𝑏𝑥(𝑡)=𝑏𝐴(𝑡)𝑦𝑏𝐴(𝑡)𝑥𝑏𝐵(𝑡)𝑦𝑏𝐵(𝑡)𝐓,𝐪(𝑡)sensor=𝑥𝑠𝐴(𝑡)𝑦𝑠𝐴(𝑡)𝑥𝑠𝐵(𝑡)𝑦𝑠𝐵(𝑡)𝐓.(8) The magnetic force (9a)–(9d) produced by an electromagnet of active magnetic bearing is linearized by Taylor series expansion which leads to the expression of the force (9a), [5]: 𝑓𝑖𝑐𝑐(𝑖𝑥,𝑖𝑦)(𝐴,𝐵)(𝑡),ğ‘žğ‘(𝑥,𝑦)(𝐴,𝐵)𝑖(𝑡)=𝑘0(𝐴,𝐵)∓𝑖𝑐(𝑖𝑥,𝑖𝑦)(𝐴,𝐵)𝑔0(𝐴,𝐵)Â±ğ‘žğ‘(𝑥,𝑦)(𝐴,𝐵)2,𝐟𝐢𝑐(𝑡),𝐪𝑏(𝑡)≈𝐊𝑆𝐪𝑏(𝑡)+𝐊𝐼𝐢𝑐(𝑡),(9a) where: 𝐊𝑆=âŽ¡âŽ¢âŽ¢âŽ¢âŽ¢âŽ¢âŽ£ğ‘˜ğ‘¥ğ´0000𝑘𝑦𝐴0000𝑘𝑥𝐵0000ğ‘˜ğ‘¦ğµâŽ¤âŽ¥âŽ¥âŽ¥âŽ¥âŽ¥âŽ¦,𝐊𝐼=âŽ¡âŽ¢âŽ¢âŽ¢âŽ¢âŽ¢âŽ£ğ‘˜ğ‘–ğ‘¥ğ´0000𝑘𝑖𝑦𝐴0000𝑘𝑖𝑥𝐵0000ğ‘˜ğ‘–ğ‘¦ğµâŽ¤âŽ¥âŽ¥âŽ¥âŽ¥âŽ¥âŽ¦,𝐢𝑐(âŽ¡âŽ¢âŽ¢âŽ¢âŽ¢âŽ¢âŽ£ğ‘–ğ‘¡)=𝑐𝑥𝐴𝑖(𝑡)𝑐𝑦𝐴𝑖(𝑡)𝑐𝑥𝐵(𝑖𝑡)ğ‘ğ‘¦ğµâŽ¤âŽ¥âŽ¥âŽ¥âŽ¥âŽ¥âŽ¦,(𝑡)(9b) where the elements of 𝐊𝑆 and 𝐊𝐼, respectively, displacement and current gains matrices are 𝑘(𝑥,𝑦)(𝐴,𝐵)=𝑖𝜕𝑓𝑐𝑐(𝑖𝑥,𝑖𝑦)(𝐴,𝐵)(𝑡),ğ‘žğ‘(𝑥,𝑦)(𝐴,𝐵)(𝑡)ğœ•ğ‘žğ‘(𝑥,𝑦)(𝐴,𝐵)||||𝑖𝑐(𝑖𝑥,𝑖𝑦)(𝐴,𝐵)ğ‘ž=0𝑏(𝑥,𝑦)(𝐴,𝐵)=0𝑖=2𝑘0(𝐴,𝐵)−𝑖𝑐(𝑖𝑥,𝑖𝑦)(𝐴,𝐵)2𝑔0(𝐴,𝐵)+ğ‘žğ‘(𝑥,𝑦)(𝐴,𝐵)𝑔0(𝐴,𝐵)+ğ‘žğ‘(𝑥,𝑦)(𝐴,𝐵)4|||||𝑖𝑐(𝑖𝑥,𝑖𝑦)(𝐴,𝐵)ğ‘ž=0𝑏(𝑥,𝑦)(𝐴,𝐵)=0=2𝑘(𝐴,𝐵)𝑖20(𝐴,𝐵)𝑔30(𝐴,𝐵),𝑘(𝑖𝑥,𝑖𝑦)(𝐴,𝐵)=𝑖𝜕𝑓𝑐(𝑖𝑥,𝑖𝑦)(𝐴,𝐵)(𝑡),ğ‘žğ‘(𝑥,𝑦)(𝐴,𝐵)(𝑡)𝜕𝑖𝑐(𝑖𝑥,𝑖𝑦)(𝐴,𝐵)||||𝑖𝑐(𝑖𝑥,𝑖𝑦)(𝐴,𝐵)ğ‘ž=0𝑏(𝑥,𝑦)(𝐴,𝐵)=0𝑖=2𝑘0(𝐴,𝐵)−𝑖𝑐(𝑖𝑥,𝑖𝑦)(𝐴,𝐵)𝑔0(𝐴,𝐵)+ğ‘žğ‘(𝑥,𝑦)(𝐴,𝐵)2|||||𝑖𝑐(𝑖𝑥,𝑖𝑦)(𝐴,𝐵)ğ‘ž=0𝑏(𝑥,𝑦)(𝐴,𝐵)=0=2𝑘(𝐴,𝐵)𝑖0(𝐴,𝐵)𝑔20(𝐴,𝐵),(9c) where 𝑘(𝐴,𝐵)=14𝜇0𝑁2(𝐴,𝐵)𝐴𝑔(𝐴,𝐵)cos𝛽(9d)is the multiplicative constant carried out from principle of virtual work and 𝛽 is the angle of two poles of one electromagnet.

In this paper, the value of displacement and current gains are taken from a real model of radial active magnetic bearing produced by SKF. The physical meaning on these gain matrices is related to the reaction force produced for the unit of displacement and current, respectively, for displacement and current gains

Due to the transformation coordinates (5), we transform (9a)–(9d) in𝐟𝐢𝑐(𝑡),𝐪𝑏(𝑡)≈𝐊𝑆𝐁Θdisp𝐁𝐓mag𝐪𝑔(𝑡)+𝐊𝐼𝐢𝑐(𝑡).(10) But we do not have to forget that mathematical model in a state space form considers the matrix product 𝐁Θmag𝐁mag before 𝐟(𝐢𝑐(𝑡),𝐪𝑏(𝑡)) as in (6).

3. Controller

In order to provide a stabilizing effect to control the position of the rotor, a suitable control system must be performed because no magnetic levitation can be stabilized without controller [1, 4]. Here we performed three different controllers or rather loop-shaping design, 𝜇-synthesis, and sub(ğ»âˆž) robust control according to the mathematical model (11) with the assumption ̇𝐱1̇𝐱(𝑡)2(=âŽ¡âŽ¢âŽ¢âŽ£ğ‘¡)ğŸŽğˆ4×4𝐌−1𝐊𝑆4×4−𝐌−1î‚€Î©ğ†î‚âŽ¤âŽ¥âŽ¥âŽ¦îƒ¬ğ±1𝐱(𝑡)2(+âŽ¡âŽ¢âŽ¢âŽ£ğ‘¡)ğŸŽğŸŽ4×4−𝐏𝑀4×4−Ω𝐌−1ğğºâŽ¤âŽ¥âŽ¥âŽ¦îƒ¬ğ®ğ‘€ğ®ğºîƒ­+âŽ¡âŽ¢âŽ¢âŽ£ğŸŽ4×4𝐌−1ğŠğ¼âŽ¤âŽ¥âŽ¥âŽ¦îƒ¬ğ³ğ®(𝑡),1𝐳(𝑡)2=⎡⎢⎢⎢⎣(𝑡)𝐌−1𝐊𝑆4×4−𝐌−1î‚€Î©ğ†î‚ğŸŽî‚€Î©ğ†î‚4×4âŽ¤âŽ¥âŽ¥âŽ¥âŽ¦îƒ¬ğ±1𝐱(𝑡)2+⎡⎢⎢⎣(𝑡)−𝐏𝑀4×4−Ω𝐌−1ğğºğŸŽğŸŽ4×4âŽ¤âŽ¥âŽ¥âŽ¦îƒ¬ğ®ğ‘€ğ®ğºîƒ­+âŽ¡âŽ¢âŽ¢âŽ£ğŒâˆ’1ğŠğ¼ğŸŽ4×4⎤⎥⎥⎦𝐮(𝑡),𝐲(𝑡)sensor=𝐁Θdisp𝐁sensor4×4ğŸŽ4×4𝐱𝐓1(𝑡)𝐱𝐓2(𝑡)𝐓,(11)𝐱1𝑥(𝑡)=𝑔(𝑡)𝛾(𝑡)𝑦𝑔(𝑡)𝜗(𝑡)𝐓,𝐱2(𝑡)=̇𝑥𝑔(𝑡)̇𝛾(𝑡)̇𝑦𝑔̇(𝑡)𝜗(𝑡)𝐓,𝐮(𝑡)=𝐢𝑐(𝐮𝑡),𝑀𝐮𝐺=îƒ¬ğš«ğ‘€ğŸŽ4×4ğŸŽ4×4𝚫𝐺𝐳1𝐳(𝑡)2.(𝑡)(12)For all kind of robust control systems performed in this paper, a state space equation in a package form is adopted for this purpose. The package form characterized by the introduction of all inputs in terms of uncertainties, disturbances, and control signal is introduced into plant as in âŽ¡âŽ¢âŽ¢âŽ¢âŽ£ğ†=𝐀𝐁1𝐁2𝐂1𝐃11𝐃12𝐂2𝐃21𝐃22⎤⎥⎥⎥⎦,(13) where âŽ¡âŽ¢âŽ¢âŽ£ğ€=ğŸŽğˆ4×4𝐌−1𝐊𝑆4×4−𝐌−1î‚€Î©ğ†î‚âŽ¤âŽ¥âŽ¥âŽ¦,𝐁1=âŽ¡âŽ¢âŽ¢âŽ£ğŸŽğŸŽ4×4−𝐏𝑀4×4−𝐌−1𝐏𝐺⎤⎥⎥⎦,𝐁2=âŽ¡âŽ¢âŽ¢âŽ£ğŸŽ4×4𝐌−1ğŠğ¼âŽ¤âŽ¥âŽ¥âŽ¦,𝐂1=âŽ¡âŽ¢âŽ¢âŽ¢âŽ£î‚€ğŒâˆ’1𝐊𝑆4×4−𝐌−1î‚€Î©ğ†î‚ğŸŽî‚€Î©ğ†î‚4×4⎤⎥⎥⎥⎦,𝐃11=âŽ¡âŽ¢âŽ¢âŽ£âˆ’ğğ‘€4×4−𝐌−1ğğºğŸŽğŸŽ4×4⎤⎥⎥⎦,𝐃12=âŽ¡âŽ¢âŽ¢âŽ£ğŒâˆ’1ğŠğ¼ğŸŽ4×4⎤⎥⎥⎦,𝐂2=𝐁Θdisp𝐁sensor4×4ğŸŽ4×4,𝐃21=ğŸŽ4×8,𝐃22=ğŸŽ4×8.(14) Usually a rotor supported by active magnetic bearing needs to reach some desired performances that are described by weighting functions. The weighting functions introduced in the plant are relative to position and control signal performances in order to impose limits in the current value and maximum displacement of each rotor’s section. The weighting functions are transfer functions in Laplace domain describing the desired dynamic behavior of the rotor’s section varying by the excitation frequency. The expressions (15) and (16) are the weighting functions for position and control input, respectively,𝑊𝑝=𝑠2+1.8𝑠+10𝑠2,+8𝑠+0.01(15)𝑊𝑢=10−5.(16)

The weighting function described in (15) is applied to the displacements of the rotor’s sections, while the control signal is scaled according to (16). The transfer functions are taken by considering that the system has to be subjected to an attenuation of disturbances. This attenuation can be developed by maintaining the frequency response of the shaped plant by the controller well below the frequency response of (15) in the frequency range of interest. The expression of (15) can be obtained by MATLab commands by specifying the upper and lower limits of range of frequency. The block scheme is shown in the Figure 2, where 𝐺mds is the nominal plant meant without uncertainties, 𝐾 is the controller, 𝑑 the disturbances, and 𝑒𝑝 and 𝑒𝑢 the outputs of weighting function with regard to the position and control signal, respectively.

The presence of weighting functions produces an increase of state vector’s variables so that the new plant is P as shown in Figure 3.

All controllers used in this paper are characterized by a common concept or rather the robustness. The robustness is meant in a double way: robust stability and robust performance. The closed-loop system achieves robust stability if it is internally stable for all possible plant models 𝐺=𝐹(𝐺mds,Δ). In the present case, this means that the system must remain stable for any value of 𝛿𝑚,𝐼𝑃,𝐼𝑃.

Since those weighting functions are introduced in order to provide some characteristic on the system’s output, the robust performance criterion (17) is introduced for all 𝐺=𝐹(𝐺mds,Δ)‖‖‖‖𝑊𝑝(𝐼+𝐺𝐾)−1𝑊𝑢𝐾(𝐼+𝐺𝐾)−1‖‖‖‖∞<1,(17)îƒ¬ğžğ‘ğžğ‘¢îƒ­=𝐖𝑝(𝐈+𝐆𝐊)−1𝐖𝑢𝐊(𝐈+𝐆𝐊)−1𝐝,(18) where Δ can be structured or unstructured uncertainty matrix; in any case it must satisfy the norm condition or rather: ‖Δ‖∞<1 for robust performances and where vector d in (18) is the disturbances on the system.

Hence, the performance criterion is that the transfer functions from 𝐝 to ğžğ‘ and ğžğ‘¢ should be small in the sense of ‖‖∞for all possible uncertain transfer matrices Δ.

4. Results, Simulations, and Discussion

The simulations are performed by considering the data contained in Table 1.

Another set of data are referred to the transfer function introduced in the plant of our system. These transfer functions are essential if a certain performance must be obtained; these performances are usually referred in the frequency domain for example. Some authors usually introduce scalar weighting function in order to describe a certain constant value they want to obtain as a particular output. In this paper, we introduce a transfer function in the output signals or rather the displacement. The introduction of transfer function in Laplace domain “s” means that the displacement must be characterized by a certain dynamic behavior according the frequency variable. This technique is commonly used above all when a flexible structure is taken into account or when some nodes are subjected to vibrations such as in this case. Figures 4, 5, and 6 show the frequency response of weighting function relative to the displacement performances in order to analyze the sensitivity function. This is made for all controlled axis of each radial active magnetic bearing. Figure 4 shows the sensitivity function for the loop-shaping controller design, and we can see that the system has a good attenuation of disturbances until a certain value of frequency equals 103 Hz for Loop-Shaping Design Procedure (LSPD) controller.

Figures 5 and 6 show that 𝜇-synthesis and Sub(ğ»âˆž) robust control are not able to maintain a good rejection of disturbances for the entire range of frequency required for the system.

In a physical sense, the μ value represents an amplification coefficient of the effect produced by an exogenous excitation on the system. In a mathematical sense, it is determined by the following expression: 𝜇𝚫,𝐆𝑧𝑤=1minî€·ğœŽ(𝚫)∶det𝐈−𝐆𝑧𝑤𝚫,(19) where ğœŽ(Δ) is the singular value of relative input-output plant as a function of Δ.

Figure 7 discusses the robust stability condition for loop-shaping design procedure and shows that system has robust stability until a value of frequency of 103 Hz. After this value of frequency, the system has a pick outside the unitary value and goes back under one when the frequency approaches to the value of 104 Hz (upper and lower bound). The loop-shaping design shows good performances until 103 Hz for upper and lower bound since 𝜇<1. The loop-shaping design shows also a good rejection of disturbances as shown in the sensitivity function in Figure 4; in fact it has maintained itself well below the transfer function (15).

Figure 8 shows that robust stability for the 𝜇-synthesis is not maintained in the range of frequency near [100, 103] Hz so that the system does not have robust stability in that range. Figure 8 shows that, in the study of robust performances, the 𝜇 value is more than one for all range of operating frequency. This means the state of uncertainties introduced in this work induces the system to have bad robust performances also for the nominal plant. The same behavior is obtained by the Sub(ğ»âˆž) control.

Obviously the previous results affect the dynamic behavior of the entire system. According to Figure 3, the study of reference’s tracking and disturbances rejection is performed. It is built for all studied controllers. In Figure 10, the disturbance rejection test is performed according to a simulation characterized by a range of time of sixty seconds and a disturbance injection built as a square wave (black line) with a period of 20 s and amplitude of 10−6 m. All three controllers implemented are capable to support the requirements to reject the disturbance in a different way. The controllers Sub(ğ»âˆž) and 𝜇-Synthesis are characterized by the same pattern and the same period of oscillation, and they are perfectly superimposed. Their dynamic behavior is typical of damped system where a certain overshoot’s value is present and different from the previous controllers. The loop-shaping controller (red line) provides a good performance for the disturbance rejection due to the short period to extinct the transient response and small overshoot’s value if compared with those offered by the previous controllers; it has a good disturbance rejection and good performance because the it has a 𝜇 value lower than one. In all cases, the controllers are able to reject the disturbances by leading all suspended section to maintain the operating point position. Figure 11 shows the reference’s tracking simulations.

The controllers are characterized by state space representation (20) with four input positions and four control current outputs:̇𝐲𝐱(𝑡)=𝐀𝐱(𝑡)+𝐁𝐮(𝑡),(𝑡)=𝐂𝐱(𝑡),𝐮(𝑡)=𝐪(𝑡)sensor,𝐲(𝑡)=𝐢𝑐(𝑡).(20)

5. Conclusions

This paper shows that a comparison of three different control systems is built for a suspended rotor by active magnetic bearings. The comparison shows that loop-shaping design procedure provides the best performance to eliminate the disturbances and to follow the reference’s tracking when a certain required performances on position and control current are necessary. The present study shows that, once the weighting functions are introduced, only the loop-shaping design procedure is able to lead the system with robust stability and robust performance. The next development is to produce the previous discussion for the flexible rotor under the assumption that sensors are not colocated and with a state of uncertainties on displacement and current control gains.