大家好，我是仿真秀专栏作者—力学AI有限元，985硕士毕业，博士在读，现从事力学仿真研究，具体研究方向为疲劳寿命分析，有6年有限元仿真经验，擅长ABAQUS模拟仿真，isight参数优化。fe-safe疲劳寿命计算和abaqus的二次开发等。当前，在仿真秀官网和APP制作并发布了10多套力学与有限元仿真视频教程，旨在帮助理工科学子和工程师能够学以致用，提升科研技术创新能力。2020年7月17日15时41分许，河北某风电场突发重大设备故障，该风场11号2000kW风机（2011年1月投产，2014年11月出质保，由风电场自主维护）在正常运行中通讯突然丢失，1分钟后附近作业人员听到倒塔声响，确认风机已倒塔报废，造成严重经济损失。事故发生时，现场小雨、风速11.2m/s，风机转速21.21rpm、负荷1827kW，各项运行参数均符合标准；事后调取该风机前7天运行数据及30天紧急变桨测试记录，未发现异常迹象。 为查明事故根源，技术团队立即开展全面失效分析，金属专家现场检查发现，128颗塔筒连接螺栓中有12颗存在疲劳失效，存在先断裂、承载力转移的情况，初步判定事故直接原因是第一节与第二节塔筒连接螺栓断裂。随后通过宏观观察、化学成分分析、扫描电镜分析、金相检验及力学性能测试等一系列专业手段，最终锁定了螺栓断裂的核心诱因。宏观观察显示，断裂部位集中在塔筒法兰连接螺栓的螺纹牙底处，断口存在明显的“海滩纹”特征，这是典型的疲劳断裂标志；扫描电镜分析进一步证实，螺纹牙底存在多处裂纹萌生源，且裂纹扩展路径上分布着加工刀痕；金相检验发现螺栓内部存在硫化物夹杂物，这些夹杂物成为应力集中的薄弱点；力学性能测试则表明，螺栓表面硬度因增碳处理出现超限情况，进一步降低了其抗疲劳能力。综合分析可知，此次螺栓断裂是多因素叠加的结果：首先，螺纹牙底的几何结构本身存在应力集中效应，在风机运行过程中承受的多方向剪切应力及拉应力作用下，极易形成裂纹源；其次，加工过程中残留的刀痕与材料内部的硫化物夹杂物，极大地促进了裂纹的形成与早期扩展；最后，表面增碳导致的硬度超限，使得螺栓韧性下降，裂纹扩展速率加快，最终在服役一段时间后发生断裂失效，进而引发风机倒塔事故。螺栓作为机械装备中不可或缺的核心连接部件，其失效往往会引发连锁反应，甚至造成重大安全事故，因此精准识别螺栓失效原因、提前规避失效风险至关重要。有限元分析（FEA）作为一种基于计算机信息化处理的“虚拟实验”技术，已成为螺栓连接设计优化与失效排查的关键手段，其核心价值体现在三个核心场景：一是结构强度评估，可精准分析螺栓在不同负载条件下的承载能力，避免局部过载现象；二是失效模式识别，能模拟螺栓在复杂工况下的力学响应，定位裂纹萌生源头与扩展路径；三是设计优化迭代，基于模拟结果可对螺栓的材料选择、几何结构、紧固方式等进行优化，在提升性能的同时降低成本。我们知道，螺栓断裂失效并非单一因素导致，而是材料缺陷、设计不足、安装不当与环境作用等多因素叠加的结果。从近期失效案例来看，主要诱因可归纳为四类：- 一是应力集中引发的疲劳断裂，如风电螺栓过渡圆角处因几何突变导致应力集中，在交变载荷下萌生裂纹并扩展；
- 二是预紧力控制不当，安装时扭矩过大导致螺栓过载，或扭矩不足引发连接松动，进而产生额外交变载荷；
- 三是材料与制造缺陷，如原材料硫化物夹杂、热处理不当产生的淬火裂纹，或表面处理引发的氢脆问题；
- 四是环境腐蚀作用，如沿海或含硫环境中的氯离子、硫化物引发应力腐蚀开裂（SCC）。
有限元分析技术能够精准捕捉上述失效诱因的力学本质，而ABAQUS的技术特性使其具备独特优势：其一，支持复杂接触关系模拟，可准确计算螺栓与螺母、被连接件之间的接触压力分布，揭示松动失效的力学机理；其二，具备完善的非线性分析功能，能模拟螺栓在复杂工况下的力学响应，定位裂纹萌生源头与扩展路径；其三，支持多物理场耦合分析，可结合腐蚀介质、温度变化等环境因素，开展应力-腐蚀耦合模拟；其四，提供丰富的单元类型与网格技术，能对螺纹、过渡圆角等关键部位进行精细化建模，精准计算应力集中系数。这些特性使ABAQUS成为螺栓连接失效机理研究与工程优化的理想工具之一。接下来我将系统梳理ABAQUS螺栓连接有限元分析的核心流程与技术要点，为工程实践提供技术参考，同时为助力技术人员快速掌握相关技能，欢迎读者朋友首次，点赞和分享到朋友圈。螺栓连接是工程结构中最常见的装配与连接形式之一，其力学行为具有显著的非线性特征。螺栓预紧力、零部件间的接触状态以及摩擦效应共同决定了连接结构的整体刚度、承载能力及局部应力分布。在有限元分析中，螺栓连接问题通常涉及多接触对、多约束方向以及预紧力加载过程，若建模策略不当，极易出现刚度异常、接触失效或应力结果失真等问题。因此，本案例以典型法兰式螺栓连接结构为对象，系统展示 Abaqus 中螺栓连接的几何建模、接触定义、预紧力施加以及结果后处理方法。模型由多个实体部件组成，包括法兰本体、螺栓、螺母以及中间垫片结构。各部件在装配阶段保持几何独立，通过接触关系建立力学联系，而非通过几何合并或绑定方式处理。在建模过程中，通过参考坐标系与基准平面定义装配方向，使螺栓轴线与结构轴线保持一致，为后续约束方向与预紧力施加提供清晰的几何参考。整体结构具有明显的周期性特征，但本案例保留完整模型，以便直观展示多螺栓连接的受力分布。针对不同接触区域，分别建立接触对，并统一纳入接触管理器中进行集中管理，便于后续检查与修改。在接触属性设置中，法向行为采用硬接触形式，保证接触界面在受压时能够有效传递法向力，同时避免非物理的拉力传递。切向行为采用库仑摩擦模型，根据不同接触界面的实际工况设置不同的摩擦系数。例如，在螺栓头与法兰接触面、螺母与垫片接触面中设置较高摩擦系数，而在结构贴合面可根据需要采用较低摩擦系数或近似无摩擦假设。通过为不同接触对分配不同的接触属性，可以更加真实地反映螺栓连接中复杂的力传递路径。在螺栓连接分析中，接触面数量多、几何复杂，主从面选择尤为重要。本案例遵循以下基本原则：所有接触对均采用 Node-to-Surface 离散方式，以保证计算稳定性并降低接触搜索成本。通过接触管理器对主从关系进行统一检查，可有效避免遗漏或定义错误的问题。为避免刚体 位移并保证结构受力合理，需要对模型施加必要的边界约束。本案例中在连接结构底部施加固定约束，用以模拟结构在实际工况中的支承状态。在约束定义过程中，明确区分全局坐标系与局部坐标系，确保约束方向与几何方向一致，避免因坐标系混淆导致约束失效或额外约束刚度的引入。Abaqus 中通过 Bolt Load 功能实现螺栓预紧力的施加。本案例在螺栓中部定义预紧截面，并沿螺栓轴向施加预紧载荷。预紧力通常在独立的分析步中施加，使螺栓产生初始拉伸状态，同时在连接件之间建立接触压力。
通过分步加载，可以有效避免接触突变引起的收敛问题，同时保证预紧力在后续分析中正确保留。螺栓、螺母等小尺寸部件采用相对较密的网格，以捕捉局部应力集中；连接件主体区域则采用相对均匀的中等密度网格。接触区域附近均进行了局部网格加密，以提高接触压力与应力结果的可靠性。在结果分析前，通过整体显示、剖切显示等方式对模型进行检查，确认装配关系、接触状态及约束设置的正确性。计算结果显示，在预紧力作用下，各螺栓呈现出明显的轴向拉应力分布，螺栓头与螺母下方的连接件区域出现局部压应力集中。多个螺栓之间的受力分布具有一定一致性，但仍受装配与接触条件影响存在细微差异。通过整体应力云图可以清晰识别连接结构的主要承载路径，而局部放大显示则有助于分析螺栓孔附近及接触边缘的应力集中情况。通过该螺栓连接分析案例，可以总结出以下关键建模经验：本案例完整展示了 Abaqus 中螺栓连接结构的接触建模与预紧力分析流程。通过合理的接触定义、主从面选择以及预紧力施加策略，可以有效模拟螺栓连接的真实受力状态。该案例可作为复杂装配结构接触分析的重要参考，也为后续疲劳、松动或多工况分析奠定了可靠的建模基础。以下结合风电法兰螺栓失效案例，详细阐述核心技术要点。螺栓连接是工程结构中最常见的装配与连接形式之一，其力学行为具有显著的非线性。几何建模阶段需在保证计算精度的前提下合理简化模型，避免冗余计算。对于螺栓连接装配体，通常包含螺栓、螺母、垫片与被连接件四个核心部件，可通过SolidWorks等CAD软件建立三维模型后导入ABAQUS。建模时需重点保留应力集中敏感区域的几何特征，如螺纹牙型、头下过渡圆角、垫片接触面等，不可随意简化。网格划分是影响计算精度的核心环节，需采用"整体粗化、局部细化"的策略。螺栓、螺母等小尺寸部件采用相对较密的网格，以捕捉局部应力集中；连接件主体区域则采用相对均匀的中等密度网格。接触区域附近均进行了局部网格加密，以提高接触压力与应力结果的可靠性。对于风电法兰螺栓这类关键部件，过渡圆角处的网格尺寸应控制在特征尺寸的1/5以内，螺纹牙型需完整建模，避免因网格粗糙导致应力计算偏差。此外，需注意螺栓与被连接件的网格兼容性，避免出现网格畸变。材料属性定义需严格依据实际螺栓材质参数，对于10.9级风电法兰螺栓常用的42CrMo合金钢，需输入弹性模量210GPa、泊松比0.28、密度7.85×10⁻⁹t/mm³，同时需提供应力-应变曲线，以准确模拟材料塑性变形行为。若涉及疲劳分析，还需输入材料的S-N曲线、疲劳强度系数等参数；若考虑环境腐蚀，需定义腐蚀介质作用下的材料性能衰减模型。接触关系设置直接影响计算结果的准确性，需合理定义接触属性与接触对。在ABAQUS中，通常将螺栓与螺母的螺纹啮合面、螺栓头部与垫片的接触面、垫片与被连接件的接触面定义为接触对。接触属性设置需遵循以下原则：切向行为采用"罚函数摩擦"模型，摩擦系数根据表面处理方式确定，如达克罗处理表面摩擦系数取0.13~0.15，发黑处理表面取0.26~0.30；法向行为采用"硬接触"模型，确保接触面上的压力传递准确，避免出现穿透现象。对于螺纹啮合接触，可采用"面-面接触"或"线-面接触"，若追求高精度，建议采用螺纹实体建模并定义面-面接触。预紧力是螺栓连接的核心力学参数，预紧力不足或过载是导致失效的主要诱因之一，因此准确施加预紧力是有限元分析的关键步骤。ABAQUS中常用两种预紧力施加方法：一是"螺栓载荷法"，通过创建螺栓截面（Section）并定义预紧力大小，适用于简化模型；二是"温度载荷法"，通过给螺栓施加温度变化，利用热胀冷缩原理模拟预紧过程，适用于复杂接触场景。以风电法兰螺栓为例，采用温度载荷法施加预紧力的步骤如下：首先在初始步中定义螺栓与被连接件的接触关系；然后创建通用静力分析步，在该步中给螺栓施加温度降，根据热膨胀系数计算所需温度变化量（ΔT=F/(α·A·E)，其中F为预紧力，α为线膨胀系数，A为螺栓横截面积，E为弹性模量）；最后在后续分析步中施加工作载荷，如法兰的轴向拉力、风力引发的交变载荷等。这种方法能够准确模拟预紧过程中的接触压力分布，以及工作载荷与预紧力的叠加效应，贴近实际安装与工作工况。求解设置需根据分析类型合理选择求解器与迭代参数。对于静力分析（如预紧力计算、过载失效分析），采用ABAQUS/Standard求解器；对于动态疲劳分析或冲击载荷下的失效分析，采用ABAQUS/Explicit求解器。求解过程中需开启接触稳定迭代，避免因接触分离导致计算不收敛。后处理阶段需重点提取与失效相关的关键指标：一是应力分布云图，重点关注过渡圆角、螺纹根部的应力集中情况，计算应力集中系数，判断是否存在疲劳裂纹萌生风险；二是接触压力分布，若接触压力出现明显不均匀或局部为零，说明存在连接松动问题；三是螺栓位移与应变数据，分析是否存在过度塑性变形；四是疲劳寿命预测，结合S-N曲线与交变载荷参数，计算螺栓的疲劳寿命，评估疲劳失效风险。对于风电法兰螺栓失效案例，通过后处理可发现，过渡圆角处的最大主应力达到520MPa，超过材料疲劳极限，且疲劳寿命仅为1.2年，与实际失效情况高度吻合。本文案例模型节选自我的《吃透Abaqus非线性核心：Abaqus接触-屈曲-材料损伤全模块工程案例实操合集》，本课程聚焦 Abaqus 接触、屈曲、材料非线性核心技术，其价值在于帮学习者打通 “理论 - 建模 - 实操 - 落地” 全链路，掌握复杂结构的力学响应分析能力，解决工程与科研中 “接触失效、结构失稳、材料损伤” 等高频难题。应用场景广泛覆盖机械、土木、材料等领域，包括螺栓连接 / 齿轮啮合等接触问题仿真、欧拉柱 / 薄板等结构屈曲稳定性分析、金属塑性变形 / 橡胶超弹性 / 裂纹扩展等材料非线性模拟，适配工程技术人员的项目设计校核、科研人员的课题研究、CAE 工程师的技能升级等核心需求，助力快速提升问题解决效率与实战竞争力。
讲师提供vip交流群/答疑/相关学习资料
《吃透Abaqus非线性核心：Abaqus接触-屈曲-材料损伤全模块工程案例实操合集》
- 球–平板压缩接触分析：压力分布与摩擦效应（插图位置：接触压力云图）
- 掌握ABAQUS中接触、屈曲与材料非线性的建模技巧；
VDI2230 螺栓连接结构标准规范设计计算高级培训（独家更新）
