CN106602952A - Flux linkage full-rank identification method for permanent magnet of PMSM - Google Patents

Flux linkage full-rank identification method for permanent magnet of PMSM Download PDF

Info

Publication number
CN106602952A
CN106602952A CN201610512256.XA CN201610512256A CN106602952A CN 106602952 A CN106602952 A CN 106602952A CN 201610512256 A CN201610512256 A CN 201610512256A CN 106602952 A CN106602952 A CN 106602952A
Authority
CN
China
Prior art keywords
pmsm
permanent magnet
identification
magnet flux
flux linkage
Prior art date
Application number
CN201610512256.XA
Other languages
Chinese (zh)
Other versions
CN106602952B (en
Inventor
陈涛
郭壮志
任鹏飞
薛鹏
程辉
肖海红
Original Assignee
河南工程学院
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 河南工程学院 filed Critical 河南工程学院
Priority to CN201610512256.XA priority Critical patent/CN106602952B/en
Publication of CN106602952A publication Critical patent/CN106602952A/en
Application granted granted Critical
Publication of CN106602952B publication Critical patent/CN106602952B/en

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

The invention discloses a flux linkage full-rank identification method for a permanent magnet of a PMSM. According to the method, an unscented-Kalman-filter-algorithm-based parameter sensitivity analysis for the flux linkage identification precision of a permanent magnet of a PMSM is carried out; a permanent magnet flux linkage full-rank identification method of different rotating speed zones of a PMSM drive system is determined. With the unscented-Kalman-filter-algorithm-based permanent-magnet flux linkage full-rank identification method of different rotating speed zones of a PMSM drive system, influences on the permanent magnet flux linkage identification precision by the system measurement noises, PMSM parameter changing, and a degenerate-rank phenomenon of an identification model can be eliminated; linear error occurrence and complicated Jacobian matrix calculation during the state identification process of the traditional extended Kalman filter algorithm can be avoided; and the algorithm implementation difficulty can be reduce. Meanwhile, the number of to-be-identified parameters can be reduced at a high-speed operating zone in a PMSM drive system; the operand of the identification algorithm can be reduced; and reasonable consideration of the identification precision and the identification speed can be realized.

Description

A kind of PMSM permanent magnet flux linkages full rank discrimination method

Technical field

The present invention relates to the technical field of permagnetic synchronous motor (PMSM), and in particular to a kind of PMSM permanent magnet flux linkages full rank Discrimination method.

Background technology

The technical advantage such as PMSM has simple structure, fault rate low and operational efficiency is high, gradually in industrial servo-drive, new The fields such as energy automobile are applied widely.However, for many applications such as electric automobile, PMSM power densities Height, radiating condition is poor, and operating condition is complicated, under jointly controlling in torque capacity electric current ratio and weak magnetic more, there is stronger electricity Pivot reacts, in addition the factor such as natural aging, permanent magnet uniformly demagnetization or local demagnetization failure easily occurs, causes motor output to turn Square is reduced and torque pulsation, directly affects the direct torque precision and operational reliability of power drive system.

Meanwhile, in the control algolithm such as MTPA controls, weak magnetic control in PMSM, it is both needed to accurate permanent magnet magnetic chain information. Therefore, in order to realize PMSM drive system permanent magnets health status monitor and high performance controller design, it is necessary to it is accurate to obtain Permanent magnet flux linkage.

In order to accurately obtain PMSM permanent magnet flux linkages, French scholar Henwood N and Shandong University Wang Song are respectively adopted dragon Burger observer realizes the observation of PMSM permanent magnet flux linkages with method of least square, yet with observed result to the quick of measurement noise Perception, limits application of the program in actual industrial system.In order to solve the problems, such as noise jamming, Tsing-Hua University Xiao Xi professors It is assumed that on the premise of other PMSM parameter constants, being estimated online to permanent magnet flux linkage using expanded Kalman filtration algorithm Meter, but affected by magnetic circuit saturation and operation temperature rise, stator resistance RsAnd d-q axle stator inductance LdAnd LqCan drive with PMSM There is different degrees of change in the change of dynamic system conditions, affects the estimated accuracy of permanent magnet flux linkage.

Need to be state variable by the PMSM parameter processings of change to ensure PMSM permanent magnet flux linkage identification precisions, realize Permanent magnet flux linkage identification under Parameters variation constraint.For this purpose, peace group's great waves is established recognize PMSM stator resistance R simultaneouslys, d-q axles Stator inductance LdAnd LqAnd permanent magnet flux linkage ψfAdaptive model;The electricity then based on PMSM q axles such as Cortes-Romero J A Pressure equation, all electromagnetic parameters of PMSM are realized using algebraically identification algorithm in the case of participating in without the need for parameter of electric machine initial value On-line identification.However, the order for realizing the PMSM state equations that permanent magnet flux linkage is recognized is 2, while two parameters can only be realized Full rank is recognized, stator resistance Rs, stator inductance LdAnd LqAnd permanent magnet flux linkage ψfWhile identification exist identification equation owe order ask Topic, the uniqueness of identification result lacks theory support, easily causes identification result to be absorbed in local optimum and even dissipates.For this purpose, northern Capital Aero-Space university professor Wang Lina proposes a kind of based on model reference adaptive calculation for face mounted permagnetic synchronous motor The step identification method of method, the method is first with d shaft voltage equation estimation PMSM armature inductance Ls, recycle the L for obtainingsIt is real Existing ψfWith RsSynchronous full rank identification.Due to adopting i face mounted permagnetic synchronous motor mored=0 control, in order to realize ψf、RsIt is same Step full rank identification, need to inject amplitude and the rational d axles current perturbation of frequency, and not account in PMSM drive system runnings LsChange is to ψfThe impact of identification precision.For this purpose, France scholar Underwood S J are based on IPM synchronous motor Rs、 Ld、Lq、ψfThe different time scales of four parameters, are divided into fast varying parameter and slow varying parameter, and using two different time chis The least-squares algorithm of degree realizes the on-line identification of aforementioned four parameter, although the deficient order that the method can solve identification model is asked Topic, but in order to ensure algorithmic statement, still need to inject amplitude and the rational current perturbation of frequency to guarantee slow time scale most in d axles The identification precision of young waiter in a wineshop or an inn's multiplication algorithm, and identification result is easily by measurement influence of noise.Artificial intelligence with evolution algorithm as representative by In with stronger Nonlinear Processing ability, there is certain application in the identification of PMSM permanent magnet flux linkages, but how to reduce it Full rank recognizes amount of calculation, is but still key technical problem urgently to be resolved hurrily.

In consideration of it, the PMSM permanent magnet flux linkages that systematic survey noise, PMSM Parameters variations, identification equation are owed under order constraint are high Precision, online, full rank identification have become the monitoring of PMSM permanent magnets health status and high performance control field pass urgently to be resolved hurrily Key problem.

The content of the invention

In order to solve above-mentioned technical problem, the present invention provides a kind of PMSM permanent magnet flux linkages full rank discrimination method, based on double Unscented kalman filtering algorithm carries out subregion joint full rank identification to PMSM permanent magnet flux linkages, realizes measurement noise, PMSM parameters PMSM permanent magnet flux linkages high accuracy, online, full rank identification under order constraint are owed in change, identification model, are to realize PMSM drivetrains The high performance control of system and the on-line monitoring of permanent magnet health status provide foundation.

In order to achieve the above object, the technical scheme is that:A kind of PMSM permanent magnet flux linkages full rank discrimination method, adopts Stator voltage u in the d-q shaftings of collection PMSMd、uq, stator current id、iqAnd PMSM drive system rotor electrical angular velocity omegase, Setting up in d-q shaftings is used for the PMSM state equations of permanent magnet flux linkage identification, realizes based on Unscented kalman filtering algorithm The parameters sensitivity analysis of PMSM permanent magnet flux linkage identification precisions, based on analysis result PMSM drive system different rotating speeds area is determined The permanent magnet flux linkage full rank discrimination method based on double Unscented kalman filtering algorithms, realize measurement noise, permagnetic synchronous motor Parameters variation and identification model owe permanent magnet flux linkage high accuracy, online, full rank identification under order constraint;Its step is as follows:

Stator voltage u in step one, the d-q shaftings of collection PMSMdAnd uq, stator current idAnd iqAnd PMSM drivetrains System rotor electrical angular velocity omegae

Step 2, the PMSM state equations set up in d-q shaftings, based on Unscented kalman filtering algorithm PMSM permanent magnetism is realized Body magnetic linkage is recognized, under analysis PMSM drive system difference operating conditions, stator resistance Rs, d-q axle stator inductance LdAnd LqChange is right The impact of permanent magnet flux linkage identification precision, realizes the parameter based on the permanent magnet flux linkage identification precision of Unscented kalman filtering algorithm Sensitivity analyses;

Step 3, according to step 2 result, in PMSM drive system low cruises area, by d-q axle stator inductance Ld、LqConnection Close identification and stator resistance Rs, permanent magnet flux linkage ψfJoint identification is combined, updated each other, under identification equation full rank state, is disappeared Except stator resistance Rs, d-q axle stator inductance LdAnd LqImpact of the change to permanent magnet flux linkage estimated accuracy;In PMSM drive systems Middle and high fast Operational Zone, using permanent magnet flux linkage ψfIdentification and d-q axle stator inductance Ld、LqJoint identification combines, in identification side Under journey full rank state, d-q axle stator inductance L are eliminatedd、LqImpact of the change to permanent magnet flux linkage identification precision.

Stator voltage u in the d-q shaftings of the PMSMdAnd uq, stator current idAnd iqAcquisition methods be:Sampling PMSM Stator line voltage uab、ubc, three-phase current ia、ib、ic, and obtained by coordinate transform, transformation matrix of coordinates is respectively:

In formula, θ is flux linkage position of the rotor angle.

Stator voltage u in the d-q shaftings of the PMSMdAnd uq, stator current idAnd iqAcquisition methods be:Directly adopt The d-q shaft voltage command values that PMSM driving system controllers are calculatedWithReplace d-q axle stator voltages udAnd uq, d-q axles Current instruction valueReplace d-q axle stator current idAnd iq

The acquisition methods of the PMSM drive systems rotor electrical angular velocity are obtained by incremental optical-electricity encoder, step Suddenly it is:

(1) in t1And t2Umber of pulse N that optical rotary encoder sends on the d-q axles of neighbouring sample instance sample PMSM1、 N2, sampling instant t1And t2Difference be sampling period T;

(2) according to rotor angular rate ωeWith optical rotary encoder impulse sampling value N1、N2And between sampling period T Relation calculate rotor angular rate ωe, its expression formula is:

In formula, M is the optical rotary encoder umber of pulse of a week, and p is permagnetic synchronous motor number of pole-pairs.

It is used to realize that the method for PMSM state equations that permanent magnet flux linkage is recognized is in the d-q shaftings:

In formula, udAnd uqD-q axle stator voltages, i are represented respectivelydAnd iqD-q axle stator currents, L are represented respectivelydAnd LqRespectively Represent d-q axle stator inductances;RsRepresent stator resistance, ψfRepresent permanent magnet flux linkage, ωeRepresent rotor electric angle frequency;

In PMSM state equations, state vector is:X=[id iq ψf]T, input quantity is:U=[ud/Ld uq/Lq]T, output Measure and be:Y=[id iq]T

It is described to realize that the method that PMSM permanent magnet flux linkages are recognized is based on Unscented kalman filtering algorithm:Permanent magnet will be realized The PMSM state equations of magnetic linkage identification are described as the measurement equation table of the general type of nonlinear system, state equation and discretization It is shown as:

In formula, x (t) be system mode vector, y (tk) it is output, f () represents systematic state transfer equation, h () Represent system measuring equation, σ (t) is the process noise for considering model uncertainty and measuring uncertainty, μ (tk) it is to consider mould Uncertain and measuring uncertainty the measurement noise of type, the variance matrix of σ (t) is Q (t), μ (tk) variance matrix be R T (), u (t) is definitiveness input vector, B represents control matrix;

Above-mentioned PMSM state equations, state vector, input vector, output vector are substituted into into Unscented kalman filtering algorithm In, realize based on the state vector recursion of Unscented kalman filtering algorithm, you can realize the on-line identification of PMSM permanent magnet flux linkages, Step is:

(1) state vector initialization

Process noise covariance matrix Q and measurement noise covariance matrix R initial values are set according to priori, shape is initialized State vectorInit state varivance matrix

(2) Sigma points are calculated

Pass through within each sampling periodCalculate Sigma Point (k=1,2, ∞) and, obtain the Sigma dot matrixs of n × (2n+1);Wherein:χk-1For Sigma dot matrixs,Pk-1It is expressed as the prediction average and covariance of a moment state vector, n is state vector dimension, λ=α2(n+ B) it is scale parameter, α is dispersion level of the Sigma points near state variable average, and b is scale coefficient;

(3) time renewal

Realize that Sigma points are transmitted by the system state equation of discretization:I.e.According to Transmission result obtains the prediction average and covariance of state vector:

Wherein,The average weight of state vector is represented,The covariance weight of state vector is represented,WithQuantitative relation be:

(4) measurement updaue

According to system measurement data, by formula It is capable of achieving the one-step prediction of state vector Update and updated with Kalman gain, obtain the optimal estimation of state vector and variance matrix;Wherein, H represents calculation matrix, KkTable Show Kalman gain;

K=k+1 is made, the iteration output of state vector is realized in repeat step (2)-(4).

It is described to realize based on the parameters sensitivity analysis of the PMSM permanent magnet flux linkage identification precisions of Unscented kalman algorithm Method is:It is determined that the PMSM stator resistance R affected by system conditions and operation temperature rises, d axle inductances Ld, q axle inductances Lq's Excursion, and the permanent magnet flux linkage identification under different system operating mode based on Unscented kalman filtering algorithm is analyzed within this range Error.

In the PMSM drive systems low cruise region, full rank identification equation is:

In the middle and high fast operation area of the PMSM drive systems, full rank identification equation is:

Beneficial effect of the present invention:Compared with prior art, the present invention can owe under order constraint, to realize PMSM in identification model The full rank identification of permanent magnet flux linkage, eliminates measurement noise, PMSM Parameters variations and identification model and owes order to permanent magnet flux linkage identification The impact of precision;Avoid simultaneously linearized stability of the conventional Extension Kalman filtering algorithm in state vector identification process and The calculating of complicated Jacobian matrix, reduces algorithm and realizes difficulty while identification precision is improved;In PMSM drive systems, High-speed cruising area, can reduce parameter R to be identified on the premise of permanent magnet flux linkage identification precision is ensureds, reduce identification and calculate Method amount of calculation.

Description of the drawings

In order to be illustrated more clearly that the embodiment of the present invention or technical scheme of the prior art, below will be to embodiment or existing The accompanying drawing to be used needed for having technology description is briefly described, it should be apparent that, drawings in the following description are only this Some embodiments of invention, for those of ordinary skill in the art, on the premise of not paying creative work, can be with Other accompanying drawings are obtained according to these accompanying drawings.

Fig. 1 is the parameter sensitivity of the PMSM permanent magnet flux linkage identification precisions based on Unscented kalman filtering algorithm of the present invention Property analysis result.

Fig. 2 is the PMSM permanent magnet flux linkages subregion joint full rank identification based on double Unscented kalman filtering algorithms of the present invention The structured flowchart of method.

Fig. 3 is the permanent magnet flux linkage identification result of the middle and high fast Operational Zones of PMSM of the present invention.

Specific embodiment

Below in conjunction with the accompanying drawing in the embodiment of the present invention, the technical scheme in the embodiment of the present invention is carried out clear, complete Site preparation is described, it is clear that described embodiment is only a part of embodiment of the invention, rather than the embodiment of whole.It is based on Embodiment in the present invention, it is every other that those of ordinary skill in the art are obtained under the premise of creative work is not paid Embodiment, belongs to the scope of protection of the invention.

A kind of PMSM permanent magnet flux linkages full rank discrimination method, step is as follows:

Step one, stator voltage u obtained in the d-q shaftings of PMSMdAnd uq, stator current idAnd iqAnd PMSM drivetrains System rotor electrical angular velocity omegae

Stator voltage u in the d-q shaftings of the PMSMdAnd uq, stator current idAnd iqAcquisition methods have following two:

(1) stator line voltage u of sampling PMSMab、ubc, three-phase current ia、ib、ic, and obtained by coordinate transform.Coordinate Transformation matrix is respectively:

In formula, θ is flux linkage position of the rotor angle.

(2) the d-q shaft voltage command values directly calculated using PMSM driving system controllersWithReplace d-q axles Stator voltage udAnd uq, d-q shaft current command valuesReplace d-q axle stator current idAnd iq

Due to adopting electric current, rotating speed double circle structure PMSM drive systems more, therefore, method (1) needs to increase voltage adopts Sample and isolation circuit, increase system hardware expense;And method (2) using command value substitute actual value, though without the need for voltage sample with Isolation circuit, but need to consider that inverter is non-linear and sample circuit time lag caused by substitute offset issue, needed if necessary Effect compensation.

The acquisition methods of the PMSM drive systems rotor electrical angular velocity are obtained by incremental optical-electricity encoder, step Suddenly it is:

(1) in t1And t2Umber of pulse N that optical rotary encoder sends on the d-q axles of neighbouring sample instance sample PMSM1、 N2, sampling instant t1And t2Difference be sampling period T.

(2) according to rotor angular rate ωeWith optical rotary encoder impulse sampling value N1、N2And between sampling period T Relation calculate rotor angular rate ωe, its expression formula is:

In formula, M is the optical rotary encoder umber of pulse of a week, and p is permagnetic synchronous motor number of pole-pairs.

Step 2, the PMSM state equations set up in d-q shaftings, based on Unscented kalman filtering algorithm PMSM permanent magnetism is realized Body magnetic linkage is recognized, under analysis PMSM drive system difference operating conditions, stator resistance Rs, stator inductance Ld、LqChange is to permanent magnet Magnetic linkage ψfThe impact of identification precision, realizes based on the permanent magnet flux linkage ψ of Unscented kalman filtering algorithmfThe parameter of identification precision is quick Perceptual analysis;

It is used to realize that the PMSM state equations that permanent magnet flux linkage is recognized are in the d-q shaftings:

Wherein, udAnd uqD-q axle stator voltages, i are represented respectivelydAnd iqD-q axle stator currents, L are represented respectivelydAnd LqRespectively Represent d-q axle stator inductances;RsRepresent stator resistance, ψfRepresent permanent magnet flux linkage, ωeRepresent rotor electric angle frequency.

For the state equation of formula (3) description, its state vector x, input quantity u and output y are respectively:

X=[id iq ψf]T, u=[ud/Ld uq/Lq]T, y=[id iq]T (4)

It is described to realize that the method that PMSM permanent magnet flux linkages are recognized is based on Unscented kalman filtering algorithm:

The PMSM state equations for realizing permanent magnet flux linkage identification that formula (3) is described are characterized as into the one of nonlinear system As the measurement equation of form, its state equation and discretization be expressed as

In formula:X (t) be system state variables, y (tk) it is output, f () represents systematic state transfer equation, h () Represent system measuring equation, σ (t), μ (tk) respectively consider the process noise of model uncertainty and measuring uncertainty and survey Amount noise, the variance matrix of σ (t) is Q (t), μ (tk) variance matrix be R (t), u (t) be definitiveness input vector, B for control Matrix processed, is constant value matrix.

PMSM state equations, state vector, input vector, output vector are substituted in Unscented kalman filtering algorithm, it is real Now it is based on the state vector recursion of Unscented kalman filtering algorithm, you can be capable of achieving PMSM permanent magnet flux linkage ψfOn-line identification, Step is:

(1) state vector initialization

Process noise covariance matrix Q and measurement noise covariance matrix R initial values are set according to priori, and is initialized State vector x and State error variance matrix P:

(2) Sigma points are calculated

Within each sampling period (k=1,2, ∞) and Sigma points are calculated according to formula (7), obtain a n × (2n + 1) Sigma dot matrixs.

In formula:χk-1For Sigma dot matrixs,Pk-1The prediction average and covariance at a moment are expressed as, n is State vector dimension, λ=α2(n+b) it is scale parameter;α is dispersion level of the Sigma points near state variable average, is determined Sigma point distribution situations, are usually taken to be interval [10-4, 1] on little positive number;B is scale coefficient, is usually taken to be 0 or 3-n.

(3) time renewal

The transmission of Sigma points should be realized by the system state equation of discretization, as shown in formula (8), according to transmission result The prediction average and its covariance of state vector are obtained, as shown in formula (9).

Wherein,The average weight of state vector is represented,Represent the covariance weight of state vector, quantitative relation For:

(4) measurement updaue

According to measurement data, it is capable of achieving one-step prediction by formula (11)-(16) and Kalman gain updates, obtains state The optimal estimation of vector and its variance matrix.

(5) k=k+1 is made, the iteration output of state vector is realized in duplication stages (2)-(4).

Formula (5) is substituted in the Unscented kalman filtering algorithm flow described by formula (6)-(16), you can realize being based on nothing The state vector recursion of mark Kalman filtering algorithm.

It is described to realize based on the PMSM permanent magnet flux linkage ψ of Unscented kalman algorithmfThe parameters sensitivity analysis of identification precision Method be:It is determined that the PMSM stator resistance R affected by operating condition and operation temperature risesAnd d-q axle stator inductance Ld、LqChange Change scope, and analyze the ψ under different system operating mode based on Unscented kalman filtering algorithm within this rangefIdentification Errors.

Step 3, according to step 2 conclusion, in PMSM drive system low cruises area, by d-q axle stator inductance Ld、 LqCarry out joint full rank identification, and with stator resistance Rs, permanent magnet flux linkage ψfThe identification of joint full rank combines, updates each other, is distinguishing Under knowing equation full rank state, stator resistance R is eliminateds, d-q axle stator inductance Ld、LqChange is to permanent magnet flux linkage ψfIdentification precision Affect, shown in full rank identification equation such as formula (16).

Further, in the middle and high fast Operational Zone of PMSM drive systems, due to stator resistance RsChange is to permanent magnet flux linkage ψf The impact of identification precision is less, therefore using permanent magnet flux linkage ψfIdentification and d-q axle stator inductance Ld、LqThe identification of joint full rank is mutually tied The method of conjunction, eliminates d-q axle stator inductance Ld、LqChange is to permanent magnet flux linkage ψfThe impact of identification precision, full rank identification equation is such as Shown in formula (17).

According to above-mentioned thinking determine considerations PMSM Parameters variations when permanent magnet flux linkage subregion combine full rank discrimination method Structured flowchart is as shown in Figure 2.In PMSM drive system Operational Zone of the rotating speed less than 100rpm, using formula (16) Suo Shi, d-q axles are determined Sub- inductance Ld、LqJoint full rank is recognized and stator resistance Rs, permanent magnet flux linkage ψfThe full rank identification side that the identification of joint full rank combines Method, eliminates Rs、Ld、LqChange is to permanent magnet flux linkage ψfThe impact of identification precision, under the constraint of PMSM Parameters variations permanent magnet is realized Magnetic linkage ψfFull rank identification.In PMSM drive system middle and high fast Operational Zone of the rotating speed more than 100rpm, using formula (17) Suo Shi, D-q stator inductance Ld、LqJoint full rank is recognized and permanent magnet flux linkage ψfThe full rank discrimination method that identification combines, eliminates d-q stators Inductance Ld、LqChange is to permanent magnet flux linkage ψfThe impact of estimated accuracy, under the constraint of PMSM Parameters variations permanent magnet flux linkage ψ is realizedf Full rank identification.

Carried out experimental verification to the method for the present invention, experiment condition is given load torque 3Nm, and take rotating speed from 900 revs/min of dynamic processes for being down to 450 revs/min.In view of under laboratory environment, parameter of electric machine change is less in the short time, this The bright checking that aforementioned high speed scheme is completed near motor parameters value.Stator resistance RsWhen being taken as the 150% of design load, D axle stator inductance Ld, q axle stator inductance LqAnd permanent magnet flux linkage ψfIdentification result is respectively such as Fig. 3 (a), 3 (b) and Fig. 3 (c) institutes Show.With PMSM design load (Ld=Lq=0.001283H, ψf=0.1278Wb) compare, the identification precision of presently disclosed method Higher, Identification Errors average can be controlled within 7%, and in the middle and high fast Operational Zone of PMSM drive systems, without the need for considering stator Resistance RsImpact of the change to permanent magnet flux linkage identification precision, reduces parameter stator resistance R to be identifieds, effectively reduce identification Algorithm amount of calculation.

The above, the only present invention preferably specific embodiment, but protection scope of the present invention is not limited thereto, Any those familiar with the art the invention discloses technical scope in, the change or replacement that can be readily occurred in, All should be included within the scope of the present invention.

Claims (9)

1. a kind of PMSM permanent magnet flux linkages full rank discrimination method, it is characterised in that the stator voltage in the d-q shaftings of collection PMSM ud、uq, stator current id、iqAnd PMSM drive system rotor electrical angular velocity omegase, setting up in d-q shaftings is used for permanent magnet flux linkage The PMSM state equations of identification, realize that the parameter based on the PMSM permanent magnet flux linkage identification precisions of Unscented kalman filtering algorithm is quick Perceptual analysis, based on analysis result determine PMSM drive system different rotating speeds area based on double Unscented kalman filtering algorithms forever Magnet magnetic linkage full rank discrimination method, realizes that measurement noise, permagnetic synchronous motor Parameters variation and identification model are owed under order constraint Permanent magnet flux linkage high accuracy, online, full rank identification;Its step is as follows:
Stator voltage u in step one, the d-q shaftings of collection PMSMdAnd uq, stator current idAnd iqAnd PMSM drive system rotors Electric angular velocity omegae
Step 2, the PMSM state equations set up in d-q shaftings, based on Unscented kalman filtering algorithm PMSM permanent magnet magnetics are realized Chain is recognized, under analysis PMSM drive system difference operating conditions, stator resistance Rs, d-q axle stator inductance LdAnd LqChange is to permanent magnetism The impact of body magnetic linkage identification precision, realizes the parameter sensitivity based on the permanent magnet flux linkage identification precision of Unscented kalman filtering algorithm Property analysis;
Step 3, according to step 2 result, in PMSM drive system low cruises area, by d-q axle stator inductance Ld、LqCombine and distinguish Know and stator resistance Rs, permanent magnet flux linkage ψfJoint identification is combined, updated each other, under identification equation full rank state, eliminates fixed Sub- resistance Rs, d-q axle stator inductance Ld、LqImpact of the change to permanent magnet flux linkage estimated accuracy;It is middle and high in PMSM drive systems Fast Operational Zone, using permanent magnet flux linkage ψfIdentification and d-q axle stator inductance Ld、LqJoint identification combines, in identification equation full rank Under state, d-q axle stator inductance L are eliminatedd、LqImpact of the change to permanent magnet flux linkage identification precision.
2. PMSM permanent magnet flux linkages full rank discrimination method according to claim 1, it is characterised in that the d-q of the PMSM Stator voltage u in shaftingdAnd uq, stator current idAnd iqAcquisition methods be:Stator line voltage u of sampling PMSMab、ubc, three-phase Electric current ia、ib、ic, and obtained by coordinate transform, transformation matrix of coordinates is respectively:
u d u q = c o s θ s i n θ - s i n θ c o s θ 2 / 3 1 / 3 0 - 3 / 3 u a b u b c ,
i d i q = c o s θ s i n θ - s i n θ c o s θ 2 / 3 - 1 / 3 - 1 / 3 0 - 3 / 3 3 / 3 i a i b i c ;
In formula, θ is flux linkage position of the rotor angle.
3. PMSM permanent magnet flux linkages full rank discrimination method according to claim 1, it is characterised in that the d-q of the PMSM Stator voltage u in shaftingdAnd uq, stator current idAnd iqAcquisition methods be:Directly calculated using PMSM driving system controllers D-q shaft voltage command values outWithReplace d-q axle stator voltages udAnd uq, d-q shaft current command valuesReplace d- Q axle stator current idAnd iq
4. PMSM permanent magnet flux linkages full rank discrimination method according to claim 1, it is characterised in that the PMSM drivetrains The acquisition methods of system rotor electrical angular velocity are obtained by incremental optical-electricity encoder, and step is:
(1) in t1And t2Umber of pulse N that optical rotary encoder sends on the d-q axles of neighbouring sample instance sample PMSM1、N2, adopt Sample moment t1And t2Difference be sampling period T;
(2) according to rotor angular rate ωeWith optical rotary encoder impulse sampling value N1、N2And the pass between sampling period T System calculates rotor angular rate ωe, its expression formula is:
ω e = N 2 - N 1 M × T × p × 2 π ;
In formula, M is the optical rotary encoder umber of pulse of a week, and p is permagnetic synchronous motor number of pole-pairs.
5. PMSM permanent magnet flux linkages full rank discrimination method according to claim 1, it is characterised in that in the d-q shaftings For realizing that the method for PMSM state equations that permanent magnet flux linkage is recognized is:
di d d t = u d L d - R s L d i d + ω e L q L d i q di q d t = u q L q - R s L q i q - ω e L d L q i d - ω e ψ f L q dψ f d t = 0 ;
In formula, udAnd uqD-q axle stator voltages, i are represented respectivelydAnd iqD-q axle stator currents, L are represented respectivelydAnd LqRepresent respectively D-q axle stator inductances;RsRepresent stator resistance, ψfRepresent permanent magnet flux linkage, ωeRepresent rotor electric angle frequency;
In PMSM state equations, state vector is:X=[id iq ψf]T, input quantity is:U=[ud/Ld uq/Lq]T, output For:Y=[id iq]T
6. PMSM permanent magnet flux linkages full rank discrimination method according to claim 1, it is characterised in that described based on without mark card Kalman Filtering algorithm realizes that the method that PMSM permanent magnet flux linkages are recognized is:The PMSM state equations that permanent magnet flux linkage is recognized will be realized The measurement equation for being described as the general type of nonlinear system, state equation and discretization is expressed as:
x · ( t ) = f [ x ( t ) ] + B u ( t ) + σ ( t ) y ( t k ) = h [ x ( t k ) ] + μ ( t k ) ;
In formula, x (t) be system mode vector, y (tk) it is output, f () represents systematic state transfer equation, and h () is represented System measuring equation, σ (t) is the process noise for considering model uncertainty and measuring uncertainty, μ (tk) it is to consider model not The measurement noise of definitiveness and measuring uncertainty, the variance matrix of σ (t) is Q (t), μ (tk) variance matrix be R (t), u T () is definitiveness input vector, B represents control matrix;
Above-mentioned PMSM state equations, state vector, input vector, output vector are substituted in Unscented kalman filtering algorithm, it is real Now it is based on the state vector recursion of Unscented kalman filtering algorithm, you can realize the on-line identification of PMSM permanent magnet flux linkages, step For:
(1) state vector initialization
Process noise covariance matrix Q and measurement noise covariance matrix R initial values are set according to priori, init state to Amount x:Init state varivance matrix P:
(2) Sigma points are calculated
Pass through within each sampling periodCalculate Sigma Point (k=1,2, ∞) and, obtain the Sigma dot matrixs of n × (2n+1);Wherein:χk-1For Sigma dot matrixs,Pk-1It is expressed as the prediction average and covariance of a moment state vector, n is state vector dimension, λ=α2(n+ B) it is scale parameter, α is dispersion level of the Sigma points near state variable average, and b is scale coefficient;
(3) time renewal
Realize that Sigma points are transmitted by the system state equation of discretization:I.e.According to transmission As a result the prediction average and covariance of state vector are obtained:
Wherein, Wi (c)Represent the average weight of state vector, Wi (m)Table Show the covariance weight of state vector, Wi (c)And Wi (m)Quantitative relation be:
W 0 ( m ) = λ / n + λ W 0 ( c ) = λ / n + λ + ( 1 - α 2 + β ) W i ( m ) = W i ( c ) = 1 / 2 ( n + λ ) , i = 1 , 2 , ... , 2 n ;
(4) measurement updaue
According to system measurement data, by formula It is capable of achieving the one-step prediction of state vector Update and updated with Kalman gain, obtain the optimal estimation of state vector and variance matrix;Wherein, H represents calculation matrix, KkTable Show Kalman gain;
K=k+1 is made, the iteration output of state vector is realized in repeat step (2)-(4).
7. PMSM permanent magnet flux linkages full rank discrimination method according to claim 1, it is characterised in that the realization is based on nothing The method of the parameters sensitivity analysis of the PMSM permanent magnet flux linkage identification precisions of mark Kalman Algorithm is:It is determined that receiving system operation work The PMSM stator resistance R that condition and operation temperature rise affects, d axle inductances Ld, q axle inductances LqExcursion, and divide within this range Permanent magnet flux linkage Identification Errors under analysis different system operating mode based on Unscented kalman filtering algorithm.
8. PMSM permanent magnet flux linkages full rank discrimination method according to claim 1, it is characterised in that the PMSM drivetrains In system low cruise region, full rank identification equation is:
di d d t = u d L d - R s L d i d + ω e L q L q i q di q d t = u q L q - R s L q i q - ω e L d L q i d - ω e ψ f L q dψ f d t = 0 dR s d t = 0 , di d d t = u d L d - R s L d i d + ω e L q L q i q di q d t = u q L q - R s L q i q - ω e L d L q i d - ω e ψ f L q dL d d t = 0 dL q d t = 0 .
9. PMSM permanent magnet flux linkages full rank discrimination method according to claim 1, it is characterised in that the PMSM drivetrains Unite in middle and high fast operation area, full rank identification equation is:
di d d t = u d L d - R s L d i d + ω e L q L q i q di q d t = u q L q - R s L q i q - ω e L d L q i d - ω e ψ f L q dψ f d t = 0 , di d d t = u d L d - R s L d i d + ω e L q L q i q di q d t = u q L q - R s L q i q - ω e L d L q i d - ω e ψ f L q dL d d t = 0 dL q d t = 0 .
CN201610512256.XA 2016-06-29 2016-06-29 A kind of PMSM permanent magnet flux linkage full rank discrimination method CN106602952B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610512256.XA CN106602952B (en) 2016-06-29 2016-06-29 A kind of PMSM permanent magnet flux linkage full rank discrimination method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610512256.XA CN106602952B (en) 2016-06-29 2016-06-29 A kind of PMSM permanent magnet flux linkage full rank discrimination method

Publications (2)

Publication Number Publication Date
CN106602952A true CN106602952A (en) 2017-04-26
CN106602952B CN106602952B (en) 2018-12-28

Family

ID=58555736

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610512256.XA CN106602952B (en) 2016-06-29 2016-06-29 A kind of PMSM permanent magnet flux linkage full rank discrimination method

Country Status (1)

Country Link
CN (1) CN106602952B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108445401A (en) * 2018-02-09 2018-08-24 深圳市鹏诚新能源科技有限公司 On-line Estimation method, electronic device and the storage medium of battery charge state SOC
CN109039219A (en) * 2018-07-06 2018-12-18 浙江零跑科技有限公司 A kind of automobile motor guard method based on rotor magnetic steel temperature

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050007044A1 (en) * 2003-07-10 2005-01-13 Ming Qiu Sensorless control method and apparatus for a motor drive system
CN102779238A (en) * 2012-08-09 2012-11-14 北京航空航天大学 Brushless DC (Direct Current) motor system identification method on basis of adaptive Kalman filter
CN103338002A (en) * 2013-06-25 2013-10-02 同济大学 Method for identifying permanent magnet flux and quadrature axis inductance of permanent magnet synchronous motor
CN103414416A (en) * 2013-07-11 2013-11-27 中国大唐集团科学技术研究院有限公司 Permanent magnet synchronous motor sensorless vector control system based on EKF
CN104034332A (en) * 2014-06-20 2014-09-10 东南大学 Kalman filtering-based method for estimating attitude angle of rescue wrecker

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050007044A1 (en) * 2003-07-10 2005-01-13 Ming Qiu Sensorless control method and apparatus for a motor drive system
CN102779238A (en) * 2012-08-09 2012-11-14 北京航空航天大学 Brushless DC (Direct Current) motor system identification method on basis of adaptive Kalman filter
CN102779238B (en) * 2012-08-09 2015-05-27 北京航空航天大学 Brushless DC (Direct Current) motor system identification method on basis of adaptive Kalman filter
CN103338002A (en) * 2013-06-25 2013-10-02 同济大学 Method for identifying permanent magnet flux and quadrature axis inductance of permanent magnet synchronous motor
CN103414416A (en) * 2013-07-11 2013-11-27 中国大唐集团科学技术研究院有限公司 Permanent magnet synchronous motor sensorless vector control system based on EKF
CN104034332A (en) * 2014-06-20 2014-09-10 东南大学 Kalman filtering-based method for estimating attitude angle of rescue wrecker

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108445401A (en) * 2018-02-09 2018-08-24 深圳市鹏诚新能源科技有限公司 On-line Estimation method, electronic device and the storage medium of battery charge state SOC
CN109039219A (en) * 2018-07-06 2018-12-18 浙江零跑科技有限公司 A kind of automobile motor guard method based on rotor magnetic steel temperature

Also Published As

Publication number Publication date
CN106602952B (en) 2018-12-28

Similar Documents

Publication Publication Date Title
Reza et al. A review of reliable and energy efficient direct torque controlled induction motor drives
CN103401501B (en) A kind of PMSM servo system control method based on fuzzy active disturbance rejection
CN103647490B (en) A kind of sliding mode control strategy of magneto
CN102208891B (en) Method for controlling PMSM (permanent magnet synchronous motor) servo system based on friction and disturbance compensation
CN103248306B (en) Online decoupling identification method of multiple parameters of PMSM (permanent magnet synchronous motor)
CN103124158B (en) Based on the automatic setting method of the permagnetic synchronous motor speed ring controling parameters of fractional order
CN100570391C (en) The real-time detection of permanent-magnetism synchronous motor permanent magnetic field aberration and analytical approach and device thereof
Hoseinnezhad et al. Calibration of resolver sensors in electromechanical braking systems: A modified recursive weighted least-squares approach
CN103296959B (en) Permagnetic synchronous motor senseless control system and method
CN105119549B (en) A kind of motor stator resistance discrimination method
CN103312244B (en) Based on the brshless DC motor Direct Torque Control of segmented sliding moding structure
CN102385342B (en) Self-adaptation dynamic sliding mode controlling method controlled by virtual axis lathe parallel connection mechanism motion
CN101917150B (en) Robust controller of permanent magnet synchronous motor based on fuzzy-neural network generalized inverse and construction method thereof
CN201910764U (en) Permanent magnet synchronous motor (PMSM) direct torque control system based on terminal sliding mode
Xu et al. Very-low speed control of PMSM based on EKF estimation with closed loop optimized parameters
CN103338003B (en) A kind of method of electric motor load torque and inertia on-line identification simultaneously
Li et al. On the rejection of internal and external disturbances in a wind energy conversion system with direct-driven PMSG
CN104242769A (en) Permanent magnet synchronous motor speed composite control method based on continuous terminal slip form technology
CN104639003B (en) A kind of method for identification of rotational inertia of AC servo
CN104953915A (en) Permanent magnet synchronous motor sliding-mode control strategy based on novel reaching law
CN102035456A (en) Direct torque control system of permanent magnet synchronous motor based on terminal sliding mode
CN103997272B (en) The load disturbance compensation device of permagnetic synchronous motor and method
Zhu et al. Estimation of winding resistance and PM flux-linkage in brushless AC machines by reduced-order extended Kalman Filter
CN104283478B (en) A kind of Over Electric Motor with PMSM current control system and control method
CN102931906B (en) Method for asynchronous motor rotor flux linkage observation and rotation speed identification

Legal Events

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