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

郭佳强 刘 旭 李珊瑚

（1. 省部共建电工装备可靠性与智能化国家重点实验室（河北工业大学） 天津 300130 2. 河北工业大学河北省电磁场与电器可靠性重点实验室 天津 300130）

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/4V FRM的拓扑结构。VFRM的电枢绕组和励磁绕组并联。在VFRM的绕组中，定子电流含有交流和直流分量。定子电流的交流分量充当电枢电流，而直流分量充当励磁电流。为了给零序电流提供回路，采用公共直流母线的双逆变器结构。三相电流的表达式为

式中，ia、ib、ic为三相电流；iac为相电流的交流分量幅值；idc为相电流的直流分量（在后文中用i0表示）；ωe为电角速度。

图1 使用零序电流励磁的6/4 VFRM的拓扑结构 Fig.1 Topology of zero-sequence current excited 6/4 VFRM

VFRM在dq0坐标系的电压方程如式（2）所示，其中，3次谐波分量为静止坐标系中的2次谐波在dq0平面上的映射[31]。

式中，ud、uq和u0分别为dq0轴电压分量；id、iq、i0分别为dq0轴电流分量；Rs为绕组的电阻；Ls为定子电感的恒定分量；Lδ为定子电感的变化分量；θe为电角度。

由式（2）可知，在正弦电流条件下，反电动势中含有谐波，而基波电流环只能输出正弦电压，不能补偿这些谐波，在实际控制中会造成电流畸变。畸变的电流不利于弱磁区中的电流控制，还会降低控制系统的稳定性，增加电机的损耗。因此，在本文中，通过谐波电流环对基波电压进行补偿，从而抑制电流的畸变[24]。通过谐波电流环，一方面对反电动势中的谐波进行补偿，抑制了谐波电流；另一方面，可以将定子电压分离为基波分量和谐波分量，有利于弱磁区中参考电流的计算。在使用谐波电流环对电压进行补偿后，在弱磁控制中推导参考电流时，可以只使用电压方程的基频分量。这时，将电压方程式（2）中的3次分量忽略，得到VFRM的简化电压方程为

VFRM的瞬时转矩为

式中，Nr为转子极数；β为超前电流角。

弱磁控制中参考电流的计算仅涉及平均转矩。因此，仅使用电磁转矩方程式的平均转矩表达式。平均转矩是式（4）中的恒定分量，即

2 弱磁区控制方法

2.1 电流和电压约束

基于逆变器或电机的温度限制，可以确定最大定子电流imax[32]。参考电流应满足

最大电压umax由直流母线电压和脉冲宽度调制（Pulse Width Modulation, PWM）策略决定[32]。在VFRM中，反电动势中的谐波会降低直流母线电压的利用率[4]。相电压的峰值大于基波幅值。因此，当相电压的峰值达到Udc时，相应的基波分量幅值被确定为最大电压，从而避免逆变器饱和。考虑零序电压后，电压约束关系为

将式（3）代入式（7）并忽略dq轴上的电阻性电压降，可以将电压约束重写为

2.2 最优电流工作点的计算

在弱磁区中，VFRM有3个电流分量需要控制。为了计算最优的电流参考值，采用拉格朗日乘子法计算id、iq和i0。为求转矩的极大值，根据转矩公式，设目标函数为

使用电压和电流限制方程作为约束条件。得到拉格朗日函数FL为

式中，λ1和λ2为拉格朗日乘子。

对拉格朗日函数中变量id、iq和i0以及λ1、λ2求偏微分可得

在方程组式（11）中，消去λ1和λ2，则得到以下方程组。该方程组的根即考虑电阻压降后弱磁区内的最优id、iq、i0值。

id和iq可以用i0表示。利用式（12）第二个等式可以得到id的表达式为

利用式（12）第三个等式和式（13）可以得到iq为

由式（12）第一个等式，设函数N0()f i为

函数N0()f i中id和iq可由式（13）和式（14）表示。对于i0，可以通过牛顿拉夫逊算法，对N0()f i= 0迭代求解得出。函数N0()f i的导数为N0()f i′ 。每次迭代后的根为i0(k+1)。当相邻两次迭代结果的差满足精度要求时，将完成迭代。迭代过程为

设置初始值为1，当所需精度为0.01A时，可以用多次迭代来完成计算。求得i0的给定值后，通过式（13）和式（14）可以求得id和iq的给定值。

dq平面中的电压极限圆方程为

电压极限圆的中心坐标为

忽略电阻压降时，将Rs设置为0。电压极限圆心的坐标为

图2为忽略Rs时，电压和电流极限圆随转速升高的变化情况。从式（19）可以看出，电压极限圆的圆心始终在d轴上。随着速度的增加，电压极限圆缩小。为了保持在弱磁区输出最大功率，根据式（12），i0将随速度增加。并且，电压极限圆的圆心向d轴的负方向移动。此外，根据式（6），随着i0的增加，电流极限圆的半径减小。电压极限圆与电流极限圆的交点即为当前工作点。

图2 忽略电阻压降的电流轨迹 Fig.2 Current trajectory with zero-sequence resistive voltage drop neglected

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

图3 考虑电阻压降的电流轨迹 Fig.3 Current trajectory with zero-sequence resistive voltage drop considered

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

图4 恒转矩区边界时工作点情况 Fig.4 Operating point at the boundary of constant torque region

2.3 最大电压计算

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

图5 相电压波形及其谐波分析 Fig.5 Phase voltage waveforms and corresponding harmonic analysis

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

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

稳态下，特定次数的谐波在相应次数下的dq分量是恒定值，但是其他次数谐波在本次数dq坐标系中的映射是交流分量。例如，基波在2次dq坐标系下表现为3次谐波，其变换矩阵为

因此，可以通过截止频率低于基频的低通滤波器来提取对应次数的dq分量。2次谐波电流环如图6所示。通过2次dq变换和低通滤波，可以获得2次谐波的dq轴分量。将调制后的二阶αβ 轴电压叠加在基频上，可以达到抑制2次谐波电流的目的。图6中，带“*”的变量代表该变量是给定值。本文后面“*”的含义也相同。

图6 2次谐波电流环框图 Fig.6 Scheme of second harmonic current loop

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

从图5a可知，式（2）计算出的相电压波形与真实值比较接近。因此可以推测相电压峰值为

式中，ud和uq采用式（2）进行计算。

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

图7 弱磁控制系统框图 Fig.7 Drive scheme with flux-weakening control

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

图8 电感随电流变化的情况 Fig.8 Inductance variation according to currents

3 实验分析

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

表1 VFRM参数 Tab.1 Parameters of VFRM

图9 实验平台 Fig.9 Picture of the experimental setup

当转速为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为考虑和忽略电阻压降时的转矩-转速曲线。由于忽略了电阻压降，所以电机在弱磁区中无法跟随参考电流，这限制了电机在高速时功率的输出。

图10 稳态电流波形（1 500r/min） Fig.10 Current waveforms at steady state (speed at 1 500r/min)

图11 VFRM在弱磁区内dq0轴电流随转速的变化 Fig.11 Currents of the VFRM varies with speed in flux-weakening region

图12 转矩-转速曲线 Fig.12 Torque-speed curves

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

图13 功率-转速曲线 Fig.13 Power-speed curves

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

图14 效率-转速曲线 Fig.14 Efficiency-speed curves

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

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

图15 功率因数-转速曲线 Fig.15 Power factor-speed curves

图16 电机升速（0～2 000r/min）时的实验结果 Fig.16 The experimental results while the speed increases from 0r/min to 2 000r/min

4 结论

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