mirror of
https://asciireactor.com/otho/phy-4660.git
synced 2024-11-23 14:45:06 +00:00
206 lines
15 KiB
TeX
206 lines
15 KiB
TeX
|
|
|
|
|
|
\documentclass[11pt,letterpaper]{article}
|
|
|
|
%\usepackage{natbib}
|
|
%\usepackage{cite}
|
|
\usepackage{biblatex}
|
|
\usepackage{graphicx}
|
|
\usepackage[margin=1.in,centering]{geometry}
|
|
\usepackage{hyperref}
|
|
\usepackage{caption}
|
|
\usepackage[export]{adjustbox}
|
|
\usepackage{float}
|
|
|
|
|
|
\bibliography{/home/caes/wmu/phy-4660/adv_lab.bib}
|
|
|
|
|
|
\begin{document}
|
|
|
|
|
|
\title{Lab 2: Chaos in a Driven Pendulum}
|
|
\author{Otho Ulrich, Eugene Kopf}
|
|
\maketitle
|
|
|
|
\begin{abstract}
|
|
The chaotic behaviour of a driven pendulum is explored. Phase space behaviours of a non-chaotic periodic response and chaotic response are generated from computational models. The form of each is compared, demonstrating the chaotic attractor as a distinguishing feature of chaos. The motion of a physical pendulum is observed under damped and damped-driven conditions. The motion is explored according to its phase space output, and determined to be chaotic.
|
|
\end{abstract}
|
|
|
|
|
|
\section{Chaos}
|
|
\label{sec:chaos}
|
|
Chaos is observed in many non-linear physical systems. It is the condition that a system's evolution is strongly sensitive to initial conditions, such that predicting the future state becomes impossible. The motion of a driven oscillator, such as a driven pendulum, becomes unpredictable as the driving frequency and natural frequency of the pendulum interact. Damping can constrain the motion, and we find that while the motion is unpredictable, it still exhibits some useful features, which will be discussed in section \ref{sec:modeling}.
|
|
|
|
For a finite period of time, chaotic behaviour isn't completely discernable from periodic behaviour, because the possibility exists that the function may repeat itself at some future time. Distinctions are seen between a periodic variable and a chaotic variable in phase space and Poincar\'{e} sections. \cite{TANGLEDTALE}
|
|
|
|
|
|
\subsection{Model of a Driven Pendulum}
|
|
\label{sec:modeling}
|
|
The driven pendulum exemplifies chaotic motion. In the angular coordinate $\theta$, the equation of motion of a driven pendulum is
|
|
|
|
\begin{center}
|
|
$\frac{d^2\theta}{dt^2} = \frac{{\omega_0}^2}{I} sin(\theta) - \frac{{\alpha}}{I} \frac{d\theta}{dt} + \frac{f}{I} cos(\omega t + \phi)$.
|
|
\end{center}
|
|
|
|
Here, $\omega_0$ represents the natural frequency of the pedulum, also its resonant frequency. $\alpha$ is a damping term -- this can take a variety of forms, and in the experiment of section \ref{sec:experiment} is produced by a neodymium magnet interacting with the metal wheel of the pendulum. $f$ is the forcing amplitude where $\omega$ is the forcing frequency, offset from the angular coordinate by a phase $\phi$. $I$ is the moment of inertia of the pendulum. The system will respond most strongly to the driver at this frequency.
|
|
|
|
|
|
\section{Phase Space}
|
|
\label{sec:phasespace}
|
|
A 2-dimensional phase space is a useful environment in which to identify the chaos in the motion in a coordinate. For an oscillator, the convenient coordinate is its angular position $\theta$, with the angular velocity response $\theta^\prime$. In figure \ref{fig:model_periodic}, observe how a periodic variable can be identified in phase space.
|
|
|
|
\begin{figure}
|
|
\includegraphics[width=6.5in]{model_periodic.png}
|
|
\caption{Computer-generated model of a periodic function in phase space. A pendulum's angular coordinate would correspond to $\theta=x$. The motion is predictable and orbits a single point.}
|
|
\label{fig:model_periodic}
|
|
\end{figure}
|
|
|
|
\subsection{Chaotic Attractor}
|
|
Figure \ref{fig:model_damped_phase} demonstrates a damped oscillator, which exhibits stable critical points where the angular velocity goes to zero. Under forced conditions, the angular velocity does not converge to zero, but the motion produces orbits about these critical points in phase space, and we call these points chaotic attractors. Poincar\'{e} first postulated that chaos would be exemplified by complicated paths that roughly follow one of these orbits (the apex) about these attractors. \cite{CHAOSAT50} Attractors are a primary identifying characteristic of chaos, and should be observable in the chaotic motion of a forced pendulum. \cite{CHAOSDYNAMICS} A driven oscillator's path orbits around these critical points but can be seen to jump between them in an unpredictable way along the position coordinate; observe figure \ref{fig:model_driven_phase}.
|
|
|
|
Damping still plays an important role in this chaotic motion; figure \ref{fig:model_driven_nodamp_phase} shows a driven pendulum with no damping, where the motion freely exhibits the effects of both the natural and forcing frequency, and becomes extremely complicated. The natural (un-driven) response of the pendulum is seen in the tall cirular strokes, which represent the pendulum's weight attempting to bring the pendulum to equilibrium. The driving frequency produces small-amplitude variations when the pendulum has a high angular speed and when the pendulum has a low angular speed, it can easily reverse the motion of the pendulum. The oscillator jumps around to many possible states, making it difficult to discern attractors. Chaotic attractors are much more clear when forcing interacts with damping in a system; figure \ref{fig:model_driven_phase} shows this very well.
|
|
|
|
Computer models were generated using the Chaos for Java program written by Brian Davies. \cite{CHAOSFORJAVA} The path is computed from the circle to the triangle.
|
|
|
|
\begin{figure}
|
|
\includegraphics[width=6.5in]{model_damped_phase.png}
|
|
\caption{Model of a damped oscillator in phase space. The critical points are stable, since the velocity approaches zero from all points within the associated region. \cite{CHAOSDYNAMICS}}
|
|
\label{fig:model_damped_phase}
|
|
\end{figure}
|
|
|
|
\begin{figure}
|
|
\includegraphics[width=6.5in]{model_driven_phase.png}
|
|
\caption{Computer-generated model of a damped and driven pendulum's angular motion. The motion jumps between critical points (called attractors in this context), and exhibits a semi-stable low-energy state around one of attractors. Arguments: $\theta=x, f=1, I=1, \omega=2/3, \alpha=0.7, k=\frac{\omega_0^2}{I}=2, \phi=0, \theta(0)=1, \theta^\prime(0)=1$, 150 time steps.}
|
|
\label{fig:model_driven_phase}
|
|
\end{figure}
|
|
|
|
\begin{figure}
|
|
\includegraphics[width=6.5in]{no_drag.png}
|
|
\caption{Computer-generated model of a driven pendulum with no damping. The path through phase space is unstable and could easily take off toward either extreme of the $\theta$ coordinate. This sensitivity is characteristic of chaotic motion. Arguments: $\theta=x, f=1, I=1, \omega=0.5, \alpha=0, k=1.5, \phi=0$, 180 time steps.}
|
|
\label{fig:model_driven_nodamp_phase}
|
|
\end{figure}
|
|
|
|
\subsection{Poincar\'{e} Map}
|
|
A Poincar\'{e} section is a 2D phase space cross-section; in the case of the driven pendulum, the cross-section at a phase $\phi$ from the forcing term. One draws a map between each successive point to create a Poincar\'e map, which is a useful representation of a system's behaviour in phase space. A Poincar\'e section of the driven pendulum model is shown in figure \ref{fig:model_driven_poincare}. An undriven and undamped oscillator will always return to the same point after one period, so a Poincar\'e map sampled using the period corresponding to the oscillator's natural frequency will consist of a single dot. The existence of multiple points at the same phase indicates chaotic motion.
|
|
|
|
\begin{figure}
|
|
\includegraphics[width=6.5in]{model_driven_poincare.png}
|
|
\caption{Poincar\'e section of the damped and driven pendulum from figure \ref{fig:model_driven_phase} at $\phi=0$. Several dots are observed at this sampling phase, so the motion can be deemed chaotic.}
|
|
\label{fig:model_driven_poincare}
|
|
\end{figure}
|
|
|
|
|
|
\section{Experimental Driven Pendulum}
|
|
\label{sec:experiment}
|
|
An experimental driven pendulum is built in order to determine whether its motion is chaotic. The apparatus is shown in figure \ref{fig:apparatus}. The springs are connected so that they both hold about equal tension while the pendulum is at its top-most position, and so that they can stretch over at least one full turn of the pendulum in either direction. The tension was not measured precisely, nor the weight of the pendulum or the damping of the pendulum. The amplitude of the driving arm was $1.45" \pm 0.1"$. A photogate is used to determine the driving period when necessary.
|
|
|
|
\begin{figure}
|
|
\includegraphics[width=6.5in]{appartus.png}
|
|
\caption{Experimental driven pendulum. This angular motion of the round metal pendulum at the top is expected to exhibit chaotic motion. \cite{ADVLABCHAOS}}
|
|
\label{fig:apparatus}
|
|
\end{figure}
|
|
|
|
\subsection{Resonant Frequency}
|
|
To determine the natural frequency of the pendulum, it was released from the top position under damped conditions and its frequency of oscillation was measured. The waveform is shown in figure \ref{fig:exp_damped}. Fourier transformations of the waveform in figure \ref{fig:exp_damped} reveal the natural frequency at where the power spectral density peaks, discounting the low-frequency peak from the flattened waveform due to damping. The natural frequency is estimated to be $0.93\pm0.6 Hz$.
|
|
|
|
\begin{figure}
|
|
\hfill
|
|
{\includegraphics[width=6in]{exp_damped_time.png}}[A]
|
|
\hfill
|
|
{\includegraphics[width=5in]{exp_damped_fourier.png}}[B]
|
|
\hfill
|
|
\caption{[A] The damped oscillations of the experimental pendulum. The fourier transforms of this function reveal the natural frequency of the pendulum, in figure B. [B] The fourier frequency-space representation of the damped pendulum. The peak near $0.9Hz$ is the natural/resonant frequency of this pendulum.}
|
|
\label{fig:exp_damped}
|
|
\end{figure}
|
|
|
|
|
|
\subsection{Periodic Motion}
|
|
To produce periodic motion, the driving arm is run with a driving period $1.45\pm0.02s$. The observed motion is plotted in figure \ref{fig:exp_periodic}. The motion still suggests the possibility of chaos by the variation in the path taken about the critical point, but the Poincar\'e section in figure \ref{fig:exp_periodic} indicates that the orbit is likely converging. The period of the observed motion is consistent with the driving period.
|
|
|
|
\begin{figure}
|
|
\hfill
|
|
[A]{
|
|
\includegraphics[width=1.8in]{exp_periodic_time.png}
|
|
}
|
|
\hfill
|
|
[B]{
|
|
\includegraphics[width=1.8in]{exp_periodic_phase.png}
|
|
}
|
|
\hfill
|
|
[C]{
|
|
\includegraphics[width=1.8in]{exp_periodic_poincare.png}
|
|
}
|
|
\hfill
|
|
\caption{[A] The observed motion when the driving arm is run with period $1.45\pm0.02s$ is periodic at $1.48\pm0.08s$. [B] The periodic nature is recognizable in phase space. [C] The Poincar\'e section shows the orbit is likely converging.}
|
|
\label{fig:exp_periodic}
|
|
\end{figure}
|
|
|
|
|
|
\subsection{Chaotic Motion}
|
|
The driving arm period was increased until the pendulum was exhibiting visibly complex behaviour. The final period was $1.13\pm0.02s$. In figure \ref{fig:exp_chaotic_A} the observed motions are plotted. Two chaotic attractors are strikingly visible in the phase diagram, and the Poincar\'e plot shows significant deviation through the phase.
|
|
|
|
\begin{figure}
|
|
\hfill
|
|
[A]{
|
|
\includegraphics[width=1.8in]{exp_chaotic_A_time.png}
|
|
}
|
|
\hfill
|
|
[B]{
|
|
\includegraphics[width=1.8in]{exp_chaotic_A_phase.png}
|
|
}
|
|
\hfill
|
|
[C]{
|
|
\includegraphics[width=1.8in]{exp_chaotic_A_poincare.png}
|
|
}
|
|
\hfill
|
|
\caption{[A] The observed motion when the driving arm is run with period $1.13\pm0.02s$ is chaotic. [B] Chaotic attractors are recognizable in phase space. [C] The Poincar\'e section shows chaos in the orbit, but that the motions are still constrained to a finite area in phase space.}
|
|
\label{fig:exp_chaotic_A}
|
|
\end{figure}
|
|
|
|
The natural frequency of the pendulum is tested next. The driving period corresponding to the natural frequency was chosen as $1.26\pm0.02s$. This driving frequency found a stronger balance between the two attractors observed at the previous frequency. Figure \ref{fig:exp_chaotic_B} shows plots of the motion. While performing this test, it became clear that the driving arm was near the resonance frequency, as the apparatus began to shake itself to tipping.
|
|
|
|
\begin{figure}
|
|
\hfill
|
|
[A]{
|
|
\includegraphics[width=1.8in]{exp_chaotic_B_time.png}
|
|
}
|
|
\hfill
|
|
[B]{
|
|
\includegraphics[width=1.8in]{exp_chaotic_B_phase.png}
|
|
}
|
|
\hfill
|
|
[C]{
|
|
\includegraphics[width=1.8in]{exp_chaotic_B_poincare.png}
|
|
}
|
|
\hfill
|
|
\caption{[A] The observed motion when the driving arm is run with period $1.26\pm0.02s$ (near the resonance frequency) is chaotic. [B] . [C] The Poincar\'e section shows the same integrated area of constraint in variation as the first test, but more evenly distributed.}
|
|
\label{fig:exp_chaotic_B}
|
|
\end{figure}
|
|
|
|
|
|
One final driving period was chosen for good measure. At $1.15\pm0.02s$, the motion appears closer to sinsuisoidal than in the pervious tests. Nonetheless, the attractors can be observed in the phase diagram and the paths through the Poincar\'e section are limited to the same integrated area.
|
|
|
|
\begin{figure}
|
|
\hfill
|
|
\includegraphics[width=1.8in]{exp_chaotic_C_time.png}
|
|
\hfill
|
|
\includegraphics[width=1.8in]{exp_chaotic_C_phase.png}
|
|
\hfill
|
|
\includegraphics[width=1.8in]{exp_chaotic_C_poincare.png}
|
|
\hfill
|
|
\caption{The observed motion when the driving arm is run with period $1.16\pm0.02s$. The time graph looks significiantly different from the previous tests, but the same attractors and Poincar\'e section are observed.}
|
|
\label{fig:exp_chaotic_C}
|
|
\end{figure}
|
|
|
|
\section{Discussion}
|
|
\label{sec:discussion}
|
|
The periodic motion is observable in the first experimental case, plotted in figure \ref{fig:exp_periodic}. Once the driving arm period was increased from that frequency to one that induced chaotic motion, the predicted spread of intercepts through the Poincar\'e section and the orbits about the chaotic attractors were visually confirmed. It would have been useful to produce the phase diagram of the damped, undriven pendulum to compare the locations of those critical points to the chaotic attractors. However, as Poincar\'e predicted so long ago, it is quite evident that constrained chaotic motion does follow complex paths about attractors, and this motion is evident in the observed motion of the damped, driven pendulum.
|
|
|
|
|
|
\printbibliography
|
|
|
|
|
|
\end{document}
|