We consider a constitutive model for polycrystalline ice, which contains delayed-elastic and viscous deformations, and a damage variable. The damage variable is coupled to the delayed-elastic deformation by a fiber bundle ansatz. We construct an isotropic theory, which can be calibrated with experimental data. Furthermore, we generalize the theory to a damage model in terms of rank-four tensors. This general model allows the evolution of anisotropic damage