三相短路故障计算

三相短路故障计算

1.原始数据：如图表示有五个节点的系统结构接线图，网络参数在图中标出。在简化的

短路电流计算中，假定发电机电势的标么值等于1，试求：（1）节点1、2、3三相短路时的各支路电流及各节点的电压；（2）在支路1——2的中点发生三相短路时的各支炉电流及各节点的电压；（3）断开支路1——2时，节点发生三相短路时的各支路电流及各节点电压。

其中： 120.25x j =，130.35x j =，230.3x j =，240.015x j =，350.03x j =，''40.02d x j =，''50.02d x j =

2.程序框图：

3.结果显示：

(1)在各节点处发生短路时的值1.2.3节点分别短路时的各值：阻抗矩阵：

0 + 0.1664i 0 + 0.0205i 0 + 0.0208i 0 + 0.0117i 0 + 0.0083i 0 + 0.0205i 0 + 0.0307i 0 + 0.0061i 0 + 0.0175i 0 + 0.0025i 0 + 0.0208i 0 + 0.0061i 0 + 0.0412i 0 + 0.0035i 0 + 0.0165i 0 + 0.0117i 0 + 0.0175i 0 + 0.0035i 0 + 0.0186i 0 + 0.0014i 0 + 0.0083i 0 + 0.0025i 0 + 0.0165i 0 + 0.0014i 0 + 0.0186i

短路点电流：

各节点电压：

各支路电流：

(2) 1-2支路中线短路

阻抗矩阵：

0.1664i 0.0205i 0.0208i 0.0117i 0.0083i 0.0934i 0.0205i 0.0307i 0.0061i 0.0175i 0.0025i 0.0256i 0.0208i 0.0061i 0.0412i 0.0035i 0.0165i 0.0135i 0.0117i 0.0175i 0.0035i 0.0186i 0.0014i 0.0146i 0.0083i 0.0025i 0.0165i 0.0014i 0.0186i 0.0054i 0.0934i 0.0256i 0.0135i 0.0146i 0.0054i 0.1220i

用上述方法可以得其他各值

短路电流I f=08.1957i

各节点电压V:

0.2341 0.7903 0.8898 0.8802 0.9559 0

各支路电流Izl:

0 - 2.5057i 0 +31.5804i 0 +16.1963i 0 + 5.1223i 0 + 4.0107i 0 -30.8900i

(3) 1-2断线

阻抗矩阵：

0 + 0.3935i 0 + 0.0045i 0 + 0.0435i 0 + 0.0026i 0 + 0.0174i 0 + 0.0045i 0 + 0.0318i 0 + 0.0045i 0 + 0.0182i 0 + 0.0018i 0 + 0.0435i 0 + 0.0045i 0 + 0.0435i 0 + 0.0026i 0 + 0.0174i 0 + 0.0026i 0 + 0.0182i 0 + 0.0026i 0 + 0.0190i 0 + 0.0010i 0 + 0.0174i 0 + 0.0018i 0 + 0.0174i 0 + 0.0010i 0 + 0.0190i 同理求取其余各值

短路节点电流If: 0 + 2.5413i

各节点电压V:

0 0.9884 0.8894 0.9934 0.9558

各支路电流Izl:

0 +20.4438i 0 -21.7822i 0 + 0.2723i 0 + 3.8119i 0

4.程序清单

%B=[1 2 0.25i 0 0;1 3 0.35i 0 0;2 3 0.3i 0 0;2 4 0.015i 0 1./0.02i;3 5 0.03i 0 1./0.02i]

%B=[1 6 1./2*0.25i 0 0;1 3 0.35i 0 0;2 3 0.3i 0 0;2 4 0.015i 0 1./0.02i;3 5 0.03i 0 1./0.02i;2 6 1./2*0.25i 0 0]

%B=[1 3 0.35i 0 0;2 3 0.3i 0 0;2 4 0.015i 0 1./0.02i;3 5 0.03i 0 1./0.02i] % ½Úµãµ¼ÄÉÕóµÄÇóÈ¡

clear

n=input('½Úµã¸öÊý:n=');

nl=input('֧·Êý:nl=');

B=input('֧·²ÎÊý:B=');

Yn=zeros(n);%YÕóÇåÁã

for i=1:nl

p=B(i,1);q=B(i,2);

y=B(i,3);s=B(i,4);

t=B(i,5);

if s==0

Yn(p,p)=Yn(p,p)+1./y;

Yn(q,q)=Yn(q,q)+1./y+t;

Yn(p,q)=Yn(p,q)-1./y;

Yn(q,p)=Yn(q,p)-1./y;

end

end%ÇóÈ¡µ¼ÄÉÕó

Zn=inv(Yn); %ÇóÈ¡×迹Õó

f=input('ÊäÈë¶Ì·½ÚµãºÅ;f=');%Çó¶Ì·µãµçÁ÷

disp('×迹¾ØÕóZn:');

disp(Zn);

Zn(f,f)=Zn(f,f);

If(f)=-1./Zn(f,f);

disp('¶Ì·½ÚµãµçÁ÷If:');

disp(If(f));

for i=1:n %Çó¸÷½Úµãµçѹ

V(i)=1-(-Zn(i,f))*If(f);

end

disp('¸÷½ÚµãµçѹV:');

disp(V);

for i=1:nl %Çó֧·µçÁ÷

p=B(i,1); q=B(i,2);

k=-(Zn(p,f)-Zn(q,f))./Zn(p,q);

Izl(i)=k*1./Zn(f,f);

end

disp('¸÷֧·µçÁ÷Izl:');

disp(Izl);