无刷直流电动机线电流的解析计算
2022-01-13
来源:小奈知识网
….堕 皇 … ! . 苎_7 ………………………………………一 … 无刷直流电动机线电流的解析计算 高庆嘉,白越,王鹤 (中国科学院长春光学精密机械与物理研究所,吉林长春130033) 摘要:为了能够方便准确地预测无刷直流电动机的机械特性,对电动机的线电流进行了解析计算。根据电动 机电路模型和电压平衡方程式,在不忽略电感和电阻的情况下对换相过程进行分析,指出了电磁时间常数和一个状 态角换相周期的比值关系是影响相电流瞬时值变化的主要因素。结合相电流瞬时值变化特点将线电流的计算模型 分为两种,在此基础上推导了不同情况下电流瞬时值解析表达式,最后借助Matlab软件得到了电动机线电流的解析 解。与两台典型无刷直流电动机实测数据进行对比验证,该方法计算精度较高计算且简便、快速,可应用于工程实际。 关键词:无刷直流电动机;线电流;解析计算;机械特性;电磁时间常数 中图分类号:TM33 文献标识码:A 文章编号:l004—7018(2011)07—00O1一o4 Analytical Calculation of Line Current for Brushless DC Motor CA0 Qing-jia,BAl Yue,WANG He (Changchun Institute of Optics,Fine Mechanics and Physics,Chinese Academy of Sciences,Changchun 130033,China) Abstract:In order to easily and accurately predict mechanical characteristics of brushless DC motors,the analytical calculation of the line current was done.Firstly,the commutation process including the case of induction and resistance was analyzed based on the motor circuit model and voltage balance equation.It was analyzed that the ratio of electromagnetic time constant and a commutation cycle period at a state angle was the main reason of the three-phase transient current vari— ations,which makes line current calculation model divided into two parts.Then,the analytical expressions of the phase tran— sient current variations were derived.At last,line current analytical solutions were given by the means of Matlab.The calcu・ lating resuhs of two typical prototypes were compared and veriifed with experimental data.111e error was less than 5%.and the method was easy and fast and had a hi【gh precision.As a result,the method in this paper can provide a reference for the motor design and application. Key words:brushless DC motor;line current;analytical calculation;mechanical characteristic;the electromagnetic time constant 0引 言 计内容的书籍看,基本都采用忽略电感的直流电动i 机分析模型 ,而在分析转矩脉动方面的文章却 无刷直流电动机具有结构简单、运行效率高、转 忽略电阻的存在 _『 J。这是由于同时求取计及电阻; 矩密度大和可靠性好等诸多优点,已广泛应用于各 和电感的电动机电路模型较困难。但元刷直流电动 种伺服控制和驱动等场合¨ 。机械特性是电动机 机的电流是由相绕组通过逆变电路按照一定次序轮 设计和应用最重要的特性之一,不同负载下线电流 流换相产生,忽略电感或电阻势必导致电流计算值 的变化表征了无刷直流电动机机械特性的变化。借 和实际值存在一定偏差,进而影响到实际应用。文j一 助测功机等仪器的试验方法可以获得较准确的线电 献[9]分析了计及绕组电感的无刷直流电动机机械;莉 流,但在某些应用场合如电动机在线测量、高速运转 特性,但平均电流积分区间不合适;文献[10—11]分 暴 以及大电流短时工作等,该方法将受到一定的限制。 析了计及绕组电感的无刷直流电动机电路模型和平;曼 利用电动机数学模型并借助仿真工具可以获得数值 均电流计算公式,但忽略了U<4E这一先决条件。 解,但人们还是期望得到解析解,这是由于解析解能 另外,上述文献中均忽视了电磁时问常数 对换相 够描述电动机内在参数与外部特性之间的函数关 电流的影响,这对于有些无槽电动机、空心杯电动机;嚣 系,为电动机设计和性能预测提供简洁的计算方法。 等较小电磁时间常数的电动机来说,将不能真实地 蟹 从目前正式出版的包含永磁无刷直流电动机设 描述换相过程和电流变化,从而导致错误的电动机i钎 线电流解析解。 i弄 收稿日期:2011~02—27 本文以常见的两相导通星型三相六状态无刷直 基金项目:中国科学院“十一五”基金(No.61501.02.03.04) 国家自然科学基金(50905174) 根据电蠹蓍 动机电路模型对换相过程进行分析,指出了_ Il …t 绥 嬲 …………………………………………………………………~ 一 .:-_.一 一:二……….: :. …. 研 微持电棚 2011年第7期 电磁时间常数和一个状态角换相周期的比值关系是 j影响相电流瞬时值变化的主要因素。根据电磁时间 常数和一个状态角换相周期的比值关系,将电流计 :算模型分为两种,并推导了相应的电流瞬时值解析 为电动机转速,P为磁极对数。不同状态下绕组及 功率管导通顺序如表1所示。 表1不同导通状态下绕组及功率管导通顺序 状态 I J Ⅱ I Ⅲ l Ⅳ V Ⅵ 导通相 A/C l B/C I B/A I C/A c,B A/B 功率管 T1T2 I T2T3 I T3T4 l T4T5 1'5T6 Tl 表达式,最后借助Matlab软件计算电动机线电流的 j解析解。通过与两台典型无刷直流电动机实测数据 进行对比验证,证明了本文分析方法和计算公式的 正确性,且计算简便、快速,可为无刷直流电动机设 由上述分析可知,只要对任意一个状态进行分 析,就可求得整个换相周期的解,下面的讨论基于电 动机从A/C相导通结束到B/C相导通换相过程,以 计和应用提供一定参考。 1电动机分析模型 本文以两相导通星型三相六状态永磁无刷直流 电动机为研究对象,电动机电路模型如图1所示,其 中,U为功率桥两端的直流电压,T1~T6为MOS— FET功率管,D1一D6为续流二极管,尺和£分别为 一相绕组的等效电阻和等效电感, ,e,i分别为各相 绕组的端电压、反电势和电流的瞬时值,电流以流向 绕组中心点为正电流,相电压以中心点为参考点。 直流电动机换相过程的换流时间通常在毫秒级,而 组的电压平衡方程式。电动机正常运行时,星型连 嚣 析 蓁 [三 ]=[蚕 兰][莩]+[§兰蔓] [萋]+[(三1l ], ■一个状态角对应的换相时间 = ‘箸= ,式中,n 及B/C相导通一个状态角内过程。 2换相过程分析 如图2所示,电动机换相过程主要包括换相区 和导通区两部分。在状态I结束、状态Ⅱ开始时,开 关 保持开通状态不变,T1关断的同时 打开, 此时C相绕组电流持续, 相绕组电流由零开始上 升,等效电路如图2a所示。由于电感存在,A相电 流不可能马上突变为零,而是通过与T4反并联的二 极管D4续流经过一段时间t,下降至零,由于二极 管阻断特性,D4关断,换相区结束,此时电动机只有 相和c相绕组导通。进入导通区,等效电路如图 2b所示,当到达 时刻一个换相状态结束,下一个 换相状态开始。 翁u (a)换相区 (b)导通区 图2电动机等效电路 下面对换相区和导通区内的电流变化分别进行 分析。 2.1换相区(0~t,时刻) 换相区内A相电流经过t 时间降为零,由图2a 可列出该时间段的电压平衡方程式: 盹+ 誓+en 一 誓 。=0 巩+£鲁 一 誓 = (2) i。+i6+i =0 其中,初始条件: ia(0)=一ic(0)=,01 (0)=0 } (3) e。=e6 一e =E J 三相电流表达式: d 。 +2E+3Ri 一=一~ 3£ di6 2(U—E+3Ri ) d£ (4) di U一4E+3Ri 一==一一 出 3 ….堕壁熏 … 苎 …由式(2)和式(3)解得: . …………………………………… …: 4a所示;另一种情况是i 在 时刻之前已经达到 一U+2E 3too+U-I-2E上 a 一■ 一+—— 广(1 . 。 (5) ,0,维持一段时间后进人下一换相状态,如图4b所 示。根据电磁时间常数的定义,i 由零增加到最大 ÷) 值所需时间为 ,当 < 时第一种情况发 ! 4E—Uc ■ 一+—— 广.U一4E一3Rto J l 生,反之,第二种情况发生。设 = ,则两种情 式中: 为电动机电磁时间常数,下=去。 况分别对应 ∈(0,1)和 ∈(1,+∞)两个区间,因 在换相结束时刻t。时,B相和C相电流值: (tl= (6) 根据式(5)可解得i。降为零和i 升到to所需 时间ta,f和tbf: ,ta,f刊n(1+ ) (7) 刊n(1+ ) (8) 实际中,由于电动机工作电压、转速以及负载的 不同, 。和 变化率可能不同,因此可分为三种情况: (1)鲁= dib,即 降为零时, 达到稳态值,0, 如图3a所示,此时ta,f=f¨,得U=4E+3矾。 (2) di.dib< , 。降为零时, 未达到稳态值,0, 如图3b所示,此时,U<4E+3mo。 (3) di.> dib, 。未降为零时, 早已达到稳态值 ,n,如图3c所示.此时.U>4E+3彤n。 图3一个换相角内三相电流变化 电磁时间常数r表征了电动机电流从零增加到 最大值的63.2%所用的时间。因此,f¨不但与电磁 时间常数.r有关,还跟一个换相状态角的时间 有 关联。通常情况下,r的数值一般在毫秒和微妙间, 而续流时间 。. 为微秒级甚至更小,因此,多数电动 机换相过程是按照图3b所描述的电流变化进行的, 本文即针对图3b情况进行分析,其他两个情况分析 过程与之类似。 2.2导通区(t -T时刻) 在导通区内,A相电路断开, 相电流继续上 升,导通区内i 变化曲线根据 和 的不同可分为 两种情况,如图4所示。一种情况是i 从t 时刻到 时刻电流持续增大,直到 时刻恰好达到,n,如图 此导通区内i 的解分两种情况分析计算。 (a) (b) 图4 i 在一个状态角内变化曲线 导通区内以A相电流降为零时间为零时刻,结 合图2b,可得B/C相在导通区的电压平衡方程: =2(Ri6+E+ dib) (9) 初始条件:ib(O)=一 。(0)=,。 当tx E(1,+∞),对应的导通过程如图4a所示, 由式(9)可得口相和G相电流: . . 。 — U一2E +—(U+2E)(R,n+2E—U)一 上 (10) 初始电流解: ^ ^一上 ,0=, ・ (11) 2一e一 式中:lr为忽略电感的电流值,, = 。 当 ∈(0,1)时,对应的导通过程如图4b所示。 B相和C相电流的解表达式: U一2E (U-I-2E)(mo+2E— ) 2R 。 2R(3R,n+U+2E) 初始电流解: ,0:lr.2—e—-—2e百---r -- (13) 2一e一 式中:t:为以t 为零时刻起点i 达到稳态值,0的时 间,t2=0 6 32一t1。 .3线电流的解析计算 线电流即流经电源母线上的平均电流值,下面 对线电流计算公式进行讨论。由于A相电流是通ii 一…t ■ {目 研 …………………………………………………一……………~ j岐持电柏 211年第7期 : 一 一:.:……0…一:: :. :…一 .过开关T2与续流二极管1)4组成一个回路,并没有 流过电源母线,因此,在讨论的B/C相导通状态角 载点转速4 468 r/rain,实测电流值为0.241 A。由 于 =9.456,利用式(16)计算,线电流0.231 A,误 差4.3%,文献[10]计算结果为0.218 A,误差 内,线电流,只与 相平均电流有关,i 的平均值即 电动机的线电流: 9.5%,本文计算公式精度更高。按式(16)计算,线 ,=Ib 寺 d 当 ∈(1,+∞)时, ,= ( ibd T-t1 d )(14) (15) 电流一转速计算值与实测值如图5a所示。 无齿槽无刷直流电动机参数: =28 V,P=3, =0.35 n,L=0.09 mH, =0.005 313 V/(r・min), J『=0.257 1 ms。该电机在输出力矩0.15 N・m时, 结合式(5)、式(10)、式(11)、式(15),利用 Maflab容易求得,,整理后,表达式: ,: 兰 二 !二!±!竺: +l 3RT _ (7I+3T一3tl一,re一 ) 2T(3Rio+ +2 ) 。(r一3 +3t1一下e一_ )Elo Tf3RIo+ +2E) 。 (÷ +E )(T一£1+7"e一_t— ) —一(16) 当 ∈(0,1)时, Jr= ( d + d + T-t1 )(17) 结合式(5)、式(10)、式(13)、式(17),利用 Matlab容易求得,表达式: ,——下~( 一志)Io+———— .(2U一2E)( r——~+ e一 ) + (r+3 一3t1一Te- ̄Tg)Vlo — 十 一T—t1 (Jr一3志+3tl—Te-- ̄-)EIo —— 一十 蟊 嘉i 3 R/o U 2E)+ + 4验证与讨论 机i 蓄 为了验证上述分析方法和计算公式的正确性, 燕;分别以定子有齿槽无刷直流电动机和无齿槽无刷直 薜i流电动机两台电动机的实测数据进行对比验证,这 群 两台电动机分别代表了 ∈(1,+o。)和 ∈(0,I)两 异 种情况。 定子有槽无刷直流电动机参数:U=329 V,P= 4,R=32 Q, =107 mH, :3.343 8 ms,Ke=0.055 3 ■v/(r.min)。该电机在输出力矩0.12 N.m时,负 负载点转速4 760 r/min,实测电流值为2.99 A,由 于 =0.581,利用式(18)计算,线电流3.08 A,误差 为3.0l%。按式(18)计算,线电流一转速计算值与 实测值如图5b所示。 7 000r 6 ooo ̄.、未计电感计算值 5 000 、\√ 4 000l 值、\ 3 000lf 计 值 2 O0 —— 广— —] IIA 胜 (a)按式(16)计算 (b)按式(18)计算 图5线电流一转速计算值与实测值 通过上述验证比较,可见两台电动机的线电流 计算值与实测值较吻合,计算精度较高,并且线电流 的整个计算过程简便、快速。计算结果与实测值间 的误差可能与电动机反电势实际波形和设定条件不 完全相同有关,此外还与电感值和电阻值的测量准 确度有关。 定子有槽无刷直流电动机在负载点转速4 468 r/rain时,实测电流值为0.241 A,不计及电感的线 电流1.28 A,是实测电流的5.3倍,两者相差较大, 此时 =9.456。无槽电动机在负载点转速4 760 r/ min时,不计及电感计算的线电流3.85 A,是实测电 流的1.3倍,两者相差相对较小,此时/z=0.581。由 此可见, 越大,电流实际值与不计电感电流计算值 相差越大。另外,由图5可以看出,空载转速时绕组 电感对电流影响较小,其计算值与实测值基本一致, 随着力矩增大,转速下降,忽略绕组电感的线电流计 算值较实测值偏差越来越大,而采用计及绕组电感 的模型计算线电流时,测试值与计算值比较吻合。 5结论 (1)无刷直流电动机的换相过程主要包括换相 区和导通区两部分,电磁时间常数与一个状态角换 相周期比值关系决定了导通区内瞬态电流具有不同 的解析表达式,进而影响线电流的解析解。 (2) ∈(1,+o。)时,瞬态电流由两段分段函数 组成, ∈(0,1)时由三段分段函数组成。忽略电感 (下转第24页) .... …………………………………………定子槽数和气堕壁 …2…011. .塑… (2)永磁体结构方面,表贴式结构具有机械对 称性好、安装方便、成本较低等优点,在高速永磁电 机的设计中较为常用,所以本设计采用了表贴式永 磁体结构。 (3)空载条件下,槽数越多,转子转矩脉动越 Application Society Annual Meeting,Conference Recond of the 1992 IEEE.1992,I:265—270. nder A,Schneider T,Klohr M.Fixation of Buried and mufface- [6] Bimounted magnets in high-speed permanent-magnet synchronous machines[J].IEEE Transactions On Industy Applications。2006。 42(4):1031-1037. 小。额定负载情况下,随着槽数的增多,转子涡流损 耗逐步减少。 [7] 王继强,王凤翔,鲍文博,等.高速永磁电机转子设计与强度分 析[J].中国电机工程学报,2005,25(15):140. [8] 王凤翔.高速电机的设计特点及相关技术研究[J】.沈阳工业 大学学报,2006,28(3):258—264. [9] 王继强.高速永磁电机的机械和电磁特性研究[D].沈阳:沈 阳工业大学,2007. J R.Miller T J E.Demgn of bmsMe ̄permanent mag- [1O] Hendemhot(4)空载条件下,转子损耗随气隙增大而减小。 同时。气隙长度达到2 mm后,再增加气隙长度,对 转子损耗的削弱效果不明显。 参考文献 [I]沈建新,费伟中,陈利根.气隙磁场波形及磁瓦充磁方式对无 刷直流电动机性能的影响[J].微特电机,2006(6):7-9. net motors[M].Magna Physics PuMishing&Clarendon Press. 1994. [2]周风争,沈建新,林瑞光.从电机设计的角度减少高速永磁电 机转子损耗[J].浙江大学学报,2007,41(9):1588—1591. [3]周凤争,沈建新,王凯.转子结构对高速无刷电机转子涡流损 耗的影响[J].浙江大学学报,2008,42(9):1588—1590. [4]余莉.高速永磁无刷直流电机性能分析与设计的研究[13].南 京:东南大学,2007. Zhu z Q,Xia Z P,Wu L J,et 1.Ianfluence fSloot and Pole Num- ber Combination on Radia】Force and Vibration Modes in Frae- tional Slot PM Brushless Machines Ha ̄ng Single—and Double- layer Windings[C]//IEEE Energy Conversion Congress and Ex. position,San Jose.2009:3316—3323. [12] 陈世坤.电机设计[M].北京:机械工业出版社,2000:10-17. [13] 周凤争,沈建新,王凯,等.高速永磁无刷直流电动机拓扑结 构的研究[J].微特电机,2007(10):1—3. [5]Hippner M,Harley R G.Looking for an optimal rotor ofr high— speed permanent—magnet synchronous machine[C]//Industry (上接第4页) [6] 叶金虎.现代无刷直流永磁电动机的原理和设计[M].北京: 科学出版社,2007. [7] 罗宏浩,吴峻,赵宏涛,等.永磁无刷直流电机换相控制研究 的电动机线电流计算值与实测值的误差随着 增大 而增大,两者近似成正比关系。 (3)借助Matlab软件对电动机线电流进行解 析计算,计算精度较高,且简便、快速,可为元刷直流 电动机的设计和应用提供一定参考。 参考文献 [1] 吴一辉,高庆嘉,白越,等.反作用飞轮驱动电机的电磁设计 [J].光学精密工程,2010,18(6):1319—1325. [J].中国电机工程学报,2008,28(24):108—122. [8] 王淑红,熊光煜.无刷直流电机换相转矩脉动减小及动态仿真 [J].电机与控制学报,20o8,12(2):169—173. [9] 王晋,陶桂林,周理兵,等.基于换相过程分析的无刷直流电动 机机械特性的研究[J].中国电机工程学报,2005,25(14):141 —145. [1O] 谭建成.无刷直流电动机换相分析和电流转矩解析表达式 [J].微电机,2010,43(5):13—18. [2] 高庆嘉,白越,吴一辉,等.姿控飞轮用永磁无刷直流电动机电 磁设计与分析[J].微特电机,2009,10:4—6. [3]张琛,直流无刷电动机原理及应用[M].北京:机械工业出版 社,2006. 李鲲鹏,胡虔生,黄允凯.计及绕组电感的永磁无刷直流电动 机电路模型及其分析[J].中国电机工程学报,2004,24(1): 76-80. [4]Hendershot Jr,J R,Miller T J E.Design of Brushless permanent -magnet motol ̄[M].Oxford University Press,1995. 作者简介:高庆嘉(1983一),男,研究实习员,主要从事永磁电动 机及飞行器设计。 [5]李钟明,刘卫国,刘景林,等.稀土永磁电机[M].北京:国防 工业出版社,]999. (上接第l2页) [3] 吴德荣,李景川.遗传算法在电机系列优化设计中的应用研究 [J].西安交通大学学报,1999(2):16-20. [4]唐任远.现代永磁电机一理论与设计[M].北京:机械工业出 版社,1997. 版社,1998. [8]李鲲鹏,胡虏生.遗传算法电机优化设计简介[J].微特电机, 2001(4):32—33. [9]周克定.遗传算法的基本理论及其在电机优化设计中的应用 [J].湖北工学院学报,1995(4):l一6. [1O]魏静微.小功率永磁电机原理、设计与应用[M].北京:机械 工业出版社,2009. [5]陈金涛,辜承林.基于遗传算法的永磁直流电机优化设计[J]. 微电机,2006,39(8):21—24,81. [6]陈世坤.电机设计(第2版)[M].北京:机械工业出版社, 一 [7]邢文训,谢金星.现代优化计算方法[M].北京:清华大学出 2Oo4. 作者简介:王广生(1986一)男,硕士研究生,研究方向为电力电 子与电力传动,电机控制。