A numerical method for introducing an arbitrary yield function into rigid–viscoplastic FEM programs