3. 柔性体离散化方法

柔性多体系统动力学联合仿真的第一步，是把连续弹性体转化为有限自由度的动力学模型。这个转化过程称为离散化，它直接决定了后续联合仿真的计算规模、数值稳定性和物理保真度。本节围绕有限元法、模态综合法、Craig-Bampton 方法以及模态截断准则展开，重点说明它们各自解决什么问题、如何实施、以及如何评估精度。

有限元法：从连续体到有限自由度

有限元法（FEM）将柔性体划分为若干单元，在每个单元内用形函数插值位移场，通过虚功原理或能量原理组装出整体动力学方程：

M q̈ + C q̇ + K q = F_ext + F_contact

其中 M、K 分别为质量矩阵与刚度矩阵，q 为节点自由度向量。对于三维实体单元，单个节点通常有 3 个平动自由度；若采用梁、壳单元，还会引入转动自由度。有限元法的优势在于几何适应性强、材料与边界条件处理成熟，能够直接给出应力、应变等局部场量。但它的代价是自由度数量巨大：一个中等复杂度的柔性体，节点数常在 10⁴～10⁶ 量级，直接用于多体动力学时域积分几乎不可行。因此，有限元法通常作为“高保真参考模型”，而不是联合仿真的在线模型。

模态综合法：用少量模态表征整体变形

模态综合法的基本思想是：柔性体的动态响应主要由低阶模态主导，高阶模态对整体运动的贡献较小。将物理坐标 q 变换到模态坐标 η：

q = Φ η,  Φ = [φ₁ φ₂ … φₙ]

代入原方程并左乘 Φᵀ，得到模态空间下的方程：

M̃ η̈ + C̃ η̇ + K̃ η = Φᵀ F_ext

当 Φ 取正则模态时，M̃ 和 K̃ 均为对角阵，各阶模态解耦，计算量大幅下降。但正则模态只适用于自由-自由或特定边界条件，对于多体系统中常见的“一端连接运动副、另一端自由”的柔性体，正则模态收敛较慢，需要保留大量模态才能准确描述连接处的局部变形。

Craig-Bampton 方法：固定界面模态与约束模态的结合

Craig-Bampton（C-B）方法针对上述问题，将自由度分为界面自由度 q_b（与其他部件连接处）和内部自由度 q_i。其变换关系为：

[q_b]   [ I      0    ] [q_b]
[   ] = [              ] [   ]
[q_i]   [Φ_ic  Φ_in   ] [η_n]

其中 Φ_ic 为约束模态（固定界面、逐个释放界面自由度得到的静变形），Φ_in 为固定界面正则模态。C-B 方法的物理意义非常清晰：界面自由度被完整保留，因此连接处的力和位移可以精确传递；内部变形则由少量固定界面模态近似。这样，即使只保留 10～30 阶内部模态，也能在连接界面处获得很高的精度。

在联合仿真中，C-B 方法生成的模型通常以模态中性文件（MNF）或状态空间形式导出，供多体动力学软件调用。界面节点对应实际连接点，其编号和坐标系必须与多体模型中的接口严格一致，否则会出现力传递错误或刚体位移漂移。

模态截断准则：保留多少阶才够用

模态截断是降阶的核心决策。保留阶数过少，会丢失高频动态特性，导致接触力振荡或能量异常；保留过多，则失去降阶意义。工程上常用的判断依据包括：

• 频率覆盖准则：最高保留模态频率应不低于激励最高频率的 2～3 倍。例如，若系统主要激励频率为 50 Hz，则至少保留到 100～150 Hz 的模态。
• 模态有效质量准则：在主要运动方向上，保留模态的有效质量之和应占总质量的 90% 以上。这一指标可通过有限元软件直接输出。
• 界面频响误差准则：在界面点施加单位力，比较降阶模型与全有限元模型在关注频段内的频响函数，误差应控制在 5% 以内。
• 能量范数收敛：逐步增加模态数，观察系统总能量或应变能的相对变化，当变化小于设定阈值时认为收敛。

需要强调的是，这些阈值是示例性参考，实际取值必须结合具体工况、激励频段和精度要求确定。对于接触碰撞问题，高频成分丰富，往往需要保留更多模态或采用接触局部加密策略。

精度评估：如何确认降阶模型可信

降阶模型建立后，不能直接投入联合仿真，必须先做精度评估。推荐的做法是：

1. 用全有限元模型计算一个基准工况（如自由衰减或简谐激励）的响应。
2. 用相同工况驱动降阶模型，比较界面点位移、速度、加速度时程。
3. 检查刚体模态频率是否接近零（自由-自由状态下前 6 阶应为刚体模态）。
4. 检查界面力与位移的传递是否满足牛顿第三定律，残差应在数值噪声量级。
5. 在联合仿真中监测能量守恒性，若总能量持续漂移，说明截断或积分参数需要调整。

总结来说，有限元法提供高保真基准，模态综合法实现降阶，Craig-Bampton 方法在界面精度与计算效率之间取得平衡，而模态截断准则与精度评估则是保证降阶模型可信的“质量关卡”。在联合仿真流程中，这四者构成一条完整的离散化技术链，缺一不可。