Method and apparatus of controlling an induction motor

Abstract: A three-phase-to-two-phase converter (10) detects two-phase currents of three-phase stator currents, converts the stator currents into a γ - δ coordinate system to derive a γ-axis current and a δ-axis current. The γ-axis current, the δ-axis current, and voltage command values outputted respectively from a γ-axis current controller (2) and a δ-axis current controller (3) are entered to a stator current and induced electromotive force estimator (9), which determines estimated values of the stator currents and estimated values of an induced electromotive force in the γ - δ coordinate system.

Description: TECHNICAL FIELD The present invention relates to a method of controlling an induction motor by bringing γ - δ coordinate axes which are reference coordinate axes of a control system configuration into agreement with d - q coordinate axes where the d-axis represents the direction of the vector of secondary magnetic fluxes of the induction motor and the q-axis represents an axis perpendicular to the d-axis. BACKGROUND ART Processes for estimating the speed of an induction motor using the induced electromotive force (induced voltage) are employed in industrial applications for the purpose of making the V/f control highly accurate and simplifying calculations for the speed-sensor-less vector control because the processes are physically obvious from the fact that an induced electromotive force is generated when the rotor rotates and also because the processes are able to calculate the rotational angular velocity of secondary magnetic fluxes with a relatively small amount of calculations. Examples of conventional processes include [1] "Study of making speed-sensor-less vector control highly responsive", 1994 Electric Society National Conference No. 655, and [2] "Speed-sensor-less vector control in view of weak field region", 1994 Electric Society National Conference No. 656. The former process pays attention to a voltage model of the q-axis of a d - q coordinate system and proposes a simple disturbance estimator for estimating a q-axis induced electromotive force using a q-axis voltage command value as an input and a q-axis current as an output. The simple disturbance estimator estimates a speed by multiplying the estimated q-axis induced electromotive force by coefficients representing a motor constant and a flux command value. The former process is an improvement over a process that uses the output of a q-axis current controller as an estimated speed. The former process achieves separation between a speed control system and a current control system by introducing the estimation of an induced electromotive force, and is successful in improving the speed control response. The latter process estimates an induced electromotive force from a measured voltage and a measured current using a voltage equation of the induction motor, estimates a power supply angular frequency from the estimated induced electromotive force and an estimated value of secondary magnetic fluxes, and subtracts a command value for a slip frequency from the estimated power supply angular frequency to estimate a speed. The above conventional processes basically use voltage and current equations in the d - q coordinate system of an induction motor, i.e., where εd: a d-axis component of the induced electromotive force, εq: a q-axis component of the induced electromotive force, Φrd: a d-axis component of the secondary magnetic fluxes, Φrq: a q-axis component of the secondary magnetic fluxes, Isd: a d-axis component of the stator current (excitation-related current), Isq: a q-axis component of the stator current (torque-related current), usd: a d-axis component of the stator voltage, usq: a q-axis component of the stator voltage, Rs: a stator resistance, Lm: mutual inductance, Lr: rotor self-inductance, Ls: stator self-inductance, δ = 1 - Lm 2 /(LsLr), ω*: a power supply angular frequency. The conventional processes estimate an induced electromotive force according to the equation (1), and estimate a speed according to the equation (2). Since a current differential term in the right-hand side of the equation (1) is omitted, the conventional processes fail to estimate a correct induced electromotive force while the current is varying. If the current differential term is not omitted, then noise contained in a measured current value is amplified. The literature [1] proposes a process of equivalently canceling a low-pass filter used to output an estimated speed and a current differential term. The literature [1] assumes an ideal vector control state, estimates only a q-axis induced electromotive force, and estimates a speed using the estimated q-axis induced electromotive force. However, there is a possibility that the ideal vector control state may collapse due to a change in the induction motor constant and errors of measured current and voltage with respect to actual current and voltage. The literature [1] is silent about how to compensate for such a collapse of the ideal vector control state. The literature [2] proposes the following equation (3) in order to compensate for a collapse of the ideal vector control state: ω* = sqn(εq)(|εq|Φrd - KΦrd*εd) In a steady state, the existence of εd signifies the existence of Φrq as can be seen from the equation (2), and indicates that the ideal vector control state has collapsed. Specifically, the γ - δ coordinate axes which are reference coordinate axes of the control system configuration and the d - q coordinate axes are displaced from each other. Axis-displacement compensation is carried out by eliminating εd in the second term of the right-hand side of the equation (3). However, the equation is contradictory in that ω* is necessary to estimate the induced electromotive force used in the equation (3) and ω* is necessary to determine ω*. In order to actually perform software-implemented calculations of the equation, ω* at a preceding time in a discrete system is used. If ω* varies quickly, then the speed response tends to be delayed. Finally, HARDAN F ET AL: 'Adaptive current regulation using extended Kalman filter for vector-controlled induction machine drive' POWER ELECTRONICS SPECIALISTS CONFERENCE, 1995. PESC '95 RECORD., 26TH ANNUAL IEEE ATLANTA, GA, USA 18-22 JUNE 1995, NEW YORK, NY, USA, IEEE, US, 18 June 1995 (1995-06-18), pages 457-463, XP010150596 ISBN: 0-7803-2730-6 discloses a model-based predictive current regulator for a vector-controlled induction motor drive. An extended Kalman filter is employed to estimate the back EMF of the machine stator circuit model used by the inner current regulator loop. This also estimates the rotor resistance for slip angular velocity calculation in the outer speed loop. DISCLOSURE OF THE INVENTION It is an object of the present invention to provide a method of controlling an induction motor according to claim 1, and an apparatus for controlling an induction motor according to claim 3. A method according to the present invention comprises estimating an induced electromotive force of an induction motor by detecting at least two of three-phase stator currents of the induction motor, converting the detected stator currents into a γ - δ coordinate system to derive a γ-axis current Isγ and a δ-axis current Isδ, constructing a state estimator using differences Isγ - Îsγ, Isδ - Îsδ between estimated values Îsγ, Îsδ of the γ-axis current and the δ-axis current in a preceding control loop and actual values thereof as corrective quantities, using stator voltage command values Usγ*, Usδ* in the γ - δ coordinate system as inputs, and using an induced electromotive force generated by rotation of a rotor as a disturbance, and estimating estimated values Îsγ, Îsδ of stator currents and estimated values ε and γ , ε and δ of the induced electromotive force in the γ - δ coordinate system. A method according to the present invention further comprises estimating a speed of an induction motor by adding the product of the difference Imrγ - Isγ between an estimated value or a command value Imry of a γ component of a secondary field current and a γ component Isγ of a stator current and a motor constant to a γ-axis induced electromotive force ε and γ in the γ - δ coordinate system which has been estimated by the above method, thereby to calculate a new induced voltage Eγest, adding the product of the difference Imrδ - Is5 between an estimated value or a command value Imrδ of a δ component of the secondary field current and a δ component Isδ of the stator current and the motor constant to a δ-axis induced electromotive force ε and δ in the γ - δ coordinate system which has been estimated by the above method, thereby to calculate a new induced voltage Eδest, estimating the magnitude of a rotor speed ω γ of the induction motor using the sum of the squares of the induced voltages Eγest, Eδest and an estimated value or a command value of secondary magnetic fluxes, determining the sign of the rotor speed with the sign of the induced voltage Eδest, and estimating an estimated value ω and γ of the rotor speed. A method according to the present invention further comprises correcting an axis displacement of an induction motor by subtracting the product of the induced voltage Eyest and a gain from the speed value estimated by the above method to adjust the estimated speed value in a direction to eliminate the induced voltage Eγest for thereby bringing γ - δ coordinate axes into agreement with d - q coordinate axes. According to the present invention, there is provided an apparatus for controlling an induction motor according to claim 3. The apparatus may further comprise : a speed controller for being supplied with the difference between a speed command value ω γ * and an estimated speed value ω and γ as a first difference and outputting a δ-axis current value Isδ*; a δ-axis current controller for being supplied with the difference between a δ-axis current command value Isδ* and a δ-axis current Isδ or a δ-axis current estimated value Îsδ and outputting a δ-axis voltage command value Usδ*; a γ-axis current controller for being supplied with the difference between a command value Imr* of a secondary field current related to the magnitude of secondary magnetic fluxes and a γ-axis current Îsγ or a γ-axis current estimated value Îsγ as a second difference and outputting a γ-axis voltage command value Usγ*; a polar coordinate converter for converting the voltage command values Usδ*, Usγ* into magnitudes and phases of the voltage command values; a slip angular frequency calculator for being supplied with the first difference and the second difference and determining a slip angular frequency command value ω sl *; an integrator for converting a rotational speed of the γ - δ coordinate axes determined by adding the estimated speed value ω and γ and the slip angular frequency command value ω sl * to each other into an angle; an inverter circuit for being supplied with an output from the polar coordinate converter and the position of the rotated γ - δ coordinate axes outputted from the integrator, and generating a PWM signal to energize the induction motor; The apparatus for controlling the induction motor according to the present invention also has a speed estimator for being supplied with the estimated values ε and γ , ε and δ of the induced electromotive force and carrying out the above method of estimating the speed. The apparatus for controlling the induction motor according to the present invention also has an axis displacement corrector for carrying out the above method of correcting the axis displacement. The method of estimating the induced electromotive force of the induction motor according to the present invention is constructed of a state estimator expressed by the following equation (4) based on the theoretical equation (1) of the induction motor: In the equation (4), an output error feedback gain for multiplying the poles of the state estimator by Isγ - Isγ, Isδ - Isδ can be set up arbitrarily. Therefore, the equation (4) can estimate an induced electromotive force in a manner insensitive to changes in the inductor motor constant, and adjust the time in which the estimated values are converged. When the equation (4) is used in an actual control process, it becomes a differential equation converted to a discrete system. Therefore, unlike the equation (1), it is not necessary to calculate a differential of the current for estimating an induced electromotive force, and measurement noise is not amplified. The state estimator is constructed on the γ - δ coordinate axes which are rotating at the angular velocity ω*. Since estimated values are DC components, changes in their errors with respect to actual values on the d - q coordinate axes are slow. If the poles are stably established, then it is not necessary to increase the poles unnecessarily to speed up the convergence. In the method of estimating the speed of the induction motor according to the present invention, when the following current equation (5) for the rotor: is substituted for the flux differential term in the first term of the right-hand side of the equation (2) to eliminate the differential term of the secondary magnetic fluxes, the following equation (6) is obtained: However, the equation is developed on the γ - δ coordinate axes which are control reference coordinate axes. If induced electromotive force estimated values are used in the equation (6), then In the equation, new induced voltages are defined as Eγest, Eδest. The new induced voltages Eγest, Eδest are calculated by the right-hand side of the first row of the equation (7). Based on the equation (7), a speed is estimated according to the following equation (8): | ω r | = Lr Lm E γ est 2 + E δ est 2 Φ r γ 2 + Φ r δ 2 The sign of the speed is determined by the sign of Eδest. Thus, ωr = sgn(Eδest)|ωr| where sgn(Eδest) represents the sign of Eδest. The speed is finally estimated according to the equation (9). For the secondary field currents Imry, Imrδ and the secondary fluxes Φrγ, Φrδ in the equations (7), (8), there is used a magnetic flux estimated value or a magnetic flux command value based on the equation (5). If a magnetic flux command value is used, then, assuming that the command value for the secondary magnetic fluxes is represented by Φr*, the speed is estimated by: ωr = sgn(Eδest)|ωr| - KcomEγest where Kcom represents a gain. If the secondary magnetic flux command value is used as described above, it is simultaneously necessary to perform a compensation process to bring the γ - δ coordinate axes and the d - q coordinate axes into agreement with each other. In the method of compensating for the axis displacement according to the present invention, since Eγest is a function of a δ component Φrδ of secondary magnetic fluxes generated when the γ - δ coordinate axes and the d - q coordinate axes are displaced from each other, as indicated by the equation (7), the speed command value is adjusted in a direction to eliminate Eγest for thereby compensating for the axis displacement. Specifically, the equation (9) is expanded into the equation (10). Φrγ = Φr*, Imrγ = Φr*/Lm, Φrδ = Imrδ = 0 In the method of estimating the speed according to the present invention, inasmuch as the speed is directly estimated using the induced electromotive force estimated values, the contradiction of the equation development according to the conventional processes is eliminated, reducing a control delay. In the method of estimating the speed according to the present invention, the product of Eγest and the gain is subtracted from the estimated speed to adjust the estimated speed value in a direction to eliminate Eγest for thereby bringing the γ - 5 coordinate axes and the d - q coordinate axes into agreement with each other. Eγest is composed of the electromotive force ε and γ estimated by the state estimator and the difference between the secondary field current and the excitation-related component of the stator current, and takes into account variations of the secondary magnetic fluxes. While the secondary magnetic fluxes are assumed to be constant according to the conventional processes, the method according to the present invention is free of such an assumption, and poses no problem even when the magnetic fluxes vary at the start of excitation when the induction motor is activated and at the time of field-weakening control. BRIEF DESCRIPTION OF THE DRAWINGS The figure is a block diagram of an apparatus for controlling an induction motor according an embodiment of the present invention. BEST MODE FOR CARRYING OUT THE INVENTION FIG. 1 is a block diagram of an apparatus for controlling an induction motor according an embodiment of the present invention. The apparatus according to the embodiment serves to control an induction motor 11, and comprises a speed controller 1, a γ-axis current controller 2, a δ-axis current controller 3, a polar coordinate converter 4, an inverter circuit 5, a slip angular frequency calculator 6, a speed estimator and axis displacement corrector 7, an integrator 8, a stator current and induced electromotive force estimator 9 on γ - δ coordinate axes, and a three-phase-to-two-phase converter 10. The speed controller 1 is supplied with the difference between a speed command value ω γ * and an estimated speed value ω and γ and outputs a δ-axis current value Isδ*. The δ-axis current controller 3 is supplied with the difference between a δ-axis current command value Isδ* and a δ-axis current Isδ and outputs a δ-axis voltage command value Usδ*. The γ-axis current controller 2 is supplied with the difference between a command value Imr* (= Imrγ = Φr*/Lm) of a secondary field current related to the magnitude of secondary magnetic fluxes and a γ-axis current estimated value Isγ and outputs a γ-axis voltage command value Usγ*. The polar coordinate converter 4 converts the voltage command values Usδ*, Usγ* into magnitudes and phases of the voltage command values. The inverter circuit 5 is supplied with an output from the polar coordinate converter 4 and the position of the rotated γ - δ coordinate axes outputted from the integrator 8, and generates a PWM signal to energize the induction motor 11. The three-phase-to-two-phase converter 10 detects three-phase currents from stator currents Isu, Isv and converting the detected three-phase currents into two-phase alternating currents. The stator current and induced electromotive force estimator 9 on the γ - δ coordinate axes is supplied with Isγ, Isδ subjected to coordinate transformation using the position of the γ - δ coordinate axes from the two-phase alternating currents and the voltage command values Usγ*, Usδ*, performs calculations thereon according to the equation (4), and outputs estimated values Îsγ, Îsδ of the stator currents and estimated values ε and γ , ε and δ of an induced electromotive force in the γ - δ coordinate system. The speed estimator and axis displacement corrector 7 is supplied with the estimated values ε and γ , ε and δ of the induced electromotive force, performs calculations according to the equations (7), (8), (10), and outputs an estimated speed value ω and γ . The slip angular frequency calculator 6 performs calculations according to the equation (12) to determine a slip angular frequency command value ω sl *. ω sl * = Rr · Is δ* Lr ·Im r * The integrator 8 converts a rotational speed of the γ - δ coordinate axes determined by adding the estimated speed value ω γ and the slip angular frequency command value ω s1 * to each other into an angle, and outputs the angle to the inverter circuit 5 and the three-phase-to-two-phase converter 10. If the present embodiment is implemented by software, then all the calculations are discretely approximated in use.

Claims: A method of controlling an induction motor having γ - δ coordinate axes which are reference coordinate axes of a control system configuration and d - q coordinate axes where the d-axis represents the direction of the vector of secondary magnetic fluxes of the induction motor and the q-axis represents an axis perpendicular to the d-axis, the method comprising the steps of: detecting at least two of three-phase stator currents of the induction motor; converting the detected stator currents into a γ - δ coordinate system to derive a γ-axis current Isγ and a δ-axis current Isδ; estimating values of stator currents Îsγ and Îsδ and of induced electromotive forces ε and γ and ε and δ in the γ - δ coordinate system with a state estimator using differences Isγ - Îsγ, Isδ - Îsδ between estimated values Îsγ, Îsδ of the γ-axis current and the δ-axis current in a previous estimation and actual values thereof as corrective quantities, stator voltage command values Usγ*, Usδ* in the γ - δ coordinate system as inputs, and induced electromotive forces ε and γ and ε and δ generated by rotation of a rotor as a disturbance; characterized in the following steps adding the product of the difference Imrγ - Isγ between an estimated value or a command value Imrγ of a γ component of a secondary field current and a γ component Isγ of a stator current and a motor constant to the estimated γ-axis induced electromotive force ε and γ in the γ - δ coordinate system, thereby to calculate a new induced voltage Eγest; adding the product of the difference Imrδ - Isδ between an estimated value or a command value Imrδ of a δ component of the secondary field current and a δ component Isδ of the stator current and the motor constant to the estimated δ-axis induced electromotive force ε and δ in the γ - δ coordinate system, thereby to calculate a new induced voltage Eδest; estimating the magnitude of a rotor speed ω γ of the induction motor using the sum of the squares of the induced voltages Eγest, Eδest and an estimated value or a command value of secondary magnetic fluxes; determining the sign of the rotor speed with the sign of the induced voltage Eδest so as to estimate an estimated value ω and γ of the rotor speed; and correcting an axis displacement by subtracting the product of the induced voltage Eyest and a gain from the estimated speed value to adjust the estimated speed value in a direction to eliminate the induced voltage Eγest, thereby bringing the γ - δ coordinate axes into agreement with the d - q coordinate axes. A method according to claim 1, further comprising the steps of: determining a δ-axis current value Isδ* with reference to the difference between a speed command value ω γ * and the estimated speed value ω and γ as a first difference; determining a δ-axis voltage command value Usδ* with reference to the difference between a δ-axis current command value Isδ* and a δ-axis current Isδ or a δ-axis current estimated value ísδ; determining a γ-axis voltage command value Usγ* with reference to the difference between a command value Imr* of a secondary field current related to the magnitude of secondary magnetic fluxes and a γ-axis current Isγ or a γ-axis current estimated value Îsγ as a second difference; converting the volage command values Usδ*, Usγ* into magnitudes and phases of the voltage command values; determining a slip angular frequency command value ω sl * with reference to the first difference and the second difference; converting a rotational speed of the γ - δ coordinate axes determined by adding the estimated speed value ω and γ and the slip angular frequency command value ω sl * to each other into an angle; generating a PWM signal to energize the induction motor with reference to the angle and the magnitudes and phases of the voltage command values; An apparatus for controlling an induction motor having γ - δ coordinate axes which are reference coordinate axes of a control system configuration and d - q coordinate axes where the d-axis represents the direction of the vector of secondary magnetic fluxes of the induction motor and the q-axis represents an axis perpendicular to the d-axis, comprising means for detecting at least two of three-phase stator currents of the induction motor; means for converting the detected stator currents into a γ - δ coordinate system to derive a γ-axis current Isγ and a δ-axis current Isδ; a state estimator for estimating values of stator currents Îsγ and Îsδ and of induced electromotive forces ^ε γ and ε and δ in the γ - δ coordinate system, adapted to use differences Isγ - Îsγ, Isδ - Îsδ between estimated values Îsγ, Îsδ of the γ-axis current and the δ-axis current in a previous estimation and actual values thereof as corrective quantities, stator voltage command values Usγ*, Usδ* in the γ - δ coordinate system as inputs, and an induced electromotive forces ε and γ and ε and δ generated by rotation of a rotor as a disturbance; means for adding the product of the difference Imrγ - Isγ between an estimated value or a command value Imrγ of a γ component of a secondary field current and a γ component Isγ of a stator current and a motor constant to the estimated γ-axis induced electromotive force ε and γ in the γ - δ coordinate system, thereby to calculate a new induced voltage Eγest; means for adding the product of the difference Imrδ - Isδ between an estimated value or a command value Imrδ of a δ component of the secondary field current and a δ component Isδ of the stator current and the motor constant to the estimated δ-axis induced electromotive force ε and δ in the γ - δ coordinate system, thereby to calculate a new induced voltage Eδest; means for estimating the magnitude of a rotor speed ω γ of the induction motor using the sum of the squares of the induced voltages Eγest, Eδest and an estimated value or a command value of secondary magnetic fluxes; means for determining the sign of the rotor speed with the sign of the induced voltage Eδest so as to estimate an estimated value ω and γ of the rotor speed; and means for correcting an axis displacement by subtracting the product of the induced voltage Eyest and a gain from the estimated speed value to adjust the estimated speed value in a direction to eliminate the induced voltage Eγest, thereby bringing the γ - δ coordinate axes into agreement with the d - q coordinate axes. An apparatus according to claim 3, further comprising: a speed controller for being supplied with the difference between a speed command value ωγ* and an estimated speed value ω and γ as a first difference and outputting a δ-axis current value Isδ*; a δ-axis current controller for being supplied with the difference between a δ-axis current command value Isδ* and a δ-axis current Isδ or a δ-axis current estimated value of Îsδ and outputting a δ-axis voltage command value Usδ*; a γ-axis current controller for being supplied with the difference between a command value Imr* of a secondary field current related to the magnitude of secondary magnetic fluxes and a γ-axis current Isγ or a γ-axis current estimated value Îsγ as a second difference and outputting a γ-axis voltage command value Usγ*; a polar coordinate converter for converting the voltage command values Usδ*, Usγ* into magnitudes and phases of the voltage command values; a slip angular frequency calculator for being supplied with the first difference and the second difference and determining a slip angular frequency command value ω sl *; an integrator for converting a rotational speed of the γ - δ coordinate axes determined by adding the estimated speed value ω and γ and the slip angular frequency command value ω sl * to each other into an angle; an inverter circuit for being supplied with an output from the polar coordinate converter and the position of the rotated γ - δ coordinate axes outputted from the integrator, and generating a PWM signal to energize the induction motor;

Inventors: IDE KOZO, MORIMOTO SHINYA, OGURO RYUICHI, IDE, KOZO, MORIMOTO, SHINYA, OGURO, RYUICHI, IDE, Kozo, K.K. Yaskawa Denki, MORIMOTO, Shinya, K.K. Yaskawa Denki, OGURO, Ryuichi, K.K. Yaskawa Denki
Applicants: YASKAWA DENKI SEISAKUSHO KK, KABUSHIKI KAISHA YASKAWA DENKI
Application year: 1998
Main classification: H02P21
IPC codes: H02P21, H02P27