We present a systematic investigation of calculating quantum dots (QDs) energy levels using finite element method in the frame of eight-band k · p method. Numerical results including piezoelectricity, electron and hole levels, as well as wave functions are achieved. In the calculation of energy levels, we do observe spurious solutions (SSs) no matter Burt Foreman or symmetrized Hamiltonians are used. Different theories are used to analyse the SSs, we find that the ellipticity theory can give a better explanation for the origin of SSs and symmetrized Hamiltonian is easier to lead to SSs. The energy levels simulated with the two Hamiltonians are compared to each other after eliminating SSs, different Hamiltonians cause a larger difference on electron energy levels than that on hole energy levels and this difference decreases with the increase of QD size.