A numerical model for the simulation of multiphase flows with free surfaces is presented. The model allows to incorporate in a unified manner several phases ranging from incompressible Newtonian flows, Oldroyd-B viscoelastic flows and neo-Hookean elastic solids deformations. We advocate a Eulerian modeling of the multiphase flows, relying on the volume fraction of liquid, describing multiple phases with those different rheologies. One advantage of the Eulerian approach is to allow for large deformations of elastic solids, and changes of topologies. The numerical framework relies on an operator splitting strategy and a two-grid method. The numerical model is validated with a numerical experiment based on the collision between two elastic bodies with free surfaces.