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来抑制车身俯仰与侧倾,但也会影响轮胎动载荷。控制器设计时需在“车身姿态舒适”与“轮胎接地性”之间折中:阻尼过大,轮胎动载荷波动加剧,附着利用率下降;阻尼过小,车身姿态恶化,载荷转移更剧烈。
五、三向耦合与算法影响
三个方向并非独立:纵向滑移率影响侧向力(摩擦椭圆),侧向加速度引起载荷转移进而改变纵向力上限,垂向动载荷又反过来影响滑移率控制精度。因此底盘控制算法原型通常采用分层结构:上层做运动规划与参考模型,中层做纵向/侧向/垂向协调分配,下层做执行器控制。理解这些方程与关键变量的耦合关系,才能在原型阶段合理简化模型、设定控制目标与门限。
六、从方程到原型:一个可执行的建模检查清单
在算法原型阶段,建议按以下顺序把上述方程落地为可运行模型,避免一开始就陷入高保真仿真:
| 步骤 | 建模对象 | 关键输入 | 关键输出 | 验证方式 |
|---|---|---|---|---|
| 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,便于控制器参数整定。
综上,纵向、侧向、垂向方程分别给出了车速、航向、姿态的演化规律,而轮胎力、载荷转移、滑移率是连接方程与执行器的三个关键变量。原型阶段不必追求轮胎模型的高阶非线性,但必须保证滑移率-力曲线、侧偏角-力曲线、载荷转移方向与摩擦椭圆约束的正确性,否则后续控制器参数整定将失去物理依据。