论文部分内容阅读
板材成形是一个影响因素众多,涉及材料、几何、接触非线性的复杂大变形过程。为了简化对板材成形问题的研究,研究者通常从线性小变形的问题开始入手,通过构造有限元单元模型对其进行数值仿真计算。由于壳单元考虑了拉伸、弯曲效应,能更精确的模拟材料的变形过程,目前普遍采用壳单元进行板材形变分析。然而,当前主流的有限元法中的壳单元每个节点都有六个自由度,特别是其中的三个旋转自由度给研究带来诸多困难,因为旋转自由度经常会带来不收敛的问题,而且仿真计算效率也会因为节点自由度数的增加而降低。同时,由于材料尺寸变大、形状变复杂的问题愈加突出,节点自由度带来的弊端愈加明显,构造出一种高效的有限元单元显得极为迫切。本文就是针对构造高效的有限元单元做出了如下工作:(1)基于无旋转自由度的思想实现了仅有平移自由度的“S3”板单元。该单元是以“片”存在的,每个主单元有三个相邻的三角形单元,单元每个节点只有3个平移自由度,它的弯曲效应是通过相邻单元的横向位移来表示的。(2)在无旋转自由度的理论基础上,结合离散的Kirchhoff理论,实现了“RFDKT”板单元。该单元是由DKT6板单元和S3单元叠加得到的,通过S3单元的相邻单元关系来去除DKT6单元三边中点的旋转值,从而形成一个有18个自由度的单元“片”。(3)借鉴RFDKT板单元的构造思想,从拟协调的方法出发,将拟协调常弯矩六参数三角形板单元与无旋转自由度的板单元结合起来提出并实现了"RFQC"板单元。该单元与RFDKT单元类似,但是列式更为简单清晰,兼有Rotation Free和拟协调技术的优点,计算效率更高。以上三种单元都是基于无旋转自由度的思想构造出来的,这些单元去除了每个节点的旋转自由度,大大缩小了单元的刚度矩阵,在保证精度的同时其计算效率的优势显著。数值算例的结果表明这三种单元均有良好的精度和收敛性,本文提出的"RFQC"单元计算效率更高。