第四节:车辆纵向动力学模型

各位工程师朋友,今天我们来聊聊车辆纵向动力学模型。说实话,这个模型是整车动力学仿真中最基础、也最常用的一块。我刚开始做车辆动力学仿真时,总觉得纵向模型太简单,不就是个牛顿第二定律吗?直到有一次在实车标定中,因为忽略了滚动阻力的非线性特性,导致仿真结果和实测数据差了将近15%。嗯,从那以后,我再也不敢小看这个"简单"模型了。

4.1 纵向力平衡方程

车辆纵向动力学,说白了就是研究车在前进方向上受的力。根据牛顿第二定律,我们可以写出最基本的平衡方程:

m * dv/dt = F_tractive - F_resistance

其中:

  • m — 整车质量(kg),包括乘客和货物
  • dv/dt — 纵向加速度(m/s²)
  • F_tractive — 驱动力(N),来自发动机/电机经传动系统传递到车轮
  • F_resistance — 总行驶阻力(N)

总行驶阻力又可以拆分为四个部分:

F_resistance = F_air + F_rolling + F_grade + F_accel

我个人习惯把阻力分成"跑起来就有的"和"特定工况才有的"两类。空气阻力和滚动阻力是常驻的,坡度阻力和加速阻力则看路况和驾驶行为。

4.2 空气阻力建模

空气阻力,也叫气动阻力。你想想看,车在空气中高速行驶,就像你把手伸出车窗一样,风会推着你的手往后。这个力跟车速的平方成正比。

标准公式:

F_air = 0.5 * ρ * Cd * A * v²

参数说明:

符号 含义 典型值 单位
ρ 空气密度 1.225(海平面,20°C) kg/m³
Cd 风阻系数 0.25~0.35(轿车) 无量纲
A 迎风面积 2.0~2.5(轿车)
v 相对车速 m/s
我的经验: 我在做某款电动SUV的续航仿真时,发现空气阻力在80km/h以上时占比急剧上升。高速工况下,空气阻力能占到总阻力的60%以上。所以电动车为什么都拼命降风阻?原因就在这里。

这里有个细节要注意:公式中的v是相对车速,不是绝对车速。如果有风,要加上风速分量。我见过有人直接用GPS车速算空气阻力,结果在逆风工况下仿真偏差很大。

4.3 滚动阻力建模

滚动阻力来自轮胎与地面的接触变形。轮胎不是刚体,压在地面上会变形,这个变形会产生能量损耗。

基本公式:

F_rolling = f * m * g * cos(θ)

其中:

  • f — 滚动阻力系数,无量纲
  • m — 整车质量(kg)
  • g — 重力加速度(9.81 m/s²)
  • θ — 道路坡度角(rad)

滚动阻力系数f不是常数,它跟车速、胎压、路面材料都有关系。我常用的经验公式是:

f = f0 + f1 * (v / 100) + f4 * (v / 100)^4

典型参数:

路面类型 f0 f1 f4
良好沥青路 0.008 0.0012 0.0002
混凝土路 0.010 0.0015 0.0003
碎石路 0.015 0.0020 0.0005
避坑指南: 我曾经在低温环境下做测试,发现滚动阻力比常温大了将近30%。后来查资料才知道,低温导致轮胎橡胶变硬,滞后损失增大。所以冬季续航打折,滚动阻力增加也是重要原因之一。

4.4 坡度阻力建模

坡度阻力,就是车在上坡时重力沿坡道方向的分量。这个力只跟坡度和车重有关,跟车速无关。

公式很简单:

F_grade = m * g * sin(θ)

上坡时F_grade为正(阻力),下坡时为负(助力)。

实际道路的坡度通常用百分比表示,比如5%的坡度意味着每100米水平距离上升5米。换算关系:

坡度百分比 = tan(θ) * 100%
θ = arctan(坡度百分比 / 100)

举个例子:5%的坡度,θ ≈ 2.86°。对于一辆2吨的车,坡度阻力约为:

F_grade = 2000 * 9.81 * sin(2.86°) ≈ 980 N

这个力不小了。我做过一个测试,在8%的坡道上,一辆2.5吨的SUV需要额外输出将近2000N的力才能维持匀速。所以山区道路的能耗比平原高很多,原因就在这里。

4.5 加速阻力建模

加速阻力,严格来说应该叫"加速惯性力"。根据达朗贝尔原理,我们把加速时产生的惯性力当作一种阻力来处理。

公式:

F_accel = δ * m * a

其中:

  • a — 纵向加速度(m/s²)
  • δ — 旋转质量换算系数,大于1

为什么要引入δ?因为车辆加速时,不仅整车质量要加速,发动机飞轮、变速箱齿轮、传动轴、车轮等旋转部件也要加速旋转。这些旋转部件的转动惯量会"消耗"一部分驱动力。

δ的估算公式:

δ = 1 + (I_wheel / (m * r²)) + (I_engine * i² * η) / (m * r²)

简化版经验值:

挡位 δ典型值
1挡 1.3~1.5
2挡 1.2~1.3
3挡 1.1~1.2
4挡及以上 1.05~1.1
我的习惯: 在做初步仿真时,我一般取δ=1.1作为常数。但如果是做精确的加速性能分析,我会根据挡位实时计算δ。特别是电动车,电机转动惯量小,δ比燃油车要低一些。

4.6 完整纵向动力学模型

把上面四个阻力合起来,完整的纵向动力学方程就是:

m * dv/dt = F_tractive - [0.5*ρ*Cd*A*v² + f*m*g*cos(θ) + m*g*sin(θ) + δ*m*a]

这个方程,就是整车纵向动力学仿真的核心。不管是做燃油车的动力性分析,还是电动车的续航仿真,都离不开它。

下面我用SVG画一张图,把整个知识体系串起来:

车辆纵向动力学模型知识体系 纵向力平衡方程 m·dv/dt = F_tractive - F_resistance 空气阻力 F_air = 0.5·ρ·Cd·A·v² 关键参数:Cd, A, ρ 与车速平方成正比 滚动阻力 F_roll = f·m·g·cos(θ) 关键参数:f, m f随车速非线性变化 坡度阻力 F_grade = m·g·sin(θ) 关键参数:θ, m 与车速无关 加速阻力 F_accel = δ·m·a 关键参数:δ, a δ考虑旋转质量效应 总阻力 = 空气阻力 + 滚动阻力 + 坡度阻力 + 加速阻力

这张图把四个阻力的关系理清楚了。你看,空气阻力和滚动阻力是"常驻"的,只要车在跑就有。坡度阻力看路况,加速阻力看驾驶行为。实际仿真时,这四个力要同时计算,然后代入平衡方程求解加速度。

实战要点总结:

  • 空气阻力:高速工况主导,注意相对风速
  • 滚动阻力:低速工况占比大,注意温度和胎压影响
  • 坡度阻力:山区道路不可忽略,用百分比坡度换算
  • 加速阻力:起步和急加速时显著,别忘了旋转质量系数δ

好了,这一节的内容就到这里。纵向动力学模型虽然基础,但它是整车仿真的基石。我建议你在做任何车辆动力学项目之前,先把这四个阻力的参数标定准确。参数不准,后面的仿真都是空中楼阁。


公众号:蓝海资料掘金营,微信deep3321