On numerical approximation of contact problems