# Patent

## Title
Motor control device

## Abstract

A motor control device that includes a first speed estimator estimating the rotation speed of the rotor of a permanent-magnet synchronous motor and that controls the motor so that a first estimated rotation speed estimated by the first speed estimator follows the specified speed value further includes a second speed estimator that estimates the rotation speed of the rotor by an estimation method different from that used by the first speed estimator. The motor control device detects synchronization failure based on a comparison between a second estimated rotation speed estimated by the second speed estimator and the first estimated rotation speed or the specified speed value.

## Background
This nonprovisional application claims priority under 35 U.S.C. §119(a) on Patent Application No. 2006-106293 filed in Japan on Apr. 7, 2006, the entire contents of which are hereby incorporated by reference.
BACKGROUND OF THE INVENTION
1. Field of the Invention
The present invention relates to a motor control device for controlling the operation of a motor, and more particularly to a motor control device with a capability of detecting synchronization failure. The present invention also relates to a motor drive system incorporating such a motor control device.
2. Description of Related Art
There have conventionally been developed motor control devices (position-sensorless control devices) that estimate the rotor position of a motor without the use of a rotor position sensor and controls the motor based on the so estimated rotor position. In motor control devices of this type, an abrupt variation in the load of the motor may cause synchronization failure. If synchronization failure occurs, the motor may come to a standstill and become uncontrollable. Thus, techniques that can surely detect synchronization failure have been sought.
There have been proposed various techniques for detecting such synchronization failure. Many of them detect synchronization failure by paying attention to a lowering in power factor, or ineffective current, resulting from a current that does not flow in normal operation.
For example, JP-A-2001-025282 discloses a technique whereby synchronization failure is detected through comparison of the period of ineffective current and the period of its voltage in case of a standstill resulting from synchronization failure.
For another example, a technique is known whereby synchronization failure is detected when the deviation between the estimated speed as estimated by a speed estimator and the specified speed has remained equal to or greater than a predetermined value for a predetermined period or longer. These conventional techniques, however, have the disadvantage that, in particular conditions, they cannot detect synchronization failure satisfactorily (for example, they simply cannot detect synchronization failure at all).
Incidentally, JP-A-H11-018499 discloses a technique whereby synchronization failure is detected by estimating and constantly comparing the γ-axis and δ-axis induction voltages that are generated as the rotor of a motor rotates. On the other hand, JP-A-2004-064902 discloses a technique whereby, to detect synchronization failure of a synchronous motor, the voltage to be specified is determined based on the motor generation coefficient and the specified excitation current, then a voltage correction amount is produced to make the specified excitation current and the actual excitation current equal, and then the voltage correction amount is compared with a predetermined value.
Thus, while, as described above, techniques that can surely detect synchronization failure have been sought in motor control devices that control a motor without the use of a rotor position sensor, none of the conventionally known techniques can satisfactorily meet the requirements.
SUMMARY OF THE INVENTION
According to one aspect of the present invention, a motor control device that includes a first speed estimator estimating the rotation speed of the rotor of a permanent-magnet synchronous motor and that controls the motor so that a first estimated rotation speed estimated by the first speed estimator follows a specified speed value further includes: a second speed estimator estimating the rotation speed of the rotor by an estimation method different from that used by the first speed estimator; and a synchronization failure detector detecting synchronization failure based on a second estimated rotation speed estimated by the second speed estimator and on the first estimated rotation speed or the specified speed value.
Specifically, for example, the synchronization failure detector may detect synchronization failure based on the ratio or deviation between the second estimated rotation speed and the first estimated rotation speed or the specified speed value.
According to another aspect of the present invention, a motor control device that includes a first speed estimator estimating the rotation speed of the rotor of a permanent-magnet synchronous motor and that controls the motor so that a first estimated rotation speed estimated by the first speed estimator follows a specified speed value further includes: a synchronization failure detector estimating, based on one or two of a first, a second, and a third state quantity, another of the first, second, and third state quantities and detecting synchronization failure based on the result of the estimation, the first state quantity being related to the motor voltage applied to the motor, the second state quantity being related to the rotation speed of the rotor, and the third state quantity being related to a motor current supplied to the motor.
Specifically, for example, let the axis parallel to the magnetic flux produced by the permanent magnet provided on the rotor be the d-axis, let the axis estimated, for control purposes, to correspond to the d-axis be the γ-axis, and let the axis estimated, for control purposes, to be perpendicular to the γ-axis be the δ-axis, then the second state quantity may be represented by the first estimated rotation speed or the specified speed value, the synchronization failure detector may estimate as the first state quantity the δ-axis component of the motor voltage based on the second state quantity alone or the second state quantity and the third state quantity represented by the measured current of the motor current or the specified current value to be followed by the measured current, and the synchronization failure detector detects synchronization failure based on the thus estimated δ-axis voltage as estimated as the first state quantity.
Moreover, for example, the synchronization failure detector may detect synchronization failure based on the ratio or deviation between the estimated δ-axis voltage and the specified δ-axis voltage value to be followed by the δ-axis component of the motor voltage.
Specifically, for example, let the axis parallel to the magnetic flux produced by the permanent magnet provided on the rotor be the d-axis, let the axis estimated, for control purposes, to correspond to the d-axis be the γ-axis, and let the axis estimated, for control purposes, to be perpendicular to the γ-axis be the δ-axis, then the first state quantity may be represented by the specified δ-axis voltage value to be followed by the δ-axis component of the motor voltage, the second state quantity may be represented by the first estimated rotation speed or the specified speed value, and the synchronization failure detector estimates as the third state quantity the δ-axis component of the motor current based on the first and second state quantities, and detects synchronization failure based on the thus estimated δ-axis current as estimated as the third state quantity.
Moreover, for example, the synchronization failure detector may detect the synchronization failure based on the ratio or deviation between the estimated δ-axis current and the δ-axis component of the measured current or the specified δ-axis current value to be followed by the δ-axis component of the measured current, or based on a comparison between the estimated δ-axis current and a predetermined value.
Specifically, for example, let the axis parallel to the magnetic flux produced by the permanent magnet provided on the rotor be the d-axis, let the axis estimated, for control purposes, to correspond to the d-axis be the γ-axis, and let the axis estimated, for control purposes, to be perpendicular to the γ-axis be the δ-axis, then the first state quantity is represented by the specified δ-axis voltage value to be followed by the δ-axis component of the motor voltage, the synchronization failure detector estimates as the second state quantity a second estimated rotation speed, separate from the first estimated rotation speed, based on the first state quantity alone or the first state quantity and the third state quantity represented by the measured current of the motor current or the specified current value to be followed by the measured current, and the synchronization failure detector detects synchronization failure based on the second estimated rotation speed.
According to still another aspect of the present invention, a motor drive system includes: a motor; an inverter driving the motor; and one of the motor control devices described above that controls the motor by controlling the inverter.



BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 is a block diagram showing an outline of the configuration of a motor drive system as a first embodiment of the present invention;
FIG. 2 is an analysis model diagram of the motor in first embodiment of the present invention;
FIG. 3 is a configuration block diagram of the motor drive system shown in FIG. 1;
FIG. 4 is a configuration block diagram of a motor drive system as a second embodiment of the present invention;
FIG. 5 is a configuration block diagram of a motor drive system as a third embodiment of the present invention;
FIG. 6 is a configuration block diagram of a motor drive system as a fourth embodiment of the present invention;
FIG. 7 is a diagram showing an example of the voltage vector locus of the superposition voltage outputted from the superposition voltage generator shown in FIG. 6;
FIG. 8 is a diagram showing the current vector locus of the superposition current that flows according to the superposition voltage shown in FIG. 7;
FIG. 9 is a waveform diagram showing the arithmetic product of the γ-axis and δ-axis components of the superposition current along with the direct-current component of the arithmetic product (assuming that the axis error is zero);
FIG. 10 is a waveform diagram showing the arithmetic product of the γ-axis and δ-axis components of the superposition current along with the direct-current component of the arithmetic product (assuming that the axis error is non-zero);
FIG. 11 is an internal block diagram of the second speed estimator shown in FIG. 6;
FIG. 12 is an analysis model diagram of the motor, illustrating a dmqm method applicable in any embodiment of the present invention;
FIG. 13 is an analysis model diagram of the motor, illustrating a dmqm method applicable in any embodiment of the present invention;
FIG. 14 is a diagram showing the current locus of the motor current, as observed when a dmqm method is adopted;
FIG. 15 is an internal block diagram of the first speed estimator shown in FIGS. 3 to 6, as observed when a dmqm method is adopted;
FIG. 16 is a diagram showing the relationship among motor parameters, as observed when a dmqm method is adopted;
FIG. 17 is a diagram illustrating the maximum torque control realized when a dmqm method is adopted; and
FIG. 18 is a vector diagram obtained in a synchronized state (in synchronized operation) when, for example, a dmqm method is adopted.



DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS
Hereinafter, the present invention will be described by way of specific embodiments with reference to the accompanying drawings. Among these drawings, such parts that appear more than once are identified with common reference designations, and in principle the description of any part, once given, will not be repeated. Likewise, throughout the drawings, the same symbols (such as θ and ω) are stuck to represent the same quantities. Moreover, to make the description simple, state quantities and the like are often referred to by their symbols alone; for example, the “estimated motor speed ωe” is sometimes referred to simply by “ωe”, and these should be understood to represent the same thing.
First Embodiment
First, a first embodiment of the present invention will be described. FIG. 1 is a block configuration diagram of a motor drive system according to the first embodiment. The reference numeral 1 represents a three-phase permanent-magnet synchronous motor (hereinafter simply the “motor 1”) having a permanent magnet on a rotor (unillustrated) and having an armature winding on a stator (unillustrated). In the following description, what are referred to simply by the “armature winding” and the “rotor” should be understood to mean the armature winding provided on the stator of the motor 1 and the rotor of the motor 1 respectively. The motor 1 is, for example, a salient-pole motor as exemplified by a interior permanent magnet synchronous motor; it may instead be a non-salient-pole motor.
The reference numeral 2 represents a PWM (pulse width modulation) inverter that, according to the rotor position of the motor 1, supplies three-phase alternating-current voltages, consisting of U-phase, V-phase, and W-phase voltages, to the armature winding of the motor 1. These voltages supplied to the armature winding of the motor 1 are collectively referred to as the motor voltage (armature voltage) Va, and the currents supplied from the inverter 2 to the armature winding of the motor 1 are collectively referred to as the motor current (armature current) Ia.
The reference numeral 3 represents a motor control device (position-sensorless control device) that, based on the motor current Ia, estimates the rotor position or the like of the motor 1, and feeds the PWM inverter 2 with a signal for rotating the motor 1 at a desired rotation speed. This desired rotation speed is fed, in the form of a specified motor speed value ω*, from an unillustrated CPU (central processing unit) or the like to a motor control device (in this embodiment, the motor control device 3).
FIG. 2 is an analysis model diagram of the motor 1. FIG. 2 shows the U-phase, V-phase, and W-phase armature winding fixed axes. The reference numeral 1a represents the permanent magnet constituting the rotor of the motor 1. In a rotating coordinate system that rotates at the same speed as the magnetic flux produced by the permanent magnet 1a, the direction of the magnetic flux produced by the permanent magnet 1a is referred to as the d-axis, and the axis estimated, for control purposes, as corresponding to the d-axis is referred to as the γ-axis. Moreover, although unillustrated, the axis having a phase leading the d-axis by an electrical angle of 90 degrees is referred to as the q-axis, and the axis estimated to have a phase leading the γ-axis by an electrical angle of 90 degrees is referred to as the δ-axis. The coordinate axes of the rotating coordinate system having the d- and q-axes as its coordinate axes are called the d-q axes (real axes). The rotating coordinate system (estimated rotating coordinate system) estimated for control purposes has the γ- and δ-axes as its coordinate axes, and its coordinate axes are called the γ-δ axes.
The d-q axes rotates, and its rotation speed (i.e., the rotation speed of the rotor of the motor 1) is called the real motor speed ω. The γ-δ axes also rotates, and its rotation speed is called the estimated motor speed ωe. With respect to the d-q axes in rotation, the phase of the d-axis at a given moment is represented, relative to the U-phase armature winding fixed axis, by θ (real rotor position θ). Likewise, with respect to the γ-δ axes in rotation, the phase of the γ-axis at that given moment is represented, relative to the U-phase armature winding fixed axis, by θe (estimated rotor position θe). Then, the axis error Δθ between the d- and γ-axes (the axis error Δθ between the d-q axes and the γ-δ axes) is given by Δθ=θ−θe. The values of ω*, ω, ωe, and ω2e (described later) are given in electrical angular velocity.
In the following description, the γ-axis, δ-axis, d-axis, and q-axis components of the motor voltage Va are referred to as the γ-axis voltage vγ, δ-axis voltage vδ, d-axis voltage vd, and q-axis voltage vq respectively; the γ-axis, δ-axis, d-axis, and q-axis components of the motor current Ia are referred to as the γ-axis current iγ, δ-axis current iδ, d-axis current id, and q-axis current iq respectively.
Moreover, in the following description, Ra represents the motor resistance (the resistance of the armature winding of the motor 1); Ld and Lq represent the d-axis and q-axis inductances (the d-axis and q-axis components of the inductance of the armature winding of the motor 1); Φa represents the armature flux linkage ascribable to the permanent magnet 1a. The values of Ld, Lq, Ra, and Φa are determined at the time of the fabrication of the motor drive system, and are used by the motor control device to perform calculations. In the formulae presented in the course of the following description, s represents the Laplace operator.
FIG. 3 is a configuration block diagram of the motor drive system, and shows the internal configuration of the motor control device 3 shown in FIG. 1. The motor control device 3 includes a current detector 11, a coordinate converter 12, a subtracter 13, a subtracter 14, a current controller 15, a magnetic flux controller 16, a speed controller 17, a coordinate converter 18, a subtracter 19, a first speed estimator 20, an integrator 21, a second speed estimator 22, and a synchronization failure checker 23. The different parts constituting the motor control device (3, and also 3a, 3b, and 3c described later) can freely use all the values produced within the motor control device (3, and also 3a, 3b, and 3c described later).
The current detector 11 is realized with, for example, Hall devices or the like, and detects the U-phase current (the current flowing through the U-phase part of the armature winding) iu and the V-phase current (the current flowing through the V-phase part of the armature winding) iv of the motor current Ia supplied from the PWM inverter 2 to the motor 1. The coordinate converter 12 receives the U-phase current iu and V-phase current iv detected by the current detector 11, and converts them into the γ-axis current iγ and δ-axis current iδ, based on the estimated rotor position θe fed from the integrator 21. The conversion here is performed according to formula (1) below.




  
    
      
        
          [
          
            
              
                
                  i
                  γ
                
              
            
            
              
                
                  i
                  δ
                
              
            
          
          ]
        
        =
        
          
            
              2
            
            ⁡
            
              [
              
                
                  
                    
                      sin
                      ⁡
                      
                        (
                        
                          
                            θ
                            e
                          
                          +
                          
                            π
                            /
                            3
                          
                        
                        )
                      
                    
                  
                  
                    
                      sin
                      ⁢
                      
                        
                      
                      ⁢
                      
                        θ
                        e
                      
                    
                  
                
                
                  
                    
                      cos
                      ⁡
                      
                        (
                        
                          
                            θ
                            e
                          
                          +
                          
                            π
                            /
                            3
                          
                        
                        )
                      
                    
                  
                  
                    
                      cos
                      ⁢
                      
                        
                      
                      ⁢
                      
                        θ
                        e
                      
                    
                  
                
              
              ]
            
          
          ⁡
          
            [
            
              
                
                  
                    i
                    u
                  
                
              
              
                
                  
                    i
                    v
                  
                
              
            
            ]
          
        
      
    
    
      
        (
        1
        )
      
    
  




The first speed estimator 20 estimates and outputs the estimated motor speed ωe (first estimated rotation speed). The method by which the estimated motor speed ωe is estimated will be described in detail later, after the description of the fourth embodiment. The integrator 21 integrates the estimated motor speed ωe estimated by the first speed estimator 20, and thereby calculates the estimated rotor position θe.
The subtracter 19 subtracts from the specified motor speed value ω* the estimated motor speed ωe fed from the first speed estimator 20, and outputs the result of the subtraction (speed error). The speed controller 17, based on the subtraction result (ω*−ωe) from the subtracter 19, produces the specified δ-axis current value iδ*. This specified δ-axis current value iδ* represents the current to be followed by the δ-axis current iδ, which is the δ-axis component of the motor current Ia. The magnetic flux controller 16 outputs the specified γ-axis current value iγ*; meanwhile, as necessary, the specified δ-axis current value iδ* and the estimated motor speed ωe are referred to. The specified γ-axis current value iγ* represents the current to be followed by the γ-axis current iγ, which is the γ-axis component of the motor current Ia.
The subtracter 13 subtracts from the specified γ-axis current value iγ* outputted from the magnetic flux controller 16 the γ-axis current iγ outputted from the coordinate converter 12, and thereby calculates the current error (iγ*−iγ). The subtracter 14 subtracts from the specified δ-axis current value iδ* outputted from the speed controller 17 the δ-axis current iδ outputted from the coordinate converter 12, and thereby calculates the current error (iδ*−iδ).
The current controller 15 receives the current errors calculated by the subtracters 13 and 14, the γ-axis current iγ and δ-axis current iδ from the coordinate converter 12, and the estimated motor speed ωe from the first speed estimator 20, and calculates and outputs the specified γ-axis voltage value vγ* to be followed by the γ-axis voltage vγ and the specified δ-axis voltage value vδ* to be followed by the δ-axis voltage vδ such that the γ-axis current iγ follows the specified γ-axis current value iγ* and in addition that the δ-axis current iδ follows the specified δ-axis current value iδ*.
Based on the estimated rotor position θe fed from the integrator 21, the coordinate converter 18 performs reverse conversion on the specified γ-axis voltage value vγ* and specified δ-axis voltage value vδ*, and thereby produces specified three-phase voltage values consisting of a specified U-phase voltage value vu*, a specified V-phase voltage value vv*, and a specified W-phase voltage value vw* representing the U-phase, V-phase, and W-phase components of the motor voltage Va. The coordinate converter 18 then outputs these specified three-phase voltage values to the PWM inverter 2. The reverse conversion here is performed according to formula (2) consisting of two equations below.




  
    
      
        
          
            [
            
              
                
                  
                    v
                    u
                    *
                  
                
              
              
                
                  
                    v
                    v
                    *
                  
                
              
            
            ]
          
          =
          
            
              
                
                  2
                  3
                
              
              ⁡
              
                [
                
                  
                    
                      
                        cos
                        ⁢
                        
                          
                        
                        ⁢
                        
                          θ
                          e
                        
                      
                    
                    
                      
                        
                          -
                          sin
                        
                        ⁢
                        
                          
                        
                        ⁢
                        
                          θ
                          e
                        
                      
                    
                  
                  
                    
                      
                        cos
                        ⁡
                        
                          (
                          
                            
                              θ
                              e
                            
                            -
                            
                              2
                              ⁢
                              
                                π
                                /
                                3
                              
                            
                          
                          )
                        
                      
                    
                    
                      
                        -
                        
                          sin
                          ⁡
                          
                            (
                            
                              
                                θ
                                e
                              
                              -
                              
                                2
                                ⁢
                                
                                  π
                                  /
                                  3
                                
                              
                            
                            )
                          
                        
                      
                    
                  
                
                ]
              
            
            ⁡
            
              [
              
                
                  
                    
                      v
                      γ
                      *
                    
                  
                
                
                  
                    
                      v
                      δ
                      *
                    
                  
                
              
              ]
            
          
        
        ⁢
        
          

        
        ⁢
        
          
            v
            w
            *
          
          =
          
            -
            
              (
              
                
                  v
                  u
                  *
                
                +
                
                  v
                  v
                  *
                
              
              )
            
          
        
      
    
    
      
        (
        2
        )
      
    
  




Based on the specified three-phase voltage values (vu*, vv* and vw*), which represent the voltages to be applied to the motor 1, the PWM inverter 2 produces pulse-width-modulated signals, and supplies the motor current Ia commensurate with the specified three-phase voltage values to the armature winding of the motor 1 to drive it.
The second speed estimator 22 estimates the rotation speed of the rotor by an estimation method different from the one by which the first speed estimator 20 estimates the rotation speed of the rotor. In other words, the second speed estimator 22 calculates (estimates) the rotation speed of the rotor according to a calculation formula different from the one according to which the first speed estimator 20 calculates (estimates) the rotation speed of the rotor. The rotation speed of the rotor estimated by the second speed estimator 22 is called the second estimated motor speed ω2e.
Formula (3) below is the voltage equation of the motor 1 (the equation of the q-axis voltage vq) as generally known. Here, p represents the differentiation operator. Ignoring the second term (a transient term) in the right side of formula (3) and then substituting vδ, iδ, and iγ on the γ-δ axes for vq, iq and id on the d-q axes gives formula (4) below, which is the equation of the δ-axis voltage vδ.

vq=Raiq+pLqiq+ωLdid+ωΦa  (3)

vδ=Raiδ+ωLdiγ+ωΦa  (4)

According to formula (4), the second speed estimator 22 calculates the second estimated motor speed ω2e. In practice, the second speed estimator 22 refers to the relevant values calculated within the motor control device 3, and calculates the second estimated motor speed ω2e according to formula (5a) below. Normally, the magnetic flux produced by the permanent magnet 1a is significantly higher than the magnetic flux produced by the γ-axis current iγ, such that Φa>>Ldiγ. In addition, at an appreciable rotation speed, the voltage drop across the motor resistance is significantly smaller than the δ-axis voltage. Hence, the second estimated motor speed ω2e may be calculated according to formula (5b) below, which is an approximation of formula (5a).




  
    
      
        
          ω
          
            2
            ⁢
            e
          
        
        =
        
          
            
              v
              δ
              *
            
            -
            
              
                R
                a
              
              ⁢
              
                i
                δ
              
            
          
          
            
              Φ
              a
            
            +
            
              
                L
                d
              
              ⁢
              
                i
                γ
              
            
          
        
      
    
    
      
        (
        
          5
          ⁢
          a
        
        )
      
    
  
  
    
      
        
          ω
          
            2
            ⁢
            e
          
        
        ≈
        
          
            v
            δ
            *
          
          
            Φ
            a
          
        
      
    
    
      
        (
        
          5
          ⁢
          b
        
        )
      
    
  




By referring to the second estimated motor speed ω2e calculated by the second speed estimator 22 and the specified motor speed value ω* (or, as will be described later, the estimated motor speed ωe), the synchronization failure checker 23 checks whether or not the rotor is in a state of synchronization failure. Here, a state of synchronization failure refers to one in which the rotor becomes unable to rotate according to the specified motor speed value ω* (for example, it comes to a standstill). Such a state occurs, for example, when an abrupt variation arises in the load while the rotor is rotating. Also referred to as a state of synchronization failure is one in which the rotor, although an attempt is being made to recover it from a standstill, cannot start rotating again for some reason (a so-called locked state).
If the rotor is not recognized to be in a state of synchronization failure, it is recognized to be operating normally, that is, in a synchronized state. In a synchronized state, the specified motor speed value ω* and the estimated motor speed ωe are equal (or substantially equal), and in addition the estimated motor speed ωe and the second estimated motor speed ω2e also are equal (or substantially equal).
For example, the synchronization failure checker 23 recognizes the rotor to be in a state of synchronization failure if formula (6) below holds, and recognizes it to be in a synchronized state if formula (6) does not hold. This is to say, with formula (6) rewritten in terms of formulae (5a) and (5b), synchronization failure is recognized if either formula (7a) or (7b) holds.




  
    
      
        
          
            ω
            *
          
          
            ω
            
              2
              ⁢
              e
            
          
        
        >
        
          k
          1
        
      
    
    
      
        (
        6
        )
      
    
  
  
    
      
        
          ω
          *
        
        >
        
          
            k
            1
          
          ⁢
          
            
              
                v
                δ
                *
              
              -
              
                
                  R
                  a
                
                ⁢
                
                  i
                  δ
                
              
            
            
              
                Φ
                a
              
              +
              
                
                  L
                  d
                
                ⁢
                
                  i
                  γ
                
              
            
          
        
      
    
    
      
        (
        
          7
          ⁢
          a
        
        )
      
    
  
  
    
      
        
          ω
          *
        
        >
        
          
            k
            1
          
          ⁢
          
            
              v
              δ
              *
            
            
              Φ
              a
            
          
        
      
    
    
      
        (
        
          7
          ⁢
          b
        
        )
      
    
  




Here, k1 has a value greater than 1 (for example, 2) that is determined previously through experiments or otherwise as necessary. Typically, k1 is a constant; however, it may instead be a variable whose value varies according to the relevant value (such as ω*, vδ*, or iδ) within the motor control device 3. When synchronization failure occurs, the rotor usually comes to a standstill. As a result, the specified δ-axis voltage value vδ*, which is commensurate with the voltage actually applied to the motor 1, becomes low, and accordingly the second estimated motor speed ω2e becomes low (see formula (5a) or (5b) above). Hence, according to formula (6) above, that is, according to formula (7a) or (7b), synchronization failure can be detected.
In formulae (5a) and (7a), iγ may be replaced with iγ*. Likewise, in formulae (5a) and (7a), iδ may be replaced with iδ*. This is because, since iγ and iδ are measured current values, they are equal (or substantially equal) to iγ* and iδ*, which are specified current values. In formulae (6), (7a), and (7b), ω* may be replaced with ωe. This is because ω*, which is a specified speed value, is equal (or substantially equal) to ωe. In formula (6), (7a), and (7b), the inequality sign “>” may be replaced with “≧”.
When formula (6) above is used, synchronization failure is detected based on the ratio of ω2e to ω* (or ωe). Alternatively, synchronization failure may be detected based on the deviation between ω2e and ω* (or ωe). Specifically, for example, instead of formula (6), formula (8) below may be used to detect synchronization failure. In that case, the rotor is recognized to be in a state of synchronization failure if formula (8) below holds, and otherwise it is recognized to be in a synchronized state.

ω2e−ω*>G1  (8)

In formula (8), G1 is a function of ω2e, ω*, or ωe, and this function is, for example, so determined that, when formula (6) holds, formula (8) also holds and in addition that, when formula (6) does not hold, formula (8) also does not hold. In this case, formula (6) is equivalent to formula (8), and hence the detection of synchronization failure using formula (6) is practically the same as the detection of synchronization failure using formula (8). The value of G1 is constantly calculated according to ω2e, ω*, or ωe, or is previously determined and is stored as table data in a memory (unillustrated).
Even when the two speed estimators adopt different estimation methods, in a synchronized state, the two methods both calculate speeds close to the real values. In a state of synchronization failure, however, neither of the methods can accurately calculate the real value, and thus the values they calculate are not close to each other. Hence, by comparing the values calculated by the two speed estimators, it is possible to detect synchronization failure with high accuracy.
The first speed estimator 20 is used to control the driving of the motor 1, and the driving of the motor 1 is so controlled that the value (ωe) estimated by the first speed estimator 20 follows the specified speed value (ω*). Thus, in case of synchronization failure, the value estimated by the first speed estimator 20 is not correct, and there is a comparatively large deviation between the estimated value and the specified speed value. This deviation is expected to be exploited for the detection of synchronization failure. Inconveniently, however, it occasionally happens that speed control effects (abnormal) energization such that the estimated value becomes equal to the specified speed value, resulting in stable control. This can be avoided by the use of the second speed estimator as in this embodiment, and it is thereby possible to surely detect synchronization failure.
Second Embodiment
Next, a second embodiment of the present invention will be described. FIG. 4 is a block configuration diagram of a motor drive system according to the second embodiment. The motor drive system shown in FIG. 4 includes a motor 1, an inverter 2, and a motor control device 3a. 
The motor control device 3a includes an applied voltage estimator 25 and a synchronization failure checker 23a instead of the second speed estimator 22 and the synchronization failure checker 23 provided in the motor control device 3 shown in FIG. 3; in other respects, the motor control device 3a and the motor drive system shown in FIG. 4 are the same as the motor control device 3 and the motor drive system shown in FIG. 3. The following description therefore places emphasis on the differences from the first embodiment, and no overlapping description of the common features will be repeated.
The applied voltage estimator 25 estimates (calculates) the δ-axis voltage vδ according to formula (4) above. The δ-axis voltage vδ estimated by the applied voltage estimator 25 is referred to as the estimated δ-axis voltage vδe. The calculated estimated δ-axis voltage vδe is fed to the synchronization failure checker 23a. 
In practice, the applied voltage estimator 25 refers to the relevant values calculated within the motor control device 3a, and calculates the estimated δ-axis voltage vδe according to formula (9a) below. Considering that iγ is usually very small, the estimated δ-axis voltage vδe may be calculated according to formula (9b), which is an approximation of formula (9a). At an appreciable rotation speed, the voltage drop across the motor resistance is significantly smaller than ω*Φa. Hence, the estimated δ-axis voltage vδe may be calculated according to formula (9c), which is an approximation of formula (9a).

vδe=Raiδ+ω*Ldiγ+ω*Φa  (9a)

vδe≈Raiδ+ω*Φa  (9b)

vδe≈ω*Φa  (9c)

By referring to the estimated δ-axis voltage vδe calculated by the applied voltage estimator 25 and the specified δ-axis voltage value vδ* calculated by the current controller 15, the synchronization failure checker 23a checks whether or not the rotor is in a state of synchronization failure. If the rotor is not recognized in a state of synchronization failure, it is recognized to be in a synchronized state.
For example, the synchronization failure checker 23a recognizes the rotor to be in a state of synchronization failure if formula (10) below holds, and recognizes it to be in a synchronized state if formula (10) does not hold. This is to say, with formula (10) rewritten in terms of formula (9a), (9b) and (9c), synchronization failure is recognized if formula (11a), (11b) or (11c) below holds. Incidentally, in a synchronized state, vδe and vδ* are equal (or substantially equal).




  
    
      
        
          
            v
            δ
            *
          
          
            v
            
              δ
              ⁢
              
                
              
              ⁢
              e
            
          
        
        <
        
          k
          2
        
      
    
    
      
        (
        10
        )
      
    
  




vδ*<k2(Raiδ+ω*Ldiγ+ω*Φa)  (11a)

vδ*<k2(Raiδ+ω*Φa)  (11b)

vδ*<k2ω*Φa  (11c)

Here, k2 has a value smaller than 1 (for example, 0.5) that is determined previously through experiments or otherwise as necessary. Typically, k2 is a constant; however, it may instead be a variable whose value varies according to the relevant value (such as ω*, vδ*, or iδ) within the motor control device 3a. When synchronization failure occurs, the rotor usually comes to a standstill. As a result, the specified δ-axis voltage value vδ*, which is commensurate with the voltage actually applied to the motor 1, becomes low; on the other hand, the estimated δ-axis voltage vδe, which is calculated on the assumption that the rotor is rotating according to ω*, has a value commensurate with ω*. Hence, according to formula (10) above, that is, according to formula (11a), (11b), or (11c), synchronization failure can be detected.
In formula (9a) and (11a), iδ may be replaced with iγ*. Likewise, in formula (9a), (9b), (11a), and (11b), iδ may be replaced with iδ*. In formula (9a), (9b), (9c), (11a), (11b), and (11c), ω* may be replaced with ωe. In formula (10), (11a), (11b), and (11c), the inequality sign “<” may be replaced with “≦”.
When formula (10) above is used, synchronization failure is detected based on the ratio of vδe to vδ*. Alternatively, synchronization failure may be detected based on the deviation between vδe and vδ*. Specifically, for example, instead of formula (10), formula (12) below may be used to detect synchronization failure. In that case, the rotor is recognized to be in a state of synchronization failure if formula (12) below holds, and otherwise it is recognized to be in a synchronized state.

vδe−vδ*>G2  (12)

In formula (12), G2 is a function of vδe or vδ*, and this function is, for example, so determined that, when formula (10) holds, formula (12) also holds and in addition that, when formula (10) does not hold, formula (12) also does not hold. In this case, formula (10) is equivalent to formula (12), and hence the detection of synchronization failure using formula (10) is practically the same as the detection of synchronization failure using formula (12). The value of G2 is constantly calculated according to vδe or vδ*, or is previously determined and is stored as table data in a memory (unillustrated).
In the first embodiment, the speed is estimated from the specified voltage value; conversely, in this embodiment, the applied voltage is estimated from the specified speed value or the estimated speed, and synchronization failure is detected based on a comparison between the estimated applied voltage and the specified voltage value.
Third Embodiment
Next, a third embodiment of the present invention will be described. FIG. 5 is a block configuration diagram of a motor drive system according to the third embodiment. The motor drive system shown in FIG. 5 includes a motor 1, an inverter 2, and a motor control device 3b. 
The motor control device 3b includes a supplied current estimator 26 and a synchronization failure checker 23b instead of the second speed estimator 22 and the synchronization failure checker 23 provided in the motor control device 3 shown in FIG. 3; in other respects, the motor control device 3b and the motor drive system shown in FIG. 5 are the same as the motor control device 3 and the motor drive system shown in FIG. 3. The following description therefore places emphasis on the differences from the first embodiment, and no overlapping description of the common features will be repeated.
The supplied current estimator 26 estimates (calculates) the δ-axis current iδ according to formula (4) above. The δ-axis current iδ estimated by the supplied current estimator 26 is referred to as the estimated δ-axis current iδe. The calculated estimated δ-axis current iδe is fed to the synchronization failure checker 23b. 
In practice, the supplied current estimator 26 refers to the relevant values calculated within the motor control device 3b, and calculates the estimated δ-axis current iδe according to formula (13a) below. Considering that iγ is usually very small, the estimated δ-axis current iδe may be calculated according to formula (13b) below, which is an approximation of formula (13a).




  
    
      
        
          i
          
            δ
            ⁢
            
              
            
            ⁢
            e
          
        
        =
        
          
            1
            
              R
              a
            
          
          ⁢
          
            (
            
              
                v
                δ
                *
              
              -
              
                (
                
                  
                    
                      ω
                      *
                    
                    ⁢
                    
                      L
                      d
                    
                    ⁢
                    
                      i
                      γ
                    
                  
                  +
                  
                    
                      ω
                      *
                    
                    ⁢
                    
                      Φ
                      a
                    
                  
                
                )
              
            
            )
          
        
      
    
    
      
        (
        
          13
          ⁢
          a
        
        )
      
    
  
  
    
      
        
          i
          
            δ
            ⁢
            
              
            
            ⁢
            e
          
        
        ≈
        
          
            1
            
              R
              a
            
          
          ⁢
          
            (
            
              
                v
                δ
                *
              
              -
              
                
                  ω
                  *
                
                ⁢
                
                  Φ
                  a
                
              
            
            )
          
        
      
    
    
      
        (
        
          13
          ⁢
          b
        
        )
      
    
  




By referring to the estimated δ-axis current iδe calculated by the supplied current estimator 26 and the specified δ-axis current value iδ* calculated by the speed controller 17 (or the δ-axis current iδ calculated by the coordinate converter 12), the synchronization failure checker 23b checks whether or not the rotor is in a state of synchronization failure. Here, the reference to iδ* (or iδ) may be omitted. If the rotor is not recognized in a state of synchronization failure, it is recognized to be in a synchronized state.
For example, the synchronization failure checker 23b recognizes the rotor to be in a state of synchronization failure if formula (14), (15), or (16) below holds, and recognizes it to be in a synchronized state if none of those formula holds. The value of iδe in formula (14), (15), and (16) is calculated according to formula (13a) or (13b) above. Incidentally, in a synchronized state, iδe and iδ* are equal (or substantially equal).




  
    
      
        
          
            i
            
              δ
              ⁢
              
                
              
              ⁢
              e
            
          
          
            i
            δ
            *
          
        
        <
        
          k
          3
        
      
    
    
      
        (
        14
        )
      
    
  




iδe−iδ*<k4  (15)

iδe<k5  (16)

Here, k3, k4, and k5 have values that are determined previously through experiments or otherwise as necessary. Typically, k3, k4, and k5 are constants; however, any of them may instead be a variable whose value varies according to the relevant value (such as ω*, vδ*, or iδ) within the motor control device 3b. When synchronization failure occurs, the rotor usually comes to a standstill. As a result, the specified δ-axis voltage value vδ*, which is commensurate with the voltage actually applied to the motor 1, becomes low, specifically smaller than ω*Φa. That is, in case of synchronization failure, iδe has a negative value. When the motor 1 is driven, iδ*>0 always holds.
Accordingly, k3 is given a value smaller than 1, for example, 0, 0.5, or a negative value close to 0; k4 is given a negative value; and k5 is set equal 0 or is given a negative value.
In formula (13a), iγ may be replace with iγ*. Likewise, in formulae (14) and (15), iδ* may be replaced with iδ. In formulae (13a) and (13b), ω* may be replaced with ωe. In formula (14), (15), and (16), the inequality sign “<” may be replaced with “≦”.
Fourth Embodiment
Next, a fourth embodiment of the present invention will be described. FIG. 6 is a block configuration diagram of a motor drive system according to the fourth embodiment. The motor drive system shown in FIG. 6 includes a motor 1, an inverter 2, and a motor control device 3c. 
As compared with the motor control device 3 shown in FIG. 3, the motor control device 3c additionally includes a superposition voltage generator 28, a LPF (low-pass filter) 29, and adder 30, and includes, instead of the second speed estimator 22 and the synchronization failure checker 23 provided in the motor control device 3 shown in FIG. 3, a second speed estimator 27 and a synchronization failure checker 23c. In other respects, the motor control device 3c and the motor drive system shown in FIG. 6 are the same as the motor control device 3 and the motor drive system shown in FIG. 3. The following description therefore places emphasis on the differences from the first embodiment, and no overlapping description of the common features will be repeated.
The superposition voltage generator 28 generates and outputs superposition voltages to be superposed on the specified U-phase, V-phase, and W-phase voltage values vu*, vv*, and vw*. These superposition voltages consist of a U-phase superposition voltage vhu for vu* (the U-phase component of the superposition voltages), a V-phase superposition voltage vhv for vv* (the V-phase component of the superposition voltages), and a W-phase superposition voltage vhw for vw* (the W-phase component of the superposition voltages). The adder 30 adds vhu, vhv and vhw respectively to vu*, vv* and vw* from the coordinate converter 18, and outputs the results of the addition to PWM inverter 2.
Thus, in this embodiment, specified three-phase voltage values consisting of (vu*+vhu), (vv*+vhv), and (vw*+vhw) are, as the voltages to be applied to the motor 1, fed to the PWM inverter 2. The PWM inverter 2 drives the motor 1 by supplying the armature winding of the motor 1 with the motor current Ia commensurate with the specified three-phase voltage values having the superposition voltages superposed thereon.
As described above, in this embodiment, superposition voltages are superposed on the drive voltages, represented by vu*, vv* and vw*, for driving the motor 1. As a result of the superposition of these superposition voltages, superposition currents commensurate with the superposition voltages are superposed on the drive currents, referred to as the specified γ-axis and δ-axis current values iγ* and iδ*, for driving the motor 1.
The superposition voltages generated by the superposition voltage generator 28 form a high-frequency rotating voltage. Here, “high-frequency” means that the frequency of the superposition voltages is significantly higher than the frequency of the drive voltages. Accordingly, the frequency of the superposition currents, which is superposed according to those superposition voltages, is significantly higher than the frequency of the drive currents. On the other hand, “rotating voltage” means that the voltage vector locus of the superposition voltages is circular on fixed coordinate axes.
Even when considered on the γ-δ axes (the γ-δ plane), the voltage vector locus of the superposition voltages generated by the superposition voltage generator 28 is circular, for example like the voltage vector locus 70 shown in FIG. 7. When the superposition voltages are three-phase balanced voltages, their voltage vector locus is, like the voltage vector locus 70, perfectly circular, lying on the γ-δ axes and centered at the origin. This rotating voltage (forming the superposition voltages) is a high-frequency voltage that is not synchronous with the motor 1, and thus its application to the motor 1 does not cause it to rotate.
When the motor 1 is an interior permanent magnet synchronous motor or the like and Ld<Lq holds, the current vector locus of the superposition currents that are caused to flow through the motor 1 by the superposition voltages having the voltage vector locus 70 is elliptic, like the current vector locus 71 shown in FIG. 8, lying on the γ-δ axes (the γ-δ plane), centered at the origin, having a major axis aligned with the γ-axis, and having a minor axis aligned with the δ-axis. Here, the current vector locus 71 is one observed when the axis error Δθ equals zero. When the axis error Δθ is non-zero, the current vector locus of the superposition currents is elliptic, like the current vector locus 72, having a major (or minor) axis departed from the γ-axis (or the δ-axis). That is, when the axis error Δθ is non-zero, the current vector locus 71 rotates about the origin on the γ-δ axes (the γ-δ plane) and is inclined like the current vector locus 72.
The γ-axis and δ-axis components of the superposition currents are referred to as the γ-axis superposition current ihγ and the δ-axis superposition current ihδ. Their arithmetic product (ihγ×ihδ) contains a direct-current component that depends on the inclination of the elliptic current vector locus 72. The arithmetic product (ihγ×ihδ) is positive in the first and third quadrant of the current vector locus, and is negative in the second and fourth quadrant thereof; it thus includes no direct-current component when the ellipse is not inclined (like the current vector locus 71) and comes to include a direct-current component as the ellipse is inclined (like the current vector locus 72). In FIG. 8, I, II, III, and IV represent the first, second, third, and fourth quadrant on the γ-δ axes (the γ-δ plane).
In FIG. 9, with the horizontal axis representing time, curves 60 and 61 respectively represent the arithmetic product (ihγ×ihδ) and its direct-current component as observed when the axis error Δθ equals zero. In FIG. 10, with the horizontal axis representing time, curves 62 and 63 respectively represent the arithmetic product (ihγ×ihδ) and its direct-current component as observed when the axis error Δθ is non-zero. As will be understood from FIGS. 9 and 10, the direct-current component of the arithmetic product (ihγ×ihδ) equals zero when Δθ=0°, and is non-zero when Δθ≠0°. Moreover, the direct-current component increases as the axis error Δθ increase (the former is substantially proportional to the latter). If control is performed such that the direct-current component converges to zero, then the axis error Δθ converges to zero.
This relationship is exploited by the second speed estimator 27, which thus estimates the second estimated motor speed ω2e such that the direct-current component of the arithmetic product (ihγ×iδ) converges to zero. This estimation method, as in the first embodiment, differs from the one adopted in the first speed estimator 20. FIG. 11 is a block diagram showing an example of the internal configuration of the second speed estimator 27. The second speed estimator 27 shown in FIG. 11 includes a BPF (band-pass filter) 31, a coordinate converter 32, a multiplier 33, a LPF (low-pass filter) 34, a proportional-plus-integral calculator 35, and an integrator 36.
Suppose now that the frequency (electrical angular velocity) of the superposition voltages on fixed coordinate axes is ωh. In the second speed estimator 27, the BPF 31 is fed with the U-phase and V-phase currents iu and iv detected by the current detector 11. The BPF 31 is also fed with the W-phase current iw (the current flowing through the W-phase part of the armature winding) detected by the current detector 11, or the W-phase current iw calculated from the U-phase and V-phase currents iu and iv.
The BPF 31 extracts the ωh frequency components from iu, iv, and iw, and outputs them. The BPF 31 is a band-pass filter that receives, as the input signals thereto, iu, iv, and iw and whose pass band includes the frequency ωh; typically, for example, the center frequency of its pass band is set at ωh. The BPF 31 eliminates the components of the frequency of the drive currents.
Based on the second estimated rotor position ω2e fed from the integrator 36, the coordinate converter 32 converts the ωh frequency components of iu, iv, and iw into two-phase currents on the γ-δ axes, namely the γ-axis superposition current ihγ and the δ-axis superposition current ihδ. The γ-axis and δ-axis superposition currents ihγ and ihδ respectively represent the γ-axis and δ-axis components of the superposition currents superposed according to the superposition voltages. Incidentally, the second estimated rotor position θ2e is calculated based on the second estimated motor speed ω2e as will be described later and is, in strict terms, different from θe; in a synchronized state, however, it is regarded as equal to θe.
The multiplier 33 calculates the arithmetic product (ihγ×ihδ) of the γ-axis and δ-axis superposition currents ihγ and ihδ. The LPF 34 eliminates high-frequency components from the arithmetic product (ihγ×ihδ) and thereby extracts the direct-current component ihD of the arithmetic product (ihγ×ihδ).
To achieve PLL (phase locked loop) control, the proportional-plus-integral calculator 35, cooperating with the integrator 36 among others, performs proportional-plus-integral control to calculate the second estimated motor speed ω2e such that the direct-current component ihD outputted from the LPF 34 converges to zero. The integrator 36 integrates the second estimated motor speed ω2e outputted from the proportional-plus-integral calculator 35 to calculate the second estimated rotor position θ2e.
In a synchronized state, the second estimated motor speed ω2e calculated as described above is equal (or substantially equal) to the estimated motor speed ωe and the specified motor speed value ω*.
The LPF 29 shown in FIG. 6 eliminates the high-frequency components attributable to the superposition voltages from vγ* and vδ* from current controller 15 and from iγ and iδ from the coordinate converter 12. Based on vγ*, vδ*, iγ, and iδ thus having the high-frequency components attributable to the superposition voltages eliminated by the LPF 29, the first speed estimator 20 calculates the estimated motor speed ωe.
The second estimated motor speed ω2e calculated by the proportional-plus-integral calculator 35 is fed to the synchronization failure checker 23c shown in FIG. 6. The synchronization failure checker 23c is the same as the synchronization failure checker 23 shown in FIG. 3, and detects synchronization failure by the same method as the synchronization failure checker 23 shown in FIG. 3.
In the second speed estimator 27 shown FIG. 11, the coordinate converter 32 performs coordinate conversion after the BPF 31 has extracted the superposition current components; alternatively, the superposition current components may be extracted after coordinate conversion has been performed. In that case, a BPF for extracting the (ωh−ω2e) frequency components is provided in the stage succeeding the coordinate converter. For example, the center frequency of the pass band of this BPF is set at (ωh−ω2e).
In FIG. 6, the superposition voltages generated by the superposition voltage generator 28 is added to the specified voltage value after the coordinate conversion by the coordinate converter 18; alternatively, the addition of the superposition voltages may be performed before the coordinate conversion. In that case, the added superposition voltages are given, for example, a voltage vector that rotates at a frequency of (ωh−ω2e) on the γ-67 axes.
The voltage vector locus of the superposition voltages generated by the superposition voltage generator 28 does not necessarily have to be a perfect circle on the γ-δ axes. So long as the voltage vector locus of the superposition voltages on the γ-δ axes (the γ-δ plane) has a shape that includes the origin and that is symmetric about the γ-axis or δ-axis, the superposition voltages may be any voltages.
Here, “includes the origin” means that the origin on the γ-δ axes is located inside the “symmetric shape” mentioned above; “symmetric about the γ-axis” means that, of the voltage vector locus on the γ-δ axes, the part located in the first and second quadrants and the part located in the third and fourth quadrants are so shaped as to be line-symmetric about the γ-axis; “symmetric about the δ-axis” means that, of the voltage vector locus on the γ-δ axes, the part located in the first and fourth quadrants and the part located in the second and third quadrants are so shaped as to be line-symmetric about the δ-axis.
For example, the voltage vector locus of the superposition voltages on the γ-δ axes may be elliptic, having a minor or major direction aligned with the γ-axis; it may be linear, aligned with the γ-axis or δ-axis (i.e., a single-phase alternating voltage may instead be used); or it may be rectangular centered at the origin.
Instead of injecting superposition currents by superposing high-frequency superposition voltages, it is also possible to provide a superposition current generator (unillustrated) so that high-frequency superposition currents are superposed directly on iγ* and iδ*. The frequency of these superposition currents is significantly higher than that of the drive currents represented by iγ* and iδ*. In that case, a γ-axis superposition current ihγ and a δ-axis superposition current ihδ are respectively superposed on the specified γ-axis current value iγ* outputted from the magnetic flux controller 16 and on the specified δ-axis current value iδ* outputted from the speed controller 17.
Moreover, in this case, from vγ* and vδ* calculated by the current controller 15, those components vhγ and vhδ in the superposition voltages which are attributable to the superposition currents are extracted with a BPF. In addition, the direct-current component vhD of the arithmetic product (vhγ×vhδ) is extracted, and proportional-plus-integral control is performed to calculate the second estimated motor speed ω2e such that the direct-current component vhD converges to zero.
The current vector locus of the superposition currents generated by the superposition current has, on the γ-δ axes, a shape that includes the origin and that is symmetric about the γ-axis or δ-axis. For example, the current vector locus on the γ-δ axes may be perfectly circular, centered at the origin; it may be elliptic, having a minor or major axis aligned with the γ-axis; it may be linear, aligned with the γ-axis or δ-axis; or it may be rectangular, centered at the origin.
In case of synchronization failure, the first speed estimator 20 used to control the driving of the motor 1 cannot calculate the correct speed. By contrast, exploiting the magnetic salient poles of the rotor, the method involving the injection of high-frequency superposition currents can detect the rotor position accurately even when the rotation speed of the rotor is low, and can detect a rotation standstill surely even in case of one resulting from synchronization failure.
Thus, by comparing the second estimated motor speed ω2e with the specified motor speed value ω*, or with the estimated motor speed ωe (described in detail later) estimated from the axis error or the like based on the induction voltages (in other words, electromotive forces), it is possible to surely detect synchronization failure.
There have conventionally been proposed various methods for estimating the rotation speed and rotor position of a rotor, and different methods have different advantages and disadvantages; for example, some methods are suitable for low speeds, while other methods are suitable for high speeds, etc. Out of this consideration, estimation methods may be switched between at comparatively low speeds and at comparatively high speeds. In that case, typically, as the rotation speed becomes high, estimation by a method suitable for low speeds is stopped, and instead estimation by a method suitable for high speeds is started. By contrast, in this embodiment (and the first embodiment), the first speed estimator for controlling the driving of the motor 1 and the second speed estimator that is not directly concerned with the driving of the motor 1 are operated concurrently, irrespective of the rotation speed, to detect synchronization failure, thereby achieving a unique detection method.
First speed estimator: Next, a detailed description will be given of the first speed estimator 20 used in the embodiments described above. The first speed estimator 20 calculates the estimated motor speed ωe based on all or part of the specified γ-axis voltage value vγ* and specified δ-axis voltage value vδ* calculated by the current controller 15 and the γ-axis current iγ and δ-axis current iδ calculated by the coordinate converter 12, with the exception that, in the fourth embodiment (FIG. 6), the first speed estimator 20 calculates ωe based on all or part of vγ*, vδ*, iγ, and iδ having the high-frequency components eliminated therefrom as fed via the LPF 29.
Various methods for estimating the rotation speed of the rotor can be adopted in the first speed estimator 20.
One example is to use formula (17) below, which is disclosed in Japanese Patent No. 3411878. Specifically, with reference to the self-estimated ωe, the axis error Δθ is calculated according to formula (17) below; then, the axis error Δθ is made to converge to zero by proportional-plus-integral control, and thereby ωe is estimated. It should be noted here that, whereas according to Japanese Patent No. 3411878, the difference between the γ-axis and the d-axis is represented by Δθ, in this embodiment the difference d-axis and the γ-axis is represented by Δθ, with the result that opposite signs appear between in the formula for calculating Δθ in Japanese Patent No. 3411878 and in formula (17).




  
    
      
        Δθ
        =
        
          
            tan
            
              -
              1
            
          
          ⁡
          
            (
            
              
                -
                
                  (
                  
                    
                      v
                      γ
                      *
                    
                    -
                    
                      
                        R
                        a
                      
                      ⁢
                      
                        i
                        γ
                      
                    
                    +
                    
                      
                        ω
                        e
                      
                      ⁢
                      
                        L
                        q
                      
                      ⁢
                      
                        i
                        δ
                      
                    
                  
                  )
                
              
              
                
                  v
                  δ
                  *
                
                -
                
                  
                    R
                    a
                  
                  ⁢
                  
                    i
                    δ
                  
                
                -
                
                  
                    ω
                    e
                  
                  ⁢
                  
                    L
                    q
                  
                  ⁢
                  
                    i
                    γ
                  
                
              
            
            )
          
        
      
    
    
      
        (
        17
        )
      
    
  




Another example is, as disclosed in JP-A-2006-067656, to estimate the δ-axis magnetic flux, which is the component of the magnetic flux parallel to the δ-axis of the permanent magnet, based on the motor current and then make the δ-axis magnetic flux converge to zero to thereby estimate ωe.
As yet another example, the method (hereinafter referred to as the dmqm method) proposed by the applicant will be described below. The dmqm method is discussed in detail in, among others, Japanese Patent Application No. 2006-043398.
First, for easy understanding of the significance of the dmqm method, a brief description will be given of a common method for achieving maximum torque control.
A common motor control device controls a motor so that the axis error Δθ converges to zero. On the other hand, a formula for calculating the d-axis current id for the purpose of maximum torque control exploiting a reluctance torque is widely known, and generally, to achieve maximum torque control, the specified γ-axis current value iγ* is calculated according to formula (18) below.




  
    
      
        
          i
          γ
          *
        
        =
        
          
            
              Φ
              a
            
            
              2
              ⁢
              
                (
                
                  
                    L
                    q
                  
                  -
                  
                    L
                    d
                  
                
                )
              
            
          
          -
          
            
              
                
                  Φ
                  a
                  2
                
                
                  4
                  ⁢
                  
                    
                      (
                      
                        
                          L
                          q
                        
                        -
                        
                          L
                          d
                        
                      
                      )
                    
                    2
                  
                
              
              +
              
                i
                δ
                
                  *
                  2
                
              
            
          
        
      
    
    
      
        (
        18
        )
      
    
  




Achieving maximum torque control based on formula (18) above presupposes that the axis error Δθ is kept equal to zero, and, to calculate the axis error Δθ, it is necessary to previously calculate parameters needed in calculation, such as the q-axis inductance Lq. Thus, to achieve maximum torque control, first, the relevant parameters need to be so adjusted that the axis error Δθ is kept equal to zero; in addition, second, the parameters used in formula (18) need to be adjusted; furthermore, third, the specified γ-axis current value iγ* needs to be calculated constantly according to formula (18).
The dmqm method applicable to any of the embodiments contributes to simplifying the adjustment of calculation parameters and reducing the amount of calculation. Now, the dmqm method will be described.
FIGS. 12 and 13 are analysis model diagrams of the motor 1, illustrating the operation of the first speed estimator 20. FIG. 12 is a more detailed version of the analysis model diagram of FIG. 2.
The rotation axis whose direction coincides with that of the current vector to be fed to the motor 1 when maximum torque control is achieved is called the qm-axis. The axis that lags behind the qm-axis by an electrical angle of 90 degrees is called the dm-axis. The coordinate axes consisting of the dm and qm-axes is called the dm-qm axes.
As is well-known, the motor current that achieves maximum torque control has a positive q-axis component and a negative d-axis component. Thus, the qm-axis leads the q-axis in phase. The lead in phase occurs counter-clockwise in FIGS. 12 and 13.
The phase (angle) of the q-axis relative to the qm-axis is represented by θm, and the phase (angle) of the qm-axis relative to the δ-axis is represented by Δθm. Then, needless to say, the phase of the d-axis relative to the dm-axis also equals θm, and the phase of the dm-axis relative to the γ-axis also equals Δθm. What θm represents is the lead angle of the qm-axis (dm-axis) relative to the q-axis (d-axis). What Δθm represents is the axis error between the qm-axis and the δ-axis (the axis error between the dm-qm axes and the γ-δ axes). The value of Δθ, which is the axis error between the d- and γ-axes, is given by Δθ=Δθm+θm.
As described above, here, it is assumed that the dm-axis leads the d-axis in phase and that θm is then negative. Likewise, when the γ-axis leads the dm-axis in phase, Δθm is negative. The vectors (such as Em) shown in FIG. 13 will be described later.
The dm-axis and qm-axis components of the motor current Ia are referred to as the dm-axis current idm and the qm-axis current iqm respectively. The dm-axis and qm-axis components of the motor voltage Va are referred to as the dm-axis voltage vdm and the qm-axis voltage vqm respectively.
According to the dmqm method, the axis error Δθm between the qm-axis (dm-axis) and the δ-axis (γ-axis) is estimated, and thereby the γ-axis, which is an estimated axis, is made to converge to dm-axis (i.e., the axis error Δθm is made to converge to zero). Then, the motor current Ia is broken down into the qm-axis current iqm parallel to the qm-axis and the dm-axis current idm parallel to the dm-axis, and thereby the motor 1 is vector-controlled.
Even in this case, as with the common method for achieving maximum torque control previously described with reference to formula (18), the parameters for estimating the axis error Δθm (for making the axis error Δθm converge to zero) need to be adjusted. Here, however, completing those adjustments simultaneously brings about the completion of the adjustment of the parameters for achieving maximum torque control. That is, since the adjustment of the parameters for estimating the axis error is shared as the adjustment of the parameters for achieving maximum torque control, quite advantageously, the necessary adjustments can be made very easily.
Moreover, as will be clear from the definition of the qm-axis, the current locus of the motor current Ia as observed in maximum torque control lies on the qm-axis, as indicated by a solid line 82 in FIG. 14. Thus, in maximum torque control, it is not necessary to calculate the specified γ-axis current value iγ* according to a complicated formula like formula (18) above; this helps alleviate the burden of calculation. Here, the specified γ-axis current value iγ* equals zero or a predetermined value close to zero, irrespective of the value of iδ.
Now, the dmqm method will be described in more detail with reference to voltage equations. The extension induction voltage (extended electromotive force) equation on the real axes are given by formula (26), and the extension induction voltage (extended electromotive force) Eex is given by formula (27). In the formulae below, p represents the differentiation operator. It should be noted that, for consistency with the formula numbers used in the specification of Japanese Patent Application No. 2006-043398, in the present specification, formula numbers (19) to (25) are vacant.




  
    
      
        
          [
          
            
              
                
                  v
                  d
                
              
            
            
              
                
                  v
                  q
                
              
            
          
          ]
        
        =
        
          
            
              [
              
                
                  
                    
                      
                        R
                        a
                      
                      +
                      
                        pL
                        d
                      
                    
                  
                  
                    
                      
                        -
                        ω
                      
                      ⁢
                      
                        
                      
                      ⁢
                      
                        L
                        d
                      
                    
                  
                
                
                  
                    
                      ω
                      ⁢
                      
                        
                      
                      ⁢
                      
                        L
                        q
                      
                    
                  
                  
                    
                      
                        R
                        a
                      
                      +
                      
                        pL
                        d
                      
                    
                  
                
              
              ]
            
            ⁡
            
              [
              
                
                  
                    
                      i
                      d
                    
                  
                
                
                  
                    
                      i
                      q
                    
                  
                
              
              ]
            
          
          +
          
            [
            
              
                
                  0
                
              
              
                
                  
                    E
                    ex
                  
                
              
            
            ]
          
        
      
    
    
      
        (
        26
        )
      
    
  




Eex=ω((Ld−Lq)id+Φa)−(Ld−Lq)(piq)  (27)

Performing coordinate conversion so that formula (26) on the real axes is converted into one on the γ-δ axes estimated for control purposes gives formula (28). Ignoring the third term in the right side of formula (28) for the sake of simplicity gives formula (29).




  
    
      
        
          [
          
            
              
                
                  v
                  γ
                
              
            
            
              
                
                  v
                  δ
                
              
            
          
          ]
        
        =
        
          
            
              [
              
                
                  
                    
                      
                        R
                        a
                      
                      +
                      
                        pL
                        d
                      
                    
                  
                  
                    
                      
                        -
                        ω
                      
                      ⁢
                      
                        
                      
                      ⁢
                      
                        L
                        q
                      
                    
                  
                
                
                  
                    
                      ω
                      ⁢
                      
                        
                      
                      ⁢
                      
                        L
                        q
                      
                    
                  
                  
                    
                      
                        R
                        a
                      
                      +
                      
                        pL
                        d
                      
                    
                  
                
              
              ]
            
            ⁡
            
              [
              
                
                  
                    
                      i
                      γ
                    
                  
                
                
                  
                    
                      i
                      δ
                    
                  
                
              
              ]
            
          
          +
          
            
              E
              ex
            
            ⁡
            
              [
              
                
                  
                    
                      
                        -
                        sin
                      
                      ⁢
                      
                        
                      
                      ⁢
                      Δ
                      ⁢
                      
                        
                      
                      ⁢
                      θ
                    
                  
                
                
                  
                    
                      cos
                      ⁢
                      
                        
                      
                      ⁢
                      Δθ
                    
                  
                
              
              ]
            
          
          -
          
            
              (
              
                p
                ⁢
                
                  
                
                ⁢
                Δ
                ⁢
                
                  
                
                ⁢
                θ
              
              )
            
            ⁢
            
              
                L
                d
              
              ⁡
              
                [
                
                  
                    
                      
                        -
                        
                          i
                          δ
                        
                      
                    
                  
                  
                    
                      
                        i
                        γ
                      
                    
                  
                
                ]
              
            
          
        
      
    
    
      
        (
        28
        )
      
    
  
  
    
      
        
          [
          
            
              
                
                  v
                  γ
                
              
            
            
              
                
                  v
                  δ
                
              
            
          
          ]
        
        =
        
          
            
              [
              
                
                  
                    
                      
                        R
                        a
                      
                      +
                      
                        pL
                        d
                      
                    
                  
                  
                    
                      
                        -
                        ω
                      
                      ⁢
                      
                        
                      
                      ⁢
                      
                        L
                        q
                      
                    
                  
                
                
                  
                    
                      ω
                      ⁢
                      
                        
                      
                      ⁢
                      
                        L
                        q
                      
                    
                  
                  
                    
                      
                        R
                        a
                      
                      +
                      
                        pL
                        d
                      
                    
                  
                
              
              ]
            
            ⁡
            
              [
              
                
                  
                    
                      i
                      γ
                    
                  
                
                
                  
                    
                      i
                      δ
                    
                  
                
              
              ]
            
          
          +
          
            
              E
              ex
            
            ⁡
            
              [
              
                
                  
                    
                      
                        -
                        sin
                      
                      ⁢
                      
                        
                      
                      ⁢
                      Δθ
                    
                  
                
                
                  
                    
                      cos
                      ⁢
                      
                        
                      
                      ⁢
                      Δθ
                    
                  
                
              
              ]
            
          
        
      
    
    
      
        (
        29
        )
      
    
  




Rewriting formula (29) with respect to the dm-qm axes gives formula (30).




  
    
      
        
          [
          
            
              
                
                  v
                  dm
                
              
            
            
              
                
                  v
                  qm
                
              
            
          
          ]
        
        =
        
          
            
              [
              
                
                  
                    
                      
                        R
                        a
                      
                      +
                      
                        pL
                        d
                      
                    
                  
                  
                    
                      
                        -
                        ω
                      
                      ⁢
                      
                        
                      
                      ⁢
                      
                        L
                        q
                      
                    
                  
                
                
                  
                    
                      ω
                      ⁢
                      
                        
                      
                      ⁢
                      
                        L
                        q
                      
                    
                  
                  
                    
                      
                        R
                        a
                      
                      +
                      
                        pL
                        d
                      
                    
                  
                
              
              ]
            
            ⁡
            
              [
              
                
                  
                    
                      i
                      dm
                    
                  
                
                
                  
                    
                      i
                      qm
                    
                  
                
              
              ]
            
          
          +
          
            
              E
              ex
            
            ⁡
            
              [
              
                
                  
                    
                      
                        -
                        sin
                      
                      ⁢
                      
                        
                      
                      ⁢
                      
                        θ
                        m
                      
                    
                  
                
                
                  
                    
                      cos
                      ⁢
                      
                        
                      
                      ⁢
                      
                        θ
                        m
                      
                    
                  
                
              
              ]
            
          
        
      
    
    
      
        (
        30
        )
      
    
  




Here, it is assumed that formula (31) holds. Then, considering that id=iqm· sin θm, formula (32) holds.

Lq1iqm=sin θm{Φa+(Ld−Lq)id}  (31)

Lq1iqm=sin θm{Φa+(Ld−Lq)id}=sin θm{Φa+(Ld−Lq)iqm sin θm}  (32)

Rearranging formula (30) based on formula (32) gives formula (33). Here, Em is given by formula (34). The symbol Lq1 represents the virtual inductance that depends on θm. Here, Lq1 is defined for the sake of convenience to handle the Eex·sin θm present in the second term in the right side of formula (30) as the voltage drop across the virtual inductance. Incidentally, Lq1 is negative.




  
    
      
        
          [
          
            
              
                
                  v
                  dm
                
              
            
            
              
                
                  v
                  qm
                
              
            
          
          ]
        
        =
        
          
            
              [
              
                
                  
                    
                      
                        R
                        a
                      
                      +
                      
                        pL
                        d
                      
                    
                  
                  
                    
                      
                        -
                        ω
                      
                      ⁢
                      
                        
                      
                      ⁢
                      
                        (
                        
                          
                            L
                            q
                          
                          +
                          
                            L
                            
                              q
                              ⁢
                              
                                
                              
                              ⁢
                              1
                            
                          
                        
                        )
                      
                    
                  
                
                
                  
                    
                      ω
                      ⁢
                      
                        
                      
                      ⁢
                      
                        L
                        q
                      
                    
                  
                  
                    
                      
                        R
                        a
                      
                      +
                      
                        pL
                        d
                      
                    
                  
                
              
              ]
            
            ⁡
            
              [
              
                
                  
                    
                      i
                      dm
                    
                  
                
                
                  
                    
                      i
                      qm
                    
                  
                
              
              ]
            
          
          +
          
            
              E
              m
            
            ⁡
            
              [
              
                
                  
                    0
                  
                
                
                  
                    1
                  
                
              
              ]
            
          
        
      
    
    
      
        (
        33
        )
      
    
  




Em=(ω((Ld−Lq)id+Φa)−(Ld−Lq)(piq))cos θm=Eex cos θm  (34)

Now, an approximation is made such that the equation Lm=Lq+Lq1 holds (since θm depends on iq and iqm, Lq1 depends on iq and iqm; moreover, under the influence of magnetic saturation, Lq also depends on iq and iqm; the dependence of Lq1 on iq and the dependence of Lq on iq are integrated into Lm so that consideration is given to the influence of iq and iqm at the time of estimation). Now, formula (33) can be rearranged into formula (35).




  
    
      
        
          [
          
            
              
                
                  v
                  dm
                
              
            
            
              
                
                  v
                  qm
                
              
            
          
          ]
        
        =
        
          
            
              [
              
                
                  
                    
                      
                        R
                        a
                      
                      +
                      
                        pL
                        d
                      
                    
                  
                  
                    
                      
                        -
                        ω
                      
                      ⁢
                      
                        
                      
                      ⁢
                      
                        L
                        m
                      
                    
                  
                
                
                  
                    
                      ω
                      ⁢
                      
                        
                      
                      ⁢
                      
                        L
                        q
                      
                    
                  
                  
                    
                      
                        R
                        a
                      
                      +
                      
                        pL
                        d
                      
                    
                  
                
              
              ]
            
            ⁡
            
              [
              
                
                  
                    
                      i
                      dm
                    
                  
                
                
                  
                    
                      i
                      qm
                    
                  
                
              
              ]
            
          
          +
          
            
              E
              m
            
            ⁡
            
              [
              
                
                  
                    0
                  
                
                
                  
                    1
                  
                
              
              ]
            
          
        
      
    
    
      
        (
        35
        )
      
    
  




Further rearranging formula (35) gives formula (36) below. Here, Eexm is given by formula (37) below.




  
    
      
        
          
            
              
                [
                
                  
                    
                      
                        v
                        dm
                      
                    
                  
                  
                    
                      
                        v
                        qm
                      
                    
                  
                
                ]
              
              =
              
              ⁢
              
                
                  
                    [
                    
                      
                        
                          
                            
                              R
                              a
                            
                            +
                            
                              pL
                              d
                            
                          
                        
                        
                          
                            
                              -
                              ω
                            
                            ⁢
                            
                              
                            
                            ⁢
                            
                              L
                              m
                            
                          
                        
                      
                      
                        
                          
                            ω
                            ⁢
                            
                              
                            
                            ⁢
                            
                              L
                              m
                            
                          
                        
                        
                          
                            
                              R
                              a
                            
                            +
                            
                              pL
                              d
                            
                          
                        
                      
                    
                    ]
                  
                  ⁡
                  
                    [
                    
                      
                        
                          
                            i
                            dm
                          
                        
                      
                      
                        
                          
                            i
                            qm
                          
                        
                      
                    
                    ]
                  
                
                +
              
            
          
        
        
          
            
              
              ⁢
              
                
                  [
                  
                    
                      
                        0
                      
                    
                    
                      
                        
                          E
                          m
                        
                      
                    
                  
                  ]
                
                +
                
                  
                    
                      ω
                      ⁡
                      
                        (
                        
                          
                            L
                            q
                          
                          -
                          
                            L
                            m
                          
                        
                        )
                      
                    
                    ⁡
                    
                      [
                      
                        
                          
                            0
                          
                          
                            0
                          
                        
                        
                          
                            1
                          
                          
                            0
                          
                        
                      
                      ]
                    
                  
                  ⁡
                  
                    [
                    
                      
                        
                          
                            i
                            dm
                          
                        
                      
                      
                        
                          
                            i
                            qm
                          
                        
                      
                    
                    ]
                  
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  
                    [
                    
                      
                        
                          
                            
                              R
                              a
                            
                            +
                            
                              pL
                              d
                            
                          
                        
                        
                          
                            
                              -
                              ω
                            
                            ⁢
                            
                              
                            
                            ⁢
                            
                              L
                              m
                            
                          
                        
                      
                      
                        
                          
                            ω
                            ⁢
                            
                              
                            
                            ⁢
                            
                              L
                              m
                            
                          
                        
                        
                          
                            
                              R
                              a
                            
                            +
                            
                              pL
                              d
                            
                          
                        
                      
                    
                    ]
                  
                  ⁡
                  
                    [
                    
                      
                        
                          
                            i
                            dm
                          
                        
                      
                      
                        
                          
                            i
                            qm
                          
                        
                      
                    
                    ]
                  
                
                +
                
                  [
                  
                    
                      
                        0
                      
                    
                    
                      
                        
                          
                            E
                            m
                          
                          +
                          
                            
                              ω
                              ⁡
                              
                                (
                                
                                  
                                    L
                                    q
                                  
                                  -
                                  
                                    L
                                    m
                                  
                                
                                )
                              
                            
                            ⁢
                            
                              i
                              dm
                            
                          
                        
                      
                    
                  
                  ]
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  
                    [
                    
                      
                        
                          
                            
                              R
                              a
                            
                            +
                            
                              pL
                              d
                            
                          
                        
                        
                          
                            
                              -
                              ω
                            
                            ⁢
                            
                              
                            
                            ⁢
                            
                              L
                              m
                            
                          
                        
                      
                      
                        
                          
                            ω
                            ⁢
                            
                              
                            
                            ⁢
                            
                              L
                              m
                            
                          
                        
                        
                          
                            
                              R
                              a
                            
                            +
                            
                              pL
                              d
                            
                          
                        
                      
                    
                    ]
                  
                  ⁡
                  
                    [
                    
                      
                        
                          
                            i
                            dm
                          
                        
                      
                      
                        
                          
                            i
                            qm
                          
                        
                      
                    
                    ]
                  
                
                +
                
                  [
                  
                    
                      
                        0
                      
                    
                    
                      
                        
                          E
                          exm
                        
                      
                    
                  
                  ]
                
              
            
          
        
      
    
    
      
        (
        36
        )
      
    
  
  
    
      
        
          
            
              
                E
                exm
              
              =
              
              ⁢
              
                (
                
                  
                    ω
                    ⁡
                    
                      (
                      
                        
                          
                            (
                            
                              
                                L
                                d
                              
                              -
                              
                                L
                                q
                              
                            
                            )
                          
                          ⁢
                          
                            i
                            d
                          
                        
                        +
                        
                          Φ
                          a
                        
                      
                      )
                    
                  
                  -
                  
                    
                      (
                      
                        
                          L
                          d
                        
                        -
                        
                          L
                          q
                        
                      
                      )
                    
                    ⁢
                    
                      (
                      
                        pi
                        q
                      
                      )
                    
                  
                
                )
              
            
          
        
        
          
            
              
              ⁢
              
                
                  cos
                  ⁢
                  
                    
                  
                  ⁢
                  
                    θ
                    m
                  
                
                +
                
                  
                    ω
                    ⁡
                    
                      (
                      
                        
                          L
                          q
                        
                        -
                        
                          L
                          m
                        
                      
                      )
                    
                  
                  ⁢
                  
                    i
                    dm
                  
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  E
                  m
                
                +
                
                  
                    ω
                    ⁡
                    
                      (
                      
                        
                          L
                          q
                        
                        -
                        
                          L
                          m
                        
                      
                      )
                    
                  
                  ⁢
                  
                    i
                    dm
                  
                
              
            
          
        
      
    
    
      
        (
        37
        )
      
    
  




Suppose that there is an axis error Δθm between the γ-δ axes and the dm-qm axes; then, formula (36) can be rearranged into formula (38) below. That is, just as formula (26) is rearranged into formula (28), performing coordinate conversion so that formula (36) on the dm-qm axes is converted into one on the γ-δ axes gives formula (38).




  
    
      
        
          [
          
            
              
                
                  v
                  γ
                
              
            
            
              
                
                  v
                  δ
                
              
            
          
          ]
        
        =
        
          
            
              [
              
                
                  
                    
                      
                        R
                        a
                      
                      +
                      
                        pL
                        d
                      
                    
                  
                  
                    
                      
                        -
                        ω
                      
                      ⁢
                      
                        
                      
                      ⁢
                      
                        L
                        m
                      
                    
                  
                
                
                  
                    
                      ω
                      ⁢
                      
                        
                      
                      ⁢
                      
                        L
                        m
                      
                    
                  
                  
                    
                      
                        R
                        a
                      
                      +
                      
                        pL
                        d
                      
                    
                  
                
              
              ]
            
            ⁡
            
              [
              
                
                  
                    
                      i
                      γ
                    
                  
                
                
                  
                    
                      i
                      δ
                    
                  
                
              
              ]
            
          
          +
          
            
              E
              exm
            
            ⁡
            
              [
              
                
                  
                    
                      
                        -
                        sin
                      
                      ⁢
                      
                        
                      
                      ⁢
                      Δ
                      ⁢
                      
                        
                      
                      ⁢
                      
                        θ
                        m
                      
                    
                  
                
                
                  
                    
                      cos
                      ⁢
                      
                        
                      
                      ⁢
                      
                        Δθ
                        m
                      
                    
                  
                
              
              ]
            
          
          -
          
            
              (
              
                p
                ⁢
                
                  
                
                ⁢
                
                  Δθ
                  m
                
              
              )
            
            ⁢
            
              
                L
                d
              
              ⁡
              
                [
                
                  
                    
                      
                        -
                        
                          i
                          δ
                        
                      
                    
                  
                  
                    
                      
                        i
                        γ
                      
                    
                  
                
                ]
              
            
          
        
      
    
    
      
        (
        38
        )
      
    
  




When approximations are made such that pΔθm≈0, idm≈0, and (Ld−Lq) (piq)≈0, then Eexm given by formula (37) is approximated by formula (39) below.




  
    
      
        
          
            
              
                E
                exm
              
              =
              
              ⁢
              
                
                  
                    (
                    
                      
                        ω
                        ⁡
                        
                          (
                          
                            
                              
                                (
                                
                                  
                                    L
                                    d
                                  
                                  -
                                  
                                    L
                                    q
                                  
                                
                                )
                              
                              ⁢
                              
                                i
                                d
                              
                            
                            +
                            
                              Φ
                              a
                            
                          
                          )
                        
                      
                      -
                      
                        
                          (
                          
                            
                              L
                              d
                            
                            -
                            
                              L
                              q
                            
                          
                          )
                        
                        ⁢
                        
                          (
                          
                            pi
                            q
                          
                          )
                        
                      
                    
                    )
                  
                  ⁢
                  cos
                  ⁢
                  
                    
                  
                  ⁢
                  
                    θ
                    m
                  
                
                +
              
            
          
        
        
          
            
              
              ⁢
              
                
                  ω
                  ⁡
                  
                    (
                    
                      
                        L
                        q
                      
                      -
                      
                        L
                        m
                      
                    
                    )
                  
                
                ⁢
                
                  i
                  dm
                
              
            
          
        
        
          
            
              ≈
              
              ⁢
              
                (
                
                  
                    ω
                    ⁡
                    
                      (
                      
                        
                          
                            (
                            
                              
                                L
                                d
                              
                              -
                              
                                L
                                q
                              
                            
                            )
                          
                          ⁢
                          
                            i
                            δ
                          
                          ⁢
                          sin
                          ⁢
                          
                            
                          
                          ⁢
                          
                            θ
                            m
                          
                        
                        +
                        
                          Φ
                          a
                        
                      
                      )
                    
                  
                  -
                  
                    
                      (
                      
                        
                          L
                          d
                        
                        -
                        
                          L
                          q
                        
                      
                      )
                    
                    ⁢
                    
                      (
                      
                        pi
                        q
                      
                      )
                    
                  
                
                )
              
            
          
        
        
          
            
              
              ⁢
              
                
                  cos
                  ⁢
                  
                    
                  
                  ⁢
                  
                    θ
                    m
                  
                
                +
                
                  
                    ω
                    ⁡
                    
                      (
                      
                        
                          L
                          q
                        
                        -
                        
                          L
                          m
                        
                      
                      )
                    
                  
                  ⁢
                  
                    i
                    dm
                  
                
              
            
          
        
        
          
            
              ≈
              
              ⁢
              
                
                  ω
                  ⁡
                  
                    (
                    
                      
                        
                          (
                          
                            
                              L
                              d
                            
                            -
                            
                              L
                              q
                            
                          
                          )
                        
                        ⁢
                        
                          i
                          δ
                        
                        ⁢
                        sin
                        ⁢
                        
                          
                        
                        ⁢
                        
                          θ
                          m
                        
                      
                      +
                      
                        Φ
                        a
                      
                    
                    )
                  
                
                ⁢
                cos
                ⁢
                
                  
                
                ⁢
                
                  θ
                  m
                
              
            
          
        
      
    
    
      
        (
        39
        )
      
    
  




Substituting Lm=Lq+Lq1 in formula (32) above, then solving the resulting formula for θm, and then assuming that iδ≈iqm gives formula (40) below. As expressed by formula (40), θm is a function of iδ, and accordingly Eexm is a function of iδ.




  
    
      
        
          θ
          m
        
        =
        
          
            sin
            
              -
              1
            
          
          (
          
            
              
                Φ
                a
              
              -
              
                
                  
                    Φ
                    a
                    2
                  
                  +
                  
                    4
                    ⁢
                    
                      (
                      
                        
                          L
                          q
                        
                        -
                        
                          L
                          m
                        
                      
                      )
                    
                    ⁢
                    
                      (
                      
                        
                          L
                          q
                        
                        -
                        
                          L
                          d
                        
                      
                      )
                    
                    ⁢
                    
                      i
                      δ
                      2
                    
                  
                
              
            
            
              2
              ⁢
              
                
                  i
                  δ
                
                ⁡
                
                  (
                  
                    
                      L
                      q
                    
                    -
                    
                      L
                      d
                    
                  
                  )
                
              
            
          
          )
        
      
    
    
      
        (
        40
        )
      
    
  




Now, with reference to FIG. 13, a description will be given of the relationship inductance Lq1.
Define Φexm=Eexm/ω; then Φexm is the sum of Φm and (Lq−Lm)idm. Thus, in the rotating coordinate system, Φexm, like Φm, is a flux linkage vector on the dm-axis. In maximum torque control, as described above, idm≈0, and hence Φexm is (substantially) equal to Φm.
FIG. 15 shows an example of the internal configuration of the first speed estimator 20, as observed when the dmqm method is adopted. The first speed estimator 20 shown in FIG. 15 includes an axis error estimator 41 and a proportional-plus-integral calculator 42.
The axis error estimator 41 calculates the axis error Δθm based on all or part of vγ*, vδ*, iγ, and iδ. To achieve PLL (phase locked loop), the proportional-plus-integral calculator 42 performs proportional-plus-integral control to calculate the estimated motor speed ωe such that the axis error Δθm calculated by the axis error estimator 41 converges to zero.
Various methods for estimating the axis error Δθm can be adopted to make the axis error estimator 41 estimate the axis error Δθm. Presented below are four examples of the calculation method adopted in the axis error estimator 41.
It should be noted that, when using any formula stated in the present specification, the axis error estimator 41 substitutes the values of vγ*, vδ*, and ωe for the values of vγ, vδ, and ω wherever they appear in the formula. It should also be noted that any description (for example, regarding how to determine the value of Lm) given in connection with a particular calculation method is applicable to all other calculation methods.
First Calculation Method: First, a first method for calculating the axis error Δθm will be described. In the first calculation method, the induction voltage Eex generated by the motor 1 is broken into an induction voltage vector on the qm-axis and an induction voltage vector on the dm-axis. Then, based on the induction voltage vector Eexm(≈Em; see FIG. 13), which is the induction voltage vector on the qm-axis, the axis error Δθm is calculated, and thereby the phase (θe) of the γ-axis, which is an axis estimated for control purposes, is calculated (i.e., the rotor position is estimated).
The γ-axis and δ-axis components of the induction voltage vector Eexm are represented by Eexmγ and Eexmδ respectively. Then, as will be clear from FIG. 13, Δθm=tan−1(−Eexmγ/Eexmδ) holds. Using a rearranged form of the first and second lines of formula (38) (a determinant), Δθm is given by formula (41) below (assuming that the third term in the right side of formula (38) is ignored). In formula (41), eventually, on the assumption that Δθm is small, an approximation is made such that tan−1(−Eexmγ/Eexmδ)≈(−Eexmγ/Eexmδ).




  
    
      
        
          
            
              
                Δθ
                m
              
              =
              
              ⁢
              
                
                  tan
                  
                    -
                    1
                  
                
                ⁢
                
                  
                    -
                    
                      E
                      
                        exm
                        ⁢
                        
                          
                        
                        ⁢
                        γ
                      
                    
                  
                  
                    E
                    
                      exm
                      ⁢
                      
                        
                      
                      ⁢
                      δ
                    
                  
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  tan
                  
                    -
                    1
                  
                
                ⁢
                
                  
                    -
                    
                      (
                      
                        
                          v
                          γ
                        
                        -
                        
                          
                            (
                            
                              
                                R
                                a
                              
                              +
                              
                                pL
                                d
                              
                            
                            )
                          
                          ⁢
                          
                            i
                            γ
                          
                        
                        +
                        
                          ω
                          ⁢
                          
                            
                          
                          ⁢
                          
                            L
                            m
                          
                          ⁢
                          
                            i
                            δ
                          
                        
                      
                      )
                    
                  
                  
                    
                      v
                      δ
                    
                    -
                    
                      
                        (
                        
                          
                            R
                            a
                          
                          +
                          
                            pL
                            d
                          
                        
                        )
                      
                      ⁢
                      
                        i
                        δ
                      
                    
                    -
                    
                      ω
                      ⁢
                      
                        
                      
                      ⁢
                      
                        L
                        m
                      
                      ⁢
                      
                        i
                        γ
                      
                    
                  
                
              
            
          
        
        
          
            
              ≈
              
              ⁢
              
                -
                
                  
                    
                      v
                      γ
                    
                    -
                    
                      
                        (
                        
                          
                            R
                            a
                          
                          +
                          
                            pL
                            d
                          
                        
                        )
                      
                      ⁢
                      
                        i
                        γ
                      
                    
                    +
                    
                      ω
                      ⁢
                      
                        
                      
                      ⁢
                      
                        L
                        m
                      
                      ⁢
                      
                        i
                        δ
                      
                    
                  
                  
                    
                      v
                      δ
                    
                    -
                    
                      
                        (
                        
                          
                            R
                            a
                          
                          +
                          
                            pL
                            d
                          
                        
                        )
                      
                      ⁢
                      
                        i
                        δ
                      
                    
                    -
                    
                      ω
                      ⁢
                      
                        
                      
                      ⁢
                      
                        L
                        m
                      
                      ⁢
                      
                        i
                        γ
                      
                    
                  
                
              
            
          
        
      
    
    
      
        (
        41
        )
      
    
  




When calculating Δθm according to formula (41), the axis error estimator 41 can ignore the differentiation terms pLdiγ and pLdiδ. Moreover, to calculate the value of Lm needed to calculate Δθm, the axis error estimator 41 uses formula (42) below. Formula (42) is obtained by substituting idm=0 and formulae (43) and (44) below in formula (32) above, and then solving the resulting formula for Lq1.




  
    
      
        
          
            
              
                L
                m
              
              =
              
                
                  L
                  q
                
                +
                
                  L
                  
                    q
                    ⁢
                    
                      
                    
                    ⁢
                    1
                  
                
              
            
          
        
        
          
            
              =
              
                
                  L
                  q
                
                +
                
                  
                    
                      i
                      d
                    
                    ⁢
                    
                      {
                      
                        
                          Φ
                          a
                        
                        +
                        
                          
                            (
                            
                              
                                L
                                d
                              
                              -
                              
                                L
                                q
                              
                            
                            )
                          
                          ⁢
                          
                            i
                            d
                          
                        
                      
                      }
                    
                  
                  
                    
                      i
                      d
                      2
                    
                    +
                    
                      i
                      q
                      2
                    
                  
                
              
            
          
        
      
    
    
      
        (
        42
        )
      
    
  
  
    
      
        
          i
          qm
        
        =
        
          
            
              i
              d
              2
            
            +
            
              i
              q
              2
            
          
        
      
    
    
      
        (
        43
        )
      
    
  




among Eex, Em, and Eem. Consider Eex, Em and Eexm as voltage vectors in a rotating coordinate system. Then, Eex can be called an extension induction voltage vector. The extension induction voltage vector Eex is an induction voltage vector on the q-axis. The extension induction voltage vector Eex can be broken into an induction voltage vector on the qm-axis and an induction voltage vector on the dm-axis. As will be understood from formula (34) above, of these vectors thus broken down, the induction voltage vector on the qm-axis is Em; and the induction voltage vector (Eex·sin θm) on the dm-axis, indicated by the reference numeral 80 in FIG. 13 is the voltage drop vector attributable to the virtual inductance Lq1.

As will be understood from a comparison between formulae (34) and (37), Eexm is the sum of Em and ω(Lq−Lm)idm. Thus, in the rotating coordinate system, Eexm, like Em, is an induction voltage vector on the qm-axis. In maximum torque control, as described above, idm≈0, and hence Eexm is (substantially) equal to Em.
Next, with reference to FIG. 13, a description will be given of the magnetic flux corresponding to Eex, Em, and Eexm. What Eex represents is an induction voltage generated by the flux linkage Φex of the motor 1 and the rotation of the motor 1. Put the other way around, Φex is calculated by dividing Eex by ω (assuming that the transient term (the second term in the right side) of Eex given by formula (27) is ignored).
Consider Φex as a flux linkage vector in a rotating coordinate system; then, the flux linkage vector Φex is a flux linkage vector on the d-axis. The flux linkage vector Φex can be broken into a flux linkage vector on the qm-axis and a flux linkage vector on the dm-axis. Of these vectors thus broken down, the flux linkage vector on the dm-axis is represented by Φm, and fulfills Φm=Em/ω; the flux linkage vector (Φex·sin θm) on the qm-axis, indicated by the reference numeral 81 in FIG. 13, is a magnetic flux vector attributable to the virtual




  
    
      
        
          sin
          ⁢
          
            
          
          ⁢
          
            θ
            m
          
        
        =
        
          
            i
            d
          
          
            
              
                i
                d
                2
              
              +
              
                i
                q
                2
              
            
          
        
      
    
    
      
        (
        44
        )
      
    
  




Rearranging formula (42) above using formula (45), which gives the d-axis current id that achieves maximum torque control, and formula (43), which gives (an approximation of) the relationship among id, iq, and iqm, gives Lm as a function of iqm (i.e., the terms of id and iq are eliminated from the calculation formula of Lm). Hence, assuming that iδ≈iqm, the axis error estimator 41 can calculate, based on iδ, the value of Lm given as a function of iqm. Then, based on the Lm thus calculated, the axis error Δθm is calculated according to formula (41).




  
    
      
        
          i
          d
        
        =
        
          
            
              Φ
              a
            
            
              2
              ⁢
              
                (
                
                  
                    L
                    q
                  
                  -
                  
                    L
                    d
                  
                
                )
              
            
          
          -
          
            
              
                
                  Φ
                  a
                  2
                
                
                  4
                  ⁢
                  
                    
                      (
                      
                        
                          L
                          q
                        
                        -
                        
                          L
                          d
                        
                      
                      )
                    
                    2
                  
                
              
              +
              
                i
                q
                2
              
            
          
        
      
    
    
      
        (
        45
        )
      
    
  




Alternatively, the value of Lm may be calculated by assuming that iδ≈iqm and using an approximation formula that gives Lm as a function of iδ; or different values of Lm corresponding to different values of iδ may be previously prepared in the form of table data so that the value of Lm is determined by referring to the table data.
FIG. 16 is a graph obtained in a numerical example, showing the dependence of Ld, Lq, and Lm on iqm (assuming that iγ*≈0). As shown in FIG. 16, the value of Lm depends on iqm, and increases as iqm increases. The graph shows that the value of Lm that achieves maximum torque control is located far closer to Ld than to Lq.
The value of Lm is so determined as to fulfill formula (46) or (47) below. This produces an intentional deviation between the d- and γ-axes, and makes iγ*≈0, thereby achieving control approximate to maximum torque control.

Ld≦Lm<Lq  (46)

Ld≦Lm<(Ld+Lq)/2  (47)

The value of Lm may be fixed; that is, Lm may be given a value that remains fixed irrespective of the value of iδ. In FIG. 17, a solid line 83 represents the relationship between the d-axis and q-axis currents id and iq as observed when Lm is given a predetermined fixed value, whereas a broken line 84 represents the relationship between the d-axis and q-axis currents id and iq as observed when maximum torque control is performed in an ideal manner. FIG. 17 shows that the curves described by the solid line 83 and the broken line 84 are very similar.
Second Calculation Method: Next, a second method for calculating the axis error Δθm will be described. Also in the second calculation method, as in the first calculation method described above, the axis error Δθm is calculated based on the induction voltage vector Eexm, and thereby the phase (θe) of the γ-axis, which is an axis estimated for control purposes, is calculated (i.e. the rotor position is estimated). In the second calculation method, however, the δ-axis component Eexmδ of the induction voltage vector Eexm is not used. Specifically, the axis error Δθm is calculated according to formula (48). In formula (48), eventually, on the assumption that Δθm is small, an approximation is made such that sin−1(−Eexmγ/Eexm)≈(−Eexmγ/Eexm).




  
    
      
        
          
            
              
                Δθ
                m
              
              =
              
              ⁢
              
                
                  sin
                  
                    -
                    1
                  
                
                ⁡
                
                  (
                  
                    
                      -
                      
                        E
                        
                          exm
                          ⁢
                          
                            
                          
                          ⁢
                          γ
                        
                      
                    
                    
                      E
                      exm
                    
                  
                  )
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  sin
                  
                    -
                    1
                  
                
                ⁢
                
                  
                    -
                    
                      (
                      
                        
                          v
                          γ
                        
                        -
                        
                          
                            (
                            
                              
                                R
                                a
                              
                              +
                              
                                pL
                                d
                              
                            
                            )
                          
                          ⁢
                          
                            i
                            γ
                          
                        
                        +
                        
                          ω
                          ⁢
                          
                            
                          
                          ⁢
                          
                            L
                            m
                          
                          ⁢
                          
                            i
                            δ
                          
                        
                      
                      )
                    
                  
                  
                    E
                    exm
                  
                
              
            
          
        
        
          
            
              ≈
              
              ⁢
              
                
                  
                    v
                    γ
                  
                  -
                  
                    
                      (
                      
                        
                          R
                          a
                        
                        +
                        
                          pL
                          d
                        
                      
                      )
                    
                    ⁢
                    
                      i
                      γ
                    
                  
                  +
                  
                    ω
                    ⁢
                    
                      
                    
                    ⁢
                    
                      L
                      m
                    
                    ⁢
                    
                      i
                      δ
                    
                  
                
                
                  E
                  exm
                
              
            
          
        
      
    
    
      
        (
        48
        )
      
    
  




When calculating Δθm using formula (48), the axis error estimator 41 can ignore the differentiation term pLdiγ. The value of Lm is determined in the same manner as in the first calculation method described above.
The value of Eexm in formula (48) is calculated according to formula (39) above. As an approximation formula for calculating Eexm, for example, formula (49), (50), or (51) below can be used. Formula (49) approximates formula (37) by making approximations such that p Δθm≈0, idm≈0, and (Ld−Lq)(piq)≈0; formula (50) further approximates formula (49) by making an approximation such that cos θm≈1; formula (51) further approximates formula (50) by making an approximation such that (Ld−Lq) iδ sin θm<<Φa. When (49), (50), or (51) is used, ωe is substituted for the value of ω.

Eexm≈ω((Ld−Lq)iδ sin θm+Φa)cos θm  (49)

Eexm≈ω((Ld−Lq)iδ sin θm+Φa)  (50)

Eexm≈ωΦa  (51)

The value of θm in formula (49) etc. is calculated according to formula (40) above. As will be understood from formula (40), θm is a function of iδ, and thus Eexm also is a function of iδ. Since calculating Eexm requires complicated calculation, it is preferable to use an appropriate approximation formula to calculate it. Alternatively, different values of Eexm corresponding to different values of iδ may be previously prepared in the form of table data so that the value of Eexm is determined by referring to the table data.
Third Calculation Method: Next, a third method for calculating the axis error Δθm will be described. In the third calculation method, the flux linkage Φex produced around the armature winding of the motor 1 is broken down into a flux linkage vector on the qm-axis and a flux linkage vector on the dm-axis. Then, based on the flux linkage vector Φexm(≈Φm; see FIG. 13), which is the flux linkage vector on the dm-axis, the axis error Δθm is calculated, and thereby the phase (θe) of the γ-axis, which is an axis estimated for control purposes, is calculated (i.e., the rotor position is estimated).
The γ-axis and δ-axis components of the flux linkage vector Φexm is represented by Φexmγ and Φexmδ respectively. Then, as will be clear from FIG. 13, Δθm=tan−1(−Φexmδ/Φexmγ) holds. Here, Φexm equals Eexm divided by ω. Accordingly, in the third calculation method, the axis error Δθm is calculated according to formula (52) below, which is a rearranged version of formula (41) used in the first calculation method. In other respects, the third calculation method is the same as the first calculation method.




  
    
      
        
          Δθ
          m
        
        =
        
          
            
              tan
              
                -
                1
              
            
            ⁢
            
              
                Φ
                
                  exm
                  ⁢
                  
                    
                  
                  ⁢
                  δ
                
              
              
                Φ
                
                  exm
                  γ
                
              
            
          
          =
          
            
              tan
              
                -
                1
              
            
            (
            
              
                
                  -
                  
                    E
                    
                      exm
                      ⁢
                      
                        
                      
                      ⁢
                      γ
                    
                  
                
                ω
              
              
                
                  E
                  
                    exm
                    ⁢
                    
                      
                    
                    ⁢
                    δ
                  
                
                ω
              
            
            )
          
        
      
    
    
      
        (
        52
        )
      
    
  




Fourth Calculation Method: Next, a fourth method for calculating the axis error Δθm will be described. In the fourth calculation method, as in the third calculation method described above, the axis error Δθm is calculated based on the flux linkage vector Φexm, and thereby the phase (θe) of the γ-axis, which is an axis estimated for control purposes, is calculated (i.e., the rotor position is estimated). In the fourth calculation method, however, the γ-axis component Φexmγ of the flux linkage vector Φexm is not used. Specifically, the axis error Δθm is calculated according to formula (53) below, which is a rearranged version of formula (48) used in the second calculation method. In other respect, the fourth calculation method is the same as the second calculation method.




  
    
      
        
          Δθ
          m
        
        =
        
          
            
              sin
              
                -
                1
              
            
            ⁡
            
              (
              
                
                  -
                  
                    Φ
                    
                      exm
                      ⁢
                      
                        
                      
                      ⁢
                      δ
                    
                  
                
                
                  Φ
                  exm
                
              
              )
            
          
          =
          
            
              sin
              
                -
                1
              
            
            (
            
              
                
                  -
                  
                    E
                    
                      exm
                      ⁢
                      
                        
                      
                      ⁢
                      γ
                    
                  
                
                ω
              
              
                
                  E
                  exm
                
                ω
              
            
            )
          
        
      
    
    
      
        (
        53
        )
      
    
  




Any calculation method other than the first to fourth calculation methods described above may be used. For example, the axis error Δθm may be calculated based on the error current between the current on the dm-qm axes (the current in the motor model) and the current on the γ-δ axes.
Next, a description will be given of how sensorless vector control is performed in case of synchronization failure. FIG. 18 is a vector diagram obtained in a synchronized state (in synchronized operation) when, for example, the dmqm method is adopted.
Here, the rotation speed (rotation angular speed) of the rotor is represented by ωr, and the rotation speed of the γ-δ axis, that is, the rotation speed (rotation angular speed) of the current vector of the motor current Ia is represented by ωs. In a synchronized state, ωr and ωs, are equal.
The motor voltage Va is given as the sum of the voltage component ωrΦex produced by the rotation of the rotor, the voltage component RaIa produced across the motor resistance, and the voltage component ωsLqiδ produced by the inductance and the motor current.
When the dmqm method is adopted, the difference between Lq and the estimation parameter Lm constantly produces a non-zero axis error Δθ between the γ-axis and the γ-axis. That is, an axis error Δθ commensurate with the difference between Lq and Lm is produced such that a synchronized state (synchronized operation) is maintained. In this state, as will be easily understood from the vector diagram of FIG. 18, formula (54) below holds.

ωsLmiδ=ωsLqiδ−ωrΦex sin Δθ  (54)

In a state of synchronization failure, ωr equals 0. Thus, as will be easily understood from the vector diagram of FIG. 18, the value of the δ-axis voltage in a state of synchronization failure is evidently smaller than in a synchronized state. Thus, as in the second embodiment, by estimating the δ-axis voltage, it is possible to detect synchronization failure.
Instead, it is also possible to divide the δ-axis voltage by Φex (or simply by Φa) to calculate ωs and detect synchronization failure based on the thus calculated ωs. This is possible because, in case of synchronization failure, the calculated ωs has an extremely small value. Examples of methods that exploit ωs are presented in the first and fourth embodiments. It is also possible to detect synchronization failure by estimating the δ-axis current, and an example of this method is presented in the third embodiment.
MODIFICATIONS AND VARIATIONS
Any feature of the embodiments described above can be applied, unless inconsistent, to any embodiment other than that in connection with which the feature is specifically described.
The synchronization failure checker 23, 23a, 23b, and 23c in FIGS. 3 to 6 each serve as a synchronization failure detector. In FIG. 3, the second speed estimator 22 can also be regarded as part of the synchronization failure detector. In FIG. 4, the applied voltage estimator 25 can also be regarded as part of the synchronization failure detector. In FIG. 5, the supplied current estimator 26 can also be regarded as part of the synchronization failure detector. In FIG. 6, the second speed estimator 27 can also be regarded as part of the synchronization failure detector, and the superposition voltage generator 28 and adder 30 can also be regarded as part of the synchronization failure detector.
In each embodiment, the coordinate converter 12 and 18, the subtracter 13 and 14, and the current controller 15 constitute a specified voltage value calculator; the magnetic flux controller 16, the speed controller 17, and subtracter 19 constitute a specified current calculator.
What the values iu, iv, iγ, and iδ fed into and out of the coordinate converter 12 represent are the values of currents as measured by the current detector (current measurer) 11, and therefore these currents can be called measured currents. The current detector 11 may be so configured as to directly detect the motor current as shown in FIG. 3 etc., or may instead be so configured as to reproduce the motor current from instantaneous currents in the supplied direct current and thereby detect motor current. Also in this case, the motor current so detected can be called as a current as measured (i.e., a measured current).
In each embodiment, the motor control device is realized, for example, with software (a program) incorporated in a general-purpose microcomputer or the like. Needless to say, the motor control device may be realized with hardware alone, instead of software (a program).
The present invention is suitable for electric devices of any kind that use a motor; for example, it is suitable for electric cars that run by exploiting the rotation of a motor, and compressors and the like used in air conditioners and the like.
Any of the different specified values (iγ*, iδ*, vγ*, vδ*, etc.) and other state quantities (ωe, ω2e, etc.) mentioned above, that is, any value that needs to be derived (or estimated) as necessary may be derived (or estimated) in any manner. That is, such values may be derived (or estimated), for example, through calculation performed within the motor control device (3, 3a, etc.), or may be derived (or estimated) from a previously set data table.

## Claims

1. A motor control device comprising a first speed estimator estimating a rotation speed of a rotor of a permanent-magnet synchronous motor, the motor control device controlling the motor so that a first estimated rotation speed estimated by the first speed estimator follows a specified speed value,
the motor control device further comprising:
a synchronization failure detector estimating, based on one or two of a first, a second, and a third state quantity, another of the first, second, and third state quantities and detecting synchronization failure based on a result of the estimation,
the first state quantity being related to a motor voltage applied to the motor,
the second state quantity being related to the rotation speed of the rotor, and
the third state quantity being related to a motor current supplied to the motor.

2. The motor control device of claim 1, wherein
let an axis parallel to a magnetic flux produced by a permanent magnet provided on the rotor be a d-axis,
let an axis estimated, for control purposes, to correspond to the d-axis be a γ-axis, and
let an axis estimated, for control purposes, to be perpendicular to the γ-axis be a δ-axis, then
the second state quantity is represented by the first estimated rotation speed or the specified speed value,
the synchronization failure detector estimates as the first state quantity a δ-axis component of the motor voltage based on
the second state quantity alone or
the second state quantity and the third state quantity represented by a measured current of the motor current or a specified current value to be followed by the measured current, and

the synchronization failure detector detects synchronization failure based on a thus estimated δ-axis voltage as estimated as the first state quantity.

3. The motor control device of claim 2,
wherein the synchronization failure detector detects synchronization failure based on a ratio or deviation between the estimated δ-axis voltage and a specified δ-axis voltage value to be followed by the δ-axis component of the motor voltage.

4. The motor control device of claim 1, wherein
let an axis parallel to a magnetic flux produced by a permanent magnet provided on the rotor be a d-axis,
let an axis estimated, for control purposes, to correspond to the d-axis be a γ-axis, and
let an axis estimated, for control purposes, to be perpendicular to the γ-axis be a δ-axis, then
the first state quantity is represented by a specified δ-axis voltage value to be followed by a δ-axis component of the motor voltage,
the second state quantity is represented by the first estimated rotation speed or the specified speed value, and
the synchronization failure detector estimates as the third state quantity a δ-axis component of the motor current based on the first and second state quantities, and detects synchronization failure based on a thus estimated δ-axis current as estimated as the third state quantity.

5. The motor control device of claim 4,
wherein the synchronization failure detector detects the synchronization failure based on
a ratio or deviation between the estimated δ-axis current and a δ-axis component of the measured current or a specified δ-axis current value to be followed by the δ-axis component of the measured current, or
a comparison between the estimated δ-axis current and a predetermined value.

6. The motor control device of claim 1, wherein
let an axis parallel to a magnetic flux produced by a permanent magnet provided on the rotor be a d-axis,
let an axis estimated, for control purposes, to correspond to the d-axis be a γ-axis, and
let an axis estimated, for control purposes, to be perpendicular to the γ-axis be a δ-axis, then
the first state quantity is represented by a specified δ-axis voltage value to be followed by a δ-axis component of the motor voltage,
the synchronization failure detector estimates as the second state quantity a second estimated rotation speed, separate from the first estimated rotation speed, based on
the first state quantity alone or
the first state quantity and the third state quantity represented by a measured current of the motor current or the specified current value to be followed by the measured current, and

the synchronization failure detector detects synchronization failure based on the second estimated rotation speed.

7. A motor drive system comprising:
a motor;
an inverter driving the motor; and
the motor control device of claim 1, the motor control device controlling the motor by controlling the inverter.