Swinging Atwood's machine

The Swinging Atwood's machine. The smaller mass, labelled m, is allowed to swing freely whereas the larger mass, M, can only move up and down. Assume the pivots to be points.

The swinging Atwood's machine (SAM) is a mechanism that resembles a simple Atwood's machine except that one of the masses is allowed to swing in a two-dimensional plane, producing a dynamical system that is chaotic for some system parameters and initial conditions.

Specifically, it comprises two masses (the pendulum, mass $m$ and counterweight, mass $M$ ) connected by an inextensible, massless string suspended on two frictionless pulleys of zero radius such that the pendulum can swing freely around its pulley without colliding with the counterweight.^[1]

The conventional Atwood's machine allows only "runaway" solutions (i.e. either the pendulum or counterweight eventually collides with its pulley), except for $M=m$ . However, the swinging Atwood's machine with $M>m$ has a large parameter space of conditions that lead to a variety of motions that can be classified as terminating or non-terminating, periodic, quasiperiodic or chaotic, bounded or unbounded, singular or non-singular^[1]^[2] due to the pendulum's reactive centrifugal force counteracting the counterweight's weight.^[1] Research on the SAM started as part of a 1982 senior thesis entitled Smiles and Teardrops (referring to the shape of some trajectories of the system) by Nicholas Tufillaro at Reed College, directed by David J. Griffiths.^[3]

Equations of motion

Motion of Swinging Atwood's Machine for M/m = 4.5

The swinging Atwood's machine is a system with two degrees of freedom. We may derive its equations of motion using either Hamiltonian mechanics or Lagrangian mechanics. Let the swinging mass be $m$ and the non-swinging mass be $M$ . The kinetic energy of the system, $T$ , is:

\begin{align} T &= \frac{1}{2} M v^2_M + \frac{1}{2} mv^2_m \\ &= \frac{1}{2}M \dot{r}^2+\frac{1}{2} m \left(\dot{r}^2+r^2\dot{\theta}^2\right) \end{align}

where $r$ is the distance of the swinging mass to its pivot, and $\theta$ is the angle of the swinging mass relative to pointing straight downwards. The potential energy $U$ is solely due to the acceleration due to gravity:

\begin{align} U &= Mgr - mgr \cos{\theta} \end{align}

We may then write down the Lagrangian, $\mathcal{L}$ , and the Hamiltonian, $\mathcal{H}$ of the system:

\begin{align} \mathcal{L} &= T-U\\ &= \frac{1}{2}M \dot{r}^2+\frac{1}{2} m \left(\dot{r}^2+r^2\dot{\theta}^2\right) - Mgr + mgr \cos{\theta}\\ \mathcal{H} &= T+U\\ &= \frac{1}{2}M \dot{r}^2+\frac{1}{2} m \left(\dot{r}^2+r^2\dot{\theta}^2\right) + Mgr - mgr \cos{\theta} \end{align}

We can then express the Hamiltonian in terms of the canonical momenta, $p_r$ , $p_\theta$ :

\begin{align} p_r &= \frac{\partial{\mathcal{L}}}{\partial \dot{r}} = \frac{\partial T}{\partial \dot{r}} = (M+m)\dot{r}\\ p_\theta &= \frac{\partial {\mathcal{L}}}{\partial \dot{\theta}} = \frac{\partial T}{\partial \dot{\theta}} = mr^2 \dot{\theta}\\ \therefore \mathcal{H} &= \frac{p_r^2}{2(M+m)} + \frac{p_\theta^2}{2mr^2} + Mgr - mgr \cos{\theta} \end{align}

Lagrange analysis can be applied to obtain two second-order coupled ordinary differential equations in $r$ and $\theta$ . First, the $\theta$ equation:

\begin{align} \frac{\partial {\mathcal{L}}}{\partial \theta} &= \frac{d}{dt} \left(\frac{\partial {\mathcal{L}}}{\partial \dot{\theta}}\right)\\ -mgr \sin{\theta} &= 2mr \dot{r}\dot{\theta} + mr^2 \ddot{\theta}\\ r\ddot{\theta} + 2\dot{r}\dot{\theta} + g\sin{\theta} &= 0 \end{align}

And the $r$ equation:

\begin{align} \frac{\partial {\mathcal{L}}}{\partial r} &= \frac{d}{dt} \left( \frac{\partial L}{\partial \dot{r}}\right)\\ mr\dot{\theta}^2 - Mg + mg\cos{\theta} &= (M+m) \ddot{r} \end{align}

We simplify the equations by defining the mass ratio $\mu = \frac{M}{m}$ . The above then becomes:

(\mu+1)\ddot{r} - r\dot{\theta}^2 + g(\mu - \cos{\theta}) = 0

Hamiltonian analysis may also be applied to determine four first order ODEs in terms of $r$ , $\theta$ and their corresponding canonical momenta $p_r$ and $p_\theta$ :

\begin{align} \dot{r}&=\frac {\partial{\mathcal{H}}} {\partial{p_r}} = \frac {p_r}{M+m} \\ \dot{p_r} &= - \frac {\partial{\mathcal{H}}} {\partial{r}} = \frac {p_\theta ^2} {mr^3} - Mg + mg\cos{\theta} \\ \dot{\theta}&=\frac {\partial{\mathcal{H}}} {\partial{p_\theta}} = \frac {p_\theta} {mr^2} \\ \dot{p_\theta} &= - \frac {\partial{\mathcal{H}}} {\partial{\theta}} = -mgr\sin{\theta} \end{align}

Notice that in both of these derivations, if one sets $\theta$ and angular velocity $\dot{\theta}$ to zero, the resulting special case is the regular non-swinging Atwood machine:

\ddot{r} = g \frac{1-\mu}{1+\mu}=g\frac{m-M}{m+M}

The swinging Atwood's machine has a four-dimensional phase space defined by $r$ , $\theta$ and their corresponding canonical momenta $p_r$ and $p_\theta$ . However, due to energy conservation, the phase space is constrained to three dimensions.

System with massive pulleys

If the pulleys in the system are taken to have moment of inertia $I$ and radius $R$ , the Hamiltonian of the SAM is then:^[4]

\mathcal{H}\left(r, \theta, \dot{r}, \dot{\theta} \right) = \underbrace{ \frac{1}{2} M_t \left( R \dot{\theta} - \dot{r} \right) ^2 + \frac{1}{2} m r^2 \dot{\theta}^2 }_{T} + \underbrace{ gr \left(M - m \cos{\theta} \right) + gR \left( m \sin{\theta} - M \theta \right)}_{U},

Where M_t is the effective total mass of the system,

M_t = M + m + \frac{I}{R^2}

This reduces to the version above when $R$ and $I$ become zero. The equations of motion are now:^[4]

\begin{align} \mu_t ( \ddot{r} - R \ddot{\theta}) & = r \dot{\theta}^2 + g (\cos {\theta} - \mu ) \\ r \ddot{\theta} & = - 2 \dot{r} \dot{\theta} + R \dot{\theta}^2 - g \sin {\theta} \\ \end{align}

where $\mu_t = M_t / m$ .

Integrability

Hamiltonian systems can be classified as integrable and nonintegrable. SAM is integrable when the mass ratio $\mu = M/m = 3$ .^[5] The system also looks pretty regular for $\mu = 4 n^2 - 1 = 3, 15, 35, ...$ , but the $\mu = 3$ case is the only integrable mass ratio found so far. For many other values of the mass ratio (and initial conditions) SAM displays chaotic motion.

Numerical studies indicate that when the orbit is singular (initial conditions: $r=0, \dot{r}=v, \theta=\theta_0, \dot{\theta}=0$ ), the pendulum executes a single symmetrical loop and returns to the origin, regardless of the value of $\theta_0$ . When $\theta_0$ is small (near vertical), the trajectory describes a "teardrop", when it is large, it describes a "heart". These trajectories can be exactly solved algebraically, which is unusual for a system with a non-linear Hamiltonian.^[6]

Trajectories

The swinging mass of the swinging Atwood's machine undergoes interesting trajectories or orbits when subject to different initial conditions, and for different mass ratios. These include periodic orbits and collision orbits.

Nonsingular orbits

For certain conditions, system exhibits complex harmonic motion.^[1] The orbit is called nonsingular if the swinging mass does not touch the pulley.

Selection of nonsingular orbits
An orbit of the swinging Atwood's machine for $\mu=2$ , $\theta_0=\frac{\pi}{2}$ , and zero initial velocity.	An orbit of the swinging Atwood's machine for $\mu=3$ , $\theta_0=\frac{\pi}{2}$ , and zero initial velocity.	An orbit of the swinging Atwood's machine for $\mu=5$ , $\theta_0=\frac{\pi}{2}$ , and zero initial velocity.	An orbit of the swinging Atwood's machine for $\mu=6$ , $\theta_0=\frac{\pi}{2}$ , and zero initial velocity.
An orbit of the swinging Atwood's machine for $\mu=16$ , $\theta_0=\frac{\pi}{2}$ , and zero initial velocity.	An orbit of the swinging Atwood's machine for $\mu=19$ , $\theta_0=\frac{\pi}{2}$ , and zero initial velocity.	An orbit of the swinging Atwood's machine for $\mu=21$ , $\theta_0=\frac{\pi}{2}$ , and zero initial velocity.	An orbit of the swinging Atwood's machine for $\mu=24$ , $\theta_0=\frac{\pi}{2}$ , and zero initial velocity.

Periodic orbits

Type A orbits for

\theta_0

ranging from 0.1 to 3.1.

When the different harmonic components in the system are in phase, the resulting trajectory is simple and periodic, such as the "smile" trajectory, which resembles that of an ordinary pendulum, and various loops.^[3]^[7] In general a periodic orbit exists when the following is satisfied:^[1]

r(t+\tau) = r(t),\, \theta(t+\tau) = \theta(t)

The simplest case of periodic orbits is the "smile" orbit, which Tufillaro termed Type A orbits in his 1984 paper.^[1]

Selection of periodic orbits
A "smile" orbit of the swinging Atwood's machine for $\mu=1.665$ , $\theta_0=\frac{\pi}{2}$ , and zero initial velocity.	An orbit of the swinging Atwood's machine for $\mu=2.394$ , $\theta_0=\frac{\pi}{2}$ , and zero initial velocity.	An orbit of the swinging Atwood's machine for $\mu=1.1727$ , $\theta_0=\frac{\pi}{2}$ , and zero initial velocity.	An orbit of the swinging Atwood's machine for $\mu=1.555$ , $\theta_0=\frac{\pi}{2}$ , and zero initial velocity.

Singular orbits

The motion is singular if at some point, the swinging mass passes through the origin. Since the system is invariant under time reversal and translation, it is equivalent to say that the pendulum starts at the origin and is fired outwards:^[1]

r(0) = 0

The region close to the pivot is singular, since $r$ is close to zero and the equations of motion require dividing by $r$ . As such, special techniques must be used to rigorously analyze these cases.^[8]

The following are plots of arbitrarily selected singular orbits.

Selection of singular orbits
An orbit of the swinging Atwood's machine for $\mu=10$ , $\theta_0=\frac{\pi}{2}$ , and zero initial velocity.	An orbit of the swinging Atwood's machine for $\mu=25$ , $\theta_0=\frac{\pi}{2}$ , and zero initial velocity.

Collision orbits

Type B orbits for

\theta_0

ranging from 0.1 to 3.1.

Collision (or terminating singular) orbits are subset of singular orbits formed when the swinging mass is ejected from the its pivot with an initial velocity, such that it returns to the pivot (i.e. it collides with the pivot):

r(\tau) = r(0) = 0, \, \tau > 0

The simplest case of collision orbits are the ones with a mass ratio of 3, which will always return symmetrically to the origin after being ejected from the origin, and were termed Type B orbits in Tufillaro's initial paper.^[1] They were also referred to as teardrop, heart, or rabbit-ear orbits because of their appearance.^[3]^[6]^[7]^[8]

When the swinging mass returns to the origin, the counterweight mass, $M$ must instantaneously change direction, causing an infinite tension in the connecting string. Thus we may consider the motion to terminate at this time.^[1]

Boundedness

For any initial position, it can be shown that the swinging mass is bounded by a curve that is a conic section.^[2] The pivot is always a focus of this bounding curve. The equation for this curve can be derived by analyzing the energy of the system, and using conservation of energy. Let us suppose that $m$ is released from rest at $r=r_0$ and $\theta=\theta_0$ . The total energy of the system is therefore:

E = \frac{1}{2}M \dot{r}^2+\frac{1}{2} m \left(\dot{r}^2+r^2\dot{\theta}^2\right) + Mgr - mgr \cos{\theta} = Mgr_0 - mgr_0 \cos{\theta_0}

However, notice that in the boundary case, the velocity of the swinging mass is zero.^[2] Hence we have:

Mgr - mgr \cos{\theta}=Mgr_0 - mgr_0 \cos{\theta_0}

To see that it is the equation of a conic section, we isolate for $r$ :

\begin{align} r&=\frac{h}{1-\frac{\cos{\theta}}{\mu}}\\ h&=r_0\left(1-\frac{\cos{\theta_0}}{\mu}\right) \end{align}

Note that the numerator is a constant dependent only on the initial position in this case, as we have assumed the initial condition to be at rest. However, the energy constant $h$ can also be calculated for nonzero initial velocity, and the equation still holds in all cases.^[2] The eccentricity of the conic section is $\frac{1}{\mu}$ . For $\mu>1$ , this is an ellipse, and the system is bounded and the swinging mass always stays within the ellipse. For $\mu=1$ , it is a parabola and for $\mu<1$ it is a hyperbola; in either of these cases, it is not bounded. As $\mu$ gets arbitrarily large, the bounding curve approaches a circle. The region enclosed by the curve is known as the Hill's region.^[2]

Recent three dimensional extension

A new integrable case for the problem of three dimensional Swinging Atwood Machine (3D-SAM) is announced recently.^[9]

References

1 2 3 4 5 6 7 8 9 Tufillaro, Nicholas B.; Abbott, Tyler A.; Griffiths, David J. (1984). "Swinging Atwood's Machine". American Journal of Physics 52 (10): 895–903. Bibcode:1984AmJPh..52..895T. doi:10.1119/1.13791. |access-date= requires |url= (help)
1 2 3 4 5 Tufillaro, Nicholas B.; Nunes, A.; Casasayas, J. (1988). "Unbounded orbits of a swinging Atwood's machine". American Journal of Physics 56: 1117. Bibcode:1988AmJPh..56.1117T. doi:10.1119/1.15774.
1 2 3 Tufillaro, Nicholas B. (1982). Smiles and Teardrops (Thesis). Reed College.
1 2 Pujol, Olivier; Perez, J.P.; Simo, C.; Simon, S.; Weil, J.A. (2010). "Swinging Atwood's Machine: Experimental and numerical results, and a theoretical study". Physica D 239 (12): 1067–1081. Bibcode:2010PhyD..239.1067P. doi:10.1016/j.physd.2010.02.017.
↑ Tufillaro, Nicholas B. (1986). "Integrable motion of a swinging Atwood's machine". American Journal of Physics 54 (2): 142. Bibcode:1986AmJPh..54..142T. doi:10.1119/1.14710.
1 2 Tufillaro, Nicholas B. (1994). "Teardrop and heart orbits of a swinging Atwoods machine,". American Journal of Physics 62 (3): 231–233. Bibcode:1994AmJPh..62..231T. doi:10.1119/1.17602.
1 2 Tufillaro, Nicholas B. (1985). "Motions of a swinging Atwood's machine". Journal de Physique 46 (9): 1495–1500. doi:10.1051/jphys:019850046090149500.
1 2 Tufillaro, Nicholas B. (1985). "Collision orbits of a swinging Atwood's machine". Journal de Physique 46: 2053–2056. doi:10.1051/jphys:0198500460120205300.
↑ Elmandouh, A. (2016). "On the integrability of the motion of 3D-Swinging Atwood machine and related problems". Physics Letters A 380: 989. doi:10.1016/j.physleta.2016.01.021.

External links

Chaos theory

Chaos theory	Anosov diffeomorphism Bifurcation theory Butterfly effect Chaos theory in organizational development Complexity Control of chaos Dynamical system Edge of chaos Fractal Predictability Quantum chaos Santa Fe Institute Synchronization of chaos Unintended consequences

Chaotic maps (list)	Arnold tongue Arnold's cat map Baker's map Complex quadratic map Complex squaring map Coupled map lattice Double pendulum Double scroll attractor Duffing equation Duffing map Dyadic transformation Dynamical billiards outer Exponential map Gauss map Gingerbreadman map Hénon map Horseshoe map Ikeda map Interval exchange map Kaplan–Yorke map Logistic map Lorenz system Multiscroll attractor Rabinovich–Fabrikant equations Rössler attractor Standard map Swinging Atwood's machine Tent map Tinkerbell map Van der Pol oscillator Zaslavskii map

Chaos systems	Bouncing ball dynamics Chua's circuit Economic bubble Tilt-A-Whirl

Chaos theorists	Michael Berry Leon O. Chua Mitchell Feigenbaum Martin Gutzwiller Brosl Hasslacher Michel Hénon Edward Norton Lorenz Aleksandr Lyapunov Benoît Mandelbrot Henri Poincaré Otto Rössler David Ruelle Oleksandr Mykolaiovych Sharkovsky Floris Takens James A. Yorke

This article is issued from Wikipedia - version of the Monday, February 15, 2016. The text is available under the Creative Commons Attribution/Share Alike but additional terms may apply for the media files.