Non-cyclic shakedown/ratcheting boundary determination - Part 2: Numerical implementation