# Sliding Mode Control for Bearingless Induction Motor Based on a Novel Load Torque Observer

Authors: Zebin Yang; Ling Wan; Xiaodong Sun; Lin Chen; Zheng Chen
Journal: Journal of Sensors
 (2016)

---

## Abstract

For the problem of low control performance of Bearingless Induction Motor (BIM) control system in the presence of large load disturbance, a novel load torque sliding mode observer is proposed on the basis of establishing sliding mode speed control system. The load observer chooses the speed and load torque of the BIM control system as the observed objects, uses the speed error to design the integral sliding mode surface, and adds the low-pass filter to reduce the torque observation error. Meanwhile, the output of the load torque is used as the feedforward compensation for the control system, which can provide the required current for load changes and reduce the adverse influence of disturbance on system performance. Besides, considering that the load changes lead to the varying rotational inertia, the integral identification method is adopted to identify the rotational inertia of BIM, and the rotational inertia can be updated to the load observer in real time. The simulation and experiment results all show that the proposed method can track load torque accurately, improve the ability to resist disturbances, and ameliorate the operation quality of BIM control system. The chattering of sliding mode also is suppressed effectively.


---

## Body

## 1. Introduction

Based on the similarity principles of magnetic bearing and alternating current (AC) motor stator structure, BIM is formed. Two sets of windings are embedded in the stator slot of BIM, which can separately produce electromagnetic torque and radial levitation force. BIM achieves the integration of rapid rotation and stable suspension of rotor by changing the currents in the windings and avoids the mechanical bearing friction, wear and tear, and lubrication. It breaks the bottleneck of traditional asynchronous motor developing towards the higher precision and higher speed direction [1–5]. BIM has many better advantages than the traditional asynchronous motor, such as simple structure, uniform air gap, high mechanical strength, high speed, and ultrahigh speed running in the corrosion or other special environments. Therefore, it shows broad development prospect in medical equipment, transportation, national defense, and so forth [6–9]. However, BIM has the characteristics of nonlinearity, multivariability, and strong coupling. The traditional PI controller cannot acquire high-performance control for BIM when the control system is disturbed by load torque [10].Sliding mode variable structure control, as a kind of special nonlinear control, can operate in accordance with the trajectory designed by people and purposefully adjust operation according to the system status, which can gain excellent control performance. Due to the fact that the sliding mode control not only can be set by people, but also does not need high precision mathematical model and has strong robustness to disturbances, it is becoming a hot research topic [11–15], and it is gradually applied in the AC servo system. In [16], a new reaching law was designed to improve the operation quality of sliding mode. At the same time, it was applied in the speed control, which effectively enhanced the robustness of permanent magnet synchronous motor (PMSM) system. In [17], the sliding mode control combining with model reference adaptive was used to obtain the speed. The results showed that it increased the estimation precision of rotor velocity for PMSM and decreased the chattering. In [18], the sliding mode control was used in a generator based on the exercise equipment with nonlinear P-V characteristic curves. The amount of generator input current harmonic is greatly reduced. In [19], the conventional sliding mode control was united with the adaptive fuzzy backstepping scheme. The simulation proved that this method improved the performance of mismatched uncertain system. In [20], the sliding mode control dealt with the difficult problem of obtaining the counterelectromotive force, and it finally implemented the direct torque control of brushless direct current motor. In [21], the sliding mode control was used to detect the speed and position for PMSM. The experimental results proved the validity of the proposed sliding mode observer. In [22], based on the nonsingular terminal sliding mode algorithm and backstepping method, the sliding mode observer and position controller were put forward, which can estimate the torque accurately and track the position quickly. In [23], the adaptive sliding mode control for uncertain singularly perturbed nonlinear system was designed. It not only reduced the effects of uncertainty, but also guaranteed the control performance. In [24, 25], the load sliding mode observers were proposed. They diminished the adverse effects of load changes on PMSM and improved the antidisturbance ability of controlled system at some level. However, they all ignored the problem that the load changes result in the different rotational inertia and the controlled system had large chattering. Hence, the system cannot achieve the best dynamic performance.A novel sliding mode observer of load torque, of which the state variables are the speed and load torque, is proposed to suppress the impacts of the load torque changes on BIM control system. A low-pass filter used in the observer reduces the observation error of torque. Moreover, the observer as feedforward compensation for the given current alleviates the output pressure of sliding mode controller (SMC). In addition, adopting the integral identification method validly identifies the rotational inertia and improves the precision of BIM. The simulation and experimental results show that the proposed method overcomes the disadvantageous effects on the speed regulation system generated by load disturbances and strengthens the antidisturbance ability of the system.


## 2. The Dynamics Model of BIM

According to the electromagnetic field theory, the radial levitation force of BIM in thed-q coordinates can be established as [6] (1)Fx=Kψ1dis2d+ψ1qis2q,Fy=Kψ1dis2q-ψ1qis2d,where K=Km+Kl, Km=πP1P2Lm1/18lrμ0N1N2, and Kl=P1N2/2rN1; Fx and Fy are the components of the radial levitation force in x and y directions; the subscript “1” represents the torque windings, the subscript “2” represents the radial levitation force windings, “s” represents the stator, and “r” represents the rotor; P1 and P2 separately represent the pole pairs of torque windings and suspension windings; is2d and is2q are the current components of the stator in levitation force windings under the d-q axis; Lm1 is mutual inductance of the levitation force windings; l is the effective length of the rotor; r is the stator inner diameter; μ0 is the permeability of vacuum; N1 and N2, respectively, show the effective number of turns of the torque windings and the levitation force windings; and ψ1d and ψ1q are the components of flux linkage for the torque winding in the d-q coordinates, respectively.With the torque windings and the levitation force windings, BIM is a nonlinear, strongly coupled, and complex system. In order to simplify the analysis of BIM, a hypothesis is given that the levitation force windings only create a rotating magnetic field. The rotor voltage equation can be described as(2)ur1d=Rr1ir1d+pψr1d-ψ1q+Lr1lir1qω1-ωr=0,ur1q=Rr1ir1q+pψr1q-ψ1d+Lr1lir1dω1-ωr=0,where ur1d and ur1q are the rotor voltages of torque windings in d-q coordinates; Rr1 is the rotor resistance; ω1 and ωr are separately the air gap field speed and rotor speed; and p is the differential operator.The flux linkage can be expressed as(3)ψ1d=is1d+ir1dLm1,ψ1q=is1q+ir1qLm1,ψs1d=ψ1d+is1dLs1l,ψs1q=ψ1q+is1qLs1l,ψr1d=ψ1d+ir1dLr1l,ψr1q=ψ1q+ir1qLr1l,where Ls1l and Lr1l are the stator leakage inductance and rotor leakage inductance of torque windings, respectively.The electromagnetic torque equation is set up as(4)Te=P1ψ1dis1q-ψ1qis1d.The equation of motion is written as(5)Te=TL+JP1pωr,where TL is the load torque and J is the rotational inertia.After coordinate transforming, the rotor flux ind-q axis can be expressed as(6)ψr1d=Lm1is1d+Lr1ir1d,ψr1q=Lm1is1q+Lr1ir1q.Making the axis of the rotating coordinates d coincide with the rotor flux linkage of torque windings, it is written as ψr1d=ψr1. Formula (6) can be simplified as(7)ir1d=ψr1-Lm1is1dLr1,ir1q=-Lm1Lr1is1q.Putting Formula (7) into Formula (2), the excitation current is1d and slip speed ωs can be obtained as follows:(8)is1d=Tr1p+1Lm1ψr1,ωs=Lm1Tr1ψr1is1q,where ωs=ω1-ωr and Tr1=Lr1/Rr1 is the time constant of rotor.The electromagnetic torque equation turns into(9)Te=P1Lm1Lr1is1qψr1,where Lr1 is the rotor self-induction. Figure 1 is the block diagram of rotor field-oriented decoupling control.Figure 1
Rotor field-oriented decoupling control.


## 3. The Speed SMC of BIM

The system state variables are defined as(10)eω1=ω∗-ω,eω2=eω1′=-ω′,where ω∗ is the given speed and ω is the actual speed.After combining with (5), Formula (10) is described as(11)eω1′=eω2=-P12ψr1JLm1Lr1is1q+P1JTL,eω2′=eω1′′=-P12ψr1JLm1Lr1is1q′.The sliding mode surface is chosen as(12)s=eω1+c1eω2.This paper chooses the reaching law [16] to weaken the inherent chattering:(13)dsdt=-εX2sgn⁡s-kX2s,limt→∞⁡X=0,a≥0,b≥0,ε>0,k>0.According to the Lyapunov stability theory, the existence condition of generalized sliding mode is written as [12](14)V=12s2<0.Differentiating (14) with respect to time, it becomes (15) by substituting (13):(15)V˙=ss˙=s-εX2sgn⁡s-kX2s=-εX2s-kX2s2<0.Formula (15) always stands up. Hence, the system can arrive at sliding mode surface in limited time.Differentiatings=eω1+c1eω2 with respect to time, it can be gained as(16)s′=eω1′+c1eω2′=eω2-c1P12ψr1Lm1JLr1is1q′.Combining Formula (13) with Formula (16) gives the following formula:(17)eω2-c1P12ψr1Lm1JLr1is1q′=-εX2sgn⁡s-kX2s.ChoosingX=eω1 to avoid the differential interference in eω2, the sliding mode controller is designed as(18)is1q=JLr1P12ψr1Lm1c1∫εeω12sat⁡s+keω12s+eω2dt.From (18), it can be seen that the current can eliminate steady-state error and improve the accuracy of system.


## 4. The Design of Novel Load Torque Sliding Mode Observer

### 4.1. The Design of Load Torque Sliding Mode Observer

Considering the high switch frequency of the controller, the load torque can be deemed to be a constant value in a control cycle. Considering the load torque as an extension, the state equation of BIM can be expressed as(19)dωdt=P12ψr1JLm1Lr1is1q-P1JTL,dTLdt=0.Based on the equation above, the extended load torque observer is written as(20)dω^dt=P12ψr1JLm1Lr1is1q-P1JT^L+V,dT^Ldt=ηV,where V=γsgn⁡(ω-ω^); γ is the sliding mode gain; η is the feedback gain; and ω^ and T^L are the estimations of electrical angular velocity and load torque, respectively.The estimation errors of speed and the load are defined as(21)x1x2=ω-ω^TL-T^L.After Formula (19) subtracts Formula (20), the observation errors of sliding mode are obtained as(22)dx1dt=-P1Jx2-V,dx2dt=-ηV.Because the torque change is expressed in the form of speed finally, the designed sliding mode surface consists of the state variablex1=ω-ω^. The sliding mode surface is established as (23) to reduce the system overshoot:(23)s=x1+c∫x1dt.Differentiating the sliding mode surface and combining with (13) and (5), it can be acquired as(24)T^L=JP1γsgn⁡ω-ω^-cω-ω^-kω-ω^2s-εω-ω^2sgn⁡s+TL0.Because the system load torque TL is an unknown variable, TL(0) is recorded as the estimation of load torque at time zero.In order to decrease the chattering in the sliding mode, the sign functionsgn⁡(s) is replaced by saturation function sat⁡(s,Δ) [25]:(25)sats,Δ=1six>-ΔsixΔ-Δ<six<Δ-1six<-Δ.Formula (24) can be expressed as(26)T^L=JP1γsat⁡ω-ω^-cω-ω^-kω-ω^2s-εω-ω^2sats,Δ+TL0.For the purpose of improving the observation precision of sliding mode observer, a low-pass filter [7] shown in (27) is added to the observer:(27)T~L=ωcs+ωcT^L.The output of load torque observer is as the feedforward disturbance compensationiq2 and the given current iq∗ is described as(28)iq∗=is1q+iq2=JLr1P12ψr1Lm1c1∫εeω12sat⁡s+keω12s+eω2dt+T~L∂,where ∂>0 is the feedforward gain of torque observation.From (18) and (28), it can be known that ɛ and k in Formula (18) need to be large enough for meeting the load disturbance, while large ɛ and k increase the amplitude of discrete magnitude and result in big chattering. However, in (28), the observed disturbance is used to provide the required current for disturbance changes and needs no large ɛ and k. As a consequence, the feedforward compensation scheme of load torque reduces the gain amplitude of sliding mode and lowers the negative impacts on the control system caused by the disturbances.


### 4.2. The Load Torque Sliding Mode Observer with the Rotational Inertia Online Identification

When a sliding mode observer is designed, the rotational inertia of system is usually regarded as a known quantity. However, in practice application, the load changes will lead to the inertia changes. Therefore, identifying the inertia online and timely updating it to the observer will greatly improve the overall control performance of BIM system. With the characteristics of high precision and strong robustness, the integral identification algorithm is used to identify the rotational inertia online. The structure diagram of BIM’s speed loop based on the integral identification is shown in Figure2.Figure 2
Diagram of inertia identification.Equation (5) is rewritten as(29)Te=J^P1dωdt+T^r,where J^ is the estimation of rotational inertia and T^r is the collection of disturbances and its specific expression is as follows:(30)T^rt=ΔJdωdt+TL,where ΔJ is the inertia error, J-J^=ΔJ.By the formula above, it can be found that the output of disturbance is in the form of torque. At the same time, observing torque can obtainΔJ. Using ΔJ and constant recursion correction based on selecting initial value of inertia, the recursion equation can be expressed as(31)J^k=J^k-1+ΔJ.Thus, the identification accuracy of inertia depends onΔJ. In order to obtain high accuracy ΔJ, this paper chooses a periodic speed signal and uses the integral to eliminate the influences of torque disturbances on rotational inertia.Because of the high sample frequency, the load torqueTL is regarded as TL(t)=TL(t+T) and the speed signal ω meets ω(t)=ω(t+T). After on both sides of Formula (30) multiplying by ω˙(t) and integrating it, the equation is written as(32)∫k-1TkTT^rtω˙tdt=ΔJ∫k-1TkTω˙2tdt+∫k-1TkTTLtω˙tdt.Due to the fact that the load torque is a constant within a cycle, the definite integral of the load torque in (32) can be calculated as follows:(33)∫k-1TkTTLtω˙tdt=TLωtk-1TkT=TLωkT-TLωk-1T=0.Hence, (32) can be simplified as(34)∫k-1TkTT^rtω˙tdt=ΔJ∫k-1TkTω˙2tdt.From (34), it can be found that the effects of load disturbance on the inertia have been solved.Equation (34) can be expressed as(35)ΔJ=∫k-1TkTT^rtω˙tdt∫k-1TkTω˙2tdt.Combining with (35), the recursive equation of inertia identification is available as(36)J^k=J^k-1+∫k-1TkTT^rtω˙tdt∫k-1TkTω˙2tdt.Inertia identification is realized and can be updated to the load observer automatically.


## 4.1. The Design of Load Torque Sliding Mode Observer

Considering the high switch frequency of the controller, the load torque can be deemed to be a constant value in a control cycle. Considering the load torque as an extension, the state equation of BIM can be expressed as(19)dωdt=P12ψr1JLm1Lr1is1q-P1JTL,dTLdt=0.Based on the equation above, the extended load torque observer is written as(20)dω^dt=P12ψr1JLm1Lr1is1q-P1JT^L+V,dT^Ldt=ηV,where V=γsgn⁡(ω-ω^); γ is the sliding mode gain; η is the feedback gain; and ω^ and T^L are the estimations of electrical angular velocity and load torque, respectively.The estimation errors of speed and the load are defined as(21)x1x2=ω-ω^TL-T^L.After Formula (19) subtracts Formula (20), the observation errors of sliding mode are obtained as(22)dx1dt=-P1Jx2-V,dx2dt=-ηV.Because the torque change is expressed in the form of speed finally, the designed sliding mode surface consists of the state variablex1=ω-ω^. The sliding mode surface is established as (23) to reduce the system overshoot:(23)s=x1+c∫x1dt.Differentiating the sliding mode surface and combining with (13) and (5), it can be acquired as(24)T^L=JP1γsgn⁡ω-ω^-cω-ω^-kω-ω^2s-εω-ω^2sgn⁡s+TL0.Because the system load torque TL is an unknown variable, TL(0) is recorded as the estimation of load torque at time zero.In order to decrease the chattering in the sliding mode, the sign functionsgn⁡(s) is replaced by saturation function sat⁡(s,Δ) [25]:(25)sats,Δ=1six>-ΔsixΔ-Δ<six<Δ-1six<-Δ.Formula (24) can be expressed as(26)T^L=JP1γsat⁡ω-ω^-cω-ω^-kω-ω^2s-εω-ω^2sats,Δ+TL0.For the purpose of improving the observation precision of sliding mode observer, a low-pass filter [7] shown in (27) is added to the observer:(27)T~L=ωcs+ωcT^L.The output of load torque observer is as the feedforward disturbance compensationiq2 and the given current iq∗ is described as(28)iq∗=is1q+iq2=JLr1P12ψr1Lm1c1∫εeω12sat⁡s+keω12s+eω2dt+T~L∂,where ∂>0 is the feedforward gain of torque observation.From (18) and (28), it can be known that ɛ and k in Formula (18) need to be large enough for meeting the load disturbance, while large ɛ and k increase the amplitude of discrete magnitude and result in big chattering. However, in (28), the observed disturbance is used to provide the required current for disturbance changes and needs no large ɛ and k. As a consequence, the feedforward compensation scheme of load torque reduces the gain amplitude of sliding mode and lowers the negative impacts on the control system caused by the disturbances.


## 4.2. The Load Torque Sliding Mode Observer with the Rotational Inertia Online Identification

When a sliding mode observer is designed, the rotational inertia of system is usually regarded as a known quantity. However, in practice application, the load changes will lead to the inertia changes. Therefore, identifying the inertia online and timely updating it to the observer will greatly improve the overall control performance of BIM system. With the characteristics of high precision and strong robustness, the integral identification algorithm is used to identify the rotational inertia online. The structure diagram of BIM’s speed loop based on the integral identification is shown in Figure2.Figure 2
Diagram of inertia identification.Equation (5) is rewritten as(29)Te=J^P1dωdt+T^r,where J^ is the estimation of rotational inertia and T^r is the collection of disturbances and its specific expression is as follows:(30)T^rt=ΔJdωdt+TL,where ΔJ is the inertia error, J-J^=ΔJ.By the formula above, it can be found that the output of disturbance is in the form of torque. At the same time, observing torque can obtainΔJ. Using ΔJ and constant recursion correction based on selecting initial value of inertia, the recursion equation can be expressed as(31)J^k=J^k-1+ΔJ.Thus, the identification accuracy of inertia depends onΔJ. In order to obtain high accuracy ΔJ, this paper chooses a periodic speed signal and uses the integral to eliminate the influences of torque disturbances on rotational inertia.Because of the high sample frequency, the load torqueTL is regarded as TL(t)=TL(t+T) and the speed signal ω meets ω(t)=ω(t+T). After on both sides of Formula (30) multiplying by ω˙(t) and integrating it, the equation is written as(32)∫k-1TkTT^rtω˙tdt=ΔJ∫k-1TkTω˙2tdt+∫k-1TkTTLtω˙tdt.Due to the fact that the load torque is a constant within a cycle, the definite integral of the load torque in (32) can be calculated as follows:(33)∫k-1TkTTLtω˙tdt=TLωtk-1TkT=TLωkT-TLωk-1T=0.Hence, (32) can be simplified as(34)∫k-1TkTT^rtω˙tdt=ΔJ∫k-1TkTω˙2tdt.From (34), it can be found that the effects of load disturbance on the inertia have been solved.Equation (34) can be expressed as(35)ΔJ=∫k-1TkTT^rtω˙tdt∫k-1TkTω˙2tdt.Combining with (35), the recursive equation of inertia identification is available as(36)J^k=J^k-1+∫k-1TkTT^rtω˙tdt∫k-1TkTω˙2tdt.Inertia identification is realized and can be updated to the load observer automatically.


## 5. Simulation and Experimental Research

### 5.1. Results and Analysis of the Simulation

To validate the effectiveness of novel load torque observer with the function of inertia identification online and feedforward compensation scheme for BIM speed regulation system, a simulation mode of control system is constructed. The control block diagram of BIM includes two parts: the rotation part and the suspension part. In rotation part, the SMC outputs the given currentiq∗ by inputting the speed error. Combining with the given air gap flux ψ1∗, the excitation component is received. With the coordinate transformation method, the two-phase excitation current is transformed into the three-phase given current. With the current regulation, the three-phase current is obtained which is used to control the rotation of rotor. In suspension part, the radial levitation force of BIM is output by PID controllers. With ψ1∗, the current in the levitation windings can be calculated. By the coordinate transformation and current regulation, the required three-phase current is gained. The steady suspension and rapid rotation of rotor are realized finally. The whole control block diagram of BIM is shown in Figure 3 and the specific parameters of BIM are shown in Table 1.Table 1
Parameters of the Bearingless Induction Motor (BIM).
Parameters Torque winding Suspension winding Rated power (Kw) 1 0.5 Rated current (A) 2.86 2.86 Stator resistance (Ω) 2.01 1.03 Rotor resistance (Ω) 11.48 0.075 Mutual inductance of stator and rotor (H) 0.15856 0.00932 Stator leakage inductance (H) 0.16310 0.01199 Rotor leakage inductance (H) 0.16778 0.01474 Rotational inertia (kg⋅m2) 0.00769 0.00769 Rotor mass (kg) 2.85 2.85 Stator inner diameter (mm) 98 98 Core length (mm) 105 105 Pole pairs 1 2Figure 3
Control system diagram of Bearingless Induction Motor (BIM).In the simulation, the rotational inertia is set to 0.00769/kg·m2 and the given speed is n=10000r/min. Based on the tracking characteristics of integral identification algorithm, a step signal is selected as the given speed whose amplitude is 10000 r/min and the sampling cycle is defined as T=0.02s. Figure 4 shows the identification waveform of rotational inertia when the system suffers the load disturbances. It can be seen that the rotational inertia J converges to the given value within a sampling period at first. When the system is attacked by the load disturbance (8 N·m), J just has small fluctuation and restores the stabilization quickly. Therefore, the integral identification algorithm has better tracking ability and shows good robustness for disturbance.Figure 4
The identification result of rotational inertia.Figures5(a) and 5(b) show the estimation load T^L and the actual load TL under load mutation. As shown in (a) and (b), a sudden load (8 N·m) is added to the system at t=0.4 s. Then, the load drops to 0 N·m at t=0.6 s. It can be seen from Figure 5 that the sliding mode observer can accurately track the load torque and has strong robustness.Figure 5
The torque waveforms under load mutation: (a) the waveform of the estimation torque; (b) the waveform of the actual torque.
(a) (b)Figure6 shows the simulations of BIM in the presence of sudden load (8 N·m). Figures 6(a) and 6(b) present the output currents from the SMC and the load torque observer, respectively. Figure 6(c) shows the speed response of BIM under the proposed control strategy in this paper and the conventional SMC. Figure 6(a) shows that the current from SMC just slightly increases and it returns back to the stable value rapidly. Figure 6(b) indicates that the compensation current from load torque observer rises quickly, which provides enough current for disturbance. Figure 6(c) demonstrates that the speed of BIM controlled by the conventional SMC has larger fluctuation and needs more time to operate at the original speed value than the proposed method when it is attacked by load mutation. Based on the above analysis, the following two conclusions can be obtained:(1)
The feedforward compensation scheme of load torque observer can provide the required current for load changes. It can reduce the output pressure of SMC and make the output of SMC almost invariant.(2)
Based on the novel load observer and feedforward compensation strategy, the speed of BIM under big disturbance has no fluctuation and can quickly converge to the original value. The method weakens the system chattering effectively and enhances the stability of system.Figure 6
The responses of current and speed under load mutation: (a) the output current of SMC; (b) the compensation current of load torque; and (c) the speed response of BIM under sudden load.
(a) (b) (c)Figures7(a) and 7(b) show the rotor radial displacement at the speed of n=10000r/min. It can be known that the rotor can arrive at the steady point rapidly with the proposed control strategy in this paper. It achieves the integration of rapid rotation and stable suspension. The system has excellent control performance.Figure 7
Waveforms of rotor radial displacement: (a) the radial displacement onx-axis; (b) the radial displacement on y-axis.
(a) (b)


### 5.2. Results and Analysis of the Experiment

In order to further verify the effectiveness of the proposed control method, an experimental prototype with two degrees of freedom is used to build experimental platform. Due to the limits of photoelectrical encoder measuring speed, the speed is set to 2000 r/min in the experiment. The air gap of motor auxiliary bearing is 0.4 mm. Moreover, the load mutation is carried out to detect antijamming performance of BIM. The experimental results are shown in Figure8.Figure 8
The experimental waveforms underload mutation: (a) the radial displacement of rotor; (b) the identification results of rotational inertia under load mutation; (c) the waveforms of torque; (d) the responses of speed and current under the conventional SMC; and (e) the responses of speed and current under the proposed method.
(a) (b) (c) (d) (e)Figure8(a) shows the radial displacement when the system is controlled by the proposed method in this paper. The rotor is running around the equilibrium point and the maximum offset value is far less than the air gap of auxiliary bearing. It indicates that the rotor is suspended steadily under the proposed control strategy. Figures 8(b)–8(e) show the responses of BIM with the load mutation. Figure 8(b) gives the identification results of rotational inertia. In view of the excellent robustness of integral identification algorithm, the rotational inertia returns back to the stable value after a slight fluctuation. The observing waveform of load torque is given in Figure 8(c). It can be seen from the waveform that the sliding mode observer tracks the load precisely. Figures 8(d) and 8(e) present the experimental results based on the control strategy of ordinary SMC and the proposed method in this paper separately. From Figure 8(d), it can be found that the current has big undulation. The speed decreases by 40 r/min and needs 20 ms to restore stability. From Figure 8(e), the load mutation does not affect the output of SMC and the speed is smooth. By comparing (d) with (e), it can be found that using the new load torque sliding mode observer and the feedforward compensation scheme can make BIM exactly and stably operate. In addition, it has low sensitivity with respect to disturbances.


## 5.1. Results and Analysis of the Simulation

To validate the effectiveness of novel load torque observer with the function of inertia identification online and feedforward compensation scheme for BIM speed regulation system, a simulation mode of control system is constructed. The control block diagram of BIM includes two parts: the rotation part and the suspension part. In rotation part, the SMC outputs the given currentiq∗ by inputting the speed error. Combining with the given air gap flux ψ1∗, the excitation component is received. With the coordinate transformation method, the two-phase excitation current is transformed into the three-phase given current. With the current regulation, the three-phase current is obtained which is used to control the rotation of rotor. In suspension part, the radial levitation force of BIM is output by PID controllers. With ψ1∗, the current in the levitation windings can be calculated. By the coordinate transformation and current regulation, the required three-phase current is gained. The steady suspension and rapid rotation of rotor are realized finally. The whole control block diagram of BIM is shown in Figure 3 and the specific parameters of BIM are shown in Table 1.Table 1
Parameters of the Bearingless Induction Motor (BIM).
Parameters Torque winding Suspension winding Rated power (Kw) 1 0.5 Rated current (A) 2.86 2.86 Stator resistance (Ω) 2.01 1.03 Rotor resistance (Ω) 11.48 0.075 Mutual inductance of stator and rotor (H) 0.15856 0.00932 Stator leakage inductance (H) 0.16310 0.01199 Rotor leakage inductance (H) 0.16778 0.01474 Rotational inertia (kg⋅m2) 0.00769 0.00769 Rotor mass (kg) 2.85 2.85 Stator inner diameter (mm) 98 98 Core length (mm) 105 105 Pole pairs 1 2Figure 3
Control system diagram of Bearingless Induction Motor (BIM).In the simulation, the rotational inertia is set to 0.00769/kg·m2 and the given speed is n=10000r/min. Based on the tracking characteristics of integral identification algorithm, a step signal is selected as the given speed whose amplitude is 10000 r/min and the sampling cycle is defined as T=0.02s. Figure 4 shows the identification waveform of rotational inertia when the system suffers the load disturbances. It can be seen that the rotational inertia J converges to the given value within a sampling period at first. When the system is attacked by the load disturbance (8 N·m), J just has small fluctuation and restores the stabilization quickly. Therefore, the integral identification algorithm has better tracking ability and shows good robustness for disturbance.Figure 4
The identification result of rotational inertia.Figures5(a) and 5(b) show the estimation load T^L and the actual load TL under load mutation. As shown in (a) and (b), a sudden load (8 N·m) is added to the system at t=0.4 s. Then, the load drops to 0 N·m at t=0.6 s. It can be seen from Figure 5 that the sliding mode observer can accurately track the load torque and has strong robustness.Figure 5
The torque waveforms under load mutation: (a) the waveform of the estimation torque; (b) the waveform of the actual torque.
(a) (b)Figure6 shows the simulations of BIM in the presence of sudden load (8 N·m). Figures 6(a) and 6(b) present the output currents from the SMC and the load torque observer, respectively. Figure 6(c) shows the speed response of BIM under the proposed control strategy in this paper and the conventional SMC. Figure 6(a) shows that the current from SMC just slightly increases and it returns back to the stable value rapidly. Figure 6(b) indicates that the compensation current from load torque observer rises quickly, which provides enough current for disturbance. Figure 6(c) demonstrates that the speed of BIM controlled by the conventional SMC has larger fluctuation and needs more time to operate at the original speed value than the proposed method when it is attacked by load mutation. Based on the above analysis, the following two conclusions can be obtained:(1)
The feedforward compensation scheme of load torque observer can provide the required current for load changes. It can reduce the output pressure of SMC and make the output of SMC almost invariant.(2)
Based on the novel load observer and feedforward compensation strategy, the speed of BIM under big disturbance has no fluctuation and can quickly converge to the original value. The method weakens the system chattering effectively and enhances the stability of system.Figure 6
The responses of current and speed under load mutation: (a) the output current of SMC; (b) the compensation current of load torque; and (c) the speed response of BIM under sudden load.
(a) (b) (c)Figures7(a) and 7(b) show the rotor radial displacement at the speed of n=10000r/min. It can be known that the rotor can arrive at the steady point rapidly with the proposed control strategy in this paper. It achieves the integration of rapid rotation and stable suspension. The system has excellent control performance.Figure 7
Waveforms of rotor radial displacement: (a) the radial displacement onx-axis; (b) the radial displacement on y-axis.
(a) (b)


## 5.2. Results and Analysis of the Experiment

In order to further verify the effectiveness of the proposed control method, an experimental prototype with two degrees of freedom is used to build experimental platform. Due to the limits of photoelectrical encoder measuring speed, the speed is set to 2000 r/min in the experiment. The air gap of motor auxiliary bearing is 0.4 mm. Moreover, the load mutation is carried out to detect antijamming performance of BIM. The experimental results are shown in Figure8.Figure 8
The experimental waveforms underload mutation: (a) the radial displacement of rotor; (b) the identification results of rotational inertia under load mutation; (c) the waveforms of torque; (d) the responses of speed and current under the conventional SMC; and (e) the responses of speed and current under the proposed method.
(a) (b) (c) (d) (e)Figure8(a) shows the radial displacement when the system is controlled by the proposed method in this paper. The rotor is running around the equilibrium point and the maximum offset value is far less than the air gap of auxiliary bearing. It indicates that the rotor is suspended steadily under the proposed control strategy. Figures 8(b)–8(e) show the responses of BIM with the load mutation. Figure 8(b) gives the identification results of rotational inertia. In view of the excellent robustness of integral identification algorithm, the rotational inertia returns back to the stable value after a slight fluctuation. The observing waveform of load torque is given in Figure 8(c). It can be seen from the waveform that the sliding mode observer tracks the load precisely. Figures 8(d) and 8(e) present the experimental results based on the control strategy of ordinary SMC and the proposed method in this paper separately. From Figure 8(d), it can be found that the current has big undulation. The speed decreases by 40 r/min and needs 20 ms to restore stability. From Figure 8(e), the load mutation does not affect the output of SMC and the speed is smooth. By comparing (d) with (e), it can be found that using the new load torque sliding mode observer and the feedforward compensation scheme can make BIM exactly and stably operate. In addition, it has low sensitivity with respect to disturbances.


## 6. Conclusions

A novel load torque sliding mode observer was proposed to eliminate the adverse impacts caused by load disturbance in BIM control system. Owning to the application of a low-pass filter in the sliding mode observer, the observation error of load torque is effectively reduced. At the same time, the output of load torque observer as disturbance compensation greatly diminishes the amplitude of discrete quantity and weakens chattering. With the integral identification method, the proposed sliding mode observer can identify the rotational inertia accurately and improve the robustness of rotational inertia for disturbances. The simulation and experimental results all show that the proposed control scheme in this paper effectively improves the dynamic and static performance of BIM control system, suppresses the system chattering, and enhances the robustness of system.



---

*Source: 8567429-2016-10-11.xml*