Element differential method for contact problems with non-conforming contact discretization

In this paper, a new strong-form numerical method, the element differential method (EDM) is employed to solve two- and three-dimensional contact problems without friction. When using EDM, one can obtain the system of equations by directly differentiating the shape functions of Lagrange isoparametric elements for characterizing physical variables and geometry without the variational principle or any integration. Non-uniform contact discretization is used to enhance contact conditions, which avoids performing identical discretization along the contact surfaces of two contact objects. Two methods for imposing contact constraints are proposed. One method imposes Neumann boundary conditions on the contact surface, whereas the other directly applies the contact constraints as collocation equations for the nodes within the contact zone. The accuracy of the two methods is similar, but the multi-point constraints method does not increase the degrees of freedom of the system equations during the iteration process. The results of four numerical examples have verified the accuracy of the proposed method.