纵横向联合激励下的动力学耦合机理

一、联合激励的物理机理

1.1 纵横向耦合的来源

车辆平面动力学中，纵向和侧向运动通过以下路径耦合：

耦合路径	机理	数学表达
运动学耦合	车身坐标系旋转	$a_x={\dot{v}}_x-r{v}_y$, $a_y={\dot{v}}_y+r{v}_x$
轮胎力耦合	摩擦椭圆约束	$F_x^2+F_y^2\le (\mu F_z{)}^2$
载荷转移耦合	联合加速度效应	$F_z=f(a_x,a_y)$
刚度耦合	载荷敏感性	$C_{\alpha }=C_{\alpha }(F_z)$

1.2 耦合的强度度量

定义耦合系数：

$${\eta }_{xy}=\frac{\text{耦合项贡献}}{\text{总响应}}=\frac{|r{v}_y|}{|{\dot{v}}_x|}+\frac{|r{v}_x|}{|{\dot{v}}_y|}$$

${\eta }_{xy}$ 范围	耦合强度	模型选择
$<0.1$	弱耦合	可解耦分析
$0.1\sim 0.3$	中等耦合	需考虑耦合项
$>0.3$	强耦合	必须联合分析

---

二、运动学耦合分析

2.1 科里奥利耦合项

车身坐标系中的加速度：

$$\begin{array}{rl}{a}_x& ={\dot{v}}_x-r{v}_y\\ a_y& ={\dot{v}}_y+r{v}_x\end{array}$$

耦合项：

• $-r{v}_y$：纵向加速度中的科里奥利项
• $+r{v}_x$：侧向加速度中的科里奥利项

2.2 科里奥利项的量级分析

典型工况下：

| 工况 | $v_x$ | $v_y$ | $r$ | $|r{v}_y|$ | $|r{v}_x|$ | |------|-------|-------|-----|-----------|-----------| | 高速直线 | $30m/s$ | $0.1m/s$ | $0.01rad/s$ | $0.001m/{\mathrm{s}}^2$ | $0.3m/{\mathrm{s}}^2$ | | 稳态转向 | $20m/s$ | $2m/s$ | $0.4rad/s$ | $0.8m/{\mathrm{s}}^2$ | $8m/{\mathrm{s}}^2$ | | 极限漂移 | $15m/s$ | $8m/s$ | $0.6rad/s$ | $4.8m/{\mathrm{s}}^2$ | $9m/{\mathrm{s}}^2$ |

物理意义

高速转向时，科里奥利项 $r{v}_x$ 对侧向加速度的贡献可能超过 ${\dot{v}}_y$。

2.3 耦合对稳定性的影响

线性化方程：

$$\begin{array}{rl}\delta a_x& =\delta {\dot{v}}_x-v_{y0}\delta r-r_0\delta v_y\\ \delta a_y& =\delta {\dot{v}}_y+v_{x0}\delta r+r_0\delta v_y\end{array}$$

耦合矩阵：

$$\mathbf{C}=\left[\begin{array}{ccc}0& -r_0& -v_{y0}\\ 0& 0& v_{x0}+r_0\\ 0& 0& 0\end{array}\right]$$

---

三、轮胎力耦合分析

3.1 摩擦椭圆约束

联合工况下，轮胎力满足：

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

等效摩擦系数：

$${\mu }_{\text{eff}}=\sqrt{{\mu }^2-{\left(\frac{F_x}{F_z}\right)}^2}$$

侧向力的有效上限：

$$|F_y|\le {\mu }_{\text{eff}}{F}_z=\sqrt{(\mu F_z{)}^2-F_x^2}$$

3.2 耦合对侧偏刚度的影响

有效侧偏刚度：

$$C_{\alpha }^{\text{eff}}={\frac{\partial F_y}{\partial \alpha }|}_{\kappa \ne 0}=C_{\alpha }\cdot \frac{{\mu }_{\text{eff}}}{\mu }$$

纵向力存在时，侧偏刚度降低：

$$\frac{C_{\alpha }^{\text{eff}}}{C_{\alpha }}=\sqrt{1-{\left(\frac{F_x}{\mu F_z}\right)}^2}$$

3.3 联合滑移的 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}}$$

---

四、载荷转移耦合分析

4.1 联合加速度下的载荷转移

四轮垂直载荷：

$$\begin{array}{rl}{F}_{z,fl}& =\frac{L_r}{2L}mg+\frac{h_{\text{CG}}}{2L}m{a}_x+{\eta }_f\frac{m{a}_y{h}_{\text{CG}}}{t_f}\\ F_{z,fr}& =\frac{L_r}{2L}mg+\frac{h_{\text{CG}}}{2L}m{a}_x-{\eta }_f\frac{m{a}_y{h}_{\text{CG}}}{t_f}\\ F_{z,rl}& =\frac{L_f}{2L}mg-\frac{h_{\text{CG}}}{2L}m{a}_x+{\eta }_r\frac{m{a}_y{h}_{\text{CG}}}{t_r}\\ F_{z,rr}& =\frac{L_f}{2L}mg-\frac{h_{\text{CG}}}{2L}m{a}_x-{\eta }_r\frac{m{a}_y{h}_{\text{CG}}}{t_r}\end{array}$$

4.2 外前轮载荷最大化

极限工况（制动 + 转向）：

$$F_{z,fl}^{\text{max}}=\frac{L_r}{2L}mg+\frac{h_{\text{CG}}}{2L}m|a_x{|}_{\text{max}}+{\eta }_f\frac{m|a_y{|}_{\text{max}}{h}_{\text{CG}}}{t_f}$$

典型值：$F_{z,fl}^{\text{max}}\approx 2\sim 3\times F_{z0}$

4.3 内后轮载荷损失

内后轮可能离地：

$$F_{z,rr}\le 0\Rightarrow \frac{L_f}{2L}mg\le \frac{h_{\text{CG}}}{2L}m|a_x|+{\eta }_r\frac{m|a_y|h_{\text{CG}}}{t_r}$$

离地临界条件：

$$|a_x|+\frac{2L{\eta }_r}{t_r}|a_y|\ge \frac{L_{f}g}{h_{\text{CG}}}$$

---

五、耦合系统的稳定性分析

5.1 线性化状态矩阵

考虑纵横向耦合的状态矩阵：

$$\mathbf{A}=\left[\begin{array}{ccc}{A}_{xx}& A_{xy}& A_{xr}\\ A_{yx}& A_{yy}& A_{yr}\\ A_{rx}& A_{ry}& A_{rr}\end{array}\right]$$

耦合元素：

• $A_{xy}=-r_0$（纵向 - 侧向耦合）
• $A_{yx}=r_0$（侧向 - 纵向耦合）
• $A_{xr},A_{yr}$（横摆耦合）

5.2 特征值轨迹

随耦合强度变化，特征值轨迹：

graph TD
    subgraph 复平面
    A[原点] --> B{耦合强度}
    B -->|弱 | C[左半平面稳定]
    B -->|中 | D[向虚轴靠近]
    B -->|强 | E[穿越虚轴不稳定]
    end

5.3 耦合导致的失稳机理

失稳模式 1：横摆发散

• 特征值穿越虚轴
• 横摆角速度指数增长

失稳模式 2：漂移振荡

• 特征值穿越虚轴为共轭复根
• 极限环振荡

---

六、数值示例

6.1 车辆参数

参数	符号	数值
整车质量	$m$	$1500\mathrm{kg}$
轴距	$L$	$2.7\mathrm{m}$
质心高度	$h_{\text{CG}}$	$0.55\mathrm{m}$
前轮距	$t_f$	$1.6\mathrm{m}$
后轮距	$t_r$	$1.55\mathrm{m}$
摩擦系数	$\mu$	$0.9$

6.2 纯转向工况（$a_x=0$）

$a_y$	$C_{\alpha }^{\text{eff}}/C_{\alpha }$	稳定性
$0.3g$	$0.95$	稳定
$0.6g$	$0.82$	稳定
$0.8g$	$0.65$	临界

6.3 联合工况（$a_x=-0.3g,a_y=0.6g$）

车轮	$F_z(\mathrm{N})$	$C_{\alpha }^{\text{eff}}/C_{\alpha 0}$
前外 (fl)	$7200$	$0.92$
前内 (fr)	$3800$	$0.78$
后外 (rl)	$4500$	$0.85$
后内 (rr)	$1900$	$0.65$

有效轴侧偏刚度：

$$\begin{array}{rl}{C}_{\alpha f}^{\text{eff}}& =0.85C_{\alpha f0}\\ C_{\alpha r}^{\text{eff}}& =0.75C_{\alpha r0}\end{array}$$

不足转向梯度变化：

$$K^{\text{eff}}=1.18K_0$$

联合工况使不足转向梯度增加 18%。

---

七、控制应用

7.1 解耦控制策略

输入 - 输出解耦：

设计控制律使：

$$\mathbf{u}={\mathbf{K}}_x\mathbf{x}+{\mathbf{K}}_d\mathbf{d}$$

闭环系统：

$$\dot{\mathbf{x}}=(\mathbf{A}-\mathbf{B}{\mathbf{K}}_x)\mathbf{x}$$

目标：使闭环矩阵对角化。

7.2 ESP 的耦合补偿

ESP 制动力分配考虑耦合：

$$T_{\text{brake},i}=T_{\text{brake},i}^0+\Delta T_{\text{coupling}}$$ $$\Delta T_{\text{coupling}}=K_c(a_x{a}_y)$$

7.3 扭矩矢量分配的解耦

独立驱动可主动利用耦合：

$$\begin{array}{rl}\Delta T_{\text{yaw}}& =\frac{t}{2}(F_{x,fr}-F_{x,fl})\\ \Delta F_x& =F_{x,fl}+F_{x,fr}\end{array}$$

---

八、相关内容

• 01-极限侧滑工况下的相平面稳定性理论
• 02-完整状态向量定义与多输入多输出方程构建
• 01-轮胎非线性特性与摩擦椭圆约束的全面引入