CN111049441B - 无刷直流电机的模型预测与反演控制方法 - Google Patents

无刷直流电机的模型预测与反演控制方法 Download PDF

Info

Publication number
CN111049441B
CN111049441B CN202010001008.5A CN202010001008A CN111049441B CN 111049441 B CN111049441 B CN 111049441B CN 202010001008 A CN202010001008 A CN 202010001008A CN 111049441 B CN111049441 B CN 111049441B
Authority
CN
China
Prior art keywords
motor
current
control
rotation
equation
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
CN202010001008.5A
Other languages
English (en)
Other versions
CN111049441A (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.)
Shanghai Xin Yue Lian Hui Electronic Technology Co ltd
Original Assignee
Shanghai Xin Yue Lian Hui Electronic Technology Co ltd
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 Shanghai Xin Yue Lian Hui Electronic Technology Co ltd filed Critical Shanghai Xin Yue Lian Hui Electronic Technology Co ltd
Priority to CN202010001008.5A priority Critical patent/CN111049441B/zh
Publication of CN111049441A publication Critical patent/CN111049441A/zh
Application granted granted Critical
Publication of CN111049441B publication Critical patent/CN111049441B/zh
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
    • H02P21/0003Control strategies in general, e.g. linear type, e.g. P, PI, PID, using robust control
    • H02P21/0017Model reference adaptation, e.g. MRAS or MRAC, useful for control or parameter estimation

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Control Of Motors That Do Not Use Commutators (AREA)
  • Control Of Electric Motors In General (AREA)

Abstract

本发明通过一种无刷直流电机的模型预测与反演控制方法,实现电机在力矩干扰及负载作用下转角、转速的高精度、高稳定度控制。在电机数字控制系统中,在电机参数已知情况下,根据电机机械特性及电气特性,求解得到电机受控转动下转角、转速的状态演变规律。在每个控制周期中,通过测量电机转角、转速及电流等状态变量,按照电机控制任务所规划出来的、电机后续期望达到的转角及转速,根据电机受控转动的规律,反演出当前控制电机转动所需的电压、电流。当电机因为外界力矩干扰或者负载作用在转角、转速上存在偏差时,通过补偿控制,实现对外界力矩干扰及负载作用的补偿。本发明的电机控制响应速度快、控制精度高、干扰及负载下控制鲁棒性好。

Description

无刷直流电机的模型预测与反演控制方法
技术领域
本发明涉及定子为三相对称星型连接绕组的无刷直流电机的一种数字控制方法,具体涉及一种无刷直流电机的模型预测与反演控制方法,可广泛应用于该类型无刷直流电机在各种力矩负载及干扰情况下的单纯转速控制以及转角、转速的同步控制。
背景技术
无刷直流电机因结构简单、运行高效可靠、便于控制而得到广泛应用。无刷直流电机的高精度、高稳定度、敏捷快速转动控制是各种电动伺服转台、数控机床、操控及服务机器人实现和达到良好性能指标的重要保障。在各种任务和负载工况下,要求无刷直流电机具有响应及时、控制超调小、控制稳定时间短、控制静差小、对力矩干扰及负载鲁棒性强的控制品质。
电压、电流、转角、转速的高精度、高速测量技术,高性能数字控制器技术以及快速、高精度、大功率的数字功放驱动电路技术使得计算机数字控制技术得以在无刷直流电机上广泛应用,先进、复杂的控制算法因而能够应用于无刷直流电机的驱动控制。
目前,无刷直流电机的PID控制、模糊控制、鲁棒控制、滑模变结构控制、自适应控制等方法各自不同程度的存在控制精度/稳定度差、对干扰和负载鲁棒性不强、控制器设计人为主观影响因素较大等问题,且通常从系统稳定性角度出发设计,在控制过程中常存在超调、震荡等问题。
本发明申请从动力学系统可控性角度出发设计,按照系统运动学、动力学规律进行状态转移的预测和控制作用的反演,类似的无刷直流电机控制技术尚未见报道。
发明内容
本发明提供一种无刷直流电机的模型预测与反演控制方法,实现无刷直流电机的转角、转速同步高精度高动态控制和误差补偿,能够实现力矩扰动下高精度的位置与速度伺服控制。
为了达到上述目的,本发明提供一种无刷直流电机的模型预测与反演控制方法,包含:
步骤A,在当前采样控制周期k,预测未来采样控制周期k+n1的电机转动状态偏差,n1≥2;
步骤B,反演k+1~k+n1-1各采样控制周期的控制电流指令序列,并对干扰及负载力矩造成的电机转动状态偏差进行补偿,利用转动状态偏差反演求解补偿电流,与计算的控制电流之和作为电流输出的指令值;
步骤C,反演k时刻的控制电压指令。根据当前k周期测量电流值和步骤B中的所得电流指令,预测k+n2周期的电流差值,n2∈[1,n1-1],反向推演k周期的指令电压,作为控制输出。在下一采样控制节拍,返回步骤A,重新进行控制电压指令解算。
可选地,所述的模型预测与反演控制方法,进一步包含以下步骤:
步骤S1,在连续时间系统中,建立电机电气特性方程、机械特性方程、转动运动学方程;
步骤S2,分别建立电机电气特性、机械特性、转动特性的离散状态转移方程;该步骤S2之后,执行所述的步骤A~步骤C;
所述的步骤S1中,在连续时间系统中,根据已知的电机反电势系数、转子转动惯量、阻尼系数、每相绕组的电阻、自感、两绕组间互感,实时测量采样得到的电机三相电流、电机转速、电机转角,以及在采样控制时刻通过电角度查表计算得出的反电势函数,以电机转角作为因变量,电机转速作为变量,得到状态空间形式的电机转动运动学方程;以电机转速作为状态变量,以三相电流作为控制输入,得到状态空间形式的电机机械特性方程;以三相电流作为状态变量,三相电压作为控制输入,得到状态空间形式的电机电气特性方程;
所述的步骤S2中,分别对电机转动运动学方程、机械特性方程、电气特性方程进行求解,得到电机转动特性、机械特性、电气特性的状态转移方程;根据状态转移方程,在离散时间系统中,建立电机转动特性、机械特性、电气特性的离散状态转移方程。
可选地,所述的步骤A中,以当前采样控制时刻作为k时刻,根据当前测量采样到的电机转角、转速、三相电流,利用电机转动特性、机械特性、电气特性的离散状态转移方程,预测未来采样控制时刻k+n1的电机转角、转速,根据k+n1时刻的转角、转速的规划值,计算得到k+n1时刻的转角、转速的偏差预测值。
可选地,所述的步骤B中,根据电机转动特性、机械特性离散状态转移方程,采用广义逆方法,计算k+1~k+n1-1各采样控制时刻的控制电流序列,得到所有的k+1~k+n1-1控制电流指令序列的集合。根据之前的电机实际转动状态的偏差信息,反向推演补偿负载电流,其与之前的控制电流的和作为指令电流值输出。
可选地,所述的步骤C中,比较k+n2周期电机电流的预测值以及指令值得到电流偏差预测值,根据电机电气特性离散状态转移方程,采用广义逆方法,计算得到k周期的控制电压指令输出给电源。
可选地,对应所述的步骤A~步骤C中,在单个控制周期中,
电流的预测,对应于电气特性离散方程
U2,k+1=G1,kU2,k+Q1,kU1,k+HkEk 式1
转速的预测,对应于机械特性离散方程
ωk+1=G2,kωk+Q2,kU2,k+PkWk 式2
转角的预测,对应于转动特性离散方程
θk+1=θk+Tωk 式3
其中,ix|x=a,b,c为三相绕组中电流
U1=[ua ub uc]T,ux|x=a,b,c为三相端电压值;U1,k是k周期时U1所示电压的实测值;
ω为电机转速;U2=[ia ib ic]T,ix|x=a,b,c为三相绕组中电流,W=[Tl],T1为负载转矩,T为一个控制周期的时长;
U2,k是k周期时U2所示电流的实测值、ωk是k周期转速的实测值;
下标k、(k+1)分别表示当前控制周期和下一控制周期;
Figure GDA0003919883650000041
Figure GDA0003919883650000042
Wk为此周期内负载;
I为单位阵,Ek为当前周期的反电势,Hk为反电势作用常值矩阵;
Figure GDA0003919883650000051
Figure GDA0003919883650000052
E=[ea eb ec]T,ex|x=a,b,c为三相反电势,
R、L、M分别为每相绕组的电阻、电感与绕组间互感;
Figure GDA0003919883650000053
J为转子转动惯量,Bf为阻尼系数,K为反电势系数,
fxe)|x=a,b,c为反电势波形函数,θe为转子电角度。
可选地,对应步骤C,在单个控制周期中,计算电压的通式为:
Figure GDA0003919883650000054
根据上式的最小二乘最小范数特解,得到控制电压并驱动电源输出:
Figure GDA0003919883650000055
其中,
Figure GDA0003919883650000056
为电流参考值,以计算得到的电流U2,k+1代入;
Figure GDA0003919883650000057
为Q1,k的广义逆,V为任意3维列向量;
Figure GDA0003919883650000058
Figure GDA0003919883650000059
为Q2,k的广义逆;
Figure GDA00039198836500000510
为根据之前转速偏差求得的补偿电流量;
kc为正值补偿系数;Δωj是j周期转速的状态偏差;
Figure GDA0003919883650000061
为给定的转速参考值,由此计算得到电流参考值
Figure GDA0003919883650000062
进而计算得到控制电压并驱动电源输出,实现电机的转速模型预测与反演控制。
可选地,对应步骤A~步骤C,通过两个控制周期,实现对电机转角转速的同步控制;其中,以电机转角、转速X=[θ ω]T作为状态变量,建立预测方程
Figure GDA0003919883650000063
式中,
Figure GDA0003919883650000064
控制电流为
Figure GDA0003919883650000065
式中,Kc=dia([kc1 kc2]),为正定对角矩阵,kc1、kc2分别为调节角度、角速度误差补偿速度而引入的正值补偿系数;
Figure GDA0003919883650000066
为k+2时刻给定的转角转速参考值,Xk为k时刻的转角转速实测值,ΔXj是j周期转角转速的状态偏差;
取控制电流中的U2,k+1输出;下一控制周期中,更新转速转角的参考值和实测值,重新计算并取相应的U2,k+1输出,代入到控制电压的算式5中作为电流参考值
Figure GDA0003919883650000067
并更新电流实测值U2,k和反电势值Ek,计算三相电压值
Figure GDA0003919883650000068
实现对电机转角转速的同步控制。
可选地,所述电机是定子为三相对称星型连接绕组的无刷直流电机。
本发明所述无刷直流电机的模型预测与反演控制方法,与现有技术相比,其优点和有益效果是:
(1)跟踪控制精度高,响应速度快,对电机参数及扰动具有强鲁棒性。
本方法基于模型及实时状态和参考状态,定量反演控制量。参数偏差及扰动最终反映于状态偏差,并依据状态偏差及扰动影响机理直接求解控制补偿量,具有较强的鲁棒性。此方法为根据状态偏差实时调整参数的精确定量控制方法,具有自适应和精确定量控制的特点。
(2)适用范围广。
根据不同的任务规划电机的转角、转速参考轨迹,可以实现高精度的位置伺服控制,如定点控制和轨迹跟踪控制。能应用于导弹舵机的高动态伺服舵机,数控机床,机器人关节,转台等高精度伺服系统中。
(3)方法简单有效,便于调节。
控制器参数由控制周期及电机本体参数确定,且随电机状态时变自适应,便于应用在数字控制系统中,适用于同种类不同参数的电机,有利于批量设计控制器。
附图说明
图1是无刷直流电机模型预测与反演控制流程图;
图2是无刷直流电机模型预测与反演控制框图。
具体实施方式
本发明公开了一种定子为三相对称星型连接绕组的无刷直流电机(文中后述电机均指该型电机)的模型预测与反演控制方法,可以实现电机在力矩干扰及负载作用下转角、转速的高精度、高稳定度控制。
本发明的构思是,根据无刷直流电机的电气特性、机械特性方程和运动学方程,得到电机转动在控制作用下的变化规律,并根据当前采样时刻的状态,预测未来一个或多个采样时刻的电机状态,并根据未来对应时刻期望达到的状态,反向推演当前时刻应定量施加的控制量,使系统在有限控制周期内快速达到期望状态(高动态)。针对电机参数偏差及扰动等因素造成的状态偏差,根据偏差量反演累加得到控制补偿输入,使伺服系统对参数和扰动具有较强的鲁棒性。控制器参数由电机本体参数和电机可测运动状态共同确定,因而所设计的控制器具有时变自适应和较强鲁棒性的特点。
以下将结合附图和实施例对本发明作进一步地描述。
本发明基于系统模型,通过控制作用对系统状态的影响机理,根据当前状态和期望状态预测未来状态偏差,反向解算当前控制周期应施加的控制作用,实现无刷直流电机的高精度控制,算法流程如图1所示。
本发明提供一种无刷直流电机的模型预测与反演控制方法中,包含:
步骤S1,在连续时间系统中,建立电机电气特性方程、机械特性方程、转动运动学方程;
步骤S2,分别建立电机电气特性、机械特性、转动特性的离散状态转移方程;该步骤S2之后,执行下述的步骤A~步骤C;
步骤A,在当前采样控制周期k,预测未来采样控制周期k+n1的电机状态偏差,n1≥2;
步骤B,反演k+1~k+n1-1各采样控制周期的控制电流指令序列,并对干扰及负载力矩造成的电机转动状态偏差进行补偿,利用转动状态偏差反演求解补偿电流,与计算的控制电流之和作为电流输出的指令值;
步骤C,反演k时刻的控制电压指令。根据当前k周期测量电流值和步骤B中的所得电流指令,预测k+n2周期的电流差值,n2∈[1,n1-1],反向推演k周期的指令电压,作为控制输出。在下一采样控制节拍,返回步骤A,重新进行控制电压指令解算。
所述的步骤S1中,在连续时间系统中,根据已知的电机反电势系数、转子转动惯量、阻尼系数、每相绕组的电阻、自感、两绕组间互感,实时测量采样得到的电机三相电流、电机转速、电机转角,以及在采样控制时刻通过电角度查表计算得出的反电势函数,以电机转角作为因变量,电机转速作为变量,得到状态空间形式的电机转动运动学方程;以电机转速作为状态变量,以三相电流作为控制输入,得到状态空间形式的电机机械特性方程;以三相电流作为状态变量,三相电压作为控制输入,得到状态空间形式的电机电气特性方程;
所述的步骤S2中,分别对电机转动运动学方程、机械特性方程、电气特性方程进行求解,得到电机转动特性、机械特性、电气特性的状态转移方程;根据状态转移方程,在离散时间系统中,建立电机转动特性、机械特性、电气特性的离散状态转移方程。
可选地,所述的步骤A中,以当前采样控制时刻作为k时刻,根据当前测量采样到的电机转角、转速、三相电流,利用电机转动特性、机械特性、电气特性的离散状态转移方程,预测未来采样控制时刻k+n1的电机转角、转速,根据k+n1时刻的转角、转速的规划值,计算得到k+n1时刻的转角、转速的偏差预测值。
可选地,所述的步骤B中,根据电机转动特性、机械特性离散状态转移方程,采用广义逆方法,计算k+1~k+n1-1各采样控制时刻的控制电流序列,得到所有的k+1~k+n1-1控制电流指令序列的集合。根据之前的电机实际转动状态的偏差信息,反向推演补偿负载电流,其与之前的控制电流的和作为指令电流值输出。
可选地,所述的步骤C中,比较k+n2周期电机电流的预测值以及指令值得到电流偏差预测值,根据电机电气特性离散状态转移方程,采用广义逆方法,计算得到k周期的控制电压指令输出给电源。
下面对预测单控制周期实施方式的特例进行介绍:
步骤S1,建立无刷直流电机的电气特性方程、机械特性方程、运动学方程。
无刷直流电机的状态空间形式的电气特性方程为
Figure GDA0003919883650000101
式中,U2=[ia ib ic]T为状态变量,ix|x=a,b,c为三相绕组中电流,U1=[ua ub uc]T为控制输入,ux|x=a,b,c为三相端电压值,
Figure GDA0003919883650000102
Figure GDA0003919883650000103
E=[ea eb ec]T,ex|x=a,b,c为三相反电势,R、L、M分别为每相绕组的电阻、电感与绕组间互感。
状态空间形式的机械特性方程为
Figure GDA0003919883650000104
式中,ω为电机转速;U2=[ia ib ic]T为控制输入,W=[Tl]为负载转矩,相关矩阵为A2=[-Bf/J],
Figure GDA0003919883650000111
J为转子转动惯量,Bf为阻尼系数,K为反电势系数,fxe)|x=a,b,c为反电势波形函数,θe为转子电角度。
电机参数K、J、Bf、R、L、M通常均为已知,三相电流、转速ω可以测量得到,反电势函数在每个控制周期中可通过电角度查表计算得出。
转动运动学方程为
Figure GDA0003919883650000113
θ为转子角位置。
步骤S2,基于步骤S1,建立数字控制系统中电机电气特性、机械特性、转动特性的离散状态转移方程。
当采样周期T足够小时,近似认为三相端电压和反电势在一个测控周期内不变,有电气特性离散方程
U2,k+1=G1,kU2,k+Q1,kU1,k+HkEk (3)
式中,
Figure GDA0003919883650000114
下标k、(k+1)分别表示当前控制周期和下一控制周期,I为单位阵,Ek为当前周期的反电势,Hk为反电势作用常值矩阵。
机械特性离散方程为
ωk+1=G2,kωk+Q2,kU2,k+PkWk (4)
式中,
Figure GDA0003919883650000112
Wk为此周期内负载。ωk、U2,k分别为测量的周期的转速值、电流值。因为A2、P为常值,因此G2,k、Pk为常值。
转动特性离散方程为
θk+1=θk+Tωk (5)
基于步骤S1、S2,执行步骤A~步骤C:
步骤A,在当前采样控制周期k,预测未来某采样控制周期k+2的电机状态偏差。在k+1周期内不加控制情况下,k+2周期的电机的转速偏差为
Figure GDA0003919883650000121
其中
Figure GDA0003919883650000122
为k+2周期的期望值,
步骤B,反演k+1周期的控制电流指令,并对干扰及负载力矩造成的电机转动状态偏差进行补偿,
Figure GDA0003919883650000123
式(6)中,
Figure GDA00039198836500001213
为Q2,k的广义逆。
Figure GDA0003919883650000124
为补偿控制电流,Δωj为之前的j周期转速偏差,kc为正值补偿系数,用以调节补偿程度。
Figure GDA0003919883650000125
为计算所得的k+1周期的电流值,作为参考指令值输入到电流控制环。
步骤C,根据步骤B中参考电流值
Figure GDA0003919883650000126
在无控下预测k+1周期电流偏差,反演得到当前k应施加的电压值
Figure GDA0003919883650000127
其通解为
Figure GDA0003919883650000128
根据上式的最小二乘最小范数特解,得到控制电压并驱动电源输出:
Figure GDA0003919883650000129
其中,U2,k为测量的k周期的电流值。
特殊地,对转角、转速进行同步控制的情况下,在步骤B中,以电机转角、转速X=[θω]T作为状态变量,建立预测方程,
Figure GDA00039198836500001210
式中,
Figure GDA00039198836500001211
T为控制周期。
控制电流为
Figure GDA00039198836500001212
式中,Kc=dia([kc1 kc2]),为正定对角矩阵,kc1、kc2分别为调节角度、角速度误差补偿速度而引入的正值补偿系数;
Figure GDA0003919883650000131
为k+2时刻给定的转角转速参考值,Xk为k时刻的转角转速实测值,ΔXj是j周期转角转速的状态偏差;
取控制电流中的U2,k+1输出;下一控制周期中,更新转速转角的参考值和实测值,重新计算并取相应的U2,k+1输出,代入到式(8)中作为电流参考值
Figure GDA0003919883650000132
并更新电流实测值U2,k和反电势值Ek,计算三相电压值
Figure GDA0003919883650000133
从而实现对电机转角转速的同步控制。
综上所述,本发明提供的无刷直流电机的模型预测与反演控制方法,在电机数字控制系统中,在电机参数已知情况下,根据电机机械特性及电气特性,求解得到电机受控转动下转角、转速的状态演变规律。在每个控制周期中,通过测量电机转角、转速及电流等状态变量,按照电机控制任务所规划出来的、电机后续期望达到的转角及转速,根据电机受控转动的规律,反演出当前控制电机转动所需的电压、电流。当电机因为外界力矩干扰或者负载作用在转角、转速上存在偏差时,通过补偿控制,实现对外界力矩干扰及负载作用的补偿。与现有技术相比,本发明具有电机控制响应速度快、控制精度高、干扰及负载下控制鲁棒性好的优点。
上述为预测单周期或两周期的特例说明,对于转角、转速和电流的多周期预测和反演控制量的方法也在本发明要求的权利范围内。尽管本发明的内容已经通过上述优选实施例作了详细介绍,但应当认识到上述的描述不应被认为是对本发明的限制。在本领域技术人员阅读了上述内容后,对于本发明的多种修改和替代都将是显而易见的。因此,本发明的保护范围应由所附的权利要求来限定。

Claims (9)

1.一种无刷直流电机的模型预测与反演控制方法,其特征在于,包含:
步骤A,在当前采样控制周期k,预测未来采样控制周期k+n1的电机转动状态偏差,n1≥2;
步骤B,反演k+1~k+n1-1各采样控制周期的控制电流指令序列,并对干扰及负载力矩造成的电机转动状态偏差进行补偿,利用转动状态偏差反演求解补偿电流,与计算的控制电流之和作为电流输出的指令值;
步骤C,反演k时刻的控制电压指令;根据当前k周期测量电流值和步骤B中的所得电流指令,预测k+n2周期的电流差值,n2∈[1,n1-1],反向推演k周期的指令电压,作为控制输出;在下一采样控制节拍,返回步骤A,重新进行控制电压指令解算。
2.如权利要求1所述的模型预测与反演控制方法,其特征在于,
进一步包含以下步骤:
步骤S1,在连续时间系统中,建立电机电气特性方程、机械特性方程、转动运动学方程;
步骤S2,分别建立电机电气特性、机械特性、转动特性的离散状态转移方程;该步骤S2之后,执行所述的步骤A~步骤C;
所述的步骤S1中,在连续时间系统中,根据已知的电机反电势系数、转子转动惯量、阻尼系数、每相绕组的电阻、自感、两绕组间互感,实时测量采样得到的电机三相电流、电机转速、电机转角,以及在采样控制时刻通过电角度查表计算得出的反电势函数,以电机转角作为因变量,电机转速作为变量,得到状态空间形式的电机转动运动学方程;以电机转速作为状态变量,以三相电流作为控制输入,得到状态空间形式的电机机械特性方程;以三相电流作为状态变量,三相电压作为控制输入,得到状态空间形式的电机电气特性方程;
所述的步骤S2中,分别对电机转动运动学方程、机械特性方程、电气特性方程进行求解,得到电机转动特性、机械特性、电气特性的状态转移方程;根据状态转移方程,在离散时间系统中,建立电机转动特性、机械特性、电气特性的离散状态转移方程。
3.如权利要求2所述的模型预测与反演控制方法,其特征在于,
所述的步骤A中,以当前采样控制时刻作为k时刻,根据当前测量采样到的电机转角、转速、三相电流,利用电机转动特性、机械特性、电气特性的离散状态转移方程,预测未来采样控制时刻k+n1的电机转角、转速,根据k+n1时刻的转角、转速的规划值,计算得到k+n1时刻的转角、转速的偏差预测值。
4.如权利要求3所述的模型预测与反演控制方法,其特征在于,
所述的步骤B中,根据电机转动特性、机械特性离散状态转移方程,采用广义逆方法,计算k+1~k+n1-1各采样控制时刻的控制电流序列,得到所有的k+1~k+n1-1控制电流指令序列的集合;根据之前的电机实际转动状态的偏差信息,反向推演补偿负载电流,其与之前的控制电流的和作为指令电流值输出。
5.如权利要求4所述的模型预测与反演控制方法,其特征在于,
所述的步骤C中,比较k+n2周期电机电流的预测值以及指令值得到电流偏差预测值,根据电机电气特性离散状态转移方程,采用广义逆方法,计算得到k周期的控制电压指令输出给电源。
6.如权利要求5所述的模型预测与反演控制方法,其特征在于,
对应所述的步骤A~步骤C中,在单个控制周期中,
电流的预测,对应于电气特性离散方程
U2,k+1=G1,kU2,k+Q1,kU1,k+HkEk 式1
转速的预测,对应于机械特性离散方程
ωk+1=G2,kωk+Q2,kU2,k+PkWk 式2
转角的预测,对应于转动特性离散方程
θk+1=θk+Tωk 式3
其中,U1=[ua ub uc]T,ux|x=a,b,c为三相端电压值;U1,k是k周期时U1所示电压的实测值;
ω为电机转速;U2=[ia ib ic]T,ix|x=a,b,c为三相绕组中电流;W=[Tl],T1为负载转矩,T为一个控制周期时长;
U2,k是k周期时U2所示电流的实测值、ωk是k周期转速的实测值;
下标k、(k+1)分别表示当前和下一控制周期;
Figure FDA0003919883640000031
Figure FDA0003919883640000041
Wk为此周期内负载;
I为单位阵,Ek为当前周期的反电势,Hk为反电势作用常值矩阵;
Figure FDA0003919883640000042
Figure FDA0003919883640000043
E=[ea eb ec]T,ex|x=a,b,c为三相反电势,
R、L、M分别为每相绕组的电阻、电感与绕组间互感;
Figure FDA0003919883640000044
J为转子转动惯量,Bf为阻尼系数,K为反电势系数,
fxe)|x=a,b,c为反电势波形函数,θe为转子电角度。
7.如权利要求6所述的模型预测与反演控制方法,其特征在于,
对应步骤C,在单个控制周期中,计算电压的通式为:
Figure FDA0003919883640000045
根据上式的最小二乘最小范数特解,得到控制电压并驱动电源输出:
Figure FDA0003919883640000046
其中,
Figure FDA0003919883640000047
为电流参考值,以计算得到的电流U2,k+1代入;
Figure FDA0003919883640000048
为Q1,k的广义逆,V为任意3维列向量;
Figure FDA0003919883640000051
Figure FDA0003919883640000052
为Q2,k的广义逆;
Figure FDA0003919883640000053
为根据之前转速偏差求得的补偿电流量;
kc为正值补偿系数;Δωj是j周期转速的状态偏差;
Figure FDA0003919883640000054
为给定的转速参考值,由此计算得到电流参考值
Figure FDA0003919883640000055
进而计算得到控制电压并驱动电源输出,实现电机的转速模型预测与反演控制。
8.如权利要求7所述的模型预测与反演控制方法,其特征在于,
对应步骤A~步骤C,通过两个控制周期,实现对电机转角转速的同步控制;其中,以电机转角、转速X=[θ ω]T作为状态变量,建立预测方程
Figure FDA0003919883640000056
式中,
Figure FDA0003919883640000057
控制电流为
Figure FDA0003919883640000058
式中,Kc=dia([kc1 kc2]),为正定对角矩阵,kc1、kc2分别为调节角度、角速度误差补偿速度而引入的正值补偿系数;
Figure FDA0003919883640000059
为k+2时刻给定的转角转速参考值,Xk为k时刻的转角转速实测值,ΔXj是j周期转角转速的状态偏差;
取控制电流中的U2,k+1输出;下一控制周期中,更新转速转角的参考值和实测值,重新计算并取相应的U2,k+1输出,代入到控制电压的算式5中作为电流参考值
Figure FDA00039198836400000510
并更新电流实测值U2,k和反电势值Ek,计算三相电压值
Figure FDA00039198836400000511
实现对电机转角转速的同步控制。
9.如权利要求1~8中任意一项所述的模型预测与反演控制方法,其特征在于,所述电机是定子为三相对称星型连接绕组的无刷直流电机。
CN202010001008.5A 2020-01-02 2020-01-02 无刷直流电机的模型预测与反演控制方法 Active CN111049441B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010001008.5A CN111049441B (zh) 2020-01-02 2020-01-02 无刷直流电机的模型预测与反演控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010001008.5A CN111049441B (zh) 2020-01-02 2020-01-02 无刷直流电机的模型预测与反演控制方法

Publications (2)

Publication Number Publication Date
CN111049441A CN111049441A (zh) 2020-04-21
CN111049441B true CN111049441B (zh) 2023-03-10

Family

ID=70243548

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010001008.5A Active CN111049441B (zh) 2020-01-02 2020-01-02 无刷直流电机的模型预测与反演控制方法

Country Status (1)

Country Link
CN (1) CN111049441B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112000014B (zh) * 2020-09-03 2022-07-15 上海航天控制技术研究所 一种操控机械臂的模型预测与反演的控制方法
CN112984701B (zh) * 2021-03-15 2022-09-20 广东积微科技有限公司 空调电机转速的调节控制方法、装置、单片机设备及介质
CN113668178B (zh) * 2021-09-08 2024-04-19 海信冰箱有限公司 洗衣机的电机控制方法、装置、计算机可读介质及洗衣机
CN113848499B (zh) * 2021-09-28 2022-11-11 珠海格力电器股份有限公司 驱动电源的参数监控方法、装置及驱动电源系统
CN115645792A (zh) * 2022-10-31 2023-01-31 徐工消防安全装备有限公司 一种泵控工作平台调平系统及云梯消防车

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20000003132A (ko) * 1998-06-26 2000-01-15 윤종용 인버터형 공기조화기의 이상상태 검출방법
CN103326780A (zh) * 2013-06-21 2013-09-25 中国科学院空间科学与应用研究中心 基于压缩感知接收机的自由空间光通信apt系统及方法
CN108233781A (zh) * 2018-01-18 2018-06-29 南京理工大学 一种基于非线性干扰观测器的直流无刷电机自适应反演滑模控制方法
CN108248605A (zh) * 2018-01-23 2018-07-06 重庆邮电大学 一种智能车辆轨迹跟随的横纵向协调控制方法
CN108437991A (zh) * 2018-04-11 2018-08-24 厦门大学 一种智能电动汽车自适应巡航控制系统及其方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20000003132A (ko) * 1998-06-26 2000-01-15 윤종용 인버터형 공기조화기의 이상상태 검출방법
CN103326780A (zh) * 2013-06-21 2013-09-25 中国科学院空间科学与应用研究中心 基于压缩感知接收机的自由空间光通信apt系统及方法
CN108233781A (zh) * 2018-01-18 2018-06-29 南京理工大学 一种基于非线性干扰观测器的直流无刷电机自适应反演滑模控制方法
CN108248605A (zh) * 2018-01-23 2018-07-06 重庆邮电大学 一种智能车辆轨迹跟随的横纵向协调控制方法
CN108437991A (zh) * 2018-04-11 2018-08-24 厦门大学 一种智能电动汽车自适应巡航控制系统及其方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Improved Dynamic Behavior in BLDC Drives Using Model Predictive Speed and Current Control;Araz Darba et al.;《IEEE Transactions on Industrial Electronics》;20150908;正文第728-740页 *
卫星姿态大角度机动的轨迹规划和模型预测与反演控制;谭天乐 等;《飞控与探测》;20190125;第2卷(第1期);正文第17-23页 *

Also Published As

Publication number Publication date
CN111049441A (zh) 2020-04-21

Similar Documents

Publication Publication Date Title
CN111049441B (zh) 无刷直流电机的模型预测与反演控制方法
CN107425769B (zh) 一种永磁同步电机调速系统的自抗扰控制方法及系统
CN111555680B (zh) 一种永磁同步电机无差拍预测电流控制方法
CN110190795B (zh) 一种永磁同步电机级联式鲁棒预测电流控制方法
CN110649845B (zh) 基于鲁棒广义预测控制的光电转台位置跟踪控制方法
CN110798116B (zh) 一种基于电枢模型独立前馈补偿的电机矢量复合控制器
CN106788059B (zh) 高动态电机伺服控制系统的延时补偿方法
CN112422004A (zh) 一种永磁同步电机弱磁控制模式下的扰动抑制方法
CN111211717A (zh) 非奇异滑模结构的ipmsm无位置传感器电机闭环结构控制方法
CN112953335A (zh) 一种永磁同步电机有限时间自适应复合控制方法和系统
CN110557069B (zh) 转子运行参数估计方法、电机控制系统和自抗扰控制器
CN116526919A (zh) 一种永磁同步电机伺服系统及其电流预测控制方法、装置
CN114567226A (zh) 考虑铁损耗的pmsm调速系统变趋近律滑模控制方法
JP2830274B2 (ja) 可変速駆動装置の電流制御方式
CN113992098A (zh) 一种永磁同步电机稳态性能的提升方法、系统
CN113794419A (zh) 一种基于增量式电机模型的无差拍电流预测鲁棒控制方法
CN110995098B (zh) 一种永磁同步电机的反演控制方法
CN109889113B (zh) 一种基于自抗扰控制的永磁电机变速扫描控制系统
CN116805849A (zh) 一种永磁同步电机的连续集模型预测控制方法
CN115051601A (zh) 变速旋转直流电机伺服系统的扰动补偿和跟踪控制方法
CN113078865B (zh) 内置式永磁同步电机无传感器控制方法
CN111181462B (zh) 基于变步长神经网络的表贴式永磁同步电机参数辨识方法
CN114826080A (zh) 一种基于积分补偿的无差拍电流预测控制方法和系统
CN114285342A (zh) 一种永磁同步电机模型预测直接速度同步控制方法
Hu et al. Research on injection molding machine drive system based on model predictive control

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