歪曲单元错误

这是个很好的问题。单元的歪曲（Distortion）是一个指标，用于定量评估单元从理想形状（比如正方形或正四面体）的变形程度。具体来说，它通过单元顶点坐标计算出的雅可比矩阵的行列式、宽高比、内角等多种几何参数来评估。例如，四面体单元歪曲度的一个指标定义如下。

本质上有两个原因。第一，数值积分精度下降。有限元分析在高斯积分点处评估应力和应变，但歪曲单元会导致积分点的物理位置偏离，积分精度大幅下降。第二，形状函数导数计算变得不稳定。当雅可比矩阵接近奇异时，其逆矩阵计算会产生数值误差爆炸，刚度矩阵无法正确计算。例如，如果某个四面体单元的雅可比矩阵行列式低于1.0e-10，即使在双精度计算中，有效数字也会几乎完全丧失，求解器会发散。

确实有差异。每个求解器都有内部阈值。例如，Abaqus/Standard对四面体和六面体单元执行"形状检查"，默认对歪曲容限设置为0.1~1.0范围，发出警告或错误。而Ansys Mechanical广泛使用"Element Quality"指标，在0.0~1.0范围内评估，其中1为理想（正四面体或直方体），一般来说0.1以下被判定为"Poor"。这些阈值在线性静力学分析和非线性、动力学分析之间是不同的。动力学分析要求更严格的标准。

基本上是对所有单元的检查。预处理器（比如Ansys Meshing或Abaqus/CAE）在网格生成时或导入时以批处理方式执行。具体来说，对每个单元类型（四面体、六面体等）计算定义的"品质指标"。对于四面体，计算体积比、边长比（宽高比）、内角的最大值和最小值等。六面体单元更复杂，例如将"歪曲参数"

Jacobian Ratio是在单元内多个位置（通常是积分点或顶点）计算的雅可比矩阵行列式的最小值与最大值的比。理想的直方体单元在所有位置的雅可比矩阵行列式相同，比值为1。当单元歪曲时，各位置的值会分散，比值会远离1。Ansys中，在评估二阶单元（有中间节点）时特别重要。例如，四面体二阶单元（SOLID187）在顶点、边的中点、体积中心等多个采样点计算雅可比矩阵行列式。允许值取决于分析类型，但在非线性分析或涉及接触的分析中，当比值超过30时会出现警告，超过100时出现错误。

二阶单元的敏感性要高得多。原因在于形状函数。一阶单元的形状函数是线性的，所以单元内的变形梯度是常数。但二阶单元的形状函数是二次的，中间节点的位置极大地影响单元内的变形映射。当中间节点远离边的中点时，单元内部的映射变得高度非线性，雅可比矩阵行列式在各处会极端变化。因此，Ansys和Abaqus在使用二阶单元时，默认启用"中间节点放在边中点"选项（Abaqus中的"Midside node parameter"）。禁用这个选项会导致即使是微小的几何凹凸也会引发大量歪曲错误。

首先，确定错误信息中列出的单元ID，然后可视化该单元在模型中的物理位置。在Ansys中使用"Select > By Element Number"，在Abaqus/CAE中使用"Query > Element"可以找到它们。然后思考"为什么在这个位置生成了这个单元"。尖角、非常细长的面、两个部件靠得很近的接触区域等是典型的发生位置。

修改几何形状是首选。在CAD中给那个尖锐肋角添加0.5mm或1mm的倒角（圆角）。实际产品上往往也存在倒角。如果无法修改CAD，就用网格控制来处理。Ansys Meshing中可以使用"Face Meshing"或"Edge Sizing"来沿肋规则排列单元。或者改变网格类型，试试能否将那个区域从四面体（Tet）改成六面体（Hex）扫掠网格。万一都不行，就整体细化网格。虽然会增加计算成本，但歪曲会缓解。

不一定需要全部修复，但需要战略性判断。首先确认这些低品质单元是否位于"预期高应力集中区域"。螺栓孔边缘或缺口根部等关键结果区域的低品质单元必须修复，即使只是警告等级。但在约束面或远处、应力几乎不作用的区域，可以忽略。不过，非线性分析（材料非线性、接触）或特征值分析（模态分析）中，低品质单元更容易引起数值不稳定，需要更严格处理。

情况而定。网格平滑化通过移动单元内部节点来改善形状，Abaqus/CAE和HyperMesh中有实现。但边界节点不动，所以根本的形状问题解决不了。过度平滑化还会导致单元互相穿刺或几何形状严重偏离原始形状。实务中，先进行局部重新网格划分或网格大小调整，实在不行再轻度平滑化作为最后手段。在Ansys Fluent等CFD求解器中网格平滑化很常见，但在有限元分析中要谨慎使用。

用户界面和默认严格程度有区别。Ansys Mechanical作为"Mesh Metric"提供"Element Quality""Skewness""Orthogonal Quality"等多角度指标，并用直方图显示分布，这是一大特色。Abaqus/CAE在"Mesh"模块的"Verify"选项卡执行"Analysis checks"，用不同颜色显示单元的"Distortion"程度。COMSOL Multiphysics在"Mesh"节点的"Statistics"中显示"Element quality"，用0~1值表示。默认警告阈值方面，Abaqus相对保守（严格），Ansys标准，COMSOL更多由用户自己设置。

差别很大。Abaqus/Standard检测到歪曲单元时，默认停止分析并输出错误信息。但如果在输入文件中添加关键字"*DIAGNOSTICS, DISTORTION CONTROL=YES"，会启动"歪曲控制"功能，自动删除或修复歪曲单元继续计算。Ansys Mechanical中，歪曲单元严重时，求解器在生成刚度矩阵时失败，通常以"A general failure occurred in the solver"这样的笼统错误停下，要花时间定位具体单元。COMSOL对"歪曲网格单元"的自动重新网格划分功能很强，易于与自适应网格细分结合。

开源求解器的预处理器与求解器联系不如商业软件紧密，用户需要自己承担网格品质管理的责任。比如CalculiX（预处理器如CGX等）在生成网格时提供品质检查命令，但标准相对简单。Code_Aster与Salome-Meca环境联动，Salome网格模块提供的品质指标（NICO、Quality）需要用户自己理解。总的来说，开源软件的错误信息比商业软件简洁，不是"歪曲单元"而是更直接的表述如"零或负体积单元"这样的错误。

几乎肯定是CAD数据的"几何缺陷"导致的。STL文件是面元数据，面与面之间可能不完全相连，存在"间隙"或"重叠"。STEP/IGES文件中也可能含有极短的边（长度0.001mm等）、微小面、表面间微小间隙。网格划分器忠实地沿着这些不完美的几何形状，生成极端细长的三角形或四边形，结果导致歪曲单元。对策是使用预处理器的"几何修复"功能（Ansys的"SCDM"、Abaqus的"Geometry Edit"）删除微小特征或缝合间隙。

这是"几何非线性"导致的典型问题。当材料大幅变形时，单元自身也被拉伸、压缩、严重歪曲。初始状态良好，但变形后雅可比矩阵行列式可能变负（单元翻转）。对策有多个。第一，审视单元公式。Abaqus的混合单元（比如C3D8H）对体积变化大的材料有效。第二，细化网格，用更小的单元追踪变形。第三，改变求解器设置。Ansys中启用大变形选项（NLGEOM,ON），进一步利用"Mesh Morphing"或"Rezoning"（Abaqus中的"Adaptive Meshing"）在变形过程中重新组织网格。这些方法计算成本高，但对继续分析是必要的。

因为接触计算直接依赖单元形状，要计算表面法线和侵入量。歪曲单元的表面法线方向急剧变化，接触判定的积分点配置不均，导致接触力计算不稳定。结果可能无法收敛，或出现异常大的接触力，进而加重单元歪曲，形成恶性循环。具体对策：(1) 接触面尽量用六面体单元或高品质四面体；(2) 接触面网格大小保持一致（大小比超过10:1容易出问题）；(3) 初始接触设置为"Adjust to Touch"等消除初始贯入；(4) 尝试把接触算法从"Augmented Lagrange"改为"Pure Penalty"（情况而定）。

检查"歪曲"以外的网格品质指标，特别是"宽高比"和"偏斜角"。细长的单元（宽高比大）会导致一个方向的变形刚度被过度高估，应力偏离真值。"偏斜角"大的单元（边与垂直方向偏离大）会在剪切变形计算中引入误差。Ansys的"Mesh Metric"可以批量检查这些。另一个重点是结果评估方法。单元内应力在高斯积分点计算，然后外推到节点。歪曲大的单元中，这种外推本身会放大误差。因此，歪曲区域的"单元解"不能信任，即使看起来"节点解（平均化）"很光滑，那也可能是误导性的平均化结果。根本解决只有改善该区域的网格品质。