Master/Slave面的不一致

这不是物理上的上下，而是数值积分中"探测点"和"目标面"的关系。通常，Slave面的高斯积分点探测Master面，检查是否有穿透。例如，Abaqus推荐将更粗的网格或刚度更高的材料面设置为Master。当刚度比超过100倍时，将硬的一侧设为Master可以改善收敛性。

计算过程更复杂。Slave点 $$ \mathbf{x}_s $$ 在Master面 $$ \Gamma_m $$ 上找到最近点 $$ \mathbf{x}_m $$，计算其法向距离（间隙 $$ g_N $$）和切向滑移。在惩罚法中，接触力为 $$ \mathbf{F}_c = \epsilon_N g_N \mathbf{n} $$（$$ \epsilon_N $$ 是惩罚刚度），当间隙为负（穿透）时会产生反力。

完全正确。大多数求解器（Ansys Mechanical、LS-DYNA）定义接触面的外向法线为正。只有当Slave点从Master面的"正"法线方向接近时，接触才有效。如果部件内向定义，接触会无法检测到，导致穿透。在网格生成时统一面的方向至关重要。

网格大小差异只是原因之一，本质是"搜索算法失败"。例如，在桶形搜索中，如果搜索距离（球形搜索半径）使用默认值，当网格尺寸差异很大时，可能找不到对应的单元。Ansys根据单元尺寸自动设置球形半径，但手动调整时，建议设为最大网格尺寸的2-3倍。

Slave点成为"孤立点"，不属于任何Master单元，接触约束无法生成。结果是该点的穿透约束不被添加到刚度矩阵中，导致部件看起来相互穿过。收敛判定会出现大的残差，导致分析停止或出现非物理的"接触突然打开"现象。

问题仍可能发生。数值上"完全一致"很少见，节点坐标通常有微小差异（1e-6 mm量级）。为处理这种"数值噪声"，接触设置有"调整容限"和"球形系数"等参数。COMSOL建议将"对的搜索容限"从"自动"改为"手动"，设置为模型总体尺寸的相对值（如1e-5）。

在创建接触面对前，进行以下预检查：1. 用前后处理软件可视化面法向，确保全部外向；2. 将接触预期区域网格尺寸比控制在3倍以内；3. 从接触面排除尖锐边缘和极度细长单元（宽高比20以上）。在Abaqus/CAE的"网格"模块中运行"验证网格"可以检测畸形单元。

改变接触定义方法。在Ansys Workbench中，"粘合"接触的"调整至接触"选项可自动修正初始间隙/过盈。或者，将细网格面（螺栓头）设为Slave，粗网格面设为Master，用"接触工具"明确指定球形搜索区域。经验上，初值可设为粗网格单元边长的1.5倍。

取决于情况。对于简单形状和少数面，手动选择更可靠。对于复杂装配体，自动接触（Ansys的"自动面对面"、Abaqus的"通用接触"）更高效，且优化了面方向和搜索算法，还能包含自接触。但计算成本会增加。使用自动接触后，应养成查看"接触工具"生成的接触对的习惯。

有很大区别。Abaqus/Standard非常严格，默认设置下不一致会输出大量警告消息（***WARNING: SLAVE NODE X IS OVERCLOSED BY Y），直接导致收敛不良。因此，初始间隙调整（"干涉配合"或"间隙"）几乎必须使用。LS-DYNA的自动接触则很宽松，有"软"选项可自动消除初始穿透。Ansys Mechanical处于中间，"程序控制"会自动尝试调整。

处理和修正方式不同。COMSOL Multiphysics的"对"型接触假定两个面法向"相对"。如果反向，有"反转"按钮简单修正。MSC Nastran的面接触（BSURF）由定义面的单元连接顺序（右手系）决定法向，需在网格生成阶段考虑。在HyperMesh中使用"定向单元"工具统一各部件法向是标准做法。

不匹配网格本身不是问题。问题在于自动网格生成器对各部件独立网格，接触边界的节点位置不一致。防止方法是使用"接触尺寸匹配"功能。在Ansys Workbench的"接触尺寸"中，对接触边设置相同的网格尺寸；或在Abaqus/CAE中用"按尺寸设置边"为接触边分配相同尺寸。但这只保证一阶单元对齐，二阶单元的中间节点位置仍可能有偏差。

首先可视化"初始接触状态"。在Abaqus中，.sta文件开头列出各接触对的初始开度/穿透量。在Ansys Mechanical中，通过"接触工具"→"初始信息"获取类似数据。如果出现意外的大穿透（例如1m模型有1mm以上），说明CAD干涉或部件位置错误。微小穿透（1e-3 mm以下）可用"调整至接触"或"界面处理"解决。

这正是"动态不一致"的典型症状。大变形或旋转导致初始Master/Slave搜索关系崩溃。对策有二：第一，使用"非线性自适应区域"（Ansys）或"*CONTACT_AUTOMATIC_SINGLE_SURFACE"（LS-DYNA）等在变形时重新定义接触面的功能。第二，动态扩大搜索半径。在Abaqus/Explicit中，将缩放因子从1.0（默认）增加到1.5或2.0，可提高变形后的接触捕捉概率。

几乎肯定是。这是"点接触"状态的证据。本应由面承受的荷载只由一个Slave节点承担。数值上，当惩罚刚度 $$ \epsilon_N $$ 过高（默认的"程序控制"不适当）时也会发生。首先，细化接触面网格至均匀。在Ansys中，将"法向刚度"从"程序控制"改为"手动"，设因子0.1试试。在Abaqus中调整"接触控制"的"刚度缩放因子"。

可以。用"对称面接触"或"基于多点约束的接触"。Ansys的"纯惩罚"算法有"对称"选项，自动反向Master/Slave计算。Abaqus/Standard的"接触对"可选"平衡Master-Slave"，或使用"通用接触"。最终方案是用"粘合"接触或"RBE3/MPC"近似约束全面，但这无法表示分离或滑动，应用范围有限。