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,⁴⁸Ca, ⁵⁶Ni and ²⁰⁸Pb. Spin-orbit splitting in Ca-isotopes, ⁵⁶Ni, ⁹⁰Zr, 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.