Patents
Search within
Search within the title, abstract, claims, or full patent document: You can restrict your search to a specific field using field names.
Use TI= to search in the title, AB= for the abstract, CL= for the claims, or TAC= for all three. For example, TI=(safety belt).
Search by Cooperative Patent Classifications (CPCs): These are commonly used to represent ideas in place of keywords, and can also be entered in a search term box. If you're searching forseat belts, you could also search for B60R22/00 to retrieve documents that mention safety belts or body harnesses. CPC=B60R22 will match documents with exactly this CPC, CPC=B60R22/low matches documents with this CPC or a child classification of this CPC.
Learn More
Title
Abstract
Claims
Full Document
Find patents
Keywords and boolean syntax (USPTO or EPO format): seat belt searches these two words, or their plurals and close synonyms. "seat belt" searches this exact phrase, in order. -seat -belt searches for documents not containing either word.
For searches using boolean logic, the default operator is AND with left associativity. Note: this means safety OR seat belt is searched as (safety OR seat) AND belt. Each word automatically includes plurals and close synonyms. Adjacent words that are implicitly ANDed together, such as (safety belt), are treated as a phrase when generating synonyms.
Learn More
Search by
Chemistry searches match terms (trade names, IUPAC names, etc. extracted from the entire document, and processed from .MOL files.)
Substructure (use SSS=) and similarity (use ~) searches are limited to one per search at the top-level AND condition. Exact searches can be used multiple times throughout the search query.
Searching by SMILES or InChi key requires no special syntax. To search by SMARTS, use SMARTS=.
To search for multiple molecules, select "Batch" in the "Type" menu. Enter multiple molecules separated by whitespace or by comma.
Learn More
Patent numbers
Search specific patents by importing a CSV or list of patent publication or application numbers.
Induction motor control system with speed and flux estimation
US6014006A
United States
- Inventor
Ross Maxwell Stuntz Vladimir Aleksandrovich Bondarko Alexander Timofeevich Zaremba - Current Assignee
- Ford Global Technologies LLC
Worldwide applications
Application US09/339,485 events
1999-06-24
2000-01-11
Application granted
2000-01-11
2000-06-23
2019-06-24
Anticipated expiration
Status
Expired - Fee Related
Description
translated from
1. Field of the Invention
The present invention relates to control systems for induction motors.
2. Description of Related Art
In recent years, there has been substantial research on advanced control systems for controlling induction motors without rotational transducers to achieve high dynamic performance of induction motor drives. A major element of these control systems has been the estimation of the rotor position or the position of the rotor flux without using a rotor position sensor. The techniques used for these purposes can be divided into two groups: a) tracking machine spatial saliencies and b) using a variety of control methods for speed and flux estimation.
The spatial saliencies tracking methods require injection of the separate high frequency carrier signal and work best at low speed. They are based on the accurate analysis of machine physical structure irregularities and have intrinsic limitations originating from the requirements of separation of fundamental and carrier frequencies and detection of low amplitude spatial information signal.
The control methods from the second group rely on different control techniques. One of the approaches is based on augmentation of the plant to include dynamics of unknown parameters and use of state dependent Riccati equation filter methodology. A particular control method employs two complementary speed and flux observers or estimators. The main high gain speed and flux estimator works in the motoring and generating region and requires the knowledge of motor speed sign. The auxiliary estimator operates in generating and braking modes but can be justified if motor speed varies sufficiently slowly. However, the switching between the two estimators can cause divergence when the induction motor is operated at low speed under large load. The high gain estimator is utilized in a reference frame related with a flux estimation vector. Rigorous closed loop analysis taking into account rotor resistance variation, and error in estimation of rotor flux and speed provides conditions under which asymptotic torque tracking is achieved.
An object of the present invention is to provide method and apparatus for controlling an induction motor using pseudo (filtered) electrical current and voltage signals and observable time derivatives thereof as inputs to a system of linear algebraic equations representing motor dynamics.
The present invention provides method and apparatus for controlling an induction motor wherein a system of linear algebraic equations is developed to define a speed observer or estimator based on motor dynamics. The speed observer equations have rotor speed as the only unknown variable and have pseudo (filtered) electrical current and voltage signals and observable derivatives thereof as inputs. Rotor speed is estimated as a rational function of the measured and observable variable stator current and voltage inputs. Rotor flux then is estimated directly from rotor speed using an induction motor equation. Rotor flux angle is determined from rotor flux for use in controlling torque of the induction motor.
The sole FIG. 1 is a schematic view of an induction motor and control system pursuant to an embodiment of the invention.
Referring to the FIGURE, a block diagram of a motor control system in accordance with an embodiment of the present invention is illustrated. In particular, an induction motor 10 is shown driven by a pulse width modulating (PWM) voltage source inverter 12, a field oriented controller 14, synchronous-to-stationary frame transformation block 16, two to three phase transformation block 17, a stator current observer 22, stator voltage observer 24, three to two phase transformation block 26, a speed observer or estimator 30, and a flux observer 32 providing estimated flux values to rotor flux calculator 34.
Polyphase current rotating in the same direction as, but faster than, the rotor 10a of the induction motor 10 will create a torque in a direction so as to aid in rotor rotation. In this operational mode, net electrical power flow is into the motor electrical terminals, and the motor 10 is considered to be in a "motoring" mode. Conversely, polyphase currents rotating slower than the motor create a torque in a direction opposing the rotor rotation. Net electrical power flow is thus out of the motor electrical terminals, and the motor 10 is considered to be in a "generating" mode.
The voltage observer 24 receives commanded stator voltage signals ux and uy calculated in the stationary frame and generates pseudo voltage signals u0 and u1 (as described below) and sends them to speed observer or estimator 30.
The speed observer 30 estimates the rotational electrical speed of rotor 10a by processing a system of algebraic equations to be described below using the pseudo stator current signals is0, is1, is2 and pseudo stator voltage signals u0, u1 as inputs to the equations.
The estimated rotor speed, ωr is sent to the flux observer 32 which estimates the rotor flux, λr, directly from rotor speed using an induction motor equation described below. The rotor flux is estimated directly from the rotor speed and stator current signals. The estimated rotor flux is used by the rotor flux angle calculator 34 to calculate the rotor flux angle Θr. The rotor flux angle is used with the torque command, Te *, and rotor flux command, λr *, and measured stator currents isx, isy to control induction motor 10 via conventional techniques such as indirect field orientation methods incorporating slip frequency calculator and current regulator. Current observer 22, voltage observer 24, speed observer 30, flux observer 32, and flux angle calculator 34 are embodied in a microprocessing system of a motor control unit.
The estimated rotor flux angle, Θr, is sent to transform block 16 for purposes of transforming voltage control signal from a synchronous frame to stator stationary frame.
A dynamic model of a three phase induction machine 10, in a conventional stationary stator reference frame (x, y) is as follows: ##EQU1##
Here the variables Θ and ω are the angular position and mechanical speed of the rotor 10a, λr, us are rotor flux, stator current and voltage in the stator frame of reference (x, y), respectively. λr T is a transposed vector for λr.
The parameters Rr, Rs are rotor and stator resistances, M is a mutual inductance, Lr=M+L1r, Ls =M+L1s are rotor and stator inductances, σ=1-M2 /Ls Lr is the leakage parameter, m is the moment of inertia of the motor, np is the number of pole pairs, and TL is an external load torque.
For vector control of induction motor 10, it is required to measure or estimate rotor mechanical speed ω. Then the necessary stator frequency command ωs * can be obtained by summing the commanded slip frequency ωs * with the measured rotor speed ω:
ω.sub.e *=n.sub.p ω+ω.sub.s * (5)
where ##EQU2## and Te *, λr * are reference (commanded) values of torque and rotor flux.
Integrating electrical frequency signal ωe * (5) provides information about orientation of rotor flux which allows separation of the problems of flux control and reference torque tracking. An object of the invention is to design the speed observer 30 for speed estimation and the flux observer 32 for flux estimation using only information about stator current signals isx, isy and stator voltage signals ux, uy.
The speed observer design assumes that all motor parameters have constant known values. To further simplify the problem, the rotor speed ω is assumed to be constant, but unknown. The second assumption can be justified by the fact that actual rotor speed varies slowly as compared to rotor flux and stator current. Rotor speed changes can be neglected in design of the speed observer 30.
The speed observer 30 may be briefly described as follows:
Due to assumption ω≈constant, the last two equations of the dynamic model (3), (4) are sufficient to describe the dynamics of induction motor 10.
These equations (3), (4) are linear and allow elimination of unobservable flux and obtainment of a differential equation having stator current and voltage signals only.
The same equations (3), (4) are valid for filtered versions of current and voltage, whose derivatives are observable.
Since rotor velocity is the only unknown variable of the equations, it may be calculated as a rational function of filtered current, voltage and their observable derivatives.
Having estimated rotor speed, we can evaluate flux directly from equation (3).
The speed observer 30 embodies equations (3), (4) rewritten as follows: ##EQU3## where is =(isx, isy), u=(ux, uy), ωr is rotor electrical speed and ωr =np ω, u is stator voltage, ##EQU4##
Differentiating equation (7) and eliminating i, by using equation (8) we get: ##EQU5## Equation (10) can be represented in the form:
a+Jb·ω.sub.r =0 (11)
where a and b are 2-vectors as follows: ##EQU6## Thus, we get an overdetermined linear algebraic system (11) with respect to ωr. Minimizing mean square residual |a+Jbωr |2 gives an estimation of the rotor speed: ##EQU7## where bt represents the transpose of the vector b. Equation (14) contains derivatives of current and voltage signals and can not be used directly to determine rotor speed. To resolve this problem, low pass filtered versions of stator current and stator voltage are introduced by current observer 22 and voltage observer 24 as pseudo current is0 and voltage u0 as follows ##EQU8## where G>0 and G is the filter gain. Equation (15) provides pseudo (filtered) current and voltage is0 and u0 from measured stator current is and stator voltage u.
Substituting (15) into (11) gives ##EQU9## This means that:
a.sub.G +Jb.sub.G ·=ε (19)
where ε is a small number and ε→0 exponentially as t (time)→∞; or, roughly speaking, we can use for speed estimation the following equation:
a.sub.G +Jb.sub.G ·ω.sub.r =0. (20)
Therefore, we get: ##EQU10## Observable first and second order derivatives of pseudo (filtered) current is0 and first order derivative of pseudo (filtered) voltage u0 are:
i.sub.s1 =di.sub.s0 |dt, i.sub.s2 =d.sup.2 i.sub.s0 |dt.sup.2, u.sub.1 =du.sub.0 |dt (24)
can be determined as a solution of equations: ##EQU11## The last equality in (25) follows from: ##EQU12##
Equations (21)-(25) define the speed observer 30. The speed is obtained as a solution (21) of the system of algebraic equation with pseudo (filtered) current and voltage signals is0, is1, is2 and u0, u1 as inputs.
The estimated rotor speed, ωr, is sent to the flux observer 32 which estimates the rotor flux, λr, using equation (3) rewritten as: ##EQU13## is an exponentially stable matrix. The flux estimation signal λr is bounded for any bounded input is.
The estimated rotor flux is used by the rotor flux angle calculator 34 to calculate the rotor flux angle, Θr, using the equation: ##EQU14## where λrx and λry are rotor flux with respect to stator reference frame (x, y).
The rotor flux angle is used with the torque command, Te *, and rotor flux command, λr *, and measured current signals isx, isy to control induction motor 10 via conventional techniques such as indirect field orientation methods incorporating slip frequency calculator and current regulator.
Speed Observer Optimized To Noise
To improve speed observer robustness with respect to noise in stator current measurements, a speed estimation is obtained as a solution of a continuous set of equations (20) over given time interval. The usual strategy is to put maximum weight on recent measurements and gradually discount previous measurements. Minimizing integral mean error: ##EQU15## where t is the time moment, s is a dummy variable of integration, and discount factor e-gs is selected as an exponential function, and leading to the following formulae: ##EQU16## where bG ' is defined by equations (23)-(25), z,900 and Φ are the filtered versions of ζ, φ, respectively, and g is the filter gain.
By comparing (30)-(34) with (21), it can be concluded that the speed observer 30 obtained from minimization of (29) contains additional filtering (32), (33) of the denominator and numerator of expression (21).
Analysis of singular point
Considering the situation when b=0 and consequently a=0 in equation (11) provides a singular case at which induction motor speed can not be determined from system (11).
The motor dynamics for b=0 are governed by the equation ##EQU17##
Substituting equation (35) into equation (7) we get:
ηMi.sub.s -Aλ.sub.r =0 (36)
From (36) it follows that at the singular point, the current and flux are related by algebraic equality. Transforming (36) to the (d, q) stator reference frame related with rotor flux, we get two equivalent equations:
Mi.sub.d =λ.sub.d (37)
ηMi.sub.q +ω.sub.r λ.sub.d =0 (38)
where ##EQU18## where id and iq are stator currents in the respective d and q stator references frames and λd is flux in the (d, q) stator reference frame and λr t is the transpose of λr.
Comparing the equation (38) with (5), we find in the singular case that the electrical excitation frequency ωe equals zero:
ω.sub.e =0. (39)
This corresponds to the special situation when the electrical field is not rotating and machine slip, s: ##EQU19## is infinite.
In the vicinity of the observer singular point, an additional speed observer (not shown) may be used if desired or needed and based on dynamic equations: ##EQU20## where
b=-i.sub.s1 +(βηM-γ)i.sub.s +δu. (42)
From equations (41) and (8), the following formulae are generated: ##EQU21## which define flux and speed changes in the area where ωe ≅0.
While the invention is described above in terms of specific embodiments thereof, it is not intended to be limited thereto but rather only to the extent set forth hereafter in the following claims.
Claims (10)
Hide Dependent
translated from
Claims (10)
Hide Dependent
translated from
1. A method of controlling torque of a multi-phase induction motor, comprising:
(a) estimating rotor speed as a function of filtered stator current and voltage signals and first and second order derivatives of the filtered stator current signal and first order derivative of the filtered stator voltage signal,
(b) estimating rotor flux from said estimated rotor speed,
(c) estimating rotor flux angle from said estimated rotor flux, and
(d) using the estimated rotor flux angle to control motor torque.
2. The method of claim 1 wherein said current signal is filtered by: ##EQU22## where is is stator current, is0 is a filtered pseudo stator current, and G is greater than 0.
3. The method of claim 2 wherein said voltage signal is filtered by: ##EQU23## where u is stator voltage, u0 is a filtered pseudo stator voltage, and G is greater than 0.
4. The method of claim 3 wherein rotor speed is estimated as a solution of equations below: ##EQU24## .
5. A method of estimating rotor speed of a multi-phase induction motor, comprising estimating rotor speed as a function of filtered stator current and voltage signals and first and second order derivatives of the filtered stator current signal and first order derivative of the filtered stator voltage signal, wherein said current signal is filtered by: ##EQU25## where is is stator current, is0 is a filtered pseudo stator current, and G is greater than 0 and wherein said voltage signal is filtered by: ##EQU26## where u is stator voltage, u0 is a filtered pseudo stator voltage, and G is greater than 0.
6. A rotor speed estimator, comprising:
(a) a current observer that filters measured stator current signals to provide pseudo current signals thereof,
(b) a voltage observer that filters measured stator voltage signals to provide pseudo voltage signals thereof,
(c) a rotor speed observer that estimates rotor speed as a function of the pseudo stator current and voltage signals and first and second order derivatives of the pseudo current signal and first order derivative of the pseudo voltage signal.
7. The estimator of claim 6 including a second rotor speed observer for a singular point of motor dynamics.
8. The estimator of claim 6 wherein said current observer filters said measured stator current signal by: ##EQU27## where is is stator current, is0 is a filtered pseudo stator current, and G is greater than 0.
9. The method of claim 6 wherein said voltage observer filters said measured voltage signal by: ##EQU28## where u is stator voltage, u0 is a filtered pseudo stator voltage, and G is greater than 0.
10. The method of claim 6 wherein rotor speed is estimated as a solution of: ##EQU29## .