螺栓连接体的线性分析
理论与物理
螺栓连接的力学
老师，螺栓连接体的FEM分析要怎么做？
螺栓连接是结构分析中最常见的连接形式，但FEM建模却出人意料地困难。螺栓的预紧力（初始张力）、被连接体的压缩、外力引起的轴力变动、接触面的滑动…全都包含非线性因素。
螺栓连接的基本力学
螺栓连接力学的基础是螺栓-被连接体的串联弹簧模型。
螺栓的轴向刚度：$k_b = E_b A_b / L_b$
被连接体的压缩刚度：$k_c$
预紧力（初始张力）$F_i$ 施加后，当外力 $F_{ext}$ 作用时：
这里 $\Phi = k_b / (k_b + k_c)$ 是内力系数（载荷导入率）。
只有外力的一部分（$\Phi$倍）作为附加载荷进入螺栓。剩余部分通过被连接体压缩的缓解而被吸收。
通常的螺栓连接中 $\Phi = 0.1 \sim 0.3$。也就是说只有外力的10~30%进入螺栓。这就是预紧力的效果，能大幅改善螺栓的疲劳寿命。
预紧力对疲劳很重要呢。
是的。没有预紧力，外力会直接作用在螺栓上，导致疲劳断裂。适当的预紧力可以大幅减小螺栓的应力波动幅度。
VDI 2230 指南
螺栓连接有设计标准吗？
VDI 2230的手算仅限于单根螺栓的轴向问题。偏心载荷或多个螺栓的相互作用则需要FEM。
FEM中的建模级别
螺栓连接的FEM建模有三个级别：
| 级别 | 建模方式 | 精度 | 成本 |
|---|---|---|---|
| Level 1 | 弹簧单元＋RBE2/RBE3 | 低（概算） | 低 |
| Level 2 | 梁单元＋接触 | 中 | 中 |
| Level 3 | 实体单元＋接触＋螺纹形状 | 高 | 高 |
Level 3要连螺纹形状都建进去吗？
如果需要精确评估螺纹的应力集中或座面压力，那是必要的。但通常Level 1~2就足够了。Level 3仅用于特殊情况（疲劳评估、高温螺栓等）。
总结
我来整理一下螺栓连接的理论。
要点：
- 串联弹簧模型 — 螺栓刚度 $k_b$ 与被连接体刚度 $k_c$ 的比率是关键
- 内力系数 $\Phi = k_b/(k_b+k_c)$ — 外力的10~30%进入螺栓
- 预紧力改善疲劳寿命 — 减小应力波动幅度
- VDI 2230是设计标准 — 手算针对单根螺栓的轴向问题
- 三个建模级别 — 根据用途选择详细程度
$\Phi$ 的值决定了螺栓设计的一切呢。
只要知道 $\Phi$ 和初始张力 $F_i$，就能计算螺栓的最大轴力、座面压力、疲劳寿命。FEM的目标也是正确评估这两个参数。
螺栓紧固理论的历史
将螺栓轴力（预紧力）概念体系化的是1938年德国VDI（德国工程师协会）发布的VDI 2230标准，这被认为是首次。该标准至今仍在修订，VDI 2230:2015作为世界标准应用于汽车、航空产业。从初版起就存在的核心“紧固系数”概念，在80年后仍未改变。
各项的物理意义
- 惯性项（质量项）：$\rho \ddot{u}$，即“质量×加速度”。您有过急刹车时身体被向前甩出去的经历吗？那种“被带走的感觉”正是惯性力。物体越重越难启动，一旦启动也越难停止。地震时建筑物摇晃，也是因为地面突然移动，而建筑物的质量“被落下”。静力分析中此项设为零，这是基于“缓慢施力，加速度可忽略”的假设。对于冲击载荷或振动问题，此项绝对不能省略。
- 刚度项（弹性恢复力）：$Ku$ 或 $\nabla \cdot \sigma$。拉弹簧时能感觉到“想要恢复原状的力”吧？那就是胡克定律 $F=kx$，也是刚度项的本质。那么提问——铁棒和橡皮筋，用相同的力拉，哪个伸得更长？当然是橡皮筋。这种“不易伸长性”就是杨氏模量 $E$，它决定了刚度。常见的误解：“刚度高=强度高”是不对的。刚度是“不易变形性”，强度是“不易破坏性”，是不同的概念。
- 外力项（载荷项）：体积力 $f_b$（重力等）和表面力 $f_s$（压力、接触力等）。可以这样想——桥上卡车的重量是“作用在整个内容物上的力”（体积力），轮胎压路面的力是“只作用在表面的力”（表面力）。风压、水压、螺栓的紧固力…全都是外力。这里容易犯的错误：弄错载荷方向。本想是“拉伸”却成了“压缩”——听起来像笑话，但在3D空间中坐标系发生旋转时，确实会发生。
- 阻尼项：瑞利阻尼 $C\dot{u} = (\alpha M + \beta K)\dot{u}$。试着弹一下吉他的弦。声音会一直持续吗？不，会逐渐变小对吧。这是因为振动能量通过空气阻力或弦的内部摩擦转化成了热能。汽车的减震器也是同样原理——特意吸收振动能量来改善乘坐舒适性。如果阻尼为零会怎样？建筑物在地震后会一直摇晃不停。实际上不会这样，所以设定适当的阻尼很重要。
假设条件与适用范围
量纲分析与单位制
| 变量 | SI单位 | 注意事项·换算备忘 |
|---|---|---|
| 位移 $u$ | m（米） | 以mm输入时，载荷、弹性模量也需统一为MPa/N系 |
| 应力 $\sigma$ | Pa（帕斯卡）= N/m² | MPa = 10⁶ Pa。与屈服应力比较时注意单位制不一致 |
| 应变 $\varepsilon$ | 无量纲（m/m） | 注意工程应变与对数应变的区别（大变形时） |
| 弹性模量 $E$ | Pa | 钢：约210 GPa，铝：约70 GPa。注意温度依赖性 |
| 密度 $\rho$ | kg/m³ | mm制时为tonne/mm³（钢为 = 10⁻⁹ tonne/mm³） |
| 力 $F$ | N（牛顿） | mm制用N，m制也用N统一 |
数值解法与实现
FEM中的预紧力设置
螺栓的预紧力在FEM中怎么设置？
各求解器都有预紧力专用功能。
Nastran
```
$ 对螺栓梁单元施加预紧力
TEMP(INIT) = 100
TEMP(LOAD) = 200
```
用温度载荷模拟预紧力的方法是经典做法。最近有SOL 400的*BOLT PRELOAD功能。
Abaqus
```
*STEP
*STATIC
*BOLT LOAD
bolt_section, bolt_mid_surface, 50000.
*END STEP
```
*BOLT LOAD是直接对螺栓截面施加初始张力。最直观。
Ansys
```
! Workbench: Bolt Pretension载荷
! APDL: PRETS179单元
SLOAD, bolt_section, PRETEN, 50000.
```
Abaqus的*BOLT LOAD最容易理解呢。切开螺栓截面施加初始张力。
Abaqus的实现中，会在螺栓的中间截面添加一个“预紧力节点”，对该节点施加轴向力。通过将力切换为锁定位移，可以在外力作用后仍保持预紧力。
2步分析
预紧力和外力是同时施加的吗？
基本做法是分2步分析：
- 对螺栓施加初始张力
- 被连接体被压缩
- 螺栓轴力 = $F_i$（设计值）
Step 2: 外力
- 将预紧力“锁定”（固定长度）
- 添加外力（内压、自重、风载荷等）
- 螺栓轴力 = $F_i + \Phi \cdot F_{ext}$
为什么要“锁定”？
如果在Step 2中保持预紧力不变，那么即使因外力变形，预紧力也保持不变（物理上不正确）。切换为锁定（位移固定）后，螺栓轴力会随着外力引起的变形而自然变动。这与VDI 2230的模型是一致的。
接触的必要性
螺栓连接的分析需要接触吗？
本质上是非线性接触问题。被连接面的夹紧力可能因外力而减小，最终可能导致分离（separation）。
相关主题
なった
詳しく
報告