电力系统分析课程设计湖北民大.docx
- 文档编号:30040918
- 上传时间:2023-08-04
- 格式:DOCX
- 页数:45
- 大小:342KB
电力系统分析课程设计湖北民大.docx
《电力系统分析课程设计湖北民大.docx》由会员分享,可在线阅读,更多相关《电力系统分析课程设计湖北民大.docx(45页珍藏版)》请在冰豆网上搜索。
电力系统分析课程设计湖北民大
2011-2012学年度下学期电力系统分析课程设计
电力系统短路故障的计算机
算法程序设计
姓名
学号03094
班级五
指导教师钟健伟
1课程设计说明与方案论证
1.1设计说明
在电力系统运行过程中,时常会发生故障,且大多是短路故障,一般指短路和断线,分为简单故障和复杂故障,一切不正常的相与相之间或相与地之间(对于中性点接地的系统)发生通路的情况即为短路。
电流剧增:
设备发热增加,若短路持续时间较长,可能使设备过热甚至损坏;由于短路电流的电动力效应,导体间还将产生很大的机械应力,致使导体变形甚至损坏。
电压大幅度下降,对用户影响很大。
当短路发生地点离电源不远而持续时间又较长时,并列运行的发电机可能失去同步,破坏系统运行的稳定性,造成大面积停电,这是短路最严重的后果。
发生不对称短路时,三相不平衡电流会在相邻的通讯线路感应出电动势,影响通讯。
短路通常分为三相短路、单相接地短路、两相短路和两相接地短路。
其中三相短路为对称短路,后三者为不对称短路。
电力运行经验指出单相接地短路占大多数,因此分析与计算不对称短路具有非常重要意义,鉴于学习过程的规律和重视基础原理,在这里我做出对三相对称短路的计算机算法。
三相短路在电力系统中虽然发生的概率很小,但是一旦发生产生的影响非常严重,电力系统中设备的动稳定和热稳定校验都要以三相短路电流为依据,因此我选择做是电力系统三相短路的计算程序,觉得更有意义。
实现了三相短路计算的计算机实现。
本程序输入含有网络拓扑信息和电器元件的参数,输出三相短路的短路电流和短路后网络各节点的电压参数和各支路的电流参数。
并以文件的形式保存,方便用户的读取,打印。
实现了设计的要求。
应用计算机进行电力系统计算,首先要建立电力系统相应的数学模型,如建立电力系统等效网络和节点方程;其次是运用合理的计算方法计算结果,如三角分解法,LDU分解法等;第三则是选择合适的计算机语言编制计算程序,实现程序的准确性,提高精确度。
本设计利用节点方程来实现电力系统三相短路计算。
首先,根据给定的电力系统运行方式制定系统的一相等值网络,计算出各元器件的参数,其次根据网络方程利用变压器和线路参数形成不含发电机和负荷的节点导纳矩阵Y;然后,利用系统形成的节点导纳矩阵求解系统阻抗矩阵Z;第三则是应用节点阻抗矩阵计算短路电流。
电力系统短路电流的实际工程计算中,许多实际问题的解决并不需要十分精确的结果。
因此采用近似计算法,在建立系统节点方程时,采用网络模型和参数做简化处理,等值电路的制定通常将发电机作为含源支路,表示为节点注入电流源和节点电抗的形式;系统中存在同步发电机时,按发电机支路处理;同时,忽略发电机、变压器和输电线路的电阻,不计线路的电容,略去变压器的励磁电流,负荷忽略不计;在标幺参数计算时,选取各级平均额定电压作为基准电压,忽略元件的额定电压和平均电压的差别,认为变压器变比的标幺值为1;此外,还假定所有的发电机的电势具有相同的相位,所有的元件仅用电抗表示,避免复数运算。
三相对称短路为金属性短路。
本课程设计的关键部分是利用系统节点导纳矩阵形成节点阻抗矩阵,本程序采用LDU分解法,然后利用节点导纳矩阵和阻抗矩阵的乘积为单位矩阵列写方程,节点导纳矩阵经LDU分解后,回代求解系统阻抗矩阵。
利用节点阻抗矩阵进行短路计算,首先利用故障前的节点注入电流和系统节点阻抗矩阵求取故障前系统的各节点电压;然后,利用金属短路时短路点残压为0,对地电阻为0,计算短路点电流;最后,利用故障前节点电压、节点阻抗和短路点电流求得短路后各节点电压和各支路电流。
1.2实现方式
matlab是一种解释性语言,集成化,专门用于和数学相关的工程计算,数值运算matlab是很简便强大,它提供了基本的数学算法,例如矩阵运算、数值分析算法,MATLAB集成了2D和3D图形功能,以完成相应数值可视化的工作,并且提供了一种交互式的高级编程语言——M语言,利用M语言可以通过编写脚本或者函数文件实现用户自己的算法。
典型的应用于:
数学和计算;算法开发;建模、模拟和原形化;数据分析,探索和可视化;科学与工程制图;应用开发,包括图形用户界面的建立等广大领域。
但是由于客观原因,对matlab的操作不是很熟练,还处于基础学习阶段,不利于设计的顺利完成与实现。
C语言是面向对象的基础语言。
更多的用于底层函数开发,软件开发,单片机控制等,C语言是一种计算机程序设计语言。
它既具有高级语言的特点,又具有汇编语言的特点。
它可以作为工作系统设计语言,编写系统应用程序,也可以作为应用程序设计语言,编写不依赖计算机硬件的应用程序。
它的应用范围广泛,具备很强的数据处理能力,不仅仅是在软件开发上,而且各类科研都需要用到C语言,适于编写系统软件,三维,二维图形和动画。
具体应用比如单片机以及嵌入式系统开发。
它把高级语言的基本结构和语句与低级语言的实用性结合起来。
C语言可以像汇编语言一样对位、字节和地址进行操作,而这三者是计算机最基本的工作单元。
结构式语言代码及数据的分隔化,即程序的各个部分除了必要的信息交流外彼此独立。
这种结构化方式可使程序层次清晰,便于使用、维护以及调试。
C语言以函数形式提供给用户,这些函数可方便的调用,并具有多种循环、条件语句控制程序流向,从而使程序完全结构化。
功能齐全。
具有各种各样的数据类型,并引入了指针概念,可使程序效率更高。
而且计算功能、逻辑判断功能也比较强大,可以实现决策目的的游戏。
适用范围大。
适合于多种操作系统,如Windows、DOS、UNIX等等;也适用于多种机型。
C语言对编写需要硬件进行操作的场合,明显优于其它高级语言,有一些大型应用软件也是用C语言编写的。
主观上分析,由于有C语言的学习基础,对C语言的熟练程度更深,和兴趣所在,决定选择C语言实现设计目的。
2总体设计
2.1功能实现
图1-1对称短路计算原理图
2.1程序流程图
3.1.1主程序流图
图1-2主流程图
3详细设计
3.1主要程序流程图
3.1.1创建新数据
图1-3创建新数据
3.1.2调取数据程序
图1-4调取数据
3.1.4改变节点子程序
图1-6改变节点子程序
3.1.3计算子程序
图1-5计算子程序
3.2程序算法及变量说明
3.2.1网络中任意一点电流计算公式:
(1.1)
在不要求精确计算的场合,可以不计负荷电流的影响。
在形成节点导纳矩阵时,所有节点的负荷都略去不计,短路前网络处于空载,各节点电压的正常分量的标幺等于1。
(1.2)
网络任意节点电压计算公式:
(1.3)
3.2.2
本程序中将电路的支路分为三类:
电抗支路、电流支路、发电机直流,数据输入时分开输入格
式如下:
电抗支路:
支路数、节点首端好、节点末端好,支路电抗值。
发电机支路:
发电机台数、发电机接入点,发电机电势、发电机之路电抗。
电容支路:
支路数、电容节点号,电容支路电抗值。
程序中的变量定义如下:
FILE*fp;//文件指针
floatIf;//短路电流
floatV0[10];//节点电压初值
floatIn[10];//注入电流
floatV[10];//节点电压
floatI[10];//支路电流
intcount[2][10];//数组用于储存支路号
floatz[10][10];//储存支路电抗值
floatyjs[10][10];//存放计算导纳矩阵中的中间值
floatye[10];//储存发电机支路导纳
floatY[10][10];//节点导纳矩阵
floatZ[10][10];//阻抗矩阵
floatdk;//电抗值
inte1[10];
floatyc[10];//节点对地电抗
floatycjs[10];
inte[10];//发电机接入点
floatE[10];
floatIg[10];
floatEmid[10];
floatze[10];//储存发电机支路导纳
电力系统分析课程设计
6
intm,n,a,f,i=0,j=0,sd,md,butn,times,timesc;/*n是节点数,a是支路数,f是短路点,i为行标,j为列
标*/
1、文件读入程序如下:
printf("请输入系统的节点数,支路数,故障节点\n");
scanf("%d%d%d",&n,&a,&f);
fp=fopen(name,"w+");
fprintf(fp,"系统参数如下所示:
\n\n");
fprintf(fp,"1.系统的节点数为%d\n2.系统的支路数为%d\n3.系统的短路点为%d\n\n",n,a,f);
fprintf(fp,"支路电抗参数如下\n");
2、计算主程序如下:
(1)、If=V0[f-1]/Z[f-1][f-1];(短路电流计算)
for(i=0;i V[i]=V0[i]-Z[i][f-1]*If; (2)、for(i=0;i { for(j=0;j Y[i][i]=Y[i][i]-yjs[i][j]; Y[i][i]=Y[i][i]-ye[i]-ycjs[i]; } for(i=0;i for(j=i+1;j { Y[i][j]=yjs[i][j]; Y[j][i]=yjs[i][j]; } 3.3程序代码及注释 2.4.1具体源程序代码及其注释见附录 3.4实例测算 【例6-3】在如图1-7所示的电力系统中分别在节点1和节点5接入发电机支路,其标幺值参数为: = =1.0, =j0.15, =j0.22。 在节点3发生三相短路,计算短路电流及网络中的电流分 布。 线路的电阻和电容略去不计,变压器的标幺变比等于1。 各元件参数的标幺值如下: =j0.105, =1.05, =j0.184, =0.03+j0.08, =0.024+j0.065, =0.018+j0.05, = =j0.02, = =j0.016, = =j0.013。 图1-7电力系统等值网络 图1-8三相短路等值网络 由以上网络图可得计算过程如下: (1)形成输入文件: 系统参数如下所示: 1.系统的节点数为5 2.系统的支路数为5 3.系统的短路点为3 支路电抗参数如下 =0.1050 =0.0650 =0.0800 =0.1840 =0.0500 发电机接入点及其电抗如下: 共二台发电机进入系统 =0.1500 =0.2200 =1.0000 =1.0000 共0个节点有接地电容 (2)计算输出如下图 计算结果如下: 1、节点导纳矩阵: -16.19059.5238 9.5238-37.408415.384612.5000 15.3846-35.384620.0000 12.500020.0000-37.93485.4348 5.4348-9.9802 2、节点阻抗矩阵: 0.11810.09580.09020.08580.0467 0.09580.16290.15330.14590.0794 0.09020.15330.18600.16110.0877 0.08580.14590.16110.17290.0941 0.04670.07940.08770.09410.1515 3、短路节点是: 3 4、短路电流: If=5.3767 5、各节点电压: V=0.5152V2=0.1758V3=0V4=0.1336V5=0.5282 6、发电机支路电流如下: Ig[1][0]=-3.2321 Ig[5][0]=-2.1445 7、各支路电流: I12=-3.2321I23=-2.7046I24=-0.5275I43=-2.6720I54=-2.1445 图1-9对称短路计算结果 图1-10改变短路节点后运算结果 【习题6-10】如图3-9所示网络,略去负荷,试用节点阻抗矩阵法求节点5发生三相短路时,短路 点的短路电流及线路L-2、L-3的电流。 已知各元件参数如下: 发电机G-1: =120MVA,x”d=0.23;G-2: =60MVA,x”d=0.14 变压器T-1: =120MVA,Vs=10.5%;T-2: =60MVA,Vs=10.5% 线路参数x1=0.4Ω/km,b1=2.8×10-6S/km 线路长度L-1: 120km,L-2: 80km,L-3: 70km 取 =120MVA, = 图1-11电力系统网络图 图1-12三相短路等效网络 (1)、形成输入文件: 系统参数如下所示: 1.系统的节点数为: 5 2.系统的支路数为: 5 3.系统的短路点为: 5 支路电抗参数如下: z35=0.29036z45=0.25406z34=0.43554z13=0.1050z24=0.2100 发电机接入点及其电抗如下: 共2台发电机进入系统: z1=0.2300z2=0.2800 e1=1.0500e2=1.0500 共3个节点有接地电容: Y3=-32.372299Y4=-34.077000Y5=-43.185200 (2)、得到的输出文件: 计算结果如下: 1、节点导纳矩阵: -13.87169.5238 -8.33334.7619 9.5238-15.23292.29603.4440 4.76192.2960-10.96463.9361 3.44403.9361-7.3569 2、节点阻抗矩阵: 0.18190.06160.16000.10790.1326 0.06160.20780.08980.15360.1242 0.16000.08980.23300.15710.1932 0.10790.15360.15710.26890.2174 0.13260.12420.19320.21740.3427 3、短路节点是: 5 4、短路电流: If=3.1263 5、各节点电压: V1=0.6471V2=0.6723V3=0.4632V4=0.3890V5=0 6、发电机支路电流如下: Ig[10]=-1.7516 Ig[20]=-1.3491 7、各支路电流: =-1.5954 =-1.5310 =-0.1705 =-1.7516 =-1.3491 【习题6-11】电力系统等值电路如图1-13所示,支路阻抗的标幺值已注明图中。 (1)形成节点导纳矩阵(或节点阻抗矩阵),并用以计算节点3的三相短路电流。 (2)另选一种方法计算短路电流,并用以验证 (1)的计算结果。 网络图 经过诺顿变换的网络 图1-13等值电路网络 由以上网络可得输入文件如下: 系统参数如下所示: 1.系统的节点数为: 3 2.系统的支路数为: 2 3.系统的短路点为: 3 支路电抗参数如下: z12=0.2000z23=0.1000 发电机接入点及其电抗如下: 共2台发电机进入系统: z1=0.2000z3=0.130435 e1=1.0500e3=0.913043 共1个节点有接地电容: Y2=-2.000000 (2)、得到的输出文件: 计算结果如下: 1、节点导纳矩阵: -10.00005.0000 5.0000-14.500010.0000 10.0000-17.6667 2、节点阻抗矩阵: 0.13940.07890.0446 0.07890.15770.0893 0.04460.08930.1071 3、短路节点是: 3 4、短路电流: If=9.1875 5、各节点电压: V1=0.6344V2=0.2188V3=0 6、发电机支路电流如下: Ig[1][0]=-2.0781 Ig[3][0]=-7.0000 7、各支路电流: =-2.0781 =-2.1875 通过计算结果与算例答案的对比,计算结果与算例标准答案一样,说明了程序计算方法和思路的准确性。 说明本次课程设计基本达到了设计的要求。 4总结 在历经多天的学习和实践之后,终于完成了本次课程设计,学习到了相当多的宝贵经验是实际操作技巧。 更加使我熟练了在电力系统设计过程中基本知识的重要性,学习应该先注重基础,再加强实践实际处理能力,才能达到能够将理论转化为实践的本领,能够胜任未来工作中的任务和需要。 另外编程能力也是一个强有力的本领应该多实践,多做事。 加强对专业软件的了解与学习,能够在发展中不断充实自己增长自己的才干,能够利用多种手段与方法解决问题难题,以便能够以最优最好的方式解决挑战性难题,增加解决困难的几率。 本系统中,只解决了电力系统中对称短路的问题,对不对称短路还没有解决,在电力系统运行过程中,时常会发生故障,且大多是短路故障。 短路通常分为三相短路、单相接地短路、两相短路和两相接地短路。 其中三相短路为对称短路,后三者为不对称短路。 电力运行经验指出单相接地短路占大多数,因此分析与计算不对称短路具有非常重要意义。 分析计算不对称短路方法很多,目前实际最常用的方法是对称分量法。 而以对称分量法为核心的计算方法又可有解析法和计算机程序算法等,本论文的主要工作即介绍这两种计算方法。 解析法,是将微分方程代数化、暂态分析稳态化、不对称转化为对称并叠加完成不对称故障的分析与计算。 计算机程序算法是在形成三个序网的节点导纳矩阵后,对其应用高斯消去法求得故障端点等值阻抗,根据故障类型选用相应公式计算各序电流、电压,进而合成三相电流、电压。 所以今后可以讲不对称短路引入设计范畴。 参考文献 (1)《电力系统分析》(上册)华中科技大学出版社何仰赞 (2)《电力系统故障的计算机辅助分析》重庆大学出版社米麟书等 (3)《电力系统故障分析》清华大学出版社周荣光 (4)《短路电流实用计算方法》电力工业出版社西安交通大学等 附录 主程序及注释(这里只给出主要计算程序): #include #include #include 声明外部函数 externvoidcreate(); externvoidload(); externvoidquit(); externvoidchange(charname[20]); externvoidcalculate(charname[20]) floatIf;//短路电流 floatV0[10];//节点电压初值 floatIn[10];//注入电流 floatV[10];//节点电压 floatI[10];//支路电流 intcount[2][10];//数组用于储存支路号 floatz[10][10];//储存支路电抗值 floatyjs[10][10];//存放计算导纳矩阵中的中间值 floatye[10];//储存发电机支路导纳 floatY[10][10];//节点导纳矩阵 floatZ[10][10];//阻抗矩阵 floatdk;//电抗值 inte1[10]; floatyc[10];//节点对地电抗 floatycjs[10]; inte[10];//发电机接入点 floatE[10]; floatIg[10]; floatEmid[10]; floatze[10];//储存发电机支路导纳 intm,n,a,f,g,i=0,j=0,sd,md,butn,times,timesc;/*n是节点数,a是支路数,f是短路点,i为行标,j为列标 for(i=0;i<10;i++) {yc[i]=0; ycjs[i]=0; E[i]=0; Emid[i]=0;} for(i=0;i<10;i++)//初始赋 { V0[i]=0; for(j=0;j<10;j++) { z[i][j]=0; Y[i][j]=0; } } /******************************读取系统参数**************************/ fp=fopen(name,"r"); fscanf(fp,"系统参数如下所示: \n\n"); fscanf(fp,"1.系统的节点数为%d\n2.系统的支路数为%d\n3.系统的短路点为%d\n\n",&n,&a,&f); /**************读电抗*************/ fscanf(fp,"支路电抗参数如下\n");
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 电力系统 分析 课程设计 湖北