In this paper, a glycolysis model subject to no-flux boundary condition is considered. First, by discussing the corresponding characteristic equation, the stability of constant steady state solution is discussed, and the Turing's instability is shown. Next, based on Lyapunov-Schmidt reduction method and singularity theory, the multiple stationary bifurcations with singularity are analyzed. In particular, under no-flux boundary condition we show the existence of nonconstant steady state solution bifurcating from a double zero eigenvalue, which is always excluded in most existing works. Also, the stability, bifurcation direction and multiplicity of the bifurcation steady state solutions are investigated by the singularity theory. Finally, the theoretical results are confirmed by numerical simulations. It is also shown that there is no Hopf bifurcation on basis of the condition (C).