异步电动机动态数学模型仿真文档格式.docx
- 文档编号:16607610
- 上传时间:2022-11-24
- 格式:DOCX
- 页数:29
- 大小:1.22MB
异步电动机动态数学模型仿真文档格式.docx
《异步电动机动态数学模型仿真文档格式.docx》由会员分享,可在线阅读,更多相关《异步电动机动态数学模型仿真文档格式.docx(29页珍藏版)》请在冰豆网上搜索。
图1-1三相异步电动机的物理模型
在图1-1中,定子三相绕组轴线A、B、C在空间是固定的,以A轴为参考坐标轴;
转子绕组轴线a、b、c随转子旋转,转子a轴和定子A轴间的电角度θ为空间角位移变量。
规定各绕组电压、电流、磁链的正方向符合电动机惯例和右手螺旋定则。
1.1.1异步电动机三相动态模型的数学表达式
异步电机的动态数学模型由下述电压方程、磁链方程、转矩方程和运动方程组成。
其中磁链方程和转矩方程为代数方程,电压方程和运动方程为微分方程。
(1)磁链方程
异步电动机每个绕组的磁链是它本身的自感磁链和其它绕组对它的互感磁链之和,因此,六个绕组的磁链可用下式表示:
(1-1)
或写成(1-1a)
其中L为6
6电感矩阵,其中对角线元素其中对角线元素
、
是各有关绕组的自感,其余各项则是绕组间的互感。
绕组之间的互感又分为两类:
①定子三相彼此之间和转子三相彼此之间位置都是固定的,故互感为常值;
②定子任一相与转子任一相之间的相对位置是变化的,互感是角位移的函数。
(2)电压方程
三相定子的电压方程可表示为:
(1-2)
方程中,
为定子三相电压;
为定子三相电流;
为定子三相绕组磁链;
为定子各相绕组电阻。
三相转子绕组折算到定子侧后的电压方程为:
(1-3)
为转子三相电压;
为转子三相电流;
为转子三相绕组磁链;
为转子各相绕组电阻。
将电压方程写成矩阵形式:
(1-4)
或写成(1-4a)
(3)转矩方程
电磁转矩方程为:
(1-5)
式中,
为电机极对数,
为角位移。
(4)运动方程
运动方程为:
(1-6)
为电磁转矩;
为负载转矩;
为电机机械角速度;
为转动惯量。
1.1.2异步电动机三相原始模型的性质
由异步电动机的三相动态模型可见,非线性耦合早电压方程、磁链方程与转矩方程中都有体现。
既存在定子和转子间的耦合,也存在三相绕组间的交叉耦合。
旋转电动势和电磁转矩中都包含变量之间的乘积,这是非线性的基本因素。
定转子间的相对运动,导致其夹角不断变化,使得互感矩阵为非线性变参数矩阵。
所有这些,都使异步电动机成为高阶、非线性、强耦合的多变量系统。
假定异步电动机三相绕组为Y无中线连接,若为Δ连接,可等效为Y连接。
可以证明:
异步电动机三相数学模型中存在一定的约束条件:
(1-7)
同理,转子绕组也存在相应的约束条件:
(1-8)
以上分析表明,对于无中性线Y/Y联结绕组的电动机,三相变量中只有两相是独立的,因此三相原始数学模型并不是物理对象最简洁的描述,完全可以而且也有必要用两项模型代替。
1.2坐标变换
异步电动机三相原始动态模型相当复杂,分析和求解这组非线性方程十分困难。
在实际中必须予以简化,简化的基本方法就是坐标变换。
异步电动机数学模型之所以复杂,关键是因为有一个复杂的电感矩阵和转矩方程,它们体现了异步电动机的电磁耦合和能量转换的复杂关系。
因此,要简化数学模型,须从电磁耦合关系入手。
1.2.1坐标变换的基本思路
如果能将交流电动机的物理模型等效地变换成类似直流电动机的模式,分析和控制就可以大大简化。
坐标变换正是按照这条思路进行的。
不同坐标系中电动机模型等效的原则是:
在不同坐标下绕组所产生的合成磁动势相等。
三相变量中只有两相为独立变量,完全可以也应该消去一相。
所以,三相绕组可以用相互独立的两相正交对称绕组等效代替,等效的原则是产生的磁动势相等。
两相绕组,通以两相平衡交流电流,也能产生旋转磁动势。
当三相绕组和两相绕组产生的旋转磁动势大小和转速都相等时,即认为两相绕组与三相绕组等效,这就是3/2变换。
两个匝数相等相互正交的绕组d、q,分别通以直流电流,产生合成磁动势F,其位置相对于绕组来说是固定的。
如果人为地让包含两个绕组在内的铁心以同步转速旋转,磁动势F自然也随之旋转起来,成为旋转磁动势。
如果旋转磁动势的大小和转速与固定的交流绕组产生的旋转磁动势相等,那么这套旋转的直流绕组也就和前面两套固定的交流绕组都等效了。
1.2.2三相-两相变换(3/2变换)
三相绕组A、B、C和两相绕组之间的变换,称作三相坐标系和两相正交坐标系间的变换,简称3/2变换。
图1-2三相坐标系和两相正交坐标系中的磁动势矢量
ABC和两个坐标系中的磁动势矢量,将两个坐标系原点重合,并使A轴和
轴重合。
按照磁动势相等的等效原则,三相合成磁动势与两相合成磁动势相等,故两套绕组磁动势在αβ轴上的投影应相等,因此
(1-9)
写成矩阵形式,得:
(1-10)
按照变换前后总功率不变,匝数比为
(1-11)
则三相坐标系变换到两相正交坐标系的变换矩阵
(1-12)
两相正交坐标系变换到三相坐标系(简称2/3变换)的变换矩阵
(1-13)
1.2.3静止两相-旋转正交变换(2s/2r变换)
从静止两相正交坐标系αβ到旋转正交坐标系dq的变换,称作静止两相-旋转正交变换,简称2s/2r变换,其中s表示静止,r表示旋转,变换的原则同样是产生的磁动势相等。
图1-3静止两相正交坐标系和旋转正交坐标系中的磁动势矢量
旋转正交变换阵为
(1-14)
静止两相正交坐标系到旋转正交坐标系的变换阵
(1-15)
旋转正交坐标系到静止两相正交坐标系的变换阵
(1-16)
定子旋转变换阵
(1-17)
转子旋转变换阵
(1-18)
电压方程
(1-19)
磁链方程
(1-20)
转矩方程
(1-21)
旋转变换是用旋转的绕组代替原来静止的定子绕组,并使等效的转子绕组与等效的定子绕组重合,且保持严格同步,等效后定、转子绕组间不存在相对运动。
旋转正交坐标系中的磁链方程和转矩方程与静止两相正交坐标系中相同,仅下标发生变化。
从表面上看来,旋转正交坐标系中的数学模型还不如静止两相正交坐标系的简单,实际上旋转正交坐标系的优点在于增加了一个输入量ω1,提高了系统控制的自由度。
2异步电动机在正交坐标系上的动态数学模型
异步电动机三相原始模型相当复杂,通过坐标变换能够简化数学模型,便于进行分析和计算。
按照从特殊到一般,首先推导静止两相正交坐标系中的数学模型,然后推广到旋转正交坐标系。
由于运动方程不随坐标变换而变化,故仅讨论电压方程、磁链方程和转矩方程。
在以下论述中,下标s表示定子,下标r表示转子。
2.1静止两相正交坐标系中的动态数学模型
异步电动机定子绕组是静止的,只要进行3/2变换就行了。
转子绕组是旋转的,必须通过3/2变换和旋转到静止的变换,才能变换到静止两相正交坐标系。
(1)定子绕组和转子绕组的3/2变换
对静止的定子三相绕组和旋转的转子三相绕组进行相同的3/2变换,变换后的定子两相正交坐标系静止,而转子两相正交坐标系以角速度逆时针旋转,如图2-1a所示。
图2-1定子、转子坐标系到静止两相正交坐标系的变换
相应的数学模型如下:
(2-1)
磁链方程
(2-2)
转矩方程
(2-3)
3/2变换将按三相绕组等效为互相垂直的两相绕组,从而消除了定子三相绕组、转子三相绕组间的相互耦合。
但定子绕组与转子绕组间仍存在相对运动,因而定、转子绕组互感阵仍是非线性的变参数阵。
输出转矩仍是定、转子电流及其定、转子夹角的函数。
与三相原始模型相比,3/2变换减少了状态变量的维数,简化了定子和转子的自感矩阵。
(2)静止两相正交坐标系中的矩阵方程
对图2-1中的转子坐标系
作旋转变换(旋转正交坐标系到静止两相正交坐标系的变换),即将
坐标系顺时针旋转
角,使其与定子
坐标系重合,且保持静止,即用静止的两相转子正交绕组等效代替原先转动的两相绕组,如图2-1b所示。
旋转变换阵为
(2-4)
变换后的电压方程为
(2-5)
磁链方程为
(2-6)
转矩方程为
(2-7)
旋转变换改变了定、转子绕组间的耦合关系,将相对运动的定、转子绕组用相对静止的等效绕组来代替,消除了定、转子绕组间夹角对磁链和转矩的影响。
旋转变换的优点在于将非线性变参数的磁链方程转化为线性定常的方程,但却加剧了电压方程中非线性耦合程度,将矛盾从磁链方程转移到电压方程中来了,并没有改变对象的非线性耦合程度。
2.2旋转正交坐标系中的动态数学模型
更广义的坐标旋转是对定子坐标系和转子坐标系同时施行旋转变换,把它们变换到同一个旋转正交坐标系dq上,dq相对于定子的旋转角速度为
,如图2-2a所示。
图2-2定子、转子坐标系到旋转正交坐标系的变换
a)定子、转子坐标系b)旋转正交坐标系
定子旋转变换阵为
(2-8)
转子旋转变换阵为
(2-9)
电压方程为
(2-10)
(2-11)
(2-12)
旋转变换是用旋转的绕组代替原来静止的定子绕组,并使等效的转子绕组与等效的定子绕组重合,且保持严格同步,等效后定、转子绕组间不存在相对运动。
故旋转正交坐标系中的磁链方程和转矩方程与静止两相正交坐标系中相同,仅下标发生变化。
从表面上看来,旋转正交坐标系中的数学模型还不如静止两相正交坐标系的简单,实际上旋转正交坐标系的优点在于增加了一个输入量ω1,提高了系统控制的自由度,磁场定向控制就是通过选择ω1而实现的。
旋转速度任意的正交坐标系无实际使用意义,常用的是同步旋转坐标系,将绕组中的交流量变为直流量,以便模拟直流电动机进行控制。
2.3异步电动机在正交坐标系上的状态方程
旋转正交坐标系上的异步电动机具有4阶电压方程和1阶运动方程,因此须选取5个状态变量。
可选的状态变量共有9个,这9个变量分为5组:
①转速;
②定子电流;
③转子电流;
④定子磁链;
⑤转子磁链。
转速作为输出变量必须选取。
其余的4组变量可以任意选取两组,定子电流可以直接检测,应当选为状态变量。
剩下的3组均不可直接检测或检测十分困难,考虑到磁链对电动机的运行很重要,可以选定子磁链或转子磁链。
状态方程为状态变量。
状态变量
输入变量
输出变量
状态方程
输出方程
转子电磁时间常数
电动机漏磁系数
根据以上公式绘制动态结构图如图:
图2-3为状态变量在dq坐标系中动态结构图
3异步电动机模型仿真
3.1ACMotor模块
根据图2-3的动态结构图,用MATLAB/SIMULINK基本模块建立在dq坐标系下异步电动机仿真模型ACMotor模块。
ACMotor模块图如图3-1。
根据图2-3计算参数为:
图3-1ACMotor模块图
3.2坐标变换模块
(1)3/2transform模块
根据静止两相正交坐标系到旋转正交坐标系的变换阵
则有Usa=0.8165*Ua-0.4082*Ub-0.4082*Uc,Usb=0.7071*Ub-0.7071*Uc
其中Ua,Ub,Uc为三相坐标系下的输入电压,Usa和Usb为静止两相正交坐标下的电压。
搭建模块如下图:
图3-23/2transform模块(a)
图3-33/2transform模块(b)
(2)2s/2rtransform模块
根据定子旋转变换阵
则有Usd=cos
Usa+sin
Usb,Usq=-sin
Usa+
Usb
其中Usa和Usb为静止两相正交坐标下的电压,Usd和Usq为两相旋转坐标系下的电压。
为d轴与a轴的夹角。
搭建模块如下图
图3-42s/2rtransform模块(a)
图3-52s/2rtransform模块(b)
(3)2r/2stransform模块
根据旋转正交坐标系到静止两相正交坐标系的变换阵
则有Isa=cos
Isd-sin
Isq,Isb=sin
Isd+cos
Isq
其中Isa和Isb为静止两相正交坐标下的电压,Isd和Isq为两相旋转坐标系下的电压。
图3-62r/2stransform模块(a)
图3-72r/2stransform模块(b)
4)2/3transform模块
两相正交坐标系变换到三相坐标系(简称2/3变换)的变换矩阵
则有Ia=0.8165Isa,Ib=-0.4082Isa+0.7071Isb,Ic=-0.4082Isa-0.7071Isb
其中Ia,Ib,Ic为三相坐标系下的输入电流,Isa和Isb为静止两相正交坐标下的电流。
搭建模块如下图:
图3-82/3transform模块(a)
图3-92/3transform模块(b)
(5)2r/2stransform模块和2s/2rtransform模块中输入量Theta根据下式得到。
(6)若由三相坐标系直接变换到两相旋转坐标系下,得到其坐标变换矩阵为:
搭建仿真模型为
图3-103/2rtransform模块
3.3仿真原理图
在进行异步电动机仿真时,以为状态变量的dq坐标系中的状态方程为内核,在外围加上坐标变换和状态变换,就可得到在dq坐标系下的仿真结果。
仿真原理图如下图所示。
图3-11仿真原理图(a)
图3-12仿真原理图(b)
其中有5个输入参数:
三相正弦交流电压Usa,Usb,Usc,同步转速W1,负载转矩T1。
三相正弦交流电压幅值均为380V,频率为100*piHZ,相角分别为0、-2*pi/3、2*pi/3,同步转速为常数100*pi,根据
所以额定输出转矩为19.76N.M,设定负载转矩为阶跃信号,阶跃时间为1s,阶跃初始值为0,终值为19.76N.M。
正弦电压源Usa参数设置如图3-13所示,且a、b、c相角互差
。
图3-13正弦电压源Usa参数设置
4仿真结果及分析
当负载转矩为额定输出转矩19.76N.M时,电磁转矩和转速结果图如下:
图4-1转矩与转速结果图
图4-2转矩与转速局部结果图
当负载转矩为额定转矩时,转速不为额定转速。
为使转速接近于额定值,调整负载转矩为15.4N.M。
此时,电磁转矩和转速结果图如下:
图4-3转矩与转速结果图
图4-4转矩与转速局部结果图
由图4-3和图4-4可知,电动机空载启动时,转速迅速上升并达到稳定值1500r/min,电磁转矩在转速上升时作衰减震荡,最后稳定值为零。
在1s时突加负载Tl=15.4N.M,转速降至1450r/min,即额定值。
三相电流结果图如下:
图4-5三相电流结果图
图4-6空载稳定三相电流
图4-7带额定负载时的三相电流
由图4-6和图4-7可知,空载运行电流幅值为4A,额定转速时运行电流幅值7A。
结论
本文详细地介绍了基于Matlab/Simulink软件下,建立异步电动机定子磁链仿真模型的原理、参数计算及动态模型框图。
三相异步电动机本身是一个高阶、非线性、强耦合的多变量系统,传统的分析方法已很难适应这样复杂系统的分析,计算机仿真技术的发展为复杂系统的分析提供了极为有利的条件。
在分析异步电动机的物理模型后,建立异步电动机的动态数学模型,然后推导出两相旋转坐标系上的状态方程和转矩方程,利用Matlab/Simulink仿真工具把数学方程转变为模型。
运行异步电动机的仿真模型,可观察到异步电动机在启动和加载的情况下,转速、电磁转矩、定子磁链和定子电流的变化曲线,同时分析各个变量之间的变化关系。
进一步了解异步电动机的运行特性。
仿真结果表明,用Simulink进行三相异步电动机仿真比较方便,且高效直观,得到的结果也是比较接近实际。
但是,在仿真的过程中,也遇到了一些问题。
比如,在设置仿真参数时,参数不合理,无法观察到正确的波形。
因此,在设置参数时,需要一定的技巧才能快速地得到满意的仿真波形。
还有就是异步电动机的参数很关键,其精确度关系到构建的异步电动机模型是否符合实际。
在建立异步电动机的状态方程时,采用的一些近似处理对模型仿真结果也有一定的影响。
因此,应尽可能得到异步电动机的精确参数来构造模块,这样针对性更强,仿真精度更高,仿真结果更可靠。
参考文献
[1]洪乃刚.电力电子和电力拖动控制系统的Matlab仿真.北京:
机械工业出版社,2006
[2]陈伯时.电力拖动自动控制系统(第三版).北京:
机械工业出版社,2003
[3]薛定宇.基于MATLAB/Simulink的系统仿真技术与应用.北京:
清华大学出版社,2002
[4]陈桂明.张明照.应用MATLAB建模与仿真[P].北京:
科学出版社,2001
[5]李夙.异步电动机直接转距控制[M].北京:
机械工业出版社,1999
[6]彭鸿才.电机原理及拖动.北京:
机械工业出版社,1996
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 异步电动机 动态 数学模型 仿真