logo资料库

基于MATLAB的异步电动机仿真.doc

第1页 / 共23页
第2页 / 共23页
第3页 / 共23页
第4页 / 共23页
第5页 / 共23页
第6页 / 共23页
第7页 / 共23页
第8页 / 共23页
资料共23页,剩余部分请下载后查看
目 录
1 引言
2 异步电动机动态数学模型
2.1异步电动机动态数学模型的性质
2.2三相异步电动机的多变量非线性数学模型
2.2.1 电压方程
2.2.2 磁链方程
2.2.3 转矩方程
2.2.4 电力拖动系统运动方程
2.2.5 三相异步电机的数学模型
3 坐标变化和变换矩阵
3.1三相--两相变换(3/2变换)
3.2三相异步电动机在两相坐标系上的数学模型
3.2.1三相异步电动机在两相坐标系上的状态方程
3.2.2两相静止坐标系中按定子磁链定向的状态方程
4 软件介绍及模型实现
4.1 Matlab/Simulink简介
4.2模型实现
4.2.1 Simulink模型设计
4.2.2模型参数设置
4.2.3仿真结果
5 结论
参考文献
目 录 1 引言................................................................................................................................ 1 2 异步电动机动态数学模型............................................................................................ 2 2.1 异步电动机动态数学模型的性质....................................................................... 2 2.2 三相异步电动机的多变量非线性数学模型....................................................... 2 2.2.1 电压方程..................................................................................................... 3 2.2.2 磁链方程..................................................................................................... 4 2.2.3 转矩方程..................................................................................................... 6 2.2.4 电力拖动系统运动方程............................................................................. 7 2.2.5 三相异步电机的数学模型......................................................................... 8 3 坐标变化和变换矩阵.................................................................................................... 9 3.1 三相--两相变换(3/2 变换)...............................................................................9 3.2 三相异步电动机在两相坐标系上的数学模型................................................. 10 3.2.1 三相异步电动机在两相坐标系上的状态方程........................................11 3.2.2 两相静止坐标系中按定子磁链定向的状态方程....................................11 4 软件介绍及模型实现.................................................................................................. 13 4.1 Matlab/Simulink 简介......................................................................................... 13 4.2 模型实现............................................................................................................. 13 4.2.1 Simulink 模型设计....................................................................................13 4.2.2 模型参数设置........................................................................................... 15 4.2.3 仿真结果................................................................................................... 18 5 结论.............................................................................................................................. 21 参考文献.......................................................................................................................... 22
课程设计说明书 1 引言 1985 年,由 Depenbrock 教授提出的直接转距控制理论将运动控制的发展向前 推进了一大步。接着 1987 年把它又推广到弱磁调速范围。不同于矢量控制技术, 它无需将交流电动机与直流电动机作比较、等效和转化,不需要模仿直流电动机的 控制,也不需要为解耦而简化交流电动机的数学模型[1]。它只是在定子坐标系下分 析交流电机的数学模型,强调对电机的转距进行直接控制,省掉了矢量旋转变换等 复杂的变换与计算。直接转距控制从一诞生,就以新颖的控制思想,简洁明了的系 统结构,优良的静、动态性能受到人们的普遍关注。 系统建模与仿真一直是各领域研究、分析和设计各种复杂系统的有力工具。建 模可以超越理想的去模拟复杂的现实物理系统;而仿真则可以对照比较各种控制策 略和方案,优化并确定系统参数。长期以来,仿真领域的研究重点是放在仿真模型 建立这一环节上,即在系统模型建立以后,设计一种算法,以使系统模型为计算机 所接受,然后再将其编制成计算机程序,并在计算机上运行。显然,为达到理想的 目的,在这一过程中编制与修改仿真程序十分耗费时间和精力,这也大大阻碍了仿 真技术的发展和应用。近年来逐渐被大家认识的 Matlab 语言则很好的解决了这个问 题。 1
课程设计说明书 2 异步电动机动态数学模型 2.1 异步电动机动态数学模型的性质 直流电动机的磁通由励磁绕组产生,可以在电枢合上电源以前建立起来而不参 与系统的动态。 过程(弱磁调速时除外)。因此,它的动态数学模型只有一个输入变量——电 枢电压和一个输出变量——转速,在控制对象中含有机电时间常数 m 和电枢回路电 磁时间常数 l ,如果电力电子变换装置也计入控制对象,则还有滞后的时间常数 s 。 在工程上能够允许的一些假定条件下,可以描述成单变量(单输入单输出)的三阶 线性系统[2],完全可以应用经典的线性控制理论和由它发展出来的工程设计方法进 行分析与设计。 但是,同样的理论和方法用来分析与设计交流调速系统时,就不那么方便了, 因为交流电机的数学模型和直流电机模型相比有着本质上的区别。 1)异步电机变压变频调速时需要进行电压(或电流)和频率的协调控制,有 电压(电流)和频率两种独立的输入变量。在输出变量中,除转速外,磁通也得算 一个独立的输出变量。因为电机只有一个三相输入电源,磁通的建立和转速的变化 是同时进行的,为了获得良好的动态性能,也希望对磁通施加某种控制,使它在动 态过程中尽量保持恒定,才能产生较大的动态转矩。 由于这些原因,异步电机是一个多变量(多输入多输出)系统,而电压(电流)、 频率、磁通、转速之间又互相都有影响,所以是强耦合的多变量系统,可以先用下 图来定性地表示。 2)在异步电机中,电流乘磁通产生转矩,转速乘磁通得到感应电动势,由于 它们都是同时变化的,在数学模型中就含有两个变量的乘积项。这样一来,即使不 考虑磁饱和等因素,数学模型也是非线性的。 3)三相异步电机定子有三个绕组,转子也可等效为三个绕组,每个绕组产生 磁通时都有自己的电磁惯性[3],再算上运动系统的机电惯性,和转速与转角的积分 关系,即使不考虑变频装置[4]的滞后因素,也是一个八阶系统。 总起来说,异步电机的动态数学模型是一个高阶、非线性、强耦合的多变量系 统。 2.2 三相异步电动机的多变量非线性数学模型 在研究异步电动机的多变量非线性数学模型时,常作如下的假设: (1)忽略空间谐波,设三相绕组对称,在空间互差 120°电角度,所产生的磁 动势沿气隙周围按正弦规律分布。 2
课程设计说明书 (2)忽略磁路饱和,各绕组的自感和互感都是恒定的。 (3)忽略铁心损耗。 (4)不考虑频率变化和温度变化对绕组电阻的影响。 异步电机的数学模型由下述电压方程、磁链方程、转矩方程和运动方程组成。 2.2.1 电压方程 u A  Ri sA 三相定子绕组的电压平衡方程为 d A  d t d B  d t d C  d t Ri sC Ri sB     u u  B C 与此相应,三相转子绕组折算到定子侧后的电压方程为 u a  Ri a r  u b  Ri b r  u c  Ri c r  d a  d t d b  d t d c  d t au , bu , cu —定子和转子相电压的瞬时值; Cu , bi , ci —定子和转子相电流的瞬时值; a , Bu , 式中 Au , ai , , , , Ci Bi Ai C , A , B , Rs, Rr—定子和转子绕组电阻 上述各量都已折算到定子侧,为了简单起见,表示折算的上角标“ ’”均省略, c —各相绕组的全磁链; b , 以下同此。 电压方程的矩阵形式 将电压方程写成矩阵形式,并以微分算子 p 代替微分符号 d /dt A B u   u   u C  u   u  u   b a c            R s 0 0 0 0 0           0 R s 0 0 0 0 0 0 R s 0 0 0 0 0 0 R r 0 0 0 0 0 0 R r 0  p   A    B   C    a   b     c           0 0 0 0 0 R r           i  A  i  B  i C  i  a  i b  i   c           3 (2-1)
课程设计说明书 或改写成 u  Ri  p 2.2.2 磁链方程 每个绕组的磁链是它本身的自感磁链和其它绕组对它的互感磁链之和,因此, 六个绕组的磁链可表达为                      L AB L BB L CB L aB L bB L cB L AC L BC L CC L aC L bC L cC L AA L BA L CA L aA L bA L cA   A    B   C    a   b     c 或改写成 Li (2-2)式中,L 是 6×6 电感矩阵,其中对角线元素 AAL , BBL , CCL , aaL , L Ab L Bb L Cb L ab L bb L cb L Ac L Bc L Cc L ac L bc L cC L Aa L Ba L Ca L aa L ba L ca i  A  i  B  i C  i  a  i b  i   c (2-2)                     bbL , ccL 是各有关绕组的自感,其余各项则是绕组间的互感。 实际上,与电机绕组交链的磁通主要只有两类:一类是穿过气隙的相间互感磁 通,另一类是只与一相绕组交链而不穿过气隙的漏磁通,前者是主要的。 电感的种类和计算如下。 定子漏感 lsL ——定子各相漏磁通所对应的电感,由于绕组的对称性,各相漏 感值均相等; 转子漏感 lrL ——转子各相漏磁通所对应的电感; 定子互感 msL ——与定子一相绕组交链的最大互感磁通; 转子互感 mrL ——与转子一相绕组交链的最大互感磁通。 由于折算后定、转子绕组匝数相等,且各绕组间互感磁通都通过气隙,磁阻相 同,故可认为 msL = mrL 。 自感表达式对于每一相绕组来说,它所交链的磁通是互感磁通与漏感磁通之 和,因此,定子各相自感为(自感=漏感+互感) AA BB  ms   L L   lL L bb L CC L  转子各相自感为 L L cc aa 两相绕组之间只有互感。互感又分为两类: (1)定子三相彼此之间和转子三相彼此之间位置都是固定的,故互感为常值; (2)定子任一相与转子任一相之间的位置是变化的,互感是角位移θ的函数。 L ms lL   s r 4
课程设计说明书 第一类固定位置绕组的互感,三相绕组轴线彼此在空间的相位差是±120°, 在假定气隙磁通为正弦分布的条件下,互感值应为 L ms cos 120  L ms cos( 120  )  1 2 L ms 于是 L AB  L BC  L CA  L BA  L CB  L AC L ab  L bc  L ca  L ba  L cb  L ac   1 L 2 ms 1 L 2 ms (2-3) 第二类变化位置绕组的互感,定、转子绕组间的互感,由于相互间位置的变化, 可分别表示为 L aA L Aa   L Bb  L bB  L Cc  L cC  L ms cos L Ab  L bA  L Bc  L cB  L Ca  L aC  L ms cos(   120 )  (2-4) L Ac  L cA  L Ba  L aB  L Cb  L bC  L ms cos(   120 )  当定、转子两相绕组轴线一致时,两者之间的互感值最大,就是每相最大互感 msL 。 整理以上各式,即得完整的磁链方程,显然这个矩阵方程是比较复杂的,为了 方便起见,可以将它写成分块矩阵的形式 AA BA L L L CA L L bA L aA cA BB L AB L L CB L L bB L aB cB BC L AC L L CC L L bC L aC cC Aa Ba L L L Ca L L ba L aa ca Ab Bb L L L Cb L L bb L ab cb Bc L Ac L L Cc L L L bc ac cC           A B i   i   i C  i   i  i   b c a              A    B   C    a   b     c           式中 T T sΨ C A B  a b  Ψ r c  T Asi i  i a r i Ti i B i C c i b 1 2 rrL           L ms 1 2 1 2   L l r  L ms  L l r L ms L ms L ms 1 2    1 2 1 2 L ms L ms         5 L ms L ms  L l r           Ψ Ψ    s r    L L    ss rs L L sr rr    i i    s r    (2-5) (2-6)
课程设计说明书 rrL           L l r L ms L ms  L ms 1 2 1 2  L rs L  sr T  msL  1 2 L ms r   L l L ms L ms 1 2 cos  120   120   cos( cos( )  )    1 2 1 2 L ms  L ms L ms         cos( L l r cos(      srL (2-7) )  )  cos( cos( 120   120   cos  )  )       120   cos  120   值得注意的是, 和 两个分块矩阵互为转置,且均与转子位置θ有关,它 们的元素都是变参数,这是系统非线性的一个根源。为了把变参数转换成常参数须 rsL 利用坐标变换,后面将详细讨论这个问题。 如果把磁链方程代入电压方程中,即得展开后的电压方程 u  Ri   Ri  Li ( p iL d d t )  Ri   d d L   i  iL d d t  iL d d t 式中,Ldi /dt 项属于电磁感应电动势中的脉变电动势(或称变压器电动势), (dL / d)i 项属于电磁感应电动势中与转速成正比的旋转电动势。 2.2.3 转矩方程 根据机电能量转换原理,在多绕组电机中,在线性电感的条件下,磁场的储能 和磁共能为 ' mW   m (电流约束为常值), ' WW m  m  1 2 ψi T  1 2 T i Li 而电磁转矩等于机械角位移变化时磁共能的变化率 且机械角位移 m =  / np ,于是 T e  ' W  m   m i  const .  ' Wn  m   p i  const . 整理上式可得 T e  1 2 n p T iLi     1 2 n p T i      0 L  rs   L  sr   0 i      又由于 6
课程设计说明书 T i  [ i T s i T r ]  [ i A i B i C i a i b i c ] 则 T e  1 2 n p   i T r  iL  rs    i T s  s iL  sr   r   转矩方程的三相坐标系形式 T sr    msL rs L      iL 1   rs  2   [( Ln ii  ms aA p ( ii ii   bA cB ( ii ii   cA aB  ii aC ii bC n T r i  p L T e T e   cos  120   120   )  )  cos( cos( cos( cos( 120   cos  120   s r  i T s  iL   sr    sin) ii ii   bB cC ) ) sin( 120    120 )] sin( )    )  )  cos( cos( 120   120   cos  )  )       (2-8) 应该指出,上述公式是在线性磁路、磁动势在空间按正弦分布的假定条件下得 出来的,但对定、转子电流对时间的波形未作任何假定,式中的 i 都是瞬时值。 因此,上述电磁转矩公式完全适用于变压变频器供电的含有电流谐波的三相异 步电机调速系统。 2.2.4 电力拖动系统运动方程 在一般情况下,电力拖动系统的运动方程式是 (2-9) T e  T L  J n p d   dt D n p   K n p TL —— 负载阻转矩; J—— 机组的转动惯量; D —— 与转速成正比的阻转矩阻尼系数; K —— 扭转弹性转矩系数。 对于恒转矩负载,D = 0 ,K = 0 ,则 T e  T L  J n p d  d t 7
分享到:
收藏