Control method of permanent magnet synchronous motor

文档序号:1523761 发布日期:2020-02-11 浏览:12次 中文

阅读说明:本技术 永磁同步电机的控制方法 (Control method of permanent magnet synchronous motor ) 是由 张瑞峰 柴璐军 王晓妮 杨高兴 詹哲军 梁海刚 贺志学 于 2019-10-29 设计创作,主要内容包括:本发明涉及交流电机的控制方法,具体为永磁同步电机的控制方法。解决现有永磁同步电机的控制方法中因使用的电机参数不准确造成电机输出的转矩精度和电机运行效率被影响的问题。该控制方法通过电机定子温度T,使用在线查表得到定子电阻<I>R</I><Sub><I>s</I></Sub><I>(T)</I>;通过电机定子温度T、电流幅值<I>I</I><Sub><I>s</I></Sub>、电流矢量角<I>β</I>,使用在线查表得到定子电感<I>L</I><Sub><I>d</I></Sub><I>(T)</I>和<I>L</I><Sub><I>q</I></Sub><I>(T)</I>;使用磁链观测模型实时计算出磁链值<I>Ψ</I><Sub><I>f</I></Sub><I>(T</I><Sub><I>r</I></Sub><I>)</I>,从而提高了电机控制和解耦的准确性;利用转矩闭环的输出结果重新分配了给定的定子电流,维持了永磁同步电机按照较优的控制轨迹运行,降低了电机的发热和损耗。实现了电机的解耦控制、增强了电机控制系统的鲁棒性。(The invention relates to a control method of an alternating current motor, in particular to a control method of a permanent magnet synchronous motor. The problem that the torque precision of motor output and the motor operation efficiency are influenced due to inaccurate motor parameters in the existing control method of the permanent magnet synchronous motor is solved. The control method obtains the stator resistance through the temperature T of the motor stator by using an online table look-up R s (T) (ii) a Through the temperature T and the current amplitude of the stator of the motor I s Current vector angle β Obtaining stator inductance using online table lookup L d (T) And L q (T) (ii) a Real-time calculation of flux linkage values using flux linkage observation models Ψ f (T r ) Therefore, the accuracy of motor control and decoupling is improved; the given stator current is redistributed by utilizing the output result of the torque closed loop, the permanent magnet synchronous motor is maintained to run according to a better control track, and the heating and the loss of the motor are reduced. The decoupling control of the motor is realized, and the robustness of a motor control system is enhanced.)

1. The control method of the permanent magnet synchronous motor is characterized in that a control block diagram of the control method comprises a rotary transformer module (1), a temperature sensor module (2), a stator resistance calculation module (3), a Clark conversion module (4), a Park conversion module (5), a stator inductance calculation table look-up module (6), a permanent magnet flux linkage calculation module (7), a torque calculation module (8), a current calculation module (9), a given current generation module (10), a robust decoupling controller module (11), a PWM modulation module (12) and a three-phase inverter bridge module (13);

1) rotary transformer module

The rotary transformer is arranged on the permanent magnet synchronous motor, the rotor position theta of the permanent magnet synchronous motor is obtained through measurement of the rotary transformer, and the rotating speed w of the permanent magnet synchronous motor is obtained through differentiation of the rotor position theta e

2) Temperature sensor module

A temperature sensor is embedded in the motor stator, and the real-time stator temperature T of the motor is obtained by the temperature sensor;

3) stator resistance calculation module

Motor stator resistance R caused by motor temperature change sTesting and drawing a temperature resistance comparison table of the motor stator, acquiring a real-time temperature value T of the stator through a temperature sensor, and acquiring a stator resistor R by inquiring the temperature resistance comparison table of the motor stator s(T);

4) Clark conversion module

Collecting two-phase stator current i a、i bObtaining stator current i through Clark conversion α、i β

5) Park conversion module

Stator current i α、i βObtaining the current i under the d-q static coordinate system through Park conversion d、i q

6) Stator inductance calculation table look-up module

Stator inductance L dAnd L qWith the motor stator temperature T and the amplitude I of the stator current SAnd phase β of stator current SThree variables are related;

firstly, the stator inductance L is manufactured at different temperature points within a certain temperature interval d、L qWith I SAnd β SA table of changes;

then using the real-time collected stator temperature T and the amplitude I of the stator current SAnd phase β of stator current SObtaining the stator inductance L by looking up the table d(T)、L q(T);

7) Permanent magnet flux linkage calculation module

The input of the permanent magnet flux linkage calculation module is as follows: current i d、i qStator resistance R s(T), stator inductance L d(T) and L q(T), stator voltage u 'of the last calculation cycle' d,u’ qPermanent magnet flux linkage psi calculated in the previous calculation cycle f’(T r) Rotational speed w e(ii) a Output Ψ of the permanent magnet flux linkage calculation Module f(T r) The calculation method is as follows:

Figure RE-FDA0002312158640000021

wherein: k p_ΨfAnd K i_ΨfIs a control parameter that is a function of,

Figure RE-FDA0002312158640000022

debugging to obtain K p_ΨfAnd K i_ΨfThe process of (2) is as follows: given K i_ΨfA smaller parameter (e.g., 0.001), the parameter K is adjusted first p_ΨfPsi making the output derived f(T r) In a state of constant amplitude oscillation, then the parameter K is adjusted i_ΨfTo make the output Ψ f(T r) Convergence, parameter K obtained at this time p_ΨfAnd K i_ΨfIs a control parameter;

Figure RE-FDA0002312158640000023

Figure RE-FDA0002312158640000025

8) torque calculation module

The input to the torque calculation module is the current i d、i qStator resistance R s(T), stator inductance L d(T) and L q(T), rotational speed w e,Ψ f(T r);

Output of torque computing module-electromagnetic torque T ebThe following calculation formula is obtained:

T ex1=1.5n pi q[(L d(T)-L q(T))i d]

in the formula, n pIs the number of pole pairs of the motor;

T ex2=1.5n pi q[(L d(T 0)-L q(T 0))i d]

in the formula, L d(T 0) And L q(T 0) The values of the stator inductances under rated operating conditions are respectively;

T ex=T ex1-T ex2

T est=1.5n pi qf(T r)+(L d(T 0)-L q(T 0))i d]

T eb=T ex+T est

9) current calculation module

Target torque T * eAfter amplitude limiting and ramp processing, the given torque T is obtained * e1

T * e1、T ebStator inductance L d(T), stator inductance L q(T)、Ψ f(T r) And n pIs the input of the current calculation module;

given torque T * e1Obtaining a given current value I through a calculation module *And current angle β;

T * e1and T ebThe difference of (a) is output as delta β through a PI regulator;

the final current angle β is obtained 1Calculated by the following formula:

β 1=β+Δβ

10) given current generation module

The input of the given current generation module is a given current value I *And current angle β 1

The calculation process is as follows:

Figure RE-FDA0002312158640000031

11) robust decoupling controller

The input parameter of the robust decoupling controller is current i d *、i q *、i d、i qStator resistor R s(T), stator inductance L d(T) and L q(T), permanent magnet flux linkage Ψ f(T r) And a rotational speed w eThe output parameter is the stator voltage u d、u q

Intermediate variable u of control algorithm d1、u q1The calculation process is as follows:

Figure RE-FDA0002312158640000032

in the formula, β xIs a control parameter, which is expressed as follows:

β x=β b×f kx

in the formula, β bIs a reference value of a control parameter, selected between 0.1 and 1.0, f kxIs the switching frequency of the inverter;

intermediate variable u of control algorithm d2、u q2The calculation process is as follows:

Figure RE-FDA0002312158640000033

in the formula, delta d、δ qThe calculation formula of (a) is as follows:

Figure RE-FDA0002312158640000034

in the formula, lambda is a control parameter, and the control parameter lambda is selected by adopting a trial and error method; i.e. i d1、i q1As intermediate variables in the algorithmic process, i d1、i q1The calculation formula of (a) is as follows:

Figure RE-FDA0002312158640000041

wherein, delta' d、δ’ qIs the variable calculated in the last calculation period (or called the last beat);

voltage u d、u qThe calculation process of (2) is as follows:

Figure RE-FDA0002312158640000042

K Ris a damping coefficient;

12) PWM modulation module

PWM modulation moduleIs a stator voltage u d、u qDc bus voltage u dcRotational speed w eAnd an angle θ; the output of the PWM modulation module is 6 paths of PWM waves, and the three-phase inverter bridge module is driven to work.

2. The control method of the permanent magnet synchronous motor according to claim 1, wherein in 6) the stator inductance calculation look-up table module:

considering i qAt 0, the program is computationally problematic, so that the denominator is given a particularly small number k α

3. The control method of the permanent magnet synchronous motor according to claim 2, wherein in 6) the stator inductance calculation look-up table module: amplitude I of stator current SAnd phase β of stator current SRespectively low-pass filtering to obtain filtered values I SLPFAnd β SLPF

4. The control method of the permanent magnet synchronous motor according to claim 2, wherein in 6) the stator inductance calculation look-up table module: the temperature interval is the interval of [ -30 ℃, 160 ℃ ] and the temperature values of each integral multiple of ten are taken as the temperature points, so that in the interval of [ -30 ℃, 160 ℃ ], twenty temperature points are present: -30 ℃, 20 ℃, 10 ℃, 0 ℃, 10 ℃, 20 ℃, … …, 150 ℃ and 160 ℃.

5. The control method of the permanent magnet synchronous motor according to claim 3, wherein in 6) the stator inductance calculation look-up table module: the temperature interval is the interval of [ -30 ℃, 160 ℃ ] and the temperature values of each integral multiple of ten are taken as the temperature points, so that in the interval of [ -30 ℃, 160 ℃ ], twenty temperature points are present: -30 ℃, 20 ℃, 10 ℃, 0 ℃, 10 ℃, 20 ℃, … …, 150 ℃ and 160 ℃.

6. The control method of the permanent magnet synchronous motor according to claim 4, wherein the process of making the table is: the tested permanent magnet synchronous motor and the counter-dragging motor are coaxial or connected through a gear box, the counter-dragging motor operates below a rated rotating speed, and different I of the permanent magnet synchronous motor is given under the condition that the rotating speed is stable and the temperature of a stator is at a certain temperature point S,I S>0 and β S,90°<β S<180 DEG, measured by high frequency injection method, of different I SAnd β SCorresponding motor stator inductance L d、L q(ii) a Changing the temperature of the stator to obtain a plurality of groups of stator inductances L at different temperature points d、L qWith I SAnd β SA table of changes.

7. The control method of a permanent magnet synchronous motor according to claim 5, wherein the process of making a table is: the tested permanent magnet synchronous motor and the counter-dragging motor are coaxial or connected through a gear box, the counter-dragging motor operates below a rated rotating speed, and different I of the permanent magnet synchronous motor is given under the condition that the rotating speed is stable and the temperature of a stator is at a certain temperature point SLPF,I SLPF>0 and β SLPF,90°<β SLPF<180 DEG, measured by high frequency injection method, of different I SLPFAnd β SLPFCorresponding motor stator inductance L d、L q(ii) a Changing the temperature of the stator to obtain a plurality of groups of stator inductances L at different temperature points d、L qWith I SLPFAnd β SLPFA table of changes.

8. The control method of the permanent magnet synchronous motor according to claim 6, wherein in 6) the stator inductance calculation look-up table module: the table lookup procedure is as follows: when the stator temperature T acquired in real time is not equal to the temperature value of any temperature point, two temperature points T closest to the stator temperature T are selected x、T x+10Stator inductance L d、L qWith I SAnd β SUsing the changed table as a look-up tableA table and a second table; when the stator temperature T acquired in real time is equal to the temperature value of a certain temperature point, selecting the temperature point T equal to the stator temperature T xStator inductance L d、L qWith I SAnd β STable of changes, temperature point T adjacent to temperature increasing side of the temperature point x+10Stator inductance L d、L qWith I SAnd β SA changed table, a first table and a second table used as table lookup; according to the amplitude I of the collected stator current SAnd phase β of stator current SLooking up L in the first table d1(T) and L q1(T), finding L in the second table d2(T) and L q2(T); l is obtained as follows d(T) and L q(T):

Figure RE-FDA0002312158640000051

9. The control method of the permanent magnet synchronous motor according to claim 7, wherein in 6) the stator inductance calculation look-up table module: the table lookup procedure is as follows: when the stator temperature T acquired in real time is not equal to the temperature value of any temperature point, two temperature points T closest to the stator temperature T are selected x、T x+10Stator inductance L d、L qWith I SLPFAnd β SLPFThe changed tables are used as a first table and a second table for table lookup; when the stator temperature T acquired in real time is equal to the temperature value of a certain temperature point, selecting the temperature point T equal to the stator temperature T xStator inductance L d、L qWith I SLPFAnd β SLPFTable of changes, temperature point T adjacent to temperature increasing side of the temperature point x+10Stator inductance L d、L qWith I SLPFAnd β SLPFA changed table, a first table and a second table used as table lookup; according to the collected I SLPFAnd β SLPFLooking up L in the first table d1(T) and L q1(T), finding L in the second table d2(T) and L q2(T); according to the followingThe following formula gives L d(T) and L q(T):

10. Method for controlling a permanent magnet synchronous machine according to any of claims 1-9, characterized in that, 9) in the current calculation module, a torque T is given * e1Obtaining a given current value I through a calculation module *And current angle β are calculated as follows:

obtaining a per unit value base value t through operation ebAnd i bxWherein i bxIs the per unit value base value of the current, passing through i bx=ψ f(T r)/(L q(T)-L d(T)) calculating; t is t ebIs a per unit value base value of torque, which can be expressed by t eb=n pψ f(T r)i bxCalculating to obtain;

given torque

Figure RE-FDA0002312158640000062

in per unit value, the relationship between torque and current of the control algorithm is expressed as:

Figure RE-FDA00023121586400000610

A given current value I *By the formula

Figure RE-FDA00023121586400000619

current angle β passing through

Figure RE-FDA00023121586400000620

in view of

Figure RE-FDA0002312158640000071

Technical Field

The invention relates to a control method of a motor, in particular to a control method of a permanent magnet synchronous motor.

Background

The permanent magnet synchronous motor has the advantages of high efficiency, energy conservation, high power density and the like, and is widely applied to the field of rail transit. For a permanent magnet synchronous motor, the most important performance requirement is to generate accurate torque and realize high-efficiency control under various environments, and one of the most important factors influencing the torque precision and the high-efficiency control of the permanent magnet motor is the change of motor parameters in the motor operation process, which causes the mismatching of a motor control algorithm and a controlled motor. The change of motor parameters is mainly caused by the magnetic saturation effect of an iron core caused by the change of the working temperature of the motor and the stator current, thereby causing the stator inductance L in the motor parameters d、L qStator resistor R sAnd a permanent magnet flux linkage Ψ fA change in (c). The traditional vector control method is widely used in permanent magnet synchronous motor control, but does not consider temperature to motor parameter R s、Ψ fAnd iron core saturation effect pair L d、L qThe accuracy of the torque output from the motor and the efficiency of the motor operation will be greatly affected.

Disclosure of Invention

The invention solves the problem that the torque precision of the motor output and the motor operation efficiency are influenced due to inaccurate motor parameters in the existing control method of the permanent magnet synchronous motor, and provides the control method of the permanent magnet synchronous motor. Meanwhile, on the basis of accurately determining the motor parameters, a robust decoupling control method for the permanent magnet synchronous motor is provided, decoupling control of the motor is achieved, and robustness of a motor control system is enhanced.

The invention is realized by adopting the following technical scheme: the control block diagram of the optimization control method of the permanent magnet synchronous motor comprises a rotary transformer module, a temperature sensor module, a stator resistance calculation module, a Clark conversion module, a Park conversion module, a stator inductance calculation table look-up module, a permanent magnet flux linkage calculation module, a torque calculation module, a current calculation module, a given current generation module, a robust decoupling controller module, a PWM (pulse width modulation) module and a three-phase inverter bridge module;

1) rotary transformer module

The rotary transformer is arranged on the permanent magnet synchronous motor, the rotor position theta of the permanent magnet synchronous motor is obtained through measurement of the rotary transformer, and the rotating speed w of the permanent magnet synchronous motor is obtained through differentiation of the rotor position theta e

2) Temperature sensor module

A temperature sensor is embedded in the motor stator, and the real-time stator temperature T of the motor is obtained by the temperature sensor;

3) stator resistance calculation module

Motor stator resistance R caused by motor temperature change sTesting and drawing a temperature resistance comparison table of the motor stator, acquiring a real-time temperature value T of the stator through a temperature sensor, and acquiring a stator resistor R by inquiring the temperature resistance comparison table of the motor stator s(T);

4) Clark conversion module

Collecting two-phase stator current i a、i bObtaining stator current i through Clark conversion α、i β

5) Park conversion module

Stator current i α、i βObtaining the current i under the d-q static coordinate system through Park conversion d、i q

6) Stator inductance calculation table look-up module

StatorInductor L dAnd L qWith the motor stator temperature T and the amplitude I of the stator current SAnd phase β of stator current SThree variables are related;

firstly, the stator inductance L is manufactured at different temperature points within a certain temperature interval d、L qWith I SAnd β SA table of changes;

then using the real-time collected stator temperature T and the amplitude I of the stator current SAnd phase β of stator current SObtaining the stator inductance L by looking up the table d(T)、L q(T);

7) Permanent magnet flux linkage calculation module

The input of the permanent magnet flux linkage calculation module is as follows: current i d、i qStator resistance R s(T), stator inductance L d(T) and L q(T), stator voltage u 'of the last calculation cycle (or title to last beat)' d,u’ qThe permanent magnet flux linkage Ψ calculated in the previous calculation cycle (or referred to as the previous beat) f’(T r) Rotational speed w e. Output Ψ of the permanent magnet flux linkage calculation Module f(T r) The calculation method is as follows:

Figure BDA0002252455400000021

wherein: k p_ΨfAnd K i_ΨfIs a control parameter that is a function of, respectively representing the estimated current values, T being the acquisition time, T rRepresents the actual operating temperature of the permanent magnet;

debugging to obtain K p_ΨfAnd K i_ΨfThe process of (2) is as follows: given K i_ΨfA smaller parameter (e.g., 0.001), the parameter K is adjusted first p_ΨfPsi making the output derived f(T r) In a state of constant amplitude oscillation, then the parameter K is adjusted i_ΨfTo make the output Ψ f(T r) Convergence, parameter K obtained at this time p_ΨfAnd K i_ΨfIs a control parameter.

Figure BDA0002252455400000023

And

Figure BDA0002252455400000024

calculated by the following formula:

Figure BDA0002252455400000031

8) torque calculation module

The input to the torque calculation module is the current i d、i qStator resistance R s(T), stator inductance L d(T) and L q(T), rotational speed w e,Ψ f(T r);

Output of torque computing module-electromagnetic torque T ebThe following calculation formula is obtained:

T ex1=1.5n pi q[(L d(T)-L q(T))i d]

in the formula, n pIs the pole pair number of the motor.

T ex2=1.5n pi q[(L d(T 0)-L q(T 0))i d]

In the formula, L d(T 0) And L q(T 0) The values of the stator inductances in the nominal operating mode are in each case.

T ex=T ex1-T ex2

T est=1.5n pi qf(T r)+(L d(T 0)-L q(T 0))i d]

T eb=T ex+T est

9) Current calculation module

Target torque T * eAfter amplitude limiting and ramp processing, the given torque T is obtained * e1

T * e1、T ebStator inductance L d(T), stator inductance L q(T)、Ψ f(T r) And n pIs the input of the current calculation module;

given torque T * e1Obtaining a given current value I through a calculation module *And current angle β (shown in FIG. 3);

given torque T * e1Obtaining a given current value I through a calculation module *And current angle β are calculated as follows:

obtaining a per unit value base value t through operation ebAnd i bxWherein i bxIs the per unit value base value of the current, passing through i bx=ψ f(T r)/(L q(T)-L d(T)) calculating; t is t ebIs a per unit value base value of torque, which can be expressed by t eb=n pψ f(T r)i bxCalculating to obtain; the per-unit value base value varies with the variation of the motor parameter.

Given torque

Figure BDA0002252455400000032

And its per unit value t enPassing through type

Figure BDA0002252455400000033

Represents; electric current

Figure BDA0002252455400000034

And its per unit value i dnCan pass through

Figure BDA0002252455400000035

Represents; electric current

Figure BDA0002252455400000036

And its per unit value i qnCan pass through

Figure BDA0002252455400000037

Represents the current of

Figure BDA0002252455400000038

And current

Figure BDA0002252455400000041

Is to calculate I *An intermediate variable of (d);

in per unit value, the relationship between torque and current of the control algorithm is expressed as:

Figure BDA0002252455400000042

by commanding a given torque

Figure BDA0002252455400000043

Becomes per unit value t enIn a format of, and then by a formula

Figure BDA0002252455400000044

Solving to obtain current

Figure BDA0002252455400000045

Per unit value of i dnAnd finally, the formula of

Figure BDA0002252455400000046

Given current can be calculated Electric current

Figure BDA0002252455400000048

Per unit value of i qnCan be represented by the formula t en=i qn(1-i dn) Is calculated to be t at this time enAnd i dnIs a known amount, and then a formula

Figure BDA0002252455400000049

Given current can be calculated

Figure BDA00022524554000000410

A given current value I *By the formula

Figure BDA00022524554000000411

Calculating to obtain;

current angle β passing through Calculating to obtain;

in view of

Figure BDA00022524554000000413

At 0, the program is computationally problematic, so that the denominator is given a particularly small number k b,k bMay be equal to 0.0000001;

T * e1and T ebThe difference is output to be delta β through a PI regulator, delta β is a compensation value of a given current angle, and the given current angle is corrected by identifying real-time torque and comparing the real-time torque with the given torque;

the final current angle β is obtained 1Calculated by the following formula:

β 1=β+Δβ

10) given current generation module

The input of the given current generation module is a given current value I *And current angle β 1

The calculation process is as follows:

11) robust decoupling controller

The input parameter of the robust decoupling controller is current i d *、i q *、i d、i qStator resistor R s(T), stator inductance L d(T) and L q(T), permanent magnet flux linkage Ψ f(T r) And a rotational speed w eThe output parameter is the stator voltage u d、u q

Intermediate variable u of control algorithm d1、u q1The calculation process is as follows:

Figure BDA0002252455400000051

in the formula, β xIs a control parameter, control parameter β xVarying with modulation strategy, control parameter β xProportional to the switching frequency, which is expressed as follows:

β x=β b×f kx

in the formula, β bIs a reference value of a control parameter, selected between 0.1 and 1.0, f kxIs the switching frequency of the inverter;

u is calculated by the above control algorithm d1、u q1The decoupling control of the permanent magnet synchronous motor is realized;

intermediate variable u of control algorithm d2、u q2The calculation process is as follows:

Figure BDA0002252455400000052

in the formula, delta d、δ qIs a key part of a robust controller, and the calculation formula is as follows:

Figure BDA0002252455400000053

in the formula, lambda is a control parameter, and the control parameter lambda is selected by a trial and error method, and can be 80; i.e. i d1、i q1As intermediate variables in the algorithmic process, i d1、i q1The calculation formula of (a) is as follows:

Figure BDA0002252455400000054

wherein, delta' d、δ’ qIs the variable calculated in the last calculation period (or called the last beat);

voltage u d、u qThe calculation process of (2) is as follows:

Figure BDA0002252455400000055

in order to enhance the stability of the control system and realize control decoupling, the term R is respectively added to the dq axis S(T)×k R×i dAnd R S(T)×k R×i qf(T r)w e,K RThe damping coefficient is 0.9;

12) PWM modulation module

The input of the PWM modulation module is a stator voltage u d、u qDc bus voltage u dcRotational speed w eAnd an angle θ; the output of the PWM modulation module is 6 paths of PWM waves, and the three-phase inverter bridge module is driven to work.

The invention has the following beneficial effects:

(1) the invention passes the temperature T and the current amplitude I of the stator of the motor sThe current vector angle β is looked up on line to obtain the accurate stator resistance R of the motor parameter s(T), stator inductance L d(T) and L q(T) calculating the flux linkage value Ψ in real time using a flux linkage observation model f(T r) Detection equipment for the temperature of the rotor is omitted, and the accuracy of motor control and decoupling is improved; and the given stator current is redistributed by utilizing the output result of the torque closed loop, so that the permanent magnet synchronous motor is kept to run in a better track, and the heating and the loss of the motor are reduced.

(2) By using a robust decoupling control method in the current loop, the anti-interference performance of system control is improved, and decoupling control of the permanent magnet synchronous motor is realized.

Drawings

FIG. 1 is a control block diagram of the control method of the present invention;

FIG. 2 is a block diagram of a flow of a stator inductance calculation look-up table module;

FIG. 3 is a control block diagram of a current calculation module;

FIG. 4 is a control block diagram of a robust decoupling controller;

fig. 5 is a schematic diagram of a segment modulation algorithm.

Detailed Description

The control block diagram (as shown in fig. 1) of the control method of the permanent magnet synchronous motor comprises a rotary transformer module 1, a temperature sensor module 2, a stator resistance calculation module 3, a Clark conversion module 4, a Park conversion module 5, a stator inductance calculation look-up table module 6, a permanent magnet flux linkage calculation module 7, a torque calculation module 8, a current calculation module 9, a given current generation module 10, a robust decoupling controller module 11, a PWM modulation module 12 and a three-phase inverter bridge module 13;

1) rotary transformer module

The rotary transformer is arranged on the permanent magnet synchronous motor, the rotor position theta of the permanent magnet synchronous motor is obtained through measurement of the rotary transformer, and the rotating speed w of the permanent magnet synchronous motor is obtained through differentiation of the rotor position theta e

2) Temperature sensor module

A temperature sensor is embedded in the motor stator, and the real-time stator temperature T of the motor is obtained by the temperature sensor;

3) stator resistance calculation module

Motor stator resistance R caused by motor temperature change sTesting and drawing a temperature resistance comparison table of the motor stator, acquiring a real-time temperature value T of the stator through a temperature sensor, and acquiring a stator resistor R by inquiring the temperature resistance comparison table of the motor stator s(T);

4) Clark conversion module

Collecting two-phase stator current i a、i bObtaining stator current i through Clark conversion α、i β

5) Park conversion module

Stator current i α、i βObtaining the current i under the d-q static coordinate system through Park conversion d、i q

6) Stator inductance calculation table look-up module

The change of the stator current can cause the magnetic saturation effect of the stator core, and the stator inductance L is changed along with the change of the d-axis current and the q-axis current d、L qWill change and the stator temperature T of the motor will also have an influence on the stator inductance. To getTo a more accurate stator inductance parameter L dAnd L qObtaining stator inductance L by using table look-up method dAnd L q

Stator inductance L dAnd L qWith the motor stator temperature T and the amplitude I of the stator current SAnd phase β of stator current SThree variables are related (as shown in fig. 2);

wherein:

Figure BDA0002252455400000071

considering i qAt 0, the program is computationally problematic, so that the denominator is given a particularly small number k α,k αMay be equal to 0.0000001;

further, the amplitude I of the stator current SAnd phase β of stator current SRespectively low-pass filtering to obtain filtered values I SLPFAnd β SLPF

Firstly, the stator inductance L is manufactured at different temperature points within a certain temperature interval d、L qWith I SAnd β S(or I) SLPFAnd β SLPF) The changed tables (the number of the tables is consistent with the number of the temperature points set in the temperature interval, and the table corresponding to each temperature point reflects the stator inductance L d、L qWith I SAnd β S(or I) SLPFAnd β SLPF) Content of the change); in specific implementation, the temperature range is-30 ℃ and 160 DEG C]Interval, and taking the temperature value of integral multiple of each ten as the temperature point, so that the temperature is between minus 30 ℃ and 160 DEG C]Interval, there are twenty temperature points: -30 ℃, -20 ℃, -10 ℃, 0 ℃, 10 ℃, 20 ℃, … …, 150 ℃, 160 ℃; the specific process of making the table is as follows: the tested permanent magnet synchronous motor and the counter-dragging motor are coaxial or connected through a gear box, the counter-dragging motor operates below the rated rotating speed, and different current instructions I of the permanent magnet synchronous motor are given under the conditions that the rotating speed is stable and the temperature of a stator is at a certain temperature point (taking-20 ℃ as an example) S,I S>0 and β S,90°<β S<180°(or I) SLPF,I SLPF>0 and β SLPF,90°<β SLPF<180 degree (Current i) * d、i * q) By high frequency injection of different I SAnd β S(or I) SLPFAnd β SLPF) Corresponding motor stator inductance L d、L q(ii) a Changing the temperature of the stator to obtain a plurality of groups of stator inductances L at different temperature points d、L qWith I SAnd β S(or I) SLPFAnd β SLPF) A table of changes.

Then using the real-time collected stator temperature T and the amplitude I of the stator current SAnd phase β of stator current S(or I) SLPFAnd β SLPF) Obtaining real-time stator inductance L by querying the table d(T)、L q(T); the table lookup procedure is as follows: when the stator temperature T acquired in real time is not equal to the temperature value of any temperature point, two temperature points T closest to the stator temperature T are selected x、T x+10Stator inductance L d、L qWith I SAnd β S(or I) SLPFAnd β SLPF) The first table and the second table used as the table look-up table (e.g., stator temperature T of 23 deg.C, selected temperature point T) x20 ℃ and T x+10Table at 30 ℃ as first and second tables); when the stator temperature T acquired in real time is equal to the temperature value of a certain temperature point, selecting the temperature point T equal to the stator temperature T xStator inductance L d、L qWith I SAnd β S(or I) SLPFAnd β SLPF) Table of changes, temperature point T adjacent to temperature increasing side of the temperature point x+10Stator inductance L d、L qWith I SAnd β S(or I) SLPFAnd β SLPF) A changed table, a first table and a second table used as table lookup; according to the amplitude I of the collected stator current SAnd phase β of stator current S(or I) SLPFAnd β SLPF) Looking up L in the first table d1(T) and L q1(T), finding L in the second table d2(T) and L q2(T);L is obtained as follows d(T) and L q(T):

Figure BDA0002252455400000081

L d(T) and L q(T) is the current operating condition (current temperature T, current I) SAnd β S(or I) SLPFAnd β SLPF) Motor stator inductance value.

7) Permanent magnet flux linkage calculation module

The permanent magnet material in the permanent magnet synchronous motor rotor is greatly influenced by temperature change, and the change relation between the permanent magnet flux linkage of the motor and the temperature can be expressed as

Figure BDA0002252455400000082

Wherein: Ψ f(T 0) Is the permanent magnet flux linkage psi under the rated working condition f(T r) Is the permanent magnet flux linkage, T, at the actual operating temperature of the motor 0May be taken at 20 ℃ T rIs the actual operating temperature of the permanent magnet and α is the temperature coefficient of remanence.

Because the permanent magnet of the motor is arranged on the rotor of the motor, the actual working temperature of the motor cannot be obtained, and the flux linkage value psi is obtained in real time by adopting a flux linkage observation model f(T r) The input of the permanent magnet flux linkage calculation module is as follows: current i d、i qStator resistance R s(T), stator inductance L d(T) and L q(T), stator voltage u 'of the last calculation cycle (or title to last beat)' d,u’ qThe permanent magnet flux linkage Ψ calculated in the previous calculation cycle (or referred to as the previous beat) f’(T r) Rotational speed w eOutput Ψ of the permanent magnet flux linkage calculation Module f(T r) The calculation method is as follows:

wherein: k p_ΨfAnd K i_ΨfIs a control parameter that is a function of, respectively representing the estimated current values, T being the acquisition time, T rRepresents the actual operating temperature of the permanent magnet;

debugging to obtain K p_ΨfAnd K i_ΨfThe process of (2) is as follows: given K i_ΨfA smaller parameter (e.g., 0.001), the parameter K is adjusted first p_ΨfPsi making the output derived f(T r) In a state of constant amplitude oscillation, then the parameter K is adjusted i_ΨfTo make the output Ψ f(T r) Convergence, parameter K obtained at this time p_ΨfAnd K i_ΨfIs a control parameter;

and

Figure BDA0002252455400000094

calculated by the following formula:

Figure BDA0002252455400000095

8) torque calculation module

The input to the torque calculation module is the current i d、i qStator resistance R s(T), stator inductance L d(T) and L q(T), rotational speed w e,Ψ f(T r);

Output of torque computing module-electromagnetic torque T ebThe following calculation formula is obtained:

T ex1=1.5n pi q[(L d(T)-L q(T))i d]

in the formula, n pIs the number of pole pairs of the motor;

T ex2=1.5n pi q[(L d(T 0)-L q(T 0))i d]

in the formula, L d(T 0) And L q(T 0) The values of the stator inductances under rated operating conditions are respectively;

T ex=T ex1-T ex2

T est=1.5n pi qf(T r)+(L d(T 0)-L q(T 0))i d]

T eb=T ex+T est

9) current calculation module

Target torque T * eAfter amplitude limiting and ramp processing, the given torque T is obtained * e1

T * e1、T ebStator inductance L d(T) and L q(T)、Ψ f(T r) And n pIs the input of the current calculation module;

given torque T * e1Obtaining a given current value I through a calculation module *And current angle β (shown in FIG. 3);

given torque T * e1Obtaining a given current value I through a calculation module *And current angle β are calculated as follows:

obtaining a per unit value base value t through operation ebAnd i bxWherein i bxIs the per unit value base value of the current, passing through i bx=ψ f(T r)/(L q(T)-L d(T)) calculating; t is t ebIs a per unit value base value of torque, which can be expressed by t eb=n pψ f(T r)i bxCalculating to obtain; the per-unit value base value varies with the variation of the motor parameter.

Given torque

Figure BDA0002252455400000101

And its per unit value t enCan pass through

Figure BDA0002252455400000102

Represents; electric current And its per unit value i dnCan pass through

Figure BDA0002252455400000104

Represents; electric current

Figure BDA0002252455400000105

And its per unit value i qnCan pass through

Figure BDA0002252455400000106

Represents the current of And current

Figure BDA0002252455400000108

Is to calculate I *An intermediate variable of (d);

in per unit value, the relationship between torque and current of the control algorithm is expressed as:

Figure BDA0002252455400000109

by applying a given torque

Figure BDA00022524554000001010

Becomes per unit value t enIn a format of, and then by a formula

Figure BDA00022524554000001011

Solving to obtain current

Figure BDA00022524554000001012

Per unit value of i dnAnd finally, the formula of

Figure BDA00022524554000001013

Given current can be calculated

Figure BDA00022524554000001014

Electric current

Figure BDA00022524554000001015

Per unit value of i qnCan be represented by the formula t en=i qn(1-i dn) Is calculated to be t at this time enAnd i dnIs a known amount, and then a formula

Figure BDA00022524554000001016

Given current can be calculated

Figure BDA00022524554000001017

A given current value I *By the formula Calculating to obtain;

current angle β passing through

Figure BDA00022524554000001019

Calculating to obtain;

in view of At 0, the program is computationally problematic, so that the denominator is given a particularly small number k b,k bMay be equal to 0.0000001;

T * e1and T ebThe difference is output to be delta β through a PI regulator, delta β is a compensation value of a given current angle, and the given current angle is corrected by identifying real-time torque and comparing the real-time torque with the given torque;

the final current angle β is obtained 1Calculated by the following formula:

β 1=β+Δβ

10) given current generation module

The input of the given current generation module is a given current value I *And current angle β 1

The calculation process is as follows:

Figure BDA0002252455400000111

11) robust decoupling controller

The robust decoupling controller realizes decoupling control of the permanent magnet synchronous motor on one hand and improves the anti-interference function of motor control on the other hand, and the input parameter of the robust decoupling controller is current i d *、i q *、i d、i qStator resistor R s(T), stator inductance L d(T) and L q(T), permanent magnet flux linkage Ψ f(T r) And a rotational speed w eWith stator voltage u as an output parameter d、u q

A robust decoupling controller is shown in fig. 4.

Intermediate variable u of control algorithm d1、u q1The calculation process is as follows:

Figure BDA0002252455400000112

in the formula, β xIs a control parameter, control parameter β xVarying with modulation strategy, control parameter β xProportional to the switching frequency, which is expressed as follows:

β x=β b×f kx

in the formula, β bIs a control parameter reference value selected from 0.1-1.0 (e.g., 0.1, 0.2, 0.3, 0.4, 0.5, 0.7, 0.9, 1.0; preferably 0.5), f kxIs the switching frequency of the inverter;

u is calculated by the above control algorithm d1、u q1The decoupling control of the permanent magnet synchronous motor is realized;

intermediate variable u of control algorithm d2、u q2The calculation process is as follows:

Figure BDA0002252455400000113

in the formula, delta d、δ qIs a key part of a robust controllerThe calculation formula is as follows:

in the formula, lambda is a control parameter, and the control parameter lambda is selected by a trial and error method, and can be 80. i.e. i d1、i q1As intermediate variables in the algorithmic process, i d1、i q1The calculation formula of (a) is as follows:

Figure BDA0002252455400000121

wherein, delta' d、δ’ qIs the variable calculated in the last calculation period (or called the last beat);

voltage u d、u qThe calculation process of (2) is as follows:

Figure BDA0002252455400000122

in order to enhance the stability of the control system and realize control decoupling, the term R is respectively added to the dq axis S(T)×k R×i dAnd R S(T)×k R×i qf(T r)w e,K RThe damping coefficient can be 0.9;

12) PWM modulation module

The input of the PWM modulation module is a stator voltage u d、u qDc bus voltage u dcRotational speed w eAnd an angle theta. The output of the PWM modulation module is 6 paths of PWM waves, and the three-phase inverter bridge module is driven to work.

Due to the restriction of conditions such as heat dissipation, the modulation algorithm of the high-power permanent magnet synchronous motor is limited by the switching frequency, and a segmented modulation strategy combining multiple modulation modes is generally adopted. A schematic diagram of the segmented modulation strategy is shown in fig. 5.

The segmented modulation is divided into asynchronous modulation and synchronous modulation, the synchronous modulation can be divided into multiple segments under the limitation of switching frequency, and finally square wave modulation is carried out. Under square wave modulation, the voltage utilization rate is high and harmonic waves are small. In the synchronous segmented modulation algorithm, available modulation algorithms include SPWM modulation, specific subharmonic elimination PWM (SHEPWM) modulation and the like, and various modulation algorithms have advantages, disadvantages and application ranges.

The control method provided by the invention can be realized in a DSP chip, the control of the permanent magnet synchronous motor is realized by two-stage interruption, a motor control algorithm (namely, steps 1-11) is operated in one-stage interruption, and the design control period is 250 us; and operating the PWM algorithm of the motor in the other stage of interruption, wherein the interruption period is related to the operation frequency of the motor according to the current modulation strategy.

18页详细技术资料下载
上一篇:一种医用注射器针头装配设备
下一篇:一种混合定子材料驱动电机的控制系统

网友询问留言

已有0条留言

还没有人留言评论。精彩留言会获得点赞!

精彩留言,会给你点赞!