3. 车辆动力学基础与底盘控制模型

底盘控制算法的原型设计,本质上是在“车辆会怎么动”和“控制器该发什么指令”之间建立可计算的桥梁。纵向、侧向、垂向三个方向的动力学方程,就是这座桥梁的三根主梁。它们分别回答三个问题:车速如何被驱动与制动改变?航向如何被转向与横摆改变?车身姿态如何被悬架与载荷转移改变?而轮胎力、载荷转移、滑移率则是贯穿三者的关键变量——控制器输出的扭矩、制动压力、转向角、阻尼力,最终都要通过轮胎与地面的接触印迹转化为车辆运动。

一、坐标系与变量约定

为避免符号混乱,本节统一采用ISO 8855车辆坐标系:x轴指向车辆前方,y轴指向驾驶员左侧,z轴垂直向上,绕z轴的横摆角速度ψ̇以逆时针为正。纵向加速度ax以前进为正,侧向加速度ay以向左为正。轮胎纵向力Fx沿轮胎平面向前为正,侧向力Fy以指向轮胎左侧为正。滑移率κ定义为κ=(ωR−vx)/max(|vx|,ε),驱动时为正、制动时为负;轮胎侧偏角α以轮胎速度方向相对轮胎朝向的夹角表示。全文与图示均遵循这一约定。

二、纵向动力学:驱动、制动与滑移率

纵向方程可写为:m·ax = Fx,f + Fx,r − Fdrag − Froll − m·g·sinθ。其中Fx,f、Fx,r分别为前、后轴轮胎纵向力之和,Fdrag=0.5ρCdAvx2为空气阻力,Froll=f·m·g·cosθ为滚动阻力,θ为道路坡度。控制器能直接改变的是驱动扭矩Td与制动压力pb,它们通过轮胎半径R和制动器增益Kb转化为轮端力矩,再经轮胎滑移特性变成Fx。

轮胎纵向力与滑移率的关系呈“先升后降”的曲线:小滑移区(约0~0.15)近似线性,Fx≈Cκ·κ,Cκ为纵向滑移刚度;峰值附着系数通常出现在κ≈0.1~0.2;超过峰值后进入不稳定区,滑移率继续增大而Fx下降。ABS与TCS的核心任务,就是把κ控制在峰值附近。例如某试验车在干沥青路面、垂直载荷4 kN条件下,Cκ约为80 kN/单位滑移率,峰值Fx约4.4 kN,对应κ≈0.12。控制器若采用滑模或PID调节制动压力,目标就是让实测κ跟踪这一参考值。

三、侧向动力学:横摆与侧偏

侧向运动常用二自由度自行车模型描述:m(v̇y+vxψ̇)=Fy,f+Fy,r,Izψ̈=a·Fy,f−b·Fy,r。其中a、b为质心到前、后轴距离,Iz为横摆转动惯量。轮胎侧向力在小侧偏角下近似Fy=Cα·α,Cα为侧偏刚度。前轮侧偏角αf=δ−arctan((vy+aψ̇)/vx),后轮αr=−arctan((vy−bψ̇)/vx)。

把线性轮胎模型代入,可得到横摆角速度对前轮转角的稳态增益:ψ̇/δ = vx/(L+Kusvx2),其中L=a+b,Kus=(m/L)(b/Cα,f−a/Cα,r)为不足转向梯度。Kus>0表示不足转向,Kus<0表示过多转向。ESC的参考横摆角速度通常由该式乘以一个安全系数得到,再与实际ψ̇比较,差值超过门限即触发单轮制动。例如某轿车m=1500 kg,a=1.2 m,b=1.4 m,Cα,f=80 kN/rad,Cα,r=90 kN/rad,可算得Kus≈0.0027 rad/(m/s²),在vx=20 m/s时稳态增益约0.31 (rad/s)/rad。

四、垂向动力学与载荷转移

垂向方程描述簧上质量运动:ms·z̈s=Fsusp−msg,悬架力Fsusp=ks(zu−zs)+cs(żu−żs)。更关键的是载荷转移:纵向加速度引起前后轴垂向载荷变化ΔFz,f=−m·ax·h/L,ΔFz,r=+m·ax·h/L;侧向加速度引起左右轮载荷转移ΔFz=m·ay·h/tw,h为质心高度,tw为轮距。载荷转移直接改变轮胎可用附着力:Fx,max=μ·Fz,因此制动时前轴载荷增加、后轴减少,若制动力分配不当,后轮先抱死会导致失稳。

悬架控制(如CDC连续阻尼控制)通过调节cs来抑制车身俯仰与侧倾,但也会影响轮胎动载荷。控制器设计时需在“车身姿态舒适”与“轮胎接地性”之间折中:阻尼过大,轮胎动载荷波动加剧,附着利用率下降;阻尼过小,车身姿态恶化,载荷转移更剧烈。

五、三向耦合与算法影响

三个方向并非独立:纵向滑移率影响侧向力(摩擦椭圆),侧向加速度引起载荷转移进而改变纵向力上限,垂向动载荷又反过来影响滑移率控制精度。因此底盘控制算法原型通常采用分层结构:上层做运动规划与参考模型,中层做纵向/侧向/垂向协调分配,下层做执行器控制。理解这些方程与关键变量的耦合关系,才能在原型阶段合理简化模型、设定控制目标与门限。

车辆三向动力学耦合关系与关键变量 纵向动力学 m·aₓ = Fₓ,f + Fₓ,r − F_drag − F_roll Fₓ = C_κ · κ(小滑移区) 峰值附着:κ ≈ 0.1~0.2 控制目标:ABS/TCS 跟踪 κ_ref 侧向动力学 m(v̇_y + vₓψ̇) = F_y,f + F_y,r I_z ψ̈ = a·F_y,f − b·F_y,r F_y = C_α · α(小侧偏角) 控制目标:ESC 跟踪 ψ̇_ref 垂向动力学 m_s·z̈_s = F_susp − m_s·g ΔF_z,f = −m·aₓ·h/L ΔF_z = m·a_y·h/t_w 控制目标:CDC 调节 c_s 摩擦椭圆 载荷转移 F_z 变化 → 附着上限变化 → 滑移率控制精度 轮胎力生成与执行器接口 驱动扭矩 T_d → 轮端力矩 → F_x(κ) | 制动压力 p_b → F_x(κ) 转向角 δ → α_f, α_r → F_y(α) | 阻尼力 c_s → F_susp → F_z 底盘控制算法原型(分层协调) 参考模型 → 纵向/侧向/垂向协调分配 → 执行器控制 关键变量: κ, α, F_z, aₓ, a_y, ψ̇ 耦合风险: 滑移率↑ → 侧向力↓

六、从方程到原型:一个可执行的建模检查清单

在算法原型阶段,建议按以下顺序把上述方程落地为可运行模型,避免一开始就陷入高保真仿真:

步骤 建模对象 关键输入 关键输出 验证方式
1 纵向单轮模型 T_d, p_b, μ, F_z κ, F_x, v_x 阶跃扭矩下κ响应是否收敛
2 二自由度侧向模型 δ, v_x, C_α ψ̇, v_y, a_y 稳态增益与理论值对比
3 载荷转移模块 a_x, a_y, h, L, t_w F_z,f, F_z,r, F_z,l, F_z,r 制动/转向工况下F_z变化方向
4 摩擦椭圆耦合 F_x, F_y, μ, F_z 可用F_y上限 联合工况下侧向力衰减
5 控制器接口 参考ψ̇, κ_ref 制动压力/扭矩指令 跟踪误差与执行器饱和

以第1步为例,若取m=1500 kg、R=0.32 m、C_κ=80 kN、F_z=4 kN,初始v_x=20 m/s,在t=0时施加T_d=800 N·m的阶跃扭矩,可观察到κ从0快速上升,约0.3 s后接近0.12,F_x达到约4.4 kN,随后若继续增大扭矩,κ越过峰值进入不稳定区,F_x反而下降。这一现象正是TCS需要快速干预的依据。原型中可用一阶惯性环节近似轮速动态,时间常数取0.05~0.1 s,便于控制器参数整定。

工程提示:在原型阶段,纵向与侧向模型可先解耦验证,但进入联合工况(如弯道制动)前,必须引入摩擦椭圆约束:√((F_x/(μF_z))² + (F_y/(μF_z))²) ≤ 1。否则控制器可能在纵向请求过大时误判侧向可用力,导致ESC介入滞后。
常见误区:把载荷转移仅当作“悬架问题”而忽略其对轮胎力的影响。实际上,制动时前轴载荷增加约m·a_x·h/L,若h=0.55 m、L=2.6 m、a_x=0.6g,前轴载荷增加约12%,后轴减少约12%,后轮附着上限同步下降,制动力分配若仍按静态比例,后轮极易先抱死。

综上,纵向、侧向、垂向方程分别给出了车速、航向、姿态的演化规律,而轮胎力、载荷转移、滑移率是连接方程与执行器的三个关键变量。原型阶段不必追求轮胎模型的高阶非线性,但必须保证滑移率-力曲线、侧偏角-力曲线、载荷转移方向与摩擦椭圆约束的正确性,否则后续控制器参数整定将失去物理依据。