通过分析接触前后雅可比矩阵的稀疏结构, 构造基于接触压力判断准则的新型自适应预处理方法。在接触压力为零时使用完全的雅可比矩阵作为Preconditioned Jacobian-free Newton-Krylov(PJFNK)方法的预处理矩阵, 在接触压力非零时使用雅可比矩阵的块对角作为预处理矩阵。数值试验表明: 本文的自适应预处理方法能够有效提升模拟效率, 求解总时间与Multiphysics Object-Oriented Simulation Environment(MOOSE)中现有方法相比减少了15.2%~33.1%。
