4. 机械多体系统建模与集成准备
多体动力学模型是虚拟样机中承载机械运动行为的核心载体。它把真实机构抽象为若干刚体(或柔性体)、约束它们的运动副,以及施加其上的力元,再通过求解器计算位移、速度、加速度与约束反力。多领域集成时,机械多体模型往往作为“被驱动对象”或“运动边界提供者”,与液压、控制、电机等模型交换接口变量。因此,本节的重点不是把模型建得越复杂越好,而是让模型在拓扑、参数和接口三个层面都具备可集成、可复用的条件。
一、多体动力学建模要点
建模前先明确三件事:坐标系约定、单位制、以及建模目的。坐标系建议统一采用全局坐标系描述位姿,局部坐标系固连于各刚体,方向定义一旦确定,全文与接口都不得随意更改。单位制推荐统一为 m–kg–s–N–rad,避免 CAD 导入后因 mm 与 m 混用导致质量、惯量量级错误。
刚体属性中,质量、质心位置和转动惯量决定动力学响应。对于形状规则的零件,可由 CAD 直接计算;对于复杂铸件,应确认材料密度是否与图纸一致。转动惯量必须相对质心定义,若软件要求输入相对某参考点,需用平行轴定理换算。例如某摆臂质量 2.4 kg,质心距转轴 0.18 m,绕质心惯量 0.012 kg·m²,则绕转轴惯量约为 0.012 + 2.4×0.18² ≈ 0.090 kg·m²,量级差异会直接影响摆动周期。
拓扑结构上,应优先构建“基座—运动链—末端执行器”的树状结构,避免出现闭环冗余约束。冗余约束会让求解器在数值上产生过约束,表现为约束反力异常增大或求解失败。若机构本身存在闭链(如四连杆、并联机构),应保留必要的闭环,但需检查每个闭环的自由度是否与预期一致。
二、CAD 导入与几何处理
CAD 导入不是简单“打开文件”,而是几何清理与轻量化的过程。导入前应删除与动力学无关的倒角、螺纹、装饰面和小孔,减少面片数量;对装配体,应确认各零件是否按装配位置导入,避免出现干涉或悬空。常用格式中,STEP 适合保留精确几何,Parasolid 适合与主流 CAD 互操作,STL 仅适合视觉参考,不宜用于质量属性计算。
导入后需检查三项:一是模型是否落在合理坐标范围内,避免因原点偏移导致数值奇异;二是各零件是否被识别为独立刚体,若被合并需重新拆分;三是质量属性是否与设计值一致。若 CAD 单位是 mm,导入后应整体缩放 0.001,并重新计算质量与惯量。对于需要柔性化的部件,可在导入后单独标记,后续替换为模态中性文件(MNF)。
三、运动副与力元定义
运动副定义的核心是“约束什么、释放什么”。旋转副约束三个平动和两个转动,只保留绕轴转动;移动副约束两个平动和三个转动,只保留沿轴移动;球副约束三个平动,保留三个转动。定义时需明确连接的两个刚体、作用点位置和轴向。作用点应尽量选在真实铰点,轴向用两点或全局坐标轴指定。
力元分为驱动与载荷两类。驱动可用位移、速度或加速度函数描述,例如旋转驱动可写为 omega = 2*pi*10*t,表示 10 转/秒的匀加速或匀速输入。载荷包括弹簧、阻尼、接触力和外部力。弹簧刚度与阻尼需与物理样机一致,接触力则需定义接触刚度、阻尼、摩擦系数和穿透深度。参数取值应来自设计值或实验,若为假设值,应在模型中标注为“待标定”。
四、输出接口标准化
多领域集成的关键是接口变量可识别、可映射、可复现。机械多体模型对外输出的常见变量包括:关节位移、速度、加速度、约束反力、驱动扭矩、接触力、以及关键点位置。输入变量则包括:外部驱动力/扭矩、来自控制模型的作动信号、来自液压模型的作用力。
标准化建议统一变量命名、单位、正方向与采样步长。例如关节角统一用 rad,角速度用 rad/s,力用 N,扭矩用 N·m;正方向按右手定则,并在接口文档中写明。接口变量应尽量使用物理量而非软件内部编号,避免集成时依赖特定工具。对于 FMI/FMU 导出,需确认输入输出端口与模型变量一一对应,并设置合理的初始值和单位。
五、集成前的检查与交付
完成建模后,建议按以下清单逐项确认,再进入多领域集成:坐标系与单位制是否统一;质量、质心、惯量是否与设计值一致;运动副自由度是否与机构预期一致;是否存在冗余约束或过约束;力元参数是否已标注来源或待标定;接口变量是否已按物理量命名并写明单位与正方向;导出 FMI/FMU 时端口是否与模型变量一一对应。