CN107276475A - 一种双电机串联缺相容错预测型直接转矩控制方法 - Google Patents
一种双电机串联缺相容错预测型直接转矩控制方法 Download PDFInfo
- Publication number
- CN107276475A CN107276475A CN201710522578.7A CN201710522578A CN107276475A CN 107276475 A CN107276475 A CN 107276475A CN 201710522578 A CN201710522578 A CN 201710522578A CN 107276475 A CN107276475 A CN 107276475A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mtd
- mtr
- msqrt
- 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
Links
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P21/00—Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P25/00—Arrangements or methods for the control of AC motors characterised by the kind of AC motor or by structural details
- H02P25/02—Arrangements 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/022—Synchronous motors
- H02P25/024—Synchronous motors controlled by supply frequency
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P2207/00—Indexing scheme relating to controlling arrangements characterised by the type of motor
- H02P2207/05—Synchronous machines, e.g. with permanent magnets or DC excitation
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Control Of Multiple Motors (AREA)
- Control Of Ac Motors In General (AREA)
- External Artificial Organs (AREA)
Abstract
本发明涉及一种双电机串联缺相容错预测型直接转矩控制方法,由逆变器、T6变换模块、旋转坐标变换模块、两台电机磁链计算模块、预测模型模块、评估函模块等环节构成。根据当前扑第k扑的电流、转子位置角、转速等信息,借助于定子磁链及电磁转矩预测模型,从损耗函数最小角度预测出下一扑第k+1扑逆变器的功率开关组合。本发明所提出的一种双电机串联缺相容错预测型直接转矩控制方法,实现了缺相后的两台电机运行相互解耦,提高了串联驱动***缺相容错运行能力,实现了两台电机电磁转矩及定子磁链的精确控制,两台电机电磁转矩及定子磁链稳态脉动大幅度降低,两台电机稳态运行更加平稳。
Description
技术领域
本发明涉及一种双电机串联缺相容错预测型直接转矩控制方法。
背景技术
单六相逆变器供电六相串联三相双永磁同步电机驱动***中两台电机之间相互解耦,而这种解耦特性的获得是依赖两台电机绕组之间的连接实现的,即把六相电机电空间对称的两相绕组尾端并联后,再与三相电机中的一相绕组串联,这样三相绕组电流均分到并联的六相电机的两相绕组中。六相电机中产生对称空间旋转磁场的电流分量不流过三相绕组;而三相电机的电流虽然流过六相绕组,但不会在六相电机中产生旋转磁场。从而实现两台电机之间的解耦控制。采用直接转矩控制策略,可以进一步提高两台电机转矩的动态控制性能,同时也可以进一步提高两台电机之间控制的可靠性。
但当六相绕组中一相绕组断路,或六相逆变桥中一个桥臂发生故障时,六相绕组中只有剩余的5相绕组可以工作,与缺相相绕组电空间对称的绕组中直接流过三相绕组相电流。显然,三相绕组电流对六相电机旋转磁场产生不利影响,如何保持两台电机之间的转矩解耦控制是期待解决的难题。
为此,本发明针对双电机串联驱动***在六相电机缺一相情况,提出一种容错预测型直接转矩控制方法。
发明内容
本发明的目的在于提供一种双电机串联缺相容错预测型直接转矩控制方法,以克服现有技术中存在的缺陷。
为实现上述目的,本发明的技术方案是:一种双电机串联缺相容错预测型直接转矩控制方法,提供一六相永磁同步电机以及一三相永磁同步电机,提供一预测型直接转矩控制***,该***包括:T6变换模块、第一旋转坐标变换模块、第二旋转坐标变换模块、六相电机定子磁链计算模块、三相电机定子磁链计算模块、预测模型模块以及评估函数模块;按照如下步骤实现:
步骤S1:所述T6变换模块获取待采样逆变器中六相功率开关管的六相电流isA~isF,并经所述T6变换模块变换后,输出αβ坐标平面中isα(k)、isβ(k)、xy坐标平面中isx(k)、isy(k)以及O1O2坐标平面中io1(k)、io2(k);
步骤S2:将isα(k)、isβ(k)、所述六相永磁同步电机的转子位置角θr1(k)传输至至所述第一旋转坐标变换模块,经坐标变换后输出d1q1坐标平面中id1(k)、iq1(k);将isx(k)、isy(k)、所述三相永磁同步电机的转角位置θr2(k)传输至所述第二旋转坐标变换模块,经坐标变换后输出d2q2坐标平面中id2(k)、iq2(k);
步骤S3:将id1(k)、iq1(k)传输至所述六相电机定子磁链计算模块,输出六相电机d1q1坐标平面中定子磁链ψd1(k)、ψq1(k);将id2(k)、iq2(k)传输至所述三相电机定子磁链计算模块,输出三相电机d2q2坐标平面中定子磁链ψd2(k)、ψq2(k);
步骤S4:把ψd1(k)、ψq1(k)、ψd2(k)、ψq2(k)、id1(k)、iq1(k)、id2(k)、iq2(k)、θr1(k)、θr2(k)、直流母线电压UDC传输至所述预测模型模块,输出两台电机电磁转矩第k+1扑预测值Te1(k+1)、Te2(k+1)以及两台电机定子磁链幅值第k+1扑预测值|ψs1(k+1)|、|ψs2(k+1)|;
步骤S5:将两台电机的给定电机电磁转矩两台电机的电机给定定子磁链两台电机电磁转矩第k+1扑预测值Te1(k+1)、Te2(k+1)、两台电机定子磁链幅值第k+1扑预测值|ψs1(k+1)|、|ψs2(k+1)|传输至所述评估函数模块,计算出对应于Sb~Sf的评估函数值cost;
步骤S6:获取cost值最小时的Sb~Sf取值,并根据Sb~Sf取值控制剩余健康五相功率开关管,从而使两台电机电磁转矩、定子磁链以最小的脉动进行闭环控制。
在本发明一实施例中,在所述步骤S1中,所述T6变换模块按照如下方式进行变换:
在本发明一实施例中,在所述步骤S2中,所述第一旋转坐标变换模块按照
如下方式进行坐标变换:
所述第二旋转坐标变换模块按照如下方式进行坐标变换:
在本发明一实施例中,在所述步骤S3中,所述六相电机d1q1坐标平面中定子磁链ψd1(k)ψq1(k)通过如下方式获取:
其中,ψf1为六相永磁同步电机的转子磁链矢量,Ld1=Lsσ1+3Lsm1+3Lrs1、Lq1=Lsσ1+3Lsm1-3Lrs1,Lsm1=(Ldm1+Lqm1)/2,Lrs1=(Ldm1-Lqm1)/2,Ldm1、Lqm1为六相电机相绕组主磁通直轴电感、交轴电感,Lsσ1为六相电机相绕组漏电感,Ld1、Lq1分别为六相电机直、交轴电感;
所述三相电机d2q2坐标平面中定子磁链ψd2(k)ψq2(k)通过如下方式获取:
其中,ψf2为三相永磁同步电机的转子磁链矢量,Ld2=Lsσ1+2Lsσ2+3Lsm2+3Lrs2,Lq2=Lsσ1+2Lsσ2+3Lsm2-3Lrs2,Lsσ2为三相电机相绕组漏电感,Lsm2=(Ldm2+Lqm2)/2,Lrs2=(Ldm2-Lqm2)/2,Ldm2、Lqm2为三相电机相绕组主磁通直轴电感、交轴电感,Ld2、Lq2分别为三相电机直、交轴电感。
在本发明一实施例中,在所述步骤S4中,还包括如下步骤:
步骤S41:将把id1(k)、iq1(k)、id2(k)、iq2(k)、θr1(k)、θr2(k)、ωr1(k)、ωr2(k)、i02(k)传输至一计算模块,并按照如下方式计算:
其中,Rs1为绕组电阻,ωr1(k)为六相电机转子旋转的电角速度,ωr2(k)为三相电机转子旋转的电角速度;
步骤S42:将一组当前Sb~Sf值、直流母线电压UDC送给电压矢量计算环节,通过如下方式获取u′sα、usβ、u′sx、usy:
步骤S43:把id1(k)、iq1(k)、id2(k)、iq2(k)、θr1(k)、θr2(k)、ωr1(k)、ωr2(k)、u′sα、usβ、u′sx、usy、ψd1(k)、ψq1(k)、ψd2(k)、ψq2(k)对应传输至两台电机定子磁链预测模块,输出第k+1扑对应Sb~Sf的两台电机定子磁链预测值ψd1(k+1)、ψq1(k+1)、ψd2(k+1)、ψq2(k+1),其计算公式如下:
步骤S44:将ψd1(k+1)ψq1(k+1)、ψd2(k+1)ψq2(k+1)传输至一磁链幅值计算模块,按照如下方式输出两台电机磁链幅值|ψs1(k+1)|、|ψs2(k+1)|:
步骤S44:把ψd1(k+1)、ψq1(k+1)、ψd2(k+1)、ψq2(k+1)传输至一电磁转矩预测模块,按照如下方式输出两台电机电磁转矩预测值Te1(k+1)、Te2(k+1):
在本发明一实施例中,在所述步骤S5中,所述评估函数值cost:
其中,k1、k2、k3、k4为损耗函数系数。
在本发明一实施例中,当两台电机的控制量为电磁转矩时,所述给定电机电磁转矩由***给定;当两台电机的控制量为转速时,所述给定电机电磁转矩由两台电机速度控制器分别输出转矩给定;当两台电机的控制量为转子位置角时,所述给定电机电磁转矩由两台电机位置控制器输出转矩给定。
在本发明一实施例中,在所述步骤S6中,把Sb~Sf值增加1,再从所述步骤S4开始计算;当Sb~Sf的取值完全后,找出对应cost值最小的Sb~Sf。
在本发明一实施例中,Sb~Sf取值为Sb~Sf=00000~11111。
相较于现有技术,本发明具有以下有益效果:
1)由于利用了考虑缺相相绕组电压的逆变器输出电压矢量对考虑缺相相绕组磁链的两台电机定子磁链及电磁转矩进行解耦控制,实现了缺相后的两台电机运行相互解耦,提高了串联驱动***缺相容错运行能力。
2)借助于电磁转矩和定子磁链的预测控制策略,加在两台电机上的电压矢量总是最优,从而实现了两台电机电磁转矩及定子磁链的精确控制,两台电机电磁转矩及定子磁链稳态脉动大幅度降低,两台电机稳态运行更加平稳。
附图说明
图1为本发明一实施例中双电机串联缺相容错预测型直接转矩控制方法的驱动***硬件结构图。
图2为本发明一实施例中六相逆变器供电六相串联三相双永磁同步电机驱动电路连接示意图。
图3(a)为本发明一实施例中六相电机机电能量转换坐标平面示意图。
图3(b)为本发明一实施例中三相电机机电能量转换坐标平面示意图。
图4为本发明一实施例中串联缺相容错预测型直接转矩控制结构示意图。
图5为本发明一实施例中预测模型模块示意图。
具体实施方式
下面结合附图,对本发明的技术方案进行具体说明。
如图1所示,本实施中双电机串联缺相容错预测型直接转矩控制方法的驱动***包括:整流电路、滤波电容、六相逆变器、偏置60度六相对称绕组永磁同步电机、三相永磁同步电机、六相绕组电流采集电路、两台电机转子位置角采集电路、隔离驱动、中央控制器、人机接口等。也可以采用合适的直流电源提供六相逆变器直流母线电压。逆变器中功率管采用IGBT或MOFET,中央控制器采用DSP或单片机。绕组电流采集电路采用霍尔电流传感器与运算放大器相结合方式构成,也可以采用绕组串功率电阻后接差分运算放大器相结合方式构成。采用霍尔方案可以有效实现控制回路与主回路的电气隔离,采用绕组串功率电阻方案可以降低驱动***成本。转子位置角采集电路可以采用旋转编码器后接电平转换电路构成,也可以采用旋转变压器后接解码电路构成,其中前者成本较低,但位置角采样精度受编码器线数限制,而后者成本较高,但位置角采样精度较高。绕组电流采集电路和转子位置角采集电路弱电压信号送到中央控制器A/D转换模块。根据取得的信号和本发明的容错预测型直接转矩控制策略计算出应发出的控制信号,经由隔离驱动去控制逆变器中的功率开关管的开关动作。
进一步的,六相逆变器供电六相串联三相双永磁同步电机驱动电路连接示意如图2所示,其中A~F为六相电机的相绕组,U~W为三相电机的相绕组。六相电机电空间对称两相绕组尾端并联后再与三相绕组串联,从而保证六相电机中实现机电能量转换的电流分量不流过三相绕组;任意一相三相电流均分到两个并联六相绕组中,在六相电机中产生的合成磁动势为零。
进一步的,本实施例中,分析中实现机电能量转换坐标平面定义如图3(a)以及图3(b)所示。图中,αβ和xy为静止坐标系,d1q1和d2q2为六相和三相电机转子同步旋转坐标系。θr1为d1轴和α轴之间的夹角即六相电机转子旋转的电角度,ωr1为六相电机转子旋转的电角速度。ψs1、ψf1、u s1、is1分别为六相电机实现机电能量转换的定子磁链矢量、转子磁链矢量、定子电压矢量以及定子电流矢量,这些矢量及驱动***输入电压、电流在d1轴、q1轴、α轴、β轴上的投影分别用下标“d1”、“q1”、“α”、“β”标示。δ1为六相电机定子磁链和转子磁链之间的夹角。同理定义三相电机中有关变量如图3(b)所示。由于驱动***有6个自由度,除了两台电机机电能量转换平面中电磁转矩和磁链共计4个自由度外,还存在非机电能量转换轴系o 1o 2上的2个自由度,本发明称之为零序轴系。机电能量转换平面中电磁转矩和磁链的脉动直接影响两台电机切向旋转的平稳性,而轴系o 1o 2上变量直接影响六相电机铜损耗、铁损耗等。A~F相逆变桥臂上下功率管互锁开关,且用一个开关变量表示功率管开关情况,即Si=1表示上桥功率管导通,下桥功率管关闭;Si=0表示上桥功率管关断,下桥功率管导通(i=A~F)。
进一步的,在本实施例中,可以根据绕组电压、电流、磁链、电感之间关系建立双电机在自然坐标系A~F内的数学模型。为了揭示双电机之间的耦合关系,利用下式(1)恒功率变换矩阵T6将自然坐标系A~F内的数学模型变换到αβ-xy-o1o 2静止坐标系。
六相电机αβ平面内的电压、磁链、电磁转矩平衡方程式如下:
Te1=p1(ψsαiβ-ψsβiα) (4)
其中,Rs1为绕组电阻,Lsσ1为每相绕组的自漏感,Lsm1=(Ldm1+Lqm1)/2,Lrs1=(Ldm1-Lqm1)/2,Ldm1、Lqm1为六相电机相绕组主磁通直轴电感、交轴电感,p1为磁极对数。
三相电机xy平面内的电压、磁链、电磁转矩平衡方程式如下:
Te2=p2(ψsxiy-ψsyix) (7)
其中,有关变量意义类似于六相电机。
零序轴系中电压平衡方程式如下:
利用θr1角,可以把式(2)-(4)变换到坐标系d1q1中得:
Te1=p1(ψd1iq1-ψq1id1) (11)
其中,Ld1=Lsσ1+3Lsm1+3Lrs1、Lq1=Lsσ1+3Lsm1-3Lrs1。
利用θr2角,可以把式(5)-(7)变换到坐标系d2q2中得:
Te2=p2(ψd2iq2-ψq2id2) (14)
其中,Ld2=Lsσ1+2Lsσ2+3Lsm2+3Lrs2、Lq2=Lsσ1+2Lsσ2+3Lsm2-3Lrs2。
当六相电机A相绕组断路或A相逆变桥故障后,六相中的A相不流过电流,只有剩余的B~F五相工作,只有四个可控自由度,无需对零序电流进行控制。为了实现两台电机的解耦控制,假设A相逆变桥始终输出六相A和三相U串联电压u0,这样六个逆变桥输出电压如下:
尽管六相中A相没有电流流过,但若考虑A相绕组后,双电机串联***输入电压之和为零特性,可以求出UNO,然后代入式(15),再利用T6把式(15)变换至αβ-xy-o 1o 2中如下:
令:
则式(16)进一步变化为:
同样可以把双电机串联***输入电流[isA isB isC isD isE isF]借助T6矩阵变换至αβ-xy-o 1o 2中,且isA=0,可以求解出:
这样,零序电压uo2变形为:
若用同步旋转坐标系中变量表示isα、isx,则零序电压uo2变形为:
其中,
其中,
把式(16)旋转变换到同步旋转坐标系中,得:
根据式(9)(10)(12)(13)可得:
连列式(23)(24)可以求解dq轴系磁链的微分表达式,并据此用离散形式表示如下:
由于求解上述4x4矩阵的逆很耗费时间,而实际***中六相电机漏感较小。为此,实际编程中可以忽略六相电机漏感对三相电机磁链的贡献。
这样式(25)进一步简化为:
Ts为数字控制周期。
则第k+1扑两台电机定子磁链幅值为:
根据式(11)(14)可以求出下一周期两台电机电磁转矩的预测值为:
本实施例控制的目的希望两台电机每一个数字控制周期所加电压矢量实现两台电机电磁转矩和定子磁链幅值误差最小,为此,本实施例预测控制中选择如下评估函数:
显然,所选择的开关组合Sb~Sf应该使得式(31)值最小。
本实施例提出一种双电机串联缺相容错预测型直接转矩控制方法,目的有两方面:一是实现双电机串联***在一相逆变桥臂故障或六相电机中一相绕组断路情况下,仍然保持两台电机高性能的转矩解耦控制,从而提高串联驱动***运行的可靠性;二是实现两台电机的电磁转矩和定子磁链的精确控制,进一步改善两台电机运行的平稳性。在构建控制方法过程,为了实现两台电机电磁转矩、定子磁链之间的解耦,考虑了缺相相绕组耦合磁链对定子磁链的贡献,即六相电机定子磁链包括缺相相绕组耦合磁链。同样,为了利用剩余5相逆变桥输出电压矢量精确控制两台电机变量,逆变桥输出电压矢量中考虑了缺相相绕组。
进一步的,尽管六相绕组缺了一相,但在构建两台电机定子磁链过程中,仍然考虑了缺相相绕组磁链;同时构建出含有缺相相绕组电压的逆变器电压矢量。借助于电磁转矩和定子磁链预测控制策略,利用最优电压矢量实现两台电机电磁转矩及定子磁链最低脉动的解耦控制,从而提高了串联驱动***的可靠性。假设当前控制处于第k扑,关键是要根据预测控制方法,找出下一扑k+1逆变器的开关组合Sb~Sf。
所提的预测型直接转矩控制***结构框图如图4所示,由逆变器、T6变换、旋转坐标变换、两台电机磁链计算环节、预测模型、评估函数、两台电机等环节构成。预测模型如图5所示。根据当前扑(第k扑)的电流、转子位置角、转速等信息,借助于定子磁链及电磁转矩预测模型,从损耗函数最小角度预测出下一扑(第k+1扑)逆变器的功率开关组合。
具体步骤如下:
步骤(1):把采样得到的逆变器输出六相电流isA~isF送给T6变换模块,输出αβ坐标平面中isα(k)isβ(k)、xy坐标平面中isx(k)isy(k)、o1o2坐标平面中io1(k)io2(k),其计算数学模型如下:
步骤(2):把isα(k)isβ(k)、六相永磁同步电机(该电机简记为PMSM6)转子位置角θr1(k)送给旋转坐标变换模块,输出d1q1坐标平面中id1(k)iq1(k);把isx(k)isy(k)、三相永磁同步电机(该电机简记为PMSM3)转子位置角θr2(k)送给旋转坐标变换模块,输出d2q2坐标平面中id2(k)iq2(k),其计算公式如下:
步骤(3):把id1(k)iq1(k)、id2(k)iq2(k)分别送给六相电机和三相电机定子磁链计算模块,分别输出六相电机d1q1坐标平面中定子磁链ψd1(k)ψq1(k)、三相电机d2q2坐标平面中定子磁链ψd2(k)ψq2(k),其计算公式如下:
步骤(4):把ψd1(k)ψq1(k)、ψd2(k)ψq2(k)、id1(k)iq1(k)、id2(k)iq2(k)、θr1(k)、θr2(k)、直流母线电压UDC送给预测模型环节,输出两台电机电磁转矩第k+1扑预测值Te1(k+1)Te2(k+1)、两台电机定子磁链幅值第k+1扑预测值ψs1(k+1)ψs2(k+1);
其中,步骤(4)进一步包括如下步骤
步骤(4.1):把id1(k)iq1(k)、id2(k)iq2(k)、θr1(k)、θr2(k)、ωr1(k)、ωr2(k)、i02(k)送给计算环节,其计算公式如下:
步骤(4.2):把一组Sb~Sf值、直流母线电压UDC送给电压矢量计算环节,输出u′sαusβu′sxusy,其计算公式如下:
步骤(4.3):把id1(k)iq1(k)、id2(k)iq2(k)、θr1(k)、θr2(k)、ωr1(k)、ωr2(k)、u′sαusβu′sxusy、ψd1(k)ψq1(k)、ψd2(k)ψq2(k)送给两台电机定子磁链预测环节,输出对应Sb~Sf的两台电机定子磁链预测值ψd1(k+1)ψq1(k+1)、ψd2(k+1)ψq2(k+1),其计算公式如下:
步骤(4.4):把ψd1(k+1)ψq1(k+1)、ψd2(k+1)ψq2(k+1)送给磁链幅值计算环节,输出两台电机磁链幅值|ψs1(k+1)|、|ψs2(k+1)|,其计算公式如下:
步骤(4.5):把ψd1(k+1)ψq1(k+1)、ψd2(k+1)ψq2(k+1)送给电磁转矩预测环节,输出两台电机电磁转矩预测值Te1(k+1)、Te2(k+1),其计算公式如下:
步骤(5):把两台电机转矩给定两台电机定子磁链给定Te1(k+1)Te2(k+1)、ψs1(k+1)ψs2(k+1)送给评估函数模块,计算出对应于Sb~Sf的评估函数值cost,其计算公式如下:
其中,k1、k2、k3、k4为损耗函数系数。
进一步的,把Sb~Sf值增加1,再从步骤(4)开始计算。当Sb~Sf的取值完全后,找出对应cost值最小的Sb~Sf,并根据该组Sb~Sf取值,通过逆变器输出电压加在串联电机驱动***上,实现k+1扑两台电机电磁转矩、定子磁链以最小的脉动闭环控制。
其中,步骤(5)中电磁转矩给定和视具体两台电机控制变量而定。若控制的是电磁转矩,则***直接给定该值;若控制的是转速,则两台电机速度控制器(例如PI控制器)分别输出转矩给定和若控制的是转子位置角,则两台电机位置控制器输出即为转矩给定和
以上是本发明的较佳实施例,凡依本发明技术方案所作的改变,所产生的功能作用未超出本发明技术方案的范围时,均属于本发明的保护范围。
Claims (9)
1.一种双电机串联缺相容错预测型直接转矩控制方法,提供一六相永磁同步电机以及一三相永磁同步电机,其特征在于,提供一预测型直接转矩控制***,该***包括:T6变换模块、第一旋转坐标变换模块、第二旋转坐标变换模块、六相电机定子磁链计算模块、三相电机定子磁链计算模块、预测模型模块以及评估函数模块;按照如下步骤实现:
步骤S1:所述T6变换模块获取待采样逆变器中六相功率开关管的六相电流isA~isF,并经所述T6变换模块变换后,输出αβ坐标平面中isα(k)、isβ(k)、xy坐标平面中isx(k)、isy(k)以及O1O2坐标平面中io1(k)、io2(k);
步骤S2:将isα(k)、isβ(k)、所述六相永磁同步电机的转子位置角θr1(k)传输至至所述第一旋转坐标变换模块,经坐标变换后输出d1q1坐标平面中id1(k)、iq1(k);将isx(k)、isy(k)、所述三相永磁同步电机的转角位置θr2(k)传输至所述第二旋转坐标变换模块,经坐标变换后输出d2q2坐标平面中id2(k)、iq2(k);
步骤S3:将id1(k)、iq1(k)传输至所述六相电机定子磁链计算模块,输出六相电机d1q1坐标平面中定子磁链ψd1(k)、ψq1(k);将id2(k)、iq2(k)传输至所述三相电机定子磁链计算模块,输出三相电机d2q2坐标平面中定子磁链ψd2(k)、ψq2(k);
步骤S4:把ψd1(k)、ψq1(k)、ψd2(k)、ψq2(k)、id1(k)、iq1(k)、id2(k)、iq2(k)、θr1(k)、θr2(k)、直流母线电压UDC传输至所述预测模型模块,输出两台电机电磁转矩第k+1扑预测值Te1(k+1)、Te2(k+1)以及两台电机定子磁链幅值第k+1扑预测值|ψs1(k+1)|、|ψs2(k+1)|;
步骤S5:将两台电机的给定电机电磁转矩两台电机的电机给定定子磁链两台电机电磁转矩第k+1扑预测值Te1(k+1)、Te2(k+1)、两台电机定子磁链幅值第k+1扑预测值|ψs1(k+1)|、|ψs2(k+1)|传输至所述评估函数模块,计算出对应于Sb~Sf的评估函数值cost;
步骤S6:获取cost值最小时的Sb~Sf取值,并根据Sb~Sf取值控制剩余健康五相功率开关管,从而使两台电机电磁转矩、定子磁链以最小的脉动进行闭环控制。
2.根据权利要求1所述的一种双电机串联缺相容错预测型直接转矩控制方法,其特征在于,在所述步骤S1中,所述T6变换模块按照如下方式进行变换:
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>i</mi>
<mrow>
<mi>s</mi>
<mi>&alpha;</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>i</mi>
<mrow>
<mi>s</mi>
<mi>&beta;</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>i</mi>
<mrow>
<mi>s</mi>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>i</mi>
<mrow>
<mi>s</mi>
<mi>y</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>i</mi>
<mrow>
<mn>01</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>i</mi>
<mrow>
<mn>02</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<msqrt>
<mn>3</mn>
</msqrt>
</mfrac>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mn>0.5</mn>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mn>0.5</mn>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mn>0.5</mn>
</mrow>
</mtd>
<mtd>
<mn>0.5</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mfrac>
<msqrt>
<mn>3</mn>
</msqrt>
<mn>2</mn>
</mfrac>
</mtd>
<mtd>
<mfrac>
<msqrt>
<mn>3</mn>
</msqrt>
<mn>2</mn>
</mfrac>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mfrac>
<msqrt>
<mn>3</mn>
</msqrt>
<mn>2</mn>
</mfrac>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mfrac>
<msqrt>
<mn>3</mn>
</msqrt>
<mn>2</mn>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mn>0.5</mn>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mn>0.5</mn>
</mrow>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mn>0.5</mn>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mn>0.5</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mfrac>
<msqrt>
<mn>3</mn>
</msqrt>
<mn>2</mn>
</mfrac>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mfrac>
<msqrt>
<mn>3</mn>
</msqrt>
<mn>2</mn>
</mfrac>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mfrac>
<msqrt>
<mn>3</mn>
</msqrt>
<mn>2</mn>
</mfrac>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mfrac>
<msqrt>
<mn>3</mn>
</msqrt>
<mn>2</mn>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mfrac>
<mn>1</mn>
<msqrt>
<mn>2</mn>
</msqrt>
</mfrac>
</mtd>
<mtd>
<mfrac>
<mn>1</mn>
<msqrt>
<mn>2</mn>
</msqrt>
</mfrac>
</mtd>
<mtd>
<mfrac>
<mn>1</mn>
<msqrt>
<mn>2</mn>
</msqrt>
</mfrac>
</mtd>
<mtd>
<mfrac>
<mn>1</mn>
<msqrt>
<mn>2</mn>
</msqrt>
</mfrac>
</mtd>
<mtd>
<mfrac>
<mn>1</mn>
<msqrt>
<mn>2</mn>
</msqrt>
</mfrac>
</mtd>
<mtd>
<mfrac>
<mn>1</mn>
<msqrt>
<mn>2</mn>
</msqrt>
</mfrac>
</mtd>
</mtr>
<mtr>
<mtd>
<mfrac>
<mn>1</mn>
<msqrt>
<mn>2</mn>
</msqrt>
</mfrac>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<msqrt>
<mn>2</mn>
</msqrt>
</mfrac>
</mrow>
</mtd>
<mtd>
<mfrac>
<mn>1</mn>
<msqrt>
<mn>2</mn>
</msqrt>
</mfrac>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<msqrt>
<mn>2</mn>
</msqrt>
</mfrac>
</mrow>
</mtd>
<mtd>
<mfrac>
<mn>1</mn>
<msqrt>
<mn>2</mn>
</msqrt>
</mfrac>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<msqrt>
<mn>2</mn>
</msqrt>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>i</mi>
<mrow>
<mi>s</mi>
<mi>A</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>i</mi>
<mrow>
<mi>s</mi>
<mi>B</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>i</mi>
<mrow>
<mi>s</mi>
<mi>C</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>i</mi>
<mrow>
<mi>s</mi>
<mi>D</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>i</mi>
<mrow>
<mi>s</mi>
<mi>E</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>i</mi>
<mrow>
<mi>s</mi>
<mi>F</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>.</mo>
</mrow>
3.根据权利要求1所述的一种双电机串联缺相容错预测型直接转矩控制方法,其特征在于,在所述步骤S2中,所述第一旋转坐标变换模块按照如下方式进行坐标变换:
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>i</mi>
<mrow>
<mi>d</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>i</mi>
<mrow>
<mi>q</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>cos&theta;</mi>
<mrow>
<mi>r</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>sin&theta;</mi>
<mrow>
<mi>r</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>sin&theta;</mi>
<mrow>
<mi>r</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>cos&theta;</mi>
<mrow>
<mi>r</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>i</mi>
<mrow>
<mi>s</mi>
<mi>&alpha;</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>i</mi>
<mrow>
<mi>s</mi>
<mi>&beta;</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
</mrow>
所述第二旋转坐标变换模块按照如下方式进行坐标变换:
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>i</mi>
<mrow>
<mi>d</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>i</mi>
<mrow>
<mi>q</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>cos&theta;</mi>
<mrow>
<mi>r</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>sin&theta;</mi>
<mrow>
<mi>r</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>sin&theta;</mi>
<mrow>
<mi>r</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>cos&theta;</mi>
<mrow>
<mi>r</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>i</mi>
<mrow>
<mi>s</mi>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>i</mi>
<mrow>
<mi>s</mi>
<mi>y</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>.</mo>
</mrow>
4.根据权利要求1所述的一种双电机串联缺相容错预测型直接转矩控制方法,其特征在于,在所述步骤S3中,所述六相电机d1q1坐标平面中定子磁链ψd1(k)、ψq1(k)通过如下方式获取:
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>&psi;</mi>
<mrow>
<mi>d</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>&psi;</mi>
<mrow>
<mi>q</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>L</mi>
<mrow>
<mi>d</mi>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<msub>
<mi>L</mi>
<mrow>
<mi>q</mi>
<mn>1</mn>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>i</mi>
<mrow>
<mi>d</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>i</mi>
<mrow>
<mi>q</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>+</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msqrt>
<mn>3</mn>
</msqrt>
<msub>
<mi>&psi;</mi>
<mrow>
<mi>f</mi>
<mn>1</mn>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
</mrow>
其中,ψf1为六相永磁同步电机的转子磁链矢量,Ld1=Lsσ1+3Lsm1+3Lrs1、Lq1=Lsσ1+3Lsm1-3Lrs1,Lsm1=(Ldm1+Lqm1)/2,Lrs1=(Ldm1-Lqm1)/2,Ldm1、Lqm1为六相电机相绕组主磁通直轴电感、交轴电感,Lsσ1为六相电机相绕组漏电感,Ld1、Lq1分别为六相电机直、交轴电感;
所述三相电机d2q2坐标平面中定子磁链ψd2(k)、ψq2(k)通过如下方式获取:
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>&psi;</mi>
<mrow>
<mi>d</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>&psi;</mi>
<mrow>
<mi>q</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>L</mi>
<mrow>
<mi>d</mi>
<mn>2</mn>
</mrow>
</msub>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<msub>
<mi>L</mi>
<mrow>
<mi>q</mi>
<mn>2</mn>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>i</mi>
<mrow>
<mi>d</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>i</mi>
<mrow>
<mi>q</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>+</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msqrt>
<mn>3</mn>
</msqrt>
<msub>
<mi>&psi;</mi>
<mrow>
<mi>f</mi>
<mn>2</mn>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
</mrow>
其中,ψf2为三相永磁同步电机的转子磁链矢量,Ld2=Lsσ1+2Lsσ2+3Lsm2+3Lrs2,Lq2=Lsσ1+2Lsσ2+3Lsm2-3Lrs2,Lsσ2为三相电机相绕组漏电感,Lsm2=(Ldm2+Lqm2)/2,Lrs2=(Ldm2-Lqm2)/2,Ldm2、Lqm2为三相电机相绕组主磁通直轴电感、交轴电感,Ld2、Lq2分别为三相电机直、交轴电感。
5.根据权利要求1所述的一种双电机串联缺相容错预测型直接转矩控制方法,其特征在于,在所述步骤S4中,还包括如下步骤:
步骤S41:将把id1(k)、iq1(k)、id2(k)、iq2(k)、θr1(k)、θr2(k)、ωr1(k)、ωr2(k)、i02(k)传输至一计算模块,并按照如下方式计算:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>u</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>o</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>s</mi>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>i</mi>
<mrow>
<mn>02</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>-</mo>
<msqrt>
<mn>2</mn>
</msqrt>
<msub>
<mi>L</mi>
<mrow>
<mi>s</mi>
<mi>&sigma;</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>&lsqb;</mo>
<mo>-</mo>
<msub>
<mi>i</mi>
<mrow>
<mi>d</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<msub>
<mi>&omega;</mi>
<mrow>
<mi>r</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mrow>
<mi>r</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>i</mi>
<mrow>
<mi>q</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<msub>
<mi>&omega;</mi>
<mrow>
<mi>r</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mrow>
<mi>r</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>i</mi>
<mrow>
<mi>d</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<msub>
<mi>&omega;</mi>
<mrow>
<mi>r</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mrow>
<mi>r</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>i</mi>
<mrow>
<mi>q</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<msub>
<mi>&omega;</mi>
<mrow>
<mi>r</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mrow>
<mi>r</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>;</mo>
</mrow>
其中,Rs1为绕组电阻,ωr1(k)为六相电机转子旋转的电角速度,ωr2(k)为三相电机转子旋转的电角速度;
步骤S42:将一组当前Sb~Sf值、直流母线电压UDC送给电压矢量计算环节,通过如下方式获取u′sα、usβ、u′sx、usy:
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msubsup>
<mi>u</mi>
<mrow>
<mi>s</mi>
<mi>&alpha;</mi>
</mrow>
<mo>&prime;</mo>
</msubsup>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>u</mi>
<mrow>
<mi>s</mi>
<mi>&beta;</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>u</mi>
<mrow>
<mi>s</mi>
<mi>x</mi>
</mrow>
<mo>&prime;</mo>
</msubsup>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>u</mi>
<mrow>
<mi>s</mi>
<mi>y</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>u</mi>
<mrow>
<mi>s</mi>
<mi>o</mi>
<mn>1</mn>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>u</mi>
<mrow>
<mi>s</mi>
<mi>o</mi>
<mn>2</mn>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mfrac>
<msqrt>
<mn>3</mn>
</msqrt>
<mn>2</mn>
</mfrac>
<msub>
<mi>U</mi>
<mrow>
<mi>D</mi>
<mi>C</mi>
</mrow>
</msub>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>S</mi>
<mi>b</mi>
</msub>
<mo>+</mo>
<msub>
<mi>S</mi>
<mi>c</mi>
</msub>
<mo>+</mo>
<msub>
<mi>S</mi>
<mi>e</mi>
</msub>
<mo>-</mo>
<msub>
<mi>S</mi>
<mi>f</mi>
</msub>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0.5</mn>
<msub>
<mi>U</mi>
<mrow>
<mi>D</mi>
<mi>C</mi>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>S</mi>
<mi>b</mi>
</msub>
<mo>+</mo>
<msub>
<mi>S</mi>
<mi>c</mi>
</msub>
<mo>-</mo>
<msub>
<mi>S</mi>
<mi>e</mi>
</msub>
<mo>-</mo>
<msub>
<mi>S</mi>
<mi>f</mi>
</msub>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mfrac>
<mn>1</mn>
<mrow>
<mn>2</mn>
<msqrt>
<mn>3</mn>
</msqrt>
</mrow>
</mfrac>
<msub>
<mi>U</mi>
<mrow>
<mi>D</mi>
<mi>C</mi>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>S</mi>
<mi>b</mi>
</msub>
<mo>-</mo>
<mn>3</mn>
<msub>
<mi>S</mi>
<mi>c</mi>
</msub>
<mo>+</mo>
<mn>4</mn>
<msub>
<mi>S</mi>
<mi>d</mi>
</msub>
<mo>-</mo>
<mn>3</mn>
<msub>
<mi>S</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<msub>
<mi>S</mi>
<mi>f</mi>
</msub>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0.5</mn>
<msub>
<mi>U</mi>
<mrow>
<mi>D</mi>
<mi>C</mi>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>S</mi>
<mi>b</mi>
</msub>
<mo>-</mo>
<msub>
<mi>S</mi>
<mi>c</mi>
</msub>
<mo>+</mo>
<msub>
<mi>S</mi>
<mi>e</mi>
</msub>
<mo>-</mo>
<msub>
<mi>S</mi>
<mi>f</mi>
</msub>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>u</mi>
<mrow>
<mi>o</mi>
<mn>2</mn>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
</mrow>
步骤S43:把id1(k)、iq1(k)、id2(k)、iq2(k)、θr1(k)、θr2(k)、ωr1(k)、ωr2(k)、u′sα、usβ、u′sx、usy、ψd1(k)、ψq1(k)、ψd2(k)、ψq2(k)对应传输至两台电机定子磁链预测模块,输出第k+1扑对应Sb~Sf的两台电机定子磁链预测值ψd1(k+1)、ψq1(k+1)、ψd2(k+1)、ψq2(k+1),其计算公式如下:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>&psi;</mi>
<mrow>
<mi>d</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>&psi;</mi>
<mrow>
<mi>q</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>&psi;</mi>
<mrow>
<mi>d</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>&psi;</mi>
<mrow>
<mi>q</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>&psi;</mi>
<mrow>
<mi>d</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>&psi;</mi>
<mrow>
<mi>q</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>&psi;</mi>
<mrow>
<mi>d</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>&psi;</mi>
<mrow>
<mi>q</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>+</mo>
<msub>
<mi>T</mi>
<mi>s</mi>
</msub>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mn>2</mn>
<mfrac>
<msub>
<mi>L</mi>
<mrow>
<mi>s</mi>
<mi>&sigma;</mi>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>L</mi>
<mrow>
<mi>q</mi>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<msup>
<mi>cos</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mrow>
<mi>r</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mn>2</mn>
<mfrac>
<msub>
<mi>L</mi>
<mrow>
<mi>s</mi>
<mi>&sigma;</mi>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>L</mi>
<mrow>
<mi>d</mi>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mrow>
<mi>r</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mrow>
<mi>r</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mn>2</mn>
<mfrac>
<msub>
<mi>L</mi>
<mrow>
<mi>s</mi>
<mi>&sigma;</mi>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>L</mi>
<mrow>
<mi>q</mi>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mrow>
<mi>r</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mrow>
<mi>r</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mn>2</mn>
<mfrac>
<msub>
<mi>L</mi>
<mrow>
<mi>s</mi>
<mi>&sigma;</mi>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>L</mi>
<mrow>
<mi>q</mi>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mrow>
<mi>r</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>&CenterDot;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>(</mo>
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>u</mi>
<mrow>
<mi>s</mi>
<mi>&alpha;</mi>
</mrow>
<mo>&prime;</mo>
</msubsup>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mrow>
<mi>r</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>u</mi>
<mrow>
<mi>s</mi>
<mi>&beta;</mi>
</mrow>
</msub>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mrow>
<mi>r</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msubsup>
<mi>u</mi>
<mrow>
<mi>s</mi>
<mi>&alpha;</mi>
</mrow>
<mo>&prime;</mo>
</msubsup>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mrow>
<mi>r</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>u</mi>
<mrow>
<mi>s</mi>
<mi>&beta;</mi>
</mrow>
</msub>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mrow>
<mi>r</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>u</mi>
<mrow>
<mi>s</mi>
<mi>x</mi>
</mrow>
<mo>&prime;</mo>
</msubsup>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mrow>
<mi>r</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>u</mi>
<mrow>
<mi>s</mi>
<mi>y</mi>
</mrow>
</msub>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mrow>
<mi>r</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msubsup>
<mi>u</mi>
<mrow>
<mi>s</mi>
<mi>x</mi>
</mrow>
<mo>&prime;</mo>
</msubsup>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mrow>
<mi>r</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>u</mi>
<mrow>
<mi>s</mi>
<mi>y</mi>
</mrow>
</msub>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mrow>
<mi>r</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>+</mo>
<msqrt>
<mn>2</mn>
</msqrt>
<msub>
<mover>
<mi>u</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>o</mi>
<mn>2</mn>
</mrow>
</msub>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mrow>
<mi>r</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mrow>
<mi>r</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mrow>
<mi>r</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mrow>
<mi>r</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>s</mi>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>i</mi>
<mrow>
<mi>d</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>&omega;</mi>
<mrow>
<mi>r</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<msub>
<mi>&psi;</mi>
<mrow>
<mi>q</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>s</mi>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>i</mi>
<mrow>
<mi>q</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>&omega;</mi>
<mrow>
<mi>r</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<msub>
<mi>&psi;</mi>
<mrow>
<mi>d</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>s</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<mn>2</mn>
<msub>
<mi>R</mi>
<mrow>
<mi>s</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>)</mo>
<msub>
<mi>i</mi>
<mrow>
<mi>d</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>&omega;</mi>
<mrow>
<mi>r</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<msub>
<mi>&psi;</mi>
<mrow>
<mi>q</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>s</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<mn>2</mn>
<msub>
<mi>R</mi>
<mrow>
<mi>s</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>)</mo>
<msub>
<mi>i</mi>
<mrow>
<mi>q</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>&omega;</mi>
<mrow>
<mi>r</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<msub>
<mi>&psi;</mi>
<mrow>
<mi>d</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>;</mo>
</mrow>
步骤S44:将ψd1(k+1)ψq1(k+1)、ψd2(k+1)ψq2(k+1)传输至一磁链幅值计算模块,按照如下方式输出两台电机磁链幅值|ψs1(k+1)|、|ψs2(k+1)|:
<mrow>
<mo>|</mo>
<msub>
<mi>&psi;</mi>
<mrow>
<mi>s</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>|</mo>
<mo>=</mo>
<msqrt>
<mrow>
<msubsup>
<mi>&psi;</mi>
<mrow>
<mi>d</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>&psi;</mi>
<mrow>
<mi>q</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
</msqrt>
<mo>,</mo>
<mo>|</mo>
<msub>
<mi>&psi;</mi>
<mrow>
<mi>s</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>|</mo>
<mo>=</mo>
<msqrt>
<mrow>
<msubsup>
<mi>&psi;</mi>
<mrow>
<mi>d</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>&psi;</mi>
<mrow>
<mi>q</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
</msqrt>
<mo>;</mo>
</mrow>
步骤S44:把ψd1(k+1)、ψq1(k+1)、ψd2(k+1)、ψq2(k+1)传输至一电磁转矩预测模块,按照如下方式输出两台电机电磁转矩预测值Te1(k+1)、Te2(k+1):
<mrow>
<msub>
<mi>T</mi>
<mrow>
<mi>e</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<mfrac>
<mrow>
<msqrt>
<mn>3</mn>
</msqrt>
<msub>
<mi>&psi;</mi>
<mrow>
<mi>f</mi>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>&psi;</mi>
<mrow>
<mi>q</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mrow>
<msub>
<mi>L</mi>
<mrow>
<mi>d</mi>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>L</mi>
<mrow>
<mi>d</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>L</mi>
<mrow>
<mi>q</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
<msub>
<mi>&psi;</mi>
<mrow>
<mi>d</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<msub>
<mi>&psi;</mi>
<mrow>
<mi>q</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>L</mi>
<mrow>
<mi>d</mi>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>L</mi>
<mrow>
<mi>q</mi>
<mn>1</mn>
</mrow>
</msub>
</mrow>
</mfrac>
</mrow>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>T</mi>
<mrow>
<mi>e</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>p</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<mfrac>
<mrow>
<msqrt>
<mn>3</mn>
</msqrt>
<msub>
<mi>&psi;</mi>
<mrow>
<mi>f</mi>
<mn>2</mn>
</mrow>
</msub>
<msub>
<mi>&psi;</mi>
<mrow>
<mi>q</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mrow>
<msub>
<mi>L</mi>
<mrow>
<mi>d</mi>
<mn>2</mn>
</mrow>
</msub>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>L</mi>
<mrow>
<mi>d</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>L</mi>
<mrow>
<mi>q</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>)</mo>
<msub>
<mi>&psi;</mi>
<mrow>
<mi>d</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<msub>
<mi>&psi;</mi>
<mrow>
<mi>q</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>L</mi>
<mrow>
<mi>d</mi>
<mn>2</mn>
</mrow>
</msub>
<msub>
<mi>L</mi>
<mrow>
<mi>q</mi>
<mn>2</mn>
</mrow>
</msub>
</mrow>
</mfrac>
</mrow>
<mo>)</mo>
</mrow>
<mo>.</mo>
</mrow>
6.根据权利要求5所述的一种双电机串联缺相容错预测型直接转矩控制方法,其特征在于,在所述步骤S5中,所述评估函数值cost:
<mrow>
<mi>C</mi>
<mi>o</mi>
<mi>s</mi>
<mi>t</mi>
<mo>=</mo>
<msub>
<mi>k</mi>
<mn>1</mn>
</msub>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>T</mi>
<mrow>
<mi>e</mi>
<mn>1</mn>
</mrow>
<mo>*</mo>
</msubsup>
<mo>-</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>e</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>k</mi>
<mn>2</mn>
</msub>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>T</mi>
<mrow>
<mi>e</mi>
<mn>2</mn>
</mrow>
<mo>*</mo>
</msubsup>
<mo>-</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>e</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>k</mi>
<mn>3</mn>
</msub>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>&psi;</mi>
<mrow>
<mi>s</mi>
<mn>1</mn>
</mrow>
<mo>*</mo>
</msubsup>
<mo>-</mo>
<msub>
<mi>&psi;</mi>
<mrow>
<mi>s</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>k</mi>
<mn>4</mn>
</msub>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>&psi;</mi>
<mrow>
<mi>s</mi>
<mn>2</mn>
</mrow>
<mo>*</mo>
</msubsup>
<mo>-</mo>
<msub>
<mi>&psi;</mi>
<mrow>
<mi>s</mi>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
其中,k1、k2、k3、k4为损耗函数系数。
7.根据权利要求1或6所述的一种双电机串联缺相容错预测型直接转矩控制方法,其特征在于,当两台电机的控制量为电磁转矩时,所述给定电机电磁转矩由***给定;当两台电机的控制量为转速时,所述给定电机电磁转矩由两台电机速度控制器分别输出转矩给定;当两台电机的控制量为转子位置角时,所述给定电机电磁转矩由两台电机位置控制器输出转矩给定。
8.根据权利要求1所述的一种双电机串联缺相容错预测型直接转矩控制方法,其特征在于,在所述步骤S6中,把Sb~Sf值增加1,再从所述步骤S4开始计算;当Sb~Sf的取值完全后,找出对应cost值最小的Sb~Sf。
9.根据权利要求1、5或8所述的一种双电机串联缺相容错预测型直接转矩控制方法,其特征在于,Sb~Sf取值为Sb~Sf=00000~11111。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710522578.7A CN107276475B (zh) | 2017-06-30 | 2017-06-30 | 一种双电机串联缺相容错预测型直接转矩控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710522578.7A CN107276475B (zh) | 2017-06-30 | 2017-06-30 | 一种双电机串联缺相容错预测型直接转矩控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107276475A true CN107276475A (zh) | 2017-10-20 |
CN107276475B CN107276475B (zh) | 2020-04-10 |
Family
ID=60070258
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710522578.7A Active CN107276475B (zh) | 2017-06-30 | 2017-06-30 | 一种双电机串联缺相容错预测型直接转矩控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107276475B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111130413A (zh) * | 2020-01-03 | 2020-05-08 | 湖南大学 | 多相储能永磁电机定子磁链预测控制方法、***及介质 |
CN113687226A (zh) * | 2021-08-20 | 2021-11-23 | 中车大连电力牵引研发中心有限公司 | 一种城轨车控牵引逆变器异步电机缺相故障的识别方法 |
WO2022088440A1 (zh) * | 2020-10-29 | 2022-05-05 | 浙江大学 | 一种双电机转矩同步***模型预测电流控制方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101931352A (zh) * | 2010-07-14 | 2010-12-29 | 中国人民解放军海军航空工程学院 | 一种单逆变器驱动的双y移30°六相永磁同步电动机双电机串联***及控制方法 |
CN105743175A (zh) * | 2016-04-07 | 2016-07-06 | 东南大学 | 一种新型集成充电机功能的电动车驱动*** |
CN105897114A (zh) * | 2016-04-28 | 2016-08-24 | 南京航空航天大学 | 双定子双凸极电机故障容错驱动***及控制方法 |
CN106487308A (zh) * | 2016-11-04 | 2017-03-08 | 福州大学 | 串联电机驱动***输入缺一相容错型直接转矩控制方法 |
EP2903155B1 (fr) * | 2014-02-04 | 2017-07-05 | Thales | Dispositif de commande de deux moteurs polyphasés |
-
2017
- 2017-06-30 CN CN201710522578.7A patent/CN107276475B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101931352A (zh) * | 2010-07-14 | 2010-12-29 | 中国人民解放军海军航空工程学院 | 一种单逆变器驱动的双y移30°六相永磁同步电动机双电机串联***及控制方法 |
EP2903155B1 (fr) * | 2014-02-04 | 2017-07-05 | Thales | Dispositif de commande de deux moteurs polyphasés |
CN105743175A (zh) * | 2016-04-07 | 2016-07-06 | 东南大学 | 一种新型集成充电机功能的电动车驱动*** |
CN105897114A (zh) * | 2016-04-28 | 2016-08-24 | 南京航空航天大学 | 双定子双凸极电机故障容错驱动***及控制方法 |
CN106487308A (zh) * | 2016-11-04 | 2017-03-08 | 福州大学 | 串联电机驱动***输入缺一相容错型直接转矩控制方法 |
Non-Patent Citations (1)
Title |
---|
刘陵顺 等: "对称六相和三相PMSM串联***反电动势谐波效应补偿控制", 《大连理工大学学报》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111130413A (zh) * | 2020-01-03 | 2020-05-08 | 湖南大学 | 多相储能永磁电机定子磁链预测控制方法、***及介质 |
CN111130413B (zh) * | 2020-01-03 | 2021-06-04 | 湖南中聚福能源科技有限公司 | 多相储能永磁电机定子磁链预测控制方法、***及介质 |
WO2022088440A1 (zh) * | 2020-10-29 | 2022-05-05 | 浙江大学 | 一种双电机转矩同步***模型预测电流控制方法 |
CN113687226A (zh) * | 2021-08-20 | 2021-11-23 | 中车大连电力牵引研发中心有限公司 | 一种城轨车控牵引逆变器异步电机缺相故障的识别方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107276475B (zh) | 2020-04-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106487308B (zh) | 串联电机驱动***输入缺一相容错型直接转矩控制方法 | |
CN107196571A (zh) | 一种双电机串联预测型直接转矩控制方法 | |
CN107222146B (zh) | 双三相永磁同步电机高负载能力的直接转矩控制方法 | |
WO2022160828A1 (zh) | 电动汽车驱动与充电一体化电路及其转矩消除控制方法 | |
CN104184380B (zh) | 偏置60度六相永磁同步电机缺一相容错型转矩控制方法 | |
CN104270063B (zh) | 六相永磁同步电机缺任意两相容错型直接转矩控制方法 | |
CN110297182B (zh) | 一种模拟开绕组永磁同步电机的电力电子负载*** | |
CN105490604B (zh) | 一种三相四开关感应电动机变频调速***的预测控制方法 | |
CN106712601A (zh) | 基于铜耗最小的轴向磁场磁通切换容错电机容错控制方法 | |
CN107645259A (zh) | 一种电动汽车用驱动电机转矩闭环控制方法 | |
CN105048910A (zh) | 基于无传感器的高压级联变频器矢量控制***及其控制方法 | |
CN106059428A (zh) | 三相四开关逆变器驱动永磁同步电机的模型预测控制方法 | |
CN112117941A (zh) | 一种基于模型预测电流控制的开绕组永磁同步电机容错控制方法 | |
CN110299882B (zh) | 混合供电型开绕组永磁同步电机三矢量模型预测控制方法 | |
CN107994832A (zh) | 一种多通道永磁同步电机混合直接转矩控制方法 | |
CN107276475A (zh) | 一种双电机串联缺相容错预测型直接转矩控制方法 | |
CN107154716A (zh) | 双定子周向错位角电励磁直流电机及其控制方法 | |
CN107453664A (zh) | 一种基于模型预测的三相四开关逆变器永磁同步电机***磁链控制方法 | |
CN103560722A (zh) | 一种永磁直线同步电机控制装置及方法 | |
CN113131816A (zh) | 混合转子双定子同步电机最大转矩电流比控制***及方法 | |
CN107134966A (zh) | 一种基于有限状态预测的开绕组永磁电机电流控制方法 | |
CN104617827A (zh) | 一种电动汽车用轴向磁场磁通切换永磁电机容错控制方法 | |
CN108054973B (zh) | 降低悬浮电流的单绕组定子永磁型磁通切换电机驱动方法 | |
CN110034719A (zh) | 一种六相单绕组无轴承磁通切换电机缺一相容错控制方法 | |
CN104852658B (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 |