CN114400942A - 一种五相双励磁同步电机最大转矩铜耗比协调控制方法 - Google Patents

一种五相双励磁同步电机最大转矩铜耗比协调控制方法 Download PDF

Info

Publication number
CN114400942A
CN114400942A CN202111668934.9A CN202111668934A CN114400942A CN 114400942 A CN114400942 A CN 114400942A CN 202111668934 A CN202111668934 A CN 202111668934A CN 114400942 A CN114400942 A CN 114400942A
Authority
CN
China
Prior art keywords
current
axis
excitation
motor
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.)
Withdrawn
Application number
CN202111668934.9A
Other languages
English (en)
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.)
Xian University of Technology
Original Assignee
Xian University of Technology
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 Xian University of Technology filed Critical Xian University of Technology
Priority to CN202111668934.9A priority Critical patent/CN114400942A/zh
Publication of CN114400942A publication Critical patent/CN114400942A/zh
Withdrawn legal-status Critical Current

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
    • 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
    • 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

Landscapes

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

Abstract

本发明公开了一种五相双励磁同步电机最大转矩铜耗比协调控制方法,通过电流矢量控制方法,结合旋转坐标系下的数学模型计算d轴电流、q轴电流和励磁电流的参考值;基于电压预测模型、参考电流及实际电流反馈值,最终实现五相混合励磁同步电机铜耗最小分区协调控制方法。本发明提出的控制方法通过定子电枢电流交直轴解耦与励磁电流的协调控制,能够实现电机在全速度范围内保持铜耗最小的高效运行策略,包括三次谐波电流的注入对损耗的影响有更全面的研究;能进行全速度区范围的的平稳切换,综合利用电枢电流和励磁电流对转矩与磁场的调节作用,提升了电机的综合性能。

Description

一种五相双励磁同步电机最大转矩铜耗比协调控制方法
技术领域
本发明属于永磁励磁同步电机控制技术领域,具体涉及一种五相双励磁同步电机最大转矩铜耗比协调控制方法。
背景技术
混合励磁同步电机是由20世纪八十年代末美国学者提出的一种新型电机结构,它结合了电励磁同步电机调磁方便和永磁同步电机效率高等优点,又克服了永磁同步电机气隙磁场难以调节的缺点,具有很好的应用价值。多相电机相比于传统的三相电机,多相结构的输出转矩脉动更小,容错能力更强,稳定持续运行等特点,更加适合电动汽车的应用场合。
五相双励磁同步电机内部有永磁和电励磁两种励磁源,其气隙磁场一般以永磁源产生的磁场为主,电励磁产生的磁场为辅,该类电机结合了混合励磁同步电机与多相电机的优点。近些年,我国电动汽车迎来了新一轮的高速发展,由于电动汽车工作环境较为苛刻,对整个系统的安全可靠性也有着较高的要求。五相双励磁同步电机内部有永磁和电励磁两种励磁源,其气隙磁场一般以永磁源产生的磁场为主,电励磁产生的磁场为辅,该类电机结合了混合励磁同步电机与多相电机的优点。因此,五相双励磁同步电机能更好的满足电动汽车驱动系统起动转矩大、过载能力强、可靠性高、调速范围宽、效率高等要求。为实现五相双励磁同步电机的可靠运行,提出一种五相双励磁同步电机最大转矩铜耗比协调控制方法,提高双励磁同步电机驱动系统的稳定性、鲁棒性、可靠性、快速响应能力,推进该类电机在电动汽车领域的应用。
发明内容
本发明的目的是提供一种五相双励磁同步电机最大转矩铜耗比协调控制方法,实现了五相双励磁同步电机在全区范围内的稳定、可靠、快速、高效运行。
本发明所采用的技术方案是,一种五相双励磁同步电机最大转矩铜耗比协调控制方法,具体按照以下步骤实施:
步骤1、从五相双励磁同步电机主电路采集五相电流iA、iB、iC、iD、iE和励磁电流if,逆变器母线电压Udc和励磁电压Uf,,对电机位置进行初始位置检测,将编码器采集到的信号送入控制器进行处理后得到电机转子位置角θ和转速n;
步骤2:将步骤1采集得到的相电流iA、iB、iC、iD、iE经A/D转换后,然后进行坐标变换分别得到两相旋转坐标系下关于基波平面与三次谐波平面对应的定子d轴电流id1与id3、q轴电流iq1与iq3
步骤3:判断电机运行区域,当电机运行于低速区时,利用拉格朗日乘数法寻找铜耗的最小值,计算得到低速区关于基波平面与三次谐波平面对应的定子d轴电流、q轴电流和励磁电流参考值;
当电机运行于中速区时,采用励磁电流为零的简单控制,根据拉格朗日方程求得铜耗的最小值,计算得到中速区关于基波平面与三次谐波平面对应的定子d轴电流、q轴电流和励磁电流参考值;
当电机运行于高速区时,建立拉格朗日方程,求得铜耗的最小值,计算得到高速区关于基波平面与三次谐波平面对应的定子d轴电流、q轴电流和励磁电流参考值;
步骤4:通过步骤3,将最大转矩铜耗比协调控制策略产生的d轴参考电流,q轴参考电流分别送入PI调节模块后进行比例积分运算,分别获得在基波平面与三次谐波平面的d轴电压与q轴电压,然后将d轴电压与q轴电压进行旋转正交-静止两相变换后,获得静止两相坐标系下α轴电压和β轴电压,再通过空间矢量脉冲宽度调制模块后输出10路脉冲宽度调制信号,从而驱动主功率变换器工作;同时步骤1得到的励磁电流if,经跟随、滤波、偏置与A/D转换后和励磁电流参考值一起送入直流励磁脉宽调制模块,然后输出4路脉冲宽度信号来驱动励磁功率变换器。
本发明的特点还在于,
步骤1中,
五相双励磁同步电机在两相旋转坐标系下的数学模型为:
磁链方程:
Figure BDA0003448930900000031
电压方程:
Figure BDA0003448930900000041
电磁转矩方程:
Figure BDA0003448930900000042
其中,脚标“1”和“3”分别代表基波平面和三次谐波平面成分,ψd1、ψd3为定子d轴磁链,ψq1、ψq3为定子q轴磁链,ψm1、ψm3为永磁磁链幅值,ψmf1、ψmf3为励磁绕组磁链幅值;Ld1、Ld3为d轴电感,Lq1、Lq3为q轴电感,Lf1、Lf3为励磁绕组电感,Mf1、Mf3为电枢与励磁绕组之间的互感;id1、id3为d轴电流,iq1、iq3为q轴电流,if1、if3为励磁绕组电流,ωe为电角速度;ud1、ud3为d轴电压,uq1、uq3为q轴电压,uf1、uf3为励磁绕组电压,Rs为电枢绕组电阻,Rf为励磁绕组电阻;Te1、Te3为电磁转矩,p为电机极对数。
步骤3中,
计算低速区d轴参考电流、q轴参考电流和励磁电流参考值的方法如下:
将电机的铜耗pcu表示为d轴电流、q轴电流和励磁电流的函数,如下所示:
Figure BDA0003448930900000051
电机运行转速小于额定转速时,将总转矩方程(3)改写为
Figure BDA0003448930900000052
同时,将铜耗方程(4)改写为:
Figure BDA0003448930900000053
电机运行在低速区时,基于式(5)的电磁转矩方程约束条件,建立如下所示拉格朗日函数:
Figure BDA0003448930900000054
将式(7)分别对iq1、iq3、if1、if3、λ1和λ3求导,可得:
Figure BDA0003448930900000055
令上式
Figure BDA0003448930900000061
Figure BDA0003448930900000062
得到如下所示的低速区参考电流:
Figure BDA0003448930900000063
其中,
Figure BDA0003448930900000064
k5=k1Mf1 2,k6=2k1ψm1Mf1+k2Mf1 2
Figure BDA0003448930900000065
k9=k5Mf3 2,k10=2k3ψm3Mf3+k4Mf3 2
Figure BDA0003448930900000066
T01=(2Te1/5p)2,T03=(2Te3/15p)2
式(9)中的最后两个方程为一元四次方程,通过求解方程中的ifref1与ifref3,带入函数表达式可解除d、q轴参考电流。
步骤3中,计算中速区d轴参考电流、q轴参考电流和励磁电流参考值的方法如下:
利用转矩电流最大比控制,中速区铜耗方程为
Figure BDA0003448930900000067
将总转矩方程(3)改写为
Figure BDA0003448930900000068
电机运行在中速区时,基于电磁转矩方程约束条件,建立如下所示拉格朗日函数:
Figure BDA0003448930900000071
将式(12)分别对id1、id3、iq1、iq3、λ1和λ3求导,可得:
Figure BDA0003448930900000072
令上式
Figure BDA0003448930900000073
Figure BDA0003448930900000074
得到如下所示的中速区参考电流:
Figure BDA0003448930900000081
步骤3中,计算高速区d轴参考电流、q轴参考电流和励磁电流参考值的方法如下:
电机的电压极限与电流极限方程如下所示:
Figure BDA0003448930900000082
式中,usN1、usN3分别为基波平面与三次谐波平面的电机额定电压,isN1、isN3分别为基波平面与三次谐波平面的电机额定电压,Udc为直流母线电压;ud1、ud3分别为基波平面与三次谐波平面的d轴电压,uq1、uq3分别为基波平面与三次谐波平面的q轴电压,
高速稳态运行时,忽略电枢绕组上的压降,式(2)重写为:
Figure BDA0003448930900000091
由式(15)和(16),可得电压极限方程为:
Figure BDA0003448930900000092
电机运行在高速区时,在此范围内,dq轴电感之间的差值产生的磁阻力矩不容忽视,所以根据转矩方程(3)和电压方程约束条件(17),建立如下拉格朗日函数:
Figure BDA0003448930900000093
将式(18)分别对id1、id3、iq1、iq3、if1、if3、λ1、λ2、λ3和λ4求导
Figure BDA0003448930900000101
令上式
Figure BDA0003448930900000102
Figure BDA0003448930900000103
电机运行在高速区时,式(19)中第8、10个方程中(Lq1iq1)2、(Lq3iq3)2较小,可以忽略,得到如下所示的高速区参考电流:
Figure BDA0003448930900000104
其中,
Figure BDA0003448930900000111
Figure BDA0003448930900000112
Figure BDA0003448930900000113
Figure BDA0003448930900000114
k7=(5k1RsMf1-2RfLd1)[(k1Ld1+Mf1)(k1(Ld1-Lq1)+Mf1)]
k8={(k1(Ld1-Lq1)+Mf1)[5k2RsMf1+(k2Ld1m1)(5k1RsMf1-2RfLd1)]
+(k1Ld1+Mf1)(5k1RsMf1-2RfLd1)[k2(Ld1-Lq1)+ψm1]}
k9={5k2RsMf1(k2Ld1m1)[k1(Ld1-Lq1)+Mf1]
+(k2(Ld1-Lq1)+ψm1)[5k2RsMf1+(k2Ld1m1)(5k1RsMf1-2RfLd1)]}
k10=5k2RsMf1(k2Ld1m1)[k2(Ld1-Lq1)+ψm1]
Figure BDA0003448930900000115
Figure BDA0003448930900000116
k13=(5k3RsMf3-2RfLd3)[(k3Ld3+Mf3)(k3(Ld3-Lq3)+Mf3)]
k14={(k3(Ld3-Lq3)+Mf3)[5k4RsMf3+(k4Ld3m3)(5k3RsMf3-2RfLd3)]
+(k3Ld3+Mf3)(5k3RsMf3-2RfLd3)[k4(Ld3-Lq3)+ψm3]}
k15={5k4RsMf3(k4Ld3m3)[k3(Ld3-Lq3)+Mf3]
+(k4(Ld3-Lq3)+ψm3)[5k4RsMf3+(k4Ld3m3)(5k3RsMf3-2RfLd3)]}
k16=5k2RsMf1(k2Ld1m1)[k2(Ld1-Lq1)+ψm1]
通过下式求出参考励磁电流ifref1、ifref3
Figure BDA0003448930900000121
本发明的有益效果是,一种五相双励磁同步电机最大转矩铜耗比协调控制方法,通过定子电枢电流交直轴解耦与励磁电流的协调控制,能够实现电机在全速度范围内保持铜耗最小的高效运行策略,包括三次谐波电流的注入对损耗的影响有更全面的研究;能进行全速度区范围的的平稳切换,综合利用电枢电流和励磁电流对转矩与磁场的调节作用,提升了电机的综合性能。多相系统具有更高的自由度,应用于电动汽车易实现输出转矩大,效率高以及调速范围宽。在整个运行区域,驱动系统都具有更强的鲁棒性、较快的动态响应。具体优点如下:
(1)转矩与转速脉动小,驱动系统运行稳定可靠,系统动态响应快,抗扰动能力与鲁棒性强;
(2)多相系统具有更高的自由度,应用于电动汽车易实现输出转矩大,效率高以及调速范围宽;
(3)综合利用电枢电流和励磁电流对转矩与磁场的调节作用,提升了电机的综合性能;
(4)充分发挥了电机与逆变器的输出能力,实现了全速度区运行的平稳切换;
(5)在全速度运行范围,相同转速条件下,注入三次谐波电流可以提高输出转矩,提高了电机带负载能力;
(6)控制方法简单,容易实现,有利于实时环境应用;
该发明提出的控制方法使五相双励磁同步电机更适合工作在低速大转矩、高效宽调速、操作灵活的变速驱动系统领域。
附图说明
图1是本发明一种五相双励磁同步电机最大转矩铜耗比协调控制方法的流程图;
图2是本发明最大转矩铜耗比协调控制策略框图;
图3是本发明一种五相双励磁同步电机最大转矩铜耗比协调控制方法所用模型的系统框图;
图4是本发明一种五相双励磁同步电机最大转矩铜耗比协调控制方法的结构框图。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
本发明一种五相双励磁同步电机最大转矩铜耗比协调控制方法的系统框图如图3所示,该控制系统由交流电源、调压器、整流电路、稳压电容、电枢与励磁驱动电路、主功率变换器、励磁功率变换器、电流和电压传感器、五相双励磁同步电机、光电编码器、DSP控制器等组成。
交流电源给整个系统供电,经过整流电路整流后,然后进行滤波、稳压送到主、励磁功率变换器,由霍尔电压传感器采集母线电压,调理后送入控制器。主、励磁功率变换器的输出端接五相双励磁同步电机,霍尔电流互感器采集相电流和励磁电流,调理后送入DSP控制器;光电编码器采集转子位置信号,处理后送入控制器计算转子位置角与角速度。控制器输出10路开关信号驱动主功率变换器,经过励磁脉宽调制模块输出4路脉冲宽度信号驱动励磁功率变换器。
一种五相双励磁同步电机最大转矩铜耗比协调控制方法,结合图1、图2、图4所示,具体按照以下步骤实施:
步骤1:从电机主电路采集五相电流iA、iB、iC、iD、iE和励磁电流if,逆变器母线电压Udc和励磁电压Uf,,对电机位置进行初始位置检测,将编码器采集到的信号送入控制器进行处理后得到电机转子位置角θ和转速n;
步骤2:将步骤1采集得到的相电流iA、iB、iC、iD、iE经A/D转换后,然后进行坐标变换得到两相旋转坐标系下关于基波平面与三次谐波平面对应的定子d轴电流id1与id3、q轴电流iq1与iq3
步骤3:判断电机运行区域,当电机运行于低速区时,根据电磁转矩方程约束条件,利用拉格朗日乘数法寻找铜耗的最小值,计算得到关于基波平面与三次谐波平面对应的定子d轴电流、q轴电流和励磁电流参考值。
当电机运行于中速区时,由于没有弱磁或增磁的需要,采用励磁电流为零的简单控制,根据拉格朗日方程求得铜耗的最小值,计算得到中速区关于基波平面与三次谐波平面对应的定子d轴电流、q轴电流和励磁电流参考值。
当电机运行于高速区时,电机反电势接近母线电压,为了同时满足负载转矩与逆变器输出电压的限制,即满足转矩方程和电压方程的约束条件下,建立拉格朗日方程,求得铜耗的最小值,计算得到关于基波平面与三次谐波平面对应的定子d轴电流、q轴电流和励磁电流参考值。
步骤3计算低速区d轴参考电流、q轴参考电流和励磁电流参考值的方法如下:
五相双励磁同步电机在dq参考坐标系下的数学模型为:
磁链方程:
Figure BDA0003448930900000151
电压方程:
Figure BDA0003448930900000152
电磁转矩方程:
Figure BDA0003448930900000153
其中,脚标“1”和“3”分别代表基波平面和三次谐波平面成分,ψd1、ψd3为定子d轴磁链,ψq1、ψq3为定子q轴磁链,ψm1、ψm3为永磁磁链幅值,ψmf1、ψmf3为励磁绕组磁链幅值;Ld1、Ld3为d轴电感,Lq1、Lq3为q轴电感,Lf1、Lf3为励磁绕组电感,Mf1、Mf3为电枢与励磁绕组之间的互感;id1、id3为d轴电流,iq1、iq3为q轴电流,if1、if3为励磁绕组电流,ωe为电角速度;ud1、ud3为d轴电压,uq1、uq3为q轴电压,uf1、uf3为励磁绕组电压,Rs为电枢绕组电阻,Rf为励磁绕组电阻;Te1、Te3为电磁转矩,p为电机极对数。
步骤3计算低速区d轴参考电流、q轴参考电流和励磁电流参考值的方法如下:
将电机的铜耗pcu表示为d轴电流、q轴电流和励磁电流的函数,如下所示:
Figure BDA0003448930900000161
电机运行转速小于额定转速时,由于dq轴电感不同产生的磁阻力矩相比于总转矩很小,所以在此速度范围内采用d轴电流为零的控制方法来实现简单有效的运行,将总转矩方程(3)改写为
Figure BDA0003448930900000162
同时,将铜耗方程(4)改写为:
Figure BDA0003448930900000163
电机运行在低速区时,基于式(5)的电磁转矩方程约束条件,建立如下所示拉格朗日函数:
Figure BDA0003448930900000164
将式(7)分别对iq1、iq3、if1、if3、λ1和λ3求导,可得:
Figure BDA0003448930900000171
令上式
Figure BDA0003448930900000172
Figure BDA0003448930900000173
得到如下所示的低速区参考电流:
Figure BDA0003448930900000174
其中,
Figure BDA0003448930900000175
k5=k1Mf1 2,k6=2k1ψm1Mf1+k2Mf1 2
Figure BDA0003448930900000176
k9=k5Mf3 2,k10=2k3ψm3Mf3+k4Mf3 2
Figure BDA0003448930900000177
T01=(2Te1/5p)2,T03=(2Te3/15p)2
式(9)中的最后两个方程为一元四次方程,通过求解方程中的ifref1与ifref3,带入函数表达式可解除d、q轴参考电流。
步骤3计算中速区d轴参考电流、q轴参考电流和励磁电流参考值的方法如下:
电机的运行转速超过额定转速而又低于弱磁基速时,若继续使用低速区的控制方式,q轴电流明显上升,铜耗无法达到最小。这时候需要充分考虑磁阻转矩,又此阶段没有增弱磁的需求,因此为简化控制采用励磁电流为零,利用转矩电流最大比控制,中速区铜耗方程为
Figure BDA0003448930900000181
将总转矩方程(3)改写为
Figure BDA0003448930900000182
电机运行在中速区时,基于电磁转矩方程约束条件,建立如下所示拉格朗日函数:
Figure BDA0003448930900000183
将式(12)分别对id1、id3、iq1、iq3、λ1和λ3求导,可得:
Figure BDA0003448930900000191
令上式
Figure BDA0003448930900000192
Figure BDA0003448930900000193
得到如下所示的中速区参考电流:
Figure BDA0003448930900000194
步骤3计算高速区d轴参考电流、q轴参考电流和励磁电流参考值的方法如下:
电机运行转速达到弱磁基速时,电机反电势接近母线电压,受电机输入电压和电流的限制,通过弱磁控制削弱电机的气隙磁场,电机的电压极限与电流极限方程如下所示:
Figure BDA0003448930900000201
式中,usN1、usN3分别为基波平面与三次谐波平面的电机额定电压,isN1、isN3分别为基波平面与三次谐波平面的电机额定电压,Udc为直流母线电压;ud1、ud3分别为基波平面与三次谐波平面的d轴电压,uq1、uq3分别为基波平面与三次谐波平面的q轴电压,
高速稳态运行时,忽略电枢绕组上的压降,式(2)重写为:
Figure BDA0003448930900000202
由式(15)和(16),可得电压极限方程为:
Figure BDA0003448930900000203
电机运行在高速区时,在此范围内,dq轴电感之间的差值产生的磁阻力矩不容忽视,所以根据转矩方程(3)和电压方程约束条件(17),建立如下拉格朗日函数:
Figure BDA0003448930900000211
将式(18)分别对id1、id3、iq1、iq3、if1、if3、λ1、λ2、λ3和λ4求导
Figure BDA0003448930900000212
令上式
Figure BDA0003448930900000213
Figure BDA0003448930900000214
电机运行在高速区时,式(19)中第8、10个方程中(Lq1iq1)2、(Lq3iq3)2较小,可以忽略,得到如下所示的高速区参考电流:
Figure BDA0003448930900000221
其中,
Figure BDA0003448930900000222
Figure BDA0003448930900000223
Figure BDA0003448930900000224
Figure BDA0003448930900000225
k7=(5k1RsMf1-2RfLd1)[(k1Ld1+Mf1)(k1(Ld1-Lq1)+Mf1)]
k8={(k1(Ld1-Lq1)+Mf1)[5k2RsMf1+(k2Ld1m1)(5k1RsMf1-2RfLd1)]
+(k1Ld1+Mf1)(5k1RsMf1-2RfLd1)[k2(Ld1-Lq1)+ψm1]}
k9={5k2RsMf1(k2Ld1m1)[k1(Ld1-Lq1)+Mf1]
+(k2(Ld1-Lq1)+ψm1)[5k2RsMf1+(k2Ld1m1)(5k1RsMf1-2RfLd1)]}
k10=5k2RsMf1(k2Ld1m1)[k2(Ld1-Lq1)+ψm1]
Figure BDA0003448930900000226
Figure BDA0003448930900000227
k13=(5k3RsMf3-2RfLd3)[(k3Ld3+Mf3)(k3(Ld3-Lq3)+Mf3)]
k14={(k3(Ld3-Lq3)+Mf3)[5k4RsMf3+(k4Ld3m3)(5k3RsMf3-2RfLd3)]
+(k3Ld3+Mf3)(5k3RsMf3-2RfLd3)[k4(Ld3-Lq3)+ψm3]}
k15={5k4RsMf3(k4Ld3m3)[k3(Ld3-Lq3)+Mf3]
+(k4(Ld3-Lq3)+ψm3)[5k4RsMf3+(k4Ld3m3)(5k3RsMf3-2RfLd3)]}
k16=5k2RsMf1(k2Ld1m1)[k2(Ld1-Lq1)+ψm1]
通过下式求出参考励磁电流ifref1、ifref3
Figure BDA0003448930900000231
步骤4:通过步骤3,将最大转矩铜耗比协调控制策略产生的d轴参考电流,q轴参考电流分别送入PI调节模块后进行比例积分运算,分别获得在基波平面与三次谐波平面的d轴电压与q轴电压,然后将d轴电压与q轴电压进行旋转正交-静止两相变换后,获得静止两相坐标系下α轴电压和β轴电压,再通过空间矢量脉冲宽度调制模块后输出10路脉冲宽度调制信号,从而驱动主功率变换器工作;同时步骤1得到的励磁电流if,经跟随、滤波、偏置与A/D转换后和励磁电流参考值一起送入直流励磁脉宽调制模块,然后输出4路脉冲宽度信号来驱动励磁功率变换器。
现有混合励磁同步电机与多相电机矢量控制与直接转矩控制方法中,由于要兼顾谐波平面的影响,相电流要注入谐波电流电机的铜耗会增加,因此控制策略有必要考虑谐波电流对于损耗的影响,又存在电机与逆变器输出能力低、高速区损耗大,等问题。本发明的有益效果是,一种五相双励磁同步电机最大转矩铜耗比协调控制方法,通过定子电枢电流交直轴解耦与励磁电流的协调控制,能够实现电机在全速度范围内保持铜耗最小的高效运行策略,包括三次谐波电流的注入对损耗的影响有更全面的研究;能进行全速度区范围的的平稳切换,综合利用电枢电流和励磁电流对转矩与磁场的调节作用,提升了电机的综合性能。多相系统具有更高的自由度,应用于电动汽车易实现输出转矩大,效率高以及调速范围宽。在整个运行区域,驱动系统都具有更强的鲁棒性、较快的动态响应。
该发明提出的控制方法使五相双励磁同步电机更适合工作在低速大转矩、高效宽调速、操作灵活的变速驱动系统领域。

Claims (5)

1.一种五相双励磁同步电机最大转矩铜耗比协调控制方法,其特征在于,具体按照以下步骤实施:
步骤1、从五相双励磁同步电机主电路采集五相电流iA、iB、iC、iD、iE和励磁电流if,逆变器母线电压Udc和励磁电压Uf,,对电机位置进行初始位置检测,将编码器采集到的信号送入控制器进行处理后得到电机转子位置角θ和转速n;
步骤2:将步骤1采集得到的相电流iA、iB、iC、iD、iE经A/D转换后,然后进行坐标变换分别得到两相旋转坐标系下关于基波平面与三次谐波平面对应的定子d轴电流id1与id3、q轴电流iq1与iq3
步骤3:判断电机运行区域,当电机运行于低速区时,利用拉格朗日乘数法寻找铜耗的最小值,计算得到低速区关于基波平面与三次谐波平面对应的定子d轴电流、q轴电流和励磁电流参考值;
当电机运行于中速区时,采用励磁电流为零的简单控制,根据拉格朗日方程求得铜耗的最小值,计算得到中速区关于基波平面与三次谐波平面对应的定子d轴电流、q轴电流和励磁电流参考值;
当电机运行于高速区时,建立拉格朗日方程,求得铜耗的最小值,计算得到高速区关于基波平面与三次谐波平面对应的定子d轴电流、q轴电流和励磁电流参考值;
步骤4:通过步骤3,将最大转矩铜耗比协调控制策略产生的d轴参考电流,q轴参考电流分别送入PI调节模块后进行比例积分运算,分别获得在基波平面与三次谐波平面的d轴电压与q轴电压,然后将d轴电压与q轴电压进行旋转正交-静止两相变换后,获得静止两相坐标系下α轴电压和β轴电压,再通过空间矢量脉冲宽度调制模块后输出10路脉冲宽度调制信号,从而驱动主功率变换器工作;同时步骤1得到的励磁电流if,经跟随、滤波、偏置与A/D转换后和励磁电流参考值一起送入直流励磁脉宽调制模块,然后输出4路脉冲宽度信号来驱动励磁功率变换器。
2.根据权利要求1所述的一种五相双励磁同步电机最大转矩铜耗比协调控制方法,其特征在于,所述步骤1中,
五相双励磁同步电机在两相旋转坐标系下的数学模型为:
磁链方程:
Figure FDA0003448930890000021
电压方程:
Figure FDA0003448930890000022
电磁转矩方程:
Figure FDA0003448930890000031
其中,脚标“1”和“3”分别代表基波平面和三次谐波平面成分,ψd1、ψd3为定子d轴磁链,ψq1、ψq3为定子q轴磁链,ψm1、ψm3为永磁磁链幅值,ψmf1、ψmf3为励磁绕组磁链幅值;Ld1、Ld3为d轴电感,Lq1、Lq3为q轴电感,Lf1、Lf3为励磁绕组电感,Mf1、Mf3为电枢与励磁绕组之间的互感;id1、id3为d轴电流,iq1、iq3为q轴电流,if1、if3为励磁绕组电流,ωe为电角速度;ud1、ud3为d轴电压,uq1、uq3为q轴电压,uf1、uf3为励磁绕组电压,Rs为电枢绕组电阻,Rf为励磁绕组电阻;Te1、Te3为电磁转矩,p为电机极对数。
3.根据权利要求2所述的一种五相双励磁同步电机最大转矩铜耗比协调控制方法,其特征在于,所述步骤3中,
计算低速区d轴参考电流、q轴参考电流和励磁电流参考值的方法如下:
将电机的铜耗pcu表示为d轴电流、q轴电流和励磁电流的函数,如下所示:
Figure FDA0003448930890000032
电机运行转速小于额定转速时,将总转矩方程(3)改写为
Figure FDA0003448930890000033
同时,将铜耗方程(4)改写为:
Figure FDA0003448930890000034
电机运行在低速区时,基于式(5)的电磁转矩方程约束条件,建立如下所示拉格朗日函数:
Figure FDA0003448930890000041
将式(7)分别对iq1、iq3、if1、if3、λ1和λ3求导,可得:
Figure FDA0003448930890000042
令上式
Figure FDA0003448930890000043
Figure FDA0003448930890000044
得到如下所示的低速区参考电流:
Figure FDA0003448930890000045
其中,
Figure FDA0003448930890000051
k5=k1Mf1 2,k6=2k1ψm1Mf1+k2Mf1 2
Figure FDA0003448930890000052
k9=k5Mf3 2,k10=2k3ψm3Mf3+k4Mf3 2
Figure FDA0003448930890000053
T01=(2Te1/5p)2,T03=(2Te3/15p)2
式(9)中的最后两个方程为一元四次方程,通过求解方程中的ifref1与ifref3,带入函数表达式可解除d、q轴参考电流。
4.根据权利要求3所述的一种五相双励磁同步电机最大转矩铜耗比协调控制方法,其特征在于,所述步骤3中,计算中速区d轴参考电流、q轴参考电流和励磁电流参考值的方法如下:
利用转矩电流最大比控制,中速区铜耗方程为
Figure FDA0003448930890000054
将总转矩方程(3)改写为
Figure FDA0003448930890000055
电机运行在中速区时,基于电磁转矩方程约束条件,建立如下所示拉格朗日函数:
Figure FDA0003448930890000056
将式(12)分别对id1、id3、iq1、iq3、λ1和λ3求导,可得:
Figure FDA0003448930890000061
令上式
Figure FDA0003448930890000062
Figure FDA0003448930890000063
得到如下所示的中速区参考电流:
Figure FDA0003448930890000064
5.根据权利要求4所述的一种五相双励磁同步电机最大转矩铜耗比协调控制方法,其特征在于,所述步骤3中,计算高速区d轴参考电流、q轴参考电流和励磁电流参考值的方法如下:
电机的电压极限与电流极限方程如下所示:
Figure FDA0003448930890000071
式中,usN1、usN3分别为基波平面与三次谐波平面的电机额定电压,isN1、isN3分别为基波平面与三次谐波平面的电机额定电压,Udc为直流母线电压;ud1、ud3分别为基波平面与三次谐波平面的d轴电压,uq1、uq3分别为基波平面与三次谐波平面的q轴电压,
高速稳态运行时,忽略电枢绕组上的压降,式(2)重写为:
Figure FDA0003448930890000072
由式(15)和(16),可得电压极限方程为:
Figure FDA0003448930890000073
电机运行在高速区时,在此范围内,dq轴电感之间的差值产生的磁阻力矩不容忽视,所以根据转矩方程(3)和电压方程约束条件(17),建立如下拉格朗日函数:
Figure FDA0003448930890000081
将式(18)分别对id1、id3、iq1、iq3、if1、if3、λ1、λ2、λ3和λ4求导
Figure FDA0003448930890000082
令上式
Figure FDA0003448930890000083
Figure FDA0003448930890000084
电机运行在高速区时,式(19)中第8、10个方程中(Lq1iq1)2、(Lq3iq3)2较小,可以忽略,得到如下所示的高速区参考电流:
Figure FDA0003448930890000091
其中,
Figure FDA0003448930890000092
Figure FDA0003448930890000093
Figure FDA0003448930890000094
Figure FDA0003448930890000095
k7=(5k1RsMf1-2RfLd1)[(k1Ld1+Mf1)(k1(Ld1-Lq1)+Mf1)]
k8={(k1(Ld1-Lq1)+Mf1)[5k2RsMf1+(k2Ld1m1)(5k1RsMf1-2RfLd1)]+(k1Ld1+Mf1)(5k1RsMf1-2RfLd1)[k2(Ld1-Lq1)+ψm1]}
k9={5k2RsMf1(k2Ld1m1)[k1(Ld1-Lq1)+Mf1]+(k2(Ld1-Lq1)+ψm1)[5k2RsMf1+(k2Ld1m1)(5k1RsMf1-2RfLd1)]}
k10=5k2RsMf1(k2Ld1m1)[k2(Ld1-Lq1)+ψm1]
Figure FDA0003448930890000096
Figure FDA0003448930890000097
k13=(5k3RsMf3-2RfLd3)[(k3Ld3+Mf3)(k3(Ld3-Lq3)+Mf3)]
k14={(k3(Ld3-Lq3)+Mf3)[5k4RsMf3+(k4Ld3m3)(5k3RsMf3-2RfLd3)]+(k3Ld3+Mf3)(5k3RsMf3-2RfLd3)[k4(Ld3-Lq3)+ψm3]}
k15={5k4RsMf3(k4Ld3m3)[k3(Ld3-Lq3)+Mf3]+(k4(Ld3-Lq3)+ψm3)[5k4RsMf3+(k4Ld3m3)(5k3RsMf3-2RfLd3)]}
k16=5k2RsMf1(k2Ld1m1)[k2(Ld1-Lq1)+ψm1]
通过下式求出参考励磁电流ifref1、ifref3
Figure FDA0003448930890000101
CN202111668934.9A 2021-12-31 2021-12-31 一种五相双励磁同步电机最大转矩铜耗比协调控制方法 Withdrawn CN114400942A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111668934.9A CN114400942A (zh) 2021-12-31 2021-12-31 一种五相双励磁同步电机最大转矩铜耗比协调控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111668934.9A CN114400942A (zh) 2021-12-31 2021-12-31 一种五相双励磁同步电机最大转矩铜耗比协调控制方法

Publications (1)

Publication Number Publication Date
CN114400942A true CN114400942A (zh) 2022-04-26

Family

ID=81228579

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111668934.9A Withdrawn CN114400942A (zh) 2021-12-31 2021-12-31 一种五相双励磁同步电机最大转矩铜耗比协调控制方法

Country Status (1)

Country Link
CN (1) CN114400942A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116191969A (zh) * 2023-02-24 2023-05-30 山东大学 基于谐波电流注入的同步电机转矩密度提升方法及系统

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116191969A (zh) * 2023-02-24 2023-05-30 山东大学 基于谐波电流注入的同步电机转矩密度提升方法及系统
CN116191969B (zh) * 2023-02-24 2023-11-28 山东大学 基于谐波电流注入的同步电机转矩密度提升方法及系统

Similar Documents

Publication Publication Date Title
CN103647489B (zh) 一种混合励磁同步电机效率最优控制方法
CN106559026B (zh) 一种电机驱动系统的控制方法、控制装置和变频空调器
CN101855825B (zh) 马达控制装置与其控制方法
CN106788081B (zh) 一种混合励磁同步电机损耗最小直接转矩控制方法
JPH1118496A (ja) 電気車の制御装置および制御方法
CN108390602B (zh) 一种混合励磁同步电机直接预测功率控制方法
CN104767455B (zh) 一种混合励磁同步电机无位置传感器直接转矩控制方法
WO2011039993A1 (ja) 電力変換装置
CN106788026A (zh) 一种空间矢量信号注入永磁同步电机最大转矩电流比控制方法
CN103997267A (zh) 一种开绕组永磁同步电机的串联补偿直接转矩控制方法
CN113992071B (zh) 一种基于运行工况的双模式绕组电机驱动系统
CN103532466B (zh) 一种用于控制永磁同步电机转矩变化率的方法及装置
CN109347386A (zh) 一种基于svpwm的五相永磁电机最大转矩电流比容错控制方法
Ellabban et al. A comparative study of different control techniques for an induction motor fed by a Z-source inverter for electric vehicles
CN102412778A (zh) 全速度域转矩最大化的感应电机矢量控制电流分配方法
CN104682806B (zh) 一种混合励磁同步电机恒磁链控制方法
CN112886880A (zh) 三电平永磁同步电机无位置传感器模型预测电流控制方法
CN112737439B (zh) 一种高压直流起动发电机电动运行弱磁控制方法
CN107947669B (zh) 一种混合励磁同步电机非线性逆推跟踪控制方法
CN108418485B (zh) 一种隐极式混合励磁电机恒功率损耗模型预测控制方法
CN106533294B (zh) 一种基于线电压调制的永磁同步电机弱磁控制方法
CN105743330B (zh) 一种双逆变器单svm简化调制算法
CN114400942A (zh) 一种五相双励磁同步电机最大转矩铜耗比协调控制方法
CN104617849B (zh) 一种混合励磁同步电机最大输出功率控制方法
JP2009183051A (ja) 同期機の制御装置

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
WW01 Invention patent application withdrawn after publication
WW01 Invention patent application withdrawn after publication

Application publication date: 20220426