同步电机的三相突然短路的MATLAB计算.doc

(完整版)同步电机的三相突然短路的MATLAB计算
同步电机的三相突然短路的MATLAB计算
姓名：林启星 专业:电力电子与电力传动 学号：100120039
1、 新建一个M文件存放主程序如下：
％＊＊*******＊*＊****＊＊＊**＊**＊*＊*＊＊*＊＊＊＊
%定义变量及变量赋值
%＊*＊**＊***＊***＊*＊*＊＊*＊***＊＊＊**＊*＊***
global A B U0
R=0.005；Rf=0。000656；RD=0。00151;RQ=0。00159；
xd=1.0;xq=0。60；xf=1。03;
xD=0。95;xQ=0。70;xad=0.85;xaq=0。45；
u0=1;i0=1；w0=0。5548；
A=［—xd，0,xad,xad，0;0，-xq,0,0，xaq；-xad，0，xf，xad,0;—xad,0,xad，xD,0;0,—xaq，0,0,xQ］；
B=[-R，xq,0,0，—xaq;—xd，—R,xad，xad，0;0，0,Rf，0，0；0，0,0,RD,0；0,0，0，0,RQ］；
％＊**＊*＊*****＊＊＊*＊**＊*****＊＊*＊＊＊*＊**＊
％求出功率角a0及ud0,uq0，id0,iq0,if0
%＊＊*＊＊*＊*＊＊＊＊****＊＊*＊＊＊**＊**＊*＊＊＊**＊
a0=atan(（u0*sin(w0）+xq*i0)/（u0＊cos(w0)+R＊i0））-w0;
ud0=u0*sin(a0）；
uq0=u0*cos（a0）;
id0=i0＊sin(a0+w0);
iq0=cos(a0+w0);
Eom=uq0+R＊iq0+xd*id0;
If0=Eom/xad;
U0=［-ud0;-uq0;0；0;0］；
%***＊**＊**＊*＊*＊＊＊＊**＊**＊***＊*＊＊＊**＊＊
％通过使用ode45来计算各绕组电流
％*＊*＊**＊＊＊***＊****＊**＊**＊＊**＊***＊***
I0=[0;0；0;0;0］;
[t，i］=ode45（'IdqfDQ'，[0，2000]，I0）；
id=i(:，1)+id0;
iq=i（:,2)+iq0;
If=i(:,3)+If0；
iD=i(：，4)；
iQ=i（：，5);
ia=id.＊cos（t+pi）-iq。＊sin(t+pi）;
%＊＊＊*＊＊＊*****＊＊***＊＊＊**＊＊＊**＊＊＊**＊**
%三相短路时的电磁转矩Te
%＊*****＊**＊**＊＊*＊*＊*＊＊*＊****＊＊**＊＊**
Wd=—xd*id+xad＊If+xad*iD;
Wq=—xq*iq+xaq＊iQ；
Te=Wd.*iq-Wq.*id；

％＊＊＊*****＊***＊＊*＊**＊**＊*＊***＊＊*＊*＊**
%计算ia及Te
％＊**＊***＊＊*****＊＊＊＊＊*＊＊********＊*＊*＊
Ten=Te(end，1）
ian=ia(end,1)
％＊*＊＊*＊＊＊**＊＊*＊＊＊**＊＊***＊＊＊**＊******
％各变量图形输出
%**＊＊＊＊＊＊＊*＊*＊*＊*＊**＊＊＊*＊**＊＊*＊＊＊＊*＊
figure(1）;
plot(t，ia）
title（'ia=f（t）'）；
xlabel(’t’);
ylabel(’ia');
figure（2）；
plot(t,If)
title（’If=f(t)’);
xlabel('t')；
ylabel('If'）;
figure（3）;
plot(t，iD)
title（'iD=f（t）’)；
xlabel（’t’）；
ylabel（’iD’）；
figure(4）；
plot（t，iQ）
title（'iQ=f（t）’）;
xlabel（’t’）；
ylabel('iQ')；
figure(5）；
plot（t,Te)
title（'Te=f（t）'）;
xlabel（’t'）;
ylabel('Te’)；
2、新建一个M文件（IdqfDQ。m）存放子程序程序如下：
function dI=IdqfDQ(t，i)
global A B U0
dI=A\U0-A\B＊i;
3、 运行结果如下：
Ten =:
0。0160
ian =
0.6706