25kV交流牵引供电与JR EH800机车Simulink仿真建模实践

25kV交流铁路牵引供电系统驱动JR EH800列车，这个课题我在接手之前一直觉得是"教科书里的东西"。真正在Simulink里面把供电系统和机车负载完整搭起来之后，才发现里面的坑远比想象中多。从牵引变电所模型怎么简化，到机车功率曲线怎么拟合，再到仿真步长怎么设才能算得快又不发散，每一步都值得单独拿出来讲。

这篇文章我尽量按照"为什么这么做"而不是"做了什么"来写，把我在建模和调试过程中验证过可行的方案、踩过的坑、以及最终确认有效的一些参数设置都整理出来。适合准备做电气化铁路牵引供电仿真、机车-电网耦合分析，或者是Simulink里搞电力电子与电力系统联合仿真的同学参考。我用的版本是MATLAB R2023b的Simulink，配合Simscape Electrical工具箱（也就是之前的SimPowerSystems），这套流程在近几个大版本里都通用。

1. 项目整体设计与思路拆解

1.1 牵引供电系统与JR EH800的核心特征

25kV交流牵引供电系统是当前最主流的大规模电气化铁路供电方案。所谓"25kV"是指牵引变电所二次侧母线电压，实际上送到接触网上的有效值通常维持在27.5kV到29kV之间，这是为了补偿线路压降。到了机车受电弓位置，正常范围一般在22.5kV到31kV之间波动，过压欠压保护都以此为基准。

从结构上看，这个系统主要包含牵引变电所（TSS）、接触网（Catenary）、钢轨回路、回流线，以及沿线分布的分区所、AT所（自耦变压器所）。国内高铁大量采用AT供电方式，即每隔10公里左右设置一台自耦变压器，将接触网和正馈线之间的电压变换为2倍关系，从而降低钢轨电位和线路阻抗，减小对邻近通信线路的干扰。JR EH800在设计上适应的供电制式是交流25kV/50Hz与25kV/60Hz双频段，这是因为它需要兼容东日本和北海道地区不同电网频率的线路。EH800是八轴货运电力机车，整备功率超过8000kW，货列牵引质量可以达到上千吨。对于这样一个"大功率移动负载"，在Simulink里建模时不能当成纯电阻看待，它的整流器、中间直流回路、逆变器和异步牵引电机的动态特性都必须反映出来，否则网压波动和谐波分析都无从谈起。

1.2 为什么选择Simulink作为仿真平台

在Simulink里做这个课题，理由很直接：它支持多物理域联合仿真，电力系统模型可以直接和电机驱动模型放在同一个仿真环境里跑。对比其他方案，比如PSCAD擅长电磁暂态，但对机车内部的电力电子变换器建模不够方便；Matlab/Simulink在这类"供电系统+机车负载"的跨界问题上有天然优势。Simscape Electrical把元件库做得很全，从三相电源、线路、变压器到IGBT、二极管、异步电机都有现成模块，而且支持连续、离散两种求解方式。

我选择搭建的是一种"大电网包小机车"的结构：外部是供电系统模型，负责输出接触网电压；内部是EH800机车的牵引传动链模型，通过受电弓模块与网压连接。这样做的好处是，后续加多列车运行场景时，只需要复制机车模型、修改速度曲线和发车间隔就能实现多车仿真。

1.3 方案选型时的几个决策点

真正动手之前有几个问题需要决策，这里重点说两个。

第一个是仿真域的选型。Simulink里既有连续时间仿真，也有离散时间仿真。牵引供电系统本身是工频系统（50Hz/60Hz），电磁暂态过程到几毫秒量级就有了，但电力电子器件PWM载波频率通常是几百赫兹到几千赫兹，如果步长太大，PWM波形细节完全失真；步长太小，整条线路在10秒仿真里要算几百万步。我最终采用了固定步长离散仿真，步长设为20微秒。这个值既能捕捉到PWM开关过程的主要特征，又不会让仿真时间完全不可接受。

第二个是机车模型的复杂度。EH800是交-直-交传动的典型代表：受电弓引入单相25kV交流电，经过主变压器降压，送入PWM整流器变成直流电，中间直流回路再给逆变器供电，逆变器输出变频变压的三相交流电驱动异步牵引电机。对"载荷研究"来说，如果只关心供电系统侧的有功和无功冲击，可以用一个受控电流源来等效机车；但如果要研究谐波、电压波动和机车内部直流母线电压的稳定性，就必须把整流器和逆变器完整搭出来。我最后选择了"整流器+直流母线+简化逆变器与电机"的方案：整流器是完整的，逆变器和电机用受控转矩负载等效。这样精度不会损失太多，但仿真速度可以快好几倍。

2. 核心细节解析与实操要点

2.1 接触网与回流系统的建模

接触网在电气上是一个沿线分布阻抗的传输结构。在Simulink里，最常用的等效方式是π型集中参数模型。导线电阻大约0.15Ω/km，电感大约1.3mH/km，我按单线20km的供电区间计算，等效串联阻抗取R=3Ω、L=26mH。如果需要反映同一馈线区段内多辆机车同时取流，最好分成多个π型段，而不是一个总集总参数。

这里有个容易被忽略的点：钢轨回流回路的阻抗。很多人建模时直接拿一根理想的接地线当回流路径，实际上钢轨对地有泄漏电阻、钢轨本身也有阻抗，而且钢轨电位是沿线的函数。在AT供电方式下，钢轨电位因为自耦变压器的钳位作用被限制在较低水平，所以我才在前面的方案里坚持用了AT模型。如果不是AT方式，普通直接供电方式下钢轨电位要达到上百伏，这对仿真结果影响很大。

2.2 JR EH800机车交-直-交传动链模型

EH800机车牵引系统的核心是四象限脉冲整流器（4QC）。这个环节是谐波和功率波动的"源头"，也是建模的重点。在Simulink里，我用Universal Bridge模块配置为单相全桥，开关器件选IGBT/Diodes，PWM由载波和调制波比较产生。调制波是电流控制器输出的电压指令，而电流控制器采用瞬态电流控制策略，包含电压外环和电流内环。

直流母线电容的取值很关键，直接关系到中间直流电压的纹波和抗负载冲击能力。EH800中间直流电压在三千伏左右。我按工程经验取电容值为10mF，能够把负载突变时的电压跌落控制在10%以内。

电机侧，我没有把异步电机本体和逆变器完全建模，而是用一个前馈关系：根据当前速度计算出牵引力，再由牵引力和速度折算成功率，将这个功率作为直流侧的扰动负载。这个简化对"载荷研究"足够用，因为载荷研究的核心是供电系统的电压、电流、功率和各节点的能量分配，而不是电机内部的磁链控制。

2.3 载荷曲线的处理

这里说的"载荷"包含两层含义：一层是线路负荷（有多少列机车在取流），另一层是机车自身的牵引载荷（牵引力、速度、坡道附加阻力等）。JR EH800作为货运机车，典型载荷场景是重载起动和持续爬坡。

我用的是经典的阻力公式：F_r = A + B·v + C·v²，其中系数按货运列车经验选取，A对应轴承摩擦阻力，B对应线性阻力，C对应空气阻力。在Simulink里用一个MATLAB Function模块实现，输入是速度，输出是阻力。然后根据加速度需求计算所需牵引力：F_t = M·a + F_r + F_g，其中F_g是坡道阻力，等于M·g·sinθ。

这些计算直接在模型里用增益模块和求和模块搭出来，或者直接写在MATLAB Function里都行。我倾向于在Simulink里把阻力公式可视化成若干子模块，这样后面做参数化扫描时只需要改系数，不用改代码。

3. 实操过程与核心环节实现

3.1 顶层模型搭建步骤

先搭一个干净的顶层模型，我建议分三步走，不要一上来就把所有模块拖进去。

第一步，搭供电系统。放一个Three-Phase Source（实际上这里应该是单相/两相考虑，需根据实际线路接法选择），因为牵引供电是单相负荷，我习惯用两个单相交流源模拟V/V接线或Scott接线的牵引变电所输出。具体接线方式根据手头数据来，但简化建模时可以直接用一个可编程交流电压源（AC Voltage Source）代表接触网端的开路电压，然后串联线路阻抗。如果要体现AT供电，则需要加入一个理想变压器模型，变比1:1，一次侧并接在接触网和正馈线之间，中间抽头接钢轨。

第二步，搭机车块。机车模型的输入是接触网电压，输出是取流电流。受电弓处我用了一个Breaker模块来控制机车是否投入，这样后面可以设置成"列车进入供电区段后才接通"，模拟列车从无电区进入带电区的过程。

第三步，搭测量与显示单元。在接触网末端（机车位置）放电压测量模块，在机车支路放电流测量模块，功率测量直接用Simscape的Power Sensor或者自己用电压电流乘积计算。示波器Scope建议用Sample based模式，配合To Workspace模块把关键信号存到工作区，方便后续做FFT分析谐波。

3.2 整流器参数与控制器整定

很多人卡在四象限整流器这里。如果你也遇到"模型跑起来了但直流母线电压稳不住"的问题，大概率是控制器参数没整定好。

以我用的瞬态电流控制为例，电压外环采用PI控制器，电流内环采用比例控制器加负载电流前馈。控制器输出的调制信号和三角载波比较，生成IGBT的PWM脉冲。三角载波频率我设为1250Hz，这是沿用实际机车变流器的常用水平。频率太低，谐波明显；频率太高，仿真时间陡增。

PI参数的整定我推荐一个比较基础的做法：先把电流内环的比例系数拉开，让电流能跟踪电压指令；再调电压外环，观察直流母线响应。仿真步长20微秒下，电压环带宽取100Hz左右，对应PI参数大约Kp=2、Ki=50；电流内环带宽取500Hz左右，比例系数大约Kp=5。当然这些数值要结合具体标幺值系统调整，直接照搬可能不收敛。

3.3 三组典型工况仿真结果

我做了三组典型工况：空载、满载牵引、电气制动回馈。

空载工况下，机车只带辅助负载，接触网电压几乎不跌落，接触网电流幅值很小，主要是无功励磁电流和谐波电流。这时候主要观察稳态电压波形和THD。

满载牵引工况是最关键的。机车在10秒内从0加速到100km/h，最大牵引功率达到8000kW。网压从27.5kV跌落到约25.8kV，原因是我把线路阻抗和变压器漏抗都计算进了压降中。接触网电流峰值超过350A，直流母线电压维持在3000V附近，纹波约为2.5%。这个结果符合工程经验范围。

电气制动回馈工况下，机车相当于一个发电单元，把动能回馈到电网。此时接触网电流相位反向，有功功率变为负值。这个工况对验证整流器四象限能力和供电系统能量吸收能力很有价值。仿真中观察到直流母线电压上升，能量回馈引起接触网末端电压小幅抬升，整个系统仍保持稳定。

3.4 仿真加速与数据管理技巧

如果觉得几千瓦功率、几十秒运行时间太长，有几个实用技巧可以提速。第一，把PWM开关管的开关频率降低到可接受范围，比如从1250Hz降为1050Hz，仿真时间能缩短不少。第二，把模型改为离散仿真并使用Dormand-Prince求解器（ode45的离散版），也可以显著减少计算量。第三，可以把机车模型在第一个工况稳定后，直接用一个查表模块代替整流器部分，这样后续多工况只需修改查表数据。

数据管理上，我习惯将每次仿真后导出的Voltage、Current、Power等信号数组，保存为一个结构体，文件名包含工况标签和日期。这样比较不同参数设置下的结果时，直接用MATLAB脚本加载批量绘图即可。

4. 常见问题与排查技巧实录

4.1 模型不收敛怎么排查

Simulink里电气仿真最常见的问题就是"Nonlinear iteration is not converging"。我的排查顺序固定如下：先检查模型是否有代数环，特别是直接电压源串联电感再接电流源的情况，这种情况下模型可能出现数值刚性。其次是检查步长，如果模型在连续时间下跑模拟量，把步长改小一个量级看是否收敛；如果是离散模型，检查是否选择了zero-crossing detection选项，有时候关掉过零检测反而更快收敛。

还有一个容易忽略的，是Simscape Electrical模块之间的物理信号和Simulink普通信号混用。比如你把电压测量输出直接接到一个普通的Gain模块，再连接到可控电流源的输入，这里就可能产生代数环。解决办法是在连接处放入Memory或Unit Delay模块，切断代数环，或者用Physical Signal to Simulink Converter做好类型转换。

4.2 直流母线电压振荡问题

我调试过程中遇到过一次直流母线电压持续振荡的情况，频率大约是50Hz，幅值约为目标电压的15%。排查下来，原因是电压外环PI参数过强，导致在工频周期内产生了次同步振荡，电流内环无法抑制。解决方法就是降低电压环PI参数，同时提高电流内环的比例增益，让电流环响应更快，把电压环的校正需求降到最低。

此外还要检查负载功率突变点。EH800的牵引力从零到满值只有几秒钟，中间直流母线电流会有剧烈的阶跃变化。如果母线电容偏小，电压纹波会更加明显。我试过把电容从10mF降到5mF，纹波率从2.5%上升到7%，效果很明显，说明至少要保证1mF/100kW这个量级。

4.3 FFT谐波分析时的窗口选择

在做接触网电流谐波分析时，FFT窗口长度要用整周期，否则频谱泄漏会把谐波成分弄得一团糟。50Hz系统下，我通常取0.2秒的数据（10个周波），并加Hanning窗。重点关注3、5、7、9次谐波。对于PWM整流器，其主要谐波集中在开关频率整数倍附近，也就是1250Hz±50Hz，以及2500Hz±50Hz附近。Simulink的Powergui模块里有FFT Analysis工具，加载To Workspace导出的波形后可以直接选周期分析。

一个关键点：在使用FFT之前要把波形数据的采样时间统一。因为Simulink里不同模块可能采用不同的采样率，To Workspace导出时如果不规定好，示波器看到的可能是混叠后的结果。我通常把示波器采样率设为仿真步长的整数倍，比如0.2ms，这样FFT窗口频率分辨率是5Hz，足够分辨主要谐波分量。

4.4 不同频率系统（50Hz/60Hz）兼容性建模

EH800需要同时兼容50Hz和60Hz，这让仿真建模多了一件事。在Simulink里，我通过一个可变频率的电压源模块实现频率切换。频率设为50Hz时，系统周期为20ms，载波比是1250/50=25；频率设为60Hz时，周期为16.67ms，载波比是1250/60≈20.83，PWM谐波分布会有所不同。建议把载波频率随系统频率同步调整，保持载波比为整数，否则谐波频谱分析会变得很复杂。

另外，在AT供电方式下，正馈线和接触网之间的电压幅值实际上是1:2关系。简化建模时，可以在正馈线端不单独建模，把自耦变压器的电压变化等效为接触网对钢轨电压的2倍关系，这样模型结构会简单很多。但要注意，这种方式无法反映钢轨电位，如果需要研究人身安全和信号干扰，必须按实际AT结构建模。

4.5 仿真时间与精度的取舍建议

关于仿真步长，我给一个比较实用的参考区间。如果只研究稳态潮流和谐波，步长可以放宽到50微秒，仿真很快。如果研究机车起动时的暂态电压跌落、直流母线动态响应，步长最好在10~20微秒，否则动态过程可能被数值滤波抹平。如果还要研究线路保护动作或故障过程，比如接触网短路，建议步长缩到5微秒以下，不过这样的仿真时长会急剧增加，要有心理准备。

我自己的经验是：主要研究目标定为"载荷"而不是"故障"，所以20微秒的固定步长足够，单次仿真10秒的运行时间大约在20分钟到1小时之间，视电脑配置而定。如果对某一个工况需要大规模参数扫描，建议把每次扫描的仿真时间缩短到5秒以内，结果差异通常可以忽略。

4.6 常见问题速查表

下表罗列了我建模过程中遇到的高频问题、可能原因和解决方案，对照参考可以省去不少重复排查时间。

5. 从建模仿真到课题结论

整个仿真跑下来，我对25kV交流牵引供电系统和JR EH800这种大功率货运机车的匹配特性有了更直观的理解。最核心的结论是：8MW级别的机车负载对25kV供电网来说是一个相当可观的扰动源，线路阻抗、供电方式、AT所分布直接决定了网压波动幅度。采用AT供电方式能够明显改善电压水平，这个在仿真中可以通过对比"有AT"和"无AT"两种模型直观看到。

另一个重要结论是，四象限整流器对电网电压的支撑能力有限。当直流母线负载突增时，如果不加限制，母线电压可能跌到规定范围以下，而此时整流器能提供的最大电流是有上限的。这个现象在仿真里很容易暴露出来，如果只做稳态核算，很难发现。

最后多说一句，仿真模型的真实性完全取决于参数的边界条件。我基于公开资料和经验数据搭建的模型，适合做趋势性研究和方案对比。如果要用在工程设计和产品开发上，建议根据实际线路数据、实际机车参数重新标定模型参数。MATLAB的Parameter Estimation工具箱可以用来做基于实测数据的参数辨识，这是把仿真模型从"大致对"变成"比较准"的一条可行路径。