Using various versions of the Skyrme force and Negele's interaction, we calculate deformation energies of nuclei by evaluating the expectation value of the many-body Hamiltonian in wave functions taken to be antisymmetrized products of single-particle functions. These single-particle functions are eigenfunctions of a phenomenological potential, here taken to be a deformed Woods-Saxon well. The method can be thought of as an extension of the Strutinsky shell-correction method, to make the connection with the two-body interaction. The method employed here is checked by comparison with Hartree-Fock (HF) results; our method is, however, much faster than the HF method, and, therefore, suitable for a wide range of problems where one tests the sensitivity of results to changes in the two-body interaction. A fairly good agreement with the HF method is obtained for ground-state energies, radii and deformations, as well as for deformations of shape isomers. The main discrepancy is that our energies tend to increase slightly too rapidly with deformation, indicating that we may not have chosen the best phenomenological well. Two-dimensional energy surfaces, which agree quite well with those from the Strutinsky method, are found for 240Pu.