孤立节点错误
孤立节点
老师，"Orphan nodes detected"是什么意思？
孤立节点错误的理论基础
孤立节点的物理意义
"孤立节点"是网格的一个常见错误，但从物理角度讲，问题到底在哪里呢？只不过是个孤零零的点，为什么分析会停止呢？
问得好。从物理角度讲，在有限元法中，"节点"是承载结构物质量和刚度的地方。孤立节点不属于任何单元，因此质量和刚度都为零，变成了一个"幽灵节点"。例如，如果你在这个点设置强制位移，由于质量为零，运动方程
质量和刚度为零意味着那个点"不存在"一样。但为什么会产生这样的点呢？即使是从CAD中导入的干净模型也会发生，原因是什么？
主要原因是"几何学不一致"。例如，两个表面"看起来"相接，但实际坐标值有微妙的偏差。CAD的合并容差（例如CATIA V5的默认值0.001毫米）小于这种隙度时，网格生成器会将其识别为不同的边。网格生成时，一个表面的边上生成节点，而另一个表面的单元可能不包含该节点，这样就留下了孤立节点。
所以微米级的隙度是原因。但那个"孤立节点"在坐标上是靠近周围节点的，不能通过自动"拼接"功能将其与周围节点合并吗？
"拼接"或"接触"设置是用来定义不同零件之间的"接触"的，不能解决网格生成前的几何问题。Abaqus/CAE的"Merge/Cut"功能或Ansys SpaceClaim的"Share Topology"是将几何一体化的操作。但是，1e-5毫米以下的极小偏差检测不到，会在网格生成阶段显现出来。这就是"看似清理过的模型仍出现错误"的原因。
孤立节点错误的数值计算方法
求解器检测机制
求解器具体用什么算法来发现孤立节点？是逐个检查每个节点吗？
大致上是通过"连接矩阵"的创建和检查。构建表示单元-节点连接关系的矩阵后，搜索任何单元都不引用的节点索引。计算复杂度为O(N)（N是节点数），但对于大规模模型不能忽视。MSC Nastran等求解器在刚度矩阵
不属于任何单元的节点在连接矩阵中行（或列）全为零。但是，如果有集中质量单元（比如MASS21）附加在孤立节点上，不会被判为错误吗？
敏锐的观察。你说得对，Ansys的MASS21或Abaqus的*MASS单元是"零维单元"，它单独就能"激活"一个节点。因此求解器不会将其视为孤立节点。但这仅限于结构分析。在热分析中，如果该节点只有温度自由度，未连接到热传导单元（例如SOLID70），仍会出现错误。自由度的一致性是关键。
自由度的一致性……也就是说，求解器还在管理"每个节点的活跃自由度"和"支撑这些自由度的单元"的映射？这些信息在哪里查看？
通常对用户是隐藏的内部数据结构。但启用调试输出就能窥探。例如，在LS-DYNA中设置`memory=2`关键字，会输出详细的内存分配信息和节点-单元连接性。Ansys的`/DEBUG`命令也一样。但是输出文件会很庞大，实际上应该在问题节点周围的子模型上试验。
孤立节点错误的实际应用
发现到修正的工作流
实际上遇到"孤立节点错误"时，首先应该做什么？仅从错误信息无法判断节点的位置。
首先从错误消息中确定节点号。Ansys会显示"Node X is not connected to any element."，Abaqus会显示"Node set
找到位置后，直接删除该节点不就行了吗？但我担心不解决为什么生成那个节点的根本原因，可能其他地方还会出现，感到不安。
你说得对。删除只是对症疗法。根本原因几乎100%在"几何"中。具体步骤是：1) 放大问题节点周围的几何（边、顶点），详细查看。2) 检查是否有重复的边或微小的面隙度。3) 回到CAD软件（SolidWorks、NX、Creo），重新定义相关部分的草图或特征，或者将"合并"容差放宽至约0.01毫米，重新导出（STEP或IGES）。
很难回到CAD的情况下，能在CAE前处理软件内修复吗？
可以。使用Ansys SpaceClaim或Abaqus/CAE的"Geometry Edit"模块，或Altair HyperMesh的"Geometry Cleanup"工具。具体来说，可以应用"Stitch"（缝合）、"Merge"（顶点合并）、"Fill Hole"（填洞）功能。但这些自动修复可能改变形状，所以应该在重要特征（如肋、凸台）远处应用，修复后还要进行尺寸检查。
在重新生成网格前，有办法预测不会产生孤立节点吗？
生成网格后、执行分析前，进行"网格质量检查"。大多数前处理器都有"查找未连接节点"或"自由边"检测功能。在HyperMesh中，打开"Tool"面板的"Check Elems"，选中"Detached"进行检查。将这作为日常工作流程的一部分，就能在执行求解器前发现并修正。
孤立节点错误的软件比较
各软件的处理方法差异
在Ansys、Abaqus、COMSOL中，孤立节点错误的处理方法有区别吗？
有显著差异。Ansys Workbench的"Geometry"单元（SpaceClaim）和"Model"单元（Mechanical）是分离的，几何修复和网格管理是独立的。孤立节点主要在"Model"单元发现，解决原因需要回到"Geometry"单元。而Abaqus/CAE是单一环境，从几何编辑到网格、分析设置都在一起，可以通过"Merge/Cut"操作直接修复。COMSOL Multiphysics有"Form Assembly"功能将组件结合，共享网格的结合界面设计，基本不会产生孤立节点。
在Ansys Mechanical内部有直接删除孤立节点的命令吗？
GUI中没有。但用APDL（Ansys Parametric Design Language）命令可以。如果孤立节点号是100，用`NDELE, 100`命令删除。但这是网格编辑操作，重新生成网格时可能复现。更有效的方法是进入`/PREP7`，用`NUMMRG, NODE`命令（节点合并）将邻近节点统合。容差用`VALUE`选项指定（例：`NUMMRG, NODE, , , 1E-4`）。
免费和开源求解器（CalculiX、Code_Aster）怎么处理？
CalculiX（CCX求解器）读取输入文件（.inp），在`.dat`文件中输出错误。孤立节点通常表现为"ERROR in e_c3d: nonpositive jacobian in element..."这样的间接错误，原因很难追踪。用Salome-Meca作前处理器时，可以在"Mesh"模块用"Clean mesh"算法执行"Remove detached nodes"。Code_Aster在`.mess`文件中输出详细错误日志，节点连接性检查在`DEFI_MODELE`阶段进行。
孤立节点错误的故障排除
常见陷阱和对策
几何修复好了，重新生成网格，又在另一个节点号出现孤立节点错误。陷入无限循环。怎么办？
这是几何拓扑有根本问题的表现。例如，用实体网格化薄壁结构时，如果面由非常细长的三角形组成，边上的节点可能被邻接面的单元排除，导致孤立。对策有两个：1) 调整网格尺寸。问题区域的局部网格尺寸设为最小边长的1/2以下（例：0.1毫米的边用0.05毫米的网格）。2) 改变网格类型。从四面体（Tet）改为六面体（Hex）支配网格，规则性网格不容易产生孤立点。
在用"接触"设置的多体模型中，接触面出现孤立节点。原因是什么？
接触设置本身不是原因。接触对象的面（Target面和Contact面）属于不同的体时，网格分别独立生成。此时，两者的网格密度相差很大，密集方的网格节点可能位于粗疏方网格单元的"外侧"，导致孤立。防止办法是用Ansys的"Contact Sizing"或Abaqus的"Mesh Constraints"，使接触面间的网格尺寸相似。
网格成功生成，求解器也开始运行，但中途因"零枢轴"或"过度变形"停止。这也可能与孤立节点有关吗？
很有可能。特别是非线性分析（接触、大变形）中，计算过程中单元严重扭曲，雅可比行列式变负。结果是单元和节点的连接性在求解器内部丢失，产生"事实上的孤立节点"，刚度矩阵变奇异。这是"零枢轴"的一个原因。对策是严格检查网格质量（纵横比、扭曲度），在预期大变形的区域应用自适应网格细分（Abaqus的ALE、LS-DYNA的Adaptive Meshing）。
最后，有办法彻底防止孤立节点错误吗？有"黄金法则"吗？
没有"黄金法则"，但有确实可行的最佳实践。1) **CAD阶段**：执行装配的"干涉检查"，零件间隙除非有意否则为0毫米。2) **导入时**：导入CAE软件时，STEP AP214（包含形状和颜色）不如STEP AP203（仅形状）可靠。3) **网格前**：在前处理器中"体积检查"确认为闭合实体。4) **网格后**：分析前必须执行"自由边"和"孤立节点"检查。养成这4步的习惯。
更详细
错误