Simulink单机无穷大系统两相接地短路暂态稳定仿真

1. 为什么单机无穷大系统是理解暂态稳定的最佳模型

1.1 从“单机无穷大”的假设说起：物理意义与建模逻辑

刚接触电力系统暂态稳定仿真的人，最容易被各种复杂模型吓住——多机系统、网络化简、动态等值，听着就头大。但几乎所有教材、几乎所有科研入门工作，都会从“单机无穷大系统”开始。这不是偷懒，而是这个模型恰好抓住了暂态稳定问题的核心矛盾。

所谓“单机无穷大系统”，就是一台同步发电机经过变压器和输电线路，连接到容量无限大的电力系统。这个“无穷大”的物理含义是：无论这台发电机怎么折腾——加速、减速、摇摆、甚至短时失步——系统侧的电压幅值和频率都纹丝不动。相当于你在一艘万吨巨轮旁边踩了一脚，巨轮不会有任何反应。这个假设把问题边界变得清爽：所有暂态行为的研究焦点，全部集中在这台发电机自身的机电暂态过程上。

用Simulink做这类仿真，首先要理解的不是Simulink怎么操作，而是这个模型背后的三套方程：发电机的转子运动方程（描述转速和功角变化）、发电机内部的电磁暂态方程（描述定子、转子绕组的电磁过程）、以及网络方程（描述发电机通过变压器和线路与系统之间的功率传输关系）。在Simulink里，你可以选择把这三套方程全部用元件模块搭出来，也可以自己写状态方程塞进一个S-Function或MATLAB Function模块里。两条路线我都走过，各有优劣，后面会详细对比。

有一点需要提前说明：Simulink里的电力系统仿真，本质上是用数值积分算法去解一个高阶微分方程组。所以你先要建立这个观念——你写的每一个参数、选的每一个求解器、设的每一步长，最终都会反映在微分方程求解的精度和稳定性上。很多新手仿真发散，不是模型错了，而是求解器设置不合理。

1.2 暂态稳定问题的本质：功角与转速的摇摆过程

要理解发电机转速变化曲线，必须先想明白一个物理过程：同步发电机在正常稳态运行时，转子以同步转速旋转，原动机输入的机械功率 ( P_m ) 与发电机输出的电磁功率 ( P_e ) 正好平衡。此时功角 ( \delta ) 保持恒定，转速 ( \omega ) 等于同步转速。

但故障发生时，一切平衡被打破。以两相接地短路为例：故障瞬间，发电机的电磁功率 ( P_e ) 急剧跌落——因为电压跌落、故障支路分流，送出去的功率大幅减少。但原动机那边的机械功率 ( P_m ) 可来不及变，汽轮机或水轮机的调速器都有惯性，短时间内机械功率基本维持原值。这一下就出现了功率差额 ( \Delta P = P_m - P_e > 0 )。发电机转子在这个不平衡力矩的驱动下开始加速，转速升高。

等故障切除后，输电线路恢复，电磁功率 ( P_e ) 又弹回来了。但如果故障期间转子已经积累了足够大的动能，功角冲得很远，那么故障切除后 ( P_e ) 可能依然大于 ( P_m )（在减速区间），转子开始减速。这个过程就像荡秋千——加速冲上去，减速荡回来，功角和转速呈现周期性的摇摆。如果摇摆的幅度不断衰减，最终稳定在新的平衡点或回到原平衡点，系统是暂态稳定的；如果功角越摆越大，转速越偏越多，最终失去同步，就是暂态失稳。

这就是你在Simulink示波器里看到的那条转速变化曲线背后的全部物理意义。所以这篇文章虽然是讲仿真实现，但核心还是帮你把物理过程和仿真现象一一对应起来，否则你仿真跑出来了也看不懂波形在说什么。

2. 两相接地短路故障：比想象中更有“学问”的扰动

2.1 两相接地短路的故障特征与复合序网

很多人做暂态稳定仿真，上手就是三相短路。这当然最简单，但现实电网里，单相接地短路占比最高，两相接地短路和两相相间短路也时有发生。两相接地短路（比如B、C两相同时接地）的物理过程比三相短路复杂得多，因为这是不对称故障。

不对称故障的分析武器是对称分量法。把故障点的电压、电流分解为正序、负序、零序三组对称分量。两相接地短路的边界条件是：故障处B、C相电压为零，A相电流为零。由这个边界条件推出的复合序网是正序、负序、零序三个网络在故障点处并联。如果你用Simulink里的Three-Phase Fault模块来设置故障类型，模块底层已经替你处理好了这些序网连接，你只需要在下拉框里选“Phase B-Ground”和“Phase C-Ground”同时选中就行。

但这里有个关键点：两相接地短路的负序和零序分量会在发电机定子绕组中产生负序电流和零序电流，这些电流会产生反向旋转的磁场和脉动转矩，作用在发电机转子上，引起二倍频的振荡分量叠加在转速响应上。所以你会在转速波形上看到，两相接地短路情况下的转速曲线，比三相短路情况的曲线多了一层高频“毛刺”或波纹。这不是数值噪声，而是负序分量真实作用的物理反映。看波形时不要把它误认为求解误差。

我在最初做这个仿真时，就犯过这个错。当时转速波形上有明显的二倍频振荡，我以为是步长太大导致数值不稳定，拼命减小步长，结果毛刺还在。后来查了文献才明白，这就是不对称故障的负序效应。所以做这类仿真，务必对“哪些现象是物理真实、哪些是数值引入”有基本判断力。

2.2 故障对发电机电磁功率的影响：等面积定则的视角

要预判转速会冲到多大、系统能不能维持稳定，最经典的工具是等面积定则。它可以在不做仿真的情况下，粗略估算出极限切除角——也就是最晚必须在功角摆到多大之前切除故障，否则系统失稳。

正常运行时，发电机送到无穷大系统的电磁功率为：

[ P_e = \frac{E'U}{X_{\Sigma}} \sin\delta ]

( E' ) 是暂态电抗后的电势，( U ) 是无穷大系统母线电压，( X_{\Sigma} ) 是发电机暂态电抗、变压器电抗、线路电抗之和。正常运行时 ( X_{\Sigma} ) 较小，( P_e ) 曲线较高。

发生两相接地短路时，故障点相当于经过一个不对称阻抗接地，复合序网中故障点的等效阻抗是 ( X_2 ) 与 ( X_0 ) 并联后再串入网络。此时从发电机看进去的等效总电抗 ( X_{\Sigma}^{(2)} ) 明显大于正常值。电抗一大，功角特性曲线的幅值就变低——同样一个功角下，能送出去的电磁功率变小了。这就是加速面积产生的原因。

有意思的是，两相接地短路时功率曲线虽然比正常低，但通常比三相短路时的功率曲线高。因为三相短路时故障点三相对称短路，从发电机看进去的等效电抗更大，功率曲线压得更低，加速面积更大，极限切除角更小。换句话说，三相短路是更严苛的故障类型，两相接地短路的暂态稳定裕度通常比三相短路大一些。所以做“两相接地短路”研究时，你会发现同样的故障切除时间下，转速摇摆的幅度比三相短路小，系统更容易保持稳定。这个认知对后续仿真参数设计很有用——你需要把故障切除时间设置得更接近临界值，才能观察到转速的显著变化。

等面积定则在Simulink里的价值不是用来精确计算，而是用来验证仿真结果。你先用等面积定则估算一个大概的极限切除角，再在仿真里扫描故障切除时间，理论上两者对应的临界切除时间是吻合的。如果差太远，要么等面积定则算错了，要么Simulink模型哪里有问题。这种交叉验证的做法，能让你的仿真实力提升一个台阶。

3. Simulink仿真实现：两种路线与选型建议

3.1 元件级建模路线：Simscape Electrical方案

第一条路线是用Simscape Electrical（早期叫SimPowerSystems）里的元件模块搭模型。这是大多数人首选的方案，因为模块化程度高，不需要手推微分方程。你需要用到的模块大致有：

• 同步发电机模块（Synchronous Machine pu Standard），选标幺制模型
• 原动机与调速器模块（比如Steam Turbine and Governor或Hydraulic Turbine and Governor，取决于你模拟的是哪种电厂）
• 变压器模块（Three-Phase Transformer Two Windings）
• 输电线路模块（Three-Phase PI Section Line），你需要用两回并联线路来模拟一个简化的输电通道
• 无穷大系统模块（Three-Phase Source，设大短路容量，相当于无穷大母线）
• 故障模块（Three-Phase Fault）
• 断路器模块（Three-Phase Breaker）或直接用故障模块本身的时序控制来模拟故障切除

这种搭法的好处是直观，物理结构跟实际系统一一对应，变压器、线路、发电机内部绕组结构全部可视化。而且模块内部自带完整的Park方程磁链模型，负序、零序效应、定子瞬变过程都会自动算出来，不需要你自己写方程，物理精度有保证。

代价是模型计算量大，仿真速度慢，而且参数设置繁琐。每个模块都有十几甚至几十个参数，每个参数都要填得合理，否则仿真结果就会很奇怪。我见过很多同学直接把同步发电机模块默认参数一放，跑出来的转速曲线跟心电图一样乱跳，根本没法用。

3.2 方程级建模路线：自建微分方程方案

第二条路线是用MATLAB Function或S-Function把发电机转子运动方程和电磁功率方程直接写出来。最核心的就是二阶转子运动方程：

[ \frac{2H}{\omega_0} \frac{d^2\delta}{dt^2} = P_m - P_e ]

或者写成两个一阶方程：

[ \frac{d\delta}{dt} = \omega - \omega_0 ]

[ \frac{2H}{\omega_0} \frac{d\omega}{dt} = P_m - P_e ]

( H ) 是发电机的惯性时间常数（单位秒），( \omega_0 ) 是同步转速（rad/s），( P_m ) 是机械功率，( P_e ) 是电磁功率。这里的 ( P_e ) 分三段按故障状态切换：正常时用正常电抗算，故障时用故障状态下的等效电抗算，故障切除后按切除后网络算。

这种做法的好处是轻量、快速、参数透明。修改一个电抗值就像改一个常量一样简单，物理过程的信息完全在你的掌控中。而且因为没有变压器、线路等模块的动态过程干扰，转速曲线的物理规律显示得特别干净，适合用来做理论验证和研究。

坏处是模型考虑的因素有限，比如负序电流引起的二倍频转矩、定子电磁暂态过程，这些在简化模型中体现不出来。而且如果你对电力系统分析理解不深，很容易在故障状态电抗计算上出错——比如两相接地短路时的等效电抗是负序和零序电抗并联后再与正序网络串联，少算一个并联步骤结果就差很远。

3.3 为什么不推荐直接拿三相短路改参数来模拟两相接地短路

顺便提醒一个常见操作误区。有人图省事，把三相短路故障模块改两相，然后其他参数不动——这是错的。两相接地短路和三相短路对应的复合序网完全不同，故障期间发电机的电磁功率曲线也不一样。如果你只是把故障类型从“三相”改成“两相接地”，而故障阻抗、接地阻抗、序网参数都没做相应调整，那你仿真出的转速曲线既不是两相接地的结果，也不能反映三相短路的特征，两头都不靠。

正确的做法是：故障模块里明确选择两相接地（B相和C相同时接地），并且设置合理的故障电阻 ( R_f ) 和接地电阻 ( R_g )。注意接地电阻的大小对零序网络影响很大，直接决定零序电流幅值，进而影响故障期间电磁功率跌落的程度。我建议一开始设置较小的故障电阻（比如0.001Ω），先把故障当作金属性短路来处理，再逐步加大故障电阻观察转速摇摆幅度的变化规律。

4. 仿真模型搭建全流程拆解

4.1 第一版：用Simscape Electrical搭出完整单机无穷大系统

如果走元件级路线，我建议按下面这个顺序搭建，这个顺序是我踩过不少坑后总结出来的，能最小化排查成本。

第一步，放置同步发电机模块。双击打开参数界面，在“Parameters”标签下选择“pu Standard”类型。额定参数设置为：额定视在功率 ( S_n = 100 )MVA，额定线电压 ( U_n = 13.8 )kV，额定频率 ( f_n = 50 )Hz。然后填暂态电抗 ( x'_d = 0.3 )pu、次暂态电抗 ( x''_d = 0.2 )pu、惯性时间常数 ( H = 3.5 )s 这些经典值。这里最关键的是惯性时间常数，它直接决定转速摇摆的周期和衰减特性—— ( H ) 越小转子越“轻”，故障时加速越快，转速偏移越大。

第二步，放置变压器。升压变压器容量100MVA，变比13.8kV/230kV，短路电抗 ( X_T = 0.1 )pu。这里要说一个细节：变压器模块中的“连接方式”选Yg-Delta还是Y-Y，决定了零序网络是否连通。做两相接地短路仿真，零序通路必须正确。一般推荐发电机侧用三角形接线，阻断零序电流进入发电机内部，而系统侧用中性点接地的星形接线，为零序电流构成通路。这个如果不设置对，故障电流和转速波形都会有问题。

第三步，放置双回线路。用两个Three-Phase PI Section Line模块并联，模拟双回线运行。每回线路长度100km，单位长度正序电阻 ( R_1 = 0.01273 )Ω/km，正序电感 ( L_1 = 0.9337 )mH/km，正序电容 ( C_1 = 0.01274 )μF/km，零序参数约为正序的3倍左右。这些参数可以直接用线路模块自带的默认值修改得到。

第四步，放置无穷大系统。用一个Three-Phase Source模块，电压230kV，频率50Hz，关键是短路容量要设得很大，比如10000MVA，相当于一个强系统。如果短路容量设置得和发电机容量一个量级，无穷大系统就不再“无穷大”，两侧都会摇摆，模型含义就变了。

第五步，放置故障模块和断路器。Three-Phase Fault模块支持编程控制，你可以直接在这个模块里设置故障开始时间和结束时间，在结束时间点就实现了故障切除。但如果你还想模拟故障后线路跳闸，就需要用Three-Phase Breaker配合Control信号来实现。

全部连好之后，在发电机模块输出端口上接一个转速测量通道。推荐直接使用发电机模块输出端的“wm”物理信号（机械转速），单位是pu，再乘以同步转速 ( \omega_0 = 2\pi f_n = 314.159 )rad/s 就能得到实际值，也可以用bus selector把转速信号引出来看。

我给一个具体的故障时序建议：0~1秒系统正常运行，让发电机达到稳定；1.0秒时在双回线其中一回的50%处发生B、C两相接地短路；1.1秒时切除故障；1.1秒后观察转速的摇摆过程直至仿真结束（比如5秒）。这样你能看到完整的“正常运行—故障扰动—故障切除—转子摇摆—趋于稳定”全过程。

4.2 第二版：自建状态方程模型（发电机组微分方程）

如果你选方程级路线，搭建过程就完全是另一种画风。打开一个新的Simulink模型，拖入一个MATLAB Function模块，在里面写转子运动方程和电磁功率计算逻辑。

我提供一个简化但物理正确的模型框架。经典二阶模型，状态变量为功角 ( \delta ) 和转速偏差 ( \Delta\omega )：

function [delta, omega, Pe] = swing_eqn(Pm, tf_type) % 单机无穷大系统二阶转子运动方程 % 输入: Pm 机械功率(pu), tf_type 故障类型(0-正常,1-两相接地,2-切除后) % 输出: delta 功角(rad), omega 机械转速(pu), Pe 电磁功率(pu) global Xd H omega0 E U delta0 omega0 % 根据故障类型确定等效电抗 switch tf_type case 0 % 正常: Xd1 + XT + XL/2 Xeq = 0.3 + 0.1 + 0.2; case 1 % 两相接地短路: 正序电抗 + 负序与零序并联 X2 = 0.2; X0 = 0.1; Xeq = 0.3 + 0.1 + 0.2 + (X2*X0/(X2+X0)); case 2 % 故障切除后单回运行 Xeq = 0.3 + 0.1 + 0.4; end % 电磁功率(暂态电抗后电势恒定模型) Pe = E * U / Xeq * sin(delta); % 转子运动方程 d_delta = omega0 * (omega - 1); % 功角导数 d_omega = (Pm - Pe) / (2*H); % 转速导数 end

当然，实际在MATLAB Function里直接写微分方程不方便，你需要配合积分器模块（Integrator）搭建一个反馈回路——把转速积分得到功角，功角反馈进函数算电磁功率，功率差再驱动求和模块，形成闭环。这里不再展开每个模块的连线，核心思想是把二阶微分方程降阶为两个一阶微分方程，用积分器串联实现。

这种模型的优势是，你可以在同一套框架下任意修改故障类型对应的电抗表达式，快速对比三相短路、两相接地短路、单相接地短路等不同故障类型的转速摇摆差异，特别适合用来做机理分析和论文里的参数扫描图。

4.3 故障时序与参数设置：暂态稳定仿真的“灵魂”

无论走哪条路线，故障时序设置都是整个仿真成败的关键。我强烈建议你建立一套“统一时序约定”，比如：

• ( t = 0 \sim 1.0 )s：系统空载或带载稳定运行。这段时间要留够，让发电机功率角、励磁电压等变量收敛到稳态。
• ( t = 1.0 )s：投入故障（两相接地短路）。
• ( t = 1.1 )s：切除故障（跳开故障线路或断开故障点）。
• ( t = 1.1 \sim 5.0 )s：观察转子摇摆过程。

故障持续时间 ( t_f ) 我在文章开头提过，它与极限切除时间 ( t_{cr} ) 的关系至关重要。你可以在保持其他条件不变的前提下，把 ( t_f ) 从0.1s逐步增加到0.3s，会看到一个有趣的临界现象：较短的故障时间下，转速经过几次衰减振荡后回到同步转速；当 ( t_f ) 超过某个阈值后，转速振荡不再衰减，甚至持续发散——这个阈值就是仿真得到的极限切除时间。

这里有一个容易被忽略的参数：仿真步长。Simulink默认的变步长求解器（ode23tb或ode15s）在电力电子类仿真中有一定优势，但暂态稳定仿真属于机电暂态，时间尺度在几百毫秒到几秒，频率成分相对较低。我个人经验是，用固定步长求解器、设置最大步长为0.1ms到1ms就足够了。如果步长取得太大（比如10ms），电磁暂态过程的分辨率不够，转速波形的高频细节会失真；如果小到微秒级，计算时间急剧增加，除非你关心的是故障瞬间的电磁暂态细节，否则没必要。

土办法是“步长敏感性测试”：先用1ms步长跑一遍，再用0.1ms跑一遍，如果转速曲线的包络线基本一致，说明1ms足够；如果两条曲线差异明显，继续减小步长直到曲线重合。

4.4 波形观测与数据导出

很多人仿真跑完就在Scope里看一眼波形，觉得差不多了。但做暂态稳定研究，光用Scope看是远远不够的。你要记录精确的转速数值、功角极值、振荡频率、衰减时间常数，这些都要靠数据导出做进一步分析。

在Simulink里，推荐用“To Workspace”模块把转速、功角、电磁功率等信号以数组形式导出到MATLAB工作区。然后再写几行脚本做后处理：

% 判断系统是否暂态稳定：检查转速偏差是否收敛 omega_pu = omega_signal; % 导出的转速信号(pu) omega_dev = abs(omega_pu - 1); [t_peak, idx_peak] = findpeaks(omega_dev); % 找振荡峰值 if t_peak(end) > t_peak(end-1) * 0.9 disp('系统暂态失稳'); else disp('系统暂态稳定'); end % 计算最大转速偏差 [max_dev, idx_max] = max(omega_dev); fprintf('最大转速偏差: %.4f pu\n', max_dev); fprintf('发生时间: %.4f s\n', t_out(idx_max));

这段脚本的价值在于把“看到波形能感觉到不稳定”变成“用指标判断是否稳定”。你可以把“最大转速偏差小于某个阈值（比如0.002pu）且振荡衰减到一定范围”作为稳定判据。学术论文里写“系统在故障切除后2秒内恢复同步转速”，就靠这个数据支撑。

5. 发电机转速变化结果解读：从波形看暂态稳定本质

5.1 典型转速响应曲线：加速、减速、再同步的过程

当你仿真跑完，把转速波形调出来，你会看到一条非常典型的曲线。正常运行时转速是一条水平直线（等于同步转速，pu值下就是1.0）。故障发生的瞬间，转速开始上升——这个上升不是瞬时的，而是有一个斜率逐渐增大的过程，因为转子有惯性，加速需要时间。故障期间转速一直上升，直到故障切除时刻。

故障切除后，转速并不会立刻掉头下降，而是会继续冲高一段距离。原因很简单：故障切除的瞬间，转子转速已经高于同步转速，转子里的动能不可能瞬间消失。功角继续向前冲，直到减速功率做的功抵消了多余动能，转速才到达峰值，开始回落。转速峰值的位置出现了一个明显的“拐点”。

然后转子进入减速阶段，转速下降。但下降过头之后，转速又会低于同步转速——因为功角冲过头，电磁功率大于机械功率，转子在减速区间被“拉过头”。于是转速在同步转速附近来回摆动，形成一条衰减振荡曲线。振荡的包络线逐渐收窄，最终转速回到同步转速，系统恢复稳定。

这就是单机无穷大系统暂态稳定的标准叙事。整个过程中，转速变化曲线的形状由三大因素决定：故障期间加速面积的大小（取决于故障类型和时长）、故障切除后的减速能力（取决于系统电抗和电压水平）、以及转子惯性（取决于H常数）。你在解读波形时，脑子里一定要有这三条线索，它们会在后续的每个参数分析场景中出现。

5.2 判断系统是否暂态稳定的三条标准

看波形不能靠“感觉”。我总结了三条可量化的稳定判据，你可以直接用MATLAB脚本计算：

第一，转速偏差是否收敛。暂态稳定系统的转速偏差应该随时间衰减，最终趋近于零。具体判据可以是：仿真结束时转速偏差的绝对值小于0.001pu，且振荡包络线单调递减。如果转速偏差越振荡越大，或者最后稳定在某个非同步转速上，说明系统失去稳定。

第二，最大功角是否不超过失稳门槛。单机无穷大系统有一个典型的失稳模式：功角在向 ( \delta > 90^\circ ) 甚至 ( \delta > 180^\circ ) 方向持续摆动，就基本宣告失步了。你可以记录每个振荡周期的功角峰值，如果峰值呈现增长趋势，且超过 ( 180^\circ ) 之后转速不再按同步规律来回穿越，系统失稳。

第三，振荡是否具有等幅或增幅特征。稳定系统的振荡是减幅振荡，相邻两个周期的振荡峰值之比小于1。你可以用信号处理的方法，从导出的转速数据中提取每个周期峰值，计算衰减率。如果衰减率为正（峰值不断增长），不管当前转速偏差有多大，系统终究要失稳。

这三个判据各有侧重，建议组合使用。单纯看转速偏差可能在某些特殊的临界切除时间下误判——比如系统处于临界失稳状态时，转速波形可能先衰减几拍，后续再重新发散。这时候只看最终值不够，必须检查振荡趋势的方向。

5.3 故障持续时间与临界切除角的关系

聊完稳定判据，再来做一个特别有价值的仿真小实验：改变故障持续时间，扫描 ( t_f ) 从0.05s到0.30s，每隔0.02s跑一次仿真，记录最大转速偏差和最终是否稳定。你会发现一个非常清晰的规律：

故障时间短时，比如0.05~0.12s，转速最大偏差随故障时间近似线性增长，但系统都能恢复稳定，最终转速回到同步转速。故障时间到了某个临界值附近，比如0.16s左右，系统反应变得敏感——故障时间微小增加，最大转速偏差却大幅跳升。超过这个临界值，转速曲线彻底发散，系统失稳。

这个临界值在物理上对应于“极限切除时间”，在功角平面上就对应于“极限切除角”。你可以把每次仿真的最大功角和故障切除功角记录下来，画在等面积定则的 ( P-\delta ) 图上，会直观看到加速面积和减速面积的关系：切除角小于极限值时减速面积大于加速面积，多余动能可以被消耗掉；切除角超过极限值，即使动用了整个减速区间也不足以消耗故障期间积累的动能，转子带着剩余动能冲过失稳点。

这个实验强烈推荐做。它不需要修改模型结构，只需要用脚本循环改变故障模块中的切除时间参数，自动批量跑仿真并记录结果。一张“故障切除时间—最大转速偏差”的关系曲线图，几乎可以当成整个暂态稳定研究的核心图表。

6. 常见问题与排查技巧实录

我把自己做这个仿真实操时遇到的问题，以及周围同学常踩的坑，整理成了一份速查表。

这里最想展开说的一点，是关于“仿真发散”这件事。Simulink的仿真发散往往不是模型数学上无解，而是数值积分算法在某个突变点（比如故障投入、切除的瞬间）遇到了不连续跳变，步长适应机制反复尝试导致步长缩到极小，或者算法在刚性方程上失效。解决办法有三个层次：

第一个层次是换求解器。Simulink默认的ode45在光滑问题上很好用，但电力系统这种强非线性+频繁开关切换的场景里，ode15s、ode23t这类刚性求解器明显更稳。第二个层次是切换故障和切除瞬间的步长设置：Simulink允许在模型中插入“Solver Configuration”或通过“Zero-Crossing Detection”选项来控制事件检测精度。第三个层次是给功率信号加上“Slew Rate Limiter”之类的限速模块，人为限制信号的跳变率，减轻数值冲击。当然，加限速模块是不得已的下策，它会稍微改变真实的物理动态，能不用就不用。

再说一个独家小技巧。如果转速波形在故障切除后出现数值振荡，但你又无法确定是物理还是数值问题，用“步长减半对比法”验证：把步长从1ms改成0.5ms，如果振荡幅度基本不变，说明是物理现象；如果振荡幅度明显下降或消失，说明是数值引入的伪振荡，需要调整求解器设置。这个方法花不了几分钟，但能帮你避免很多误判。

7. 扩展玩法：从单机到多机，从短路到重合闸

7.1 加入自动重合闸时序

两相接地短路仿真做完，我建议你再往前走一步：把故障切除改成“故障切除—无流时间—自动重合闸”的完整时序。实际电网输电线路绝大多数都配置自动重合闸，因为瞬时性故障（特别是单相接地）占了绝大多数，重合闸能显著提高供电可靠性。

在Simulink里实现这个时序很简单：用两个断路器模块叠加控制信号。第一个断路器在故障后0.1s跳开，隔离故障；等待0.5s的无流时间（让故障点电弧去游离）；第二个断路器（或同一个断路器再合闸）在跳开0.5s后重合。如果是瞬时性故障，重合成功，系统恢复；如果是永久性故障，重合后再次跳开，系统可能失稳。

你会在转速波形上看到比单次切除更复杂的摇摆过程：第一次跳闸引起一波振荡，重合闸瞬间又叠加一次冲击。有时候单次切除系统是稳定的，但重合闸后的二次冲击反而导致失稳。这个现象叫“重合闸冲击下的暂态稳定问题”，在工程实际中非常重要。做个自动重合闸仿真，你的单机无穷大模型就从一个教学示例变成了一个小型工程研究平台。

7.2 励磁系统与调速器的协同影响

我前面建议用的发电机模型是经典二阶模型——暂态电抗 ( x'_d ) 后电势恒定。这种模型忽略了励磁系统的动态调节作用。现实中发电机都装自动电压调节器（AVR），故障时电压跌落会触发强励，提高暂态电势，从而增大电磁功率上限，改善暂态稳定。

如果你想更贴近实际，可以在Simulink里给发电机加上励磁系统模块（比如IEEE Type I Excitation System）和调速器模块（比如Hydraulic Turbine and Governor）。你会发现两个有趣的现象：加了励磁系统后，同样的故障情况下，转速振荡幅度变小，系统极限切除时间变长——稳定裕度提高；而调速器虽然主要管频率恢复，但也在故障过程中缓慢调节机械功率，影响摇摆曲线的末期衰减特性。

这里有个建模上的提醒：加了励磁和调速器后，模型时间常数跨度变大（从毫秒级到秒级），数值求解的刚性增强。你要特别注意求解器的选择，建议用ode23t或ode15s，并适当放宽仿真步长，否则计算时间会明显增长。

7.3 向多机系统扩展的思路

单机无穷大系统的一切分析都建立在“无穷大系统纹丝不动”的假设上。但真实电网是很多台发电机互相摇摆、互相影响的。多机系统的暂态稳定分析复杂度大增，因为每一台发电机都有自己的转子运动方程，彼此通过网络方程耦合，功角之间的相对运动取代了绝对功角成为观测量。

在Simulink里做多机扩展，最简洁的方式是在现有单机模型基础上，把无穷大系统模块替换成另一台同步发电机模型（带自己的原动机、励磁），两台发电机之间用线路连接。然后你可以在某台发电机附近设置两相接地短路故障，观察两台发电机的转速偏差曲线——你会发现它们不是一起摇摆，而是相对摇摆，功角差来回振荡。两个机群的相对功角才是判断多机系统暂态稳定的关键指标。

这个扩展方向做成之后，研究问题的深度和广度会完全不同，因为两机系统是所有多机稳定性的最小原型，两机搞懂了，多机系统只是规模和网络拓扑的进步。你现在的单机模型完全可以作为搭建多机模型的基础素材——变压器、线路、故障模块直接复用，只需要新增一台发电机模块及其控制信号就行。

我个人在实际操作中最大的体会是：单机无穷大系统之所以经典，不是因为它简单，而是因为它把“什么因素影响暂态稳定”这个问题拆解得足够干净。做完两相接地短路仿真，如果你能不看示波器、只靠等面积定则估算出大概的临界切除时间，再看仿真波形验证自己的判断，那你对这个领域的理解已经上了一个台阶。

还有个小建议：做这类研究时，养成“每改一个参数就单独导出数据并做记录”的好习惯。不要只在Scope上看一眼觉得像了就过——把转速峰值、振荡周期、衰减率这些量化指标存下来，建立自己的仿真数据表。后面写报告、写论文、做对比分析，这套数据就是你最可靠的素材。