电力系统下课程设计报告

老詢理N*穿广州学院GuangzhouCollegeofSouthChinaUniversityofTechnology《电力系统分析》课程MATCH_ word word文档格式规范word作业纸小票打印word模板word简历模板免费word简历 _1717202180289_2 报告 软件系统测试报告下载sgs报告如何下载关于路面塌陷情况报告535n,sgs报告怎么下载竣工报告下载 题目：3G9bus短路电流计算系别电气工程学院专业班级10级电气四班学生姓名学号指导教师提交日期2012年12月10日目录TOC\o"1-5"\h\zHYPERLINK\l"bookmark6"\o"CurrentDocument"一、设计目的3HYPERLINK\l"bookmark8"\o"CurrentDocument"二、短路电流计算的基本原理和 方法 快递客服问题件处理详细方法山木方法pdf计算方法pdf华与华方法下载八字理论方法下载 3HYPERLINK\l"bookmark10"\o"CurrentDocument"2.1电力系统节点方程的建立3HYPERLINK\l"bookmark20"\o"CurrentDocument"2.2利用节点阻抗矩阵计算短路电流4HYPERLINK\l"bookmark22"\o"CurrentDocument"三、3G9bus短路电流在计算机的编程6HYPERLINK\l"bookmark24"\o"CurrentDocument"3.1、三机九节点系统6HYPERLINK\l"bookmark28"\o"CurrentDocument"3.3输出并计算结果13四•总结15、设计目的掌握电力系统短路计算的基本原理；掌握并能熟练运用一门计算机语言（MATLA语言或FORTRA或C语言或C+h语言）；3.采用计算机语言对短路计算进行计算机编程计算。、短路电流计算的基本原理和方法2.1电力系统节点方程的建立利用节点方程作故障计算，需要形成系统的节点导纳（或阻抗）矩阵。一般短路电流计算以前要作电力系统的潮流计算，假定潮流计算的节点导纳矩阵已经形成，在此基础上通过追加支路的方式形成电力短路电流计算的节点导纳矩阵YN1）对发电机节点在每一发电机节点增加接地有源支路Ei与Z^RijXi串联求短路稳态解：Ei=EQj乙=R•jXqi求短路起始次暂态电流解：E=E「乙=R+jXi般情况下发电机定子绕组电阻忽略掉，并将巳与乙=R•jXi的有源支路转化成电流源1Y=G+Bi=———h=巳忆与导纳R+jXi并联的形式2）负荷节点的处理负荷节点在短路计一算中一般作为节点的接地支路，并用恒定阻抗表示，其数值由短路前瞬间的负荷功率和节点实际电压算出，即首先根据给定的电力系统运行方式制订系统的等值电路，并进行各元件标么值参数的计算，然后利用变压器和线路的参数形成不含发电机和负荷的节点导纳矩阵YN。ZLDk=RLDk+jXLDkVk2S_DkYLDk-GLDk'jBLDk二I¥1fi-1境电机和负荷曙负电跡的接人图c-m对称址跻井折3）如果短路故障发生在输电线路内，在短路点要增加新节点，将一条输电线分成两段，并形成短路电流计算的节点导纳矩阵Y,最后形成包括所有发电机支路和负荷支路的节点方程如下：YV=1一般，节点导纳矩阵丫阵与YN阶次相同。其差别只在于YN阵不含发电机和负荷的阻抗；节点注入电流向量I中只有发电机端节点的电流不为零。有非零电流源注人的节点称为有源节点。系统中的同步调相机可按发电机处理。任进行起始次暂态电流计算时，大型同步电动机、感应电动机以及以电动机为主要成分的综合负荷，特别是在短路点近处的这些负荷，必要时也可以用有源支路表示，并仿照发电机进行处理。必须指出，在计算机已普遍应用的情况下，如果有必要的话，只要能提供短路计算所需的准确的原始数据，对短路进行更精确的计算并不存在什么障碍。2.2利用节点阻抗矩阵计算短路电流假定系统中的节点f经过渡阻抗zf发生短路。过渡阻抗渡阻抗zf,不参与形成网络的节点导纳（或阻抗）矩阵。图6-3中方框内的有源网络代表系统正常状态的单相等值网络。现在我们保持故障处的边界条件不变，把网络的原有部分同故障支路分开（见图6-3）容易看出，对于正常状态的网络而言，发生短路相当于在故障节点f增加了一个注人电流一If（短路电流以流出故障点为正，节点电流则以注入为正）。因此，网络中任一节点iVi=HZjlj—Zifllf（i=1,2,..,f,.,n）的电压可表示为：jG（6-3）由式（6-3）可见，任一节点i的电压郁由两项叠加而成。第一项表示当lf=0时由网络内所有电源在节点i产生的电压，也就是短路前瞬间正常运行状态一F的节点电压，（0）这是节点电压的正常分量，记为Vi是电网的潮流解。第二项是当网络中所有电流源都断开，仅仅由短路电流If在节点i产生的电压，这就是节点电压的故障分童。上述两个分量的叠加，就等于发生短路后节点，的实际电压，即V=V⑼—ZifIlf（i=1,2,..,仁.,n）Vf=Vf（0）—ZifIf（6-4）（o）Z式中，i是短路前故障点的正常电压；当i=f时，"是故障节点f的自阻抗，也称V_zIV—V（0）zI输入阻抗。Vf如f（边界条件方程）带入Vf=Vf-Zif1f可得：IlfVf⑼ZffZfV—「VVf⑼(i=1,2,..,f,.,n)Si行潇曲端网皓对于非变压器支路，令图的5玄路电流计誅ZffZfkVp-Vqpq:Zpqk=1。从计算 公式 小学单位换算公式大全免费下载公式下载行测公式大全下载excel公式下载逻辑回归公式下载 （6-7）和（6-8）可以看到，式中所用到的阻抗矩阵元素都带有列标f。这就是说，如果网络在正常状态下的节点电压为已知，为了进行短路计算，只须利用节点阻抗矩阵中与故障点f对应的一列元素。因此，尽管是采用了阻抗型的节点方程，但是并不需要作出全部阻抗矩阵。在短路的实际计一算中，一般只需形成网络的节点导纳矩阵•并根据具体要求，用第四章所讲的方法求出阻抗矩阵的某一列或某几列元素即可。在应用节点阻抗矩阵进行短路计算时，我们都将采用这种算法。1）解潮流计算，呼（i二1,2，..，f，.，n）2）修正解潮流的YN形成丫3）指定短路点fY2iY2Y22Yn2(-第f元素,(k=1,2,川,n)出\\\5)计算短路电流：IfNzf6)计算节点电压：Vi二Vi(0)—ZfIlf(i=1,2,..,f,.,n)IkVp-Vq7)计算支路电流：1pq_zzpq8)输出计算结果三、3G9bus短路电流在计算机的编程3.1、三机九节点系统27Q图1三机九节点系统3表1九节点系统支路参数支路R(p.u.)X(p.u.)B/2(p.u.)1~400.05761.02~700.06251.03~900.05861.04~50.010.0850.0884~60.0170.0920.0795~70.0320.1610.1536~90.0390.170.1797~80.00850.0720.07458~90.01190.10080.1045表2九节点系统发电机参数发电机编号节点类型Pg(p.u.)Vg(p.u.)X；(p.u.)E“(p.u.)1V1.040.31.1372PV1.631.0250.31.2113PV0.851.0250.31.043表3九节点系统负荷参数节点编号节点类型Pi(p.u.)Qi(p.u.)4PQ005PQ1.250.56PQ0.90.37PQ008PQ10.359PQ00主程序3.2程序设计主函数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_KVA2/Trans(i).fbBase_KVA2;Trans(i).X=tempx*Trans(i).X;Trans(i).R=tempx*Trans(i).R;end%N=0%Trans(1)%Trans(2)%forY=G+jBmatrix[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);GBB2pause[JP,JQ]=FormJPQmatrix(Bus,B,B2,busnumber);JPiJP=-inv(JP)JQiJQ=-inv(JQ)pause%maxiteration=0fori=1:busnumberNodeV(i)=Bus(i).V;Nodea(i)=Bus(i).angle;VX(i)=Bus(i).V*cos(Bus(i).angle);VY(i)=Bus(i).V*sin(Bus(i).angle);dQGQL(i)=Bus(i).QG-Bus(i).QL;dPGPL(i)=Bus(i).PG-Bus(i).PL;endNodeV=NodeV'Nodea=Nodea'%VX=VX'%VY=VY'dQGQL=dQGQL'dPGPL=dPGPL'pause%fornointer=1:1OmaxdP=1.;maxdQ=1.;epsilon=0.000001;noiteration=0;while(maxdP>epsilon)&(maxdP>epsilon)[deltaP,deltaQ,maxdP,maxdQ]=FormdPQvector(Bus,NodeV,Nodea,dQGQL,dPGPL,B,G,busnumber);deltaP;deltaQ;maxdP;maxdQ;da=iJP*deltaP;dV=iJQ*deltaQ;Nodea=Nodea+da;NodeV=NodeV+dV;noiteration=noiteration+1;ifnoiteration>20breakendendfori=1:busnumberBus(i).V=NodeV(i);NodeV(i)=NodeV(i)*Bus(i).Base_KV;Bus(i).angle=Nodea(i);Nodea(i)=Nodea(i)*180/pi;endnoiterationNodea=Nodea'NodeV=NodeV'Clear子函数%生成GB矩阵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)+Bus(i).pb*j;endfori=1:aclinenumberf=Acline(i).fbno;t=Acline(i).tbno;Y(f,f)=Y(f,f)+Acline(i).hB*j+1/(Acline(i).R+Acline(i).X*j);Y(t,t)=Y(t,t)+Acline(i).hB*j+1/(Acline(i).R+Acline(i).X*j);Y(f,t)=Y(f,t)-1/(Acline(i).R+Acline(i).X*j);Y(t,f)=Y(t,f)-1/(Acline(i).R+Acline(i).X*j);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+Trans(i).X*j);Y(t,t)=Y(t,t)+1/(Trans(i).R+Trans(i).X*j)/Trans(i)."2;Y(f,t)=Y(f,t)-1/(Trans(i).R+Trans(i).X*j)/Trans(i).k;Y(t,f)=Y(t,f)-1/(Trans(i).R+Trans(i).X*j)/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%生成JP、JQ矩阵function[JP,JQ]=FormJPQmatrix(Bus,B,B2,busnumber)JP=B;JQ=B2;fori=1:busnumberifBus(i).type==1fork=1:busnumberJQ(i,k)=0.;JQ(k,i)=0.;JP(i,k)=0.;JP(k,i)=0.;endJQ(i,i)=1.;JP(i,i)=1.;endifBus(i).type==3fork=1:busnumberJQ(i,k)=0.;JQ(k,i)=0.;endJQ(i,i)=1.；endendendfunction[deltaP,deltaQ,maxdP,maxdQ]=FormdPQvector(Bus,NodeV,Nodea,dQGQL,dPGPL,B,G,busnumber)deltaQ=dQGQL;deltaP=dPGPL;maxdP=0.;maxdQ=0.;fori=1:busnumberifBus(i).type==1deltaQ(i)=0.;deltaP(i)=0.;endifBus(i).type==3deltaQ(i)=0.;%y1=0;%y2=0;y3=0;fork=1:busnumberif(B(i,k)~=0|G(i,k)~=0)%y仁y1+(G(i,k)*VX(k)-B(i,k)*VY(k));%y2=y2+(G(i,k)*VY(k)+B(i,k)*VX(k));y3=y3+NodeV(k)*(G(i,k)*cos(Nodea(i)-Nodea(k))+B(i,k)*sin(Nodea(i)-Nodea(k)));endenddeltaP(i)=deltaP(i)-y3*NodeV(i);%deltaP2(i)=(deltaP2(i)-(y1*VX(i)+y2*VY(i)))/Bus(i).V;endifBus(i).type==2%y1=0;%y2=0;y3=0;y4=0;fork=1:busnumberif(B(i,k)~=0|G(i,k)~=0)%y1=y1+(G(i,k)*VX(k)-B(i,k)*VY(k));%y2=y2+(G(i,k)*VY(k)+B(i,k)*VX(k));y3=y3+NodeV(k)*(G(i,k)*cos(Nodea(i)-Nodea(k))+B(i,k)*sin(Nodea(i)-Nodea(k)));y4=y4+NodeV(k)*(G(i,k)*sin(Nodea(i)-Nodea(k))-B(i,k)*cos(Nodea(i)-Nodea(k)));endenddeltaP(i)=deltaP(i)-y3*NodeV(i);%deltaP2(i)=(deltaP2(i)-(y1*VX(i)+y2*VY(i)))/Bus(i).V;deltaQ(i)=deltaQ(i)-y4*NodeV(i);%deltaQ2(i)=(deltaQ2(i)-(y1*VY(i)-y2*VX(i)))/Bus(i).V;endifmaxdP