无刷直流电机数学模型完整版

1、电机数学模型以二相导通星形三相六状态为例，分析BLDC 的数学模型及电磁转矩等特性。为了便于分析，假定：a)三相绕组完全对称，气隙磁场为方波，定子电流、转子磁场分布皆对称；b)忽略齿槽、换相过程和电枢反应等的影响；c)电枢绕组在定子内表面均匀连续分布；d)磁路不饱和，不计涡流和磁滞损耗。则三相绕组的电压平衡方程可表示为：(1)式中：为定子相绕组电压(V) ；为定子相绕组电流(A) ；为定子相绕组电动势 (V) ；L 为每相绕组的自感 (H) ；M 为每相绕组间的互感 (H)；p 为微分算子 p=d/dt。三相绕组为星形连接，且没有中线，则有(2)(3)得到最终电压方程：(4)rL-Miaeai

2、bicr L-MrL-Mebec图 .无刷直流电机的等效电路无刷直流电机的电磁转矩方程与普通直流电动机相似， 其电磁转矩大小与磁通和电流幅值成正比(5)所以控制逆变器输出方波电流的幅值即可以控制 BLDC 电机的转矩。为产生恒定的电磁转矩， 要求定子电流为方波， 反电动势为梯形波， 且在每半个周期内，方波电流的持续时间为 120°电角度，梯形波反电动势的平顶部分也为 120° 电角度，两者应严格同步。由于在任何时刻，定子只有两相导通，则：电磁功率可表示为：(6)电磁转矩又可表示为：(7)无刷直流电机的运动方程为：(8)其中为电磁转矩；为负载转矩； B 为阻尼系数；为电机机械

3、转速；J为电机的转动惯量。传递函数：无刷直流电机的运行特性和传统直流电机基本相同， 其动态结构图可以采用直流电机通用的动态结构图，如图所示：+I(s)+U(s)TC (s)N(s)1/RCt365/(GD2s)-T L (s)Ce图2. 无刷直流电机动态结构图由无刷直流电机动态结构图可求得其传递函数为:式中：K1为电动势传递系数，Ce 为电动势系数；K2为转矩传递函数，R 为电动机内阻， Ct 为转矩系数；Tm 为电机时间常数，G 为转子重量， D 为转子直径。基于 MATLAB 的 BLDC 系统模型的建立在 Matlab 中进行 BLDC 建模仿真方法的研究已受到广泛关注，已有提出采用节点

4、电流法对电机控制系统进行分析，通过列写m 文件，建立 BLDC 仿真模型，这种方法实质上是一种整体分析法，因而这一模型基础上修改控制算法或添加、删除闭环就显得很不方便；为了克服这一不足，提出在Matlab/Simulink 中构造独立的功能模块，通过模块组合进行 BLDC 建模，这一方法可观性好，在原有建模的基础上添加、 删除闭环或改变控制策略都十分便捷， 但该方法采用快速傅立叶变换 (FFT)方法求取反电动势，使得仿真速度受限制。本文提出了一种新型的 BLDC 建模方法，将控制单元模块化，在 Matlab/Simulink 建立独立的功能模块：BLDC 本体模块、电流滞环控制模块、速度控制模

5、块、参考电流模块、转矩计算模块和电压逆变模块， 对这些功能模块进行有机整合， 即可搭建出无刷直流电机系统的仿真模型。 在建模过程中， 梯形波反电动势的求取方法一直是较难解决的问题 27,28，本文采用分段线性法成功地化解了这一难点， 克服了建模方法存在的不足。Matlab6.5 针对电气传动控制领域所设计的工具箱SimPowerSystemToolbox2.3已提供了 PMSM 的电机模型，但没有给出BLDC 的电机模型。因此，本文在分析无刷直流电机数学模型的基础上，借助于Matlab 强大的仿真建模能力，在Matlab/Simulink 中建立了 BLDC 控制系统的仿真模型。BLDC 建模

6、仿真系统采用双闭环控制方案：下即为 BLDC 建模的整体控制框图， 其中主要包括： BLDC 本体模块、电流滞环控制模块、速度控制模块、参考电流模块、转矩计算模块和电压逆变模块。BLDC 本体结构(1) BLDCM 本体模块在整个控制系统的仿真模型中，BLDCM 本体模块是最重要的部分，该模块根据BLDC 电压方程式 (4)求取 BLDC 三相相电流，结构框图如图所示图 .BLDCM本体模块结构框图及其封装形式在整个控制系统的仿真模型中，BLDC 本体模块是最重要的部分，该模块根据BLDC 电压方程式 (2-4)求取 BLDC 三相相电流，而要获得三相相电流信号ia，ib，ic，必需首先求得三

7、相反电动势信号 ea，eb，ec 控制框图如图 2-11 所示。而 BLDC 建模过程中，梯形波反电动势的求取方法一直是较难解决的问题， 反电动势波形不理想会造成转矩脉动增大、 相电流波形不理想等问题， 严重时会导致换相失败，电机失控。因此，获得理想的反电动势波形是 BLDC 仿真建模的关键问题之一。本文采用了分段线性法，如图 2-12 所示，将一个运行周期 0° 360°分为 6 个阶段，每 60°为一个换相阶段，每一相的每一个运行阶段都可用一段直线进行表示，根据某一时刻的转子位置和转速信号，确定该时刻各相所处的运行状态，通过直线方程即可求得反电动势波形。 分段

8、线性法简单易行， 且精度较高， 能够较好的满足建模仿真的设计要求。 因而，本文采用分段线性法建立梯形波反电动势波形。理想情况下，二相导通星形三相六状态的BLDC定子三相反电动势的波形如图2-12 所示。图中，根据转子位置将运行周期分为6 个阶段： 0 /3， /32 /3，2 /3，4 /3， 4 /35 /3，5 /32。以第一阶段 0 /3为例， A 相反电动势处于正向最大值 Em， B 相反电动势处于负向最大值 -Em，C 相反电动势处于换相阶段，由正的最大值 Em 沿斜线规律变化到负的最大值 -Em。根据转子位置和转速信号，就可以求出各相反电动势变化轨迹的直线方程，其它5 个阶段，也是

9、如此。据此规律，可以推得转子位置和反电动势之间的线性关系，如表2-1所示，从而采用分段线性法，解决了在BLDC 本体模块中梯形波反电动势的求取问题。Emea-EmposEmebpos-EmEmecpos-Em0/32/34/35/32图 . 三相反电动势波形转子位置和反电动势之间的线性关系表转子位置eabcee0/3K*w- K*wK*w*(per-pos)/( /6)+1)K*w*(pos- /6/32/3K*w- K*w-per)/( -1)/6)K*w*(per+2* /32/3K*w- K*w-pos)/( /6)+1)4/3- K*wK*wK*w*(pos- -per)/( /6)-

10、1)4/35/3- K*wK*w*(per+4* /3K*w-pos)/( /6)+1)5/32K*w*(pos- 5* /3K*w- K*w-per)/( -1)/6)表中： K 为反电动势系数(V/(r/min) ， pos 为角度信号， w 为转速信号，转数per=fix(pos/(2*pi)*2*pi ， fix 函数是实现取整功能。根据上式，用 M 文件编写反电势系数的S 函数如下：反电动势S 函数 (emf.m)%= %BLDCM模型中反电动势函数%= function sys,x0,str,ts =emf(t,x,u,flag) switch flagcase 0,%初始化设置s

11、ys,x0,str,ts=mdlInitializeSizes;case 3,%输出量计算sys = mdlOutputs(t,x,u);case 1,2,4,9%未定义标志sys = ;otherwise%错误处理error('unhandled flag = ',num2str(flag);end%=%mdlInitializeSizes进行初始化，设置系统变量的大小%= function sys,x0,str,ts=mdlInitializeSizes()sizes = simsizes;%取系统默认设置sizes.NumContStates= 0;sizes.NumDi

12、scStates= 0;sizes.NumOutputs= 3;sizes.NumInputs= 2;sizes.DirFeedthrough= 1;sizes.NumSampleTimes= 1;sys = simsizes(sizes);x0 = ;str = ;ts = -1 0;%=%mdlOutputs计算系统输出%= function sys=mdlOutputs(t,x,u)globalk;globalPos;globalw;k=0.060;% V/(r/min)反电动势系数w=u(1);% 转速 (rad/s)Pos=u(2);% 角度 (rad)ifPos>=0 &am

13、p; Pos<=pi/3sys=k*w,-k*w,k*w*(-Pos)/(pi/6)+1);elseifPos>=pi/3 & Pos<=2*pi/3sys=k*w,k*w*(Pos-pi/3)/(pi/6)-1),-k*w;elseifPos>=2*pi/3 & Pos<=pisys=k*w*(2*pi/3-Pos)/(pi/6)+1),k*w,-k*w;elseifPos>=pi & Pos<=4*pi/3sys=-k*w,k*w,k*w*(Pos-pi)/(pi/6)-1);elseifPos>=4*pi/3 &am

14、p; Pos<=5*pi/3sys=-k*w,k*w*(4*pi/3-Pos)/(pi/6)+1),k*w;elsePos>=5*pi/3 & Pos<=2*pisys=k*w*(Pos-5*pi/3)/(pi/6)-1),-k*w,k*w;end转矩计算模块根据 BLDC 数学模型中的电磁转矩方程式， 可以建立图 5.7 所示的转矩计算模块，模块输入为三相相电流与三相反电动势，通过加、乘模块即可求得电磁转矩信号Te 。转矩计算模块结构框图及其封装形式转速计算模块根据运动方程式 (2.4)，由电磁转矩、负载转矩以及摩擦转矩，通过加乘、积分环节即可得到转速信号，求得的转

15、速信号经过积分就可得到电机转角信号，如图转速计算模块结构框图及其封装形式电流滞环控制模块在这个仿真模块中采用滞环控制原理来实现电流的调节， 使得实际电流随跟定电流的变化。模块结构框图如图 5.10 所示 40 ，输入为三相参考电流和三相实际电流，输出为 PWM 逆变器控制信号。电流滞环控制模块结构框图及其封装参考电流模块参考电流模块的作用是根据电流幅值信号 Is 和位置信号给出三相参考电流，输出的三相参考电流直接输入电流滞环控制模块， 用于与实际电流比较进行电流滞环控制。转子位置和三相参考电流之间的对应关系如表所示， 参考电流模块的这一功能可通过 S 函数编程实现，程序如下参考电流 S 函数

16、(mod.m)function sys,x0,str,ts =mod(t,x,u,flag)switch flagcase 0,sys,x0,str,ts=mdlInitializeSizes;case 3,sys = mdlOutputs(t,x,u);case 2,sys = ;case 9,sys = ;otherwiseerror('unhandled flag = ',num2str(flag);endfunction sys,x0,str,ts=mdlInitializeSizes()sizes = simsizes;sizes.NumContStates= 0;s

17、izes.NumDiscStates= 0;sizes.NumOutputs= 1;sizes.NumInputs= 1;sizes.DirFeedthrough = 1;sizes.NumSampleTimes = 1;sys = simsizes(sizes);x0= ;str = ;ts= -1 0;function sys=mdlOutputs(t,x,u)global Pos;global w;global Theta;Theta=u;b=fix(Theta/(2*pi);% 取整if Theta=0sys=0;else if (Theta/(2*pi)=bsys=2*pi;else

18、sys=Theta-b*2*pi;endendPos=sys; %位置表 5.2 转子位置和三相参考电流之间的对应关系表转子位置 I_ar I_br I_cr0/3Is - Is0/32/3Is 0 - Is 2/3 0 Is - Is4/3- Is Is04/35/3 - Is0Is5/32 0 - Is Is位置计算模块电机转角信号到电机位置信号的转换可通过S 函数编程实现，程序如下位置计算S 函数 (is.m)function sys,x0,str,ts =is(t,x,u,flag)switch flagcase 0,sys,x0,str,ts=mdlInitializeSizes;c

19、ase 3,sys = mdlOutputs(t,x,u);case 2,sys = ;case 9,sys = ;otherwiseerror('unhandled flag = ',num2str(flag);endfunction sys,x0,str,ts=mdlInitializeSizes()sizes = simsizes;sizes.NumContStates sizes.NumDiscStates= 0;= 0;sizes.NumOutputs= 3;sizes.NumInputssizes.DirFeedthrough = 1;sizes.NumSample

20、Times = 1;sys = simsizes(sizes);x0= ;str = ;ts= -1 0;= 2;function sys=mdlOutputs(t,x,u)global Is;global Pos;Is=u(1); %电流Pos=u(2);%位置if Pos>=0& Pos<=pi/3sys=Is,-Is,0;elseif Pos>=pi/3& Pos<=2*pi/3sys=Is,0,-Is;elseif Pos>=2*pi/3& Pos<=pisys=0,Is,-Is;elseif Pos>=pi& Pos<=4*pi/3sys=-Is,Is,0;elseif Pos>=4*pi/3& Pos<=5