Simple lattice dynamical models have been fitted to the experimental phonon dispersion curves of niobium, molybdenum, and their alloys. The 'best-fit' parameters are used to compute the phonon distribution functions, g(v), and several properties which depend on these functions.