这次我们来看电力系统三相短路的 MATLAB Simulink 仿真怎么搭。短路暂态分析是继电保护整定、开关设备选型、电力系统稳定研究的基础而三相短路又是最严重、最常被拿来做首轮研究的故障形式。很多做毕设、课程设计或者刚接触电力系统仿真的人都会卡在同一个问题上Simulink 模型到底怎么建故障模块在哪里参数该怎么填仿真出来的波形怎么解读。这篇文章会从零开始带你在 Simulink 中完成一个包含三相电源、输电线路、短路故障模块、测量环节的暂态仿真模型然后演示怎么样把短路电流波形录出来、读取冲击电流值、对比不同故障类型的结果并给出一个用 MATLAB 脚本批量扫描故障时刻的示例方便你做参数分析和后续扩展。全文主要面向电气工程学生、继电保护方向研究生、刚转电力系统仿真的一线工程师以及任何需要在 MATLAB 里快速验证短路暂态过程的人。不需要先具备 Simulink 高级建模经验能看懂常用模块和基础连线即可但需要有 Simulink 与 Simscape Electrical 相关工具箱的 MATLAB 环境。1. 核心能力速览能力项说明项目类型电力系统暂态仿真模型搭建与验证使用工具MATLAB Simulink Simscape ElectricalSpecialized Power Systems主要功能三相短路暂态过程仿真、短路电流波形记录、冲击电流读取、不同类型短路对比模型复杂度中等适合教学、课程设计、科研预研硬件要求普通台式机 / 笔记本即可CPU 为主内存建议 8GB 以上工具箱要求Simulink、Simscape、Simscape Electrical需安装完整许可输出效果Scope 波形、To Workspace 数据、MATLAB 后处理图表接口与批量可用 sim 命令批量仿真可通过 set_param 修改参数实现参数扫描适合场景继电保护整定分析、暂态特性教学、开关设备选型前期研究2. 适用场景与使用边界这个仿真模型解决的核心问题是在没有实际实验设备的情况下用数字仿真还原三相短路发生时电力系统的暂态行为。它适合做继电保护原理验证、短路电流衰减特性分析、不同中性点接地方式对短路电流的影响、以及保护装置的动模试验前校验。不过要特别注意使用边界。Simulink 里用集中参数线路和理想故障模块建立的模型不等同于实际电网。它忽略了很多工程细节比如线路的分布参数特性、电弧动态特性、变压器饱和、实测故障电阻变化等。因此仿真结果可以用来做相对趋势分析但不能直接替代实际系统的短路计算更不能直接把仿真得到的短路电流数值拿去做现场保护定值整定。涉及真实电力系统数据时参数应脱敏处理并且最终结果必须由专业工程师结合实测和业界标准计算工具校验。在法律合规方面使用 MATLAB Simulink 做常规电力系统仿真不涉及敏感内容但如果你的模型基于实际变电站或电网数据请务必确保数据来源合法不要公开未脱敏的网络拓扑、设备铭牌和运行参数。仿真模型本身可以用来做保护逻辑测试但严禁将模型输出未经校验地用于真实设备投运决策。3. 环境准备与前置条件开始搭建之前先确认你的 MATLAB 环境。比较稳妥的配置是 R2020a 及以上版本因为 Simscape Electrical 在较新版本中刀具更完整模块搜索也更方便。如果你还在用 R2018 或者更早的版本也能做但部分模块路径会不同需要手动查找。打开 MATLAB 后在命令行窗口输入ver查看已安装的工具箱列表中是否有下列内容SimulinkSimscapeSimscape ElectricalMATLAB Report Generator非必需但有助于导出仿真报告如果没有 Simscape Electrical最直接的解法是重新运行 MATLAB 安装程序勾选“Simscape Electrical”组件或者联系管理员开通许可证。注意这个工具箱早期叫 SimPowerSystems后来并入 Simscape所以你在较老教程里看到“SimPowerSystems”字样指的就是同一套电力系统仿真库。环境检查建议清单检查项要求校验方式MATLAB 版本R2020a 及以上越低版本兼容性越差version命令Simscape Electrical必须安装ver查看工具箱列表内存建议 8GB 以上模型不大时 4GB 可运行系统任务管理器磁盘空间至少预留 5GBMATLAB 运行时需要缓存系统磁盘属性路径权限模型文件和脚本所在目录可写新项目放在独立文件夹中4. Simulink 三相短路仿真模型搭建下面进入核心部分。目标是建立一个10kV三相电源经过一段输电线路后发生三相短路的模型故障由Three-Phase Fault模块控制波形用Scope显示同时用To Workspace把电流数据导出到 MATLAB 工作区方便后续计算冲击电流和批量扫描。4.1 新建模型并添加模块在 Simulink 起始页点击“Blank Model”新建模型然后打开库浏览器Library Browser搜索并添加以下模块模块名称数量作用Three-Phase Source1提供三相交流电源Three-Phase Series RLC Branch1模拟输电线路阻抗Three-Phase Fault1控制三相短路发生和切除Three-Phase V-I Measurement1测量母线三相电压和电流Scope2显示电压和电流波形To Workspace2把电压电流导出到工作区Powergui1电力系统模型求解器设置与初始化把模块拖到画布后按图 1 的思路连线Three-Phase Source 输出端接到 Three-Phase V-I Measurement 的输入侧V-I Measurement 的输出端串联 Three-Phase Series RLC Branch线路输出端接到 Three-Phase Fault 模块的输入端和输出端电测量模块的测量输出端接 Scope 和 To Workspace。注意交流电源、线路、故障模块之间要使用 Simscape Electrical 电气连接端口波形显示和 To Workspace 需要先经过信号转换节点一般是 V-I Measurement 输出的信号端口。4.2 电源与线路参数配置双击Three-Phase Source按下表思路设置Phase-to-phase RMS voltage10e3VFrequency50HzInternal connectionYg星形接地Source short-circuit level可以作为系统等值阻抗的简化体现不设置时按理想电压源处理双击Three-Phase Series RLC Branch时需要先确认使用方式是串联 RLCResistance R1ΩInductance L10e-3HCapacitance Cinf线路对地电容暂不考虑时设为 inf相当于开路这几个参数决定了短路发生后的时间常数。线路的 X/R 比值直接影响非周期分量的衰减速率X/R 越大短路电流直流分量衰减越慢。你后续做批量分析时可以替换成不同 L/R 组合观察波形差异。4.3 故障模块配置双击Three-Phase Fault关键参数直接决定仿真行为参数名建议值说明Fault typeABC 或 A-G、A-B 等选择短路类型三相短路选 ABCTransition Times[0.2 0.3]第一个是故障开始时刻第二个是故障切除时刻Fault resistance0.001Ω故障电阻越小短路电流越接近金属性短路Ground resistance0.001Ω接地故障时的接地电阻三相短路也可保留默认这里0.001欧姆的故障电阻既能模拟接近金属性短路的特性又能避免纯零电阻导致部分求解器出现数值振荡。注意故障模块的 Transition Times 是两个时间点格式是类似[0.2 0.3]的数组表示故障从 0.2 秒发生、0.3 秒切除。如果只填一个时间点表示故障在该时刻发生并永久保持。4.4 测量与波形显示Three-Phase V-I Measurement的测量端输出一组电压和电流信号直接接到两个 Scope 上就能看到三相电压和三相电流波形。为了让数据进入 MATLAB 工作区同时接两个To Workspace模块Variable name分别设为Vs_abc和Is_abcSave format选TimeseriesSample time设为0继承输入信号步长这里建议优先用 Timeseries 格式因为后续在 MATLAB 脚本里处理起来最方便可以直接拿到时间和数据值。4.5 Powergui 求解器设置电力系统模型里必须加入Powergui模块否则仿真会报错或结果异常。双击 Powergui 进入设置Simulation type 选ContinuousSolver type 选Variable-stepSolver 选ode23tb或ode15s电力系统短路暂态属于典型的刚性系统短路瞬间电流变化非常快非周期分量缓慢衰减这两段时间尺度差异很大。ode23tb和ode15s都是 MATLAB 里处理刚性问题的常见求解器能有效避免步长过小导致仿真卡死。5. 功能测试与效果验证模型搭好后不要直接一上来就跑短路工况。更合理的顺序是先验证无故障稳态再叠加短路这样才能判断波形异常到底是模型本身问题还是故障引入的。5.1 无故障稳态校验先把Three-Phase Fault模块的 Transition Times 临时改成一个远超仿真结束时间的值比如[999 1000]相当于整个仿真过程中故障都不发生。设置仿真停止时间为0.2秒点击运行。预期结果是三相电压和三相电流均为标准正弦波三相电流幅值由线路阻抗决定幅值基本恒定。这一步验证了电源、线路、测量链路是通的。如果波形畸变明显或者电流幅值异常大先回去检查线路阻抗与电源电压是否匹配比如线电压 10kV 时1Ω 线路电阻加 10mH 电感下稳态电流应该在几百安到上千安级别如果看到上百万安培的电流大概率是电源内阻没设置并且线路阻抗设成了接近 0。5.2 三相短路暂态波形恢复故障模块 Transition Times 为[0.2 0.3]仿真停止时间设为1秒。再次运行。这里重点看短路点附近的电流也就是Is_abc。预期你可以看到0.2 秒前三相电流正常正弦。0.2 秒时刻故障相电流突然激增出现很大的冲击峰值。第一个峰值后电流不是立刻变成对称正弦而是叠加了一个按指数规律衰减的直流分量导致波形的正负半轴不对称整体包络线逐渐回落到稳态短路电流水平。0.3 秒故障切除后电流快速恢复或过渡到新的稳态。这就是典型的短路暂态过程。观察非周期分量衰减快慢时可以对比线路电阻 R 和电感 L 的不同组合R 越大衰减越快反映到波形上就是不对称很快消失R 很小、L 很大时直流分量能维持好几个周波。5.3 用 MATLAB 脚本读取冲击电流波形只能靠“看”要量化分析就需要把数据拿出来算。如果你在模型中设置了To Workspace运行后工作区里已经有Is_abc这个 Timeseries 对象。在命令行窗口执行下面的示例脚本即可提取故障期间的最大电流绝对值也就是近似冲击电流严格意义上冲击电流还跟衰减时间常数有关但作为仿真学习足够% 读取短路电流数据 t Is_abc.Time; Ia Is_abc.Data(:, 1); % A 相电流列号视电压电流顺序调整 Ib Is_abc.Data(:, 2); Ic Is_abc.Data(:, 3); % 提取故障时间段 fault_start 0.2; fault_end 0.3; idx t fault_start t fault_end; % 计算三相电流峰值 peak_a max(abs(Ia(idx))); peak_b max(abs(Ib(idx))); peak_c max(abs(Ic(idx))); disp([A 相故障期间最大电流绝对值: , num2str(peak_a), A]); disp([B 相故障期间最大电流绝对值: , num2str(peak_b), A]); disp([C 相故障期间最大电流绝对值: , num2str(peak_c), A]);注意Is_abc.Data(:, 1)的具体索引取决于 V-I Measurement 输出的信号顺序和 To Workspace 模块的接线顺序。如果运行后报数组越界直接在命令行执行size(Is_abc.Data)查看实际维度再调整列号。把算出来的峰值除以故障前的稳态电流幅值就能得到冲击倍数。电力系统教材里常提的这个倍数通常与系统 X/R 比值有关你的仿真结果应该能直观反映这种关联。5.4 不同短路类型对比三相短路数值最大但实际电力系统中单相接地短路发生概率更高。回头双击Three-Phase Fault把 Fault type 改成A-GA 相接地再仿真对比 A 相电流和三相短路时的差异。如果改成A-B两相短路可以看到非故障相 C 相电流明显小于故障相并且故障相的相位关系也不一样。这个对比实验对理解对称短路和不对称短路的暂态差异非常有帮助是理解后续基于序分量保护原理的重要基础。6. MATLAB Simulink 求解器与性能观察Simulink 仿真电力系统时性能表现跟求解器选择关系很大。很多初学者遇到仿真卡住或波形毛刺多第一反应是模型错了但实际上问题往往出在求解器配置。求解器选择的通用原则场景推荐求解器说明普通连续系统ode45速度快精度可接受电力电子 / 短路暂态ode23tb 或 ode15s刚性系统能承受故障瞬间的快速变化离散系统 / 实时仿真固定步长离散求解器需要先通过 Powergui 进行离散化设置单纯看稳态结果加大步长 ode45可以提高计算速度但暂态细节会丢失如果短路瞬间波形产生明显的高频振荡毛刺优先检查仿真步长上限。在 Configuration Parameters 里最大步长不要设得太大理论上建议不大于系统最小时间常数的四分之一。10mH 电感和 1Ω 电阻对应的电气时间常数为 10ms最大步长建议设在 0.001 秒级别太大会把衰减过程直接抹掉。仿真速度的观察方法也很简单。在 MATLAB 命令行运行脚本时用tic和toc包裹sim命令tic; out sim(mdl, StopTime, 1); toc;如果一次 1 秒仿真时间要跑几十秒甚至几分钟需要看模型是否有大量Scope在实时更新画面。Scope 实时绘制对 CPU 开销很大做批量仿真时不要开着好几个 Scope 窗口跑改用To Workspace记录数据仿真结束后统一画图。另外一个隐藏开销是 Powergui 里如果开启了FFT Analysis或Impedance Measurement仿真前必须关闭或断开否则每次积分都要消耗额外资源。内存占用方面常规的小规模三相短路模型不会把机器内存吃满。真正吃内存的是你把仿真数据全部存成 Timeseries并且仿真时间很长、步长很小时返回的数组可能非常庞大。实际编程时建议只输出关键通道或者在 To Workspace 里设置降采样 Sample time例如设为0.0001秒而不是0。7. 批量仿真与数据导出单次仿真只能得到一个工况的结果做课程设计或写论文时通常需要分析不同故障时刻、不同故障电阻、不同线路阻抗对短路暂态的影响。手动改参数再点运行效率低而且容易漏掉组合。正确的做法是写一个 MATLAB 脚本用set_param修改模型参数再用sim命令循环运行模型。以批量扫描故障发生时刻为例。假设故障模块名是Three-Phase Fault模型名是three_phase_fault_demo脚本可以这样写% 选择模型 mdl three_phase_fault_demo; open_system(mdl); % 需要扫描的故障发生时刻 fault_times [0.05 0.1 0.2 0.3 0.5]; fault_duration 0.1; peak_current zeros(size(fault_times)); for k 1:numel(fault_times) t_start fault_times(k); t_end t_start fault_duration; % 修改故障模块的切换时间 set_param([mdl /Three-Phase Fault], ... SwitchTimes, sprintf([%g %g], t_start, t_end)); % 运行仿真 simOut sim(mdl, StopTime, 1); % 提取 A 相电流并计算故障期间峰值 Is_abc simOut.get(Is_abc); t Is_abc.Time; Ia Is_abc.Data(:, 1); idx t t_start t t_end; peak_current(k) max(abs(Ia(idx))); end % 汇总结果 results_table table(fault_times, peak_current, ... VariableNames, {FaultTime, PeakCurrent}); disp(results_table);这段脚本里simOut.get(Is_abc)的变量名必须和模型中To Workspace设置的变量名一致。如果你没有勾选模型配置里的“Single simulation output”选项sim命令返回的结果结构会有一点差别最稳妥的情况是调用前用get_param(mdl, SaveOutput)查看。报错时在命令行执行help sim确认当前 MATLAB 版本的输出格式。另外一个容易踩的坑是set_param修改参数后如果模型处于运行状态或者自动保存状态没刷新第二次循环可能沿用旧参数。建议循环开头先sim一下小步长测试同步。如果你有 Parallel Computing Toolbox还可以把for改成parfor。不过并行跑 Simulink 模型时要注意每个 worker 必须能访问到当前模型通常需要先将模型保存到当前路径并且每个 worker 都要能打开该模型。并行仿真的优势在模型本身很大、单次仿真时间超过几十秒时才会明显显现。批量运行结束后除了表格还可以画一张“故障发生时间对冲击电流的影响”曲线figure; plot(fault_times, peak_current, o-, LineWidth, 1.5); xlabel(故障发生时刻 / s); ylabel(A 相峰值电流 / A); title(不同故障时刻下三相短路电流峰值); grid on;8. 常见问题与排查方法问题现象可能原因排查方式解决方案找不到 Three-Phase Fault 模块没有安装 Simscape Electrical 工具箱ver查看工具箱列表重装 MATLAB 并勾选 Simscape Electrical双击 Powergui 报错或无法设置MATLAB 版本太旧工具箱冲突查看错误日志which powergui升级到 R2020a 及以上仿真一直卡住不完成求解器选择不当或步长上限过大查看仿真时间、尝试减小最大步长切换到 ode23tb / ode15s最大步长设置为 0.001短路瞬间波形毛刺严重故障电阻为 0 或求解器处理突变能力有限检查 Fault resistance 是否为 0将故障电阻改为 0.001 ΩScope 中波形全部为一条线或异常信号类型不匹配或数据未更新查看 Scope 连接信号是否有离散/连续混合使用 V-I Measurement 输出信号检查采样时间MATLAB R2022b 启动报 Error 9安装目录权限异常或许可证服务未正常启动Windows 事件查看器查找 MATLAB 日志以管理员身份重新激活或修复安装To Workspace 导出后数据维度不对Save format 设置的不是 Timeseries双击 To Workspace 查看 Save format建议设为 Timeseries列号以实际输出为准批量仿真第一次跑通第二次报错模型状态污染或临时文件残留循环开始前执行bdclose(all)在批量循环前重新打开模型保存一次修改故障类型后波形没有变化参数没有生效或模型未保存检查 set_param 修改的模块路径是否正确先手动双击故障模块确认参数已变针对“MATLAB R2022b error 9 错误”这类启动或运行级问题多数情况与 MATLAB 本身的安装、许可证激活、路径权限有关与 Simulink 模型关系不大。遇到这类报错时先看命令行或 Windows 事件查看器里的完整错误文本确认是不是工具箱许可证检查失败。若确认是许可证问题在 MathWorks 账户中心重新激活许可证通常能解决。9. 最佳实践与使用建议把模型搭通只是第一步真正让仿真结果可用、可复现、可用于论文或报告需要建立一套规范化的使用习惯。第一运行顺序固定为“先稳态、再故障、再恢复”。不要一上来就把故障模块设置在 0 时刻因为仿真起始瞬间本身就有系统初始化的暂态过程和短路暂态混在一起会很难分析。建议仿真时间设置 1 秒故障放在 0.2 秒之后留给系统足够的稳定时间。第二每个参数都采用有依据的数值。比如线路阻抗的 X/R 比值是决定短路电流衰减速率的关键不要随手填一个 1 mH 的电感建议先查教材或实际工程参考值明确自己的模型是研究“大电抗系统”还是“低阻抗系统”。第三模型命名和输出变量做到可读。建议模型名采用three_phase_fault_demo这类带下划线的风格。To Workspace的变量名不要用simout、out1这种默认名建议直接使用Is_abc、Vs_abc这类带物理含义的名字批量脚本会好写很多。第四所有数据、波形图、脚本、模型文件分目录保存。可以用下面的结构组织project/ ├── model/ │ └── three_phase_fault_demo.slx ├── script/ │ ├── run_fault_test.m │ └── plot_results.m ├── data/ │ └── sim_20250101.mat └── figure/ └── fault_waveform.png第五涉及数据安全和合规不要在任何公开场合直接贴出未经脱敏的实际电网拓扑和运行参数。即使论文需要也要做等效处理例如把实际变电站参数转换为标幺值或统一基准值下的模型参数既不影响分析趋势也保护了数据边界。10. 总结与下一步这个三相短路 Simulink 仿真模型最值得上手的地方在于它用很少的模块就完整呈现了电力系统暂态分析中的核心现象短路冲击、非周期分量衰减、故障切除后的恢复过程。你最先要验证的不是复杂的参数扫描而是跑通“稳态—短路—切除”三段波形能正确读出故障期间最大电流。最容易踩的坑有两个一是没有安装 Simscape Electrical 导致找不到故障模块二是求解器选错导致仿真时间过长或波形毛刺明显这两点都能通过本文第 3 章和第 8 章的检查表快速解决。后续扩展方向很明确把故障类型改成单相接地和两相短路对比不同短路类型下的故障相与非故障相电流增加一个 Three-Phase Breaker 模块配合保护逻辑模拟保护动作切除故障线路接入三相变压器模型分析短路电流穿越变压器后的变化甚至可以用 Simulink Coder 把保护逻辑导出为嵌入式代码与外部硬件联动实验。对整个电力系统仿真体系来说三相短路模型不是最终目标它更像是你进入暂态分析世界的第一块敲门砖。建议把当前模型和脚本保存好后续做故障分析和保护整定研究时稍作改造就能继续用。
电力系统三相短路Simulink仿真模型搭建与波形分析
这次我们来看电力系统三相短路的 MATLAB Simulink 仿真怎么搭。短路暂态分析是继电保护整定、开关设备选型、电力系统稳定研究的基础而三相短路又是最严重、最常被拿来做首轮研究的故障形式。很多做毕设、课程设计或者刚接触电力系统仿真的人都会卡在同一个问题上Simulink 模型到底怎么建故障模块在哪里参数该怎么填仿真出来的波形怎么解读。这篇文章会从零开始带你在 Simulink 中完成一个包含三相电源、输电线路、短路故障模块、测量环节的暂态仿真模型然后演示怎么样把短路电流波形录出来、读取冲击电流值、对比不同故障类型的结果并给出一个用 MATLAB 脚本批量扫描故障时刻的示例方便你做参数分析和后续扩展。全文主要面向电气工程学生、继电保护方向研究生、刚转电力系统仿真的一线工程师以及任何需要在 MATLAB 里快速验证短路暂态过程的人。不需要先具备 Simulink 高级建模经验能看懂常用模块和基础连线即可但需要有 Simulink 与 Simscape Electrical 相关工具箱的 MATLAB 环境。1. 核心能力速览能力项说明项目类型电力系统暂态仿真模型搭建与验证使用工具MATLAB Simulink Simscape ElectricalSpecialized Power Systems主要功能三相短路暂态过程仿真、短路电流波形记录、冲击电流读取、不同类型短路对比模型复杂度中等适合教学、课程设计、科研预研硬件要求普通台式机 / 笔记本即可CPU 为主内存建议 8GB 以上工具箱要求Simulink、Simscape、Simscape Electrical需安装完整许可输出效果Scope 波形、To Workspace 数据、MATLAB 后处理图表接口与批量可用 sim 命令批量仿真可通过 set_param 修改参数实现参数扫描适合场景继电保护整定分析、暂态特性教学、开关设备选型前期研究2. 适用场景与使用边界这个仿真模型解决的核心问题是在没有实际实验设备的情况下用数字仿真还原三相短路发生时电力系统的暂态行为。它适合做继电保护原理验证、短路电流衰减特性分析、不同中性点接地方式对短路电流的影响、以及保护装置的动模试验前校验。不过要特别注意使用边界。Simulink 里用集中参数线路和理想故障模块建立的模型不等同于实际电网。它忽略了很多工程细节比如线路的分布参数特性、电弧动态特性、变压器饱和、实测故障电阻变化等。因此仿真结果可以用来做相对趋势分析但不能直接替代实际系统的短路计算更不能直接把仿真得到的短路电流数值拿去做现场保护定值整定。涉及真实电力系统数据时参数应脱敏处理并且最终结果必须由专业工程师结合实测和业界标准计算工具校验。在法律合规方面使用 MATLAB Simulink 做常规电力系统仿真不涉及敏感内容但如果你的模型基于实际变电站或电网数据请务必确保数据来源合法不要公开未脱敏的网络拓扑、设备铭牌和运行参数。仿真模型本身可以用来做保护逻辑测试但严禁将模型输出未经校验地用于真实设备投运决策。3. 环境准备与前置条件开始搭建之前先确认你的 MATLAB 环境。比较稳妥的配置是 R2020a 及以上版本因为 Simscape Electrical 在较新版本中刀具更完整模块搜索也更方便。如果你还在用 R2018 或者更早的版本也能做但部分模块路径会不同需要手动查找。打开 MATLAB 后在命令行窗口输入ver查看已安装的工具箱列表中是否有下列内容SimulinkSimscapeSimscape ElectricalMATLAB Report Generator非必需但有助于导出仿真报告如果没有 Simscape Electrical最直接的解法是重新运行 MATLAB 安装程序勾选“Simscape Electrical”组件或者联系管理员开通许可证。注意这个工具箱早期叫 SimPowerSystems后来并入 Simscape所以你在较老教程里看到“SimPowerSystems”字样指的就是同一套电力系统仿真库。环境检查建议清单检查项要求校验方式MATLAB 版本R2020a 及以上越低版本兼容性越差version命令Simscape Electrical必须安装ver查看工具箱列表内存建议 8GB 以上模型不大时 4GB 可运行系统任务管理器磁盘空间至少预留 5GBMATLAB 运行时需要缓存系统磁盘属性路径权限模型文件和脚本所在目录可写新项目放在独立文件夹中4. Simulink 三相短路仿真模型搭建下面进入核心部分。目标是建立一个10kV三相电源经过一段输电线路后发生三相短路的模型故障由Three-Phase Fault模块控制波形用Scope显示同时用To Workspace把电流数据导出到 MATLAB 工作区方便后续计算冲击电流和批量扫描。4.1 新建模型并添加模块在 Simulink 起始页点击“Blank Model”新建模型然后打开库浏览器Library Browser搜索并添加以下模块模块名称数量作用Three-Phase Source1提供三相交流电源Three-Phase Series RLC Branch1模拟输电线路阻抗Three-Phase Fault1控制三相短路发生和切除Three-Phase V-I Measurement1测量母线三相电压和电流Scope2显示电压和电流波形To Workspace2把电压电流导出到工作区Powergui1电力系统模型求解器设置与初始化把模块拖到画布后按图 1 的思路连线Three-Phase Source 输出端接到 Three-Phase V-I Measurement 的输入侧V-I Measurement 的输出端串联 Three-Phase Series RLC Branch线路输出端接到 Three-Phase Fault 模块的输入端和输出端电测量模块的测量输出端接 Scope 和 To Workspace。注意交流电源、线路、故障模块之间要使用 Simscape Electrical 电气连接端口波形显示和 To Workspace 需要先经过信号转换节点一般是 V-I Measurement 输出的信号端口。4.2 电源与线路参数配置双击Three-Phase Source按下表思路设置Phase-to-phase RMS voltage10e3VFrequency50HzInternal connectionYg星形接地Source short-circuit level可以作为系统等值阻抗的简化体现不设置时按理想电压源处理双击Three-Phase Series RLC Branch时需要先确认使用方式是串联 RLCResistance R1ΩInductance L10e-3HCapacitance Cinf线路对地电容暂不考虑时设为 inf相当于开路这几个参数决定了短路发生后的时间常数。线路的 X/R 比值直接影响非周期分量的衰减速率X/R 越大短路电流直流分量衰减越慢。你后续做批量分析时可以替换成不同 L/R 组合观察波形差异。4.3 故障模块配置双击Three-Phase Fault关键参数直接决定仿真行为参数名建议值说明Fault typeABC 或 A-G、A-B 等选择短路类型三相短路选 ABCTransition Times[0.2 0.3]第一个是故障开始时刻第二个是故障切除时刻Fault resistance0.001Ω故障电阻越小短路电流越接近金属性短路Ground resistance0.001Ω接地故障时的接地电阻三相短路也可保留默认这里0.001欧姆的故障电阻既能模拟接近金属性短路的特性又能避免纯零电阻导致部分求解器出现数值振荡。注意故障模块的 Transition Times 是两个时间点格式是类似[0.2 0.3]的数组表示故障从 0.2 秒发生、0.3 秒切除。如果只填一个时间点表示故障在该时刻发生并永久保持。4.4 测量与波形显示Three-Phase V-I Measurement的测量端输出一组电压和电流信号直接接到两个 Scope 上就能看到三相电压和三相电流波形。为了让数据进入 MATLAB 工作区同时接两个To Workspace模块Variable name分别设为Vs_abc和Is_abcSave format选TimeseriesSample time设为0继承输入信号步长这里建议优先用 Timeseries 格式因为后续在 MATLAB 脚本里处理起来最方便可以直接拿到时间和数据值。4.5 Powergui 求解器设置电力系统模型里必须加入Powergui模块否则仿真会报错或结果异常。双击 Powergui 进入设置Simulation type 选ContinuousSolver type 选Variable-stepSolver 选ode23tb或ode15s电力系统短路暂态属于典型的刚性系统短路瞬间电流变化非常快非周期分量缓慢衰减这两段时间尺度差异很大。ode23tb和ode15s都是 MATLAB 里处理刚性问题的常见求解器能有效避免步长过小导致仿真卡死。5. 功能测试与效果验证模型搭好后不要直接一上来就跑短路工况。更合理的顺序是先验证无故障稳态再叠加短路这样才能判断波形异常到底是模型本身问题还是故障引入的。5.1 无故障稳态校验先把Three-Phase Fault模块的 Transition Times 临时改成一个远超仿真结束时间的值比如[999 1000]相当于整个仿真过程中故障都不发生。设置仿真停止时间为0.2秒点击运行。预期结果是三相电压和三相电流均为标准正弦波三相电流幅值由线路阻抗决定幅值基本恒定。这一步验证了电源、线路、测量链路是通的。如果波形畸变明显或者电流幅值异常大先回去检查线路阻抗与电源电压是否匹配比如线电压 10kV 时1Ω 线路电阻加 10mH 电感下稳态电流应该在几百安到上千安级别如果看到上百万安培的电流大概率是电源内阻没设置并且线路阻抗设成了接近 0。5.2 三相短路暂态波形恢复故障模块 Transition Times 为[0.2 0.3]仿真停止时间设为1秒。再次运行。这里重点看短路点附近的电流也就是Is_abc。预期你可以看到0.2 秒前三相电流正常正弦。0.2 秒时刻故障相电流突然激增出现很大的冲击峰值。第一个峰值后电流不是立刻变成对称正弦而是叠加了一个按指数规律衰减的直流分量导致波形的正负半轴不对称整体包络线逐渐回落到稳态短路电流水平。0.3 秒故障切除后电流快速恢复或过渡到新的稳态。这就是典型的短路暂态过程。观察非周期分量衰减快慢时可以对比线路电阻 R 和电感 L 的不同组合R 越大衰减越快反映到波形上就是不对称很快消失R 很小、L 很大时直流分量能维持好几个周波。5.3 用 MATLAB 脚本读取冲击电流波形只能靠“看”要量化分析就需要把数据拿出来算。如果你在模型中设置了To Workspace运行后工作区里已经有Is_abc这个 Timeseries 对象。在命令行窗口执行下面的示例脚本即可提取故障期间的最大电流绝对值也就是近似冲击电流严格意义上冲击电流还跟衰减时间常数有关但作为仿真学习足够% 读取短路电流数据 t Is_abc.Time; Ia Is_abc.Data(:, 1); % A 相电流列号视电压电流顺序调整 Ib Is_abc.Data(:, 2); Ic Is_abc.Data(:, 3); % 提取故障时间段 fault_start 0.2; fault_end 0.3; idx t fault_start t fault_end; % 计算三相电流峰值 peak_a max(abs(Ia(idx))); peak_b max(abs(Ib(idx))); peak_c max(abs(Ic(idx))); disp([A 相故障期间最大电流绝对值: , num2str(peak_a), A]); disp([B 相故障期间最大电流绝对值: , num2str(peak_b), A]); disp([C 相故障期间最大电流绝对值: , num2str(peak_c), A]);注意Is_abc.Data(:, 1)的具体索引取决于 V-I Measurement 输出的信号顺序和 To Workspace 模块的接线顺序。如果运行后报数组越界直接在命令行执行size(Is_abc.Data)查看实际维度再调整列号。把算出来的峰值除以故障前的稳态电流幅值就能得到冲击倍数。电力系统教材里常提的这个倍数通常与系统 X/R 比值有关你的仿真结果应该能直观反映这种关联。5.4 不同短路类型对比三相短路数值最大但实际电力系统中单相接地短路发生概率更高。回头双击Three-Phase Fault把 Fault type 改成A-GA 相接地再仿真对比 A 相电流和三相短路时的差异。如果改成A-B两相短路可以看到非故障相 C 相电流明显小于故障相并且故障相的相位关系也不一样。这个对比实验对理解对称短路和不对称短路的暂态差异非常有帮助是理解后续基于序分量保护原理的重要基础。6. MATLAB Simulink 求解器与性能观察Simulink 仿真电力系统时性能表现跟求解器选择关系很大。很多初学者遇到仿真卡住或波形毛刺多第一反应是模型错了但实际上问题往往出在求解器配置。求解器选择的通用原则场景推荐求解器说明普通连续系统ode45速度快精度可接受电力电子 / 短路暂态ode23tb 或 ode15s刚性系统能承受故障瞬间的快速变化离散系统 / 实时仿真固定步长离散求解器需要先通过 Powergui 进行离散化设置单纯看稳态结果加大步长 ode45可以提高计算速度但暂态细节会丢失如果短路瞬间波形产生明显的高频振荡毛刺优先检查仿真步长上限。在 Configuration Parameters 里最大步长不要设得太大理论上建议不大于系统最小时间常数的四分之一。10mH 电感和 1Ω 电阻对应的电气时间常数为 10ms最大步长建议设在 0.001 秒级别太大会把衰减过程直接抹掉。仿真速度的观察方法也很简单。在 MATLAB 命令行运行脚本时用tic和toc包裹sim命令tic; out sim(mdl, StopTime, 1); toc;如果一次 1 秒仿真时间要跑几十秒甚至几分钟需要看模型是否有大量Scope在实时更新画面。Scope 实时绘制对 CPU 开销很大做批量仿真时不要开着好几个 Scope 窗口跑改用To Workspace记录数据仿真结束后统一画图。另外一个隐藏开销是 Powergui 里如果开启了FFT Analysis或Impedance Measurement仿真前必须关闭或断开否则每次积分都要消耗额外资源。内存占用方面常规的小规模三相短路模型不会把机器内存吃满。真正吃内存的是你把仿真数据全部存成 Timeseries并且仿真时间很长、步长很小时返回的数组可能非常庞大。实际编程时建议只输出关键通道或者在 To Workspace 里设置降采样 Sample time例如设为0.0001秒而不是0。7. 批量仿真与数据导出单次仿真只能得到一个工况的结果做课程设计或写论文时通常需要分析不同故障时刻、不同故障电阻、不同线路阻抗对短路暂态的影响。手动改参数再点运行效率低而且容易漏掉组合。正确的做法是写一个 MATLAB 脚本用set_param修改模型参数再用sim命令循环运行模型。以批量扫描故障发生时刻为例。假设故障模块名是Three-Phase Fault模型名是three_phase_fault_demo脚本可以这样写% 选择模型 mdl three_phase_fault_demo; open_system(mdl); % 需要扫描的故障发生时刻 fault_times [0.05 0.1 0.2 0.3 0.5]; fault_duration 0.1; peak_current zeros(size(fault_times)); for k 1:numel(fault_times) t_start fault_times(k); t_end t_start fault_duration; % 修改故障模块的切换时间 set_param([mdl /Three-Phase Fault], ... SwitchTimes, sprintf([%g %g], t_start, t_end)); % 运行仿真 simOut sim(mdl, StopTime, 1); % 提取 A 相电流并计算故障期间峰值 Is_abc simOut.get(Is_abc); t Is_abc.Time; Ia Is_abc.Data(:, 1); idx t t_start t t_end; peak_current(k) max(abs(Ia(idx))); end % 汇总结果 results_table table(fault_times, peak_current, ... VariableNames, {FaultTime, PeakCurrent}); disp(results_table);这段脚本里simOut.get(Is_abc)的变量名必须和模型中To Workspace设置的变量名一致。如果你没有勾选模型配置里的“Single simulation output”选项sim命令返回的结果结构会有一点差别最稳妥的情况是调用前用get_param(mdl, SaveOutput)查看。报错时在命令行执行help sim确认当前 MATLAB 版本的输出格式。另外一个容易踩的坑是set_param修改参数后如果模型处于运行状态或者自动保存状态没刷新第二次循环可能沿用旧参数。建议循环开头先sim一下小步长测试同步。如果你有 Parallel Computing Toolbox还可以把for改成parfor。不过并行跑 Simulink 模型时要注意每个 worker 必须能访问到当前模型通常需要先将模型保存到当前路径并且每个 worker 都要能打开该模型。并行仿真的优势在模型本身很大、单次仿真时间超过几十秒时才会明显显现。批量运行结束后除了表格还可以画一张“故障发生时间对冲击电流的影响”曲线figure; plot(fault_times, peak_current, o-, LineWidth, 1.5); xlabel(故障发生时刻 / s); ylabel(A 相峰值电流 / A); title(不同故障时刻下三相短路电流峰值); grid on;8. 常见问题与排查方法问题现象可能原因排查方式解决方案找不到 Three-Phase Fault 模块没有安装 Simscape Electrical 工具箱ver查看工具箱列表重装 MATLAB 并勾选 Simscape Electrical双击 Powergui 报错或无法设置MATLAB 版本太旧工具箱冲突查看错误日志which powergui升级到 R2020a 及以上仿真一直卡住不完成求解器选择不当或步长上限过大查看仿真时间、尝试减小最大步长切换到 ode23tb / ode15s最大步长设置为 0.001短路瞬间波形毛刺严重故障电阻为 0 或求解器处理突变能力有限检查 Fault resistance 是否为 0将故障电阻改为 0.001 ΩScope 中波形全部为一条线或异常信号类型不匹配或数据未更新查看 Scope 连接信号是否有离散/连续混合使用 V-I Measurement 输出信号检查采样时间MATLAB R2022b 启动报 Error 9安装目录权限异常或许可证服务未正常启动Windows 事件查看器查找 MATLAB 日志以管理员身份重新激活或修复安装To Workspace 导出后数据维度不对Save format 设置的不是 Timeseries双击 To Workspace 查看 Save format建议设为 Timeseries列号以实际输出为准批量仿真第一次跑通第二次报错模型状态污染或临时文件残留循环开始前执行bdclose(all)在批量循环前重新打开模型保存一次修改故障类型后波形没有变化参数没有生效或模型未保存检查 set_param 修改的模块路径是否正确先手动双击故障模块确认参数已变针对“MATLAB R2022b error 9 错误”这类启动或运行级问题多数情况与 MATLAB 本身的安装、许可证激活、路径权限有关与 Simulink 模型关系不大。遇到这类报错时先看命令行或 Windows 事件查看器里的完整错误文本确认是不是工具箱许可证检查失败。若确认是许可证问题在 MathWorks 账户中心重新激活许可证通常能解决。9. 最佳实践与使用建议把模型搭通只是第一步真正让仿真结果可用、可复现、可用于论文或报告需要建立一套规范化的使用习惯。第一运行顺序固定为“先稳态、再故障、再恢复”。不要一上来就把故障模块设置在 0 时刻因为仿真起始瞬间本身就有系统初始化的暂态过程和短路暂态混在一起会很难分析。建议仿真时间设置 1 秒故障放在 0.2 秒之后留给系统足够的稳定时间。第二每个参数都采用有依据的数值。比如线路阻抗的 X/R 比值是决定短路电流衰减速率的关键不要随手填一个 1 mH 的电感建议先查教材或实际工程参考值明确自己的模型是研究“大电抗系统”还是“低阻抗系统”。第三模型命名和输出变量做到可读。建议模型名采用three_phase_fault_demo这类带下划线的风格。To Workspace的变量名不要用simout、out1这种默认名建议直接使用Is_abc、Vs_abc这类带物理含义的名字批量脚本会好写很多。第四所有数据、波形图、脚本、模型文件分目录保存。可以用下面的结构组织project/ ├── model/ │ └── three_phase_fault_demo.slx ├── script/ │ ├── run_fault_test.m │ └── plot_results.m ├── data/ │ └── sim_20250101.mat └── figure/ └── fault_waveform.png第五涉及数据安全和合规不要在任何公开场合直接贴出未经脱敏的实际电网拓扑和运行参数。即使论文需要也要做等效处理例如把实际变电站参数转换为标幺值或统一基准值下的模型参数既不影响分析趋势也保护了数据边界。10. 总结与下一步这个三相短路 Simulink 仿真模型最值得上手的地方在于它用很少的模块就完整呈现了电力系统暂态分析中的核心现象短路冲击、非周期分量衰减、故障切除后的恢复过程。你最先要验证的不是复杂的参数扫描而是跑通“稳态—短路—切除”三段波形能正确读出故障期间最大电流。最容易踩的坑有两个一是没有安装 Simscape Electrical 导致找不到故障模块二是求解器选错导致仿真时间过长或波形毛刺明显这两点都能通过本文第 3 章和第 8 章的检查表快速解决。后续扩展方向很明确把故障类型改成单相接地和两相短路对比不同短路类型下的故障相与非故障相电流增加一个 Three-Phase Breaker 模块配合保护逻辑模拟保护动作切除故障线路接入三相变压器模型分析短路电流穿越变压器后的变化甚至可以用 Simulink Coder 把保护逻辑导出为嵌入式代码与外部硬件联动实验。对整个电力系统仿真体系来说三相短路模型不是最终目标它更像是你进入暂态分析世界的第一块敲门砖。建议把当前模型和脚本保存好后续做故障分析和保护整定研究时稍作改造就能继续用。
