双轨至单轨模型的渐近降阶理论

一、模型降阶的动机

1.1 高保真模型的局限性

双轨模型（7 自由度）虽然精确，但存在以下局限：

局限性	说明	后果
计算复杂度高	非线性方程 + 轮胎查表	实时仿真困难
控制器设计困难	MIMO 非线性系统	难以应用线性控制理论
参数辨识复杂	参数数量多	实验成本高
解析分析困难	耦合方程无解析解	物理洞察受限

1.2 降阶的目标

模型降阶（Model Order Reduction, MOR）的目标：

• 保留原系统的核心动态特性
• 降低状态变量数量
• 简化方程形式
• 便于控制器设计和分析

1.3 降阶的层次

模型	自由度数	适用场景
双轨完整模型	7+	高保真仿真、硬件在环
双轨简化模型	4	操纵稳定性分析
单轨动力学模型	2	控制器设计、稳定性分析
单轨运动学模型	2	低速路径规划

---

二、降阶的数学方法

2.1 奇异摄动理论

奇异摄动（Singular Perturbation）方法适用于多时间尺度系统。

系统标准形式：

$$\begin{array}{rl}\dot{\mathbf{x}}& =\mathbf{f}(\mathbf{x},\mathbf{z},\mathbf{u})\\ \epsilon \dot{\mathbf{z}}& =\mathbf{g}(\mathbf{x},\mathbf{z},\mathbf{u})\end{array}$$

其中：

• $\mathbf{x}$：慢状态
• $\mathbf{z}$：快状态
• $\epsilon \ll 1$：小参数

2.2 准稳态近似

当 $\epsilon \to 0$ 时：

$$\mathbf{g}(\mathbf{x},\mathbf{z},\mathbf{u})=0$$

解得快状态的准稳态解：

$${\mathbf{z}}_{\text{ss}}=\mathbf{h}(\mathbf{x},\mathbf{u})$$

代入慢状态方程：

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

2.3 车轮动力学的准稳态近似

车轮转动方程：

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

时间尺度分析：

动态	时间常数	相对速度
车轮转动	${\tau }_w\approx 0.05\mathrm{s}$	快
车身侧向	${\tau }_y\approx 0.3\mathrm{s}$	中
车身纵向	${\tau }_x\approx 1.0\mathrm{s}$	慢

车轮动力学的准稳态近似：

$$0=T_{\text{drive},i}-T_{\text{brake},i}-R_e{F}_{x,i}$$ $$F_{x,i}=\frac{T_{\text{drive},i}-T_{\text{brake},i}}{R_e}$$

---

三、双轨到单轨的渐近降阶

3.1 轮距效应的渐近分析

双轨模型横摆方程：

$$I_z\dot{r}=L_f{F}_{y,f}-L_r{F}_{y,r}+\frac{t_f}{2}\Delta F_{y,f}+\frac{t_r}{2}\Delta F_{y,r}$$

轮距项与主项的比值：

$${\epsilon }_t=\frac{t\Delta F_y}{2L{F}_y}$$

典型值分析：

工况	${\epsilon }_t$
小侧向加速度	$<0.05$
中等侧向加速度	$0.05\sim 0.15$
极限工况	$>0.15$

3.2 左右轮力的平均化

定义平均侧向力：

$${\bar{F}}_{y,f}=\frac{F_{y,fl}+F_{y,fr}}{2},{\bar{F}}_{y,r}=\frac{F_{y,rl}+F_{y,rr}}{2}$$

双轨模型简化为：

$$\begin{array}{rl}m({\dot{v}}_y+r{v}_x)& =2{\bar{F}}_{y,f}+2{\bar{F}}_{y,r}\\ I_z\dot{r}& =2L_f{\bar{F}}_{y,f}-2L_r{\bar{F}}_{y,r}\end{array}$$

3.3 等效单轮参数

等效前轮侧偏刚度：

$$C_{\alpha f}^{\text{eq}}=2{\bar{C}}_{\alpha f}$$

等效后轮侧偏刚度：

$$C_{\alpha r}^{\text{eq}}=2{\bar{C}}_{\alpha r}$$

单轨模型方程恢复：

$$\begin{array}{rl}m({\dot{v}}_y+r{v}_x)& =-C_{\alpha f}^{\text{eq}}{\alpha }_f-C_{\alpha r}^{\text{eq}}{\alpha }_r\\ I_z\dot{r}& =-L_f{C}_{\alpha f}^{\text{eq}}{\alpha }_f+L_r{C}_{\alpha r}^{\text{eq}}{\alpha }_r\end{array}$$

---

四、降阶误差分析

4.1 轮距忽略误差

轮距项引起的横摆力矩误差：

$$\Delta M_z=\frac{t_f}{2}(F_{y,fl}-F_{y,fr})+\frac{t_r}{2}(F_{y,rl}-F_{y,rr})$$

相对误差：

$${\epsilon }_M=\frac{|\Delta M_z|}{|L_f{F}_{y,f}-L_r{F}_{y,r}|}$$

4.2 载荷转移忽略误差

双轨模型考虑载荷转移：

$$C_{\alpha }^{\text{eff}}=C_{\alpha }(F_z^{\text{avg}})\cdot (1-{\epsilon }_{\text{load}})$$

载荷转移因子：

$${\epsilon }_{\text{load}}=\frac{n(1-n)}{2}{\left(\frac{\Delta F_z}{F_{z0}}\right)}^2$$

典型值：${\epsilon }_{\text{load}}\approx 0.05\sim 0.15$（$0.5g$ 转向）

4.3 降阶误差总结

误差源	符号	典型范围	建议处理
轮距忽略	${\epsilon }_t$	$0.02\sim 0.1$	可忽略（线性区）
载荷转移	${\epsilon }_{\text{load}}$	$0.05\sim 0.2$	等效刚度修正
左右力差	${\epsilon }_F$	$0.01\sim 0.05$	可忽略
耦合项	${\epsilon }_c$	$0.05\sim 0.2$	增益调度

---

五、改进的单轨模型

5.1 等效参数修正

载荷转移修正：

$$C_{\alpha }^{\text{eff}}=C_{\alpha 0}\left[1-k_{\text{load}}{\left(\frac{a_y}{g}\right)}^2\right]$$

典型值：$k_{\text{load}}\approx 0.1\sim 0.2$

轮距修正：

$$K^{\text{eff}}=K_0+\Delta K_t$$ $$\Delta K_t=\frac{m}{L}\left(\frac{t_f^2}{8C_{\alpha f}}+\frac{t_r^2}{8C_{\alpha r}}\right)$$

5.2 增益调度单轨模型

车速调度的参数：

$$\begin{array}{rl}{C}_{\alpha f}(v_x)& =C_{\alpha f0}(1-k_v{v}_x^2)\\ C_{\alpha r}(v_x)& =C_{\alpha r0}(1-k_v{v}_x^2)\end{array}$$

侧向加速度调度：

$$C_{\alpha }(a_y)=C_{\alpha 0}\left[1-k_a{\left(\frac{a_y}{g}\right)}^2\right]$$

5.3 非线性单轨模型

保留轮胎非线性但忽略左右差异：

$$\begin{array}{rl}{F}_{y,f}& =F_{y,f}({\alpha }_f)\text{(魔术公式)}\\ F_{y,r}& =F_{y,r}({\alpha }_r)\text{(魔术公式)}\end{array}$$

方程形式：

$$\begin{array}{rl}m({\dot{v}}_y+r{v}_x)& =F_{y,f}({\alpha }_f)+F_{y,r}({\alpha }_r)\\ I_z\dot{r}& =L_f{F}_{y,f}({\alpha }_f)-L_r{F}_{y,r}({\alpha }_r)\end{array}$$

---

六、数值验证

6.1 仿真工况

工况	描述	验证重点
阶跃转向	$\delta =5^{\circ }$, $v_x=80km/h$	瞬态响应
正弦转向	$f=0.5\mathrm{Hz}$, ${\delta }_0=3^{\circ }$	频响特性
极限转向	$\delta =10^{\circ }$, $v_x=60km/h$	非线性精度

6.2 响应比较

横摆角速度响应（阶跃转向）：

模型	$r_{\text{peak}}(rad/s)$	$t_p(\mathrm{s})$	稳态误差
双轨完整	0.52	0.71	-
双轨简化	0.51	0.72	$-2\%$
单轨线性	0.48	0.75	$-8\%$
单轨非线性	0.50	0.73	$-4\%$

6.3 相轨迹比较

graph TD
    subgraph β-r 相平面
    A[双轨模型轨迹] --> B[收敛到平衡点]
    C[单轨模型轨迹] --> D[收敛到平衡点]
    B --> E{稳态误差}
    D --> E
    E --> |< 5%| F[可接受]
    E --> |> 5%| G[需修正]
    end

---

七、降阶模型的应用

7.1 控制器设计

单轨模型适用于：

• LQR 横摆稳定性控制
• MPC 轨迹跟踪控制
• 滑模控制设计

7.2 稳定性分析

单轨模型可用于：

• 特征值分析
• 根轨迹分析
• 相平面稳定性分析

7.3 实时仿真

计算效率对比：

模型	计算时间/步长	实时性
双轨完整	$500\mu \mathrm{s}$	×（需 $<100\mu \mathrm{s}$）
单轨非线性	$100\mu \mathrm{s}$	✓
单轨线性	$20\mu \mathrm{s}$	✓✓

---

八、相关内容

• 02-平面动力学向高阶自由度扩展的理论框架
• 01-极限侧滑工况下的相平面稳定性理论
• 02-纵横向联合激励下的动力学耦合机理