CN102291079A - Speed sensor-free control algorithm for direct-drive permanent magnet synchronous wind power generation system - Google Patents

Speed sensor-free control algorithm for direct-drive permanent magnet synchronous wind power generation system Download PDF

Info

Publication number
CN102291079A
CN102291079A CN2011101861095A CN201110186109A CN102291079A CN 102291079 A CN102291079 A CN 102291079A CN 2011101861095 A CN2011101861095 A CN 2011101861095A CN 201110186109 A CN201110186109 A CN 201110186109A CN 102291079 A CN102291079 A CN 102291079A
Authority
CN
China
Prior art keywords
current
rotor
stator
speed
rotating speed
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.)
Pending
Application number
CN2011101861095A
Other languages
Chinese (zh)
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.)
Southeast University
Original Assignee
Southeast University
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 Southeast University filed Critical Southeast University
Priority to CN2011101861095A priority Critical patent/CN102291079A/en
Publication of CN102291079A publication Critical patent/CN102291079A/en
Pending legal-status Critical Current

Links

Images

Abstract

The invention discloses a control algorithm for estimating real-time stator resistance and speed values of a wind generator by detecting current and voltage and controlling the permanent magnet synchronous wind driven generator by taking the speed estimated by the algorithm as speed feedback. The speed value is estimated by the following steps of: online recognizing the stator resistance value of the wind driven generator based on a model reference adaptation system (MRAS) according to detected stator current and voltage values of the generator, obtaining a computational formula for an error between an estimated electric angle and an actual electric angle of a rotor based on a counter electromotive force and estimating the rotating speed of a fan by eliminating the error. The current and voltage values in the computational formula for the electric angle error are obtained according to the detection, the variable stator resistance value is recognized according to the MRAS in the running of the generator, a differential term of the stator current is realized by adopting a tracking differentiator, thereby more accurately and rapidly estimating the rotating speed.

Description

The synchronous wind power system Speedless sensor of direct-drive permanent-magnetism control algolithm
Technical field
The present invention is the Speedless sensor control algolithm of the synchronous wind power system of a kind of direct-drive permanent-magnetism, adopt the method for stator resistance identification to improve Speedless sensor control algolithm velocity estimation precision and vector control performance, belong to the technical field of wind power generation control appliance.
Background technology
In the synchronous wind power system of direct-drive permanent-magnetism, the vector control of magneto alternator (PMSG) needs the speed and the positional information of rotor, the most frequently used method is a sensor installation on armature spindle, but when the installation of transducer has reduced system reliability and increased cost, brought problem also for the maintenance of system.And general sensor for temperature, humidity, electromagnetic environment also all have requirement, and wind energy conversion system is arranged on the abominable field of environment usually, makes that general working sensor is unreliable.And PMSG rotor mechanical angular speed is very low in the more powerful directly driven wind-powered system, and traditional transducer is difficult to accurately locate electrical degree, thereby is difficult to reach fine control performance.But the PMSG number of pole-pairs is more, and motor electric angle speed can be very not low when the rotor low cruise, thereby stator voltage can be very not low yet, and the influence that the electric angle Speed identification was subjected to when mechanical separator speed was low is less; And do not require Speed identification under the zero-speed state in that system is in service.Therefore,, reduce the maintenance cost of system, in directly driven wind-powered system, generally all save encoder uniform velocity transducer, adopt the Speedless sensor control technology in order to put forward the reliability of wind-powered electricity generation unit operation.Permagnetic synchronous motor Speedless sensor control technology mainly can be divided into following 3 classes at present: a class is based on the open loop computational methods of motor ideal model, the open-loop algorithm that is based on motor mathematical model that adopts, estimation algorithm that changes as direct computing method, based on inductance and back electromotive force integration method etc.This class methods computational process is simple, can access rotor position estimation result comparatively accurately under the correct prerequisite of the parameter of electric machine.But because it need measure the current/voltage value exactly, and has very strong dependence to the parameter of electric machine, add that it is a kind of open loop Calculation Method, can't guarantee that motor is being subjected to noise jamming or still can obtains correct result when parameter changes; The another kind of closed loop algorithm that is based on various observer model, Chinese scholars has been studied the closed loop algorithm based on various observers, current application has extended Kalman filter, sliding mode observer, model reference adaptive algorithm and other adaptive algorithm etc. comparatively widely, and the permagnetic synchronous motor Speedless sensor is controlled on the stability of speed accuracy of observation and system and all is greatly improved.All these are applicable to that the no transducer control technology of high speed operation all is directly or indirectly based on extracting position signal from back-emf, and the voltage equation by motor calculates the electromotive force of being responded to and carries out the estimation of rotor-position.Because the back-emf amplitude is directly proportional with electric angle speed, when the electric angle rotating speed very low in addition when the zero-speed signal to noise ratio of back electromotive force little, other disturbing factor in addition, this class algorithm can not accurately be estimated spinner velocity and position, but more powerful direct wind-driven generator number of pole-pairs is more, even electric angle speed is not low yet when low rate start, and operates in certain range of speeds, the influence that the rotating speed identification is subjected to is less, and therefore this method can be applied in the directly driven wind-powered system; Be to serve as the algorithm based on the motor ideal characterisitics of typical case's representative with the high frequency injection method at last, in order to solve the speed observation problem under the low speed, 1993, M. Corley and R. D. Lorenz have proposed the high-frequency signal method for implanting, measure the variation of motor inductance parameters by injecting high-frequency current signal or high-frequency voltage signal, estimate the low speed and the zero-speed rotor-position of permanent magnet synchronous motor, this method utilizes the salient pole nature of salient-pole machine self or non-salient pole machine showed under high-frequency signal salient pole nature to estimate motor rotor position and rotating speed, is that current permagnetic synchronous motor does not have the effective method of mechanical transducer low cruise.From latter stage the 90's of 20th century so far, had increasing people to study the no transducer control technology of low speed and zero-speed permanent magnet synchronous motor both at home and abroad, and constantly had new method to occur.Avoided direct calculating counter electromotive force of motor based on the method for motor non-ideal characteristic, but obtained rotor position information, all obtained result preferably at motor low speed even zero-speed with the excitation of characteristic frequency.But motor needs the excitation that continues in running, reduced the voltage utilization of inverter, and wind power generator control system does not start when low speed and zero-speed.Because this class methods signal processing is comparatively complicated, dynamic property is not very good in addition.And the torque pulsation that high-frequency signal produced can cause damage to the power transmission shaft of PMSG, and is inapplicable in directly driven wind-powered system.
Summary of the invention
Variation causes the coarse problem of rotating speed identification at parameters such as motor stator resistance, the invention provides the synchronous wind power system Speedless sensor of a kind of direct-drive permanent-magnetism control algolithm.It is to contain the equation that remains to estimate the stator resistance parameter as adjustable model that the identification of MRAS stator resistance gets main thought, and as the reference model, two models have the output variable of identical physical significance with the equation that do not contain unknown parameter.Two models are worked simultaneously, and the stator resistance parameter of utilizing the difference of its output variable to come the real-time regulated adjustable model according to adaptive rate, follow the tracks of the purpose of actual resistance to reach the stator resistance identifier.Adaptive rate need satisfy stability of a system principle, and the asymptotic convergence of system and speed is guaranteed by the Popov superstability.The present invention has utilized the PMSG model parameter so that control effect and control cost and reach comparatively ideal coordination, adopts stator resistance parameter identification algorithm based on MRAS to solve parameter such as motor stator resistance and changes the coarse problem of rotating speed identification that causes.The stator resistance parameter that utilization picks out obtains rotor based on back electromotive force and estimates electrical degree and actual electrical angle calculation formula, and the stator current differential term in the formula has then adopted tracking differentiator (TD) to realize.Estimate that by eliminating rotor electrical degree and rotor actual electrical angular error estimate the wind energy conversion system rotating speed, utilize the estimation rotating speed wind-driven generator to be realized closed-loop vector control as feedback.
Technical solution of the present invention can be divided into following step and realize:
(1) current/voltage coordinate transform, the aerogenerator stator current value and magnitude of voltage and the estimated motor electrical degree that obtain according to detection are that the benchmark of rotational coordinates carries out the dq rotating coordinate transformation to stator current and voltage with the axis of rotor permanent magnetic potential.
(2) stator resistance identification utilizes the current/voltage value behind the rotating coordinate transformation, based on MRAS the motor stator resistance parameter of real-time change is picked out.
(3) rotating speed is estimated, be used to receive the motor stator resistance parameter of current value after the coordinate transform and magnitude of voltage, on-line identification, obtain rotor based on back electromotive force and estimate electrical degree and actual electrical angle calculation formula, the current/voltage value obtains according to detection in the computing formula, identification obtains the stator resistance value according to MRAS, and the stator current differential term has then adopted tracking differentiator (TD) to realize.Estimate that by eliminating rotor electrical degree and rotor actual electrical angular error estimate the wind energy conversion system rotating speed, utilize the estimation rotating speed wind-driven generator to be realized closed-loop vector control as feedback.
Advantage of the present invention: if the rotating speed discrimination method does not carry out on-line identification to the parameter of electric machine, then can rotating speed occur because of variations in temperature and the identification of rotor space position is inaccurate, thereby influence the closed-loop drive performance of PMSG.The present invention proposes a kind of based on MRAS permanent-magnetic synchronous motor stator resistance identification scheme, be applied in the synchronous wind power control system of high-power direct-drive permanent-magnetism, for raising the efficiency, reference model selects motor body, adjustable model selects voltage model, the stator resistance value adopts the PI adaptive rate to reach the accurate estimation of stator resistance as the adjustable parameter of adjustable model.After stator resistance on-line identification goes out, calculate two components of back electromotive force under the dq rotational coordinates based on motor mathematical model, calculate estimation electrical degree and actual electrical angular error based on the parameter of electric machine of measuring and estimate, eliminate this error by the PI control law and extract rotor-position and tach signal.The present invention adopts the stator resistance parameter identification algorithm based on MRAS to solve parameter such as motor stator resistance to change the coarse problem of rotating speed identification that causes.Thereby weakened the influence of image parameter variation to a great extent to the rotating speed identification.The stability and the reliability of Speedless sensor control algolithm have been improved.
Description of drawings
Fig. 1 stator resistance identification algorithm block diagram
Fig. 2 Speedless sensor algorithm block diagram
Among the figure 1 is magneto alternator, the 2nd, electric current coordinate converter, the 3rd, voltage coordinate converter, the 4th, stator current adjustable model, the 5th, adaptive control laws, the 6th, back electromotive force estimation, the 7th, error angle estimation, the 8th, PI controller, the 9th, low pass filter, the 10th, integrator.
Symbol and variable declaration among the figure
Figure 920437DEST_PATH_IMAGE001
, ,
Figure 747808DEST_PATH_IMAGE003
The stator three-phase voltage;
Figure 143017DEST_PATH_IMAGE004
,
Figure 895072DEST_PATH_IMAGE005
, The stator three-phase current;
Figure DEST_PATH_IMAGE007A
,
Figure 296153DEST_PATH_IMAGE008
Stator voltage on the dq coordinate system;
,
Figure 619687DEST_PATH_IMAGE010
Stator current on the dq coordinate system;
Figure DEST_PATH_IMAGE011A
,
Figure 413199DEST_PATH_IMAGE012
The stator dq electric current estimated value of adjustable model;
The stator resistance estimated value of adjustable model;
Figure DEST_PATH_IMAGE014AA
The estimation rotating speed;
Figure DEST_PATH_IMAGE015A
The rotor estimation error;
Figure 223209DEST_PATH_IMAGE016
The estimated rotor position.
Embodiment
Analyze with the non-salient pole permagnetic synchronous motor, the stator current Mathematical Modeling under rotating coordinate system is
Figure DEST_PATH_IMAGE017A
(1)
The estimation stator current Mathematical Modeling that constructing variable is adjustable is
(2)
Figure DEST_PATH_IMAGE007AA
, Be the stator voltage on the dq coordinate system;
Figure DEST_PATH_IMAGE009AA
,
Figure 62036DEST_PATH_IMAGE010
Be the stator current on the dq coordinate system;
Figure DEST_PATH_IMAGE011AA
,
Figure 429564DEST_PATH_IMAGE019
Stator dq electric current estimated value for adjustable model;
Figure DEST_PATH_IMAGE007AAA
Be the stator equivalent inductance on the dq coordinate system;
Be stator resistance;
Figure DEST_PATH_IMAGE021A
Stator resistance estimated value for adjustable model;
Be the rotor permanent magnet magnetic linkage;
Figure DEST_PATH_IMAGE023A
Be differential operator;
Figure DEST_PATH_IMAGE024AA
Be rotor electric angle speed;
Definition generalized error vector is
Figure DEST_PATH_IMAGE025A
Figure 294358DEST_PATH_IMAGE026
Figure 569482DEST_PATH_IMAGE027
Formula (2) deducts formula (1),
Figure DEST_PATH_IMAGE028AAA
(3)
In the formula
Figure DEST_PATH_IMAGE029A
,
Get D= I, then v= De= e
By the overstable theory of Popov, system stability need be satisfied:
(1) linear element Transitive Matrices G( s)= D( SI AM) -1Just real for strictness, this proves easily.
(2) the nonlinear time-varying link satisfies the Popov integral inequality, needs to select the suitable parameters adaptive rate to satisfy.
Figure 803246DEST_PATH_IMAGE031
(4)
Figure 719118DEST_PATH_IMAGE032
Be a limited positive constant, do not rely on
Figure 797933DEST_PATH_IMAGE033
Will v, ω substitution integral inequality,
Figure DEST_PATH_IMAGE034A
(5)
Press the common configuration of model reference adaptive rate,
Figure DEST_PATH_IMAGE035A
Get following proportional integral form
Figure DEST_PATH_IMAGE036A
(6)
Formula (6) substitution formula (5) can get
Figure DEST_PATH_IMAGE037A
(7)
Figure DEST_PATH_IMAGE038A
(8)
Figure DEST_PATH_IMAGE039A
(9)
Want
Figure DEST_PATH_IMAGE040A
Set up, can select
Figure 732129DEST_PATH_IMAGE041
, ,
Figure 579048DEST_PATH_IMAGE043
,
Figure 461553DEST_PATH_IMAGE044
For with
Figure 988DEST_PATH_IMAGE045
The limited positive number that same nature is arranged.
For formula (8), get
Figure DEST_PATH_IMAGE046A
,?
Figure 62310DEST_PATH_IMAGE047
>0 (10)
Then
Figure DEST_PATH_IMAGE049A
(11)
To formula (10) both sides differentiate:
Figure 983179DEST_PATH_IMAGE048
>0 (12)
To formula (9), desirable ,
Figure 580382DEST_PATH_IMAGE052
0, make
Figure 880782DEST_PATH_IMAGE042
Set up.
According to selected adaptive law, can guarantee
Figure DEST_PATH_IMAGE040AA
Set up, can get
Figure 930647DEST_PATH_IMAGE053
Identification algorithm is
Figure DEST_PATH_IMAGE054A
Figure 279589DEST_PATH_IMAGE055
Figure 278769DEST_PATH_IMAGE056
(13)
Stator resistance identification algorithm block diagram such as Fig. 1.
Because the rotor-position of control system without position sensor is not to obtain by actual detected, and is based on the result of stator resistance identifier
Figure 563120DEST_PATH_IMAGE057
Construct that motor mathematical model estimates, so in the actual rotor position
Figure 162597DEST_PATH_IMAGE058
And estimated position
Figure 190596DEST_PATH_IMAGE016
Between can have error
Figure 309862DEST_PATH_IMAGE059
Suppose that motor rotor position and the physical location estimated differ
Figure 686485DEST_PATH_IMAGE060
, at this moment, the back electromotive force that is produced by rotor permanent magnet has produced in the coordinate with the estimated rotor position orientation
Figure DEST_PATH_IMAGE061A
Two components of coordinate system, promptly
Figure 523991DEST_PATH_IMAGE062
With If the directed coordinate of actual rotor position is
Figure DEST_PATH_IMAGE064A
Coordinate system, and based on estimated position
Figure 709214DEST_PATH_IMAGE016
Be
Figure 69788DEST_PATH_IMAGE065
Coordinate system.
Can obtain based on the estimation rotor-position by coordinate transform With the estimation rotating speed The permagnetic synchronous motor Mathematical Modeling
Figure DEST_PATH_IMAGE068AAA
(14)
Figure DEST_PATH_IMAGE061AA
Back electromotive force under the coordinate system is:
Figure DEST_PATH_IMAGE069A
(15)
(16)
Figure DEST_PATH_IMAGE071A
It is back EMF coefficient.Suppose to estimate that rotating speed and actual speed error are enough little, then back electromotive force is exactly a site error function.Calculate the rotor estimation error
Figure DEST_PATH_IMAGE015AA
:
(17)
In the rotor estimation error
Figure DEST_PATH_IMAGE015AAA
The differential term that pair electric current is arranged in the computing formula (17) as can be seen.Differential is realized relatively difficulty in practice.Differential term can be ignored in invariablenes turning speed and load variations that does not have and disturbance, but change in reference rotation velocity and to cause current value to change, perhaps when changing, applied load also can cause current value to change, this moment, dq shaft current differential term had to consider, otherwise dynamic property and disturbance rejection ability that rotating speed is estimated will descend.And in wind power control system, change of wind velocity is frequent, thereby the blower fan torque that wind-force forms is all changing constantly, simultaneously blower fan also need according to wind speed adjust dq shaft current value to change wind speed so that the wind energy of catching reaches best, so the speed dynamic of wind-driven generator is estimated very important.To this, adopt and follow the tracks of the differential that differentiator is realized dq shaft current item.
The main thought that d axle or q shaft current are followed the tracks of differentiator (TD) is: with d axle or q shaft current
Figure DEST_PATH_IMAGE073A
Be input signal, it will export 2 signals
Figure DEST_PATH_IMAGE074A
With
Figure 841357DEST_PATH_IMAGE075
, wherein Follow the tracks of
Figure 466690DEST_PATH_IMAGE077
, and
Figure 75526DEST_PATH_IMAGE078
=
Figure 435969DEST_PATH_IMAGE079
Thereby, As " approximate differential ".Appropriate structuring is followed the tracks of the internal structure of differentiator, signal
Figure DEST_PATH_IMAGE074AA
,
Figure 573056DEST_PATH_IMAGE075
Can be at input end signal discontinuous or extract continuous signal and differential signal thereof preferably when comprising certain noise, and speed is fast, the precision height.The dq shaft current is followed the tracks of differentiator and is adopted the way of realization of discrete system time-optimal control comprehensive function algorithm as follows:
(18)
Figure 336798DEST_PATH_IMAGE083
(19)
In formula (18), (19)
Figure DEST_PATH_IMAGE084A
Be the sampling period,
Figure 600289DEST_PATH_IMAGE085
Be that the speed density factor is followed the tracks of in reflection.The tracking differentiator of this form all has effect well at aspects such as tracking performance, differential quality and elimination electric shocks.Can obtain dq shaft current differential signal with the method.In the high-power wind turbine system, blower fan is in the permanent environment that becomes of a wind speed, current value is not only very big, the current transients value also may be very big, and change very frequent, dq shaft current differential signal molecule denominator proportion in formula (17) arctan function can not be ignored, and accurately follows the tracks of the rotation speed of fan variation fast and and must obtain dq shaft current differential signal, has realized dq shaft current differential function well and follow the tracks of differentiator.
The rotor estimation error
Figure DEST_PATH_IMAGE015AAAA
All can obtain after the estimation of each parameter by detection or preceding method in the formula, thereby can estimate error angle
Figure DEST_PATH_IMAGE015AAAAA
, by right
Figure 246034DEST_PATH_IMAGE086
Carry out obtaining estimating rotating speed after the proportional integral adjusting , to estimating that thereby rotating speed carries out integration and obtains the estimated rotor position
Figure DEST_PATH_IMAGE066AA
Bigger by the fluctuation of speed that calculates, therefore output adds the one-level low pass filter usually.Speedless sensor algorithm structure such as Fig. 2.
The transfer function that can derive this algorithm by its principle is:
Figure DEST_PATH_IMAGE087A
(20)
This is a second-order system, can draw the PI parameter according to response characteristic.

Claims (4)

1. the synchronous wind power system Speedless sensor of direct-drive permanent-magnetism control algolithm may further comprise the steps:
(1) current/voltage coordinate transform, the aerogenerator stator current value and magnitude of voltage and the estimated motor electrical degree that obtain according to detection are that the benchmark of rotational coordinates carries out the dq rotating coordinate transformation to stator current and voltage with the axis of rotor permanent magnetic potential.
(2) stator resistance identification utilizes the current/voltage value behind the rotating coordinate transformation, based on MRAS the motor stator resistance parameter of real-time change is picked out.
(3) rotating speed is estimated, be used to receive the motor stator resistance parameter of current value after the coordinate transform and magnitude of voltage, on-line identification, obtain rotor based on back electromotive force and estimate electrical degree and actual electrical angle calculation formula, the current/voltage value obtains according to detection in the computing formula, identification obtains the stator resistance value according to MRAS, and the stator current differential term has then adopted tracking differentiator (TD) to realize.Estimate that by eliminating rotor electrical degree and rotor actual electrical angular error estimate the wind energy conversion system rotating speed, utilize the estimation rotating speed wind-driven generator to be realized closed-loop vector control as feedback.
2. the synchronous wind power system Speedless sensor of direct-drive permanent-magnetism according to claim 1 control algolithm, it is characterized in that described step (1) is for sending power-on command when wind power system, beginning was from detecting after system received the current transformer power-on command, start the net side converter, temporary transient starter motor side converter not, generator unit stator has only magnitude of voltage, according to the magnitude of voltage that records by phase-locked electrical degree and the electric angle speed that calculates wind-driven generator.According to electrical degree and the electric angle speed calculated, the generator converter controller adopts Speedless sensor closed-loop vector control method starter motor side converter, begin to detect simultaneously the input as coordinate converter of stator voltage and electric current after the startup, the value after the conversion is again as the input of MRAS stator resistance identifier and rotating speed estimator.
3. the synchronous wind power system Speedless sensor of direct-drive permanent-magnetism according to claim 1 control algolithm is characterized in that beginning to start based on MRAS permanent-magnetic synchronous motor stator resistance identification scheme after described step (2) is for the wind-driven generator operation.Reference model selects motor body, and adjustable model selects voltage model, and the stator resistance value adopts the PI adaptive rate to reach the accurate estimation of stator resistance as the adjustable parameter of adjustable model.Analyze with the non-salient pole permagnetic synchronous motor, the stator current Mathematical Modeling under rotating coordinate system is
p i d i q = - R s L s ω - ω - R s L s i d i q + u d u q - ψ f · ω - - - ( 1 )
The estimation stator current Mathematical Modeling that constructing variable is adjustable is
p i ^ d i ^ q = - R ^ s L s ω - ω - R ^ s L s i ^ d i ^ q + u d u q - ψ f · ω - - - ( 2 )
u d, u qBe the stator voltage on the dq coordinate system;
i d, i qBe the stator current on the dq coordinate system;
Figure FDA0000073735770000013
Stator dq electric current estimated value for adjustable model;
L sBe the stator equivalent inductance on the dq coordinate system;
R sBe stator resistance;
Figure FDA0000073735770000021
Stator resistance estimated value for adjustable model;
ψ fBe the rotor permanent magnet magnetic linkage;
P is a differential operator;
ω is a rotor electric angle speed.
Definition generalized error vector is
e = i - i ^ = i d - i ^ d i q - i ^ q , i = i d i q , i ^ = i ^ d i ^ q
Press the common configuration of model reference adaptive rate, select adaptive law for use according to the overstable theory of Popov,
Figure FDA0000073735770000025
Identification algorithm is
R ^ s = k i ∫ 0 t ( i ^ d 2 - i d i ^ d + i ^ q 2 - i q i ^ q ) dτ + k p ( i ^ d 2 - i d i ^ d + i ^ q 2 - i q i ^ q ) + R ^ s ( 0 ) - - - ( 3 )
4. the synchronous wind power system Speedless sensor of the direct-drive permanent-magnetism control algolithm based on the identification of stator resistance according to claim 1, it is characterized in that described step (3) rotating speed estimates to utilize the motor stator resistance parameter of current value after the coordinate transform and magnitude of voltage, on-line identification, estimate that by eliminating rotor electrical degree and rotor actual electrical angular error estimate the wind energy conversion system rotating speed.Suppose the rotor estimated position
Figure FDA0000073735770000027
θ differs θ with the actual rotor position e, the back electromotive force that is produced by rotor permanent magnet has produced in the coordinate with the estimated rotor position orientation
Figure FDA0000073735770000028
Two components of coordinate system With
Figure FDA00000737357700000210
If the directed coordinate of actual rotor position is d qCoordinate system, and based on estimated position
Figure FDA00000737357700000211
Be Coordinate system.Can obtain based on the estimation rotor-position by coordinate transform
Figure FDA00000737357700000213
With the estimation rotating speed
Figure FDA00000737357700000214
The permagnetic synchronous motor Mathematical Modeling
u ^ d = R ^ s i ^ d + L s p i ^ d - ω ^ L s i ^ q + e ^ d u ^ q = R ^ s i ^ q + L s p i ^ q + ω ^ L s i ^ d + e ^ q - - - ( 4 )
Figure FDA00000737357700000216
Back electromotive force under the coordinate system is:
e ^ d = - E ex sin θ e - ( ω ^ - ω ) L s i q e ^ q = E ex cos θ e + ( ω ^ - ω ) L s i d - - - ( 5 )
E ex=ω·ψ f (6)
E ExIt is back EMF coefficient.Suppose to estimate that rotating speed and actual speed error are enough little, then back electromotive force is exactly a site error function.Calculate the rotor estimation error
Figure FDA0000073735770000031
θ ^ e = tan - 1 ( - e ^ d e ^ q )
= tan - 1 - ( u ^ d - R ^ s i ^ d + L s ω ^ i ^ q - L s p i ^ d ) u ^ q - R ^ s i ^ q - L s ω ^ i ^ d - L s p i ^ q - - - ( 7 )
In the rotor estimation error
Figure FDA0000073735770000034
Computing formula has the current differential item in (7).In wind power control system, change in reference rotation velocity and to cause current value to change, perhaps also can cause current value to change when applied load changes, this moment, dq shaft current differential term can not be ignored, otherwise dynamic property and disturbance rejection ability that rotating speed is estimated will descend.Adopt and follow the tracks of the differential that differentiator is realized dq shaft current item.The structure that d axle or q shaft current are followed the tracks of differentiator (TD) is: d axle or q shaft current i (t) are input signal, and it will export 2 signal x 1And x 2, x wherein 1Follow the tracks of i (t), and
Figure FDA0000073735770000035
Thereby x 2" approximate differential " as i (t).Appropriate structuring is followed the tracks of the internal structure of differentiator, signal x 1, x 2Can be at input end signal discontinuous or extract continuous signal and differential signal thereof preferably when comprising certain noise, and speed is fast, the precision height.The dq shaft current is followed the tracks of differentiator and is adopted the way of realization of discrete system time-optimal control comprehensive function algorithm as follows:
x 1 ( k + 1 ) = x 1 ( k ) + hx 2 ( k ) x 2 ( k + 1 ) = x 2 ( k ) + hf , | f | ≤ r f = fhan ( x 1 ( k ) , x 2 ( k ) , r , h ) - - - ( 8 )
d = rh d 0 = hd y = x 1 + h x 2 d 0 = d 2 + 8 r | y | a = x 2 + ( a 0 - d ) 2 · sign ( y ) , | y | > d 0 x 2 + y h , | y | ≤ d 0 fhan = r · sign ( a ) , | a | > d r a d , | a | ≤ d - - - ( 9 )
H is the sampling period in formula (8), (9), and r is that the speed density factor is followed the tracks of in reflection.Obtain dq shaft current differential signal with the method.
The rotor estimation error
Figure FDA0000073735770000038
Obtain after the estimation of each parameter by detection or preceding method in the formula, thereby estimate error angle
Figure FDA0000073735770000039
By right
Figure FDA00000737357700000310
Carry out obtaining the initial estimation rotating speed after the PI adjusting, finally estimated rotating speed through the one-level low pass filter
Figure FDA00000737357700000311
To estimating that thereby rotating speed carries out integration and obtains the estimated rotor position
Figure FDA0000073735770000041
CN2011101861095A 2011-07-05 2011-07-05 Speed sensor-free control algorithm for direct-drive permanent magnet synchronous wind power generation system Pending CN102291079A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2011101861095A CN102291079A (en) 2011-07-05 2011-07-05 Speed sensor-free control algorithm for direct-drive permanent magnet synchronous wind power generation system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2011101861095A CN102291079A (en) 2011-07-05 2011-07-05 Speed sensor-free control algorithm for direct-drive permanent magnet synchronous wind power generation system

Publications (1)

Publication Number Publication Date
CN102291079A true CN102291079A (en) 2011-12-21

Family

ID=45337235

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2011101861095A Pending CN102291079A (en) 2011-07-05 2011-07-05 Speed sensor-free control algorithm for direct-drive permanent magnet synchronous wind power generation system

Country Status (1)

Country Link
CN (1) CN102291079A (en)

Cited By (37)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102420561A (en) * 2011-12-01 2012-04-18 国电南京自动化股份有限公司 Speed sensorless vector control method on basis of cascaded high voltage inverter
CN102857161A (en) * 2012-09-24 2013-01-02 海信(山东)空调有限公司 Air-conditioning apparatus and permanent magnet synchronous motor control method and device of air-conditioning apparatus
CN103051280A (en) * 2012-12-19 2013-04-17 厦门理工学院 Method for estimating rotating speed and rotor position of low-speed directly-driven permanent magnet synchronous generator
CN103259481A (en) * 2012-02-15 2013-08-21 通用汽车环球科技运作有限责任公司 Method and system for estimating electrical angular speed of permanent magnet machine
CN103439657A (en) * 2013-07-23 2013-12-11 南京工程学院 AC servo motor drive parameter detection method and application of AC servo motor drive parameter detection method in fault detection
CN103701388A (en) * 2013-11-06 2014-04-02 国家电网公司 Selection method for speed identification strategy of permanent magnet direct drive wind power system
CN104009696A (en) * 2014-05-08 2014-08-27 昆明理工大学 Interactive model reference adaptive speed and stator resistance identification method based on sliding-mode control
CN104639008A (en) * 2013-11-08 2015-05-20 Lg电子株式会社 Motor driving apparatus and laundry treatment machine including the same
CN104682810A (en) * 2013-11-29 2015-06-03 长沙市比亚迪汽车有限公司 Method for compensating current waveform of permanent magnet synchronous motor
CN105553370A (en) * 2015-12-16 2016-05-04 四川长虹电器股份有限公司 Motor control method and electronic equipment
CN106059435A (en) * 2016-07-04 2016-10-26 南京航空航天大学 Method for improving estimation precision of rotor position of permanent magnet synchronous motor
CN106059436A (en) * 2016-07-04 2016-10-26 南京航空航天大学 Method of reducing errors in estimating position and rotating speed of permanent magnet synchronous motor rotor
CN106096298A (en) * 2016-06-20 2016-11-09 重庆大学 Transducer parameters on-line identification method based on observer
US9708746B2 (en) 2013-11-08 2017-07-18 Lg Electronics Inc. Motor driving apparatus and laundry treatment machine including the same
US9708745B2 (en) 2013-11-08 2017-07-18 Lg Electronics Inc. Motor driving apparatus and laundry treatment machine including the same
US9765468B2 (en) 2013-11-08 2017-09-19 Lg Electronics Inc. Seoul National University R&Db Foundation Motor driving apparatus and laundry treatment machine including the same
CN107769655A (en) * 2017-11-06 2018-03-06 广东美的暖通设备有限公司 Permagnetic synchronous motor turn count method, apparatus, computing device and storage medium
CN108540031A (en) * 2018-04-24 2018-09-14 南京信息职业技术学院 The method for estimating rotating speed and moment controlling system of bearingless synchronous reluctance motor
CN109596241A (en) * 2018-11-28 2019-04-09 同济大学 The online Temperature fusion estimation method of rotor and system
CN109742992A (en) * 2019-01-22 2019-05-10 南通大学 A kind of new type composite excitation synchronous motor control method of position-sensor-free
CN109802604A (en) * 2019-01-30 2019-05-24 深圳市泰道精密机电有限公司 A kind of magneto driving method and system
CN109981018A (en) * 2017-12-27 2019-07-05 上海大郡动力控制技术有限公司 The starting of position-sensor-free permanent magnet synchronous motor and vector control method
CN110380656A (en) * 2019-06-03 2019-10-25 中国人民解放军陆军装甲兵学院 A kind of permanent magnet synchronous motor method for controlling position-less sensor based on simplified MRAS method
CN110417320A (en) * 2019-07-04 2019-11-05 东南大学 A kind of vertical operation magnetic flux switching permanent-magnetism linear motor uplink and downlink position-sensorless control method
CN110635738A (en) * 2019-10-12 2019-12-31 东南大学 Real-time identification method for stator resistance and motor temperature of permanent magnet synchronous motor
CN110995093A (en) * 2019-12-05 2020-04-10 北京动力机械研究所 Rotor position and rotating speed detection method based on back emf open loop estimation
CN111543003A (en) * 2018-01-12 2020-08-14 三菱电机株式会社 Control device for rotary machine
CN111711398A (en) * 2020-07-13 2020-09-25 华中科技大学 Dynamic performance improvement method for permanent magnet synchronous motor position sensorless control system
CN112366989A (en) * 2020-11-19 2021-02-12 北京信息科技大学 Brushless direct current motor control method based on parameter identification
CN112886885A (en) * 2021-03-26 2021-06-01 上海电气风电集团股份有限公司 Permanent magnet synchronous generator control method, device and system of wind generating set
CN112968644A (en) * 2021-03-08 2021-06-15 上海交通大学 Permanent magnet synchronous motor parameter online identification method, system, terminal and medium
CN113346810A (en) * 2021-06-18 2021-09-03 湖南科技大学 Speed and current double closed-loop fuzzy control PMSM sensorless control method
CN113965126A (en) * 2021-10-25 2022-01-21 深圳市振邦智能科技股份有限公司 Rotor initial position detection method based on current differential method
CN114204867A (en) * 2021-12-16 2022-03-18 上海氢恒汽车电子有限公司 Rotor position determination method of permanent magnet synchronous motor
CN114337402A (en) * 2021-12-08 2022-04-12 燕山大学 Synchronous reluctance motor inductance parameter identification and rotor position hybrid estimation method
CN115864929A (en) * 2022-12-27 2023-03-28 江西安百川电气有限公司 Method and device for detecting rotor step loss of permanent magnet synchronous motor
CN116317770A (en) * 2023-02-03 2023-06-23 北京中科昊芯科技有限公司 Method for identifying motor stator resistance offline, method for controlling motor and medium

Cited By (63)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102420561B (en) * 2011-12-01 2013-07-10 国电南京自动化股份有限公司 Speed sensorless vector control method on basis of cascaded high voltage inverter
CN102420561A (en) * 2011-12-01 2012-04-18 国电南京自动化股份有限公司 Speed sensorless vector control method on basis of cascaded high voltage inverter
CN103259481A (en) * 2012-02-15 2013-08-21 通用汽车环球科技运作有限责任公司 Method and system for estimating electrical angular speed of permanent magnet machine
CN103259481B (en) * 2012-02-15 2016-07-06 通用汽车环球科技运作有限责任公司 The method and system of the angular rate of estimation permanent magnet motor
CN102857161A (en) * 2012-09-24 2013-01-02 海信(山东)空调有限公司 Air-conditioning apparatus and permanent magnet synchronous motor control method and device of air-conditioning apparatus
CN103051280A (en) * 2012-12-19 2013-04-17 厦门理工学院 Method for estimating rotating speed and rotor position of low-speed directly-driven permanent magnet synchronous generator
CN103439657B (en) * 2013-07-23 2016-05-11 南京康尼机电股份有限公司 AC servo motor transmission parameter detection method and the application in fault detect thereof
CN103439657A (en) * 2013-07-23 2013-12-11 南京工程学院 AC servo motor drive parameter detection method and application of AC servo motor drive parameter detection method in fault detection
CN103701388A (en) * 2013-11-06 2014-04-02 国家电网公司 Selection method for speed identification strategy of permanent magnet direct drive wind power system
CN103701388B (en) * 2013-11-06 2017-01-11 国家电网公司 Selection method for speed identification strategy of permanent magnet direct drive wind power system
CN104639008A (en) * 2013-11-08 2015-05-20 Lg电子株式会社 Motor driving apparatus and laundry treatment machine including the same
US9765468B2 (en) 2013-11-08 2017-09-19 Lg Electronics Inc. Seoul National University R&Db Foundation Motor driving apparatus and laundry treatment machine including the same
CN104639008B (en) * 2013-11-08 2017-06-20 Lg电子株式会社 Motor driver and the washings processing equipment with the motor driver
US9708746B2 (en) 2013-11-08 2017-07-18 Lg Electronics Inc. Motor driving apparatus and laundry treatment machine including the same
US9708744B2 (en) 2013-11-08 2017-07-18 Lg Electronics Inc. Motor driving apparatus and laundry treatment machine including the same
US9708745B2 (en) 2013-11-08 2017-07-18 Lg Electronics Inc. Motor driving apparatus and laundry treatment machine including the same
CN104682810A (en) * 2013-11-29 2015-06-03 长沙市比亚迪汽车有限公司 Method for compensating current waveform of permanent magnet synchronous motor
CN104682810B (en) * 2013-11-29 2018-09-07 长沙市比亚迪汽车有限公司 A method of for compensating permanent magnet synchronous motor current waveform
CN104009696A (en) * 2014-05-08 2014-08-27 昆明理工大学 Interactive model reference adaptive speed and stator resistance identification method based on sliding-mode control
CN104009696B (en) * 2014-05-08 2017-02-08 昆明理工大学 Interactive model reference adaptive speed and stator resistance identification method based on sliding-mode control
CN105553370A (en) * 2015-12-16 2016-05-04 四川长虹电器股份有限公司 Motor control method and electronic equipment
CN105553370B (en) * 2015-12-16 2018-12-25 四川长虹电器股份有限公司 A kind of motor control method and electronic equipment
CN106096298B (en) * 2016-06-20 2019-02-15 重庆大学 Transducer parameters on-line identification method based on observer
CN106096298A (en) * 2016-06-20 2016-11-09 重庆大学 Transducer parameters on-line identification method based on observer
CN106059436B (en) * 2016-07-04 2018-08-28 南京航空航天大学 A method of reducing permanent-magnet synchronous motor rotor position and error for rotating speed estimation
CN106059435B (en) * 2016-07-04 2018-12-07 南京航空航天大学 A method of improving permanent-magnet synchronous motor rotor position estimated accuracy
CN106059436A (en) * 2016-07-04 2016-10-26 南京航空航天大学 Method of reducing errors in estimating position and rotating speed of permanent magnet synchronous motor rotor
CN106059435A (en) * 2016-07-04 2016-10-26 南京航空航天大学 Method for improving estimation precision of rotor position of permanent magnet synchronous motor
CN107769655A (en) * 2017-11-06 2018-03-06 广东美的暖通设备有限公司 Permagnetic synchronous motor turn count method, apparatus, computing device and storage medium
CN107769655B (en) * 2017-11-06 2020-02-28 广东美的暖通设备有限公司 Method and device for estimating rotating speed of permanent magnet synchronous motor, computing equipment and storage medium
CN109981018B (en) * 2017-12-27 2021-01-29 上海大郡动力控制技术有限公司 Starting and vector control method of permanent magnet synchronous motor without position sensor
CN109981018A (en) * 2017-12-27 2019-07-05 上海大郡动力控制技术有限公司 The starting of position-sensor-free permanent magnet synchronous motor and vector control method
CN111543003B (en) * 2018-01-12 2023-12-12 三菱电机株式会社 Control device for rotary machine
CN111543003A (en) * 2018-01-12 2020-08-14 三菱电机株式会社 Control device for rotary machine
CN108540031A (en) * 2018-04-24 2018-09-14 南京信息职业技术学院 The method for estimating rotating speed and moment controlling system of bearingless synchronous reluctance motor
CN109596241A (en) * 2018-11-28 2019-04-09 同济大学 The online Temperature fusion estimation method of rotor and system
CN109596241B (en) * 2018-11-28 2020-07-14 同济大学 Online temperature fusion estimation method for motor rotor
CN109742992A (en) * 2019-01-22 2019-05-10 南通大学 A kind of new type composite excitation synchronous motor control method of position-sensor-free
CN109742992B (en) * 2019-01-22 2022-06-17 南通大学 Novel hybrid excitation synchronous motor control method without position sensor
CN109802604A (en) * 2019-01-30 2019-05-24 深圳市泰道精密机电有限公司 A kind of magneto driving method and system
CN110380656A (en) * 2019-06-03 2019-10-25 中国人民解放军陆军装甲兵学院 A kind of permanent magnet synchronous motor method for controlling position-less sensor based on simplified MRAS method
CN110417320A (en) * 2019-07-04 2019-11-05 东南大学 A kind of vertical operation magnetic flux switching permanent-magnetism linear motor uplink and downlink position-sensorless control method
CN110417320B (en) * 2019-07-04 2021-06-01 东南大学 Up-down no-position control method for vertical operation magnetic flux switching permanent magnet linear motor
CN110635738B (en) * 2019-10-12 2020-12-18 东南大学溧阳研究院 Real-time identification method for stator resistance and motor temperature of permanent magnet synchronous motor
CN110635738A (en) * 2019-10-12 2019-12-31 东南大学 Real-time identification method for stator resistance and motor temperature of permanent magnet synchronous motor
CN110995093A (en) * 2019-12-05 2020-04-10 北京动力机械研究所 Rotor position and rotating speed detection method based on back emf open loop estimation
CN111711398A (en) * 2020-07-13 2020-09-25 华中科技大学 Dynamic performance improvement method for permanent magnet synchronous motor position sensorless control system
CN111711398B (en) * 2020-07-13 2022-02-18 华中科技大学 Dynamic performance improvement method for permanent magnet synchronous motor position sensorless control system
CN112366989A (en) * 2020-11-19 2021-02-12 北京信息科技大学 Brushless direct current motor control method based on parameter identification
CN112968644A (en) * 2021-03-08 2021-06-15 上海交通大学 Permanent magnet synchronous motor parameter online identification method, system, terminal and medium
CN112886885B (en) * 2021-03-26 2022-07-29 上海电气风电集团股份有限公司 Permanent magnet synchronous generator control method, device and system of wind generating set
CN112886885A (en) * 2021-03-26 2021-06-01 上海电气风电集团股份有限公司 Permanent magnet synchronous generator control method, device and system of wind generating set
CN113346810A (en) * 2021-06-18 2021-09-03 湖南科技大学 Speed and current double closed-loop fuzzy control PMSM sensorless control method
CN113346810B (en) * 2021-06-18 2022-11-11 湖南科技大学 Speed and current double closed-loop fuzzy control PMSM sensorless control method
CN113965126A (en) * 2021-10-25 2022-01-21 深圳市振邦智能科技股份有限公司 Rotor initial position detection method based on current differential method
CN113965126B (en) * 2021-10-25 2024-02-20 深圳市振邦智能科技股份有限公司 Rotor initial position detection method based on current differentiation method
CN114337402A (en) * 2021-12-08 2022-04-12 燕山大学 Synchronous reluctance motor inductance parameter identification and rotor position hybrid estimation method
CN114337402B (en) * 2021-12-08 2023-09-22 燕山大学 Synchronous reluctance motor inductance parameter identification and rotor position hybrid estimation method
CN114204867A (en) * 2021-12-16 2022-03-18 上海氢恒汽车电子有限公司 Rotor position determination method of permanent magnet synchronous motor
CN115864929A (en) * 2022-12-27 2023-03-28 江西安百川电气有限公司 Method and device for detecting rotor step loss of permanent magnet synchronous motor
CN115864929B (en) * 2022-12-27 2023-06-06 江西安百川电气有限公司 Method and device for detecting rotor out-of-step of permanent magnet synchronous motor
CN116317770A (en) * 2023-02-03 2023-06-23 北京中科昊芯科技有限公司 Method for identifying motor stator resistance offline, method for controlling motor and medium
CN116317770B (en) * 2023-02-03 2024-01-26 北京中科昊芯科技有限公司 Method for identifying motor stator resistance offline, method for controlling motor and medium

Similar Documents

Publication Publication Date Title
CN102291079A (en) Speed sensor-free control algorithm for direct-drive permanent magnet synchronous wind power generation system
CN101964624B (en) Sensorless control system of permanent magnet synchronous motor
CN104009697B (en) Substation inspection robot uses the method for mixing observation device detection positional information
CN106655942B (en) Permanent magnet synchronous motor method for controlling position-less sensor
CN110350835A (en) A kind of permanent magnet synchronous motor method for controlling position-less sensor
CN103346726B (en) Based on the permanent-magnetic synchronous motor stator Flux Observation Method of expansion flux observer
CN103414423A (en) Surface-mounted permanent magnet synchronous motor sensorless direct torque control method
CN103138672A (en) Active disturbance rejection control method of direct-driven permanent magnet synchronization wind power system
CN107579690B (en) A kind of ultrahigh speed permanent magnet synchronous motor method for estimating rotating speed based on sliding formwork observation
CN105450106B (en) Magneto alternator rotor position estimation method
CN103338002A (en) Method for identifying permanent magnet flux and quadrature axis inductance of permanent magnet synchronous motor
Zhou et al. Sensorless direct torque control for electrically excited synchronous motor based on injecting high-frequency ripple current into rotor winding
CN102208894A (en) Rotating speed evaluation method used for speed-sensorless induction motor vector controlling
CN107911057A (en) A kind of rotor-position robust observation procedure for flywheel energy storage system
CN105227010A (en) A kind of permagnetic synchronous motor position-sensor-free position detection error harmonic pulse removing method
Karthikeyan et al. A versatile rotor position computation algorithm for the power control of a grid-connected doubly fed induction generator
CN103904973A (en) Method for controlling salient pole permanent magnet synchronous machine in speed sensorless mode
CN106208872A (en) PMSM high frequency voltage Square wave injection method for estimating rotating speed based on immune algorithm
CN104300861A (en) Method for controlling three-phase permanent magnet synchronous motor
CN102647134A (en) Efficiency optimization control method without angle sensor for permanent magnet synchronous motor
CN108377117A (en) Permanent magnet synchronous motor recombination current control system based on PREDICTIVE CONTROL and method
CN106788049A (en) Speedless sensor moment controlling system and method based on cascade sliding mode observer
CN105071736B (en) A kind of wind turbine permanent magnet synchronous motor is without sensorless rotor position detection method
CN101282102A (en) Control method for megawatt permanent magnetism synchronous directly-driving-wind motor group non-speed sensor
Zhang et al. Sensorless control of salient PMSM with EKF of speed and rotor position

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C02 Deemed withdrawal of patent application after publication (patent law 2001)
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20111221