CN103647489B - 一种混合励磁同步电机效率最优控制方法 - Google Patents

Abstract

本发明公开了一种混合励磁同步电机效率最优控制方法，在满足电机输出转矩与转速情况下，控制d轴、q轴和励磁绕组电流，使混合励磁同步电机铁耗和铜耗最小，达到混合励磁同步电机效率最优，提高了电机效率。电机运行于低速区时，根据负载大小，通过d轴、q轴和励磁绕组电流协调控制，使铁耗与铜耗之和最小。电机运行于高速区时，利用d轴电流与励磁绕组电流共同弱磁，通过d轴、q轴和励磁绕组电流协调控制，使电机铜耗和铁耗最小，达到混合励磁同步电机效率最优。混合励磁同步电机效率最优控制方法，减小了电机损耗，提高了能源使用效率，达到节约电能效果。

Description

一种混合励磁同步电机效率最优控制方法

技术领域

本发明属于电气传动技术领域，涉及一种电流分配策略，特别是涉及一种混合励磁同步电机控制方法。

背景技术

混合励磁同步电机是在永磁同步与电励磁同步电机的基础上发展起来的一种宽调速电机，其主要目的是为了解决永磁同步电机气隙磁场难以调节的问题。混合励磁同步电机具有两种励磁源，一种是永磁体，另一种是电励磁，永磁体产生的磁势为主磁势，励磁绕组产生的磁势为辅磁势。这种电机结合了永磁同步与电励磁同步电机的优点，两种励磁源在电机气隙中相互作用产生主磁通，当电励磁线圈通入正向的励磁电流时，产生正向电磁转矩而增大了电机转矩；反之，当电励磁线圈通入反向励磁电流时，则产生反向磁场削弱气隙磁场达到弱磁升速的目的，从而拓宽了电机的调速范围。

目前，对于混合励磁电机控制方法及驱动***研究还较少，相关文献不是很多，基本可以分为两种：

（1）id=0的控制方法

id=0控制是混合励磁同步电机矢量控制中最为简单且使用最为广泛的一种控制算法，计算量小、控制方便；缺点是没有考虑磁阻转矩，定子电流没有达到最优状态，效率没有最优化。

（2）铜耗最小的控制方法

铜耗最小控制方法是在输出所要求的转矩下，保持铜耗最小。相对于id=0控制，铜耗最小控制提高了混合励磁同步电机的效率，但是没有考虑电机的铁耗，而混合励磁同步电机的铁耗又是损耗中较大的一部分，所以在这种控制方法中，电机的定子电流还是没有达到最优状态，电机效率仍然没有最优化。

发明内容

技术问题：本发明针对现有技术之不足，在分析现有混合励磁同步电机控制方法的基础上，提出了一种混合励磁同步电机效率最优控制方法。

技术方案：本发明的混合励磁同步电机效率最优控制方法，包括以下步骤：

（1）从电机主电路采集相电流i_a、i_b和励磁电流i_f，对电机进行准确初始位置检测，从电机编码器上采集信号，送入控制器进行处理，得出转速n和转子位置角θ；

（2）将采集的相电流i_a、i_b经信号调理和A/D转换，然后进行帕克变换，得到两相旋转坐标系下的d轴电流i_d和q轴电流i_q；

（3）将实测转速n与给定转速n^*比较后得到转速偏差Δn，将转速偏差Δn输入速度调节器经比例积分运算后得到转矩参考值将转矩参考值实测转速n和给定转速n^*输入电流分配器，判断实际转速是否小于弱磁基速，如是，则电机运行于低速区，进入步骤4），否则，电机运行于高速区，进入步骤5）；

（4）根据下列方程组求解计算d轴电流参考值i_dref、q轴电流参考值i_qref和励磁电流参考值i_fref：

$\left\{\begin{array}{c}{i}_{\mathrm{dref}}={\mathrm{ai}}_{\mathrm{fref}}+b\\ i_{\mathrm{qref}}=\sqrt{{\mathrm{ci}}_{\mathrm{fref}}^2+{\mathrm{di}}_{\mathrm{fref}}+e}\\ T_e^{*}=\frac{3}{2}p[{\mathrm{\ψ}}_f+({\mathrm{ai}}_{\mathrm{fref}}+b)(L_d-L_q)+M_{\mathrm{sf}}{i}_{\mathrm{fref}}]\sqrt{{\mathrm{ci}}_{\mathrm{fref}}^2+{\mathrm{di}}_{\mathrm{fref}}+e}\end{array}\right.$

式中，系数a、b、c、d、e别由下列式子求得：

$a=\frac{[\frac{{2R}_f(L_d-L_1)}{M_{\mathrm{sf}}}+]{2c}_{\mathrm{Fe}}{f}^{\mathrm{\β}}{M}_{\mathrm{sf}}({2L}_d-L_q)}{[{2c}_{\mathrm{Fe}}{f}^{\mathrm{\β}}{L}_d(L_q-{2L}_d)-{3R}_s]},$ $b=\frac{{2c}_{\mathrm{Fe}}{f}^{\mathrm{\β}}({2L}_d-L_q){\mathrm{\ψ}}_f}{[{2c}_{\mathrm{Fe}}{f}^{\mathrm{\β}}{L}_d(L_q-{2L}_d)-{3R}_s]}$

$c=\frac{2}{M_{\mathrm{sf}}({3R}_s+{2c}_{\mathrm{Fe}}{f}^{\mathrm{\β}}{L}_q^2)}[a(L_d-L_q)({2R}_f+{2c}_{\mathrm{Fe}}{f}^{\mathrm{\β}}{\mathrm{aM}}_{\mathrm{sf}}{L}_d+{2c}_{\mathrm{Fe}}{f}^{\mathrm{\β}}{M}_{\mathrm{sf}}^2),$

$+M_{\mathrm{sf}}({2R}_f+{2c}_{\mathrm{Fe}}{f}^{\mathrm{\β}}{\mathrm{aM}}_{\mathrm{sf}}{L}_d+{2c}_{\mathrm{Fe}}{f}^{\mathrm{\β}}{M}_{\mathrm{sf}}^2)]$

$d=\frac{1}{M_{\mathrm{sf}}({3R}_s+{2c}_{\mathrm{Fe}}{f}^{\mathrm{\β}}{L}_q^2)}\left[{\mathrm{\ψ}}_f({2R}_f+{2c}_{\mathrm{Fe}}{f}^{\mathrm{\β}}{\mathrm{aM}}_{\mathrm{sf}}{L}_d+{2c}_{\mathrm{Fe}}{f}^{\mathrm{\β}}{M}_{\mathrm{sf}}^2)\right]$

$+b(L_d-L_q)({2R}_f+{2c}_{\mathrm{Fe}}{f}^{\mathrm{\β}}{\mathrm{aM}}_{\mathrm{sf}}{L}_d+{2c}_{\mathrm{Fe}}{f}^{\mathrm{\β}}{M}_{\mathrm{sf}}^2),$

$+a(L_d-L_q)({2c}_{\mathrm{Fe}}{f}^{\mathrm{\β}}{\mathrm{bM}}_{\mathrm{sf}}{L}_d+{2c}_{\mathrm{Fe}}{f}^{\mathrm{\β}}{M}_{\mathrm{sf}}{\mathrm{\ψ}}_f)$

$+M_{\mathrm{sf}}({2c}_{\mathrm{Fe}}{f}^{\mathrm{\β}}{\mathrm{bM}}_{\mathrm{sf}}{L}_d+{2c}_{\mathrm{Fe}}{f}^{\mathrm{\β}}{M}_{\mathrm{sf}}{\mathrm{\ψ}}_f)]$

$e=\frac{[{\mathrm{\ψ}}_f({2c}_{\mathrm{Fe}}{f}^{\mathrm{\β}}{\mathrm{bM}}_{\mathrm{sf}}{L}_d+{2c}_{\mathrm{Fe}}{f}^{\mathrm{\β}}{M}_{\mathrm{sf}}{\mathrm{\ψ}}_f)+b(L_d-L_q)({2c}_{\mathrm{Fe}}{f}^{\mathrm{\β}}{\mathrm{bM}}_{\mathrm{sf}}{L}_d+{2c}_{\mathrm{Fe}}{f}^{\mathrm{\β}}{M}_{\mathrm{sf}}{\mathrm{\ψ}}_f)]}{M_{\mathrm{sf}}({3R}_s+{2c}_{\mathrm{Fe}}{f}^{\mathrm{\β}}{L}_q^2)}$

i_dref为d轴电流参考值，i_qref为q轴电流参考值，i_fref为励磁绕组电流参考值；L_d为d轴电感，L_q为q轴电感，M_sf为电枢绕组与励磁绕组之间的互感，ψ_f为永磁体磁链，p为电机极对数，R_s为定子相绕组电阻，R_f为励磁绕组电阻，c_Fe为铁耗系数，β为电机铁耗修正参数，取值1.5～2之间，f为电机运行频率；

（5）根据下列方程组求解计算d轴电流参考值i_dref、q轴电流参考值i_qref与励磁电流参考值i_fref：

$\left\{\begin{array}{c}{i}_{\mathrm{dref}}=\frac{{2R}_f{L}_d{\mathrm{\ψ}}_f}{{3R}_s{M}_{\mathrm{sf}}^2+{2R}_f{L}_d^2}(\frac{n_{\mathrm{Bdec}}}{n_r}-1)\\ i_{\mathrm{qref}}={\mathrm{kT}}_e^{*}\\ i_{\mathrm{fref}}=\frac{{3R}_s{M}_{\mathrm{sf}}{\mathrm{\ψ}}_f}{{3R}_s{M}_{\mathrm{sf}}^2+{2R}_f{L}_d^2}(\frac{n_{\mathrm{Bdec}}}{n_r}-1)\end{array}\right.$

其中，n_r为电机速度，n_Bdec为电机弱磁基速，k为速度调节器比例系数，L_d为d轴电感，M_sf为电枢绕组与励磁绕组之间的互感，ψ_f为永磁体磁链，R_s为定子相绕组电阻，R_f为励磁绕组电阻；

（6）将d轴电流参考值i_dref和q轴电流i_qref分别与步骤（2）中的d轴电流i_d和q轴电流i_q比较后得到d轴电流偏差Δi_d和q轴电流偏差Δi_q，将d轴电流偏差Δi_d输入d轴电流调节器进行比例积分运算，得到d轴电压u_d，将q轴电流偏差Δi_q输入q轴电流调节器进行比例积分运算，得到q轴电压u_q，然后对d轴电压u_d和q轴电压u_q进行旋转正交-静止两相变换后，得到静止两相坐标系下α轴电压u_α和β轴电压u_β，将α轴电压u_α和β轴电压u_β输入脉冲宽度调制模块，运算输出6路脉冲宽度调制信号，驱动主功率变换器；

同时将步骤（1）中采集的励磁电流i_f，经信号调理与A/D转换后和励磁电流参考值i_fref一起送入直流励磁脉宽调制模块，运算输出4路脉冲宽度调制信号来驱动励磁功率变换器。

本发明方法的一种优选方案中，步骤6）中的脉冲宽度调制模块为空间矢量脉冲宽度调制模块。

有益效果：现有混合励磁同步电机控制方法主要有id=0控制和铜耗最小控制。其中id=0控制的思想是：在转矩给定情况下，保持d轴电流恒等于1，转矩与定子电流呈线性关系，控制简单，但是对于磁阻率较大的混合励磁同步电机，id=0控制没有充分利用磁阻转矩，定子电流没有达到最优，所以电机效率没有达到最优。铜耗最小控制的思想是：在转矩给定情况下，保持混合励磁同步电机铜耗最小，充分利用了磁阻转矩，优化了电流，提高了电机效率，但是铜耗最小控制没有考虑电机的铁耗，电机效率还没有达到最优情况。所以本发明提出了一种基于铜耗和铁耗最小的混合励磁同步电机效率最优控制方法。本发明通过步骤4）和步骤5）的效率最优控制方法，使混合励磁同步电机无论运行在低速区还是高速区，都使得电机铜耗和铁耗之和最小，达到效率最优。本发明相对现有控制方法具有以下优点：

（1）效率最优控制提高了电机效率；

（2）效率最优控制提高了能源使用率；

（3）效率最优控制达到节约电能效果。

附图说明

图1是本发明方法的逻辑流程框图；

图2是本发明方法的***框图；

图3是实现本发明方法的结构框图；

图4是电流分配结构框图。

具体实施方式

图3为实现本发明混合励磁同步电机效率最优控制方法的***框图，该控制***由交流电源、整流器、母线电容、DSP控制器、主功率变换器、辅功率变换器、传感器、混合励磁同步电机、光电编码器等组成。

交流电源给整个***供电，经过整流器整流后，滤波、稳压，送给主、辅功率变换器，霍尔电压传感器采集母线电压，调理后送入控制器。主、辅功率变换器的输出端接混合励磁同步电机，霍尔电流互感器采集相电流和励磁电流，调理后送入控制器，编码器信号采集转速与转子位置信号，处理后送入控制器计算转子位置角与转速。控制器输出10路PWM信号分别驱动主、励磁功率变换器。

本发明的混合励磁同步电机效率最优控制方法，图3所示，具体包括以下步骤：

（1）三个霍尔电流传感器分别从电机主电路采集相电流i_a、i_b和励磁电流i_f，将采集到的信号经电压跟随、滤波、偏置及过压保护等信号调理后送入控制器，对电机进行准确初始位置检测，从电机编码器上采集信号，处理送入控制器计算得出转速n和转子位置角θ；

（2）将送入控制器的相电流i_a、i_b进行A/D转换，经过三相坐标系到两相旋转坐标系的帕克变换得到两相旋转坐标系下的d轴电流i_d和q轴电流i_q；

（3）将编码器实测转速n与给定转速n^*比较后得到转速偏差Δn，转速偏差Δn进入速度调节器后得到转矩参考值将转矩参考值实测转速n和给定转速n^*送入电流分配器，判断实际转速n是否小于弱磁基速n_Bdec，如是，则电机运行于低速区，进入步骤4），否则，电机运行于高速区，进入步骤5）；

（4）下面分析混合励磁同步电机的效率最优控制原理，根据矢量控制原理，在d-q坐标系中，得出混合励磁同步电机的数学模型。

磁链方程：

$\left[\begin{array}{c}{\mathrm{\ψ}}_d\\ {\mathrm{\ψ}}_q\\ {\mathrm{\ψ}}_f\end{array}\right]=\left[\begin{array}{ccc}{L}_d& 0& M_{\mathrm{sf}}\\ 0& L_q& 0\\ 3/2M_{\mathrm{sf}}& 0& L_f\end{array}\right]\left[\begin{array}{c}{i}_d\\ i_q\\ i_f\end{array}\right]+\left[\begin{array}{c}{\mathrm{\ψ}}_f\\ 0\\ {\mathrm{\ψ}}_{\mathrm{fm}}\end{array}\right]---\left(1\right)$

电压方程：

$\left\{\begin{array}{c}{u}_d=R_s{i}_d+\frac{{\mathrm{d\ψ}}_d}{\mathrm{dt}}-{\mathrm{\ω}}_e{\mathrm{\ψ}}_q\\ u_q=R_s{i}_q+\frac{{\mathrm{d\ψ}}_q}{\mathrm{dt}}+{\mathrm{\ω}}_e{\mathrm{\ψ}}_d\\ u_f=R_f{i}_f+\frac{{\mathrm{d\ψ}}_f}{\mathrm{dt}}\end{array}\right.---\left(2\right)$

转矩方程：

$T_e=\frac{3}{2}{\mathrm{pi}}_q[{\mathrm{\ψ}}_f+i_d(L_d-L_q)+M_{\mathrm{sf}}{i}_f]---\left(3\right)$

其中，i_d、i_q分别为d轴与q轴电流，i_f为励磁绕组电流；L_d、L_q分别为d轴与q轴电感，M_sf为电枢与励磁绕组之间的互感；ω_e为电角速度；ψ_f为永磁体磁链，p为电机极对数，u_d、u_q分别为d轴与q轴的电压，u_f为励磁绕组电压；R_s为电枢绕组电阻，R_f为励磁绕组电阻；ψ_d、ψ_q、ψ_f分别d轴、q轴与励磁绕组磁链；ψ_f为永磁体磁链幅值，ψ_fm为永磁体穿过励磁绕组的磁链。

混合励磁同步电机总的损耗表达式为：

$P_{\mathrm{Loss}}=P_{\mathrm{Cu}}+P_{\mathrm{Fe}}=\frac{3}{2}{R}_s(i_d^2+i_q^2)+R_f{i}_f^2+c_{\mathrm{Fe}}{\mathrm{\ω}}_e^{\mathrm{\β}}[{(L_d{i}_d+{\mathrm{\ψ}}_f+M_{\mathrm{sf}}{i}_f)}^2+{\left(L_q{i}_q\right)}^2]---\left(4\right)$

其中， $P_{\mathrm{Cu}}=\frac{3}{2}{R}_s(i_d^2+i_q^2)+R_f{i}_f^2$ 铜耗表达式；

$P_{\mathrm{Fe}}=c_{\mathrm{Fe}}{\mathrm{\ω}}_e^{\mathrm{\β}}[{(L_d{i}_d+{\mathrm{\ψ}}_f+M_{\mathrm{sf}}{i}_f)}^2+{\left(L_q{i}_q\right)}^2]$ 铁耗表达式。

根据式（3）与式（4）建立如下拉格朗日函数，λ为拉格朗日乘子。

$L(i_d,i_q,i_f,\mathrm{\λ})=\frac{3}{2}{R}_s(i_d^2+i_q^2)+R_f{i}_f^2+c_{\mathrm{Fe}}{\mathrm{\ω}}_e^{\mathrm{\β}}[{(L_d{i}_d+{\mathrm{\ψ}}_f+M_{\mathrm{sf}}{i}_f)}^2+{\left(L_q{i}_q\right)}^2]---\left(5\right)$

$+\mathrm{\λ}\left\{\frac{3}{2}{\mathrm{pi}}_q\right[{\mathrm{\ψ}}_f+(L_d-L_q)i_d+M_{\mathrm{sf}}{i}_f]-T_e\}$

上式分别对i_d、i_q、i_f求导，得

$\left\{\begin{array}{c}\frac{\∂L}{{\∂i}_d}={3R}_s{i}_d+{2c}_{\mathrm{Fe}}{\mathrm{\ω}}_e^{\mathrm{\β}}{L}_d(L_d{i}_d+{\mathrm{\ψ}}_f+M_{\mathrm{sf}}{i}_f)-\frac{3}{2}\mathrm{\λp}(L_d-L_q)i_q\\ \frac{\∂L}{\∂i_q}={3R}_s{i}_q+{2c}_{\mathrm{Fe}}{\mathrm{\ω}}_e^{\mathrm{\β}}{L}_q^2{i}_q+\frac{3}{2}\mathrm{\λp}[{\mathrm{\ψ}}_f+(L_d-L_q)i_d+M_{\mathrm{sf}}{i}_f]\\ \frac{\∂L}{{\∂i}_f}={2R}_f{i}_f+{2c}_{\mathrm{Fe}}{\mathrm{\ω}}_e^{\mathrm{\β}}{M}_{\mathrm{sf}}(L_d{i}_d+{\mathrm{\ψ}}_f+M_{\mathrm{sf}}{i}_f)+\frac{3}{2}\mathrm{\λp}{M}_{\mathrm{sf}}{i}_q\end{array}\right.---\left(6\right)$

令得式（7），根据式（7）计算d轴电流参考值i_dref、q轴电流参考值i_qref与励磁电流参考值i_fref。

$\left\{\begin{array}{c}{i}_{\mathrm{dref}}={\mathrm{ai}}_{\mathrm{fref}}+b\\ i_{\mathrm{qref}}=\sqrt{{\mathrm{ci}}_{\mathrm{fref}}^2+{\mathrm{di}}_{\mathrm{fref}}+e}\\ T_e^{*}=\frac{3}{2}p[{\mathrm{\ψ}}_f+({\mathrm{ai}}_{\mathrm{fref}}+b)(L_d-L_q)+M_{\mathrm{sf}}{i}_{\mathrm{fref}}]\sqrt{{\mathrm{ci}}_{\mathrm{fref}}^2+{\mathrm{di}}_{\mathrm{fref}}+e}\end{array}\right.---\left(7\right)$

式中，系数a、b、c、d、e别由下列式子求得：

$a=\frac{[\frac{{2R}_f(L_d-L_q)}{M_{\mathrm{sf}}}+{2c}_{\mathrm{Fe}}{f}^{\mathrm{\β}}{M}_{\mathrm{sf}}({2L}_d-L_q)]}{\left[{2c}_{\mathrm{Fe}}{f}^{\mathrm{\β}}{L}_d(L_q-{2L}_d){-3R}_s\right]},$ $b=\frac{{2c}_{\mathrm{Fe}}{f}^{\mathrm{\β}}({2L}_d-L_q){\mathrm{\ψ}}_f}{[{2c}_{\mathrm{Fe}}{f}^{\mathrm{\β}}{L}_d(L_q-{2L}_d)-{3R}_s]}$

$c=\frac{1}{M_{\mathrm{sf}}({3R}_s+{2c}_{\mathrm{Fe}}{f}^{\mathrm{\β}}{L}_q^2)}[a(L_d-L_q)({2R}_f+{2c}_{\mathrm{Fe}}{f}^{\mathrm{\β}}{\mathrm{aM}}_{\mathrm{sf}}{L}_d+{2c}_{\mathrm{Fe}}{f}^{\mathrm{\β}}{M}_{\mathrm{sf}}^2),$

$+M_{\mathrm{sf}}({2R}_f+{2c}_{\mathrm{Fe}}{f}^{\mathrm{\β}}{\mathrm{aM}}_{\mathrm{sf}}{L}_d+{2c}_{\mathrm{Fe}}{f}^{\mathrm{\β}}{M}_{\mathrm{sf}}^2)]$

$d=\frac{1}{M_{\mathrm{sf}}({3R}_s+{2c}_{\mathrm{Fe}}{f}^{\mathrm{\β}}{L}_q^2)}[{\mathrm{\ψ}}_f({2R}_f+{2c}_{\mathrm{Fe}}{f}^{\mathrm{\β}}{\mathrm{aM}}_{\mathrm{sf}}{L}_d+{2c}_{\mathrm{Fe}}{f}^{\mathrm{\β}}{M}_{\mathrm{sf}}^2)$

$+b(L_d-L_q)({2R}_f+{2c}_{\mathrm{Fe}}{f}^{\mathrm{\β}}{\mathrm{aM}}_{\mathrm{sf}}{L}_d+{2c}_{\mathrm{Fe}}{f}^{\mathrm{\β}}{M}_{\mathrm{sf}}^2),$

$+a(L_d-L_q)({2c}_{\mathrm{Fe}}{f}^{\mathrm{\β}}{\mathrm{bM}}_{\mathrm{sf}}{L}_d+{2c}_{\mathrm{Fe}}{f}^{\mathrm{\β}}{M}_{\mathrm{sf}}{\mathrm{\ψ}}_f)$

$+M_{\mathrm{sf}}({2c}_{\mathrm{Fe}}{f}^{\mathrm{\β}}{\mathrm{bM}}_{\mathrm{sf}}{L}_d+{2c}_{\mathrm{Fe}}{f}^{\mathrm{\β}}{M}_{\mathrm{sf}}{\mathrm{\ψ}}_f)]$

$e=\frac{[{\mathrm{\ψ}}_f({2c}_{\mathrm{Fe}}{f}^{\mathrm{\β}}{\mathrm{bM}}_{\mathrm{sf}}{L}_d+{2c}_{\mathrm{Fe}}{f}^{\mathrm{\β}}{M}_{\mathrm{sf}}{\mathrm{\ψ}}_f)+b(L_d-L_q)({2c}_{\mathrm{Fe}}{f}^{\mathrm{\β}}{\mathrm{bM}}_{\mathrm{sf}}{L}_d+{2c}_{\mathrm{Fe}}{f}^{\mathrm{\β}}{M}_{\mathrm{sf}}{\mathrm{\ψ}}_f)]}{M_{\mathrm{sf}}({3R}_s+{2c}_{\mathrm{Fe}}{f}^{\mathrm{\β}}{L}_q^2)}$

i_dref为d轴电流参考值，i_qref为q轴电流参考值，i_fref为励磁绕组电流参考值；L_d为d轴电感，L_q为q轴电感，M_sf为电枢绕组与励磁绕组之间的互感，ψ_f为永磁体磁链，p为电机极对数，R_s为定子相绕组电阻，R_f为励磁绕组电阻，c_Fe为铁耗系数，β为电机铁耗修正参数，取值1.5～2之间，f为电机运行频率；

（5）建立如下拉格朗日函数，λ为拉格朗日乘子。

$L(i_d,i_f,\mathrm{\λ})=\frac{3}{2}{R}_s{i}_d^2+R_f{i}_f^2+c_{\mathrm{Fe}}{\mathrm{\ω}}_e^{\mathrm{\β}}{(L_d{i}_d+{\mathrm{\ψ}}_f+M_{\mathrm{sf}}{i}_f)}^2+\mathrm{\λ}[L_d{i}_d+M_{\mathrm{sf}}{i}_f-{\mathrm{\ψ}}_f(\frac{n_{\mathrm{Bdec}}}{n_r}-1)]---\left(8\right)$

上式分别对i_d、i_f求导，得

$\left\{\begin{array}{c}\frac{\∂L}{\∂i_d}={3R}_s{i}_d+{2c}_{\mathrm{Fe}}{\mathrm{\ω}}_e^{\mathrm{\β}}{L}_d(L_d{i}_d+{\mathrm{\ψ}}_f+M_{\mathrm{sf}}{i}_f)+\mathrm{\λ}{L}_d\\ \frac{\∂L}{\∂i_f}={2R}_f{i}_f+{2c}_{\mathrm{Fe}}{\mathrm{\ω}}_e^{\mathrm{\β}}{M}_{\mathrm{sf}}(L_d{i}_d+{\mathrm{\ψ}}_f+M_{\mathrm{sf}}{i}_f)+\mathrm{\λ}{M}_{\mathrm{sf}}\\ L_d{i}_d+M_{\mathrm{sf}}{i}_f={\mathrm{\ψ}}_f(\frac{n_{\mathrm{Bdec}}}{n_r}-1)\end{array}\right.---\left(9\right)$

令得式（10），根据式（10）计算d轴电流参考值i_dref、q轴电流参考值i_qref与励磁电流参考值i_fref。

$\left\{\begin{array}{c}{i}_{\mathrm{dref}}=\frac{2R_f{L}_d{\mathrm{\ψ}}_f}{{3R}_s{M}_{\mathrm{sf}}^2+{2R}_f{L}_d^2}(\frac{n_{\mathrm{Bdec}}}{n_r}-1)\\ i_{\mathrm{qref}}={\mathrm{kT}}_e^{*}\\ i_{\mathrm{fref}}=\frac{{3R}_s{M}_{\mathrm{sf}}{\mathrm{\ψ}}_f}{{3R}_s{M}_{\mathrm{sf}}^2+{2R}_f{L}_d^2}(\frac{n_{\mathrm{Bdec}}}{n_r}-1)\end{array}\right.---\left(10\right)$

其中，n_r为电机速度，n_Bdec为电机弱磁基速，k为速度调节器比例系数；

（6）将d轴电流参考值i_dref和q轴电流i_qref分别与步骤（2）中的d轴电流i_d和q轴电流i_q比较后得d轴电流偏差Δi_d和q轴电流偏差Δi_q，将Δi_d与Δi_q分别送入d轴电流调节器与q轴电流调节器，得到d轴电压u_d和q轴电压u_q，进行旋转正交-静止两相变换后得到静止两相坐标系下的电压信号u_α与u_β，送入空间矢量脉冲宽度调制模块后输出6路脉冲宽度调制信号，驱动主功率变换器；同时将步骤（1）中采集的励磁电流i_f，经信号调理与A/D转换后和励磁电流参考值i_fref一起送入直流励磁脉宽调制模块，运算输出4路脉冲宽度调制信号来驱动励磁功率变换器。

Claims (2)

1.一种混合励磁同步电机效率最优控制方法，其特征在于，该方法包括以下步骤：

(1)从电机主电路采集相电流i_a、i_b和励磁电流i_f，对电机进行准确初始位置检测，从电机编码器上采集信号，送入控制器进行处理，得出转速n和转子位置角θ；

(2)将采集的相电流i_a、i_b经信号调理和A/D转换，然后进行帕克变换，得到两相旋转坐标系下的d轴电流i_d和q轴电流i_q；

(3)将实测转速n与给定转速n^*比较后得到转速偏差Δn，将所述转速偏差Δn输入速度调节器经比例积分运算后得到转矩参考值将转矩参考值实测转速n和给定转速n^*输入电流分配器，判断实际转速是否小于弱磁基速，如是，则电机运行于低速区，进入步骤(4)，否则，电机运行于高速区，进入步骤(5)；

(4)根据下列方程组求解计算d轴电流参考值i_dref、q轴电流参考值i_qref和励磁电流参考值i_fref：

$\left\{\begin{array}{c}{i}_{dref}=a{i}_{fref}+b\\ i_{qref}=\sqrt{{\mathrm{ci}}_{fref}^2+{\mathrm{di}}_{fref}+e}\\ T_e^{*}=\frac{3}{2}p\[{\mathrm{\ψ}}_f+({\mathrm{ai}}_{fref}+b)(L_d-L_q)+M_{sf}{i}_{fref}\]\sqrt{{\mathrm{ci}}_{fref}^2+{\mathrm{di}}_{fref}+e}\end{array}\right.$

式中，系数a、b、c、d、e别由下列式子求得：

$a=\frac{\[\frac{2R_f(L_d-L_q)}{M_{sf}}+2c_{Fe}{f}^{\mathrm{\β}}{M}_{sf}(2L_d-L_q)\]}{\[2c_{Fe}{f}^{\mathrm{\β}}{L}_d(L_q-2L_d)-3R_s\]},b=\frac{2c_{Fe}{f}^{\mathrm{\β}}(2L_d-L_q){\mathrm{\ψ}}_f}{\[2c_{Fe}{f}^{\mathrm{\β}}{L}_d(L_q-2L_d)-3R_s\]}$

$\begin{array}{c}c=\frac{1}{M_{sf}(3R_s+2c_{Fe}{f}^{\mathrm{\β}}{L}_q^2)}\[a(L_d-L_q)(2R_f+2c_{Fe}{f}^{\mathrm{\β}}{\mathrm{aM}}_{sf}{L}_d+2c_{Fe}{f}^{\mathrm{\β}}{M}_{sf}^2)\\ +M_{sf}(2R_f+2c_{Fe}{f}^{\mathrm{\β}}{\mathrm{aM}}_{sf}{L}_d+2c_{Fe}{f}^{\mathrm{\β}}{M}_{sf}^2)\end{array}$

$\begin{array}{c}d=\frac{1}{M_{sf}(3R_s+2c_{Fe}{f}^{\mathrm{\β}}{L}_q^2)}{\mathrm{\ψ}}_f(2R_f+2c_{Fe}{f}^{\mathrm{\β}}{\mathrm{aM}}_{sf}{L}_d+2c_{Fe}{f}^{\mathrm{\β}}{M}_{sf}^2)\\ +b(L_d-L_q)(2R_f+2c_{Fe}{f}^{\mathrm{\β}}{\mathrm{aM}}_{sf}{L}_d+2c_{Fe}{f}^{\mathrm{\β}}{M}_{sf}^2)\\ +a(L_d-L_q)(2c_{Fe}{f}^{\mathrm{\β}}{\mathrm{bM}}_{sf}{L}_d+2c_{Fe}{f}^{\mathrm{\β}}{M}_{sf}{\mathrm{\ψ}}_f)\\ +M_{sf}(2c_{Fe}{f}^{\mathrm{\β}}{\mathrm{bM}}_{sf}{L}_d+2c_{Fe}{f}^{\mathrm{\β}}{M}_{sf}{\mathrm{\ψ}}_f)\end{array}$

$e=\frac{\[{\mathrm{\ψ}}_f(2c_{Fe}{f}^{\mathrm{\β}}{\mathrm{bM}}_{sf}{L}_d+2c_{Fe}{f}^{\mathrm{\β}}{M}_{sf}{\mathrm{\ψ}}_f)+b(L_d-L_q)(2c_{Fe}{f}^{\mathrm{\β}}{\mathrm{bM}}_{sf}{L}_d+2c_{Fe}{f}^{\mathrm{\β}}{M}_{sf}{\mathrm{\ψ}}_f)\]}{M_{sf}(3R_s+2c_{Fe}{f}^{\mathrm{\β}}{L}_q^2)}$

i_dref为d轴电流参考值，i_qref为q轴电流参考值，i_fref为励磁绕组电流参考值；L_d为d轴电感，L_q为q轴电感，M_sf为电枢绕组与励磁绕组之间的互感，ψ_f为永磁体磁链，p为电机极对数，R_s为定子相绕组电阻，R_f为励磁绕组电阻，c_Fe为铁耗系数，β为电机铁耗修正参数，取值1.5～2之间，f为电机运行频率；

(5)根据下列方程组求解计算d轴电流参考值i_dref、q轴电流参考值i_qref与励磁电流参考值i_fref：

$\left\{\begin{array}{c}{i}_{dref}=\frac{2R_f{L}_d{\mathrm{\ψ}}_f}{3R_S{M}_{sf}^2+2R_f{L}_d^2}(\frac{n_{Bdec}}{n_r}-1)\\ i_{\mathrm{qr}ef}={\mathrm{kT}}_e^{*}\\ i_{fref}=\frac{2R_s{M}_{sf}{\mathrm{\ψ}}_f}{3R_S{M}_{sf}^2+2R_f{L}_d^2}(\frac{n_{Bdec}}{n_r}-1)\end{array}\right.$

其中，n_r为电机速度，n_Bdec为电机弱磁基速，k为速度调节器比例系数，L_d为d轴电感，M_sf为电枢绕组与励磁绕组之间的互感，ψ_f为永磁体磁链，R_s为定子相绕组电阻，R_f为励磁绕组电阻；

(6)将d轴电流参考值i_dref和q轴电流i_qref分别与所述步骤(2)中的d轴电流i_d和q轴电流i_q比较后得到d轴电流偏差Δi_d和q轴电流偏差Δi_q，将所述d轴电流偏差Δi_d输入d轴电流调节器进行比例积分运算，得到d轴电压u_d，将q轴电流偏差Δi_q输入q轴电流调节器进行比例积分运算，得到q轴电压u_q，然后对所述d轴电压u_d和q轴电压u_q进行旋转正交-静止两相变换后，得到静止两相坐标系下α轴电压u_α和β轴电压u_β，将所述α轴电压u_α和β轴电压u_β输入脉冲宽度调制模块，运算输出6路脉冲宽度调制信号，驱动主功率变换器；

同时将步骤(1)中采集的励磁电流i_f，经信号调理与A/D转换后和励磁电流参考值i_fref一起送入直流励磁脉宽调制模块，运算输出4路脉冲宽度调制信号来驱动励磁功率变换器。

2.根据权利要求1所述的混合励磁同步电机效率最优控制方法，所述步骤(6)中的脉冲宽度调制模块为空间矢量脉冲宽度调制模块。