1、附录 A（流程图）输入数据根据 fault 值选择故障类型用节点导纳阵求逆得到节点阻抗矩阵选择短路点f计算各序电流的故障分量计算故障点各序电流计算故障点各序电压的故障分量计算故障后各节点电压值形成节点导纳矩阵确定 fault 值开始计算支路电流计算计及发电机端点的各节点电压结束【流程图一.短路计算】否否否是是是开始清屏设初值：Xd_all,Xdd_all,Xq,Xq_all,Eq0,UG0,Xe,Xd,XddKe=0.1Ke_result=Ke,delta_result=delta, K1_result= K1,Kemax_result=KemaxP_result=Eq*U*sin(delta
2、*pi/180)/Xd_all,delta=delta+0.01delta=49.00求 Kemax,Kemin求 Eq,Eqq求 UGd,UGq,UG求 K1K6Ke0)&(Ke=Kemin)&(Ke=Kemax)结束打印 Ke_result，delta_result, K1_result,Kemax_result【流程图二.暂态稳定性】是否开始结束赋初值（delta,omega, 步长h,时段长度Duration等）等输入故障切除时间 CutTimei=0i=round(CutTime/h)求时段初变化率求时段末估计值求时段末变化率求平均变化率求时段末计算值i=i+1显示此时功角同理求故障切除后的功角变化过程画出摇摆曲线【流程图三.静态稳定性】
