The Active Fractional Order Control for Maglev Suspension System
Maglev suspension system is the core part of maglev train. In the practical application, the load uncertainties, inherent nonlinearity, and misalignment between sensors and actuators are the main issues that should be solved carefully. In order to design a suitable controller, the attention is paid to the fractional order controller. Firstly, the mathematical model of a single electromagnetic suspension unit is derived. Then, considering the limitation of the traditional PD controller adaptation, the fractional order controller is developed to obtain more excellent suspension specifications and robust performance. In reality, the nonlinearity affects the structure and the precision of the model after linearization, which will degrade the dynamic performance. So, a fractional order controller is addressed to eliminate the disturbance by adjusting the parameters which are added by the fractional order controller. Furthermore, the controller based on LQR is employed to compare with the fractional order controller. Finally, the performance of them is discussed by simulation. The results illustrated the validity of the fractional order controller.
The developing of maglev train is under a rapid speed around the world. Thanks to no physic contact with guide-way, the maglev system gets the advantage of lower noise, less costly maintenance, less exhaust fumes emission, and so on, which is suitable for the urban transportation. Until now, two commercial maglev routes and several test routes have been built around the world.
Suspension control is one of the core technologies in maglev train [1, 2]. The practical model of suspension system is a nonlinear and inherent unstable system due to the fact that the electromagnetic force produced by a constant current is inversely proportional to the square of the levitation gap. Besides, the parameters of the levitation system are varying with the difference of the operating condition. So, it is hard to choose a practicable control strategy to design an ideal controller that can keep satisfied dynamic performance when the train is running under the actual environment. Considering the complicated working situation, many researchers have investigated many new controllers with different theories for the maglev suspension system such as adaptive control , nonlinear control , fuzzy control , and robust control .
The quick development of fractional calculus attracts extensive attentions due to its nonlinear feature and unique memory characteristic. In electricity, properties between resistance and capacitance are intermediated by a concept of fractance . Fractional calculus has been applied in the modeling and control of various kinds of physical systems and many real systems are modeled or fitted by fractional order systems. So, in recent years, with the progress of fractional order integral theory, many papers have discussed the design and application of controller design based on fractional order theory.
Fractional order controller attracts increasing attention because of the higher freedom degree provided [8–11]. In paper  a controller is proposed; this control scheme can improve the robustness of system against disturbances. In paper , Oustaloup discussed the controller. In paper , Podlubny focused on the controller; this method is the milestone of the ; it increased two parameters for the controller, so, with the controller, the system can get desired performance. Since then, many works have been done to apply these new methods to different nonlinear systems.
The maglev train, as a typical nonlinear system, is an ideal model for fractional order controller design. This essay paid attention to discussing the usage of fractional order control theory and application of the controller on solving the system inherent nonlinearity and parameters uncertain problems of the suspension system. The rest of the paper is organized as follows. The mathematical model of the basic maglev suspension system is developed in Section 2. In Section 3, the design theory and procedure of controller based on the simplified model are discussed. In Section 4, simulation is proposed, and we compared the performance with controller based on . Finally, conclusions are presented in Section 5.
A whole maglev train includes the air springs, actively controlled electromagnets, and the cabin. An overall view of the maglev vehicle is shown in Figure 1. The maglev train is made up of ten levitation modules, and each module can be decomposed into two single electromagnetic suspension units which are shown in Figure 2. To derive the dynamic model of the single electromagnetic suspension unit, from the opinion of engineering, several hypotheses were made. Firstly, it is assumed that the mass of cabin is transferred to the bogie only in the longitudinal direction. Secondly, neglecting the flux leakage of the electromagnetic due to the operating levitation gap is 9 mm. Finally, considering the working gap and engineering analysis, the magnetic resistances of the ferrite core and rail are neglected; namely, the magnetic potential falls on the air gap evenly.
The fundamental principle of the electromagnet suspension system (shown in Figure 2) is described as the following: according to the desired suspension gap and status information gotten by position sensor and accelerometer, the controller adjusts the current supply so that the suspension system can produce suitable force that can overcome the gravitation and keep the system stable in the expected operating position. A single electromagnet suspension unit can be simplified as shown in Figure 2.
In Figure 2, is the total mass of suspension system. denotes the number of turns of a single electromagnet. is the area of the magnetic pole. is the magnetic resistance. is the magnetic permeability of atmosphere. is the gravity acceleration. The magnetic inductance is defined as . is the force produced by electromagnetic for suspension. is the measured position value from position sensor. means the value of measured acceleration from accelerometers. The voltage supplied to magnetic suspension system is expressed by . The symbol is the current of magnetic suspension system. is the flux density of the gap. is the disturbance of the load.
There are three main relations in the suspension system; there are voltage () to current (), acceleration () to electromagnetic force (), and current () to electromagnetic force (). Based on the electromagnetic theory and dynamics analysis of electromagnet , the following equations can be built:where , , and .
There is a need to pay attention to the fact that this model did not put the guide-way irregularities into consideration. So, with , is equal to . From (1), it is obvious that the system is internally nonlinear, and it will be discussed later. Furthermore, simplifying the above equations, the following equation can be deduced:
From (2), it is obvious that electromagnetic force is dependent on the current and gap ; the relationship between them is nonlinear. That means there is nonlinear relationship between them and it is the source of the system nonlinearity.
Traditionally, the model is linearized in the general working point , although this will bring the loss of precision, but it is still an effective way to analyse the main feature of the system. Certainly, it is a good idea to design the controller based on the nonlinear model, but it is necessary to consider the difficulty brought for the model building and difficulty of controller designing.
The linearized system can be described aswhere , . Replacing , , and with , , and and defining that , , and , the state-space representation of the above equation iswherewhere includes the disturbance of from track, cabin, and other environment factors. From (4), it is obvious that the value of parameters is changing according to the working position . Furthermore, the disturbance is varying with the working environment at the same time. Facing so much uncertainty, as a good controller for the suspension system, it not only should keep the system stable in various working positions, but also has to reject the disturbance from environment.
3. Procedure of Controller Design
In this section, firstly, the classical controller design based on the method of is recommended . Then, the procedure of fractional order controller design is introduced based on the principle of designing a fractional order controller and the suspension system model built in Section 2.
3.1. Classical Controller Based on
There are many ways to design classical controllers. This part focused on the design of controller.
For a state-space equation,
is a state matrix. is a control matrix. is a output matrix. , , and are , , and time varying matrices. The quadratic performance index of the equation can be expressed aswhere , are nonnegative definite matrix, , are weighting matrix of state-variable and input-variable separately, is positive definite matrix, is a fixed time, and is free. The aim of the method is to be sure that the quadratic performance index should be the smallest, which means, in the process of control, the error and the consumption of energy are the least. The feedback control value is
Firstly, a Hamilton function is established:
Then, making a derivation to , the smallest value can be calculated as
Following the above steps, a classical controller can be designed with desired performance.
3.2. Fractional Order Controller Design
Define a system transfer function and its fractional order controller and is the open-loop transfer function. For the purpose of designing a fractional order controller, there are three specifications needed to be met [14–16].(a)Phase margin specification: is the gain crossover frequency interested, is the phase margin, and is the open-loop transfer function of the system.(b)Flat phase specification: The phase Bode is flat at the gain crossover frequency. It means that the system is insensitive to the gain changes in the crossover frequency.(c)Gain crossover frequency specification:
3.2.1. Preliminary Knowledge
For a fractional order controllerwhere and and are the relevant parameters needed to be calculated. The phase and gain of the controller can be expressed as
According to (4), the suspension system transform function is given as
Assuming , , . Then it is similar to that the phase and gain of are
3.2.2. Design of Fractional Order Controller
The open-loop transform function is
By applying the three specifications introduced above to (19), we get
Considering principle (a) that , then
Considering principle (b),
So, combining (22) and (23), it is easy to deduce and . It is also a good idea to get and by plotting a picture of the relationship of and . At the same time, considering principle (c), can be calculated as
3.3. Approximation of the Fractional Order Operator
Due to the fractional order differentiator, the fractional order controller is infinite dimensional. In practical application, it is crucial for us to get the finite dimensional approximation of the fractional order system. Many researchers have investigated many ways for this problem, in which the Oustaloup Recursive Algorithm  is a good way for this approximation. Choosing the interested frequency range, (), then the fractional order operator can be approximated with Oustaloup Algorithm as where the parameters in (25) can be evaluated by the following equations:
By using the Oustaloup Recursive Algorithm, a lot of convenience is brought.
4. Simulation and Analysis
4.1. Controller Designs
In order to discuss the fractional order controller design based on Section 3, this section investigated the performance of fractional order controller by simulation.
This paper chose a single standard maglev bogie as the model for study. The values of the basic parameters are listed in Table 1.
In practical, the input of the closed-loop system is the expected gap of the electromagnet. So the block diagram of the system can be shown in Figure 3.
By putting these parameters in the model, the mathematical model of suspension system in the working position can be expressed asand the transfer function of the system can be expressed as
The controller and fractional order controller can be designed according to Section 3.
The controller based on , by choosing , , and , and the feedback matrix can be calculated: . Certainly, when using the parameters, there is a need to do some fine adjustments.
The controller design based on fractional order theory is according to Section 3. And the designed controller can be expressed as . Figure 4 is the closed-loop system Bode diagram with fractional order controller.
From Figure 4, it is easy to know that the overshoot of the closed-loop system is very small. Bandwidth is about 15 Hz, which can meet the requirement of the system. So, for these principles, the controller design based on fractional order theory can satisfy the basic demand of the maglev system.
4.2. Simulation Comparisons
This part compared the performance of the closed-loop system with different controllers under different working conditions.
At first, this part compared the step response of the two closed-loop systems. In the step response, the overshoot and rising time are the crucial factors we were concerned about. Figures 5 and 6 showed the step responses of the closed-loop system with controller based on and controller based on fractional order theory, respectively.
Secondly, considering that the load of maglev suspension system is varying when operating, we hope the system can get stable performance when load is varying seriously. So, this section compared the step responses of the systems while load is changing. Figures 7 and 8 show the step response of different controllers when the load is varying.
At last, when maglev train is running, the working gap is not constant. It is easy to know that when working position is varying, nearly all the parameters of the system will change. This phenomenon will bring the uncertain to the system and degrade the performance of the controller. So, much attention is paid to this problem. Figures 9 and 10 showed the step response of different controllers under different work positions.
In Figure 5, the rising time is about 0.18 s and the overshoot is about 5%. In Figure 6, the rising time is 0.05 s, much shorter than in Figure 5, and the overshoot can nearly be neglected. Comparing Figures 5 and 6, obviously, the system with fractional order controller has better performance, not only in overshoot, but also in rising time and settling time. The fractional order controller, with the characteristic of nonlinear, improved the system performance of rising time.
In Figure 7, when the load changes from 350 kg to 1250 kg, not only the rising time, but also the overshoot of the step response varies seriously. And the heavier the load is, the longer the rising time is and the larger the overshoot is. When kg, compared with kg, the overshoot is more than 20%. The rising time, compared with kg, is about 0.08 s longer. Similar variables are applied to the system with fractional order controller and Figure 8 shows the response of it.
In Figure 8, the overshoots are smaller than those in Figure 7. When kg, the overshoot is less than 10%, and the rising time is less than 0.05 s. When changed to 653 kg, the overshoot and rising time are varying, but they are smaller than the response controlled by controller all the same.
Figures 7 and 8 present the step response when the load is varying. In contrast, in Figure 8, although the load is varying as in Figure 7, the overshoot and rising time of the step response are changing much smaller than in Figure 7. That means the system with fractional order controller is robust to load varying.
There is also a need to pay attention to the fact that, in Figure 8, when is varying, the difference of the response exists.
In Figure 9, when the working position changes from 0.006 m to 0.012 m, the rising time and the overshoot of the step response varied with the changing of working position.
In Figure 10, the overshoot and rising time remain constant when working position is varying. On the other hand, in Figure 9, the rising time and overshoot are unstable and changing seriously. The system performance with controller is sensitive to the varying of working position. But, in Figure 10, it is obvious that system with fractional order controller is not sensitive to the change of working position and is robust to this kind of parameter change.
The conclusion can be concluded that the fractional order controller can suppress the disturbance brought by working position varying remarkably.
(1)This paper applied a fractional order controller on the maglev suspension system and introduced the procedure of the controller design. Three specifications are discussed for the selection of parameters of the controller.(2)Thanks to its nonlinear feature and higher degree free, the fractional order controller, compared with classical controller based on , can achieve better dynamic performance and robustness for the system, especially to the uncertain load and operating gap. Simulation results show the advantages of the fractional order controller.(3)In the next stage, the main work should be focused on the realization of the controller on hardware for the usage in practical suspension system and the building of a reliable fractional order model for the suspension system.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
The authors gratefully acknowledge financial support supported by the National Natural Science Foundation of China (no. 11202230 and no. 11302252).
P. K. Sinha, Electromagnetic Suspension: Dynamics & Control, Short Run Press, London, UK, 1987.
L. Yung, D. Zhang, and C. Hu, “The design and simulation of an adaptive maglev control algorithm based on oscillation observation,” in Proceedings of the 18th International Conference on Magnetically Levitated Systems and Linear Drives (MAGLEV '04), vol. 2, pp. 984–990, October 2004.View at: Google Scholar
Three-parameter tunable tilt integral derivative (TID) controller, US Patent US5371670, 1994.
I. Podlubny, Fractional Differential Equations, Academic Press, 1999.View at: MathSciNet