二维接触问题的互补算法及工程应用

接触问题是一种常见的非线性问题,如何能够很好地模拟接触面的变形及受力特性,以及实现对变形体间的接触问题的真实模拟是该领域研究的难点问题。基于二维接触问题的实际物理意义,分别在法向和切向建立等价的互补模型。用非线性互补函数(NCP)中的Fischer-Burmeister(FB)函数将互补函数模型转化为非光滑方程组表达,用常规的Newton法求解。同时,基于高斯积分法可以用较少的积分点达到较高的精度,为了进一步提高求解精度,改善不连续的通病,对面-面接触模型在高斯点上对接触面进行处理,可通过调节积分点数目对求解精度进行控制,方法易于理解,实现方便。在此基础上建立二维接触有限元模型,通过一系列工程算例验证该方法的可行性与有效性。结果表明,与ABAQUS有限元的计算结果相比,该方法有着较高的精度,更真实地反映问题的实际。