# Control System Design for 16/6/8 Double-Stator Bearingless Switched Reluctance Motor

Authors: Chuanyu Sun; Hang Yang; Shangke Han; Hongchang Ding; Jiaqing Li; Ning Han
Journal: Mathematical Problems in Engineering
 (2021)

---

## Abstract

The 16/6/8 double-stator bearingless switched reluctance motor (DSBSRM) is used as the object of study in this paper. To solve the problem of torque and levitation force ripples in this motor, a control system direct force control (DFC) and direct instantaneous torque control (DITC) based on the torque sharing function (TSF) are proposed. With the strong nonlinearity and approximation capability of radial basis function neural networks, the torque and levitation force observer are designed. The observed torque and levitation forces are used as feedback for the internal loop control, which is combined with the external loop control to make a double closed-loop control. In order to further improve the output torque and system robustness and suppress the torque ripple in steady-state process, the motor winding method is optimized and a set of switching angles is added on the basis of TSF. The simulation results verify the effectiveness and superiority of the proposed control method. It effectively suppresses speed ripple and reduces torque and levitation force fluctuations and rotor radial displacement jitter.


---

## Body

## 1. Introduction

With the development of modern industry, the performance requirements of motors are becoming higher and higher. The traditional motor uses mechanical bearings, so the bearing has high mechanical wear, works less efficiently, and has a low life at high speeds. In order to overcome the disadvantages of mechanical bearing, scholars proposed the bearingless switched reluctance motor (BSRM) [1, 2]. It not only has the advantages of simple structure, good speed regulation performance, and high mechanical strength of switched reluctance motor (SRM), but also solves the problem of mechanical bearings. Due to the longer motor life, the absence of mechanical friction, and the lack of lubrication, it is used in a wide range of applications [3, 4]. However, BSRM has strong coupling and fluctuation problems. To solve these problems, scholars at home and abroad have conducted a lot of research mainly in two aspects: structure and control algorithm.Based on the structural characteristics, four types can be distinguished: mixed stator tooth structure, wide rotor tooth structure, double-stator structure, and composite rotor structure [5–11]. The dual stator structure, with the outer stator winding providing the torque and the inner stator winding providing the levitation, provides better decoupling and easier control by separating the torque and levitation force magnetic circuits [10, 12–14]. Current square wave control is the main control method, where the winding current is directly regulated by a current hysteresis controller to produce the required torque and levitation force [13, 15]. However, this method requires complex mathematical model calculations and has large fluctuation of levitation force and torque. At present, there is a hot topic of research in applying the control methods from SRM to BSRM [15–20]. The DTC (direct torque control) and DFC can better suppress torque and levitation force ripples by setting the space voltage vector table reasonably well. They can also improve control accuracy and motor robustness by combining neural networks or superhelix algorithms [15–18]. It was found that DITC is more effective than DTC in suppressing torque ripples. With a limited switching frequency, the DITC improves the operating efficiency of the motor by enabling smoother output torque during phase change [19, 20].In order to suppress torque and levitation force torque pulsation for 16/6/8 DSBSRM [21], this paper proposes a DTC and DFC control system based on RBF. In addition, to increase the output torque, the motor winding is optimized and a set of switching angles is added to the original TSF. Finally simulation verifies the effectiveness and superiority of the control system.


## 2. Motor Structure, Operating Principle, and Coupling Analysis

### 2.1. Motor Structure and Operating Principle

The structure of the 16/6/8 DSBSRM is shown in Figure1(a), with the outer stator, rotor, and inner stator from outside to inside. The outer stator is a torque stator with 16 teeth, and each tooth is wound with two sets of coils. There are 8 phases and the coils A1+, A1−, A2+, and A2− form phase A in series, as shown in Figure 1(b). The rotor consists of 6 torque teeth, a magnetic isolation frame, and a toroidal core. The torque teeth are set as “U” to increase the output torque. The inner stator is a levitation stator, consisting of 4 permanent magnets and 4 levitation teeth. The permanent magnets are arranged at 45°, 135°, 225°, and 315° to provide magnetic pull; the levitation teeth are arranged at 0°, 90°, 180°, and 270°, and a set of coils are wound. M+ and M− coils in series form the control winding of one phase. The control winding and the permanent magnets work together to achieve the stable levitation of the motor.Figure 1
The structure of the DSBSRM: (a) exploded view of 16/8/6 DSBSRM; (b) DSBSRM winding mode.
(a)(b)Figure2 shows the operating principle of 16/8/6 DSBSRM. As shown in Figure 2(a), when the A-phase winding is energized, two excitation loops are formed between the outer stator teeth and the rotor teeth. Based on the “principle of minimum magnetoresistance”, two magnetic pulls which drive the rotor counterclockwise are generated. If the A⟶D⟶G⟶B ⟶E⟶H⟶C⟶F-phase winding is energized in turn, the rotor rotates counterclockwise, and if the F⟶C⟶H⟶E⟶B⟶G⟶D⟶A-phase winding is energized, the rotor rotates clockwise.Figure 2
Schematic diagram of motor operation: (a) equivalent magnetic circuit model of torque winding and permanent magnet; (b) force analysis diagram of rotor; (c) equivalent magnetic circuit model of control winding.
(a)(b)(c)As shown in Figures2(a) and 2(b), when the rotor is levitated in the ideal center, the FA1 and FA2 combined torque force is 0, so the levitation force is provided only by the permanent magnet. Since the four permanent magnets are structured symmetrically, levitation forces provided by permanent magnets F1, F2, F3, and F4 will be of the same size and 90° apart in direction. When the rotor is eccentric in the −X direction, the air gap reluctance corresponding to each magnetic circuit changes, which causes FA1, FA2F1, F2, F3, and F4 combined force to be nonzero. Therefore, the inner stator has a levitation force on the rotor in the −X direction when the M-phase winding is not energized. As the M-phase winding current gradually increases, the magnetic flux density increases in the +X direction and decreases in −X direction, as shown in Figures 2(a) and 2(c). When the magnetic flux density is greater in the +X direction than in the −X direction, suspending force gradually changes from negative value to positive value, and levitation force is generated in the +X direction. Similarly, levitation force is generated in the +Y direction when N+ and −N are energized positively. By adjusting the current of the four windings appropriately, the levitation forces can be provided in any direction to satisfy levitation requirements.


### 2.2. Mathematical Analytical Model

In order to achieve stable levitation and rotation of a bearingless motor, the relationship between torque and levitation force, rotor position angle, and winding current needs to be expressed to clarify the control laws for the phase windings. The virtual-work method is adopted to derive the mathematical model of torque and levitation forces.As derived in [21], the torque can be expressed as(1)TA=2N2i2μ0hrc1δ1+4c2πrθ−δ1.Magnetic leakage factorc1 = 0.2 and c2 = −3, when θ is at [6°–7.5°]. When θ is at [7.5°–13.5°], c1 = 0.925 and c2 = 0.756. δ1 is the thickness of the air gap between the outer stator and the rotor.Similarly, ignoring magnetic leakage conditions, we can expressX direction suspending force Fx as(2)Fx=ΦPM2+μ0Aδ2Ni2−ΦPM1−μ0Aδ2Ni22μ0A,where δ2 is the thickness of the air gap between the inner stator and the rotor. x is the eccentric displacement in the X direction. ΦPM1 is the upper permanent magnet bias magnetic flux, and ΦPM2 is the lower permanent magnet bias magnetic flux. ΦPM1 and ΦPM2 can be expressed as(3)ΦPM1=2δ2μ0AFPMμ0ARPM+δ−xcos45°2δ2−x+2δ2δ2−x,ΦPM2=2δ2μ0AFPMμ0ARPM+δ2+xcos45°2δ2+x+2δ2δ2+x.


### 2.3. Coupling Analysis

To verify the decoupling performance of the 16/6/8 DSBSRM, a simulation model was built in Maxwell with the motor specific parameters shown in Table1.Table 1
Parameters of 16/6/8 DSBSRM.
Parameter nameParameter valueOuter stator outer diameter80 mmOuter stator inner diameter53 mmRotor outer diameter52 mmRotor inner diameter24 mmTorque tooth height7 mmOuter diameter of inner stator23 mmInner stator inner diameter5 mmMotor axial length80 mmExternal stator tooth pole arc7.5°Protruding pole arc of torque tooth7.5°Torque tooth pole arc30°Length of permanent magnet3 mmFigure3 shows the torque and permanent magnet flux density diagram. When the A-phase torque winding is energized with 10 A (IA = 10 A), the magnetic separator separates the magnetic field lines generated by the torque winding from the magnetic field lines generated by the permanent magnet, so that the two magnetic field lines do not intersect. It shows that the decoupling effect of torque and levitation force is very good. As shown in Figures 4(b) and 4(c), when the levitation winding of phase M is energized with 6 A (IM = 6 A), the magnetic flux density is greater in the +X direction than in the −X direction. As the M-phase winding current gradually increases, in the +X direction, the magnetic flux density X-B increases and the levitation force Fx gradually increases. However, in the +Y direction, the magnetic flux density Y-B and the levitation force Fy remain largely unchanged.Figure 3
Torque and permanent magnet flux density diagram.Figure 4
Levitation force coupling analysis diagram: (a) magnetic flux density diagram for levitation force; (b) levitation force and magnetic flux density graphs.
(a)(b)The analysis of Figures3 and 4 shows that the decoupling between torque and levitation forces works well, as well as the equally good decoupling of levitation forces in X and Y. As a result, torque and levitation forces can be controlled separately, and the levitation forces in the X and Y directions can also be controlled separately.


## 2.1. Motor Structure and Operating Principle

The structure of the 16/6/8 DSBSRM is shown in Figure1(a), with the outer stator, rotor, and inner stator from outside to inside. The outer stator is a torque stator with 16 teeth, and each tooth is wound with two sets of coils. There are 8 phases and the coils A1+, A1−, A2+, and A2− form phase A in series, as shown in Figure 1(b). The rotor consists of 6 torque teeth, a magnetic isolation frame, and a toroidal core. The torque teeth are set as “U” to increase the output torque. The inner stator is a levitation stator, consisting of 4 permanent magnets and 4 levitation teeth. The permanent magnets are arranged at 45°, 135°, 225°, and 315° to provide magnetic pull; the levitation teeth are arranged at 0°, 90°, 180°, and 270°, and a set of coils are wound. M+ and M− coils in series form the control winding of one phase. The control winding and the permanent magnets work together to achieve the stable levitation of the motor.Figure 1
The structure of the DSBSRM: (a) exploded view of 16/8/6 DSBSRM; (b) DSBSRM winding mode.
(a)(b)Figure2 shows the operating principle of 16/8/6 DSBSRM. As shown in Figure 2(a), when the A-phase winding is energized, two excitation loops are formed between the outer stator teeth and the rotor teeth. Based on the “principle of minimum magnetoresistance”, two magnetic pulls which drive the rotor counterclockwise are generated. If the A⟶D⟶G⟶B ⟶E⟶H⟶C⟶F-phase winding is energized in turn, the rotor rotates counterclockwise, and if the F⟶C⟶H⟶E⟶B⟶G⟶D⟶A-phase winding is energized, the rotor rotates clockwise.Figure 2
Schematic diagram of motor operation: (a) equivalent magnetic circuit model of torque winding and permanent magnet; (b) force analysis diagram of rotor; (c) equivalent magnetic circuit model of control winding.
(a)(b)(c)As shown in Figures2(a) and 2(b), when the rotor is levitated in the ideal center, the FA1 and FA2 combined torque force is 0, so the levitation force is provided only by the permanent magnet. Since the four permanent magnets are structured symmetrically, levitation forces provided by permanent magnets F1, F2, F3, and F4 will be of the same size and 90° apart in direction. When the rotor is eccentric in the −X direction, the air gap reluctance corresponding to each magnetic circuit changes, which causes FA1, FA2F1, F2, F3, and F4 combined force to be nonzero. Therefore, the inner stator has a levitation force on the rotor in the −X direction when the M-phase winding is not energized. As the M-phase winding current gradually increases, the magnetic flux density increases in the +X direction and decreases in −X direction, as shown in Figures 2(a) and 2(c). When the magnetic flux density is greater in the +X direction than in the −X direction, suspending force gradually changes from negative value to positive value, and levitation force is generated in the +X direction. Similarly, levitation force is generated in the +Y direction when N+ and −N are energized positively. By adjusting the current of the four windings appropriately, the levitation forces can be provided in any direction to satisfy levitation requirements.


## 2.2. Mathematical Analytical Model

In order to achieve stable levitation and rotation of a bearingless motor, the relationship between torque and levitation force, rotor position angle, and winding current needs to be expressed to clarify the control laws for the phase windings. The virtual-work method is adopted to derive the mathematical model of torque and levitation forces.As derived in [21], the torque can be expressed as(1)TA=2N2i2μ0hrc1δ1+4c2πrθ−δ1.Magnetic leakage factorc1 = 0.2 and c2 = −3, when θ is at [6°–7.5°]. When θ is at [7.5°–13.5°], c1 = 0.925 and c2 = 0.756. δ1 is the thickness of the air gap between the outer stator and the rotor.Similarly, ignoring magnetic leakage conditions, we can expressX direction suspending force Fx as(2)Fx=ΦPM2+μ0Aδ2Ni2−ΦPM1−μ0Aδ2Ni22μ0A,where δ2 is the thickness of the air gap between the inner stator and the rotor. x is the eccentric displacement in the X direction. ΦPM1 is the upper permanent magnet bias magnetic flux, and ΦPM2 is the lower permanent magnet bias magnetic flux. ΦPM1 and ΦPM2 can be expressed as(3)ΦPM1=2δ2μ0AFPMμ0ARPM+δ−xcos45°2δ2−x+2δ2δ2−x,ΦPM2=2δ2μ0AFPMμ0ARPM+δ2+xcos45°2δ2+x+2δ2δ2+x.


## 2.3. Coupling Analysis

To verify the decoupling performance of the 16/6/8 DSBSRM, a simulation model was built in Maxwell with the motor specific parameters shown in Table1.Table 1
Parameters of 16/6/8 DSBSRM.
Parameter nameParameter valueOuter stator outer diameter80 mmOuter stator inner diameter53 mmRotor outer diameter52 mmRotor inner diameter24 mmTorque tooth height7 mmOuter diameter of inner stator23 mmInner stator inner diameter5 mmMotor axial length80 mmExternal stator tooth pole arc7.5°Protruding pole arc of torque tooth7.5°Torque tooth pole arc30°Length of permanent magnet3 mmFigure3 shows the torque and permanent magnet flux density diagram. When the A-phase torque winding is energized with 10 A (IA = 10 A), the magnetic separator separates the magnetic field lines generated by the torque winding from the magnetic field lines generated by the permanent magnet, so that the two magnetic field lines do not intersect. It shows that the decoupling effect of torque and levitation force is very good. As shown in Figures 4(b) and 4(c), when the levitation winding of phase M is energized with 6 A (IM = 6 A), the magnetic flux density is greater in the +X direction than in the −X direction. As the M-phase winding current gradually increases, in the +X direction, the magnetic flux density X-B increases and the levitation force Fx gradually increases. However, in the +Y direction, the magnetic flux density Y-B and the levitation force Fy remain largely unchanged.Figure 3
Torque and permanent magnet flux density diagram.Figure 4
Levitation force coupling analysis diagram: (a) magnetic flux density diagram for levitation force; (b) levitation force and magnetic flux density graphs.
(a)(b)The analysis of Figures3 and 4 shows that the decoupling between torque and levitation forces works well, as well as the equally good decoupling of levitation forces in X and Y. As a result, torque and levitation forces can be controlled separately, and the levitation forces in the X and Y directions can also be controlled separately.


## 3. Design of Control System

According to Section2.3, the decoupling effect of torque and levitation is very good. Therefore, when the control system of 16/8/6 DSBSRM is designed, the torque control system and levitation control system are designed, respectively.


### 3.1. Design of Torque Control System

The torque characteristic is analyzed by Maxwell simulation within a rotor period of 60°, as shown in Figure5. It is found that the torque tends to oscillate and gradually decrease with increasing rotor position angle. Moreover, there are two positive torque zones: 0–15° and 27–35°. To better understand the two positive torque zones, at the same rotor cycle, the A-phase winding and the H-phase winding are enabled separately with a current of 10 A. As shown in Figure 6, when the rotor angle is right at 0–15°, the first positive torque region of A-phase intersects with the second positive torque region of H phase. Consequently, when both phase windings are on at the same time, the torque generated by the H-phase winding is superimposed on the torque generated by the A-phase winding, resulting in an increase in output torque.Figure 5
Torque characteristic curve.Figure 6
Torque superposition curve.In order to make full use of the two positive torque zones of the motor and to improve the operating efficiency of the motor, a TSF-based DITC control method was used. The method can output a specific reference torque by setting the switching angle reasonably and adjusting the torque change rate, so as to achieve smooth commutation and effectively reduce the torque fluctuation [19, 20, 22].The TSF includes linear, sine, cosine, cubic, and exponential types. The cosine function is more suitable for nonlinear systems than the linear type and has a smaller copper loss compared to other functions [22]. Therefore, it can suppress torque fluctuations well, and the cosine function is shown as follows:(4)fkθ=12−12cosπθ−θon1θov,θon≤θ<θon+θov,1,θon+θov≤θ<θoff,12+12cosπθ−θoff1θov,θoff≤θ<θoff+θov,0,else,where θ is the rotor position angle, θon is the turn-on angle, θoff is the turn-off angle, and θov is the phase current commutation overlap angle. The waveform is shown in Figure 7.Figure 7
Cosine TSF waveform.Based on the analysis of the torque characteristic curve of the motor, a set of switching angles is added to the traditional cosine function to regulates the on/off of the second positive torque region. The improved TSF is shown as follows:(5)fkθ=12−12cosπθ−θon1θov,θon1≤θ<θon1+θov,1,θon1+θov≤θ<θoff1,12+12cosπθ−θoff1θov,θoff1≤θ<θoff1+θov,12−12cosπθ−θon2θov,θon2≤θ<θon2+θov,1,θon2+θov≤θ<θoff2,12+12cosπθ−θoff2θov,θoff2≤θ<θoff2+θov,0,else,where θon1 and θoff1 are the first set of switching angles, and θon2 and θoff2 are the second set of switching angles. By setting the size of θon2 and θoff2 reasonably, it makes the H phase on when the A-phase is on at the same time. The TSF is made to have a constant sum of 1 at any point in time, and the ratio of the maximum torque in the first positive torque region to the maximum torque in the second positive torque region is about 4 : 1 when the current is 10 A. Therefore, a coefficient of 0.8 is added in the interval θon1,θoff1+θov, and a coefficient of 0.2 is added in the interval [θon2,θoff2+θov, as shown in equation (6). A schematic diagram of the waveform diagram is shown in Figure 8.(6)fθ=0.8fkθ,θon1≤θ≤θoff1+θov,0.2fkθ,θon2≤θ≤θoff2+θov,0,else.Figure 8
Waveform of improved TSF.To meet the requirements, the motor winding method is optimized. The double windings are wound in the same direction on the same outer stator teeth, as shown in Figure9. When A and H are on at the same time, the output torque increases due to the superposition of the magnetic fields. If the rotor rotates clockwise, the energizing sequence is A (H)⟶D (C)⟶G (F)⟶B (A)⟶E (D)⟶H (G)⟶C (B)⟶F (E), and if the rotor rotates counterclockwise, the energizing sequence is A (H)⟶F (E)⟶ C (B)⟶H (G)⟶E (D)⟶B (A)⟶G (F)⟶D (E).Figure 9
Winding optimization diagram.


### 3.2. Design of Levitation Force Control System

DFC uses the levitation force as the control object and adds the closed-loop buoyancy on the basis of the traditional square wave control to form a double closed-loop control. At the same time, the two-point levitation force regulator is adopted, as shown in Figure10, to realize the high dynamics of motor rotor levitation control.Figure 10
Levitation force control schematic.For hysteresis control of the levitation force, combined with Figure10, taking M-phase as an example, we assume that the rotor produces an eccentricity in the −X direction, so the levitation force in the +X direction needs to be increased. At this point, the actual levitation force F is less than the reference levitation force Fref, and the levitation force deviation is greater than the limit of the internal hysteresis loop ΔF (Fref − F > ΔF). In order to increase the levitation force, a positive voltage should be applied across the winding M to increase the current, and the voltage symbol “1” is chosen. The difference between Fxref and Fx decreases as the current increases. Similarly when the actual levitation force F is greater than the reference levitation force Fxref, the levitation force deviation is greater than the internal hysteresis loop limit ΔF (Fref − F > ΔF). In order to reduce the levitation force, a negative voltage should be applied to the winding M to increase the current; the voltage symbol “−1” is selected. This creates a cycle to achieve stable control of the levitation forces. Similarly, the regulation of the levitation force in the Y direction can be achieved by controlling the voltage sign of the levitation winding N, selected in the same way as in the X direction.


## 3.1. Design of Torque Control System

The torque characteristic is analyzed by Maxwell simulation within a rotor period of 60°, as shown in Figure5. It is found that the torque tends to oscillate and gradually decrease with increasing rotor position angle. Moreover, there are two positive torque zones: 0–15° and 27–35°. To better understand the two positive torque zones, at the same rotor cycle, the A-phase winding and the H-phase winding are enabled separately with a current of 10 A. As shown in Figure 6, when the rotor angle is right at 0–15°, the first positive torque region of A-phase intersects with the second positive torque region of H phase. Consequently, when both phase windings are on at the same time, the torque generated by the H-phase winding is superimposed on the torque generated by the A-phase winding, resulting in an increase in output torque.Figure 5
Torque characteristic curve.Figure 6
Torque superposition curve.In order to make full use of the two positive torque zones of the motor and to improve the operating efficiency of the motor, a TSF-based DITC control method was used. The method can output a specific reference torque by setting the switching angle reasonably and adjusting the torque change rate, so as to achieve smooth commutation and effectively reduce the torque fluctuation [19, 20, 22].The TSF includes linear, sine, cosine, cubic, and exponential types. The cosine function is more suitable for nonlinear systems than the linear type and has a smaller copper loss compared to other functions [22]. Therefore, it can suppress torque fluctuations well, and the cosine function is shown as follows:(4)fkθ=12−12cosπθ−θon1θov,θon≤θ<θon+θov,1,θon+θov≤θ<θoff,12+12cosπθ−θoff1θov,θoff≤θ<θoff+θov,0,else,where θ is the rotor position angle, θon is the turn-on angle, θoff is the turn-off angle, and θov is the phase current commutation overlap angle. The waveform is shown in Figure 7.Figure 7
Cosine TSF waveform.Based on the analysis of the torque characteristic curve of the motor, a set of switching angles is added to the traditional cosine function to regulates the on/off of the second positive torque region. The improved TSF is shown as follows:(5)fkθ=12−12cosπθ−θon1θov,θon1≤θ<θon1+θov,1,θon1+θov≤θ<θoff1,12+12cosπθ−θoff1θov,θoff1≤θ<θoff1+θov,12−12cosπθ−θon2θov,θon2≤θ<θon2+θov,1,θon2+θov≤θ<θoff2,12+12cosπθ−θoff2θov,θoff2≤θ<θoff2+θov,0,else,where θon1 and θoff1 are the first set of switching angles, and θon2 and θoff2 are the second set of switching angles. By setting the size of θon2 and θoff2 reasonably, it makes the H phase on when the A-phase is on at the same time. The TSF is made to have a constant sum of 1 at any point in time, and the ratio of the maximum torque in the first positive torque region to the maximum torque in the second positive torque region is about 4 : 1 when the current is 10 A. Therefore, a coefficient of 0.8 is added in the interval θon1,θoff1+θov, and a coefficient of 0.2 is added in the interval [θon2,θoff2+θov, as shown in equation (6). A schematic diagram of the waveform diagram is shown in Figure 8.(6)fθ=0.8fkθ,θon1≤θ≤θoff1+θov,0.2fkθ,θon2≤θ≤θoff2+θov,0,else.Figure 8
Waveform of improved TSF.To meet the requirements, the motor winding method is optimized. The double windings are wound in the same direction on the same outer stator teeth, as shown in Figure9. When A and H are on at the same time, the output torque increases due to the superposition of the magnetic fields. If the rotor rotates clockwise, the energizing sequence is A (H)⟶D (C)⟶G (F)⟶B (A)⟶E (D)⟶H (G)⟶C (B)⟶F (E), and if the rotor rotates counterclockwise, the energizing sequence is A (H)⟶F (E)⟶ C (B)⟶H (G)⟶E (D)⟶B (A)⟶G (F)⟶D (E).Figure 9
Winding optimization diagram.


## 3.2. Design of Levitation Force Control System

DFC uses the levitation force as the control object and adds the closed-loop buoyancy on the basis of the traditional square wave control to form a double closed-loop control. At the same time, the two-point levitation force regulator is adopted, as shown in Figure10, to realize the high dynamics of motor rotor levitation control.Figure 10
Levitation force control schematic.For hysteresis control of the levitation force, combined with Figure10, taking M-phase as an example, we assume that the rotor produces an eccentricity in the −X direction, so the levitation force in the +X direction needs to be increased. At this point, the actual levitation force F is less than the reference levitation force Fref, and the levitation force deviation is greater than the limit of the internal hysteresis loop ΔF (Fref − F > ΔF). In order to increase the levitation force, a positive voltage should be applied across the winding M to increase the current, and the voltage symbol “1” is chosen. The difference between Fxref and Fx decreases as the current increases. Similarly when the actual levitation force F is greater than the reference levitation force Fxref, the levitation force deviation is greater than the internal hysteresis loop limit ΔF (Fref − F > ΔF). In order to reduce the levitation force, a negative voltage should be applied to the winding M to increase the current; the voltage symbol “−1” is selected. This creates a cycle to achieve stable control of the levitation forces. Similarly, the regulation of the levitation force in the Y direction can be achieved by controlling the voltage sign of the levitation winding N, selected in the same way as in the X direction.


## 4. Control System Based on RBF Observer

In the proposed control system, the torque and levitation forces are used as direct control objects, so the torque and levitation forces need to be estimated. From (1) and (2), the method ignores the magnetic saturation and leakage of the motor and is too complex to derive a torque model for the full angle [21]. The neural network observer is more effective and convenient than the mathematical model approach.


### 4.1. RBF Observer

Compared with other neural networks, RBF neural networks have the advantages of simple structure, good approximation and generalization ability, and fast learning speed. They are widely used in nonlinear control systems, and their structure is shown in Figure11 [23]. There are three layers: input layer, hidden layer, and output layer, and there are no links between the same layers.Figure 11
Structure of RBF neural network.The input vector isX=x1,x2T, the hidden layer is H=h1,h2,…,hjT, the weight of the implied layer to the output layer is W=W1,W2,…,WjT, and the output is Y=WTH. The hj is the output of the j neuron of the hidden layer, as shown in the following equation:(7)hj=exp−X−cj2bj2,where cj is central vector and bj is the width of the Gaussian basis function. The relationship between phase current, rotor position angle, and torque, as shown in Figure 12(a), and the data relationship between levitation current, rotor eccentric displacement, and levitation force, as shown in Figure 13(a), were obtained by Maxwell simulation.Figure 12
Torque RBF training results: (a) simulation torque value; (b) torque training result; (c) error of simulation value and training result.
(a)(b)(c)Figure 13
RBF training results of levitation force: (a) simulation levitation force value; (b) result of levitation force training; (c) error of simulation value and training result.
(a)(b)(c)In the torque observer, the phase current and rotor position angle are used as inputs and the torque as output. The phase currents were taken in the range of 0–60 A, one value for every 2 A; the rotor position angles were taken in the range of 0–60°, one value for every 2°, for a total of 900 samples; and the input quantities were normalized. The training results are shown in Figure12(b). The levitation force observer takes levitation current and rotor eccentric displacement as input and levitation force as output. The levitation current takes values from 0 to 6 A, one value for every 0.2 A, and the eccentric displacement takes values from 0 to 0.4 mm, one value for every 0.01 mm. 1200 training samples were taken, and the training results are shown in Figure 13(b). The smaller the target error, the higher the observation accuracy, but the longer the training time, the more the neurons. Therefore, in order to improve the training efficiency, the number of neurons is reduced while ensuring the training accuracy. After repeated training, the training target errors are determined as 0.0001 and 0.035.As shown in Figures12(c) and 13(c), the simulation values of torque and levitation force are plotted against the training results. The errors are kept within ±0.05 Nm and ±0.2 Nm, respectively. The training results are very close to the simulation values, and the errors are very small, achieving the observation of torque and levitation forces.


### 4.2. Control Block Diagram

According to the above analysis results, the control block diagram of 16/6/8 DSBSRM is constructed as shown in Figure14. The asymmetric half-bridge circuit is selected as the power converter with 8 phases for torque control and 4 phases for levitation force control.Figure 14
Control block diagram.In the torque control system, the required torqueTrefA∼TrefH can be obtained from the regulation of speed error by PID and TSF. Meanwhile, the actual torque value TA ∼ TH is obtained according to the RBF torque observer. TrefA∼TrefH and TA∼TH are converted into control signals by the hysteresis loop controller, which controls the on/off of the power converter to turn the torque winding on/off for stable rotation.In the levitation control system, the required levitation forcesFxref, Fyref are obtained from the regulation of displacement errors by the PID controller. At the same time, the actual levitation force values Fx, Fy are obtained according to the RBF levitation force observer. Fxref, Fyref and Fx, Fy torque are converted into control signals by the hysteresis loop controller, which controls the on/off of the power converter to achieve stable levitation.


## 4.1. RBF Observer

Compared with other neural networks, RBF neural networks have the advantages of simple structure, good approximation and generalization ability, and fast learning speed. They are widely used in nonlinear control systems, and their structure is shown in Figure11 [23]. There are three layers: input layer, hidden layer, and output layer, and there are no links between the same layers.Figure 11
Structure of RBF neural network.The input vector isX=x1,x2T, the hidden layer is H=h1,h2,…,hjT, the weight of the implied layer to the output layer is W=W1,W2,…,WjT, and the output is Y=WTH. The hj is the output of the j neuron of the hidden layer, as shown in the following equation:(7)hj=exp−X−cj2bj2,where cj is central vector and bj is the width of the Gaussian basis function. The relationship between phase current, rotor position angle, and torque, as shown in Figure 12(a), and the data relationship between levitation current, rotor eccentric displacement, and levitation force, as shown in Figure 13(a), were obtained by Maxwell simulation.Figure 12
Torque RBF training results: (a) simulation torque value; (b) torque training result; (c) error of simulation value and training result.
(a)(b)(c)Figure 13
RBF training results of levitation force: (a) simulation levitation force value; (b) result of levitation force training; (c) error of simulation value and training result.
(a)(b)(c)In the torque observer, the phase current and rotor position angle are used as inputs and the torque as output. The phase currents were taken in the range of 0–60 A, one value for every 2 A; the rotor position angles were taken in the range of 0–60°, one value for every 2°, for a total of 900 samples; and the input quantities were normalized. The training results are shown in Figure12(b). The levitation force observer takes levitation current and rotor eccentric displacement as input and levitation force as output. The levitation current takes values from 0 to 6 A, one value for every 0.2 A, and the eccentric displacement takes values from 0 to 0.4 mm, one value for every 0.01 mm. 1200 training samples were taken, and the training results are shown in Figure 13(b). The smaller the target error, the higher the observation accuracy, but the longer the training time, the more the neurons. Therefore, in order to improve the training efficiency, the number of neurons is reduced while ensuring the training accuracy. After repeated training, the training target errors are determined as 0.0001 and 0.035.As shown in Figures12(c) and 13(c), the simulation values of torque and levitation force are plotted against the training results. The errors are kept within ±0.05 Nm and ±0.2 Nm, respectively. The training results are very close to the simulation values, and the errors are very small, achieving the observation of torque and levitation forces.


## 4.2. Control Block Diagram

According to the above analysis results, the control block diagram of 16/6/8 DSBSRM is constructed as shown in Figure14. The asymmetric half-bridge circuit is selected as the power converter with 8 phases for torque control and 4 phases for levitation force control.Figure 14
Control block diagram.In the torque control system, the required torqueTrefA∼TrefH can be obtained from the regulation of speed error by PID and TSF. Meanwhile, the actual torque value TA ∼ TH is obtained according to the RBF torque observer. TrefA∼TrefH and TA∼TH are converted into control signals by the hysteresis loop controller, which controls the on/off of the power converter to turn the torque winding on/off for stable rotation.In the levitation control system, the required levitation forcesFxref, Fyref are obtained from the regulation of displacement errors by the PID controller. At the same time, the actual levitation force values Fx, Fy are obtained according to the RBF levitation force observer. Fxref, Fyref and Fx, Fy torque are converted into control signals by the hysteresis loop controller, which controls the on/off of the power converter to achieve stable levitation.


## 5. Simulation Validation

In order to verify the feasibility and effectiveness of the above control method, a 16/6/8 DSBSRM control simulation model is built based on MATLAB/Simulink simulation platform.n∗=600rpm; x∗=y∗=0mm; θon1=2°; θoff1=10°; θon2=24°; θoff2=31°; θov=2°; the load is 1.3 Nm.Figures15(a) and 15(b) show the speed waveforms of current square wave control, traditional DITC, and modified DITC. The n∗ increases to 1200 rpm at 2 s and decreases to 600 rpm at 5 s. The speed reaches steady state at 0.23 and 2.54 s under current control and at 0.21 s and 2.31 s under traditional DITC, which is 8.6% and 42.6% improvement compared to the current control. The speed reaches steady state at 0.2 s and 2.22 s under the improved DITC. The improved DITC is a further shortening of the traditional DITC with improvement of 4.7% and 29%. At low speeds, the overshoot of the improved DITC is greater than that of the conventional DITC, due to its higher output torque. However, at high speeds, the difference is not significant.Figure 15
Dynamic response of speed: (a) speed change; (b) enlarged view of speed change; (c) steady speed of 600 rpm; (d) steady speed of 1200 rpm.
(a)(b)(c)(d)To research the output speed shaking and accuracy, Figures15(c) and 15(d) show the amplification curves of the motor speed stabilized at 600 rpm and 1200 rpm. Table 2 shows the comparison of speed fluctuation at 600 rpm and 1200 rpm for the three control methods.Table 2
Comparison of speed fluctuation.
Speed (r/min)MethodWave range (r/min)Maximum speed error (r/min)600Current control596～600.24Traditional DITC599.2～600.30.8Improved DITC599.7～600.30.31200Current control1191.4～11998.6Traditional DITC1197.8～1199.92.2Improved DITC1198.5～1199.91.5As known from Table2, at 600 rpm, the traditional DITC reduces the maximum speed error by 80% and the fluctuation amplitude by 73.8% compared to the current control. Under the improved DITC, the error and speed fluctuation are further reduced based on the conventional DITC. The maximum speed error is reduced by 62.5% and the fluctuation amplitude is reduced by 45.5%. At 1200 rpm, the improved DITC reduces the maximum speed error by 31.8% and the fluctuation amplitude by 33.3% compared with the current control. Analyzing Figures 15(c) and 15(d) and Table 2, we find that the improved DITC can further suppress speed shake and improve speed accuracy.Figures16(a) and 16(b) show the torque waveforms at 1200 rpm for both control methods. Compared to the traditional DITC, the torque fluctuation peak is smaller under the improved DITC, which effectively reduces the torque fluctuation. At the same time, when the A-phase winding provides torque, the H-phase winding provides torque at the same time. The energization law of the winding is A (H)⟶D (C)⟶G (F)⟶B (A)⟶E (D)⟶H (G)⟶C (B)⟶F (E), which is consistent with the theoretical analysis in Section 3.1.Figure 16
Various torque curves: (a) various torque curves of traditional DITC; (b) various torque curves of improved DITC.
(a)(b)Figure17 shows the torque waveform curves under the three control methods. In Figure 17(a), steady-state torque ranges are 0.05–1.62 Nm and 0.32–2.54 Nm, and pulsation amplitudes are 1.57 Nm and 2.26 Nm for 600 rpm and 1200 rpm under control. In Figure 17(b), steady-state torque ranges are 0.32–1.22 Nm and 0.65–2.45 Nm under traditional DITC, which reduces torque pulsations by 42.6% and 20.35% compared to current control. In Figure 17(c), steady-state torque ranges are 0.34–0.86 Nm and 0.79–2.18 Nm under the improved DITC. The torque pulsations of the improved DITC are reduced by 42.2% and 22.8% compared to the traditional DITC. Analyzing Figures 18(a) and 18(b), we find that the RBF torque observer has an observation error under the three control methods. With all three control methods, the observer can observe the torque values well and within a manageable margin of error. The torque fluctuation of the improved DITC is the smallest, so the error is the smallest and fluctuates within ±0.02 Nm in steady state. Analyzing Figures 17 and 18, we infer that the improved torque control system not only observes torque well, but also suppresses torque fluctuations well and improves the robustness of the motor.Figure 17
Torque ripple and observation: (a) current control; (b) traditional DITC; (c) improving the DITC.
(a)(b)(c)Figure 18
Torque observation error: (a) torque observation error at 600 rpm; (b) torque observation error at 1200 rpm.
(a)(b)Figure19 shows the steady-state displacement waveforms in the X direction controlled by current control and DFC. Under current control, the radial displacement of the rotor fluctuates at −0.107 to 0.139 mm. The fluctuation amplitude is 0.248 mm, and the maximum displacement error is 0.139 mm. Under DFC, the radial displacement of the rotor fluctuates at -0.035–0.044 mm. The fluctuation amplitude is 0.079 mm, and the maximum displacement error is 0.044 mm. The radial displacement of the rotor was reduced by 68.1%, and the maximum displacement error was reduced by 68.3%. This shows that DFC can better suppress rotor jitter and improve levitation accuracy.Figure 19
Rotor displacement inX direction.Figures20(a) and 20(b) show the waveform of the levitation force in X direction controlled by current control and DFC. The radial levitation force fluctuates at −74–70 Nm with a fluctuation of 144 Nm under current control. Under DFC, the levitation force is within −22 to 34 Nm with a fluctuation of 56 Nm. The radial levitation force fluctuation of the rotor is reduced by 61.13%. Analyzing Figure 21, we find that the RBF levitation force observers have an observation error under the two control methods. The levitation force under the DFC is smaller, so the observation error minimum is within −0.1–0.4 Nm. Analyzing Figures 20 and 21, we conclude that the improved levitation force control system not only observes levitation force well, but also suppresses levitation force fluctuations.Figure 20
Levitation force ripple and observation: (a) current control; (b) DFC.
(a)(b)Figure 21
Levitation force observation error.


## 6. Conclusions

This paper first analyzes the operating principle and decoupling effect of the 16/6/8 DSBSRM, and the control system of DITC and DFC based on TSF is proposed. Meanwhile, RBF neural network is used to build a torque and levitation force observer to form a double closed-loop control. Then, torque ripple is reduced by optimizing the motor winding method and proposing a double switching angle DITC. Finally, the simulation is verified by the MATLAB/Simulink simulation platform. The results show that the control system not only is able to observe torque and levitation forces well, but also effectively suppresses speed jitter. At the same time, torque and levitation force ripples are reduced, and rotor radial displacement jitter is decreased. This increases the dynamic response time and robustness of the system.



---

*Source: 4727917-2021-08-31.xml*