螺栓连接的有限元分析

1 概述螺栓是机载设备设计中常用的联接件之一。其具有结构简单, 拆装方便,调整容易等优点, 被广泛应用于航空、航天、汽车以及各种工程结构之中。在航空机载环境下，由于振动冲击的影响，设备往往产生较大的过载，对作为紧固件的螺栓带来强度高要求。螺栓是否满足强度要求，关系到机载设备的稳定性和安全性。

传统力学的解析方法对螺栓进行强度校核，主要是运用力的分解和平移原理，解力学平衡方程，借助理论和经验公式，理想化和公式化。没有考虑到连接部件整体性、力的传递途径、部件的局部细节( 如应力集中、应力分布) 等等。通过有限元法，整体建模，局部细化，可以弥补传统力学解析的缺陷。用有限元分析软件

MSC.Patran/MSC.Nastran 提供的特殊单元来模拟螺栓连接，过程更方便，计算更精确，结果更可靠。因此，有限元在螺栓强度校核中的应用越来越广泛。

2 有限元模型的建立

对于螺栓的模拟，有多种模拟方法，如多点约束单元法和梁元法等。

多点约束单元法(MPC)即采用特殊单元RBE2来模拟螺栓连接。在螺栓连接处，设置其中一节点为从节点(Dependent) ，另外一个节点为主节点(Independent) 。主从节点之间位移约束关系使得从节点跟随主节点位移变化。比例因子选为1, 使从节点和主节点位移变化协调一致，从而模拟实际工作状态下，螺栓对法兰的连接紧固作用。

梁元法模拟即采用两节点梁单元Beam其能承受拉伸、剪切、扭转。通过参数设置，使梁元与螺栓几何属性一致。

本文分别用算例来说明这两种方法的可行性。

2.1 几何模型

如图 1 所示组合装配体，底部约束。两圆筒连接法兰通过8 颗螺栓固定。端面受联合载荷作用。

图1三维几何模型

2.2单元及网格

抽取圆筒壁中性面建模，采用四节点壳元（shell）,设置壳元厚度等于实际壁厚。

法兰处的过渡圆弧处网格节点设置密一些，其它可以相对稀疏。

在法兰上下两节点之间建立多点约束单元（RBE2,算例1,图3）或梁元（Beam,算例2,图4）来模拟该位置处的螺栓连接。

图3算例1（多点约束单元法）连接网格

图4算例2（梁元法）连接网格在圆筒端面中心建立不属于结构模型的参考节点，通过加权平均约束单元RBE3 建立端面节点与参考点的主从约束关系。外加载荷施加在参考点上，然后被均匀分配到端面节点。

这里，对于多个面的网格划分，应当注意在各几何连接面法矢量的一致性。这样划分网格时，才能保证shell单元法矢量的一致性。图2显示了各面的法矢量方向是一致的。

图2面法向量方向图

对于复杂曲面模型，还应当注意连接面接缝处网格协调；网格划分结束，必须用Equivale nee合并相同节点。

图5整体模型有限元网格

2.3材料属性、边界约束及载荷

计算中所使用的材料参数如下：

圆筒：E=70 GPa 卩=0.3

螺栓：E=184GPa 卩=0.3

底部法兰在8处螺栓处约束，在独立节点处施加联合载荷。

3有限元结果

3.1应力云图从图6图7看出，两种模拟方法，结构整体应力分布相当

图7算例2（梁元法）应力云图 3.2螺栓强度核算 在两算例中，可以在F06结果文件中得到螺栓对应的节点编号和节点载荷。从结 果文件可以看出，模拟螺栓的两对应节点载荷大小相等、方向相反。所以，只需 取其中一个节点分析即可。下表1、表2以8个上法

Fhrqp

>14kl Msk 酮 Un 2.5

口囲 mWi

2K4+W 2 SSi-C^X

图6算例1（多点约束单元法）应力云图

3.湍4■遡

?aut(K 圍

?4i 心

1.7V+QW UtldJi Filr^e

Max361'3K! :^4F»r 爭 Mkn 負醫T 殴

3

61+0<2.

■-I i 9

Li* T*r

”'

异 5 -rt-

fr i

!曲ErOJr

爲I 匚〜“ xK, 5i|lh>- Tfl

产 I L 诂1- •・.Z72、

伽ODQ

l.W+<

1 Wi

1.221-1

話曲z 沁*

SCWtWl

2 56*0C€

兰节点为例，各节点载荷分量即为单个螺栓所受的载荷，载荷单位No 表1算例1（多点约束单元法）螺栓连接处节点载荷

表2算例2（梁元法）螺栓连接处节点载荷

沖：沏连按螺性的常向載荷・贬值表示螺栓受.

由表可以看出，Fy为连接螺栓的轴向载荷，正值表示螺栓受拉，负值表示螺栓受压缩载荷。而实际工作状况下，连接螺栓是不会受压。表中负值的出现，是由构成单元的两节点之间位移约束特性所决定，这里应当舍负取正。

表1、2中各对应节点Fy值近似相等，Fx和Fz值有所差异。

为了计算方便，以表1（算例1多点约束单元法）为例，分别选取螺栓最大拉伸载荷和螺栓最大剪切载荷计算其相关强度，计算结果偏保守。

螺栓材料1Cr18Ni9Ti，M6

螺栓拉伸载荷：Fy=4194 N

螺栓剪切载荷：

巴=+ 心13+ 602’ =3373 N

螺栓拉伸:

说明螺栓强度满足要求。

4分析与结论 由上分析可知，在有限元分析时，

多点约束单元法和梁元法均可以对装配

螺栓剪切:

3373

3.14 x 6亠

4

=1195/?3 根据第4强度理论:

-V148- + 3X1192

-254MPa

螺栓剩余强度系数:

4194

D~