The same driven damped pendulum

PHYSICS · Multi-part · Fermi University Physics · Mechanics & Dynamical Systems · Nonlinear dynamics and chaos · Nonlinear dynamics and chaos

Computational Physics, Ch When trajectories separate: the driven pendulum and the logistic map

Problem

The same driven damped pendulum, θ¨+γθ˙+ω02sin⁡θ=Adcos⁡(ωdt)\ddot\theta+\gamma\dot\theta+\omega_0^2\sin\theta=A_d\cos(\omega_d t) with ω0=1 s−1\omega_0=1\,\mathrm{s^{-1}}, γ=0.5 s−1\gamma=0.5\,\mathrm{s^{-1}}, ωd=2/3 s−1\omega_d=2/3\,\mathrm{s^{-1}}, is now driven at Ad=1.2 s−2A_d=1.2\,\mathrm{s^{-2}}. Air resistance beyond the stated linear damping is ignored. ((a)) The pendulum is released from θ=0.2 rad\theta=0.2\,\mathrm{rad} at rest and integrated by fourth-order Runge-Kutta with 942942 steps to the drive period, Δt=0.0100 s\Delta t=0.0100\,\mathrm{s}. After discarding a transient of 100100 drive periods, the state is sampled once per drive period for ten further periods, with the angle wrapped to (−π,π](-\pi,\pi], giving θ\theta (rad): 1.1611.161, −0.954-0.954, −0.116-0.116, 0.4290.429, −2.512-2.512, 1.9021.902, 0.2440.244, 2.7752.775, 0.2580.258, 1.1541.154, with no two values repeating and no obvious cycle among them. These ten numbers belong to that one run: a hundred drive periods is 942 s942\,\mathrm{s}, far beyond this pendulum's predictability horizon, so a different step or a different starting angle gives ten different numbers, and it is their pattern rather than their values that the question is about. Repeating this for 10510^5 drive periods gives 10510^5 points that never coincide but stay confined to a bounded region of the (θ,θ˙)(\theta,\dot\theta) plane. What do these sampled points show about the periodicity of the motion? ((b)) Describe the shape traced out in the (θ,θ˙)(\theta,\dot\theta) phase portrait once the transient is discarded, and contrast it with the phase portraits of period-1 and period-2 motion.

Browse the problem library