电力系统潮流及短路电流计算程序汇总.docx
- 文档编号:11689982
- 上传时间:2023-03-30
- 格式:DOCX
- 页数:13
- 大小:103.18KB
电力系统潮流及短路电流计算程序汇总.docx
《电力系统潮流及短路电流计算程序汇总.docx》由会员分享,可在线阅读,更多相关《电力系统潮流及短路电流计算程序汇总.docx(13页珍藏版)》请在冰豆网上搜索。
电力系统潮流及短路电流计算程序汇总
班级:
姓名:
学号:
一、作业要求
编写程序计算图1所示算例系统的潮流及三相短路电流。
潮流计算:
方法不限,计算系统的节点电压和相角。
短路电流:
4号母线发生金属性三相短路时(zf=0),分别按照精确算法和近似算法计算短路电流、系统中各节点电压以及网络中各支路的电流分布,并对两种情况下的计算结果进行比较。
二、电路图及参数
图13机9节点系统表19节点系统支路参数
表29节点系统发电机参数
表39节点系统负荷参数
三、计算步骤
(1)进行系统正常运行状态的潮流计算,求得Ui
(2)形成不含发电机和负荷的节点导纳矩阵YN;
(0)
)的并联组合;)和导纳yi(1/jxdi(3)将发电机表示为电流源Ii(Ei/jxdi节点负荷用恒阻
抗的接地支路表示;形成包括所有发电机支路和负荷支路的节点导纳矩阵Y,即在YN中的发电机节点和负荷节点的自导纳上分别增加发电机导纳yi和负荷导纳yLD,i
(yLD,i
SLD,iPLDijQLDi
);Vi2Vi2
1
*
(4)利用ZY
,计算节点阻抗矩阵,从而得到阻抗矩阵中的第f列;
(5)利用公式(6-7)或(6-10)计算短路电流;
(6)利用公式(6-8)或(6-11)计算系统中各节点电压;
(7)利用公式(6-9)计算变压器支路的电流;对输电线路利用П型等值电路计算支路电流。
四、计算结果
节点导纳矩阵Yn:
Columns1through5
0-17.3611i000+17.3611i000-16.0000i000000-17.0648i000+17.3611i003.3074-39.3089i-1.3652+11.6041i000-1.3652+11.6041i2.5528-17.3382i000-1.9422+10.5107i000+16.0000i00-1.1876+5.9751i
00000000+17.0648i00
Columns6through9
0000
00+16.0000i00
0000+17.0648i
-1.9422+10.5107i000
0-1.1876+5.9751i00
3.2242-15.8409i00-1.2820+5.5882i
02.8047-35.4456i-1.6171+13.6980i0
0-1.6171+13.6980i2.7722-23.3032i-1.1551+9.7843i
-1.2820+5.5882i0-1.1551+9.7843i2.4371-32.1539i
电压幅值:
1.04001.02501.02501.02580.99561.01271.02581.01591.0324
电压相角:
00.16200.0814-0.0387-0.0696-0.06440.06490.01270.0343
节点有功:
0.71641.63000.85000.0000-1.2500-0.9000-0.0000-1.0000-0.0000
节点无功:
0.27050.0665-0.10860.0000-0.5000-0.3000-0.0000-0.3500-0.0000
修正后的节点导纳矩阵Y:
Columns1through5
0-20.6944i000+17.3611i000-19.3333i000000-20.3982i000+17.3611i003.3074-39.3089i-1.3652+11.6041i000-1.3652+11.6041i3.8716-17.6627i000-1.9422+10.5107i000+16.0000i00-1.1876+5.9751i
00000000+17.0648i0
Columns6through9
0000
00+16.0000i00
0000+17.0648i
-1.9422+10.5107i000
0-1.1876+5.9751i00
4.1321-16.0184i00-1.2820+5.5882i
02.8047-35.4456i-1.6171+13.6980i0
0-1.6171+13.6980i3.7323-23.6669i-1.1551+9.7843i
-1.2820+5.5882i0-1.1551+9.7843i2.4371-32.1539i
节点阻抗矩阵Z的第4列:
0.0463+0.1252i
0.0329+0.0693i
0.0316+0.0707i
0.0552+0.1493i
0.0589+0.1204i
0.0562+0.1226i
0.0397+0.0838i
0.0416+0.0814i
0.0378+0.0845i
精确计算结果:
短路电流:
模值:
6.4459
相角:
-71.9365
节点电压模值:
0.18310.56870.54270.00000.14660.15060.45370.4463
400.4495
支路电流:
ijIij
140.5779-3.1264i271.3702-1.4433i390.64294-1.4808i45-0.77968+1.5248i46-0.6411+1.477i57-0.89528+1.6436i69-0.73353+1.5487i780.50734+0.10234i890.062766+0.056451i近似计算结果:
短路电流:
模值:
6.2838
相角:
-69.7198
节点电压模值:
0.16110.52140.5157
0.00000.182750.16750.42270.43480.4217
五、程序流程图
六、程序及输入文件input_data.xls文件:
powerflow_cal.m文件:
l=9;%支路数n=9;%节点数m=6;%PQ节点数
Yn=zeros(n);%初始化节点导纳矩阵Y
DATA1=xlsread('input_data.xls',1);%计算节点导纳矩阵Yfork=1:
l
i=DATA1(k,1);j=DATA1(k,2);R=DATA1(k,3);X=DATA1(k,4);B2=DATA1(k,5);
Yn(i,i)=Yn(i,i)+1i*B2+1/(R+1i*X);Yn(j,j)=Yn(j,j)+1i*B2+1/(R+1i*X);Yn(i,j)=Yn(i,j)-1/(R+1i*X);Yn(j,i)=Yn(j,i)-1/(R+1i*X);
end
disp('节点导纳矩阵Yn:
');
disp(Yn);
G=real(Yn);
B=imag(Yn);
DATA2=xlsread('input_data.xls',2);
P=zeros(1,n);
Q=zeros(1,n);
U=ones(1,n);
P(2:
n)=DATA2(2:
n,3);
Q(4:
n)=DATA2(4:
n,4);
U(1:
3)=DATA2(1:
3,5);
%设置节点电压初值
e
(1)=DATA2(1,5);
e(2:
n)=1.0;
f(1:
n)=0.0;
%设置迭代次数
t=0;
tmax=10;
whilet<=tmax
%计算f(x)
a(1:
n)=0.0;
c(1:
n)=0.0;
fori=2:
n
forj=1:
n
a(i)=a(i)+G(i,j)*e(j)-B(i,j)*f(j);c(i)=c(i)+G(i,j)*f(j)+B(i,j)*e(j);end
end
fori=2:
n
deltaP(i)=P(i)-e(i)*a(i)-f(i)*c(i);end
forj=4:
n
deltaQ(j)=Q(j)-f(j)*a(j)+e(j)*c(j);end
fork=2:
3
deltaU2(k)=U(k)*U(k)-e(k)*e(k)-f(k)*f(k);end
fx=[deltaP(2:
n)deltaQ(4:
n)deltaU2(2:
3)]';%计算雅克比矩阵J
fori=2:
n
forj=2:
n
ifi~=j
H(i,j)=-(G(i,j)*e(i)+B(i,j)*f(i));8
N(i,j)=B(i,j)*e(i)-G(i,j)*f(i);
else
H(i,j)=-a(i)-(G(i,i)*e(i)+B(i,i)*f(i));
N(i,j)=-c(i)+(B(i,i)*e(i)-G(i,i)*f(i));
end
end
end
fori=4:
n
forj=2:
n
ifi~=j
M(i,j)=B(i,j)*e(i)-G(i,j)*f(i);
L(i,j)=G(i,j)*e(i)+B(i,j)*f(i);
else
M(i,j)=c(i)+(B(i,i)*e(i)-G(i,i)*f(i));
L(i,j)=-a(i)+(G(i,i)*e(i)+B(i,i)*f(i));
end
end
end
fori=2:
3
forj=2:
n
ifi~=j
R(i,j)=0;
S(i,j)=0;
else
R(i,j)=-2*e(i);
S(i,j)=-2*f(i);
end
end
end
J=[H(2:
n,2:
n)N(2:
n,2:
n);M(4:
n,2:
n)L(4:
n,2:
n);R(2:
3,2:
n)S(2:
3,2:
n)];
ifmax(abs(fx))<0.0001
%输出结果
break;
else
%求解修正方程获得dx
dx=-J^(-1)*fx;
dx=dx';
e(2:
n)=e(2:
n)+dx(1:
n-1);
f(2:
n)=f(2:
n)+dx(n:
2*(n-1));
t=t+1;
end
end
ift>tmax
str='潮流计算不收敛';
disp(str);
else
a(1:
n)=0.0;
c(1:
n)=0.0;
fori=1:
n
forj=1:
n
a(i)=a(i)+G(i,j)*e(j)-B(i,j)*f(j);c(i)=c(i)+G(i,j)*f(j)+B(i,j)*e(j);end
end
fori=1:
n
U(i)=e(i)+1i*f(i);
amp(i)=abs(U(i));
arg(i)=angle(U(i));
P(i)=e(i)*a(i)+f(i)*c(i);
Q(i)=f(i)*a(i)-e(i)*c(i);
end
disp('电压幅值:
');
disp(amp);
disp('电压相角:
');
disp(arg);
disp('节点有功:
');
disp(P);
disp('节点无功:
');
disp(Q);
end
%计算短路电流
f=4;
zf=0.0;
%修正节点导纳矩阵
Xd=DATA2(1:
3,6);
E=DATA2(1:
3,7);
fori=1:
3
Ii(i)=E(i)/(1i*Xd(i));
end
Y=Yn;
fori=1:
3
Y(i,i)=Y(i,i)+1/(1i*Xd(i));
end
forj=4:
n
Y(j,j)=Y(j,j)+(-P(j)+1i*Q(j))/(U(j)*U(j));end
disp('修正后的节点导纳矩阵Y:
');
disp(Y);
Z=Y^(-1);
disp('节点阻抗矩阵Z的第4列:
');
disp(Z(:
4));
%精确计算
disp('精确计算结果:
');
U0=U;
If=U0(f)/(Z(f,f)+zf);
amp=abs(If);
arg=atand(imag(If)/real(If));
disp('短路电流:
');
disp('模值:
');
disp(amp);
disp('相角:
');
disp(arg);
fori=1:
n
U(i)=U0(i)-Z(i,f)*If;
amp=abs(U);
end
disp('节点电压模值:
');
disp(amp);
disp('支路电流:
');
str=['i''j''Iij'];disp(str);
fork=1:
l
i=DATA1(k,1);
j=DATA1(k,2);
r=DATA1(k,3);
x=DATA1(k,4);
z=r+1i*x;
I=(U(i)-U(j))/z;
str=[num2str(i)''num2str(j)disp(str);
end
%近似计算
disp('近似计算结果:
');
U0(1:
n)=1.0;
If=U0(f)/(Z(f,f)+zf);
amp=abs(If);
arg=atand(imag(If)/real(If));
disp('短路电流:
');
disp('模值:
');
disp(amp);
11num2str(I)];''
disp('相角:
');
disp(arg);
fori=1:
n
U(i)=U0(i)-Z(i,f)*If;amp=abs(U);
end
disp('节点电压模值:
');disp(amp);
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 电力系统 潮流 短路 电流 计算 程序 汇总