This study investigates forced nonlinear vibrations of a simply supported Euler-Bernoulli beam on a nonlinear elastic foundation with quadratic and cubic nonlinearities. Applying the homotopy analysis method (HAM) to the spatially discretized governing equation, we derive novel analytical solutions and discuss their convergence to present nonlinear frequency responses with varying contributions of the nonlinearity coefficients. A comparison with numerical solutions is conducted and nonlinear time responses and phase planes are compared to the results from linear beam theory. The study demonstrates that apart from nonlinear problems of free vibrations, HAM is equally capable of solving strongly nonlinear problems of forced vibrations.