短路电流计算

PAGE华南理工大学广州学院电气工程学院专业课程设计报告PAGEPAGE8课程设计报告题目：短路电流计算 所在学院专业班级学生姓名学生学号指导教师提交日期华南理工大学广州学院电气工程学院电力系统分析课程设计报告PAGE0目录一、计算原理： 1二、计算流程： 1三、程序分析： 3四、短路电流计算主程序输入文件清单： 4五、导纳矩阵程序输入文件清单： 7六、整理计算结果： 8七、短路电流计算的要点： 12八、体会 12附件一： 13一、计算原理：基于节点阻抗矩阵用MATLAB语言编写程序计算图1所示的算例系统的三相短路电流，当4号母线发生金属性三相短路时（zf=0），分别按照精确算法和近似算法计算短路电流、系统中各节点电压以及网络中各支路的电流分布，并对两种情况下的计算结果进行比较。近似算法是指：形成节点导纳矩阵时，所有节点的负荷都略去不计，短路前网络处于空载状态，短路前各节点电压均为1.0。图1三机九节点系统二、计算流程：将发电机表示为电流源()和导纳()的并联组合；节点负荷用恒阻抗的接地支路表示；形成包括所有发电机支路和负荷支路的节点导纳矩阵Y，即在YN中的发电机节点和负荷节点的自导纳上分别增加发电机导纳将发电机表示为电流源()和导纳()的并联组合；节点负荷用恒阻抗的接地支路表示；形成包括所有发电机支路和负荷支路的节点导纳矩阵Y，即在YN中的发电机节点和负荷节点的自导纳上分别增加发电机导纳和负荷导纳；形成不含发电机和负荷的节点导纳矩阵YN；进行系统正常运行状态的潮流计算，求得利用利用，计算节点阻抗矩阵，从而得到阻抗矩阵中的第f列计算短路电流计算系统中各节点电压计算各支路的短路电流乘以IB得到有名值解毕注：计算短路电流以及各端电压请参见附录一。三、程序分析：为了计算短路电流，那么必须先得出计算导纳矩阵的程序，其中我们需要根据以下导纳矩阵的公式变压器：线路：接地支路：四、短路电流计算主程序输入文件清单：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)fori=1:busnumberBus(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;endfid=fopen('Aclinedata.txt');A=textscan(fid,'%s%s%f%f%f%f');fclose(fid);aclinenumber=size(A{1},1)fori=1:aclinenumberAcline(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);fork=1:busnumberifstrcmp(Acline(i).fbname,Bus(k).name)Acline(i).fbno=Bus(k).no;endifstrcmp(Acline(i).tbname,Bus(k).name)Acline(i).tbno=Bus(k).no;endendendfid=fopen('Transdata.txt');T=textscan(fid,'%s%f%f%s%f%f%f%f');fclose(fid);tansnumber=size(T{1},1)fori=1:tansnumberTrans(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);fork=1:busnumberifstrcmp(Trans(i).fbname,Bus(k).name)Trans(i).fbno=Bus(k).no;endifstrcmp(Trans(i).tbname,Bus(k).name)Trans(i).tbno=Bus(k).no;endendTrans(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不包括发电机导纳和负荷导纳的节点导纳矩阵pausey=1.0/(j*0.15)%每一台发电机的次暂态电抗0.15fori=1:3YN(i,i)=YN(i,i)+y;%将发电机导纳加入节点导纳矩阵end%荷导节点开路YN%包括发电机导纳节点导纳矩阵pauseZ=inv(YN)%Z节点阻抗矩阵%4号节点三相短路If=(Bus(4).V*cos(Bus(4).angle)+j*Bus(4).V*sin(Bus(4).angle))/Z(4,4)%母线对地短路电流标幺值/(0-9.1239)i/puIB=Sbase_MVA/sqrt(3)/230%0.251kAIfymzh=abs(If*IB)%母线对地短路电流有名值2.2903kAfori=1:busnumberV0(i)=Bus(i).V*cos(Bus(i).angle)+j*Bus(i).V*sin(Bus(i).angle);Vf(i)=V0(i)-Z(i,4)*If;I(i)=Vf(i)/Z(i,4);I(i)=abs(I(i)*IB);endpauseV0%111111111Vf%0.27750.67690.67710.00000.18740.19350.54220.54590.5509Clear五、导纳矩阵程序输入文件清单：function[G,B,X]=FormYmatrix(Bus,busnumber,Acline,aclinenumber,Trans,tansnumber)Y=zeros(busnumber);X=zeros(busnumber);fori=1:busnumberY(i,i)=Y(i,i)+1/j*Bus(i).pb;endfori=1:aclinenumberf=Acline(i).fbno;t=Acline(i).tbno;Y(f,f)=Y(f,f)+j*Acline(i).hB+1/(Acline(i).R+j*Acline(i).X);Y(t,t)=Y(t,t)+j*Acline(i).hB+1/(Acline(i).R+j*Acline(i).X);Y(f,t)=Y(f,t)-1/(Acline(i).R+j*Acline(i).X);Y(t,f)=Y(t,f)-1/(Acline(i).R+j*Acline(i).X);X(f,f)=X(f,f)-1/Acline(i).X;X(t,t)=X(t,t)-1/Acline(i).X;X(f,t)=1/Acline(i).X;X(t,f)=1/Acline(i).X;endfori=1:tansnumberf=Trans(i).fbno;t=Trans(i).tbno;Y(f,f)=Y(f,f)+1/(Trans(i).R+j*Trans(i).X);Y(t,t)=Y(t,t)+1/[(Trans(i).R+j*Trans(i).X)/Trans(i).k^2];Y(f,t)=Y(f,t)-1/[(Trans(i).R+j*Trans(i).X)/Trans(i).k];Y(t,f)=Y(t,f)-1/[(Trans(i).R+j*Trans(i).X)/Trans(i).k];X(f,f)=X(f,f)-1/Trans(i).X;X(t,t)=X(t,t)-1/Trans(i).X;X(f,t)=1/Trans(i).X;X(t,f)=1/Trans(i).X;endG=real(Y);B=imag(Y);end六、整理计算结果：Sbase_MVA=100busnumber=9aclinenumber=6tansnumber=3YN=Columns1through50-17.3611i000+17.3611i000-16.0000i000000-17.0648i000+17.3611i000-39.9954i0+11.7647i0000+11.7647i0-17.9759i0000+10.8696i000+16.0000i000+6.2112i00000000+17.0648i00Columns6through9000000+16.0000i000000+17.0648i0+10.8696i00000+6.2112i000-16.7519i000+5.8824i00-36.1001i0+13.8889i000+13.8889i0-23.8095i0+9.9206i0+5.8824i00+9.9206i0-32.8678iy=0-6.6667iYN=Columns1through50-24.0278i000+17.3611i000-22.6667i000000-23.7315i000+17.3611i000-39.9954i0+11.7647i0000+11.7647i0-17.9759i0000+10.8696i000+16.0000i000+6.2112i00000000+17.0648i00Columns6through9000000+16.0000i000000+17.0648i0+10.8696i00000+6.2112i000-16.7519i000+5.8824i00-36.1001i0+13.8889i000+13.8889i0-23.8095i0+9.9206i0+5.8824i00+9.9206i0-32.8678iZ=Columns1through50+0.0988i0+0.0256i0+0.0256i0+0.0792i0+0.0644i0+0.0256i0+0.0963i0+0.0281i0+0.0354i0+0.0487i0+0.0256i0+0.0281i0+0.0963i0+0.0354i0+0.0369i0+0.0792i0+0.0354i0+0.0354i0+0.1096i0+0.0891i0+0.0644i0+0.0487i0+0.0369i0+0.0891i0+0.1378i0+0.0639i0+0.0367i0+0.0494i0+0.0884i0+0.0758i0+0.0363i0+0.0739i0+0.0398i0+0.0502i0+0.0690i0+0.0360i0+0.0594i0+0.0546i0+0.0498i0+0.0617i0+0.0356i0+0.0391i0+0.0753i0+0.0492i0+0.0514iColumns6through90+0.0639i0+0.0363i0+0.0360i0+0.0356i0+0.0367i0+0.0739i0+0.0594i0+0.0391i0+0.0494i0+0.0398i0+0.0546i0+0.0753i0+0.0884i0+0.0502i0+0.0498i0+0.0492i0+0.0758i0+0.0690i0+0.0617i0+0.0514i0+0.1412i0+0.0520i0+0.0590i0+0.0687i0+0.0520i0+0.1047i0+0.0842i0+0.0554i0+0.0590i0+0.0842i0+0.1227i0+0.0760i0+0.0687i0+0.0554i0+0.0760i0+0.1048iIf=0-9.1239iIB=0.2510Ifymzh=2.2903Z=Columns1through50+0.0988i0+0.0256i0+0.0256i0+0.0792i0+0.0644i0+0.0256i0+0.0963i0+0.0281i0+0.0354i0+0.0487i0+0.0256i0+0.0281i0+0.0963i0+0.0354i0+0.0369i0+0.0792i0+0.0354i0+0.0354i0+0.1096i0+0.0891i0+0.0644i0+0.0487i0+0.0369i0+0.0891i0+0.1378i0+0.0639i0+0.0367i0+0.0494i0+0.0884i0+0.0758i0+0.0363i0+0.0739i0+0.0398i0+0.0502i0+0.0690i0+0.0360i0+0.0594i0+0.0546i0+0.0498i0+0.0617i0+0.0356i0+0.0391i0+0.0753i0+0.0492i0+0.0514iColumns6through90+0.0639i0+0.0363i0+0.0360i0+0.0356i0+0.0367i0+0.0739i0+0.0594i0+0.0391i0+0.0494i0+0.0398i0+0.0546i