Table of Contents Author Guidelines Submit a Manuscript
Modelling and Simulation in Engineering
Volumeย 2011, Article IDย 414286, 10 pages
Research Article

Performances Comparison for a Rotating Shaft Suspended by 4-Axis Radial Active Magnetic Bearings via ๐œ‡-Synthesis, Loop-Shaping Design, and Sub(๐ป)โˆžwith Uncertainties

Dipartimento di Ingegneria Chimica, Gestionale, Informatica e Meccanica, Universita' di Palermo, Viale delle Scienze, 90128 Palermo, Italy

Received 3 March 2011; Revised 8 April 2011; Accepted 9 May 2011

Academic Editor: Jing-songย Hong

Copyright ยฉ 2011 G. Barbaraci and G. Virzi' Mariotti. 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.


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.

Figure 1: Schematic view of 4-axis rotating shaft supported by two radial active magnetic bearings with sensors.

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.

Figure 2: Block scheme of plant with the introduction of weighting functions as further outputs.

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.

Figure 3: Block scheme of plant showing the new plant.

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.

Table 1: Data for simulation.

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.

Figure 4: Sensitivity function describing how the disturbance affects the displacement with LSDP controller.
Figure 5: Sensitivity function describing how the disturbance affects the displacement with ๐œ‡-synthesis controller.
Figure 6: Sensitivity function describing how the disturbance affects the displacement with Sub(๐ปโˆž) 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 7: Robust stability, nominal and robust performance function described by ๐œ‡ value with LSDP controller.

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.

Figure 8: Robust stability, nominal and robust performance function described by ๐œ‡-value with ๐œ‡-Synthesis.
Figure 9: Robust stability, nominal and robust performance function described by ๐œ‡-value with Sub(๐ปโˆž).

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.

Figure 10: Comparison via simulation for disturbance rejection test.
Figure 11: Comparison via simulation for reference tracking test.

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.


  1. G. Schweitzer, H. Bleuler, and A. Traxler, Active Magnetic Bearings, Vdf Hochschulverlag AG, Zurich, Switzerland, 1994.
  2. G. Genta, Dynamics of Rotating Systems, Mechanical Engineering Series, Springer, Torino, Italy, 2005.
  3. G. Barbaraci and G. Virzì Mariotti, โ€œControllo subottimo per un albero rotante in levitazione magnetica attiva,โ€ in Proceedings of the 38th Convegno Nazionale AIAS, Torino, Italy, Settembre 2009.
  4. G. Schweitzer and E. H. Maslen, Magnetic Bearings: Theory, Design and Application to Rotating Machinery, Springer, Zurich, Switzerland, 2009.
  5. G. Barbaraci, A. H. Pesch, and J. T. Sawicki, โ€œExperimental investigations of minimum power consumption optimal control for variable speed amb rotor,โ€ in Proceedings of the ASME 2010 International Mechanical Engineering Congress & Exposition (IMECE '10), pp. 1โ€“10, Vancouver, Canada.
  6. D. Bugajski, D. Enns, and A. Tannenbaum, โ€œSynthesis methods for robust nonlinear control,โ€ in Proceedings of the American Control Conference, pp. 531โ€“535, June 1993. View at Scopus
  7. C. R. Knospe, R. L. Fittro, and L. S. Stephens, โ€œControl of a high speed machining spindle via μ-synthesis,โ€ in Proceedings of the 1997 IEEE International Conference on Control Applications, pp. 912โ€“917, Hartford, Conn, USA, October 1997. View at Scopus
  8. J. D. Wu and J. H. Lin, โ€œImplementation of an active vibration controller for gear-set shaft using μ-analysis,โ€ Journal of Sound and Vibration, vol. 281, no. 3–5, pp. 1037โ€“1055, 2005. View at Publisher ยท View at Google Scholar ยท View at Scopus
  9. S. Kaitwanidvilai and M. Parnichkun, โ€œDesign of structured controller satisfying H loop shaping using evolutionary optimization: application to a pneumatic robot arm,โ€ Engineering Letters, vol. 16, no. 2, 2003. View at Google Scholar
  10. D. McFarlane and K. Glover, โ€œA loop-shaping design procedure using H synthesis,โ€ IEEE Transactions on Automatic Control, vol. 37, no. 6, pp. 759โ€“769, 1992. View at Publisher ยท View at Google Scholar ยท View at MathSciNet ยท View at Scopus