“…Another key problem in contact analysis is how to enforce the contact constraint conditions. The penalty function method [1,27], Lagrange multiplier method [5,24,28], augmented Lagrange multiplier method [29] and mathematical programming method [3,4,22,30] are often used to treat contact constraint conditions. The penalty function method is an approximate method, which allows the contact bodies to embed with each other.…”