2. 多体系统建模基础
多体系统动力学的核心任务，是把一台真实机械装置翻译成计算机能够求解的数学模型。这个翻译过程不是简单地把零件画进软件，而是要在物理对象与数学变量之间建立一套严格的对应关系。本节围绕六个基本要素展开：刚体、铰链、约束、力元、坐标系与广义坐标，并说明它们如何共同构成从物理系统到数学模型的完整链条。
一、刚体：被抽象掉变形的质量载体
刚体是多体系统中最基本的运动单元。所谓“刚”，是指物体内部任意两点间距离在运动过程中保持不变，即忽略弹性变形。对于大多数机构动力学问题，构件本身的变形远小于整体运动幅度，因此刚体假设是合理且高效的。
一个刚体在数学上需要描述三类属性：质量 m、质心位置以及惯性张量 J。惯性张量通常先在随体坐标系下定义，再通过方向余弦矩阵转换到全局坐标系。以一根均质矩形杆为例，长 L=0.5 m、截面 0.04 m×0.02 m、质量 2 kg，其绕质心垂直于杆长方向的转动惯量约为：
J = m·L² / 12 = 2 × 0.5² / 12 ≈ 0.0417 kg·m²这个数值直接决定角加速度对力矩的响应。若把杆件误设为点质量，转动惯量为零，动力学响应会完全失真。因此刚体建模的第一件事，是确认质量、质心与惯性参数是否与真实构件一致。
二、铰链与约束：限制相对运动的规则
铰链（也称运动副）描述两个刚体之间允许的相对运动形式。旋转铰只允许绕某一轴转动，剩余 5 个自由度被限制；滑移铰只允许沿某一轴平移；球铰允许三个转动自由度；万向节则允许两个转动自由度。约束的数学表达通常写成位置约束方程 Φ(q,t)=0，其时间导数给出速度约束，二阶导数给出加速度约束。
约束分为完整约束与非完整约束。完整约束可写成坐标的代数方程，例如旋转铰的轴对齐条件；非完整约束只能写成速度的不可积关系，例如纯滚动轮。多体动力学中绝大多数机械铰链属于完整约束，可用拉格朗日乘子法统一处理。
约束数量与自由度之间存在直接关系：若系统有 n 个广义坐标、m 个独立完整约束，则系统自由度为 n−m。例如平面四杆机构有 3 个活动构件，每个构件 3 个自由度，共 9 个广义坐标；4 个旋转铰各提供 2 个平面约束，共 8 个约束，自由度 = 9−8 = 1，与曲柄摇杆机构的实际单自由度一致。
三、力元：驱动与耗散的来源
力元描述系统内部或外部的作用力，包括弹簧、阻尼器、驱动器、重力、接触力等。力元不改变系统自由度，但决定运动方程右端项。常见力元包括：线性弹簧 F=k(x−x₀)，线性阻尼 F=c·v，以及电机/液压缸输出的驱动力或驱动力矩。
力元可以作用在刚体质心，也可以作用在任意标记点。作用点位置会影响等效力矩，因此建模时必须明确力元的作用点与方向。例如在旋转铰处施加驱动力矩 T=5 N·m，若转动惯量为 0.0417 kg·m²，忽略阻尼时角加速度约为 120 rad/s²；加入阻尼 c=0.1 N·m·s/rad 后，稳态角速度会趋于 T/c=50 rad/s。
四、坐标系与广义坐标：描述运动的语言
多体系统至少涉及两类坐标系：全局坐标系（惯性系，固定于地面）和随体坐标系（固定于各刚体，随刚体运动）。全局坐标系用于统一度量位置与方向，随体坐标系用于定义刚体自身的惯性属性与连接点。
广义坐标是描述系统位形的独立变量集合。常见选择包括：笛卡尔坐标加欧拉角、自然坐标（质心坐标加方向余弦）、相对坐标（铰链转角）等。广义坐标不唯一，但必须能唯一确定系统位形，且数量等于系统自由度。
以平面曲柄滑块机构为例，若采用相对坐标，可取曲柄转角 θ 为广义坐标，滑块位移 x 由几何约束 x = r·cosθ + √(l²−r²sin²θ) 确定。这样只需一个变量即可描述整个系统，方程规模最小。若采用笛卡尔坐标，则需要曲柄质心、曲柄转角、连杆质心、连杆转角、滑块位移共 5 个坐标，再附加 4 个约束方程，方程规模更大但约束处理更规整。
五、从物理系统到数学模型的完整链条
把上述要素串起来，多体系统建模遵循一条清晰的映射链：物理构件 → 刚体（质量、质心、惯性张量）；物理连接 → 铰链与约束方程；物理作用 → 力元；物理位形 → 广义坐标；物理运动 → 约束方程与动力学方程。最终得到形如 M(q)q̈ + Φ_qᵀλ = Q(q,q̇,t) 的微分代数方程组，其中 M 为质量矩阵，Φ_q 为约束雅可比，λ 为拉格朗日乘子，Q 为广义力。
六、建模实践中的关键判断
在实际建模中，最容易出错的不是方程推导，而是要素选择与参数一致性。以下几点值得特别注意：
- 惯性参数必须与随体坐标系匹配。若惯性张量在质心坐标系下给出，而建模时把参考点设在铰链处，需要用平行轴定理平移：J_O = J_c + m·d²。
- 约束不能重复也不能遗漏。同一方向若被两个铰链同时限制，会导致约束雅可比行相关，求解器报奇异。建模后应检查约束雅可比的秩。
- 力元方向要与坐标系定义一致。弹簧力沿两标记点连线方向时，需明确该方向在全局系下的分量；若方向写反，系统会发散。
- 广义坐标数量应等于自由度。若采用非独立坐标，必须同时给出约束方程，不能只给坐标不给约束。
总结来看，多体系统建模的本质是建立一套一一对应的映射：每个物理构件对应一个刚体，每个物理连接对应一组约束方程，每个物理作用对应一个力元，每个独立运动对应一个广义坐标。坐标系提供度量基准，广义坐标提供变量语言，约束方程提供运动规则，力元提供驱动与耗散。六者齐备，物理系统才真正转化为可求解的数学模型。