舰船科学技术  2020, Vol. 42 Issue (8): 140-144    DOI: 10.3404/j.issn.1672-7649.2020.08.026   PDF    
六相双Y30°绕组感应电机建模与控制技术研究
王贤明, 程晗, 何露, 吴建峰     
中国船舶集团公司 第七二二研究所,湖北 武汉 430205
摘要: 六相感应电机具有转矩脉动低、电机损耗小、电机极限容量大、能量密度高等一系列优点。本文基于三相感应电机αβ两相静止坐标系下的数学模型表达式,将六相双Y30°绕组感应电机空间磁场分布等效成2套三相电机绕组磁势在空间上的合成,推导出空间坐标6/2矢量变换表达式,并设计出六相双Y30°绕组感应电机矢量控制系统,采用速度外环、电流内环双环控制系统。通过理论分析及仿真,验证了本文电机数学模型、坐标变换、控制策略的有效性。
关键词: 感应电机     六相     转子磁链定向     矢量控制    
Research on modeling and control technology of six-phase induction motor with double Y-connected 3-phase symmetrical windings displaced in turn by 30°
WANG Xian-ming, CHENG Han, HE Lu, WU Jian-feng     
The 722 Research Institute of CSSC, Wuhan 430205, China
Abstract: The six-phase induction motor has a series of advantages such as low torque ripple, low motor loss, large motor limit capacity, and high energy density. Based on the mathematical model expression of αβ two-phase stationary coordinate system of three-phase induction motor, the spatial magnetic field distribution of six-phase double Y30° winding induction motor is equivalent to the spatial synthesis of two sets of three-phase motor windings. Coordinates 6/2 vector transformation expression, and designed a six-phase double Y30° winding induction motor vector control system are given. The double loop control system of speed outer loop and current inner loop is adopted. Through theoretical analysis and simulation, the effectiveness of the mathematical model, coordinate transformation and control strategy of the motor is verified.
Key words: induction motor     six-phase     rotor flux-oriented     vector control    
0 引 言

六相感应电机具有转矩脉动低、电机损耗小、电机极限容量大、能量密度高等一系列优点,因此成为国内外学者近年来研究的热点[18]。六相双Y30°绕组感应电机的定子绕组由空间相差30°电角度2套结构完全一样的三相绕组所组成,转子结构与三相感应电机笼型结构相同。在定子绕组缺相时可以降低负载功率继续运行,从而提升六相感应电机运行的可靠性,在舰船动力推进、航空航天、电动汽车等领域的研究及应用日益广泛。

本文基于三相感应电机矢量控制策略扩展到六相感应电机中,通过推导出空间坐标6/2矢量变换表达式,设计出六相双Y30°绕组感应电机矢量控制系统。通过理论分析及仿真,验证了本文电机数学模型、坐标变换、控制策略的有效性。

1 六相双Y30°绕组感应电机模型

为了更好地研究及分析六相双Y30°绕组感应电机的特性,实现六相双Y30°绕组电机高性能矢量控制,需要从空间物理结构及数学模型等方面展开相应的研究。

1.1 物理结构模型

六相双Y30°绕组感应电机的绕组物理结构如图1所示。在空间上具有2套三相互差1200度对称绕组,2套对称绕组之间相关30°的电角度。六相双Y30°绕组感应电机空间磁场分布可以等效成2套三相电机绕组磁势在空间上的合成。下面从数学模型的角度进一步分析六相绕组之间的相互关系,并推导出相应的坐标变换表达式。

图 1 六相双Y30°绕组感应电机绕组物理结构 Fig. 1 Winding physical structure of six-phase winding induction motor with double Y30°
1.2 三相感应电机数学模型

三相感应电机静止αβ坐标系磁链方程:

$\left\{ \begin{aligned} & {\psi _{s\alpha }} = {L_s}{i_{s\alpha }} + {L_m}{i_{r\alpha }}\text{,} \\ & {\psi _{s\beta }} = {L_s}{i_{s\beta }} + {L_m}{i_{r\beta }}\text{,} \\ & {\psi _{r\alpha }} = {L_r}{i_{r\alpha }} + {L_m}{i_{s\alpha }}\text{,} \\ & {\psi _{r\beta }} = {L_r}{i_{r\beta }} + {L_m}{i_{s\beta }}\text{;} \\ \end{aligned} \right.$ (1)

三相感应电机静止αβ坐标系电压方程:

$\left\{ \begin{aligned} & {u_{s\alpha }} = {R_s}{i_{s\alpha }} + p{\psi _{s\alpha }} = {R_s}{i_{s\alpha }} + {L_s}p{i_{s\alpha }} + {L_m}p{i_{r\alpha }} \text{,} \\ & {u_{s\beta }} = {R_s}{i_{s\beta }} + p{\psi _{s\beta }} = {R_s}{i_{s\beta }} + {L_s}p{i_{s\beta }} + {L_m}p{i_{r\beta }} \text{,}\\ & 0 = {R_r}{i_{r\alpha }} + {\omega _r}{\psi _{r\beta }} + p{\psi _{r\alpha }} \text{,} \\ & 0 = {R_r}{i_{r\beta }} - {\omega _r}{\psi _{r\alpha }} + p{\psi _{r\beta }}{\rm{ }} \text{。} \\ \end{aligned} \right.$ (2)

三相感应电机αβ两相静止坐标系下的数学模型表达式如下:

$\small \begin{array}{*{20}{l}} \begin{gathered} \left[ {\begin{array}{*{20}{l}} {p{i_{s\alpha }}} \\ {p{i_{s\beta }}} \\ {p{\psi _{r\alpha }}} \\ {p{\psi _{r\beta }}} \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} { - \displaystyle\frac{a}{b}}&{{\text{0}}}&{\displaystyle\frac{{{R_r}{L_m}}}{{{L_r}b}}}&{\displaystyle\frac{{{\omega _r}{L_m}}}{b}} \\ 0& {- \displaystyle\frac{a}{b}}&{ - \displaystyle\frac{{{\omega _r}{L_m}}}{b}}&{\displaystyle\frac{{{R_r}{L_m}}}{b}} \\ {\displaystyle\frac{{{R_r}{L_m}}}{{{L_r}}}}&{{\text{0}}}&{ - \displaystyle\frac{{{R_r}}}{{{L_r}}}}&{ - {\omega _r}} \\ 0&{\displaystyle\frac{{{R_r}{L_m}}}{{{L_r}}}}&{{\omega _r}}&{ - \displaystyle\frac{{{R_r}}}{{{L_r}}}} \end{array}} \right]\left[ {\begin{array}{*{20}{l}} {{i_{s\alpha }}} \\ {{i_{s\beta }}} \\ {{\psi _{r\alpha }}} \\ {{\psi _{r\beta }}} \end{array}} \right] + \hfill \\ \quad \quad \quad \quad \quad \left[ {\begin{array}{*{20}{l}} {\displaystyle\frac{{{L_r}}}{b}}&{{\text{0}}} \\ 0&{\displaystyle\frac{{{L_r}}}{b}} \\ 0&{{\text{0}}} \\ 0&{0} \end{array}} \right]\left[ {\begin{array}{*{20}{l}} {{u_{s\alpha }}} \\ {{u_{s\beta }}} \end{array}} \right]\text{。} \hfill \\[-10pt] \end{gathered} \end{array} $ (3)

式中:Rs为定子绕组电阻;Rr为转子绕组电阻;Ls为定子绕组自感;Lr为转子绕组自感;Lm为绕组激磁电感;ωr为电机电角速度;i为定子绕组α轴电流;isβ为定子绕组β轴电流;irα为转子绕组α轴电流;irβ为转子绕组β轴电流;usα为定子绕组α轴电压;usβ为定子绕组β轴电压;ψsα为定子α轴磁链;ψsβ为定子β轴磁链;ψrα为转子α轴磁链;ψrβ为转子β轴磁链。

1.3 六相双Y30°绕组感应电机数学模型

六相双Y30°绕组感应电机的电压方程为:

${u_{6s}} = {R_{6s}}{i_{6s}} + p{\psi _{6s}}\text{,}$ (4)

$\small\left[ \begin{array}{l} {u_A} \\ {u_B} \\ {u_C} \\ {u_D} \\ {u_E} \\ {u_F} \\ \end{array} \right] = R\left[ \begin{array}{l} 1\;\;\;0\;\;\;0\;\;\;0\;\;\;0\;\;\;0 \\ 0\;\;\;1\;\;\;0\;\;\;0\;\;\;0\;\;\;0 \\ 0\;\;\;0\;\;\;1\;\;\;0\;\;\;0\;\;\;0 \\ 0\;\;\;0\;\;\;0\;\;\;1\;\;\;0\;\;\;0 \\ 0\;\;\;0\;\;\;0\;\;\;0\;\;\;1\;\;\;0 \\ 0\;\;\;0\;\;\;0\;\;\;0\;\;\;0\;\;\;1 \\ \end{array} \right]\left[ \begin{array}{l} {i_A} \\ {i_B} \\ {i_C} \\ {i_D} \\ {i_E} \\ {i_F} \\ \end{array} \right] + p\left[ \begin{array}{l} {\psi _A} \\ {\psi _B} \\ {\psi _C} \\ {\psi _D} \\ {\psi _E} \\ {\psi _F} \\ \end{array} \right]\text{。}$ (5)

图2为双Y30°绕组六相静止到两相静止6 s/2 s变换的坐标关系。为了方便起见,令双Y30°绕组的A轴与两相的α轴重合。假设磁势波形正弦分布,只计其基波分量,当二者的旋转磁势完全等效时,合成磁势沿相同轴向的分量一定相同。也就是说,双Y30°绕组与两相绕组沿αβ轴瞬时磁势的投影相等,可以推出:

图 2 六相静止与两相静止坐标系示意图 Fig. 2 Schematic diagram of six-phase stationary and two-phase stationary coordinate system
$\small\left\{ \begin{aligned} & {F_\alpha } = {F_A} + {F_B}\cos ({\text π} /6) + {F_C}\cos (2{\text π} /3) +\\ & \quad \quad\; {F_D}\cos (5{\text π}/6) + {F_E}\cos (4{\text π} /3) \text{,} \\ & {F_\beta } = {F_B}\sin ({\text π} /6) + {F_C}\sin (2{\text π} /3) + \\ & {\rm{ }}{F_D}\sin (5{\text π} /6) + {F_E}\sin (4{\text π} /3) - {F_F} \text{。} \\ \end{aligned} \right.$ (6)

写成矩阵形式为:

$ \left[\!\! \begin{array}{l} {F_\alpha }\\ {F_\beta } \end{array} \!\! \right] = \left[ \!\!{\begin{array}{*{20}{l}} 1&{\sqrt 3 /2}&{ - {\rm{1/2}}}&{ - \sqrt 3 /2}&{ - {\rm{1/2}}}&{\rm{0}}\\ 0&{1/2}&{\sqrt 3 /2}&{{\rm{1/2}}}&{ - \sqrt 3 /2}&{ - {\rm{1}}} \end{array}} \!\!\right]{F_{6s}} \text{。} $ (7)

虽然零序分量没有物理意义,但为了将上式中的矩阵用单位正交阵表示,需要补充定义4个零序分量:

$ \begin{aligned} & {C_{6s/2s}} = \frac{1}{{\sqrt 3 }} \times \\ & \left[ {\begin{array}{*{20}{l}} 1&{\sqrt 3 /2}&{ - {\text{1/2}}}&{ - \sqrt 3 /2}&{ - {\text{1/2}}}&{{\text{0}}} \\ 0&{1/2}&{\sqrt 3 /2}&{{\text{1/2}}}&{ - \sqrt 3 /2}&{ - {\text{1}}} \\ 1&{ - \sqrt 3 /2}&{ - {\text{1/2}}}&{\sqrt 3 /2}&{ - {\text{1/2}}}&{{\text{0}}} \\ 0&{1/2}&{ - \sqrt 3 /2}&{{\text{1/2}}}&{\sqrt 3 /2}&{ - {\text{1}}} \\ 1&{0}&{1}&{0}&{1}&{0} \\ 0&{1}&{0}&{1}&{0}&{1} \end{array}} \right] \text{,} \end{aligned} $ (8)

为了验证上面坐标变换矩阵的正交性,得出 ${C_{6s/2s}}$ 的转置矩阵如下:

$ \begin{aligned} & {C_{2s/6s}} = C_{6s/2s}^T = \frac{1}{{\sqrt 3 }} \times \\ & \left[ {\begin{array}{*{20}{l}} 1&{0}&{1}&{0}&{1}&{0} \\ {\sqrt 3 /2}&{1/2}&{ - \sqrt 3 /2}&{{\text{1/2}}}&{0}&{1} \\ { - {\text{1/2}}}&{\sqrt 3 /2}&{ - {\text{1/2}}}&{ - \sqrt 3 /2}&{1}&{0} \\ { - \sqrt 3 /2}&{1/2}&{\sqrt 3 /2}&{{\text{1/2}}}&{0}&{1} \\ { - {\text{1/2}}}&{ - \sqrt 3 /2}&{ - {\text{1/2}}}&{\sqrt 3 /2}&{1}&{0} \\ 0&{ - 1}&{0}&{ - 1}&{0}&{1} \end{array}} \right] \text{,} \end{aligned} $ (9)

将两矩阵 ${C_{6s/2s}}$ ${C_{2s/6s}}$ 相乘得到单位矩阵,从而说明 ${C_{6s/2s}}$ 具有正交性。

${C_{6s/2s}} \times {C_{2s/6s}} = \left[ \begin{array}{l} 1\;\;\;0\;\;\;0\;\;\;0\;\;\;0\;\;\;0 \\ 0\;\;\;1\;\;\;0\;\;\;0\;\;\;0\;\;\;0 \\ 0\;\;\;0\;\;\;1\;\;\;0\;\;\;0\;\;\;0 \\ 0\;\;\;0\;\;\;0\;\;\;1\;\;\;0\;\;\;0 \\ 0\;\;\;0\;\;\;0\;\;\;0\;\;\;1\;\;\;0 \\ 0\;\;\;0\;\;\;0\;\;\;0\;\;\;0\;\;\;1 \\ \end{array} \right] \text{。} $ (10)

图3为两相静止坐标系与两相旋转坐标系示意图。使M轴与α轴重合,M轴滞后T轴π/2,旋转方向为逆时针。根据2个轴系形成的旋转磁场等效的原则,消去磁动势中的匝数,直接用电流来表示,则

图 3 两相静止坐标系与两相旋转坐标系示意图 Fig. 3 Schematic diagram of two-phase stationary and two-phase rotating coordinate system
$ \left[ \begin{array}{l} {i_d}\\ {i_q} \end{array} \right] = \left[ {\begin{array}{*{20}{c}} {\cos (\theta )}&{{\rm{sin(}}\theta {\rm{)}}}\\ { - \sin (\theta )}&{{\rm{cos(}}\theta {\rm{)}}} \end{array}} \right]\left[ \begin{array}{l} {i_\alpha }\\ {i_\beta } \end{array} \right] \text{,} $ (11)

将上述变换矩阵扩展成6*6方阵,由于零序电流不形成旋转磁场,不用转换,只需在主对角线上增加4个1,使矩阵增加4行4列:

$ {C_{2s/2r}} = \left[ {\begin{array}{*{20}{l}} {{\rm{cos(}}\theta {\rm{)}}}&{{\rm{sin(}}\theta {\rm{)}}}&0&0&0&0\\ { - {\rm{sin(}}\theta {\rm{)}}}&{\cos (\theta )}&0&0&0&0\\ 0&0&1&0&0&0\\ 0&0&0&1&0&0\\ 0&0&0&0&1&0\\ 0&0&0&0&0&1 \end{array}} \right] \text{。} $ (12)

六相静止坐标系与两相旋转坐标系之间的变换如下:

$\small\begin{split} &{C_{6s/2r}} = {C_{2s/2r}} \times {C_{6s/2s}} = \displaystyle\frac{1}{{\sqrt 3 }}\times \\ & \left[{\begin{array}{*{20}{l}} {{\rm{cos(}}\theta {\rm{)}}}&{\cos (\theta - \displaystyle\frac{\text π }{6})}&{{\rm{cos(}}\theta - \displaystyle\frac{{2{\text π}}}{3}{\rm{)}}}\\ {{\rm{ - sin(}}\theta {\rm{)}}}&{{\rm{ - sin}}(\theta - \displaystyle\frac{\text π}{6})}&{{\rm{ - sin(}}\theta {\rm{ - }}\displaystyle\frac{{2{\text π}}}{3}{\rm{)}}}\\ 1&{{\rm{ - }}\sqrt 3 /2}&{ - {\rm{1/2}}}\\ 0&{1/2}&{ - \sqrt 3 /2}\\ 1&0&1\\ 0&1&0 \end{array}}\right.\\ &\left. {\begin{array}{*{20}{l}} {{\rm{cos(}}\theta - \displaystyle\frac{{5{\text π}}}{6}{\rm{)}}}&{{\rm{cos(}}\theta - \displaystyle\frac{{4{\text π}}}{3}{\rm{)}}}&{{\rm{cos(}}\theta - \displaystyle\frac{{3 {\text π}}}{2}{\rm{)}}}\\ {{\rm{ - sin(}}\theta {\rm{ - }}\displaystyle\frac{{5{\text π}}}{6}{\rm{)}}}&{{\rm{ - sin(}}\theta {\rm{ - }}\displaystyle\frac{{4{\text π}}}{3}{\rm{)}}}&{{\rm{ - sin(}}\theta {\rm{ - }}\displaystyle\frac{{3{\text π}}}{2}{\rm{)}}}\\ {\sqrt 3 /2}&{ - {\rm{1/2}}}&{\rm{0}}\\ {{\rm{1/2}}}&{\sqrt 3 /2}&{ - {\rm{1}}}\\ 0&1&0\\ 1&0&1 \end{array}}\right]\text{。} \end{split}$ (13)

通过坐标变换,将六相双Y30°绕组变换到静止αβ坐标系,得到的电压、磁链方程与三相感应电机αβ坐标系下的电压、磁链方程相同,如式(1)和式(2)所示,六相双Y30°绕组感应电机αβ两相静止坐标系下的数学模型表达式如式(3)。

2 六相双Y30°绕组感应电机仿真 2.1 仿真模型构建

图4为六相感应电机控制系统主电路在电力电子仿真软件(PSIM)中搭建的仿真模型,图5为基于六相双Y30°绕组感应电机αβ两相静止坐标系下的数学表达式(3)建立的仿真模型,与上面三相感应电机相同。只是三相感应电机需进行3/2变换,而六相感应电机是基于6/2变换,静止αβ坐标系输入电压表达式如下式:

图 4 六相感应电机控制系统主电路 Fig. 4 Control system main circuit of six-phase induction motor

图 5 六相双Y30°绕组感应电机仿真模型 Fig. 5 Simulation model of six-phase winding induction motor with double Y30°
$\small \begin{gathered} \left[ {\begin{array}{*{20}{l}} {{u_{s\alpha }}} \\ {{u_{s\beta }}} \end{array}} \right] = \frac{1}{{\sqrt 3 }}\times \hfill \\ \left[ {\begin{array}{*{20}{l}} 1&{\sqrt 3 /2}&{ - {\text{1/2}}}&{ - \sqrt 3 /2}&{ - {\text{1/2}}}&{{\text{0}}} \\ 0&{1/2}&{\sqrt 3 /2}&{{\text{1/2}}}&{ - \sqrt 3 /2}&{ - {\text{1}}} \end{array}} \right]\left[ {\begin{array}{*{20}{l}} {{u_A}} \\ {{u_B}} \\ {{u_C}} \\ {{u_D}} \\ {{u_E}} \\ {{u_F}} \end{array}} \right] \text{,} \end{gathered} $ (14)
$\small \left[ \begin{array}{l} {i_{sa}}\\ {i_{ib}}\\ {i_{sc}}\\ {i_{sd}}\\ {i_{se}}\\ {i_{sf}} \end{array} \right] = \frac{1}{{\sqrt 3 }}\left[ {\begin{array}{*{20}{l}} 1&0\\ {\sqrt 3 /2}&{1/2}\\ { - {\rm{1/2}}}&{\sqrt 3 /2}\\ { - \sqrt 3 /2}&{1/2}\\ { - {\rm{1/2}}}&{ - \sqrt 3 /2}\\ 0&{ - 1} \end{array}} \right]\left[ \begin{array}{l} {i_{s\alpha }}\\ {i_{s\beta }} \end{array} \right] \text{。} $ (15)

图6为六相双Y30°绕组感应电机矢量控制原理框图。采用双闭环控制系统,外环为速度环PI控制,内环为电流环PI控制,包含励磁电流iM和转矩电流iT调节器。基于同步旋转dq坐标系实现电压、电流坐标变换,并采用转子间接磁场定向控制。基于表达式(14)将六相静止坐标系下绕组电压变换到两相旋转坐标系,送入闭环控制系统。控制系统为了能快速实现电流跟踪,基于解耦电路实现电压前馈补偿控制。

图 6 六相双Y30°绕组感应电机矢量控制原理框图 Fig. 6 Vector control block diagram of six-phase winding induction motor with double Y30°
2.2 系统仿真结果

图7为六相双Y30°绕组感应电机在转子间接磁场定向控制下MT轴电流参考值与实际值波形。从波形中可以看出,能实现较好的静、动态跟踪,励磁电流及转矩电流静态误差较小。图7(a)给出了线电压uabude波形,图7(b)给出了电机角速度及六相绕组电流随时间变化曲线,图7(c)图7(b)在时间轴局部放大图。从图中可以看出,2套三相绕组电流分别1200对称,B相绕组电流滞后A相绕组300电角度,跟电机实际空间结构相符,从而说明六相双Y30°绕组感应电机数学模型及坐标变换公式的正确性,矢量闭环控制的有效性。

图 7 六相双Y30°绕组感应电机矢量控制仿真波形 Fig. 7 Vector control simulation waveforms of six-phase winding induction motor with double Y30°
3 结 语

本文基于三相感应电机αβ两相静止坐标系下的数学模型及六相双Y30°绕组磁势的等效关系,推导出六相双Y30°绕组感应电机的数学模型及坐标变换表达式,并通过电力电子仿真软件PSIM构建系统仿真模型,仿真结果证明建模方法的正确性及控制策略的有效性。

参考文献
[1]
LYRA R O C, LIPO T A. Torque density improvement in a six-phase induction motor with third harmonic current injection[J]. IEEE Trans. on IAS, 2002, 38(5): 1351-1360.
[2]
BOJOI R, LAZZARI M, PROFUMO F, et al. Digital field oriented control for dual three phase induction motor drives[J]. IEEE Trans Ind App, 2002(IA-02): 818-815.
[3]
BOJOI R, GRIVA G, PROFUMO F. Field oriented control of dual three- phase induction motor drives using a luenberger flux observer[J]. IEEE Trans Ind App, 2006, IA-06: 1253-1260.
[4]
杨金波, 杨贵杰, 李铁才. 六相电压源逆变器PWM算法[J]. 电工技术学报, 2012, 27(7): 205-211.
[5]
郑剑, 肖蕙蕙, 李山. 六相感应电机直接转矩控制系统[J]. 电机与控制应用, 2010, 37(11): 15-19.
ZHENG Jian, XIAO Hui-hui, LI Shan. Direct torque control system of 6-pahxe induction motor[J]. Electric Machines and Control Application, 2010, 37(11): 15-19. DOI:10.3969/j.issn.1673-6540.2010.11.004
[6]
黄守道, 王铭, 钱鋆. 双Y移30度六相PMSM的空间矢量控制[J]. 控制工程, 2011, 12(6): 202-205.
[7]
艾永乐, 王玉梅, KAMPER Marrten. 梯形波相电流驱动六相感应电机建模与控制[J]. 电机与控制学报, 2012, 16(3): 72-77.
AI Yong-le, WANG Yu-mei, KAMPER Marrten. Modeling and controlling of six-phase induction machine driven trapezoidal phase current waveform[J]. Electric machines and control, 2012, 16(3): 72-77.
[8]
艾永乐, 王玉梅, KAMPER Marrten. 基于梯形波相电流驱动的六相交流感应电动机性能[J]. 电工技术学报, 2009, 24(9): 59-63.
AI Yong-le, WANG Yu-mei, KAMPER Marrten. Performances of six-phase induction motor with trapezia phase current waveforms[J]. Transactions of china electro technical society, 2009, 24(9): 59-63. DOI:10.3321/j.issn:1000-6753.2009.09.009