22节点初步短路计算

最新推荐文章于 2021-12-21 17:10:03 发布

expect_bug

最新推荐文章于 2021-12-21 17:10:03 发布

阅读量905

点赞数

分类专栏：作业文章标签： matlab

本文链接：https://blog.csdn.net/expect_bug/article/details/95915937

版权

使用MATLAB进行22节点电力系统的短路计算，包括单相接地、两相、两相接地和三相短路的故障电流、电压及支路电流的详细计算过程。通过读取数据文件，构建节点导纳矩阵，并输出不同故障类型的电气参数。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

clear all;
fn=‘22节点.txt’; %数据文件

sys=DataRead(fn); %读取数据文件
Y1=sys.Y1; %正序
Y2=sys.Y2; %负序
Y0=sys.Y0; %零序
Vf0=1; %故障前短路点电压
pos=11; %短路点位置
Z1=full(inv(Y1)); %结点导纳矩阵：正序
Z2=full(inv(Y2)); %负序
Z0=full(inv(Y0)); %零序
Z1n=Z1(pos,pos); %短路点戴维宁阻抗：正序，负序，零序
Z2n=Z2(pos,pos);
Z0n=Z0(pos,pos);
%% 单相接地短路
Ifn=Vf0/(Z1n+Z2n+Z0n); %故障点故障电流：正序，负序，零序
Ifn1=Ifn;
Ifn2=Ifn;
Ifn0=Ifn;

Vfn1=Vf0-Z1nIfn1; %故障点故障电压：正序，负序，零序
Vfn2=-Z2nIfn2;
Vfn0=-Z0nIfn0;
V1=1+Z1(:,pos)(-Ifn1); %各节点电压：正序，负序，零序
V2=Z2(:,pos)(-Ifn2);
V0=Z0(:,pos)(-Ifn0);
a=exp(1j2pi/3);
T=[1 1 1;a^2 a 1;a a^2 1];
Vabc=[T*[V1,V2,V0]’]’; %结点电压

%支路电流
sys=dlmread(fn);
flag=find(sys(:,1)==0);
N=sys(1,1);
Branch=sys(flag(1)+1:flag(2)-1,1:6);

%计算各支路电流：正序，负序，零序
for K=1:size(Branch,1)
I=Branch(K,3);
J=Branch(K,4);
x1=Branch(K,5);
x0=Branch(K,6);
IS1=(V1(I)-V1(J))/(Z1(I,J));
IS2=(V2(I)-V2(J))/(Z2(I,J));
if x0~=0
IS0=(V0(I)-V0(J))/(1j*x0);
else
IS0=0;

end
BSC(K,:)=full([K,I,J,IS1,IS2,IS0]);

end

%输出结果
fprintf(’ 单相接地短路 \n\n’);
fprintf(’=\n’);
fprintf(’ 接地点故障分量 |\n’);
fprintf(’=\n’);
fprintf(’ 正序电流 | 负序电流 | 零序电流 |\n’);
fprintf(’=\n’);
fprintf(’ % 4.4f+j% 4.4f % 4.4f+j%4.4f % 4.4f+j% 4.4f \n’,…
real(Ifn1),imag(Ifn1),real(Ifn2),imag(Ifn2),real(Ifn0),imag(Ifn0));
fprintf(’=\n’);
fprintf(’ 正序电压 | 负序电压 | 零序电压 |\n’);
fprintf(’=\n’);
fprintf(’ % 4.4f+j% 4.4f % 4.4f+j% 4.4f % 4.4f+j% 4.4f \n’,…
real(Vfn1),imag(Vfn1),real(Vfn2),imag(Vfn2),real(Vfn0),imag(Vfn0));
fprintf(’=\n’);
fprintf(’\n\n\n’);
fprintf(’ 各支路电流 \n\n’);
fprintf(’=====================================================================================================\n’);
fprintf(’ 序号 | From BUS | To BUS | 正序 | 负序 | 零序 |\n’);
fprintf(’----------

最低0.47元/天解锁文章