电力系统暂态上机计算课程设计报告.doc

课程设计报告( 20102011年度第二学期)名 称：电力系统暂态上机 院 系：电气与电子工程学院 班 级： 学 号： 学生姓名： 指导教师： 齐郑 设计周数： 两周 成 绩： 日期： 2011年 7月8日 课程 电力系统暂态上机课程设计报告一、课程设计的目的与要求巩固电力系统暂态分析理论知识，使学生掌握采用计算机对电力系统电磁暂态过程和机电暂态过程进行计算的基本方法，并进一步巩固计算机编程能力，为将来从事相关的技术工作打下必要的基础。二、设计正文1 短路计算过程流程图根据系统参数计算节点阻抗矩阵故障类型判断根据故障类型所对应的复合序网计算故障点各序电流分量由故障点各序短路电流计算电流故障分量，并计算各节点电压故障分量故障后各节点电压计算各支路电流计算发电机机端电压计算流过发电机支路的电流短路计算结果：三相短路ABC绝对值节点1 Uabc0.4507-0.2254 - 0.3903i-0.2254 + 0.3903i0.4507节点2 Uabc0.5352-0.2676 - 0.4635i-0.2676 + 0.4635i0.5352节点3 Uabc0000节点4 Uabc0.5489+0.3169i-0.0000-0.6338i-0.5489+0.3169i0.6338节点5 Uabc0.5977 +0.3451i0-0.6901i-0.5977+0.3451i0.6901节点1和2 Iabc0+0.8451i0.7319-0.4225i-0.7319-0.4225i0.8451节点1和3 Iabc0-4.5070i-3.9032+2.2535i3.9032+2.2535i4.5070节点2和3 Iabc0-5.3521i-4.6351+2.6761i4.6351+2.6761i5.3521节点4和1 Iabc0-3.6620i-3.1714+1.8310i3.1714+1.8310i3.6620节点5和2 Iabc0-6.1972i-5.3669+3.0986i5.3669+3.0986i6.1972 经A相接地ABC绝对值节点1 Uabc0.5822-0.3983-0.8660i-0.3983+0.8660i0.5822/0.9532/0.9532节点2 Uabc0.6610-0.3994-0.8660i-0.3994+0.8660i0.6610/0.9537/0.9537节点3 Uabc0.0000-0.5543-0.8660i-0.5543+0.8660i0.0000/1.0282/1.0282节点4 Uabc0.6661+0.5000i0.0000-1.0000i-0.6661+0.5000i0.8329/1.0000/0.8329节点5 Uabc0.6968+0.5000i0.0000-1.0000i-0.6968+0.5000i0.8577/1.0000/0.8577节点1和2 Iabc0+0.6608i0+0.6608i-0.0000-0.1384i0.6608/0.1384/0.1384节点1和3 Iabc0-4.3317i0.0000-0.0692i-0.0000-0.0692i4.3317/0.0692/0.0692节点2和3 Iabc0-4.9924i0+0.0692i0+0.0692i4.9924/0.0692/0.0692节点4和1 Iabc0-3.6709i0.0000-0.2077i-0.0000-0.2077i3.6709/0.2077/0.2077节点5和2 Iabc0-5.6532i-0.0000+0.2077i0.0000+0.2077i5.6532/0.2077/0.2077A相经10欧姆电阻接地短路ABC绝对值节点1 Uabc0.6456-0.1499i-0.4138-0.8296i-0.4138+0.9025i0.6628/0.9270/0.9928节点2 Uabc0.7125-0.1216i-0.4147-0.8299i-0.4147+0.9021i0.7228/0.9278/0.9929节点3 Uabc0.1518-0.3588i-0.5460-0.8855i-0.5460+0.84660.3896/1.0403/1.0074节点4 Uabc0.6964+0.4283i0.0000-1.0000i-0.6964+0.5717i0.8176/1.0000/0.9011节点5 Uabc0.7225+0.4393i-0.0000-1.0000i-0.7225+0.5607i0.8456/1.0000/0.9146节点1和2 Iabc-0.2371+0.5605i0.0497-0.1174i0.0497-0.1174i0.6085/0.1275/0.1275节点1和3 Iabc1.5542-3.6742i0.0248-0.0587i0.0248-0.0587i3.9894/0.0638/0.0638节点2和3 Iabc1.7913-4.2346i-0.0248+0.0587i-0.0248+0.0587i4.5979/0.0638/0.0638节点4和1 Iabc1.3172-3.1137i0.0745-0.1761i0.0745-0.1761i3.3808/0.1913/0.1913节点5和2 Iabc2.0284-4.7951i-0.0745+0.1761i-0.0745+0.1761i5.2065/0.1913/0.1913 bc两相相间短路ABC绝对值节点1 Uabc1.0000-0.5000-0.3903i-0.5000+0.3903i1.0000/0.6343/0.6343节点2 Uabc1.0000-0.5000-0.4635i-0.5000+0.4635i1.0000/0.6818/0.6818节点3 Uabc1.0000-0.5000-0.50001.0000/0.5000/0.5000节点4 Uabc0.8660+0.3169i-0.0000-0.6338i-0.8660+0.3169i0.9222/0.6338/0.9222节点5 Uabc0.8660+0.3451i-0.0000-0.6901i-0.8660+0.3451i0.9322/0.6901/0.9322节点1和2 Iabc0+0.0000i0.7319-0.0000i-0.7319-0.0000i0.0000/0.7319/0.7319节点1和3 Iabc0+0.0000i-3.9032-0.0000i3.9032+0.0000i0.0000/3.9032/3.9032节点2和3 Iabc1.7913-4.2346i-0.0248+0.0587i-0.0248+0.0587i4.5979/0.0638/0.0638节点4和1 Iabc1.3172-3.1137i0.0745-0.1761i0.0745-0.1761i3.3808/0.1913/0.1913节点5和2 Iabc2.0284-4.7951i-0.0745+0.1761i-0.0745+0.1761i5.2065/0.1913/0.1913bc两相短路接地ABC绝对值节点1 Uabc0.9036-0.3500-0.3903i-0.3500+0.3903i0.9036/0.5243/0.5243节点2 Uabc0.9046-0.3870-0.4635i-0.3870+0.4635i0.9046/0.6038/0.6038节点3 Uabc1.05150.0000-0.0000i0.0000+0.0000i1.0515/0.0000/0.0000节点4 Uabc0.7712+0.3169i-0.0000-0.6338i-0.7712+0.3169i0.8338/0.6338/0.8338节点5 Uabc0.7858+0.3451i-0.0000-0.6901i-0.7858+0.3451i0.8582/0.6901/0.8582节点1和2 Iabc0+0.1313i0.7319-0.2477i-0.7319-0.2477i0.1313/0.7726/0.7726节点1和3 Iabc0+0.0657i-3.9032+2.0872i3.9032+2.0872i0.0657/4.4262/4.4262节点2和3 Iabc0-0.0657i-4.6351+2.3349i4.6351+2.3349i0.0657/5.1899/5.1899节点4和1 Iabc0+0.1970i-3.1714+1.8394i3.1714+1.8394i0.1970/3.6662/3.6662节点5和2 Iabc0-0.1970i-5.3669+2.5826i5.3669+2.5826i0.1970/5.9560/5.95602 静态稳定性计算过程流程图输入系统参数，设定初值Pmax=0，Ke=0.1While (KePmax,将P，Ke，delta值赋给Pmax，Ke0，deltamax输出Pmax，Ke0，deltamaxKe=Ke+0.1结束静态稳定性计算结果：3Ke=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_best=2.100000, delta_best=97, P_max=1.623379储备系数Kp=0.623379.3.暂态稳定性计算过程流程图根据系统正常运行方式计算所需初值。i=0判断系统有无扰动，根据扰动性质选择功率特性曲线。D_Delta1=(Omega(i)-1)*360*f D_Omega1=(PT-P2M*sin(Delta(i)*pi/180)/TjDelta0(i+1)=Delta(i)+D_Delta1*hOmega0(i+1)=Omega(i)+D_Omega1*h D_Delta2=(Omega0(i+1)-1)*360*fD_Omega2=(PT-P2M*sin(Delta0(i+1)*pi/180)/Tj Delta(i+1)=Delta(i)+(D_Delta1+D_Delta2)*0.5*hOmega(i+1)=Omega(i)+(D_Omega1+D_Omega2)*0.5*h Omega(i+1)=Omega(i)+(D_