第44卷2011年第4期4月M ICR OM OTOR SV ol 44.N o 4A pr 2011收稿日期:2010-04-07基金项目:西北工业大学研究生创业种子基金项目作者简介:李延升(1983),男,博士研究生,研究方向为电机与电器。
E-m a i:lliyanchao mm@yahoo .co 窦满峰(1967),男,教授,博导,研究方向为电机与电器。
对转式永磁无刷直流电机的建模与仿真李延升,窦满峰,雷金莉(西北工业大学,西安 710072)摘 要:该文根据对转式与普通永磁无刷直流电机区别,建立了对转永磁无刷直流电机的数学模型,采用M atlab /S i m u li nk 仿真软件建立了电机的仿真模型,并对电机带螺旋桨负载进行仿真分析。
仿真结果表明:仿真波形与理论分析基本一致,验证该模型的有效性,为对转式永磁无刷直流电机的控制算法研究提供了工具。
关键词:对转式;无刷直流电机;建模;仿真中图分类号:TM 36+1 文献标志码:A 文章编号:1001-6848(2011)04-0019-04M odeli ng and Si m ul ati on of the Contra -rotati ng BLDC M otor Control Syste mLI Yansheng ,DOU M anfeng ,LE I Jinli(N ort h w estern P oly technical University ,X i an 710072,China )Abst ract :Contra -rotati n g per m anentm agnet br ush less DC m otor uses per m anentm agnet as the ou ter rotor ,the ar m ature w inding as the i n ner rotor ,both inner and outer rotor i n teracts on the reverse ro tation by m eans of t h e m agne tic force .Based on the ana l y sis of the m athe m atica lm odel o f contra -rotating BLCDM,the mode l of BLDC M w as estab lished by the m odu lar design in M atlab /S i m ulink ,and the si m ulati o n experi m ent w as acco m p li s hed w ith a pr ope ller loads .The si m ulati o n resu lts are consistentw ith t h e theory analysis ,and the m ethod is va li d .The para m eter of th ism ethod is suitable for verif y ing the reasonability o f other contr o l algo -rit h m s and provides a ne w w ay fo r further research o f the con tra -rotati n g BLDC M.K ey w ords :contra -rotati n g ;BLCDM;m odeli n g ;si m u lati o n0 引 言对转式无刷直流电机直接驱动对转螺旋桨,在水下航行器中广泛应用[1]。
它与普通永磁无刷直流电机比较,除永磁体部分可以旋转,电枢部分也相对静止部分旋转,即电磁转矩驱动两个转子朝相反的方向旋转。
以电枢部分为参照系来观察永磁体部分的旋转行为,可以发现对转式永磁无刷直流电机与普通的永磁无刷直流电动机的电流方程、电压平衡方程一致,数学模型中仅仅多了一个运动方程[2]。
根据这一思路,本文根据对转永磁无刷直流电机的数学模型,在S i m u li n k 软件中建立仿真模型,并对其进行仿真分析。
1 对转式BLDC M 数学模型无刷直流电机的基本物理量有电磁转矩、电枢电流、反电动势和转速等[3],这些物理量的计算与电机的气隙磁场分布、绕组形式有十分密切的关系。
对于稀土永磁无刷直流电动机,其气隙磁场波形可以为方波也可以为正弦波或梯形波,这与选用电机的磁路结构和永磁体的形状有关。
本文研究的对转式永磁无刷直流电机,其气隙磁场波形为方波,绕组中感应电动势为梯形波,采用方波电流驱动。
在分析和仿真控制系统时,可直接利用电机原有的相变量来建立数学模型,既方便,又能获得准确结果。
假定永磁无刷直流电机工作在二相导通星形三相六状态下,工作过程中磁路不饱和,不计涡流和磁滞损耗,三相绕组完全对称,那么三相绕组的电压平衡方程式为:U a U b U c=R 000R 000Ri a i b i c+L -M 000L -M00L -Md i ad td i bd td i c d t+e ae b e c(1)44卷当三相绕组星形连接,没有中线时,其电流方程为:i a +i b +i c =0(2)式中,U a 、U b 、U c 为绕组相电压(V );i a 、i b 、i c 为绕组相电流(A );e a 、e b 、e c 为绕组反电动势(V );R 为每相绕组的电阻( );L 为每相绕组的自感(H );M 为两相绕组间互感(H )。
根据作用力与反作用力定律,对转式永磁无刷直流电机的内外转子受到大小相等、方向相反的磁场力[4]。
所以,电磁转矩表达式可写为:T em =e a i a +e b i b +e c i cr 1+ r2(3)其中, r 1、 r 3为转子的角速度(1表示电枢转子,2表示永磁体转子)。
三相绕组的反电势波形是梯形波,大小可表示为:e a =k ef a ( r1+ r 2)( r 1+ r2)e b =k e f b ( r 1+ r2)( r 1+ r 2)e c =k e f c ( r1+ r 2)( r 1+ r2)(4)k e 为相反电势系数; r1、 r 2为转子转角。
f a ( r1+ r2)、f b ( r 1+ r2)、f c ( r1+ r 2)表示反电势的波形,它们是关于 r1和 r2的函数,其最大值为1、最小值为-1。
如图1所示:图1 反电势波形图增加电枢转子的机械运动方程,则对转式永磁无刷直流电机的运动方程为:T em -T 01-B 1 2r1=J 1d r 1d t T em -T 02-B 2 2r2=J 2d r 2d t(5)T 01、T 02为转子的摩擦转矩;B 1 2r 1、B 2 2r2为螺旋桨负载转矩;J 1、J 2为转子的转动惯量。
在对转式永磁无刷直流电机运行过程中中,内外转子在任何时候的电磁转矩都相同,那么假设:(a)摩擦转矩T 01=T 02:由电机的运动方程可以得出,当螺旋桨负载的系数B 1=B 2时,由于内外转子转动惯量J 1 J 2,所以电机起动时,两转子的角速度变化量不同,电机稳定时,两转子速度大小相等,方向相反;当系数B 1 B 2时,两转子角速度变化量不同,转速也不同。
(b)摩擦转矩T 01 T 02:由于电机的摩擦转矩是定值,与螺旋桨负载转矩相比较,可以忽略,所以这种情况和上述a 假设基本相同。
2 仿真模型2 1 电机电压方程模型根据对转式电机无刷直流电机的电压平衡方程式(1)可以得出电机等效模型如图2所示,图2 电压平衡方程模快2 2 反电势计算模块由电压平衡方程式(1)可知,要获得三相电流信号i a 、i b 、i c ,必需首先求得三相反电动势信号e a 、e b 、e c 。
为获得恒定的电磁转矩,将永磁无刷直流电机的反电动势设计为梯形波,其平顶宽度大于120 ,梯形波的幅值与电机转速成正比,计算公式如式(4)所示。
对转式永磁无刷直流电机在运行过程中,内转子与外转子的相对位置不同,电机的反电动势的方向不同。
以二相导通星形三相六状态为例,反电势计算可利用分段线形法模拟实现[5]。
根据两个转子的相对位置,将一个运行周期0-360 分为6个阶段,每60 为一个换向阶段,每一相的每一个运行阶段都可用一段直线进行表示,根据某一时刻的两转子之间夹角和转速信号,确定该时刻各相所处的运行状态,通过直线方程即可求得反电动势波形。
对转式无刷直流电机反电势计算模块如图3所示。
图3 反电势计算模块204期李延升等:对转式永磁无刷直流电机的建模与仿真2 3 电磁转矩测量模块根据电磁转矩方程(3)可知,电磁转矩与三相相电流、反电动势以及内外转子转速有关,所以在S i m ulink 中建立转矩模块时,模块输入信号为三相相电流、反电动势以及转速,通过加法和乘法计算,即可得到电磁转矩,如图4所示。
图4 电磁转矩计算模块2 4 转速计算模块根据式(5),即可搭建出如图5所示的转速模块。
电磁转矩和负载转矩经过加法、乘法和积分运算后,得到转子角速度信号,转子角速度乘以30/ 即为转速信号;对转子角速度进行积分运算即可得到转子位置信号。
图5 转速计算模块电压平衡方程模块、反电势计算模块、转矩计算模块和转速计算模块集合在一起就构成了永磁无刷直流电机的本体模型[7]。
3 实例仿真3 1 仿真对象为了验证模型的有效性,本文在M atlab /Si m u -li n k 中建立了数字控制系统的仿真模型,如图6所示。
仿真实验参数如下:电机工作状态为三相六状态,极对数p =5,电压U =270V,转速1000r/m in ,相电阻R =0 464 ,相电感L -M =0 0015H,电势系数k e =0 6V /rad s ,两转子转动惯量分别为J 1=0 01kg m 2、J 2=0 015kg m 2(不计螺旋桨转动惯量)。
摩擦转矩和负载转矩见表1:图6 实例电机仿真表1 三种负载情况对比表摩擦转矩Nm螺旋桨转矩系数Nm s/radT 01T 02B 1B 2图号111 50E-031 50E -03图7(a)112 00E-031 50E -03图7(b)20 51 50E-031 50E -03图7(c)3 2 仿真结果分析根据表1中负载数据,对模型进行仿真,结果如图7所示:2144卷图7 仿真结果曲线图7(a)-(c)是对应于表1中三种负载下的转速曲线,可以看出对转式永磁无刷直流电机起动过程中,内转子转速n 1上升快,外转子转速n2上升慢,这是由于两个转子的转动惯量不同,导致转子转速的变化量不同。