短路电流计算

课程设计报告

题目：短路电流计算

所在学院

专业班级

学生姓名

学生学号

指导教师

提交日期

目录

一、计算原理： (1)

二、计算流程： (1)

三、程序分析： (3)

四、短路电流计算主程序输入文件清单： (4)

五、导纳矩阵程序输入文件清单： (7)

六、整理计算结果： (8)

七、短路电流计算的要点： (12)

八、体会 (12)

附件一： (13)

一、计算原理： 基于节点阻抗矩阵用MATLAB 语言编写程序计算图1所示的算例系统的三相短路电流，当4号母线发生金属性三相短路时（z f =0），分别按照精确算法和近似算法计算短

路电流、系统中各节点电压以及网络中各支路的电流分布，并对两种情况下的计算结果进行比较。

近似算法是指：形成节点导纳矩阵时，所有节点的负荷都略去不计，短路前网络处于空载状态，短路前各节点电压均为1.0。

图1 三机九节点系统

二、计算流程：

计算流程框图(0)i V

i I (/i di

E x ''''=阻抗的接地支路表示；形成包括所有发电机支路和负荷支路的节点导纳矩阵中的发电机节点和负荷节点的自导纳上分别增加发电机导纳

利用，计算节点阻抗矩阵，从而得到阻抗矩阵中的第f列

计算短路电流

计算系统中各节点电压

计算各支路的短路电流

乘以IB得到有名值

解毕

注：计算短路电流以及各端电压请参见附录一。

三、程序分析：

为了计算短路电流，那么必须先得出计算导纳矩阵的程序，其中我们需要根据以下导纳矩阵的公式

变压器：

线路：

接地支路：

211;11;i ,i i ,i j ,j j ,j i ,j i ,j j ,i j ,i Y Y Y Y R jX k (R jX )Y Y Y Y k(R jX )k(R jX )⎧⎫←+←+⎪⎪++⎪⎪⎨⎬⎪⎪←-←-⎪⎪++⎩⎭

11;2211;i ,i i ,i j ,j j ,j i ,j i ,j j ,i j ,i B B Y Y j Y Y j R jX R jX Y Y Y Y R jX R jX ⎧⎫←++←++⎪⎪++⎪⎪⎨⎬⎪⎪←-←-⎪⎪++⎩⎭1i ,i i ,i Y Y R jX ⎧⎫←+⎨⎬+⎩

⎭

四、短路电流计算主程序输入文件清单：

Sbase_MVA=100.

fid=fopen('Nodedata.txt');

N=textscan(fid, '%s %u %d %f %f %f %f %f %f');

fclose(fid);

busnumber=size(N{1},1)

for i=1:busnumber

Bus(i).name=N{1}(i);

Bus(i).type=N{2}(i);

Bus(i).no=i;

Bus(i).Base_KV=N{3}(i);

Bus(i).PG=N{4}(i);

Bus(i).QG=N{5}(i);

Bus(i).PL=N{6}(i);

Bus(i).QL=N{7}(i);

Bus(i).pb=N{8}(i); %接地电抗（或电容）导纳标幺值Bus(i).V=1.0;

Bus(i).angle=0;

end

fid=fopen('Aclinedata.txt');

A=textscan(fid, '%s %s %f %f %f %f');

fclose(fid);

aclinenumber=size(A{1},1)

for i=1:aclinenumber

Acline(i).fbname=A{1}(i);

Acline(i).tbname=A{2}(i);

Acline(i).Base_KV=A{3}(i);

Acline(i).R=A{4}(i);

Acline(i).X=A{5}(i);

Acline(i).hB=A{6}(i);

for k=1:busnumber

if strcmp(Acline(i).fbname, Bus(k).name)

Acline(i).fbno=Bus(k).no;

end

if strcmp(Acline(i).tbname, Bus(k).name)

Acline(i).tbno=Bus(k).no;

end

end

end

fid=fopen('Transdata.txt');

T=textscan(fid, '%s %f %f %s %f %f %f %f');

fclose(fid);

tansnumber=size(T{1},1)

for i=1:tansnumber

Trans(i).fbname=T{1}(i);

Trans(i).fbBase_KV=T{2}(i);

Trans(i).fbrated_KV=T{3}(i);

Trans(i).tbname=T{4}(i);

Trans(i).tbBase_KV=T{5}(i);

Trans(i).tbrated_KV=T{6}(i);

Trans(i).R=T{7}(i);

Trans(i).X=T{8}(i);

for k=1:busnumber

if strcmp(Trans(i).fbname, Bus(k).name)

Trans(i).fbno=Bus(k).no;

end

if strcmp(Trans(i).tbname, Bus(k).name)

Trans(i).tbno=Bus(k).no;

end

end

Trans(i).k=Trans(i).tbrated_KV*Trans(i).fbBase_KV/Trans(i).fbrated_KV/Trans(i).tbBase_KV; tempx=Trans(i).fbrated_KV^2/Trans(i).fbBase_KV^2;

Trans(i).X=tempx*Trans(i).X;

Trans(i).R=tempx*Trans(i).R;

end

[G,B,B2]=FormYmatrix(Bus,busnumber,Acline,aclinenumber,Trans,tansnumber); %B:=B';B2:=B"

%dlmwrite('Gmatrix.txt', G, 'delimiter', '\t','precision', 6);

%dlmwrite('Bmatrix.txt', B, 'delimiter', '\t','precision', 6);

YN=j*B2 %按近似算法所得的节点导纳矩阵，YN不包括发电机导纳和负荷导纳的节点导纳矩阵pause

y=1.0/(j*0.15) %每一台发电机的次暂态电抗0.15

for i=1:3