1、课设报告信息工程系电力系统分析课程设计报告书题目: 电力系统短路故障的计算机算法程序设计 专 业: 电气工程及其自动化 班 级: 学 号: 学生姓名: 指导教师: 钟 健 伟 2012年3月6日 信息工程学院课程设计任务书学生姓名 学 号 成 绩设计题目电力系统短路故障的计算机算法程序设计设计内容电力系统故障的计算程序设计及编制和调试。设计要求1在对称短路计算、简单不对称短路计算中任选一种计算作为计算作业。2计算机语言自选。3设计、编制、调试出相关的通用计算程序。4输入输出数据一律以文件格式形成。(1)对称短路计算的输入输出数据 (供参考)输入数据可考虑为以下两个部分:数据1(可用变量表示):
2、节点数、支路数、故障节点数据2(可用二维数组表示):支路参数即 支路首端号、支路末端号、支路电抗值输出数据(可用二维数组表示)为:节点导纳矩阵或阻抗矩阵、故障节点、故障点电流、支路电流、各节点电压。(2)简单不对称短路计算的输入输出数据 (供参考)输入数据可考虑为以下两个部分:数据1(可用变量表示):节点数、支路数、故障节点数据2(可用二维数组表示):支路参数即支路首端号、支路末端号、支路正序电抗值、负序电抗值、零序电抗值输出数据(可用二维数组表示)为:故障节点、故障点电流、各支路电流、各节点电压的序分量和相分量。时间安排2月273月1号确认题目并找资料3月2号3月9号整理修改打印参考资料何仰
3、赞,电力系统分析(上),华中科技大学出版社,2001 1.课程设计目的及意义: 根据所给的电力系统,编制短路电流计算程序,通过计算机进行调试,最后完成一个切实可行的电力系统计算应用程序。通过自己设计电力系统计算程序使同学对电力系统分析有进一步理解,同时加强计算机实际应用能力的训练。所谓短路,是指电力系统正常运行情况以外的相与相之间或相与地之间的“短接”。在电力系统正常运行时,除中性点外,相与相或相与地之间是绝缘的。如果由于某种原因使其绝缘破坏而构成了通路,我们就称电力系统是发生了短路故障。在三相系统中,短路故障可分为两大类:即对称短路(三相短路)和不对称短路(两相短路、两相接地短路、单相接地短
4、路)。其中三相短路虽然发生的机会较少,但情况严重,又是研究其它短路的基础。所以我们先研究最简单的三相短路电流的暂态变化规律。2. 课程设计内容电力系统故障的计算程序设计及编制和调试。3.计算机计算的原理我使用的是第四代计算机语言的MATLAB,MATLAB求解下列方程,并求解矩阵A的特征值。Ax=b,其中:A= 32 13 45 67 23 79 85 12 43 23 54 65 98 34 71 35b= 1 2 3 4解为:x=Ab;设A的特征值组成的向量e,e=eig(A)。可见,MATLAB的程序极其简短。更为难能可贵的是,MATLAB甚至具有一定的智能水平,比如上面的解方程,MAT
5、LAB会根据矩阵的特性选择方程的求解方法,所以用户根本不用怀疑MATLAB的准确性。2)运算符丰富。由于MATLAB是用C语言编写的,MATLAB提供了和C语言几乎一样多的运算符,灵活使用MATLAB的运算符将使程序变得极为简短。3)MATLAB既具有结构化的控制语句(如for循环,while循环,break语句和if语句),又有面向对象编程的特性。4)程序限制不严格,程序设计自由度大。例如,在MATLAB里,用户无需对矩阵预定义就可使用。5)程序的可移植性很好,基本上不做修改就可以在各种型号的计算机和操作系统上运行。6)MATLAB的图形功能强大。在FORTRAN和C语言里,绘图都很不容易,
6、但在MATLAB里,数据的可视化非常简单。MATLAB还具有较强的编辑图形界面的能力。7)MATLAB的缺点是,它和其他高级程序相比,程序的执行速度较慢。由于MATLAB的程序不用编译等预处理,也不生成可执行文件,程序为解释执行,所以速度较慢。8)功能强大的工具箱是MATLAB的另一特色。MATLAB包含两个部分:核心部分和各种可选的工具箱。核心部分中有数百个核心内部函数。其工具箱又分为两类:功能性工具箱和学科性工具箱。功能性工具箱主要用来扩充其符号计算功能,图示建模仿真功能,文字处理功能以及与硬件实时交互功能。功能性工具箱用于多种学科。而学科性工具箱是专业性比较强的,如control,too
7、lbox,signl proceessing toolbox,commumnication toolbox等。这些工具箱都是由该领域内学术水平很高的专家编写的,所以用户无需编写自己学科范围内的基础程序,而直接进行高,精,尖的研究。9)源程序的开放性。开放性也许是MATLAB最受人们欢迎的特点。除内部函数以外,所有MATLAB的核心文件和工具箱文件都是可读可改的源文件,用户可通过对源文件的修改以及加入自己的文件构成新的工具箱。4.短路电流的构成及求解A. (1).周期分流电路求解:(2).非周期分量的求解:B.电机的过渡过程理论及短路电流的求解5例题的求解:例6-3的数据,线路的电阻和电容略去不
8、计,电压器的标么变比等于参数:z12=J0.105, z45=J0.184,z24=J0.08,z23=J0.065,z34=J0.05Y11=-J16.1905,Y12=J905238,Y21=J905238,Y22-J37.4084,Y23=J15.3846,Y24=J12.500,Y32=J15.3846,Y33=-J35.3846,Y34=J20.000,Y42=J15.5000,Y43=-J37.9348,Y45=J5.4348,Y54=J5.4348,Y55=-J9.9802 (1)电路图:(2)流程图: 输入数据形成节点导纳矩阵选择故障点If输入数据计算节点阻抗矩阵If列元素用公
9、式(6-10)计算短路电流If If=1/(Zff+zf)z电流If用公式(6-11)计算各点电压 Vi=1-Zif/( Zff+zf)用公式(6-9)计算指定支路的电路 Ipq=(Kvp-Vq)/zpq输出结果Y=0 -16.1905i 0 + 9.5238i 0 0 0 0 + 9.5238i 0 -37.4084i 0 +15.3846i 0 +12.5000i 0 0 0 +15.3846i 0 -35.3846i 0 +20.0000i 0 0 0 +12.5000i 0 +20.0000i 0 -37.9348i 0 + 5.4348i 0 0 0 0 + 5.4348i 0 -
10、9.9802i 采用4-3节的方法Z=0 + 0.1181i 0 + 0.0958i 0 + 0.0902i 0 + 0.0858i 0 + 0.0467i 0 + 0.0958i 0 + 0.1629i 0 + 0.1533i 0 + 0.1459i 0 + 0.0794i 0 + 0.0902i 0 + 0.1533i 0 + 0.1860i 0 + 0.1611i 0 + 0.0877i 0 + 0.0858i 0 + 0.1459i 0 + 0.1611i 0 + 0.1729i 0 + 0.0941i 0 + 0.0467i 0 + 0.0794i 0 + 0.0877i 0 + 0
11、.0941i 0 + 0.1515i 计算电流及网络中的电流分布If=Vi/Zii=1/Zii V1=V2=V3=V4=V5=1 (注意:Vi=0)V1=V1-Z1i*IfV2=V2-Z2i*IfV3=V3-Z3i*IfV4=V4-Z4i*IfV5=V5-Z5i*IfI54=(V5-V4)/z45I43=(V4-V3)/z43I23=(V2-V3)/z23I12=(V1-V2)/z12I24=(V2-V4)/z426程序及说明程序如下:jd=input(please input jd=);if jd=1Y=-16.1905j, 9.5238j, 0, 0 , 0; 9.5238j, 37.40
12、84j, 15.3846j, 12.5000j, 0; 0, 15.3846j, -35.3846j, 20.000j, 0; 0, 12.5000j, 20.000j, -37.9348j, 5.4348j; 0,0,0,5.4348j,-9.9802jinv(Y)I=1/0.1181iV1=0V2=1-0.0958i*IV3=1-0.0902i*IV4=1-0.0858i*IV5=1-0.0467i*II54=(V5-V4)/0.184iI43=(V4-V3)/0.05iI23=(V2-V3)/0.065iI12=(V1-V2)/0.105iI24=(V2-V4)/0.08ielseif
13、jd=2Y=-16.1905j, 9.5238j, 0, 0 , 0; 9.5238j, 37.4084j, 15.3846j, 12.5000j, 0; 0, 15.3846j, -35.3846j, 20.000j, 0; 0, 12.5000j, 20.000j, -37.9348j, 5.4348j; 0,0,0,5.4348j,-9.9802jinv(Y)I=1/0.1629iV1=1-0.0958i*IV2=0V3=1-0.1533i*IV4=1-0.1459i*IV5=1-0.0794i*II54=(V5-V4)/0.184iI43=(V4-V3)/0.05iI23=(V2-V3
14、)/0.065iI12=(V1-V2)/0.105iI24=(V2-V4)/0.08ielseif jd=3Y=-16.1905j, 9.5238j, 0, 0 , 0; 9.5238j, 37.4084j, 15.3846j, 12.5000j, 0; 0, 15.3846j, -35.3846j, 20.000j, 0; 0, 12.5000j, 20.000j, -37.9348j, 5.4348j; 0,0,0,5.4348j,-9.9802jinv(Y)I=1/0.1860iV1=1-0.0902i*IV2=1-0.1533i*IV3=0V4=1-0.1611i*IV5=1-0.08
15、77i*II54=(V5-V4)/0.184iI43=(V4-V3)/0.05iI23=(V2-V3)/0.065iI12=(V1-V2)/0.105iI24=(V2-V4)/0.08ielseif jd=4Y=-16.1905j, 9.5238j, 0, 0 , 0; 9.5238j, 37.4084j, 15.3846j, 12.5000j, 0; 0, 15.3846j, -35.3846j, 20.000j, 0; 0, 12.5000j, 20.000j, -37.9348j, 5.4348j; 0,0,0,5.4348j,-9.9802jinv(Y)I=1/0.1729iV1=1-
16、0.0858i*IV2=1-0.1459i*IV3=1-0.1611i*IV4=0V5=1-0.0941i*II54=(V5-V4)/0.184iI43=(V4-V3)/0.05iI23=(V2-V3)/0.065iI12=(V1-V2)/0.105iI24=(V2-V4)/0.08ielseif jd=5Y=-16.1905j, 9.5238j, 0, 0 , 0; 9.5238j, 37.4084j, 15.3846j, 12.5000j, 0; 0, 15.3846j, -35.3846j, 20.000j, 0; 0, 12.5000j, 20.000j, -37.9348j, 5.4
17、348j; 0,0,0,5.4348j,-9.9802jinv(Y)I=1/0.1515iV1=1-0.0467i*IV2=1-0.0794i*IV3=1-0.0877i*IV4=1-0.0941i*IV5=0I54=(V5-V4)/0.184iI43=(V4-V3)/0.05iI23=(V2-V3)/0.065iI12=(V1-V2)/0.105iI24=(V2-V4)/0.08i else c=errorEnd7.程序结果及分析节点1:please input jd=1 y11 =0 -16.1905i y12 =0 + 9.5238i y13 = 0y14 = 0y15 =0y21 =0
18、 + 9.5238iy22 =0 -37.4084iy23 =0 +15.3846iy24 =0 +12.5000iy25 =0y31 =0y32 =0 +15.3846iy33 =0 -35.3846iy34 =0 +20.0000iy35 =0y41 =0y42 =0 +12.5000iy43 =0 +20.0000iy44 =0 -37.9348iy45 =0 + 5.4348iy51 =0y52 =0y53 =0y54 =0 + 5.4348iy55 =0 - 9.9802iy0 = Columns 1 through5 0 -16.1905i 0 + 9.5238i 0 0 0 0
19、+ 9.5238i 0 -37.4084i 0 +15.3846i 0 +12.5000i 0 0 0 +15.3846i 0 -35.3846i 0 +20.0000i 0 0 0 +12.5000i 0 +20.0000i 0 -37.9348i 0 + 5.4348i 0 0 0 0 + 5.4348i 0 - 9.9802ians =Columns 1 through 4 0 + 0.1181i 0 + 0.0958i 0 + 0.0902i 0 + 0.0858i 0 + 0.0958i 0 + 0.1629i 0 + 0.1533i 0 + 0.1459i 0 + 0.0902i
20、0 + 0.1533i 0 + 0.1860i 0 + 0.1611i 0 + 0.0858i 0 + 0.1459i 0 + 0.1611i 0 + 0.1729i 0 + 0.0467i 0 + 0.0794i 0 + 0.0877i 0 + 0.0941iColumn 5 0 + 0.0467i 0 + 0.0794i 0 + 0.0877i 0 + 0.0941i 0 + 0.1515ii =0 - 8.4674iv1 =0v2 =0.1888v3 =0.2362v4 =0.2735v5 =0.6046i54 =0 - 1.7993ii43 =0 - 0.7451ii23 =0 + 0
21、.7295ii12 =0 + 1.7983ii24 =0 + 1.0584i节点2: please input jd=2 i =0 - 6.1387iv1 =0.4119v2 =0v3 =0.0589v4 =0.1044v5 =0.5126i54 =0 - 2.2186ii43 =0 - 0.9085ii23 =0 + 0.9066ii12 =0 - 3.9229ii24 =0 + 1.3045i节点3:please input jd=3 i =0 - 5.3763iv1 =0.5151v2 =0.1758v3 =0v4 =0.1339v5 =0.5285i54 =0 - 2.1447ii43
22、 =0 - 2.6774ii23 =0 - 2.7047ii12 =0 - 3.2309ii24 =0 - 0.5242i节点4:please input jd=4 i =0 - 5.7837iv1 =0.5038v2 =0.1562v3 =0.0682v4 =0v5 =0.4558i54 =0 - 2.4769ii43 =0 + 1.3650ii23 =0 - 1.3525ii12 =0 - 3.3105ii24 =0 - 1.9520i节点5:please input jd=5 i =0 - 6.6007iv1 =0.6917v2 =0.4759v3 =0.4211v4 =0.3789v5
23、 =0i54 = 0 + 2.0591ii43 =0 + 0.8449ii23 =0 - 0.8429ii12 = 0 - 2.0556ii24 =0 - 1.2129i8 .设计感悟: 在这次电力系统分析课程设计(电力系统短路故障的计算机算法程序设计)中我学到了很多。明白短路是电力系统中危害最严重的一种故障,短路计算是一项基础而又重要的工作。更进一步的熟悉了MATLAB的操作与使用,结合课本和图书馆查阅资料,对于电力系统短路故障的计算有了清晰的认识,无论是公式还是一些算法,在程序编写中都发挥的淋漓尽致,节点之间短路电流,节点电压,过渡电阻,自电阻等有了了解。感谢所有对我有过帮助的老师和同学!
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1