1. 二自由度车辆模型基础概念
在车辆动力学研究中，二自由度模型是最基础也是最重要的分析工具之一。这个简化模型抓住了车辆横向运动的核心特性，特别适合分析转向工况下的动态响应。所谓"二自由度"，指的是模型主要考虑车辆的横向运动和横摆运动两个自由度。
质心侧偏角（β）是车辆纵向轴线与速度矢量之间的夹角，它反映了车辆实际运动方向与指向方向之间的偏差。当车辆转弯时，这个角度会明显显现。横摆角速度（γ）则是车辆绕垂直轴旋转的角速度，也就是我们常说的"转得快不快"。
这两个参数构成了车辆横向动力学分析的基础。在Simulink中建立这个模型时，我们需要基于以下基本方程：
code复制m(v̇ + uγ) = Fyf + Fyr
Izγ̇ = aFyf - bFyr
其中：
- m 为车辆质量
- u 为纵向车速（假设恒定）
- v 为横向速度
- Iz 为绕z轴的转动惯量
- a,b 分别为质心到前、后轴的距离
- Fyf,Fyr 为前、后轮胎侧向力
2. Simulink建模环境准备
2.1 基础模块配置
开始建模前，需要确保Simulink环境准备就绪。建议使用MATLAB R2020b或更新版本，这些版本对车辆动力学工具箱的支持更完善。新建一个空白模型后，我们需要准备以下基本模块：
- 输入模块：使用Constant或Signal Builder模块作为前轮转角输入
- 运算模块：主要使用Math Operations库中的加减乘除模块
- 积分模块：使用Continuous库中的Integrator模块
- 显示模块：使用Scope或To Workspace模块记录结果
提示：建模前建议先保存空白模型，Simulink有时会出现意外崩溃。可以设置自动保存功能：File > Preferences > General > AutoSave。
2.2 参数初始化设置
在模型初始化脚本中（通常为Model Properties > Callbacks > InitFcn），需要定义以下基础参数：
matlab复制m = 1500; % 车辆质量(kg)
Iz = 2500; % 横摆转动惯量(kg·m^2)
a = 1.2; % 前轴到质心距离(m)
b = 1.5; % 后轴到质心距离(m)
u = 20; % 纵向车速(m/s)
Cf = 80000; % 前轮侧偏刚度(N/rad)
Cr = 100000; % 后轮侧偏刚度(N/rad)
这些参数可以根据具体车型进行调整，是影响模型响应的关键因素。
3. 核心模型搭建过程
3.1 运动方程实现
根据二自由度模型方程，我们需要在Simulink中搭建对应的计算结构。以下是具体步骤：
- 创建前轮转角输入δ_f（使用Step或Ramp信号模拟转向输入）
- 计算前后轮侧向力：
code复制Fyf = Cf * (δ_f - (v+aγ)/u) Fyr = Cr * (- (v-bγ)/u) - 实现横向加速度计算：
code复制
v̇ = (Fyf + Fyr)/m - uγ - 实现横摆角加速度计算：
code复制γ̇ = (aFyf - bFyr)/Iz
在Simulink中，这些方程需要使用Sum、Product和Gain模块组合实现。特别注意符号的正负关系，这是模型正确性的关键。
3.2 积分环节设置
对v̇和γ̇进行积分得到v和γ：
- 添加两个Integrator模块，分别命名为'v_Integrator'和'gamma_Integrator'
- 设置初始条件为0（静止状态开始）
- 将输出分别连接到v和γ的信号线
质心侧偏角β通过v和u计算得到：
code复制β = atan(v/u)
在Simulink中使用Trigonometric Function模块实现arctan计算。
4. 相平面分析实现
4.1 相平面理论基础
相平面分析是研究动态系统稳定性的有力工具。对于车辆二自由度模型，我们通常绘制：
- 质心侧偏角β vs 横摆角速度γ
- 质心侧偏角β vs 其导数β̇
- 横摆角速度γ vs 其导数γ̇
这些相图可以直观展示车辆在不同工况下的稳定性特征。例如，相轨迹收敛表示系统稳定，发散则表示不稳定。
4.2 Simulink实现方法
在模型中添加以下模块实现相平面分析：
- 使用To Workspace模块记录β和γ的时间序列
- 在MATLAB脚本中处理数据：
matlab复制% 计算β的导数
beta_dot = gradient(beta, t);
% 绘制β-γ相平面
figure;
plot(beta, gamma);
xlabel('质心侧偏角β(rad)');
ylabel('横摆角速度γ(rad/s)');
grid on;
% 绘制β-β̇相平面
figure;
plot(beta, beta_dot);
xlabel('质心侧偏角β(rad)');
ylabel('质心侧偏角变化率β̇(rad/s)');
grid on;
4.3 典型工况测试
建议测试以下几种典型输入：
- 阶跃转向输入（Step信号）
- 正弦扫频转向输入（Sine Wave模块）
- 脉冲转向输入（Pulse Generator模块）
每种工况下的相图会呈现不同特征，反映了车辆在不同动态条件下的响应特性。
5. 参数影响分析与调试技巧
5.1 关键参数影响
通过修改以下参数，观察相图变化：
| 参数 | 影响效果 | 典型调整范围 |
|---|---|---|
| 车速u | 高速时相轨迹可能发散 | 10-40 m/s |
| 前轮侧偏刚度Cf | 影响不足转向特性 | 50000-120000 N/rad |
| 质量m | 影响响应速度 | 1000-2500 kg |
| 轴距分配(a/b) | 影响转向特性 | 0.7-1.5比例 |
5.2 常见问题排查
-
模型不收敛：
- 检查积分器设置，尝试减小步长
- 确认没有代数环问题（使用Unit Delay模块切断）
-
相图异常：
- 检查参数单位是否一致
- 确认arctan计算范围正确
-
数值不稳定：
- 尝试使用ode15s求解器
- 增加相对容差(RelTol)到1e-4
注意：当车速很高时，线性轮胎模型可能不再适用，此时相图会显示发散趋势，这实际反映了车辆真实的失稳现象。
6. 模型扩展与进阶应用
基础模型搭建完成后，可以考虑以下扩展方向：
-
非线性轮胎模型：
替换线性侧偏刚度，使用Pacejka魔术公式：code复制Fy = D*sin(C*atan(B*α - E*(B*α - atan(B*α)))) -
道路坡度影响：
在运动方程中加入重力分量：code复制m(v̇ + uγ) = Fyf + Fyr + mg*sinφ -
联合仿真：
将Simulink模型与CarSim或TruckSim等专业软件联合，获取更真实的车辆参数。 -
控制器设计：
基于相平面分析设计ESP或AEB控制器，观察相轨迹如何被控制到稳定区域。
在实际操作中，我发现一个很有用的技巧：将常用参数组保存为MAT文件，测试不同车型时可以直接加载。例如：
matlab复制% 保存参数组
sedan_params = struct('m',1500,'Iz',2500,...);
save('vehicle_params.mat','sedan_params');
% 加载参数组
load('vehicle_params.mat');
这样能大大提高工作效率，特别是在需要对比多种配置时。另外，建议为每个测试工况创建单独的模型副本，避免意外覆盖之前的成果。