An arbitrary Lagrangian Eulerian finite element method based numerical code for viscoelastic fluids using well-known stabilization techniques (SUPG, DEVSS, log-conformation) is adapted to perform a 3D study of soft systems (drops, elastic particles) suspended in Newtonian and viscoelastic fluids under unbounded shear flow. Since the interface between the suspended objects and the matrix needs to be tracked, a finite element method with SUPG stabilization and second-order time discretization is defined on the interface, with the normal velocity of the interface equal to the normal component of the fluid velocity and a tangential velocity such that the elements on the interface are evenly distributed. This allows the mesh to get rid of the ta...