4. 板料成形仿真单元技术

板料成形仿真的精度，往往不是被材料模型“卡住”，而是被单元技术“拖住”。同一套材料参数、同一副模具型面，仅改变单元类型、积分方案或厚度方向积分点数量，回弹预测值就可能相差 20%～40%，破裂位置也可能从圆角根部“跳”到侧壁中部。本节围绕壳/实体单元选择、积分与沙漏、厚度方向积分点、单元尺寸四条主线，说明它们各自影响什么、如何配置、以及结果如何判读。

一、壳单元与实体单元：先看厚度方向应力是否需要显式求解

板料厚度通常为 0.6～2.0 mm，而面内尺寸可达数百毫米，厚度/面内比约 1/100～1/500，属于典型薄壁结构。壳单元把厚度方向应力做简化假设，用中面位移加转动自由度描述变形，计算量小、适合覆盖件整体成形；实体单元在厚度方向布置多层单元，可显式求解厚度应力与剪切，适合翻边、包边、局部镦厚、厚板成形等场景。

选择判据可归纳为三点：其一，若厚度方向正应力与横向剪应力对结果影响小（如普通覆盖件拉延），优先壳单元；其二，若存在明显厚度压缩、局部增厚或剪切主导变形（如厚板翻边、包边、镦粗类工序），应使用实体单元或壳-实体混合；其三，若只关心回弹趋势与破裂风险区，壳单元配合足够厚度积分点通常已够用。工程上常见做法是：拉延工序用壳单元，翻边/整形等局部工序切换到实体单元，并在同一模型中保持材料与接触参数一致，避免“换单元就换结论”。

二、积分方案与沙漏控制：减缩积分的代价必须补回来

壳单元常用积分方案有两类：完全积分与减缩积分。完全积分在厚度方向与面内均取足够积分点，刚度偏硬，不易出现沙漏，但计算量大，且对剪切锁死较敏感；减缩积分每个方向少取一个积分点，计算效率高、抗剪切锁死好，但会引入零能模式，即沙漏。沙漏表现为单元呈棋盘状、无物理意义的变形，能量不守恒，回弹与破裂预测都会失真。

控制沙漏的常用手段包括：采用沙漏控制算法（如刚度型或黏性型）、适当增加面内积分点、细化网格、或在关键区域改用完全积分。判断沙漏是否被有效抑制，可观察沙漏能占总内能的比例，工程经验上该比例宜控制在 5% 以内，超过 10% 时结果可信度明显下降。若沙漏能持续偏高，优先检查网格是否过粗、减缩积分是否用于弯曲主导区，而不是盲目加大沙漏控制系数——过强的沙漏控制会人为增加刚度，使回弹被低估。

三、厚度方向积分点：回弹预测的“分辨率”

壳单元在厚度方向用积分点近似应力沿厚度的分布。积分点太少（如 3 点），厚度方向弯曲应力梯度被“抹平”，回弹驱动力矩被低估，回弹预测偏小；积分点增加到 5～7 点，厚度方向应力分布可分辨，回弹预测趋于稳定。工程经验是：普通覆盖件回弹分析建议 5 点起步，高强钢或回弹敏感件用 7 点；继续增加到 9 点以上，精度提升有限，计算量却明显上升。

需要区分的是：厚度方向积分点数量影响的是“厚度方向应力分辨率”，而面内积分方案影响的是“面内刚度与沙漏”。两者不能互相替代。若回弹预测始终偏小，先检查厚度积分点是否足够，再检查材料硬化模型与卸载路径，而不是一味加密面内网格。

四、单元尺寸：破裂预测看圆角，回弹预测看整体

单元尺寸对两类预测的影响机制不同。破裂预测对局部应变梯度敏感：凹模圆角、凸模圆角、拉延筋入口处应变集中，若单元尺寸大于圆角半径的 1/3～1/2，峰值应变会被平均掉，破裂风险被低估。因此圆角区单元尺寸通常控制在 0.5～1.0 mm，且圆角处至少布置 4～6 层单元。回弹预测对整体弯曲刚度与残余应力分布敏感：单元过粗会高估弯曲刚度、低估回弹；单元过细则计算量上升，且接触搜索成本增加。覆盖件回弹分析常用 2～5 mm 单元，并在拉延筋、圆角、切边线附近局部细化。

一个实用的检查方法是做单元尺寸收敛性验证：以 4 mm、2 mm、1 mm 三档网格分别计算同一零件，观察回弹角与最大减薄率的变化。若从 2 mm 到 1 mm 结果变化小于 5%，可认为 2 mm 已收敛；若变化仍大于 10%，说明当前网格尚未捕捉到关键梯度，需继续细化或检查圆角区单元质量。

综合来看，单元技术配置应形成一条清晰的判断链：先判断厚度方向应力是否重要，决定壳/实体；再根据弯曲主导程度选择积分方案并监控沙漏能；然后按回弹敏感度确定厚度方向积分点；最后按破裂风险区与回弹整体变形分别确定单元尺寸，并用收敛性验证确认结果稳定。四者协同，才能让回弹与破裂预测同时可信。