从0.02s看起的短路冲击无穷大功率电源供电系统三相短路Simulink仿真全解析先说个结论三相短路仿真在电气工程里不算稀奇但真正把模型调明白、波形看透的人并不多。随便搜一下Matlab/Simulink短路仿真能找到大批“带跑通代码”的资源可一问到“为什么0.02s这个时刻短路”“变压器参数为什么这样设”“为什么仿真一跑就发散”能讲清楚的人就少了大半。这篇就用无穷大功率电源给变压器供电、0.02s时低压侧发生三相短路这个场景把Simulink建模、参数整定、结果解读到坑点排查完整捋一遍希望给正在做课程设计、毕业设计或者打牢短路分析基本功的同学一点实质参考。无穷大功率电源这个概念很多人觉得它就是“电压恒定不变”的代名词。实际仿真里没有真正“无穷大”的电源只有“等效内阻小到可以忽略”的电源模型。之所以用无穷大电源是因为短路发生在变压器低压侧时系统侧等值阻抗远小于变压器短路阻抗此时短路电流主要受变压器漏抗限制电源侧的电压幅值和频率在故障前后几乎不变化。这个前提让短路电流的解析计算变得非常简单也方便跟仿真波形互相验证——这正是这个仿真的最大价值用仿真结果反过来验证书本上的解析公式。我会按模型设计思路、模块搭建细节、理论计算对照、波形解读和问题排查五个部分来讲全程用Simulink标准库里的模块不涉及自定义小工具照着搭就行。1. 场景设计与数学原理先搞清楚“在算什么”1.1 无穷大功率电源是什么为什么仿真里还要设内阻“无穷大功率电源”是一种理想化建模。它的定义是电源内部阻抗为零端电压在任何扰动下保持恒定频率恒定不变。现实中并不存在但电网中很多场合可以近似为无穷大系统。比如一个几十万千伏安的220kV变电站对于它下游一台几千甚至几万千伏安的配电变压器而言当这台变压器低压侧发生短路时220kV系统自身的电压降极其微小把它当作恒定电压源完全够用。仿真里有个细节要注意Three-Phase Source模块参数里如果填内阻为0会在部分求解器下引起数值奇异出现仿真发散或者初始潮流计算失败。所以即使做“无穷大”模型依然要填一个极小的阻抗比如0.0001Ω或更小模拟“近似无穷大”。这里不要较真说“无穷大为什么有内阻”这是数值计算的妥协也是行业惯例——大多数教科书例题、工程短路计算也会给一个“系统短路容量无穷大”或“系统等值阻抗忽略不计”的前提仿真里的极小内阻就对应着这个假设。1.2 三相短路的暂态过程周期分量与直流分量三相短路是电力系统中最严重的故障类型之一因为三相同时短接故障点电压几乎为零短路电流是三相对称的没有负序和零序分量所以只用正序等值电路就能计算。短路发生瞬间的电流和稳态短路电流往往是两码事。稳态短路电流由系统等值阻抗决定大小恒定但短路暂态过程里由于电感中磁链不能突变电流里会叠加一个衰减的直流分量。这个直流分量的大小取决于短路瞬间的电压相位角——电压过零附近短路时直流分量最大电压峰值附近短路时直流分量最小。直流分量随时间按指数规律衰减衰减快慢取决于短路回路的X/R值电抗与电阻之比X/R越大衰减越慢。这就是为什么仿真波形里短路瞬间会看到一个很大的冲击电流尖峰过几个周波后才慢慢回落到稳态值。0.02s这个时间点对应中国电网工频50Hz的一个完整周波也就是电压相位经过360度回到起点。选择这个时刻做短路投切波形观测起来最规整也方便和理论计算的冲击系数对应。1.3 冲击电流与最大有效值电流的计算逻辑理论分析中冲击电流即短路后半个周波时出现的最大瞬时值的计算公式是[ i_{sh} K_{sh} \times \sqrt{2} \times I_{sc} ]其中 (I_{sc}) 是稳态短路电流的有效值(K_{sh}) 是冲击系数(K_{sh}1e^{-\pi / (X/R)})。如果X/R很大Ksh接近2意味着冲击电流峰值接近稳态短路电流幅值的两倍如果X/R较小比如接近纯电阻Ksh接近1冲击就不太明显。这一点在仿真验证中很重要仿真波形上读出来的短路电流最大瞬时值应当落在 ( \sqrt{2}I_{sc} ) 和 ( K_{sh}\sqrt{2}I_{sc} ) 之间并且半个周波0.01s后达到最大。如果你的波形峰值远高于理论冲击值说明模型参数有误——大概率是没等效对系统阻抗或者故障电阻设成0导致极端短路。2. Simulink模型搭建模块选型与参数设置2.1 用到的模块清单与模型结构Simulink里搭建这个仿真不需要额外安装工具箱Simscape Electrical早前叫SimPowerSystems在标准版Matlab里自带。核心用到的模块如下Three-Phase Source三相电压源模拟无穷大系统。这里要设置幅值、频率、内阻抗。Three-Phase Transformer (Two Windings)双绕组三相变压器Y/Y或Y/Δ连接均可需要填容量、频率、两侧电压、短路阻抗百分比、铜耗/铁耗参数。Three-Phase Fault三相故障模块可以在0.02s时刻将三相同时接地短接。Three-Phase V-I Measurement电压电流测量模块接在变压器低压侧与故障模块之间用来抓短路电流波形。Three-Phase Series RLC Branch如果系统需要等值阻抗来近似短路容量可以用这个模块串在电源和变压器之间。Scope示波器看波形就靠它。整体连接顺序是Three-Phase Source → 串联阻抗如果需要→ 变压器高压侧 → 变压器低压侧 → V-I测量 → Three-Phase Fault → 负载可选。短路故障时Three-Phase Fault内部开关闭合相当于低压侧三相直接短接。2.2 变压器参数怎么填标幺值是关键变压器模块参数设定是整个模型里最需要动脑子的地方。以容量100MVA、电压等级110kV/10.5kV的变压器为例Nominal power额定容量100e6 VANominal frequency额定频率50 HzWinding 1 (HV) connectionYg高压侧中性点接地或直接Y接Winding 2 (LV) connectionDelta低压侧三角接这是配电变压器最常见的接法Winding 1 voltage (V1)110e3 V线电压有效值Winding 2 voltage (V2)10.5e3 V线电压有效值Magnetization resistance Rm 和 reactance Xm通常在100倍标幺值附近默认值即可Short-circuit resistance 和 leakage reactance这两个参数要和短路阻抗百分比Uk%对应起来关于短路阻抗变压器铭牌上一般标Uk%10.5%意思是在额定频率下把一侧短路另一侧施加电压达到额定电流时的电压标幺值为0.105。换算成标幺阻抗就是0.105 pu。其中电阻分量很小通常电抗占比高电阻/电抗比在1/10到1/20之间。如果只填了Uk%但没填电阻电抗分别多少模块里可以用“resistance and reactance”选项分开填或者简化填“resistance and leakage reactance, in pu”并直接给电阻值比如0.01 pu和电抗值比如0.1044 pu算出 ( \sqrt{0.01^20.1044^2} \approx 0.1049 ) pu非常接近0.105。注意变压器短路阻抗的正确定义是决定短路电流大小的核心。短路电流计算时变压器低压侧三相短路等效到高压侧看进去回路的阻抗就是系统阻抗加变压器漏抗。2.3 电源等效内阻的数值处理三相对称短路时无穷大电源的内阻可以忽略不计但模块上得填一个有限值。我习惯填1e-4欧姆的电阻串联1e-3亨的电感这样数值计算不会报错且对短路电流的影响远小于变压器阻抗仿真结果与理想无穷大电源之间的误差在0.5%以内完全可以忽略。也有人直接在Three-Phase Source里填0模型也能跑通但遇到某些求解器在高频暂态时容易出数值振荡。为了少踩坑还是填个极小阻抗稳妥。2.4 Fault模块的配置细节Three-Phase Fault模块是触发短路的关键。双击打开后关键设置如下Fault type故障类型选Three-phase faultABC这就是三相短路。External control of fault timing如果要在仿真中精确控制0.02s时刻用内部timing即可。Transition status初始状态选open不故障0.02s时切换为closed故障。Transition time 1s填0.02。Fault resistance R fault (ohm)一般填0.001或更小的过渡电阻。理论上三相短路金属性短路的过渡电阻为零但设为0可能引起某些求解器的数值振荡建议填一个极小值如1e-3欧姆。Snubber resistance默认1e6欧姆保持默认。2.5 求解器设置与步长选择Simulink模型搭好之后求解器设置直接决定仿真能不能收敛、波形准不准。电气系统是典型的刚性系统包含电磁暂态过程中功率电子开关和电感的快动态与机械慢动态的混合建议Solver选Variable-step变步长Solver type选ode23tb或ode15s。电气仿真中ode23tb往往比ode45更稳定尤其开关动作频繁或系统阻抗比很大的场合用ode45经常会出现“仿真发散”的报错核心原因就是刚性太强。Max step size建议设为1e-4秒或更小。虽然变步长求解器会自动调整步长但如果限制最大步长过大的话在0.02s短路时刻会产生很大的局部截断误差波形峰值和相位的精度都会下降。Relative tolerance保持1e-4默认即可对绝大多数电气仿真已经足够如果波形有高频毛刺可以适当缩小到1e-5。3. 短路电流理论计算把公式和仿真值对上3.1 用标幺值法快速算出稳态短路电流在做仿真之前先把理论值算出来这样拿到仿真波形之后可以直接对比立刻能判断模型对不对。以100MVA、110kV/10.5kV变压器Uk%10.5%系统无穷大容量为例变压器低压侧额定电流为[ I_{2N} \frac{S_N}{\sqrt{3} \times U_{2N}} \frac{100 \times 10^6}{\sqrt{3} \times 10.5 \times 10^3} 5499 \text{ A} ]低压侧三相短路时由于系统侧阻抗忽略短路电流主要受变压器漏抗限制。用标幺值算更直观短路电流标幺值等于1除以短路阻抗标幺值即 ( I_{sc}^* 1 / 0.105 9.52 ) pu。那么实际短路电流为[ I_{sc} 9.52 \times 5499 52351 \text{ A} \approx 52.4 \text{ kA} ]这个52.4kA就是稳态短路电流的有效值。如果变压器高压侧还有一段线路阻抗或系统等值阻抗那这个值会更低一些。3.2 估算冲击电流的峰值量级要估算冲击电流峰值先算回路X/R。假如变压器短路阻抗里电阻分量R0.01pu、电抗X0.1044pu那么X/R约为10.44。冲击系数为[ K_{sh} 1 e^{-\pi / 10.44} 1 e^{-0.301} 1 0.74 1.74 ]冲击电流峰值为[ i_{sh} K_{sh} \times \sqrt{2} \times I_{sc} 1.74 \times 1.414 \times 52351 \approx 128.8 \text{ kA} ]看到这个数值不要慌它是短路后半个周波0.02s 0.01s 0.03s附近电流的瞬时最大值。稳态幅值大约为 ( \sqrt{2} \times 52.4 74.1 ) kA冲击电流比稳态幅值高了将近一倍这就是直流分量的“加持”。仿真波形上如果峰值落在125kA到130kA之间说明模型设得很准。注意冲击电流是短路电流的最大瞬时值也是校验电气设备动稳定性的关键参数。断路器开断能力、母线机械强度校验用的都是它。仿真里读出峰值后可以顺手算一下“峰值系数”是否在1.8附近这能快速验证模型是否合理。3.3 为什么说0.02s是“最不利”的短路时机的近似不是随便选的0.02s对应50Hz工频的一个整周波。在这个时刻投入短路A相电压相位与0时刻相同但三相中总有一相电压接近过零点。而直流分量最大的情况正是发生在电压过零附近短路时。所以用0.02s作为故障时刻虽然不是数学上让某相直流分量取最大值的精确角度但已经足够观测到直流分量对冲击电流的明显抬升。如果把短路时刻改到电压峰值附近例如0.025s对应5/4周波你会在波形上看到冲击电流明显减小直流分量很小几乎看不出暂态凸起。很多教材喜欢强调“短路瞬间电压相位对暂态过程的影响”这个仿真模型的优势就在于可以通过改动“Transition time”一个参数直观看到同一系统在不同短路相位下电流波形的差异——这是静态计算给不了的。4. 波形仿真结果解读Scope里看到的到底是什么4.1 全时段观测短路前、短路瞬间、暂态衰减仿真时长我习惯设为0.1s0.02s短路之后留0.08s观察暂态衰减和稳态。Scope里同时拉出三相电流和低压侧线电压波形。短路前如果是空载状态电流为零如果是带负载状态电流为一个较小的正弦波。0.02s时三条电流曲线突然窜高出现明显的暂态冲击随后在几个周波内逐渐“收敛”到对称的正弦波幅值稳定在52.4kA附近。具体看电流波形有这几个特征值得盯着观察短路瞬间电流有“毛刺”和较大尖峰这是正常现象不是仿真错误它表示直流分量叠加在交流分量上。三相对称短路电流波形互相保持120度相位差稳态时三相对称暂态过程中三相的直流分量大小不同导致三相波形呈现“不对称”的偏移。冲击电流峰值出现的时刻应该在短路后半个周波附近但仿真中由于直流分量衰减、系统参数不完全理想等因素峰值可能略有偏移这是正常的。4.2 用测量工具读取关键数值Scope模块自带光标测量功能仿真结束后点开Scope在波形窗口的工具栏里选择“cursor measurement”就能看到当前游标位置的数值。如果觉得分辨率不够可以在Scope左上角设置里减小时间轴显示范围把0.02s到0.035s这段放大再用光标读取。我实测过一组参数110kV/10.5kV、100MVA、Uk%10.5%、电源内阻0.01Ωj0.1Ω仿真得到的稳态短路电流有效值约51.8kA冲击峰值约127kA与理论值52.4kA/128.8kA的误差在2%以内。这个精度已经足够验证模型的正确性。4.3 电压波形低压侧母线塌陷与高压侧维持变压器低压侧三相短路后低压侧电压几乎降为零。Simulink的测量模块会显示一个接近零的电压波形在故障瞬间存在一个急剧跌落的过程。而高压侧电压因为无穷大电源的钳制基本保持不变——这个特征完美体现了“无穷大系统的电压支撑能力”。如果你把高压侧电压也拉进Scope就会看到故障前后高压侧电压幅值几乎是一条平行线只在故障瞬间因为数值计算产生微小的抖动抖动幅度通常在1%以内。这个现象在实际工程中对应的概念是“母线残压”。低压侧短路后母线电压几乎为零该母线上其他负荷全部失压这也是短路故障对系统最直接的影响之一。4.4 参数扫描看看X/R大小对波形的影响模型搭好以后我特别建议做一组参数扫描修改变压器短路阻抗中的电阻与电抗比值观察冲击电流的变化。方法是对同一个模型复制几个副本分别设X/R比为2、10、20运行仿真把Scope重叠对比。你会很直观地看到X/R越大短路电流暂态过程越长、冲击峰值越高X/R越小冲击峰越低、直流分量衰减越快。这个练习比单纯读一百遍教材的“阻尼作用”都有效。因为X/R比越大回路时间常数 ( L/R X/(\omega R) ) 越大直流分量衰减越慢冲击系数越接近2。在模型里改几个数字、跑一次仿真这个物理规律就刻在脑子里了。5. 常见问题与排查技巧仿真发散不是玄学5.1 “仿真发散”的常见原因和解决方法“仿真发散”是电气仿真里出现频率最高的报错我自己初学阶段至少被它折磨过几十次。在无穷大电源三相短路这个模型里仿真发散的原因主要有这么几类求解器选错。使用ode45或ode23求解刚性电气系统时步长会变得非常小仿真时间却推不进甚至直接报数值奇异。换ode23tb或ode15s后大多数情况下立刻好转。最大步长设得太大。设成1s甚至默认值会让短路瞬间的开关切换产生巨大的局部误差。建议最大步长不高于1e-4秒如果电脑配置好可以压到5e-5秒。电源内阻为0或变压器短路阻抗为0。阻抗为零会引入退化方程导致数值矩阵奇异。任何阻抗都不要填0哪怕填一个1e-6的小数都行。故障电阻设为0时某些求解器也会出问题解法是填0.001Ω或更小的正数。提醒如果反复调参仍然发散先检查模型连线是否正确尤其是Three-Phase V-I Measurement模块的方向箭头是否朝向负载/故障一侧。方向接反会出现某一路电流为负这本身不导致发散但会引起后续分析混淆。5.2 波形上出现大量高频毛刺这种情况一般不是模型错了而是数值刚性引起的振荡。可能原因包括最大步长太小、相对误差限太小、求解器不适合高频分量。处理办法是换ode23tb适当放宽相对误差到1e-3或者打开Simulink配置里的“zero-crossing detection”为Enabled。还有一种情况是Scope采样时间太短导致看起来像毛刺这种情况下调小Scope的采样间隔或者改用decimation采样即可。5.3 短路后电流反而下降或者波形不对称这通常是变压器磁饱和模型或励磁支路参数设置不当引起的但更常见的原因是故障模块设置了故障相不正确比如选了A/B相而漏了C相。三相短路必须选ABC三相不然波形必然不对称这与题目设定不符合。另有一种情况是变压器本身Y/Δ接法导致低压侧绕组电流和线电流同相位不一致Scope里三根电流线看起来相位和幅值有点差别这可能不是故障而是接线组别的影响。判断办法是看低压侧线电流Line current而不是相电流Phase current两者在Δ接法下有根号三倍的关系别弄混。5.4 想要看不对称故障怎么办这是一个很自然的扩展方向把Three-Phase Fault的故障类型改成“A相接地”或“AB相间短路”仿真立刻就能看到不同故障类型的波形特征。三相短路与不对称短路最直观的区别是三相短路波形三相对称幅值最大单相接地故障有零序分量故障相电流未必最大相间故障两相电流升高、另一相为零。这些差异用同一套模型改一个选项就可以看到强烈建议把几种故障类型都跑一遍。收尾的一点体会这个模型我调过很多次最深刻的体会是书本上那些公式和系数只有在亲手拉过波形、量过峰值之后才真正内化成自己的东西。无穷大电源供电系统三相短路仿真看似简单但它把短路电流计算、暂态过程分析、电磁暂态仿真调试这些基本功全串起来了。建议所有刚接触电力系统仿真的同学都不要跳过理论计算这一步先把52.4kA和128.8kA算出来再跑仿真你会发现自己对“仿真结果对不对”的判断力直接上一个台阶。最后再分享一个小技巧保存波形时别只存Scope截图用“To Workspace”模块把电流时间序列导到MATLAB工作区然后就可以用plot命令自定义绘图了。做课程设计或写论文时自己控制线宽、坐标标签、图例的图比Scope截图专业很多而且还能顺手算一下FFT谐波含量——这又是另一篇可以展开的内容了。