# A Bearingless Induction Motor Direct Torque Control and Suspension Force Control Based on Sliding Mode Variable Structure

Authors: Zebin Yang; Lin Chen; Xiaodong Sun; Weiming Sun; Dan Zhang
Journal: Mathematical Problems in Engineering
 (2017)

---

## Abstract

Aiming at the problems of the large torque ripple and unstable suspension performance in traditional direct torque control (DTC) for a bearingless induction motor (BIM), a new method of DTC is proposed based on sliding mode variable structure (SMVS). The sliding mode switching surface of the torque and flux linkage controller are constructed by torque error and flux error, and the exponential reaching law is used to design the SMVS direct torque controller. On the basis of the radial suspension force mathematical model of the BIM, a radial suspension force closed-loop control method is proposed by utilizing the inverse system theory and SMVS. The simulation models of traditional DTC and the new DTC method based on SMVS of the BIM are set up in the MATLAB/Simulink toolbox. On this basis, the experiments are carried out. Simulation and experiment results showed that the stable suspension operation of the BIM can be achieved with small torque ripple and flux ripple. Besides, the dynamic response and suspension performance of the motor are improved by the proposed method.


---

## Body

## 1. Introduction

Bearingless motors inherited the characteristics of the traditional magnetic bearing motors, such as being without lubrication, having no wear, and having no mechanical noise. Compared with the traditional magnetic bearing motor, the space utilization rate, electromagnetic efficiency, and other aspects have been improved. It is of great significance to the development of the field of the biochemical medicine and industrial field such as the heart blood pump, the turbo molecular pump, and the high speed flywheel energy storage system [1, 2]. The levitation windings of bearingless motors, breaking the balance of the original state of the air gap magnetic field, changing the distribution of the air gap magnetic field, and producing the Maxwell force which have an effect on the rotor were to realize the stable suspension. Many scholars at home and abroad have studied the bearingless permanent magnet motor. However, the bearingless induction motor (BIM) has the advantages of low loss, low cost, simple structure, and withstanding high temperature, and it is easy to realize the open-loop speed control by the voltage source inverter, which people are increasingly concerned with [3].The BIM is a strongly coupled nonlinear system. In order to realize the decoupling between the radial suspension force and the torque, the rotor flux oriented vector control is introduced into the BIM in [4]. However, the use of the rotor flux instead of the air gap flux in this method has bad influence on the rotor stable suspension. Reference [5] overcame the limitation of rotor flux orientation, introduced the air gap flux oriented of the BIM, and deduced the suspension force expression of the air gap flux oriented control, but the control algorithm is much more complicated and highly nonlinear, and the torque response will be instable, which makes the application of this method limited.Direct torque control (DTC) is another way of AC speed regulation after field oriented vector control. The method is simple and does not need the complex vector transform and current control, which improves the dynamic response of the system. From the retrieved papers, there are not many literatures researching on the DTC of bearingless motor. The DTC is introduced into the bearingless permanent magnet synchronous motor and the bearingless permanent magnet slice motor in [6, 7]. The DTC is applied to the BIM in [8, 9] and the DTC based on space vector pulse width modulation is proposed. The experimental results showed that the method can realize the stable suspension of the motor. However, the traditional DTC has the problems of high torque ripple and nonfixed switching frequency. Because the stator flux estimation directly affects the performance of DTC, a large number of papers have been studied on the stator flux estimation [10–12]. Some papers committed to improving the structure of the controller [13–16], but these controllers have increased the complexity of the system, which make the control system of the BIM more complicated and influence the motor stable suspension operation.In order to reduce the torque ripple and improve the performance of the traditional DTC of the BIM, sliding mode variable structure (SMVS) control is introduced to the DTC system. Two hysteresis controllers in the traditional DTC are replaced by the sliding mode controller. The switching table is no longer used to select the space voltage vector; instead, the voltage vectors of the two-phase static coordinate system are output by the sliding mode controller, which is constructed according to the torque error and flux linkage error. The choice of inverter switching value based on space voltage pulse width modulation (SVPWM) reduces the torque ripple fundamentally.The BIM control system includes torque control and suspension force control. For the part of suspension force control, the traditional method is by detecting the radial displacement of the rotor, and the given value of suspension force is obtained by the PID modulation, which can realize the closed-loop control of the suspension force [17]. The method is relevant to air gap flux identification, and its dynamic response and the anti-interference performance are not good. To overcome this problem, the SMVS suspension force control method based on inverse system is proposed. On the basis of the radial suspension force mathematical model of the BIM, the reversibility of the inverse system is analyzed; the SMVS control is introduced into the inverse system and realized suspension force control. The BIM DTC and suspension force control system based on SMVS are built in MATLAB/Simulink toolbox. The dynamic performance of the torque response and the rotor suspension performance are compared and analyzed for two cases that are based on traditional DTC and SMVS-DTC, respectively. On this basis, the experiments were carried out. Simulation and experiment results showed that the torque ripple has greatly reduced and the performance of the rotor suspension is improved.


## 2. Sliding Mode Variable Structure Direct Torque Control of the BIM

### 2.1. Principle of SMVS

SMVS control is a special kind of nonlinear control method, which is a kind of control strategy of variable structure control system. The difference between this control strategy and the traditional control strategy is that the control is not continuous. It is a kind of switching characteristic that makes the system structure change with time. The control strategy makes the system pass through the prescribed state locus back and forth in high frequency and small amplitude. Because the state locus can be designed, and not related to the system parameters and outside disturbance, it has a strong robustness [18].Sliding mode variable structure control has three basic problems: firstly, the sliding mode should exist; secondly, the reach-ability condition should be satisfied; finally, the stability of the sliding mode should be analyzed [19–21].


### 2.2. Design of SMVS-DTC Controller

SMVS controller is used to replace the traditional hysteresis comparators. The torque winding reference voltage vector should be calculated by torque error and flux error, and thus the controller includes flux controller and torque controller.The mathematical model of the BIM in the two-phase stationaryα-βcoordinate system and with the torque winding stator current and stator flux as state variables is as follows:(1)ddti1sα=-1σR1sL1s+R1rL1ri1sα-ωri1sβ+R1rσL1sL1rψ1sα+ωrσL1sψ1sβ+1σL1su1sα,ddti1sβ=-1σR1sL1s+R1rL1ri1sβ+ωri1sα+R1rσL1sL1rψ1sβ-ωrσL1sψ1sα+1σL1su1sβ,ddtψ1sα=u1sα-R1si1sα,ddtψ1sβ=u1sβ-R1si1sβ,where σ=1-L1m2/L1sL1r; i1sα and i1sβ, respectively, represent two components of the torque winding stator current and rotor current in coordinate system. ψ1rα and ψ1rβ, respectively, represent two components of the torque winding stator flux and rotor flux in coordinate system. L1s and L1r are stator and rotor self-inductance of torque winding. R1s and R1r are stator and rotor resistance of torque winding; L1m is mutual inductance between stator and rotor of torque winding. u1sα and u1sβ, respectively, represent two components of the torque winding stator voltage and rotor voltage in coordinate system. ωr is the rotor angular velocity.The electromagnetic torque equation of the BIM is as follows:(2)Te=32P1ψ1sαi1sβ-ψ1sβi1sα,where P1 is the number of the pole-pairs of torque winding.The voltage-current model is used to estimate the stator flux of the torque winding in the two-phase stationary coordinate system (α-β).(3)ψ1sα=∫u1sα-R1si1sαdt,ψ1sβ=∫u1sβ-R1si1sβdt.In order to track the desired trajectory of the torque and flux, the sliding mode switching surface is selected as follows:(4)S1=Te∗-Te,S2=ψ1s2∗-ψ1s2,S=S1S2T,where Te∗ and ψ1s2∗ are the given values of torque and the square of flux. Te and ψ1s2 are the calculated values of torque and the square of flux.To make the system have a good dynamic performance, the reaching law is selected as the exponential reaching law:(5)S˙=-εsgn⁡S-KS,where ε and K are positive constants. When the value of ε is very small and the value of K is great, the speed of the reaching law away from the switching surface is fast and near the switching surface is slow. This can effectively reduce the chattering and shorten the transition time [22].Calculating the derivation ofS1 and S2,(6)S˙1=-32P11σR1s/L1s+R1r/L1rψ1sβi1sα-ψ1sαi1sβ+ωrψ1sαi1sα+ψ1sβi1sβ-ωrσL1sψ1s2+i1sβ-ψ1sβσL1su1sα-i1sα-ψ1sβσL1su1sβ=C1+D1u,where ψ1s2=ψ1sα2+ψ1sβ2.(7)S˙2=-2ψ1sαu1sα-2ψ1sβu1sβ+2Rsψ1sαi1sα+ψ1sβi1sβ=C2+D2u,in whichC1,D1,C2,D2, anduare as follows:(8)C1=-32P11σR1sL1s+R1rL1rψ1sβi1sα-ψ1sαi1sβ+ωrψ1sαi1sα+ψ1sβi1sβ-ωrσL1sψ1s2,D1=-32P1i1sβ-ψ1sβσL1s32P1i1sα-ψ1sασL1s,C2=2Rsψ1sαi1sα+ψ1sβi1sβ,D2=-2ψ1sα-2ψ1sβ,u=u1sαu1sβT.Combining (5) and (6) one can obtain(9)S˙=S˙1S˙2=C1C2+D1D2u1sαu1sβ.Substituting (4) into (7) one can obtain the controller formula:(10)u=-D-1C1+ε1sgnS1+K1S1C2+ε2sgnS2+K2S2and D-1=d1d2d3d4, where(11)d1=-2ψ1sβ3p1ψ1sβi1sβ+ψ1sαi1sα-3p1ψ1sβ2+ψ1sα2/σL1s;d2=-1.5p1i1sα-ψ1sα/σL1s3p1ψ1sβi1sβ+ψ1sαi1sα-3p1ψ1sβ2+ψ1sα2/σL1s;d3=2ψ1sα3p1ψ1sβi1sβ+ψ1sαi1sα-3p1ψ1sβ2+ψ1sα2/σL1s;d4=-1.5p1i1sβ-ψ1sβ/σL1s3p1ψ1sβi1sβ+ψ1sαi1sα-3p1ψ1sβ2+ψ1sα2/σL1s.


### 2.3. Analyses of the Reach-Ability Condition of the Sliding Mode

If the sliding mode exists, the moving point beyond the switching surface will reach the switching surface in finite time and the sliding mode motion is stable. The following reach-ability condition is satisfied:(12)limS→0+S˙<0,limS→0-S˙>0.The global reach-condition isSS˙<0. From the exponential reaching law, it can be obtained that(13)SS˙=S-εsgnS-KS=-εS-KS2.Becauseε and K are positive constant, SS˙<0. The sliding mode satisfies the reach-ability condition.


### 2.4. Analyses of the Stability of SMVS-DTC Controller

Selecting Lyapunov function:V=SST/2, the derivation of it can be calculated:(14)V˙=S1S2-ε1sgnS1-K1S1-ε2sgnS2-K2S2=-S1ε1sgnS1+K1S1-S2ε2sgnS2+K2S2=-ε1S1-K1S12-ε2S2-K2S22,in whichε1, K1, ε2, and K2 are positive constant, so it can be obtained that(15)V˙=-ε1S1-K1S12-ε2S2-K2S22<0.It can be concluded that this SMVS-DTC controller satisfies the stability conditions.


## 2.1. Principle of SMVS

SMVS control is a special kind of nonlinear control method, which is a kind of control strategy of variable structure control system. The difference between this control strategy and the traditional control strategy is that the control is not continuous. It is a kind of switching characteristic that makes the system structure change with time. The control strategy makes the system pass through the prescribed state locus back and forth in high frequency and small amplitude. Because the state locus can be designed, and not related to the system parameters and outside disturbance, it has a strong robustness [18].Sliding mode variable structure control has three basic problems: firstly, the sliding mode should exist; secondly, the reach-ability condition should be satisfied; finally, the stability of the sliding mode should be analyzed [19–21].


## 2.2. Design of SMVS-DTC Controller

SMVS controller is used to replace the traditional hysteresis comparators. The torque winding reference voltage vector should be calculated by torque error and flux error, and thus the controller includes flux controller and torque controller.The mathematical model of the BIM in the two-phase stationaryα-βcoordinate system and with the torque winding stator current and stator flux as state variables is as follows:(1)ddti1sα=-1σR1sL1s+R1rL1ri1sα-ωri1sβ+R1rσL1sL1rψ1sα+ωrσL1sψ1sβ+1σL1su1sα,ddti1sβ=-1σR1sL1s+R1rL1ri1sβ+ωri1sα+R1rσL1sL1rψ1sβ-ωrσL1sψ1sα+1σL1su1sβ,ddtψ1sα=u1sα-R1si1sα,ddtψ1sβ=u1sβ-R1si1sβ,where σ=1-L1m2/L1sL1r; i1sα and i1sβ, respectively, represent two components of the torque winding stator current and rotor current in coordinate system. ψ1rα and ψ1rβ, respectively, represent two components of the torque winding stator flux and rotor flux in coordinate system. L1s and L1r are stator and rotor self-inductance of torque winding. R1s and R1r are stator and rotor resistance of torque winding; L1m is mutual inductance between stator and rotor of torque winding. u1sα and u1sβ, respectively, represent two components of the torque winding stator voltage and rotor voltage in coordinate system. ωr is the rotor angular velocity.The electromagnetic torque equation of the BIM is as follows:(2)Te=32P1ψ1sαi1sβ-ψ1sβi1sα,where P1 is the number of the pole-pairs of torque winding.The voltage-current model is used to estimate the stator flux of the torque winding in the two-phase stationary coordinate system (α-β).(3)ψ1sα=∫u1sα-R1si1sαdt,ψ1sβ=∫u1sβ-R1si1sβdt.In order to track the desired trajectory of the torque and flux, the sliding mode switching surface is selected as follows:(4)S1=Te∗-Te,S2=ψ1s2∗-ψ1s2,S=S1S2T,where Te∗ and ψ1s2∗ are the given values of torque and the square of flux. Te and ψ1s2 are the calculated values of torque and the square of flux.To make the system have a good dynamic performance, the reaching law is selected as the exponential reaching law:(5)S˙=-εsgn⁡S-KS,where ε and K are positive constants. When the value of ε is very small and the value of K is great, the speed of the reaching law away from the switching surface is fast and near the switching surface is slow. This can effectively reduce the chattering and shorten the transition time [22].Calculating the derivation ofS1 and S2,(6)S˙1=-32P11σR1s/L1s+R1r/L1rψ1sβi1sα-ψ1sαi1sβ+ωrψ1sαi1sα+ψ1sβi1sβ-ωrσL1sψ1s2+i1sβ-ψ1sβσL1su1sα-i1sα-ψ1sβσL1su1sβ=C1+D1u,where ψ1s2=ψ1sα2+ψ1sβ2.(7)S˙2=-2ψ1sαu1sα-2ψ1sβu1sβ+2Rsψ1sαi1sα+ψ1sβi1sβ=C2+D2u,in whichC1,D1,C2,D2, anduare as follows:(8)C1=-32P11σR1sL1s+R1rL1rψ1sβi1sα-ψ1sαi1sβ+ωrψ1sαi1sα+ψ1sβi1sβ-ωrσL1sψ1s2,D1=-32P1i1sβ-ψ1sβσL1s32P1i1sα-ψ1sασL1s,C2=2Rsψ1sαi1sα+ψ1sβi1sβ,D2=-2ψ1sα-2ψ1sβ,u=u1sαu1sβT.Combining (5) and (6) one can obtain(9)S˙=S˙1S˙2=C1C2+D1D2u1sαu1sβ.Substituting (4) into (7) one can obtain the controller formula:(10)u=-D-1C1+ε1sgnS1+K1S1C2+ε2sgnS2+K2S2and D-1=d1d2d3d4, where(11)d1=-2ψ1sβ3p1ψ1sβi1sβ+ψ1sαi1sα-3p1ψ1sβ2+ψ1sα2/σL1s;d2=-1.5p1i1sα-ψ1sα/σL1s3p1ψ1sβi1sβ+ψ1sαi1sα-3p1ψ1sβ2+ψ1sα2/σL1s;d3=2ψ1sα3p1ψ1sβi1sβ+ψ1sαi1sα-3p1ψ1sβ2+ψ1sα2/σL1s;d4=-1.5p1i1sβ-ψ1sβ/σL1s3p1ψ1sβi1sβ+ψ1sαi1sα-3p1ψ1sβ2+ψ1sα2/σL1s.


## 2.3. Analyses of the Reach-Ability Condition of the Sliding Mode

If the sliding mode exists, the moving point beyond the switching surface will reach the switching surface in finite time and the sliding mode motion is stable. The following reach-ability condition is satisfied:(12)limS→0+S˙<0,limS→0-S˙>0.The global reach-condition isSS˙<0. From the exponential reaching law, it can be obtained that(13)SS˙=S-εsgnS-KS=-εS-KS2.Becauseε and K are positive constant, SS˙<0. The sliding mode satisfies the reach-ability condition.


## 2.4. Analyses of the Stability of SMVS-DTC Controller

Selecting Lyapunov function:V=SST/2, the derivation of it can be calculated:(14)V˙=S1S2-ε1sgnS1-K1S1-ε2sgnS2-K2S2=-S1ε1sgnS1+K1S1-S2ε2sgnS2+K2S2=-ε1S1-K1S12-ε2S2-K2S22,in whichε1, K1, ε2, and K2 are positive constant, so it can be obtained that(15)V˙=-ε1S1-K1S12-ε2S2-K2S22<0.It can be concluded that this SMVS-DTC controller satisfies the stability conditions.


## 3. SMVS Suspension Force Control of the BIM

### 3.1. Radial Suspension Force Mathematical Model of the BIM

According to the relationship between the suspension force and the suspension winding stator current of the BIM in the two-phase stationary coordinate system (α-β) [23],(16)Fx=Kmi2sαψ1α+i2sβψ1β,Fy=Kmi2sβψ1α-i2sαψ1β,where Km=πP1P2L2m/18μ0lrW1W2. P1 and P2 are the numbers of pole-pairs of torque winding and suspension winding, respectively. ψ1α and ψ1β are the component of air gap flux in α-β axis, respectively. Fx and Fy are the suspension force of radial x and y. L2m is the mutual inductance between stator and rotor of suspension winding. μ0 is the permeability of vacuum. l is the effective length of the rotor. r is the outside diameter of rotor. W1 and W2 are the torque winding turns and suspension winding turns.When the rotor deviates from the center of the motor stator, it will cause the motor flux distribution. Then the Maxwell force is not zero, and the direction of its action and the direction of the rotor eccentricity are consistent, pointing to the direction of minimum air gap; it is called unbalanced magnetic force and is showed as(17)Fsx=ksx,Fsy=ksy,where ks=kπrlB12/2μ0g0 and Fsx and Fsy are the unbalanced magnetic force in the direction of x and y, respectively. x and y are the displacement in the direction of x and y, respectively, and g0 is the air gap length.The motion equations of the rotor are expressed as(18)mx¨=Fx-Fsx,my¨=Fy-Fsy-mg,where m is the mass of the rotor and g is the gravity acceleration. By (16)–(18), it can be obtained that(19)x¨=Kmi2sαψ1α+i2sβψ1β-ksxm,y¨=Kmi2sβψ1α-i2sαψ1β-ksy-mgm.


### 3.2. SMVS Radial Suspension Force Control of the BIM Based on Inverse System

The basic principle of inverse system method is to compensate the controlled object as a system with a linear transfer relationship by using the inverse system of the controlled object, then, to synthesize the system according to the linear system theory, and to realize the decoupling in nonlinear systems [24, 25].Theorem 1.
For a system ofp dimensional input, u(t)=(u1,u2,…,up)T, and, for q dimensional output, y(t)=(y1,y2,…,yq)T, with a set of initial states x(t0)=x0, and it can be expressed as the following state equation:(20)x˙=fx,u,y=hx,u,xt0=x0.The necessary and sufficient condition for the system to be reversible in the neighborhood of (x0,u0) has vector relative order α=(α1,α2,…,αq)T in this neighborhood.Select the state variables in (19).(21)X=x1x2x3x4T=xyx˙y˙T.The input variable is(22)U=u1u2T=i2sαi2sβT.The output variable is(23)Y=y1y2T=xyT.The state equation of the composite controlled object is obtained by (19) and(21):(24)x˙1=x3,x˙2=x4,x˙3=Kmmu1ψ1α+u2ψ1β-ksmx1,x˙4=Kmmu2ψ1α-u1ψ1β-ksmx2-g.Taking the derivative of the output function until each component inYq=(y1α1,…,yqαq) include input u explicit. (25)y˙1=x3,y¨1=Kmmu1ψ1α+u2ψ1β-ksmx1,y˙2=x4,y¨2=Kmmu2ψ1α-u1ψ1β-ksmx2-g.The Jacobi matrix of the inputu can be expressed as(26)A=∂y¨1∂u1∂y¨1∂u2∂y¨2∂u1∂y¨2∂u2=Kmmψ1αψ1β-ψ1βψ1α.BecauseDetA=Km(ψ1α2+ψ1β2)/m≠0, according to Theorem 1, the relative order of the system is α=(α1,a2)T=2,2T, α1+a2=4=n, so the system is reversible.AssumingY¨=y¨1y¨2T=ϕ1ϕ2T, (25) can be calculated as(27)u1=mKmψ1α2+ψ1β2ψ1αϕ1-ψ1βϕ2+ksmx1ψ1α-ksmx2ψ1β-gψ1β,u2=mKmψ1α2+ψ1β2ψ1βϕ1+ψ1αϕ2+ksmx1ψ1β+ksmx2ψ1α+gψ1α.In order to realize closed-loop displacement control and improve the performance of anti-interference, SMVS is introduced into the inverse system. The sliding mode switching surface of the sliding mode controller is defined as follows:(28)S3=c1x∗-x+x˙∗-x˙,S4=c2y∗-y+y˙∗-y˙,where x∗ and y∗ are the given displacement of x and y direction, x and y are the measured displacement, and c1 and c2 are constants.Exponential reaching law is adopted in this paper and the expression is as follows:(29)S˙3=-ε3sgnS3-K3S3,S˙4=-ε4sgnS4-K4S4.According to (14)-(15), the sliding mode controller is stable. It can be obtained by (28)-(29) that (30)x¨=ϕ1=c1x˙∗-x˙+x¨∗+ε3sgnS3+K3S3,y¨=ϕ1=c1y˙∗-y˙+y¨∗+ε4sgnS4+K4S4.Torque winding and suspension winding of the BIM are coupled by air gap flux of torque winding, and torque winding can be obtained by torque winding stator flux subtracting stator leakage inductance.(31)ψ1α=ψ1sα-L1sli1sα,ψ1β=ψ1sβ-L1sli1sβ.It can be known from (3) that the stator flux identification in the DTC has the pure integral part. The initial value of integration, the accumulated error, and other factors will make the flux identification inaccurate. Therefore, the traditional suspension force control method cannot meet the requirements of the stable suspension of the BIM. Because of the introduction of SMVS to the suspension force control, the robustness of the controller is improved, the dependence of the accuracy of stator flux estimation is reduced, and the performance of stator stable suspension is improved.Through the separate control of torque and suspension force, the SMVS-DTC system of the BIM can be constructed, and the structure diagram of the system is showed as Figure1.Figure 1
The structure diagram of the BIM SMVS-DTC system.The control block diagram of the BIM includes two parts: torque winding control part and suspension winding control part. In torque winding control part, the voltage and current of the stator are obtained by the coordinate transformation, and then the estimated value of the stator flux is calculated by the stator flux calculation block. The given value of the torque is obtained by the speed regulator, and the estimation value of torque is calculated by torque calculation block. The torque error signal is obtained by given value and the estimation value of torque, and then the voltage control signal of inverter is obtained by SMVS-DTC block. In suspension winding control part, through the SMVS block, the component of levitation force can be calculated by position error signal. The air gap flux is calculated by the stator flux of torque winding; through the inverse system block, the stator current of suspension winding can be calculated by the component of levitation force and the air gap flux, so as to realize radial stable suspension of the rotor.


## 3.1. Radial Suspension Force Mathematical Model of the BIM

According to the relationship between the suspension force and the suspension winding stator current of the BIM in the two-phase stationary coordinate system (α-β) [23],(16)Fx=Kmi2sαψ1α+i2sβψ1β,Fy=Kmi2sβψ1α-i2sαψ1β,where Km=πP1P2L2m/18μ0lrW1W2. P1 and P2 are the numbers of pole-pairs of torque winding and suspension winding, respectively. ψ1α and ψ1β are the component of air gap flux in α-β axis, respectively. Fx and Fy are the suspension force of radial x and y. L2m is the mutual inductance between stator and rotor of suspension winding. μ0 is the permeability of vacuum. l is the effective length of the rotor. r is the outside diameter of rotor. W1 and W2 are the torque winding turns and suspension winding turns.When the rotor deviates from the center of the motor stator, it will cause the motor flux distribution. Then the Maxwell force is not zero, and the direction of its action and the direction of the rotor eccentricity are consistent, pointing to the direction of minimum air gap; it is called unbalanced magnetic force and is showed as(17)Fsx=ksx,Fsy=ksy,where ks=kπrlB12/2μ0g0 and Fsx and Fsy are the unbalanced magnetic force in the direction of x and y, respectively. x and y are the displacement in the direction of x and y, respectively, and g0 is the air gap length.The motion equations of the rotor are expressed as(18)mx¨=Fx-Fsx,my¨=Fy-Fsy-mg,where m is the mass of the rotor and g is the gravity acceleration. By (16)–(18), it can be obtained that(19)x¨=Kmi2sαψ1α+i2sβψ1β-ksxm,y¨=Kmi2sβψ1α-i2sαψ1β-ksy-mgm.


## 3.2. SMVS Radial Suspension Force Control of the BIM Based on Inverse System

The basic principle of inverse system method is to compensate the controlled object as a system with a linear transfer relationship by using the inverse system of the controlled object, then, to synthesize the system according to the linear system theory, and to realize the decoupling in nonlinear systems [24, 25].Theorem 1.
For a system ofp dimensional input, u(t)=(u1,u2,…,up)T, and, for q dimensional output, y(t)=(y1,y2,…,yq)T, with a set of initial states x(t0)=x0, and it can be expressed as the following state equation:(20)x˙=fx,u,y=hx,u,xt0=x0.The necessary and sufficient condition for the system to be reversible in the neighborhood of (x0,u0) has vector relative order α=(α1,α2,…,αq)T in this neighborhood.Select the state variables in (19).(21)X=x1x2x3x4T=xyx˙y˙T.The input variable is(22)U=u1u2T=i2sαi2sβT.The output variable is(23)Y=y1y2T=xyT.The state equation of the composite controlled object is obtained by (19) and(21):(24)x˙1=x3,x˙2=x4,x˙3=Kmmu1ψ1α+u2ψ1β-ksmx1,x˙4=Kmmu2ψ1α-u1ψ1β-ksmx2-g.Taking the derivative of the output function until each component inYq=(y1α1,…,yqαq) include input u explicit. (25)y˙1=x3,y¨1=Kmmu1ψ1α+u2ψ1β-ksmx1,y˙2=x4,y¨2=Kmmu2ψ1α-u1ψ1β-ksmx2-g.The Jacobi matrix of the inputu can be expressed as(26)A=∂y¨1∂u1∂y¨1∂u2∂y¨2∂u1∂y¨2∂u2=Kmmψ1αψ1β-ψ1βψ1α.BecauseDetA=Km(ψ1α2+ψ1β2)/m≠0, according to Theorem 1, the relative order of the system is α=(α1,a2)T=2,2T, α1+a2=4=n, so the system is reversible.AssumingY¨=y¨1y¨2T=ϕ1ϕ2T, (25) can be calculated as(27)u1=mKmψ1α2+ψ1β2ψ1αϕ1-ψ1βϕ2+ksmx1ψ1α-ksmx2ψ1β-gψ1β,u2=mKmψ1α2+ψ1β2ψ1βϕ1+ψ1αϕ2+ksmx1ψ1β+ksmx2ψ1α+gψ1α.In order to realize closed-loop displacement control and improve the performance of anti-interference, SMVS is introduced into the inverse system. The sliding mode switching surface of the sliding mode controller is defined as follows:(28)S3=c1x∗-x+x˙∗-x˙,S4=c2y∗-y+y˙∗-y˙,where x∗ and y∗ are the given displacement of x and y direction, x and y are the measured displacement, and c1 and c2 are constants.Exponential reaching law is adopted in this paper and the expression is as follows:(29)S˙3=-ε3sgnS3-K3S3,S˙4=-ε4sgnS4-K4S4.According to (14)-(15), the sliding mode controller is stable. It can be obtained by (28)-(29) that (30)x¨=ϕ1=c1x˙∗-x˙+x¨∗+ε3sgnS3+K3S3,y¨=ϕ1=c1y˙∗-y˙+y¨∗+ε4sgnS4+K4S4.Torque winding and suspension winding of the BIM are coupled by air gap flux of torque winding, and torque winding can be obtained by torque winding stator flux subtracting stator leakage inductance.(31)ψ1α=ψ1sα-L1sli1sα,ψ1β=ψ1sβ-L1sli1sβ.It can be known from (3) that the stator flux identification in the DTC has the pure integral part. The initial value of integration, the accumulated error, and other factors will make the flux identification inaccurate. Therefore, the traditional suspension force control method cannot meet the requirements of the stable suspension of the BIM. Because of the introduction of SMVS to the suspension force control, the robustness of the controller is improved, the dependence of the accuracy of stator flux estimation is reduced, and the performance of stator stable suspension is improved.Through the separate control of torque and suspension force, the SMVS-DTC system of the BIM can be constructed, and the structure diagram of the system is showed as Figure1.Figure 1
The structure diagram of the BIM SMVS-DTC system.The control block diagram of the BIM includes two parts: torque winding control part and suspension winding control part. In torque winding control part, the voltage and current of the stator are obtained by the coordinate transformation, and then the estimated value of the stator flux is calculated by the stator flux calculation block. The given value of the torque is obtained by the speed regulator, and the estimation value of torque is calculated by torque calculation block. The torque error signal is obtained by given value and the estimation value of torque, and then the voltage control signal of inverter is obtained by SMVS-DTC block. In suspension winding control part, through the SMVS block, the component of levitation force can be calculated by position error signal. The air gap flux is calculated by the stator flux of torque winding; through the inverse system block, the stator current of suspension winding can be calculated by the component of levitation force and the air gap flux, so as to realize radial stable suspension of the rotor.


## 4. Simulation and Experimental Study of the BIM SMVS-DTC and Suspension Force Control System

### 4.1. Results and Analysis of Simulation

In order to verify the feasibility of the control method, the simulation model of the BIM SMVS-DTC and suspension force control system is established. The simulation results of SMVS-DTC control method proposed and the traditional DTC control method are compared and analyzed. The simulation parameters are in Table1.Table 1
Parameters of the BIM.
Motor parameters
Torque winding
Suspension winding
P (W)
1500
500
n (r/min)
6000
6000
R
s (Ω)
2.01
1.03
R
r (Ω)
11.48
0.075
L
m (H)
0.15856
0.00932
L
1
s
l
(
H
)
0.00454
0.00267
L
1
r
l
(
H
)
0.00922
0.00542
J (kg⋅m2)
0.00769
0.00769
m (kg)
2.85
2.85
p
1
2The simulation parameters are as follows:ε1=ε2=0.03, K1=K2=1500, ε3=ε4=0.01, and K3=K4=3000. The given speed is 6000 r/min, and the simulation time is 0.6 seconds. The load torque is added to 2 N·m at t = 0.3 s. To verify the superiority of the proposed method, SMVS-DTC with SMVS suspension force control method and the traditional DTC with traditional suspension force control method are compared and analyzed by simulation. The simulation results are in Figures 2, 3, and 4.Figure 2
Trajectory of stator flux.
(a)
DTC
(b)
SMVS-DTCFigure 3
Torque response.
(a)
DTC
(b)
SMVS-DTCFigure 4
Speed response.The stator flux trajectory, torque response, and speed response are showed in Figures2, 3, and 4, respectively. It can be seen from Figure 2(a) that the amplitude of stator flux oscillations is large and the dynamic response is not good by traditional DTC method. From Figure 2(b) we can see that the amplitude of stator flux oscillations is smaller by SMVS-DTC method. Compared with Figures 3(a) and 3(b), we can see that the torque ripple of SMVS-DTC is reduced greatly than traditional DTC. And the dynamic performance of the motor is also improved when the motor is started. When t = 0.3 s the load torque is added. The transition time of the proposed method is reduced. It can be found in Figure 3(b) that there is chattering at steady state due to the sliding mode variable structure control and it can not be avoided, but the chattering is small, which does not affect the stable operation of the motor. From Figure 4, it can be found that the speed rise time has been reduced, and the dynamic response has been improved when using SMVS-DTC method compared to traditional DTC method. In addition, when the load is added, the speed fluctuation is obviously reduced when using SMVS-DTC method.The simulation results of the rotor motion trajectory and rotor radial displacement by traditional suspension force control method and the proposed method are showed in Figures5 and 6, respectively. We can see that the proposed method can effectively reduce the radial displacement of the rotor and the displacement fluctuation time, which can realize stable suspension operation of the BIM. In order to verify the performance of antidisturbance of the suspension force control method, the given radial displacement was set to 0.05 mm when t = 0.3, and the simulation waveforms are shown  in Figure 7.Figure 5
Rotor trajectory.
(a)
DTC
(b)
SMVS-DTCFigure 6
Rotor radial displacement.
(a)
DTC
(b)
SMVS-DTCFigure 7
Rotor radial displacement.As can be seen from the figure, the radial displacement fluctuation is only 0.015 mm, when the given displacement suddenly changed, and the fluctuation time is less than 0.05 seconds, so it has good performance of robustness.


### 4.2. Results and Analysis of Experiment

In order to verify the effect of the proposed method, experimental research on the prototype of the BIM has been carried out. We analyzed and compared this method with the traditional DTC method. The TI company’s 32-bit fixed point DSP chip F2812 is used as the control chip of the control system. The EVA and EVB on this chip can produce six complementary PWM waveforms; therefore, the torque winding and suspension winding of the BIM can be independently controlled. The experimental parameters are similar to the simulation parameters, and rotor speed is set to 2000 r/min, started with no-load, and then the load torque of 2 N·m is added to the motor. The torque responses of these two kinds of DTC methods are compared with the experimental results. The experimental block diagram of the control system is shown in Figure8 and the experimental results are shown in Figures 9 and 10.Figure 8
The experiment block diagram of the control system.Figure 9
Rotor radial displacement.Figure 10
Speed and torque response.
(a)
DTC
(b)
SMVS-DTCThe experimental block diagram of the control system includes two parts: torque winding control part and suspension winding control part. In torque winding control part, the speed signal is obtained by the photoelectric encoder and the DSP (Digital Signal Processor) as the control chip of the control system to achieve the SMVS-DTC control. The sample of busbar current and voltage was achieved by the current and voltage sensor. The PWM signal output by DSP was put into IPM (Intelligent Power Module). In suspension winding control part, the radial displacement signal of the suspension winding is obtained by the position sensor, and this signal was put into DSP to realize the closed-loop control of radial displacement. The power drive part of the suspension winding is the same as the torque winding.Figure9 is the displacement response of the rotor with SMVS suspension force control method based on inverse system. Figure 10 is the speed and torque response of traditional DTC and SMVS-DTC methods. Figure 9 shows that the rotor displacement in x and y direction is small with the suspension force control method, and it can realize the motor stable suspension operation. We can see from Figure 10(a) that the torque ripple is relatively large when the motor operates with no load and with load, and the torque ripple is large and the dynamic performance is not good when the load is suddenly added; at the same time, the rotor speed is significantly fluctuated. In addition, the BIM is vibrated when the motor starts, accompanied by a sharp noise. It can be known from Figures 10(a) and 10(b) that the torque ripple can be reduced effectively by the proposed method, the torque transition process is faster, and the speed is not fluctuated obviously when the load is suddenly added. The operation noise of experimental prototype is reduced and the starting process is stable with no violent vibration.


## 4.1. Results and Analysis of Simulation

In order to verify the feasibility of the control method, the simulation model of the BIM SMVS-DTC and suspension force control system is established. The simulation results of SMVS-DTC control method proposed and the traditional DTC control method are compared and analyzed. The simulation parameters are in Table1.Table 1
Parameters of the BIM.
Motor parameters
Torque winding
Suspension winding
P (W)
1500
500
n (r/min)
6000
6000
R
s (Ω)
2.01
1.03
R
r (Ω)
11.48
0.075
L
m (H)
0.15856
0.00932
L
1
s
l
(
H
)
0.00454
0.00267
L
1
r
l
(
H
)
0.00922
0.00542
J (kg⋅m2)
0.00769
0.00769
m (kg)
2.85
2.85
p
1
2The simulation parameters are as follows:ε1=ε2=0.03, K1=K2=1500, ε3=ε4=0.01, and K3=K4=3000. The given speed is 6000 r/min, and the simulation time is 0.6 seconds. The load torque is added to 2 N·m at t = 0.3 s. To verify the superiority of the proposed method, SMVS-DTC with SMVS suspension force control method and the traditional DTC with traditional suspension force control method are compared and analyzed by simulation. The simulation results are in Figures 2, 3, and 4.Figure 2
Trajectory of stator flux.
(a)
DTC
(b)
SMVS-DTCFigure 3
Torque response.
(a)
DTC
(b)
SMVS-DTCFigure 4
Speed response.The stator flux trajectory, torque response, and speed response are showed in Figures2, 3, and 4, respectively. It can be seen from Figure 2(a) that the amplitude of stator flux oscillations is large and the dynamic response is not good by traditional DTC method. From Figure 2(b) we can see that the amplitude of stator flux oscillations is smaller by SMVS-DTC method. Compared with Figures 3(a) and 3(b), we can see that the torque ripple of SMVS-DTC is reduced greatly than traditional DTC. And the dynamic performance of the motor is also improved when the motor is started. When t = 0.3 s the load torque is added. The transition time of the proposed method is reduced. It can be found in Figure 3(b) that there is chattering at steady state due to the sliding mode variable structure control and it can not be avoided, but the chattering is small, which does not affect the stable operation of the motor. From Figure 4, it can be found that the speed rise time has been reduced, and the dynamic response has been improved when using SMVS-DTC method compared to traditional DTC method. In addition, when the load is added, the speed fluctuation is obviously reduced when using SMVS-DTC method.The simulation results of the rotor motion trajectory and rotor radial displacement by traditional suspension force control method and the proposed method are showed in Figures5 and 6, respectively. We can see that the proposed method can effectively reduce the radial displacement of the rotor and the displacement fluctuation time, which can realize stable suspension operation of the BIM. In order to verify the performance of antidisturbance of the suspension force control method, the given radial displacement was set to 0.05 mm when t = 0.3, and the simulation waveforms are shown  in Figure 7.Figure 5
Rotor trajectory.
(a)
DTC
(b)
SMVS-DTCFigure 6
Rotor radial displacement.
(a)
DTC
(b)
SMVS-DTCFigure 7
Rotor radial displacement.As can be seen from the figure, the radial displacement fluctuation is only 0.015 mm, when the given displacement suddenly changed, and the fluctuation time is less than 0.05 seconds, so it has good performance of robustness.


## 4.2. Results and Analysis of Experiment

In order to verify the effect of the proposed method, experimental research on the prototype of the BIM has been carried out. We analyzed and compared this method with the traditional DTC method. The TI company’s 32-bit fixed point DSP chip F2812 is used as the control chip of the control system. The EVA and EVB on this chip can produce six complementary PWM waveforms; therefore, the torque winding and suspension winding of the BIM can be independently controlled. The experimental parameters are similar to the simulation parameters, and rotor speed is set to 2000 r/min, started with no-load, and then the load torque of 2 N·m is added to the motor. The torque responses of these two kinds of DTC methods are compared with the experimental results. The experimental block diagram of the control system is shown in Figure8 and the experimental results are shown in Figures 9 and 10.Figure 8
The experiment block diagram of the control system.Figure 9
Rotor radial displacement.Figure 10
Speed and torque response.
(a)
DTC
(b)
SMVS-DTCThe experimental block diagram of the control system includes two parts: torque winding control part and suspension winding control part. In torque winding control part, the speed signal is obtained by the photoelectric encoder and the DSP (Digital Signal Processor) as the control chip of the control system to achieve the SMVS-DTC control. The sample of busbar current and voltage was achieved by the current and voltage sensor. The PWM signal output by DSP was put into IPM (Intelligent Power Module). In suspension winding control part, the radial displacement signal of the suspension winding is obtained by the position sensor, and this signal was put into DSP to realize the closed-loop control of radial displacement. The power drive part of the suspension winding is the same as the torque winding.Figure9 is the displacement response of the rotor with SMVS suspension force control method based on inverse system. Figure 10 is the speed and torque response of traditional DTC and SMVS-DTC methods. Figure 9 shows that the rotor displacement in x and y direction is small with the suspension force control method, and it can realize the motor stable suspension operation. We can see from Figure 10(a) that the torque ripple is relatively large when the motor operates with no load and with load, and the torque ripple is large and the dynamic performance is not good when the load is suddenly added; at the same time, the rotor speed is significantly fluctuated. In addition, the BIM is vibrated when the motor starts, accompanied by a sharp noise. It can be known from Figures 10(a) and 10(b) that the torque ripple can be reduced effectively by the proposed method, the torque transition process is faster, and the speed is not fluctuated obviously when the load is suddenly added. The operation noise of experimental prototype is reduced and the starting process is stable with no violent vibration.


## 5. Conclusion

Aimed at the problems of the large torque ripple and unstable suspension performance in traditional DTC and suspension force control of the BIM, the traditional method is improved in this paper and the DTC and suspension force control method based on SMVS are proposed. The performance of torque control and suspension force control of the BIM is better with the proposed method. From the simulation and experiment results of the control system of the BIM based on the traditional DTC and SMVS-DTC, we can draw the following conclusions:(1)
The torque ripple can be reduced effectively and the dynamic performance and the anti-interference ability are improved by SMVS-DTC method.(2)
The SMVS suspension force control method not only reduced the rotor radial displacement fluctuation and fluctuation time but also realized the motor stable suspension operation.



---

*Source: 2409179-2017-08-01.xml*