螺栓连接仿真，90%的错误出在第一步
上次我让AI独立建一个螺栓连接模型，它交了一个光杆圆柱——没有螺栓头、没有螺母、没有预紧力。模型跑通了，报告465MPa，看起来像模像样。
但结果完全不能用。
不是因为网格不够细，也不是因为接触没设对——而是因为模型从一开始就没有表达正确的传力机制。
这件事让我意识到一个更大的问题：不只是AI，很多工程师在建螺栓连接模型时，也没有想清楚一个最基本的问题——
你这个螺栓连接，到底靠什么传力？
一个被忽略的根本问题
大部分螺栓连接教程一上来就教你怎么设接触对、怎么加预紧力、怎么划网格。这些当然重要。但在动手之前，有一个问题必须先回答：
你的连接设计，在工况载荷下，主要靠哪种机制传递载荷？
这个问题决定了你的建模策略、网格重点、甚至什么简化可以接受、什么简化绝对不行。回答错了，后面每一步都可能走偏。
螺栓连接的传力机制不是单一的，而是三种机制在不同载荷阶段依次登场——而且同一个连接，在加载过程中会从一种机制过渡到另一种。
图1：螺栓连接的三种传力机制
三种传力机制
机制一：摩擦传力（预紧力主导阶段）
原理： 螺栓头和螺母从两侧夹紧板，预紧力产生法向压力，法向压力产生摩擦力。当外载荷小于板间最大摩擦力时，两块板不发生相对滑移——载荷完全靠摩擦传递，螺栓杆不受剪。
类比： 两个人用力夹住一本书，你拽书的边缘——只要拽力小于两个人夹持的摩擦力，书不会滑动。
什么时候主导： 外载荷小于滑移载荷F_slide = μ × F_preload × n（μ=摩擦系数，F_preload=预紧力，n=接触面数）时。
为什么螺杆-孔壁接触可以放松？ 因为孔径比螺杆大（本例1mm间隙），在摩擦阶段板没有滑移，间隙没有被消除——螺杆和孔壁之间实际上不接触。间隙的消除需要板的滑移，而滑移需要外载荷超过摩擦力。 这是一个耦合关系：间隙→滑移→接触→承压，只有前一步发生了，后一步才会启动。
建模含义：
螺栓头和螺母的几何必须建出来——没有它们就没有夹紧力
预紧力必须施加——没有预紧力就没有摩擦
板间接触的法向行为必须正确——法向压力决定摩擦力大小
螺杆-孔壁接触在模型中仍要定义（因为可能滑移），但网格精度可以适当放松
网格重点在板间接触面和螺栓头/螺母的接触边缘
机制二：承压传力（滑移后阶段）
原理： 当外载荷超过滑移载荷，板开始相对滑移。间隙被消除后，螺杆开始接触孔壁，载荷通过孔壁承压传递。此时螺杆承受剪切，孔壁承受挤压。
类比： 那本书终于从两个人手里滑出来了，但书中间有一根针穿过——现在靠针和纸张之间的挤压来抵抗拽力。
为什么承压型需要更多厚度方向网格？ 承压应力沿板厚度方向不是均匀分布的。靠近螺杆接触面的一侧应力最大，沿厚度方向逐渐衰减——这是因为载荷从孔壁传入板后，通过剪切扩散到整个截面。如果厚度方向只有2-3层单元，这个应力梯度会被严重抹平，导致你既看不到真实的峰值，也看不到正确的衰减趋势。4-5层是捕捉这个梯度的最低要求，6-8层才能做网格独立性验证。
建模含义：
孔周网格必须够细——承压应力梯度沿周向和厚度方向都很大
螺杆和孔壁的接触必须精确定义——这是主要传力路径
厚度方向网格至少4-5层——承压应力沿厚度方向有梯度
螺栓头和螺母仍然需要（提供残余夹紧力），但精度要求降低
机制三：螺栓抗拉（轴向受载）
原理： 当外载荷沿螺栓轴向拉伸时，螺杆承受轴向拉力，板间压力减小。当外载荷等于预紧力时，板间压力降为零，连接开始分离。
建模含义：
螺杆的轴向刚度必须准确——不能用梁单元简化（梁单元无法捕捉弯曲和应力集中）
预紧力衰减过程必须捕捉——这是设计的关注点
螺栓头/螺母与板的接触面可能分离——接触状态变化剧烈
需要关注撬力效应——被连接件的弯曲会在螺栓上产生额外拉力
三种机制的建模策略对比
图2：三种传力机制的建模重点差异
怎么判断你的连接处于哪种机制
图3：螺栓连接建模前必须做的判断流程
在动手建模之前，做一个简单的手算：
第一步：计算滑移载荷
F_slide = μ × F_preload × n
比如：μ=0.3，预紧力F_preload=10kN，单剪（n=1）→ F_slide = 0.3 × 10000 × 1 = 3000N
第二步：对比工况载荷
工况载荷5000N > 3000N → 连接会滑移 → 进入承压阶段。
工况载荷2000N < 3000N → 连接不滑移 → 摩擦主导。
第三步：根据判断结果决定建模策略
摩擦主导 → 重点在预紧力和板间接触。
承压主导 → 重点在孔周网格和螺杆-孔壁接触。
很多人跳过了这三步，直接开始建模。 然后发现结果不对，也不知道哪里出了问题。
软件实操：一个正确的螺栓连接模型
为了把理论落地，我用Abaqus演示完整的建模流程。参数如下：
几何、材料、网格
几何关键： 螺栓不是一根杆——必须包含螺栓头、螺杆、螺母三个部分。
网格关键：
单元类型用六面体缩减积分（C3D8R）。网格独立性验证：至少做3级加密，最大应力变化≤5%才算收敛。
接触定义
共7对接触面——除了常规的板-板和螺杆-孔壁，最容易遗漏的是螺栓头/螺母与板面的接触：
摩擦属性：法向Hard接触，切向罚函数μ=0.3。
预紧力施加
两个分析步：
- Step 1（预紧）： Bolt Load = 10kN（"Apply force"）
- Step 2（拉伸）： 固定螺栓长度（"Fix at current length"）+ 施加0.1mm拉伸位移
最常见错误： Step 2继续施加10kN力而不是固定长度——这会导致预紧力不随外载荷变化，不符合物理。
仿真结果深度分析
两个Step各6个增量步，全部收敛。下面不是简单报数字——而是逐项分析每个数据说明了什么物理。
应力分布：四个部件的应力排序
图4：预紧力施加后的Mises应力剖面（Step-1，中心线剖面，螺栓头朝上）
这些数字说明了什么？
① 应力排序始终是：螺栓 > 螺母 > 上板 > 下板。 这不是随机的——它精确反映了传力路径。预紧力从螺栓头→上板→下板→螺母形成闭合回路，离预紧力源越远的部件受力越小。如果你的仿真结果里板的应力比螺栓高，说明预紧力施加有问题。
② 螺栓应力在拉伸后翻倍（245→487 MPa）。 为什么？外载荷方向是X（横向），但螺栓主要承受Z方向拉伸。翻倍不是因为螺栓直接承受了X方向的力——而是板间滑移导致载荷传递路径改变。当板开始滑移，螺栓杆开始接触孔壁，螺栓同时承受轴向拉伸（预紧力残余）和横向剪切（承压传力），Mises应力自然飙升。
③ 上板应力始终大于下板（119 vs 75, 153 vs 109）。 上板6mm比下板8mm薄——同样的力，截面积小应力就大。但更重要的是上板靠近螺栓头——预紧力从螺栓头直接压入上板，应力集中最严重。如果你的结果反过来（下板应力>上板），检查螺栓头和螺母的接触定义是否搞反了。
图5：拉伸载荷后的Mises应力剖面（Step-2，能直接看到螺栓、板内部应力分布和孔周应力集中）
板间滑移：连接传力状态的直接证据
这是整组数据里最有信息量的一行。
Step-1（预紧后）：完全没滑移。 8253个接触点，没有一个点的滑移量超过0.001mm。这说明预紧力产生的摩擦力完全"锁住"了两块板——这正是摩擦型连接的健康状态。
Step-2（拉伸后）：21%的点开始滑移。 最大滑移量0.1mm（等于施加的位移量），1766个点超过0.001mm。这意味着什么？
连接已经从纯摩擦型过渡到了摩擦+承压混合状态。 21%的点在滑移说明搭接区的边缘已经失去了摩擦锁紧——外载荷超过了该区域的摩擦力。但另外79%的点仍然被摩擦锁住。连接处于过渡区，不是纯承压。
这个数据直接验证了建模前手算的预测：0.1mm位移产生的载荷超过了F_slide = 3000N，连接确实开始滑移了。如果你建模型前没做这个手算，你不会知道结果里的0.1mm滑移是"正常的"还是"连接已经失效了"。
接触压力分布：预紧力是否正确施加的验尸报告
接触压力的三个信号：
① 最大CPRESS = 66.6 MPa（预紧后）。 这个数值能反算预紧力是否生效。螺栓头面积≈π×(9²-5²)=201mm²，10kN/201mm²≈49.8MPa。仿真值66.6比理论值高34%——因为应力集中在螺栓头边缘（不是均匀分布），边缘压力大于平均值是正常的。如果你的CPRESS远低于50MPa，说明预紧力没生效。
② Step-2后CPRESS增加50%（66.6→99.9 MPa）。 外载荷让接触压力增加——这意味着板间被进一步压紧。什么情况下会这样？当板被拉伸时，泊松效应导致板变窄，但螺栓夹紧力还在——板变窄让接触面积减小，单位面积压力增大。如果你的CPRESS在加载后反而下降，说明连接可能开始分离——这是危险的信号。
③ 有效接触点从2762降到2602（-5.8%）。 少了160个接触点——这是板间开始分离的边缘效应。分离从搭接区最外端开始（那里离螺栓最远，夹紧力最弱）。如果有效接触点下降超过20%，说明连接大面积分离——结果已经不可信了。
变形模式：每个部件在做什么
变形模式说明：
- 下板被拉开0.1mm，上板不动——载荷正确传递
- 螺栓和螺母只移动了0.005mm（5%的位移量）——螺栓几乎没动
- 0.1mm位移大部分被板间滑移和板弹性变形吸收了
如果你的仿真里螺栓位移很大（跟板一样），说明螺栓头/螺母与板的接触没生效——螺栓没有被"卡住"，模型变成了自由浮动的销轴。
从数据反推：这个连接处于什么状态
把所有证据放在一起，给这个连接做一次"诊断"：
综合判断：连接处于摩擦→承压的过渡区。 搭接区边缘已开始滑移，但中心区域仍被摩擦锁住。螺栓同时承受轴向残余预紧力和横向承压剪切力。如果继续加大载荷，滑移区会扩大到整个搭接面，连接完全进入承压阶段——螺栓应力会进一步飙升，直到剪切失效或孔壁承压失效。
这个诊断过程，就是你拿到仿真结果后应该做的分析——不是只看Max Mises是多少，而是逐项检查每个指标说明了什么物理。
判断结果可信度的检查清单
上次AI犯错复盘
AI建的模型——光杆圆柱、没有螺栓头和螺母、没有预紧力——用三种传力机制框架分析：
- 摩擦传力：完全缺失。 没有夹紧力=没有摩擦。μ=0.3设了也白设。
- 承压传力：部分表达。 螺杆-孔壁接触定义了，但没有预紧力。
- 螺栓抗拉：不适用。
AI的模型只表达了承压传力的一部分。465MPa反映的是无夹紧状态下的孔壁挤压应力，跟真实螺栓连接相差甚远。
用本例的数据对比：有预紧力的模型，螺栓应力487 MPa（拉伸+剪切耦合）；AI的模型没有预紧力，465 MPa只反映孔壁挤压。数字看起来差不多，但物理含义完全不同。
这就是"不知道自己在分析哪种传力机制"的后果——两个数字接近，但回答的是两个不同的问题。
你在做螺栓连接仿真时，有没有遇到过"模型跑通了但结果不对"的情况？
评论区聊聊——你最常踩的螺栓建模坑是什么？你拿到结果后第一个检查的是什么？
