第二节 基础物理模型：列车动力学模型、基本阻力与附加阻力、牵引力与制动力特性

各位同学，咱们今天聊聊列车运行的“老本行”——物理模型。说白了，就是搞清楚列车到底是怎么动起来的，又是什么在拖它的后腿。

我个人习惯，做任何节能曲线设计之前，先把物理底子摸透。你想想看，连车怎么跑、阻力从哪来都搞不清楚，后面谈什么优化节能？那不成了空中楼阁嘛。

2.1 列车动力学模型

列车跑起来，本质上就是一个质点的运动。嗯，这里要注意，虽然列车很长，但在宏观运行分析中，我们通常把它简化成一个单质点模型。

核心公式其实就一个——牛顿第二定律的变体：

m · a = F_t - F_b - R(v) - G(x)

其中：

• m —— 列车质量（含旋转质量系数，我一般取1.06~1.08）
• a —— 加速度
• F_t —— 牵引力
• F_b —— 制动力
• R(v) —— 基本阻力（速度的函数）
• G(x) —— 附加阻力（坡道、曲线等，位置的函数）

我在项目中遇到过一件事：有同事直接用列车总质量代入计算，结果仿真出来的加速曲线总是偏快。后来一查，忘了考虑旋转质量系数。别小看这6%~8%，累积下来误差不小。

2.2 基本阻力与附加阻力

阻力这东西，看不见摸不着，但实实在在影响能耗。我把它分成两类：

2.2.1 基本阻力

基本阻力是列车在平直轨道上运行时始终存在的阻力。它跟速度关系密切，通常用戴维斯公式拟合：

R_b(v) = A + B·v + C·v²

参数含义：

参数	物理意义	典型值（地铁列车）
A	机械阻力（轴承、车轮等）	0.5~1.5 kN
B	与速度成正比的阻力	0.01~0.03 kN/(km/h)
C	空气阻力系数	0.0002~0.0005 kN/(km/h)²

我曾经吃过一次亏。在某条线上做节能测试，仿真结果和实测数据总是对不上。排查了三天，最后发现是空气阻力系数C取值偏小了。那条线路隧道断面小，活塞效应明显，空气阻力比标准值大了将近30%。

2.2.2 附加阻力

附加阻力是“额外”的，遇到特定线路条件才会出现。主要有三种：

• 坡道阻力： R_g = m·g·sin(θ)。上坡为正阻力，下坡为负阻力（其实是助力）。
• 曲线阻力： 经验公式 R_c = 600 / R（R为曲线半径，单位米）。半径越小，阻力越大。
• 隧道阻力： 跟隧道断面、长度、列车外形都有关。一般取基本阻力的10%~30%。

你想想看，如果一条线路既有大坡道又有小半径曲线，附加阻力可能比基本阻力还大。这时候节能策略就得调整——比如上坡前提前加速冲坡，下坡时充分利用重力势能。

2.3 牵引力与制动力特性

牵引力和制动力，是列车运行的“油门”和“刹车”。它们的特性直接决定了列车能跑多快、能停多稳。

2.3.1 牵引力特性

牵引力不是恒定的。它随速度变化，典型特性分三段：

1. 恒力矩区（低速段）： 牵引力恒定，约0~10 km/h。电机输出最大转矩。
2. 恒功率区（中速段）： 牵引力随速度增加而下降，功率恒定。约10~40 km/h。
3. 自然特性区（高速段）： 牵引力继续下降，电机特性限制。约40 km/h以上。

用代码表示大概是这样：

def calc_tractive_force(v, F_max, P_max):
    if v < v_base:
        return F_max  # 恒力矩区
    else:
        return P_max / v  # 恒功率区

嗯，这里要注意，实际牵引特性曲线不是这么理想化的。电机有弱磁控制、有温升限制，实际输出会略有不同。我建议做高精度仿真时，直接用厂家提供的牵引特性曲线数据表，别自己拟合。

2.3.2 制动力特性

制动力分两种：

• 电制动（再生制动）： 电机反转，把动能转化为电能回馈电网。节能！但低速时制动力会衰减。
• 空气制动（摩擦制动）： 闸瓦抱紧车轮。可靠！但磨损大，有粉尘。

现代列车通常采用“电制动优先、空气制动补充”的策略。我参与过一个项目，为了最大化节能，把电制动的退出速度从10 km/h降到了5 km/h。结果确实多回收了一些能量，但低速段电制动扭矩波动大，乘客舒适性受了影响。后来还是调回了8 km/h——节能和舒适之间得找个平衡。

好了，物理模型这块就讲这么多。下一节咱们聊聊怎么把这些模型用起来，构建节能优化的目标函数。有什么问题随时问我。