CN113489409A - 一种基于旋转综合矢量的脉振高频注入的永磁电机的位置估计方法 - Google Patents

一种基于旋转综合矢量的脉振高频注入的永磁电机的位置估计方法 Download PDF

Info

Publication number
CN113489409A
CN113489409A CN202110675485.4A CN202110675485A CN113489409A CN 113489409 A CN113489409 A CN 113489409A CN 202110675485 A CN202110675485 A CN 202110675485A CN 113489409 A CN113489409 A CN 113489409A
Authority
CN
China
Prior art keywords
axis
motor
frequency
generated
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
CN202110675485.4A
Other languages
English (en)
Other versions
CN113489409B (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.)
Nanjing University of Aeronautics and Astronautics
Original Assignee
Nanjing University of Aeronautics and Astronautics
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 Nanjing University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN202110675485.4A priority Critical patent/CN113489409B/zh
Publication of CN113489409A publication Critical patent/CN113489409A/zh
Application granted granted Critical
Publication of CN113489409B publication Critical patent/CN113489409B/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/0003Control strategies in general, e.g. linear type, e.g. P, PI, PID, using robust control
    • 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/05Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation specially adapted for damping motor oscillations, e.g. for reducing hunting
    • 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/24Vector control not involving the use of rotor position or rotor speed sensors
    • 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
    • H02P25/00Arrangements or methods for the control of AC motors characterised by the kind of AC motor or by structural details
    • H02P25/02Arrangements or methods for the control of AC motors characterised by the kind of AC motor or by structural details characterised by the kind of motor
    • H02P25/022Synchronous motors
    • 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
    • H02P2203/00Indexing scheme relating to controlling arrangements characterised by the means for detecting the position of the rotor
    • H02P2203/11Determination or estimation of the rotor position or other motor parameters based on the analysis of high-frequency signals
    • 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
    • H02P2207/00Indexing scheme relating to controlling arrangements characterised by the type of motor
    • H02P2207/05Synchronous machines, e.g. with permanent magnets or DC excitation

Landscapes

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

Abstract

本发明公开了一种基于旋转综合矢量的脉振高频注入的永磁电机的位置估计方法,首先设置旋转电压矢量,在估计的d轴和q轴上注入该旋转电压矢量的分量,然后求取各自的零序电压并求和;对求和电压进行变换后输入至低通滤波器;最后将滤波结果输入至PI环节和积分环节,即可获得电机的估计位置;本发明通过在估计的d轴和q轴上注入高频信号,使得进入PI调节器的信号中不存在六倍频干扰信号,解决了传统高频注入方法中存在的六倍频干扰问题;其次,高频信号产生的旋转矢量是长度固定的且相对d轴旋转频率为ωh,解决了逆二倍频坐标算法中存在的转矩脉动谐波多、存在三倍频转矩脉动、转速波动大的问题。同时提高了转子位置辨识的精度,减小了转速的波动。

Description

一种基于旋转综合矢量的脉振高频注入的永磁电机的位置估 计方法
技术领域
本发明涉及电机控制技术领域,主要涉及一种基于旋转综合矢量的脉振高频注入的永磁电机的位置估计方法。
背景技术
永磁同步电机驱动系统广泛应用于航空航天、交通运输、家用电器等领域。而采用无位置技术的的永磁电机驱动系统可以减小系统的体积与成本。
经典的永磁同步电机的无位置传感器技术主要有滑模观测器法和高频注入法。滑模观测器法是通过构造一个滑模面,然后设计滑模控制律,进而得到扩展反电动势,实现位置信息的观测。滑模观测器设计的根本目标是得到扩展反电动势,从而得到位置信息。电机在零速情况下没有扩展反电动势,低速情况下反电动势谐波成分含量大,因此滑模等一系列通过反电动势来观测位置信息的方法,在零低速情况下,均无法准确观测出电机的位置信息。
高频信号注入法主要利用电机自身的非理想特性对电机的转速和位置信号进行估计。主要原理是在电机定子侧注入高频电压或者高频电流信号,注入的高频信号在电机内会形成高频磁场。由于电机具有非理想特性,例如转子结构凸极性、转子饱和凸极性等,这种非理想特性会对高频磁场产生调制作用,使得电机定子侧产生和电机位置和转速相关的高频电压或电流信号,通过提取有效信号和解调制等方式得到电机的转速和位置信息。由于没有用到反电动势等受到速度约束的物理量,因此在零低速情况下也有较好的性能。脉振高频注入法通过在同步旋转坐标系的直轴d轴中注入高频正弦电压信号,注入的信号会在静止坐标系中形成一个高频的脉振电压信号,通过对交轴高频电流信号进行幅值调制后,可以提取出其中与转子位置有关的信息,以此方式来获得转子的位置速度信息。
但是,对于脉振高频注入法,在计算转子位置的实际值和估计值的误差的过程中,会引入一个六倍频扰动量(该扰动量的频率为电机电角频率的六倍),从而影响了转子位置的辨识精度。
为了解决六倍频扰动量的干扰问题,文献《P.L.Xu,and Z.Q.Zhu,Novel CarrierSignal Injection Method Using Zero-Sequence Voltage for Sensorless Control ofPMSM Drives,IEEE TRANSACTIONS ON INDUSTRIAL ELECTRONICS,VOL.63,NO.4,APRIL,2016》提出了一种逆二倍频同步坐标的脉振高频注入法,该方法将高频信号注入在了逆二倍频同步坐标的q轴上,可以使得六倍频扰动量得到很好的抑制。然而,该方法逆二倍频同步坐标与电机实际的旋转坐标逆向旋转,导致注入的高频信号在电机中产生了两种不同频率的转矩脉动,给转矩脉动的抑制和转速控制器的设计带来了麻烦,增加了转矩脉动和转速波动。
发明内容
发明目的:本发明提供了一种综合矢量旋转式的脉振高频注入的永磁电机的位置估计方法,在抑制了六倍频扰动量的干扰问题的同时,减少了转矩脉动的频率数量,解决了现有逆二倍频同步坐标的脉振高频注入法中存在多种频率的转矩脉动、转矩脉动大、转速波动大的问题。
技术方案:为实现上述目的,本发明采用的技术方案为:
一种基于旋转综合矢量的脉振高频注入的永磁电机的位置估计方法,所述永磁电机的转速控制方法如下:
将给定电机电角频率ω*与估计电角频率
Figure BDA00031208362500000223
的差值
Figure BDA0003120836250000021
输入至第一PI环节,获得电机的给定q轴电流
Figure BDA0003120836250000022
通过电流传感器获取电机的三相电流ia,ib和ic,对ia,ib和ic进行abc/dq变换,获得实际q轴电流iq和实际d轴电流id;将给定q轴电流
Figure BDA0003120836250000023
和实际q轴电流iq的差值
Figure BDA0003120836250000024
输入至第二PI环节,得到
Figure BDA0003120836250000025
设定给定d轴电流
Figure BDA0003120836250000026
将给定d轴电流
Figure BDA0003120836250000027
和实际d轴电流id的差值
Figure BDA0003120836250000028
输入至第二PI环节,得到
Figure BDA0003120836250000029
Figure BDA00031208362500000210
与高频电压Uhsinωht相加获得
Figure BDA00031208362500000211
Figure BDA00031208362500000212
与高频电压Uhcosωht相加得到
Figure BDA00031208362500000213
分别进行dq/abc变换,得到对应逆变器三相桥臂的占空比;其特征在于,采用基于旋转综合矢量的脉振高频注入法估计电机转子位置,具体步骤如下:
步骤S1、设置旋转电压矢量
Figure BDA00031208362500000214
其中Uh为旋转电压矢量的长度,
Figure BDA00031208362500000215
为旋转电压矢量的相位角,ωh为旋转电压矢量的角频率,
Figure BDA00031208362500000216
为估计的转子位置角;
步骤S2、电机的α轴为电机A相绕组的轴线方向,将估计的d轴记为
Figure BDA00031208362500000217
轴,且
Figure BDA00031208362500000218
轴与α轴的夹角为
Figure BDA00031208362500000219
将估计的q轴记为
Figure BDA00031208362500000220
轴,
Figure BDA00031208362500000221
轴与α轴的夹角为
Figure BDA00031208362500000222
所述旋转电压矢量
Figure BDA0003120836250000031
Figure BDA0003120836250000032
轴上的分量为Uh cosωht,在
Figure BDA0003120836250000033
轴上的分量为Uhsinωht
步骤S3、将Uh cosωht注入到电机的
Figure BDA0003120836250000034
轴上,求取Uh cosωht产生的零序电压:
Figure BDA0003120836250000035
其中,L2为电机自感的二次谐波的幅值,L0为电机自感的平均值,M0为电机互感的平均值,Ld为电机d轴电感,Lq为电机q轴电感;
将Uh sinωht注入到电机的
Figure BDA0003120836250000036
轴上,求取Uh sinωht产生的零序电压:
Figure BDA0003120836250000037
其中
Figure BDA0003120836250000038
Figure BDA0003120836250000039
轴与电机的实际d轴夹角;
将Uhcosωht产生的零序电压和Uh sinωht产生的零序电压求和得到:
Figure BDA00031208362500000310
步骤S4、对求和电压进行变换如下:
Figure BDA00031208362500000311
将变换结果输入至低通滤波器进行滤波可得:
Figure BDA00031208362500000312
其中k为滤波系数;
步骤S5、将步骤S4所述低通滤波结果输入至PI环节,获得电机的估计电角频率,将所述估计电角频率输入至积分环节即可获得电机的估计位置。
进一步地,所述步骤S3中Uh cosωht产生的零序电压求取过程如下:
步骤S3.1、记
Figure BDA00031208362500000313
轴与电机的实际d轴夹角为
Figure BDA00031208362500000314
电机的实际d轴与α轴的夹角为θ,则
Figure BDA00031208362500000315
Uh cosωht产生的dq轴电流的变化率为:
Figure BDA0003120836250000041
步骤S3.2,分别计算Uh cosωht产生的三相电流的变化率为:
Figure BDA0003120836250000042
Figure BDA0003120836250000043
Figure BDA0003120836250000044
步骤S3.3,分别计算Uh cosωht产生的三相电压为:
Figure BDA0003120836250000045
Figure BDA0003120836250000046
Figure BDA0003120836250000047
步骤S3.4,计算Uh cosωht产生的零序电压如下:
Figure BDA0003120836250000048
Uh sinωht产生的零序电压求取过程如下:
步骤L3.1、记
Figure BDA0003120836250000049
轴与电机的实际d轴夹角为
Figure BDA00031208362500000410
电机的实际d轴与α轴的夹角为θ,则
Figure BDA0003120836250000051
Uh sinωht产生的dq轴电流的变化率为:
Figure BDA0003120836250000052
步骤L3.2、分别计算Uh sinωht产生的三相电流的变化率如下:
Figure BDA0003120836250000053
Figure BDA0003120836250000054
Figure BDA0003120836250000055
步骤L3.3、分别计算Uh sinωht产生的三相电压为:
Figure BDA0003120836250000056
Figure BDA0003120836250000057
Figure BDA0003120836250000058
步骤L3.4、计算Uh sinωht产生的零序电压为:
Figure BDA0003120836250000059
有益效果:
(1)、本发明提出的方案通过在估计的d轴和q轴上注入高频信号,使得进入PI调节器的信号中不存在六倍频干扰信号,解决了传统高频注入方法中存在的六倍频干扰问题;
(2)、本发明提出的方案中,高频信号产生的旋转矢量是长度固定的且相对d轴旋转频率为ωh,有效地解决了逆二倍频坐标算法中存在的转矩脉动谐波多、存在三倍频转矩脉动、转速波动大的问题;
(3)、本发明提出的方案可以有效提高转子位置辨识的精度,减小了转子转速的波动。
附图说明
图1为本发明永磁同步电机的转速控制框图;
图2为传统的基于脉振高频注入的无位置方法控制框图;
图3为传统的基于脉振高频注入的无位置方法的关键波形之一——转子位置及误差;
图4为传统的基于脉振高频注入的无位置方法的关键波形之二——转子位置误差的频谱;
图5为传统的基于脉振高频注入的无位置方法的关键波形之三——转速;
图6为传统的基于脉振高频注入的无位置方法的关键波形之四——转速误差的频谱;
图7为文献“P.L.Xu and Z.Q.Zhu,Novel Carrier Signal Injection MethodUsing Zero-Sequence Voltage for Sensorless Control of PMSM Drives,IEEETRANSACTIONS ON INDUSTRIAL ELECTRONICS,63(4),pp:2053-2061,2016.”提出的基于逆二倍频坐标的脉振高频注入的无位置方法的控制框图;
图8为基于逆二倍频坐标的脉振高频注入的无位置方法的关键波形之一——转子位置及误差;
图9为基于逆二倍频坐标的脉振高频注入的无位置方法的关键波形之二——转子位置误差的频谱;
图10为基于逆二倍频坐标的脉振高频注入的无位置方法的关键波形之三——转矩;
图11为基于逆二倍频坐标的脉振高频注入的无位置方法的关键波形之四——转速;
图12为基于逆二倍频坐标的脉振高频注入的无位置方法的关键波形之五——转速误差的频谱;
图13为基于逆二倍频坐标的脉振高频注入的无位置方法的关键波形之六——高频信号所产生的转矩的频谱;
图14为本发明提出的基于旋转综合矢量的脉振高频注入的永磁电机的位置估计方法的控制框图;
图15为本发明提出的基于旋转综合矢量的脉振高频注入的永磁电机的位置估计方法的关键波形之一——转子位置及误差;
图16为本发明提出的基于旋转综合矢量的脉振高频注入的永磁电机的位置估计方法的关键波形之二——转子位置误差的频谱;
图17为本发明提出的基于旋转综合矢量的脉振高频注入的永磁电机的位置估计方法的关键波形之三——转矩;
图18为本发明提出的基于旋转综合矢量的脉振高频注入的永磁电机的位置估计方法的关键波形之四——转速;
图19为本发明提出的基于旋转综合矢量的脉振高频注入的永磁电机的位置估计方法的关键波形之五——转速误差的频谱;
图20为本发明提出的基于旋转综合矢量的脉振高频注入的永磁电机的位置估计方法的关键波形之六——高频信号所产生的转矩的频谱。
具体实施方式
下面结合附图对本发明作更进一步的说明。
如图1所示的永磁同步电机的转速控制方法,将给定电机电角频率ω*与估计电角频率
Figure BDA0003120836250000071
的差值
Figure BDA0003120836250000072
输入至第一PI环节,获得电机的给定q轴电流
Figure BDA0003120836250000073
通过电流传感器获取电机的三相电流ia,ib和ic,对ia,ib和ic进行abc/dq变换,获得实际q轴电流iq和实际d轴电流id;将给定q轴电流
Figure BDA0003120836250000074
和实际q轴电流iq的差值
Figure BDA0003120836250000075
输入至第二PI环节,得到
Figure BDA0003120836250000076
设定给定d轴电流
Figure BDA0003120836250000077
将给定d轴电流
Figure BDA0003120836250000078
和实际d轴电流id的差值
Figure BDA0003120836250000079
输入至第二PI环节,得到
Figure BDA00031208362500000710
Figure BDA00031208362500000711
与高频电压Uhsinωht相加获得
Figure BDA00031208362500000712
Figure BDA00031208362500000713
与高频电压Uhcosωht相加得到
Figure BDA00031208362500000714
分别进行dq/abc变换,得到对应逆变器三相桥臂的占空比。
传统基于脉振高频注入的无位置方法如图2所示,它的核心思想是在估计的d轴上注入Uh cosωht注入到电机的
Figure BDA0003120836250000081
轴上,Uh cosωht产生的零序电压求取过程如下:记
Figure BDA0003120836250000082
轴与电机的实际d轴夹角为
Figure BDA0003120836250000083
Figure BDA0003120836250000084
电机的实际d轴与α轴的夹角为θ,那么Uh cosωht产生的dq轴电流的变化率为:
Figure BDA0003120836250000085
分别计算Uh cosωht产生的三相电流的变化率如下:
Figure BDA0003120836250000086
Figure BDA0003120836250000087
Figure BDA0003120836250000088
分别计算Uh cosωht产生的三相电压如下:
Figure BDA0003120836250000089
Figure BDA00031208362500000810
Figure BDA00031208362500000811
Uh cosωht产生的零序电压为:
Figure BDA0003120836250000091
Figure BDA0003120836250000092
作如下变换:
Figure BDA0003120836250000093
Figure BDA0003120836250000094
经过低通滤波器得到
Figure BDA0003120836250000095
其中k为滤波系数;
将上述滤波结果经过积化和差,得到:
Figure BDA0003120836250000096
Figure BDA0003120836250000097
依次通过PI调节器和积分环节得到转子的辨识位置。
从以上推导可以看出,对于传统的基于脉振高频注入的无位置方法,它的核心思想是在估计的d轴上注入高频信号,该信号在估计的d轴上是一个脉振信号,它对应的综合矢量不是一个旋转量。在这样的注入方式下,
Figure BDA0003120836250000098
会伴随着
Figure BDA0003120836250000099
一同进入PI调节器,稳态运行时,
Figure BDA00031208362500000910
那么
Figure BDA00031208362500000911
所以在辨识的转子位置中会有一个六倍频的干扰信号(该干扰信号的频率为电机转角频率的六倍),如图3所示。图4为转子位置误差的频谱分析,可以看出,六次谐波的含量较高。转速的脉动的峰峰值达到了5rpm,如图5所示。
文献“P.L.Xu and Z.Q.Zhu,Novel Carrier Signal Injection Method UsingZero-Sequence Voltage for Sensorless Control of PMSM Drives,IEEE TRANSACTIONSON INDUSTRIAL ELECTRONICS,63(4),pp:2053-2061,2016.”提出的基于逆二倍频坐标的脉振高频注入的无位置方法的控制框图如7所示,它的主要步骤如下:
设置一个旋转的
Figure BDA0003120836250000101
坐标,该坐标中,
Figure BDA0003120836250000102
轴超前
Figure BDA0003120836250000103
轴90°,
Figure BDA0003120836250000104
轴的旋转频率为
Figure BDA0003120836250000105
轴旋转频率的2倍,且旋转方向相反,在
Figure BDA0003120836250000106
轴上注入信号Uh cosωht,Uh cosωht产生的零序电压求取过程如下:记
Figure BDA0003120836250000107
轴与电机α轴的夹角为
Figure BDA0003120836250000108
那么
Figure BDA0003120836250000109
轴与电机α轴的夹角为
Figure BDA00031208362500001010
电机的实际d轴与α轴的夹角为θ,那么Uh cosωht产生的dq轴电流的变化率为:
Figure BDA00031208362500001011
计算Uh cosωht产生的三相电流的变化率为:
Figure BDA00031208362500001012
Figure BDA00031208362500001013
Figure BDA00031208362500001014
计算Uh cosωht产生的三相电压为:
Figure BDA0003120836250000111
Figure BDA0003120836250000112
Figure BDA0003120836250000113
计算Uh cosωht产生的零序电压为:
Figure BDA0003120836250000114
Figure BDA0003120836250000115
乘以cosωht得到
Figure BDA0003120836250000116
Figure BDA0003120836250000117
经过低通滤波器得到
Figure BDA0003120836250000118
其中k为滤波系数;
Figure BDA0003120836250000119
依次通过PI调节器和积分环节得到转子的辨识位置。
可以看出,对于文献提出的算法,进入PI调节器的只有
Figure BDA00031208362500001110
没有六倍频的干扰,如图8所示。图9为图8中转子位置误差的频谱,可以看出,与图4相比,六次谐波的含量大大减小。
因此与图2中的传统方法相比,基于逆二倍频坐标的脉振高频注入的无位置方法提高了转子辨识精度。
但是,在该算法中,高频信号注入的轴的旋转速度与电机的d轴旋转不同步,所以注入的高频信号在电机q轴产生的电流如下所示:
Figure BDA00031208362500001111
该iq产生的转矩脉动如下所示:
Figure BDA0003120836250000121
其中,Pr为电机的转子极对数,ψpm为电机的永磁磁链。
从式中可以看出,转矩脉动中既包含了角频率为ωh的高频脉动转矩,也包含了角频率为
Figure BDA0003120836250000122
的低频脉动转矩,如图10所示,转矩波形为包络线性质,脉动频率的数量较多,且频率数值相差较大,导致转速调节器难以优化设计,使得转速波动较大,如图11所示,转速波动的峰峰值为3rpm。角频率
Figure BDA0003120836250000123
为电机同步角频率的3倍,所以转速误差和高频信号产生的转矩都包含了较高的3次谐波分量,如图12与13所示。
本发明提出了一种基于旋转综合矢量的脉振高频注入的永磁电机的位置估计方法,如图14所示,具体步骤如下:
步骤S1、设置旋转电压矢量
Figure BDA0003120836250000124
其中Uh为旋转电压矢量的长度,
Figure BDA0003120836250000125
为旋转电压矢量的相位角,ωh为旋转电压矢量的角频率,
Figure BDA0003120836250000126
为估计的转子位置角;
步骤S2、电机的α轴为电机A相绕组的轴线方向,将估计的d轴记为
Figure BDA0003120836250000127
轴,且
Figure BDA0003120836250000128
轴与α轴的夹角为
Figure BDA0003120836250000129
将估计的q轴记为
Figure BDA00031208362500001210
轴,
Figure BDA00031208362500001211
轴与α轴的夹角为
Figure BDA00031208362500001212
所述旋转电压矢量
Figure BDA00031208362500001213
Figure BDA00031208362500001214
轴上的分量为Uh cosωht,在
Figure BDA00031208362500001215
轴上的分量为Uhsinωht;
步骤S3、将Uh cosωht注入到电机的
Figure BDA00031208362500001216
轴上,求取Uh cosωht产生的零序电压。具体地,
步骤S3.1、记
Figure BDA00031208362500001217
轴与电机的实际d轴夹角为
Figure BDA00031208362500001218
电机的实际d轴与α轴的夹角为θ,则
Figure BDA00031208362500001219
Uh cosωht产生的dq轴电流的变化率为:
Figure BDA00031208362500001220
步骤S3.2,分别计算Uh cosωht产生的三相电流的变化率为:
Figure BDA0003120836250000131
Figure BDA0003120836250000132
Figure BDA0003120836250000133
步骤S3.3,分别计算Uh cosωht产生的三相电压为:
Figure BDA0003120836250000134
Figure BDA0003120836250000135
Figure BDA0003120836250000136
步骤S3.4,计算Uh cosωht产生的零序电压如下:
Figure BDA0003120836250000137
将Uh sinωht注入到电机的
Figure BDA0003120836250000138
轴上,求取Uh sinωht产生的零序电压。具体如下:
步骤L3.1、记
Figure BDA0003120836250000139
轴与电机的实际d轴夹角为
Figure BDA00031208362500001310
电机的实际d轴与α轴的夹角为θ,则
Figure BDA00031208362500001311
Uh sinωht产生的dq轴电流的变化率为:
Figure BDA0003120836250000141
步骤L3.2、分别计算Uh sinωht产生的三相电流的变化率如下:
Figure BDA0003120836250000142
Figure BDA0003120836250000143
Figure BDA0003120836250000144
步骤L3.3、分别计算Uh sinωht产生的三相电压为:
Figure BDA0003120836250000145
Figure BDA0003120836250000146
Figure BDA0003120836250000147
步骤L3.4、计算Uh sinωht产生的零序电压为:
Figure BDA0003120836250000148
其中
Figure BDA0003120836250000149
Figure BDA00031208362500001410
轴与电机的实际d轴夹角;
将Uhcosωht产生的零序电压和Uh sinωht产生的零序电压求和得到:
Figure BDA0003120836250000151
步骤S4、对求和电压进行变换如下:
Figure BDA0003120836250000152
将变换结果输入至低通滤波器进行滤波可得:
Figure BDA0003120836250000153
其中k为滤波系数;
步骤S5、将步骤S4所述低通滤波结果输入至PI环节,获得电机的估计电角频率,将所述估计电角频率输入至积分环节即可获得电机的估计位置。
可以看出,对于本发明的算法,进入PI调节器的只有
Figure BDA0003120836250000154
没有六倍频的干扰,如图15所示,因此与图2中的传统方法相比,本发明提出的方案提高了转子辨识精度。图16为本发明提出的方案的转子位置误差的频谱,可以看出,六次谐波的含量大大降低。
同时,本发明提出的方案中,高频信号注入在
Figure BDA0003120836250000155
轴和
Figure BDA0003120836250000156
轴上,且在
Figure BDA0003120836250000157
轴和
Figure BDA0003120836250000158
轴上的两个高频信号相差90°,注入的高频信号形成的综合矢量相对于d轴是一个稳定的旋转矢量,该旋转矢量相对于d轴的旋转角频率为ωh,该旋转矢量在q轴上产生的电流如下所示:
Figure BDA0003120836250000159
该电流产生的转矩脉动如下所示:
Figure BDA00031208362500001510
可以看出,式中只有频率为ωh的一个脉动转矩,如图17所示,相对于基于逆二倍频坐标的脉振高频注入的无位置方法,本发明产生的转速波动减小,转速波动的峰峰值只有1.4rpm,如图18所示。
同时可以看出,现有文献和本发明提出的方案都解决了六倍频的干扰问题。但是文献注入的信号还引入了3倍频的转矩脉动,体现为图12和图13中的3次谐波的含量相对较高,而本发明不存在这样的问题,体现为图中19和图20中3次谐波的含量减小。
综上所述,首先,本发明提出的方案通过在估计的d轴和q轴上注入高频信号,使得进入PI调节器的信号中不存在六倍频干扰信号,解决了传统高频注入方法中存在的六倍频干扰问题;其次,本发明提出的方案中,高频信号产生的旋转矢量是长度固定的且相对d轴旋转频率为ωh,有效地解决了逆二倍频坐标算法中存在的转矩脉动谐波多、存在三倍频转矩脉动、转速波动大的问题。本发明同时提高了转子位置辨识的精度和减小了转速的波动。
以上所述仅是本发明的优选实施方式,应当指出:对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。

Claims (2)

1.一种基于旋转综合矢量的脉振高频注入的永磁电机的位置估计方法,所述永磁电机的转速控制方法如下:
将给定电机电角频率ω*与估计电角频率
Figure FDA0003120836240000011
的差值
Figure FDA0003120836240000012
输入至第一PI环节,获得电机的给定q轴电流
Figure FDA0003120836240000013
通过电流传感器获取电机的三相电流ia,ib和ic,对ia,ib和ic进行abc/dq变换,获得实际q轴电流iq和实际d轴电流id;将给定q轴电流
Figure FDA0003120836240000014
和实际q轴电流iq的差值
Figure FDA0003120836240000015
输入至第二PI环节,得到
Figure FDA0003120836240000016
设定给定d轴电流
Figure FDA0003120836240000017
将给定d轴电流
Figure FDA0003120836240000018
和实际d轴电流id的差值
Figure FDA0003120836240000019
输入至第二PI环节,得到
Figure FDA00031208362400000110
Figure FDA00031208362400000111
与高频电压Uhsinωht相加获得
Figure FDA00031208362400000112
Figure FDA00031208362400000113
与高频电压Uhcosωht相加得到
Figure FDA00031208362400000114
分别进行dq/abc变换,得到对应逆变器三相桥臂的占空比;其特征在于,采用基于旋转综合矢量的脉振高频注入法估计电机转子位置,具体步骤如下:
步骤S1、设置旋转电压矢量
Figure FDA00031208362400000115
其中Uh为旋转电压矢量的长度,
Figure FDA00031208362400000116
为旋转电压矢量的相位角,ωh为旋转电压矢量的角频率,
Figure FDA00031208362400000117
为估计的转子位置角;
步骤S2、电机的α轴为电机A相绕组的轴线方向,将估计的d轴记为
Figure FDA00031208362400000118
轴,且
Figure FDA00031208362400000119
轴与α轴的夹角为
Figure FDA00031208362400000120
将估计的q轴记为
Figure FDA00031208362400000121
轴,
Figure FDA00031208362400000122
轴与α轴的夹角为
Figure FDA00031208362400000123
所述旋转电压矢量
Figure FDA00031208362400000124
Figure FDA00031208362400000125
轴上的分量为Uhcosωht,在
Figure FDA00031208362400000126
轴上的分量为Uhsinωht;
步骤S3、将Uhcosωht注入到电机的
Figure FDA00031208362400000127
轴上,求取Uhcosωht产生的零序电压:
Figure FDA00031208362400000128
其中,L2为电机自感的二次谐波的幅值,L0为电机自感的平均值,M0为电机互感的平均值,Ld为电机d轴电感,Lq为电机q轴电感;
将Uhsinωht注入到电机的
Figure FDA00031208362400000129
轴上,求取Uhsinωht产生的零序电压:
Figure FDA00031208362400000130
其中
Figure FDA00031208362400000131
Figure FDA00031208362400000132
轴与电机的实际d轴夹角,
Figure FDA00031208362400000133
将Uhcosωht产生的零序电压和Uhsinωht产生的零序电压求和得到:
Figure FDA0003120836240000021
步骤S4、对求和电压进行变换如下:
Figure FDA0003120836240000022
将变换结果输入至低通滤波器进行滤波可得:
Figure FDA0003120836240000023
其中k为滤波系数;
步骤S5、将步骤S4所述低通滤波结果输入至PI环节,获得电机的估计电角频率,将所述估计电角频率输入至积分环节即可获得电机的估计位置。
2.根据权利要求1所述的一种基于旋转综合矢量的脉振高频注入的永磁电机的位置估计方法,其特征在于,所述步骤S3中Uhcosωht产生的零序电压求取过程如下:
步骤S3.1、记
Figure FDA0003120836240000024
轴与电机的实际d轴夹角为
Figure FDA0003120836240000025
电机的实际d轴与α轴的夹角为θ,则
Figure FDA0003120836240000026
Uhcosωht产生的dq轴电流的变化率为:
Figure FDA0003120836240000027
步骤S3.2,分别计算Uhcosωht产生的三相电流的变化率为:
Figure FDA0003120836240000028
Figure FDA0003120836240000031
步骤S3.3,分别计算Uhcosωht产生的三相电压为:
Figure FDA0003120836240000032
步骤S3.4,计算Uhcosωht产生的零序电压如下:
Figure FDA0003120836240000033
Uhsinωht产生的零序电压求取过程如下:
步骤L3.1、记
Figure FDA0003120836240000034
轴与电机的实际d轴夹角为
Figure FDA0003120836240000035
电机的实际d轴与α轴的夹角为θ,则
Figure FDA0003120836240000036
Uhsinωht产生的dq轴电流的变化率为:
Figure FDA0003120836240000037
步骤L3.2、分别计算Uhsinωht产生的三相电流的变化率如下:
Figure FDA0003120836240000038
Figure FDA0003120836240000041
步骤L3.3、分别计算Uhsinωht产生的三相电压为:
Figure FDA0003120836240000042
步骤L3.4、计算Uhsinωht产生的零序电压为:
Figure FDA0003120836240000043
CN202110675485.4A 2021-06-18 2021-06-18 基于旋转综合矢量的永磁电机高频注入位置估计方法 Active CN113489409B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110675485.4A CN113489409B (zh) 2021-06-18 2021-06-18 基于旋转综合矢量的永磁电机高频注入位置估计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110675485.4A CN113489409B (zh) 2021-06-18 2021-06-18 基于旋转综合矢量的永磁电机高频注入位置估计方法

Publications (2)

Publication Number Publication Date
CN113489409A true CN113489409A (zh) 2021-10-08
CN113489409B CN113489409B (zh) 2023-06-09

Family

ID=77935673

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110675485.4A Active CN113489409B (zh) 2021-06-18 2021-06-18 基于旋转综合矢量的永磁电机高频注入位置估计方法

Country Status (1)

Country Link
CN (1) CN113489409B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114123901A (zh) * 2021-11-24 2022-03-01 南京航空航天大学 一种抑制转矩脉动的逆二倍频同步坐标脉振高频注入的六相永磁电机无位置传感器控制方法
CN114301356A (zh) * 2021-12-14 2022-04-08 南京航空航天大学 基于旋转综合矢量脉振高频电压双套绕组反向注入的无位置传感器控制方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030062870A1 (en) * 2001-09-28 2003-04-03 Semyon Royak Stator and rotor resistance identifier using high frequency injection
US20030169015A1 (en) * 2002-03-05 2003-09-11 Semyon Royak Flux position identifier using high frequency injection with the presence of a rich harmonic spectrum in a responding signal
CN101958686A (zh) * 2009-07-15 2011-01-26 通用汽车环球科技运作公司 无传感器转子角位置估计的方法、系统和设备
CN103956955A (zh) * 2014-04-11 2014-07-30 浙江大学 一种单边可控的共母线开绕组永磁电机系统及其零序电流的抑制方法
CN104579085A (zh) * 2014-12-23 2015-04-29 广东威灵电机制造有限公司 电机转子位置的检测方法和检测装置
CN109546904A (zh) * 2018-11-20 2019-03-29 上海交通大学 双三相永磁同步电机的转子位置检测方法
CN111181456A (zh) * 2020-01-08 2020-05-19 武汉船用电力推进装置研究所(中国船舶重工集团公司第七一二研究所) 一种无位置传感器永磁同步电机振动控制方法
CN112737450A (zh) * 2020-12-24 2021-04-30 上海大学 一种用于spmsm转子位置估计的高频注入补偿方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030062870A1 (en) * 2001-09-28 2003-04-03 Semyon Royak Stator and rotor resistance identifier using high frequency injection
US20030169015A1 (en) * 2002-03-05 2003-09-11 Semyon Royak Flux position identifier using high frequency injection with the presence of a rich harmonic spectrum in a responding signal
CN101958686A (zh) * 2009-07-15 2011-01-26 通用汽车环球科技运作公司 无传感器转子角位置估计的方法、系统和设备
CN103956955A (zh) * 2014-04-11 2014-07-30 浙江大学 一种单边可控的共母线开绕组永磁电机系统及其零序电流的抑制方法
CN104579085A (zh) * 2014-12-23 2015-04-29 广东威灵电机制造有限公司 电机转子位置的检测方法和检测装置
CN109546904A (zh) * 2018-11-20 2019-03-29 上海交通大学 双三相永磁同步电机的转子位置检测方法
CN111181456A (zh) * 2020-01-08 2020-05-19 武汉船用电力推进装置研究所(中国船舶重工集团公司第七一二研究所) 一种无位置传感器永磁同步电机振动控制方法
CN112737450A (zh) * 2020-12-24 2021-04-30 上海大学 一种用于spmsm转子位置估计的高频注入补偿方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
GAOLIN WANG ET AL.: "Low-Frequency Pulse Voltage Injection Scheme-Based Sensorless Control of IPMSM Drives for Audible Noise Reduction", 《IEEE TRANSACTIONS ON INDUSTRIAL ELECTRONICS》 *
魏海峰等: "考虑逆变器输出谐波的永磁同步电机转子位置估计误差分析", 《电工技术学报》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114123901A (zh) * 2021-11-24 2022-03-01 南京航空航天大学 一种抑制转矩脉动的逆二倍频同步坐标脉振高频注入的六相永磁电机无位置传感器控制方法
CN114123901B (zh) * 2021-11-24 2023-12-01 南京航空航天大学 一种抑制转矩脉动的逆二倍频同步坐标脉振高频注入的六相永磁电机无位置传感器控制方法
CN114301356A (zh) * 2021-12-14 2022-04-08 南京航空航天大学 基于旋转综合矢量脉振高频电压双套绕组反向注入的无位置传感器控制方法
CN114301356B (zh) * 2021-12-14 2023-12-01 南京航空航天大学 基于旋转综合矢量脉振高频电压双套绕组反向注入的无位置传感器控制方法

Also Published As

Publication number Publication date
CN113489409B (zh) 2023-06-09

Similar Documents

Publication Publication Date Title
Almarhoon et al. Improved pulsating signal injection using zero-sequence carrier voltage for sensorless control of dual three-phase PMSM
Sun et al. Speed sensorless control for permanent magnet synchronous motors based on finite position set
Jang et al. Sensorless drive of surface-mounted permanent-magnet motor by high-frequency signal injection based on magnetic saliency
Benjak et al. Review of position estimation methods for IPMSM drives without a position sensor part I: Nonadaptive methods
Xu et al. Novel square-wave signal injection method using zero-sequence voltage for sensorless control of PMSM drives
CN110198150B (zh) 一种永磁同步电机多参数在线辨识方法
Bolognani et al. Sensorless control of IPM motors in the low-speed range and at standstill by HF injection and DFT processing
JP4674516B2 (ja) 同期モータの磁極位置推定方法
CN110048655B (zh) 反电势基波提取的永磁同步电机无位置传感器控制系统
CN109981001B (zh) 一种低噪音的永磁同步电机转子初始位置检测方法
Moghadam et al. Sensorless control of PMSMs with tolerance for delays and stator resistance uncertainties
Shuang et al. Improved cross-coupling effect compensation method for sensorless control of IPMSM with high frequency voltage injection
Wang et al. Improved fast method of initial rotor position estimation for interior permanent magnet synchronous motor by symmetric pulse voltage injection
CN113489409A (zh) 一种基于旋转综合矢量的脉振高频注入的永磁电机的位置估计方法
CN106026831B (zh) 一种无位置传感器控制下的电流解耦方法
JP3687590B2 (ja) Pmモータの制御方法、および制御装置
CN115173774B (zh) 一种永磁同步电机无位置传感器控制方法及系统
Wu et al. An optimized asymmetric pulsewidth modulation for sensorless control of permanent magnet synchronous machines
Mamo et al. Novel rotor position extraction based on carrier frequency component method (CFCM) using two reference frames for IPM drives
Qiang et al. Position estimation of AC machines at all frequencies using only space vector PWM based excitation
CN113676106B (zh) 用于六相永磁电机的双套绕组高频注入的无位置方法
CN114362622B (zh) 一种消除观测值波动误差的永磁同步电机无位置传感器控制方法、系统及装置
Lim et al. Sensorless control using a full-order observer based on a novel flux model of high power interior permanent magnet synchronous motor
CN111371363B (zh) 基于信号注入的电机最大转矩电流比控制方法及其装置
Li et al. High frequency response current self-demodulation method for sensorless control of interior permanent magnet synchronous motor

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