平面动力学向高阶自由度扩展的理论框架

一、高阶扩展的动机

1.1 平面模型的假设与局限

平面动力学模型（2DOF/4DOF）基于以下假设：

假设	内容	局限
平面运动假设	忽略垂向运动	无法分析路面不平激励
无侧倾/俯仰	$\varphi =\theta =0$	忽略载荷转移动态
左右对称	左右轮特性相同	无法模拟路面μ分裂
恒定车速	$v_x=\text{const}$	忽略纵向动力学耦合

1.2 需要高阶扩展的场景

场景	需要增加的 DOF	物理由来
平顺性分析	垂向 $z$、侧倾 $\varphi$、俯仰 $\theta$	悬架弹性
极限操纵	侧倾 $\varphi$、侧倾角速度 $p$	载荷转移动态
主动悬架控制	四轮悬架行程 $z_i$	作动器动态
制动稳定性	纵向 $v_x$、俯仰 $\theta$	载荷前移
翻滚稳定性	侧倾 $\varphi$、侧倾角速度 $p$	高质心车辆

---

二、7DOF 完整车辆模型

2.1 自由度的定义

完整 7DOF 模型包含：

$$\mathbf{q}={\left[\begin{array}{ccccccc}X& Y& Z& \varphi & \theta & \psi & \dots \end{array}\right]}^{\top }$$

或车身坐标系：

$$\mathbf{x}={\left[\begin{array}{cccccc}{v}_x& v_y& v_z& p& q& r\end{array}\right]}^{\top }$$

自由度	符号	含义	单位
纵向平动	$v_x$	质心纵向速度	$m/s$
侧向平动	$v_y$	质心侧向速度	$m/s$
垂向平动	$v_z$	质心垂向速度	$m/s$
侧倾转动	$p$	绕 $x$ 轴角速度	$rad/s$
俯仰转动	$q$	绕 $y$ 轴角速度	$rad/s$
横摆转动	$r$	绕 $z$ 轴角速度	$rad/s$
（车轮转动）	${\omega }_i$	四轮角速度	$rad/s$

2.2 牛顿 - 欧拉方程的完整形式

平动方程：

$$m\left[\begin{array}{c}{\dot{v}}_x-r{v}_y+q{v}_z\\ {\dot{v}}_y-p{v}_z+r{v}_x\\ {\dot{v}}_z-q{v}_x+p{v}_y\end{array}\right]=\left[\begin{array}{c}{F}_x\\ F_y\\ F_z\end{array}\right]$$

转动方程：

$$\left[\begin{array}{ccc}{I}_x& 0& 0\\ 0& I_y& 0\\ 0& 0& I_z\end{array}\right]\left[\begin{array}{c}\dot{p}\\ \dot{q}\\ \dot{r}\end{array}\right]+\left[\begin{array}{c}p\\ q\\ r\end{array}\right]\times \mathbf{I}\left[\begin{array}{c}p\\ q\\ r\end{array}\right]=\left[\begin{array}{c}{M}_x\\ M_y\\ M_z\end{array}\right]$$

2.3 耦合项的分析

科里奥利耦合：

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

耦合项	来源	量级
$-r{v}_y$	横摆 - 侧向	中等
$q{v}_z$	俯仰 - 垂向	小
$-p{v}_z$	侧倾 - 垂向	小
$r{v}_x$	横摆 - 纵向	大
$-q{v}_x$	俯仰 - 纵向	中等
$p{v}_y$	侧倾 - 侧向	小

---

三、侧倾动力学的引入

3.1 侧倾运动的物理机理

侧倾角 $\varphi$：车身绕 $x$ 轴的转动。

侧倾力矩平衡：

$$I_x\dot{p}+c_{\varphi }p+k_{\varphi }\varphi =m{a}_y{h}_{\text{CG}}+mg{h}_{\text{CG}}\sin \varphi$$

其中：

• $I_x$：侧倾转动惯量
• $c_{\varphi }$：侧倾阻尼系数
• $k_{\varphi }$：侧倾刚度
• $h_{\text{CG}}$：质心到侧倾轴的距离

3.2 侧倾动力学的状态方程

状态变量：$\varphi$（侧倾角）、$p$（侧倾角速度）

$$\begin{array}{rl}\dot{\varphi }& =p\\ \dot{p}& =\frac{1}{I_x}\left(m{a}_y{h}_{\text{CG}}+mg{h}_{\text{CG}}\sin \varphi -c_{\varphi }p-k_{\varphi }\varphi \right)\end{array}$$

小角度近似：

$$\begin{array}{rl}\dot{\varphi }& =p\\ \dot{p}& =\frac{m{h}_{\text{CG}}}{I_x}{a}_y-\frac{c_{\varphi }}{I_x}p-\frac{k_{\varphi }-mg{h}_{\text{CG}}}{I_x}\varphi \end{array}$$

3.3 侧倾对载荷转移的动态影响

动态载荷转移：

$$\Delta F_{z,f}(t)={\eta }_f\frac{1}{t_f}\left(m{a}_y(t)h_{\text{CG}}+k_{\varphi ,f}\varphi (t)+c_{\varphi ,f}p(t)\right)$$ $$\Delta F_{z,r}(t)={\eta }_r\frac{1}{t_r}\left(m{a}_y(t)h_{\text{CG}}+k_{\varphi ,r}\varphi (t)+c_{\varphi ,r}p(t)\right)$$

3.4 侧倾动力学的固有频率

特征方程：

$$s^2+\frac{c_{\varphi }}{I_x}s+\frac{k_{\varphi }-mg{h}_{\text{CG}}}{I_x}=0$$

自然频率：

$${\omega }_{n,\varphi }=\sqrt{\frac{k_{\varphi }-mg{h}_{\text{CG}}}{I_x}}$$

典型值：${\omega }_{n,\varphi }\approx 1.5\sim 2.5rad/s$（$0.25\sim 0.4\mathrm{Hz}$）

---

四、俯仰动力学的引入

4.1 俯仰运动的物理机理

俯仰角 $\theta$：车身绕 $y$ 轴的转动。

俯仰力矩平衡：

$$I_y\dot{q}+c_{\theta }q+k_{\theta }\theta =-m{a}_x{h}_{\text{CG}}-mg{x}_{\text{CG}}\sin \theta$$

其中：

• $I_y$：俯仰转动惯量
• $c_{\theta }$：俯仰阻尼系数
• $k_{\theta }$：俯仰刚度

4.2 俯仰动力学的状态方程

$$\begin{array}{rl}\dot{\theta }& =q\\ \dot{q}& =\frac{1}{I_y}\left(-m{a}_x{h}_{\text{CG}}-c_{\theta }q-k_{\theta }\theta \right)\end{array}$$

（忽略重力项，因质心接近俯仰轴）

4.3 俯仰对纵向载荷转移的影响

动态载荷转移：

$$\Delta F_{z,\text{long}}(t)=\frac{1}{L}\left(m{a}_x(t)h_{\text{CG}}+k_{\theta }\theta (t)+c_{\theta }q(t)\right)$$

---

五、垂向动力学的引入

5.1 簧载质量的垂向运动

垂向运动方程：

$$m_s{\ddot{z}}_s+c_z{\dot{z}}_s+k_z{z}_s=F_{z,\text{ext}}$$

其中：

• $m_s$：簧载质量
• $c_z$：垂向阻尼
• $k_z$：垂向刚度

5.2 四轮独立悬架模型

每轮悬架动态：

$$\begin{array}{rl}{m}_s{\ddot{z}}_s+\sum _i{c}_i({\dot{z}}_s-{\dot{z}}_{u,i})+\sum _i{k}_i(z_s-z_{u,i})& =0\\ m_{u,i}{\ddot{z}}_{u,i}+c_i({\dot{z}}_{u,i}-{\dot{z}}_s)+k_i(z_{u,i}-z_s)+k_{t,i}(z_{u,i}-z_{r,i})& =0\end{array}$$

其中：

• $z_s$：车身垂向位移
• $z_{u,i}$：非簧载质量位移
• $z_{r,i}$：路面输入
• $k_{t,i}$：轮胎垂向刚度

5.3 状态空间形式

状态变量（10DOF 四轮悬架）：

$$\mathbf{x}={\left[\begin{array}{ccccccc}{z}_s& {\dot{z}}_s& z_{u,fl}& {\dot{z}}_{u,fl}& \dots & z_{u,rr}& {\dot{z}}_{u,rr}\end{array}\right]}^{\top }$$

共 10 个状态。

---

六、完整车辆模型的耦合结构

6.1 耦合关系图

graph TD
    subgraph 车身运动
    A[纵向 v_x] --> B[侧向 v_y]
    B --> C[横摆 r]
    C --> D[侧倾 φ]
    C --> E[俯仰 θ]
    B --> D
    A --> E
    D --> F[垂向 z]
    E --> F
    end
    
    subgraph 轮胎力
    G[F_x] --> A
    H[F_y] --> B
    I[F_z] --> D
    I --> E
    I --> F
    end
    
    subgraph 载荷转移
    D --> H
    E --> G
    F --> G
    F --> H
    end

6.2 耦合强度矩阵

定义耦合强度矩阵 $\mathbf{C}$：

$$C_{ij}=\frac{|\text{第 }j\text{ DOF 对第 }i\text{ DOF 的贡献}|}{|\text{第 }i\text{ DOF 的自贡献}|}$$

典型值（乘用车）：

$$\mathbf{C}=\left[\begin{array}{cccccc}{v}_x& v_y& r& \varphi & \theta & z\\ 1& 0.05& 0.02& 0.01& 0.1& 0.01\\ 0.05& 1& 0.1& 0.2& 0.02& 0.05\\ 0.02& 0.1& 1& 0.15& 0.05& 0.02\\ 0.01& 0.2& 0.15& 1& 0.05& 0.3\\ 0.1& 0.02& 0.05& 0.05& 1& 0.1\\ 0.01& 0.05& 0.02& 0.3& 0.1& 1\end{array}\right]$$

---

七、模型选择指南

7.1 应用场景与模型复杂度

应用场景	推荐模型	自由度数	理由
路径规划（低速）	运动学模型	2	侧滑可忽略
路径规划（高速）	单轨动力学	2	侧滑显著
稳定性控制设计	单轨 + 侧倾	4	载荷转移重要
平顺性分析	完整 + 悬架	10+	垂向动态关键
极限操纵仿真	完整 + 非线性轮胎	14+	全耦合必需
翻滚分析（SUV）	完整 + 侧倾	6+	侧倾稳定性关键

7.2 模型复杂度与计算成本

模型	状态数	计算时间/步长	实时性
运动学	2	$5\mu \mathrm{s}$	✓✓✓
单轨动力学	2	$10\mu \mathrm{s}$	✓✓✓
单轨 + 侧倾	4	$30\mu \mathrm{s}$	✓✓
双轨动力学	7	$100\mu \mathrm{s}$	✓
完整 7DOF	7	$200\mu \mathrm{s}$	✓
完整 + 悬架	17	$500\mu \mathrm{s}$	×

---

八、数值示例

8.1 侧倾动力学响应

参数	符号	数值
侧倾刚度	$k_{\varphi }$	$1400\mathrm{N}\cdot m/rad$
侧倾阻尼	$c_{\varphi }$	$400\mathrm{N}\cdot \mathrm{m}\cdot s/rad$
侧倾惯量	$I_x$	$500\mathrm{kg}\cdot {\mathrm{m}}^2$
质心侧偏距	$h_{\text{CG}}$	$0.55\mathrm{m}$

阶跃侧向加速度 $a_y=0.5g$ 下的侧倾角：

$${\varphi }_{\text{ss}}=\frac{m{a}_y{h}_{\text{CG}}}{k_{\varphi }-mg{h}_{\text{CG}}}\approx 4.2^{\circ }$$

8.2 俯仰动力学响应

参数	符号	数值
俯仰刚度	$k_{\theta }$	$1200\mathrm{N}\cdot m/rad$
俯仰阻尼	$c_{\theta }$	$350\mathrm{N}\cdot \mathrm{m}\cdot s/rad$
俯仰惯量	$I_y$	$600\mathrm{kg}\cdot {\mathrm{m}}^2$

制动 $a_x=-0.5g$ 下的俯仰角：

$${\theta }_{\text{ss}}=\frac{-m{a}_x{h}_{\text{CG}}}{k_{\theta }}\approx 2.8^{\circ }$$

---

九、相关内容

• 01-双轨至单轨模型的渐近降阶理论
• 01-极限侧滑工况下的相平面稳定性理论
• 02-纵横向联合激励下的动力学耦合机理