Meshless solutions of 2D contact problems by subdomain variational inequality and MLPG method with radial basis functions

A subdomain variational inequality and its meshless linear complementary formulation are developed in the present paper for solving two-dimensional contact problems. The subdomain variational inequality will be defined in detail. The meshless method is based on a local weighted residual method with the Heaviside step function as the weighting function over a local subdomain and radial basis functions as trial functions for interpolation. Three different radial basis functions (RBFs), i.e. Multiquadrics (MQ), Gaussian (EXP) and Thin Plate Splines (TPS) are examined and the selection of their shape parameters is studied based on 2D solid stress problems with closed-form solutions. The developed meshless/linear complementary method is applied to solve two frictionless contact problems. For the RBFs, it has been found that the TPS shape parameter is not sensitive to nodal distance and a value of 4 is found as a good choice for TPS from this research.