3. 仿真模型构建

在 LS-DYNA 仿真与试验对标中，模型构建质量直接决定后续对标结论是否可信。几何清理、网格划分、单元类型选择与对标目标区域确定，这四件事并非彼此独立的“建模步骤”，而是一条相互约束的决策链：几何清理决定网格能否顺利生成，网格尺寸与质量决定单元类型能否发挥预期精度，而单元类型又反过来影响对标目标区域上可提取的物理量。本节围绕这条决策链展开，重点说明每一步“实际怎么做”以及“做到什么程度可以进入下一步”。

一、几何清理：先明确“为谁清理”

几何清理的目标不是把 CAD 修得漂亮，而是让几何适合显式动力学网格划分与接触定义。导入 LS-DYNA 前处理环境（如 LS-PrePost、HyperMesh、ANSA）后，建议按以下顺序处理：

• 缝合自由边与重复面：容差通常取 0.01～0.1 mm，视零件最小特征尺寸而定。若最小圆角半径为 0.5 mm，容差不应超过 0.05 mm，否则会误缝合真实间隙。
• 删除微小特征：直径小于 2 倍目标网格尺寸的孔、宽度小于 1 倍网格尺寸的凸台与倒角，通常简化为平面或忽略。例如目标网格 3 mm 时，直径 5 mm 以下的孔可填平。
• 中面抽取：薄壁件（厚度/跨度 < 1/10）优先抽中面，用壳单元建模；厚实体件保留实体几何。
• 修复拓扑错误：检查并消除尖角、重叠面、非流形边，避免网格划分时出现负体积或畸变单元。

清理完成后应做一次“几何健康检查”：自由边数量、重复节点、最小面尺寸。若仍有大量自由边，说明缝合不彻底，继续划分只会把问题传递到网格阶段。

二、网格划分：尺寸、质量与局部加密

网格尺寸由对标目标区域的特征尺度和关心的最高频率共同决定。显式动力学中，稳定时间步长近似为 Δt ≈ L/c，其中 L 为最小单元特征长度，c 为材料声速。以钢为例，c ≈ 5000 m/s，若最小单元 2 mm，则 Δt ≈ 4×10⁻⁷ s。若希望总计算时间控制在合理范围，网格不能无限细化。

实际划分时建议遵循：

• 全局尺寸：取结构最小关注特征尺寸的 1/3～1/5。例如关注 15 mm 宽的筋板，全局网格取 3～5 mm。
• 局部加密：在对标目标区域（如碰撞接触区、应力集中区、应变片粘贴区）加密到全局尺寸的 1/2～1/3，并保证过渡区单元尺寸变化不超过 2 倍，避免刚度突变。
• 质量指标：壳单元长宽比 < 5，翘曲角 < 15°，雅可比 > 0.6；实体单元雅可比 > 0.5，最小内角 > 20°。出现负体积必须立即修正。

网格质量检查应作为独立环节执行，而不是“划完就用”。建议输出一份质量统计：单元总数、最小/最大尺寸、雅可比分布、负体积数量。若负体积不为零，后续计算必然报错。

三、单元类型选择：壳、实体与积分方式

LS-DYNA 中常用单元类型及适用场景如下：

选择原则：薄壁件优先壳单元，厚向尺寸接近面内尺寸时用实体单元。壳单元需注意厚向积分点数量，3 点可捕捉弯曲，5 点更适合大变形。实体单元若出现剪切自锁，可改用 ELFORM=-1 或 -2，但计算成本上升。对标目标区域若涉及厚度方向应力梯度，必须用实体单元或壳的厚向多点积分。

四、对标目标区域确定：从试验测点反推模型关注区

对标目标区域不是“整个模型”，而是试验中真正有测点、有数据、能比较的区域。确定方法：先列出试验测点（应变片位置、加速度计安装点、位移传感器方向），再在模型上标出对应位置，最后检查该区域网格是否足够细、单元类型是否能输出对应物理量。

例如试验在 B 柱中部贴了应变片，模型就应在该处加密壳网格，并确保输出 *DATABASE_ELOUT 中对应单元的应变。若试验测的是加速度，则目标区域应包含加速度计安装点附近的节点，并检查该处是否因网格过粗导致局部模态失真。目标区域确定后，应形成一份“对标映射表”：试验测点编号、模型对应单元/节点、输出变量、允许误差范围。

五、构建完成后的自检清单

在提交计算前，建议逐项确认：几何自由边是否为零；网格负体积是否为零；目标区域网格尺寸是否达到试验测点分辨率要求；单元类型是否与厚度方向应力梯度匹配；对标目标区域是否已建立命名 Set 并确认输出变量。任何一项不满足，都应回到对应环节修正，而不是带着问题进入求解。模型构建阶段多花一小时检查，往往能避免后续数小时的无效计算与错误对标结论。