Using Skyrme’s density dependent interaction the evolution of nuclear shells has been studied in Hartree-Fock formalism. Optimization of the strength of tensor interaction has been done in reproducing the observed splitting of shell model states of 40,48Ca, 56Ni and 208Pb. Spin-orbit splitting in Ca-isotopes, 56Ni, 90Zr, N = 82 isotones, Sn-isotopes and evolution of gaps in Z, N=8, 20 have been reanalyzed with the inclusion of tensor interaction. For doubly shell closed nuclei it has been observed that tensor interaction is sensitive to spin saturation of nuclear shells.