仿真电机变频转速转矩

1、异步电动机全速度范围内的矩速特性仿真The simulation of Asynchronous motor full speed range of torque speed characteristic姓名：胡珍，班级：09电牵（3）班，学号：20090210470402摘要：此论文主要是运用MATLAB编写m函数程序，绘制出牵引电动机低频到高频时，恒磁通,恒转矩起动，恒功率运行的矩速曲线。根据以往所学电机学的知识可以知道电磁转矩正比于气隙磁通的平方。若在调节时维持Es/fs不变，亦即Im为常数，则电磁转矩完全由转子的转差频率fsl所决定，而与定子频率fs无关。关键字：异步电动机 MATLA

2、B仿真 电动机变频 恒磁通 恒转矩 恒功率KEY WORDS: asynchronous motor Matlab simulation Motor frequency conversion Constant flux Constant torque constant powerABSTRACE: This thesis mainly is to use MATLAB m function program，Draw the traction motor low frequency to high frequency, constant flux, constant torque startin

3、g, constant power operation of the moment speed curve. According to the previous learning motor learning knowledge can know electromagnetic torque is proportional to the square of the air gap flux. If in the regulation when maintaining Es/fs is changeless, i.e. Im for constant, the electromagnetic t

4、orque completely by the rotor of the slip frequency FSL decided, and stator frequency fs have nothing to do.正文：由变频调速异步电动机的稳态等值电路可知通过空气隙传入转子的电磁功率和电磁转矩的一般表达式，为获得不变的气隙磁通，需要对电动机按照比值Eg/fs不变进行调节。可得转矩的表达式如下示：转子的临界转差频率：fm=Rr2Llr电动机的最大转矩为：Tm=mnp2(Egfs)214Llr运用matlab编写m函数文件，绘制出电动机低频到高频的矩速曲线，其m函数如下：function f=

5、fi(f)a=2;m=3;np=2;Rs=2.497;Rr=5.503;Ls=0.0217;Lr=0.0217;A=m*np/(2*pi);s=0.02:0.001:1;fs=40;n0=60*fs/np;n=n0-s.*n0;B=2*pi*Lr*fs;T=A*a2*(Rr*fs*s./(Rr2+B2*s.2);T1=fs*s./(0.012+2*s.2);plot(T,n) hold onfs=60;A=m*np/(2*pi);s=0.02:0.001:1;n0=60*fs/np;n=n0-s.*n0;B=2*pi*Lr*fs;T=A*a2*(Rr*fs*s./(Rr2+B2*s.2);T1

6、=fs*s./(0.012+2*s.2);plot(T,n)hold onfs=80;A=m*np/(2*pi);s=0.02:0.001:1;n0=60*fs/np;n=n0-s.*n0;B=2*pi*Lr*fs;T=A*a2*(Rr*fs*s./(Rr2+B2*s.2);T1=fs*s./(0.012+2*s.2);plot(T,n)hold onfs=100;A=m*np/(2*pi);s=0.02:0.001:1;n0=60*fs/np;n=n0-s.*n0;B=2*pi*Lr*fs;T=A*a2*(Rr*fs*s./(Rr2+B2*s.2);T1=fs*s./(0.012+2*s.2

7、);plot(T,n)hold onfs=120;A=m*np/(2*pi);s=0.02:0.001:1;n0=60*fs/np;n=n0-s.*n0;B=2*pi*Lr*fs;T=A*a2*(Rr*fs*s./(Rr2+B2*s.2);T1=fs*s./(0.012+2*s.2);plot(T,n)hold onfs=160;A=m*np/(2*pi);s=0.02:0.001:1;n0=60*fs/np;n=n0-s.*n0;B=2*pi*Lr*fs;T=A*a2*(Rr*fs*s./(Rr2+B2*s.2);T1=fs*s./(0.012+2*s.2);plot(T,n)hold onfs=180;A=m*np/(2*pi);s=0.02:0.001:1;n0=60*fs/np;n=n0-s.*