Digital rotor flux observer

Information

  • Patent Grant
  • 6509711
  • Patent Number
    6,509,711
  • Date Filed
    Wednesday, April 26, 2000
    24 years ago
  • Date Issued
    Tuesday, January 21, 2003
    21 years ago
Abstract
A system (10) for controlling the torque of an induction motor (12) utilizes a flux observer (14). The flux observer (14) receives stator current inputs (16, 18) and a rotor speed estimate (80) and then outputs a rotor flux estimate (20) that provides increased motor control stability at all speeds.
Description




TECHNICAL FIELD




The present invention relates to a system for controlling an electric induction motor and, more particularly, to rotor flux observer design for controlling an electric induction motor.




BACKGROUND OF THE PRESENT INVENTION




Typical induction motors include a rotor mounted inside a stator for rotation about a rotation axis. A common rotor design includes a “squirrel cage winding” in which axial conductive bars are connected at either end by shorting rings to form a generally cylindrical structure. The stator is formed by a plurality of winding which surround the rotor and are typically arranged in three separate phases. Time varying voltage is applied across the stator windings which generate an electromotive force (emf) and an associated stator magnetic field which rotates around the stator at a stator field frequency primarily within the space defined by the rotor.




As the stator field rotates about the rotor, relative motion between the stator field flux and the rotor bars induces voltages in the rotor at a slip frequency. Slip is the difference between the stator frequency. Slip is the difference between the stator field frequency and the rotor electrical speed. The voltages induced in the rotor cause rotor bar currents, which in turn generate a rotor magnetic field. The stator and rotor fields are stationary with respect to each other, but are separated by a changeable rotor angle (a). The two fields interact to produce torque, which causes rotor rotation. The amount of motor torque for an induction motor will vary as the voltages induced in the rotor change.




It is known that for efficient control of an induction motor, the rotor position must be measured or estimated. A rotor position sensor is typically employed to sense the position of the rotor and feed it directly to a controller or other device that is controlling the operation of the motor. Based on a reading from the position sensor, in field oriented control of an induction motor, the orientation of rotor flux is determined and the stator voltage signal is calculated. While position sensors have been widely used for sensing the rotor position, they can be unreliable, are relatively noisy, and are relatively expensive.




Various control schemes for torque control of induction motors (IM) without the use of a position sensor have been employed. These IM torque control schemes include as a part the flux estimation procedure, the rotor flux observer, which is based on the following known rotor flux dynamic equation:














λ
^

r




t


=



(



-


R
r


L
r




I

+


ω
r


J


)




λ
^

r


+



R
r


L
r



M






i
s







(
1
)





















Where







λ
^

r


=

[



λ
^

rx



λ
^

ry


]


,


i
s

=

[


i
xs


i
ys


]


,

I
=

[



1


0




0


1



]


,

J
=

[



0



-
1





1


0



]






(
2
)













In the equations above, {circumflex over (λ)}


r


is the rotor flux estimation, i


s


is the stator current, and ω


r


is the rotor electrical speed and R


r


, L


r


, and M are the rotor resistance, rotor inductance, and machine mutual inductance, respectively.




For higher rotor speeds, the discrete implementation of equation (1), shown below:











λ
^


rk
+
1


=



λ
^

rk

+

T


[



(



-


R
r


L
r




I

+


ω
rk


J


)




λ
^

rk


+



R
r


L
r



M






i
sk



]







(
3
)













where T is a sampling period, can be numerically unstable. The solution of (1) is oscillatory with frequency ω


r


, and the simple Euler integration method implemented in equation (3) is not adequate for high rotor speeds. Thus, an improved solution for controlling the motor torque that operates reliably at all speeds would be advantageous.




SUMMARY OF THE INVENTION




It is an object of the present invention to provide a discrete time flux observer for torque control of induction motors.




It is another object of the present invention to provide a flux observer that provides significantly improved performance at high rotor speeds.




It is yet another object of the present invention to provide a flux observer that provides an accurate flux estimation that provides increased endurance and ease of maintenance preserving high dynamic performance.




In accordance with the above and other objects of the present invention, a system for controlling the torque of an induction motor is provided. The system includes a controller that receives a signal which is representative of the motor stator current. The controller also receives information on the rotor speed. The actual rotor speed or the estimated rotor speed is then fed to a flux observer. The flux observer calculates the rotor flux which is necessary to calculate the stator voltage to control the output torque of the induction motor. The observer is a discrete time flux observer that utilizes analytical solutions of the rotor flux dynamic equation over a time sampling interval and thus provides higher accuracy and improved performance over standard Euler or Tustin approximations.




These and other features of the present invention will become apparent from the following description of the invention, when viewed in accordance with the accompanied drawings and appended claims.











BRIEF DESCRIPTION OF THE DRAWINGS





FIG. 1

is a block diagram of a sensorless field orientation control scheme in accordance with a preferred embodiment of the present invention;





FIG. 2

is a block diagram of a discrete time flux observer for torque control of an induction motor in accordance with a preferred embodiment of the present invention;





FIG. 3

is a graph illustrating an estimation of one component of the rotor flux using the rotor speed estimated value in accordance with a preferred embodiment of the present invention; and





FIG. 4

is a graph illustrating an estimation of the rotor flux with a fifty percent (50%) variation in the rotor time constant in accordance with a preferred embodiment of the present invention.











BEST MODE(S) FOR CARRYING OUT THE INVENTION




Referring now to

FIG. 1

, which illustrates a sensorless field orientation control system


10


in accordance with the present invention. The components of the control system


10


, which are used to control the torque of an induction motor


12


, are briefly and generally described below, as the basic configuration is well known. It should be understood that the control system


10


is merely exemplary and a variety of other components and configurations may be utilized. Further, the induction motor


12


can be incorporated into a variety of applications, including both automotive and non-automotive. An example of an automotive application in which the induction motor


12


can be used is in an electronic power assisted steering system.




For example, as shown in

FIG. 1

, a flux observer


14


is provided with two current inputs. The first current input (i


sx


)


16


is the x component of the stator current while the second current input (i


sy


)


18


is the y component of the stator current. The flux observer


14


outputs the rotor flux (λ


r


) which is determined as set forth below and is generally indicated by line


20


. Thereafter, the rotor flux angle (θ


r


) is calculated, as generally represented by reference number


22


, and is output as generally indicated by line


24


. The output


24


, indicative of the rotor flux angle, is fed to a field oriented controller


26


and is also used to assist in the synchronous to stator transformation, which is generally indicated by reference number


28


.




The field oriented controller


26


determines the stator voltage in the frame related to the rotor magnetic field and controls the (d, q) components of stator current according to standard techniques which are well known in the art. The field oriented controller


24


is provided with various inputs, including the reference value of the rotor flux (λ


r


*), as generally indicated by line


21


, the motor torque reference signal (T


e


*), as generally indicated by line


33


, the stator x current input (i


sx


)


16


, and the stator y current input (i


sy


)


18


. The field oriented controller


26


outputs a first voltage component output (u


d


), as generally indicated by reference number


34


, and a second voltage component output (u


q


), as generally indicated by reference number


36


. These voltage outputs


34


,


36


are transformed from the synchronous to stator frame in block


28


.




The synchronous to stator transformation, as indicated by reference number


28


, outputs an x component of the voltage (u


x


) as generally indicated by reference number


38


and a y component of the voltage (u


y


) is generally indicated by reference number


40


. The voltage components


38


and


40


are then fed into a two phase to three phase converter


42


. The two phase to three phase converter transforms the x component of voltage (u


x


)


38


and the y component of voltage (u


y


)


40


into three voltage phases represented by (u


a


) as indicated by reference number


44


, (u


b


) as indicated by reference number


46


, and (u


s


) as indicated by reference number


48


. The three voltage phases (u


a


), (u


b


) and (u


s


), are input into a phase wave modulator voltage source inverter


50


, which inputs the three voltage phases (u


a


), (u


b


) and (u


s


), into corresponding stator windings and the three stator current signals (i


sa


, i


sb


, i


sc


) are measured and generally indicated by reference numbers


52


,


54


and


56


, respectively.




Additionally, the x component of voltage (u


x


)


38


and the y component of voltage (u


y


)


40


are also fed to a voltage observer


58


, which then outputs two voltage components (u


o


) as indicated by reference number


60


and (u


l


) as indicated by reference number


62


, which are then fed to the speed observer


64


. The current signals from the stator windings (i


sa


, i


sb


and i


sc


) are then input into a three phase to two phase converter


66


which generates two current outputs (i


sx


)


68


and (i


sy


)


70


. The current output


68


,


70


are then fed into a current observer


72


. The current observer generates current output (i


so


), as indicated by reference number


74


, (i


s1


) as indicated by reference number


76


, and (i


s2


) as indicated by reference number


78


. The current outputs


74


,


76


,


78


are then fed into the speed observer


64


. The speed observer


64


outputs a signal


80


which represents the rotor speed estimation ({circumflex over (ω)}


r


). The rotor speed estimation ({circumflex over (ω)}


r


) is then fed into the flux observer


14


.




Referring now to

FIG. 2

, which schematically illustrates the digital rotor flux observer


14


in accordance with the present invention. The digital rotor flux observer


14


is designed based upon a zero order hold (ZOH) approach and analytical solutions of the flux dynamic equation. In order to overcome the instability experienced by prior flux observers that are based on standard approximations, equation (1) set forth above is transformed, in accordance with the present invention, to discrete time using ZOH approximation.




Straightforward implementation of ZOH approximation requires a set of discrete models which are calculated off line over the range of the rotor speed. This requires high memory resources and essential work on preliminary observer calibration. To avoid this, an analytical solution for the flux observer design is suggested. The design of the flux observer


14


is based on the scalar form of the rotor dynamic equation using complex variables for flux and current. Thus, rewriting equation (1) in scalar form with complex variables for flux and current results in the following discrete time equation:














λ
r




t


=



(



-
η

+

j





ω


,

)



λ
r


+

η





M






i
s






where











λ
r

=


λ
rx

+

j






λ
ry




,


i
s

=


i
sx

+

j






i
sy




,

η
=




R
r


L
r







and





j

=



(

-
1

)


.








(
4
)













By using the analytical solution of equation (4) above, the flux observer


14


is preferably designed in accordance with the following equation over the sampling interval [kT,(k+1)T]:










λ

k
+
1


=



λ
k







kT


(

k
+
1

)


T





[

η
+

j






ω


(
s
)




]








s





+

η





M




kT


(

k
+
1

)


T





i


(
t
)








kT
l




[

η
+

j






ω


(
s
)




]








s











t









(
5
)













In accordance with the present invention, it is assumed that the rotor speed ω(s) is constant over the sampling interval set forth in equation (5). Thus, the rotor speed estimation ({circumflex over (ω)}


r


)


80


is input into the flux observer


14


. Given that the speed dependant exponent can be written as a constant σ=η−jω the exponential in equation (5) can thus be expressed as e


−στ


, as shown in FIG.


2


. Equation (5) can therefore be simplified, as the rotor speed is assumed to be constant over a sampling interval, and expressed in the form of equation (6) below:











λ
r



(
t
)


=



λ
ro





-

σ


(

t
-

t
o


)





+

η





M




to
t






σ


(

τ
-
t

)





i
s








τ









(
6
)













Accordingly, by using equation (6), the discrete form of the flux observer, which includes the speed as a parameter, is then calculated as follows:










λ

rk
+
1


=



λ
k






-
σ






T



+

η





M






i
sk


1

-





-
σ






T


σ






(
7
)













The discrete form of the flux observer is represented generally by reference number


84


. The variable T is the sampling period. After recovering the scalar components of the rotor flux from equation (7), the analytical digital form of the flux observer


14


is calculated, as set forth below:











λ

rx


(

k
+
1

)



=



(

1
-
α

)



λ
rxk


-

βλ
ryk

+



η





M



η
2

+

ω
r
2





[



(

ηα
+


ω
r


β


)



i
sx


-


(



ω
r


α

-
ηβ

)



i
sy



]











λ

ry


(

k
+
1

)



=



(

1
-
α

)



λ
ryk


-

βλ
rxk

+



η





M



η
2

+

ω
r
2





[



(

ηα
+


ω
r


β


)



i
sy


-


(



ω
r


α

-
ηβ

)



i
sx



]











where





α

=


1
-





-
η






T



cos






ω
r


T





β


=





-
η






T



sin






ω
r


T







(
8
)













In accordance with the present invention, equations (8) are the discrete ZOH approximation of the continuous equation (1). The rotor speed ({circumflex over (ω)}


r


)


80


and the stator current (i


s


)


16


,


18


are parameters that are inputs into the flux observer


14


to arrive at the discrete time instance. The rotor flux angle ({circumflex over (θ)}


r


) is preferably determined at


22


in accordance with the equation λ


ry


/({square root over (λ


rx





2





ry





2


)})=cos{circumflex over (θ)}


r


and then input to the field oriented controller


26


.




The plots in

FIGS. 3 and 4

show the result of the estimation of the rotor flux {circumflex over (λ)}


r


for the modeling of engine idle active damping using an induction motor (IM) with parameters given in Table 1 below:




















R


s


(ohm)




R


r


(ohm)




L


r


(H)




L


s


(H)




M (H)











η


(

1
sec

)


















.15




.0187




.0011




.0011




.00081




17















FIG. 3

plots the flux on the y-axis versus time in seconds on the x-axis.

FIG. 3

illustrates the estimation of one component of the rotor flux as determined by the flux observer


14


in accordance with equations (8) and a sampling frequency of IKHZ using the rotor speed estimated value. The flux observer of the present invention that operates in accordance with equation (8) demonstrate good performance, while the estimation based on the standard algorithm (3) is unstable.





FIG. 4

shows that the rotor flux observer


14


of the present invention is robust to motor parameter variations. For example,

FIG. 4

illustrates a flux estimator for a fifty percent (50%) variation of the rotor time constant. As shown, the flux estimate is higher during transient. However, the flux phase, which is important for the control, is not affected by variations in the rotor time constant




The flux observer


14


of equations of the present invention, thus provides an accurate flux estimator at all speeds, including at high speeds where solution is oscillatory and prior estimations failed. The disclosed flux observer


14


is robust with respect to motor parameter variations. For example, varying the rotor time constant changes the amplitude of flux estimate during transient, but the phase of the flux estimate and the rotor angle estimate remain essentially unchanged. Thus, implementation of the disclosed flux observer


14


eliminates the need for a position sensor, reduces cost and improves reliability while preserving the system high dynamic performance.




Having now fully described the invention, it will be apparent to one of ordinary skill in the art that many changes and modifications can be made thereto without departing from the spirit or scope of the invention as set forth herein.



Claims
  • 1. A method of utilizing a rotor flux observer to control the torque of an induction motor, having a rotor and a stator, comprising:receiving a signal indicative of the stator current at the rotor flux observer; receiving a signal representative of the rotor speed at the rotor flux observer; calculating an exponent of said rotor speed signal; estimating the rotor flux based on the assumption that the rotor speed is constant over a sampling interval determining a discrete form of the rotor flux; estimating a rotor flux angle; calculating a voltage to be used in controlling the motor torque, based on a field orientation scheme.
  • 2. A method of utilizing a rotor flux observer to control the torque of an induction motor, having a rotor and a stator, comprising:receiving a signal indicative of the stator current at the rotor flux observer; receiving a signal representative of the rotor speed at the rotor flux observer; calculating an exponent of said rotor speed signal; estimating the rotor flux based on the assumption that the rotor speed is constant over a sampling interval; determining a discrete form of the rotor flux; recovering scalar components of the rotor flux to obtain a flux estimated volume, wherein said scalar components of the rotor flux are determined according to the following equations: λrx⁢(k+1)=(1-α)⁢λrxk-βλryk+η⁢ ⁢Mη2+ωr2⁡[(ηa+ωr⁢β)⁢isx-(ωr⁢a-ηβ)⁢isy]λry⁢(k+1)=(1-α)⁢λryk+βλrxk+η⁢ ⁢Mη2+ωr2⁡[(ηa+ωr⁢β)⁢isy-(ωr⁢a-ηβ)⁢isx]estimating a rotor flux angle; calculating a voltage to be used in controlling the motor torque, based on a field orientation scheme.
  • 3. The method of claim 1, wherein said rotor flux angle is estimated according to the following equation: cos⁢ ⁢θ^r=λryλry2+λry2
  • 4. The flux observer of claim 2, wherein said rotor speed signal is received from a speed observer.
  • 5. The flux observer of claim 4, wherein said speed observer receives current inputs and voltage inputs from a current observer and a voltage observer respectively.
  • 6. The flux observer of claim 2, wherein said rotor flux angle is determined according to the following equation: cos⁢ ⁢θ^r=λryλry2+λry2
  • 7. A discrete time flux observer for torque control of an induction motor, the induction motor having a rotor and a stator, comprising:a rotor speed input for receiving a signal indicative of the rotor speed; a stator current input for receiving a signal indicative of the stator current; a rotor flux estimating mechanism that estimates the rotor flux based on the assumption that the rotor speed is constant over a sampling interval; a rotor flux output for outputting said estimated rotor flux; a flux angle calculation for receiving said estimated rotor flux and determining a rotor flux angle; and a field oriented controller for receiving said determined rotor flux angle and determining a stator voltage; wherein the flux observer is used to control an induction motor that is incorporated into a vehicle electronic power assisted steering system.
  • 8. The flux observer of claim 7, wherein a discrete form of said rotor flux observer is expressed as follows: λrk+1=λk⁢e-σ⁢ ⁢T+η⁢ ⁢M⁢ ⁢isk⁢i-e-σ⁢ ⁢Tσ
  • 9. The flux observer of claim 8, wherein the rotor flux observer analytical digital form is determined as follows: λrx⁡(k+1)=(1-α)⁢λrxk-βλryk+η⁢ ⁢Mη2+ωr2⁡[(ηα+ωr⁢β)⁢isx-(ωr⁢α-ηβ)⁢isy]⁢⁢λry⁡(k+1)=(1-α)⁢λryk-βλrxk+η⁢ ⁢Mη2+ωr2⁡[(ηα+ωr⁢β)⁢isy-(ωr⁢α-ηβ)⁢isx].(8)
  • 10. The control scheme of claim 9, wherein the flux observer receives a pair of current inputs with the first current input being representative of the x component of the stator current and the second current input being representative of the y component of the stator current.
  • 11. The control scheme of claim 9, wherein said flux observer output or estimated rotor flux which is then used to determine a rotor flux angle.
  • 12. The control scheme of claim 11, wherein a discrete form of said rotor flux is expressed in accordance with the following equation: λrk+1=λk⁢ⅇ-σ⁢ ⁢T+η⁢ ⁢M⁢ ⁢isk⁢1-ⅇ-σ⁢ ⁢Tσ.
  • 13. The control scheme of claim 12, wherein an analytical digital form of said flux observer is determined as follows: λrx⁡(k+1)=(1-α)⁢λrxk-βλryk+η⁢ ⁢Mη2+ωr2⁡[(ηα+ωr⁢β)⁢isx-(ωr⁢α-ηβ)⁢isy]⁢⁢λry⁡(k+1)=(1-α)⁢λryk-βλrxk+η⁢ ⁢Mη2+ωr2⁡[(ηα+ωr⁢β)⁢isy-(ωr⁢α-ηβ)⁢isx].(8)
  • 14. A sensorless field orientation control scheme for controlling an induction motor, comprising:a flux observer for determining a rotor flux; a speed observer in communication with said flux observer to provide said flux observer with a signal indicative of the rotor speed; a voltage observer in communication with said speed observer to provide said speed observer with a signal representative of the filtered motor voltage; a current observer to provide said speed observer with a signal representative of the filtered motor current; and a field oriented controller in communication with said flux observer to control the torque of the motor; whereby said flux observer estimates said motor flux based on the assumption that the motor speed is constant over a sampling interval and wherein the inductive motor is incorporated into a vehicle electronic power assisted steering system.
  • 15. The control scheme of claim 11, wherein said rotor flux angle is determined according to the following equation: cos⁢ ⁢θ^r=λryλry2+λry2.
US Referenced Citations (16)
Number Name Date Kind
4707651 Schauder Nov 1987 A
4777422 Slicker et al. Oct 1988 A
5296793 Lang Mar 1994 A
5296794 Lang et al. Mar 1994 A
5365158 Tanaka et al. Nov 1994 A
5541488 Bansai et al. Jul 1996 A
5585709 Jansen Dec 1996 A
5708346 Schöb Jan 1998 A
5729113 Jansen et al. Mar 1998 A
5796235 Schrodl et al. Aug 1998 A
5796236 Royak Aug 1998 A
5973474 Yamamoto Oct 1999 A
6014006 Stuntz et al. Jan 2000 A
6069467 Jansen May 2000 A
6137258 Jansen Oct 2000 A
6163127 Patel et al. Dec 2000 A