CN102931906B - Method for asynchronous motor rotor flux linkage observation and rotation speed identification - Google Patents

Method for asynchronous motor rotor flux linkage observation and rotation speed identification Download PDF

Info

Publication number
CN102931906B
CN102931906B CN201210486855.0A CN201210486855A CN102931906B CN 102931906 B CN102931906 B CN 102931906B CN 201210486855 A CN201210486855 A CN 201210486855A CN 102931906 B CN102931906 B CN 102931906B
Authority
CN
China
Prior art keywords
centerdot
psi
beta
alpha
rotor flux
Prior art date
Application number
CN201210486855.0A
Other languages
Chinese (zh)
Other versions
CN102931906A (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
Priority to CN201210417887 priority Critical
Priority to CN201210417887.5 priority
Application filed by 浙江大学 filed Critical 浙江大学
Priority to CN201210486855.0A priority patent/CN102931906B/en
Publication of CN102931906A publication Critical patent/CN102931906A/en
Application granted granted Critical
Publication of CN102931906B publication Critical patent/CN102931906B/en

Links

Abstract

The invention discloses a method for asynchronous motor rotor flux linkage observation and rotation speed identification, which comprises the following steps of: according to stator voltage us and current signals is measured in real time, observing to obtain a rotor flux linkage Psi through a sliding mode observer based on a super-twisting theory and using the rotor flux linkage as a reference value; using a rotor flux linkage Psi~r calculated according to a rotor flux linkage current model as an adjustable value; constructing a model reference self-adaptive system (MRAS) based on the rotor flux linkages; and identifying the asynchronous motor rotation speed omega~r through the self-adaptive law. The rotation speed identification method provided by the invention has high robustness for the variation of stator resistance, and can provide accurate rotor flux linkages and rotation speed values, thereby being especially suitable for an asynchronous motor vector control system.

Description

The method of asynchronous machine rotor flux observation and Speed Identification

Technical field

The present invention discloses a kind of discrimination method of the parameter of electric machine, the particularly very strong rotor flux observer of a kind of robustness, and the Speedless sensor Vector Control System of Induction Motor scheme based on this observer.

Background technology

At present, it is ripe that the control mode of threephase asynchronous machine has become, and vector control and direct torque control can meet most of operating mode demand.Which kind of control mode no matter, rotating speed is all a very important controlled quentity controlled variable.But velocity transducer is installed difficulty in some cases, or sometimes in order to save this part cost, people only wish that, according to the phase current easily recording on frequency converter, phase voltage data, real-time identification rotating speed, realizes Speedless sensor and control.

Speed Sensorless Induction Motor principle can be divided into two types: a generic request rotor has asymmetry, as rotor Harmonic Method, high-frequency signal injection, these class methods need to be analyzed the frequency spectrum of signal, and program wastes time and energy, very harsh to hardware requirement during high speed; Another kind of method is the Mathematical Modeling based on asynchronous machine, with certain mathematical method identification rotating speed wherein, as the method based on state observer, relates to the method for artificial intelligence, and the method based on model reference adaptive principle (MRAS).

Nowadays, researcher has been developed a lot of Speedless sensor algorithms based on asynchronous machine Mathematical Modeling.Belong to having of state observer category: full Order Observers and reduced-order state observer, extended Kalman filter (EKF), and sliding mode observer (SMO).State observer method, to parameter of electric machine sensitive, makes algorithm complicated in order to meet Existence of Global Stable; EKF calculation of complex, a large amount of random parameters will be debugged and obtain; SMO robustness is stronger, but intrinsic shake is harmful to motor low cruise.The method that relates to artificial intelligence is one of study hotspot of the industry always, is just limited to hardware, with a certain distance from practical, also has.The explicit physical meaning of tradition MRAS, algorithm is simpler, and stable state accuracy is relatively good, but is subject to parameter of electric machine variable effect, such as stator resistance R s.Along with the further investigation to MRAS algorithm, it is found that and select different reference model and adjustable model, can develop and the MRAS identification algorithm of different structure, as based on rotor flux, based on back-emf, based on instantaneous reactive power etc.In addition, apply sliding formwork principle, fuzzy control principle etc. under MRAS structure, can develop and a lot of different structures, the leeway of Improvement is very large.

Judging a kind of quality of Speed Identification algorithm, is mainly to see that can this algorithm keep identification precision in a wider speed adjustable range.The Mathematical Modeling of motor is always accurate not, and some parameters also can and change along with motor operation, thereby greatly affects the accuracy of identification.Under low speed, stator resistance R svariation to rotor flux ψ rhaving the greatest impact of observation, and then the identification precision of the Speedless sensor of impact based on flux observer.For R simpact, have scholar in observer, to add R sadaptive Identification improve low-speed performance, but this method has relatively high expectations to the stability Design of observer, versatility is not strong.

Summary of the invention

In order to solve the problems of the technologies described above, the present invention proposes the irrelevant rotor flux observation procedure of a kind of and stator resistance, rotor flux that can online observation asynchronous machine, Speed Identification.The method changes insensitive to stator resistance, and has stronger robustness to disturbing.

The method of this asynchronous machine rotor flux observation and Speed Identification, according to the stator voltage u recording in real time s, stator current i s, by based on the theoretical sliding mode observer of ultimate attainment distortion (Super-Twisting), observe and obtain rotor flux ψ rand as reference value, the rotor flux of being tried to achieve by rotor flux current model and as adjustable value, form the model reference adaptive system (Model Reference AdaptiveSystem, MRAS) based on rotor flux, by adaptive law, carry out identification Rotational Speed of Asynchronous Motor

By the theoretical sliding mode observer of described ultimate attainment distortion (Super-Twisting), observation obtains rotor flux ψ rprocess as follows:

First, threephase asynchronous machine is write to following form in the Mathematical Modeling under rest frame:

i · sα = - k 1 · i sα + k 2 · ( b · ψ rα + ω r · ψ rβ ) + k 3 · u sα i · sβ = - k 1 · i sβ + k 2 · ( b · ψ rβ - ω r · ψ rα ) + k 3 · u sβ ψ · rα = a · i sα - b · ψ rα - ω r · ψ rβ ψ · rβ = a · i sβ - b · ψ rβ + ω r · ψ rα

In above formula, i s αand i s βfor stator current i scomponent, u s αand u s βfor stator voltage u scomponent, ψ r αand ψ r βfor rotor flux ψ rcomponent;

R sfor stator resistance, L mfor the steady-state value of magnetizing inductance, L sand L rbe respectively stator, inductor rotor, T rfor rotor time constant, ω rfor rotor speed, other scales that occur in above formula are shown:

σ = 1 - L m 2 / L s L r

k 1 = R s / σL s + L m 2 / σL s L r T r

k 2=L m/σL sL r

k 3=1/σL s

a=L m/T r

b=1/T r

Secondly, the conversion of above-mentioned model being done based on ultimate attainment distortion theory is as follows:

z 1 = i sα z 2 = i sβ z 3 = b · ψ rα + ω r · ψ rβ z 4 = b · ψ rβ - ω r · ψ rα

Z wherein 1~ z 4for mathematics intermediate variable;

After conversion, the Mathematical Modeling under threephase asynchronous machine rest frame is written as:

z · 1 = - k 1 · z 1 + k 2 · z 3 + k 3 · u sα z · 2 = - k 1 · z 2 + k 2 · z 4 + k 3 · u sβ

So, according to the form of ultimate attainment distortion theory, build rotor flux observer as follows:

z ^ · 1 = - k 1 · z 1 + k 2 · z ~ 3 + k 3 · u sα + λ 1 · | e 1 | 0.5 · sgn ( e 1 ) z ~ · 3 = δ 1 · sgn ( e 1 ) z ^ · 2 = - k 1 · z 2 + k 2 · z ~ 4 + k 3 · u sβ + λ 2 · | e 2 | 0.5 · sgn ( e 2 ) z ~ · 4 = δ 2 · sgn ( e 2 )

Wherein, e 1and e 2for error amount, have sgn () is-symbol function, λ 1, λ 2and δ 1, δ 2for predefined sliding formwork gain;

Finally, the measured value obtaining from above formula with they and rotor flux have following relation:

ψ rα = ∫ 0 t ( - z ~ 3 + L m T r · z ^ 1 ) · dt

ψ rβ = ∫ 0 t ( - z ~ 4 + L m T r · z ^ 2 ) · dt

In above formula, L mfor magnetizing inductance, T rfor rotor time constant.

The rotor flux being obtained by above-mentioned observer is expressed as [ψ under rest frame r αψ r β] t, as the reference value of model reference adaptive method.

Described rotor flux current model, as adjustable model, carrys out identification Rotational Speed of Asynchronous Motor by adaptive law process as follows:

First, build the observer based on rotor flux current model,

ψ ^ rα = 1 1 + T r · p · ( L m · i sα - ω ^ r · T r · ψ ^ rβ )

ψ ^ rβ = 1 1 + T r · p · ( L m · i sβ + ω ^ r · T r · ψ ^ rα )

In above formula, ψ ^ rα ψ ^ rβ T Be the rotor flux value of estimating, p is differential operator, the Rotational Speed of Asynchronous Motor that identification obtains, T rfor rotor time constant;

Secondly, by the output valve [ψ based on the theoretical sliding mode observer of ultimate attainment distortion (Super-Twisting) r αψ r β] tas reference value, ψ ^ rα ψ ^ rβ T As adjustable value, according to model reference adaptive (MRAS) principle, adaptive error value is tried to achieve by following formula:

ϵ = ψ ^ rα · ψ rβ - ψ ^ rβ · ψ rα

Finally, according to the adaptive law of stability principle design speed, with the equivalence of PI controller, so the output valve of PI controller is Identification of Rotational Speed of Asynchronous Motor result

Beneficial effect of the present invention is, has adopted the sliding mode observer based on ultimate attainment distortion theory, and this observer does not need stator resistance R svalue, can in full speed range, guarantee rotor flux ψ raccuracy of observation.Using its output valve as with reference to value, and the MRAS Speedless sensor of structure based on rotor flux, can effectively improve the identification precision of motor speed in low speed situation.

Accompanying drawing explanation

Fig. 1 Vector Control System for Asynchronous Machine schematic diagram;

The structural representation of Fig. 2 asynchronous machine rotor flux observation and Speed Identification method;

The structural representation of the rotor flux sliding mode observer of Fig. 3 based on ultimate attainment distortion theory;

Fig. 4 flux observation and Speed Identification result figure-rotating speed contrast;

Fig. 5 flux observation and Speed Identification result figure-magnetic linkage contrast.

Embodiment

Below in conjunction with drawings and Examples, the invention will be further elaborated.

Referring to Fig. 1, strong power part, three-phase alternating-current supply obtains DC bus-bar voltage U through not controlling rectification dc, service voltage source type inverter, then obtain supplying with the three phase mains of asynchronous machine.

Weak current part, adopt vector control mode, comprise voltage, current sensor, 3 phase/2 phase static coordinate conversion modules, 2 mutually static/2 same pace coordinate transformation modules, Type New Observer module, rotor flux amplitude judge module, speed ring PI controller module, electric current loop PI controller module, 2 same pace/2 phase static coordinate conversion modules, Realization of pulse width modulation based on voltage space vector module.

The present invention relates generally to Type New Observer module, and other modules are the required functional module of Vector Control System of Induction Motor, is general knowledge known in this field.

The workflow of whole system is described below, to introduce the annexation of each module.

1. by transducer, recorded each phase current and the voltage of threephase asynchronous machine, input " 3 phase/2 phase static coordinate conversion module ", obtains stator current i scomponent i s αand i s β, stator voltage u scomponent u s αand u s β.

2. utilize stator voltage, current signal, by Type New Observer of the present invention, obtain real-time rotate speed and synchronous speed angle.Type New Observer comprises: (a) the sliding mode observer module based on ultimate attainment distortion theory, (b) rotor flux space bit angle setting computing module, (c) rotor flux current model block, (d) error amount computing module, (e) rotating speed adaptation module.Concrete details as shown in Figure 2.

(a) first, the Mathematical Modeling under threephase asynchronous machine rest frame is write to following form:

i · sα = - k 1 · i sα + k 2 · ( b · ψ rα + ω r · ψ rβ ) + k 3 · u sα i · sβ = - k 1 · i sβ + k 2 · ( b · ψ rβ - ω r · ψ rα ) + k 3 · u sβ ψ · rα = a · i sα - b · ψ rα - ω r · ψ rβ ψ · rβ = a · i sβ - b · ψ rβ + ω r · ψ rα

In above formula, under rest frame, i s αand i s βfor stator current i scomponent, u s αand u s βfor stator voltage u scomponent, ψ r αand ψ r βfor rotor flux ψ rcomponent;

R sfor stator resistance, L mfor the steady-state value of magnetizing inductance, L sand L rbe respectively stator, inductor rotor, T rfor rotor time constant, ω rfor rotor speed, other scales that occur in above formula are shown:

σ = 1 - L m 2 / L s L r

k 1 = R s / σL s + L m 2 / σL s L r T r

k 2=L m/σL sL r

k 3=1/σL s

a=L m/T r

b=1/T r

Secondly, to above-mentioned model, do based on the theoretical conversion of ultimate attainment distortion (Super-Twisting) as follows:

z 1 = i sα z 2 = i sβ z 3 = b · ψ rα + ω r · ψ rβ z 4 = b · ψ rβ - ω r · ψ rα

Z wherein 1~ z 4for mathematics intermediate variable;

After conversion, the Mathematical Modeling under threephase asynchronous machine rest frame is written as:

z · 1 = - k 1 · z 1 + k 2 · z 3 + k 3 · u sα z · 2 = - k 1 · z 2 + k 2 · z 4 + k 3 · u sβ

So the form theoretical according to ultimate attainment distortion (Super-Twisting) builds rotor flux observer:

z ^ · 1 = - k 1 · z 1 + k 2 · z ~ 3 + k 3 · u sα + λ 1 · | e 1 | 0.5 · sgn ( e 1 ) z ~ · 3 = δ 1 · sgn ( e 1 ) z ^ · 2 = - k 1 · z 2 + k 2 · z ~ 4 + k 3 · u sβ + λ 2 · | e 2 | 0.5 · sgn ( e 2 ) z ~ · 4 = δ 2 · sgn ( e 2 )

Wherein, e 1and e 2for error amount, have sgn () is-symbol function, λ 1, λ 2and δ 1, δ 2for predefined sliding formwork gain.

By above formula discretization, obtain following formula statement:

z ^ 1 ( k ) = z ^ 1 ( k - 1 ) + T · [ - k 1 · z 1 ( k - 1 ) + k 2 · z ~ 3 ( k - 1 ) + k 3 · u sα ( k - 1 ) + λ 1 · | e 1 ( k - 1 ) | 0.5 · sgn ( e 1 ( k - 1 ) ) ] z ~ 3 ( k ) = z ~ 3 ( k - 1 ) + T · δ 1 · sgn ( e 1 ( k - 1 ) ) z ^ 2 ( k ) = z ^ 2 ( k - 1 ) + T · [ - k 1 · z 2 ( k - 1 ) + k 2 · z ~ 4 ( k - 1 ) + k 3 · u sβ ( k - 1 ) + λ 2 · | | e 2 ( k - 1 ) | | 0.5 · sgn ( e 2 ( k - 1 ) ) ] z ~ 4 ( k ) = z ~ 4 ( k - 1 ) + T · δ 2 · sgn ( e 2 ( k - 1 ) )

Wherein T represents system treatment cycle, and k represents certain once-through operation.

Finally, above formula is followed the tracks of stator current component i s αand i s β, output observed result with they and rotor flux have following relation:

ψ rα = ∫ 0 t ( - z ~ 3 + L m T r · z ^ 1 ) · dt

ψ rβ = ∫ 0 t ( - z ~ 4 + L m T r · z ^ 2 ) · dt

Wherein, L mfor magnetizing inductance, T rfor rotor time constant;

In the implementation procedure of above formula, there are two point values to illustrate.

The first, in above formula, comprise pure integral element.Pure integral element is affected by initial value for integral and drift, has direct current biasing and initial phase problem.For addressing this problem, with low-pass first order filter, replace pure integral element, and suitable amplitude and phase compensation in addition.In addition, according to ultimate attainment distortion (Super-Twisting), although the rotor flux observer that theoretical form builds also contains pure integral element, it has applied sliding formwork simultaneously, so the meeting that affects of direct current biasing is corrected by sliding formwork effect.

The second, in above formula, contain rotor time constant T r.When motor moves, T rvalue can be because temperature rise and slip change the kelvin effect unusual fluctuation causing, and change.But, in real system ratio large two orders of magnitude, so T rvariation little on the observed result impact of rotor flux.

During concrete enforcement, use the low-pass first order filter with amplitude and phase compensation, replace pure integral element, its flow chart as shown in Figure 3.Suppose that current stator frequency is ω e, the cut-off frequency ω of low-pass first order filter c, compensating gain COEFFICIENT K should be K=ω c/ ω e.During actual use, can keep K is normal value 3, and cut-off frequency ω cgiven according to motor running condition.When whole system has just brought into operation, by cut-off frequency ω cbe made as a suitable fixed value (as 30rad/s); Treat following space bit angle setting θ 1after stable, ω cby θ 1the stator angular frequency that converts and obtain through differential edetermine.

The rotor flux being obtained by above-mentioned observer is expressed as [ψ under rest frame r αψ r β] t, as the reference value of model reference adaptive method.

(b) by the rotor flux value [ψ of above-mentioned observation r αψ r β] t, calculate its space bit angle setting θ 1.

θ 1 = a tan ( ψ rα ψ rβ )

The foundation that this space bit angle setting is rotor flux linkage orientation, i.e. the space bit angle setting of synchronous speed coordinate system.

(c), using rotor flux current model as adjustable model, by adaptive law, come the process of identification Rotational Speed of Asynchronous Motor as follows.

First, build the observer based on rotor flux current model:

ψ ^ rα = 1 1 + T r · p · ( L m · i sα - ω ^ r · T r · ψ ^ rβ )

ψ ^ rβ = 1 1 + T r · p · ( L m · i sβ + ω ^ r · T r · ψ ^ rα )

In above formula, ψ ^ rα ψ ^ rβ T Be the rotor flux value of estimating, p is differential operator, the Rotational Speed of Asynchronous Motor that identification obtains, L mfor magnetizing inductance, T rfor rotor time constant.

Secondly, by the output valve [ψ based on the theoretical sliding mode observer of ultimate attainment distortion (Super-Twisting) r αψ r β] tas reference value, ψ ^ rα ψ ^ rβ T As adjustable value, according to model reference adaptive (MRAS) principle, under rest frame, adaptive error value is tried to achieve by following formula:

ϵ = ψ ^ rα · ψ rβ - ψ ^ rβ · ψ rα

Finally, according to the adaptive law of stability principle design speed, with the equivalence of PI controller, therefore the output valve of " rotating speed adaptation module " PI controller is Identification of Rotational Speed of Asynchronous Motor result

3. by the amount under rest frame, and the rotor flux space bit angle setting θ being obtained by observer 1, input " 2 mutually static/2 same pace coordinate transformation modules ", obtains the stator current i under synchronous speed coordinate system scomponent i sdand i sq.

4. rotor flux amplitude is controlled: input speed set-point by " rotor flux amplitude judge module ", judge whether motor operates under weak magnetic state.If motor is in permanent magnetic flux running status, output stator electric current d axle component set-point for normal value.If motor weak magnetic field operation, obtains according to look-up table will as set-point, i sdas value of feedback, by " electric current loop PI controller module ", the set-point of output stator voltage d axle component

5. speed and torque are controlled: outside appointment as rotary speed setting value, by identification result as value of feedback, by " speed ring PI controller module ", output stator electric current q axle component set-point will as set-point, i sqas value of feedback, by " electric current loop PI controller module ", the set-point of output stator voltage q axle component

6. inverter control signal produces: by above-mentioned stator voltage set-point of trying to achieve with and rotor flux space bit angle setting θ 1, input " 2 same pace/2 phase static coordinate conversion module ", obtains the component of stator voltage under rest frame with above-mentioned component calculates inverter trigger impulse via Realization of pulse width modulation based on voltage space vector module, is defeated by inverter, can control each phase voltage of supplying with threephase asynchronous machine.

Embodiment 1

Provide a specific embodiment below, the parameter of dispatching from the factory of a threephase asynchronous machine is as following table:

Parameter of electric machine table

Apply the method for asynchronous machine rotor flux observation of the present invention and Speed Identification, controlling motor running condition stops for startup-40r/min-500r/min-1000r/min-500r/min-40r/min-, move altogether approximately 40 seconds, the duration of each operation conditions and result as shown in Figure 4, Figure 5, two different typesettings that figure is same result.

In Fig. 4: A is that rotating speed, the B of algorithm identification of the present invention is that rotating speed, C that photoelectric code disk reads are that poor, the D that two kinds of methods obtain rotating speeds is as the flux observation result (being close to coincidences) with reference to model and adjustable model, and E is the waveform of black surround place amplification in figure.

In Fig. 5: A is that poor, C that Speed Identification is the flux observation of reference model and adjustable model with code-disc rotating speed (be close to and overlap), B are to be as with reference to model flux observation result as adjustable model flux observation result, D, and E is the waveform that in figure, amplify at black surround place.

Above-mentioned Fig. 4, Fig. 5 are exported by control system digital-to-analogue conversion, with the DL750 of YOKOGAWA company oscillograph recording, obtain.

Under 40r/min operation conditions, the rotor flux observer based on ultimate attainment distortion theory, its observed result is more accurate than conventional art, can be as with reference to model, so the present invention's Speed Identification knot is under the low speed comparatively accurate.

Claims (2)

1. a method for asynchronous machine rotor flux observation and Speed Identification, is characterized in that, according to the stator voltage u of the asynchronous machine recording in real time s, stator current i s, by the sliding mode observer based on ultimate attainment distortion theory, observation obtains rotor flux ψ rand as reference value, the rotor flux of being tried to achieve by rotor flux current model as adjustable value, form the model reference adaptive system based on rotor flux, by adaptive law, carry out identification Rotational Speed of Asynchronous Motor
By the described sliding mode observer based on ultimate attainment distortion theory, observation obtains rotor flux ψ rprocess as follows:
First, threephase asynchronous machine is write to following form in the Mathematical Modeling under rest frame:
i · sα = - k 1 · i sα + k 2 · ( b · ψ rα + ω r · ψ rβ ) + k 3 · u sα i · sβ = - k 1 · i sβ + k 2 · ( b · ψ rβ - ω r · ψ rα ) + k 3 · u sβ ψ · rα = a · i sα - b · ψ rα - ω r · ψ rβ ψ · rβ = a · i sβ - b · ψ rβ + ω r · Ψ rα
In above formula, i s αand i s βfor stator current i scomponent, u s αand u s βfor stator voltage u scomponent, ψ r αand ψ r βfor rotor flux ψ rcomponent;
R sfor stator resistance, L mfor the steady-state value of magnetizing inductance, L sand L rbe respectively stator, inductor rotor, T rfor rotor time constant, ω rfor rotor speed, other scales that occur in above formula are shown:
σ = 1 - L m 2 / L s L r
k 1 = R s / σL s + L m 2 / σL s L r T r
k 2=L m/σL sL r
k 3=1/σL s
a=L m/T r
b=1/T r
Secondly, the conversion of above-mentioned model being done based on ultimate attainment distortion theory is as follows:
z 1 = i sα z 2 = i sβ z 3 = b · ψ rα + ω r · ψ rβ z 4 = b · ψ rβ - ω r · ψ rα
Z wherein 1~z 4for mathematics intermediate variable;
After conversion, the Mathematical Modeling under threephase asynchronous machine rest frame is written as:
z · 1 = - k 1 · z 1 + k 2 · z 3 + k 3 · u sα z · 2 = - k 1 · z 2 + k 2 · z 4 + k 3 · u sβ
So, according to the form of ultimate attainment distortion theory, build rotor flux observer:
z ^ · 1 = - k 1 · z 1 + k 2 · z ~ 3 + k 3 · u sα + λ 1 · | e 1 | 0.5 · sgn ( e 1 ) z ~ · 3 = δ 1 · sgn ( e 1 ) z ^ · 2 = - k 1 · z 2 + k 2 · z ~ 4 + k 3 · u sβ + λ 2 · | e 2 | 0.5 · sgn ( e 2 ) z ~ · 4 = δ 2 · sgn ( e 2 )
Wherein, e 1and e 2for error amount, have sgn () is-symbol function, λ 1, λ 2and δ 1, δ 2for predefined sliding formwork gain;
Finally, the measured value obtaining from above formula with they and rotor flux have following relation:
ψ rα = ∫ 0 t ( - z ~ 3 + L m T r · z ^ 1 ) · dt
ψ rβ = ∫ 0 t ( - z ~ 4 + L m T r · z ^ 2 ) · dt
In above formula, L mfor magnetizing inductance, T rfor rotor time constant;
The rotor flux being obtained by above-mentioned observer is expressed as [ψ under rest frame r αψ r β] t, as the reference value of model reference adaptive method.
2. the method for rotor flux observation according to claim 1 and Speed Identification, is characterized in that, described rotor flux current model is as adjustable model, by adaptive law identification Rotational Speed of Asynchronous Motor process as follows:
First, build the observer based on rotor flux current model,
ψ ^ rα = 1 1 + T r · p · ( L m · i sα - ω ^ r · T r · ψ ^ rβ )
ψ ^ rβ = 1 1 + T r · p · ( L m · i sβ + ω ^ r · T r · ψ ^ rα )
In above formula, be the rotor flux value of observation, p is differential operator, the Rotational Speed of Asynchronous Motor that identification obtains, T rfor rotor time constant;
Secondly, by the output valve [ψ of the sliding mode observer based on ultimate attainment distortion theory r αψ r β] tas reference value, as adjustable value, according to model reference adaptive method, adaptive error value is tried to achieve by following formula:
ϵ = ψ ^ rα · ψ rβ - ψ ^ rβ · ψ rα
Finally, according to the adaptive law of stability principle design speed, with the equivalence of PI controller, so the output valve of PI controller is Identification of Rotational Speed of Asynchronous Motor result
CN201210486855.0A 2012-10-26 2012-11-23 Method for asynchronous motor rotor flux linkage observation and rotation speed identification CN102931906B (en)

Priority Applications (3)

Application Number Priority Date Filing Date Title
CN201210417887 2012-10-26
CN201210417887.5 2012-10-26
CN201210486855.0A CN102931906B (en) 2012-10-26 2012-11-23 Method for asynchronous motor rotor flux linkage observation and rotation speed identification

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210486855.0A CN102931906B (en) 2012-10-26 2012-11-23 Method for asynchronous motor rotor flux linkage observation and rotation speed identification

Publications (2)

Publication Number Publication Date
CN102931906A CN102931906A (en) 2013-02-13
CN102931906B true CN102931906B (en) 2014-11-12

Family

ID=47646636

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210486855.0A CN102931906B (en) 2012-10-26 2012-11-23 Method for asynchronous motor rotor flux linkage observation and rotation speed identification

Country Status (1)

Country Link
CN (1) CN102931906B (en)

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103715962B (en) * 2013-12-25 2016-10-05 西安理工大学 The permagnetic synchronous motor sliding-mode speed observer that dual stage matrix converter drives
CN104009696B (en) * 2014-05-08 2017-02-08 昆明理工大学 Interactive model reference adaptive speed and stator resistance identification method based on sliding-mode control
CN104092378A (en) * 2014-06-16 2014-10-08 西北工业大学 Control method of robust high order sliding mode for Flyback convertor
CN104065320B (en) * 2014-06-17 2017-08-01 合肥学院 Double-fed wind power generator group Speed Sensorless Control Method
CN105207561B (en) * 2015-09-08 2018-06-01 湖南工程学院 A kind of Flux Observation Method and equipment for motor
CN105634367A (en) * 2016-02-23 2016-06-01 江苏方天电力技术有限公司 MRAS-based high-voltage asynchronous motor control method
CN105610369B (en) * 2016-03-04 2018-05-29 合肥工业大学 A kind of asynchronous machine Flux Observation Method based on sliding mode observer
CN105717453B (en) * 2016-03-25 2018-10-23 珠海格力节能环保制冷技术研究中心有限公司 A kind of detection method of motor magnetic linkage and the detection method of electrical fault
CN105932924B (en) * 2016-06-07 2018-06-29 合肥工业大学 Asynchronous machine rotor angular rate observation procedure
CN106788051B (en) * 2016-11-23 2019-02-05 江苏大学 A kind of induction-type bearingless motor method for estimating rotating speed
CN106849799B (en) * 2017-03-20 2019-03-05 江苏大学 A kind of method of induction-type bearingless motor flux linkage observation and Speed Identification
CN107370436A (en) * 2017-06-27 2017-11-21 上海电机学院 A kind of direct Torque Control based on Super twisting sliding formworks

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1889358A (en) * 2006-07-11 2007-01-03 合肥工业大学 Permasyn electric machine control system based on adaptive sliding moding structure

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1889358A (en) * 2006-07-11 2007-01-03 合肥工业大学 Permasyn electric machine control system based on adaptive sliding moding structure

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
"交流电机无速度传感器矢量控制系统变结构模型参考自适应转速辨识";王庆龙等;《中国电机工程学报》;20070531;第27卷(第15期);第70-74页 *
S. SOLVAR等."Sensorless second order sliding mode observer for induction motor".《2010 IEEE International Conference on Control Applications,Part of 2010 IEEE Multi-Conference on Systems and Control》.2010,1933-1938. *
S. SOLVAR等."Sensorless second order sliding mode observer for induction motor".《2010 IEEE International Conference on Control Applications,Part of 2010 IEEE Multi-Conference on Systems and Control》.2010,1933-1938. *
王庆龙等."交流电机无速度传感器矢量控制系统变结构模型参考自适应转速辨识".《中国电机工程学报》.2007,第27卷(第15期),第70-74页. *

Also Published As

Publication number Publication date
CN102931906A (en) 2013-02-13

Similar Documents

Publication Publication Date Title
Benlaloui et al. Implementation of a new MRAS speed sensorless vector control of induction machine
CN103647490B (en) A kind of sliding mode control strategy of magneto
CN103401501B (en) A kind of PMSM servo system control method based on fuzzy active disturbance rejection
Yin et al. Research on robust performance of speed-sensorless vector control for the induction motor using an interfacing multiple-model extended Kalman filter
Rehman et al. A new current model flux observer for wide speed range sensorless control of an induction machine
CN103414427B (en) Brushless direct current motor control method
CN103066910B (en) Direct torque control method capable of hindering torque impulse
CN103427746B (en) A kind of face mounted permagnetic synchronous motor full speed range method for controlling position-less sensor
Lascu et al. A class of speed-sensorless sliding-mode observers for high-performance induction motor drives
US6014006A (en) Induction motor control system with speed and flux estimation
AU2013306382B2 (en) System and method for error correction in angular position sensors
KR101244299B1 (en) Method and system for sensorless torque control for ironless permanent magnet machine
CN100583620C (en) Rotor position presuming method and apparatus, motor control method and compressor
CN103828221B (en) The method and system of motor is controlled with variable switching frequency under variable running speed
CN100477483C (en) Motor controller
CN100530935C (en) Vector controller for a permanent magnet synchronous motor
CN101388639B (en) Non-position sensor vector control method for double-feed wind power generator
JP5492192B2 (en) AC motor control device
CN104335476B (en) The control device of motor and the control method of motor
Xu et al. Very-low speed control of PMSM based on EKF estimation with closed loop optimized parameters
CN102420561B (en) Speed sensorless vector control method on basis of cascaded high voltage inverter
CN1964182B (en) Electrical equipment
CN103412484B (en) A kind of control moment gyro framework disturbing moment suppressing method
CN101247104B (en) Conjecture method of stator magnetic linkage oriented AC motor rotation velocity and rotor position
Marchesoni et al. A simple approach to flux and speed observation in induction motor drives

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant