CN112688616B - 开关磁阻电机相电流估测方法 - Google Patents

开关磁阻电机相电流估测方法 Download PDF

Info

Publication number
CN112688616B
CN112688616B CN202110107841.2A CN202110107841A CN112688616B CN 112688616 B CN112688616 B CN 112688616B CN 202110107841 A CN202110107841 A CN 202110107841A CN 112688616 B CN112688616 B CN 112688616B
Authority
CN
China
Prior art keywords
phase
current
reluctance motor
phase current
acc
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
CN202110107841.2A
Other languages
English (en)
Other versions
CN112688616A (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.)
RONGCHENG RONGJIA POWER CO Ltd
China University of Mining and Technology CUMT
Nanjing University of Aeronautics and Astronautics
Jiangsu University
Original Assignee
RONGCHENG RONGJIA POWER CO Ltd
China University of Mining and Technology CUMT
Nanjing University of Aeronautics and Astronautics
Jiangsu 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 RONGCHENG RONGJIA POWER CO Ltd, China University of Mining and Technology CUMT, Nanjing University of Aeronautics and Astronautics, Jiangsu University filed Critical RONGCHENG RONGJIA POWER CO Ltd
Priority to CN202110107841.2A priority Critical patent/CN112688616B/zh
Publication of CN112688616A publication Critical patent/CN112688616A/zh
Application granted granted Critical
Publication of CN112688616B publication Critical patent/CN112688616B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Control Of Electric Motors In General (AREA)

Abstract

由于目前开关磁阻电机不存在可观测的电流模型,因此无法采用闭环观测方法估测电流。本发明提出了一种开关磁阻电机相电流估测方法。该方法建立的相电流模型不包含微分方程,因此不需要闭环校正。该相电流模型具有结构简单和精度高的特点,保证了相电流的计算实时性和估测精度。由于相电流估测方法减少了开关磁阻电机系统中电流传感器的数量,因此降低了系统成本并提高了系统可靠性。

Description

开关磁阻电机相电流估测方法
技术领域
本发明涉及开关磁阻电机传感器领域。
背景技术
在控制系统中,传感器是获得时变参数的常用手段。然而,由于传感器的使用会导致系统成本的上升以及可靠性的降低。各种代替传感器的时变参数估测方法不断被提出。然而,到目前为止,对于开关磁阻电机,关于电流的可观测模型尚未被找到。因此,传统的闭环观测器方法无法用于开关磁阻电机的电流估测。为降低开关磁阻电机系统的成本及提高其可靠性,有必要研究一种适合于开关磁阻电机的电流估测方法。
发明内容
本发明提出的开关磁阻电机相电流估测方法如下。
一种开关磁阻电机相电流估测方法,其特征在于:
该开关磁阻电机相电流估测方法对应的开关磁阻电机上安装有一个测量相电流总和的电流传感器。开关磁阻电机的功率变换器采用不对称半桥功率变换器或其它适用于开关磁阻电机的功率变换器。所述相电流估测方法依次包括以下步骤:
1、建立相电流的计算模型;
2、建立开关磁阻电机的电磁转矩解析表达式;
3、对所述步骤2中的电磁转矩解析表达式进行误差补偿;
4、确定所述相电流计算模型的在线计算方式。
所述的相电流计算模型如下:
Figure GDA0004097180610000011
其中,T(ic1,θc1)为前一相导通绕组电流产生的电磁转矩,T(ic2,θc2)为后一相导通绕组电流产生的电磁转矩,ic1为前一相导通绕组电流值,ic2为后一相导通绕组电流值,θc1为前一导通相的相对转子位置,θc2为后一导通相的相对转子位置,相对转子位置为转子相对定子的位置,某一相的相对转子位置即转子相对该相定子的位置。TL为负载转矩,J为转动惯量,D为粘性摩擦系数,acc为电机转子的角加速度,ω为电机转子的角速度,is为所述电流传感器的电流测量值,该电流测量值为相电流总和。
采用分段建模方法建立所述电磁转矩解析表达式,依相对转子位置进行分段,在每个相对转子位置区间的电磁转矩解析表达式如下:
Figure GDA0004097180610000012
其中,T(i,θ)为电磁转矩,i为相电流值,θ为相对转子位置,Dx(θ)为关于相对转子位置的给定表达式,Cx(i)为关于相电流的待定表达式。Cx(i)的求解方法如下:
Figure GDA0004097180610000013
其中,θy为某一固定的相对转子位置,Ty(i,θy)为θy处的电磁转矩表达式。式(4)满足如下两个条件:①(Dx1),Dx2)…Dxn))为线性无关向量组;②m=n。求解式(4),得到Cx(i)的表达式,由此得到电磁转矩解析表达式。
由于所述的电磁转矩解析表达式存在误差,因此
T(im1,θc1)+T(im2,θc2)-TL≠J×acc+D×ω (5)
其中,im1为前一相导通绕组测量电流值,im2为后一相导通绕组测量电流值。设E(acc,ω)为误差补偿函数,使
T(im1,θc1)+T(im2,θc2)-TL=J×acc+D×ω+E(acc,ω) (6)
定义参数comp=T(im1,θc1)+T(im2,θc2)-TL=J×acc+D×ω+E(acc,ω)。在不同的acc和不同的ω下计算T(im1,θc1)+T(im2,θc2)-TL获得comp的值,建立表格comp—acc-ω,在线计算中,通过acc和ω查找对应的comp,comp包含E(acc,ω),因此对电磁转矩解析模型具有补偿作用。此时,待求解的相电流计算模型的式(1)转化为:
T(ic1,θc1)+T(ic2,θc2)-TL=comp (7)
将式(2)代入式(7),得
T(ic1,θc1)+T(is-ic1,θc2)-TL=comp (8)
当只有一相绕组被导通时,该相绕组的电流值为所述电流传感器的电流测量值。当有两相绕组被导通时,采用牛顿迭代法求解式(8),得到ic1的数值解,再根据式(2)求得ic2。至此相电流的估测完成。
其中,Dx(θ)选择为θx-1、sin(x×θ)、cos((x-1)×θ)、tan(x×θ)、或e(x-1)×θ其中之一。
其中,Ty(i,θy)通过在θy处的T-i测量数据拟合得到,T为电磁转矩值。拟合函数采用多项式函数。
当所述Ty(i,θy)的拟合多项式函数的最高次项的次数不大于3时,式(8)的求解方法还可以采用公式法。
其中,牛顿迭代法的迭代终止条件为下列二者之一:1、设ic1(n)为ic1的第n次迭代值,ic1(n+1)为ic1的第n+1次迭代值,ε为给定正数,当|ic1(n+1)-ic1(n)|≤ε时,停止迭代;2、设e(n)=|ic1(n)-im1|,给定误差阈值ethreshold,当|e(n+1)-e(n)|<ethreshold时,停止迭代,计算一个转子周期的n的平均值,将该平均值向上取整后作为在线计算的预设迭代次数。
给定电流估测阈值,当所述ic1的数值解小于电流估测阈值时,认为前一相的相电流下降为零,此时只有一相绕组被导通。
本发明的相电流估测方法减少了开关磁阻电机系统中电流传感器的数量,降低了系统成本并提高了系统可靠性。
附图说明
图1为本发明所述的开关磁阻电机相电流估测方法的步骤。
图2为三相12/8开关磁阻电机的功率变换器以及电流传感器的安装位置。
图3为在电流滞环控制方式下,按照本发明所述的开关磁阻电机相电流估测方法的估测结果。
图4为在电压斩波控制方式下,按照本发明所述的开关磁阻电机相电流估测方法的估测结果。
具体实施方式
本实施例的电机为三相12/8开关磁阻电机。其功率变换器采用不对称半桥功率变换器,电流传感器采用LEM传感器,三相绕组的引出线穿过LEM传感器中间的孔,如图1所示。
定义某一相最小电感处的相对转子位置为0°,那么该相最大电感处的相对转子位置为22.5°。根据在相对转子位置0°、1°、2°、3°、4°、5°、6°、7°、8°、9°、10°、11°、12°、13°、14°、15°、16°、17°、18°、19°、20°、21°、22°、22.5°、23°、24°、25°、26°、27°、28°、29°、30°、31°、32°、33°、34°、35°、36°、37°、38°、39°、40°、41°、42°、43°、44°、45°处的T-i测量数据分别拟合得到对应的Ty(i,θy)。Ty(i,θy)的拟合函数采用三次多项式函数。Dx(θ)选择为θx-1。利用式(4)分别在相对转子位置区间(0°,1°)、(1°,2°)、(2°,3°)、(3°,4°)、(4°,5°)、(5°,6°)、(6°,7°)、(7°,8°)、(8°,9°)、(9°,10°)、(10°,11°)、(11°,12°)、(12°,13°)、(13°,14°)、(14°,15°)、(15°,16°)、(16°,17°)、(17°,18°)、(18°,19°)、(19°,20°)、(20°,21°)、(21°,22°)、(22°,22.5°)、(22.5°,23°)(23°,24°)、(24°,25°)、(25°,26°)、(26°,27°)、(27°,28°)、(28°,29°)、(29°,30°)、(30°,31°)、(31°,32°)、(32°,33°)、(33°,34°)、(34°,35°)、(35°,36°)、(36°,37°)、(37°,38°)、(38°,39°)、(39°,40°)、(40°,41°)、(41°,42°)、(42°,43°)、(43°,44°)、(44°,45°)上求解Cx(i)并得到其电磁转矩表达式。由此得到的电磁转矩表达式为
T(i,θ)=f1(θ)×i3+f2(θ)×i2+f3(θ)×i+f4(θ) (9)
其中,f1(θ)、f2(θ)、f3(θ)、f4(θ)均为关于θ的函数。
当只有一相绕组被导通时,该相绕组的电流值为传感器的电流测量值。当有两相绕组被导通时,采用牛顿迭代法求解式(8),得到ic1的数值解,设
F(ic1)=(f1c1)-f1c2))×ic1 3+(f2c1)+f2c2)+3×f1c2)×is)×ic1 2+(-3×f1c2)×is 2-2×f2c2)×is+f3c1)-f3c2))×ic1+f1c2)×is 3+f2c2)×is 2+f3c2)×is+f4c1)+f4c2)-TL—comp (10)其迭代方程为
ic1(n+1)=ic1(n)-F(ic1)/F′(ic1) (11)
其中,F′(ic1)为F(ic1)的导函数。迭代终止条件选择为所述迭代终止条件2。给定误差阈值ethreshold设为0.1A。由此得到的预设迭代次数为2。其迭代结果为ic1的数值解,再根据式(2)求得ic2
给定电流估测阈值设为0.1A。
给定目标转速为600转/min,给定负载转矩为0.2N·m,分别在电流滞环控制方式和电压斩波控制方式下验证了本发明的相电流估测效果,见说明书附图3-4。

Claims (6)

1.一种开关磁阻电机相电流估测方法,其特征在于:该开关磁阻电机相电流估测方法对应的开关磁阻电机上安装有一个测量相电流总和的电流传感器;所述相电流估测方法依次包括以下步骤:
步骤1、建立相电流的计算模型;
步骤2、建立开关磁阻电机的电磁转矩解析表达式;
步骤3、对所述步骤2中的电磁转矩解析表达式进行误差补偿;
步骤4、确定所述相电流计算模型的在线计算方式;
所述的相电流计算模型为:
T(ic1,θc1)+T(ic2,θc2)-TL=J×acc+D×ω (1)
ic1+ic2=is (2)
其中,T(ic1,θc1)为前一相导通绕组电流产生的电磁转矩,T(ic2,θc2)为后一相导通绕组电流产生的电磁转矩,ic1为前一相导通绕组电流值,ic2为后一相导通绕组电流值,θc1为前一导通相的相对转子位置,θc2为后一导通相的相对转子位置,相对转子位置为转子相对定子的位置,某一相的相对转子位置即转子相对该相定子的位置;TL为负载转矩,J为转动惯量,D为粘性摩擦系数,acc为电机转子的角加速度,ω为电机转子的角速度,is为所述电流传感器的电流测量值,该电流测量值为相电流总和;
所述电磁转矩解析表达式采用分段建模方法建立,依相对转子位置进行分段,在每个相对转子位置区间的电磁转矩解析表达式如下:
Figure FDA0004151095160000011
其中,T(i,θ)为电磁转矩,i为相电流值,θ为相对转子位置,Dx(θ)为关于相对转子位置的给定表达式,Cx(i)为关于相电流的待定表达式;Cx(i)的求解方法如下:
Figure FDA0004151095160000012
其中,θy为某一固定的相对转子位置,Ty(i,θy)为θy处的电磁转矩表达式;式(4)满足如下两个条件:①(Dx1),Dx2)...Dxn))为线性无关向量组;②m=n;求解式(4),得到Cx(i)的表达式,由此得到电磁转矩解析表达式;
由于所述的电磁转矩解析表达式存在误差,因此
T(im1,θc1)+T(im2,θc2)-TL≠J×acc+D×ω (5)
其中,im1为前一相导通绕组测量电流值,im2为后一相导通绕组测量电流值;设E(acc,ω)为误差补偿函数,使
T(im1,θc1)+T(im2,θc2)-TL=J×acc+D×ω+E(acc,ω) (6)
定义参数comp=T(im1,θc1)+T(im2,θc2)-TL=J×acc+D×ω+E(acc,ω);在不同的acc和不同的ω下计算T(im1,θc1)+T(im2,θc2)-TL获得comp的值,建立表格comp-acc-ω,在线计算中,通过acc和ω查找对应的comp,comp包含E(acc,ω),因此对电磁转矩解析模型具有补偿作用;此时,待求解的相电流计算模型的式(1)转化为:
T(ic1,θc1)+T(ic2,θc2)-TL=comp (7)
将式(2)代入式(7),得
T(ic1,θc1)+T(is-ic1,θc2)-TL=comp (8)
当只有一相绕组被导通时,该相绕组的电流值为所述电流传感器的电流测量值;当有两相绕组被导通时,采用牛顿迭代法求解式(8),得到ic1的数值解,再根据式(2)求得ic2
2.根据权利要求1所述的开关磁阻电机相电流估测方法,其中,Dx(θ)选择为θx-1、sin(x×θ)、cos((x-1)×θ)、tan(x×θ)、或e(x-1)×θ其中之一。
3.根据权利要求2所述的开关磁阻电机相电流估测方法,其中,Ty(i,θy)通过在θy处的T-i测量数据拟合得到,T为电磁转矩值,拟合函数采用多项式函数。
4.根据权利要求3所述的开关磁阻电机相电流估测方法,其中,当所述Ty(i,θy)的拟合多项式函数的最高次项的次数不大于3时,式(8)的求解方法还可以采用公式法。
5.根据权利要求1-4之一所述的开关磁阻电机相电流估测方法,其中,牛顿迭代法的迭代终止条件为下列二者之一:1、设ic1(n)为ic1的第n次迭代值,ic1(n+1)为ic1的第n+1次迭代值,ε为给定正数,当|ic1(n+1)-ic1(n)|≤ε时,停止迭代;2、设e(n)=|ic1(n)-im1|,给定误差阈值ethreshold,当|e(n+1)-e(n)|<ethreshold时,停止迭代,计算一个转子周期的n的平均值,将该平均值向上取整后作为在线计算的预设迭代次数。
6.根据权利要求5所述的开关磁阻电机相电流估测方法,给定电流估测阈值,当所述ic1的数值解小于电流估测阈值时,认为前一相的相电流下降为零,此时只有一相绕组被导通。
CN202110107841.2A 2021-01-27 2021-01-27 开关磁阻电机相电流估测方法 Active CN112688616B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110107841.2A CN112688616B (zh) 2021-01-27 2021-01-27 开关磁阻电机相电流估测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110107841.2A CN112688616B (zh) 2021-01-27 2021-01-27 开关磁阻电机相电流估测方法

Publications (2)

Publication Number Publication Date
CN112688616A CN112688616A (zh) 2021-04-20
CN112688616B true CN112688616B (zh) 2023-07-04

Family

ID=75459319

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110107841.2A Active CN112688616B (zh) 2021-01-27 2021-01-27 开关磁阻电机相电流估测方法

Country Status (1)

Country Link
CN (1) CN112688616B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105811849A (zh) * 2016-05-06 2016-07-27 桂林电子科技大学 一种电流非线性补偿的开关磁阻电机转矩控制方法与系统
CN107453672A (zh) * 2017-08-01 2017-12-08 合肥工业大学 分段式谐波电流抑制开关磁阻电机转矩脉动的控制方法
CN107979255A (zh) * 2016-10-24 2018-05-01 南京理工大学 一种大转矩低脉动的双定子轴向磁通开关磁阻电机
CN108429509A (zh) * 2018-03-15 2018-08-21 合肥工业大学 开关磁阻电机分段pwm占空比控制方法
CN110417328A (zh) * 2019-06-26 2019-11-05 湘潭大学 一种基于ckmtoa的srm转矩分配函数优化补偿控制策略

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105811849A (zh) * 2016-05-06 2016-07-27 桂林电子科技大学 一种电流非线性补偿的开关磁阻电机转矩控制方法与系统
CN107979255A (zh) * 2016-10-24 2018-05-01 南京理工大学 一种大转矩低脉动的双定子轴向磁通开关磁阻电机
CN107453672A (zh) * 2017-08-01 2017-12-08 合肥工业大学 分段式谐波电流抑制开关磁阻电机转矩脉动的控制方法
CN108429509A (zh) * 2018-03-15 2018-08-21 合肥工业大学 开关磁阻电机分段pwm占空比控制方法
CN110417328A (zh) * 2019-06-26 2019-11-05 湘潭大学 一种基于ckmtoa的srm转矩分配函数优化补偿控制策略

Also Published As

Publication number Publication date
CN112688616A (zh) 2021-04-20

Similar Documents

Publication Publication Date Title
CN104218863B (zh) 使用状态观测器确定同步电机中内部温度的方法和系统
JP6333772B2 (ja) 同期電動機の温度推定装置
CN103825525B (zh) 一种改进的无传感器永磁同步电机速度估测方法
EP2258043B1 (en) Sensorless control of salient-pole machines
JP2014204668A (ja) 誘導電動機のセンサレスベクトル制御装置
CN110165960B (zh) 偏差角度估计方法、系统、计算机装置及介质
US6509711B1 (en) Digital rotor flux observer
CN111293932A (zh) 控制永磁体同步电机的方法以及电机电路
JP2016518105A (ja) 多相回転電機の回転子の角度位置を推定するための方法、および、そのような機械のための多相インバータの制御への応用
JP6557452B2 (ja) 状態観測器を使用して同期電動装置の中の回転子の位置および速度を求める方法
JP3896047B2 (ja) モータ駆動制御装置
US7774148B2 (en) Torque estimator for IPM motors
KR101121474B1 (ko) Dc 서보모터 시스템의 관성모멘트와 마찰계수를 구하는 방법
CN112688616B (zh) 开关磁阻电机相电流估测方法
JP2002051595A (ja) 電動機の制御装置
CN109699197B (zh) 用于估算机动车辆的交流电机的转子的位置和速度的方法及相应系统
EP3477846B1 (en) Method for determining a measuring offset of a rotor position sensor, controller unit for an electric machine and electric machine for a vehicle
CN108649849A (zh) 一种简单的无传感器永磁同步电机速度估测方法
JP6108114B2 (ja) 永久磁石形同期電動機の制御装置
JP4061446B2 (ja) 同期電動機の抵抗値同定方法とその制御装置
CN114070146A (zh) 一种故障检测方法、装置、设备及存储介质
JP4281054B2 (ja) 誘導電動機の速度センサレスベクトル制御方法及び装置
US11936316B2 (en) Position estimation for permanent magnet synchronous machines via rotor flux space vector identification
JP2023169905A (ja) Pmモータの温度推定装置
CN108521245A (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