电力系统分析课设-C语言.doc
《电力系统分析课设-C语言.doc》由会员分享,可在线阅读,更多相关《电力系统分析课设-C语言.doc(25页珍藏版)》请在咨信网上搜索。
1、电力系统分析课程设计目录1课程设计说明12编程语言的选择及理由13程序设计23。1主体流程图23。2详细流程图33。2。1创建系统(create)流程图33。2。2加载系统函数(load)流程图33.2。3计算子函数(calculate)流程图43.2.4改变短路节点(change)流程图43.3数据及变量说明53.4程序代码及注释63.5测试算例64设计体会13参考文献13附录13i1 课程设计说明本课程设计根据电力系统分析课程中所讲的基于节点方程的三相短路计算的原理和方法。采用c语言编程 。实现了三相短路计算的计算机实现。本程序输入含有网络拓扑信息和电器元件的参数,输出三相短路的短路电流和
2、短路后网络各节点的电压参数和各支路的电流参数.并以文件的形式保存,方便用户的读取,打印.实现了设计的要求。三相短路在电力系统中虽然发生的概率很小,但一旦发生产生的影响非常严重,电力系统中设备的动稳定和热稳定校验都要以三相短路电流为依据,因此我选择做是电力系统三相短路的计算程序,觉得更有意义.应用计算机进行电力系统计算,首先要建立电力系统相应的数学模型,如建立电力系统等效网络和节点方程;其次是运用合理的计算方法计算结果,如三角分解法,LDU分解法等;第三则是选择合适的计算机语言编制计算程序,实现程序的准确性,提高精确度.本设计利用节点方程来实现电力系统三相短路计算。首先,根据给定的电力系统运行方
3、式制定系统的一相等值网络,计算出各元器件的参数,其次根据网络方程利用变压器和线路参数形成不含发电机和负荷的节点导纳矩阵Y;然后,利用系统形成的节点导纳矩阵求解系统阻抗矩阵Z;第三则是应用节点阻抗矩阵计算短路电流。电力系统短路电流的实际工程计算中,许多实际问题的解决并不需要十分精确的结果.因此采用近似计算法,在建立系统节点方程时,采用网络模型和参数做简化处理,等值电路的制定通常将发电机作为含源支路,表示为节点注入电流源和节点电抗的形式;系统中存在同步发电机时,按发电机支路处理;同时,忽略发电机、变压器和输电线路的电阻,不计线路的电容,略去变压器的励磁电流,负荷忽略不计;在标幺参数计算时,选取各级
4、平均额定电压作为基准电压,忽略元件的额定电压和平均电压的差别,认为变压器变比的标幺值为1;此外,还假定所有的发电机的电势具有相同的相位,所有的元件仅用电抗表示,避免复数运算。三相对称短路为金属性短路.本课程设计的关键部分是利用系统节点导纳矩阵形成节点阻抗矩阵,本程序采用LDU分解法,然后利用节点导纳矩阵和阻抗矩阵的乘积为单位矩阵列写方程,节点导纳矩阵经LDU分解后,回代求解系统阻抗矩阵.利用节点阻抗矩阵进行短路计算,首先利用故障前的节点注入电流和系统节点阻抗矩阵求取故障前系统的各节点电压;然后,利用金属短路时短路点残压为0,对地电阻为0,计算短路点电流;最后,利用故障前节点电压、节点阻抗和短路
5、点电流求得短路后各节点电压和各支路电流.2 编程语言的选择及理由能够实现程序编写的语言有多种,本课程设计中采用C语言进行电力系统对称短路计算程序的编写.个人对c语言和metlab相对较熟悉。之所以选择C语言编程,从主观上说,是因为我学过C语言课程,相对于其他的编程语言而言,C语言的应用更为熟悉和了解。另外,从客观上说,C语言是目前世界上最为流行、使用最为广泛的高级程序设计语言.,C的效率MATLAB比高多了,MATLAB虽然语法简单,但执行效率不高.MATLAB是用于特定的方面的,比如说矩阵运算方面、DSP就很有优势,而C语言不同.很多软件的底层都是C编写的。MATLAB提供的是一种基于解释的
6、语言,虽然也是一种高级语言,但是还是很简单的,它的目标在于用户方便,开发简单。而实际上这些方便也是有代价的,编译器承受了这些代价,程序执行效率方面也有损失。另一方面MATLAB这么做也是基于其专用于科学计算的目的来考虑的,无法应用于通用的程序设计。c语言这么做保持了其编译的速度,便于大型程序设计,执行效率也很高,因此我选择用c语言编写程序 。3 程序设计3.1 主体流程图开始欢迎使用用户选择操作,按键123进入创建系统子函数(create)进入加载系统子函数(load)进入退出子函数进入计算子函数(calculate)改短路点子函数(change)是否改路点?结束是图31 主体流程图3.2 详
7、细流程图3.2.1 创建系统(create)流程图输入文件名输入电抗支路参数结束输入节点数、支路数、故障节点输入发电机支路参数输入电容支路参数图3-2 create函数流程图3.2.2 加载系统函数(load)流程图输入系统文件名输出查看系统参数选择进入下一个子函数图33 节点导纳矩阵形成流程图3.2.3 计算子函数(calculate)流程图 文件读入电路参数求取LDU分解的因子表矩阵下三角求取DUZj向量求解完毕?否是结束求取节点阻抗矩阵列向量根据支路参数建立节点导纳矩阵求取短路点电流求取节点电压求取支路电流建立文件存储计算结果结束图34 计算子函数流程图3.2.4 改变短路节点(chan
8、ge)流程图输入短路点从新读取文件中的参数选择进入下一个子函数计算图35 change函数流程图3.3 数据及变量说明本程序中将电路的支路分为三类:电抗支路、电流支路、发电机直流,数据输入时分开输入格式如下:电抗支路:支路数、节点首端好、节点末端好,支路电抗值.发电机支路:发电机台数、发电机接入点,发电机电势、发电机之路电抗。电容支路:支路数、电容节点号,电容支路电抗值.程序中的变量定义如下: FILE fp; /文件指针 float If; /短路电流 float V010; /节点电压初值 float In10; /注入电流 float V10; /节点电压 float I10; /支路电
9、流int count210; /数组用于储存支路号float z1010; /储存支路电抗值float yjs1010; /存放计算导纳矩阵中的中间值float ye10; /储存发电机支路导纳float Y1010; /节点导纳矩阵float Z1010; /阻抗矩阵float dk; /电抗值int e110;float yc10; /节点对地电抗float ycjs10;int e10; /发电机接入点float E10; float Ig10;float Emid10;float ze10; /储存发电机支路导纳int m,n,a,f,i=0,j=0,sd,md,butn,times,
10、timesc;/n是节点数,a是支路数,f是短路点,i为行标,j为列标*/1、文件读入程序如下: printf(请输入系统的节点数,支路数,故障节点n); scanf(”%d%dd,&n,a,f); fp=fopen(name,”w+); fprintf(fp,”系统参数如下所示:nn”); fprintf(fp,1.系统的节点数为dn2。系统的支路数为%dn3。系统的短路点为%dnn,n,a,f); fprintf(fp,”支路电抗参数如下n);2、计算主程序如下: (1)、 If=V0f-1/Zf-1f1;(短路电流计算) for(i=0;in;i+) Vi=V0i-Zif1If; (2)
11、、for(i=0;in;i+)(导纳矩阵计算) for(j=0;jn;j+) Yii=Yiiyjsij; Yii=Yiiyeiycjsi; for(i=0;in;i+) for(j=i+1;jn;j+) Yij=yjsij; Yji=yjsij; 3.4 程序代码及注释 具体源程序代码及其注释见附录。3.5 测试算例【例6-3】在如图35所示的电力系统中分别在节点1和节点5接入发电机支路,其标幺值参数为:。在节点3发生三相短路,计算短路电流及网络中的电流分布。线路的电阻和电容略去不计,变压器的标幺变比等于1。各元件参数的标幺值如下:图3-6 电力系统等值网络图图37 三相短路时的等值网络图由以
12、上网络图可得计算过程如下:(1)、形成输入文件:系统参数如下所示:1。系统的节点数为52。系统的支路数为53。系统的短路点为3支路电抗参数如下z12=0.1050z23=0。0650z24=0.0800z43=0.0500z54=0.1840发电机接入点及其电抗如下:共2台发电机进入系统z1=0。1500z5=0。2200e1=1.0000e5=1.0000共0个节点有接地电容(2)、计算后得到的输出文件:计算结果如下:1、节点导纳矩阵:-16。1905 9.5238 9.5238 -37.4084 15.3846 12.5000 15。3846 -35.3846 20。0000 12。500
13、0 20.0000 -37。9348 5.4348 5。4348 -9。9802 2、节点阻抗矩阵: 0.1181 0.0958 0.0902 0。0858 0.0467 0.0958 0.1629 0.1533 0.1459 0。0794 0。0902 0.1533 0.1860 0.1611 0.0877 0.0858 0。1459 0。1611 0。1729 0.0941 0.0467 0。0794 0.0877 0.0941 0。1515 3、短路节点是:34、短路电流:If=5.37675、各节点电压:V1=0.5152V2=0.1758V3=0V4=0。1336V5=0.52826
14、、发电机支路电流如下:Ig10=-3。2321Ig50=-2.14457、各支路电流:I12=-3。2321I23=-2。7046I24=-0.5275I43=2。6720I54=-2。1445图38 对称短路计算结果图38 改变短路点后计算结果【习题6-10】如图3-9所示网络,略去负荷,试用节点阻抗矩阵法求节点5发生三相短路时,短路点的短路电流及线路L2、L3的电流。已知各元件参数如下:发电机G1: SN = 120 MVA, x”d = 0。23; G2: SN = 60 MVA, x”d = 0。14变压器T-1: SN = 120 MVA, Vs = 10。5%; T2: SN =
15、60 MVA, Vs = 10。5%线路参数 x1 = 0。4/km, b1 = 2.810-6 S/km线路长度 L1: 120 km, L2: 80 km, L3: 70 km取SB = 120 MVA, VB = VavGGG-1G-2T-1T-2L-1L-2L-312345图3-9 电力系统网络图图3-10 三相短路时的等值网络图(1)、形成输入文件:系统参数如下所示:1。系统的节点数为:52.系统的支路数为:53。系统的短路点为:5支路电抗参数如下:z3 5=0。29036 z4 5=0.25406 z3 4=0.43554z1 3=0。1050z2 4=0.2100发电机接入点及其
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 电力系统 分析 语言
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【人****来】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【人****来】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。