CN105450124A - Method and apparatus for obtaining motor parameters - Google Patents

Method and apparatus for obtaining motor parameters Download PDF

Info

Publication number
CN105450124A
CN105450124A CN201510944054.8A CN201510944054A CN105450124A CN 105450124 A CN105450124 A CN 105450124A CN 201510944054 A CN201510944054 A CN 201510944054A CN 105450124 A CN105450124 A CN 105450124A
Authority
CN
China
Prior art keywords
prime
calculation
shaft current
result
take turns
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201510944054.8A
Other languages
Chinese (zh)
Other versions
CN105450124B (en
Inventor
陈跃
刘启武
唐婷婷
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Sichuan Changhong Electric Co Ltd
Original Assignee
Sichuan Changhong Electric Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Sichuan Changhong Electric Co Ltd filed Critical Sichuan Changhong Electric Co Ltd
Priority to CN201510944054.8A priority Critical patent/CN105450124B/en
Publication of CN105450124A publication Critical patent/CN105450124A/en
Application granted granted Critical
Publication of CN105450124B publication Critical patent/CN105450124B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P21/00Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
    • H02P21/14Estimation or adaptation of machine parameters, e.g. flux, current or voltage

Abstract

Embodiments of the invention provide a method and an apparatus for obtaining motor parameters, and aim to achieve a technical effect of obtaining actual parameters of a motor. The method comprises the steps of obtaining (n-1)th actual q axis current in the (n-1)th period of the motor, wherein n is a positive integer; obtaining (n-1)th putative q axis current output by an adjustable model corresponding to the motor, wherein the (n-1)th putative q axis current is the q axis current output by the adjustable model based on a (n-1)th-round calculation result; obtaining an nth-round calculation result according to the (n-1)th actual q axis current and the (n-1)th putative q axis current, wherein the nth-round calculation result comprises the nth motor parameters or the nth motor parameter correlative; judging whether the nth-round calculation result and the (n-1)th-round calculation result satisfy a preset relation or not; and if so, determining the nth motor parameters corresponding to the nth-round calculation result as the actual parameters of the motor.

Description

A kind of method and apparatus obtaining the parameter of electric machine
Technical field
The present invention relates to electronic technology field, particularly relate to a kind of method and apparatus obtaining the parameter of electric machine.
Background technology
At present, when controlling motor, the parameter of electric machine adopted is generally the specifications parameter that manufacturer provides, and such as motor inductances is with the discrete data of curent change, electric motor resistance and back electromotive force constant etc.But some parameters of electric machine are not thick-and-thin, but have relation with motor actual operating state, such as resistance can increase with the electric current of input motor and increase.Meanwhile, the environment of test specification parameter and the different of motor actual application environment, also may cause there is either large or small error between the specifications parameter that tests out and motor actual parameter.And the error of specifications parameter and actual parameter, cause Electric Machine Control inaccurate, operational efficiency is low and rotate the technical problems such as noise is large.
Summary of the invention
The embodiment of the present application provides a kind of method and apparatus obtaining the parameter of electric machine, for realizing the technique effect obtaining motor actual parameter.
First aspect, this application provides a kind of method obtaining the parameter of electric machine, comprising:
Obtain the (n-1)th actual q shaft current of motor in the (n-1)th cycle; N is positive integer;
Obtain that the adjustable model corresponding with described motor export (n-1)th estimates q shaft current; Described (n-1)th presumption q shaft current is that described adjustable model takes turns the q shaft current of result of calculation output based on (n-1)th;
According to described (n-1)th actual q shaft current and described (n-1)th presumption q shaft current, obtain n-th and take turns result of calculation; Described n-th takes turns result of calculation comprises n-th parameter of electric machine or the n-th parameter of electric machine correlative;
Judge that described n-th takes turns result of calculation and described (n-1)th and take turns result of calculation and whether meet preset relation;
When described n-th take turns result of calculation and described (n-1)th take turns result of calculation meet described preset relation time, determine that described n-th to take turns in result of calculation the actual parameter that corresponding described n-th parameter of electric machine is described motor.
Optionally, judging that described n-th takes turns result of calculation and described (n-1)th and take turns after whether result of calculation meet preset relation, also comprise:
When described n-th take turns result of calculation and described (n-1)th take turns result of calculation do not meet described preset relation time, take turns result of calculation by described n-th and input described adjustable model, with make described adjustable model based on described n-th take turns result of calculation the (n+1)th cycle export and described (n-1)th estimate q shaft current different n-th estimate q shaft current.
Optionally, the described parameter of electric machine comprises electric motor resistance, d axle inductance, q axle inductance and back electromotive force constant, according to described (n-1)th actual q shaft current and described (n-1)th presumption q shaft current, obtains n-th and takes turns result of calculation, comprising:
Obtain described n-th according to following formula and take turns result of calculation:
a n ′ = a 0 - Σ i = 0 n - 1 k 1 ( i q i - i q i ′ ) i q i ′ Δ T b n ′ = b 0 + Σ i = 0 n - 1 k 2 ( i q i - i q i ′ ) u q i ′ Δ T c n ′ = c 0 - Σ i = 0 n - 1 k 3 ( i q i - i q i ′ ) ω i Δ T d n ′ = d 0 - Σ i = 0 n - 1 k 4 ( i q i - i q i ′ ) ω i i d i ′ Δ T
Or
a n ′ = a 0 - Σ i = 0 n - 1 k 1 ( i q i - i q i ′ ) i q i ′ Δ T b n ′ = b 0 + Σ i = 0 n - 1 k 2 ( i q i - i q i ′ ) u q i ′ Δ T c n ′ = c 0 - Σ i = 0 n - 1 k 3 ( i q i - i q i ′ ) ω i Δ T d n ′ = d 0 - Σ i = 0 n - 1 k 5 ( i q i - i q i ′ ) ω i i d i ′ Δ T
Wherein, a n ′ = R n ′ L q n ′ , b n ′ = 1 L q n ′ , c n ′ = K e n ′ L q n ′ , d n ′ = L d n ′ L q n ′ , R n' be the n-th electric motor resistance, L dn' be the n-th d axle inductance, L qn' be the n-th q axle inductance, K en' be the n-th back electromotive force constant, r 0for the initial motor resistance of described adjustable model, L q0for the initial q axle inductance of described adjustable model, K e0for the initial back electromotive force constant of described adjustable model, L d0for the initial d axle inductance of described adjustable model, i qibe the i-th actual q shaft current, i dibe the i-th actual d shaft current, i qi' be the i-th presumption q shaft current, i di' be the i-th presumption d shaft current, Δ T is the cycle, ω ibe the i-th motor angular velocity, u qibe the i-th q shaft voltage, k 1, k 2, k 3, k 4and k 5be positive number.
Optionally, when described n-th result of calculation specifically comprises described n-th parameter of electric machine, judge that described n-th takes turns result of calculation and described (n-1)th and take turns result of calculation and whether meet preset relation, comprising:
Judge the (n-1)th electric motor resistance R n-1' and R n' ratio, described (n-1)th q shaft current L qn-1' and L qn' ratio, described (n-1)th back electromotive force constant K en-1' and K en' ratio and the (n-1)th d shaft current L dn-1' and L dn' ratio whether all in preset range;
Work as R n-1' and R n' ratio, L qn-1' and L qn' ratio, K en-1' and K en' ratio and L dn-1' and L dn' ratio all in described preset range, represent that described n-th takes turns result of calculation and described (n-1)th and take turns result of calculation and meet described preset relation.
Optionally, when described n-th take turns result specifically comprise described n-th parameter of electric machine correlative time, judge that described n-th takes turns result of calculation and described (n-1)th and take turns result of calculation and whether meet preset relation, comprising:
Judge a n-1' and a n' difference, b n-1' and b n' difference, c n-1' and c n' difference and d n-1' and d n' the absolute value of difference whether be all less than threshold value; Work as a n-1' and a n' difference, b n-1' and b n' difference, c n-1' and c n' difference and d n-1' and d n' the absolute value of difference when being all less than described threshold value, represent that described n-th takes turns result of calculation and described (n-1)th and take turns result of calculation and meet described preset relation.
Optionally, described adjustable model according to the following formula, is taken turns result of calculation based on described n-th and is exported described n-th presumption q shaft current in the (n+1)th cycle:
i q n ′ - i q n - 1 ′ Δ T = - R n ′ L q n ′ i q n ′ - L d n ′ L q n ′ ω n i d n ′ + u q n - K ′ e n ω n L ′ q n + k 1 ( i q n - i q n ′ )
Wherein, i qn-1' be described (n-1)th presumption q shaft current, i qn' be described n-th presumption q shaft current, i dn' be described n-th presumption d shaft current, i qnfor described n-th actual q shaft current, ω nbe the n-th motor angular velocity, u qnit is the n-th q shaft voltage.
On the other hand, this application provides a kind of device obtaining the parameter of electric machine, comprising:
Actual current obtains unit, for obtaining the (n-1)th actual q shaft current of motor in the (n-1)th cycle; N is positive integer;
Presumption electric current obtains unit, and export for obtaining the adjustable model corresponding with described motor (n-1)th estimates q shaft current; Described (n-1)th presumption q shaft current is that described adjustable model takes turns the q shaft current of result of calculation output based on (n-1)th;
Computing unit, for according to described (n-1)th actual q shaft current and described (n-1)th presumption q shaft current, obtains n-th and takes turns result of calculation; Described n-th takes turns result of calculation comprises n-th parameter of electric machine or the n-th parameter of electric machine correlative;
Judging unit, for judging that described n-th takes turns result of calculation and described (n-1)th and take turns result of calculation and whether meet preset relation;
Determining unit, for when described n-th take turns result of calculation and described (n-1)th take turns result of calculation meet described preset relation time, determine that described n-th to take turns in result of calculation the actual parameter that corresponding described n-th parameter of electric machine is described motor.
Optionally, described device also comprises:
Input unit, for judging that described n-th takes turns result of calculation and described (n-1)th and take turns after whether result of calculation meet preset relation, when described n-th take turns result of calculation and described (n-1)th take turns result of calculation do not meet described preset relation time, take turns result of calculation by described n-th and input described adjustable model, with make described adjustable model based on described n-th take turns result of calculation the (n+1)th cycle export and described (n-1)th estimate q shaft current different n-th estimate q shaft current.
Optionally, the described parameter of electric machine comprises electric motor resistance, d axle inductance, q axle inductance and back electromotive force constant, and described computing unit is used for obtaining described n-th according to following formula and takes turns result of calculation:
a n ′ = a 0 - Σ i = 0 n - 1 k 1 ( i q i - i q i ′ ) i q i ′ Δ T b n ′ = b 0 + Σ i = 0 n - 1 k 2 ( i q i - i q i ′ ) u q i ′ Δ T c n ′ = c 0 - Σ i = 0 n - 1 k 3 ( i q i - i q i ′ ) ω i Δ T d n ′ = d 0 - Σ i = 0 n - 1 k 4 ( i q i - i q i ′ ) ω i i d i ′ Δ T
Or
a n ′ = a 0 - Σ i = 0 n - 1 k 1 ( i q i - i q i ′ ) i q i ′ Δ T b n ′ = b 0 + Σ i = 0 n - 1 k 2 ( i q i - i q i ′ ) u q i ′ Δ T c n ′ = c 0 - Σ i = 0 n - 1 k 3 ( i q i - i q i ′ ) ω i Δ T d n ′ = d 0 - Σ i = 0 n - 1 k 5 ( i q i - i q i ′ ) ω i i d i ′ Δ T
Wherein, a n ′ = R n ′ L q n ′ , b n ′ = 1 L q n ′ , c n ′ = K e n ′ L q n ′ , d n ′ = L d n ′ L q n ′ , R n' be the n-th electric motor resistance, L dn' be the n-th d axle inductance, L qn' be the n-th q axle inductance, K en' be the n-th back electromotive force constant, r 0for the initial motor resistance of described adjustable model, L q0for the initial q axle inductance of described adjustable model, K e0for the initial back electromotive force constant of described adjustable model, L d0for the initial d axle inductance of described adjustable model, i qibe the i-th actual q shaft current, i dibe the i-th actual d shaft current, i qi' be the i-th presumption q shaft current, i di' be the i-th presumption d shaft current, Δ T is the cycle, ω ibe the i-th motor angular velocity, u qibe the i-th q shaft voltage, k 1, k 2, k 3, k 4and k 5be positive number.
Optionally, when described n-th result of calculation specifically comprises described n-th parameter of electric machine, described judging unit is for judging the (n-1)th electric motor resistance R n-1' and R n' ratio, described (n-1)th q shaft current L qn-1' and L qn' ratio, described (n-1)th back electromotive force constant K en-1' and K en' ratio and the (n-1)th d shaft current L dn-1' and L dn' ratio whether all in preset range;
Work as R n-1' and R n' ratio, L qn-1' and L qn' ratio, K en-1' and K en' ratio and L dn-1' and L dn' ratio all in described preset range, represent that described n-th takes turns result of calculation and described (n-1)th and take turns result of calculation and meet described preset relation.
Optionally, when described n-th take turns result specifically comprise described n-th parameter of electric machine correlative time, described judging unit is for judging a n-1' and a n' difference, b n-1' and b n' difference, c n-1' and c n' difference and d n-1' and d n' the absolute value of difference whether be all less than threshold value; Work as a n-1' and a n' difference, b n-1' and b n' difference, c n-1' and c n' difference and d n-1' and d n' the absolute value of difference when being all less than described threshold value, represent that described n-th takes turns result of calculation and described (n-1)th and take turns result of calculation and meet described preset relation.
Optionally, described adjustable model according to the following formula, is taken turns result of calculation based on described n-th and is exported described n-th presumption q shaft current in the (n+1)th cycle:
i q n ′ - i q n - 1 ′ Δ T = - R n ′ L q n ′ i q n ′ - L d n ′ L q n ′ ω n i d n ′ + u q n - K ′ e n ω n L ′ q n + k 1 ( i q n - i q n ′ )
Wherein, i qn-1' be described (n-1)th presumption q shaft current, i qn' be described n-th presumption q shaft current, i dn' be described n-th presumption d shaft current, i qnfor described n-th actual q shaft current, ω nbe the n-th motor angular velocity, u qnit is the n-th q shaft voltage.
Above-mentioned one or more technical scheme in the embodiment of the present application, at least has one or more technique effects following:
In the technical scheme of the embodiment of the present application, first, the (n-1)th actual q shaft current of motor in the (n-1)th cycle is obtained, and the (n-1)th presumption q shaft current that the adjustable model corresponding with motor exports.(n-1)th presumption q shaft current is that adjustable model takes turns the q shaft current of result of calculation output based on (n-1)th; N is positive integer.Then, according to the (n-1)th actual q shaft current and the (n-1)th presumption q shaft current, obtain n-th and take turns result of calculation, wherein, n-th takes turns result of calculation comprises n-th parameter of electric machine or the n-th parameter of electric machine correlative.Then, judge that n-th takes turns result of calculation and (n-1)th and take turns result of calculation and whether meet preset relation, when meeting preset relation, because adjustable model is corresponding with motor, so now the parameter of adjustable model is equal or close to the actual parameter of motor, so determine n-th, to take turns n-th parameter of electric machine corresponding to result of calculation be the actual parameter of motor.So, just obtained the actual parameter of motor by the technical scheme in the embodiment of the present application, so make to control motor more accurate.
Accompanying drawing explanation
Fig. 1 is the method flow diagram obtaining the parameter of electric machine in the embodiment of the present application;
Fig. 2 is the schematic diagram of the model framework obtaining the parameter of electric machine in the embodiment of the present application;
Fig. 3 is the apparatus structure schematic diagram obtaining the parameter of electric machine in the embodiment of the present application.
Embodiment
The embodiment of the present application provides a kind of method and apparatus obtaining the parameter of electric machine, for realizing the technique effect obtaining motor actual parameter.
In order to solve the problems of the technologies described above, the technical scheme general thought that the application provides is as follows:
In the technical scheme of the embodiment of the present application, first, the (n-1)th actual q shaft current of motor in the (n-1)th cycle is obtained, and the (n-1)th presumption q shaft current that the adjustable model corresponding with motor exports.(n-1)th presumption q shaft current is that adjustable model takes turns the q shaft current of result of calculation output based on (n-1)th; N is positive integer.Then, according to the (n-1)th actual q shaft current and the (n-1)th presumption q shaft current, obtain n-th and take turns result of calculation, wherein, n-th takes turns result of calculation comprises n-th parameter of electric machine or the n-th parameter of electric machine correlative.Then, judge that n-th takes turns result of calculation and (n-1)th and take turns result of calculation and whether meet preset relation, when meeting preset relation, because adjustable model is corresponding with motor, so now the parameter of adjustable model is equal or close to the actual parameter of motor, so determine n-th, to take turns n-th parameter of electric machine corresponding to result of calculation be the actual parameter of motor.So, just obtained the actual parameter of motor by the technical scheme in the embodiment of the present application, so make to control motor more accurate.
Below by accompanying drawing and specific embodiment, technical solution of the present invention is described in detail, the specific features being to be understood that in the embodiment of the present application and embodiment is the detailed description to technical scheme, instead of the restriction to technical scheme, when not conflicting, the technical characteristic in the embodiment of the present application and embodiment can combine mutually.
Term "and/or" herein, being only a kind of incidence relation describing affiliated partner, can there are three kinds of relations in expression, and such as, A and/or B, can represent: individualism A, exists A and B simultaneously, these three kinds of situations of individualism B.In addition, character "/" herein, general expression forward-backward correlation is to the relation liking a kind of "or".
The application's first aspect provides a kind of method obtaining the parameter of electric machine, please refer to Fig. 1, comprises the steps:
S101: obtain the (n-1)th actual q shaft current of motor in the (n-1)th cycle.
S102: obtain that the adjustable model corresponding with described motor export (n-1)th estimates q shaft current.
S103: according to described (n-1)th actual q shaft current and described (n-1)th presumption q shaft current, obtain n-th and take turns result of calculation.
S104: judge that described n-th takes turns result of calculation and described (n-1)th and take turns result of calculation and whether meet preset relation.
S105: when described n-th take turns result of calculation and described (n-1)th take turns result of calculation meet described preset relation time, determine that described n-th to take turns in result of calculation the actual parameter that corresponding described n-th parameter of electric machine is described motor.
In S101, obtain the (n-1)th actual q shaft current of motor in the (n-1)th cycle, detect specifically by current sensing means or control software design acquisition.In the embodiment of the present application, the cycle refers to the cycle of often taking turns calculating, and Cycle Length Δ T can be consistent with the cycle of pulse-modulated signal, such as 1.4 × 10 -4s is to 2.5 × 10 -4arbitrary value in s.N is positive integer, such as 1,10 or 112 etc., and the application does not limit.In addition, n is positive integer, such as 3,4 or 10 etc.
In S102, obtain the (n-1)th presumption q shaft current that the adjustable model corresponding with motor exports.Specifically, motor is in operation, and can obtain the realistic model of motor, i.e. Transmission function of motor.In the embodiment of the present application, the transfer function of adjustable model is consistent with the realistic model form of motor, and the function of Parameter adjustable.So adjustable model is corresponding with motor.The presumption q shaft current of adjustable model is according to the transfer function of adjustable model and input and the q shaft current obtained.For example, the realistic model of such as motor is
di d d t = - R L d i d + L q L d ωi q + u d L d di q d t = - R L q i q - L d L q ωi d + u q - K e ω L q , Formula (1)
Wherein, i dfor the actual d shaft current of motor, represent i dto t differentiate.I qfor the actual q shaft current of motor, represent i qto t differentiate.R is real electrical machinery resistance, L dfor the actual d axle inductance of motor, L qfor the actual q axle inductance of motor, ω is the actual angular speed of motor, u dfor the actual d shaft voltage of motor, u qfor the actual q shaft voltage of motor, K efor the back electromotive force constant of motor.Above-mentioned formula includes a large amount of parameter of electric machine in (1), and in the embodiment of the present application, only get second equation, namely formula (2) is as actual q shaft model.
di q d t = - R L q i q - L d L q ωi d + u q - K e ω L q Formula (2)
Order a = R L q , b = 1 L q , c = K e L q , d = L d L q , Formula (2) can be reduced to
di q d t = - ai q - dωi d + bu q - c ω . Formula (3)
So adjustable model is consistent with motor realistic model form
di d ′ d t = - R ′ L d ′ i d ′ + L q ′ L d ′ ωi q ′ + u d L d ′ di q ′ d t = - R ′ L q ′ i q ′ - L d ′ L q ′ ωi d ′ + u q - K e ′ ω L q ′ . Formula (4)
Wherein, i d' be the d shaft current of adjustable model, represent i d' to t differentiate.I q' be the q shaft current of adjustable model, represent i q' to t differentiate.R' is the electric motor resistance of adjustable model, L d' be the d axle inductance of adjustable model, L q' be the q axle inductance of adjustable model.ω is the angular speed of adjustable model, u dfor the d shaft voltage of adjustable model, u qfor the q shaft voltage of adjustable model, ω, u dand u qidentical with actual q shaft voltage with the actual angular speed of motor, actual d shaft voltage.K e' be the back electromotive force constant of adjustable model.
In the embodiment of the present application, the q shaft model of adjustable model is corresponding with actual q axle
di q ′ d t = - a ′ i q ′ - d ′ ωi d ′ + b ′ u q - c ′ ω , Formula (5)
Wherein, a ′ = R ′ L q ′ , b ′ = 1 L q ′ , c ′ = K e ′ L q ′ , d ′ = L d ′ L q ′ .
Further, in order to ensure that adjustable model is restrained, there is negative characteristic root, in the embodiment of the present application, introducing error minus feedback, to make the closed function of adjustable model being
di q ′ d t = - a ′ i q ′ - d ′ ωi d ′ + b ′ u q - c ′ ω + k 1 ( i q - i q ′ ) , Formula (6)
Wherein, k1 is positive number.
(n-1)th presumption q shaft current be adjustable model based on last round of result of calculation, namely (n-1)th take turns result of calculation and export q shaft current.Specifically, adjustable model takes turns result of calculation adjustment parameter based on (n-1)th, the electric motor resistance of such as adjustable model, and the d axle inductance of adjustable model and/or q axle inductance etc., then according to exporting the (n-1)th presumption q shaft current based on closed function and input.
In specific implementation process, for obtain motor actual parameter and carry out n-th take turns calculating time, S101 can perform prior to S102, and S102 also can perform prior to S101, can also perform S101 and S102 simultaneously, and the application does not do concrete restriction.
After obtaining the (n-1)th actual q shaft current and the (n-1)th presumption q shaft current, in S103, carry out n-th based on the (n-1)th actual q shaft current and the (n-1)th presumption q shaft current and take turns calculating, and then result of calculation is taken turns in acquisition n-th.As shown in Figure 2, in the embodiment of the present application, u dand u qit is input.Each takes turns calculating, from the actual q shaft current i that motor obtains q, and the presumption q shaft current i that adjustable model exports q' be all input in computing module and calculate.Computing module according to the algorithm process preset, and then exports epicycle result of calculation.
Specifically, n-th in the embodiment of the present application is taken turns result of calculation and is comprised n-th parameter of electric machine or parameter of electric machine correlative.In the embodiment of the present application, n-th parameter of electric machine is specially the n-th electric motor resistance R n', the n-th d axle inductance L dn', the n-th q axle inductance L qn' and the n-th back electromotive force constant K en'.N-th parameter of electric machine correlative is R n', L dn', L qn' and K en' dependent variable, be specially so n-th in the embodiment of the present application takes turns result of calculation at least can comprise R n', L dn', L qn' and K en', or at least comprise a n', b n', c n' and d n'.
Obtain n-th take turns result of calculation after, in S104, judge that n-th takes turns result of calculation and (n-1)th and take turns result of calculation and whether meet preset relation.In the embodiment of the present application, preset relation represents that the parameter of adjustable model is consistent with the parameter of motor realistic model or close.Further, when n-th take turns result of calculation and (n-1)th take turns result of calculation meet preset relation time, determine in S105 that n-th n-th parameter of electric machine of taking turns in result of calculation is the actual parameter of motor, or n-th n-th parameter of electric machine of taking turns the n-th motor correlative in result of calculation corresponding is the actual parameter of motor.
Specifically, because the functional form of adjustable model is consistent with the form of motor realistic model, and motor realistic model represents the virtual condition of motor, so, when n-th take turns result of calculation and (n-1)th take turns result of calculation meet preset relation time, represent that the parameter adjusted of adjustable model is to consistent with the parameter of motor realistic model or close, and then the adjustable model representing now is identical with motor realistic model or close to identical.Therefore, now n-th of adjustable model take turns n-th parameter of electric machine in result of calculation or the parameter of electric machine corresponding to the n-th parameter of electric machine correlative is just identical with the actual parameter of motor or close, therefore can as the actual parameter of motor.
So, seen from the above description, when n-th take turns result of calculation and (n-1)th take turns result of calculation meet preset relation time, adjustable model is equivalent to motor realistic model, n-th parameter of electric machine of adjustable model or the parameter of electric machine corresponding to the n-th parameter of electric machine correlative are just equivalent to the actual parameter of motor, so, using n-th parameter of electric machine or the parameter of electric machine corresponding to the n-th parameter of electric machine correlative as actual parameter, just obtain the actual parameter of motor.
In another embodiment of the application, after S104, can also comprise:
When described n-th take turns result of calculation and described (n-1)th take turns result of calculation do not meet described preset relation time, take turns result of calculation by described n-th and input described adjustable model, with make described adjustable model based on described n-th take turns result of calculation the (n+1)th cycle export and described (n-1)th estimate q shaft current different n-th estimate q shaft current.
Specifically, when judge in S104 n-th take turns result of calculation and (n-1)th take turns result of calculation do not meet preset relation time, represent that the parameter of now adjustable model is not also adjusted to consistent with motor or close consistent, so adjustable model still needs to continue adjustment parameter.So exported by computing module further n-th takes turns result of calculation input adjustable model, and then makes adjustable model adjust parameter, thus the (n+1)th cycle export and (n-1)th estimate q shaft current different n-th estimate q shaft current.
In the embodiment of the present application, if n-th takes turns result of calculation and (n-1)th take turns result of calculation when not meeting preset relation, result of calculation input adjustable model is taken turns by n-th, and then perform S101 to S104, until epicycle result of calculation and last round of result of calculation meet preset relation, determine that the actual parameter of motor just terminates.
Next, take turns result of calculation to how obtaining n-th describe in detail.
In the embodiment of the present application, according to following formula, obtain n-th and take turns result of calculation:
a n ′ = a 0 - Σ i = 0 n - 1 k 1 ( i q i - i q i ′ ) i q i ′ Δ T b n ′ = b 0 + Σ i = 0 n - 1 k 2 ( i q i - i q i ′ ) u q i ′ Δ T c n ′ = c 0 - Σ i = 0 n - 1 k 3 ( i q i - i q i ′ ) ω i Δ T d n ′ = d 0 - Σ i = 0 n - 1 k 4 ( i q i - i q i ′ ) ω i i d i ′ Δ T Formula (7)
Or
a n ′ = a 0 - Σ i = 0 n - 1 k 1 ( i q i - i q i ′ ) i q i ′ Δ T b n ′ = b 0 + Σ i = 0 n - 1 k 2 ( i q i - i q i ′ ) u q i ′ Δ T c n ′ = c 0 - Σ i = 0 n - 1 k 3 ( i q i - i q i ′ ) ω i Δ T d n ′ = d 0 - Σ i = 0 n - 1 k 5 ( i q i - i q i ′ ) ω i i d i ′ Δ T . Formula (8)
Specifically, r 0for the initial motor resistance of adjustable model, L q0for the initial q axle inductance of adjustable model, K e0for the initial back electromotive force constant of adjustable model, L d0for the initial d axle inductance of adjustable model.R 0, L q0, K e0and L d0arranged arbitrarily by those skilled in the art.Wherein, due to L q0for denominator, thus when arranging L q0specifically can be set to any non-zero real, the application does not do concrete restriction.I qibe the i-th actual q shaft current, i dibe the i-th actual d shaft current, i.e. the actual q shaft current of the i-th cycle motor output and d shaft current.I qi' be the i-th presumption q shaft current, i di' be the i-th presumption d shaft current, i.e. the q shaft current of the i-th cycle adjustable model output and d shaft current.Δ T is the cycle.ω ibe the i-th motor angular velocity, i.e. the actual angular speed of motor in the i-th cycle.U qibe the i-th q shaft voltage, i.e. the q shaft voltage of the i-th periodical input motor and adjustable model.K 1, k 2, k 3, k 4and k 5be positive number.In specific implementation process, k 1, k 2, k 3, k 4and k 5can identical also can not be identical, the application does not do concrete restriction.
At initial time, will as initial value, calculate a n', b n', c n' and d n'.To calculate a n' be example:
Work as n=1, namely carry out the 1st when taking turns calculating:
a 1 ′ = a 0 - Σ i = 0 0 k 1 ( i q 0 - i q 0 ′ ) i q 0 ′ Δ T = a 0 - k 1 ( i q 0 - i q 0 ′ ) i q 0 ′ Δ T .
Work as n=2, namely carry out the 2nd when taking turns calculating:
a 2 ′ = a 0 - Σ i = 0 1 k 1 ( i q 0 - i q 0 ′ ) i q 0 ′ Δ T = a 0 - k 1 ( i q 0 - i q 0 ′ ) i q 0 ′ Δ T - k 1 ( i q 1 - i q 1 ′ ) i q 1 ′ Δ T .
When n >=3, the like, just repeat no longer one by one here.So, the value from 1 of the n in the embodiment of the present application.
In specific implementation process, the application those of ordinary skill in the field can select formula (7) or formula (8) to obtain a n', b n', c n' and d n', the application does not do concrete restriction.
In the embodiment of the present application, n-th takes turns result of calculation specifically can comprise n-th parameter of electric machine, and n-th parameter of electric machine comprises R n', L dn', L qn' and K en'.Or n-th takes turns result of calculation also specifically can comprise the n-th parameter of electric machine correlative, and the n-th parameter of electric machine correlative at least comprises a n', b n', c n' and d n', pass through a n', b n', c n' and d n' acquisition R can be calculated further n', L dn', L qn' and K en'.
In the embodiment of the present application, motor realistic model is for specifically to get formula (3), and then adjustable model is specially formula (6).In specific implementation process, obtain a based on the (n-1)th actual q shaft current and the (n-1)th presumption q shaft current n', b n', c n' and d n' method have multiple, such as, based on liapunov function, Popov's inequality etc., the application does not do concrete restriction.In the embodiment of the present application, the derivation of above-mentioned formula (7) and formula (8) will be obtained in detail for Popov function.
First, formula (3) deducts formula (6),
e q · = ( - ai q + a ′ i ′ q + ai ′ q - ai ′ q ) - dωi d + d ′ ωi ′ d + ( b - b ′ ) u q + ( - c + c ′ ) ω - k 1 ( i q - i q ′ ) = ( - ai q + ai ′ q + a ′ i ′ q - ai ′ q ) - dωi d + d ′ ωi ′ d + ( b - b ′ ) u q + ( - c + c ′ ) ω - k 1 ( i q - i q ′ ) = - ae q - ( a - a ′ ) i q ′ - dωi d + d ′ ωi ′ d + ( b - b ′ ) u q + ( - c + c ′ ) ω - k 1 ( i q - i q ′ ) = - ae q - ( a - a ′ ) i q ′ - dωi d + dωi d ′ - dωi d ′ + d ′ ω i ′ d + ( b - b ′ ) u q + ( - c + c ′ ) ω - k 1 ( i q - i q ′ ) = - ae q - ( a - a ′ ) i q ′ - dωe d + ( - d ω + d ′ ω ) i d ′ + ( b - b ′ ) u q + ( - c + c ′ ) ω - k 1 ( i q - i q ′ ) = - ( a + k 1 ) e q - ( a - a ′ ) i q ′ - dωe d + ( - d ω + d ′ ω ) i d ′ + ( b - b ′ ) u q + ( - c + c ′ ) ω ,
Formula (9)
Wherein e d=i d-i d '', e q=i q-i q '', " " represents for differentiating operator.Due to when the parameter adjustment of adjustable model is to consistent with motor model or close to time consistent, a'=a, b'=b, c'=c, d'=d, i can be thought d=i' d, and then e d=i d-i d ''=0.So d ω e can be made d=0, and then formula (9) becomes
e q · = - ( a + k 1 ) e q - ( a - a ′ ) i q ′ + ( - d ω + d ′ ω ) i d ′ + ( b - b ′ ) u q + ( - c + c ′ ) ω . Formula (10)
Because liapunov function form is wherein p and Γ is symmetric positive definite matrix arbitrarily.So, in the embodiment of the present application, make A=[-a], negative feedback G=[-k1], φ t=[-a+a' ,-d+d', b-b' ,-c+c']=[a1, a2, a3, a4] and then: e · = ( A + G ) e + φ T s , Wherein
Structure liapunov function
V ( e , t ) · = 1 2 ( e T · p e + e T p e · + φ T · Γ φ + φ T Γ φ · ) . Formula (11)
Select P=[1], even P is unit matrix, then have:
e T · p e + e T p e · = [ ( A + G ) e + φ T s ] T p e + e T p ( ( A + G ) e + φ T s ) = e T ( A + G ) T p e + s T φ p e + e T p ( A + G ) e + e T pφ T s , Formula (12)
e T · p e + e T p e · = e T ( ( A + G ) T p + p ( A + G ) ) e + s T φ p e + e T pφ T s
And according to Lyapunov stability condition, liapunov function be made to stablize, need negative definite.In the embodiment of the present application, P=[1], then have: (A+G) tp+p (A+G)=-2 (a+k1), due to R > 0, L q> 0, so so as k1 > 0, (A+G) tp+p (A+G) negative definite.
Further, due to:
In addition, make in the embodiment of the present application Γ = g 1 , 0 , 0 , 0 0 , g 2 , 0 , 0 0 , 0 , g 3 , 0 0 , 0 , 0 , g 4 , G 1, g 2, g 3and g 4for positive number, therefore
φ T · Γ φ = [ a 1 · , a 2 · , a 3 · , a 4 · ] g 1 , 0 , 0 , 0 0 , g 2 , 0 , 0 0 , 0 , g 3 , 0 0 , 0 , 0 , g 4 [ a 1 , a 2 , a 3 , a 4 ] T = g 1 a 1 · a 1 + g 2 a · 2 a 2 + g 3 a 3 · a 3 + g 4 a · 4 a 4 . Formula (14)
φ T Γ φ · = g 1 a 1 · a 1 + g 2 a · 2 a 2 + g 3 a 3 · a 3 + g 4 a · 4 a 4
So, according to formula (12), formula (13) and formula (14), by formula (11) arrangement be
V ( e · , t ) = 1 2 ( e T · p e + e T p e · + φ T · Γ φ + φ T Γ φ · ) = 1 2 e T ( ( A + G ) T p + p ( A + G ) ) e + [ a 1 i ′ q e q + a 2 ωi d ′ e q + a 3 u q e q + a 4 ωe q ] + [ g 1 a 1 · a 1 + g 2 a · 2 a 2 + g 3 a 3 · a 3 + g 4 a 4 · a 4 ] = 1 2 ( e T ( A T p + p A ) e ) + a 1 ( g 1 a 1 · + i ′ q e q ) + a 2 ( g 2 a · 2 + ωi d ′ e q ) + a 3 ( g 3 a 3 · + u q e q ) + a 4 ( g 4 a 4 · + ωe q ) . Formula (15)
As can be seen here, when
g 1 a 1 · + i ′ q e q = 0 g 2 a · 2 + ωi d ′ e q = 0 ( g 3 a 3 · + u q e q ) = 0 g 4 a 4 · + ωe q = 0 When formula (16) is set up, negative definite.
Wherein, [-a+a' ,-d+d', b-b' ,-c+c']=[a1, a2, a3, a4], i.e. a1=-a+a', a2=-d+d', a3=b-b', the a4=-c+c' of formula (16).
And r, L q, L dand K ebeing all actual parameter, is definite value, so the derivative of a1, a2, a3 and a4 is 0.And then formula (16) can abbreviation be further
g 1 a ′ · + i ′ q e q = 0 g 2 d · ′ + ωi d ′ e q = 0 - g 3 b · ′ + u q e q = 0 g 4 c · ′ + ωe q = 0 . Formula (17)
By formula (17) transposition also integration, can obtain,
a ′ = a 0 - ∫ 0 t k 1 i q ′ e q d t b ′ = b 0 + ∫ 0 t k 2 u q e q d t c ′ = c 0 - ∫ 0 t 1 k 3 ωe q d t d ′ = d 0 - ∫ 0 t k 4 ωi d ′ e q d t , Formula (18)
Wherein, k 1 = 1 g 1 , k 2 = 1 g 3 , k 3 = 1 g 4 , k 4 = 1 g 2 .
In addition, due to when the parameter adjustment of adjustable model is to consistent with motor model or close to time consistent, a'=a, b'=b, c'=c, d'=d, i can be thought d=i' d.So can also be deformed into
d ′ = d 0 - ∫ 0 t k 5 ωi d e q d t . Formula (19)
Wherein k 5for positive number.
Therefore, except formula (18), formula (20) also can calculate a', b', c' and d'.
a ′ = a 0 - ∫ 0 t k 1 i q ′ e q d t b ′ = b 0 + ∫ 0 t k 2 u q e q d t c ′ = c 0 - ∫ 0 t 1 k 3 ωe q d t d ′ = d 0 - ∫ 0 t k 5 ωi d e q d t Formula (20)
Finally, for ease of computer disposal, by above-mentioned formula (18) discretization, obtain
a n ′ = a 0 - Σ i = 0 n - 1 k 1 ( i q i - i q i ′ ) i q i ′ Δ T b n ′ = b 0 + Σ i = 0 n - 1 k 2 ( i q i - i q i ′ ) u q i ′ Δ T c n ′ = c 0 - Σ i = 0 n - 1 k 3 ( i q i - i q i ′ ) ω i Δ T d n ′ = d 0 - Σ i = 0 n - 1 k 4 ( i q i - i q i ′ ) ω i i d i ′ Δ T Formula (7)
Or
a n ′ = a 0 - Σ i = 0 n - 1 k 1 ( i q i - i q i ′ ) i q i ′ Δ T b n ′ = b 0 + Σ i = 0 n - 1 k 2 ( i q i - i q i ′ ) u q i ′ Δ T c n ′ = c 0 - Σ i = 0 n - 1 k 3 ( i q i - i q i ′ ) ω i Δ T d n ′ = d 0 - Σ i = 0 n - 1 k 5 ( i q i - i q i ′ ) ω i i d i Δ T . Formula (8)
In the embodiment of the present application, n-th takes turns result of calculation comprises n-th parameter of electric machine, and wherein n-th parameter of electric machine is specially R n', L dn', L qn' and K en'.Or n-th takes turns result of calculation comprises the n-th parameter of electric machine correlative, and the n-th parameter of electric machine correlative at least comprises a n', b n', c n' and d n'.Due to n-th take turns result of calculation have two kinds may, so judge in S104 that n-th takes turns result of calculation and whether the (n-1)th result of calculation meets preset relation, also have two kinds of specific implementations.
The first:
When the n-th result of calculation specifically comprises n-th parameter of electric machine, and n-th parameter of electric machine at least comprises R n', L qn', K en' and L dn' time, S104 specifically comprises:
Judge the (n-1)th electric motor resistance R n-1' and R n' ratio, described (n-1)th q shaft current L qn-1' and L qn' ratio, described (n-1)th back electromotive force constant K en-1' and K en' ratio and the (n-1)th d shaft current L dn-1' and L dn' ratio whether all in preset range.
In the first implementation, for judging that n-th takes turns result of calculation and (n-1)th and take turns result of calculation and whether meet preset relation, (n-1)th parameter of electric machine corresponding to result of calculation is taken turns in (n-1)th of acquisition, comprises the (n-1)th electric motor resistance R n-1', the (n-1)th q axle inductance L qn-1', the (n-1)th back electromotive force constant K en-1' and the (n-1)th d axle inductance L dn-1'.Then R is obtained further n-1' and R n' ratio, L qn-1' and L qn' ratio, K en-1' and K en' ratio and L dn-1' and L dn' ratio be δ 1 = R n - 1 ′ R n ′ , δ 2 = L q n - 1 ′ L q n ′ , δ 3 = K e n - 1 ′ K e n ′ With δ 4 = L d n - 1 ′ L d n ′ .
In the first implementation, judge that n-th takes turns result of calculation and (n-1)th and take turns result of calculation and whether meet preset relation, be specially judgement δ 1 = R n - 1 ′ R n ′ , δ 2 = L q n - 1 ′ L q n ′ , δ 3 = K e n - 1 ′ K e n ′ With δ 4 = L d n - 1 ′ L d n ′ Whether all in preset range.In the embodiment of the present application, preset range be near 1 among a small circle, such as [0.995,1.005], or [0.998,1.002] etc.In specific implementation process, preset range is less, finally as n-th parameter of electric machine of real electrical machinery parameter more close to real electrical machinery parameter.The application those of ordinary skill in the field can be arranged according to actual, and the application does not do concrete restriction.
Specifically, when δ 1 = R n - 1 ′ R n ′ , δ 2 = L q n - 1 ′ L q n ′ , δ 3 = K e n - 1 ′ K e n ′ With δ 4 = L d n - 1 ′ L d n ′ Time all in preset range, represent δ 1 = R n - 1 ′ R n ′ , δ 2 = L q n - 1 ′ L q n ′ , δ 3 = K e n - 1 ′ K e n ′ With δ 4 = L d n - 1 ′ L d n ′ Equal 1 or close to 1.And δ 1 = R n - 1 ′ R n ′ , δ 2 = L q n - 1 ′ L q n ′ , δ 3 = K e n - 1 ′ K e n ′ With δ 4 = L d n - 1 ′ L d n ′ Equal 1 or close to 1, show i dn-1' equal or close to i dn-1, and i qn-1' equal or close to i qn-1.Because motor is consistent with the input of adjustable model, function model form is consistent, and output is now relevant or close relevant, so show that the parameter of now adjustable model is equal or close to motor realistic model, so can by R n', L qn', K en' and L dn' be considered as motor actual parameter R, L q, K eand L d.So, when δ 1 = R n - 1 ′ R n ′ , δ 2 = L q n - 1 ′ L q n ′ , δ 3 = K e n - 1 ′ K e n ′ With δ 4 = L d n - 1 ′ L d n ′ Time all in preset range, represent that n-th takes turns result of calculation and (n-1)th and take turns result of calculation and meet preset relation.
Further, in the first implementation, due to n-th, what take turns that result of calculation comprises is exactly n-th parameter of electric machine itself, so directly determine that in S105 n-th parameter of electric machine is the actual parameter of motor.
The second:
When n-th takes turns result of calculation and specifically comprise the n-th parameter of electric machine correlative, and the n-th motor correlative at least comprises a n', b n', c n' and d n' time, S104 specifically comprises:
Judge a n-1' and a n' difference, b n-1' and b n' difference, c n-1' and c n' difference and d n-1' and d n' the absolute value of difference whether be all less than threshold value.
When n-th take turns result of calculation specifically comprise the n-th parameter of electric machine correlative time, (n-1)th need be obtained and take turns the (n-1)th parameter of electric machine correlative a in result of calculation n-1', b n-1', c n-1' and d n-1'.Then a is obtained further n-1' and a n' difference, b n-1' and b n' difference, c n-1' and c n' difference and d n-1' and d n' the absolute value of difference.
In the second implementation, judge that n-th takes turns result of calculation and (n-1)th and take turns result of calculation and whether meet preset relation, be specially and judge a n-1' and a n' difference, b n-1' and b n' difference, c n-1' and c n' difference and d n-1' and d n' the absolute value of difference whether be all less than threshold value.In the embodiment of the present application, threshold value be 0 or close to 0 number, such as 0,0.05, or 0.1 etc.In specific implementation process, threshold value is less, more close to 0, finally as n-th parameter of electric machine of real electrical machinery parameter more close to real electrical machinery parameter.The application those of ordinary skill in the field can be arranged according to actual, and the application does not do concrete restriction.
Specifically, a n-1' and a n' difference, b n-1' and b n' difference, c n-1' and c n' difference and d n-1' and d n' the absolute value of difference when being all less than threshold value, show i dn-1' equal or close to i dn-1, i qn-1' equal or close to i qn-1, and then, show that the parameter of now adjustable model is equal or close to motor realistic model, so can by R n', L qn', K en' and L dn' be considered as motor actual parameter R, L q, K eand L d.So, work as a n-1' and a n' difference, b n-1' and b n' difference, c n-1' and c n' difference and d n-1' and d n' the absolute value of difference when being all less than threshold value, represent that n-th takes turns result of calculation and (n-1)th and take turns result of calculation and meet preset relation.
Further, what take turns that result of calculation specifically comprises is the n-th parameter of electric machine correlative due to n-th so in the second implementation, n-th parameter of electric machine corresponding to the n-th parameter of electric machine correlative determined in S105 is just with L d n ′ = d n ′ b n ′ .
Further, when n-th take turns result of calculation and (n-1)th take turns result of calculation do not meet pre-conditioned time, take turns result of calculation input adjustable model by n-th, according to formula (6), adjustable model exports the n-th presumption q shaft current according to formula (21) in the (n+1)th cycle.
di q n ′ d t = - a n ′ i q n ′ - d n ′ ω n i d n ′ + b n ′ u q n - c n ′ ω n + k 1 ( i q n - i q n ′ ) Formula (21)
Wherein, i dn' be the n-th presumption d shaft current, i qn' be the n-th presumption q shaft current, ω nbe the actual angle of motor in the n-th cycle, u qnbe the q shaft voltage of the n-th periodical input adjustable model and motor, i qnit is the n-th actual q shaft current.
Further, in specific implementation process, if utilize computer calculate, then computer is specially to utilize and is calculated by the formula after discretization.So the formula (21) after discretization is
i q n ′ - i q n - 1 ′ Δ T = - R n ′ L q n ′ i q n ′ - L d n ′ L q n ′ ω n i d n ′ + u q n - K ′ e n ω n L ′ q n + k 1 ( i q n - i q n ′ ) . Formula (22)
Based on the inventive concept same with the method obtaining the parameter of electric machine in previous embodiment, the application's second aspect also provides a kind of device obtaining the parameter of electric machine, as shown in Figure 3, comprising:
Actual current obtains unit 301, for obtaining the (n-1)th actual q shaft current of motor in the (n-1)th cycle; N is positive integer;
Presumption electric current obtains unit 302, and export for obtaining the adjustable model corresponding with described motor (n-1)th estimates q shaft current; Described (n-1)th presumption q shaft current is that described adjustable model takes turns the q shaft current of result of calculation output based on (n-1)th;
Computing unit 303, for according to described (n-1)th actual q shaft current and described (n-1)th presumption q shaft current, obtains n-th and takes turns result of calculation; Described n-th takes turns result of calculation comprises n-th parameter of electric machine or the n-th parameter of electric machine correlative;
Judging unit 304, for judging that described n-th takes turns result of calculation and described (n-1)th and take turns result of calculation and whether meet preset relation;
Determining unit 305, for when described n-th take turns result of calculation and described (n-1)th take turns result of calculation meet described preset relation time, determine that described n-th to take turns in result of calculation the actual parameter that corresponding described n-th parameter of electric machine is described motor.Further, the device in the embodiment of the present application can further include:
Input unit, for judging that described n-th takes turns result of calculation and described (n-1)th and take turns after whether result of calculation meet preset relation, when described n-th take turns result of calculation and described (n-1)th take turns result of calculation do not meet described preset relation time, take turns result of calculation by described n-th and input described adjustable model, with make described adjustable model based on described n-th take turns result of calculation the (n+1)th cycle export and described (n-1)th estimate q shaft current different n-th estimate q shaft current.
Specifically, the parameter of electric machine in the embodiment of the present application comprises electric motor resistance, d axle inductance, q axle inductance and back electromotive force constant, and computing unit 303 takes turns result of calculation for obtaining described n-th according to following formula:
a n ′ = a 0 - Σ i = 0 n - 1 k 1 ( i q i - i q i ′ ) i q i ′ Δ T b n ′ = b 0 + Σ i = 0 n - 1 k 2 ( i q i - i q i ′ ) u q i ′ Δ T c n ′ = c 0 - Σ i = 0 n - 1 k 3 ( i q i - i q i ′ ) ω i Δ T d n ′ = d 0 - Σ i = 0 n - 1 k 4 ( i q i - i q i ′ ) ω i i d i ′ Δ T
Or
a n ′ = a 0 - Σ i = 0 n - 1 k 1 ( i q i - i q i ′ ) i q i ′ Δ T b n ′ = b 0 + Σ i = 0 n - 1 k 2 ( i q i - i q i ′ ) u q i ′ Δ T c n ′ = c 0 - Σ i = 0 n - 1 k 3 ( i q i - i q i ′ ) ω i Δ T d n ′ = d 0 - Σ i = 0 n - 1 k 5 ( i q i - i q i ′ ) ω i i d i Δ T
Wherein, a n ′ = R n ′ L q n ′ , b n ′ = 1 L q n ′ , c n ′ = K e n ′ L q n ′ , d n ′ = L d n ′ L q n ′ , R n' be the n-th electric motor resistance, L dn' be the n-th d axle inductance, L qn' be the n-th q axle inductance, K en' be the n-th back electromotive force constant, a 0 = R 0 L q 0 , b 0 = 1 L q 0 , c 0 = K e 0 L q 0 , d 0 = L d 0 L q 0 , R 0for the initial motor resistance of described adjustable model, L q0for the initial q axle inductance of described adjustable model, K e0for the initial back electromotive force constant of described adjustable model, L d0for the initial d axle inductance of described adjustable model, i qibe the i-th actual q shaft current, i dibe the i-th actual d shaft current, i qi' be the i-th presumption q shaft current, i di' be the i-th presumption d shaft current, Δ T is the cycle, ω ibe the i-th motor angular velocity, u qibe the i-th q shaft voltage, k 1, k 2, k 3, k 4and k 5be positive number.
Optionally, when the n-th result of calculation specifically comprises n-th parameter of electric machine, judging unit 304 is for judging the (n-1)th electric motor resistance R n-1' and R n' ratio, the (n-1)th q shaft current L qn-1' and L qn' ratio, the (n-1)th back electromotive force constant K en-1' and K en' ratio and the (n-1)th d shaft current L dn-1' and L dn' ratio whether all in preset range;
Work as R n-1' and R n' ratio, L qn-1' and L qn' ratio, K en-1' and K en' ratio and L dn-1' and L dn' ratio all in preset range, represent that n-th takes turns result of calculation and (n-1)th and take turns result of calculation and meet preset relation.
Or, when n-th take turns result specifically comprise the n-th parameter of electric machine correlative time, judging unit 304 is for judging a n-1' and a n' difference, b n-1' and b n' difference, c n-1' and c n' difference and d n-1' and d n' the absolute value of difference whether be all less than threshold value; Work as a n-1' and a n' difference, b n-1' and b n' difference, c n-1' and c n' difference and d n-1' and d n' the absolute value of difference when being all less than threshold value, represent that n-th takes turns result of calculation and (n-1)th and take turns result of calculation and meet preset relation.
Further, adjustable model according to the following formula, is taken turns result of calculation based on n-th and is exported the n-th presumption q shaft current in the (n+1)th cycle:
i q n ′ - i q n - 1 ′ Δ T = - R n ′ L q n ′ i q n ′ - L d n ′ L q n ′ ω n i d n ′ + u q n - K ′ e n ω n L ′ q n + k 1 ( i q n - i q n ′ )
Wherein, i qn-1' be described (n-1)th presumption q shaft current, i qn' be described n-th presumption q shaft current, i dn' be described n-th presumption d shaft current, i qnfor described n-th actual q shaft current, ω nbe the n-th motor angular velocity, u qnit is the n-th q shaft voltage.
The various variation pattern of the method for the parameter of the acquisition motor in earlier figures 1-Fig. 2 embodiment and instantiation are equally applicable to the device of the parameter of the acquisition motor of the present embodiment, by the detailed description of the method for the aforementioned parameter to acquisition motor, those skilled in the art clearly can know in the present embodiment the implementation method of the device of the parameter obtaining motor, so succinct in order to specification, be not described in detail in this.
Above-mentioned one or more technical scheme in the embodiment of the present application, at least has one or more technique effects following:
In the technical scheme of the embodiment of the present application, first, the (n-1)th actual q shaft current of motor in the (n-1)th cycle is obtained, and the (n-1)th presumption q shaft current that the adjustable model corresponding with motor exports.(n-1)th presumption q shaft current is that adjustable model takes turns the q shaft current of result of calculation output based on (n-1)th; N is positive integer.Then, according to the (n-1)th actual q shaft current and the (n-1)th presumption q shaft current, obtain n-th and take turns result of calculation, wherein, n-th takes turns result of calculation comprises n-th parameter of electric machine or the n-th parameter of electric machine correlative.Then, judge that n-th takes turns result of calculation and (n-1)th and take turns result of calculation and whether meet preset relation, when meeting preset relation, because adjustable model is corresponding with motor, so now the parameter of adjustable model is equal or close to the actual parameter of motor, so determine n-th, to take turns n-th parameter of electric machine corresponding to result of calculation be the actual parameter of motor.So, just obtained the actual parameter of motor by the technical scheme in the embodiment of the present application, so make to control motor more accurate.
Obviously, those skilled in the art can carry out various change and modification to the present invention and not depart from the spirit and scope of the present invention.Like this, if these amendments of the present invention and modification belong within the scope of the claims in the present invention and equivalent technologies thereof, then the present invention is also intended to comprise these change and modification.

Claims (12)

1. obtain a method for the parameter of electric machine, it is characterized in that, comprising:
Obtain the (n-1)th actual q shaft current of motor in the (n-1)th cycle; N is positive integer;
Obtain that the adjustable model corresponding with described motor export (n-1)th estimates q shaft current; Described (n-1)th presumption q shaft current is that described adjustable model takes turns the q shaft current of result of calculation output based on (n-1)th;
According to described (n-1)th actual q shaft current and described (n-1)th presumption q shaft current, obtain n-th and take turns result of calculation; Described n-th takes turns result of calculation comprises n-th parameter of electric machine or the n-th parameter of electric machine correlative;
Judge that described n-th takes turns result of calculation and described (n-1)th and take turns result of calculation and whether meet preset relation;
When described n-th take turns result of calculation and described (n-1)th take turns result of calculation meet described preset relation time, determine that described n-th to take turns in result of calculation the actual parameter that corresponding described n-th parameter of electric machine is described motor.
2. the method for claim 1, is characterized in that, judging that described n-th takes turns result of calculation and described (n-1)th and take turns after whether result of calculation meet preset relation, also comprises:
When described n-th take turns result of calculation and described (n-1)th take turns result of calculation do not meet described preset relation time, take turns result of calculation by described n-th and input described adjustable model, with make described adjustable model based on described n-th take turns result of calculation the (n+1)th cycle export and described (n-1)th estimate q shaft current different n-th estimate q shaft current.
3. method as claimed in claim 2, it is characterized in that, the described parameter of electric machine comprises electric motor resistance, d axle inductance, q axle inductance and back electromotive force constant, according to described (n-1)th actual q shaft current and described (n-1)th presumption q shaft current, obtain n-th and take turns result of calculation, comprising:
Obtain described n-th according to following formula and take turns result of calculation:
a n ′ = a 0 - Σ i = 0 n - 1 k 1 ( i q i - i q i ′ ) i q i ′ Δ T b n ′ = b 0 + Σ i = 0 n - 1 k 2 ( i q i - i q i ′ ) u q i ′ Δ T c n ′ = c 0 - Σ i = 0 n - 1 k 3 ( i q i - i q i ′ ) ω i Δ T d n ′ = d 0 - Σ i = 0 n - 1 k 4 ( i q i - i q i ′ ) ω i i d i ′ Δ T
Or
a n ′ = a 0 - Σ i = 0 n - 1 k 1 ( i q i - i q i ′ ) i q i ′ Δ T b n ′ = b 0 + Σ i = 0 n - 1 k 2 ( i q i - i q i ′ ) u q i ′ Δ T c n ′ = c 0 - Σ i = 0 n - 1 k 3 ( i q i - i q i ′ ) ω i Δ T d n ′ = d 0 - Σ i = 0 n - 1 k 5 ( i q i - i q i ′ ) ω i i d i Δ T
Wherein, a n ′ = R n ′ L q n ′ , b n ′ = 1 L q n ′ , c n ′ = K e n ′ L q n ′ , d n ′ = L d n ′ L q n ′ , R n' be the n-th electric motor resistance, L dn' be the n-th d axle inductance, L qn' be the n-th q axle inductance, K en' be the n-th back electromotive force constant, r 0for the initial motor resistance of described adjustable model, L q0for the initial q axle inductance of described adjustable model, K e0for the initial back electromotive force constant of described adjustable model, L d0for the initial d axle inductance of described adjustable model, i qibe the i-th actual q shaft current, i dibe the i-th actual d shaft current, i qi' be the i-th presumption q shaft current, i di' be the i-th presumption d shaft current, Δ T is the cycle, ω ibe the i-th motor angular velocity, u qibe the i-th q shaft voltage, k 1, k 2, k 3, k 4and k 5be positive number.
4. method as claimed in claim 3, is characterized in that, when described n-th result of calculation specifically comprises described n-th parameter of electric machine, judges that described n-th takes turns result of calculation and described (n-1)th and take turns result of calculation and whether meet preset relation, comprising:
Judge the (n-1)th electric motor resistance R n-1' and R n' ratio, described (n-1)th q shaft current L qn-1' and L qn' ratio, described (n-1)th back electromotive force constant K en-1' and K en' ratio and the (n-1)th d shaft current L dn-1' and L dn' ratio whether all in preset range;
Work as R n-1' and R n' ratio, L qn-1' and L qn' ratio, K en-1' and K en' ratio and L dn-1' and L dn' ratio all in described preset range, represent that described n-th takes turns result of calculation and described (n-1)th and take turns result of calculation and meet described preset relation.
5. method as claimed in claim 3, is characterized in that, when described n-th take turns result specifically comprise described n-th parameter of electric machine correlative time, judge that described n-th takes turns result of calculation and described (n-1)th and take turns result of calculation and whether meet preset relation, comprising:
Judge a n-1' and a n' difference, b n-1' and b n' difference, c n-1' and c n' difference and d n-1' and d n' the absolute value of difference whether be all less than threshold value; Work as a n-1' and a n' difference, b n-1' and b n' difference, c n-1' and c n' difference and d n-1' and d n' the absolute value of difference when being all less than described threshold value, represent that described n-th takes turns result of calculation and described (n-1)th and take turns result of calculation and meet described preset relation.
6. the method as described in claim 4 or 5, is characterized in that, described adjustable model according to the following formula, is taken turns result of calculation based on described n-th and exported described n-th presumption q shaft current in the (n+1)th cycle:
i q n ′ - i q n - 1 ′ Δ T = - R n ′ L q n ′ i q n ′ - L d n ′ L q n ′ ω n i d n ′ + u q n - K ′ e n ω n L ′ q n + k 1 ( i q n - i q n ′ )
Wherein, i qn-1' be described (n-1)th presumption q shaft current, i qn' be described n-th presumption q shaft current, i dn' be described n-th presumption d shaft current, i qnfor described n-th actual q shaft current, ω nbe the n-th motor angular velocity, u qnit is the n-th q shaft voltage.
7. obtain a device for the parameter of electric machine, it is characterized in that, comprising:
Actual current obtains unit, for obtaining the (n-1)th actual q shaft current of motor in the (n-1)th cycle; N is positive integer;
Presumption electric current obtains unit, and export for obtaining the adjustable model corresponding with described motor (n-1)th estimates q shaft current; Described (n-1)th presumption q shaft current is that described adjustable model takes turns the q shaft current of result of calculation output based on (n-1)th;
Computing unit, for according to described (n-1)th actual q shaft current and described (n-1)th presumption q shaft current, obtains n-th and takes turns result of calculation; Described n-th takes turns result of calculation comprises n-th parameter of electric machine or the n-th parameter of electric machine correlative;
Judging unit, for judging that described n-th takes turns result of calculation and described (n-1)th and take turns result of calculation and whether meet preset relation;
Determining unit, for when described n-th take turns result of calculation and described (n-1)th take turns result of calculation meet described preset relation time, determine that described n-th to take turns in result of calculation the actual parameter that corresponding described n-th parameter of electric machine is described motor.
8. device as claimed in claim 7, it is characterized in that, described device also comprises:
Input unit, for judging that described n-th takes turns result of calculation and described (n-1)th and take turns after whether result of calculation meet preset relation, when described n-th take turns result of calculation and described (n-1)th take turns result of calculation do not meet described preset relation time, take turns result of calculation by described n-th and input described adjustable model, with make described adjustable model based on described n-th take turns result of calculation the (n+1)th cycle export and described (n-1)th estimate q shaft current different n-th estimate q shaft current.
9. device as claimed in claim 8, it is characterized in that, the described parameter of electric machine comprises electric motor resistance, d axle inductance, q axle inductance and back electromotive force constant, and described computing unit is used for obtaining described n-th according to following formula and takes turns result of calculation:
a n ′ = a 0 - Σ i = 0 n - 1 k 1 ( i q i - i q i ′ ) i q i ′ Δ T b n ′ = b 0 + Σ i = 0 n - 1 k 2 ( i q i - i q i ′ ) u q i ′ Δ T c n ′ = c 0 - Σ i = 0 n - 1 k 3 ( i q i - i q i ′ ) ω i Δ T d n ′ = d 0 - Σ i = 0 n - 1 k 4 ( i q i - i q i ′ ) ω i i d i ′ Δ T
Or
a n ′ = a 0 - Σ i = 0 n - 1 k 1 ( i q i - i q i ′ ) i q i ′ Δ T b n ′ = b 0 + Σ i = 0 n - 1 k 2 ( i q i - i q i ′ ) u q i ′ Δ T c n ′ = c 0 - Σ i = 0 n - 1 k 3 ( i q i - i q i ′ ) ω i Δ T d n ′ = d 0 - Σ i = 0 n - 1 k 5 ( i q i - i q i ′ ) ω i i d i Δ T
Wherein, a n ′ = R n ′ L q n ′ , b n ′ = 1 L q n ′ , c n ′ = K e n ′ L q n ′ , d n ′ = L d n ′ L q n ′ , R n' be the n-th electric motor resistance, L dn' be the n-th d axle inductance, L qn' be the n-th q axle inductance, K en' be the n-th back electromotive force constant, r 0for the initial motor resistance of described adjustable model, L q0for the initial q axle inductance of described adjustable model, K e0for the initial back electromotive force constant of described adjustable model, L d0for the initial d axle inductance of described adjustable model, i qibe the i-th actual q shaft current, i dibe the i-th actual d shaft current, i qi' be the i-th presumption q shaft current, i di' be the i-th presumption d shaft current, Δ T is the cycle, ω ibe the i-th motor angular velocity, u qibe the i-th q shaft voltage, k 1, k 2, k 3, k 4and k 5be positive number.
10. device as claimed in claim 9, it is characterized in that, when described n-th result of calculation specifically comprises described n-th parameter of electric machine, described judging unit is for judging the (n-1)th electric motor resistance R n-1' and R n' ratio, described (n-1)th q shaft current L qn-1' and L qn' ratio, described (n-1)th back electromotive force constant K en-1' and K en' ratio and the (n-1)th d shaft current L dn-1' and L dn' ratio whether all in preset range;
Work as R n-1' and R n' ratio, L qn-1' and L qn' ratio, K en-1' and K en' ratio and L dn-1' and L dn' ratio all in described preset range, represent that described n-th takes turns result of calculation and described (n-1)th and take turns result of calculation and meet described preset relation.
11. devices as claimed in claim 9, is characterized in that, when described n-th take turns result specifically comprise described n-th parameter of electric machine correlative time, described judging unit is for judging a n-1' and a n' difference, b n-1' and b n' difference, c n-1' and c n' difference and d n-1' and d n' the absolute value of difference whether be all less than threshold value; Work as a n-1' and a n' difference, b n-1' and b n' difference, c n-1' and c n' difference and d n-1' and d n' the absolute value of difference when being all less than described threshold value, represent that described n-th takes turns result of calculation and described (n-1)th and take turns result of calculation and meet described preset relation.
12. devices as described in claim 10 or 11, is characterized in that, described adjustable model according to the following formula, is taken turns result of calculation based on described n-th and exported described n-th presumption q shaft current in the (n+1)th cycle:
i q n ′ - i q n - 1 ′ Δ T = - R n ′ L q n ′ i q n ′ - L d n ′ L q n ′ ω n i d n ′ + u q n - K ′ e n ω n L ′ q n + k 1 ( i q n - i q n ′ )
Wherein, i qn-1' be described (n-1)th presumption q shaft current, i qn' be described n-th presumption q shaft current, i dn' be described n-th presumption d shaft current, i qnfor described n-th actual q shaft current, ω nbe the n-th motor angular velocity, u qnit is the n-th q shaft voltage.
CN201510944054.8A 2015-12-16 2015-12-16 A kind of method and apparatus for obtaining the parameter of electric machine Active CN105450124B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510944054.8A CN105450124B (en) 2015-12-16 2015-12-16 A kind of method and apparatus for obtaining the parameter of electric machine

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510944054.8A CN105450124B (en) 2015-12-16 2015-12-16 A kind of method and apparatus for obtaining the parameter of electric machine

Publications (2)

Publication Number Publication Date
CN105450124A true CN105450124A (en) 2016-03-30
CN105450124B CN105450124B (en) 2018-05-01

Family

ID=55559997

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510944054.8A Active CN105450124B (en) 2015-12-16 2015-12-16 A kind of method and apparatus for obtaining the parameter of electric machine

Country Status (1)

Country Link
CN (1) CN105450124B (en)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2013215064A (en) * 2012-04-04 2013-10-17 Nsk Ltd Motor control device and electric power steering device mounted with the same
CN103501150A (en) * 2013-10-12 2014-01-08 上海联孚新能源科技有限公司 Embedded permanent magnet synchronous motor parameter identification device and method
US20140167659A1 (en) * 2011-09-28 2014-06-19 Denso Corporation Motor control device and motor control method
CN104539211A (en) * 2014-12-18 2015-04-22 珠海格力节能环保制冷技术研究中心有限公司 Motor parameter identification method and device

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140167659A1 (en) * 2011-09-28 2014-06-19 Denso Corporation Motor control device and motor control method
JP2013215064A (en) * 2012-04-04 2013-10-17 Nsk Ltd Motor control device and electric power steering device mounted with the same
CN103501150A (en) * 2013-10-12 2014-01-08 上海联孚新能源科技有限公司 Embedded permanent magnet synchronous motor parameter identification device and method
CN104539211A (en) * 2014-12-18 2015-04-22 珠海格力节能环保制冷技术研究中心有限公司 Motor parameter identification method and device

Also Published As

Publication number Publication date
CN105450124B (en) 2018-05-01

Similar Documents

Publication Publication Date Title
CN106909074B (en) Control device
CN107121932A (en) Motor servo system error symbol integrates Robust Adaptive Control method
CN103684183B (en) Rotational Speed of Asynchronous Motor method of estimation
CN105424276A (en) Method and device for obtaining rotational inertia of motor
CN105353789A (en) Continuous vibration signal time history replication control method
CN106094701A (en) Control device and control method, control system, control program and record medium
CN106169897A (en) A kind of motor speed real-time anti-saturation PID control method and device
Oliva et al. An engineering method for the power flow assessment in servo-actuated automated machinery: Mechatronic modeling and experimental evaluation
CN103942380B (en) Graphical control system design and simulation tool
CN109426143A (en) Load torque evaluation method, system, Mechatronic control system, method and motor
CN105450124A (en) Method and apparatus for obtaining motor parameters
CN104238376A (en) Low-temperature and low-air-pressure environment generator set running simulating system based on data
Burkus et al. A validation procedure to identify joint friction, reductor self-locking and gear backlash parameters
CN106998161A (en) The method for identification of rotational inertia and device of a kind of servo-drive system
CN104991445A (en) Globally-stable motor servo system self-adaptive output feedback robustness control method
CN105553371A (en) Method and device for obtaining motor parameters
CN106526299A (en) Current detection method for power converter based on non-smooth observation technology
CN105424277A (en) Method and device for obtaining rotational inertia of motor
CN105429544A (en) Method and apparatus for obtaining motor parameters
CN106208844B (en) A kind of motor servo system output feedback robust control method of Existence of Global Stable
CN104251784A (en) Reliability accelerated testing method of combined stress of integrated mechanical and electrical product
Baser et al. Transparency improvement in haptic devices with a torque compensator using motor current
CN105553370A (en) Motor control method and electronic equipment
Dahale et al. Design and Implementation of Controller using MPC Toolbox
KR102570839B1 (en) Method and apparatus of generating pulse for estimating a position of rotator in motor for phil simulator, and motor control system using the same

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant