4. 单元类型与网格策略对失效模拟的影响

在LS-DYNA中做失效与断裂分析，单元类型和网格策略往往比材料参数本身更能决定结果的可信度。同一组MAT_ADD_EROSION或MAT_JOHNSON_COOK参数，换一种单元公式或网格尺寸，裂纹起裂位置、扩展路径甚至是否起裂都可能完全不同。本节围绕实体单元、壳单元、梁单元、SPH与EFG五类常用离散方式，说明它们各自适合的失效场景，并讨论网格尺寸敏感性、自适应网格与网格删除技术的关键控制点。

一、五类单元的失效适用性

实体单元（solid element）是三维断裂分析的主力。LS-DYNA中常用的ELFORM=1（常应力单点积分）计算最快，但容易出现沙漏，需要配合IHQ=4等沙漏控制；ELFORM=2为全积分，无沙漏但易发生体积锁死，在近不可压材料中会显著偏硬；ELFORM=-1（单点积分+物理沙漏控制）和ELFORM=-2（全积分+假定应变）是冲击与失效分析中较稳妥的选择。实体单元能自然表现三轴应力状态，因此对剪切主导的韧性断裂、孔洞扩张和层裂都较合适。

壳单元（ELFORM=2 Belytschko-Tsay、ELFORM=16全积分等）适合薄壁结构的撕裂、褶皱和冲压开裂。壳单元厚度方向只有一个积分点（NIP默认5），厚度方向应力近似为零，因此对以厚度减薄为判据的成形极限图（FLD）类失效很有效，但对需要厚度方向应力梯度的层裂问题不适用。梁单元（ELFORM=1 Hughes-Liu、ELFORM=6等）主要用于框架、钢筋、索缆的拉断与屈曲失效，其截面积分点（如NIP=5的矩形截面）决定塑性铰和断裂能否被捕捉。

SPH（*SECTION_SPH）是无网格粒子法，适合大变形、碎片飞散、射流侵彻等网格会严重畸变的问题，粒子间距相当于网格尺寸，失效通过粒子删除或损伤累积实现。EFG（*SECTION_EFG）为无网格伽辽金法，适合裂纹扩展路径未知、需要避免网格重划的场合，但计算成本高，接触处理也更复杂。工程上常见做法是：主体结构用实体或壳单元，局部高风险区嵌入SPH或EFG，通过*CONTACT_TIED_SURFACE_TO_SURFACE或*CONSTRAINED_INTERPOLATION连接。

二、网格尺寸敏感性：为什么结果会随网格变化

失效模拟对网格尺寸的敏感性来自两方面。其一是应变局部化：当单元特征尺寸 le 大于断裂过程区尺寸时，软化段无法被正确分辨，裂纹会被“抹平”到整个单元带，导致能量耗散被高估或低估。其二是删除判据的离散性：以等效塑性应变EPS为判据时，粗网格单元积分点少，达到临界应变的时刻偏晚，峰值力偏高、吸能偏低。工程经验是，韧性断裂的 le 宜控制在0.2～0.5 mm量级（视板厚与材料而定），并至少做粗、中、细三级网格对比，观察力-位移曲线的峰值与软化段是否趋于稳定。

为降低敏感性，可采用正则化手段：在MAT_ADD_EROSION中同时给出MXEPS与断裂能Gf，或使用*MAT_GURSON、*MAT_JOHNSON_COOK配合*MAT_ADD_DAMAGE_GISSMO，让损伤演化与单元尺寸挂钩。GISSMO中通过*MAT_ADD_DAMAGE_GISSMO的DMGEXP和FADEXP控制损伤累积，并用*DEFINE_CURVE定义失效应变随应力三轴度的变化，能在一定程度上补偿网格依赖性。

三、自适应网格与网格删除的配合

自适应网格（*CONTROL_ADAPTIVE）通过在计算中动态细分高梯度区域，把计算资源集中到裂纹尖端或剪切带。关键参数包括ADPFREQ（细化检查频率）、ADPTOL（细化判据容差）和ADPASS。需要特别注意的是，自适应细化会引入历史变量插值误差，对损伤类材料可能造成损伤值跳变，因此建议只在弹性或小损伤阶段启用，或在细化后重置损伤变量。对于已经进入软化段的区域，继续细化反而可能加速删除、造成能量不守恒。

网格删除是模拟裂纹扩展最直接的手段，但必须谨慎。删除单元后，原本由该单元承担的载荷会瞬间转移，若周围单元未及时进入塑性，会出现“人工卸载波”，表现为力-位移曲线上的异常抖动。建议在*CONTROL_TERMINATION中设置DTMIN，避免时间步因小单元删除而崩溃；同时用*DATABASE_GLSTAT监控总能量，若沙漏能或删除能占比超过5%，说明删除策略需要调整。删除判据不宜只用单一应变阈值，可组合SIGVM（等效应力）、MAXPS（最大主应力）与MXEPS，并设置不同的失效顺序。

四、一个可复用的网格策略检查清单

最后强调一点：单元类型和网格策略不是独立选项，而是与材料模型、接触定义、加载速率共同构成一个系统。失效分析中，先确定物理机制（韧性、脆性、层裂、撕裂），再选单元与网格，最后用能量与收敛性验证，才能让仿真结果真正支撑工程判断。