二自由度车辆模型相平面分析与Matlab仿真

1. 项目背景与核心概念解析

在车辆动力学研究中，二自由度模型是最基础也最经典的车辆稳定性分析工具。这个看似简单的模型，却能揭示车辆在极限工况下的失稳机理。相平面分析法作为非线性系统研究的利器，通过将车辆状态变量（质心侧偏角和横摆角速度）的关系可视化，为工程师提供了直观的稳定性判据。

相平面上的每个点都代表车辆的瞬态运动状态，而临界轨迹则是区分稳定域和不稳定域的分界线。鞍点作为相平面中的特殊平衡点，对系统稳定性具有决定性影响。通过Matlab仿真绘制这些关键要素，我们能够：

• 量化评估车辆在不同工况下的稳定裕度
• 预测车辆失稳的临界条件
• 验证控制策略的有效性

提示：二自由度模型虽然简化了轮胎力、悬架等复杂因素，但抓住了车辆横向动力学最本质的特征，是理解ESC（电子稳定控制系统）工作原理的基础。

2. 二自由度车辆模型构建

2.1 基本假设与坐标系定义

建立模型前需要明确几个关键假设：

1. 忽略悬架运动，车辆简化为平行于地面的刚体
2. 仅考虑横向和横摆两个自由度
3. 小角度假设（侧偏角<5°时线性轮胎模型有效）
4. 匀速行驶工况（纵向速度恒定）

采用车辆坐标系：

• 原点位于质心
• x轴指向车辆前进方向
• y轴指向驾驶员左侧
• z轴垂直向上（右手定则）

2.2 动力学方程推导

通过牛顿-欧拉法建立运动方程：

横向动力学方程：
m(v̇ + ur) = F_yf + F_yr
其中：

• m：整车质量
• v：质心侧向速度
• u：纵向速度（恒定）
• r：横摆角速度
• F_yf, F_yr：前/后轴侧向力

横摆动力学方程：
I_z ṙ = aF_yf - bF_yr
其中：

• I_z：绕z轴的转动惯量
• a,b：质心到前/后轴距离

线性轮胎模型：
F_yf = -C_f α_f
F_yr = -C_r α_r
其中：

• C_f, C_r：前/后轮侧偏刚度
• α_f, α_r：前/后轮侧偏角

运动学关系：
α_f ≈ δ - (v + ar)/u
α_r ≈ -(v - br)/u
（δ为前轮转角）

2.3 状态空间表达

选择状态变量x = [v; r]，输入u = δ，整理得到标准状态方程：

ẋ = Ax + Bu