




已阅读5页,还剩18页未读, 继续免费阅读
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
课程 电力系统暂态上机课程设计报告课程设计报告( 20122013年度第二学期)名 称: 院 系: 班 级: 学 号: 学生姓名: 指导教师: 设计周数: 两 周 成 绩: 日期:2013年 7月 10 日21一、课程设计的目的与要求巩固电力系统暂态分析理论知识,使学生掌握采用计算机对电力系统电磁暂态过程和机电暂态过程进行计算的基本方法,并进一步巩固计算机编程能力,为将来从事相关的技术工作打下必要的基础。二、设计正文(详细内容见附录,用A4纸,页数不限)1 对称短路计算过程流程图和计算结果2 不对称短路计算过程流程图和计算结果3 静态稳定性计算过程流程图和计算结果4 暂态稳定性计算过程流程图和计算结果5 思考题三、课程设计总结或结论本次电力系统暂态上机主要讨论了计算各种类型短路故障下系统网络中的电压电流分布以及电力系统的静态稳定、暂态稳定问题。通过本次课程设计,本人对电力系统故障分析有了更深刻的理解,电力系统的故障时,大部分电磁量将随时间变化,描述其特性的是微分方程,这给分析计算带来一定困难。在分析过程中通常尽量避免对微分方程直接求解,而是采用一定的工具和假设使问题得以简化,即把“微分方程代数化,暂态分析稳态化”。在分析不对称故障时,各相之间电磁量的耦合使问题的分析更为复杂,此时常用的分析方法是采用对称分量法将不对称问题转化为对称问题来求解。同时我对用来分析电力系统静态稳定的试探法,用来分析电力系统暂态稳定的改进欧拉法有了一些使用心得。这与手算系统短路电流时使用的网络化简方法大大不同。另外,通过此次课程设计,我学习并掌握了一定的Matlab知识,对简单的选择语句、循环语句、数组的使用等有了一些了解,为今后学习、工作中再次学习使用此软件打下了基础,获益匪浅。四、参考文献1. 电力系统暂态分析,李光琦,中国电力出版社,2006年,第三版;2. 电力系统分析(上、下),何仰赞,华中科技大学出版社,1998年,第二版。3. 电力系统故障的计算机辅助分析 重庆大学出版社 米麟书等 4. 电力系统潮流计算 天津大学出版社 宋文南等 5. 电力系统故障分析 清华大学出版社 周荣光6. 短路电流实用计算方法 电力工业出版社 西安交通大学等7. 精通Matlab6.5北京航空航天大学出版社,张志涌附录(设计流程图、计算结果、思考题答案)1 对称短路计算过程流程图和计算结果流程图输入数据选择故障类型用节点导纳阵求逆得到节点阻抗矩阵计算各节点的序电压、相电压根据故障类型计算故障点序电流、相电流计算各支路的序电流、相电流计算发电机节点的序电压、相电压形成节点导纳矩阵选择短路点开始计算发电机支路的序电流、相电流结束计算结果:(1) 导纳矩阵:Y1 =Y2 = Y0 = YY1=YY2 = (2) 节点3发生三相短路故障a. 故障点三序电流、三相电流:I1= 0 - 9.8592i I2=0 I3=0 ; Ia= 0 - 9.8592i Ib= -8.5383 + 4.9296i Ic= 8.5383 + 4.9296ib.各节点三序电压、三相电压:节点U1U2U0UaUbUc10.45070.00000.00000.4507-0.2254 - 0.3903i -0.2254 + 0.3903i20.53520.00000.00000.5352-0.2676 - 0.4635i-0.2676 + 0.4635i30.00000.00000.00000.00000.00000.000040.63380.00000.00000.5489 + 0.3169i0.0000 - 0.6338i-0.5489 + 0.3169i50.69010.00000.00000.5977 + 0.3451i0.0000 - 0.6901i-0.5977 + 0.3451ic.各支路三序电流、三相电流:支路I1I2I0IaIbIc1-20 + 0.8451i0.00000.00000 + 0.8451i0.7319 - 0.4225i-0.7319 - 0.4225i1-30 - 4.5070i0.00000.00000 - 4.5070i-3.9032 + 2.2535i3.9032 + 2.2535i2-30 - 5.3521i0.00000.00000 - 5.3521i-4.6351 + 2.6761i4.6351 + 2.6761i1-40 + 3.6620i0.00000.00000 + 3.6620i3.1714 - 1.8310i-3.1714 - 1.8310i2-50 + 6.1972i0.00000.00000 + 6.1972i5.3669 - 3.0986i-5.3669 - 3.0986i发电机1-1.8310 + 3.1714i00-1.8310 + 3.1714i3.6620 + 0.0000i-1.8310 - 3.1714i发电机23-3.0986 + 5.3669i00-3.0986 + 5.3669i6.1972 + 0.0000i-3.0986 - 5.3669i2不对称短路计算过程流程图和计算结果计算结果:(1) 节点3发生A相短路接地故障a. 故障点三序电流、三相电流:I1= 0 - 3.1080i I2= 0 - 3.1080i I3= 0 0 - 3.1080i ; Ia= 0 - 9.3241i Ib= 0 - 0.0000i Ic= 0 - 0.0000ib. 各节点三序电压、三相电压:节点U1U2U0UaUbUc10.8268-0.1732-0.07150.5822-0.3983 - 0.8660i -0.3983 + 0.8660i20.8535-0.1465-0.0459 0.6610-0.3994 - 0.8660i-0.3994 + 0.8660i30.6848-0.3152-0.36950.0000-0.5543 - 0.8660i-0.5543 + 0.8660i40.8846-0.11540.00000.6661 + 0.5000i0.0000 - 1.0000i-0.6661 + 0.5000i50.9023-0.09770.00000.6968 + 0.5000i0.0000 - 1.0000i-0.6968 + 0.5000ic. 各支路三序电流、三相电流:支路I1I2I0IaIbIc1-20 + 0.2664i0 + 0.2664i0 + 0.1280i0 + 0.6608i0.0000 - 0.1384i-0.0000 - 0.1384i1-30 - 1.4208i0 - 1.4208i0 - 1.4900i0 - 4.3317i0.0000 - 0.0692i-0.0000 - 0.0692i2-30 - 1.6872i0 - 1.6872i0 - 1.6180i0 - 4.9924i0 + 0.0692i0 + 0.0692i1-40 + 1.1544i0 + 1.1544i0 + 1.3621i0 + 3.6709i0.0000 + 0.2077i-0.0000 + 0.2077i2-50 + 1.9536i0 + 1.9536i0 + 1.7459i0 + 5.6532i0.0000 - 0.2077i-0.0000 - 0.2077i发电机1-0.5772 + 0.9997i0.5772 + 0.9997i0-0.0000 + 1.9995i0.0000 + 0.0000i-0.0000 - 1.9995i发电机23-0.9768 + 1.6919i0.9768 + 1.6919i0 -0.0000 + 3.3838i0.0000 + 0.0000i-0.0000 - 3.3838i(2) 节点3发生A相经10电阻接地故障a. 故障点三序电流、三相电流:I1= 0.0333 - 0.0004i I2= 0.0333 - 0.0004i I3= 0.0333 - 0.0004i ; Ia= 0.1000 - 0.0011i Ib= 0 Ic= 0b. 各节点三序电压、三相电压:节点U1U2U0UaUbUc10.8531 - 0.0621i-0.1469 - 0.0621i-0.0607 - 0.0257i0.6456 - 0.1499i-0.4138 - 0.8296i-0.4138 + 0.9025i20.8757 - 0.0526i-0.1243 - 0.0526i-0.0389 - 0.0165i0.7125 - 0.1216i-0.4147 - 0.8299i-0.4147 + 0.9021i30.7326 - 0.1131i-0.2674 - 0.1131i-0.3134 - 0.1326i 0.1518 - 0.3588i-0.5460 - 0.8855i-0.5460 + 0.8466i40.9021 - 0.0414i-0.0979 - 0.0414i00.6964 + 0.4283i0.0000 - 1.0000i-0.6964 + 0.5717i50.9171 - 0.0350i-0.0829 - 0.0350i00.7225 + 0.4393i0.0000 - 1.0000i-0.7225 + 0.5607ic. 各支路三序电流、三相电流:支路I1I2I0IaIbIc1-2-0.0956 + 0.2260i-0.0956 + 0.2260i-0.0459 + 0.1085i-0.2371 + 0.5605i0.0497 - 0.1174i0.0497 - 0.1174i1-30.5098 - 1.2052i0.5098 - 1.2052i0.5346 - 1.2639i1.5542 - 3.6742i0.0248 - 0.0587i0.0248 - 0.0587i2-30.6054 - 1.4311i0.6054 - 1.4311i0.5806 - 1.3724i1.7913 - 4.2346i-0.0248 + 0.0587i-0.0248 + 0.0587i1-4-0.4142 + 0.9792i-0.4142 + 0.9792i-0.4887 + 1.1553i-1.3172 + 3.1137i-0.0745 + 0.1761i-0.0745 + 0.1761i2-5-0.7010 + 1.6571i-0.7010 + 1.6571i-0.6265 + 1.4809i-2.0284 + 4.7951i0.0745 - 0.1761i0.0745 - 0.1761i发电机1-0.8483 + 0.6409i0.1309 + 1.0551i0-0.7174 + 1.6960i0.0000 - 0.0000i0.7174 - 1.6960i发电机23-1.4356 + 1.0846i0.2215 + 1.7856i0 -1.2141 + 2.8702i-0.0000 - 0.0000i1.2141 - 2.8702i(3) 节点3发生b、c两相短路故障a. 故障点三序电流、三相电流:I1= 0 - 4.9296i I2= 0 + 4.9296i I3= 0 ; Ia= 0 Ib= -8.5383 - 0.0000i Ic= 8.5383 + 0.0000ib.各节点三序电压、三相电压:节点U1U2U0UaUbUc10.72540.27460.00001.0000-0.5000 - 0.3903i-0.5000 + 0.3903i20.76760.23240.00001.0000-0.5000 - 0.4635i-0.5000 + 0.4635i30.50000.50000.00001.0000-0.5000 - 0.0000i-0.5000 + 0.0000i40.81690.18310.00000.8660 + 0.3169i0.0000 - 0.6338i-0.8660 + 0.3169i50.84510.15490.00000.8660 + 0.3451i0.0000 - 0.6901i-0.8660 + 0.3451ic.各支路三序电流、三相电流:支路I1I2I0IaIbIc1-20 + 0.4225i0 - 0.4225i0.00000 + 0.0000i0.7319 - 0.0000i-0.7319 - 0.0000i1-30 - 2.2535i0 + 2.2535i0.00000 - 0.0000i-3.90323.9032 + 0.0000i2-30 - 2.6761i0 + 2.6761i0.00000 - 0.0000i-4.6351 + 0.0000i4.6351 + 0.0000i1-40 + 1.8310i0 - 1.8310i0.00000 + 0.0000i3.1714 - 0.0000i-3.1714 - 0.0000i2-50 + 3.0986i0 - 3.0986i0.000005.3669 + 0.0000i-5.3669 - 0.0000i发电机1-0.9155 + 1.5857i-0.9155 - 1.5857i0-1.8310 + 0.0000i3.6620 + 0.0000i-1.8310 - 0.0000i发电机23-1.5493 + 2.6835i-1.5493 - 2.6835i0-3.0986 6.1972 + 0.0000i-3.0986 - 0.0000i(4) 节点3发生b、c两相短路接地故障a. 故障点三序电流、三相电流:I1= 0 - 6.4036i I2=0 + 3.4556i I3= 0 + 2.9480i ; Ia= 0 Ib= -8.5383 + 4.4220i Ic= 8.5383 + 4.4220ib. 各节点三序电压、三相电压:节点U1U2U0UaUbUc10.64320.19250.06780.9036-0.3500 - 0.3903i-0.3500 + 0.3903i20.69810.16290.04360.9046-0.3870 - 0.4635i-0.3870 + 0.4635i30.35050.35050.35051.05150040.76220.128300.7712 + 0.3169i0.0000 - 0.6338i-0.7712 + 0.3169i50.79870.108600.7858 + 0.3451i0.0000 - 0.6901i-0.7858 + 0.3451ic. 各支路三序电流、三相电流:支路I1I2I0IaIbIc1-20 + 0.5489i0 - 0.2962i0 - 0.1214i 0 + 0.1313i0.7319 - 0.2477i-0.7319 - 0.2477i1-30 - 2.9274i0 + 1.5797i0 + 1.4133i0 + 0.0657i-3.9032 + 2.0872i3.9032 + 2.0872i2-30 - 3.4762i0 + 1.8759i0 + 1.5347i0 - 0.0657i-4.6351 + 2.3349i4.6351 + 2.3349i1-40 + 2.3785i0 - 1.2835i0 - 1.2920i0 - 0.1970i3.1714 - 1.8394i-3.1714 - 1.8394i2-50 + 4.0251i0 - 2.1721i0 - 1.6561i0 + 0.1970i5.3669 - 2.5826i-5.3669 - 2.5826i发电机1-1.1892 + 2.0598i-0.6417 - 1.1115i0-1.8310 + 0.9483i3.6620 + 0.0000i-1.8310 - 0.9483i发电机23-2.0126 + 3.4858i-1.0860 - 1.8811i0-3.0986 + 1.6048i6.1972 + 0.0000i-3.0986 - 1.6048i2 静态稳定性计算过程流程图和计算结果流程图 否否否是是是开始清屏设初值:Xd_all,Xdd_all,Xq,Xq_all,Eq0,UG0,Xe,Xd,XddKe=0.1Ke_result=Ke,delta_result=delta,P_result=Eq*U*sin(delta*pi/180)/Xd_all, delta=delta+1delta=49求Kemax,Kemin求Eq,Eqq求UGd,UGq,UG求K1K6Ke5.7delta=110(V10)&(V20)&(V30)结束打印Ke_result,delta_result, P_resultKe=ke+0.1计算结果:Ke=0.100000, delta=91, P=1.345117Ke=0.200000, delta=92, P=1.363000Ke=0.300000, delta=93, P=1.380693Ke=0.400000, delta=95, P=1.398409Ke=0.500000, delta=96, P=1.415747Ke=0.600000, delta=97, P=1.432842Ke=0.700000, delta=98, P=1.449679Ke=0.800000, delta=99, P=1.466244Ke=0.900000, delta=100, P=1.482520Ke=1.000000, delta=101, P=1.498487Ke=1.100000, delta=101, P=1.513985Ke=1.200000, delta=102, P=1.529380Ke=1.300000, delta=103, P=1.544429Ke=1.400000, delta=103, P=1.558878Ke=1.500000, delta=104, P=1.573336Ke=1.600000, delta=105, P=1.587402Ke=1.700000, delta=105, P=1.600902Ke=1.800000, delta=105, P=1.613924Ke=1.900000, delta=102, P=1.621480Ke=2.000000, delta=99, P=1.622734Ke=2.100000, delta=97, P=1.623379Ke=2.200000, delta=94, P=1.614349Ke=2.300000, delta=92, P=1.607801Ke=2.400000, delta=89, P=1.589402Ke=2.500000, delta=87, P=1.576380Ke=2.600000, delta=85, P=1.560855Ke=2.700000, delta=83, P=1.542968Ke=2.800000, delta=81, P=1.522859Ke=2.900000, delta=80, P=1.514020Ke=3.000000, delta=78, P=1.490598Ke=3.100000, delta=76, P=1.465314Ke=3.200000, delta=75, P=1.453440Ke=3.300000, delta=73, P=1.425420Ke=3.400000, delta=72, P=1.411962Ke=3.500000, delta=70, P=1.381541Ke=3.600000, delta=69, P=1.366704Ke=3.700000, delta=68, P=1.351397Ke=3.800000, delta=66, P=1.318171Ke=3.900000, delta=65, P=1.301745Ke=4.000000, delta=64, P=1.284942Ke=4.100000, delta=63, P=1.267785Ke=4.200000, delta=62, P=1.250299Ke=4.300000, delta=61, P=1.232503Ke=4.400000, delta=60, P=1.214420Ke=4.500000, delta=59, P=1.196069Ke=4.600000, delta=58, P=1.177469Ke=4.700000, delta=57, P=1.158638Ke=4.800000, delta=56, P=1.139594Ke=4.900000, delta=55, P=1.120353Ke=5.000000, delta=54, P=1.100930Ke=5.100000, delta=53, P=1.081340Ke=5.200000, delta=52, P=1.061599Ke=5.300000, delta=51, P=1.041718Ke=5.400000, delta=51, P=1.041792Ke=5.500000, delta=50, P=1.021746Ke=5.600000, delta=49, P=1.001589Ke=5.600000, delta=49, P=1.001589最终选择放大倍数Ke=2.14暂态稳定性计算过程流程图和计算结果 流程图.暂态稳定性是否开始结束输入系统数据,赋初值(delta,omega,步长h,时段长度)Duration等)等输入故障切除时间CutTimei=0i=round(CutTime/h)求时段初功角、角速度变化率求时段末功角、角速度估计值求时段末功角、角速度变化率求平均功角、角速度变化率求时段末功角、角速度计算值i=i+1显示此时功角同理求故障切除后的功角变化过程画出摇摆曲线结果: 0.15s时切除故障的摇摆曲线 0.25s时切除故障的摇摆曲线根据摇摆曲线判断:0.15s时切除故障系统暂态稳定,0.25s时切除故障系统失稳。另据试探,为保证系统暂态稳定,切除时间最大值约为0.2s。思考题1计算短路电流,书中给出的手算方法与计算机编程方法有何区别?答:手算短路电流时,要将用各序等值电抗表示的复杂网络简化为一个电抗值,并根据边界条件计算出三序网络的关系做出复合序网图,然后计算出正负零各序电流,之后转化为abc三相短路电流。而用计算机时可以不用简化网络,求出各序导纳矩阵后利用I=YU来计算各节点正负零各序电流后再转换成三相短路电流。2用计算机方法,进行电力系统潮流计算和短路电流计算,有哪些区别?答:潮流计算时所求的导纳矩阵为正序的导纳矩阵,而短路计算时要求出正负零三序的导纳矩阵。潮流计算时需要迭代计算出结果,而短路计算一次就可以计算出短路电流。潮流计算需要选定PQ和PV 节点,短路计算不需要区分节点类型。3如果交给你一个任务,请你用已有的短路电流计算软件计算某地区电网某点发生短路后的短路电流,你应该做哪些工作?(收集哪些数据,如何整理,计算结果如何分析)答:1.网络中各个发电机的容量、额定电压、功率因数和正负序阻抗;变压器的容量、额定电压、额定变比和短路电压百分数;各个线路的电抗。2.计算各元件等值电抗,求出系统的等值电路;3.短路前运行状况的分析计算;4.修改程序中的参数,计算故障分量。4静态稳定的概念,暂态稳定的概念?答:静态稳定:指电力系统受到小干扰后,不发生非周期性失步或自发振荡,自动恢复到初始运行状态的能力。暂态稳定:指电力系统受到大干扰后,各发电机组保持同步运行并过渡到新的或恢复到原来稳定运行状态的能力。5提高静态稳定、暂态稳定的措施?答:提高静态稳定的措施:一、 采用自动调节励磁装置二、 减小元件的电抗1、 采用分裂导线2、 采用串联电容补偿三、 提高线路标称等级四、 改善系统的结构和采用中间补偿设备提高暂态稳定的措施:一、改变制动功率1、故障快速切除和应用自动重合闸2、对发电机进行强行励磁3、电气制动4、变压器中性点经小电阻接地5、输电线路设置开关站6、输电线路采用强行串联电容补偿二、改变原动功率1、快速的自动调速系统或快速关闭汽门2、联锁切除部分发电机3、合理选择远距离输电系统的运行接线三、系统失去稳定后的措施1、设置解列点2、短期异步运行和再同步3、做好系统“黑启动”方案程序1clc;clear;%-输入已知条件-%bus_Num1=3; %节点数bus_Num2=5; %包括发电机节点的节点数branch_Num1=3; %线路数branch_Num2=5; %包括发电机支路的支路数branch1_No1=1,1,2;branch1_No2=2,3,3;branch2_No1=1,1,2,1,2;branch2_No2=2,3,3,4,5;%输入支路各序阻抗,z1_branch表示支路正序阻抗,z2_branch表示支路负序阻抗,z0_branch表示支路零序阻抗z1_branch(1)=j*0.1;z1_branch(2)=j*0.1;z1_branch(3)=j*0.1;z2_branch=z1_branch;z0_branch(1)=j*0.2;z0_branch(2)=j*0.2;z0_branch(3)=j*0.2;%-第一步:不考虑发电机节点计算节点导纳矩阵-%节点导纳矩阵,Y1表示不计发电机节点的正序网络节点导纳阵,Y2表示不计发电机节点的负序网络节点导纳阵,Y0表示不计发电机节点的零序网络节点导纳阵,Y1=zeros(bus_Num1);Y1(1,1)=1/(j*0.15);Y1(2,2)=1/(j*0.075);%请同学们求正序节点导纳矩阵for m=1:branch_Num1 No1=branch1_No1(m); No2=branch1_No2(m); Y1(No1,No1)=Y1(No1,No1)+1/z1_branch(m); Y1(No2,No2)=Y1(No2,No2)+1/z1_branch(m); Y1(No1,No2)=-1/z1_branch(m); Y1(No2,No1)=-1/z1_branch(m);endstr=正序阻抗Y11=Y1Y2=Y1;%负序等于正序str=负序阻抗 Y22=Y2Y0=zeros(bus_Num1);Y0(1,1)=1/(j*0.0525);Y0(2,2)=1/(j*0.0263);%请同学们求零序节点导纳矩阵for m=1:branch_Num1 No1=branch1_No1(m); No2=branch1_No2(m); Y0(No1,No1)=Y0(No1,No1)+1/z0_branch(m); Y0(No2,No2)=Y0(No2,No2)+1/z0_branch(m); Y0(No1,No2)=-1/z0_branch(m); Y0(No2,No1)=-1/z0_branch(m);endstr=零序阻抗 Y00=Y0%-第二步:考虑发电机节点计算节点导纳矩阵-%节点导纳矩阵,YY1表示计及发电机节点的正序网络节点导纳阵,YY2表示计及发电机节点的负序网络节点导纳阵z1_branch(4)=j*0.05;z1_branch(5)=j*0.025;z2_branch=z1_branch;YY1=zeros(bus_Num1);YY1(4,4)=1/(j*0.1);YY1(5,5)=1/(j*0.05);%请同学们求正序节点导纳矩阵for m=1:bus_Num2 No1=branch2_No1(m); No2=branch2_No2(m); YY1(No1,No1)=YY1(No1,No1)+1/z1_branch(m); YY1(No2,No2)=YY1(No2,No2)+1/z1_branch(m); YY1(No1,No2)=-1/z1_branch(m); YY1(No2,No1)=-1/z1_branch(m);endYY2=YY1;%负序等于正序str=正负序阻抗 YY12=YY1%-第三步:计算节点阻抗矩阵-%Z1=inv(Y1);Z2=inv(Y2);Z0=inv(Y0);程序2clc;clear;%数据来源于教材电力系统暂态分析P77例(3-4),P143例(5-7)Matrix %计算节点导纳矩阵、节点阻抗矩阵,形成全局变量Fault_Node=input(输入短路点编号 ;n Fault_Node=);Fault_Type=input(输入短路类型 ;n(1)Fault_Type=0为三相短路;n(2)Fault_Type=1为a相接地短路;n(3)Fault_Type=2为a相经10欧姆电阻接地短路n(4)Fault_Type=3为bc两相相间短路n(5)Fault_Type=4为bc两相短路接地nFault_Type=)a=-0.5+j*sqrt(3)/2;T=1 1 1 % T为对称分量法的变换矩阵,见P87公式(4-4) a2 a 1 a a2 1;%-第一步:计算短路点的序电流,相电流-%计根据故障类型选择不同的计算公式,计算故障点各序电流if Fault_Type=0 I_Fault1=1/Z1(Fault_Node,Fault_Node); %三相短路只有正序电流 I_Fault2=0; I_Fault0=0;elseif Fault_Type=1 I_Fault1=1/(Z1(Fault_Node,Fault_Node)+Z2(Fault_Node,Fault_Node)+Z0(Fault_Node,Fault_Node); %3个序网阻抗串联 I_Fault2=I_Fault1; I_Fault0=I_Fault1; elseif Fault_Type=2 Zf=10*60/(1152); I_Fault1=1/(Z1(Fault_Node,Fault_Node)+Z2(Fault_Node,Fault_Node)+Z0(Fault_Node,Fault_Node)+3*Zf); %3个序阻抗和电阻串联 I_Fault2=I_Fault1; I_Fault0=I_Fault1;elseif Fault_Type=3 I_Fault1=1/(Z1(Fault_Node,Fault_Node)+Z2(Fault_Node,Fault_Node); %正序阻抗和负序阻抗并联 I_Fault2=-I_Fault1; I_Fault0=0;elseif Fault_Type=4 %三个序阻抗并联 I_Fault1=1/(Z1(Fault_Node,Fault_Node)+(Z2(Fault_Node,Fault_Node)*Z0(Fault_Node,Fault_Node)/(Z2(Fault_Node,Fault_Node)+Z0(Fault_Node,Fault_Node); I_Fault2=-I_Fault1*Z0(Fault_Node,Fault_Node)/(Z2(Fault_Node,Fault_Node)+Z0(Fault_Node,Fault_Node); I_Fault0=-I_Fault1*Z2(Fault_Node,Fault_Node)/(Z2(Fault_Node,Fault_Node)+Z0(Fault_Node,Fault_Node);end;str=短路点序电流;strI1=I_Fault1I2=I_Fault2I0=I_Fault0str=短路点相电流;strIabc=T*I_Fault1 I_Fault2 I_Fault0. %相量Iabc_effective=abs(Iabc) %有效值%-第二步:计算各个节点的序电压,相电压-%故障分量 Ifault_node1=zeros(bus_Num1,1);Ifault_node2=zeros(bus_Num1,1);Ifault_node0=zeros(bus_Num1,1);for m=1:bus_Num1 if m=Fault_Node Ifault_node1(m,1)=-1*I_Fault1; Ifault_node2(m,1)=-1*I_Fault2; Ifault_node0(m,1)=-1*I_Fault0; else Ifault_node1(m,1)=0; Ifault_node2(m,1)=0; Ifault_node0(m,1)=0; endendUfault_node1=zeros(bus_Num1,1);Ufault_node2=zeros(bus_Num1,1);Ufault_node0=zeros(bus_Num1,1);Ufault_node1=Y1Ifault_node1;Ufault_node2=Y2Ifault_node2;Ufault_node0=Y0Ifault_node0;%正常分量 Unormal_node1=ones(bus_Num1,1);Unormal_node2=zeros(bus_Num1,1);Unormal_node0=zeros(bus_Num1,1);U_node1=zeros(bus_Num1,1);U_node2=zeros(bus_Num1,1);U_node0=zeros(bus_Num1,1);U_node1=Ufault_node1+Unormal_node1; %正常分量加故障分量U_node2=Ufault_node2+Unormal_node2;U_node0=Ufault_node0+Unormal_node0;str=1.2.3节点电压U1=U_node1 U2=U_node2U0=U_node0Uabc123=T*U_node1 U_node2 U_node0.%-第三步:计算各个支路的序电流,相电流-%I_branch1=zeros(bus_Num1,1);I_branch2=zeros(bus_Num1,1);I_branch0=zeros(bus_Num1,1);for i=1:bus_Num1 m=branch1_No1(i); n=branch1_No2(i); I_branch1(i,1)=(U_node1(m,1)-U_node1(n,1)/z1_branch(i); I_branch2(i,1)=(U_node2(m,1)-U_node2(n,1)/z2_branch(i); I_branch0(i,1)=(U_node0(m,1)-U_node0(n,1)/z0_branch(i);endstr=各支路序电流;strI1= I_branch1I2= I_branch2I0= I_branch0str=各支路相电流;Iabc1=T*I_branch1 I_branch2 I_branch0.Iabc1_effective=abs(Iabc1)%-第四步:计算发电机节点的序电压,相电压-%I2fault_node1=zeros(bus_Num2,1);I2fault_node2=zeros(bus_Num2,1);I2fault_node0=zeros(bus_Num2,1);for m=1:bus_Num2 if m=Fault_Node I2fault_node1(m,1)=-1*I_Fault1; I2fault_node2(m,1)=-1*I_Fault2; I2fault_node0(m,1)=-1*I_Fault0; else I2fault_node1(m,1)=0; I2fault_node2(m,1)=0; I2fault_node0(m,1)=0; endendU2fault_node1=zeros(bus_Num2,1);U2fault_node2=zeros(bus_Num2,1);U2fault_n
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 二零二五年度高品质实木地板专业施工合同
- 二零二五版智能社区房地产开发经营合同规范
- 二零二五年度矿业权抵押担保合作协议范本
- 二零二五年度茶叶电商平台技术开发合同
- 二零二五年度房产项目股权转让与代建服务一体化协议
- 2025版新能源工程师岗位责任及绩效合同样本
- 二零二五年度餐饮企业员工绩效考核合同
- 2025防火门节能环保型安装工程合同
- 2025年海水淡化及水处理设备项目合作计划书
- 稚子弄冰说课课件
- 云端漫步云端飞车创新创业项目商业计划书
- 2025年中国工程质量检测行业市场前景预测及投资价值评估分析报告
- 宁夏资环技术有限公司招聘考试真题2024
- 高职院校与企业合作中的资源整合与共享
- 2025至2030中国烫金箔行业发展趋势分析与未来投资战略咨询研究报告
- 2025云南省初中学业水平考试数学
- 税务稽查程序培训
- 2025年 长春中医药大学附属医院招聘考试笔试试题附答案
- 地理●甘肃卷丨2024年甘肃省普通高中学业水平等级性考试高考地理真题试卷及答案
- 全国公开课一等奖七年级历史统编版上册《第4课夏商西周王朝的更替》课件(内嵌视频)
- T/SHEPEA 003-202435 kV架空导线绝缘护套技术规范
评论
0/150
提交评论