同步电机三相短路电流和转矩计算.docx
《同步电机三相短路电流和转矩计算.docx》由会员分享,可在线阅读,更多相关《同步电机三相短路电流和转矩计算.docx(12页珍藏版)》请在冰豆网上搜索。
同步电机三相短路电流和转矩计算
同步电机三相短路电流与电磁转矩计算
编写佘名寰
本文就是按照陈珩教授所著得‘同步电机运行基本理论与计算机算法’一书介绍得算法与例题计算同步电机得三相短路电流。
计算程序用MATLAB语言编写,计算结果与书中结果基本一致。
本文可供电力系统电气技术人员与大专院校电力专业学生参考。
1.计算方法
1、1初始数据计算
由短路前得机端电压u[0],定子绕组电流i[o],与功率因数角φ[0]求得短路前得功率角
δ0=
从而得u[0],i[0]得正、交轴分量
ud[0]=u[0]sinδ0
uq[0]=u[0]cosδ0
id[0]=i[0]sin(δ0+φ[0])
iq[0]=i[0]cos(δ0+φ[0])
短路前得空载电势就是
Eq[0]=uq[0]+riq[0]+xdid[0]
励磁电流为
if[0]=Eq[0]/xaf
式中xd为同步电机正轴同步电抗
xq同步电机交轴同步电抗
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)
%k=0;
forn=2:
10
a
(1)=h/2;
a
(2)=a
(1);
a(3)=h;
a(4)=h;
fork=1:
3
fori=1:
5
c(i)=b(i)+a(k)*dy(i);
y(i)=y(i)+a(k+1)*dy(i)/3、0;
end
dy=z2*(c)、'+i1;
end
fori=1:
5
y(i)=y(i)+h*dy(i)/6、0;
end
t=t+h;
fori=1:
5
b(i)=y(i);
end
dy=z2*(y)、'+i1;
%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));
end
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=
Columns1through7
-0、7966-0、9644241603207350、97802、4496686213438025、6777204819841818、4172515171420659、427882627330178
Columns8through10
8、3622653501447395、9184325579004333、227399551484802
横轴短路电流
qi
qi=
Columns1through7
0、63681、94522、66602、4495933972144211、51170、3252-1、494184215750709
Columns8through10
-2、7523-3、3834-3、227468137315201
励磁绕组电流
fi
fi=
Columns1through7
2、28842、5178174245469513、2548509777309844、0544、69255、7904、941547215637447
Columns8through10
4、4986640266012453、81283、3557
纵轴阻尼绕组电流
zdi
zdi=
Columns1through7
00、8597117550616702、1774129105476173、59164、7227836673887345、2743848475096735、1786
Columns8through10
4、2841844368405483、89831、679680254626384
横轴阻力尼绕组电流
zqi
zqi=
Columns1through7
00、86861、29281、30910、578725790607257-0、3501-1、352763638980671
Columns8through10
-2、2395-2、5611-2、4590
电磁转矩
tm
tm=
Columns1through7
0、85512、9771766601552444、8382496751246415、5892558297652774、5397477838432851、849594513339010-1、3946
Columns8through10
-3、827143999454925-4、6864-4、4002
三相短路电流与转矩波形图
A相短路电流图(故障前额定负载,θ0=3、1416)
转矩
A相短路电流(短路前空载,转子角θ0=π)
转矩
A相短路电流(短路前空载,转子角θ0=π/2)
转矩
参考文献:
1.同步电机运行