1. 这项研究到底在研究什么突然三相短路为什么值得单独拎出来做仿真做电力系统方向的人几乎都绕不开同步发电机的突然三相短路暂态过程。我最早接触这个课题的时候也觉得不就是短路嘛稳态短路算个电流不就完事了。真把Simulink模型搭起来、波形拉出来看之后才发现突然短路和稳态短路完全是两码事这里面的暂态分量、衰减时间常数、冲击电流峰值每一环都藏着大量工程细节。先说清楚一个容易被新手忽略的概念所谓“突然三相短路”指的是发电机在正常运行状态下机端突然发生三相对称短路。这个时候定子绕组里会出现远远大于稳态短路电流的冲击电流峰值可能达到额定电流的十几倍甚至更高。这个冲击电流会对绕组产生巨大的电动力同时电磁转矩剧烈波动对轴系造成冲击。因此这项研究不是为了看“短路以后最终电流是多少”而是要看“短路的头几个周期内电流怎么变、峰值多大、按什么规律衰减”。从物理本质上讲突然三相短路的暂态过程由两个基本定律支配。一个是电磁惯性定律也就是磁链不能突变另一个是超导回路磁链守恒原理。需要注意这里说的“超导回路”是等效概念因为在分析暂态过程的初始瞬间我们可以忽略所有绕组电阻把各个绕组看成理想的超导闭合回路。这样一来短路前定子、转子各绕组的磁链会保持守恒迫使短路后在各绕组中感应出附加的电流分量以此来维持磁链不变。很多资料直接给出结论短路电流中包含基波交流分量、非周期分量以及二次谐波分量。基波分量由转子励磁磁势和电枢反应共同决定非周期分量是为了维持短路瞬间定子磁链初始值而产生的直流分量二次谐波分量则是转子不对称性直轴和交轴磁阻不同造成的。理解这三个分量的来源直接关系到后面仿真时怎么看波形、怎么对应理论公式。这个课题的实验条件在物理实验室比较难做一来发电机容量不小二来真的去制造一次突然短路对设备冲击很大风险也高。所以Simulink仿真是非常合适的研究手段既能把电磁暂态过程还原出来又能随意改变故障合闸角、切除时间、负载条件等变量很划算。2. Simulink仿真建模从元件选型到主电路拓扑2.1 选哪款发电机模型同步电机模块的两种用法Simulink的Simscape Electrical早期叫SimPowerSystems里提供了同步电机模型一般用“Synchronous Machine”模块。标准库里有两种典型配置一种是基本型输入端口有Pm机械功率和Ef励磁电压输出是定子三相电流、转速、电磁转矩等测量信号。这种模型适合我们这种把励磁系统简化处理的研究。另一种是带励磁系统细节的需要外接励磁调节器。我做这个研究时直接用的基本型因为重点是短路暂态励磁调节器在暂态初期还没来得及动作对第一个峰值的贡献可忽略但如果你想研究暂态后期电压恢复那就另说。选型上还有一个关键点要把电感的饱和效应设置为考虑或者不考虑。同步电机模块里有一个参数“Rotor type”一般选“Salient-pole”凸极还是“Round”隐极这决定了直轴和交轴电抗是否相同。汽轮发电机一般是隐极直轴和交轴电抗相近水轮发电机是凸极交轴电抗明显大于直轴。做突然三相短路研究时为了简化理论分析常假设发电机为隐极机也就是xd等于xq。如果你用的参数表里xd和xq不相等在理论对比的时候就要注意公式会带上二次谐波分量。2.2 主电路拓扑用无穷大电网还是机端直接短路做同步发电机突然短路仿真通常有两种接线方案。第一种是发电机单独带一个负载通过一个三相断路器直接接到地或者短路点。这个方案贴近“发电机机端短路”的研究场景。第二种是发电机经变压器和双回线接入无穷大系统然后在机端或者线路某点设置短路故障。这更贴近实际电力系统但在分析暂态过程时容易把电网侧的影响混进来不容易聚焦发电机本体特性。我做这个课题时采用的方案是同步发电机出口接一个三相断路器断路器另一侧直接三相对地短路同时保留一个可选的负载支路模拟短路前带载状态。断路器模块用“Three-Phase Breaker”Simulink里需要注意的是断路器有一个“Breaker resistance Ron”和“Snubber circuits”参数。很多人忽略了Snubber缓冲电路的设置导致仿真中出现非物理的电压尖峰甚至数值振荡。我的做法是保留默认的缓冲电阻和电容但把缓冲电阻值调大到1e6欧姆让它在稳态时几乎不起作用只在数值上帮助解算器收敛。拓扑结构大致这样走同步电机端口A,B,C接三相互感或直接接测量模块Three-Phase V-I Measurement测量模块输出接示波器和“RMS”模块读取电压电流有效值电机端口再串一个Three-Phase Breaker初始状态为闭合断路器另一侧连接一个“Three-Phase Fault”模块或者直接接地这里我建议直接用Three-Phase Fault模块而不用断路器原因是Three-Phase Fault支持指定故障开始时间和结束时间还能方便地选择故障相并且内部带有接地电阻和故障电阻参数。断路器需要设计外部控制信号去切换状态麻烦一点。用Three-Phase Fault时把它的“Transition times”设置为比如[5 8]表示第5秒发生短路第8秒切除。为了模拟“突然”短路故障发生时间必须在发电机稳定运行之后比如从t0开始发电机逐步建压到t5s时基本稳定这时再触发故障。2.3 发电机参数设置用标幺值还是国际单位同步电机模块的参数输入方式有两种选择一种是国际单位SI另一种是标幺值pu。我看到不少人直接用电机铭牌参数换算成SI填进去结果仿真出来的波形怎么看怎么不对原因是SI参数换算涉及一堆绕组折算关系很容易出错。我建议直接用pu方式配合基准值设置这样既直观又不容易出错。具体参数填写实例以一台模拟的100MVA、10.5kV发电机为例额定视在功率Pn: 100e6 VA额定线电压Vn: 10.5e3 V额定频率fn: 50 Hz定子电阻Rs: 0.0025 pu定子漏感Ll: 0.2 pu这个值实际对应漏抗直轴同步电抗xd: 1.5 pu直轴暂态电抗xd: 0.3 pu直轴次暂态电抗xd: 0.21 pu交轴同步电抗xq: 1.5 pu隐极机取与xd相等交轴暂态电抗xq: 0.3 pu交轴次暂态电抗xq: 0.21 pu直轴暂态时间常数Td: 0.85 s直轴次暂态时间常数Td: 0.05 s交轴次暂态时间常数Tq: 0.07 s惯性时间常数H: 3.5 s极对数p: 2对应1500 rpm at 50Hz在这些参数里xd、xd、xd三个值非常关键。从电机学理论可知突然短路初瞬间定子电流大小主要由次暂态电抗xd限制随后次暂态分量很快衰减由Td决定过渡到由暂态电抗xd限制的阶段最后经过Td的衰减到达由同步电抗xd决定的稳态短路电流。你如果在仿真结果里看到电流峰值比用xd算出来的稳态值大很多那就是因为初始阶段实际限制电流的是xd。填写参数时要注意单位的坑Simulink的pu基准是额定电压和额定功率阻抗基准是 (Vn^2)/Pn。比如10.5kV、100MVA阻抗基准就是 10.5^2/100 1.1025 欧姆你填入的pu电抗乘以这个基准就是实际电抗值。但模块内部会自动换算不需要你手动算只要保证pu值正确就行。3. 核心细节暂态过程的理论公式与Simulink中的物理对应3.1 定子短路电流公式长什么样在忽略电阻、假设空载短路且A相磁链初始值为零时A相短路电流的解析表达式可以写成这里我直接用工程习惯的标幺形式[ i_a(t) \sqrt{2} E_q \left[ \left(\frac{1}{x_d} - \frac{1}{x_d}\right) e^{-t/T_d} \left(\frac{1}{x_d} - \frac{1}{x_d}\right) e^{-t/T_d} \frac{1}{x_d} \right] \cos(\omega t)\frac{\sqrt{2} E_q}{x_d} e^{-t/T_a} \cos(\theta_0)\text{二次谐波项} ]式中 (E_q) 为空载电动势(\theta_0) 是短路瞬间A相轴线与励磁磁场轴线的夹角(T_a) 是定子非周期分量的衰减时间常数。别看公式长每一项你都可以在Simulink波形里找到对应第一项方括号内对应基波周期分量的包络它由三个指数衰减项叠加而成分别对应次暂态、暂态、稳态三个阶段的电流幅值。短路后的第一周期内电流主要看第一项与第二项的叠加。第二项为非周期分量它使短路电流波形整体偏移时间轴表现为上半波大、下半波小。第三项二次谐波分量一般比较小在凸极机中更明显若你把xq和xd设成相等这项就消失了。在仿真中查看A相电流波形时你会发现短路后的最初几个周期波形明显偏离纯正弦波形一边高一边低这就是非周期分量的作用。非周期分量衰减到零后波形逐渐变成对称的稳态短路电流。3.2 为什么次暂态和暂态会有区别磁链守恒的直观理解为了讲清这个我用一个生活类比。转子励磁绕组和阻尼绕组就像两个“储物间”定子突然短路时外部磁场要挤进定子绕组但磁链守恒定律不允许定子绕组的磁链突变所以定子绕组里会立刻感应出电流来“顶住”外部磁场的变化。这个感应电流的路径同时受到定子漏感、转子励磁绕组、阻尼绕组的影响。在最初瞬间阻尼绕组和励磁绕组都相当于“超导”状态它们会感应出自己的电流来维持自己的磁链不变这导致定子电流的初始路径对应的等效电抗非常小也就是xd。随后阻尼绕组电流衰减Td很小通常0.03~0.15s它就相当于从“超导”变回普通绕组定子电流的限制电抗变为xd再之后励磁绕组电流衰减Td通常0.8~2s限制电抗慢慢过渡到xd。整个过程就是“短路电流先很大、然后快速掉一截、再慢慢掉到稳态”的三阶段。在Simulink里如果你用电机模块的内部模型这些时间常数都体现在参数Td、Td中不需要你额外加入任何惯性环节。但要注意Simulink的同步电机模型是基于Park方程建立的默认考虑了d、q轴之间的耦合所以仿真结果能够反映上述暂态过程。前提是你要确认自己填的Td和Td是“开路时间常数”而不是“短路时间常数”。电机学里给出的是开路时间常数也就是转子绕组开路时的时间常数通常用Td0、Td0表示。但从转子看进去的等效阻抗和定子侧短路状态有关实际短路时的衰减时间常数Td和Td是从开路时间常数换算来的。Simulink模块参数里写的是“d-axis transient open-circuit time constant Td0”还是“Td”这个要看清。如果你填的是开路时间常数仿真结果中的短路电流衰减会偏慢还是偏快需要对上。我用的模块版本里参数名是“d-axis transient open-circuit time constant Td0(s)”我填0.85s是开路值。如果你想验证短路时间常数可以用公式 Td ≈ Td0 * (xd/xd)按上面的参数就是 0.85*0.3/1.50.17s。在波形上你可以从次暂态衰减结束后到稳态之间的时间尺度来估算这个值能对上就说明模型合理。3.3 励磁绕组和阻尼绕组电流的变化趋势突然短路不只影响定子侧转子侧同样精彩。短路瞬间定子电流激增电枢反应突然增强为了维持转子励磁绕组磁链不变励磁绕组里会感应出一个很大的附加直流分量方向与原有励磁电流相同如果短路发生在空载且励磁电流为正值的时候所以励磁电流会突然猛增出现一个尖峰随后按Td和Td衰减回原值。阻尼绕组电流则在短路初瞬间最大随后按Td快速衰减到零。在Simulink中你可以通过电机的“m”端口输出信号找到励磁电流If和阻尼绕组电流。我实测的波形非常典型短路发生的瞬间If从1.0 pu猛跳到2.5 pu左右具体看参数然后快速衰减约0.2s后回落随后再缓慢衰减。这个现象在理论分析中常被提及但真正从仿真里看到冲击包络线时我对“磁链守恒”的理解才落到了实处。3.4 短路电流的非周期分量与合闸角的关系另一个让人容易忽略的细节是短路发生的相位角。定子绕组的磁链初始值与转子位置有关如果短路瞬间A相绕组的磁链正好为零那么A相非周期分量的初始值也为零A相电流波形上下对称如果短路瞬间磁链最大非周期分量初始值最大波形严重偏移。在Simulink的Three-Phase Fault模块里你可以设置故障发生时间但时间与相位的关系取决于发电机的初始相位和转速。如果你想专门研究合闸角的影响可以给发电机设置初始状态时调整转子初始位置。但更方便的做法是设置一个小的“Phase-angle shift”或者通过改变故障发生时间比如在50Hz系统中改变1/50s的时间对应改变360度的相位来观察不同合闸角。我一般通过把故障时间从5.000s逐步增加到5.020s每次增加0.005s相当于90度电角度就能看到A相非周期分量从最大到最小的全部变化。这个实验我强烈建议做因为理论课本上只会给一个包络图亲手改变相位看到波形偏移的过程对理解“磁链不能突变”会有质的提升。4. 仿真结果怎么看波形解读与理论验证4.1 定子三相短路电流波形搭好模型后我习惯先跑一个5秒的稳态过程然后在t5s时刻投入故障仿真持续到8s。第一眼看A相电流波形时你会发现短路前电流幅值很小空载短路后瞬间电流幅值暴增第一个峰值甚至接近额定电流的十倍。这个峰值大小和理论计算基本一致。以我用的参数为例额定电流基准 (I_n 100e6 / (\sqrt{3} * 10.5e3) \approx 5499A)。空载电动势大约1.0 pu因为励磁设置了1.0 pu对应空载额定电压。短路瞬间A相周期分量幅值为 (E_q / x_d 1.0 / 0.21 4.76 pu)考虑非周期分量的叠加第一个半波峰值可能达到5倍pu对应实际电流约3万安培左右。你从示波器里看到的第一个周波峰值应该在5pu附近而不是稳态短路电流的1/1.50.667pu。这就直观说明了次暂态电抗的重要性。用Scope读数时注意Scope显示的是瞬时值还是标幺值取决于你测量模块的设置。当电机模块输出定子电流时Simulink默认输出是pu你需要乘以电流基准才能得到安培数。建议在多路测量模块后面接一个“Gain”模块乘以基准电流或者直接改用“Machines”输出端口的is_abc信号然后自备标幺说明。不然你看到波形数值只有5容易误以为只比额定大5倍。4.2 电磁转矩与转速波动短路瞬间电磁转矩也会出现剧烈冲击。初始时转矩可能出现一个负向大尖峰对应发电机减速随后振荡衰减。在Simulink中你可以在发电机模块的m端口输出Te电磁转矩。你会看到转矩波形在短路瞬间出现一个接近额定转矩数倍的负向脉冲然后呈振荡衰减趋势。这个现象背后是定子电流非周期分量与转子励磁磁场共同作用产生的脉动转矩。如果只分析周期分量转矩通常稳定在某一负值但叠加了非周期分量后会产生一个按基波频率脉动的单轴转矩它正是导致发电机轴系扭振的主要激励。做机电暂态的同学可能更关注这个。如果以后往“次同步振荡”方向扩展这个仿真模型加一段轴系多质量块就能用。转速方面短路后原动机输入机械功率来不及调整发电机输出电磁功率骤降短路时电功率很小因此剩余功率会加速转子转速会缓慢上升。由于惯性时间常数H较大3.5s在0.5s的仿真窗口内转速变化可能只有百分之几但你能看到转速从1500rpm开始逐渐升高的趋势。如果原动机模型是恒功率的话。如果你用的是恒转速理想原动机转速就恒定不变那就看不到了。我这里为了让现象更真实用的是“Hydraulic Turbine and Governor”模块或者简化的常数机械功率Pm保持转矩不变。4.3 用FFT分析谐波构成不少人做完仿真就直接截图收工其实还可以用FFT工具进一步验证短路电流的谐波构成。Simulink的Powergui模块里自带FFT分析工具把A相电流信号加载进去选择第一个或第二个周波进行FFT你能清楚看到频谱上存在50Hz基波、100Hz二次谐波和直流分量。如果选择的FFT窗口在短路后0.02s以内直流分量和二次谐波都很明显如果在0.5s以后再看直流分量基本消失二次谐波也几乎没了基波幅值趋近于稳态值。这比单纯看波形更能说清楚“暂态分量”的含义。FFT操作上有个小提示Powergui中“FFT Analysis”要勾选“Display”和“Export”选项信号来源选“Scope”或“Structure with time”。推荐在故障测量模块后面接一个“To Workspace”模块存成structure格式这样FFT可以直接读取。如果你直接用Scope需要先把波形输出保存成变量操作路径有点绕。我一般是用“To Workspace”存变量一劳永逸。5. 实操中的常见问题与排查记录这个课题仿真本身不复杂但模型搭好之后往往会遇到各种奇怪现象。以下是我和身边人踩过的一些坑整理成表格方便你排查。现象可能原因解决办法仿真一开始就发散错误提示“Solver failed”同步电机初始化失败xd或Td等参数超出合理范围或使用了不合适的解算器改用“ode23tb”或“ode15s”并与Powergui一致检查参数pu值是否合理尤其是电抗介于0和1之间同步电抗可以大于1但漏感不能太大短路电流波形出现高频毛刺或锯齿解算步长太大断路器故障模块导致事件触发时数值突变把最大步长设为1e-4或更小在Powergui中勾选“Continuous”仿真类型不要用离散短路前电流不为零甚至很大你可能没有设置电机初始功率或负载空载仿真时电机端口接了负载或并联电容如果要做空载短路确保断路器之前没有负载故障模块在触发前应处于断开状态短路后电流波形不对称但没有衰减到零非周期分量衰减时间常数Ta被设得过大或定子电阻太小导致衰减太慢检查同步电机参数里的“stator resistance Rs”如果设得太小非周期分量衰减非常慢波形长时间偏移这是正常的但不是故障转速持续下降或上升发散了原动机输入功率与电功率不平衡或者调速器参数不合适先用恒定机械功率Pm额定值保证空载稳态时转速为同步速后面再加调速器模型Scope中看到多个波形重叠难以分辨各波形幅值差异太大电流几千A电压上万V分多个Scope显示或者用“Multimeter”加不同输出端合理设置坐标范围仿真结果与理论公式差1-2倍你可能用线电压、线电流的基准值但电机模块内部用的是相电压、相电流pu值基准选取要一致理论计算时用额定相电压和额定相电流作为基准然后再对比波形三相短路电流的包络线不是平滑衰减而是周期波动这其实不是错误因为非周期分量在每个相的衰减不同且转子的不对称也会导致包络线波动检查你是否设定了xd不等于xq如果是隐极机确认xdxq包络线会更平滑除了这些我还遇到过一类比较隐蔽的问题Simulink的Three-Phase Fault模块默认情况下接地电阻Rg为0.01欧姆。理想三相短路应该是阻抗为零但Rg为0.01欧会导致短路电流比理论值略小。如果要做精确对比把接地电阻改成1e-4甚至0。不过要注意直接把Rg设为0在某些求解器下可能引起数值奇异一般设为1e-6就可以忽略不计。另外有一个我经常用的小技巧在故障模块之前接一个“Three-Phase V-I Measurement”并用“Multimeter”选择测量类型为“Currents”和“Voltages”然后在Powergui的“Tools Measurements”中勾选并导出一系列电量。这比直接接Scope更能灵活记录不同位置的电气量。特别是记录定子A相电流的同时还想记录励磁电流和转矩就需要把机器信号和电路信号分别处理最后统一导入时间向量。6. 如何扩展这个仿真从教科书案例走向实际应用做完基础版本的突然三相短路仿真我个人觉得有四个方向特别值得继续往下做每个方向都能把模型价值放大不少。第一个方向是考虑励磁系统的影响。基本模型用的是恒定Ef短路暂态初期励磁调节器确实来不及动作但几个周期之后自动电压调节器AVR的输出就开始追踪电压偏差增大励磁从而影响短路电流的稳态值和衰减趋势。你可以把同步电机的Ef端口接到一个简化的AVR模型上用一个比例积分环节模拟。这样短路电流的稳态值不再由xd单独决定而是由励磁调节特性决定更接近实际电力系统。第二个方向是研究不同故障类型的影响。把所有故障类型都试一遍单相接地、两相短路、两相接地短路。你会发现单相接地短路时由于零序通路和序分量叠加电流特性与三相短路差别很大。此时需要用序分量理论来解释也可以在Simulink中观察各序电流波形。这个方向能帮你把“对称分量法”从理论带到实践。第三个方向是联网运行下的短路暂态。把同步发电机接一个升压变压器再接到无穷大电网在变压器高压侧设置短路故障。这时你需要加入变压器模型、合理设置基准值发电机与电网之间会发生功率交换短路电流分析会更复杂也会涉及保护定值配合。这是与实际工程最贴近的一个扩展。第四个方向是机电暂态与轴系扭振分析。前面提到短路瞬间电磁转矩冲击很大如果你把同步电机模块扩展到多质量块轴系比如高压缸、低压缸、发电机、励磁机多质量块弹簧模型就能研究短路对轴系的扭振激励。这会用到电力系统机电暂态与转子动力学交叉的知识属于比较进阶的方向。我个人在实际操作中的体会是Simulink这个模型的价值不在于它能不能还原教科书上的波形而在于它逼着你去确认每个参数的物理意义。如果你只是填一组默认参数、跑一条曲线那就失去了这个课题的意义。把xd、xd、xd、Td、Td挨个改一遍观察短路电流包络线怎么变比看十篇论文都管用。最后再分享一个小技巧仿真时间设置上不要只跑到短路后0.1s就停。即使你只关心第一个峰值也建议跑到故障后2s左右让波形充分衰减到稳态。这样你能一次性在图上看到“次暂态—暂态—稳态”三个阶段的界限后续做FFT和参数提取也方便。真正看到那三段衰减曲线的时候突然三相短路这件事才算真正吃透了。
