目录1.设计题目2.设计要求3.基础资料用于仿真程序验证4.设计内容三、设计内容1.仿真程序流程图3.仿真程序使用说明四、设计总结1.设计过程中遇到的问题及解决方法2.设计体会3.对设计的建议五、参考文献1.设计题目两相短路故障短路电流仿真程序设计2.设计要求采用计算机算法计算故障点不对称故障电流、短路功率等参数。需绘制计算机算法的流程图并基于流程图编写相应计算程序。3.基础资料用于仿真程序验证电力系统图1所示输电线路、电力变压器、系统、发电机及负载的参数均标注于图中系统图f点发生不对称短路b、c两相短路分析与计算短路电流与短路功率不对称故障的短路功率为故障电流有效值。变压器连接组标号均为YNd11。各元件三序参数可参考后面附表。电动机三序参数可参考教材表10.1或者按照祝淑萍的书135页略去不计。距离较远电动机电势对短路电流的影响可忽略。本题复合序网只有正序和负序故计算元件参数时可只计算正序和负序参数。关于线路的几何均距参考后面表格自己选择即可。4.设计内容1设计两相短路故障计算程序流程图。2编写基于节点导纳矩阵的系统两相短路相关参数计算程序该程序应实现以下主要功能①基于网络各序标幺值等值参数采用导纳矩阵求逆法求取系统的正、负序节点阻抗矩阵②针对任意给定的故障点计算其在发生任意两相短路时的故障点正、负序电流、故障点短路电流以及短路功率。③在故障点发生两相短路的条件下计算任意指定支路的各相电流和任意指定节点的各相电压。3以3基础资料所述系统为例以100 MVA为基准容量采用平均额定电压作为基准电压计算各元件的正、负序标幺值等值参数据此绘制标注各节点序号的系统正、负序等值电路(请自行对各节点进行编号)采用所编写程序对其仿真分析验证程序正确性。三、设计内容1.仿真程序流程图2.仿真程序源代码%%节点阻抗矩阵函数 % B Z 支路参数矩阵 n1*6 % Y 节点导纳矩阵 function [Z]T3(B,X) nllength(B(:,1)); nlength(X(:,1)); Yzeros(n,n); for i1:nl if(B(i,6)1) LB(i,1); RB(i,2); else LB(i,2); RB(i,1); end Y(L,L)Y(L,L)1/B(i,3)*B(i,5)^2B(i,4)/2; Y(R,R)Y(R,R)1/B(i,3)B(i,4)/2; Y(L,R)Y(L,R)-1/B(i,3)*B(i,5); Y(R,L)Y(L,R); end for i1:n Y(i,i)Y(i,i)X(i,2); end Zinv(Y); end B[1 2 0.058i 0 0.9565 1; 2 3 0.329i 0 1 0; 3 4 0.1665i 0 0.9565 0; 3 6 0.161i 0 1 0; 6 5 1.05i 0 1.00435 0]; X1[1 1/0.204i; 2 0; 3 1/6.321i; 4 1/0.208i; 5 1/8.5i; 6 0]; X2[1 1/0.204i; 2 0; 3 1/3.321i; 4 1/0.208i; 5 1/8.5i; 6 0]; %正序阻抗矩阵 Z1T3(B,X1); %负序阻抗矩阵 Z2T3(B,X2); %正负序电流 I_11.12/(Z1(3,3)Z2(3,3)) I_2-I_1 %短路电流标幺值 有名值 If_1abs(I_1)*sqrt(3) IfIf_1*100/sqrt(3)/115*1000 %短路功率 P100*If_1 %输入节点位置 i3; j2; %%各节点电压 %正序 U_01.12; U_i_1U_0-Z1(i,3)*I_1 Zjk_10.0887i; U_j_1U_0-Z1(j,3)*I_1 %负序 U_i_2-Z2(i,3)*I_2 Zjk_20.0862i; U_j_2-Z2(j,3)*I_2 %%各支路电流正序 %正序 I_ij_1(U_i_1-U_j_1)/0.329i %负序 I_ij_2(U_i_2-U_j_2)/0.329i %%各相电流 a-0.51i*sqrt(3)/2; I_aI_ij_1I_ij_2 I_ba^2*I_ij_1a*I_ij_2 I_ca*I_ij_1a^2*I_ij_2 %%各相电压 U_aU_i_1U_i_2 U_ba^2*U_i_1a*U_i_2 U_ca*U_i_1a^2*U_i_2 %%判断支路是否存在变压器 %t是判断条件 t0; % 角度转弧度 theta_pos deg2rad(30); % 正序30°YNd11 theta_neg deg2rad(-30); % 负序-30° if(t1) disp(支路存在变压器从高压侧到低压侧); %正序参数相位偏移 U_i_1U_i_1*exp(1i*theta_pos) U_j_1U_j_1*exp(1i*theta_pos) I_ij_1I_ij_1*exp(1i*theta_pos) %负序参数相位偏移 U_i_2U_i_2*exp(1i*theta_neg) U_j_2U_j_2*exp(1i*theta_neg) I_ij_2I_ij_2*exp(1i*theta_neg) %%各相电流 a-0.50.866i; I_aI_ij_1I_ij_2 I_ba^2*I_ij_1a*I_ij_2 I_ca*I_ij_1a^2*I_ij_2 %%各相电压 U_aU_i_1U_i_2 U_ba^2*U_i_1a*U_i_2 U_ca*U_i_1a^2*U_i_2 end3.仿真程序使用说明(一)、程序功能该MATLAB仿真程序用于电力系统不对称BC两相短路故障计算依托自定义T3函数实现导纳矩阵求逆法求解序网节点阻抗矩阵。程序先录入电网支路、节点对地序参数分别搭建正序网络、负序网络节点导纳矩阵对导纳矩阵求逆得到正负序节点阻抗矩阵再输入故障节点编号根据两相短路故障边界条件计算故障点正序、负序电流标幺值推导故障短路电流、短路容量接着遍历全部网络节点计算所有节点的正序、负序序电压提取任意支路的正负序支路电流最后利用对称分量变换公式将序电压、序电流转换为实际A、B、C三相电压与三相电流完整实现含变压器、发电机、线路的多节点电网不对称短路数值仿真可用于课程设计故障数值验算与电网故障特性分析。(二)、输入参数、输入矩阵及各元素含义1.支路参数矩阵B该矩阵为核心拓扑输入矩阵每一行对应电网一条电气支路固定六列参数。第一列代表支路起始节点编号第二列代表支路终止节点编号用来确定电网各元件的拓扑连接关系第三列是支路串联电抗的复数标幺值本次仿真无电阻仅保留虚数电抗第四列为线路对地充电电纳本次模型不考虑电容全部填充0第五列是变压器变比标幺值普通输电线路无变压器统一填1变压器支路填入折算后的标准变比第六列为接地标志数字数值1代表该支路连接接地元件数值0代表支路两端均不接地。2.节点对地正序导纳矩阵X1矩阵为两列多行结构用来录入每个节点对地的正序电气参数。每行第一列是节点编号第二列是该节点对地元件的正序导纳标幺值发电机节点填入发电机正序电抗倒数无接地元件的线路、负荷节点直接填0用于函数生成完整正序导纳矩阵。3.节点对地负序导纳矩阵X2矩阵格式与X1完全一致仅第二列为节点对地负序导纳标幺值发电机、电动机正负序电抗数值不同对应节点导纳数值会发生变化无源节点依旧填充0用于构建负序网络导纳矩阵。4.单独输入常量参数包含系统基准容量、基准电压、发电机空载电动势标幺值、故障节点编号以及对称分量变换计算所需的复数算子a若网络含YNd11变压器还需输入相位校正系数用于修正变压器带来的序分量相位偏移。(三)、输出参数、输出矩阵及各元素含义1.正序节点阻抗矩阵Z1、负序节点阻抗矩阵Z2二者均为和网络节点数量一致的方阵矩阵行号、列号均对应电网节点编号。矩阵对角线上的元素Z(k,k)为第k个节点的序自阻抗物理含义是仅在k节点通入单位电流时k节点自身产生的序电压非对角线元素Z(m,n)为m节点与n节点之间的序互阻抗代表n节点通入单位电流时m节点产生的序电压完整反映网络各节点之间的电气耦合关系。2.故障点电气量输出参数输出故障点正序电流、负序电流复数标幺值BC两相短路工况下二者幅值相等输出故障短路电流标幺值与有名值、短路容量标幺值直观反映短路故障严重程度。3.节点、支路序电气量输出参数输出全网每个节点对应的正序电压、负序电压复数标幺值输出任意指定支路流过的正序支路电流、负序支路电流复数标幺值复数同时包含电气量幅值与相位信息。4.三相电气量输出参数依托对称分量公式换算得到A相、B相、C相三相电压、三相电流复数可直接提取幅值与相位直观体现不对称短路带来的三相电压跌落、三相电流不平衡问题存在变压器时额外输出相位校正后的修正版三相电压、三相电流消除变压器绕组相位带来的计算误差。(四)、自定义函数T3输入输出补充说明函数输入第一个输入B为支路参数矩阵第二个输入X为对应序网的节点对地导纳矩阵X1或X2函数输出Z为对应序网完整节点阻抗方阵矩阵行列对应系统各节点编号。函数内部逻辑先搭建节点导纳矩阵逐条录入支路的自导、互导叠加对地电容与发电机对地导纳最后对导纳矩阵求逆得到节点阻抗矩阵输出。4.仿真程序验证四、设计总结1.设计过程中遇到的问题及解决方法1、在绘制标注各节点序号的系统正负序等值电路时有问题。解决办法通过和老师及交流和查询课本例题找到错误。2、仿真时节点导纳矩阵等于零写入的数据有误矩阵X第二列接地导纳值输入有问题。解决办法灵活运用Chatgpt 、deepseek和Gemini等等以及和老师交流查询课本资料找到错误并修改。3、计算相关变压器支路时没考虑到相位偏移。解决办法通过查询教材上公式和资料得知正负序电压电流经过星三角形变压器的相位偏移情况。2.设计体会本次两相短路故障仿真程序设计是我第一次完整经历“理论推导→参数计算→程序编写→结果校验”的电力系统工程实践收获远超预期。最初我以为只是把课本公式敲进MATLAB实际却踩了不少“隐形坑”。比如绘制正、负序等值电路时我一开始忽略了变压器YNd11接线对相位的影响导致支路电流计算结果与理论值有偏差后来在搭建节点导纳矩阵时又误将节点对地导纳填为0导致矩阵奇异、程序报错反复排查后才发现是发电机次暂态电抗的导纳值漏写。这些问题让我深刻意识到电力系统计算容不得半点马虎一个参数的疏漏就会导致整个结果失效。最让我有成就感的是对“仿真结果与理论边界条件差异”的分析。一开始我发现短路点a相电压接近0和课本中“两相短路a相电压不为零”的结论矛盾一度以为是程序错误。后来通过查阅《电力系统故障分析》才明白课本结论是“理想不接地系统”的假设而实际工程中110kV系统多为中性点接地短路点电位会被钳位到地电位。这次“纠错”让我真正打通了“理论-仿真-工程”的逻辑链条仿真不是对理论的简单复刻而是要反映真实系统的复杂性。此外我也体会到模块化编程的优势。把导纳矩阵构建、阻抗矩阵求逆、序分量合成拆分成独立模块后调试效率提升了不少也更理解了计算机算法在电力系统大规模计算中的必要性——如果用手算这个6节点系统至少需要半天而程序几秒就能出结果。这次设计不仅巩固了不对称故障分析的知识更培养了我的工程严谨性和问题解决能力为后续的毕业设计和工程实践打下了坚实基础。3.对设计的建议1、 增加“理论-仿真差异”的引导分析本次设计中“短路点a相电压接近0”的问题困扰了我很久后来才明白是理想模型与实际系统的差异。建议在教学中专门增设一节“仿真结果合理性分析”引导学生区分“教科书理想假设”和“工程实际场景”比如讲解中性点接地方式对短路电压的影响避免学生把“数值差异”简单归为“程序错误”培养工程思维。2、 优化参数查询的便利性本次设计需要查询线路电抗、变压器阻抗等参数附录中的表格排版较乱查找效率低。建议将常用参数整理成电子表格按电压等级、设备类型分类索引或提供在线参数查询工具减少学生在非核心任务上的时间消耗。3、增加“变压器相位偏移”的专项练习YNd11变压器的±30°相位偏移是不对称故障计算的重点和难点本次设计中我差点遗漏这一环节。建议在设计任务中明确要求对比“考虑相位偏移”和“不考虑相位偏移”的结果差异加深学生对变压器接线组别影响的理解。4、开放“扩展故障类型”的选做任务本次设计仅要求实现两相短路仿真对学有余力的同学可以开放单相接地、两相接地短路等扩展任务鼓励对比不同故障类型的序网连接方式和计算结果进一步提升综合能力。五、参考文献[1] 何仰赞,温增银.电力系统分析[M].4版.武汉:华中科技大学出版社,2016.[2] 陈珩.电力系统稳态分析[M].4版.北京:中国电力出版社,2018.[3] 李光琦.电力系统暂态分析[M].3版.北京:中国电力出版社,2012.[6] 刘万顺.电力系统故障分析[M].3版.北京:中国电力出版社,2010.[7] DL/T 5429-2009 电力系统设计技术规程[Z].北京:中国电力出版社,2009.[9] 电力工业部电力规划设计总院.电力系统设计手册[Z].北京:中国电力出版社,1998.[10] 孟祥萍高燕.电力系统分析.3版. 北京:高等教育出版社2021.5