A data-driven harmonic approach to constructing anisotropic damage models with a minimum number of internal variables