CN113644852A - 用于永磁同步电动机的鲁棒三矢量模型预测磁链控制方法 - Google Patents
用于永磁同步电动机的鲁棒三矢量模型预测磁链控制方法 Download PDFInfo
- Publication number
- CN113644852A CN113644852A CN202110565963.6A CN202110565963A CN113644852A CN 113644852 A CN113644852 A CN 113644852A CN 202110565963 A CN202110565963 A CN 202110565963A CN 113644852 A CN113644852 A CN 113644852A
- Authority
- CN
- China
- Prior art keywords
- flux linkage
- permanent magnet
- magnet synchronous
- synchronous motor
- psi
- 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.)
- Granted
Links
- 230000004907 flux Effects 0.000 title claims abstract description 164
- 230000001360 synchronised effect Effects 0.000 title claims abstract description 74
- 239000013598 vector Substances 0.000 title claims abstract description 63
- 238000000034 method Methods 0.000 title claims abstract description 47
- 238000013178 mathematical model Methods 0.000 claims abstract description 11
- 230000009471 action Effects 0.000 claims description 19
- 230000008569 process Effects 0.000 claims description 13
- 238000005070 sampling Methods 0.000 claims description 4
- 230000008859 change Effects 0.000 abstract description 8
- 230000035945 sensitivity Effects 0.000 abstract description 2
- 230000006870 function Effects 0.000 description 16
- 238000004364 calculation method Methods 0.000 description 9
- 238000006243 chemical reaction Methods 0.000 description 9
- 238000001514 detection method Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 6
- 238000005457 optimization Methods 0.000 description 5
- 238000004088 simulation Methods 0.000 description 4
- 230000015572 biosynthetic process Effects 0.000 description 3
- 238000011217 control strategy Methods 0.000 description 3
- 238000003786 synthesis reaction Methods 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 238000009795 derivation Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 229910052761 rare earth metal Inorganic materials 0.000 description 1
- 150000002910 rare earth metals Chemical class 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P21/00—Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
- H02P21/14—Estimation or adaptation of machine parameters, e.g. flux, current or voltage
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P21/00—Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
- H02P21/06—Rotor flux based control involving the use of rotor position or rotor speed sensors
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P21/00—Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
- H02P21/13—Observer control, e.g. using Luenberger observers or Kalman filters
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P25/00—Arrangements or methods for the control of AC motors characterised by the kind of AC motor or by structural details
- H02P25/02—Arrangements 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/022—Synchronous motors
- H02P25/024—Synchronous motors controlled by supply frequency
- H02P25/026—Synchronous motors controlled by supply frequency thereby detecting the rotor position
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P27/00—Arrangements or methods for the control of AC motors characterised by the kind of supply voltage
- H02P27/04—Arrangements 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/06—Arrangements 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/08—Arrangements 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
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P6/00—Arrangements for controlling synchronous motors or other dynamo-electric motors using electronic commutation dependent on the rotor position; Electronic commutators therefor
- H02P6/34—Modelling or simulation for control purposes
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P2207/00—Indexing scheme relating to controlling arrangements characterised by the type of motor
- H02P2207/05—Synchronous 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
本发明公开了用于永磁同步电动机的鲁棒三矢量模型预测磁链控制方法,具体按照以下步骤实施:步骤1、建立永磁同步电动机的数学模型;步骤2、选取一拍延时补偿控制下新的磁链控制变量,计算新的磁链控制变量控制条件下永磁同步电动机的磁链预测方程;步骤3、分析参数扰动及未建模动态扰动对永磁同步电动机的磁链预测方程准确性影响;步骤4、根据步骤3得到的准确性影响,设计龙贝格扰动观测器以及龙贝格扰动观测器稳定条件;步骤5、将龙贝格观测的集总扰动补偿到永磁同步电动机的磁链预测方程中,实现对永磁同步电动机的准确控制;解决了模型预测磁链控制对参数变化敏感及未考虑未建模动态扰动,带来的鲁棒性问题。
Description
技术领域
本发明属于电机控制技术领域,具体涉及用于永磁同步电动机的鲁棒三 矢量模型预测磁链控制方法。
背景技术
永磁同步电动机(Permanent Magnet Synchronous Motor,PMSM)利用 稀土永磁体励磁,具有结构简单、可靠性高、效率高、功率因数高等优点, 广泛应用于交流电机传动领域,如工业机器人、航空航天、电动汽车等领域。 随着上述领域的不断发展,对永磁同步电动机的控制也有了更高的要求。
模型预测控制是20世纪80年代初次应用在功率变换器和传动装置的新 兴控制策略,具有概念简单、可考虑多变量、易于包含约束条件和非线性因 素等优点。其主要思想是根据***当前时刻状态值,依据被控对象数学模型 预测出***下一时刻的状态,通过构建的价值函数进行在线寻优,确定出最 优控制序列,并将这个控制序列作用到控制***中。模型预测控制根据控制 目标的不同可分为模型预测电流控制、模型预测速度控制、模型预测转矩控 制。模型预测磁链控制是针对模型预测转矩中权重系数难以整定的一种改进 方法,由转矩、磁链和负载角的关系推导得到等效定子磁链矢量,以等效的 定子磁链矢量作为控制目标,从而无需权重系数。然而在模型预测磁链控制 方法中,等效定子磁链的给定值为永磁体磁链值,负载角的计算含有电机参 数,预测方程也基于***模型,因此对电机参数依赖性较大。在实际应用中, 电机参数实际值会因为温度变化、磁路饱和、***测量不准确等因素与电机 铭牌值不一致,并且模型预测磁链控制未考虑***的未建模动态扰动,从而 引起磁链给定值与实际值之间的稳态误差,最终会影响***的性能,导致系 统鲁棒性低。
目前,为了提高预测控制鲁棒性的方法主要包括:参数辨识、加入预测 误差矫正项、采用具有鲁棒性强的控制器、超局部无模型预测控制和采用观 测器观测集总扰动进行补偿等方法。其中采用观测器观测集总扰动可以实现 对参数扰动和未建模动态扰动的集总补偿,并且具有原理实现简单等优点, 从而得到了广泛的应用。
由于模型预测磁链控制不仅在预测模型中含有电机参数,在其给定值和 负载角计算中也含有电机参数,并且预测过程中没有考虑未建模的动态扰动 因此有必要进行模型预测磁链控制鲁棒性提升的研究。
发明内容
本发明的目的是提供用于永磁同步电动机的鲁棒三矢量模型预测磁链 控制方法,解决了模型预测磁链控制对参数变化敏感及未考虑未建模动态扰 动,带来的鲁棒性问题。
本发明所采用的技术方案是,用于永磁同步电动机的鲁棒三矢量模型预 测磁链控制方法,具体按照以下步骤实施:
步骤1、建立永磁同步电动机的数学模型;
步骤2、选取一拍延时补偿控制下新的磁链控制变量,计算新的磁链控 制变量控制条件下永磁同步电动机的磁链预测方程;
步骤3、分析参数扰动及未建模动态扰动对永磁同步电动机的磁链预测 方程准确性影响;
步骤4、根据步骤3得到的准确性影响,设计龙贝格扰动观测器以及龙 贝格扰动观测器稳定条件;
步骤5、将龙贝格观测的集总扰动补偿到永磁同步电动机的磁链预测方 程中,实现对永磁同步电动机的准确控制。
本发明的特点还在于:
步骤1中永磁同步电动机的数学模型为:
ψd=Ldid+ψf (5)
ψq=Lqiq (6)
式(1)-(6)中,ud为定子电压的直轴分量,uq为定子电压的交轴分量;Rs为定子电阻;id为定子电流的直轴分量,iq为定子电流的交轴分量;ωre是转 子电角速度;ψd为定子磁链的直轴分量,ψq为定子磁链的交轴分量;Ld为 定子电感的直轴分量,Lq为定子电感的交轴分量;ψf是永磁体磁链;Te为 电磁转矩;pn为电机极对数;J为转动惯量;B为摩擦粘滞系数;TL为负载 转矩;t为时间;其中,表贴式永磁同步电动机中Ld=Lq=Ls。
步骤2具体过程为:
根据永磁同步电动机的数学模型,令ψds(k)=Ldid(k),采用后向欧拉法 离散化可得:
式中:ψd(k+1)为(k+1)时刻直轴磁链的预测值;ψq(k)为k时刻交轴磁链的 预测值;ud(k)为定子电压在k时刻的直轴分量;
将式(7)中对直轴磁链ψd的控制等效为对变量ψds的控制,采用新的控 制变量ψds、ψq对永磁同步电机的模型预测磁链控制,重新构造价值函数为:
C=|ψds *-ψds(k+1)|2+|ψq *-ψq(k+1)|2 (8)
式中:ψds *为变量ψds的给定值;
对永磁同步电机控制采用一拍延时补偿,则式(8)中的ψds(k+1)和 ψq(k+1)替换为定子磁链在(k+2)时刻的直交轴分量ψds(k+2)和交轴分量 ψq(k+2),则修正后的价值函数表示为:
C=|ψds *-ψds(k+2)|2+|ψq *-ψq(k+2)|2 (9)
新的磁链控制变量控制条件下永磁同步电动机的磁链预测方程为:
式中:ud(k+1)为定子电压在(k+1)时刻的直轴分量,uq(k+1)为定子电 压在(k+1)时刻的交轴分量。
步骤3具体过程为:
将电机参数实际值与该电机铭牌标称值无偏差时的预测磁链表示为ψ′ds、 ψ′q,将电机参数实际值与铭牌标称值有偏差时的预测磁链表示为ψds、ψq, 用Δψf表示永磁体磁链实际值与铭牌标称值之间的误差,ΔLs表示定子电感 实际值与铭牌标称值之间的误差,ΔRs表示定子电阻实际值与铭牌标称值之 间的误差,用ε表示未建模动态的不确定扰动;
考虑电机参数实际值与铭牌标称值存在偏差以及未建模扰动,预测磁链 表示为:
根据式(13)获得参数扰动及未建模动态对永磁同步电动机的磁链预测 方程准确性影响。
步骤4具体过程为:
根据式(13)可得:
式中:fd直轴预测磁链的总扰动,fq是交轴预测磁链的总扰动;
离散后的龙贝格观测器为:
忽略-ωreTsj项,根据式(18)可得龙贝格观测器的特征方程为:
即得:
D(z)=c0+c1z1+c2z2=0 (21)
根据朱利判据可知,观测器稳定的充要条件是:
将式(21)中的c0、c1、c2代入式(22)得:
根据式(23)得观测器稳定需要满足的最终条件为:
步骤5具体过程为:
加入龙贝格观测器估计的扰动量后,根据式(14)并考虑一拍延迟补偿, 得到鲁棒三矢量模型预测磁链k+2时刻的表达式:
将鲁棒三矢量模型预测磁链k+2时刻的表达式带入修正后的价值函数, 选取修正后的价值函数取值最小时,对应的一组电压矢量及其作用时间,作 用于逆变器,实现对永磁同步电动机的准确控制。
本发明的有益效果是:
(1)本发明是的控制方法摒弃了模型预测磁链控制通过分析参考磁链 幅值、参考转矩和负载角之间的数学关系,得到等效的参考磁链矢量的方法, 重新选取了磁链控制变量,实现永磁同步电动机的模型预测磁链控制,可使 磁链的给定值不包含电机参数,且无需通过含有电机参数的等效负载角进行 转矩和磁链的转换,具有算法简单,参数鲁棒性强的优点。
(2)本发明通过分析参数变化以及不确定动态扰动对预测磁链方程的 影响,根据预测磁链方程,构造了龙贝格扰动观测器,用来观测***包括电 机参数变化在内的集总扰动。在电机参数发生变化或者存在外部扰动时,使 得磁链和电流的实际值可以更加准确地跟随给定值,以减小磁链和电流的稳 态误差,进一步提高了***的鲁棒性。
附图说明
图1是本发明的控制方法应用的***原理图;
图2是图1的***在定子电阻、定子磁链及定子电感同时变化时的转速 波形图;
图3是图1的***在定子电阻、定子磁链及定子电感同时变化时的ψds波 形图;
图4是图1的***在定子电阻、定子磁链及定子电感同时变化时的ψq波 形图;
图5是图1的***在定子电阻、定子磁链及定子电感同时变化时的a相 电流波形图;
图中,1.三相逆变器,2.电流检测电路,3.永磁同步电动机,4.旋转编码 器,5.Clark变换模块,6.Park变换模块,7.k时刻磁链值估计模块,8.龙贝格 扰动观测器模块,9.一拍延迟补偿模块,10.k+2时刻磁链预测模块,11.作用 时间计算模块,12.期望矢量合成模块,13.价值函数优化模块。
具体实施方式
下面结合附图及具体实施方式对本发明进行详细说明。
本发明方法主要用于永磁同步电动机,其控制***如图1所示,包括信 号检测电路、主电路和控制电路;主电路包括三相逆变器1,主要用于驱动 永磁同步电动机3;信号检测电路包括电流检测电路2和旋转编码器4,主 要用于检测永磁同步电动机3的电流和转子位置信号;控制电路包括Clark 变换模块5、Park变换模块6、k时刻磁链值估计模块7、龙贝格扰动观测器 模块8、一拍延迟补偿模块9、k+2时刻磁链预测模块10、作用时间计算模 块11、期望矢量合成模块12、价值函数优化模块13,主要用于对信号检测 电路得到的信号就行处理,得到控制主电路的控制信号。
其中,控制电路将电流检测电路2检测永磁同步电动机3后得到的三相 电流ia、ib和ic经过Clark变换模块5处理后得到iα和iβ;旋转编码器4检测 永磁同步电动机3的转子位置角θ,求导得到电机机械角速度ωm;iα、iβ与θ 经过Park变换模块6处理后得到两相旋转坐标系下的直轴和交轴电流反馈 值id和iq;id和iq经过k时刻磁链值估计模块7得到k时刻的磁链值ψds(k)、 ψq(k);k时刻的磁链值ψds(k)、ψq(k)经过龙贝格扰动观测器模块8得到估 计的集总扰动电机机械角速度ωm、估计的集总扰动 和k时刻的磁链值ψds(k)、ψq(k)经过一拍延迟补偿模块9得到k+1时 刻的磁链预测值ψds(k+1)、ψq(k+1);k+1时刻的磁链预测值ψds(k+1)、ψq(k+1) 和估计的集总扰动经过k+2时刻磁链预测模块10得到k+2时 刻的磁链预测值ψds(k+2)、ψq(k+2);电机机械角速度给定值和电机机械 角速度ωm作差后经过PI控制器处理后得到q轴磁链给定值由估计的集总扰动及k+1时刻的磁链预测值ψds(k+1)、ψq(k+1)经过作用时间计算模块11得到两个有效矢量和一个零矢量的作用时间ti,j,z; 矢量作用时间ti,j,z经过期望矢量合成模块12可以得到6组期望的电压矢量 uevvⅠ~Ⅵ;将磁链给定值期望电压矢量uevvⅠ~Ⅵ及k+2时刻磁链预测值 ψds(k+2)、ψq(k+2)经过价值函数优化模块13,选择出最佳的一组电压矢量 及作用时间作用于逆变器用以控制永磁同步电动机3。
由此可知,通过获得最佳的一组电压矢量及作用时间便可提高永磁同步 电动机的鲁棒性。
本发明用于永磁同步电动机的鲁棒三矢量模型预测磁链控制方法,具体 按照以下步骤实施:
步骤1、建立永磁同步电动机的数学模型;
永磁同步电动机的数学模型为:
ψd=Ldid+ψf (5)
ψq=Lqiq (6)
式(1)-(6)中,ud为定子电压的直轴分量,uq为定子电压的交轴分 量;Rs为定子电阻;id为定子电流的直轴分量,iq为定子电流的交轴分量;ωre是转子电角速度;ψd为定子磁链的直轴分量,ψq为定子磁链的交轴分量;Ld为定子电感的直轴分量,Lq为定子电感的交轴分量;ψf是永磁体磁链;Te为 电磁转矩;pn为电机极对数;J为转动惯量;B为摩擦粘滞系数;TL为负载 转矩;t为时间;其中,表贴式永磁同步电动机中Ld=Lq=Ls。
步骤2、选取一拍延时补偿控制下新的磁链控制变量,计算新的磁链控 制变量控制条件下永磁同步电动机的磁链预测方程;
步骤2具体过程为:
根据永磁同步电动机的数学模型,选取新的控制变量进行模型预测磁链 控制:由式(3)可以看出交轴磁链与转矩有关,因此可由转速环PI调节器 直接得到。令ψds(k)=Ldid(k),采用后向欧拉法离散化可得:
式中:ψd(k+1)为(k+1)时刻直轴磁链的预测值;ψq(k)为k时刻交轴磁 链的预测值;ud(k)为定子电压在k时刻的直轴分量;
将式(7)中对直轴磁链ψd的控制等效为对变量ψds的控制,采用新的控 制变量ψds、ψq对永磁同步电机的模型预测磁链控制,重新构造价值函数为:
C=|ψds *-ψds(k+1)|2+|ψq *-ψq(k+1)|2 (8)
式中:式中:ψds *为变量ψds的给定值,根据id=0控制,ψds *=0。因此可以 避免由于永磁体磁链值无法准确获取而带来的直轴磁链给定值不准确,无法 满足控制***要求的问题。
由于数字控制***存在保持和量化等环节,存在数字延时问题,对永磁 同步电机控制采用一拍延时补偿,经过一拍延迟补偿,即将式(8)所示的 价值函数中的ψds(k+1)和ψq(k+1)替换为定子磁链在(k+2)时刻的直交轴分 量ψds(k+2)和ψq(k+2),则修正后的价值函数表示为:
C=|ψds *-ψds(k+2)|2+|ψq *-ψq(k+2)|2 (9)
新的磁链控制变量控制条件下永磁同步电动机的磁链预测方程为:
式中:ud(k+1)为定子电压在(k+1)时刻的直轴分量,uq(k+1)为定子电压 在(k+1)时刻的交轴分量。
步骤3、分析参数扰动及未建模动态扰动对永磁同步电动机的磁链预测 方程准确性影响;步骤3具体过程为:
为了便于分析参数等各种扰动带来的预测误差,将电机参数实际值与该 电机铭牌标称值无偏差时的预测磁链表示为ψ′ds、ψ′q,将电机参数实际值与 铭牌标称值有偏差时的预测磁链表示为ψds、ψq,用Δψf表示永磁体磁链实 际值与铭牌标称值之间的误差,ΔLs表示定子电感实际值与铭牌标称值之间 的误差,ΔRs表示定子电阻实际值与铭牌标称值之间的误差,用ε表示未建模 动态的不确定扰动;
理想状态下,电机参数实际值就是铭牌标称值,预测磁链如式(12)所 示,
考虑电机参数实际值与铭牌标称值存在偏差以及未建模扰动,预测磁链 表示为:
由于温度和磁饱和等因素会使得***中存在参数不匹配现象,由式(12) 和式(13)可知,定子电阻、定子电感及磁链的不准确和***存在的未建模 动态不确定扰动,会影响磁链预测的准确性。
根据式(13)获得参数扰动及未建模动态对永磁同步电动机的磁链预测 方程准确性影响。
步骤4、根据步骤3得到的准确性影响,设计龙贝格扰动观测器以及龙 贝格扰动观测器稳定条件;步骤4具体过程为:
采用龙贝格观测器观测参数不匹配引起的扰动和未建模动态不确定扰 动,并在预测模型中进行补偿,以减小三矢量模型预测磁链控制策略中控制 变量的稳态误差。考虑到参数不匹配和未建模动态不确定扰动ε带来的影响, 式(13)可写作:
式中:fd直轴预测磁链的总扰动,fq是交轴预测磁链的总扰动;包括定子电 感、永磁体磁链和定子电阻不匹配引起的扰动以及未建模动态扰动。
离散后的龙贝格观测器为:
u(k)=ud(k)+j(uq(k)-ωreψf) (19)
龙贝格观测器的稳定性证明如下:朱利判据是现代控制***中可以在Z 域内应用的稳定判据,即可以根据***在Z域的特征值是否在单位圆内,从 而判断该***能否稳定。因此,先求取式(18)的特征值,实际控制***中, 由于Ts较小,忽略-ωreTsj项,根据式(18)可得龙贝格观测器的特征方程为:
即得:
D(z)=c0+c1z1+c2z2=0 (21)
根据朱利判据可知,观测器稳定的充要条件是:
将式(21)中的c0、c1、c2代入式(22)得:
根据式(23)得观测器稳定需要满足的最终条件为:
步骤5、将龙贝格观测的集总扰动补偿到永磁同步电动机的磁链预测方 程中,实现对永磁同步电动机的准确控制;
步骤5具体过程为:
加入龙贝格观测器估计的扰动量后,根据式(14)并考虑一拍延迟补偿, 得到鲁棒三矢量模型预测磁链k+2时刻的表达式:
将鲁棒三矢量模型预测磁链k+2时刻的表达式带入修正后的价值函数, 选取修正后的价值函数取值最小时,对应的一组电压矢量及其作用时间,作 用于逆变器,实现对永磁同步电动机的准确控制。
其作用时间计算采用现有的方法,具体为:
三矢量模型预测磁链控制策略中三个电压矢量的作用时间是根据直交轴 磁链无差拍原理计算的。推导出鲁棒三矢量模型预测磁链控制零矢量作用时 直交轴磁链斜率sd0和sq0的计算公式,如式(26)和(27)所示:
因此两个相邻有效电压矢量ui、uj作用时直交轴磁链斜率sdi、sqi、sdj和 sqj公式分别为:
式中:udi、uqi表示有效矢量ui在直、交轴的电压分量;
udj、uqj表示有效矢量uj在直、交轴的电压分量。
令定子磁链的预测值在一个控制周期结束时等于等效参考定子磁链矢量, 故直、交轴定子磁链预测公式可改写为:
ψds(k+2)=ψds(k+1)+sditi+sdjtj+sd0tz=ψds * (32)
ψq(k+2)=ψq(k+1)+sqiti+sqjtj+sq0tz=ψq * (33)
式中:ti表示有效矢量ui的作用时间;
tj表示有效矢量uj的作用时间;
tz表示零矢量的作用时间。
三个电压矢量的作用时间之和刚好等于一个控制周期Ts,即
Ts=ti+tj+tz (34)
将式(26)~(34)联立可解得:
tz=Ts-ti-tj (37)
由于两电平电压源逆变器共有6个有效电压矢量和2个零矢量,因此根 据步骤6计算得到的矢量作用时间,由两个相邻的基本电压矢量和零矢量共 可以合成6个幅值和方向都可变的有效矢量如式(38)、(39)所示:
将式(38)、(39)合成的电压矢量作用于修正后的价值函数(11)中, 选择出使得价值函数值最小的一组电压矢量作用于逆变器,即可控制永磁同 步电机。
实施例
本发明提供了永磁同步电动机一种鲁棒三矢量模型预测磁链控制方法, 该方法首先重新选取了模型预测磁链控制的磁链变量,使磁链的给定值不受 永磁体磁链值的影响,无须电机参数的负载角计算,具有计算简单,且鲁棒 性强的优点;同时分析了模型预测磁链控制由于参数扰动和未建模动态扰动 造成的预测误差,研究了龙贝格扰动观测器观测三矢量模型预测磁链控制中 存在的集总扰动,提高了***的鲁棒性。为了验证本发明方法的有效性,利 用MATLAB/SIMULINK进行仿真验证。
仿真模型中Clark变换模块5、Park变换模块6、k时刻磁链值估计模块 7、龙贝格扰动观测器模块8、一拍延迟补偿模块9、k+2时刻磁链预测模块 10、作用时间计算模块11、期望矢量合成模块12、价值函数优化模块13均 采用S-function Builder功能模块和C语言编程实现,采样频率均为10kHz。 在上述仿真模型中,永磁同步电动机参数设置为:永磁体磁链为0.303Wb、 定子电感为11.57mH、额定电压为380V、额定电流为4.4A、定子电阻为 1.64Ω、额定转速为2430(r/min)、极对数为4、转子惯量为0.0011(kg·m2)、 阻尼系数为0.001、额定负载转矩为9.6(N·m);仿真模型参数设置为:PI例 系数为0.001、积分系数为0.00003、龙贝格扰动观测器系数kT1=0.286、 kT2=-58.6。
图2为电机永磁体磁链变为120%ψf,定子电阻变为120%Rs,定子电感 变为80%Ls,空载起动至1000r/min;0.1s转速给定由1000r/min突变为 1500r/min;0.2s突加9.6N·m负载;0.3s负载转矩由9.6N·m突减至空载的 转速响应波形图,图3为ψds波形图;图4为ψq波形图;图5为a相电流 波形图;
从图2可以看出,本发明控制方法控制永磁同步电动机控制***上,在 电机的磁链、定子电阻及定子电感实际值与铭牌值不一致时,其转速能够准 确地跟随给定,转速不受电机参数变化及未建模动态的影响。
从图3和图4可以看出,本发明方法控制在电机的磁链、定子电阻及定 子电感实际值与铭牌值不一致时,ψds和ψq的实际值仍然可以准确地跟随给 定值,不存在由于参数变化和未建模动态造成的稳态误差,能够实现更加准 确地控制,提高了***的鲁棒性。
从图5可以看出本发明永磁同步电动机一种鲁棒三矢量模型预测磁链控 制在参数变化时相电流仍可以保持较好地正弦度。
上述实施例进一步说明本发明永磁同步电动机一种鲁棒三矢量模型预 测磁链控制方法的有效性。因此本发明用于永磁同步电动机的鲁棒三矢量模 型预测磁链控制方法可以使得电机在参数不匹配以及存在为建模动态扰动 时仍具有准确地控制,提高了模型预测磁链的鲁棒性。
Claims (6)
1.用于永磁同步电动机的鲁棒三矢量模型预测磁链控制方法,其特征在于,具体按照以下步骤实施:
步骤1、建立永磁同步电动机的数学模型;
步骤2、选取一拍延时补偿控制下新的磁链控制变量,计算新的磁链控制变量控制条件下永磁同步电动机的磁链预测方程;
步骤3、分析参数扰动及未建模动态扰动对永磁同步电动机的磁链预测方程准确性影响;
步骤4、根据步骤3得到的准确性影响,设计龙贝格扰动观测器以及龙贝格扰动观测器稳定条件;
步骤5、将龙贝格观测的集总扰动补偿到永磁同步电动机的磁链预测方程中,实现对永磁同步电动机的准确控制。
2.根据权利要求1所述用于永磁同步电动机的鲁棒三矢量模型预测磁链控制方法,其特征在于,步骤1中所述永磁同步电动机的数学模型为:
ψd=Ldid+ψf (5)
ψq=Lqiq (6)
式(1)-(6)中,ud为定子电压的直轴分量,uq为定子电压的交轴分量;Rs为定子电阻;id为定子电流的直轴分量,iq为定子电流的交轴分量;ωre是转子电角速度;ψd为定子磁链的直轴分量,ψq为定子磁链的交轴分量;Ld为定子电感的直轴分量,Lq为定子电感的交轴分量;ψf是永磁体磁链;Te为电磁转矩;pn为电机极对数;J为转动惯量;B为摩擦粘滞系数;TL为负载转矩;t为时间;其中,表贴式永磁同步电动机中Ld=Lq=Ls。
3.根据权利要求2所述用于永磁同步电动机的鲁棒三矢量模型预测磁链控制方法,其特征在于,步骤2具体过程为:
根据永磁同步电动机的数学模型,令ψds(k)=Ldid(k),采用后向欧拉法离散化可得:
式中:ψd(k+1)为(k+1)时刻直轴磁链的预测值;ψq(k)为k时刻交轴磁链的预测值;ud(k)为定子电压在k时刻的直轴分量;
将式(7)中对直轴磁链ψd的控制等效为对变量ψds的控制,采用新的控制变量ψds、ψq对永磁同步电机的模型预测磁链控制,重新构造价值函数为:
C=|ψds *-ψds(k+1)|2+|ψq *-ψq(k+1)|2 (8)
式中:ψds *为变量ψds的给定值;
对永磁同步电机控制采用一拍延时补偿,则式(8)中的ψds(k+1)和ψq(k+1)替换为定子磁链在(k+2)时刻的直交轴分量ψds(k+2)和交轴分量ψq(k+2),则修正后的价值函数表示为:
C=|ψds *-ψds(k+2)|2+|ψq *-ψq(k+2)|2 (9)
新的磁链控制变量控制条件下永磁同步电动机的磁链预测方程为:
式中:ud(k+1)为定子电压在(k+1)时刻的直轴分量,uq(k+1)为定子电压在(k+1)时刻的交轴分量。
5.根据权利要求4所述用于永磁同步电动机的鲁棒三矢量模型预测磁链控制方法,其特征在于,步骤4具体过程为:
根据式(13)可得:
式中:fd直轴预测磁链的总扰动,fq是交轴预测磁链的总扰动;
根据式(14),以ψds,ψq为状态变量,设计龙贝格观测器为:
离散后的龙贝格观测器为:
u(k)=ud(k)+j(uq(k)-ωreψf) (19)
忽略-ωreTsj项,根据式(18)可得龙贝格观测器的特征方程为:
即得:
D(z)=c0+c1z1+c2z2=0 (21)
根据朱利判据可知,观测器稳定的充要条件是:
将式(21)中的c0、c1、c2代入式(22)得:
根据式(23)得观测器稳定需要满足的最终条件为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110565963.6A CN113644852B (zh) | 2021-05-24 | 2021-05-24 | 用于永磁同步电动机的鲁棒三矢量模型预测磁链控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110565963.6A CN113644852B (zh) | 2021-05-24 | 2021-05-24 | 用于永磁同步电动机的鲁棒三矢量模型预测磁链控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113644852A true CN113644852A (zh) | 2021-11-12 |
CN113644852B CN113644852B (zh) | 2024-02-06 |
Family
ID=78415809
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110565963.6A Active CN113644852B (zh) | 2021-05-24 | 2021-05-24 | 用于永磁同步电动机的鲁棒三矢量模型预测磁链控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113644852B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114614712A (zh) * | 2022-03-27 | 2022-06-10 | 西北工业大学 | 基于模型预测-观测器的飞机电推进***扰动抑制方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110299882A (zh) * | 2019-06-21 | 2019-10-01 | 浙江大学 | 基于混合供电型开绕组永磁同步电机的快速三矢量模型预测控制方法 |
US10819264B1 (en) * | 2020-08-04 | 2020-10-27 | Wolong Electric Group Co. Ltd. | Robust starting system and method for interior permanent magnet synchronous motor control |
CN111987943A (zh) * | 2020-07-23 | 2020-11-24 | 西安理工大学 | 一种永磁同步电机模型预测控制的鲁棒性能提升方法 |
CN112737452A (zh) * | 2020-12-25 | 2021-04-30 | 珠海格力电器股份有限公司 | 三电平永磁同步电机的控制方法及装置 |
-
2021
- 2021-05-24 CN CN202110565963.6A patent/CN113644852B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110299882A (zh) * | 2019-06-21 | 2019-10-01 | 浙江大学 | 基于混合供电型开绕组永磁同步电机的快速三矢量模型预测控制方法 |
CN111987943A (zh) * | 2020-07-23 | 2020-11-24 | 西安理工大学 | 一种永磁同步电机模型预测控制的鲁棒性能提升方法 |
US10819264B1 (en) * | 2020-08-04 | 2020-10-27 | Wolong Electric Group Co. Ltd. | Robust starting system and method for interior permanent magnet synchronous motor control |
CN112737452A (zh) * | 2020-12-25 | 2021-04-30 | 珠海格力电器股份有限公司 | 三电平永磁同步电机的控制方法及装置 |
Non-Patent Citations (2)
Title |
---|
EMRAH ZERDALI ET AL.: "Speed-Sensorless Finite Control Set Model Predictive Control of PMSM with Permanent Magnet Flux Linkage Estimation", 《2020 2ND GLOBAL POWER, ENERGY AND COMMUNICATION CONFERENCE (IEEE GPECOM2020)》, pages 1 - 6 * |
沈攀 等: "改进的永磁同步电机双矢量模型 预测转矩控制", 《新型工业化》, vol. 9, no. 1, pages 1 - 7 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114614712A (zh) * | 2022-03-27 | 2022-06-10 | 西北工业大学 | 基于模型预测-观测器的飞机电推进***扰动抑制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113644852B (zh) | 2024-02-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107154763B (zh) | 永磁同步电机无差拍直接转矩控制***及控制方法 | |
CN109768753B (zh) | 新型滑模观测器的无位置传感器永磁同步电机模型预测控制方法 | |
CN112701968B (zh) | 一种永磁同步电机模型预测控制鲁棒性能提升方法 | |
CN106788054B (zh) | 一种基于旋转高频注入法的无速度传感器控制方法 | |
CN103780187B (zh) | 永磁同步电机高动态响应电流控制方法及*** | |
CN112271970B (zh) | 永磁同步电机矢量控制方法、设备及存储介质 | |
CN104579083A (zh) | 永磁同步电机矢量控制方法及*** | |
CN110739893A (zh) | 一种改进自适应无轨迹卡尔曼滤波转动惯量辨识方法 | |
CN112217436B (zh) | 一种抑制电流测量误差引起的永磁同步电机转速脉动方法 | |
CN106385216B (zh) | 一种永磁同步电机电流预测控制稳态误差消除方法及*** | |
CN110912480A (zh) | 基于扩张状态观测器的永磁同步电机无模型预测控制方法 | |
CN111656674B (zh) | 用于电力转换装置的控制装置、控制方法、以及电动机驱动*** | |
CN104779873B (zh) | 一种用于pmsm伺服***的预测函数控制方法 | |
CN115864928A (zh) | 一种基于校正电流预测的pmsm模型参考自适应转速估算方法 | |
CN116526919A (zh) | 一种永磁同步电机伺服***及其电流预测控制方法、装置 | |
CN103986392B (zh) | 一种低速直驱式交流伺服***的控制方法 | |
CN113644852B (zh) | 用于永磁同步电动机的鲁棒三矢量模型预测磁链控制方法 | |
CN106712629B (zh) | 一种永磁同步电机的电流控制方法 | |
Horch et al. | Nonlinear integral backstepping control for induction motor drive with adaptive speed observer using super twisting strategy | |
CN110649851B (zh) | 异步电机多参数解耦在线辨识方法 | |
CN111835251A (zh) | 一种基于无速度传感的永磁同步电机高性能控制方法 | |
CN113965129B (zh) | 一种永磁同步电机控制***电流测量偏移误差的补偿方法 | |
CN112271972B (zh) | 带有电流误差校正的永磁同步电机直接转矩控制方法 | |
Bai et al. | Position and speed detection method based on adaptive extended moving-window linear regression for traction machine drives | |
CN113691187A (zh) | 一种无差拍预测转速伺服永磁电机控制方法及*** |
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 |