Summary: Numerical simulations in three-dimensional for porous media models in soil mechanics are a difficult task for the engineering modelling as well as for numerical realization. Here, we present a general numerical scheme for the simulation of two-phase models in combination with an abstract material model via the stress response with a specialized parallel saddle point solver. Therefore, we give a brief introduction into the theoretical background of the theory of porous media and constitute a two-phase model consisting of a porous solid skeleton saturated by a viscous pore-fluid. The material behaviour of the skeleton is assumed to be elasto-viscoplastic. The governing equations are transfered to a weak formulation suitable for the application of finite element method. Introducing an abstract formulation in terms of the stress response, we define a clear interface between the assembling process and the parallel solver modules. We demonstrate the efficiency of this approach by challenging numerical experiments realized on Linux Cluster in Chemnitz.