2. 虚拟样机建模基础

虚拟样机的可信度，首先取决于建模方法是否与所研究的物理问题匹配。多体系统动力学、有限元方法和系统建模语言分别面向三类典型问题：机构运动与载荷传递、结构变形与应力分布、跨能量域的动态耦合。理解它们的适用边界，比记住软件操作更重要。

多体系统动力学：处理大范围运动与约束

多体系统动力学（MBD）把机械系统抽象为若干刚体或柔性体，通过运动副（转动副、移动副、球铰等）和力元（弹簧、阻尼、驱动、接触）连接。其核心是建立广义坐标下的运动方程：

M(q)q̈ + C(q,q̇)q̇ + Φ_q^T λ = Q_e
Φ(q,t) = 0

其中 q 为广义坐标，M 为质量矩阵，C 包含科氏力与离心力项，Φ 为约束方程，λ 为拉格朗日乘子，Q_e 为广义外力。求解器通常采用变步长隐式积分（如 BDF、Radau），因为机械系统方程常呈刚性。

以一个曲柄滑块机构为例：曲柄长 0.1 m，连杆长 0.3 m，转速 600 r/min，滑块质量 2 kg。建模时先定义转动副轴线方向，再施加驱动转速。仿真输出滑块位移、速度、加速度以及转动副反力。若加速度曲线在死点附近出现尖峰，说明约束或接触刚度设置过刚，需要检查是否引入柔性体或调整接触参数。判断结果是否合理，可先看能量守恒：无阻尼时机械能波动应小于 1%。

有限元方法：处理变形、应力与模态

有限元方法（FEM）把连续体离散为单元，通过形函数插值位移场，得到刚度矩阵 K、质量矩阵 M 和阻尼矩阵 C，求解：

M ü + C u̇ + K u = F(t)

静力分析退化为 K u = F；模态分析求解广义特征值问题 K φ = ω² M φ。工程上最需要判断的是网格收敛性：对同一结构逐步加密网格，若最大应力变化小于 5%，可认为网格足够。例如一块 200 mm × 50 mm × 5 mm 的悬臂梁，自由端受 100 N 集中力，理论挠度 δ = FL³/(3EI)。取 E = 210 GPa，I = 5×5³/12 mm⁴ ≈ 52.08 mm⁴，计算得 δ ≈ 0.61 mm。若有限元结果偏离超过 10%，应检查约束是否过约束、单元类型是否合适。

模态分析中，前六阶刚体模态应接近 0 Hz；若不为零，说明约束不足或存在未连接节点。第一阶弹性模态频率应与解析解或试验值对比，偏差通常控制在 5% 以内。

系统建模语言：跨能量域的统一表达

系统建模语言（如 Modelica）采用非因果、基于方程的建模方式，适合液压、气动、热、电与控制系统的耦合。其基本单元是“组件 + 连接”，连接器传递势变量和流变量。例如一个液压缸，一端连接泵出口压力，另一端连接负载质量，通过流量连续和力平衡自动组装方程。

与 MBD 和 FEM 相比，系统建模语言不强调几何细节，而强调能量守恒和接口一致性。它常与 MBD 联合仿真：MBD 提供机械运动与载荷，系统建模语言提供驱动与控制的动态响应。

模型简化原则与精度—成本权衡

模型简化不是随意删减，而是保留主导动态、忽略次要因素。常用原则包括：

• 频率截断：若关注 0～50 Hz 的动态，可忽略 200 Hz 以上的局部模态，用集中质量代替。
• 对称与周期简化：对称结构取一半或一个扇区，施加对称边界。
• 小变形假设：当变形小于特征尺寸的 1% 时，可线性化。
• 刚度等效：用弹簧—阻尼—质量系统替代复杂柔性体，但需保证低频模态一致。

精度与计算成本并非线性关系。网格加密一倍，自由度约增加 8 倍，计算时间可能增加 10 倍以上。因此工程上常采用“先粗后细、先刚后柔”的策略：先用刚体 MBD 快速扫描工况，再对关键部件做 FEM 校核，最后用联合仿真验证系统级响应。

实际建模时，建议先明确分析目标：若关注机构运动与驱动载荷，优先用 MBD；若关注局部应力或疲劳，必须用 FEM；若涉及液压、气动与控制耦合，用系统建模语言搭建回路，再通过接口与 MBD/FEM 联合。三者不是替代关系，而是按问题尺度分层使用。建模精度应服务于决策需求，而不是追求无限细化。