Polygonal type wear is often seen on a circumference of a rotating viscoelastic disk under contact load. This phenomenon is caused by a combination of physical and mechanical features, such as a dynamic behavior of the material, a heat generation due to loss energy, a heat conduction and various modes of vibrations in the viscoelastic disk including the loading system. The authors had already discussed about the dominant factor of the peculiar phenomenon using an experimental-theoretical hybrid type simulation. In this study, numerical analyses were attempted by using four simplified numerical models, namely a multistage modeling, for a finite element analysis in order to investigate a dominant factor and an opportunity of the periodic stress and temperature distribution of the rotating viscoelastic disk. The criterion of the occurrence of the periodic stress and temperature distribution in a rotating viscoelastic disk was investigated accurately, comparing the results of the experimental/theoretical hybrid simulations and FE analysis.