A stochastic projection method is employed to obtain the probability distribution of pitch angle of an airfoil in pitch and plunge subject to probabilistic uncertainty in both the initial pitch angle and the cubic spring coefficient of the restoring pitch force. Historically, the selected basis for the stochastic projection method has been orthogonal polynomials, referred to as the polynomial chaos. Such polynomials, however, result in unacceptable computational expense for applications involving oscillatory motion, and a new basis, the Fourier chaos, is introduced for computing limit-cycle oscillations. Unlike the polynomial chaos expansions, which cannot predict limit-cycle oscillations, the Fourier chaos expansions predict both subcritical and supercritical responses even with low-order expansions and high-order nonlinearities. Bifurcation diagrams generated with this new approximate method compare well to Monte Carlo simulations.