2. 多体动力学建模数学基础

传动系统多体动力学仿真，本质上是在计算机里“重演”齿轮、轴、轴承、壳体等部件在力与运动作用下的相互作用。要让计算机能算，必须先把机械系统翻译成一套数学语言：用什么量描述位置，用什么方程描述连接，用什么规律描述运动与受力的因果关系，最后用什么算法把方程解出来。本节围绕这四件事展开：广义坐标、约束方程、拉格朗日方程、数值求解方法。

一、广义坐标：把机械系统的位形压缩成一组数

描述一个刚体在空间中的位置和姿态，最直观的是笛卡尔坐标加欧拉角，共 6 个量。若传动系统有 n 个刚体，就需要 6n 个坐标。但机械系统里大量部件被铰链、轴承、齿轮啮合约束，真正独立的自由度远小于 6n。广义坐标就是用来描述系统位形的一组独立或非独立的参数，记作 q = [q₁, q₂, …, q_m]^T。

在传动系统建模中，常见选择是：轴系用转角 θ 作为广义坐标，齿轮副用主动轮转角与从动轮转角，行星排用太阳轮、齿圈、行星架的转角，轴承支撑处用平移位移 x、y、z 与转动角。选择原则是“能少不多、能独立不冗余”：能用转角描述转动就不要引入完整欧拉角，能用一个啮合约束消去一个坐标就不要保留两个独立坐标。广义坐标选得好，方程规模小、数值条件好；选得差，约束方程复杂、求解容易病态。

二、约束方程：把“不能动的方向”写成等式

约束方程描述系统中各坐标之间必须满足的几何或运动关系，一般写成 Φ(q, t) = 0。传动系统里最常见的约束有三类。

• 转动副/轴承约束：限制两个部件之间除绕轴转动外的所有相对运动，例如轴与轴承座之间只允许绕 z 轴相对转动。
• 齿轮啮合约束：用节圆半径比建立两齿轮转角关系，如 r₁θ₁ = r₂θ₂（外啮合时考虑转向相反）。
• 驱动约束：给定输入轴转速或转角随时间变化，如 θ_in = ωt，把系统运动“驱动”起来。

约束方程对时间求导可得速度约束 Φ_qq̇ = −Φ_t，再求导得加速度约束 Φ_qq̈ = γ，其中 Φ_q 是约束雅可比矩阵。雅可比矩阵的每一行代表一个约束对各广义坐标的敏感程度，是后续动力学方程组装的核心数据。若约束之间线性相关，雅可比矩阵秩亏，系统会出现多余约束，仿真中表现为求解失败或结果跳变，因此建模时要检查约束独立性。

三、拉格朗日方程：从能量角度建立动力学方程

牛顿-欧拉法需要逐部件画受力图，对传动系统这种多部件、多约束的系统非常繁琐。拉格朗日方程从能量出发，自动消去理想约束力，更适合多体系统。其一般形式为：

d/dt(∂T/∂q̇) − ∂T/∂q = Q + Φ_qᵀλ

其中 T 是系统动能，Q 是广义力，Φ_q^Tλ 是约束力项，λ 是拉格朗日乘子，物理意义是约束反力。对传动系统，动能 T 通常写成 ½q̇^TM(q)q̇，M 是质量矩阵，包含各轴转动惯量、齿轮惯量、行星排等效惯量。广义力 Q 包括驱动力矩、负载力矩、齿轮啮合力、轴承阻尼力等。

以一对直齿轮为例，设主动轮转角 θ₁、从动轮转角 θ₂，转动惯量 J₁、J₂，啮合约束 r₁θ₁ = r₂θ₂。动能 T = ½J₁θ̇₁² + ½J₂θ̇₂²，广义力 Q₁ = T_in，Q₂ = −T_load。代入拉格朗日方程并加入约束乘子，可得：

J₁θ̈₁ + r₁λ = T_in
J₂θ̈₂ − r₂λ = −T_load
r₁θ₁ − r₂θ₂ = 0

这就是齿轮传动最简化的动力学方程。λ 对应啮合力，正负号取决于约束方程写法。实际建模中还要加入时变啮合刚度、齿侧间隙、轴承支撑刚度，方程会扩展为更复杂的微分-代数方程组。

四、动力学方程数值求解：从方程到曲线

拉格朗日方程最终形成的是微分-代数方程组（DAE）：

M q̈ + Φ_qᵀλ = Q
Φ(q, t) = 0

这类方程不能直接解析求解，必须用数值方法逐步推进。传动系统仿真中常用的方法有三类。

• 显式积分：如四阶龙格-库塔（RK4），每步计算 q̈ 后直接更新 q 和 q̇。实现简单、单步计算量小，但稳定性受步长限制，适合非刚性、频率不高的系统。
• 隐式积分：如 BDF、Newmark-β，需要求解非线性方程组，每步计算量大，但稳定性好，适合含高刚度轴承、齿轮啮合刚度的刚性系统。
• 约束稳定化：数值积分会累积约束误差，导致 Φ(q,t) 逐渐偏离零，出现“约束漂移”。常用 Baumgarte 稳定化把 Φ 和 Φ̇ 反馈到加速度约束中，或每步做位置/速度投影，把 q 拉回约束流形。

实际仿真中，步长选择是关键。以齿轮啮合为例，啮合刚度可能达到 10⁸ N/m 量级，系统最高关注频率若为 5 kHz，则步长通常取最高频率周期的 1/10～1/20，即 10～20 μs。若用显式积分，步长还要更小；若用隐式积分，可适当放大，但需检查能量误差和约束漂移。判断结果是否可信，可以观察三条曲线：约束方程残差是否保持在 10⁻⁶ 量级以内、能量是否在无阻尼工况下近似守恒、啮合力是否出现非物理的高频振荡。若残差持续增大或啮合力剧烈跳变，应减小步长、改用隐式方法或检查约束独立性。

把以上四步串起来：先选广义坐标把系统位形参数化，再用约束方程描述部件之间的连接关系，然后用拉格朗日方程组装出动力学方程，最后用合适的数值方法求解并检查约束残差与能量行为。这套数学基础决定了后续齿轮、轴系、行星排建模的精度和效率，也是理解仿真软件内部求解逻辑的关键。