Nonlinear dynamic responses of an Euler–Bernoulli beam attached to a rotating rigid hub with a constant angular velocity under the gravity load are investigated. The slope angle of the centroid line of the beam is used to describe its motion, and the nonlinear integro-partial differential equation that governs the motion of the rotating hub-beam system is derived using Hamilton's principle. Spatially discretized governing equations are derived using Lagrange's equations based on discretized expressions of kinetic and potential energies of the system, yielding a set of second-order nonlinear ordinary differential equations with combined parametric and forced harmonic excitations due to the gravity load. The incremental harmonic balance (IHB) method is used to solve for periodic responses of a high-dimensional model of the system for which convergence is reached and its period-doubling bifurcations. The multivariable Floquet theory along with the precise Hsu's method is used to investigate the stability of the periodic responses. Phase portraits and bifurcation points obtained from the IHB method agree very well with those from numerical integration.