延性破坏 — CAE术语解说

最大的区别是，破坏前材料是否会发生「塑性变形」。在延性破坏中，例如软钢（S45C）被拉伸时，会局部出现缩颈（Necking），断面积大幅减少。断裂面呈纤维状，肉眼可见变形。而脆性破坏（铸铁等）几乎无变形就突然断裂，断裂面呈粒状且平坦。

问得好。在单轴拉伸中，公称应力下降，但真应力却不断增加。然而，材料的「加工硬化」能力无法跟上局部断面积减少导致的应力集中，缩颈就由此开始。这被称为Considère条件。在真应力-真应变曲线上，满足

许多金属经历「空洞的发生、成长、合并」这个微观过程。在材料中的夹杂物（MnS等）或第二相粒子界面处，微小空洞（void）发生，在大的塑性应变下这些空洞成长，相邻空洞连接（合并），最终导致断裂。在电子显微镜下观察断裂面，会看到称为小窝（dimple）的特征图案。

单纯的塑性应变阈值是不够的。最典型的方法是基于「损伤力学」的模型。例如Lemaître模型或Gurson模型。Gurson模型将空洞体积分率

完全正确。这就是「应力三轴度」的影响。应力三轴度

「单元删除法」是早期的做法，但在动力分析中会产生问题因为应力会不连续地变为零。目前通常用「粘性正则化」或「弹性刚度降低法」。例如，Abaqus的延性损伤模型在达到断裂应变后，在指定的应变范围（如0.05）内逐渐衰减应力，使单元刚度平稳地趋于零。这样可以防止过度的收敛问题和不稳定的能量释放。

你说对了，这是最大的难点。最起码需要不同应力三轴度下的实验数据。实际中，结合光滑材拉伸试验（η≈0.33）、带槽圆棒试验（η≈0.6-1.0）、纯剪切试验（η≈0）等，收集各种条件下的断裂应变数据。例如，在汽车钢板成形限界图（FLD）制作中，会按照JIS Z 2254标准进行Nakazima试验，收集各应变路径下的断裂数据。

由于要追踪局部应变集中，网格依赖性是固有的问题。对策主要有两种。一是「根据网格尺寸调整参数」。例如，使断裂能量密度保持恒定，用网格特征长度

断裂位置吻合是必要条件但不充分。荷载-位移曲线的峰值后衰减行为、到达断裂时的能量吸收量，以及在多种不同试验形式（例如平面应变和平面应力状态）下的断裂行为都能正确重现，才能说明模型信度高。对汽车碰撞分析，需要在部件单体三点弯曲试验和实车碰撞试验两种情况下都保持一致性。

区别很大。Abaqus/Explicit内置了「延性损伤」模型，可以定义考虑应力三轴度和应变速率的断裂应变。还可以使用Gurson-Tvergaard-Needleman（GTN）模型，直接模拟空洞成长。而Ansys Mechanical（除APDL/LS-DYNA外）的标准塑性模型没有破坏判定，一般需要用用户子程序（UPF）自己实现或与LS-DYNA联合。

LS-DYNA的强项是破坏模型丰富且有大量应用实绩。特别是「GISSMO（Generalized Incremental Stress State dependent damage Model）」，同时考虑应力三轴度和Lode角，可以精确重现钢板的复杂破坏模式（如剪切破坏）。还有「MAT_ADD_EROSION」（单元删除条件），可以设置塑性应变、应力、时间等多种简单阈值，用户可根据目的和数据量选择合适的方法。

COMSOL Multiphysics的「结构力学模块」中有「损伤」接口可以处理延性损伤。默认实现了Lemaître损伤模型，可以定义损伤变量D的演化法则。用户还可以用偏微分方程（PDE）接口实现自己的损伤演化律，灵活性高。不过，与Abaqus和LS-DYNA相比，面向汽车碰撞的、预先用实验数据拟合的模型较少，需要用户自己根据理论和实验数据来构建模型。

主要有两个原因。第一，单元删除或刚度降低导致模型质量急剧减少或刚度矩阵病态。第二，破坏尖端出现极大应变速率，严格化了陆朗根据条件（CFL条件）。对策是调整「粘性正则化」参数使破坏平稳化，同时应用「质量缩放」来维持破坏区域以外单元的时间步长。

这是「损伤局部化带」集中在单个单元厚度、不向相邻单元传播的典型问题。非局部化模型是解决方案。如果使用的求解器没有这个功能，可以在用户子程序中实现一个过滤，即计算破坏参数（如应力三轴度）时不用单元中心值，而用周围数个单元的平均值。也可以选择保持网格相对粗糙，正确散逸破坏能量的「网格尺寸调整」方法。

首先确认输入参数是否「真的相同」。例如，断裂应变是真应变还是工程应变，应力三轴度的定义（η = p/q 或反过来等），单元积分方法的默认值（完全积分/减缩积分）导致的应力应变计算精度差异等。信任的准则是：用简单验证模型（如带槽拉伸试验）与实验结果最吻合的求解器。分析时一定要查阅求解器手册，确认所用模型的理论细节，比较时不能把求解器当黑盒。