简介面向地铁牵引供电系统设计、运维技术人员及高校电力系统仿真方向师生的技术文档围绕Matlab/Simulink平台展开地铁直流牵引供电系统的完整建模与故障分析。内容涵盖外部电源、24脉波整流变电站、直流母线、接触网PI分布参数模型及钢轨大地回路并延伸到牵引传动系统的SVPWM逆变控制与故障工况仿真对成都地铁实际参数进行复现验证仿真结果与现场故障录波对比可用于理解建模关键技术、掌握故障特性分析方法。资源包为1个docx文档约49KB内含可运行的复现代码及逐段中文注释涉及24脉波整流器、接触网特征阻抗计算、空间矢量PWM扇区判断与系统级主仿真脚本等模块兼顾理论推导与工程落地。已有183人学习适合需要快速搭建仿真框架、对照代码理解电力电子系统建模细节的读者参考。1. 直流牵引供电仿真的误差大头往往不在整流器很多人拿到“地铁牵引供电系统建模与仿真”这个题目第一反应是把 24 脉波整流器的二极管桥搭出来看着直流母线纹波小就以为模型成了。真做过现场对录波的人会知道空载电压对得上、满载压降对不上问题往往出在钢轨—大地回流路径和接触网分布参数上。直流牵引是走行轨回流体系电流从接触网去、经钢轨回钢轨对地并非理想绝缘泄漏电流在地层里扩散这条路径的阻抗随频率和土壤分层变化直接决定短路电流峰值和馈线保护定值。本篇按外部电源、整流变电站、直流母线、接触网、列车负载、回流网络的顺序把各子系统在 Simulink 里怎么落地、参数从哪来、故障工况怎么注入讲清楚。适合做供电系统设计、继保整定以及拿这条线做毕业设计的同学。2. 24 脉波整流变电站的移相变压器与直流母线建模2.1 为什么是 24 脉波而不是 12 脉波城市轨道交通牵引变电所普遍采用 24 脉波整流机组两台机组并联运行。原理是用两段相位差 15° 的阀侧绕组各自接三相全桥直流侧并联输出。12 脉波的特征谐波是 11、13 次24 脉波把特征谐波推到 23、25 次网侧电流 THD 从 8% 左右压到 3% 以内这是城市电网对地铁这种大功率非线性负荷的基本准入条件。理解这一点才能明白为什么模型里“移相”这一步绝对不能省。建模时最容易踩的坑是以为在 Simulink 里放两个 Universal Bridge 再并联就算 24 脉波了。没有移相变压器两个桥的输入电压同相并联后谐波是叠加而不是抵消直流纹波反而更差。判断移相是否真正生效不看模型结构看直流侧纹波的主频。2.2 整流机组的关键参数参数典型值说明网侧电压35 kV部分线路 110 kV 直降阀侧电压1180 V空载直流约 1590 V移相角±7.5° 或 0°/15°两组之间差 15°联结组Dyn11 / Dy0配合实现 15° 相位差理想空载直流电压√2·1180·(3√3/π) ≈ 1593 V理论值额定直流电压1500 V牵引网标称值参数表看起来简单但每一行都对应一个仿真里的可调量。阀侧电压偏高空载电压就对不上移相角给错纹波主频就错。整定这些参数的原则是先让空载工况的直流电压和纹波对上现场再谈负载和故障。2.3 移相变压器组的搭建代码function tf_group build_phase_shift_transformers() % 网侧 35kV两个阀侧绕组相位相差 15 度 base struct(PrimaryVoltage, 35e3, ... SecondaryVoltage, 1180, ... VectorGroup, Dyn11, ... PhaseShift, 0); tf1 ThreePhaseTransformer(base); base.PhaseShift 15; % 第二组移相 15 度 base.VectorGroup Dy0; tf2 ThreePhaseTransformer(base); tf_group [tf1, tf2]; % 说明移相通过阀侧绕组接法星/三角配合实现 % 两组阀侧分别接三相全桥直流侧串联或并联 end逻辑说明移相角在 Simulink 的 Three-Phase Transformer (Two Windings) 里由 Winding 2 的 connection 参数决定Dyn11 与 Dy0 组合天然产生 30° 相位差要细分到 15° 需要延边三角形或 Zigzag 绕组。参数说明PhaseShift 只是工程标记真正决定相位的是绕组接法和匝数比不要把它当成仿真里会直接生效的物理量否则跑出来的纹波会骗人。2.4 直流母线 RLC 与纹波校验整流器输出到直流母线之间有平波电抗器和母线电容。验证模型是否正确的第一关是空载纹波第二关才是满载压降。% 直流母线纹波快速校验 Vdc0 sqrt(2)*1180*(3*sqrt(3)/pi); % 理想空载直流电压 ≈1593V L_smooth 2e-3; % 平波电抗 2 mH C_bus 10e-3; % 母线电容 10 mF f_line 50; % 网侧频率 f_ripple 24*f_line; % 24 脉波纹波基频 1200 Hz % 简单 LC 二阶梯波衰减估算 atten 1 / (2*pi*f_ripple)^2 / (L_smooth*C_bus 1e-12); fprintf(空载纹波衰减系数 ≈ %.2e\n, atten);逻辑说明纹波基频取 24×501200 Hz不是 300 Hz这是判断模型是否真的 24 脉波的关键。参数说明L_smooth 和 C_bus 的乘积决定纹波衰减母线电容取值要在纹波和故障电流上升率之间折中取值太大短路电流 di/dt 会被拉低和现场录波对不上。如果仿真出来的纹波主频落在 300 Hz回头检查变压器接法几乎可以确定移相没生效。3. 接触网 PI 链式电路与钢轨—大地回流阻抗建模3.1 集中参数为什么不成立接触网随线路长度变化一段馈线供电半径常见 13 km工频和纹波频段下必须用分布参数。PI 型链式电路是最省算力的做法把馈线切成 N 段每段用串联 RL 和两端各半个并联 CG。这样既保留了行波特性又把状态变量数量控制在可控范围。参数成都地铁取值单位R per km0.02Ω/kmL per km1.2e-3H/kmC per km12e-9F/kmG per km1e-6S/km这些参数不是拍脑袋来的要按实际接触网型号、悬挂方式和回流轨截面换算。同一批参数下分段数不同算出来的短路峰值能差 5% 以上所以参数和分段数必须一起标定。3.2 特征阻抗与分段数的选择function [Zc, seg_num] contact_line_model(freq, length_km) R 0.02; L 1.2e-3; C 12e-9; G 1e-6; Zs R 1j*2*pi*freq*L; % 串联阻抗 Yp G 1j*2*pi*freq*C; % 并联导纳 Zc sqrt(Zs/Yp); % 特征阻抗 % 分段数每段长度不超过最高关注频率对应波长的 1/10 v 1/sqrt(L*C); % 行波速度 lambda v / freq; seg_num max(4, ceil(length_km / (lambda/10))); fprintf(特征阻抗 %.3f∠%.2f° Ω建议分段 %d\n, ... abs(Zc), angle(Zc)*180/pi, seg_num); end逻辑说明特征阻抗决定馈线末端短路时反射系数的量级分段数决定仿真能否捕捉行波过程。参数说明freq 取关注的最高频率做保护整定关注 50 Hz 足够做纹波或暂态要按 1200 Hz 甚至更高算。分段数不是越多越好每段都引入状态变量段数翻倍仿真步长要跟着减半算力成本平方级上升。3.3 钢轨对地不是绝缘四层大地模型走行轨通过扣件、道床对地存在过渡电阻回流电流一部分经钢轨回一部分经大地回。Carson 方程给出地回路阻抗但完整形式要算 Bessel 函数的无穷积分Simulink 里跑不动工程上通常做分层加权简化。function Z_earth carson_equation(freq, depth, rho) % depth: 各层深度(m)最后一层用 Inf % rho : 各层电阻率(Ω·m) mu0 4*pi*1e-7; omega 2*pi*freq; % 分层土壤等效电阻率简化按深度加权 rho_eq 0; for k 1:numel(rho) if k numel(rho) w 0.1; % 底层权重 else w min(depth(k), 5) / 5; % 浅层权重更高 end rho_eq rho_eq w*rho(k); end Z_earth (mu0*omega/8) 1j*omega*(mu0/(2*pi))*log(2.0/1.0) ... rho_eq*freq*1e-6; % 经验修正项 end逻辑说明这里给的是工程可用的简化形式用一个加权等效电阻率替代分层积分的复杂计算。参数说明depth 按地层分层给成都地区常见表层黏土、下层砂卵石电阻率从 100 降到 30 Ω·m最后两层权重自动降低。钢轨自身取 R0.02 Ω/km、L1.5e-3 H/km两部分相加得到总回流阻抗。这个模型在 50 Hz 到 1 kHz 的误差可以接受再高的频段需要回到精细 Carson 展开式。提示如果只关心直流侧稳态压降钢轨—大地阻抗可以先用常数 0.05 Ω/km 顶替等整流和逆变调通之后再换成频变模型调试效率会高很多。4. SVPWM 牵引逆变与列车负载的联合仿真4.1 矢量控制与 SVPWM 的分工牵引传动做的是把直流 1500 V 逆变成变频变压的三相交流驱动异步或永磁电机。外层是转子磁场定向的电流环内层是 SVPWM 输出六路驱动。两层在 Simulink 里不要混在一个 S-Function 里调试时很难定位是控制参数不对还是扇区判断出错。分开成两个 Function-Call Subsystem中间用电压矢量传递出问题一眼能看出是外层还是内层。4.2 电流环参数表参数取值说明d 轴电感 Ld8 mH表贴式 PMSM 时 Ld≈Lqq 轴电感 Lq12 mH凸极率约 1.5定子电阻 Rs0.2 Ω影响低速转矩精度电流环 Kp_d0.5按带宽取开关频率的 1/10电流环 Ki_d10消除稳态误差开关频率2 kHzIGBT 常用这些值对应的是 1.5 MW 级牵引电机换功率等级时电阻电感要按比例缩放。Kp 的整定原则是先按带宽 1/10 开关频率推算再在仿真里阶跃响应观察超调最终以调整时间小于 10 ms 为准。4.3 SVPWM 扇区判断与作用时间function duty svpwm_control(Vref, Vdc) % Vref: [Valpha, Vbeta] Valpha Vref(1); Vbeta Vref(2); theta atan2(Vbeta, Valpha); theta mod(theta, 2*pi); % 归一到 [0, 2pi) sector floor(theta/(pi/3)) 1; % 扇区 1..6 % 相邻基本矢量的作用时间 T1 sqrt(3)*abs(Vref)*sin(sector*pi/3 - theta)/Vdc; T2 sqrt(3)*abs(Vref)*sin(theta - (sector-1)*pi/3)/Vdc; T0 max(0, 1 - T1 - T2); % 零矢量时间 duty struct(sector, sector, ... T1, T1, T2, T2, T0, T0, ... saturate, (T1T2) 1); % 过调制标志 end逻辑说明T1、T2 是两个相邻基本电压矢量的作用时间T0 是零矢量时间三段加起来等于一个开关周期。参数说明saturate 标志很重要当 T1T21 说明参考矢量超出内切圆需要做过调制处理否则输出波形顶部被削平电机电流出现明显低次谐波。Vdc 用实时母线电压而不是额定 1500 VPWM 才能在网压波动时自动补偿。4.4 列车负载曲线的接入% 单车运行阻力与加速度曲线 v_max 80/3.6; % 最高速度 m/s m_train 280000; % 列车质量 kg R_run (v) 5.2 0.045*v 0.00086*v.^2; % 运行阻力 N a_cmd [0 1.0; 10 0.8; 20 0; 30 -0.6; 40 0]; % 加速度指令(m/s^2) t_cmd a_cmd(:,1); a_val a_cmd(:,2); F_trac (t,v) m_train*interp1(t_cmd, a_val, t, linear, 0) R_run(v);逻辑说明牵引力由加速度指令加运行阻力构成负载功率再乘速度得到接到直流母线上作为受控电流源。参数说明运行阻力公式的二次项来自空气阻力速度单位必须是 m/s直接填 km/h 算出的功率偏小一个数量级这是新手最常犯的错误。多车同时启动时用 this 累加母线电压跌落幅度是校核馈线保护的核心指标。5. 短路故障注入与求解器级校核技巧5.1 短路故障的参数化注入function set_fault(model, t_fault, duration, R_fault) % 正极对地短路注入 set_param([model /Fault_Breaker], ... SwitchTimes, mat2str([t_fault, t_faultduration]), ... InitialState, 0); % 0断开 set_param([model /Fault_Resistance], ... Resistance, num2str(R_fault)); % 过渡电阻 Ω end逻辑说明用 Breaker 的 SwitchTimes 做时序开关比在模块里用 ClockRelational Operator 触发更稳变步长求解器下不会因过零检测丢事件。参数说明过渡电阻 R_fault 取 0.01 Ω 对应金属性短路取 0.52 Ω 对应经电弧或道床电阻的非金属性短路两者峰值电流差两三倍整定保护时都要跑一遍。5.2 校核指标与实测对比场景实测仿真允许误差空载母线电压1650 V1635 V 2%满载母线压降1420 V1405 V 2%正极接地短路峰值12.3 kA12.1 kA 5%短路电流上升率 di/dt按录波仿真值 10%四个指标里最容易超标的是 di/dt它受母线电容和接触网分段数共同影响调纹波时电容加大di/dt 就偏小。两个指标互相拉扯需要在参数上做折中一般先满足 short-circuit 峰值再回头微调分段数把 di/dt 拉回来。5.3 求解器与步长的取舍Simulink 里电力电子系统常用 ode23tb 或 ode15s。变步长下若故障事件时刻被跨过短路电流峰值会偏低。稳妥做法是把 Max step size 固定到开关周期的 1/202 kHz 开关对应 25 μs故障时刻再在求解器里加一个 Output times 断点。这样峰值误差能压到 3% 以内。注意功率器件用 ideal switch 时步长太大会出现数值振荡母线电压波形上看到寄生的高频毛刺那不是模型物理特性是求解器伪影不要当成谐振去解释。5.4 多区间仿真的一条提速技巧整条线路多个区间一起仿真时接触网链式电路的状态变量数量是瓶颈。把不关注的远端区间用戴维南等效压缩成一个电压源加阻抗只在故障点附近保留完整 PI 分段单次仿真时间通常能降 60% 以上对故障点近端波形的影响可忽略。这是把“能跑通”推进到“能迭代参数”的关键一步很多时候模型精度不是不够而是跑一次太慢导致没法调参。本文还有配套的精品资源点击获取