This work is devoted to the development and manufacturing of realistic benchmark phantoms to evaluate the performance of microwave imaging devices. The 3D (3 dimensional) printed phantoms contain several cavities, designed to be filled with liquid solutions that mimic biological tissues in terms of complex permittivity over a wide frequency range. Numerical versions (stereolithography (STL) format files) of these phantoms were used to perform simulations to investigate experimental parameters. The purpose of this paper is two-fold. First, a general methodology for the development of a biological phantom is presented. Second, this approach is applied to the particular case of the experimental device developed by the Department of Electronics and Telecommunications at Politecnico di Torino (POLITO) that currently uses a homogeneous version of the head phantom considered in this paper. Numerical versions of the introduced inhomogeneous head phantoms were used to evaluate the effect of various parameters related to their development, such as the permittivity of the equivalent biological tissue, coupling medium, thickness and nature of the phantom walls, and number of compartments. To shed light on the effects of blood circulation on the recognition of a randomly shaped stroke, a numerical brain model including blood vessels was considered.