基于MATLAB和VisualBasic的电力系统课程设计报告书

《基于MATLAB和VisualBasic的电力系统课程设计报告书》由会员分享，可在线阅读，更多相关《基于MATLAB和VisualBasic的电力系统课程设计报告书（26页珍藏版）》请在装配图网上搜索。

1、 . . . 信息工程学院 电力系统分析课程设计报告书 题目:电力系统短路故障的计算机算法程序设计 专 业： 班 级： 学 号： 学生： 指导 年 月 日 信息工程学院课程设计任务书学生学 号成 绩设计题目 电力系统短路故障的计算机算法程序设计设计容电力系统发生短路故障造成的危害性是最大的。作为电力系统三大计算之一,分析短路故障的参数更为重要。通过课程设计, 使学生巩固电力系统三相短路计算的基本原理与方法,掌握短路电流的数值求解方法,开发系统短路故障电流的计算程序。同时,通过软件开发,也有助于计算机操作能力和软件开发能力的提高。设计要求1、运用合理的计算方法,就是要求所选用的计算方法能快速准确

2、地得出正确结果,同时还应要求在解算过程中占用存少,以利提高计算机的解题规模。2、选择合适的语言编写程序,就是首先确定用什么计算机语言来编制程序；其次是作出计算的流程图；第三根据流程图用选择的语言编写计算程序。3、考虑影响问题的主要因素,而忽略一些次要因素,使数学模型既能正确地反映实际问题,又使计算不过于复杂。时间安排 1、先熟悉编程语言一周时间 2、复习电力系统分析相关容一周时间 3、建立模型,确定算法,编程,完成任务参考资料 1电力系统分析上册 华中科技大学 何仰赞 2电力系统故障的计算机辅助分析 大学 米麟书等 3电力系统故障分析 清华大学 周荣光 4短路电流实用计算方法 电力工业 交通大

3、学等目 录1 任务提出与方案论证11.1 任务提出11.2 方案论证12设计思路22.1 项目分析23 详细设计43.1 VB界面制作详细程序43.2 测试算例94 总 结18参 考 文 献19附 录.2021 / 261 任务提出与方案论证1.1 任务提出电力系统的故障有很多种,而短路则是最为严重的一种,包括一相短路接地,两相短路,两相短路接地,三相短路,其中三相短路在电力系统中虽然发生的概率很小,但一旦发生产生的影响非常严重,电力系统中设备的动稳定和热稳定校验都要以三相短路电流为依据,因此我选择做电力系统三相短路的计算程序更有意义。电力系统发生短路故障造成的危害性是最大的。作为电力系统三大

4、计算之一,分析短路故障的参数更为重要。通过课程设计, 使学生巩固电力系统三相短路计算的基本原理与方法,掌握短路电流的数值求解方法,开发系统短路故障电流的计算程序。同时,通过软件开发,也有助于计算机操作能力和软件开发能力的提高。1.2 方案论证本课程设计根据电力系统分析课程中所讲的基于节点方程的三相短路计算的原理和方法,编程实现了三相短路计算。本程序输入含有网络拓扑信息和电器元件的参数,输出三相短路的短路电流和短路后网络各节点的电压参数和各支路的电流参数,并以文件的形式保存,方便用户的读取、打印,实现了设计的要求。对于编程语言的选择有很多种,目前比较流行的是C语言和MATLAB两种,因为这两种语

5、言相对较简单,而MATLAB作为专门的数学软件又更具简介性和直观性,是一种基于解释的语言,目标在于用户方便,开发简单,但编译器则承担了更大重任,因此程序执行效率方面也有损失。VB语言是目前世界上最为流行使用最为广泛的面向对象设计语言用他来设计界面,很多软件的界面都使用VB编写的,能起到美化界面的作用,就我个人而言,我学过VB课程,相对于其他语言来说,VB的运用更为熟悉和了解,因此,在本课程设计中我选择VB语言进行界面的美化,MATLAB在后台调用函数对电力系统对称短路计算,从而实现VB和MATLAB联合编程。2设计思路2.1 项目分析电力系统短路电流的实际工程计算中,许多实际问题的解决并不需要

6、十分精确地结果,因此采用近似计算法,在建立系统节点方程时,采用网络模型和参数做简化处理,等值电路的制定通常将发电机作为含源支路,表示为节点注入电流源和节点电抗的形式,系统中存在同步发电机时,按发电机支路处理,同时忽略发电机、变压器和输电线路的电阻,不计线路的电容,略去变压器的励磁电流,负荷忽略不计,在标幺参数计算时,选取各级平均额定电压作为基准电压,忽略元件的额定电压和平局电压的差别,认为变压器变比的标幺值为1,假定所有的发电机的电势具有相同的相位,所有的元件仅用电抗表示,避免复数运算。本课程设计的关键部分是利用系统节点导纳矩阵形成节点阻抗矩阵,本程序采用LDU分解法,然后利用节点导纳矩阵和阻

7、抗矩阵的乘积为单位矩阵列些方程,节点导纳矩阵经LDU分解后,回代求解系统阻抗矩阵。利用节点阻抗矩阵进行短路计算,首先利用故障前的节点注入电流和系统节点阻抗矩阵求取故障前系统的各节点电压,然后利用金属短路时短路点残压为0、对地点电阻为0,计算短路点电流,最后利用故障前节点电压、节点阻抗和短路点电流求得短路后各节点电压和各支路电流。应用计算机进行电力系统计算,首先要建立电力系统相应的数学模型,如建立电力系统等效网络和节点方程,其次是运用合理的计算方法计算结果,如三角分解发,LDU分解法等,第三则是选择合适的计算机语言编制计算程序,实现程序的准确性,提高精确度。本设计利用节点方程来实现电力系统三相短

8、路计算。首先,根据给定的电力系统运行方式制定系统的一相等值网络,计算出各元器件的参数,其次根据网络方程利用变压器和线路参数形成不含发电机和负荷的节点导纳矩阵Y,然后利用系统形成的节点导纳矩阵求解系统阻抗矩阵Z,最后应用节点阻抗矩阵计算短路电流。2.2 总流程图程序的总流程图如图2-1所示。输入文件名输入节点数、支路数、故障输入电抗支路参数输入发电机支路参数输入电容支路参数结束 3 详细设计3.1 VB界面制作详细程序Dim n%Dim f%Dim z_data As DoubleDimt%Private Sub cmmd1_Clickn = Valf = ValIf n 0 And f 0 T

9、hencmmd1.Enabled = FalseCmmd3.Enabled = TrueLab1.Caption = 总共 & Str & 个节点 & Space & 当前输入 & Str & 个End IfReDim z_data 重定义大小End SubPrivate Sub cmmd2_Clickn = 0f = 0Text1.Text = Text2.Text = Text3.Text = cmmd1.Enabled = TrueCmmd3.Enabled = FalseCmmd4.Visible = FalseLab1.Caption = Pic.ClsPic.ClsEnd SubP

10、rivate Sub Cmmd3_ClickDim I%If cnt Then For I = 0 To n - 1 z_data = ValSplit NextEnd Ifcnt =t + 1Lab1.Caption = 总共 & Str & 个节点 & Space & 当前输入 & Str & 个If Thencmmd1.Enabled = Falsecmmd2.Enabled = FalseCmmd4.Visible = Truecnt = 0Text3.Text = End IfEnd SubPrivate Sub Cmmd4_ClickDim Y As VariantDim I1 A

11、s VariantDim I As VariantDim X As IntegerDim j As IntegerDim qiu As New current_c.currentDim temp As DoubleCall qiu.curren_catPic.PrintPic.Print 导纳矩阵：Dim strout As StringDim s As String For X = 1 To n For j = 1 To n temp = RoundY, 4 s = Str If Abs 1 And temp 0 Then s = StrAbs If temp 0 Then s = 0 &

12、Trim Else s = -0 & Trim End If End If Pic.Print s, Next j Pic.Print s = Next XPic.PrintPic.Print 各支路电流：Pic.Print For X = 1 To n For j = 1 To n temp = RoundI, 4 s = Str If Abs 1 And temp 0 Then s = StrAbs If temp 0 Then s = 0 & Trim Else s = -0 & Trim End If End If Pic.Print s, Next j Pic.Print s = N

13、ext X Pic.Print strout = Chr & StrRound Pic.Print 短路点电流：, strout Cmmd4.Enabled = False cmmd2.Enabled = TrueEnd SubPrivate Sub Command1_ClickDim stra As StringDim Y%Dim d%Dim temp As StringDim tdata As DoubleWith CommonDialog1 .CancelError = True .InitDir = App.Path 予设存档路径 .Filter = 文档文件|*.txt .ShowS

14、ave 或commonDialog1.Action = 2 保存文件 End WithOpen CommonDialog1.FileName For Input As #1stra = StrConvInputB$LOF, #1, vbUnicodeClose #1temp = SplitReDim tdataUBoundFor Y = 0 To UBoundtdata = ValtempNextn = tdataf = tdataText1.Text = nText2.Text = fFor Y = 2 To UBoundtdata = tdataNextReDim z_dataFor Y

15、= 0 To n - 1 For d = 0 To n - 1 z_data = tdata NextNextCmmd4.Enabled = TrueCmmd4.Visible = Truecmmd1.Enabled = FalseEnd SubPrivate Sub Form_LoadCmmd3.Enabled = FalseCmmd4.Visible = Falsecnt = 0End SubPrivate Sub Text1_KeyPress If KeyAscii Asc Or KeyAscii Asc Then KeyAscii = 0 End IfEnd SubPrivate Su

16、b Text2_KeyPress If KeyAscii 57 Then KeyAscii = 0 End IfEnd Sub3.2 测试算例例6-3在如图3-5所示的电力系统中分别在节点1和节点5接入发电机支路,其标幺值参数为：。在节点3发生三相短路,计算短路电流及网络中的电流分布。线路的电阻和电容略去不计,变压器的标幺变比等于1。各元件参数的标幺值如下：图3-5 电力系统等值网络图其三相短路时的等值网络图如图3-6所示。图3-6 三相短路时的等值网络图由以上网络图可得计算过程如下：1、形成输入文件：系统参数如下所示:1.系统的节点数为52.系统的支路数为53.系统的短路点为3支路电抗参数如

17、下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.5000 20.0000 -37.9348 5.4348 5.4348 -9.9802 2、节点阻抗矩阵: 0.1181 0.0958 0.0902

18、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、发电机支路电流如下:Ig10=-3.2321Ig50=-2.14457、各支路电流:I12=-3.2321I23=-2.7046I24=-0.

19、5275I43=-2.6720I54=-2.1445程序运行结果如图3-7所示。图3-7 程序运行结果 习题6-10如图3-8所示网络,略去负荷,试用节点阻抗矩阵法求节点5发生三相短路时,短路点的短路电流及线路L-2、L-3的电流。已知各元件参数如下：发电机G-1: SN = 120 MVA, xd = 0.23; G-2: SN = 60 MVA, xd = 0.14变压器T-1: SN = 120 MVA, Vs = 10.5%; T-2: SN = 60 MVA, Vs = 10.5%线路参数 x1 = 0.4/km, b1 = 2.810-6 S/km线路长度 L-1: 120 km,

20、 L-2: 80 km, L-3: 70 km取SB = 120 MVA, VB = VavGGG-1G-2T-1T-2L-1L-2L-312345图3-8 电力系统网络图其三相短路时的等值网络图如图3-9所示。图3-9 三相短路时的等值网络图由以上网络图可得计算过程如下：1、形成输入文件：系统参数如下所示:1.系统的节点数为：52.系统的支路数为：53.系统的短路点为：5支路电抗参数如下：z3 5=0.29036z4 5=0.25406z3 4=0.43554z1 3=0.1050z2 4=0.2100发电机接入点及其电抗如下：共2台发电机进入系统：z1=0.2300z2=0.2800e1=

21、1.0500e2=1.0500共3个节点有接地电容：Y3=-32.372299Y4=-34.077000Y5=-43.1852002、得到的输出文件：计算结果如下：1、节点导纳矩阵:-13.8716 9.5238 -8.3333 4.7619 9.5238 -15.2329 2.2960 3.4440 4.7619 2.2960 -10.9646 3.9361 3.4440 3.9361 -7.3569 2、节点阻抗矩阵: 0.1819 0.0616 0.1600 0.1079 0.1326 0.0616 0.2078 0.0898 0.1536 0.1242 0.1600 0.0898 0.

22、2330 0.1571 0.1932 0.1079 0.1536 0.1571 0.2689 0.2174 0.1326 0.1242 0.1932 0.2174 0.3427 3、短路节点是:54、短路电流：If=3.12635、各节点电压:V1=0.6471V2=0.6723V3=0.4632V4=0.3890V5=06、发电机支路电流如下:Ig10=-1.7516Ig20=-1.34917、各支路电流:I35=-1.5954I45=-1.5310I34=-0.1705I13=-1.7516I24=-1.3491程序运行结果如图3-10所示。图3-10 程序运行结果习题6-11电力系统等值

23、电路如图3-11所示,支路阻抗的标幺值已注明图中。1形成节点导纳矩阵或节点阻抗矩阵,并用以计算节点3的三相短路电流。2另选一种方法计算短路电流,并用以验证1的计算结果。图3-11 电力系统等值电路图由以上网络图可得计算过程如下：1、形成输入文件：系统参数如下所示:1.系统的节点数为：32.系统的支路数为：23.系统的短路点为：3支路电抗参数如下：z1 2=0.2000z2 3=0.1000发电机接入点及其电抗如下：共2台发电机进入系统：z1=0.2000z3=0.130435e1=1.0500e3=0.913043共1个节点有接地电容：Y2=-2.0000002、得到的输出文件：计算结果如下：

24、1、节点导纳矩阵:-10.0000 5.0000 5.0000 -14.5000 10.0000 10.0000 -17.6667 2、节点阻抗矩阵: 0.1394 0.0789 0.0446 0.0789 0.1577 0.0893 0.0446 0.0893 0.1071 3、短路节点是:34、短路电流：If=9.18755、各节点电压:V1=0.6344V2=0.2188V3=06、发电机支路电流如下:Ig10=-2.0781Ig30=-7.00007、各支路电流:I12=-2.0781I23=-2.1875运行结果如图3-12所示。图3-12 程序运行结果程序运行结果与算例标准答案一样

25、,说明本程序的计算方法和思路的准确性,程序基本符合设计的要求。4 总 结本次课程设计使我感受最深的是计算机在电力系统领域里的巨大应用价值。在这次设计中,用手工计算是相当麻烦的,短路电流所涉及的参数众多,计算公式也是相当繁冗,手工计算起来十分地乏味,且需要有相当的耐心和精力,稍有不慎,即会出错。而通过计算机来进行计算,我们只需要把计算的流程搞懂了,设计出相关的程序,简简单单地敲击键盘就可以得出答案来了,且答案与实际值相差无几。我想,计算机不仅仅在这一次设计中体现出其价值,在今后的学习里,我们还一定会用到计算机编程来解决问题的情况。因此,我觉得学好计算机,掌握好一门计算机语言是非常有必要的。电力系

26、统短路电流故障是电力系统三大计算之一,学好它对于学好电力系统这门课程,对以后的毕业设计和实际工作具有十分重要的意义。而利用计算机来计算和分析其大小、流向是十分方便的。因此,在今后的学习中,我会继续努力学好电力系统分析这门课程,并学习另外一些相关的课外知识,开阔自己的视野。该课程设计给出了基于组合故障模型的电力系统不对称短路故障通用算法公式的推导过程a一方面解决了任意过渡阻抗的故障计算问题a另一方面拓宽了故障的计算围a计及了组合型故障发生的可能a因此涵盖了短路故障所有特殊情形O 该算法针对不对称短路故障而言具有严格意义上的通用性a大大简化了计算机的编程实现a有效地实现了任意基本故障型短路及组合故

27、障型短路的统一计算a具有应用价值。参 考 文 献1何仰赞. 电力系统分析上册. ：华中科技大学,20022米麟书等. 电力系统故障的计算机辅助分析. ：大学,19923周作仁等. 电网计算与程序设计. ：科学技术,19834Visual Basic 程序设计教程. ：人民邮电,2007.25康华光,大钦. 电子技术基础模拟部分第五版M.：高等教育,20056皮文兵. 一种宽输入围的Gillbert模拟乘法器设计J. 电子设计应用.2007.137 焦春生.新型绿色能效D类音频放大器设计应用EB/OL8 筑生. 微分半动力系统的不变集 D. : 大学数学系数学研究所, 19839 MATLAB

28、7.0从入门到精通修订版. ：人民邮电 ,2010.5附 录主要数据变量说明如下：n ：总节点数f ：短路点z ：系统阻抗矩阵 ：节点i和节点j之间的阻抗,如无连接,则0Z ：系统的阻抗Z矩阵Y ：系统的导纳矩阵If ：短路电流 ：短路后节点i电压的向量 ：短路后节点i和j之间的支路电流计算主程序如下：load impedance.dat; 载入阻抗矩阵load point.dat; 载入短路节点,节点总数 n=point; 赋节点总数值f=point; 赋短路节点值z=1j.*impedance;Y=zeros;I=zeros; 初始定义Y,I阵为零矩阵for i=1:n for j=1:n

29、 if z=0 Y=Y+1/z; 求出Y阵中的对角线上元素,自导纳 end end endfor i=1:n for j=i+1:n if z=0 Y=0; else Y=-1/z; 求出Y阵中非对角线上元素,互导纳 end Y=Y; endendZ=inv; 导纳矩阵求逆,得出阻抗矩阵If=1/Z; 计算短路电流,公式来源6-10V=ones; 产生长度为n的全1行数组,定义各节点电压为额定制for i=1:n 求解短路后各节点电压,公式来源6-8 V=V-If*Z;endfor i=1:n for j=1:n I=V-V*Y; 求出各支路电流,公式来源6-9 endendY=imag; 取

30、出导纳矩阵的虚部为导纳矩阵Z=imag;If=imag;I=imag;save daona.dat Y -ascii 保存Y阵结果存于daona文件中save zukang.dat Z -ascii 保存Z阵结果存于zukang文件中save current.dat I -ascii 保存支路电流存于currrent文件中save voltage.dat V -ascii 保存各节点电压存于voltage文件中save d-currentt.dat If -ascii 保存短路电流信息存于d-current文件中save zong.m Y Z I V If -ascii 保存所有信息于总文件中