1. 金属成形仿真概览

金属成形工艺——冲压、锻造、挤压、拉延——本质上都是"大变形、强非线性、短时程"的力学过程。板料在几十毫秒内被拉深成杯形件，坯料在锻锤下瞬间充满模腔，金属在挤压筒内经历剧烈塑性流动。这类问题的共同特征是：材料进入塑性大变形阶段、接触边界不断变化、几何构型持续更新，且惯性效应与应变率效应不可忽略。这正是显式动力学方法的主场。

为什么是显式动力学

隐式方法求解的是平衡方程 K·u = F，每一步都要组装并求逆整体刚度矩阵。当成形过程出现网格严重畸变、接触状态频繁切换时，刚度矩阵病态甚至奇异，迭代收敛困难。显式方法则直接对运动方程做时间积分：

M·a = F_ext - F_int
v_(n+1/2) = v_(n-1/2) + a_n · Δt
u_(n+1)   = u_n + v_(n+1/2) · Δt

其中 M 为集中质量矩阵，求逆只是逐节点除法，无需迭代。代价是条件稳定：时间步长必须小于临界值 Δt_crit ≈ L_min / c，L_min 是最小单元特征长度，c 是材料声速。对钢材料，c ≈ 5000 m/s，若最小单元为 1 mm，则 Δt_crit ≈ 0.2 μs。冲压过程持续 50 ms，意味着约 25 万步——这正是显式方法"以步数换稳定性"的典型代价，也是它适合短时程、高非线性问题的原因。

四类典型工艺的仿真定位

冲压/拉延：关注板料减薄率、起皱、破裂与回弹。显式求解成形阶段，回弹常需切换到隐式或采用显式-隐式联合流程。锻造：关注充型完整性、折叠缺陷、模具应力与锻件温度场，常需耦合热-力分析。挤压：关注金属流动均匀性、挤压力行程曲线与模具磨损，网格重划分需求高。拉延：与冲压高度重叠，但更强调压边力、拉延筋与材料各向异性对成形极限的影响。

LS-DYNA求解流程的实际操作

一个完整的板料拉延仿真，通常按以下顺序推进。首先是几何与网格：板料用壳单元（如 *SECTION_SHELL，厚度方向 5 个积分点），模具用刚性壳或实体单元。板料网格在圆角区域需细化到 1–2 mm，否则无法捕捉弯曲变形。接着是材料定义：常用 *MAT_PIECEWISE_LINEAR_PLASTICITY（24 号材料），输入弹性模量 210 GPa、泊松比 0.3、密度 7.85e-9 t/mm³，以及真实应力-塑性应变曲线。若考虑应变率效应，可改用 123 号或 124 号材料模型。

接触设置是成形仿真的成败关键。板料与模具之间用 *CONTACT_FORMING_ONE_WAY_SURFACE_TO_SURFACE，静摩擦系数取 0.12、动摩擦系数取 0.10 是冷冲压的常见起点。压边圈通过 *BOUNDARY_PRESCRIBED_MOTION_RIGID 施加恒定压边力，凸模则用速度曲线控制下行。

求解控制方面，*CONTROL_TIMESTEP 中的 DT2MS 设为负值可启用质量缩放，把最小时间步"抬"到目标值。例如目标步长 1 μs，而实际最小单元对应 0.3 μs，则通过缩放将步长提升，代价是引入附加质量。质量增加应控制在总质量的 1% 以内，否则惯性效应会污染结果。*CONTROL_TERMINATION 设定终止时间，*DATABASE_BINARY_D3PLOT 控制输出帧间隔。

结果判读与常见异常

求解完成后，先看 d3hsp 文件中的能量摘要：沙漏能应小于内能的 5%，滑移能不应出现负值，总能量应基本守恒。再看 d3plot 云图：板料减薄率超过 25% 通常预示破裂风险，厚度增加区域可能起皱。若出现单元负体积或接触穿透，需回到网格与接触参数排查。

文件体系与工作习惯

LS-DYNA 的输入是纯文本 K 文件，可用文本编辑器直接修改，也可由 ANSA、HyperMesh、LS-PrePost 等前处理器生成。一个良好的工作习惯是：主 K 文件只保留 *INCLUDE 语句，把网格、材料、接触、控制卡分文件管理。这样修改材料参数时不必翻动上万行网格数据，也便于版本对比与团队协作。