In this paper, the novel model of fluid-conveying imperfect pipe supported at both ends is established by considering the geometric imperfection and the geometric nonlinearity induced by mid-plane stretching. The integral-partial differential equation is discretized by the Galerkin method and solved by a fourth-order Runge–Kutta integration algorithm. Compared with the supercritical pitchfork bifurcation of the perfect pipe conveying fluid, the results show that the cusp bifurcation occurs in the imperfect pipe when increasing the flow velocity. Excellent agreement is observed between the numerical results and the analytical results. The two stable asymmetry bifurcation branches bring interesting phenomena in the post-buckling state. The global nonlinear dynamic behaviors of the imperfect pipe are studied by establishing the bifurcation diagrams. The influence of the geometric imperfection amplitude on the nonlinear response is leading to cusp bifurcation comparing with pitchfork bifurcation of the perfect pipe. When pulsation frequency is set as the bifurcation parameter, there are clear nonresonance ranges, low energy resonance ranges and high energy resonance ranges. In the high energy resonance ranges, the first mode vibration coexisting with the sub-harmonic resonance and combination resonance occurs. As the mean flow velocity and pulsation amplitude are set as bifurcation parameters, the vibration of the imperfect pipe becomes more and more complicated. The vibration exhibits far richer dynamic behaviors including periodic, multi-periodic, quasi-periodic, and chaotic motions. The viscoelastic damping can effectively suppress the vibration response and transfer the high energy resonance to the low energy resonance state. The improved model and corresponding results provide useful information for further studying the dynamic behaviors of fluid-conveying pipe with geometric imperfections.