基于Simulink的轮胎动力学模型（魔术公式）探索

基于simulink的轮胎动力学模型(魔术公式) 以路面附着系数、垂直载荷、车轮侧偏角为输入，纵向力 侧向力 和 附着椭圆为输出。

在车辆动力学的研究领域，轮胎动力学模型扮演着举足轻重的角色，其中魔术公式模型更是因其强大的描述能力而备受青睐。今天咱就唠唠基于Simulink搭建以路面附着系数、垂直载荷、车轮侧偏角为输入，纵向力、侧向力和附着椭圆为输出的轮胎动力学魔术公式模型。

魔术公式简介

魔术公式，正式名称是“Magic Formula”，由 Pacejka 教授提出。它能够以简洁的数学表达式精准地描述轮胎的力学特性。其基本形式如下：

\[y(x) = D \sin \{ C \arctan [B x - E (B x - \arctan (B x))] \}\]

这里的 \(x\) 是输入变量（如侧偏角、纵向滑移率等）， \(y(x)\) 是对应的输出力（纵向力或侧向力）。 \(B\)、\(C\)、\(D\)、\(E\) 是模型参数，它们决定了轮胎特性曲线的形状和位置，并且这些参数通常是垂直载荷 \(F_z\) 和路面附着系数 \(\mu\) 的函数。

Simulink 模型搭建

输入模块

咱先在Simulink里把输入搞定。创建三个输入端口，分别对应路面附着系数 \(\mu\)、垂直载荷 \(F_z\) 和车轮侧偏角 \(\alpha\) 。就像这样简单设置下：

% 假设用MATLAB脚本生成初始值作为测试 mu = 0.8; % 路面附着系数假设值 F_z = 5000; % 垂直载荷假设值，单位 N alpha = 5; % 车轮侧偏角假设值，单位 度

这几个值后续在Simulink里就作为输入参数的初始设定，当然实际应用中可能从其他模块实时获取。

参数计算模块

前面提到魔术公式的 \(B\)、\(C\)、\(D\)、\(E\) 参数和 \(\mu\)、\(Fz\) 有关。咱得根据经验公式或者试验数据拟合出这些关系。例如简单假设 \(D\) 与垂直载荷 \(Fz\) 成正比关系，\(D = kD * Fz\)，这里 \(kD\) 是比例系数，假设 \(kD = 0.1\)。

k_D = 0.1; D = k_D * F_z;

对于其他参数 \(B\)、\(C\)、\(E\) 类似地根据其与 \(\mu\)、\(F_z\) 的关系计算。这部分计算在Simulink里可以用Math Function模块来实现，把上述计算逻辑写成函数形式就行。

纵向力与侧向力计算模块

计算纵向力和侧向力就直接套用魔术公式。以侧向力 \(F_y\) 计算为例（假设侧偏角 \(\alpha\) 已经换算成弧度）：

% 假设已经计算好B、C、D、E参数 alpha_rad = deg2rad(alpha); F_y = D * sin(C * atan(B * alpha_rad - E * (B * alpha_rad - atan(B * alpha_rad))));

在Simulink里用一系列的数学运算模块搭建出这个公式的计算逻辑，输入就是前面算好的参数和侧偏角。纵向力 \(F_x\) 的计算类似，只不过输入变量可能是纵向滑移率（这里假设纵向滑移率计算好作为输入）。

附着椭圆计算模块

附着椭圆描述了轮胎纵向力和侧向力之间的相互关系，基于轮胎的附着特性。其方程为：

基于simulink的轮胎动力学模型(魔术公式) 以路面附着系数、垂直载荷、车轮侧偏角为输入，纵向力 侧向力 和 附着椭圆为输出。

\[\left( \frac{Fx}{F{xmax}} \right)^2 + \left( \frac{Fy}{F{ymax}} \right)^2 = 1\]

这里 \(F{xmax} = \mu Fz\)，\(F{ymax} = \mu Fz\)。在Simulink里，先计算出 \(F{xmax}\) 和 \(F{ymax}\) ，然后利用这个方程画出附着椭圆。这部分可能用到一些图形绘制相关的模块或者自定义函数，比如在MATLAB Function模块里写如下代码来绘制简单的附着椭圆数据点：

function [x_data, y_data] = draw_ellipse(mu, F_z, F_x, F_y) F_xmax = mu * F_z; F_ymax = mu * F_z; theta = 0:0.01:2*pi; x_data = F_xmax * cos(theta); y_data = F_ymax * sin(theta); % 这里简单把实际计算出的纵向力和侧向力也绘制进去 x_data = [x_data, F_x]; y_data = [y_data, F_y]; end

模型测试与结果分析

搭建好模型后，咱给不同的输入参数，比如改变路面附着系数 \(\mu\) 从 0.4 到 0.9 ，看看纵向力、侧向力和附着椭圆如何变化。要是 \(\mu\) 增大，从公式和实际物理意义都能知道， \(F{xmax}\) 和 \(F{ymax}\) 会增大，反映在附着椭圆上就是椭圆变大。纵向力和侧向力在相同侧偏角和滑移率下也会增大。

通过基于Simulink搭建这个轮胎动力学魔术公式模型，我们能很直观地研究轮胎在不同工况下的力学特性，为车辆动力学的进一步研究和整车性能优化提供有力支持。