The ground state properties of the closed-shell 90Zr and closed-subshell 96,106Zr isotopes were investigated using a self-consistent spherical Hartree-Fock HF method and Random Phase Approximation RPA with five different Skyrme type effective nucleon-nucleon interaction: KDE0v1, BsK1, SIII, SVII and SGOI. Having a large number of Skyrme-force parameterizations requires a continuous search for the best for describing the experimental data. The presented approach attempts to accurately describe the structure of the Zirconium isotopes 90, 96, 106Zr. The calculated binding energies, root-mean squares, charge distributions, single particle energies were compared with the experimental data. Moreover, the transition strength function S(E) of the isoscalar giant monopole resonance (ISGMR), Jπ = 0+, T=0, the isoscalar giant dipole resonance (ISGDR), Jπ = 1-, T=0, isoscalar giant quadrupole resonance (ISGQR), Jπ = 2+, T=0, and isoscalar giant octupole resonance (ISGOR), Jπ = 3-, T=0 were calculated for the investigated isotopes and compared with the available experimental data.