无网格局部彼得洛夫伽辽金法在大变形问题中的应用
【摘要】:
由于网格依赖性,有限元法在计算一些大变形和移动边界问题中遇到了许多困难。一些发展比较成熟的伽辽金型无网格法可以很好地避免这些困难,然而它们需要使用背景网格进行积分,对于大规模问题计算效率较低。本文采用的无网格局部彼得洛夫-伽辽金法(MLPG)具有完全不需要网格、计算速度快和精度高等优点。本文在这种方法的框架内提出了计算非线性材料大变形问题的计算方案,并得到了算例验证,主要工作如下:
1.对MLPG法的计算效率和数值稳定性进行了分析。基于对结点搜索、矩阵存储和求解方案的计算效率进行的分析,提出了优化方案,把MLPG法的算法复杂度降低到了O(N)量级;使用谱分解法对MLPG方法的数值稳定性进行了分析,确定了MLPG方法保持数值稳定的参数范围。
2.提出了一种基于拉格朗日坐标计算路径无关材料大变形问题的MLPG法计算方案,经过通用的线性化处理得到增量求解格式。以超弹性材料为例对本文提出的计算方案进行了验证,算例表明,非线性MLPG法不需要再做其它处理就可以减轻体积锁死问题,能够比有限元法模拟更大的材料变形,并且可以避免有限元由于网格畸变导致精度下降的问题。
3.提出了一种计算路径相关材料大变形问题的MLPG法计算方案。为避免对试探子域和检验子域的更新,基本变量在当前构形的空间坐标中表示,而数值积分在初始构形的物质坐标中进行,对空间坐标的导数通过张量变换转换到物质坐标。为模拟更大的变形采用了极分解的乘法超弹塑性本构模型。作为数值检验,使用这种方法模拟了颈缩、剪切带形成等典型的应变局部化现象。数值算例表明:提出的MLPG计算方案能较精确模拟塑性大变形,可以避免塑性变形体积不变导致体积锁死问题和大变形导致的网格畸变问题;同时,由于无网格形函数可以构造高阶光滑的位移场,MLPG法可以减轻由于应变局部化导致的位移场不连续带来的数值困难。