[发明专利]一种杆单元接触计算的处理方法

权利要求书

1.一种杆单元接触计算的处理方法，已知杆单元、实体单元以及相应边界条件的有限元模型，其特征在于包含以下内容：1)在已知的有限元模型内，确定杆单元与实体单元发生接触的区域，该区域为处理区域；2)在确定的处理区域内，沿每个杆单元的轴向，建立壳单元，壳单元与对应的杆单元共节点；3)将步骤2)的壳单元加载进已知的有限元模型，形成第一轮有限元模型；4)在第一轮有限元模型的基础上设定壳单元的厚度并赋予壳单元对应杆单元的材料属性，形成第二轮有限元模型；5)在第二轮有限元模型基础上，建立壳单元与实体单元的接触条件，形成第三轮有限元模型；6)将第三轮有限元模型提交计算，输出该有限元模型的计算结果，包含每个杆单元的应力值；7)对比所有杆单元之间的应力值之差是否在预定的阈值范围内，如果是，完成计算；否则，减小步骤4)中壳单元厚度尺寸，重新形成第二轮有限元模型；再重复步骤5)和步骤6)，直至所有杆单元之间的应力值之差在预定的阈值范围内为止。

2.如权利要求1所述的杆单元接触计算的处理方法，其特征在于，所述的杆单元的轴向两侧分别建立一个相互对称的壳单元。

3.如权利要求1或2所述的杆单元接触计算的处理方法，其特征在于，所述的壳单元应尽可能平行于该杆单元与实体单元的接触面。

4.如权利要求3所述的杆单元接触计算的处理方法，其特征在于，所述的壳单元宽度不超过杆单元长度的十分之一。

5.如权利要求1所述的杆单元接触计算的处理方法，其特征在于，所述的壳单元初始厚度不超过0.01mm。

6.如权利要求1所述的杆单元接触计算的处理方法，其特征在于，所述的杆单元最大应力值和最小应力值差的阈值范围不大于杆单元最小应力值的5％。