CN111564999B - 一种基于mras算法的电机低速测速方法 - Google Patents

一种基于mras算法的电机低速测速方法 Download PDF

Info

Publication number
CN111564999B
CN111564999B CN201911393819.8A CN201911393819A CN111564999B CN 111564999 B CN111564999 B CN 111564999B CN 201911393819 A CN201911393819 A CN 201911393819A CN 111564999 B CN111564999 B CN 111564999B
Authority
CN
China
Prior art keywords
speed
model
equation
time
permanent magnet
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
CN201911393819.8A
Other languages
English (en)
Other versions
CN111564999A (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.)
Shenzhen Graduate School Harbin Institute of Technology
Original Assignee
Shenzhen Graduate School Harbin Institute of 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 Shenzhen Graduate School Harbin Institute of Technology filed Critical Shenzhen Graduate School Harbin Institute of Technology
Priority to CN201911393819.8A priority Critical patent/CN111564999B/zh
Publication of CN111564999A publication Critical patent/CN111564999A/zh
Application granted granted Critical
Publication of CN111564999B publication Critical patent/CN111564999B/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/14Estimation or adaptation of machine parameters, e.g. flux, current or voltage
    • H02P21/18Estimation of position or speed
    • 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/13Observer control, e.g. using Luenberger observers or Kalman filters

Landscapes

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

Abstract

本发明涉及一种基于MRAS算法的电机低速测速方法,其采用永磁同步电机的运动学方程,结合M/T法以及模型参考自适应法进行低速转速测量,将M/T法求得的速度信息代入MRAS算法的参数输入侧,进而修正MRAS算法的速度观测数据,使得低速测量值比单独使用某一种测速方法更为精确。

Description

一种基于MRAS算法的电机低速测速方法
技术领域
本发明涉及一种测速方法,特别是指一种基于MRAS算法的电机低速测速方法。
背景技术
众所周知,永磁同步电机凭借其优良的力矩性能,具有更高的速度控制和位置控制性能,目前已经被广泛应用在了如机械臂关节、云台、相机镜头等需要更高精度和低速稳定性的场合。
永磁同步电机一般需要安装位置传感器,其用以采集电机转子位置信息用于电流坐标变换,转速控制或位置控制。常用的转子位置传感器为增量式光电编码器。在低速情况下,光电编码器所得到的角度存在一定的误差。同时,传统的测速方法测频(M)法低速精度不佳,测时(T)法及时频同测(M/T)法低速测量时间过长影响动态性能。传统的改变控制器结构或者前馈补偿的方法均引入了电流进行计算,经过派克变换后的电流存在由于转子位置不精确以及电流采样不精确而造成的误差。这在一定程度上影响了控制器在低转速条件下的控制性能,而此为传统技术的主要缺点。
发明内容
本发明所采用的技术方案为:一种基于MRAS算法的电机低速测速方法,采用永磁同步电机的运动学方程,结合M/T法以及模型参考自适应法进行低速转速测量,将M/T法求得的速度信息代入MRAS算法的参数输入侧,进而修正 MRAS算法的速度观测数据,使得低速测量值比单独使用某一种测速方法更为精确。
将斯特里贝克摩擦转矩模型代替运动学方程中的粘滞阻尼系数项,修正 MRAS算法的低速速度观测器模型,使得所提改进的低速测速方法更加精确。
一种基于MRAS算法的电机低速测速方法,M/T法中在一定时间间隔后,在测量时间内对编码器脉冲进行计数,同时时钟脉冲计数器进行计数,达到规定的测量时间后,下一个编码器脉冲上升沿到来时停止时钟脉冲计数器计数,并由此计数器值m2来确定检测时间。
M/T法测得的转速为:
Figure BDA0002345719550000021
M/T法测量时间由两部分组成,分别为M法测量所需时间Tg以及下一个编码器脉冲间隔时间ΔT,其检测时间为:
Figure BDA0002345719550000022
M/T法测速的最大误差εmax为:
Figure BDA0002345719550000023
计数器值m2与最大误差εmax的关系式:
Figure BDA0002345719550000031
模型参考自适应系统(MRAS),由参考模型,可调模型及自适应律组成,可调模型以及参考模型具有相同的输入,将参考模型的输出与可调模型的输出求差,将差值通过自适应律修改可调模型参数,使可调模型逼近参考模型。
永磁同步电机其电流方程为:
Figure BDA0002345719550000032
定义:
Figure BDA0002345719550000033
则将电流方程写成如下状态空间表达式:
Figure BDA0002345719550000041
状态矩阵A包含待辨识的转速信息,因此将此式作为可调模型,永磁同步电机本体作为参考模型,具有相同的输入与输出参数。
利用永磁同步电机的运动学方程,分别代入上述M/T法、模型参考自适应法测得的转速。
采用斯特里贝克摩擦转矩模型代替运动学方程中的粘滞阻尼系数项,转动惯量项用模型参考自适应法估算速度计算:
Figure BDA0002345719550000044
式中Teˊ为电磁转矩计算值(N·m);ωE为当前M/T法测速结果 (rad/s);ωe为模型参考自适应法当前测速结果(rad/s)。
电机的转动惯量,最大静摩擦力,库伦摩擦力为固定参数,不随转速变化而改变,负载转矩恒定。
利用此时的电磁转矩计算值根据永磁同步电机的运动学方程可得:
Figure BDA0002345719550000042
式中
Figure BDA0002345719550000043
为转速估计值(rad/s)。
将公式与运动学方程相减得电磁转矩误差表达式:
Figure BDA0002345719550000051
式中Te *为电磁转矩真实值(N·m);
Figure BDA0002345719550000053
为转速真实值(rad/s)。
设置观测器增益为K,则在上式中引入电磁转矩误差修正项,则观测器方程为:
Figure BDA0002345719550000054
式中K′T为系数行向量,K′T=[1-K-1-1];
T为力矩列向量,
Figure BDA0002345719550000055
采用永磁同步电机的运动学方程,结合M/T法以及模型参考自适应法进行低速转速测量,将M/T法求得的速度信息代入MRAS算法的参数输入侧,进而修正MRAS算法的速度观测数据。
本发明的有益效果为:本发明依据电机的运动学方程设计转速观测器,提出了一种融合M/T法与模型参考自适应法(MRAS)的永磁同步电动机低速测速方法。本发明的技术能够有效的改善低速时的测速精确度。
附图说明
图1为本发明M/T法测速示意图。
图2为本发明配合观测器结构框图。
具体实施方式
一种基于MRAS算法的电机低速测速方法,其依据电机的运动学方程设计转速观测器,提出了一种融合M/T法与模型参考自适应法(MRAS)的永磁同步电动机低速测速方法。本发明的技术能够有效的改善低速时的测速精确度。
M/T法基本原理与误差分析。
如图1所示,M/T法即规定一定时间间隔后,在测量时间内对编码器脉冲进行计数,同时时钟脉冲计数器进行计数。达到规定的测量时间后,下一个编码器脉冲上升沿到来时停止时钟脉冲计数器计数,并由此计数器值m2来确定检测时间。
如图1所示,M/T法测得的转速为:
Figure BDA0002345719550000061
M/T法测量时间由两部分组成,分别是M法测量所需时间Tg,下一个编码器脉冲间隔时间ΔT,所以其检测时间为:
Figure BDA0002345719550000062
M/T法测速的最大误差为:
Figure BDA0002345719550000071
综合前面公式,得m2和εmax的关系式:
Figure BDA0002345719550000072
从上式可以看出,M/T法的测量误差与计时脉冲数有关,所以其具有较大的测量范围。低速情况下m2较大,具有更高的精度。有利于拓宽伺服系统的调速范围。在低速情况下,ΔT较大从而使测量时间较大,降低转速的采样频率,影响控制器性能。在高精度应用场合,由此带来控制上的延迟以及更大的控制量波动,会使伺服系统的性能达不到要求。
模型参考自适应算法(MRAS)设计。
模型参考自适应系统(Model Reference Adaptive System,MRAS),由参考模型(实际系统),可调模型及自适应律组成。可调模型和参考模型具有相同的输入,将参考模型的输出与可调模型的输出求差,将差值通过自适应律修改可调模型参数,使可调模型逼近参考模型。
以表贴式永磁同步电机为例,其电流方程:
Figure BDA0002345719550000073
定义:
Figure BDA0002345719550000081
则可以将电流方程写成如下状态空间表达式:
Figure BDA0002345719550000082
状态矩阵A包含待辨识的转速信息,因此将此式作为可调模型。永磁同步电机本体作为参考模型,具有相同的输入与输出参数。
改进基于MRAS算法的电机低速测速方法:
M/T法在测量时刻时能得出较高精度的转速计算结果,误差来源于计时脉冲数、计数脉冲数及编码器机械制造误差,而模型参考自适应法测量结果误差来源于电流采样及电机模型本身误差。两种方法的速度计算相互独立,不存在误差传递。
在测速结果计算时间上,M/T法要短于模型参考自适应法,并且模型参考自适应法因为含有积分环节而具有一定的结果延迟。模型参考自适应法的采样频率要高于M/T法。模型参考自适应法的误差来源多样,直接计算其误差较为困难,从而无法直接比较当前两种方法的当前时刻测速结果的优劣。为了避免电流测量误差影响测速方法优劣的判断,采用永磁同步电机的运动学方程进行分析:根据电机运动学方程,分别代入之前M/T法、模型参考自适应法测得的转速。在低速下,斯特里贝克效应对摩擦转矩的影响较大,为了方便分析,采用斯特里贝克摩擦转矩模型代替运动学方程中的粘滞阻尼系数项,由于在M/T 法采样时间内M/T法的计算结果不可导,所以转动惯量项用模型参考自适应法估算速度计算:
Figure BDA0002345719550000091
式中Teˊ——电磁转矩计算值(N·m);
ωE——当前M/T法测速结果(rad/s);
ωe——模型参考自适应法当前测速结果(rad/s)。
电机的转动惯量,最大静摩擦力,库伦摩擦力为固定参数,不随转速变化而改变,负载转矩恒定。
利用此时的电磁转矩计算值根据永磁同步电机的运动学方程可得:
Figure BDA0002345719550000101
式中
Figure BDA0002345719550000102
——转速估计值(rad/s)。
将公式与运动学方程相减得电磁转矩误差表达式:
Figure BDA0002345719550000103
式中Te *——电磁转矩真实值(N·m);
Figure BDA0002345719550000105
——转速真实值(rad/s)。
设置观测器增益为K,则在上式中引入电磁转矩误差修正项,则观测器方程为:
Figure BDA0002345719550000106
式中K′T——系数行向量,K′T=[1-K-1-1];
T——力矩列向量,
Figure BDA0002345719550000107
转速观测器结构框图如图2所示,本专利提出采用永磁同步电机的运动学方程,融合M/T法与模型参考自适应法进行低速转速测量,并且考虑了低转速情况下电机运动学方程与传统电机运动学方程的差异,用M/T法求得的速度信息代入MRAS的参数输入侧,修正MRAS的速度观测数据,使得低速测量值比单独使用某一种测速方法更加精确。
本专利的技术内容归纳如下:为了避免电流测量误差影响测速方法优劣的判断,本专利提出采用永磁同步电机的运动学方程,融合M/T法与模型参考自适应法进行低速转速测量,用M/T法求得的速度信息代入MRAS的参数输入侧,修正MRAS的速度观测数据,使得低速测量值比单独使用某一种测速方法更加精确。
考虑了低转速情况下电机运动学方程与传统电机运动学方程的差异,提出采用斯特里贝克摩擦转矩模型代替运动学方程中的粘滞阻尼系数项,修正 MRAS的低速速度观测器模型,使得所提改进的低速测速方法更加精确。

Claims (1)

1.一种基于MRAS算法的电机低速测速方法,其特征在于:
M/T法中在一定时间间隔后,在测量时间内对编码器脉冲进行计数,同时时钟脉冲计数器进行计数,达到规定的测量时间后,下一个编码器脉冲上升沿到来时停止时钟脉冲计数器计数,并由此计数器值m2来确定检测时间,M/T法测得的转速为:
Figure FDA0003278735270000011
M/T法测量时间由两部分组成,分别为M法测量所需时间Tg以及下一个编码器脉冲间隔时间ΔT,其检测时间为:
Figure FDA0003278735270000012
M/T法测速的最大误差εmax为:
Figure FDA0003278735270000013
计数器值m2与最大误差εmax的关系式:
Figure FDA0003278735270000014
模型参考自适应系统,由参考模型,可调模型及自适应律组成,可调模型以及参考模型具有相同的输入,将参考模型的输出与可调模型的输出求差,将差值通过自适应律修改可调模型参数,使可调模型逼近参考模型,
永磁同步电机其电流方程为:
Figure FDA0003278735270000021
定义:
Figure FDA0003278735270000022
则将电流方程写成如下状态空间表达式:
Figure FDA0003278735270000023
状态矩阵A包含待辨识的转速信息,因此将此式作为可调模型,永磁同步电机本体作为参考模型,具有相同的输入与输出参数,
利用永磁同步电机的运动学方程,分别代入上述M/T法、模型参考自适应法测得的转速,
采用斯特里贝克摩擦转矩模型代替运动学方程中的粘滞阻尼系数项,转动惯量项用模型参考自适应法估算速度计算:
Figure FDA0003278735270000031
式中Teˊ为电磁转矩计算值;ωE为当前M/T法测速结果;ωe为模型参考自适应法当前测速结果,
电机的转动惯量,最大静摩擦力,库伦摩擦力为固定参数,不随转速变化而改变,负载转矩恒定,
利用此时的电磁转矩计算值根据永磁同步电机的运动学方程可得:
Figure FDA0003278735270000032
式中
Figure FDA0003278735270000033
为转速估计值,
将公式与运动学方程相减得电磁转矩误差表达式:
Figure FDA0003278735270000041
式中Te *为电磁转矩真实值;
Figure FDA0003278735270000042
为转速真实值,
设置观测器增益为K,则在上式中引入电磁转矩误差修正项,则观测器方程为:
Figure FDA0003278735270000043
式中K′T为系数行向量,K′T=[1 -K -1 -1];
T为力矩列向量,
Figure FDA0003278735270000044
采用永磁同步电机的运动学方程,结合M/T法以及模型参考自适应法进行低速转速测量,将M/T法求得的速度信息代入MRAS算法的参数输入侧,进而修正MRAS算法的速度观测数据。
CN201911393819.8A 2019-12-30 2019-12-30 一种基于mras算法的电机低速测速方法 Active CN111564999B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911393819.8A CN111564999B (zh) 2019-12-30 2019-12-30 一种基于mras算法的电机低速测速方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911393819.8A CN111564999B (zh) 2019-12-30 2019-12-30 一种基于mras算法的电机低速测速方法

Publications (2)

Publication Number Publication Date
CN111564999A CN111564999A (zh) 2020-08-21
CN111564999B true CN111564999B (zh) 2021-12-14

Family

ID=72072922

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911393819.8A Active CN111564999B (zh) 2019-12-30 2019-12-30 一种基于mras算法的电机低速测速方法

Country Status (1)

Country Link
CN (1) CN111564999B (zh)

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5663847A (en) * 1995-03-27 1997-09-02 Abramovitch; Daniel Y. Rejection of disturbances on a disk drive by use of an accelerometer
CN101005263A (zh) * 2007-01-25 2007-07-25 上海交通大学 交流电机伺服系统速度控制方法
CN101800505A (zh) * 2010-03-12 2010-08-11 北京航空航天大学 一种磁悬浮飞轮转速控制方法
CN103684182A (zh) * 2013-11-14 2014-03-26 西安理工大学 一种永磁同步电机参数辨识方法
CN104393815A (zh) * 2014-12-16 2015-03-04 电子科技大学 一种基于复合估计转速的永磁同步电机容错控制装置
CN105680746A (zh) * 2015-11-20 2016-06-15 宋正强 一种利用在线粒子群优化算法来设计永磁同步电动机电流和速度控制器pi参数的方法
CN106026822A (zh) * 2016-06-13 2016-10-12 上海电气集团股份有限公司 伺服电机驱动系统的惯量在线辨识方法及系统
CN106646225A (zh) * 2016-10-10 2017-05-10 杭州戈虎达科技有限公司 一种便携式电机效率现场检测装置及方法
CN106998164A (zh) * 2017-04-12 2017-08-01 武汉理工大学 一种永磁同步电机故障诊断与容错控制系统及方法
CN107181438A (zh) * 2017-06-06 2017-09-19 哈尔滨工业大学深圳研究生院 基于改进型q‑mras的异步电机无速度传感器控制方法
CN109412486A (zh) * 2018-10-23 2019-03-01 深圳市深信创联智能科技有限责任公司 直线感应电机的速度观测方法

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5663847A (en) * 1995-03-27 1997-09-02 Abramovitch; Daniel Y. Rejection of disturbances on a disk drive by use of an accelerometer
CN101005263A (zh) * 2007-01-25 2007-07-25 上海交通大学 交流电机伺服系统速度控制方法
CN101800505A (zh) * 2010-03-12 2010-08-11 北京航空航天大学 一种磁悬浮飞轮转速控制方法
CN103684182A (zh) * 2013-11-14 2014-03-26 西安理工大学 一种永磁同步电机参数辨识方法
CN104393815A (zh) * 2014-12-16 2015-03-04 电子科技大学 一种基于复合估计转速的永磁同步电机容错控制装置
CN105680746A (zh) * 2015-11-20 2016-06-15 宋正强 一种利用在线粒子群优化算法来设计永磁同步电动机电流和速度控制器pi参数的方法
CN106026822A (zh) * 2016-06-13 2016-10-12 上海电气集团股份有限公司 伺服电机驱动系统的惯量在线辨识方法及系统
CN106646225A (zh) * 2016-10-10 2017-05-10 杭州戈虎达科技有限公司 一种便携式电机效率现场检测装置及方法
CN106998164A (zh) * 2017-04-12 2017-08-01 武汉理工大学 一种永磁同步电机故障诊断与容错控制系统及方法
CN107181438A (zh) * 2017-06-06 2017-09-19 哈尔滨工业大学深圳研究生院 基于改进型q‑mras的异步电机无速度传感器控制方法
CN109412486A (zh) * 2018-10-23 2019-03-01 深圳市深信创联智能科技有限责任公司 直线感应电机的速度观测方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
《光电编码器在伺服控制系统中的应用及误差分析》;赵尚武;《第十三届全国敏感元件与传感器学术会议论文集 下 STC》;20141231;1068-1071 *

Also Published As

Publication number Publication date
CN111564999A (zh) 2020-08-21

Similar Documents

Publication Publication Date Title
CN109945819B (zh) 一种永磁同步电机转子位置测量方法
CN109323711A (zh) 一种陀螺仪模态反转零位自校正方法及系统
CN111238535B (zh) 一种基于因子图的imu误差在线标定方法
CN101902187A (zh) 一种用于永磁电机低速运转的控制系统
JP6091523B2 (ja) サーボ制御装置
CN109067274A (zh) 基于干扰补偿的永磁同步电机调速的积分滑模控制方法
CN108959183B (zh) 角度传感器和旋转编码器的前向插值方法
CN110426033A (zh) 基于松耦合imu阵列导航系统的时间同步算法
CN111564999B (zh) 一种基于mras算法的电机低速测速方法
CN113702705B (zh) 一种双轴精密离心机相位同步测量系统及测量方法
US6677722B2 (en) Servo control method
CN112923966A (zh) 一种基于卡尔曼滤波的双读数头光电编码器角度估计方法
CN116499696B (zh) 一种提高风洞模型试验模型姿态角动态精度的方法
CN102608912B (zh) 主动式半捷联惯性测量装置驱动系统的精确控制方法
EP0450084B1 (en) Closed loop servo motor control method
CN109959391B (zh) 一种用于光纤陀螺动态性能评测的信号同步方法
CN112504527A (zh) 一种基于力平衡的小型稳定平台摩擦力矩测试方法
JPS6130770A (ja) 車速検出装置
CN107508512B (zh) 闭环步进电机的超低速预估器控制算法
CN104678301B (zh) 一种基于脉冲反馈检测电机运动学参数的方法
CN110657800B (zh) 一种位置测量组合导航系统的时间同步方法
CN110243462B (zh) 运输测振仪振动检测方法
Boggarpu et al. New learning algorithm for high-quality velocity measurement from low-cost optical encoders
CN109048907A (zh) 消除惯量变化的机器人控制方法
CN107271162B (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