做牵引供电系统仿真这个方向我踩了不少坑也积累了一些可以复用的经验。这次把整套25kV交流电铁路牵引供电系统在Matlab/Simulink里完整搭了一遍从牵引变电所、接触网到EH800机车的变压-整流-逆变链路再到1kV直流母线载荷特性分析把之前零散的知识点串成了闭环。这套模型我用在了JR EH800机车载荷研究上仿真出来的趋势和实际线路运行数据能对上所以把整个建模过程、参数计算和踩过的坑整理出来给做铁路电气化、电力电子变流和控制方向的朋友一个可以直接参考的落地模板。1. 系统结构拆解与建模思路1.1 25kV交流牵引供电系统从变电站到受电弓先捋一下电气化铁路的供电路径。发电厂发出的电通过公网输电线路传输到牵引变电所牵引变电所里的牵引变压器把电压从110kV或220kV降到27.5kV工程上标称25kV再通过接触网就是铁路正上方的悬链线结构输送到机车顶部的受电弓。电流经过机车负载后通过钢轨和回流线回到牵引变电所形成完整的回流回路。这个结构和普通电力用户有一个本质差异牵引负荷是单相大功率动态负荷而且机车是不断运动的所以会带来接触网电压波动、负序电流、谐波等一系列问题。在做仿真建模的时候我刻意没有把供电侧做成理想电压源而是保留了牵引变压器的漏感和短路阻抗这样才能看到机车加载时网侧电压的真实跌落情况。这一点特别重要如果你的模型里供电侧是理想源那你后面做的所有载荷分析都会偏乐观接触网压降对机车功率输出上限的影响根本体现不出来。1.2 EH800机车与1kV中间直流环节EH800是JR北海道的货运电力机车采用交-直-交传动方式最大功率6720kW。受电弓把25kV单相交流引入机车后先经车载主变压器降压然后通过四象限脉冲整流器PWM整流器变成直流这个直流母线就是中间直流环节。标题里提到的1kV我理解为这个中间直流母线的电压等级在仿真模型里按DC 1000V来做整个中间回路。直流母线之后是三相两电平逆变器把直流再逆变成三相变频交流电驱动异步牵引电机。整个链路可以分成三段来研究供电侧25kV交流、中间侧1kV直流、传动侧变频交流。载荷研究的关键就是看这三段在不同机车工况下各自的电压、电流、功率怎么动态变化。这里补充一个背景EH800是日本货运机车中首款能适应50Hz/60Hz两种电网频率的交流电力机车所以模型里电源频率这个参数建议做成可配置的方便扩展研究。2. Matlab/Simulink建模方案选型2.1 为什么选Simscape Electrical而不是纯数学模型建模方式上我试过两条路。第一条是用Simulink里的基础数学模块把供电系统、整流器、逆变器、电机全部用微分方程和传递函数搭出来。这条路模型轻、跑得快但有个致命问题——谐波分析做不精确因为开关器件被简化成了理想受控源无法真实反映PWM调制的脉动过程。第二条就是用Simscape Electrical老名字叫SimPowerSystems里面有现成的三相变压器、IGBT、二极管、异步电机模型能直接模拟开关过程网侧电流谐波、直流纹波这些细节都能看到。我的选择是主电路用Simscape Electrical控制系统用标准Simulink模块搭。这样既有物理模型的真实性控制算法又能自由修改。使用Simscape Electrical时要注意模型里必须放一个Powergui模块它负责配置求解器和分析工具。这个模块不是摆设仿真步长设置、FFT分析、阻抗测量都要靠它。很多新手上来就搭模型忘了放Powergui一仿真就报错这个坑我后面会详细说。2.2 分块建模五层结构我的模型按物理层级分成五块每块单独封装成子系统Subsystem模块之间通过电气端口和信号线连接。这个分层方式虽然看起来笨但确实是排查问题最有效的手段调试的时候可以单独看每一级的波形不用从头查到尾。电网与牵引变压器模块三相电压源、单相牵引变压器、接触网等效阻抗车载主变压器模块EH800机车内的高压变压器含漏抗四象限整流器模块IGBT桥臂、PWM控制器、中间直流电容逆变器与牵引电机模块三相两电平逆变器、异步电机、负载转矩输入控制逻辑模块整流器双闭环控制、逆变器矢量控制每一层封装好以后用总线信号把控制指令在层间传递。比如整流器控制层需要采集直流母线电压和网侧电流这些信号通过Goto/From标签跨子系统传递减少连线杂乱。模型整体跑通之后想改某一层的参数直接进子系统改就行不影响其他部分。2.3 版本和工具箱说明我用的环境是Matlab R2023a需要安装Simscape Electrical工具箱在Add-On Explorer里搜Simscape Electrical。如果你用的是老版本模块库的名字会不一样R2018a之前叫SimPowerSystemsR2018b到R2021a叫Simscape Electrical Specialized Power Systems。命名变化不影响模块使用方法但从网上搜参考案例时要注意对应版本。有一个小技巧Simulink的模块库浏览器里可以按CtrlF搜索模块名称。比如要找单相变压器直接搜Single-Phase Transformer比在电气库的菜单树里一层层翻快得多。如果你要用的模块在当前版本里找不到多半是名字变了或者被合并到别的库了去MathWorks官网文档里查一下模块别名就行。3. 关键参数计算与模块实现3.1 牵引变电所参数与容量选择牵引变压器的容量不能拍脑袋决定。EH800最大功率是6720kW考虑机车效率约0.88-0.9从电网侧取的电功率大概在7500kW左右。再考虑功率因数按0.95算单台机车所需容量大约8000kVA。实际牵引变电所还要同时带多趟车所以我的模型里设一个25MVA左右的单相牵引变压器短路阻抗取10.5%这是牵引变电所比较常见的配置。Matlab里用Simscape Electrical的Single-Phase TransformerTwo Windings模块参数设置如下参数项数值说明额定容量25MVA满足多车同时取电需求频率50Hz可选50/60Hz切换一次侧电压110kV公网侧电压二次侧电压27.5kV接触网标称25kV对应的实际空载电压短路阻抗百分比10.5%影响电压跌落和谐波特性如果要做接触网电压降落分析还需要在变压器二次侧串联一段线路阻抗。接触网的单位阻抗大概是0.15Ω/km左右我按10km接触网长度设置了等效阻抗。铁路牵引供电系统对电压波动有要求通过这个阻抗能模拟出机车加载时27.5kV母线电压跌落的情况。想验证这个阻抗设置是否合理可以在空载和满载两种情况下测量接触网末端电压对比电压跌落比例是否在5%到10%的合理范围内。3.2 车载主变压器与整流器主电路参数车载主变压器的作用是把接触网来的25kV降到适合整流器工作的电压等级。单相PWM整流器是升压型变流器直流母线电压U_dc必须高于交流侧电压峰值一般取1.2到1.5倍裕量。对于U_dc1000V的直流母线理想情况下交流侧电压有效值要低于1000/√2≈707V。考虑到IGBT的耐压和调制比限制我取车载变压器二次侧电压有效值为650V对应峰值约920V调制比留了约8%的裕量。整流器主电路是两桥臂的IGBT单相全桥结构交流侧串联进线电感L。电感取值对电流纹波影响很大经验公式是L U_dc / (8 × f_sw × ΔI_max)其中f_sw是开关频率取2000HzΔI_max是允许电流纹波峰值按额定电流的15%-20%估算。EH800单单元额定功率约3360kW直流电压1000V所以额定直流电流约3360AΔI_max按500A算L 1000 / (8×2000×500) ≈ 0.125mH。实际取0.15mH留一点裕量。进线电感的另外一个作用是滤波它和整流器交流侧电压共同决定了电流的谐波含量。如果电感太小网侧电流里的高次谐波会明显增大如果电感太大电流内环的响应速度会变慢动态性能下降。0.15mH这个值在谐波抑制和动态响应之间比较平衡。3.3 直流母线电容计算与二倍频脉动问题单相整流器有一个固有特性瞬时功率以两倍电网频率脉动所以直流电压必然会有100Hz的纹波。如果不加处理这个纹波会传导到逆变器侧直接影响电机控制精度。实际工程中EH800这样的机车会在中间直流环节加二次谐振滤波支路LC谐振在100Hz我在模型里也采用了这个方案。先说支撑电容的估算。设脉动功率幅值为P对单相PWM整流器脉动功率幅值等于平均功率直流母线电压U_dc1000V允许纹波幅值ΔU20V即纹波峰峰值40V对应2%的电压波动。根据能量守恒关系电容储能的变化要能吸收脉动能量C P / (ω × ΔU × U_dc)其中ω2π×50314.16rad/sP3.36MW计算得到C ≈ 0.535F也就是535000μF。这个电容值非常大实际工程中不可能用这么大容量的电容去做支撑这正是需要谐振支路的原因。二次谐振支路的设计方法如下让LC串联支路谐振在100Hz即满足LC 1 / (2π × 100)^2 ≈ 2.53e-6取L0.5mH则C 2.53e-6 / 0.5e-3 ≈ 5.07mF。这个谐振支路专门吸收100Hz的纹波电流支撑电容只用取4-5mF就够了体积和成本大幅下降。我在模型里加了谐振支路后100Hz纹波从±80V压到了±15V左右效果非常明显。3.4 四象限整流器双闭环控制系统实现四象限整流器的控制是整个仿真中最费时间的部分。控制目标有两个把直流电压稳定在1000V让网侧电流正弦化且功率因数接近1。系统结构上用经典的电压外环加电流内环双闭环。电压外环检测直流母线电压U_dc与给定值1000V做差经过PI调节器后输出网侧电流指令的幅值I_ref。这个I_ref要和电网电压同步通过锁相环PLL得到电网电压的相位θ然后用cosθ做同步信号生成电流指令i_ref I_ref × cosθ。电流内环检测实际网侧电流i_s与指令i_ref做差经过PI调节器后输出调制信号。由于交流侧存在电感L系统方程是U_s L×di/dt U_ab其中U_ab为整流器交流侧电压为了让电流准确跟踪指令调制信号里还要加入电网电压前馈分量。PI参数整定我用的是工程设计法。电流内环的截止频率取得比开关频率低5到10倍开关频率2kHz电流环带宽取200Hz左右。这时Kp ≈ L×ω_c 0.00015×2π×200 ≈ 0.19Ki ≈ R×ω_c进线电阻取0.01ΩKi约12.6。实际调的时候Kp取0.15Ki取10波形就比较干净了。注意这个计算值是起点不是终点仿真里要微调。电压外环的带宽不能太高否则整流器从电网抽走的电流幅值指令会震荡。我取电压环带宽为电流环的五分之一左右即30-40Hz。P参数先给Kp2Ki30然后边看波形边调。如果直流电压出现持续振荡大概率是电压环Kp给大了。3.5 逆变器与异步牵引电机模型中间直流1000V经过三相两电平逆变器用SVPWM调制输出三相交流驱动异步牵引电机。EH800每个单元配4台三相异步牵引电机每台额定功率约840kW总功率6720kW除以8台电机。我在模型里用一台等效电机代表一个单元的4台电机也就是额定功率3360kW的等效电机模型这样能减少计算量同时保证载荷特性研究的准确性。异步电机的等效电路参数在Simulink的Asynchronous Machine SI Units模块里设置需要填定子电阻、转子电阻、定子漏感、转子漏感、励磁电感、转动惯量等。我没有EH800的实测电机参数用的是同等级大功率牵引电机参数的近似值这在仿真是常见做法。关键是要保证额定功率和额定电压匹配否则转矩-转速特性曲线会失真。逆变器控制我用的是转子磁场定向矢量控制FOC把定子电流分解成励磁分量和转矩分量分别用两个PI控制器控制。再上层是转速外环和转矩限幅。载荷研究中我给电机加一个负载转矩曲线来模拟不同的牵引工况。矢量控制的两个电流内环PI参数可以参考异步电机参数来设计励磁电流环带宽和转矩电流环带宽可以取一样都取200Hz左右。4. 载荷研究仿真流程与结果分析4.1 典型工况设计与仿真设置载荷研究的重点不在稳态波形而在动态过程中各环节怎么响应。我设计了四个典型工况工况A满载起动加速从0加速到60km/h工况B恒速巡航过程中负载转矩阶跃增加工况C坡道载荷模拟长大坡道持续大功率输出工况D再生制动能量从电机回馈到直流母线仿真参数方面求解器在Powergui里选离散Discrete模式步长100μs。仿真总时长30s。整流器开关频率2kHz逆变器开关频率1kHz。离散模式比连续模式跑得快很多因为Simscape Electrical的连续模式需要极小步长才能精确捕捉开关动作计算量大到难以接受。离散模式下把步长设为开关周期的二十分之一到五十分之一既能捕捉PWM细节运行速度也可控。4.2 供电侧响应网压跌落与谐波机车加速时接触网电流会从几十安快速冲到数百安。由于接触网有阻抗网压会有明显跌落。在我的模型里机车主变压器一次侧电流峰值能到250A左右对应功率约6720kW、25kV电压、功率因数0.95的情况。接触网电压跌落幅度在4%到5%之间这个数值在铁路牵引供电系统允许范围内。谐波方面加速过程中网侧电流的THD大约在5%到8%之间主要取决于整流器PWM的调制比和进线电感的大小。如果THD超标优先调整进线电感其次是提高开关频率。把进线电感从0.15mH加大到0.2mHTHD能从8%降到5%左右但代价是电流环响应变慢这个要权衡。4.3 直流母线动态响应这是载荷研究中最有意思的部分。当转矩从50%阶跃到100%时直流母线电压会先跌下去十几伏然后经过大约0.5秒回到1000V。这个过程的过渡时间和电压跌落幅值直接反映电压外环PI参数整定得合不合理。如果跌落超过50V说明电压环响应太慢可以适当加大Kp如果出现振荡说明Kp太大。稳态下的直流电压纹波主要来自两方面单相整流器的二倍频脉动100Hz和逆变器的开关谐波。我在模型里加了谐振支路后100Hz纹波从±80V压到了±15V。这里要强调一个细节谐振支路的电容初始电压要设成0否则仿真开始时会有一个很大的电流冲击。4.4 电机侧响应与载荷突变矢量控制下转矩电流和励磁电流可以实现解耦控制。加速阶段电流内环快速跟随电机转速平稳上升。载荷突变时转速跌落不超过2%大约3秒内恢复到设定值。这些动态指标能说明整套模型的控制参数基本合理。再生制动工况下电机的机械能通过逆变器回馈到直流母线直流电压会先上升然后整流器进入逆变状态把能量回馈到电网。这时候能看到网侧电流的相位反转180度功率倒流这是四象限整流器区别于不可控整流器的主要特征。在做能量回馈研究时直流母线电压的上升幅度是重要指标如果超过1100V说明整流器回馈功率受限需要增大电流内环的限幅。5. 常见问题与排查技巧实录5.1 仿真速度慢到无法忍受Simscape Electrical的模型本身运行就很慢如果不做任何配置仿真1秒可能要跑十几分钟。我的做法是Powergui里选离散求解器步长按PWM开关频率的二十分之一到五十分之一设置。开关频率2kHz时步长50-100μs比较合适。步长太大会丢PWM细节太小则仿真时间成倍增长。另外一个有效手段是降低逆变器开关频率。逆变器开关频率从2kHz降到1kHz对系统级载荷研究结果影响很小但仿真速度能提升将近一倍。如果只是研究供电侧载荷特性可以干脆把逆变器和电机简化成受控电流源省掉整个传动侧的开关模型这样速度还能再快一个量级。5.2 整流器电流环发散新手调四象限整流器最容易遇到电流发散。表现是网侧电流波形变成锯齿状且幅度激增继续跑仿真直接崩掉。原因基本是电流环PI参数给大了或者进线电感设置太小。我的排查思路是先把电压外环断开只用固定占空比的PWM调制看整流器能不能稳住然后确认进线电感值在合理范围用前面那个经验公式算一下最后电流环Kp先给计算值的三分之一Ki先给0跑通了再加积分。这个顺序能帮你快速定位是参数问题还是模型接线问题。5.3 中间直流电压纹波大如果直流电压纹波超过20V优先检查谐振支路是否装配正确。谐振频率要调到100Hz即LC乘积为2.53e-6。拿L0.5mH、C5.07mF的组合来说仿真波形里应该能看到电感电流里有一个明显的100Hz交流分量这个分量正是用来吸收脉动功率的。如果谐振支路参数正确但纹波依旧大检查是不是支撑电容容值太小了。有些情况下电容的等效串联电阻ESR没有设置也会导致高频纹波异常。在Simscape Electrical的电容模块里可以设置串联电阻和并联电导不建议全部设成理想值否则仿真结果会和实际偏离很大。5.4 FFT分析结果异常谐波分析容易在两种情况下出问题一是仿真没进入稳态就去取数据FFT窗口取了开关瞬间的突变二是FFT窗口时间长度不是信号周期的整数倍导致频谱泄漏。我惯用的做法是取仿真最后2秒的数据用Powergui的FFT Analysis Tool把FFT窗口设置成0.1秒也就是5个工频周期基波频率填50Hz最大频率填2000Hz。这样出来的谐波柱状图才可靠。另外FFT分析前要确认波形数据已经稳定可以先把波形放大看一眼确认没有明显的动态过渡过程后再做分析。5.5 常见问题速查表问题现象可能原因处理方案仿真启动就报错缺少Powergui模块在模型任意位置放一个Powergui模块直流电压起不来整流器触发脉冲未使能检查PWM信号是否接入IGBT门极仿真速度极慢求解器用了连续模式Powergui里切换为离散模式步长100μs网侧电流THD超标进线电感偏小加大L按0.15-0.2mH范围调整直流母线电压振荡电压环Kp偏大减小Kp增大Ki带宽降到30Hz电机转速超调大转速环PI参数不当先整定电流环再整定转速环6. 后续扩展方向这套模型搭完之后我已经在上面做了不少工况仿真但还有很多方向可以继续挖。比如多列车同时运行场景下的供电系统协同仿真这种研究对牵引变电所的容量规划和调度策略优化非常有价值。再比如低次谐波电流的补偿策略研究牵引供电系统的谐波治理一直是个热点问题。另外如果对控制算法感兴趣可以试试在四象限整流器上改用模型预测控制MPC替代传统双闭环PI控制动态响应速度会有明显提升但控制器的计算量会翻好几倍仿真步长也需要压得更小。这个方向适合做控制算法研究的读者进一步探索。我个人在实际调试中的体会是牵引供电系统仿真最需要耐心的地方在控制参数整定而整定参数的前提是理解每一级之间的耦合关系。整流器的电压环调不好直流母线波动会影响逆变器输出最终在电机转速上看到奇怪的振荡。不要一上来就想着把整个系统一次调通先分级调试——把整流器单独调稳了再接上逆变器把逆变器调好了再接入电机负载。每加一级都确认波形正常再进入下一级。这个过程虽然慢但能帮你积累很多直觉以后拿到任何一套需要仿真的系统都知道该从哪里下手。