![电力系统下课程设计-短路电流计算_第1页](http://file3.renrendoc.com/fileroot_temp3/2022-2/23/5508b848-80d7-4c7c-8e8c-762d09184419/5508b848-80d7-4c7c-8e8c-762d091844191.gif)
![电力系统下课程设计-短路电流计算_第2页](http://file3.renrendoc.com/fileroot_temp3/2022-2/23/5508b848-80d7-4c7c-8e8c-762d09184419/5508b848-80d7-4c7c-8e8c-762d091844192.gif)
![电力系统下课程设计-短路电流计算_第3页](http://file3.renrendoc.com/fileroot_temp3/2022-2/23/5508b848-80d7-4c7c-8e8c-762d09184419/5508b848-80d7-4c7c-8e8c-762d091844193.gif)
![电力系统下课程设计-短路电流计算_第4页](http://file3.renrendoc.com/fileroot_temp3/2022-2/23/5508b848-80d7-4c7c-8e8c-762d09184419/5508b848-80d7-4c7c-8e8c-762d091844194.gif)
![电力系统下课程设计-短路电流计算_第5页](http://file3.renrendoc.com/fileroot_temp3/2022-2/23/5508b848-80d7-4c7c-8e8c-762d09184419/5508b848-80d7-4c7c-8e8c-762d091844195.gif)
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、电力系统分析课程设计报告题目: 3G9bus短路电流计算 系 别 电气工程学院 专业班级 10级电气四班 学生姓名 学 号 指导教师 提交日期 2012年12月10日 目录一、设计目的2二、短路电流计算的基本原理和方法22.1电力系统节点方程的建立22.2利用节点阻抗矩阵计算短路电流3三、3G9bus短路电流在计算机的编程53.1、三机九节点系统53.3输出并计算结果12四总结14一、设计目的1.掌握电力系统短路计算的基本原理;2.掌握并能熟练运用一门计算机语言(MATLAB语言或FORTRAN或C语言或C+语言);3.采用计算机语言对短路计算进行计算机编程计算。二、短路电流计算的基本原理和方
2、法2.1电力系统节点方程的建立利用节点方程作故障计算,需要形成系统的节点导纳(或阻抗)矩阵。一般短路电流计算以前要作电力系统的潮流计算,假定潮流计算的节点导纳矩阵已经形成,在此基础上通过追加支路的方式形成电力短路电流计算的节点导纳矩阵YN。1)对发电机节点在每一发电机节点增加接地有源支路 与串联求短路稳态解: 求短路起始次暂态电流解: 一般情况下发电机定子绕组电阻忽略掉,并将与的有源支路转化成电流源与导纳并联的形式2)负荷节点的处理负荷节点在短路计一算中一般作为节点的接地支路,并用恒定阻抗表示,其数值由短路前瞬间的负荷功率和节点实际电压算出,即首先根据给定的电力系统运行方式制订系统的等值电路,
3、并进行各元件标么值参数的计算,然后利用变压器和线路的参数形成不含发电机和负荷的节点导纳矩阵 YN 。 3)如果短路故障发生在输电线路内,在短路点要增加新节点,将一条输电线分成两段,并形成短路电流计算的节点导纳矩阵Y,最后形成包括所有发电机支路和负荷支路的节点方程如下:一般,节点导纳矩阵 Y 阵与 YN 阶次相同。其差别只在于 YN 阵不含发电机和负荷的阻抗;节点注入电流向量 I中只有发电机端节点的电流不为零。有非零电流源注人的节点称为有源节点。系统中的同步调相机可按发电机处理。任进行起始次暂态电流计算时,大型同步电动机、感应电动机以及以电动机为主要成分的综合负荷,特别是在短路点近处的这些负荷,
4、必要时也可以用有源支路表示,并仿照发电机进行处理。必须指出,在计算机已普遍应用的情况下,如果有必要的话,只要能提供短路计算所需的准确的原始数据,对短路进行更精确的计算并不存在什么障碍。2.2利用节点阻抗矩阵计算短路电流假定系统中的节点f经过渡阻抗zf发生短路。过渡阻抗渡阻抗zf,不参与形成网络的节点导纳(或阻抗)矩阵。图6-3中方框内的有源网络代表系统正常状态的单相等值网络。现在我们保持故障处的边界条件不变,把网络的原有部分同故障支路分开(见图6-3)。容易看出,对于正常状态的网络而言,发生短路相当于在故障节点f增加了一个注人电流一 If (短路电流以流出故障点为正,节点电流则以注入为正)。因
5、此,网络中任一节点 i 的电压可表示为: (6-3)由式(6-3)可见,任一节点 i 的电压郁由两项叠加而成。第一项表示当 If=0时由网络内所有电源在节点i产生的电压,也就是短路前瞬间正常运行状态一 F 的节点电压,这是节点电压的正常分量,记为是电网的潮流解。第二项是当网络中所有电流源都断开,仅仅由短路电流If在节点i产生的电压,这就是节点电压的故障分童。上述两个分量的叠加,就等于发生短路后节点,的实际电压,即 (6-4)式中,是短路前故障点的正常电压;当i=f时,是故障节点 f 的自阻抗,也称输入阻抗。(边界条件方程)带入可得: 对于非变压器支路,令 k=1。从计算公式(6-7)和(6-8
6、)可以看到,式中所用到的阻抗矩阵元素都带有列标f。这就是说,如果网络在正常状态下的节点电压为已知,为了进行短路计算,只须利用节点阻抗矩阵中与故障点f对应的一列元素。因此,尽管是采用了阻抗型的节点方程,但是并不需要作出全部阻抗矩阵。在短路的实际计一算中,一般只需形成网络的节点导纳矩阵并根据具体要求,用第四章所讲的方法求出阻抗矩阵的某一列或某几列元素即可。在应用节点阻抗矩阵进行短路计算时,我们都将采用这种算法。1) 解潮流计算,2) 修正解潮流的YN形成Y3) 指定短路点f4) 计算节点阻抗矩阵第k列5) 计算短路电流:6) 计算节点电压: 7) 计算支路电流:8) 输出计算结果三、3G9bus短
7、路电流在计算机的编程3.1、三机九节点系统图1 三机九节点系统表1 九节点系统支路参数支路R(p.u.)X(p.u.)B/2(p.u.)1400.05761.02700.06251.03900.05861.0450.010.0850.088460.0170.0920.079570.0320.1610.153690.0390.170.179780.00850.0720.0745890.01190.10080.1045表2 九节点系统发电机参数发电机编号节点类型PG(p.u.)VG(p.u.)(p.u.)(p.u.)1V1.040.31.1372PV1.631.0250.31.2113PV0.85
8、1.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(N1,1)for i=1:busnumber Bus(i).name=N1(i); Bus(i).type=N2(i); Bus(i)
9、.no=i; Bus(i).Base_KV=N3(i); Bus(i).PG=N4(i); Bus(i).QG=N5(i); Bus(i).PL=N6(i); Bus(i).QL=N7(i); Bus(i).pb=N8(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(A1,1)for i=1:aclinenumber Acline(i).fbname=A
10、1(i); Acline(i).tbname=A2(i); Acline(i).Base_KV=A3(i); Acline(i).R=A4(i); Acline(i).X=A5(i); Acline(i).hB=A6(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 endfid=fopen('Trans
11、data.txt');T=textscan(fid, '%s %f %f %s %f %f %f %f')fclose(fid);tansnumber=size(T1,1)for i=1:tansnumber Trans(i).fbname=T1(i); Trans(i).fbBase_KV=T2(i); Trans(i).fbrated_KV=T3(i); Trans(i).tbname=T4(i); Trans(i).tbBase_KV=T5(i); Trans(i).tbrated_KV=T6(i); Trans(i).R=T7(i); Trans(i).X=T8
12、(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_KV2/Trans(i).fbBase_KV
13、2; Trans(i).X=tempx*Trans(i).X; Trans(i).R=tempx*Trans(i).R;end%N=0%Trans(1)%Trans(2)% for Y=G+jB matrixG,B,B2=FormYmatrix(Bus,busnumber,Acline,aclinenumber,Trans,tansnumber); %B:=B'B2:=B"dlmwrite('Gmatrix.txt', G, 'delimiter', 't','precision', 6);dlmwrite(
14、39;Bmatrix.txt', B, 'delimiter', 't','precision', 6);GBB2pauseJP,JQ=FormJPQmatrix(Bus,B,B2,busnumber);JPiJP=-inv(JP)JQiJQ=-inv(JQ)pause%maxiteration=0for i=1:busnumber NodeV(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).an
15、gle); 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%for nointer=1:10maxdP=1.;maxdQ=1.;epsilon=0.000001;noiteration=0;while (maxdP>epsilon)&(maxdP>epsilon)deltaP,deltaQ,maxdP,maxdQ=F
16、ormdPQvector(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; if noiteration>20 break endendfor i=1:busnumber Bus(i).V=NodeV(i); NodeV(i)=NodeV(i)*Bus(i).Base_KV; Bus(i).angle=Nodea(i); Nodea(
17、i)=Nodea(i)*180/pi;endnoiterationNodea=Nodea'NodeV=NodeV'Clear子函数%生成G、B矩阵function G,B,X=FormYmatrix(Bus,busnumber,Acline,aclinenumber,Trans,tansnumber)Y=zeros(busnumber);X=zeros(busnumber);for i=1:busnumberY(i,i)=Y(i,i)+Bus(i).pb*j;endfor i=1:aclinenumber f=Acline(i).fbno; t=Acline(i).tbno;
18、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/A
19、cline(i).X; endfor i=1:tansnumber f=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).k2; 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
20、; 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;for i=1:busnumber if Bus(i).type=1 for k=1:busnumber JQ(i,k)=0.; JQ(k,i)=0.; JP(i,k)=0.; JP(k,i)=0.; end JQ(i,i)=1.; JP(i,i)=1.; end
21、if Bus(i).type=3 for k=1:busnumber JQ(i,k)=0.; JQ(k,i)=0.; end JQ(i,i)=1.; end endend%计算偏节点PQ差量function deltaP,deltaQ,maxdP,maxdQ=FormdPQvector(Bus,NodeV,Nodea,dQGQL,dPGPL,B,G,busnumber) deltaQ=dQGQL; deltaP=dPGPL; maxdP=0.; maxdQ=0.; for i=1:busnumber if Bus(i).type=1 deltaQ(i)=0.; deltaP(i)=0.; en
22、d if Bus(i).type=3 deltaQ(i)=0.; %y1=0; %y2=0; y3=0; for k=1:busnumber if (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); end end deltaP(i)=deltaP(i)-y3*NodeV(i); %deltaP2(i)=(deltaP2
23、(i)-(y1*VX(i)+y2*VY(i)/Bus(i).V; end if Bus(i).type=2 %y1=0; %y2=0; y3=0; y4=0; for k=1:busnumber if (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
24、)-Nodea(k)-B(i,k)*cos(Nodea(i)-Nodea(k); end end deltaP(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; end if maxdP<abs(deltaP(i); maxdP=abs(deltaP(i); end if maxdQ<abs(deltaQ(i); maxdQ=abs(deltaQ(i); end deltaP(i)=deltaP(i)/NodeV(i); deltaQ(i)=deltaQ(i)/NodeV(i);endend3.3输出并计算结果由以上程序可以得出:(1)进行系统正常运行状态的潮流计算,求得:己知公式1-(Zif/(Zff+Zf))因为金属性短路时Zf0公式为1-(Zif/Zff)再根据公式可得,V1&(0) =0.2774 ; V2&(0) 0.6770; V3&(0)0.6770; V4&(0) 0;V5&(0)
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 稀土永磁项目可行性研究报告-高端稀土永磁行业放量可期
- 2024年越南纽扣称重传感器行业现状及前景分析2024-2030
- 2024年欧洲生物气溶胶实时在线监测系统市场主要企业市场占有率及排名
- 2024-2030年墨西哥可吸收止血剂市场前景分析
- 集成电路设计(第4版) 试题及答案 作业8
- 公司年终总结发言稿范文大全10篇完整版
- 端午假期安全教育
- 微机原理与接口技术(第3版)课件 第9章 可编程定时计数器8253-part1
- 框架剪力墙施工组织
- 浙教版八年级下册《6.1 反比例函数》同步练习
- 2024年典型事故案例警示教育手册15例
- 2024年江苏省高考生物试卷(解析版)
- 一年级下册竖式计算题大全(打印版)-
- 高中语文课文中常见的偏义复词
- JJF(苏)243-2021分布光度计校准规范-(现行有效)
- 开笔礼主持词(精选4篇)_开笔礼仪式主持词
- 紧急放行控制管理程序
- 网络的利与弊辩论赛正方材料
- 平桂区山林权属纠纷的成因及对策
- 药物制剂设备教学大纲
- 中国政区空白地图
评论
0/150
提交评论