学院: 电气工程与信息工程学院
专业: 电气工程
姓名。学号。
指导教师: 教授
摘要:根据异步电机在两相同步旋转坐标系下的数学模型, 利用simulink软件包中的s函数模块和基本模块, 构建出电机的**模型, 并通过实例验证, 该建模方法求解精度高、简洁实用。
关键词:异步电动机;s函数;**;matlab
0引言。异步电机的动态模型是高阶、非线性、强耦合的多变量系统,通过坐标变换的方法对其进行简化后,模型简单得多,但其非线性、多变量的本质并未改变。描述电机的仍是一组高阶、变系数的微分方程,用传统的方法对其进行**分析并非易事。
将电机在两相静止坐标系下的数学模型做了一定的变换后,利用matlab/simulink**软件成功构建了电机的**模型。为使模型的建立更加简洁、明了,本文根据异步电机在两相同步旋转坐标系下的数学模型,充分利用matlab强大的矩阵运算功能及simulink提供的s函数模块,建立了普通异步电动机的**模型,并对实际电机进行了**。
1 电机两相同步旋转坐标系下的数学模型。
1.1 _ 三相静止到两相同步旋转坐标系的变换。
由三相静止坐标系变换到m -t坐标系要先将三相静止电流先变换到两相静止电流再将其变换到m -t坐标系下。
以下同坐标系下产生的磁动势完全一致为原则,两相同步旋转坐标系上的直流绕组可以等效于静止的三相对称交流绕组。其坐标变换矩阵为:
1 为同步角速度;ω2 为转子角速度;ω3为旋转坐标轴相对转子的角速度。
1.2 电压方程。
异步电机在两相同步旋转坐标系下的电压矩阵方程为:
式中,下标s、r表示定、转子分量;d、q表示同步坐标系的坐标轴;usd,usq;urd,urq分别为定、转子在d、q轴上的电压分量;isd、isq、ird、irq分别为定、转子在d、q轴上的电流分量;rs、rr分别为定、转子电阻;ls,lr分别为定、转子自感;lm为定、转子绕组间互感;p为微分算子。
1.3 转矩与运动方程。
转矩方程:运动方程;
式中,np为极对数;j为转动惯量;tl为负载转矩。
2 _ 电机**模块。
2.1 _ **模块。
matlab具有强大的矩阵运算能力和简便的绘图功能,simulink是matlab中用来对动态系统进行建模、**和综合分析的集成软件包。本文用simulink**软件提供的基本模块实现电机的坐标转换方程和转矩与运动方程,用s函数实现电机的电压方程,从而构建电机**图,如图1、图2、图3所示。
图1 电机**图。
图2 坐标转换图。
图3 转矩与运动模块。
2.2 电压方程的s函数描述。
s函数是用户借以自建simulink模块所必需的、具有特殊调用格式的函数文件。matlab提供了s函数m文件的标准模板,用户可以根据需要对其进行适当地改编,生成自己的s函数。s函数一旦被正确地嵌入simulink标准模块库中的s-function框架模块中,就可以像其他标准模块一样,与simulink的方程解算器solver交互,实现其功能。
由电压方程确定s函数模块的输入u(1),u(2),u(3)分别为坐标变换后的定子电压usd、usq及转子转速_;输出x(1),x(2),x(3),x(4)分别为定、转子d、q轴电流isd、isq、ird、irq。下面给出**模型中调用的v-equations_m函数中的主要程序。
定义用户s函数。
function [sys,x0,str,ts]=equations(t,x,u,flag,ls,lr,lm,rr,rs,w1)
switch flag,case 0,[sys,x0,str,ts]=mdlinitializesizes;
case 1,sys= mdlderivatives(t,x,u,ls,lr,lm,rr,rs,w1);
case 2,sys=mdlupdate(t,x,u);
case 3,sys= mdloutputs(t,x,u);
case 9,sys=mdlterminate(t,x,u);
otherwise
error(['unhandled flag = num2str(flag)])
end模块初始化子程序。
function [sys,x0,str,ts]=mdlinitializesizes
sizes = simsizes;
sizes. numcontstates = 4;
sizes. numdiscstates = 0;
sizes. numoutputs = 4;
sizes. numinputs = 3;
sizes. dirfeedthrough = 1;
sizes. numsampletimes = 1;
sys = simsizes(sizes);
x0 = 0; 0; 0; 0] ;
str =
ts = 0,0];
计算模块导数子程序。
function sys=mdlderivatives(t,x,u,ls,lr,lm,rr,rs,w1);
a= [ls,0,lm,0 ; 0,ls,0,lm ; lm,0,lr,0 ; 0,lm,0,lr] ;
b= [u(1) ;u(2) ;0 ; 0] ;
c= [rs,-w1*ls,0,-w1*lm ; w1*ls,rs,w1*lm,0 ; 0,(u(3)-w1)*lm,rr,(u(3)-w1)*lr;
(w1-u(3))*lm,0,(w1- u(3))*lr,rr] ;
x= [x(1);x(2);x(3);x(4)];
dx= inv(a)*(b-c*x);
sys=dx;
计算模块输出子程序。
function sys=mdloutputs(t,x,u)
sys=[x(1);x(2);x(3);x(4)];
function sys=mdlupdate(t,x,u)
sys=function sys=mdlterminate(t,x,u)
sys=3 电机**实验。
以1台三相四极星形联结的感应电动机为例,电机参数为:额定功率pn=1.1kw;极对数np=2;转动惯量j=0.
008kg·m2;定子电阻rs=3.433_;转子电阻rr=5.533_;定子电感ls=0.
203h;转子电感lr=0.207h;定、转子互感lm=0.193h。
给电机直接加入380v、50hz额定电压,让其空载起动,并在0.4s时加入3n·m的负载转矩,**结果如图4、图5所示。由输出波形可以看出,模型成功地**了电机的起动及加载运行过程。
图4 转子转速**曲线。
图5 电磁转矩**曲线。
4结论。应用matlab/simulink**时,可以利用simulink提供的各种标准模块,根据信号流程来搭建电机模型。这种方法方便、直观,但在求解高阶非线性微分方程组时,模型的搭建将会非常繁琐且容易出错。
s函数是matlab中对一般动态系统的计算语言描述,可以用matlab的m文件或c语言编写。通过灵活地编写s函数m文件或c语言编写。通过灵活地编写s函数m文件的标准模板,能够准确、快捷地实现高阶非线性微分方程组的求解。
本文介绍了异步感应电机在两相同步旋转坐标系下数学模型的建立,并结合实例介绍了一种应用simulink进行电机实时**的方法。该建模方法充分利用s函数,并灵活结合基本模块的搭建,使得**模型简洁易懂、通用性强,可以推广至其他类型电机的**。
地图数学模型原理结课作业
一 在空间数据库中,把大比例尺图形数据缩编成小比例尺,图形数据是按要素分层的,各要素应采用什么模型确定选取指标?答 地图缩编通常会将图形数据按要素分层,如居民地 河流 道路等图层,各种要素选取指标的确定通常又有多种模型可以选择,下面采用一元回归数学模型说明居民地和河流的指标选取过程。1.居民地选取指...
地图数学模型原理结课作业
武汉大学研究生结课作业。地图数学模型原理与分析。院 系 名称 资源与环境科学学院。专业名称 土地资源管理。学生姓名 秦坤。学生学号 2014202050103 指导教师 何宗宜。二 一五年六月。一 在空间数据库中,把大比例尺图形数据缩编成小比例尺,图形数据是按要素分层的,各要素应采用什么模型确定选取...
研究生课程紊流数学模型结课作业
解 时均方程的定义 紊流运动的实验研究表明,虽然紊流结构十分复杂,但它仍然遵循连续介质的一般动力学规律,因此,雷诺在1886年提出用时均值概念来研究紊流运动。他认为,紊流中任何物理量虽然都随时间和空间而变化,但是任一瞬时的运动仍然符合连续介质流动的特征,流场中任一空间点上应该适用黏性流体运动的基本方...