In this paper, we study the stochastic P-bifurcation problem for axially moving of a bistable viscoelastic beam with fractional derivatives of high order nonlinear terms under Gaussian white noise excitation. First, using the principle for minimum mean square error, we show that the fractional derivative term is equivalent to a linear combination of the damping force and restoring force, so that the original system can be simplified to an equivalent system. Second, we obtain the stationary Probability Density Function (PDF) of the system’s amplitude by stochastic averaging, and using singularity theory, we find the critical parametric condition for stochastic P-bifurcation of amplitude of the system. Finally, we analyze the types of the stationary PDF curves of the system qualitatively by choosing parameters corresponding to each region within the transition set curve. We verify the theoretical analysis and calculation of the transition set by showing the consistency of the numerical results obtained by Monte Carlo simulation with the analytical results. The method used in this paper directly guides the design of the fractional order viscoelastic material model to adjust the response of the system.