打印部分

学号：01203490308电力系统分析课程设计题目不对称短路电流计算学院自动化学院专业电气工程及其自动化班级电气1203班姓名廖超指导教师唐爱红2015年6月24日课程设计任务书学生姓名：廖超专业班级：电气1203班指导教师：唐爱红工作单位：武汉理工大学题目:不对称短路电流计算初始条件：如图所示电力系统，线路和变压器都并联运行。各线路的电阻可以忽略，架空线路为有钢质避雷线路。架空线路和电缆的正序电抗分别取0.4Ω/km和0.08Ω/km，计算图示各点的短路电流。要求完成的主要任务:（包括课程设计工作量及其技术要求，以及说明书撰写等具体要求）1.写出不对称短路电流的计算方法；2.用C或FORTRAN语言实现不对称短路电流的计算程序，调试程序并求出计算结果；3.用某种电力系统仿真工具搭建仿真模型，将结果与程序计算结果比较分析；4.撰写设计说明书；5.提供计算程序代码。时间安排：5月19日：领取任务书，分小组学习设计指导书；5月20日：分小组学习潮流计算、短路电流计算和稳定性计算方法及计算步骤；5月21-23日：分小组学习MATLAB、PSCAD/EMTDC、BPA、PSS/E、PSASP电力系统仿真软件；5月24日-27日：搭建仿真模型并调试；5月28日：对仿真结果进行分析，修改仿真模型；5月29日：撰写设计说明书；5月30日：答辩。指导教师签名：年月日系主任（或责任教师）签名：年月日本科生课程设计成绩评定表姓名性别专业、班级课程设计题目：课程设计答辩或质疑记录：成绩评定依据：序号项目满分得分备注1学习态度端正102基本概念、原理清晰203书写格式规范、设计内容条理清晰204能综合运用所学知识205设计内容正确、创新106阐述设计内容流畅、正确回答答辩问题20合计最终评定成绩（以优、良、中、及格、不及格评定）指导教师签字：年月日5计算步骤及实现各部分功能的程序5.1原始数据的的输入根据算法流程图，此程序应该首先进行变量的定义和初始数据的输入，这里的原始数据包括故障信息和基本的数据信息。#includestdio.h#defineSB100#defineUav115main(){intK,j,n,d,S;floatX0,X1,X2,X3,X4,X5,X6,X7,X8,X9,X10,X11,X12,Y[7][7],U[7],I[7],UK,SN,SK,P,Q,X,b,L,Y1,Y2,Y3,K1,K2,K3,IK,IM,SD,PN,Xd;K1=K2=1.05;K3=0.96;clear();prinf(pleaseinputthepatamentofgenerater:Xd,PN,cosa\n);scanf(%f,%f,%f,&Xd,&PN,&X);X1=Xd*SN/(PN/X);scanf(%f,%f,%f,&Xd,&PN,&X);X10=Xd*SN/(PN/X);prinf(X1=%f,X10=%f\n,X0,X10);prinf(pleaseinputtheputamenteroftramsfer:UKandSN\n1:);scanf(%f,%f,,&UK,&SN);X2=UK*SN/SN;prinf(\n2:\n);scanf(%f,%f,,&UK,&SN);X3=UK*SN/SN;prinf(\n3:\n);scanf(%f,%f,,&UK,&SN);X4=UK*SN/SN;prinf(X2=%f,X3=%f,X4=%f\n,X2,X3,X4);prinf(pleaseinputtheputamenterofload:P,Q\n);scanf(%f,%f,,&P,&Q);X8=1/((Q/100)*Q/100+(P/100)*(P/100))*Q/100;X11=0.35;X12=1.2;prinf(pleaseinputtheMofload:P,cosa\n);scanf(%f,%f,,&P,&X);X9=1/6.5*100/(P/X);prinf(X8=%f,X9=%f\n,X8,X9);prinf(pleaseinputtheputamenterofline:\n);prinf(forexample:X=0.39,b=02.92e-6,i=100);prinf(\n1:);scanf(%f,%f,,&X,&b,&L);X5=X*L*SN/(Uav*Uav);Y1=(double)1/2*b*L*(Uav*Uav)/100;prinf(\n2:);scanf(%f,%f,,&X,&b,&L);X6=X*L*SN/(Uav*Uav);Y2=(double)1/2*b*L*(Uav*Uav)/100;prinf(\n3:);scanf(%f,%f,,&X,&b,&L);X7=X*L*SN/(Uav*Uav);Y3=(double)1/2*b*L*(Uav*Uav)/100;prinf(X5=%f,X6=%f,X7=%f\n,X5,X6,X7);prinf(Y1=%f,Y2=%f,Y3=%f\n,Y1,Y2,Y3);prinf(Y1=Y1,Y2=Y2,Y3=Y3\n,Y1,Y2,Y3);prinf(pleaseinputtheK(7)toseeYarray:\n);scanf(&d,&K);for(n=0;nK;n++)for(j=0;jK;j++)Y[n][j]=0;Y[0][0]=1/X1+1/X2;Y[0][1]=-1/(X2*K1);Y[1][0]=Y[0][1];Y[1][1]=1/X5+Y1+1/(X2*K1*K1);Y[1][2]=-1/X5;Y[2][1]=Y[1][2];Y[2][2]=Y1+1/X5+Y2+1/X7+(1/2*X3*K2*K2);Y[2][3]=-1/(1/2*X3*K2);Y[2][4]=-1/X6;Y[2][5]=-1/X7;Y[3][2]=Y[2][3];Y[3][3]=1/(1/2*X3)+(1/2*X10);Y[4][2]=Y[2][4];Y[4][4]=Y2+1/X8+1/X6;Y[5][5]=Y3+1/X7+1/(X4*K3*K3);Y[5][6]=-1/(X4*K3);Y[6][5]=Y[5][6];Y[6][6]=1/X4;for(n=0;nK;n++)for(j=0;jK;j++){prinf(j%,2f\t,Y[n][j]);if(j==6)prinf(\n);}for(n=0;nK;n++)U[n]=0.0;prinf(pleaseinputthenumberofshortcuicit:K=I(i7)\n);scanf(&f,&d);U[d]=1;I[0]=0.0;for(n=0;nK;n++)for(j=0;jK;j++)I[n]=Y[n][j]*U[j];n=0;prinf(I[0]=%f\n,I[n++]);prinf(I[1]=%f\n,I[n++]);prinf(I[2]=%f\n,I[n++]);prinf(I[3]=%f\n,I[n++]);prinf(I[4]=%f\n,I[n++]);prinf(I[5]=%f\n,I[n++]);prinf(I[6]=%f\n,I[n++]);IK=I[d]*SN/(sqrt(3)*Uav);IM=2.55*IK;SK=1.732*IK*Uav;prinf(IK=%f,SD=%f,IM=%f,IK,SD,IM);prinf(pleaseinputtheK(7)toseeYarray:\n);scanf(&d,&K);for(n=0;nK;n++)for(j=0;jK;j++)Y[n][j]=0;Y[0][0]=1/X1+1/X2;Y[0][1]=-1/(X2*K1);Y[1][0]=Y[0][1];Y[1][1]=1/X5+Y1+1/(X2*K1*K1);Y[1][2]=-1/X5;Y[2][1]=Y[1][2];Y[2][2]=Y1+1/X5+Y2+1/X7+(1/2*X3*K2*K2);Y[2][3]=-1/(1/2*X3*K2);Y[2][4]=-1/X6;Y[2][5]=-1/X7;Y[3][2]=Y[2][3];Y[3][3]=1/(1/2*X3)+(1/2*X10);Y[4][2]=Y[2][4];Y[4][4]=Y2+1/X8+1/X6;Y[5][5]=Y3+1/X7+1/(X4*K3*K3);Y[5][6]=-1/(X4*K3);Y[6][5]=Y[5][6];Y[6][6]=1/X4;for(n=0;nK;n++)for(j=0;jK;j++){prinf(j%,2f\t,Y[n][j]);if(j==6)prinf(\n);}scanf(&d,&K);for(n=0;nK;n++)for(j=0;jK;j++)Y[n][j]=0;Y[0][0]=1/(X2*K1*K1)+1/(3*X5);Y[0][1]=-1/(3*X5);Y[1][0]=Y[0][1];Y[1][1]=1/(3*X5)+1/(3*X6)+1/(3*X7)+1/(1/2*X3*K2*K2);Y[1][2]=-1/(3*X6);Y[2][1]=Y[1][2];Y[2][2]=1/(3*X6)+Y2+1/X12;Y[1][3]=Y[3][1];Y[3][3]=1/(X4*K3*K3)+1/(3*X7);Y[3][1]=-1/(3*X7);Y[4][3]=Y[3][4];Y[4][4]=1/X4;for(n=0;nK;n++)for(j=0;jK;j++){prinf(j%,2f\t,Y[n][j]);if(j==6)prinf(\n);}}6MATLAB仿真6.1仿真模型图图6.1三相不对称短路仿真模型6.2仿真结果图6.2故障点A相电压正序分量波形图图6.2故障点A相电压负序分量波形图图6.2故障点A相电压零序分量波形图