In this study, the homotopy perturbation method (HPM) is applied to free vibration analysis of beam on elastic foundation. This numerical method is applied on three different axially loaded cases, namely: 1) one end fixed, the other end simply supported; 2) both ends fixed and 3) both ends simply supported cases. Analytical solutions and frequency factors are evaluated for different ratios of axial load N acting on the beam to Euler buckling load, $N_r$. The application of HPM for the particular problem in this study gives results which are in excellent agreement with both analytical solutions and the variational iteration method (VIM) solutions for all the cases considered in this study and the differential transform method (DTM) results available in the literature for the fixed-pinned case.