5.1 功角稳定性：同步的拉锯战

5.1 功角稳定性：同步的拉锯战

本节摘要：功角稳定研究发电机受扰后能否保持同步。静态稳定由功角特性曲线的上升段裕度决定，暂态稳定由等面积定则判定。本节用欧拉积分复现一次完整的暂态稳定仿真，找出临界切除时间。

摆回来的，和摆丢的

回到 2.1 节建立的图像：并网发电机之间靠转子磁场的"咬合"保持同步，功角是这种咬合的松紧度量。正常运行时机械功率与电磁功率平衡，功角停在某个稳定值。短路发生的瞬间，发电机端电压跌落、送出的电磁功率骤减，可原动机的蒸汽推力不会立刻变——多出来的功率全部变成转子加速的动能，功角开始增大。故障切除后，端电压恢复，功角摆到比原来大的位置时电磁功率反过来超过机械功率，转子减速。接下来就是一场拉锯：动能与位置量的交换，像弹簧上的小球。摆幅逐次衰减回到原点，叫暂态稳定；一摆冲过临界角度再也回不来，发电机与系统失步，叫暂态失稳——失步后的发电机平均发出异步功率、巨大电流冲击机组，必须立即解列。

它在体系中的位置：本节是第五章三条边界的第一条，向上直接使用 2.1 的功角概念与 4.2 的故障场景，向下的"切除越快越稳"是第六章保护速动性的物理出处，而断面稳定限额则是 9.4 节安全标准的计算内核之一。

功角特性：一条正弦曲线上的两个世界

单机对无穷大系统的情形最能看到本质：送出有功等于两电势幅值乘积除以联系电抗，再乘功角的正弦。这条正弦曲线被三十度到九十度的分界切成两个世界。三十度附近（曲线上升段，斜率为正）：功角被扰动推大一点，电磁功率随之增大，把功角拉回去——像碗底的小球，这是静态稳定的运行区。九十度（曲线顶点）：电磁功率达到最大，是数学上的静稳极限；九十度以后（下降段）：功角增大功率反而减小，加速度更大，失稳雪崩式发展——碗沿翻了过来。工程上当然不会贴着九十度运行，实际功角常压在三四十度以内，留出承受扰动的裕度；所谓"稳定限额"本质就是"输送功率加上扰动还够不着顶点"的那条线。

从这条曲线还能读出两个工程手段的原理：串联补偿减小联系电抗，直接抬高整条曲线的顶点，同样功率对应的功角更小——这就是 3.2 节串补与 9.2 节灵活直流的共同原理"缩短电气距离"；快速励磁强行增大发电机电势，也抬曲线、扩裕度。

静稳与暂稳是同一座山的两条上山路径，值得把关系挑明。静态稳定回答"小扰动下能否回到原点"——运行点在曲线上升段且离顶点有距离就行；暂态稳定回答"大扰动后能否被拉回来"——能量法则说了算。二者共享同一条功角曲线，但问的问题不同：前者是局部的线性化问题（曲线斜率即同步转矩系数），后者是全局的能量问题（加速与减速面积的积分）。工程上"静稳储备系数"（输送功率距静稳极限的百分比）与"暂态稳定裕度"（临界切除时间与实际切除时间的差）分别量化两条边界，方式单上"按暂稳限额控制断面"的条款，背后就是这两个指标的计算结果——第五章的知识到运行规程里，就浓缩成这几个数字。

曲线上的拉锯：等面积定则

等面积定则：不用解方程的判据

暂态稳定本质要解摇摆方程（二阶非线性微分方程），但首摆稳定性有个优雅的几何判据。在功率-功角平面上：故障期间运行点沿被压低的故障曲线走，机械功率线与故障曲线围成的面积是转子吸收的加速能量；故障切除后运行点跳回正常曲线，正常曲线与机械功率线围成的是可释放的减速能量上限。首摆稳定的充要条件：加速面积小于最大可用减速面积。一眼就能看出两个规律：故障切除越晚，功角被推得越远，加速面积越大；输送功率越重，初始功角越大、可用减速面积越小——重载方式稳定裕度小，这就是方式安排要控制断面功率的第一个原因。

亲手复现一次暂态仿真

等面积定则给判据，数值积分给全过程。取单机无穷大系统：机械功率零点九标幺，正常功角特性幅值一点八（初始功角恰为三十度），故障期间幅值被压到零点四；惯性常数五秒，阻尼系数二。用欧拉法积分摇摆方程，比较不同切除时间的结局：

import math

Pm, Pmax_n, Pmax_f = 0.9, 1.8, 0.4      # 机械功率、正常/故障期间幅值（标幺）
H, D, ws = 5.0, 2.0, 2 * math.pi * 50   # 惯性常数、阻尼、同步角速度
delta0 = math.asin(Pm / Pmax_n)          # 初始功角：三十度

def simulate(tc, T=3.0, dt=0.0005):
    d, s, t, peak = delta0, 0.0, 0.0, delta0
    while t < T:
        Pe = (Pmax_f if t < tc else Pmax_n) * math.sin(d)
        s += (Pm - Pe - D * s) * dt / (2 * H)   # 转差（标幺）
        d += ws * s * dt                          # 功角（弧度）
        peak = max(peak, d)
        if d > math.pi:
            return "失稳", math.degrees(peak)
        t += dt
    return "稳定", math.degrees(peak)

for tc in (0.10, 0.20, 0.30, 0.35):
    state, peak = simulate(tc)
    print(f"切除时间 {tc:.2f} 秒 -> {state}，最大功角 {peak:.0f} 度")
​

运行结果：切除零点一秒，最大功角只摆到四十九度即被拉回；零点二秒摆到七十四度；零点三秒摆到一百一十四度，仍在回摆的边缘惊险稳住；零点三五秒则一摆冲过一百八十度，判定失稳。临界切除时间落在零点三到零点三五秒之间——把这个数字与第六章保护的动作时间（二到五个周波，即四十到一百毫秒）对照，你会看到保护速动性的全部意义：保护每快一个周波，稳定裕度就厚一分。改参数再试也有收获：把输送功率加到一点二（初始功角约四十二度），临界切除时间显著缩短；投入强力励磁抬高故障后曲线，临界时间延长——教材里的定性结论，在这段代码里都变成可复现的定量曲线。

这段仿真是工程稳定计算的"细胞"。真实的暂态稳定程序无非是三个方向的扩展：机组从一台到上千台（摇摆方程组耦合求解）、模型从经典二阶到包含调速励磁的详细模型（每台机十几个状态量）、场景从单一故障到数千个预想故障的批量扫描。扫描的目的正是找"最紧的边界"——哪个方式、哪个故障、哪台机组最接近失稳，稳定限额就从那里反推出来。你在方式单上看到的每一条断面限额，背后都是一晚上机器轰鸣的批量仿真；而读懂那些结果报告的能力，起点就是本节这一百行不到的循环。

⚠️ 工程与教科书的差距：真实系统是多机系统，"一台机对无穷大"的图像只是入门。多机稳定的失步形态复杂得多（同调机群分群、区间振荡），工程仿真用逐步积分法算全网的机-网耦合方程，单机版的等面积直觉仍是读懂仿真结果的底色。

要点打包带走

• 两个世界：功角特性上升段稳定、下降段发散，九十度是数学极限，工程运行压在三四十度留裕度。
• 等面积定则：加速面积对减速面积的赛跑，几何上解释了"快切除"与"轻方式"为何都增稳。
• 临界切除时间：稳定与否的分水岭，由仿真扫描得出，是保护整定与稳定限额的共同依据。
• 增稳组合：串补与直流缩短电气距离、强励抬曲线、快关汽门减机械功率，全部作用在两张面积的消长上。
• 静稳与暂稳：一个问小扰动能否回原点、一个问大扰动能否被拉回，同一曲线上的两类问题、两套指标。

第二条边界看的不是转子而是电压：无功耗尽时，电压如何沿着鼻形曲线滑向崩溃，是下一节的内容。