Numerical implementation of a mathematical model for finite–strain elastoplastic consolidation