Abstract:
A method and a system in connection with a speed and position sensorless permanent magnet synchronous machine equipped with an output filter and driven by an inverter, which method comprises the steps of forming a speed-adaptive full-order observer based on the dynamic model of the combination of the permanent magnet synchronous machine (PMSM) and the output filter, measuring the inverter output current (i A ), estimating inverter output current (î A ), determining the estimate for the electrical angular speed ({circumflex over (ω)} m ) using the estimated and measured inverter output currents (î A , i A ) in an adaptation law, injecting a voltage signal (u c ) into the inverter voltage reference (u A,ref 0 ) to obtain a modified voltage reference (u A,ref ), detecting an error signal (ε) from the measured inverter output current (i A ), and calculating the speed correction term (ω ε ) used in the adaptation of the observer from the error signal (ε).

Description:
BACKGROUND OF THE INVENTION 
   The invention relates to a method and system in connection with permanent magnet synchronous machines (PMSM). More specifically, the invention relates to a method of determining the angular speed and the rotor position of a speed and position sensorless PMSM drive equipped with an output filter and to an apparatus for carrying out the method. 
   Problems may be encountered in AC motor drives due to the non-sinusoidal voltage produced by a pulse-width modulated (PWM) inverter. The high rate of change of the voltage (i.e. high du/dt) may cause excessive voltage stresses in the stator winding insulations. It may also excite the parasitic capacitances of the stator winding and produce bearing currents. Lower-order harmonics cause acoustic noise and power losses; the losses caused by eddy currents are a special concern in high-speed solid-rotor motors. 
   A common approach to overcome these problems is to use an inverter output filter [1]-[4]. An LC filter, having the resonance frequency below the switching frequency, is a typical choice for the filter topology if a nearly sinusoidal output voltage is required. If a conventional scalar control method is used, the heavy filtering of the LC filter does not complicate the drive control. When better dynamic properties are demanded, a vector control method must be used. For vector control, the filter dynamics should be taken into account in the control design. 
   Various methods have been proposed for the vector control of variable-speed drives equipped with an LC filter [2]-[8]. Methods based on a feed-forward action and a sliding mode control are proposed for compensating for the effects of the filter in a speed-sensorless permanent magnet synchronous motor (PMSM) drive [2]. A full-order observer for a PMSM is implemented in the stator reference frame for estimating the rotor position in [3]. A feed-forward controller is used in a high-speed synchronous reluctance motor drive with an LC filter in [4]. In these methods, stator current or stator voltage measurements are needed. Vector control methods for induction motor drives with an LC filter are proposed in [5, 6] also require measurements from the motor side of the filter. 
   Due to the LC filter, the electrical quantities of the motor input differ from those of the inverter output. Frequency converters are equipped with measurements of the inverter output electrical quantities, but the motor terminal quantities are needed for vector control. In [2]-[6], the motor voltages or currents are measured by additional sensors, requiring hardware modifications in the motor drive. If the motor quantities are estimated instead, as proposed in [7, 8] for induction motor drives, the additional measurements are avoided and a filter can be added to an existing drive. 
   BRIEF DESCRIPTION OF THE INVENTION 
   It is an object of the present invention to provide a method that avoids the above-mentioned drawbacks and enables determining the rotor speed and position without additional measurements in a PMSM drive equipped with an LC filter. This object is achieved by a method and a system of the invention that are characterized by what is stated in independent claims  1  and  6 . 
   The invention is based on the idea that the speed and position of the rotor can be accurately determined even at low speeds by using an adaptive full-order observer in combination with pulsating high-frequency signal injection [9]. The signal injection method can be surprisingly used in connection with an LC output filter, although the LC filter is situated on the current path between the inverter and the permanent magnet synchronous machine. 
   By using the method of the invention, the PMSM can be controlled without any speed or position sensors using vector control, thus eliminating the need of mechanical vulnerable instruments and additional wiring. Further, the method does not require any additional voltage or current measurements. Only the voltage of the intermediate circuit of the frequency converter and the output current of the inverter need to be measured. The method of the invention provides accurate information on the rotor position and speed and thus enables accurate control of PMSM drives equipped with an LC filter. 

   
     BRIEF DESCRIPTION OF THE FIGURES 
     The invention will now be described in more detail by means of preferred embodiments and with reference to the attached drawings, in which 
       FIG. 1  shows a PMSM drive system equipped with a three-phase LC filter; 
       FIG. 2  shows a simplified block diagram of a control system used in connection with the PMSM drive; 
       FIG. 3  shows amplitude and phase responses of the PMSM from d-axis voltage to d-axis current; and 
       FIG. 4  shows amplitude and phase responses of the PMSM from d-axis voltage to q-axis current. 
   

   DETAILED DESCRIPTION OF THE INVENTION 
   In the following subsections, first a filter and machine model will be briefly discussed. Then an example of a control system suitable to be used in connection with the invention is described, after which the structure of the speed-adaptive full-order observer and the high-frequency signal injection used in the method of the invention is described. 
   Filter and Motor Models 
     FIG. 1  shows a PMSM drive system equipped with an LC filter  14 . The inverter output voltage u A  is filtered by the LC filter consisting of filter inductances L f  and filter capacitances C f , resulting in a nearly sinusoidal stator voltage u s . The inverter output current i A  and the dc-link voltage u dc  are the only measured quantities. The dc-link voltage is used in the control of the inverter itself. The dc-link voltage is formed by using a diode bridge  11  and inverter  12  is used to form alternating voltage from the dc-link voltage. The speed reference signal ω m,ref  is given as an input to the control  13  of the drive. 
   In the d-q reference frame fixed to the rotor, the model of the three-phase LC filter and PMSM can be written as
 
 {dot over (x)}=Ax+B[u   A ψ pm ] T   (1)
 
i A =Cx  (2)
 
where x=[i A  u s  ψ s ] T  is the state vector consisting of the inverter output current i A =[i Ad  i Aq ] T , the stator voltage u s =[u sd  u sq ] T , and the stator flux linkage ψ s =[ψ sd  ψ sq ] T . The inverter output voltage u A =[u Ad  u Aq ] T  and the permanent magnet flux ψ pm =[ψ pm  0] T  are considered as inputs to the system. The matrix transpose is denoted by superscript T. The system matrices in (1) and (2) are
 
                 A   =     [               -     R   Lf       ⁢     L   f     -   1       ⁢   I     -       ω   m     ⁢   J               -     L   f     -   1         ⁢   I         0               C   f     -   1       ⁢   I             -     ω   m       ⁢   J             -     C   f     -   1         ⁢     L   s     -   1                 0       I             -     R   s       ⁢     L   s     -   1         -       ω   m     ⁢   J             ]             (   3   )               B   =     [             L   f     -   1       ⁢   I         0           0           C   f     -   1       ⁢     L   s     -   1                 0           R   s     ⁢     L   s     -   1               ]             (   4   )               C   =     [         I       0       0         ]             (   5   )               
where L f  is the inductance and R Lf  is the series resistance of the filter inductor, C f  is the filter capacitance, R S  is the stator resistance, ω m  is the electrical angular speed of the rotor, and
 
             I   =     [         1       0           0       1         ]       ,     J   =     [         0         -   1             1       0         ]             
The inductance matrix
 
             L   s     =     [           L   d         0           0         L   q           ]           
consists of the direct-axis inductance L d  and quadrature-axis inductance L q .
 
Control System
 
     FIG. 2  shows a simplified block diagram of the control system using the estimates obtained with the invention. In  FIG. 2 , the estimated quantities are marked by the symbol ^. The cascade control and speed-adaptive full-order observer  25  are implemented in the estimated rotor reference frame. The estimated rotor position {circumflex over (θ)} m  is obtained by integrating the estimated rotor angular speed {circumflex over (ω)} m . The inverter current i A , the stator voltage u s , and the stator current i s  are controlled by PI controllers  21 ,  22 ,  23  respectively, and cross-couplings due to the rotating reference frame are compensated for. A maximum torque per current method is used for calculating the stator current reference i s,ref . The rotor speed is governed by a PI controller  24  with active damping. 
   Space vectors on the left-hand side of coordinate transformations  27 ,  28  are in the estimated rotor reference frame and on the right-hand side in the stator reference frame. Double lines in  FIG. 2  indicate vector quantities whereas single lines indicate scalar quantities. The coordinate transformations are carried out by rotating the vector values with an estimated rotor position angle {circumflex over (θ)} m . This estimated value is obtained from the observer system and its correctness is an important factor in the accuracy of the control. The control system of  FIG. 2  will be further referred to in the following paragraphs. The dc-voltage measurement, as shown in  FIG. 1  and needed for the inverter control, is not shown in  FIG. 2 . 
   Observer Structure 
   A speed-adaptive full-order observer has been successfully used in a sensorless induction motor drive equipped with an LC filter [7, 8]. A similar observer structure is constructed for the PMSM drive in the following. For stable operation at low speeds, the observer is augmented with an HF signal injection technique similarly as in [10]. 
   HF Signal Injection 
   The HF signal injection method is based on [9]. A carrier excitation signal 
                   u   c     =         u   ^     c     ⁢       cos   ⁡     (       ω   c     ⁢   t     )       ⁡     [         1           0         ]                 (   6   )               
having amplitude û c  and angular frequency ω c  is superimposed on the inverter voltage reference u A,ref 0  in the estimated rotor reference frame as shown in  FIG. 2 . The inverter output current i A   s  is measured, an HF current response is detected on the q-axis of the estimated rotor reference frame, and the amplitude is modulated by the rotor position estimation error. Thus the larger the error of the estimate of the rotor position is, the larger the amplitude of the current response is on the q-axis. If the error is zero and the estimate of the rotor position corresponds to the actual rotor position, the injected voltage does not produce any current to q-axis. The measured current is transformed to rotor reference frame with a transformation block  28 , which multiplies the measured current by e −J{circumflex over (θ)}     m   . The q-axis current signal is then bandpass filtered, demodulated, and low-pass filtered in an error signal generation block  30  to obtain error signal
   ε≈K   ε sin(2{tilde over (θ)} m )  (7) 
where K ε  is the signal injection gain and {tilde over (θ)} m  is the estimated error of the rotor position. Without the inverter output LC filter, the signal injection gain would be
 
   
     
       
         
           
             
               
                 
                   K 
                   ɛ 
                 
                 = 
                 
                   
                     
                       
                         u 
                         ^ 
                       
                       c 
                     
                     
                       ω 
                       c 
                     
                   
                   ⁢ 
                   
                     
                       
                         
                           L 
                           q 
                         
                         - 
                         
                           L 
                           d 
                         
                       
                       
                         4 
                         ⁢ 
                         
                           L 
                           q 
                         
                         ⁢ 
                         
                           L 
                           d 
                         
                       
                     
                     . 
                   
                 
               
             
             
               
                 ( 
                 8 
                 ) 
               
             
           
         
       
     
   
   The error signal (7) is used as a correction in the adaptive full-order observer described in the following subsection. 
   The LC filter affects the signal injection gain (8). To illustrate the effect, amplitude and phase responses from the inverter voltage to the inverter current were calculated by Control System Toolbox of the MATLAB® software. 
   Parameters given in Table 1 were used for this example. The response from the d-axis voltage to the d-axis current is shown in  FIG. 3  for the rotor position estimation error {tilde over (θ)} m =10°. The amplitude response has a notch at the resonance point of the filter capacitor and the d-axis inductance of the PMSM, and a peak at the filter resonance point. Above f=500 Hz, the LC filter amplifies the response compared to the response obtained by the PMSM only. 
     FIG. 4  shows the amplitude and phase response from the d-axis voltage to the q-axis current for the rotor position estimation error {tilde over (θ)} m =10°. Compared to  FIG. 3 , the parallel resonance of the filter capacitance and the motor inductance is not visible. Instead, the peak exists in the amplitude response. After the peak, the amplitude response decays rapidly when the frequency increases. It is reasonable to select the frequency of the HF excitation voltage below the resonant frequency of the filter. If the margin between the signal injection frequency and the LC filter resonant frequency is small, the HF current amplitude, and hence also the signal injection gain, are increased by the LC filter. 
   Speed-Adaptive Full-Order Observer 
   The adaptive full order observer is based on the dynamic model of the system, the inverter current being the measured feedback signal for the observer. The electrical angular speed of the rotor is estimated using an adaptation mechanism. The observer is defined by
 
 {circumflex over ({dot over (x)}=Â{circumflex over (x)}+{circumflex over (B)}[u   A {circumflex over (ψ)} pm ] T   +K ( i   A   −î   A )  (9)
 
where the system matrix and the observer gain matrix are
 
                   A   ^     =       [             -       R   ^     Lf       ⁢       L   ^     f     -   1       ⁢   I             -       L   ^     f     -   1         ⁢   I         0                 C   ^     f     -   1       ⁢   I         0           -       C   ^     f     -   1         ⁢       L   ^     s     -   1                 0       I           -       R   ^     s       ⁢       L   ^     s     -   1               ]     -       (         ω   ^     m     -     ω   ɛ       )     ⁡     [         J       0       0           0       J       0           0       0       J         ]                 (   10   )                 B   ^     =     [               L   ^     f     -   1       ⁢   I         0           0             C   ^     f     -   1       ⁢       L   ^     s     -   1                 0             R   ^     s     ⁢       L   ^     s     -   1               ]             (   11   )               K   =     [               k     1   ⁢   d       ⁢   I     +       k     1   ⁢   q       ⁢   J                     k     2   ⁢   d       ⁢   I     +       k     2   ⁢   q       ⁢   J                     k     3   ⁢   d       ⁢   I     +       k     3   ⁢   q       ⁢   J             ]             (   12   )               
and ^ denotes estimated quantity. The factors k id  and k iq  (i=1, 2, 3) are scalar gain parameters, of which k 1d  has to be a large positive constant (several p.u.) to ensure stability. The observer defined in (9) produces estimates for the states of inverter output current vector î A , stator voltage vector û s  and stator flux linkage vector {circumflex over (ψ)} s  these being the state variables in vector {circumflex over (x)}. As seen in  FIG. 2 , the observer system  25  outputs estimates for stator voltage vector û s , stator current vector is, electrical angular speed of the rotor {circumflex over (ω)} m  and rotor position {circumflex over (θ)} m . As described earlier, the estimates for rotor speed, stator current and stator voltage are used as feedback signals in the control of the system in the embodiment of  FIG. 2 .
 
   The observer defined in (9) receives as inputs measured inverter output current i A , inverter output voltage u A  and estimated permanent magnet flux {circumflex over (ψ)} pm . The inverter output voltage u A  is usually determined directly as being the voltage reference that is used to control the inverter. Inverters are usually so fast in operation that the voltage reference given to them is implemented very accurately. In the embodiment of  FIG. 2 , the observer  25  receives inverter voltage reference u A,ref 0  as the input u A . The permanent magnet flux {circumflex over (ψ)} pm  of equation (9) is a constant value estimated before the start of the drive of during the commissioning of the drive for example. 
   The measured inverter output current i A  is used directly in the observer defined in equation (9) as a feedback signal providing correction to the observer. The estimated inverter output current is subtracted from the measured one to provide current error. This error is multiplied by the matrix K as defined in (9). 
   As seen from  FIG. 2 , the estimated flux linkage vector {circumflex over (ψ)} s  is not used in the control as such. However, the estimated flux linkage vector is used in the calculation of an estimate of stator current vector î s . This estimate is calculated from the flux estimates with equation is î s ={circumflex over (L)} s   −1 ({circumflex over (ψ)} s −{circumflex over (ψ)} pm ). The control system could also be constructed in such a way that it would utilize estimated flux linkage directly. Also, it is possible to construct the observer system such that it would use stator current as a state variable. 
   The adaptation law used to correct the adaptive system to correspond with the real system is
 
{circumflex over (ω)} m   =−K   p ( i   Aq   −î   Aq )− K   i ∫( i   Aq   −î   Aq ) dt   (13)
 
where K p  and K i  are nonnegative adaptation gains. Equation (13) thus produces an estimate for the electrical angular speed of the rotor based on the estimated and measured value of the inverter output current. Equation (13) is basically a PI controller having the error between the estimated and the measured inverter output current as the error that is to be minimized.
 
   The speed correction term ω ε  is obtained by a PI mechanism
 
ω ε =ν p ε+ν i   ∫εdt   (14)
 
from the error signal ε obtained from signal injection structure, where ν p  and ν i  are nonnegative gains. The speed correction term ω ε  and the electrical angular speed {circumflex over (ω)} m  are used in the system matrix Â to adapt the observer system with the controlled real system and to correct the angular speed estimate to correspond with the actual value. The feedback and corrections made to both equations (9) and (10) affect the observed values obtained from the observer. Since the angular speed is corrected by using signal injection technique and feedback is obtained from the measured current, the observed values and the speed and position information are accurate and the whole control of the drive is reliable even at low speeds despite the LC filter.
 
   The change in the signal injection gain (8) caused by the LC filter should be taken into account when selecting gains for equation (14). These gains can be determined roughly from the responses of  FIGS. 3 and 4 , and more specifically by simulations or experiments. 
   The digital implementation of the adaptive full-order observer can be based on a simple symmetric Euler method [11], for example. 
   It is obvious to a person skilled in the art that the basic idea of the invention can be implemented in many different ways. The control system described above is only presented to give an example how the method and system according to the present invention can be taken into use. The structure of the speed adaptive observer can also be implemented in a plurality of ways. The invention and its embodiments are thus not restricted to the examples described above, but can vary within the scope of the claims. 
   
     
       
             
           
             
             
             
             
           
             
             
             
           
             
           
             
             
             
             
           
         
             
               TABLE 1 
             
             
                 
             
           
           
             
               Motor Parameters 
             
           
        
         
             
                 
               Stator resistance R s   
               3.59 
               Ω 
             
             
                 
               Direct-axis inductance L d   
               36.0 
               mH 
             
             
                 
               Quadrature-axis inductance L q   
               51.0 
               mH 
             
             
                 
               Permanent magnet flux ψ pm   
               0.545 
               Vs 
             
             
                 
               Nominal power 
               2.2 
               kW 
             
             
                 
               Nominal current I N   
               4.3 
               A 
             
             
                 
               Nominal frequency f N   
               75 
               Hz 
             
           
        
         
             
                 
               Number of pole pairs p 
               3 
             
           
        
         
             
               Filter Parameters 
             
           
        
         
             
                 
               Inductance L f   
               5.1 
               mH 
             
             
                 
               Capacitance C f   
               6.8 
               μF 
             
             
                 
               Series resistance R Lf   
               0.1 
               Ω 
             
             
                 
                 
             
           
        
       
     
   
   REFERENCES 
   
       
       [1] Y. Murai, T. Kubota, and Y. Kawase, “Leakage current reduction for a high-frequency carrier inverter feeding an induction motor,”  IEEE Trans. Ind. Applicat ., vol. 28, no. 4, pp. 858-863, July/August 1992 
       [2] M. Carpita, D. Colombo, A. Monti, and A. Fradilli, “Power converter filtering techniques design for very high speed drive systems,” in  Proc. EPE&#39; 01, Graz, Austria, August 2001 
       [3] T. D. Batzel and K. Y. Lee, “Electric propulsion with sensorless permanent magnet synchronous motor: implementation and performance,”  IEEE Trans. Energy Conversion , vol. 20, no. 3, pp. 575-583, September 2005 
       [4] J.-D. Park, C. Khalizadeh, and H. Hofmann, “Design and control of high-speed solid-rotor synchronous reluctance drive with three-phase LC filter,” in  Conf Rec. IEEE - IAS Annu. Meeting , Hong Kong, China, October 2005, pp. 715-722 
       [5] R. Seliga and W. Koczara, “Multiloop feedback control strategy in sinewave voltage inverter for an adjustable speed cage induction motor drive system,” in  Proc. EPE&#39; 01, Graz, Austria, August 2001, CD-ROM 
       [6] M. Kojima, K. Hirabayashi, Y. Kawabata, E. C. Ejiogu, and T. Kawabata, “Novel vector control system using deadbeat-controlled PWM inverter with output LC filter,”  IEEE Trans. Ind. Applicat ., vol. 40, no. 1, pp. 162-169, January/February 2004 
       [7] J. Salomäki, M. Hinkkanen, and J. Luomi, “Sensorless vector control of an induction motor fed by a PWM inverter through an output LC filter,”  Trans. IEEJ , vol. 126-D, no. 4, pp. 430-437, April 2006 
       [8] - - - , “Sensorless control of induction motor drives with inverter output filter,” in  Proc. IEEE IEMDC&#39; 05, San Antonio, Tex., May 2005, pp. 332-339 
       [9] M. Corley and R. D. Lorenz, “Rotor position and velocity estimation for a salient-pole permanent magnet synchronous machine at standstill and high speeds,”  IEEE Trans. Ind. Applicat ., vol. 43, no 4, pp 784-789, July/August 1998 
       [10] A. Piippo and J. Luomi, “Adaptive observer combined with HF signal injection for sensorless control of PMSM drives,” in  Proc. IEEE IEMDC&#39; 05, San Antonio, Tex., May 2005, pp. 674-681. 
       [11] J. Niiranen, “Fast and accurate symmetric Euler algorithm for electromechanical simulations,” in  Proc. Electrimacs&#39; 99, vol. 1, Lisboa, Portugal, September 1999, pp. 71-78.