论文部分内容阅读
SSOR—PCG方法对于大型对称正定问题具有很高的求解效率，但采用求解静动力接触问题的Lagrange乘子法导致结构刚度矩阵对应Lagrange乘子的对角元为零，不满足传统的SSOR—PCG方法的应用条件．为此通过建立联系Lagrange乘子的罚函数矩阵，提出了SSOR—PCG罚函数方法，并通过主动自由度和被动自由度的关系，提出了SSOR—PCG变量替换法．数值例题证明SSOR—PCG变量替换法具有良好的精度和效率．