单机无穷大系统发生各类(三相短路，单相接地，两相接地，两相相间短路)等短路故障，各类(单相断线...
单机无穷大系统发生各类(三相短路，单相接地，两相接地，两相相间短路)等短路故障，各类(单相断线，两相断线，三相断线)等断线故障，暂态稳定仿真分析
电力系统仿真这玩意儿说难不难，但碰到短路断线这些幺蛾子总能让人抓狂。今天咱们直接上手单机无穷大系统，拿Python搞点硬核的暂态稳定仿真。别慌，代码我帮你踩过坑了。
先说短路故障这块，三相短路最要命。系统导纳矩阵突变相当于电网"心脏骤停"，这时候发电机的转子摇摆方程直接进入狂暴模式：
def swing_equation(delta, t, P_m, E, V, Xd, H): d_delta = delta[1] P_e = (E * V) / Xd * np.sin(delta[0]) d_omega = (P_m - P_e) / (2 * H) return [d_delta, d_omega] Xd_fault = Xd / 2 sol = odeint(swing_equation, [delta0, 0], t, args=(Pm, E, V, Xd_fault, H))注意看第5行的Pe计算，故障时电抗Xd突然变小会导致电磁功率暴增。这时候如果机械功率Pm跟不上，转子加速度直接起飞。仿真结果里功角曲线要是突破120度，系统基本就凉了。
接地故障更复杂点，得用对称分量法。比如单相接地得玩零序网络：
% 正序网络保持原样 Z1 = 0.1; % 零序网络参数突变 Z0 = 0.3 + 0.4i; I_fault = 3 * V_prefault / (2*Z1 + Z0 + 3*Zf);这段MATLAB代码里的3倍系数是精髓，不对称故障会产生零序电流，保护装置能不能及时切除就看这个电流值了。仿真时记得故障持续时间设短点，0.2秒是临界值。
断线故障仿真更带劲。比如两相断线，系统直接变不对称运行。这时候得用节点阻抗矩阵重构：
# 原始导纳矩阵 Ybus = np.array([[1/(0.1j), -1/(0.1j)], [-1/(0.1j), 1/(0.1j)]]) # B相断线时修改矩阵 Ybus[1,1] -= 1/(0.5j) # 断线相导纳归零注意看第5行的操作，断线相当于该相导纳变成无穷大（阻抗无穷小），但实际操作中要避免数值计算问题，改成加个极大值更稳。仿真时观察三相电压波形，断线相电压会飙升到1.5倍以上。
暂态稳定判据别死记教科书那套，实战中看这三个指标：
- 功角曲线最大摆幅（别超120度）
- 转速偏差恢复情况（±0.5Hz内）
- 母线电压恢复时间（故障切除后20周波内）
搞仿真的兄弟注意了，故障时序设置是门艺术。比如两相接地转单相断线这种组合故障，时间步长设1ms以下才抓得到暂态细节。用变步长算法时，记得锁住故障发生前后1个周期的时间步。
最后说个骚操作：三相断线时别急着重启动，先看发电机有没有进入孤岛运行模式。这时候调频特性决定生死，仿真代码里调频系数设太小的话，频率能给你飙到55Hz去。所以备用容量参数千万别瞎设，血泪教训啊！
仿真这玩意儿说到底就是和微分方程斗智斗勇，下次碰到仿真发散别急着砸电脑，先把故障持续时间减半试试。稳住，咱们能赢！