Dance education in colleges and universities is the most important means of inheriting dance skills, cultivating dance talents, and promoting the prosperity and development of dance art. In the new era, the country’s emphasis on “aesthetic education” has provided fertile policy soil for the development of dance majors in universities. Based on the spiritual and cultural needs of the people and the development needs of the national dance art, it is of great urgency for colleges and universities to explore the future oriented Chinese dance higher education and dance creation. Dance education and dance creation are closely linked and interdependent. In the process of Dance education, dance creation inspiration is stimulated. Dance creation and innovation inject new soul into Dance education. College Dance education should combine the two organically to promote the high-quality development of Chinese dance art. The classroom teaching quality evaluation in dance aesthetic education is classical multiple-attributes decision-making (MADM). The probabilistic hesitancy fuzzy sets (PHFSs) are used as a tool for characterizing uncertain information during the classroom teaching quality evaluation in dance aesthetic education. In this paper, we extend the classical grey relational analysis (GRA) method to the probabilistic hesitancy fuzzy MADM with unknown weight information. Firstly, the basic concept, comparative formula and Hamming distance of PHFSs are introduced. Then, the information entropy is used to compute the attribute weights based on the expected values and deviation degree. Then, probabilistic hesitancy fuzzy GRA (PHF-GRA) method is built for MADM under PHFSs. Finally, a practical case study for classroom teaching quality evaluation in dance aesthetic education is designed to validate the proposed method and some comparative studies are also designed to verify the applicability.