Mathematical Problems in Engineering

Volume 2013 (2013), Article ID 768736, 12 pages

http://dx.doi.org/10.1155/2013/768736

Research Article

A Two-stage Kalman Filter for Sensorless Direct Torque Controlled PM Synchronous Motor Drive

Boyu Yi, Longyun Kang, Kai Jiang, and Yujian Lin

Electric Power College, South China University of Technology, Guangzhou, Guangdong 510640, China

Received 5 December 2012; Accepted 15 March 2013

Academic Editor: Jun-Juh Yan

Copyright © 2013 Boyu Yi et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

This paper presents an optimal two-stage extended Kalman filter (OTSEKF) for closed-loop flux, torque, and speed estimation of a permanent magnet synchronous motor (PMSM) to achieve sensorless DTC-SVPWM operation of drive system. The novel observer is obtained by using the same transformation as in a linear Kalman observer, which is proposed by C.-S. Hsieh and F.-C. Chen in 1999. The OTSEKF is an effective implementation of the extended Kalman filter (EKF) and provides a recursive optimum state estimation for PMSMs using terminal signals that may be polluted by noise. Compared to a conventional EKF, the OTSEKF reduces the number of arithmetic operations. Simulation and experimental results verify the effectiveness of the proposed OTSEKF observer for DTC of PMSMs.

1. Introduction

Owing to their characteristics of high efficiency, high power density, and reliability, AC machines and more recently especially permanent magnet synchronous machines (PMSMs) have obtained dominance [1]. To control the torque and flux levels of PMSM rotor flux oriented control is seen as an industry standard [2, 3]. Alternatives however exist, and direct torque control (DTC) can provide accurate fast torque control. Moreover, DTC uses no current controller and no motor parameters other than the stator resistance, which yields a faster torque response and lower parameter dependence than with field oriented control [4–6].

To implement the idea of the DTC on PMSM, the controller must obtain some real-time parameters of PMSM, including stator flux, position, and rotor speed. So the mechanical sensor which can measure rotor position and speed is necessary for the realization of DTC. But the application of mechanical sensor increases the complexity of the control system and degrades the performance of control system when encountering adverse environmental conditions. Moreover, it increases the system cost and maintenance requirements. So sensorless DTC technique has become the hot issue in research and drawn many researchers’ and engineer s’ attention.

Several parameter estimation techniques have been reported in the literature [7–13]. One of major methods is based on the linear or nonlinear state observers, such as extended Luenberger observer and extended Kalman filter (EKF). Some works [10–13] have shown the applicability of the EKF to estimate parameters for electric machines. These observers have been shown to be the best computer algorithms for processing noisy discrete measurements and obtaining high-accuracy estimates of dynamic system states (stator flux, rotor position, and speed). The EKF-based algorithms described in [10–13] have good performance, but they appear to be very complex because of the high order of the mathematical models. The applicability of the extended Kalman filter to real-time signal processing problems is generally limited by the complex mathematical operations.

Therefore, this paper presents an optimal two-stage extended Kalman filter (OTSEKF) for flux, position, and speed estimation of a direct torque controlled (DTC) PMSM drive. The proposed observer is an effective implementation of extended Kalman filters. Following the same approach as given in [14], the two-stage structure employed in OTSEKF can decouple the conventional EKF into two parallel observers, which are called the full order filter and the augmented state filter. Compared to the conventional EKF, OTSEKF can reduce the computational burden. To facilitate the understanding, the complete equations of this filter are presented and compared to a straight implementation of the conventional EKF equations.

The paper is organized in six sections. In Section 2, the continuous and the discrete model of the PMSM are recalled. In Section 3, the DTC-SVPWM strategy of PMSMs is introduced briefly. In Section 4, based on a conventional EKF algorithm, optimal two-stage extended Kalman filter equations are developed, and theoretical operation requirements of both OTSEKF and EKF are calculated. In Section 5, experimental and simulation results are discussed. Finally, a conclusion wraps up the paper.

2. The Model of PMSM

As elaborated in [10], the machine equations in the rotor (dq) reference frame are as follows:

where , , , and are the stator voltages and currents in the () reference frame, , are stator flux linkages in the reference frame, and are the machine axes inductances, is the stator winding resistance, and is the flux produced by the magnets. The angular velocity is measured in electrical radians per second. is the electrical position.

In order to remove constant term in (1), the state vector is chosen to be , estimated parameter vector to be , input vector to be , and output to be . The PMSM model is described by the general nonlinear state-space model [15, 16]:

with

For digital implementation of the observer, the discretized machine equations are required, provided that the input vector is nearly constant during a sampling period . These equations can be obtained from (2):

Tolerating a small discretization error, a first-order series expansion of the matrix exponential is used:

with

3. Principle of Sensorless DTC-SVPWM

The basic idea of DTC technique is to calculate and control stator flux linkage and torque of PMSM directly to achieve high dynamic performance. DTC technique involves statorflux vector, torque estimators, hysteresis controllers, and switching tables in order to determine directly an inverter switching state [17, 18]. The hysteresis controller minimizes the flux and torque errors [17, 18]. However, in the conventional DTC system, the switchover among the basic vectors is discontinuous because the universal voltage inverter has only eight available basic space vectors. In a control period, only one voltage space vector can be selected so the flux and torque ripples are unavoidable.

To reduce the ripples of the electromagnetic torque and flux linkage in PMSM drives, a modified DTC using Space Vector Pulse Width Modulation (SVPWM) method (called DTC-SVPWM) is proposed in this paper. SVPWM techniques have several advantages that are offering better DC bus utilization, lower torque ripples, lower total harmonic distortion in the AC motor current, and lower switching loss. The same flux and torque estimators as for basic DTC are also used in the DTC-SVPWM method. The main difference is that DTC-SVPWM has two PI controllers and a Reference Flux Vector Calculator (RFVC) instead of hysteresis controllers and the switching table [19–21]. The system structure of DTC-SVPWM can be built and shown in Figure 1. This system uses three-way closed-loop control of speed, flux linkage, and torque. Adopting speed deviation as input value, outer loop Proportional Integral (PI) controller outputs torque reference input . Then taking torque deviation as input value, torque loop PI controller outputs , which is deviation of rotating speed between stator flux linkage and rotor flux linkage. Therefore, the reference speed of stator flux can be obtained by adding the rotor flux linkage speeds and . RFVC is used to determine the reference stator flux linkage vector for the next sample time. The inputted to RFVC is amplitude of given value of the stator flux linkage.

Figure 1: System diagram of the SVPWM-DTC scheme.

The vector relationships between stator flux linkage vector and rotor flux linkage can be drawn in the rotor flux , stator flux , and stationary frames as shown in Figure 2. , which is the phase angle of stator flux linkage, can be obtained by the flux estimator. The increment of stator flux angle in the next sampling time can be acquired by .

Figure 2: The vector diagram of PMSM.

Define the flux deviations between and as; then

where and are components of in the () reference frame. and are amplitudes of and .

In order to make up for flux deviations and , reference stator voltages and should be applied on the motor which can be calculated by

Substituting (7) into (8), (9) can be acquired:

Based on stator voltage components and , voltage vector selection, duration time, and switch signal of inverter can be obtained through SVPWM module.

4. Design of OTSEKF Observer

To enhance the reliability and robustness of the overall system, the sensorless control of DTC-SVPWM is desirable. So a state observer is necessary to achieve sensorless operation of a PMSM drive. The linear algorithm proposed by Hsieh and Chen [14], which is named the optimal two-stage Kalman estimator (OTSKE), is extended to the nonlinear estimation case. Then, the optimal two-stage extended Kalman filter (OTSEKF) is proposed to estimate the stator flux, position and rotor speed needed for DTC-SVPWM. This algorithm can effectively save computation cost compared to conventional EKF.

4.1. Conventional EKF

Based on discretized machine equations, an EKF is constructed to estimate the stator flux, torque, and speed of the PMSM. The state vector is chosen to be . , and , are chosen as input and output vectors because these quantities can be easily obtained from measurements of stator phase currents and voltage construction using DC link voltage and switching status. Considering the noise and parameter errors, the state space model in the rotor reference frame is described by

with

where and are zero-mean noise with covariances and , respectively, and are independent from the system state . The system noise takes into account system disturbance and model inaccuracies, while represents the measurement noise. The noise covariance matrices are defined as follows:

The overall structure of the EKF is well known by employing a two-step prediction and correction algorithm [10]. Hence, the filter is given by

with

The EKF mentioned above takes into account the system and measurement noise and exhibits excellent robustness to model inaccuracies, measurement noise, and system uncertainties. But it is obviously that four-order matrix operations are necessary to complete EKF computational operation. Considering Pulse Width Modulation (PWM) period is very small, only high performance DSP or FPGA can qualify for this work.

4.2. OTSEKF Algorithm

Using the same processing method as in [14], the OTSEKF can be obtained by making coordinate transformation. So it is necessary to define a transformation matrix ; the is specified as follows:

The main advantage of using the transformation is that the inverse transformation involves only a change of sign. Two blending matrices and are defined, respectively, by and . The transformation operation can be achieved by two transformation matrices and so that the variance-covariance matrices in new base are block diagonal:

Using the two transformation matrices defined above, overlined expressions that correspond to vectors and matrices in the new base can be obtained:

where

Considering characteristic of matrix , (17) become

To decouple the EKF into a full order filter and an augmented state filter, the following two-step iterative substitution method is used.

Step 1. Substituting (13) into (19), we have

Step 2. Substituting (17) into the right-hand side of (20), the following equations are written:

Supposing variance-covariance matrices are block diagonal, the following relations are obtained by using (22) and (25):

where

The above equations lead to

The equations of full order filter are acquired by the next steps.

From (21), we give the prediction equations of and as

where

From (22), we have

Then using (27), (29), and (32), the above equation can be written as

where

From (23) and (28), we have

Then

From (24), (28), and (30), we have

Then

where

From (25), (28), and (30), we have

Then

Besides prediction function (32), the rest of augmented state filter is obtained by expanding (21)–(25) and using (26)–(28):

where and can be calculated by (31) and (35).

Finally, based on the above analysis, the OTSEKF algorithm can be organized by the next two parts [15, 16]. The first part of OTSEKF for state and parameter prediction is as follows:

The second part for state and parameter correction is as follows:

Using (17), the original state can be obtained by the sum of the state with the augmented state :

The initial conditions of this OTSEKF are established with the initial conditions of a conventional EKF (, , , , ), so that

According to variables of full order filter ( and ), stator flux linkage in machine rotor reference frame are obtained as follows:

The estimated flux linkages , are transformed into the stationary reference frame via

The stator flux linkage and electromagnetic torque estimators are then given by

where is the electromagnetic torque and is the motor pole pairs. The estimated speed obtained from the OTSEKF observer is used to close the speed loop to achieve sensorless operation. So the sensorless DTC-SVPWM scheme is illustrated in Figure 3.

Figure 3: Block diagram of the proposed sensorless DTC-SVPWM scheme.

4.3. Algorithm Complexity

Tables 1 and 2 show the number of arithmetic operations required at each time sample by the standard EKF algorithm and the OTSEKF, respectively (where rough matrix-based implementation is used). The state vector , measurement , input vector , and parameter dimensions are, respectively, , , , and . To compute standard EKF, the total number of multiplications is 358, and the total number of additions is 289, whereas the number of arithmetic operations required by OTSEKF algorithm is 284 multiplications and 215 additions. This means that the number of multiplications is reduced by 20.67% and the number of additions is reduced by 25.61%.

Table 1: EKF arithmetic operation requirements.

Number of multiplications

(, , )
Number of additions

(, , )

, ,
Function of system (22)
Function of system (5)

(24)
(20)

(128)
(112)

(16)
(16)

(72)
(56)

(96)
(80)

Total
358
289

Table 2: OTSEKF arithmetic operation requirements.

Number of multiplications

(, , , )
Number of additions

(, , , )

, , , , , ,
Function of system (20)
Function of system (3)

(16)
(12)

(16)
(12)

(16)
(16)

(24)
(24)

(16)
(12)

(4)
(2)

(20)
(18)

(8)
(8)

(8)
(8)

(8)
(4)

(16)
(12)

(32)
(20)

(16)
(12)

(12)
(12)

(16)
(8)

(16)
(12)

(12)
(16)

(8)
(4)

Total
284
215

5. Simulation and Experimental Results

5.1. Simulation Results

To verify validity and feasibility of OTSEKF algorithm, the sensorless DTC-SVPWM technique described in the previous sections is implemented in MATLAB/SIMULINK environment. As stated in [14], the EKF and OTSEKF are mathematically equivalent without requiring system constraints. So a comparison is made between outputs of two algorithms. Real-time parameter estimations observed by OTSEKF are used to formulate the closed loop, such as stator flux, rotor speed, and position. The estimations obtained by conventional EKF algorithm are not included in the DTC-SVPWM and only evaluated in open loop. The block diagram of evaluation module is shown in Figure 4. Figure 5 is the simulation results of parameter estimation based on OTSEKF. The speed dynamic response of the sensorless drive is shown in Figure 5(a). The machine is accelerated from 0 rpm to 600 rpm. During transients, OTSEKF has good tracking performance (speed and position). In Figures 5(g) and 5(h), the torque and flux responses show that ripples are significantly suppressed due to the SVPWM modulation scheme. In Figures 5(c) and 5(f), difference in speed and position estimations between the two observers is very little. It proves equivalence of both observers. The difference between two observers is caused by accuracy loss in OTSEKF, which uses more calculation steps.

Figure 4: Open-loop and closed-loop simulation structure.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

Figure 5: Simulation results for parameters estimation. (a) Actual rotor speed and estimated rotor speed (OTSEKF). (b) OTSEKF speed error. (c) Difference of rotor speed estimation between EKF and OTSEKF. (d) Actual rotor position and estimated rotor position (OTSEKF). (e) OTSEKF position error. (f) Difference of rotor position estimation between EKF and OTSEKF. (g) Estimated torque. (h) Estimated stator flux linkage.

5.2. Experimental Results

The drive system consists of a three-phase, eight-pole, 1.2 kW, PMSM driven from a three-phase Insulated Gate Bipolar Transistor (IGBT) inverter. The PMSM is mechanically coupled to a magnetic clutch, which provides rated torque, even at very low speeds. The drive system is implemented by an Expert 3 control system of Myway Company. The core DSP processor in Expert 3 control system is TMS320C6713 with a clock speed of 225 MHz. The actual rotor position and speed are obtained from an incremental encoder with 10000 pulses per revolution. The inverter space-vector PWM control, current regulator, and observer algorithm operate with a 100 us sampling/integration time step. The currents flowing in the stator windings are measured with two hall-effect current sensors.

The experiment test is conducted to testify the tracking performance of OTSEKF and to demonstrate that the OTSEKF is mathematically equivalent to EKF. Figure 6 presents the corresponding experimental results. Figures 6(a) and 6(b) show the rotor speed and angle dynamic response of OTSEKF and the errors between the estimated and real values when the command speed is 600 rpm. We can observe that the OTSEKF tracks the actual speed and angle very well, and the steady state errors are very small. Figure 6(e) presents trajectory of the estimated stator flux components, and the stator flux trajectory is almost circular. Figure 6(f) shows the response of torque is fast and smooth. The differences between the two observers are represented in Figures 6(c) and 6(d). As we can expect, the differences in speed and position estimations are still small, and this proves that the two observers are mathematically equivalent.

(a)

(b)

(c)

(d)

(e)

(f)

Figure 6: Experimental results for parameters estimation. (a) Real speed, estimated speed (OTSEKF), and estimated speed error. (b) Real rotor position, estimated rotor position (OTSEKF), and estimated position error. (c) Difference of speed estimation between EKF and OTSEKF. (d) Difference of rotor position estimation between EKF and OTSEKF. (e) Estimated torque. (f) Estimated stator flux linkage.

6. Conclusion

The real-time digital implementation of extended Kalman filters algorithm for senseless DTC-SVPWM requires a very fast signal processor specialized and optimized to perform complex mathematical calculations and manipulate a large amount of data. This is generally considered as a major disadvantage to apply Kalman filter algorithm to industrial field. So this paper proposes a novel extended Kalman filter called optimal two-stage extended Kalman filter. OTSEKF can reduce computational cost effectively compared to conventional EKF. The simulation and experimental results show that OTSEKF has good tracking performance. Equivalence is verified by little difference in estimations between OTSEKF and EKF. The OTSEKF algorithm provides a new idea to solve computational complexion and allows us to use cheaper DSP in practical application.

References

T. Sebastian, G. R. Slemon, and M. A. Rahman, “Modelling of permanent magnet synchronous motors,” IEEE Transactions on Magnetics, vol. 22, no. 5, pp. 1069–1071, 1986.

M. W. Degner, J. M. Guerrero, and F. Briz, “Slip-gain estimation in field-orientation-controlled induction machines using the system transient response,” IEEE Transactions on Industry Applications, vol. 42, no. 3, pp. 702–711, 2006.

M. Rashed, P. F. A. MacConnell, A. F. Stronach, and P. Acarnley, “Sensorless indirect-rotor-field-orientation speed control of a permanent-magnet synchronous motor with stator-resistance estimation,” IEEE Transactions on Industrial Electronics, vol. 54, no. 3, pp. 1664–1675, 2007.

L. Zhong, M. F. Rahman, W. Y. Hu, and K. W. Lim, “Analysis of direct torque control in permanent magnet synchronous motor drives,” IEEE Transactions on Power Electronics, vol. 12, no. 3, pp. 528–536, 1997.

L. Zhong, M. F. Rahman, W. Y. Hu, K. W. Lim, and M. A. Rahman, “A direct torque controller for permanent magnet synchronous motor drives,” IEEE Transactions on Energy Conversion, vol. 14, no. 3, pp. 637–642, 1999.

Y. Wang and Z. Deng, “Improved stator flux estimation method for direct torque linear control of parallel hybrid excitation switched-flux generator,” IEEE Transactions on Energy Conversion, vol. 27, no. 3, pp. 747–756, 2012.

K. B. Lee and F. Blaabjerg, “An improved DTC-SVM method for sensorless matrix converter drives using an overmodulation strategy and a simple nonlinearity compensation,” IEEE Transactions on Industrial Electronics, vol. 54, no. 6, pp. 3155–3166, 2007.

Y. Yan, J. G. Zhu, and Y. G. Guo, “Initial rotor position estimation and sensorless direct torque control of surface-mounted permanent magnet synchronous motors considering saturation saliency,” IET Electric Power Applications, vol. 2, no. 1, pp. 42–48, 2008.

Y. Wang and Z.-Q. Deng, “A position sensorless method for direct torque control with space vector modulation of hybrid excitation flux-switching generator,” IEEE Transactions on Energy Conversion, vol. 27, no. 4, pp. 912–921, 2012.

T. J. Vyncke, R. K. Boel, and J. A. A. Melkebeek, “On the stator flux linkage estimation of an PMSM with extended Kalman filters,” in Proceedings of the 5th IET International Conference on Power Electronics, Machines and Drives (PEMD '10), pp. 1–6, April 2010.

G. Foo, S. Sayeef, and M. F. Rahman, “SVM direct torque controlled interior permanent magnet synchronous motor drive using an extended Kalman filter,” in Proceedings of the 4th IET International Conference on Power Electronics, Machines and Drives (PEMD '08), pp. 712–716, April 2008.

M. Barut and R. Demir, “Bi input-extended Kalman filter based speed-sensorless direct torque control of IMs,” in Proceedings of the 19th International Conference on Electrical Machines (ICEM '10), pp. 1–5, September 2010.

I. M. Alsofyani, N. Idris, T. Sutikno, and Y. A. Alamri, “An optimized Extended Kalman Filter for speed sensorless direct troque control of an induction motor,” in Proceedings of the IEEE International Conference on Power and Energy (PECon '12), pp. 319–324, December 2012.

C.-S. Hsieh and F.-C. Chen, “Optimal solution of the two-stage Kalman estimator,” IEEE Transactions on Automatic Control, vol. 44, no. 1, pp. 194–199, 1999.

M. Hilairet, F. Auger, and E. Berthelot, “Speed and rotor flux estimation of induction machines using a two-stage extended Kalman filter,” Automatica, vol. 45, no. 8, pp. 1819–1827, 2009.

A. Akrad, M. Hilairet, and D. Diallo, “A sensorless PMSM drive using a two stage extended Kalman estimator,” in Proceedings of the 34th Annual Conference of the IEEE Industrial Electronics Society (IECON '08), pp. 2776–2781, November 2008.

M. Depenbrock, “Direct self-control (DSC) of inverter-fed induction machine,” IEEE Transactions on Power Electronics, vol. 3, no. 4, pp. 420–429, 1988.

I. Takahashi and T. Noguchi, “A new quick-response and high-efficiency control strategy of an induction motor,” IEEE Transactions on Industry Applications, vol. IA-22, no. 5, pp. 820–827, 1986.

D. Swierczynski, M. P. Kazmierkowski, and F. Blaabjerg, “DSP based direct torque control of permanent magnet synchronous motor (PMSM) using space vector modulation (DTC-SVM),” in Proceedings of the IEEE International Symposium on Industrial Electronics (ISIE '02), vol. 3, pp. 723–727, 2002.

S. Sayeef and M. F. Rahman, “Comparison of proportional+integral control and variable structure control of interior permanent magnet synchronous motor drives,” in Proceedings of the 38th IEEE Annual Power Electronics Specialists Conference (PESC '07), pp. 1645–1650, June 2007.

J. Siahbalaee, S. Vaez-Zadeh, and F. Tahami, “A new loss minimization approach with flux and torque ripples reduction of direct torque controlled permanent magnet synchronous motors,” in Proceedings of the 13th European Conference on Power Electronics and Applications (EPE '09), pp. 1–8, September 2009.