A numerical method for coupled electro-mechanical contact problems