同步电机三相短路电流和转矩计算Word下载.docx
《同步电机三相短路电流和转矩计算Word下载.docx》由会员分享,可在线阅读,更多相关《同步电机三相短路电流和转矩计算Word下载.docx(12页珍藏版)》请在冰豆网上搜索。
xaf定子绕组与劢磁绕组间得互感电抗
r定子绕组电阻
1、2电流变化量得状态空间方程式
同步电机突然短路时各绕组电流得变化量得计算可运用以派克分量表示得状态空间方程式
方程中各下标变量得含义为
d---纵轴,q---横轴,f----励磁绕组,D---纵轴阻尼绕组,Q---横轴阻尼绕组,a---定子绕组
上式可简化为
化作电流变化量得常系数一阶微分方程组形式
在三相短路时若励磁电压不可调,则
由于电流不能突变,t=0瞬间电流变化量得初值
将电压变化量与电流变化量得初值代入微分方程,用数值计算得龙格---库塔法即可求出
t=0+Δh时刻得各电流变化量,反复计算则可求得各个时刻得,叠加短路前绕组电流
可得短路时电流全量
用派克逆变换可得定子三相电流,以a相为例
ia=
短路t=0时转子位置角
2.、同步电机三相短路电流计算例题与程序
电机参数
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
xaf=0、85,xaD=0、85,xfD=0、85,xaQ=0、45
原始运行条件为额定负载
U[0]=1,i[0]=1,φ[0]=0、5548(单位为弧度,相对于cosφ=0、8)
短路时得转子位置角
θ0=3、1416
三相短路计算程序:
CMSHORT3、M
%part1
ra=0、005;
rf=、000656;
rzd=、00151;
rzq=、00159;
xd=1、0;
xq=、60;
xf=1、03;
xzd=、95;
xzq=、70;
xaf=、85;
xazd=、85;
xfzd=、85;
xazq=、450;
u0=1、0;
i0=1、0;
phas=、5548;
cita0=3、1416;
p=31、4160;
h=、5236;
x1=[-1、0,0、0,0、85,0、85,0、0;
0、0,-0、60,0、0,0、0,0、45;
-0、85,0、0,1、03,0、85,0、0;
-0、85,0、0,0、85,0、95,0、0;
0、0,-0、45,0、0,0、0,0、70];
z1=[-0、005,0、6,0、0,0、0,-0、45;
-1、0,-0、005,0、85,0、85,0、0;
0、0,0、0,0、000656,0、0,0、0;
0、0,0、0,0、0,0、00151,0、0;
0、0,0、0,0、0,0、0,0、00159];
g0=(u0*sin(phas)+xq*i0)/(u0*cos(phas)+ra*i0);
g0=atan(g0)-phas;
ud0=u0*sin(g0);
uq0=u0*cos(g0);
di0=i0*sin(g0+phas);
qi0=i0*cos(g0+phas);
eq0=uq0+ra*qi0+xd*di0;
fi0=eq0/xaf;
du=[-ud0,-uq0,0、0,0、0,0、0];
x2=inv(x1);
z2=-x2*z1;
i1=x2*du'
;
y=[0、0,0、0,0、0,0、0,0、0];
%part2
t=0、0;
fori=1:
5
b(i)=y(i);
end
dy=z2*(y)、'
+i1;
tt
(1)=t;
di
(1)=y
(1)+di0;
qi
(1)=y
(2)+qi0;
fi
(1)=y(3)+fi0;
zdi
(1)=y(4);
zqi
(1)=y(5);
ai
(1)=di
(1)*cos(t+cita0)-qi
(1)*sin(t+cita0);
tm
(1)=qi
(1)*(-xd*di
(1)+xaf*fi
(1)+xazd*zdi
(1))-di
(1)*(-xq*qi
(1)+xazq*zqi
(1));
%while(tt(n)<
p)
%k=0;
forn=2:
10
a
(1)=h/2;
a
(2)=a
(1);
a(3)=h;
a(4)=h;
fork=1:
3
c(i)=b(i)+a(k)*dy(i);
y(i)=y(i)+a(k+1)*dy(i)/3、0;
end
dy=z2*(c)、'
y(i)=y(i)+h*dy(i)/6、0;
t=t+h;
%k=k+1;
%part3
tt(n)=t;
di(n)=y
(1)+di0;
qi(n)=y
(2)+qi0;
fi(n)=y(3)+fi0;
zdi(n)=y(4);
zqi(n)=y(5);
ai(n)=di(n)*cos(t+cita0)-qi(n)*sin(t+cita0);
tm(n)=qi(n)*(-xd*di(n)+xaf*fi(n)+xazd*zdi(n))-di(n)*(-xq*qi(n)+xazq*zqi(n));
plot(tt,ai);
xlabel('
t(rad)'
);
ylabel('
ia'
title('
3-phasefaultofsynchronousmachine'
grid
3.程序说明:
X1为状态空间方程得Xdq0矩阵
Z1为状态空间方程得Zdq0矩阵
X2=
I1=
程序%part1第一部分输入初始数据,建立状态空间方程;
程序%part2第二部分用定步长四阶龙格---库塔法解一阶微分方程组,计算tn+1=tn+h时刻各绕组电流得变化量。
时间步长h用标么值,单位为弧度,h取1/6π
定步长四阶龙格---库塔法解一阶微分方程组得基本公式:
y1’=f1(t,y1,y2,…,ym),y1(t0)=y10
y2’=f2(t,y1,y2,…,ym),y2(t0)=y20
……、
Ym’=fm(t,y1,y2,…,ym),ym(t0)=ym0
微分方程组由tj积分一步到tj+1=tj+h
Yi,j+1=yij+(k1i+2k2i+2k3i+k4i)h/6,i=1,2,…,m
k1i=fi(tj,y1j,y2j,…,ymj),i=1,2,…,m
k2i=fi(tj+h/2,y1j+k11h/2,…,ymj+k1mh/2),i=1,2,…,m
k3i=fi(tj+h/2,y1j+k21h/2,…,ymj+k2mh/2),i=1,2,…,m
k4i=fi(tj+h,y1j+k31h,…,ymj+k3mh),i=1,2,…,m
程序%part3第三部分输出计算结果,tm为转矩
Tm=iqψd+idψq
ψd=-xdid+xafif+xaDiD
ψq=-xqiq+xaQiQ
本程序计算时间周期为t=0—10X1/6π,绘画取20π
4、程序计算结果:
T=0--10X1/6π,取10个点
纵轴短路电流
di
di=
Columns1through7
0、7966190047267142、23674、4478626375923196、8263978366408118、7369322382886849、68989、427860674567469
Columns8through10
8、05455、98053、734356462313833
A相短路电流
ai
ai=
-0、7966-0、9644241603207350、97802、4496686213438025、6777204819841818、4172515171420659、427882627330178
8、3622653501447395、9184325579004333、227399551484802
横轴短路电流
qi
qi=
0、63681、94522、66602、4495933972144211、51170、3252-1、494184215750709
-2、7523-3、3834-3、227468137315201
励磁绕组电流
fi
fi=
2、28842、5178174245469513、2548509777309844、0544、69255、7904、941547215637447
4、4986640266012453、81283、3557
纵轴阻尼绕组电流
zdi
zdi=
00、8597117550616702、1774129105476173、59164、7227836673887345、2743848475096735、1786
4、2841844368405483、89831、679680254626384
横轴阻力尼绕组电流
zqi
zqi=
00、86861、29281、30910、578725790607257-0、3501-1、352763638980671
-2、2395-2、5611-2、4590
电磁转矩
tm
tm=
0、85512、9771766601552444、8382496751246415、5892558297652774、5397477838432851、849594513339010-1、3946
-3、827143999454925-4、6864-4、4002
三相短路电流与转矩波形图
A相短路电流图(故障前额定负载,θ0=3、1416)
转矩
A相短路电流(短路前空载,转子角θ0=π)
A相短路电流(短路前空载,转子角θ0=π/2)
参考文献:
1.同步电机运行