CN101769998A - 开关磁阻电机磁链特性测量方法 - Google Patents

开关磁阻电机磁链特性测量方法 Download PDF

Info

Publication number
CN101769998A
CN101769998A CN201010018146A CN201010018146A CN101769998A CN 101769998 A CN101769998 A CN 101769998A CN 201010018146 A CN201010018146 A CN 201010018146A CN 201010018146 A CN201010018146 A CN 201010018146A CN 101769998 A CN101769998 A CN 101769998A
Authority
CN
China
Prior art keywords
current
phase current
phase
angle
voltage
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
Application number
CN201010018146A
Other languages
English (en)
Other versions
CN101769998B (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.)
Nantong University
Original Assignee
Nantong University
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 Nantong University filed Critical Nantong University
Priority to CN2010100181460A priority Critical patent/CN101769998B/zh
Publication of CN101769998A publication Critical patent/CN101769998A/zh
Application granted granted Critical
Publication of CN101769998B publication Critical patent/CN101769998B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Control Of Electric Motors In General (AREA)

Abstract

本发明公开了一种开关磁阻电机磁链特性测量方法,包括:校准位置传感器,确定测量电路的系数;移除负载,运行开关磁阻电机,调节开关角度,使电机的相电流近似为平顶波形状;减小直流电压使得相电流的峰值尽量的小,采样和存储数据矩阵;增大直流电压,使相电流的峰值缓慢增大;加上负载,改变开通角或者关断角,检测需要的数据组,采样和存储这些数据;使电流峰值逐渐增大,直到电机容许的电流上限值,记录下不同电流情况下的数据向量,直到最大的相电流达到需要电流峰值的最大电流imax。本发明整个过程只用到常规的开关磁阻电机调速系统和存储式数字示波器,不需要专门的转子机械定位和夹紧装置,测量设备很容易得到。

Description

开关磁阻电机磁链特性测量方法
技术领域:
本发明涉及一种开关磁阻电机磁链特性测量方法。
背景技术:
磁链特性是开关磁阻电机建模、仿真以及基于模型的控制等研究的基础。每相磁链关于相电流和位置角度的函数可表示为ψ(i,θ),由于开关磁阻电机的双凸极结构和磁路饱和性质,导致其磁链呈现高度的非线性特性,准确描述ψ(i,θ)函数也比较困难。
文献提出了多种获取开关磁阻电机磁链特性的方法。这些方法可分为有限元法(FEA)、直接测量方法和间接测量方法。
(1)有限元法:磁链可以通过专业的软件和详细的开关磁阻电机设计数据计算得到。然而专业的软件很昂贵,而且开关磁阻电机设计数据用户通常无法得到。另一方面,有限元的计算过程很耗时,并且计算结果精度在很大程度上依赖于对有限元网格的具体划分。
(2)直接测量方法:直接测量方法是借助安装在电机内部的磁传感器测量磁通密度来获得磁链,由于此种方法不方便而且价格昂贵,在实际中使用较少。
(3)间接测量方法:间接测量方法的思想就是通过其它容易测量的变量来计算磁链。这些变量可以是电压、电流、电阻或者转矩。由于这些用于计算磁链的变量相对容易测量得到,所以间接测量方法被人们普遍使用。间接测量方法可以被归纳为测试线圈法、电压脉冲注入法和转矩法。大部分的间接检测方法都需要专门的转子机械定位和夹紧装置,并且在测量过程中当电磁转矩很大时,转子会发生一定程度的转动,影响了测量结果。
发明内容:
本发明的目的在于提供一种不需要专门的转子机械定位和夹紧装置,测量设备很容易得到,并且方法比较简单的开关磁阻电机磁链特性测量方法。
本发明的技术解决方案是:
一种开关磁阻电机磁链特性测量方法,包括测量用于计算磁链的变量,其特征是:测量用于计算磁链的变量依次包括下列步骤:
(1)校准位置传感器,确定测量电路的系数;
(2)移除负载,运行开关磁阻电机,调节开关角度,使电机稳定运行于中速区,并且相电流近似为平顶波形状;减小直流电压使得相电流的峰值尽量的小,测量相电压u、相电流i、转子位置角度θ,并将测量得到的相电压u、相电流i、转子位置角度θ与采样时间t以
Figure G2010100181460D00021
的矩阵形式储存在存储式数字示波器中,记为S1,其中
Figure G2010100181460D00022
Figure G2010100181460D00023
都是列向量;
(3)增大直流电压,加负载,调节负载大小和开关角度,使相电流的峰值增加一个电流步长(电流步长影响测量结果精度,可选择为1A或更小),使电机稳定运行于中速区,且相电流近似为平顶波形状,检测相电压u、相电流i、转子位置角度θ,并将测量得到的相电压u、相电流i、转子位置角度θ与采样时间t以
Figure G2010100181460D00031
的矩阵形式储存在存储式数字示波器中,采样和存储这些数据记为S2
(4)重复步骤(3),使电流峰值逐渐增大,直到电机容许的电流上限值imax,记录下不同电流情况下的数据向量,标记为S3、S4…Sn
本发明的测量磁链的方法,整个过程只用到常规的开关磁阻电机调速系统(SRD)和数字存储示波器(DSO),不需要专门的转子机械定位和夹紧装置,测量设备很容易得到,并且测量和计算过程比较简单。
附图说明:
下面结合附图和实施例对本发明作进一步说明。
图1是本发明方法框架图。
图2是相电感和检测信号(u-相电压,i-相电流,θ-转子位置)关系示意图。
图3是ψ-i平面上的磁链轨迹和imax *图。
图4是磁链计算的流程图。
图5是用于磁链计算的典型测量波形图(5V/div,32V/div,0.8A/div,1ms/div),图中:1-转子位置,2-相电压,3-相电流。
图6是实验与FEA方法所得磁链的比较图。
图7是APC模式下试验与仿真结果对比图,其中θon=-27.5°,θoff=-14°,TL=2.7N-m。
图8是CCC模式下试验与仿真结果对比图,其中:θon=-27°,θoff=-6.5°,TL=1.8N-m
具体实施方式:
本方法的框架图如图1所示,由常用的开关磁阻电机调速系统和数字存储示波器组成。为了得到磁链值,需使SRD运行在稳态,然后测量其中某相的电压u、电流i、位置角θ。由于电机每相的参数是相同的,所以任意相都可作为测量相。测量得到的u、i、θ和采样时间t以
Figure G2010100181460D00041
的矩阵形式储存在DSO中,其中
Figure G2010100181460D00042
Figure G2010100181460D00043
都是是列向量。通过调节SRD的运行参数可以测得和存储不同的实验需求属性的样本数据组。在获得这些数据以后,磁链就可以在这些数据基础上根据下面提到的算法在PC机上计算得到。
电感L和测得的信号(θ,u,i)之间的关系如图2所示,其中τr是转子周期,θu是转子非对齐位置,θa是转子对齐位置,θon是主开关器件的开通角,θ*是转子介于θu和θa之间的任意位置角度。对应于某个具体的电流i*的位置角度θ*处的磁链可以由公式(1)计算得到。
ψ ( i * , θ * ) = ∫ 0 t * [ u ( t ) - ri ( t ) ] dt + ψ ( 0 , θ on ) - - - ( 1 )
其中t*是对应转子位置θ*的瞬时时刻,r是相电阻,ψ(0,θon)是当电流i=0也就是在开通角θon处的磁链值。特别注意,对于任何θon由于对应的电流都是0,使得ψ(0,θon)=0。
由于采样存储的数据是离散的,要计算离散的磁链ψ(i*,θ*),公式(1)的连续积分形式就需转化成如公式(2)所示的离散积分形式。
ψ ( p ) = ψ ( p - 1 ) + Δt 2 { u ( p ) + u ( p - 1 ) - R [ i ( p ) - i ( p - 1 ) ] } ψ ( 1 ) = 0 , p = 2,3,4 . . . . . . - - - ( 2 )
其中Δt是DSO的采样周期,采样点数1对应的时间是t=0时刻,p与第t*时刻相对应。在实际测量中,变量u(p)和i(p)如果不正好是样本矩阵S中的数据,则可以通过下面的方法计算得到:
由于数据矩阵S是在SRD运行在稳态时得到的,所以可以假设在[tu,ta]区间内,转子速度是常数,tu和ta是分别对应转子位置θu和θa的瞬时时刻,这些时刻可以通过在S矩阵中寻找信号θ的跳变边沿得到。所以t*可以由公式(3)计算得到。
t * = t u + θ * - θ u t r / 2 ( t a - t u ) - - - ( 3 )
那么u(t*)和i(t*)相当于公式(2)中的u(p)和i(p),可以根据t*由数值插值计算得到。
通过调节SRD的运行参数,可以得到位置θ*处与不同电流ij *相对应的一系列样本数据Sj,其中1≤j≤n,n是Sj总组数,那么对应于第j组离散电流ij *在位置角θ*处的磁链ψ(ij *,θ*)可以由公式(2)求得。
如果将在位置角θ*出的不同电流和与其对应的磁链组合成矩阵
Figure G2010100181460D00052
然后根据ij *大小对Z*进行排序,那么在Z*中对应于位置角θ*的最小电流值imin *和最大电流值imax *就可以得到。在位置角θ*处的对应电流范围为
Figure G2010100181460D00053
内的磁链值ψ(i,θ*)就可以由数值插值方法计算得到。由于i=0时,ψ(i,θ*)=0。所以在电流范围为内的磁链值ψ(i,θ*)也可以由数值插值方法计算得到。从而位置角度θ*处对应的整个电流范围为内的ψ(i,θ*)都可以求得。电流步长Δi*越小,插值精度越高。由于θ*是介于θu和θa之间的任意位置角度,那么在角度范围为θu≤θ≤θa,电流范围为
Figure G2010100181460D00061
的ψ(i,θ)都可以根据上述的方法求得。
实验得到的角度位置θ*处的最大电流imax *总是在小于最大电流imax的范围内变动。不同角度θ*对应的最大电流imax *在ψ-i平面图上的典型形状如图3所示,图3中的阴影部分表示电流值大于imax *处的磁链值没有得到。但是,在实际应用中所需的磁链范围只需要如图3所示的磁链轨迹范围内的磁链值,磁链轨迹内电流的范围都没有超过imax *。即使需要的磁链范围对应的电流最大值大于imax *,也可以通过在测量过程中统一增大最大电流imax *的峰值得到图3阴影部分的磁链值。
为了有效地实施本文所述方法,在测试过程中可遵循如下所列原则:
(1)使开关磁阻电机运行在稳态,转速不能太低,例如:大于1000rev/min(1000转/分)时有利于提高磁链计算精度。
(2)使得相电流形状接近于平顶电流,使得电流的平顶部分最大范围地介于θu和θa区域。这样的好处是使对于某个位置角每组样本的相电流不一样,以至于样本数据能够最佳地计算磁链值。
测量方法可以归纳为以下步骤:
以三相12/8极开关磁阻电机为例,其转子周期为45°,由于磁链特性的对称性,只需测量0°~22.5°范围内各离散位置的磁链特性。设所需建立的磁链特性数据表格形式为ψ(i,θ),其中i为电流,范围从0A到8A,步长为1A,θ为转子角度,范围从0°到22.5°,步长1.5°,测量步骤如下:
(1)校准位置传感器,确定测量电路的系数;
(2)移除负载,运行开关磁阻电机,调节开关角度,使电机稳定运行于中速区,并且相电流近似为平顶波形状;减小直流电压使得相电流的峰值尽量的小,测量相电压u、相电流i、转子位置角度θ,并将测量得到的相电压u、相电流i、转子位置角度θ与采样时间t以
Figure G2010100181460D00071
的矩阵形式储存在存储式数字示波器中,记为S1,其中
Figure G2010100181460D00072
Figure G2010100181460D00073
都是列向量;
(3)增大直流电压,加负载,调节负载大小和开关角度,使相电流的峰值增加一个电流步长(电流步长影响测量结果精度,可选择为1A或更小),并使电机稳定运行于中速区,且相电流近似为平顶波形状,检测相电压u、相电流i、转子位置角度θ,并将测量得到的相电压u、相电流i、转子位置角度θ与采样时间t以
Figure G2010100181460D00074
的矩阵形式储存在存储式数字示波器中,采样和存储这些数据记为S2
(4)重复步骤(3),使电流峰值逐渐增大,直到电机容许的电流上限值imax,记录下不同电流情况下的数据向量,标记为S3、S4...Sn
根据存储的n组数据矩阵Sj,其中1≤j≤n,磁链可以由前文所描述的方法和图4所示的计算过程得到。
实验测量与验证:
实验环境由一台1kW三相12/8极开关磁阻电机调速系统和数字存储示波器构成,如图5所示。电机转子位置检测使用传统的光电位置传感器,相电流、电压检测使用霍尔传感器,涡流测功机作为负载,安捷伦6014A型DSO用于检测和存储数据。
A实验测量与有限元方法结果的比较
实验测量步骤按照前文提到的方法进行。典型的用于磁链计算的转子位置,相电压,相电流波形如图5所示。由图5可以看出相电流接近于平顶电流,如前文所述。
本次实验的最大电流设为imax=8.5A,采样的样本组数为25组。磁链根据采样的数据由前文所述的方法计算得到,为了验证试验得到的磁链的精确度,将实验结果与FEA方法得到的磁链值做比较,两种方法得到的磁链比较结果如图6所示。从图6可以看出本文的方法计算的磁链值与FEA方法得到的磁链很接近。并且每个离散位置对应的最大电流imax *也在图6上,他们形成的基本形状跟图3预测的基本相似。
B仿真与实验结果的对比
为了进一步验证本方法的可行性,将实验与仿真的结果作对比。使用本方法计算所得的磁链进行动态仿真,得到一组电压和电流波形;设置同样的控制参数进行实验,得到一组电压和电流波形。如果仿真波形和实验结果吻合,则表示所测量的磁链特性是可靠的。
仿真使用Simulink软件,仿真系统的建立进行了下列假设:
1)直流电压是常数;
2)功率变换器是理想的;
3)忽略开关磁阻电机的互感,涡流和磁滞效应。
开关磁阻电机的控制方式有角度位置控制(APC)和电流斩波控制(CCC)两种,在APC和CCC模式下的仿真结果与实验结果分别如图7和图8所示。从图7和图8可以看出,实验结果与仿真结果有少许的差别,这个现象是由于下列因素引起的:
1)实验所得磁链值的误差;
2)建立仿真系统进行了一定的理想化假设;
3)实验中开通角与关断角的控制不可能与仿真值完全一致;
4)电压和电流的测量误差。
在工程实际中这些因素都是不可避免的。由图7和图8所示结果可见,这些因素的影响已被控制在较小范围之内。

Claims (1)

1.一种开关磁阻电机磁链特性测量方法,包括测量用于计算磁链的变量,其特征是:测量用于计算磁链的变量依次包括下列步骤:
(1)校准位置传感器,确定测量电路的系数;
(2)移除负载,运行开关磁阻电机,调节开关角度,使电机稳定运行于中速区,并且相电流近似为平顶波形状;减小直流电压使得相电流的峰值尽量的小,测量相电压u、相电流i、转子位置角度θ,并将测量得到的相电压u、相电流i、转子位置角度θ与采样时间t以的矩阵形式储存在存储式数字示波器中,记为S1,其中
Figure F2010100181460C00012
Figure F2010100181460C00013
都是列向量;
(3)增大直流电压,加负载,调节负载大小和开关角度,使相电流的峰值增加一个电流步长,并使电机稳定运行于中速区,且相电流近似为平顶波形状,检测相电压u、相电流i、转子位置角度θ,并将测量得到的相电压u、相电流i、转子位置角度θ与采样时间t以的矩阵形式储存在存储式数字示波器中,采样和存储这些数据记为S2
(4)重复步骤(3),使电流峰值逐渐增大,直到电机容许的电流上限值imax,记录下不同电流情况下的数据向量,标记为S3、S4...Sn
CN2010100181460A 2010-01-18 2010-01-18 开关磁阻电机磁链特性测量方法 Expired - Fee Related CN101769998B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2010100181460A CN101769998B (zh) 2010-01-18 2010-01-18 开关磁阻电机磁链特性测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2010100181460A CN101769998B (zh) 2010-01-18 2010-01-18 开关磁阻电机磁链特性测量方法

Publications (2)

Publication Number Publication Date
CN101769998A true CN101769998A (zh) 2010-07-07
CN101769998B CN101769998B (zh) 2012-07-18

Family

ID=42502979

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2010100181460A Expired - Fee Related CN101769998B (zh) 2010-01-18 2010-01-18 开关磁阻电机磁链特性测量方法

Country Status (1)

Country Link
CN (1) CN101769998B (zh)

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102368674A (zh) * 2011-06-22 2012-03-07 南通大学 一种开关磁阻电机转子位置解算方法及系统
CN102621492A (zh) * 2012-04-02 2012-08-01 湘潭大学 开关磁阻发电机检测装置
CN103869245A (zh) * 2014-01-20 2014-06-18 哈尔滨理工大学 基于磁共能修正的开关磁阻电机磁链曲线测试方法及装置
CN104014076A (zh) * 2014-04-21 2014-09-03 广州三瑞医疗器械有限公司 一种镇痛波形检测装置及其方法
CN104330755A (zh) * 2014-08-26 2015-02-04 南通大学 一种开关磁阻电机磁链特性在线测量方法
CN104539212A (zh) * 2015-01-20 2015-04-22 中国矿业大学 一种开关磁阻电机磁链特性在线测量方法
CN106441664A (zh) * 2016-11-21 2017-02-22 南通大学 一种开关磁阻电机转矩测量仪
US10079566B1 (en) 2017-07-14 2018-09-18 Caterpillar Inc. Pulse injection current pattern for rotor position estimation
CN109309462A (zh) * 2018-05-10 2019-02-05 烟台仙崴机电有限公司 一种能源互联网新型发电机系统效率优化控制方法
CN111190128A (zh) * 2018-11-15 2020-05-22 北京自动化控制设备研究所 一种磁阻电机铁磁材料bh特性检测算法
CN113676107A (zh) * 2021-08-24 2021-11-19 山东交通学院 一种开关磁阻电机转子位置检测装置原点校准方法及系统

Cited By (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102368674B (zh) * 2011-06-22 2013-11-06 南通大学 一种开关磁阻电机转子位置解算方法及系统
CN102368674A (zh) * 2011-06-22 2012-03-07 南通大学 一种开关磁阻电机转子位置解算方法及系统
CN102621492A (zh) * 2012-04-02 2012-08-01 湘潭大学 开关磁阻发电机检测装置
CN103869245B (zh) * 2014-01-20 2017-02-01 哈尔滨理工大学 基于磁共能修正的开关磁阻电机磁链曲线测试方法
CN103869245A (zh) * 2014-01-20 2014-06-18 哈尔滨理工大学 基于磁共能修正的开关磁阻电机磁链曲线测试方法及装置
CN104014076A (zh) * 2014-04-21 2014-09-03 广州三瑞医疗器械有限公司 一种镇痛波形检测装置及其方法
CN104014076B (zh) * 2014-04-21 2016-02-24 广州三瑞医疗器械有限公司 一种镇痛波形检测装置及其方法
CN104330755A (zh) * 2014-08-26 2015-02-04 南通大学 一种开关磁阻电机磁链特性在线测量方法
CN104330755B (zh) * 2014-08-26 2017-03-29 南通大学 一种开关磁阻电机磁链特性在线测量方法
CN104539212A (zh) * 2015-01-20 2015-04-22 中国矿业大学 一种开关磁阻电机磁链特性在线测量方法
CN106441664A (zh) * 2016-11-21 2017-02-22 南通大学 一种开关磁阻电机转矩测量仪
CN106441664B (zh) * 2016-11-21 2019-04-19 南通大学 一种开关磁阻电机转矩测量仪
US10079566B1 (en) 2017-07-14 2018-09-18 Caterpillar Inc. Pulse injection current pattern for rotor position estimation
CN109309462A (zh) * 2018-05-10 2019-02-05 烟台仙崴机电有限公司 一种能源互联网新型发电机系统效率优化控制方法
CN111190128A (zh) * 2018-11-15 2020-05-22 北京自动化控制设备研究所 一种磁阻电机铁磁材料bh特性检测算法
CN111190128B (zh) * 2018-11-15 2022-10-18 北京自动化控制设备研究所 一种磁阻电机铁磁材料bh特性检测算法
CN113676107A (zh) * 2021-08-24 2021-11-19 山东交通学院 一种开关磁阻电机转子位置检测装置原点校准方法及系统
CN113676107B (zh) * 2021-08-24 2023-06-20 山东交通学院 一种开关磁阻电机转子位置检测装置原点校准方法及系统

Also Published As

Publication number Publication date
CN101769998B (zh) 2012-07-18

Similar Documents

Publication Publication Date Title
CN101769998B (zh) 开关磁阻电机磁链特性测量方法
CN105915141B (zh) 永磁同步电机永磁体磁链在线测量系统与测量方法
CN106655941B (zh) 一种内嵌式永磁同步电机的参数估计方法及装置
CN102931903B (zh) 一种永磁同步电动机转子初始位置检测方法
CN103018577B (zh) 一种测量永磁同步电机参数的方法
CN101604946B (zh) 一种适合高速的开关磁阻电机无位置传感器控制方法
CN104836506A (zh) 一种永磁同步电机转子零位校正系统及方法
CN104579083A (zh) 永磁同步电机矢量控制方法及系统
CN101630938A (zh) 无位置传感器永磁同步电机转子初始位置辨识方法
CN107659237A (zh) 一种永磁同步电机的无模型无差拍电流预测控制装置及其方法
CN103155399B (zh) Pm电动机的电流控制增益调整方法、电流控制方法以及控制装置
CN103869170A (zh) 电机中电阻的估算
CN110336513B (zh) 一种无位置传感器开关磁阻电机初始定位方法
CN104316876B (zh) 一种考虑互感耦合的快速获取三相12/8极srm磁链特性的方法
CN106441664B (zh) 一种开关磁阻电机转矩测量仪
CN101526823B (zh) 开关磁阻电机恒定转矩的控制方法
CN101533071A (zh) 一种开关磁阻电机全域非线性电感测试方法
Cai et al. An inductive position sensor with switched reluctance motor structure
CN110912482B (zh) 一种永磁同步电机磁链及温度估计方法
CN113067505B (zh) 在永磁同步电机控制过程中对电压矢量进行补偿的方法
CN104330755B (zh) 一种开关磁阻电机磁链特性在线测量方法
CN101707466A (zh) 一种步进电机、其丢步检测结构及其检测丢步的方法
CN110798111B (zh) 永磁同步电机旋转变压器零位检测方法、装置
CN106169895A (zh) 一种永磁直线同步电机电气参数测量方法
CN105406777A (zh) 一种永磁同步电机定子磁链的检测装置及检测方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
C17 Cessation of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20120718

Termination date: 20140118