# Research on the Improvement of Feedback Linearization Control in Suspension System Countering Inductance Variation

Authors: Liwei Zhang; Yue Zhang; Chao Zhang; He Zhao
Journal: Mathematical Problems in Engineering
 (2019)

---

## Abstract

The safety of the magnetic levitation (maglev) train is closely related to the control performance of the suspension module. However, during operation, the working conditions vary and are vulnerable to the external disturbances. In this work, a large-scale variation of the inductance of the magnetic levitation operation under different air gap conditions is considered, where the transfer function of the system changes nonlinearly. On the basis of the classical feedback linearization method, the algorithm of the first-order derivative for a single equilibrium point is improved, and then a multiequilibrium point feedback linearization method subject to the variation of the inductance is derived. The proposed linearization method can decouple the inductance from the air gap dynamics in any state of levitation, thus, reducing the model error. Using a general linear controller, the closed-loop control performance of the nonlinear hybrid excitation suspension system is run in MATLAB®. The simulation results show that the proposed method achieves good dynamic performance under various operating conditions and it improves the robust performance of the system.


---

## Body

## 1. Introduction

For urban low-speed magnetic levitation traffic, the hybrid excitation suspension electromagnet combines the permanent magnet excitation and the electric excitation and is connected in series in the magnetic circuit to generate a control voltage with a varying amplitude and direction through the winding to realize the magnetic circuit. The electric excitation control plays a role in adjusting the magnetic flux density, changing the magnetic attraction force, and adjusting the motion state of the system. Compared with the conventional electric excitation suspension magnet system, the most significant advantage of a hybrid excitation suspension system is energy savings, and the power required to the control system is small [1]. Studying hybrid excitation suspension system has important practical significance for improving the economic benefit and practical value of magnetic suspension technology.In recent years, many scholars have conducted in-depth researches and discussions on the hybrid excitation magnetic levitation technology. Among the hybrid excitation control methods at the present, the Taylor series expansion about an equilibrium point is commonly used in the modeling approach [1–3]. However, at the positions that lie far away from the equilibrium point, especially during the course of events such as levitation and landing, the relative stability of the control system drops significantly. With the large-scale application of differential geometry in electrical control theory, nonlinear control methods have emerged. References [3–6] apply the state feedback linearization to magnetic levitation systems, and the control performance is improved compared with Taylor’s linearization method. Among them, references [3, 4] use the current as the control input for being easily measured, but it does not represent the electromagnetic and dynamic changes inside the magnetic levitation system well. References [5, 6] use the air gap magnetic flux density as the internal control input that clearly describes the inherent law of the system. However, on this basis, it is necessary to install a magnetic density measuring device which is very rare in the field to form internal feedback of the system. Through the improvement of the above theoretical basis, references [7, 8] propose a feedback linearization control method based on the rated operating point, which is a first-order derivative transformation process for the nonlinear system at the rated operating point. The simulation showed that the method is more robust to the uncertain systems.Based on the feedback linearization method, references [9–13] robustly optimize the controller according to system characteristics and control requirements. In the design process, the control idea fully considers the uncertainties that appear inside and outside of the system, which include the system parameter perturbation, external noise, and high frequency interference. Combining the nominal and uncertain models with the robust controllers, the design optimization can be implemented for a given class of models. Compared with only considering the nominal model, this control method improves the conservativeness in the nominal design, but the system performance can be well-preserved under the parameter perturbation.For systems that contain nonlinearities or potential bounded time-varying uncertainties, or systems that do not satisfy the global matching conditions, [9] propose an adaptive robust control method with state constraints which is suitable for the nonlinear maglev train suspension control. The method consists of a three-step state transition that converts the maglev train into an interconnected uncertainty system. And then, based on the proposed robust control of transformed system, the overall uncertainty of the adaptive law simulation system is constructed.For the classes of fractional and integer-order systems with mismatched perturbations, a sliding mode control method based on a new fractional-order disturbance observer is proposed in [10]. This method can effectively handle the mismatched disturbance between quadcopters and the magnetic suspension systems, with better control performance, faster response, and reduced overshoot and jitter.Reference [11] proposes a real-time adaptive control method suitable for large-scale mass change of magnetic levitation system. The method uses a fast internal current loop combined with adaptive control to manipulate the suspension control system, and even when the load weight varies widely, the system is still in a stable state.In order to deal with the large-scale parameter perturbation changes, [12] performs inter-partition processing on the parameters. According to the experimental method, the weight functions of each segment are optimized separately to reduce the system conservativeness.In [13], a Takagi-Sugeno (T-S) fuzzy controller based on the improved form of the piecewise Lyapunov function is used to relax the stability condition of the model, and the proposed method is combined with the Parallel Distributed Compensation (PDC) controller to stabilize the position of the magnetic ball in the magnetic suspension system under the disturbance. In addition, this technology ensures that the system is robust, the disturbance tolerance is maximized, the performance is improved, and the H∞ performance metrics γ is reduced.Based on [5, 6], this paper establishes a mathematical model of hybrid excitation single magnet. The differential geometry transformation in [6, 7] is used to obtain the classical Brunovsky’s Canonical Form, and further first-order derivative operating is performed at the suspension equilibrium point to obtain a robust feedback linearization system. Based on the above linearization system, in the process of combining system air gap variation and the variation of the inductance value, the decoupling of system inductance value and the air gap change is proposed to improve the multiequilibrium point-based feedback linearization control method. The method can effectively cope with the model error caused by the variation of the system inductance and can improve the dynamic performance by real-time correction. Combining the linearization method proposed in this paper and the fundamental linear PID controller, the closed-loop control of the hybrid excitation suspension system is studied and verified through the MATLAB® simulation. In this simulation, the commands of levitation, square wave follower, and sudden load are given, and the common drift system parameter drift is simulated to investigate and compare the dynamic performance under different commands of each system. The simulation results demonstrate the effectiveness of the proposed linearization method.


## 2. Analytical Model

The research object of this paper is a hybrid excitation suspension module, which is mainly composed of long straight guide rail, electromagnet, and permanent magnet. The structure diagram is shown in Figure1. In recent years, the low-speed maglev trains designed and manufactured in China use U-shaped magnet structures, which can be decoupled from the linear motor propulsion structure and more flexible in train control [1].Figure 1
Schematic diagram of hybrid excitation.According to the working principle of the electromagnetic attraction, the air gap flux formula of the magnetic circuit of the hybrid excitation suspension module in Figure1 is(1)Φz,i=FmRm=μ0ANit+Hchpm2zt+hpm/μrIn Formula (1):  Φ is air gap flux; Fm is the total magnetomotive force;  Rm is the total reluctance; μ0 is the vacuum permeability; μr is the relative permeability;  A is the electromagnetic attraction cross-sectional area; z(t) is the variable air gap; i(t) is the electric excitation current; Hc is the permanent magnet nominal coercivity; hpm is the total thickness of the permanent magnet; N is the number of winding turns.The formula for the electromagnetic attraction force without considering the magnetic flux leakage effect when the magnetic circuit passes through the air gap is obtained as(2)Ft=Bt2Aμ0=μ0ANit+Hchpm2zt+hpm/μr2The system is controlled by the voltage input signalu(t). Assuming the circuit impedance isR, the circuit formula is obtained as(3)ut=Rit+ddtNΦz,i=RNμpmBt·2μrz+hpm-Hchpmμpm+NAB˙In Formula (3), μpm=μ0.μr is the permeability of the permanent magnet.The dynamics of the suspension module is mainly affected by the electromagnetic attraction forceF(z,i), the external interference  fd, and the vehicle weight Mg, where M is the vehicle mass and g is Earth’s gravity:(4)Md2ztdt2=-Fz,i+Mg+fdIn this paper, the external disturbance forcefd is mainly composed of normal attraction Fd. The permanent magnet synchronous linear motor above the suspension module will produce the interference effect under static or working conditions [14], so the suppression of this force is one of focal research points of this paper.Figure2 is the system open-loop control block diagram, obtained from Formulas (1) to (4).Figure 2
Open-loop control block diagram of the hybrid excitation suspension single iron system.


## 3. Feedback Linearization

The state-space variablesx=x1x2x3T=ztvtBtT are defined as the air gap, the vertical motion velocity, and the air gap flux density, respectively.The suspension module described by Formulas (1)-(4) is a nonlinear system, and the system can be described as a mathematical model of a standard affine nonlinear system in the form of the following formulas: (5)x˙=fx+gxuy=hx(6)fx=x2g-Aμ0Mx32RAμpmN2Hchpmμpm-x3·2μrx1+hpmgx=001NAThx=x1u=utSince the relative order of this system is 3, which is equal to the system dimension, therefore, the system satisfies the controllability and the matching conditions and can be accurately linearized. By solving the Lie derivatives according to Formula (6), a new control variable v  is obtained after the coordinate transformation. The control input voltage uc of the original nonlinear system can be expressed as the following formula:(7)uc=αcx+βcxv=-Lf3hx+vLgLf2hxNote thatv  is a new control variable and is not the same as vertical motion velocity v(t). αc(x) and βc(x) are defined in Formula (8), representing the constant term and the coefficient term, respectively, for the classical feedback linearization. In Formula (8), the other ϕc(x) is defined as the state variable vector group.(8)αcx=-Lf3hxLgLf2hx=-RμpmNHchpmμpm-x3·2μrx1+hpmβcx=1LgLf2hx=-μ0NM2x3ϕcx=x1x2g-Aμ0Mx32TThe original nonlinear system (5) can be converted to the Brunovsky canonical form, and its affine transformed state variable  ϕc is expressed as(9)z˙=010001000z+001vy=100zz=ϕcx=hxLfhxLf2hxT=x1x2g-Aμ0Mx32T.Note thatz is the transformed coordinate and is not the same as the air gap z(t). The nominal transfer function gs1 of system (9) can be obtained as (10)gs1=csI-A-1B=1s3.


## 4. Feedback Linearization Based on Rated Operating Point

During the implementation of the above control method, the nonlinear system is converted by the feedback linearization into a linear representation of the Brunovsky canonical form (9), which is a general result of any third-order SISO nonlinear system. In the subsequent model reference controller design, much of system information is hidden in the linearization transformation. Therefore, the feedback linearized model containing the physical information of the system must be obtained by means of the first-order derivative of the equilibrium operating point.The Brunovsky canonical form (9) is expressed as(11)z˙=Az+Bu.where the state matricesAr and  Br are the first-order derivatives of the system at the rated equilibrium point x0=[Z0 V0 B0]T:(12)Ar≜∂f∂xx=x0,Br≜gx0.Then the system control inputu is(13)ux,v=αx+βxvwhere (14)αx≜αcx+βcxYT-1ϕcxβx≜βcxH-1ϕx≜T-1ϕcxIn Formula (14), α,  β,  ϕ represent the constant term of the robust feedback linearization, the controller coefficient, term and the state variable vector group, respectively, and the variables are as follows:(15)G≜LgLf2hx0Y≜-G∂xαcx0T≜∂xϕcx0H≜G-1where G is the second-order Lie derivative, Y is the partial derivative of αc at the equilibrium point and the product of-G, T is the partial derivative of the state variable vector ϕc at the equilibrium point, and H is the reciprocal of G.According to the original nonlinear system, the feedback linearization results based on the equilibrium point x0=[Z0 V0 B0]T are given by(16)Ar=∂xfx0=01001-2AB0μ0M-2μrB0RAμpmN20-2μrZ0+hpmRAμpmN2,Br=001NAAccording to Formula (13), the actual control input for the system, Formulas (17) and (18) can be derived:(17)αx=αcx+βcxYT-1ϕcx=-RμpmNHchpmμpm-x3·2μrx1+hpm-2RB02μ0Nx1x3+μ0RM2μpmNA1x3g-Aμ0Mx322μrZ0+hpmβx=βcxH-1=-B0x3(18)ϕx=T-1ϕcx=x1x2-μ0M2AB0g-Aμ0Mx32According to the system state Formula (18), the transfer function gs2 is derived as (19)(19)gs2=CrSI-Ar-1Br+d=-2B0AμpmNμpmμ0MN2As3+RMμ02μrZ0+hpms2-4AB02RμrNow, the transfer functiongs2 containing the rich physical information of the system can be obtained, and the controller can be designed by the robust optimization method.


## 5. Multiequilibrium Point Feedback Linearization Method for Variation of Inductance Value of Suspension System

During the motion of the suspension system, the above feedback linearization system is based on the rated levitating air gap. Compared with the decoupled linearization system of the conventional motor, the inductance value changes nonlinearly in response to the levitating air gap, as shown in Formula (20) and Figure 3:(20)L=N2Aμ02zt+hpm/μrFigure 3
Relationship between single magnet suspension inductance and air gap.During the system levitation process, the air gap of the single magnet model is reduced from 0.025 m to the rated operating point of 0.01 m, and its inductance value is also changed from 0.17H to 0.03H. In other words, the suspension system is a strongly nonlinear system compared to the motor control system. The physical model of the system changes drastically throughout the start-up process. If the system only uses the transfer function based on the equilibrium point parameters, as Formula (19) shows, the deviation of the system characteristics will be generated at other operating points. Considering the above shortcomings, an improved multiequilibrium point feedback linearization method for the variation of the inductance value of the suspension system is obtained. This allows Formula (17) to be rewritten as the following formula:(21)αLx=-RμpmNHchpmμpm-x3·2μrx1+hpm-2RB02μ0Nx1x3+μ0RM2μpmNA1x3g-Aμ0Mx322μrzt+hpmβx=-B0x3The method replaces the rated operating point air gapZ0 in the  α term with a measurable air gap z(t). That is, at any time, ignoring the influence of nonlinear factors such as leakage inductance, the feedback linearization method can reproduce better dynamic characteristics of the original single-magnet nonlinear suspension system. The transfer function gs3 of the new feedback linearization method is shown in the following formula:(22)gs3=CrSI-Ar-1Br+d=-2B0AμpmNμpmμ0MN2As3+RMμ02μrzt+hpms2-4AB02RμrThe system’s fixed value ofB0 is only related to the load weight, independent of the system air gap, and the inductance is completely decoupled from the levitating state and position. Based on this method, a linear controller is designed, which can be added into the total controller of the nonlinear hybrid excitation suspension system. This particular closed-loop control block diagram is shown in Figure 4.Figure 4
Improved robust feedback linearization method suspension system control block diagram.


## 6. System Stability Analysis

In order to verify the stability of the feedback linearization system combined with the nonlinear magnetic levitation system, this paper carries out the following studies.In Figure4, the nonlinear magnetic levitation system processed by a robust feedback linearization method with inductance value decoupling can be expressed as a linear system, as shown in the following formula:(23)x˙=Arx+BruThe parameters of the maglev module are shown in Table1.Table 1
The parameters of the maglev module.
Parameter  Value Overall mass of the suspension M (Kg) 650 Electric excitation winding number N 200 Permanent magnet coercivity H C (Am) 8.8×105 Vacuum conductivity μ 0 (H/m) 4π×10−7 Permanent magnet relative permeability μ r 1.05 Electric field winding resistance R (Ω) 0.5 Permanent magnet radial length h p m (m) 0.02 Permanent magnet cross-sectional area A (m2) 0.0235 Gravity acceleration g (m/s2) 9.8According the Formula (16) and the parameters shown in Table 1, the specific expression of Ar,Br can be obtained.(24)Ar=∂xfx0=01001-33.62-494.650-16.54,Br=004.7The state feedback control is represented byk1,k2, andk3, representing the feedback gains of the system state space variable system air gapx1, velocityx2, and accelerationx3, respectively:(25)u=k1k2k3x1x2x3=k1x1+k2x2+k3x3Let the closed-loop system characteristic formula equal to 0.(26)det⁡Ar+Bru=0The following conditions can be gotten(27)k1>-2324.855k2>0k3<-77.738That is, if the state feedback control parameter satisfies Formula (27), the closed loop system shown as Figure 4 is stable.When the proposed feedback linearization method decouples the air gap inductance, the linearized system equation of the system will change monotonously with the air gap, and the parameter stabilization condition of the state feedback controller will change slightly. The control parameters used in this paper meet the stability requirements within the range of the model physical model.The system is also analyzed by the Linear Analysis module of the MATLAB®/Simulink® software. The system transfer functions are shown in Table2. The original nonlinear system with the equilibrium points of 0.01m, 0.0075m, and 0.0125m and the improved feedback linearization system Bode diagram are shown in Figures 5 and 6. The results show that the improved feedback linearization method can reproduce better frequency domain characteristics of the original nonlinear system.Table 2
System transfer function.
Desired Air gap (m) Original system Improved feedback linearization system 0.01 -7.145s3+16.53s2-1.659e4 - 7.15 s 3 + 16.54 s 2 - 1.661 e 4 0.0075 -7.145s3+14.41s2-1.659e4 - 7.15 s 3 + 14.42 s 2 - 1.661 e 4 0.0125 -7.145s3+18.65s2-1.659e4 - 7.15 s 3 + 18.67 s 2 - 1.661 e 4Figure 5
Original nonlinear system Bode diagram.Figure 6
Improved feedback linearization system Bode diagram.


## 7. Simulation

Based on the MATLAB®/Simulink® simulation platform, this paper simulates the working process of a magnetic hybrid excitation suspension module. The parameters of the system are shown in Table1. The parameters of the PID controller are 10000, 10000, and 1500, respectively. The acceleration state feedback gain is -2600. To validate the algorithm, the same state feedback controller is used for the following four control systems:ORI: the original nonlinear control system,CFB: the classical feedback linearized control system,RFB: the equilibrium point-based feedback linearized control system,LRFB: the feedback linearized control system based on the improvement of the inductance value.The control commands include the system response from the levitation to the rated operating point, the tracking of the operating point, and the sudden load, as shown in Figures7, 8, and 9.Figure 7
Levitation of the suspension system to the rated operating point.Figure 8
The suspension system tracks the rated working points.Figure 9
Response of a sudden load of the suspension system.In Figure7, the system levitates from the initial position to the equilibrium point of 0.01 m.RFB andLRFB reach the steady-state with a small overshoot after 1 s, whileORI andCFB have a large overshoot (0.0075 meters) and reach the steady-state at 2.5 s after five damping oscillations.In Figure8, all four systems track the setpoints command from the equilibrium position (0.01 m) with a 1 Hz, a peak of 0.0125 meters, and a valley of 0.0075 m for motion. The responses ofORI,RFB, andLRFB follow the dynamic similarity of the system, leading to the amplitude of the levitation action decreasing in turn. Under the action ofCFB, that is, the classical feedback linearization method, larger overshoots are generated during the even-number of overshoots, which is distinguished from the dynamic characteristics of the other three systems. The reason for the explanationCFB result is thatORI,RFB, andLRFB have similar system models, whileCFB transforms the original system into gs=1/s3, which weakens the dynamic control performance of the system to some extents.Figure9 shows the response of the system to a sudden load of 0.2 kg under equilibrium. Similar to the above simulation,ORI,RFB, andLRFB have similar disturbance rejection characteristics. After 1s of transient time,ORI,RFB, andLRFB return to the equilibrium position, and the overshooting characteristic ofLRFB is optimal (which means the lowest). After passing 2 s from 3.2 s,CFB achieves the maximum overshoot that is 1.5 times of the maximum overshoot ofLRFB. However, CFB oscillates many times (only 2.5 oscillations are shown in Figure 9) before readjusting back to the equilibrium point.The reason for the explanation of CFB’s dynamical behavior is that when there is a sudden load, the model’s system transfer function changes. Since the linearized models under the action ofRFB andLRFB are closer to the nonlinear model after the load changes, then this explains the dynamical behaviors ofRFB andLRFB. Similar to theRFB, theLRFB method can deal with the system disturbance that is caused by the mismatch of the model, better than the others. The mismatch happens when the magnetic levitation attraction is caused by the load variation, the external disturbance of gravity direction, or the leakage inductance. This verifies the system’s robustness.In order simulate the various conditions encountered in actual operation, a combination of parameter variations involving different types and states are applied. This helps to further verify the robust performance of the system. The simulation uses the sampling method for uncertain parameter combination that is specified by human operators. This method combines the uncertain parameter in the traversal way. Next, the simulation experiment is conducted in Simulink® repeatedly. Then, the results are summarized in the same observation coordinate system, and thereby, achieving the purpose of simulated random sampling. There are three parameters considered: the armature winding resistance R, which is greatly affected by the temperature; the coercive force Hc of the permanent magnetic physical patch inside the hybrid excitation structure; and the levitated overall mass M.Each parameter is taken one nominal quantity and two deviation quantities. Thus, there are a total of 27 parametric combinations of 33 variations.The performances of the four different feedback linearization methods are analyzed after summarizing the effects of various possibilities on the dynamic responses of the magnetic levitation control system. In order to distinguish the variations of the suspension mass significantly up a certain degree, the specified ‘nominal’ massM1 is labelled in Blue line, the total massM2 (heavier thanM1 by 20%) is labelled in Red line, and the total massM3 (lighter thanM1 by 20%) is labelled in Black line (see Figure 10). The specific variation parameters are shown in Table 3.Table 3
The parameters of maglev module.
Parameter  value Overall mass of suspension M1 (Kg) 650 Overall mass of suspension M2 (Kg) 780 Overall mass of suspension M3 (Kg) 520 Electric field winding resistance R1 (Ω) 0.5 Electric field winding resistance R2 (Ω) 0.525 Electric field winding resistance R3 (Ω) 0.55 Permanent magnet coercivity Hc1 (Am) 8.8×105 Permanent magnet coercivity Hc2 (Am) 8.36×105 Permanent magnet coercivity Hc3 (Am) 7.92×105Figure 10
Dynamic response of system floating and step disturbance under multiple operating conditions.Figure10 shows the dynamic response curves of the four linearization methods that are combined with the same linear controller for levitation and sudden load under multiple operating conditions. It can be seen that bothRFB andLRFB maintain the consistency of the dynamic curves, and their overshoots are relatively small. However,CFB shows large inconsistencies of dynamical responses under different system parameter combinations, and their transient times are relatively longer. The maximum undershoot ofCFB is about 25%. By directly controlling the suspension system with a linear controller, only two model curves of the specified mass achieve the nominal performances, and the others produce large undershoot.Figure11 shows the dynamic responses of the suspension system subjecting to the square wave command input under multiple operating conditions. In these cases,CFB also exhibits poor dynamic performance with the largest overshoot and the longest transient period. Comparatively,ORI, RFB, andLRFB perform better in tracking the square wave signal, but among them, dynamic responses ofLRFB are the most definitive. The reason is that the upper and lower amplitudes of the square wave signal do not lie at the rated operating point, which enables theLRFB to change the feedback coefficient according to the actual position of the system. This linearization method can reflect better system characteristics at the nonrated points. In particular, the inductance of the levitated magnet is decoupled to avoid modeling errors caused by different positions of the levitation.Figure 11
Dynamic response of system following square wave under multiple operating conditions.


## 8. Conclusion

This paper proposes a robust feedback linearization method based on the improved variations of inductance for the hybrid excitation magnetic levitation system. The proposed method is built on the basis of various control methods of feedback linearization in the reference. This method considers three things, namely, the levitation excitation of the hybrid excitation magnet, the different operating conditions, and the sensitivity of the system model to the inductance parameter error. This method can further improve the linearization of the system, making it closer to the dynamics of the original nonlinear system, therefore, enabling a more robust linear controller to be designed.On the general basis of the linear controller and the same system modeling conditions, three different linearization methods are used to simulate the original nonlinear system. The results show that even in the presence of the internal system parameter drift and the external disturbances, the linearization method combined with the basic linear controller proposed in this paper can still maintain a nominal performance of the system. In other words, the proposed method has a certain degree of robust optimization effect on the linear controller of the front end. In this paper, the proposed method is beneficial to the enhancement of the system model details through the high precision feedback measurement of the state variables and the regulation of nominal system performance requirements. In addition, if the total mass of the suspension can be accurately measured, the dynamic performance of the system under heavy load or sudden load can be better guaranteed.



---

*Source: 5747812-2019-07-01.xml*