2. 车辆动力学模型：单轮模型与运动方程推导

图中几个关键量：车体质量 m、车轮转动惯量 J、轮胎有效滚动半径 R。制动时，制动力矩 T_b 作用在车轮上，地面制动力 F_x 让车减速。

2.2 车辆运动方程（车体平动）

根据牛顿第二定律，车体平动方程很简单：

m · dv/dt = -Fx

这里 v 是车速，F_x 是地面制动力。负号表示制动力让车速下降。

那 F_x 怎么来的？它取决于轮胎与路面的附着系数 μ 和垂向载荷 F_z：

Fx = μ(λ) · Fz

其中 λ 就是滑移率。你看，滑移率通过附着系数 μ 直接影响制动力，进而影响车速——这就是为什么滑移率估算如此重要。

2.3 车轮旋转动力学方程

车轮转起来，得用转动定律。绕车轮中心的力矩平衡：

J · dω/dt = -Tb + Fx · R

解释一下：

• J · dω/dt —— 车轮角加速度，转动惯量乘角加速度
• -T_b —— 制动力矩，让车轮减速（负号）
• F_x · R —— 地面制动力产生的力矩，让车轮继续转（正号）

制动时，T_b 大于 F_x·R，所以车轮角速度下降。如果 T_b 太大，车轮可能抱死——ω 变成 0，但车还在滑行。

2.4 两个方程联立，构成完整模型

把平动和转动方程放一起：

这三个方程就是单轮模型的全部。你想想看，输入只有 T_b（制动压力调节的输出），状态变量是 v 和 ω，通过滑移率 λ 把两者耦合起来。

2.5 模型怎么用？一个仿真例子

我习惯在 MATLAB/Simulink 里搭这个模型。代码很简单：

% 单轮模型参数
m   = 400;      % 四分之一车重 [kg]
J   = 1.2;      % 车轮等效转动惯量 [kg·m²]
R   = 0.32;     % 轮胎有效滚动半径 [m]
Fz  = m * 9.81; % 垂向载荷 [N]

% 初始状态
v0  = 30;       % 初始车速 [m/s] ≈ 108 km/h
w0  = v0 / R;   % 初始轮速 [rad/s]

% 仿真参数
dt  = 0.001;    % 步长 1ms
t_end = 3;      % 仿真时长 3s

% 制动压力输入（阶跃）
Tb = 1200;      % 制动力矩 [Nm]

% 主循环（简化示意）
for k = 1:t_end/dt
    % 计算滑移率
    lambda = (v - w*R) / max(v, 0.01);
    
    % 查附着系数曲线（魔术公式简化）
    mu = mu_lambda(lambda);
    
    % 地面制动力
    Fx = mu * Fz;
    
    % 更新车速
    v = v - (Fx/m) * dt;
    
    % 更新轮速
    w = w + (-Tb + Fx*R)/J * dt;
end

2.6 模型的局限与工程取舍

单轮模型虽然简单，但有几个地方你得心里有数：

• 忽略了载荷转移 —— 制动时前后轮垂向力会变，单轮模型没法体现
• 没有侧向动力学 —— 转弯制动时，侧向力会影响纵向附着
• 轮胎模型简化 —— 实际 μ-λ 曲线受路面、温度、磨损影响很大

但话说回来，做滑移率估算和制动压力调节的算法验证，单轮模型是业界标配。我参与过的三个量产项目，算法前期验证全是用这个模型。等算法成熟了，再往七自由度甚至整车模型里移植。