A macro-scale model of spark plasma sintering (SPS) that couples electrical, thermal, stress-strain and densification components is presented. The continuum theory of sintering is incorporated enabling the evolution of the densification based on local conditions, thus a true spatial density distribution could be obtained. Specimen behavior is described through a non-linear viscous constitutive relation. The simulation is based on an FEM computer code. Several examples are shown and results are compared with experimental data available.