Abaqus网格大变形之CEL（欧拉法）

模型穿透在Abaqus仿真中非常常见且关键的问题。模型穿透通常发生在接触分析中，具体表现为一个物体（从面）的节点穿透了另一个物体（主面）的几何表面。

这不仅仅是视觉效果问题，它意味着你的接触定义失效了，会导致计算结果严重失真、不收敛，甚至得到完全错误的结果。

以下是导致模型穿透的主要原因以及相应的解决方案，我将从最常见到较少见的顺序进行说明。

一、 接触定义设置问题（最主要的原因）

这是绝大多数穿透问题的根源。

1. 接触属性定义不当

 · 法向行为： 默认的“硬”接触通常是最好的，但要确保其生效。

 · 切向行为： 如果摩擦系数设置过小或不设置摩擦，物体更容易发生滑动和穿透。

 · 解决方案： 检查接触属性，确保“硬”接触被使用，并根据实际情况设置合理的摩擦系数。

2. 接触对调整不足

 · 初始过盈： 如果模型在分析一开始就存在穿透（例如，装配时两个零件重叠了），Abaqus在第一步计算中会尝试解决它，但如果过盈量太大，可能会导致求解失败或无法完全修正穿透。

 · 解决方案：

  · 在Interaction模块中： 定义接触时，在Edit Contact Pair对话框中选择Adjust only to remove overclosure。这会在分析开始前自动调整节点位置以消除初始穿透。

  · 在Load模块中： 使用BC或Predefined Field中的Adjust for contact功能。

  · 修改几何模型： 从根本上修正CAD模型，确保没有初始重叠。

3. 接触主从面选择错误

 · 原则： 从面应该是网格更细、材料更软或几何更复杂的表面。

 · 错误示例： 如果将网格粗糙的面设为从面，其节点可能会在网格精细的主面单元之间“掉进去”，造成穿透。

 · 解决方案： 严格遵守主从面选择原则。当两个接触面网格密度相差较大时，将网格更细的面作为从面。

二、 网格问题

网格质量直接影响接触计算的精度。

1. 网格过于粗糙

 · 原因： 粗糙的网格无法精确描述接触表面的几何曲率，使得接触点计算不准确，节点很容易从一个单元“跳”到另一个单元后面，造成穿透。

 · 解决方案： 在接触区域进行网格细化。这是解决穿透问题最有效的方法之一。

2. 单元类型不当

 · 原因： 对于接触问题，使用线性减缩积分单元（如C3D8R）有时会因为沙漏模式或本身的计算特点而导致接触力计算不精确。

 · 解决方案： 在接触区域考虑使用线性完全积分单元（如C3D8） 或二次单元（如C3D20）。二次单元因为中间节点的存在，能更好地描述曲面，接触精度更高。

三、 求解控制与步长问题

1. 时间增量步过大

 · 原因： Abaqus/Explicit中，过大的时间步长可能导致在一个增量步内，物体移动的距离超过了接触单元的尺寸，从而“跳过”了接触检测，直接发生穿透。在Abaqus/Standard中，过大的增量步会导致收敛困难，也可能在迭代过程中出现穿透。

 · 解决方案：

  · Explicit： 使用质量缩放 来稳定地增大稳定时间极限，但需谨慎使用，避免引入过多的人工惯性效应。让Abaqus自动控制时间步长通常是最佳选择。

  · Standard： 在Step模块中，减小初始增量步大小和最小增量步大小，并增加最大增量步数，给求解器更多机会去收敛。

2. 接触约束算法问题

 · 原因： Abaqus/Standard默认的接触算法是“拉格朗日乘子法”或“罚函数法”的变体。有时默认设置对某些剧烈变化的接触不够“强硬”。

 · 解决方案： 在Interaction模块中，编辑接触对，进入Contact Controls。可以尝试创建更严格的接触控制，例如：

  · 将初始接触和严重不连续迭代的许可度收紧。

  · 对于大滑移问题，确保设置了正确的滑移公式。

四、 模型本身的问题

1. 刚体设置

 · 原因： 如果一个物体被定义为刚体，其表面不会发生变形。如果另一个可变形体与之接触，在计算不收敛或接触定义不当时，可变形体的节点可能会穿透刚体表面。

 · 解决方案： 确保刚体表面的网格足够精细，并且正确设置了主从面。

总结与排查步骤

当遇到穿透问题时，建议按以下顺序排查：

1. 可视化检查： 首先在Visualization模块中，打开Symbols和Deformed Shape Plot，将变形比例设置为1。查看穿透是发生在分析过程中，还是在初始状态就已存在。

 · 初始穿透： 重点检查接触对调整和几何模型。

 · 分析中穿透： 继续以下步骤。

2. 检查接触定义：

 · 确认主从面选择是否正确。

 · 确认接触属性（硬接触、摩擦系数）是否合理。

3. 细化网格： 这是解决大多数动态穿透问题的首选方法。显著细化接触区域的网格，然后重新计算。

4. 调整求解参数：

 · Explicit： 检查时间增量步是否合理，考虑质量缩放。

 · Standard： 减小时间增量步，放宽收敛准则（如果可能），或使用自动稳定。

5. 检查单元类型： 在关键接触区域尝试使用更高阶的单元。

6. 使用接触诊断工具： 在Interaction模块中，使用Tool -> Query下的接触诊断工具，检查初始状态下的接触间隙和过盈量。

一句话总结： 模型穿透的本质是接触力未能及时、足够地阻止物体的相对运动。解决思路就是优化接触定义、细化接触网格、控制求解步长，三者结合，缺一不可。

免责声明：本文系网络转载或改编，未找到原创作者，版权归原作者所有。如涉及版权，请联系删