退化单元错误

很好的问题。典型的情况是四面体单元的1条边被压扁，导致4个节点几乎在同一平面上。例如，理想四面体的体积计算公式为

这是自动网格划分算法的固有问题。特别是当用"四面体网格"填充复杂的CAD几何体时，在狭窄的间隙或尖锐边缘附近，单元会被压扁。此外，网格"扫掠"操作在扫掠路径急剧变化的地方也容易出现这种问题。Ansys Meshing的"Patch Conforming"算法对这个问题比较敏感，为了忠实追踪几何体的细微特征，在微小的面附近会生成极端扁平的单元。

因为数值积分的精度会致命性地降低。有限元法通过高斯积分点来计算单元内的应力和应变。退化的单元会导致形状函数的雅可比矩阵

两个阶段都会检查，但检查标准和严重程度不同。前处理网格划分（如Altair HyperMesh的"Quality Index"检查）会从形状比、歪斜度、最小角、雅可比数等多个角度检测"低品质"单元，并给出警告或尝试自动修复。而求解器（如MSC Nastran）在组装刚性矩阵时会尝试计算单元刚性矩阵，如果在数值上判定为奇异，就会输出致命错误（如FATAL MSG 3009）并停止求解。

这在求解器内部的单元库例程中进行。例如，对于等参数单元，在每个积分点计算雅可比矩阵行列式

有很大的区别。二次单元（有中间节点）由于形状函数是高次的，可以用弯曲的边来表示，但如果中间节点位置不当，很容易导致雅可比数为负的"反转单元"。相比之下，一次单元的雅可比数只取决于节点位置，检测相对简单，但一旦退化，积分点本身的值就会出错。因此，在非线性分析中，当出现大变形、单元被压扁时，二次单元往往会更早产生错误。

首先要确定错误信息中记载的单元ID。大多数求解器（如Ansys Mechanical的"Solution Information"）会输出导致错误的单元号。将该单元可视化，确认它是否靠近某个几何特征（尖角、极小的面、薄壁部分等）。这是根本原因分析的第一步。

有的，这被称为"几何理想化"或"特征抑制"。例如，Siemens NX Simcenter 3D或CATIA的CAE模块提供了"Virtual Topology"功能，可以将不必要的小边或小面并入相邻的大面。具体来说，直径小于1mm的凸台或加强筋、小于0.5mm的圆角通常对结构强度没有贡献，是可以抑制的候选。但是，应力集中位置的重要圆角应该保留。

有几种方法。第一，调整网格大小。在问题小面附近设置"局部网格大小"，将其设为略大于该面的尺寸（例如面为0.5mm，网格设为1.0mm），这样网格划分器可能会忽略这个面而直接通过。第二，改变网格类型。从四面体改为以六面体为主的网格。扫掠网格或映射网格更容易控制单元的纵横比，不容易产生退化。但对于复杂形状可能无法使用。

有业界经验法则。对于线性静力分析，最小雅可比比（相对于理想形状的比值）0.5以上还可以，0.7以上算是良好的。但对于非线性分析或接触分析需要更严格的标准，目标是0.8以上。关于歪斜度（Skewness），四面体的目标是0.9以下，六面体是0.8以下。这些值在NASA的"结构分析参考指南"和汽车业界的内部最佳实践指南中都有记载。

网格修复的理念和自动化程度上有明显差异。Ansys Workbench的"Mesh Metrics"可以详细展示网格质量指标，并提供"Patch Independent"网格方法，这种方法可以在一定程度上忽略小的几何特征。Abaqus/CAE的网格划分器对几何体忠实度高，但"Mesh Edit"工具包有强大的功能，可以手动移动或分割单个单元，适合熟练用户进行精准修复。COMSOL Multiphysics在"自由四面体"网格中内置了"避免极度细长单元"的算法，可以事先进行一定程度的预防，但一旦发生问题，诊断信息较为简单。

正是如此。开源网格划分器（如Netgen、Gmsh）在算法上很优秀，但缺乏商用软件那样鲁棒的"几何修复"功能。因此，读入稍微"脏"一点的STEP或IGES文件时，容易直接生成退化单元。另外，CalculiX等求解器的错误信息很简洁，只输出"negative jacobian"，要找出问题单元比较费时。商用软件的优势在于对"肮脏现实"数据的前处理耐性和丰富的诊断信息。

有两个原因。第一，它们以"网格质量管理"为最优先，在UI和功能上都围绕这个中心设计。例如，ANSA的"Morphing"工具箱可以直接拖动网格节点修改形状，同时实时监控品质指标（雅可比数、纵横比等）的变化。第二，其中融入了行业特有的最佳实践。面向汽车业的"半自动中面网格"功能能高效处理薄板零件，将其转换为中面，内置的算法可以预防薄板退化单元。

原因相似但数值状态不同。"Negative Jacobian"表示单元节点顺序错误，行列式值实际为负。这在大变形分析中单元被反转或网格生成时节点编号有问题时发生。"Zero or Negative Jacobian"除了为负外，还包括接近零的情况。后者更普遍，单纯单元被压扁时也会出现。无论如何，求解器都判断无法正确计算该单元的刚性矩阵而输出错误。

有最后的手段，但结果的可信度会降低，需要谨慎。第一，改变求解器设置。Abaqus/Explicit可以改变默认的"Element deletion"设置，让被压扁的单元被删除后继续求解。第二，改变单元类型。如果Abaqus/Standard中C3D4（线性四面体）不行，改用C3D4H（线性四面体、混合压力公式）有时能缓解体积锁定，增强对形状歪斜的容错性。第三，极端粗化网格。选择问题局部区域，将网格大小改为原来的2倍或3倍，使小特征被网格跨越。

这是由于"大变形导致单元破坏"。初始状态形状比良好的单元在求解过程中发生大变形或屈曲，最终被物理上压扁了。这在橡胶材料的超弹性分析或金属大塑性变形分析中特别常见。对策有两个。一是使用网格自适应（重新网格划分）功能。Abaqus的"ALE adaptive meshing"或LS-DYNA的"Adaptive Remeshing"可以在求解中自动为变形大的区域重新生成网格。二是有意地将问题区域的初始网格设得更细，分散因变形引起的形状歪斜。

可能是间接原因。特别是用"罚函数法"进行接触时，接触面的法向方向会插入一个极大的虚拟弹簧。如果接触面网格不当（例如从属面网格太粗），接触反力会集中在个别节点上，这些反力会使周围的单元在数值上扭曲，恶化雅可比数。对策是使接触面网格大小均匀，从属面应比主面细（经验法则是细1.2至2倍）。另外，Ansys中可以试着从罚函数法切换到"Normal Lagrange"法，这样可以在某些情况下缓解罚函数法引起的过度反力集中。