论文部分内容阅读
由于材料不可压缩(泊松比趋于0.5)而引起的体积锁定问题一直是有限元计算中的难题之一。在阅读国内外文献的基础上,本文详细介绍了目前通用的一些解决体积锁定的数值方法以及该领域的一些最新进展,以及目前较为广泛研究的光滑有限元法(SFEM),并通过模态分析法分析了常见单元(基于FEM的三角元Tri3、基于FEM的四边元Q4以及基于NS-FEM的三角元Tri3-NS)的不锁定模态、物理及非物理锁定模态。研究表明,Q4、Tri3单元的体积锁定是由于在体积不可压缩条件下,某些基本变形模态(例如:“沙漏模态”)无法准确描述所致。因此,本文受无网格法富集节点域思想的启发,构建了一种新型的5节点四边形单元(Q5),使其在不可压缩条件下能够准确描述基本变形模态。通过Inf-Sup检验发现,目前商用有限元软件通用的选择性降阶积分法(Q4-SRI)无法解决体积锁定引起的“压力波动”问题,主要原因在于单元之间的应变是不连续的,为此基于Q5单元提出了两种光滑积分算法:基于边的面积权重应变光滑(EAW)以及基于节点的面积权重应变光滑(NAW),并导出了一系列基于Q5单元的免于体积锁定的光滑积分策略。随后,基于Q5单元气泡形函数富集节点域的思想,提出了一系列气泡形富集节点域单元,包括:9节点六面体单元(H9单元,适用于三维问题)、4节点三角元(Tri4单元,适用于二维平面应变问题)以及5节点四面体单元(Tet5单元,适用于三维问题),并给出了相应的免于体积锁定的光滑积分策略。大量典型算例结果表明,这种新型气泡形函数富集单元及衍生的一系列光滑积分策略在不可压缩条件下能成功解锁,本文提出的Q5-pEAW/NAW以及H9-Gi/NVW模型能够有效地解决体积自锁问题并通过数值Inf-Sup检验。本文主要工作如下:(1)基于无网格法气泡形函数富集节点域的思路,提出了一种Q5单元,使其能够在高斯点上满足不可压缩条件,解决不可压缩条件下的“沙漏模态”锁定问题。(2)将Q5单元与光滑有限元法(SFEM)相结合,提出了一系列解决体积锁定问题的积分策略:包括基于子单元域的应变光滑法(Q5-CS-SC4)、基于边界域光滑的面积权重应变光滑法(Q5-EAW)、基于边界域光滑的可变面积权重应变光滑法(Q5-pEAW)基于节点域光滑的可变面积权重选择性积分法(Q5-pEAW/NAW)。Q5-CS-SC4法继承了Q5单元采用高斯积分的优点,且具有“软化”刚度矩阵的性质。Q5-EAW法继承了Q5-CS-SC4法解决“沙漏”锁定的特性,同时还兼具光滑单元间应变的特性,可以解决体积锁定导致的压力波动问题。Q5-pEAW法提高了解的准确性。Q5-pEAW/NAW模型兼具Q5-pEAW法以及NAW法的优点,相比于Q5-EAW法更加准确,且更易于推广至三维问题。并推导出了Q5-pEAW/NAW模型的几何非线性算法构造及迭代步骤。(3)基于Q5单元气泡形函数富集节点域的思想,提出了一系列气泡形函数富集节点域单元,包括:H9单元、Tri4单元以及Tet5单元,并给出了相应的免于体积锁定的光滑积分策略。