并行有限元计算中的接触算法

摘要

针对动态接触问题的有限元并行计算,提出了一种新的接触算法.新算法引入局部拉氏乘子技术来计算接触力.由于同时考虑了无穿透的接触约束条件和相邻接触对的相互影响,较之广泛使用的罚参数法,新算法使接触约束条件和系统平衡方程得到更充分的满足.虽然为提高接触计算精度而在局部采用了迭代技术,但算法仍然具有较高的效率,且与显式时间积分方案完全相容.此外,通过构造专门的区域分解方案,实现了将现有为串行程序开发的搜索算法平滑移植到并行环境的目标.数值算例表明,所提出的接触算法具有很好的并行性,在保证了接触问题并行计算精度的同时,取得了满意的并行效率.