# Adaptive Control of Electromagnetic Suspension System by HOPF Bifurcation

Authors: Aming Hao; Xiaolong Li; Longhua She
Journal: Mathematical Problems in Engineering
 (2013)

---

## Abstract

EMS-type maglev system is essentially nonlinear and unstable. It is complicated  to design a stable controller for maglev system which is under large-scale disturbance and parameter variance. Theory analysis expresses that this phenomenon corresponds to a HOPF bifurcation in mathematical model. An adaptive control law which adjusts the PID control parameters is given in this paper according to HOPF bifurcation theory. Through identification of the levitated mass, the controller adjusts the feedback coefficient to make the system far from the HOPF bifurcation point and maintain the stability of the maglev system. Simulation result indicates that adjusting proportion gain parameter using this method can extend the state stability range of maglev system and avoid the self-excited vibration efficiently.


---

## Body

## 1. Introduction

With the success of the Shanghai Pudong International Airport link in China, magnetic levitation vehicle (maglev vehicle) through electromagnetic suspension (EMS) technology has come into our life [1, 2]. The EMS provides nonmechanical contacting suspension force by means of electromagnets with an air gap controller realized through position, velocity, and acceleration feedback. Classical control theories have provided much of the design rules for many EMS including maglev vehicle [3, 4] and magnetic suspension and balance system [5]. However, because the constant linear control strategies deteriorate rapidly with increasing deviations from the nominal operating point caused by significant changes in the suspended load, some approaches to the problem of ensuring consistent performance while operating point changes have been reported. A model reference adaptive controller to compensate for payload variations and external force disturbances has been given [6]. A gain scheduling method has been designed [7], where the nonlinear force/current/air-gap relationship of the magnetic suspension is successively linearized at various operating points with a suitable controller is designed for each of these operating points. Some other kinds of method also have been developed by [8–11].However, the above papers have not discussed how to deal with this case in which the mass levitated changes abruptly. Because the EMS is inherently unstable, nonlinear dynamic system, the above methods may not guarantee its stability with abrupt changing mass levitated sometimes. When the EMS becomes unstable by an abrupt disturbance, it may produce vibration firstly corresponding to HOPF bifurcation of a nonlinear state space equation [12]. That is to say, when designing the EMS controller parameters, if we consider an additional constraint which makes the controller parameters always far away from the HOPF bifurcation points while the mass levitated varies, a larger state stable range will be attained. Then, the anti-abrupt-disturbance ability of the system will be enhanced. This paper studies how to design parameter self-tuning adaptive controller according to the above principle.


## 2. One Degree of Freedom EMS Model

The one degree of freedom EMS is shown schematically in Figure1. It is composed of an electromagnet, an air gap sensor, a DSP controller, a current driver, and a levitated object. Through proper control algorithm in DSP, the object will be levitated stably with certain gap.Figure 1
Single degree of freedom EMS.Set the levitation air gap with the direction down ass, the expected air gap is s0, the velocity of air gap changing is y, the gravity acceleration is g, the mass of the object levitated is m,  the current in electromagnet is I, the resistance of the electromagnet is R,  and the electromagnet voltage on the system equilibrium point is U0.Set the turns of the electromagnet coil asN,  pole area of the electromagnet is A,  and permeability of vacuum is μ0. Then, set k=(μ0N2A)/4  as the constant of electromagnet determined by N, A, and μ0.Suppose that the following dual-loop controller acts on the system:(1)Uk=KP(c-c0)+KDc˙+KII,
where KP and KD are adjustable gap loop feedback parameters. KP is the proportional feedback coefficient, and KD is the differential feedback coefficient.K
I is the current loop feedback coefficient. By adopting the proper current feedback, the time constant of electromagnet can be greatly reduced. So in the band of Ems system, the transfer function from voltage to the current of electromagnet can be seen as proportional component and the controller design becomes more convenient [4].Then the the system model can be written as follows [13]:
(2)c˙=y,y˙=g-kI2mc2,I˙=-(c2kR-yc)I+c2k[U0+KII+KP(c-c0)+KDy].


## 3. System Analysis

Settings˙=0, y˙=0, and I˙=0, the equilibrium point of the state variables can be obtained:
(3)s0=I0kmg,y0=0,I0=U0R-KI.SetR′=R-KI.The Jacobian matrix of the system at the equilibrium point is(4)A(KP,KD)=[0102gs00-2s0kgmKPs02ks0KD2k+mgk-R′s02k].
The eigen polynomial of the Jacobian matrix of the system at the equilibrium point is
(5)|λI-A|=λ3+R′s02kλ2+KDgkmλ+KPgkm-gkR′.
The corresponding Routh table can be computed as in Table 1.Table 1
Routh table of EMS Jacobian matrix.
λ
3
c
1,1
=
1
c
1,2
=
K
D
g
k
m
λ
2
c
2,1
=
R
s
0
2
k
c
2,2
=
K
P
g
k
m
-
g
k
R
λ
1
c
3,1
=
2
s
0
m
g
k
(
m
g
k
+
K
D
s
0
2
k
-
K
P
R
)
λ
0
c
4,1
=
K
P
g
k
m
-
g
k
RIt is easy to know according to Routh table that whenR′mg/k<KP<(R′s0/2k)KD+R′mg/k, the system (2) is stable at the equilibrium point.Furthermore, there is a bifurcation whenc4,1=0 or c3,1=0. In fact,  KP>R′mg/k, so c4,1=0 is impossible. This means that there is a bifurcation in the system only when c3,1=0. Next, we will show the bifurcation is HOPF bifurcation.Whenc3,1=0, there is a pair of pure imaginary roots and a negative real root. The pure imaginary roots are [14]
(6)λ=α(KP,KD)±ω(KP,KD)j,
where  α(KP,KD)=0 and ω(KP,KD)=KDg/mk.The other root can be achieved through computation [14]:
(7)α′(KP,KD)=Es02g(R′KP-1),
where E=[R′2(1+ω2)+KP2]-1.In practice,KP≫R′; then α′(KP,KD)<0.So there is a HOPF bifurcation in EMS whenc3,1=0. That is to say, the system has periodical solutions which means that there is a self-excited vibration at this point.We define the scope ofKP in which the system is stable as
(8)W=(KPL,KPR),
where KPL is the left pole of proportional feedback coefficient and KPR is the right pole, and
(9)KPL=R′mgk,(10)KPR=R′mgk+R′s02kKD.


## 4. Parameters Design of Adaptive Controller

The design of adaptive controller parameters is to choose the proper levitation air gap coefficientKP and the differential feedback coefficient KD, so that they satisfy the following conditions:(1)
KP∈W;(2)
KP and KPR keep a certain far distance.Condition(1) guarantees the stability. Condition (2) guarantees that the system has no HOPF bifurcation.In EMS, the mass levitated is not always fixed. According to (8), (9), and (10), W is a function of the mass levitated, and therefore it will slide with the mass varying. If KP keeps fixed, it will move beyond Wor approach to KPR while the mass levitated varying. The first condition broken will cause the system unstability, and the second condition broken will cause self-excited vibration. To avoid this, the proportional feedback coefficient should be adjusted with the varying of the mass levitated, confined in W, and kept far away from KPR.In order to provide enough stiffness and to avoid self-excited vibration, in this paper, the proportional feedback coefficient is placed in the middle of intervalW. That is, KP is determined by the following equation:
(11)KP=KPL+12(KPR-KPL)=R′mgk+R′s04kKD.So it is necessary to know the levitated mass continuously and automatically; then we can adjust the proportional feedback coefficient so that both condition(1) and condition (2) can be satisfied together. But in fact, the levitated mass cannot be measured directly easily. To overcome this problem, a system variable must be found such out that it can reflect the variance of the levitated mass and can be real time measured.As we know, when the system reaches the stable state, the current in the electromagnet, the electromagnet voltage, and the levitation mass satisfy the following equation:(12)I0=s0mgk,I0=U0R,(13)U=U0+LdIdt,
where L=2kN/s0 is the equivalent induction of electromagnet. Then, the following equation can be put forward:
(14)U=Rs0mgk+LdIdt.In (13), the acceleration of gravity g and the constant of electromagnet k are fixed. The levitation air gap s0, the resistance of electromagnet R, and the induction of electromagnet L are also constant. The voltage of electromagnet U and the current  in the electromagnet  I can be measured by sensor. dI/dt, that is, the velocity of current changing, can be calculated by using the current I.So the levitated mass can be calculated by the following equation:(15)m=(U-L(dI/dt)Rs0)2kg.That is to say, we can get the variance of the levitated mass from the electromagnet voltage and current and get the variance scope ofW from the variance of the levitated mass. Then, the proportional feedback coefficient can be adjusted timely so that it can satisfy both condition (1) and condition (2) together.


## 5. Simulation Testing

We have an EMS experiment equipment as shown in Figure2. Table 2 includes its parameters.Table 2
Parameters of Figure2.
k
m
(
kg
)
R
(
Ω
)
I
0
(
A
)
K
D
K
I
s
0 (m)
7
.
5599
×
10
-
4
750
2
34.3
200
−20
0.01Figure 2
EMS experiment equipment of NUDT.Then, according to (8), we have
(16)W=(6863297733).Simulation 1. Set KP0 in the middle of W; that is, KP0=83183; after the system runs for 40 seconds, the levitated mass decreases by 15% abruptly, and the step response of the system can be obtained as in Figure 3.Figure 3
Step response of Simulation 1.From Figure3 we know that due to the abrupt decrease of the levitated mass at time 40 s, the system, through a transitory regulating process, returns to a special state in which there is an evident self-excited vibration.According to (8), (9), and (10), due to the decrease of the levitated mass, W slides to 5356. And because of the fixed KP0, the distance between KP0 and KPR decreases to 5356, which results in KP0 approaching KPR. Then, the system produces a self-excited vibration as shown in Figure 3.Simulation 2.After the system runs for 40 seconds, the levitated mass decreases by 15% abruptly. The estimated levitation mass use (15) is shown in Figure 4 and the proportional feedback coefficient while the levitated mass varies according to HOPF bifurcation theory described in Section 4 is shown in Figure 5. The step response can be obtained as in Figure 6 in this case.Figure 4
Curve of the estimated levitation mass.Figure 5
Curve of the proportional feedback coefficient.Figure 6
Step response of Simulation 2.From Figure6, we know that after the levitated mass abruptly decreases by 15%, the system returns to the stable state, and there is no self-excited vibration.According to the former computation, when the levitated mass decreases by 15%,W slides to left 5356. According to (12), the electromagnet voltage decreases by 7.8%, and according to (8), KP slides to left 5356. So there is no change of the relative distance between KP and W. Therefore, the dynamic characteristics of the system have not been changed.


## 6. Conclusion

By analyzing the HOPF bifurcation of the nonlinear EMS type maglev system, a self-tuning adaptive control method is put forward, which uses an additional constraint as the adaptive rule.  This constraint makes the controller parameter always far away from the HOPF bifurcation points while the mass levitated varies. Proved by simulation, adaptive control method in this paper can make electromagnetic suspension system still stable after an abrupt changing of the mass levitated and avoid its self-excited vibration.The stability of the adaptive controller needs to be proven theoretically furthermore. And the adaptive relationship needs to be formalized furthermore.



---

*Source: 928719-2013-11-07.xml*