Multi‐variant experimental investigations and multi‐physics microstructural modeling of the spark plasma sintering process of metallic powders have been performed up to a relative density of approximately 80 %. In comparison, the effect of sintering temperature, pressure and particle size on the interparticle contact area growth and axial shrinkage of cylindrical specimens of copper and nickel particles is measured in laboratory scaled tests. In the present research, for the first time all relevant for sintering phenomena were considered simultaneously: the fully coupled thermo‐electro‐mechanical modeling of the spark plasma sintering processes, additionally taking into account for lattice, grain boundary, surface diffusion, electromigration and thermomigration has been carried out. The computational analysis of various physical phenomena allows to identify dominant and insignificant mechanisms. The two‐level numerical simulation includes the modeling of the sintering setup at the macroscopic level and the neck formation process in particle chain systems at the microscopic level. The results of the numerical simulations show a very good agreement with the experimental data. Therefore, the impact of electrical and mechanical loads as well as of particle size on microscopic distribution of temperature, inelastic strain and on densification has been studied by the finite element simulations.This article is protected by copyright. All rights reserved.