2. 多体动力学建模基础
多体动力学仿真的本质，是把一台真实机械系统翻译成计算机能求解的数学对象。这个翻译过程不是“画几个方块加箭头”，而是需要明确五类基本要素：刚体、约束、力元、坐标系、自由度，再由它们组装出运动学与动力学方程。理解这套映射规则，是后续台架对标能否成立的前提——如果模型里的约束和力元与物理系统不一致，仿真曲线再漂亮也没有对标价值。
一、五类基本要素及其物理含义
刚体是建模的最小惯性单元。每个刚体需要定义质量 m、质心位置、以及相对于质心坐标系的惯性张量 I（3×3 对称矩阵，含 Ixx、Iyy、Izz 及三个惯性积）。例如一根 1.2 m 长、外径 60 mm、壁厚 4 mm 的钢制摆臂，质量约 6.7 kg，若按简化杆件估算，绕质心横轴惯量约 0.8 kg·m²。刚体本身不产生运动，它只是承载质量和惯量的“容器”。
约束限制刚体之间的相对运动，是自由度的“减法器”。常见类型包括：旋转铰（释放 1 个转动自由度）、平移铰（释放 1 个移动自由度）、球铰（释放 3 个转动自由度）、万向节、圆柱铰、固定副、齿轮副与耦合副。约束的数学形式是代数方程 Φ(q,t)=0，其雅可比矩阵 Φq 决定约束方向。
力元是自由度的“激励源”，包括弹簧、阻尼、衬套（六向刚度阻尼）、接触力、轮胎力、驱动与负载力矩、重力等。力元不直接约束运动，而是通过广义力 Q 进入动力学方程右端。
坐标系提供描述基准。工程上通常区分：全局惯性系（地面）、局部体坐标系（固连于刚体）、标记点坐标系（用于定义铰点、力作用点、传感器测点）。所有铰点位置、力矢量方向都必须明确声明在哪个坐标系下表达，否则正负号会出错。
自由度是系统独立运动参数的数量。对 N 个刚体、每个 6 个空间自由度，共 6N 个广义坐标；若施加 m 个独立约束，则系统自由度 DOF = 6N − m。例如双横臂悬架：上摆臂、下摆臂、转向节、轮毂共 4 个刚体（24 自由度），若用 4 个旋转铰加 1 个球铰约束，剩余自由度约 2～3 个，正好对应轮跳与转向两个独立运动。
二、运动学与动力学方程
运动学层面，位置约束 Φ(q,t)=0 对时间求导得到速度约束 Φq·q̇ = −Φt，再求导得到加速度约束 Φq·q̈ = γ。给定驱动后，可先解位置、再解速度、最后解加速度，这就是运动学分析的“三步走”。
动力学层面，常用拉格朗日方程或第一类拉格朗日方程（含拉格朗日乘子 λ）：
M(q)·q̈ + Φqᵀ·λ = Q(q, q̇, t)
Φ(q, t) = 0其中 M 为质量矩阵，Q 为广义力，Φqᵀ·λ 为约束反力。求解器（如 Adams、RecurDyn、Simpack、Simscape Multibody）内部通常把上式转为微分代数方程（DAE）或通过坐标缩减转为常微分方程（ODE）后积分。理解这一点，就能明白为什么“约束冗余”会导致求解失败——约束方程线性相关时，Φq 秩亏，λ 不唯一。
三、从物理系统到仿真模型的映射方法
映射不是一步到位，而是按“拓扑—参数—验证”三层推进：
- 拓扑映射：先画物理系统的机构简图，标出所有刚体、铰点、力作用点，确定哪些是运动件、哪些是接地件。这一步决定约束类型和数量。
- 参数映射：把 CAD 或实测得到的质量、惯量、铰点坐标、衬套刚度阻尼、弹簧预载填入模型。参数来源必须可追溯，建议建立“参数—来源—单位—置信度”对照表。
- 验证映射：先做静平衡或运动学检查（如轮跳时轮心轨迹是否合理），再做动力学试算，最后才进入台架对标。
四、一个可复用的建模检查清单
| 检查项 | 判断依据 | 常见错误 |
|---|---|---|
| 刚体数量与惯量 | 运动件是否全部覆盖，惯量单位是否为 kg·m² | 漏掉小质量件导致高频模态缺失 |
| 约束类型与数量 | DOF = 6N − m 是否与物理一致 | 过约束导致求解器报冗余 |
| 铰点坐标 | 是否声明坐标系，正负号是否统一 | 左右侧坐标混用 |
| 力元方向 | 弹簧/阻尼/衬套的作用线是否与物理一致 | 衬套局部坐标系装反 |
| 驱动与负载 | 驱动函数、负载力矩是否可复现 | 驱动加在错误自由度上 |
下面这张图把上述映射关系压缩成一张“物理系统 → 仿真模型”的结构图，左侧是物理对象，中间是五类建模要素，右侧是求解输出，便于在建模时逐项核对。
把这张图当作建模时的“对照表”使用：每建立一个刚体、每添加一个约束或力元，都回到图中确认它对应物理系统的哪一部分、进入哪一类方程、最终影响哪一项输出。当仿真与台架出现偏差时，也能按图索骥——先查约束与自由度是否一致，再查力元参数，最后才怀疑积分器设置。这样，多体动力学建模才不是“调参玄学”，而是一条可追溯、可复现的工程链路。