企业官网建设流程全解析

做电力系统仿真的人，几乎都会撞上同一个题目：同步发电机突然三相短路。我第一次拿到这个题时，心里还嘀咕：“突然”两个字是不是多余的，难道还有“逐渐短路”吗？后来真正把Simulink模型搭起来、把波形拉出来之后才明白，“突然”强调的是短路发生前电机处于稳态，而短路瞬间定子电流几乎是在几个毫秒内跳变上去的，这个跳变光靠欧姆定律解释不通，必须把磁链守恒、转子绕组感应电流和一系列时间常数全部拉进来，才能把过程说清楚。这篇文章就把“同步发电机突然三相短路暂态过程研究”的完整路线写下来，从物理原理、模型搭建、波形解读到调试踩坑，尽量按我自己做过的流程来还原。适合电气工程方向的学生、做发变组保护整定的工程师，以及所有想用Simulink把课本上的公式变成眼见为实波形的人。

1. 突然三相短路：为什么它是发电机绕组的一场“生存危机”

1.1 电抗突减，电流自然暴增

同步发电机稳态运行时，机端突然发生三相短路，外接负载阻抗被短路点几乎“短路”掉，发电机定子绕组面对的等效阻抗瞬间从同步电抗变成更小的暂态电抗、次暂态电抗。可以这样类比：自来水管前面装了一个大阻力阀门，水流本来不大；结果阀门突然被拆掉，水源压力基本没变，管道里的流速立刻冲上去。定子电流的大小由电动势和电抗决定，电抗从同步电抗 (X_d) 骤降到次暂态电抗 (X_d'')，而 (X_d'') 往往只有 (X_d) 的十分之一甚至更低，短路电流必然在很短时间内冲到很高的数值。

很多教材上用“冲击系数”来描述最恶劣情况下的电流峰值，冲击电流可能达到周期分量起始幅值的1.8倍到2倍。为什么能这么高？因为定子电流不仅含有基频周期分量，还叠加了一个按指数衰减的直流分量，两者在同一时刻方向相同，瞬时值自然被顶得很高。这种冲击电流对发电机绕组端部产生巨大的电动力，对绕组绝缘和机械结构都是严重考验，所以电力系统设计、保护整定、断路器开断能力校核都必须研究这个过程。

1.2 磁链守恒背后的“补偿电流”

真正让突然短路区别于稳态短路的核心原因是绕组的磁链不能突变。无论是定子绕组还是转子绕组，只要它是闭合回路，在短路瞬间都要保持自身磁链不变。外接阻抗突然减少，定子电流趋势上是跳变的，但磁链的连续性不允许定子电流按照简单的欧姆定律立刻跳变，于是绕组内部会感应出额外的自由电流，用磁场的变化去“抵消”外部短路带来的冲击。

这个思路也是后续分析所有暂态分量的钥匙。定子绕组为了保持磁链不变，会产生直流分量和非周期性质的自由电流；转子上的励磁绕组和阻尼绕组为了保持自身磁链不变，也会感应出直流电流，这些转子电流反过来又在定子侧感应出基频交流分量。转子绕组里的电流不能一直存在，因为绕组有电阻，自由电流会按各自的时间常数衰减，衰减的快慢直接决定了定子短路电流周期分量的衰减速度。

1.3 时间常数决定了故障发展的“三幕剧”

从工程角度观察，突然三相短路后的电流并不是瞬间稳定下来，而是分三个阶段变化：

这三个阶段对应三个电抗：(X_d'')、(X_d')、(X_d)。短路瞬间首先由 (X_d'') 决定电流幅值，因为阻尼绕组和励磁绕组都在起作用，等效电抗最小，电流最大；阻尼绕组作用衰减掉以后，只剩下励磁绕组在维持磁链，等效电抗变成 (X_d')；等到励磁绕组的自由电流也衰减干净，才会过渡到 (X_d) 决定的稳态短路电流。很多同学仿真出来的波形“看不出次暂态”，本质上是因为时间常数太小、仿真步长太大或者观测窗口不对，把这个过程漏掉了。

1.4 仿真研究解决了什么问题

用Simulink做这个题目，一方面是把复杂的电磁暂态过程可视化，让“次暂态”“暂态”“直流衰减”这些抽象概念直接在波形上显示出来；另一方面，工程中需要知道短路电流的起始峰值、周期分量衰减到多少、什么时候达到稳态，这些数据可以直接用于断路器开断能力校验、发电机绕组受力分析和继电保护定值配合。单纯解析计算很繁琐，尤其是带阻尼绕组、饱和效应、多机互联的情况下，手算几乎不可能，仿真就成了最实用的工具。

2. 建模之前，先在纸上把暂态电流的变化规律走一遍

2.1 定子电流其实由两个分量“拼”出来

很多初学者一拿到Simulink就开始拖模块，结果波形出来后完全看不懂。我建议先别急着头文件，拿纸笔把电流组成捋清楚。

定子任意一相短路电流的瞬时值，大体由两部分组成：

• 基频周期分量：它本身也在随时间变化，起始幅值由 (X_d'') 决定，然后按照次暂态、暂态两个时间常数衰减到稳态幅值；
• 非周期分量：这是短路瞬间为了满足定子电流不能突变而产生的直流偏移量，按定子非周期时间常数 (T_a) 衰减，典型值0.1s~0.3s。

为什么会有直流分量？看一个简单类比：交流电流在正常情况下是正弦曲线，短路瞬间如果按新的电抗计算，交流分量的起始值与故障前瞬间的电流值往往不相等，而电感电流又不能突变，中间出现了“缺口”。为了让电流在短路瞬间保持连续，就只能叠加一个方向合适、大小等于缺口的直流分量，让整体电流在故障瞬间“接得上”，然后这个直流分量再慢慢衰减掉。

所以仿真中短路发生时刻不同，直流分量的大小完全不同。如果短路瞬间恰逢某相电流过零，那么这个缺口可能很小；如果恰逢电流峰值，缺口可能很大。三相短路时三相的直流偏移也不一样，这就是为什么我们用示波器看三相电流时，会发现某一相或两相明显偏向时间轴一侧。

2.2 转子绕组如何“拖住”定子电流

定子周期分量的衰减，本质上是转子侧自由电流衰减的结果。短路瞬间，定子电枢反应产生的磁通突然变化，转子励磁绕组和阻尼绕组都会感应出直流电流，去维持扰动前的磁链。这些直流电流建立的磁通会在定子侧感应出基频电动势，相当于让定子短路电流暂时保持一个很高的水平。励磁绕组电阻小、电感大，时间常数长，所以它维持的“暂态分量”衰减慢；阻尼绕组电阻相对大一点、时间常数短，所以它维持的“次暂态分量”衰减快。

这里也解释了为什么没有阻尼绕组的发电机，短路电流起始阶段看不到明显的次暂态过程，电流直接从暂态幅值开始衰减。Simulink中的同步电机模型如果设置了阻尼绕组参数，次暂态阶段就能体现出来；如果参数不设置或设成零，波形就只有暂态和稳态两部分。

2.3 用dq0坐标看更清爽

三相静止坐标系里各种交变分量混在一起，很难直接读出时间常数。实际分析时都会通过Park变换把定子量变到转子旋转坐标系，也就是dq0坐标。在这个坐标系里，定子电流的基频周期分量变成了近似直流量，非周期分量反而变成了交流量，而且会分解出直轴分量 (i_d) 和交轴分量 (i_q)。

对于理想的三相突然短路，直轴电流 (i_d) 的变化最能说明问题。短路前 (i_d) 是某个稳态值；短路的瞬间跳到一个很高的次暂态幅值，随后按 (T_d'') 快速衰减到暂态幅值，再按 (T_d') 衰减到稳态短路值。把这条曲线拉出来，用指数拟合就能实测出时间常数，这是仿真研究最有价值的地方，也是我在后面章节会详细展开的操作。

2.4 简化表达式，心里有数

工程上常把直轴电流的变化写成三段式：

[ i_d(t) = I_{\infty} + \left(I_d' - I_{\infty}\right)e^{-t/T_d'} + \left(I_d'' - I_d'\right)e^{-t/T_d''} ]

其中 (I_{\infty}) 是对应同步电抗的稳态短路电流，(I_d') 是对应暂态电抗的起始电流，(I_d'') 是对应次暂态电抗的起始电流。这个表达式不是严格从电机五阶方程推出来的，但用来理解仿真波形的包络线足够。实际波形还会叠加直流分量和由于转子不对称引起的2倍频分量，这些在仿真波形里也能观察到，但并不影响我们抓住主趋势。

3. Simulink中的建模步骤：从同步电机到故障投切

3.1 先确定模型拓扑

我用的仿真拓扑很基础：一台同步发电机经过一段线路和负载连接，机端并联一个三相短路故障模块。短路时刻由故障模块控制，短路前发电机带负载稳态运行，短路后三相同时接地。不需要接无穷大系统，因为我们要观察的是发电机本身对短路点提供的暂态电流；当然，如果后面要研究系统侧馈入短路电流，可以再加入电源和变压器，不过初学者先不要急着扩大模型。

在Simulink库中，我建议搜索关键字定位模块，不同版本库路径略有差异，大致在：

• 同步电机：Simscape > Electrical > Specialized Power Systems > Machines > Synchronous Machine
• 三相故障：Simscape > Electrical > Specialized Power Systems > Elements > Three-Phase Fault
• 电压电流测量：Simscape > Electrical > Specialized Power Systems > Measurements > Three-Phase V-I Measurement
• 电源、负载、线路：Simscape > Electrical > Specialized Power Systems下相应元件
• Powergui模块：Simscape > Electrical > Specialized Power Systems > Powergui

Powergui这个模块很多人容易漏掉，凡是仿真模型里含有同步电机、变压器或电力电子开关的，都必须放一个Powergui。它是电力系统仿真的求解器管理器，也是稳态初始化工具。

3.2 同步电机模型怎么选

Synchronous Machine模块有几种变体，比如SI Fundamental、pu Fundamental、pu Standard等。我的选择是pu Fundamental，因为参数直接使用标幺值，和电机铭牌、暂态电抗参数手册对应起来非常方便。Standard模型通常更复杂，考虑更多细节，但对突然三相短路这个题目来说，Fundamental模型已经足够清晰，跑起来也更快。

如果是做更深入的研究，比如分析谐波、饱和、多质量轴系扭振，再考虑换更细的模型。先用简单模型把物理过程跑通，再逐步加复杂度，这个思路适用于绝大多数仿真课题。

3.3 发电机参数怎么填

我建议先用一组典型参数做参考，熟悉之后再换成你自己研究的实际电机参数。以下是一组100MVA级汽轮发电机的标幺值示例：

注意不同版本模块可能要求的是开路时间常数 (T_{d0}')、(T_{d0}'')，填写前先看对话框里的说明，用哪种就填哪种。电机参数用pu模型时，基准值由模块内部按额定功率、额定电压自动生成，不需要额外换算成国际单位，这也是我推荐pu模型的原因。

3.4 原动机输入和励磁怎么设置

同步电机模块的输入端口通常有两个：机械功率 (P_m) 和励磁电压 (V_f)。最简单的做法是用两个常数模块分别给0.8 pu和1.0 pu。0.8 pu代表原动机输入机械功率大约是额定值的80%，对应带负载运行；励磁电压给1.0 pu后，让发电机在初始化时自动整定稳态响应。

如果你没有连接调速器和励磁调节器，短路后发电机转速会缓慢上升，因为短路时机端电压接近零、电磁功率很小，而原动机还在持续注入机械功率。这个现象是正常的，正好也能帮助我们理解短路对发电机功角稳定性的影响。如果只关心电磁暂态，这个简单输入方式足够了。

3.5 故障模块和测量环节

Three-Phase Fault模块可以设置故障相别和动作时间。我这里做三相短路，就选中A、B、C三相都接地；故障时刻设为0.5s，结束时刻设为2s；故障电阻我习惯设置成0.001Ω而不是0Ω，原因后面会在踩坑部分详细说。

测量方面，我建议至少加两组：

1. 在发电机机端串联一个Three-Phase V-I Measurement，把三相电压、电流引到Scope里；
2. 从同步电机模块的测量输出端口引出一根线，接Bus Selector，挑选你关心的信号，比如电磁转矩 (T_e)、转速 (w_m)、励磁电流 (i_f)、直轴电流 (i_d)、交轴电流 (i_q)。

这样既能看定子瞬时电流，又能从dq坐标角度观察暂态衰减过程。

3.6 Powergui初始化和求解器参数

双击Powergui，选择“Machines and Load Flow Initialization”。在这个面板里选择同步电机，并点击执行初始化，软件会根据当前拓扑和给定功率、电压自动计算出初始励磁、转子角和各绕组初始电流。做完这一步，仿真才能从稳态直接开始，否则一开机就是剧烈的启动过渡，和短路的暂态过程混在一起，没法看。

求解器方面，我习惯在Simulink仿真参数里选变步长隐式求解器，例如ode23t或ode15s，最大步长设为1e-4s。因为次暂态时间常数只有几十毫秒，最大步长太大，仿真直接就“跳过”了次暂态阶段，波形只剩暂态衰减，看起来就像纯电抗模型算出来的结果。仿真总时长可以先设2s，够观察到暂态时间常数衰减到稳态的过程。

4. 仿真结果怎么读：包络线、直流偏移与时间常数

4.1 定子三相电流波形：第一眼看什么

仿真跑通后，先看Scope里的三相定子电流。短路前是三相对称正弦波形，幅值比较小；短路瞬间，电流猛然增大，而且三相波形明显不再对称于横轴。这里的不对称不是指幅值不同，而是某一相或两相的波形整体向上或向下偏移，这就是前面说的非周期直流分量在起作用。

随着时间推移，直流偏移逐渐衰减到零，三相电流重新变成对称的正弦波，但幅值仍然高于短路前，然后继续缓慢下降，最后稳定在由同步电抗决定的稳态短路电流值。整个过程看起来就像一个很宽大的包络线慢慢“收窄”。

要提取冲击电流，就把短路后第一个周波附近的波形放大，找最高峰值。以50Hz系统为例，冲击电流通常出现在短路后半个周波，也就是10ms左右。拿仿真得到的瞬时电流峰值除以周期分量起始幅值，能得到冲击系数，正常情况应该在1.8附近。如果偏得太多，先检查是不是短路时刻恰好打在相角峰值上，或者直流分量衰减太快。

4.2 从dq轴分量看衰减的三段性

三相电流混合了基频和直流，读起来费劲。这时候一定要看 (i_d) 的波形。理想情况下，短路前 (i_d) 近似恒定；短路瞬间跳到一个很高的值，然后出现两个明显速度不同的衰减段：

• 前50ms内， (i_d) 快速下降，这是次暂态分量在衰减；
• 之后下降速度明显变慢，这是暂态分量在衰减；
• 最终平缓下来，进入稳态短路电流。

用我在3.3节给的参数估算，(i_d) 起始幅值接近 (1.0 / X_d'' = 5) pu，经过次暂态衰减后回到 (1.0 / X_d' = 3.3) pu，再经过暂态衰减后回到 (1.0 / X_d = 0.56) pu。这个“5 → 3.3 → 0.56”的数字，在仿真波形里是可以清楚看出来的。如果仿真步长太大，第一段断崖式下降会直接变成一条直线，看不到中间的3.3 pu平台。

同时可以看励磁电流 (i_f)。短路瞬间 (i_f) 会有明显突增，随后逐渐回落，还可能叠加交流振荡分量。这是因为定子非周期分量产生的静止磁场相对转子是旋转的，会在转子绕组里感应出交流量。看到这些现象不要慌，这是电磁暂态的固有特征。

4.3 电磁转矩和转速的“异常”反而有价值

短路瞬间电磁转矩会出现很大的冲击振荡，随后因为机端电压很低，电磁功率几乎为零，而原动机机械功率还在输入，发电机转子开始加速。如果你只跑2s，转速变化可能不大，但趋势明确；如果把仿真时间延长，转速上升会更明显。这个现象直接引出了暂态稳定问题：短路故障如果不能及时切除，发电机持续加速，功角不断增大，最终可能失步。

所以这个仿真不只是为了看电流大小，它同时展示了电磁暂态与机电暂态的耦合。短路瞬间先看到电气的快速变化，随后才看到转子转速和功角慢慢变化，这两种时间尺度差异在同一个模型里体现得非常直观。

4.4 用Matlab拟合实测时间常数

仿真不只是看个热闹，还可以把数据导出来做定量验证。通过To Workspace把 (i_d) 和仿真时间导出到工作区，然后用Matlab的lsqcurvefit拟合指数曲线。下面是一个示例脚本：

% 假设 Simulink 导出的变量为 t_out（时间）和 id_out（直轴电流） idx = t_out >= 0.5; % 取故障发生后的数据 t = t_out(idx) - 0.5; % 时间归零 y = id_out(idx); % 拟合模型: i_d(t) = c0 + c1 * exp(-t/T1) + c2 * exp(-t/T2) model = @(p, t) p(1) + p(2) * exp(-t / p(3)) + p(4) * exp(-t / p(5)); p0 = [0.56, 1.6, 0.5, 2.0, 0.03]; % 初值需根据波形大致调整 options = optimoptions('lsqcurvefit', 'Display', 'off'); pfit = lsqcurvefit(model, p0, t, y, [], [], options); fprintf('稳态值 = %.3f pu\n', pfit(1)); fprintf('暂态时间常数 Td1 = %.3f s\n', pfit(3)); fprintf('次暂态时间常数 Td2 = %.3f s\n', pfit(5));

拟合前最好先截掉第一个周波内的数据，或者只拟合直流分量已基本衰减到零之后的包络线，否则非周期分量会干扰结果。实际做下来，只要模型参数和求解器设置合理，拟合出的 (T_d')、(T_d'') 能和设置值对应上，误差来源主要在于非周期分量和高频分量的残留。

5. 调试中踩过的坑：初始化、故障电阻和求解器这三件事

5.1 初始条件没设置好，看到的是“假短路”

最开始我做这个仿真时，直接拖好模型就点运行，结果短路前发电机电流就在疯狂振荡，转速从零开始爬升，励磁电流也是从头建立。等到0.5s短路发生，前面根本还没进入稳态，得到的结果完全无法和理论对应。

问题就出在没有做稳态初始化。解决办法是在Powergui里执行Machines and Load Flow Initialization。这一步会基于你给的机械功率、励磁电压和负载，反算出电机应该有的初始状态。做完初始化再运行，短路前波形就是一条平整的正弦曲线，短路后的波形才能真正反映暂态过程。

如果模型加了变压器、线路、调整了负载，每次改动参数都建议重新做一次初始化。我见过不少同学改了负载电阻后没重新初始化，仿真结果越来越离谱，还以为是模型搭错了。

5.2 故障电阻不是越小越好

理论上三相金属性短路就是故障点阻抗为零，所以很多人把Three-Phase Fault的故障电阻填成0。但Simulink的电力电子/机电仿真模型在零过渡电阻时容易出现数值刚性问题，表现为仿真变慢、报错，或者波形上出现密集的毛刺。

我建议填0.001Ω，这个值和金属性短路在工程上几乎没有区别，但数值稳定性会好很多。另外，Three-Phase Fault模块里如果选了带外部控制断口，还需要检查外部控制信号是否在正确时间切换；如果用内部定时，故障时间矩阵要写成两列，每一行代表一个开关的动作区间，例如[0.5 2]表示0.5s发生短路并一直保持到2s。

5.3 最大步长和采样时间决定你能看到什么

次暂态时间常数可能只有30ms，如果最大步长设为0.01s，每个周期只有几个采样点，(i_d) 从5pu降到3.3pu的过程可能就被一个大步长“一步跨过去”，波形看起来直接从5跳到了3.3，中间的指数过渡完全丢失。更严重点说，冲击电流峰值根本采不到。

所以最大步长至少要小于最快时间常数的五分之一到十分之一，我通常设1e-4s。Scope或者To Workspace的采样时间也要同步设置，如果是连续示波器，通过Decimation控制采样率；如果用离散采样模块，建议采样时间设1e-4s，否则可能出现“波形稀疏、峰值丢失”的现象。

5.4 常见问题和排查方法

5.5 仿真时间和观测窗口怎么折中

观察次暂态需要很小的步长，观察暂态和稳态需要很长的时间窗，这二者对计算量是矛盾的。我的习惯是把仿真分为两个阶段看：

• 先跑0到0.1s，重点观察短路瞬间的冲击电流和次暂态衰减；
• 再跑完整2s，重点考察暂态时间常数和转速变化。

如果模型很大、仿真很慢，可以用两个模型分别做短时间和长时间仿真，不必在同一个模型里追求所有细节都兼顾。这种做法听着简单，实际能省掉大量调参时间。

6. 从单机突然短路说开去：还能再接再厉做哪些扩展

6.1 从三相短路到不对称故障

突然三相短路是分析基础，但实际电网中单相接地短路、两相短路更常见。Simulink里改起来很简单，Three-Phase Fault模块里选不同的故障相别即可。不过分析波形时会发现，不对称故障除了出现电流增大，还伴随着负序和零序分量，需要借助对称分量法来理解。单机模型的短路电流中会出现负序磁场，在转子里感应出100Hz分量，所以波形振荡特征和三相短路很不一样。

如果要做不对称接地故障，还要注意故障模块的接地方式、变压器中性点接地方式等细节，这些都会影响零序通路。建议先把三相短路彻底搞明白，再玩不对称故障，否则很容易被波形搞懵。

6.2 把结果用在校验和整定上

发电机端的突然三相短路电流是继电保护整定计算的重要输入。保护用电流互感器看到的电流中包含衰减的周期分量和非周期分量，这对保护的启动、速断、延时配合都会产生影响。把Simulink输出的暂态电流数据导入保护逻辑模型，可以验证差动保护、过流保护在暂态过程中的动作行为。

断路器开断能力校核也需要关注短路电流中的直流分量，因为直流分量不会过零，断路器开断电流时如果直流分量过大，熄弧会更加困难。仿真中得到的冲击电流、衰减时间常数，可以作为动稳定、热稳定校核的输入。把一个“课程仿真题”变成工程计算工具，其实只差这一步。

6.3 接入多机系统和暂态稳定分析

单机模型的下一步自然是多机系统。两台发电机经线路并联，某台机端短路时，不但故障机本身经历电磁暂态，另一台机组也会受到冲击，功角在两机之间摆动。此时需要关注的物理过程从“定子电流冲击”扩展到“暂态稳定”，也就是短路故障能否及时切除，发电机组是否保持同步。

在多机系统中，短路时刻的选择、故障持续时间、切除时间对结果影响非常大。同一个故障，0.2s切除和0.4s切除，系统可能一个是稳定的、一个是失步的。沿着这个方向，可以做极限切除时间CCT的计算，这也是电力系统暂态稳定课程里特别经典的课题。Simulink模型不用推倒重来，只要把原来的单机增加成两机或多机，加上输电线路和断路器，故障切除改用时间控制的算法即可。

6.4 给做课程设计和毕业设计的人一点建议

如果这个题目是你的课程设计或毕业设计，我的建议是不要贪大。先把单机带载模型跑稳，把三相短路波形、时间常数、冲击电流这几个核心结果整理出来，再考虑扩展。汇报时最打动人的不是你建了多复杂的系统，而是你能不能清楚地解释：为什么短路瞬间电流那么大，为什么衰减分几段，为什么短路时刻不同直流分量不同。这些用Simulink波形配合理论公式讲清楚，已经是一份很扎实的工作。

我在实际带着别人跑这个模型时，最大的体会是：真正重要的不是把模块拖出来，而是知道每个参数、每个模块、每个初始条件在物理上对应什么。初始化和最大步长这两个细节，几乎决定了仿真结果能不能看，但很多教程一句话带过，导致新手反复踩坑。建议你跑通一次之后，故意把自己的参数改到“错误”状态，比如不初始化、步长调大、故障电阻设零，再对照正确配置看波形差异。这几组对比做下来，对暂态过程的理解会比只看一遍理论深刻得多。