CN108075704B - 电压模型定子磁链估计方法 - Google Patents
电压模型定子磁链估计方法 Download PDFInfo
- Publication number
- CN108075704B CN108075704B CN201711362389.4A CN201711362389A CN108075704B CN 108075704 B CN108075704 B CN 108075704B CN 201711362389 A CN201711362389 A CN 201711362389A CN 108075704 B CN108075704 B CN 108075704B
- Authority
- CN
- China
- Prior art keywords
- stator
- flux linkage
- coordinate system
- direct current
- 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
Links
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/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
- 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
- H02P21/141—Flux estimation
-
- 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/01—Asynchronous machines
Abstract
本发明公开了一种新颖的改进电压模型定子磁链估计方法。该估计方法根据采集并运算得到的静止坐标系下的异步电机定子反感应电动势和定子角频率,然后经过直流偏置消除模块得到不含直流量的定子反感应电动势ysα、ysβ;最后通过磁链计算模块得到估计的α轴定子磁链估计的β轴定子磁链相对于传统的技术方案,本发明方案具有如下优点:1、消除电压模型磁链观测器的直流偏置和稳态误差问题;2、磁链估计的动态响应速度较快。
Description
技术领域
本发明涉及一种磁链观测器,具体说是一种新颖的改进电压模型定子磁链估计方法。
背景技术
异步电机磁链的准确观测一直是电机高性能控制的关键课题。利用观测磁链可以实现电机速度、转矩等参数的估计,进而实现电机无速度、转矩闭环控制。基于电压纯积分的磁链观测方案,因其结构简单、易于实现等优点而受到人们青睐。然而在实际应用中,电压模型固有的纯积分问题限制了磁链观测的精度。针对这一问题,国内外学者提出了一系列改进方案。
题为“定子磁链观测器低通滤波器的改进”(何志明,廖勇,向大为,中国电机工程学报,2008年第28卷第18期61-65页)的文章提出一种基于补偿加低通滤波器的方案。该方法虽然可以有效提高磁链观测速度,但是无法完全消除估计磁链中的直流偏置。2015年IEEE文献“Improved stator flux estimator for speed sensorless induction motordrives”(“基于改进磁链观测器的异步电机无速度控制”——2015年IEEE电力电子期刊),提出一种变低通滤波器截止频率的改进方案。该方法虽然可以进一步提高磁链观测速度、降低估计磁链中的高频噪声,但是仍然无法完全消除估计磁链中的直流偏置。题为“一种能消除直流偏置和稳态误差的电压型磁链观测器”(张星,瞿文龙,陆海峰,电工电能新技术,2006年第25卷第1期39-43页)的文章提出一种带通滤波器串联补偿网络的方案。虽然该方案可以消除积分初值与直流漂移问题,但是该方案磁链观测的动态响应速度仍然有待提升。
发明内容
本发明的目的是为了克服现有改进电压模型磁链观测方法无法兼顾磁链观测速度和精度的问题,提出一种新颖的改进方案。本发明首先采集异步电机的定子频率ωe,静止αβ坐标系下的定子电压usα、usβ和定子电流isα、isβ,在计算得到定子反电动势后,依次经过直流偏置消除模块、磁链计算模块,最终实现定子磁链的准确观测。
本发明的技术方案按照以下步骤实现。
步骤1:采集异步电机的定子频率ωe,静止αβ坐标系下的定子电压分量usα、usβ和静止αβ坐标系下的定子电流分量isα、isβ,根据公式(1)计算得到定子反感应电动势esa、esβ;
其中Rs为定子电阻;
步骤2:根据步骤1得到的定子反感应电动势esa、esβ,经过直流偏置消除模块得到不含直流量的定子反感应电动势ysα、ysβ;
所述直流偏置消除模块,包括一个补偿网络计算方程和一个高通滤波器计算方程,其表达式分别为:
补偿网络计算方程:
高通滤波器计算方程:
其中cα、cβ为中间变量,kH为高通滤波器补偿系数,ysα、ysβ为不含直流量的定子反感应电动势,ysα、ysβ即直流偏置消除模块的输出,||代表绝对值运算,s为拉普拉斯算子,sgn()为符号函数;
步骤3:根据步骤2得到的不含直流量的定子反感应电动势ysα、ysβ,通过磁链计算模块得到估计的α轴定子磁链估计的β轴定子磁链
所述磁链计算模块,包括一个补偿网络计算方程和一个低通滤波器计算方程,其表达式分别为:
补偿网络计算方程:
低通滤波器计算方程:
其中dα、dβ为中间变量,kL为低通滤波器补偿系数,||代表绝对值运算,s为拉普拉斯算子,sgn()为符号函数。
优选地,步骤1中所述静止αβ坐标系下的定子电压分量usα、usβ的采集方式包括以下两种:
第一种,采样得到实时异步电机线电压Uab、Ubc,经过公式(2)计算得到静止αβ坐标系下的定子电压分量usα、usβ:
第二种,直接采用电机控制器运算单元计算出逆变器调制电压来获得静止αβ坐标系下的定子电压分量usα、usβ。
优选地,步骤1中所述静止αβ坐标系下的定子电流分量isα、isβ的采集步骤如下:
1)采样得到异步电机三相定子电流isa、isb、isc;
2)根据公式(3)计算得到异步电机静止αβ坐标系下的定子电流isα、isβ:
本发明相对于现有技术的有益效果是:
1、与基于低通滤波器的改进电压模型相比,本发明可以完全消除估计磁链中的直流偏置。
2、与基于带通滤波器的改进电压模型相比,本发明可以提高磁链观测动态响应速度。
附图说明
图1为本发明执行的流程图;
图2为本发明提出算法的实现框图;
图3为电机定子频率变化图;
图4为电机定子频率按图3变化时,本发明方案和传统方案A的估计定子磁链的αβ轴误差对比波形。
图5为电机定子频率按图3变化时,本发明方案和传统方案B的估计定子磁链的αβ轴误差对比波形。
具体实施方式
下面结合附图对本发明的技术方案作进一步详细的描述。
图1为本发明执行的流程图,图2为本发明提出算法的实现框图,由图1、图2所示,本发明实施例包括以下步骤:
步骤1:采集异步电机的定子频率ωe,静止αβ坐标系下的定子电压分量usα、usβ和静止αβ坐标系下的定子电流分量isα、isβ,根据公式(1)计算得到定子反感应电动势esa、esβ;
其中Rs为定子电阻。
所述静止αβ坐标系下的定子电压分量usα、usβ的采集方式包括以下两种:
第一种,采样得到实时异步电机线电压Uab,Ubc,经过公式(2)计算得到静止αβ坐标系下的定子电压分量usα、usβ:
具体的,首先利用霍尔电压传感器采样定子A、B两相之间的线电压和定子B、C两相之间的线电压,然后输入带有滤波器功能的采样调理电路,实现电压信号的采样,获得实时异步电机线电压Uab、Ubc,并在数字信号处理芯片中通过编程对获得的定子电压信号Uab和Ubc进行如公式(2)所示的坐标变换,计算得到异步电机静止αβ坐标系下的定子电流usα、usβ。
第二种,直接采用电机控制器运算单元计算出逆变器调制电压来获得静止αβ坐标系下的定子电压分量usα、usβ。
所述静止αβ坐标系下的定子电流分量isα、isβ的采集步骤如下:
1)利用霍尔电流传感器采样得到异步电机三相定子电流isa、isb、isc;
2)将异步电机三相定子电流isa、isb、isc输入电流采样通道;在数字信号处理芯片中通过编程对获得的电机三相电流isa、isb、isc,根据公式(3)计算得到异步电机静止αβ坐标系下的定子电流isα、isβ:
步骤2:根据步骤1得到的定子反感应电动势esa、esβ,经过直流偏置消除模块得到不含直流量的定子反感应电动势ysα、ysβ。
所述直流偏置消除模块,包括一个补偿网络计算方程和一个高通滤波器计算方程,其表达式分别为:
补偿网络计算方程:
高通滤波器计算方程:
其中cα、cβ为中间变量,kH为高通滤波器补偿系数,ysα、ysβ为不含直流量的定子反感应电动势,ysα、ysβ即直流偏置消除模块的输出,||代表绝对值运算,s为拉普拉斯算子,sgn()为符号函数。
步骤3:根据步骤2得到的不含直流量的定子反感应电动势ysα、ysβ,通过磁链计算模块得到估计的α轴定子磁链估计的β轴定子磁链
所述磁链计算模块,包括一个补偿网络计算方程和一个低通滤波器计算方程,其表达式分别为:
补偿网络计算方程:
低通滤波器计算方程:
其中dα、dβ为中间变量,kL为低通滤波器补偿系数,||代表绝对值运算,s为拉普拉斯算子,sgn()为符号函数。
本实施例中的仿真电机参数:额定功率10kW,额定电压220V,额定电流32A,定子电阻0.076Ω,极对数2,额定频率100Hz,高通滤波器补偿系数kH和低通滤波器补偿系数kL分别为:
kL=10kH (9)
当电机定子频率按照图3变化时,图4给出了本发明方案和传统方案A估计定子磁链的αβ轴误差对比波形,图5给出了本发明方案和传统方案B估计定子磁链的αβ轴误差对比波形,其中,磁链误差值由估计定子磁链减去仿真中电机模型输出的参考磁链获取。
从图4和5可以看出:本发明提出的方案在电机全频率区间内都可以实现定子磁链的准确估计,同时在定子频率变化时,本发明方案具有更快的磁链观测速度。传统方案A为2015年IEEE文献“Improved stator flux estimator for speed sensorless inductionmotor drives”(“基于改进磁链观测器的异步电机无速度控制”——2015年IEEE电力电子期刊)提出的一种可变截止频率的低通滤波器方案。传统方案B为题为“一种能消除直流偏置和稳态误差的电压型磁链观测器”(张星,瞿文龙,陆海峰,电工电能新技术,2006年第25卷第1期39-43页)的文章提出一种带通滤波器串联补偿网络的方案。
Claims (3)
1.一种电压模型定子磁链估计方法,其特征在于,包括如下步骤:
步骤1:采集异步电机的定子频率ωe,静止αβ坐标系下的定子电压分量usα、usβ和静止αβ坐标系下的定子电流分量isα、isβ,根据公式(1)计算得到定子反感应电动势esa、esβ;
其中Rs为定子电阻;
步骤2:根据步骤1得到的定子反感应电动势esa、esβ,经过直流偏置消除模块得到不含直流量的定子反感应电动势ysα、ysβ;
所述直流偏置消除模块,包括一个补偿网络计算方程和一个高通滤波器计算方程,其表达式分别为:
补偿网络计算方程:
高通滤波器计算方程:
其中cα、cβ为中间变量,kH为高通滤波器补偿系数,ysα、ysβ为不含直流量的定子反感应电动势,ysα、ysβ即直流偏置消除模块的输出,||代表绝对值运算,s为拉普拉斯算子,sgn()为符号函数;
步骤3:根据步骤2得到的不含直流量的定子反感应电动势ysα、ysβ,通过磁链计算模块得到估计的α轴定子磁链估计的β轴定子磁链
所述磁链计算模块,包括一个补偿网络计算方程和一个低通滤波器计算方程,其表达式分别为:
补偿网络计算方程:
低通滤波器计算方程:
其中dα、dβ为中间变量,kL为低通滤波器补偿系数,||代表绝对值运算,s为拉普拉斯算子,sgn()为符号函数。
2.根据权利要求1所述的一种电压模型定子磁链估计方法,其特征在于,步骤1中所述静止αβ坐标系下的定子电压分量usα、usβ的采集方式包括以下两种:
第一种,采样得到实时异步电机线电压Uab、Ubc,经过公式(2)计算得到静止αβ坐标系下的定子电压分量usα、usβ:
第二种,直接采用电机控制器运算单元计算出逆变器调制电压来获得静止αβ坐标系下的定子电压分量usα、usβ。
3.根据权利要求1所述的一种电压模型定子磁链估计方法,其特征在于,步骤1中所述静止αβ坐标系下的定子电流分量isα、isβ的采集步骤如下:
1)采样得到异步电机三相定子电流isa、isb、isc;
2)根据公式(3)计算得到异步电机静止αβ坐标系下的定子电流isα、isβ:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711362389.4A CN108075704B (zh) | 2017-12-18 | 2017-12-18 | 电压模型定子磁链估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711362389.4A CN108075704B (zh) | 2017-12-18 | 2017-12-18 | 电压模型定子磁链估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108075704A CN108075704A (zh) | 2018-05-25 |
CN108075704B true CN108075704B (zh) | 2019-12-27 |
Family
ID=62158500
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711362389.4A Active CN108075704B (zh) | 2017-12-18 | 2017-12-18 | 电压模型定子磁链估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108075704B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111953242A (zh) * | 2019-04-29 | 2020-11-17 | 博世力士乐(西安)电子传动与控制有限公司 | 电动机工作状态的估计方法与装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101795109A (zh) * | 2010-02-11 | 2010-08-04 | 中国船舶重工集团公司第七一二研究所 | 一种异步电机矢量控制器 |
CN102780444A (zh) * | 2012-08-06 | 2012-11-14 | 无锡艾柯威科技有限公司 | 一种电机定子磁链估算方法 |
CN103414424A (zh) * | 2013-06-19 | 2013-11-27 | 山西潞安环保能源开发股份有限公司 | 交流电机定子磁链估计方法 |
CN104362930A (zh) * | 2014-12-09 | 2015-02-18 | 南京国电南自新能源科技有限公司 | 基于定子反电动势的同步电机转速的快速算法 |
CN105846748A (zh) * | 2016-05-17 | 2016-08-10 | 华中科技大学 | 一种基于矢量变换与信号滤波的定子磁链计算方法 |
-
2017
- 2017-12-18 CN CN201711362389.4A patent/CN108075704B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101795109A (zh) * | 2010-02-11 | 2010-08-04 | 中国船舶重工集团公司第七一二研究所 | 一种异步电机矢量控制器 |
CN102780444A (zh) * | 2012-08-06 | 2012-11-14 | 无锡艾柯威科技有限公司 | 一种电机定子磁链估算方法 |
CN103414424A (zh) * | 2013-06-19 | 2013-11-27 | 山西潞安环保能源开发股份有限公司 | 交流电机定子磁链估计方法 |
CN104362930A (zh) * | 2014-12-09 | 2015-02-18 | 南京国电南自新能源科技有限公司 | 基于定子反电动势的同步电机转速的快速算法 |
CN105846748A (zh) * | 2016-05-17 | 2016-08-10 | 华中科技大学 | 一种基于矢量变换与信号滤波的定子磁链计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108075704A (zh) | 2018-05-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107659237B (zh) | 一种永磁同步电机的无模型无差拍电流预测控制装置及其方法 | |
CN105610369B (zh) | 一种基于滑模观测器的异步电机磁链观测方法 | |
CN110022106B (zh) | 一种基于高频信号注入的永磁同步电机无位置传感器控制方法 | |
CN107482982B (zh) | 一种基于铁损模型的异步电机矢量控制方法 | |
CN106788049B (zh) | 基于级联滑模观测器的无速度传感器转矩控制系统及方法 | |
CN112910329A (zh) | 一种基于移动均值滤波器的pmsm无位置传感器控制方法 | |
CN110808703B (zh) | 考虑铁损电阻的永磁同步电机转速及转子位置估计方法 | |
CN107222135A (zh) | 一种直流无刷电机无位置传感器控制系统换相控制方法 | |
CN111211717B (zh) | 非奇异滑模结构的ipmsm无位置传感器电机闭环控制方法 | |
CN111239661B (zh) | 基于固定点采样的三相电流传感器误差校正系统及方法 | |
CN103281030A (zh) | 一种混合励磁电机无位置传感器矢量控制方法 | |
CN112117943B (zh) | 一种新型ipmsm高频方波注入无位置传感器控制 | |
CN114665772B (zh) | 永磁同步电机的控制方法 | |
CN108075704B (zh) | 电压模型定子磁链估计方法 | |
CN114944801A (zh) | 一种基于新息自适应扩展卡尔曼的pmsm无位置传感器控制方法 | |
CN113364375B (zh) | 变结构电流调节器的pmsm驱动系统无传感器控制方法 | |
CN108649849B (zh) | 一种简单的无传感器永磁同步电机速度估测方法 | |
CN116191964B (zh) | 基于矢量谐振控制器的正弦型电励磁双凸极电机控制方法 | |
CN108063573B (zh) | 基于直流偏置双闭环补偿的改进电压模型磁链观测方法 | |
CN108683376A (zh) | 一种基于滤波原理的交流电机电流测量偏移误差补偿方法 | |
CN116232154A (zh) | 基于复矢量离散滤波器的电机参数估计及控制方法和系统 | |
CN107994829A (zh) | 一种基于带通滤波器的改进电压模型定子磁链估计系统及方法 | |
CN114696695A (zh) | 永磁同步直线电机无传感器控制性能的提升方法 | |
CN110112977B (zh) | 一种基于永磁同步电机相电压测量的无位置传感器控制方法 | |
CN113258841A (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 |