Formation and incidence of shell galaxies in the Illustris simulation