1. 课程导入与金属成形仿真概览

金属成形工艺在汽车覆盖件、航空结构件、五金件和精密锻件生产中占据核心地位。冲压、拉延、锻造、挤压等工艺的共同特征是:板料或坯料在模具作用下发生大位移、大转动和大塑性变形,局部区域可能经历减薄、增厚、回弹、起皱甚至开裂。传统试模依赖经验与反复修模,周期长、成本高。LS-DYNA 作为显式动力学有限元程序的代表,为这类强非线性、瞬态接触问题提供了高效的数值实验手段。

理解 LS-DYNA 在金属成形中的定位,关键在于把握显式时间积分与隐式求解的本质差异。显式方法采用中心差分格式,无需组装和求解整体刚度矩阵,每个时间步的计算量近似与单元数成正比,且天然适合处理接触、材料失效和大变形。其代价是时间步长受 Courant 条件限制,稳定步长约为最小单元特征长度除以材料声速。隐式方法则通过 Newton-Raphson 迭代求解平衡方程,时间步可取得较大,但每次迭代需要求解大型线性方程组,在接触状态频繁变化、材料高度非线性时收敛困难。金属成形中板料与模具的接触边界不断迁移,网格畸变严重,这正是显式方法占据优势的场景。

金属成形显式仿真:工艺对象 → 求解链路 → 结果判读 成形工艺对象 冲压 / 拉延 锻造 / 挤压 板料 / 坯料 大位移 · 大转动 · 大塑性 LS-DYNA 显式求解链路 中心差分时间积分 Δt ≤ L_min / c 接触界面搜索 罚函数 / 对称罚 材料本构更新 弹塑性 / 各向异性 单元应力回映 壳 / 实体 / 自适应 无需整体刚度矩阵 · 每步计算量线性增长 结果判读与决策 减薄率 / 开裂风险 回弹量 / 补偿方向 起皱 / 压边力调整 工艺窗口收敛判断 显式 vs 隐式:金属成形适用性对比 显式(LS-DYNA) 小步长 · 无迭代收敛 接触切换稳定 适合瞬态大变形 隐式(如 ANSYS) 大步长 · 需迭代收敛 接触突变易发散 适合静力 / 回弹 成形主过程用显式,回弹与 springback 可用隐式或显式-隐式联合 显式时间步典型量级:1e-7 ~ 1e-6 s(视网格与声速) 质量缩放可提升步长,但需控制惯性效应误差 工程应用场景与学习路径 覆盖件冲压:开裂 / 起皱 锻造成形:充型 / 折叠 挤压:流速 / 温度耦合 回弹补偿:模具型面修正 学习路径:关键字卡片 → 单单元验证 → 板料拉延 → 接触与压边 → 回弹补偿 每阶段以可复现算例闭环,逐步建立参数-结果映射直觉

从工程应用看,LS-DYNA 显式仿真主要覆盖四类场景。第一类是覆盖件冲压与拉延,关注减薄率、开裂、起皱和压边力窗口;第二类是锻造与挤压,关注金属充型、折叠缺陷、模具应力和温度耦合;第三类是回弹预测与模具型面补偿,通常先显式完成成形,再隐式或显式求解回弹;第四类是高速冲击、旋压、电磁成形等瞬态问题,显式几乎是唯一可行选择。

学习路径建议按“关键字卡片 → 单单元验证 → 板料拉延 → 接触与压边 → 回弹补偿”推进。关键字卡片阶段要能独立写出 *CONTROL_TIMESTEP、*MAT_PIECEWISE_LINEAR_PLASTICITY、*CONTACT_FORMING_SURFACE_TO_SURFACE 等核心卡片;单单元验证用于确认材料曲线、单位制和沙漏控制是否正确;板料拉延算例则把网格、接触、压边力和时间步串起来。每阶段都应以可复现算例闭环,而不是只看理论。

显式仿真的时间步由最小单元尺寸和材料声速共同决定。以钢为例,声速约 5000 m/s,若最小单元尺寸为 1 mm,则稳定步长约 2e-7 s。若把步长放大到 1e-6 s,就必须通过质量缩放或增大最小单元尺寸实现,但质量缩放会引入附加惯性,需检查动能与内能之比是否仍在可接受范围。

一个典型的板料拉延教学算例可以这样设置:板料 300 mm × 200 mm × 1 mm,采用 Belytschko-Tsay 壳单元,网格尺寸 2 mm;材料为低碳钢,密度 7.85e-9 t/mm³,弹性模量 210 GPa,泊松比 0.3,屈服强度 180 MPa,硬化曲线用 *DEFINE_CURVE 给出;模具和压边圈设为刚体,压边力 200 kN,冲头速度 5 m/s。求解后重点观察厚度云图中减薄率是否超过 20%、板料边缘是否出现起皱波,以及冲头力-位移曲线是否出现异常抖动。若减薄率超标,可调整压边力、拉延筋或坯料形状;若起皱,则需提高压边力或优化材料流动。