[go: up one dir, main page]

CN116094400A - A MTPA predictive control method for permanent magnet synchronous motor - Google Patents

A MTPA predictive control method for permanent magnet synchronous motor Download PDF

Info

Publication number
CN116094400A
CN116094400A CN202310250531.5A CN202310250531A CN116094400A CN 116094400 A CN116094400 A CN 116094400A CN 202310250531 A CN202310250531 A CN 202310250531A CN 116094400 A CN116094400 A CN 116094400A
Authority
CN
China
Prior art keywords
current
axis
mtpa
permanent magnet
magnet synchronous
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.)
Withdrawn
Application number
CN202310250531.5A
Other languages
Chinese (zh)
Inventor
高丽华
王琰
马凯
朱丽琴
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing Polytechnic Institute
Original Assignee
Nanjing Polytechnic Institute
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 Nanjing Polytechnic Institute filed Critical Nanjing Polytechnic Institute
Priority to CN202310250531.5A priority Critical patent/CN116094400A/en
Publication of CN116094400A publication Critical patent/CN116094400A/en
Withdrawn legal-status Critical Current

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
    • H02P21/22Current control, e.g. using a current control loop
    • 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/0003Control strategies in general, e.g. linear type, e.g. P, PI, PID, using robust control
    • 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
    • H02P25/00Arrangements or methods for the control of AC motors characterised by the kind of AC motor or by structural details
    • H02P25/02Arrangements or methods for the control of AC motors characterised by the kind of AC motor or by structural details characterised by the kind of motor
    • H02P25/022Synchronous motors
    • H02P25/024Synchronous motors controlled by supply frequency
    • 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
    • H02P27/00Arrangements or methods for the control of AC motors characterised by the kind of supply voltage
    • H02P27/04Arrangements or methods for the control of AC motors characterised by the kind of supply voltage using variable-frequency supply voltage, e.g. inverter or converter supply voltage
    • H02P27/06Arrangements or methods for the control of AC motors characterised by the kind of supply voltage using variable-frequency supply voltage, e.g. inverter or converter supply voltage using DC to AC converters or inverters
    • H02P27/08Arrangements or methods for the control of AC motors characterised by the kind of supply voltage using variable-frequency supply voltage, e.g. inverter or converter supply voltage using DC to AC converters or inverters with pulse width modulation
    • 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
    • H02P2207/00Indexing scheme relating to controlling arrangements characterised by the type of motor
    • H02P2207/05Synchronous machines, e.g. with permanent magnets or DC excitation

Landscapes

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

Abstract

本发明公开了一种永磁同步电机MTPA预测控制方法,包括:建立永磁同步电机的基于MTPA的电流预测的多项目价值函数;将当前时刻的采样的永磁同步电机电流与电压值带入价值函数求解得到当前时刻满足电机最优控制的变量,并根据该最优控制的变量中的电流值计算出下一时刻即k时刻的电流值;将电机转速参考值和得到电流值输入PI调节器后获得该时刻的电磁转矩参考值;将获得的电磁转矩参考值、预测到的下一时刻的d轴预测电流作为MTPA模块的输入,MTPA模块输出其工作时的q轴电流以及d轴电流。本发明在电流预测控制的基础上,通过上一个周期预测的当前周期的d轴电流近似获取MTPA工作点电流,有效简化公式法中电流计算方法。

Figure 202310250531

The invention discloses a MTPA predictive control method for a permanent magnet synchronous motor, which includes: establishing a multi-item value function based on MTPA current prediction of the permanent magnet synchronous motor; bringing the current and voltage values of the permanent magnet synchronous motor sampled at the current moment into Solve the value function to obtain the variable that satisfies the optimal control of the motor at the current moment, and calculate the current value at the next moment, that is, the k time according to the current value in the variable of the optimal control; input the motor speed reference value and the obtained current value into the PI adjustment The electromagnetic torque reference value at this moment is obtained after the controller; the obtained electromagnetic torque reference value and the predicted d-axis predicted current at the next moment are used as the input of the MTPA module, and the MTPA module outputs the q-axis current and d shaft current. On the basis of current predictive control, the present invention approximates the current of the MTPA working point through the d-axis current of the current cycle predicted in the previous cycle, and effectively simplifies the current calculation method in the formula method.

Figure 202310250531

Description

一种永磁同步电机MTPA预测控制方法A predictive control method for MTPA of permanent magnet synchronous motor

技术领域Technical Field

本发明涉及交流伺服控制技术领域,具体涉及一种永磁同步电机MTPA预测控制方法。The present invention relates to the technical field of AC servo control, and in particular to a permanent magnet synchronous motor MTPA predictive control method.

背景技术Background Art

永磁同步电机因其结构简单、功率密度高、重量轻、低损耗和效率高等性能,在汽车工业等多个领域得到广泛的应用。随着控制理论和计算机技术的发展,越来越多先进的控制方法被应用到永磁同步电机的控制中,常用的如改进PI控制、滑模控制、MTPA控制、预测控制等。Permanent magnet synchronous motors are widely used in many fields such as the automotive industry due to their simple structure, high power density, light weight, low loss and high efficiency. With the development of control theory and computer technology, more and more advanced control methods are applied to the control of permanent magnet synchronous motors, such as improved PI control, sliding mode control, MTPA control, predictive control, etc.

最大转矩电流比MTPA(Maximum Torque per Ampere)控制是一种提升电机运行效率的控制方法,它使得电流一定的情况下输出最大电磁转矩。MTPA控制主要分为三类:查表法,信号注入法和公式法。公式法是MTPA控制中常用的一种方法,公式法MTPA通过建立电磁转矩和定子电流约束条件,并通过极值求解得到电磁转矩与交轴、直轴电流之间的关系,进而得到直轴和交轴电流给定值。现有的公式法如附图2所示中,不仅包括原有公式法MTPA、还包括超前角公式法MTPA;MTPA法直观易懂,但其中通过转矩获得电流的公式中含有高次多项式、反三角函数等运算,计算复杂且量大。Maximum torque per current ratio MTPA (Maximum Torque per Ampere) control is a control method to improve the operating efficiency of the motor. It enables the maximum electromagnetic torque to be output when the current is constant. MTPA control is mainly divided into three categories: table lookup method, signal injection method and formula method. The formula method is a commonly used method in MTPA control. The formula method MTPA establishes electromagnetic torque and stator current constraints, and obtains the relationship between the electromagnetic torque and the quadrature axis and direct axis currents through extreme value solution, and then obtains the direct axis and quadrature axis current given values. The existing formula method is shown in Figure 2, which includes not only the original formula method MTPA, but also the lead angle formula method MTPA; the MTPA method is intuitive and easy to understand, but the formula for obtaining current through torque contains high-order polynomials, inverse trigonometric functions and other operations, and the calculation is complex and large.

模型预测控制MPC(Model Predictive Control)是在建立在系统数学模型基础上,实现在当前时刻预测下一时刻的状态。而永磁同步电机连续电流预测控制是将参考电流与运行电流之间的误差进行约束求解,然后将下一个控制周期的最优控制电压用于电机控制,动态性能好。因此需要更好的方法来解决如何在永磁同步电机电流预测控制的基础上研究一种计算更为简单的MTPA控制方法。Model Predictive Control (MPC) is based on the mathematical model of the system to predict the state of the next moment at the current moment. The continuous current predictive control of permanent magnet synchronous motors solves the error between the reference current and the operating current, and then uses the optimal control voltage of the next control cycle for motor control, which has good dynamic performance. Therefore, a better method is needed to solve how to study a MTPA control method with simpler calculations based on the current predictive control of permanent magnet synchronous motors.

发明内容Summary of the invention

1.所要解决的技术问题:1. Technical problems to be solved:

针对上述技术问题,本发明提供一种永磁同步电机MTPA预测控制方法,本方法能够简化MTPA工作点求解方法,同时在电流控制中引入基于MTPA的预测控制方案,确保电流轨迹能够按照最大转矩电流比曲线运行,提升电机运行效率。In view of the above technical problems, the present invention provides a permanent magnet synchronous motor MTPA predictive control method, which can simplify the MTPA operating point solution method and introduce a MTPA-based predictive control scheme into current control to ensure that the current trajectory can operate according to the maximum torque current ratio curve and improve the motor operation efficiency.

2.技术方案:2. Technical solution:

一种永磁同步电机MTPA预测控制方法,其特征在于:包括以下步骤:A permanent magnet synchronous motor MTPA predictive control method, characterized in that it comprises the following steps:

步骤一:在d-q旋转坐标系下建立永磁同步电机的基于MTPA的电流预测的多项目价值函数;所述多项目价值函数包括电流误差优化、最大转矩电流比优化以及电压误差优化;Step 1: Establish a multi-item value function for current prediction based on MTPA of a permanent magnet synchronous motor in a d-q rotating coordinate system; the multi-item value function includes current error optimization, maximum torque current ratio optimization, and voltage error optimization;

步骤二:通过价值函数求出当前k时刻电机控制最优电压,并根据该最优控制电压和永磁同步电机得电流预测模型预测下一个k+1时刻的d轴电流值;Step 2: Calculate the optimal voltage of the motor control at the current k moment through the value function, and predict the d-axis current value at the next k+1 moment based on the optimal control voltage and the current prediction model of the permanent magnet synchronous motor;

步骤三:将电机转速参考值和电机当前运行转速值输入PI调节器后获得该时刻的电磁转矩参考值Te *(k);Step 3: Input the motor speed reference value and the current running speed value of the motor into the PI regulator to obtain the electromagnetic torque reference value T e * (k) at that moment;

步骤四:将步骤三获得的电磁转矩参考值、根据在k-1时刻预测到的k时刻的d轴预测电流作为MTPA模块的输入,MTPA模块输出其工作时的q轴电流以及d轴电流。Step 4: The electromagnetic torque reference value obtained in step 3 and the predicted d-axis current at time k predicted at time k-1 are used as inputs of the MTPA module, and the MTPA module outputs the q-axis current and d-axis current during its operation.

进一步地,所述步骤一具体包括:Furthermore, the step 1 specifically includes:

S11:永磁同步电机电流方程为:S11: The current equation of the permanent magnet synchronous motor is:

Figure BDA0004127752170000021
Figure BDA0004127752170000021

Figure BDA0004127752170000022
Figure BDA0004127752170000022

公式(1)中id为d轴电流,i q为q轴电流;Rs为电阻;Ld为d轴磁导、Lq为q轴磁导;P为电机极对数;ud为d轴电压,uq为q轴电压;ωm为机械转速;φf为磁链;In formula (1), i d is the d-axis current, i q is the q-axis current; R s is the resistance; L d is the d-axis permeance, Lq is the q-axis permeance; P is the number of motor pole pairs; ud is the d-axis voltage, u q is the q-axis voltage; ω m is the mechanical speed; φ f is the flux linkage;

S12:将永磁同步电机电流方程离散化;得到永磁同步电机的电流离散模型如下式(2):S12: Discretize the permanent magnet synchronous motor current equation; obtain the current discrete model of the permanent magnet synchronous motor as follows (2):

Figure BDA0004127752170000023
Figure BDA0004127752170000023

Figure BDA0004127752170000024
Figure BDA0004127752170000024

(2)式中,T为采样周期,k表示k个采样时刻;(2) In the formula, T is the sampling period, k represents k sampling moments;

S13:将公式(2)改写为状态空间模型得到如下式(3)的状态空间方程:S13: Rewrite formula (2) into a state space model to obtain the state space equation of formula (3):

x(k+1)=Ax(k)+Bu(k)x(k+1)=Ax(k)+Bu(k)

y(k)=Cx(k) (3);y(k)=Cx(k) (3);

(3)式中,

Figure BDA0004127752170000025
(3) In the formula,
Figure BDA0004127752170000025

Figure BDA0004127752170000031
Figure BDA0004127752170000031

Figure BDA0004127752170000032
Figure BDA0004127752170000032

S14:永磁同步电机采用MTPA模式运行时轨迹为MTPA控制曲线,即d、q轴电流满足公式:

Figure BDA0004127752170000033
基于该MTPA控制曲线的原理建立如下式(4)的基于MTPA电流预测的多项目价值函数;其中多项目价值优化包括d、q轴电流误差优化,d、q轴电压误差优化以及MTPA曲线误差优化;所述d、q轴电流误差为经过MTPA模块得到的电机d、q轴参考电流与电机当前运行d、q轴电流误差,通过该项控制电机电流能够跟踪MTPA控制工作点的参考电流运行;所述d、q轴电压误差为永磁同步电机控制电压增量误差达到最小;所述MTPA曲线误差为保证电流控制按照预设的MTPA曲线运行;S14: When the permanent magnet synchronous motor is operated in MTPA mode, the trajectory is the MTPA control curve, that is, the d-axis and q-axis currents satisfy the formula:
Figure BDA0004127752170000033
Based on the principle of the MTPA control curve, a multi-project value function based on MTPA current prediction is established as shown in the following formula (4); wherein the multi-project value optimization includes d and q axis current error optimization, d and q axis voltage error optimization and MTPA curve error optimization; the d and q axis current errors are the motor d and q axis reference currents obtained by the MTPA module and the motor's current d and q axis current errors, and the motor current can track the reference current operation of the MTPA control working point through this control; the d and q axis voltage errors are the permanent magnet synchronous motor control voltage increment errors that are minimized; the MTPA curve error is to ensure that the current control runs according to the preset MTPA curve;

Figure BDA0004127752170000034
Figure BDA0004127752170000034

(4)式中,kid为当前d轴电流与d轴参考电流之

Figure BDA0004127752170000035
间的误差权重;kiq为q轴电流与q轴参考电流之
Figure BDA0004127752170000036
间的误差权重;kmtpa为当电机采用MTPA运行时d、q轴电流关系误差的权重;kud和kuq分别为d、q轴电压在k时刻增量的权重。(4) where kid is the difference between the current d-axis current and the d-axis reference current.
Figure BDA0004127752170000035
The error weight between q-axis current and q-axis reference current; kiq is the error between q-axis current and q-axis reference current
Figure BDA0004127752170000036
kmtpa is the weight of the error between the d-axis and q-axis currents when the motor adopts MTPA operation; kud and kuq are the weights of the increments of the d-axis and q-axis voltages at time k, respectively.

进一步地,步骤二中通过对价值函数分别进行对△ud和△uq求偏导数,即:Furthermore, in step 2, the partial derivatives of the value function are calculated for △u d and △u q , namely:

Figure BDA0004127752170000037
Figure BDA0004127752170000037

Figure BDA0004127752170000038
Figure BDA0004127752170000038

从而得到当前k时刻的最优控制电压增量△ud和△uq;然后根据公式(5)获得k时刻最优电压,并将求得的该电压用于永磁同步电机控制;Thus, the optimal control voltage increments △u d and △u q at the current time k are obtained; then the optimal voltage at the time k is obtained according to formula (5), and the obtained voltage is used for permanent magnet synchronous motor control;

ud(k)=ud(k-1)+△ud(k)u d (k)=u d (k-1)+△u d (k)

uq(k)=uq(k-1)+△uq(k) (5);u q (k)=u q (k-1)+Δu q (k) (5);

将上式中获取的d、q轴电压带入公式(6),得到k+1时刻的d轴电流idmiddle(k+1),作为MTPA工作点求解的中间变量:Substituting the d-axis and q-axis voltages obtained in the above formula into formula (6), we can obtain the d-axis current i dmiddle (k+1) at time k+1, which is used as the intermediate variable for solving the MTPA operating point:

Figure BDA0004127752170000041
Figure BDA0004127752170000041

进一步地,步骤三具体为:Furthermore, step three is specifically as follows:

永磁同步电机速度公式如下:The speed formula of permanent magnet synchronous motor is as follows:

Figure BDA0004127752170000042
Figure BDA0004127752170000042

利用上式公式,将电机速度参考值和当前时刻电机的速度反馈值输入PI调节器并得到电机的电磁转矩参考值Te *(k)。By using the above formula, the motor speed reference value and the speed feedback value of the motor at the current moment are input into the PI regulator to obtain the electromagnetic torque reference value Te * (k) of the motor.

进一步地,将步骤三获得的电磁转矩参考值Te *(k)、步骤一中的d轴预测电流idpre(k+1)作为MTPA算法模块的输入求解电机MTPA控制时的工作点,具体如下:Furthermore, the electromagnetic torque reference value Te * (k) obtained in step 3 and the d-axis predicted current i dpre (k+1) in step 1 are used as inputs of the MTPA algorithm module to solve the working point of the motor during MTPA control, as follows:

S41:永磁同步电机电磁转矩方程如(8)所示,S41: The electromagnetic torque equation of the permanent magnet synchronous motor is shown in (8),

Figure BDA0004127752170000043
Figure BDA0004127752170000043

将电磁转矩Te *(k)、预测电流idmiddle(k)代入公式(8),可以得到MTPA工作时的q轴电流:Substituting the electromagnetic torque Te * (k) and the predicted current i dmiddle (k) into formula (8), the q-axis current when MTPA is working can be obtained:

Figure BDA0004127752170000044
Figure BDA0004127752170000044

S42:由于当永磁同步电机采用MTPA模式运行时,其d、q轴电流满足下式(10):S42: When the permanent magnet synchronous motor is operated in MTPA mode, its d-axis and q-axis currents satisfy the following equation (10):

Figure BDA0004127752170000045
Figure BDA0004127752170000045

因此将得到的中iqmtpa(k)带入公式(10),即得到MTPA模板工作时的d轴工作电流,如下式(11):Therefore, the obtained i qmtpa (k) is substituted into formula (10), and the d-axis working current of the MTPA template when working is obtained, which is as follows (11):

Figure BDA0004127752170000046
Figure BDA0004127752170000046

3.有益效果:3. Beneficial effects:

(1)本发明提供了的一种MTPA工作点电流的计算方法,在电流预测控制的基础上,通过上一个周期预测的当前周期的d轴电流近似获取MTPA工作点电流,附图2中公开的是现有技术中的公式法、超前角公式法以及本申请中公开的改进MTPA的计算公式,从中不难看出本申请公开的方法能够简化公式法中电流计算方法,消除了传统的公式法中的开方等运算。(1) The present invention provides a method for calculating the MTPA operating point current. On the basis of current prediction control, the MTPA operating point current is approximately obtained by predicting the d-axis current of the current cycle from the previous cycle. Attached Figure 2 discloses the formula method in the prior art, the lead angle formula method and the calculation formula of the improved MTPA disclosed in the present application. It is not difficult to see that the method disclosed in the present application can simplify the current calculation method in the formula method and eliminate the square root operation and the like in the traditional formula method.

(2)本发明还提供了一种MTPA电流预测控制方法,在电流预测价值函数中除了电流误差、电压误差外,增加了MTPA工作轨迹误差,使得永磁同步电机d轴、q轴电流尽可能沿着最大转矩电流比轨迹运行。同时,电流预测中获得的下一时刻的d轴电流作为MTPA工作点求解的变量,用于简化MTPA工作点求解的公式。(2) The present invention also provides an MTPA current prediction control method. In addition to the current error and voltage error, the MTPA working trajectory error is added to the current prediction value function, so that the d-axis and q-axis currents of the permanent magnet synchronous motor run along the maximum torque current ratio trajectory as much as possible. At the same time, the d-axis current at the next moment obtained in the current prediction is used as a variable for solving the MTPA working point, which is used to simplify the formula for solving the MTPA working point.

附图说明BRIEF DESCRIPTION OF THE DRAWINGS

图1为本发明的系统控制框图;FIG1 is a system control block diagram of the present invention;

图2为本发明中MTPA公式与现有技术中的MTPA公式。FIG. 2 shows the MTPA formula in the present invention and the MTPA formula in the prior art.

具体实施方式DETAILED DESCRIPTION

下面结合附图对本发明进行具体的说明。The present invention will be described in detail below with reference to the accompanying drawings.

现有的技术方案,如公式法、查表法等,在准确性、计算量、实施速度等方面具有一定的不足。Existing technical solutions, such as formula method, table lookup method, etc., have certain shortcomings in terms of accuracy, calculation amount, implementation speed, etc.

公式求解方法是应用广泛的一种MTPA工作点电流的求解方法。现有的MTPA公式法主要有电流求解法、超前角求解法等,如附图2中所示。电流求解法以速度环输出电磁转矩,并根据公式求解出d、q轴电流,用于MTPA控制,该公式较为计算复杂;超前角电流法相对于直接公式法公式形式简单,但设计到三角函数,计算也较为复杂。本申请提出的改进的MTPA控制算法首先采用预测控制保持系统性能的稳定,在此基础上利用预测控制的预测电流作为MTPA工作点计算的输入,简化了直接公式法(1)中的工作点q轴电流的求解。The formula solving method is a widely used method for solving the MTPA working point current. The existing MTPA formula methods mainly include the current solving method, the lead angle solving method, etc., as shown in Figure 2. The current solving method uses the speed loop to output the electromagnetic torque, and solves the d and q axis currents according to the formula for MTPA control. The formula is relatively complex to calculate; the lead angle current method is simpler than the direct formula method, but it is designed with trigonometric functions and the calculation is also relatively complex. The improved MTPA control algorithm proposed in this application first adopts predictive control to maintain the stability of system performance, and on this basis uses the predicted current of predictive control as the input for MTPA working point calculation, which simplifies the solution of the working point q axis current in the direct formula method (1).

如附图1所示,一种永磁同步电机MTPA预测控制方法,其特征在于:包括以下步骤:As shown in FIG1 , a permanent magnet synchronous motor MTPA predictive control method is characterized by comprising the following steps:

步骤一:在d-q旋转坐标系下建立永磁同步电机的基于MTPA的电流预测的多项目价值函数;所述多项目价值函数包括电流误差优化、最大转矩电流比优化以及电压误差优化;Step 1: Establish a multi-item value function for current prediction based on MTPA of a permanent magnet synchronous motor in a d-q rotating coordinate system; the multi-item value function includes current error optimization, maximum torque current ratio optimization, and voltage error optimization;

步骤二:通过价值函数求出当前k时刻电机控制最优电压,并根据该最优控制电压和永磁同步电机得电流预测模型预测下一个k+1时刻的d轴电流值;Step 2: Calculate the optimal voltage of the motor control at the current k moment through the value function, and predict the d-axis current value at the next k+1 moment based on the optimal control voltage and the current prediction model of the permanent magnet synchronous motor;

步骤三:将电机转速参考值和电机当前运行转速值输入PI调节器后获得该时刻的电磁转矩参考值Te *(k);Step 3: Input the motor speed reference value and the current running speed value of the motor into the PI regulator to obtain the electromagnetic torque reference value T e * (k) at that moment;

步骤四:将步骤三获得的电磁转矩参考值、根据在k-1时刻预测到的k时刻的d轴预测电流作为MTPA模块的输入,MTPA模块输出其工作时的q轴电流以及d轴电流。Step 4: The electromagnetic torque reference value obtained in step 3 and the predicted d-axis current at time k predicted at time k-1 are used as inputs of the MTPA module, and the MTPA module outputs the q-axis current and d-axis current during its operation.

进一步地,所述步骤一具体包括:Furthermore, the step 1 specifically includes:

S11:永磁同步电机电流方程为:S11: The current equation of the permanent magnet synchronous motor is:

Figure BDA0004127752170000061
Figure BDA0004127752170000061

Figure BDA0004127752170000062
Figure BDA0004127752170000062

公式(1)中id为d轴电流,i q为q轴电流;Rs为电阻;Ld为d轴磁导、Lq为q轴磁导;P为电机极对数;ud为d轴电压,uq为q轴电压;ωm为机械转速;φf为磁链;In formula (1), i d is the d-axis current, i q is the q-axis current; R s is the resistance; L d is the d-axis permeance, Lq is the q-axis permeance; P is the number of motor pole pairs; ud is the d-axis voltage, u q is the q-axis voltage; ω m is the mechanical speed; φ f is the flux linkage;

S12:将永磁同步电机电流方程离散化;得到永磁同步电机的电流离散模型如下式(2):S12: Discretize the permanent magnet synchronous motor current equation; obtain the current discrete model of the permanent magnet synchronous motor as follows (2):

Figure BDA0004127752170000063
Figure BDA0004127752170000063

Figure BDA0004127752170000064
Figure BDA0004127752170000064

(2)式中,T为采样周期,k表示k个采样时刻;(2) In the formula, T is the sampling period, k represents k sampling moments;

S13:将公式(2)改写为状态空间模型得到如下式(3)的状态空间方程:S13: Rewrite formula (2) into a state space model to obtain the state space equation of formula (3):

x(k+1)=Ax(k)+Bu(k)x(k+1)=Ax(k)+Bu(k)

y(k)=Cx(k) (3)y(k)=Cx(k) (3)

(3)式中,

Figure BDA0004127752170000065
(3) In the formula,
Figure BDA0004127752170000065

Figure BDA0004127752170000066
Figure BDA0004127752170000066

Figure BDA0004127752170000067
Figure BDA0004127752170000067

S14:永磁同步电机采用MTPA模式运行时轨迹为MTPA控制曲线,即d、q轴电流满足公式:

Figure BDA0004127752170000071
基于该MTPA控制曲线的原理建立如下式(4)的基于MTPA电流预测的多项目价值函数;其中多项目价值优化包括d、q轴电流误差优化,d、q轴电压误差优化以及MTPA曲线误差优化;所述d、q轴电流误差为经过MTPA模块得到的电机d、q轴参考电流与电机当前运行d、q轴电流误差,通过该项控制电机电流能够跟踪MTPA控制工作点的参考电流运行;所述d、q轴电压误差为永磁同步电机控制电压增量误差达到最小;所述MTPA曲线误差为保证电流控制按照预设的MTPA曲线运行;S14: When the permanent magnet synchronous motor is operated in MTPA mode, the trajectory is the MTPA control curve, that is, the d-axis and q-axis currents satisfy the formula:
Figure BDA0004127752170000071
Based on the principle of the MTPA control curve, a multi-project value function based on MTPA current prediction is established as shown in the following formula (4); wherein the multi-project value optimization includes d and q axis current error optimization, d and q axis voltage error optimization and MTPA curve error optimization; the d and q axis current errors are the motor d and q axis reference currents obtained by the MTPA module and the motor's current d and q axis current errors, and the motor current can track the reference current operation of the MTPA control working point through this control; the d and q axis voltage errors are the permanent magnet synchronous motor control voltage increment errors that are minimized; the MTPA curve error is to ensure that the current control runs according to the preset MTPA curve;

Figure BDA0004127752170000072
Figure BDA0004127752170000072

(4)式中,kid为当前d轴电流与d轴参考电流之

Figure BDA0004127752170000073
间的误差权重;kiq为q轴电流与q轴参考电流之
Figure BDA0004127752170000074
间的误差权重;kmtpa为当电机采用MTPA运行时d、q轴电流关系误差的权重;kud和kuq分别为d、q轴电压在k时刻增量的权重。(4) where kid is the difference between the current d-axis current and the d-axis reference current.
Figure BDA0004127752170000073
The error weight between q-axis current and q-axis reference current; kiq is the error between q-axis current and q-axis reference current
Figure BDA0004127752170000074
kmtpa is the weight of the error between the d-axis and q-axis currents when the motor adopts MTPA operation; kud and kuq are the weights of the increments of the d-axis and q-axis voltages at time k, respectively.

进一步地,步骤二中通过对价值函数分别进行对△ud和△uq求偏导数,即:Furthermore, in step 2, the partial derivatives of the value function are calculated for △u d and △u q , namely:

Figure BDA0004127752170000075
Figure BDA0004127752170000075

Figure BDA0004127752170000076
Figure BDA0004127752170000076

从而得到当前k时刻的最优控制电压增量△ud和△uq;然后根据公式(5)获得k时刻最优电压,并将求得的该电压用于永磁同步电机控制;Thus, the optimal control voltage increments △u d and △u q at the current time k are obtained; then the optimal voltage at the time k is obtained according to formula (5), and the obtained voltage is used for permanent magnet synchronous motor control;

ud(k)=ud(k-1)+△ud(k)u d (k)=u d (k-1)+△u d (k)

uq(k)=uq(k-1)+△uq(k) (5);u q (k)=u q (k-1)+Δu q (k) (5);

将上式中获取的d、q轴电压带入公式(6),得到k+1时刻的d轴电流idmiddle(k+1),作为MTPA工作点求解的中间变量:Substituting the d-axis and q-axis voltages obtained in the above formula into formula (6), we can obtain the d-axis current i dmiddle (k+1) at time k+1, which is used as the intermediate variable for solving the MTPA operating point:

Figure BDA0004127752170000077
Figure BDA0004127752170000077

进一步地,步骤三具体为:Furthermore, step three is specifically as follows:

永磁同步电机速度公式如下:The speed formula of permanent magnet synchronous motor is as follows:

Figure BDA0004127752170000081
Figure BDA0004127752170000081

利用上式公式,将电机速度参考值和当前时刻电机的速度反馈值输入PI调节器并得到电机的电磁转矩参考值Te *(k)。By using the above formula, the motor speed reference value and the speed feedback value of the motor at the current moment are input into the PI regulator to obtain the electromagnetic torque reference value Te * (k) of the motor.

进一步地,将步骤三获得的电磁转矩参考值Te *(k)、步骤一中的d轴预测电流idpre(k+1)作为MTPA算法模块的输入求解电机MTPA控制时的工作点,具体如下:Furthermore, the electromagnetic torque reference value Te * (k) obtained in step 3 and the d-axis predicted current i dpre (k+1) in step 1 are used as inputs of the MTPA algorithm module to solve the working point of the motor during MTPA control, as follows:

S41:永磁同步电机电磁转矩方程如(8)所示,S41: The electromagnetic torque equation of the permanent magnet synchronous motor is shown in (8),

Figure BDA0004127752170000082
Figure BDA0004127752170000082

将电磁转矩Te *(k)、预测电流idmiddle(k)代入公式(8),可以得到MTPA工作时的q轴电流:Substituting the electromagnetic torque Te * (k) and the predicted current i dmiddle (k) into formula (8), the q-axis current when MTPA is working can be obtained:

Figure BDA0004127752170000083
Figure BDA0004127752170000083

S42:由于当永磁同步电机采用MTPA模式运行时,其d、q轴电流满足下式(10):S42: When the permanent magnet synchronous motor is operated in MTPA mode, its d-axis and q-axis currents satisfy the following equation (10):

Figure BDA0004127752170000084
Figure BDA0004127752170000084

因此将得到的中iqmtpa(k)带入公式(10),即得到MTPA模板工作时的d轴工作电流,如下式(11):Therefore, the obtained i qmtpa (k) is substituted into formula (10), and the d-axis working current of the MTPA template when working is obtained, which is as follows (11):

Figure BDA0004127752170000085
Figure BDA0004127752170000085

虽然本发明已以较佳实施例公开如上,但它们并不是用来限定本发明的,任何熟习此技艺者,在不脱离本发明之精神和范围内,自当可作各种变化或润饰,因此本发明的保护范围应当以本申请的权利要求保护范围所界定的为准。Although the present invention has been disclosed as above in terms of preferred embodiments, they are not intended to limit the present invention. Anyone skilled in the art can make various changes or modifications without departing from the spirit and scope of the present invention. Therefore, the scope of protection of the present invention should be based on the scope of protection defined by the claims of this application.

Claims (5)

1.一种永磁同步电机MTPA预测控制方法,其特征在于:包括以下步骤:1. A permanent magnet synchronous motor MTPA predictive control method, characterized in that it comprises the following steps: 步骤一:在d-q旋转坐标系下建立永磁同步电机的基于MTPA的电流预测的多项目价值函数;所述多项目价值函数包括电流误差优化、最大转矩电流比优化以及电压误差优化;Step 1: Establish a multi-item value function for current prediction based on MTPA of a permanent magnet synchronous motor in a d-q rotating coordinate system; the multi-item value function includes current error optimization, maximum torque current ratio optimization, and voltage error optimization; 步骤二:通过价值函数求出当前k时刻电机控制最优电压,并根据该最优控制电压和永磁同步电机得电流预测模型预测下一个k+1时刻的d轴电流值;Step 2: Calculate the optimal voltage of the motor control at the current k moment through the value function, and predict the d-axis current value at the next k+1 moment based on the optimal control voltage and the current prediction model of the permanent magnet synchronous motor; 步骤三:将电机转速参考值和电机当前运行转速值输入PI调节器后获得该时刻的电磁转矩参考值Te *(k);Step 3: Input the motor speed reference value and the current running speed value of the motor into the PI regulator to obtain the electromagnetic torque reference value T e * (k) at that moment; 步骤四:将步骤三获得的电磁转矩参考值、根据在k-1时刻预测到的k时刻的d轴预测电流作为MTPA模块的输入,MTPA模块输出其工作时的q轴电流以及d轴电流。Step 4: The electromagnetic torque reference value obtained in step 3 and the predicted d-axis current at time k predicted at time k-1 are used as inputs of the MTPA module, and the MTPA module outputs the q-axis current and d-axis current during its operation. 2.根据权利要求1所述的一种永磁同步电机MTPA预测控制方法,其特征在于:所述步骤一具体包括:2. A permanent magnet synchronous motor MTPA predictive control method according to claim 1, characterized in that: the step 1 specifically comprises: S11:永磁同步电机电流方程为:S11: The current equation of the permanent magnet synchronous motor is:
Figure FDA0004127752160000011
Figure FDA0004127752160000011
Figure FDA0004127752160000012
Figure FDA0004127752160000012
公式(1)中id为d轴电流,iq为q轴电流;Rs为电阻;Ld为d轴磁导、Lq为q轴磁导;P为电机极对数;ud为d轴电压,uqq轴电压;ωm为机械转速;ff为磁链;In formula (1), i d is the d-axis current, i q is the q-axis current; R s is the resistance; L d is the d-axis permeance, Lq is the q-axis permeance; P is the number of motor pole pairs; ud is the d-axis voltage, u q is the q -axis voltage; ω m is the mechanical speed; f f is the magnetic flux; S12:将永磁同步电机电流方程离散化;得到永磁同步电机的电流离散模型如下式(2):S12: Discretize the permanent magnet synchronous motor current equation; obtain the current discrete model of the permanent magnet synchronous motor as follows (2):
Figure FDA0004127752160000013
Figure FDA0004127752160000013
(2)式中,T为采样周期,k表示k个采样时刻;(2) In the formula, T is the sampling period, k represents k sampling moments; S13:将公式(2)改写为状态空间模型得到如下式(3)的状态空间方程:S13: Rewrite formula (2) into a state space model to obtain the state space equation of formula (3):
Figure FDA0004127752160000014
Figure FDA0004127752160000014
(3)式中,
Figure FDA0004127752160000015
(3) In the formula,
Figure FDA0004127752160000015
Figure FDA0004127752160000021
Figure FDA0004127752160000021
Figure FDA0004127752160000022
Figure FDA0004127752160000022
S14:永磁同步电机采用MTPA模式运行时轨迹为MTPA控制曲线,即d、q轴电流满足公式:
Figure FDA0004127752160000023
基于该MTPA控制曲线的原理建立如下式(4)的基于MTPA电流预测的多项目价值函数;其中多项目价值优化包括d、q轴电流误差优化,d、q轴电压误差优化以及MTPA曲线误差优化;所述d、q轴电流误差为经过MTPA模块得到的电机d、q轴参考电流与电机当前运行d、q轴电流误差,通过该项控制电机电流能够跟踪MTPA控制工作点的参考电流运行;所述d、q轴电压误差为永磁同步电机控制电压增量误差达到最小;所述MTPA曲线误差为保证电流控制按照预设的MTPA曲线运行;
S14: When the permanent magnet synchronous motor is operated in MTPA mode, the trajectory is the MTPA control curve, that is, the d-axis and q-axis currents satisfy the formula:
Figure FDA0004127752160000023
Based on the principle of the MTPA control curve, a multi-project value function based on MTPA current prediction is established as shown in the following formula (4); wherein the multi-project value optimization includes d and q axis current error optimization, d and q axis voltage error optimization and MTPA curve error optimization; the d and q axis current errors are the motor d and q axis reference currents obtained by the MTPA module and the motor's current d and q axis current errors, and the motor current can track the reference current operation of the MTPA control working point through this control; the d and q axis voltage errors are the permanent magnet synchronous motor control voltage increment errors that are minimized; the MTPA curve error is to ensure that the current control runs according to the preset MTPA curve;
Figure FDA0004127752160000024
Figure FDA0004127752160000024
(4)式中,kid为当前d轴电流与d轴参考电流之
Figure FDA0004127752160000025
间的误差权重;kiq为q轴电流与q轴参考电流之
Figure FDA0004127752160000026
间的误差权重;kmtpa为当电机采用MTPA运行时d、q轴电流关系误差的权重;kud和kuq分别为d、q轴电压在k时刻增量的权重。
(4) where kid is the difference between the current d-axis current and the d-axis reference current.
Figure FDA0004127752160000025
The error weight between q-axis current and q-axis reference current; kiq is the error between q-axis current and q-axis reference current
Figure FDA0004127752160000026
kmtpa is the weight of the error between the d-axis and q-axis currents when the motor adopts MTPA operation; kud and kuq are the weights of the increments of the d-axis and q-axis voltages at time k, respectively.
3.根据权利要求2所述的一种永磁同步电机MTPA预测控制方法,其特征在于:步骤二中通过对价值函数分别进行对△ud和△uq求偏导数,即:3. A permanent magnet synchronous motor MTPA predictive control method according to claim 2, characterized in that: in step 2, partial derivatives of the cost function are obtained for Δu d and Δu q respectively, that is:
Figure FDA0004127752160000027
Figure FDA0004127752160000027
Figure FDA0004127752160000028
Figure FDA0004127752160000028
从而得到当前k时刻的最优控制电压增量△ud和△uq;然后根据公式(5)获得k时刻最优电压,并将求得的该电压用于永磁同步电机控制;Thus, the optimal control voltage increments △u d and △u q at the current time k are obtained; then the optimal voltage at the time k is obtained according to formula (5), and the obtained voltage is used for permanent magnet synchronous motor control;
Figure FDA0004127752160000031
Figure FDA0004127752160000031
将上式中获取的d、q轴电压带入公式(6),得到k+1时刻的d轴电流idmiddle(k+1),作为MTPA工作点求解的中间变量:Substituting the d-axis and q-axis voltages obtained in the above formula into formula (6), we can obtain the d-axis current i dmiddle (k+1) at time k+1, which is used as the intermediate variable for solving the MTPA operating point:
Figure FDA0004127752160000032
Figure FDA0004127752160000032
4.根据权利要求3所述的一种永磁同步电机MTPA预测控制方法,其特征在于:步骤三具体为:永磁同步电机速度公式如下:4. A permanent magnet synchronous motor MTPA predictive control method according to claim 3, characterized in that: step three is specifically: the permanent magnet synchronous motor speed formula is as follows:
Figure FDA0004127752160000033
Figure FDA0004127752160000033
利用上式公式,将电机速度参考值和当前时刻电机的速度反馈值输入PI调节器并得到电机的电磁转矩参考值
Figure FDA0004127752160000034
Using the above formula, the motor speed reference value and the current motor speed feedback value are input into the PI regulator to obtain the motor electromagnetic torque reference value.
Figure FDA0004127752160000034
5.根据权利要求4所述的一种永磁同步电机MTPA预测控制方法,其特征在于:将步骤三获得的电磁转矩参考值
Figure FDA0004127752160000035
步骤一中的d轴预测电流idpre(k+1)作为MTPA算法模块的输入求解电机MTPA控制时的工作点,具体如下:
5. A permanent magnet synchronous motor MTPA predictive control method according to claim 4, characterized in that: the electromagnetic torque reference value obtained in step 3 is
Figure FDA0004127752160000035
The d-axis predicted current i dpre (k+1) in step 1 is used as the input of the MTPA algorithm module to solve the working point of the motor during MTPA control, as follows:
S41:永磁同步电机电磁转矩方程如(8)所示,S41: The electromagnetic torque equation of the permanent magnet synchronous motor is shown in (8),
Figure FDA0004127752160000036
Figure FDA0004127752160000036
将电磁转矩
Figure FDA0004127752160000037
预测电流idmiddle(k)代入公式(8),可以得到MTPA工作时的q轴电流:
The electromagnetic torque
Figure FDA0004127752160000037
Substituting the predicted current i dmiddle (k) into formula (8), we can obtain the q-axis current when MTPA is working:
Figure FDA0004127752160000038
Figure FDA0004127752160000038
S42:由于当永磁同步电机采用MTPA模式运行时,其d、q轴电流满足下式(10):S42: When the permanent magnet synchronous motor is operated in MTPA mode, its d-axis and q-axis currents satisfy the following equation (10):
Figure FDA0004127752160000039
Figure FDA0004127752160000039
因此将得到的中iqmtpa(k)带入公式(10),即得到MTPA模板工作时的d轴工作电流,如下式(11):Therefore, the obtained i qmtpa (k) is substituted into formula (10), and the d-axis working current of the MTPA template when working is obtained, which is as follows (11):
Figure FDA00041277521600000310
Figure FDA00041277521600000310
CN202310250531.5A 2023-03-15 2023-03-15 A MTPA predictive control method for permanent magnet synchronous motor Withdrawn CN116094400A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310250531.5A CN116094400A (en) 2023-03-15 2023-03-15 A MTPA predictive control method for permanent magnet synchronous motor

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310250531.5A CN116094400A (en) 2023-03-15 2023-03-15 A MTPA predictive control method for permanent magnet synchronous motor

Publications (1)

Publication Number Publication Date
CN116094400A true CN116094400A (en) 2023-05-09

Family

ID=86214239

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310250531.5A Withdrawn CN116094400A (en) 2023-03-15 2023-03-15 A MTPA predictive control method for permanent magnet synchronous motor

Country Status (1)

Country Link
CN (1) CN116094400A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117978030A (en) * 2024-02-05 2024-05-03 南京工业大学 A method and system for predictive control of MTPA of permanent magnet synchronous motor with constant parameters

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117978030A (en) * 2024-02-05 2024-05-03 南京工业大学 A method and system for predictive control of MTPA of permanent magnet synchronous motor with constant parameters
CN117978030B (en) * 2024-02-05 2025-03-07 南京工业大学 Permanent magnet synchronous motor constant parameter MTPA predictive control method and system

Similar Documents

Publication Publication Date Title
CN110336501B (en) A Model Predictive Control Method for Embedded Permanent Magnet Synchronous Motors
CN112701968B (en) A Robust Performance Improvement Method for Model Predictive Control of Permanent Magnet Synchronous Motors
CN107317532B (en) Predictive current control method and system for permanent magnet synchronous motor based on sliding mode
CN106788075B (en) Rapid vector screening prediction method for controlling torque based on improved Euler method
CN108092567B (en) A speed control system and method for a permanent magnet synchronous motor
CN106788045B (en) A kind of permanent magnet synchronous motor model prediction PI changeable weight concurrency control method
CN106655938B (en) Control system for permanent-magnet synchronous motor and control method based on High-Order Sliding Mode method
CN112422004B (en) Disturbance suppression method for permanent magnet synchronous motor in weak magnetic control mode
CN110829902B (en) Model-free predictive speed control method and system for permanent magnet synchronous motor for electric vehicle
CN111987943B (en) A Robust Performance Improvement Method for Model Predictive Control of Permanent Magnet Synchronous Motors
CN111555684B (en) Multi-step Model Predictive Torque Control Method for Permanent Magnet Synchronous Motors with Variable Switch Points
CN107154763A (en) Permagnetic synchronous motor dead beat direct Torque Control and control method
CN114679095A (en) A finite set model predictive current control method for permanent magnet motors based on disturbance compensation
CN111092583A (en) A three-phase permanent magnet synchronous motor drive system current loop delay compensation method
CN116961512B (en) Model prediction-based current control method, device and storage medium
CN110061671A (en) Permanent magnet synchronous motor control method and system based on speed change approach rate
CN112701970B (en) PMSM (permanent magnet synchronous motor) accurate discrete active-disturbance-rejection control method under low carrier ratio
CN110649854A (en) A Continuous Set Model Predictive Torque Control Method for Surface Mount Permanent Magnet Synchronous Motors
CN116094400A (en) A MTPA predictive control method for permanent magnet synchronous motor
Liu et al. Speed estimation and Parameters Identification simultaneously of PMSM based on MRAS
CN114157193B (en) A search and interpolation synchronous motor torque ripple suppression control method and system
CN117724342A (en) Combined control system of solar wing driving mechanism
CN113824364A (en) Method, system, equipment and medium for flux weakening control of permanent magnet synchronous motor
CN107800334B (en) PMSM pre-synchronization coaxial operation control method and system
Shi et al. Improved deadbeat predictive current control for SPMSM based on multi-step 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
WW01 Invention patent application withdrawn after publication
WW01 Invention patent application withdrawn after publication

Application publication date: 20230509