轮胎非线性特性与摩擦椭圆约束的全面引入

一、轮胎非线性的来源

1.1 线性模型的局限性

线性轮胎模型：

$$F_y=-C_{\alpha }\alpha$$

适用范围：

• 侧偏角 $\alpha <4^{\circ }$
• 侧向加速度 $a_y<0.3g$
• 小滑移率 $|\kappa |<0.05$

局限性：

• 无法描述轮胎饱和特性
• 无法处理联合滑移工况
• 无法预测极限操纵行为

1.2 非线性的物理机理

非线性源	机理	数学表现
橡胶材料非线性	粘弹性滞后	$F(\alpha )$ 饱和曲线
接触patch 变形	剪切应力分布	渐进滑移
摩擦极限	库伦摩擦	$
载荷敏感性	接触面积变化	$C_{\alpha }(F_z)$ 非线性

1.3 典型的侧向力 - 侧偏角曲线

graph LR
    A[侧偏角 α] --> B{α 大小判断}
    B -->|α < 3°| C[线性区 F_y ≈ -C_α α]
    B -->|3° < α < 8°| D[过渡区 二阶项显著]
    B -->|α > 8°| E[饱和区 F_y → μF_z]

---

二、魔术公式轮胎模型

2.1 魔术公式的标准形式

Pacejka 魔术公式（Magic Formula）：

$$F(\alpha )=D\sin \left[C\arctan \left(B\alpha -E(B\alpha -\arctan (B\alpha ))\right)\right]$$

其中：

• $B$：刚度因子（Stiffness factor）
• $C$：形状因子（Shape factor）
• $D$：峰值因子（Peak factor）
• $E$：曲率因子（Curvature factor）

2.2 参数的物理意义

参数	符号	物理意义	典型值
峰值因子	$D$	最大侧向力 $\mu F_z$	$0.8\sim 1.2F_z$
形状因子	$C$	曲线整体形状	$1.3\sim 1.7$
刚度因子	$B$	原点斜率相关	$C_{\alpha }/(CD)$
曲率因子	$E$	过渡区曲率	$-2\sim 2$

2.3 侧偏刚度的关系

原点刚度：

$$C_{\alpha }={\frac{d{F}_y}{d\alpha }|}_{\alpha =0}=BCD$$

因此：

$$B=\frac{C_{\alpha }}{CD}$$

2.4 载荷依赖性

魔术公式参数随垂直载荷变化：

$$\begin{array}{rl}D& =\mu F_z(a_1+a_2\Delta F_z)\\ C& =b_1+b_2\Delta F_z\\ E& =c_1+c_2\Delta F_z\\ BCD& =C_{\alpha 0}(1-d_1|\Delta F_z|)\end{array}$$

其中 $\Delta F_z=(F_z-F_{z0})/F_{z0}$。

---

三、摩擦椭圆约束

3.1 联合滑移的物理机理

轮胎同时承受纵向滑移率 $\kappa$ 和侧偏角 $\alpha$ 时，接触 patch 内的摩擦应力同时产生纵向力 $F_x$ 和侧向力 $F_y$。

摩擦极限条件：

$$\sqrt{F_x^2+F_y^2}\le \mu F_z$$

3.2 摩擦椭圆方程

摩擦椭圆（Friction Ellipse）描述了联合滑移下的力约束：

$${\left(\frac{F_x}{{\mu }_x{F}_z}\right)}^2+{\left(\frac{F_y}{{\mu }_y{F}_z}\right)}^2\le 1$$

或等效形式：

$${\left(\frac{F_x}{F_{x,\text{max}}}\right)}^2+{\left(\frac{F_y}{F_{y,\text{max}}}\right)}^2\le 1$$

其中：

• $F_{x,\text{max}}={\mu }_x{F}_z$：纯纵向滑移最大力
• $F_{y,\text{max}}={\mu }_y{F}_z$：纯侧偏最大力

graph TD
    subgraph 摩擦椭圆
    A[F_x 轴] --> B[椭圆边界]
    C[F_y 轴] --> B
    D[工作点 F_x, F_y] --> E{是否在椭圆内？}
    E -->|是 | F[轮胎未饱和]
    E -->|否 | G[超出摩擦极限]
    end

3.3 联合滑移下的力计算

简化模型（弹性摩擦圆）：

$$\begin{array}{rl}{F}_x& =F_{x0}(\kappa )\cdot \frac{{\mu }_{\text{eff}}}{\mu }\\ F_y& =F_{y0}(\alpha )\cdot \frac{{\mu }_{\text{eff}}}{\mu }\end{array}$$

其中有效摩擦系数：

$${\mu }_{\text{eff}}=\mu \sqrt{1-{\left(\frac{F_{x0}}{\mu F_z}\right)}^2-{\left(\frac{F_{y0}}{\mu F_z}\right)}^2}$$

3.4 Dugoff 轮胎模型

Dugoff 模型提供了摩擦椭圆的解析形式：

$$\begin{array}{rl}{F}_x& =\frac{C_{\kappa }\kappa }{1+\kappa }f(\lambda )\\ F_y& =\frac{C_{\alpha }\tan \alpha }{1+\kappa }f(\lambda )\end{array}$$

其中：

$$\lambda =\frac{\mu F_z(1+\kappa )}{2\sqrt{(C_{\kappa }\kappa {)}^2+(C_{\alpha }\tan \alpha {)}^2}}$$ $$f(\lambda )=\{\begin{array}{ll}\lambda (2-\lambda ),& \lambda <1\\ 1,& \lambda \ge 1\end{array}$$

物理意义

$f(\lambda )$ 是摩擦限制因子：

• $\lambda \ge 1$：轮胎未饱和，线性区
• $\lambda <1$：轮胎部分或完全饱和

---

四、双轨模型的非线性动力学方程

4.1 状态向量的扩展

双轨非线性模型的状态向量：

$$\mathbf{x}={\left[\begin{array}{ccccccc}{v}_x& v_y& r& {\omega }_{fl}& {\omega }_{fr}& {\omega }_{rl}& {\omega }_{rr}\end{array}\right]}^{\top }$$

共 7 个状态变量。

4.2 纵向动力学方程

$$\begin{array}{rl}m({\dot{v}}_x-r{v}_y)=& \sum _{i=fl,fr}(F_{x,i}\cos {\delta }_i-F_{y,i}\sin {\delta }_i)\\ & +\sum _{i=rl,rr}{F}_{x,i}-F_{\text{drag}}\end{array}$$

其中空气阻力：

$$F_{\text{drag}}=\frac{1}{2}\rho C_{d}A{v}_x^2$$

4.3 侧向动力学方程

$$\begin{array}{rl}m({\dot{v}}_y+r{v}_x)=& \sum _{i=fl,fr}(F_{x,i}\sin {\delta }_i+F_{y,i}\cos {\delta }_i)\\ & +\sum _{i=rl,rr}{F}_{y,i}\end{array}$$

4.4 横摆动力学方程

$$\begin{array}{rl}{I}_z\dot{r}=& \sum _{i=fl,fr}\left[L_f(F_{x,i}\sin {\delta }_i+F_{y,i}\cos {\delta }_i)\pm \frac{t_f}{2}(F_{x,i}\cos {\delta }_i-F_{y,i}\sin {\delta }_i)\right]\\ & +\sum _{i=rl,rr}\left[-L_r{F}_{y,i}\pm \frac{t_r}{2}{F}_{x,i}\right]\end{array}$$

4.5 车轮转动动力学

每个车轮的转动方程：

$$I_w{\dot{\omega }}_i=T_{\text{drive},i}-T_{\text{brake},i}-R_e{F}_{x,i}$$

其中：

• $I_w$：车轮转动惯量
• $T_{\text{drive},i}$：驱动力矩
• $T_{\text{brake},i}$：制动力矩
• $R_e$：轮胎有效半径

---

五、非线性方程的数值求解

5.1 状态空间的标准形式

非线性状态空间方程：

$$\dot{\mathbf{x}}=\mathbf{f}(\mathbf{x},\mathbf{u})$$

其中输入向量：

$$\mathbf{u}={\left[\begin{array}{ccc}\delta & T_{\text{drive}}& T_{\text{brake},i}\end{array}\right]}^{\top }$$

5.2 数值积分方法

四阶 Runge-Kutta 法（RK4）：

$$\begin{array}{rl}{k}_1& =\mathbf{f}({\mathbf{x}}_n,{\mathbf{u}}_n)\\ k_2& =\mathbf{f}({\mathbf{x}}_n+\frac{h}{2}{k}_1,{\mathbf{u}}_{n+1/2})\\ k_3& =\mathbf{f}({\mathbf{x}}_n+\frac{h}{2}{k}_2,{\mathbf{u}}_{n+1/2})\\ k_4& =\mathbf{f}({\mathbf{x}}_n+h{k}_3,{\mathbf{u}}_{n+1})\\ {\mathbf{x}}_{n+1}& ={\mathbf{x}}_n+\frac{h}{6}(k_1+2k_2+2k_3+k_4)\end{array}$$

时间步长选择：

工况	推荐步长	理由
稳态转向	$10\mathrm{ms}$	低频动态
阶跃输入	$5\mathrm{ms}$	中频响应
ABS/ESP 干预	$1\mathrm{ms}$	高频切换

5.3 雅可比矩阵的解析计算

对于隐式积分或稳定性分析，需要雅可比矩阵：

$$\mathbf{J}=\frac{\partial \mathbf{f}}{\partial \mathbf{x}}$$

元素形式：

$$J_{ij}=\frac{\partial f_i}{\partial x_j}$$

由于轮胎力的非线性，雅可比矩阵需要数值计算或查表插值。

---

六、极限工况下的动力学特性

6.1 轮胎饱和的影响

当轮胎进入饱和区（$\alpha >8^{\circ }$ 或 $|\kappa |>0.1$）：

效应	线性模型预测	非线性模型预测
侧向力	随 $\alpha$ 线性增长	趋于 $\mu F_z$ 饱和
横摆阻尼	恒定	显著降低
稳定性边界	固定	随工况变化
极限环	不存在	可能出现

6.2 不足转向的极限变化

线性模型预测的不足转向梯度 $K$ 为常数。

非线性模型中，有效 $K$ 随侧向加速度变化：

$$K_{\text{eff}}(a_y)=\frac{m}{L}\left(\frac{L_r}{C_{\alpha r}^{\text{eff}}(a_y)}-\frac{L_f}{C_{\alpha f}^{\text{eff}}(a_y)}\right)$$

通常 $K_{\text{eff}}$ 随 $a_y$ 增大而增大（更趋向不足转向）。

6.3 横摆角速度的极限

最大稳态横摆角速度受摩擦极限限制：

$$|r_{\text{max}}|=\frac{\mu g}{v_x}$$

超过此值，车辆将进入侧滑状态。

---

七、数值示例

7.1 车辆与轮胎参数

参数	符号	数值
整车质量	$m$	$1500\mathrm{kg}$
前轮侧偏刚度	$C_{\alpha f}$	$60000N/rad$
后轮侧偏刚度	$C_{\alpha r}$	$55000N/rad$
摩擦系数	$\mu$	$0.9$（干沥青）
峰值因子	$D$	$0.9F_z$
形状因子	$C$	$1.5$
曲率因子	$E$	$-0.5$

7.2 线性与非线性响应对比

在 $\delta =5^{\circ }$、$v_x=80km/h$ 的稳态转向工况下：

响应	线性模型	非线性模型	差异
稳态横摆角速度	$0.46rad/s$	$0.43rad/s$	$-6.5\%$
侧偏角 ${\alpha }_f$	$3.2^{\circ }$	$3.4^{\circ }$	$+6\%$
侧向力 $F_{y,f}$	$3420\mathrm{N}$	$3280\mathrm{N}$	$-4\%$

7.3 极限侧向加速度

模型	$a_{y,\text{max}}$
线性（外推）	$>1.0g$（无意义）
非线性	$\approx 0.85g$

---

八、相关内容

• 02-完整状态向量定义与多输入多输出方程构建
• 02-准静态平面载荷转移理论推导
• 01-极限侧滑工况下的相平面稳定性理论