简介基于MATLAB的三相短路分析资料包适合电力系统专业学生与工程技术人员用于掌握故障建模、仿真与结果分析的方法。资源从三相短路基本概念出发结合Simulink平台完整演示了电力系统模型搭建、故障类型与位置设置、仿真运行及数据解读的全过程并给出可复用的test.m脚本和详细课设报告。压缩包共12个文件大小仅1.09MB包含m脚本、PDF报告、4个txt数据文件与6张PNG结果图txt文件保存了电路原始数据、节点阻抗/导纳矩阵及输出结果PNG图清晰呈现发电机等效前后的矩阵变化便于对照报告逐项理解。目前已有2576人学习尤其适合正在完成电力系统课设或MATLAB仿真作业的读者。通过研读代码与报告可以快速掌握用Simulink分析三相短路电流、验证保护设备动作特性并能直接修改参数扩展至其他故障类型是一份实用的电力仿真学习资料。 做短路计算校验的时候三相短路几乎是每个电力系统从业者都躲不开的基本功。我最早接触三相短路分析是在学暂态过程的时候当时全靠手算公式、查冲击系数曲线算一条辐射线路的短路电流要折腾大半天后来真正参与设备选型和保护整定发现短路电流校核几乎贯穿所有项目才开始系统性用MATLAB来搭仿真模型。MATLAB里的Simscape Electrical模块老版本叫SimPowerSystems支持直接拖拽电源、变压器、线路和故障模块设置一个故障时间就能看到三相短路后的电压电流波形比手算直观太多。这篇把整套流程和踩过的坑整理出来给准备用MATLAB做三相短路分析的同学和工程师一个可以照做的参考。1. 三相短路分析的本质算的是暂态电流不是稳态电流1.1 为什么会短路短路电流为什么可怕三相短路是指供电系统中三相之间发生金属性或非金属性的直接连接相当于负载被短路掉回路阻抗骤减系统电压被瞬间拉低短路电流可以高达额定电流的数倍甚至数十倍。为什么都拿三相短路作为最严重的故障工况来校核因为对称短路时正序阻抗最小短路电流通常最大设备动稳定、热稳定、断路器遮断容量、保护定值整定几乎都要围绕这个最大工况来做。真正让三相短路分析变得棘手的是暂态过程。短路瞬间电流不能突变为了维持磁链守恒回路里会感应出一个非周期分量电流这个直流分量和周期分量叠加导致短路电流波形严重偏置。用公式来表达故障相的短路电流大约为i(t) sqrt(2) I * [sin(ωt α - φ) - sin(α - φ) * e^(-t/Ta)]其中I是次暂态电流有效值α是故障瞬间的电压相位角φ是短路回路阻抗角Ta L/R为直流分量衰减时间常数。初始的冲击电流峰值在半个周波附近出现50Hz系统约0.01s其幅值可以近似为 Kimp * sqrt(2) * IKimp就是冲击系数取值通常在1.8左右。如果不理解这个暂态过程仿真跑出来的波形就只是“看到了一个波形”没办法从里面提取设备校核需要的关键参数。1.2 MATLAB做短路分析的两条路线用MATLAB做三相短路分析我试过两条路线各有适用场景。第一条是纯脚本编程把系统等效成R-L阻抗用解析公式直接算短路电流周期分量、非周期分量和冲击电流。好处是快速、透明适合简单辐射网和批量参数扫描比如我要比较不同线路长度对短路电流的影响写一个循环就出来了。缺点是处理复杂拓扑和多节点网络时手写节点阻抗矩阵非常费劲很容易出错。第二条是用Simulink里的Simscape Electrical模块搭建数学模型这也是我后期项目里用得最多的一条路。它把电源、变压器、线路、负载、故障模块做成可视化元件搭完模型就能看到各点电压电流波形还能方便地切换故障类型、故障时刻和故障位置。缺点是模型自由度多求解器、步长、阻抗参数设置不当结果可能歪到离谱。我个人的建议是做教学演示或验证理论公式优先用脚本做工程系统分析、保护整定或需要给同事看波形时用Simulink建模。两者结合是最好的先用脚本算一个基准值再拿Simulink去验证和扩展。2. 建模前的参数整理从单线图到Simulink元件参数2.1 系统阻抗与标幺值换算不管用什么工具三相短路分析的输入基础都是系统阻抗。模型里的每个元件都要有明确的R和X这些参数可以从设备铭牌、线路单位阻抗表和系统短路容量中获取。Simulink里可以直接填欧姆Ω和亨H但我建议在做理论校核时先换算成标幺值这样对比手算结果更方便。基准值的选取遵循 S_B 100 MVA或10 MVA基准电压取电网额定电压例如10kV系统选U_B 10.5kV基准阻抗为 Z_B U_B^2 / S_B。举个例子我常用的一个验证模型是10kV辐射网系统电源侧短路容量为200MVA折算到10.5kV基准下的电源等效阻抗约 Z_source U_B^2 / S_k 10.5^2 / 200 0.551Ω。变压器容量2MVA短路电压百分比 U_k% 6%折算到10.5kV侧变压器阻抗约为 Z_T U_k% * U_B^2 / (100 * S_T) 6 * 10.5^2 / (100 * 2) 3.307Ω。线路用10km的10kV架空线每公里电抗0.35Ω电阻0.17Ω线路阻抗 Z_L 1.7 j3.5Ω。这些值最后填到Simulink元件里就是有名值。2.2 在Simscape Electrical库中选型与搭建新版MATLAB的Simulink库浏览器里路径一般是Simscape Electrical Specialized Power Systems Fundamental Blocks。老版本直接叫SimPowerSystems位置类似。需要用到的基本元件有Three-Phase Source三相电压源设置相电压/线电压有效值、频率、内阻和电感。Three-Phase Transformer (Two Windings)双绕组三相变压器设置额定容量、频率、一次二次额定电压、短路阻抗百分比。Three-Phase Series RLC Branch线路阻抗可以直接填R、L。Three-Phase Fault三相故障模块可选故障类型和切换时间。Three-Phase V-I Measurement测量三相电压电流输出到Scope或Workspace。Powergui电力系统图形用户界面模块选择连续模式或相量模式。我一般选Continuous才能看到暂态波形。搭建顺序不复杂三相电压源接到变压器一次侧变压器二次侧经过线路阻抗到故障模块再接到负载。故障模块放在短路点测量模块放在电源侧和故障点两侧各一个方便对比电压跌落和电流变化。别忘了在模型中放置Powergui模块否则整个模型会报错。把上节计算出的参数填入元件对话框Three-Phase Source的Phase-to-phase RMS voltage填10.5kV如果电源侧视为无穷大系统内阻为0但为了仿真接近实际建议填0.5Ω和0.01H变压器一次10.5kV、二次0.4kV也可以不降压用于等值短路计算线路RLC填1.7Ω和0.35H对应10km注意电抗换算到电感需要除以2πf。2.3 故障模块参数设置细节Three-Phase Fault模块是所有仿真里的核心开关。双击模块Fault type选Three-phase faultABC这时三相同时短路。Fault resistance填0.01Ω不要填0数值上可以更稳定如果做金属性短路填一个很小的值即可。切换方式我习惯用外部控制External control用Step模块从0跳到1触发故障模块导通也可以用内部定时直接设置Transition time为0.1s故障初始状态选Open经过0.1s后切换到Closed。这里有个容易被忽略的点Transition time是故障开始时间不是故障持续时间。如果希望故障持续0.2s后恢复需要设置两个Transition time比如0.1和0.3s状态依次为Open→Closed→Open。只填一个时间故障会一直保持到仿真结束。对于短路分析我们通常只关心故障后前几个周波所以故障持续0.2s足够。我还习惯在故障模块的External Fault Breaker控制端口接一个Signal Generator或者Step这样可以精确控制故障触发相位。想研究短路初始角的影响就通过Step的延迟时间调整短路发生的合闸角。3. 仿真运行与波形解读短路的真正信息藏在暂态过程里3.1 从波形中读出周期分量和非周期分量模型搭好设置仿真时长0.4s最大步长0.1ms运行后双击Scope看到故障点电流波形。稳定情况下看到故障前是规则正弦故障发生后电流跳变前半段波形上下严重不对称包络线一高一低随后逐渐恢复到对称正弦。这个不对称就是非周期分量在衰减而对称部分就是短路电流的周期分量。我在实际读波形时会把数据导到MATLAB工作区再画图方便标注峰值。比如用To Workspace模块输出三相电流然后在命令窗口提取数据t out.tout; iA out.iabc(:,1); plot(t, iA); xline(0.1, --); grid on;如果仿真的故障开始时间是0.1s波形图里0.1s之前的正弦就是负荷电流或空载电流0.1s之后才是短路电流。用数据游标读出故障后第一个周期波峰大约15kA再往后几个周波峰值回落到约8kA前者包含了直流分量后者才接近周期分量幅值。3.2 冲击电流与衰减时间常数的实测方法冲击电流是短路后第一个峰值是设备动稳定校验的核心参数。我在模型中读取A相电流最大值为15kA左右与理论估算的冲击电流对照。理论值怎么算先由短路回路总阻抗算出次暂态周期分量有效值I然后乘以sqrt(2)再乘以冲击系数Kimp。Kimp可以用1.8的典型值也可以按Kimp 1 e^(-0.01/Ta)计算。T_a L/R (总电抗/2πf)/总电阻。如果L/R算下来是0.03sKimp 1 e^(-0.333) ≈ 1.717理论峰值约为 sqrt(2) * I * 1.717。衰减时间常数还能直接从波形里测出来把故障后的电流包络线提取出来直流分量衰减到初始值的36.8%所需的时间就是Ta。由于非周期分量在示波器里表现为包络线的不对称差可以取上包络线减下包络线再除以2得到直流衰减曲线用对数坐标拟合得出时间常数。这个方法适合验证实际模型中的回路时间常数是否与手算一致。3.3 不同故障时刻对波形的影响短路初始角对冲击电流影响很大。这是因为公式里sin(α-φ)决定了非周期分量的初始幅值。如果短路瞬间电压过零α≈0回路阻抗角φ约80度sin(α-φ)很大直流分量几乎以最大初始值出现冲击电流最大如果短路瞬间电压刚好在峰值附近非周期分量初始值很小波形基本对称冲击系数接近1。用MATLAB模拟这个现象很简单在Three-Phase Source模块里设置Phase angle of phase A或者用外部Step控制故障触发时间把故障时间从0.1s改成0.105s等观察第一个峰值的不同。我实测发现同一个模型只改故障时刻冲击电流能从14kA变到16kA差别超过了10%。所以在做保护整定时不能光看一组仿真结果要找最严重的合闸角。4. 理论手算与仿真结果交叉验证差距从哪来4.1 手算短路电流的简化模型对于辐射网手算三相短路电流的简化步骤是忽略负荷把电源等效为次暂态电势E短路点电压在故障前按额定电压考虑短路电流有效值I_sc c U_N / (sqrt(3) Z_total)。电压系数c在高压系统通常取1.1低压系统取1.05。Z_total是短路点到电源的全部等值阻抗之和包括电源内阻抗、变压器阻抗和线路阻抗。我用前面10kV模型来算Z_source 0.551ΩZ_T 3.307ΩZ_L 1.7 j3.5Ω。要注意电抗和电阻要分别加再求模。总电阻R 0.5 0 1.7 2.2Ω电源电阻很小忽略或按前面填的0.5Ω总电抗X 0.551 3.307 3.5 7.358Ω假设电源内阻电抗0.551。Z_total sqrt(2.2^2 7.358^2) ≈ 7.68Ω。取c 1.1短路电流有效值 I_sc 1.1 * 10.5k / (sqrt(3) * 7.68) ≈ 0.868kA 868A。这里因为模型用了10kV/0.4kV变压器和0.4kV侧故障阻抗折算会不同手算时应注意基准侧选择为了简单这里是在10kV侧等值短路点。把仿真结果的短路电流周期分量有效值从波形里提炼出来两者通常能对上偏差在5%以内。如果偏差超过10%不要急着说仿真不对先去检查参数折算和基准值。4.2 仿真和手算不一致的典型原因我在几次对比中总结出几个最常见的不一致来源变压器模型是否包含励磁支路和饱和非线性。手算基本忽略励磁电流Simulink的变压器模型如果选择饱和饱和特性短路电流会略有偏差但对短路电流影响不大如果变压器模型没有选择模拟饱和则和手算接近。故障电阻。手算假设金属性短路R_fault0仿真为了数值收敛填了0.01Ω多了一个微小的电阻测出来的电流会偏小一点可以忽略。负荷的影响。手算短路电流通常忽略负荷但Simulink模型里如果接了额定负载故障前的电流不是0故障瞬间的初始相位条件会变进而影响非周期分量但周期分量有效值基本不受负荷影响。线路的分布电容。手算用集中阻抗Simulink的分布式参数线路模型Distributed Parameters Line在短路瞬间会有暂态行波看起来波形高频振荡这和简化手算差别很大。对于短线路建议用Series RLC Branch或PI型线路模型波形更干净结果也更接近手算。4.3 用脚本批量对比参数我做参数校核时会用脚本循环修改线路长度把仿真后的周期分量有效值和手算值放在一起画曲线这样能直观看出趋势和偏差。做法是把模型参数用变量L_km替代在仿真脚本里用set_param改参数并运行sim再把结果存到数组。这里提供一段简化脚本思路L_vec 2:2:20; I_sim zeros(size(L_vec)); for i 1:length(L_vec) set_param(short_circuit/Line RLC, R, num2str(0.17*L_vec(i))); set_param(short_circuit/Line RLC, L, num2str(0.35*L_vec(i)/(2*pi*50))); sim(short_circuit); % 提取故障后电流稳态峰值 I_sim(i) max(abs(out.iabc(300:end,1))) / sqrt(2); end这里set_param只改阻抗模型会自动按新参数重新仿真非常方便扫参。注意L的单位是亨利需要把线路电抗除以314。5. 实操中踩过的坑与调试经验5.1 仿真步长与求解器设置三相短路仿真最怕步长太大把冲击电流峰值给抹平了。默认的变步长求解器在快速暂态过程中会自动减步长但也可能不够密。我通常在模型配置参数里把最大步长固定为1e-4s0.1ms求解器选ode23tb或ode15s这两种对电力电子和电磁暂态问题比较友好。如果只关心波形宏观趋势可以放开到0.5ms但读冲击电流时必须用0.1ms左右否则读出的峰值可能偏小5%以上。另一个坑是仿真总时长。故障时间设0.1s仿真总时长至少要0.2s才能看到故障后的暂态衰减过程我一般设0.4s如果还想看故障清除后的恢复过程就设0.6s并在故障模块里配置第二个Transition time。5.2 测量环节的读数陷阱用Three-Phase V-I Measurement输出到工作区后数据是结构体不是数组。我一开始直接plot出错了后来用out.iabc.signals.values提取才行。如果测量模块里选择的是Use a label那么会被命名成变量需要在模型里用Signal Label把标签对应好。Scope里看到的电流是瞬时值最大值是峰值不是有效值。很多新手直接把Scope里的峰值当作短路电流有效值去套公式结果对不上。要取有效值应该把故障后达到稳态的波形做RMS计算或者用峰值除以sqrt(2)。还要注意区分是故障相电流还是线电流三相故障时三相电流等大但如果是单相故障读故障相和非故障相意义完全不同。另外故障模块如果设置了外部控制Step模块从0跳到1之后要确保Step的Final value是1若Final value是0故障模块可能一直在导通或断开状态检查仿真波形有无异常通常一眼能看出来。5.3 模型扩展从三相短路到其他故障类型三相短路分析是基础但实际系统中单相接地、两相短路也很常见。Simulink里切换故障类型很快Three-Phase Fault模块选择Phase A fault、Phase B fault或AB fault即可。如果是接地故障记得把Fault resistance和接地方案配好要观察零序电流和三相不对称波形。故障位置也能随便挪我把故障模块从线路末端移到母线侧、变压器出口、发电机出口测不同位置的短路电流能直观理解短路点越靠近电源阻抗越小短路电流越大。这种扩展对一个校核项目来说很有价值比如断路器能不能装在当前母线位置就看该点的短路电流是否在遮断容量以内。5.4 关于模型收敛的最后提醒有时候模型运行报错显示“Solver encountered an error”十有八九是初始条件或参数突变造成的。常见解决办法把故障电阻从0.01Ω提高到0.1Ω看能否运行把Powergui的Simulation type切到PhasorPhasor模式不看暂态波形但能快速算稳态短路电流把负荷从模型中断开用纯阻抗网络测试。我在调试时习惯把故障模块的Transition time设置成0.1s但先不接入线路直接测空载模型是否稳定再加故障一层层排查能省很多时间。总的来说基于MATLAB做三相短路分析本质上不是学会拖模块而是能把系统阻抗、暂态过程、冲击电流这些物理概念用数字曲线表达出来并反复验证。这个流程跑通之后后面不管换什么系统接线我都是先用脚本估一个数再搭Simulink模型精细看波形最后回到手算公式解释偏差。整个过程里最有价值的一点就是不要盲信仿真波形任何一组数据都要能回到物理原理上讲得通才算真可靠。本文还有配套的精品资源点击获取
尧图网络