This paper investigates some nonlinear dynamical behaviors about domains of attraction, bifurcations, and chaos in an axially accelerating viscoelastic beam under a time-dependent tension and a time-dependent speed. The axial speed and the axial tension are coupled to each other on the basis of a harmonic variation over constant initial values. The transverse motion of the moving beam is governed by nonlinear integro-partial-differential equations with the rheological model of the Kelvin–Voigt energy dissipation mechanism, in which the material derivative is applied to the viscoelastic constitutive relation. The fourth-order Galerkin truncation is employed to transform the governing equation to a set of nonlinear ordinary differential equations. The nonlinear phenomena of the system are numerically determined by applying the fourth-order Runge–Kutta algorithm. The tristable and bistable domains of attraction on the stable steady state solution with a three-to-one internal resonance are analyzed emphatically by means of the fourth-order Galerkin truncation and the differential quadrature method, respectively. The system parameters on the bifurcation diagrams and the maximum Lyapunov exponent diagram are demonstrated by some numerical results of the displacement and speed of the moving beam. Furthermore, chaotic motion is identified in the forms of time histories, phase-plane portraits, fast Fourier transforms, and Poincaré sections.
Two Embedded Pairs of Runge-Kutta Type Methods for Direct Solution of Special Fourth-Order Ordinary Differential Equations
