Abstract

In the real applications, we found that it is difficult to achieve good control performance through manually tuning proportional–integral (PI) parameters of phase locked loop (PLL) and speed-loop of Luenberger observer (LO) for the PMSM sensorless control system. Therefore, this paper is to use the particle swarm optimization (PSO) algorithm to optimize the PI parameters of PLL and speed-loop of Luenberger observer of the system. Firstly, the ranges of PLL parameters are obtained by analyzing the PLL subsystem stability. Then, the ranges of PI parameters of PLL and speed-loop are set based on theoretical estimation and empirical values. The control system model is realized in MATLAB/Simulink that considers the constraints such as the saturation. The integral time absolute error is the objective function, and the PSO with different topologies is used to optimize the PI parameters. The simulation and experimental results show that the proposed method is feasible, and the optimized parameters can effectively improve the precision of position estimation and speed estimation. Moreover, the simulations and experiments are carried out to verify the robustness of the proposed method, and the results show that the optimized system can achieve good performance when there are uncertainties or disturbances.

1. Introduction

The advantages of permanent magnet synchronous motor (PMSM) include high efficiency, high power density, and small size [14]. With the development of science and technology, customers, who need higher quality productions and many products, put forward higher requirements for motor control and energy efficiency, which cause the widespread adoption of PMSM in industries like home appliance, industrial engineering, and autoindustry. The field-oriented control (FOC) technology is widely used to control PMSM. In the FOC system, it is necessary to obtain the information of rotor position and speed. In order to obtain the rotor position, some position sensors (such as Hall sensor, photoelectric encoder, and rotary converter) are used to detect the rotor position of the motor. However, the sensor increases the complexity and cost of control systems and reduces the reliability of the system [4]. In order to solve the series of problems caused the position sensors, researchers have proposed the sensorless control technology for PMSM. The research and application of this technology is of great significance and has gradually become one of the hottest fields of motor control, and this kind of control method is called sensorless control [5].

In general, the main methods of sensorless control of PMSM include high frequency injection (HFI), model reference adaptive system (MRAS), flux linkage observer (FLO), extended Kalman filter (EKF), sliding mode observer (SMO), and Luenberger observer (LO). HFI has good control performance on the sensorless low-speed operation [68]. However, this method has a poor effect on the recognition of the hidden pole motors besides causing additional noise [9]. MRAS has been used in references [10, 11], but it is dependent on motor parameters excessively [12]. FLO was used in [13, 14], but it has issues of DC offset and harmonics. EKF is good at dynamic performance and anti-interference ability [15, 16], but this algorithm is complex [1, 9]. SMO is simple, and it has good robustness [1720], but it has discontinuous switch control, which can cause system chattering [9]. Especially when it runs at low speed, SMO can cause huge torque ripple. LO uses the classical theory of the geometric property of linear (and stationary) systems known as observability [2123] to estimate the internal states of the system based on the input and output values of the system. It is simple and not chattering and has good robustness [21], but there is division operation in the arctangent function, which can cause larger estimation error of rotor position and speed [19]. It was verified that PLL algorithm can efficiently improve estimation accuracy of rotor position and speed without division operation [2, 20], but it is difficult to tune the PI parameters of PLL [1].

To achieve good control performance of PMSM sensorless control, this paper proposes to use PSO algorithm to optimize PI parameters of PLL and speed-loop. This method can reduce the parameter tuning difficulty and the deviation of PLL and speed-loop effectively, besides improving the speed performance of PMSM sensorless control. The main contributions of this paper are (1)this study proposed to use PSO algorithm to optimize the PI parameters of PLL and speed-loop of the Luenberger observer-based PMSM sensorless control system, which can reduce the parameters tuning challenges(2)the optimized PI parameters based on the simulation model of PMSM sensorless control system have been applied to the physical system, and it shows that they can improve the performance of position estimation and speed estimation(3)simulations and experiments have been used to verify the robustness of the proposed method, and the results show that the optimized system has good performance when there are uncertainties or disturbances(4)it is found that the random topology of PSO can achieve better performance than the ring topology and global topology

The rest of this paper is arranged as follows. Section 2 gives a brief description of the Luenberger observer-based PMSM sensorless control. A new application for the PSO-optimized PMSM sensorless control system is described in Section 3. In Section 4, the system simulation and experiment are given to verify the performance of the proposed method, and the results are investigated based on some comparisons. Finally, Section 5 gives the conclusion of this study.

2. Luenberger Observer Based on PMSM Sensorless Control System

The PMSM sensorless control system using Luenberger observer is shown in Figure 1. The whole control system is comprised of space vector pulse width modulation (SVPWM), three-phase full bridge inverter circuit, Luenberger observer, and speed and position estimator. The details of Figure 1 are described in the following subsections.

2.1. Mathematical Model of PMSM

This paper uses the stationary stator coordinate system of PMSM as shown in Figure 2. The axes , , and refer to three-phase stationary stator coordinates. Each electric angle is 120°. The axes and refer to two-phase stationary stator coordinates. The vector diagram between two-phase stationary stator coordinates and two-phase rotating rotor coordinates is shown as Figure 3. The axes and refer to two-phase rotating rotor coordinates. According to rotor field-orientated control, is the axis where the rotor locates, and the axis gets ahead 90° of the axis in an anticlockwise direction. The axis and rotate at rotor angular speed in an anticlockwise direction.

According to Figures 2 and 3, it comes up with the voltage equations that PMSM works under ABC three-phase coordinates.

Here, is the stator winding resistance; , , and are three-phase current of ; , and are three-phase voltage of ; , , and are three-phase stator flux linkage of .

The three-phase stator flux linkage equations of PMSM are

Here, is the stator winding inductance; , , and are the flux linkages that the permanent magnet makes on the stator winding.

The stator voltage equations that PMSM works under the stationary α-β coordinates are

The back EMF equations that PMSM works under the stationary coordinates are

Here, , , , and are the stator current and voltage of - axis; is the flux linkage of rotor permanent magnet; is the electrical angle speed of the rotor;is the electrical angle of the rotor.

The voltage equations that PMSM works under and coordinates are as follows: where and are the stator voltage on the axis and , respectively, and are the stator current on the axis and , respectively, and and are the stator flux linkage component on the axis and , respectively.

The stator flux linkage component on the axis and can be described by

The equation of electromagnetic torque is

Here, is the pole-pairs number of the motor.

2.2. Mathematical Model of Luenberger Observer of PMSM

The basic structure of Luenberger observer is shown in Figure 4. According to the mathematical model of PMSM, we usually use the voltage and current of the PMSM, which could be measured easily, as the inputs to restructure its states. In Figure 4, a feedback control is introduced, and the feedback signal is the comparison of the estimated current with the actual one. It makes the estimated value be close to the actual value infinitely.

In Figure 4, , , , and are the stator current and voltage of - axis; and are the estimated currents of - axis, respectively; and are the estimated back EMFs of the - axis, respectively.

In the block diagram of Figure 4, , , and .

According to Figure 4, we can obtain the following equations:

Substituting , , and into equations (8) and (9), we can design the PMSM mathematical model of the Luenberger observer as follows:

Here, and are the gains of the Luenberger observer.

After back EMF of Luenberger observer being estimated, the position angle of the rotor can be obtained directly by using the function below.

Here, refers to the estimated rotor electrical angle.

According to equations (4) and (11), the following function can be obtained.

Here, refers to the estimated electrical angle speed.

According to equation (11), the arctangent function is used to calculate rotor electrical angle. In practice, there are lots of harmonics and noise in the control system due to the nonlinearity of the inverter. Therefore, the division operation in the arctangent function causes larger estimation error of rotor position and speed [19].

Another method for rotor position and speed estimation is using PLL, which can efficiently avoid the noise made by arctangent function, and the specific block diagram is shown in Figure 5 (we call this method as “2021 LO”). The whole control system is comprised of SVPWM, three-phase full bridge inverter circuit, Luenberger observer, and phase locked loop.

PLL is used to obtain the position and speed due to the excellent tracking performance on frequency and phase of PLL. The generic PLL block diagram is shown in Figure 6.

According to Figure 6, the following function can be obtained

In Figure 6, refers to the estimated rotor speed by the PLL, and refers to the estimated rotor position angle by the PLL.

When , the following function can be used.

According to equation (14), Figure 6 can be equivalent to Figure 7.

According to Figure 7, the closed loop transfer function and error transfer function of phase locked loop can be obtained, respectively, as below:

Here, . Since is in the outside of the closed-loop system, it does not affect the stability of the PLL.

Since , (15) and (16) are both minimum phase system. So, the stability of (15) and (16) is determined by their poles.

The characteristic equation of (15) and (16) is

Hence, the poles of (15) and (16) are

Since and , if , and the real parts of are negative if .

Hence, for all , (15) and (16) are stable.

As rotor position is a ramp function, the steady-state error of rotor position and speed estimation system by phase locked loop can be obtained as

Equation (19) shows that actual rotor position can be estimated by the PLL, and the PLL has a low pass filtering characteristic, which can perform secondary filtering on the back EMF to obtain the rotor position. Because the performance of phase-locked loop is determined by and , we found that it is very complex to calculate the value of and based on the analytical solution. In order to improve the estimating performance of position and speed, PSO can be used to optimize PI parameters of PLL and speed-loop in PMSM sensorless control system based on the 2021 LO. It can also improve estimated precision of position and speed and performance of dynamic response.

3. Optimizing PMSM Sensorless Control System

In practice, the load is often fluctuant when the PMSM is used, and it can cause maladjustment or overshooting of rotation speed if the parameters of PI controllers are not properly tuned. It is hard to meet the requirements in high precision rotation speed control. Hence, it is necessary to optimize the PI values of speed-loop in PMSM sensorless control system and improve the performance of speed control. The specific block diagram of the proposed system is shown in Figure 8. The following two subsections show how to apply PSO to the Luenberger observer-based PMSM sensorless control system.

3.1. Introduction of PSO

PSO was put forward by Oubelaid et al. [24]. The information of every particle in the particle swarm is iterated and updated constantly. In every iteration, it can calculate the optimal values of individual particle and all the particles ( and ). After calculating the and , every particle adjusts its flight speed and location. The formulas of speed and location updating of particles are

In the equations, refers to velocity of particle in dimension after times’ iteration; refers to location of particle in dimension after times’ iteration; and are learning factors in particle swarm; and refer to two random numbers in (0, 1).

The topology of the particles plays important role for the optimization performance, and there are several kinds of topologies such as ring topology, global topology (also called full connection topology), and random topology (it is the default set of the standard PSO 2007 [25]). The ring topology means that each particle only has two fixed neighbors, the global topology is the above described one and its formula is (20), and the random topology means that the topology of particles are not fixed and is dynamically changed according to preset condition(s).

3.2. PSO-Based PMSM Sensorless Control System

In this paper, the optimization of the PI parameters of PLL and speed-loop is realized using the standard PSO 2007 [25]. Firstly, we need to select the objective function. In order to get a good control performance, the absolute value of speed difference between the target setting value and estimating value of observer should be small enough. The objective function is

Here, is the speed difference between the target setting value and estimating value of observer at time , and is the target setting speed in time .

In consideration of the system’s tuning time should be small, we use the ITAE (integral time absolute error) as the evaluation function. Its continuous state function is (22), and its discrete state function is (23).

Since the fixed-step solver ode3 is used to simulate the control system, we can use the following equation as the objective function which is similar to (22).

The variables to be optimized include 4 parameters, and of PLL and and of speed-loop. We label them as , , , and .

4. Simulation of System and Experimental Verification

In order to verify the feasibility and performance of the proposed control system, simulations and experiments are carried out. The parameters of PMSM, PSO, the 2021 LO, and the LO are shown in Tables 14.

4.1. Construct Simulation Model

According to the above analysis, the simulation model of PSO-based PMSM sensorless control system is shown in Figure 9 realized in MATLAB/Simulink. The whole simulation model is comprised of PMSM, SVPWM, three-phase full bridge inverter circuit, Luenberger observer, phase locked loop, and PSO. Here, PSO is used to optimize the whole system by calling the simulation file.

4.2. Comparative Analysis of PSO with Different Topologies

The parameters of PMSM in this simulation are shown in Table 1. According to foregoing analysis, the PLL is stable when both the value of and are greater than zero, but the range is too big. According to theoretical estimation [26] and experience in engineering application [27], the variable search spaces and PSO parameter set are given in Table 2, and the flowchart of PSO algorithm optimization is given in Figure 10. The whole simulation time is 0.1 second, and the target speed is 1000 RPM. When the time is at 0.05 second, we increase the load to 2 Nm. After 98 iterations of PSO, we got the optimal parameters of PI shown in Table 5. is convergent and close to zero. The value of evaluation function is 3.7386e+06.

In order to find which topology can achieve better optimization performance, we did comparisons among the ring topology, global topology, and the standard PSO 2007 (random topology). In the same condition, we simulated 10 times for each topology, and the results are shown in Table 6. From Table 6, we can conclude that the ITAE of standard PSO 2007 (the random topology) is better or smaller than that of the PSO with ring topology and global topology, but the computational time of ring topology is better than standard PSO and global topology, since there are fewer neighbors and less comparisons for the ring topology.

4.3. Simulation and Analysis of PMSM Performance

According to the simulation result of Subsection 4.2, we used the optimal parameters of PI by standard PSO 2007 (shown in Table 5) to verify the quantitative performance of the proposed method. LO [23] and 2021 LO [2] are compared with the proposed method based on the simulations. In the simulations, the same PMSM model and same simulation conditions were used. (1)Comparison of speed curve

The simulation results are shown in Figure 11. We can conclude that PSO has obvious advantages in speed dynamic performance. (2)Comparison of control performance

The simulation comparisons are based on the maximum speed overshoot, transient time and ITAE, and they are shown in Table 7. From Table 7, we can conclude that the performance of PMSM has been improved significantly using the proposed method.

4.4. Comparative Analysis of Simulation Results of PMSM Sensorless Control System

According to the simulation results of Subsection 4.2, the PI parameters optimized by PSO were applied to the simulation model of PMSM sensorless control system. The performance of PMSM was compared with the LO and the 2021 LO. The simulation condition is same as Subsection 4.2.

Figures 1214 show the back EMF waveforms using the LO, the 2021 LO and the proposed method, respectively. These waveforms are perfect sine waves of three back EMF. The start and operation processes of PMSM are quite stable, which build a good foundation for the accurate estimations of the rotor position and speed. This proves that there is no obvious difference on back EMF performances of these three observers under this simulation condition.

Figures 1517 show the estimated position and the actual position waveforms for the LO, the 2021 LO, and the proposed method. From Figures 1517, it can be found that the position angle deviation of proposed method between estimation and actual position is smaller than other two methods.

Figures 1820 show the position estimation errors using the LO, the 2021 LO, and the proposed method, respectively. Comparing these waveforms, we can get the following: (a)The position estimation errors of the LO, the 2021 LO, and the proposed method are 2.9° (about 0.051 rad), 2.18° (about 0.038 rad), and 1.58° (about 0.0275 rad), respectively, when start up(b)The position estimation errors are 2.5° (about 0.0436 rad), 1.8° (about 0.0314 rad), and 1.2° (about 0.0209 rad), respectively, in steady state

It shows that the position estimation error of PMSM is smaller using the proposed method than that of other two methods.

Figures 2123 show the estimated speed and actual speed about the LO, the 2021 LO, and the proposed method, respectively. It showed that compared with the proposed method, the speed overshoot of PMSM is bigger by using the LO and the 2021 LO, and the time from start to steady state is longer.

Figures 2426 show the speed estimation error using the LO, the 2021 LO, and the proposed method, respectively. Comparing these waveforms, it shows that the proposed method is better than the LO and the 2021 LO of the error of speed estimation when start-up, and the speed error is about ±0.1 RPM when in steady state.

Table 8 lists the simulation data on the accuracy of estimated position and estimated speed, and it can be found that the proposed method can achieve better performance than other two methods.

4.5. Robustness Analysis of PMSM Sensorless Control Based on Simulations

In practice, the robustness is important for a control system. In order to verify the robustness of the optimized PMSM system, the work performance on no-load and full-load is to be verified based on two cases.

Case 1. The whole simulation time is 0.1 second. Firstly, we set the PMSM running at 1000 RPM without load, and then, the 2 Nm load is applied to the system at second. The simulation results are shown in Figure 26, which shows the waveform of speed appeared fluctuations, but it goes stably immediately.

Case 2. The whole simulation time is 0.1 second. Firstly, we set the PMSM running at the 1000 RPM with full load 2 Nm. When second, the load and speed are changed to 1 Nm (50% of load torque changes) and 500 RPM (50% of speed changes), respectively. Simulation results of speed and speed estimation error are shown in Figures 27 and 28, respectively. The simulation results showed that the speed fluctuates slightly and then be stabilized quickly when the speed and torque are changed.

According to the simulation results above, it proves that the control system has good robustness and dynamic response.

4.6. Experimental Results and Analysis

It is important to verify the proposed PSO optimization method using physical system. Here, using STM32F030 MCU as its main control chip and writing the software code in C language, a PMSM sensorless control system is designed. The parameters of PMSM are same as Table 1. Experimental platform is shown in Figure 29.

In order to verify the accuracy of speed, 2 Nm loads are applied to the PMSM by hysteresis dynamometer when speed running at 1000 RPM, and the phase current waveform, torque, and speed are shown in Figure 30. The waveforms show that the phase current presents a perfect sine wave, and the motor speed is stabilized at about 1000.05 RPM. It shows that the speed results of experiment are consistent with the simulation results.

To verify the robustness of optimized PMSM system, two experimental cases are carried out.

Case 1. In order to verify the robustness of the PMSM physical system, 2 Nm loads are applied to PMSM by hysteresis dynamometer after the speed stabilizes at 1000 RPM without load, and the test data is shown in Figure 31. These phase currents are perfect sine waves. The current changes very smoothly, without abnormal fluctuations when the load is applied to the PMSM system.

Case 2. In order to verify the robustness of PMSM physical system under full load condition, the load and speed are changed to 1 Nm and 500 RPM, respectively, after PMSM system becomes stable with a speed of 1000 RPM and a 2 Nm load. The test data is shown in Figure 32. The waveform shows that the current changes very smoothly, without abnormal fluctuations when the test condition changes.

According to the above experiments, it proves that experimental results are consistent with the simulation results, and the optimized PMSM system has good stability and robustness. Because of the limited resources, the rotor position angle and speed overshoot are not obtained in the experiment. According to equations (4) and (13), the speed is proportional to the position angle. Combined with the speed testing data above, it can be deduced that the position angle is consistent to the simulation data.

5. Conclusion

This paper proposed to use the PSO algorithm to optimize the PI parameters of PLL and the speed-loop of the Luenberger observer-based PMSM sensorless control system. The simulation results in MATLAB/Simulink showed that the PSO algorithm can reduce the parameter tuning difficulty of PLL and speed loop effectively. Then, we applied the optimized PI parameters to the Simulink model of PMSM sensorless control system and physical system. The simulation results and experimental results showed that the proposed method can improve performance of position estimation and speed estimation. Moreover, the PSO with different topologies was tested, and it was found that the random topology of PSO can achieve better performance than the ring topology and global topology. Finally, simulations and experiments have been used to verify the robustness of the proposed method, and the results show that the optimized system has good performance when there are uncertainties or disturbances. Since the system model is based on MATLAB/Simulink, it is difficult to convert it to the codes, which can be run in the microcontroller of the experimental system. In the future, we will simplify the system model and the PSO algorithm, so the proposed method can real-time tune the system parameters.

Abbreviations

PMSM:Permanent magnet synchronous motor
PSO:Particle swarm optimization
PI:Proportional–integral
PLL:Phase locked loop
LO:Luenberger observer
FOC:Field-oriented control
HFI:High-frequency injection
MRAS:Model reference adaptive system
FLO:Flux linkage observer
EKF:Extended Kalman filter
SMO:Sliding mode observer
SVPWM:Space vector pulse width modulation
EMF:Electromotive force
ITAE:Integral time absolute error.

Data Availability

No data were used to support this study.

Conflicts of Interest

The author declares that they have no conflicts of interest.

Acknowledgments

This research is supported partially by the South African National Research Foundation Grants (Nos. 141951, 132159, 137951, and 132797).