# Robust Output Feedback Stabilization of a Field-Sensed Magnetic Suspension System

Authors: Jen-Hsing Li; Juing-Shian Chiou
Journal: Mathematical Problems in Engineering
 (2013)

---

## Abstract

The magnetic suspension system (MSS) is very important in many engineering applications. This paper proposes the dynamic output feedback control of a field-sensed MSS (FSMSS). Subsequently, the mathematical model of the MSS is described by discrete-time systems. Ideally, the coefficients of a nominal polynomial can precisely determine the Schur stability. But in reality, the coefficients may contain uncertainties due to reasons such as computational errors. Therefore, there is a need to address the problem of robust stability for discrete-time systems. In this paper, the size of allowable perturbation in polynomial coefficient space was estimated for the output feedback control of the MSS. Theℓ∞-norm and a lower bound for the size of the Schur stability hypercube are provided in this paper.


---

## Body

## 1. Introduction

Magnetic suspension technique has numerous applications such as levitation of high speed trains, frictionless bearings, magnetically suspended wind tunnels, rocket-guiding projects, and vibration isolation tables in semiconductor manufacturing.The classroom demonstration device was first proposed by Wong [1] in 1986 with the phase-lead compensation and root locus analysis shown for the undergraduate course. In addition, a photoresistor was utilized to sense the position signal. In 1993, the sliding mode control was successfully extended to the magnetic suspension system (MSS) by Cho et al. [2]. The sliding mode control is used to stabilize the nonlinear systems in the global sense and considered robustness issues in modeling uncertainties and disturbances. In a sequel, the MSS was examined by μ-synthesis in [3, 4]. The unstructured modeling uncertainties [3] and the real parametric uncertainties [4] have been issued for the MSS device. Moreover, the μ-synthesis algorithm was utilized to design a robust controller to achieve robust performance. Oliveira et al. [5] proposed the robust control material applied to the MSS. The utilization of developed techniques in robust control such as H∞, μ analysis, and synthesis was issued. Oliveira et al. focused on education using the MSS to enhance the content of control courses. Likewise, previous studies by Oliveira et al. [6] proposed a digital control instructional material applicable to the MSS to reinforce control education. Shiakolas et al. [7, 8] proposed a hardware-in-the-loop environment in which the MSS was demonstrated. The emphasis was on digital control in [7] and neural network feed-forward control in [8]. MATLAB was used as the primary realization tool.Hurley and Wölfle [9] concentrated on electromagnet design of the MSS. Their study made use of a metallic sphere as the suspended object and showed that the electromagnet could be optimized for a given sphere. A set of design rules was proposed which enabled users to design an optimum electromagnet for a given size ball. Subsequently, Hurley et al. [10] concentrated on PWM (pulse width modulation) drive of the MSS. They analyzed the averaged circuit models and the small-signal circuit models of the PWM drive and confirmed that the PWM drive was efficient for the MSS. J. H. Li and T. H.  S. Li [11] concluded that thyristor drive was successful to stabilize the MSS earlier in 1995. Li [12] later extended the PWM drive to a field-sensed MSS (FSMSS) and guaranteed the stability of this FSMSS.The robust nonlinear compensation algorithm based on geometric feedback linearization techniques was developed in [13], and the MSS was exploited to guarantee the stability of the proposed scheme. The modeling and nonlinear control of the MSS was also investigated in [14]. Hajjaji and Ouladsine [14] focused on modeling and identification of the MSS, and afterwards the feedback linearization approach was employed to stabilize the MSS.Marsden [15] provided the FSMSS project kit, and Lundberg [16] used inexpensive experiments to teach analysis and design in undergraduate control courses. Referring to [15], the analog control of the FSMSS is a class project at MIT (Massachusetts Institute of Technology) on feedback systems. While in [17–19], the current driver is implemented by the PWM technique. In this paper, the FSMSS is extended to a digital control system, and the problem of robust stability on discrete-time systems is addressed.This paper is organized as follows. Section2 reviews the discrete model and the PD control of a MSS. Section 3 describes the output feedback control in detail. Section 4 addresses the robust stability analysis. The apparatus and experiments with the results are delineated in Section 5, and a brief conclusion is made in Section 6.


## 2. Review of Discrete Model and PD Control

The general mathematical model of an MSS is given as follows [1, 11, 12, 17–19]:
(1)md2xdt2=mg-Ci2x2,
where x is the distance between the electromagnet and the steel ball, m is the mass of the steel ball, g is the gravitational acceleration, i is the coil current, and C is the force constant. Using the Taylor series expansion and neglecting all higher-order terms, the piece-wise linearized equation is
(2)d2dt2Δx=-2Ci0mxo2Δi+2Ci02mxo3Δx.The equilibrium condition is(3)mg=Ci02xo2.Taking the Laplace transform of (2), the transfer function of MSS is
(4)G(s)=ΔX(s)ΔI(s)=-2Ci0/mxo2s2-2Ci02/mxo3,
where ΔX(s) is the Laplace transform of Δx and ΔI(s) is the Laplace transform of Δi. The discrete-time model of (4) can be calculated by the following formula [20]:
(5)G(z)=∑atpolesofG(α)[residuesofG(α)11-z-1eTα],
where T is the sampling period and can be expressed as the following equation:
(6)G(z)=ΔX(z)ΔI(z)=γ(zz-1/β-zz-β)=-zγ(β2-1)/β(z-β)(z-1/β).The parameters are given as follows:(7)β=eT2Ci02/mxo3>1,γ=C2mxo.The parameterγ~=γρ(β2-1)/β is a positive real number, where ρ is the linear factor of the position sensor circuit and its unit is volt/meter. R(z) is the reference input of the closed loop system. If the stable PD controller [17–19]
(8)GC(z)=Kz-1(z+φ)
is exploited as shown in Figure 1, the z-1 term is designed to cancel the z term of G(z) in (6). We introduced a zero z=-φ and assigned the gain K to move the roots into the unit circle. The resulting transfer function ΔX~(z)/R(z) is
(9)ΔX~(z)R(z)=(z+φ)Kγ~(z-β)(z-1/β)+(z+φ)Kγ~,
while the characteristic polynomial is
(10)Q(z)=z2+z[Kγ~-(β+1β)]+1+Kγ~φ.Figure 1
Block diagram of the PD control system.By the Jury stability criterion [20], the stable conditions for K and φ are
(11)(β-1)γρ(β+1)(1+φ)<K<(β+1)γρ(β-1)(1-φ),(12)|1+Kγ~φ|<1.Becauseγ~>0 and K>0, the inequality (12) implies
(13)-2Kγ~<φ<0.From inequalities (11) and (13), the stable PD controller is in certain of closed range. Hence we first assign the zero z=-φ then increase the gain K slowly to stabilize the MSS. If this controller is not satisfied, we can reassign another zero and gain until the controller is satisfied. This design procedure is the rule of thumb in related experiments.However, the current command of equilibrium point is the summation of the bias current (3) and the linear control current. The current command is given as follows:
(14)i(k)=io(k)+Δi(k)=io(k)+KΔx~(k)+KφΔx~(k-1),
where io(k) is the bias current of the equilibrium point, Δi(k) is the linear control current, Δx~(k) is the present position measurement perturbation, and Δx~(k-1) is the previous position measurement perturbation.


## 3. Output Feedback Control

From Figure1, the transfer function which is from ΔI(z) to ΔX~(z) is given as follows:
(15)ΔX~(z)ΔI(z)=zγ~z2-β~z+1,
where parameters are γ~=γρ(β2-1)/β and β~=β+β-1. The controllable canonical form of state space representation is (Figure 2)
(16)[x1(k+1)x2(k+1)]=[01-1β~][x1(k)x2(k)]+[01]Δi(k),Δx(k)=[0γ~][x1(k)x2(k)].Figure 2
Block diagram representation of (16).The dynamic output feedback control is(17)Δi(k)=K1Δx(k)+K2Δx(k-1)=K1γ~x2(k)+K2γ~x1(k).Hence, the closed loop system is(18)[x1(k+1)x2(k+1)]=[01-1+K2γ~β~+K1γ~][x1(k)x2(k)].The characteristic polynomial is(19)Q(z)=z2-(β~+K1γ~)z+(1-K2γ~).Comparing the characteristic polynomial of the PD type control (8) with (19), the controller parameters are
(20)K1=-K,K2=-Kφ.If the parameters are perturbed to beβ~+δβ~ and γ~+δγ~, then the characteristic polynomial is perturbed to be
(21)Q(z)+δQ(z)=z2-(β~+K1γ~)z+(1-K2γ~)-(δβ~+K1δγ~)z-K2δγ,~
where
(22)δQ(z)=-(δβ~+K1δγ~)z-K2δγ~=δ1z+δ0.


## 4. Robust Stability Analysis

Consider the nominal Schur polynomial of a discrete system [21]
(23)Q(z)=anzn+an-1zn-1+⋯+a1z+a0=an∏j=1n(z-zj),an>0
which is subject to the perturbation
(24)δQ(z)=δnzn+δn-1zn-1+⋯+δ1z+δ0.To measure this size, we used theℓ∞-norm; that is, for δQ⃑=[δ0,δ1,…,δn]T, ∥δQ⃑∥∞=maxj|δj|. The perturbed polynomial is P(z)=Q(z)+δQ(z). We introduced an (n-1)×(n-1) real matrix formed from the coefficients of Q(z) as follows:
(25)Qn-1C=Qn-1R-Qn-1L,
where
(26)Qn-1R=[a0a1a2⋯an-3an-20a0a1⋯an-4an-300a0⋯an-5an-4⋮⋮⋮⋮⋮⋮00000a0],Qn-1L=[a2a3⋯an-2an-1ana3a4⋯an-1an0a4a5⋯an00⋮⋮⋮⋮⋮⋮an0⋯000],detQn-1C=(-1)n(n-1)/2ann-1∏j,k=1,j<kn(1-zjzk),QinvC={qjk}=(Qn-1C)-1.Then we can obtain the following result, and the proof was provided in detail in [21].Theorem 1.
If∥δQ⃑∥∞<min(ℓ∞,ℓ1∞,ℓ-1∞), then the perturbed polynomial P(z) is Schur stable,  where
(27)ℓ∞=(2∑j,k=1n|qjk|)-1,ℓ1∞=Q(1)n+1,ℓ-1∞=(-1)nQ(-1)n+1.A simulation example [18] is given as follows. Let the discrete-time plant be G(z)=-z/(z-2)(z-0.5) and ρ=1. And let the controller be GC(z)=Kz-1(z-0.3). From (11), the closed loop system is stable for 0.714<K<3.46. Assume that the gain is set to be K=2 and the output feedback controller is K1=-2, K2=0.6. The controlled nominal polynomial is given as follows:
(28)Q(z)=z2-0.5z+0.4.From Theorem1, ℓ∞=0.3, ℓ1∞=0.3, and ℓ-1∞=0.6333. Hence the lower bound for the size of the maximal Schur stability hypercube is ∥δQ⃑∥∞<min(ℓ∞,ℓ1∞,ℓ-1∞)=0.3.If the controller is redesigned asGC(z)=Kz-1(z-0.7), the stable gain is 1.66<K<2.65. Assume that the gain is set to be K=2 and the output feedback controller is K1=-2, K2=1.4. The nominal controlled polynomial is given as follows:
(29)Q(z)=z2-0.5z+1.From Theorem1, ℓ∞=0.7, ℓ1∞=0.0333, and ℓ-1∞=0.3667. Hence the lower bound for the size of the maximal Schur stability hypercube is ∥δQ⃑∥∞<min(ℓ∞,ℓ1∞,ℓ-1∞)=0.0333.


## 5. Experiment and Result

The photograph of a FSMSS is shown in Figure3, while the magnet position sensor circuit is shown in Figure 4. The SS495A is a sensor which has a ratiometric output voltage, set by the supply voltage. It varies in proportion to the strength of the magnetic field. The inverting amplifier is adjusted, so that the magnet position signal varies from 0 Volt to 3 Volts. When the magnet is placed on the top of SS495A, the output voltage can be measured as 3 Volts. While a decreasing output voltage is manifested as the magnet is pulled away from the SS495A. So we can measure the position of the magnet by this circuit. From Figure 5, the bottom of the magnet is exactly xo=1 cm away from the sensor, then the measured voltage is 1.5 Volts, and the distance between the electromagnet and the top of the magnet is 4 cm. The bias current required to balance the gravitational force is 0.47 Amps.Figure 3
A photograph of FSMSS.Figure 4
Magnet position sensor circuit.Figure 5
The motion direction of the FSMSS.The current driver circuit is shown in Figure6. The isolated gate bipolar transistor (IGBT) is a power device which serves as a switch in power circuits. The TLP250 is a photocoupler which consists of a light emitting diode and an integrated photodetector. The TLP250 is suitable for the gate driving circuit of an IGBT. If the IGBT is on, the applied voltage is 12 Volts, and the coil is charged. If the IGBT is off, the applied voltage is zero, and the coil is discharged via the flywheel diode. Further, the average voltage is determined by duty cycle. In this PWM circuit, the modulation frequency is 18 kHz. The electromagnetic coil with low pass property can attenuate the high frequency harmonics making the coil current smooth. Subsequently, the rectangle with a hole (LA55-P) is a coil current sensor. The linear factor of the current sensor circuit is set to be 3 V/A by the noninverting amplifier.Figure 6
The PWM current driver circuit and current sensor circuit.The study made use of the eZdsp F2812 [22] embedded control circuit board. The pack of the eZdsp F2812 provides a Code Composer Studio (CCS) software that includes a text-editor, a C-complier, an assembler, a debugger, and a download tool. Users can implement control theories with the use of C-algorithms.The function block diagram is shown in Figure7. It includes two control loops in the form of the coil current control loop (inner loop) and the magnet position control loop (outer loop). For the two-loop control system, the inner loop is designed at first. The coil current is first designed by the proportional and integral (PI) controller [17–19]. Furthermore, Figure 8 shows the step response of the coil current from one Amp to two Amps. And base from Figure 8, the settling time of the coil current is about 20 ms. The dynamics of the coil current is faster than that of position; hence the gain can be approximated as unity. Afterwards, the outer control loop is designed using the discrete PD control to stabilize the FSMSS. The experiment is the set point stabilization by use of the PD control. The step response is shown in Figure 9. It can be observed that the FSMSS is stable.Figure 7
Robust output feedback control of a FSMSS.Figure 8
The step response of the coil current.Figure 9
The magnet position response.In response to the robust output feedback control, we designed additional experiments to demonstrate its effectiveness. Four kinds of suspended objects are used to examine its robustness. The first suspended object is shown in Figure10(a). The middle magnet cube was replaced by a 1 cm white foam tape. Subsequently, the weight of this suspended object is lighter than the three magnet cubes in Figure 3. The weight difference varies the system parameters β~ and γ~. Although changes have been made on the system parameters, the system is still able to maintain stability and guarantee robustness.Various kinds of suspended objects.
(a)
Two cubes and 1 cm foam tape
(b)
The plastic stick with magnets
(c)
The plastic stick with 2 extra magnets
(d)
The plastic stick with 3 extra magnetsFigures10(b)–10(d) exhibit three experiments examined using a plastic stick with magnets. The plastic stick is red in color with two small magnets embedded at both ends. In Figure 10(b), no extra magnets are added. However, Figures 10(c) and 10(d) had two and three extra magnets added at both ends, respectively. Subsequently, two magnets were added at the bottom, and a magnet is added at the top of the plastic stick. Results of the experiments showed that the three suspended objects are all balanced and stable in the air. These experiments demonstrate the robust stability of the proposed methodology in this paper.


## 6. Conclusion

This paper deals with the FSMSS, and a discrete robustness analysis of the FSMSS is issued. Also the PD control and dynamic output feedback control are both considered in this paper. Because the discrete model of the FSMSS is a second-order system, it is more convenient to perform design and analysis. The system parameters are always unknown or partially unknown, so the proposed discrete PD control is a good strategy to stabilize the MSS. Moreover, the proposed output feedback control can be applied to achieve better performance. The robustness analysis of this paper can measure the robust stability margin of the designed controller. The simulations and experiments are also provided in this paper.



---

*Source: 192431-2013-05-19.xml*