简介面向电气工程专业学生、科研人员及电力系统设计维护人员的MATLAB供配电系统短路故障仿真资料重点解决利用SimPowerSystem搭建无穷大功率电源供电系统模型、模拟三相短路并分析暂态过程的问题。压缩包内为1个PDF文件大小116KB内容包含课程设计任务书、技术参数、仿真要求及工作计划便于按步骤学习和打印查阅。已有92人浏览学习。内容围绕110kV无穷大系统在0.02s发生三相对称短路的场景完整给出变压器容量、短路电压、线路电阻电抗等关键元件的参数取值与计算方法并对比无穷大系统与同步发电机突然三相短路时电流变化的异同帮助读者深入理解短路暂态特性。同时覆盖从理论复习、SimPowerSystem建模调试、仿真结果分析到社会、安全、环境等因素思考的完整流程既适合课程设计参考也可作为工程项目实施与教学训练素材。1. 短路故障仿真第一步把无穷大电源的暂态过程读懂做供配电系统课程设计时很多人一上来就打开 Simulink 拖元件结果波形出来了却讲不清每个拐点对应什么物理过程。电力系统短路故障仿真里最容易被忽视、也最值得先搞懂的概念是三相短路瞬间电流最大值并不是故障刚发生的那个时刻而是大约半个周波之后。这个现象由非周期分量直流分量叠加造成是冲击电流计算的物理基础。本文以 MATLAB/SimPowerSystems现为 Simscape Electrical Specialized Power Systems为工具完整走一遍无穷大功率电源供电系统三相对称短路的理论推导、元件参数解析计算、仿真建模到结果分析的全流程适合电气工程专业学生做课程设计、毕业设计也适合刚接触电力系统仿真的工程技术人员对照排查自己的建模方法。2. 无穷大电源三相短路的暂态过程周期分量与非周期分量的分离2.1 无穷大电源的数学模型无穷大功率电源是电力系统短路计算中的理想化概念指电源容量远大于被研究系统、内阻抗为零、短路前后母线电压幅值和频率均保持恒定的电源。实际系统中 110kV 及以上电网的相对容量足够大时可近似按无穷大电源处理。设电源相电压为 $u_a U_m \sin(\omega t \alpha)$其中 $\alpha$ 为短路瞬间 A 相电压的初相角。短路前系统带纯阻感负载A 相电流表达式为$$i_a I_m \sin(\omega t \alpha - \varphi)$$其中 $I_m U_m / \sqrt{R^2 X^2}$$\varphi \arctan(X/R)$。这里 $R$、$X$ 是从电源到短路点之间的等值电阻与电抗。关键点在于短路瞬间电感中的磁链不能突变但短路后的电路阻抗突然变小稳态电流幅值会跳变到另一个值。为了满足磁链守恒约束回路中必然出现一个自由分量电流这个分量按指数规律衰减其初始值恰好抵消周期分量在短路瞬间的突变部分。2.2 短路电流的完整表达式与冲击电流对三相短路取 A 相分析短路后的电流可以写成强制分量与自由分量之和$$i_a I_{pm} \sin(\omega t \alpha - \varphi_k) \left[ I_m \sin(\alpha - \varphi) - I_{pm} \sin(\alpha - \varphi_k) \right] e^{-t/T_a}$$式中 $I_{pm}$ 是短路电流周期分量的幅值$\varphi_k$ 是短路回路阻抗角$T_a L_k / R_k$ 是短路回路非周期分量衰减时间常数。方括号里那一项就是非周期分量也叫直流分量它的大小取决于短路时刻的电压初相角 $\alpha$。当 $\alpha - \varphi \approx \pi/2$ 时非周期分量初始值最大合成电流在最严重情况下会在短路后约半个周期50Hz 系统为 0.01s达到峰值这个峰值称为冲击电流$$i_{sh} I_{pm} I_{pm} e^{-0.01/T_a} K_{sh} \sqrt{2} I$$$K_{sh} 1 e^{-0.01/T_a}$ 为冲击系数工程上高压系统一般取 1.61.9。这一结论直接指导仿真中的故障时间设置0.02s 设置短路则波形约在 0.028~0.03s 处出现第一个最高峰分析波形时应对照这个时刻验证模型是否正确。2.3 为什么 0.02s 是一个合理的短路时刻设定课程设计任务书中要求 0.02s 发生三相对称短路这个时间不是随便给的。电力系统继电保护的主保护动作时间通常在 2040ms 量级0.02s 短路对应的是保护尚未切除故障的最严重工况。仿真时还要设置一个故障切除时间比如 0.3s让故障模块在分析完暂态全过程后再断开目的是观察短路电流从暂态恢复到稳态的完整曲线。另外要注意任务书把短路点设在变压器低压侧出口。这意味着短路回路中包含线路阻抗和变压器短路阻抗两部分电源电压等级为 110kV短路点电压等级为 11kV做理论计算时必须用标幺值法或折算到同一电压等级后进行否则阻抗直接相加会出现数量级错误。3. 元件参数解析计算从变压器铭牌到标幺值阻抗3.1 基准值选取与线路参数计算理论计算采用标幺值法基准容量取 $S_B 100\text{MVA}$电压基准取各电压等级的平均额定电压110kV 侧取 $U_{B1} 115\text{kV}$11kV 侧取 $U_{B2} 10.5\text{kV}$。这样选取的动机是平均额定电压更接近实际运行电压且各电压等级间通过标准变比自然解耦。线路参数按任务书给定$L 50\text{km}$$x_1 0.4\Omega/\text{km}$$r_1 0.17\Omega/\text{km}$则线路总电抗 $X_L 20\Omega$、总电阻 $R_L 8.5\Omega$。高压侧基准阻抗为$$Z_{B1} \frac{U_{B1}^2}{S_B} \frac{115^2}{100} 132.25\Omega$$线路标幺值$$X_{L}^{} \frac{20}{132.25} 0.1512, \quad R_{L}^{} \frac{8.5}{132.25} 0.0643$$这部分计算的意义在于后续 Simulink 建模中如果直接填有名值需要把线路阻抗换算到 11kV 侧而用标幺值计算短路电流时则完全不受电压等级换算困扰。两套体系各有利弊建议理论分析用标幺值仿真建模直接在线路模块里填 50km、$r_1 0.17\Omega/\text{km}$、$x_1 0.4\Omega/\text{km}$ 这三个原始参数让 SimPowerSystems 自动处理。3.2 变压器短路阻抗与损耗参数折算变压器铭牌参数为 $S_N 20\text{MVA}$$U_k% 10.5$$\Delta P_k 135\text{kW}$$\Delta P_0 22\text{kW}$$I_0% 0.8$变比 $110/11\text{kV}$。先计算以变压器额定容量为基准的短路阻抗标幺值短路电抗标幺值 $X_T^* 10.5% 0.105$短路电阻标幺值 $R_T^* \Delta P_k / S_N 135\text{kW} / 20\text{MVA} 0.00675$短路阻抗模值 $Z_T^* \sqrt{0.105^2 0.00675^2} \approx 0.1052$可见电阻占比极小当统一换算到 $S_B 100\text{MVA}$ 基准时变压器电抗标幺值变为$$X_T^{*} 0.105 \times \frac{100}{20} 0.525$$这个数值远大于线路标幺电抗 0.1512说明在 110kV 供电的短路计算中变压器短路阻抗是决定短路电流大小的主因素。线路电阻和变压器电阻对短路电流幅值影响很小但计算冲击系数 $K_{sh}$ 时不能忽略因为 $T_a X / (\omega R)$ 对电阻敏感。3.3 短路电流理论值手算主要计算过程如下% 供配电系统三相短路电流理论计算 % 基准值 SB 100e6; % 基准容量 100MVA UB1 115e3; % 110kV侧平均额定电压 UB2 10.5e3; % 10.5kV侧平均额定电压 IB2 SB / (sqrt(3) * UB2); % 低压侧基准电流约 5499A % 元件标幺值SB基准 xL_pu (0.4 * 50) / (UB1^2 / SB); % 线路电抗标幺值 0.1512 rL_pu (0.17 * 50) / (UB1^2 / SB); % 线路电阻标幺值 0.0643 xT_pu 0.105 * (SB / 20e6); % 变压器电抗标幺值 0.525 rT_pu (135e3 / 20e6) * (SB / 20e6);% 变压器电阻标幺值 0.003375 X_sum xL_pu xT_pu; % 总电抗标幺值 0.6762 R_sum rL_pu rT_pu; % 总电阻标幺值 0.0677 I_pu 1 / sqrt(X_sum^2 R_sum^2); % 短路电流标幺值 1.4758 Ik I_pu * IB2; % 短路电流有名值约 8114A Ta X_sum / (2 * pi * 50 * R_sum); % 非周期分量时间常数约 0.0318s Ksh 1 exp(-0.01 / Ta); % 冲击系数约 1.73 ish Ksh * sqrt(2) * Ik; % 冲击电流约 19.8kA代码逻辑说明电压基准取 115kV 和 10.5kV 两级短路点基准电流用 10.5kV 侧计算线路阻抗标幺值需要除以高压侧基准阻抗变压器阻抗标幺值需要按容量折算不能直接用铭牌值 0.105最后短路电流标幺值等于 1 除以短路回路总阻抗模值。其中Ksh用 50Hz 的半个周期时间 0.01s 计算反映的是最严重短路时刻的冲击系数。3.4 励磁支路参数对短路计算的影响变压器空载损耗 $\Delta P_0 22\text{kW}$ 和空载电流百分数 $I_0% 0.8$ 用于计算励磁支路参数。励磁电阻 $R_m^* 1 / (\Delta P_0 / S_N) 909$励磁电抗 $X_m^* 1 / \sqrt{(I_0/100)^2 - (\Delta P_0/S_N)^2} \approx 126$。由于这两个标幺值远大于短路回路阻抗短路时流过励磁支路的电流极小因此手算短路电流时可以忽略励磁支路。但 Simulink 建模时需要填入这些参数否则变压器空载损耗特性和仿真启动瞬态与实际不符这也是理论计算与仿真调参过程中容易被忽视的差异点。4. SimPowerSystems 建模实操三相短路故障的参数设置与调试4.1 模型拓扑与元件清单SimPowerSystems 库在较新版本 MATLAB 中的路径为 Simscape Electrical Specialized Power Systems核心求解模块 Powergui 必须拖入模型否则仿真报错。模型拓扑按电源到短路点的功率流向布置三相电源 - 三相串联 RLC 支路模拟 50km 线路- 双绕组三相变压器 - 三相故障模块 - 接地。元件库路径本例参数设置Three-Phase SourceElectrical Sources相电压 110kV/sqrt(3)频率 50Hz内阻抗设极小值Three-Phase Series RLC BranchElementsR0.1750L0.450/(2pi50)Y 连接Three-Phase Transformer (Two Windings)Elements20MVA110/11kVUk%10.5Pk135kWP022kWI0%0.8Three-Phase FaultElements0.02s 投入0.3s 退出三相短路接地电阻 0.001ΩThree-Phase V-I MeasurementMeasurements接在故障点前端测三相电压电流Powergui主界面连续仿真 Continuous50Hz连线顺序上故障模块的 Com 端必须接地否则短路回路不闭合。测量模块放在变压器低压侧与故障点之间这样记录的电流就是短路电流本体而不是电源侧电流。4.2 变压器参数填写的关键换算Simulink 变压器对话框中的 R、L 若选 pu 单位是以变压器自身额定值为基准的因此直接填 $R_T^* 0.00675$$X_T^* 0.1048$扣除电阻后的电抗分量不需要按 100MVA 重新折算。但如果选择填写有名值则需要把低压侧参数按变比折算到高压侧相对繁琐。用脚本生成参数最稳妥% 由铭牌参数换算 Simulink 变压器对话框所需标幺值 SN 20e6; Pk 135e3; P0 22e3; I0_pct 0.8; Uk_pct 10.5; RT_pu Pk / SN; % 0.00675 ZT_pu Uk_pct / 100; % 0.105 XT_pu sqrt(ZT_pu^2 - RT_pu^2); % 0.10478 % 填到 Simulink 的 R pu 和 L pu 中 fprintf(R pu %.5f\n, RT_pu); fprintf(L pu %.5f (电抗标幺值, 对话框中填 L 时按 LX/(2*pi*f))\n, XT_pu);参数说明Simulink 变压器对话框中的 L 单位是亨利当选择 pu 时填入的是标幺值但内部按 $L X_{pu}/(2\pi f)$ 换算成实际电感。励磁支路的 Rm 填入 909、Lm 填入 126 的标幺值。注意变压器连接方式要选 Y-Y与任务书要求的高低压绕组均为 Y 形连接一致。4.3 故障模块与求解器设置Three-Phase Fault 模块是建模中直接影响波形形态的元件。Transition 参数组中 Transitional 1 填[0.02 0.3]表示 0.02s 故障合闸、0.3s 故障切除Fault resistance 填0.001欧姆模拟金属性短路。Switch 组中把三相 Fault A、Fault B、Fault C 全部勾选即为三相对称短路。若只勾选单相或两相故障类型就变成不对称短路波形分析逻辑完全不同。仿真器设置建议Stop time 填 0.4s求解器选变步长 ode23tbMax step size 设为1e-4。短路暂态过程中非周期分量衰减时间常数只有几十毫秒0.4s 足够观察从短路发生到稳态的全过程最大步长限制到 0.1ms 是为了保证短路瞬间电流跳变的波形采集足够密否则冲击电流峰值可能因为采样点稀疏而被低估。提示Three-Phase Source 内部阻抗不要填 0SimPowerSystems 在部分版本中对纯理想电压源会报数值求解错误。一般填 0.001Ω 电阻加 0.001H 电感即可数值上远小于线路阻抗不影响短路计算结果。4.4 调试中常见问题与排查仿真最常见的报错是 Algebraic loop 或数值不收敛。遇到这种情况优先检查 Powergui 是否在模型中、故障模块是否接地、变压器连接组别是否选择正确。第二个常见问题是波形启动阶段出现高频振荡这通常是变压器励磁涌流与短路暂态叠加造成的属于正常现象分析时取 0.02s 之后的数据即可。如果短路后电流快速衰减到零说明故障模块的 Fault resistance 填得过大或者故障端口与测量模块接反。5. 结果对比分析无穷大系统与同步发电机短路的差异5.1 波形读取与冲击电流验证仿真运行后用 Scope 观察变压器低压侧三相电流波形。0.02s 前电流接近零空载0.02s 短路瞬间电流突增到数倍额定值0.03s 附近出现第一个也是最大的峰值这就是冲击电流。把理论值 19.8kA 与仿真波形最大值对比误差一般在 5% 以内说明模型参数设置正确。若需要量化对比用 To Workspace 模块把电流数据导出或用 MATLAB 脚本直接提取% 从仿真输出结构体提取电流波形并计算冲击电流 t out.simout.Time; % 时间序列 ia out.simout.Data(:, 1); % A相电流 [ish_sim, idx] max(ia); % 仿真冲击电流 t_sh t(idx); % 冲击电流出现时刻 fprintf(仿真冲击电流: %.2f kA, 出现在 %.3f s\n, ish_sim/1e3, t_sh);这段脚本用max检索短路后波形最大值输出冲击电流幅值和对应时刻。注意短路前空载电流接近零max的结果一定落在故障时段内不会误判。若波形中含高频噪声可先用medfilt1(ia, 5)做中值滤波再取峰值。5.2 无穷大系统与同步发电机短路电流特性对比对比项无穷大电源系统同步发电机机端短路周期分量幅值恒定不变次暂态、暂态、稳态三段衰减非周期分量按 TaL/R 衰减衰减更快受定子回路时间常数控制冲击电流倍数1.61.9可达 2 倍左右与 Xd 相关短路电流计算式IE/XΣIE/Xd稳态短路电流等于短路瞬间周期分量远小于初始值受励磁调节影响从对比可以看出无穷大系统短路电流分析的关键是正确计算短路回路总阻抗而同步发电机短路分析的核心是电机内部次暂态电抗 Xd 和暂态电抗 Xd 的分阶段衰减。若课程设计第二阶段要求扩展对比可以在同一模型中用 Synchronous Machine 模块替换无穷大电源对发电机机端设置三相短路对比两种工况下电流衰减形态的差异。5.3 仿真参数复用的通用化技巧将变压器参数、线路长度、短路时间等全部提取为脚本变量用set_param批量修改模型参数后反复仿真可以快速完成多组工况对比。例如把线路长度从 50km 改为 80km只需更新xL_pu并重新计算理论值Simulink 模型中线路模块的电阻和电感按长度线性更新。这种参数化思路也适用于后续做短路电流与保护整定配合分析将计算出的短路电流和冲击电流直接作为继电保护整定计算中电流速断保护和过流保护的输入值整个课程设计的工程完整性会提升一个档次。本文还有配套的精品资源点击获取 
长河君 · 2026/9/17 16:55:34 · 阅读 1.2k 无穷大电源三相短路故障仿真：从暂态理论到Simulink冲击电流验证
简介面向电气工程专业学生、科研人员及电力系统设计维护人员的MATLAB供配电系统短路故障仿真资料重点解决利用SimPowerSystem搭建无穷大功率电源供电系统模型、模拟三相短路并分析暂态过程的问题。压缩包内为1个PDF文件大小116KB内容包含课程设计任务书、技术参数、仿真要求及工作计划便于按步骤学习和打印查阅。已有92人浏览学习。内容围绕110kV无穷大系统在0.02s发生三相对称短路的场景完整给出变压器容量、短路电压、线路电阻电抗等关键元件的参数取值与计算方法并对比无穷大系统与同步发电机突然三相短路时电流变化的异同帮助读者深入理解短路暂态特性。同时覆盖从理论复习、SimPowerSystem建模调试、仿真结果分析到社会、安全、环境等因素思考的完整流程既适合课程设计参考也可作为工程项目实施与教学训练素材。1. 短路故障仿真第一步把无穷大电源的暂态过程读懂做供配电系统课程设计时很多人一上来就打开 Simulink 拖元件结果波形出来了却讲不清每个拐点对应什么物理过程。电力系统短路故障仿真里最容易被忽视、也最值得先搞懂的概念是三相短路瞬间电流最大值并不是故障刚发生的那个时刻而是大约半个周波之后。这个现象由非周期分量直流分量叠加造成是冲击电流计算的物理基础。本文以 MATLAB/SimPowerSystems现为 Simscape Electrical Specialized Power Systems为工具完整走一遍无穷大功率电源供电系统三相对称短路的理论推导、元件参数解析计算、仿真建模到结果分析的全流程适合电气工程专业学生做课程设计、毕业设计也适合刚接触电力系统仿真的工程技术人员对照排查自己的建模方法。2. 无穷大电源三相短路的暂态过程周期分量与非周期分量的分离2.1 无穷大电源的数学模型无穷大功率电源是电力系统短路计算中的理想化概念指电源容量远大于被研究系统、内阻抗为零、短路前后母线电压幅值和频率均保持恒定的电源。实际系统中 110kV 及以上电网的相对容量足够大时可近似按无穷大电源处理。设电源相电压为 $u_a U_m \sin(\omega t \alpha)$其中 $\alpha$ 为短路瞬间 A 相电压的初相角。短路前系统带纯阻感负载A 相电流表达式为$$i_a I_m \sin(\omega t \alpha - \varphi)$$其中 $I_m U_m / \sqrt{R^2 X^2}$$\varphi \arctan(X/R)$。这里 $R$、$X$ 是从电源到短路点之间的等值电阻与电抗。关键点在于短路瞬间电感中的磁链不能突变但短路后的电路阻抗突然变小稳态电流幅值会跳变到另一个值。为了满足磁链守恒约束回路中必然出现一个自由分量电流这个分量按指数规律衰减其初始值恰好抵消周期分量在短路瞬间的突变部分。2.2 短路电流的完整表达式与冲击电流对三相短路取 A 相分析短路后的电流可以写成强制分量与自由分量之和$$i_a I_{pm} \sin(\omega t \alpha - \varphi_k) \left[ I_m \sin(\alpha - \varphi) - I_{pm} \sin(\alpha - \varphi_k) \right] e^{-t/T_a}$$式中 $I_{pm}$ 是短路电流周期分量的幅值$\varphi_k$ 是短路回路阻抗角$T_a L_k / R_k$ 是短路回路非周期分量衰减时间常数。方括号里那一项就是非周期分量也叫直流分量它的大小取决于短路时刻的电压初相角 $\alpha$。当 $\alpha - \varphi \approx \pi/2$ 时非周期分量初始值最大合成电流在最严重情况下会在短路后约半个周期50Hz 系统为 0.01s达到峰值这个峰值称为冲击电流$$i_{sh} I_{pm} I_{pm} e^{-0.01/T_a} K_{sh} \sqrt{2} I$$$K_{sh} 1 e^{-0.01/T_a}$ 为冲击系数工程上高压系统一般取 1.61.9。这一结论直接指导仿真中的故障时间设置0.02s 设置短路则波形约在 0.028~0.03s 处出现第一个最高峰分析波形时应对照这个时刻验证模型是否正确。2.3 为什么 0.02s 是一个合理的短路时刻设定课程设计任务书中要求 0.02s 发生三相对称短路这个时间不是随便给的。电力系统继电保护的主保护动作时间通常在 2040ms 量级0.02s 短路对应的是保护尚未切除故障的最严重工况。仿真时还要设置一个故障切除时间比如 0.3s让故障模块在分析完暂态全过程后再断开目的是观察短路电流从暂态恢复到稳态的完整曲线。另外要注意任务书把短路点设在变压器低压侧出口。这意味着短路回路中包含线路阻抗和变压器短路阻抗两部分电源电压等级为 110kV短路点电压等级为 11kV做理论计算时必须用标幺值法或折算到同一电压等级后进行否则阻抗直接相加会出现数量级错误。3. 元件参数解析计算从变压器铭牌到标幺值阻抗3.1 基准值选取与线路参数计算理论计算采用标幺值法基准容量取 $S_B 100\text{MVA}$电压基准取各电压等级的平均额定电压110kV 侧取 $U_{B1} 115\text{kV}$11kV 侧取 $U_{B2} 10.5\text{kV}$。这样选取的动机是平均额定电压更接近实际运行电压且各电压等级间通过标准变比自然解耦。线路参数按任务书给定$L 50\text{km}$$x_1 0.4\Omega/\text{km}$$r_1 0.17\Omega/\text{km}$则线路总电抗 $X_L 20\Omega$、总电阻 $R_L 8.5\Omega$。高压侧基准阻抗为$$Z_{B1} \frac{U_{B1}^2}{S_B} \frac{115^2}{100} 132.25\Omega$$线路标幺值$$X_{L}^{} \frac{20}{132.25} 0.1512, \quad R_{L}^{} \frac{8.5}{132.25} 0.0643$$这部分计算的意义在于后续 Simulink 建模中如果直接填有名值需要把线路阻抗换算到 11kV 侧而用标幺值计算短路电流时则完全不受电压等级换算困扰。两套体系各有利弊建议理论分析用标幺值仿真建模直接在线路模块里填 50km、$r_1 0.17\Omega/\text{km}$、$x_1 0.4\Omega/\text{km}$ 这三个原始参数让 SimPowerSystems 自动处理。3.2 变压器短路阻抗与损耗参数折算变压器铭牌参数为 $S_N 20\text{MVA}$$U_k% 10.5$$\Delta P_k 135\text{kW}$$\Delta P_0 22\text{kW}$$I_0% 0.8$变比 $110/11\text{kV}$。先计算以变压器额定容量为基准的短路阻抗标幺值短路电抗标幺值 $X_T^* 10.5% 0.105$短路电阻标幺值 $R_T^* \Delta P_k / S_N 135\text{kW} / 20\text{MVA} 0.00675$短路阻抗模值 $Z_T^* \sqrt{0.105^2 0.00675^2} \approx 0.1052$可见电阻占比极小当统一换算到 $S_B 100\text{MVA}$ 基准时变压器电抗标幺值变为$$X_T^{*} 0.105 \times \frac{100}{20} 0.525$$这个数值远大于线路标幺电抗 0.1512说明在 110kV 供电的短路计算中变压器短路阻抗是决定短路电流大小的主因素。线路电阻和变压器电阻对短路电流幅值影响很小但计算冲击系数 $K_{sh}$ 时不能忽略因为 $T_a X / (\omega R)$ 对电阻敏感。3.3 短路电流理论值手算主要计算过程如下% 供配电系统三相短路电流理论计算 % 基准值 SB 100e6; % 基准容量 100MVA UB1 115e3; % 110kV侧平均额定电压 UB2 10.5e3; % 10.5kV侧平均额定电压 IB2 SB / (sqrt(3) * UB2); % 低压侧基准电流约 5499A % 元件标幺值SB基准 xL_pu (0.4 * 50) / (UB1^2 / SB); % 线路电抗标幺值 0.1512 rL_pu (0.17 * 50) / (UB1^2 / SB); % 线路电阻标幺值 0.0643 xT_pu 0.105 * (SB / 20e6); % 变压器电抗标幺值 0.525 rT_pu (135e3 / 20e6) * (SB / 20e6);% 变压器电阻标幺值 0.003375 X_sum xL_pu xT_pu; % 总电抗标幺值 0.6762 R_sum rL_pu rT_pu; % 总电阻标幺值 0.0677 I_pu 1 / sqrt(X_sum^2 R_sum^2); % 短路电流标幺值 1.4758 Ik I_pu * IB2; % 短路电流有名值约 8114A Ta X_sum / (2 * pi * 50 * R_sum); % 非周期分量时间常数约 0.0318s Ksh 1 exp(-0.01 / Ta); % 冲击系数约 1.73 ish Ksh * sqrt(2) * Ik; % 冲击电流约 19.8kA代码逻辑说明电压基准取 115kV 和 10.5kV 两级短路点基准电流用 10.5kV 侧计算线路阻抗标幺值需要除以高压侧基准阻抗变压器阻抗标幺值需要按容量折算不能直接用铭牌值 0.105最后短路电流标幺值等于 1 除以短路回路总阻抗模值。其中Ksh用 50Hz 的半个周期时间 0.01s 计算反映的是最严重短路时刻的冲击系数。3.4 励磁支路参数对短路计算的影响变压器空载损耗 $\Delta P_0 22\text{kW}$ 和空载电流百分数 $I_0% 0.8$ 用于计算励磁支路参数。励磁电阻 $R_m^* 1 / (\Delta P_0 / S_N) 909$励磁电抗 $X_m^* 1 / \sqrt{(I_0/100)^2 - (\Delta P_0/S_N)^2} \approx 126$。由于这两个标幺值远大于短路回路阻抗短路时流过励磁支路的电流极小因此手算短路电流时可以忽略励磁支路。但 Simulink 建模时需要填入这些参数否则变压器空载损耗特性和仿真启动瞬态与实际不符这也是理论计算与仿真调参过程中容易被忽视的差异点。4. SimPowerSystems 建模实操三相短路故障的参数设置与调试4.1 模型拓扑与元件清单SimPowerSystems 库在较新版本 MATLAB 中的路径为 Simscape Electrical Specialized Power Systems核心求解模块 Powergui 必须拖入模型否则仿真报错。模型拓扑按电源到短路点的功率流向布置三相电源 - 三相串联 RLC 支路模拟 50km 线路- 双绕组三相变压器 - 三相故障模块 - 接地。元件库路径本例参数设置Three-Phase SourceElectrical Sources相电压 110kV/sqrt(3)频率 50Hz内阻抗设极小值Three-Phase Series RLC BranchElementsR0.1750L0.450/(2pi50)Y 连接Three-Phase Transformer (Two Windings)Elements20MVA110/11kVUk%10.5Pk135kWP022kWI0%0.8Three-Phase FaultElements0.02s 投入0.3s 退出三相短路接地电阻 0.001ΩThree-Phase V-I MeasurementMeasurements接在故障点前端测三相电压电流Powergui主界面连续仿真 Continuous50Hz连线顺序上故障模块的 Com 端必须接地否则短路回路不闭合。测量模块放在变压器低压侧与故障点之间这样记录的电流就是短路电流本体而不是电源侧电流。4.2 变压器参数填写的关键换算Simulink 变压器对话框中的 R、L 若选 pu 单位是以变压器自身额定值为基准的因此直接填 $R_T^* 0.00675$$X_T^* 0.1048$扣除电阻后的电抗分量不需要按 100MVA 重新折算。但如果选择填写有名值则需要把低压侧参数按变比折算到高压侧相对繁琐。用脚本生成参数最稳妥% 由铭牌参数换算 Simulink 变压器对话框所需标幺值 SN 20e6; Pk 135e3; P0 22e3; I0_pct 0.8; Uk_pct 10.5; RT_pu Pk / SN; % 0.00675 ZT_pu Uk_pct / 100; % 0.105 XT_pu sqrt(ZT_pu^2 - RT_pu^2); % 0.10478 % 填到 Simulink 的 R pu 和 L pu 中 fprintf(R pu %.5f\n, RT_pu); fprintf(L pu %.5f (电抗标幺值, 对话框中填 L 时按 LX/(2*pi*f))\n, XT_pu);参数说明Simulink 变压器对话框中的 L 单位是亨利当选择 pu 时填入的是标幺值但内部按 $L X_{pu}/(2\pi f)$ 换算成实际电感。励磁支路的 Rm 填入 909、Lm 填入 126 的标幺值。注意变压器连接方式要选 Y-Y与任务书要求的高低压绕组均为 Y 形连接一致。4.3 故障模块与求解器设置Three-Phase Fault 模块是建模中直接影响波形形态的元件。Transition 参数组中 Transitional 1 填[0.02 0.3]表示 0.02s 故障合闸、0.3s 故障切除Fault resistance 填0.001欧姆模拟金属性短路。Switch 组中把三相 Fault A、Fault B、Fault C 全部勾选即为三相对称短路。若只勾选单相或两相故障类型就变成不对称短路波形分析逻辑完全不同。仿真器设置建议Stop time 填 0.4s求解器选变步长 ode23tbMax step size 设为1e-4。短路暂态过程中非周期分量衰减时间常数只有几十毫秒0.4s 足够观察从短路发生到稳态的全过程最大步长限制到 0.1ms 是为了保证短路瞬间电流跳变的波形采集足够密否则冲击电流峰值可能因为采样点稀疏而被低估。提示Three-Phase Source 内部阻抗不要填 0SimPowerSystems 在部分版本中对纯理想电压源会报数值求解错误。一般填 0.001Ω 电阻加 0.001H 电感即可数值上远小于线路阻抗不影响短路计算结果。4.4 调试中常见问题与排查仿真最常见的报错是 Algebraic loop 或数值不收敛。遇到这种情况优先检查 Powergui 是否在模型中、故障模块是否接地、变压器连接组别是否选择正确。第二个常见问题是波形启动阶段出现高频振荡这通常是变压器励磁涌流与短路暂态叠加造成的属于正常现象分析时取 0.02s 之后的数据即可。如果短路后电流快速衰减到零说明故障模块的 Fault resistance 填得过大或者故障端口与测量模块接反。5. 结果对比分析无穷大系统与同步发电机短路的差异5.1 波形读取与冲击电流验证仿真运行后用 Scope 观察变压器低压侧三相电流波形。0.02s 前电流接近零空载0.02s 短路瞬间电流突增到数倍额定值0.03s 附近出现第一个也是最大的峰值这就是冲击电流。把理论值 19.8kA 与仿真波形最大值对比误差一般在 5% 以内说明模型参数设置正确。若需要量化对比用 To Workspace 模块把电流数据导出或用 MATLAB 脚本直接提取% 从仿真输出结构体提取电流波形并计算冲击电流 t out.simout.Time; % 时间序列 ia out.simout.Data(:, 1); % A相电流 [ish_sim, idx] max(ia); % 仿真冲击电流 t_sh t(idx); % 冲击电流出现时刻 fprintf(仿真冲击电流: %.2f kA, 出现在 %.3f s\n, ish_sim/1e3, t_sh);这段脚本用max检索短路后波形最大值输出冲击电流幅值和对应时刻。注意短路前空载电流接近零max的结果一定落在故障时段内不会误判。若波形中含高频噪声可先用medfilt1(ia, 5)做中值滤波再取峰值。5.2 无穷大系统与同步发电机短路电流特性对比对比项无穷大电源系统同步发电机机端短路周期分量幅值恒定不变次暂态、暂态、稳态三段衰减非周期分量按 TaL/R 衰减衰减更快受定子回路时间常数控制冲击电流倍数1.61.9可达 2 倍左右与 Xd 相关短路电流计算式IE/XΣIE/Xd稳态短路电流等于短路瞬间周期分量远小于初始值受励磁调节影响从对比可以看出无穷大系统短路电流分析的关键是正确计算短路回路总阻抗而同步发电机短路分析的核心是电机内部次暂态电抗 Xd 和暂态电抗 Xd 的分阶段衰减。若课程设计第二阶段要求扩展对比可以在同一模型中用 Synchronous Machine 模块替换无穷大电源对发电机机端设置三相短路对比两种工况下电流衰减形态的差异。5.3 仿真参数复用的通用化技巧将变压器参数、线路长度、短路时间等全部提取为脚本变量用set_param批量修改模型参数后反复仿真可以快速完成多组工况对比。例如把线路长度从 50km 改为 80km只需更新xL_pu并重新计算理论值Simulink 模型中线路模块的电阻和电感按长度线性更新。这种参数化思路也适用于后续做短路电流与保护整定配合分析将计算出的短路电流和冲击电流直接作为继电保护整定计算中电流速断保护和过流保护的输入值整个课程设计的工程完整性会提升一个档次。本文还有配套的精品资源点击获取