计及电枢和零序回路电阻压降的可变磁通磁阻电机弱磁控制

　　摘要 在可变磁通磁阻电机的弱磁控制中，忽略电枢和零序回路上的电阻压降会影响弱磁区内参考电流的计算，这会导致电流环饱和，电机电流不能跟随给定值。在这种情况下，电机在弱磁区的输出功率会受到限制。为了解决此问题，提出一种计及电枢和零序回路电阻压降的弱磁控制方法。首先，该方法基于电压方程以及电压和电流的约束，利用拉格朗日极值法求取弱磁区内考虑电阻压降后的最优电流参考值的公式。然后，利用牛顿拉夫逊迭代法根据当前电机的状态计算出电流参考值。考虑电阻压降后，弱磁区内参考电流的计算准确度得到提高，有利于可变磁通磁阻电机在弱磁区内的稳定运行。最后，通过实验结果验证了该方法的有效性。

　　本文源自郭佳强; 刘旭; 李珊瑚， 电工技术学报 发表时间：2021-07-05

　　关键词：可变磁通磁阻电机 弱磁控制 电阻压降 参考电流

　　0 引言

　　可变磁通磁阻电机(Variable Flux Reluctance Machine, VFRM)是一种定子励磁的无刷同步电机[1]，其转子上没有绕组或永磁体。VFRM 的散热性能良好、调速范围宽、性能稳定、成本低[2-3]，具有良好的应用前景。在 VFRM中使用集成励磁和电枢绕组，绕组中注入含有交流和直流分量的定子电流。定子电流的交流分量作为电枢电流，直流分量作为励磁电流[4]。根据零序分量的定义[5-11]，定子电流的直流分量是零序电流。为了给零序电流提供回路，在 VFRM 的控制系统中使用共直流母线的双逆变器结构[12-13]。

　　VFRM 的反电动势幅值与转速成正比，因此在额定负载下，当电机的速度超过基速时，逆变器的输出电压将无法满足需求[14]。为了拓宽调速范围，需要通过弱磁控制将 d 轴电流调节为负值从而削弱反电动势[15]。弱磁控制通常分为基于反馈或基于前馈的方法[16-25]。在基于反馈的方法中，通常通过电流环的输出电压与最大电压之间的偏差来调制 d 轴电流。在前馈方法中，通常通过等转矩曲线，电流或电压极限圆的交点来计算弱磁区中的工作点。在基于前馈的方法中，为了简化电流参考值的计算，通常忽略电阻压降。文献[26]分析了忽略电阻压降对弱磁控制的影响：当忽略绕组电阻压降时，参考电流将高于电机的实际需求，使得电流环饱和，从而导致逆变器饱和，这不仅增加了电流谐波，还将使电机失控甚至损坏。考虑电阻压降后计算出的电流参考值可以使电机在弱磁区域内沿着最佳的电流轨迹稳定运行。文献[27-28]提出基于插值的参考电流计算方法，该方法考虑了弱磁区域的 dq 轴上的电阻压降。文献[29]分析了绕组和逆变器中的电阻压降，并用来补偿 q 轴电压，从而使电机沿着最佳电流轨迹运行。文献[30]提出一种基于牛顿拉夫逊算法的参考电流在线计算的方法。然而，上述研究针对的都是永磁同步电机中 dq 轴(电枢绕组)上的电阻压降。与永磁同步电机相比，VFRM 除了电枢绕组上的电阻压降外，还有零序回路中的电阻压降，如果要准确计算弱磁区中的参考电流，需要同时考虑电枢和零序回路上的电阻压降。

　　为了解决上述问题，本文对 VFRM 在弱磁区内考虑电阻压降的电流参考值计算进行分析。基于电压和电流的约束方程，采用拉格朗日乘数法计算出弱磁区中 VFRM 的最优参考电流方程，并用牛顿拉夫逊算法进行迭代求解。首先介绍 VFRM 的拓扑和数学模型。然后分析了考虑电阻压降的最优工作点轨迹的计算方法。最后将该方法在一台 6/4VFRM 上实现，通过实验结果对该方法进行验证。

　　1 VFRM 的模型

　　图 1 为使用零序电流励磁的 6/4VFRM 的拓扑结构。VFRM 的电枢绕组和励磁绕组并联。在 VFRM 的绕组中，定子电流含有交流和直流分量。定子电流的交流分量充当电枢电流，而直流分量充当励磁电流。为了给零序电流提供回路，采用公共直流母线的双逆变器结构。三相电流的表达式为  a a c e d c b a c e d c c a c e d c c o s 2 π c o s 3 2 π c o s 3 i i t i i i t i i i t i                        (1)式中，ia、ib、ic 为三相电流;iac 为相电流的交流分量幅值;idc 为相电流的直流分量(在后文中用 i0 表示);e 为电角速度。

　　完整的 VFRM 电压方程如式(2)所示，其中含有 3 次变化分量[4]，这些 3 次分量是相 2 次谐波在 dq 平面上的映射[31]，有 d δ δ s e e δ d d q δ δ q s q e s e 0 0 δ 0 s d c o s 3 s in 3 2 2 d d s in 3 c o s 3 0 2 2 d d 0 2 d L L i L L t u i L L i u R i L t u i L i L t                                                                   δ e s δ e d e s δ e δ e δ q 0 s in 3 c o s 3 0 c o s 3 s in 3 0 0 0 L L L i L L L L i i                       (2)式中，ud、uq 和 u0 分别为 dq0 轴电压分量;id、iq、 i0 分别为 dq0 轴电流分量;Rs 为绕组的电阻;Ls 为定子电感的恒定分量;L 为定子电感的变化分量;e 为电角度。

　　由式(2)可知，在正弦电流条件下，反电动势中含有谐波，而基波电流环只能输出正弦电压，不能补偿这些谐波，在实际控制中会造成电流畸变。畸变的电流不利于弱磁区中的电流控制，还会降低控制系统的稳定性，增加电机的损耗。因此，在本文中，通过谐波电流环对基波电压进行补偿，从而抑制电流的畸变[24]。通过谐波电流环，一方面对反电动势中的谐波进行补偿，抑制了谐波电流;另一方面，可以将定子电压分离为基波分量和谐波分量，有利于弱磁区中参考电流的计算。在使用谐波电流环对电压进行补偿后，在弱磁控制中推导参考电流时，可以只使用电压方程的基频分量。这时，将电压方程式(2)中的 3 次分量忽略，得到 VFRM 的简化电压方程为 d 0 d s d e s q s δ q q s q e s d e δ 0 s δ d 0 0 s 0 s d d d d d d d d 2 d d i i u R i L i L L t t i u R i L i L i L t L i i u R i L t t              (3) VFRM 的瞬时转矩为    r 2 2 e δ 0 q δ d q e 3 1 s in 3 2 2 4 N T L i i L i i          (4)式中，Nr 为转子极数; 为超前电流角。

　　弱磁控制中参考电流的计算仅涉及平均转矩。因此，仅使用电磁转矩方程式的平均转矩表达式。平均转矩是式(4)中的恒定分量，即 e r δ q 0 3 2 T N L i i  (5)

　　2 弱磁区控制方法

　　2.1 电流和电压约束

　　基于逆变器或电机的温度限制，可以确定 imax (最大定子电流)[33]。参考电流应满足 2 2 d q 2 2 0 m a x 2 2 i i   i i ≤ (6)最大电压 umax 由直流母线电压和脉冲宽度调制(Pulse Width Modulation, PWM)策略决定[33]。在 VFRM 中，反电动势中的谐波会降低直流母线电压的利用率[4]。相电压的峰值大于基波幅值。因此，当相电压的峰值达到 Udc 时，相应的基波分量幅值被确定为最大电压，从而避免逆变器饱和。考虑零序电压后，电压约束关系为  2 2 2 d q m a x 0 u u u u   ≤ (7)将式(3)代入式(7)并忽略 dq 轴上的电阻性电压降，可以将电压约束重写为    2 2 2 2 m a x s 0 s s q s d δ 0 e u R i u L i L i L i         ≤ (8)

　　2.2 最优电流工作点的计算

　　在弱磁区中，VFRM 有 3 个电流分量需要控制。为了计算最优的电流参考值，采用拉格朗日乘子法计算 id、iq 和 i0。为求转矩的极大值，根据转矩公式，设目标函数为 f i i i i L q 0 q 0  ,   (9)使用电压和电流限制方程作为约束条件。得到拉格朗日函数 FL 为  2 2 L q 0 1 2 s d e 2 d q 2 2 0 s s q 2 2 i i F i i R i L i    i i                  2 2 s d e s q e δ 0 m a x s 0 R i L i L i u R i      (10)式中，1 和2 为拉格朗日乘子。

　　对拉格朗日函数中变量 id、iq 和 i0 以及1、2 求偏微分可得函数 N 0 f i ( ) 中 id 和 iq 可由式(13)和式(14)表示。对于 i0，可以通过牛顿拉夫逊算法，对 N 0 f i ( ) = 0 迭代求解得出。函数 N 0 f i ( ) 的导数为 N 0 f i  ( ) 。每次迭代后的根为 i0(k+1)。当相邻两次迭代结果的差满足精度要求时，将完成迭代。迭代过程为   N 0 0 0 N 0 ( ) ( 1) ( ) ( ) f i k i k i k f i k    (16)设置初始值为 1，当所需精度为 0.01A 时，可以用多次迭代来完成计算。求得 i0 的给定值后，通过式(13)和式(14)可以求得 id 和 iq 的给定值。

　　dq 平面中的电压极限圆方程为下方程组。该方程组的根即考虑电阻压降后弱磁区内的最优 id、iq、i0 值。        2 2 s d e s q s q e s d e δ 0 2 m a 2 2 2 2 2 2 δ s 0 d δ s 0 δ s s 0 e s d s 0 s x s 2 2 s 0 d q 2 2 0 0 2 2 2 3 0 3 0 2 2 L L i i L L i L L R i L i R i L i L i u i i u i R i R R i i i i i                             (12) id 和 iq 可以用 i0 表示。利用式(12)第二个等式可以得到 id 的表达式为  2 2 2 δ s e 0 s d 2 2 2 2 2 δ 0 e s 0 e s 0 s s 2 2 2 3 L L i i i L i L i R i R u     (13)利用式(12)第三个等式和式(13)可以得到 iq 为 电压极限圆的中心坐标为 δ s 0 δ s 0 e 2 2 2 2 s s s s 2 2 e e , L R i L L i R R L L             (18)忽略电阻压降时，将 Rs 设置为 0。电压极限圆心的坐标为 δ 0 s , 0 L i L      (19)图 2 为忽略 Rs 时，电压和电流极限圆随转速升高的变化情况。从式(19)可以看出，电压极限圆的圆心始终在 d 轴上。随着速度的增加，电压极限圆缩小。为了保持在弱磁区输出最大功率，根据式(12)，i0 将随速度增加。并且，电压极限圆的圆心向 d 轴的负方向移动。此外，根据式(6)，随着 i0 的增加，电流极限圆的半径减小。电压极限圆与电流极限圆的交点即为当前工作点。

　　图 3 为考虑电阻压降时电压和电流的极限圆随速度的增加而变化的示意图。电机从零开始加速，负载增加到定子电流达到最大值。在此过程中，电机在 A 点处工作。当转速为基速时，电机运行到达恒定转矩区域的边界。如果转速继续上升，需要进行弱磁控制。在弱磁区，工作点将移动到 dq 电流平面的第二象限，其轨迹为 AB，如图 3 所示。与 A 和 B 所在电压圆对应的圆心分别为 O1 和 O3。

　　图 4 为当转速达到恒定转矩区域的边界时考虑和忽略电阻压降的工作点分析。图 4 中的电压圆 1 和 2 分别是忽略和考虑电阻压降的情况。A 点是实际工作点。如果忽略电阻压降，根据电压圆 1，电机尚在恒转矩区域，而根据电压圆 2，则达到了逆变器的最大输出电压。如果转速继续提高，在忽略电阻压降时，由于控制系统认为电机仍在恒转矩区，所以仍将根据恒定转矩区域的规则计算参考电流。由于电压和电流的限制，电机电流将无法跟踪该参考电流值，在这种情况下，电流调节器会饱和，从而导致电机失去控制。为了使电机稳定运行，工作点保持在电压极限圆内，必须减小负载。这样，输出功率将降低。相对应的，考虑电阻压降后，可以准确获取电机工作所需的电流参考值。

　　2.3 最大电压计算

　　根据前文分析可知，VFRM 的定子电压获取对判断是否进入弱磁区和参考电流的计算非常重要。然而，定子电压的峰值与基波分量的幅值不一致，这给电压约束计算带来了困难。因此本节中对最大电压进行分析。

　　根据电压方程式(2)可知，VFRM 的 dq 轴电压中含有 3 次谐波变化分量。这部分 3 次谐波变化分量是相电压的 2 次谐波在 dq 平面中的映射[27]。在转速 400r/min、iq 为 2A、id=0A 的条件下，根据式(2)计算出 dq 轴电压的波形后通过 Park 反变换，可以得到相电压的波形，如图 5 所示[34]。为了抑制相电流中的谐波，本文使用了指定谐波次数抑制法[31]。在指定谐波次数抑制法中，需要抑制的谐波电流都有对应的谐波电流环。

　　以 2 次谐波电流为例。为了抑制 2 次谐波电流，构建了 2 次谐波电流环，并将参考值设为零。2 次谐波的旋转坐标变换矩阵为

　　稳态下，特定次数的谐波在相应次数下的 dq 分量是恒定值，但是其他次数谐波在本次数 dq 坐标系中的映射是交流分量。例如，基波在 2 次 dq 坐标系下表现为 3 次谐波，其变换矩阵为 e e dq2 dq1 e e c o s (3 ) s in (3 ) 0 s in (3 ) c o s (3 ) 0 0 0 1           C (21)因此，可以通过截止频率低于基频的低通滤波器来提取对应次数的 dq 分量。2 次谐波电流环如图 6 所示。通过 2 次 dq 变换和低通滤波，可以获得 2 次谐波的 dq 轴分量。将调制后的二阶 轴电压叠加在基频上，可以达到抑制 2 次谐波电流的目的。图 6 中，带“*”的变量代表该变量是给定值。本文后面“*”的含义也相同。

　　通过谐波电流环，一方面对反电动势中的谐波进行补偿，抑制了谐波电流;另一方面，将相电压分为了基波分量和 2 次谐波分量。对于弱磁控制，需要基波的电压幅值，这样，基波相电压的 dq 轴分量被分离出，更容易求取定子电压的基波分量。

　　从图 5a 可知，公式计算出的相电压波形与真实值比较接近。因此可以推测相电压峰值为 2 2 s d q 0 u u u u    (22)式中，ud 和 uq 采用式(2)进行计算。

　　最大电压 umax=36V 情况下，根据式(22)计算出的定子电压在转速 0～2 000r/min 范围内峰值的最大值约为 62V。考虑 5%的电压裕量后，65V 的直流 母 线 电 压 可 以 使 电 机 的 最 大 基 波 定 子 电 压 为 36V。

　　整个弱磁控制系统框图如图 7 所示。参考电流由 MTPA 和弱磁模块计算。共直流母线的双逆变器用来为零序电流提供回路。电压补偿模块用于对反电动势中特定次数的谐波进行补偿。另外，通过使用有限元法生成的电感表来考虑不同工作条件下的电感变化[35]。如图 8 所示为电感随电流变化的有限元仿真结果。电流的变化范围是 id(1.2～0A)，iq (1～2A)，i0 为 1A 和 2A。可见，电感随电流变化的范围很小，可以认为电机没有进入饱和区。

　　3 实验分析

　　为了验证考虑电阻压降的弱磁控制方法，搭建了基于 6/4 VFRM 的双逆变器控制系统。表 1 是实验中使用电机的参数。实验平台如图 9 所示。实验中使用的控制器是 dSPACE;VFRM 由双逆变器驱动;IGBT 的开关频率设置为 10kHz;电机的输出转矩由奇石乐 4 502A 传感器测量;本实验使用的示波器是横河 DL850W;功率、效率和功率因数等参数的测量采用横河 WT3000E 功率分析仪。

　　当转速为 100r/min 时，增加负载转矩使定子电流达到最大值。然后，增加转速并记录电流、速度、转矩等变量。图 10 为在稳态下测得的 dq0 轴电流以及 A 相电流的波形。此时的转速为 1 500r/min，转矩为 0.32N·m。VFRM 这时处于弱磁区，id 为负值。图 11 为随着转速的增加，测得的 id、iq、i0 的变化。在 0～930r/min 阶段，电机处于恒转矩区，负载转矩约为 0.41N·m，其中，id=0，iq 和 i0 保持固定比例。在 930～2 000r/min 阶段，电机处于弱磁区，负载转矩随着速度的增加而减小。当电阻压降被忽略并且定子电压被低估时，参考电流将继续根据恒转矩区的规则进行计算。但是，由于电压和电流的限制，电机在转速提高后会无法跟踪参考电流。为了使 VFRM 可控，必须降低输出转矩。考虑电阻压降后， VFRM 可以沿着弱磁区的最佳电流轨迹运行。图 12 为考虑和忽略电阻压降时的转矩-转速曲线。由于忽略了电阻压降，所以电机在弱磁区中无法跟随参考电流，这限制了电机在高速时功率的输出。

　　如图 13 所示为实验测得的考虑和忽略电阻压降时的功率-转速曲线。采用文中控制后，可变磁通磁阻电机可以准确判断并进入弱磁区。因此，采用考虑电阻压降的弱磁控制后，样机可以在弱磁区内沿着最优电流轨迹运行，并输出最大功率。而忽略电阻压降时，对于本实验中使用的样机，由于控制系统不能准确判断其进入弱磁控制的时刻，因此在高速时的输出功率下降。

　　图 14 为实验测得的采用文中方法前后的电机效率。不考虑电阻压降时，电机在高速时定子电流幅值减小，因此铜耗也减小。在输入输出功率均减小的情况下，其效率反而比考虑电阻压降时的高。但是在其输出功率减小的情况下，效率高没有意义。

　　图 15 为实验测得的不同转速下功率因数变化情况。在恒转矩区中，有功功率中的铜耗较大，因此功率因数在低速区很高。随着转速的升高，无功功率上升，铜耗占比减小，功率因数也逐渐减小。在弱磁区中，考虑电阻压降控制下，由于样机能沿着最优电流轨迹运行，电流和电压的有效值均不变，因此视在功率不变。在这种情况下，功率因数与有功功率的趋势有关。电流有效值不变时，铜耗相对变化较小。因此功率因数在弱磁区中与输出功率的变化趋势一致，均是先增加后趋于平稳。而在未考虑电阻压降控制下，高速时电机的输入输出功率均减小。

　　图 16 为在 0～2 000r/min 的过程中转速、转矩和电流的变化情况。当给出转速指令时，为了使电机加速，控制系统会增加 iq 和 i0，以使电机的输出转矩迅速增加，达到升速的目的。同时，在速度增加到基速后，电机进入弱磁区。为了削弱反电动势，id 被调节为负值。在速度达到 2 000r/min 后，电机重新恢复稳定状态。

　　4 结论

　　本文提出考虑电阻压降的 VFRM的弱磁控制策略。基于电机的电压方程和电压、电流约束，使用拉格朗日乘数法和牛顿拉夫逊法来计算考虑电枢和零序回路电阻压降的最优参考电流。另外，分析了定子电压最大值的计算。最后在一台 6/4 VFRM 上实现了本文的控制策略并验证了其有效性。实验结果表明，如果忽略电枢和零序回路上的电阻压降，则弱磁区中电流参考值的计算会受到影响。在这种情况下，逆变器将会饱和，电机中的电流将不能跟随给定值。因此，忽略电阻压降会造成电机在弱磁区的输出功率下降。而考虑了电枢和零序回路电阻压降后，VFRM 可以在弱磁区沿着最优电流轨迹稳定运行并保持最大功率输出。