A formalism is developed which can be used to calculate nuclear level densities with realistic interactions. The grand partition function is constructed and the level density is obtained by Laplace inversion. The formalism is a generalization of the technique used when the only residual interaction is the pairing force. It takes into account both Hartree-Fock and pairing-type effects. Illustrative examples are given.