4. 列车运行仿真模型构建入门

要分析CBTC系统的运营能力，仅靠静态的间隔公式往往不够，因为列车在区间内的速度变化、制动响应和前后车追踪关系都是动态过程。本节从移动闭塞原理出发，搭建一个可运行的列车运行仿真模型骨架，包含列车动力学、制动模型与安全防护距离三部分，并用一个简单区间案例验证模型输出是否合理。

一、坐标系与建模约定

建模前必须先统一参照约定，否则后续计算极易出现符号错误。本节采用如下约定：以线路起点为里程零点，列车位置 s 沿运行方向递增，单位为米；速度 v 沿运行方向为正，单位为 m/s；加速度 a 以牵引为正、制动为负。仿真步长取 dt = 0.1 s，采用欧拉法逐步积分：v(k+1) = v(k) + a·dt，s(k+1) = s(k) + v(k)·dt。步长越小精度越高，但计算量增大，0.1 s 对运营能力分析已足够。

二、列车动力学参数与牵引制动特性

列车纵向运动遵循牛顿第二定律，单位质量合力为：a = (F_t - F_b - R(v)) / M，其中 F_t 为牵引力，F_b 为制动力，R(v) 为基本阻力，M 为列车质量。基本阻力常用戴维斯公式描述：R(v) = A + B·v + C·v²，A 反映滚动与轴承阻力，B 反映机械摩擦，C 反映空气阻力。以某城市轨道交通 6 辆编组列车为例（示例假设参数）：M = 3.2×10⁵ kg，A = 2.5 kN，B = 0.06 kN·s/m，C = 0.0035 kN·s²/m²。牵引力与制动力并非恒定，通常随速度升高而衰减，工程上可用分段折线近似：低速段保持恒牵引力（如 300 kN），达到转折速度后按恒功率（如 4 MW）反比衰减，即 F = P/v。

实现时，把牵引特性写成速度的函数表，仿真每一步根据当前速度插值取力，再叠加阻力求加速度。这样列车在区间的加速、惰行、制动全过程都能被逐步还原，而不是用一个平均速度粗略代替。

三、制动模型与安全防护距离

CBTC 采用移动闭塞，后车追踪的目标是前车尾部加上一段安全防护距离。制动模型通常采用分段模型：先是系统反应时间（含通信延迟与车载设备响应，示例取 t_r = 2.0 s），此阶段列车仍以原速走行；随后制动率线性上升至最大常用制动率（示例取 β = 1.0 m/s²）。因此后车从当前位置到完全停住所需距离为：

d_brake = v·t_r + v² / (2·β_eff)
d_safe  = d_brake + d_margin

其中 β_eff 为考虑制动率爬升后的等效减速度，d_margin 为安全余量（示例取 50 m，覆盖测距误差与停车点偏差）。移动闭塞下，后车的移动授权终点即为前车尾部位置减去 d_safe，后车据此生成制动曲线。需要强调的是，d_safe 随速度平方增长，这正是高速段追踪间隔难以压缩的根本原因。

四、简单区间案例验证

设区间长度 1200 m，前车在 t=0 时位于 1000 m 处并以 15 m/s 匀速运行，后车从 0 m 处静止起动，采用最大牵引加速至 20 m/s 后匀速，接近授权终点时按制动模型减速。仿真中每一步计算后车移动授权终点 s_auth = s_lead - d_safe(v)，若 s_rear + v·dt > s_auth 则触发制动。运行后可观察三个关键量：后车是否在授权终点前停住、最小追踪间距是否始终大于零、以及后车从起动到停稳的总时间。

按上述示例参数估算：后车以 20 m/s 运行时，d_brake = 20×2.0 + 400/(2×1.0) = 40 + 200 = 240 m，加 50 m 余量得 d_safe = 290 m。若前车尾部在 1000 m，则授权终点约在 710 m，后车须在 710 m 前完成制动。若仿真结果显示后车越过该点，说明制动触发逻辑或参数设置有误，应检查反应时间、制动率与步长是否匹配。反之，若后车过早制动、追踪间距远大于 d_safe，则说明触发阈值过于保守，会人为拉大追踪间隔，影响能力评估的真实性。

五、建模要点小结

本节模型的核心是把“移动闭塞”翻译成可计算的约束：前车位置决定授权终点，授权终点减去安全防护距离得到后车可用的制动起点，列车动力学决定后车能否在该点前停住。参数取值应来自车辆与信号系统的实际特性，示例中的数值仅用于教学演示。模型跑通后，就可以通过改变站间距、列车性能、反应时间等参数，观察追踪间隔的变化，为后续的行车间隔优化提供定量基础。