多相感应电机任意非正弦气隙磁密的磁动势分析与应用
来源:叨叨游戏网
第34卷第2期 2016年4月 轻工枕械 Light Industry Machinery Vo1.34 NO.2 Apt.20I6 [自控·检测] DOI:10.3969/j.issn.1005-2895.2016.02.01 1 多相感应电机任意非正弦气隙磁密的 磁动势分析与应用 梅泽挺,蔡卓剑,赵荣祥 (浙江大学电气工程学院,浙江杭州摘310027) 要:多相感应电机转矩密度的提高依赖于非正弦气隙磁密行波的合理构建,由于传统的电流相位、幅值控制策略无 法给出期望的非正弦磁密行波,因此文章研究了控制多相感应电机定子相电流中基波和谐波幅值、相位、频率的大小来 构建期望的气隙磁密波形的方法。首先基于单根裁流导体产生的气隙磁动势波形,推导了任意相数、不对称、包含任意 高次谐波电流波形激励的具有典型绕组结构的多相感应电机的磁动势分布表达式。根据以上磁动势分析结论,多相电 机可通过变电流相序变极。针对气隙磁密行波逼近方波波形的案例,根据所提出的磁动势分析方法分析和预测了理想 空栽励磁电流的波形 将预测波形和方波供电下的励磁电流波形进行了比较,从FFT分析结果来看。仍需微调各次谐波 电压相位,从而实现各次电流相位完全一致。 关键词:多相电机;磁动势分析;任意波形的磁场行波;气隙磁密波形 文献标志码:A 文章编号:1005—2895(2016)02-0049-07 中图分类号:TM343.2 MMF Analysis and Application for Arbitrary Air-Gap Flux Density in Multiphase Induction Machines MEI Zeting,CAI Zhuojian,ZHA0 Rongxiang (College of Electrical Engineering,Zhejiang University,Hangzhou 3 l0027,China) Abstract:Torque density enhancement of muhiphase induction machines(M1Ms)depends on the reasonable construction of non—sinusoidal air—gap flux densitY travelling waveform。Due to the defect of traditional strategy with amplitude and phase shitfs control of phase currents which can t yield the desirable flux density travelling waveform.the method was found to control the phase shifts,amplitude,fi'equencies of phase currents.MMF distribution expressions of MIMs With typical winding stucture and arbitrary phase number based on the MMF wayefrorm produced by one current— carrying conductor was derived.Phase currents were unsymmetrical and contain arbitrary higher harmonics.According to the eonclusions from MMF analysis of MIMs above,the pole palrs of MIMs could be eontrolled by phase shift between adjacent phases.In the case study of the air—gap flux density wave approaching square wave,the magnetizing current waye was analyzed and predicted in the light of the method of MMF analysis.It is found that phase shitfs of voltage harmonies should be finely tuned.SO that the EHrrent harmonies are in phase wih each otther in terms of FFT analyzed results from the comparison of predicted wave with the magnetizing current wave under square wave voltage power supply. Key words:muhiphase induetion machines;magnetic motive force(MMF)analysis;travelling wave with arbitrary waveform;lfux densi ̄wave in the air gap 收稿日期:2015-08-26;修回日期:2015-09-22 基金项目:高等学校博i学科点々项科研基金博导类rJ主课题(20120101110113)。 作者简介:梅泽挺(1991),男,湖北红安人,硕十研究隹,主要研究方向为电机建模机理分析、感应电机控制、电力电予。E—Iliail: meizeta@zju.edu.en通信作者:赵荣祥(1962),男,教授,主要从事变频涮速系统、电力电子技术在储能巾的应 等方面的研 究。E—mail:rongxiang@zju.edu.cn ·50· 轻工机械LightIndustryMachinery 2016年第2期 自1980年以来,多相感应电机的研究便成为电机 方向的一个研究热点。随着功率电子学、微控制器和 电机基本理论的快速发展,这一研究趋势得以加强。 F( )=(n( )一(n( ))) (t)=N( )i(t), (1) 特别是现在电力电子变换技术可以给电机提供可变频 率、幅值、相位的供电电源,因此电机的相数可以不受 (n( )) 1T O dln( )d 。 (2) 式中:i(t)为电流随时间变化的函数; 为沿着圆周的 机械角度,且 E(0,2叮T);N( )被称为绕组函数;n 三相电网的约束。研究表明,多相电机相对于三相电 机,优势在于:拥有更小的脉动转矩幅值、更高的转矩 脉动频率、更好的容错能力和更高的工作效率。据报 道,多相电机已经在舰船推进,机车牵引,更加“电气 ( )为匝数函数,可以用数导体的方法来确定¨ 。 当然,可以将绕组函数用来确定任意节距的载流 线圈产生的磁动势波形F :( ),即 化”的航天器方面找到了用武之地¨ 。 2O世纪90年代,LIPO等学者发现多相电机除了 使气隙磁密逼近正弦波形以外,还可以利用低次谐波 电流注入使气隙磁场波形为非正弦,从而提高电机的 转矩密度。这一现象引起了很多电机学者的兴趣。之 后,TOLIYAT等学者对这一发生在集中整距绕组电机 中的现象进行了深入分析和研究 。 在文献[5—7]中,绕组函数作为电机磁动势分析和 电机参数计算的强有力工具。早期电机研究学者们将 这种方法用于电机绕组的优化设计以及噪声与振动分 析 加]。 目前多相感应电机磁动势的分析工作的局限在于 励磁电流为具体的电流波形(基波+三次谐波,相与 相之间的相位差固定),没有给出多相电机磁动势分 析的完整分析过程以及量化结果。本文从最基础的单 根载流导体产生的磁动势波形开始分析,直至得出具 有典型绕组结构的感应电机所有相绕组合成的总磁动 势表达式。在磁路B—H曲线为非线性的情况下,磁动 势波形和磁密波形并不完全相似,传统的方法并没有给 出气隙磁密为非正弦行波时相电流的基波、谐波幅值和 相位之间的量化关系。文章在研究了任意形状的行波 的基础上,提出了控制多相电机定子各相相电流的基波 和谐波幅值、相位和频率的大小来构建期望的气隙磁密 波形的方法。最后,基于上述分析结果,预测了气隙磁 密逼近方波时电流的波形,并将其与方波供电下的励磁 电流波形进行了比较。从Fn'分析结果来看,方波供电 下的励磁电流波形已经比较接近理想电流波形。 1 多相感应电机的磁动势分析 1.1 多相感应电机磁动势分析基础 在对多相感应电机进行磁动势分析之前,首先要做 一些基本假设 。基于磁场的安培定律,高斯定理, D.W.Novotny推导出了沿着定、转子圆周任意分布的 载流绕组(如图1(a)所示,其中R ,R 为定子内半径和 转子外半径)。在电机气隙中产生的磁动势波形,有 ( )= 1 sin( s Jf ( )(1一 ),细 0≤ < , (3) 【一Nil(t) ot, < <2订。 式中: 表示单个线圈边所包含的匝数, 为线圈的节 距,用弧度表示。公式(3)的Fourier级数和原函数在不 同的情况下可以选择性使用,其原函数如图1(b)所示。 根据绕组函数,可以得到单根载流导体产生的气 隙磁动势波形: F(qS):坐 in 。 (4) 公式(4)的原函数如图1(C)所示。 (a)沿转子外圆周分布的任意节距的线圈 一 f(0(1 1/2腿 f) 0 -U ̄i( -1/2N ̄i(t) 机械角度 ̄/rad 机械角度 ̄/rad (c)单根导体产生的气隙 (b)单个线圈产生的气隙磁动势分布函数 磁动势分布函数 图1 电机圆周任意节距栽流线圈分布图 及其磁动势分布波形 Figure 1 Distribution diagram of arbitrary ^ pitch current—carrying coil around motor periphery and MMF distribution wave [自控·检溯] 梅泽挺,等:多相感应电机任意非正弦气隙磁密的磁动势分析与应用 1.2只考虑基波电流作用时的磁动势分析 由于多相电机的磁动势分析方法是建立在深人挖 掘三相电机绕组结构的数学特点的基础之上,在对多 相电机的磁动势进行分析之前,先熟悉一下三相电机 的典型绕组结构。如图2所示,三相电机的绕组结构 是典型的双层、短距、分布式绕组结构。表1所示为定 子槽号的分配情况。在图2中,定子拥有24槽,4极, 分布式线圈个数q=2,短距线圈都拥有 =5槽(整距 =6), 和。表示电流流人导体的方向是指向纸面 内还是纸面外。绕组圆图相对于绕组平面展开图的优 势是可以看出绕组结构的数学特点,便于做抽象的数 学分析。槽中的每根导体含有Ⅳ 匝,具有典型绕组结 构的感应电机的定子槽数应该满足 Z=2npq。 (5) 式中:z代表定子槽数目;n代表定子相数;p代表极对 数;q代表分布式线圈的个数。 图2 具有双层短距分布式绕组结构的 三相感应电机 Figure 2 Three phase induction machines with double-·layer and short--pitched distributed windings 表1 定子槽号分配表 Tabl 1 DistribUti0n tabl nf stator slots number 磁动势分析将会首先从线圈组的q个载流线圈开 始,有 ( )= g 。=I 1V 警)×£ (6) COS V( 一(i一1) 1。 式中:Hi( )表示一相绕组线圈组中的第i(i:1,2,… q)个载流线圈产生的磁场强度;OL 表示机械短距角 度; =2 ̄r/Z为槽距角;g 表示等效气隙长度。由式 (6)可得线圈组的总磁动势为 HN/s(ci2):Hl( )+月 (ci2)+…+Ha( )= × COS ( 一(q一1)÷ )。 (7) 式中H ( )代表了一相的一个线圈组所产生的总磁 动势。他的物理意义是代表着一相载流绕组产生一个 等效的N/S极。 当涉及到一相绕组的磁动势分析时,引入变量 取代 一÷(q一1) ,则一相绕组的磁动势H ( ) 的分布表达式为 =2q]V if (t) 1 in( 孚)翼× [COS B—COS ( 一 )+COS (卢一2a )一…一 COS v(/3一(2p一1)OL )]。 (8) 式中: =1T/p是整距线圈的机械角度。重新整理式 (8),简化可得 一 ( )si 警)× 望 sin v(/3-(p-1 psin v% 一 n ))。(。 9) 从相绕组的磁动势表达式可以看出: 1)在相绕组的磁动势中只有P的奇数倍的极对 数存在; 2)短距系数和分布系数都存在于相绕组的磁动 势中。 最后,分析由定子所有的相绕组产生的总的磁动 势分布。假设电流为只有对称的正序基波分量,则有 H = 等 1× [sin(∞ + )cos +sin( f一1n 2订+ )c。s ( 一 轻工机械2 丌 + LightIndustryMachinery 2016年第2期 二 一。 + i ( 扣 一 n 耵+ )。。 ( 一 (10) 势分布。基频激励的双层集中整距绕组结构,不同相 数的感应电机的磁动势分布如表2所示。在表2中, 不同转速的磁动势时空谐波的幅值和速度大小都选取 相对于基速旋转的磁动势的幅值和速度大小。最低次 式中:r=1,3,…;,表不相电流的有效值; 表不电流 的角频率; 表示相电流的初始相。 进一步简化式(10),则有 的磁动势谐波的极对数为P,幅值为2√2qN In/'rrg ,它 是引起电机旋转的工作极对数。电机所谓的变极方法 就是要改变幅值最大的工作极对数;正弦磁密的电机 : 等,翌1× (詈+ ))+ 绕组的设计目标就是要削弱其他高次磁动势空间谐波 幅值相对于它的大小。表2中,n表示定子相数,取 n=3,5,7,9,更大的n值没有继续在表中反映。由表 2可以看出:定子中的基频电流可以产生多种极对数 的磁动势谐波,这些磁动势谐波的旋转速度是不同的; 转子中的电流将由于定子电流产生的磁动势谐波在电 距分布式线圈的情况下磁动势谐波含量随着相数的增 [ sin( t+ + 一 t=n sin( ̄oot+ 一 + (詈一 ))]。( ) 机气隙圆周中旋转而感应出来;多相电机在不采取短 (12) 只有满足以下约束条件的极对数存在 =p(2kn±1),k:0,1,2,…。 由式(10)和(12)可得多相电.机基频激励的磁动 加而明显减少,从而磁密波形可以更加逼近正弦。 表2基频激励产生的多相感应电机磁动势行波谐波分布(相移=2'rr/n) Table 2 MMF travelling wave harmonics distribution produced by fundamental frequency current excitation(phase shit=2zfr/n) P 5p 7p 1 1/5 一1/7 ∞ ,p —1/5 1/7 P 9p 1 1/9 ∞ ,p P l ∞ P 1 m ,p 1/9 l/1l 13p 15D 1/13 一1/l5 —1/13 1/15 17p 19p 35p 37p I/l7 一1/19 一1/35 l/37 —1/I7 1/19 —1/35 1/37 11D —l/ll 1lp13p —-1/II1/13 —-1/1I 1/13 19p 21D 一1/19 I/21 1/19 1/21 27p 29p —l/27 1/29 一l/27 1/29 : 1.3高频电流激励的磁动势分析 邻相的相电流的相位差在3,5倍基频电流激励时都为 6"rr/n,这时磁动势谐波分布如表3—4所示。 表3 3倍基频激励产生的多相感应电机 采用简化的时域、空间域双重傅里叶级数来分析 高频电流激励的磁动势,则有: ,£)=∑g ( ) ( ), "=l Ⅱ=∞ 磁动势行波谐波分布 Table 3 MMF travelling wave harmonics distribution gm( )=∑Avsin( + ), im( )=∑A sin( t+ )。 式中:g ( )表示第m相相绕组的空间分布函数; i (t)表示第m相相绕组的相电流随时间变化的函数; A 是周期函数傅里叶级数第//'项对应的系数和相 位;4 , 是周期函数傅里叶级数第 项对应的系数和 produced by 3 times fundamental frequency current excitation 相位。先考虑3,5倍基频电流激励不同相数的整距集 中式绕组结构的感应电机的磁动势分布。首先假设相 [自控·检测] 梅泽挺,等:多相感应电机任意非正弦气隙磁密的磁动势分析与应用 ·53· 表4 5倍基频激励产生的多相感应电机 磁动势行波谐波分布 Table 4 MMF travelling wave harmonics distribution …在公式(13)中,a=exp(J×2"rr/n),n个不对称相 量 , ,…,, 的第h+1个对称分量组为J『 口 ,, a Js 。当非对称相量被转换成了对称相量,与 produced by 5 times fundamental frequency current excitation 前面相同的分析过程将同样适用于不对称多相电流存 在的情形。 ● 『 = ● : - : ● L_ . 在表3~4中,磁动势谐波的幅值和速度大小都选 取相对于工作极对数的磁动势的幅值和速度大小,最 (13) 当一台n相电机的定子相电流不对称时,与工作极 低次的磁动势的极对数为3p,它是工作极对数,幅值 对数不等的低次磁动势谐波将会出现并且会产生脉动 转矩。假若只有零序相量组,那么将没有磁动势行波产 生,只有驻波存在。当存在第h个对称分量 时,意味 着第h个对称分量组, n , ^,…,a ’, ^一定存 分别为2 qNc13n/3wg ,2√2q ,5/3wg。。13,,5分别 是3,5次正弦电流谐波的有效值。表3~4表明,磁动 势谐波的旋转速度等于电流频率/极对数,可以通过控 制相电流的频率和工作极对数来控制磁动势谐波的旋 在。这时,只有满足以下约束条件的磁动势行波存在: =P(kn±h) k=0,1,2,…;h:0,1,2,…,n一1。 转速度,同时,磁动势谐波的极对数只与相数以及相邻 相电流之间的相位差有关。 1.4不对称多相电流情形下的磁动势分析 (14) 除了要满足式(14)外,还必须是正奇数。显然, 当相电流的基波分量不对称时,意味着相邻相电 流的相位差不全是2 ̄h/n,h=0,1,2,…,n一1。相电 当相数越多时,对称分量的数目也越多,这就为多相电 机的变极提供了另外一种方法:改变相邻相电流之间 的相位差而不是改变电机端部接线方式。整理表1~ 3可得表5。 流幅值不全相等,这时可以使用多相对称变换来解决 电流不对称的问题 12-13]。 表5 多相感应电机相电流不同相序时的磁动势分布 Table 5 MMF travelling wave harmonics distribution produced by phase currents with different phase sequences 从表5中可以看出,5相电机中当相电流之间的 相位差从2 ̄r/5变为6-rr/5时,工作极对数将从P对极 变为3p对极,但是磁动势谐波含量加重。同样的分析 适用于更高相数的多相电机。 以来,众多电机研究学者就将其注意力聚焦在它上面。 非正弦磁密行波波形在多相电机中实现,是由于多相 电机中电流对称分量数目变多,改变相邻相电流之间 的相位差就可以改变多相电机的工作极对数,当电流 中存在基波、高次谐波并且它们各自的相邻相电流之 间的相位差不一致时,就可以实现不同工作极对数的 磁动势具有相同的旋转速度以及相位。当气隙磁密波 形不是正弦波,而是平顶波、阶梯波时,理论分析表明 2任意非正弦气隙磁密行波波形的实现 2.1 非正弦气隙磁密行波波形的特点 自从发现多相电机可以通过注入电流谐波,实现 气隙磁密为非正弦行波,从而提高转矩密度这个事实 轻工机械LightIndustryMachinery 2016年第2期 这样的气隙磁密行波的确可以提升转矩密度。然而, 些高次的磁动势行波谐波只会产生脉动的转矩,增加 功率损耗。鉴于以上因素,我们可以利用的电流高次 谐波次数并不是很多。 磁通密度,r t2)f 非正弦的行波波形是通过实验来优化各次谐波电流大 小 】引,因此,有必要研究任意的一维行波的数学特点, 从而在理论上对实际如何优化磁密波形给予支持。 图3所示为一种典型的磁密行波波形厂( )。假 设行波的角速度为09 ,通过观察发现,行波在不同 的时刻t 和£:时,它们的波形是相同的。在t:时刻的 、、 = 波形可以看成t。时刻的波形向右移动了b=∞ /p(t 一t.)。这种特点在数学上可以表示为: 厂-( )=∑(0 c0s +bksin kp,p)= ∑Aksin( (f)+ ); (15) 一b)= ∑(akcos ( 一6)+bksin ( 一6))= ∑A sin( ( 一6)+ ); (16) 厂((f)一 P f): =1 akcos印( 一鲁 p )+bksin ( 一警 p ,)) 一∑A sin( ((£,。t—p )一 )。 (17) 式中 )的周期T:A=2,n-/p,人为磁场行波波形的 波长,代表的是机械角度,仅由电机的绕组结构决定; 代表不同极对数时正弦行波的相位;a ,b ,A 为傅 里叶系数。 根据式(17),可以看出非正弦形状的任意行波波 形其本质是一系列具有相同转速、不同极对数的正弦 行波的组成,且有: sin( ( 一p )一 ) 詈{ sin( (㈠) ).c。s( ( 一 (㈠))))。 (18) 式中:k#cn,C=1,2,3,…。 联立式(11),式(14)~(18)就可以知道如何控制 定子电流逼近任意形状的行波。由式(17)一(18)可 知,要想构造非正弦的磁场行波,高次电流时间谐波必 须注入。在现实中必须承认任意行波波形只能逼近, 无法合成,这是由于定子有限的相数,由逆变器控制输 出的电流波形不可能没有高次谐波,一种频率的电流 可以产生多种高次的磁动势行波谐波(表1~4)。这 一 L /2 \ l/., ;一 图3 电机磁动势谐波分析中典型的磁密行波波形 Figure 3 Typical flux density travelling wave in MMF harmonics analysis of electrical machines 2.2非正弦气隙磁密导致的电流过冲现象 仅从改善气隙磁密这一点来看,需要磁密波形逼 近方波如图4所示。但是,从磁动势分析的结论可知, 逼近方波要求集中整距绕组感应电机中的电流谐波幅 值满足,l=一13=15=…, 1= 3= 5=…=0。这 样的话,电流在零点附近会有尖峰出现。这一点既可 以从图5的理论分析结果中看出,也得到了图6¨ 实 验波形的验证。 磁通密度/1r 图4方波波形及其傅里叶级数 Figure 4 Square wave and its Fourier series 6, 4, 蛋21 虐0 羹.2/ .4, 一61 图5 空载时定子电流的预测波形 Figure 5 Predictive wave of stator current at no load 图6(b)所示是110 kW九相集中整距绕组电机在 空载、外施加电压为900 V,50 Hz的图6(a)中方波电 压下电流的波形,和图5有一定差别的原因是实验电