CN112583315B - Three-vector model prediction torque control method for three-level permanent magnet synchronous motor - Google Patents

Three-vector model prediction torque control method for three-level permanent magnet synchronous motor Download PDF

Info

Publication number
CN112583315B
CN112583315B CN202011352040.4A CN202011352040A CN112583315B CN 112583315 B CN112583315 B CN 112583315B CN 202011352040 A CN202011352040 A CN 202011352040A CN 112583315 B CN112583315 B CN 112583315B
Authority
CN
China
Prior art keywords
vector
voltage vector
optimal voltage
optimal
action
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.)
Active
Application number
CN202011352040.4A
Other languages
Chinese (zh)
Other versions
CN112583315A (en
Inventor
涂建军
李众
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Jiangsu University of Science and Technology
Original Assignee
Jiangsu University of Science and Technology
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 Jiangsu University of Science and Technology filed Critical Jiangsu University of Science and Technology
Priority to CN202011352040.4A priority Critical patent/CN112583315B/en
Publication of CN112583315A publication Critical patent/CN112583315A/en
Application granted granted Critical
Publication of CN112583315B publication Critical patent/CN112583315B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P21/00Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P21/00Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
    • H02P21/14Estimation or adaptation of machine parameters, e.g. flux, current or voltage
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P21/00Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
    • H02P21/14Estimation or adaptation of machine parameters, e.g. flux, current or voltage
    • H02P21/20Estimation of torque
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P6/00Arrangements for controlling synchronous motors or other dynamo-electric motors using electronic commutation dependent on the rotor position; Electronic commutators therefor
    • H02P6/34Modelling or simulation for control purposes
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T10/00Road transport of goods or passengers
    • Y02T10/60Other road transportation technologies with climate change mitigation effect
    • Y02T10/64Electric machine technologies in electromobility

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Control Of Ac Motors In General (AREA)

Abstract

本发明公开了一种三电平永磁同步电动机三矢量模型预测转矩控制方法,采用矢量分区的方法选择电压矢量,通过将三电平逆变器备选矢量中的大矢量代入到价值函数中寻优,从而确定第一最优电压矢量,通过第一最优电压矢量的位置确定扇区,第二最优电压矢量从第一最优电压矢量所在扇区中的三个非零矢量代入价值函数中寻优得到,结合第一最优电压矢量、第二最优电压矢量和零矢量的作用时间控制三电平逆变器动作,这样第一最优电压矢量从原来的寻优27次降低到了6次,这样的寻优过程既考虑到了两个最优电压矢量的选择,又结合了其作用时间,保证所选的两个最优电压矢量能够保证中点电位平衡和较低的开关频率。

Figure 202011352040

The invention discloses a three-level permanent magnet synchronous motor three-vector model prediction torque control method. The vector partition method is adopted to select the voltage vector, and the large vector in the three-level inverter candidate vector is substituted into the value function. The first optimal voltage vector is determined by searching for the optimal voltage vector, the sector is determined by the position of the first optimal voltage vector, and the second optimal voltage vector is substituted from the three non-zero vectors in the sector where the first optimal voltage vector is located. The value function is optimized, and the action time of the first optimal voltage vector, the second optimal voltage vector and the zero vector is combined to control the action of the three-level inverter, so that the first optimal voltage vector is optimized from the original 27 times. It is reduced to 6 times. This optimization process not only considers the selection of the two optimal voltage vectors, but also combines its action time to ensure that the selected two optimal voltage vectors can ensure the balance of the mid-point potential and lower switching. frequency.

Figure 202011352040

Description

一种三电平永磁同步电动机三矢量模型预测转矩控制方法A three-level permanent magnet synchronous motor three-vector model predictive torque control method

技术领域technical field

本发明涉及永磁同步电动机转矩控制,具体是涉及一种三电平永磁同步电动机三矢量模型预测转矩控制方法。The invention relates to torque control of permanent magnet synchronous motors, in particular to a three-level permanent magnet synchronous motor three-vector model predictive torque control method.

背景技术Background technique

由于永磁同步电机结构简单、体积小、功率密度高、效率高、运行性能优良等优点,近年来获得了广泛的应用,例如:交通领域、极爱单领域、车船的驱动领域,以及作为泵类、压缩机、高精度伺服系统等方面的转动机构。因此对永磁同步电机控制的研究也越来越多。Due to the advantages of simple structure, small size, high power density, high efficiency, and excellent operating performance, permanent magnet synchronous motors have been widely used in recent years, such as in the field of transportation, the field of love, the driving field of vehicles and ships, and as a pump Class, compressor, high-precision servo system and other aspects of the rotating mechanism. Therefore, there are more and more studies on the control of permanent magnet synchronous motors.

现有的永磁同步电机控制多用模型预测转矩控制,模型预测转矩控制是有限集预测算法中的一种经典控制方法,它非常适用于永磁同步电机控制系统,具有控制策略简单,动态性能好的优点,更适用于多标量、多约束条件、非线性系统。传统的模型预测转矩控制通过得到系统的预测模型,遍历可知数目的状态矢量,再根据预期控制目标设定合适的价值函数,通过价值函数来比较,即可得到功率器件未来需要执行的开关序列。Existing permanent magnet synchronous motor control uses model predictive torque control. Model predictive torque control is a classic control method in finite set prediction algorithms. It is very suitable for permanent magnet synchronous motor control systems. It has simple control strategies and dynamic The advantages of good performance are more suitable for multi-scalar, multi-constraint, and nonlinear systems. The traditional model predictive torque control obtains the predictive model of the system, traverses a known number of state vectors, and then sets an appropriate value function according to the expected control target. By comparing the value function, the switching sequence that the power device needs to execute in the future can be obtained. .

三电平逆变器相较于两电平逆变器来说有更小的开关频率和更广的矢量选择范围,然而在实际运用中需要考虑到中点电位平衡的问题。现有的模型预测控制策略中将中点电位平衡项引入到价值函数中,在寻优过程中考虑中点电位平衡问题,但也增加了控制策略的实行难度。Compared with the two-level inverter, the three-level inverter has a smaller switching frequency and a wider range of vector selection. However, in practical applications, the balance of the mid-point potential needs to be considered. In the existing model predictive control strategy, the midpoint potential balance term is introduced into the value function, and the midpoint potential balance problem is considered in the optimization process, but it also increases the difficulty of implementing the control strategy.

现有技术缺点:Disadvantages of the existing technology:

(1)已有的模型预测控制策略里在每个控制周期里有且只能得到一个大小固定,方向固定的电压矢量,但实际需要跟踪的参考电压矢量是不定的,因此输出的电压矢量不能有效的使得反馈转矩追踪参考转矩,直接导致永磁同步电机磁通脉动大,输出转矩脉动不平稳,电机稳态性能不佳。(1) In the existing model predictive control strategy, only one voltage vector with a fixed size and a fixed direction can be obtained in each control cycle, but the actual reference voltage vector to be tracked is uncertain, so the output voltage vector cannot be Effectively makes the feedback torque track the reference torque, which directly leads to the large magnetic flux pulsation of the permanent magnet synchronous motor, the unstable output torque pulsation, and the poor steady-state performance of the motor.

(2)三电平或更多电平拓扑结构的器件开关状态数目过多,计算成本过大,且存在复杂的冗余矢量,不易确定开关序列组合;(2) The number of switch states of devices with three-level or more-level topology is too large, the calculation cost is too large, and there are complex redundant vectors, so it is difficult to determine the switch sequence combination;

(3)价值函数中含有中点电位控制项,权重因子的整定变得更加的复杂,导致控制的精度不高,计算复杂;(3) The value function contains the mid-point potential control term, and the setting of the weight factor becomes more complicated, resulting in low control accuracy and complicated calculation;

(4)多矢量预测控制造成较高的开关频率,增加了开关损耗。(4) Multi-vector predictive control causes higher switching frequency and increases switching loss.

发明内容SUMMARY OF THE INVENTION

发明目的:针对以上缺点,本发明提供一种快速筛选矢量、降低开关损耗的三电平永磁同步电动机三矢量模型预测转矩控制方法。Purpose of the invention: In view of the above shortcomings, the present invention provides a three-vector model predictive torque control method for a three-level permanent magnet synchronous motor that can quickly screen vectors and reduce switching losses.

技术方案:为解决上述问题,本发明采用一种三电平永磁同步电动机三矢量模型预测转矩控制方法,包括以下步骤:Technical solution: In order to solve the above problems, the present invention adopts a three-level permanent magnet synchronous motor three-vector model predictive torque control method, which includes the following steps:

(1)根据永磁同步电动的状态空间模型,构建电流和磁通的状态方程,根据电流和磁通的状态方程得到转矩和磁通的预测模型,根据三电平逆变器直流侧电容电压模型得到第一电容电压预测模型;根据第一电容电压预测模型构建第一价值函数;(1) According to the state space model of permanent magnet synchronous motor, the state equations of current and magnetic flux are constructed, and the prediction model of torque and magnetic flux is obtained according to the state equations of current and magnetic flux. According to the DC side capacitance of the three-level inverter obtaining a first capacitor voltage prediction model from the voltage model; constructing a first value function according to the first capacitor voltage prediction model;

(2)采样得到k时刻的电流;(2) Sampling to obtain the current at time k;

(3)将采样取得的数据代入到转矩和磁通的预测模型中,得到下一时刻k+1时刻的转矩和磁链;(3) Substitute the data obtained by sampling into the prediction model of torque and magnetic flux, and obtain the torque and flux linkage at the next moment k+1;

(4)将三电平逆变器空间矢量图中的大矢量作为第一最优电压矢量的备选矢量代入到第一价值函数中,选出第一最优电压矢量;(4) Substitute the large vector in the three-level inverter space vector diagram into the first value function as the candidate vector of the first optimal voltage vector, and select the first optimal voltage vector;

(5)将第一最优电压矢量对应的扇区中的三个非零矢量作为第二最优电压矢量的备选矢量,零矢量作为第三电压矢量;(5) Three non-zero vectors in the sector corresponding to the first optimal voltage vector are used as candidate vectors for the second optimal voltage vector, and the zero vector is used as the third voltage vector;

(6)依次选取三个备选矢量作为第二最优电压矢量,计算第一最优电压矢量、第二最优电压矢量和第三电压矢量的作用时间,并求出占空比,根据占空比调制出新电压矢量;(6) Select three candidate vectors in turn as the second optimal voltage vector, calculate the action time of the first optimal voltage vector, the second optimal voltage vector and the third voltage vector, and obtain the duty ratio. The empty ratio modulates a new voltage vector;

(7)根据第一最优电压矢量、第二最优电压矢量和第三电压矢量的作用时间建立第二电容电压预测模型,并根据第二电容电压预测模型构建第二价值函数;(7) establishing a second capacitance voltage prediction model according to the action time of the first optimal voltage vector, the second optimal voltage vector and the third voltage vector, and building a second cost function according to the second capacitance voltage prediction model;

(8)将步骤6中,调制出的三个新电压矢量代入到第二价值函数中选出最优的第二最优电压矢量;(8) Substitute the three new voltage vectors modulated in step 6 into the second value function to select the optimal second optimal voltage vector;

(9)将所选的第一最优电压矢量、第二最优电压矢量和第三电压矢量对应的开关序列及对应的作用时间输入脉冲调制器控制三电平逆变器动作。(9) Input the switching sequence corresponding to the selected first optimal voltage vector, the second optimal voltage vector and the third voltage vector and the corresponding action time into the pulse modulator to control the action of the three-level inverter.

进一步的,所述步骤1中电流和磁通的状态方程为:Further, the state equations of the current and the magnetic flux in the step 1 are:

Figure BDA0002801587600000021
Figure BDA0002801587600000021

Figure BDA0002801587600000022
Figure BDA0002801587600000022

其中,iα、iβ分别表示定子电流矢量在α、β轴的分量;Ls表示定子电感;uα、uβ分别表示定子电压矢量在α、β轴下的分量;Rs表示定子电阻;ωr表示电动机的转子电角速度;ψα、ψβ分别表示电子磁链矢量在α、β轴下的分量。Among them, i α and i β represent the components of the stator current vector on the α and β axes respectively; L s represents the stator inductance; u α and u β represent the components of the stator voltage vector on the α and β axes respectively; R s represents the stator resistance ; ω r represents the rotor electrical angular velocity of the motor; ψ α , ψ β represent the components of the electron flux linkage vector under the α and β axes, respectively.

进一步的,所述步骤1中,转矩和磁通的预测模型为:Further, in the step 1, the prediction model of torque and magnetic flux is:

Figure BDA0002801587600000031
Figure BDA0002801587600000031

ψs(k+1)=ψs(k)+Ts[us(k)-Rsis(k)]ψ s (k+1)=ψ s (k)+T s [u s (k)-R s i s (k)]

Figure BDA0002801587600000032
Figure BDA0002801587600000032

其中,Ts为采样周期;is(k+1)、ψs(k+1)、Te(k+1)分别为k+1时刻定子电流矢量、定子磁链和转矩矢量;is(k)、ψs(k)、us(k)、ωr(k)分别为k时刻定子电流矢量、定子磁链、定子电压矢量和转子的电角速度;j为虚数单位,p为极对数,

Figure BDA0002801587600000033
为叉乘符号。Among them, T s is the sampling period; i s (k+1), ψ s (k+1), and Te (k+1) are the stator current vector, stator flux linkage and torque vector at time k+1, respectively; i s (k), ψ s (k), u s (k), and ω r (k) are the stator current vector, stator flux linkage, stator voltage vector and the electrical angular velocity of the rotor at time k, respectively; j is an imaginary unit, and p is Number of pole pairs,
Figure BDA0002801587600000033
is the cross product symbol.

进一步的,所述步骤1中,电容电压预测模型为:Further, in the step 1, the capacitance voltage prediction model is:

Figure BDA0002801587600000034
Figure BDA0002801587600000034

其中,vc1(k)、vc2(k)为k时刻采样的电容电压;vc1(k+1)、vc2(k+1)为k+1时刻预测的电容电压;ic1(k)和ic2(k)是根据k时刻逆变器的开关状态和输出电流的定义值,C1、C2为直流侧两电容。Among them, v c1 (k) and v c2 (k) are the capacitor voltages sampled at time k; v c1 (k+1) and v c2 (k+1) are the predicted capacitor voltages at time k+1; i c1 (k ) and i c2 (k) are defined values according to the switching state and output current of the inverter at time k, C 1 and C 2 are the two capacitors on the DC side.

进一步的,所述步骤4中,第一价值函数为:Further, in the step 4, the first value function is:

Figure BDA0002801587600000035
Figure BDA0002801587600000035

其中,

Figure BDA0002801587600000036
为转矩给定值,
Figure BDA0002801587600000037
为定子磁通给定值,λ1、λ2、λ3为权重系数,Fsw为一个控制周期内开关管发生开通和关断的动作变换次数。in,
Figure BDA0002801587600000036
is the torque given value,
Figure BDA0002801587600000037
is the given value of the stator magnetic flux, λ 1 , λ 2 , λ 3 are the weight coefficients, and F sw is the switching times of switching on and off in a control cycle.

进一步的,所述步骤6中,第一最优电压矢量、第二最优电压矢量和第三电压矢量的作用时间分别为t1、t2和t0Further, in the step 6, the action times of the first optimal voltage vector, the second optimal voltage vector and the third voltage vector are respectively t 1 , t 2 and t 0 :

t0=Ts-t1-t2 t 0 =T s -t 1 -t 2

Figure BDA0002801587600000038
Figure BDA0002801587600000038

其中,

Figure BDA0002801587600000039
S为一个量,且
Figure BDA0002801587600000041
Figure BDA0002801587600000042
为转矩和磁通在第一最优电压矢量作用下的变化率;
Figure BDA0002801587600000043
Figure BDA0002801587600000044
为转矩和磁通在第二最优电压矢量作用下的变化率;
Figure BDA0002801587600000045
为转矩和磁通在零矢量作用下的变化率。in,
Figure BDA0002801587600000039
S is a quantity, and
Figure BDA0002801587600000041
Figure BDA0002801587600000042
is the rate of change of torque and magnetic flux under the action of the first optimal voltage vector;
Figure BDA0002801587600000043
Figure BDA0002801587600000044
is the rate of change of torque and magnetic flux under the action of the second optimal voltage vector;
Figure BDA0002801587600000045
is the rate of change of torque and flux under the action of zero vector.

进一步的,t1、t2、t0在0~Ts范围内,且为保证t1、t2、t0在0~Ts范围内,需作如下处理:Further, t 1 , t 2 , and t 0 are in the range of 0 to T s , and to ensure that t 1 , t 2 , and t 0 are in the range of 0 to T s , the following processing is required:

(6.1)当作用时间t1、t2和t0中任意两个作用时间小于0,且另一个作用时间大于Ts,则令小于0的作用时间等于0,大于Ts的作用时间等于Ts(6.1) When any two of the action times t 1 , t 2 and t 0 are less than 0, and the other action time is greater than T s , then the action time less than 0 is equal to 0, and the action time greater than T s is equal to T s ;

(6.2)当作用时间t1+t2≤Ts,则各矢量的作用时间不变;(6.2) When the action time t 1 +t 2 ≤T s , the action time of each vector remains unchanged;

(6.3)当作用时间t1+t2≥Ts,则需要将三个矢量的作用时间代入如下公式进行修正:(6.3) When the action time t 1 +t 2 ≥T s , the action time of the three vectors needs to be substituted into the following formula for correction:

Figure BDA0002801587600000046
Figure BDA0002801587600000046

其中:t1'、t'2和t'0分别表示修正后的第一最优电压矢量、第二最优电压矢量和第三电压矢量的作用时间。Wherein: t 1 ', t' 2 and t' 0 respectively represent the action time of the corrected first optimal voltage vector, the second optimal voltage vector and the third voltage vector.

进一步的,所述步骤6中,新电压矢量的调制公式:Further, in the step 6, the modulation formula of the new voltage vector is:

Figure BDA0002801587600000047
Figure BDA0002801587600000047

其中,uα2、uβ2为新电压矢量的α、β轴分量,uopt1_α、uopt1_β为第一最优电压矢量的α、β轴分量,uopt2_n_α、uopt2_n_β为第二最优电压矢量集合中第n个电压矢量的α、β轴分量,第一最优电压矢量uopt1、第二最优电压矢量uopt2_n与零矢量u0的占空比分别ηopt1、ηopt2_n、η0Among them, u α2 and u β2 are the α and β axis components of the new voltage vector, u opt1_α and u opt1_β are the α and β axis components of the first optimal voltage vector, and u opt2_n_α and u opt2_n_β are the second optimal voltage vector set In the α and β axis components of the nth voltage vector, the duty ratios of the first optimal voltage vector u opt1 , the second optimal voltage vector u opt2_n and the zero vector u 0 are η opt1 , η opt2_n , and η 0 , respectively.

进一步的,所述步骤7中第二电容电压预测模型为:Further, the second capacitor voltage prediction model in the step 7 is:

Figure BDA0002801587600000048
Figure BDA0002801587600000048

其中,v'c1(k+1)、v'c2(k+1)为第一最优电压矢量、第二最优电压矢量和第三电压矢量的作用时间下k+1时刻的直流侧两电容电压,v'c1(k)、v'c2(k)为第一最优电压矢量、第二最优电压矢量和第三电压矢量的作用时间下k时刻两电容电压,ic11(k)、ic21(k)、ic31(k)、ic12(k)、ic22(k)、ic32(k)分别为第一最优电压矢量、第二最优电压矢量和第三电压矢量作用时的输出电流值。Among them, v' c1 (k+1), v' c2 (k+1) are the two voltages on the DC side at time k+1 under the action time of the first optimal voltage vector, the second optimal voltage vector and the third voltage vector Capacitor voltage, v' c1 (k), v' c2 (k) are the two capacitor voltages at time k under the action time of the first optimal voltage vector, the second optimal voltage vector and the third voltage vector, i c11 (k) , ic21 (k), ic31 (k), ic12 (k), ic22 (k), ic32 (k) are the first optimal voltage vector, the second optimal voltage vector and the third voltage vector respectively Output current value when active.

进一步的,所述步骤7中第二价值函数gn为:Further, in the step 7, the second value function g n is:

Figure BDA0002801587600000051
Figure BDA0002801587600000051

其中,gn为第二价值函数。Among them, g n is the second value function.

有益效果:本发明相对于现有技术,其显著优点是通过第一最优电压矢量的位置确定扇区从而选出第二最优电压矢量,降低了寻优次数,并结合两个最优矢量和零矢量的作用时间,保证了中点电位平衡和较低的开关频率;使用三矢量模型预测控制,合成的矢量不仅在方向可以变化,还可以改变其大小,使其达到的控制效果最优。Beneficial effect: Compared with the prior art, the present invention has the significant advantage that the sector is determined by the position of the first optimal voltage vector to select the second optimal voltage vector, which reduces the number of optimizations, and combines the two optimal voltage vectors and the action time of the zero vector to ensure the neutral point potential balance and lower switching frequency; using the three-vector model predictive control, the synthesized vector can not only change in direction, but also change its size to achieve the best control effect. .

附图说明Description of drawings

图1是本发明NPC三电平逆变器拓扑图;1 is a topology diagram of an NPC three-level inverter of the present invention;

图2是本发明提供的矢量分区的三电平基本空间矢量分布图;Fig. 2 is the three-level basic space vector distribution diagram of the vector partition provided by the present invention;

图3是本发明三电平逆变器驱动永磁同步电机三矢量模型预测控制的框图;Fig. 3 is the block diagram of the three-vector model predictive control of the three-level inverter-driven permanent magnet synchronous motor of the present invention;

图4是本发明矢量选择的流程图。Figure 4 is a flow chart of the vector selection of the present invention.

具体实施方式Detailed ways

本实施例中一种三电平永磁同步电动机三矢量模型预测转矩控制方法,其中In this embodiment, a three-level permanent magnet synchronous motor three-vector model predicted torque control method, wherein

根据永磁同步电动的状态空间模型,构建电流和磁通的状态方程为:According to the state space model of permanent magnet synchronous electromotive force, the state equations of current and magnetic flux are constructed as:

Figure BDA0002801587600000052
Figure BDA0002801587600000052

Figure BDA0002801587600000053
Figure BDA0002801587600000053

其中,iα、iβ分别表示定子电流矢量在α、β轴的分量;Ls表示定子电感;uα、uβ分别表示定子电压矢量在α、β轴下的分量;Rs表示定子电阻;ωr表示电动机的转子电角速度;ψα、ψβ分别表示电子磁链矢量在α、β轴下的分量。Among them, i α and i β represent the components of the stator current vector on the α and β axes respectively; L s represents the stator inductance; u α and u β represent the components of the stator voltage vector on the α and β axes respectively; R s represents the stator resistance ; ω r represents the rotor electrical angular velocity of the motor; ψ α , ψ β represent the components of the electron flux linkage vector under the α and β axes, respectively.

根据永磁同步电动机电流和磁通的状态方程,对其中的微分项进行离散化处理,计算得到永磁同步电动机转矩和磁通的预测模型为:According to the state equation of the permanent magnet synchronous motor current and magnetic flux, the differential term is discretized, and the prediction model of the permanent magnet synchronous motor torque and magnetic flux is calculated as follows:

Figure BDA0002801587600000054
Figure BDA0002801587600000054

ψs(k+1)=ψs(k)+Ts[us(k)-Rsis(k)]ψ s (k+1)=ψ s (k)+T s [u s (k)-R s i s (k)]

Figure BDA0002801587600000061
Figure BDA0002801587600000061

其中,Ts为采样周期;is(k+1)、ψs(k+1)、Te(k+1)分别为k+1时刻定子电流矢量、定子磁链和转矩矢量;is(k)、ψs(k)、us(k)、ωr(k)分别为k时刻定子电流矢量、定子磁链、定子电压矢量和转子的电角速度;j为虚数单位,p为极对数,

Figure BDA0002801587600000062
为叉乘符号。Among them, T s is the sampling period; i s (k+1), ψ s (k+1), and Te (k+1) are the stator current vector, stator flux linkage and torque vector at time k+1, respectively; i s (k), ψ s (k), u s (k), and ω r (k) are the stator current vector, stator flux linkage, stator voltage vector and the electrical angular velocity of the rotor at time k, respectively; j is an imaginary unit, and p is Number of pole pairs,
Figure BDA0002801587600000062
is the cross product symbol.

如图1所示,三电平逆变器有三个桥臂,为A、B、C三相,每一相上有四个开关管。在正常工作时,每相有三种开关状态,对于A相,“P”状态对应的是Sa1和Sa2导通,Sa3和Sa4截止;“O”状态对应Sa2和Sa3导通,Sa1和Sa4截止;“N”状态对应Sa3和Sa4导通,Sa1和Sa2截止。因此三相桥臂有27种开关状态,有19种不同的电压矢量。不同的电压矢量对于中点电位的影响是不同的,因此在控制系统中需要设置中点电位平衡控制。根据三电平逆变器直流侧电容电压模型得到第一电容电压预测模型为:As shown in Figure 1, the three-level inverter has three bridge arms, which are three-phase A, B, and C, and each phase has four switching tubes. During normal operation, each phase has three switching states. For phase A, the "P" state corresponds to the conduction of Sa1 and Sa2, and the conduction of Sa3 and Sa4; the state of "O" corresponds to the conduction of Sa2 and Sa3, and the conduction of Sa1 and Sa4; The "N" state corresponds to that Sa3 and Sa4 are turned on, and Sa1 and Sa2 are turned off. Therefore, the three-phase bridge arm has 27 switching states and 19 different voltage vectors. Different voltage vectors have different effects on the midpoint potential, so it is necessary to set the midpoint potential balance control in the control system. According to the DC side capacitor voltage model of the three-level inverter, the first capacitor voltage prediction model is obtained as:

Figure BDA0002801587600000063
Figure BDA0002801587600000063

Figure BDA0002801587600000064
Figure BDA0002801587600000064

其中,vc1(k)、vc2(k)为k时刻采样的电容电压;vc1(k+1)、vc2(k+1)为k+1时刻预测的电容电压;ic1(k)和ic2(k)是根据k时刻逆变器的开关状态和输出电流的定义值,C1、C2为直流侧两电容,ia(k)、ib(k)、ic(k)为三相定子电流,idc(k)为直流母线电流,H1a、H1b、H1c、H2a、H2b、H2c为三相状态值,对A相来说,“P”状态时H1a为1,其余为0;“O”状态时H1b为1,其余为0;“N”状态时H1c为1,其余为0。Among them, v c1 (k) and v c2 (k) are the capacitor voltages sampled at time k; v c1 (k+1) and v c2 (k+1) are the predicted capacitor voltages at time k+1; i c1 (k ) and i c2 (k) are defined values according to the switching state and output current of the inverter at time k, C 1 and C 2 are the two capacitors on the DC side, i a (k), i b (k), ic ( k) is the three-phase stator current, i dc (k) is the DC bus current, H 1a , H 1b , H 1c , H 2a , H 2b , H 2c are the three-phase state values, for phase A, "P" In the state, H 1a is 1 and the rest are 0; in the "O" state, H 1b is 1 and the rest are 0; in the "N" state, H 1c is 1 and the rest are 0.

以永磁同步电动机的电磁转矩和磁通预测值与参考值之差的绝对值与逆变器直流侧电容电压之差绝对值及开关频率限制构建目标第一价值函数g:The target first value function g is constructed by the absolute value of the difference between the electromagnetic torque and magnetic flux predicted value of the permanent magnet synchronous motor and the reference value, the absolute value of the difference between the DC side capacitor voltage of the inverter and the switching frequency limit:

Figure BDA0002801587600000065
Figure BDA0002801587600000065

其中,

Figure BDA0002801587600000066
为转矩给定值,
Figure BDA0002801587600000067
为定子磁通给定值,λ1、λ2、λ3为权重系数,Fsw为一个控制周期内开关管发生开通和关断的动作变换次数,具体计算如下:in,
Figure BDA0002801587600000066
is the torque given value,
Figure BDA0002801587600000067
is the given value of the stator magnetic flux, λ 1 , λ 2 , and λ 3 are the weight coefficients, and F sw is the switching times of switching on and off in a control cycle. The specific calculation is as follows:

Fsw=|Sa(k)-Sa(k-1)|+|Sb(k)-Sb(k-1)|+|Sc(k)-Sc(k-1)|F sw =|S a (k)-S a (k-1)|+|S b (k)-S b (k-1)|+|S c (k)-S c (k-1)|

其中Sa(k)、Sb(k)、Sc(k)为k时刻逆变器三相的开关状态;Sa(k-1)、Sb(k-1)、Sc(k-1)为上一时刻k-1时刻逆变器三相的开关状态。Among them, S a (k), S b (k), S c (k) are the switching states of the three-phase inverter at time k; S a (k-1), S b (k-1), S c (k -1) is the switching state of the three-phase inverter at the last moment k-1.

如图2所示,三电平逆变器的空间矢量分别为3个零矢量,12个正负冗余小矢量,6个中矢量和6个大矢量,其中大矢量对中点电位无影响,因此以大矢量为中心将空间矢量图分为六个扇区。As shown in Figure 2, the space vectors of the three-level inverter are 3 zero vectors, 12 positive and negative redundant small vectors, 6 medium vectors and 6 large vectors, of which the large vectors have no effect on the midpoint potential , so the space vector is divided into six sectors centered on the large vector.

将三电平逆变器空间矢量图中的大矢量作为第一最优电压矢量的备选矢量代入到第一价值函数g中,选取使结果最小的电压矢量作为第一最优电压矢量uopt1;将与第一最优电压矢量uopt1在同一扇区的非零矢量纳入第二最优电压矢量备选集合中,每个扇区中有三个备选电压矢量;并将与第一最优电压矢量uopt1在同一扇区的零矢量作为第三电压矢量。Substitute the large vector in the three-level inverter space vector diagram as the candidate vector of the first optimal voltage vector into the first value function g, and select the voltage vector that minimizes the result as the first optimal voltage vector u opt1 ; Include the non-zero vectors in the same sector as the first optimal voltage vector u opt1 into the second optimal voltage vector candidate set, and there are three candidate voltage vectors in each sector; The zero vector of the voltage vector u opt1 in the same sector is used as the third voltage vector.

依次选取三个备选电压矢量作为第二最优电压矢量,计算第一最优电压矢量uopt1、第二最优电压矢量uopt2和第三电压矢量(零矢量)u0的作用时间,根据转矩和磁链无差拍控制原则,在一个采样周期结束后,转矩和磁链预测值等于给定值,即:Select three candidate voltage vectors in turn as the second optimal voltage vector, calculate the action time of the first optimal voltage vector u opt1 , the second optimal voltage vector u opt2 and the third voltage vector (zero vector) u 0 , according to According to the deadbeat control principle of torque and flux linkage, after a sampling period ends, the predicted value of torque and flux linkage is equal to the given value, namely:

Figure BDA0002801587600000071
Figure BDA0002801587600000071

Figure BDA0002801587600000072
Figure BDA0002801587600000072

其中

Figure BDA0002801587600000073
为转矩和磁通在第一最优电压矢量作用下的变化率;
Figure BDA0002801587600000074
为转矩和磁通在第二最优电压矢量作用下的变化率;
Figure BDA0002801587600000075
Figure BDA0002801587600000076
为转矩和磁通在零矢量作用下的变化率;
Figure BDA0002801587600000077
为第一最优电压矢量uopt1作用下的转矩和磁链的预测值;Te(k+1)opt2_n、|ψs(k+1)opt2_n|为第二最优电压矢量集合中第n个电压矢量uopt2_n作用下的转矩和磁链的预测值;Te(k+1)_0、|ψS(k+1)_0|为第三电压矢量(零矢量)作用下转矩和磁链的预测值。in
Figure BDA0002801587600000073
is the rate of change of torque and magnetic flux under the action of the first optimal voltage vector;
Figure BDA0002801587600000074
is the rate of change of torque and magnetic flux under the action of the second optimal voltage vector;
Figure BDA0002801587600000075
Figure BDA0002801587600000076
is the rate of change of torque and magnetic flux under the action of zero vector;
Figure BDA0002801587600000077
are the predicted values of torque and flux linkage under the action of the first optimal voltage vector u opt1 ; T e (k+1) opt2_n , |ψ s (k+1) opt2_n | Predicted values of torque and flux linkage under the action of n voltage vectors u opt2_n ; Te (k+1) _0 , |ψ S (k+1) _0 | are the torque under the action of the third voltage vector (zero vector) and the predicted value of the flux linkage.

第一最优电压矢量、第二最优电压矢量和第三电压矢量的作用时间分别为t1、t2和t0The action times of the first optimal voltage vector, the second optimal voltage vector and the third voltage vector are respectively t 1 , t 2 and t 0 :

t0=Ts-t1-t2 t 0 =T s -t 1 -t 2

Figure BDA0002801587600000081
Figure BDA0002801587600000081

其中,

Figure BDA0002801587600000082
S为一个量,且in,
Figure BDA0002801587600000082
S is a quantity, and

Figure BDA0002801587600000083
Figure BDA0002801587600000083

计算出t1、t2、t0后,为保证t1、t2、t0在0~Ts范围内,需作如下处理:After calculating t 1 , t 2 , and t 0 , in order to ensure that t 1 , t 2 , and t 0 are within the range of 0 to T s , the following processing is required:

(6.1)当作用时间t1、t2和t0中任意两个作用时间小于0,且另一个作用时间大于Ts,则令小于0的作用时间等于0,大于Ts的作用时间等于Ts(6.1) When any two of the action times t 1 , t 2 and t 0 are less than 0, and the other action time is greater than T s , then the action time less than 0 is equal to 0, and the action time greater than T s is equal to T s ;

(6.2)当作用时间t1+t2≤Ts,则各矢量的作用时间不变;(6.2) When the action time t 1 +t 2 ≤T s , the action time of each vector remains unchanged;

(6.3)当作用时间t1+t2≥Ts,则需要将三个矢量的作用时间代入如下公式进行修正:(6.3) When the action time t 1 +t 2 ≥T s , the action time of the three vectors needs to be substituted into the following formula for correction:

Figure BDA0002801587600000084
Figure BDA0002801587600000084

其中:t1'、t'2和t'0分别表示修正后的第一最优电压矢量、第二最优电压矢量和第三电压矢量的作用时间。Wherein: t 1 ', t' 2 and t' 0 respectively represent the action time of the corrected first optimal voltage vector, the second optimal voltage vector and the third voltage vector.

第一最优电压矢量uopt1、第二最优电压矢量uopt2_n与零矢量u0的占空比分别ηopt1、ηopt2_n、η0,根据占空比调制出新电压矢量:The duty cycles of the first optimal voltage vector u opt1 , the second optimal voltage vector u opt2_n and the zero vector u 0 are η opt1 , η opt2_n , and η 0 respectively, and a new voltage vector is modulated according to the duty cycle:

Figure BDA0002801587600000085
Figure BDA0002801587600000085

其中,uα2、uβ2为新电压矢量的α、β轴分量,uopt1_α、uopt1_β为第一最优电压矢量的α、β轴分量,uopt2_n_α、uopt2_n_β为第二最优电压矢量集合中第n个电压矢量的α、β轴分量。Among them, u α2 and u β2 are the α and β axis components of the new voltage vector, u opt1_α and u opt1_β are the α and β axis components of the first optimal voltage vector, and u opt2_n_α and u opt2_n_β are the second optimal voltage vector set The α and β axis components of the nth voltage vector in .

在三电平逆变器中,第一最优电压矢量、第二最优电压矢量和第三电压矢量对中点电位有不同的影响,因此根据三个电压矢量的作用时间重新建立直流侧两电容的预测电压,第二电容电压预测模型为:In a three-level inverter, the first optimal voltage vector, the second optimal voltage vector and the third voltage vector have different effects on the midpoint potential, so the two voltages on the DC side are re-established according to the action time of the three voltage vectors. The predicted voltage of the capacitor, the second capacitor voltage prediction model is:

Figure BDA0002801587600000091
Figure BDA0002801587600000091

其中,v'c1(k+1)、v'c2(k+1)为第一最优电压矢量、第二最优电压矢量和第三电压矢量的作用时间下k+1时刻的直流侧两电容电压,v'c1(k)、v'c2(k)为第一最优电压矢量、第二最优电压矢量和第三电压矢量的作用时间下k时刻两电容电压,ic11(k)、ic21(k)、ic31(k)、ic12(k)、ic22(k)、ic32(k)分别为第一最优电压矢量、第二最优电压矢量和第三电压矢量作用时的输出电流值。Among them, v' c1 (k+1), v' c2 (k+1) are the two voltages on the DC side at time k+1 under the action time of the first optimal voltage vector, the second optimal voltage vector and the third voltage vector Capacitor voltage, v' c1 (k), v' c2 (k) are the two capacitor voltages at time k under the action time of the first optimal voltage vector, the second optimal voltage vector and the third voltage vector, i c11 (k) , ic21 (k), ic31 (k), ic12 (k), ic22 (k), ic32 (k) are the first optimal voltage vector, the second optimal voltage vector and the third voltage vector respectively Output current value when active.

根据第二电容电压预测模型构建第二价值函数gn为:The second value function g n constructed according to the second capacitor voltage prediction model is:

Figure BDA0002801587600000092
Figure BDA0002801587600000092

将依次选取三个备选电压矢量作为第二最优电压矢量调制出的三个新电压矢量代入第二价值函数gn中,选取使结果最小的电压矢量作为第二最优电压矢量,最后将选择的第一最优电压矢量、第二最优电压矢量和第三电压矢量(零矢量)对应的开关序列及其作用时间输入脉冲调制器来控制三电平逆变器动作。The three new voltage vectors modulated by sequentially selecting three candidate voltage vectors as the second optimal voltage vector are substituted into the second value function g n , and the voltage vector with the smallest result is selected as the second optimal voltage vector, and finally the The selected switching sequences corresponding to the first optimal voltage vector, the second optimal voltage vector and the third voltage vector (zero vector) and their action time are input to the pulse modulator to control the action of the three-level inverter.

Claims (10)

1.一种三电平永磁同步电动机三矢量模型预测转矩控制方法,其特征在于,包括以下步骤:1. a three-level permanent magnet synchronous motor three-vector model predictive torque control method, is characterized in that, comprises the following steps: (1)根据永磁同步电动机的状态空间模型,构建电流和磁通的状态方程,根据电流和磁通的状态方程得到转矩和磁通的预测模型,根据三电平逆变器直流侧电容电压模型得到第一电容电压预测模型;根据第一电容电压预测模型构建第一价值函数;(1) According to the state space model of the permanent magnet synchronous motor, the state equations of current and magnetic flux are constructed, and the prediction model of torque and magnetic flux is obtained according to the state equations of current and magnetic flux. According to the DC side capacitance of the three-level inverter obtaining a first capacitor voltage prediction model from the voltage model; constructing a first value function according to the first capacitor voltage prediction model; (2)采样得到k时刻的电流;(2) Sampling to obtain the current at time k; (3)将采样取得的数据代入到转矩和磁通的预测模型中,得到下一时刻k+1时刻的转矩和磁链;(3) Substitute the data obtained by sampling into the prediction model of torque and magnetic flux, and obtain the torque and flux linkage at the next moment k+1; (4)将三电平逆变器空间矢量图中的大矢量作为第一最优电压矢量的备选矢量代入到第一价值函数中,选出第一最优电压矢量;(4) Substitute the large vector in the three-level inverter space vector diagram into the first value function as the candidate vector of the first optimal voltage vector, and select the first optimal voltage vector; (5)将第一最优电压矢量对应的扇区中的三个非零矢量作为第二最优电压矢量的备选矢量,零矢量作为第三电压矢量;(5) Three non-zero vectors in the sector corresponding to the first optimal voltage vector are used as candidate vectors for the second optimal voltage vector, and the zero vector is used as the third voltage vector; (6)依次选取三个备选矢量作为第二最优电压矢量,计算第一最优电压矢量、第二最优电压矢量和第三电压矢量的作用时间,并求出占空比,根据占空比调制出新电压矢量;(6) Select three candidate vectors in turn as the second optimal voltage vector, calculate the action time of the first optimal voltage vector, the second optimal voltage vector and the third voltage vector, and obtain the duty ratio. The empty ratio modulates a new voltage vector; (7)根据第一最优电压矢量、第二最优电压矢量和第三电压矢量的作用时间建立第二电容电压预测模型,并根据第二电容电压预测模型构建第二价值函数;(7) establishing a second capacitance voltage prediction model according to the action time of the first optimal voltage vector, the second optimal voltage vector and the third voltage vector, and building a second cost function according to the second capacitance voltage prediction model; (8)将步骤6中,调制出的三个新电压矢量代入到第二价值函数中选出最优的第二最优电压矢量;(8) Substitute the three new voltage vectors modulated in step 6 into the second value function to select the optimal second optimal voltage vector; (9)将步骤(4)所选的第一最优电压矢量、步骤(8)所选的第二最优电压矢量和步骤(5)所选的第三电压矢量对应的开关序列及对应的作用时间输入脉冲调制器控制三电平逆变器动作。(9) Compare the switching sequence corresponding to the first optimal voltage vector selected in step (4), the second optimal voltage vector selected in step (8) and the third voltage vector selected in step (5) and the corresponding The action time input pulse modulator controls the action of the three-level inverter. 2.根据权利要求1所述的三矢量模型预测转矩控制方法,其特征在于,所述步骤1中电流和磁通的状态方程为:2. The three-vector model predictive torque control method according to claim 1, wherein the state equation of the current and the magnetic flux in the step 1 is:
Figure FDA0003565004000000011
Figure FDA0003565004000000011
Figure FDA0003565004000000012
Figure FDA0003565004000000012
其中,iα、iβ分别表示定子电流矢量在α、β轴的分量;Ls表示定子电感;uα、uβ分别表示定子电压矢量在α、β轴下的分量;Rs表示定子电阻;ωr表示电动机的转子电角速度;ψα、ψβ分别表示电机磁链矢量在α、β轴下的分量。Among them, i α and i β represent the components of the stator current vector on the α and β axes respectively; L s represents the stator inductance; u α and u β represent the components of the stator voltage vector on the α and β axes respectively; R s represents the stator resistance ; ω r represents the rotor electrical angular velocity of the motor; ψ α , ψ β represent the components of the motor flux linkage vector under the α and β axes, respectively.
3.根据权利要求2所述的三矢量模型预测转矩控制方法,其特征在于,所述步骤1中,转矩和磁通的预测模型为:3. The three-vector model predictive torque control method according to claim 2, wherein in the step 1, the predictive models of torque and magnetic flux are:
Figure FDA0003565004000000021
Figure FDA0003565004000000021
ψs(k+1)=ψs(k)+Ts[us(k)-Rsis(k)]ψ s (k+1)=ψ s (k)+T s [u s (k)-R s i s (k)]
Figure FDA0003565004000000022
Figure FDA0003565004000000022
其中,Ts为采样周期;is(k+1)、ψs(k+1)、Te(k+1)分别为k+1时刻定子电流矢量、定子磁链和转矩矢量;is(k)、ψs(k)、us(k)、ωr(k)分别为k时刻定子电流矢量、定子磁链、定子电压矢量和转子的电角速度;j为虚数单位,p为极对数,
Figure FDA0003565004000000023
为叉乘符号。
Among them, T s is the sampling period; i s (k+1), ψ s (k+1), and Te (k+1) are the stator current vector, stator flux linkage and torque vector at time k+1, respectively; i s (k), ψ s (k), u s (k), and ω r (k) are the stator current vector, stator flux linkage, stator voltage vector and the electrical angular velocity of the rotor at time k, respectively; j is an imaginary unit, and p is Number of pole pairs,
Figure FDA0003565004000000023
is the cross product symbol.
4.根据权利要求3所述的三矢量模型预测转矩控制方法,其特征在于,所述步骤1中,电容电压预测模型为:4. The three-vector model predicted torque control method according to claim 3, wherein in the step 1, the capacitor voltage prediction model is:
Figure FDA0003565004000000024
Figure FDA0003565004000000024
其中,vc1(k)、vc2(k)为k时刻采样的电容电压;vc1(k+1)、vc2(k+1)为k+1时刻预测的电容电压;ic1(k)和ic2(k)是根据k时刻逆变器的开关状态和输出电流的定义值,C1、C2为直流侧两电容。Among them, v c1 (k) and v c2 (k) are the capacitor voltages sampled at time k; v c1 (k+1) and v c2 (k+1) are the predicted capacitor voltages at time k+1; i c1 (k ) and i c2 (k) are defined values according to the switching state and output current of the inverter at time k, C 1 and C 2 are the two capacitors on the DC side.
5.根据权利要求4所述的三矢量模型预测转矩控制方法,其特征在于,所述步骤4中,第一价值函数为:5. The three-vector model predictive torque control method according to claim 4, wherein in the step 4, the first value function is:
Figure FDA0003565004000000025
Figure FDA0003565004000000025
其中,
Figure FDA0003565004000000026
为转矩给定值,
Figure FDA0003565004000000027
为定子磁通给定值,λ1、λ2、λ3为权重系数,Fsw为一个控制周期内开关管发生开通和关断的动作变换次数。
in,
Figure FDA0003565004000000026
is the torque given value,
Figure FDA0003565004000000027
is the given value of the stator magnetic flux, λ 1 , λ 2 , λ 3 are the weight coefficients, and F sw is the switching times of switching on and off in a control cycle.
6.根据权利要求5所述的三矢量模型预测转矩控制方法,其特征在于,所述步骤6中,第一最优电压矢量、第二最优电压矢量和第三电压矢量的作用时间分别为t1、t2和t06 . The three-vector model predictive torque control method according to claim 5 , wherein in the step 6, the action times of the first optimal voltage vector, the second optimal voltage vector and the third voltage vector are respectively 6 . for t 1 , t 2 and t 0 : t0=Ts-t1-t2 t 0 =T s -t 1 -t 2
Figure FDA0003565004000000031
Figure FDA0003565004000000031
其中,
Figure FDA0003565004000000032
S为一个量,且
Figure FDA0003565004000000033
Figure FDA0003565004000000034
为转矩和磁通在第一最优电压矢量作用下的变化率;
Figure FDA0003565004000000035
Figure FDA0003565004000000036
为转矩和磁通在第二最优电压矢量作用下的变化率;
Figure FDA0003565004000000037
为转矩和磁通在零矢量作用下的变化率。
in,
Figure FDA0003565004000000032
S is a quantity, and
Figure FDA0003565004000000033
Figure FDA0003565004000000034
is the rate of change of torque and magnetic flux under the action of the first optimal voltage vector;
Figure FDA0003565004000000035
Figure FDA0003565004000000036
is the rate of change of torque and magnetic flux under the action of the second optimal voltage vector;
Figure FDA0003565004000000037
is the rate of change of torque and flux under the action of zero vector.
7.根据权利要求6所述的三矢量模型预测转矩控制方法,其特征在于,t1、t2、t0在0~Ts范围内,且为保证t1、t2、t0在0~Ts范围内,需作如下处理:7 . The three-vector model predictive torque control method according to claim 6 , wherein t 1 , t 2 , and t 0 are within the range of 0 to T s , and in order to ensure that t 1 , t 2 , and t 0 are within the range of t 1 , t 2 , and t 0 . In the range of 0~T s , the following processing is required: (6.1)当作用时间t1、t2和t0中任意两个作用时间小于0,且另一个作用时间大于Ts,则令小于0的作用时间等于0,大于Ts的作用时间等于Ts(6.1) When any two of the action times t 1 , t 2 and t 0 are less than 0, and the other action time is greater than T s , then the action time less than 0 is equal to 0, and the action time greater than T s is equal to T s ; (6.2)当作用时间t1+t2≤Ts,则各矢量的作用时间不变;(6.2) When the action time t 1 +t 2 ≤T s , the action time of each vector remains unchanged; (6.3)当作用时间t1+t2≥Ts,则需要将三个矢量的作用时间代入如下公式进行修正:(6.3) When the action time t 1 +t 2 ≥T s , the action time of the three vectors needs to be substituted into the following formula for correction:
Figure FDA0003565004000000038
Figure FDA0003565004000000038
其中:t'1、t'2和t'0分别表示修正后的第一最优电压矢量、第二最优电压矢量和第三电压矢量的作用时间。Wherein: t' 1 , t' 2 and t' 0 respectively represent the action time of the corrected first optimal voltage vector, the second optimal voltage vector and the third voltage vector.
8.根据权利要求7所述的三矢量模型预测转矩控制方法,其特征在于,所述步骤6中,新电压矢量的调制公式:8. The three-vector model predictive torque control method according to claim 7, wherein in the step 6, the modulation formula of the new voltage vector:
Figure FDA0003565004000000039
Figure FDA0003565004000000039
其中,uα2、uβ2为新电压矢量的α、β轴分量,uopt1_α、uopt1_β为第一最优电压矢量的α、β轴分量,uopt2_n_α、uopt2_n_β为第二最优电压矢量集合中第n个电压矢量的α、β轴分量,第一最优电压矢量uopt1、第二最优电压矢量uopt2_n与零矢量u0的占空比分别ηopt1、ηopt2_n、η0Among them, u α2 and u β2 are the α and β axis components of the new voltage vector, u opt1_α and u opt1_β are the α and β axis components of the first optimal voltage vector, and u opt2_n_α and u opt2_n_β are the second optimal voltage vector set In the α and β axis components of the nth voltage vector, the duty ratios of the first optimal voltage vector u opt1 , the second optimal voltage vector u opt2_n and the zero vector u 0 are η opt1 , η opt2_n , and η 0 , respectively.
9.根据权利要求8所述的三矢量模型预测转矩控制方法,其特征在于,所述步骤7中第二电容电压预测模型为:9. The three-vector model predicted torque control method according to claim 8, wherein the second capacitor voltage prediction model in the step 7 is:
Figure FDA0003565004000000041
Figure FDA0003565004000000041
其中,v'c1(k+1)、v'c2(k+1)为第一最优电压矢量、第二最优电压矢量和第三电压矢量的作用时间下k+1时刻的直流侧两电容电压,v'c1(k)、v'c2(k)为第一最优电压矢量、第二最优电压矢量和第三电压矢量的作用时间下k时刻两电容电压,ic11(k)、ic21(k)、ic31(k)、ic12(k)、ic22(k)、ic32(k)分别为第一最优电压矢量、第二最优电压矢量和第三电压矢量作用时直流侧两电容的输出电流值。Among them, v' c1 (k+1), v' c2 (k+1) are the two voltages on the DC side at time k+1 under the action time of the first optimal voltage vector, the second optimal voltage vector and the third voltage vector Capacitor voltage, v' c1 (k), v' c2 (k) are the two capacitor voltages at time k under the action time of the first optimal voltage vector, the second optimal voltage vector and the third voltage vector, i c11 (k) , ic21 (k), ic31 (k), ic12 (k), ic22 (k), ic32 (k) are the first optimal voltage vector, the second optimal voltage vector and the third voltage vector respectively The output current value of the two capacitors on the DC side when acting.
10.根据权利要求9所述的三矢量模型预测转矩控制方法,其特征在于,所述步骤7中第二价值函数为:10. The three-vector model predictive torque control method according to claim 9, wherein the second value function in the step 7 is:
Figure FDA0003565004000000042
Figure FDA0003565004000000042
其中,gn为第二价值函数。Among them, g n is the second value function.
CN202011352040.4A 2020-11-27 2020-11-27 Three-vector model prediction torque control method for three-level permanent magnet synchronous motor Active CN112583315B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011352040.4A CN112583315B (en) 2020-11-27 2020-11-27 Three-vector model prediction torque control method for three-level permanent magnet synchronous motor

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011352040.4A CN112583315B (en) 2020-11-27 2020-11-27 Three-vector model prediction torque control method for three-level permanent magnet synchronous motor

Publications (2)

Publication Number Publication Date
CN112583315A CN112583315A (en) 2021-03-30
CN112583315B true CN112583315B (en) 2022-05-17

Family

ID=75124006

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011352040.4A Active CN112583315B (en) 2020-11-27 2020-11-27 Three-vector model prediction torque control method for three-level permanent magnet synchronous motor

Country Status (1)

Country Link
CN (1) CN112583315B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112953341A (en) * 2021-04-22 2021-06-11 哈尔滨理工大学 Three-vector model predictive control algorithm for zero common-mode voltage of permanent magnet synchronous motor
CN113179065B (en) * 2021-05-12 2022-08-19 西南交通大学 Permanent magnet synchronous motor model prediction pulse sequence control method
CN113708688B (en) * 2021-08-31 2023-06-30 西南交通大学 Permanent magnet motor vector-reduction model predictive control method

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2546979A1 (en) * 2011-07-15 2013-01-16 ABB Research Ltd. Method for controlling harmonics and resonances in an inverter
CN106936356B (en) * 2017-04-24 2023-03-28 东南大学盐城新能源汽车研究院 Vector screening and duty ratio combined motor model prediction control system and method
CN110350568B (en) * 2019-06-26 2022-12-16 武汉大学 Universal model prediction optimization control method under unbalanced power grid condition
CN110707978B (en) * 2019-10-12 2021-06-11 南通大学 Three-level permanent magnet synchronous motor model prediction control method considering vector partition
CN110867887B (en) * 2019-10-16 2022-09-06 南京理工大学 NPC three-level inverter optimal sequence model prediction control device and method

Also Published As

Publication number Publication date
CN112583315A (en) 2021-03-30

Similar Documents

Publication Publication Date Title
CN112583315B (en) Three-vector model prediction torque control method for three-level permanent magnet synchronous motor
CN106936356B (en) Vector screening and duty ratio combined motor model prediction control system and method
CN109660170B (en) A high-reliability current predictive control method and system for a permanent magnet synchronous motor
CN107565865B (en) A kind of fault-tolerant double vector prediction control method and device of six-phase permanent-magnet motor
CN110557074B (en) Predictive Current Control Method for Cascaded H-Bridge Inverter Permanent Magnet Motor Systems
CN112737444B (en) Dual three-phase permanent magnet synchronous motor control method with alternate execution of sampling and control programs
CN113708688B (en) Permanent magnet motor vector-reduction model predictive control method
CN109274303A (en) Asynchronous motor model predictive control method based on extended control set
CN111555684B (en) Multi-step Model Predictive Torque Control Method for Permanent Magnet Synchronous Motors with Variable Switch Points
CN110492821B (en) Permanent magnet motor direct flux linkage control method based on unfixed vector action time
CN110829940A (en) A Torque Distribution Model Predictive Control Method for Switched Reluctance Motor
CN113179065B (en) Permanent magnet synchronous motor model prediction pulse sequence control method
CN114584040A (en) Predictive torque control method of permanent magnet synchronous motor based on discrete space vector modulation
CN111130419A (en) Predictive flux linkage control method for permanent magnet motor based on extended step size and variable action time
CN110401393A (en) Tolerant layered sequence model predictive control method, device and storage medium
CN111969916A (en) Multi-current-difference-updating model-free prediction current control method for permanent magnet synchronous motor
CN112821816B (en) Predictive current control method based on PMSM model of NPC three-level inverter
CN115622040A (en) Multi-step model predictive control method and DSP controller based on hysteresis control
CN117526790A (en) Five-level model predictive control method for open winding permanent magnet synchronous motor
CN117375476A (en) Three-level model-free prediction current control method for permanent magnet synchronous motor
CN117175988A (en) Global predictive control method for quasi-Z source inverter-permanent magnet synchronous motor system without weight coefficient
CN108649852B (en) Permanent magnet synchronous motor control method for improving current loop
CN111106773A (en) Permanent magnet synchronous motor model prediction control method based on optimized duty ratio
CN116667732A (en) A model predictive current control method for permanent magnet synchronous motor with three-level inverter
CN113285634B (en) High-speed field weakening control method and system for permanent magnet synchronous motor based on multi-step zero-delay model prediction

Legal Events

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