# Patent

## Title
Speed sensorless control of an induction machine using a PWM inverter with output LC filter

## Abstract

A method of controlling of an induction machine using an inverter is disclosed with output LC filter. The method includes determining an inverter output current vector (iA), determining an inverter output voltage vector (uA), and forming a full-order observer having a system matrix (Â) and gain vector (K). The observer produces an estimated rotor flux linkage vector ({circumflex over (ψ)}R), an estimated stator current vector (îs), an estimated stator voltage vector (ûs) and an estimated inverter output current vector (îA). An estimation error (iA−îA) of the inverter output current vector is determined, and a speed adaptation law is formed based on the estimation error of the inverter output current vector for determining an estimate for electrical angular speed ({circumflex over (ω)}m) of the induction machine. The induction machine is controlled based on the produced estimates and a measured inverter output current.

## Background
This application claims priority under 35 USC § 119 to EP Patent Application No. 04104340.7, the disclosure of which is hereby incorporated by reference in its entirety.
FIELD OF THE INVENTION
This invention relates to a speed sensorless control method for an induction motor that is supplied by a PWM inverter through an output LC filter.
BACKGROUND OF THE INVENTION
A voltage source PWM inverter enables stepless speed and torque control of AC motors, allowing reduced energy consumption and increased control performance. The use of a PWM inverter, however, not only brings advantages but also causes unwanted effects in the motor. The output voltage of the inverter consists of sharp-edged voltage pulses, producing bearing currents and high voltage stresses in motor insulations [1], [2]. The oscillation at the switching frequency causes additional losses and acoustic noise. These phenomena can be eliminated by adding an LC filter to the output of the PWM inverter. In addition, the EMI shielding of the motor cable may be avoided if the voltage is nearly sinusoidal.
Adding an LC filter to a variable speed drive makes the motor control more difficult. Usually, a simple volts-per-hertz control method is chosen. Better control performance is achieved by using vector control, i.e. field oriented control. However, there are only few publications that deal with the vector control of a motor fed through an LC filter [3]–[5]. In these papers, an extra current or voltage measurement was necessary, and a speed encoder was used. In order to obtain cost savings and reliability improvements, a full-order observer was proposed in [6], making additional current or voltage measurements unnecessary.
Recently, the speed sensorless control of ac motors has become popular. Promising estimation methods for speed sensorless induction motor drives are speed-adaptive full-order observers [7] combined with improvements in regeneration mode operation [8], [9]. However, a speed sensorless control methods for induction motor supplied through an LC filter are not yet published.
A problem associated with the prior art control systems with an output filter is the need for either the measurements of current or voltage from the motor or the use of speed encoders. Both of the above solutions can increase costs both for installation and maintenance on the system.
BRIEF DESCRIPTION OF THE INVENTION
An exemplary method disclosed herein is based on the use of an adaptive full-order observer, and no additional voltage, current or speed measurements are needed for the vector control of the motor. The rotor speed adaptation is based on the estimation error of the inverter current. The rotor speed adaptation can be based on the measured inverter current due to the surprisingly noticed fact, that the quadrature components of the inverter current and the stator current are almost identical, i.e. the applied LC filter does not distort the q component of the current.
An exemplary advantage of the method is that the induction motor can be controlled without any additional measurements even when LC filter is used in the inverter output. With the method of the invention the benefits achieved with the inverter output LC filter can be utilized in connection with a drive that uses no additional measurement or feedback signals.
This can be achieved by using a speed-adaptive observer, which is extended for the induction motor drive equipped with an LC filter, resulting in a drive where only the inverter output current and dc-link voltage are measured. A simple observer gain is used, and a speed-adaptation law basing on the estimation error of the inverter current is employed. The regeneration mode operation at low speeds is further stabilized by modifying the speed adaptation law. The vector control of the motor described in the specification is based on nested control loops.



BRIEF DESCRIPTION OF THE DRAWINGS
In the following the invention will be described in greater detail by means of preferred embodiments with reference to the attached drawings, in which
FIG. 1 shows a principle of the exemplary control system,
FIG. 2 is an exemplary signal flow diagram of the cascade control,
FIGS. 3a, 3b and 4 show loci of a current estimation error,
FIG. 5 is an exemplary signal flow diagram of linearized model,
FIGS. 6a and 6b show observer poles,
FIG. 7 represents an exemplary simulation result,
FIG. 8 shows voltage and current waveforms from the simulation of FIG. 7,
FIG. 9 shows the exemplary experimental setup, and
FIGS. 10 and 11 represent exemplary experimental results.



DETAILED DESCRIPTION OF THE INVENTION
In the following first the system model and control are explained. Then attention is given to dynamic analysis of the system and simulation and experimental results are also described.
System Model and Control
A principle of the control system is shown in FIG. 1. The inverter output voltage uA is filtered by an LC filter, and the induction motor (IM) is fed by the filtered voltage us. The inverter output current iA and the dc link voltage udc are the only measured quantities, whereas the stator voltage us, stator current is, and the electrical angular speed ωm of the rotor are estimated by an observer (the estimated quantities being marked by ‘^’). The system is controlled by nested control loops in estimated rotor flux reference frame. It should be noted that the control arrangement presented is only one example of suitable control.
A. Filter and Motor Models
In a reference frame rotating at angular frequency ωs, the equations for the LC filter are




  
    
      
        
          
            d
            
              
                i
                _
              
              A
            
          
          
            d
            t
          
        
        =
        
          
            
              -
              j
            
            ⁢
            
              
            
            ⁢
            
              ω
              s
            
            ⁢
            
              
                i
                _
              
              A
            
          
          -
          
            
              
                R
                Lf
              
              
                L
                f
              
            
            ⁢
            
              
                i
                _
              
              A
            
          
          +
          
            
              1
              
                L
                f
              
            
            ⁢
            
              (
              
                
                  
                    u
                    _
                  
                  A
                
                -
                
                  
                    u
                    _
                  
                  s
                
              
              )
            
          
        
      
    
    
      
        (
        1
        )
      
    
  








  
    
      
        
          
            d
            
              
                u
                _
              
              s
            
          
          
            d
            t
          
        
        =
        
          
            
              -
              j
            
            ⁢
            
              
            
            ⁢
            
              ω
              s
            
            ⁢
            
              
                u
                _
              
              s
            
          
          +
          
            
              1
              
                C
                f
              
            
            ⁢
            
              (
              
                
                  
                    i
                    _
                  
                  A
                
                -
                
                  
                    i
                    _
                  
                  s
                
              
              )
            
          
        
      
    
    
      
        (
        2
        )
      
    
  




where Lf is the inductance and RLf the series resistance of the inductor, and Cf is the capacitance of the filter.
The motor model is based on the inverse-Γ model [10] of the induction motor. The stator and rotor voltage equations are




  
    
      
        
          
            u
            _
          
          s
        
        =
        
          
            
              R
              s
            
            ⁢
            
              
                i
                _
              
              s
            
          
          +
          
            
              d
              
                
                  ψ
                  _
                
                s
              
            
            
              d
              t
            
          
          +
          
            j
            ⁢
            
              
            
            ⁢
            
              ω
              s
            
            ⁢
            
              
                ψ
                _
              
              s
            
          
        
      
    
    
      
        (
        3
        )
      
    
  








  
    
      
        0
        =
        
          
            
              R
              R
            
            ⁢
            
              
                i
                _
              
              R
            
          
          +
          
            
              d
              
                
                  ψ
                  _
                
                R
              
            
            
              d
              t
            
          
          +
          
            
              j
              ⁡
              
                (
                
                  
                    ω
                    s
                  
                  -
                  
                    ω
                    m
                  
                
                )
              
            
            ⁢
            
              
                ψ
                _
              
              R
            
          
        
      
    
    
      
        (
        4
        )
      
    
  




respectively, where Rs and RR are the stator and rotor resistances, respectively, and iR is the rotor current. The stator and rotor flux linkages are

ψs=(L′s+LM)iS+LMiR  (5)

ψR=LM(is+iR)  (6)

respectively, where L′s denotes the stator transient inductance and LM is the magnetizing inductance. Based on (1)–(6), the state-space representation of the system can be written as shown in (7) and (8).




  
    
      
        
          
            x
            _
          
          .
        
        =
        
          
            
              
                [
                
                  
                
                ⁢
                
                  
                    
                      
                        
                          -
                          
                            
                              R
                              Lf
                            
                            
                              L
                              f
                            
                          
                        
                        -
                        
                          j
                          ⁢
                          
                            
                          
                          ⁢
                          
                            ω
                            s
                          
                        
                      
                    
                    
                      
                        -
                        
                          1
                          
                            L
                            f
                          
                        
                      
                    
                    
                      0
                    
                    
                      0
                    
                  
                  
                    
                      
                        1
                        
                          C
                          f
                        
                      
                    
                    
                      
                        
                          -
                          j
                        
                        ⁢
                        
                          
                        
                        ⁢
                        
                          ω
                          s
                        
                      
                    
                    
                      
                        -
                        
                          1
                          
                            C
                            f
                          
                        
                      
                    
                    
                      0
                    
                  
                  
                    
                      0
                    
                    
                      
                        1
                        
                          
                            L
                            ′
                          
                          s
                        
                      
                    
                    
                      
                        
                          -
                          
                            1
                            
                              
                                τ
                                ′
                              
                              σ
                            
                          
                        
                        -
                        
                          j
                          ⁢
                          
                            
                          
                          ⁢
                          
                            ω
                            s
                          
                        
                      
                    
                    
                      
                        
                          1
                          
                            
                              L
                              ′
                            
                            s
                          
                        
                        ⁢
                        
                          (
                          
                            
                              1
                              
                                τ
                                r
                              
                            
                            -
                            
                              j
                              ⁢
                              
                                
                              
                              ⁢
                              
                                ω
                                m
                              
                            
                          
                          )
                        
                      
                    
                  
                  
                    
                      0
                    
                    
                      0
                    
                    
                      
                        R
                        R
                      
                    
                    
                      
                        
                          -
                          
                            1
                            
                              τ
                              r
                            
                          
                        
                        -
                        
                          j
                          ⁡
                          
                            (
                            
                              
                                ω
                                s
                              
                              -
                              
                                ω
                                m
                              
                            
                            )
                          
                        
                      
                    
                  
                
                ⁢
                
                  
                
                ]
              
              
                }
                
                  A
                  _
                
              
            
            ⁢
            
              x
              _
            
          
          +
          
            
              
                (
                
                  
                
                ⁢
                
                  
                    
                      
                        1
                        
                          L
                          f
                        
                      
                    
                  
                  
                    
                      0
                    
                  
                  
                    
                      0
                    
                  
                  
                    
                      0
                    
                  
                
                ⁢
                
                  
                
                )
              
              
                }
                B
              
            
            ⁢
            
              
                u
                _
              
              A
            
          
        
      
    
    
      
        (
        7
        )
      
    
  








  
    
      
        
          
            i
            _
          
          A
        
        =
        
          
            
              [
              
                1
                ⁢
                
                  
                
                ⁢
                0
                ⁢
                
                  
                
                ⁢
                0
                ⁢
                
                  
                
                ⁢
                0
              
              ]
            
            
              }
              C
            
          
          ⁢
          
            x
            _
          
        
      
    
    
      
        (
        8
        )
      
    
  




The state vector is x=[iA uS is ψR]T, and the two time constants are defined as τ′π=L′s/(Rs+RR) and πr=LM/RR.
B. Cascade Control
FIG. 2 illustrates the proposed cascade control of the system in the estimated rotor flux reference frame. In the LC filter control, the innermost control loop governs the inverter current iA by means of a PI controller, and the stator voltage us is governed by a P-type controller in the next control loop. In both control loops, decoupling terms are used to compensate the cross-couplings caused by the rotating reference frame.
The motor control forms the two outermost control loops. The stator current is is controlled by a PI-type controller with cross-coupling compensation, and the rotor speed is governed by a PI-controller. In addition, a PI-type rotor flux controller is used. It should be noted, however, that the presented control system is provided only as an example.
C. Adaptive Full-Order Observer
The system states are estimated by means of a full-order observer. The electrical angular speed of the rotor, included in the state space representation (7), is estimated using an adaptation mechanism. The observer is implemented in the estimated rotor flux reference frame, i.e., in a reference frame where {circumflex over (ψ)}R={circumflex over (ψ)}R+j0. The observer is given as

{circumflex over ({dot over (x)}=Â{circumflex over (x)}+BuA+K(iA−îA)  (9)

The system matrix and the observer gain vector in (9) are




  
    
      
        
          
            A
            ^
          
          _
        
        =
        
          [
          
            
          
          ⁢
          
            
              
                
                  
                    -
                    
                      
                        R
                        Lf
                      
                      
                        L
                        f
                      
                    
                  
                  -
                  
                    j
                    ⁢
                    
                      
                    
                    ⁢
                    
                      ω
                      s
                    
                  
                
              
              
                
                  -
                  
                    1
                    
                      L
                      f
                    
                  
                
              
              
                0
              
              
                0
              
            
            
              
                
                  1
                  
                    C
                    f
                  
                
              
              
                
                  
                    -
                    j
                  
                  ⁢
                  
                    
                  
                  ⁢
                  
                    ω
                    s
                  
                
              
              
                
                  -
                  
                    1
                    
                      C
                      f
                    
                  
                
              
              
                0
              
            
            
              
                0
              
              
                
                  1
                  
                    
                      L
                      ′
                    
                    s
                  
                
              
              
                
                  
                    -
                    
                      1
                      
                        
                          τ
                          ′
                        
                        σ
                      
                    
                  
                  -
                  
                    j
                    ⁢
                    
                      
                    
                    ⁢
                    
                      ω
                      s
                    
                  
                
              
              
                
                  
                    1
                    
                      
                        L
                        ′
                      
                      s
                    
                  
                  ⁢
                  
                    (
                    
                      
                        1
                        
                          τ
                          r
                        
                      
                      -
                      
                        j
                        ⁢
                        
                          
                        
                        ⁢
                        
                          
                            ω
                            ^
                          
                          m
                        
                      
                    
                    )
                  
                
              
            
            
              
                0
              
              
                0
              
              
                
                  R
                  R
                
              
              
                
                  
                    -
                    
                      1
                      
                        τ
                        r
                      
                    
                  
                  -
                  
                    j
                    ⁡
                    
                      (
                      
                        
                          ω
                          s
                        
                        -
                        
                          
                            ω
                            ^
                          
                          m
                        
                      
                      )
                    
                  
                
              
            
          
          ⁢
          
            
          
          ]
        
      
    
    
      
        (
        10
        )
      
    
  




K=[k1 k2 k3 k4]T  (11)

where the estimated states are marked by the symbol ‘^’.
The conventional speed adaptation law for the induction motor [7] is modified for the case where an LC filter is used. The estimation error of the inverter current is used for the speed adaptation, instead of the estimation error of the stator current as in the prior art systems. In order to stabilize the regeneration mode at low speeds, the idea of a rotated current estimation error [9], [11] is adopted.
The speed-adaptation law in the estimated rotor flux reference frame is

{circumflex over (ω)}m=−KpIm{(iA−îA)e−jφ}−Ki∫Im{(iA−îA)e−jφ}dt  (12)

where Kp and Ki are real adaptation gains, and the angle φ changes the direction of the error projection. The digital implementation of the adaptive full-order observer is based on a simple symmetric Euler method [12].
According to the method of the invention, the inverter output current vector iA and the inverter output voltage vector uA are determined. These determinations are normal current and voltage measurements. In a practice, the inverter output voltage vector uA is replaced in equation (9) by its reference uA,ref. Typically only two phase currents are measured to obtain the current vector. The output voltage of the inverter can be determined together from information of the states of the output switches and from voltage of the intermediate circuit.
Method of the invention further comprises a step of forming a full-order observer having a system matrix Â and gain vector K as explained above. The observer produces the estimated rotor flux linkage vector {circumflex over (ψ)}R, the estimated stator current vector îs, the estimated stator voltage vector ûs and the estimated inverter output current vector îA. These estimates can be used in the motor control in normal manner.
Since the inverter output current is both estimated and determined, the difference of them can be used in a speed adaptation loop, which produces an estimate for the electrical angular speed {circumflex over (ω)}m of the induction machine. Basically the speed adaptation law corrects the estimate of the angular speed so that the determined and estimated current vectors are similar.
Surprisingly the q components of estimated stator current vector and inverter output current vector are almost identical, which makes it possible to use the inverter output current instead of stator current in the speed adaptation as explained later.
With the method of the invention, all required information is gathered in order to control the induction machine.
Steady-State Analysis
The dynamics of the estimation error {tilde over (x)}=x−{circumflex over (x)} is obtained from (7) and (9):

{tilde over ({dot over (x)}=(A−KC){tilde over (x)}+(A−Â){circumflex over (x)}  (13)

where the difference between system matrices is




  
    
      
        
          
            A
            _
          
          -
          
            
              A
              _
            
            ^
          
        
        =
        
          
            
              [
              
                
              
              ⁢
              
                
                  
                    0
                  
                  
                    0
                  
                  
                    0
                  
                  
                    0
                  
                
                
                  
                    0
                  
                  
                    0
                  
                  
                    0
                  
                  
                    0
                  
                
                
                  
                    0
                  
                  
                    0
                  
                  
                    0
                  
                  
                    
                      
                        -
                        j
                      
                      ⁢
                      
                        1
                        
                          L
                          
                            ′
                            s
                          
                        
                      
                    
                  
                
                
                  
                    0
                  
                  
                    0
                  
                  
                    0
                  
                  
                    j
                  
                
              
              ⁢
              
                
              
              ]
            
            
              }
              
                M
                _
              
            
          
          ⁢
          
            (
            
              
                ω
                m
              
              -
              
                
                  ω
                  ^
                
                m
              
            
            )
          
        
      
    
    
      
        (
        14
        )
      
    
  




For the steady-state analysis, the derivative of estimation error (13) is set to zero. The operation point is determined by the synchronous angular frequency ωs, the slip angular frequency ωr=ωs−ωm, and the estimated rotor flux {circumflex over (ψ)}R. The example values of a 2.2-kW four-pole induction motor (400 V, 50 Hz), shown in Table I, were used for the following analysis.
FIG. 3 illustrates the current estimation error as the slip is varied from the negative rated slip to the positive rated slip, for various values of the rotor speed estimation error {tilde over (ω)}m=ωm−{circumflex over (ω)}m between −0.005 p.u. and 0.005 p.u. The estimated rotor flux is constant and the synchronous frequency ωs is 0.1 p.u. in FIG. 3a and 0.8 p.u. in FIG. 3b. The observer gain is K=[3000s−1 0 0 0]T, and the base value of the angular frequency is 2π50 rad/s, and the rated slip frequency is ωrN=0.05 p.u.
If the angle correction is not used in the adaptation law (12), the rotor speed estimate is calculated according to the imaginary part of the current estimation error. This kind of adaptation law works well in the motoring mode (where ωsωr>0), but at low synchronous speeds in the regeneration mode (ωsωr<0), the imaginary part of the current estimation error changes its sign at a certain slip angular frequency as can be seen in FIG. 3(a). Beyond this point, the estimated rotor speed is corrected to the wrong direction, leading to unstable operation. This problem is not encountered at high synchronous speeds, as can seen in FIG. 3(b).
The unstable operation can be avoided if the real part of the current estimation error is also taken into account in the speed adaptation. Correspondingly, the current estimation error is rotated by a factor e−jφ. The angle φ is selected as [9]




  
    
      
        φ
        =
        
          {
          
            
              
                
                  
                    
                      φ
                      max
                    
                    ⁢
                    
                      
                    
                    ⁢
                    sign
                    ⁢
                    
                      
                    
                    ⁢
                    
                      (
                      
                        ω
                        s
                      
                      )
                    
                    ⁢
                    
                      (
                      
                        1
                        -
                        
                          
                            
                            
                              ω
                              s
                            
                            
                          
                          
                            ω
                            φ
                          
                        
                      
                      )
                    
                    ⁢
                    
                      
                    
                    ⁢
                    if
                    ⁢
                    
                      
                    
                    ⁢
                    
                      
                      
                        ω
                        s
                      
                      
                    
                  
                  <
                  
                    
                      ω
                      φ
                    
                    ⁢
                    
                      
                    
                    ⁢
                    and
                    ⁢
                    
                      
                    
                    ⁢
                    
                      ω
                      s
                    
                    ⁢
                    
                      
                        ω
                        ^
                      
                      r
                    
                  
                  <
                  0
                
              
            
            
              
                
                  
                    0
                    ⁢
                    
                      
                    
                    ⁢
                    otherwise
                  
                  ⁢
                  
                    
                  
                
              
            
          
        
      
    
    
      
        (
        15
        )
      
    
  




where φmax is the maximum correction angle and ωφ is the limit for the synchronous angular frequency after which the correction is not used.
The influence of the angle correction is shown in FIG. 4. The parameter values φmax=0.467 and ωφ=0.75 p.u. are used in this example. In the regeneration-mode operation, the current error vector is rotated clockwise and the estimated rotor speed is corrected to the right direction.
FIG. 4 shows loci of current estimation error as slip varies from negative rated slip to positive rated slip for various rotor speed estimation error values between −0.005 p.u. and 0.005 p.u. The synchronous angular frequency 0.1 p.u. The current error is rotated at negative slip angular frequencies (dashed curves).
Dynamic Analysis
The dynamic behavior of the speed-adaptive observer can be analyzed via linearization. The operating point is set by the equilibrium quantities: the rotor angular speed ωm0 and the synchronous angular frequency ωs0, and the rotor flux ψR0. The linearized estimation error is

{tilde over ({dot over (x)}=(A0−KC){tilde over (x)}+Mx0(ωm−{circumflex over (ω)}m)  (16)

The transfer function from the speed estimation error to the inverter current estimation error obtained from (16) is




  
    
      
        
          
            G
            _
          
          ⁡
          
            (
            s
            )
          
        
        =
        
          
            
              
                
                  
                    i
                    _
                  
                  A
                
                ⁡
                
                  (
                  s
                  )
                
              
              -
              
                
                  
                    
                      i
                      ^
                    
                    _
                  
                  A
                
                ⁡
                
                  (
                  s
                  )
                
              
            
            
              
                
                  ω
                  m
                
                ⁡
                
                  (
                  s
                  )
                
              
              -
              
                
                  
                    ω
                    ^
                  
                  m
                
                ⁡
                
                  (
                  s
                  )
                
              
            
          
          =
          
            
              
                C
                ⁡
                
                  (
                  
                    sI
                    -
                    
                      
                        A
                        _
                      
                      0
                    
                    +
                    
                      
                        K
                        _
                      
                      ⁢
                      C
                    
                  
                  )
                
              
              
                -
                l
              
            
            ⁢
            
              
                Mx
                _
              
              o
            
          
        
      
    
    
      
        (
        17
        )
      
    
  




Based on (12), the transfer function from the imaginary part of the rotated inverter current error, Im{(iA−îA)e−jφ}, to the speed estimate, {circumflex over (ω)}m, is




  
    
      
        
          K
          ⁡
          
            (
            s
            )
          
        
        =
        
          
            -
            
              K
              p
            
          
          -
          
            
              K
              i
            
            s
          
        
      
    
    
      
        (
        18
        )
      
    
  




The resulting linearized system model for dynamic analysis is shown in FIG. 5. This model is used for investigating the pole locations of the linearized system at different operating points.
The observer gain K affects the stability of the system. In an induction motor drive without an LC filter, the adaptive observer is stable in the motoring mode even with zero gain [8]. However, zero gain cannot be used when an LC filter is present. FIG. 6(a) shows the observer poles as the synchronous angular frequency ωs is varied from −1 p.u. to 1 p.u. and the slip is rated. The adaptive observer with zero gain is unstable in the motoring mode, corresponding to the poles in the right half-plane.
In order to obtain a simple observer structure, the observer gain

K=[k, 0 0 0]T  (19)

is proposed. The poles obtained are shown in FIG. 6(b). All poles stay in the left half-plane in the whole inspected operation region.
Simulation Results
The system was investigated by computer simulation with Matlab/Simulink software. The data of a 2.2-kW induction motor, given in Table I, were used for the simulations. The LC filter was designed according to the design rules in [13], [14]. The sampling frequency was equal to the switching frequency of 5 kHz. The bandwidths of the controllers were 500 Hz for the inverter current, 250 Hz for the stator voltage, 150 Hz for the stator current, 15 Hz for the rotor speed, and 1.5 Hz for the rotor flux. The speed estimate was filtered using a low-pass filter having the bandwidth of 40 Hz. The reference voltage uA,ref was used in the observer instead of the actual inverter output voltage uA. The observer gain was K=[3000s−1 0 0 0]T, and the adaptation gains were chosen as Kp=10(As)−1 and Ki=20000 (As2)−1.
FIG. 7 shows simulation results obtained for a sequence consisting of a fast acceleration from zero speed to 0.8 p.u., a rated load step, a slow speed reversal, and a stepwise load reversal to the negative rated load. The motor was in the regenerating mode between t=7.6 s and t=12 s. During the rest of the sequence, the motor was in the motoring mode with the exception of a short time in the plugging mode during the speed reversal. The first subplot shows the rotor speed (solid) and its estimate (dashed). The second subplot shows the q component of the inverter current (solid) and its estimate (dashed). The third subplot shows the q component of the stator current (solid) and its estimate (dashed). It should be noted, that the solid and dashed plots are similar, meaning that the estimates are accurate.
Slow speed reversals at load are difficult for sensorless induction motor drives. Although the exact motor and filter parameters are used in the simulation, small vibrations appear in iAq and isq at t=7.8 s when the synchronous frequency is zero. More problems are encountered if slower speed reversals are needed or parameter estimates are inaccurate. It is to be noted that the q components of the inverter and stator currents are nearly equal, which makes it possible to use iAq in the speed adaptation law instead of isq as in the prior art systems. The voltage and current waveforms before and after the LC filter are illustrated in detail in FIG. 8. The stator voltage and current are nearly sinusoidal. The first subplot shows the inverter output voltage (phase-to-phase) and the stator voltage (phase-to-phase). The second subplot shows the inverter current and the stator current.
EXPERIMENTAL RESULTS
The experimental setup is illustrated in FIG. 9. The 2.2-kW four-pole induction motor was fed by a frequency converter controlled by a dSPACE DS1103 PPC/DSP board. The parameters of the experimental setup correspond to those given in Table I. In the LC filter three 3.3-μF capacitors were used in delta connection, giving the per-phase capacitance value of 9.9 μF. The dc-link voltage was measured, and the reference voltage obtained from the inverter current controller was used for the observer. The rotor speed and the shaft torque were measured only for monitoring. A simple current feedforward compensation for dead times and power device voltage drops was applied [15]. A permanent magnet servo motor was used as a load machine.
FIG. 10 presents the experimental results corresponding to the simulations shown in FIG. 7. The measured performance corresponds well to the simulation results, but the vibrations in iAq and isq are pronounced at zero synchronous frequency. The vibrations are due to inaccurate motor parameter estimates and inverter nonidealities. The explanations of FIG. 10 are as in FIG. 7.
The vibration can be attenuated by lowering the observer gain at low speeds. The observer gain is selected as




  
    
      
        
          k
          1
        
        =
        
          {
          
            
              
                
                  
                    
                      k
                      
                        1
                        ⁢
                        I
                      
                    
                    +
                    
                      
                        
                          
                            k
                            
                              1
                              ⁢
                              h
                            
                          
                          -
                          
                            k
                            
                              1
                              ⁢
                              I
                            
                          
                        
                        
                          ω
                          d
                        
                      
                      ⁢
                      
                        
                      
                      ⁢
                      
                        
                        
                          
                            ω
                            ^
                          
                          m
                        
                        
                      
                      ⁢
                      
                        
                      
                      ⁢
                      if
                      ⁢
                      
                        
                      
                      ⁢
                      
                        
                        
                          
                            ω
                            ^
                          
                          m
                        
                        
                      
                    
                  
                  <
                  
                    ω
                    d
                  
                
              
            
            
              
                
                  
                    
                  
                  ⁢
                  
                    
                      k
                      
                        1
                        ⁢
                        h
                      
                    
                    ⁢
                    
                      
                    
                    ⁢
                    otherwise
                  
                  ⁢
                  
                    
                  
                
              
            
          
        
      
    
    
      
        (
        20
        )
      
    
  




where k1l and k1h are the minimum and the maximum gains. The speed limit after which the maximum gain is used is ωd. The limit angular frequencies ωφ for the angle correction in (15) must also be changed because of the lower gain at low speeds. In order to avoid stepwise changes, the limit is selected as




  
    
      
        
          ω
          φ
        
        =
        
          {
          
            
              
                
                  
                    
                      ω
                      
                        φ
                        ⁢
                        
                          
                        
                        ⁢
                        l
                      
                    
                    +
                    
                      
                        
                          
                            ω
                            
                              φ
                              ⁢
                              
                                
                              
                              ⁢
                              h
                            
                          
                          -
                          
                            ω
                            
                              φ
                              ⁢
                              
                                
                              
                              ⁢
                              l
                            
                          
                        
                        
                          ω
                          d
                        
                      
                      ⁢
                      
                        
                        
                          
                            ω
                            ^
                          
                          m
                        
                        
                      
                      ⁢
                      
                        
                      
                      ⁢
                      if
                      ⁢
                      
                        
                      
                      ⁢
                      
                        
                        
                          
                            ω
                            ^
                          
                          m
                        
                        
                      
                    
                  
                  <
                  
                    ω
                    d
                  
                
              
            
            
              
                
                  
                    ω
                    
                      φ
                      ⁢
                      
                        
                      
                      ⁢
                      h
                    
                  
                  ⁢
                  
                    
                  
                  ⁢
                  otherwise
                
              
            
          
        
      
    
    
      
        (
        21
        )
      
    
  




where ωφl and ωφh are the limit angular frequencies for the angle correction at low and high speeds, respectively.

Experimental results obtained for the modified observer gain and angle correction are shown in FIG. 11. The parameters in (20) and (21) were k11=1000 s−1, k1h=3000 s−1, ωφl=0.4 p.u., ωφh=0.75 p.u., and ωd=0.38 p.u. The explanations of FIG. 11 are as in FIG. 7.
The vibrations at the zero synchronous speed are reduced significantly, but not totally removed.
It should be noted, that the described control system is only one possible system for controlling an induction machine based on the method of the invention. It will be obvious to a person skilled in the art that, as the technology advances, the inventive concept can be implemented in various ways. The invention and its embodiments are not limited to the examples described above but may vary within the scope of the claims.







TABLE I





PARAMETERS OF THE MOTOR AND THE LC FILTER


















Motor Parameters




Stator resistance RS
3.67 Ω



Rotor resistance RR
1.65 Ω



Stator transient inductance L′S
0.0209 H



Magnetizing inductance LM
0.264 H



Total moment of inertia J
0.0155 kgm2



Rated speed nN
1430 r/min



Rated (base) current IN
5.0 A



Rated torque TN
14.6 Nm



LC Filter Parameters



Inductance Lf
8 mH



Capacitance Cf
9.9 μF



Series resistance RLf
0.1 Ω










REFERENCES


    [1] D. F. Busse, J. M. Erdman, R. J. Kerkman, D. W. Schlegel, and G. L. Skibinski, “The effects of PWM voltage source inverters on the mechanical performance of rolling bearings,” IEEE Trans. Ind. Applicat., vol. 33, March/April 1997, pp. 567–576.
    [2] A. von Jouanne, P. Enjeti, and W. Gray, “The effect of long motor leads on PWM inverter fed AC motor drive systems,” in Proc. IEEE APEC, vol. 2, Dallas, Tex., March 1995, pp. 592–597.
    [3] 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, January/February 2004, pp. 162–169.
    [4] A. Nabae, H. Nakano, and Y. Okamura, “A novel control strategy of the inverter with sinusoidal voltage and current outputs,” in Proc. IEEE PESC'94, vol. 1, Taipei, Taiwan, June 1994, pp. 154–159.
    [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 2001, Graz, Austria, August 2001, CD-ROM.
    [6] J. Salomaki and J. Luomi, “Vector control of an induction motor fed by a PWM inverter with output LC filter,” in Proc. NORPIE/12004, Trondheim, Norway, June 2004, CD-ROM.
    [7] H. Kubota, K. Matsuse, and T. Nakano, “DSP-based speed adaptive flux observer of induction motor,” IEEE Trans. Ind. Applicat., vol. 29, March/April 1993, pp. 344–348.
    [8] H. Kubota, I. Sato, Y. Tamura, K. Matsuse, H. Ohta and Y. Hori, “Regenerating-mode low-speed operation of sensorless induction motor drive with adaptive observer,” IEEE Trans. Ind. Applicat., vol. 38, July/August 2002, pp. 1081–1086.
    [9] M. Hinkkanen and J. Luomi, “Stabilization of the regenerating mode of full-order flux observers for sensorless induction motors,” in Proc. IEEE IEMDC'03, Madison, Wis., June 2003.
    [10] G. R. Slemon, “Modelling of induction machines for electric drives,” IEEE Trans. Ind. Applicat., vol. 25, November/December 1989, pp. 1126–1131.
    [11] M. Hasegawa and K. Matsui, “Robust adaptive full-order observer design with novel adaptive scheme for speed sensorless vector controlled induction motors,” in Proc. IEEE IECON '02, Sevilla, Spain, November 2002.
    [12] J. Niiranen, “Fast and accurate symmetric Euler algorithm for electromechanical simulations,” in Proc. Electrimacs'99, Lisboa, Portugal, September 1999, pp. 71–78.
    [13] J. Steinke, C. Stulz, and P. Pohjalainen, “Use of a LC filter to achieve a motor friendly performance of the PWM voltage source inverter,” in Proc. IEEE IEMDC'97, Milwaukee, Wis., May 1997, pp. TA2/4.1–TA2/4.3.
    [14] C. Xiyou, Y. Bin, and G. Yu, “The engineering design and the optimization of inverter output RLC filter in AC motor drive system,” in Proc. IEEE IECON'02, vol. 1, Sevilla, Spain, November 2002, pp. 175–180.
    [15] J. K. Pedersen, F. Blaabjerg, J. W. Jensen, and P. Thogersen, “An ideal PWM-VSI inverter with feedforward and feedback compensation,” in Proc. EPE'93, vol. 4, Brighton, U.K., September 1993, pp. 312–318.

## Claims

1. A method of controlling of an induction machine using an inverter with an output LC filter, comprising:
determining an inverter output current vector (iA)
determining an inverter output voltage vector (uA),
forming a full-order observer having a system matrix (Â) and gain vector (K), the observer producing an estimated rotor flux linkage vector ({circumflex over (ψ)}R), an estimated stator current vector (îs), an estimated stator voltage vector (û{circumflex over (us)}) and an estimated inverter output current vector (îA)
determining an estimation error (iA−îA) of the inverter output current vector,
forming a speed adaptation law based on the estimation error of the inverter output current vector for determining an estimate for electrical angular speed ({circumflex over (ω)}m) of the induction machine, and
controlling the induction machine based on at least one of the estimated rotor flux linkage vector, estimated stator current vector, estimated stator voltage vector and estimated inverter output current vector and a measured inverter output current.

2. A method according to claim 1, comprising:
defining a correction angle (φ) and
changing a direction of the estimation error (iA−îA) according to the correction angle (φ) before forming the speed adaptation law to keep the observer stable in a regenerating mode.

3. A method according to claim 1, wherein determining of the inverter output voltage vector (uA) comprises:
using an inverter output voltage reference uA,ref.

4. A method according to claim 1, wherein the controlling is based on the estimated rotor flux linkage vector, estimated stator current vector, estimated stator voltage vector and estimated inverter output current vector.