完整word版直线电机本体建模.docx

1、完整word版直线电机本体建模基于 Simulink 的直线电机本体建模电磁发射课题组2015年 10月 29日1直线感应电动机的等效电路直线电机在结构上可看作是沿径向剖开并将圆周展为直线的旋转电机，如图 1 所示。直线感应电动机的稳态特性近似计算方法基本可以沿用旋转感应电动机的等效电路 1 。图 1 旋转电机演变为直线电机示意图对于旋转异步电机而言， 与电机绕组交链的磁通主要有两类： 一类是穿过气隙的相间互感磁通； 另一类是只与一相绕组交链而不穿过气隙的漏磁通， 前者是主要的。 定子各相漏磁通所对应的电感称作定子漏感 Lls ，由于绕组的对称性，各相漏感值均相等。同样，转子各相漏磁通则对应于

2、转子漏感 Llr 。对于每一相绕组来说， 它所交链的磁通是互感磁通和漏感磁通之和，因此，定子各相自感为：LAA LBB LCC Lms Lls（1）转子各相自感为：LaaLbb Lcc Lmr Llr Lms Llr（2）两相绕组之间只有互感，互感又分为两类：1） 定子三相彼此之间和转子三相彼此之间位置都是固定的，故互感为常值；2） 定子任一相与转子任一相之间的位置是变化的，互感是角位移 的函数。由于三相绕组轴线彼此在空间的相位差为120 ，因此互感为：Lms cos 120Lms cos1201 Lms（3）2于是：LABLBCLCALBALCBLAC1Lms（4）2LabLbcLcaLba

3、LcbLac1Lmr1Lms（5）22定转子绕组间的互感由于相互间的位置的变化，为：LAaLaALBbLbBLcCLCcLms cosLAbLbALBcLcBLaCLCaLms cos120LAcLcALBaLaBLbCLCbLms cos120（6）（7）（8）以上是针对旋转异步电机的参数的推到过程，而对于直线电机，文献 3 中作者给出了圆筒形直线感应电动机的等效电路，如所示：图 2 圆筒形直线电机的等效电路图 2 中， Rs 和 X s 分别代表初级绕组的电阻和漏抗； Rm 代表励磁电阻； X m 代表励磁电抗； r20 代表次级表面电阻； x20 代表次级表面电抗； Red 代表边端效应

4、影响纵向边电功率产生的损耗折算成的等效电阻； Rr 代表在次级铜层中的折算的电阻值。在该文献 3 中次级使用的是导电层和导磁层所构成的复合材料。至于图 2 中的相关参数的计算过程， 在该文献中都有详细的说明， 不再赘述。文献 1 中给出了计及边端效应的等效电路，如所示：图 3 计及边端效应的等效电路图 3 中， b0 为励磁电纳（ ）； r1 为初级绕组电阻； x1为初级绕组漏电抗； r2 为次级导体电阻折算到初级的换算值， Re 为边端效应消耗功率的等效电阻折算到初级的换算值。2直线感应电机的数学模型（ 1）电压方程参看海军工程大学鲁军勇在文献 4 中给出的电压方程，即：vds Rsids

5、Dvds Rsids D1Rr idr D1Rr iqr DdsVeqsdsVeqs（9）VeVdrqrqrVeVdr式中： Rs 为通电段定子绕组电阻； Rr 为通电段定子绕组电阻； Ve为同步速度， V 为动子实际速度； D 为微分算子； 。注释：对上式进行简要的推导：利用三相静止坐标系到两相任意旋转坐标系间的转换矩阵 C3s 2r可将三相静止坐标系下的定子电压方程转换到任意旋转坐标系 dq0坐标系下，即：uARs00i AdAuB0 Rs0iBdtBuC00Rs iCCuARs00i AdC3 s/2 ruBC3 s/2 r0Rs0dds1dsiBC3 s/2 rC3 s/2 ruC00

6、Rsdtqsdtqsi CudsRsidsdds01dsuqsiqsdtqs10qsudsRs idsDds1qsuqsRs iqsDqs1ds（ 10）对于转子电压方程的推导过程类似， 只是转子坐标系转换矩阵与定子坐标系的转换矩阵不一样，即：cosC3 s 2r2sin312rcos( r2)cos( r2 )33sin( r2)sin( r2（11）r3)31122利用该转换矩阵将转子电压方程由三相静止坐标系转换到两相任意旋转坐标系下，即：uaRr00iadaub0Rr0ibdtbuc00RriccuaRr00iadd C3 s/2 r C3 s/2 rC3 s/2 rubC3s/2 r0

7、Rr0ibdr1druc00RricdtqrdtqrudrRridrddr0sdruqriqrdtqrs0qrudrRr idrDdrsqruqrRr iqrDqrsdr（12）综上，将电压方程归结为：- （ 13）考虑角速度与速度间的关系，即：R vvv（ 14）R将式（ 14）带入到式（ 13）中可得：uds Rsi ds Duqs Rsi qs D0 Rr idr D0 Rr iqr DdsVeqsqsVeds（ 15）VeVdrqrqrVeVdr注意：式（13）中的 1 是电角速度， s 1 中的 为次级折算的旋v转角速度（机械量），折算关系是：，而式（ 15）中的速度 Ve是同步速度

8、（定子磁场的速度），V 是动子实际的运动速度 （机械运动速度）。（ 2）磁链方程鲁军勇在文献 4中给出的直线电机的磁链方程为：dsLlsLu1LmidsLm idrqsLlsLu1LmiqsLm iqrdrLmidsLlrLmidrqrLmiqsLlrLmiqr（ 3）电磁推力方程文献 4 中给出的直线电机电磁推力方程为：Fe3Lmqr ids2dr iqsLr注释：对上式（ 17）进行简要的推导：从电磁功率的角度入手，则：（16）（17）FevrPeTeTe mTe m（ 18）PeFenpnp mvrnp m因此，电磁推力与电磁转矩的关系为：Te（ 19）Fenp而我们知道 对于旋转异步电

9、机 而言，其电磁转矩的表达式为：Tenp Lm (is ir is ir )（ 20）结合磁链方程将式（ 20）中的转子电流分量消掉，则：irirLmisLm isLr irLrr（ 21）LmisLr irLm isrirrLr将式（ 21）带入到式（ 20）中可得：Tenp Lm (isq irdisdi rq )np LmisqrdLmisdisdrqLm isqLrLrnpLmisqi sd（ 22）rdrqLrFeTeLmisq rdisdrqnpLr疑问：式（ 17）中的系数如何理解？个人认为应该是转换矩阵的不同带来的这个系数， 因为在上面的分析中采用的都是恒功率转换矩阵， 而在鲁

10、军勇的文献中所使用的转换矩阵是恒幅值转换矩阵，下面我们验证这种猜测：由文献 7 可知恒功率转换矩阵 C3s 2r 和 C3 s1 2 r 分别为：2coscos120cos120C3 s 2 rsinsin120sin120 （23）3111222cossin12C3 s12 r2cos120sin1201（ 24）32cos120sin12012恒幅值转换矩阵为：coscos120cos120C3 s 2 r2sin120sin120（25）sin3111222cossin1C3 s12 rcos120sin1201（ 26）cos120sin1201仍然借助旋转异步电机的电磁转矩来推导电磁

11、推力，将式（ 25）和（ 26）带入到文献 7 中给出的电磁转矩表达式中，即：Te np Lms i A iai BibiC i c siniA ib iB iciC ia sin( 120 )iA iciB iaiC ibsin(120 )（27）利用恒幅值转换矩阵将ABC 坐标系上的定、转子电流转换到 dq0坐标系，即：i Acos 1sin11isdi Bcos1120sin11201isqiCcos1120sin11201is 0i acos rsinr1irdi bcosr120sinr1201irqi ccosr120sinr1201ir 0由于推导过程相当复杂， 但是我们发现在文

12、献 7 中作者指出：在化简过程中的零轴分量完全抵消了，所以对比两种情况的转换矩阵，可做如下的推导：当使用恒功率转换矩阵时：Te np Lms i Ai aiB ibiC ic sini Aib i B ic iC i a sin(120 )i Aici B iai Ci bsin(120 )（ 28）2np Lms K当使用恒幅值转换矩阵时：Te np Lms i Ai aiB ibiC ic sini Aib i B i c iC i a sin(120 )i Aici B iai Ci bsin(120 )（ 29）np Lms K因此，采用恒幅值转换矩阵运算时的电磁转矩为采用恒功率转换矩

13、阵运算时电磁转矩的Te3 np Lm (isqirdisd irq )1.5 倍，即2。将其带入到式（ 19）中可得电磁推力为：FeTe3Lm ( isqirdisd irq )3Lm ( isq rd i sd rq ) （30）np22Lr（ 4）运动方程文献 4 中给出的直线电机在发射阶段的运动方程为：M mdvdtFeBv2M m g（ 31）式中： M 风摩系数；为负载质量； m 为动子本体质量；为滑动摩擦系数。Fe 为电磁推力；B 为注释：个人认为如果按照式（ 23）来编写状态方程时，较难列写出速度的状态方程，因为我们知道状态方程的形式为： x& Ax Bu ，考虑是否能将运动方程

14、简化为这种容易列写状态方程的形式， 为此，参看文献56 中给出的运动方程的形式，即：dvM m Fe FL Bv v （ 32）式中： FL 负载阻力；v 机械运动速度；Bv 与速度有关的阻尼系数；将电磁推力的表达式带入到式（ 32）中，得：dvBv vLmdr iqsqr idsFL（ 33）dtM 总Lr M 总M 总文献 6 中作者将粘滞阻尼系数取：Bv0.2 N s m 。3状态方程推导状态方程是指刻画系统输入和状态关系的表达式。 状态向量所满足的向量常微分方程称为控制系统的状态方程， 状态方程是控制系统数学模型的重要组成部分。对于线性系统而言，我们知道其状态方程的形式为：x& A t

15、 x B t u（ 33）y C t x D t u状态变量的选取：直线电机作为异步电机的一种，同样具有 4 阶电压方程和 1 阶运动方程，因此其状态方程也应该是 5 阶的，因此必须选取 5 个状态变量 7 。在旋转异步电机中可选的变量共有 9 个，即转速 、4个电流变量 i sd 、 i sq 、 i rd 、 i rq 和 4 个磁链变量 sd 、 sq 、 rd 、 rq 。个人认为针对直线电机而言，将其中的转速换成速度 v ，另外，转子电流 i rd 、 irq 是不可测的，因此不宜作为状态变量，故只能选择定子电流 i sd 和 isq ，另两个状态变量必须是转子磁链 rd 、 rq

16、，或定子磁链 sd 、 sq 。为了推导出状态方程，需要结合电压方程（ 15）和磁链方程，现将两个方程重新列出：电压方程：uds Rsids Duqs Rsiqs D0 Rr idr D0 Rr iqr Dds1qsqs1ds（34）drsqrqrsdr磁链方程：sdLsi sdLmirdsqLsi sqLmirqLmisd（35）rdLr irdrqLmi sqLri rq式（ 34）中： s1；注意到如果采用转子磁链定向，则有rqLmisqLrirq0，即irqLmi sqLmisdLrd。r。由式（ 35）的第 3 式可得： irdLr将磁链方程代入到电压方程中，消去其中的sd 、 sq

17、 、 i rd 、 irq ：uR isdp L isdL ird L isqL irqsdssmdqssmusqRsisqpLsi sqLmirqdqsLsi sdLmird0 R 1 L ip rLrrdmsdrddqsrq0 R 1 L ip rLrrqmsqrqdqsrd由上式的第 3、4 式可得：d rdRrLm isdRr1 rddqs rqdtLrLrd rqRrLm isqRr1 rqdqs rddtLrLrd iRsL2rRr L2m i iLmLm dtsdLsL2rsddqs sqLsLrTrrdLsLrrqd i iRsL2rRr L2m iLm Lmdtsqdqs s

18、dLsL2rsqLsLrrdLs Lr TrrqL2m 1其中： 电机漏磁系数，LsLrTrLrTrRr转子电磁时间常数，dqs1usdLsusqLs将上述推导出的状态方程写成矩阵的形式，则：Rs L2r Rr L2mLmLm 2Ls LrdqsLs Lr TrLs Lr10i sdRs L2r Rr Lm2Lm Lmi sdLsdqsdi sqLs L2rLs LrLs Lr Trisq1usd0drdLm1rdL susqt0dqs TrTr00rqrq0Lmdqs100TrTr上式与式（33）即 dvBv vLmdr iqsqri dsFL组成dtM 总Lr M总M 总直线电机的状态方程。从上述的状态方程中可知，状态变量为：Xv, sd , sq,i sd ,i sqT（ 36）输入变量为：TUusd , usq , 1, FL（37）如果在推导状态方程时使用的是鲁军勇文献 4 中给出的电压方程和磁链方程，则只需对上述的状态方程做如下的修改：4 S-Function 的编写4.1 S 函数的原理Simulink 中的大部分模块都具有一个输入向量 u 、一个输出向量 y 和一个状态向量 x ，如所示：引入 S 函数的引入 S 函数的目的是为了使 Simulink 有能力构作一般的仿真框图，去处理如下各种系统的仿真： 连续系统、离散系