CN110808702B - 嵌入式永磁同步电机的控制方法 - Google Patents

嵌入式永磁同步电机的控制方法 Download PDF

Info

Publication number
CN110808702B
CN110808702B CN201911036906.8A CN201911036906A CN110808702B CN 110808702 B CN110808702 B CN 110808702B CN 201911036906 A CN201911036906 A CN 201911036906A CN 110808702 B CN110808702 B CN 110808702B
Authority
CN
China
Prior art keywords
stator
temperature
current
permanent magnet
module
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
CN201911036906.8A
Other languages
English (en)
Other versions
CN110808702A (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.)
CRRC Yongji Electric Co Ltd
Original Assignee
CRRC Yongji Electric Co Ltd
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 CRRC Yongji Electric Co Ltd filed Critical CRRC Yongji Electric Co Ltd
Priority to CN201911036906.8A priority Critical patent/CN110808702B/zh
Publication of CN110808702A publication Critical patent/CN110808702A/zh
Application granted granted Critical
Publication of CN110808702B publication Critical patent/CN110808702B/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
    • 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/141Flux estimation
    • 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/20Estimation of torque
    • 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
    • H02P25/024Synchronous motors controlled by supply frequency
    • H02P25/026Synchronous motors controlled by supply frequency thereby detecting the rotor position
    • 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
    • H02P27/00Arrangements or methods for the control of AC motors characterised by the kind of supply voltage
    • H02P27/04Arrangements or methods for the control of AC motors characterised by the kind of supply voltage using variable-frequency supply voltage, e.g. inverter or converter supply voltage
    • H02P27/06Arrangements or methods for the control of AC motors characterised by the kind of supply voltage using variable-frequency supply voltage, e.g. inverter or converter supply voltage using dc to ac converters or inverters
    • H02P27/08Arrangements or methods for the control of AC motors characterised by the kind of supply voltage using variable-frequency supply voltage, e.g. inverter or converter supply voltage using dc to ac converters or inverters with pulse width modulation
    • 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
    • H02P29/00Arrangements for regulating or controlling electric motors, appropriate for both AC and DC motors
    • H02P29/60Controlling or determining the temperature of the motor or of the drive
    • 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
    • H02P2205/00Indexing scheme relating to controlling arrangements characterised by the control loops
    • H02P2205/01Current loop, i.e. comparison of the motor current with a current reference
    • 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
    • H02P2205/00Indexing scheme relating to controlling arrangements characterised by the control loops
    • H02P2205/05Torque loop, i.e. comparison of the motor torque with a torque reference
    • 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

Abstract

本发明涉及电机的控制方法,具体为嵌入式永磁同步电机的控制方法。解决现有永磁同步电机的控制方法没有实时而准确使用电机参数的问题。该嵌入式永磁同步电机的控制方法在线考虑了温度变化和电机饱和效应对电机参数的影响。通过电机定子温度T、电流幅值I s 、电流矢量角β在线查表得到定子电阻R s (T),定子电感L d (T)L q (T),使用磁链观测模型实时计算出磁链值Ψ f (T r ),提高了电机控制和解耦的准确性;且利用转矩闭环的输出结果重新分配了给定的励磁电流i * d 和转矩电流i * q ,维持了永磁同步电机按照MTPA轨迹运行,降低了电机的发热和损耗。实现了电机的解耦控制、增强了电机控制系统的鲁棒性。

Description

嵌入式永磁同步电机的控制方法
技术领域
本发明涉及电机的控制方法,具体为嵌入式永磁同步电机的控制方法。
背景技术
永磁同步牵引系统在轨道交通领域中运用越来越广泛,特别是嵌入式永磁同步电机需求量很大。对于嵌入式永磁同步电机而言,最重要的性能要求是在广泛的环境下产生准确和高效率的控制性能,而影响永磁电机控制精度的一个最重要因素是电机参数定子电感Ld、Lq、定子电阻Rs和永磁体磁链Ψf的变化。而电机参数的变化主要由电机工作温度和定子电流变化引起的铁芯磁饱和效应引起。传统的矢量控制方法在永磁同步电机控制中被广泛使用,其没有考虑温度对电机参数Rs、Ψf,以及铁芯饱和效应对Ld、Lq的影响,电机输出的转矩精度和电机运行效率将很大程度上被影响。
发明内容
本发明解决现有永磁同步电机控制方法中没有实时而准确使用电机参数,因而电机输出的转矩精度和电机运行效率被影响的问题,提供一种嵌入式永磁同步电机的控制方法。该控制方法实时而准确的确定电机参数,从而降低电机参数不准确对电机控制的影响。同时,在准确的确定电机参数的基础上,提供一种永磁同步电机鲁棒解耦控制方法,实现了电机的解耦控制、增强了电机控制系统的鲁棒性。
本发明是采用如下技术方案实现的:嵌入式永磁同步电机的控制方法,其控制器包括旋转变压器模块、温度传感器模块、定子电阻计算模块、Clark变换模块、Park变换模块、定子电感计算查表模块、永磁体磁链计算模块、转矩计算模块、电流I*计算模块、电流角度计算模块、鲁棒解耦控制器模块、PWM调制模块、三相逆变桥模块;
1)旋转变压器模块
旋转变压器安装在嵌入式永磁同步电机上,通过旋转变压器测量得到永磁同步电机的转子位置θ,转子位置θ经过微分,可以得到永磁同步电机的转速we
2)温度传感器模块
电机定子中埋有温度传感器,由温度传感器得到电机的实时定子温度T;
3)定子电阻计算模块
电机温度变化导致电机定子电阻Rs的变化,测试并绘制电机定子的温度阻值对照表,通过温度传感器获取定子实时温度值T,并通过查询电机定子的温度阻值对照表得到定子电阻Rs(T);
4)Clark变换模块
采集两相定子电流ia、ib,经过Clark变换得到定子电流iα、iβ
5)Park变换模块
定子电流iα、iβ经过Park变换得到d-q静止坐标系下的电流id、iq
6)定子电感计算查表模块
定子电感Ld和Lq与电机定子温度T、定子电流的幅值IS和定子电流的相位βS三个变量相关;
先制作在一定温度区间内的不同温度点下,定子电感Ld、Lq随IS和βS变化的表格;
随后用实时采集到的定子温度T、定子电流的幅值IS和定子电流的相位βS,通过查询表格,得到定子电感Ld(T)、Lq(T);
7)永磁体磁链计算模块
永磁体磁链计算模块的输入为:电流id、iq,定子电阻Rs(T),定子电感Ld(T)和Lq(T),上一计算周期(或称上一拍)的定子电压u’d,u’q,上一计算周期(或称上一拍)计算得到的永磁体磁链Ψ’f(Tr),转速we。永磁体磁链计算模块的输出Ψf(Tr)计算方式如下:
Figure GDA0002939984400000021
其中:Kp_Ψf和Ki_Ψf是控制参数,
Figure GDA0002939984400000022
分别表示估计的电流值,t为采集时刻,Tr代表永磁体的实际工作温度;
调试得到Kp_Ψf和Ki_Ψf的过程如下:给定Ki_Ψf一个较小的参数(例如0.001),先调节参数Kp_Ψf使输出得到的Ψf(Tr)处于等幅震荡的状态,这时再调节参数Ki_Ψf,使输出Ψf(Tr)收敛,这时得到的参数Kp_Ψf和Ki_Ψf是控制参数。
Figure GDA0002939984400000023
Figure GDA0002939984400000024
通过以下公式计算得到:
Figure GDA0002939984400000031
8)转矩计算模块
转矩计算模块的输入为电流id、iq,定子电阻Rs(T),定子电感Ld(T)和Lq(T),转速we,Ψf(Tr);
转矩计算模块的输出——电磁转矩Teb按如下计算公式得到:
Tex1=1.5npiq[(Ld(T)-Lq(T))id]
式中,np是电机的极对数。
Tex2=1.5npiq[(Ld(T0)-Lq(T0))id]
式中,Ld(T0)和Lq(T0)分别是在额定工况下的定子电感的值。
Tex=Tex1-Tex2
Test=1.5npiqf(Tr)+(Ld(T0)-Lq(T0))id]
Teb=Tex+Test
9)电流I*计算模块
目标转矩T* e经过限幅和斜坡处理后,得到给定转矩T* e1
T* e1为电流I*计算模块的输入;
给定转矩T* e1经过计算模块,得到给定电流值I1 *;给定转矩T* e1经过计算模块,得到给定电流值I1 *的计算过程如下:
根据永磁同步电机在额定工况下的电机参数得到运算中用到的标幺值基值teb和ibx,其中ibx是电流的标幺值基值,通过ibx=ψf(T0)/(Lq(T0)-Ld(T0))计算得到,Ld(T0)和Lq(T0)分别是在额定工况下定子电感的值,Ψf(T0)是在额定工况下永磁体磁链的值,以上三个值都是固定不变的;teb是转矩的标幺值基值,可以通过teb=npψf(T0)ibx计算得到;
给定转矩
Figure GDA0002939984400000032
与其标幺值ten可以通过式
Figure GDA0002939984400000033
表示;电流
Figure GDA0002939984400000034
与其标幺值idn可以通过式
Figure GDA0002939984400000035
表示;电流
Figure GDA0002939984400000036
与其标幺值iqn可以通过式
Figure GDA0002939984400000037
表示,电流
Figure GDA0002939984400000038
和电流
Figure GDA0002939984400000041
是计算I1 *的中间变量;
在标幺值的形式下,控制算法的转矩和电流的关系可以表示为:
Figure GDA0002939984400000042
通过将给定转矩指令
Figure GDA0002939984400000043
变为标幺值ten的格式,再通过公式
Figure GDA0002939984400000044
求解得到电流
Figure GDA0002939984400000045
的标幺值idn,最后再通过式
Figure GDA0002939984400000046
可计算得到给定电流
Figure GDA0002939984400000047
电流
Figure GDA0002939984400000048
的标幺值iqn可通过公式ten=iqn(1-idn)计算得到,此时ten和idn是已知量,再通过式
Figure GDA0002939984400000049
可计算得到给定电流
Figure GDA00029399844000000410
给定电流值I1 *通过公式
Figure GDA00029399844000000411
计算得到;
T* e1与Teb的差经过PI调节器,输出为ΔI,ΔI是对电流给定的补偿值,是电机参数不准确时,通过辨识实时转矩,对给定电流进行的补偿量;
最终得到电流I*计算模块的输出是电流给定I*,通过如下公式计算得到:
I*=I1 *+ΔI
10)电流角度计算模块
电流角度计算模块的输入为电机参数Ψf(Tr)、Ld(T)、Lq(T)和电流给定I*
电流角度计算模块的输出为电流角度β1
Figure GDA00029399844000000412
电流幅值I*和β1用来生成给定电流i* d、i* q,计算过程如下:
Figure GDA00029399844000000413
11)鲁棒解耦控制器
鲁棒解耦控制器的输入参数有电流id *、iq *、id、iq、定子电阻Rs(T)、定子电感Ld(T)和Lq(T)、永磁体磁链Ψf(Tr)和转速we,其输出参数为定子电压ud、uq
控制算法中间变量ud1、uq1计算过程如下:
Figure GDA0002939984400000051
式中,βx是控制参数,控制参数βx随着调制策略的不同而变化,控制参数βx与开关频率成一定的比例,其表示如下:
βx=βb×fkx
式中,βb是控制参数基准值,在0.1-1.0之间选择,fkx是逆变器的开关频率;
通过以上的控制算法计算得到ud1、uq1,实现了永磁同步电机的解耦控制;
控制算法中间变量ud2、uq2计算过程如下:
Figure GDA0002939984400000052
式中,δd、δq是鲁棒控制器的关键部分,其计算公式如下:
Figure GDA0002939984400000053
式中,λ为控制参数,控制参数λ采用试凑法选取,可取80;id1、iq1为算法过程中的中间变量,id1、iq1的计算公式如下:
Figure GDA0002939984400000054
其中,δ’d、δ’q是上一计算周期(或称上一拍)计算得到的变量δd、δq
电压ud、uq的计算过程如下:
Figure GDA0002939984400000055
这里为了增强控制系统的稳定性,实现控制解耦,在dq轴上分别加入项RS(T)×kR×id和RS(T)×kR×iqf(Tr)we,KR为阻尼系数,取为0.9;
12)PWM调制模块
PWM调制模块的输入为定子电压ud、uq,直流母线电压udc,转速we和角度θ;PWM调制模块的输出为6路PWM波,驱动三相逆变桥模块工作。
本发明带来的有益效果:
(1)本发明同时在线考虑了温度变化和电机饱和效应对电机参数的影响,提高了永磁同步电机的控制精度。
(2)本发明通过电机定子温度T、电流幅值Is、电流矢量角β在线查表得到准确的电机参数定子电阻Rs(T),定子电感Ld(T)和Lq(T),使用磁链观测模型实时计算出磁链值Ψf(Tr),省去了对转子温度的检测设备,提高了电机控制和解耦的准确性;且利用转矩闭环的输出结果重新分配了给定的励磁电流i* d和转矩电流i* q,维持了永磁同步电机按照较优的轨迹运行,降低了电机的发热和损耗。
(3)通过使用鲁棒解耦控制方法,提高了系统控制的抗干扰性能。
附图说明
图1为本发明所述控制方法的控制框图;
图2为定子电感计算查表模块的流程框图;
图3为电流I*计算模块的控制框图;
图4为鲁棒解耦控制器的控制框图;
图5为分段调制算法示意图。
具体实施方式
嵌入式永磁同步电机的控制方法,其控制器(如图1所示)包括旋转变压器模块1、温度传感器模块2、定子电阻计算模块3、Clark变换模块4、Park变换模块5、定子电感计算查表模块6、永磁体磁链计算模块7、转矩计算模块8、电流I*计算模块9、电流角度计算模块10、鲁棒解耦控制器模块11、PWM调制模块12、三相逆变桥模块13;
1)旋转变压器模块
旋转变压器安装在嵌入式永磁同步电机上,通过旋转变压器可测量得到永磁同步电机的转子位置θ,转子位置θ经过微分,可以得到永磁同步电机的转速we
2)温度传感器模块
电机定子中埋有温度传感器,由温度传感器得到电机的实时定子温度T;
3)定子电阻计算模块
电机温度变化导致电机定子电阻Rs的变化,测试并绘制电机定子的温度阻值对照表,通过温度传感器获取定子实时温度值T,并通过查询电机定子的温度阻值对照表得到定子电阻Rs(T);
4)Clark变换模块
采集两相定子电流ia、ib,经过Clark变换得到定子电流iα、iβ
5)Park变换模块
定子电流iα、iβ经过Park变换得到d-q静止坐标系下的电流id、iq
6)定子电感计算查表模块
定子电流的变化会引起定子铁芯的磁饱和效应,随着d、q轴电流的变化,定子电感Ld、Lq都会发生变化,同时电机的定子温度T也会对定子电感产生影响。为了得到较为准确的定子电感参数Ld和Lq,采用查表法来得到定子电感Ld和Lq
定子电感Ld和Lq与电机定子温度T、定子电流的幅值IS和定子电流的相位βS三个变量相关(如图2所示);
其中:
Figure GDA0002939984400000071
考虑到iq为0时,程序计算时会出现问题,因此给分母加一个特别小的数字kα,kα可以等于0.0000001;
进一步地,定子电流的幅值IS和定子电流的相位βS分别经过低通滤波,得到滤波后的值ISLPF和βSLPF
先制作在一定温度区间内的不同温度点下,定子电感Ld、Lq随IS和βS(或ISLPF和βSLPF)变化的表格(表格数量与温度区间内所设温度点的数量一致,与每个温度点对应的表格反映定子电感Ld、Lq随IS和βS(或ISLPF和βSLPF)变化的内容);具体实施时,温度区间为[-30℃,160℃]区间,并以每个十的整数倍的温度值作为温度点,这样,在[-30℃,160℃]区间,有二十个温度点:-30℃、-20℃、-10℃、0℃、10℃、20℃、……、150℃、160℃;具体的制作表格的过程为:测试的永磁同步电机和对拖电机同轴或通过齿轮箱连接,将对拖电机运行在额定转速以下,在转速稳定且定子温度在某一温度点(以-20℃为例)的情况下,给定永磁同步电机不同的电流指令IS(IS>0)和βS(90°<βS<180°)(或ISLPF(ISLPF>0)和βSLPF(90°<βSLPF<180°))(电流i* d、i* q),用高频注入法测得与不同的IS和βS(或ISLPF和βSLPF)对应的电机定子电感Ld、Lq;改变定子温度,从而获得多组不同温度点下的定子电感Ld、Lq随IS和βS(或ISLPF和βSLPF)变化的表格。
随后用实时采集到的定子温度T、定子电流的幅值IS和定子电流的相位βS(或ISLPF和βSLPF),通过查询表格,得到实时的定子电感Ld(T)、Lq(T);查表的过程如下:实时采集的定子温度T不等于任何温度点的温度值时,选择离定子温度T最近的两个温度点Tx、Tx+10的定子电感Ld、Lq随IS和βS(或ISLPF和βSLPF)变化的表格作为查表使用的第一表格和第二表格(如,定子温度T为23℃,选择温度点Tx=20℃和Tx+10=30℃的表格作为第一和第二表格);实时采集的定子温度T等于某温度点的温度值时,选择与定子温度T相等的温度点Tx的定子电感Ld、Lq随IS和βS(或ISLPF和βSLPF)变化的表格,和该温度点增温侧相邻的温度点Tx+10的定子电感Ld、Lq随IS和βS(或ISLPF和βSLPF)变化的表格,作为查表使用的第一表格和第二表格;根据采集的定子电流的幅值IS和定子电流的相位βS(或ISLPF和βSLPF),在第一表格里查到Ld1(T)和Lq1(T),在第二表格里查到Ld2(T)和Lq2(T);按如下公式得到Ld(T)和Lq(T):
Figure GDA0002939984400000081
Ld(T)和Lq(T)为当前工况下(当前温度T,当前IS和βS(或ISLPF和βSLPF))的电机定子电感值。
7)永磁体磁链计算模块
永磁同步电机转子中永磁材料受温度变化影响大,电机永磁体磁链与温度的变化关系可以表示为
Figure GDA0002939984400000082
其中:Ψf(T0)是额定工况下的永磁体磁链,Ψf(Tr)是电机实际运行温度下的永磁体磁链,T0可取为20℃,Tr是永磁体实际工作温度,α是剩磁的温度系数。
由于电机的永磁体安装在电机的转子上,无法得到其实际工作温度,因此采用磁链观测模型实时得到磁链值Ψf(Tr),永磁体磁链计算模块的输入为:电流id、iq,定子电阻Rs(T),定子电感Ld(T)和Lq(T),上一计算周期(或称上一拍)的定子电压u’d,u’q,上一计算周期(或称上一拍)计算得到的永磁体磁链Ψf’(Tr),转速we,永磁体磁链计算模块的输出Ψf(Tr)计算方式如下:
Figure GDA0002939984400000091
其中:Kp_Ψf和Ki_Ψf是控制参数,
Figure GDA0002939984400000092
分别表示估计的电流值,t为采集时刻,Tr代表永磁体的实际工作温度;
调试得到Kp_Ψf和Ki_Ψf的过程如下:给定Ki_Ψf一个较小的参数(例如0.001),先调节参数Kp_Ψf使输出得到的Ψf(Tr)处于等幅震荡的状态,这时再调节参数Ki_Ψf,使输出Ψf(Tr)收敛,这时得到的参数Kp_Ψf和Ki_Ψf是控制参数;
Figure GDA0002939984400000093
Figure GDA0002939984400000094
通过以下公式计算得到:
Figure GDA0002939984400000095
8)转矩计算模块
转矩计算模块的输入为电流id、iq,定子电阻Rs(T),定子电感Ld(T)和Lq(T),转速we,Ψf(Tr);
转矩计算模块的输出——电磁转矩Teb按如下计算公式得到:
Tex1=1.5npiq[(Ld(T)-Lq(T))id]
式中,np是电机的极对数;
Tex2=1.5npiq[(Ld(T0)-Lq(T0))id]
式中,Ld(T0)和Lq(T0)分别是在额定工况下的定子电感的值;
Tex=Tex1-Tex2
Test=1.5npiqf(Tr)+(Ld(T0)-Lq(T0))id]
Teb=Tex+Test
9)电流I*计算模块
目标转矩T* e经过限幅和斜坡处理后,得到给定转矩T* e1
T* e1为电流I*计算模块的输入;
给定转矩T* e1经过计算模块,得到给定电流值I1 *(如图3所示);给定转矩T* e1经过计算模块,得到给定电流值I1 *的计算过程如下:
根据永磁同步电机在额定工况下的电机参数得到运算中用到的标幺值基值teb和ibx,其中ibx是电流的标幺值基值,通过ibx=ψf(T0)/(Lq(T0)-Ld(T0))计算得到,Ld(T0)和Lq(T0)分别是在额定工况下定子电感的值,Ψf(T0)是在额定工况下永磁体磁链的值,以上三个值都是固定不变的;teb是转矩的标幺值基值,可以通过teb=npψf(T0)ibx计算得到;
转矩指令
Figure GDA0002939984400000101
与其标幺值ten可以通过式
Figure GDA0002939984400000102
表示;电流
Figure GDA0002939984400000103
与其标幺值idn可以通过式
Figure GDA0002939984400000104
表示;电流
Figure GDA0002939984400000105
与其标幺值iqn可以通过式
Figure GDA0002939984400000106
表示,电流
Figure GDA0002939984400000107
和电流
Figure GDA0002939984400000108
是计算I1 *的中间变量;
在标幺值的形式下,控制算法的转矩和电流的关系表示为:
Figure GDA0002939984400000109
通过将给定转矩指令
Figure GDA00029399844000001010
变为标幺值ten的格式,再通过公式
Figure GDA00029399844000001011
求解得到电流
Figure GDA00029399844000001012
的标幺值idn,最后再通过式
Figure GDA00029399844000001013
可计算得到给定电流
Figure GDA00029399844000001014
电流
Figure GDA00029399844000001015
的标幺值iqn可通过公式ten=iqn(1-idn)计算得到,此时ten和idn是已知量,再通过式
Figure GDA00029399844000001016
可计算得到给定电流
Figure GDA00029399844000001017
给定电流值I1 *通过公式
Figure GDA00029399844000001018
计算得到;
T* e1与Teb的差经过PI调节器,输出为ΔI,ΔI是对电流给定的补偿值,是电机参数不准确时,通过辨识实时转矩,对给定电流进行的补偿量;
最终得到电流I*计算模块的输出是电流给定I*,通过如下公式计算得到:
I*=I1 *+ΔI
10)电流角度计算模块
电流角度计算模块的输入为电机参数Ψf(Tr)、Ld(T)、Lq(T)和电流给定I*
电流角度计算模块的输出为电流角度β1
Figure GDA0002939984400000111
电流幅值I*和β1用来生成给定电流i* d、i* q,计算过程如下:
Figure GDA0002939984400000112
11)鲁棒解耦控制器
鲁棒解耦控制器一方面实现永磁同步电机的解耦控制,另一方面提高电机控制的抗干扰功能,其输入参数有电流id *、iq *、id、iq、定子电阻Rs(T)、定子电感Ld(T)和Lq(T)、永磁体磁链Ψf(Tr)和转速we,其输出参数有定子电压ud、uq
鲁棒解耦控制器如图4所示。
控制算法中间变量ud1、uq1计算过程如下:
Figure GDA0002939984400000113
式中,βx是控制参数,控制参数βx随着调制策略的不同而变化,控制参数βx与开关频率成一定的比例,其表示如下:
βx=βb×fkx
式中,βb是控制参数基准值,在0.1-1.0之间选择(如,0.1、0.2、0.3、0.4、0.5、0.7、0.9、1.0;其中优选0.5),fkx是逆变器的开关频率;
通过以上的控制算法计算得到ud1、uq1,实现了永磁同步电机的解耦控制;
控制算法中间变量ud2、uq2计算过程如下:
Figure GDA0002939984400000114
式中,δd、δq是鲁棒控制器的关键部分,其计算公式如下:
Figure GDA0002939984400000121
式中,λ为控制参数,控制参数λ采用试凑法选取,可取80。id1、iq1为算法过程中的中间变量,id1、iq1的计算公式如下:
Figure GDA0002939984400000122
其中,δ’d、δ’q是上一计算周期(或称上一拍)计算得到的变量δd、δq
电压ud、uq的计算过程如下:
Figure GDA0002939984400000123
这里为了增强控制系统的稳定性,实现控制解耦,在dq轴上分别加入项RS(T)×kR×id和RS(T)×kR×iqf(Tr)we,KR为阻尼系数,可以取为0.9;
12)PWM调制模块
PWM调制模块的输入为定子电压ud、uq,直流母线电压udc,转速we和角度θ。PWM调制模块的输出为6路PWM波,驱动三相逆变桥模块工作。
因散热等条件的制约,大功率永磁同步电机的调制算法受到开关频率的限制,一般采用多种调制方式相结合的分段调制策略。分段调制策略的示意图如图5所示。
分段调制分为异步调制和同步调制,同步调制受到开关频率的限制可分为多段,最终进入方波调制。方波调制下,电压利用率高、谐波小。同步分段调制算法中,可采用的调制算法有SPWM调制、特定次谐波消除PWM(SHEPWM)调制等,各种调制算法有其优缺点和适用范围。
本发明中提出的控制方法,可以在DSP芯片中实现,永磁同步电机控制由两级中断实现,一级中断中运行电机控制算法(即步骤1-11),设计控制周期为250us;另一级中断中运行电机PWM调制算法,其中断周期与当前的调制策略与电机运行频率相关。

Claims (10)

1.嵌入式永磁同步电机的控制方法,其特征在于,其控制器包括旋转变压器模块(1)、温度传感器模块(2)、定子电阻计算模块(3)、Clark变换模块(4)、Park变换模块(5)、定子电感计算查表模块(6)、永磁体磁链计算模块(7)、转矩计算模块(8)、电流I*计算模块(9)、电流角度计算模块(10)、鲁棒解耦控制器模块(11)、PWM调制模块(12)、三相逆变桥模块(13);
1)旋转变压器模块
旋转变压器安装在嵌入式永磁同步电机上,通过旋转变压器测量得到永磁同步电机的转子位置θ,转子位置θ经过微分,得到永磁同步电机的转速we
2)温度传感器模块
电机定子中埋有温度传感器,由温度传感器得到电机的实时定子温度T;
3)定子电阻计算模块
电机温度变化导致电机定子电阻Rs的变化,测试并绘制电机定子的温度阻值对照表,通过温度传感器获取定子实时温度值T,并通过查询电机定子的温度阻值对照表得到定子电阻Rs(T);
4)Clark变换模块
采集两相定子电流ia、ib,经过Clark变换得到定子电流iα、iβ
5)Park变换模块
定子电流iα、iβ经过Park变换得到d-q静止坐标系下的电流id、iq
6)定子电感计算查表模块
定子电感Ld和Lq与电机定子温度T、定子电流的幅值IS和定子电流的相位βS三个变量相关;
先制作在一定温度区间内的不同温度点下,定子电感Ld、Lq随IS和βS变化的表格;
随后用实时采集到的定子温度T、定子电流的幅值IS和定子电流的相位βS,通过查询表格,得到定子电感Ld(T)、Lq(T);
7)永磁体磁链计算模块
永磁体磁链计算模块的输入为:电流id、iq,定子电阻Rs(T),定子电感Ld(T)和Lq(T),上一计算周期的定子电压u’d,u’q,上一计算周期计算得到的永磁体磁链Ψ’f(Tr),转速we,永磁体磁链计算模块的输出Ψf(Tr)计算方式如下:
Figure FDA0002939984390000021
其中:Kp_Ψf和Ki_Ψf是控制参数,
Figure FDA0002939984390000022
分别表示估计的电流值,t为采集时刻,Tr代表永磁体的实际工作温度;
调试得到Kp_Ψf和Ki_Ψf的过程如下:给定Ki_Ψf一个较小的参数,先调节参数Kp_Ψf使输出得到的Ψf(Tr)处于等幅震荡的状态,这时再调节参数Ki_Ψf,使输出Ψf(Tr)收敛,这时得到的参数Kp_Ψf和Ki_Ψf是控制参数;
Figure FDA0002939984390000023
Figure FDA0002939984390000024
通过以下公式计算得到:
Figure FDA0002939984390000025
8)转矩计算模块
转矩计算模块的输入为电流id、iq,定子电阻Rs(T),定子电感Ld(T)和Lq(T),转速we,Ψf(Tr);
转矩计算模块的输出——电磁转矩Teb按如下计算公式得到:
Tex1=1.5npiq[(Ld(T)-Lq(T))id]
式中,np是电机的极对数,
Tex2=1.5npiq[(Ld(T0)-Lq(T0))id]
式中,Ld(T0)和Lq(T0)分别是在额定工况下的定子电感的值;
Tex=Tex1-Tex2
Test=1.5npiqf(Tr)+(Ld(T0)-Lq(T0))id]
Teb=Tex+Test
9)电流I*计算模块
目标转矩T* e经过限幅和斜坡处理后,得到给定转矩T* e1
T* e1为电流I*计算模块的输入;
给定转矩T* e1经过计算模块,得到给定电流值I1 *
T* e1与Teb的差经过PI调节器,输出为ΔI;
最终得到电流I*计算模块的输出是电流给定I*,通过如下公式计算得到:
I*=I1 *+ΔI
10)电流角度计算模块
电流角度计算模块的输入为电机参数Ψf(Tr)、Ld(T)、Lq(T)和电流给定I*
电流角度计算模块的输出为电流角度β1
Figure FDA0002939984390000031
电流幅值I*和β1用来生成给定电流i* d、i* q,计算过程如下:
Figure FDA0002939984390000032
11)鲁棒解耦控制器
鲁棒解耦控制器的输入参数有电流id *、iq *、id、iq、定子电阻Rs(T)、定子电感Ld(T)和Lq(T)、永磁体磁链Ψf(Tr)和转速we,其输出参数为定子电压ud、uq
控制算法中间变量ud1、uq1计算过程如下:
Figure FDA0002939984390000033
式中,βx是控制参数,其表示如下:
βx=βb×fkx
式中,βb是控制参数基准值,在0.1-1.0之间选择,fkx是逆变器的开关频率;
控制算法中间变量ud2、uq2计算过程如下:
Figure FDA0002939984390000034
式中,δd、δq的计算公式如下:
Figure FDA0002939984390000041
式中,λ为控制参数,控制参数λ采用试凑法选取;
id1、iq1为算法过程中的中间变量,id1、iq1的计算公式如下:
Figure FDA0002939984390000042
其中,δ’d、δ’q是上一计算周期计算得到的变量δd、δq
电压ud、uq的计算过程如下:
Figure FDA0002939984390000043
KR为阻尼系数;
12)PWM调制模块
PWM调制模块的输入为定子电压ud、uq,直流母线电压udc,转速we和角度θ;PWM调制模块的输出为6路PWM波,驱动三相逆变桥模块工作。
2.根据权利要求1所述的嵌入式永磁同步电机的控制方法,其特征在于,在6)定子电感计算查表模块中:
Figure FDA0002939984390000044
考虑到iq为0时,程序计算时会出现问题,因此给分母加一个特别小的数字kα
3.根据权利要求2所述的嵌入式永磁同步电机的控制方法,其特征在于,在6)定子电感计算查表模块中:定子电流的幅值IS和定子电流的相位βS分别经过低通滤波,得到滤波后的值ISLPF和βSLPF
4.根据权利要求2所述的嵌入式永磁同步电机的控制方法,其特征在于,在6)定子电感计算查表模块中:温度区间为[-30℃,160℃]区间,并以每个十的整数倍的温度值作为温度点,这样,在[-30℃,160℃]区间,有二十个温度点:-30℃、-20℃、-10℃、0℃、10℃、20℃、……、150℃、160℃。
5.根据权利要求3所述的嵌入式永磁同步电机的控制方法,其特征在于,在6)定子电感计算查表模块中:温度区间为[-30℃,160℃]区间,并以每个十的整数倍的温度值作为温度点,这样,在[-30℃,160℃]区间,有二十个温度点:-30℃、-20℃、-10℃、0℃、10℃、20℃、……、150℃、160℃。
6.根据权利要求4所述的嵌入式永磁同步电机的控制方法,其特征在于,制作表格的过程为:测试的永磁同步电机和对拖电机同轴或通过齿轮箱连接,将对拖电机运行在额定转速以下,在转速稳定且定子温度在某一温度点的情况下,给定永磁同步电机不同的IS,IS>0和βS,90°<βS<180°,用高频注入法测得与不同的IS和βS对应的电机定子电感Ld、Lq;改变定子温度,从而获得多组不同温度点下的定子电感Ld、Lq随IS和βS变化的表格。
7.根据权利要求5所述的嵌入式永磁同步电机的控制方法,其特征在于,制作表格的过程为:测试的永磁同步电机和对拖电机同轴或通过齿轮箱连接,将对拖电机运行在额定转速以下,在转速稳定且定子温度在某一温度点的情况下,给定永磁同步电机不同的ISLPF,ISLPF>0和βSLPF,90°<βSLPF<180°,用高频注入法测得与不同的ISLPF和βSLPF对应的电机定子电感Ld、Lq;改变定子温度,从而获得多组不同温度点下的定子电感Ld、Lq随ISLPF和βSLPF变化的表格。
8.根据权利要求6所述的嵌入式永磁同步电机的控制方法,其特征在于,在6)定子电感计算查表模块中:查表的过程如下:实时采集的定子温度T不等于任何温度点的温度值时,选择离定子温度T最近的两个温度点Tx、Tx+10的定子电感Ld、Lq随IS和βS变化的表格作为查表使用的第一表格和第二表格;实时采集的定子温度T等于某温度点的温度值时,选择与定子温度T相等的温度点Tx的定子电感Ld、Lq随IS和βS变化的表格,和该温度点增温侧相邻的温度点Tx+10的定子电感Ld、Lq随IS和βS变化的表格,作为查表使用的第一表格和第二表格;根据采集的定子电流的幅值IS和定子电流的相位βS,在第一表格里查到Ld1(T)和Lq1(T),在第二表格里查到Ld2(T)和Lq2(T);按如下公式得到Ld(T)和Lq(T):
Figure FDA0002939984390000051
9.根据权利要求7所述的嵌入式永磁同步电机的控制方法,其特征在于,在6)定子电感计算查表模块中:查表的过程如下:实时采集的定子温度T不等于任何温度点的温度值时,选择离定子温度T最近的两个温度点Tx、Tx+10的定子电感Ld、Lq随ISLPF和βSLPF变化的表格作为查表使用的第一表格和第二表格;实时采集的定子温度T等于某温度点的温度值时,选择与定子温度T相等的温度点Tx的定子电感Ld、Lq随ISLPF和βSLPF变化的表格,和该温度点增温侧相邻的温度点Tx+10的定子电感Ld、Lq随ISLPF和βSLPF变化的表格,作为查表使用的第一表格和第二表格;根据采集的ISLPF和βSLPF,在第一表格里查到Ld1(T)和Lq1(T),在第二表格里查到Ld2(T)和Lq2(T);按如下公式得到Ld(T)和Lq(T):
Figure FDA0002939984390000061
10.根据权利要求1-9任一个所述的嵌入式永磁同步电机的控制方法,其特征在于,9)电流I*计算模块中,给定转矩T* e1经过计算模块,得到给定电流值I1 *的计算过程如下:
根据永磁同步电机在额定工况下的电机参数得到运算中用到的标幺值基值teb和ibx,其中ibx是电流的标幺值基值,通过ibx=ψf(T0)/(Lq(T0)-Ld(T0))计算得到,Ld(T0)和Lq(T0)分别是在额定工况下定子电感的值,Ψf(T0)是在额定工况下永磁体磁链的值,以上三个值都是固定不变的;teb是转矩的标幺值基值,可以通过teb=npψf(T0)ibx计算得到;
转矩指令
Figure FDA0002939984390000062
与其标幺值ten可以通过式
Figure FDA0002939984390000063
表示;电流
Figure FDA0002939984390000064
与其标幺值idn可以通过式
Figure FDA0002939984390000065
表示;电流
Figure FDA0002939984390000066
与其标幺值iqn可以通过式
Figure FDA0002939984390000067
表示,电流
Figure FDA0002939984390000068
和电流
Figure FDA0002939984390000069
是计算I1 *的中间变量;
在标幺值的形式下,控制算法的转矩和电流的关系可以表示为:
Figure FDA00029399843900000610
通过将给定转矩
Figure FDA00029399843900000611
变为标幺值ten的格式,再通过公式
Figure FDA00029399843900000612
求解得到电流
Figure FDA00029399843900000613
的标幺值idn,最后再通过式
Figure FDA00029399843900000614
可计算得到给定电流
Figure FDA00029399843900000615
电流
Figure FDA00029399843900000616
的标幺值iqn可通过公式ten=iqn(1-idn)计算得到,此时ten和idn是已知量,再通过式
Figure FDA0002939984390000071
可计算得到给定电流
Figure FDA0002939984390000072
给定电流值I1 *通过公式
Figure FDA0002939984390000073
计算得到。
CN201911036906.8A 2019-10-29 2019-10-29 嵌入式永磁同步电机的控制方法 Active CN110808702B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911036906.8A CN110808702B (zh) 2019-10-29 2019-10-29 嵌入式永磁同步电机的控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911036906.8A CN110808702B (zh) 2019-10-29 2019-10-29 嵌入式永磁同步电机的控制方法

Publications (2)

Publication Number Publication Date
CN110808702A CN110808702A (zh) 2020-02-18
CN110808702B true CN110808702B (zh) 2021-03-26

Family

ID=69489422

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911036906.8A Active CN110808702B (zh) 2019-10-29 2019-10-29 嵌入式永磁同步电机的控制方法

Country Status (1)

Country Link
CN (1) CN110808702B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112737441B (zh) * 2020-12-25 2022-11-22 中车永济电机有限公司 一种永磁辅助同步磁阻电机的控制方法
CN113315432B (zh) * 2021-03-31 2023-09-12 金风科技股份有限公司 风力发电机的参数辨识方法以及参数辨识装置
CN113489398B (zh) * 2021-06-03 2022-04-01 西安理工大学 一种内置式永磁同步电机无位置传感器参数误差补偿策略
CN114499331A (zh) * 2021-12-28 2022-05-13 中车永济电机有限公司 一种永磁同步电机控制方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101295953A (zh) * 2007-04-25 2008-10-29 株式会社日立制作所 永磁电动机的弱磁控制装置及应用它的电动动力转向器
CN106655941A (zh) * 2017-01-24 2017-05-10 广州汽车集团股份有限公司 一种内嵌式永磁同步电机的参数估计方法及装置
CN107645259A (zh) * 2017-07-26 2018-01-30 中国第汽车股份有限公司 一种电动汽车用驱动电机转矩闭环控制方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2019009970A (ja) * 2017-06-28 2019-01-17 ルネサスエレクトロニクス株式会社 モータ駆動装置およびモータシステム
JP6647262B2 (ja) * 2017-10-04 2020-02-14 キヤノン株式会社 モータ制御装置、シート搬送装置、原稿読取装置及び画像形成装置

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101295953A (zh) * 2007-04-25 2008-10-29 株式会社日立制作所 永磁电动机的弱磁控制装置及应用它的电动动力转向器
CN106655941A (zh) * 2017-01-24 2017-05-10 广州汽车集团股份有限公司 一种内嵌式永磁同步电机的参数估计方法及装置
CN107645259A (zh) * 2017-07-26 2018-01-30 中国第汽车股份有限公司 一种电动汽车用驱动电机转矩闭环控制方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
A Torque Compensation Control Scheme of PMSM Considering Wide Variation of Permanent Magnet Temperature;Su-Yeon Cho 等;《IEEE TRANSACTIONS ON MAGNETICS》;20190228;第55卷(第2期);全文 *

Also Published As

Publication number Publication date
CN110808702A (zh) 2020-02-18

Similar Documents

Publication Publication Date Title
CN110808702B (zh) 嵌入式永磁同步电机的控制方法
CN107645259B (zh) 一种电动汽车用驱动电机转矩闭环控制方法
CN109194229B (zh) 一种基于转矩闭环的永磁同步电机mtpa控制系统及方法
CN110784141B (zh) 永磁同步电机的控制方法
Cheok et al. A new torque and flux control method for switched reluctance motor drives
CN103931096B (zh) 用温度补偿控制电动机的方法和系统
CN103187919B (zh) 一种永磁同步电机弱磁调速的系统和方法
CN110868120B (zh) 一种内置式永磁同步电机控制方法
CN104579083A (zh) 永磁同步电机矢量控制方法及系统
CN108521243B (zh) 一种基于空间矢量调制的高速永磁同步电机直接功率控制方法
WO2022134772A1 (zh) 一种永磁辅助同步磁阻电机的控制方法
CN103178769A (zh) 永磁同步电机静止状态下参数离线辨识方法
CN112740537B (zh) 永磁同步电机的mtpa控制方法、装置、系统及设备
CN107294447B (zh) 永磁同步电机的自适应最大转矩电流比控制装置
JP5595835B2 (ja) 電動機の駆動装置
CN111245328B (zh) 查表法结合调节器的永磁同步电机控制方法
CN109782173B (zh) 异步电机励磁互感曲线测量系统及其测量方法
Korlinchak et al. Sensorless field orientation of an induction motor drive using a time-varying observer
CN104767457B (zh) 直流变频压缩机运行过程中参数自适应的方法
CN110112974A (zh) 电机控制方法、控制器、存储介质及电机驱动系统
CN110649851B (zh) 异步电机多参数解耦在线辨识方法
Zhu et al. Phase current reconstruction error suppression method for single DC-link shunt PMSM drives at low-speed region
CN107276440A (zh) 一种逆变器的非线性补偿装置、系统及控制方法
CN112821834B (zh) 一种永磁同步电机的在线参数辨识方法和装置
CN111082726B (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