FEA仿真分析 — 螺栓连接非线性接触的收敛调试实战

为什么一个”简单”的螺栓连接分析跑了三天还不收敛

在FEA仿真分析中，螺栓连接几乎出现在每一个装配体分析中。表面上看起来很简单——一个螺栓穿过两个零件的孔，拧紧螺母产生预紧力。但在有限元模型中，这涉及到接触非线性（螺栓头与法兰面、螺母与法兰面、螺纹与孔壁、法兰面之间——至少四处接触对）、几何非线性（如果有大变形）和材料非线性（如果在螺栓中发生屈服）。这三重非线性叠加使得收敛问题成为螺栓分析的头号障碍。

我们的案例是一个法兰连接——6组M20螺栓将两个DN200法兰连接在一起，螺栓预紧力60 kN，法兰承受内压和弯矩。初始模型用了Abaqus Standard（隐式求解器），可是计算一直在0.1秒附近反复迭代不收敛——耗时三天仍没有走出第一步。

本文复盘这个看似简单却处处陷阱的FEA仿真分析案例，重点分享接触设置策略和收敛控制的实战经验。

接触对设置：你以为的接触可能根本不是接触

螺栓连接中接触对的数量和类型直接影响收敛性。对于M20螺栓，螺栓头下表面与法兰顶面的接触、螺纹区域与法兰孔壁的接触——这些是显而易见的。但法兰面之间的接触——两个法兰端面在预紧力未施加时可能处于分离状态，在预紧力施加后发生闭合——这涉及”初始间隙”的处理，是收敛问题的重灾区。

我们的模型在使用默认的”Adjust only to remove overclosure”选项时，由于法兰端面的初始间隙（约0.02毫米），接触对在主分析步开始时是未闭合的。在预紧力施加过程中，接触突然从”开”变为”闭”——这种状态跳跃在隐式求解器中极易导致收敛失败。

解决方法是在接触对定义中使用”Interference fit”选项或施加一个小的”稳定化”阻尼。我们选择了后者——添加了一个与接触间隙比例相关的粘性阻尼项，在接触闭合过程中自动衰减。这个技术相当于给接触状态的变化添加了一个”缓冲”——避免跳跃式变化，使牛顿-拉弗森迭代更平滑。参数设置为默认阻尼系数的1%，足以稳定收敛而不影响最终结果的精度。

预紧力施加：Bolt Load的”正确”使用方式

Abaqus中的Bolt Load功能可以在螺栓截面上施加轴向力或缩短量来模拟预紧效果。使用Bolt Load时有三个关键步骤和两个常见错误。

三个关键步骤：第一步，在螺栓上创建一个截面（Surface），Bolt Load作用于此截面。截面必须选择螺栓杆的横截面（垂直于螺栓轴线），且必须在螺栓杆处于自由状态的位置（不要在螺纹或螺栓头过渡区创建截面）。第二步，在初始分析步中施加一个非常小的预紧力（如10 N）来”固定”接触状态——这相当于给接触一个初始值，帮助后续的牛顿迭代找到正确的解方向。第三步，在后续分析步中将预紧力增加到目标值（60 kN），然后”Fix at current length”——将螺栓的伸长量冻结，后续载荷（内压、弯矩）作用时螺栓保持定长而非定力（后者更符合物理实际）。

两个常见错误：错误一，跳过”小预紧力”步骤直接施加目标预紧力。在接触尚未建立的情况下施加60 kN的轴向力，螺栓头和螺母会”穿透”法兰面——导致极端的单元畸变和立即的不收敛。错误二，在预紧力施加完成后没有”Fix at current length”——如果保持定力状态，后续的内部压力和弯矩载荷会导致螺栓载荷异常增加或减小，与实际物理行为不符。

收敛控制：自动稳定化与人工干预的平衡

当常规的牛顿-拉弗森迭代不收敛时，Abaqus提供两种自动干预手段：自动时间步缩减（Automatic Stabilization）和接触稳定化阻尼。关键是不滥用它们。

自动时间步缩减通过将当前增量步一分为二（甚至更小）来挽救不收敛的迭代。这在接触状态突变的初期是有帮助的——比如螺栓头和法兰面发生第一次接触时。但过度依赖自动时间步缩减会导致增量步大小不断萎缩——从0.01秒减到10⁻⁵秒，分析时间呈指数增长。有一个实用技巧：设置最小增量步大小为初始值的10⁻⁴，一旦触发这个下限就停止分析并检查模型——而不是让分析无限地”挣扎”在极度细碎的时间步中。

人工稳定化（Damping Factor）是另一个双刃剑。在接触分析中添加一个小数值的人工阻尼（如1×10⁻⁷），阻尼力与节点速度成正比，能有效抑制高频振动改善收敛。但阻尼太大（>1×10⁻⁴）会导致非物理的能量耗散——阻尼力做的功与外力、应变能相比较变得显著。一个验证阻尼合理性的方法是比较ALLSD（稳定化耗散能）与ALLIE（总应变能）的比值——如果ALLSD/ALLIE > 5%，需要减小阻尼系数。在我们的案例中，ALLSD/ALLIE稳定在约1.5%，说明阻尼设置合理且对结果精度影响可忽略。

结果解读：不要只看应力云图

螺栓FEA仿真分析完成后，很多工程师只看Von Mises应力云图——”螺栓最大应力不超过屈服强度，结构安全”——然后就关闭文件了。这是严重的不足。螺栓连接的关键性能指标不是螺栓本体的应力，而是连接面之间的接触压力分布和可能的分离状态。

在我们的法兰案例中，Von Mises应力一切正常——螺栓最大应力约550 MPa，低于8.8级螺栓的屈服强度640 MPa。但检查接触状态后发现，在弯矩最大的工况下，法兰面在远离螺栓的位置发生了约0.03毫米的分离——虽然这个数值极小，但对于密封性能要求严格的法兰连接来说，微小的面分离意味着潜在的泄漏风险。

另一个容易被忽视的结果是螺栓的残余夹紧力。在加载和卸载后，螺栓中的力是否恢复到初始值？如果存在显著的夹紧力损失（超过初始值的10%），说明发生了局部塑性变形或接触面嵌入——这会影响连接的长期可靠性和疲劳寿命。在FEA仿真分析中，螺栓连接的”隐性失败”往往藏在接触压力和残余夹紧力这些不那么直观的结果量中。

——————————————————————————————

📐 科研计算代算服务

我们提供专业的第一性原理计算、分子动力学模拟、有限元仿真分析及材料表征测试服务。
团队由博士学历的计算化学、计算材料学专家组成，覆盖VASP、LAMMPS、GROMACS、ANSYS、COMSOL等主流工具链。

  🔗 立即咨询计算方案